Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Geophysics for exploration”

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.

At least 415 records · Page 23
All sourcesNASA NTRSDOE OSTI
NASA NTRS2015

Revisiting the West Clearwater Lake Impact Structure, Canada

The West and East Clearwater Lake impact structures are two of the most distinctive and recognizable impact structures on Earth. Known regionally as the "Clearwater Lake Complex", these structures are located in northern Quebec, Canada (56 deg 10 N, 74 deg 20 W) approximately 125 km east of Hudson Bay. The currently accepted diameters are 36 km and 26 km for the West and East structures, respectively. Long thought to represent a rare example of a double impact, recent age dating has called this into question with ages of approximately 286 Ma and approximately 460-470 Ma being proposed for the West and East structures, respectively. Relatively little is known about the East Clearwater Lake structure. There is no surface exposure and what information there is comes from geophysics and two drill cores obtained in the 1960s. In contrast, the West Clearwater Lake structure is relatively well preserved with large ring of islands in the approximately 30 km diameter lake. Much of the work done on West Clearwater stems from field investigations carried out in 1977 driven by the Apollo program, with a focus on the impact melt rocks and other impactites, which are well exposed on the ring of islands. To our knowledge, the Clearwater Lake impact structures have not been the focus of detailed impact geology field investigations since the 1977 expedition and the only geological map that exists is from the 1960s and is at the reconnaissance level. Our knowledge of impact cratering processes have increased substantially since this time, as have the analytical techniques available for samples. This provided the motivation for a joint Canadian-US-UK expedition to the West Clearwater Lake impact structure in August and September 2015, under the auspices of the FINESSE (Field Investigations to Enable Solar System Science and Exploration) project, part of NASA's Solar System Exploration Research Virtual Institute (SSERVI). We focus here on the impactites of the West Clearwater Lake impact structure. Other ongoing studies, also presented at this conference, focus on central uplift formation, the impact-generated hydrothermal system, xxxx and using WCIS as an analog test site for crew studies of sampling protocols].

West Clearwater Lake↗
NASA NTRS2019

The Lunar Geophysical Network Mission

Overarching Principles: Must be better than Apollo (coverage, duration, instrument performance); Learn from the Apollo experience. Lunar Geophysical Network (LGN) New Frontiers (NF)-class mission, as part of the NF-5 call. “This mission consists of several identical landers distributed across the lunar surface, each carrying geophysical instrumentation. The primary science objectives are to characterize the Moon’s internal structure, seismic activity, global heat flow budget, bulk composition, & magnetic field.” Global distribution of multiple stations. Each station should contain a seismometer, heat flow probe, electromagnetic sounder, laser retroreflector (lunar nearside). Each station must be long-lived (e.g., approximately10 years)to allow other stations (from other countries?) to be integrated with the anchor nodes to form the International Lunar Network. Why LGN? Planetary Science: Moon represents an end-member in planetary evolution (large small body, small rocky planet); Primary planetary differentiation preserved; Key to understanding terrestrial planet initial differentiation. Lunar Science: Heat flow probes yield crustal heat budget estimates; Combined with EMS (ElectroMagnetic Sounding), the temperature profile of the deep interior can be modeled along with mineralogy; Seismic and LLR (Lunar Laser Ranging) data also yield structure and compositional information of the lunar interior; High fidelity data from LGN would enhance the usefulness of the GRAIL (Gravity Recovery and Interior Laboratory) and SELENE (Selenological and Engineering Explorer) gravity data. Human Exploration: LGN must be established prior to renewed human lunar activity - we do not know the exact locations or causes of the shallow moonquakes (SMQs) - the largest magnitude seismic events recorded by Apollo (1 event per year of magnitude greater than or equal to 5); Establishing surface infrastructure near SMQ epicenters must be avoided.

Neal, C. R.↗
DOE OSTI2023

High-pressure behavior of 3.65 Å phase: Insights from Raman spectroscopy

