Seis2Rock: A Data-Driven Approach to Direct Petrophysical Inversion of Pre-Stack Seismic Data

The inversion of petrophysical parameters from seismic data represents a fundamental step in the process of characterizing the subsurface. We propose a novel, data-driven approach named Seis2Rock that utilizes optimal basis functions learned from well log information to directly link band-limited petrophysical reflectivities to pre-stack seismic data. Seis2Rock is composed of two stages: training and inference. During training, a set of optimal basis functions are identified by performing singular value decomposition on one or more synthetic AVO gathers created from measured or rock-physics synthesized elastic well-logs. In inference, seismic pre-stack data are first projected into a set of band-limited petrophysical properties using the previously computed basis functions; this is followed by regularized post-stack seismic inversion of the individual properties. In this work, we apply the Seis2Rock methodology to a synthetic dataset based on the Smeaheia reservoir model and the open Volve field dataset. Numerical results reveal the ability of the proposed method in recovering accurate porosity, shale content, and water saturation models. Finally, the proposed methodology is applied in the context of reservoir monitoring to invert time-lapse, pre-stack seismic data for water saturation changes

Modeling Lost-Circulation into Fractured Formation in Rock Drilling Operations

Loss of circulation while drilling is a challenging problem that may interrupt drilling operations, reduce efficiency, and increases cost. When a drilled borehole intercepts conductive faults or fractures, lost circulation manifests as a partial or total escape of drilling, workover, or cementing fluids into the surrounding rock formations. Studying drilling fluid loss into a fractured system has been investigated using laboratory experiments, analytical modeling, and numerical simulations. Analytical modeling of fluid flow is a tool that can be quickly deployed to assess lost circulation and perform diagnostics, including leakage rate decline and fracture conductivity. In this chapter, various analytical methods developed to model the flow of non-Newtonian drilling fluid in a fractured medium are discussed. The solution methods are applicable for yield-power-law, including shear-thinning, shear-thickening, and Bingham plastic fluids. Numerical solutions of the Cauchy equation are used to verify the analytical solutions. Type-curves are also described using dimensionless groups. The solution methods are used to estimate the range of fracture conductivity and time-dependent fluid loss rate, and the ultimate total volume of lost fluid. The applicability of the proposed models is demonstrated for several field cases encountering lost circulations

Effect of methyl orange on wettability of sandstone formations: Implications for enhanced oil recovery

With the increasing global population, fossil fuel resources still represent a main contributor to the energy supply, despite the progress made in the field of renewable energies. Large quantities of residual oil from mature reservoirs cannot be produced through primary and secondary recovery methods. Among alternative recovery techniques, chemically enhanced oil recovery methods are attracting considerable interest to increase the hydrocarbon recovery from oil-bearing geological formations. The wettability of any particular formation can be used to predict the oil recovery factor of a reservoir based on its wetting state. However, due to the complex nature of geological porous media, special treatments are required to control the wetting characteristics for improving the oil recovery. In this work, methyl orange (MO), a hazardous pollutant widely discharged in industrial wastewater, was used as a chemical agent for the purpose of altering the wettability. Initially, quartz substrates were aged with 10 − 2 mol/L n-decane/stearic acid solution to mimic natural geological conditions; then, stearic acid-aged quartz substrates were treated in various concentrations of MO (10, 25, 50, 75, and 100 mg/L) for 7 days at 50 °C, followed by advancing and receding contact angle measurements at various physico-thermal geological conditions (temperature 25, 50 °C, pressure 10, 15, 20 MPa, and brine salinity 0 – 0.3 M). Our results demonstrate that increasing the temperature, pressure, and salinity of quartz aged with stearic acid has a negative effect on the wettability (resulting in a higher hydrophobicity). However, at any constant physio-thermal condition, MO significantly alters the wettability of the organic-aged quartz substrates from oil-wet to water-wet conditions, thus improving oil recovery. The concentration of MO plays a critical role, with increasing concentrations favouring the water-wet conditions. Quartz aged with MO at a concentration of 100 mg/L shows water-wet behaviour, with the lowest advancing and receding contact angles of 31 ° and 29 °, respectively, at 25 °C, 20 MPa, and 0.3 M salinity. The findings of this study provide new insights that can be useful for disposing MO in deep underground reservoirs rather than discharging into the hydrosphere, thus mitigating climate change. In addition, the present data can be helpful for improving the oil productivity from sandstone reservoirs

