Search NASA⌕ Search

Engineering topics

Fu, Pengcheng

Publications and source records attributed to Fu, Pengcheng.

At least 55 records · Page 3

EGS Collab Experiment 1: Microseismic Monitoring

The U.S. Department of Energy's Enhanced Geothermal System (EGS) Collab project aims to improve our understanding of hydraulic stimulations in crystalline rock for enhanced geothermal energy production through execution of intensely monitored meso-scale experiments. The first experiment is being performed at the 4850 ft level of the Sanford Underground Research Facility (SURF), approximately 1.5 km below the surface at Lead, South Dakota. Here we report on microseismic monitoring of repeated stimulation experiments and subsequent flow tests between two boreholes in the Poorman Formation. Stimulations were performed at several locations in the designated injection borehole at flow rates from 0.1 to 5 L/min over temporal durations from minutes to hours. Microseismic monitoring was performed using a dense 3D sensor array including two cemented hydrophone strings with 12 sensors at 1.75 m spacing accompanied by 18 3-C accelerometers, deployed in 6 monitoring boreholes, completely surrounding the stimulation region. Continuous records were obtained over a two-month period using a novel dual recording system consisting of a conventional 96 channel exploration seismograph and a high-performance 64 channel digitizer sampling sensors at 4 and 100 kHz respectively. Using a standard STA/LTA triggering algorithm, we detected thousands of microseismic events with recorded energy in a frequency range generally above 3 kHz and up to 40 kHz. The locations of these events are consistent with creation of a hydraulic fracture and additional reactivation of pre-existing structures. Using manual pick refinement and double-difference relocation we are able to track the fracture growth to high precision. We estimate the times and locations of the fracture intersecting a monitoring and the production borehole using microseismic events. They are in excellent agreement with independent measurements using distributed temperature sensing, in-situ strain observations and measurements of conductivity changes. This submission includes a microearthquake catalog, raw event files, a subset of the continuous microseismic monitoring data collected during stimulations and flow test activity on 05/22/2018, 05/23/2018, 05/24/2018, 05/25/2018, 06/25/2018, 07/19/2018, 07/20/2018, 12/7/2018, 12/20/2018, and 12/21/2018 (in binary format), and a binary file interpreter to read the continuous microseismic monitoring data. A Stanford Geothermal Workshop paper is also included to describe microseismic monitoring activities at SURF during these periods.

15 GEOTHERMAL ENERGY↗

EGS Collab Experiment 2: Core Logs

Core logs and photos from the EGS Collab project Experiment 2 for the Top Vertical well (TV4100) and the Top Horizontal well (TV 4100) on the 4100 Level of SURF (the Sanford Underground Research Facility). The core logs are stored in a single PDF file with 5-ft run intervals. In the monitoring well IDs, "O" indicates that the well is orthogonal to the anticipated fracture plane, "T" refers to top, and "H" refers to horizontal. A core log CT scan for TV4100 and a layout image of the 4100 wells are included as well. Logs include: experiment number; borehole ID; depth interval; run number; final packed core box number; scribe line (yes/no; red-on-right convention); logging dates; logger initials; as well as sketches of core foliation, folding, and fracturing with additional details and notes on other features of interest. Shift reports include: date, location, personnel, summary of site activity, and field notes.

15 GEOTHERMAL ENERGY↗

Fabric characteristics and processes influencing the liquefaction and re-liquefaction of sand

Although the density of sand has long been considered to be the predominant factor for liquefaction resistance, observations during earthquake sequences and laboratory tests have suggested that other factors may also cause significant anomalies in the relationship between density and liquefaction resistance. This study pinpoints the fabric characteristics and processes influencing the liquefaction resistance of sand. Numerical DEM tests on 2D polydispersed circular granular material are conducted on specimens with various stress-strain and liquefaction histories, showing that changes in micromechanical fabric characteristics of sand can indeed influence its liquefaction resistance, even overpowering the influence of overall density under certain conditions. The liquefaction resistance of initially isotropic sand is strongly correlated with the amount of increase in mean neighboring particle distance needed for it to reach liquefaction (ΔMNPD liq ). Smaller ΔMNPD liq results in weaker liquefaction resistance. Furthermore, fabric anisotropy also significantly affects the liquefaction resistance of sand, due to its influence on dilatancy. Greater fabric anisotropy intensity results in smaller liquefaction resistance. The findings in this study enhances the understanding of liquefaction behavior of sand by establishing a link between the changes in quantifiable fabric states and liquefaction resistance.

