Prepulse combined hydraulic fracturing facilitates the development of fracture networks by integrating prepulse hydraulic loading with conventional hydraulic fracturing.The formation mechanisms of fracture networks be...Prepulse combined hydraulic fracturing facilitates the development of fracture networks by integrating prepulse hydraulic loading with conventional hydraulic fracturing.The formation mechanisms of fracture networks between hydraulic and pre-existing fractures under different prepulse loading parameters remain unclear.This research investigates the impact of prepulse loading parameters,including the prepulse loading number ratio(C),prepulse loading stress ratio(S),and prepulse loading frequency(f),on the formation of fracture networks between hydraulic and pre-existing fractures,using both experimental and numerical methods.The results suggest that low prepulse loading stress ratios and high prepulse loading number ratios are advantageous loading modes.Multiple hydraulic fractures are generated in the specimen under the advantageous loading modes,facilitating the development of a complex fracture network.Fatigue damage occurs in the specimen at the prepulse loading stage.The high water pressure at the secondary conventional hydraulic fracturing promotes the growth of hydraulic fractures along the damage zones.This allows the hydraulic fractures to propagate deeply and interact with pre-existing fractures.Under advantageous loading conditions,multiple hydraulic fractures can extend to pre-existing fractures,and these hydraulic fractures penetrate or propagate along pre-existing fractures.Especially when the approach angle is large,the damage range in the specimen during the prepulse loading stage increases,resulting in the formation of more hydraulic fractures.展开更多
In order to identify the development characteristics of fracture network in tight conglomerate reservoir of Mahu after hydraulic fracturing,a hydraulic fracturing test site was set up in the second and third members o...In order to identify the development characteristics of fracture network in tight conglomerate reservoir of Mahu after hydraulic fracturing,a hydraulic fracturing test site was set up in the second and third members of Triassic Baikouquan Formation(T1b2 and T1b3)in Ma-131 well area,which learned from the successful experience of hydraulic fracturing test sites in North America(HFTS-1).Twelve horizontal wells and a high-angle coring well MaJ02 were drilled.The orientation,connection,propagation law and major controlling factors of hydraulic fractures were analyzed by comparing results of CT scans,imaging logs,direct observation of cores from Well MaJ02,and combined with tracer monitoring data.Results indicate that:(1)Two types of fractures have developed by hydraulic fracturing,i.e.tensile fractures and shear fractures.Tensile fractures are approximately parallel to the direction of the maximum horizontal principal stress,and propagate less than 50 m from perforation clusters.Shear fractures are distributed among tensile fractures and mainly in the strike-slip mode due to the induced stress field among tensile fractures,and some of them are in conjugated pairs.Overall,tensile fractures alternate with shear fractures,with shear fractures dominated and activated after tensile ones.(2)Tracer monitoring results indicate that communication between wells was prevalent in the early stage of production,and the static pressure in the fracture gradually decreased and the connectivity between wells reduced as production progressed.(3)Density of hydraulic fractures is mainly affected by the lithology and fracturing parameters,which is smaller in the mudstone than the conglomerate.Larger fracturing scale and smaller cluster spacing lead to a higher fracture density,which are important directions to improve the well productivity.展开更多
Multi-stage and multi-cluster fracturing(MMF)is a crucial technology in unconventional oil and gas development,aiming to enhance production by creating extensive fracture networks.However,achieving uniform expansion o...Multi-stage and multi-cluster fracturing(MMF)is a crucial technology in unconventional oil and gas development,aiming to enhance production by creating extensive fracture networks.However,achieving uniform expansion of multi-cluster hydraulic fractures(HFs)in MMF remains a significant challenge.Field practice has shown that the use of temporary plugging and diversion fracturing(TPDF)can promote the balanced expansion of multi-cluster HFs.This study conducted TPDF experiments using a true triaxial fracturing simulation system setting a horizontal well completion with multi-cluster jetting perforations to investigate the equilibrium initiation and extension of multi-cluster fractures.The influence of key parameters,including cluster spacing,fracturing fluid viscosity,differential stress,and fracturing fluid injection rate,on fracture initiation and propagation was systematically examined.The results indicate that while close-spaced multi-cluster fracturing significantly increases the number of HFs,it also leads to uneven extension of HFs in their propagation.In contrast,TPDF demonstrates effectiveness in mitigating uneven HF extension,increasing the number of HFs,and creating a larger stimulated reservoir volume,ultimately leading to improved oil and gas well productivity.Moreover,under conditions of high differential stress,the differential stress within the formation exerts a stronger guiding effect in HFs,which are more closely aligned with the minimum principal stress.Low-viscosity fluids facilitate rapid and extensive fracture propagation within the rock formation.High-volume fluid injection,on the other hand,more comprehensively fills the formation.Therefore,employing lowviscosity and high-volume fracturing is advantageous for the initiation and extension of multi-cluster HFs.展开更多
The core sampling experiments were conducted after hydraulic fracturing in the three-dimensional development zone of Fuling shale gas.Six coring wells of different well types were systematically designed.Based on the ...The core sampling experiments were conducted after hydraulic fracturing in the three-dimensional development zone of Fuling shale gas.Six coring wells of different well types were systematically designed.Based on the integrated engineering technology of post-fracturing drilling,coring and monitoring of shale and the analysis of fracture source tracing,the evaluation of the fracture network after fracturing in the three-dimensional development of shale gas was conducted.The data of core fractures after fracturing indicate that three major types of fractures are formed after fracturing:natural fractures,hydraulic fractures,and fractures induced by external mechanical force,which are further classified into six subcategories:natural structural fractures,natural bedding fractures,hydraulic fractures,hydraulically activated fractures,drilling induced fractures,and fractures induced by core transportation.The forms of the artificial fracture network after fracturing are complex.Hydraulic fractures and hydraulically activated fractures interweave with each other,presenting eight forms of artificial fracture networks,among which the“一”-shaped fracture is the most common,accounting for approximately 70%of the total fractures.When the distance to the fractured wellbore is less than 35 m,the density of the artificial fracture network is relatively high;when it is 35–100 m,the density is lower;and when it is beyond 100 m,the density gradually increases.The results of the fracture tracing in the core sampling area confirm that the current fracturing technology can essentially achieve the differential transformation of the reservoir in the main area of Jiaoshiba block in Fuling.The three-layer three-dimensional development model can efficiently utilize shale gas reserves,although there is still room for improvement in the complexity and propagation uniformity of fractures.It is necessary to further optimize technologies such as close-cutting combined with temporary blocking and deflection within fractures or at fracture mouths,as well as limited flow perforation techniques,to promote the balanced initiation and extension of fractures.展开更多
The real-time monitoring of fracture propagation during hydraulic fracturing is crucial for obtaining a deeper understanding of fracture morphology and optimizing hydraulic fracture designs.Accurate measurements of ke...The real-time monitoring of fracture propagation during hydraulic fracturing is crucial for obtaining a deeper understanding of fracture morphology and optimizing hydraulic fracture designs.Accurate measurements of key fracture parameters,such as the fracture height and width,are particularly important to ensure efficient oilfield development and precise fracture diagnosis.This study utilized the optical frequency domain reflectometer(OFDR)technique in physical simulation experiments to monitor fractures during indoor true triaxial hydraulic fracturing experiments.The results indicate that the distributed fiber optic strain monitoring technology can efficiently capture the initiation and expansion of fractures.In horizontal well monitoring,the fiber strain waterfall plot can be used to interpret the fracture width,initiation location,and expansion speed.The fiber response can be divided into three stages:strain contraction convergence,strain band formation,and postshutdown strain rate reversal.When the fracture does not contact the fiber,a dual peak strain phenomenon occurs in the fiber and gradually converges as the fracture approaches.During vertical well monitoring in adjacent wells,within the effective monitoring range of the fiber,the axial strain produced by the fiber can represent the fracture height with an accuracy of 95.6%relative to the actual fracture height.This study provides a new perspective on real-time fracture monitoring.The response patterns of fiber-induced strain due to fractures can help us better understand and assess the dynamic fracture behavior,offering significant value for the optimization of oilfield development and fracture diagnostic techniques.展开更多
With the increasing demand for energy,traditional oil resources are facing depletion and insufficient supply.Many countries are rapidly turning to the development of unconventional oil and gas resources.Among them,sha...With the increasing demand for energy,traditional oil resources are facing depletion and insufficient supply.Many countries are rapidly turning to the development of unconventional oil and gas resources.Among them,shale oil and gas reservoirs have become the focus of unconventional oil and gas resources exploration and development.Based on the characteristics of shale oil and gas reservoirs,supercritical CO_(2) fracturing is more conducive to improving oil recovery than other fracturing technologies.In this paper,the mechanism of fracture initiation and propagation of supercritical CO_(2) in shale is analyzed,including viscosity effect,surface tension effect,permeation diffusion effect of supercritical CO_(2),and dissolution-adsorption effect between CO_(2) and shale.The effects of natural factors,such as shale properties,bedding plane and natural fractures,and controllable factors,proppant,temperature,pressure,CO_(2) concentration and injection rate on fracture initiation and propagation are clarified.The methods of supercritical CO_(2) fracturing process,thickener and proppant optimization to improve the efficiency of supercritical CO_(2) fracturing are discussed.In addition,some new technologies of supercritical CO_(2) fracturing are introduced.The challenges and prospects in the current research are also summarized.For example,supercritical CO_(2) is prone to filtration when passing through porous media,and it is difficult to form a stable flow state.Therefore,in order to achieve stable fracturing fluid suspension and effectively support fractu res,it is urge nt to explo re new fracturing fluid additives or improve fracturing fluid formulations combined with the research of new proppants.This paper is of great significance for understanding the behavior mechanism of supercritical CO_(2) in shale and optimizing fracturing technology.展开更多
A precise diagnosis of the complex post-fracturing characteristics and parameter variations in tight gas reservoirs is essential for optimizing fracturing technology,enhancing treatment effectiveness,and assessing pos...A precise diagnosis of the complex post-fracturing characteristics and parameter variations in tight gas reservoirs is essential for optimizing fracturing technology,enhancing treatment effectiveness,and assessing post-fracturing production capacity.Tight gas reservoirs face challenges due to the interaction between natural fractures and induced fractures.To address these issues,a theoretical model for diagnosing fractures under varying leak-off mechanisms has been developed,incorporating the closure behavior of natural fractures.This model,grounded in material balance theory,also accounts for shut-in pressure.The study derived and plotted typical G-function charts,which capture fracture behavior during closure.By superimposing the G-function in the closure phase of natural fractures with pressure derivative curves,the study explored how fracture parameters—including leak-off coefficient,fracture area,closure pressure,and closure time—impact these diagnostic charts.Findings show that variations in natural fracture flexibility,fracture area,and controlling factors influence the superimposed G-function pressure derivative curve,resulting in distinctive“concave”or“convex”patterns.Field data from Well Y in a specific tight gas reservoir were used to validate the model,confirming both its reliability and practicality.展开更多
This study aims to understand the effect of injection rate on injection-induced fracture activation in granite.We performed water injection-induced slip tests on samples containing either a smooth or a rough fracture ...This study aims to understand the effect of injection rate on injection-induced fracture activation in granite.We performed water injection-induced slip tests on samples containing either a smooth or a rough fracture at four different injection rates under undrained conditions and monitored the acoustic emission(AE)signals during the tests.Experimental results reveal that the critical activation fluid pressure is related to the injection rate,pressure diffusion rate,stress state,and fracture roughness.For the smooth fracture,as the injection rate increases,the critical activation fluid pressure increases significantly,while the injection rate has little effect on the critical activation fluid pressure of the rough fracture.The quasi-static slip distance of fractures decreases as the injection rate increases,with rough fractures exhibiting a greater overall slip distance compared to smooth fractures.The number of AE events per unit sliding distance increases with the injection rate,while the global b value decreases.These results indicate that higher injection rates produce more large-magnitude AE events and more severe slip instability and asperity damage.We established a linkage between fluid injection volume,injection rate,and AE events using the seismogenic index(Σ).The smooth fracture exhibits a steadily increasingΣwith the elapse of injection time,and the rate of increase is higher at higher injection rates;while the rough fracture is featured by a fluctuatingΣ,signifying the intermittent occurrence of large-magnitude AE events associated with the damage of larger fracture asperities.Our results highlight the importance of fracture surface heterogeneity on injection-induced fracture activation and slip.展开更多
Hydraulic fracturing is a crucial technique for efficient development of coal reservoirs.Coal rocks typically contain a high density of natural fractures,which serve as conduits for fracturing fluid.Upon injection,the...Hydraulic fracturing is a crucial technique for efficient development of coal reservoirs.Coal rocks typically contain a high density of natural fractures,which serve as conduits for fracturing fluid.Upon injection,the fluid infiltrates these natural fractures and leaks out,resulting in complex fracture morphology.The prediction of hydraulic fracture network propagation for coal reservoirs has important practical significance for evaluating hydraulic fracturing.This study proposes a novel inversion method for predicting fracture networks in coal reservoirs,explicitly considering the distribution of natural fractures.The method incorporates three distinct natural fracture opening modes and employs a fractal probability function to constrain fracture propagationmorphology.Based on thismethod,the study compares hydraulic fracture networkmorphologies in coal reservoirs with andwithout the presence of natural fractures.Theresults showthatwhile both reservoir types exhibitmulti-branch fracture networks,reservoirs containing natural fractures demonstrate greater branching and a larger stimulated reservoir volume(SRV).Additionally,the study employs a fractal dimension calculation method to quantitatively describe the geometric distribution characteristics of fractures.The analysis reveals that the geometry and distribution of natural fractures,as well as reservoir geological parameters,significantly influence the fracture networkmorphology and fractal dimension.The contact angle between natural and hydraulic fractures affects propagation direction;specifically,when the contact angle isπ/2,the fractal dimension of the hydraulic fracture network is maximized.Moreover,smaller lengths and spacings of natural fracture led to higher fractal dimensions,which can significantly increase the SRV.The proposed method offers an effective tool for evaluating the hydraulic fracturing of coal reservoirs.展开更多
Hydraulic fracturing technology has played an important role in the exploitation of unconventional oil and gas resources,however,its application to gas hydrate reservoirs has been rarely studied.Currently,there is sti...Hydraulic fracturing technology has played an important role in the exploitation of unconventional oil and gas resources,however,its application to gas hydrate reservoirs has been rarely studied.Currently,there is still limited understanding of the propagation and extension of fractures around the wellbore during the fracturing process of horizontal wells in hydrate reservoirs,as well as the stress interference patterns between fractures.This study simulates hydraulic fracturing processes in hydrate reservoirs using a fluidsolid coupling discrete element method(DEM),and analyzes the impacts of hydrate saturation and geological and engineering factors on fracture extension and stress disturbance.The results show that hydraulic fracturing is more effective when hydrate saturation exceeds 30%and that fracture pressure increases with saturation.The increase in horizontal stress differential enhances the directionality of fracture propagation and reduces stress disturbance.The distribution uniformity index(DUI)reveals that injection pressure is directly proportional to the number of main fractures and inversely proportional to fracturing time,with fracturing efficiency depending on the spacing between injection points and the distance between wells.This work may provide reference for the commercial exploitation of natural gas hydrates.展开更多
Understanding the anisotropic fracture behavior and the characteristics of the fracture process zone(FPZ)under size effects in laminated rocks,as well as its role in rock fracturing,is crucial for various engineering ...Understanding the anisotropic fracture behavior and the characteristics of the fracture process zone(FPZ)under size effects in laminated rocks,as well as its role in rock fracturing,is crucial for various engineering applications.In this study,three-point bending tests were conducted on shale specimens with varying bedding angles and sizes.The anisotropic characteristics and size effects of fracture parameters were revealed.A comparative analysis was performed on the evolutions of FPZs computed using size effect theory,digital image correlation(DIC),and linear elastic fracture mechanics.The results divulged that:(i)With increasing bedding angles,there is a noticeable decrease in apparent fracture toughness(KICA),apparent fracture energy(GICA),and nominal strength(σ_(Nu)).When the bedding angle of shale is less than 45°,the crack propagation and fracture parameters are mainly influenced by the matrix.Contrary,shale with bedding angles greater than 60°,the crack propagation and fracture parameters are mainly controlled by the bedding.When the bedding angle is between 45°and 60°,the fracture propagation evolves from permeating the matrix to extending along the bedding;(ii)The fracture parameters exhibit significant size dependent behavior,as KICA and GICA rise with increasing specimen size,butσNu falls with increasing specimen sizes.The fracture parameters align with the theoretical predictions of Bažant size effect law;and(iii)The lengths of DIC-based FPZ,effective FPZ,and inelastic zone follow W-shape variations with bedding angle.The dimensionless sizes of FPZ and inelastic zone decrease with specimen size,indicating a size effect.Furthermore,there is a negative relation between KICA and the dimensionless size of the FPZ,whileσNu is positively correlated to the dimensionless size of the FPZ.This highlights the essential role of the FPZ in the size effect of rock fracture.The bedding angle exerts an influence on the FPZ,subsequently affecting the anisotropic fracture and size-dependent behavior of shale.展开更多
Coupled thermo-hydro-mechanical(THM)processes in fractured rock are playing a crucial role in geoscience and geoengineering applications.Diverse and conceptually distinct approaches have emerged over the past decades ...Coupled thermo-hydro-mechanical(THM)processes in fractured rock are playing a crucial role in geoscience and geoengineering applications.Diverse and conceptually distinct approaches have emerged over the past decades in both continuum and discontinuum perspectives leading to significant progress in their comprehending and modeling.This review paper offers an integrated perspective on existing modeling methodologies providing guidance for model selection based on the initial and boundary conditions.By comparing various models,one can better assess the uncertainties in predictions,particularly those related to the conceptual models.The review explores how these methodologies have significantlyenhanced the fundamental understanding of how fractures respond to fluid injection and production,and improved predictive capabilities pertaining to coupled processes within fractured systems.It emphasizes the importance of utilizing advanced computational technologies and thoroughly considering fundamental theories and principles established through past experimental evidence and practical experience.The selection and calibration of model parameters should be based on typical ranges and applied to the specificconditions of applications.The challenges arising from inherent heterogeneity and uncertainties,nonlinear THM coupled processes,scale dependence,and computational limitations in representing fieldscale fractures are discussed.Realizing potential advances on computational capacity calls for methodical conceptualization,mathematical modeling,selection of numerical solution strategies,implementation,and calibration to foster simulation outcomes that intricately reflectthe nuanced complexities of geological phenomena.Future research efforts should focus on innovative approaches to tackle the hurdles and advance the state-of-the-art in this critical fieldof study.展开更多
Based on the finite element-discrete element numerical method,a numerical model of fracture propagation in deflagration fracturing was established by considering the impact of stress wave,quasi-static pressure of expl...Based on the finite element-discrete element numerical method,a numerical model of fracture propagation in deflagration fracturing was established by considering the impact of stress wave,quasi-static pressure of explosive gas,and reflection of stress wave.The model was validated against the results of physical experiments.Taking the shale reservoirs of Silurian Longmaxi Formation in Luzhou area of the Sichuan Basin as an example,the effects of in-situ stress difference,natural fracture parameters,branch wellbore spacing,delay detonation time,and angle between branch wellbore and main wellbore on fracture propagation were identified.The results show that the fracture propagation morphology in deflagration fracturing is less affected by the in-situ stress difference when it is 5-15 MPa,and the tendency of fracture intersection between branch wellbores is significantly weakened when the in-situ stress difference reaches 20 MPa.The increase of natural fracture length promotes the fracture propagation along the natural fracture direction,while the increase of volumetric natural fracture density and angle limits the fracture propagation area and reduces the probability of fracture intersection between branch wells.The larger the branch wellbore spacing,the less probability of the fracture intersection between branch wells,allowing for the fracture propagation in multiple directions.Increasing the delay detonation time decreases the fracture spacing between branch wellbores.When the angle between the branch wellbore and the main wellbore is 45°and 90°,there is a tendency of fracture intersection between branch wellbores.展开更多
Shale reservoirs are characterized by numerous geological discontinuities,such as bedding planes,and exhibit pronounced heterogeneity across rock layers separated by these planes.Bedding planes often possess distinct ...Shale reservoirs are characterized by numerous geological discontinuities,such as bedding planes,and exhibit pronounced heterogeneity across rock layers separated by these planes.Bedding planes often possess distinct mechanical properties compared to the surrounding rock matrix,particularly in terms of damage and fracture behavior.Consequently,vertical propagation of hydraulic fractures is influenced by both bedding planes and the heterogeneity.In this study,a numerical investigation into the height growth of hydraulic fractures was conducted using the finite element method,incorporating zero-thickness cohesive elements.The analysis explored the effects of bedding planes,toughness contrasts between layers,and variations in in-situ stress across different strata.The results reveal that hydraulic fractures are more likely to propagate along bedding planes instead of traversing them and extending vertically into barrier layers when(1)bedding strength is low,(2)stress contrast between layers is high,and(3)toughness contrast is significant.Furthermore,for a given bedding strength,increased stress contrast or higher toughness contrast between layers elevate hydraulic fracture extension pressure.When a substantial stress difference exists between layers(Lc 0.4),hydraulic=fractures preferentially propagate along bedding planes.Conversely,as bedding strength increases,the propagation distance along bedding planes decreases,accompanied by an amplified horizontal compressive stress field.Notably,when the stress difference is sufficiently small(SD a phenomenon termed“stress rolling”emerges,wherein<-0.2),hydraulic fractures deviate from vertical growth and instead extend along a near-horizontal trajectory.展开更多
Based on continuum-discontinuum element method,the numerical simulation of fracture propagation during deflagration-hydraulic composite fracturing was constructed by considering deflagration stress impact induced frac...Based on continuum-discontinuum element method,the numerical simulation of fracture propagation during deflagration-hydraulic composite fracturing was constructed by considering deflagration stress impact induced fracture creation,deflagrating gas driven fracture propagation,and hydraulic fracture propagation,exploring the effects of in-situ stress difference,deflagration peak pressure,deflagration pressurization rate,hydraulic fracturing displacement and hydraulic fracturing fluid viscosity on fracture propagation in deflagration-hydraulic composite fracturing.The deflagration-hydraulic composite fracturing combines the advantages of deflagration fracturing in creating complex fractures near wells and the deep penetration of hydraulic fracturing at the far-field region,which can form multiple deep penetrating long fractures with better stimulation effects.With the increase of in-situ stress difference,the stimulated area of deflagration-hydraulic composite fracturing is reduced,and the deflagration-hydraulic composite fracturing is more suitable for reservoirs with small in-situ stress difference.Higher peak pressure and pressurization rate are conducive to increasing the maximum fracture length and burst degree of the deflagration fractures,which in turn increases the stimulated area of deflagration-hydraulic composite fracturing and improves the stimulation effect.Increasing the displacement and viscosity of hydraulic fracturing fluid can enhance the net pressure within the fractures,activate the deflagration fractures,increase the turning radius of the fractures,generate more long fractures,and effectively increase the stimulated reservoir area.The stimulated reservoir area is not completely positively correlated with the hydraulic fracturing displacement and fracturing fluid viscosity,and there is a critical value.When the critical value is exceeded,the stimulated area decreases.展开更多
Stereoscopic particle image velocimetry technology was employed to investigate the planar three-dimensional velocity field and the process of proppant entry into branch fractures in a fracture configuration of“vertic...Stereoscopic particle image velocimetry technology was employed to investigate the planar three-dimensional velocity field and the process of proppant entry into branch fractures in a fracture configuration of“vertical main fracture-vertical branch fracture”intersecting at a 90°angle.This study analyzed the effects of pumping rate,fracturing fluid viscosity,proppant particle size,and fracture width on the transport behavior of proppant into branch fractures.Based on the deflection behavior of proppant,the main fractures can be divided into five regions:pre-entry transition,pre-entry stabilization,deflection entry at the fracture mouth,rear absorption entry,and movement away from the fracture mouth.Proppant primarily deflects into the branch fracture at the fracture mouth,with a small portion drawn in from the rear of the intersection.Increasing the pumping rate,reducing the proppant particle size,and widening the branch fracture are conducive to promoting proppant deflection into the branch.With increasing fracturing fluid viscosity,the ability of proppant to enter the branch fracture first improves and then declines,indicating that excessively high viscosity is unfavorable for proppant entry into the branch.During field operations,a high pumping rate and micro-to small-sized proppant can be used in the early stage to ensure effective placement in the branch fractures,followed by medium-to large-sized proppant to ensure adequate placement in the main fracture and enhance the overall conductivity of the fracture network.展开更多
Accurate prediction of hydraulic fracture propagation is vital for Enhanced Geothermal System(EGS)design.We study the first hydraulic fracturing job at the GR1 well in the Gonghe Basin using field data,where the overa...Accurate prediction of hydraulic fracture propagation is vital for Enhanced Geothermal System(EGS)design.We study the first hydraulic fracturing job at the GR1 well in the Gonghe Basin using field data,where the overall direction of hydraulic fractures does not show a delineated shape parallel to the maximum principal stress orientation.A field-scale numerical model based on the distinct element method is set up to carry out a fully coupled hydromechanical simulation,with the explicit representation of natural fractures via the discrete fracture network(DFN)approach.The effects of injection parameters and in situ stress on hydraulic fracture patterns are then quantitatively assessed.The study reveals that shear-induced deformation primarily governs the fracturing morphology in the GR1 well,driven by smaller injection rates and viscosities that promote massive activation of natural fractures,ultimately dominating the direction of hydraulic fracturing.Furthermore,the increase of in situ differential stress may promote shear damage of natural fracture surfaces,with the exact influence pattern depending on the combination of specific discontinuity properties and in situ stress state.Finally,we provide recommendations for EGS fracturing based on the influence characteristics of multiple parameters.This study can serve as an effective basis and reference for the design and optimization of EGS in the Gonghe basin and other sites.展开更多
Controllable shock wave fracturing is an innovative engineering technique used for shale reservoir fracturing and reformation.Understanding the anisotropic fracture mechanism of shale under impact loading is vital for...Controllable shock wave fracturing is an innovative engineering technique used for shale reservoir fracturing and reformation.Understanding the anisotropic fracture mechanism of shale under impact loading is vital for optimizing shock wave fracturing equipment and enhancing shale oil production.In this study,using the well-known notched semi-circular bend(NSCB)sample and the novel double-edge notched flattened Brazilian disc(DNFBD)sample combined with a split Hopkinson pressure bar(SHPB),various dynamic anisotropic fracture properties of Lushan shale,including failure characteristics,fracture toughness,energy dissipation and crack propagation velocity,are comprehensively compared and discussed under mode Ⅰ and mode Ⅱ fracture scenarios.First,using a newly modified fracture criterion considering the strength anisotropy of shale,the DNFBD specimen is predicted to be a robust method for true mode Ⅱ fracture of anisotropic shale rocks.Our experimental results show that the dynamic mode Ⅱ fracture of shale induces a rougher and more complex fracture morphology and performs a higher fracture toughness or fracture energy compared to dynamic mode Ⅰ fracture.The minimal fracture toughness or fracture energy occurs in the Short-transverse orientation,while the maximal ones occur in the Divider orientation.In addition,it is interesting to find that the mode Ⅱ fracture toughness anisotropy index decreases more slowly than that in the mode Ⅰ fracture scenario.These results provide significant insights for understanding the different dynamic fracture mechanisms of anisotropic shale rocks under impact loading and have some beneficial implications for the controllable shock wave fracturing technique.展开更多
The linear elastic hydraulic fracture criterion is not applicable to deep reservoirs when nonlinear behavior is present over an extensive zone at the fracture tip.This study aims to develop a criterion for nonlinear h...The linear elastic hydraulic fracture criterion is not applicable to deep reservoirs when nonlinear behavior is present over an extensive zone at the fracture tip.This study aims to develop a criterion for nonlinear hydraulic fracture considering the fracture process zone(FPZ)and seeks to reveal the causes of nonlinearity during fracture propagation in deep reservoirs.A closing stress profile considering the in-situ stress was established by using the cohesive zone model(CZM)to describe the FPZ at the fracture tip.An analytical model for the FPZ length was derived,while the criterion for nonlinear fracture propagation was proposed.The FPZ fully developed and the fracture began to propagate when the apparent stress intensity at the fracture tip reached the apparent fracture toughness or when the in-situ stress intensity reached the in-situ fracture toughness.The proposed criterion can clearly determine the length of the FPZ,accurately predict the breakdown pressure during fracturing operations,and establish a relationship between these two parameters.It addresses the inherent limitations of conventional linear elastic fracture mechanics(LEFM),which often underestimates fracture toughness and neglects the effects of the FPZ.This research is expected to enhance the fracturing design in deep reservoirs.展开更多
Heterogeneity in carbonate formations due to discontinuities(e.g.,fractures and cavities)will bring about distinctive acid stimulation effects.However,the differences in fracturing behavior between homogeneous and het...Heterogeneity in carbonate formations due to discontinuities(e.g.,fractures and cavities)will bring about distinctive acid stimulation effects.However,the differences in fracturing behavior between homogeneous and heterogeneous carbonate formations remain unclear,complicating the optimization of acid fracturing strategies.In this paper,full-diameter carbonate rock samples with different degrees of discontinuity development are selected to investigate the fracturing behavior under different fluid types and injection schemes.Advanced techniques,including 3D CT scanning and 3D laser scanning,are employed to analyze fracture morphology and etching characteristics,respectively.Experimental results show that the coupled hydraulic-chemical effects play different roles in fracture induction between fracture-cavity developed and undeveloped carbonate rocks.Acid-fracturing stimulation consistently induces multiple types of complex fractures in fracture-cavity carbonate rocks,whereas it results in a single artificial fracture in less fracture-cavity carbonate rocks.Furthermore,localized etching patterns are prevalent in most fracture-cavity carbonate rocks,whereas homogeneous carbonate rocks exhibit regional or global etching characteristics.In both carbonate rocks,the stimulation effect of guar fluid is inferior to that of gelled acid but comparable to self-generating acid.Further findings are that alternating fracturing with guar and acid fluids in fracture-cavity carbonate rocks can sustain or even increase the injection pressure,facilitating the formation of new or depth-penetrating fractures.This phenomenon,however,is not observed in fracture-cavity undeveloped carbonate rocks.Potential interaction modes between induced fracture and natural fractures/cavities under different injection conditions are also identified.Finally,preferred fracturing schemes applicable to different carbonate formations are recommended based on the area,number and roughness of the induced fractures.展开更多
基金financially supported by,the Fundamental Research Funds for the Central Universities(Grant No.2023QN1064)the China Postdoctoral Science Foundation(Grant No.2023M733772)Jiangsu Funding Program for Excellent Postdoctoral Talent(Grant No.2023ZB847)。
文摘Prepulse combined hydraulic fracturing facilitates the development of fracture networks by integrating prepulse hydraulic loading with conventional hydraulic fracturing.The formation mechanisms of fracture networks between hydraulic and pre-existing fractures under different prepulse loading parameters remain unclear.This research investigates the impact of prepulse loading parameters,including the prepulse loading number ratio(C),prepulse loading stress ratio(S),and prepulse loading frequency(f),on the formation of fracture networks between hydraulic and pre-existing fractures,using both experimental and numerical methods.The results suggest that low prepulse loading stress ratios and high prepulse loading number ratios are advantageous loading modes.Multiple hydraulic fractures are generated in the specimen under the advantageous loading modes,facilitating the development of a complex fracture network.Fatigue damage occurs in the specimen at the prepulse loading stage.The high water pressure at the secondary conventional hydraulic fracturing promotes the growth of hydraulic fractures along the damage zones.This allows the hydraulic fractures to propagate deeply and interact with pre-existing fractures.Under advantageous loading conditions,multiple hydraulic fractures can extend to pre-existing fractures,and these hydraulic fractures penetrate or propagate along pre-existing fractures.Especially when the approach angle is large,the damage range in the specimen during the prepulse loading stage increases,resulting in the formation of more hydraulic fractures.
基金Supported by the National Natural Science Foundation of China(52274051)CNPC-China University of Petroleum(Beijing)Strategic Cooperative Project(ZLZX2020-01).
文摘In order to identify the development characteristics of fracture network in tight conglomerate reservoir of Mahu after hydraulic fracturing,a hydraulic fracturing test site was set up in the second and third members of Triassic Baikouquan Formation(T1b2 and T1b3)in Ma-131 well area,which learned from the successful experience of hydraulic fracturing test sites in North America(HFTS-1).Twelve horizontal wells and a high-angle coring well MaJ02 were drilled.The orientation,connection,propagation law and major controlling factors of hydraulic fractures were analyzed by comparing results of CT scans,imaging logs,direct observation of cores from Well MaJ02,and combined with tracer monitoring data.Results indicate that:(1)Two types of fractures have developed by hydraulic fracturing,i.e.tensile fractures and shear fractures.Tensile fractures are approximately parallel to the direction of the maximum horizontal principal stress,and propagate less than 50 m from perforation clusters.Shear fractures are distributed among tensile fractures and mainly in the strike-slip mode due to the induced stress field among tensile fractures,and some of them are in conjugated pairs.Overall,tensile fractures alternate with shear fractures,with shear fractures dominated and activated after tensile ones.(2)Tracer monitoring results indicate that communication between wells was prevalent in the early stage of production,and the static pressure in the fracture gradually decreased and the connectivity between wells reduced as production progressed.(3)Density of hydraulic fractures is mainly affected by the lithology and fracturing parameters,which is smaller in the mudstone than the conglomerate.Larger fracturing scale and smaller cluster spacing lead to a higher fracture density,which are important directions to improve the well productivity.
基金funded by the National Natural Science Foundation of China(52104046).
文摘Multi-stage and multi-cluster fracturing(MMF)is a crucial technology in unconventional oil and gas development,aiming to enhance production by creating extensive fracture networks.However,achieving uniform expansion of multi-cluster hydraulic fractures(HFs)in MMF remains a significant challenge.Field practice has shown that the use of temporary plugging and diversion fracturing(TPDF)can promote the balanced expansion of multi-cluster HFs.This study conducted TPDF experiments using a true triaxial fracturing simulation system setting a horizontal well completion with multi-cluster jetting perforations to investigate the equilibrium initiation and extension of multi-cluster fractures.The influence of key parameters,including cluster spacing,fracturing fluid viscosity,differential stress,and fracturing fluid injection rate,on fracture initiation and propagation was systematically examined.The results indicate that while close-spaced multi-cluster fracturing significantly increases the number of HFs,it also leads to uneven extension of HFs in their propagation.In contrast,TPDF demonstrates effectiveness in mitigating uneven HF extension,increasing the number of HFs,and creating a larger stimulated reservoir volume,ultimately leading to improved oil and gas well productivity.Moreover,under conditions of high differential stress,the differential stress within the formation exerts a stronger guiding effect in HFs,which are more closely aligned with the minimum principal stress.Low-viscosity fluids facilitate rapid and extensive fracture propagation within the rock formation.High-volume fluid injection,on the other hand,more comprehensively fills the formation.Therefore,employing lowviscosity and high-volume fracturing is advantageous for the initiation and extension of multi-cluster HFs.
文摘The core sampling experiments were conducted after hydraulic fracturing in the three-dimensional development zone of Fuling shale gas.Six coring wells of different well types were systematically designed.Based on the integrated engineering technology of post-fracturing drilling,coring and monitoring of shale and the analysis of fracture source tracing,the evaluation of the fracture network after fracturing in the three-dimensional development of shale gas was conducted.The data of core fractures after fracturing indicate that three major types of fractures are formed after fracturing:natural fractures,hydraulic fractures,and fractures induced by external mechanical force,which are further classified into six subcategories:natural structural fractures,natural bedding fractures,hydraulic fractures,hydraulically activated fractures,drilling induced fractures,and fractures induced by core transportation.The forms of the artificial fracture network after fracturing are complex.Hydraulic fractures and hydraulically activated fractures interweave with each other,presenting eight forms of artificial fracture networks,among which the“一”-shaped fracture is the most common,accounting for approximately 70%of the total fractures.When the distance to the fractured wellbore is less than 35 m,the density of the artificial fracture network is relatively high;when it is 35–100 m,the density is lower;and when it is beyond 100 m,the density gradually increases.The results of the fracture tracing in the core sampling area confirm that the current fracturing technology can essentially achieve the differential transformation of the reservoir in the main area of Jiaoshiba block in Fuling.The three-layer three-dimensional development model can efficiently utilize shale gas reserves,although there is still room for improvement in the complexity and propagation uniformity of fractures.It is necessary to further optimize technologies such as close-cutting combined with temporary blocking and deflection within fractures or at fracture mouths,as well as limited flow perforation techniques,to promote the balanced initiation and extension of fractures.
基金supported by the National Natural Science Foundation of China(Grant No.52104060)the Natural Science Foundation of Shandong Province,China(Grant No.ZR2021QE015).
文摘The real-time monitoring of fracture propagation during hydraulic fracturing is crucial for obtaining a deeper understanding of fracture morphology and optimizing hydraulic fracture designs.Accurate measurements of key fracture parameters,such as the fracture height and width,are particularly important to ensure efficient oilfield development and precise fracture diagnosis.This study utilized the optical frequency domain reflectometer(OFDR)technique in physical simulation experiments to monitor fractures during indoor true triaxial hydraulic fracturing experiments.The results indicate that the distributed fiber optic strain monitoring technology can efficiently capture the initiation and expansion of fractures.In horizontal well monitoring,the fiber strain waterfall plot can be used to interpret the fracture width,initiation location,and expansion speed.The fiber response can be divided into three stages:strain contraction convergence,strain band formation,and postshutdown strain rate reversal.When the fracture does not contact the fiber,a dual peak strain phenomenon occurs in the fiber and gradually converges as the fracture approaches.During vertical well monitoring in adjacent wells,within the effective monitoring range of the fiber,the axial strain produced by the fiber can represent the fracture height with an accuracy of 95.6%relative to the actual fracture height.This study provides a new perspective on real-time fracture monitoring.The response patterns of fiber-induced strain due to fractures can help us better understand and assess the dynamic fracture behavior,offering significant value for the optimization of oilfield development and fracture diagnostic techniques.
文摘With the increasing demand for energy,traditional oil resources are facing depletion and insufficient supply.Many countries are rapidly turning to the development of unconventional oil and gas resources.Among them,shale oil and gas reservoirs have become the focus of unconventional oil and gas resources exploration and development.Based on the characteristics of shale oil and gas reservoirs,supercritical CO_(2) fracturing is more conducive to improving oil recovery than other fracturing technologies.In this paper,the mechanism of fracture initiation and propagation of supercritical CO_(2) in shale is analyzed,including viscosity effect,surface tension effect,permeation diffusion effect of supercritical CO_(2),and dissolution-adsorption effect between CO_(2) and shale.The effects of natural factors,such as shale properties,bedding plane and natural fractures,and controllable factors,proppant,temperature,pressure,CO_(2) concentration and injection rate on fracture initiation and propagation are clarified.The methods of supercritical CO_(2) fracturing process,thickener and proppant optimization to improve the efficiency of supercritical CO_(2) fracturing are discussed.In addition,some new technologies of supercritical CO_(2) fracturing are introduced.The challenges and prospects in the current research are also summarized.For example,supercritical CO_(2) is prone to filtration when passing through porous media,and it is difficult to form a stable flow state.Therefore,in order to achieve stable fracturing fluid suspension and effectively support fractu res,it is urge nt to explo re new fracturing fluid additives or improve fracturing fluid formulations combined with the research of new proppants.This paper is of great significance for understanding the behavior mechanism of supercritical CO_(2) in shale and optimizing fracturing technology.
文摘A precise diagnosis of the complex post-fracturing characteristics and parameter variations in tight gas reservoirs is essential for optimizing fracturing technology,enhancing treatment effectiveness,and assessing post-fracturing production capacity.Tight gas reservoirs face challenges due to the interaction between natural fractures and induced fractures.To address these issues,a theoretical model for diagnosing fractures under varying leak-off mechanisms has been developed,incorporating the closure behavior of natural fractures.This model,grounded in material balance theory,also accounts for shut-in pressure.The study derived and plotted typical G-function charts,which capture fracture behavior during closure.By superimposing the G-function in the closure phase of natural fractures with pressure derivative curves,the study explored how fracture parameters—including leak-off coefficient,fracture area,closure pressure,and closure time—impact these diagnostic charts.Findings show that variations in natural fracture flexibility,fracture area,and controlling factors influence the superimposed G-function pressure derivative curve,resulting in distinctive“concave”or“convex”patterns.Field data from Well Y in a specific tight gas reservoir were used to validate the model,confirming both its reliability and practicality.
基金supported by the National Key Research and Development Program of China(Grant No.2023YFB2390300)the National Natural Science Foundation of China(Grant No.42172292)Yinlin Ji is grateful for the support by the Helmholtz Association's Initiative and Networking Fund for the Helmholtz Young Investigator Group ARES(Contract No.VH-NG-1516).
文摘This study aims to understand the effect of injection rate on injection-induced fracture activation in granite.We performed water injection-induced slip tests on samples containing either a smooth or a rough fracture at four different injection rates under undrained conditions and monitored the acoustic emission(AE)signals during the tests.Experimental results reveal that the critical activation fluid pressure is related to the injection rate,pressure diffusion rate,stress state,and fracture roughness.For the smooth fracture,as the injection rate increases,the critical activation fluid pressure increases significantly,while the injection rate has little effect on the critical activation fluid pressure of the rough fracture.The quasi-static slip distance of fractures decreases as the injection rate increases,with rough fractures exhibiting a greater overall slip distance compared to smooth fractures.The number of AE events per unit sliding distance increases with the injection rate,while the global b value decreases.These results indicate that higher injection rates produce more large-magnitude AE events and more severe slip instability and asperity damage.We established a linkage between fluid injection volume,injection rate,and AE events using the seismogenic index(Σ).The smooth fracture exhibits a steadily increasingΣwith the elapse of injection time,and the rate of increase is higher at higher injection rates;while the rough fracture is featured by a fluctuatingΣ,signifying the intermittent occurrence of large-magnitude AE events associated with the damage of larger fracture asperities.Our results highlight the importance of fracture surface heterogeneity on injection-induced fracture activation and slip.
基金This study was supported by the Foundation of National Engineering Laboratory for Exploration and Development of Low-Permeability Oil and Gas Fields(2023-015).
文摘Hydraulic fracturing is a crucial technique for efficient development of coal reservoirs.Coal rocks typically contain a high density of natural fractures,which serve as conduits for fracturing fluid.Upon injection,the fluid infiltrates these natural fractures and leaks out,resulting in complex fracture morphology.The prediction of hydraulic fracture network propagation for coal reservoirs has important practical significance for evaluating hydraulic fracturing.This study proposes a novel inversion method for predicting fracture networks in coal reservoirs,explicitly considering the distribution of natural fractures.The method incorporates three distinct natural fracture opening modes and employs a fractal probability function to constrain fracture propagationmorphology.Based on thismethod,the study compares hydraulic fracture networkmorphologies in coal reservoirs with andwithout the presence of natural fractures.Theresults showthatwhile both reservoir types exhibitmulti-branch fracture networks,reservoirs containing natural fractures demonstrate greater branching and a larger stimulated reservoir volume(SRV).Additionally,the study employs a fractal dimension calculation method to quantitatively describe the geometric distribution characteristics of fractures.The analysis reveals that the geometry and distribution of natural fractures,as well as reservoir geological parameters,significantly influence the fracture networkmorphology and fractal dimension.The contact angle between natural and hydraulic fractures affects propagation direction;specifically,when the contact angle isπ/2,the fractal dimension of the hydraulic fracture network is maximized.Moreover,smaller lengths and spacings of natural fracture led to higher fractal dimensions,which can significantly increase the SRV.The proposed method offers an effective tool for evaluating the hydraulic fracturing of coal reservoirs.
基金financially supported by the National Key Research and Development Plan(2023YFC2811001)the National Natural Science Foundation of China(42206233)the Taishan Scholars Program(tsqn202312280,tsqn202306297)。
文摘Hydraulic fracturing technology has played an important role in the exploitation of unconventional oil and gas resources,however,its application to gas hydrate reservoirs has been rarely studied.Currently,there is still limited understanding of the propagation and extension of fractures around the wellbore during the fracturing process of horizontal wells in hydrate reservoirs,as well as the stress interference patterns between fractures.This study simulates hydraulic fracturing processes in hydrate reservoirs using a fluidsolid coupling discrete element method(DEM),and analyzes the impacts of hydrate saturation and geological and engineering factors on fracture extension and stress disturbance.The results show that hydraulic fracturing is more effective when hydrate saturation exceeds 30%and that fracture pressure increases with saturation.The increase in horizontal stress differential enhances the directionality of fracture propagation and reduces stress disturbance.The distribution uniformity index(DUI)reveals that injection pressure is directly proportional to the number of main fractures and inversely proportional to fracturing time,with fracturing efficiency depending on the spacing between injection points and the distance between wells.This work may provide reference for the commercial exploitation of natural gas hydrates.
基金funded by the National Natural Science Foundation of China(Grant Nos.12172230,U22A20166)the Department of Science and Technology of Guangdong Province(Grant No.2019ZT08G315).
文摘Understanding the anisotropic fracture behavior and the characteristics of the fracture process zone(FPZ)under size effects in laminated rocks,as well as its role in rock fracturing,is crucial for various engineering applications.In this study,three-point bending tests were conducted on shale specimens with varying bedding angles and sizes.The anisotropic characteristics and size effects of fracture parameters were revealed.A comparative analysis was performed on the evolutions of FPZs computed using size effect theory,digital image correlation(DIC),and linear elastic fracture mechanics.The results divulged that:(i)With increasing bedding angles,there is a noticeable decrease in apparent fracture toughness(KICA),apparent fracture energy(GICA),and nominal strength(σ_(Nu)).When the bedding angle of shale is less than 45°,the crack propagation and fracture parameters are mainly influenced by the matrix.Contrary,shale with bedding angles greater than 60°,the crack propagation and fracture parameters are mainly controlled by the bedding.When the bedding angle is between 45°and 60°,the fracture propagation evolves from permeating the matrix to extending along the bedding;(ii)The fracture parameters exhibit significant size dependent behavior,as KICA and GICA rise with increasing specimen size,butσNu falls with increasing specimen sizes.The fracture parameters align with the theoretical predictions of Bažant size effect law;and(iii)The lengths of DIC-based FPZ,effective FPZ,and inelastic zone follow W-shape variations with bedding angle.The dimensionless sizes of FPZ and inelastic zone decrease with specimen size,indicating a size effect.Furthermore,there is a negative relation between KICA and the dimensionless size of the FPZ,whileσNu is positively correlated to the dimensionless size of the FPZ.This highlights the essential role of the FPZ in the size effect of rock fracture.The bedding angle exerts an influence on the FPZ,subsequently affecting the anisotropic fracture and size-dependent behavior of shale.
基金funding from the European Research Council(ERC)under the European Union’s Horizon 2020 Research and Innovation Program through the Starting Grant GEoREST(grant agreement No.801809)support by MICIU/AEI/10.13039/501100011033 and by"European Union Next Generation EU/PRTR"through the‘Ramón y Cajal’fellowship(reference RYC2021-032780-I)+9 种基金funding by MICIU/AEI/10.13039/501100011033 and by“ERDF,EU”through the‘HydroPoreII’project(reference PID2022-137652NBC44)support by the Institute for Korea Spent Nuclear Fuel(iKSNF)National Research Foundation of Korea(NRF)grant funded by the Korea government(Ministry of Science and ICT,MSIT)(2021M2E1A1085196)support by the Swedish Radiation Safety(SSM),Swedish Transport Administration(Trafikverket),Swedish Rock Engineering Foundation(BeFo),and Nordic Energy Research(Grant 187658)supported by the US Department of Energy(DOE),the Officeof Nuclear Energy,Spent Fuel and Waste Science and Technology Campaign,and by the US Department of Energy(DOE),the Office of Basic Energy Sciences,Chemical Sciences,Geosciences,and Biosciences Division both under Contract Number DE-AC02-05CH11231 with Lawrence Berkeley National Laboratorysupport from the US National Science Foundation(grant CMMI-2239630)funding from the European Research Council(ERC)under the European Union’s Horizon 2020 research and innovation programme(grant agreement No.101002507)the UK Natural Environment Research Council(NERC)for funding SeisGreen Project(Grant No.NE/W009293/1)which supported this workthe Royal Society UK for supporting this research through fellowship UF160443IMEDEA is an accredited"Maria de Maeztu Excellence Unit"(Grant CEX2021-001198,funded by MICIU/AEI/10.13039/501100011033).
文摘Coupled thermo-hydro-mechanical(THM)processes in fractured rock are playing a crucial role in geoscience and geoengineering applications.Diverse and conceptually distinct approaches have emerged over the past decades in both continuum and discontinuum perspectives leading to significant progress in their comprehending and modeling.This review paper offers an integrated perspective on existing modeling methodologies providing guidance for model selection based on the initial and boundary conditions.By comparing various models,one can better assess the uncertainties in predictions,particularly those related to the conceptual models.The review explores how these methodologies have significantlyenhanced the fundamental understanding of how fractures respond to fluid injection and production,and improved predictive capabilities pertaining to coupled processes within fractured systems.It emphasizes the importance of utilizing advanced computational technologies and thoroughly considering fundamental theories and principles established through past experimental evidence and practical experience.The selection and calibration of model parameters should be based on typical ranges and applied to the specificconditions of applications.The challenges arising from inherent heterogeneity and uncertainties,nonlinear THM coupled processes,scale dependence,and computational limitations in representing fieldscale fractures are discussed.Realizing potential advances on computational capacity calls for methodical conceptualization,mathematical modeling,selection of numerical solution strategies,implementation,and calibration to foster simulation outcomes that intricately reflectthe nuanced complexities of geological phenomena.Future research efforts should focus on innovative approaches to tackle the hurdles and advance the state-of-the-art in this critical fieldof study.
基金Supported by the National Natural Science Foundation of China(52374004)National Key R&D Program of China(2023YFF0614102,2023YFE0110900).
文摘Based on the finite element-discrete element numerical method,a numerical model of fracture propagation in deflagration fracturing was established by considering the impact of stress wave,quasi-static pressure of explosive gas,and reflection of stress wave.The model was validated against the results of physical experiments.Taking the shale reservoirs of Silurian Longmaxi Formation in Luzhou area of the Sichuan Basin as an example,the effects of in-situ stress difference,natural fracture parameters,branch wellbore spacing,delay detonation time,and angle between branch wellbore and main wellbore on fracture propagation were identified.The results show that the fracture propagation morphology in deflagration fracturing is less affected by the in-situ stress difference when it is 5-15 MPa,and the tendency of fracture intersection between branch wellbores is significantly weakened when the in-situ stress difference reaches 20 MPa.The increase of natural fracture length promotes the fracture propagation along the natural fracture direction,while the increase of volumetric natural fracture density and angle limits the fracture propagation area and reduces the probability of fracture intersection between branch wells.The larger the branch wellbore spacing,the less probability of the fracture intersection between branch wells,allowing for the fracture propagation in multiple directions.Increasing the delay detonation time decreases the fracture spacing between branch wellbores.When the angle between the branch wellbore and the main wellbore is 45°and 90°,there is a tendency of fracture intersection between branch wellbores.
基金the National Natural Science Foundation of China(No.52204052)the NationalNatural Science Foundation of China(No.U23B20156)the Sichuan Science and Technology Program(No.2023NSFSC0933).
文摘Shale reservoirs are characterized by numerous geological discontinuities,such as bedding planes,and exhibit pronounced heterogeneity across rock layers separated by these planes.Bedding planes often possess distinct mechanical properties compared to the surrounding rock matrix,particularly in terms of damage and fracture behavior.Consequently,vertical propagation of hydraulic fractures is influenced by both bedding planes and the heterogeneity.In this study,a numerical investigation into the height growth of hydraulic fractures was conducted using the finite element method,incorporating zero-thickness cohesive elements.The analysis explored the effects of bedding planes,toughness contrasts between layers,and variations in in-situ stress across different strata.The results reveal that hydraulic fractures are more likely to propagate along bedding planes instead of traversing them and extending vertically into barrier layers when(1)bedding strength is low,(2)stress contrast between layers is high,and(3)toughness contrast is significant.Furthermore,for a given bedding strength,increased stress contrast or higher toughness contrast between layers elevate hydraulic fracture extension pressure.When a substantial stress difference exists between layers(Lc 0.4),hydraulic=fractures preferentially propagate along bedding planes.Conversely,as bedding strength increases,the propagation distance along bedding planes decreases,accompanied by an amplified horizontal compressive stress field.Notably,when the stress difference is sufficiently small(SD a phenomenon termed“stress rolling”emerges,wherein<-0.2),hydraulic fractures deviate from vertical growth and instead extend along a near-horizontal trajectory.
基金Supported by the Basic Science Center Project of the National Natural Science Foundation of China(52288101).
文摘Based on continuum-discontinuum element method,the numerical simulation of fracture propagation during deflagration-hydraulic composite fracturing was constructed by considering deflagration stress impact induced fracture creation,deflagrating gas driven fracture propagation,and hydraulic fracture propagation,exploring the effects of in-situ stress difference,deflagration peak pressure,deflagration pressurization rate,hydraulic fracturing displacement and hydraulic fracturing fluid viscosity on fracture propagation in deflagration-hydraulic composite fracturing.The deflagration-hydraulic composite fracturing combines the advantages of deflagration fracturing in creating complex fractures near wells and the deep penetration of hydraulic fracturing at the far-field region,which can form multiple deep penetrating long fractures with better stimulation effects.With the increase of in-situ stress difference,the stimulated area of deflagration-hydraulic composite fracturing is reduced,and the deflagration-hydraulic composite fracturing is more suitable for reservoirs with small in-situ stress difference.Higher peak pressure and pressurization rate are conducive to increasing the maximum fracture length and burst degree of the deflagration fractures,which in turn increases the stimulated area of deflagration-hydraulic composite fracturing and improves the stimulation effect.Increasing the displacement and viscosity of hydraulic fracturing fluid can enhance the net pressure within the fractures,activate the deflagration fractures,increase the turning radius of the fractures,generate more long fractures,and effectively increase the stimulated reservoir area.The stimulated reservoir area is not completely positively correlated with the hydraulic fracturing displacement and fracturing fluid viscosity,and there is a critical value.When the critical value is exceeded,the stimulated area decreases.
基金Supported by Joint Funds of the National Natural Science Foundation of China(U23B6004).
文摘Stereoscopic particle image velocimetry technology was employed to investigate the planar three-dimensional velocity field and the process of proppant entry into branch fractures in a fracture configuration of“vertical main fracture-vertical branch fracture”intersecting at a 90°angle.This study analyzed the effects of pumping rate,fracturing fluid viscosity,proppant particle size,and fracture width on the transport behavior of proppant into branch fractures.Based on the deflection behavior of proppant,the main fractures can be divided into five regions:pre-entry transition,pre-entry stabilization,deflection entry at the fracture mouth,rear absorption entry,and movement away from the fracture mouth.Proppant primarily deflects into the branch fracture at the fracture mouth,with a small portion drawn in from the rear of the intersection.Increasing the pumping rate,reducing the proppant particle size,and widening the branch fracture are conducive to promoting proppant deflection into the branch.With increasing fracturing fluid viscosity,the ability of proppant to enter the branch fracture first improves and then declines,indicating that excessively high viscosity is unfavorable for proppant entry into the branch.During field operations,a high pumping rate and micro-to small-sized proppant can be used in the early stage to ensure effective placement in the branch fractures,followed by medium-to large-sized proppant to ensure adequate placement in the main fracture and enhance the overall conductivity of the fracture network.
基金support from the National Natural Science Foundation of China(Grant Nos.42320104003,42177175,and 42077247)the Fundamental Research Funds for the Central Universities.
文摘Accurate prediction of hydraulic fracture propagation is vital for Enhanced Geothermal System(EGS)design.We study the first hydraulic fracturing job at the GR1 well in the Gonghe Basin using field data,where the overall direction of hydraulic fractures does not show a delineated shape parallel to the maximum principal stress orientation.A field-scale numerical model based on the distinct element method is set up to carry out a fully coupled hydromechanical simulation,with the explicit representation of natural fractures via the discrete fracture network(DFN)approach.The effects of injection parameters and in situ stress on hydraulic fracture patterns are then quantitatively assessed.The study reveals that shear-induced deformation primarily governs the fracturing morphology in the GR1 well,driven by smaller injection rates and viscosities that promote massive activation of natural fractures,ultimately dominating the direction of hydraulic fracturing.Furthermore,the increase of in situ differential stress may promote shear damage of natural fracture surfaces,with the exact influence pattern depending on the combination of specific discontinuity properties and in situ stress state.Finally,we provide recommendations for EGS fracturing based on the influence characteristics of multiple parameters.This study can serve as an effective basis and reference for the design and optimization of EGS in the Gonghe basin and other sites.
基金supported by the National Natural Science Foundation of China(Grant No.12302500)the National Key Research and Development Program of China(Grant No.2020YFA0710503)Postdoctoral Fellowship Program(Grade B)of China Postdoctoral Science Foundation(Grant No.GBZ20230022).
文摘Controllable shock wave fracturing is an innovative engineering technique used for shale reservoir fracturing and reformation.Understanding the anisotropic fracture mechanism of shale under impact loading is vital for optimizing shock wave fracturing equipment and enhancing shale oil production.In this study,using the well-known notched semi-circular bend(NSCB)sample and the novel double-edge notched flattened Brazilian disc(DNFBD)sample combined with a split Hopkinson pressure bar(SHPB),various dynamic anisotropic fracture properties of Lushan shale,including failure characteristics,fracture toughness,energy dissipation and crack propagation velocity,are comprehensively compared and discussed under mode Ⅰ and mode Ⅱ fracture scenarios.First,using a newly modified fracture criterion considering the strength anisotropy of shale,the DNFBD specimen is predicted to be a robust method for true mode Ⅱ fracture of anisotropic shale rocks.Our experimental results show that the dynamic mode Ⅱ fracture of shale induces a rougher and more complex fracture morphology and performs a higher fracture toughness or fracture energy compared to dynamic mode Ⅰ fracture.The minimal fracture toughness or fracture energy occurs in the Short-transverse orientation,while the maximal ones occur in the Divider orientation.In addition,it is interesting to find that the mode Ⅱ fracture toughness anisotropy index decreases more slowly than that in the mode Ⅰ fracture scenario.These results provide significant insights for understanding the different dynamic fracture mechanisms of anisotropic shale rocks under impact loading and have some beneficial implications for the controllable shock wave fracturing technique.
基金supported by the National Natural Science Foun-dation of China(No.52434001).
文摘The linear elastic hydraulic fracture criterion is not applicable to deep reservoirs when nonlinear behavior is present over an extensive zone at the fracture tip.This study aims to develop a criterion for nonlinear hydraulic fracture considering the fracture process zone(FPZ)and seeks to reveal the causes of nonlinearity during fracture propagation in deep reservoirs.A closing stress profile considering the in-situ stress was established by using the cohesive zone model(CZM)to describe the FPZ at the fracture tip.An analytical model for the FPZ length was derived,while the criterion for nonlinear fracture propagation was proposed.The FPZ fully developed and the fracture began to propagate when the apparent stress intensity at the fracture tip reached the apparent fracture toughness or when the in-situ stress intensity reached the in-situ fracture toughness.The proposed criterion can clearly determine the length of the FPZ,accurately predict the breakdown pressure during fracturing operations,and establish a relationship between these two parameters.It addresses the inherent limitations of conventional linear elastic fracture mechanics(LEFM),which often underestimates fracture toughness and neglects the effects of the FPZ.This research is expected to enhance the fracturing design in deep reservoirs.
基金supported by the National Natural Science Foundation of China(Grant No.52404020 and Grant No.U1762216).
文摘Heterogeneity in carbonate formations due to discontinuities(e.g.,fractures and cavities)will bring about distinctive acid stimulation effects.However,the differences in fracturing behavior between homogeneous and heterogeneous carbonate formations remain unclear,complicating the optimization of acid fracturing strategies.In this paper,full-diameter carbonate rock samples with different degrees of discontinuity development are selected to investigate the fracturing behavior under different fluid types and injection schemes.Advanced techniques,including 3D CT scanning and 3D laser scanning,are employed to analyze fracture morphology and etching characteristics,respectively.Experimental results show that the coupled hydraulic-chemical effects play different roles in fracture induction between fracture-cavity developed and undeveloped carbonate rocks.Acid-fracturing stimulation consistently induces multiple types of complex fractures in fracture-cavity carbonate rocks,whereas it results in a single artificial fracture in less fracture-cavity carbonate rocks.Furthermore,localized etching patterns are prevalent in most fracture-cavity carbonate rocks,whereas homogeneous carbonate rocks exhibit regional or global etching characteristics.In both carbonate rocks,the stimulation effect of guar fluid is inferior to that of gelled acid but comparable to self-generating acid.Further findings are that alternating fracturing with guar and acid fluids in fracture-cavity carbonate rocks can sustain or even increase the injection pressure,facilitating the formation of new or depth-penetrating fractures.This phenomenon,however,is not observed in fracture-cavity undeveloped carbonate rocks.Potential interaction modes between induced fracture and natural fractures/cavities under different injection conditions are also identified.Finally,preferred fracturing schemes applicable to different carbonate formations are recommended based on the area,number and roughness of the induced fractures.
