Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “3D Geologic Model”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 217 records · Page 12

Real Time 3D Observations of Portland Cement Carbonation at CO2 Storage Conditions

Depleted oil reservoirs are considered a viable solution to the global challenge of CO 2 storage. A key concern is whether the wells can be suitably sealed with cement to hinder the escape of CO 2 . Under reservoir conditions, CO 2 is in its supercritical state, and the high pressures and temperatures involved make real-time microscopic observations of cement degradation experimentally challenging. Here, we present an in situ 3D dynamic X-ray micro computed tomography (μ-CT) study of well cement carbonation at realistic reservoir stress, pore-pressure, and temperature conditions. The high-resolution time-lapse 3D images allow monitoring the progress of reaction fronts in Portland cement, including density changes, sample deformation, and mineral precipitation and dissolution. By switching between flow and nonflow conditions of CO 2 -saturated water through cement, we were able to delineate regimes dominated by calcium carbonate precipitation and dissolution. For the first time, we demonstrate experimentally the impact of the flow history on CO 2 leakage risk for cement plugging. In-situ μ-CT experiments combined with geochemical modeling provide unique insight into the interactions between CO 2 and cement, potentially helping in assessing the risks of CO 2 storage in geological reservoirs.

58 GEOSCIENCES↗

Bayesian seismic inversion based on rock-physics prior modeling for the joint estimation of acoustic impedance, porosity and lithofacies

We propose a Bayesian approach for seismic inversion to estimate acoustic impedance, porosity and lithofacies within the reservoir conditioned to post-stack seismic and well data. The link between elastic and petrophysical properties is given by a joint prior distribution for the logarithm of impedance and porosity, based on a rock-physics model. The well conditioning is performed through a background model obtained by well log interpolation. Two different approaches are presented: in the first approach, the prior is defined by a single Gaussian distribution, whereas in the second approach it is defined by a Gaussian mixture to represent the well data multimodal distribution and link the Gaussian components to different geological lithofacies. The forward model is based on a linearized convolutional model. For the single Gaussian case, we obtain an analytical expression for the posterior distribution, resulting in a fast algorithm to compute the solution of the inverse problem, i.e. the posterior distribution of acoustic impedance and porosity as well as the facies probability given the observed data. For the Gaussian mixture prior, it is not possible to obtain the distributions analytically, hence we propose a Gibbs algorithm to perform the posterior sampling and obtain several reservoir model realizations, allowing an uncertainty analysis of the estimated properties and lithofacies. Both methodologies are applied to a real seismic dataset with three wells to obtain 3D models of acoustic impedance, porosity and lithofacies. The methodologies are validated through a blind well test and compared to a standard Bayesian inversion approach. Using the probability of the reservoir lithofacies, we also compute a 3D isosurface probability model of the main oil reservoir in the studied field.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Reconstituted Three-Dimensional Interactive Imaging

A method combines two-dimensional images, enhancing the images as well as rendering a 3D, enhanced, interactive computer image or visual model. Any advanced compiler can be used in conjunction with any graphics library package for this method, which is intended to take digitized images and virtually stack them so that they can be interactively viewed as a set of slices. This innovation can take multiple image sources (film or digital) and create a "transparent" image with higher densities in the image being less transparent. The images are then stacked such that an apparent 3D object is created in virtual space for interactive review of the set of images. This innovation can be used with any application where 3D images are taken as slices of a larger object. These could include machines, materials for inspection, geological objects, or human scanning. Illuminous values were stacked into planes with different transparency levels of tissues. These transparency levels can use multiple energy levels, such as density of CT scans or radioactive density. A desktop computer with enough video memory to produce the image is capable of this work. The memory changes with the size and resolution of the desired images to be stacked and viewed.

Hamilton, Joseph↗

Quality Analysis of Baseline Time-Lapse CSEM Data CarbonSAFE Project, North Dakota

