Results for “3D Geologic Model” — Page 14

NASA NTRS2019

Astrobiology Survey of a Lava Cave at Lava Beds National Monument by a Rover Carrying a Remote Sensing Instrument Payload

We report here on a survey of a lava tube cave by a rover that is instrumented for astrobiology missions. The NASA Ames testbed rover, CaveR, was deployed in Valentine Cave in Lava Beds National Monument (N. CA, USA) during August of 2018. The rover carried an instrument package consisting of Near Infrared and Visible Spectrometer System (NIRVSS) a point spectrometer operating in 1590-3400 nm range, sensitive to H2O and -OH bearing minerals, pyroxenes, and carbonates (Roush, et al 2018); the bore sighted Drill Operations Camera (DOC), a monochrome imager illuminated by LEDs at 410, 540, 640, 740, 905 and 940 nm; a Realsense ™ depth sensor system for 3D model generation; and a high resolution DSLR stereo camera. The payload was mounted on a tiltable instrument platform attached to the left side of the rover. The rover was driven manually in the cave by field operators, following instructions from a remote science operations team, and simulating a mission concept with science-guided autonomy. A simulated mission took place for 3 days with a team of 3 scientists selecting targets and interpreting data from the payload. To begin the mission, the rover drove along one wall of the cave imaging continuously with the Realsense in 20 m cave segments, three total. At the start of each day, the images from a 20m segment and a panorama stitched from them were provided to the science team to examine. The science team used these data to prioritize specific points along the cave wall for the collection of NIRVSS, DOC, and DSLR data. The objective of the data collection was to identify and study putative biological and mineralogical features in the cave. The data were delivered in xGDS, a customized mapping, planning, and data base management software developed at NASA Ames (Lee, et al 2013). Once the targets for further observations were selected, a plan for collecting the observations (positions in the cave and pointing for each requested observation) was constructed using xGDS and delivered to a rover team to execute the science data collection plan. Acquired data were delivered back to the science team for analysis. Preliminary results from the experiment illustrate the utility of the system (rover plus payload) to study the cave geology and mineralogy and its potential for identifying biomineral features.

Wong, Uland↗

DOE OSTI2014

Utilization of the St. Peter Sandstone in the Illinois Basin for CO2 Sequestration

This project is part of a larger project co-funded by the United States Department of Energy (US DOE) under cooperative agreement DE-FE0002068 from 12/08/2009 through 9/31/2014. The study is to evaluate the potential of formations within the Cambro-Ordovician strata above the Mt. Simon Sandstone as potential targets for carbon dioxide (CO2) sequestration in the Illinois and Michigan Basins. This report evaluates the potential injectivity of the Ordovician St. Peter Sandstone. The evaluation of this formation was accomplished using wireline data, core data, pressure data, and seismic data acquired through funding in this project as well as existing data from two additional, separately funded projects: the US DOE funded Illinois Basin – Decatur Project (IBDP) being conducted by the Midwest Geological Sequestration Consortium (MGSC) in Macon County, Illinois, and the Illinois Industrial Carbon Capture and Sequestration (ICCS) Project funded through the American Recovery and Reinvestment Act (ARRA), which received a phase two award from DOE. This study addresses the question of whether or not the St. Peter Sandstone may serve as a suitable target for CO2 sequestration at locations within the Illinois Basin where it lies at greater depths (below the underground source of drinking water (USDW)) than at the IBDP site. The work performed included numerous improvements to the existing St. Peter reservoir model created in 2010. Model size and spatial resolution were increased resulting in a 3 fold increase in the number of model cells. Seismic data was utilized to inform spatial porosity distribution and an extensive core database was used to develop porosity-permeability relationships. The analysis involved a Base Model representative of the St. Peter at “in-situ” conditions, followed by the creation of two hypothetical models at in-situ + 1,000 feet (ft.) (300 m) and in-situ + 2,000 ft. (600 m) depths through systematic depthdependent adjustment of the Base Model parameters. Properties for the depth shifted models were based on porosity versus depth relationship extracted from the core database followed by application of the porosity-permeability relationship. Each of the three resulting models were used as input to dynamic simulations with the single well injection target of 3.2 million tons per annum (MTPA) for 30 years using an appropriate fracture gradient based bottom hole pressure limit for each injection level. Modeling results are presented in terms of well bottomhole pressure (BHP), injection rate profiles, and three-dimensional (3D) saturation and differential pressure volumes at selected simulation times. Results suggest that the target CO2 injection rate of 3.2 MTPA may be achieved in the St. Peter Sandstone at in-situ conditions and at the in-situ +1,000 ft. (300 m) depth using a single injector well. In the latter case the target injection rate is achieved after a ramp up period which is caused by multi-phase flow effects and thus subject to increased modeling uncertainty. Results confirm that the target rate may not be achieved at the in-situ +2,000 ft. (600 m) level even with multiple wells. These new modeling results for the in-situ case are more optimistic than previous modeling results. This difference is attributed to the difference in methods and data used to develop model permeability distributions. Recommendations for further work include restriction of modeling activity to the in-situ +1,000 ft. (300 m) and shallower depth interval, sensitivity and uncertainty analysis, and refinement of porosity and permeability estimates through depth and area selective querying of the available core database. It is also suggested that further modeling efforts include scope for evaluating project performance in terms of metrics directly related to the Environmental Protection Agency (EPA) Class VI permit requirements for the area of review (AoR) definition and post injection site closure monitoring.

Will, Robert↗

DOE OSTI2022

A Bayesian Approach for In-Situ Stress Prediction and Uncertainty Quantification for Subsurface Engineering

Many subsurface engineering applications require accurate knowledge of the in-situ state of stress for their safe design and operation. Existing methods to meet this need primarily include field measurements for estimating one or more of the principal stresses from a borehole, or optimization methods for constructing a 3D geomechanical model in terms of geophysical measurements. These methods, however, often contain considerable uncertainty in estimating the state of stress. Here, in this paper, we build on a Bayesian approach to quantify uncertainty in stress estimations for subsurface engineering applications. This approach can provide an estimate of the 3D distribution of stress throughout the volume of interest and provide an estimate of the uncertainty arising from the stress measurement, the rheology parameters, and a paucity of measurements. The value of this approach is demonstrated using stress measurements from the In Salah carbon storage site, which was one of the world’s first industrial carbon capture and storage projects. This demonstration shows the application of this Bayesian approach for estimating the initial state of stress for In Salah and quantifying the uncertainty in the estimated stress. Also, an assessment of a maximum injection pressure to prevent geomechanical risks from CO 2 injection pressures is provided in terms of the probability distribution of the minimum principal stress quantified by the approach. With the In Salah case study, this paper demonstrates that using the Bayesian approach can provide additional insights for site explorations and/or project operations to make informed-site decisions for subsurface engineering applications.

