Engineering Papers⌕ Search

Engineering topics

Fu, Pengcheng

Publications and source records attributed to Fu, Pengcheng.

At least 37 records · Page 2
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Developing Upscaling Approach for Swarming Hydraulic Fractures Observed at Hydraulic Fracturing Test Site through Multiscale Simulations

We show that this work aims to address a challenge posed by recent observations of tightly spaced hydraulic fractures in core samples from the hydraulic fracturing test site (HFTS) in the Middle Wolfcamp Formation. Many fractures in retrieved cores have subfoot spacing, which is at odds with conventional models in which usually one hydraulic fracture is initiated per cluster. Models assuming a single fracture at each cluster, although a common practice, often predict excessive fracture propagation that is inconsistent with microseismic observation. Here, we aim to develop a numerical approach to effectively account for densely spaced hydraulic fractures in field-scale simulations. Because it is impractical to explicitly model all aforementioned fractures, we develop a new upscaling law that enables existing simulation tools to predict reservoir response to fracture swarms. The upscaling law is derived based on an energy equivalence argument and validated through multiscale simulations using a high-fidelity code, GEOS. The swarming fractures are first modeled with a spacing that is much smaller than the cluster spacing; these fractures are then approximated by an upscaled, single fracture based on the proposed upscaling law. The upscaled fracture is shown to successfully match the energy input rate and produce the total fracture aperture and average propagation length of the explicitly simulated swarm. Afterward, the upscaling approach is further implemented in 3D field-scale simulations and validated against the HFTS microseismic data of a horizontal well. Our results show that hydraulic fracture swarming can significantly affect fracture propagation behaviors compared with the propagation of single fractures as assumed by conventional modeling approaches. Under the considered situations, the conventional treatment yields fast propagation speed that far exceeds that indicated by the microseismic data. We also illustrate that this discrepancy can be reduced readily through the implementation of the upscaling law. Our results demonstrate the importance of accounting for the fracture swarming effect in field-scale simulations and the efficacy of this approach to enable realistic predictions of reservoir responses to fracture swarms, without the need to model tightly spaced fractures individually.

02 PETROLEUM↗
DOE OSTI2020

Characterization of flow and transport in a fracture network at the EGS Collab field experiment through stochastic modeling of tracer recovery

Energy extraction from subsurface reservoirs is important for addressing the increasing energy demand and environmental concerns such as global warming. However, the characterization of subsurface reservoirs, particularly reservoirs dominated by fracture networks remains a challenge due to the lack of means to directly observe subsurface processes. This study explores the feasibility and efficacy of characterizing fracture flow and transport processes in an enhanced geothermal system (EGS) testbed through stochastic tracer modeling. There are two enabling factors that allow application of stochastic modeling to characterize a subsurface reservoir. First, an abundance of geological and geophysical measurements enables the development of a high-fidelity and well-constrained fracture network model. Second, high-performance computing (HPC) allows running massive realizations efficiently. Six conservative tracer tests were stochastically modeled and produced satisfactory realizations that successfully reproduce field tracer recovery data from each tracer test. The evolution of flow and transport processes in the fracture network was then analyzed from these satisfactory realizations. This work demonstrates that stochastic tracer modeling on a high-fidelity fracture network model is feasible and can provide important insights regarding flow and transport characteristics in subsurface fractured reservoirs.

58 GEOSCIENCES↗
DOE OSTI2020

Fracture Caging to Limit Induced Seismicity

Geothermal resources offer a stable low-carbon energy source. However, geothermal sites can collocate with the hypocenters of large-magnitude seismic events. Large seismic events pose a risk to public safety and are therefore a liability for efforts to develop geothermal resources. In this work, we propose “fracture caging” to limit induced seismic event magnitudes and present evidence from numerical model predictions, laboratory experiments, and field observations. Fracture caging involves drilling tactical production wells around a geothermal injection zone to contain fluids in fracture-dominated flow systems. Prior to our work, the effect of small wells on the growth of large fractures and on flow through fractures was subject to debate. Our work shows that production wells can impede fracture growth and contain high-pressure fluids in fracture-dominated rocks. This containment offers a mechanism to limit induced seismicity.

