Abstract
Multiscale simulation is a promising approach to facilitate direct simulation of large and complex grid models for highly heterogeneous petroleum reservoirs. Unlike traditional simulation, approaches based on upscaling/downscaling, multiscale methods seek to solve the full flow problem by incorporating subscale heterogeneities into local discrete approximation spaces. We consider a multiscale formulation based on a hierarchical grid approach, where basis functions with subgrid resolution are computed numerically to correctly and accurately account for subscale variations from an underlying (fine-scale) geomodel when solving the global flow equations on a coarse grid. By using multiscale basis functions to discretise the global flow equations on a (moderately sized) coarse grid, one can retain the efficiency of an upscaling method and, at the same time, produce detailed and conservative velocity fields on the underlying fine grid. For pressure equations, the multiscale mixed finite-element method (MsMFEM) has been shown to be a particularly versatile approach. In this paper, we extend the method to corner-point grids, which is the industry standard for modelling complex reservoir geology. To implement MsMFEM, one needs a discretisation method for solving local flow problems on the underlying fine grids. In principle, any stable and conservative method can be used. Here, we use a mimetic discretisation, which is a generalisation of mixed finite elements that gives a discrete inner product, allows for polyhedral elements, and can (easily) be extended to curved grid faces. The coarse grid can, in principle, be any partition of the subgrid, where each coarse block is a connected collection of subgrid cells. However, we argue that, when generating coarse grids, one should follow certain simple guidelines to achieve improved accuracy. We discuss partitioning in both index space and physical space and suggest simple processing techniques. The versatility and accuracy of the new multiscale mixed methodology is demonstrated on two corner-point models: a small Y-shaped sector model and a complex model of a layered sedimentary bed. A variety of coarse grids, both violating and obeying the above mentioned guidelines, are employed. The MsMFEM solutions are compared with a reference solution obtained by direct simulation on the subgrid.
Similar content being viewed by others
References
Aarnes, J.E.: On the use of a mixed multiscale finite element method for greater flexibility and increased speed or improved accuracy in reservoir simulation. Multiscale Model. Simul. 2(3), 421–439 (2004)
Aarnes, J.E., Kippe, V., Lie, K.A.: Mixed multiscale finite elements and streamline methods for reservoir simulation of large geomodels. Adv. Water Resour. 28(3), 257–271 (2005)
Aarnes, J.E., Krogstad, S., Lie, K.A.: A hierarchical multiscale method for two-phase flow based upon mixed finite elements and nonuniform coarse grids. Multiscale Model. Simul. 5(2), 337–363 (2006)
Aarnes, J.E., Krogstad, S., Lie, K.A., Laptev, V.: Multiscale mixed methods on corner-point grids: mimetic versus mixed subgrid solvers. Technical report, SINTEF (2006)
Aavatsmark, I., Barkve, T., Bøe, Ø., Mannseth, T.: Discretization on unstructured grids for inhomogeneous, anisotropic media. Part I: derivation of the methods. SIAM J. Sci. Comput. 19(5), 1700–1716 (1998)
Aavatsmark, I., Barkve, T., Bøe, Ø., Mannseth, T.: Discretization on unstructured grids for inhomogeneous, anisotropic media. Part II: discussion and numerical results. SIAM J. Sci. Comput. 19(5), 1717–1736 (1998)
Arbogast, T.: Numerical subgrid upscaling of two-phase flow in porous media. In: Chen, Z., Ewing, R., Shi, Z.C. (eds.) Lecture Notes in Physics, pp. 35–49. Springer, Berlin (2000)
Arbogast, T.: Analysis of a two-scale, locally conservative subgrid upscaling for elliptic problems. SIAM J. Numer. Anal. 42(2), 576–598 (2004)
Arbogast, T., Boyd, K.: Subgrid upscaling and mixed multiscale finite elements. SIAM J. Numer. Anal. 44(3), 1150–1171 (2006)
Arbogast, T., Minkoff, S.E., Keenan, P.T.: An operator-based approach to upscaling the pressure equation. In: Burganos, V.N. et al. (eds.) Computational Methods in Water Resources XII, Vol. 1: Computational Methods in Contamination and Remediation of Water Resources, pp. 405–412. Computational Mechanics Publications, Southampton (1998)
Arnold, D.N., Brezzi, F.: Mixed and nonconforming finite element methods: implementation, postprocessing and error estimates. RAIRO Model. Math. Anal. Numer. 19(1), 7–32 (1985)
Brezzi, F., Lipnikov, K., Shashkov, M.: Convergence of mimetic finite difference methods for diffusion problems on polyhedral meshes. SIAM J. Numer. Anal. 43, 1872–1895 (2005)
Brezzi, F., Lipnikov, K., Shashkov, M.: Convergence of mimetic finite difference methods for diffusion problems on polyhedral meshes with curved faces. Math. Models Methods Appl. Sci. 16(2), 275–297 (2006)
Brezzi, F., Lipnikov, K., Shashkov, M., Simoncini, V.: A new discretization methodology for diffusion problems on generalized polyhedral meshes. Comput. Methods Appl. Mech. Eng. 196(37–40), 3682–3692 (2007)
Brezzi, F., Lipnikov, K., Simoncini, V.: A family of mimetic finite difference methods on polygonial and polyhedral meshes. Math. Models Methods Appl. Sci. 15, 1533–1553 (2005)
Chen, Z., Hou, T.: A mixed multiscale finite element method for elliptic problems with oscillating coefficients. Math. Comput. 72, 541–576 (2003)
Edwards, M., Lazarov, R., Yotov, I. (eds.): Special issue on locally conservative numerical methods for flows in porous media. Comput. Geosci. 6(3–4), 225–579 (2002)
Gautier, Y., Blunt, M.J., Christie, M.A.: Nested gridding and streamline-based simulation for fast reservoir performance prediction. Comput. Geosci. 3, 295–320 (1999)
Hou, T., Wu, X.H.: A multiscale finite element method for elliptic problems in composite materials and porous media. J. Comput. Phys. 134, 169–189 (1997)
Hughes, T., Feijoo, G., Mazzei, L., Quincy, J.B.: The variational multiscale method—a paradigm for computational mechanics. Comput. Methods. Appl. Mech. Eng. 166, 3–24 (1998)
Jenny, P., Lee, S.H., Tchelepi, H.A.: Multi-scale finite-volume method for elliptic problems in subsurface flow simulation. J. Comput. Phys. 187, 47–67 (2003)
Jenny, P., Lee, S.H., Tchelepi, H.A.: Multi-scale finite-volume method for elliptic problems in subsurface flow simulation. Multiscale Model. Simul. 3(1), 50–64 (2004)
Juanes, R.: A variational multiscale method for the simulation of porous media flow in highly heterogeneous formations. In: Binning, P.J., Engesgaard, P.K., Dahle, H.K., Pinder, G.F., Gray, W.G. (eds.) Proceedings of the XVI International Conference on Computational Methods in Water Resources (2006). http://proceedings.cmwr-xvi.org
Kippe, V., Aarnes, J., Lie, K.A.: A comparison of multiscale methods for elliptic problems in porous media flow. Comput. Geosci. (2008). doi:10.1007/s10596-007-9074-6
Målqvist, A.: Adaptive variational multiscale methods. Ph.D. Thesis, Chalmers University of Technology (2005)
Ponting, D.K.: Corner point geometry in reservoir simulation. In: King, P. (ed.) Proceedings of the 1st European Conference on Mathematics of Oil Recovery, Cambridge, 1989, pp. 45–65. Clarendon, Oxford (1989)
Raviart, P., Thomas, J.: A mixed finite element method for second order elliptic equations. In: Galligani, I., Magenes, E. (eds.) Mathematical Aspects of Finite Element Methods, pp. 292–315. Springer, Berlin Heidelberg New York (1977)
Author information
Authors and Affiliations
Corresponding author
Rights and permissions
About this article
Cite this article
Aarnes, J.E., Krogstad, S. & Lie, KA. Multiscale mixed/mimetic methods on corner-point grids. Comput Geosci 12, 297–315 (2008). https://doi.org/10.1007/s10596-007-9072-8
Received:
Accepted:
Published:
Issue Date:
DOI: https://doi.org/10.1007/s10596-007-9072-8