Hydrogen (H₂) is an attractive energy carrier to move, store, and deliver energy in a form that can be easily used. Field proven technology for underground hydrogen storage (UHS) is essential for a successful hydrogen economy. Options for this are manmade caverns, salt domes/caverns, saline aquifers, and depleted oil/gas fields, where large quantities of gaseous hydrogen have been stored in caverns for many years. The key requirements intrinsic of a porous rock formation for seasonal storage of hydrogen are: adequate capacity, ability to contain H₂, capability to inject/extract high volumes of H₂, and a reliable caprock to prevent leakage. We have carefully evaluated a commercial non-isothermal compositional gas reservoir simulator and its suitability for hydrogen storage and withdrawal from saline aquifers and depleted oil/gas reservoirs. We have successfully calibrated the gas equation of state model against published laboratory H₂ density and viscosity data as a function of pressure and temperature. Comparisons between the H₂, natural gas and CO₂ storage in real field models were also performed. Our numerical models demonstrated more lateral spread of the H₂ when compared to CO₂ and natural gas with a need for special containment in H₂ projects. It was also observed that the experience with CO₂ and natural gas storage cannot be simply replicated with H₂.

Summary: Hydrogen (H2) is an attractive energy carrier, and its true potential is in decarbonizing industries, such as for providing heat for buildings and being a reliable fuel for trains, buses, and heavy trucks. Industry is already making tremendous progress in cutting costs and improving the efficiency of hydrogen infrastructure. Currently, heating is primarily provided by using natural gas and transportation by gasoline with a large carbon footprint. Hydrogen has a similarly high energy density, but there are technical challenges preventing its large-scale use as an energy carrier. Underground geologic storage of hydrogen in porous media (aquifers and hydrocarbon reservoirs) could offer substantial storage capacity at low cost as well as buffer capacity to meet changing seasonal electricity demands or possible disruptions in power supply. Underground geologic storage must have adequate capacity and ability to inject/extract high volumes with a reliable caprock. A thorough study is essential for a large number of site surveys to locate and fully characterize the subsurface geological storage sites both onshore and offshore. An isothermal compositional reservoir simulator was used to evaluate hydrogen storage and withdrawal from saline aquifers and depleted oil/gas reservoirs. The phase behavior, fluid properties, and petrophysical models were all calibrated against published laboratory data for density, viscosity, relative permeability, and capillary pressure for a given site. History-matched dynamic models of two CO2 injection field projects in saline aquifers and one natural gas storage in a depleted oil reservoir were considered as hypothetical hydrogen seasonal storage sites. A wide range of working gas volume/cushion gas volume ratios was observed, meaning that careful optimization is required for a successful storage project. For the aquifer cases, the range was 0.292 to 1.883 and a range of 1.045 to 4.4 was observed for the depleted hydrocarbon reservoir scenarios. For the saline aquifer cases, a higher injection rate, longer injection/withdrawal (I/W) cycles, and the use of pump wells to control the hydrogen plume spreading were all beneficial for improving the working gas/cushion gas ratio and the working gas volume. Plume control was important for storage in the oil reservoir in which changes in the well length location and orientation showed high sensitivity in the working and cushion gas volumes. Sensitivities to the initial gas saturation in the depleted gas reservoir scenarios suggested that both cushion and working gas volumes increased with the initial gas saturation while the ratio of working to cushion gas volumes decreased with the initial gas saturation. Finally, when comparing the ratios of working to cushion gas volumes, it was the highest for the depleted oil reservoir, followed by the depleted gas reservoir, and the aquifer.

Three different reservoir simulators that utilize both two-phase and three-phase microemulsion phase behavior models are used to model surfactant–polymer flooding to determine and compare their results. Different models are used in each simulator to describe the physical behavior of injected chemicals into the reservoir, which raises the need to benchmark their results. The physical behavior models of polymer and surfactant were constructed and verified on a 1D scale reservoir model and further verified in a 3D model. Finally, simulations were conducted in a field-scale reservoir containing 680,400 grids, where results were compared and analyzed. The 1D and 3D model results suggest an excellent match between the different simulators in modeling surfactant–polymer floods. In the case of the field-scale model, the simulators matched in terms of oil recovery and total volumes produced and injected, while having similar reservoir pressure profiles but with significant discrepancies in terms of injected and produced chemicals. These results indicate that despite the differences in the calculated injected and produced chemicals due to the different models in the simulators, the effect of surfactant–polymer floods on oil recovery, total injected and produced fluids, and average pressure profiles can be comparably modeled in all of the three simulators.

