In this work, we present a kinetic simulation model for gas hydrates in porous media using the Operator-Based Linearization (OBL) technique. The OBL approach introduces algebraic operators that represent the physical terms in the mass and energy balance equations. Operators are calculated only in supporting points comprising the discretized parameter space, and operator values and partial derivatives for linear system assembly are readily obtained through (multi-)linear interpolation. Taking advantage of this setup, the implementation of advanced thermodynamic models for hydrate formation and dissociation under kinetic assumptions is simplified. We test the assumptions for thermodynamic modelling by analysing the Gibbs energy surfaces of the fluid and hydrate phases and demonstrate that, in the limit, the thermodynamic equilibrium for both kinetic and equilibrium reaction models is equivalent. We compare the simulation results with the published experimental results for CH₄-hydrates and extend the assessment to a CO₂-hydrate formation experiment in a semi-batch, constant-pressure configuration. The model reproduces the main pressure–temperature transients and hydrate evolution for both CH₄- and CO₂-systems. We demonstrate applicability at core scale for hydrate formation and, at field scale, for gas production from CH₄-hydrates by thermal stimulation and depressurization. The interaction of thermal-compositional phenomena (phase changes, adiabatic expansion, kinetic rates, and reaction enthalpy) gives rise to highly nonlinear physics that an appropriate OBL discretization resolves. Overall, the patterns of hydrate formation and dissociation are highly sensitive to the kinetic-rate inputs; hence, the appropriate choice of the reaction model remains a key consideration from both physical and numerical perspectives.

Hydrocarbon fuels are widely recognized as significant contributors to climate change and the rising levels of CO2 in the atmosphere. As a result, it is crucial to reduce the net carbon intensity of energy derived from these fuels. This study explores the feasibility of using dimethyl ether (DME), produced through the hydrogenation of CO2, as a low-carbon method for generating electricity from hydrocarbon fuels. The proposed approach involves capturing the emitted CO2 during combustion and utilizing it to produce the necessary DME in a closed cycle. It is shown that for a mature reservoir in the Middle East, this method can mitigate approximately 75% of the CO2 emissions released from burning the produced oil. By incorporating zero-carbon electricity throughout the process, the total abatement of CO2 can reach 85%. Furthermore, the study highlights the importance of improving the DME utilization factor (bbl-oil/tDME). By optimizing this factor, high abatement rates can be achieved. However, it is important to note that implementing this method comes with a high exergetic cost. During a certain period in the field’s lifetime, the invested energy exceeds the energy produced. The stages with the highest exergy consumption are CO2 capture and hydrogen production.

To achieve effective long-term CO2 storage in saline aquifers, it is essential to understand and monitor CO2 distribution and trapping mechanisms, which are significantly influenced by groundwater flow. This study investigates the impact of background flow velocity and direction on CO2 plume behavior and different trapping mechanisms (residual and solubility) using numerical analysis. The results of simulation show that in the flat (0° dip) model, increasing background flow velocity significantly extends the plume migration distance, enhancing both solubility and residual trapping through a larger CO2-water contact area and increased pore space occupation. The analysis is further extended to a dipping aquifer scenario to assess the role of groundwater flow direction. In the co-current flow case, where water and CO2 move in the same direction, the plume attains its maximum lateral extension, resulting in the highest storage efficiency. Conversely, in the counter-current flow scenario, where CO2 and water move in opposite directions, lower CO2 trapping is observed because, particularly at high velocity, the drag force exerted by water overcomes buoyancy force and limits further plume extension.

This study uses the concept of exergy-return on exergy-investment (ERoEI) to evaluate the life-cycle exergetic efficiency and CO₂ intensity (grams CO₂ per MJ of electricity) of (diabatic and adiabatic) compressed air energy storage (CAES) systems. Several CAES configurations are assessed under defined system boundaries, including diabatic systems powered by methane (CH₄) or hydrogen (H2), and adiabatic system with a thermal energy storage (TES) facility.