58 GEOSCIENCES↗

Recovering the microseismic response from a geomechanical simulation

The assessment of microseismic events is one of a few tools available for characterizing processes that occur within the subsurface. Here, we have developed a method for modeling microseismic activity in the subsurface that allows us to model microseismic events at the reservoir scale. By embedding this method in a fully coupled hydromechanical numerical code, we simulate the development of hydraulic fractures in an unconventional reservoir with an explicitly represented discrete fracture network for a relatively simple synthetic model and a more complicated model based upon field observations. The findings demonstrate that the model can be effectively calibrated against measured microseismic activity, and it can be used to make predictions regarding the timing, location, and amplitude of events in the subsurface.

58 GEOSCIENCES↗

Validation of Numerical Models of the Offshore Wind Turbine from the Alpha Ventus Wind Farm Against Full-Scale Measurements Within OC5 Phase III: Preprint

The main objective of the Offshore Code Comparison Collaboration Continuation, with Correlation (OC5) project, is validation of aero-hydro-servo-elastic simulation tools for offshore wind turbines (OWTs) through comparison of simulated results to the response data of physical systems. Phase III of the OC5 project validates OWT models against the measurements recorded on Senvion 5M wind turbine supported by the OWEC Quattropod from the alpha ventus offshore wind farm. The following operating conditions of the wind turbine were chosen for the validation: (1) Idling below the cut-in wind speed; (2) Rotornacelle assembly (RNA) rotation maneuver below the cut-in wind speed; (3) Power production below and above the rated wind speed; and (4) Shut-down. A number of validation load cases were defined based on these operating conditions. The following measurements were used for validation: (1) Strains and accelerations recorded on the support structure; (2) Pitch, yaw, and azimuth angles, generator speed, and electrical power recorded from the RNA. Strains were not directly available from the majority of the OWT simulation tools. Therefore, strains were calculated based on out-of-plane bending moments, axial forces, and geometrical properties of the structural members. Also, a number of issues arose during the validation: (1) The need for a thorough quality check of sensor measurements; (2) The sensitivity of the turbine loads to the controller and airfoil properties, which were only approximated in the modeling approach; (3) The importance of estimating and applying an appropriate damping value for the structure; and (4) The importance of wind characteristics beyond turbulence on the loads. The simulation results and measurements were compared in terms of time series, discrete Fourier transforms (DFTs), power spectral densities (PSDs), probability density functions (PDFs) of strains and accelerometers. A good match was achieved between the measurements and models set up by OC5 Phase III participants.

17 WIND ENERGY↗

EGS Collab Experiment 1: Circulation Testing

These data and test descriptions comprise a chilled circulation test conducted at the 164' fracture in the EGS Collab Experiment 1 testbed on the 4850 ft level of the Sanford Underground Research Facility. Descriptions of the meta data, design drawings for the flow testing system, and evaluation of the thermistor data are provided here. The test ran from April 2019 through early March of 2020, when testing was concluded at the experiment 1 site. These data are are complementary to the stimulation data provided in another submission which is linked below (i.e. stimulation at the 164' notch). More information about the test itself as well as the rationale and process of data processing is available on the EGS Collab Experiment 1 Long Term Circulation Test wiki page which is also linked below.

15 GEOTHERMAL ENERGY↗

The Necessity for Iteration in the Application of Numerical Simulation to EGS: Examples from the EGS Collab Test Bed 1

