Full text loading...
-
Coupling Vertex Centered Based Flow Elements With Poromechanical Finite Elements Using Unstructured Grids
- Publisher: European Association of Geoscientists & Engineers
- Source: Conference Proceedings, ECMOR XVI - 16th European Conference on the Mathematics of Oil Recovery, Sep 2018, Volume 2018, p.1 - 23
Abstract
Significant work has been done for generating unstructured grids. Coupled geomechanics simulation and hydraulic fracture flow for gas shale simulation have given a new impulse for unstructured gridding.
The objective is coupling flow and geomechanics using unstructured grid models, demonstrate the ability to apply a more efficient pressure coupling using discretization at vertices on both sides (mechanical and flow equations).
Coupled equations are discretized and solved on an unstructured flow grid and a geomechanical finite element grid which are composed of various types of elements such as tetrahedrons and hexahedrons. On the flow side, a recent multi-point flux method, the Vertex Approximate Gradient (VAG) (Eymard 2012) is investigated for solving the reservoir equations on such unstructured grid. SPEJ paper 173309 (Samier, Masson Apr. 2017) presented the implementation of VAG scheme inside a next generation reservoir simulator designed for handling unstructured grids. This paper proposes an iterative coupling scheme with full pressure coupling at vertices. The geo-mechanics equations fully coupled to a single phase flow are solved using global pressure.
Then the resulting deformations are iteratively coupled to the multi-phase flow simulator.
Since most geo-mechanic simulators propose fully coupled single phase flow features, the main advantage of this method is the ability to use a full pressure coupling method with industrial simulators.
The convergence of this new scheme is discussed and results are presented for two cases described below. The first case is a validation case used by other SPE papers. The second case is a synthetic faulted and stress sensitive black oil reservoir simulation. Faults are modeled using specific cohesive elements. Material properties are nonlinear according to a Camclay elastoplastic model. Results are compared with standard loose iterative coupling method using TPFA cell centered elements.