Building PEBI Grids Honoring 3D Geological Features And Flow Patterns Using Centroidal Voronoi Tessellations

in: Proc. 32nd Gocad Meeting, Nancy, pages 12, ASGA

Abstract

In this paper, several constraints are put together to build unstructured reservoir grids for flow simulation: (1) refinement in high flow/vorticity areas, (2) cell orientation according to flow pattern, (3) well conformity, (4) faults and horizons conformity, (5) cell aspect ratio. The generated grids are fully PErpandicular BIsector (Voronoi) and globally honor each of these constraints. An iterative process minimizes a function that is a weighted sum of the constraints. As this function is not proven C 1 , an adaptive steepest gradient solver is used. Flow simulation properties are upscaled in the reservoir grid from the geological grid using a global upscaling technique. A fast implementation is proposed, that takes advantage of the Voronoi pattern. The methodology is applied in 3D. 2D flow simulation results using a Two-point-Flux Approximation scheme are presented.

Download / Links

BibTeX Reference

@INPROCEEDINGS{,
    author = { Merland, Romain and Levy, Bruno and Caumon, Guillaume },
     title = { Building PEBI Grids Honoring 3D Geological Features And Flow Patterns Using Centroidal Voronoi Tessellations },
     month = { "sep" },
 booktitle = { Proc. 32nd Gocad Meeting, Nancy },
      year = { 2012 },
     pages = { 12 },
 publisher = { ASGA },
  abstract = { In this paper, several constraints are put together to build unstructured reservoir grids for flow simulation: (1) refinement in high flow/vorticity areas, (2) cell orientation according to flow
pattern, (3) well conformity, (4) faults and horizons conformity, (5) cell aspect ratio. The generated
grids are fully PErpandicular BIsector (Voronoi) and globally honor each of these constraints. An
iterative process minimizes a function that is a weighted sum of the constraints. As this function
is not proven C 1 , an adaptive steepest gradient solver is used. Flow simulation properties are
upscaled in the reservoir grid from the geological grid using a global upscaling technique. A
fast implementation is proposed, that takes advantage of the Voronoi pattern. The methodology
is applied in 3D. 2D flow simulation results using a Two-point-Flux Approximation scheme are
presented. }
}