Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “subsurface storage”

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 235 records · Page 13
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Design, construction, and validation of an in-situ groundwater trace element analyzer with applications in carbon storage

Abstract It is estimated that carbon emissions should reach net-zero by 2050 to meet important climate targets. Carbon capture is likely necessary to reach these targets, requiring a long-term storage solution such as geological carbon sequestration. However, as with any subsurface activity, leakage can occur, potentially impacting groundwater quality near the storage site. Rapid detection is essential to mitigate damage to this resource. Since CO 2 will acidify groundwater, the concentrations of acid soluble minerals and associated cations will increase. Thus, an in-situ, real-time element analysis system based on laser-induced breakdown spectroscopy (LIBS) is under development to monitor these elements. The system splits the traditional LIBS system into a miniature, all-optical sensor head built around a passively Q-switch laser fiber coupled to a control unit. Previous work has validated the LIBS technique for use at high pressure as well as the split system design. In this work, a fieldable prototype sensor is developed and tested in an onsite monitoring well where trace elements concentrations (approx. 0–3 ppm) were tracked over 20 days. These concentrations varied in response to local rainfall, diluting with increased rain, demonstrating the ability of a LIBS-based sensor to track trace elements under real-world conditions.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Basic Energy Sciences Roundtable: Foundational Science for Carbon Dioxide Removal Technologies (Brochure)

Basic research to identify and understand the fundamental principles governing carbon dioxide removal processes of capture, conversion, and storage is essential for achieving zero carbon emissions by 2050. In March 2022, the Office of Basic Energy Sciences (BES)—in coordination with the US Department of Energy (DOE) technology Offices of Energy Efficiency and Renewable Energy, Fossil Energy and Carbon Management, and Nuclear Energy—held a roundtable titled, “Foundational Science for Carbon Dioxide Removal Technologies,” to discuss the scientific and technical barriers for CO 2 capture, conversion, and storage. Five priority research opportunities were identified to address these scientific and technical challenges and accelerate progress toward the realization of zero carbon emissions. The PROs are as follows: 1) Master Interfacial Processes of CO2 Transport and Reactivity Across Multiple Length and Time Scales; 2) Create Materials that Simultaneously Exhibit Multiple Properties for CO 2 Capture and Release or Conversion; 3) Discover Unconventional Pathways for Energy-Efficient CO 2 Capture, Release, and Conversion; 4) Control Multiphase Interactions Required for CO 2 Conversion into Minerals and Materials; and 5) Achieve Predictive Understanding of Coupled Processes in Complex Subsurface Geologic Systems for Secure Carbon Storage.

36 MATERIALS SCIENCE↗
DOE OSTI2022

A Project Lifetime Approach to the Management of Induced Seismicity Risk at Geologic Carbon Storage Sites

The geologic storage of carbon dioxide (CO 2 ) is one method that can help reduce atmospheric CO 2 by sequestering it into the subsurface. Large-scale deployment of geologic carbon storage, however, may be accompanied by induced seismicity. We present a project lifetime approach to address the induced seismicity risk at these geologic storage sites. This approach encompasses both technical and nontechnical stakeholder issues related to induced seismicity and spans the time period from the initial consideration phase to postclosure. These recommendations are envisioned to serve as general guidelines, setting expectations for operators, regulators, and the public. They contain a set of seven actionable focus areas, the purpose of which are to deal proactively with induced seismicity issues. Although each geologic carbon storage site will be unique and will require a custom approach, these general best practice recommendations can be used as a starting point to any site-specific plan for how to systematically evaluate, communicate about, and mitigate induced seismicity at a particular reservoir.

58 GEOSCIENCES↗
DOE OSTI2024

Impacts of Topography‐Driven Water Redistribution on Terrestrial Water Storage Change in California Through Ecosystem Responses

