Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “geothermal fields”

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.

337 records · Page 19

Inaugural Molten Salt Technologies Workshop Powering the Future

The inaugural “Molten Salt Technologies – Powering the Future” Workshop marked a significant convergence of minds from diverse industries, each harnessing molten-salt technologies in innovative ways. Participants from sectors such as solar, geothermal, and advanced nuclear energy, as well as those involved in cutting-edge applications like thermal transport and rare-earth metals extraction, gathered to discuss their shared technical challenges and opportunities. Industry leaders in metal extraction and recycling, alongside experts in high-temperature sensor technology and advanced material manufacturing, also brought their unique perspectives to the table. This workshop served as a crucial platform for these varied industries to delve into the engineering intricacies that molten salt technologies entail. Common challenges such as understanding the thermophysical properties of salts, tackling corrosion mechanisms in harsh environments, and enhancing material resilience under extreme conditions were at the forefront of discussions. These technical sessions highlighted the critical need for cross-industry collaboration to address issues like salt life-cycle process engineering, impurity mitigation, and the development of durable, high-performance materials. By bringing together academia, industry, and representatives from national laboratories and the U.S. Department of Energy (DOE), the workshop facilitated a rich exchange of knowledge and experiences. This interaction not only fostered new partnerships but also strengthened the network among existing collaborators, setting the stage for joint solutions to the complex problems faced by all sectors using molten salt technologies. The event underscored the importance of collaborative efforts in overcoming common engineering challenges and advancing the application of molten-salt technologies across various industries. The workshop not only provided an essential forum for networking and idea exchange but also highlighted the collective drive towards innovative solutions that could benefit multiple fields.

Department of Energy↗

The First X-ray Diffraction Patterns of Clay Minerals from Gale Crater

The Mars Science Laboratory (MSL) Rover, Curiosity spent approx 150 sols at Yellowknife Bay (YKB) studying a section of fluvio-lacustrine sedimentary rocks (with potential indications of volcanic influence), informally known as the Yellowknife Bay formation. YKB lies in a distal region of the Peace Vallis alluvial fan, which extends from the northern rim of Gale Crater toward the dune field at the base of Mt Sharp. Sedimentological and stratigraphic observations are consistent with the Yellowknife Bay formation being part of a distal fan deposit, which could be as young as middle Hesperian to even early Amazonian in age (approx 3.5 to 2.5 Ga). The Yellowknife Bay formation hosts a unit of mudstone called the Sheepbed member. Curiosity obtained powdered rock samples from two drill holes in the Sheepbed Member, named John Klein and Cumberland, and delivered them to instruments in Curiosity. Data from CheMin, a combined X-ray diffraction (XRD)/X-ray fluorescence instrument (XRF), has allowed detailed mineralogical analysis of mudstone powders revealing a clay mineral component of approx 20 wt.% in each sample. The clay minerals are important indicators of paleoenvironmental conditions and sensitive recorders of post-depositional alteration processes. The XRD pattern of John Klein reveals a 021 band consistent with a trioctahedral phyllosilicate. A broad peak at approx 10A with a slight inflexion at approx 12A indicates the presence of 2:1 type clay minerals in the John Klein sample. The trioctahedral nature of the clay minerals, breadth of the basal reflection, and presence of a minor component with larger basal spacing suggests that John Klein contains a trioctahedral smectite (probably saponite), whose interlayer is largely collapsed because of the low-humidity conditions. The XRD patterns show no evidence of corrensite (mixed-layer chlorite/smectite) or chlorite, which are typical diagenetic products of trioctahedral smectites when subjected to burial and heating >60degC in the presence of water. Given estimated geothermal gradients on Mars temperatures <60 degC might still be consistent with (but do not require) moderate burial. However, our ability to identify interstratified minerals is greatly limited by the lack of access to traditional treatments methods used in the lab (e.g., ethylene glycol solvation). Our preferred explanation for the origin of trioctahedral smectites in Sheepbed mudstone is in situ production via reaction of olivine, water and Si-bearing amorphous material, an important mudstone component detected by XRD. Elevated levels of magnetite in the Sheepbed and the trioctahedral monomineralic nature of the clay minerals support this model. These observations, combined with previous studies of olivine stability, support the persistence of circum-neutral hydrous conditions for thousands of years at YKB.

Bristow, Thomas↗