The United States Department of Energy, Geothermal Technologies Office (GTO) is funding a collaborative investigation of enhanced geothermal systems (EGS) processes at the meso-scale. This study, referred to as the EGS Collab project, is a unique opportunity for scientists and engineers to investigate the creation of fracture networks and circulation of fluids across those networks under in-situ stress conditions. The EGS Collab project is envisioned to comprise three experiments and the site for the first experiment is on the 4850 Level in phyllite of the Precambrian Poorman formation, at the Sanford Underground Research Facility, located at the former Homestake Gold Mine, in Lead, South Dakota. Principal objectives of the project are to develop a number of intermediate-scale field sites and to conduct well-controlled in situ experiments focused on rock fracture behavior and permeability enhancement. Data generated during these experiments will be compared against predictions of a suite of computer codes specifically designed to solve problems involving coupled thermal, hydrological, geomechanical, and geochemical processes. Comparisons between experimental and numerical simulation results will provide code developers with direction for improvements and verification of process models, build confidence in the suite of available numerical tools, and ultimately identify critical future development needs for the geothermal modeling community. Moreover, conducting thorough comparisons of models, modelling approaches, measurement approaches and measured data, via the EGS Collab project, will serve to identify techniques that are most likely to succeed at the Frontier Observatory for Research in Geothermal Energy (FORGE), the GTO’s flagship EGS research effort. As noted, outcomes from the EGS Collab project experiments will serve as benchmarks for computer code verification, but numerical simulation additionally plays an essential role in designing these meso-scale experiments. This paper reviews specific numerical simulations supporting the design of experiments within Test Bed 1, a volume of phyllite rock under in-situ stress conditions off the western side of the West Access Drift on the 4850 Level, near Governor’s Corner. Numerical simulations were executed prior to the start of hydraulic stimulation activities within Test Bed 1 following standard practices, using best estimates of principal stress conditions, thermal conditions, and the rock petrophysical properties, including geomechanical properties. These simulations indicated notching of the borehole would promote the initiation of transverse fractures, seismic magnitudes during the hydraulic stimulation would be below 0.1 magnitude on the Richter scale, a rock thermal conductivity of 5.0 W/m K yielded agreement with kISMET borehole temperature logs, back pressure on the production borehole would increase circulation across the test bed, and the production borehole would serve to halt propagation of the hydraulic fracture to the drift. Pre-existing natural fractures, heterogeneities in the rock properties, monitoring boreholes, and overlooked mine elements have prompted a second look at numerically modeling stimulation, fluid circulation, tracer migration, and thermal breakthrough. Numerical simulation is an invaluable tool for providing insight and understanding to complex physical processes. The success of simulations, however, often depends on including all of the salient features of the system in the founding conceptual model. This paper takes a retrospective look at examples where the conceptual model and simulation was sufficient to provide accurate forecasts and those where elements were missing, necessitating rethinking of the simulation.

numerical simulation, EGS Collab, enhanced geother↗

The Necessity for Iteration in the Application of Numerical Simulation to EGS: Examples from the EGS Collab Test Bed 1

The United States Department of Energy, Geothermal Technologies Office (GTO) is funding a collaborative investigation of enhanced geothermal systems (EGS) processes at the meso-scale. This study, referred to as the EGS Collab project, is a unique opportunity for scientists and engineers to investigate the creation of fracture networks and circulation of fluids across those networks under in-situ stress conditions. The EGS Collab project is envisioned to comprise three experiments and the site for the first experiment is on the 4850 Level in phyllite of the Precambrian Poorman formation, at the Sanford Underground Research Facility, located at the former Homestake Gold Mine, in Lead, South Dakota. Principal objectives of the project are to develop a number of intermediate-scale field sites and to conduct well-controlled in situ experiments focused on rock fracture behavior and permeability enhancement. Data generated during these experiments will be compared against predictions of a suite of computer codes specifically designed to solve problems involving coupled thermal, hydrological, geomechanical, and geochemical processes. Comparisons between experimental and numerical simulation results will provide code developers with direction for improvements and verification of process models, build confidence in the suite of available numerical tools, and ultimately identify critical future development needs for the geothermal modeling community. Moreover, conducting thorough comparisons of models, modelling approaches, measurement approaches and measured data, via the EGS Collab project, will serve to identify techniques that are most likely to succeed at the Frontier Observatory for Research in Geothermal Energy (FORGE), the GTO's flagship EGS research effort. As noted, outcomes from the EGS Collab project experiments will serve as benchmarks for computer code verification, but numerical simulation additionally plays an essential role in designing these meso-scale experiments. This paper reviews specific numerical simulations supporting the design of experiments within Test Bed 1, a volume of phyllite rock under in-situ stress conditions off the western side of the West Access Drift on the 4850 Level, near Governor's Corner. Numerical simulations were executed prior to the start of hydraulic stimulation activities within Test Bed 1 following standard practices, using best estimates of principal stress conditions, thermal conditions, and the rock petrophysical properties, including geomechanical properties. These simulations indicated notching of the borehole would promote the initiation of transverse fractures, seismic magnitudes during the hydraulic stimulation would be below 0.1 magnitude on the Richter scale, a rock thermal conductivity of 5.0 W/m K yielded agreement with kISMET borehole temperature logs, back pressure on the production borehole would increase circulation across the test bed, and the production borehole would serve to halt propagation of the hydraulic fracture to the drift. Pre-existing natural fractures, heterogeneities in the rock properties, monitoring boreholes, and overlooked mine elements have prompted a second look at numerically modeling stimulation, fluid circulation, tracer migration, and thermal breakthrough. Numerical simulation is an invaluable tool for providing insight and understanding to complex physical processes. The success of simulations, however, often depends on including all of the salient features of the system in the founding conceptual model. This paper takes a retrospective look at examples where the conceptual model and simulation was sufficient to provide accurate forecasts and those where elements were missing, necessitating rethinking of the simulation.

