To investigate the height growth of multi-cluster fractures during variable fluid-viscosity fracturing in a layered shale oil reservoir,a two-dimensional finite element method(FEM)-discrete fracture network(DFN)model ...To investigate the height growth of multi-cluster fractures during variable fluid-viscosity fracturing in a layered shale oil reservoir,a two-dimensional finite element method(FEM)-discrete fracture network(DFN)model coupled with flow,stress and damage is proposed.A traction-separation law is used to describe the mixed-mode response of the damaged adhesive fractures,and the cubic law is used to describe the fluid flow within the fractures.The rock deformation is controlled by the in-situ stress,fracture cohesion and fluid pressure on the hydraulic fracture surface.The coupled finite element equations are solved by the explicit time difference method.The effects of the fracturing treatment parameters including fluid viscosity,pumping rate and cluster spacing on the geometries of multifractures are investigated.The results show that variable fluid-viscosity injection can improve the complexity of the fracture network and height of the main fractures simultaneously.The pumping rate of15 m^(3)/min,variable fluid-viscosity of 3-9-21-36-45 mPa s with a cluster spacing of 7.5 m is the ideal treatment strategy.The field application shows that the peak daily production of the application well with the optimized injection procedu re of variable fluid-viscosity fracturing is 171 tons(about 2.85 times that of the adjacent well),which is the highest daily production record of a single shale oil well in China,marking a strategic breakthrough of commercial shale oil production in the Jiyang Depression,Shengli Oilfield.The variable fluid-viscosity fracturing technique is proved to be very effective for improving shale oil production.展开更多
The continental shale reservoirs of Jurassic Lianggaoshan Formation in Sichuan Basin contain thin lamina,which is characterized by strong plasticity and developed longitudinal shell limestone interlayer.To improve the...The continental shale reservoirs of Jurassic Lianggaoshan Formation in Sichuan Basin contain thin lamina,which is characterized by strong plasticity and developed longitudinal shell limestone interlayer.To improve the production efficiency of reservoirs by multi-cluster fracturing,it is necessary to consider the unbalanced propagation of hydraulic fractures and the penetration effect of fractures.This paper constructed a numerical model of multi-fracture propagation and penetration based on the finite element coupling cohesive zone method;considering the construction cluster spacing,pump rate,lamina strength and other parameters studied the influencing factors of multi-cluster fracture interaction propagation;combined with AE energy data and fracture mode reconstruction method,quantitatively characterized the comprehensive impact of the strength of thin interlayer rock interfaces on the initiation and propagation of fractures that penetrate layers,and accurately predicted the propagation pattern of hydraulic fractures through laminated shale oil reservoirs.Simulation results revealed that in the process of multi-cluster fracturing,the proportion of shear damage is low,and mainly occurs in bedding fractures activated by outer fractures.Reducing the cluster spacing enhances the fracture system's penetration ability,though it lowers the activation efficiency of lamina.The high plasticity of the limestone interlayer may impact the vertical propagation distance of the main fracture.Improving the interface strength is beneficial to the reconstruction of the fracture height,but the interface communication effect is limited.Reasonable selection of layers with moderate lamina strength for fracturing stimulation,increasing the pump rate during fracturing and setting the cluster spacing reasonably are beneficial to improve the effect of reservoir stimulation.展开更多
This study conducted temporary plugging and diversion fracturing(TPDF)experiments using a true triaxial fracturing simulation system within a laboratory setting that replicated a lab-based horizontal well completion w...This study conducted temporary plugging and diversion fracturing(TPDF)experiments using a true triaxial fracturing simulation system within a laboratory setting that replicated a lab-based horizontal well completion with multi-cluster sand jetting perforation.The effects of temporary plugging agent(TPA)particle size,TPA concentration,single-cluster perforation number and cluster number on plugging pressure,multi-fracture diversion pattern and distribution of TPAs were investigated.A combination of TPAs with small particle sizes within the fracture and large particle sizes within the segment is conducive to increasing the plugging pressure and promoting the diversion of multi-fractures.The addition of fibers can quickly achieve ultra-high pressure,but it may lead to longitudinal fractures extending along the wellbore.The temporary plugging peak pressure increases with an increase in the concentration of the TPA,reaching a peak at a certain concentration,and further increases do not significantly improve the temporary plugging peak pressure.The breaking pressure and temporary plugging peak pressure show a decreasing trend with an increase in single-cluster perforation number.A lower number of single-cluster perforations is beneficial for increasing the breaking pressure and temporary plugging peak pressure,and it has a more significant control on the propagation of multi-cluster fractures.A lower number of clusters is not conducive to increasing the total number and complexity of artificial fractures,while a higher number of clusters makes it difficult to achieve effective plugging.The TPAs within the fracture is mainly concentrated in the complex fracture areas,especially at the intersections of fractures.Meanwhile,the TPAs within the segment are primarily distributed near the perforation cluster apertures which initiated complex fractures.展开更多
The exploitation of shale gas is promising due to depletion of the conventional energy and intensification of the greenhouse effect.In this paper,we proposed a heat-fluid-solid coupling damage model of supercritical C...The exploitation of shale gas is promising due to depletion of the conventional energy and intensification of the greenhouse effect.In this paper,we proposed a heat-fluid-solid coupling damage model of supercritical CO_(2)(SC-CO_(2))compound fracturing which is expected to be an efficient and environmentally friendly way to develop shale gas.The coupling model is solved by the finite element method,and the results are in good agreement with the analytical solutions and fracturing experiments.Based on this model,the fracture propagation characteristics at the two stages of compound fracturing are studied and the influence of pressurization rate,in situ stress,bedding angle,and other factors are considered.The results show that at the SC-CO_(2)fracturing stage,a lower pressurization rate is conducive to formation of the branches around main fractures,while a higher pressurization rate inhibits formation of the branches around main fractures and promotes formation of the main fractures.Both bedding and in situ stress play a dominant role in the fracture propagation.When the in situ stress ratio(δ_(x)/δ_(y))is 1,the presence of bedding can reduce the initiation pressure and failure pressure.Nevertheless,it will cause the fracture to propagate along the bedding direction,reducing the fracture complexity.In rocks without bedding,hydraulic fracturing has the lengthening and widening effects for SC-CO_(2)induced fracture.In shale,fractures induced at the hydraulic fracturing stage are more likely to be dominated by in situ stresses and have a shorter reorientation radius.Therefore,fracture branches propagating along the maximum principal stress direction may be generated around the main fractures induced by SC-CO_(2)at the hydraulic fracturing stage.When the branches converge with the main fractures,fracture zones are easily formed,and thus the fracture complexity and damage area can be significantly increased.The results are instructive for the design and application of SC-CO_(2)compound fracturing.展开更多
The flow of fluid through the porous matrix of a reservoir rock applies a seepage force to the solid rock matrix.Although the seepage force exerted by fluid flow through the porous matrix of a reservoir rock has a not...The flow of fluid through the porous matrix of a reservoir rock applies a seepage force to the solid rock matrix.Although the seepage force exerted by fluid flow through the porous matrix of a reservoir rock has a notable influence on rock deformation and failure,its effect on hydraulic fracture(HF)propagation remains ambiguous.Therefore,in this study,we improved a traditional fluid–solid coupling method by incorporating the role of seepage force during the fracturing fluid seepage,using the discrete element method.First,we validated the simulation results of the improved method by comparing them with an analytical solution of the seepage force and published experimental results.Next,we conducted numerical simulations in both homogeneous and heterogeneous sandstone formations to investigate the influence of seepage force on HF propagation.Our results indicate that fluid viscosity has a greater impact on the magnitude and extent of seepage force compared to injection rate,and that lower viscosity and injection rate correspond to shorter hydraulic fracture lengths.Furthermore,seepage force influences the direction of HF propagation,causing HFs to deflect towards the side of the reservoir with weaker cementation and higher permeability.展开更多
Deep shale reservoirs are characterized by elevated breakdown pressures,diminished fracture complexity,and reduced modified volumes compared to medium and shallow reservoirs.Therefore,it is urgent to investigate parti...Deep shale reservoirs are characterized by elevated breakdown pressures,diminished fracture complexity,and reduced modified volumes compared to medium and shallow reservoirs.Therefore,it is urgent to investigate particular injection strategies that can optimize breakdown pressure and fracturing efficiency to address the increasing demands for deep shale reservoir stimulation.In this study,the efficiency of various stimulation strategies,including multi-cluster simultaneous fracturing,modified alternating fracturing,alternating shut-in fracturing,and cyclic alternating fracturing,was evaluated.Subsequently,the sensitivity of factors such as the cycle index,shut-in time,cluster spacing,and horizontal permeability was investigated.Additionally,the flow distribution effect within the wellbore was discussed.The results indicate that relative to multi-cluster simultaneous fracturing,modified alternating fracturing exhibits reduced susceptibility to the stress shadow effect,which results in earlier breakdown,extended hydraulic fracture lengths,and more consistent propagation despite an increase in breakdown pressure.The alternating shut-in fracturing benefits the increase of fracture length,which is closely related to the shut-in time.Furthermore,cyclic alternating fracturing markedly lowers breakdown pressure and contributes to uniform fracture propagation,in which the cycle count plays an important role.Modified alternating fracturing demonstrates insensitivity to variations in cluster spacing,whereas horizontal permeability is a critical factor affecting fracture length.The wellbore effect restrains the accumulation of pressure and flow near the perforation,delaying the initiation of hydraulic fractures.The simulation results can provide valuable numerical insights for optimizing injection strategies for deep shale hydraulic fracturing.展开更多
Indirect fracturing in the roof of broken soft coal seams has been demonstrated to be a feasible technology.In this work,the No.5 coal seam in the Hancheng block was taken as the research object.Based on the findings ...Indirect fracturing in the roof of broken soft coal seams has been demonstrated to be a feasible technology.In this work,the No.5 coal seam in the Hancheng block was taken as the research object.Based on the findings of true triaxial hydraulic fracturing experiments and field pilot under this technology and the cohesive element method,a 3D numerical model of indirect fracturing in the roof of broken soft coal seams was established,the fracture morphology propagation and evolution law under different conditions was investigated,and analysis of main controlling factors of fracture parameters was conducted with the combination weight method,which was based on grey incidence,analytic hierarchy process and entropy weight method.The results show that“士”-shaped fractures,T-shaped fractures,cross fractures,H-shaped fractures,and“干”-shaped fractures dominated by horizontal fractures were formed.Different parameter combinations can form different fracture morphologies.When the coal seam permeability is lower and the minimum horizontal principal stress difference between layers and fracturing fluid injection rate are both larger,it tends to form“士”-shaped fractures.When the coal seam permeability and minimum horizontal principal stress between layers and perforation position are moderate,cross fractures are easily generated.Different fracture parameters have different main controlling factors.Engineering factors of perforation location,fracturing fluid injection rate and viscosity are the dominant factors of hydraulic fracture shape parameters.This study can provide a reference for the design of indirect fracturing in the roof of broken soft coal seams.展开更多
Multistage multi-cluster hydraulic fracturing has enabled the economic exploitation of shale reservoirs,but the interpretation of hydraulic fracture parameters is challenging.The pressure signals after pump shutdown a...Multistage multi-cluster hydraulic fracturing has enabled the economic exploitation of shale reservoirs,but the interpretation of hydraulic fracture parameters is challenging.The pressure signals after pump shutdown are influenced by hydraulic fractures,which can reflect the geometric features of hydraulic fracture.The shutdown pressure can be used to interpret the hydraulic fracture parameters in a real-time and cost-effective manner.In this paper,a mathematical model for shutdown pressure evolution is developed considering the effects of wellbore friction,perforation friction and fluid loss in fractures.An efficient numerical simulation method is established by using the method of characteristics.Based on this method,the impacts of fracture half-length,fracture height,opened cluster and perforation number,and filtration coefficient on the evolution of shutdown pressure are analyzed.The results indicate that a larger fracture half-length may hasten the decay of shutdown pressure,while a larger fracture height can slow down the decay of shutdown pressure.A smaller number of opened clusters and perforations can significantly increase the perforation friction and decrease the overall level of shutdown pressure.A larger filtration coefficient may accelerate the fluid filtration in the fracture and hasten the drop of the shutdown pressure.The simulation method of shutdown pressure,as well as the analysis results,has important implications for the interpretation of hydraulic fracture parameters.展开更多
Karst fracture-cavity carbonate reservoirs,in which natural cavities are connected by natural fractures to form cavity clusters in many circumstances,have become significant fields of oil and gas exploration and explo...Karst fracture-cavity carbonate reservoirs,in which natural cavities are connected by natural fractures to form cavity clusters in many circumstances,have become significant fields of oil and gas exploration and exploitation.Proppant fracturing is considered as the best method for exploiting carbonate reservoirs;however,previous studies primarily focused on the effects of individual types of geological formations,such as natural fractures or cavities,on fracture propagation.In this study,true-triaxial physical simulation experiments were systematically performed under four types of stress difference conditions after the accurate prefabrication of four types of different fracture-cavity distributions in artificial samples.Subsequently,the interaction mechanism between the hydraulic fractures and fracture-cavity structures was systematically analyzed in combination with the stress distribution,cross-sectional morphology of the main propagation path,and three-dimensional visualization of the overall fracture network.It was found that the propagation of hydraulic fractures near the cavity was inhibited by the stress concentration surrounding the cavity.In contrast,a natural fracture with a smaller approach angle(0°and 30°)around the cavity can alleviate the stress concentration and significantly facilitate the connection with the cavity.In addition,the hydraulic fracture crossed the natural fracture at the 45°approach angle and bypassed the cavity under higher stress difference conditions.A new stimulation effectiveness evaluation index was established based on the stimulated reservoir area(SRA),tortuosity of the hydraulic fractures(T),and connectivity index(CI)of the cavities.These findings provide new insights into the fracturing design of carbonate reservoirs.展开更多
In this paper,a viscoelasticity-plastic damage constitutive equation for naturally fractured shale is deduced,coupling nonlinear tensile-shear mixed fracture mode.Dynamic perforation-erosion on fluid re-distribution a...In this paper,a viscoelasticity-plastic damage constitutive equation for naturally fractured shale is deduced,coupling nonlinear tensile-shear mixed fracture mode.Dynamic perforation-erosion on fluid re-distribution among multi-clusters are considered as well.DFN-FEM(discrete fracture network combined with finite element method)was developed to simulate the multi-cluster complex fractures propagation within temporary plugging fracturing(TPF).Numerical results are matched with field injection and micro-seismic monitoring data.Based on geomechanical characteristics of Weiyuan deep shale gas reservoir in Sichuan Basin,SW China,a multi-cluster complex fractures propagation model is built for TPF.To study complex fractures propagation and the permeability-enhanced region evolution,intersecting and competition mechanisms between the fractures before and after TPF treatment are revealed.Simulation results show that:fracture from middle cluster is restricted by the fractures from side-clusters,and side-clusters plugging is benefit for multi fractures propagation in uniformity;optimized TPF timing should be delayed within a higher density or strike of natural fractures;Within a reservoir-featured natural fractures distribution,optimized TPF timing for most clustered method is 2/3 of total fluid injection time as the optimal plugging time under different clustering modes.展开更多
This paper presents an integrated study from fracture propagation modeling to gas flow modeling and a correlation analysis to explore the key controlling factors of intensive volume fracturing.The fracture propagation...This paper presents an integrated study from fracture propagation modeling to gas flow modeling and a correlation analysis to explore the key controlling factors of intensive volume fracturing.The fracture propagation model takes into account the interaction between hydraulic fracture and natural fracture by means of the displacement discontinuity method(DDM)and the Picard iterative method.The shale gas flow considers multiple transport mechanisms,and the flow in the fracture network is handled by the embedded discrete fracture model(EDFM).A series of numerical simulations are conducted to analyze the effects of the cluster number,stage spacing,stress difference coefficient,and natural fracture distribution on the stimulated fracture area,fractal dimension,and cumulative gas production,and their correlation coefficients are obtained.The results show that the most influential factors to the stimulated fracture area are the stress difference ratio,stage spacing,and natural fracture density,while those to the cumulative gas production are the stress difference ratio,natural fracture density,and cluster number.This indicates that the stress condition dominates the gas production,and employing intensive volume fracturing(by properly increasing the cluster number)is beneficial for improving the final cumulative gas production.展开更多
Hydraulic fracturing is widely used in geothermal resource exploitation, and many natural fractures exist in hot dry rock reservoirs due to in-situ stress and faults. However, the infuence of natural fractures on hydr...Hydraulic fracturing is widely used in geothermal resource exploitation, and many natural fractures exist in hot dry rock reservoirs due to in-situ stress and faults. However, the infuence of natural fractures on hydraulic fracture propagation is not considered in the current study. In this paper, based on the phase feld model, a thermo-hydro-mechanical coupled hydraulic fracture propagation model was established to reveal the infuence of injection time, fracturing method, injection fow rate, and natural fracture distribution on the fracture propagation mechanism. The results show that fracture complexity increases with an increase in injection time. The stress disturbance causes the fracture initiation pressure of the second cluster signifcantly higher than that of the frst and third clusters. The zipper-type fracturing method can reduce the degree of stress disturbance and increase fracture complexity by 7.2% compared to simultaneous hydraulic fracturing. Both low and high injection fow rate lead to a decrease in fracture propagation time, which is not conducive to an increase in fracture complexity. An increase in the natural fracture angle leads to hydraulic fracture crossing natural fracture, but has a lesser efect on fracture complexity. In this paper, we analyzed the infuence of diferent factors on initiation pressure and fracture complexity, providing valuable guidance for the exploitation of geothermal resources.展开更多
This paper establishes a 3D multi-well pad fracturing numerical model coupled with fracture propagation and proppant migration based on the displacement discontinuity method and Eulerian-Eulerian frameworks,and the fr...This paper establishes a 3D multi-well pad fracturing numerical model coupled with fracture propagation and proppant migration based on the displacement discontinuity method and Eulerian-Eulerian frameworks,and the fracture propagation and proppant distribution during multi-well fracturing are investigated by taking the actual multi-well pad parameters as an example.Fracture initiation and propagation during multi-well pad fracturing are jointly affected by a variety of stress interference mechanisms such as inter-cluster,inter-stage,and inter-well,and the fracture extension is unbalanced among clusters,asymmetric on both wings,and dipping at heels.Due to the significant influence of fracture morphology and width on the migration capacity of proppant in the fracture,proppant is mainly placed in the area near the wellbore with large fracture width,while a high-concentration sandwash may easily occur in the area with narrow fracture width as a result of quick bridging.On the whole,the proppant placement range is limited.Increasing the well-spacing can reduce the stress interference of adjacent wells and promote the uniform distribution of fractures and proppant on both wings.The maximum stimulated reservoir volume or multi-fracture uniform propagation can be achieved by optimizing the well spacing.Although reducing the perforation-cluster spacing also can improve the stimulated reservoir area,a too low cluster spacing is not conducive to effectively increasing the propped fracture area.Since increasing the stage time lag is beneficial to reduce inter-stage stress interference,zipper fracturing produces more uniform fracture propagation and proppant distribution.展开更多
A fracture propagation model of radial well fracturing is established based on the finite element-meshless method.The model considers the coupling effect of fracturing fluid flow and rock matrix deformation.The fractu...A fracture propagation model of radial well fracturing is established based on the finite element-meshless method.The model considers the coupling effect of fracturing fluid flow and rock matrix deformation.The fracture geometries of radial well fracturing are simulated,the induction effect of radial well on the fracture is quantitatively characterized,and the influences of azimuth,horizontal principle stress difference,and reservoir matrix permeability on the fracture geometries are revealed.The radial wells can induce the fractures to extend parallel to their axes when two radial wells in the same layer are fractured.When the radial wells are symmetrically distributed along the direction of the minimum horizontal principle stress with the azimuth greater than 15,the extrusion effect reduces the fracture length of radial wells.When the radial wells are symmetrically distributed along the direction of the maximum horizontal principal stress,the extrusion increases the fracture length of the radial wells.The fracture geometries are controlled by the rectification of radial borehole,the extrusion between radial wells in the same layer,and the deflection of the maximum horizontal principal stress.When the radial wells are distributed along the minimum horizontal principal stress symmetrically,the fracture length induced by the radial well decreases with the increase of azimuth;in contrast,when the radial wells are distributed along the maximum horizontal principal stress symmetrically,the fracture length induced by the radial well first decreases and then increases with the increase of azimuth.The fracture length induced by the radial well decreases with the increase of horizontal principal stress difference.The increase of rock matrix permeability and pore pressure of the matrix around radial wells makes the inducing effect of the radial well on fractures increase.展开更多
How to simulate fracture mode and crack propagation path in a plate with multiple cracks is an attractive but difficult issue in fracture mechanics.Peridynamics is a recently developed nonlocal continuum formulation t...How to simulate fracture mode and crack propagation path in a plate with multiple cracks is an attractive but difficult issue in fracture mechanics.Peridynamics is a recently developed nonlocal continuum formulation that can spontaneously predict the crack nucleation,branch and propagation in materials and structures through a meshfree discrete technique.In this paper,the peridynamic motion equation with boundary traction is improved by simplifying the boundary transfer functions.We calculate the critical cracking load and the fracture angles of the plate with multiple cracks under uniaxial tension.The results are consistent with those predicted by classical fracture mechanics.The fracture mode and crack propagation path are also determined.The calculation shows that the brittle fracture process of the plate with multiple cracks can be conveniently and correctly simulated by the peridynamic motion equation with boundary conditions.展开更多
In this study, an orthogonal array experiment is conducted by using a transparent fracture network replica. Image processing and theoretical analysis are performed to investigate the model sealing efficiency(SE), fact...In this study, an orthogonal array experiment is conducted by using a transparent fracture network replica. Image processing and theoretical analysis are performed to investigate the model sealing efficiency(SE), factors influencing SE, and the effect of flowing water on propagation. The results show that grout propagation can be classified into three patterns in the fracture network: sealing off, partial sealing,and major erosion. The factors controlling the SE in a descending order of the amount of influence are the initial water flow speed, fracture aperture, grout take, and gel time. An optimal value for the combination of the gel time and grout take(artificial factors) can result in a good SE. The grouting and seepage pressures are measured, and the results reveal that their variations can indicate the SE to some extent. The SE is good when the seepage pressure at each point increases overall;the frequent fluctuations in the seepage pressure indicate a moderately poor SE, and an overall decline in the seepage pressure indicates a major erosion type. The deflection effect of grouting shows an approximately elliptical propagation with the long axis expanding along the wider fracture opening, demonstrating further application in grouting design.展开更多
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.展开更多
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.展开更多
A novel laboratory simulation method for modeling multi-staged fracturing in a horizontal well was established based on a true tri-axial hydraulic fracturing simulation system. Using this method, the influences of net...A novel laboratory simulation method for modeling multi-staged fracturing in a horizontal well was established based on a true tri-axial hydraulic fracturing simulation system. Using this method, the influences of net pressure in hydraulic fracture, stage spacing, perforation parameter, horizontal stress bias and well cementation quality on the propagation geometry of multiple fractures in a tight sandstone formation were studied in detail. The specimen splitting and analogy analysis of fracturing curve patterns reveals: Multiple fractures tend to merge under the condition of high horizontal stress bias and short stage spacing with pre-existing hydraulic fractures under critical closure situation, and the propagation of subsequent fractures is possibly suppressed because of high net pressure in pre-created fractures and asymmetric distribution of fracture width. And the subsequently created fractures are situated in the induced stress decreasing zone due to long stage spacing, leading to weak stress interference, and perforation with intense density and deep penetration facilitates the decrease of initiation fracture pressure. The deflection angle of subsequent fracture and horizontal stress variation tend to be amplified under low horizontal bias with constant net pressure in fractures. The longitudinal fracture is likely to be initiated at the interface of wellbore and concrete sample with poor cementation quality. The initiation fracture pressure of the different stages increases in turn, with the largest increase of 30%. Pressure quickly declines after initiation with low propagation pressure when the transverse hydraulic fracture is formed. The pressure reduces with fluctuation after the initiation of fracture when the fracture deflects, the extension pressure is high, and the fracture formed is tortuous and narrow. There is a violently fluctuant rise of pressure with multiple peak values when longitudinal fracture created, and it is hard to distinguish the features between the initiation stage and propagation stage.展开更多
Objective As the most widely used and effective technique in reservoir reconstruction of unconventional natural gas,hydraulic fracturing has been achieved good effect in CBM development.It is important to note that co...Objective As the most widely used and effective technique in reservoir reconstruction of unconventional natural gas,hydraulic fracturing has been achieved good effect in CBM development.It is important to note that coal seam is both source rock and reservoir,展开更多
基金funded by the National Natural Science Foundation of China(Nos.52192622,51874253,U20A202)
文摘To investigate the height growth of multi-cluster fractures during variable fluid-viscosity fracturing in a layered shale oil reservoir,a two-dimensional finite element method(FEM)-discrete fracture network(DFN)model coupled with flow,stress and damage is proposed.A traction-separation law is used to describe the mixed-mode response of the damaged adhesive fractures,and the cubic law is used to describe the fluid flow within the fractures.The rock deformation is controlled by the in-situ stress,fracture cohesion and fluid pressure on the hydraulic fracture surface.The coupled finite element equations are solved by the explicit time difference method.The effects of the fracturing treatment parameters including fluid viscosity,pumping rate and cluster spacing on the geometries of multifractures are investigated.The results show that variable fluid-viscosity injection can improve the complexity of the fracture network and height of the main fractures simultaneously.The pumping rate of15 m^(3)/min,variable fluid-viscosity of 3-9-21-36-45 mPa s with a cluster spacing of 7.5 m is the ideal treatment strategy.The field application shows that the peak daily production of the application well with the optimized injection procedu re of variable fluid-viscosity fracturing is 171 tons(about 2.85 times that of the adjacent well),which is the highest daily production record of a single shale oil well in China,marking a strategic breakthrough of commercial shale oil production in the Jiyang Depression,Shengli Oilfield.The variable fluid-viscosity fracturing technique is proved to be very effective for improving shale oil production.
基金financial support by the National Key Research and Development Program of China (No.2022YFE0129800)the National Natural Science Foundation of China (No.52074311)。
文摘The continental shale reservoirs of Jurassic Lianggaoshan Formation in Sichuan Basin contain thin lamina,which is characterized by strong plasticity and developed longitudinal shell limestone interlayer.To improve the production efficiency of reservoirs by multi-cluster fracturing,it is necessary to consider the unbalanced propagation of hydraulic fractures and the penetration effect of fractures.This paper constructed a numerical model of multi-fracture propagation and penetration based on the finite element coupling cohesive zone method;considering the construction cluster spacing,pump rate,lamina strength and other parameters studied the influencing factors of multi-cluster fracture interaction propagation;combined with AE energy data and fracture mode reconstruction method,quantitatively characterized the comprehensive impact of the strength of thin interlayer rock interfaces on the initiation and propagation of fractures that penetrate layers,and accurately predicted the propagation pattern of hydraulic fractures through laminated shale oil reservoirs.Simulation results revealed that in the process of multi-cluster fracturing,the proportion of shear damage is low,and mainly occurs in bedding fractures activated by outer fractures.Reducing the cluster spacing enhances the fracture system's penetration ability,though it lowers the activation efficiency of lamina.The high plasticity of the limestone interlayer may impact the vertical propagation distance of the main fracture.Improving the interface strength is beneficial to the reconstruction of the fracture height,but the interface communication effect is limited.Reasonable selection of layers with moderate lamina strength for fracturing stimulation,increasing the pump rate during fracturing and setting the cluster spacing reasonably are beneficial to improve the effect of reservoir stimulation.
基金Supported by the National Natural Science Foundation of China(51974332).
文摘This study conducted temporary plugging and diversion fracturing(TPDF)experiments using a true triaxial fracturing simulation system within a laboratory setting that replicated a lab-based horizontal well completion with multi-cluster sand jetting perforation.The effects of temporary plugging agent(TPA)particle size,TPA concentration,single-cluster perforation number and cluster number on plugging pressure,multi-fracture diversion pattern and distribution of TPAs were investigated.A combination of TPAs with small particle sizes within the fracture and large particle sizes within the segment is conducive to increasing the plugging pressure and promoting the diversion of multi-fractures.The addition of fibers can quickly achieve ultra-high pressure,but it may lead to longitudinal fractures extending along the wellbore.The temporary plugging peak pressure increases with an increase in the concentration of the TPA,reaching a peak at a certain concentration,and further increases do not significantly improve the temporary plugging peak pressure.The breaking pressure and temporary plugging peak pressure show a decreasing trend with an increase in single-cluster perforation number.A lower number of single-cluster perforations is beneficial for increasing the breaking pressure and temporary plugging peak pressure,and it has a more significant control on the propagation of multi-cluster fractures.A lower number of clusters is not conducive to increasing the total number and complexity of artificial fractures,while a higher number of clusters makes it difficult to achieve effective plugging.The TPAs within the fracture is mainly concentrated in the complex fracture areas,especially at the intersections of fractures.Meanwhile,the TPAs within the segment are primarily distributed near the perforation cluster apertures which initiated complex fractures.
基金the funding support from the National Natural Science Foundation of China(Nos.52274093 and 52004236).
文摘The exploitation of shale gas is promising due to depletion of the conventional energy and intensification of the greenhouse effect.In this paper,we proposed a heat-fluid-solid coupling damage model of supercritical CO_(2)(SC-CO_(2))compound fracturing which is expected to be an efficient and environmentally friendly way to develop shale gas.The coupling model is solved by the finite element method,and the results are in good agreement with the analytical solutions and fracturing experiments.Based on this model,the fracture propagation characteristics at the two stages of compound fracturing are studied and the influence of pressurization rate,in situ stress,bedding angle,and other factors are considered.The results show that at the SC-CO_(2)fracturing stage,a lower pressurization rate is conducive to formation of the branches around main fractures,while a higher pressurization rate inhibits formation of the branches around main fractures and promotes formation of the main fractures.Both bedding and in situ stress play a dominant role in the fracture propagation.When the in situ stress ratio(δ_(x)/δ_(y))is 1,the presence of bedding can reduce the initiation pressure and failure pressure.Nevertheless,it will cause the fracture to propagate along the bedding direction,reducing the fracture complexity.In rocks without bedding,hydraulic fracturing has the lengthening and widening effects for SC-CO_(2)induced fracture.In shale,fractures induced at the hydraulic fracturing stage are more likely to be dominated by in situ stresses and have a shorter reorientation radius.Therefore,fracture branches propagating along the maximum principal stress direction may be generated around the main fractures induced by SC-CO_(2)at the hydraulic fracturing stage.When the branches converge with the main fractures,fracture zones are easily formed,and thus the fracture complexity and damage area can be significantly increased.The results are instructive for the design and application of SC-CO_(2)compound fracturing.
基金National Natural Science Foundation of China(51934005,U23B2089)Shaanxi Provincial Natural Science Basic Research Program Project(2024JC-YBQN-0554).
文摘The flow of fluid through the porous matrix of a reservoir rock applies a seepage force to the solid rock matrix.Although the seepage force exerted by fluid flow through the porous matrix of a reservoir rock has a notable influence on rock deformation and failure,its effect on hydraulic fracture(HF)propagation remains ambiguous.Therefore,in this study,we improved a traditional fluid–solid coupling method by incorporating the role of seepage force during the fracturing fluid seepage,using the discrete element method.First,we validated the simulation results of the improved method by comparing them with an analytical solution of the seepage force and published experimental results.Next,we conducted numerical simulations in both homogeneous and heterogeneous sandstone formations to investigate the influence of seepage force on HF propagation.Our results indicate that fluid viscosity has a greater impact on the magnitude and extent of seepage force compared to injection rate,and that lower viscosity and injection rate correspond to shorter hydraulic fracture lengths.Furthermore,seepage force influences the direction of HF propagation,causing HFs to deflect towards the side of the reservoir with weaker cementation and higher permeability.
基金supported by the National Natural Science Foundation of China(NSFC)(Grant Nos.42377156,42077251 and 42202305).
文摘Deep shale reservoirs are characterized by elevated breakdown pressures,diminished fracture complexity,and reduced modified volumes compared to medium and shallow reservoirs.Therefore,it is urgent to investigate particular injection strategies that can optimize breakdown pressure and fracturing efficiency to address the increasing demands for deep shale reservoir stimulation.In this study,the efficiency of various stimulation strategies,including multi-cluster simultaneous fracturing,modified alternating fracturing,alternating shut-in fracturing,and cyclic alternating fracturing,was evaluated.Subsequently,the sensitivity of factors such as the cycle index,shut-in time,cluster spacing,and horizontal permeability was investigated.Additionally,the flow distribution effect within the wellbore was discussed.The results indicate that relative to multi-cluster simultaneous fracturing,modified alternating fracturing exhibits reduced susceptibility to the stress shadow effect,which results in earlier breakdown,extended hydraulic fracture lengths,and more consistent propagation despite an increase in breakdown pressure.The alternating shut-in fracturing benefits the increase of fracture length,which is closely related to the shut-in time.Furthermore,cyclic alternating fracturing markedly lowers breakdown pressure and contributes to uniform fracture propagation,in which the cycle count plays an important role.Modified alternating fracturing demonstrates insensitivity to variations in cluster spacing,whereas horizontal permeability is a critical factor affecting fracture length.The wellbore effect restrains the accumulation of pressure and flow near the perforation,delaying the initiation of hydraulic fractures.The simulation results can provide valuable numerical insights for optimizing injection strategies for deep shale hydraulic fracturing.
基金National Natural Science Foundation of China(11672333).
文摘Indirect fracturing in the roof of broken soft coal seams has been demonstrated to be a feasible technology.In this work,the No.5 coal seam in the Hancheng block was taken as the research object.Based on the findings of true triaxial hydraulic fracturing experiments and field pilot under this technology and the cohesive element method,a 3D numerical model of indirect fracturing in the roof of broken soft coal seams was established,the fracture morphology propagation and evolution law under different conditions was investigated,and analysis of main controlling factors of fracture parameters was conducted with the combination weight method,which was based on grey incidence,analytic hierarchy process and entropy weight method.The results show that“士”-shaped fractures,T-shaped fractures,cross fractures,H-shaped fractures,and“干”-shaped fractures dominated by horizontal fractures were formed.Different parameter combinations can form different fracture morphologies.When the coal seam permeability is lower and the minimum horizontal principal stress difference between layers and fracturing fluid injection rate are both larger,it tends to form“士”-shaped fractures.When the coal seam permeability and minimum horizontal principal stress between layers and perforation position are moderate,cross fractures are easily generated.Different fracture parameters have different main controlling factors.Engineering factors of perforation location,fracturing fluid injection rate and viscosity are the dominant factors of hydraulic fracture shape parameters.This study can provide a reference for the design of indirect fracturing in the roof of broken soft coal seams.
基金The work is supported by the Sub-Project of“Research on Key Technologies and Equipment of Reservoir Stimulation”of China National Petroleum Corporation Post–14th Five-Year Plan Forward-Looking Major Science and Technology Project“Research on New Technology of Monitoring and Diagnosis of Horizontal Well Hydraulic Fracture Network Distribution Pattern”(2021DJ4502).
文摘Multistage multi-cluster hydraulic fracturing has enabled the economic exploitation of shale reservoirs,but the interpretation of hydraulic fracture parameters is challenging.The pressure signals after pump shutdown are influenced by hydraulic fractures,which can reflect the geometric features of hydraulic fracture.The shutdown pressure can be used to interpret the hydraulic fracture parameters in a real-time and cost-effective manner.In this paper,a mathematical model for shutdown pressure evolution is developed considering the effects of wellbore friction,perforation friction and fluid loss in fractures.An efficient numerical simulation method is established by using the method of characteristics.Based on this method,the impacts of fracture half-length,fracture height,opened cluster and perforation number,and filtration coefficient on the evolution of shutdown pressure are analyzed.The results indicate that a larger fracture half-length may hasten the decay of shutdown pressure,while a larger fracture height can slow down the decay of shutdown pressure.A smaller number of opened clusters and perforations can significantly increase the perforation friction and decrease the overall level of shutdown pressure.A larger filtration coefficient may accelerate the fluid filtration in the fracture and hasten the drop of the shutdown pressure.The simulation method of shutdown pressure,as well as the analysis results,has important implications for the interpretation of hydraulic fracture parameters.
基金sponsored by the National Natural Science Foundation of China(Grants Nos.52104046 and 52104010).
文摘Karst fracture-cavity carbonate reservoirs,in which natural cavities are connected by natural fractures to form cavity clusters in many circumstances,have become significant fields of oil and gas exploration and exploitation.Proppant fracturing is considered as the best method for exploiting carbonate reservoirs;however,previous studies primarily focused on the effects of individual types of geological formations,such as natural fractures or cavities,on fracture propagation.In this study,true-triaxial physical simulation experiments were systematically performed under four types of stress difference conditions after the accurate prefabrication of four types of different fracture-cavity distributions in artificial samples.Subsequently,the interaction mechanism between the hydraulic fractures and fracture-cavity structures was systematically analyzed in combination with the stress distribution,cross-sectional morphology of the main propagation path,and three-dimensional visualization of the overall fracture network.It was found that the propagation of hydraulic fractures near the cavity was inhibited by the stress concentration surrounding the cavity.In contrast,a natural fracture with a smaller approach angle(0°and 30°)around the cavity can alleviate the stress concentration and significantly facilitate the connection with the cavity.In addition,the hydraulic fracture crossed the natural fracture at the 45°approach angle and bypassed the cavity under higher stress difference conditions.A new stimulation effectiveness evaluation index was established based on the stimulated reservoir area(SRA),tortuosity of the hydraulic fractures(T),and connectivity index(CI)of the cavities.These findings provide new insights into the fracturing design of carbonate reservoirs.
基金Supported by the National Natural Science Foundation of China(52192622,52204005,U20A20265)Sichuan Outstanding Young Scientific and Technological Talents Project(2022JDJQ0007).
文摘In this paper,a viscoelasticity-plastic damage constitutive equation for naturally fractured shale is deduced,coupling nonlinear tensile-shear mixed fracture mode.Dynamic perforation-erosion on fluid re-distribution among multi-clusters are considered as well.DFN-FEM(discrete fracture network combined with finite element method)was developed to simulate the multi-cluster complex fractures propagation within temporary plugging fracturing(TPF).Numerical results are matched with field injection and micro-seismic monitoring data.Based on geomechanical characteristics of Weiyuan deep shale gas reservoir in Sichuan Basin,SW China,a multi-cluster complex fractures propagation model is built for TPF.To study complex fractures propagation and the permeability-enhanced region evolution,intersecting and competition mechanisms between the fractures before and after TPF treatment are revealed.Simulation results show that:fracture from middle cluster is restricted by the fractures from side-clusters,and side-clusters plugging is benefit for multi fractures propagation in uniformity;optimized TPF timing should be delayed within a higher density or strike of natural fractures;Within a reservoir-featured natural fractures distribution,optimized TPF timing for most clustered method is 2/3 of total fluid injection time as the optimal plugging time under different clustering modes.
基金supported by the National Natural Science Foundation of China(Nos.52274038,5203401042174143)+1 种基金the Taishan Scholars Project(No.tsqnz20221140)the Open Fund of State Key Laboratory of Oil and Gas Reservoir Geology and Exploitation(Southwest Petroleum University)of China(No.PLN2020-5)。
文摘This paper presents an integrated study from fracture propagation modeling to gas flow modeling and a correlation analysis to explore the key controlling factors of intensive volume fracturing.The fracture propagation model takes into account the interaction between hydraulic fracture and natural fracture by means of the displacement discontinuity method(DDM)and the Picard iterative method.The shale gas flow considers multiple transport mechanisms,and the flow in the fracture network is handled by the embedded discrete fracture model(EDFM).A series of numerical simulations are conducted to analyze the effects of the cluster number,stage spacing,stress difference coefficient,and natural fracture distribution on the stimulated fracture area,fractal dimension,and cumulative gas production,and their correlation coefficients are obtained.The results show that the most influential factors to the stimulated fracture area are the stress difference ratio,stage spacing,and natural fracture density,while those to the cumulative gas production are the stress difference ratio,natural fracture density,and cluster number.This indicates that the stress condition dominates the gas production,and employing intensive volume fracturing(by properly increasing the cluster number)is beneficial for improving the final cumulative gas production.
基金supported by the National Natural Science Foundation of China(52174024).
文摘Hydraulic fracturing is widely used in geothermal resource exploitation, and many natural fractures exist in hot dry rock reservoirs due to in-situ stress and faults. However, the infuence of natural fractures on hydraulic fracture propagation is not considered in the current study. In this paper, based on the phase feld model, a thermo-hydro-mechanical coupled hydraulic fracture propagation model was established to reveal the infuence of injection time, fracturing method, injection fow rate, and natural fracture distribution on the fracture propagation mechanism. The results show that fracture complexity increases with an increase in injection time. The stress disturbance causes the fracture initiation pressure of the second cluster signifcantly higher than that of the frst and third clusters. The zipper-type fracturing method can reduce the degree of stress disturbance and increase fracture complexity by 7.2% compared to simultaneous hydraulic fracturing. Both low and high injection fow rate lead to a decrease in fracture propagation time, which is not conducive to an increase in fracture complexity. An increase in the natural fracture angle leads to hydraulic fracture crossing natural fracture, but has a lesser efect on fracture complexity. In this paper, we analyzed the infuence of diferent factors on initiation pressure and fracture complexity, providing valuable guidance for the exploitation of geothermal resources.
基金Supported by National Natural Science Foundation of China(51974332)Strategic Cooperation Project Between PetroChina and China University of Petroleum(Beijing)(ZLZX2020-07).
文摘This paper establishes a 3D multi-well pad fracturing numerical model coupled with fracture propagation and proppant migration based on the displacement discontinuity method and Eulerian-Eulerian frameworks,and the fracture propagation and proppant distribution during multi-well fracturing are investigated by taking the actual multi-well pad parameters as an example.Fracture initiation and propagation during multi-well pad fracturing are jointly affected by a variety of stress interference mechanisms such as inter-cluster,inter-stage,and inter-well,and the fracture extension is unbalanced among clusters,asymmetric on both wings,and dipping at heels.Due to the significant influence of fracture morphology and width on the migration capacity of proppant in the fracture,proppant is mainly placed in the area near the wellbore with large fracture width,while a high-concentration sandwash may easily occur in the area with narrow fracture width as a result of quick bridging.On the whole,the proppant placement range is limited.Increasing the well-spacing can reduce the stress interference of adjacent wells and promote the uniform distribution of fractures and proppant on both wings.The maximum stimulated reservoir volume or multi-fracture uniform propagation can be achieved by optimizing the well spacing.Although reducing the perforation-cluster spacing also can improve the stimulated reservoir area,a too low cluster spacing is not conducive to effectively increasing the propped fracture area.Since increasing the stage time lag is beneficial to reduce inter-stage stress interference,zipper fracturing produces more uniform fracture propagation and proppant distribution.
基金Supported by the National Natural Science Foundation of China(51827804)CNPC Strategic Cooperation Science and Technology Major Project(ZLZX2020-01-05)Open Fund of State Key Laboratory of Rock Mechanics and Engineering(SKLGME021024).
文摘A fracture propagation model of radial well fracturing is established based on the finite element-meshless method.The model considers the coupling effect of fracturing fluid flow and rock matrix deformation.The fracture geometries of radial well fracturing are simulated,the induction effect of radial well on the fracture is quantitatively characterized,and the influences of azimuth,horizontal principle stress difference,and reservoir matrix permeability on the fracture geometries are revealed.The radial wells can induce the fractures to extend parallel to their axes when two radial wells in the same layer are fractured.When the radial wells are symmetrically distributed along the direction of the minimum horizontal principle stress with the azimuth greater than 15,the extrusion effect reduces the fracture length of radial wells.When the radial wells are symmetrically distributed along the direction of the maximum horizontal principal stress,the extrusion increases the fracture length of the radial wells.The fracture geometries are controlled by the rectification of radial borehole,the extrusion between radial wells in the same layer,and the deflection of the maximum horizontal principal stress.When the radial wells are distributed along the minimum horizontal principal stress symmetrically,the fracture length induced by the radial well decreases with the increase of azimuth;in contrast,when the radial wells are distributed along the maximum horizontal principal stress symmetrically,the fracture length induced by the radial well first decreases and then increases with the increase of azimuth.The fracture length induced by the radial well decreases with the increase of horizontal principal stress difference.The increase of rock matrix permeability and pore pressure of the matrix around radial wells makes the inducing effect of the radial well on fractures increase.
基金The work was supported by the National Nature Science Foundation of China through the Grant Nos.12072145 and 11672129.
文摘How to simulate fracture mode and crack propagation path in a plate with multiple cracks is an attractive but difficult issue in fracture mechanics.Peridynamics is a recently developed nonlocal continuum formulation that can spontaneously predict the crack nucleation,branch and propagation in materials and structures through a meshfree discrete technique.In this paper,the peridynamic motion equation with boundary traction is improved by simplifying the boundary transfer functions.We calculate the critical cracking load and the fracture angles of the plate with multiple cracks under uniaxial tension.The results are consistent with those predicted by classical fracture mechanics.The fracture mode and crack propagation path are also determined.The calculation shows that the brittle fracture process of the plate with multiple cracks can be conveniently and correctly simulated by the peridynamic motion equation with boundary conditions.
基金supported by the Natural Science Foundation of China under (Nos. 42172293, 4190020747, and 41472268)。
文摘In this study, an orthogonal array experiment is conducted by using a transparent fracture network replica. Image processing and theoretical analysis are performed to investigate the model sealing efficiency(SE), factors influencing SE, and the effect of flowing water on propagation. The results show that grout propagation can be classified into three patterns in the fracture network: sealing off, partial sealing,and major erosion. The factors controlling the SE in a descending order of the amount of influence are the initial water flow speed, fracture aperture, grout take, and gel time. An optimal value for the combination of the gel time and grout take(artificial factors) can result in a good SE. The grouting and seepage pressures are measured, and the results reveal that their variations can indicate the SE to some extent. The SE is good when the seepage pressure at each point increases overall;the frequent fluctuations in the seepage pressure indicate a moderately poor SE, and an overall decline in the seepage pressure indicates a major erosion type. The deflection effect of grouting shows an approximately elliptical propagation with the long axis expanding along the wider fracture opening, demonstrating further application in grouting design.
基金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.
基金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.
文摘A novel laboratory simulation method for modeling multi-staged fracturing in a horizontal well was established based on a true tri-axial hydraulic fracturing simulation system. Using this method, the influences of net pressure in hydraulic fracture, stage spacing, perforation parameter, horizontal stress bias and well cementation quality on the propagation geometry of multiple fractures in a tight sandstone formation were studied in detail. The specimen splitting and analogy analysis of fracturing curve patterns reveals: Multiple fractures tend to merge under the condition of high horizontal stress bias and short stage spacing with pre-existing hydraulic fractures under critical closure situation, and the propagation of subsequent fractures is possibly suppressed because of high net pressure in pre-created fractures and asymmetric distribution of fracture width. And the subsequently created fractures are situated in the induced stress decreasing zone due to long stage spacing, leading to weak stress interference, and perforation with intense density and deep penetration facilitates the decrease of initiation fracture pressure. The deflection angle of subsequent fracture and horizontal stress variation tend to be amplified under low horizontal bias with constant net pressure in fractures. The longitudinal fracture is likely to be initiated at the interface of wellbore and concrete sample with poor cementation quality. The initiation fracture pressure of the different stages increases in turn, with the largest increase of 30%. Pressure quickly declines after initiation with low propagation pressure when the transverse hydraulic fracture is formed. The pressure reduces with fluctuation after the initiation of fracture when the fracture deflects, the extension pressure is high, and the fracture formed is tortuous and narrow. There is a violently fluctuant rise of pressure with multiple peak values when longitudinal fracture created, and it is hard to distinguish the features between the initiation stage and propagation stage.
基金supported by the National Natural Science Foundation of China(grant No.41572140)the National Major Special Project of Science and Technology of China(grant No.2016ZX05044-001)+1 种基金the Fundamental Research Funds for the Central Universities(grant No.2015XKZD07)the Qing Lan Project
文摘Objective As the most widely used and effective technique in reservoir reconstruction of unconventional natural gas,hydraulic fracturing has been achieved good effect in CBM development.It is important to note that coal seam is both source rock and reservoir,