Conference presentation at International Meeting for Applied Geoscience & Energy (IMAGE), Houston, TX, August 28 – September 2, 2022. When a series of time-lapse CSEM surveys are designed to measure the often-small variations in signal seen in CCUS projects, the primary concern must be collecting the highest-quality data and reducing as much noise as possible. This includes optimization during planning and feasibility, careful and consistent quality checks in the field, and transparent and repeatable postprocessing steps. We apply this practice to a baseline (prior to CO 2 injection) CSEM data set collected for a time-lapse survey in Center, North Dakota, as a part of the North Dakota CarbonSAFE project, and describe the rigorous quality control and assessment methods used, including noise removal and data validation with 1D and 3D inversion to tie results to borehole logs. The final result is an accurate and representative CSEM data set and information that stakeholders can use to inform future time-lapse survey costs and designs.

20 FOSSIL-FUELED POWER PLANTS↗

Plasma IMS Composition Measurements for Europa and the Other Galilean Moons

NASA and ESA are planning the joint Europa Jupiter System Mission (EJSM) to the Jupiter system with specific emphasis to Europa and Ganymede, respectively. The Japanese Space Agency is also planning an orbiter mission to explore Jupiter's magnetosphere and the Galilean satellites. For NASA's Jupiter Europa Orbiter (JEO) we are developing the 3D Ion Mass Spectrometer (IMS) with two main goals which can also be applied to the other Galilean moons, 1) measure the plasma interaction between Europa and Jupiter's magnetosphere and 2) infer the 4 pi surface composition to trace elemental and significant isotopic levels. The first goal supports the magnetometer (MAG) measurements, primarily directed at detection of Europa's sub-surface ocean, while the second gives information about transfer of material between the Galilean moons, and between the moon surfaces and subsurface layers putatively including oceans. The measurement of the interactions for all the Galilean moons can be used to trace the in situ ion measurements of pickup ions back to either Europa's or Ganymede's surface from the respectively orbiting spacecraft. The IMS instrument, being developed under NASA's Astrobiology Instrument Development Program, would maximally achieve plasma measurement requirements for JEO and EJSM while moving forward our knowledge of Jupiter system composition and source processes to far higher levels than previously envisaged. The composition of the global surfaces of Europa and Ganymede can be inferred from the measurement of ejected neutrals and pick-up ions using at minimum an in situ payload including MAG and IMS also fully capable of meeting Level 1 mission requirements for ocean detection and survey. Elemental and isotopic analysis of potentially extruded oceanic materials at the moon surfaces would further support the ocean objectives. These measurements should be made from a polar orbiting spacecraft about Europa or Ganymede at height 100 km. The ejecta produced by sputtering of the surfaces of Europa and Ganymede has been shown to be representative of the surface composition. Level 2 science on surface geology and composition can then be further enhanced by addition of the following: 3D Ion Neutral Mass Spectrometer (INNS), 3D plasma electron spectrometer (ELS), and hot plasma energetic particle instrument. The measurement approach is to alternate between times measuring pickup ions and times measuring plasma and magnetic field parameters along the spacecraft trajectory. By measuring the pickup ion energy, arrival direction and mass-per-charge, the ion can be traced back along the ejection trajectory to the approximate area of origin if the 3-D electric field and magnetic field are known. In situ observations of plasma flow velocities and vector magnetic fields can be used to determine the local convective electric field (E = -VXB) along the spacecraft trajectory. By combining this information with models of the magnetospheric interaction with Europa, one can generate 3D maps of the electric and magnetic field and compute the trajectories of the pickup ions back to the surface or exospheric points of origin. In the case of Ganymede there is the additional complexity of its own internal dipole magnetic field, while Io's volcanic activity introduces the complexity of a highly structured denser atmosphere. Callisto with its less globally extended exosphere will have a simpler interaction than for Europa (i.e., more like our moon). We will discuss these differences in light of the above proposed technique. Finally, the INNS observations and neutral exosphere models are needed to estimate production rates of pickup ions. The hot plasma measurements are needed to correct for sputtering rates which can be time dependent and electron plasma observations for electron impact ionization rates. Instrument characteristics, field-of-view requirements, modes of operation and effects of radiation on instrument functionality will be discussed.

Sittler, Edward↗

Newmark-Type Pseudo-Three-Dimensional Back-Analysis of Co-Seismic Landslides in Egkremnoi, Lefkada, Greece