15 GEOTHERMAL ENERGY↗

Deformation of Granular Material under Continuous Rotation of Stress Principal Axes

The mechanical response of granular material to the rotation of stress principal axes is an issue of both practical and theoretical importance in soil mechanics. This study used a two-dimensional (2D) discrete-element method (DEM) simulation to investigate the deformation of granular material under the continuous rotation of stress principal axes while maintaining fixed principal stress values. The results show that under such rotation, the deformation can exceed that caused by fixed principal stress axes cyclic biaxial compression with a maximum deviatoric stress ratio that equals the fixed stress ratio during rotation. The volumetric strain was found to be contractive overall while oscillating within each load cycle for specimens with significant fabric anisotropy. Here, the initial fabric-anisotropy orientation was found to have little influence on the overall evolution of volumetric strain, although it was found to affect the development of shear strain. During the cycles of rotation of initial stress principal axes, the dilatancy of granular material, in the contractive direction, was found to be dominated by the evolving orientation of the major principal stress axis in reference to the initial major principal stress axis during anisotropic consolidation, and the influence of the initial fabric-anisotropy orientation was found to be trivial. As stress rotation continued, the influence of the relative orientation between stress and fabric became prominent, causing the material to dilate when the major principal stress axis rotated from the normal of the bedding plane to being perpendicular to it and causing contraction during the other half of the cycle.

58 GEOSCIENCES↗

Hydraulic-Fracture-Height Growth Under the Combined Influence of Stress Barriers and Natural Fractures

A fully coupled finite-element/finite-volume code is used here to model 3D hydraulically driven fractures under the influence of strong vertical variations in closure stress interacting with natural fractures. Previously unknown 3D interaction mechanisms on fracture-height growth are revealed. Slipping of a natural fracture, triggered by elevated fluid pressure from an intersecting hydraulic fracture, can induce both increases and decreases of normal stress in the minimum-horizontal-stress direction, toward the center and tip of the natural fracture, respectively. Consequently, natural fractures are expected to be able to both encourage and inhibit the progress of hydraulic fractures propagating through stress barriers, depending on the relative locations between the intersecting fractures. Once the hydraulic fracture propagates above the stress barrier through the weakened segment near a favorably located natural fracture, a configuration consisting of two opposing fractures cuts the stress barrier from above and below. The fluid pressure required to break the stress barrier under such opposing-fracture configurations is substantially lower than that required by a fracture penetrating the same barrier from one side. Sensitivity studies of geologic conditions and operational parameters have also been performed to explore the feasibility of controlled fracture height. The interactions between hydraulic fractures, natural fractures, and geologic factors such as stress barriers in three dimensions are shown to be much more complex than in two dimensions. Although it is impossible to exhaust all the possible configurations, the ability of a 3D, fully coupled numerical model to naturally capture these processes is well-demonstrated.

42 ENGINEERING↗

Evaluating a Simple Fracturing Criterion for a Hydraulic Fracture Crossing Stress and Stiffness Contrasts