Paleoenvironmental Implications of Clay Minerals at Yellowknife Bay, Gale Crater, Mars

The Mars Science Laboratory (MSL) Rover, Curiosity spent approx 150 sols at Yellowknife Bay (YKB) studying a section of fluvio-lacustrine sedimentary rocks (with potential indications of volcanic influence), informally known as the Yellowknife Bay formation. YKB lies in a distal region of the Peace Vallis alluvial fan, which extends from the northern rim of Gale Crater toward the dune field at the base of Mt Sharp. Sedimentological and stratigraphic observations are consistent with the Yellowknife Bay formation being part of a distal fan deposit, which could be as young as middle Hesperian to even early Amazonian in age (approx. 3.5 to 2.5 Ga). The Yellowknife Bay formation hosts a unit of mudstone called the Sheepbed member. Curiosity obtained powdered rock samples from two drill holes in the Sheepbed Member, named John Klein and Cumberland, and delivered them to instruments in Curiosity. Data from CheMin, a combined X-ray diffraction (XRD)/X-ray fluorescence instrument (XRF), has allowed detailed mineralogical analysis of mudstone powders revealing a clay mineral component of approx. 20 wt.% in each sample. The clay minerals are important indicators of paleoenvironmental conditions and sensitive recorders of post-depositional alteration processes. The XRD pattern of John Klein reveals a 02l band consistent with a trioctahedral phyllosilicate. A broad peak at approx. 10A with a slight inflexion at approx. 12A indicates the presence of 2:1 type clay minerals in the John Klein sample. The trioctahedral nature of the clay minerals, breadth of the basal reflection, and presence of a minor component with larger basal spacing suggests that John Klein contains a trioctahedral smectite (probably saponite), whose interlayer is largely collapsed because of the low-humidity conditions. The XRD patterns show no evidence of corrensite (mixed-layer chlorite/smectite) or chlorite, which are typical diagenetic products of trioctahedral smectites when subjected to burial and heating >60 C in the presence of water. Given estimated geothermal gradients on Mars temperatures <60 C might still be consistent with (but do not require) moderate burial. However, our ability to identify interstratified minerals is greatly limited by the lack of access to traditional treatments methods used in the lab (e.g., ethylene glycol solvation). Our preferred explanation for the origin of trioctahedral smectites in Sheepbed mudstone is in situ production via reaction of olivine, water and Si-bearing amorphous material, an important mudstone component detected by XRD. Elevated levels of magnetite in the Sheepbed and the trioctahedral monomineralic nature of the clay minerals support this model. These observations, combined with previous studies of olivine stability, support the persistence of circum-neutral hydrous conditions for thousands of years at YKB.

MSL↗

Modeling supercritical CO 2 flow and mineralization in reactive host rocks with PFLOTRAN v7.0

Understanding the flow and reactivity of CO 2 injected into geological reservoirs is important for many subsurface applications including secure geologic carbon storage (GCS), critical mineral extraction, enhanced geothermal systems (EGS), and enhanced oil recovery (EOR). Traditionally, subsurface CO 2 injection for GCS applications has focused on geologic formations with favorable subsurface configurations for CO 2 migration and trapping through non-reactive mechanisms such as structural, solubility, and petrophysical trapping. Recently, CO 2 -reactive rocks such as mafic and ultramafic basalts have been investigated for their potential to react with injected CO 2 in situ to simultaneously dissolve host rock minerals and mineralize CO 2 as carbonates. Engineering rapid CO 2 mineralization in the subsurface is attractive because of the increased density of stored CO 2 , the additional safety factors associated with solidification, and the potential to extract valuable critical minerals. Here we present recent developments in the parallel flow and reactive transport simulator PFLOTRAN to model coupled CO 2 -brine flow and reactive transport for a wide range of injection and production applications involving reactive CO 2 -brine systems. These developments are based on the well established and trusted CO 2 flow capabilities in the STOMP-CO 2 simulator. New capabilities added to PFLOTRAN include new CO 2 -brine equations of state with optional thermal coupling, several new constitutive relationships like capillary pressure smoothing and scanning path hysteresis, a fully implicit well model, and native linkage with PFLOTRAN's well-established reactive transport libraries. A series of benchmarks between PFLOTRAN and STOMP-CO 2 verify the newly developed CO 2 -brine flow capabilities. Demonstrations of coupled CO 2 -brine flow modeling and reactive transport show how CO 2 mineralization can be engineered in reactive host rocks. Finally, an example use case involving copper leaching by CO 2 and critical mineral extraction is presented to showcase the strengths of this new implementation. Several limitations still remain, including limited availability of field data to parameterize models. Future work should constrain the evolution of mineral surface area during mineralization and the temperature and/or pH dependence of geochemical reactions for specific systems of interest.