Abstract The 3.65 Å phase [MgSi(OH)6] is a hydrous phase that is predicted to be stable in a simplified MgO-SiO2-H2O (MSH) ternary system at pressures exceeding 9 GPa. Along cold subduction zones, it is likely to transport water, bound in its crystalline lattice, into the Earth’s interior. The 3.65 Å phase consists of Mg and Si octahedral sites attached to the hydroxyl group that forms a hydrogen bond and is predicted to undergo pressure-induced symmetrization of the hydrogen bond. Therefore, in this study, we investigate the high-pressure behavior of the 3.65 Å phase using Raman spectroscopy. We have conducted five distinct compressions up to ~60 GPa using two different pressure-transmitting media—alcohol mixture and neon. At ambient conditions, we identified vibrational modes using complementary first-principles simulations based on density functional perturbation theory. Upon compression, we note that the first derivative of the vibrational modes in the lattice region stiffens, i.e., b1lattice > 0. In contrast, the hydroxyl region softens, i.e., b1OH > 0. This is indicative of the strengthening of hydrogen bonding upon compression. We noticed a significant broadening of vibrational modes related to hydroxyl groups that are indicative of proton disorder. However, within the maximum pressures explored in this study, we did not find evidence for pressure-induced symmetrization of the hydrogen bonds. We used the pressure derivative of the vibrational modes to determine the ratio of the bulk moduli and their pressure derivative. We note that the smaller bulk moduli of hydrous phases compared to the major mantle phases are compensated by significantly larger pressure derivatives of the bulk moduli for the hydrous phases. This leads to a significant reduction in the elasticity contrast between hydrous and major mantle phases. Consequently, the detection of the degree of mantle hydration is likely to be challenging at greater depths.

Geochemistry & Geophysics↗
NASA NTRS2002

Initial Results on the Meteoritic Component of new Sediment Cores Containing Deposits of the Eltanin Impact Event

The late Pliocene impact of the Eltanin asteroid is the only known asteroid impact in a deep- ocean (-5 km) basin . This was first discovered in 1981 as an Ir anomaly in sediment cores collected by the USNS Eltanin in 1965. In 1995, Polarstern expedition ANT XII/4 made the first geological survey of the suspected impact region. Three sediment cores sampled around the San Martin seamounts (approx. 57.5 S, 91 W) contained well-preserved impact deposits that include disturbed ocean sediments and meteoritic impact ejecta. The latter is composed of shock-melted asteroidal materials and unmelted meteorites. In 2001, the FS Polarstern returned to the impact area during expedition ANT XVIIU5a. At least 16 cores were recovered that contain ejecta deposits. These cores and geophysical data from the expedition can be used to map the effects of the impact over a region of about 80,000 square km. To date we have measured Ir concentrations in sediments from seven of the new cores and preliminary data should be available for a few more by the time of the meeting. Our initial interpretation of these data is that there is a region in the vicinity of the San Martin Seamounts comprising at least 20,000 square km in which the average amount of meteoritic material deposited was more than 1 g per square cm. This alone is enough material to support a 500 m asteroid. Beyond this is a region of about 60,000 square km, mostly to the north and west, where the amount of ejecta probably averages about 0.2 g per square cm. Another 400 km to the east, USNS Eltanin core E10-2 has about 0.05 g per square cm, so we know that ejecta probably occurs across more than a million square km of ocean floor. A key to future exploration of this impact is to find evidence of the ejecta at more sites distant from the seamounts. We currently have almost no data from regions to the west or south of the San Martin seamounts.

Kyte, Frank T.↗
NASA NTRS2011

Cryogenic Behavior of the High Temperature Crystal Oscillator PX-570

Microprocessors, data-acquisition systems, and electronic controllers usually require timing signals for proper and accurate operation. These signals are, in most cases, provided by circuits that utilize crystal oscillators due to availability, cost, ease of operation, and accuracy. Stability of these oscillators, i.e. crystal characteristics, is usually governed, amongst other things, by the ambient temperature. Operation of these devices under extreme temperatures requires, therefore, the implementation of some temperature-compensation mechanism either through the manufacturing process of the oscillator part or in the design of the circuit to maintain stability as well as accuracy. NASA future missions into deep space and planetary exploration necessitate operation of electronic instruments and systems in environments where extreme temperatures along with wide-range thermal swings are countered. Most of the commercial devices are very limited in terms of their specified operational temperature while very few custom-made and military-grade parts have the ability to operate in a slightly wider range of temperature. Thus, it is becomes mandatory to design and develop circuits that are capable of operation efficiently and reliably under the space harsh conditions. This report presents the results obtained on the evaluation of a new (COTS) commercial-off-the-shelf crystal oscillator under extreme temperatures. The device selected for evaluation comprised of a 10 MHz, PX-570-series crystal oscillator. This type of device was recently introduced by Vectron International and is designed as high temperature oscillator [1]. These parts are fabricated using proprietary manufacturing processes designed specifically for high temperature and harsh environment applications [1]. The oscillators have a wide continuous operating temperature range; making them ideal for use in military and aerospace industry, industrial process control, geophysical fields, avionics, and engine control. They exhibit low jitter and phase noise, consume little power, and are suited for high shock and vibration applications. The unique package design of these crystal oscillators offers a small ceramic package footprint, as well as providing both through-hole mounting and surface mount options.

