1887

Abstract

Summary

Grids based on Voronoi diagrams, comprise the dual of Delaunay triangulations (DTs), and remain predominant in reservoir simulation. Voronoi-based grids are locally orthogonal, i.e., PErpendicular BIsectional(PEBI), and permit consistent two-point flux approximation if the permeability field is isotropic, or if the grid is generated to be K-orthogonal. In addition, Voronoi grids can be made to honor classical key constraints, minimizing discretization error, thereby aligning associated control volumes with solid walls, well-trajectories, and geological features: layers, shale barriers, fractures, faults, and pinch-outs.

Typically, in reservoirs formed by deposition, the directional trend in the horizontal plane is not very distinct, whereas across the layers, rock properties may jump by orders of magnitude. The PEBI property associated with DTs is of major significance, and can only be exploited provided the circumcenter is used as the approximation point. This requires that the grids generated be boundary aligned, and comprised of entirely acute simplexes. The control-volume centroid is commonly used as the approximation point, because a geological feature honored acute triangulation cannot be guaranteed in the general case; especially in the presence of complex geometries and geological constraints. Development of a geological feature-based acute DT technique is presented. A boundary-aligned grid generation method is augmented with a mesh reconstruction technique, which can ensure circumcenter containment of a DT. To honor the geological feature idea of protection-circle is used. In the mesh reconstruction technique, each mesh point is optimized iteratively, using a local-advancing front method incorporating the length of opposite edges of the set of simplexes sharing it. The methods presented generate boundary-aligned acute DT, where previously proposed methods fail to ensure the acute DT property. Details of the method will be presented, together with results for a number of test cases that verify consistency of the two-point flux on the resulting boundary-aligned acute Voronoi grids.

Loading

Article metrics loading...

/content/papers/10.3997/2214-4609.201802273
2018-09-03
2020-04-01
Loading full text...

Full text loading...