A pseudo-three-dimensional (pseudo-3D) methodology to back-analyze coseismic landslides was developed and applied to 68mapped landslides, which occurred over approximately 1 km2area in Egkremnoi, Greece, during the 2015 Mw6.5 Lefkada earthquake. The methodology is based on a one-dimensional (1D) Newmark-type sliding block model to assess instability and a spatial projection in 3D topography to derive landslide geometry. The strength parameters for modeled landslides that best match the landslide location, area, and volume were derived through an iterative scheme that optimizes the match using predefined matching criteria. The range of different-sized landslides produced shear strength estimates from ~10 to 300 kPa and led to the derivation of a regionally averaged strength envelope characterized by a cohesion of 6 kPa and a friction angle of 53°for the highly fractured limestones that are encountered in this area. Compared to previous full 3Dslope stability analyses in this area, the friction angle using this methodology was found to be generally consistent, but the cohesion was lower. The presented methodology can provide a computationally efficient method to estimate the average shear strength of a geologic unit over large areas, especially where extensive field and laboratory tests on the materials are unavailable or difficult to conduct.

Weibing Gong↗

Hawai‘i Play Fairway (Final Report)

Most of Hawai'i's geothermal resources are blind—their manifestations, such as hot springs and steam vents, do not appear on the ground surface because the heated water flows far below. With the exception of $K\bar{i}lauea East Rift Zone$, in most areas of Hawai'i, high lateral permeability in the first kilometer below ground surface prevents surface thermal features from developing. As a methodology for discovering these blind resources, Play Fairway Analysis (PFA) involves finding potential locations of blind hydrothermal systems and describing potential geothermal sources in rift-zone settings. Using the PFA to find Hawai'i's geothermal resources, the University of Hawai'i (UH) conducted the Hawai‘i Play Fairway Project, Hawai'i's first statewide geothermal resource assessment since 1985. Sponsored by the U.S. Department of Energy, the Hawai'i Play Fairway Project provided an updated resource assessment, a roadmap for additional exploration activities, and the identification of areas for further exploration. Benefitting from UH's core competency in earth sciences and experienced geothermal researchers, the project comprised three phases. During the first phase, the team identified, compiled, and ranked existing geologic, groundwater, and geophysical datasets relevant to subsurface heat, fluid and permeability. Using a Bayesian statistical approach, the team developed a statistical methodology to integrate these data into a resource probability map. The team evaluated the confidence in the probability value and considered development viability of areas with geothermal resources. With these analyses, the team identified 10 locations in the Hawaiian Islands for exploration activities. For the second phase, the team collected new groundwater data in 10 locations across the state and new geophysical data on $L\bar{a}na‘i, Maui$, and central Hawai'i Island and modeled topographically induced stress to better characterize subsurface permeability. Analyzing the subsurface stresses, the team evaluated the potential for fracture-induced permeability. The team inverted the MT and gravity data to produce 3D models of resistivity and density, respectively, on $L\bar{a}na‘i$, across $Haleakal\bar{a}'s$ SW rift (Maui), and surrounding Mauna Kea (Hawai‘i Island). The team developed and applied a new method for incorporating depth information about resistivity, density, and potential for fracture-induced permeability into the statistical method for computing resource probability in these three focus areas. The team incorporated the new groundwater results with the new geophysical results and the calculations of potential for fracture-induced permeability to produce updated maps of resource probability and confidence. Through combining data from the first and second phases, the team determined locations for further exploration during the third phase. For MT and gravity surveys, the team recommended $Kaua'i's$ $L\bar{i}hu'e$ $Basin$, the east rift of $Maui's$ $Haleakal\bar{a}$ volcano, and the southwest rift of Hawai'i Island's Mauna Loa volcano. The MT and gravity surveys aimed to enable improved confidence in the resource potential in these locations. For drilling deep groundwater well(s), the team recommended Southeast Mauna Kea and $L\bar{a}na's$ $P\bar{a}l\bar{a}wai$ $Basin$. During the third phase, further exploration involved drilling a groundwater well in $L\bar{a}na's$ $P\bar{a}l\bar{a}wai$ $Basin$ and performing more geophysical surveys. We deepened an existing water well proximal to our target area on $L\bar{a}na'i$ due to funding constraints that precluded us from spudding a new well that would exceed 1km depth. Drilling was preceded by a number of substantial elements including: writing an Environmental Assessment and the subsequent legal process, performance of deviation logging, lowering a camera down the well, coordinating site preparation with $P\bar{u}lama$ $L\bar{a}na'i$, shipping the UH-owned rig interisland, procuring supplies, and leading 3 community meetings on $L\bar{a}na'i$. Drilling occurred 24/7 the entire month of June 2019 over which time $L\bar{a}na'i$ $Well$ 10 was deepened from 427 m to 1057 m, with continuous core collected. We measured a roughly linear temperature gradient averaging 42°C/km and a maximum bottom hole temperature of 66°C. This gradient is more than twice the background for Hawai'i and within a range of gradients measured in this depth range for some exploration wells within KERZ. We consider these results encouraging for $L\bar{a}na'i's$ resource potential and recommend following with a slim hole within $L\bar{a}na'i's$ caldera (our target zone) to ~ 2 km. Further, the positive implications such results have for the island of O‘ahu are substantial - the shield stage of O'ahu's volcanoes ended 1-2 My earlier. However, O'ahu uses more electricity than the rest of the islands combined, and the utility recently called for 500-700MW of firm, dispatchable renewable electricity on O'ahu by 2033. In Phase 3, we also collected limited new encouraging groundwater data, and updated our thoughts on the probabilities of fluid and permeability at resource depths (Pr F = 1; Pr P = mostly unconstrained). Ultimately, we advocate for using our final probability of heat, and confidence in this probability, to drive the next phase of exploration. We contend further development of geothermal in Hawai‘i will enable the state to achieve its 100% renewable policy objective and Hawai'i to transition off of fossil fuels through geothermal discovery and development. The project not only produced a large amount of data and expanded the existing knowledge of Hawai'i's geothermal resources, but also produced publications, theses, presentations, core photos, datasets, media reports, television interviews, community events, and a blog. Students and new professionals benefitted from the project's hands-on research experiences and educational opportunities and earned awards and recognition.

