Hydraulic Fracture Network Research Articles

The hydraulic fracture propagation pattern induced by multi-stage temporary plugging and diverting fracturing in reservoirs with various lithologies: An experimental investigation

Multi-stage temporary plugging and diverting fracturing (TPDF) is an effective method for generating hydraulic fracture (HF) networks. This study investigates various lithological reservoirs in the Xinjiang region, obtaining downhole full-diameter cores for experimental analysis using true triaxial TPDF. The characteristics of HF morphology are quantitatively assessed by employing computed tomography (CT) scanning. The findings are summarized as follows: (1) Initial hydraulic fracturing of specimens with different lithologies results in σH-direction double-wing HF, while the first TPDF generates a single-wing HF along the σh direction, and the second TPDF produces a single-wing HF along the σH direction. (2) The volume and area of HFs in the first TPDF of volcanic rock increased by over 30%. The first TPDF effect is more pronounced in conglomerate rock, with HF volume over 25% and surface area increasing by more than 35%. (3) During multi-stage TPDF, volcanic rock transitions from initial HF to the formation of new HF, sandstone diverts from the wellbore to create new HF, and conglomerate generates new HF through multi-point initiation in the wellbore and HF. Each TPDF process involves the propagation of existing HFs and the generation of new ones. (4) The breakdown pressure in multi-stage TPDF increased by 46.5% and 51.6% in volcanic rock, while in sandstone, the first TPDF increased by 90.6%. In conglomerate rock, multi-stage TPDF saw increases of 51.2% and 41.9%, respectively. These findings offer theoretical insights for optimizing TPDF design in diverse lithological reservoirs.

Utilizing well-reservoir pseudo-connections for multi-stage hydraulic fracturing modeling in tight gas saturated formations

Purpose. Research is aimed at integrating multi-stage hydraulic fracturing in horizontal wells with hydrodynamic simulation as a mandatory part of planning the mining of any shale oil or gas reservoir. Methods. Geological and hydrodynamic reservoir modeling is part of the research. The properties and geometries of the hydraulic fracture network and its representation in the dynamic reservoir model were assessed. The comparative characterization was carried out based on the two methods of fracture modeling: cell dimension reduction for explicit fracture modeling (LGR – local grid refinement) and implicit fracture modeling method, presented in this paper, with additional pseudo-connections between well and reservoir. Findings. A hydrodynamic model for low-permeable reservoir, produced by horizontal well, hydraulically fractured with 5 stages, has been generated. This model is calibrated to the production history and flowing bottom hole pressure by applying two methods of fracture modeling. Modeling results show that it is possible to replicate historical well production by using both methods. However, the proposed method with pseudo connections has several advantages compared to the generally accepted, local grid refinement (LGR) method. Originality. For the first time, a system of pseudo connections between well and reservoir was constructed to model a multi-stage hydraulic fracturing for a hydrodynamic model of tight reservoir. Hydrodynamic simulation results were refined and calibrated to the history of hydrocarbon production and flowing bottom hole pressure data using the pseudo-connections and LGR methods. The similarity of the results by applying LGR and pseudo-connections methods was revealed. Practical implications. The use of pseudo connections for hydraulic fracturing modeling can reduce simulation run time for cases where multi-stage hydraulic fracturing has already been carried out or is planned in the future. Additionally, the use of this method allows testing a larger number of realizations and scenarios, including hydraulic fracturing design (number of stages, size and conductivity of resulted fracture systems, fracture orientation, etc.), well placement and fracture growth relative to well trajectory. Also, there is no need to rebuild a model every time for each realization, as is the case with the LGR method.

Effects of microstructure on hydraulic fracturing and gas–water production in coal reservoirs: A case study of the Dahebian coalbed methane block in Western Guizhou, China

