Multi-thread band matrix solver for well system fluid flow modeling

Patent metadata

Field	Value
Publication number	US-9284820-B2
Application number	US-201314144034-A
Country	US
Kind code	B2
Filing date	Dec 30, 2013
Priority date	Aug 27, 2013
Publication date	Mar 15, 2016
Grant date	Mar 15, 2016

Abstract

In some aspects, techniques and systems for operating a subterranean region model are described. A band matrix that represents well system fluid flow in a subterranean region is accessed. The band matrix includes flow variable coefficients based on governing flow equations for the well system fluid flow. A first elimination is performed on an upper part of the band matrix. The upper part of the band matrix includes a first subset of the flow variable coefficients corresponding to a first subset of flow variables. In parallel with the first elimination, a second elimination is performed on a lower part of the band matrix. The lower part of the band matrix includes a second subset of flow variable coefficients corresponding to a second subset of the flow variables. An intermediate part of the band matrix can be solved. The intermediate part of the band matrix includes other flow variable coefficients.

First claim

The invention claimed is:

1. A method of operating a subterranean region model, the method comprising:
accessing a band matrix that represents well system fluid flow in a subterranean region, the band matrix including flow variable coefficients based on governing flow equations for the well system fluid flow;
performing, by operation of one or more computers, a first elimination on an upper part of the band matrix, the upper part of the band matrix including a first subset of the flow variable coefficients corresponding to a first subset of flow variables in the governing flow equations;
performing, in parallel with the first elimination, a second elimination on a lower part of the band matrix, the lower part of the band matrix including a second subset of flow variable coefficients corresponding to a second subset of the flow variables in the governing flow equations; and
solving an intermediate part of the band matrix, the intermediate part of the band matrix including other flow variable coefficients of the band matrix, wherein the intermediate part comprises at least ρ rows, ρ being a bandwidth of the band matrix.

2. The method of claim 1 , comprising:
solving the upper part of the band matrix based on a solution of the intermediate part; and
solving the lower part of the band matrix based on the solution of the intermediate part.

3. The method of claim 1 , comprising monitoring a distance between the upper part of the band matrix and the lower part of the band matrix.

4. The method of claim 1 , wherein performing the first elimination on the upper part of the band matrix comprises converting the upper part of the band matrix to an upper triangular representation by forward elimination.

5. The method of claim 1 , wherein performing the second elimination on the lower part of the band matrix comprises converting the lower part of the band matrix to a lower triangular representation by backward elimination.

6. The method of claim 1 , wherein:
performing the first elimination on the upper part of the band matrix comprises performing, by a first thread, the first elimination on the upper part of the band matrix; and
performing the second elimination on the lower part of the band matrix comprises performing, in parallel with the first elimination, by a second thread, the second elimination on the lower part of the band matrix.

7. The method of claim 6 , comprising using a shared memory between the first thread and the second thread.

8. The method of claim 6 , comprising asynchronously communicating between the first thread and the second thread.

9. The method of claim 6 , comprising solving an intermediate part of the band matrix, the intermediate part of the band matrix including other flow variable coefficients of the band matrix, and wherein a matrix position of the intermediate part depends on respective processing powers of the first and second threads.

10. A non-transitory computer-readable medium storing instructions that, when executed by data processing apparatus, perform operations comprising:
performing a first elimination on an upper part of a band matrix, the band matrix including flow variable coefficients based on governing flow equations representing well system fluid flow in a subterranean region, the upper part of the band matrix including a first subset of the flow variable coefficients corresponding to a first subset of flow variables in the governing flow equations;
performing, in parallel with the first elimination, a second elimination on a lower part of the band matrix, the lower part of the band matrix including a second subset of flow variable coefficients corresponding to a second subset of the flow variables in the governing flow equations; and
solving an intermediate part of the band matrix, the intermediate part of the band matrix including other flow variable coefficients of the band matrix, wherein the intermediate part comprises at least ρ rows, ρ being a bandwidth of the band matrix.

11. The computer-readable medium of claim 10 , the operations comprising:
solving the upper part of the band matrix based on a solution of the intermediate part; and
solving the lower part of the band matrix based on the solution of the intermediate part.

12. The computer-readable medium of claim 10 , the operations comprising monitoring a distance between the upper part of the band matrix and the lower part of the band matrix.

13. The computer-readable medium of claim 10 , wherein performing the first elimination on the upper part of the band matrix comprises converting the upper part of the band matrix to an upper triangular representation by forward elimination.

14. The computer-readable medium of claim 10 , wherein performing the second elimination on the lower part of the band matrix comprises converting the lower part of the band matrix to a lower triangular representation by backward elimination.

15. The computer-readable medium of claim 10 , wherein:
performing the first elimination on the upper part of the band matrix comprises performing, by a first thread, the first elimination on the upper part of the band matrix; and
performing the second elimination on the lower part of the band matrix comprises performing, in parallel with the first elimination, by a second thread, the second elimination on the lower part of the band matrix.

16. The computer-readable medium of claim 15 , the operations comprising asynchronously communicating between the first thread and the second thread.

17. The computer-readable medium of claim 15 ,
wherein a matrix position of the intermediate part depends on respective processing powers of the first and second threads.

18. A flow modeling system comprising one or more computers that include:
memory operable to store a band matrix that includes flow variable coefficients based on governing flow equations representing well system fluid flow in a subterranean region; and
data processing apparatus operable to:
perform a first elimination on an upper part of the band matrix, the upper part of the band matrix including a first subset of the flow variable coefficients corresponding to a first subset of flow variables in the governing flow equations;
perform, in parallel with the first elimination, a second elimination on a lower part of the band matrix, the lower part of the band matrix including a second subset of flow variable coefficients corresponding to a second subset of the flow variables in the governing flow equations; and

solve an intermediate part of the band matrix, the intermediate part of the band matrix including other flow variable coefficients of the band matrix, wherein the intermediate part comprises at least ρ rows, ρ being a bandwidth of the band matrix.

19. The flow modeling system of claim 18 , the data processing apparatus being operable to:
solve the upper part of the band matrix based on a solution of the intermediate part; and
solve the lower part of the band matrix based on the solution of the intermediate part.

20. The flow modeling system of claim 18 , the data processing apparatus operable to monitor a distance between the upper part of the band matrix and the lower part of the band matrix.

21. The flow modeling system of claim 18 , the data processing apparatus being operable to execute:
a first thread to perform the first elimination on the upper part of the band matrix; and
a second thread to perform the second elimination on the lower part of the band matrix in parallel with the first elimina

Assignees

• Halliburton Energy Services Inc

Inventors

• Shetty Dinesh Ananda
• Lin Avi

Classifications

• G06F17/10

  Complex mathematical operations {(function generation by table look-up G06F1/03; evaluation of elementary functions by calculation G06F7/544)} · CPC title

• G06F17/16

  Matrix or vector computation {, e.g. matrix-matrix or matrix-vector multiplication, matrix factorization (matrix transposition G06F7/78)} · CPC title

• G06F17/11

  for solving equations {, e.g. nonlinear equations, general mathematical optimization problems (optimization specially adapted for a specific administrative, business or logistic context G06Q10/04)} · CPC title

• E21B43/00

  Methods or apparatus for obtaining oil, gas, water, soluble or meltable materials or a slurry of minerals from wells (applicable only to water E03B) · CPC title

• G06F30/20

  Design optimisation, verification or simulation (optimisation, verification or simulation of circuit designs G06F30/30) · CPC title