Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “BOREHOLES”

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 271 records · Page 15

Microseismic Monitoring at the Farnsworth CO 2 -EOR Field

The Farnsworth Unit in northern Texas is a field site for studying geologic carbon storage during enhanced oil recovery (EOR) using CO 2 . Microseismic monitoring is essential for risk assessment by detecting fluid leakage and fractures. We analyzed borehole microseismic data acquired during CO 2 injection and migration, including data denoising, event detection, event location, magnitude estimation, moment tensor inversion, and stress field inversion. We detected and located two shallow clusters, which occurred during increasing injection pressure. The two shallow clusters were also featured by large b values and tensile cracking moment tensors that are obtained based on a newly developed moment tensor inversion method using single-borehole data. The inverted stress fields at the two clusters showed large deviations from the regional stress field. The results provide evidence for microseismic responses to CO 2 /fluid injection and migration.

02 PETROLEUM↗

Autonomous Inversion of In Situ Deformation Measurement Data for Injection-Induced Stress Change

Geologic carbon storage (GCS) is likely to play a key part of the global effort to dramatically reduce CO2 emissions and perhaps even reduce atmospheric CO2 concentrations through carbon negative operations. A critical part of effort to commercialize and widely deploy this technology is developing the capability to rapidly assimilate real-time monitoring data into a form that will enable site operators to make decisions to manage the safe and efficient operations. Two of the risks associate with GCS are the risk of inducing fractures in the sealing formations that can create leakage pathways and the risk of inducing earthquakes of sufficient magnitude to cause public concern, property damage, or safety risks. To properly manage these risks the site operator needs to know the initial state of stress, the change in stress induced by injection, and the relationship between operational parameters such as injection rate and pressure and the change in stress. Current methods of estimating the change in stress require choosing the type of constitutive model and the model parameters based on core, log, and geophysical data during the characterization phase, with little feedback from operational observations to validate or refine these choices. These characterization methods interrogate the geologic formations using length scales, loading rates or magnitudes that are quite different from those encountered by the actual storage system. It is shown that errors in the assumed constitutive response, even when informed by laboratory tests on core samples, are likely to be common, large, and underestimate the magnitude of stress change caused by injection. Recent advances in borehole-based strain instruments and borehole and surface-based tilt and displacement instruments have now enabled monitoring of the deformation of the storage system throughout its operational lifespan. This data can enable validation and refinement of the knowledge of the geomechanical properties and state of the system, but brings with it a challenge to transform the raw data into actionable knowledge. We demonstrate a method that uses automatic differentiation and a finite-element based geomechanical model perform a gradient-based deterministic inversion of geomechanical monitoring data. This approach allows autonomous integration of the instrument data without the need for time consuming manual interpretation and selection of updated model parameters. Furthermore, only isotropic linear elasticity is considered in this paper, the approach presented is very flexible as to what type of geomechanical constitutive response can be used. The approach is easily adaptable to nonlinear physics-based constitutive models to account for common rock behaviors such as creep and plasticity. The approach also enables training of machine learning-based constitutive models by allowing back propagation of errors through the finite element calculations. This enables strongly enforcing known physics, such as conservation of momentum and continuity, while allowing data-driven models to learn the truly unknown physics such as the constitutive or petrophysical responses.

Burghardt, Jeffrey A.↗

Electrical Resistivity Tomography based monitoring of stress perturbations to optimize placement of high-precision strain meters

The Center for Understanding Subsurface Signals and Permeability is a new U.S. Department of Energy Earthshot Center focused on understanding and predicting the long-term evolution of permeability in enhanced geothermal systems. The center will use a highly instrumented testbed within the Sanford Underground Research Facility to conduct field scale experiments that elucidate and test capabilities to simulate geochemical-geomechanical interactions and permeability evolution. Here we demonstrate initial developments using previously collected electrical resistivity tomography (ERT) monitoring data with high-performance multi-physics modelling advancements to inform the optimal location of two new monitoring boreholes. Specifically, ERT monitoring data collected during shear stimulation testing shows marked responses to changes in stress during borehole pressurization. We demonstrate how the same response is being simulated, ultimately to train a machine-learning algorithm to estimate rock properties and enable enhanced prediction of stress and strain responses anticipated during future testing campaigns.