Critical Minerals↗

Ultra High-Temperature Magnetic Bearing System for s-CO 2 Turbines/Expanders. Phase II Final Report

The goal of this project is to design and construct of a prototype of ultra-high temperature permanent magnet (UHT–PM), which provides the force to pull the spinning shaft towards the target position in the magnetic bearing (MB) clearance circle. This project will develop a PM biased MB system, that is actively cooled to the 550°C maximum operating temperature limit of PMs. The MB test environment will include 700°C, 40krpm and 4,000 psi requirements to advance the supercritical carbon dioxide (sCO 2 ) and other types of Brayton power cycle systems. The proposed novel MBs combine PMs for supplying the bias field and electromagnets for supplying the control field, both in UHT conditions. The choice of material is the UHT Sm-Co PMs (US patent 6,451,132) due to their significantly better thermal stability. The 6-pole homopolar design permits uninterrupted force delivery via a redundant control capability (US patent 7,429,811) even if 3 poles were to fail. The Phase I efforts was to focus on development and demonstration of component technologies and detailed design of full scale components and a demonstrator test rig for use in Phase II. UHT–PMs were produced and used in the MB actuator. The operating environment of the EEC-T550 Sm-Co magnet material could be extended by implementing active cooling. Electromagnetic coils imbedded in the MB provide forces to stabilize the shaft position at the target via sensing and control stages. Shaft position sensors enable to identify the present position of the spinning shaft relative to its desired target position in the MB clearance circle (air gap). Auxiliary bearings will be installed to provide a backup support system in the event that power to the MB fails or there is a fault in a device providing current to the MB coils. For control, reliability and safety reasons the sensors and auxiliary bearings are typically located as close as possible to the MBs, which exposes them to the same environmental conditions. Development of only the actuator for the UHT environment is an inadequate approach to complete integration of the MB system in the turbine/expander. Therefore the UHT sensor and auxiliary bearing were designed and developed along with the actively cooled, UHT-PM actuator development. Although shaft position sensors are commercially available (Kaman) they may cost up to $10k per channel, and a full 5 axis magnetic suspension is best implemented with up to 20 sensors. Once successfully completed, this project will significantly extend the operating environment capability for turbomachinery applications such as nuclear power, concentrated solar thermal, fossil fuel, geothermal, and shipboard propulsion. The further advancement of UHT-PM-MB will also strengthen the US position in PM technology, which is presently heavily dominated by China. The UHT Sm-Co magnet material that will be utilized in the proposed work has been developed exclusively by the project proposer. Many actuators, motors, generators and other MB applications will greatly benefit from the proposed work.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Monitoring Fracture Hydromechanical Evolution in the Lab and Field Using Unsupervised Metric Learning

Fractures evolve in time through thermal‐hydraulic‐mechanical‐chemical (THMC) processes that alter their long‐range hydraulic transport properties and modify subsurface behavior and activities. The location of subsurface fractures makes it necessary to use remote sensing techniques such as passive or active seismic monitoring for fracture characterization. In this paper, we develop a machine learning approach to monitor the evolution of fracture properties using passive seismic sources in a laboratory setting and using active seismic monitoring from the Sanford Underground Research Facility in Lead, South Dakota, at a depth of 1.25 km in amphibolite rock during stimulation of natural fractures as well as during induced fracturing. The unsupervised metric learning technique applies tandem neural networks (twin (Siamese) or triplet) with contrastive loss and adaptive margins to track slowly varying systems for which class or similarity labels are not available. The approach adopts locality‐sensitive hashing to divide time‐ordered contiguous data into an arbitrary number of pseudo‐classes. Contrastive‐loss training with many hash bins generates an evolving latent‐space trajectory. This approach enables unsupervised metric learning for seismic data stacks under the condition of contiguous state sampling and slowly varying fracture properties. The displacement discontinuity theory provides a mechanistic foundation for the fracture‐dependent trajectories that are related to relaxation of fractures with time‐dependent specific stiffness responding to changes in stress or fluid saturation.