Abstract Lateral subsurface flow plays an essential role in sustaining the terrestrial ecosystem, but it is not explicitly represented in most Earth System Models. In this study, we implemented an explicit lateral saturated flow model into the E3SM land model (ELM). The model explicitly describes lateral flow in the saturated zone by representing, for each model grid, an idealized hillslope consisting of five hydrologically connected soil columns. We conducted three model experiments driven by 0.125° atmospheric forcing data during 1980–2015 over California using models of the default ELM, a modified version of ELM to enhance infiltration, and the model with the lateral saturated flow model. The simulated runoff, evapotranspiration, and terrestrial water storage anomaly (TWSA) from the three simulations were evaluated against available observations, and the model explicitly representing lateral flow performs best. The new model produces greater gridcell‐averaged evapotranspiration especially over the mountainous regions with moderate relief and seasonally dry climates. Most importantly, it improves the modeled seasonal variations, interannual variabilities, and the recent decadal decline of TWSA. Many of these improvements can be attributed to the enhanced ecosystem resilience to droughts as demonstrated by transpiration increases caused by lateral flow. Model sensitivity experiments suggest that subsurface runoff is most sensitive to the ratio between horizontal and vertical saturated hydraulic conductivity, followed by hillslope planforms (convergent, divergent, and uniform), number of columns, and lower boundary conditions. Future work should effectively characterize hillslopes in global models and explore the long‐term influences of lateral water movement on modeled biogeochemical cycle.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

The Kimberlina synthetic multiphysics dataset for CO 2 monitoring investigations

Abstract We present a synthetic multi‐scale, multi‐physics dataset constructed from the Kimberlina 1.2 CO 2 reservoir model based on a potential CO 2 storage site in the Southern San Joaquin Basin of California. Among 300 models, one selected reservoir‐simulation scenario produces hydrologic‐state models at the onset and after 20 years of CO 2 injection. Subsequently, these models were transformed into geophysical properties, including P‐ and S‐wave seismic velocities, saturated density where the saturating fluid can be a combination of brine and supercritical CO 2 , and electrical resistivity using established empirical petrophysical relationships. From these 3D distributions of geophysical properties, we have generated synthetic time‐lapse seismic, gravity and electromagnetic responses with acquisition geometries that mimic realistic monitoring surveys and are achievable in actual field situations. We have also created a series of synthetic well logs of CO 2 saturation, acoustic velocity, density and induction resistivity in the injection well and three monitoring wells. These were constructed by combining the low‐frequency trend of the geophysical models with the high‐frequency variations of actual well logs collected at the potential storage site. In addition, to better calibrate our datasets, measurements of permeability and pore connectivity have been made on cores of Vedder Sandstone, which forms the primary reservoir unit. These measurements provide the range of scales in the otherwise synthetic dataset to be as close to a real‐world situation as possible. This dataset consisting of the reservoir models, geophysical models, simulated time‐lapse geophysical responses and well logs forms a multi‐scale, multi‐physics testbed for designing and testing geophysical CO 2 monitoring systems as well as for imaging and characterization algorithms. The suite of numerical models and data have been made publicly available for downloading on the National Energy Technology Laboratory's (NETL) Energy Data Exchange (EDX) website.

58 GEOSCIENCES↗
DOE OSTI2020

Impacts on mechanical strength of chemical reactions induced by hydrous supercritical CO2 in Boise Sandstone

Geomechanics experiments were used to assess mechanical alteration of Boise Sandstone promoted by reactions with supercritical carbon dioxide (scCO 2 ) and water vapor. During geologic carbon storage, scCO 2 is injected into subsurface reservoirs, forming buoyant plumes. At brine-plume interfaces, scCO 2 can dissolve into native brines, and water from brines can partition into scCO 2 , forming hydrous scCO 2 . This study investigates the effect of hydrous scCO 2 on the strength of Boise Sandstone. Samples are first exposed to recirculating hydrous scCO 2 for 24 h at 70 °C and 13.8 MPa scCO 2 pressure. Samples are reacted with scCO 2 with added water contents up to 500 mL. After scCO 2 exposure, samples are deformed at room temperature under confining pressures of 3.4, 6.9, and 10.3 MPa. The results demonstrate that hydrous scCO 2 induces chemical reactions in Boise Sandstone, with ions migrating from the solid into the hydrous scCO 2 phase. At the longer time-scales, these reactions could lead to mechanical weakening in the samples; however, on the scale of our experiments, the strength changes are within sample variability. Because the solubility of water in scCO 2 is extremely low (0.008 mol H 2 O per 1 mol CO 2 ), the mineral dissolution of Boise Sandstone was under 0.002 wt.%. Additionally, mineral grains and pore throats in Boise Sandstone are cemented with quartz, which is not susceptible to dissolution at these conditions. Our results indicate that humidity in scCO 2 plumes is unlikely to sustain chemical reactions and induce long term strength changes in quartz cemented sandstones due to resistant mineralogies and low water solubility.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2026