Stress, EGS, CUSSP, 3D Electrical Imaging↗

Insights into the hydro-mechanical behavior of a decimeter-scale fracture using the mini-SIMFIP probe

ABSTRACT: Understanding hydro-mechanical couplings in fractured rocks is essential for predicting the rock mass response during high-pressure fluid injection, including the stimulation of enhanced geothermal systems. However, fluid-driven fracture dislocations are challenging to measure due to the need for local displacement data at high fluid pressures. In this study, fluid-driven displacement across a decimeter-scale laboratory fracture was investigated using the mini-SIMFIP (step rate injection method for fracture in-situ properties) probe, which is a smaller version of the SIMFIP tool (Guglielmi et al., 2014). The mini-SIMFIP probe is able to resolve the full 3D displacement vector of a fracture. The probe was installed in one of two boreholes across a decimeter-scale saw-cut granite fracture. Two pressure step injection tests were conducted under the same isotropic stress conditions. By varying injection between the boreholes, we estimated the aperture profile across the fracture. The results show a consistent pressure-dependent opening as long as steady-state flow was maintained. At a certain pressure step, steady-state conditions were no longer achievable and the fracture opened rapidly. The pressure-opening relationship diverged between the two tests and indicated a homogenization of the aperture profile and an increase in the non-linearity of the flow regime.

Osten, J↗

System for detecting substructure microfractures and method therefore

Bursts of signals at different frequencies are induced into substructure, adjacent to a borehole. The return signals from each burst of signals are normalized to compensate for the attenuation, experienced by more distant return signals. The peak amplitudes of return signals, above a selected level, are cut off, and an average signal is produced from the normalized amplitude-limited return signals of each burst. The averaged signals of the return signals of all the signal bursts at the different frequencies are processed to provide a combined signal, whose amplitude is related to the microfracture density of the substructure adjacent to the borehole.

Parthasarathy, S. P.↗

Geological assessment probe

A probe is described which can be installed in a side hole that extends from a bore hole in the Earth, to assess the permeability of the strata surrounding the borehole. The probe is elongated and has a plurality of seals spaced therealong and sealed to the walls of the side hole to form a plurality of chambers sealed from one another. A tracer fluid injector on the probe can inject a tracer fluid into one of the chambers, while a tracer fluid detector located in another chamber can detect the tracer fluid, to thereby sense the permeability of the strata surrounding the side hole. The probe can include a train of modules, with each module having an inflatable packer which is inflated by the difference between the borehole pressure and the strata pressure.

Collins, E. R.↗

Mechanical and hydrologic basis for the rapid motion of a large tidewater glacier. 1: Observations

Measurements of glacier flow velocity and basal water pressure at two sites on Columbia Glacier, Alaska, are combined with meteorological and hydrologic data to provide an observational basis for assessing the role of water storage and basal water pressure in the rapid movement of this large glacier. During the period from July 5 to August 31, 1987, coordinated observations were made of glacier surface motion and of water level in five boreholes drilled to (or in one case near to) the glacier bed at two sites, 5 and 12 km from the terminus. Glacier velocities increased downglacier in this reach from about 4 m/d to about 7 m/d. Three types of time variation in velocity and other variables were revealed: (1) Diurnal fluctuation in water input/output, borehole water level, and ice velocity (fluctuation amplitude 5 to 8%); (2) Speed-up events in glacier motion (15-30% speed-up), lasting about three days, and ocurring at times of enhanced input of water, in some cases from rain and in others from ice ablation enhanced by strong, warm winds; (3) 'Extra-slowdown' events, in which, after a speed-up event, the ice velocity decreased in about 3 days to a level consistently lower than that prior to the speed-up event. All of the time variations were due, directly or indirectly, to variations in water input to the glacier.

Meier, Mark↗

Osmium Isotope Constraints on the Proportion of Bolide Component in Chicxulub Impact Melts