Patterson, Richard↗
DOE OSTI2025

Nonlinear Ensemble Filtering with Diffusion Models: Application to the Surface Quasigeostrophic Dynamics

The intersection between classical data assimilation methods and novel machine learning techniques has attracted significant interest in recent years. Here, we explore another promising solution in which diffusion models are used to formulate a robust nonlinear ensemble filter for sequential data assimilation. Unlike standard machine learning methods, the proposed ensemble score filter (EnSF) is completely training free and can efficiently generate a set of analysis ensemble members. Here, in this study, we apply the EnSF to a surface quasigeostrophic model and compare its performance against the popular local ensemble transform Kalman filter (LETKF), which makes Gaussian assumptions in the analysis step. Numerical tests demonstrate that EnSF maintains stable performance in the absence of localization and for a variety of experimental settings. We find that while LETKF maintains optimal performance in the case of linear observations of the entire state and a perfect model, EnSF shows improvements over LETKF when nonlinear observations are assimilated and the system is subject to unexpected model errors. A spectral decomposition of the analysis results in this nonlinear observation regime shows that the largest improvements over LETKF occur at large scales (small wavenumbers), where LETKF lacks sufficient ensemble spread. Overall, this initial application of EnSF to a geophysical model of intermediate complexity motivates further development of the algorithm for more realistic problems.

Artificial intelligence↗
DOE OSTI2022

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

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

58 GEOSCIENCES↗
NASA NTRS2011

A Pragmatic Path to Investigating Europa's Habitability

Assessment of Europa's habitability, as an overarching science goal, will progress via a comprehensive investigation of Europa's subsurface ocean, chemical composition, and internal dynamical processes, The National Research Council's Planetary Decadal Survey placed an extremely high priority on Europa science but noted that the budget profile for the Jupiter Europa Orbiter (1EO) mission concept is incompatible with NASA's projected planetary science budget Thus, NASA enlisted a small Europa Science Definition Team (ESDT) to consider more pragmatic Europa mission options, In its preliminary findings (May, 2011), the ESDT embraces a science scope and instrument complement comparable to the science "floor" for JEO, but with a radically different mission implementation. The ESDT is studying a two-element mission architecture, in which two relatively low-cost spacecraft would fulfill the Europa science objectives, An envisioned Europa orbital element would carry only a very small geophysics payload, addressing those investigations that are best carried out from Europa orbit An envisioned separate multiple Europa flyby element (in orbit about Jupiter) would emphasize remote sensing, This mission architecture would provide for a subset of radiation-shielded instruments (all relatively low mass, power, and data rate) to be delivered into Europa orbit by a modest spacecraft, saving on propellant and other spacecraft resources, More resource-intensive remote sensing instruments would achieve their science objectives through a conservative multiple-flyby approach, that is better situated to handle larger masses and higher data volumes, and which aims to limit radiation exposure, Separation of the payload into two spacecraft elements, phased in time, would permit costs to be spread more uniformly over mUltiple years, avoiding an excessively high peak in the funding profile, Implementation of each spacecraft would be greatly simplified compared to previous Europa mission concepts, minimizing new development while achieving the key Europa science objectives. We will report on the status of this evolving concept, and will solicit community feedback, as we pursue an innovative and low-cost ways to explore Europa and investigate its habitability.

Pappalardo↗
NASA NTRS2021

Determining Planetary Boundary Layer Height by Micro-pulse Lidar with Validation by UAV Measurements

