Enhanced-Velocity-Multiblock Method for Coupled Flow and Reactive-Species Transport Through Porous Media: Applications to Bioremediation and Carbon Sequestration
- Sunil G. Thomas (Chevron) | Mary F. Wheeler (University of Texas at Austin)
- Document ID
- Society of Petroleum Engineers
- SPE Journal
- Publication Date
- September 2012
- Document Type
- Journal Paper
- 794 - 804
- 2012. Society of Petroleum Engineers
- 5.10.1 CO2 Capture and Sequestration, 5.5 Reservoir Simulation, 4.1.2 Separation and Treating, 5.3.2 Multiphase Flow
- 1 in the last 30 days
- 262 since 2007
- Show more detail
- View rights & permissions
|SPE Member Price:||USD 10.00|
|SPE Non-Member Price:||USD 30.00|
This paper presents a multiblock-discretization method--the enhanced-velocity mixed-?nite-element method (EVMFEM) (Wheeler et al. 2002--for coupled multiphase flow and reactive-species-transport modeling in porous-media applications. The method provides local mass balance and a continuous approximation of fluxes across interfaces of elements and subdomains. It can treat nonmatching grids, allowing for a flexible choice of grid refinements. Further, by distributing the blocks among processors such that each block has approximately the same number of elements, this method can be implemented efficiently in parallel, thereby offering further reductions in computational cost.
The paper also presents recent application of EVMFEM to challenging problems such as compositional flow simulations of CO2 sequestration. Tests with EVMFEM suggest that it is advantageous to apply grid refinements around wells and to areas in which dynamics of chemical-species concentration is highest. Allowing for variable grid refinements greatly reduces the simulation cost, while preserving overall accuracy of the solution. For completeness, a few significant analytic results on convergence of the method are stated and referenced, omitting proof.
This work is significant in advancing the discretization and application of EVMFEMs in reservoir-simulation development. Problems such as transport of chemical species in multiphase flow and CO2 sequestration have begun to assume significant importance in decisions regarding the preservation of our environment and in the safe and reliable means of delivering energy. This paper offers useful methods and some innovative future directions to address the huge computational costs involved in solving such complex problems.
|File Size||6 MB||Number of Pages||11|
Arbogast, T., Cowsar, L.C., Wheeler, M.F., and Yotov, I. 2000. Mixed FiniteElement Methods on Nonmatching Multiblock Grids. SIAM J. Numer. Anal. 37 (4): 1295-1315. http://dx.doi.org/10.1137/s0036142996308447.
Arbogast, T., Pencheva, G., Wheeler, M.F., and Yotov, I. 2007. A MultiscaleMortar Mixed Finite Element Method. Multiscale Modeling & Simulation 6 (1): 319-346. http://dx.doi.org/10.1137/060662587.
Arbogast, T., Wheeler, M.F., and Yotov, I. 1997. Mixed Finite Elements forElliptic Problems with Tensor Coefficients as Cell-Centered Finite Differences.SIAM J. Numer. Anal. 34 (2): 828-852. http://dx.doi.org/10.1137/s0036142994262585.
Bell, J.B., Dawson, C.N., and Shubin, G.R. 1988. An unsplit, higher ordergodunov method for scalar conservation laws in multiple dimensions. J.Comput. Phys. 74 (1): 1-24. http://dx.doi.org/10.1016/0021-9991(88)90065-4.
Brezzi, F. and Fortin, M. 1991. Mixed and Hybrid Finite ElementMethods, Vol. 15. New York: Springer Series in Computational Mathematics,Springer-Verlag.
Brooks, R.H. and Corey, A.T. 1964. Hydraulic properties of porous media.Hydrology Paper No. 3, Colorado State University, Fort Collins, Colorado,22-27.
Chavent, G., Cohen, G., Jaffre, J., Dupuy, M., and Ribera, I. 1984.Simulation of Two-Dimensional Waterflooding By Using Mixed Finite Elements.SPE J. 24 (4): 382-390. SPE-10502-PA. http://dx.doi.org/10.2118/10502-PA.
Chen, Z., Ma, Y., and Chen, G. 2007. A sequential numerical chemicalcompositional simulator. Transport Porous Media 68 (3):389-411. http://dx.doi.org/10.1007/s11242-006-9050-y.
Christie, M.A. and Blunt, M.J. 2001. Tenth SPE Comparative Solution Project:A Comparison of Upscaling Techniques. SPE Res Eval & Eng 4 (4): 308-317. SPE-72469-PA. http://dx.doi.org/10.2118/72469-PA.
Class, H., Ebigbo, A., Helmig, R. et al. 2009. A benchmark study on problemsrelated to CO2 storage in geologic formations. Comput.Geosci. 13 (4): 409-434. http://dx.doi.org/10.1007/s10596-009-9146-x.
Dawson, C.R. 1993. Godunov-Mixed Methods for Advection-Diffusion Equationsin Multidimensions. SIAM J. Numer. Anal. 30 (5): 1315-1332.http://dx.doi.org/10.1137/0730068.
Delshad, M., Thomas, S.G., and Wheeler, M.F. 2011. Parallel NumericalReservoir Simulations of Nonisothermal Compositional Flow and Chemistry. SPEJ. 16 (2): 239-248. SPE-118847-PA. http://dx.doi.org/10.2118/118847-PA.
Durlofsky, L.J. and Chien, M.C.H. 1993. Development of a MixedFinite-Element-Based Compositional Reservoir Simulator. Paper SPE 25253presented at the SPE Reservoir Simulation Symposium, New Orleans, 28 February-3March. http://dx.doi.org/10.2118/25253-MS.
Ewing, R.E., Russell, T.F., and Wheeler, M.F. 1983. Simulation of MiscibleDisplacement Using Mixed Methods and a Modified Method of Characteristics.Paper SPE 12241 presented at the SPE Reservoir Simulation Symposium, SanFrancisco, 15-18 November. http://dx.doi.org/10.2118/12241-MS.
Eymard, R., Gallouët, T., and Herbin, R. 2000. Finite volume methods. InHandbook of Numerical Analysis VII, ed. P.G. Ciarlet and J.L. Lions,713-1020. Amsterdam, The Netherlands: Elsevier B.V.
Firoozabadi, A. and Hoteit, H. 2006. Compositional Modeling by the CombinedDiscontinuous Galerkin and Mixed Methods. SPE J. 11 (1):19-34. SPE-90276-PA. http://dx.doi.org/10.2118/90276-PA.
Nedelec, J.C. 1980. Mixed finite elements in R3. NumerischeMathematik 35 (3): 315-341. http://dx.doi.org/10.1007/bf01396415.
Peng, D.-Y. and Robinson, D.B. 1976. A New Two-Constant Equation of State.Industrial & Engineering Chemistry Fundamentals 15 (1):59-64. http://dx.doi.org/10.1021/i160057a011.
Peszynska, M. and Sun, S. 2001. Multiphase reactive transport module(TRCHEM) in IPARS. Technical Report ICES 01-32, Texas Institute forComputational and Applied Mathematics, The University of Texas at Austin,Austin, Texas (October 2001).
Raviart, P. and Thomas, J. 1977. A mixed finite element method for 2nd orderelliptic problems. In Mathematical Aspects of Finite Element Methods,ed. I. Galligani and E. Magenes, Vol. 606, 292-315. Berlin: Lecture Notes inMathematics, Springer Verlag.
Thomas, S. and Wheeler, M. 2011. Enhanced velocity mixed finite elementmethods for modeling coupled flow and transport on non-matching multiblockgrids. Comput. Geosci. 15 (4): 605-625. http://dx.doi.org/10.1007/s10596-011-9227-5.
Thomas, S.G. 2009. On some problems in the simulation of flow andtransport through porous media. PhD dissertation, University of Texas atAustin, Austin, Texas (August 2009).
Universität Stuttgart. 2008. Workshop on Numerical Models for Carbon DioxideStorage in Geological Formations, Stuttgart, Germany, 2-4 April, http://www.iws.uni-stuttgart.de/co2-workshop.
Wheeler, J.A. 2001. Integrated parallel and accurate reservoir simulator.Technical Report TICAM01-25, Institute of Computational Engineering andSciences, The University of Texas at Austin, Austin, Texas. [author added]
Wheeler, J.A., Wheeler, M.F., and Yotov, I. 2002. Enhanced Velocity MixedFinite Element Methods for Flow in Multiblock Domains. Comput. Geosci. 6 (3): 315-332. http://dx.doi.org/10.1023/a:1021270509932.