Injection-rate histories and well placement influence fault yielding during CO2 geological storage
Abstract
Safe CO ₂ storage requires injection designs that manage both when faults yield and how yielding spreads, yet temporal schedules and well layouts are often evaluated separately. We developed a two-dimensional quasi-static hydromechanical model of immiscible CO ₂ –brine flow around an offset normal fault and compared constant and variable-rate histories, well-placement configurations, and well–fault distances over 30 d. Coulomb failure function (CFF) threshold attainment with frictional weakening defined modeled onset; failure-area ratio and maximum shear plastic strain served as within-mesh comparative indicators of domain-scale localization. The simulations reveal a hierarchy of controls: injection rate primarily governs pressure buildup and onset timing; rate history reorganizes transient stress paths and localization; and well placement and distance redistribute pressure and yielding across the fault. Higher rates increased the within-mesh failure-area ratio, whereas variable histories produced distinct responses that cannot be attributed to cycling alone because peak rate and shut-in transitions were confounded. Local threshold indicators and domain-scale metrics were complementary rather than interchangeable: temporal onset, spatial extent and deformation intensity were distinct response dimensions. These findings motivate joint optimization of injection schedules and well layouts against both temporal onset and spatial localization, while providing model-scale guidance for early quasi-static risk screening rather than predictions of seismicity, dynamic rupture or long-term stability.
Related
- An AI-assisted conditioning and geological interpretation workflow for usage in implicit geological modeling
- Learning to Guide Ontology-Driven Geological Interpretation with Graph Neural Networks
- Emirates Exploration Imager (EXI) Reveals Crater and Boulder Distribution on Deimos: Insights Into Its Dynamical and Geological Evolution
- Escaping Geological Generative Prior Distributions in Bayesian Travel Time Tomography
- Two-stage assessment of regional geological hazard susceptibility under mining disturbance
BibTeX
@article{cui2026injectionrate,
title = {{Injection-rate histories and well placement influence fault yielding during CO2 geological storage}},
author = {Zhendong Cui and Dongxu Liu and Shouding Li and Xiao Li and Ruidong Peng and Xiaojie Jiao and Jiayi Li},
journal = {Scientific Reports},
year = {2026},
month = oct,
doi = {10.1038/s41598-026-74920-7},
url = {https://www.nature.com/articles/s41598-026-74920-7},
note = {CC BY 4.0},
}
RIS
TY - JOUR
TI - Injection-rate histories and well placement influence fault yielding during CO2 geological storage
AU - Zhendong Cui
AU - Dongxu Liu
AU - Shouding Li
AU - Xiao Li
AU - Ruidong Peng
AU - Xiaojie Jiao
AU - Jiayi Li
PY - 2026
DA - 2026/10/08
JO - Scientific Reports
DO - 10.1038/s41598-026-74920-7
UR - https://www.nature.com/articles/s41598-026-74920-7
AB - Safe CO ₂ storage requires injection designs that manage both when faults yield and how yielding spreads, yet temporal schedules and well layouts are often evaluated separately. We developed a two-dimensional quasi-static hydromechanical model of immiscible CO ₂ –brine flow around an offset normal fault and compared constant and variable-rate histories, well-placement configurations, and well–fault distances over 30 d. Coulomb failure function (CFF) threshold attainment with frictional weakening defined modeled onset; failure-area ratio and maximum shear plastic strain served as within-mesh comparative indicators of domain-scale localization. The simulations reveal a hierarchy of controls: injection rate primarily governs pressure buildup and onset timing; rate history reorganizes transient stress paths and localization; and well placement and distance redistribute pressure and yielding across the fault. Higher rates increased the within-mesh failure-area ratio, whereas variable histories produced distinct responses that cannot be attributed to cycling alone because peak rate and shut-in transitions were confounded. Local threshold indicators and domain-scale metrics were complementary rather than interchangeable: temporal onset, spatial extent and deformation intensity were distinct response dimensions. These findings motivate joint optimization of injection schedules and well layouts against both temporal onset and spatial localization, while providing model-scale guidance for early quasi-static risk screening rather than predictions of seismicity, dynamic rupture or long-term stability.
SN - 2045-2322
ER -