Recent drilling has produced the samples of impact melt breccias examined here. All samples examined here are from the YAX-1 borehole [7]. Two main types of impact breccias have been studied. The first type is a green altered impact meltrock found in the lower portion of the impact sequence. The texture of the rock is microcrystalline and is composed of pyroxene, plagioclase, and alkali feldspars. Its composition is consistent with continental margin rocks. It is generally massive with some flow structure. The rock was brecciated and altered after solidfication and contains small amounts of both shocked and unshocked clasts of the impacted lithologies. These lithologies include lithic quartzite, and isolated feldspar crystals. The compositions of these rocks are similar to those seen in meltrocks sampled by the Yucatan-6 borehole [8-9]. Our study includes samples YAX-1_861.4, YAX-1_863.51, and YAX- 1_876.46, which represent both the top and lower portion of the green impact meltrock. The middle sample in the sequence has the least amount of (mineralogical) alteration [10]. The second type of melt breccia under study is a brown altered impact meltrock. It also has a microcrystalline texture and both shocked and unshocked clasts of the target material. Even though this rock type has been altered, remnant schleiren, metaquartzite, and micritic calcite have been identified. Sample YAX-1_841.32 is representative of this type of rock. It was recovered from a polymict breccia in the middle of the impact sequence.

Gelinas, A.↗

Construction and Resource Utilization Explorer (CRUX): Implementing Instrument Suite Data Fusion to Characterize Regolith Hydrogen Resources

CRUX is a modular suite of geophysical and borehole instruments combined with display and decision support system (MapperDSS) tools to characterize regolith resources, surface conditions, and geotechnical properties. CRUX is a NASA-funded Technology Maturation Program effort to provide enabling technology for Lunar and Planetary Surface Operations (LPSO). The MapperDSS uses data fusion methods with CRUX instruments, and other available data and models, to provide regolith properties information needed for LPSO that cannot be determined otherwise. We demonstrate the data fusion method by showing how it might be applied to characterize the distribution and form of hydrogen using a selection of CRUX instruments: Borehole Neutron Probe and Thermal Evolved Gas Analyzer data as a function of depth help interpret Surface Neutron Probe data to generate 3D information. Secondary information from other instruments along with physical models improves the hydrogen distribution characterization, enabling information products for operational decision-making.

in situ resource utilization (ISRU)↗

Subsurface Ice Probe

The subsurface ice probe (SIPR) is a proposed apparatus that would bore into ice to depths as great as hundreds of meters by melting the ice and pumping the samples of meltwater to the surface. Originally intended for use in exploration of subsurface ice on Mars and other remote planets, the SIPR could also be used on Earth as an alternative to coring, drilling, and melting apparatuses heretofore used to sample Arctic and Antarctic ice sheets. The SIPR would include an assembly of instrumentation and electronic control equipment at the surface, connected via a tether to a compact assembly of boring, sampling, and sensor equipment in the borehole (see figure). Placing as much equipment as possible at the surface would help to attain primary objectives of minimizing power consumption, sampling with high depth resolution, and unobstructed imaging of the borehole wall. To the degree to which these requirements would be satisfied, the SIPR would offer advantages over the aforementioned ice-probing systems.

Hecht, Michael↗

Repurposing Drilling Control Diagnostics for Subsurface Edge Detection and Boundary Advisement During Planetary Drilling