02 PETROLEUM↗

UT-GOM2-2 Preliminary Report: Terrebonne Basin Northern Gulf of Mexico, 30 July-28 September 2023

In the summer and fall of 2023, the Gulf of Mexico Deepwater Hydrate Coring Expedition (UT-GOM2-2) drilled, cored, made downhole measurements, and analyzed samples from the seafloor to the base of the gas hydrate stability zone in one location (Site H, WR313) in the Terrebonne basin, deepwater Gulf of Mexico. Analyses of data and samples from the expedition will inform biological, geochemical, and geomechanical models to constrain the role of gas hydrates in the carbon cycle and the potential for gas hydrates as an energy resource. Pressure and conventional cores were collected continuously to a depth of 155.1 meters below the seafloor (mbsf). At deeper depths, cores were taken periodically from hydrate-bearing sands and their bounding muds to a total depth of 861.3 mbsf. 162.6 m of conventional core and 54.8 m of pressure core were obtained. Twelve temperature measurements were made between 27.1 and 144.5 mbsf to determine the geothermal gradient. At the seafloor, more than 4 m of sandy silt of unknown origin was encountered. Beneath this sand, to a depth of ~200 mbsf, the section was composed of interbedded mud and biogenic carbonate ooze. The biogenic ooze correlated to low density and high porosity intervals observed in the previously acquired logging while drilling (LWD) data and as measured. Calcareous nannofossil biostratigraphy constrains the entire record to the Pleistocene (< 0.91 million years) with a pronounced increase in sedimentation rate with depth. Beneath 200 mbsf, the section was predominantly composed of mud with two thicker, hydrate-bearing coarse-grained intervals, which are commonly known as the Blue and Orange sands. The dissolved gas concentration was quantified from pressure cores. In the shallow section, dissolved methane concentration increased below the sulfate-methane transition zone (SMTZ) and reaches saturation (the limit of solubility for methane) at 147 mbsf. Gas expansion was very common in conventional and depressurized pressure (conventionalized) cores below the SMTZ. At deeper depths, the methane concentration within muds bounding the Blue and Orange reservoirs was generally found to be less than saturation. The dissolved and hydrate gas composition is consistent with a microbial source, containing greater than 99.99% methane and only trace concentrations of ethane, propane, and butane. The methane to ethane ratio (C 1 /C 2 ) and the methane to ethane plus propane (C 1 /(C 2 +C 3 )) decrease with depth down to at least 678 mbsf, mainly driven by the increase in ethane with depth. It is unclear if this trend continues through the Orange sand interval. The δ 13 C isotopic signature of methane ranges between -69.9 and -78.5 ‰ Vienna Pee Dee Belemnite (VPDB). Pressure core recovery of all sandy intervals was poor. However, pressure core logs of the Orange sand show intervals of low density and high velocity, which are indicative of high hydrate saturation. One core from within the Orange sand was composed of interbedded graded sandy silt and mud. The sandy silts from this core are composed of mainly quartz and feldspar with some lithics. Most of the recovered pressure core samples are maintained at near in-situ pressure and temperature (within the hydrate stability field) at the University of Texas Pressure Core Center awaiting analysis. In the shallow section, samples will be used to determine the flux of organic carbon through the basin system, find the rate at which that carbon was consumed, and understand the microbial population responsible for these processes. In the deeper section, samples from in and around the hydrate reservoirs will be used to determine the petrophysical properties of the reservoir and bounding seals in these systems.

03 NATURAL GAS↗

Utah FORGE: Phase 3 Magnetotelluric (MT) Data

New high-quality tensor MT data at 122 sites, including the vertical magnetic field and utilizing ultra-remote referencing, have been acquired over the Utah FORGE project area. The results will be used to delineate the densities of faults and fractures in crystalline basement rocks, to define the heat sources, and to derive baseline 3D resistivity structure for later MT monitoring of temporal changes in resistivity structure following well stimulation in the EGS reservoir. There are three files here related to Utah FORGE magnetotelluric (MT) data acquisition and processing. The FORGE MT EDIs zip file contains the observed MT responses in industry-standard EDI format. For each site, there is an EDI response file that utilized a local independent reference for noise cancellation, and a file that utilized a distant reference cancelling noise associated with the DC transmission line of the Delta IPP passing down the west side of Milford Valley. These two site files could be merged as appropriate. The FORGE Model Cell Center file contains the model volume of the 3D Forge MT inversion for characterizing the resistivity structure in the project area. It was derived using finite element inversion methodology described in Wannamaker et al, in the attached FORGE Phase3 Geophysics paper, from the MT observation EDI files. It is ASCII format (.dat) and entries are defined at top of the file in a simple x-y-z-Rho listing in UTM coordinates. The element layers drape the topography so the Rho value layers are not purely horizontal slices. This greatly simplifies the listing.