58 GEOSCIENCES↗

DOE OSTI2021

2021 Summer SPE Project

SPE is a project in order to develop new, more physics-based, seismic models of explosions (see Nelson et al. 2013). One key component of this effort is numerical modeling enabled by modern state-of-art code/software. Accurate modeling of the shape and amplitude of seismic waves from their generation to their propagation to remote monitoring seismic stations is important to our ability to determine the origin and strength of the source from remote recording. In this project, the modeling is performed by coupling two codes HOSS and SPECFEM3D. HOSS models the dynamic nonlinear processes happening near the explosion. SPECFEM3D computes the propagation of seismic waves as they travel through 3D complex Earth models. SPECFEM3D was modified in order to be driven by a set of time-series calculated by HOSS in lieu of a seismic source. The goal of this summer project is the investigation of several questions pertinent to the establishment of a full end-to-end modeling capability from the high-rate strain regime area to remote distances where seismic station record seismic waves generated by explosions. The investigated questions are: (1) How to perform proper filtering? Direct-solution modeling only sustains a limited range of frequency depending on the grid size. As we go from one modeling domain to the other via coupling, the mesh size is getting coarser to allow modeling at large scale but also to account for the fact that high-frequency waves do not physically travel to large distances. So filtering of the time-series generated by the near-field hydrodynamic modeling is a current practice often employed but its effects on the modeled waveforms has to be investigated. (2) Quantitative assessment of the efficiency of attenuation to remove high-frequency content of the wavefield. This assessment will allow to create meshes with a grid size appropriate to the actual physics of wave propagation for a given explosion. (3) Checking that the coupling process and the two codes respect the cylindrical symmetry that is expected in the case of a pure explosion in a half-space. (4) The effect of the state-of-stress in the near-source area on the modeled seismic waveforms. These questions will be investigated through the modeling of SPE-4P, the fourth explosion of this series because it was designed to have little interaction with the geologic surrounding and the free surface of the Earth so that it is the most explosion-like experiment, with the most symmetries to be verified.

58 GEOSCIENCES↗

DOE OSTI2024

Applicability of 1D site response analysis to shallow sedimentary basins: A critical evaluation through physics-based 3D ground motion simulations

One-dimensional site response analysis (1D SRA) remains the standard practice in considering the effect of local soil deposits and predicting site-specific ground motions, although its range of applicability to realistic seismic wavefields is still in question. In this 1D approach, horizontal and vertical ground shaking are assumed to be induced by vertically propagating shear and compressional waves, respectively. A recent study based on analytical two-dimensional (2D) plane waves and simple point source earthquake simulations has shown two mechanistic limitations in this 1D modelling technique for general inclined seismic waves, that is, systematic over-prediction of the vertical motion and wave trapping in the 1D soil column. In this article, we evaluate in detail the applicability of this 1D modelling approach to realistic three-dimensional (3D) simulated seismic wavefields in shallow sedimentary basins. Linear-viscoelastic 1D SRA predictions using two types of input motions that are commonly used in practice—rock outcrop and in-column motions, are compared with the reference true site response results from 3D earthquake simulations in terms of various measures in the frequency and time domain. It is shown that the horizontal motion in the 3D seismic wavefield exhibits dominant shear wave propagation phenomenon, while the vertical motion is a combined effect of compressional and shear waves and can be over-predicted by the 1D approach when the incident seismic waves are inclined. Direct evidence of the wave refraction process that leads to the vertical motion over-prediction is provided. 1D SRA with in-column inputs can yield motions that have significantly longer duration compared to the true 3D site response solution due to trapped waves, casting in doubt the frequent need for increased soil damping in existing site studies to compensate for wave attenuation due to scattering alone. Sensitivity investigation on the increase of soil profile damping by a multiplier D mul shows D mul values compatible with those found in the literature for both horizontal and vertical motions. It is shown that the level of D mul optimized for a best match of the spectral acceleration is dependent on the characteristic of the input motion and a larger D mul is typically required for the vertical component. In contrast, 1D SRA with outcrop motions predicts motions with shorter significant duration due to its inability to capture the basin-edge generated surface waves. A suite of ground motion simulations was performed to assess the sensitivity of the observations to the basin geologic structure including the velocity gradient, rock-basin impedance contrast and basin depth. The analysis results show that the accuracy of the simplified 1D procedure is dependent on the wavefield composition of both the input motions and the true 3D site response solution. While the horizontal motions in shallow sedimentary basins can, to the first order, be reasonably captured by the simplified 1D approach, 1D SRA for the vertical component is in general not reliable and contributions from inclined shear waves should be accounted for in site-specific evaluation of the vertical design ground motion.

1D assumption↗

DOE OSTI2021

Impacts of Sea‐Level Rise on Hypoxia and Phytoplankton Production in Chesapeake Bay: Model Prediction and Assessment

In this study, the influence of sea-level rise (SLR) on seasonal hypoxia and phytoplankton production in Chesapeake Bay is investigated using a 3D unstructured grid model. Three SLR scenarios (0.17, 0.5, and 1.0 m) were conducted from 1991 to 1995. Results show that the summer hypoxic volume (HV) increases about 2%, 8%, and 16%, respectively, for these three scenarios, compared with Base Scenario. The contributions of physical and biological processes on the increase in the HV were analyzed. With the projected SLR, enhanced gravitational circulation transports more oxygen-rich water in the bottom layer from the mouth. However, the pycnocline moves upwards along with increasing water depth, which largely prolongs the time for dissolved oxygen (DO) to be transported to the bottom. The altered physical processes contribute greatly to a larger HV bay-wide. Besides, SLR increases the whole Bay phytoplankton production, with a larger increase in shallow areas (e.g., 53% in areas with depth <1 m under SLR of 0.5 m). Enhanced light availability is suggested to be the major driver of blooming phytoplankton under SLR in shallow areas. While increased DO production over the euphotic zone is mostly released to the atmosphere and transported downstream, the increase in settled organic matter greatly promotes DO consumption in the water column. The increased respiration is another major cause of the HV increase besides the physical contributions.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2021

