Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “fracturing fluid”

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 325 records · Page 18

Inorganic Alterations in Unconventional Shale Reservoirs: Importance of Additive and Base Fluid Chemistry

The effective development of unconventional petroleum systems requires the use of significant water resources. In an effort to reduce the consumption of freshwater resources for hydraulic fracturing, highly saline produced waters are increasingly recycled for use as a base fluid. However, there are significant knowledge gaps regarding potential water–rock interactions resulting from the introduction of produced waters and associated additives into shale reservoirs such as formation and deposition of mineral scale, which can negatively affect hydrocarbon production through wellbore restriction and damage to hydraulically generated fractures. To assess the impacts of field stimulation practices in the subsurface, a series of laboratory experiments were completed using (a) three distinct sedimentary rock formations of the Midland Basin (Texas, USA) and (b) additives with two different base fluids: municipal fresh water and clean brine. The experimental approach used relevant injection sequences and mixing ratios in specialized reactors for 3 weeks. Static pressurized experiments and nonpressurized time-resolved experiments were undertaken. The resulting solids and liquids were analyzed by using a variety of laboratory- and synchrotron-based techniques. The use of an acid spearhead (15% HCl) resulted in texturing of both clay-rich and calcareous shales, which can temporarily enhance porosity but subsequently result in mineral scale deposition. The primary matrix scale was Fe(III)-bearing phases, which occurred in all experiments regardless of base fluid chemistry. Additionally, strontium sulfate (SrSO 4 ) precipitated on shale surfaces when clean brines were used. It was concluded that clean brine was the main source of Sr 2+ species, while persulfate breaker degradation and oxidation of pyrite were the sources of SO 4 2– . Sulfate scaling was more pronounced in clay-rich shales, suggesting that Sr sorption is important for promoting celestite formation. This work demonstrates that mineral scale deposition is a complex phenomenon, whereby the type and proportions of various mineral phases are determined from reservoir alteration processes and coprecipitation of constituents from injection fluids. In conclusion, the experimental results shown here should be considered when evaluating different base fluids and additives in order to mitigate mineral precipitation in unconventional shale reservoirs, which could result in reservoir degradation.