58 GEOSCIENCES↗
DOE OSTI2020

CO2 injection with fracturing in geomechanically protected caprock: Task 6 of LLNL's Research Activities to Support DOE's Carbon Storage Program (FWP-FEW0191) (Final Technical Report)

Conventional principles of the design and operation of geologic carbon storage (GCS) require injecting CO 2 below the caprock fracturing pressure to ensure the integrity of the storage complex. In non-ideal storage reservoirs with relatively low permeability, modest injection rates can lead to pressure buildup and hydraulic fracturing of the reservoir and caprock. While the GCS community has generally viewed hydraulic fractures as a key risk to storage integrity, the actual behavior of a caprock hydraulic fracture, particularly from a geomechanical perspective, has not been thoroughly studied. The prevailing method of treating a hydraulic fracture as a high-permeability wing of the storage reservoir might have resulted in erroneous understanding of caprock hydraulic fracture behavior. Comprehensive analyses of monitoring data from the In Salah project suggested that one or more hydraulic fractures may have been created in the reservoir and lower caprock system during injection operations, but did not cause detectable leakage of CO 2 out of the storage complex. This observation coincides with a well-known phenomenon in the oil and gas industry: sedimentary rock formations have many inherent features that naturally protect the formations from unbounded vertical growth of hydraulic fractures.

58 GEOSCIENCES↗
DOE OSTI2020

EGS Collab Experiment 1 Stimulation Data

Stimulation data from Experiment 1 of EGS Collab, which occurred on the 4850 ft level of the Sanford Underground Research Facility (SURF). A detailed description of the stimulation data is provided in the StimulationDataNotes.docx and is also available on the EGS Collab Wiki. A Meta Data Cheat Sheet, which describes all of the channels in the Raw CSV files, is available as well. Note that this cheat sheet is a comprehensive meta data descriptor and channels were added as the experiment evolved. This means that some columns may not be populated in early data. Additionally, we have included the chat logs from these experiments. The experiments were broadcast over teleconferencing software and real-time data displays were available to remote observers. The logs contain important observations from those personnel performing the experiment and the remote contributors. Finally, we have included summary and individual plots of all of the data for the user to compare to.

15 GEOTHERMAL ENERGY↗
DOE OSTI2020