AbstractThe development of coalbed methane (CBM) in China is susceptible to the influence of microstructure. Therefore, it is crucial to understand the extension laws of hydraulic fractures and engineering responses in coal reservoirs affected by microstructure development. Utilizing the Dahebian block in western Guizhou as the study area, this investigation examines the coal reservoir characteristics, fracturing, and drainage engineering analysis of the well DC1 group in the region. The aim is to discuss the spatial distribution characteristics of hydraulic fractures, geological controlling factors influenced by microstructure, and their corresponding engineering responses. The results indicate that, for coal reservoirs unaffected by microstructure, the extension laws of the fracture network in both longitudinal and planar directions are influenced by burial depth and the regional stress field. In the microstructural belt, tectonic stress dominates, causing changes in the ground stress field. Consequently, the hydraulic fracture network deviates from the direction of the maximum principal stress during the extension process. When a secondary fracture is nearby, the hydraulic fracture network extends towards the shortest path radial secondary fracture direction, leading to a rapid increase in fracture width per unit length until it intersects with the secondary fracture. Additionally, the presence of secondary joints near the microfault structure decreases fracturing pressure and results in a dense distribution of the fracture network. This promotes the formation of a complex fracture network, favorable for fracturing. The extension of the fracture network in complex structural development areas is influenced by the microfault structure between wells, which is reflected in the fracturing construction pressure and fluid output. This accounts for the significant variations in the early drainage performance of CBM wells.

Effects of saturated fluids on petrophysical properties of hot dry rock at high temperatures: An experimental study

Hydraulic fracturing is a crucial technology for efficient heat extraction from hot dry rock (HDR). The fracturing process involves the injection of large volumes of cold fluids, which can significantly affect the petrophysical properties of the rock, particularly when temperature differences and fluid-rock interactions are taken into account. However, there is a lack of research on the evolution of HDR petrophysical properties under the influence of different fluids at high temperatures. In this paper, we utilized a high-temperature high-pressure tri-axial loading system, rock spontaneous imbibition instrument, porosity-permeability tester, and thermal conductivity meter, along with comprehensive analysis techniques such as ultrasonic wave velocity monitoring, CT scanning, and nuclear magnetic resonance (NMR), to systematically investigate the effects of distilled water and nano emulsion on the physical and mechanical properties of the HDR at high temperatures. Through comparison with the experimental results of dry samples, comprehensive characteristics of physical parameters for the HDR under the influence of distilled water and nano emulsion are obtained, such as spontaneous imbibition mass and permeability, thermal physical properties, and rock mechanical parameters. The experimental results show that temperature increase promotes the formation of massive thermal damage fractures in the HDR samples, and the fluid weakens the mechanical strength of the HDR. When the effects of temperature and saturated fluids are combined, the hydraulic fracture network in the HDR samples exhibits increased complexity, contrasting with the fracture patterns observed in the heat-treated dry samples. The distilled water has a more pronounced effect on weakening the mechanical strength of the samples than the nano emulsion. Moreover, the analyses of permeability, thermal conductivity, ultrasonic wave velocity, and NMR results are consistent with the results of mechanical experiments. As the temperature rises, the wave velocity of saturated distilled water samples exceeds that of both saturated nano emulsion samples and dry samples, the velocity of the saturated fluid front in the samples increases, and the range of the fluid in the pores expands. Overall, these findings provide valuable experimental evidence and theoretical guidance for better understanding the hydro-thermal-mechanical behaviors of the HDR, which are crucial for the efficient development and utilization of geothermal energy.

Study on geothermal energy self-recycling extraction and ScCO2 storage in the horizontal well annuli with fracture network system