Hydraulic fracture height containment is a critical issue in the development of unconventional reservoirs. The extent of fracture height growth depends on a variety of factors, particularly stress and stiffness contrasts between adjacent layers. Accurate simulation of fracture growth and containment requires a reliable fracturing criterion. The virtual crack closure technique (VCCT) is a widely used method for computing energy release rate. However, it is based on the assumption that a small crack extension does not significantly alter the state of the crack tip, which is generally not the case when a fracture crosses strong stress and/or stiffness contrasts. In this work, we assess the applicability and accuracy of a modified virtual crack closure technique (MVCCT) for a fluid-driven fracture in breaking through interfaces with significant stress and/or stiffness contrasts, through comparisons with analytical and reference numerical solutions. Here, the results show that significant error could occur when the fracture tip is very near or at stress/stiffness interfaces. However, this error is localized to the interface and proves to be inconsequential to the predicted penetration depth into the rock layer beyond the interface. Finally, this study validates the applicability of MVCCT in 3D hydraulic fracturing simulation in strongly heterogeneous rock formations.

36 MATERIALS SCIENCE↗

Toward better hydraulic fracturing fluids and their application in energy production: A review of sustainable technologies and reduction of potential environmental impacts

Recent advances in hydraulic fracturing, in conjunction with horizontal drilling, have enabled large-scale extraction of natural gas and oil from shale formations. Despite its advances and enormous economic benefits, opportunities remain to increase hydraulic fracturing efficiency and minimize potential environmental impacts. This review specifically examines three key themes associated with development and utilization of hydraulic fracturing fluids: 1) characteristics and behavior of fracturing fluids, 2) understanding and predicting migration and fate of fracturing fluids, 3) technologies to reduce environmental impact of fracturing fluids. Furthermore, the paper discusses key and new techniques and findings on rheology of hydrogel-based fluids, high fidelity simulation of propagation transport, potential environmental impacts, geosynthetics in mitigating contamination, and greener fracturing fluids. It is indicated that future development relies on advances in understanding of physical processes, modeling capabilities, and monitoring techniques.

02 PETROLEUM↗

Verification of a Numerical Model of the Offshore Wind Turbine from the Alpha Ventus Wind Farm Within OC5 Phase III: Preprint

The main objective of the Offshore Code Comparison Collaboration Continuation, with Correlation (OC5) project, is validation of aero-hydro-servo-elastic simulation tools for offshore wind turbines (OWTs) through comparison of simulated results to the response data of physical systems. Phase III of the OC5 project analyzes the Senvion 5M wind turbine supported by the OWEC Quattropod from the alpha ventus offshore wind farm. This paper shows results of the verification of the OWT models (code-to-code comparison). A subsequent publication will focus on their validation (comparison of simulated results to measured physical system response data). Based on the available data, the participants of Phase III set up numerical models of the OWT in their simulation tools. It was necessary to verify and to tune these models. The verification and tuning were performed against an OWT model available at the University of Stuttgart - Stuttgart Wind Energy (SWE) and documentation provided by Senvion and OWEC Tower. A very good match was achieved between the results from the reference SWE model and models set up by OC5 Phase III participants.

17 WIND ENERGY↗

Verification of a Numerical Model of the Offshore Wind Turbine From the Alpha Ventus Wind Farm Within OC5 Phase III

The main objective of the Offshore Code Comparison Collaboration Continuation, with Correlation (OC5) project, is validation of aero-hydro-servo-elastic simulation tools for offshore wind turbines (OWTs) through comparison of simulated results to the response data of physical systems. Phase III of the OC5 project analyzes the Senvion 5M wind turbine supported by the OWEC Quattropod from the alpha ventus offshore wind farm. This paper shows results of the verification of the OWT models (code-to-code comparison). A subsequent publication will focus on their validation (comparison of simulated results to measured physical system response data). Based on the available data, the participants of Phase III set up numerical models of the OWT in their simulation tools. It was necessary to verify and to tune these models. The verification and tuning were performed against an OWT model available at the University of Stuttgart - Stuttgart Wind Energy (SWE) and documentation provided by Senvion and OWEC Tower. A very good match was achieved between the results from the reference SWE model and models set up by OC5 Phase III participants.

17 WIND ENERGY↗

Inferring Geothermal Reservoir Processes at the Raft River Geothermal Field, Idaho, USA, Through Modeling InSAR-Measured Surface Deformation