15 GEOTHERMAL ENERGY↗

Reservoir Thermal Energy Storage Benchmarking (Rev. 3)

A benchmarking analysis of RTES research funded by GTO through the Beyond Batteries projects was conducted against the ESGC to see where they fit within the identified ESGC Use Cases. The projects were found to advance knowledge in multiple ESGC use cases, either directly or in some cases, indirectly as enabling technologies. This analysis is helpful to understand where RTES and associated research fits into the larger discussion around energy storage technologies. Also, a retrospective analysis of the Beyond Batteries projects was conducted to evaluate what the projects learned and how the results can be applied to advance the value of RTES. Major results of each of the studies are summarized in Table 2. Additionally, a comparative metrics analysis for RTES was completed to understand where RTES lies within the energy storage industry. Metrics for evaluation of RTES and its comparison to other storage technologies were selected and ranges of their values compiled. The selected metrics – LCOE (levelized cost of energy), capital costs, roundtrip efficiency, energy storage capacity, and storage time – were chosen based on data availability and have a particularly strong influence on the potential deployment of a storage technology. Charts which compare the metrics are presented in section 4.3 and show ranges for each of the 10 selected technologies. However, due to a lack of domestic operational facilities, values for RTES and for portions of the remaining technologies are based on theoretical modeling and studies of best-case scenarios. LCOE estimates for RTES fall within the lower reaches of Figure 15, but nevertheless amount to 2 – 5 times the ESGC Roadmap goal for LCOE, for example in the Facilitating and Evolving Grid Use Case. Capital costs for RTES sit on the higher end (Figure 16) but are expected to decrease as new projects are developed and the technology is refined. The theoretical roundtrip efficiency reported for RTES varies from mid to high percentages (Figure 17) with efficiencies upwards of 93% in modeled scenarios in the Portland Basin (Bershaw et al.,2020). RTES is also expected to have the largest energy storage capacities and longest storage times, likely matched only by lower efficiency hydrogen storage. To better assess the role that RTES could play in energy storage we examined it’s potential in the U.S. The potential depends on many factors. Recently, many researchers have started looking at deep sedimentary basins, depleted oil and gas fields, and basalt formations as potential targets for RTES development. The United States Geological Survey (USGS) has analyzed various cities and shown substantial RTES potential in the cooling sector (Pepin et al., 2021). By modeling RTES in low-quality groundwater (e.g., brackish), it is shown to be favorable across the U.S. with particular suitability in the Illinois Basin, Coastal Plains, and Basin and Range regions. Seasonal RTES operations have also been modeled in the Portland Basin by those at the USGS and Portland State University to simulate an RTES system supplying heating loads needed for the Oregon Health and Science University. Simulations suggest that high conductive heat loss in the initial years exists but tends to decrease with increasing time and development of the resource due to self-insulating nature of the basalts (Burns et al., 2020). Other national laboratory efforts are taking a close look at many of the technical issues involved with RTES (McLing et al., 2019, McLing et al., 2022). These include difficulties in understanding geochemical, hydrogeological, mechanical, and microbiological changes at such elevated temperatures and operational scenarios. Major gaps in research are identified and suggested for future work. With this increased focus to understand how to make RTES successful in the U.S., this technology could be a potential solution to many of the nation’s energy storage problems. For the energy independence of this country, the DOE should prioritize de-risking this technology by making future investments in pilot-scale demonstrations to attract potential investors.

15 GEOTHERMAL ENERGY↗

Expedition UT-GOM2-2 Site H