A robust framework for frictional fault contact in geological formations using a stabilized augmented Lagrangian approach

Numerical simulations are essential to evaluate the performance and safety of engineered subsurface systems such as geological carbon storage sites, enhanced geothermal fields, and oil and gas reservoirs. A key challenge lies in accurately modeling the frictional contact behavior along fault surfaces. This problem involves inequality constraints that arise from the physics of frictional slip, requiring specialized numerical methods to handle the resulting highly nonlinear and path-dependent behavior. Here, in this work, we address this challenge using an Augmented Lagrangian Method (ALM) implemented via the Uzawa algorithm. The formulation employs mixed finite element spaces, combining low-order piecewise linear displacements within the 3D domain cells with piecewise constant tractions defined on the fault surfaces. Furthermore, to ensure stability and satisfy the inf-sup condition, the discrete displacement space is enriched with face bubble functions on both sides of the contact interfaces. This approach offers several advantages over other stabilization techniques that rely on additional terms, and it integrates naturally in the Uzawa framework.

58 GEOSCIENCES↗
DOE OSTI2022

Time-lapse seismic data inversion for estimating reservoir parameters using deep learning

Geologic carbon sequestration involves the injection of captured carbon dioxide ([Formula: see text]) into subsurface formations for long-term storage. The movement and fate of the injected [Formula: see text] plume is of great concern to regulators because monitoring helps to identify potential leakage zones and determines the possibility of safe long-term storage. To address this concern, we design a deep-learning framework for [Formula: see text] saturation monitoring to determine the geologic controls on the storage of the injected [Formula: see text]. We use different combinations of porosities and permeabilities for a given reservoir to generate saturation and velocity models. We train the deep-learning model with a few time-lapse seismic images and their corresponding changes in saturation values for a particular [Formula: see text] injection site. The deep-learning model learns the mapping from the change in the time-lapse seismic response to the change in [Formula: see text] saturation during the training phase. We then apply the trained model to data sets comprising different time-lapse seismic image slices (corresponding to different time instances) generated using different porosity and permeability distributions that are not part of the training to estimate the [Formula: see text] saturation values along with the plume extent. Our algorithm provides a deep-learning assisted framework for the direct estimation of [Formula: see text] saturation values and plume migration in heterogeneous formations using the time-lapse seismic data. Our method improves the efficiency of time-lapse inversion by streamlining the large number of intermediate steps in the conventional time-lapse inversion workflow. This method also helps to incorporate the geologic uncertainty for a given reservoir by accounting for the statistical distribution of porosity and permeability during the training phase. Tests on different examples verify the effectiveness of our approach.

Geochemistry & Geophysics↗
DOE OSTI2022

Task 5 Field Scale Stress Modeling Scripts and Data

Project data and Matlab scripts for Task 5 of the project “A Non-Invasive Approach for Elucidating the Spatial Distribution of In Situ Stress in Deep Subsurface Geologic Formations Considered for CO2 Storage”. Project includes the static earth models, grids and other data used for simulation and calibration of the two case studies addressed in the project, referred to herein as the “FutureGen2” model (Illinois) and the “Perch” model (Otsego county, Michigan). Large Dataset, please contact EDXSupport@netl.doe.gov

FutureGen2↗
DOE OSTI2026

Illinois Storage Corridor - CarbonSAFE Phase III: Policy, Regulatory, Legal and Permitting Characteristics; Subtask 5.5