With the extensive exploitation of fossil fuels, a substantial number of production and injection wells have been or are on the verge of being abandoned. However, many of these wells present viable opportunities for geothermal energy extraction, as they are situated in geothermal reservoirs with favourable temperatures. In alignment with carbon neutrality policies, supercritical CO2 (ScCO2) has emerged as the preferred working fluid for geothermal extraction, facilitating subsequent storage. This study has developed three transient models to characterize the heat and mass transfer processes of ScCO2 circulating within a horizontal wellbore and the reservoir's fracture network. Finite element method and Comsol Multiphysics 6.1 software are employed to solve these models. Results show that ScCO2 efficiently flows from the annuli into the hydraulic fracture network and then circulates back through the annuli completing a self-recycling heat extraction process. The findings reveal several key insights: (a) Lower injection mass flow rates and flow rates lead to more efficient extraction of thermal energy within horizontal wellbores. (b) Low-temperature ScCO2 enhance thermal energy recovery. (c) The larger the coverage area of the fracture network, the more effective the heat extraction. (d) The efficiency of heat extraction and ScCO2 storage in the horizontal well-fracture network self-circulation system is 3 to 5 times and 6 to 7 times higher than that of the horizontal well self-recycling and injection-production systems. (e) Periodic thermal extraction proves to be more efficient than continuous methods. (f) The cumulative heat production of the vertical dual well injection-production system is the highest, reaching 2.45 × 1011 kJ. But it's still lower than the self-circulating system. The heat extraction effect of the diagonal dual well injection-production system is the lowest, which is 1.83 × 1011 kJ. (g) Three developed transient models facilitate efficient heat extraction and ScCO2 storage reservoir screening and evaluation. In conclusion, the self-recycling system has achieved efficient heat extraction and carbon dioxide storage, providing a theoretical basis for future energy conversion and management.

Numerical Simulation of Fracturing Fluid Storage in Shale Reservoirs Based on Experimental Measurements of Stress Sensitivity of Hydraulic Fracture Network Conductivity

Hydraulic fracturing is used in shale reservoir production, with low flowback rates and a large amount of fracturing fluid retained inside the reservoir. In this study, a stress sensitivity analysis experiment on the fracture inflow capacity was implemented to investigate the relationship between the hydraulic fracture (HF) and natural fracture (NF) inflow capacities and effective stress. A three-dimensional shale reservoir model was also constructed to couple the experimentally obtained laws with the numerical model to investigate the effects of the connection and closure of the fracture network on the retention of the fracturing fluid. The results show that the stress sensitivity of natural fractures is two orders of magnitude higher than that of hydraulic fractures. The seepage-absorption effect of capillary forces is not the whole reason for the large amount of fracturing fluid retention. The closure of the fracture network formed by natural and hydraulic fractures during the production process led to the storage of a large amount of fracturing fluid, and this process maintained the stability of the water production rate during the steady water production period.

Numerical study of multiple hydraulic fractures propagation in poroelastic media based on energy decomposition phase field methods

This paper proposes a novel phase-filed model for simulating hydraulic fracturing in poroelastic media under complex stress conditions. The main theoretical contribution lies in the fact that a generalized strain energy density decomposition is incorporated to hydraulic fracturing phase field model, which unifies hydrostatic-deviatoric and spectral decompositions and predicts material failure through classical Mohr-Coulomb failure criterion. Based on Biot’s poroelastic theory, an indicator function is used to distinguish between fluid flow in fractures and in pores. The fully coupled equations are solved by combining staggered schemes with Newton-Raphson methods. The proposed model is validated by comparing laboratory experimental results with numerical results. In order to demonstrate the influence of introducing hydrostatic-deviatoric and spectral energy decompositions, we simulated propagation of multiple hydraulic fractures and coalescence of hydraulic fractures with natural fractures. Finally, the effects of stress ratios (σx/σy) on multiple hydraulic fractures are analysed. It is found that when the stress ratio is 2, a more complex hydraulic fracture network is formed in the porous media.

Experimental Investigation into the Process of Hydraulic Fracture Propagation and the Response of Acoustic Emissions in Fracture–Cavity Carbonate Reservoirs

