XFVM modelling of fracture aperture induced by shear and tensile opening.

In reservoir simulation, it is important to understand the mechanical behaviour of fractured rocks and the effect of shear and tensile displacements of fractures on their aperture. Tensile opening directly enhances the fracture aperture, whereas shear of a preexisting rough-walled fracture creates aperture changes dependent on the local stress state. Since fracture dilatation increases reservoir permeability, both processes must be included in a realistic and consistent manner into the mechanical reservoir simulation model. Here, we use the extended finite volume method (XFVM) to conduct flow and geomechanics simulations. In XFVM, fractures are embedded in a poroelastic matrix and are modelled with discontinuous basis functions. On each fracture segment the tractions and compressive forces are calculated, and one extra degree of freedom is added for both the shear and tensile displacement. In this particular XFVM implementation we assume that linear elasticity and steady state fluid pressure adequately constrain the effective stress. In this paper, shear dilation is not calculated a posteriori, but it enters the equations such that aperture changes directly affect the stress state. This is accomplished by adding shear dilation to the displacement gradients and therefore ascertains a consistent representation in the stress-strain relations and force balances. We illustrate and discuss the influence of this extra term in two simple test cases and in a realistic layer-restricted two-dimensional fracture network subjected to plausible in situ stress and pore pressure conditions.

Numerical study of the effects of particle shape and polydispersity on permeability.

We study through numerical simulations the dependence of the hydraulic permeability of granular materials on the particle shape and the grain size distribution. Several models of sand are constructed by simulating the settling under gravity of the grains; the friction coefficient is varied to construct packs of different porosity. The size distribution and shapes of the grains mimic real sands. Fluid flow is simulated in the resulting packs using a finite element method and the permeability of the packs is successfully compared with available experimental data. Packs of nonspherical particles are less permeable than sphere packs of the same porosity. Our results indicate that the details of grain shape and size distribution have only a small effect on the permeabilty of the systems studied.

Transitions between nanomechanical and continuum mechanical contacts: new insights from liquid structure.

The use of continuum mechanics to describe contacts involving nanoscale and atomic interactions has been one of the key controversies in nanoscience, tribology, and petrophysical and geological studies. By applying a novel nonequilibrium molecular dynamics scheme to wet quartz contacts, this study revealed the key transitions between continuum electrostatic, nanomechanical and Hertzian contact behaviors at around one nm of surface separation, which results in critical contact pressure fluctuations between -30 and 100 MPa. Using a novel liquid-structure analysis scheme based on the spatial distribution of water molecules, the nanomechanical behavior was found to originate from the collapse and localization of layers of water molecules. Moreover, the role of surface curvature on this effect was also quantified and explained based on a new topological descriptor. The findings of this study enrich our understanding of wet contacts and have a wide range of applications from the nanoscale to macroscale.