Evaluation of cubic, PC-SAFT, and GERG2008 equations of state for accurate calculations of thermophysical properties of hydrogen-blend mixtures

Hydrogen (H2) is a clean fuel and key enabler of energy transition into green renewable sources and a method of achieving net-zero emissions by 2050. Underground H2 storage (UHS) is a prominent method offering a permanent solution for a low-carbon economy to meet the global energy demand. However, UHS is a complex procedure where containment security, pore-scale scattering, and large-scale storage capacity can be influenced by H2 contamination due to mixing with cushion gases and reservoir fluids. The literature lacks comprehensive investigations of existing thermodynamic models in calculating the accurate transport properties of H2-blend mixtures essential to the efficient design of various H2 storage processes. This work benchmarks cubic equations of state (EoSs), namely Peng–Robinson (PR) and Soave Redlich–Kwong (SRK) and their modifications by Boston–Mathias (PR-BM) and Schwartzentruber–Renon (SR-RK), for their reliability in predicting the thermophysical properties of binary and ternary H2-blend mixtures, including CH4, C2H6, C3H8, H2S, H2O, CO2, CO, and N2, in addition to Helmholtz-energy-based EoSs (i.e., PC-SAFT and GERG2008). The benchmarked models are regressed against the experimental data for vapor–liquid equilibrium (VLE) that covers a wide range of pressures (0.01 to 101 MPa), temperatures (92 K to 367 K), and mole fractions (0.001 to 0.90) of H2. The novelty of this work is in benchmarking and optimizing the parameters of the mentioned EoSs to study VLE envelopes, densities, and other critical transport properties, such as heat capacity and the Joule–Thomson coefficient of H2 mixtures in a wide range of associated conditions. The results highlight the significant effect of the temperature-dependent binary interaction parameters on the calculations of thermophysical properties. The SR-RK EoS demonstrated the highest agreement with VLE data among the cubic EoSs with a low root mean square error and absolute average deviation. The PC-SAFT VLE models demonstrated results comparable to the SR-RK. The sensitivity analysis highlighted the high influence of impurity on changing the thermophysical behavior of H2-blend streams during the H2 storage process

Thermodynamic modeling of hydrogen–water systems with gas impurity at various conditions using cubic and PC-SAFT equations of state

Hydrogen (H2) has emerged as a viable solution for energy storage of renewable sources, supplying off-seasonal demand. Hydrogen contamination due to undesired mixing with other fluids during operations is a significant problem. Water contamination is a regular occurrence; therefore, an accurate prediction of H2-water thermodynamics is crucial for the design of efficient storage and water removal processes. In thermodynamic modeling, the Peng–Robinson (PR) and Soave Redlich–Kwong (SRK) equations of state (EoSs) are widely applied. However, both EoSs fail to predict the vapor-liquid equilibrium (VLE) accurately for H2-blend mixtures with or without fine-tuning binary interaction parameters due to the polarity of the components. This work investigates the accuracy of two advanced EoSs: the Schwartzentruber and Renon modified Redlich–Kwong cubic EoS (SR-RK) and perturbed-chain statistical associating fluid theory (SAFT) in predicting VLE and solubility properties of H2 and water. The SR-RK involves the introduction of polar parameters and a volume translation term. The proposed workflow is based on optimizing the binary interaction coefficients using regression against experimental data that cover a wide range of pressure (0.34 to 101.23 MPa), temperature (273.2 to 588.7 K), and H2 mole fraction (0.0004 to 0.9670) values. A flash liberation model is developed to calculate the H2 solubility and water vaporization at different temperature and pressure conditions. The model captures the influence of H2-gas (CO2) impurity on VLE. The results agreed well with the experimental data, demonstrating the model\u27s capability of predicting the VLE of hydrogen-water mixtures for a broad range of pressures and temperatures. Optimized coefficients of binary interaction parameters for both EoSs are provided. The sensitivity analysis indicates an increase in H2 solubility with temperature and pressure and a decrease in water vaporization. Moreover, the work demonstrates the capability of SR-RK in modeling the influence of gas impurity (i.e., H2–CO2 mixture) on the H2 solubility and water vaporization, indicating a significant influence over a wide range of H2–CO2 mixtures. Increasing the CO2 ratio from 20% to 80% exhibited almost the opposite behavior of H2 solubility compared to the pure hydrogen feed solubility. Finally, the work emphasizes the critical selection of proper EoSs for calculating thermodynamic properties and the solubility of gaseous H2 and water vaporization for the efficient design of H2 storage and fuel cells

