Integrating Improved Neural Operators and Graph Convolutional Networks for Scalable Geological Carbon Storage Modeling
Explore the source record for details and available documents.
SEARCH · Engineering Papers
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.
Explore the source record for details and available documents.
Explore the source record for details and available documents.
Crustal motions in Alaska are computed by calculating the finite element solution for an elastic spherical shell problem. Very long baseline interferometry and geologic data are used and seem to be in rough agreement.
Explore the source record for details and available documents.
Explore the source record for details and available documents.
These resources describe the 3D geophysical inversion modeling of gravity data at the FORGE site near Milford, Utah. FORGE is the Frontier Observatory for Research in Geothermal Energy and the site in Utah has been selected by the U.S. Dept. of Energy for a 5-year R&D program to test technologies for the development of Engineered Geothermal Systems (EGS). 3D modelling of gravity data at the FORGE site is to help characterize the subsurface geologic framework. Specifically, modelling of gravity data in 3D, used in conjunction with rock density measurements and other subsurface geologic information can provide an independent test of an existing 3D geologic model (e.g. Witter et al., 2018). Such an exercise can be useful for reducing uncertainty in 3D geologic models (Witter et al, 2019).
Hydrogeological field studies rely often on a single conceptual representation of the subsurface. This is problematic since the impact of a poorly chosen conceptual model on predictions might be significantly larger than the one caused by parameter uncertainty. Furthermore, conceptual models often need to incorporate geological concepts and patterns in order to provide meaningful uncertainty quantification and predictions. Consequently, several geologically realistic conceptual models should ideally be considered and evaluated in terms of their relative merits. Here, we propose a full Bayesian methodology based on Markov chain Monte Carlo to enable model selection among 2-D conceptual models that are sampled using training images and concepts from multiple-point statistics. More precisely, power posteriors for the different conceptual subsurface models are sampled using sequential geostatistical resampling and Graph Cuts. To demonstrate the methodology, we compare and rank five alternative conceptual geological models that have been proposed in the literature to describe aquifer heterogeneity at the MAcroDispersion Experiment site in Mississippi, USA. We consider a small-scale tracer test for which the spatial distribution of hydraulic conductivity impacts multilevel solute concentration data observed along a 2-D transect. The thermodynamic integration and the stepping-stone sampling methods were used to compute the evidence and associated Bayes factors using the computed power posteriors. We find that both methods are compatible with multiple-point statistics-based inversions and provide a consistent ranking of the competing conceptual models considered.
In many hydrothermal systems, fracture permeability along faults provides pathways for groundwater to transport heat from depth. Faulting generates a range of deformation styles that cross-cut heterogeneous geology, resulting in complex patterns of permeability, porosity, and hydraulic conductivity. Vertical connectivity (a through going network of permeable areas that allows advection of heat from depth to the shallow subsurface) is rare and is confined to relatively small volumes that have highly variable spatial distribution. This local compartmentalization of connectivity represents a significant challenge to understanding hydrothermal circulation and for exploring, developing, and managing hydrothermal resources. Here, we present an evaluation of the geologic characteristics that control this compartmentalization in hydrothermal systems through 3-D analysis of the Brady geothermal field in western Nevada. A published 3-D geologic map of the Brady area is used as a basis to develop structural and geological variables that are hypothesized to control or effect permeability or connectivity. The 3-D distribution of these variables is compared to the distribution of productive and non-productive fluid flow intervals along production wells and non-productive wells via principal component analysis (PCA). This comparison elucidates which geologic and structural variables are most closely associated with productive fluid flow intervals. Results indicate that production intervals at Brady are located: (1) within or near to known and stress-loaded macro-scale faults, and (2) in areas of high fault and fracture density. This submission includes the published journal article detailing this work, the published 3-D geologic map of the Brady Geothermal Area used as a basis to develop structural and geological variables that are hypothesized to control or effect permeability or connectivity, 3-D well data, along which geologic data were sampled for PCA analyses, and associated metadata file. This work was done using existing R programs.
The distribution of lunar geologic units in space and time and their mode of origin were considered since they provided significant data which bear on a number of current problems in lunar geophysics. Observations and problems were discussed which deal with the characterization of the upper 25 km of the lunar crust, the tectonic style of the crust, the formation of mascons within major basins, analysis of lunar magnetic anomalies, and the history of the lunar crust.
Colorado School of Mines conducted research and training in the development and validation of an advanced CO 2 GS (Geological Sequestration) probabilistic simulation and risk assessment model. CO 2 GS simulation and risk assessment is used to develop advanced numerical simulation models of the subsurface to forecast CO 2 behavior and transport; optimize site operational practices; ensure site safety; and refine site monitoring, verification, and accounting efforts. As simulation models are refined with new data, the uncertainty surrounding the identified risks decrease, thereby providing more accurate risk assessment. The models considered the full coupling of multiple physical processes (geomechanical and fluid flow) and describe the effects of stochastic hydro-mechanical (H-M) parameters on the modeling of CO 2 flow and transport in fractured porous rocks. Graduate students were involved in the development and validation of the model that can be used to predict the fate, movement, and storage of CO 2 in subsurface formations, and to evaluate the risk of potential leakage to the atmosphere and underground aquifers. The main major contributions from the project include the development of: 1) an improved procedure to rigorously couple the simulations of hydro-thermo-mechanical (H-M) processes involved in CO 2 GS; 2) models for the hydro-mechanical behavior of fractured porous rocks with random fracture patterns; and 3) probabilistic methods to account for the effects of stochastic fluid flow and geomechanical properties on flow, transport, storage and leakage associated with CO 2 GS. The research project provided the means to educate and train graduate students in the science and technology of CO 2 GS, with a focus on geologic storage. Specifically, the training included the investigation of an advanced CO 2 GS simulation and risk assessment model that can be used to predict the fate, movement, and storage of CO 2 in underground formations, and the evaluation of the risk of potential CO 2 leakage to the atmosphere and underground aquifers.
This project conducted a modeling study on seismic invisibility cloaks that render geologic targets invisible to seismic waves, using the concept of seismic metamaterials. We present a parametric numerical study on the behaviors of seismic waves through cloaks with different design parameters as well as degrees of geologic heterogeneity. In addition, a seismic cloaking strategy is proposed for a future field-scale experiment at a real-world test bed. This feasibility study will guide future field experiment designs and ultimately allow us to conduct systematic field-scale tests employing Sandia’s existing resources and field expertise. The ultimate goal is to develop methods and design parameters of seismic invisibility cloaks to protect against natural and man-made seismic waves. Seismic cloaking has potential applications in several areas of national security, energy, and natural hazard reduction.
In Rock Valley, Nevada, the USA, past earthquake sequences and shallow faulting in Paleozoic carbonates remain poorly understood. Carbonates are typically more ductile rocks that do not experience large stress drops or fracture coalescence, which contradicts previous observations in the region. As part of the Source Physics Experiment, an effort has been made to experimentally characterize the petrophysical and geomechanical properties of carbonates from Rock Valley. Well core and outcrop samples were used to determine the difference between shallow-buried Tertiary limestones and deeply buried Paleozoic limestones and dolostones. The carbonates possessed porosities between 3 and 9%, with V P and V S in the Tertiary carbonates around 2700 and 1800 m/s, respectively, and 6000 and 3100 m/s, respectively, in the Paleozoic carbonates. Microstructural characterization revealed that the Paleozoic carbonates contained significantly more pre-existing damage and alteration than the Tertiary carbonates, though the deformation varies from localized to diffuse. Unconfined Brazilian tests and triaxial tests with confining pressures between 0 and 50 MPa showed that, although samples all experienced failure, the dolostones typically failed at greater stresses and possessed greater E/ν ratios than the limestones. Furthermore, the Hoek–Brown failure criterion was used to construct failure envelopes with the compressive and tensile failure tests. Velocity, density, and porosity measurements were used to construct a hypothetical velocity-depth profile and compare the results with the theoretical petrophysical measurements. Using brittleness analysis, the likely failure conditions were determined to be low-porosity dolomitic rock for fault nucleation in the Paleozoic basement.
A reservoir modeling study was conducted to assess the feasibility of storing 162.5 million tonnes (8.125 million tonnes annually) of industrially sourced carbon dioxide (CO 2 ) in the St. Peter – Everton and Knox storage complexes at the Prairie State Generating Company’s (PSGC) site near Marissa, Washington County, Illinois (PSGC site). Two separate models were constructed for the St. Peter – Everton and the Knox storage complexes. The St. Peter and Everton sandstones are the target storage units for the St. Peter – Everton storage complex. The Knox Group formations are the target storage units for the Knox storage complex. The Maquoketa Shale is the primary confining unit for both storage complexes. The static reservoir models used as input to the dynamic reservoir models of the storage complexes were developed from interpretations of logs (e.g., gamma ray, resistivity, porosity, photoelectric, and sonic), structure surfaces using data from 36 wells, thickness maps, well test data, seismic data, and permeability data in PetrelTM.
This Storage Field Development plan presents the Storage Complex characterization results, construction, monitoring, operational plans, and costs associated with the proposed Prairie State Energy Campus - Carbon Capture and Storage (PSEC-CCS) project in Washington County, Illinois, near Marissa. The proposed storage complex, known as the St. Peter-Everton Storage Complex, comprises the St. Peter Sandstone and Everton sandstone reservoirs, with the Maquoketa Shale serving as the primary confining unit and the Knox Group as lower confining unit. The lowermost Underground Source of Drinking Water (USDW) identified for the site is the Pennsylvanian section.
Explore the source record for details and available documents.
Deep learning-based surrogate models have been widely applied in geological carbon storage (GCS) problems to accelerate the prediction of reservoir pressure and CO2 plume migration. Large amounts of data from physics-based numerical simulators are required to train a model to accurately predict the complex physical behaviors associated with this process. In practice, the available training data are always limited in large-scale 3D problems due to the high computational cost. Therefore, we propose to use a multi-fidelity Fourier neural operator (FNO) to solve large-scale GCS problems with more affordable multi-fidelity training datasets. FNO has a desirable grid-invariant property, which simplifies the transfer learning procedure between datasets with different discretization. Here, we first test the model efficacy on a GCS reservoir model being discretized into 110 k grid cells. The multi-fidelity model can predict with accuracy comparable to a high-fidelity model trained with the same amount of high-fidelity data with 81% less data generation costs. We further test the generalizability of the multi-fidelity model on a same reservoir model with a finer discretization of 1 million grid cells. This case was made more challenging by employing high-fidelity and low-fidelity datasets generated by different geostatistical models and reservoir simulators. We observe that the multi-fidelity FNO model can predict pressure fields with reasonable accuracy even when the high-fidelity data are extremely limited. The findings of this study can help for better understanding of the transferability of multi-fidelity deep learning surrogate models.