Seismic monitoring of the growth of a hydraulic fracture zone at Fenton Hill, New Mexico

Geophysics ◽  
1998 ◽  
Vol 63 (1) ◽  
pp. 120-131 ◽  
Author(s):  
Yingping Li ◽  
Chuen H. Cheng ◽  
M. Nafi Toksöz

The hydraulic fracturing technique is an important method for enhancing hydrocarbon recovery, geothermal energy extraction, and solid waste disposal. Determination of the geometry and growth process of a hydraulic fracture zone is important for monitoring and assessing subsurface fractures. A relative‐source‐location approach, based on a waveform correlation and a grid search method, has been developed to estimate relative hypocenter locations for a cluster of 157 microearthquakes induced by hydraulic fracturing at the Los Alamos Hot Dry Rock (HDR) geothermal site. Among the 157 events, 147 microearthquakes occurred in a tight cluster with a dimension of 40 m, roughly defining a vertical hydraulic fracture zone with an orientation of N40°W. The length, height, and width of the hydraulic fracture zone are estimated to be 40, 35, and 5 m, respectively. Analysis of the spatial‐temporal pattern of the induced microearthquakes reveals that the fracture zone grew significantly, averaging 0.2 m/minute in a two‐hour period toward the northwest along the fracture zone strike.

Geofluids ◽  
2021 ◽  
Vol 2021 ◽  
pp. 1-24
Author(s):  
Xin Zhang ◽  
Yuqi Zhang ◽  
Bingxiang Huang

Hydraulic fracturing applications have shown a stress disturbance effect during hydraulic fracture propagation, which is often ignored. Using laboratory and discrete element numerical simulation tests, hydraulic fracture propagation under this stress disturbance is systematically studied. The results show that during hydraulic fracturing, the bedding plane is damaged by the stress disturbance, forming a bedding fracture zone (BFZ). The nonlinear fracture characteristics of the formation process of the disturbed fracture zone are revealed, and two indexes (the number of fractures in the disturbed fracture zone and the size of the disturbed fracture zone) are proposed to evaluate the fracturing effect of the stress disturbance. Based on these indexes, multifactor sensitivity tests are conducted under different geological conditions and operational factors. When the principal stress ( σ 1 ) difference is large, the number of hydraulic fractures gradually decreases from many to one, and the direction of the hydraulic fractures gradually approaches the vertical direction of σ 3 , but the change in the in situ stress condition has no obvious effect on the stress disturbance effect. The weaker the bonding strength of the bedding plane, the more significant the stress disturbance effect is, and the easier it is for the fractures to expand along the bedding plane. With increasing injection rate, the stress disturbance effect first increases and then decreases, and the hydraulic fracture propagates from along the bedding plane to cross the bedding plane. With increasing relative distance between the injection hole and bedding plane, the stress disturbance effect presents a linearly increasing trend, and the hydraulic fractures along the bedding planes extend. Based on the experimental results, the relationship between the fracturing effect of the stress disturbance and the extension mode of the hydraulic fracture is determined, and an optimization method for hydraulic fracturing in composite rock reservoirs is given. The research results can provide a theoretical basis for controlling the formation of complex fracture networks in composite rock reservoirs.


1998 ◽  
Vol 38 (1) ◽  
pp. 830 ◽  
Author(s):  
S.P. Narayan ◽  
D. Naseby ◽  
Z. Yang ◽  
S.S. Rahman

The Cooper Basin is the largest gas-producing basin in Australia and hosts a huge volume of natural gas in low permeability (known as 'tight gas') sandstone formations. Hydraulic fracture treatments, based on 'opening mode' fracture mechanics, are routinely carried out to unlock tight gas and to accelerate its recovery. Information regarding insitu stresses and natural fractures is required for successful fracture treatments. However, hydraulic fracturing is still often problematic, in part due to the relatively high insitu stresses and temperatures in the region. A vast amount of Hot Dry Rock (HDR) geothermal energy resources exists in granites below the sedimentary rocks in the Cooper Basin. Exploitation of HDR requires the same drilling and completion technologies as used in the petroleum industry. Hydraulic fracturing is also necessary for HDR reservoir creation, and requires characterisation of insitu stresses and natural fractures, as does tight gas production. It has been realised that the mechanism for reservoir stimulation in granitic rocks is proppant free shear dilation that is related to 'sliding mode' fracture mechanics. Furthermore, seismic imaging of hydraulic fracture propagation is well established in the HDR industry. These two technologies, developed in HDR, may have potential application to the petroleum industry for tight gas production. The geographic proximity of tight gas and HDR geothermal energy in the Cooper Basin and common exploitation technologies should justify close collaboration between the petroleum industry and HDR researchers.


