Method For Solving Implicit Reservoir Simulation Matrix

• US 20070255779A1
• Filed: 04/13/2005
• Published: 11/01/2007
• Est. Priority Date: 06/07/2004
• Status: Active Grant

• Current Assignee: James Watts
• Original Assignee: James Watts
• Inventors: Watts, James
• Granted Patent: US 7,672,818 B2
• Time in Patent Office: Days
• US Class Current: 708/520
• CPC Class Codes: E21B 43/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 Class Codes: G06F 17/12 Simultaneous equations , e.g. systems of linear equations
• CPC Class Codes: G06F 17/16 Matrix or vector computation , e.g. matrix-matrix or matrix-vector multiplication, matrix factorization matrix transposition G06F7/78
• CPC Class Codes: G06F 30/20 Design optimisation, verification or simulation optimisation, verification or simulation of circuit designs G06F30/30

Abstract

A method for solving a matrix equation AX=B, wherein A represents a block sparse matrix, B represents a right hand side block vector and X represents a solution block vector. In one embodiment, the method includes receiving the block sparse matrix and the right hand side block vector, constructing a reduced transformed block sparse matrix from the block sparse matrix, constructing a reduced transformed residual block vector from the block sparse matrix and the right hand side block vector, and solving for the solution block vector using the reduced transformed block sparse matrix and the reduced transformed residual block vector.

Claims

1. A method for solving a matrix equation AX=B, wherein A represents a block sparse matrix, B represents a right hand side block vector and X represents a solution block vector, the method comprising:

receiving the block sparse matrix and the right hand side block vector;

constructing a reduced transformed block sparse matrix from the block sparse matrix;

constructing a reduced transformed residual block vector from the block sparse matrix and the right hand side block vector;

wherein constructing the reduced transformed block sparse matrix comprises;

assembling one or more coefficients of one or more mass change terms from a column of the block sparse matrix into a column matrix;

performing a singular value decomposition on the column matrix to generate a left matrix, a diagonal matrix and a right matrix; and

discarding one or more singular values that are less than a predetermined threshold value to generate a reduced diagonal matrix; and

solving for the solution block vector using the reduced transformed block sparse matrix and the reduced transformed residual block vector.

2. The method of claim 1, wherein the matrix equation represents a system of fluid flow equations in one or more dimensions having one or more pressure change terms and the one or more mass change terms, wherein the block sparse matrix contains one or more coefficients of the pressure change terms and the one or more coefficients of the mass change terms.

3. The method of claim 1, wherein constructing the reduced transformed block sparse matrix further comprises constructing a transformed block sparse matrix having the same block structure and submatrix form as the block sparse matrix.

4. The method of claim 1, wherein the matrix equation represents a system of fluid flow equations in one or more dimensions having one or more pressure change terms and the one or more mass change terms, wherein the block sparse matrix contains one or more coefficients of the pressure change terms and the one or more coefficients of the mass change terms, and wherein the diagonal matrix comprises one or more singular values arranged in ascending order.

5. The method of claim 1, wherein the matrix equation represents a system of fluid flow equations in one or more dimensions having one or more pressure change terms and the one or more mass change terms, wherein the block sparse matrix contains one or more coefficients of the pressure change terms and the one or more coefficients of the mass change terms, and wherein constructing the reduced transformed block sparse matrix comprises multiplying the left matrix by the reduced diagonal matrix to generate a temporary column matrix.

6. The method of claim 1, wherein the matrix equation represents a system of fluid flow equations in one or more dimensions having one or more pressure change terms and the one or more mass change terms, wherein the block sparse matrix contains one or more coefficients of the pressure change terms and the one or more coefficients of the mass change terms, and wherein constructing the reduced transformed block sparse matrix comprises:

multiplying the left matrix by the reduced diagonal matrix to generate a temporary column matrix having a plurality of mass change terms coefficients; and

assembling the mass change terms coefficients of the temporary column matrix into one or more temporary mass change terms coefficient submatrices that correspond to the mass change terms coefficients of the column matrix.

7. The method of claim 6, wherein constructing the reduced transformed block sparse matrix comprises premultiplying each temporary mass change terms coefficient submatrix by the transpose of the right matrix to generate a transformed temporary mass change terms coefficient submatrix.

8. The method of claim 7, wherein constructing the reduced transformed block sparse matrix comprises premultiplying each subvector containing the pressure change terms coefficients in the block sparse matrix by the transpose of the right matrix to generate a transformed pressure change terms coefficients subvector.

9. The method of claim 8, wherein constructing the reduced transformed block sparse matrix comprises postmultiplying each fluid partial volumes subvector in the block sparse matrix by the right matrix to generate a transformed fluid partial volumes vector.

10. The method of claim 9, wherein constructing the reduced transformed block sparse matrix comprises constructing a transformed block sparse matrix from one or more of the transformed temporary mass change terms coefficient submatrix, the transformed pressure change terms coefficients subvector and the transformed fluid partial volumes vector.