Fracture–cavity carbonate reservoirs account for a considerable proportion of oil and gas resources. Because of the complicated relationships between cavities, fractures and pores in these reservoirs, which are defined as cavity clusters, fracturing technology is employed to enhance their hydrocarbon productivity. However, almost all previous studies have just considered the effect of a single natural cavity or fracture on the propagation of a hydraulic fracture; therefore, the mechanism by which a hydraulic fracture interacts with a cavity cluster needs to be clarified. In this study, cavity clusters with different distributions were accurately prefabricated in synthetically made samples, and large-scale simulation equipment was employed to systematically perform fracturing experiments considering different horizontal differential stress levels. Meanwhile, the hydraulic fracture propagation behaviors were comprehensively analyzed through fracture morphology, fracturing curves, the complexity of the fracture network and acoustic emission monitoring. It was found that a natural fracture with a smaller approach angle is favorable in guiding a hydraulic fracture to a cavity. The fracturing curves were divided into the following four types: frequent fluctuations with “step-like” shapes, great fluctuations with slightly lower closure pressure, fluctuations with obviously lower closure pressure, and little fluctuations with obviously lower closure pressure. And different cavity cluster distributions play a dominant role in the complexity of generated hydraulic fracture networks. In addition, AE energy was used to judge the ease of crossing the cavity. The above findings indicated that for the actual exploration and exploitation of carbonate reservoirs, the geological exploration of different fracture–cavity structures in reservoirs would be required, and targeted fracturing engineering designs need to be carried out for different fracture–cavity carbonate reservoirs.

Characterization of methane hydrate extraction influenced by hydraulic fractures using a coupled thermo-hydro-mechanical-chemical model

The low permeability of the methane hydrate-bearing sediment limits the methane gas extraction. To enhance methane hydrate extraction, hydraulic fracturing can be a promising approach to improve the hydrate reservoir permeability by creating a fracture network in the reservoir. In this study, a coupled thermo-hydro-mechanical-chemical mathematical model and its numerical implementation based on finite element technology are introduced to analyze the methane hydrate extraction through fractured methane hydrate-bearing sediment considering methane hydrates dissociation, gas-water two-phase flow, heat transfer, dynamic changes of the sediment permeability, and deformation of both sediment matrix and fractures as well as capturing the interplay between them. The coupled thermo-hydro-mechanical-chemical numerical model is verified by reproducing a methane hydrates dissociation laboratory test. Finally, we conduct a series of simulations for the methane gas depressurization extraction through the sediments with the DFNs assigned as diverse geometrical characteristics. The influence of hydraulic fracture network geometrical and hydraulic characteristics on methane hydrate extraction are discussed. The results can offer a reference for enhancing the methane hydrate extraction efficiency.

Effect of natural fracture on the heat mining of enhanced geothermal system based on the thermo-hydro-mechanical coupling model

ABSTRACT Hydraulic fracture, natural fracture and fracture networks play a vital role for heat mining of Enhanced Geothermal System (EGS). In this work, a two-dimensional thermo-hydro-mechanical (THM) EGS model was proposed to investigate the effect of natural fracture on heat extraction performance. Based on the local thermal non-equilibrium theory and discrete natural fracture network, the spatial and temporal distribution of temperature, pore pressure, effective stress, rock matrix porosity and permeability are investigated. The results show that the natural fracture has a significant negative effect on the production temperature. After the 40-year prediction, the production temperature of the base model is 45.5°C lower than that of the without natural fracture model. In the initial operation stage, the low-temperature region is mainly distributed in the region of artificial fracture and injection well, then it expands significantly to the natural fracture due to the stress redistribution. The key parameters of quantity and orientation of natural fracture can lead to a decline in production temperature and earlier thermal breakthroughs within the thermal reservoir.

Shale brittleness evaluation under high temperature and high confining stress based on energy evolution

Deep and ultra-deep shale gas reservoirs are crucial targets for future natural gas development. Extensive research has shown the shale brittleness is closely associated with the hydraulic fracture network quality. However, the shale brittleness may be altered in deep reservoirs due to the high temperature and high geo-stress. At present, commonly used shale brittleness evaluation methods cannot distinguish the differences in energy evolution caused by different mineral components in shale rock. In order to distinguish this difference, this study innovatively proposes a shale brittleness evaluation method that considers mineral composition based on energy evolution. A high-temperature triaxial compression experiment was carried out on the Longmaxi Formation shale, its fracture mode was analyzed, and the brittleness of the shale was evaluated based on this innovative method. The results indicate that with increasing confining stress, the shale fracture mode transitions from brittle tensile-shear multi-crack penetrating fracture to brittle-plastic shear single-crack fracture, and the brittleness of shale shows a decreasing trend. Under low confining stress, the shale brittleness decreases with temperature. Under high confining stress, temperature is the non-main controlling factor of shale brittleness. This study provides valuable guidance for the design and optimization of hydraulic fracturing in deep shale reservoirs.

