Adaptive Multiscale Finite-Volume Framework for Reservoir Simulation
- Hamdi A. Tchelepi (Stanford University) | Patrick Jenny (Institute of Fluid Dynamics, ETH) | Seong Hee Lee (Chevron ETC) | Christian Wolfsteiner (Chevron ETC)
- Document ID
- Society of Petroleum Engineers
- SPE Journal
- Publication Date
- June 2007
- Document Type
- Journal Paper
- 188 - 195
- 2007. Society of Petroleum Engineers
- 5.1 Reservoir Characterisation, 5.3.2 Multiphase Flow, 4.1.2 Separation and Treating, 5.3.1 Flow in Porous Media, 5.5 Reservoir Simulation, 2.2.2 Perforating, 5.5.3 Scaling Methods, 4.3.4 Scale, 4.6 Natural Gas, 5.6.5 Tracers
- 2 in the last 30 days
- 707 since 2007
- Show more detail
- View rights & permissions
|SPE Member Price:||USD 10.00|
|SPE Non-Member Price:||USD 30.00|
A multiscale finite-volume (MSFV) framework for reservoir simulation is described. This adaptive MSFV formulation is locally conservative and yields accurate results of both flow and transport in large-scale highly heterogeneous reservoir models. IMPES and sequential implicit formulations are described. The algorithms are sensitive to the specific characteristics of flow (i.e., pressure and total velocity) and transport (i.e., saturation). To compute the fine-scale flow field, two sets of basis functions - dual and primal - are constructed. The dual basis functions, which are associated with the dual coarse grid, are used to calculate the coarse scale transmissibilities. The fine-scale pressure field is computed from the coarse grid pressure via superposition of the dual basis functions. Having a locally conservative fine scale velocity field is essential for accurate solution of the saturation equations (i.e., transport). The primal basis functions, which are associated with the primal coarse grid, are constructed for that purpose. The dual basis functions serve as boundary conditions to the primal basis functions. To resolve the fine-scale flow field in and around wells, a special well basis function is devised. As with the other basis functions, we ensure that the support for the well basis is local.
Our MSFV framework is designed for adaptive computation of both flow and transport in the course of a simulation run. Adaptive computation of the flow field is based on the time change of the total mobility field, which triggers the selective updating of basis functions. The key to achieving scalable (efficient for large problems) adaptive computation of flow and transport is the use of high fidelity basis functions with local support. We demonstrate the robustness and computational efficiency of the MSFV simulator using a variety of large heterogeneous reservoir models, including the SPE 10 comparative solution problem.
|File Size||2 MB||Number of Pages||8|
Arbogast, T. 2002. Implementation of a locally conservative numericalsubgrid upscaling scheme for two phase Darcy flow. ComputationalGeosciences 6: 453-481.
Arbogast, T. and Bryant, S.L. 2002. A Two-Scale Numerical SubgridTechnique for Waterflood Simulations. SPEJ 7: 446-457.SPE-81909-PA. DOI: 10.2118/81909-PA.
Aziz, K. and Settari, A. 1979. Petroleum Reservoir Simulation.Applied Science Publishers, London, England.
Chen, Z. and Hou, T.Y. 2003. A Mixed Multiscale Finite Element Method forElliptic Problems with Oscillating Coefficients. Math. Comput.72: 541-576.
Chien, M.C.H., Lee, S.T., and Chen, W.H. 1985. A New Fully Implicit CompositionalSimulator. Paper SPE 13385 presented at the SPE Symposium on ReservoirSimulation, Dallas. DOI: 10.2118/13385-MS.
Christie, M.A. and Blunt, M.J. 2001. Tenth SPE comparative solutionproject: A comparison of upscaling techniques. SPEREE 4:308-317. SPE-72469-PA. DOI: 10.2118/72469-PA.
Coats, K.H. 2000. A Note onIMPES and Some IMPES-Based Simulation Models. SPEJ 5:245-251. SPE-65092-PA. DOI: 10.2118/65092-PA.
Coats, K.H., George, W.D., Chu, C., and Marcum, B.E. 1974. Three-Dimensional Simulation ofSteamflooding. SPEJ 14 (6): 573-592. SPE-4500-PA. DOI:10.2118/4500-PA.
Dagan, G. 1989. Flow and Transport in Porous Formations.Springer-Verlag, New York, NY.
Gautier, Y., Blunt, M.J., and Christie, M.A. 1999. Nested Gridding andStreamline-Based Simulation for Fast Reservoir Performance Prediction.Computational Geosciences 3: 295-320.
Hoeksema, R.J. and Kitanidis, P. 1985. Analysis of the spatial structureproperties of selected aquifers. Water Resour. Res. 21:563-572.
Hou, T. and Wu, X.H. 1997. A Multiscale Finite Element Method for EllipticProblems in Composite Materials and Porous Media. J. Comp. Phys.134: 169-189.
Jenny, P., Lee, S.H., and Tchelepi, H.A. 2003. Multi-Scale Finite-VolumeMethod for Elliptic Problems in Subsurface Flow Simulation. J. Comp.Phys. 187: 47-67.
Lee, S.H., Tchelepi, H.A., Jenny, P., and Dechant, L. 2002a. Implementation of a Flux ContinuousFinite-Difference Method for Stratigraphic, Hexahedron Grids. SPEJ7 (3): 267-277. SPE-80117-PA. DOI: 10.2118/80117-PA.
Lee, S.H., Jenny, P., and Tchelepi, H.A. 2002b. A Finite-Volume Method withHexahedral Multiblock Grids for Modeling Flow in Porous Media. ComputationalGeosciences 6: 353-379.
Lunati, E. and Jenny, P. 2006. Multi-Scale Finite-Volume Method for DensityDriven Flow in Porous Media. Computational Geosciences. Submitted.
Spillette, A.G., Hillestad, J.G., and Stone, H.L. 1973. A High-Stability Sequential SolutionApproach to Reservoir Simulation. Paper SPE 4542 presented at the 48thAnnual Fall Meeting of SPE of AIME, Las Vegas, Nevada, 30 September-3 October.DOI: 10.2118/4542-MS.
Stone, H. and Garder, A.O.J. 1961. Analysis of Gas-Cap or Dissolved-GasDrive Reservoirs. Trans. AIME, 222.
Watts, J.W. 1986. ACompositional Formulation of the Pressure and Saturation Equations.SPERE 1 (3): 243-252. SPE-12244-PA. DOI: 10.2118/12244-PA.
Wolfsteiner, C., Lee, S.H., and Tchelepi, H.A. 2006. Well Modeling in theMultiscale Finite Volume Method for Subsurface Flow Simulation. MultiscaleModel. Simul. 5: 900 - 917.