15 GEOTHERMAL ENERGY↗

The effects of well damage and completion designs on geoelectrical responses in mature wellbore environments

Well integrity is one of the major concerns in long-term geologic storage sites due to the potential risk of well leakage and groundwater contamination. Evaluating changes in electrical responses due to energized steel-cased wells has the potential to quantify and predict possible wellbore failures because any kind of breakage or corrosion along highly conductive well casings will have an impact on the distribution of the subsurface electrical potential. However, realistic wellbore-geoelectrical models that can fully capture fine-scale details of well completion design and the state of well damage at the field scale require extensive computational effort, or they can even be intractable to simulate. To overcome this computational burden while still keeping the model realistic, we have used the hierarchical finite-element method that represents electrical conductivity at each dimensional component (1D edges, 2D planes, and 3D cells) of a tetrahedral mesh. This allows well completion designs with real-life geometric scales and well systems with realistic, detailed, progressive corrosion and damage in our models. We have developed a comparison of possible discretization approaches of a multicasing completion design in the finite-element model. The effects of the surface casing and the coupling between concentric well casings as well as the effects of the degree and the location of well damage on the electrical responses are also examined. As a result, we analyze real surface electric field data to detect wellbore integrity failure associated with damage.

58 GEOSCIENCES↗

Final Project Report: Imaging Fault Zones Using a Novel Elastic Reverse-Time Migration Imaging Technique

Imaging fault zones and fractures is crucial for geothermal operators, providing important information for reservoir evaluation and management strategies. However, there are no existing techniques available for directly and clearly imaging fault zones, particularly for steeply dipping faults and fracture zones. In this project, we developed novel acoustic- and elastic-waveform inversion methods for high-resolution velocity model building. In addition, we developed acoustic and elastic reverse-time migration methods for high-resolution subsurface imaging of complex subsurface structures and steeply-dipping fault/fracture zones. We first evaluated and verified the improved capabilities of our newly developed seismic inversion and migration imaging methods using synthetic seismic data. Our numerical tests verified that our new methods directly image subsurface fracture/fault zones using surface seismic reflection data. We then applied our novel seismic inversion and migration imaging methods to a field 3D surface seismic dataset acquired at the Soda Lake geothermal field using Vibroseis sources. Our migration images of the Soda Lake geothermal field obtained using our seismic inversion and migration imaging algorithms revealed several possible fault/fracture zones. AltaRock Energy, Inc. is working with Cyrq Energy, Inc. to refine the geologic interpretation at the Soda Lake geothermal field. Trenton Cladouhos, Senior Vice President R&D of AltaRock, was very interested in our imaging results of 3D surface seismic data from the Soda Lake geothermal field. He planed to perform detailed interpretation of our images in collaboration with James Faulds and Holly McLachlan of University of Nevada at Reno. Using our high-resolution seismic inversion and migration imaging results can help determine the optimal locations to drill wells for geothermal energy production and reduce the risk of geothermal exploration.