Field Application of a Novel Multiresolution Multiwell Unconventional Reservoir Simulation: History Matching and Parameter Identification

Summary This study focuses on developing an efficient workflow by integrating a multiresolution simulation model and a multiobjective evolutionary algorithm (MOEA) for application to multiwell unconventional reservoirs. In this approach, hydraulic fractures are represented using a dual porosity, dual permeability system facilitated by an embedded discrete fracture model (EDFM). A novel fast-marching simulation method is used to cut down on computational expenses by an order of magnitude, greatly accelerating the history-matching process. A variety of integrated monitoring technologies were implemented to map out the hydraulic fracture network. Insights into hydraulic fracture locations were gleaned from warm-back analysis of distributed temperature sensing data, and these locations were then assimilated into the simulation model as embedded discrete fractures. For the simulation, a fast-marching-based multiresolution model was used to partition the reservoir into local and shared domains guided by the diffusive-time-of-flight (DTOF) principle. The local domain maintained the original 3D grids near the wells while transforming the remaining area into 1D grids to accelerate the simulation process. Before history matching, a thorough sensitivity analysis was conducted to pinpoint the most impactful parameters. Subsequently, the model was fine-tuned using production data through an MOEA. The most sensitive parameters in history matching were identified as fracture geometry and conductivity, fluid saturations, and rock compressibility in the stimulated reservoir volume (SRV) areas. After history matching, there was a noteworthy reduction in the uncertainty of these tuning parameters. The calibrated parameters are valuable to evaluate the effectiveness of the well completion design. Overall, this work emphasizes the innovative combination of techniques applied, the efficiency gains in the history-matching process, and the scalability of the approach to other oilfield applications.

Dynamic Pressure Analysis of Shale Gas Wells Considering Three-Dimensional Distribution and Properties of the Hydraulic Fracture Network

A major challenge in transient pressure analysis for shale gas wells is their complex transient flow behavior and fracturing parameters. While numerical simulations offer high accuracy, analytical models are attractive for transient pressure analysis due to their high computational efficiency and broad applicability. However, traditional analytical models are often oversimplified, making it difficult to capture the complex seepage system, and three-dimensional fracture characteristics are seldom considered. To address these limitations, this study presents a comprehensive hybrid model that characterizes the transient flow behavior and analyzes the pressure response of a fractured shale gas well with a three-dimensional discrete fracture. To achieve this, the hydraulic fracture is discretized into several panels, and the transient flow equation is numerically solved using the finite difference method. Based on the Langmuir adsorption isotherm and the pseudo-steady diffusion in matrix and Darcy flow in the network of micro-fractures, a reservoir model is established, and the Laplace transformation is adopted to solve the model analytically. The transient responses are obtained by dynamically coupling the flow in the reservoir and the discrete fracture. The precision of the proposed model is validated using the commercial numerical simulator, Eclipse. A series of transient pressure dynamic curves are drawn to make a precise observation of different flow regimes, and the effects of several parameters on transient pressure response are also examined. The results show that the shale gas well testing interpretation curves comprise nine flow stages. The pressure drop of shale gas reservoirs is lower than that of conventional gas reservoirs due to the replenishment of desorbed gas. The artificial fracture flow capacity, fracture length, and height are the main engineering factors affecting the pressure responses of shale gas wells. Maximizing the degree and scope of reconstruction can enhance the gas well production capacity during fracturing construction. The research results also indicate that our model is a reliable semi-analytical model for well test interpretations in real case studies.

Model Shows Computational Gains, Preserves Accuracy in Tight Rock EOR

