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.
This report compiles the results of geologic characterization of Task 3 (late-stage reef), Task 4 (active EOR reefs), and Task 5 (new EOR reefs) to demonstrate developed methodologies, geologic variability, and reservoir potential.
The overall objective of this project is to develop a framework for reliable characterization and prediction of the state of stress in the overburden and underburden (including the basement) in CO 2 storage reservoirs using machine learning and integrated geomechanics and geophysical methods. Specifically, we propose to develop workflow encompassing of technologies and/or methods to predict stress and pressure changes due to CO 2 injection in an active tertiary recovery site and their impacts on subtle fault activation, fractures and occurrence of microseismic events and compare responses to field observations. In this project, we anticipate using dataset from the Farnsworth field Unit (FWU) which is operated by Purdure Petroleum. A novel elastic-waveform VSP inversion technique will be used to estimate high-resolution spatial and temporal changes of elastic moduli in CO 2 storage reservoirs, which will be combined with velocity-stress relationship derived from laboratory tests to obtain subsurface pressure and stress. Clustered microseismic data will be jointly inverted for improved focal mechanisms. Least-squares reverse-time migration of microseismic waveform data will be performed to directly image fracture/fault zones. Additionally, a deep neural network machine learning technique with convolutional and recurrent layers will be used for learning the spectro-temporal structures in microseismic waveforms. The results of this geotechnical data analysis will be integrated to develop a high-resolution 3D mechanical earth model extending from the overburden sealing formations to the underburden including the basement. Mechanical properties will be derived through integration of mechanical logs, tests, available results from chemo-mechanical laboratory tests, and elastic inversion of seismic data using a combination of Bayesian and stochastic methods as well as machine learning technique. Failure features (faults/fractures) will be represented and/or modeled based on seismic and core data analysis. A transient hydrodynamic-geomechanical model will be developed through coupling with the calibrated FWU reservoir simulation model. The full physics coupled model will be used to train a reduced order proxy model using machine learning algorithm for estimating stress which will then be used with appropriate constitutive relationships and forward seismological models to simulate pressure changes and induced microseismicity. An advanced optimization framework will be developed to perform a history match to minimize error between field observations and simulated. The history matched proxy model will be verified against the full-physics equivalent. The field observations that will be used in the coupled model calibration process include pressure/stress inverted from VSP, moment magnitude from microseismic analysis, real time downhole pressure measurements, production and injection data. Parameter sensitivity and uncertainty analysis will be performed to characterize the impact of model parameter uncertainty on stress estimates. The proposed project will have significant impact on future field implementation of the proposed technology. Because the project field site is an ongoing CO 2 EOR development, the value of the new technology will be demonstrated in an operational context and evaluated as a viable risk mitigation strategy. Cost/benefit will be evaluated together with the various commercial incentives for CO 2 sequestration available to oil and gas operators. The extensive available dataset and ongoing data acquisition under the SWP Phase III work plan provides flexibility for investigation of multiple approaches and reduces technical risk.
The pore structure of mudrocks is a key characteristic to evaluate flow behavior through these rocks. Although significant advances have been made to resolve pore characteristics, porosity evolution, or pore connectivity, there is still insufficient knowledge linking porosity evolution to flow and transport in mudrocks. To better understand these links, we conducted very small-angle (VSANS) and small-angle neutron scattering (SANS) experiments on 13 sets of mudrocks from global locations, characterized by differences in composition, maturity, and depositional environment. Our results indicate that a homogeneous pore structure reflects a stable and low-energy depositional environment. Organic-lean mudrocks mainly contain pores <10 nm, which restricts flow and is beneficial for, e.g., caprock seal integrity. Increasing maturity as a result of increasing burial depth, and therefore, temperature leads to porosity evolution. In organic-rich mudrocks (ORMs), this happens through the generation of organic matter pores and pore preservation around diagenetically formed carbonates. ORMs mainly contain macropores (>50 nm), favoring production from, e.g., unconventional reservoirs. SANS-derived pore size distributions will help to better understand fluid imbibition and flow properties in mudrocks through comprehensive quantitative characterization of the nano- to micron scale pore network conduits and their relationships to burial diagenesis.
This report presents recommendations for quantitatively assessing and managing leakage risk and effective containment assurance at geologic carbon storage (GCS) sites, amidst uncertainty in site characteristics and performance. These recommended practices are a product of the U.S. Department of Energy’s (DOE) Office of Fossil Energy and Carbon Management’s National Risk Assessment Partnership (NRAP)—a multi-year collaborative research effort of Lawrence Berkeley National Laboratory, Lawrence Livermore National Laboratory, Los Alamos National Laboratory, National Energy Technology Laboratory, and Pacific Northwest National Laboratory. These recommended practices are organized as a set of conceptual workflows that correspond with the stages of the GCS project lifecycle. They consider aspects of site characterization, quantification of potential leakage, assessing the spatial extent of potential environmental impacts, optimizing monitoring design for detection of potential leakage, evaluation of concordance between monitored and forecasted reservoir response, assessing conformance with expected/accepted performance, and determination of the timeline over which risk of impactful leakage can be expected to abate. These conceptual workflows are intended to aid stakeholders in making risk management decisions at GCS sites and to support development of sound, risk-based justifications for characterizing the subsurface, defining area of potential leakage impact, designing site monitoring plans, selecting leakage risk management/mitigation alternatives, and estimating site closure timelines.
The gas hydrate petroleum system (GHPS) approach, which has been used to characterize gas hydrates in nature, uses three distinct components: a methane source, a methane migration pathway, and a reservoir that not only contains gas hydrate, but also acts as a seal to prevent methane loss. Unlike GHPS, a traditional petroleum system (PS) approach further distinguishes between the reservoir, a unit with generally coarser sediment grains, and a separate overlying seal unit with generally finer sediment grains. Adopting this traditional PS distinction in the GHPS approach facilitates assessments of reservoir growth and production potential. The significance of the seal for the formation of a gas hydrate reservoir as well as for efficiency in methane extraction from the reservoir as an energy resource is evident in findings from recent offshore field expeditions, such as India’s second National Gas Hydrate Program expedition (NGHP-02). In regard to gas hydrate-bearing reservoir formations, the NGHP-02 gas chemistry data indicate a primarily microbial methane source. Fine-grained seal sediment in contact with coarser grained reservoir sediment can facilitate that microbial methane production. Logging-while-drilling and sediment core data also indicate that the overlying fine-grained seal sediment is less permeable than the underlying, highly gas hydrate-saturated reservoir sediment. The overlying seal’s capacity to act as a low-permeability boundary is important not only for preventing methane migration out of the reservoir over time, but also for preventing water invasion into the reservoir during methane extraction from the reservoir. Ultimately, the presence of an overlying, fine-grained, low-permeability “seal” influences how gas hydrate initially forms in a coarse-grained reservoir and dictates how efficiently methane can be extracted as an energy resource from the gas hydrate reservoir via depressurization.
Cumulative volcano volume curves are presented as evidence for steady-state behavior at certain volcanoes and to develop a model of steady-state volcanism. A minimum criteria of five eruptions over a year was chosen to characterize a steady-state volcano. The subsequent model features a constant head of magmatic pressure from a reservoir supplied from depth, a sawtooth curve produced by the magma arrivals or discharge from the subvolcanic reservoir, large volume eruptions with long repose periods, and conditions of nonsupply of magma. The behavior of Mts. Etna, Nyamuragira, and Kilauea are described and show continuous levels of plasma output resulting in cumulative volume increases. Further discussion is made of steady-state andesitic and dacitic volcanism, long term patterns of the steady state, and magma storage, and the lack of a sufficient number of steady-state volcanoes in the world is taken as evidence that further data is required for a comprehensive model.
On terrestrial planets and moons of our solar system cores reveal details about a geological structure's formation, content, and history. The strategy for the search for life is focused first on finding water which serves as a universal solvent, and identifying the rocks which such solvent act upon to release the constituent salts, minerals, ferrites, and organic compounds and chemicals necessary for life. Dielectric spectroscopy measures the dielectric properties of a medium as a function of frequency. Reflection measurements in the frequency range from 300 kHz to 300 MHz were carried out using RF and microwave network analyzers interrogating SansEC Sensors placed on clean geological core samples. These were conducted to prove the concept feasibility of a new geology instrument useful in the field and laboratory. The results show that unique complex frequency spectra can be acquired for a variety of rock core samples. Using a combination of dielectric spectroscopy and computer simulation techniques the magnitude and phase information of the frequency spectra can be converted to dielectric spectra. These low-frequency dielectric properties of natural rock are unique, easily determined, and useful in characterizing geology. TIPPER is an Electro-Magnetic Passive-Source Geophysical Method for Detecting and Mapping Geothermal Reservoirs and Mineral Resources. This geophysical method uses distant lightning and solar wind activity as its energy source. The most interesting deflections are caused by the funneling of electrons into more electrically conductive areas like mineralized faults, water or geothermal reservoirs. We propose TIPPER to be used with SansEC for determining terrain/ocean chemistry, ocean depth, geomorphology of fracture structures, and other subsurface topography characteristics below the ice crust of Jovian moons. NASA envisions lander concepts for exploration of these extraterrestrial icy surfaces and the oceans beneath. One such concept would use a nuclear powered heated tip for melting through the ice sheath of Europa and inserting a down hole SansEC with TIPPER interface. NASA's Juno space probe already on the way to Jupiter as part of the Exploration New Frontiers Program and the planned Europa mission will conduct detailed reconnaissance of Jupiter's moon Europa and investigate whether the icy moon could harbor conditions suitable for life. It has already been observed that Jovian moons have auroras that may serve as naturally occurring active energy sources for a TIPPER instrument.
Carbon dioxide is important in its role as coupler of the terrestrial biosphere to inorganic chemical processes and as the principal greenhouse gas controlling Earth's surface temperature. The hypothesis that atmospheric CO2 levels have diminished with time, with the resulting cooling effect offsetting an increase in the solar constant, seems firmly established, and it is shown that feedback mechanisms exist which can maintain the terrestrial surface in a relatively narrow temperature range over geological time. Of the factors involved in such CO2 variation, the oceanic reservoir appears the most important. Surface waters are probably in approximate equilibrium with regard to CO2 exchange with the ambient atmosphere in most regions, but data from deep-ocean water sampling indicates that such waters are somewhat undersaturated in the sense that they would tend to absorb CO2 from the atmosphere if brought to the surface without change in composition or temperature. If major impacts into the ocean can result in loss of a substantial portion of the atmospheric CO2 reservoir, then any such future event could imperil the continuation of most higher forms of life on Earth. The most likely candidate for an inverse Nyos global event in previous Earth history is the Cretaceous-Tertiary terminal extinction event. The Cretaceous was characterized by warm, equable temperatures presumably indicative of relatively high CO2 levels and an intense greenhouse heating. Cooling of the oceans in absence of massive transfer of CO2 to the oceanic reservoir in itself would promote a condition of CO2 undersaturation in abyssal waters, and this is made even more extreme by the pattern of ocean water circulation. It is possible to envision a situation in which deep ocean waters were at least occasionally profoundly undersaturated with regard to CO2. Turnover of a major fraction of such an ocean would then remove, on a very short time scale, as much as 90 percent of the atmospheric CO2 inventory.
This Storage Field Development plan presents the Storage Complex characterization results, construction, monitoring, and operational plans, and costs associated with the proposed One Earth Sequestration Carbon Capture and Storage (OES-CCS) site in McLean County, Illinois, near Gibson City. The proposed storage complex, known as the Mt. Simon Storage Complex, comprises the Cambrian Mt. Simon Sandstone reservoir and the primary seal, the Cambrian Eau Claire Formation. The lowermost Underground Source of Drinking Water (USDW) identified for the site is the Ordovician St. Peter Sandstone. Geologic characterization of the Mt. Simon Storage Complex at the OES-CCS site was performed by the Illinois Storage Corridor CarbonSAFE Phase III project, which also prepared and submitted three UIC Class VI applications to construct three injection wells; the permit applications were submitted and are in the federal EPA review process. A characterization well, OEE #1, was drilled to collect site-specific data. These data were analyzed and used to develop the UIC Class VI applications. The OEE #1 well will be converted to an in-zone monitoring (IZM) well for the injection phase. The proposed buildout for the OES-CCS site includes (1) three injection wells (OES #1, OES #2, and OES #3), (2) two IZM wells, (3) two above confining zone (ACZ) monitoring wells, one of which will be used to monitor the lowermost USDW, (4) capture and compression facilities, and (5) transportation facilities, i.e., pipelines. A pre-operational testing program was proposed in the Class VI permit application and will be employed at the site pending approval. Additional pre-injection (baseline), syn-injection, and post-injection monitoring and site care procedures will be followed by OES to ensure that injection activities are protective of human health and the environment. Injection is scheduled to begin in 2025, distributed across the three injection wells in accordance with the permit operating conditions. One Earth Sequestration intends to inject up to 90 million tonnes of CO 2 over a period of approximately 20 years. Injection will begin at approximately 0.5 million tonnes of CO 2 annually and ramp up to a maximum of 4.5 million tonnes annually. Daily injection rates are expected to range from 1,400 to 1,500 tonnes per day initially and reach a maximum of approximately 4,225 tonnes per day, depending on site geology and injectivity at each injection well location, and CO 2 availability. The costs associated with the OES-CCS project include pre-operational costs (e. g. additional seismic data acquisition and well drilling), capture and transportation facility and equipment costs, predicted field operating expenditures (OpEx), and decommissioning and post-injection site care (PISC) costs. The risks associated with project activities, such as site construction, injection operations, and verification of secure storage were evaluated, and mitigation strategies proposed to alleviate those risks.
Multidimensional, experimental and modeling evaluation of permeability evolution, the Caney Shale Field Lab, OK, USA
Human health is affected by indoor air quality. One distinctive aspect of the indoor environment is its very large surface area that acts as a poorly characterized sink and source of gas-phase chemicals. In this work, air-surface interactions of 19 common indoor air contaminants with diverse properties and sources were monitored in a house using fast-response, on-line mass spectrometric and spectroscopic methods. Enhanced-ventilation experiments demonstrate that most of the contaminants reside in the surface reservoirs and not, as expected, in the gas phase. They participate in rapid air-surface partitioning that is much faster than air exchange. Phase distribution calculations are consistent with the observations when assuming simultaneous equilibria between air and large weakly polar and polar absorptive surface reservoirs, with acid-base dissociation in the polar reservoir. Chemical exposure assessments must account for the finding that contaminants that are fully volatile under outdoor air conditions instead behave as semivolatile compounds indoors.
Microbial conversion of coal to methane is a promising technology during transition of coal-based energy from conventional coal mining to natural gas recovery from coal. Significant research advances have been made towards engineering ideal microbial communities and nutrients for bio-stimulation of coal. However, actual field applications require geomechanical and flow behavior characterization of coal during the bioconversion process as well as gas production over the life of the created “biogenic gas reservoirs”. This work presents the results of an experimental investigation to estimate gas/methane production using bioconversion of coal and analyses of variations in bulk modulus, strain and permeability. Bulk modulus of coal demonstrated time-dependent behavior with continued bioconversion of solid coal. The phenomenon was modeled using mass balancing in a closed environment and the logistic equation based biogenic gas production with time. The modeled gas production and changes in modulus of coal showed excellent agreement with the corresponding experimental results. Next, numerical simulation of biogenic conversion of coal in a constant stress in-situ condition, replicating the Huff ‘n Puff method under field conditions, was carried out. The results showed that bioconversion can lead to decrease in effective stress, increase in permeability and decrease in the modulus of coal with time. Repeated recharge of coal with nutrients and recovering the produced gas, that is, sequential cycles of Huff n’ Puff, showed further increase in permeability and decrease in coal strength, potentially leading to coal failure in-situ, further increasing the permeability, thus enhancing the prospect of field application of the technology.
Anaerobic fungi (class Neocallimastigomycetes) thrive as low-abundance members of the herbivore digestive tract. The genomes of anaerobic gut fungi are poorly characterized and have not been extensively mined for the biosynthetic enzymes of natural products such as antibiotics. Here, we investigate the potential of anaerobic gut fungi to synthesize natural products that could regulate membership within the gut microbiome. Here, complementary 'omics' approaches were combined to catalog the natural products of anaerobic gut fungi from four different representative species: Anaeromyces robustus (A. robustus), Caecomyces churrovis (C. churrovis), Neocallimastix californiae (N. californiae), and Piromyces finnis (P. finnis). In total, 146 genes were identified that encode biosynthetic enzymes for diverse types of natural products, including nonribosomal peptide synthetases and polyketide synthases.
Estimating the distance from hydraulic fracture tip to monitor well can be very useful for fracture characterization, well spacing optimization, and preventing parent-child well interference. Heart-shape signal is referred to as the extensional precursor of fracture hit recorded by cross-well strain measurements and can be served as a vital tool to make such estimation. This study incorporates the 3D Displacement Discontinuity Method to understand the impact of fracture geometry and monitor well offset on the heart-shape signal’s characteristics. Results from numerical simulation and analytical solutions reveal a strong linear correlation between the spatial extent of heart shape signal and fracture tip distance. This relationship was further developed to predict tip distance using field data from the Hydraulic Fracture Test Site 2. A reasonable approximation result from field data further validates the methodology. In addition, it is worth noting that the estimation accuracy depends on the ratio between fracture dimension and tip distance. The findings of this study offer a novel approach for real-time monitoring and characterizing hydraulic fracture propagation, which can be further used for well spacing optimization in unconventional and Enhanced Geothermal System reservoir development, as well as cap rock integrity monitoring for carbon sequestration projects.
Abilities to accurately characterize the subsurface in a geothermal setting is key to assess and support production. An important element of geothermal reservoir monitoring is also the ability to investigate fluid transport within fracture network. This report focuses on improving subsurface imaging and monitoring in geothermal settings using full waveform inversion based on the adjoint method and time-lapse imaging. To assess our method, we rely on a dense seismic dataset collected in 2016 at the Brady Hot Springs geothermal site in Nevada for the DOE-funded project Poroelastic Tomography by Adjoint Inverse Modeling of Data from Seismology, Geodesy, and Hydrology. This dataset captures subsurface changes across four stages of geothermal power plant operations, which involve varying rates of fluid injection and extraction. Two velocity models were previously derived from this dataset using different methods: one based on travel times and another on sweep interferometry. Our first step is to refine these models using adjoint tomography, which has been applied successfully at global and regional-scales but is less common at the reservoir-scale. Two approaches are then explored for time-lapse analysis: directly comparing refined tomographic models from different stages or backpropagating waveform differences relative to a baseline tomographic model. The main take away is that both approaches highlight similar reservoir behaviors, but the latter approach is more computationally effective in capturing small-scale changes in subsurface properties. For this work, we leverage the use of Salvus (www.mondaic.com), an end-to-end seismic imaging solution, relying on the spectral element method to compute forward and adjoint simulations, and developed by Mondaic Ltd. It includes integrated workflow management that handles waveform and metadata, launches simulations, computes waveform misfits and adjoint sources, and iterates for model updates by nonlinear optimization.
The proposed paper will present a numerical investigation of the flow characteristics and boundary layer development in the nozzles of high enthalpy shock tunnel facilities used for hypersonic propulsion testing. The computed flow will be validated against existing experimental data. Pitot pressure data obtained at the entrance of the test cabin will be used to validate the numerical simulations. It is necessary to accurately model the facility nozzles in order to characterize the test article flow conditions. Initially the axisymmetric nozzle flow will be computed using a Navier Stokes solver for a range of reservoir conditions. The calculated solutions will be compared and calibrated against available experimental data from the DLR HEG piston-driven shock tunnel and the 16-inch shock tunnel at NASA Ames Research Center. The Reynolds number is assumed to be high enough at the throat that the boundary layer flow is assumed turbulent at this point downstream. The real gas affects will be examined. In high Mach number facilities the boundary layer is thick. Attempts will be made to correlate the boundary layer displacement thickness. The displacement thickness correlation will be used to calibrate the quasi-1D codes NENZF and LSENS in order to provide fast and efficient tools of characterizing the facility nozzles. The calibrated quasi-1D codes will be implemented to study the effects of chemistry and the flow condition variations at the test section due to small variations in the driver gas conditions.
The Cypress Sandstone is the youngest and shallowest unit in the Illinois Basin that was featured in the United States Carbon Utilization and Storage Atlas IV as a target for saline carbon storage with an estimated 0.2 to 2.3 GT of storage potential. Additional research on a residual oil zone (ROZ) developed within the Cypress Sandstone has delineated 27 prospects with approximately 290.8 million m3 (1.8 billion barrels) of oil in place. 21 to 31 million m 3 (144 to 196 million barrels) of oil is estimated to be recoverable using carbon dioxide enhanced oil recovery (CO 2 -EOR). Storage of CO 2 associated with EOR in these ROZ prospects alone, not accounting for associated main pay zones (MPZs), underlying brine formation, or intervals adjacent to or between prospects, is estimated to be up to 10.4 billion tonnes. The Cypress Sandstone is thus well understood to be a CO 2 injection target, both for EOR and associated storage. However, the seal system overlying the Cypress Sandstone is poorly understood. Unlike deeper formations such as the Mt. Simon Sandstone or the St. Peter Sandstone which are either in use as a CO 2 sink or being characterized for prospective storage, respectively, the Cypress is not overlain by hundreds of feet of impermeable shale. Rather, the Cypress is overlain by a lithologically variable interval that is composed generally of shales and limestones with some sandstone in the part of the Basin where the Cypress is deep enough to facilitate CO 2 storage. Also, due to its status as one of the shallowest and most prolific oil reservoirs in the Basin, the seal system overlying the Cypress Sandstone has a relatively high number of legacy well penetrations. The purpose of this report is to characterize the Cypress Sandstone seal system using well logs and available core. Gross thickness, lithology (facies), and mineralogy of seals is described and mapped across the Basin. The column height of CO 2 that can be held is calculated using capillary pressure data from a representative core.
As the earliest-dated solids in our Solar System, calcium-aluminum-rich inclusions (CAIs) provide a record of their formation environment near the young Sun and hold clues to the formation of planetary-scale isotopic reservoirs in the solar protoplanetary disk. Although CAIs from several CV, CK, CM, CO, and ordinary chondrites have been analyzed previously for their Ti isotopic compositions, CAIs from just three CV chondrites have been analyzed for their Cr isotopic compositions, and only a handful of CAIs have been measured for both their Ti and Cr isotopic compositions. We report mass-independent Ti and Cr isotopic anomalies in several CAIs from CV and CK chondrites; this is the first report of the Cr isotopic composition of a CAI from a CK chondrite. With this data, we aim to better constrain the compositional range of CAIs in ε 50 Ti versus ε 54 Cr space, thereby facilitating the isotopic characterization of the material inherited by the solar protoplanetary disk and the role of CAIs in the formation of distinct planetary-scale isotopic reservoirs in our early Solar System. The narrow range in isotopic anomalies in CAIs when compared to other early-formed refractory inclusions such as platy hibonite crystals (PLACs) and spinel-hibonite inclusions (SHIBs) suggests that CAIs record the mixing of these precursor materials and the averaging of their larger isotopic anomalies. The isotopic composition of CAIs is therefore likely the result of a combination of factors, including mixing of material inherited from their formation region, heterogeneous carrier phase distribution, and thermal processing in the disk. The ε 50 Ti and ε 54 Cr isotopic compositions of CAIs are not correlated, further demonstrating that these isotopic anomalies have different carrier phases. The Ti and Cr isotopic compositions of CAIs additionally show that CAIs alone cannot be responsible for the compositional difference between the non-carbonaceous chondritic (NC) and carbonaceous chondritic (CC) isotopic reservoirs but nevertheless do play a role in the formation of these large-scale isotopic reservoirs in the early Solar System.