Thresholding Analysis and Feature Extraction from 3D Ground Penetrating Radar Data for Noninvasive Assessment of Peanut Yield

This study explores the efficacy of utilizing a novel ground penetrating radar (GPR) acquisition platform and data analysis methods to quantify peanut yield for breeding selection, agronomic research, and producer management and harvest applications. Sixty plots comprising different peanut market types were scanned with a multichannel, air-launched GPR antenna. Image thresholding analysis was performed on 3D GPR data from four of the channels to extract features that were correlated to peanut yield with the objective of developing a noninvasive high-throughput peanut phenotyping and yield-monitoring methodology. Plot-level GPR data were summarized using mean, standard deviation, sum, and the number of nonzero values (counts) below or above different percentile threshold values. Best results were obtained for data below the percentile threshold for mean, standard deviation and sum. Data both below and above the percentile threshold generated good correlations for count. Correlating individual GPR features to yield generated correlations of up to 39% explained variability, while combining GPR features in multiple linear regression models generated up to 51% explained variability. The correlations increased when regression models were developed separately for each peanut type. This research demonstrates that a systematic search of thresholding range, analysis window size, and data summary statistics is necessary for successful application of this type of analysis. The results also establish that thresholding analysis of GPR data is an appropriate methodology for noninvasive assessment of peanut yield, which could be further developed for high-throughput phenotyping and yield-monitoring, adding a new sensor and new capabilities to the growing set of digital agriculture technologies.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2021

Tightly-coupled camera/LiDAR integration for point cloud generation from GNSS/INS-assisted UAV mapping systems

Unmanned aerial vehicles (UAVs) equipped with integrated global navigation satellite systems/inertial navigation systems (GNSS/INS) together with cameras and/or LiDAR sensors are being widely used for topographic mapping in a variety of applications such as precision agriculture, coastal monitoring, and archaeological documentation. Integration of image-based and LiDAR point clouds can provide a comprehensive 3D model of the area of interest. For such integration, ensuring a good alignment between data from the different sources is critical. Although many works have been conducted on this topic, there is still a need for a rigorous integration approach that minimizes the discrepancy between camera and LiDAR data caused by inaccurate system calibration parameters and/or trajectory artifacts. This study proposes an automated tightly-coupled camera/LiDAR integration workflow for GNSS/INS-assisted UAV systems. The proposed strategy is conducted in three main steps. First, an image-based point cloud is generated using a LiDAR/GNSS/INS-assisted structure from motion (SfM) strategy. Then, feature correspondences between image-based and LiDAR point clouds are automatically identified. Finally, an integrated-bundle adjustment procedure including image points, LiDAR raw measurements, and GNSS/INS information is conducted to minimize the discrepancy between point clouds from different sensors while estimating system calibration parameters and refining the trajectory information. The proposed SfM strategy and integration framework are evaluated using five datasets. The SfM results show that using LiDAR data can facilitate feature matching and further increase the number of reconstructed 3D points. The experimental results also illustrate that the developed automated camera/LiDAR integration strategy is capable of accurately estimating system calibration parameters to achieve good alignment among camera/LiDAR data from single/multiple systems. Finally, an absolute accuracy in the range of 3–5 cm is achieved for the image/LiDAR point clouds after the integration process.

42 ENGINEERING↗

DOE OSTI2014

Modeling CO 2 Sequestration in Saline Aquifer and Depleted Oil Reservoir To Evaluate Regional CO 2 Sequestration Potential of Ozark Plateau Aquifer System, South-Central Kansas

1. Drilled, cored, and logged three wells to the basement and collecting more than 2,700 ft of conventional core; obtained 20 m i2 of multicomponent 3D seismic imaging and merged and reprocessed more than 125 mi 2 of existing 3D seismic data for use in modeling CO 2 - EOR oil recovery and CO 2 storage in five oil fields in southern Kansas. 2. Determined the technical feasibility of injecting and sequestering CO 2 in a set of four depleted oil reservoirs in the Cutter, Pleasant Prairie South, Eubank, and Shuck fields in southwest Kansas; of concurrently recovering oil from those fields; and of quantifying the volumes of CO2 sequestered and oil recovered during the process. 3. Formed a consortium of six oil operating companies, five of which own and operate the four fields. The consortium became part of the Southwest Kansas CO 2 -EOR Initiative for the purpose of sharing data, knowledge, and interest in understanding the potential for CO 2 -EOR in Kansas. 4. Built a regional well database covering 30,000 mi 2 and containing stratigraphic tops from ~90,000 wells; correlated 30 major stratigraphic horizons; digitized key wells, including wireline logs and sample logs; and analyzed more than 3,000 drill stem tests to establish that fluid levels in deep aquifers below the Permian evaporites are not connected to the surface and therefore pressures are not hydrostatic. Connectivity with the surface aquifers is lacking because shale aquitards and impermeable evaporite layers consist of both halite and anhydrite. 5. Developed extensive web applications and an interactive mapping system that do the following: a. Facilitate access to a wide array of data obtained in the study, including core descriptions and analyses, sample logs, digital (LAS) well logs, seismic data, gravity and magnetics maps, structural and stratigraphic maps, inferred fault traces, earthquakes, Class I and II disposal wells, and surface lineaments. b. Provide real-time analysis of the project dataset, including automated integration and viewing of well logs, core, core analyses, brine chemistry, and stratigraphy using the Java Profile app. A cross-section app allows for the display of log data for up to four wells at a time. 6. Integrated interpretations from the project’s interactive web-based mapping system to gain insights to aid in assessing the efficacy of geologic CO 2 storage in Kansas and insights toward understanding recent seismicity to aid in evaluating induced vs. naturally occurring earthquakes. 7. Developed a digital type-log system, including web-based software to modify and refine stratigraphic nomenclature to provide stakeholders a common means for communication about the subsurface. 8. Contracted use of a nuclear magnetic resonance (NMR) log and ran it slowly to capture response and characterize larger pores common for carbonate reservoirs. Used NMR to extend core analyses to apply permeability, relative permeability to CO 2 , and capillary pressure to the major rock types, each uniquely expressed as a reservoir quality index (RQI), present in the Mississippian and Arbuckle rocks. 9. Characterized and evaluated the possible role of microbes in dense brines. Used microbes to compliment H/O stable isotopes to fingerprint brine systems. Used perforation/swabbing to obtain samples from multiple hydrostratigraphic units and confirmed equivalent results using less expensive drill stem tests (DST). 10. Used an integrated approach from whole core, logs, tests, and seismic to verify and quantify properties of vuggy, brecciated, and fractured carbonate intervals. 11. Used complex geocellular static and dynamic models to evaluate regional storage capacity using large parallel processing. 12. Carbonates are complex reservoirs and CO 2 -EOR needs to move to the next generation to increase effectiveness of CO 2 and efficiency and safety of the injection.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2014

