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.展开更多
Multistage fracturing technology has been used to enhance tight hydrocarbon resource recovery.Determining the proper well spacing and fracturing strategy is crucial for generating a complex fracture network that facil...Multistage fracturing technology has been used to enhance tight hydrocarbon resource recovery.Determining the proper well spacing and fracturing strategy is crucial for generating a complex fracture network that facilitates oil and gas flow in reservoirs.The stress-shadow effect that occurs between multiple wells significantly affects the development of fracture networks in reservoirs.However,the quantification of the stress-shadow effect and its influence on fracture networks has not been satisfactorily resolved because of the difficulties in detecting and identifying fracture propagation and reorientation in reservoirs.In this study,based on the geological information from the Shengli oilfield,we applied a hybrid finite element-discrete element method to analyze engineering-scale three-dimensional fracture propagation and reorientation by altering well spacings and fracturing strategies.The results indicate that the fracturing area generated by the synchronous fracturing scheme is much smaller than those generated by the sequential and alternative schemes.An alternative hydrofracturing scheme is optimal with respect to fracturing area.The stress-blind area was defined to quantify the mechanical disturbance between adjacent wells.Our study improves the understanding of the effect of fracturing schemes on fracture networks and the impact of independent factors contributing to stress-shadow effects.展开更多
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.展开更多
Deep geological repository is typically situated at depths ranging from several hundred to 1000 m below ground,making bentonite engineered barrier potentially vulnerable to high water pressure and even inducing hydrau...Deep geological repository is typically situated at depths ranging from several hundred to 1000 m below ground,making bentonite engineered barrier potentially vulnerable to high water pressure and even inducing hydraulic fracturing.This study conducted injection tests on compacted GMZ(Gaomiaozi)bentonite with a self-developed visualization set-up.The objective was to unveil the roles of dry density,water content,and pressurization rate in hydraulic fracturing from the perspective of fracturing macromorphological dynamics and breakthrough characteristics.Moreover,the relationships between breakthrough characteristics and microstructure were examined by MIP(mercury intrusion porosimetry)analysis.Results showed that the fracturing dynamics were characterized by three stages:hydration,cracking,and fracturing stages.Compared to water content and pressurization rate,dry density exerted more pronounced effects on these stages.Increasing dry density can lead to an expansion of circular hydration zone,a more complex cracking network,and a change in fracturing patterns from long and clear to short and fuzzy.In terms of breakthrough characteristics,the breakthrough pressure was positively correlated with dry density and negatively correlated with water content.Interestingly,there is a good and unique logarithmic correlation between the breakthrough pressure and the ratio eM/em of inter-aggregate void ratio and intra-aggregate void ratio,regardless of dry density and water content.Within a certain range(i.e.200-50 kPa/min),breakthrough pressure showed slight dependency on pressurization rate.Nevertheless,an extremely low pressurization rate of 20 kPa/min caused a transition for the specimen from quasi-brittle to plastic state owning to more water infiltration,thereby hindering fracture initiation and propagation.展开更多
Medium-high maturity continental shale oil is one of the hydrocarbon resources with the most potential for successful development in China.Nevertheless,the unique geological conditions of a multi-lithologic superposit...Medium-high maturity continental shale oil is one of the hydrocarbon resources with the most potential for successful development in China.Nevertheless,the unique geological conditions of a multi-lithologic superposition shield the vertical propagation of hydraulic fractures and limit the longitudinal reconstruction in reservoirs,posing a great challenge for large-scale volumetric fracturing.Radial wellbore crosslayer fracturing,which transforms the interaction between the hydraulic fractures and lithologic interface into longitudinal multilayer competitive initiation,could provide a potential solution for this engineering challenge.To determine the longitudinal propagation behaviors of fractures guided by radial wellbores,true triaxial fracturing experiments were performed on multilayer shale-sandstone samples,with a focus on the injection pressure response,fracture morphology,and cross-layer pattern.The effects of the radial borehole length L,vertical stress difference K_(v),injection rate Q,and viscosity m of the fracturing fluid were analyzed.The results indicate that radial wellbores can greatly facilitate fracture initiation and cross-layer propagation.Unlike conventional hydraulic fracturing,there are two distinct fracture propagation patterns in radial wellbore fracturing:cross-layering and skip-layering.The fracture height guided by a radial wellbore is positively correlated with K_(v),Q,and m.Increasing these parameters causes a shift in the fracture initiation from a single root to an asynchronous root/toe end and can improve the cross-layer propagation capacity.Critical parameter thresholds exist for fracture propagation through and across interlayers under the guidance of radial boreholes.A parameter combination of critical cross-layering/skip-layering or alternating displacement/viscosity is recommended to simultaneously improve the fracture height and degree of lateral activation.The degree of correlation of different parameters with the vertical fracture height can be written as L>Q/m>K_(v).Increasing the radial wellbore length can effectively facilitate fracture cross-/skip-layer propagation and reduce the critical threshold of injection parameters,which is conducive to maximizing the stimulated reservoir volume.展开更多
This paper introduces a novel approach combining radial borehole fracturing with Water-Alternating-Gas(WAG)injection,enabling simultaneous WAG injection and shale oil production in a single vertical well.A numerical r...This paper introduces a novel approach combining radial borehole fracturing with Water-Alternating-Gas(WAG)injection,enabling simultaneous WAG injection and shale oil production in a single vertical well.A numerical reservoir model incorporating the modified exponential non-Darcy law,stress sensitivity,and diffusion is established.The spatial distribution of permeability reduction shows that stress sensitivity enhances the non-Darcy effect,with apparent permeability decreasing to 0-92.1%of the initial value,highlighting the importance of maintaining reservoir pressure.Continuous CO_(2) flooding leads to early gas breakthrough,while continuous water flooding has less displacement efficiency.A 30%water-to-gas injection time ratio improves oil production and delays gas breakthrough compared to continuous CO_(2) injection.Optimal conditions for effective recovery are identified as an initial production period of 100 d and a well vertical spacing of 30 m.This study compares the production capacity of WAG operations under radial borehole fracturing and horizontal well fracturing.When the number of wells is two for both cases,the production capacity of radial borehole fracturing is comparable to that of five-stage horizontal well fracturing,indicating that radial borehole fracturing can serve as an alternative or supplement to horizontal well fracturing when the reservoir volume is limited.This study offers a new method and theoretical basis for the efficient development of shale oil.展开更多
Two sets of alloys,Mg-Zn-Ca-xNi(0≤x≤5),have been developed with tunable corrosion and mechanical properties,optimized for fracturing materials.High-zinc artificial aged(T6)Mg-12Zn-0.5Ca-x Ni(0≤x≤5)series,featuring...Two sets of alloys,Mg-Zn-Ca-xNi(0≤x≤5),have been developed with tunable corrosion and mechanical properties,optimized for fracturing materials.High-zinc artificial aged(T6)Mg-12Zn-0.5Ca-x Ni(0≤x≤5)series,featuring a straightforward preparation method and the potential for manufacturing large-scale components,exhibit notable corrosion rates up to 29 mg cm^(-2)h^(-1)at 25℃ and 643 mg cm^(-2)h^(-1)at 93℃.The high corrosion rate is primary due to the Ni–containing second phases,which intensify the galvanic corrosion that overwhelms their corrosion barrier effect.Low-zinc rolled Mg-1.5Zn-0.2Ca-x Ni(0≤x≤5)series,characterizing excellent deformability with an elongation to failure of~26%,present accelerated corrosion rates up to 34 mg cm^(-2)h^(-1)at 25℃ and 942 mg cm^(-2)h^(-1)at 93℃.The elimination of corrosion barrier effect via deformation contributes to the further increase of corrosion rate compared to the T6 series.Additionally,Mg-Zn-Ca-xNi(0≤x≤5)alloys exhibit tunable ultimate tensile strengths ranging from~190 to~237 MPa,depending on their specific composition.The adjustable corrosion rate and mechanical properties render the Mg-Zn-Ca-x Ni(0≤x≤5)alloys suitable for fracturing materials.展开更多
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.展开更多
Through a case analysis,this study examines the spatiotemporal evolution of microseismic(MS)events,energy characteristics,volumetric features,and fracture network development in surface well hydraulic fracturing.A tot...Through a case analysis,this study examines the spatiotemporal evolution of microseismic(MS)events,energy characteristics,volumetric features,and fracture network development in surface well hydraulic fracturing.A total of 349 MS events were analyzed across different fracturing sections,revealing significant heterogeneity in fracture propagation.Energy scanning results showed that cumulative energy values ranged from 240 to 1060 J across the sections,indicating notable differences.Stimulated reservoir volume(SRV)analysis demonstrated well-developed fracture networks in certain sections,with a total SRV exceeding 1540000 m^(3).The hydraulic fracture network analysis revealed that during the midfracturing stage,the density and spatial extent of MS events significantly increased,indicating rapid fracture propagation and the formation of complex networks.In the later stage,the number of secondary fractures near fracture edges decreased,and the fracture network stabilized.By comparing the branching index,fracture length,width,height,and SRV values across different fracturing sections,Sections No.1 and No.8 showed the best performance,with high MS event densities,extensive fracture networks,and significant energy release.However,Sections No.4 and No.5 exhibited sparse MS activity and poor fracture connectivity,indicating suboptimal stimulation effectiveness.展开更多
To more accurately describe the coal damage and fracture evolution law during liquid nitrogen(LN_(2))fracturing under true triaxial stress,a thermal-hydraulic-mechanical-damage(THMD)coupling model for LN_(2) fracturin...To more accurately describe the coal damage and fracture evolution law during liquid nitrogen(LN_(2))fracturing under true triaxial stress,a thermal-hydraulic-mechanical-damage(THMD)coupling model for LN_(2) fracturing coal was developed,considering the coal heterogeneity and thermophysical parameters of nitrogen.The accuracy and applicability of model were verified by comparing with LN_(2) injection pre-cooling and fracturing experimental data.The effects of different pre-cooling times and horizontal stress ratios on coal damage evolution,permeability,temperature distribution,and fracture characteristics were analyzed.The results show that the permeability and damage of the coal increase exponentially,while the temperature decreases exponentially during the fracturing process.As the pre-cooling time increases,the damage range of the coal expands,and the fracture propagation becomes more pronounced.The initiation pressure and rupture pressure decrease and tend to stabilize with longer precooling times.As the horizontal stress ratio increases,fractures preferentially extend along the direction of maximum horizontal principal stress,leading to a significant decrease in both initiation and rupture pressures.At a horizontal stress ratio of 3,the initiation pressure drops by 48.07%,and the rupture pressure decreases by 41.36%.The results provide a theoretical basis for optimizing LN_(2) fracturing techniques and improving coal seam modification.展开更多
Discrete fracture network(DFN)commonly existing in natural rock masses plays an important role in geological complexity which can influence rock fracturing behaviour during fluid injection.This paper simulated the hyd...Discrete fracture network(DFN)commonly existing in natural rock masses plays an important role in geological complexity which can influence rock fracturing behaviour during fluid injection.This paper simulated the hydraulic fracturing process in lab-scale coal samples with DFNs and the induced seismic activities by the discrete element method(DEM).The effects of DFNs on hydraulic fracturing,induced seismicity and elastic property changes have been concluded.Denser DFNs can comprehensively decrease the peak injection pressure and injection duration.The proportion of strong seismic events increases first and then decreases with increasing DFN density.In addition,the relative modulus of the rock mass is derived innovatively from breakdown pressure,breakdown fracture length and the related initiation time.Increasing DFN densities among large(35–60 degrees)and small(0–30 degrees)fracture dip angles show opposite evolution trends in relative modulus.The transitional point(dip angle)for the opposite trends is also proportionally affected by the friction angle of the rock mass.The modelling results have much practical meaning to infer the density and geometry of pre-existing fractures and the elastic property of rock mass in the field,simply based on the hydraulic fracturing and induced seismicity monitoring data.展开更多
Supercritical CO_(2)(SC-CO_(2))fracturing stands out a promising waterless stimulation technique in the development of unconventional resources.While numerous studies have delved into the inducedfracture mechanism of ...Supercritical CO_(2)(SC-CO_(2))fracturing stands out a promising waterless stimulation technique in the development of unconventional resources.While numerous studies have delved into the inducedfracture mechanism of SC-CO_(2),the small scale of rock samples and synthetic materials used in many studies have limited a comprehensive understanding of fracture propagation in unconventional formations.In this study,cubic tight sandstone samples with dimensions of 300 mm were employed to conduct SC-CO_(2)fractu ring experiments under true-triaxial stre ss conditions.The spatial morphology and quantitative attributes of fracture induced by water and SC-CO_(2)fracturing were compared,while the impact of in-situ stress on fracture propagation was also investigated.The results indicate that the SCCO_(2)fracturing takes approximately ten times longer than water fracturing.Furthermore,under identical stress condition,the breakdown pressure(BP)for SC-CO_(2)fracturing is nearly 25%lower than that for water fracturing.A quantitative analysis of fracture morphology reveals that water fracturing typically produces relatively simple fracture pattern,with the primary fracture distribution predominantly controlled by bedding planes.In contrast,SC-CO_(2)fracturing results in a more complex fracture morphology.As the differential of horizontal principal stress increases,the BP for SC-CO_(2)fractured rock exhibits a downward trend,and the induced fracture morphology becomes more simplified.Moreover,the presence of abnormal in-situ stress leads to a further increase in the BP for SC-CO_(2)fracturing,simultaneously enhancing the development of a more conductive fracture network.These findings provide critical insights into the efficiency and behavior of SC-CO_(2)fracturing in comparison to traditional water-based fracturing,offering valuable implication for its potential applications in unconventional reservoirs.展开更多
Source properties and stress fields are critical to understand fundamental mechanisms for fluid-induced earthquakes.In this study,we identify the focal mechanism solutions(FMSs)of 360 earthquakes with local magnitude ...Source properties and stress fields are critical to understand fundamental mechanisms for fluid-induced earthquakes.In this study,we identify the focal mechanism solutions(FMSs)of 360 earthquakes with local magnitude M_(L)≥1.5 in the Changning shale gas field from January 2016 to May 2017 by fitting three-component waveforms.We then constrain the directions of the maximum horizontal stress(σ_(H_(max)))for four dense earthquake clusters using the stress tensor inversion method.The stress drops of 121 earthquakes with M_(L)≥1.5 are calculated using the spectral ratio method.We examine the spatiotemporal heterogeneity of stress field,and discuss the cause of non-double-couple(non-DC)components in seismicity clusters.Following the Mohr-Coulomb criterion,we estimate the fluid overpressure thresholds from FMS for different seismic clusters,providing insights into potential physical mechanisms for induced seismicity.The FMS results indicate that shallow reverse earthquakes,with steep dip angles,characterize most events.The source mechanisms of earthquakes with M_(L)≥1.5 are dominated by DC components(>70%),but several earthquakes with M_(L)>3.0 and the microseismic events nearby during injection period display significant non-DC components(>30%).Stress inversion results reveal that the σ_(H_(max)) direction ranges from 120°to 128°.Stress drops of earthquakes range between 0.10 and 64.49 MPa,with high values occurring on reverse faults situated at a greater distance from the shale layer,accompanied by a moderate rotation(≤25°)in the trend of σ_(H_(max)).The seismic clusters close to the shale layer exhibit low fluid overpressure thresholds,prone to being triggered by high pore-pressure fluid.The integrated results suggest that the diffusion of high pore pressures is likely to be the primary factor for observed earthquakes.The present results are expected to offer valuable insights into the origin of anomalous seismicity near the shale gas sites.展开更多
A coupled thermal-hydro-mechanical cohesive phase-field model for hydraulic fracturing in deep coal seams is presented.Heat exchange between the cold fluid and the hot rock is considered,and the thermal contribution t...A coupled thermal-hydro-mechanical cohesive phase-field model for hydraulic fracturing in deep coal seams is presented.Heat exchange between the cold fluid and the hot rock is considered,and the thermal contribution terms between the cold fluid and the hot rock are derived.Heat transfer obeys Fourier's law,and porosity is used to relate the thermodynamic parameters of the fracture and matrix domains.The net pressure difference between the fracture and the matrix is neglected,and thus the fluid flow is modeled by the unified fluid-governing equations.The evolution equations of porosity and Biot's coefficient during hydraulic fracturing are derived from their definitions.The effect of coal cleats is considered and modeled by Voronoi polygons,and this approach is shown to have high accuracy.The accuracy of the proposed model is verified by two sets of fracturing experiments in multilayer coal seams.Subsequently,the differences in fracture morphology,fluid pressure response,and fluid pressure distribution between direct fracturing of coal seams and indirect fracturing of shale interlayers are explored,and the effects of the cluster number and cluster spacing on fracture morphology for multi-cluster fracturing are also examined.The numerical results show that the proposed model is expected to be a powerful tool for the fracturing design and optimization of deep coalbed methane.展开更多
Aiming at mitigating the high risks associated with conventional explosive blasting,this study developed a safe directional fracturing technique,i.e.instantaneous expansion with a single fracture(IESF),using a coal-ba...Aiming at mitigating the high risks associated with conventional explosive blasting,this study developed a safe directional fracturing technique,i.e.instantaneous expansion with a single fracture(IESF),using a coal-based solid waste expanding agent.First,the mechanism of directional fracturing blasting by the IESF was analyzed,and the criterion of directional crack initiation was established.On this basis,laboratory experiments and numerical simulations were conducted to systematically evaluate the directional fracturing blasting performance of the IESF.The results indicate that the IESF presents an excellent directional fracturing effect,with average surface undulation differences ranging from 8.1 mm to 22.7 mm on the fracture surfaces.Moreover,during concrete fracturing tests,the stresses and strains in the fracturing direction are measured to be 2.16-3.71 times and 8 times larger than those in the nonfracturing direction,respectively.Finally,the IESF technique was implemented for no-pillar mining with gob-side entry retaining through roof cutting and pressure relief in an underground coal mine.The IESF technique effectively created directional cracks in the roof without causing severe roadway deformation,achieving an average cutting rate and maximum roadway deformation of 94%and 197 mm,respectively.These on-site test results verified its excellent directional rock fracturing performance.The IESF technique,which is safe,efficient,and green,has considerable application prospects in the field of rock mechanics and engineering.展开更多
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.展开更多
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.展开更多
To accurately analyze proppant transport in rough intersecting fractures and elucidate the interaction mechanisms among liquid,particles,and rough walls,this study reconstructed a numerical model of fractures in inhom...To accurately analyze proppant transport in rough intersecting fractures and elucidate the interaction mechanisms among liquid,particles,and rough walls,this study reconstructed a numerical model of fractures in inhomogeneous reservoirs with varying brittleness index(BI).Various auto-correlation Gaussian rough fracture models were created using Matlab to assess roughness through the fractal dimension method.This research innovatively combined Boolean operations to establish three-dimensional rough fracture models,incorporating(Computational Fluid Dynamics)CFD-DEM(Discrete Element Method)with a bidirectional method for cosimulation.The proppant transport in fractures was categorized into three zones based on the difference in the turbulent kinetic energy.Artificially induced fracture roughness increases fluid retention and turbulence,causing plugging effects and limiting proppant flow into branch fractures.Additionally,compared with the superior deposition and significant support effects of the spherical proppant,the low-sphericity proppant traveled farther under fracturing fluid,inducing more pronounced plugging near curved fracture intersections;the variation in fracture intersection angles primarily impacted the wall shear stress within the flow field,indicating smaller angles led to higher shear energy at the intersection.Compared with the intersection angle of 30°,the height and area deposited in the 90 branch fracture increased by 52.25%and 65.33%,respectively:notably,injecting proppant from smaller to larger particles(S:M:L)and a low velocity effectively ensured fracture conductivity near the wellbore at joint roughness coefficient(JRC)≥46 while achieving satis-factory placement in the branch fracture,making it a recommended approach.展开更多
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.展开更多
The Carter model is used to characterize the dynamic behaviors of fracture growth and fracturing fluid leakoff.A thermo-fluid coupling temperature response forward model is built considering the fluid flow and heat tr...The Carter model is used to characterize the dynamic behaviors of fracture growth and fracturing fluid leakoff.A thermo-fluid coupling temperature response forward model is built considering the fluid flow and heat transfer in wellbore,fracture and reservoir.The influences of fracturing parameters and fracture parameters on the responses of distributed temperature sensing(DTS)are analyzed,and a diagnosis method of fracture parameters is presented based on the simulated annealing algorithm.A field case study is introduced to verify the model’s reliability.Typical V-shaped characteristics can be observed from the DTS responses in the multi-cluster fracturing process,with locations corresponding to the hydraulic fractures.The V-shape depth is shallower for a higher injection rate and longer fracturing and shut-in time.Also,the V-shape is wider for a higher fracture-surface leakoff coefficient,longer fracturing time and smaller fracture width.Additionally,the cooling effect near the wellbore continues to spread into the reservoir during the shut-in period,causing the DTS temperature to decrease instead of rise.Real-time monitoring and interpretation of DTS temperature data can help understand the fracture propagation during fracturing operation,so that immediate measures can be taken to improve the fracturing performance.展开更多
文摘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.
基金supported in part by the National Key Research and Development Project of China(No.2022YFC3004602)in part by the National Natural Science Foundation of China(Nos.52121003 and 52342403).
文摘Multistage fracturing technology has been used to enhance tight hydrocarbon resource recovery.Determining the proper well spacing and fracturing strategy is crucial for generating a complex fracture network that facilitates oil and gas flow in reservoirs.The stress-shadow effect that occurs between multiple wells significantly affects the development of fracture networks in reservoirs.However,the quantification of the stress-shadow effect and its influence on fracture networks has not been satisfactorily resolved because of the difficulties in detecting and identifying fracture propagation and reorientation in reservoirs.In this study,based on the geological information from the Shengli oilfield,we applied a hybrid finite element-discrete element method to analyze engineering-scale three-dimensional fracture propagation and reorientation by altering well spacings and fracturing strategies.The results indicate that the fracturing area generated by the synchronous fracturing scheme is much smaller than those generated by the sequential and alternative schemes.An alternative hydrofracturing scheme is optimal with respect to fracturing area.The stress-blind area was defined to quantify the mechanical disturbance between adjacent wells.Our study improves the understanding of the effect of fracturing schemes on fracture networks and the impact of independent factors contributing to stress-shadow effects.
基金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 the National Natural Science Foundation of China(Grant Nos.42430713 and 42125701)Innovation Program of Shanghai Municipal Education Commission(Grant No.2023ZKZD26)。
文摘Deep geological repository is typically situated at depths ranging from several hundred to 1000 m below ground,making bentonite engineered barrier potentially vulnerable to high water pressure and even inducing hydraulic fracturing.This study conducted injection tests on compacted GMZ(Gaomiaozi)bentonite with a self-developed visualization set-up.The objective was to unveil the roles of dry density,water content,and pressurization rate in hydraulic fracturing from the perspective of fracturing macromorphological dynamics and breakthrough characteristics.Moreover,the relationships between breakthrough characteristics and microstructure were examined by MIP(mercury intrusion porosimetry)analysis.Results showed that the fracturing dynamics were characterized by three stages:hydration,cracking,and fracturing stages.Compared to water content and pressurization rate,dry density exerted more pronounced effects on these stages.Increasing dry density can lead to an expansion of circular hydration zone,a more complex cracking network,and a change in fracturing patterns from long and clear to short and fuzzy.In terms of breakthrough characteristics,the breakthrough pressure was positively correlated with dry density and negatively correlated with water content.Interestingly,there is a good and unique logarithmic correlation between the breakthrough pressure and the ratio eM/em of inter-aggregate void ratio and intra-aggregate void ratio,regardless of dry density and water content.Within a certain range(i.e.200-50 kPa/min),breakthrough pressure showed slight dependency on pressurization rate.Nevertheless,an extremely low pressurization rate of 20 kPa/min caused a transition for the specimen from quasi-brittle to plastic state owning to more water infiltration,thereby hindering fracture initiation and propagation.
基金supported by the National Natural Science Foun-dation of China(52421002,U24B6001,52204019,and 52192624)the Open Foundation of the Shanxi Key Laboratory of Carbon Dioxide Sequestration and Enhanced Oil Recovery。
文摘Medium-high maturity continental shale oil is one of the hydrocarbon resources with the most potential for successful development in China.Nevertheless,the unique geological conditions of a multi-lithologic superposition shield the vertical propagation of hydraulic fractures and limit the longitudinal reconstruction in reservoirs,posing a great challenge for large-scale volumetric fracturing.Radial wellbore crosslayer fracturing,which transforms the interaction between the hydraulic fractures and lithologic interface into longitudinal multilayer competitive initiation,could provide a potential solution for this engineering challenge.To determine the longitudinal propagation behaviors of fractures guided by radial wellbores,true triaxial fracturing experiments were performed on multilayer shale-sandstone samples,with a focus on the injection pressure response,fracture morphology,and cross-layer pattern.The effects of the radial borehole length L,vertical stress difference K_(v),injection rate Q,and viscosity m of the fracturing fluid were analyzed.The results indicate that radial wellbores can greatly facilitate fracture initiation and cross-layer propagation.Unlike conventional hydraulic fracturing,there are two distinct fracture propagation patterns in radial wellbore fracturing:cross-layering and skip-layering.The fracture height guided by a radial wellbore is positively correlated with K_(v),Q,and m.Increasing these parameters causes a shift in the fracture initiation from a single root to an asynchronous root/toe end and can improve the cross-layer propagation capacity.Critical parameter thresholds exist for fracture propagation through and across interlayers under the guidance of radial boreholes.A parameter combination of critical cross-layering/skip-layering or alternating displacement/viscosity is recommended to simultaneously improve the fracture height and degree of lateral activation.The degree of correlation of different parameters with the vertical fracture height can be written as L>Q/m>K_(v).Increasing the radial wellbore length can effectively facilitate fracture cross-/skip-layer propagation and reduce the critical threshold of injection parameters,which is conducive to maximizing the stimulated reservoir volume.
基金the Young Scientists Fund of the National Natural Science Foundation of China(52204063)the Key Laboratory of Shale Gas Exploration,Ministry of Natural Resources(Chongqing Institute of Geology and Mineral Resources),Chongqing,China(KLSGE-202202).
文摘This paper introduces a novel approach combining radial borehole fracturing with Water-Alternating-Gas(WAG)injection,enabling simultaneous WAG injection and shale oil production in a single vertical well.A numerical reservoir model incorporating the modified exponential non-Darcy law,stress sensitivity,and diffusion is established.The spatial distribution of permeability reduction shows that stress sensitivity enhances the non-Darcy effect,with apparent permeability decreasing to 0-92.1%of the initial value,highlighting the importance of maintaining reservoir pressure.Continuous CO_(2) flooding leads to early gas breakthrough,while continuous water flooding has less displacement efficiency.A 30%water-to-gas injection time ratio improves oil production and delays gas breakthrough compared to continuous CO_(2) injection.Optimal conditions for effective recovery are identified as an initial production period of 100 d and a well vertical spacing of 30 m.This study compares the production capacity of WAG operations under radial borehole fracturing and horizontal well fracturing.When the number of wells is two for both cases,the production capacity of radial borehole fracturing is comparable to that of five-stage horizontal well fracturing,indicating that radial borehole fracturing can serve as an alternative or supplement to horizontal well fracturing when the reservoir volume is limited.This study offers a new method and theoretical basis for the efficient development of shale oil.
基金supported by the National Key Research and Development Program(No.2022YFE0122000)National Natural Science Foundation of China under Grant Nos.52234009,52274383,52222409,and 52201113。
文摘Two sets of alloys,Mg-Zn-Ca-xNi(0≤x≤5),have been developed with tunable corrosion and mechanical properties,optimized for fracturing materials.High-zinc artificial aged(T6)Mg-12Zn-0.5Ca-x Ni(0≤x≤5)series,featuring a straightforward preparation method and the potential for manufacturing large-scale components,exhibit notable corrosion rates up to 29 mg cm^(-2)h^(-1)at 25℃ and 643 mg cm^(-2)h^(-1)at 93℃.The high corrosion rate is primary due to the Ni–containing second phases,which intensify the galvanic corrosion that overwhelms their corrosion barrier effect.Low-zinc rolled Mg-1.5Zn-0.2Ca-x Ni(0≤x≤5)series,characterizing excellent deformability with an elongation to failure of~26%,present accelerated corrosion rates up to 34 mg cm^(-2)h^(-1)at 25℃ and 942 mg cm^(-2)h^(-1)at 93℃.The elimination of corrosion barrier effect via deformation contributes to the further increase of corrosion rate compared to the T6 series.Additionally,Mg-Zn-Ca-xNi(0≤x≤5)alloys exhibit tunable ultimate tensile strengths ranging from~190 to~237 MPa,depending on their specific composition.The adjustable corrosion rate and mechanical properties render the Mg-Zn-Ca-x Ni(0≤x≤5)alloys suitable for fracturing materials.
基金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 Yunlong Lake Laboratory of Deep Underground Science and Engineering Project(No.104024008)the National Natural Science Foundation of China(Nos.52274241 and 52474261)the Natural Science Foundation of Jiangsu Province(No.BK20240207).
文摘Through a case analysis,this study examines the spatiotemporal evolution of microseismic(MS)events,energy characteristics,volumetric features,and fracture network development in surface well hydraulic fracturing.A total of 349 MS events were analyzed across different fracturing sections,revealing significant heterogeneity in fracture propagation.Energy scanning results showed that cumulative energy values ranged from 240 to 1060 J across the sections,indicating notable differences.Stimulated reservoir volume(SRV)analysis demonstrated well-developed fracture networks in certain sections,with a total SRV exceeding 1540000 m^(3).The hydraulic fracture network analysis revealed that during the midfracturing stage,the density and spatial extent of MS events significantly increased,indicating rapid fracture propagation and the formation of complex networks.In the later stage,the number of secondary fractures near fracture edges decreased,and the fracture network stabilized.By comparing the branching index,fracture length,width,height,and SRV values across different fracturing sections,Sections No.1 and No.8 showed the best performance,with high MS event densities,extensive fracture networks,and significant energy release.However,Sections No.4 and No.5 exhibited sparse MS activity and poor fracture connectivity,indicating suboptimal stimulation effectiveness.
基金financially supported by the National Natural Science Foundation of China(Nos.51874236 and 52174207)Shaanxi Science and Technology Innovation Team(No.2022TD02)Henan University of Science and Technology PhD Funded Projects(No.B2025-9)。
文摘To more accurately describe the coal damage and fracture evolution law during liquid nitrogen(LN_(2))fracturing under true triaxial stress,a thermal-hydraulic-mechanical-damage(THMD)coupling model for LN_(2) fracturing coal was developed,considering the coal heterogeneity and thermophysical parameters of nitrogen.The accuracy and applicability of model were verified by comparing with LN_(2) injection pre-cooling and fracturing experimental data.The effects of different pre-cooling times and horizontal stress ratios on coal damage evolution,permeability,temperature distribution,and fracture characteristics were analyzed.The results show that the permeability and damage of the coal increase exponentially,while the temperature decreases exponentially during the fracturing process.As the pre-cooling time increases,the damage range of the coal expands,and the fracture propagation becomes more pronounced.The initiation pressure and rupture pressure decrease and tend to stabilize with longer precooling times.As the horizontal stress ratio increases,fractures preferentially extend along the direction of maximum horizontal principal stress,leading to a significant decrease in both initiation and rupture pressures.At a horizontal stress ratio of 3,the initiation pressure drops by 48.07%,and the rupture pressure decreases by 41.36%.The results provide a theoretical basis for optimizing LN_(2) fracturing techniques and improving coal seam modification.
基金Australian Research Council Linkage Program(LP200301404)for sponsoring this researchthe financial support provided by the State Key Laboratory of Geohazard Prevention and Geoenvironment Protection(Chengdu University of Technology,SKLGP2021K002)National Natural Science Foundation of China(52374101,32111530138).
文摘Discrete fracture network(DFN)commonly existing in natural rock masses plays an important role in geological complexity which can influence rock fracturing behaviour during fluid injection.This paper simulated the hydraulic fracturing process in lab-scale coal samples with DFNs and the induced seismic activities by the discrete element method(DEM).The effects of DFNs on hydraulic fracturing,induced seismicity and elastic property changes have been concluded.Denser DFNs can comprehensively decrease the peak injection pressure and injection duration.The proportion of strong seismic events increases first and then decreases with increasing DFN density.In addition,the relative modulus of the rock mass is derived innovatively from breakdown pressure,breakdown fracture length and the related initiation time.Increasing DFN densities among large(35–60 degrees)and small(0–30 degrees)fracture dip angles show opposite evolution trends in relative modulus.The transitional point(dip angle)for the opposite trends is also proportionally affected by the friction angle of the rock mass.The modelling results have much practical meaning to infer the density and geometry of pre-existing fractures and the elastic property of rock mass in the field,simply based on the hydraulic fracturing and induced seismicity monitoring data.
基金funded by the National Natural Scientific Foundation of China(Nos.52304008,52404038,52474043)the China Postdoctoral Science Foundation(No.2023MD734223)+1 种基金the Key Laboratory of Well Stability and Fluid&Rock Mechanics in Oil and Gas Reservoir of Shaanxi Province(No.23JS047)the Youth Talent Lifting Program of Xi'an Science and Technology Association(No.959202413078)。
文摘Supercritical CO_(2)(SC-CO_(2))fracturing stands out a promising waterless stimulation technique in the development of unconventional resources.While numerous studies have delved into the inducedfracture mechanism of SC-CO_(2),the small scale of rock samples and synthetic materials used in many studies have limited a comprehensive understanding of fracture propagation in unconventional formations.In this study,cubic tight sandstone samples with dimensions of 300 mm were employed to conduct SC-CO_(2)fractu ring experiments under true-triaxial stre ss conditions.The spatial morphology and quantitative attributes of fracture induced by water and SC-CO_(2)fracturing were compared,while the impact of in-situ stress on fracture propagation was also investigated.The results indicate that the SCCO_(2)fracturing takes approximately ten times longer than water fracturing.Furthermore,under identical stress condition,the breakdown pressure(BP)for SC-CO_(2)fracturing is nearly 25%lower than that for water fracturing.A quantitative analysis of fracture morphology reveals that water fracturing typically produces relatively simple fracture pattern,with the primary fracture distribution predominantly controlled by bedding planes.In contrast,SC-CO_(2)fracturing results in a more complex fracture morphology.As the differential of horizontal principal stress increases,the BP for SC-CO_(2)fractured rock exhibits a downward trend,and the induced fracture morphology becomes more simplified.Moreover,the presence of abnormal in-situ stress leads to a further increase in the BP for SC-CO_(2)fracturing,simultaneously enhancing the development of a more conductive fracture network.These findings provide critical insights into the efficiency and behavior of SC-CO_(2)fracturing in comparison to traditional water-based fracturing,offering valuable implication for its potential applications in unconventional reservoirs.
基金financially supported by the National Natural Science Foundation of China(Grant Nos.U20A20266 and 12302503)Scientific and technological research projects in Sichuan province(Grant No.2024NSFSC0973).
文摘Source properties and stress fields are critical to understand fundamental mechanisms for fluid-induced earthquakes.In this study,we identify the focal mechanism solutions(FMSs)of 360 earthquakes with local magnitude M_(L)≥1.5 in the Changning shale gas field from January 2016 to May 2017 by fitting three-component waveforms.We then constrain the directions of the maximum horizontal stress(σ_(H_(max)))for four dense earthquake clusters using the stress tensor inversion method.The stress drops of 121 earthquakes with M_(L)≥1.5 are calculated using the spectral ratio method.We examine the spatiotemporal heterogeneity of stress field,and discuss the cause of non-double-couple(non-DC)components in seismicity clusters.Following the Mohr-Coulomb criterion,we estimate the fluid overpressure thresholds from FMS for different seismic clusters,providing insights into potential physical mechanisms for induced seismicity.The FMS results indicate that shallow reverse earthquakes,with steep dip angles,characterize most events.The source mechanisms of earthquakes with M_(L)≥1.5 are dominated by DC components(>70%),but several earthquakes with M_(L)>3.0 and the microseismic events nearby during injection period display significant non-DC components(>30%).Stress inversion results reveal that the σ_(H_(max)) direction ranges from 120°to 128°.Stress drops of earthquakes range between 0.10 and 64.49 MPa,with high values occurring on reverse faults situated at a greater distance from the shale layer,accompanied by a moderate rotation(≤25°)in the trend of σ_(H_(max)).The seismic clusters close to the shale layer exhibit low fluid overpressure thresholds,prone to being triggered by high pore-pressure fluid.The integrated results suggest that the diffusion of high pore pressures is likely to be the primary factor for observed earthquakes.The present results are expected to offer valuable insights into the origin of anomalous seismicity near the shale gas sites.
基金Project supported by the National Natural Science Foundation of China(No.42202314)。
文摘A coupled thermal-hydro-mechanical cohesive phase-field model for hydraulic fracturing in deep coal seams is presented.Heat exchange between the cold fluid and the hot rock is considered,and the thermal contribution terms between the cold fluid and the hot rock are derived.Heat transfer obeys Fourier's law,and porosity is used to relate the thermodynamic parameters of the fracture and matrix domains.The net pressure difference between the fracture and the matrix is neglected,and thus the fluid flow is modeled by the unified fluid-governing equations.The evolution equations of porosity and Biot's coefficient during hydraulic fracturing are derived from their definitions.The effect of coal cleats is considered and modeled by Voronoi polygons,and this approach is shown to have high accuracy.The accuracy of the proposed model is verified by two sets of fracturing experiments in multilayer coal seams.Subsequently,the differences in fracture morphology,fluid pressure response,and fluid pressure distribution between direct fracturing of coal seams and indirect fracturing of shale interlayers are explored,and the effects of the cluster number and cluster spacing on fracture morphology for multi-cluster fracturing are also examined.The numerical results show that the proposed model is expected to be a powerful tool for the fracturing design and optimization of deep coalbed methane.
基金supported by the National Natural Science Foundation of China(Grant No.52404155)State Key Laboratory of Mining Disaster Prevention and Control(Shandong University of Science and Technology)+1 种基金Ministry of Education(Grant No.JMDPC202402)supported by the opening project of State Key Laboratory of Explosion Science and Safety Protection(Beijing Institute of Technology).The opening project number is KFJJ24-20M.
文摘Aiming at mitigating the high risks associated with conventional explosive blasting,this study developed a safe directional fracturing technique,i.e.instantaneous expansion with a single fracture(IESF),using a coal-based solid waste expanding agent.First,the mechanism of directional fracturing blasting by the IESF was analyzed,and the criterion of directional crack initiation was established.On this basis,laboratory experiments and numerical simulations were conducted to systematically evaluate the directional fracturing blasting performance of the IESF.The results indicate that the IESF presents an excellent directional fracturing effect,with average surface undulation differences ranging from 8.1 mm to 22.7 mm on the fracture surfaces.Moreover,during concrete fracturing tests,the stresses and strains in the fracturing direction are measured to be 2.16-3.71 times and 8 times larger than those in the nonfracturing direction,respectively.Finally,the IESF technique was implemented for no-pillar mining with gob-side entry retaining through roof cutting and pressure relief in an underground coal mine.The IESF technique effectively created directional cracks in the roof without causing severe roadway deformation,achieving an average cutting rate and maximum roadway deformation of 94%and 197 mm,respectively.These on-site test results verified its excellent directional rock fracturing performance.The IESF technique,which is safe,efficient,and green,has considerable application prospects in the field of rock mechanics and engineering.
基金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.
文摘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.
文摘To accurately analyze proppant transport in rough intersecting fractures and elucidate the interaction mechanisms among liquid,particles,and rough walls,this study reconstructed a numerical model of fractures in inhomogeneous reservoirs with varying brittleness index(BI).Various auto-correlation Gaussian rough fracture models were created using Matlab to assess roughness through the fractal dimension method.This research innovatively combined Boolean operations to establish three-dimensional rough fracture models,incorporating(Computational Fluid Dynamics)CFD-DEM(Discrete Element Method)with a bidirectional method for cosimulation.The proppant transport in fractures was categorized into three zones based on the difference in the turbulent kinetic energy.Artificially induced fracture roughness increases fluid retention and turbulence,causing plugging effects and limiting proppant flow into branch fractures.Additionally,compared with the superior deposition and significant support effects of the spherical proppant,the low-sphericity proppant traveled farther under fracturing fluid,inducing more pronounced plugging near curved fracture intersections;the variation in fracture intersection angles primarily impacted the wall shear stress within the flow field,indicating smaller angles led to higher shear energy at the intersection.Compared with the intersection angle of 30°,the height and area deposited in the 90 branch fracture increased by 52.25%and 65.33%,respectively:notably,injecting proppant from smaller to larger particles(S:M:L)and a low velocity effectively ensured fracture conductivity near the wellbore at joint roughness coefficient(JRC)≥46 while achieving satis-factory placement in the branch fracture,making it a recommended approach.
基金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.
基金Supported by the National High-Tech Research Project(GJSCB-HFGDY-2024-004)National Natural Science Foundation of China(12402305)+2 种基金Postdoctoral Fellowship Program of CPSF(GZC20232200)China Postdoctoral Science Foundation(2024M762703)Sichuan Science and Technology Program(2025ZNSFSC1352)。
文摘The Carter model is used to characterize the dynamic behaviors of fracture growth and fracturing fluid leakoff.A thermo-fluid coupling temperature response forward model is built considering the fluid flow and heat transfer in wellbore,fracture and reservoir.The influences of fracturing parameters and fracture parameters on the responses of distributed temperature sensing(DTS)are analyzed,and a diagnosis method of fracture parameters is presented based on the simulated annealing algorithm.A field case study is introduced to verify the model’s reliability.Typical V-shaped characteristics can be observed from the DTS responses in the multi-cluster fracturing process,with locations corresponding to the hydraulic fractures.The V-shape depth is shallower for a higher injection rate and longer fracturing and shut-in time.Also,the V-shape is wider for a higher fracture-surface leakoff coefficient,longer fracturing time and smaller fracture width.Additionally,the cooling effect near the wellbore continues to spread into the reservoir during the shut-in period,causing the DTS temperature to decrease instead of rise.Real-time monitoring and interpretation of DTS temperature data can help understand the fracture propagation during fracturing operation,so that immediate measures can be taken to improve the fracturing performance.
