Luke P. Frash

Understanding hydraulic fracturing: a multi-scale problem

Despite the impact that hydraulic fracturing has had on the energy sector, the physical mechanisms that control its efficiency and environmental impacts remain poorly understood in part because the length scales involved range from nanometres to kilometres. We characterize flow and transport in shale formations across and between these scales using integrated computational, theoretical and experimental efforts/methods. At the field scale, we use discrete fracture network modelling to simulate production of a hydraulically fractured well from a fracture network that is based on the site characterization of a shale gas reservoir. At the core scale, we use triaxial fracture experiments and a finite-discrete element model to study dynamic fracture/crack propagation in low permeability shale. We use lattice Boltzmann pore-scale simulations and microfluidic experiments in both synthetic and shale rock micromodels to study pore-scale flow and transport phenomena, including multi-phase flow and fluids mixing. A mechanistic description and integration of these multiple scales is required for accurate predictions of production and the eventual optimization of hydrocarbon extraction from unconventional reservoirs. Finally, we discuss the potential of CO2 as an alternative working fluid, both in fracturing and re-stimulating activities, beyond its environmental advantages. This article is part of the themed issue ‘Energy and the subsurface’.

High-stress triaxial direct-shear fracturing of Utica shale and in situ X-ray microtomography with permeability measurement: SHALE FRACTURE, μCT, AND PERMEABILITY

The challenge of characterizing subsurface fluid flow has motivated extensive laboratory studies, yet fluid-flow through rock specimens in which fractures are created and maintained at high-stress conditions remains under-investigated at this time. The studies of this type that do exist do not include in situ fracture geometry measurements acquired at stressed conditions, which would be beneficial for interpreting the flow behavior. Therefore, this study investigates the apparent permeability induced by direct-shear fracture stimulation through Utica shale (a shale gas resource and potential caprock material) at high triaxial-stress confinement and for the first time relates these values to simultaneously acquired in situ X-ray radiography and microtomography images. Change in fracture geometry and apparent permeability was also investigated at additional reduced stress states. Finite element and combined finite discrete element modeling were used to evaluate the in situ observed fracturing process. Results from this study indicate that the increase in apparent permeability through fractures created at high-stress (22.2 MPa) was minimal relative to the intact rock (less than 1 order of magnitude increase) while fractures created at low stress (3.4 MPa) were significantly more permeable (2 to 4 orders of magnitude increase). This study demonstrates the benefit of in situ X-ray observation coupled with apparent permeability measurement to analyze fracture creation in the subsurface. Our results show that the permeability induced by fractures through shale at high stress can be minor and therefore favorable in application to CO2 sequestration caprock integrity but unfavorable for hydrocarbon recovery from unconventional reservoirs.

A non-locking composite tetrahedron element for the combined finite discrete element method

Purpose In order to avoid the problem of volumetric locking often encountered when using constant strain tetrahedral finite elements, the purpose of this paper is to present a new composite tetrahedron element which is especially designed for the combined finite-discrete element method (FDEM). Design/methodology/approach A ten-noded composite tetrahedral (COMPTet) finite element, composed of eight four-noded low order tetrahedrons, has been implemented based on Munjiza’s multiplicative decomposition approach. This approach naturally decomposes deformation into translation, rotation, plastic stretches, elastic stretches, volumetric stretches, shear stretches, etc. The problem of volumetric locking is avoided via a selective integration approach that allows for different constitutive components to be evaluated at different integration points. Findings A number of validation cases considering different loading and boundary conditions and different materials for the proposed element are presented. A practical application of the use of the COMPTet finite element is presented by quantitative comparison of numerical model results against simple theoretical estimates and results from acrylic fracturing experiments. All of these examples clearly show the capability of the composite element in eliminating volumetric locking. Originality/value For this tetrahedral element, the combination of “composite” and “low order sub-element” properties are good choices for FDEM dynamic fracture propagation simulations: in order to eliminate the volumetric locking, only the information from the sub-elements of the composite element are needed which is especially convenient for cases where re-meshing is necessary, and the low order sub-elements will enable robust contact interaction algorithms, which maintains both relatively high computational efficiency and accuracy.

Flow-through compression cell for small-angle and ultra-small-angle neutron scattering measurements

In situ measurements of geological materials under compression and with hydrostatic fluid pressure are important in understanding their behavior under field conditions, which in turn provides critical information for application-driven research. In particular, understanding the role of nano- to micro-scale porosity in the subsurface liquid and gas flow is critical for the high-fidelity characterization of the transport and more efficient extraction of the associated energy resources. In other applications, where parts are produced by the consolidation of powders by compression, the resulting porosity and crystallite orientation (texture) may affect its in-use characteristics. Small-angle neutron scattering (SANS) and ultra SANS are ideal probes for characterization of these porous structures over the nano to micro length scales. Here we show the design, realization, and performance of a novel neutron scattering sample environment, a specially designed compression cell, which provides compressive stress and hydrostatic pressures with effective stress up to 60 MPa, using the neutron beam to probe the effects of stress vectors parallel to the neutron beam. We demonstrate that the neutron optics is suitable for the experimental objectives and that the system is highly stable to the stress and pressure conditions of the measurements.