Multi-physics and multi-scale characterization of shale anisotropy

Shales are the most abundant sedimentary rock type in the Earth's shallow crust. In the past decade or so, they have attracted increased attention from the petroleum industry as reservoirs, as well as more traditionally for their sealing capacity for hydrocarbon/CO2 traps or underground waste repositories. The effectiveness of both fundamental and applied shale research is currently limited by (i) the extreme variability of physical, mechanical and chemical properties observed for these rocks, and by (ii) the scarce data currently available. The variability in observed properties is poorly understood due to many factors that are often irrelevant for other sedimentary rocks. The relationships between these properties and the petrophysical measurements performed at the field and laboratory scales are not straightforward, translating to a scale dependency typical of shale behaviour. In addition, the complex and often anisotropic micro-/meso-structures of shales give rise to a directional dependency of some of the measured physical properties that are tensorial by nature such as permeability or elastic stiffness. Currently, fundamental understanding of the parameters controlling the directional and scale dependency of shale properties is far from complete. Selected results of a multi-physics laboratory investigation of the directional and scale dependency of some critical shale properties are reported. In particular, anisotropic features of shale micro-/meso-structures are related to the directional-dependency of elastic and fluid transport properties: - Micro-/meso-structure (μm to cm scale) characterization by electron microscopy and X-ray tomography; - Estimation of elastic anisotropy parameters on a single specimen using elastic wave propagation (cm scale); - Estimation of the permeability tensor using the steady-state method on orthogonal specimens (cm scale); - Estimation of the low-frequency diffusivity tensor using NMR method on orthogonal specimens (<μm scale). For each of the above properties, leading-edge experimental techniques have been associated with novel interpretation tools. In this contribution, these experimental and interpretation methods are described. Relationships between the measured properties and the corresponding micro-/meso-structural features are discussed. For example, P-wave velocity was measured along 100 different propagation paths on a single cylindrical shale specimen using miniature ultrasonic transducers. Assuming that (i) the elastic tensor of this shale is transversely isotropic; and (i) the sample has been cored perfectly perpendicular to the bedding plane (symmetry plane is horizontal), Thomsen's anisotropy parameters inverted from the measured velocities are: - P-wave velocity along the symmetry axis (perpendicular to the bedding plane) αo=3.45km/s; - P-wave anisotropy ɛ=0.12; - Parameter controlling the wave front geometry δ=0.058. A novel inversion algorithm allows for recovering these parameters without assuming a priori a horizontal bedding (symmetry) plane. The inversion of the same data set using this algorithm yields (i) αo=3.23km/s, ɛ=0.25 and δ=0.18, and (ii) the elastic symmetry axis is inclined of ω=30° with respect to the specimen's axis. Such difference can have strong impact on field applications (AVO, ray tracing, tomography).