Pressure and conventional cores were collected at Site H of the Walker Ridge Protracted Area Block 313 in the Terrebonne Basin, deepwater Gulf of America (Gulf of Mexico) during the University of Texas (UT) Deepwater Hydrate Coring Expedition (UT-GOM2-2). Pressure and conventional cores were collected continuously to a depth of 155.1 meters below the seafloor (mbsf). At deeper depths, cores were taken periodically from hydrate-bearing sands and their bounding muds to a total depth of 861.3 mbsf. 162.6 m of conventional core and 54.8 m of pressure core were recovered. Twelve temperature measurements were made between 27.1 and 144.5 mbsf to determine the geothermal gradient. At the seafloor, more than 4 m of sandy silt of unknown origin was encountered. Beneath this sand, to a depth of about 200 mbsf, the section was composed of interbedded mud and biogenic carbonate ooze. The ooze correlated to low density and high porosity intervals observed in the previously acquired logging while drilling (LWD) data and as measured. These ooze intervals also correspond to lighter sediment color, increased Ca content based on X-ray florescence (XRF) core scanning, and increased calcareous nannofossil abundance. Calcareous nannofossil biostratigraphy constrains the entire record to the Pleistocene (< 0.91 million years), with a pronounced increase in sedimentation rate with depth. Below 200 mbsf, the section was predominantly composed of mud with two thicker, hydrate-bearing coarse-grained intervals, which are commonly known as the Blue and Orange sands. The dissolved gas concentration was quantified from pressure cores. In the shallow section, dissolved methane concentration increased below the sulfate-methane transition zone (SMTZ) and reaches saturation (the limit of solubility for methane) at 147 mbsf. Gas expansion was very common in conventional and depressurized pressure (conventionalized) cores below the SMTZ. At deeper depths, the methane concentration within muds bounding the Blue and Orange reservoirs was generally found to be less than saturation. The dissolved and hydrate gas composition is consistent with a microbial source, containing greater than 99.99% methane and only trace concentrations of ethane, propane, and butane. The methane to ethane ratio (C1 /C2 ) and the methane to ethane plus propane (C1 /(C2 +C3 )) decrease with depth down to at least 678 mbsf, mainly driven by the increase in ethane with depth. It is unclear if this trend continues through the Orange sand interval. The δ13C isotopic signature of methane ranges between -69.9 and -78.5 ‰ relative to the Vienna Pee Dee Belemnite (VPDB) standard. Pressure core recovery in sandy intervals was poor. However, pressure core logs of the Orange sand show intervals of low density and high velocity, which are indicative of high hydrate saturation. One pressure core was degassed and the average hydrate saturation in the core was determined to be 24%. One core from within the Orange sand was composed of interbedded graded sandy silt and mud. The sandy silts from this core are composed of mainly quartz and feldspar with some lithics. Most of the recovered pressure core samples are maintained at near in-situ pressure and temperature (within the hydrate stability field) at the University of Texas Pressure Core Center awaiting analysis. In the shallow section, samples will be used to determine the flux of organic carbon through the basin system, find the rate at which that carbon was consumed, and understand the microbial population responsible for these processes. In the deeper section, samples from in and around the hydrate reservoirs will be used to determine the petrophysical properties of the reservoir and bounding seals in these systems.

03 NATURAL GAS↗

Verification of Coupled Hydraulic Fracturing Simulators Using Laboratory-Scale Experiments

In this work, we aim to verify the predictions of the numerical simulators, which are used for designing field-scale hydraulic stimulation experiments. Although a strong theoretical understanding of this process has been gained over the past few decades, numerical predictions of fracture propagation in low-permeability rocks still remains a challenge. Against this background, we performed controlled laboratory-scale hydraulic fracturing experiments in granite samples, which not only provides high-quality experimental data but also a well-characterized experimental set-up. Using the experimental pressure responses and the final fracture sizes as benchmark, we compared the numerical predictions of two coupled hydraulic fracturing simulators—CSMP and GEOS. Both the simulators reproduced the experimental pressure behavior by implementing the physics of Linear Elastic Fracture Mechanics (LEFM) and lubrication theory within a reasonable degree of accuracy. The simulation results indicate that even in the very low-porosity (1–2 %) and low-permeability (10 -18 m 2 - 10 -19 m 2 ) crystalline rocks, which are usually the target of EGS, fluid-loss into the matrix and unsaturated flow impacts the formation breakdown pressure and the post-breakdown pressure trends. Therefore, underestimation of such parameters in numerical modeling can lead to significant underestimation of breakdown pressure. The simulation results also indicate the importance of implementing wellbore solvers for considering the effect of system compressibility and pressure drop due to friction in the injection line. The varying injection rate as a result of decompression at the instant of fracture initiation affects the fracture size, while the entry friction at the connection between the well and the initial notch may cause an increase in the measured breakdown pressure.