Modeling CO 2 Sequestration in Saline Aquifer and Depleted Oil Reservoir To Evaluate Regional CO 2 Sequestration Potential of Ozark Plateau Aquifer System, South-Central Kansas

1. Drilled, cored, and logged three wells to the basement and collecting more than 2,700 ft of conventional core; obtained 20 m i2 of multicomponent 3D seismic imaging and merged and reprocessed more than 125 mi 2 of existing 3D seismic data for use in modeling CO 2 - EOR oil recovery and CO 2 storage in five oil fields in southern Kansas. 2. Determined the technical feasibility of injecting and sequestering CO 2 in a set of four depleted oil reservoirs in the Cutter, Pleasant Prairie South, Eubank, and Shuck fields in southwest Kansas; of concurrently recovering oil from those fields; and of quantifying the volumes of CO2 sequestered and oil recovered during the process. 3. Formed a consortium of six oil operating companies, five of which own and operate the four fields. The consortium became part of the Southwest Kansas CO 2 -EOR Initiative for the purpose of sharing data, knowledge, and interest in understanding the potential for CO 2 -EOR in Kansas. 4. Built a regional well database covering 30,000 mi 2 and containing stratigraphic tops from ~90,000 wells; correlated 30 major stratigraphic horizons; digitized key wells, including wireline logs and sample logs; and analyzed more than 3,000 drill stem tests to establish that fluid levels in deep aquifers below the Permian evaporites are not connected to the surface and therefore pressures are not hydrostatic. Connectivity with the surface aquifers is lacking because shale aquitards and impermeable evaporite layers consist of both halite and anhydrite. 5. Developed extensive web applications and an interactive mapping system that do the following: a. Facilitate access to a wide array of data obtained in the study, including core descriptions and analyses, sample logs, digital (LAS) well logs, seismic data, gravity and magnetics maps, structural and stratigraphic maps, inferred fault traces, earthquakes, Class I and II disposal wells, and surface lineaments. b. Provide real-time analysis of the project dataset, including automated integration and viewing of well logs, core, core analyses, brine chemistry, and stratigraphy using the Java Profile app. A cross-section app allows for the display of log data for up to four wells at a time. 6. Integrated interpretations from the project’s interactive web-based mapping system to gain insights to aid in assessing the efficacy of geologic CO 2 storage in Kansas and insights toward understanding recent seismicity to aid in evaluating induced vs. naturally occurring earthquakes. 7. Developed a digital type-log system, including web-based software to modify and refine stratigraphic nomenclature to provide stakeholders a common means for communication about the subsurface. 8. Contracted use of a nuclear magnetic resonance (NMR) log and ran it slowly to capture response and characterize larger pores common for carbonate reservoirs. Used NMR to extend core analyses to apply permeability, relative permeability to CO 2 , and capillary pressure to the major rock types, each uniquely expressed as a reservoir quality index (RQI), present in the Mississippian and Arbuckle rocks. 9. Characterized and evaluated the possible role of microbes in dense brines. Used microbes to compliment H/O stable isotopes to fingerprint brine systems. Used perforation/swabbing to obtain samples from multiple hydrostratigraphic units and confirmed equivalent results using less expensive drill stem tests (DST). 10. Used an integrated approach from whole core, logs, tests, and seismic to verify and quantify properties of vuggy, brecciated, and fractured carbonate intervals. 11. Used complex geocellular static and dynamic models to evaluate regional storage capacity using large parallel processing. 12. Carbonates are complex reservoirs and CO 2 -EOR needs to move to the next generation to increase effectiveness of CO 2 and efficiency and safety of the injection.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2014

Modeling CO 2 Sequestration in Saline Aquifer and Depleted Oil Reservoir To Evaluate Regional CO 2 Sequestration Potential of Ozark Plateau Aquifer System, South-Central Kansas