Validation of Numerical Models of the Offshore Wind Turbine From the Alpha Ventus Wind Farm Against Full-Scale Measurements 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 validates OWT models against the measurements recorded on a 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) rotor-nacelle assembly (RNA) rotation maneuver below the cut-in wind speed, (3) power production below and above the rated wind speed, and (4) shutdown. 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 and (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 cross-sectional properties of the structural members. The simulation results and measurements were compared in terms of time series, discrete Fourier transforms, power spectral densities, and probability density functions of strains and accelerometers. A good match was achieved between the measurements and models setup by OC5 Phase III participants.

17 WIND ENERGY↗
DOE OSTI2020

Portable interactive visualization of large-scale simulations in geotechnical engineering using Unity3D

Development in large-scale geotechnical engineering simulation places tremendous demand for efficient visualization of such simulation data. This study presents a lightweight software tool, i.e. Geotechnical Interactive Visualization (GIV), as a solution to this challenge, which achieves efficient interactive visualization of large-scale simulations data in geotechnical engineering. Visualization data flow and algorithms specifically optimized for common geotechnical engineering applications are implemented in GIV. GIV can visualize geotechnical structure models with time-varying attributes attached to mesh with fixed topology, and also models with time-varying mesh topologies but no attributes attached, the two most common visualization tasks in geotechnical engineering. Furthermore, challenges for large-scale simulation data visualization, including parallel simulation data redundancy, massive data size, dynamic user interaction, and portability are overcome via specifically designed algorithms for simulation data preprocessing and optimized visualization modules using the powerful 3D rendering and interactive game engine Unity3D. Comparison of GIV with several widely used visualization tools for the visualization of large-scale idealized datasets and realistic geotechnical simulations highlights the visualization efficiency, smooth interactivity, and lightweight features of GIV.

42 ENGINEERING↗
DOE OSTI2020

Gas Production from Hot Water Circulation through Hydraulic Fractures in Methane Hydrate-Bearing Sediments: THC-Coupled Simulation of Production Mechanisms

Methane hydrates, widely found in permafrost and deep marine sediments, have great potential as a future energy 6 source. Conventional production schemes perform poorly for challenging hydrate reservoirs with low permeability. We propose an 7 efficient production scheme by combining hydraulic fracturing from horizontal wells and hot water circulation through fractures. A 8 fully coupled thermo-hydro-chemical (THC) model is developed to simulate the key physical processes during gas production from 9 a hydrate reservoir representative of typical geological settings in Shenhu, South China Sea. We found that the gas production 10 process has two distinct stages divided by thermal breakthrough: a relatively short prebreakthrough stage and a postbreakthrough 11 stage yielding stable gas production. Heat advection along and near the hydraulic fracture dominates the prebreakthrough stage, 12 whereas conduction-driven thermal recovery in the volume around fractures dominates the postbreakthrough stage. We identified 13 that the steady-state injection temperature has a strong effect on the performance of the proposed scheme while the fluid mass 14 circulation rate has a moderate impact beyond a threshold. The proposed scheme proves to be efficient and robust over a range of 15 reservoir conditions with respect to initial hydrate saturation and intrinsic permeability, including their spatial heterogeneities, 16 thereby offering a promising solution for challenging reservoir conditions.

gas hydrate, fracture flow, gas production, numeri↗
DOE OSTI2020

The EGS Collab Project: Learnings from Experiment 1

The primary objective of the EGS Collab Project sponsored by DOE is to increase the understanding needed to efficiently implement enhanced geothermal systems (EGS). One goal of the EGS Collab project is to create a collaborative research environment in which to study stimulation of crystalline rock at the 10 meter scale. Key to this effort is the collection of high quality data to allow comparison to numerical coupled process models in an effort to build confidence in the codes and modeling techniques used. In response to this, the EGS Collab team has created an underground test bed at the Sanford Underground Research Facility (SURF) in Lead SD at a depth of approximately 1.5 km to examine hydraulic fracturing (Experiment 1). We are currently designing a second test bed aimed at investigating shear stimulation (Experiment 2). At the Experiment 1 location, we have characterized our host rock using laboratory testing and numerous field-based geophysical and geological techniques, and created a well-instrumented test bed to allow us to carefully monitor stimulation events and flow tests. In addition to the installed geophysical sensors, we have used tracer tests, differences in the ambient microbial communities at flow collection locations, and cold water injection to inform us about dynamic flow pathways. In Experiment 1, we have hydraulically stimulated the host rock in a number times at several locations in one well, creating new fractures that connect to existing fractures between the injection and production boreholes. We have performed long-term ambient and chilled water injection tests as an analog to EGS, and have monitored system changes resulting from these water injections through geophysical monitoring, flow and pressure measurements, tracer tests, and microbiology. Here, we summarize the tests performed, issues identified including poroelastic and thermoelastic effects, Joule-Thomson effects, restarting effects, indications of flow channeling, and the primary learnings from Experiment 1.

Enhanced Geothermal Systems, EGS Collab, stimulati↗
DOE OSTI2020

Application of an Embedded Fracture and Borehole Modeling Approach to the Understanding of EGS Collab Experiment 1

EGS Collab is a series of meso-scale experiments and associated numerical simulation activities being funded by the United States Department of Energy, Geothermal Technologies Office (GTO) to investigate enhanced geothermal system processes under in-situ stress and slightly elevated temperature conditions. This project is designed to provide scientists and engineers with immediate access to impermeable rock at scales larger than possible in the laboratory, but generally smaller than those for commercial production. Immediate access to rock is provided via the existing drifts of the former Homestake Gold Mine, now operated as the Sanford Underground Research Facility in Lead, South Dakota. The objectives of the EGS Collab project are to develop well controlled fracture networks between injector and producer boreholes using normal and shear stimulation for permeability enhancement. The first experimental site was located off the West Access drift on the 4850 Level (4850 feet below ground surface) in phyllite of the Precambrian Poorman formation, and involved the creation of a fracture network comprising a combination of hydraulic and natural fractures. A second experimental site is now being considered near the battery alcove on the 4100 Level in amphibolite and rhyolite of the Yates Unit. 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. Experiment 1 has comprised a series of successful tests, including a long-term chilled-water circulation test, but the testbed has two atypical EGS elements. Active ventilation in the adjacent drift over a 50-year period cooled the testbed rock mass and hydraulic fracturing intersected monitoring boreholes making them conduits for fluid flow. Numerical simulations executed in support of the design of the EGS Collab Experiment 1, computed a radial temperature distribution and stress gradient orthogonal to the drift axis, resulting in the forecast of an oblong hydraulic fracture geometry, extended in the direction of the drift from the stimulation borehole. This paper describes the numerical simulation of long-term chilled-water circulation test with an embedded fracture and borehole modeling approach. The principal objective of the simulation work is show agreement between experimental observations in terms of production fluid temperatures, monitoring borehole temperatures, injection and production pressures, and tracer recoveries using a single conceptual model for the fracture network, using characterization and monitoring data generally available for EGS. A secondary objective is to improve the numerical simulation result comparisons with additional data available from the broader monitoring equipment within the testbed versus deeper, hotter, and more remote EGS. One complicating factor for comparisons between the numerical simulation results and experimental observations, is the Joule-Thomson heating associated with large pressure drops across the fracture network.

numerical simulation, embedded fractures and boreh↗
DOE OSTI2020

A spatial-statistical investigation of surface expressions associated with cyclic steaming in the Midway-Sunset Oil Field, California

In the Midway Sunset Oil Field in Central California, operators inject steam into the shallow diatomite formation to enhance heavy oil recovery through imbibition, wettability alteration, and viscosity reduction, among other mechanisms. The injected steam, however, does not always remain in the reservoir or return through the wells. In two zones in the study area, the steam comes out at the surface, creating sinkholes, seeps, and steam outlets. These phenomena, called “surface expressions,” pose safety and environmental hazards. Even though these surface expressions are a widespread problem in Central California, they are not well documented and understood. Possible causes of the surface expressions include: high injection pressure, structurally controlled flow patterns, leakage of steam through old improperly abandoned wells, high injection volumes, or flow along naturally occurring faults, among other possible factors. This work examines attributes of the zones with surface expressions in order to determine factors that may contribute to their occurrence. Spatial statistical analysis using logistic regression, random forests, and classification trees is used to explore the relationship between the surface expressions and geological and production-related attributes. The results point to a significant spatial correlation between the surface expressions and two predictors: concentration of plugged wells and geologic seal thickness. The results guide follow-up studies to further investigate the role of well abandonment and seal thickness in the occurrence of surface expressions.

02 PETROLEUM↗
DOE OSTI2019

Validation of Numerical Models of the Offshore Wind Turbine From the Alpha Ventus Wind Farm Against Full-Scale Measurements 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 validates OWT models against the measurements recorded on a 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) Rotor-nacelle assembly (RNA) rotation maneuver below the cut-in wind speed; (3) Power production below and above the rated wind speed; and (4) Shutdown. 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 cross-sectional 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, power spectral densities, probability density functions of strains and accelerometers. A good match was achieved between the measurements and models set up by OC5 Phase III participants.