Jew, Adam D. [SLAC National Accelerator Laboratory↗

HOSS: Hybrid Optimization Software Suite

HOSS software, developed at LANL, is a powerful, flexible, and parallelized package that enables the user to predict deformation and failure of materials in a variety of situations. LANL’s HOSS software is an R&D 100 finalist that integrates solid mechanics and computational fluid dynamics into a multi-physics platform able to handle millions of interacting and fracturing solids and/or particle systems. This nextgeneration simulation software integrates solid mechanics and computational fluid dynamics into a multi-physics platform able to provide state-of-the-art simulation solutions for a wide field of applications including oil and gas, construction, mining, defense, material sciences, and biomechanics.

36 MATERIALS SCIENCE↗

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↗

Heterogeneous multiphase flow properties of volcanic rocks and implications for noble gas transport from underground nuclear explosions

Of interest to the Underground Nuclear Explosion Signatures Experiment are patterns and timing of explosion-generated noble gases that reach the land surface. The impact of potentially simultaneous flow of water and gas on noble gas transport in heterogeneous fractured rock is a current scientific knowledge gap. This article presents field and laboratory data to constrain and justify a triple continua conceptual model with multimodal multiphase fluid flow constitutive equations that represents host rock matrix, natural fractures, and induced fractures from past underground nuclear explosions (UNEs) at Aqueduct and Pahute Mesas, Nevada National Security Site, Nevada, USA. Capillary pressure from mercury intrusion and direct air–water measurements on volcanic tuff core samples exhibit extreme spatial heterogeneity (i.e., variation over multiple orders of magnitude). Petrographic observations indicate that heterogeneity derives from multimodal pore structures in ash-flow tuff components and post-depositional alteration processes. Comparisons of pre- and post-UNE samples reveal different pore size distributions that are due in part to microfractures. Capillary pressure relationships require a multimodal van Genuchten (VG) constitutive model to best fit the data. Relative permeability estimations based on unimodal VG fits to capillary pressure can be different from those based on bimodal VG fits, implying the choice of unimodal vs. bimodal fits may greatly affect flow and transport predictions of noble gas signatures. The range in measured capillary pressure and predicted relative permeability curves for a given lithology and between lithologies highlights the need for future modeling to consider spatially distributed properties.

45 MILITARY TECHNOLOGY, WEAPONRY, AND NATIONAL DEF↗

Hydraulic characterization of a fault zone from fracture distribution

A quantitative assessment of how faults control the migration of geofluids is critical in many areas of geosciences. Here we integrated geological fieldwork, quantitative analysis of the fractures distribution and numerical modeling to build a geometrical representation of a fault zone and to characterize its hydraulic properties. Our target is a fault located in the Majella Mountain (Italy). We collected 21 scan lines across the fault profile in order to characterize its architecture. The numerical modeling of the fracture network of the damage zones and their hydraulic parameters was performed using both commercial (Move ® ) and open source software (dfnWorks and PFLOTRAN). Move ® was used to build a representative model of the fault zone using fracture spacing as a proxy, and to model the hydraulic parameters of the different fault domains. dfnWorks and PFLOTRAN were employed to infer the hydraulic parameters of the damage zones of the fault and then upscale these properties to an equivalent continuum domain, suitable for fluid flow simulations through the whole fault zone. Our findings show how even in a relatively small area it is possible to describe changes in terms of hydraulic properties of a fault zone and to build models capable to represent these variations.

58 GEOSCIENCES↗

Hawthorne Nevada Deep Direct-Use Feasibility Study - Data Used for Geothermal Resource Conceptual Modeling and Power Capacity Estimates

This data submission includes several data components that were used to develop a conceptual model and power capacity-estimates of two low-temperature geothermal resources that define geothermal prospect A at Hawthorne, Nevada. Data are sourced from a combination of legacy publicly-available data and more recent data acquisition conducted by the US Navy Geothermal Program Office (2008-2013) and the Great Basin Center for Geothermal Energy at the University of Nevada, Reno (2008-2010). Data sets include compiled fluid geochemistry data, down-hole temperature logs for wells in the vicinity of prospect A, 2 meter temperature survey data, temperature-spinner logs acquired in well HWAAD-2A, fracture picks from image log data acquired in wells HWAAD-2 and HWAAD-3, and X-Ray Diffraction (XRD) analyses on cuttings from wells HWAAD-2A and HWAAD-3. These data have been reviewed for errors and inconsistencies, but it is possible that few errors could still remain. The resource conceptual model and power capacity estimates are included in the final report to the US Department of Energy, and are presented in a manuscript by Ayling and Hinz. A link to the manuscript published in Geothermics is linked below in this submission.

15 GEOTHERMAL ENERGY↗

Experimental and numerical investigation of proppant embedment and conductivity reduction within a fracture in the Caney Shale, Southern Oklahoma, USA

The current worldwide energy supply is insufficient to meet the rising demand. As a result, the energy prices are expected to keep soaring despite the recent increases in a variety of renewable energy resources. Although not renewable, shale oil and gas — “unconventional” hydrocarbon resources are relatively clean forms of energy resources, which still hold a vast share of the energy market. For many oil and gas companies, meeting profitable production goals from shale reservoirs is sometimes challenging, due to the loss of fracture conductivity and premature declines in the production. In this paper we investigate the stress-dependent changes in the hydraulic conductivity of proppant-filled fractures and mechanical fracture–proppant interactions in Caney Shale, a calcareous, organic-rich mudrock, through laboratory experiments and numerical modeling. American Petroleum Institute (API) fracture conductivity tests were conducted using 2% KCl on five locations within the Caney Shale that consisted of selecting three brittle (reservoir) zones and two ductile zones. Confining pressures ranged from 1,000 psi to 12,000 psi at 210 °F. Conductivity, permeability as well as embedment were measured during the test. Also, an additional, laboratory in-situ visualization test was conducted to examine the detailed proppant-shale matrix interaction under elevated stress (3,920 psi effective stress) and temperature (252 °F), with a synthetic reservoir fluid. Our experimental results have confirmed that improved fracture conductivity is attributed to proppant size, and that the increase in porosity of the proppant pack, closure pressure changes and the reduction in fracture conductivity are a function of many factors such as fracture closure stress.

04 OIL SHALES AND TAR SANDS↗

Evaluating proxies for the drivers of natural gas productivity using machine-learning models

We report the extensive development of unconventional reservoirs using horizontal drilling and multistage hydraulic fracturing has generated large volumes of reservoir characterization and production data. The analysis of this abundant data using statistical methods and advanced machine-learning (ML) techniques can provide data-driven insights into well performance. Most predictive modeling studies have focused on the impact that different well completion and stimulation strategies have on well production but have not fully exploited the available in situ rock property data to determine its role in reservoir productivity. We have used machine-learning techniques to rank rock mechanical properties, microseismic attributes, and stimulation parameters in the order of their significance for predicting natural gas production from an unconventional reservoir. The data for this study came from a hydraulically fractured well in the Marcellus Shale in Monongalia County, West Virginia. The data classes included measurements aggregated by well completion stage that included (1) gas production, (2) well-log-derived measurements including bulk density, elastic moduli, shear impedance, compressional impedance, brittleness, and gamma measurements, (3) microseismic attributes, (4) long-period long-duration (LPLD) event counts, (5) fracture counts, and (6) stimulation parameters that included the fluid injection volume and average pumping pressure. To identify observable proxies for the drivers of gas production, we evaluated five commonly used ML approaches including multivariate adaptive regression spline, Gaussian mixture model, random forest, gradient boosting, and neural network. We selected five variables including LPLD event count, seismogenic b-value, hydraulic diffusivity, cumulative moment, and fluid volume as the features most likely to impact gas productivity at the stage level in the study area. The data-driven selection of these parameters for their importance in determining gas production can help reservoir engineers design more effective hydraulic-fracture treatments in the Marcellus Shale and other similar unconventional reservoirs. Plain language summary: We use machine-learning methods and data-driven selection of reservoir parameters to rank and better understand their importance in determining gas production, which can help reservoir engineers design more effective hydraulic-fracture treatments in the Marcellus Shale and other similar unconventional reservoirs.

58 GEOSCIENCES↗

Coupled Thermo–Hydro–Mechanical Modeling of Hydro-Shearing Stimulation in an Enhanced Geothermal System in the Raft River Geothermal Field, USA

The Raft River Enhanced Geothermal System (EGS) demonstration project aims to improve the geothermal production by enhancing the reservoir permeability via shear stimulation in well RRG-9. In this paper, we performed a series of 3D thermo– hydro–mechanical (THM) simulations to investigate the influence of hydraulic and thermal effects on the development of this EGS. The model includes synthetic fracture populations based on borehole televiewer images and in-situ stress measurements from well RRG-9. Fracture permeability evolution is determined using an empirical permeability law developed from laboratory experiments. The model was calibrated by comparing the hydraulic response to field observations, including wellhead pressure, injection rate, and well injectivity. Particularly, we analyzed the enhancement of reservoir permeability and the spatial extent of the stimulation zone for the given injection schedule. Our results indicate that the permeability enhancement of fractured geothermal reservoir is caused by the combined effects of injection-induced cooling and fluid pressure increase. The decrease of temperature plays a dominant role in reactivation of natural fractures under the hydro-shearing mechanism, while the higher injection pressure promotes shear failure and enlarges the stimulation zone. Finally, for the specified extensional stress state, the model favors greatest permeability enhancement along the maximum horizontal principal stress, a moderate vertical enhancement, and a smaller gain along the minimum horizontal principal stress.

58 GEOSCIENCES↗

3-D Geologic Controls of Hydrothermal Fluid Flow at Brady geothermal field, Nevada, USA

In many hydrothermal systems, fracture permeability along faults provides pathways for groundwater to transport heat from depth. Faulting generates a range of deformation styles that cross-cut heterogeneous geology, resulting in complex patterns of permeability, porosity, and hydraulic conductivity. Vertical connectivity (a throughgoing network of permeable areas that allows advection of heat from depth to the shallow subsurface) is rare and is confined to relatively small volumes that have highly variable spatial distribution. This local compartmentalization of connectivity represents a significant challenge to understanding hydrothermal circulation and for exploring, developing, and managing hydrothermal resources. Here, we present an evaluation of the geologic characteristics that control this compartmentalization in hydrothermal systems through 3-D analysis of the Brady geothermal field in western Nevada. A published 3-D geologic map of the Brady area is used as a basis to develop structural and geological variables that are hypothesized to control or effect permeability or connectivity. The 3-D distribution of these variables is compared to the distribution of productive and non-productive fluid flow intervals along production wells and non-productive wells via principal component analysis (PCA). This comparison elucidates which geologic and structural variables are most closely associated with productive fluid flow intervals. Results indicate that production intervals at Brady are located: (1) within or near to known and stress-loaded macro-scale faults, and (2) in areas of high fault and fracture density.

58 GEOSCIENCES↗

Continuous surface-to-distributed acoustic sensor snapshots explain reactivation of individual natural fractures during an unconventional reservoir stimulation

ABSTRACT Fiber-optic sensing technologies allow petroleum engineering teams to detect hydraulic fracture interaction with boreholes during unconventional reservoir stimulation. In combination with high-repeatability seismic sources, the same distributed acoustic sensors (DASs) enable vertical seismic profiling (VSP) of the fracture evolution away from the boreholes. We discovered clear signatures of seismic scattering on activated fractures during nine days of continuous seismic monitoring of the fracturing stages at the Austin Chalk/Eagle Ford Field Laboratory. The present study applies a novel approach for quantitative analysis of the scattering events in terms of the evolution of the geometry and elastic stiffness of individual fractures. Our characterization strategy sequentially refines the fracture models: from a stack of 1D soft layers to 3D rectangular inclusions. First, we estimate the number of fracture locations and reflectivity using a modified sparse-spike deconvolution of the stacked VSP traces. The fracture set consists of five fractures spaced by 15–30 m with a reflectivity of approximately 1%. Then, we develop a scattering integral method to refine these estimates along with an inversion of the fracture top and bottom for each monitoring vintage. We find that, initially, some of the fractures are located above the monitoring fiber with the height of approximately 100 m. Then we integrate the seismic interpretation with the low-frequency DAS and pressure and microseismic monitoring to reconstruct the activation process of the fractures. Most likely, some of the natural fractures slowly grew downward to the monitoring fiber as a result of fluid injections in the stimulated well. This led to bright strain anomalies but did not trigger seismicity. The top of the fractures remained almost constant and were limited by a lithologic boundary/stress barrier. To our knowledge, this is the first time VSP data enabled tracking of the fracture evolution with such high spatial and temporal resolution, which was previously only available for crosswell surveys and at a much smaller scale.

Glubokovskikh, Stanislav↗

Influence of extreme fracture flow channels on the thermal performance of open-loop geothermal systems at commercial scale

Abstract Adequate stewardship of geothermal resources requires accurate forecasting of long-term thermal performance. In enhanced geothermal systems and other fracture-dominated reservoirs, predictive models commonly assume constant-aperture fractures, although spatial variations in aperture can greatly affect reservoir permeability, fluid flow distribution, and heat transport. Whereas previous authors have investigated the effects of theoretical random aperture distributions on thermal performance, here we further explore the influence of permeability heterogeneity considering field-constrained aperture distributions from a meso-scale field site in northern New York, USA. Using numerical models of coupled fluid flow and heat transport, we conduct thermal–hydraulic simulations for a hypothetical reservoir consisting of a relatively impervious porous matrix and a single, horizontal fracture. Our results indicate that in highly channelized fields, most well design configurations and operating conditions result in extreme rates of thermal drawdown (e.g., 50% drop in production well temperatures in under 2 years). However, some other scenarios that account for the risks of short-circuiting can potentially enhance heat extraction when mass flow rate is not excessively high, and the direction of geothermal extraction is not aligned with the most permeable features in the reservoir. Through a parametric approach, we illustrate that well separation distance and relative positioning play a major role in the long-term performance of highly channelized fields, and both can be used to help mitigate premature thermal breakthrough.

Rangel-Jurado, Nicolás↗

Revealing complex subsurface dynamics with continuous seismic monitoring: Observations using distributed acoustic sensing and surface orbital vibrators during hydraulic fracturing

Understanding hydraulic fracturing is crucial to improving the stimulation of unconventional reservoirs and increasing fluid production. This study develops a novel seismic monitoring technology using distributed acoustic sensing (DAS) and surface orbital vibrators (SOV) to capture fracture seismic response and mechanical properties at high temporal intervals. We analyze continuous time-lapse vertical seismic profiling (VSP) data acquired every hour during the first nine days of treatment of an unconventional reservoir in the Austin Chalk/Eagle Field Laboratory. The VSP data contain clear seismic signals scattered from the activated fractures. The spatiotemporal changes of the fracture reflectivity revealed by the SOV/DAS data correlate well with the observations of fracture locations inferred from low-frequency DAS data. These results capture the fracture opening and closure processes, as well as highlight potential prestage activations of the fractures due to hydraulic connectivity with preexisting fracture systems. Therefore, analysis of the presented data set provides a unique opportunity to understand fracture initiation and subsequent evolution, not only in the context of unconventional resources but also in enhanced geothermal systems.

Correa, Julia↗

Detecting and Characterizing Fracture Zones Using a Convolutional Neural Network

This project directly supports the Geothermal Technologies Office (GTO) objectives outlined in the Multi-Year Program Plan (MYPP) by advancing two key research areas: “Exploration and Characterization” and “Data, Modeling, and Analysis.” This project has successfully demonstrated a pre-drilling ability to image and characterize the distribution and connectivity of subsurface faults and fractures, key parameters for identifying permeable pathways that enable geothermal fluids to circulate and produce energy. Specifically, we developed and implemented innovative machine learning methodologies to enhance geothermal exploration. Large-scale faults were detected using a Convolutional Neural Network (CNN), while small-scale fractures were characterized using a novel Double-Beam Neural Network (DBNN). These tools have proven both technically effective and cost-efficient by reducing reliance on expensive exploratory drilling. Through collaboration with our geothermal industry partner, this research has significantly advanced techniques for identifying hidden geothermal systems and extending the productive lifespan of existing geothermal fields. We applied our methods to two geothermal fields—Soda Lake (Nevada) and Lightning Dock (New Mexico)—to identify shallow steam-charged fracture zones and characterize deep faults at depths of 1.5-2 km. The steam zone identified at the Soda Lake geothermal field showed excellent agreement with prior drilling data, validating the effectiveness of our approaches. In addition, the analysis revealed three new prospective drilling targets for further development and verification. The outcomes of this project improve our scientific understanding of geothermal reservoir behavior, enhance exploration efficiency, extend the economic life of existing geothermal plants. Ultimately, these advancements contribute to GTO’s goal of achieving more sustainable, affordable, and data-driven geothermal energy development across the United States.

15 GEOTHERMAL ENERGY↗

Well-Controlled Hydro-thermo-chemical Stimulation of a Coso Geothermal Well to Enhance Productivity

The goal of this project was to demonstrate controlled fracturing and permeability increase to enhance resource production and power generation at Coso. It was planned to stimulate at least two intervals in the Well 73-18RD with zonal isolation using advanced high-temperature packer technology. The anticipation was that a hydraulic fracture would form in one zone and in the second zone natural fractures would reopen, shear and extend, thereby increasing fluid communication in the field. In both cases, the intent was to carefully design the stimulation to take advantage of the relative roles of hydraulic pressure, thermal stress and chemical (“hydro-thermo-chemical”) dissolution. The project had planned to use the Welltec Annular Barrier (WAB) system developed as part of another DOE-funded project which was on schedule to be ready for field deployment. If successful, the Coso Operating Company would incorporate this approach into their ongoing well maintenance in other wells increasing the total MWs produced. The project had multiple tasks and the objective of Task 1 was to develop a workover plan, to build a geologic model, to design a stimulation program, and to work with the operating and services companies to develop a detailed budget and obtain the necessary permits and sundry notices. In Phase 2 the workover plan would be implemented, the designed stimulation would be carried out, and subsequently the well would logged and tested to measure the impact of stimulation. Per sponsor request, a Pre-Task 1 activity with a Go/No Go decision point was scheduled in Task 1.

02 PETROLEUM↗

Accessing the Subsurface Biosphere Within Rocks Undergoing Active Low-Temperature Serpentinization in the Samail Ophiolite (Oman Drilling Project)

The Oman Drilling Project established an “Active Alteration” multi-borehole observatory in peridotites undergoing low-temperature serpentinization in the Samail Ophiolite. The highly serpentinized rocks are in contact with strongly reducing fluids. Distinct hydrological regimes, governed by differences in rock porosity and fracture density, give rise to steep redox (Eh +200 to –750 mV) and pH (pH range 8.5–11.2) gradients within the 300–400 m deep boreholes. The serpentinites and fluids host an active subsurface ecosystem. Microbial cell abundances in serpentinite vary at least six orders of magnitude, from ≤3.5 × 10 1 to 2.9 × 10 7 cells/g. Low levels of biological sulfate reduction (2–1,000 fmol/cm 3 /day) can be detected in rock cores, particularly in rocks in contact with reduced groundwaters with pH < 10.5. Thermodesulfovibrio is the predominant sulfate reducer identified via metagenomic sequencing of adjacent groundwater communities. Here, we infer that transport and reaction of microbially generated sulfide with the serpentine and brucite assemblages gives rise to optical darkening and sulfide overprinting, including the formation of tochilinite-vallerite group minerals, potentially serving as an indicator that this system is inhabited by microbial life. Olivine mesh-cores replaced with ferroan brucite and minor awaruite, abundant veins containing hydroandradite garnet and polyhedral serpentine, and late-stage carbonate veins are suggested as targets for future spatially resolved life-detection investigations. The high-quality whole-round core samples that have been preserved can be further probed to define how life distributes itself and functions within a system where chemical disequilibria are sustained by low-temperature water/rock interaction, and how biosignatures of in situ microbial activity are generated.

58 GEOSCIENCES↗

Long-term permeability evolution of shale seal rocks with argon and scCO 2

In this work, we conducted a suite of experiments to evaluate the long-term permeability evolution of shale rocks under constant effective stress, before, during, and after interaction with supercritical carbon dioxide (scCO 2 ). To do so, we measured the time-dependent evolution of argon permeabilities for the pre- and post-reacted samples (i.e., before and after long-term interaction with scCO 2 ). In addition, we obtained permeability evolution during long-term interaction with scCO 2 . The samples showed either relatively constant permeabilities or a moderate decrease during pre reaction long-term argon tests. The permeability evolution during long-term CO 2 tests showed continuous increase, continuous decrease, or cycles of increase/decrease in permeability. The long term response of the samples to CO 2 included phenomena such as 20 (i) salt precipitation, (ii) swelling-induced cracks, and (iii) carbonate dissolution. While it is obvious that salt precipitation and swelling-induced cracks decrease and increase the permeability, respectively, the sample response to carbonate dissolution proved to be more complex and may increase or decrease the permeability. The permeability evolution during post-reaction long-term argon injection is also affected by the contribution of each of these three phenomena, during long-term interaction with CO 2 . We observe increase, decrease, and constant permeability evolution during post-reaction argon tests. Our experiments reveal that the initial permeability of the samples plays a significant role on the long-term permeability response of shales in the presence argon and scCO 2 fluids. This study shows, when shales are hydraulically fractured with CO 2 , their initial permeability has a more significant role than their permeability evolution over time.

58 GEOSCIENCES↗

Coupling Subsurface and Above-Surface Models for Optimizing the Design of Borefields and District Heating and Cooling Systems

Accurate dynamic energy simulation is important for the design and sizing of district heating and cooling systems with geothermal heat exchange for seasonal energy storage. Current modeling approaches in building and district energy simulation tools typically consider heat conduction through the ground between boreholes without flowing groundwater. While detailed simulation tools for subsurface heat and mass transfer exist, these fall short in simulating above-surface energy systems. To support the design and operation of such systems, the study developed a coupled model including a software package for building and district energy simulation, and software for detailed heat and mass transfer in the subsurface. For the first, it uses 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, it uses the TOUGH simulator. The TOUGH family of codes 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. The study validated the coupled modeling approach by comparing the simulation results with one from the g-function based ground response model. It then looked into effects when the water table and the regional groundwater flow are considered in the ground, from the perspective of heat exchange between borehole and ground, and the electrical consumption of the district heating and cooling systems. To access the simulation models, please find the links in the submission: -- For coupled approach validation: see model Buildings.Fluid.Geothermal.Borefields.Examples.BorefieldsWithTough and Buildings.Examples.DistrictReservoirNetworks.Examples.Reservoir3Variable_TOUGH from the "Modelica Building Library" resource, branch issue1495_tough_interface, commit a2667c0. -- For the study of the effect of water table: see model Buildings.Examples.DistrictReservoirNetworks.Examples.Reservoir3Variable_TOUGH from he "Modelica Building Library" resource, branch issue1495_tough_interface_moreIO, commit 760de49. -- For the study of the effect of regional groundwater flow: see Buildings.Examples.DistrictReservoirNetworks.Examples.Reservoir3Variable_TOUGH from he "Modelica Building Library" resource, branch issue1495_tough_interface_moreIO_3D, commit c2a2d2a. The coupling interface script "GrounResponse.py" can be found from the above links in the folder Buildings/Resources/Python-Sources. Also, the needed files for TOUGH simulation are in the folder Buildings/Resources/Python-Sources/ToughFiles that can be accessed through the above links. A brief description of these files is given below; detailed specifications for the first three files may be found in the TOUGH3 Users Guide (Jung et al., 2018) https://tough.lbl.gov/documentation/tough-manuals/. (1) INCON - initial conditions for each grid block (2) INFILE - main input file with material properties and control parameters (3) MESH - description of the computational grid (4) readsave - Modelica/TOUGH interface program: read the final output of TOUGH simulation after TOUGH time step and prepare for transfer to Modelica for next Modelica time step (5) readsave.inp - input parameters for program readsave (6) writeincon - Modelica/TOUGH interface program: write the output of Modelica after Modelica time step and prepare for transfer to TOUGH as initial conditions for the next TOUGH step (7) writeincon.inp - input parameters for program writeincon

15 GEOTHERMAL ENERGY↗