Geophysics ◽  
1994 ◽  
Vol 59 (6) ◽  
pp. 881-888 ◽  
Author(s):  
Said Attia al Hagrey

Electric and electromagnetic methods have been applied for mapping subsurface fractures and the directional dependence of in‐situ electric parameters at the hot dry rock site at Falkenberg, Germany. This study includes the determination of several anisotropy parameters like the mean, longitudinal and transverse resistivity components ([Formula: see text], [Formula: see text] and [Formula: see text], respectively), the anisotropy coefficient λ, and the strike angle Θ. Terrain conductivity measurements using the technique of frequency‐domain electromagnetic induction reveal a dominant anomaly strike of east‐south‐east— west‐north‐west, nearly parallel to the fracturing strike of N110°. With increasing distance from the central borehole, the mise à la masse potential differences exhibit a transition from a direct to a paradoxical relationship to the resistivity anisotropy induced by the fracturing. These observations are explained using a model for an ellipsoidal fracture. The qualitative interpretation of the sounding data of Schlumberger and crossed‐square arrays clearly shows the anisotropy paradox related to the N110° strike. The crossed‐square method applied over a range of electrode spacings of 2.1–447 m yields apparent anisotropy values of [Formula: see text], [Formula: see text], [Formula: see text], λ = 1.15–1.34 and Θ ≈ 110°. The quantitative interpretation of the sounding data shows a conductive fracture zone (ρ ≈ 300–700 ω ⋅ m) embedded in the resistive granitic basement (ρ > 1200 ω ⋅ m) at a depth of 55–85 m. The anisotropy parameters for the fracture zone are ρm = 700 ω ⋅ m, [Formula: see text], ρt = 1225 ω ⋅ m and λ = 1.75. The fact that the anisotropy coefficient of the conductive zone is higher than the apparent λ is in agreement with the conclusion that the anisotropy at the study site is related mainly to fracturing.


2021 ◽  
pp. 014459872198899
Author(s):  
Weiyong Lu ◽  
Changchun He

Directional rupture is one of the most important and most common problems related to rock breaking. The goal of directional rock breaking can be effectively achieved via multi-hole linear co-directional hydraulic fracturing. In this paper, the XSite software was utilized to verify the experimental results of multi-hole linear co-directional hydraulic fracturing., and its basic law is studied. The results indicate that the process of multi-hole linear co-directional hydraulic fracturing can be divided into four stages: water injection boost, hydraulic fracture initiation, and the unstable and stable propagation of hydraulic fracture. The stable expansion stage lasts longer and produces more microcracks than the unstable expansion stage. Due to the existence of the borehole-sealing device, the three-dimensional hydraulic fracture first initiates and expands along the axial direction in the bare borehole section, then extends along the axial direction in the non-bare hole section and finally expands along the axial direction in the rock mass without the borehole. The network formed by hydraulic fracture in rock is not a pure plane, but rather a curved spatial surface. The curved spatial surface passes through both the centre of the borehole and the axial direction relative to the borehole. Due to the boundary effect, the curved spatial surface goes toward the plane in which the maximum principal stress occurs. The local ground stress field is changed due to the initiation and propagation of hydraulic fractures. The propagation direction of the fractures between the fracturing boreholes will be deflected. A fracture propagation pressure that is greater than the minimum principle stress and a tension field that is induced in the leading edge of the fracture end, will aid to fracture intersection; as a result, the possibility of connecting the boreholes will increase.


2020 ◽  
Vol 10 (8) ◽  
pp. 3333-3345
Author(s):  
Ali Al-Rubaie ◽  
Hisham Khaled Ben Mahmud

