1887

Abstract

Summary

Development of boundary aligned grid generation is presented together with comparative performance of cell-vertex versus cell-centred CVD-MPFA finite-volume formulations using equivalent degrees of freedom. When generating structured or unstructured grids for reservoir simulation, classical key constraints involve boundary aligned grid generation with control-volume boundaries aligned with solid walls and geological features such as layers, shale barriers, fractures, faults, pinchouts and multilateral wells.

The schemes used are control-volume distributed (CVD) with flow variables and rock properties sharing the same control-volume location and are comprised of a multipoint flux family formulation (CVD-MPFA). Consequently a natural choice is for primal grid cells to act as control-volumes, then grid generation can be performed with primal grid cell boundaries being aligned with key interior constraint boundaries. This naturally leads to cell-centred approximation, where flow variables and rock properties are associated with grid cell centres.

The alternative is to employ cell-vertex approximation which uses far fewer approximation points on a given unstructured grid. In this case control-volumes are constructed around primal grid vertices. The grid generation process is less straight forward since control-volumes must be constrained to satisfy interior boundary alignment. A novel grid generation procedure is proposed that automates control-volume boundary alignment and yields a Voronoi mesh. The actual grid is then generated such that dual boundaries are aligned with key internal constraint boundaries and cell-vertex approximation becomes the natural choice. In this case flow variables and rock properties are associated with grid cell vertices and their dual control-volumes.

The relative benefits of both types of approximation is made clear in terms of flow resolution and degrees of freedom required.

Loading

Article metrics loading...

/content/papers/10.3997/2214-4609.20141858
2014-09-08
2020-04-02
Loading full text...

Full text loading...