17 WIND ENERGY↗
DOE OSTI2019

EGS Collab Project: Accomplishments and Plan

To better understand stimulation and inter-well flow in crystalline rock with the goal of developing and improving subsurface heat exchangers for enhanced geothermal energy recovery, the EGS (Enhanced Geothermal Systems) Collab project is developing highly monitored, well- characterized intermediate scale (~10-20 m) field test beds, stimulating, and performing flow tests. These tests provide the ability to evaluate concepts, data to test and validate models, and tools to be used for EGS development, which may be evaluated at DOE’s Frontier Observatory for Research in Geothermal Energy (FORGE) project. Numerical simulation is key to the project, for experiment design and interpretation of results. Testbed characterization included optical and acoustic televiewers, full waveform seismic, electromagnetic, gamma, and temperature. The test block was further characterized using seismic tomography (p- and s-), electrical resistance tomography (ERT), and hydrologic characterization including tracer tests. Coupled process models are constrained by the detailed site characterization data. Installed monitoring systems include: 1) passive seismic/acoustic emissions; 2) continuous active source seismic; 3) ERT, with dynamic imaging using high contrast fluids; and 4) distributed acoustic, temperature, and strain fiber optic sensors. During flow and stimulation tests, fracture aperture strain monitoring is performed using the Step-rate Injection Method for Fracture In-situ Properties (SIMFIP) tool. Pressure and flow in the injection and production boreholes are monitored, and tracer tests and wavefield imaging and inversion are performed to provide additional constraints. Stimulations have been performed at several locations, connecting our injection and production wells. An overview of the EGS Collab project is presented to provide a framework for other EGS Collab presentations, along with results from numerical simulations and experiments performed at the Sanford Underground Research Facility in South Dakota.

