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.展开更多
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 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.展开更多
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.展开更多
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.展开更多
Due to the presence of ice and unfrozen water in pores of frozen rock,the rock fracture behaviors are susceptible to temperature.In this study,the potential thawing-induced softening effects on the fracture behaviors ...Due to the presence of ice and unfrozen water in pores of frozen rock,the rock fracture behaviors are susceptible to temperature.In this study,the potential thawing-induced softening effects on the fracture behaviors of frozen rock is evaluated by testing the tension fracture toughness(KIC)of frozen rock at different temperatures(i.e.-20℃,-15℃,-12℃,-10℃,-8℃,-6℃,-4℃,-2℃,and 0℃).Acoustic emission(AE)and digital image correlation(DIC)methods are utilized to analyze the microcrack propagation during fracturing.The melting of pore ice is measured using nuclear magnetic resonance(NMR)method.The results indicate that:(1)The KIC of frozen rock decreases moderately between-20℃ and-4℃,and rapidly between-4℃ and 0℃.(2)At-20℃ to-4℃,the fracturing process,deduced from the DIC results at the notch tip,exhibits three stages:elastic deformation,microcrack propagation and microcrack coalescence.However,at-4℃e0℃,only the latter two stages are observed.(3)At-4℃e0℃,the AE activities during fracturing are less than that at-20℃ to-4℃,while more small events are reported.(4)The NMR results demonstrate a reverse variation trend in pore ice content with increasing temperature,that is,a moderate decrease is followed by a sharp decrease and-4℃ is exactly the critical temperature.Next,we interpret the thawing-induced softening effect by linking the evolution in microscopic structure of frozen rock with its macroscopic fracture behaviors as follow:from-20℃ to-4℃,the thickening of the unfrozen water film diminishes the cementation strength between ice and rock skeleton,leading to the decrease in fracture parameters.From-4℃ to 0℃,the cementation effect of ice almost vanishes,and the filling effect of pore ice is reduced significantly,which facilitates microcrack propagation and thus the easier fracture of frozen rocks.展开更多
As a calculation method based on the Galerkin variation,the numerical manifold method(NMM)adopts a double covering system,which can easily deal with discontinuous deformation problems and has a high calculation accura...As a calculation method based on the Galerkin variation,the numerical manifold method(NMM)adopts a double covering system,which can easily deal with discontinuous deformation problems and has a high calculation accuracy.Aiming at the thermo-mechanical(TM)coupling problem of fractured rock masses,this study uses the NMM to simulate the processes of crack initiation and propagation in a rock mass under the influence of temperature field,deduces related system equations,and proposes a penalty function method to deal with boundary conditions.Numerical examples are employed to confirm the effectiveness and high accuracy of this method.By the thermal stress analysis of a thick-walled cylinder(TWC),the simulation of cracking in the TWC under heating and cooling conditions,and the simulation of thermal cracking of the SwedishÄspöPillar Stability Experiment(APSE)rock column,the thermal stress,and TM coupling are obtained.The numerical simulation results are in good agreement with the test data and other numerical results,thus verifying the effectiveness of the NMM in dealing with thermal stress and crack propagation problems of fractured rock masses.展开更多
Radial borehole fracturing that combines radial boreholes with hydraulic fracturing is anticipated to improve the output of tight oil and gas reservoirs.This paper aims to investigate fracture propagation and pressure...Radial borehole fracturing that combines radial boreholes with hydraulic fracturing is anticipated to improve the output of tight oil and gas reservoirs.This paper aims to investigate fracture propagation and pressure characteristics of radial borehole fracturing in multiple layers.A series of laboratory experiments with artificial rock samples(395 mm×395 mm×395 mm)was conducted using a true triaxial fracturing device.Three crucial factors corresponding to the vertical distance of adjacent radial borehole layers(vertical distance),the azimuth and diameter of the radial borehole are examined.Experimental results show that radial borehole fracturing in multiple layers generates diverse fracture geometries.Four types of fractures are identified based on the connectivity between hydraulic fractures and radial boreholes.The vertical distance significantly influences fracture propagation perpendicular to the radial borehole axis.An increase in the vertical distance impedes fracture connection across multiple radial borehole layers and reduces the fracture propagation distance along the radial borehole axis.The azimuth also influences fracture propagation along the radial borehole axis.Increasing the azimuth reduces the guiding ability of radial boreholes,which makes the fracture quickly curve to the maximum horizontal stress direction.The breakdown pressure correlates with diverse fracture geometries observed.When the fractures connect multi-layer radial boreholes,increasing the vertical distance decreases the breakdown pressure.Decreasing the azimuth and increasing the diameter also decrease the breakdown pressure.The extrusion force exists between the adjacent fractures generated in radial boreholes in multiple rows,which plays a crucial role in enhancing the guiding ability of radial boreholes and results in higher breakdown pressure.The research provides valuable theoretical insights for the field application of radial borehole fracturing technology in tight oil and gas reservoirs.展开更多
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 order to study and analyze the stability of engineering rock mass under non-uniform triaxial stress and obtain the evolution mechanism of the whole process of fracture,a series of conventional triaxial compression ...In order to study and analyze the stability of engineering rock mass under non-uniform triaxial stress and obtain the evolution mechanism of the whole process of fracture,a series of conventional triaxial compression tests and three-dimensional numerical simulation tests were carried out on hollow granite specimens with different diameters.The bearing capacity of hollow cylindrical specimen is analyzed based on elasticity.The results show that:1)Under low confining pressure,the tensile strain near the hole of the hollow cylindrical specimen is obvious,and the specimen deformation near the hole is significant.At the initial stage of loading,the compressive stress and compressive strain of the specimen are widely distributed.With the progress of loading,the number of microelements subjected to tensile strain gradually increases,and even spreads throughout the specimen;2)Under conventional triaxial compression,the cracking position of hollow cylinder specimens is concentrated in the upper and lower parts,and the final fracture mode is generally compressive shear failure.The final fracture mode of complete specimen is generally tensile fracture.Under high confining pressure,the tensile cracks of the sample are concentrated in the upper and lower parts and are not connected,while the cracks of the upper and lower parts of the intact sample will expand and connect to form a fracture surface;3)In addition,the tensile crack widths of intact and hollow cylindrical specimens under low confining pressure are larger than those under high confining pressure.展开更多
基金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.
基金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.
基金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.
基金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.
基金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.
基金We acknowledge the funding support from the National Natural Science Foundation of China(Grant No.42271148).
文摘Due to the presence of ice and unfrozen water in pores of frozen rock,the rock fracture behaviors are susceptible to temperature.In this study,the potential thawing-induced softening effects on the fracture behaviors of frozen rock is evaluated by testing the tension fracture toughness(KIC)of frozen rock at different temperatures(i.e.-20℃,-15℃,-12℃,-10℃,-8℃,-6℃,-4℃,-2℃,and 0℃).Acoustic emission(AE)and digital image correlation(DIC)methods are utilized to analyze the microcrack propagation during fracturing.The melting of pore ice is measured using nuclear magnetic resonance(NMR)method.The results indicate that:(1)The KIC of frozen rock decreases moderately between-20℃ and-4℃,and rapidly between-4℃ and 0℃.(2)At-20℃ to-4℃,the fracturing process,deduced from the DIC results at the notch tip,exhibits three stages:elastic deformation,microcrack propagation and microcrack coalescence.However,at-4℃e0℃,only the latter two stages are observed.(3)At-4℃e0℃,the AE activities during fracturing are less than that at-20℃ to-4℃,while more small events are reported.(4)The NMR results demonstrate a reverse variation trend in pore ice content with increasing temperature,that is,a moderate decrease is followed by a sharp decrease and-4℃ is exactly the critical temperature.Next,we interpret the thawing-induced softening effect by linking the evolution in microscopic structure of frozen rock with its macroscopic fracture behaviors as follow:from-20℃ to-4℃,the thickening of the unfrozen water film diminishes the cementation strength between ice and rock skeleton,leading to the decrease in fracture parameters.From-4℃ to 0℃,the cementation effect of ice almost vanishes,and the filling effect of pore ice is reduced significantly,which facilitates microcrack propagation and thus the easier fracture of frozen rocks.
基金supported by the National Natural Science Foundation of China(Grant No.42277165)the Fundamental Research Funds for the Central Universities,China University of Geosciences(Wuhan)(Grant No.CUGCJ1821)the National Overseas Study Fund(Grant No.202106410040).
文摘As a calculation method based on the Galerkin variation,the numerical manifold method(NMM)adopts a double covering system,which can easily deal with discontinuous deformation problems and has a high calculation accuracy.Aiming at the thermo-mechanical(TM)coupling problem of fractured rock masses,this study uses the NMM to simulate the processes of crack initiation and propagation in a rock mass under the influence of temperature field,deduces related system equations,and proposes a penalty function method to deal with boundary conditions.Numerical examples are employed to confirm the effectiveness and high accuracy of this method.By the thermal stress analysis of a thick-walled cylinder(TWC),the simulation of cracking in the TWC under heating and cooling conditions,and the simulation of thermal cracking of the SwedishÄspöPillar Stability Experiment(APSE)rock column,the thermal stress,and TM coupling are obtained.The numerical simulation results are in good agreement with the test data and other numerical results,thus verifying the effectiveness of the NMM in dealing with thermal stress and crack propagation problems of fractured rock masses.
基金The authors gratefully acknowledge the financial support received from the Strategic Cooperation Technology Projects of CNPC and CUPB(No.ZLZX2020-01)the National Key Scientific Research Instrument Research Project of NSFC(No.51827804).
文摘Radial borehole fracturing that combines radial boreholes with hydraulic fracturing is anticipated to improve the output of tight oil and gas reservoirs.This paper aims to investigate fracture propagation and pressure characteristics of radial borehole fracturing in multiple layers.A series of laboratory experiments with artificial rock samples(395 mm×395 mm×395 mm)was conducted using a true triaxial fracturing device.Three crucial factors corresponding to the vertical distance of adjacent radial borehole layers(vertical distance),the azimuth and diameter of the radial borehole are examined.Experimental results show that radial borehole fracturing in multiple layers generates diverse fracture geometries.Four types of fractures are identified based on the connectivity between hydraulic fractures and radial boreholes.The vertical distance significantly influences fracture propagation perpendicular to the radial borehole axis.An increase in the vertical distance impedes fracture connection across multiple radial borehole layers and reduces the fracture propagation distance along the radial borehole axis.The azimuth also influences fracture propagation along the radial borehole axis.Increasing the azimuth reduces the guiding ability of radial boreholes,which makes the fracture quickly curve to the maximum horizontal stress direction.The breakdown pressure correlates with diverse fracture geometries observed.When the fractures connect multi-layer radial boreholes,increasing the vertical distance decreases the breakdown pressure.Decreasing the azimuth and increasing the diameter also decrease the breakdown pressure.The extrusion force exists between the adjacent fractures generated in radial boreholes in multiple rows,which plays a crucial role in enhancing the guiding ability of radial boreholes and results in higher breakdown pressure.The research provides valuable theoretical insights for the field application of radial borehole fracturing technology in tight oil and gas reservoirs.
基金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.
基金Projects(52074116,51804113)supported by the National Natural Science Foundation of China。
文摘In order to study and analyze the stability of engineering rock mass under non-uniform triaxial stress and obtain the evolution mechanism of the whole process of fracture,a series of conventional triaxial compression tests and three-dimensional numerical simulation tests were carried out on hollow granite specimens with different diameters.The bearing capacity of hollow cylindrical specimen is analyzed based on elasticity.The results show that:1)Under low confining pressure,the tensile strain near the hole of the hollow cylindrical specimen is obvious,and the specimen deformation near the hole is significant.At the initial stage of loading,the compressive stress and compressive strain of the specimen are widely distributed.With the progress of loading,the number of microelements subjected to tensile strain gradually increases,and even spreads throughout the specimen;2)Under conventional triaxial compression,the cracking position of hollow cylinder specimens is concentrated in the upper and lower parts,and the final fracture mode is generally compressive shear failure.The final fracture mode of complete specimen is generally tensile fracture.Under high confining pressure,the tensile cracks of the sample are concentrated in the upper and lower parts and are not connected,while the cracks of the upper and lower parts of the intact sample will expand and connect to form a fracture surface;3)In addition,the tensile crack widths of intact and hollow cylindrical specimens under low confining pressure are larger than those under high confining pressure.