15 GEOTHERMAL ENERGY↗

Deep Basalt Aquifers in Orcus Patera, Elysium Basin Mars: Perspectives for Exobiology Exploration

Direct indicators of shorelines, spillways, and terraces allowed to determine the extent of the Elysium Paleolake between the contour-lines 1000 and 500 m below the Martian datum. The Elysium Paleolake is bordered north by Orcus Patera (14N/181W), which lies west of the Tartarus Montes and Tartarus Colles. The Orcus Patera displays an ellipse-shaped collapsed caldera of 360-km long and 100-km wide. Viking topographic data show that the bottom of the caldera is located at 2500 below the Martian datum, and surrounded by a steep-walled ram art which crest is located at about 0 m elevation. Considering the localization of Orcus Patera in the Elysium paleolake, its altimetry, and the magmatic origin of this caldera, we propose the existence of a paleolake in Orcus Patera generated (a) by juvenile water from magma during the Noachian period, and (b) by intermittent influx of the Elysium Basin from Hesperian to Amazonian. Results are encouraging to consider this site as a potential high-energy source environment for microbial communities. are circumscribed by a 50-km wide lava field mapped as Noachian material. The structure of Orcus Patera represents the record of material erupted from a magmatic reservoir. The caldera is enclosed by steep inner walls (25% measured from topographic data), values which could be in agreement with the presence of a deep magmatic reservoir, as suggested by the typology of Crumpler et.al. The depth of the caldera might be due to the collapse of the magma reservoir, and the release of gases accompanying the magma thermal evolution. Origins of water for the paleolake(s): The water that generated a paleolake in Orcus Patera may have come from two origins: (1) Juvenile water: Plescia and Crips estimated a magma H20 content by weight between 0.5% and 1.5% using for the first value a comparison with terrestrial basalt, and for the second values from a Martian meteorite. The amount of H20 can be estimated by the volume of erupted lava, and the lava content of the caldera. In this study, we adopt a water content of 1%. The total volume of magma that has been contained in the caldera, and the volume of lava contained in the observed lava field is about 110 x 10(exp 6) cubic km, that gives a total volume of 1.10 x 10(exp 6) cubic km of water. The juvenile water expelled by the overpressure within the magma chamber charged with desolved water-vapor may have moved into the crust. The decrease in overburden pressure led to bubble formation. The ascent of these bubbles generated a pressurization of the magma, which was sufficient to fracture the overlaying magma layer, (2) Water from Elysium paleolake. During the Amazonian, the rise of the Elysium paleolake level generated an overspilling that supplied the caldera with water. The southern portion of the crest shows a deep gap 12-km wide at -1500 m elevation, locating the gap between 500 to 1000 in below the assumed water of Elysium paleolake, thus facilitating the influx of Elysium paleolake water into Orcus Patera. Bathymetric calculations give a floor area of 25,500 sq km at -2000 m elevation, and a water volume of 42,000 cubic km, with a lake-level at -1500 m. A substantial amount of water may have percolated through the fractured lava, and part of the volume may have overspilled the northern crest of Orcus Patera to debouch in the Tartarus Montes region. We envision the formation of a subsurface aqueous environment in basaltic rocks at the contact of the two water-source origins, possibly the percolating surface lake water, and more likely the juvenile water. Similarly to terrestrial calderas, Orcus Patera might be surrounded by ring-fractures caused by the collapse of the magma chamber that followed the release of gases. These ring-fractures may have been covered later by sedimentation in the caldera (lacustrine, aeolian, and volcanic), and by mass wasting. The detumescence of the magma in the caldera, and the vesiculation of the juvenile water may have operated simultaneously. Comparatively to terrestrial melts, Martian iron-rich melts are denser. This greater density implies greater effusion rates (eight-times terrestrial values), and larger fissuration widths (two-times terrestrial ones). With increasing vesiculation of magma, the bubbles interact with one-another because there are of similar pressure. They make a magma froth at the contact with the caldera surface, and on the walls of the fractures. In the saturated magma, froth, where the volume ratio of gases-to-liquid is about 4:1, the bubbles form a huge surface area of interconnected spaces. Bubbles near the caldera surface disrupt the magma, and fragmentation takes place, which moves downward through the magma column. On Earth, the bubbles are likely to grow between 1 and 50 mm in diameter due to the difference between the magma surface tension, and the bubble supersaturation pressure. The Martian low-pressure at surface level is likely to accelerate the expansion of the bubbles, and increase their final diameter and number, creating more voids in the magma. The strong magma froth with enclosed juvenile water bubbles interconnected with exsolved gas bubbles constitute a potential geothermal environment for geochemical energy production from basalt and water that does not require excessive temperatures. This process can start at +20C. Similar types of environments have been shown on Earth as potential energy sources for microbial metabolism, and could have provided deep aqueous basaltic niches for possible Martian microorganisms, even geologically recently. During the Amazonian, combination of volcanism and water activity still existed on Mars. Moreover, this type of potential niches open ways for investigation of possible oases of extinct or extant life, not only on paleolakes, and surface hydrothermalism spring areas, but also all large systems of fossae, which combine hydrologic and volcanic activities, and which provide an energy source, and an underground shelter to prevent surface UV bombardment. Additional information contained in the original.