1. Drilled, cored, and logged three wells to the basement and collecting more than 2,700 ft of conventional core; obtained 20 m i2 of multicomponent 3D seismic imaging and merged and reprocessed more than 125 mi 2 of existing 3D seismic data for use in modeling CO 2 - EOR oil recovery and CO 2 storage in five oil fields in southern Kansas. 2. Determined the technical feasibility of injecting and sequestering CO 2 in a set of four depleted oil reservoirs in the Cutter, Pleasant Prairie South, Eubank, and Shuck fields in southwest Kansas; of concurrently recovering oil from those fields; and of quantifying the volumes of CO2 sequestered and oil recovered during the process. 3. Formed a consortium of six oil operating companies, five of which own and operate the four fields. The consortium became part of the Southwest Kansas CO 2 -EOR Initiative for the purpose of sharing data, knowledge, and interest in understanding the potential for CO 2 -EOR in Kansas. 4. Built a regional well database covering 30,000 mi 2 and containing stratigraphic tops from ~90,000 wells; correlated 30 major stratigraphic horizons; digitized key wells, including wireline logs and sample logs; and analyzed more than 3,000 drill stem tests to establish that fluid levels in deep aquifers below the Permian evaporites are not connected to the surface and therefore pressures are not hydrostatic. Connectivity with the surface aquifers is lacking because shale aquitards and impermeable evaporite layers consist of both halite and anhydrite. 5. Developed extensive web applications and an interactive mapping system that do the following: a. Facilitate access to a wide array of data obtained in the study, including core descriptions and analyses, sample logs, digital (LAS) well logs, seismic data, gravity and magnetics maps, structural and stratigraphic maps, inferred fault traces, earthquakes, Class I and II disposal wells, and surface lineaments. b. Provide real-time analysis of the project dataset, including automated integration and viewing of well logs, core, core analyses, brine chemistry, and stratigraphy using the Java Profile app. A cross-section app allows for the display of log data for up to four wells at a time. 6. Integrated interpretations from the project’s interactive web-based mapping system to gain insights to aid in assessing the efficacy of geologic CO 2 storage in Kansas and insights toward understanding recent seismicity to aid in evaluating induced vs. naturally occurring earthquakes. 7. Developed a digital type-log system, including web-based software to modify and refine stratigraphic nomenclature to provide stakeholders a common means for communication about the subsurface. 8. Contracted use of a nuclear magnetic resonance (NMR) log and ran it slowly to capture response and characterize larger pores common for carbonate reservoirs. Used NMR to extend core analyses to apply permeability, relative permeability to CO 2 , and capillary pressure to the major rock types, each uniquely expressed as a reservoir quality index (RQI), present in the Mississippian and Arbuckle rocks. 9. Characterized and evaluated the possible role of microbes in dense brines. Used microbes to compliment H/O stable isotopes to fingerprint brine systems. Used perforation/swabbing to obtain samples from multiple hydrostratigraphic units and confirmed equivalent results using less expensive drill stem tests (DST). 10. Used an integrated approach from whole core, logs, tests, and seismic to verify and quantify properties of vuggy, brecciated, and fractured carbonate intervals. 11. Used complex geocellular static and dynamic models to evaluate regional storage capacity using large parallel processing. 12. Carbonates are complex reservoirs and CO 2 -EOR needs to move to the next generation to increase effectiveness of CO 2 and efficiency and safety of the injection.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2014

Modeling CO 2 Sequestration in Saline Aquifer and Depleted Oil Reservoir To Evaluate Regional CO 2 Sequestration Potential of Ozark Plateau Aquifer System, South-Central Kansas

1. Drilled, cored, and logged three wells to the basement and collecting more than 2,700 ft of conventional core; obtained 20 m i2 of multicomponent 3D seismic imaging and merged and reprocessed more than 125 mi 2 of existing 3D seismic data for use in modeling CO 2 - EOR oil recovery and CO 2 storage in five oil fields in southern Kansas. 2. Determined the technical feasibility of injecting and sequestering CO 2 in a set of four depleted oil reservoirs in the Cutter, Pleasant Prairie South, Eubank, and Shuck fields in southwest Kansas; of concurrently recovering oil from those fields; and of quantifying the volumes of CO2 sequestered and oil recovered during the process. 3. Formed a consortium of six oil operating companies, five of which own and operate the four fields. The consortium became part of the Southwest Kansas CO 2 -EOR Initiative for the purpose of sharing data, knowledge, and interest in understanding the potential for CO 2 -EOR in Kansas. 4. Built a regional well database covering 30,000 mi 2 and containing stratigraphic tops from ~90,000 wells; correlated 30 major stratigraphic horizons; digitized key wells, including wireline logs and sample logs; and analyzed more than 3,000 drill stem tests to establish that fluid levels in deep aquifers below the Permian evaporites are not connected to the surface and therefore pressures are not hydrostatic. Connectivity with the surface aquifers is lacking because shale aquitards and impermeable evaporite layers consist of both halite and anhydrite. 5. Developed extensive web applications and an interactive mapping system that do the following: a. Facilitate access to a wide array of data obtained in the study, including core descriptions and analyses, sample logs, digital (LAS) well logs, seismic data, gravity and magnetics maps, structural and stratigraphic maps, inferred fault traces, earthquakes, Class I and II disposal wells, and surface lineaments. b. Provide real-time analysis of the project dataset, including automated integration and viewing of well logs, core, core analyses, brine chemistry, and stratigraphy using the Java Profile app. A cross-section app allows for the display of log data for up to four wells at a time. 6. Integrated interpretations from the project’s interactive web-based mapping system to gain insights to aid in assessing the efficacy of geologic CO 2 storage in Kansas and insights toward understanding recent seismicity to aid in evaluating induced vs. naturally occurring earthquakes. 7. Developed a digital type-log system, including web-based software to modify and refine stratigraphic nomenclature to provide stakeholders a common means for communication about the subsurface. 8. Contracted use of a nuclear magnetic resonance (NMR) log and ran it slowly to capture response and characterize larger pores common for carbonate reservoirs. Used NMR to extend core analyses to apply permeability, relative permeability to CO 2 , and capillary pressure to the major rock types, each uniquely expressed as a reservoir quality index (RQI), present in the Mississippian and Arbuckle rocks. 9. Characterized and evaluated the possible role of microbes in dense brines. Used microbes to compliment H/O stable isotopes to fingerprint brine systems. Used perforation/swabbing to obtain samples from multiple hydrostratigraphic units and confirmed equivalent results using less expensive drill stem tests (DST). 10. Used an integrated approach from whole core, logs, tests, and seismic to verify and quantify properties of vuggy, brecciated, and fractured carbonate intervals. 11. Used complex geocellular static and dynamic models to evaluate regional storage capacity using large parallel processing. 12. Carbonates are complex reservoirs and CO 2 -EOR needs to move to the next generation to increase effectiveness of CO 2 and efficiency and safety of the injection.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2018

Assessing the carbon sequestration potential of basalt using X-ray micro-CT and rock mechanics

Mineral carbonation in basaltic rock provides a permanent storage solution for the mitigation of anthropogenic CO 2 emissions in the atmosphere. 3D X-ray micro-CT (XCT) image analysis is applied to a core sample from the main basaltic reservoir of the CarbFix site in Iceland, which obtained a connected porosity of 2.05–8.76%, a reactive surface area of 0.10–0.33 mm -1 and a larger vertical permeability (2.07 × 10 -10 m 2 ) compared to horizontal permeability (5.10 × 10 -11 m 2 ). The calculations suggest a CO 2 storage capacity of 0.33 Gigatonnes at the CarbFix pilot site. The XCT results were compared to those obtained from a hydromechanical test applied to the same sample, during which permeability, electrical resistivity and volumetric deformation were measured under realistic reservoir pressure conditions. It was found that permeability is highly stress sensitive, dropping by two orders of magnitude for a -0.02% volumetric deformation change, equivalent to a mean pore diameter reduction of 5 μm. This pore contraction was insufficient to explain such a permeability reduction according to the XCT analysis, unless combined with the effects of clay swelling and secondary mineral pore clogging. Here the findings provide important benchmark data for the future upscaling and optimisation of CO 2 storage in basalt formations.

