Reservoir simulation predominately uses structured grids, but they have difficulties for field-scale simulation with maximum reservoir contact (MRC) wells. In this work, the grid resolution for numerical convergence using structured grids is studied. Results showed that sufficiently fine grid is needed to obtain converged solution. In common practice, the geocellular models are simply upscaled and used in history matching. This study highlights the need for grid resolution to resolve flow dynamics in reservoir simulation and improve well inflow performance calculation. Otherwise, the near-well flow may not have converged and the fidelity of the models for performance prediction is in question.

Improving the fidelity of simulation results by adequately modeling the complex multiphasic flow dynamics near complex wells in full-field simulation is of paramount importance nowadays. This work introduces a full-field unstructured gridding method which can optimally place unstructured grid cells where the resolution is needed. The method produces a consistent discretization that is efficient to compute by using a parallel unstructured reservoir simulator. For a giant Middle-East carbonate reservoir that was developed primarily using complex MRC wells, unstructured grid models were used to improve the accuracy of near-well flow and to better represent well inflow performances. The unstructured grid models are compared against the original structured grid models that show computational cost saving and require few grid cells.

Simulation results demonstrate an unstructured workflow that is practical and can be used to validate near-well modeling accuracy for the existing structured grid simulation results. The method is particular attractive for situation with densely-spaced complex wells where a structured local grid refinement (LGR) method will be ineffective. An unstructured grid is well suited to honor the near-well flow geometry and to focus grid resolution where it is needed. Perpendicular bisection (PEBI) grids are orthogonal by construction. This reduces computational complexity because two-point flux approximation (TPFA) can be applied. In field-scale simulation, this results in a significant improvement to accuracy and computational cost savings.


The current state-of-the-art simulation practice in the oil industry is to use structured grids (corner-point grids or variable-depth-thickness Cartesian grids). While unstructured grid methods such as perpendicular bisection1 (PEBI) grids or Control-Volume Finite-Element2 (CVFE) grids have been active research areas since the late eighties and commercial packages are in service, field-scale simulation applications of unstructured grids are rare. Some of the reasons cited for this include:

  1. the lack of an effective, easy-to-use full-field complex well gridding tool;

  2. the lack of supporting unstructured workflow for full-cycle reservoir simulation;

  3. the cost of unstructured grid simulation; and

  4. the availability of post analysis and visualization tools for unstructured grid simulation.

This content is only available via PDF.
You can access this article if you purchase or spend a download.