The Illinois Storage Corridor (ISC) project aims to advance the commercial development and implementation of carbon capture and storage (CCS) technologies within a region in Illinois of suitable geology for carbon dioxide (CO 2 ) storage in deep saline aquifers. The project partners have annual emissions greater than 6.5 million tonnes per year; storage hubs are being explored at two sites—one near the One Earth Energy facility in east-central Illinois, and one near the Prairie State Generating Company campus in southwest Illinois. To implement the technology, legal, policy, and economic considerations must be explored. The United States Environmental Protection Agency (US EPA) administers the Underground Injection Control (UIC) program and is responsible for issuing Class VI permits to construct and operate CO 2 injection wells, i.e., a UIC Class VI well permit is required to inject CO 2 into the subsurface for geologic sequestration. The largest consideration of the Class VI well requirements is to protect underground sources of drinking water (USDWs). Additionally discussed are considerations for Class VI permits relating to public engagement. Besides permitting, property rights to storage sites, subsurface pore spaces, and areas for pipeline transportation must be secured. Pore space rights is still a novel concept being explored and not yet addressed by the Illinois legislature. It is believed that surface property rights are required for the subsurface pore space below, so long as there is not a separated mineral estate in the subsurface. Illinois legislature has addressed securing rights-of-way for CO 2 transportation, allowing for easements and the exercise of eminent domain to secure such rights. Economically, the incentives for CCS are ever-expanding. Recently, the federal government broadened the availability and increased the dollar-amounts for the § 45Q tax credits for geological storage of carbon oxides. Congress has also authorized the Advanced Industrial Facilities Development Program which allocates billions of dollars in funds to installing technology at industrial facilities to reduce greenhouse gas emissions. On the state level, Illinois has had its own incentives for CCS projects since 2009 and expanded its emission-related goals again in 2021. The law enacted in 2021 specifically creates a commission to explore implementing CCS at Prairie State Generating Company, a partner on the Illinois Storage Corridor project.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Exploring the Feasibility of a Carbon Dioxide Storage Hub in Western North Dakota

The University of North Dakota Energy & Environmental Research Center (EERC) and project partner ONEOK, Inc. (ONEOK) are investigating the feasibility of establishing a CO2 (carbon dioxide) storage hub in western North Dakota—in the heart of the Williston Basin. The conceived Roughrider Carbon Storage Hub would store CO2 captured from six gas-processing plants owned and operated by project partner ONEOK and a planned gas-to-liquids plant. This 2-year U.S. Department of Energy-sponsored Carbon Storage Assurance Facility Enterprise (CarbonSAFE) Phase II feasibility study is evaluating the aggregation of the CO2 captured from these seven sources for injection into stacked geologic storage complexes. The proposed hub includes several aspects that make it a highly qualified candidate for a feasibility study with a notably reduced project risk profile. These include 1) a project partner (ONEOK) with a committed goal to reduce greenhouse gas emissions; 2) prior subsurface data analysis supporting a potential stacked storage configuration with adequate CO2 storage resource; 3) commitment from local, regional, and state-level stakeholders; and 4) a state with U.S. Environmental Protection Agency underground injection control Class VI primacy. ONEOK’s assets in the Williston Basin provide significant environmental benefits by capturing and processing natural gas that may otherwise be flared or vented. Storing CO2 from these gas-processing facilities will reduce overall CO2 emissions in the basin while providing essential services to producers there and contributing to continued energy independence in the domestic markets.

03 NATURAL GAS↗
DOE OSTI2024

Editorial: Subsurface microbiology within hydrocarbon resources or stored gases

A Research Topic on the microbiology of hydrocarbon and gas storage reservoirs has far reaching industrial applications. In recent decades, there has been a growing interest in understanding microbial communities in subsurface energy reservoirs, such as coal, oil, and shale beds. This area of research has broadened to include gas storage reservoirs for hydrogen and CO 2 . Scientists are beginning to unravel the unexpected impact microorganisms have on these systems, through changing the fluid geochemistry, the gas content, and even the permeability. By recognizing the influence of these tiny organisms on our engineered environments, we can develop better risk assessments, target mitigation strategies, expand energy production, and refine operational guidance, ultimately contributing to a more sustainable energy future.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2020

Wellbore injectivity response to step-rate CO 2 injection: Coupled thermo-poro-elastic analysis in a vertically heterogeneous formation

