Xavier Raynaud

Odd Andersen, Halvor M. Nilsen, Xavier Raynaud

In this paper we study the use of Virtual Element method for geomechanics. Our emphasis is on applications to reservoir simulations. The physical processes that form the reservoirs, such as sedimentation, erosion and faulting, lead to complex geometrical structures. A minimal representation, with respect to the physical parameters of the system, then naturally leads to general polyhedral grids. Numerical methods which can directly handle this representation will be highly favorable, in particular in the setting of advanced work-flows. The Virtual Element method is a promising candidate to solve the linear elasticity equations on such models. In this paper, we investigate some of the limits of the VEM method when used on reservoir models. First, we demonstrate that care must be taken to make the method robust for highly elongated cells, which is common in these applications, and show the importance of calculating forces in terms of traction on the boundary of the elements for elongated distorted cells. Second, we study the effect of triangulations on the surfaces of curved faces, which also naturally occur in subsurface models. We also demonstrate how a more stable regularization term for reservoir application can be derived.

Halvor Møll Nilsen, Idar Larsen, Xavier Raynaud

Simulation of fracturing processes in porous rocks can be divided into two main branches: (i) modeling the rock as a continuum which is enhanced with special features to account for fractures, or (ii) modeling the rock by a discrete (or discontinuous) approach that describes the material directly as a collection of separate blocks or particles, e.g., as in the discrete element method (DEM). In the modified discrete element (MDEM) method, the effective forces between virtual particles are modified in all regions, without failing elements, so that they reproduce the discretization of a first order finite element method (FEM) for linear elasticity. This provides an expression of the virtual forces in terms of general Hook's macro-parameters. Previously, MDEM has been formulated through an analogy with linear elements for FEM. We show the connection between MDEM and the virtual element method (VEM), which is a generalization of FEM to polyhedral grids. Unlike standard FEM, which computes strain-states in a reference space, MDEM and VEM compute stress-states directly in real space. This connection leads us to a new derivation of the MDEM method. Moreover, it gives the basis for coupling (M)DEM to domain with linear elasticity described by polyhedral grids, which makes it easier to apply realistic boundary conditions in hydraulic-fracturing simulations. This approach also makes it possible to combine fine-scale (M)DEM behavior near the fracturing region with linear elasticity on complex reservoir grids in the far-field region without regridding. To demonstrate the simulation of hydraulic fracturing, the coupled (M)DEM-VEM method is implemented using the Matlab Reservoir Simulation Toolbox (MRST) and linked to an industry-standard reservoir simulator.

David Cohen, Xavier Raynaud

We propose a fully discretised numerical scheme for the hyperelastic rod wave equation on the line. The convergence of the method is established. Moreover, the scheme can handle the blow-up of the derivative which naturally occurs for this equation. By using a time splitting integrator which preserves the invariants of the problem, we can also show that the scheme preserves the positivity of the energy density.