The results show that conventional (diabatic) CAES system powered by natural gas has the lower exergetic efficiency and higher CO2 intensity compared to adiabatic CAES due to the heat dissipation during compression stage and additional fuel requirements for reheating the air during expansion. Integrating carbon capture and storage (CCS) plant with conventional diabatic CAES can nearly halve the CO₂ intensity for electricity generation although the additional exergy investment for the CCS process reduces the exergetic efficiency of the system. Transitioning to green H2 (produced from low-carbon electricity) as the primary turbine fuel in the diabatic CAES results in a 65–76 % reduction in CO₂ intensity. However, the average exergetic efficiency of system decreases by around 10 %, mainly due to the substantial exergy investment associated with hydrogen production. It is also found that the adiabatic CAES system integrated with TES demonstrates the highest thermodynamic and environmental performance. When 100 % of compression heat is captured and reused during discharge phase, the system reaches ERoEI values up to 61 % with CO2 intensity of 12–26 g CO₂ per MJe.

Disclaimer: The results and performance metrics presented in this study are based on modelled scenarios and literature-derived parameters under defined system boundaries. Actual performance of CAES systems may vary depending on site-specific conditions, technology maturity, and operational configurations. All efficiency values, CO₂ intensity estimates, and comparative assessments should be interpreted within the context of the assumptions and limitations described herein. This study does not constitute a commercial endorsement or performance guarantee. The authors have made every effort to ensure accuracy but accept no liability for decisions made based on this analysis.

The gas displacement in porous media is a crucial process with extensive industrial and environmental applications. A notable example is underground hydrogen storage, where it is important to understand hydrogen mixing with cushion gas. The current paper explores anomalies in dispersion behaviour of gas mixtures under opposing flow directions (injection and production) from a modelling perspective. Due to the gaseous nature of the system, it presents significant complexities due to non-ideal mixing, compressibility, and higher diffusivity compared to Newtonian fluid transport. The findings reveal distinct dispersion behaviour during injection and production, where augmenting the mixture non-ideality enhanced the non-unique behaviour. In contrast to the dispersivity seen in Newtonian fluid flow in porous media, our research identifies that dispersivity in gas displacement depends not only on the porous medium but also on the gaseous components’ properties.

Halite precipitation during CO2 injection can significantly reduce injectivity and impact long-term storage in saline aquifers and depleted reservoirs. However, the impact of geological porous media on salt precipitation and brine movement is not fully undersood. This study explores salt precipitation dynamics from the pore to core scale using microfluidic experiments and core-flooding tests. Microfluidic results reveal three distinct phases of salt deposition: slow evaporation, rapid evaporation, and complete dry-out. Heterogeneous pore structures retain more initial water, leading to localized salt accumulation due to capillary effects. Core-scale experiments show that permeability strongly influences salt penetration and porosity reduction. In high-permeability cores, salt fronts extend deeper into the rock, causing up to 70% porosity reduction. In contrast, heterogeneous cores experience limited salt penetration but increased surface accumulation, suggesting that capillary pressure and brine redistribution control final deposition patterns. These findings highlight the complex interactions between fluid flow, rock properties, and salt crystallization, which provides valuable insights for predicting injectivity loss and optimizing CO2 storage strategies. Understanding these mechanisms is essential for improving reservoir management and ensuring the long-term stability of geological CO2 sequestration projects.