Numerical Reliability and CPU Time for the Mixed Methods applied to Flow Problems in Porous Media

This work is devoted to the numerical reliability and time requirements of the Mixed Finite Element (MFE) and Mixed-Hybrid Finite Element (MHFE) methods. The behavior of these methods is investigated under the influence of two factors: the mesh discretization and the medium heterogeneity. We show that, unlike the MFE, the MHFE "suffers" with the presence of flatted triangular elements. A numerical reliability analyzing software (Aquarels) is used to detect the instability of the matrix-inversion code generated by MAPLE which is used in the MHFE code. We also show that the spectral condition number of the algebraic systems furnished by both methods in heterogeneous media grows up linearly according to the smoothness of the hydraulic conductivity. Furthermore, it is found that the MHFE could accumulate numerical errors if the conductivity varies abruptly in space. Finally, we compare running-times for both algorithms by giving various numerical experiments

Hydrogen, carbon dioxide, and methane adsorption potential on Jordanian organic-rich source rocks: Implications for underground H2 storage and retrieval

Hydrogen (H2) storage in geological formations offers a potential large-scale solution suitable for an industrial-scale hydrogen economy. However, the presence of organic residuals can significantly influence the H2 storage efficiency, as well as cushion gas performance, such as CO2 and CH4, injected to maintain healthy reservoir pressure. Thus, the H2 storage efficiency and cushion gas selectivity were thoroughly investigated in this work based on H2, CO2, and CH4 adsorption measurements using, for the first time, actual organic-rich carbonate-rich Jordanian source rock samples (TOC = 13 % to 18 %), measured at 60 °C temperature and a wide range of pressure (0.1 – 10.0 MPa). Initially, the samples were characterized using various analytical methods. Results demonstrated that H2 adsorption capacities reached up to 0.47 mol/kg at 9.0 MPa. The measured adsorption of CO2 was four times higher than H2. An increase in TOC significantly decreased H2 adsorption compared to CO2 and CH4. Additionally, CO2 demonstrated preferential behavior as a cushion gas compared to CH4, attributed mainly to the calcite content and presence of carboxyl and sulfonyl groups. This study provides fundamental data for understanding H2 potential storage issues in an organic-rich rock formation and thus aids in the industrial implementation of an H2 supply chain

Residual trapping of CO2, N2, and a CO2-N2 mixture in Indiana limestone using robust NMR coreflooding: Implications for CO2 geological storage

Carbon capture and sequestration (CCS) in geological formations is a prominent solution for reducing anthropogenic carbon emissions and mitigating climate change. The capillary trapping of CO2 is a primary trapping mechanism governed by the pressure difference between the wetting and nonwetting phases in a porous rock, making the latter a key input parameter for dynamic simulation models. During the CCS operational process, however, the CO2 is prone to contamination by impurities from various sources such as surfaces (e.g., pipelines and tanks) and the subsurface (e.g., existing natural gas). Such contamination can strongly influence the overall CO2 wettability, storage capacity, and containment security. Hence, the present study uses the nuclear magnetic resonance (NMR) core flooding technique to investigate and compare the residual saturations of pure CO2, pure N2, and a 50:50 CO2/N2 mixture in an Indiana limestone. The longitudinal and transverse relaxation times (T1 and T2) are measured to examine the displacement process of the pore network, and the trapping mechanism is evaluated at the pore scale as a determinant of the field-scale flow behavior. The NMR T1-T2 and 2D maps are used to observe the fluid configurations in the pore network, and the T1/T2 ratios are used to evaluate the microscopic wettability of the limestone grains by the pore-space fluids following each drainage/imbibition process step. The results indicate substantial residual gas trapping in the rock for the CO2-brine, N2-brine, and CO2/N2-brine systems, corresponding to gas saturations of 25%, 27%, and 26%, respectively. In the CO2-brine system, the intermolecular interplay between the CO2-enriched brine and limestone grains results in a higher T1/T2 ratio and significantly reduces the hydrophilicity of the limestone. Furthermore, the NMR T2 distribution reveals the occurrence of preferential water displacement into the large pores (r \u3e 1 m) and from the intermediate pores (0.03 m \u3c r \u3c 1 m), whereas water remains immobile in the smaller pores (r \u3c 0.03 m). The insignificant difference in residual trapping saturation between pure CO2 and the CO2-N2 mixture indicates the potential to allow for impurities in the CO2 phase in CCS without reducing the residual trapping capacity. Thus, the present work provides comprehensive information on the impact of gas injection on residual gas trapping in subsurface geological formations at the pore scale, thereby aiding in the development of CCS and other potential applications in enhanced oil recovery (EOR)

