Operator-Based Multiscale Method for Compressible Flow
- Hui Zhou (Stanford University) | Hamdi A. Tchelepi (Stanford University)
- Document ID
- Society of Petroleum Engineers
- SPE Journal
- Publication Date
- June 2008
- Document Type
- Journal Paper
- 267 - 273
- 2008. Society of Petroleum Engineers
- 4.3.4 Scale, 5.5.1 Simulator Development, 5.5 Reservoir Simulation, 4.1.2 Separation and Treating, 5.3.2 Multiphase Flow, 5.5.8 History Matching, 5.2.1 Phase Behavior and PVT Measurements
- 2 in the last 30 days
- 395 since 2007
- Show more detail
- View rights & permissions
|SPE Member Price:||USD 12.00|
|SPE Non-Member Price:||USD 35.00|
Multiscale methods have been developed for accurate and efficient numerical solution of flow problems in large-scale heterogeneous reservoirs. A scalable and extendible Operator-Based Multiscale Method (OBMM) is described here. OBMM is cast as a general algebraic framework. It is natural and convenient to incorporate more physics in OBMM for multiscale computation. In OBMM, two operators are constructed: prolongation and restriction. The prolongation operator is constructed by assembling the multiscale basis functions. The specific form of the restriction operator depends on the coarse-scale discretization formulation (e.g., finitevolume or finite-element). The coarse-scale pressure equation is obtained algebraically by applying the prolongation and restriction operators to the fine-scale flow equations. Solving the coarse-scale equation results in a high-quality coarse-scale pressure. The finescale pressure can be reconstructed by applying the prolongation operator to the coarse-scale pressure. A conservative fine-scale velocity field is then reconstructed to solve the transport (saturation) equation. We describe the OBMM approach for multiscale modeling of compressible multiphase flow. We show that extension from incompressible to compressible flows is straightforward. No special treatment for compressibility is required. The efficiency of multiscale formulations over standard fine-scale methods is retained by OBMM. The accuracy of OBMM is demonstrated using several numerical examples including a challenging depletion problem in a strongly heterogeneous permeability field (SPE 10).
The accuracy of simulating subsurface flow relies strongly on the detailed geologic description of the porous formation. Formation properties such as porosity and permeability typically vary over many scales. As a result, it is not unusual for a detailed geologic description to require 107-108 grid cells. However, this level of resolution is far beyond the computational capability of state-of-the-art reservoir simulators (106 grid cells). Moreover, in many applications, large numbers of reservoir simulations are performed (e.g., history matching, sensitivity analysis and stochastic simulation). Thus, it is necessary to have an efficient and accurate computational method to study these highly detailed models.
Multiscale formulations are very promising due to their ability to resolve fine-scale information accurately without direct solution of the global fine-scale equations. Recently, there has been increasing interest in multiscale methods. Hou and Wu (1997) proposed a multiscale finite-element method (MsFEM) that captures the fine-scale information by constructing special basis functions within each element. However, the reconstructed fine-scale velocity is not conservative. Later, Chen and Hou (2003) proposed a conservative mixed finite-element multiscale method. Another multiscale mixed finite element method was presented by Arbogast (2002) and Arbogast and Bryant (2002). Numerical Green functions were used to resolve the fine-scale information, which are then coupled with coarse-scale operators to obtain the global solution. Aarnes (2004) proposed a modified mixed finite-element method, which constructs special basis functions sensitive to the nature of the elliptic problem. Chen et al. (2003) developed a local-global upscaling method by extracting local boundary conditions from a global solution, and then constructing coarse-scale system from local solutions. All these methods considered incompressible flow in heterogeneous porous media where the pressure equation is elliptic.
A multiscale finite-volume method (MsFVM) was proposed by Jenny et al. (2003, 2004, 2006) for heterogeneous elliptic problems. They employed two sets of basis functions--dual and primal. The dual basis functions are identical to those of Hou and Wu (1997), while the primal basis functions are obtained by solving local elliptic problems with Neumann boundary conditions calculated from the dual basis functions.
Existing multiscale methods (Aarnes 2004; Arbogast 2002; Chen and Hou 2003; Hou and Wu 1997; Jenny et al. 2003) deal with the incompressible flow problem only. However, compressibility will be significant if a gas phase is present. Gas has a large compressibility, which is a strong function of pressure. Therefore, there can be significant spatial compressibility variations in the reservoir, and this is a challenge for multiscale modeling. Very recently, Lunati and Jenny (2006) considered compressible multiphase flow in the framework of MsFVM. They proposed three models to account for the effects of compressibility. Using those models, compressibility effects were represented in the coarse-scale equations and the reconstructed fine-scale fluxes according to the magnitude of compressibility.
Motivated to construct a flexible algebraic multiscale framework that can deal with compressible multiphase flow in highly detailed heterogeneous models, we developed an operator-based multiscale method (OBMM). The OBMM algorithm is composed of four steps: (1) constructing the prolongation and restriction operators, (2) assembling and solving the coarse-scale pressure equations, (3) reconstructing the fine-scale pressure and velocity fields, and (4) solving the fine-scale transport equations.
OBMM is a general algebraic multiscale framework for compressible multiphase flow. This algebraic framework can also be extended naturally from structured to unstructured grid. Moreover, the OBMM approach may be used to employ multiscale solution strategies in existing simulators with a relatively small investment.
|File Size||282 KB||Number of Pages||7|
Aarnes, J. E. 2004. On the useof a mixed multiscale finite element method for greater flexibility andincreased speed or improved accuracy in reservoir simulation. MultiscaleModeling and Simulation 2 (3): 421-439. doi: 10.1137/030600655.
Arbogast, T. 2002. Implementation of a locallyconservative numerical subgrid upscaling scheme for two-phase Darcy flow.Computational Geosciences 6 (3-4): 453-481. doi:10.1023/A:1021295215383.
Arbogast, T. and Bryant, S. L. 2002. A Two-Scale Numerical SubgridTechnique for Waterflood Simulations. SPEJ 7 (4): 446-457.SPE-81909-PA doi: 10.2118/81909-PA
Chen, Y., Durlofsky, L. J., Gerritsen, M., and Wen, X. H. 2003. A coupled local-globalupscaling approach for simulating flow in highly heterogeneous formations.Advances in Water Resources 26 (10): 1041-1060. doi:10.1016/S0309-1708(03)00101-5.
Chen, Z. and Hou, T. 2003. A mixed finite elementmethod for elliptic problems with rapidly oscillating coefficients.Mathematic Computation 72 (242): 541-576. doi:10.1090/S0025-5718-02-01441-2.
Christie, M. A. and Blunt, M. J. 2001. Tenth SPE Comparative SolutionProject: A Comparison of Upscaling Techniques. SPEREE 4 (4):308-317. SPE-72469-PA doi: 10.2118/72469-PA
Hou, T. and Wu, X. H. 1997. A multiscale finite elementmethod for elliptic problems in composite materials and porous media.Journal of Computational Physics 134 (1): 169-189. doi:10.1006/jcph.1997.5682.
Jenny, P., Lee, S. H., and Tchelepi, H. A. 2003. Multiscale finite volumemethod for elliptic problems in subsurface flow simulation. Journal ofComputational Physics 187 (1): 47-67. doi:10.1016/S0021-9991(03)00075-5.
Jenny, P., Lee, S. H., and Tchelepi, H. A. 2004. Adaptive multiscale finite-volumemethod for multiphase flow and transport in porous media. MultiscaleModeling and Simulation 3 (1): 50-64. doi: 10.1137/030600795.
Jenny, P., Lee, S. H., and Tchelepi, H. A. 2006. Adaptive fully implicitmultiscale finite-volume method for multi-phase flow and transport inheterogeneous porous media. Journal of Computational Physics217 (2): 627-641. doi: 10.1016/j.jcp.2006.01.028.
Lunati, I. and Jenny, P. 2006. Multiscale finite-volumemethod for compressible multiphase flow in porous media. Journal ofComputational Physics 216 (2): 616-636. doi:10.1016/j.jcp.2006.01.001.
Tchelepi, H. A., Jenny, P., Lee, S. H., and Wolfsteiner, C. 2007. Adaptive Multiscale Finite-VolumeFramework for Reservoir Simulation. SPEJ 12 (2): 188-195.SPE-93395-PA doi: 10.2118/93395-PA
Zhou, H. 2006. Operator based multiscale method for compressible Flow. MSthesis, Stanford, California: Stanford University.