1887

Abstract

Summary

A large number of multiscale methods have been developed based on the same basic concept: Solve localized flow problems to estimate the local effects of fine-scale petrophysical properties. Use the resulting multiscale basis functions to pose a global flow problem a coarser grid. Reconstruct conservative fine-scale approximations from the coarse-scale solution. By extending the basic concept with iteration cycles and additional local stages, one can systematically drive the fine-scale residual towards machine precision. Posed algebraically, this can be seen as a set of restriction operators for computing a reduced global problem and a set of prolongation operators for constructing conservative fine-scale approximations.

Such multiscale finite-volume methods have been extensively developed for Cartesian grids in the literature. The industry, however, uses very complex with unstructured connections and degenerate cell geometries to represent realistic structural frameworks and stratigraphic architectures. A successful multiscale method should therefore be able to handle unstructured polyhedral grids, both on the fine and coarse scale, and be as flexible as possible to enable automatic coarse partitionings that adapt to wells and geological features in a way that ensures optimal accuracy for a chosen level of coarsening.

Herein, we will discuss a compare a set of prolongation operators that can be combined with finite-element or finite-volume restriction operators to form different multiscale finite-volume methods. We consider the MsFV prolongation operator (developed on a dual coarse grid with unitary at coarse block vertices), the more recent MsTPFA operator (developed on primal grid with unitary flux across coarse block faces), as well as a simplified constant prolongation operator. The methods will be compared on a variety of test cases ranging from simple synthetic grids to highly complex, real-world, field models. Our discussion will focus on flexibility wrt (coarse) grids and tendency of creating oscillatory approximations. In addition, we will look at various methods for improving the methods’ convergence properties when used as preconditioners, as well as for generating novel prolongation operators.

This is relevant for oil recovery because:

  • Multiscale methods may provide a way to significantly speed up reservoir simulation and make previously intractable problems possible to solve.
  • The extension of such methods to industry standard grids used for reservoir modelling enables the evaluation of the methods on real world models
  • The construction of basis functions for multiscale methods may have direct connections to the process of upscaling rock derived properties such as transmissibility

Loading

Article metrics loading...

/content/papers/10.3997/2214-4609.20141775
2014-09-08
2024-03-28
Loading full text...

Full text loading...

