Search

Natural fractures may not be developed in hot dry rock reservoirs, and tectonism or fracturing can induce the shear deformation of fractures and result in large-scale rough-walled fractures. Previous studies usually ignore the effect of rough walls on production performance. This study constructs a 3D multiwell enhanced geothermal system (EGS) model with two rough-walled fractures after shear deformation, systematically investigates the response characteristics of each physical field under the thermo–hydro–mechanical coupling, and determines the effect of the distribution of two rough-walled fractures and the relationship between the well layout scheme and the shear deformation direction on the production performance. The obtained results indicate that the distribution and evolution of the temperature and seepage fields are controlled by the fractures. The stress has a clear impact on the fracture permeability, up to approximately four times the initial permeability. The water preferentially extracts the heat from the reservoir rock between fractures. When the intersecting angles of fractures are 10° and 30°, there exists a sufficient heat exchange domain and a shorter average distance between fractures, which can enhance the production efficiency. The well layout perpendicular to the shear deformation direction is conducive to delaying the lifespan of the EGS and extracting more heat over a longer exploitation time. These key results can provide reasonable suggestions for the optimal development strategy in EGSs. [ABSTRACT FROM AUTHOR]

The mid-deep geothermal energy is extracted from the fracture network created by hydraulic fracturing. Generally, the hydraulic fracture created in the deep thermal reservoir shows the planar shape but not the network shape, and this significantly impacts the geothermal extraction efficiency. In this study, a geothermal extraction technology by radial well assisted hydraulic fracturing, which connects hydraulic fractures to each other, was proposed. An operation technique of radial well-assisted hydraulic fracturing in geothermal reservoirs without well-developed NFs was established. The model considering evolution of fracture permeability was resolved by thermo-hydromechanical (THM) coupling numerical simulation. The results show that the radial wells should be placed vertically within a large range. The thermal production increases by 20%, about 0.25 MW as the radial well length increases from 50 m to 200 m. As many fractures as possible are needed in thermal exploitation, and the fracture number should be at least 7. The simulation demonstrates that the radial well-assisted hydraulic fracturing heat recovery method is feasible in simulation. In this thesis, a more accurate numerical simulation method of geothermal production capacity is proposed, which has engineering guiding significance for the actual development and prediction of geothermal energy. [ABSTRACT FROM AUTHOR]

Natural gas hydrate (NGH) is a prospective clean energy. Conventional development method of the NGH reservoir by pressure drop in vertical well faces the problems of limited stimulated reservoir volume, low production efficiency, and low conductivity. At present, pressure drop production in various well types and stimulation treatments such as hydraulic fracturing in fracable NGH reservoirs have become the focus. The effects of pressure drop production in vertical well, horizontal well, radial wells, and hydraulic fracturing in vertical and horizontal wells in NGH reservoirs were simulated with HydrateResSim. We simulated the productivity of NGH reservoir for 1000 days with several types of development mode. The productivity is highest in horizontal well staged fracturing (5 stages), followed by pressure drop production in horizontal well, fractured vertical well (single fracture), and radial well (16 holes). The worst effect is obtained through production by pressure drop in single vertical well. The optimal effects were obtained with the fracture spacing of 58 m, the fracture half-length of 100 m and the conductivity of 30 μm2 cm. The horizontal well multi-cluster fracturing significantly enlarges the stimulated area and provides channels for gas migration and enhances productivity. [ABSTRACT FROM AUTHOR]

A proppant transport simulator coupled with multi-planar 3D (multi-PL3D) fracture propagation has been developed to examine the proppant distribution among multiple hydraulic fractures during multi-cluster fracturing in a horizontal well. The multi-PL3D fracture model considers wellbore friction, multi-fracture stress interaction, fluid leak-off, and multi-scale propagation regimes. The proppant transport is described by the two-phase (slurry/proppant) flow equations that consider proppant settling, jamming and flow regime transition. A high-resolution weighted essentially non-oscillatory (WENO) finite difference (FD) scheme is adopted to solve the nonlinear proppant transport equations. An efficient time-stepping scheme is developed to solve the solid/fluid coupling equations and moving boundaries for the multi-PL3D model. The proppant transport model and multi-PL3D model are both validated against previously published results. Using the model, we examine the proppant distributions under different injection schedules, proppant sizes, proppant density, and fluid viscosity. Results show that proppant distribution among multiple fractures is different as the flow rate and fracture width distribution vary due to multi-fracture stress interaction. The proppant in the middle cluster settles remarkably as the flow rate is lowest among the multiple clusters. The proppant is usually jammed at the pinch point, where the fracture width reduces sharply. Proppant adding schedule has a significant effect on the proppant distribution. A constant-concentration results in a proppant stack at the fracture front. In contrast, an increasing concentration favors the prop of the near-wellbore fracture. The proppant distribution area ratio (defined as the proppant distribution area divided by the fracture area) is only 20% for 20/40 mesh proppant, while the ratio is 45% for 100 mesh proppant. Slick water can increase the fracture area but not favor promoting the proppant distribution area ratio. The results can be helpful for proppant design for multi-cluster fracturing in a horizontal well. [ABSTRACT FROM AUTHOR]

Shale reservoirs are usually buried deep, and contain lots of complex natural fractures (NFs) under the influence of tectonic stress and faults. Previous studies usually focus on small-scale reservoir models and there is also a lack of research that completely focuses on NF characteristics. Using the continuous-discontinuous element method (CDEM), a Hydro-Mechanical (HM) coupling model was constructed to explore the propagation mechanism of hydraulic fracture (HF) and the effects of NF characteristics on HF propagation. The results show that as the HF approaches the NFs, the NFs preferentially tend to undergo shear failure, resulting in a decelerated HF propagation rate and an increased injection pressure. After water enters the NFs, the HF rapidly propagates, so the fracture area swiftly increases and the HF has a more tortuous fracture morphology. The NF orientation (α) and the internal friction angle (φ) exhibit a pronounced and distinct effect on HF propagation. Specifically, when α is 30°∼60°, the HF propagation is more tortuous and complex, so it is easier to form a larger transverse reconstruction range. When the NF density (η) is no less than 0.12, the HF tends to open the NFs and form a more tortuous morphology, but a higher η does not necessarily imply a higher fracture complexity. A longer NF (l > 12 m) benefits to improve the longitudinal reservoir reconstruction scope, but has little impact on the overall reconstruction effect. In the presence of stress shadow effect, it is improbable for the HF to completely open the intersecting NFs; instead, it typically opens only a segment or one of them. Consequently, it may be challenging to establish a fracture network with the simultaneous development of primary and branching fractures. Additionally, some water leakage into the inactive NFs gives rise to a decline in both the fracture length and fracture area. The findings can offer theoretical insights for refining the fracturing design in shale reservoirs with NFs. [ABSTRACT FROM AUTHOR]

