Efficient block preconditioners for the coupled equations of pressure and deformation in highly discontinuous media

SUMMARY Large-scale simulations of flow in deformable porous media require efficient iterative methods for solving the involved systems of linear algebraic equations. Construction of efficient iterative methods is particularly challenging in problems with large jumps in material properties, which is often the case in geological applications, such as basin evolution at regional scales. The success of iterative methods for this type of problems depends strongly on finding effective preconditioners. This paper investigates how the block-structured matrix system arising from single-phase flow in elastic porous media should be preconditioned, in particular for highly discontinuous permeability and significant jumps in elastic properties. The most promising preconditioner combines algebraic multigrid with a Schur complement-based exact block decomposition. The paper compares numerous block preconditioners with the aim of providing guidelines on how to formulate efficient preconditioners. Copyright © 2010 John Wiley & Sons, Ltd.

[1]  K. Toh,et al.  An efficient diagonal preconditioner for finite element solution of Biot's consolidation equations , 2002 .

[2]  C. Neuzil.  How permeable are clays and shales , 1994 .

[3]  K. Toh,et al.  Block preconditioners for symmetric indefinite linear systems , 2004 .

[4]  Ray S. Tuminaro,et al.  Parallel Smoothed Aggregation Multigrid : Aggregation Strategies on Massively Parallel Machines , 2000, ACM/IEEE SC 2000 Conference (SC'00).

[5]  M. Hestenes,et al.  Methods of conjugate gradients for solving linear systems , 1952 .

[6]  Mark F. Adams,et al.  Evaluation of three unstructured multigrid methods on 3D finite element problems in solid mechanics , 2000 .

[7]  C. E. Neuzil,et al.  Hydromechanical coupling in geologic processes , 2003 .

[8]  Marte Gutierrez,et al.  Coupling of Fluid Flow and Deformation in Underground Formations , 2002 .

[9]  William Pao,et al.  Fully coupled modeling of seabed subsidence and reservoir compaction of North Sea oil fields , 2003 .

[10]  J. Piau,et al.  Contribution of poroelasticity to reservoir engineering: Lab experiments, application to core decompression and implication in HP-HT reservoirs depletion , 1994 .

[11]  J. Shewchuk.  An Introduction to the Conjugate Gradient Method Without the Agonizing Pain , 1994 .

[12]  Aslak Tveito,et al.  Systems of PDEs and Block Preconditioning , 2003 .

[13]  M. Benzi.  Preconditioning techniques for large linear systems: a survey , 2002 .

[14]  Bjørn Fredrik Nielsen,et al.  Finite Element Discretizations of Elliptic Problems in the Presence of Arbitrarily Small Ellipticity: An Error Analysis , 1999 .

[15]  Ulrike Meier Yang,et al.  Parallel Algebraic Multigrid Methods — High Performance Preconditioners , 2006 .

[16]  Takumi Washio,et al.  Ordering strategies and related techniques to overcome the trade-off between parallelism and convergence in incomplete factorizations , 1999, Parallel Comput..

[17]  J. Pasciak,et al.  A preconditioning technique for indefinite systems resulting from mixed approximations of elliptic problems , 1988 .

[18]  Tamara G. Kolda,et al.  An overview of the Trilinos project , 2005, TOMS.

[19]  A. George,et al.  On the Complexity of Sparse $QR$ and $LU$ Factorization of Finite-Element Matrices , 1988 .

[20]  Thomas A. Manteuffel,et al.  Adaptive Smoothed Aggregation (αSA) , 2004, SIAM J. Sci. Comput..

[21]  M. Biot.  General Theory of Three‐Dimensional Consolidation , 1941 .

[22]  J. W. Ruge,et al.  4. Algebraic Multigrid , 1987 .

[23]  Mark F. Adams,et al.  Ultrascalable Implicit Finite Element Analyses in Solid Mechanics with over a Half a Billion Degrees of Freedom , 2004, Proceedings of the ACM/IEEE SC2004 Conference.

[24]  Henk A. van der Vorst,et al.  Bi-CGSTAB: A Fast and Smoothly Converging Variant of Bi-CG for the Solution of Nonsymmetric Linear Systems , 1992, SIAM J. Sci. Comput..

[25]  Mary F. Wheeler,et al.  Overcoming the problem of locking in linear elasticity and poroelasticity: an heuristic approach , 2009 .