Parallel numerical modeling of hybrid-dimensional compositional non-isothermal Darcy flows in fractured porous media

This paper introduces a new discrete fracture model accounting for non-isothermal compositional multiphase Darcy flows and complex networks of fractures with intersecting, immersed and non immersed fractures. The so called hybrid-dimensional model using a 2D model in the fractures coupled with a 3D model in the matrix is first derived rigorously starting from the equi-dimensional matrix fracture model. Then, it is dis-cretized using a fully implicit time integration combined with the Vertex Approximate Gradient (VAG) finite volume scheme which is adapted to polyhedral meshes and anisotropic heterogeneous media. The fully coupled systems are assembled and solved in parallel using the Single Program Multiple Data (SPMD) paradigm with one layer of ghost cells. This strategy allows for a local assembly of the discrete systems. An efficient preconditioner is implemented to solve the linear systems at each time step and each Newton type iteration of the simulation. The numerical efficiency of our approach is assessed on different meshes, fracture networks, and physical settings in terms of parallel scalability, nonlinear convergence and linear convergence

An exploration of the IGA method for efficient reservoir simulation

Novel numerical methods present exciting opportunities to improve the efficiency of reservoir simulators. Because potentially significant gains to computational speed and accuracy may be obtained, it is worthwhile explore alternative computational algorithms for both general and case-by-case application to the discretization of the equations of porous media flow, fluid-structure interaction, and/or production. In the present work, the fairly new concept of isogeometric analysis (IGA) is evaluated for its suitability to reservoir simulation via direct comparison with the industry standard finite difference (FD) method and 1st order standard finite element method (SFEM). To this end, two main studies are carried out to observe IGA’s performance with regards to geometrical modeling and ability to capture steep saturation fronts. The first study explores IGA’s ability to model complex reservoir geometries, observing L2 error convergence rates under a variety of refinement schemes. The numerical experimental setup includes an 'S' shaped line sink of varying curvature from which water is produced in a 2D homogenous domain. The accompanying study simplifies the domain to 1D, but adds in multiphase physics that traditionally introduce difficulties associated with modeling of a moving saturation front. Results overall demonstrate promise for the IGA method to be a particularly effective tool in handling geometrically difficult features while also managing typically challenging numerical phenomena.Petroleum and Geosystems Engineerin

Multiscale Finite Element Methods for Nonlinear Problems and their Applications

In this paper we propose a generalization of multiscale finite element methods (Ms-FEM) to nonlinear problems. We study the convergence of the proposed method for nonlinear elliptic equations and propose an oversampling technique. Numerical examples demonstrate that the over-sampling technique greatly reduces the error. The application of MsFEM to porous media flows is considered. Finally, we describe further generalizations of MsFEM to nonlinear time-dependent equations and discuss the convergence of the method for various kinds of heterogeneities

Effects of finite strains in fully coupled 3D geomechanical simulations

Numerical modeling of geomechanical phenomena and geo-engineering problems often involves complex issues related to several variables and corresponding coupling effects. Under certain circumstances, both soil and rock may experience a nonlinear material response caused by, for example, plastic, viscous, or damage behavior or even a nonlinear geometric response due to large deformations or displacements of the solid. Furthermore, the presence of one or more fluids (water, oil, gas, etc.) within the skeleton must be taken into account when evaluating the interaction between the different phases of the continuum body. A multiphase three-dimensional (3D) coupled model of finite strains, suitable for dealing with solid-displacement and fluid-diffusion problems, is described for assumed elastoplastic behavior of the solid phase. Particularly, a 3D mixed finite element was implemented to fulfill stability requirements of the adopted formulation, and a permeability tensor dependent on deformation is introduced. A consolidation scenario induced by silo filling was investigated, and the effects of the adoption of finite strains are discusse

Finite element modelling of interaction between surface and Darcy flow regimes through soils

The present work deals with the impact of surface flow on hydrodynamic conditions in saturated underground domains. A three dimensional finite element analysis of water flow has been used to obtain the required simulations. The results clearly show the effects of the surface flow on the hydrodynamic conditions of the subsurface porous regions. This analysis is an important prerequisite for the prediction of contaminant mobility in soils and hence provides a convenient tool for the prediction of environmentally important subsurface flow processes. For low permeability cases, considered here, governing equations consist of water continuity and Darcy equations. These equations are solved using a robust and reliable finite element procedure