Informed decision-making during lunar drilling and sampling missions will require data monitoring tools and specialized ground data systems. Accurate and updated situational awareness, with ongoing data monitoring, is critical for timely responses by to incoming science data. Traverse plans and scheduled activities may need to be flexibly changed in order to react to unexpected data or situations. Unlike (for example) Mars missions, the relative lightspeed closeness of the Moon allows for near-real-time ground processing of incoming mission and instrument data. An Apollo-class lunar regolith drill will in a sense “travel” a meter or two vertically at a given subsurface characterization site. As the drill penetrates into lunar regolith, it is likely to encounter a range of material densities, orientations, fracture toughness, and (perhaps) ice percentages. Lunar drill telemetry can provide science teams with a valuable first look into the subsurface structure, the regolith bulk properties, and constituents at each drilled site. Real-time AI-based recognition and reaction to downhole situations has been developed for automated deeper drilling on Mars and beyond. We can leverage the same knowledge bases and pattern-matching as areal-time interpreter of the subsurface, a situational awareness tool during drilling operations. We recently (Sept. 2019) demonstrated this AI drilling monitoring and analysis capability, in control of in-situ drilling and sampling operations, mounted on a KREX-2 rover in Chile’s Atacama Desert. Terrestrial automated drilling log analyses in oil exploration have used similar machine learning techniques in classifying and identifying features in drilling logs –but these typically are designed assuming a drilling fluid influencing downhole measurements and data (permeability, resistivity). Drilling models and existing AI software designed to detect and respond to drilling faults and hard materials can be repurposed, for near-real-time (ground-based) interpretation of drilling telemetry –a potentially valuable advisory tool for strata boundaries and changes in drilling parameters. On the Moon, this approach could be used to study the structure and to some extent the composition of lunar regolith vs. borehole depth, based on recognizable variations in fracture hardness, drilling energy and penetration rates while actively drilling. Since the early 2000s, a series of increasingly-capable real-time drilling telemetry interpretation and characterization software tools have been developed. These subsurface models and software tools have monitored the real-time drilling data received, and automatically identified changes in drill behavior (e.g., encountering a harder target layer, bit inclusions, drill choking due to infall downhole, and others) correlating these with subsurface structures and features. We discuss the mappings between drill borehole parameters, faults or events detected, and modeled changes in rock layer boundaries, in examples drawn from field testing at analog sites in an Arctic impact crater, Rio Tinto, and Chile’s Atacama Desert. These demonstrate how subsurface structural boundaries led to fault detections and responses by the software.

drilling advisor↗

TRIDENT Drill Validation at Mars and Lunar Analog Field Sites