Hydraulic fracturing, a production enhancement technique, is widely used in the development of unconventional oil and gas reservoirs. The formation of a complex network of hydraulic fractures that connect with natural fractures is crucial for hydraulic fracturing in unconventional reservoirs. However, the current understanding of vertical fracture propagation behavior under the influence of variations in the mechanical properties of interbedded rock is insufficient to meet the requirements for simulations of unconventional reservoirs under complex geological conditions. In this study, a three-dimensional discrete grid method was employed to establish a model of three-dimensional fracture propagation. Numerical simulations were conducted to investigate the vertical growth of fractures while considering the influence of bedding planes. The effects of formation factors (stress, Young's modulus) and bedding planes (cohesion, density) on the height growth of hydraulic fractures were explored. The results indicated that the interbedded stress contrast, Young's modulus contrast, and bedding planes collectively controlled the height growth of hydraulic fractures. The height of hydraulic fractures decreased with increasing minimum horizontal principal stress of adjacent layers. When the minimum horizontal principal stress of adjacent layers exceeded the vertical stress, hydraulic fractures gradually deflected into the horizontal plane. Adjacent layers with large values of Young's modulus promoted the height growth of hydraulic fractures, while adjacent layers with small values of Young's modulus inhibited the height growth of hydraulic fractures. The presence of bedding planes further suppressed height growth, and the degree of suppression was related to the cohesion and density of the bedding planes. Weaker cohesion and higher density resulted in greater suppression. The results of this study provide a reference for the design and optimization of hydraulic fracturing treatments in unconventional oil and gas reservoirs. [ABSTRACT FROM AUTHOR]

As a green, low-carbon and renewable energy, the geothermal energy has attracted great attentions. Abundant geothermal resources are hosted in several depleted oilfields. Geothermal development in abandoned oilfields has its advantages that operators are familiar with the geographic environment, geological conditions, and drilling and development techniques in the oilfield. In this study, a thermo-hydro- mechanical multi-field coupling mathematical model was established with data from Liubei Buried Hill reservoir in North China Oilfield. Then, the geothermal development in abandoned high-temperature oilfields was simulated, and the sensitivity of the geothermal development effect to well patterns, fracturing conditions, and working media was analyzed. The results show that the well pattern of one-injector and four-producer leads to the largest heat-exchange area. The hydraulic fracture model causes earlier thermal breakthrough but produces more heat than the non-fractured model. CO 2 has the better working fluid performance and leads to the better geothermal development effect than water. This study provides a theoretical basis for the high-efficiency development of geothermal energy in abandoned oilfields, and a guidance for research on geothermal development model and fracturing design. [Display omitted] [ABSTRACT FROM AUTHOR]

In unconventional reservoir hydraulic fracturing, using a pump schedule with proppant concentration constant (PCC) instead of the proppant concentration stepwise increasing (PCSI) is considered an innovative approach to reduce operation costs and risks. However, the large-scale application of PCC lacks sufficient theoretical support. This study investigates the impact of two pump schedules on proppant transport behavior through visualized proppant transport experiments and the model that couples fracture propagation with proppant transport. The experiments demonstrate that PCC facilitates proppant transport and reduces flow fluctuations in complex fractures, while PCSI enhances the filling of near-well fractures. The proppant concentration (the ratio of proppant volume to fracturing fluid volume) reaches 20 % or when a combination of 40/70-mesh and 20/40-mesh quartz sands is used for injection, the pump schedule influence becomes negligible. An increase in fracture complexity, proppant size, and proppant concentration all lead to a reduced influence of the pump schedule on the dune shape. Numerical simulation reveal that the fracture area and average fracture conductivity is similar for different pump schedules, but PCC yields a slightly larger propped fracture area (ranging from 1 % to 4 %) and results in a more uniform conductivity. PCC could achieve similar fracture filling when the proppant concentration is more than 15 % or there is no significant demand for near-well conductivity. Overall, both two pump schedules have a relatively insignificant impact on the transport process of proppants and the final sand dune shape. By comparing pressure curves during hydraulic fracturing, it was evident that PCC resulted in fewer fluctuations. Analyzing production data from 93 wells with PCSI and 59 wells with PCC, it was observed that PCC enhances production and reduces proppant flowback post-fracturing. This research provides a theoretical foundation for the application of PCC pump schedule in hydraulic fracturing. • Compared the effects of different pump schedule on proppant transport. • Proppant concentration constant and stepwise injection have little difference. • Demonstrated the feasibility of proppant concentration constant injection. • Field tests show proppant concentration constant injection performs equally or better. [ABSTRACT FROM AUTHOR]

The deep geothermal energy produced from Enhanced Geothermal System (EGS) has a great development prospect because of enormous potential and environmental friendliness. EGS process involves a complex thermal-hydraulic process, and fractures in EGS are main channels for fluid flow and heat transfer, the understanding of which is crucial to the sustainable utilization of geothermal reservoirs. In this paper, a 3D thermal-hydraulic coupled numerical model is proposed to describe the interaction of fluid flow and heat transfer. Besides, the EGS with multiple fracturing horizontal wells (MFHW) is adopted to evaluate the effect of multiple hydraulic fractures on geothermal energy extraction performance. The MFHW with multiple stimulated fractures could increase fluid flow path and heat exchange area significantly, thereby enhance the heat recovery ability. Firstly, we analyzed the evolution of temperature and flow fields in EGS and compared the MFHW EGS with conventional vertical EGS. Secondly, the effects of fracturing parameters, including the fracture number, fracture length, and fracture conductivity, on heat extraction performance were investigated. Finally, the cost for drilling and hydraulic fracturing in MFHW EGS was calculated. The results indicate that MFHW EGS has a higher cumulative thermal production and a better heat extraction performance than that of conventional vertical EGS. For the optimization of hydraulic fracture parameters, the cumulative thermal production firstly increases and then decreases as the fracture number increases, the cumulative thermal production curve exists an inflection point of fracture number. Longer fracture length and higher fracture conductivity could enhance the cumulative thermal production, but the output growth slows down gradually. Considering economic cost, the best fracture parameters for MFHW EGS in this paper are the fracture number of 7, the fracture length of 300 m, and the fracture conductivity of 350 μm2•cm, respectively. The research provides a better study for multiple fracturing horizontal wells (MFHW) EGS and helps to optimize fracture parameters and geothermal reservoir management, which is conductive to improve the geothermal energy efficiency. • An EGS with multiple fracturing horizontal well (MFHW) was put forward. • A 3D MFHW thermal-hydraulic coupled model was established. • The MFHW EGS was compared with conventional vertical EGS. • The optimal fracture parameters for MFHW EGS were concluded. [ABSTRACT FROM AUTHOR]

Methane deflagration fracturing is a new reservoir stimulation method that serves the efficient development of shale gas reservoirs. However, the propagation law of deflagration fractures is still unclear. In this paper, a numerical model considering the effect of stress wave impact and gas drive of deflagration fracturing was established based on the continuum–discontinuum element method (CDEM). The correctness of the numerical model was verified by comparing it with a laboratory experiment, the steady and unsteady analytical solutions of gas flow, and the approximate solution of fracture propagation. Then, numerical simulations of methane deflagration fracturing in vertical wells and horizontal wells under different factors were carried out to analyze the fracture mechanism. The results indicate that deflagration fracturing in vertical wells can break through the stress concentration around the borehole; the initial radial fractures are formed under the action of stress wave impact and then propagate substantially under the driving action of high-pressure gas. The in-situ stress difference affects the deflagration fracture propagation and makes the half-fracture length in the direction of maximum principal stress larger than that in the direction of minimum principal stress. The more significant the stress difference is, the more noticeable this deviation will be. When the deflagration peak pressure is high, the reservoir burst degree is large, which is conducive to enlarging the stimulation range of deflagration fracturing. Staged deflagration fracturing in horizontal wells can form 5–8 obvious fractures perpendicular to the horizontal borehole in each explosion section. A large cluster spacing and explosion section length are conducive to expanding the stimulation scope. Moreover, the propagation of deflagration fractures will be induced by the natural fractures, and the natural fracture with a considerable length or a slight angle between the dip angle and the propagation direction of deflagration fractures is more likely to be activated. • The deflagration fracturing for shale gas reservoir stimulation is proposed. • The numerical model considers the effects of stress wave and gas drive. • The propagation law of deflagration fracture is revealed by numerical simulation. • Deflagration fracturing can break the stress concentration to form complex fracture. [ABSTRACT FROM AUTHOR]

Abstract The SRV fracturing is the core of commercial exploitation of shale gas, and the design philosophy of the SRV fracturing is gradually introduced into the stimulation of conventional low permeability reservoir and unconventional reservoir. However, in reservoirs with poorly developed natural fractures or high horizontal principal stress difference, it is difficult to achieve the multi-fracture or fracture network by conventional hydraulic fracturing, and effectively stimulating reservoir needs other assisting methods. This paper aims to present a method of guiding the multi-directional hydraulic fractures by multiple radial slim holes drilled in spatial positions, and creating multi-fractures and even complex fracture network. In this paper, a 3D hydraulic fracturing extended finite element numerical model in the vertical well assisted by planar multi-radial slim holes was established, and the effect of radial slim hole number, diameter, length, vertical density, azimuth and phase angle, Poisson's ratio, horizontal principal stress difference and rock Young's modulus, fracturing fluid viscosity and injection rate, reservoir permeability, etc. on guidance of radial slim holes were investigated. Then, we conducted the true tri-axial hydraulic fracturing test with the artificial core, and the feasibility of the technology was verified. The results show that decreasing azimuth and phase angle and increasing hole diameter, length and vertical density strengthen the guiding effect of planar multi-radial slim holes, and increasing Poisson's ratio, reservoir permeability, and fracturing fluid injection rate and viscosity enhance the guidance effect of radial slim holes. Increasing horizontal principal stress difference and reservoir Young's modulus weaken the radial slim holes guiding effect. The physical modeling experiment proves that controlling generation of complex multi-fractures is feasible by rationally arranging the spatial positions of radial slim holes. The spatial arrangement of four planar radial slim holes with phase angle of 90°, azimuth of 45° and vertical density not less than 2 holes/m in the vertical well is a scientific scheme. With the reasonable treatment scheme, creating 6 main fractures is possible. The research results provide theoretical support for multi-fractured stimulation assisted by planar multi-radial slim holes, and help the design of well completion parameters and fracturing treatment parameters of the technology. Highlights • The new multi-fractured stimulation technique for hydraulic fracturing was put forward. • Stimulation experiments of radial slim holes-assisted hydraulic fracturing were carried out. • The effects of multiple factors on the guiding strength of radial slim holes were studied. • The best technical solution for radial slim holes-assisted hydraulic fracturing was proposed. [ABSTRACT FROM AUTHOR]

Abstract Until now, hydraulic fracturing has played prominent role in increasing production of hydrogen energy such as natural gas, geothermal energy, natural gas hydrate. Hydrogen energy stimulation is realized by many methods. Such as multiple-clusters in staged fracturing, and fracture propagation guided by combined radial boreholes in different azimuths. In order to extend lateral stimulation area, the zipper fracturing is investigated continuously. However, the report on the zipper fracturing is limited to the post-frac productivity. In this study, model of simultaneous and zipper fracturing was established in ABAQUS to investigate the effect of six factors on the fracture propagation. The results showed that in simultaneous fracturing, the fracture network is formed spontaneously. With same brittleness of reservoir, fracture in zipper fracturing always propagates slightly longer than that in simultaneous fracturing. The results provide theoretically support for both fracturing modes, which helps design of well completion and fracturing operation parameters. Highlights • The effects of multiple factors on fracture propagation of the new crafts on hydraulic fracturing was studied. • Simulation experiments of zipper fracturing were in-depth carried out. • Simulation of fracturing cracking stress was studied. • Difference of fracture propagation between simultaneous fracturing and zipper fracturing was discussed. [ABSTRACT FROM AUTHOR]

Under the influence of tectonic stress and faults, there are many natural fractures (NFs) in hot dry rock (HDR) reservoirs. However, the effect of NFs on the propagation of hydraulic fractures (HFs) has often not been considered in previous studies. A coupled thermo-hydro-mechanical-damage (THMD) model by the finite element method is constructed to investigate hydraulic fracturing and production performance in a fractured geothermal reservoir. The results show that thermal stress and injection pressure jointly promote rock initiation and propagation. HF can propagate tortuously and deviate from the maximum horizontal stress direction under the influence of NFs. HF and open NFs constitute the main channel of heat transfer, which determines the heat extraction performance. The fluid viscosity (μ) and injection flow rate (q in) are the main factors affecting the fracturing and production performance. A higher μ and q in can significantly increase HF length and activate the NFs along the propagation path, so the fracture area and production temperature have a distant ascension. The horizontal stress difference (Δ σ) and NF number (n) are secondary factors affecting the fracturing and production performances. A higher or lower Δ σ and n are not conducive to forming a long HF and enhancing the production temperature. A larger fracture area or higher fracture complexity does not necessarily lead to better production performance. For the development mode of HDR using two injection wells and one production well, it is beneficial to enhance the production temperature to give priority to ensuring the uniform propagation of HF. The results can provide a theoretical basis for the optimal design of an enhanced geothermal system in a fractured geothermal reservoir. • A coupled THMD model for the fracturing and production is constructed and verified • HF propagation under the influence of NFs and the production performance are studied • The main controlling factors affecting the fracturing and production are determined • The uniform propagation of HF is beneficial to improve the production performance [ABSTRACT FROM AUTHOR]

The conventional hydraulic fracturing fails in the target oil/gas zone (remaining oil, closed reservoir, etc.) which is not located in the azimuth of maximum horizontal stress of available wellbores. The technology of directional propagation of hydraulic fracture guided by vertical multi-radial slim holes is innovatively developed. In order to verify this technology, lots of true triaxial hydraulic fracturing simulation experiments were carried out with artificial cores, aiming at investigating the influence of in-situ stress, fracturing fluid displacement, and azimuth, diameter, number, and spacing of radial slim holes on fracture propagation. The results show that directional propagation of hydraulic fractures guided by vertical multi-radial slim holes is feasible. Regardless of varied radial slim hole and in-situ stress parameters, hydraulic fracture always initiates in the heel of radial slim hole. For hydraulic fracturing assisted by single radial slim hole, smaller horizontal stress difference (3 MPa, σ H /σ h  = 1.25) and smaller radial slim hole azimuth (15°) may guide propagation of hydraulic fractures along radial slim hole, and larger horizontal stress difference (6 MPa, σ H /σ h  = 1.67) can sharply reduce guidance of radial slim hole. Affected by mutual interference from guidance of radial slim holes and controlling of horizontal in-situ stress, fractures tend to distort along fracture length and fracture height. For hydraulic fracturing assisted by vertical multi-radial slim holes, horizontal stress difference is one of key factors influencing the directional propagation effect of hydraulic fracture. Larger horizontal in-situ stress difference (>6 MPa, σ H /σ h  > 1.67) hardly results in directional propagation of hydraulic fracture along radial slim hole row, and creates bigger diversion angle of fracture, and larger azimuth (>45°) of radial slim hole row reduces the guidance strength, resulting small fracture diversion radius. In conditions of larger angle between target zone and maximum horizontal stress, and requirement of large azimuth of radial slim hole row, the effective hydraulic fracture propagation along radial slim holes and ideal fracture height will be achieved through human intervention, e.g. increment of the number and diameter of radial slim holes, reduction of radial slim holes spacing, increment of fracturing fluid displacement, etc. The simulation results provide methods for designing directional propagation of hydraulic fractures guided by vertical multi-radial slim holes. [ABSTRACT FROM AUTHOR]

The clear mechanism of hydraulic fracture propagation in glutenite reservoirs with high heterogeneity is still not obtained, thus it is difficult to carry out the design of fracturing plan effectively. Based on the characteristics of the glutenite reservoirs, a coupled flow-stress-damage (FSD) model of hydraulic fracture propagation with gravels is established. This model is experimentally verified and the research on the influence of rock physical parameters and gravel property on the hydraulic fracture propagation is conducted. It is shown that as the gravel tensile strength increases, the hydraulic fracture is prone to propagate around the gravel, where the fracture deflection always occurs; as the gravel Young's modulus increases, there is high probability that hydraulic fracture propagates around the gravel, with more obvious fracture deflection; the matrix permeability influences fracture propagating morphology when encountering gravel and total fracture length; the horizontal geostress difference seriously impacts the fracture deflection; as the fracturing fluids injection displacement increases, the fracture is prone to deflect when encountering gravel; the low viscosity fracturing fluids result in the shorter fracture; the larger gravel increases the possibility of fracture deflection; in case of smaller gravel sizes, the increasing gravel content has a big influence on fracture deflection, and the increasing content of large gravel complicates the fracture morphology, resulting in the fine branched fractures; for the well rounded gravel, the fracture propagation around the gravel is prone to occur, and the fracture is not prone to deflect. Compared with the conventional sandstone reservoir, the glutenite reservoirs have higher breakdown and extension pressures, which fluctuate due to the gravel; the larger gravel size results in higher extension pressure. In this paper, a simulation method of hydraulic fracture propagation in the glutenite reservoirs is introduced, and the result provides the theoretical support for prediction of fracture propagation morphology and plan design of hydraulic fracturing in the glutenite reservoirs. [ABSTRACT FROM AUTHOR]

The conventional hydraulic fracturing fails in the target oil development zone (remaining oil or gas, closed reservoir, etc.) which is not located in the azimuth of maximum horizontal in situ stress of available wellbores. The technology of directional propagation of hydraulic fracture guided by vertical multi-radial boreholes is innovatively developed. The effects of in situ stress, wellbore internal pressure and fracturing fluid percolation effect on geostress field distribution are taken into account, a mechanical model of two radial boreholes (basic research unit) is established, and the distribution and change rule of the maximum principal stress on the various parameters have been studied. The results show that as the radial borehole azimuth increases, the preferential rock tensile fracturing in the axial plane of radial boreholes becomes increasingly difficult. When the radial borehole azimuth increases to a certain extent, the maximum principal stress no longer appears in the azimuth of the radial boreholes, but will go to other orientations outside the axial plane of radial boreholes and the maximum horizontal stress orientation. Therefore, by reducing the ratio between the distance of the radial boreholes and increasing the diameter of the radial boreholes can enhance the guiding strength. In the axial plane of the radical boreholes, particularly in the radial hole wall, position closer to the radial boreholes is more prone to rock tensile destruction. Even in the case of large radial borehole azimuth, rock still preferentially ruptures in this position. The more the position is perpendicularly far from the axis of the wellbore, the lesser it will be affected by wellbore, and the lesser the tensile stress of each point. Meanwhile, at a certain depth, due to the decrease in the impact of the wellbore and the impact of the two radial boreholes increases accordingly, at the further position from the wellbore axis, the tensile fracture is the most prone to occur and it will be closer to the axial plane of the two radial boreholes. The study provides theoretical support for the technology of directional propagation of hydraulic fracture promoted by radial borehole, which is helpful for planning well-completion parameters in technology of hydraulic fracturing promoted by radial borehole. [ABSTRACT FROM AUTHOR]

In order to study the effects of different parameters in hydraulic fracturing on formation temperature after low-temperature fracturing fluid is injected into high-temperature reservoir, in this paper, considering (1) osmotic pressure and capillary pressure (2) microthermal effect (heat conduction, heat convection, heat expansion, viscous dissipation) (3) the coupling processes of temperature (T), hydraulic, (H) mechanical (M) and chemical (C) in hydraulic fracturing process, and the coupled temperature prediction model of oil-water two-phase THMC based on discrete fractures is established. The effects of fracturing fluid temperature, reservoir temperature, dimensionless conductivity, Young's modulus, injection rate, cluster spacing, and proportion of branch fracture area on reservoir and fracture temperature during the pumping and shut-in stages of shale oil reservoir were studied. The results show (1) For formations with different reservoir temperatures, the temperature drop at the fracture was approximately 97% from the start of the operation to the end of the 30-day shut-in. The saturation change caused by fracturing fluid filtration is much larger than the temperature change caused by heat transfer, and it is found that the temperature and saturation change mainly occur in the early shut-in period (2) Imbibition accelerates the filtration loss of fracturing fluid. During pumping, the bottom hole temperature with imbibition considered is lower than that without imbibition considered. The temperature curve with imbibition considered begins to deviate about 10 min after pumping, but the final temperature tends to be the same, (3) When the fracturing fluid is pumped at high pressure, the fracture and matrix will be deformed, which will improve the porosity and permeability of the near well zone, thus speeding up the heat transfer rate. The high Young's modulus inhibits the increase of permeability. Compared with the young's modulus of 40 GPa and 20 GPa, the increase of permeability decreases by 9% (4) When the well was shut in for 30 days, the cluster spacing was greater than 20 m, and the temperature change around the fracture did not affect the clusters on both sides. The spacing between the clusters was 5 m, and the perforation clusters on both sides had a significant effect on the temperature. The temperature of the reservoir between the two adjacent clusters dropped from 393 K to 385 K. When there is a branch fracture, the small cluster spacing is equivalent to increasing the complexity of the fracture and accelerating the initial heat transfer rate. (5) High injection temperature, high dimensionless conductivity, low Young's modulus and low pumping rate are conducive to the heat transfer of fracturing fluid in the formation. • Based on the discrete fracture model, the temperature prediction model of oil-water two-phase flow was established. Considering the capillary force and osmotic pressure, the temperature (T) -hydraulic (H) -mechanical (M) -chemical (C) was realized. The finite element method is used to solve the model. • The effects of multiple factors on reservoir and fracture temperature during fracturing pumping and shut-in stages are analyzed. [ABSTRACT FROM AUTHOR]

The temperature of a deep shale reservoir may reach more than 100°C, and the effect of thermal shock on shale hydraulic fracturing has rarely not been considered in previous studies. Based on mesoscopic damage mechanics and the finite element method, a thermo-hydro-mechanical-damage (THMD) coupling model considering temperature, seepage, stress, and damage fields was constructed to investigate the effects of reservoir temperature, convective heat transfer coefficient (h), in-situ stress difference and bedding plane angle (α θ) on shale hydraulic fracturing. The results show that multiple hydraulic fractures (HFs) can occur under thermal shock and that HFs control the distribution of seepage, temperature, and stress fields. Reservoir temperature, in-situ stress difference and α θ are primary factors affecting hydraulic fracturing, whereas h is a secondary factor. When the reservoir temperature rises from 50°C to 150°C, the initiation and breakdown pressures decrease by 65.5% and 16.7%, respectively. HFs cross the bedding plane more easily, and fracture complexity is obviously enhanced. A higher h is favourable for slightly reducing the initiation and breakdown pressures, but it has little influence on the fracture complexity. Once the in-situ stress difference is low, there is a high fracture complexity, but HFs are more easily captured by bedding planes to limit the propagation of fracture height. When the in-situ stress difference is high, HFs are more likely to form bi-wing fractures. Whether α θ is too large or small, it is not conducive to improving the fracture complexity. In this study, when α θ is 30°, HFs and bedding planes intersect to form a fracture network. Essentially, thermal shock plays a key role in reducing the initiation pressure and forming multiple HFs during the fracturing process, and fracture propagation mainly depends on the injection pressure. The results can serve as reasonable suggestions for the optimization of shale hydraulic fracturing. [ABSTRACT FROM AUTHOR]

An effective evaluation on fracture network forming capacity is key to the whole process of shale gas exploration. At present, neither a clear standard nor a generally accepted evaluation method exist. In this study, a novel “Soundless Cracking Agent (SCA) fracturing evaluation method” was developed. The fractures were characterized quantitatively using fractal dimension of the trace on the core surface and areal density. Acoustic emission (AE) location was used for dynamic monitoring and analysis. The results show that the fractal dimension can be used for quantitative evaluation of complexity of fracture network. The higher the rock hardness, the smaller the fracture density after fracturing is; the higher the brittleness, the larger the fracture density after fracturing is. The development degree of natural fracture systems and sedimentary bedding is a key factor to control the propagating morphology of fractures. The number of AE events for sandstones with low clay content (<25%) is huge, and there are obvious take-off spots for cumulative curves and frequency distribution curves. The AE events for sandstones are distributed along the main fractures, with simple planar fractures clearly present after fracturing. But for shale, the number of AE events is less, with no obvious take-off spots, and AE events are scatteredly distributed. The higher the clay content and the lower the quartz content, the smaller the number of AE events is, and the smaller the frequency and the sound source amplitude are. For sandstone, the number of AE events decreased by about 75% due to the increase of clay content by 20%. The new method enables a comprehensive reflection of the characteristics of rock hardness, brittleness and natural fractures system. This work is valuable for the evaluation of hydraulic fracturing effects in unconventional oil and gas reservoirs in the future. [ABSTRACT FROM AUTHOR]

There is a wealth of experimental evidence that fracture initiation and breakdown pressures differ depending on in-situ stress status, rock properties, and injection conditions. However, the mechanism is not fully understood from a theoretical modeling perspective. In this study, a fully coupled plain-strain fracture model is proposed to interpret the mechanism of fracture initiation and breakdown pressures. The fracture model consists of fracture initiation and propagation governed by linear elastic fracture mechanics. The effects of wellbore compliance (or compressibility), solid-fluid coupling, and fracture multiscale propagation behavior are fully considered. The solid-fluid coupling equations are solved using the Newton–Raphson iterative method. The explicit time marching method is used to capture the fracture initiation process. An implicit time-stepping with the fracture tip asymptotic solution is used to capture fracture propagation fronts. The model is validated against the analytical solutions of the plane-strain model. Sensitivity analysis demonstrates that the initiation pressure mainly depends on rock properties, especially the fracture toughness. The peak pressure (breakdown pressure) is related to rock properties and injection conditions and usually occurs before the peak of the pressurization rate is reached. It increases with the injection rate, fluid viscosity, Young's modulus, and fracture toughness. The dimensionless inlet flux into the fracture can be used to determine the fracture initiation pressure. The pressurization rate during the fracture initiation stage is constant and can be used to assess wellbore compliance. Using a low injection rate and a low-viscosity fluid is beneficial to capturing the fracture initiation pressure. This study can help understand fracture initiation and propagation and interpret hydraulic fracture initiation and breakdown pressures. • The wellbore compressibility, solid-fluid coupling, and multiscale propagation behavior are considered. • Breakdown pressure is dependent on injection conditions and rock properties. • Fracture initiation pressure and wellbore compliance can be determined by inlet flux and pressurization rate. • Low-viscosity fluid and low injection rate favor investigating fracture initiation process. [ABSTRACT FROM AUTHOR]

Highlights: [•] Hydraulic fracturing simulation experiments of shale outcrops were first carried out. [•] Fracture morphology was observed for the first time by high-energy CT scanning. [•] The effects of multiple factors on fractures propagating in shale play were studied. [•] CT scanning images were combined with internal fractures photographs for analysis. [•] Hydraulic fracturing of horizontal well was simulated for shale specimens. [Copyright &y& Elsevier]

Compared with conventional hydrofracturing fluid, SC-CO 2 has the merits of no water sensitivity and effective carbon dioxide storage. At present, the research on the fracture initiation mechanism of SC-CO 2 fracturing is still at the exploring phase. This paper conducts a detailed study on the fracture initiation mechanism of SC-CO 2 fracturing, and establishes a fracture initiation and propagation model of flow-stress-damage coupling according to damage mechanics. A series of numerical simulation studies have been carried out on the rock physical parameters, fracturing construction parameters, natural fracture distribution and other factors that affect the fracture propagation of the reservoir to provide a reference for the SC-CO 2 fracturing construction design. The findings suggest: the lower formation permeability, the less SC-CO 2 filtration loss, and the larger total fracture length and width; SC-CO 2 is easier to improve the reservoir continuity and result in a complex fracture network; The lower the viscosity of SC-CO 2 , the larger total length of fractures, and the more complex the fracture network; As the delivery of pump increases, the width of fracture continues to increase, and the total fracture length decreases; SC-CO 2 has a better effect on communicating natural fractures and porous media, and the fracture network is more complex. • Fracture propagation model with flow-stress-damage (FSD) coupling is developed. • A simulation method of hydraulic fracture propagation in the glutenite reservoirs is introduced. • The model accuracy is validated through comparison between results of physical experiment and numerical simulation. • The effects of multiple factors on fracture propagating in the glutenite reservoirs are analyzed. [ABSTRACT FROM AUTHOR]

Abstract: It is critical to optimize the proppant size for an effective sand control in weakly consolidated formations and to obtain high productivity of medium-to-high permeability formations. Compared with the gravel pack completion, the frac pack has a completely different sand control mechanism and pressure field on the proppant. To achieve effective sand control and stimulation, a treatment should be designed to obtain the maximum conductivity and optimize the proppant size for sand production. At present, there are no unified criteria for proppant-size selection for the frac pack; neither is the equipment for conducting such experiments available. In order to experimentally study the frac pack, we designed a special equipment called the “Frac Pack Sand Control Model” for the experiments in which the proppant size was matched with the formation particle size. By measuring the fracture conductivity, sand production rate, and the size of the particles produced, we investigated the effect of fluid viscosity, flow rate, median grain diameter of formation sand, and proppant concentration on sand control and conductivity. We also studied the proppant embedment in unconsolidated formations and the effects of combinations of proppants of different sizes on sand control and conductivity. Experimental results show that the optimal proppant size for the frac pack is 6.8–9.7 times that of formation sands. Compared with proppant of a fixed size, a combination of proppants of variable sizes is preferred for a better sand control or higher conductivity, depending on the conditions of a particular formation. In a weakly consolidated formation, proppant embedment can be very serious, and becomes even more serious with an increasing proppant size. While ensuring the optimal size of proppant, i.e. D 50=(6.8–9.7)d 50, a high proppant concentration is needed to obtain effective sand control in a frac pack. The proppant concentration not less than 15kg/m2 is usually recommended in a frac pack. The finding in the research is significant for both theoretical study and field engineering. [Copyright &y& Elsevier]

The stimulation method of the marine natural gas hydrate (NGH) reservoir through hydraulic fracturing has been proposed to resolve the problem of the low production capacity in the conventional development method of pressure drawdown. Nevertheless, due to the strong plasticity and high argillaceous siltstone content of the marine NGH reservoir, conventional small-particle-size proppant cannot form effective support for fractures after fracturing because of serious embedding in the reservoir. To solve this problem, the large-size phase change proppants were developed in this study. First, an epoxy resin curing system that can reduce curing time to 40 min in low temperature and humid environment was developed. Then, the epoxy resin and curing system was emulsified, and through the optimization of the emulsification process, the particle size of the proppant can be controlled in 0.5–4.5 mm and the cementation between the proppant particles during the curing process can be prevented. Finally, the proppant performances were evaluated. The performance evaluation shows that the cured proppants have regular structure and good compressive strength, and the emulsion proppants have good transport capacity. Their large sizes provide effective propping effects for fractures generated in weakly cemented clayey silt marine NGH reservoirs. [ABSTRACT FROM AUTHOR]

• A coupled THMD model for thermal reservoir fracturing propagating simulation was built. • The model accuracy was validated through case study, theoretical model and experiment. • The effect of different parameters on HDR hydraulic fracturing under thermal stress was explored. The hot dry rock (HDR) hydraulic fracturing is a complex physical process coupling the effects of seepage, stress, temperature, and damage. The high temperature and brittleness of the HDR leads to the great thermal stress, and the rock is possibly thermally damaged, thus promoting hydraulic fracture (HF) extension and significantly improving the permeability around the HF. In this paper, a thermo-hydro-mechanical-damage (THMD) coupling model is established based on elastic thermodynamics, Biot's classic seepage mechanics and mesoscopic damage mechanics, and its accuracy is evaluated through case study and verification with theoretical models and experiments. The evolution of multi-physics during hydraulic fracturing of HDR is studied, and the effects of rock thermophysical parameters, temperature difference, rock heterogeneity, Young's modulus, permeability, and injection rate on HF extension in the HDR are investigated. The results show that initially, due to the severe temperature variation near the borehole, the higher thermal expansion coefficient leads to the greater thermal tensile stress and facilitates rock damage, thus reducing the fracture pressure. The research results provide theoretical basis and technical support for fracturing design of geothermal system. [ABSTRACT FROM AUTHOR]

• An improved thermal-hydraulic-mechanical coupling model was established. • Evaluation index system for geothermal energy development effect was put forward. • Evaluation method for EGS based on AHP and FCE was proposed. • Well arrangement patter and development parameters were optimized. In the paper, the thermal-hydraulic-mechanical (THM) coupling model is improved by considering more factors, and a comprehensive evaluation method is proposed to optimize the heat extraction way. The thermal extraction characteristics in the Enhanced Geothermal System (EGS) and their influence factors were modeled by several researchers, and various evaluation and optimization methods were established. However, by comparing various thermal-hydraulic-mechanical (THM) models, we find that these factors of THM models are not considered comprehensively, and current EGS evaluation focuses on the individual index, the multi-objective evaluation index system of EGS has not been proposed, and a comprehensive evaluation system considering multiple indexes has not been established yet. In this study, a thermal-hydraulic-mechanical coupling model, considering matrix, fractures and changes in physical parameters of rock and fluid, was established for optimizing the EGS. The evaluation index system, including the operation life, the thermal breakthrough time, the average heat production rate, the total heat recovery factor and the re-injection factor, were presented. And the comprehensive optimization method based on analytic hierarchy process (AHP) and fuzzy comprehensive evaluation (FCE) of EGS were proposed. The model is verified by two analytical solutions. Moreover, the model was solved by finite element method (FEM) to optimize various programs of developing the EGS. In this paper, the performances of different geothermal extraction methods are compared and the different extraction parameters are optimized. The results show that the optimized geothermal mining method has better effect than that before optimization. The orthogonal test shows the optimal case of 5 fractures, the fracture permeability of 10D, the well spacing is 400 m and the production pressure draw-down of 20 MPa. The research results provide an accurate and effective method to optimize the development mode of geothermal resources, which is conductive to improve the geothermal energy efficiency. [ABSTRACT FROM AUTHOR]

Hydraulic fracturing technology with horizontal wells is a highly efficient approach to stimulate the unconventional reservoir by creating lots of fractures to increase the oil and gas recovery rate. However, some field monitoring data indicate that multiple fractures could not propagate uniformly due to stress interference, which leads to a lower stimulation volume and production capacity. Therefore, investigating the stress interference mechanism between multiple fractures is important to obtain a higher production capacity for unconventional reservoirs. In this study, a two-dimensional XFEM (expanded finite element method) model was established to investigate the influence of stress interference on the fracture propagation of multiple fractures in horizontal wells. The different injection parameters, completion parameters, and fracturing patterns were investigated by sixteen simulation cases. The results show that first the influencing degree of stress interference from different parameters is fracture spacing, injection flow rate, number of fractures, and fluid viscosity, respectively. Second, compared to the simultaneous fracturing and sequential fracturing, the zip fracturing pattern (sides fractures created first and then mid-fracture) could create more uniform fractures in the reservoir. Third, a small fracture spacing and a higher injection flow rate could lead to a faster stress interference, but the fluid viscosity and the number of fractures have little effect on the occur timing of stress interference. This study could provide some meaningful perspective on the field fracturing design in horizontal wells. [ABSTRACT FROM AUTHOR]

Formation of enhanced geothermal system provides flow and heat transfer paths for heat extraction in HDR (hot dry rock). Based on the local thermal non-equilibrium theory, the THM(thermal-hydraulic-mechanical) coupling is established to describe the interaction of fluid flow, heat transfer and rock deformation in heat mining process. Besides, the randomly generated fractures are adopted to descript the various fracture network. Firstly, the heat mining processes of variably fractured HDR with different fracture network characteristics are researched. Secondly, heat mining performance of multi-well patterns are compared with that of doublet-well pattern. Finally, effect of the key parameters on heat mining are also investigated. The results indicate that the increase of fracture density enhances flow around and retards thermal drawdown. When the angle between fracture orientation and inlet-outlet connection is 45°, the optimal heat mining performance can be achieved. The addition of inlets in multi-well pattern conduces to large-scale utilization, and a reasonable well arrangement can avoid the thermal breakthrough. This study provides the support for EGS (enhanced geothermal system) stimulation and well arrangement to obtain better heat mining performance. • Based on the local thermal non-equilibrium theory the THM coupling is established. • Various fracture network morphologies of EGS are characterized by randomly generated fractures. • Effect of fracture network density and orientation on heat mining are studied. • Multi-well patterns are compared with doublet-well pattern in EGS exploitation. • Long-time outlet temperature, outlet mass flow rate and heat mining rate are researched. [ABSTRACT FROM AUTHOR]

Abstract Conventional hydraulic fracturing fails to develop low permeability reservoirs of fossil hydrogen energy that are not located in the direction of maximum principal in-situ stress. A new technology of fracture propagation guided by radial holes is proposed, which can realize directional propagation of hydraulic fracture along radial holes in fossil hydrogen energy development. In order to verify this new technology, a model of radial holes combined with hydraulic fracturing is established by the ABAQUS extended finite element method. Simulation results show that radial holes play a guiding role in fractures propagation. The influence extent of seven factors on the directional propagation of hydraulic fracture is listed as follows (from strong to weak): azimuth of radial holes > horizontal in-situ stress difference of fossil hydrogen reservoir > injection rate of fracturing fluid > Young's modulus of rock > permeability of fossil hydrogen reservoir > Poisson ratio of rock > viscosity of fracturing fluid. True tri-axial experiment is carried out to verify the accuracy of numerical simulation, and the result is consistent with numerical model, which indicates that numerical simulation is reliable. Highlights • A three-dimension stress coupled seepage model of radial holes is established. • Numerical model proves that radial holes play a guiding role in the directional propagation of hydraulic fractures. • The influence of seven factors on directional propagation of hydraulic fracture guided by radial holes is analyzed. • Tri-axial hydraulic fracturing experiment is carried out to verify the accuracy of numerical simulation results. [ABSTRACT FROM AUTHOR]

The simulation of hydraulic fracturing by the conventional ABAQUS cohesive finite element method requires a preset fracture propagation path, which restricts its application to the hydraulic fracturing simulation of a naturally fractured reservoir under full coupling. Based on the further development of a cohesive finite element, a new dual-attribute element of pore fluid/stress element and cohesive element (PFS-Cohesive) method for a rock matrix is put forward to realize the simulation of an artificial fracture propagating along the arbitrary path. The effect of a single spontaneous fracture, two intersected natural fractures, and multiple intersected spontaneous fractures on the expansion of an artificial fracture is analyzed by this method. Numerical simulation results show that the in situ stress, approaching angle between the artificial fracture and natural fracture, and natural fracture cementation strength have a significant influence on the propagation morphology of the fracture. When two intersected natural fractures exist, the second one will inhibit the propagation of artificial fractures along the small angle of the first natural fractures. Under different in situ stress differences, the length as well as aperture of the hydraulic fracture in a rock matrix increases with the development of cementation superiority of natural fractures. And with the increasing of in situ horizontal stress differences, the length of the artificial fracture in a rock matrix decreases, while the aperture increases. The numerical simulation result of the influence of a single natural fracture on the propagation of an artificial fracture is in agreement with that of the experiment, which proves the accuracy of the PFS-Cohesive FEM for simulating hydraulic fracturing in shale gas reservoirs. [ABSTRACT FROM AUTHOR]

Hydraulic fracturing of radial wells have become new technologies to effectively develop low permeability reservoir, thin reservoir, fractured reservoir, “dead area” after water injection etc. Influenced by ground stress and reservoir heterogeneity, hydraulically created fracture will probably not extend smoothly to connect the target area along the direction of radial wells, the reservoir stimulation effect is not ideal. Based on the pressure analysis of fracturing initiation of multiple radial holes and the theory of plasticity district, criterion of multiple radial wells orientating directional fracture propagation in the condition of ground stress is derived. In this condition, multiple radial wells will produce a continuous plasticity district in the reservoirs, while maintaining the fracture extending in plastic zone directionally without the effect of ground stress, which guarantees that fracture is interconnected between the radial wells during fracturing operation, to form the main crack along the direction of radial wells’ axis plane. In consideration of economy, maximum spacing of radial wells are optimized, in order to provide a reliable scientific basis for effectively implementation of hydraulic fracturing technology in radial wells. [ABSTRACT FROM AUTHOR]

The Stokes diameter is commonly used to characterize the size of non-spherical particles, but its application is limited to Newtonian fluids. In view of this case, an equivalent hydrodynamic diameter, defined as the effective settlement diameter, has been proposed from force balance for particles settling in non-Newtonian fluids, and an iterative calculation procedure was presented. In the present work, proppant effective settlement diameters falling in non-Newtonian fracturing fluids were studied. Settling velocities of three kinds of proppants falling in glycerol solutions, uncross-linked guar solutions and slick-waters were measured, numerous random repeated experiments were conducted, the corresponding proppant effective settlement diameters were calculated and the relationships between proppant effective settlement diameter and their sieve diameter falling in different fluids were studied. Results indicate that the proppant effective settlement diameter falling in uncross-linked guar solutions is bigger than that in glycerol, while it is smaller for slick-waters. When proppants settle in glycerol solutions, the proppants effective settling diameters are constantly 1.24 times bigger than their corresponding sieve diameter. For proppants settling in uncross-linked guar solutions and slick-waters respectively, two sigmoid curves of diameter ratio are obtained followed by reaching an asymptotic value, 1. This work presented a new calculation procedure of equivalent hydrodynamic diameter for non-spherical particles settling in non-Newtonian fluids. The effective settlement diameter extends the application conditions of Stokes diameter to non-Newtonian fluids and can be used in wider engineering applications. [ABSTRACT FROM AUTHOR]

Hydraulic fracturing (HF) has been proved to be a promising technology to achieve economic production of hydrate. However, the research on HF of hydrate reservoirs is at its initial stage with limited information being available. In particular, the mechanism of fracture propagation in hydrate reservoirs is not well understood and requires further investigations. In this study, a new coupled thermo-hydrologic-mechanical-damage (THMD) model for HF simulations in hydrate reservoirs is proposed. Damage mechanics theory is used as the criterion of hydraulic fracture initiation and propagation, and the variation of hydrate properties due to hydrate phase transformation is considered in this THMD coupled model. The influence of hydrate saturation, reservoir permeability, fracturing fluid viscosity and fluid injection rate on HF were analyzed, and the mechanism of fracture initiation and propagation during HF in hydrate reservoir was revealed for the first time. The results showed that fracturing fluid destroys phase equilibria and causes hydrate dissociation, which in turn, releases the pore spaces occupied by hydrate, resulting in the increase of reservoir permeability near fracture surface. Besides, the hydrate dissociation reduces its cementation on sediment particles, causing the decrease of cohesion near the fracture surface. The fracture initiates and propagates perpendicularly to the direction of minimum principle field stress. Hydrate dissociation leads to heterogeneity of reservoir physico-mechanical properties, resulting in a irregular fracture surface morphology. Enhancing the fracturing fluid viscosity and injection rate promotes to improve the reservoir pore pressure and inhibit methane hydrate dissociation, which is beneficial to increase the fracture length. The hydrate reservoir with high methane hydrate saturation and low permeability is conductive to form long fractures during HF. • A new thermo-hydrologic-mechanical-damage (THMD) coupling model is proposed. • The dynamic property of reservoir caused by hydrate decomposition and damage is considered. • The mechanism of fracture initiation and propagation during hydraulic fracturing is revealed. [ABSTRACT FROM AUTHOR]

As an important part in the process of hydraulic fracturing, fracturing fluid flowback after fracturing will directly affect the productivity of the fractured wells. It is of great significance to carry out optimization research on the flowback working system of fracturing fluid (the time of flowback beginning and the size of flowback nozzles) to enhance the efficiency of fracturing. At present, researches on optimization model of fracturing fluid flowback have taken into account the effects of fracturing fluid loss and proppant reflux, while ignoring the effects of proppant sedimentation. A new calculation model of proppant regurgitation velocity has been established in recent years, but it is rarely used in the research of flowback working system. This study innovatively integrated the calculation model of modified proppant sedimentation velocity obtained through experiments, the new calculation model of critical proppant regurgitation velocity, and the wellhead pressure calculation model to form a comprehensive flowback model of fracturing fluid. Then coupled with the flowback criteria of quickly flowback, reducing proppant reflux, and maintaining proppant settlement degree (ratio of proppant settlement distance to crack height) less than 60%, an optimization model of fracturing fluid flowback is constructed and realized by MATLAB software programs. The optimization model were verified by four fractured wells in DaQing Oilfield. The effects of reservoir permeability, reservoir porosity, fracturing fluid density, and proppant particle size on wellhead pressure drop, fracturing fluid flowback rate, and proppant settlement degree were calculated and analyzed by using the above flowback model. The verification results show that the simulated wellhead pressure and flowback flowrate are in good agreement with the actual dataset, with the average errors of 11.9% and 11.6% respectively, which proves the accuracy of the flowback optimization model. Parametric analysis shows that, with the increase of reservoir permeability and reservoir porosity, the wellhead pressure drop velocity increases. When the permeability increases from 1mD to 5mD, the flowback rate of fracturing fluid and proppant settlement degree are reduced by 12.3% and 69.8% respectively. When the porosity increases from 9% to 36%, the proppant settlement degree fall by 25.2%. Moreover, the flowback rate of fracturing fluid decreases linearly and the proppant settlement degree increases linearly with the increase of fracturing fluid density; and the proppant settlement degree increases by 3.5 times when the proppant particle size increases from 0.25 mm to 1.65 mm. Therefore, choosing smaller fracturing fluid density and smaller proppant particle size can ensure higher fluid flowback efficiency and lower proppant settlement degree. • A new flowback model of fracturing fluid is constructed to optimize the flowback working system, which takes the influence of proppant settling into account. • The optimization model is realized by programming with MATLAB software and its accuracy is verified by four fracturing wells of Daqing Oilfield. • Then the influences of reservoir permeability, reservoir porosity, fracturing fluid density and proppant particle size on the wellhead pressure drop, fracturing fluid flowback rate and proppant settlement degree are clarified. [ABSTRACT FROM AUTHOR]

To realize the effective extraction of high-temperature thermal energy in enhanced geothermal system (EGS), it is necessary to explore how to efficiently construct fracture network in hot dry rock (HDR). Considering the cryogenic induced thermal stress, we investigate the cracking mechanism when using supercritical CO 2 (SCCO 2) and hydraulic fracturing to stimulate HDR. Firstly, the established three-dimensional THMD coupling model is verified by conducting high temperature granite fracturing experiments. Then, characteristics of the heat and mass transfer and fracture propagation when using SCCO 2 and H 2 O in HDR fracturing are explored. Finally, the cracking mechanism of hydraulic and SCCO 2 fracturing in HDR under thermal stress is revealed. The results indicate that the cryogenic induced thermal stress can reduce the cracking pressure and tend to form branch fractures under the cooperation of injection pressure. Adopting SCCO 2 in HDR fracturing more micro fractures can be generated in near-well zone. Hydraulic fracturing has better cooling efficiency and preferably capability of extending fractures. For HDR fracturing, the scope of fracture network tends to increase first and then decrease with the raise of injection mass flux. The lower viscosity and higher specific heat capacity of fracturing fluid can promote the formation of fracture network. • True-triaxial fracturing tests on granites are conducted under different temperature. • Three-dimensional THMD coupling model is established to simulate hot rock fracturing. • Damage evolution of hydraulic and SCCO 2 fracturing in HDR stimulation is investigated. • Effect of injection mass flux on fracture morphology in HDR fracturing is researched. • Lower viscosity and higher heat capacity of fracturing fluid can promote the formation of EGS. [ABSTRACT FROM AUTHOR]

Thermal fluid injection is proven to be one of the effective methods to dissociate hydrate. The injection of thermal fluid can destroy the hydrate phase equilibrium, resulting in hydrate decomposition and hydrate saturation decline, which further cause reservoir damage. Previous works separated the study of hydrate saturation variation caused by thermal fluid injection and the reservoir damage caused by hydrate saturation variation. In this paper, a novel thermo-hydrologic-mechanical (THM) coupling model is proposed to simulate the thermal fluid injection in hydrate reservoir, and its effects on reservoir damage were analyzed. The results showed that thermal fluid injection causes the decomposition of hydrate, resulting in the decrease of hydrate saturation. The solid phase hydrates decompose into methane and water, which increases the pore pressure of hydrate reservoir and further affects the in-situ stress field. The change of effective stress induces pore deformation, and the hydrate decomposition also releases the space occupied by hydrate in reservoir pores, both affecting the reservoir permeability. Hydrate phase transition weakens the cementation of hydrate on rock particles and reduces the reservoir rock cohesion. In addition, the elastic modulus of hydrate decreases due to hydrate decomposition. The temperature and injection rate of thermal fluid are the main external factors affecting the hydrate phase transition. With the increase of thermal fluid temperature, hydrate saturation decreases due to hydrate decomposition, resulting in variation of reservoir parameters. Enhancing fluid injection rate can increase the reservoir pore pressure, which inhibits the hydrate decomposition. • A new thermo-hydrologic-mechanical (THM) coupling model is proposed. • The dynamic property of reservoir caused by thermal fluid injection is considered. • The mechanism of hydrate phase transition during thermal fluid injection is analyzed. • Increasing injection rate and reducing fluid temperature can reduce reservoir damage. [ABSTRACT FROM AUTHOR]

• In fluence of stress shadow on multi-clusters fracturing was analyzed. • The proppant transport and placement were considered in the model. • Multiple cluster hydraulic fracturing in multiple horizontal wells was simulated. Multiple cluster hydraulic fracturing is widely used in the development of shale reservoir. The morphology of multiple fractures propagating from multiple cluster is complex due to the stress shadow effect. The design and operation of multiple cluster hydraulic fracturing in horizontal wells require adequate knowledge of the effect of different factors, including rock properties, in-situ stresses and fluid properties on the fracture morphology. In this paper, lattice simulation, a new particle based computational method was used to investigate the multiple cluster hydraulic fracturing in shale. The results showed that stress anisotropy and cluster spacing play an important role on geometry of the propagating fracture. The middle hydraulic fracture is restricted to propagate when cluster spacing is decreased. Simultaneous, two-step and sequential fracturing scenarios in a single horizontal well were simulated. The simultaneous and sequential fracturing showed to mainly affect the fracture propagation morphology with little effect on fracture length, while the middle fracture is shorter than its two side fractures in two-step fracturing. In two horizontal models, the simultaneous and sequential fractures showed similar morphology during multiple cluster hydraulic fracturing. In term of stimulated reservoir volume (SRV), zipper fracturing showed the largest SRV during multiple cluster hydraulic fracturing in two horizontal wells. [ABSTRACT FROM AUTHOR]

• A new model is proposed to study the multi-cluster fracturing. • The influence of different factors on the fracture propagation is studied. • The accuracy of model is verified by experiment. Multi-cluster fracturing in horizontal wells is a key technology for successful development of ultra-low permeability reservoir. The propagation of hydraulic fracture during multi-cluster fracturing is complicated, especially in shale reservoir with multiple natural fractures. The design and operation of multi-cluster fracturing requires adequate understanding of influence of different factors on hydraulic fracture propagation. Up to now, many scholars have studied the hydraulic fracture morphology in multi-cluster fracturing, but few have analyzed the effect of natural fractures on hydraulic fracture propagation during multi-cluster fracturing. In this paper, a new computationally version of the particle-based model is established by Xsite to study the fracture propagation in multi-cluster fracturing with natural fractures. The tensile strength and rock toughness are calculated, and tri-axial experiments are performed to verify the accuracy of model. Simulation results show that cluster spacing and in-situ stress difference have a significant influence on the length of the hydraulic fracture and the morphology of fracture. The length of middle fracture increases with the increase of the cluster spacing, but decreases with the increase of the in-situ stress difference during multi-cluster fracturing with three natural fractures. The enhancing of cluster spacing can reduce the deflection of left and right fractures, and the increase of the in-situ stress difference can improve the ability of middle fracturing penetrating the natural fracture. Three fracturing sequence of synchronous fracturing, two-step fracturing and sequential fracturing is simulated. The left and right fractures can always penetrate the natural fracture with different fracturing sequence. But the middle fracture shows different morphology of arresting by natural fracture (during synchronous fracturing), partially penetrating the natural fracture (during two-step fracturing) and directly penetrating the natural fracture (during sequential fracturing). Different cement strengths of natural fracture are analyzed. The increase of strength of natural fracture can enhance the ability of hydraulic fracture penetrating the natural fracture. Hydraulic fracture is arrested by weak natural fracture (shear strength of 0.5 MPa), resulting in natural fracture opening. As hydraulic fracture intersecting with strong natural fracture (shear strength of 20 MPa), the hydraulic fracture can penetrate the natural fracture directly without natural fracture opening. [ABSTRACT FROM AUTHOR]

Water cut is generally rising with the increase of water injection during water flooding development of fossil hydrogen energy reservoir. Conventional chemical profile control agents fail to work for a long time because they are easily affected by the temperature and salinity of fossil hydrogen energy reservoir to form crossflow channel. In this paper, a new dual-disperse profile control agent and flooding system (DDS) is proposed, which uses polymers and microspheres to plug polymer crossflow and improve oil recovery in fossil hydrogen energy development. The viscosity, suspension property, thermal stability and microsphere size expansion of the DDS are tested, and the oil displacement performance is evaluated. The results show that the viscosity of polymer solution is not affected by the addition of microspheres, and the DDS system has good suspension and thermal stability. Due to the adsorption of microspheres in polymers, the particle size of microspheres tends to increase with time. But the expansion slows down gradually, which is beneficial to deep migration. In terms of plugging efficiency, the optimum concentration of DDS is 1500 mg/L polymer and 500 mg/L microsphere, and the optimum injection volume is 0.5Pv. Compared with single polymer flooding system, the new dual dispersion system has better washing resistance and enhanced oil recovery in fossil hydrogen energy development. • A new DDS dual dispersion profile control system is proposed. • The physical property of DDS is tested. • The oil displacement performance of DDS is evaluated. • DDS can increase oil recovery significantly. [ABSTRACT FROM AUTHOR]

The existence of solid hydrate in porous media always greatly reduces the relative permeability, which limits the development of hydrate reservoirs. Therefore, to recover the gas, it is necessary to carry out reservoir reconstruction measures that induce fractures around the production wells to increase the gas production of hydrate reservoirs. At present, there is limited quantitative research on the promotion of gas production by induced fractures in hydrate reservoirs. In this paper, numerical models of Class Ⅲ hydrate reservoirs with induced fractures are established. By comparing the gas production with and without induced fractures, the effectiveness of induced fractures in promoting the depressurization production of natural gas hydrates is verified. An index of improved efficiency is established and calculated to measure the promotion effect of induced fractures during the fracture promotion stage. The multiple regression calculation formulas of improved efficiency with reservoir geological parameters are obtained by integrating nonlinear regression and linear regression. The calculation model is helpful to quickly and directly measure the rationality and economy of the induced fracture depressurization method to promote hydrate development. • Numerical models of Class Ⅲ hydrate reservoirs with induced fractures are established. • An index of Improved Efficiency is established to measure the promotion effect of induced fractures. • The improved efficiency of induced fractures under different geological parameters is calculated and compared. • The multiple calculation formulas of improved efficiency with geological parameters are obtained by stepwise regression. [ABSTRACT FROM AUTHOR]

Highlights • An effective method for simulating wormhole in fractured formation was proposed. • The effects of three typical fractures on wormhole propagation were studied. • The local control domain around the fracture was found. • A step-by-step calculation method was proposed in this paper. Abstract Wormholes can effectively connect the fracture system as the dominant channels when carbonate matrix is acidified, thereby increasing the distance of acidification. It is critical to clarify the wormhole propagation law in carbonate reservoirs so as to optimize the acidification design. Previous studies highlight the influences of matrix as well as single simple fracture on wormhole propagation. However, the works characterizing impacts of single and multiple fractures with different geometry are insufficient. In this paper, the two-scale continuum model and pseudo-fracture model were combined to study the wormhole with different fractures. It was found that the fractures parallel to injection direction can concentrate fluids and thus accelerate the acid penetrating the formation, while the perpendicular ones will disperse fluids, thereby slowing down the breakthrough. The effect of straight fractures and arc fractures on wormholes are similar. When the inclination angle is less than 60°, straight fracture and arc fracture can be regarded as parallel fracture of corresponding length. When the inclination angle is greater than 60°, they can be regarded as superposition of parallel fracture and perpendicular fracture. As for the circular and polygonal fractures, only the parts near the outlet impact the subsequent wormhole propagation trajectories. The analysis of flow field shows that there is a control domain of about 3· L ·cos(θ) around the fractures. When there is no intersection of control domains of multiple fractures, there is no interference between them. There is no need to consider all fractures in the formation. Based on this, the step-by-step calculation method proposed in this paper is suitable for acidizing large-scale formation with complex fracture network. [ABSTRACT FROM AUTHOR]

Abstract At present, the proportion of tight oil in non-renewable hydrogen energy is increasing. According to an initial exploration and attemptable practice on the exploration of tight oil, it is found that the cost can be controlled effectively and positive effects are achieved. But this technique cannot make sure the proppants filled uniformly in the long fracture. Several researches on the proppants migration experiment devices and factors influencing on proppant setting are reviewed and a new set of experimental device to simulate the laws of proppants setting in long fracture is developed. This device can simulate the main factors influencing proppants setting performance. It analyzes several factors such as wall filtration, construction displacement, sand concentration, proppant size and density, viscosity of fracturing fluid is used to rank the influencing degree of every factor. Considering the effects of mutual interference between proppants, width of fracture, rough fracture surface and fracture surface filtration during the proppants setting progress, the mathematical model of proppant setting is modified by adding sand concentration correction factor, wall effect correction factor and filtration correction factor. The experimental data verify the accuracy of the settlement model is established using the data getting from experiment. Highlights • Simulated the factors influencing proppant placement, including wall filtration. • The mathematical model of proppant placement is modified by adding filtration correction factor. • The influence of different factors on proppant placement is quantitatively analyzed. [ABSTRACT FROM AUTHOR]

• Water-gas relative permeability in hydrate-bearing porous media is studied by LBM. • Relationship between relative permeability and fluid distribution is characterized. • Effect of hydrate on water-gas relative permeability is revealed at a pore scale. The water-gas relative permeability is an important parameter to characterize multiphase flow in sediments. To study the water-gas relative permeability of hydrate-bearing porous media, multiphase flow simulations were carried out at the pore scale using the lattice Boltzmann method. The effects of hydrate saturation and hydrate-growth habits on the water-gas relative permeability, which is scaled by the relative permeability considering the hydrate only, were evaluated in a two-dimensional porous medium. Results show that the increase of hydrate saturation causes the decrease of water-gas effective permeability as expected. However, the effect of hydrate saturation on the water-gas relative permeability is different from that of hydrate saturation on the water-gas effective permeability. The water-gas relative permeability increases with the increase of hydrate saturation in the pore-filling case. The water-gas relative permeability decreases with the increase of hydrate saturation in the grain-coating case. The wettability of solid phase has a different effect on the relative permeability of wetting phase and nonwetting phase. The Jamin effect (phase blocking) was observed and may exist in the production of gas from natural gas hydrate reservoirs. This seriously affects the multiphase flow characteristics. The changes of microscale fluid distribution effect the changes of water-gas relative permeability. The relationship between the water-gas relative permeability and the characterization parameters of microscale fluid distribution was analyzed. [ABSTRACT FROM AUTHOR]

To understand the fracture behavior for variable density in-plane perforations on a horizontal well, true triaxial hydraulic fracturing experiments were performed with acoustic emission monitoring. The experimental results indicate that the variable density in-plane perforations can not only create transverse fractures but also arrest fracture propagation toward undesirable zones. The in situ stress and treatment parameters play significant roles in the fracture morphology and breakdown pressure. Under a large horizontal principal stress difference, there is a high possibility for simple transverse fracture creation, which results in a lower breakdown pressure. Moreover, the decrease in fracturing fluid viscosity can increase the fracture complexity and increase in pump rate can increase breakdown pressure and induce the complex fracture geomtetry. The acoustic emission (AE) characteristics and cutting fracture morphology demonstrate that not all perforations can be successfully initiated during the fracturing process. Because of the stress interaction from side perforations, fracture initiation and propagation from the middle perforation are strongly suppressed. Although the middle fracture is not fully initiated and propagated, it is still beneficial for the side hydraulic fracture connection and propagation in one plane. Secondary and axial fractures can still be observed, primarily due to the stress interaction between neighboring perforations and misalignment of the perforation tunnel orientation with the principal stress direction. In addition, the creation of complex fractures tends to occur using a high pump rate, which relates to tense fluid distribution competition from each perforation tunnel. Fracture initiation may follow an order but mostly tends to initiate from side perforations. Because of the stress interaction, hydraulic fractures from the side perforation can sometimes deviate from the perforation tunnel direction, twist and kink, and finally, generate nonplanar fractures. A large horizontal principal stress difference is recommended and beneficial for simple, straight and planar fracture creation during hydraulic fracturing stimulation with variable density in-plane completion for the horizontal well. • The application of local in-plane perforations can control the fracture path and complexity. • Tri-axial fracturing experiments were conducted on tight rock samples with local in-plane perforations. • The change in fracturing fluid viscosity and pump rate can alter the fracture path between perforations. • The connectivity of neighboring fractures can lower the breakdown pressure. [ABSTRACT FROM AUTHOR]

• The reason for the poor effect of acid fracturing with retarded acid in dolomite reservoir was found. • The correlation between the acid-rock reaction control mode and etching morphology was analyzed. • The effects of different factors on the conductivity of fractures were studied. • Two uneven acid injection approaches were discussed. Retarded acids are often used for acid fracturing in carbonate reservoirs to reduce acid leakoff. Compared with limestone reservoirs, dolomite reservoirs are seen with lower fracture conductivity using retarded acid. This study focuses on the acid etching morphology of fracture, and investigates influential factors and improvement methods of fracture conductivity. Our analysis shows the etching morphology has no direct dependency upon the controlling regime of acid-rock reaction. It is determined by reaction rates. Low conductivity is caused by the slow reaction between dolomite and retarded acid, and difficulties in creating irregular etched grooves on fracture surfaces. A small amount of dolomite in reservoirs does not impact conductivity, and yet dolomitic contents over 80% will lead to tremendous reduction of conductivity. It is feasible to form the irregular etching morphology by reducing the flow rate and the thickening agent concentration while alternately injecting two types of acid solutions into reservoirs. [ABSTRACT FROM AUTHOR]

Slickwater fracturing can create complex fracture networks in shale. A uniform proppant distribution in the network is preferred. However, proppant transport mechanism in the fracture network is still uncertain, which restricts the optimization of sand addition schemes. In this study, slot flow experiments are conducted to analyze the proppant placement in the complex fracture system. Dense discrete phase method is used to track the particle trajectories to study the transport mechanism into the branch. The effects of the pumping rate, sand ratio, sand size, and branch angle and location are discussed in detail. Results demonstrate that: (1) under a low pumping rate or coarse proppant conditions, the dune development in the branch depends on the dune geometry in the primary fracture, and a high proportion of sand can transport into the branch; (2) using a high pumping rate or fine proppants is beneficial to the uniform placement in the fracture system; (3) sand ratio dominates the proppant placement in the branch and passing-intersection fraction of a primary fracture; (4) more proppants may settle in the near-inlet and large-angle branch due to the size limit. Decreasing the pumping rate can contribute to a uniform proppant distribution in the secondary fracture. This study provides some guidance for the optimization of proppant addition scheme in the slickwater fracturing in unconventional resources. [ABSTRACT FROM AUTHOR]

Temporarily plugging staged fracturing (TPSF) can effectively increase the complexity of fracture network in naturally fractured reservoirs. In different geological and engineering conditions, how to evaluate the stress interference of multiple fractures and interaction of natural fractures should be investigated further. To realize the temporarily plugging, new perforation elements were developed, and modified pore pressure cohesive zone (PPCZ) model was combined to simulate the interaction between hydraulic fractures and natural fractures. All simulated results are validated by comparing with analytical model. Results show that lateral fractures can really break through the suppression from previous fractures and enhance the complexity of fracture network when using TPSF. Then, case studies are presented to investigate the influence of cluster spacing, horizontal stress difference and number of perforation clusters. The results reveal that, when increasing the cluster spacing, fractures propagate uniformly and influence less; when the spacing reaches at a certain degree (40 m), the complexity of fracture network will not change. Higher stress difference results in simpler fracture network, which owns a straighter propagation trajectory. Further, more perforation clusters lead to stronger interference among multiple fractures, and more divergence and coalescence will take place. [ABSTRACT FROM AUTHOR]