Adaptation of the CPR Preconditioner for Efficient Solution of the Adjoint Equation
- Choongyong Han (Chevron Energy Technology Company) | John Wallis (Wallis Consulting Incorporated) | Pallav Sarma (Chevron Energy Technology Company) | Gary Li (Chevron Energy Technology Company) | Mark Schrader (Chevron Energy Technology Company) | Wen H Chen (Chevron Energy Technology Company)
- Document ID
- Society of Petroleum Engineers
- SPE Journal
- Publication Date
- January 2013
- Document Type
- Journal Paper
- 207 - 213
- 2013. Society of Petroleum Engineers
- 5.5 Reservoir Simulation, 5.5.8 History Matching
- 1 in the last 30 days
- 336 since 2007
- Show more detail
- View rights & permissions
|SPE Member Price:||USD 10.00|
|SPE Non-Member Price:||USD 30.00|
It is well known that the adjoint approach is the most efficient approach for gradient calculation, and it can be used with gradient-based optimization techniques to solve various optimization problems, such as the production-optimization problem and the history-matching problem. The adjoint equation to be solved in the approach is a linear equation formed with the "transpose" of the Jacobian matrix from a fully implicit reservoir simulator. For a large and/or complex reservoir model, generalized preconditioners often prove impractical for solving the adjoint equation. Preconditioners specialized for reservoir simulation, such as constrained pressure residual (CPR), exploit properties of the Jacobian matrix to accelerate convergence, so they cannot be applied directly to the adjoint equation. To overcome this challenge, we have developed a new two-stage preconditioner for efficient solution of the adjoint equation by adaptation of the CPR preconditioner (named CPRA: CPR preconditioner for adjoint equation).
The CPRA preconditioner has been coupled with an algebraic multigrid (AMG) linear solver and implemented in Chevron's extended applications reservoir simulator (CHEARS(R)). The AMG solver is well known for its outstanding capability to solve the pressure equation of complex reservoir models; solving the linear system with the "transpose" of the pressure matrix is one of the two stages of construction of the CPRA preconditioner.
Through test cases, we have confirmed that the CPRA/AMG solver with generalized minimal residual (GMRES) acceleration solves the adjoint equation very efficiently with a reasonable number of linear-solver iterations. Adjoint simulations to calculate the gradients with the CPRA/AMG solver take approximately the same amount of time (at most) as do the corresponding CPR/AMG forward simulations. Accuracy of the solutions has also been confirmed by verifying the gradients against solutions with a direct solver. A production-optimization case study for a real field using the CPRA/AMG solver has further validated its accuracy, efficiency, and the capability to perform long-term optimization for large, complex reservoir models at low computational cost.
|File Size||664 KB||Number of Pages||7|
Appleyard, J.R. 1983. Nested Factorization. Presented at the SPE ReservoirSimulation Symposium, San Francisco, California, USA, 15-18 November.SPE-12264-MS. http://dx.doi.org/10.2118/12264-MS.
Brouwer, D.R. and Jansen, J.-D. 2004. Dynamic Optmization of Water FloodingWith Smart Wells Using Optimal Control Theory. SPE J. 9(4): 391-402. SPE-78278-PA. http://dx.doi.org/10.2118/78278-PA.
Cao, H., Tchelepi, H.A., Wallis, J. et al. 2005. Parallel ScalableUnstructured CPR-Type Linear Solver for Reservoir Simulation. Presented at theSPE Annual Technical Conference and Exhibition, Dallas, 9-12 October.SPE-96809-MS. http://dx.doi.org/10.2118/96809-MS.
Chaudhri, M.M., Phale, H.A., and Liu, N. 2009. An Improved Approach ForEnsemble Based Production Optimization: Application To A Field Case. Presentedat the EUROPEC/EAGE Conference and Exhibition, Amsterdam, 8-11 June.SPE-121307-MS. http://dx.doi.org/10.2118/121307-MS.
Chen, Y., Oliver, D.S., and Zhang, D. 2009. Efficient Ensemble-BasedClosed-Loop Production Optimization. SPE J. 14 (4):634-645. SPE-112873-PA. http://dx.doi.org/10.2118/112873-PA.
Golub, G.H., Stoll, M., and Wathen, A. 2008. Approximation of the ScatteringAmplitude and Linear Systems. Electron. Trans. Numer. Anal. 31(2008): 178-203.
Li, R., Reynolds, A.C., and Oliver, D.S. 2003. History Matching ofThree-Phase Flow Production Data. SPE J. 8 (4): 328-340.SPE-87336-PA. http://dx.doi.org/10.2118/87336-PA.
Rao, S.S. 2009. Engineering Optimization: Theory and Practice, fourthedition. New York: John Wiley & Sons.
Saad, Y. and Schultz, M.H. 1986. GMRES: A Generalized Minimal ResidualAlgorithm for Solving Nonsymmetric Linear Systems. SIAM J. Sci. and Stat.Comput. 7 (3): 856-869. http://dx.doi.org/10.1137/0907058.
Sarma, P., Aziz, K., and Durlofsky, L.J. 2005. Implementation of AdjointSolution for Optimal Control of Smart Wells. Presented at the SPE ReservoirSimulation Symposium, The Woodlands, Texas, USA, 31 January-2 February.SPE-92864-MS. http://dx.doi.org/10.2118/92864-MS.
Sarma, P., Chen, W.H., Durlofsky, L.J. et al. 2008. Production OptimizationWith Adjoint Models Under Nonlinear Control-State Path Inequality Constraints.SPE Res Eval & Eng 11 (2): 326-339. SPE-99959-PA. http://dx.doi.org/10.2118/99959-PA.
Schneider, R. and Jimack, P.K. 2005. Efficient preconditioning of thediscrete adjoint equations for the incompressible Navier-Stokes equations.Int. J. Numer. Methods Fluids 47 (10-11): 1277-1283. http://dx.doi.org/10.1002/fld.874.
Stengel, R.F. 1994. Optimal Control and Estimation. Mineola, NewYork: Dover Books on Advanced Mathematics, Dover Publications.
Wallis, J., Tchelepi, H.A., and Cao, H. 2009. Apparatus, method and systemfor improved reservoir simulation using a multiplicative overlapping Schwarzpreconditioning for adaptive implicit linear systems. US Patent No.7,516,056.
Wallis, J.R. 1983. Incomplete Gaussian Elimination as a Preconditioning forGeneralized Conjugate Gradient Acceleration. Presented at the SPE ReservoirSimulation Symposium, San Francisco, California, USA, 15-18 November.SPE-12265-MS. http://dx.doi.org/10.2118/12265-MS.
Wallis, J.R., Kendall, R.P., and Little, T.E. 1985. Constrained ResidualAcceleration of Conjugate Residual Methods. Presented at the SPE ReservoirSimulation Symposium, Dallas, 10-13 February. SPE-13536-MS. http://dx.doi.org/10.2118/13536-MS.