Using the visualized experimental device of temporary plugging in hydraulic fractures, the plugging behaviors of temporary plugging particles with different sizes and concentrations in hydraulic fractures were experim...Using the visualized experimental device of temporary plugging in hydraulic fractures, the plugging behaviors of temporary plugging particles with different sizes and concentrations in hydraulic fractures were experimentally analyzed under the conditions of different carrier fluid displacements and viscosities. The results show that the greater the carrier fluid viscosity and displacement, the more difficult it is to form a plugging layer, and that the larger the size and concentration of the temporary plugging particle, the less difficult it is to form a plugging layer. When the ratio of particle size to fracture width is 0.45, the formation of the plugging layer is mainly controlled by the mass concentration of the temporary plugging particle and the viscosity of the carrier fluid, and a stable plugging layer cannot form if the mass concentration of the temporary plugging particle is less than 20 kg/m^(3)or the viscosity of the carrier fluid is greater than 3 mPa·s. When the ratio of particle size to fracture width is 0.60, the formation of the plugging layer is mainly controlled by the mass concentration of the temporary plugging particle, and a stable plugging layer cannot form if the mass concentration of the temporary plugging particle is less than 10 kg/m^(3). When the ratio of particle size to fracture width is 0.75, the formation of the plugging layer is basically not affected by other parameters, and a stable plugging layer can form within the experimental conditions. The formation process of plugging layer includes two stages and four modes. The main controlling factors affecting the formation mode are the ratio of particle size to fracture width, carrier fluid displacement and carrier fluid viscosity.展开更多
Hydraulic fracturing is considered the main stimulation method to develop shale gas reservoirs. Due to its strong heterogeneity, the shale gas formation is typically embedded with geological discontinuities such as be...Hydraulic fracturing is considered the main stimulation method to develop shale gas reservoirs. Due to its strong heterogeneity, the shale gas formation is typically embedded with geological discontinuities such as bedding planes and natural fractures. Many researchers realized that the interaction between natural fractures and hydraulic fractures plays a crucial role in generating a complex fracture network. In this paper, true tri-axial hydraulic fracturing tests were performed on polymethyl methacrylate (PMMA), on which pre-existing fracture was pre-manufactured to simulate natural fracture. A cohesive model has been developed to verify the results of the experimental tests. The key findings demonstrate that the experimental results agreed well with the numerical simulation outcomes where three main interaction modes were observed: crossing;being arrested by opening the pre-existing fracture;being arrested without dilating the pre-existing fracture. Crossing behavior is more likely to occur with the approaching angle, horizontal stress difference, and injection rate increase. Furthermore, the higher flow rate might assist in reactivating the natural fractures where both sides of the pre-existing fractures were reactivated as the injection rate increased from 5 to 20 mL/min. Additionally, hydraulic fractures show a tendency to extend vertically rather than along the direction of maximum horizontal stress when they are first terminated at the interface. This research may contribute to the field application of hydraulic fracturing in shale gas formation.展开更多
The effective plugging of artificial fractures is key to the success of temporary plugging and diverting fracturing technology,which is one of the most promising ways to improve the heat recovery efficiency of hot dry...The effective plugging of artificial fractures is key to the success of temporary plugging and diverting fracturing technology,which is one of the most promising ways to improve the heat recovery efficiency of hot dry rock.At present,how temporary plugging agents plug artificial fractures under high temperature remains unclear.In this paper,by establishing an improved experimental system for the evaluation of temporary plugging performance at high temperature,we clarified the effects of high temperature,injection rate,and fracture width on the pressure response and plugging efficiency of the fracture.The results revealed that the temporary plugging process of artificial fractures in hot dry rock can be divided into four main stages:the initial stage of temporary plugging,the bridging stage of the particles,the plugging formation stage,and the high-pressure dense plugging stage.As the temperature increases,the distribution distance of the temporary plugging agent,the number of pressure fluctuations,and the time required for crack plugging increases.Particularly,when the temperature increases by 100℃,the complete plugging time increases by 90.7%.展开更多
Depletion-induced stress change causes the redistribution of stress field in reservoirs,which can lead to the reorientation of principal stresses.Stress reorientation has a direct impact on fracture propagation of inf...Depletion-induced stress change causes the redistribution of stress field in reservoirs,which can lead to the reorientation of principal stresses.Stress reorientation has a direct impact on fracture propagation of infill wells.To understand the effect of stress reorientation on the propagation of infill well’s fractures,an integrated simulation workflow that combines the reservoir flow calculation and the infill well hydraulic fracturing modeling is adopted.The reservoir simulation is computed to examine the relationship between the extent of stress reversal region and reservoir properties.Then,the hydraulic fracturing model considering the altered stress field for production is built to characterize the stress evolution of secondary fracturing.Numerical simulations show that stress reorientation may occur due to the decreasing of the horizontal stresses in an elliptical region around the parent well.Also,the initial stress difference is the driving factor for stress reorientation.However,the bottom hole pressure,permeability and other properties connected with fluid flow control timing of the stress reorientation.The decrease of the horizontal stresses around the parent well lead to asymmetrical propagation of a hydraulic fracture of the infill well.The study provides insights on understanding the influence of stress reorientation to the infill well fracturing treatment and interference between parent and infill wells.展开更多
Pre-driven longwall retracement roadway(PLRR)is commonly used in large mine shaft.The support crushing disasters occur frequently during the retracement,and roof management is necessary.Taking the 31107 panel as resea...Pre-driven longwall retracement roadway(PLRR)is commonly used in large mine shaft.The support crushing disasters occur frequently during the retracement,and roof management is necessary.Taking the 31107 panel as research background,the roof breaking structure of PLRR is analyzed.It is concluded that the roof cutting with vertical hydraulic fracture(HF)at a specified position,that is,fixed-length roof cutting,can reduce support load and keep immediate roof intact.The extended finite element method(XFEM)is applied to simulate hydraulic fracturing.The results show that both the axial and transverse hydraulic fracturing cannot effectively create vertical HFs.Therefore,a novel construction method of vertical HF based on the stress shadow effect(SSE)is proposed.The stress reversal region and HF orientation caused by the prefabricated hydraulic fracture(PF)are verified in simulation.The sub-vertical HFs are obtained between two PFs,the vertical extension range of which is much larger than that of directional hydraulic fracturing.The new construction method was used to determine the field plan for fixed-length roof cutting.The roof formed a stable suspended structure and deformation of the main PLRR was improved after hydraulic fracturing.展开更多
Hydraulic fracture is important in unconventional oil and gas exploration.During the propagation of the hydraulic fracture,the crack tip is blunted due to the development of the process zone in the near-tip area.In th...Hydraulic fracture is important in unconventional oil and gas exploration.During the propagation of the hydraulic fracture,the crack tip is blunted due to the development of the process zone in the near-tip area.In this study,the blunting of the hydraulic fracture in polymethyl methacrylate specimens due to multi-timescale stress concentration is investigated.The ratio of the initiation toughness to the arrest toughness of the blunted hydraulic fracture is measured using both the dynamic and the static methods.Results show that a hydraulic fracture can be blunted with the time span of stress concentration from 1 ms to 600 s.It is also shown that the blunting of hydraulic fracture is a highly localized process.The morphology of the blunted crack depends on the stress distribution in the vicinity of the crack tip.展开更多
This study applies the Lindenmayer system based on fractal theory to generate synthetic fracture networks in hydraulically fractured wells.The applied flow model is based on complex analysis methods,which can quantify...This study applies the Lindenmayer system based on fractal theory to generate synthetic fracture networks in hydraulically fractured wells.The applied flow model is based on complex analysis methods,which can quantify the flow near the fractures,and being gridless,is computationally faster than traditional discrete volume simulations.The representation of hydraulic fractures as fractals is a more realistic representation than planar bi-wing fractures used in most reservoir models.Fluid withdrawal from the reservoir with evenly spaced hydraulic fractures may leave dead zones between planar fractures.Complex fractal networks will drain the reservoir matrix more effectively,due to the mitigation of stagnation flow zones.The flow velocities,pressure response,and drained rock volume(DRV)are visualized for a variety of fractal fracture networks in a single-fracture treatment stage.The major advancement of this study is the improved representation of hydraulic fractures as complex fractals rather than restricting to planar fracture geometries.Our models indicate that when the complexity of hydraulic fracture networks increases,this will suppress the occurrence of dead flow zones.In order to increase the DRV and improve ultimate recovery,our flow models suggest that fracture treatment programs must find ways to create more complex fracture networks.展开更多
Hydraulic fracture(HF) propagation behavior is significant when building enhanced geothermal systems(EGS). HF geometry is closely related to the structural planes(SPs) in hot dry rock(HDR), such as natural fractures(N...Hydraulic fracture(HF) propagation behavior is significant when building enhanced geothermal systems(EGS). HF geometry is closely related to the structural planes(SPs) in hot dry rock(HDR), such as natural fractures(NFs), quartz veins(QVs) and lithologic interfaces(LIs). However, the HF behaviors in HDR have not been well understood, especially the influence of multiple SPs on the HF geometry. To clarify this mechanism, several groups of physical simulation experiments of hydraulic fracturing were conducted to investigate the intersection relationship between the HFs and the SPs. Results show that the HF geometry shows great differences when intersecting with different SPs. In summary, the HF geometry displays four basic patterns, namely, propagation along the SPs, branching, capture, penetration/non-dilation. The fluctuation degree of the pressure-time curve and the HF complexity show a positive correlation. The cementing strength of the SP and their different mechanical properties from rock matrix influence the HF behaviors significantly. Therefore, the HF shows diverse geometries when intersecting with the NFs and LIs, while propagating along the QV when intersecting with it. For the complex networks, it is favorable for the HF to penetrate through and dilate several SPs, rather than simply cross or propagate along the SP.展开更多
Vertical height growth of hydraulic fractures(HFs)can unexpectedly penetrate a stratigraphic interface and propagate into neighboring layers,thereby resulting in low gas-production efficiency and high risk of groundwa...Vertical height growth of hydraulic fractures(HFs)can unexpectedly penetrate a stratigraphic interface and propagate into neighboring layers,thereby resulting in low gas-production efficiency and high risk of groundwater contamination or fault reactivation.Understanding of hydraulic fracture behavior at the interface is of pivotal importance for the successful development of layered reservoirs.In this paper,a twodimensional analytical model was developed to examine HF penetration and termination behavior at an orthogonal interface between two dissimilar materials.This model involves changes in the stress singularity ahead of the HF tip,which may alter at the formation interface due to material heterogeneity.Three critical stress conditions were considered to assess possible fracture behavior(i.e.,crossing,slippage,and opening)at the interface.Then,this model was verified by comparing its theoretical predictions to numerical simulations and three independent experiments.Good agreement with the simulation results and experimental data was observed,which shows the validity and reliability of this model.Finally,a parametric study was conducted to investigate the effects of key formation parameters(elastic modulus,Poisson’s ratio,and fracture toughness)between adjacent layers.These results indicate that the variation in the introduced parameters can limit or promote vertical HF growth by redistributing the induced normal and shear stresses at the interface.Among the three studied parameters,Poisson’s ratio has the least influence on the formation interface.When the fracture toughness and elastic modulus of the bounding layer are larger than those of the pay zone layer,the influence of fracture toughness will dominate the HF behavior at the interface;otherwise,the HF behavior will more likely be influenced by elastic modulus.展开更多
The fracture-cavity carbonate reservoirs in the Tahe Oilfield in China are mainly exploited by fracturing.We need the hydraulic fractures to communicate with caves to create a flow channel.However,due to the existence...The fracture-cavity carbonate reservoirs in the Tahe Oilfield in China are mainly exploited by fracturing.We need the hydraulic fractures to communicate with caves to create a flow channel.However,due to the existence of the fracture-cavity systems,the hydraulic fracture propagation morphology is complicated,while the propagation characteristics are not clear.To analyze the hydraulic fracture propagation in fracture-cavity carbonate formations,based on the discontinuous discrete fracture model,we developed a solid-seepage-freeflow coupled fracturing model for fracture-cavity formations,which can simulate the complex interaction behavior of fractures and caves.Based on the simulation results,we found the interaction rule between hydraulic fractures and fracture-cavity systems:the stress concentration around caves is the main factor that determines the fracture propagation path.Deflection due to stress concentration is usually not conducive to communication,while natural fractures distributed around caves could break the rejection action.Increasing the hydraulic energy in the hydraulic fracture can make fracture propagate directly and reduce the influence of deflection.The steering fracture formed by temporary plugging is beneficial to the communication of fracture-cavity systems in the non-principal stress direction.According to the simulation results of different fracture-cavity characteristics,we raised four optimization communication modes for fracture-cavity carbonate formation to provide references for fracturing optimization design and parameter optimization.展开更多
The existing acoustic logging methods for evaluating the hydraulic fracturing effectiveness usually use the fracture density to evaluate the fracture volume, and the results often cannot accurately reflect the actual ...The existing acoustic logging methods for evaluating the hydraulic fracturing effectiveness usually use the fracture density to evaluate the fracture volume, and the results often cannot accurately reflect the actual productivity. This paper studies the dynamic fluid flow through hydraulic fractures and its effect on borehole acoustic waves. Firstly, based on the fractal characteristics of fractures observed in hydraulic fracturing experiments, a permeability model of complex fracture network is established. Combining the dynamic fluid flow response of the model with the Biot-Rosenbaum theory that describes the acoustic wave propagation in permeable formations, the influence of hydraulic fractures on the velocity dispersion of borehole Stoneley-wave is then calculated and analyzed, whereby a novel hydraulic fracture fluid transport property evaluation method is proposed. The results show that the Stoneley-wave velocity dispersion characteristics caused by complex fractures can be equivalent to those of the plane fracture model, provided that the average permeability of the complex fracture model is equal to the permeability of the plane fracture. In addition, for fractures under high-permeability(fracture width 10~100 μm, permeability ~100 μm^(2)) and reduced permeability(1~10 μm, ~10 μm^(2), as in fracture closure) conditions, the Stoneley-wave velocity dispersion characteristics are significantly different. The field application shows that this fluid transport property evaluation method is practical to assess the permeability and the connectivity of hydraulic fractures.展开更多
In this study,we use the extended finite element method(XFEM)with a consideration of junction enrichment functions to investigate the mechanics of hydraulic fractures related to naturally cemented fractures.In the pro...In this study,we use the extended finite element method(XFEM)with a consideration of junction enrichment functions to investigate the mechanics of hydraulic fractures related to naturally cemented fractures.In the proposed numerical model,the lubrication equation is adopted to describe the fluid flow within fractures.The fluid-solid coupling systems of the hydraulic fracturing problem are solved using the Newton-Raphson method.The energy release rate criterion is used to determine the cross/arrest behavior between a hydraulic fracture(HF)and a cemented natural fracture(NF).The failure patterns and mechanisms of crack propagation at the intersection of natural fractures are discussed.Simulation results show that after crossing an NF,the failure mode along the cemented NF path may change from the tensile regime to the shear or mixed-mode regime.When an advancing HF kinks back toward the matrix,the failure mode may gradually switch back to the tensile-dominated regime.Key factors,including the length of the upper/lower portion of the cemented NF,horizontal stress anisotropy,and the intersection angle of the crack propagation are investigated in detail.An uncemented or partially cemented NF will form a more complex fracture network than a cemented NF.This study provides insight into the formation mechanism of fracture networks in formations that contain cemented NF.展开更多
The distribution of proppant injected in hydraulic fractures significantly affects the fracture conductivity and well performance.The proppant transport in thin fracturing fluid used during hydraulic fracturing in the...The distribution of proppant injected in hydraulic fractures significantly affects the fracture conductivity and well performance.The proppant transport in thin fracturing fluid used during hydraulic fracturing in the unconventional reservoirs is considerably different from fracturing fluids in the conventional reservoir due to the very low viscosity and quick deposition of the proppants.This paper presents the development of a three-dimensional Computational Fluid Dynamics(CFD)modelling technique for the prediction of proppant-fluid multiphase flow in hydraulic fractures.The proposed model also simulates the fluid leak-off behaviour from the fracture wall.The Euler-Granular and CFD-Discrete Element Method(CFD-DEM)multiphase modelling approach has been applied,and the equations defining the fluid-proppant and inter-proppant interaction have been solved using the finite volume technique.The proppant transport in hydraulic fractures has been studied comprehensively,and the computational modelling results of proppant distribution and other flow properties are in good agreement with the published experimental study.The parametric study is performed to investigate the effect of variation in proppant size,fluid viscosity and fracture width on the proppant transport.Smaller proppants can be injected early,followed by larger proppants to maintain high propping efficiency.This study has enhanced the understanding of the complex flow phenomenon between proppant and fracturing fluid and can play a vital role in hydraulic fracturing design.展开更多
During the stimulating unconventional reservoirs, the vertical propagation of hydraulic fractures is crucial for enlarging the stimulated reservoir volume, especially in multi-layers of sandstone, mudstone and shale(s...During the stimulating unconventional reservoirs, the vertical propagation of hydraulic fractures is crucial for enlarging the stimulated reservoir volume, especially in multi-layers of sandstone, mudstone and shale(sand-mud-shale). To investigate the effects of lithological interface and fracturing fluid viscosity on the fracture propagation vertically in the multi-layers, hydraulic fracturing experiments in laboratory were performed on the outcrop samples of 30 cm × 30 cm × 30 cm collected from Yanchang Formation in Ordos Basin. The results show that hydraulic fractures are multi-branched and zig-zagged when they initiate in shale, simple when they commence in sandstone or mudstone. Hydraulic fractures created with low-viscosity fracturing fluid can only cross sandstone from mudstone, but those induced by high-viscosity fracturing fluid can cross the sand-mud-shale layers. Furthermore, the high-viscosity fracturing fluid reduces the fractures complexity in shale, facilitating vertical fracture propagation. The injection pressure fluctuates slightly as the hydraulic fracture extends from shale to sandstone or mudstone, otherwise it fluctuates significantly. From the laboratory investigation, a hydraulic fracturing scheme for Chang 7 Member was proposed, with its feasibility proved in field tests.展开更多
Small-scale true triaxial sand fracturing experiments are conducted on thin interbedded shale samples made from cores of Permian Lucaogou Formation shale oil reservoir in Jimsar sag, Junggar Basin, NW China. Combined ...Small-scale true triaxial sand fracturing experiments are conducted on thin interbedded shale samples made from cores of Permian Lucaogou Formation shale oil reservoir in Jimsar sag, Junggar Basin, NW China. Combined with high-precision CT scanning digital core model reconstruction technology, hydraulic fracture geometry and proppant distribution in thin interbedded shale oil reservoirs are studied. The research shows that: In thin interbedded shale oil reservoir, the interlayer difference of rock mechanics and the interlayer interface near the wellbore cannot restrain the growth of fracture height effectively, but has a significant impact on the fracture width distribution in the fracture height direction. Hydraulic fractures in these reservoirs tend to penetrate into the adjacent layer in “step-like” form, but have a smaller width at the interface deflection, which hinders the transport of proppant in vertical direction, resulting in a poor effect of layer-crossing growth. In shale layers with dense laminae, hydraulic fractures tend to form “丰” or “井” shapes. If the perforated interval is large in rock strength and high in breakdown pressure, the main fracture is fully developed initially, large in width, and supported by enough sand. In contrast, if the perforated interval is low in strength and rich in laminae, the fracturing fluid filtration loss is large, the breakdown pressure is low, the main fracture will not open wide initially, and likely to have sand plugging. Proppant is mainly concentrated in the main hydraulic fractures with large width near the perforated layer, activated laminae, branch fractures and fractures in adjacent layers contain only a small amount of(or zero) proppant. The proppant is placed in a limited range on the whole. The limit width of fracture that proppant can enter is about 2.7 times the proppant particle size.展开更多
Field data suggests that carbonate reservoirs contain abundant natural fractures and cavities.The propagation mechanisms of hydraulic fractures in fracture-cavity reservoirs are different from conventional reservoirs ...Field data suggests that carbonate reservoirs contain abundant natural fractures and cavities.The propagation mechanisms of hydraulic fractures in fracture-cavity reservoirs are different from conventional reservoirs on account of the stress concentration surrounding cavities.In this paper,we develop a fully coupled numerical model using the extended finite element method(XFEM)to investigate the behaviors and propagation mechanisms of hydraulic fractures in fracture-cavity reservoirs.Simulation results show that a higher lateral stress coefficient can enhance the influence of the natural cavity,causing a more curved fracture path.However,lower confining stress or smaller in-situ stress difference can reduce this influence,and thus contributes to the penetration of the hydraulic fracture towards the cavity.Higher fluid viscosity and high fluid pumping rate are both able to attenuate the effect of the cavity.The frictional natural fracture connected to the cavity can significantly change the stress distribution around the cavity,thus dramatically deviates the hydraulic fracture from its original propagation direction.It is also found that the natural cavity existing between two adjacent fracturing stages will significantly influence the stress distribution between fractures and is more likely to result in irregular propagation paths compared to the case without a cavity.展开更多
This study extends an integrated field characterization in Eagle Ford by optimizing the numerical reservoir simulation of highly representative complex fractured systems through embedded discrete fracture modeling(EDF...This study extends an integrated field characterization in Eagle Ford by optimizing the numerical reservoir simulation of highly representative complex fractured systems through embedded discrete fracture modeling(EDFM). The bottom-hole flowing pressure was history-matched and the field production was forecasted after screening complex fracture scenarios with more than 100 000 fracture planes based on their propped-type. This work provided a greater understanding of the impact of complex-fractures proppant efficiency on the production. After compaction tables were included for each propped-type fracture group, the estimated pressure depletion showed that the effective drainage area can be smaller than the complex fracture network if modeled and screened by the EDFM method rather than unstructured gridding technique. The essential novel value of this work is the capability to couple EDFM with third-party fracture propagation simulation automatically, considering proppant intensity variation along the complex fractured systems. Thus, this work is pioneer to model complex fracture propagation and well interference accurately from fracture diagnostics and pseudo 3 D fracture propagation outcomes for multiple full wellbores to capture well completion effectiveness after myriads of sharper field simulation cases with EDFM.展开更多
The main area of the Jiaoshiba anticline of the Fuling shale gas field was taken as the research object,laboratory rock mechanical experiments and direct shear experiments were conducted to clarify the mechanical anis...The main area of the Jiaoshiba anticline of the Fuling shale gas field was taken as the research object,laboratory rock mechanical experiments and direct shear experiments were conducted to clarify the mechanical anisotropy characteristics and parameters of rock samples with rich beddings.Based on the experimental results,a 3D fracture propagation model of the target reservoir taking mechanical anisotropy,weak bedding plane and vertical stress difference into account was established by the discrete element method to analyze distribution patterns of hydraulic fractures under different bedding densities,mechanical properties,and fracturing engineering parameters(including perforation clusters,injection rates and fracturing fluid viscosity).The research results show that considering the influence of the weak bedding plane and longitudinal stress difference,the interlayer stress difference 3–4 MPa in the study area can control the fracture height within the zone of stress barrier,and the fracture height is less than 40 m.If the influence of the weak bedding plane is not considered,the simulation result of fracture height is obviously higher.Although the opening of high-density bedding fractures increases the complexity of hydraulic fractures,it significantly limited the propagation of fracture height.By reducing the number of clusters,increasing the injection rate,and increasing the volume and proportion of high-viscosity fracturing fluid in the pad stage,the restriction on fracture height due to the bedding plane and vertical stress difference can be reduced,and the longitudinal propagation of fractures can be promoted.The fracture propagation model was used to simulate one stage of Well A in Fuling shale gas field,and the simulation results were consistent with the micro-seismic monitoring results.展开更多
This paper describes numerical simulation of hydraulic fracturing using fracture-based continuum modeling(FBCM)of coupled geomechanical-hydrological processes to evaluate a technique for high-density fracturing and fr...This paper describes numerical simulation of hydraulic fracturing using fracture-based continuum modeling(FBCM)of coupled geomechanical-hydrological processes to evaluate a technique for high-density fracturing and fracture caging.The simulations are innovative because of modeling discrete fractures explicitly in continuum analysis.A key advantage of FBCM is that fracture initiation and propagation are modeled explicitly without changing the domain grid(i.e.no re-meshing).Further,multiple realizations of a preexisting fracture distribution can be analyzed using the same domain grid.The simulated hydraulic fracturing technique consists of pressurizing multiple wells simultaneously:initially without permeating fluids into the rock,to seed fractures uniformly and at high density in the wall rock of the wells;followed by fluid injection to propagate the seeded fracture density hydraulically.FBCM combines the ease of continuum modeling with the potential accuracy of modeling discrete fractures and fracturing explicitly.Fractures are modeled as piecewise planar based on intersections with domain elements;fracture geometry stored as continuum properties is used to calculate parameters needed to model individual fractures;and rock behavior is modeled through tensorial aggregation of the behavior of discrete fractures and unfractured rock.Simulations are presented for previously unfractured rock and for rock with preexisting fractures of horizontal,shallow-dipping,steeply dipping,or vertical orientation.Simulations of a single-well model are used to determine the pattern and spacing for a multiple-well design.The results illustrate high-density fracturing and fracture caging through simultaneous fluid injection in multiple wells:for previously unfractured rock or rock with preexisting shallow-dipping or horizontal fractures,and in situ vertical compressive stress greater than horizontal.If preexisting fractures are steeply dipping or vertical,and considering the same in situ stress condition,well pressurization without fluid permeation appears to be the only practical way to induce new fractures and contain fracturing within the target domain.展开更多
The hydraulic testing of pre-existing fractures(HTPF)is one of the most promising in situ stress measurement methods,particularly for three-dimensional stress tensor determination.However,the stress tensor determinati...The hydraulic testing of pre-existing fractures(HTPF)is one of the most promising in situ stress measurement methods,particularly for three-dimensional stress tensor determination.However,the stress tensor determination based on the HTPF method requires at least six tests or a minimum of 14-15 tests(under different conditions)for reliable results.In this study,we modified the HTPF method by considering the shear stress on each pre-existing fracture,which increased the number of equations for the stress tensor determination and decreased the number of tests required.Different shear stresses were attributed to different fractures by random sampling;therefore,the stress tensors were obtained by searching for the optimal solution using the least squares criterion based on the Monte Carlo method.Thereafter,we constrained the stress tensor based on the tensile strength criterion,compressive strength criterion,and vertical stress constraints.The inverted stress tensors were presented and analyzed based on the tensorial nature of the stress using the Euclidean mean stress tensor.Two stress-measurement campaigns in Weifang(Shandong Province,China)and Mercantour road tunnel(France)were implemented to highlight the validity and efficiency of the modified HTPF(M-HTPF)method.The results showed that the M-HTPF method can be applied for stress tensor inversion using only three to four tests on pre-existing fractures,neglecting the stress gradient.The inversion results were confined to relatively small distribution dispersions and were significantly reliable and stable due to the shear stresses on the fractures and the stress constraints employed.The M-HTPF method is highly feasible and efficient for complete stress tensor determination in a single borehole.展开更多
基金Supported by National Natural Science Foundation of China (U21A20105)Science and Technology Innovation Fund of PetroChina (2020D-5007-0208)。
文摘Using the visualized experimental device of temporary plugging in hydraulic fractures, the plugging behaviors of temporary plugging particles with different sizes and concentrations in hydraulic fractures were experimentally analyzed under the conditions of different carrier fluid displacements and viscosities. The results show that the greater the carrier fluid viscosity and displacement, the more difficult it is to form a plugging layer, and that the larger the size and concentration of the temporary plugging particle, the less difficult it is to form a plugging layer. When the ratio of particle size to fracture width is 0.45, the formation of the plugging layer is mainly controlled by the mass concentration of the temporary plugging particle and the viscosity of the carrier fluid, and a stable plugging layer cannot form if the mass concentration of the temporary plugging particle is less than 20 kg/m^(3)or the viscosity of the carrier fluid is greater than 3 mPa·s. When the ratio of particle size to fracture width is 0.60, the formation of the plugging layer is mainly controlled by the mass concentration of the temporary plugging particle, and a stable plugging layer cannot form if the mass concentration of the temporary plugging particle is less than 10 kg/m^(3). When the ratio of particle size to fracture width is 0.75, the formation of the plugging layer is basically not affected by other parameters, and a stable plugging layer can form within the experimental conditions. The formation process of plugging layer includes two stages and four modes. The main controlling factors affecting the formation mode are the ratio of particle size to fracture width, carrier fluid displacement and carrier fluid viscosity.
基金support from Major Program of the National Natural Science Foundation of China(Grant No.52192621)the National Natural Science Foundation of China for Major International(Regional)Joint Research Project(Grant No.52020105001)+1 种基金Major Science and Technology Project of Yunnan Province(Grant No.202302AF080001)Beijing Outstanding Young Scientist Program(Grant No.BJJWZYJH01201911414038).
文摘Hydraulic fracturing is considered the main stimulation method to develop shale gas reservoirs. Due to its strong heterogeneity, the shale gas formation is typically embedded with geological discontinuities such as bedding planes and natural fractures. Many researchers realized that the interaction between natural fractures and hydraulic fractures plays a crucial role in generating a complex fracture network. In this paper, true tri-axial hydraulic fracturing tests were performed on polymethyl methacrylate (PMMA), on which pre-existing fracture was pre-manufactured to simulate natural fracture. A cohesive model has been developed to verify the results of the experimental tests. The key findings demonstrate that the experimental results agreed well with the numerical simulation outcomes where three main interaction modes were observed: crossing;being arrested by opening the pre-existing fracture;being arrested without dilating the pre-existing fracture. Crossing behavior is more likely to occur with the approaching angle, horizontal stress difference, and injection rate increase. Furthermore, the higher flow rate might assist in reactivating the natural fractures where both sides of the pre-existing fractures were reactivated as the injection rate increased from 5 to 20 mL/min. Additionally, hydraulic fractures show a tendency to extend vertically rather than along the direction of maximum horizontal stress when they are first terminated at the interface. This research may contribute to the field application of hydraulic fracturing in shale gas formation.
基金supported financially by the Beijing Natural Science Foundation Project(No.3222030)the National Natural Science Foundation of China(No.51936001,No.52274002 and No.52192622)+1 种基金the PetroChina Science and Technology Innovation Foundation Project(2021DQ02–0201)Award Cultivation Foundation from Beijing Institute of Petrochemical Technology(No.BIPTACF-002).
文摘The effective plugging of artificial fractures is key to the success of temporary plugging and diverting fracturing technology,which is one of the most promising ways to improve the heat recovery efficiency of hot dry rock.At present,how temporary plugging agents plug artificial fractures under high temperature remains unclear.In this paper,by establishing an improved experimental system for the evaluation of temporary plugging performance at high temperature,we clarified the effects of high temperature,injection rate,and fracture width on the pressure response and plugging efficiency of the fracture.The results revealed that the temporary plugging process of artificial fractures in hot dry rock can be divided into four main stages:the initial stage of temporary plugging,the bridging stage of the particles,the plugging formation stage,and the high-pressure dense plugging stage.As the temperature increases,the distribution distance of the temporary plugging agent,the number of pressure fluctuations,and the time required for crack plugging increases.Particularly,when the temperature increases by 100℃,the complete plugging time increases by 90.7%.
基金the support provided by the Scientific Research and Technology Development Project of CNPC(Grant No.kt2017-19-01-1)the National Natural Science Foundation of China(Grant No.41772286,No.42077247 and No.42002271)+2 种基金Petro China Innovation Foundation(Grant No.2018D-5007-0202)Project funded by China Postdoctoral Science Foundation(Grant No.2021T140514)Fund of State Key Laboratory of Geomechanics and Geotechnical Engineering,Institute of Rock and Soil Mechanics,Chinese Academy of Sciences(Grant No.Z020009)。
文摘Depletion-induced stress change causes the redistribution of stress field in reservoirs,which can lead to the reorientation of principal stresses.Stress reorientation has a direct impact on fracture propagation of infill wells.To understand the effect of stress reorientation on the propagation of infill well’s fractures,an integrated simulation workflow that combines the reservoir flow calculation and the infill well hydraulic fracturing modeling is adopted.The reservoir simulation is computed to examine the relationship between the extent of stress reversal region and reservoir properties.Then,the hydraulic fracturing model considering the altered stress field for production is built to characterize the stress evolution of secondary fracturing.Numerical simulations show that stress reorientation may occur due to the decreasing of the horizontal stresses in an elliptical region around the parent well.Also,the initial stress difference is the driving factor for stress reorientation.However,the bottom hole pressure,permeability and other properties connected with fluid flow control timing of the stress reorientation.The decrease of the horizontal stresses around the parent well lead to asymmetrical propagation of a hydraulic fracture of the infill well.The study provides insights on understanding the influence of stress reorientation to the infill well fracturing treatment and interference between parent and infill wells.
基金financially supported by the Postgraduate Research&Practice Innovation Program of Jiangsu Province(KYCX21_2358)the National Key Research and Development Program of China(2020YFB1314204)National Natural Science Foundation of China(No.52074239)。
文摘Pre-driven longwall retracement roadway(PLRR)is commonly used in large mine shaft.The support crushing disasters occur frequently during the retracement,and roof management is necessary.Taking the 31107 panel as research background,the roof breaking structure of PLRR is analyzed.It is concluded that the roof cutting with vertical hydraulic fracture(HF)at a specified position,that is,fixed-length roof cutting,can reduce support load and keep immediate roof intact.The extended finite element method(XFEM)is applied to simulate hydraulic fracturing.The results show that both the axial and transverse hydraulic fracturing cannot effectively create vertical HFs.Therefore,a novel construction method of vertical HF based on the stress shadow effect(SSE)is proposed.The stress reversal region and HF orientation caused by the prefabricated hydraulic fracture(PF)are verified in simulation.The sub-vertical HFs are obtained between two PFs,the vertical extension range of which is much larger than that of directional hydraulic fracturing.The new construction method was used to determine the field plan for fixed-length roof cutting.The roof formed a stable suspended structure and deformation of the main PLRR was improved after hydraulic fracturing.
基金the support from the China National Science and Technology Major Project"Changning-Weiyuan shale gas development demonstration project"(2016ZX05062)the support from the China National Petroleum Corporation:"Research on Influencing Factors of Gas Hydrate Sand Production and Experimental Design"(No.CPETQ201921)。
文摘Hydraulic fracture is important in unconventional oil and gas exploration.During the propagation of the hydraulic fracture,the crack tip is blunted due to the development of the process zone in the near-tip area.In this study,the blunting of the hydraulic fracture in polymethyl methacrylate specimens due to multi-timescale stress concentration is investigated.The ratio of the initiation toughness to the arrest toughness of the blunted hydraulic fracture is measured using both the dynamic and the static methods.Results show that a hydraulic fracture can be blunted with the time span of stress concentration from 1 ms to 600 s.It is also shown that the blunting of hydraulic fracture is a highly localized process.The morphology of the blunted crack depends on the stress distribution in the vicinity of the crack tip.
文摘This study applies the Lindenmayer system based on fractal theory to generate synthetic fracture networks in hydraulically fractured wells.The applied flow model is based on complex analysis methods,which can quantify the flow near the fractures,and being gridless,is computationally faster than traditional discrete volume simulations.The representation of hydraulic fractures as fractals is a more realistic representation than planar bi-wing fractures used in most reservoir models.Fluid withdrawal from the reservoir with evenly spaced hydraulic fractures may leave dead zones between planar fractures.Complex fractal networks will drain the reservoir matrix more effectively,due to the mitigation of stagnation flow zones.The flow velocities,pressure response,and drained rock volume(DRV)are visualized for a variety of fractal fracture networks in a single-fracture treatment stage.The major advancement of this study is the improved representation of hydraulic fractures as complex fractals rather than restricting to planar fracture geometries.Our models indicate that when the complexity of hydraulic fracture networks increases,this will suppress the occurrence of dead flow zones.In order to increase the DRV and improve ultimate recovery,our flow models suggest that fracture treatment programs must find ways to create more complex fracture networks.
基金support from the National Key Research and Development Project in China(Grant nos.2018YFB1501801,2020YFE020130005)the National Natural Science Foundation of China(Grant nos.42004115,42102353)the Open Research Fund Program of Key Laboratory of Metallogenic Prediction of Nonferrous Metals and Geological Environment Monitoring(Central South University),Ministry of Education(Grant No.2021YSJS08)。
文摘Hydraulic fracture(HF) propagation behavior is significant when building enhanced geothermal systems(EGS). HF geometry is closely related to the structural planes(SPs) in hot dry rock(HDR), such as natural fractures(NFs), quartz veins(QVs) and lithologic interfaces(LIs). However, the HF behaviors in HDR have not been well understood, especially the influence of multiple SPs on the HF geometry. To clarify this mechanism, several groups of physical simulation experiments of hydraulic fracturing were conducted to investigate the intersection relationship between the HFs and the SPs. Results show that the HF geometry shows great differences when intersecting with different SPs. In summary, the HF geometry displays four basic patterns, namely, propagation along the SPs, branching, capture, penetration/non-dilation. The fluctuation degree of the pressure-time curve and the HF complexity show a positive correlation. The cementing strength of the SP and their different mechanical properties from rock matrix influence the HF behaviors significantly. Therefore, the HF shows diverse geometries when intersecting with the NFs and LIs, while propagating along the QV when intersecting with it. For the complex networks, it is favorable for the HF to penetrate through and dilate several SPs, rather than simply cross or propagate along the SP.
基金supported by the National Natural Science Foundation of China(No.52064006,52164001 and 52004072)the Guizhou Provincial Science and Technology Foundation(No.[2020]4Y044,No.[2021]292,No.GCC[2022]005 and[2021]N404)the China Scholarship Council program(202006050112)
文摘Vertical height growth of hydraulic fractures(HFs)can unexpectedly penetrate a stratigraphic interface and propagate into neighboring layers,thereby resulting in low gas-production efficiency and high risk of groundwater contamination or fault reactivation.Understanding of hydraulic fracture behavior at the interface is of pivotal importance for the successful development of layered reservoirs.In this paper,a twodimensional analytical model was developed to examine HF penetration and termination behavior at an orthogonal interface between two dissimilar materials.This model involves changes in the stress singularity ahead of the HF tip,which may alter at the formation interface due to material heterogeneity.Three critical stress conditions were considered to assess possible fracture behavior(i.e.,crossing,slippage,and opening)at the interface.Then,this model was verified by comparing its theoretical predictions to numerical simulations and three independent experiments.Good agreement with the simulation results and experimental data was observed,which shows the validity and reliability of this model.Finally,a parametric study was conducted to investigate the effects of key formation parameters(elastic modulus,Poisson’s ratio,and fracture toughness)between adjacent layers.These results indicate that the variation in the introduced parameters can limit or promote vertical HF growth by redistributing the induced normal and shear stresses at the interface.Among the three studied parameters,Poisson’s ratio has the least influence on the formation interface.When the fracture toughness and elastic modulus of the bounding layer are larger than those of the pay zone layer,the influence of fracture toughness will dominate the HF behavior at the interface;otherwise,the HF behavior will more likely be influenced by elastic modulus.
基金the National Natural Science Foundation Program(No.51874321)。
文摘The fracture-cavity carbonate reservoirs in the Tahe Oilfield in China are mainly exploited by fracturing.We need the hydraulic fractures to communicate with caves to create a flow channel.However,due to the existence of the fracture-cavity systems,the hydraulic fracture propagation morphology is complicated,while the propagation characteristics are not clear.To analyze the hydraulic fracture propagation in fracture-cavity carbonate formations,based on the discontinuous discrete fracture model,we developed a solid-seepage-freeflow coupled fracturing model for fracture-cavity formations,which can simulate the complex interaction behavior of fractures and caves.Based on the simulation results,we found the interaction rule between hydraulic fractures and fracture-cavity systems:the stress concentration around caves is the main factor that determines the fracture propagation path.Deflection due to stress concentration is usually not conducive to communication,while natural fractures distributed around caves could break the rejection action.Increasing the hydraulic energy in the hydraulic fracture can make fracture propagate directly and reduce the influence of deflection.The steering fracture formed by temporary plugging is beneficial to the communication of fracture-cavity systems in the non-principal stress direction.According to the simulation results of different fracture-cavity characteristics,we raised four optimization communication modes for fracture-cavity carbonate formation to provide references for fracturing optimization design and parameter optimization.
基金Supported by the National Natural Science Foundation of China (41821002,42174145)PetroChina Science and Technology Major Project (ZD2019-183-004)China University of Petroleum (East China) Graduate Student Innovation Project (YCX2019001)。
文摘The existing acoustic logging methods for evaluating the hydraulic fracturing effectiveness usually use the fracture density to evaluate the fracture volume, and the results often cannot accurately reflect the actual productivity. This paper studies the dynamic fluid flow through hydraulic fractures and its effect on borehole acoustic waves. Firstly, based on the fractal characteristics of fractures observed in hydraulic fracturing experiments, a permeability model of complex fracture network is established. Combining the dynamic fluid flow response of the model with the Biot-Rosenbaum theory that describes the acoustic wave propagation in permeable formations, the influence of hydraulic fractures on the velocity dispersion of borehole Stoneley-wave is then calculated and analyzed, whereby a novel hydraulic fracture fluid transport property evaluation method is proposed. The results show that the Stoneley-wave velocity dispersion characteristics caused by complex fractures can be equivalent to those of the plane fracture model, provided that the average permeability of the complex fracture model is equal to the permeability of the plane fracture. In addition, for fractures under high-permeability(fracture width 10~100 μm, permeability ~100 μm^(2)) and reduced permeability(1~10 μm, ~10 μm^(2), as in fracture closure) conditions, the Stoneley-wave velocity dispersion characteristics are significantly different. The field application shows that this fluid transport property evaluation method is practical to assess the permeability and the connectivity of hydraulic fractures.
基金financially supported by the National Science Foundation of China(Grant Nos.51804033 and 51936001)Natural Science Foundation of Jiangsu Province(Grant No.BK20170457)+3 种基金Program of Great Wall Scholar(Grant No.CIT&TCD20180313)Jointly Projects of Beijing Natural Science FoundationBeijing Municipal Education Commission(Grant No.KZ201810017023)Beijing Youth Talent Support Program(CIT&TCD201804037).
文摘In this study,we use the extended finite element method(XFEM)with a consideration of junction enrichment functions to investigate the mechanics of hydraulic fractures related to naturally cemented fractures.In the proposed numerical model,the lubrication equation is adopted to describe the fluid flow within fractures.The fluid-solid coupling systems of the hydraulic fracturing problem are solved using the Newton-Raphson method.The energy release rate criterion is used to determine the cross/arrest behavior between a hydraulic fracture(HF)and a cemented natural fracture(NF).The failure patterns and mechanisms of crack propagation at the intersection of natural fractures are discussed.Simulation results show that after crossing an NF,the failure mode along the cemented NF path may change from the tensile regime to the shear or mixed-mode regime.When an advancing HF kinks back toward the matrix,the failure mode may gradually switch back to the tensile-dominated regime.Key factors,including the length of the upper/lower portion of the cemented NF,horizontal stress anisotropy,and the intersection angle of the crack propagation are investigated in detail.An uncemented or partially cemented NF will form a more complex fracture network than a cemented NF.This study provides insight into the formation mechanism of fracture networks in formations that contain cemented NF.
文摘The distribution of proppant injected in hydraulic fractures significantly affects the fracture conductivity and well performance.The proppant transport in thin fracturing fluid used during hydraulic fracturing in the unconventional reservoirs is considerably different from fracturing fluids in the conventional reservoir due to the very low viscosity and quick deposition of the proppants.This paper presents the development of a three-dimensional Computational Fluid Dynamics(CFD)modelling technique for the prediction of proppant-fluid multiphase flow in hydraulic fractures.The proposed model also simulates the fluid leak-off behaviour from the fracture wall.The Euler-Granular and CFD-Discrete Element Method(CFD-DEM)multiphase modelling approach has been applied,and the equations defining the fluid-proppant and inter-proppant interaction have been solved using the finite volume technique.The proppant transport in hydraulic fractures has been studied comprehensively,and the computational modelling results of proppant distribution and other flow properties are in good agreement with the published experimental study.The parametric study is performed to investigate the effect of variation in proppant size,fluid viscosity and fracture width on the proppant transport.Smaller proppants can be injected early,followed by larger proppants to maintain high propping efficiency.This study has enhanced the understanding of the complex flow phenomenon between proppant and fracturing fluid and can play a vital role in hydraulic fracturing design.
基金sponsored by the Strategic Cooperation Technology Projects of CNPC and CUPB (ZLZX2020-02)the National Science Fund for Distinguished Young Scholars (Grant No.51925405)the National Natural Science Foundation of China(Grant no. 51774299)。
文摘During the stimulating unconventional reservoirs, the vertical propagation of hydraulic fractures is crucial for enlarging the stimulated reservoir volume, especially in multi-layers of sandstone, mudstone and shale(sand-mud-shale). To investigate the effects of lithological interface and fracturing fluid viscosity on the fracture propagation vertically in the multi-layers, hydraulic fracturing experiments in laboratory were performed on the outcrop samples of 30 cm × 30 cm × 30 cm collected from Yanchang Formation in Ordos Basin. The results show that hydraulic fractures are multi-branched and zig-zagged when they initiate in shale, simple when they commence in sandstone or mudstone. Hydraulic fractures created with low-viscosity fracturing fluid can only cross sandstone from mudstone, but those induced by high-viscosity fracturing fluid can cross the sand-mud-shale layers. Furthermore, the high-viscosity fracturing fluid reduces the fractures complexity in shale, facilitating vertical fracture propagation. The injection pressure fluctuates slightly as the hydraulic fracture extends from shale to sandstone or mudstone, otherwise it fluctuates significantly. From the laboratory investigation, a hydraulic fracturing scheme for Chang 7 Member was proposed, with its feasibility proved in field tests.
基金National Natural Science Foundation of China(NO.51974332)Strategic Cooperation Project Between PetroChina and China University of Petroleum(Beijing)(NO.ZLZX2020-07)。
文摘Small-scale true triaxial sand fracturing experiments are conducted on thin interbedded shale samples made from cores of Permian Lucaogou Formation shale oil reservoir in Jimsar sag, Junggar Basin, NW China. Combined with high-precision CT scanning digital core model reconstruction technology, hydraulic fracture geometry and proppant distribution in thin interbedded shale oil reservoirs are studied. The research shows that: In thin interbedded shale oil reservoir, the interlayer difference of rock mechanics and the interlayer interface near the wellbore cannot restrain the growth of fracture height effectively, but has a significant impact on the fracture width distribution in the fracture height direction. Hydraulic fractures in these reservoirs tend to penetrate into the adjacent layer in “step-like” form, but have a smaller width at the interface deflection, which hinders the transport of proppant in vertical direction, resulting in a poor effect of layer-crossing growth. In shale layers with dense laminae, hydraulic fractures tend to form “丰” or “井” shapes. If the perforated interval is large in rock strength and high in breakdown pressure, the main fracture is fully developed initially, large in width, and supported by enough sand. In contrast, if the perforated interval is low in strength and rich in laminae, the fracturing fluid filtration loss is large, the breakdown pressure is low, the main fracture will not open wide initially, and likely to have sand plugging. Proppant is mainly concentrated in the main hydraulic fractures with large width near the perforated layer, activated laminae, branch fractures and fractures in adjacent layers contain only a small amount of(or zero) proppant. The proppant is placed in a limited range on the whole. The limit width of fracture that proppant can enter is about 2.7 times the proppant particle size.
基金This research was jointly funded by the National Natural Science Foundation of China(No.51904111)the Natural Science Foundation of Jiangsu Province(No.BK20170457)+1 种基金the Open Fund for Jiangsu Key Laboratory of Advanced Manufacturing Technology(No.HGAMTL-1712)the Natural Science Research of Institution of Higher Education of Jiangsu Province(No.17KJA460003).
文摘Field data suggests that carbonate reservoirs contain abundant natural fractures and cavities.The propagation mechanisms of hydraulic fractures in fracture-cavity reservoirs are different from conventional reservoirs on account of the stress concentration surrounding cavities.In this paper,we develop a fully coupled numerical model using the extended finite element method(XFEM)to investigate the behaviors and propagation mechanisms of hydraulic fractures in fracture-cavity reservoirs.Simulation results show that a higher lateral stress coefficient can enhance the influence of the natural cavity,causing a more curved fracture path.However,lower confining stress or smaller in-situ stress difference can reduce this influence,and thus contributes to the penetration of the hydraulic fracture towards the cavity.Higher fluid viscosity and high fluid pumping rate are both able to attenuate the effect of the cavity.The frictional natural fracture connected to the cavity can significantly change the stress distribution around the cavity,thus dramatically deviates the hydraulic fracture from its original propagation direction.It is also found that the natural cavity existing between two adjacent fracturing stages will significantly influence the stress distribution between fractures and is more likely to result in irregular propagation paths compared to the case without a cavity.
文摘This study extends an integrated field characterization in Eagle Ford by optimizing the numerical reservoir simulation of highly representative complex fractured systems through embedded discrete fracture modeling(EDFM). The bottom-hole flowing pressure was history-matched and the field production was forecasted after screening complex fracture scenarios with more than 100 000 fracture planes based on their propped-type. This work provided a greater understanding of the impact of complex-fractures proppant efficiency on the production. After compaction tables were included for each propped-type fracture group, the estimated pressure depletion showed that the effective drainage area can be smaller than the complex fracture network if modeled and screened by the EDFM method rather than unstructured gridding technique. The essential novel value of this work is the capability to couple EDFM with third-party fracture propagation simulation automatically, considering proppant intensity variation along the complex fractured systems. Thus, this work is pioneer to model complex fracture propagation and well interference accurately from fracture diagnostics and pseudo 3 D fracture propagation outcomes for multiple full wellbores to capture well completion effectiveness after myriads of sharper field simulation cases with EDFM.
基金Supported by the China National Science and Technology Major Project(2016ZX05060001-032)
文摘The main area of the Jiaoshiba anticline of the Fuling shale gas field was taken as the research object,laboratory rock mechanical experiments and direct shear experiments were conducted to clarify the mechanical anisotropy characteristics and parameters of rock samples with rich beddings.Based on the experimental results,a 3D fracture propagation model of the target reservoir taking mechanical anisotropy,weak bedding plane and vertical stress difference into account was established by the discrete element method to analyze distribution patterns of hydraulic fractures under different bedding densities,mechanical properties,and fracturing engineering parameters(including perforation clusters,injection rates and fracturing fluid viscosity).The research results show that considering the influence of the weak bedding plane and longitudinal stress difference,the interlayer stress difference 3–4 MPa in the study area can control the fracture height within the zone of stress barrier,and the fracture height is less than 40 m.If the influence of the weak bedding plane is not considered,the simulation result of fracture height is obviously higher.Although the opening of high-density bedding fractures increases the complexity of hydraulic fractures,it significantly limited the propagation of fracture height.By reducing the number of clusters,increasing the injection rate,and increasing the volume and proportion of high-viscosity fracturing fluid in the pad stage,the restriction on fracture height due to the bedding plane and vertical stress difference can be reduced,and the longitudinal propagation of fractures can be promoted.The fracture propagation model was used to simulate one stage of Well A in Fuling shale gas field,and the simulation results were consistent with the micro-seismic monitoring results.
文摘This paper describes numerical simulation of hydraulic fracturing using fracture-based continuum modeling(FBCM)of coupled geomechanical-hydrological processes to evaluate a technique for high-density fracturing and fracture caging.The simulations are innovative because of modeling discrete fractures explicitly in continuum analysis.A key advantage of FBCM is that fracture initiation and propagation are modeled explicitly without changing the domain grid(i.e.no re-meshing).Further,multiple realizations of a preexisting fracture distribution can be analyzed using the same domain grid.The simulated hydraulic fracturing technique consists of pressurizing multiple wells simultaneously:initially without permeating fluids into the rock,to seed fractures uniformly and at high density in the wall rock of the wells;followed by fluid injection to propagate the seeded fracture density hydraulically.FBCM combines the ease of continuum modeling with the potential accuracy of modeling discrete fractures and fracturing explicitly.Fractures are modeled as piecewise planar based on intersections with domain elements;fracture geometry stored as continuum properties is used to calculate parameters needed to model individual fractures;and rock behavior is modeled through tensorial aggregation of the behavior of discrete fractures and unfractured rock.Simulations are presented for previously unfractured rock and for rock with preexisting fractures of horizontal,shallow-dipping,steeply dipping,or vertical orientation.Simulations of a single-well model are used to determine the pattern and spacing for a multiple-well design.The results illustrate high-density fracturing and fracture caging through simultaneous fluid injection in multiple wells:for previously unfractured rock or rock with preexisting shallow-dipping or horizontal fractures,and in situ vertical compressive stress greater than horizontal.If preexisting fractures are steeply dipping or vertical,and considering the same in situ stress condition,well pressurization without fluid permeation appears to be the only practical way to induce new fractures and contain fracturing within the target domain.
基金supported by the National Natural Science Foundation of China(Grant No.42174118)a research grant(Grant No.ZDJ 2020-7)from the National Institute of Natural Hazards,Ministry of Emergency Management of China.
文摘The hydraulic testing of pre-existing fractures(HTPF)is one of the most promising in situ stress measurement methods,particularly for three-dimensional stress tensor determination.However,the stress tensor determination based on the HTPF method requires at least six tests or a minimum of 14-15 tests(under different conditions)for reliable results.In this study,we modified the HTPF method by considering the shear stress on each pre-existing fracture,which increased the number of equations for the stress tensor determination and decreased the number of tests required.Different shear stresses were attributed to different fractures by random sampling;therefore,the stress tensors were obtained by searching for the optimal solution using the least squares criterion based on the Monte Carlo method.Thereafter,we constrained the stress tensor based on the tensile strength criterion,compressive strength criterion,and vertical stress constraints.The inverted stress tensors were presented and analyzed based on the tensorial nature of the stress using the Euclidean mean stress tensor.Two stress-measurement campaigns in Weifang(Shandong Province,China)and Mercantour road tunnel(France)were implemented to highlight the validity and efficiency of the modified HTPF(M-HTPF)method.The results showed that the M-HTPF method can be applied for stress tensor inversion using only three to four tests on pre-existing fractures,neglecting the stress gradient.The inversion results were confined to relatively small distribution dispersions and were significantly reliable and stable due to the shear stresses on the fractures and the stress constraints employed.The M-HTPF method is highly feasible and efficient for complete stress tensor determination in a single borehole.
