Search NASA⌕ Search

SEARCH · Search NASA

Results for “permeability evolution”

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 37 records · Page 2

Coupled Thermo-Hydrological-Mechanical-Chemical Behavior of Anisotropic Granite for Geologic Disposal of High-Level Radioactive Waste: A Core-Scale Laboratory Investigation

The coupled thermo-hydrological-mechanical-chemical (THMC) behavior of rock within an Excavation Damaged Zone (EDZ) is critical for the safety and long-term performance of a geological repository for high-level radioactive wastes. While many laboratory experiments have been conducted to investigate EDZ rocks, the flow and deformation characteristics resulting from anisotropic rock textures and microcrack distribution under triaxial loading and elevated temperatures remain poorly understood. Particularly, cracks at various scales serve as fast paths for fluid flow and solute transport and present as focal points of mechanical weakness, which complicate the coupled THMC processes in anisotropic EDZ rocks and challenge modeling predictions. Here, in this study, a series of core-scale experiments was conducted on three granite samples under repository-relevant conditions. These rock samples were obtained from the Grimsel Underground Research Laboratory (URL), featured by anisotropic minerals (represented by bedding layers) and microcrack distributions and coarse cm-scale grain sizes. During the experiments, samples were subjected to an elevated temperature at 90 °C and different triaxial loading conditions either by radial (normal to bedding layers) or axial (parallel to bedding layers) compaction. Water was injected into the samples, and the rock permeability evolutions and effluent water chemistry were monitored closely. For intact samples, thermal expansion of minerals at 90 °C resulted in a large, 75% irreversible permeability reduction and rock strengthening under radial compaction, while thermal impact was limited to a 15% permeability reduction under axial compaction. In contrast, for a sample containing open cracks, the growth of fractures during the experiment resulted in an abrupt permeability increase and fast failure at 90 °C. The effluent water chemistry indicates much more considerable mineral dissolution from large shear sliding than that in rocks dominated by mechanical compaction. These results helped better understand the coupled THMC processes in anisotropic rocks containing cracks, evaluate the behaviors of EDZ rocks, and predict the long-term evolution of EDZ for the performance of the repositories.

Coupled THMC processes↗

Fluid migration in low-permeability faults driven by decoupling of fault slip and opening

Understanding the response of faults to the injection of high-pressure fluids is important for several subsurface applications, for example, geologic carbon sequestration or energy storage. Lab-based experiments suggest that fluid injection can activate fault slip and that this slip can lead to increased fluid transmission along low-permeability faults. Here we present in situ observations from a cross-borehole fluid-injection experiment in a low-permeability shale-bearing fault, which show fault displacement occurring before fluid-pressure build-up. Comparing these observations with numerical models with differing permeability evolution histories, we find that the observed variation in fluid pressure is best explained by a change in permeability only after the fault fails and slips beyond the pressurized area. Once fluid migration occurs along the fault as a result of slip-induced permeability increase, the fault experiences further opening due to a decrease in the effective normal stress. We then suggest that decoupling of fault slip and opening, leading to a rapid increase in fluid pressurization following the initial fault slip, could be an efficient driver for fluid migration in low-permeability faults.

58 GEOSCIENCES↗

Fault Slip and Fluid Flow: Seismic Source Analysis to Assess Role of Multiple Slip Patches in Fault Permeability

The relationship between fault reactivation, microearthquakes (MEQs), and permeability evolution during fluid injection plays a critical role in energy harvesting and waste disposal. Recent studies have demonstrated the possibility of predicting fault permeability using cumulative seismic moments of MEQs quantitatively. To understand the underlying physical processes, we conduct fault reactivation experiments using Utah FORGE granitoid and analyze acoustic emission (AE) signals generated during stepwise increases in fluid injection pressure. Frequency analysis of thousands of calibrated AE signals reveals that fault reactivation produces multiple AE source patches with millimeter-scale radii—smaller than the sample fault radius. The cumulative area of the reactivated patches covers the fault multiple times over (∼10x–50x area) for each pressure step. These findings provide mechanistic insight that measured permeability enhancement is not driven by a single large slip event, but by the sequential and interacting activation of multiple slip patches that create a continuous flow pathway.