Drilling on Earth is typically a human-intensive activity. Drilling on other planets is further complicated by the lack of prior local field surveys of their target area, hence blindly drilling into uncertain target rocks. Field conditions on the Moon or Mars are also different than for shallow drilling on Earth: lower temperatures and pressures, less power available, low masses (hence less weight-on-bit). Given the cost of transport from Earth, no drilling muds or working fluids are likely to be available to carry away cuttings. And impact-gardened regolith and dust vary mechanically and texturally from most terrestrial soils. The Regolith and Ice Drill for Exploration of New Terrains (TRIDENT) is a rotary-percussive 1m-class drill from Honeybee Robotics. It is low-power (rotary and percussive actuators are 200 W each) and lightweight (<20 kg) with the maximum weight on bit limited to 200 N. TRIDENT has been manifested for the Volatiles Investigating Polar Exploration Rover (VIPER) and PRIME-1 lunar south pole missions in 2024, has previously been field tested at a hot, dry analog site in the Atacama Desert, and in lunar conditions in thermal vacuum chamber tests. TRIDENT was also part of the 2019 Icebreaker Mars Discovery proposal, as well as in the Mars Life Explorer concept. During ARADS tests, drill control and fault recovery automation software enabled hands-off operations of a rover-mounted TRIDENT drill. TRIDENT Drill Analog Site Validation: Past TRIDENT tests in thermal vacuum (TVAC) chambers targeted containers of manufactured lunar simulants with added volatiles. 2022 TRIDENT ambient testing at NASA Ames drilled into cemented lunar simulant materials. Low cuttings-permeability led to cuttings buildup, and drill choking and binding was observed. The Atacama analog site in ARADS had desiccated unconsolidated sediments that did not challenge the TRIDENT design. However, lunar polar regolith is expected to be diverse and heterogeneous with varying clast sizes, with abundant impactites and perhaps subsurface ice deposits. Neither the simulants nor Atacama testing had completely covered the TRIDENT-targeted field characteristics, motivating further analog tests prior to the planned lunar missions. To gain more insight into the behavior of the TRIDENT hardware in diverse impactites and subsurface ice, and to verify the software automation in that environment, in August 2023 TRIDENT was brought to Haughton Crater, a field analog site in the Canadian Arctic. In September 2023 the same drill was brought to the Bishop Tuff in southern California to verify whether drilling binding behaviors previously seen in lunar simulant testing would be observed in naturally occurring fine-grained massive layers. The ~22 Ma Haughton Crater impact structure is located at 75 ̊22’ N, 89 ̊41’ W, on northwestern Devon Island, Nunavut, Canada. Numerous deposits of pale-grey crater-fill polymictic impact-melt breccia are found within the crater with a typical thickness reaching ~125 m or greater and covering ~60 km -2 . An approximately 600m-thick permafrost layer is also present with ice typically found within 0.5-0.6m of the surface. The volcanic tableland north of Bishop, CA exposes densely welded tuff laid down during the eruption that created the Long Valley Caldera at approximately 0.76 Ma. Extensional faults and the Owens River gorge expose cross-sections across the plateau. The area is viewed as an analog site for Mars features believed to be of pyroclastic origin. Results: Haughton Crater.Drilling tests were conducted 8-13 August 2023 at a previously undisturbed area separated 5-10 m from past years’ Drill Hill test sites (75.4208, -089.7613). In six days, TRIDENT drilled 8 holes to nearly 1 m depth each, totaling 7.80 m. The active layer/ice boundary was at ~67 cm depth, with a total of approximately 2.4 m drilled into ice or ice-cemented impact breccia. During drilling, five drill fault states were observed and successfully recovered. Holes 23-1, 4 and 7 were drilled under manual control, using Honeybee’s Thorax user interface. Holes 23-5, 6, and 8 were drilled with the Ames IBexec automated drilling control software. TRIDENT was observed to have little difficulty in the thawed uncemented impact breccia above the active layer boundary, but required percussion to make slower headway in the ice-cemented breccia. In Hole 23-7 (Fig.1), drilling slowed down in a massive unit just above the active-layer boundary (perhaps a large rock extending into the ice-cementation?), with only 7cm progress made in 27 minutes of high auger torque and constant percussion, leading to a choking fault and then a binding fault. A similar pattern had been observed in TRIDENT Rio Tinto test data from 2017 [6] as well as in the 2022laboratory tests.Bishop Tuff.A team from NASA Ames and the US Geological Survey deployed the same TRIDENT drill to Bishop Sites 1B and 1C (37.4203, -118.4289; 37.4265, -118.4215) on 13-16 September 2023, on the Bishop Tuff plateau. A third drill site was used 17-18 September 2023(37.4598, -118.3667) in an abandoned pumice mine. Four holes (totaling 2.5m depth) were drilled into the fine-grained, meters-thick tuff units at Sites 1B and 1C, and a further two boreholes (totaling 1.98m depth) were made at the pumice site. Drill behavior in the tuff below 10 cm depth was similar to that seen at 65-74 cm depth in Haughton Hole 23-7 (Fig. 1) and that seen in the 2022 lab simulant drilling. Drill safety torque limits were exceeded multiple times resulting in drill stops downhole. These freezes then required external added torques (with a pipe wrench) to resume rotation, to unstick the drill for withdrawal. To prevent this choking and binding behavior we found that more-frequent cuttings removal was necessary, e.g., reducing the “drill bite” size from the nominal 10 cm to 2 cm per bite --bringing the auger up to the surface more frequently, as seen in Bishop Site 1C Hole 2 (Fig.2). This permitted slow progress without drill binding and without external interventions. Conversely, TRIDENT drilling in the more porous pumice target material showed no cuttings buildup issue, and single bites as large as 40 cm were demonstrated. Discussion: We observed that TRIDENT easily penetrated unconsolidated heterogeneous soils (both above the active layer boundary at Haughton and previously in the Atacama). Cemented or consolidated targets that were cuttings-permeable (icy impact breccia, pumice) required more energy applied and percussion. However, in non-cuttings-permeable targets (welded microporous tuff, cemented simulants, boulder) TRIDENT was observed to be prone to excessive cuttings accumulation leading to choking/binding faults and stalling. The wedge cutting bit, used by TRIDENT in field tests and in its flight versions, pulverizes the target rock and creates fine cuttings that ideally are transported up the auger spirals for removal. In porous, fractured or vesicular target materials (such as at the Bishop pumice site) a significant portion of the cuttings are pushed aside, but for non-fractured, microporous targets the cuttings remain in the borehole and accumulate. Rock powder is relatively incompressible as a working fluid at only 100-200N downward force (TRIDENT limits) and hence eventual drilling progress slows or stops. Our recommended strategy for improving TRIDENT cuttings removal in massive target units with low cuttings-permeability is to reduce TRIDENT bite sizes when encountering these units, from 10cm to as little as 1-2cm, to effectively bail the accumulating cuttings. This approach was demonstrated to reduce choking and allowed slow progress to continue in cuttings-impermeable microporous target units (viz. the Bishop tuff in our September 2023 tests or cemented simulants in 2022 ambient tests).