58 GEOSCIENCES↗

NASA NTRS2022

Simulating Mars: Enabling Testing of the Perseverance Rover Sampling and Caching Subsystem on Earth

The development of the Sampling and Caching Subsystem (SCS) on the JPL Perseverance Rover lies at the intersection of testing, robotics, and geology. The SCS team established three primary system test campaigns and venues to aid in the development of SCS through verification and validation testing – Qualification Model Dirty Testing (QMDT) to provide a venue for testing in a Martian environment, Vehicle System Testbed (VSTB) for testing while integrated with the mobility subsystem on Martian-like terrain, and the Flight Software Testbed (FSWTB) for conducting tests using the flight motor controllers and software system on a hexapod which had the ability to simulate rover tilt. Each venue contributed a vital piece to the SCS building blocks. However, the QMDT venue operating within a 10-ft diameter Thermal Vacuum chamber to simulate Martian environment provided a sui generis opportunity to fine tune the entire sampling and caching process while building the team’s knowledge base about rock drillability, system life, and target selection. On Earth, because Martian rocks are not readily available, the development team must utilize geoanalogs to the rocks and regolith on Mars. Geologists on the team helped establish a set of standard rock types to use for Mars missions, like Basalt, Sandstone, Mudstone, Gypsum, and other related geoanalogs. These geoanalogs are characterized with a standard suite of tests for density, compressibility, and other characteristics to categorize potential drillability. This concept of drillability is what links the geoanalogs on Earth to the samples we collect on Mars. With the simulant characteristics defined, these geoanalog rocks are ready to be drilled into as we do on the Martian surface. A key aspect of interacting with the surface on Mars is rock target identification and selection. The Perseverance robotic system uses the on-board cameras, instrumentation, and software to collect enough information to identify potential scientific targets. With the targets identified, SCS can place the Corer and abrade the surface or collect a sample. For a ground test activity like QMDT, the test team did not have all of the camera and instrumentation systems that the rover does, so the team developed ground test equivalents to process a rock, build a target map, and define the target. The team constructed a Rock Scanning Station to build a 3D point cloud of the rock. This point cloud was then processed and evaluated with predefined and programmed criteria in a Target Downselect Tool. A primary output of the Target Downselect Tool is a defined target that can be uploaded directly to the robotic software system to simulate and build the robotic sequences used in tests. With these insights and programmatic definition of targets, the QMDT test team was able to make the same decisions that the Perseverance surface operations team does. In addition, valuable lessons learned from developing the target selection ground tools and using them were implemented into the tools used for surface operations.

Kim, Junggon↗

DOE OSTI2017

Preface to the Special Issue on TOUGH Symposium 2015

The TOUGH Symposium 2015 was held in Berkeley, California, September 28-30, 2015. The TOUGH family of codes, developed at the Energy Geosciences Division of Lawrence Berkeley National Laboratory (LBNL), is a suite of computer programs for the simulation of multiphase and multicomponent fluid and heat flows in porous and fractured media with applications in many geosciences fields, such as geothermal reservoir engineering, nuclear waste disposal, geological carbon sequestration, oil and gas reservoirs, gas hydrate research, vadose zone hydrology and environmental remediation. Since the first release in the 1980s, many modifications and enhancements have been continuously made to TOUGH and its various descendants (iTOUGH2, TOUGH+, TOUGH-MP, TOUGHREACT, TOUGH+HYDRATE, TMVOC...), at LBNL and elsewhere. Today, these codes are used worldwide in academia, government organizations and private companies in problems involving coupled hydrological, thermal, biogeochemical and geomechanical processes. The Symposia, organized every 2-3 years, bring together developers and users for an open exchange on recent code enhancements and applications. In 2015, the Symposium was attended by one hundred participants, representing thirty-four nationalities. This Special Issue in Computers & Geosciences gathers extended versions of selected Symposium proceedings related to (i) recent enhancements to the TOUGH family of codes and (ii) coupled flow and geomechanics processes modeling. In this special issue, a paper by Jung et al. provides a new base version of TOUGH, TOUGH3. This latest version addresses the computational robustness and efficiency needed for the increasing complexity of the simulated processes and the growing size of model domains. Moreover, it consolidates the serial and parallel versions into a single code, hence, inconsistencies between different versions are less likely to occur, and code maintenance and development become more efficient. Finsterle et al. present a review of the capabilities of iTOUGH2 and illustrate the workflow of an integrated simulation-optimization framework. Liu et al. present an approach to reduce the computational cost of performing Bayesian inverse modeling and demonstrate their approach using a synthetic experiment simulated with TOUGH2. Regarding the mesh generation process, the paper by Doughty presents a generalization of the grid-generation algorithm embedded in TOUGH2 (MESHMAKER) to create one-column grids with fractal dimension, useful for representing flow through fractures or highly heterogeneous media with fractal geometry, and also for inverse methods. Sentís and Gable present a module to interface the open source code Lagrit with TOUGH2, that assists in the processes of mesh generation and model setup, and that can be particularly useful for complex geometries and large models. Bonduà et al. present a pre-processing tool (VORO2MESH) that features a full 3D Voronoi approach for structured and unstructured grids, and an enhanced post-processor (TOUGH2Viewer) for the visualization of 3D Voronoi grids. The literature review shows that coupled fluid flow and geomechanics has gained increasing attention in the last years in many geosciences fields. In this context, several codes of the TOUGH family have been successfully coupled with geomechanics software. The paper by Rutqvist presents a review of some of these TOUGH-based geomechanics models, of which the TOUGH-FLAC simulator – presented for the first time in the TOUGH Symposium 2003 – is well-known worldwide and has been used in a large variety of applications, encompassing the application domains of the TOUGH codes. Blanco-Martín et al. present an extension of TOUGH-FLAC to the finite strain framework, motivated by the existence of large strains associated with rock salt deformation in the context of heat-generating nuclear waste disposal. Following a paper published in the Special Issue of the TOUGH Symposium 2012 in this journal, Kim et al. present the current status of the TOUGH-RSBN simulator for coupled hydro-mechanical modeling of hydraulic fracturing using a discrete fracture network approach, and present some validation and demonstration cases. Rinaldi and Nespoli present TOUGH-SEED, an injection-induced seismicity modeling approach using TOUGH2 and a geomechanics stochastic model, and apply it to the Basel EGS project. Finally, Rinaldi et al. use the PEST protocol implemented in iTOUGH2, combined with TOUGH-FLAC as the forward model, to perform inverse modeling for the In Salah CO2 demonstration site in Algeria. Using pore pressures and InSAR displacements as observations, this approach allows for a multi-physics joint inversion and takes advantage of the inverse capabilities of iTOUGH2 (parameter estimation, sensitivity analysis, uncertainty estimation …), applied to coupled fluid flow and geomechanics. On behalf of the members of the TOUGH Symposium 2015 Organizing Committee, I hope that this collection of papers will provide the readers of Computers & Geosciences with a recent view of the capabilities of the TOUGH family of codes, as well as with current modeling research domains in the geosciences. Last but not least, I would like to express my gratitude to all the reviewers who kindly contributed to this Special Issue, as well as the Journal staff who have made this publication possible.