Abstract All reservoirs are fractured to some degree. Depending on the density, dimension, orientation and the cementation of natural fractures and the location where the hydraulic fracturing is done, preexisting natural fractures can impact hydraulic fracture propagation and the associated flow capacity. Understanding the interactions between hydraulic fracture and natural fractures is crucial in estimating fracture complexity, stimulated reservoir volume, drained reservoir volume and completion efficiency. However, because of the presence of natural fractures with diffuse penetration and different orientations, the operation is complicated in naturally fractured gas reservoirs. For this purpose, two numerical methods are proposed for simulating the hydraulic fracture in a naturally fractured gas reservoir. However, what hydraulic fracture looks like in the subsurface, especially in unconventional reservoirs, remain elusive, and many times, field observations contradict our common beliefs. In this study, the hydraulic fracture model is considered in terms of the state of tensions, on the interaction between the hydraulic fracture and the natural fracture (45°), and the effect of length and height of hydraulic fracture developed and how to distribute induced stress around the well. In order to determine the direction in which the hydraulic fracture is formed strikethrough, the finite difference method and the individual element for numerical solution are used and simulated. The results indicate that the optimum hydraulic fracture time was when the hydraulic fracture is able to connect natural fractures with large streams and connected to the well, and there is a fundamental difference between the tensile and shear opening. The analysis indicates that the growing hydraulic fracture, the tensile and shear stresses applied to the natural fracture.


2021 ◽  
pp. 014459872098153
Author(s):  
Yanzhi Hu ◽  
Xiao Li ◽  
Zhaobin Zhang ◽  
Jianming He ◽  
Guanfang Li

Hydraulic fracturing is one of the most important technologies for shale gas production. Complex hydraulic fracture networks can be stimulated in shale reservoirs due to the existence of numerous natural fractures. The prediction of the complex fracture network remains a difficult and challenging problem. This paper presents a fully coupled hydromechanical model for complex hydraulic fracture network propagation based on the discontinuous deformation analysis (DDA) method. In the proposed model, the fracture propagation and rock mass deformation are simulated under the framework of DDA, and the fluid flow within fractures is simulated using lubrication theory. In particular, the natural fracture network is considered by using the discrete fracture network (DFN) model. The proposed model is widely verified against several analytical and experimental results. All the numerical results show good agreement. Then, this model is applied to field-scale modeling of hydraulic fracturing in naturally fractured shale reservoirs. The simulation results show that the proposed model can capture the evolution process of complex hydraulic fracture networks. This work offers a feasible numerical tool for investigating hydraulic fracturing processes, which may be useful for optimizing the fracturing design of shale gas reservoirs.


2021 ◽  
Author(s):  
A. Kirby Nicholson ◽  
Robert C. Bachman ◽  
R. Yvonne Scherz ◽  
Robert V. Hawkes

Abstract Pressure and stage volume are the least expensive and most readily available data for diagnostic analysis of hydraulic fracturing operations. Case history data from the Midland Basin is used to demonstrate how high-quality, time-synchronized pressure measurements at a treatment and an offsetting shut-in producing well can provide the necessary input to calculate fracture geometries at both wells and estimate perforation cluster efficiency at the treatment well. No special wellbore monitoring equipment is required. In summary, the methods outlined in this paper quantifies fracture geometries as compared to the more general observations of Daneshy (2020) and Haustveit et al. (2020). Pressures collected in Diagnostic Fracture Injection Tests (DFITs), select toe-stage full-scale fracture treatments, and offset observation wells are used to demonstrate a simple workflow. The pressure data combined with Volume to First Response (Vfr) at the observation well is used to create a geometry model of fracture length, width, and height estimates at the treatment well as illustrated in Figure 1. The producing fracture length of the observation well is also determined. Pressure Transient Analysis (PTA) techniques, a Perkins-Kern-Nordgren (PKN) fracture propagation model and offset well Fracture Driven Interaction (FDI) pressures are used to quantify hydraulic fracture dimensions. The PTA-derived Farfield Fracture Extension Pressure, FFEP, concept was introduced in Nicholson et al. (2019) and is summarized in Appendix B of this paper. FFEP replaces Instantaneous Shut-In Pressure, ISIP, for use in net pressure calculations. FFEP is determined and utilized in both DFITs and full-scale fracture inter-stage fall-off data. The use of the Primary Pressure Derivative (PPD) to accurately identify FFEP simplifies and speeds up the analysis, allowing for real time treatment decisions. This new technique is called Rapid-PTA. Additionally, the plotted shape and gradient of the observation-well pressure response can identify whether FDI's are hydraulic or poroelastic before a fracture stage is completed and may be used to change stage volume on the fly. Figure 1Fracture Geometry Model with FDI Pressure Matching Case studies are presented showing the full workflow required to generate the fracture geometry model. The component inputs for the model are presented including a toe-stage DFIT, inter-stage pressure fall-off, and the FDI pressure build-up. We discuss how to optimize these hydraulic fractures in hindsight (look-back) and what might have been done in real time during the completion operations given this workflow and field-ready advanced data-handling capability. Hydraulic fracturing operations can be optimized in real time using new Rapid-PTA techniques for high quality pressure data collected on treating and observation wells. This process opens the door for more advanced geometry modeling and for rapid design changes to save costs and improve well productivity and ultimate recovery.