Enhanced Geothermal Systems, EGS Collab, stimulati↗
DOE OSTI2019

Fabric evolution and dilatancy within anisotropic critical state theory guided and validated by DEM

Fabric, expressed by means of an evolving deviatoric fabric tensor F, plays a very important role in the anisotropic mechanical response of granular materials. The Anisotropic Critical State Theory (ACST) addresses fabric anisotropy by rendering dilatancy a function of F, in addition to other state variables. In this paper, 3D DEM is used to guide the specific grain-level definition of F, the formulation of its continuum evolution equation and its effect on anisotropic dilatancy within ACST. DEM provides stress-ratio and shear strain variations as input for ACST analytical calculations of evolving fabric tensor and dilatancy, which are then favourably compared with totally independent direct measurements of these quantities by DEM. Furthermore, dilatancy is shown to be strongly affected by the contact normal-based fabric tensor Fc, whose evolution is best described by a continuum equation within ACST that includes dilatancy and a quantity related to particle orientation-based fabric tensor Fp. The aforementioned favourable comparison of the results for fabric tensor and dilatancy obtained independently by ACST and DEM, confirms the validity of the core framework of ACST irrespective of any constitutive model that addresses the deviatoric stress-strain relations.

58 GEOSCIENCES↗
DOE OSTI2019

ICAT: A numerical scheme to minimize numerical diffusion in advection-dispersion modeling and its application in identifying flow channeling

The advection-dispersion equation for scalar transport is essential for the numerical modeling of many fluid dynamics problems. However, solutions from numerical schemes always suffer from numerical diffusion and/or oscillation. In this study, we develop an Intra-Cell Advection Tracking (ICAT) scheme to minimize numerical diffusion and preserve monotonicity for advection-dispersion modeling. The key idea is to introduce “queues” in each discretized cell, and using a sequential transport rule and a flow distribution mechanism to track the scalar transport in these queues temporally and spatially. The capability and limitations of ICAT are first investigated through three test cases. Compared with the results obtained from other numerical schemes, the results from ICAT show substantially reduced numerical diffusion and agree better with analytical solutions. We also employ ICAT to simulate the transport process of a conservative tracer in a fracture with a highly heterogeneous aperture distribution. We report discrete flow channels in the fracture are better discerned by ICAT than by other numerical schemes, indicating the suitability of ICAT for modeling tracer transport in channelized flow fields.

15 GEOTHERMAL ENERGY↗
DOE OSTI2019

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

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

15 GEOTHERMAL ENERGY↗