Fluid dispersion directly influences the transport, mixing, and efficiency of hydrogen storage in depleted gas reservoirs. Pore structure parameters, such as pore size, throat geometry, and connectivity, influence the complexity of flow pathways and the interplay between advective and diffusive transport mechanisms. Hence, these factors are critical for predicting and controlling flow behavior in the reservoirs. Despite its importance, the relationship between pore structure and dispersion remains poorly quantified, particularly under elevated flow conditions. To address this gap, this study employs pore network modeling (PNM) to investigate the influence of sandstone and carbonate structures on fluid flow properties at the micro-scale. Eleven rock samples, comprising seven sandstone and four carbonate, were analyzed. Pore network extraction from CT images was used to obtain detailed pore structure parameters and their statistical measures. Pore-scale simulations were conducted across 60 scenarios with varying average interstitial velocities and water as the injected fluid. Effluent hydrogen concentrations were measured to generate elution curves as a function of injected pore volumes (PV). This approach enables the assessment of the relationship between the dispersion coefficient and pore structure parameters across all rock samples at consistent average interstitial velocities. Additionally, dispersivity and n-exponent values were calculated and correlated with pore structure parameters.

The global energy transition requires novel carbon utilization methods to enable integrated and optimized low-carbon energy production. Coupling CO₂-based geothermal energy extraction with CO₂-enhanced oil recovery (EOR) represents a promising yet largely unexplored approach for improving resource efficiency and carbon sequestration. This study investigates the integration of CO₂-Plume Geothermal (CPG) energy production with CO₂-EOR in mature oil reservoirs using numerical simulations of conceptual heterogeneous reservoir models. The interplay between EOR and CPG performance in terms of energy production and CO₂ storage is evaluated and compared to understand the geotechnical implications of this integration. The analysis highlights that initiating CPG operations after EOR significantly benefits from the established CO₂ plume, facilitating immediate and efficient geothermal energy extraction. Results show that integrating CPG with EOR increases total energy recovery by 20%–50% relative to the energy produced by EOR alone, yielding CPG thermal power outputs ranging from 13 to 23 MW_th/km². Continued CO₂ injection during CPG operations further increases total CO₂ storage by 80%–280%, driven primarily by improved volumetric sweep of previously unswept reservoir volumes and enhanced CO₂ density resulting from reservoir cooling. While reservoir heterogeneity strongly influences oil recovery during EOR, its effect on CPG thermal output is less pronounced, since native reservoir fluids (oil and brine) have already been largely displaced during the EOR stage, and the CO₂ plume gradually stabilizes over time. These findings demonstrate the viability and advantages of integrated CO₂-EOR and CPG systems, offering insights into novel methods essential for sustainable subsurface resource management and climate-change mitigation.

Gas hydrates are crystalline compounds of water and small guest molecules, relevant both as a hazard in hydrocarbon production and CO2 sequestration, and as a potential energy resource in natural reservoirs. This work presents a kinetic simulation model for hydrate formation and dissociation in porous media, implemented using the Operator-Based Linearization (OBL) technique. We verify thermodynamic assumptions through Gibbs energy analysis, showing consistency between kinetic and equilibrium reaction models. The framework is validated against literature on methane hydrates and can be extended to CO2 systems. Applications are demonstrated at core and field scales, including gas production by depressurization and thermal stimulation. Results highlight the strong influence of kinetic parameters on hydrate behavior, underscoring the importance of selecting appropriate reaction models for accurate physical and numerical predictions.

Joule-Thomson cooling during CO2 injection into low-pressure fields can lead to injectivity impairment due to hydrate formation. This paper presents axial-symmetric flow model, which can be used to predict propagation of temperature and CO2 fronts during CO2 injection into porous formations accounting for Joule-Thomson cooling and unsteady-state delayed heat exchange between the reservoir and the adjacent formations. The solution of the 1D flow is validated by comparing with the quasi 2D analytical heat-conductivity solution. The non-steady state heat exchange results in a temperature front that propagates without limit into the reservoir with time. The temperature profiles exhibit a temperature decrease from the injected temperature to a minimum value, followed by a sharp increase to initial reservoir temperature on the temperature front. The solution allows plotting temperature-pressure (T-P) profiles at fixed moments in the CO2-water phase diagram. By changing injection parameters such as injection rate, the T-P trajectories allow for assessment of hydrate formation.