2021 ◽  
Author(s):  
Seyhan Emre Gorucu ◽  
Vijay Shrivastava ◽  
Long X. Nghiem

Abstract An existing equation-of-state compositional simulator is extended to include proppant transport. The simulator determines the final location of the proppant after fracture closure, which allows the computation of the permeability along the hydraulic fracture. The simulation then continues until the end of the production. During hydraulic fracturing, proppant is injected in the reservoir along with water and additives like polymers. Hydraulic fracture gets created due to change in stress caused by the high injection pressure. Once the fracture opens, the bulk slurry moves along the hydraulic fracture. Proppant moves at a different speed than the bulk slurry and sinks down by gravity. While the proppant flows along the fracture, some of the slurry leaks off into the matrix. As the fracture closes after injection stops, the proppant becomes immobile. The immobilized proppant prevents the fracture from closing and thus keeps the permeability of the fracture high. All the above phenomena are modelled effectively in this new implementation. Coupled geomechanics simulation is used to model opening and closure of the fracture following geomechanics criteria. Proppant retardation, gravitational settling and fluid leak-off are modeled with the appropriate equations. The propped fracture permeability is a function of the concentration of immobilized proppant. The developed proppant simulation feature is computationally stable and efficient. The time step size during the settling adapts to the settling velocity of the proppants. It is found that the final location of the proppants is highly dependent on its volumetric concentration and slurry viscosity due to retardation and settling effects. As the location and the concentration of the proppants determine the final fracture permeability, the additional feature is expected to correctly identify the stimulated region. In this paper, the theory and the model formulation are presented along with a few key examples. The simulation can be used to design and optimize the amount of proppant and additives, injection timing, pressure, and well parameters required for successful hydraulic fracturing.


Energies ◽  
2018 ◽  
Vol 11 (10) ◽  
pp. 2647
Author(s):  
Gang Wang ◽  
Cheng Fan ◽  
Hao Xu ◽  
Xuelin Liu ◽  
Rui Wang

Accurately determining the height of the gas-guiding fracture zone in the overlying strata of the goaf is the key to find the height of the long horizontal borehole in the roof. In order to determine the height, in this study we chose the 6306 working face of Tangkou Coal Mine in China as a research example and used both the theoretical model and discrete element method (DEM) numerical simulation to find the height of the gas-guiding fracture zone and applied the height to drill a long horizontal borehole in the roof of the 6303 working face. Furthermore, the borehole was utilized to deep into the roof for coalbed methane drainage and the results were compared with conventional gas drainage measures from other aspects. The height of the gas-guiding fracture zone was found to be 48.57 m in theoretical model based on the bulk coefficient and the void ratio and to be 51.19 m in the DEM numerical simulation according to the temporal and spatial variation characteristics of porosity. Taking both the results of theoretical analysis and numerical simulation into consideration, we determined that gas-guiding fracture zone is 49.88 m high and applied it to drill a long horizontal borehole deep into the roof in the 6303 working face field. Compared with conventional gas drainage measures, we found that the long horizontal borehole has the high stability, high efficiency and strong adaptability for methane drainage.