Planetary boundary layer height (PBLH) is often used to characterize the structure of the lower atmosphere. Aerosol lidar, a ground-based remote sensing method, provides the vertical distribution of aerosol at a high temporal resolution observation data, from which, the PBL structure and the position of the PBL top can be comprehensively studied. PBLH determination with lidar data depends primarily on the characteristic turbulent motions in the atmosphere and the geophysical location. However, lidar determination of PBLH over densely populated subtropical locations has rarely been discussed; thus, developing retrieval techniques suitable to these areas is necessary. In this study, four PBLH determination methods (Gradient, δ–threshold, Haar wavelet transform, and hybrid image processing) are applied to estimate the PBLH from lidar observations over an urban area in East Asia, and one—the Gradient method—relied on potential temperature measurements from an unmanned aerial vehicle (UAV) flights to validate our results. Our results indicate that a combination of the gradient method and δ-threshold method can provide better results, in terms of diurnal pattern, than using either method individually. Furthermore, the Haar wavelet and the Hybrid image processing can detect the PBL development comparably well, but both methods are dependent on their initial conditions and optimized algorithm settings. In addition, the accompanying UAV observations are conclusively shown to have a high degree of efficacy for validating the lidar data. This research highlights that a combination of PBLH determination methods can better describe the PBLH evolution throughout a day in some cases, while in others less common determination methods are proving useful, and a suite of retrieval methods should still be explored for precisely mapping the PBL in densely populated subtropical areas.

Yueh-Chen Wang↗
DOE OSTI2021

Data Arrays for Microearthquake (MEQ) Monitoring using Deep Learning for the Newberry EGS Sites

The 'Machine Learning Approaches to Predicting Induced Seismicity and Imaging Geothermal Reservoir Properties' project looks to apply machine learning (ML) methods to Microearthquake (MEQ) data for imaging geothermal reservoir properties and forecasting seismic events, in order to advance geothermal exploration and safe geothermal energy production. As part of the project, this submission provides data arrays for 149 microearthquakes between the year 2012 and 2013 at the Newberry EGS Site for use with the Deep Learning Algorithm that has been developed. The data provided includes raw waveform data, location data, normalized waveform data, and processed waveform data. Penn State Geothermal Team has shared the following files from the project: - 149 microearthquakes (MEQs) between 2012 and 2013 at Newberry EGS sites, 'Normalized Waveform Inputs.npz' are normalized waveforms. - labels of 149 MEQs: Processed Waveform Inputs.npz - location labels of 149 MEQs: Location Data.npz Note: .npz is the python file format by NumPy that provides storage of array data.

15 GEOTHERMAL ENERGY↗
NASA NTRS2018

Development of a Mars Lidar (MARLI) for Measuring Wind and Aerosol Profiles from Orbit

Our understanding of the Mars atmosphere and the coupled atmospheric processes that drive its seasonal cycles is limited by a lack of observation data, particularly measurements that capture diurnal and seasonal variations on a global scale. As outlined in the 2011 Planetary Science Decadal Survey and the recent Mars Exploration Program Analysis Group(MEPAG) Goals Document, near-polar-orbital measurements of height-resolved aerosol backscatter and wind profiles area high-priority for the scientific community and would be valuable science products as part of a next-generation orbital science package. To address these needs, we have designed and tested a breadboard version of a direct detection atmospheric wind lidar for Mars orbit. It uses a single-frequency, seeded Nd:YAG laser ring oscillator operating at 1064nm (4 kHz repetition rate), with a 30-ns pulse duration amplified to 4 mJ pulse energy. The receiver uses a Fabry-Perotetalon as part of a dual-edge optical discrimination technique to isolate the Doppler-induced frequency shift of the back scattered photons. To detect weak aerosol backscatter profiles, the instrument uses a 4x4 photon-counting HgCdTeAPD detector with a 7 MHz bandwidth and < 0.4 fW/Hz(exp 1/2) noise equivalent power. With the MARLI lidar breadboard instrument, we were able to measure Doppler shifts continuously between 1 and 30 m/s by using a rotating chopper wheel to impart a Doppler shift to incident laser pulses. We then coupled the transmitter and receiver systems to a laser ranging telescope at the Goddard Geophysical and Astronomical Observatory (GGAO) to measure backscatter and Doppler wind profiles in the atmosphere from the ground. We measured a 5.3 ± 0.8 m/s wind speed from clouds in the planetary boundary layer at a range of 4 to 6 km. This measurement was confirmed with a range-over-time measurement to the same clouds as well as compared to EMC meteorological models. Here we describe the lidar approach and the breadboard instrument, and report some early results from ongoing field experiments.