CO2 storage in deep saline aquifers is an effective strategy for reducing greenhouse gas emission. However, salt precipitation triggered by evaporation of water into injected dry CO2 causes injectivity reduction. Predicting the distribution of precipitated salts and their impact on near-well permeability remains challenging. Therefore, a detailed investigation of the interactions between salt precipitation and porous domain is essential for of revealing the mechanisms of pore blockage due to salt crystallization. Through series of microfluidic experiments, direct observations, coupled with detailed imaging processing, form the basis for explaining these phenomena and provide a relationship between water and salt saturations, highlighting the critical roles played by local capillary-driven flow and water film along grains in influencing water relocation. The results reveal two distinct types of salt crystallization: occurring inside the brine with smooth edges and at the CO2-brine interface with rough edges. Furthermore, the impact of local heterogeneity and surface wettability on salt precipitation patterns is discussed. The transition region between the porous domains and inlet/outlet channels exhibits brine backflow and a larger amount of salt accumulation. This paper presents a comprehensive analysis of the dynamic process of salt dry-out occurring during CO2 injection at the pore scale.

Depleted gas reservoirs are viable choices for large-scale CO₂ storage and to displace remaining methane volumes to further increase the storage capacity (EGR). However, deployment of such projects depends on an informed knowledge of the magnitude of mixing of the miscible gases, efficiency in displacing in-situ methane by CO₂, composition of the produced gas, and CO₂ storage capacity. This study focuses on the fundamental analysis of mixing during CO₂-EGR using a numerical approach. We propose to conduct very fine grid compositional simulations to provide insights into the mixing of CO₂ and methane in a gas reservoir at different reservoir and operational conditions. We first analyze a stratified layer model to understand the basic mechanisms of scale-dependency of dispersion and the significance of reservoir heterogeneity on fluid mixing. To consider more realistic reservoir heterogeneity, a two-dimensional stochastic reservoir model is analyzed to estimate dispersivity generated as fluids flow in porous media at different scales. Reservoir heterogeneity is represented by the Dykstra Parsons coefficient (V_DP) and autocorrelation length, and fluid properties are modeled depending on pressure and temperature conditions. Field-scale simulation is also performed to discuss the way dispersion is modeled in reservoir simulation affects simulated gas recovery. Our study shows that the variance of permeability and convective spreading are the primary causes of fluid mixing at any scale. In addition, molecular diffusion is not always negligible in gas mixing even in large-scale heterogeneous reservoirs since gas has much larger diffusivity than liquid. Furthermore, the mechanism of fluid mixing during CO₂-EGR is complex with the interplay between convective spreading, transverse dispersion (including molecular diffusion), and gravity segregation. Although geoscientists often assume numerical dispersion can represent physical dispersion, our study indicates this is an oversimplification and could cause significant errors in calculated gas recovery. Permeability heterogeneity is essential for the dispersion growth process and the final displacement behavior. Reservoir heterogeneity should be modeled with high-resolution grid models to analyze mixing behaviors more accurately.

Foam is utilized in enhanced oil recovery and CO₂ sequestration. Surfactant-alternating-gas (SAG) is a preferred approach for placing foam into reservoirs, due to it enhances gas injection and minimizes corrosion in facilities. Our previous studies with similar permeability cores show that during SAG injection, several banks occupy the area near the well where fluid exhibits distinct behaviour. However, underground reservoirs are heterogeneous, often layered. It is crucial to understand the effect of permeability on fluid behaviour and injectivity in a SAG process. In this work, coreflood experiments are conducted in cores with permeabilities ranging from 16 to 2300 mD. We observe the same sequence of banks in cores with different permeabilities. However, the speed at which banks propagate and their overall mobility can vary depending on permeability. At higher permeabilities, the gas-dissolution bank and the forced-imbibition bank progress more rapidly during liquid injection. The total mobilities of both banks decrease with permeability. By utilizing a bank-propagation model, we scale up our experimental findings and compare them to results obtained using the Peaceman equation. Our findings reveal that the liquid injectivity in a SAG foam process is misestimated by conventional simulators based on the Peaceman equation. The lower the formation permeability, the greater the error.

