Effects of reservoir mechanical properties on induced seismicity during subsurface hydrogen storage.

The intermittent storage of hydrogen in subsurface porous media such as depleted gas fields could be pivotal to a successful energy transition. Numerical simulations investigate the intermittent storage of hydrogen in a porous, depleted subsurface reservoir. Various parametric studies are performed to assess the effect of mechanical properties of the reservoir (i.e. Young's modulus, Poisson's ratio, Biot coefficient and permeability) on the induced fault slip of a single through-going fault that transverses the entire reservoir. Simulations are run using a three-dimensional, finite element, fully coupled hydromechanical code with explicit representations of layers and faults. The effect of the domain mesh refinement and fault mesh refinement on the fault slip versus operation time solution is investigated. The fault is observed to slip in two distinct events, one during the second injection period and one in the third injection period. The fault is not observed to slip during the storage or withdrawal periods. It is found that in order to minimize seismic risk, a reservoir rock with high Young's modulus (>40 GPa), high Poisson's ratio (>0.30) and high Biot coefficient (>0.65) would be preferable for hydrogen storage. Reservoir rocks of low Young's modulus (10-30 GPa), intermediate Poisson's ratio (0.00-0.30) and low-to-intermediate Biot coefficient (0.25-0.65), at high injection rates, were found to have higher potential of inducing large seismic events.This article is part of the theme issue 'Induced seismicity in coupled subsurface systems'.

How tough is brittle bone? Investigating osteogenesis imperfecta in mouse bone.

The multiscale hierarchical structure of bone is naturally optimized to resist fractures. In osteogenesis imperfecta, or brittle bone disease, genetic mutations affect the quality and/or quantity of collagen, dramatically increasing bone fracture risk. Here we reveal how the collagen defect results in bone fragility in a mouse model of osteogenesis imperfecta (oim), which has homotrimeric α1(I) collagen. At the molecular level, we attribute the loss in toughness to a decrease in the stabilizing enzymatic cross-links and an increase in nonenzymatic cross-links, which may break prematurely, inhibiting plasticity. At the tissue level, high vascular canal density reduces the stable crack growth, and extensive woven bone limits the crack-deflection toughening during crack growth. This demonstrates how modifications at the bone molecular level have ramifications at larger length scales affecting the overall mechanical integrity of the bone; thus, treatment strategies have to address multiscale properties in order to regain bone toughness. In this regard, findings from the heterozygous oim bone, where defective as well as normal collagen are present, suggest that increasing the quantity of healthy collagen in these bones helps to recover toughness at the multiple length scales.

Role of geomechanically grown fractures on dispersive transport in heterogeneous geological formations.

A second order in space accurate implicit scheme for time-dependent advection-dispersion equations and a discrete fracture propagation model are employed to model solute transport in porous media. We study the impact of the fractures on mass transport and dispersion. To model flow and transport, pressure and transport equations are integrated using a finite-element, node-centered finite-volume approach. Fracture geometries are incrementally developed from a random distributions of material flaws using an adoptive geomechanical finite-element model that also produces fracture aperture distributions. This quasistatic propagation assumes a linear elastic rock matrix, and crack propagation is governed by a subcritical crack growth failure criterion. Fracture propagation, intersection, and closure are handled geometrically. The flow and transport simulations are separately conducted for a range of fracture densities that are generated by the geomechanical finite-element model. These computations show that the most influential parameters for solute transport in fractured porous media are as follows: fracture density and fracture-matrix flux ratio that is influenced by matrix permeability. Using an equivalent fracture aperture size, computed on the basis of equivalent permeability of the system, we also obtain an acceptable prediction of the macrodispersion of poorly interconnected fracture networks. The results hold for fractures at relatively low density.