Cremons, Daniel R.↗
DOE OSTI2024

Tracking Seismic Velocity Perturbations at Ridgecrest Using Ballistic Correlation Functions

Abstract We present results based on data of a dense nodal array composed of 147 stations, deployed in 2022 near the epicenter of the 2019 Mw 7.1 Ridgecrest earthquake to investigate characteristics of the seismic wavefields. Through array analyses, we identified two primary components. First, we observed far-field P waves dominating the 0.5–1.2 Hz frequency range, which are likely primarily generated by wind-driven oceanic swell activity. Second, we detected near-field body waves resulting from anthropogenic activities in the frequency range 2–8 Hz. We examined noise correlation functions derived from data of the dense deployment and regional stations to explore fault-zone seismic velocity changes using ballistic arrivals, with a focus on velocity perturbation shortly before and after the Ridgecrest earthquake sequence. Our findings exhibit distinct behavior compared to results obtained through standard coda-wave interferometry. Particularly, we observed a decrease in P-wave travel time on certain station pairs prior to the 2019 earthquake sequence. Supported by detailed investigation of the local seismic wavefields, we interpret the decreasing P-wave travel time as likely caused by a velocity increase away from the fault, possibly related to fluid migration. However, additional information is necessary to verify this hypothesis.

Geochemistry & Geophysics↗
NASA NTRS2021

The Climates of Earth's Next Supercontinent: Effects of Tectonics, Rotation Rate, and Insolation

We explore two possible Earth climate scenarios, 200 and 250 million years into the future, using projections of the evolution of plate tectonics, solar luminosity, and rotation rate. In one scenario, a supercontinent forms at low latitudes, whereas in the other it forms at high northern latitudes with an Antarctic subcontinent remaining at the south pole. The climates between these two end points are quite stark, with differences in mean surface temperatures approaching several degrees. The main factor in these differences is related to the topographic height of the high latitude supercontinents where higher elevations promote snowfall and subsequent higher planetary albedos. These results demonstrate the need to consider multiple boundary conditions when simulating Earth-like exoplanetary climates.

insolation↗
NASA NTRS1987

Investigation of the enhanced spatial density of submicron lunar ejecta between L values 1.2 and 3.0 in the earth's magnetosphere: Theory

Initial results from the measurement conducted by the dust particle experiment on the lunar orbiting satellite Lunar Explorer 35 (LE 35) were reported with the data interpreted as indicating that the moon is a significant source of micrometeroids. Primary sporadic and stream meteoroids impacting the surface of the moon at hypervelocity was proposed as the source of micron and submicron particles that leave the lunar craters with velocities sufficient to escape the moon's gravitational sphere of influence. No enhanced flux of lunar ejecta with masses greater than a nanogram was detected by LE 35 or the Lunar Orbiters. Hypervelocity meteoroid simulation experiments concentrating on ejecta production combined with extensive analyses of the orbital dynamics of micron and submicron lunar ejecta in selenocentric, cislunar, and geocentric space have shown that a pulse of these lunar ejecta, with a time correlation relative to the position of the moon relative to the earth, intercepts the earth's magnetopause surface (EMPs). As shown, a strong reason exists for expecting a significant enhancement of submicron dust particles in the region of the magnetosphere between L values of 1.2 and 3.0. This is the basis for the proposal of a series of experiments to investigate the enhancement or even trapping of submicron lunar ejecta in this region. The subsequent interaction of this mass with the upper-lower atmosphere of the earth and possible geophysical effects can then be studied.

Alexander, W. M.↗
NASA NTRS2011

Imaging the Moon's Core with Seismology