47 OTHER INSTRUMENTATION↗

Geomechanics-Based Stochastic Analysis of Injection- Induced Seismicity

The production of geothermal energy from dry and low permeability reservoirs is achieved by water circulation in natural and/or man-made fractures, and is referred to as enhanced or engineered geothermal systems (EGS). Often, the permeable zones have to be created by stimulation, a process which involves fracture initiation and/or activation of discontinuities such as faults and joints due to pore pressure and the in-situ stress perturbations. The stimulation of a rock mass is often accompanied by multiple microseismic events. Micro-seismic events associated with rock failure in shear, and shear slip on new or pre-existing fracture planes and possibly their propagations. The microseismic signals contain information about the sources of energy that can be used for understanding the hydraulic fracturing process and the created reservoir properties. Detection and interpretation of microseismic events is useful for estimating the stimulated zone, created reservoir permeability and fracture growth, and geometry of the geological structures and the in-situ stress state. The process commonly is referred to as seismicity-based reservoir characterization (SBRC). Although, progress has been made by scientific & geothermal communities for quantitative and qualitative analysis of reservoir stimulation using SBRC several key questions remain unresolved in the analysis of micro-seismicity namely, variation of seismic activity with injection rate, delayed micro-seismicity, and the relation of stimulated zone to the injected volume and its rate, and the resulting reservoir permeability. In addition, the current approach to SBRC does not consider the full range of relevant poroelastic and thermoelastic phenomena and neglects the uncertainty in rock properties and in-situ stress in the data inversion process. The objective of this research and technology developments was to develop a 3D SBRC model that addresses these shortcomings by taking into account hydro-thermo-poro-mechanical mechanisms associated with injection and utilizing a state-of-the-art stochastic inversion procedure. The approach proposed herein is innovative and significantly improves the existing SBCR technology (e.g., Shapiro et al. 2003) for geothermal reservoirs in several ways. First, the current scope of the SBRC is limited with respect to the physical processes considered and the rock properties used. Usually, the geomechanics analyses within SBRC is limited to the pore pressure diffusion in the rock mass, which is modeled using a time-dependent parabolic equation and solved using a finite element algorithm with either a line or a point source. However, water injection induces both poroelastic and thermoelastic stresses in the rock mass which affect the stress state. In fact, it has been suggested that thermoelastic stresses can play a dominant role in reservoir seismicity (Ghassemi et al., 2007). We include these important effects by using a fully-coupled poro-thermoelastic constitutive equations for the rock mass which will be solved using a 3D finite element model with more realistic injection geometries such as multiple injection/extraction sources (and in fractures), uncertainty in the material parameters and the in-situ stress distribution to better reflect the pore pressure and stress distributions. In addition, we developed a 3D stochastic fracture network model to study MEQ generation in fracture rocks. The model was verified using laboratory experiments, and calibrated and applied to Newberry EGS stimulation. In previous SBRC approaches, the triggering of micro-seismicity is modeled base on the assumption that the prior stochastic criticality model of the rock mass is a valid and adequate description. However, this assumption often does not hold in the field. Thus, we improved upon the current SBRC approach by using the micro-seismic responses to estimate the hydraulic diffusivity as well as the criticality distribution itself within the field. In this way, instead of relying on our a priori knowledge of criticality distribution, we combine an initial probabilistic description of criticality with the information contained in microseismic measurements to arrive at criticality solutions that are conditioned on both field data and our prior knowledge. Previous SBRC have relied upon a deterministic inversion approach to estimate the permeability, and the extent of the stimulated zone, whereas a stochastic inversion algorithm that recognizes and quantifies the uncertainties in the prior model, the time evolution of pore pressure distributions (modeling errors), and the observed seismic events is developed and used herein to realistically assess the quality of the solution. Finally, we developed a technique for processing discrete MEQ data to estimate fracture network properties such as dip and dip directions. The approach was successfully applied to the Fenton Hill HRD experiment and the Newberry EGS with results in good agreement with field observations.