This study develops a conceptual framework for characterizing reservoir architecture in multi-component, discrete systems using pressure transient analysis (PTA), aimed at calibrating inflow geometry prior to full-field dynamic simulation for subsurface gas storage applications such as CO₂ and hydrogen. A secondary objective is to identify variations in permeability over time by analyzing flow capacity trends and evaluating the dynamic influence of faults and fractures. The analysis is based on a gas-condensate field comprising seven wells and four zones (A, B, C, D), using integrated dynamic datasets including extended well tests (EWTs), mud loss, production logs, and production data. Detailed interpretation of PX-1’s EWT indicated delayed re-pressurization and persistent under-pressure, suggesting a compartmentalized or transient system with limited gas-in-place connectivity. Four reservoir architecture concepts were developed: (1) lithology-dominated inflow, (2) structurally controlled inflow, (3) discrete, weakly connected compartments, and (4) transient-dominated systems with tight matrix GIIP. These concepts informed four reservoir models: matrix-only (M), areal heterogeneity (A), sparse bodies (B), and sparse networks (S). Application of these models across other wells revealed consistent localized KH (permeability–thickness product) behavior, with all models fitting short-duration data comparably. However, only sparse drainage models (B/S) adequately matched PX-1’s EWT response. PTA results confirm that well tests constrain KH locally but provide limited insight into large-scale reservoir architecture. EWTs may reach ~1 km, while shorter tests are confined to ~200–400 m, typically within one to two simulation grid blocks. This study demonstrates how integrating PTA with multi-scale data improves characterization of naturally fractured, tight carbonate reservoirs and supports reservoir simulation and history matching for hydrogen storage evaluation. Based on reservoir simulations, this study concluded that naturally fractured carbonate gas reservoirs can provide significant storage and injection capacities for underground hydrogen storage. This study exemplifies how to characterize the naturally fractured tight carbonate reservoirs by integrating multi-scale and multi-dimensional data such as PTA. Furthermore, this study assists in gridding for full-field reservoir models, for history matching and quantifying the potential of hydrogen storage in these complex reservoirs. The proposed workflow provides an uncertainty-bounded reservoir characterization framework and should not be interpreted as a complete field-design methodology for hydrogen storage. The modeling does not explicitly couple geomechanical fracture growth, hydrogen diffusion, long-term geochemical reactions, or caprock integrity degradation. Therefore, the presented storage scenarios represent technically feasible cases under defined assumptions. Comprehensive site-specific geomechanical and containment assessments are required prior to field-scale implementation.

Summary: The study aims to quantitatively assess the risk of hydrate formation within the porous formation and its consequences on injectivity during storage of CO2 in depleted gas reservoirs considering low temperatures caused by the Joule-Thomson (JT) effect and hydrate kinetics. Hydrates formed during CO2 storage operation can occupy porous spaces in the reservoir rock, reducing the rock’s permeability and thus becoming a hindrance to the storage project. The aim was to understand which mechanisms can mitigate or prevent the formation of hydrates. The key mechanisms we studied included water dry-out, heat exchange with surrounding rock formation, and capillary pressure. A semicompositional thermal reservoir simulator is used to model the fluid and heat flow of CO2 through a reservoir initially composed of brine and methane. The simulator can model the formation and dissociation of both methane and CO2 hydrates using kinetic reactions. This approach has the advantage of computing the amount of hydrate deposited and estimating its effects on the porosity and permeability alteration. Sensitivity analyses are also carried out to investigate the impact of different parameters and mechanisms on the deposition of hydrates and the injectivity of CO2. Simulation results for a simplified model were verified with results from the literature. The key results of this work are as follows: (1) The JT effect strongly depends on the reservoir permeability and initial pressure and could lead to the formation of hydrates within the porous media even when the injected CO2 temperature was higher than the hydrate equilibrium temperature; (2) the heat gain from underburden and overburden rock formations could prevent hydrates formed at late time; (3) permeability reduction increased the formation of hydrates due to an increased JT cooling; and (4) water dry-out near the wellbore did not prevent hydrate formation. Finally, the role of capillary pressure was quite complex, as it reduced the formation of hydrates in certain cases and increased in other cases. Simulating this process with heat flow and hydrate reactions was also shown to present severe numerical issues. It was critical to select convergence criteria and linear system tolerances to avoid large material balance and numerical errors.