Constraining the structure of the lunar core is necessary to improve our understanding of the present-day thermal structure of the interior and the history of a lunar dynamo, as well as the origin and thermal and compositional evolution of the Moon. We analyze Apollo deep moonquake seismograms using terrestrial array processing methods to search for the presence of reflected and converted energy from the lunar core. Although moonquake fault parameters are not constrained, we first explore a suite of theoretical focal spheres to verify that fault planes exist that can produce favorable core reflection amplitudes relative to direct up-going energy at the Apollo stations. Beginning with stacks of event seismograms from the known distribution of deep moonquake clusters, we apply a polarization filter to account for the effects of seismic scattering that (a) partitions energy away from expected components of ground motion, and (b) obscures all but the main P- and S-wave arrivals. The filtered traces are then shifted to the predicted arrival time of a core phase (e.g. PcP) and stacked to enhance subtle arrivals associated with the Moon s core. This combination of filtering and array processing is well suited for detecting deep lunar seismic reflections, since we do not expect scattered wave energy from near surface (or deeper) structure recorded at varying epicentral distances and stations from varying moonquakes at varying depths to stack coherently. Our results indicate the presence of a solid inner and fluid outer core, overlain by a partial-melt-containing boundary layer (Table 1). These layers are consistently observed among stacks from four classes of reflections: P-to-P, S-to-P, P-to-S, and S-to-S, and are consistent with current indirect geophysical estimates of core and deep mantle properties, including mass, moment of inertia, lunar laser ranging, and electromagnetic induction. Future refinements are expected following the successful launch of the GRAIL lunar orbiter and SELENE 2 lunar lander missions.

Weber, Renee C.↗
NASA NTRS2003

An Autonomous Instrument Package for Providing 'Pathfinder' Network Measurements on the Surface of Mars

The investigations of the interior and atmosphere of Mars have been identified as high scientific priorities in most planetary exploration strategy document since the time of Viking. Most recently, the National Academy of Sciences has recommended a long-lived Mars network mission as its second highest scientific priority for Mars (after sample return) for the purpose of performing seismological investigations of the interior and studying the activity and composition of the atmosphere. Despite consistent recommendations by advisory groups, Mars network missions (MESUR, Marsnet, InterMarsnet, NetLander/MSR 05, NetLander/Premier 07, NetLander/?? 09) have undergone a strikingly consistent 'Phoenix' cycle of death and rebirth over the past 15 years, and there are still no confirmed plans to address the interior and atmosphere of Mars. The latest attempt is the NetLander mission. The objective of NetLander is to place a network of four landers on Mars to perform detailed measurements of the seismicity and atmospheric pressure, temperature, wind, humidity, and opacity (as well as provide images, subsurface radar sounding profiles, and electric/magnetic field measurements). However, this mission has recently encountered major programmatic difficulties within CNES and NASA. NASA has already cancelled its participation and the mission itself is facing imminent cancellation if CNES cannot solve programmatic issues associated with launching the mission in 2009. In this presentation we will describe an approach that could move us closer to realizing the goals of a Mars network mission and will secure at least one geophysical and meteorological observatory in 2009.

W B Banerdt↗
DOE OSTI2024

Techno-Economic Analysis and Market Potential of Geological Thermal Energy Storage (GeoTES) Charged With Solar Thermal and Heat Pumps