_ This article, written by JPT Technology Editor Chris Carpenter, contains highlights of paper URTeC 3834855, “Advanced Dual-Porosity and Dual-Permeability Model for Tight Rock Integrated Primary Depletion and Enhanced Oil Recovery Simulation,” by Daegil Yang, SPE, Tyler Do, and Changdong Yang, SPE, Chevron, et al. The paper has not been peer reviewed. _ Discrete fracture or dual-porosity modeling has been widely adopted for tight rock simulation. To represent the complex fracture heterogeneity and simplify its characterization, the authors present in the complete paper a new simulation model to holistically optimize the tight rock reservoir completion, primary depletion, and enhanced oil recovery modeling and retain the injection mechanism and simplify the complex fracture characterization to enhance computational efficiency. The model enables users to characterize different properties in multiple regions. This tool will have broad applications in supporting asset-development decisions. Introduction Traditional dual-porosity, dual-permeability (DPDK) simulation incorporates microscale fracture characterization through the shape factor (matrix/fracture coupling), which offers a significant contact area for the gas/oil mixture to represent the field response reliably. However, its homogeneous fracture network cannot be used to optimize the completion design. The authors have developed what they refer to as the advanced DPDK (A-DPDK) simulation work flow to maintain heterogeneous fracture characterization and to simplify microfracture characterization. The work flow contains the following features: - A newly derived shape-factor algorithm that can represent the scale and heterogeneity of the fracture network - Local grid refinement (LGR) that accelerates the speed of computation (five-times-faster acceleration) and improves accuracy - A fracture-characterization feature that serves as a helpful diagnostic tool to assign a flexible cutoff value for both propped fractures (PF) and unpropped fractures to enhance simulation speed and robustness - A connected PF feature that can label the effective fracture region to simplify completion design work flow and prioritize reservoir simulation The authors validated the new features implemented in the work flow systematically and demonstrated the benefits of the proposed work flow. The authors also applied the work flow in one gas-injection pilot and successfully history-matched the primary depletion and the gas-injection response. Modeling Approach The A-DPDK model design is aimed at modeling a complex fracture network using a structured grid system while maintaining the background natural fracture. In the modeling process, the main fracture network is screened and identified and its transmissibility calculated. Then, the microfracture is calculated to matrix mass transfer. This can be numerically modeled using a DPDK framework. The geometry of the hydraulic fracture network is preserved explicitly as grid cells in the fracture domain. Also, in the matrix domain, the enhanced permeability region (EPR), which represents the microfracture, is included. The mass transfer between matrix cells and the EPR is modeled as a grid-cell connection (transmissibility), and the mass transfer between the hydraulic fracture cells and the microfractures is modeled as a grid-cell connection as well. A new way of calculating the coupling term that models the mass transfer between the hydraulic fracture and the EPR and the matrix is introduced.

A numerical study of complex fractures propagation path in naturally fractured gas shale: Reorientation, deflection and convergence

Shale gas reservoir is characterized as the naturally fractured formation, which is benefit for complex hydraulic fracture network generation. The propagation path is crucial to determine the reservoir stimulation result, and evaluate the quality of complex hydraulic fracture network. To investigate the propagation path of complex hydraulic fracture in gas shale, the hybrid combined discrete fracture network-finite element method (DFN-FEM) model is involved. Combining a linear hardening and a non-linear softening process, a non-linear traction-separation law is proposed to describe the fracture damage evolution. Validated with open-reported laboratory experiments, the hydraulic-natural fracture (HF–NF) interaction criterion is depicted under four geomechanical and injection factors. Based on a field-scale model of Sichuan basin shale gas reservoir, the complex hydraulic fractures propagation process is simulated with rock deformation and stress shadow. Several propagation path indexes are defined with fracture morphology to discuss the effect of representative geomechanical, injection and multi-clustering parameters. The results show that the tension fractures are predominant to hydraulic fractures propagation path and shear the natural fractures, due to the combined influence of natural fractures network and interference between clusters. The horizontal stress contrast and the injection fluid viscosity is governing to the propagation path and morphology, rather than the friction coefficient and injection rate. Multi-clustered hydraulic fractures converge into a dominant fracture due to the attraction with a short clusters-spacing, while the exclusion is observed along the propagation path with the increase of clusters number.