1985 ◽  
Vol 25 (01) ◽  
pp. 78-88 ◽  
Author(s):  
T.K. Perkins ◽  
J.A. Gonzalez

Abstract When a cool fluid such as water is injected into a hot reservoir, a growing region of cooled rock is established around the injection well. The rock matrix within the cooled region contracts, and a thermoelastic stress field is induced around the well. For typical waterflooding of a moderately deep reservoir, horizontal earth stresses may be reduced by several hundred psi. If the injection pressure is too high or if suspended solids in the water plug the formation face at the perforations, the formation will be fractured hydraulically. As the fracture grows, the flow system evolves from an essentially circular geometry in the plan view to one characterized more nearly as elliptical. This paper considers thermoelastic stresses that would result from cooled regions of fixed thickness and of elliptical cross section. The stresses for an infinitely thick reservoir have been deduced from information available in public literature. A numerical method has been developed to calculate thermoelastic stresses induced within elliptically shaped regions of finite thickness. Results of these two approaches were combined, and empirical equations were developed to give an approximate but convenient, explicit method for estimating induced stresses. An example problem is given that shows how this theory can be applied to calculate the fracture lengths, bottomhole pressures (BHP's), and elliptical shapes of the flood front as the injection process progresses. Introduction When fluids are injected into a well, such as during waterflooding or other secondary or tertiary recovery processes, the temperatures of the injected fluids are typically cooler than the in-situ reservoir temperatures. A region of cooled rock forms around each injection well, and this region grows as additional fluid is injected. Formation rock within the cooled region contracts, and this leads to a decrease in horizontal earth stress near the injection well. In Ref. 1, the magnitude of the reduction in horizontal earth stress was given for the case of a radially symmetrical cooled region. Another factor, which may occur simultaneously, is the plugging of formation rock by injected solids. There is extensive literature indicating that waters normally available for injection contain suspended solids. Laboratory tests demonstrate that these waters, when injected into formation rocks, can plug the face of the rock or severely limit injectivity. In field operations, injection often simply continues at a BHP that is high enough to initiate and extend hydraulic fractures." The injected fluid then can leak off readily through the large fracture face area. Because of the lowering of horizontal earth stresses that results from cold fluid injection, hydraulic fracturing pressures can be much lower than would be expected for an ordinary low-leakoff hydraulic fracturing treatment. For this reason, the well operator may not be aware that injected fluid is being distributed through an extensive hydraulic fracture. If injection conditions are such that a hydraulic fracture is created, then the flow system will evolve from an essentially circular geometry in the plan view to one characterized more nearly as elliptical. In this paper, thermoelastic stresses for cooled regions of fixed thickness and of elliptical cross section are determined, and a theory of hydraulic fracturing of injection wells is developed. Conditions under which secondary fractures (perpendicular to the primary, main fracture) will open also are discussed. Finally, an example problem is given to illustrate how this theory can be applied to calculate fracture lengths, BHP'S, and elliptical shapes of the flood front as the injection process progresses. Thermoelastic Stresses in Regions of Elliptical Cross Section If fluid of constant viscosity is injected into a line crack (representing a two-wing, vertical hydraulic fracture), the flood front will progress outward. so its outer boundary at any time can be described approximately as an ellipse that is confocal with the line crack. If the injected fluid is at a temperature different from the formation temperature, a region of changed rock temperature with fairly sharply defined boundaries will progress outward from the injection well but lag behind the flood front. The outer boundary of the region of changed temperature also will be elliptical in its plan view and confocal with the line crack (see Fig. 1). Stresses within the region of altered temperature, as well as stress in the surrounding rock, which remains at its initial temperature, will be changed because of the expansion or contraction of the rock within the region of altered temperature. The thermoelastic stresses within an infinitely tall cylinder of elliptical cross section can be determined from information available in the literature. 10 The interior thermoelastic stresses perpendicular and parallel to the major axes of the ellipse are given by Eqs. 1 and 2, respectively. SPEJ P. 78^


Sign in / Sign up

Export Citation Format

Share Document