This paper discusses axi-symmetric flow during CO₂ injection into a non-adiabatic reservoir accounting for Joule-Thomson cooling and steady-state heat exchange between the reservoir and the adjacent layers by Newton's law. An exact solution for this 1D problem is derived and a new method for model validation by comparison with quasi 2D analytical heat-conductivity solution is developed. The temperature profile obtained by the analytical solution shows a temperature decrease to a minimum value, followed by a sharp increase to initial reservoir temperature on the temperature front. The temperature distribution head of the front is determined by the initial reservoir temperature, while the solution behind the front is determined by the temperature of injected CO₂. The analytical model exhibits stabilisation of the temperature profile and the cooled zone. The explicit formula for temperature distributions allows determining the maximum injection rate that avoids hydrate formation.

Injection of high-pressure CO2 into depleted gas reservoirs can lead to low temperatures promoting formation of hydrate in the near wellbore area resulting in reduced injection rates. The design of effective mitigation methods requires an understanding of the impact of crucial parameters on the formation and dissociation of CO2 hydrate within the porous medium under flowing conditions. This study investigates the influence of water saturation (ranging from 20 % to 40 %) on the saturation and kinetics of CO2 hydrate during continuous CO2 injection. The experiments were conducted under a medical X-ray computed tomography (CT) to monitor the dynamics of hydrate growth inside the core and to calculate the hydrate saturation profile. The experimental data reveal increase in CO2 hydrate saturation with increasing water saturation levels. The extent of permeability reduction is strongly dependent on the initial water saturation: beyond a certain water saturation the core is fully blocked. For water saturations representative of the depleted gas fields, although the amount of generated hydrate is not sufficient to fully block the CO2 flow path, a significant reduction in permeability (approximately 80 %) is measured. It is also observed that the volume of water + hydrate phases increases during hydrate formation, indicating a lower-than-water density for CO2 hydrate. Having a history of hydrate at the same water saturation leads to an increase in CO2 consumption compared to the primary formation of hydrate, confirming the existence of the water memory effect in porous media.

Hydrogen plays a crucial role in the transition to low-carbon energy systems, especially when integrated into energy storage applications. In this study, the concept of exergy-return on exergy-investment (ERoEI) is applied to investigate the exergetic efficiency (defined as the ratio of useful exergy output to invested exergy input) and CO2 equivalent intensity of the hydrogen supply chain, with a specific focus on the underground hydrogen storage process. Our findings reveal that the overall exergetic efficiency of the electricity-to-hydrogen-to-electricity conversion process can reach up to 25 %. Among the hydrogen production methods, green hydrogen, produced via electrolysis powered by renewable energy, exhibits the lowest CO2 equivalent intensity. Blue hydrogen, produced from natural gas with carbon capture, can significantly reduce the carbon footprint of electricity generation, though this advantage comes at the expense of decreased exergetic efficiency. Analysis further indicates that the exergetic efficiency of underground storage components ranges from 72 % to 92 %. A substantial fraction of the exergy is lost during compression and injection of the stored hydrogen. Nevertheless, the subsurface operations contribute a minimal CO2 emission, between 1.46–4.56 grams of equivalent CO2 per megajoule (gr-CO2eq/MJ) when powered by low-carbon energy sources. Furthermore, it is found that hydrogen loss in the reservoir, along with methane and hydrogen leak during surface operations, notably affects the overall efficiency of the storage process.