99 GENERAL AND MISCELLANEOUS↗

DOE OSTI2014

Influence of local capillary trapping on containment system effectiveness

Immobilization of CO 2 injected into deep subsurface storage reservoirs is a critical component of risk assessment for geologic CO 2 storage (GCS). Local capillary trapping (LCT) is a recently established mode of immobilization that arises when CO 2 migrates due to buoyancy through heterogeneous storage reservoirs. This project sought to assess the amount and extent of LCT expected in storage formations under a range of injection conditions, and to confirm the persistence of LCT if the seal overlying the reservoir were to lose its integrity. Numerical simulation using commercial reservoir simulation software was conducted to assess the influence of injection. Laboratory experiments, modeling and numerical simulation were conducted to assess the effect of compromised seal integrity. Bench-scale (0.6 m by 0.6 m by 0.03 m) experiments with surrogate fluids provided the first empirical confirmation of the key concepts underlying LCT: accumulation of buoyant nonwetting phase at above residual saturations beneath capillary barriers in a variety of structures, which remains immobile under normal capillary pressure gradients. Immobilization of above-residual saturations is a critical distinction between LCT and the more familiar “residual saturation trapping.” To estimate the possible extent of LCT in a storage reservoir an algorithm was developed to identify all potential local traps, given the spatial distribution of capillary entry pressure in the reservoir. The algorithm assumes that the driving force for CO 2 migration can be represented as a single value of “critical capillary entry pressure” P c,entry crit , such that cells with capillary entry pressure greater/less than P c,entry crit act as barriers/potential traps during CO 2 migration. At intermediate values of P c,entry crit , the barrier regions become more laterally extensive in the reservoir, approaching a percolation threshold while non-barrier regions remain numerous. The maximum possible extent of LCT thus occurs at P c,entry crit near this threshold. Testing predictions of this simple algorithm against full-physics simulations of buoyancy-driven CO 2 migration support the concept of critical capillary entry pressure. However, further research is needed to determine whether a single value of critical capillary entry pressure always applies and how that value can be determined a priori. Simulations of injection into high-resolution (cells 0.3 m on a side) 2D and 3D heterogeneous domains show two characteristic behaviors. At small gravity numbers (vertical flow velocity much less than horizontal flow velocity) the CO 2 fills local traps as well as regions that would act as local barriers if CO 2 were moving only due to buoyancy. When injection ceases, the CO 2 migrates vertically to establish large saturations within local traps and residual saturation elsewhere. At large gravity numbers, the CO 2 invades a smaller portion of the perforated interval. Within this smaller swept zone the local barriers are not invaded, but local traps are filled to large saturation during injection and remain during post-injection gravity-driven migration. The small gravity number behavior is expected in the region within 100 m of a vertical injection well at anticipated rates of injection for commercial GCS. Simulations of leakage scenarios (through-going region of large permeability imposed in overlying seal) indicate that LCT persists (i.e. CO 2 remains held in a large fraction of the local iv traps) and the persistence is independent of injection rate during storage. Simulations of leakage for the limiting case of CO 2 migrating vertically from an areally extensive emplacement in the lower portion of a reservoir showed similar strong persistence of LCT. This research has two broad implications for GCS. The first is that LCT can retain a significant fraction of the CO 2 stored in a reservoir – above and beyond the residual saturation -- if the overlying seal were to fail. Thus frameworks for risk assessment should be extended to account for LCT. The second implication is that compared to pressure driven flow in reservoirs, CO 2 migration and trapping behave in a qualitatively different manner in heterogeneous reservoirs when buoyancy is the dominant driving force for flow. Thus simulations of GCS that neglect capillary heterogeneity will fail to capture important features of the CO 2 plume. While commercial reservoir simulation software can account for fine scale capillary heterogeneity, it has not been designed to work efficiently with such domains, and no simulators can handle fine-scale resolution throughout the reservoir. A possible way to upscale the migration and trapping is to apply an “effective residual saturation” to coarse-scale grids. While the extent of overall immobilization can be correlated in this way, all coarser grids failed to capture the distance traveled by the migrating CO 2 for large gravity number. Thus it remains unclear how best to account for LCT in the routine simulation work-flow that will be needed for large-scale GCS. Alternatives meriting investigation include streamline methods, reduced-physics proxies (e.g. particle tracking), and biased invasion percolation algorithms, which are based on precisely the capillary heterogeneity essential for LCT.

58 GEOSCIENCES↗

NASA NTRSDate not supplied

VERNE: Revealing the Mysteries and Histories of Venus