robotic drilling↗

Slip Tendency Analysis of Fracture Networks to Determine Suitability of Candidate Testbeds for the EGS Collab Hydroshear Experiment

Experiment 2 of the EGS Collab project is aimed at testing stimulation by hydro-shearing of existing natural fractures, versus Experiment 1, which is focused on hydraulic fracturing a rock mass to enhance permeability. The main criterion for the testbed selection in Experiment 2 is the presence of an interconnected network of fractures, at least 10 meters in extent, and with orientations such that shear slip can be induced at injection pressures less than the minimum horizontal stress (Shmin). The feasibility analysis for this experiment requires a well-constrained stress state along with well-characterized fracture networks. The fracture systems at two candidate locations in the Sanford Underground Research Facility (SURF), the 4850 Level (number refers to depth below ground surface in feet) and the 4100 Level, have been characterized to different extents. The 4850 Level has a well-characterized Discrete Fracture Network (DFN) from borehole (drilled for an experimental test bed) and drift observations. In contrast, the 4100 level has a fracture network characterized only by observations from the drift wall, which provides little constraint on fracture extents. This paper will present assessments for the interpreted slip potential for the natural fractures at the two locations. Data uncertainties are addressed by performing a probabilistic analysis that takes into account the uncertainty in the stress state, uncertainty in the fracture properties, and preliminary borehole locations. There is strong evidence of shear stimulation of some natural fractures in Experiment 1 on the 4850 level. We use examples of natural fractures strongly linked to shear stimulation in Experiment 1 from multiple monitoring indicators to test the consistency of the stress model and guide the uncertainty interpretation. This analysis is intended to guide the site selection process for Experiment 2 by highlighting the fracture orientations that are likely to be shear-stimulated in a majority of modeled realizations.

40 EE - Geothermal Technologies Office (EE-4G)↗

Thermal Stimulation and Injectivity Testing at Raft River, ID EGS Site

The Raft River geothermal field is the site of an innovative Department of Energy Enhanced Geothermal System (EGS) project to determine the viability of using combined thermal and hydraulic stimulation techniques to improve energy production. Well RRG-9 is currently undergoing a stimulation program using injectate from the US Geothermal Raft River Power Plant and cold water from a cooling tower make-up water well. The stimulation began on 13 June 2013 with injection from the power plant at a temperature of about 39 °C and a pressure of 275 psig. Next, two positive displacement plunger type pumps were used to increase the injection pressure and flow rate for about one month. The highest rate achieved was 258 gpm at a pressure of 741 psig. During this time, fluid from the cooler water well was injected for about 2 weeks at various pressures. Then, the pumps were removed and plant injection resumed on 25 September. Plant injection will continue until the spring of 2014, when a high pressure hydraulic stimulation will be conducted. A series of seismic monitoring stations deployed around the well are providing data on seismic events occurring at the site. Over the past year, 51 microseismic events have been recorded, all less than Magnitude 1. During injection, several diagnostic tests were conducted to gain a better understanding of the well and reservoir. A step-rate test was performed on 22 August to measure the in-situ stress and aid in modeling in-situ fractures. A tracer was injected into the well on 9 September. No tracer was detected in adjacent production wells after several months. A second borehole televiewer survey was conducted for comparison to pre-stimulation images. Here, a third borehole televiewer survey is planned after the high pressure stimulation. Injection test data is evaluated in real time. A modified Hall plot analysis indicates the effective permeability is increasing. The injectivity index supports the results of the modified Hall plot analysis. As the thermal stimulation has continued, the injectivity index has consistently followed an upward trend from 0.1 gpm/psi to 0.53 gpm/psi.