The success of CO2 storage projects largely depends on addressing formation damage, such as salt precipitation, hydrate formation, and fines migration. While analytical models for reservoir behaviour during CO2 storage in aquifers and depleted gas fields are widely available, models addressing formation damage and injectivity decline are scarce. This work aims to develop an analytical model for CO2 injection in a layer-cake reservoir, considering permeability damage. We extend Dietz’s model for gravity-dominant flows by incorporating an abrupt permeability decrease upon the gas-water interface arrival in each layer. The exact Buckley-Leverett solution of the averaged quasi-2D (x, z) problem provides explicit formulae for sweep efficiency, well impedance, and skin factor of the injection well. Our findings reveal that despite the induced permeability decline and subsequent well impedance increase, reservoir sweep efficiency improves, enhancing storage capacity by involving a larger rock volume in CO2 sequestration. The formation damage factor d, representing the ratio between damaged and initial permeabilities, varies from 0.016 in highly damaged rock to 1 in undamaged rock, resulting in a sweep efficiency enhancement from 1–3% to 50–53%. The developed analytical model was applied to predict CO2 injection into a depleted gas field.

Foam is a promising means to assist in the permanent, safe subsurface sequestration of CO₂, whether in aquifers or as part of an enhanced-oil-recovery (EOR) process. Here we review the advantages demonstrated for foam that would assist CO₂ sequestration, in particular sweep efficiency and residual trapping, and the challenges yet to be overcome. We also review the research and field-trial literature on CO₂ foam sweep efficiency, capillary gas trapping in foam, issues involved in surfactant selection for CO₂ foam applications, foam field trials, and the state of the art from laboratory and modelling research on CO₂ foam properties, in order to present the prospects and challenges for foam-assisted CO₂ sequestration. Challenges to foam-assisted CO₂ sequestration include the following: 1) verifying the advantages indicated by laboratory research at the field scale 2) optimizing surfactant performance, while further reducing cost and adsorption if possible 3) long-term chemical stability of surfactant, and dilution of surfactant in the foam bank by flow of water. Residual gas must reside in place for decades, even if surfactant degrades or is diluted. 4) optimizing injectivity and sweep efficiency in the field-design strategy.

The study aims to quantitatively assess the risk of hydrate formation within the porous formation and its consequences to injectivity during storage of CO₂ in depleted gas reservoirs considering low temperatures caused by the Joule Thomson (JT) effect and hydrate kinetics. 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 compositional thermal reservoir simulator is used to model the fluid and heat flow of CO₂ through a reservoir initially composed of brine and methane. The simulator can model the formation and dissociation of both methane and CO₂ 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 CO₂. Simulation results for a simplified model were verified with results from the literature. The key results of this work are: (1) The Joule-Thomson 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 CO₂ 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, where 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.

Injection of high-pressure CO₂ into depleted gas reservoirs can lead to low temperatures promoting formation of hydrate in the near wellbore area resulting in reduced injection rates. The design of effective mitigation methods requires an understanding of the impact of crucial parameters on the formation and dissociation of CO₂ hydrate within the porous medium under flowing conditions. This study investigates the influence of water saturation (ranging from 20% to 40%) on the saturation and kinetics of CO₂ hydrate during continuous CO₂ injection. The experiments were conducted under a medical X-ray computed tomography (CT) to monitor the dynamics of hydrate growth inside the core and to calculate the hydrate saturation profile. The experimental data reveal increase in CO₂ hydrate saturation with increasing water saturation levels. The extent of permeability reduction is strongly dependent on the initial water saturation: beyond a certain water saturation the core is fully blocked. For water saturations representative of the depleted gas fields, although the amount of generated hydrate is not sufficient to fully block the CO₂ flow path, a significant reduction in permeability (approximately 80%) is measured. It is also observed that the volume of water+hydrate phases increases during hydrate formation, indicating a lower-than-water density for CO₂ hydrate. Having a history of hydrate at the same water saturation leads to an increase in CO₂ consumption compared to the primary formation of hydrate, confirming the existence of the water memory effect in porous media.