SEARCH · Search NASA
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.
Numerical modeling of fluid and proppant transport. Shows impact of transport rate and proppant size on the extent of fracture coverage by proppant.
This study presents the results from in-well hydraulic fracture monitoring within a horizontal well in an unconventional reservoir utilizing Distributed Fiber Optic Sensing (DFOS). An in-house-developed Brillouin-based Distributed Strain Sensing (DSS) interrogator was deployed to obtain strain measurements, complemented by a commercial Raman-based Distributed Temperature Sensing (DTS) interrogator for temperature measurements and a commercial Rayleigh-based Low-Frequency Distributed Acoustic Sensing (LF-DAS) interrogator for strain-rate measurements. Examined over a ten-day period, the spatio-temporal distribution of temperature-compensated strain obtained from DSS and DTS revealed distinct signatures of the multi-stage hydraulic fracturing process. These signatures were analyzed with respect to fracture width growth and closure, residual strain effects, and fracture conductivity near the wellbore. Fracture widths within the fracture zone were estimated for individual stages. The findings were assessed with LF-DAS measurements for further evaluation. This work integrates DFOS-measured strain, temperature, and strain-rate data for monitoring in-well hydraulic fracturing, with the goal of supporting future studies in interpreting DFOS measurements for improved understanding of hydraulic fracturing in unconventional reservoirs.
The objective of this study is to understand how microseismic events are generated during hydraulic fracturing, as well as the role of geomechanical conditions (i.e., stress and mechanical stratigraphy) in this process. In the industry, microseismic event clouds have been generally used as an “outer-boundary” of the “stimulated reservoir volume” (SRV). However, by comparing with other surveillance data (Low frequency Distributed Acoustic Sensing, or LF-DAS strain) in the Hydraulic Fracturing Test Site (HFTS) 2 experiment, we show that this assumption is fundamentally flawed. The HFTS 2 data has three unique observations that have not been commonly observed in other datasets: 1. Due to influence of offset pad depletion, microseismic data shows that hydraulic fractures from the child well can propagate over 3000 feet into the depleted low stress zone. 2. By comparing microseismic and horizontal fiber LF-DAS strain data, we observe that microseismic event cloud does not necessarily reflect the created hydraulic fracture volume. Particularly, the extent of microseismic event clouds near heel stages are much shorter than what is shown with LF-DAS strain data. 3. Microseismic event magnitudes are larger in the depleted regions. Through geomechanical analysis, we demonstrate that the “bedding-plane-slip” model is likely the mechanism for microseismic generation during hydraulic fracturing. Furthermore, this model successfully explains the above field observations from HFTS 2 experiment. We also provide a quantitative relationship connecting the microseismic event magnitude with fracture width increment and layer mechanical property contrast.
We investigated the relationship between stratigraphy, stress, and microseismicity at the Hydraulic Fracture Test Site-1. The site comprises two sets of horizontal wells in the Wolfcamp shale and a deviated well drilled after hydraulic fracturing. Regional stresses indicate normal/strike-slip faulting with E-W compression. Stress measurements in vertical and horizontal wells show that the minimum principal stress varies with depth. Strata with high clay and organic content show high values of the least compressive stress, consistent with the theory of viscous stress relaxation. By integrating data from core, logs, and the hydraulic fracturing stages, we constructed a stress profile for the Wolfcamp sequence, which predicts how much pressure is required for hydraulic fracture growth. We applied the results to fracture orientation data from image logs to determine the population of pre-existing faults that are expected to slip during stimulation. We also determined microseismic focal plane mechanisms and found slip on steeply dipping planes striking NW, consistent with the orientations of potentially active faults predicted by the stress model. This case study represents a general approach for integrating stress measurements and rock properties to predict hydraulic fracture growth and the characteristics of injection-induced microseismicity.
Our team has conducted electromagnetic (EM) surveys for the past six years to monitor hydraulic-fracture behavior at the Devine Fracture Pilot Site (DFPS). The sub-horizontal orientation of a shallow hydraulic fracture at the DFPS provides uniform access to the fracture area for interrogation and data collection. Ahmadian et al. (2023) suggested a possible correlation between spatiotemporal changes in the flow rate, bottomhole pressure (BHP), and the observed surface recorded electric field at the DFPS. In this paper, we present the development of poroelastic forward models and pressure transient analyses (PTAs) to support the development of a multiphysics inverse model for these EM surveys. First, we conducted PTAs of the shut-in periods after six injections out of 10 to determine the fracture closure pressure (FCP) or the overburden pressure used in a poroelastic fracture reopening model. Second, we developed a finite-element poroelastic model throughout five injection cycles to include the effect of the cumulative injected volumes due to the previous injections on current fracture dilation in the presence of highly permeable unpropped and propped zones adjacent to the cohesive layer that models fracture reopening. Fracture reopening in this poroelastic model is based on a calibrated traction-separation response using the bottomhole pressure collected in two injection campaigns in 2020 and 2022. We used the outcomes of a previous simulation study of the primary hydraulic-fracturing stimulation to define the dimension of an unpropped fracture zone ahead of the propped fracture area. The PTAs led to FCPs consistent with those obtained using the injection data collected at the DFPS in 2020. Further, these analyses showed that at later injections, the fracture closure occurred at a later time with respect to the shut-in time, inferring the effect of cumulative injected volumes in previous injections. The simulation results show that considering the propped and unpropped fracture zones improves our poroelastic model in predicting the injection-well BHP. The numerical simulation results demonstrate a significant excess pore pressure near the fracture because of the preceding formation loadings by the previous injections. The obtained fracture dilation area and fluid pressure distribution provide a basis to improve the development of a multiphysics inverse model. Furthermore, in an iteratively coupled scheme, this pressure distribution can be introduced into EM models to render a holistic view of the causative mechanisms for the surface signal anomalies.
Enhanced geothermal systems can provide a substantial share of the global energy demand. There exist several hurdles in the engineering implementations of such geothermal systems. One such hurdle is the accurate monitoring of the fracture networks created in subsurface through hydraulic stimulation of these systems. Micro seismicity associated with the stimulation is the primary means to locate the event hypocenters for estimating the stimulated rock volume. Existing methods for location the hypocenters are restricted to only the highest amplitude impulsive signals that are simultaneously detected on several sensors. Consequently, a large portion (usually ~99%) of the measurements are left unused. In this paper, an unsupervised manifold-approximation followed by clustering of 3-component accelerometer data is used to analyze the seismicity recorded on a monitoring well. With this method, a larger portion of the measured signal is used for the monitoring of the hydraulic fracture network. We analyze the EGS Collab experiment 1 microseismic data, recorded at the Sanford Underground Research Facility, South Dakota. Using the data from a single three-component accelerometer, the polarization features viz. Azimuth, incidence, rectilinearity, and planarity are used as inputs for the unsupervised manifold approximation followed by clustering. Our study shows that density-based clusters in the projected 3D space correspond to distinct types of hydraulically fractured zones around the injection point. Finally, we show that the temporal evolution of these clusters can be used to track fracture creation and propagation.
This study investigates the mineralization occurring at the interface between injected hydraulic fracturing fluids and steel production casing during the shut-in period of unconventional shale wells. Mild steel coupons were reacted with synthetic hydraulic fracturing fluids with varying base fluid and dissolved gas concentrations. Experiments were carried out in an autoclave system set at 50 °C and 2000 PSI for 2 and 14 days. Physical and mineralogical changes to the steel coupons’ surface were characterized using an FEI Quanta 600 field emission gun (FEG) scanning electron microscope (SEM) and a PANalytical XPert Pro diffractometer. Barite precipitation occurred ubiquitously during the first 48 hours in all experimental scenarios. Sulfate radicals created from degradation and reduction of ammonium persulfate breakers facilitate barite formation. Concomitant barite and corrosion related precipitates on the sample surfaces suggest iron hydroxy salts act as nucleation sites for mineral scale growth. These results suggest hydraulic fracturing fluids may have deleterious effects on the wellbore before production of the well begins.
Southwest Research Institute® (SwRI ® ), Schlumberger Technology Corporation (SLB), and Chevron Corporation ® (Chevron) developed a novel, optimized, and lightweight process for natural gas (NG) to replace water as a low-cost fracturing medium with a low environmental impact. Hydraulic fracturing is used to increase oil and NG production by injecting high-pressure fluid, primarily water, into a rock formation, which fractures the rock and releases trapped oil and NG. This method was developed to increase yield and make feasible production areas that would not otherwise be viable for large-scale oil and NG extraction using traditional drilling technologies. Since the fracturing fluid is composed of approximately 90% water, one of the principal drawbacks to hydraulic fracturing is its excessive water use and associated large environmental footprint. According to recent data collected at fracture sites within the United States, fracturing applications in North America can consume as much as 9 million gallons of water per well. During the fracturing process, some of the fracturing fluid is permanently lost and the portion that is recovered is contaminated by both fracturing chemicals and dissolved solids from the formation. The recovered water or flow-back represents a significant environmental challenge, as it must be treated before it can be reintroduced into the natural water system. Although there is some recycling for future fracturing, the majority of the flow-back water is hauled from the well site to a treatment facility or to an injection well for permanent underground disposal. To mitigate these issues, an optimized, lightweight and modular surface process using NG to replace a majority of the water was developed as a cost-effective and environmentally clean fracturing fluid. Using NG will result in significantly less consumption since the gas that is injected as a fracturing fluid will be mixed with the formation gas and extracted as if it were from the formation itself. This process will minimize the collection, waste, and treatment of large amounts of water and will reduce the environmental impact of transporting and storing the fracturing fluid. There are two major steps involved in utilizing NG as the primary fracturing medium: (1) increasing the supply pressure of NG to wellhead pressures suitable for fracturing and (2) mixing the required chemicals and proppant needed for the fracturing process at these elevated pressures. The second step (NG-proppant mixing at elevated pressures) still requires demonstration; but very similar processes have been demonstrated in the field with other gases such as nitrogen (N 2 ) and carbon dioxide (CO 2 ). However, the first step (a compact, on-site unit for generating high-pressure NG at costs feasible for fracturing) has not been developed and is currently not commercially available. Due to the inherent compressibility of NG, more energy is required to compress the gas than what is required to pump water (or other incompressible liquids) to the very high pressure required for downhole injection. This project aimed to develop a novel, hybrid method to overcome this challenge. The project accomplishments and findings are discussed in this report. Ultimately, the research and development efforts described herein demonstrate that fracturing with NG foam is a feasible alternative to using water.
Understanding fracture propagation behavior is essential for optimizing hydraulic fracturing in unconventional reservoirs. This study demonstrates the value of integrating Low-Frequency Distributed Acoustic Sensing (LF-DAS) and microseismic data, which together provide a more complete picture of fracture growth. Using data from Hydraulic Fracturing Test Site 2 (HFTS2), we identify stress changes in depletion zones induced by parent wells as a key factor influencing fracture propagation. This result is shown by new measurements of in-situ fracture propagation velocity and fracture-hit volume (fluid volume at fracture hit?) from LF-DAS and event density from microseismic. These findings highlight the importance of considering parent well effects, well spacing, and stimulation sequencing in completion design to improve reservoir development and production efficiency.
Understanding fracture height growth can be of great significance to optimizing field development and improving recovery. The Hydraulic Fracturing Test Site 2 (HFTS-2) has provided a unique opportunity and an advanced dataset to allow us to observe and understand fracture geometries rigorously. Low frequency distributed acoustic sensing (LF-DAS) data from a vertical well in HFTS-2 showed three key observations: (i) excessive upward height growth (>1000 ft) and limited downward growth of the hydraulic fractures during pumping, (ii) considerable additional upward fracture height growth (~300 ft) after well shut in, and (iii) very complex LF-DAS strain rate patterns for a small fiber-to-stage offset. Advanced geomechanical modeling was performed to simulate the hydraulic fracture propagation and the resulting strain responses in the vertical direction. The modeling results demonstrated asymmetric upward and downward fracture height growths as observed in HFTS-2 with a similar upward height growth rate. Simulated waterfall plots of vertical strain rate showed distinct patterns for different fiber-to-fracture distances. The upward-growing fracture tip can be clearly identified by the interfaces between compressing and extending zones. It was also found that the complex strain rate patterns observed in HFTS-2 for small fiber-stage offsets were not caused by the mechanical layering but possibly result from the simultaneous propagation of multiple hydraulic fractures at different rates. Furthermore, the simulation results improved the understating of the HFTS-2 LF-DAS data and the simulated strain rate patterns could also serve as templates for fracture height interpretation from LF-DAS data in future.
Hydraulic fractures (HF) play a dominant role in oil and gas production from unconventional reservoirs, fracture characteristics, geometry, and extent are generally not well known and are typically inferred from various monitoring methods or subsurface models. The Hydraulic Fracturing Test Site-2 (HFTS-2) is a field-based research experiment performed in the Wolfcamp formation of the Permian (Delaware) Basin. HFTS-2, a unique joint industry program (JIP) that brings together government, academia, and industry, focused on field data acquisition and interpretation. The project has built a comprehensive HF research dataset that includes microseismic (MS), 1,500 ft of core, diagnostic formation injection tests (DFIT), fluid characterization, advanced logs, geochemistry, downhole gauges, varying completions, and permanent fiber optic (FO) cable in vertical and horizontal wells. These data were integrated and used to characterize HF geometry/dimensions and resulting depletion profiles, gain insights for completion design, and achieve a better understanding of parent well depletion effects on child well performance. The analysis and integration of the comprehensive HFTS-2 diagnostics dataset provide insights for optimal well placement and completions and improves understanding about the effects of parent well depletion. Furthermore, this dataset can be used to blind test and develop subsurface workflows for unconventional resources. Key learnings from HFTS-2 and the new technologies tested can easily be transferred to other unconventional plays.
Slickwater hydraulic fracturing is becoming a prevalent approach to economically recovering shale hydrocarbon. It is very important to understand the proppant’s transport behavior during slickwater hydraulic fracturing treatment for effective creation of a desired propped fracture geometry. The currently available models are either oversimplified or have been performed at limited length scales to avoid high computational requirements. Another limitation is that the currently available hydraulic fracturing simulators are developed using only single-sized proppant particles. Motivated by this, in this work, a computationally efficient, three-dimensional, multiphase particle-in-cell (MP-PIC) model was employed to simulate the multi-size proppant transport in a field-scale geometry using the Eulerian–Lagrangian framework. Instead of tracking each particle, groups of particles (called parcels) are tracked, which allows one to simulate the proppant transport in field-scale geometries at an affordable computational cost. Then, we found from our sensitivity study that pumping schedules significantly affect propped fracture surface area and average fracture conductivity, thereby influencing shale gas production. Motivated by these results, we propose an optimization framework using the MP-PIC model to design the multi-size proppant pumping schedule that maximizes shale gas production from unconventional reservoirs for given fracturing resources.
Hydraulic Fracturing Test Site-2 (HFTS-2) is a field-based research experiment performed in the Permian (Delaware) Basin. The unique aspect of this program was the acquisition of a unique, comprehensive, diagnostic dataset. Additionally, shorter parent wells drilled three years before the child wells offered clear distinction between the stages influenced by parent-child effects and the stages without any effects. The goal of this study was to analyze and integrate this comprehensive diagnostic dataset to understand the areal and vertical extent of hydraulic fractures (HF). The paper also provides insights on the effects of parent wells’ depletion on child well HF geometry based on various monitoring methods and subsurface models. Areal and vertical coverage for all HFTS-2 wells during stimulation and depletion was estimated based on analysis and interpretation of diagnostics and advanced modeling results. HFTS-2 diagnostics included microseismic (MS), pre- and post-stimulation logs and cores, bottomhole gauges, and fiber optic (FO) data. The diagnostics results (MS, FO, image logs) were integrated and used to calibrate subsurface models. Additional field tests were designed and implemented for depletion monitoring. The tailored program for monitoring depletion included vertical and slant well pressures, interference testing, and a vertical strain depletion trial. Areal Coverage: Conventional MS (and FO MS) were used to compute HF dimensions, which were compared with diagnostics (FO strain, gauge, image logs) observations and calibrated subsurface models. A post-production interference test did not show offset well communication. Vertical Coverage: Vertical coverage during stimulation was monitored using a vertical monitoring well. The stronger mechanical strain signals showed good correlation with MS event intensities, geomechanical properties, and gauge inferences. Vertical depletion was estimated based on vertical/slant well gauges and strain depletion tests. Parent-Child Effects: Diagnostics and calibrated subsurface models show asymmetry in child well fracture geometries for stages that overlap parent wells. Child well image logs serve as a good indicator for parent Downloaded from http://onepetro.org/URTECONF/proceedings-pdf/21URTC/2-21URTC/D021S031R004/2477423/urtec-2021-5241-ms.pdf/1 by Carol Worster on 28 February 2022 URTeC 5241 well HF tracking. Child well MS events had an eastward bias, in line with pre-stimulation image logs, and was confirmed by parent well frac hits. Novel/Additive Information: The dataset presents a unique, over-constrained problem space to compare independent techniques to arrive at HF metrics (i.e., stimulation height and/or half-length), unlike a single source dataset, in which calibration is done using available data to guide predictions. Here, the asymmetry in HF geometry seen in the stages influenced by parent-child effects offers unique insights into well spacing and landing, which are key capital decisions the unconventional resources industry is seeking to optimize.
This paper describes a new modeling framework for microscopic to reservoir-scale simulations of hydraulic fracturing and production. The approach builds upon a fusion of two existing high-performance simulators for reservoir-scale behavior: the GEOS code for hydromechanical evolution during stimulation and the TOUGH+ code for multi-phase flow during production. The reservoir-scale simulations are informed by experimental and modeling studies at the laboratory scale to incorporate important micro-scale mechanical processes and chemical reactions occurring within the fractures, the shale matrix, and at the fracture-fluid interfaces. These processes include, among others, changes in stimulated fracture permeability as a result of proppant behavior rearrangement or embedment, or mineral scale precipitation within pores and microfractures, at µm to cm scales. In our new modeling framework, such micro-scale testing and modeling provides upscaled hydromechanical parameters for the reservoir scale models. We are currently testing the new modeling framework using field data and core samples from the Hydraulic Fracturing Field Test (HFTS), a recent field-based joint research experiment with intense monitoring of hydraulic fracturing and shale production in the Wolfcamp Formation in the Permian Basin (USA). Below, we present our approach coupling the reservoir simulators GEOS and TOUGH+ informed by upscaled parameters from micro-scale experiments and modeling. We provide a brief overview of the HFTS and the available field data, and then discuss the ongoing application of our new workflow to the HFTS data set.
Hydraulic fracturing of deep shale formations generates large volumes of wastewater that must be managed through treatment, reuse, or disposal. Produced wastewater liberates formation-derived radionuclides and contains previously uncharacterized organohalides thought to be generated within the shale well, both posing unknown toxicity to human and ecological health. Here, we assess the toxicity of 42 input media and produced fluid samples collected from four wells in the Utica formation and Marcellus Shale using two distinct endpoint screening assays. Broad spectrum acute toxicity was assessed using a bioluminescence inhibition assay employing the halotolerant bacterium Aliivibrio fischeri, while predictive mammalian cytotoxicity was evaluated using a N-acetylcysteine (NAC) thiol reactivity assay. The acute toxicity and thiol reactivity of early-stage flowback was higher than later produced fluids, with levels diminishing through time as the natural gas wells matured. Acute toxicity of early stage flowback and drilling muds were on par with the positive control, 3,5-dichlorophenol (6.8 mg L -1 ). Differences in both acute toxicity and thiol reactivity between paired natural gas well samples were associated with specific chemical additives. Samples from wells containing a larger diversity and concentration of organic additives resulted in higher acute toxicity, while samples from a well applying a higher composition of ammonium persulfate, a strong oxidizer, showed greater thiol reactivity, predictive of higher mammalian toxicity. Both acute toxicity and thiol reactivity are consistently detected in produced waters, in some cases present up to nine months after hydraulic fracturing. These results support that specific chemical additives, the reactions generated by the additives, or the constituents liberated from the formation by the additives contribute to the toxicity of hydraulic fracturing produced waters and reinforces the need for careful consideration of early produced fluid management.
Abstract We present a multi-resolution approach for constructing model-based simulations of hydraulic fracturing, wherein flow through porous media is coupled with fluid-driven fracture. The approach consists of a hybrid scheme that couples a discrete crack representation in a global domain to a phase-field representation in a local subdomain near the crack tip. The multi-resolution approach addresses issues such as the computational expense of accurate hydraulic fracture simulations and the difficulties associated with reconstructing crack apertures from diffuse fracture representations. In the global domain, a coupled system of equations for displacements and pressures is considered. The crack geometry is assumed to be fixed and the displacement field is enriched with discontinuous functions. Around the crack tips in the local subdomains, phase-field sub-problems are instantiated on the fly to propagate fractures in arbitrary, mesh independent directions. The governing equations and fields in the global and local domains are approximated using a combination of finite-volume and finite element discretizations. The efficacy of the method is illustrated through various benchmark problems in hydraulic fracturing, as well as a new study of fluid-driven crack growth around a stiff inclusion.
In this work, we present a domain-based algorithm to simulate the propagation of a plane-strain hydraulic fracture in a zero-toughness permeable elastic medium. The algorithm utilizes a domain-based method to solve the elasticity equation and integrates a multi-scale tip asymptote, which is particular to hydraulic fractures, into this framework. This integration is key to accurately model the energy dissipation and the fluid leak-off in the fracture tip region. The algorithm combines a 2D finite volume method (FVM) for solving the elasticity equation with a 1D FVM for solving the nonlinear lubrication equation. Incorporating the far-field asymptotics and using a moving-mesh scheme reduces the computational burden while improving the accuracy of the scheme. The paper concludes with an analysis of the numerical results. This study demonstrates the potential of this domain-based approach for modeling hydraulic fractures in poroelastic media.
Bacteria remodel their plasma membrane lipidome to maintain key biophysical attributes in response to ecological disturbances. For Halanaerobium and other anaerobic halotolerant taxa that persist in hydraulically fractured deep subsurface shale reservoirs, salinity, and hydraulic retention time (HRT) are important perturbants of cell membrane structure, yet their effects remain poorly understood. Membrane-linked activities underlie in situ microbial growth kinetics and physiologies which drive biogeochemical reactions in engineered subsurface systems. Hence, we used gas chromatography–mass spectrometry (GC–MS) to investigate the effects of salinity and HRT on the phospholipid fatty acid composition of H. congolense WG10 and mixed enrichment cultures from hydraulically fractured shale wells. We also coupled acyl chain remodeling to membrane mechanics by measuring bilayer elasticity using atomic force microscopy (AFM). For these experiments, cultures were grown in a chemostat vessel operated in continuous flow mode under strict anoxia and constant stirring. Our findings show that salinity and HRT induce significant changes in membrane fatty acid chemistry of H. congolense WG10 in distinct and complementary ways. Notably, under nonoptimal salt concentrations (7% and 20% NaCl), H. congolense WG10 elevates the portion of polyunsaturated fatty acids (PUFAs) in its membrane, and this results in an apparent increase in fluidity (homeoviscous adaptation principle) and thickness. Double bond index (DBI) and mean chain length (MCL) were used as proxies for membrane fluidity and thickness, respectively. These results provide new insight into our understanding of how environmental and engineered factors might disrupt the physical and biogeochemical equilibria of fractured shale by inducing physiologically relevant changes in the membrane fatty acid chemistry of persistent microbial taxa. GRAPHICAL ABSTRACT Salinity significantly alters membrane bilayer fluidity and thickness in Halanaerobium congolense WG10.