15 GEOTHERMAL ENERGY↗

DEEPEN 3D PFA Index Models for Exploration Datasets at Newberry Volcano

DEEPEN stands for DE-risking Exploration of geothermal Plays in magmatic ENvironments. As part of the development of the DEEPEN 3D play fairway analysis (PFA) methodology for magmatic plays (conventional hydrothermal, superhot EGS, and supercritical), index models needed to be developed to map values in geoscientific exploration datasets to favorability index values. This GDR submission includes those index models. Index models were created by binning values in exploration datasets into chunks based on their favorability, and then applying a number between 0 and 5 to each chunk, where 0 represents very unfavorable data values and 5 represents very favorable data values. To account for differences in how exploration methods are used to detect each play component, separate index models are produced for each exploration method for each component of each play type. Index models were created using histograms of the distributions of each exploration dataset in combination with literature and input from experts about what combinations of geophysical, geological, and geochemical signatures are considered favorable at Newberry. This is in attempt to create similar sized bins based on the current understanding of how different anomalies map to favorable areas for the different types of geothermal plays (i.e., conventional hydrothermal, superhot EGS, and supercritical). For example, an area of partial melt would likely appear as an area of low density, high conductivity, low vp, and high vp/vs. This means that these target anomalies would be given high (4 or 5) index values for the purpose of imaging the heat source. To account for differences in how exploration methods are used to detect each play component, separate index models are produced for each exploration method for each component of each play type. Index models were produced for the following datasets: - Geologic model - Alteration model - vp/vs - vp - vs - Temperature model - Seismicity (density*magnitude) - Density - Resistivity - Fault distance - Earthquake cutoff depth model

15 GEOTHERMAL ENERGY↗

Pisgah Lava Cave Communication Test: Science Case Study for the Networked Constellations Initiative

As part of the science case study for the Networked Constellations initiative, a team of JPL scientists explore the possibility of a mission to study the lava caves on Mars. Natural caves on Mars and the Moon present a unique opportunity to learn about the planetary geology and to provide a shelter for human explorers. Due to power and communication challenges, a network of assets has significant advantages over a single asset sent inside a cave. However, communication between the assets and the data downlink present significant difficulties due to the presence of rough walls, boulders, and other obstacles with unknown dielectric constant inside a typical cave, disturbing the propagation of the radio waves. A detailed study is needed to establish the limitations of the current communication technologies and to develop requirements for the new communication technology applicable to the cave environment. On May 4 of 2017, Konstantin Belov, Doug Ellison, and Abby Fraeman visited a lava cave in Pisgah, CA. The purpose of the visit was to build a 3D map of the cave, which could be used to create a model of radio wave propagation, and to conduct a series of communication tests using off-the-shelf equipment to verify the in-cave communication challenges. This experiment should be considered as a simple 'proof of concept' and is the subject of this report.

Belov, K.↗

Preliminary result of teleseismic double-difference relocation of earthquakes in the Molucca collision zone with a 3D velocity model

We have relocated hypocenters of earthquakes occurring in the Molucca collision zone and surrounding region taken from the BMKG catalog using teleseismic double-difference relocation algorithm (teletomoDD). We used P-wave arrival times of local, regional, and teleseismic events recorded at 304 recording stations. Over 7,000 earthquakes were recorded by the BMKG seismographicnetworkin the study region from April, 2009 toJune, 2014. We used a 3D regional-global nested velocity modelresulting fromprevious global tomographystudy. In this study, the3D seismic velocity model was appliedto theIndonesian region, whilethe1D seismicvelocity model (ak135)wasused for regions outside of Indonesia. Our relocation results show a better improvement in travel-time RMS residuals comparedto those of the BMKG catalog.Ourresultsalso show that relocation shifts were dominated intheeast-west direction, whichmaybeinfluenced by theexistingvelocity anomaly related to the reversed V-shaped slabbeneaththestudy region. Our eventrelocation results refine the geometry of slabs beneath the Halmahera and Sangihe arcs.

58 GEOSCIENCES↗

Study of seismic diffractions caused by a fracture zone at In Salah carbon dioxide storage project