References

  1. Aavatsmark, I.
    [2002] An Introduction to Multipoint Flux Approximations for Quadrilateral Grids. Computational Geosciences6(3),405–432.
    [Google Scholar]
  2. Aavatsmark, I., BarkveT., Boe, O. and MannsethM.
    [1998] Discretization on Unstructured Grids for Inhomogeneous, Anisotropic Media. Part I: Derivation of the Methods. SIAM Journal on Scientific Computing19: 1716–.
    [Google Scholar]
  3. ArnoldN.D.
    [1990] Mixed finite element methods for elliptic problems. Comput. Meth. Appl. Mech. Engrg, 281–300.
    [Google Scholar]
  4. Barth, T.
    [1994] Aspects of Unstructured Grids and Finite Volume Solvers for the Euler and Navier-Stokes Equations. Computational Fluid Dynamics No. 1994–04 in Lecture Series, von Karman Institute for Fluid Dynamics
    [Google Scholar]
  5. Branets, L., Ghai, S. S., Lyons, S. L., and Wu, X. H.
    [2009]. Efficient and Accurate Reservoir Modeling Using Adaptive Gridding with Global Scale Up. Society of Petroleum Engineers, SPE Reservoir Simulation Symposium, 2–4 February, The Woodlands, Texas(SPE-118946-MS).
    [Google Scholar]
  6. Branets, L. V., Kartasheva, E., Krasnogorov, I., Kubyak, V. and Wu, X.H.
    [2012] Method and apparatus for reservoir modeling and simulation. EP Patent App. EP20,100,830,349.
    [Google Scholar]
  7. Castellini,A.
    Flow Based Grids for Reservoir Simulation MS thesis, Stanford University, Stanford, CA USA.
    [Google Scholar]
  8. Courrioux, G., Nullans, S., Guillen, A., Boissonnat, J.D., Repusseau, P., Renaud, X. and Thibaut, M.
    [2001] 3D volumetric modelling of Cadomian terranes (Northern Brittany,France): an automatic method using Voronoi diagrams. Tectonophysics, 331, 196–.
    [Google Scholar]
  9. Du, Q., Faber, V. and Gunzburger, M.
    [1999] Centroidal Voronoi Tessellations: Applications and Algorithms. Society for Industrial and Applied Mathematics (SIAM) Rev., 41(4), 637–676.
    [Google Scholar]
  10. Diskin, B., Thomas, J., Nielsen, E., White, J. and Nishikawa, H.
    [2009] Comparison of Node-Centered and Cell-Centered Unstructured Finite-Volume Discretizations Part I: Viscous Fluxes. 47th AIAA Aerospace Sciences Meeting including The New Horizons Forum and Aerospace Exposition. 5–8 JanuaryOrlando Florida.
    [Google Scholar]
  11. Edwards, M.G. and Rogers, C.F.
    [1998] Finite volume discretization with imposed flux continuity for the general tensor pressure equation. Computational Geosciences, 2, 290–.
    [Google Scholar]
  12. Edwards, M.G.
    [2002] Unstructured, Control-Volume Distributed, Full-Tensor Finite Volume Schemes With Flow Based Grids. Computational Geosciences, 6, 433 – 452.
    [Google Scholar]
  13. Edwards, M.G. and Zheng, H.
    [2008] A quasi-positive family of continuous Darcy-flux finite-volume schemes with full pressure support. Journal of Computational Physics, 227, 9364–.
    [Google Scholar]
  14. [2010] Double-families of quasi-positive Darcy-flux approximations with highly anisotropic tensors on structured and unstructured grids. Journal of Computational Physics, 229, 625–.
    [Google Scholar]
  15. EdwardsM G, ZhengH., LamineS. and PalM.
    [2011] Continuous Elliptic and Multi-dimensional Hyperbolic Darcy-Flux Finite-volume Methods. Computers and Fluids, 46 (1), 12–22
    [Google Scholar]
  16. Fortune, S.
    [1997]. In: Goodman, J.E., and O’Rourke, J.(Eds.) Voronoi diagrams and Delaunay triangulations – Handbook of discrete and computational geometry, 377–388.
    [Google Scholar]
  17. Friis, H.A., Edwards, M.G. and Mykkeltveit, J.
    [2008] Symmetric Positive Definite Flux-Continuous Full-Tensor Finite-Volume Schemes on Unstructured Cell-Centered Triangular Grids. SIAM Journal on Scientific Computing, 31, 1220–.
    [Google Scholar]
  18. Fung, L.S.K., Ding, X.Y., and Dogru, A.H.
    [2013] An Unstructured Gridding Method for Densely-Spaced Complex Wells in Full-Field Reservoir Simulation. SPE Reservoir Simulation Symposium (The Woodlands, Texas, USA), 995–1007.
    [Google Scholar]
  19. Friis, H.A. and Edwards, M.G.
    [2011] Family of MPFA finite-volume schemes with full pressure support for the general tensor pressure equation on cell-centered triangular grids. em Journal of Computational Physics, 230(1), 205–231.
    [Google Scholar]
  20. George, P.L., Borouchaki, H.
    [1998] In: HERMES(Ed) Delaunay Triangulation – Delaunay Triangulation and Meshing Applications to Finite Elements, 37,238.
    [Google Scholar]
  21. Guibas, L., and Stolfi, J.
    [1985] Primitives for the manipulation of general subdivisions and the computation of Voronoi. ACM Trans. Graph., 4(2), 74–123.
    [Google Scholar]
  22. Green, P.J., and Sibson, R.
    [1978] Computing Dirichlet Tessellations in the Plane. The Computer Journal, 21(2), 168–173.
    [Google Scholar]
  23. Hægland, H., Assteerawatt, A., Dahle, H.K., Eigestad, G.T. and Helmig, R.
    [2009] Comparison of cell- and vertex-centered discretization methods for flow in a two-dimensional discrete-fracture-matrix system. Advances in Water Resources, 32(12), 1740–1755.
    [Google Scholar]
  24. Jenkins, C., Ouenes, A., Zellou, A. and Wingard, J.
    [2009] Quantifying and predicting naturally fractured reservoir behavior with continuous fracture models. The American Association of Petroleum Geologists, 93(11), 1597–1608.
    [Google Scholar]
  25. Karimi-Fard, M.
    [2004] Growing Region Technique Applied to Grid Generation of Complex Fractured Porous Media. European Association of Geoscientists & Engineers – 9th European Conference on the Mathematics of Oil Recovery
    [Google Scholar]
  26. KeilegavlenE. and AavatsmarkI.
    [2011] Monotonicity for MPFA methods on triangular grids.Computational Geosciences,15(1), 3–16.
    [Google Scholar]
  27. Liebling, T.M., and Pournin, L.
    [2012] Voronoi Diagrams and Delaunay Triangulations: Ubiquitous Siamese Twins. In: Martin GrotschelM.(Ed.)Journal of Documenta Mathematica – Optimization Stories, 419–431.
    [Google Scholar]
  28. Lee, D.T., and Schachter, B.J.
    [1980] Two algorithms for constructing a Delaunay triangulation.International Journal of Computer & Information Sciences, 9(3), 219–242.
    [Google Scholar]
  29. Lee, S.H., L. J.Durlofsky, M. F.Lough and W. H.Chen
    . [1998] Finite difference simulation of geologically complex reservoirs with tensor permeabilities.SPERE, 1, 567–574.
    [Google Scholar]
  30. Lawson, C.L.
    [1986] Properties of n-dimensional triangulations.Computer Aided Geometric Design, 3(4), 231–246.
    [Google Scholar]
  31. Lipnikov, K., Svyatskiy, D. and Vassilevski, Y.
    [2009] Interpolation-free monotone finite volume method for diffusion equations on polygonal meshes. Journal of Computational Physics, 228(3), 703–716.
    [Google Scholar]
  32. Löner, R. and Parikh, P.
    [1988] Three dimensional grid generation by advancing front method. International Journal for Numerical Methods in Fluids, 8, 1135–1149.
    [Google Scholar]
  33. Mahani, H. and Evazi, M.
    [2010] Vorticity-Based Perpendicular-Bisector Grids for Improved Upscaling of Two-Phase Flow. SPEJ, Dec 2010, 995–1008.
    [Google Scholar]
  34. Marina, G. and Jon, R.
    [1996] An efficient algorithm for construction of the power diagram from the voronoi diagram in the plane. International Journal of Computer Mathematics, 61(1–2), 49–61.
    [Google Scholar]
  35. Marcum, D.L., and Weatherill, N.P.
    [1995] Unstructured grid generation using iterative point insertion and local reconnection. AIAA Journal, 33, 1619–1625.
    [Google Scholar]
  36. Mavriplis, D.J.
    [1995] An Advancing Front Delaunay Triangulation Algorithm Designed for Robustness. Journal of Computational Physics, 117, 90–101.
    [Google Scholar]
  37. Merland, R., Lévy, B. and CaumonG.
    [2011] Building PEBI Grids Conforming to 3D Geological Features Using Centroidal Voronoi Tessellations. In: Marschallinger, R. and Zolb, R. (Eds.) Mathematical Geosciences at the Crossroads of Theory and Practice – Proceedings ofIAMG conference. Slazburg, 254–265.
    [Google Scholar]
  38. Mlacnik, M., Durlofsky, L.J., and Heinemann, Z.E.
    [2007] Sequentially Adapted Flow-Based PEBI Grids for Reservoir Simulation. Society of Petroleum Engineers, 11(3),317–327.
    [Google Scholar]
  39. Müller, J. D., Roe, P. L., and Deconinck, H.
    [1993]A frontal approach for internal node generation in Delaunay triangulations. International Journal for Numerical Methods in Fluids, 17(3), 241–255.
    [Google Scholar]
  40. Hassan, O., Probert, E.J., Morgan, K. and Weatherill, N.P.
    [2000] Unsteady flow simulation using unstructured meshes. Computer Methods in Applied Mechanics and Engineering,189(4),1247–1275.
    [Google Scholar]
  41. Palagi, C.L. and Aziz, K.
    [2007] Use of Voronoi Grid in Reservoir Simulation.Society of Petroleum Engineers, 2(2) 69–77.
    [Google Scholar]
  42. Pal, M. and Edwards, M.G.
    [2010] q-Families of CVD(MPFA) Schemes on General Elements: Numerical Convergence and the Maximum Principle. Archives of Computational Methods in Engineering, 17(2), 137–189.
    [Google Scholar]
  43. Rebay, S.
    [1993] Efficient Unstructured Mesh Generation by Means of Delaunay Triangulation and Bowyer-Watson Algorithm. Journal of Computational Physics, 106(1) 125–138.
    [Google Scholar]
  44. RussellT. and WheelerM.F.
    [1984] Finite element and finite difference methods for continuous flows in porous media. The Mathematics of Reservoir Simulation, R. E.Ewing, ed., Frontiers in Applied Mathematics, SIAM, 35–106
    [Google Scholar]
  45. Sahimi, M., Darvishi, R., Haghighi, M. and Rasaei, M.R.
    [2010] Upscaled Unstructured Computational Grids for Efficient Simulation of Flow in Fractured Porous Media. Transport in Porous Media, 83, 218–.
    [Google Scholar]
  46. Sazonov, I., Wang, D., Hassan, O., Morgan, K. and Weatherill, N.
    [2006] A stitching method for the generation of unstructured meshes for use with co-volume solution techniques. Computer Methods in Applied Mechanics and Engineering, 195, 1845–.
    [Google Scholar]
  47. Sibson, R.
    [1978] Locally equiangular triangulations, The Computer Journal, 21(3), 243–245.
    [Google Scholar]
  48. Sloan, S.W.
    [1993] A fast algorithm for generating constrained delaunay triangulations. Computers & Structures, 47(3), 441–450.
    [Google Scholar]
  49. [1987] A fast algorithm for constructing Delaunay triangulations in the plane. Advances in Engineering Software,9(1), 34–55.
    [Google Scholar]
  50. Su, P. and Scot, R., L.
    [1997] A comparison of sequential Delaunay triangulation algorithms. Computational Geometry–11th ACM Symposium on Computational Geometry, 7(5–6), 361–385.
    [Google Scholar]
  51. Sloan, S.W., and Houlsby, G.T.
    [1984] An implementation of Watson’s algorithm for computing 2-dimensional delaunay triangulations. Advances in Engineering Software, 6(4),192–197.
    [Google Scholar]
  52. Thompson, J.F. and Weatherill, N.P.
    [1998]. In: Thompson, J.F., Soni, B, K. and Weatherill, N.P.(Eds.)Fundamental Concepts and Approaches – Handbook of Grid Generation, 419–431.
    [Google Scholar]
  53. Verma, S. K.
    [1996] Flexible grids for reservoir simulation. PhD thesis, Stanford University.
    [Google Scholar]
  54. Watson, D. F.
    [1981] Computing the n-dimensional Delaunay tessellation with application to Voronoi polytopes. The Computer Journal, 24(2), 167–172.
    [Google Scholar]
  55. WenX. H., L.J.Durlofsky and M.G.Edwards
    [2003] Upscaling of Channel Systems in Two Dimensions Using Flow Based Grids. Transport in Porous Media., 51343–366.
    [Google Scholar]
  56. WheelerM. F. and YotovI.
    [2006] A multipoint flux mixed finite element method. SIAM J. Numer. Anal., 44:52082–2106..
    [Google Scholar]
http://instance.metastore.ingenta.com/content/papers/10.3997/2214-4609.20141858
Loading
/content/papers/10.3997/2214-4609.20141858
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