1887

Abstract

Numerical simulations of two-phase Darcy flows in heterogeneous porous media requires choosing an appropriate set of primary unknowns, which may be challenging, especially when dealing with very flat capillary pressure curves, dry regions or saturation jumps at the rock type interfaces. The classical approaches fail to cope with all of those difficulties. In particular the two-pressure formulation allows handling saturation jumps, but beaks down if the capillary pressure doesn’t depend on saturation. It also lacks robustness when dealing with nearly residual water saturations. On the other hand, for homogeneous medium, the pressure – saturation formulation is known to be robust when dealing with dry media and can handle vanishing or constant capillary pressure curves. Unfortunately it is not always possible to extend it to the case of discontinuous capillary pressure curves. In this paper, a new formulation based on parametrization techniques for the capillary pressure monotone graph extension is proposed which handles all the above mentioned difficulties while still using only two unknowns by degree of freedom. We illustrate the efficiency of our approach by numerous numerical experiments dealing with water gas flow in fractured tight gas reservoirs using the data set presented in [2]. Following [1], the fractures are modelled as interfaces of codimension one with continuous pressure at the matrix fracture interfaces. During the injection phase of the simulation, water penetrates only a few tens of centimetres deep in the matrix rock. Therefore, in order to obtain an accurate numerical approximation, an anisotropic refinement of the mesh is used in the neighbourhood of the fractures using prismatic elements. The connection with the surrounding tetrahedral mesh in the matrix domain is achieved using pyramids. Following [1], the model is discretized using the Vertex Approximate Gradient scheme which allows for general polyhedral cells. The numerical performance of the new approach is evaluated for various choices of capillary pressures curves. The comparison with classical formulations shows that the new approach is more efficient both in terms of Newton iterations and CPU time. [1] K.Brenner, M.Groza, C.Guichard, R.Masson: Vertex Approximate Gradient Scheme for Hybrid Dimensional Two-Phase Darcy Flows in Fractured Porous Media, M2AN, pp.49 2, 303-330, 2015. [2] D.Y.Ding, H.Langouet, L.Jeannin: Simulation of Fracturing Induced Formation Damage and Gas from Fractured Wells in Tight Gas Reservoirs, SPE 153255, 2012.

Loading

Article metrics loading...

/content/papers/10.3997/2214-4609.201601787
2016-08-29
2024-03-28
Loading full text...

Full text loading...

http://instance.metastore.ingenta.com/content/papers/10.3997/2214-4609.201601787
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