Influence of pressure, temperature and organic surface concentration on hydrogen wettability of caprock, implications for hydrogen geo-storage

Hydrogen (H2) as a cleaner fuel has been suggested as a viable method of achieving the de-carbonization objectives and meeting increasing global energy demand. However, successful implementation of a full-scale hydrogen economy requires large-scale hydrogen storage (as hydrogen is highly compressible). A potential solution to this challenge is injecting hydrogen into geologic formations from where it can be withdrawn again at later stages for utilization purposes. The geo-storage capacity of a porous formation is a function of its wetting characteristics, which strongly influence residual saturations, fluid flow, rate of injection, rate of withdrawal, and containment security. However, literature severely lacks information on hydrogen wettability in realistic geological and caprock formations, which contain organic matter (due to the prevailing reducing atmosphere). We, therefore, measured advancing (θa) and receding (θr) contact angles of mica substrates at various representative thermo-physical conditions (pressures 0.1-25 MPa, temperatures 308–343 K, and stearic acid concentrations of 10−9 - 10−2 mol/L). The mica exhibited an increasing tendency to become weakly water-wet at higher temperatures, lower pressures, and very low stearic acid concentration. However, it turned intermediate-wet at higher pressures, lower temperatures, and increasing stearic acid concentrations. The study suggests that the structural H2 trapping capacities in geological formations and sealing potentials of caprock highly depend on the specific thermo-physical condition. Thus, this novel data provides a significant advancement in literature and will aid in the implementation of hydrogen geo-storage at an industrial scale

Influence of organic molecules on wetting characteristics of mica/H2/brine systems: Implications for hydrogen structural trapping capacities

Hypothesis: Actualization of the hydrogen (H2) economy and decarbonization goals can be achieved with feasible large-scale H2 geo-storage. Geological formations are heterogeneous, and their wetting characteristics play a crucial role in the presence of H2, which controls the pore-scale distribution of the fluids and sealing capacities of caprocks. Organic acids are readily available in geo-storage formations in minute quantities, but they highly tend to increase the hydrophobicity of storage formations. However, there is a paucity of data on the effects of organic acid concentrations and types on the H2-wettability of caprock-representative minerals and their attendant structural trapping capacities. Experiment: Geological formations contain organic acids in minute concentrations, with the alkyl chain length ranging from C4 to C26. To fully understand the wetting characteristics of H2 in a natural geological picture, we aged mica mineral surfaces as a representative of the caprock in varying concentrations of organic molecules (with varying numbers of carbon atoms, lignoceric acid C24, lauric acid C12, and hexanoic acid C6) for 7 days. To comprehend the wettability of the mica/H2/brine system, we employed a contact-angle procedure similar to that in natural geo-storage environments (25, 15, and 0.1 MPa and 323 K). Findings: At the highest investigated pressure (25 MPa) and the highest concentration of lignoceric acid (10−2 mol/L), the mica surface became completely H2 wet with advancing (θa= 106.2°) and receding (θr=97.3°) contact angles. The order of increasing θa and θr with increasing organic acid contaminations is as follows: lignoceric acid \u3e lauric acid \u3e hexanoic acid. The results suggest that H2 gas leakage through the caprock is possible in the presence of organic acids at higher physio-thermal conditions. The influence of organic contamination inherent at realistic geo-storage conditions should be considered to avoid the overprediction of structural trapping capacities and H2 containment security