Introduction: The three Venus missions that were recently selected for upcoming flight (VERITAS, DAVINCI+, and EnVision) will be incredibly valuable to our understanding of Venus’ history, geology, and atmosphere. However, even once completed, key gaps in our knowledge of Venus, and more generally the formation and active processes on rocky, Earth-like planets, will still persist. Remaining questions include 1) how global intrinsic magnetic fields might be maintained on rocky worlds, and how they could then go extinct, and 2) what role atmospheric sulfur chemistry plays in climates of Earth-like planets, which is an increasingly timely subject as Earth’s own atmospheric sulfur content is climbing due to human activity. These questions require in-situ observations from Venus’ cloud deck, at the altitudes at which the UV absorber exists. The Venus Environment Research and Novel Exploration (VERNE) mission will address these questions with an aerial platform that will drift around the equatorial region of the planet for 9 days. VERNE will collect data to determine the identity of the mysterious UV absorber, while also taking magnetic field measurements over the tesserae, the regions on Venus that are most likely to retain remanent crustal magnetization, in order to understand the potential role of a past intrinsically-generated global magnetic field on Venus. Mission Objectives: The two major science objectives that drive the VERNE mission are: 1) Determine the identity of the Venusian unknown ultraviolet absorber(s). First observed approximately a century ago [1], the composition of Venus’ ultraviolet (UV) absorber is one of the oldest mysteries in Venus atmospheric chemistry [2,3]. While several candidate UV absorbers (mostly sulfur species) have been proposed, no consensus has been reached on its composition and its specific interactions with the atmosphere. Determining the identity of the UV absorber will aid climate models by showing how and where incident solar energy is absorbed by the atmosphere and make progress towards understanding the chemical and energetic processes taking place above Venus’ upper cloud deck [4]. 2) Determine if Venus retains evidence of a past, internally-generated magnetic field. While Venus does not currently have an intrinsically-generated magnetic field, evidence for the existence of a past field on Venus and a timeline of its decay will fill in a more holistic picture of the evolution of Venus’ geological record and atmosphere. As the oldest geologic units on the surface, Venus’ tesserae may still have remanent crustal magnetization signatures within the rocky composition [5]. In any case, the signatures detected will provide insight into Venus’ past geological and core dynamo activity. Mission Summary: VERNE includes a 3-part flight system made up of 1) an entry system with a HEEET (Heatshield for Extreme Entry Environment Technology) aeroshell, 2) an orbiter for data relay, and 3) an in-situ balloon and gondola. After entry into Venus’ atmosphere, the balloon will be deployed within the upper cloud deck, at an altitude of 62 km, over a tesserae region. The in-situ data collection will last for the duration of 2 full circumnavigations of the planet, which will take ~9 days. Instrumentation Suite: The four instruments that comprise VERNE’s instrument payload will enable the identification of the unknown UV absorber and the detection of remanent crustal magnetization if it exists in the tesserae. The proposed instrument suite cycle during mission operations is shown in Fig. 1. The four instruments are described below: 1) Adams (ion neutral mass spectrometer) has a range of 18-257 AMU and will detect and distinguish the mixing ratios of O₂, H₂O, H₂SO₄, S, S₂, S₃, S₄, S₅, S₆, S₇, S₈, SO, SO₂, OSSO, SO₃, Cl₂, FeCl₃ and other trace sulfur and organic species. The spatial (longitudinal) and temporal (day/night) variations will be observed throughout 2 circumnavigations with a sample cadence of 12 minutes. 2) Shelley (nephelometer) will determine the size distribution of the aerosols (0.4 to 36 um) in the atmosphere. With a size resolution of <0.7 um, it can determine which mode of H2SO4 is present and characterize the large (>30 um) organic particles previously detected by the Venera and Galileo missions [6,7]. 3) Herbert (UV imager) will measure UV radiance at 283 nm (the wavelength of SO2 absorption) and 365 nm (the unknown part of the absorber). UV images will be taken concurrently with the INMS and nephelometer to correlate UV absorption with abundances of the UV absorber species. 4) Vonnegut (magnetometer) has a range of >600 nT and a precision and accuracy of 1 nT. If magnetized by a past field, the crust may have retained a magnetization of up to 3 A/m2 [5]. With a noise floor of 10 nT, the magnetometer will be able to detect RCM from an altitude of 62 km, even if the thickness of the magnetized crust is just 1 km (Fig. 2). Mission Concept Design: VERNE will be launched with a mass of 3300 kg in an intermediate-high performance class vehicle with a 4-m fairing. The 475-day mission includes 466 days for the cruise, coasting, and orbit initialization phases before entry, descent, and balloon deployment. During the 9-day science phase, the aerial platform will make 2 circumnavigations of the planet at an altitude of 62 km. The INMS and the nephelometer will acquire data for 2 hours during the daytime and nighttime during each circumnavigation, while the UV imager will be on for the duration of the daytime, and the magnetometer will be operational throughout the entirety of the science phase. Data will be stored and processed with the JPL-designed Sphinx command and data handling system. Data will be sent from the balloon to the orbiter using an S-band relay link, stored on the orbiter, and then forwarded to Earth where it will be received by the DSN. Conclusion: VERNE will fill key gaps in our understanding of the history and ongoing processes related to the geology and atmosphere of Venus and rocky worlds in general. Even with adequate flight system contingencies and expected costs below the $900M New Frontiers cost cap, VERNE is not without its risks and challenges. Further trade spaces to explore include 1) using solely battery power vs. including solar panels to increase the mission duration and 2) investigating the use of lightweight materials and 3D-printed structures to reduce the gondola mass, among others. Acknowledgments: We would like to thank the JPL Planetary Science Summer School, especially our mentors Troy Hudson, Karl Mitchell, and Leslie Lowes, as well as our Team-X study lead Al Nash and the members of Team-X. We’d additionally like to thank our review panel for asking insightful questions and providing valuable feedback. References: [1] Ross, F. E. (1928) Astrophysical J., 68, 57-92. [2] Rossow, W. B. et al. (1980) J. Geophysical Research, 85, 8107-8128. [3] Pinto, J. P. et al. (2021) Nature Communications, 12, 175. [4] Titov, D. V. et al. (2007) Cosmic Research, 21, 401. [5] O’Rourke, J. et al. (2019) Geophysical Research Lett., 46, 5768–5777. [6] Limaye, S. S. et al. (2018) Astrobiology, 18(9), 1181-1198. [7] Grinspoon, D. H. et al. (2013) Planetary and Space Sci., 41(7), 515-542. [8] Parker, R. L. (2003) J. Geophysical Research, 108, 5006. *The cost information contained in this document is of a budgetary and planning nature and is intended for informational purposes only. It does not constitute a commitment on the part of JPL and/or Caltech.

H Alpert↗