enhanced geothermal systems↗

The weathered bedrock vadose zone: A hidden control on water availability in the western United States (Final Technical Report)

This report summarizes the research completed during the project entitled "The weathered bedrock vadose zone: A hidden control on water availability in the western United States" (Award: DE-SC0018039). The report covers work completed between 9/15/17 and 9/14/20. The objective of the exploratory project work was to evaluate the overarching hypothesis that the subsurface underlying actively eroding landscapes is systematically structured as a result of a balance between weathering and erosion and that this structure leads to dynamic moisture storage within weathered rock. To accomplish project objectives, we characterized hillslope weathering profile structure and conducted subsurface moisture monitoring across four study sites with comparable bedrock types but differing erosion rate and hydroclimate: The Angelo Coast Range Reserve (ACRR), Sagehorn-Russell Ranch (SRR), Antelope Valley Ranch (AVR), and the LBNL Watershed Function Scientific Focus Area (SFA). New characterization and monitoring infrastructure at the AVR and SFA sites was established through this project and enables long-term monitoring of hydrologic and carbon fluxes at the sites for comparison with long-term monitoring in other sites. Deep drilling, downhole logging, and characterization of chemistry and porosity of rock recovered from drilling reveals two weathering fronts across the sites. A shallower front is marked by the extent of pervasive fracturing and oxidation of pyrite and organic carbon. A deeper weathering front is marked by the appearance of open fractures and discoloration. At the AVR site, multiple ridge tops of different hillslope length were investigated to establish the relationship between the weathering front and topography. Drilling reveals that the fraction of the hillslope relief that is weathered scales with hillslope length, consistent with the theoretical framework that guided the project relating groundwater drainage with landscape evolution. With the establishment of new boreholes, we have documented significant water storage in weathered bedrock using low-field borehole nuclear magnetic resonance (NMR) across all four sites. In the upper 4 m of bedrock, the average water contents at ACRR, AVR, SRR, and SFA are approximately 25%, 16%, 25%, 23% respectively. At the ACRR site, where long-term monitoring via neutron probe surveys and time-domain transmission sensors are available, we have demonstrated that NMR reliably records moisture content and importantly, moisture content dynamics in weathered bedrock. Across all sites, our NMR monitoring has also revealed that a considerable fraction of the water storage occurs in fractures or pores that are significantly larger than the fine-grained bedrock matrix of the parent rock. To identify the proportion of water storage occurring in fractures and large pores, we analyzed the NMR signal by evaluating both the sum of echoes and an inversion of the NMR signal for a distribution of T2 relaxation times. Both methods reveal that at least 20% of the water storage occurs in fractures and larger pores across all sites. To constrain the interpretation of the in-situ field-scale NMR measurements, we compare the NMR response in the field to independent pore-structure information obtained via helium pycnometry. Results support the inference that dynamic seasonal water storage is dominantly restricted to the fracture network, which can reach 7% of the total volume. This dynamic range in water content is comparable to that of some soils, underscoring the need to understand how this region of the Earth System is structured, and functions as a hydrologic reservoir.

58 GEOSCIENCES↗

Coupling subsurface and above-surface models for optimizing the design of borefields and district heating and cooling systems in the presence of varying water-table depth. In: Proceedings, 46th Workshop on Geothermal Reservoir Engineering