Grin, E. A.↗

Quarterly Research Performance Progress Report (Q8)

As part of Task 1, we have started by testing our modeling capabilities by reproducing isothermal DFIT simulations presented in the literature. Once satisfied with the results we have started by targeting the modeling of the DFITs at conducted at well 58-32. We have a identified a specific test (cycle 4 in zone 2) as the most interesting to be model with GEOS hydraulic fracturing module. Thus, we have first produced results with an isothermal model and adjusted model parameters to get a satisfying match with field pressure data. The, we have added thermal effects and compared the modeling results with and without thermal effects to estimate how thermal effects may influence test interpretation. Models seem to suggest that, for small volumes of fluid, thermal effects are moderate. In Task 2, we have adapted GEOS phase-field formulation to be able to simulate near-wellbore hydraulic fracture nucleation and propagation. We have devised a novel formulation that, compared to other existing ones, incorporates rock strengths. We have submitted a journal publication about our work. We are currently employing this phase-field formulation to model the experiments taking place at U Pitt and help us understand the effect of various parameters. In Task 3, we have built a model of the region surround well 16A and have started modeling stage 3 stimulation because of its simpler planar geometry. After calibrating simulation parameters using known analytical solutions, we have simulated the stage 3 stimulation using our isothermal hydraulic fracturing module, varying the permeability field, the stress conditions including different physics to get a better understanding of the numerical challenges and of the effects of varying these parameters on the simulation results. In Task 4 laboratory experimentation, a set of specialized drilling and injection tools has been customized and constructed to accommodate an inclined well with an orientation of up to 30 degrees relative to material anisotropy or principal stress axes. These inclined samples have also undergone thermal stress and hydraulic fracturing at a temperature of 190 degrees Celsius. Furthermore, both vertical and deviated sampling testing setups enable an extended analysis of post-peak pressure behaviors, facilitating post-test pressure analyses such as the G-function, step rate, and fracture reopening measurements. Thus, the key components of in-situ stress estimation can be extracted and validated through our experiment, providing a solid foundation for validating existing in-situ stress estimation theories or proposing new ones. Simultaneously, we are integrating computer vision techniques with traditional experimental fracture observation methods such as multi-overcore/slicing and water-penetration fracture observation. This combination will prove beneficial in populating the hydraulic fracture patterns database, generated under challenging EGS conditions. This approach aims to deepen our understanding of the complexities in EGS reservoirs and pave the way for future data-driven investigations. Additionally, PITT has also equipped the ELE International compression machine, which is now prepared for conducting indirect tensile and fracture toughness tests. These tests will aid in characterizing how rock fabrics influence the resulting fracture patterns. Additionally, we have completed the required personnel training and gained access to Scanning Electron Microscopy (SEM) and Energy Dispersive Spectroscopy (EDS) for conducting more detailed characterization and analysis of rock fabrics, as well as the examination of thermal and hydraulically induced fracture patterns. Thus, the PITT team has effectively demonstrated the capabilities of our experimental apparatuses in exploring the thermal effects, well deviation angles, material anisotropy, and operational choices (such as circulation rate, injection fluid viscosity, and injection rate) and their impact on pressure responses and fracture trajectories under the Utah FORGE conditions.

15 GEOTHERMAL ENERGY↗