11. The method of claim 10, wherein constructing the reduced transformed block sparse matrix comprises eliminating one or more transformed fluid partial volumes that correspond to the discarded singular values in the transformed block sparse matrix to generate the reduced transformed block sparse matrix.

12. The method of claim 1, wherein constructing the reduced transformed block sparse matrix comprises:

determining a column matrix norm for the column matrix;

skipping the singular value decomposition if the column matrix norm is equal to or less than the predetermined threshold value;

setting the singular values to zero if the column matrix norm is equal to or less than the predetermined threshold value; and

setting the right matrix to the identity matrix if the column matrix norm is equal to or less than the predetermined threshold value.

13. The method of claim 12, wherein determining the column matrix norm comprises calculating the Frobenius norm of the column matrix.

14. The method of claim 1, wherein the reduced transformed block sparse matrix comprises one or more reduced transformed diagonal submatrices and one or more reduced transformed off-diagonal submatrices, wherein each reduced transformed diagonal submatrix comprises mass change terms coefficients and pressure change terms coefficients only within the bottom r_i+1 rows and right most r_i+1 columns of each transformed diagonal submatrix.

15. The method of claim 14, wherein each reduced transformed off-diagonal submatrix comprises mass change terms coefficients and pressure change terms coefficients only within the bottom r_i+1 rows and the right most r_j+1 columns of each transformed off-diagonal submatrix.

16. The method of claim 1, wherein constructing the reduced transformed residual block vector comprises constructing a transformed residual block vector.

17. The method of claim 1, wherein constructing the transformed residual block vector comprises constructing a transformed residual block vector having a transformed mass balance residual subvector and a transformed volume constraint residual subvector.

18. The method of claim 1, wherein constructing the reduced transformed residual block vector comprises premultiplying a mass balance residual subvector by the transpose of the right matrix to generate a transformed mass balance residual subvector.

19. The method of claim 1, wherein solving for the solution block vector comprises solving for a reduced transformed solution change block vector using the reduced transformed block sparse matrix and the reduced transformed residual block vector.

20. The method of claim 1, wherein solving for the solution block vector comprises:

solving for a reduced transformed solution change block vector using the reduced transformed block sparse matrix and the reduced transformed residual block vector; and

converting the reduced transformed solution change block vector to a solution change block vector.

21. The method of claim 1, wherein solving for the solution block vector comprises

adding the solution change block vector to a current estimate of the solution block vector to update the solution block vector.

22. The method of claim 1, wherein solving for the solution block vector comprises:

solving for a reduced transformed solution change block vector using the reduced transformed block sparse matrix and the reduced transformed residual block vector; and

converting the reduced transformed solution change block vector to a solution change block vector having one or more changes in mass unknowns and one or more changes in pressure unknowns.

23. A method for solving a matrix equation AX=B, wherein A represents a block sparse matrix, B represents a right hand side block vector and X represents a solution block vector, the method comprising:

constructing a reduced transformed block sparse matrix from the block sparse matrix;

constructing a reduced transformed residual block vector from the block sparse matrix and the right hand side block vector;

solving for a reduced transformed solution change block vector using the reduced transformed block sparse matrix and the reduced transformed residual block vector;

converting the reduced transformed solution change block vector to a solution change block vector having one or more changes in mass unknowns and one or more changes in pressure unknowns; and

adding the solution change block vector to a current estimate of the solution block vector to update the solution block vector.

24. A method for solving a matrix equation AX=B, wherein A represents a block sparse matrix, B represents a first block vector and X represents a solution block vector, the method comprising:

receiving the block sparse matrix and the first block vector;

constructing a reduced transformed block sparse matrix from the block sparse matrix;

constructing a reduced transformed residual block vector from the block sparse matrix and the first block vector;

wherein constructing the reduced transformed block sparse matrix comprises;

assembling at least one coefficient associated with the at least one mass change term from a column of the block sparse matrix into a column matrix;

performing a singular value decomposition on the column matrix to generate a first matrix, a second matrix and a third matrix; and

discarding each singular value less than a predetermined threshold value to generate a reduced diagonal matrix; and

solving for the solution block vector using the reduced transformed block sparse matrix and the reduced transformed residual block vector.

25. The method of claim 24, wherein the matrix equation represents fluid flow equations in a reservoir having at least one pressure change term and at least one mass change term, wherein the block sparse matrix has at least one coefficient associated with the at least one pressure change term and at least one coefficient associated with the at least one mass change term.

26. The method of claim 24, wherein constructing the reduced transformed block sparse matrix comprises:

determining a column matrix norm for the column matrix;

skipping the singular value decomposition if the column matrix norm is equal to or less than a predetermined threshold value;

setting the singular values to zero if the column matrix norm is equal to or less than the predetermined threshold value; and

setting the third matrix to the identity matrix if the column matrix norm is equal to or less than the predetermined matrix threshold.

27. The method of claim 26, wherein determining the column matrix norm comprises calculating the Frobenius norm of the column matrix.