The In Salah CO 2 storage project in Algeria has injected over 3 million tonnes of carbon dioxide into a water filled tight sand formation. Interferometric Synthetic Aperture Radar (InSAR) range change data revealed a double lobbed pattern of deformation that has been modeled as the opening of a sub-vertical fracture, or damage, zone. The location and geometry of the linear feature were subsequently verified using a seismic reflection survey. The elastic properties of the fracture zone, including anisotropic Poisson ratio ( ν ), Young's ( E ) and shear ( G ) moduli, were estimated from coupled geomechanical and hydrological modeling of surface deformation and pressure variations in the injection well. The elastic moduli reflect the fracture properties after CO 2 flow through the fracture zone. Thus, the seismic signature of the fracture zone could be used for monitoring the CO 2 plume. Using the estimated fracture model, we built two and three dimensional models consisting of an anisotropic fracture zone embedded within an isotropic background. Finite-difference modeling of seismic shot gathers allows us to estimate the effects of scattering from the fracture zone, potentially further constraining the geomechanical model. From the seismic modeling results, we find diffracted waves, induced by the fracture zone, which behave similar to point source diffractions. This modeling is intended to guide a search for diffraction events in the 3D surface seismic field data. The modeling results indicate that using the moduli estimated from geomechanical modeling, fracture scattered events would be 100 times lower amplitude than the interface reflections, and thus would be hard to detect. While diffracted waves are observed in the field data, which may imply the need for revision of the fracture model, including shape and elastic moduli, we are not able to match the field observation with our modeled events. This article presents a frontier study on the integration of geomechanical and geophysical methods at the In Salah site as a means to test the estimate of the subsurface CO 2 flooded fracture properties.

58 GEOSCIENCES↗

Machine learning-based inversion for acoustic impedance with large synthetic training data: Workflow and data characterization

Where wells are sparse or training data are difficult to label with high-quality wireline-derived impedance logs, machine learning (ML)-based inversion of acoustic impedance typically depends on small training data sets, leading to biased prediction. We have advanced a novel workflow that applies large synthetic seismic training data to reduce facies-related bias. Using a geologically realistic model as the truth model, we randomly select sparse seed wells to perform sequential Gaussian simulation (SGS) for impedance models of the same geometry and simulate facies variability. We implement random forest regression on 30 features extracted from the synthetic volume. We observe that more seed wells tend to reduce facies-induced bias by sampling more types of facies, resulting in a better prediction. We then focus on the responses of SGS models to facies changes, the number of seed wells necessary for a useful synthetic model, and how much a synthetic model can help ML-based inversion. Here, we observe that the SGS synthetic training model outperforms well-direct training in general. For modeled clastic shore-zone systems in Miocene Gulf of Mexico, two or more seed wells are necessary for a significant reduction of root-mean-square error and outliners, and improvement of facies imaging. In a field-data test, we apply a similar workflow to quantitatively predict acoustic impedance, which is then converted to a sand-volume map at a high-frequency sequence (10–100 m), revealing detailed facies and sandstone patterns. Such results are valuable in many geologic and engineering applications, such as hydrocarbon and CO 2 reservoir prospecting, reserve estimation, simulation, etc.

3D seismic↗

Monitoring Offshore CO 2 Sequestration Using Marine CSEM Methods; Constraints Inferred from Field- and Laboratory-Based Gas Hydrate Studies

Offshore geological sequestration of CO 2 offers a viable approach for reducing greenhouse gas emissions into the atmosphere. Strategies include injection of CO 2 into the deep-ocean or ocean-floor sediments, whereby depending on pressure–temperature conditions, CO 2 can be trapped physically, gravitationally, or converted to CO 2 hydrate. Energy-driven research continues to also advance CO 2 -for-CH 4 replacement strategies in the gas hydrate stability zone (GHSZ), producing methane for natural gas needs while sequestering CO 2 . In all cases, safe storage of CO 2 requires reliable monitoring of the targeted CO 2 injection sites and the integrity of the repository over time, including possible leakage. Electromagnetic technologies used for oil and gas exploration, sensitive to electrical conductivity, have long been considered an optimal monitoring method, as CO 2 , similar to hydrocarbons, typically exhibits lower conductivity than the surrounding medium. We apply 3D controlled-source electromagnetic (CSEM) forward modeling code to simulate an evolving CO 2 reservoir in deep-ocean sediments, demonstrating sufficient sensitivity and resolution of CSEM data to detect reservoir changes even before sophisticated inversion of data. Laboratory measurements place further constraints on evaluating certain systems within the GHSZ; notably, CO 2 hydrate is measurably weaker than methane hydrate, and >1 order of magnitude more conductive, properties that may affect site selection, stability, and modeling considerations.