Dynamic energy simulation is important for the design and sizing of district heating and cooling systems with geothermal heat exchange. Current modeling approaches in building and district energy simulation tools typically consider heat conduction through the ground between boreholes, without flow of groundwater. On the other hand, detailed simulation tools for subsurface heat and mass transfer exist, but these fall short in simulating above-surface energy systems. To support the design and operation of such systems, we have developed a coupled model including a software package for building and district energy simulation, and software for detailed heat and mass transfer in the ground. For the first, we use the open-source Modelica Buildings Library, which includes dynamic simulation models for building and district energy and control systems. For the heat and mass transfer in the soil, we use the TOUGH simulator. TOUGH can model heat and multi-phase, multi-component mass transport for a variety of fluid systems, as well as chemical reactions, in fractured porous media. In previous work, we described the coupling of these software packages, including how time-dependent boundary conditions for the borehole walls are synchronized for use in Modelica and TOUGH. We verified that the coupled Modelica/TOUGH code produced consistent results with the original Modelica code for an idealized problem in which heat transfer was purely by conduction in a uniform geologic medium. Here, we examine less idealized problems for which TOUGH’s advanced capabilities for modeling fluid flow are required. The first problem has a shallow vadose zone, and the second problem has a thicker vadose zone with a water-table depth that varies in time, which requires a fine vertical grid discretization for the TOUGH model.

Doughtry, Christine↗

Experimental Evaluation of Thermal Storage Performance of a Dual-Purpose Underground Thermal Battery

A geothermal heat pump (GHP) system is an energy-efficient building heating and cooling technology with great potential for reducing energy consumption and decarbonization. However, applications of GHP are still limited due to the high cost, of which 30% is related to the cost of installing the conventional vertical bore ground heat exchangers, which are usually installed in boreholes 60 meters deep. A dual-purpose underground thermal battery (DPUTB) has been developed to offer a low-cost ground heat exchanger with a built-in thermal storage capacity. The DPUTB innovatively integrates a shallow-bore ground heat exchanger (the outer tank), which can be installed in a borehole less than 6 m deep, with thermal energy storage (TES) (the inner tank). DPUTB has the potential to reduce the cost of a ground source heat pump system while allowing shifting the electric demand of the building served by the GHP system from peak to off-peak hours of the electric grid by charging and discharging the thermal storage. A lab-scale (1:125 in volume) DPUTB prototype was built. Phase change material (PCM) was added to increase the thermal storage capacity and maintain the supply water temperature from the TES within the desired range for direct cooling operation during the discharge period. As PCMs are critical to the TES performance of the DPUTB, this study compared the influence of different PCMs (including salt hydrate and organic PCMs) on the discharge performance of the DPUTB. The thermal State of Charge (SoC) of the DPUTB was used to compare the performance resulting from using different PCMs. Test results indicate that the organic PCM (Methyl Laurate) outperforms salt-hydrate PCMs due to a lower melting temperature and a narrower melting temperature range during the phase change process. The results of this study provide a guide for PCM selection and the optimal design of DPUTB.

Wang, Lingshi↗

Development of New g-function Data for Simulating a Novel Shallow Bore Ground Heat Exchanger

Geothermal heat pump (GHP), which is also referred as ground source heat pump, is the most energy-efficient technology for space heating and cooling. However, the application of GHP is hindered by its high initial cost, of which approximately 30% is for the ground heat exchanger (GHE). In recent years, researchers have developed different types of shallow bore ground heat exchangers (SBGHE) intending to reduce the cost of drilling. The ability to predict the thermal response of an SBGHE is critical for sizing SBGHE. While g-functions have been commonly used for predicting the thermal response of conventional vertical bore ground heat exchangers (VBGHE), they cannot be directly used for predicting the performance of SBGHE because they did not account for the impacts of the seasonal variation of the soil temperature along with the depth of an SBGHE. In addition, an SBGHE has a larger thermal mass within the borehole than the VBGHE due to the larger borehole diameter of SBGHE. This study develops new g-function data pairs for a new design of SBGHE, which is named Underground Thermal Battery (UTB). Impacts of the seasonal variation of soil temperature on the thermal response of UTB were accounted for by superposing a time-dependent soil temperature onto the g-functions calculated with a numerical model that assumes constant undisturbed soil temperature. The TOUGH program was used to predict the thermal response of several configurations involving multiple UTBs. The results indicate that the proposed methodology is appropriate to generate g-functions for the UTB, and the g-function value of UTB is much lower than that of VBGHE in the time range of 15 min to 1 year due to the large thermal mass and convection heat transfer within the UTB.

Shi, Liang↗