Safe and permanent carbon geological storages require an elaborate analysis of geomechanical stability. Subsurface injection of CO 2 changes the local temperature and pore pressure and, further, alters the stress due to thermo-poro-elastic responses. A permanent CO 2 storage testing was conducted in the water leg of Cranfield reservoir in Mississippi, USA, where hosted CO 2 enhanced-oil recovery activities. During CO 2 injection, the injection rate was ramped up twice but the bottom-hole pressure did not increase with the imposed injection rates as expected. This unexpected field observation suggests the possibility of an open-mode fracture development at the injector. However, the injector response and potential fracture development have not been rigorously interpreted with detailed near-wellbore temperature and pore pressure change upon the CO 2 injection. In this study, we performed history matching of CO 2 injection using coupled thermo-poro-elastic reservoir simulation to examine the possibility of fracturing. We built a reservoir model including vertical heterogeneity in both petrophysical and geomechanical properties estimated from well-logging analysis and laboratory experiments. The simulation results show that CO 2 injection changes stresses and support the hypothesis of development of an open-mode fracture at the injector during the Cranfield test. The near-injector region exhibits a large temperature reduction up to 55 °C with ensuing effective horizontal stress reduction up to 9.1 MPa. However, a caprock integrity issue is unlikely because of the horizontal stress contrast within layers and low hydraulic communication with the injection zone. This work indicates that the injectant temperature should be considered in the design of high-rate CO 2 injector.

58 GEOSCIENCES↗
DOE OSTI2022

Quantitative Evaluation of CO 2 Storage Potential in the Offshore Atlantic Lower Cretaceous Strata, Southeastern United States

The geological storage of CO 2 in the Earth’s subsurface has the potential to significantly offset greenhouse gas emissions for safe, economical, and acceptable public use. Due to legal advantages and vast resource capacity, offshore CO 2 storage provides an attractive alternative to onshore options. Although offshore Lower Cretaceous reservoirs have a vast expected storage capacity, there is a limited quantitative assessment of the offshore storage resource in the southeastern United States. This work is part of the Southeast Offshore Storage Resource Assessment (SOSRA) project, which presents a high-quality potential geological repository for CO 2 in the Mid- and South Atlantic Planning Areas. This is the first comprehensive investigation and quantitative assessment of CO 2 storage potential for the Lower Cretaceous section of the outer continental shelf that includes the Southeast Georgia Embayment and most of the Blake Plateau. An interpretation of 200,000 km of legacy industrial 2D seismic reflection profiles and geophysical well logs (i.e., TRANSCO 1005-1-1, COST GE-1, and EXXON 564-1) were utilized to create structure and thickness maps for the potential reservoirs and seals. We identified and assessed three target reservoirs isolated by seals based on their effective porosity values. The CO 2 storage capacity of these reservoirs was theoretically calculated using the DOE-NETL equation for saline formations. The prospective storage resources are estimated between 450 and 4700 Mt of CO 2 , with an offshore geological efficiency factor of dolomite between 2% and 3.6% at the formation scale.

03 NATURAL GAS↗
DOE OSTI2023

Thermo-Hydrological Modeling of Thermal Energy Storage in a Depleted Oil Reservoir

Thermal energy storage in oil and gas reservoirs leverages the existing surface and subsurface infrastructure, which can pave the way for economic production of geothermal energy. Existing studies on geothermal energy storage are focused mostly on the use of aquifers with more homogeneous rock and fluid properties. Coupling of heat and fluid flow in a multiphase-multicomponent system, such as an oil reservoir, is imperative especially if existing oil field assets need to be repurposed as required for a sustainable energy transition. The objective is to model the subsurface thermo-hydrological processes associated with reservoir performance and operational sustainability. The model evaluates formation pressure and temperature within the reservoir and at the injection/production wells during multiple charge and discharge cycles. Hot water (approximately 200 degrees C) heated by Concentrating Solar Power (CSP) at high pressure is injected into the existing oil reservoir for storage and produced as thermal energy for power generation, which will be accompanied by enhanced oil recovery. To demonstrate the coupled fluid and heat flow during the injection/production cycle in the subsurface reservoir, TOUGH3 (developed by Berkeley Lab) is used to simulate the thermo-hydrological (TH) processes in a multiphase, multicomponent system. Two well geometries are considered within the reservoir grid: 1) a single-well huff-n-puff system (same well is used for injection and production), and 2) an isolated injection-production well doublet. Seasonal charge and discharge cycling are implemented based on the scheduling specified in the model input file. The model reports pressure, temperature, enthalpy, liquid fluxes, heat fluxes, pore velocities, and changes in porosity & permeability due to temperature and pressure variations during the cyclic Reservoir Thermal Energy Storage (RTES) operations. The results from the simulations can be used to optimize the operational parameters (such as well spacing and injection/production rates) and round-trip efficiency for surface power-plants coupled with thermal energy storage over time. They can also serve as important inputs for levelized cost of storage estimations. The research will help to design and integrate surface renewable energy sources, such as concentrating solar power (CSP), with RTES to help balance out power supply and demand on the grid.