Ground surface deformations detected with Interferometric Synthetic Aperture Radar (InSAR) provide valuable information for inferring subsurface reservoir processes that are difficult to observe directly. This study aims at building a reservoir model that honors the available geological, hydrological and geo-mechanical data and also produces ground surface deformation consistent with InSAR measurements at the Raft River Geothermal Field. In the thermo-hydro-mechanical (THM) coupled model developed, the reservoir deforms as a result of the rock’s poroelastic response to changes in hydrologic pressure and thermal expansion/contraction. The results indicate that the observed deformation is the result of pressure decrease in the deep production reservoir and pressure increase in a shallower reservoir that accommodates the injected fluid (likely in the Salt Lake Formation). The combination of the uplift around injection wells with subsidence around the production wells, affected by the Bridge Fault as a flow barrier, creates a complex pattern of surface deformation in which the center of subtle subsidence significantly deviates from the location of the production wells. A parametric study suggests that (1) the Bridge Fault Zone is likely a flow barrier, (2) the surface deformation appears to be insensitive to the presence of the Narrows Structure, and (3) additional flow barriers likely exist to shape the flow system. This case study demonstrates the utility of a high-fidelity forward model that honors available known information and THM coupled processes in understanding geothermal reservoir characteristics.

58 GEOSCIENCES↗

Numerical Simulation Applications in the Design of EGS Collab Experiment 1

The United States Department of Energy, Geothermal Technologies Office (GTO) is funding a collaborative investigation of enhanced geothermal systems (EGS) processes at the meso-scale. This study, referred to as the EGS Collab project, is a unique opportunity for scientists and engineers to investigate the creation of fracture networks and circulation of fluids across those networks under in-situ stress conditions. The EGS Collab project is envisioned to comprise three experiments and the site for the first experiment is on the 4850 Level (4,850 feet below ground surface) in phyllite of the Precambrian Poorman formation, at the Sanford Underground Research Facility, located at the former Homestake Gold Mine, in Lead, South Dakota. Principal objectives of the project are to develop a number of intermediate-scale field sites and to conduct well-controlled in situ experiments focused on rock fracture behavior and permeability enhancement. Data generated during these experiments will be compared against predictions of a suite of computer codes specifically designed to solve problems involving coupled thermal, hydrological, geomechanical, and geochemical processes. Comparisons between experimental and numerical simulation results will provide code developers with direction for improvements and verification of process models, build confidence in the suite of available numerical tools, and ultimately identify critical future development needs for the geothermal modeling community. Moreover, conducting thorough comparisons of models, modelling approaches, measurement approaches and measured data, via the EGS Collab project, will serve to identify techniques that are most likely to succeed at the Frontier Observatory for Research in Geothermal Energy (FORGE), the GTO's flagship EGS research effort. As noted, outcomes from the EGS Collab project experiments will serve as benchmarks for computer code verification, but numerical simulation additionally plays an essential role in designing these meso-scale experiments. This paper describes specific numerical simulations supporting the design of Experiment 1, a field test involving hydraulic stimulation of two fractures from notched sections of the injection borehole and fluid circulation between sub-horizontal injection and production boreholes in each fracture individually and collectively, including the circulation of chilled water. Whereas the mine drift allows for accurate and close placement of monitoring instrumentation to the developed fractures, active ventilation in the drift cooled the rock mass within the experimental volume. Numerical simulations were executed to predict seismic events and magnitudes during stimulation, initial fracture orientations for smooth horizontal wellbores, pressure requirements for fracture initiation from notched wellbores, fracture propagation during stimulation between the injection and production boreholes, tracer travel times between the injection and production boreholes, produced fluid temperatures with chilled water injections, pressure limits on fluid circulation to avoid fracture growth, temperature environment surrounding the 4850 Level drift, and fracture propagation within a stress field altered by drift excavation, ventilation cooling, and dewatering.

15 GEOTHERMAL ENERGY↗

Outcomes from a Collaborative Approach to a Code Comparison Study for Enhanced Geothermal Systems

