Coupling surface flow and subsurface flow in complex soil structures using mimetic finite differences
Not Available
54 ENVIRONMENTAL SCIENCES↗
SEARCH · Search NASA
Results for “Complex stratigraphy”
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.
The Role of Stratigraphy and Loading History in Generating Complex Compaction Bands in Idealized Field-Scale Settings
The Buckskin Gulch locality in Utah is a landmark example of compaction localization. The outcrop of this locality involves distinct stratigraphic heterogeneity and was exposed to complex loading history. It features multiple sets of deformation bands with different kinematics and orientation. Similar formations were seen in the Valley of Fire, Nevada, and the Orange quarry, France, among other localities. The formation of such complex structures, their propagation mechanisms, and frequency is affected by numerous local and ambient factors whose impacts are not yet fully understood. The simulation of the above-mentioned localities is not feasible because of the limited amount of available information. This work, instead, investigates from a geomechanics standpoint how the interplay among material nonlinearity, outcrop stratigraphy, and loading history interconnects with specific spatiotemporal patterns of compaction band propagation. Our study shows that the system stratigraphy can be responsible for the emergence of coexisting compaction bands with different inclination and kinematics. Specifically, we show that stiffness contrasts induce nonlocal stress changes which may favor the initiation of secondary structures with different compaction localization characteristics. Furthermore, systems of inclined compaction bands induced by burial increase display secondary, noncontemporaneous sets of vertical compaction bands under the effects of postburial shortening. Our results indicate that stages of intermediate burial decrease prior to tectonic shortening can promote the formation of such complex systems. Despite the simplifications involved in our analyses, these findings show how geomechanics computations complement field observations and could provide a mechanics-based validation of site-specific reconstruction hypothesis.
58 GEOSCIENCES↗
Stratigraphy‐Induced Localization of Microseismicity During CO 2 Injection in Illinois Basin
Abstract Subsurface fluid injection stimulates complex hydromechanical interaction, necessitating the integration of geomechanical data across spatial and temporal scales to consider the sophisticated behavior. Induced seismic response is usually associated with the complex reservoir architecture and pre‐existing features that are three‐dimensional, such as local stratigraphy, fractures, faults, and other discontinuities. This study encompasses laboratory characterization of the coupled hydromechanical response of cores extracted from rock formations in Illinois Basin: reservoir ‐ Mt. Simon sandstone, basal seal ‐ Argenta sandstone, and crystalline basement ‐ Precambrian rhyolite. High‐resolution numerical modeling allows considering the three‐dimensional complexity of the Illinois Basin Decatur Project with spatial resolution comparable to one of the active seismic surveys. A detailed reconstruction of the evolving state of stress in formations lacking direct stress measurements is achieved by numerical modeling that integrated laboratory‐derived hydromechanical properties, a porosity‐permeability relationship, active seismic data, and an inverted three‐dimensional porosity distribution. It appears that the microseismic clusters, mainly observed in the crystalline basement during the injection, are linked to zones experiencing more critically stressed conditions prior to injection. These zones have a potential for reactivation during the injection and are attributed to the specific local stratigraphy of the injection site, as well as transfer of triggering perturbations during the injection.
Bondarenko, N. [University of Illinois Urbana‐Cham↗
AEOLIAN INTERDUNE FACIES OF THE NAVAJO SANDSTONE, UTAH
During the Early Jurassic an extensive desert environment covered a large area of the Western Interior of the USA. The Navajo Sandstone erg is thought to be one of the largest dunefields preserved in the rock record, estimated to have covered well over 500,000 km2, though the preserved extent is less due to erosion. The Navajo Sandstone outcrops extensively in eastern and southern Utah, forming spectacular domes and arches, including those in Arches National Park. Navajo Sandstone aeolian dunes consist of stacked, tabular sandstone bodies, displaying steeply dipping crossbedding, and are on average 15 m thick, with high porosity and permeability. In addition to these extensive paleodunes, interdune deposits are also documented, consisting of lensoid sandstone or pale grey limestone, of approximately 1-7 m thick. This research forms part of a broader investigation into the Navajo Sandstone in Utah, focusing on its potential as a CO2 reservoir. The internal stratigraphy of the Navajo Sandstone can be complex, with multiple types of aeolian surfaces present including discontinuous inter-dune contacts, and inclined bedding planes associated with large climbing dunes. Interdune deposits may be associated with more laterally extensive stratigraphic surfaces, highlighting their importance in the internal stratigraphic framework of the Navajo Sandstone erg. Here we interpret facies and depositional processes of aeolian interdune deposits of the Navajo Sandstone throughout Utah, using well log, outcrop, and thin sections. Interdune deposits range in type from short lived, small ponds of predominantly reworked aeolian sediment, to comparatively long-lived oases consisting of lakes with extensive ecosystems. This project aims to increase our understanding of an important potential CO2 reservoir, and enhance understanding of stratigraphic complexity within one of the world's most significant aeolian systems.
Mahon, Elizabeth (ORCID:0000000252692454)↗
Timing of iceberg scours and massive ice-rafting events in the subtropical North Atlantic
Abstract High resolution seafloor mapping shows extraordinary evidence that massive (>300 m thick) icebergs once drifted >5,000 km south along the eastern United States, with >700 iceberg scours now identified south of Cape Hatteras. Here we report on sediment cores collected from several buried scours that show multiple plow marks align with Heinrich Event 3 (H3), ~31,000 years ago. Numerical glacial iceberg simulations indicate that the transport of icebergs to these sites occurs during massive, but short-lived, periods of elevated meltwater discharge. Transport of icebergs to the subtropics, away from deep water formation sites, may explain why H3 was associated with only a modest increase in ice-rafting across the subpolar North Atlantic, and implies a complex relationship between freshwater forcing and climate change. Stratigraphy from subbottom data across the scour marks shows there are additional features that are both older and younger, and may align with other periods of elevated meltwater discharge.
54 ENVIRONMENTAL SCIENCES↗
INTERDUNE DEPOSITS AND STRATIGRAPHIC FRAMEWORK OF THE NAVAJO SANDSTONE IN THE LAKE POWELL AREA, UTAH, USA
The Western Interior of the USA was dominated by desert environments for much of the Jurassic. Of these desert environments, the most extensive by far was the Navajo Sandstone erg, which, in conjunction with its correlative Nugget and Aztec Sandstones, is thought to have covered well over 500,000 km2, and is one of the largest ergs preserved in the rock record. The Navajo Sandstone is also one of a few geologic units which show promise as a potential CO2 sequestration reservoir in Utah. It outcrops along the eastern and southern margin of the state, and subcrops throughout most of the central area. The Navajo Sandstone consists of stacked eolian deposits, composed of well sorted, clean sandstone, deposited in commonly tabular, high angle cross bedded paleodunes, which are on average 15 m thick. Interspersed between these paleodunes are interdune deposits. Interdune deposits consist of sandy to calcareous, predominantly planar bedded, lensoid units, ranging from <1 m up to 7 m thick. In some cases, these interdune deposits occur along more extensive stratigraphic bounding surfaces. Dating within the Navajo Sandstone has historically been very difficult, as a result of its predominantly clean, aeolian sandstone composition. In addition, complexity of internal stratigraphic surfaces coupled with a low structural dip and generally flat topography has resulted in a poor understanding of the internal stratigraphy of the Navajo Sandstone. Exceptional three-dimensional exposures of the Navajo Sandstone in the upper portions of Lake Powell permit us to relate interdune carbonates to a hierarchical bounding surfaces model. Here we use this combined approach to interpret the internal stratigraphy of the Navajo Sandstone in this region. This work forms part of a larger study of the Navajo Sandstone across Utah and its potential as a CO2 reservoir, by developing a more comprehensive stratigraphic framework, and increasing understanding of stratigraphic complexity within one of the most significant eolian systems in the world.
Mahon, Elizabeth↗
The dynamics and stratigraphy of distributary channel networks (Final Technical Report)
Sedimentary reservoirs are frequently composed of the deposits of distributary channel networks. These networks are associated with coarse sediments (and high permeability) within the channel network, and relatively fine sediments (and low permeability) outside the channel network. The shape of the network is complex, and a product of the environmental boundary conditions under which it grows (Figure 1). These systems control heterogeneity at small and large scales within sedimentary basins, and understanding this heterogeneity and the links between scales represents the important scientific frontier addressed by this project.
58 GEOSCIENCES↗
Evaluating the Thickness and Stratigraphy of Ejecta Materials at the Chang’e-4 Landing Site
The landing site of China’s Chang’e-4 (CE-4) probe is located on the mare basalts on the floor of the Von Kármán crater on the lunar far side. The Von Kármán crater is inside the ancient and highly cratered South Pole–Aitken basin, which has experienced complex emplacement sequences of both near and distant ejecta materials. These issues complicated the interpretation of the CE-4 surface in situ measurements of the visible and near-infrared spectrometer and the lunar penetrating radar (LPR) onboard Yutu-2 rover. To evaluate the sources and amounts of all principle foreign materials at the CE-4 landing site, we thoroughly examine the ejecta delivered by crater-forming events that occurred later than the formation of the mare basalts at the CE-4 landing site. We found a total of 16 craters that may have delivered ejecta thicker than 10 cm level superposed on the mare basalts at the CE-4 landing site. Crater Finsen, Von Kármán L, Von Kármán L′, and Maksutov are the top four major foreign material sources, and each of them contributed ejecta thicker than 1 m. Our surveys confirm that the ejecta from Finsen crater are the most dominant foreign materials in the uppermost few meters at the CE-4 landing site and the total impact ejecta deposited upon the mare basalts at the landing site is estimated to be thinner than 30 m. We found that the estimations from Pike’s model are the most consistent with the Yutu-2 LPR observations.
58 GEOSCIENCES↗
Mudstone Baffles and Barriers in Lower Cretaceous Strata at a Proposed CO 2 Storage Hub in Kemper County, Mississippi, United States
The Cretaceous and Tertiary deposits in Mississippi, Alabama, and the adjacent continental shelf constitute a widespread succession of sandstone, mudstone, and carbonate that has proven to be an important target for geologic CO 2 storage in the onshore Gulf of Mexico basin. Integrated analysis of stratigraphy, sedimentology, and reservoir properties based on cores and geophysical well logs indicates that the Paluxy Formation and Washita- Fredericksburg interval present gigatonne-class storage opportunities. The distribution, geometry, and composition of the area is a direct reflection of the original depositional environments, and understanding these factors is essential for understanding the geologic storage potential of the Paluxy Formation and Washita-Fredericksburg interval at the Kemper County energy facility in Mississippi. Geologic characterization of the Mississippi Embayment at the energy facility focused primarily on characterizing the confinement potential of the storage complex. Integration of core analyses and geophysical well logs has yielded a high-resolution stratigraphic analysis of storage reservoirs, baffles, barriers, and seals. Scanning electron microscopy (SEM) coupled with energy dispersive X-ray spectroscopy (EDS), and quantitative X-ray diffraction (XRD) was used to characterize microfabric, pore types, and mineralogy within mudstone of the east-central Mississippi Embayment at the Kemper County energy facility. Mudstone beds in the Paluxy Formation and Washita-Fredericksburg interval have variable thickness and continuity. High water saturation in the Cretaceous mudstone units influences swelling smectite clay and mudrock permeability based on pulse decay analysis is 1–96 nD. These low permeability values indicate that the mudstone units are effective baffles, barriers, and confining intervals that make significant migration of injected CO 2 out of the storage complex unlikely. The numerous baffles and barriers within the target reservoir intervals, moreover, favor the retention of multiple CO 2 in plumes within the abundant stacked sandstone bodies.
03 NATURAL GAS↗
Joint Inversion of Surface Electrical Resistivity Tomography and Seismic Refraction Data between the 200 Areas
Geologic stratigraphy on the Hanford Site influences groundwater and contaminant migration through the aquifer system and the vadose zone. The current geologic framework model (GFM) relies heavily on a sparse distribution of borehole data in some locations to map geologic contacts and hydrologic properties in the subsurface. Non-invasive geophysical methods such as electrical resistivity tomography (ERT), transient electromagnetic surveying, and seismic imaging are being used at Hanford to map subsurface structure in areas with limited well observations. This is to develop and mature the capability of geophysical methods to aid in GFM refinement, to identify regions of subsurface complexity, and for optimal well siting. A joint inversion of co-located seismic refraction and ERT data was carried out for data collected on a ~2.3-km profile between the 200 Areas on the Hanford Site. While ERT and seismic refraction images have sensitivity to overlapping physical properties (porosity, moisture content, lithology), the resolution and physics used to acquire each of these datasets are different and therefore information can be different or mutually complementary. Performing a joint inversion provides a reasonable option for a coherent, coupled interpretation for mutually complementary datasets. Between the 200 Areas, there are few boreholes to interpret the geologic framework model, and these data sets were obtained to provide a first line of evidence toward identifying stratigraphic structure. The seismic refraction and ERT data were independently inverted during fiscal year 2022 and broadly showed a two-layer structure with a trough-like feature that is ~1 km wide and upwards of 150 m deep. The depth of the trough feature was greater in the ERT image compared to the seismic image, which indicated a maximum depth of approximately 110 m. The objective of the joint inversion described in this report was to invert the seismic refraction and ERT data together while constraining the ERT image to be structurally similar to the seismic refraction image. The approach was applied using the geophysical inverse modeling program E4D, which has the capability to invert first-arrival times from seismic refraction data and ERT resistances using a “cross-gradient” constraint. The application of cross-gradient constraints with different weights produces ERT models that show a high degree of similarity within the upper 100 m (above ~120 m elevation). None of the ERT models show an improved structural similarity to the seismic result; therefore, it is recommended that further attempts to jointly interpret these models focus on petrophysics and image resolution. Petrophysical measurements of core samples would improve knowledge of what drives the ERT response in this region and, along with downhole geophysical measurements, could be used to “ground truth” the surface-based geophysical results. Image resolution studies would provide insight into which regions of the inverted images are reliable and which regions are poorly constrained.
58 GEOSCIENCES↗
Timescales of mafic magmatic fractionation documented by paleosecular variation in basaltic drill core, Snake River Plain volcanic province, Idaho, USA
Abstract The timescales over which fractional crystallization and recharge work in mafic volcano-plutonic provinces is subject to great uncertainty. Currently modeled processes are subject to the scale of measurement: monogenetic basaltic fields accumulate over hundreds of thousands of years, consistent with U-Th-Ra isotopic variations that imply 50% crystallization of basic magmas on timescales of 100,000 years or more, whereas crystal diffusion modeling implies phenocryst residence times of ~1–1000 years. Monogenetic basalts of the Snake River Plain in southern Idaho, USA, are up to 2 km thick and postdate passage over the Yellowstone–Snake River Plain hotspot. Detailed lithologic and geophysical logging of core from deep drill holes, along with chemical stratigraphy and high-resolution paleomagnetic inclination measurements, document individual eruptive units, compound lava flows, and basaltic flow groups that accumulated over 1–6 m.y. Hiatuses are commonly marked by loess or fluvial interbeds that vary from ~0.1 m thick to 20 m thick. Radiometric (40Ar-39Ar, detrital zircon U-Pb) and paleomagnetic timescale ages show that the deepest hole (Kimama drill hole, 1912 m total depth) accumulated over ~6 m.y. Cycles of fractional crystallization and recharge are recognized in the chemical stratigraphy as up-section shifts in major and trace elements; these fractionation cycles commonly represent 40%–50% fractionation. Individual fractionation cycles may comprise 20–40 eruptive units (8–17 lava flows) with little to no change in paleomagnetic inclination (0°–1°), whereas adjacent cycles may differ by several degrees from one another or reflect changes in polarity. Rates of paleosecular variation in Holocene lavas and sediments dated using 14C document significant shifts in magnetic inclination over short timescales, ranging from ~0.05° to 2°/decade, with an average of ~0.5°/decade and a minimum rate of 0.05°/decade. This implies that fractionation cycles with ≤1° variation in magnetic inclination formed on timescales of a few decades up to a few centuries (20–200 years). Thus, the lavas collectively represent only a few thousand years of eruptive activity, with major flow groups separated in time by tens to hundreds of thousands of years. We suggest that the rates defined by paleosecular variation capture the timescales of magmatic chamber evolution (fractionation/recharge) in the seismically imaged mid-crustal sill complex; in contrast, we suggest that crystal diffusion modeling captures the residence times in shallow subvolcanic magmatic chambers that underlie individual monogenetic volcanoes.
Geology↗
Geology of the One Earth Energy Site
The One Earth Energy site is one of two sites in the Illinois Storage Corridor (ISC) project. The objectives of the ISC project is to accelerate commercial deployment of carbon capture utilization and storage at two individual sites and receive approvals for Underground Injection Control (UIC) Class VI permits for construction at each site. At the One Earth Energy site, an extensive data collection program was undertaken, which included the drilling of a test well (One Earth Energy #1 [OEE #1]), four 2D seismic lines, and a small 3D seismic survey. The OEE #1 well was drilled in 2022 and acquired extensive core, log, and testing data to characterize the subsurface geology of the site. Coring was focused on the storage interval, the Mt. Simon Sandstone, and the confining interval, the Eau Claire Formation. The core and log data were used to evaluate the sedimentology and sequence stratigraphy, as well as to develop the conceptual geologic model. This report includes the geological summaries of the Mt. Simon Sandstone and the Eau Claire Formation. The extensive analysis of the log data is included in the petrophysical section, showing ranges of porosity, estimated pore size, and the mineral content of selected zones in the well. The separate petrographic technical report entitled “Petrographic and Advanced Geologic Characterization Report on One Earth Energy #1 (API# 1211325373)”, report number DOE-UIUC-0031892-04, details thin section point-counting analysis that includes mineralogical and pore space analysis, including grain size analysis, annotated thin section photomicrographs, scanning electron microscopy (SEM) with energy dispersive X-ray spectroscopy (EDS), and statistics of grain size analysis on Mt. Simon thin sections from OEE #1. The final OEE #1 well data to be included in this geology report is the routine core analysis of both whole core plugs and rotary sidewall core plugs. In addition to the OEE #1 well, four 2D seismic lines and a small 3D survey were acquired as part of the overall subsurface geological characterization. This geology report references the seismic interpretation report, entitled “One Earth Energy Site Seismic Interpretation Task 5.0”, report number DOE-UIUC-0031892-07. This report details the stratigraphic and structural interpretation of the 2D and 3D seismic data acquired at the One Earth Energy site. The 2D seismic data was acquired in 2019 and 2021, and the 3D survey was acquired in 2022. The objectives of the seismic programs were to contribute to the subsurface characterization of the Mt. Simon-Eau Claire Storage Complex by evaluating the continuity of potential storage reservoirs and containment intervals across the project area, and to determine if any geologic features are present that would increase containment risk to the proposed carbon storage project.
09 BIOMASS FUELS↗
Dense Seismic Array Study of a Legacy Underground Nuclear Test at the Nevada National Security Site
The complex postdetonation geologic structures that form after an underground nuclear explosion are difficult to constrain because increased heterogeneity around the damage zone affects seismic waves that propagate through the explosion site. Generally, a vertical rubble-filled structure known as a chimney is formed after an underground nuclear explosion that is composed of debris that falls into the subsurface cavity generated by the explosion. Compared with chimneys that collapse fully, leaving a surface crater, partially collapsed chimneys can have remnant subsurface cavities left in place above collapsed rubble. The 1964 nuclear test HADDOCK, conducted at the Nevada test site (now the Nevada National Security Site), formed a partially collapsed chimney with no surface crater. Understanding the subsurface structure of these features has significant national security applications, such as aiding the study of suspected underground nuclear explosions under a treaty verification. In this study, we investigated the subsurface architecture of the HADDOCK legacy nuclear test using hybrid 2D–3D active source seismic reflection and refraction data. The seismic data were acquired using 275 survey shots from the Seismic Hammer (a 13,000 kg weight drop) and 65 survey shots from a smaller accelerated weight drop, both recorded by ~1000 three-component 5 Hz geophones. First-arrival, P-wave tomographic modeling shows a low-velocity anomaly at ~200 m depth, likely an air-filled cavity caused by partial collapse of the rock column into the temporary postdetonation cavity. A high-velocity anomaly between 20 and 60 m depth represents spall-related compaction of the shallow alluvium. Hints of low velocities are also present near the burial depth (~364 m). The reflection seismic data show a prominent subhorizontal reflector at ~300 m depth, a short-curved reflector at ~200 m, and a high-amplitude reflector at ~50 m depth. Comparisons of the reflection sections to synthetic data and borehole stratigraphy suggest that these features correspond to the alluvium–tuff contact, the partial collapse cavity, and the spalled layer, respectively.
58 GEOSCIENCES↗
Mechanism of Microseismic Generation During Hydraulic Fracturing - With Evidence From HFTS 2 Observations
The objective of this study is to understand how microseismic events are generated during hydraulic fracturing, as well as the role of geomechanical conditions (i.e., stress and mechanical stratigraphy) in this process. In the industry, microseismic event clouds have been generally used as an “outer-boundary” of the “stimulated reservoir volume” (SRV). However, by comparing with other surveillance data (Low frequency Distributed Acoustic Sensing, or LF-DAS strain) in the Hydraulic Fracturing Test Site (HFTS) 2 experiment, we show that this assumption is fundamentally flawed. The HFTS 2 data has three unique observations that have not been commonly observed in other datasets: 1. Due to influence of offset pad depletion, microseismic data shows that hydraulic fractures from the child well can propagate over 3000 feet into the depleted low stress zone. 2. By comparing microseismic and horizontal fiber LF-DAS strain data, we observe that microseismic event cloud does not necessarily reflect the created hydraulic fracture volume. Particularly, the extent of microseismic event clouds near heel stages are much shorter than what is shown with LF-DAS strain data. 3. Microseismic event magnitudes are larger in the depleted regions. Through geomechanical analysis, we demonstrate that the “bedding-plane-slip” model is likely the mechanism for microseismic generation during hydraulic fracturing. Furthermore, this model successfully explains the above field observations from HFTS 2 experiment. We also provide a quantitative relationship connecting the microseismic event magnitude with fracture width increment and layer mechanical property contrast.
58 GEOSCIENCES↗