Production induced fracture closure of deep shale gas well under thermo-hydro-mechanical conditions

Deep shale gas reservoirs have geological characteristics of high temperature, high pressure, high stress, and inferior ability to pass through fluids. The multi-stage fractured horizontal well is the key to exploiting the deep shale gas reservoir. However, during the production process, the effectiveness of the hydraulic fracture network decreases with the closure of fractures, which accelerates the decline of shale gas production. In this paper, we addressed the problems of unclear fracture closure mechanisms and low accuracy of shale gas production prediction during deep shale gas production. Then we established the fluid–solid–heat coupled model coupling the deformation and fluid flow among the fracture surface, proppant and the shale matrix. When the fluid–solid–heat coupled model was applied to the fracture network, it was well solved by our numerical method named discontinuous discrete fracture method. Compared with the conventional discrete fracture method, the discontinuous discrete fracture method can describe the three-dimensional morphology of the fracture while considering the effect of the change of fracture surface permeation coefficient on the coupled fracture–matrix flow and describing the displacement discontinuity across the fracture. Numerical simulations revealed that the degree of fracture closure increases as the production time proceeds, and the degree of closure of the secondary fractures is higher than that of the primary fractures. Shale creep and proppant embedment both increase the degree of fracture closure. The reduction in fracture surface permeability due to proppant embedment reduces the rate of fluid transfer between matrix and fracture, which has often been overlooked in the past. However, it significantly impacts shale gas production, with calculations showing a 24.7% cumulative three-year yield reduction. This study is helpful to understand the mechanism of hydraulic fracture closure. Therefore, it provides the theoretical guidance for maintaining the long-term effectiveness of hydraulic fractures.

Fracture Initiation and Propagation in the Hot Dry Rock Subject to the Cyclic Injection Hydraulic Fracturing Treatment

Cyclic injection hydraulic fracturing is a promising way for the geothermal energy exploitation by reactivating the fractures in geothermal reservoir. However, fracture initiation and growth induced by cyclic injection schemes have been inadequately studied for hot dry rock (HDR), and the cyclic injection fracturing optimized often by experience. For this reason, the initiation and propagation of hydraulic fractures in the HDR under different cyclic injection methods were determined by experiment research for hydraulic fracturing. The results show that the cyclic frequency and injection rate play different roles in the stimulation of HDR. The cyclic injection with low frequency-low pressure can create more branched fractures, forming a short but complex hydraulic fracture network. However, when high flow-high frequency injection method is subjected, the branch fractures formed are significantly reduced, but each branch fracture can be fully expanded. To fully exploit the advantages of different injection methods, a numerical model that contains a fracture network was established with PFC software, and an alternating cyclic injection scheme with synergistic control of the cyclic frequency and injection rate was proposed. The comparison results indicated that the alternating cyclic injection method can effectively improve the fracturing effect in the HDR. The stimulation area of the alternating cyclic injection method is about 2.3 times and 2.7 times that of the low flow-low frequency and high flow-high frequency injection methods, respectively. The method presented here can be adopted to optimize the fracture growth regime and provide a scientific basis for EGS hydraulic fracturing design.

Model and Analysis of Pump-Stopping Pressure Drop with Consideration of Hydraulic Fracture Network in Tight Oil Reservoirs

