Search

Underground gas storage is a versatile tool for managing energy resources and addressing pressing environmental concerns. While natural gas is stored in geological formations since the beginning of the 20th century, hydrogen has recently been considered as a potential candidate toward a more flexible and sustainable energy infrastructure. Furthermore, these formations can also be used to sequester environmentally harmful gases such as CO2. When such operations are implemented in faulted basins, however, safety concerns may arise due to the possible reactivation of pre-existing faults, which could result in (micro)-seismicity events. In the Netherlands, it has been recently noted that fault reactivation can occur "unexpectedly" during the life of an underground gas storage (UGS) site, even when stress conditions are not expected to cause a failure. The present two-part work aims to develop a modeling framework to investigate the physical mechanisms causing such occurrences and define a safe operational bandwidth for pore pressure variation for UGS operations in the faulted reservoirs of the Rotliegend formation, the Netherlands. In this follow-up paper, we investigate in detail the mechanisms and crucial factors that result in fault reactivation at various stages of a UGS. The mathematical and numerical model described in Part I is used, also accounting for the effect of geochemical dissolution on reservoir and caprock weakening. TThe study investigates the risks of fault activation caused by the storage of different fluids for various purposes, such as long-term CO2 sequestration, CH4 and N2 injection and extraction cycles, and N2 permanent storage. The results show how geomechanical properties and reservoir operating conditions may increase the risk of fault reactivation at various UGS stages. Finally, operational guidelines for improving secure storage operations are presented., Comment: 31 pages, 24 figures, 4 tables

Underground gas storage (UGS) is a worldwide well-established technology that is becoming even more important to cope with seasonal peaks of gas consumption due to the growing uncertainties of the energy market. Safety issues concerning the reactivation of pre-existing faults might arise if the target reservoir is located in a faulted basin, where human activities can trigger (micro-)seismicity events. In the Netherlands, it has been observed that fault activation can occur somehow "unexpectedly" after the primary production (PP), i.e., during cushion gas injection (CGI) and UGS cycles, when the stress regime should be in the unloading/reloading path. To understand the physical mechanisms responsible for such occurrences, a 3D mathematical model coupling frictional contact mechanics in faulted porous rocks with fluid flow is developed, implemented and tested. The final aim of this two-part work is to define a safe operational bandwidth for the pore pressure range for UGS activities in the faulted reservoirs of the Rotliegend formation. Part I of this work concerns the development of the mathematical and numerical model of frictional contact mechanics and flow in faulted porous rocks. A mixed discretization of the governing PDEs under frictional contact constraints along the faults is used. A slip-weakening constitutive law governing the fault macroscopic behavior is also presented. The model is tested in the setting of an ideal reservoir located in the Rotliegend formation. The analyses point out how fault reactivation during PP can lead to a stress redistribution, giving rise to a new equilibrium configuration. When the fault is reloaded in the opposite direction during the CGI and/or UGS stages, further activation events can occur even if the stress range does not exceed either the undisturbed initial value or the maximum strength ever experienced by the formation., Comment: 22 pages, 20 figures

Understanding the heterogeneity arising from the complex architecture of sedimentary sequences in alluvial fans is challenging. This paper develops a statistical inverse framework in a multi-zone transition probability approach for characterizing the heterogeneity in alluvial fans. An analytical solution of the transition probability matrix is used to define the statistical relationships among different hydrofacies and their mean lengths, integral scales, and volumetric proportions. A statistical inversion is conducted to identify the multi-zone transition probability models and estimate the optimal statistical parameters using the modified Gauss-Newton-Levenberg-Marquardt method. The Jacobian matrix is computed by the sensitivity equation method, which results in an accurate inverse solution with quantification of parameter uncertainty. We use the Chaobai River alluvial fan in the Beijing Plain, China, as an example for elucidating the methodology of alluvial fan characterization. The alluvial fan is divided into three sediment zones. In each zone, the explicit mathematical formulations of the transition probability models are constructed with optimized different integral scales and volumetric proportions. The hydrofacies distributions in the three zones are simulated sequentially by the multi-zone transition probability-based indicator simulations. The result of this study provides the heterogeneous structure of the alluvial fan for further study of flow and transport simulations., Comment: 29 pages, 7 figures, and 3 tables