References

  1. Aarnes, J.E., Kippe, V. and Lie, K.A.
    [2005] Mixed multiscale finite elements and streamline methods for reservoir simulation of large geomodels. Adv. Water Resour., 28(3), 257–271, doi:10.1016/j.advwatres.2004.10.007.
    https://doi.org/10.1016/j.advwatres.2004.10.007 [Google Scholar]
  2. Aarnes, J.E., Krogstad, S. and Lie, K.A.
    [2006] A hierarchical multiscale method for two-phase flow based upon mixed finite elements and nonuniform coarse grids. Multiscale Model. Simul., 5(2), 337–363, ISSN 1540-3459, doi:10.1137/050634566.
    https://doi.org/10.1137/050634566 [Google Scholar]
  3. [2008] Multiscale mixed/mimetic methods on corner-point grids. Comput. Geosci., 12(3), 297–315, ISSN 1420-0597, doi:10.1007/s10596‑007‑9072‑8.
    https://doi.org/10.1007/s10596-007-9072-8 [Google Scholar]
  4. Alpak, F.O., Pal, M. and Lie, K.A.
    [2012] A multiscale method for modeling flow in stratigraphically complex reservoirs. SPE J., 17(4), 1056–1070, doi:10.2118/140403‑PA.
    https://doi.org/10.2118/140403-PA [Google Scholar]
  5. Chen, Z. and Hou, T.
    [2003] A mixed multiscale finite element method for elliptic problems with oscillating coefficients. Math. Comp., 72, 541–576, doi:10.1090/S0025‑5718‑02‑01441‑2.
    https://doi.org/10.1090/S0025-5718-02-01441-2 [Google Scholar]
  6. Christie, M.A. and Blunt, M.J.
    [2001] Tenth SPE comparative solution project: A comparison of upscaling techniques. SPE Reservoir Eval. Eng., 4, 308–317, doi:10.2118/72469‑PA, url: http://www.spe.org/csp/.
    https://doi.org/10.2118/72469-PA [Google Scholar]
  7. Hajibeygi, H.
    [2011] Iterative multiscale finite volume method for multiphase flow in porous media with complex physics. Ph.D. thesis, ETH Zurich, doi:10.3929/ethz‑a‑006696714.
    https://doi.org/10.3929/ethz-a-006696714 [Google Scholar]
  8. Hajibeygi, H., Bonfigli, G., Hesse, M.A. and Jenny, P.
    [2008] Iterative multiscale finite-volume method. J. Comput. Phys, 227(19), 8604–8621, doi:10.1016/j.jcp.2008.06.013.
    https://doi.org/10.1016/j.jcp.2008.06.013 [Google Scholar]
  9. Hou, T. and Wu, X.H.
    [1997] A multiscale finite element method for elliptic problems in composite materials and porous media. J. Comput. Phys., 134, 169–189, doi:10.1006/jcph.1997.5682.
    https://doi.org/10.1006/jcph.1997.5682 [Google Scholar]
  10. IO Center, NTNU
    [2012] The Norne benchmark case. url: http://www.ipt.ntnu.no/∼norne/wiki/doku.php.
    [Google Scholar]
  11. Jenny, P., Lee, S.H. and Tchelepi, H.A.
    [2003] Multi-scale finite-volume method for elliptic problems in subsurface flow simulation. J. Comput. Phys., 187, 47–67, doi:10.1016/S0021‑9991(03)00075‑5.
    https://doi.org/10.1016/S0021-9991(03)00075-5 [Google Scholar]
  12. [2004] Adaptive multiscale finite-volume method for multiphase flow and transport in porous media. Multiscale Model. Simul., 3(1), 50–64, doi:10.1137/030600795.
    https://doi.org/10.1137/030600795 [Google Scholar]
  13. Karypis, G. and Kumar, V.
    [1998] A fast and high quality multilevel scheme for partitioning irregular graphs. SIAM Journal on scientific Computing, 20(1), 359–392.
    [Google Scholar]
  14. Krogstad, S., Lie, K.A. and Skaflestad, B.
    [2012] Mixed multiscale methods for compressible flow. Proceedings of ECMOR XIII–13th European Conference on the Mathematics of Oil Recovery, EAGE, Biarritz, France.
    [Google Scholar]
  15. Lie, K.A., Krogstad, S., Ligaarden, I.S., Natvig, J.R., Nilsen, H. and Skaflestad, B.
    [2012] Open-source MAT-LAB implementation of consistent discretisations on complex grids. Comput. Geosci., 16, 297–322, doi: 10.1007/s10596‑011‑9244‑4.
    https://doi.org/10.1007/s10596-011-9244-4 [Google Scholar]
  16. Lunati, I. and Jenny, P.
    [2006] Multiscale finite-volume method for compressible multiphase flow in porous media. J. Comput. Phys., 216(2), 616–636, doi:10.1016/j.jcp.2006.01.001.
    https://doi.org/10.1016/j.jcp.2006.01.001 [Google Scholar]
  17. Lunati, I. and Lee, S.H.
    [2009] An operator formulation of the multiscale finite-volume method with correction function. Multiscale Model. Simul., 8(1), 96–109, doi:10.1137/080742117.
    https://doi.org/10.1137/080742117 [Google Scholar]
  18. Lunati, I., Tyagi, M. and Lee, S.H.
    [2011] An iterative multiscale finite volume algorithm converging to the exact solution. J. Comput. Phys., 230(5), 1849–1864, doi:10.1016/j.jcp.2010.11.036.
    https://doi.org/10.1016/j.jcp.2010.11.036 [Google Scholar]
  19. Metis website
    [2012] Metis – serial graph partitioning and fill-reducing matrix ordering. Url: http://glaros.dtc.umn.edu/gkhome/views/metis.
    [Google Scholar]
  20. Møyner, O.
    [2012] Multiscale finite-volume methods on unstructured grids. Master’s thesis, Norwegian University of Science and Technology, Trondheim.
    [Google Scholar]
  21. Møyner, O. and Lie, K.A.
    [2014a] The multiscale finite-volume method on stratigraphic grids. SPE J., doi:10.2118/163649‑PA.
    https://doi.org/10.2118/163649-PA [Google Scholar]
  22. [2014b] A multiscale two-point flux-approximation method. J. Comput. Phys.
    [Google Scholar]
  23. MRST
    [2014] The MATLAB Reservoir Simulation Toolbox, version 2014a. http://www.sintef.no/MRST/.
    [Google Scholar]
  24. Nordbotten, J.M. and Bjørstad, P.
    [2008] On the relationship between the multiscale finite-volume method and domain decomposition preconditioners. Comput. Geosci., 12(3), 367–376.
    [Google Scholar]
  25. Nordbotten, J.M., Keilegavlen, E. and Sandvin, A.
    [2012] Mass conservative domain decomposition for porous media flow. In: Petrova, R. (Ed.)Finite Volume Method-Powerful Means of Engineering Design. InTech Europe, Rijeka, Croatia, 235–256, doi:10.5772/38700.
    https://doi.org/10.5772/38700 [Google Scholar]
  26. Sandvin, A., Keilegavlen, E. and Nordbotten, J.M.
    [2013] Auxiliary variables for 3d multiscale simulations in heterogeneous porous media. J. Comput. Phys, 238, 141–153, doi:10.1016/j.jcp.2012.12.016.
    https://doi.org/10.1016/j.jcp.2012.12.016 [Google Scholar]
  27. Sandvin, A., Nordbotten, J. and Aavatsmark, I.
    [2011] Multiscale mass conservative domain decomposition preconditioners for elliptic problems on irregular grids. Comput. Geosci., 15, 587–602, doi:10.1007/s10596‑011‑9226‑6.
    https://doi.org/10.1007/s10596-011-9226-6 [Google Scholar]
  28. Wang, Y., Hajibeygi, H. and Tchelepi, H.A.
    [2012] Algebraic multiscale linear solver for heterogeneous elliptic problems. ECMOR XIII - 13th European Conference on the Mathematics of Oil Recovery, EAGE, Biarritz, France.
    [Google Scholar]
  29. [2014] Algebraic multiscale solver for flow in heterogeneous porous media. J. Comput. Phys., 259, 284–303, doi:10.1016/j.jcp.2013.11.024.
    https://doi.org/10.1016/j.jcp.2013.11.024 [Google Scholar]
  30. Zhou, H. and Tchelepi, H.A.
    [2012] Two-stage algebraic multiscale linear solver for highly heterogeneous reservoir models. SPE J., 17(2), 523–539, doi:10.2118/141473‑PA.
    https://doi.org/10.2118/141473-PA [Google Scholar]
http://instance.metastore.ingenta.com/content/papers/10.3997/2214-4609.20141775
Loading
/content/papers/10.3997/2214-4609.20141775
Loading

Data & Media loading...

This is a required field
Please enter a valid email address
Approval was a Success
Invalid data
An Error Occurred
Approval was partially successful, following selected items could not be processed due to error