Block matrix solver for well system fluid flow modeling

Patent metadata

Field	Value
Publication number	US-9206671-B2
Application number	US-201314143958-A
Country	US
Kind code	B2
Filing date	Dec 30, 2013
Priority date	Aug 27, 2013
Publication date	Dec 8, 2015
Grant date	Dec 8, 2015

Abstract

In some aspects, techniques and systems for operating a subterranean region model are described. A block matrix is identified. The block matrix includes a banded portion and an augmented portion. The banded portion includes flow variable coefficients based on governing flow equations of a model representing well system fluid flow in a subterranean region. In a row of the block matrix, decomposition of a pivot element in the row is performed. Elements along the row are modified based on an inverse of the decomposition of the pivot element. A flow variable of the model is solved based on the block matrix that contains the modified elements.

First claim

The invention claimed is:

1. A method of operating a subterranean region model, the method comprising:
identifying a block matrix, the block matrix including a banded portion and an augmented portion, the banded portion including flow variable coefficients based on governing flow equations of a model representing well system fluid flow in a subterranean region;
in a row of the block matrix, performing, by operation of one or more computers, decomposition of a pivot element in the row;
modifying, by operation of the one or more computers, elements along the row based on an inverse of the decomposition of the pivot element; and
solving, by operation of the one or more computers, a flow variable of the model based on the block matrix comprising the modified elements.

2. The method of claim 1 , wherein each element of the block matrix is a sub-matrix, and the decomposition is computed without computing the inverse of the pivot element.

3. The method of claim 1 , wherein performing decomposition of the pivot element comprises performing LU-decomposition of the pivot element.

4. The method of claim 1 , wherein the pivot element resides in a main diagonal of the block matrix.

5. The method of claim 1 , wherein modifying the elements along the row based on the inverse of the decomposition of the pivot element comprises multiplying other elements along the row by the inverse of the decomposition of the pivot element.

6. The method of claim 1 , comprising eliminating non-zero elements below the pivot element by row operations.

7. The method of claim 6 , wherein eliminating non-zero elements below the pivot element comprises eliminating k non-zero elements below the pivot element, k being the smaller of a bandwidth of the banded portion of the block matrix and a distance between the row and a last row of the banded portion of the block matrix.

8. The method of claim 1 , wherein the augmented portion comprises more than one column.

9. The method of claim 1 , wherein each row of the block matrix is stored in a row-major format.

10. The method of claim 1 , wherein each block of the block matrix is stored in a column-major format.

11. The method of claim 1 , wherein the block matrix further comprises off-band elements, the off-band elements representing fluid flow interactions in a fracture network of the subterranean region, and solving the flow variable comprises accounting for the off-band elements.

12. A non-transitory computer-readable medium storing instructions that, when executed by data processing apparatus, perform operations comprising:
identifying a block matrix, the block matrix including a banded portion and an augmented portion, the banded portion including flow variable coefficients based on governing flow equations of a model representing well system fluid flow in a subterranean region;
in a row of the block matrix, performing decomposition of a pivot element in the row;
modifying elements along the row based on an inverse of the decomposition of the pivot element; and
solving a flow variable of the model based on the block matrix comprising the modified elements.

13. The computer-readable medium of claim 12 , wherein performing decomposition of the pivot element comprises performing LU-decomposition of the pivot element.

14. The computer-readable medium of claim 12 , wherein modifying the elements along the row based on the inverse of the decomposition of the pivot element comprises multiplying other elements along the row by the inverse of the decomposition of the pivot element.

15. The computer-readable medium of claim 12 , the operations comprising eliminating non-zero elements below the pivot element by row operations.

16. The computer-readable medium of claim 12 , wherein the augmented portion comprises more than one column.

17. A flow modeling system comprising one or more computers that include:
memory operable to store a block matrix including a banded portion and an augmented portion, the banded portion including flow variable coefficients based on governing flow equations of a model representing well system fluid flow in a subterranean region; and
data processing apparatus operable to:
perform, in a row of the block matrix, decomposition of a pivot element in the row;
modify elements along the row based on an inverse of the decomposition of the pivot element; and
solve a flow variable of the model based on the block matrix comprising the modified elements.

18. The flow modeling system of claim 17 , wherein performing decomposition of the pivot element comprises performing LU-decomposition of the pivot element.

19. The flow modeling system of claim 17 , wherein modifying the elements along the row based on the inverse of the decomposition of the pivot element comprises multiplying other elements along the row by the inverse of the decomposition of the pivot element.

20. The flow modeling system of claim 17 , the data processing apparatus being operable to eliminate non-zero elements below the pivot element by elementary row operations.

Assignees

• Halliburton Energy Services Inc

Inventors

• Shetty Dinesh Ananda
• Lin Avi

Classifications

• G06G7/57

  for fluid flow (G06G7/50 takes precedence){; for distribution networks} · 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

• 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

• G06F30/20

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