Nurshal, M. E. M. [Pennsylvania State University, ↗

Rapid Mineral Precipitation During Shear Fracturing of Carbonate-Rich Shales

Target subsurface reservoirs for emerging low-carbon energy technologies and geologic carbon sequestration typically have low permeability and thus rely heavily on fluid transport through natural and induced fracture networks. Sustainable development of these systems requires deeper understanding of how geochemically mediated deformation impacts fracture microstructure and permeability evolution, particularly with respect to geochemical reactions between pore fluids and the host rock. Here, a series of triaxial direct shear experiments was designed to evaluate how fractures generated at subsurface conditions respond to penetration of reactive fluids with a focus on the role of mineral precipitation. Calcite-rich shale cores were directly sheared under 3.5 MPa confining pressure using BaCl 2 -rich solutions as a working fluid. Experiments were conducted within an X-ray computed tomography (xCT) scanner to capture 4-D evolution of fracture geometry and precipitate growth. Three shear tests evidenced nonuniform precipitation of barium carbonates (BaCO 3 ) along through-going fractures, where the extent of precipitation increased with increasing calcite content. Precipitates were strongly localized within fracture networks due to mineral, geochemical, and structural heterogeneities and generally concentrated in smaller apertures where rock:water ratios were highest. The combination of elevated fluid saturation and reactive surface area created in freshly activated fractures drove near-immediate mineral precipitation that led to an 80% permeability reduction and significant flow obstruction in the most reactive core. While most previous studies have focused on mixing-induced precipitation, this work demonstrates that fluid–rock interactions can trigger precipitation-induced permeability alterations that can initiate or mitigate risks associated with subsurface energy systems.

58 GEOSCIENCES↗

Feasibility of using reactive silicate particles with temperature-responsive coatings to enhance the security of geologic carbon storage

The large-scale deployment of geologic carbon storage will require minimizing the risk of CO 2 leakage. Here we report on a novel strategy to seal leakage pathways by injecting functionalized reactive calcium silicate (CaSiO 3 ) particles. The particles consist of a mineral silicate core within a temperature-sensitive Poly(N,N-dimethylaminoethyl methacrylate) coating that limits contact between the core and the surrounding CO 2 (aq) above a switching temperature making them relatively unreactive deep in the subsurface. If the CO 2 leaks to shallower depths and cooler temperatures, the coatings would extend and enable the silicate core to react with CO 2(aq) to form solid carbonate precipitates that reduce the permeability of the fracture or flow path. A combined laboratory- and simulation-based investigation was conducted to explore key early-stage research questions related to the feasibility of this leakage management strategy. Sand columns were injected with either coated or uncoated CaSiO 3 particles and the permeability was measured before and after exposure to CO 2 . Experiments were conducted at 35 °C and 70 °C, a range that spans the switching temperature for the coatings. Results show that at the cooler temperature the coating effectively enabled reaction and permeability reduction. At higher temperatures, the coated particles limited mineral carbonation and the permeability decrease was minimal compared to columns with uncoated particles. Here, the morphology and location of the precipitates provides insight about permeability evolution. Model simulations of 2D radial particle transport from an injection point provided insight into the effects of particle size, coating properties, and concentrations of particles in several representative rock formations.

54 ENVIRONMENTAL SCIENCES↗

Initial Simulations of Empty Room Collapse and Reconsolidation at the Waste Isolation Pilot Plant

The Waste Isolation Pilot Plant (WIPP) is a geologic repository for defense-related nuclear waste. If left undisturbed, the virtually impermeable rock salt surrounding the repository will isolate the nuclear waste from the biosphere. If humans accidentally intrude into the repository in the future, then the likelihood of a radionuclide release to the biosphere will depend significantly on the porosity and permeability of the repository itself. Room ceilings and walls at the WIPP tend to collapse over time, causing rubble piles to form on floors of empty rooms. The surrounding rock formation will gradually compact these rubble piles until they eventually become solid salt, but the length of time for a rubble pile to reach a certain porosity and permeability is unknown. This report details the first efforts to build models to predict the porosity and permeability evolution of an empty room as it closes. Conventional geomechanical numerical methods would struggle to model empty room collapse and rubble pile consolidation, so three different meshless methods, the Immersed Isogeometric Analysis Meshfree, Reproducing Kernel Particle Method (RKPM), and the Conformal Reproducing Kernel method, were assessed. First, the meshless methods and the finite element method each simulated gradual room closure, without ceiling or wall collapse. All three methods produced equivalent room closure predictions with comparable computational speed. Second, the Immersed Isogeometric Analysis Meshfree method and RKPM simulated two-dimensional empty room collapse and rubble pile consolidation. Both methods successfully simulated large viscoplastic deformations, fracture, and rubble pile rearrangement to produce qualitatively realistic results. In addition to geomechanical simulations, the flow channels in damaged salt and crushed salt were measured using micro-computed tomography, and input into a computational fluid dynamics simulation to predict the salt's permeability. Although room for improvement exists, the current simulation approaches appear promising.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Computational thermomechanics for crystalline rock. Part II: Chemo-damage-plasticity and healing in strongly anisotropic polycrystals

We present a thermal–mechanical–chemical-phase field model that captures the multi-physical coupling effects of precipitation creeping, crystal plasticity, anisotropic fracture, and crack healing in polycrystalline rock at various temperature and strain-rate regimes. This model is solved via a fast Fourier transfer solver with an operator-split algorithm to update displacement, temperature and phase field, and chemical concentration incrementally. In nuclear waste disposal in salt formation, brine inside the crystal salt may migrate along the grain boundary and cracks due to the gradient of interfacial energy and pressure. This migration has a significant implication on the permeability evolution, creep deformation, and crack healing within rock salt but is difficult to incorporate implicitly via effective medium theories compared with computational homogenization. As such, we introduce a thermodynamic framework and a corresponding computational implementation that explicitly captures the brine diffusion along the grain boundary and crack at the grain scale. Meanwhile, the anisotropic fracture and healing are captured via a high-order phase field that represents the regularized crack region in which a newly derived non-monotonic driving force is used to capture the fracture and healing due to the solution–precipitation. Numerical examples are presented to demonstrate the capacity of the thermodynamic framework to capture the multiphysics material behaviors of rock salt.

42 ENGINEERING↗

Impact of cement composition, brine concentration, diffusion rate, reaction rate and boundary condition on self-sealing predictions for cement-CO 2 systems

Geological CO 2 storage (GCS) plays an important role in curbing CO 2 emissions by reducing the carbon footprint of difficult to decarbonize operations and achieve negative CO 2 emissions through activities like Bioenergy with Carbon Capture and Storage (BECCS) and Direct Air Carbon Capture and Storage (DACCS). Leakage of CO 2 through wells is an important concern when it comes to deployment of large-scale GCS. Existing wells in sites that are otherwise suitable for GCS can act as conduits for stored CO 2 to escape the reservoir. There is broad consensus that the main risk of leakage through wellbores is via fractures/damaged pathways. The results from several studies evaluating the permeability evolution of cement fractures in wells upon leakage of CO 2 agree that smaller fracture apertures, slower brine velocities and higher brine residence times promote self-sealing of fractures by mineral precipitation. Quantitatively, however, the differences in sealing conditions are significant and are typically attributed to differences in experimental conditions or model assumptions. Here we examine the sensitivity of our model, describing CO 2 leakage through wellbores, to cement composition, brine concentration, diffusion rates, and reaction rates. We also evaluate the impact of the boundary condition to allow comparisons between observations from experiments performed at constant flow rate and model predictions made at constant pressure conditions. Our results show that diffusion and reactions rates have the most impact on the self-sealing criteria for cement-CO 2 systems. In addition, conditions associated with self-sealing of fractures at constant flow rate require longer fractures, smaller fracture apertures, and slower velocities than under constant pressure.

58 GEOSCIENCES↗

Reservoir-scale model of geologic hydrogen production from serpentinization: Cyclic injection in a dual-permeability fracture-matrix system

Geologic hydrogen (GeoH 2 ) from serpentinization is a promising low-carbon resource, but its reservoir-scale behavior remains poorly understood. We develop a dual-permeability reactive transport model for an injector–producer well pair in ultramafic rock that couples multiphase flow, heat transfer, geochemistry, and porosity–permeability evolution. Here, the model is calibrated to olivine flow-through experiments and upscaled to two-year long simulations with continuous injection and cyclic injection with shut-in-to-injection ratios (SIR = 1, 0.5, 0.1). Olivine reacts along high-flux pathways to form lizardite and magnetite, increasing pH; H 2 (aq) and H 2 (g) peak early and then decline as exsolution and advective export outpace local generation. Continuous injection yields the highest cumulative H 2 but the lowest water-use efficiency (1.138 x 10 –6 mol/kgw). Cyclic injection increases this ratio to 1.35 x 10 –6 , 1.377 x 10 –6 , and 1.182 x 10 –6 mol/kgw for SIR = 1, 0.5, and 0.1, respectively; SIR = 0.5 provides a ~20% improvement over continuous injection and the best compromise for pilot design.

08 HYDROGEN↗

Wide-blocky veins explained by dependency of crystal growth rate on fracture surface type: Insights from phase-field modeling

Vein microstructures contain a wealth of information on coupled chemical and mechanical processes of fracturing, fluid transport, and crystal growth. Numerical simulations have been used for exploring the factors controlling the development of vein microstructures; however, they have not been quantitatively validated against natural veins. Here we combined phase-field modeling with microtextural analysis of previously unexplained wide-blocky calcite veins in natural limestone and of the fresh fracture surface in this limestone. Results show that the wide-blocky vein textures can only be reproduced if ~10%–20% of crystals grow faster than the rest. This fraction corresponds to the amount of transgranularly broken grains that were observed on the experimental fracture surfaces, which are dominantly intergranular. We hypothesize that transgranular fractures allow faster growth of vein minerals due to the lack of clay coatings and other nucleation discontinuities that are common along intergranular cracks. Our simulation results show remarkable similarity to the natural veins and reproduce the nonlinear relationship between vein crystal width and vein aperture. This allows accurate simulations of crystal growth processes and related permeability evolution in fractured rocks.

Geology↗

DECOVALEX-2023: Task G, SAFENET Final Report

DECOVALEX Task G deals with fracture mechanics at several scales using a combined approach of experimental work, related model development, benchmarking and experimental analysis. Fig. 0.1 provides a graphical abstract for Task G. The experimental basis for the related mechanical (M), hydro-mechanical (HM), and thermo-mechanical (TM) processes comes from the rock mechanics laboratories of Universities of Freiberg (TUBAF) and Edinburgh, as well as the Korea Institute of Civil Engineering and Building Technology (KICT). The experimental work in the rock laboratories is closely linked to the underground research laboratories (URLs) Reiche Zeche (Germany), KURT (Korea), and Mont Terri (Switzerland). The scientific key questions are related to fracture permeability evolution under THM conditions, anisotropy effects on fracturing processes, and thermal fracture slip, which are being addressed in specific steps of the task. A large variety of numerical methods have been developed and applied for experimental analysis, ranging from continuum to discontinuum approaches. A detailed comparison of mechanical and hydro-mechanical processes is given in the section 3 via the benchmarking exercises. As a result of DECOVALEX-2023 Task G we further improved our understanding and predictability of fracturing processes under THM conditions. This was based on robust numerical simulation methods and in-depth experimental analysis.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Prediction and Analysis of Utah FORGE Injection Activities using a Coupled Thermo-hydro-mechanical and Earthquake (THM+E) Modeling Workflow

A coupled thermo-hydro-mechanical (THM) numerical workflow that is capable of modeling seismic slip is critical for the successful development of enhanced geothermal systems (EGS). By integrating key physical processes, this workflow enables accurate simulation of temperature and pressure diffusions, stress changes, and induced seismicity. As a result, it serves as a vital tool for predicting induced seismicity and optimizing reservoir stimulation strategies. The Utah FORGE (Frontier Observatory for Research in Geothermal Energy) project, located near Milford, Utah, is a U.S. Department of Energy initiative aimed at advancing EGS technology. In April 2024, eight new stimulation stages (Stages 3R-10) were conducted in well 16A (injection well) subsequent to the first series of stimulation (Stages 1-3) performed in April, 2022. To monitor the induced seismicity, geophones were deployed in wells 58-32, 56-32, and 78B-32, while fiber optic cables were also installed in wells 16B, 78-32, and 78B-32 to collect microseismic data and detect frac hits Preliminary analyses of microseismic catalogs and fiber optic data suggest that the stimulated fractures in Stages 3R–6 closely align with that generated during Stage 3, indicating that the new stimulations were likely reactivating the previously stimulated fracture. To better understand the underlying process, a comprehensive modeling approach that can accurately capture thermal, hydrological, mechanical, and seismic responses is essential. In this work, we propose and utilize a coupled thermo-hydro-mechanical and earthquake (THM+E) simulation workflow to numerically investigate the stimulation activities on well 16A. The specific objective is to confirm whether the new stimulation stages (Stages 3R–6) reactivated fractures previously stimulated during Stage 3. For this purpose, we perform THM+E simulations individually for Stages 3, 3R, 4, and 5, incorporating the discrete fracture networks (DFNs) created by the plane-fitting technique based on the microseismic catalogs. The simulation workflow consists of two separate models: a THM model and an earthquake model, coupled in a one-way manner. Detailed descriptions of the workflow are provided in Section 3. Simulation results are presented in terms of injection pressure, permeability evolution, and predicted seismic catalogs, which are then compared with field data for further analyses. This report is structured as follows. In Section 2, we present detailed analyses of the field data and propose the hypothesis that the new stimulation stages (Stages 3R–6) were probably reactivating the previously stimulated fractures in Stage 3. In Section 3, we introduce the coupled THM+E workflow and the problem setup to validate our hypothesis, followed by the simulation results for each stage in Section 4. Meanwhile, discussions are included to analyze the model predictions and their comparison with field data. Lastly, we conclude the report and outline future plans in Section 5.

15 GEOTHERMAL ENERGY↗

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↗

Machine Learning Approaches to Predicting Induced Seismicity and Imaging Geothermal Reservoir Properties

This project developed machine learning (ML) methods, lab data sets, and field data to advance geothermal exploration and geothermal energy production. The work had three focus areas. One involved the development of ML methods to use microearthquakes (MEQs) for imaging geothermal reservoir properties and improving subsurface characterization – most importantly the evolution of permeability within the evolving reservoir. This part of the work included development of ML approaches for automated MEQ location, focal mechanism determination and identification of earthquake precursors. The second area focused on using MEQ signals generated by geothermal exploration and production to predict the relationship between fluid injection and seismicity. Here, we extended to reservoir scale our success in using ML to predict laboratory earthquakes and fault zone stress state. The third focus area was on lab experiments. Here, we developed new ML models for lab earthquake prediction and identification of precursors to failure to improve earthquake forecasting and early warning in geothermal settings. Major outcomes of our work include ML models that learn from MEQ signals during geothermal exploration and production to predict induced seismicity. MEQs occur naturally in connection with drilling and energy production. We developed ML methods to use the seismic waves from these events to characterize the elastic, hydraulic and poromechanical properties of reservoirs. Our work illuminated fracture geometry and the evolution of fracture permeability by incorporating seismic coda wave analysis and ML methods to relate fluid injection and seismicity. We significantly expanded laboratory earthquake prediction to include methods that use both passive measurements of microearthquakes within the lab fault zones and also active source acoustic measurements of fault zone elastic properties. These methods can now predict fault zone stress state, time to failure and the magnitude of lab earthquakes. Our work showed that repetitive stick- slip failure events during frictional sliding (the lab equivalent of earthquakes) are preceded by a cascade of micro-failure events that radiate energy in a manner that foretells unstable failure – manifest as laboratory MEQs. We documented a mapping between fracture properties and statistical attributes of elastic radiation. We extended existing works to geothermal reservoir scale and developed ML methods to determine reservoir permeability, fracture properties, and their evolution during geothermal energy production. An attractive feature of ML algorithms is their ability to handle big datasets and reveal patterns and correlations that may remain invisible to conventional analyses. Our work connected data from field, laboratory and intermediate scales to study permeability, stress, strength, fracture stiffness and geometry. At the field scale we used data from the Newberry Volcano field site, UtahFORGE, EGS Collab, and also the Bedretto underground research lab in Switzerland. These data sets are bridging the gap between the lab scale, theory, and reservoir scale. Our work produced plain language summaries to improve public understanding of DOE research. We also developed openly distributed ML and seismicity datasets for use by all researchers and we published connections between induced seismicity in geothermal areas and reservoir properties including permeability, fracture properties, and stress state. Our models are designed for the large data sets of induced seismicity typically associated with geothermal sites. We produced labeled event catalogs and used them on geothermal data to assess how ML can facilitate geothermal production and exploration. All datasets are available on the GDR Productivity: The project produced 32 publications in peer reviewed journals (two are in review). It supported the work of 6 PhD students, 40 conference presentations, 6 keynote talks at national meetings, and mentoring and professional development for 4 postdoctoral fellows.

15 GEOTHERMAL ENERGY↗

Evolution of Pore Structure and Permeability of Rocks Under Hydrothermal Conditions (Final Report)

The physical and transport properties of porous rocks can be altered by a variety of diagenetic, metamorphic, and tectonic processes, and the changes that result are of critical importance to such industrial applications as resource recovery, carbon dioxide sequestration, and waste isolation in geologic formations. The associated interrelationships between rocks, pore fluids and deformation are also key to understanding many natural processes, including dynamic metamorphism, fault mechanics, fault stability and pressure solution creep. In this project, we investigated the changes of permeability and pore geometry owing to inelastic deformation by solution-transfer, brittle fracturing and dislocation creep. In particular, we studied the coupling between pore fluids and deformation in fluid-filled quartz, calcite, mudstones and ultramafic rocks and examine the effects of coupled mechanical and chemical processes on the evolution of porosity and permeability under hydrothermal conditions. The investigations used a combination of techniques, that included triaxial laboratory experiments; uniaxial compressive loading with in situ observations of microstructure; permeability under hydrostatic compression; observations of microstructure using optical and electron microscopes, and micro-CT imaging; and numerical calculations. Laboratory experiments were designed to provide mechanical and transport data under conditions that isolate the particular mechanisms responsible for the changes. The data obtained were used to quantify changes in surface roughness, porosity, pore dimensions and their spatial fluctuations. The results of the experiments and data from image analyses were compared to the results of network, finite-difference and other numerical models to verify the validity of experimentally established relations between permeability and other rock properties. A bibliography for the period from 2015-2019 is appended below. New results obtained in the final year and the unfunded extension period are detailed below.

58 GEOSCIENCES↗

Utah FORGE 5-2557: Fluid and Temperature in Fracture Mechanics and Coupled THMC Processes - Workshop Presentation

This is a presentation on the Role of Fluid and Temperature in Fracture Mechanics and Coupled Thermo-Hydro-Mechanical-Chemical (THMC) Processes for Enhanced Geothermal Systems project by Purdue University, presented by Distinguished Professor of Physics & Astronomy, Laura J. Pyrak-Nolte. The project's objective was to develop and validate a macroscopic model that accounts for local deformation/frictional behavior, seismic/aseismic behavior, chemical reactions, and determine the adequacy of classic Coulomb failure vs. rate-and-state friction. This presentation was featured in the Utah FORGE R&D Annual Workshop on September 8, 2023. The workshop provided a valuable opportunity to explore the progress made in each of the 17 Research and Development projects funded under Solicitation 2020-1 which aim to enhance our understanding of the crucial factors influencing the development of Enhanced Geothermal Systems (EGS) reservoirs and resources.

15 GEOTHERMAL ENERGY↗