References

  1. Arbogast, T., Douglas, J.Jr. and Hornung, U.
    [1990] Derivation of the Double Porosity Model of Single Phase Flow via Homogenization Theory.SIAM Journal on Mathematical Analysis, 21(4), 823–836.
    [Google Scholar]
  2. Borouchaki, H. and George, P.
    [1997] Aspects of 2-D Delaunay mesh generation.Int. J. Numer. Meth. Engng, 40, 1957–1975.
    [Google Scholar]
  3. Branets, L., Ghai, S., Lyons, S. and Wu, X.
    [2009] Efficient and Accurate Reservoir Modeling Using Adaptive Gridding with Global Scale Up.SPE Reservoir Simulation Symposium, Society of Petroleum Engineers, Texas.
    [Google Scholar]
  4. 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]
  5. Edwards, M.
    [2006] Higher-resolution hyperbolic-coupled-elliptic flux-continuous CVD schemes on structured and unstructured grids in 2-D.International Journal for Numerical Methods in Fluids, 51(9–10), 1059–1077.
    [Google Scholar]
  6. Edwards, M. and Rogers, C.
    [1998] Finite volume discretization with imposed flux continuity for the general tensor pressure equation.Computational Geosciences, 2, 259–290.
    [Google Scholar]
  7. Eigestad, G. and Klausen, R.
    [2005] On the convergence of the multi-point flux approximation O-method: Numerical experiments for discontinuous permeability.Numerical Methods for Partial Differential Equations, 21(6), 1079 – 1098.
    [Google Scholar]
  8. Frey, P., Borouchaki, H. and George, P.
    [1998] 3D Delaunay mesh generation coupled with an advancing-front approach.Computer Methods in Applied Mechanics and Engineering, 157(1–2), 115–131.
    [Google Scholar]
  9. Friis, H., Edwards, M. 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, 1192–1220.
    [Google Scholar]
  10. Fung, L., Ding, X. and Dogru, A.
    [2013] An Unstructured Gridding Method for Densely-Spaced Complex Wells in Full-Field Reservoir Simulation.SPE Reservoir Simulation Symposium, Society of Petroleum Engineers, 995–1007.
    [Google Scholar]
  11. George, P. and Borouchaki, H.
    [1998] Delaunay Triangulation and Meshing Applications to Finite Elements, chap. Delaunay Triangulation. 37–71.
    [Google Scholar]
  12. George, P., Borouchaki, H., Frey, P., Laug, P. and Saltel, E.
    [2004] Mesh Generation and Mesh Adaptivity: Theory and Techniques.John Wiley & Sons, Ltd.
    [Google Scholar]
  13. 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]
  14. Kocabas, I. and Maier, F.
    [2013] Analytical and Numerical Modeling of Tracer Flow in Oil Reservoirs Containing High Permeability Streaks., In: ,. Middle East Oil and Gas Show and Conference(SPE),Manama, Bahrain.
    [Google Scholar]
  15. Loner, 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]
  16. Manzoor, S., Edwards, M., Dogru, A. and Al-Shaalan, T.
    [2016] 3D Geological Feature Honored Cell-centered and Vertex-centered Unstructured Grid Generation, and CVD-MPFA Performance. In: EC-MOR XV - 15th European Conference on the Mathematics of Oil Recovery.
    [Google Scholar]
  17. Manzoor, S., Edwards, M., Dogru, A. and Shaalan, T.
    [2015] Boundary Aligned Grid Generation in Three Dimensions and CVD-MPFA Discretization. In: Society of Petroleum Engineers (SPE-RSS).
    [Google Scholar]
  18. Manzoor, S., Edwards, M.G., Dogru, A.H. and Al-Shaalan, T.M.
    [2018] Interior boundary-aligned unstructured grid generation and cell-centered versus vertex-centered CVD-MPFA performance.Computational Geosciences, 22(1), 195–230.
    [Google Scholar]
  19. Marcum, D. and Weatherill, N.
    [1995] Unstructured grid generation using iterative point insertion and local reconnection.AIAA Journal, 33, 1619–1625.
    [Google Scholar]
  20. Mavriplis, D.
    [1995] An Advancing Front Delaunay Triangulation Algorithm Designed for Robustness.Journal of Computational Physics, 117, 90–101.
    [Google Scholar]
  21. Mlacnik, M., Durlofsky, L. and Heinemann, Z.
    [2007] Sequentially Adapted Flow-Based PEBI Grids for Reservoir Simulation.Society of Petroleum Engineers, 11(3), 317–327.
    [Google Scholar]
  22. Müller, J., Roe, P. 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]
  23. Pal, M., Edwards, M. and Lamb, A.
    [2006] Convergence study of a family of flux-continuous finite-volume schemes for the general tensor pressure equation.Int. J. Numer Meth, Fluids, 51(1177–1203).
    [Google Scholar]
  24. Palagi, C. and Aziz, K.
    [2007] Use of Voronoi Grid in Reservoir Simulation.Society of Petroleum Engineers, 2(2), 69–77.
    [Google Scholar]
  25. Ponting, D.
    [1989] Corner Point Grid Geometry in Reservoir Simulation. In: ECMOR I - 1st European Conference on the Mathematics of Oil Recovery.
    [Google Scholar]
  26. Renka, R.
    [2015] Mesh improvement by minimizing a weighted sum of squared element volumes.Int. J. Numer. Meth. Engng, 101, 870–886.
    [Google Scholar]
  27. 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, 1826–1845.
    [Google Scholar]
  28. Watson, D.
    [1981] Computing the n-dimensional Delaunay tessellation with application to Voronoi poly-topes.The Computer Journal, 24(2), 167–172.
    [Google Scholar]
  29. Wen, X., Durlofsky, L. and Edwards, M.
    [2003] Upscaling of Channel Systems in Two Dimensions Using Flow Based Grids.Transport in Porous Media, 51 (3), 343–366.
    [Google Scholar]
http://instance.metastore.ingenta.com/content/papers/10.3997/2214-4609.201802273
Loading
/content/papers/10.3997/2214-4609.201802273
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