In this project, we developed a techno-economic analysis (TEA) model that can be used to evaluate the viability of a proposed Geological Thermal Energy Storage (GeoTES) design. This MATLAB-based model integrates distinct subsystem models for the reservoir, wells, power cycle, and solar field to capture their distinct characteristics. It applies this approach in simulating GeoTES storage and dispatch operations for durations ranging from hourly to seasonal. Using cases studies based on GeoTES designs provided by industry partners - Premier Resource Management (PRM) and EarthBridge Energy - we validated the TEA model estimations of system performance and costs (such as thermal and electrical power/energy inflow and outflow, capital costs, and levelized costs of energy and storage) for both concentrating solar thermal (CST) and Carnot Battery (CB) pairings with GeoTES (CST-GeoTES and CB-GeoTES). For the CST-GeoTES case, the model was validated against the proposed system designed by PRM. It showed good agreement with PRM's estimations when well and pump costs derived from PRM's estimations were used. When GETEM-based costs were used, there was a slight overprediction due to GETEM's project/site agnostic assumption of these costs. From a sensitivity analysis perspective, the levelized cost of electricity (LCOE) of the CST-GeoTES case was most sensitive to well flow rate and the charging temperature. An optimal design scenario resulted in an LCOE of 0.11 $\$$/kWhe. CST-GeoTES can also provide a source of heat to meet seasonal demands. With 12-hour and 24-hour levelized cost of heat (LCOH) of 0.018 $\$$/kWhth and 0.022 $\$$/kWhth, respectively, CST-GeoTES could be competitive in the California market with an average industrial price of natural gas in California between 0.041-0.047 $\$$/kWhth. The levelized cost of storage (LCOS) for CST-GeoTES depends on the energy storage duration. Although the LCOS is relatively higher for shorter durations (e.g., ~0.50 $\$$/kWhe for 1 hour of storage), it is an order of magnitude lower (0.06 $\$$/kWhe) for longer storage durations and competitive with lithium-ion batteries (beyond 12 hours of storage) and molten-salt thermal energy storage (beyond 32 hours). Energy. Three options were explored and applied to the EarthBridge case study: (1) A Carnot Battery design using R125 working fluid with both hot and cold storage; (2) A Carnot Battery design using R125 working fluid with only hot storage; (3) A Carnot Battery using a commercially available heat pump with carbon dioxide (CO2) working fluid and hot storage only. The CB-GeoTES with cold storage only had a slight (round-trip) efficiency advantage over the system without (43.4% vs. 42.8%). This is because the cold storage is limited by the freezing point of water, so the cold storage is not much colder than the environment. The system using commercially available technologies was the least efficient - partly because different cycles were used in the heat pump (CO2) and heat engine (binary cycle) which leads to some inefficiencies. Using the commercially available design, the levelized cost of energy (LCOS) from the model (0.10 $\$$/kWhe) was higher than that estimated by EarthBridge (0.068 $\$$/kWhe). This is because of the low round-trip (38.7%) efficiency of the commercially available design. Sensitivity analysis reveals that the model is most sensitive to electricity price. Including electricity price in the TEA for CB-GeoTES leads to an increase in LCOS from the base value to 0.25 $\$$/kWhe. To determine storage sites suitable for GeoTES, we gathered and analyzed geological, petrophysical, and geophysical data of oil and gas reservoir and aquifers in California and Texas. We down-selected possible sites based on cut-off values for site characteristics (e.g., reservoir temperature, formation thickness, permeability, porosity, depth, and brine salinity) and preliminary costs. Using this approach, the Carrizo-Wilcox, Yegua-Jackson, and Dockum brackish aquifers in Texas were identified as having the highest suitability. Similarly, in the central California region, the White Wolf, Belridge South Tulare, and Belridge South Reef Ridge were the most suitable. Going further, we assessed the storage potential in the selected sites. To do this we developed distributions of reservoir characteristic data and applied a Monte Carlo-based analysis to account for intrinsic uncertainty in the acquired data. The analysis revealed that the Carrizo-Wilcox aquifer had the highest storage potential with a mean capacity of 554 TWhth (i.e., 63 TWhe). The estimated capacity serves as an upper limit of storage potential given that not all fields in the basin will be developed. We participated in multiple outreach activities including conference presentations, panel session discussions, and the facilitation of a GeoTES workshop at the NREL Golden campus.

15 GEOTHERMAL ENERGY↗
DOE OSTI2022

Probabilistic inversion of circular phase spectra: application to two-station phase-velocity dispersion estimation in western Canada

SUMMARY Periodic directional and temporal measurements are common in seismology, and necessitate specific statistical analyses that are appropriate for circular quantities. In this work, we explore the use of a von Mises distribution as a representation of errors on circular seismological observations. Specifically, we automate the estimation of surface-wave phase-velocity dispersion for the teleseismic two-station method, which generally suffers from a 2π phase ambiguity. The use of Bayesian inverse techniques, which aim to rigorously quantify model parameter uncertainty, have become widespread throughout seismology over the last decade. Here, we apply Bayesian inversion to measurements of surface-wave phase spectra in order to estimate 1-D, path-averaged Earth structure between station pairs. The dispersion curve and associated uncertainties are additional results of the inversion, which can then be used as input for subsequent analyses (e.g. tomography). We demonstrate this technique through application to surface-wave recordings from long-running seismic stations throughout western Canada. Our results for over 10 000 station pairs reveal first-order tectonic features consistent with previous studies, which provides confidence in our approach as well as an initial step towards resolving a full 3-D seismic velocity model for the region. This work also presents a foundation for the inversion of surface-wave phase spectra to estimate 3-D Earth structure directly. Finally, the ideas presented in this work are not limited to the inversion of surface-wave phase spectra, but can also be considered for Bayesian geophysical inversion of any circular quantities.

Gosselin, Jeremy M. (ORCID:0000000203754102)↗