A code comparison study, supported by the U.S. Department of Energy, Geothermal Technologies Office, has recently been completed whose stated purpose was to test, diagnose differences, and demonstrate the state of modeling capabilities developed at national laboratories and academic institutes within the United States for modeling enhanced geothermal systems. This three-year study comprised two sets of problems; those classified as benchmark problems and those classified as challenge problems. The benchmark problems were structured to test the ability of the collection of numerical simulators to solve various combinations of coupled thermal, hydrologic, geomechanical, and geochemical processes. This class of problems was strictly defined in terms of properties, driving forces, initial conditions, and boundary conditions. Study participants submitted solutions to problems for which their simulation tools were deemed capable or nearly capable. Some participating codes were originally developed for EGS applications whereas some others were designed for different applications but can simulate processes similar to those in EGS. Solution submissions from both were encouraged. In some cases participants made small incremental changes to their numerical simulation codes to address specific elements of the problem, and in other cases participants submitted solutions with existing simulation tools, acknowledging the limitations of the code. The challenge problems were based on the enhanced geothermal systems research conducted in hot dry rock at Fenton Hill, near Los Alamos, New Mexico, between 1974 and 1995, covering two research stimulation, development and circulation phases in two separate reservoirs. Both challenge problems have specific questions to be answered via numerical simulation in three topical areas: 1) reservoir creation/stimulation, 2) reactive and passive transport, and 3) thermal recovery. Whereas the benchmark class of problems were designed to test capabilities for modeling coupled processes under strictly specified conditions, the stated objective for the challenge class of problems was to demonstrate what new understanding of the Fenton Hill experiments could be realized via the application of modern numerical simulation tools by recognized expert practitioners. Critical observations and data from the Fenton Hill experiments were varied and scattered among a number of individual experiments, and numerical simulation solutions were sought that satisfied as many observations at possible. Achieving agreement in multiple experimental observations from Fenton Hill often required participants to critically think about stimulation mechanisms involving natural and hydraulic fractures and re-evaluate their conceptual models and numerical solution approaches. This process has yielded new insights to Fenton Hill reservoirs and direction for future EGS research. The keys to the success of this study were the high quality numerical simulation tools available, the expertise of the developers and practitioners, and the collaborative approach to reporting preliminary results among the study participants. This paper describes results submitted to the challenge problems and the outcomes of a collaborative approach to conducting a code comparison study.

Enhanced Geothermal Systems, numerical simulation,↗

Numerical Simulation Applications in the Design of EGS Collab Experiment 1

The United States Department of Energy, Geothermal Technologies Office (GTO) is funding a collaborative investigation of enhanced geothermal systems (EGS) processes at the meso-scale. This study, referred to as the EGS Collab project, is a unique opportunity for scientists and engineers to investigate the creation of fracture networks and circulation of fluids across those networks under in-situ stress conditions. The EGS Collab project is envisioned to comprise three experiments and the site for the first experiment is on the 4850 Level in phyllite of the Precambrian Poorman formation, at the Sanford Underground Research Facility, located at the former Homestake Gold Mine, in Lead, South Dakota. Principal objectives of the project are to develop a number of intermediate-scale field sites and to conduct well-controlled in situ experiments focused on rock fracture behavior and permeability enhancement. Data generated during these experiments will be compared against predictions of a suite of computer codes specifically designed to solve problems involving coupled thermal, hydrological, geomechanical, and geochemical processes. Comparisons between experimental and numerical simulation results will provide code developers with direction for improvements and verification of process models, build confidence in the suite of available numerical tools, and ultimately identify critical future development needs for the geothermal modeling community. Moreover, conducting thorough comparisons of models, modelling approaches, measurement approaches and measured data, via the EGS Collab project, will serve to identify techniques that are most likely to succeed at the Frontier Observatory for Research in Geothermal Energy (FORGE), the GTO’s flagship EGS research effort. As noted, outcomes from the EGS Collab project experiments will serve as benchmarks for computer code verification, but numerical simulation additionally plays an essential role in designing these meso-scale experiments. This paper describes specific numerical simulations supporting the design of Experiment 1, a field test involving hydraulic stimulation of two fractures from notched sections of the injection borehole and fluid circulation between sub-horizontal injection and production boreholes in each fracture individually and collectively, including the circulation of chilled water. Whereas the mine drift allows for accurate and close placement of monitoring instrumentation to the developed fractures, active ventilation in the drift cooled the rock mass within the experimental volume. Numerical simulations were executed to predict seismic events and magnitudes during stimulation, initial fracture orientations for smooth horizontal wellbores, pressure requirements for fracture initiation from notched wellbores, fracture propagation during stimulation between the injection and production boreholes, tracer travel times between the injection and production boreholes, and produced fluid temperatures with chilled water injections, and pressure limits on fluid circulation to avoid fracture growth.

Numerical Simulation, EGS Collab, EGS, Sanford Und↗