CSP↗
DOE OSTI2026

The future of subsurface monitoring: AEC’s breakthroughs in CCS technology

Carbon capture and storage (CCS) has emerged as a key solution in the fight against climate change. However, for CCS to succeed, it is crucial to ensure that the sequestered CO2 stays safely trapped underground. The U.S. Department of Energy (DOE) has emphasized the need for advancements in subsurface monitoring, measurement, reporting, and verification. Aside from caprock integrity failure, the other primary failure points usually involve defective cement in the casing annulus of wellbores or plugged and abandoned wells. In addition, many energy producers (e.g., oil and gas, geothermal) and storage and disposal operators (e.g., H2 and water) must deal with the same issue. Poorly placed or degraded cement can create pathways for gas or fluid to escape from casing annuli and in plugged and abandoned or orphan wells, posing environmental risks. Yet, a reliable and cost-effective way to monitor cement and well integrity over multiple decades is still unavailable. Traditional geophysical methods like 4D seismic imaging and surface-based electromagnetic monitoring lack the resolution and accuracy for detecting these types of failures (Vasco et al., 2022; Fawad and Mondol, 2021). Wireline logging is expensive to run continuously and is obtrusive to the operation. While fiber optics can potentially be a solution, its bulkiness can significantly compromise the cement's integrity. To address these challenges, the Advanced Energy Consortium (AEC) at The University of Texas at Austin’s Bureau of Economic Geology (the Bureau) has been pioneering research in subsurface monitoring using its portfolio of distributed autonomous microfabricated sensors for harsh subsurface environments since 2008. A class of these microsensors [System on a Chip (SoC)] can be mixed in cement and permanently placed without compromising the cement column; the sensors would then communicate with each other or a data acquisition (DAQ) master node. Another class of the AEC microsensors can be fully autonomous, with rechargeable micro-batteries capable of exceeding 100°C, flash memory, and, currently, a pressure and temperature sensor. They are designed to circulate in mud, geothermal fluids, U-loops, or pipelines. They can log data into memory and are unobtrusive to operations. Our team has been working on a multi-year DOE-funded project (DE-FE0031856)—supported by $2.95M in federal funding and $0.75M in cost-matching from the AEC—to demonstrate SoC sensor utility for CO2 leakage monitoring in CCS applications. This multi-institutional collaboration developed a novel sensing architecture utilizing radiofrequency (RF) microsensors embedded within the cement sheath. These sensors detect CO2 migration and are interrogated via a Smart Casing Collar (SCC).

58 GEOSCIENCES↗
DOE OSTI2021

An HPC-Based Hydrothermal Finite Element Simulator for Modeling Underground Response to Community-Scale Geothermal Energy Production

Geothermal heat, as renewable energy, shows great advantage with respect to its environmental impact due to its significantly lower CO2 emissions than conventional fossil fuel. Open and closed-loop geothermal heat pumps, which utilize shallow geothermal systems, are an efficient technology for cooling and heating buildings, especially in urban areas. Integrated use of geothermal energy technologies for district heating, cooling, and thermal energy storage can be applied to optimize the subsurface for communities to provide them with multiple sustainable energy and community resilience benefits. The utilization of the subsurface resources may lead to a variation in the underground environment, which might further impact local environmental conditions. However, very few simulators can handle such a highly complex set of coupled computations on a regional or city scale. We have developed high-performance computing (HPC) based hydrothermal finite element (FE) simulator that can simulate the subsurface and its hydrothermal conditions at a scale of tens of km. The HPC simulator enables us to investigate the subsurface thermal and hydrologic response to the built underground environment (such as basements and subways) at the community scale. In this study, a coupled hydrothermal simulator is developed based on the open-source finite element library deal.II. The HPC simulator was validated by comparing the results of a benchmark case study against COMSOL Multiphysics, in which Aquifer Thermal Energy Storage (ATES) is modeled and a process of heat injection into ATES is simulated. The use of an energy pile system at the Treasure Island redevelopment site (San Francisco, CA, USA) was selected as a case study to demonstrate the HPC capability of the developed simulator. The simulator is capable of modeling multiple city-scale geothermal scenarios in a reasonable amount of time.

Sun, Xiang↗