The existing pump-stopping pressure drop models for the hydraulic fracturing operation of tight oil reservoirs only consider the main hydraulic fracture and the single-phase flow of fracturing fluid. In this paper, a new pump-stopping pressure drop model for fracturing operation based on coupling calculation of the secondary fracture and oil-water two-phase flow is proposed. The physical model includes the horizontal wellbore, the fracture network and the tight oil reservoir. Through the numerical simulation and calculation, the wellbore afterflow performance, the crossflow performance between the main hydraulic fracture and the secondary fracture, the fracturing fluid leakoff and the oil-water replacement after termination of pumping are obtained. The pressure drop characteristic curve is drawn out by the bottom-hole flow pressure calculated through the numerical simulation, and a series of analyses are carried out on the calculated pressure drop curve, which is helpful to diagnose the -oil-water two-phase flow state and the fracture closure performance under the control of the fracture network after hydraulic fracturing pumping. Finally, taking a multi-stage fractured horizontal well in a tight oil reservoir in the Junggar basin, China as an example, the pump-stopping pressure drop data of each stage after hydraulic fracturing are analyzed. Through the history fitting of the pressure drop characteristic curve, the key parameters such as fracture network parameters, which include the half-length of main hydraulic fracture, the conductivity of main hydraulic fracture and the density of secondary fracture, the fracture closure pressure are obtained by inversion, thus, the hydraulic fracturing effect of fractured horizontal well in tight oil reservoirs is further quantified.

An efficient 3D cell-based discrete fracture-matrix flow model for digitally captured fracture networks.

Complex hydraulic fracture networks are critical for enhancing permeability in unconventional reservoirs and mining industries. However, accurately simulating the fluid flow in realistic fracture networks (compared to the statistical fracture networks) is still challenging due to the fracture complexity and computational burden. This work proposes a simple yet efficient numerical framework for the flow simulation in fractured porous media obtained by 3D high-resolution images, aiming at both computational accuracy and efficiency. The fractured rock with complex fracture geometries is numerically constructed with a cell-based discrete fracture-matrix model (DFM) having implicit fracture apertures. The flow in the complex fractured porous media (including matrix flow, fracture flow, as well as exchange flow) is simulated with a pipe-based cell-centered finite volume method. The performance of this model is validated against analytical/numerical solutions. Then a lab-scale true triaxial hydraulically fractured shale sample is reconstructed, and the fluid flow in this realistic fracture network is simulated. Results suggest that the proposed method achieves a good balance between computational efficiency and accuracy. The complex fracture networks control the fluid flow process, and the opened natural fractures behave as primary fluid pathways. Heterogeneous and anisotropic features of fluid flow are well captured with the present model.

Unstable propagation of hydraulic fractures under varying well spacing and initiation sequence of multiple horizontal wells considering shear stress shadows

PurposeMulti-well hydrofracturing is an important technology to create new fractures and expand existing fractures to increase reservoir permeability. The propagation morphology of the fracture network is affected by the disturbance between the fractures initiation sequences and spacings between adjacent wells. However, it remains unclear how well spacing and initiation sequences lead to fracture propagation, deflection and connection.Design/methodology/approachIn this study, the thermal-hydro-mechanical coupling effect in the hydrofracturing process was considered, to establish a finite element-discrete element model of multistage hydrofracturing in a horizontal well. Using typical cases, the unstable propagation of hydraulic fractures in multiple horizontal wells was investigated under varying well spacing and initiation sequences. Combined with the shear stress shadow caused by in situ stress disturbed by fracture tip propagation, the quantitative indexes of fracture propagation such as length, volume, displacement vector, deflection and unstable propagation behavior of the hydrofracturing fracture network were analyzed.FindingsThe results show that the shear stress disturbance caused by multiple hydraulic fractures is a significant factor in multi-well hydrofracturing. Reducing the spacing between multiple wells increases the stress shadow area and aggravates the mutual disturbance and deflection between the fractures. The quantitative analysis results show that a decrease of well spacing reduces the total length of hydraulic fractures but increases the total volume of the fracture; compared with sequential and simultaneous fracturing, alternate fracturing can effectively reduce stress shadow area, alleviate fracture disturbance and generate larger fracture propagation length and volume.Originality/valueThe numerical models and results of the unstable propagation and stress evolution of the hydraulic fracture network under thermal-hydro-mechanical coupling obtained in this study can provide useful guidance for the evaluation and design of rock mass fracture networks in deep unconventional oil and gas reservoirs.