58 GEOSCIENCES↗

DECOVALEX-2023: Task F2 Salt Final Report

The subject of Task F of DECOVALEX-2023 concerns performance assessment modelling of radioactive waste disposal in deep mined repositories. The primary objectives of Task F are to build confidence in the models, methods, and software used for performance assessment (PA) of deep geologic nuclear waste repositories, and/or to bring to the fore additional research and development needed to improve PA methodologies. In Task F2- (salt), these objectives have been accomplished through staged development and comparison of the models and methods used by participating teams in their PA frameworks. Coupled-process submodels and deterministic simulations of the entire PA model for a reference scenario for waste disposal in domal salt have been conducted. The task specification has been updated continuously since the initiation of the project to reflect the staged development of the conceptual repository model and performance metrics. Thermal, hydrological, mechanical, and chemical properties of individual components of the engineered and natural system were chosen for relevance by participating teams. The salt reference case system was characterized using data and measurements collected at relevant underground research laboratories (URLs), field sites, and simulation results from teams with specialized modelling capability. Participating teams made a wide range of model assumptions from compartmentalized networks to full 3D models of the salt formation. No single contributed model includes full-fidelity representation of all the features, events, and processes (FEPs) detailed in the task specification, but almost all features and processes are represented in at least one model. Despite differences in the modelling strategies developed by participating teams, all models indicate that salt compaction and radionuclide diffusion are key processes in the repository, and for the FEPs and model scenario considered, little of the disposed radionuclides will migrate beyond the repository seal over the 100,000 year simulations. In general, the model output quantities have the largest differences over the short term and near the waste. The models tend to be more similar further from waste and at later time. Disparities between the models are believed to be due to differing simplifications from the task specification, some of which are chosen simplifications to reduce complexity, and some are restrictions imposed by the modelling tools. A second round of this task has been accepted for DECOVALEX-2027 in conjunction with Task F1 on crystalline PA modelling. The future round includes waste package heating, improved modelling of salt creep closure, additional comparisons of coupled-process sub-models, and the impact of repository engineering design on radionuclide migration in the repository. Participants will also propose and finalize a set of uncertain inputs for the reference case simulations, propagate these uncertainties in a set of realizations, and conduct sensitivity analyses on the simulation results.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Machine learning-assisted upscaling analysis of reservoir rock core properties based on micro-computed tomography imagery

Optimum solutions for geologic modeling and reservoir simulation in industries such as oil and gas recovery and carbon capture and storage require accurate characterization of reservoir properties, which are often heterogeneous. In this study, high-quality micro-computed tomography (CT) images (1.475-μm/pixel resolution) of a sandstone core acquired from the Bell Creek oil field, USA, were used to provide nondestructive analysis of pore- and core-scale heterogeneity across measurement scales of 94–566 μm. In addition to characterizing the as-received sample, the core sample was flooded with brine to evaluate the capacity of the core sample to receive injected fluids. The micro-CT images were systematically segmented into pore spaces and grains via machine learning (ML) steps including image preprocessing, label creation using a traditional ML method based on limited manual image annotation, and finally U-Net segmentation. The segmented image stacks were reconstructed into digital cubes of various scales of voxel lengths. The 3D porosity values were calculated for all the digital cubes, and the fractal dimensions of the cubes were estimated using a box-counting method. The results showed that smaller cubes had greater heterogeneity and that the porosity values could be accurately estimated by fractal dimension and voxel lengths using ML models. For the core sample with brine flooding, the ratio of pores filled by brine to the total pore space was related to the porosity and could also be accurately estimated by porosity, fractal dimension, and voxel lengths using ML models. In conclusion, the results of this study demonstrate that the concept of fractal dimension can be a useful vector to perform upscaling analysis of sandstone rock heterogeneity from the pore to core scale and that fractal dimensions can be used to estimate porosity values and pore space-filling capacity across those scales.

58 GEOSCIENCES↗