A Variable-Switching Method for Mass-Variable-Based Reservoir Simulators
- Faruk O. Alpak (Shell International Exploration and Production Inc.) | Jeroen C. Vink (Shell Global Solutions International B.V.)
- Document ID
- Society of Petroleum Engineers
- SPE Journal
- Publication Date
- October 2018
- Document Type
- Journal Paper
- 1,469 - 1,495
- 2018.Society of Petroleum Engineers
- Newton-Raphson method, variable switching, natural variable, simulator formulation, mass variable
- 9 in the last 30 days
- 92 since 2007
- Show more detail
- View rights & permissions
|SPE Member Price:||USD 10.00|
|SPE Non-Member Price:||USD 30.00|
Field-scale simulations of complex processes often suffer from long simulation times. One of the main reasons is that the Newton-Raphson (NR) process used to solve each simulation timestep requires many iterations and small timestep sizes to converge. Because the selection of solution variables affects the nonlinearity of the equations, it is attractive to have a practical method to rapidly explore the use of alternative primary variables in general-purpose reservoir simulators.
Many reservoir simulators use pressure, saturations, and temperature in each gridblock as the primary solution variables, which are referred to as natural variables. There is also a class of reservoir simulators that uses pressure, total component masses (or moles), and internal energy in each gridblock as primary variables. These simulators are referred to as mass-variable-based reservoir simulators. For a given choice of primary variables, most simulators have dedicated, highly optimized procedures to compute the required derivatives and chain rules required to build the Jacobian matrix. Hence, it is usually not possible to switch between mass and natural variables. In this work, however, we establish a link at the numerical-solution level between natural- and mass-variable formulations and design a novel (nonlinear) block-local method that transforms mass-variable shifts (computed at each NR iteration) into equivalent natural-variable shifts.
We demonstrate on a number of simulation models of varying complexity that by use of the proposed approach, a mass-variable-based flow simulator can still effectively use natural variables, where the change of variables can be made locally per gridblock. Results indicate that in some models the total number of NR iterations, linear-solver (LS) iterations, and timestep-size cuts (caused by the nonconvergence of the NR procedure, also known as backups) are reduced when using natural variables instead of mass variables. However, the improvement is relatively modest and not generally observed. Our findings also signify that depending on the specific characteristics of the simulation problem at hand, mass-variable-based simulators may perform comparably or outperform natural-variable-based simulators.
The proposed variable-switching method can be used effectively to evaluate the effect of using different primary solution variables on problem nonlinearity and solver efficiency. With this method, the effect of interchanging primary solution variables on problem nonlinearity can be rapidly evaluated.
|File Size||2 MB||Number of Pages||27|
Aavatsmark, I., Barkve, T., and Mannseth, T. 1998. Control-Volume Discretization Methods for 3D Quadrilateral Grids in Inhomogeneous, Anisotropic Reservoirs. SPE J. 3 (2): 146–154. SPE-38000-PA. https://doi.org/10.2118/38000-PA.
Alpak, F. O., Vink, J. C., Gao, G. et al. 2013. Techniques for Effective Simulation, Optimization, and Uncertainty Quantification of the In-Situ Upgrading Process. J. Unconven. Oil Gas Resour. 3–4 (December): 1–14. https://doi.org/10.1016/j.juogr.2013.09.001.
Appleyard, J. R., Cheshire, I. M., and Pollard, R. K. 1981. Special Techniques for Fully-Implicit Simulators. Technical Report HL81/2039, Atomic Energy Research Establishment, Harwell, Oxfordshire, UK.
Aziz, K. and Wong, T. W. 1989. Considerations in the Development of Multipurpose Reservoir Simulation Models. Oral presentation given at the 2nd International Forum on Reservoir Simulation, Alpbach, Austria, 4–8 September.
Aziz, K., Ramesh, A. B., and Woo, P. T. 1987. Fourth SPE Comparative Solution Project: Comparison of Steam Injection Simulators. J Pet Technol 39 (12): 1576–1584. SPE-13510-PA. https://doi.org/10.2118/13510-PA.
Behie, A. 1985. Comparison of Nested Factorization, Constrained Pressure Residual, and Incomplete Factorization Preconditionings. Presented at the SPE Reservoir Simulation Symposium, Dallas, 10–13 February. SPE-13531-MS. https://doi.org/10.2118/13531-MS.
Cao, H. 2002. Development of Techniques for General Purpose Simulators. PhD dissertation, Stanford University, Stanford, California.
Cao, H., Crumpton, P. I., and Schrader, M. L. 2009. Efficient General Formulation Approach for Modeling Complex Physics. Presented at the SPE Reservoir Simulation Symposium, The Woodlands, Texas, 2–4 February. SPE-119165-MS. https://doi.org/10.2118/119165-MS.
Cao, H., Tchelepi, H. A., Wallis, J. R. et al. 2005. Parallel Scalable Unstructured CPR-Type Linear Solver for Reservoir Simulation. Presented at the SPE Annual Technical Conference and Exhibition, Dallas, 9–12 October. SPE-96809-MS. https://doi.org/10.2118/96809-MS.
Cao, H., Zaydullin, R., and Obi, E. 2017. Nonlinear Convergence for Near-Miscible Problem: A Mystery Unveiled for Natural Variable Simulator. Presented at the SPE Reservoir Simulation Conference, Montgomery, Texas, 20–22 February. SPE-182633-MS. https://doi.org/10.2118/182633-MS.
Chien, M. C. H., Lee, S. T., and Chen, W. H. 1985. A New Fully Implicit Compositional Simulator. Presented at the SPE Reservoir Simulation Symposium, Dallas, 10–13 February. SPE-13385-MS. https://doi.org/10.2118/13385-MS.
Chien, M. C. H., Yardumian, H. E., Chung, E. Y. et al. 1989. The Formulation of a Thermal Simulation Model in a Vectorized, General Purpose Reservoir Simulator. Presented at the SPE Symposium on Reservoir Simulation, Houston, 6–8 February. SPE-18418-MS. https://doi.org/10.2118/18418-MS.
Coats, K. H. 1980. An Equation of State Compositional Model. SPE J. 20 (5): 363–376. SPE-8284-PA. https://doi.org/10.2118/8284-PA.
Concus, P., Golub, G. H., and Meurant, G. 1985. Block Preconditioning for the Conjugate Gradient Method. SIAM J. Sci. Stat. Comput. 6 (1): 120–252. https://doi.org/10.1137/0906018.
Collins, D. A., Nghiem, L. X., Li, Y.-K. et al. 1992. An Efficient Approach to Adaptive-Implicit Compositional Simulation With an Equation of State. SPE Res Eval & Eng 7 (2): 259–264. SPE-15133-PA. https://doi.org/10.2118/15133-PA.
Courant, R., Friedrichs, K., and Lewy, H. 1928. Über die Partiellen Differenzengleichungen der Mathematischen Physik (in German). Mathematische Annalen 100 (1): 32–74.
Das, S. K. 2007. Application of Thermal Processes in Heavy Oil Carbonate Reservoirs. Presented at the SPE Middle East Oil and Gas Show and Conference, Manama, Bahrain, 11–14 March. SPE-105392-MS. https://doi.org/10.2118/105392-MS.
Dreher, K. D., Kenyon, D. E., and Iwere, F. O. 1986. Heat Flow During Steam Injection Into a Fractured Carbonate Reservoir. Presented at the SPE Enhanced Oil Recovery Symposium, Tulsa, 20–23 April. SPE-14902-MS. https://doi.org/10.2118/14902-MS.
Fussell, L. T. and Fussell, D. D. 1979. An Iterative Technique for Compositional Reservoir Models. SPE J. 19 (4): 211–220. SPE-6891-PA. https://doi.org/10.2118/6891-PA.
Haugen, K. B. and Beckner, B. L. 2015. A General Flow Equation Framework. Presented at the SPE Reservoir Simulation Symposium, Houston, 23–25 February. SPE-173266-MS. https://doi.org/10.2118/173266-MS.
Killough, J. E. and Kossack, C. A. 1987. Fifth Comparative Solution Project: Evaluation of Miscible Flood Simulators. Presented at the SPE Reservoir Simulation Symposium, San Antonio, Texas, 1–4 February. SPE-16000-MS. https://doi.org/10.2118/16000-MS.
Maliska, C. R., da Silva, A. F. C., Czesnat, A. O. et al. 1997. Three-Dimensional Multiphase Flow Simulation in Petroleum Reservoirs Using the Mass Fractions as Dependent Variables. Presented at the Latin American and Caribbean Petroleum Engineering Conference, Rio de Janeiro, 30 August–3 September. SPE-39067-MS. https://doi.org/10.2118/39067-MS.
Por, G. J., Boerrigter, P., Maas, J. G. et al. 1989. A Fractured Reservoir Simulator Capable of Modeling Block-Block Interaction. Presented at the SPE Annual Technical Conference and Exhibition, San Antonio, Texas, 8–11 October. SPE-19807-MS. https://doi.org/10.2118/19807-MS.
Regtien, J. M. M., Por, G. J. A., van Stiphout, M. T. et al. 1995. Interactive Reservoir Simulation. Presented at the SPE Annual Technical Conference and Exhibition, San Antonio, Texas, 12–15 February. SPE-29146-MS. https://doi.org/10.2118/29146-MS.
Saad, Y. 2003. Iterative Methods for Sparse Linear Systems, second edition. Philadelphia, Pennsylvania: Society for Industrial and Applied Mathematics.
Saad, Y. and Schultz, M. H. 1986. GMRES: A Generalized Minimal Residual Algorithm for Solving Nonsymmetric Linear Systems. SIAM J. Sci. Stat. Comput. 7 (3): 856–869. https://doi.org/10.1137/0907058.
Sahuquet, B. C. and Ferrier, J. J. 1982. Steam-Drive Pilot in a Fractured Carbonated Reservoir: Lacq Superieur Field. J Pet Technol 34 (4): 873–880. SPE-9453-PA. https://doi.org/10.2118/9453-PA.
Schlumberger. 2003. Eclipse Technical Description. Houston: Schlumberger.
Schmall, L., Varavei, A., and Sepehrnoori, K. 2013. A Comparison of Various Formulations for Compositional Reservoir Simulation. Presented at the SPE Reservoir Simulation Symposium, The Woodlands, Texas, 18–20 February. SPE-163630-MS. https://doi.org/10.2118/163630-MS.
Shahin, G. T., Moosa, R., Al-Kharusi, B. S. et al. 2006. The Physics of Steam Injection in Fractured Carbonate Reservoirs: Engineering Development Options That Minimize Risk. Presented at the SPE Annual Technical Conference and Exhibition, San Antonio, Texas, 24–27 September. SPE-102186-MS. https://doi.org/10.2118/102186-MS.
Shi, X., Chang, Y.-B., Muller, M. et al. 2009. A General Unstructured Grid, Parallel, Fully Implicit Thermal Simulator and Its Application for Large Scale Thermal Models. Presented at the SPE Reservoir Simulation Symposium, The Woodlands, Texas, 2–4 February. SPE-119172-MS. https://doi.org/10.2118/119172-MS.
Tismenetsky, M. 1991. A Note on a Block Preconditioner. Appl. Math. Lett. 4 (2): 75–79. https://doi.org/10.1016/0893-9659(91)90173-S.
van der Vorst, H. A. 1992. Bi-CGSTAB: A Fast and Smoothly Converging Variant of Bi-CG for the Solution of Nonsymmetric Linear Systems. SIAM J. Sci. Stat. Comput. 13 (2): 631–644. https://doi.org/10.1137/0913035.
van Wunnik, J. N. M. and Wit, K. 1992. Improvement of Gravity Drainage by Steam Injection Into a Fractured Reservoir: An Analytical Evaluation. SPE Res Eng 7 (1): 59–66. SPE-20251-PA. https://doi.org/10.2118/20251-PA.
Vinsome, P. K.W. 1976. Orthomin, an Iterative Method for Solving Sparse Banded Sets of Simultaneous Linear Equations. Presented at the SPE Symposium on Numerical Simulation of Reservoir Performance, Los Angeles, 19–20 February. SPE-5729-MS. https://doi.org/10.2118/5729-MS.
Voskov, D. V. and Tchelepi, H. A. 2011. Compositional Nonlinear Solver Based on Trust Regions of the Flux Function Along Key Tie-Lines. Presented at the SPE Reservoir Simulation Symposium, The Woodlands, Texas, 21–23 February. SPE-141743-MS. https://doi.org/10.2118/141743-MS.
Voskov, D. V. and Tchelepi, H. A. 2012. Comparison of Nonlinear Formulations for Two-Phase Multi-Component EoS Based Simulation. J. Pet. Sci. Eng. 82–83 (February–March): 101–111. https://doi.org/10.1016/j.petrol.2011.10.012.
Voskov, D. V., Tchelepi, H. A., and Younis, R. M. 2009. General Nonlinear Solution Strategies for Multiphase Multicomponent EoS Based Simulation. Presented at the SPE Reservoir Simulation Symposium, The Woodlands, Texas, 2–4 February. SPE-118996-MS. https://doi.org/10.2118/118996-MS.
Wang, P., Yotov, I., Wheeler, M. et al. 1997. A New Generation EoS Compositional Reservoir Simulator: Part I—Formulation and Discretization. Presented at the SPE Reservoir Simulation Symposium, Dallas, 8–11 June. SPE-37979-MS. https://doi.org/10.2118/37979-MS.
Yanosik, J. L. and McCracken, T. A. 1979. A Nine-Point, Finite-Difference Reservoir Simulator for Realistic Prediction of Adverse Mobility Ratio Displacements. SPE J. 19 (4): 253–262. SPE-5734-PA. https://doi.org/10.2118/5734-PA.
Young, L. C. and Stephenson, R. E. 1983. A Generalized Compositional Approach for Reservoir Simulation. SPE J. 23 (5): 727–742. SPE-10516-PA. https://doi.org/10.2118/10516-PA.
Zaydullin, R., Voskov, D., and Tchelepi, H. A. 2012. Nonlinear Formulation Based on an Equation-of-State Free Method for Compositional Flow Simulation. SPE J. 18 (2): 264–273. SPE-146989-PA. https://doi.org/10.2118/146989-PA.