An Efficient Algebraic Multigrid Solver Strategy for Adaptive Implicit Methods in Oil-Reservoir Simulation
- Tanja Clees (Fraunhofer Institute for Algorithms and Scientific Computing) | Leonhard Ganzer (Clausthal University of Technology)
- Document ID
- Society of Petroleum Engineers
- SPE Journal
- Publication Date
- September 2010
- Document Type
- Journal Paper
- 670 - 681
- 2010. Society of Petroleum Engineers
- 5.5 Reservoir Simulation
- systems of partial differential equations, adaptive implicit methods (AIM), algebraic multigrid methods (AMG), automatic solver switching, constraint pressure residual (CPR)
- 0 in the last 30 days
- 502 since 2007
- Show more detail
- View rights & permissions
|SPE Member Price:||USD 12.00|
|SPE Non-Member Price:||USD 35.00|
We propose a new, efficient, adaptive algebraic multigrid (AMG) solver strategy for the discrete systems of partial-differential equations (PDEs) arising from structured or unstructured grid models in reservoir simulation. The proposed strategy has been particularly tailored to linear systems of equations arising in adaptive implicit methods (AIMs). The coarsening process of the AMG method designed automatically employs information on the physical structure of the models; as a smoother, an adaptive incomplete LU factorization with thresholding (ILUT) method is employed, taking care of an efficient solution of the hyperbolic parts while providing adequately smooth errors for the elliptic parts. To achieve a good compromise of high efficiency and robustness for a variety of problem classes--ranging from simple, small black-oil to challenging, large compositional models--an automatic, adaptive ILUT parameter and AMG solver switching strategy, a-SAMG, has been developed. Its efficiency is demonstrated for eight industrial benchmark cases by comparison against standard one-level and AMG solvers, including constraint pressure residual (CPR), as well as the pure one-level variant of the proposed new strategy. In addition, very promising results of first parallel runs are shown.
|File Size||788 KB||Number of Pages||12|
Amado, L.C.N., Ganzer, L.J., and Heinemann, Z.E. 1994. Finite VolumeDiscretization of the Fluid Flow Equations on General Perpendicular BisectionGrids. Paper presented at the Fifth International Forum on ReservoirSimulation, Muscat, Oman, 10-14 December.
Brand, C.W. and Ganzer, L. 1996. Iterative Solvers for Dynamically ImplicitReservoir Flow Equations on Irregular Grids. Paper presented at the 5thEuropean Conference on the Mathematics of Oil Recovery (ECMOR), Leoben,Austria, 3-6 September.
Cao, H., Tchelepi, H.A., Wallis, J., and Yardumian, H. 2005. Parallel Scalable UnstructuredCPR-Type Linear Solver for Reservoir Simulation. Paper SPE 96809 presentedat the SPE Annual Technical Conference and Exhibition, Dallas, 9-12 October.doi: 10.2118/96809-MS.
Clees, T. 2005. AMG Strategies for PDE Systems with Applications inIndustrial Semiconductor Simulation. PhD thesis, University of Cologne(Universität zu Köln), Cologne, Germany (November 2004).
Fung, L.S.K. and Dogru, A.H. 2008. Parallel Unstructured-Solver Methodsfor Simulation of Complex Giant Reservoirs. SPE J. 13(4): 440-446. SPE-106237-PA. doi: 10.2118/106237-PA.
Klie, H. 1996. Krylov-Secant Methods for Solving Large Scale Systems ofCoupled Nonlinear Parabolic Equations. PhD thesis TR96-28, Dept. ofComputational and Applied Mathematics, Rice University, Houston, Texas(September 1996).
Klie, H., Wheeler, M.F., Clees, T., and Stüben, K. 2007. Deflation AMG Solvers for HighlyIll-Conditioned Reservoir Simulation Problems. Paper SPE 105820 presentedat the SPE Reservoir Simulation Symposium, Houston, 28-30 February. doi:10.2118/105820-MS.
Krechel, A. and Stüben, K. 2005. SAMGp User's Manual, Release 21z (October2005). St. Augustin, Germany: Fraunhofer SCAI. http://www.scai.fraunhofer.de/samg.html.
Papadopoulos, A.T. and Tchelepi, H. 2003. Block Smoothed Aggregation AMGPreconditioning for Oil Reservoir Simulation Systems. Research Report NA-03/04,Numerical Analysis Group, Oxford University Computing Laboratory, Oxford, UK(September 2003).
Saad, Y. 1994. ILUT: Adual threshold incomplete ILU factorization. Numerical Linear Algebrawith Applications 1 (4): 387-402.doi:10.1002/nla.1680010405.
Saad, Y. 2003. Iterative Methods for Sparse Linear Systems, secondedition. Philadelphia, Pennsylvania: SIAM Publications.
Saad, Y. and Schultz, M.H. 1986. GMRES: A Generalized Minimal ResidualAlgorithm for Solving Nonsymmetric Linear Systems. SIAM J. Sci. Stat.Comp. 7 (3): 856-869. doi:10.1137/0907058.
Schenk, O. and Gärtner, K. 2006. On fast factorization pivoting methods forsymmetric indefinite systems. Electronic Transactions on Numerical Analysis(ETNA) 23: 158-179.
Schenk, O. and Gärtner, K. 2004. Solving unsymmetricsparse systems of linear equations with PARDISO. Future GenerationComputer Systems 20 (3): 475-487.doi:10.1016/j.future.2003.07.011.
Stüben, K. 2001a. A review of algebraicmultigrid. Journal of Computational and Applied Mathematics 128 (1-2): 281-309. doi:10.1016/S0377-0427(00)00516-1.
Stüben, K. 2001b. An Introduction to Algebraic Multigrid. InMultigrid, ed. U. Trottenberg, C.W. Ooosterlee, and A. Schüller,Appendix, 413-532. London: Academic Press.
Stüben, K. and Clees, T. 2005. SAMG User's Manual, Release 22c (June 2005).St. Augustin, Germany: Fraunhofer SCAI. http://www.scai.fraunhofer.de/samg.html.
Stüben, K., Clees, T., Klie, H., Lou, B., and Wheeler, M.F. 2007. Algebraic Multigrid Methods (AMG)for the Efficient Solution of Fully Implicit Formulations in ReservoirSimulation. Paper SPE 105832 presented at the SPE Reservoir SimulationSymposium, Houston, 28-30 February. doi: 10.2118/105832-MS.
Stüben, K., Delaney, P., and Chmakov, S. 2003. Algebraic Multigrid (AMG) forGround Water Flow and Oil Reservoir Simulation. Proc., MODFLOW and More 2003:Understanding through Modeling, International Ground Water Modeling Center(IGWMC), Golden, Colorado, USA, 17-19 September.
Thomas, G.W. and Thurnau, D.H. 1983. Reservoir Simulation Using anAdaptive Implicit Method. SPE J. 23 (5): 759-768.SPE-10120-PA. doi: 10.2118/10120-PA.
Trottenberg, U. and Clees, T. 2009. Multigrid Software for IndustrialApplications--From MG00 to SAMG. In 100 Volumes of 'Notes on Numerical FluidMechanics': 40 Years of Numerical Fluid Mechanics and Aeorodynamics inRetrospect, E.H. Hirschel and E. Krause ed., NNFM Vol. 100, 423. Berlin,Germany: Notes on Numerical Fluid Mechanics and Multidisciplinary Design,Springer-Verlag.
Trottenberg, U., Oosterlee, C.W., and Schüller, T. 2001. Multigrid.London: Academic Press.
van der Vorst, H.A. 1992. Bi-CGSTAB: A Fast and SmoothlyConverging Variant of Bi-CG for the Solution of Nonsymmetric LinearSystems. SIAM J. Sci. and Stat. Comput. 13 (2):631-644. doi:10.1137/0913035.
Vanek, P., Mandel, J., and Brezina, M. 1996. Algebraic multigrid by smoothedaggregation for second and fourth order problems. Computing 56 (3): 179-196. doi:10.1007/BF02238511.
Vinesome, P.K.W. 1976. Orthomin,an Iterative Method for Solving Sparse Banded Sets of Simultaneous LinearEquations. Paper SPE 5729 presented at the SPE Symposium on NumericalSimulation of Reservoir Performance, Los Angeles, 19-20 February. doi:10.2118/5729-MS.
Wallis, J.R., Kendall, R.P., and Little, T.E. 1985. Constrained Residual Acceleration ofConjugate Residual Methods. Paper SPE 13563 presented at the SPE Symposiumon Reservoir Simulation, Dallas, 10-13 February. doi: 10.2118/13563-MS.