Multidimensional, experimental and modeling evaluation of permeability evolution, the Caney Shale Field Lab, OK, USA
Multidimensional, experimental and modeling evaluation of permeability evolution, the Caney Shale Field Lab, OK, USA
Katende, A.↗
SEARCH · Engineering Papers
Results for “Subsurface data models”
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.
Play Fairway analysis of geothermal resources across the State of Hawai‘i: 4. Updates with new groundwater chemistry, subsurface stress analysis, and focused geophysical surveys
This paper is the fourth in a series on a play fairway analysis of geothermal resources across the State of Hawai‘i. Here we describe recent exploration activities that include groundwater sampling in ten locations statewide, as well as geophysical surveys on Lana‘i, across the SW rift of Haleakala Volcano (Maui), and surrounding Mauna Kea Volcano (Hawai‘i Island). We derive a first-order method for computing topographic stress using Green’s functions, finding that topography can induce appreciable crustal stress. We develop a new method for incorporating depth information about resistivity, density, and topographic stresses into our previously published equations to calculate resource probability and confidence. We incorporate newly collected and modeled data into our calculations to update statewide maps of probability and confidence. Lana‘i Island and southeast Mauna Kea are identified as the top targets for exploratory drilling. The east rift of Haleakala Volcano, the southern-most region of Mauna Loa’s SW rift (Hawai‘i Island), and central Kauai are identified as targets for geophysical surveying.
15 GEOTHERMAL ENERGY↗
Deep Learning for Simultaneous Inference of Hydraulic and Transport Properties
Abstract Identification of a heterogeneous conductivity field and reconstruction of a contaminant release history are key aspects of subsurface remediation. These two goals are achieved by combining model predictions with sparse and noisy hydraulic head and concentration measurements. Solution of this inverse problem is notoriously difficult due to, in part, high dimensionality of the parameter space and high computational cost of repeated forward solves. We use a convolutional adversarial autoencoder (CAAE) to parameterize a heterogeneous non‐Gaussian conductivity field via a low‐dimensional latent representation. A three‐dimensional dense convolutional encoder‐decoder (DenseED) network serves as a forward surrogate of the flow and transport model. The CAAE‐DenseED surrogate is fed into the ensemble smoother with multiple data assimilation (ESMDA) algorithm to sample from the Bayesian posterior distribution of the unknown parameters, forming a CAAE‐DenseED‐ESMDA inversion framework. The resulting CAAE‐DenseED‐ESMDA inversion strategy is used to identify a three‐dimensional contaminant source and conductivity field. A comparison of the inversion results from CAAE‐ESMDA with physical flow and transport simulator and from CAAE‐DenseED‐ESMDA shows that the latter yields accurate reconstruction results at the fraction of the computational cost of the former.
Zhou, Zitong↗
Relaxation behavior in low-frequency complex conductivity of sands caused by bacterial growth and biofilm formation by Shewanella oneidensis under a high-salinity condition
Complex electrical conductivity is increasingly used to monitor subsurface processes associated with microbial activities because microbial cells mostly have surface charges and thus electrical double layers. Although highly saline environments are frequently encountered in coastal and marine sediments, there are limited data available on the complex conductivity associated with microbial activities under a high-salinity condition. Therefore, we have developed the spectral responses of complex conductivity of sand associated with bacterial growth and biofilm formation under a highly saline condition of approximately 1% salinity and approximately 2 S/m pore water conductivity with an emphasis on relaxation behavior. A column test is performed, in which the model bacteria Shewanella oneidensis MR-1 are stimulated for cell growth and biofilm formation in a sand pack, whereas the complex conductivity is monitored from 0.01 Hz to 10 kHz. The test results indicate that the real conductivity increases in the early stage due to the microbial metabolites and the increased surface conduction with cell growth but soon begin to decrease because of the reduction of charge passages due to bioclogging. However, the imaginary conductivity significantly increases with time, and clear bell-shaped relaxation behaviors are observed with the peak frequency of 0.1–1 Hz, associated with the double-layer polarization of cells and electrically conductive pili and biofilms. The Cole-Cole relaxation model appears to capture such relaxation behaviors well, and the modeling results indicate gradual increases in normalized chargeability and decreases in relaxation time during bacterial growth and biofilm formation in the highly saline condition. Comparison with previous literature confirms that the high-salinity condition further increases the normalized chargeability, whereas it suppresses the phase shift and thus the imaginary conductivity. Our results suggest that the complex conductivity can effectively capture microbial biomass formation in sands under a highly saline condition.
58 GEOSCIENCES↗
The effect of sub-surface volume scattering on the accuracy of ice-sheet altimeter retracking algorithms
The NASA and ESA retracking algorithms are compared with an algorithm based upon a combined surface and volume (S/V) scattering model. First, the S/V, NASA, and ESA algorithms were used to retrack over 400,000 altimeter return waveforms from the Greenland and Antarctic ice sheets. The surface elevations from the S/V algorithm were compared with the elevations produced by the NASA and ESA algorithms to determine the relative accuracy of these algorithms when subsurface volume-scattering occurs. The results show that the NASA algorithm produced surface elevations within 35 to 50 cm of the S/V algorithm, while the performance of the ESA algorithm was slightly worse. Next, by analyzing several thousand satellite crossover points from the Antarctic data set, we determined the retracking algorithm that produced the most repeatable surface elevations. The elevations derived from the S/V algorithm had the smallest RMS error for the region of the East Antarctic plateau examined here. The ESA algorithm produced erroneous estimates of elevation change when seasonal variations were present; it measured 0.7 to 1.6-m change in elevation over a 6-month period on the East Antarctic plateau where accumulation rates are only 10 cm/year.
Davis, Curt H.↗
Permafrost thermal conditions are sensitive to shifts in snow timing
Abstract Changes in snow precipitation at high latitudes can significantly affect permafrost thermal conditions and thaw depth, potentially exposing more carbon-laden soil to microbial decomposition. A fully coupled process-based surface/subsurface thermal hydrology model with surface energy balance is used to analyze the impact of intra-annual variability in snow on permafrost thermal regime and the active layer thickness. In the four numerical scenarios considered, simulations were forced by the same meteorological data, except the snow precipitation, which was systematically altered to change timing of snowfall. The scenarios represent subtle shifts in snow timing, but the snow onset/melt days, the end of winter snowpack depth, and total annual snow precipitation are unchanged among scenarios. The simulations show that small shifts in the timing of snow accumulation can have significant effects on subsurface thermal conditions leading to active layer deepening and even talik formation when snowfall arrives earlier in the winter. The shifts in snow timing have a stronger impact on wetter regions, especially soil underneath small ponds, as compared to drained regions. This study highlights the importance of understanding potential changes in winter precipitation patterns for reliable projections of active-layer thickness in a changing Arctic climate.
54 ENVIRONMENTAL SCIENCES↗
Planetary Geochemistry Using Active Neutron and Gamma Ray Instrumentation
The Pulsed Neutron Generator-Gamma Ray And Neutron Detector (PNG-GRAND) experiment is an innovative application of the active neutron-gamma ray technology so successfully used in oil field well logging and mineral exploration on Earth, The objective of our active neutron-gamma ray technology program at NASA Goddard Space Flight Center (NASA/GSFC) is to bring the PNG-GRAND instrument to the point where it can be flown on a variety of surface lander or rover missions to the Moon, Mars, Venus, asterOIds, comets and the satellites of the outer planets, Gamma-Ray Spectrometers have been incorporated into numerous orbital planetary science missions and, especially in the case of Mars Odyssey, have contributed detailed maps of the elemental composition over the entire surface of Mars, Neutron detectors have also been placed onboard orbital missions such as the Lunar Reconnaissance Orbiter and Lunar Prospector to measure the hydrogen content of the surface of the moon, The DAN in situ experiment on the Mars Science Laboratory not only includes neutron detectors, but also has its own neutron generator, However, no one has ever combined the three into one instrument PNG-GRAND combines a pulsed neutron generator (PNG) with gamma ray and neutron detectors to produce a landed instrument that can determine subsurface elemental composition without drilling. We are testing PNG-GRAND at a unique outdoor neutron instrumentation test facility recently constructed at NASA/GSFC that consists of a 2 m x 2 m x 1 m granite structure in an empty field, We will present data from the operation of PNG-GRAND in various experimental configurations on a known sample in a geometry that is identical to that which can be achieved on a planetary surface. We will also compare the material composition results inferred from our experiments to both an independent laboratory elemental composition analysis and MCNPX computer modeling results,
Parsons, A.↗
Geothermal play fairway analysis, part 1: Example from the Snake River Plain, Idaho
The Snake River Plain (SRP) volcanic province overlies the track of the Yellowstone hotspot, a thermal anomaly that extends deep into the mantle. Most of the area is underlain by a basaltic volcanic province that overlies a mid-crustal intrusive complex, which in turn provides the long-term heat flux needed to sustain geothermal systems. Previous studies have identified several known geothermal resource areas within the SRP. For the geothermal study presented herein, our goals were to: (1) adapt the methodology of Play Fairway Analysis (PFA) for geothermal exploration to create a formal basis for its application to geothermal systems, (2) assemble relevant data for the SRP from publicly available and private sources, and (3) build a geothermal PFA model for the SRP and identify the most promising plays, using GIS-based software tools that are standard in the petroleum industry. The study focused on identifying three critical resource parameters for exploitable hydrothermal systems in the SRP: heat source, reservoir and recharge permeability, and cap or seal. Data included in the compilation for heat source were heat flow, distribution and ages of volcanic vents, groundwater temperatures, thermal springs and wells, helium isotope anomalies, and reservoir temperatures estimated using geothermometry. Reservoir and recharge permeability was inferred from the analysis of stress orientations and magnitudes, post-Miocene faults, and subsurface structural lineaments based on magnetics and gravity data. Data for cap or seal included the distribution of impermeable lake sediments and clay-seal associated with hydrothermal alteration below the regional aquifer. These data were used to compile Common Risk Segment maps for heat, permeability, and seal, which were combined to create a Composite Common Risk Segment map for all southern Idaho that reflects the risk associated with geothermal resource exploration and identifies favorable resource tracks. Our regional assessment indicated that undiscovered geothermal resources may be located in several areas of the SRP. Two of these areas, the western SRP and Camas Prairie, were selected for more detailed assessment, during which heat, permeability, and seal were evaluated using newly collected field data and smaller grid parameters to refine the location of potential resources. These higher resolution assessments illustrate the flexibility of our approach over a range of scales.
54 ENVIRONMENTAL SCIENCES↗
Hillslope-Channel Transitions and the Role of Water Tracks in a Changing Permafrost Landscape
The Arctic is experiencing rapid climate change, and the effect on hydrologic processes and resulting geomorphic changes to hillslopes and channels is unclear because we lack quantitative models and theory for rapid changes resulting from thawing permafrost. Here, the presence of permafrost modulates water flow and the stability of soil-mantled slopes, implying that there should be a signature of permafrost processes, including warming-driven disturbance, in channel network extent. To inform understanding of hillslope-channel dynamics under changing climates, we examined soil-mantled hillslopes within a ~300 km 2 area of the Seward Peninsula, western Alaska, where discontinuous permafrost is particularly susceptible to thaw and rapid landscape change. In this study, we pair high-resolution topographic and satellite data to multi-annual observations of InSAR-derived surface displacement over a 5-year period to quantify spatial variations in topographic change across an upland landscape. We find that neither the basin slope nor the presence of knickzones controls the magnitude of recent surface displacements within the study basin, as may be expected under conceptual models of temperate hillslope evolution. Rather, the highest displacement magnitudes tended to occur at the broad hillslope-channel transition zone. In this study area, this zone is occupied by water tracks, which are zero-order ecogeomorphic features that concentrate surface and subsurface flow paths. Our results suggest that water tracks, which appear to occupy hillslope positions between saturation and incision thresholds, are vulnerable to warming-induced subsidence and incision. We hypothesize that gullying within water tracks will outpace infilling by hillslope processes, resulting in the growth of the channel network under future warming.
58 GEOSCIENCES↗
Technoeconomic Analysis of Infrastructure Buildout Scenarios
The success of CCUS deployment in the southeast will depend on the optimized pipeline network from CO 2 point sources to geologic sinks for CO 2 storage. An optimal pipeline framework will reduce both environmental impacts and pipeline building cost. However, finding an optimal pipeline/transport infrastructure design is a non-trivial task and requires simultaneous usage of large volumes of data, computational resources, and a state-of-the-art simulator. Such a transport infrastructure model must consider point sources for CO 2 capture and associated volumes, sinks for CO 2 storage, and transportation from source to sink via pipeline networks. These considerations must be addressed simultaneously and systemically in an optimization process, in which a defined objective function (i.e., capital, variable CO 2 capture, transport, and storage cost function) is required to be minimized with the consideration of practical constraints (e.g., CO 2 flow through pipelines is less than the maximum capacity). Although optimization has been widely used in subsurface resources production and CO 2 storage, its application in large scale CCUS infrastructure design is rarely reported in the literature. SimCCS, developed by Los Alamos National Laboratory, is an open-source CCUS pipeline infrastructure design toolset that facilitates the optimization of pipeline infrastructure networks.
99 GENERAL AND MISCELLANEOUS↗
Optimizing Alabama’s CO 2 Storage in Shelby County (Project OASIS) - Deliverable Task 7: Risk Assessment
The “Optimizing Alabama’s CO 2 Storage in Shelby County: Project OASIS” CarbonSAFE Phase II Project seeks to build on regional data sets that demonstrate that the subsurface within Shelby County, Alabama has the potential to store commercial volumes of CO 2 safely, permanently, and economically. The primary target CO 2 Storage Complex is the deep Ketona Dolomite located within a 140 square mile area of the Valley and Ridge Region of Alabama. This deep saline reservoir is beneath a confining system encompassing at least 6,500 ft of shales and other low permeability sediments. Project OASIS drilled a deep stratigraphic test well to confirm the geological properties of the confining system and storage reservoir(s) within the Storage Complex. The geological data was incorporated into numerical models to establish the areal extent of the CO 2 plume and help design the storage site and its monitoring system. The project is managed by the Southern States Energy Board (SSEB), an interstate compact organization consisting of governors and state legislative leaders from sixteen southern states, Puerto Rico, and the U.S. Virgin Islands, as well as an appointee by the President of the United States. The organizational compact provides it with access to state government organizations and legislatures. The Board also maintains an Associate Members program comprised of energy resource companies, utilities, trade groups, academic R&D science and technology experts and energy consultants. Further, the SSEB staff is experienced in managing and coordinating complex energy and environmental programs, from research programs to full-scale design and demonstrations of new and innovative technologies. Project OASIS is a public-private partnership of six entities with multiple principal investigators (PIs). SSEB’s Lead PI and Co-PI are responsible for all aspects of project performance in accordance with the DOE-NETL Cooperative Agreement. SSEB has issued subgrants to Advanced Resources International, Inc., Alabama A&M University, Auburn University, Crescent Resource Innovation, and Oklahoma State University. Advanced Resources International, Inc., issued subgrants to Baker Hughes and Loudon Technical Services for field services
20 FOSSIL-FUELED POWER PLANTS↗
Detrital U-Pb zircon and 40 Ar/ 39 Ar muscovite geochronology from Middle Pennsylvanian strata in the Anadarko Basin, Texas Panhandle, USA
The Late Mississippian-Permian Anadarko Basin formed in Texas and Oklahoma, USA as the result of inversion of Neoproterozic and Cambrian rift structures. Subsidence was driven by flexural loading of the Amarillo-Wichita Uplift, and this uplift may represent the easternmost element of the Ancestral Rocky Mountains system. The northwestern part of this basin has generally been interpreted to have been filled by sediment derived from the Ancestral Front Range Uplift, ~475 km to the northwest during the early stages of basin filling. We test this model using U-Pb detrital zircon and 40 Ar/ 39 Ar detrital muscovite results from three subsurface samples of the Morrow B sandstone in the northwestern part of the Anadarko Basin. We provide a new maximum depositional age of 310.9 ± 4.9 Ma that indicates the age of the Morrow B to be no older than late Atokan to early Desmoinesian Age (Moscovian), ~10 Myr younger than previously interpreted. In contrast to some previous interpretations, we propose that the most likely source for the sediment in the Morrow B is the Amarillo Uplift to the south. Detrital zircon and detrital muscovite data have age peaks at 900–1300 Ma, 1370 Ma and 1600–1800 Ma corresponding to derivation from Grenville, Granite-Rhyolite and Yavapai-Mazatzal basement provinces, respectively. A dominant detrital zircon peak at 1370 Ma suggests that Mesoproterozoic granites in the Amarillo Uplift were exposed by Middle Pennsylvanian time, and detritus eroded from the Amarillo Uplift dominated the lowstand sediment budget of the Texas Panhandle during this time; small volumes of sediment were likely sourced from the Ancestral Front Range. In conclusion, this study presents the first detrital geochronology data from the subsurface Anadarko Basin and the first detrital muscovite data from late Paleozoic southwestern Laurentia. The results presented here highlight the interpretive power of combined detrital zircon and muscovite datasets.
58 GEOSCIENCES↗
3D Geologic Framework Modelling of the Los Alamos National Laboratory Site and Pajarito Plateau: Integrating a realistic 3D fault network and modelling subsurface relationships in a sparsely sampled and complex geologic region
The subsurface geology beneath the Pajarito Plateau is critical to understanding the seismic hazard of the Pajarito Fault System, yet our understanding of this geology is relatively poor. While previous 3D geologic framework models of the area have been created for the purposes of understanding hydrogeologic flow, they are inadequate for the purposes of understanding the Pajarito Fault System. The specific challenges of using oil and gas software for this purpose include: (1) the geologic complexities resulting from volcanism and tectonism; (2) a need for a high level of stratigraphic detail over a large area; (3) a near complete lack of seismic data; and (4) sparse wellbore data. Presented here is a workflow that handles these challenges of adapting commercially available software used by the oil and gas industries to this seismic hazard problem.
58 GEOSCIENCES↗
Reduced ENSO Variability at the LGM Revealed by an Isotope-Enabled Earth System Model
Studying the El Nino Southern Oscillation (ENSO) in the past can help us better understand its dynamics and improve its future projections. However, both paleoclimate reconstructions and model simulations of ENSO strength at the Last Glacial Maximum (LGM; 21 ka B.P.) have led to contradicting results. Here we perform model simulations using the recently developed water isotope-enabled Community Earth System Model (iCESM). For the first time, model-simulated oxygen isotopes are directly compared with those from ENSO reconstructions using the individual foraminifera analysis (IFA). We find that the LGM ENSO is most likely weaker comparing with the preindustrial. The iCESM suggests that total variance of the IFA records may only reflect changes in the annual cycle instead of ENSO variability as previously assumed. Furthermore, the interpretation of subsurface IFA records can be substantially complicated by the habitat depth of thermocline-dwelling foraminifera and their vertical migration with a temporally varying thermocline.
El Niño–Southern Oscillation;isotope-enabled si↗
Field and Lab Testing With Trident Drill to Help Prepare for Future Missions
Field work with instrumentation and technologies that are planned for flight missions is an important means to gain understanding that improves mission performance. This paper reports on results from field workin lunar analog volcanic terrain using an engineering model of the Honeybee Robotics TRIDENT (The Regolith and Ice Drill for Exploration of New Terrains)drill[1]. TRIDENT is a rotary percussive 1-meter class drill that is carried on the PRIME-1and VIPER (Volatiles Investigating Polar Exploration Rover)[2] missions to the moon scheduled to launch in 2024. A similar drilling system was planned for the proposed IcebreakerDiscovery class mission to Mars[3]and the Mars Life Explorer mission recommended by the 2020 Decadal Survey of planetary science[4].The field work objectives were (1) to use data collected by the drill for operational purposes as a probe of subsurface material properties in formations that are analogous to those that may be encountered on planetary surfaces; (2) to correlate subsurface structures deduced by drilling with those inferred from interpretations ofGround Penetrating Radar (GPR)data
C R Stoker↗
Big Data Meets Geothermal Exploration (CRADA Final Report)
As part of the Cyclotron Road program, Zanskar Geothermal & Minerals, Inc. investigated the application of micro-earthquake and ambient noise seismology methods to imaging and characterizing the structural characteristics and hydrothermal flux of subsurface faults. Significant advances in what could be resolved were enabled by two major developments in seismology: 1) the availability of large-n arrays of low-cost seismometers, and 2) the availability of increased computational power and semi-automated data reduction algorithms. In tandem, these advances may improve the signal-to-noise ratio and spatial precision of the data collected and enable higher-resolution characterization of subsurface fracture systems and their spatio-temporal evolution. These tools supported efforts to reduce dry-hole risk and to improve wellfield productivity for geothermal resource development. In particular, two applications of these advances were evaluated: 1) fracture-seismic imaging, which was used to detect ambient emissions from fluid-filled fractures, and 2) reservoir tomography, which used information about travel paths, source locations, and source parameters of micro-earthquakes to identify areas of enhanced permeability. Integration of these methods provided guidance for siting wells and served as prior constraints for reservoir models, informing forecasts of power potential and production and injection strategies aimed at minimizing temperature decline and improving overall resource productivity.
15 GEOTHERMAL ENERGY↗
Insights from a coupled thermo-hydro-mechanical analysis of a layered high-temperature thermal energy storage reservoir
Coupled thermal-hydraulic-mechanical (THM) modeling is applied to investigate the performance of a seasonal high-temperature aquifer thermal energy storage operation based on data and conditions from current site investigations at the Geostorage Forsthaus pilot project in Bern (Switzerland). The model includes subhorizontal sand lenses of various lengths and dips that are embedded in a low permeability clay matrix. Thermal energy storage is simulated by seasonal injection and withdrawal of hot (up to 90 °C) water from a main well, with reservoir pressure regulated by two auxiliary wells at a distance of about 70 m from the main well. The results show how targeted injection into deeper permeable storage formations, along with active deep well pressure control, can effectively minimize geomechanical impact and the potential risk of damaging subsurface storage and sealing formations, or even surface facilities. With such pressure control, the subsurface mechanical responses are dominated by thermal strain and stress, which can be monitored with subsurface fiber optics. The study demonstrates how coupled THM modeling can be applied for the design of a safe and efficient thermal energy storage operation, and how subsurface fiber optic monitoring can be applied for performance confirmation, allowing for more confident operational forecasting.
Rutqvist, Jonny↗
Data for Soil Oxygen Dynamics: A Key Mediator of Tile Drainage Impacts on Coupled Hydrological, Biogeochemical, and Crop Systems
Tile drainage removes excess water and is an essential, widely adopted management practice to enhance crop productivity in the US Midwest and throughout the world. Tile drainage has been shown to significantly change hydrological and biogeochemical cycles by lowering the water table and reducing the residence time of soil water, although examining the complex interactions and feedbacks in an integrated hydrology–biogeochemistry–crop system remains elusive. Oxygen dynamics are critical to unraveling these interactions and have been ignored or oversimplified in existing models. Understanding these impacts is essential, particularly so because tile drainage has been highlighted as an adaptation under projected wetter springs and drier summers in the changing climate in the US Midwest. We used the ecosys model that uniquely incorporates first-principle soil oxygen dynamics and crop oxygen uptake mechanisms to quantify the impacts of tile drainage on hydrological and biogeochemical cycles and crop growth in corn–soybean rotation fields. The model was validated with data from a multi-treatment, multi-year experiment in Washington, IA. The relative root mean square error (rRMSE) for the corn and soybean yield in validation is 5.66 % and 12.57 %, respectively. The Pearson coefficient (r) of the monthly tile flow during the growing season is 0.78. Plant oxygen stress turns out as an emergent property of the equilibrium between the soil oxygen supply and biological demand. The impact of tile drainage on the system is achieved through a series of coupled feedback mechanisms. The model results show that tile drainage reduces the soil water content and enhances soil oxygenation. It additionally increases the subsurface discharge and elevates inorganic nitrogen leaching, with seasonal variations influenced by climate and crop phenology. The improved aerobic condition alleviates crop oxygen stress during wet springs, thereby promoting crop root growth during the early growth stage. The development of greater root density, in turn, mitigates water stress during dry summers, leading to an overall increase in the crop yield by ∼6 %. These functions indicate the potential of tile drainage in bolstering crop resilience to climate change and the use of this modeling tool for large-scale assessments of tile drainage. The model reveals the underlying causal mechanisms that drive the agroecosystem response to drainage on the coupled hydrology, biogeochemistry, and crop system dynamics.
Modeling↗