Hydraulic fracturing technology has achieved remarkable results in improving the production of tight gas reservoirs,but its effectiveness is under the joint action of multiple factors of complexity.Traditional analysi...Hydraulic fracturing technology has achieved remarkable results in improving the production of tight gas reservoirs,but its effectiveness is under the joint action of multiple factors of complexity.Traditional analysis methods have limitations in dealing with these complex and interrelated factors,and it is difficult to fully reveal the actual contribution of each factor to the production.Machine learning-based methods explore the complex mapping relationships between large amounts of data to provide datadriven insights into the key factors driving production.In this study,a data-driven PCA-RF-VIM(Principal Component Analysis-Random Forest-Variable Importance Measures)approach of analyzing the importance of features is proposed to identify the key factors driving post-fracturing production.Four types of parameters,including log parameters,geological and reservoir physical parameters,hydraulic fracturing design parameters,and reservoir stimulation parameters,were inputted into the PCA-RF-VIM model.The model was trained using 6-fold cross-validation and grid search,and the relative importance ranking of each factor was finally obtained.In order to verify the validity of the PCA-RF-VIM model,a consolidation model that uses three other independent data-driven methods(Pearson correlation coefficient,RF feature significance analysis method,and XGboost feature significance analysis method)are applied to compare with the PCA-RF-VIM model.A comparison the two models shows that they contain almost the same parameters in the top ten,with only minor differences in one parameter.In combination with the reservoir characteristics,the reasonableness of the PCA-RF-VIM model is verified,and the importance ranking of the parameters by this method is more consistent with the reservoir characteristics of the study area.Ultimately,the ten parameters are selected as the controlling factors that have the potential to influence post-fracturing gas production,as the combined importance of these top ten parameters is 91.95%on driving natural gas production.Analyzing and obtaining these ten controlling factors provides engineers with a new insight into the reservoir selection for fracturing stimulation and fracturing parameter optimization to improve fracturing efficiency and productivity.展开更多
To analyze the differences in the transport and distribution of different types of proppants and to address issues such as the short effective support of proppant and poor placement in hydraulically intersecting fract...To analyze the differences in the transport and distribution of different types of proppants and to address issues such as the short effective support of proppant and poor placement in hydraulically intersecting fractures,this study considered the combined impact of geological-engineering factors on conductivity.Using reservoir production parameters and the discrete elementmethod,multispherical proppants were constructed.Additionally,a 3D fracture model,based on the specified conditions of the L block,employed coupled(Computational Fluid Dynamics)CFD-DEM(Discrete ElementMethod)for joint simulations to quantitatively analyze the transport and placement patterns of multispherical proppants in intersecting fractures.Results indicate that turbulent kinetic energy is an intrinsic factor affecting proppant transport.Moreover,the efficiency of placement and migration distance of low-sphericity quartz sand constructed by the DEM in the main fracture are significantly reduced compared to spherical ceramic proppants,with a 27.7%decrease in the volume fraction of the fracture surface,subsequently affecting the placement concentration and damaging fracture conductivity.Compared to small-angle fractures,controlling artificial and natural fractures to expand at angles of 45°to 60°increases the effective support length by approximately 20.6%.During hydraulic fracturing of gas wells,ensuring the fracture support area and post-closure conductivity can be achieved by controlling the sphericity of proppants and adjusting the perforation direction to control the direction of artificial fractures.展开更多
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.展开更多
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.展开更多
The prediction of the fracture plane orientation in fatigue is a scientific topic and remains relevant for every type of material. However, in this work, we compared the orientation of the fracture plane obtained expe...The prediction of the fracture plane orientation in fatigue is a scientific topic and remains relevant for every type of material. However, in this work, we compared the orientation of the fracture plane obtained experimentally through tests on specimens under multiaxial loading with that calculated by the variance method. In the statistical approach criteria, several methods have been developed but we have presented only one method, namely the variance method using the equivalent stress. She assumes that the fracture plane orientation is the one on which the variance of the equivalent stress is maximum. Three types of equivalent stress are defined for this method [1]: normal stress, shear stress and combined normal and shear stress. The results obtained were compared with experimental results for multiaxial cyclic stress states, and it emerges that the variance method for the case of combined loading is conservative as it gives a better prediction of the fracture plane.展开更多
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.展开更多
Hydraulic fracturing,an effective method for enhancing coal seam productivity,largely determines coalbed methane(CBM)production,which is significantly influenced by geological and engineering factors.This study focuse...Hydraulic fracturing,an effective method for enhancing coal seam productivity,largely determines coalbed methane(CBM)production,which is significantly influenced by geological and engineering factors.This study focuses on the L block to investigate the mechanisms influencing efficient fracture propagation and enhanced stimulated reservoir volume(SRV)in fracturing.To explore the mechanisms influencing effective fracture propagation and enhanced SRV,the L block was selected as the research object,with a comprehensive consideration of geological background,reservoir properties,and dynamic production data.By combining the discrete lattice method with numer-ical analysis and true triaxial experimental simulation,the fracture morphology of a single cluster and the propagation patterns of multiple clusters of complex fractures were obtained.Additionally,the optimization of temporary plugging timing and the fracture map under multiple factors were innovatively proposed.Results indicate that greater flow rate and viscosity can effectively overcome the stress shadow effect of the outermost fractures(1st and 6th clusters),increasing the fracture pressure of the single cluster and the equilibrium degree of multiple fracture propagation,thus forming a more complex fracture network.Moreover,when viscosity exceeds 45 pressure concentrates at fracture mPa⋅s,tips,promoting discontinuous propagation and reducing flow resistance.Conversely,increased gangue thickness and spacing between horizontal wells increase the vertical propagation pressure,suppressing fracture growth and reducing central flow velocity.This study provides a multi-cluster fracture propagation map for optimizing volumetric fracturing in coal seams and suggests that the optimal temporary plugging time significantly enhances the SRV.展开更多
Hydraulic fracturing plays a critical role in enhancing shale gas production in deep shale reservoirs.Conventional hydraulic fracturing simulation methods rely on prefabricated grids,which can be hindered by the chall...Hydraulic fracturing plays a critical role in enhancing shale gas production in deep shale reservoirs.Conventional hydraulic fracturing simulation methods rely on prefabricated grids,which can be hindered by the challenge of being computationally overpowered.This study proposes an efficient fracturing simulator to analyze fracture morphology during hydraulic fracturing processes in deep shale gas reservoirs.The simulator integrates the boundary element displacement discontinuity method and the finite volume method to model the fluid-solid coupling process by employing a pseudo-3D fracture model to calculate the fracture height.In particular,the Broyden iteration method was introduced to improve the computational efficiency and model robustness;it achieved a 46.6%reduction in computation time compared to the Newton-Raphson method.The influences of horizontal stress differences,natural fracture density,and natural fracture angle on the modified zone of the reservoir were simulated,and the following results were observed.(1)High stress difference reservoirs have smaller stimulated reservoir area than low stress difference reservoirs.(2)A higher natural fracture angle resulted in larger modification zones at low stress differences,while the effect of a natural fracture angle at high stress differences was not significant.(3)High-density and long natural fracture zones played a significant role in enhancing the stimulated reservoir area.These findings are critical for comprehending the impact of geological parameters on deep shale 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.展开更多
This paper presents a novel element differential method for modeling cracks in piezoelectric materials,aiming to simulate fracture behaviors and predict the fracture parameter known as the J-integral accurately.The me...This paper presents a novel element differential method for modeling cracks in piezoelectric materials,aiming to simulate fracture behaviors and predict the fracture parameter known as the J-integral accurately.The method leverages an efficient collocation technique to satisfy traction and electric charge equilibrium on the crack surface,aligning internal nodes with piezoelectric governing equations without needing integration or variational principles.It combines the strengths of the strong form collocation and finite element methods.The J-integral is derived analytically using the equivalent domain integral method,employing Green's formula and Gauss's divergence theorem to transform line integrals into area integrals for solving two-dimensional piezoelectric material problems.The accuracy of the method is validated through comparison with three typical examples,and it offers fracture prevention strategies for engineering piezoelectric structures under different electrical loading patterns.展开更多
Wellbore breakout is one of the critical issues in drilling due to the fact that the related problems result in additional costs and impact the drilling scheme severely.However,the majority of such wellbore breakout a...Wellbore breakout is one of the critical issues in drilling due to the fact that the related problems result in additional costs and impact the drilling scheme severely.However,the majority of such wellbore breakout analyses were based on continuum mechanics.In addition to failure in intact rocks,wellbore breakouts can also be initiated along natural discontinuities,e.g.weak planes and fractures.Furthermore,the conventional models in wellbore breakouts with uniform distribution fractures could not reflect the real drilling situation.This paper presents a fully coupled hydro-mechanical model of the SB-X well in the Tarim Basin,China for evaluating wellbore breakouts in heavily fractured rocks under anisotropic stress states using the distinct element method(DEM)and the discrete fracture network(DFN).The developed model was validated against caliper log measurement,and its stability study was carried out by stress and displacement analyses.A parametric study was performed to investigate the effects of the characteristics of fracture distribution(orientation and length)on borehole stability by sensitivity studies.Simulation results demonstrate that the increase of the standard deviation of orientation when the fracture direction aligns parallel or perpendicular to the principal stress direction aggravates borehole instability.Moreover,an elevation in the average fracture length causes the borehole failure to change from the direction of the minimum in-situ horizontal principal stress(i.e.the direction of wellbore breakouts)towards alternative directions,ultimately leading to the whole wellbore failure.These findings provide theoretical insights for predicting wellbore breakouts in heavily fractured rocks.展开更多
The glutenite reservoir is strongly heterogeneous due to the random distribution of gravels, making it challenging to perform hydraulic fracturing effectively. To solve this issue, it is essential to study interaction...The glutenite reservoir is strongly heterogeneous due to the random distribution of gravels, making it challenging to perform hydraulic fracturing effectively. To solve this issue, it is essential to study interaction behavior between hydraulic fractures(HFs) and gravels. A coupled hydro-mechanical model is proposed for HF propagation in glutenite using a grain-based discrete element method. This paper first investigates the dynamic evolution of HFs in glutenite, then analyzes the influences of various factors such as horizontal stress difference(Δσ), minimum horizontal stress(σh), gravel content(Vg), gravel size(dg), and stiffness ratio of gravel to matrix(Rs) on HF propagation geometries. Results show that penetrating the gravel is the primary HF-gravel interaction behavior, which follows sequential and staggered initiation modes. Bypassing the gravel is the secondary behavior, which obeys the sequential initiation mode and occurs when the orientation of the gravel boundary is inclined to the maximum horizontal stress(σH). An offset along the gravel boundary is usually formed while penetrating gravels, and the offsets may cause fracture widths to decrease by 37.8%-84.4%. Even if stress dominates the direction of HF propagation, HFs still tend to deflect within gravels. The deviation angle from σH decreases with rising Δσand increases with the increase of dgand Rs. Additionally, intra-gravel shear HFs(IGS-HFs) are prone to be generated in coarse-grained glutenite under high Δσ, while more gravel-bypassing shear HFs(GBSHFs) tend to be created in argillaceous glutenite with high Rsthan in sandy glutenite with low Rs. The findings above prompt the emergence of a novel HF propagation pattern in glutenite, which helps to understand the real HF geometries and to provide theoretical guidance for treatments in the field.展开更多
Reservoirs with a group of vertical fractures in a vertical transversely isotropic(VTI)background are considered as orthorhombic(ORT)medium.However,fracture detection in ORT medium using seismic inversion methods rema...Reservoirs with a group of vertical fractures in a vertical transversely isotropic(VTI)background are considered as orthorhombic(ORT)medium.However,fracture detection in ORT medium using seismic inversion methods remains challenging,as it requires the estimation of more than eight parameters.Assuming the reservoir to be a weakly anisotropic ORT medium with small contrasts in the background elastic parameters,a new azimuthal elastic impedance equation was first derived using parameter combinations and mathematical approximations.This equation exhibited almost the same accuracy as the original equation and contained only six model parameters:the compression modulus,anisotropic shear modulus,anisotropic compression modulus,density,normal fracture weakness,and tangential fracture weakness.Subsequently,a stepwise inversion method using second-order derivatives of the elastic impedance was developed to estimate these parameters.Moreover,the Thomsen anisotropy parameter,epsilon,was estimated from the inversion results using the ratio of the anisotropic compression modulus to the compression modulus.Synthetic examples with moderate noise and field data examples confirm the feasibility and effectiveness of the inversion method.The proposed method exhibited accuracy similar to that of previous inversion strategies and could predict richer vertical fracture information.Ultimately,the method was applied to a three-dimensional work area,and the predictions were consistent with logging and geological a priori information,confirming the effectiveness of this method.Summarily,the proposed stepwise inversion method can alleviate the uncertainty of multi-parameter inversion in ORT medium,thereby improving the reliability of fracture detection.展开更多
Horizontal well drilling and multi-stage hydraulic fracturing are key technologies for the development of shale gas reservoirs.Instantaneous acquisition of hydraulic fracture parameters is crucial for evaluating fract...Horizontal well drilling and multi-stage hydraulic fracturing are key technologies for the development of shale gas reservoirs.Instantaneous acquisition of hydraulic fracture parameters is crucial for evaluating fracturing effectiveness,optimizing processes,and predicting gas productivity.This paper establishes a transient flow model for shale gas wells based on the boundary element method,achieving the characterization of stimulated reservoir volume for a single stage.By integrating pressure monitoring data following the pumping shut-in period of hydraulic fracturing for well testing interpretation,a workflow for inverting fracture parameters of shale gas wells is established.This new method eliminates the need for prolonged production testing and can interpret parameters of individual hydraulic fracture segments,offering significant advantages over the conventional pressure transient analysismethod.The practical application of thismethodology was conducted on 10 shale gaswellswithin the Changning shale gas block of Sichuan,China.The results show a high correlation between the interpreted single-stage total length and surface area of hydraulic fractures and the outcomes of gas production profile tests.Additionally,significant correlations are observed between these parameters and cluster number,horizontal stress difference,and natural fracture density.This demonstrates the effectiveness of the proposed fracture parameter inversion method and the feasibility of field application.The findings of this study aim to provide solutions and references for the inversion of fracture parameters in shale gas wells.展开更多
Ground hydraulic fracturing plays a crucial role in controlling the far-field hard roof,making it imperative to identify the most suitable target stratum for effective control.Physical experiments are conducted based ...Ground hydraulic fracturing plays a crucial role in controlling the far-field hard roof,making it imperative to identify the most suitable target stratum for effective control.Physical experiments are conducted based on engineering properties to simulate the gradual collapse of the roof during longwall top coal caving(LTCC).A numerical model is established using the material point method(MPM)and the strain-softening damage constitutive model according to the structure of the physical model.Numerical simulations are conducted to analyze the LTCC process under different hard roofs for ground hydraulic fracturing.The results show that ground hydraulic fracturing releases the energy and stress of the target stratum,resulting in a substantial lag in the fracturing of the overburden before collapse occurs in the hydraulic fracturing stratum.Ground hydraulic fracturing of a low hard roof reduces the lag effect of hydraulic fractures,dissipates the energy consumed by the fracture of the hard roof,and reduces the abutment stress.Therefore,it is advisable to prioritize the selection of the lower hard roof as the target stratum.展开更多
The exploitation of shale gas is promising due to depletion of the conventional energy and intensification of the greenhouse effect.In this paper,we proposed a heat-fluid-solid coupling damage model of supercritical C...The exploitation of shale gas is promising due to depletion of the conventional energy and intensification of the greenhouse effect.In this paper,we proposed a heat-fluid-solid coupling damage model of supercritical CO_(2)(SC-CO_(2))compound fracturing which is expected to be an efficient and environmentally friendly way to develop shale gas.The coupling model is solved by the finite element method,and the results are in good agreement with the analytical solutions and fracturing experiments.Based on this model,the fracture propagation characteristics at the two stages of compound fracturing are studied and the influence of pressurization rate,in situ stress,bedding angle,and other factors are considered.The results show that at the SC-CO_(2)fracturing stage,a lower pressurization rate is conducive to formation of the branches around main fractures,while a higher pressurization rate inhibits formation of the branches around main fractures and promotes formation of the main fractures.Both bedding and in situ stress play a dominant role in the fracture propagation.When the in situ stress ratio(δ_(x)/δ_(y))is 1,the presence of bedding can reduce the initiation pressure and failure pressure.Nevertheless,it will cause the fracture to propagate along the bedding direction,reducing the fracture complexity.In rocks without bedding,hydraulic fracturing has the lengthening and widening effects for SC-CO_(2)induced fracture.In shale,fractures induced at the hydraulic fracturing stage are more likely to be dominated by in situ stresses and have a shorter reorientation radius.Therefore,fracture branches propagating along the maximum principal stress direction may be generated around the main fractures induced by SC-CO_(2)at the hydraulic fracturing stage.When the branches converge with the main fractures,fracture zones are easily formed,and thus the fracture complexity and damage area can be significantly increased.The results are instructive for the design and application of SC-CO_(2)compound fracturing.展开更多
As a calculation method based on the Galerkin variation,the numerical manifold method(NMM)adopts a double covering system,which can easily deal with discontinuous deformation problems and has a high calculation accura...As a calculation method based on the Galerkin variation,the numerical manifold method(NMM)adopts a double covering system,which can easily deal with discontinuous deformation problems and has a high calculation accuracy.Aiming at the thermo-mechanical(TM)coupling problem of fractured rock masses,this study uses the NMM to simulate the processes of crack initiation and propagation in a rock mass under the influence of temperature field,deduces related system equations,and proposes a penalty function method to deal with boundary conditions.Numerical examples are employed to confirm the effectiveness and high accuracy of this method.By the thermal stress analysis of a thick-walled cylinder(TWC),the simulation of cracking in the TWC under heating and cooling conditions,and the simulation of thermal cracking of the SwedishÄspöPillar Stability Experiment(APSE)rock column,the thermal stress,and TM coupling are obtained.The numerical simulation results are in good agreement with the test data and other numerical results,thus verifying the effectiveness of the NMM in dealing with thermal stress and crack propagation problems of fractured rock masses.展开更多
In this study,a fully coupled hydromechanical model within the extended finite element method(XFEM)-based cohesive zone method(CZM)is employed to investigate the simultaneous height growth behavior of multi-cluster hy...In this study,a fully coupled hydromechanical model within the extended finite element method(XFEM)-based cohesive zone method(CZM)is employed to investigate the simultaneous height growth behavior of multi-cluster hydraulic fractures in layered porous reservoirs with modulus contrast.The coupled hydromechanical model is first verified against an analytical solution and a laboratory experiment.Then,the fracture geometry(e.g.height,aperture,and area)and fluid pressure evolutions of multiple hydraulic fractures placed in a porous reservoir interbedded with alternating stiff and soft layers are investigated using the model.The stress and pore pressure distributions within the layered reservoir during fluid injection are also presented.The simulation results reveal that stress umbrellas are easily to form among multiple hydraulic fractures’tips when propagating in soft layers,which impedes the simultaneous height growth.It is also observed that the impediment effect of soft layer is much more significant in the fractures suppressed by the preferential growth of adjoining fractures.After that,the combined effect of in situ stress ratio and fracturing spacing on the multi-fracture height growth is presented,and the results elucidate the influence of in situ stress ratio on the height growth behavior depending on the fracture spacing.Finally,it is found that the inclusion of soft layers changes the aperture distribution of outmost and interior hydraulic fractures.The results obtained from this study may provide some insights on the understanding of hydraulic fracture height containment observed in filed.展开更多
During the implementation of CO_(2) fracturing for oil and gas development,the force transfer effect caused by the unsteady flow of high-pressure CO_(2) fluid can lead to forced vibration of the tubing and ensuing str...During the implementation of CO_(2) fracturing for oil and gas development,the force transfer effect caused by the unsteady flow of high-pressure CO_(2) fluid can lead to forced vibration of the tubing and ensuing structural fatigue.In this study,a forced vibration analysis of tubing under CO_(2) fracturing conditions is carried out by taking into account the fluid-structure coupling and related interaction forces by means of the method of characteristics(MOC).The results show that for every 1 m^(3)/min increase in pumping displacement,the fluid flow rate increases up to 3.67 m/s.The flow pressure in the pipe tends to be consistent with the pumping pressure at the initial stage and then decreases with an increase in the pump starting time.When the pumping pressure increases by 10 MPa,the additional stress in the tubing increases by 11.8%,and the peak value of the additional stress at the bottom of the well is the largest.The temperature in the tubing grows with well depth,which causes a phase change in CO_(2) due to heat absorption.At this time the pressure in the tubing decreases,the fluid flow rate increases by about 1.12 m/s,and the additional stress grows by about 1.5 MPa.展开更多
To analyze and predict the mechanical behaviors of deep hard rocks,some key issues concerning rock fracturing mechanics for deep hard rock excavations are discussed.First,a series of apparatuses and methods have been ...To analyze and predict the mechanical behaviors of deep hard rocks,some key issues concerning rock fracturing mechanics for deep hard rock excavations are discussed.First,a series of apparatuses and methods have been developed to test the mechanical properties and fracturing behaviors of hard rocks under high true triaxial stress paths.Evolution mechanisms of stress-induced disasters in deep hard rock excavations,such as spalling,deep cracking,massive roof collapse,large deformation and rockbursts,have been recognized.The analytical theory for the fracturing process of hard rock masses,including the three-dimensional failure criterion,stress-induced mechanical model,fracturing degree index,energy release index and numerical method,has been established.The cracking-restraint method is developed for mitigating or controlling rock spalling,deep cracking and massive collapse of deep hard rocks.An energy-controlled method is also proposed for the prevention of rockbursts.Finally,two typical cases are used to illustrate the application of the proposed methodology in the Baihetan caverns and Bayu tunnels of China.展开更多
基金funded by the Key Research and Development Program of Shaanxi,China(No.2024GX-YBXM-503)the National Natural Science Foundation of China(No.51974254)。
文摘Hydraulic fracturing technology has achieved remarkable results in improving the production of tight gas reservoirs,but its effectiveness is under the joint action of multiple factors of complexity.Traditional analysis methods have limitations in dealing with these complex and interrelated factors,and it is difficult to fully reveal the actual contribution of each factor to the production.Machine learning-based methods explore the complex mapping relationships between large amounts of data to provide datadriven insights into the key factors driving production.In this study,a data-driven PCA-RF-VIM(Principal Component Analysis-Random Forest-Variable Importance Measures)approach of analyzing the importance of features is proposed to identify the key factors driving post-fracturing production.Four types of parameters,including log parameters,geological and reservoir physical parameters,hydraulic fracturing design parameters,and reservoir stimulation parameters,were inputted into the PCA-RF-VIM model.The model was trained using 6-fold cross-validation and grid search,and the relative importance ranking of each factor was finally obtained.In order to verify the validity of the PCA-RF-VIM model,a consolidation model that uses three other independent data-driven methods(Pearson correlation coefficient,RF feature significance analysis method,and XGboost feature significance analysis method)are applied to compare with the PCA-RF-VIM model.A comparison the two models shows that they contain almost the same parameters in the top ten,with only minor differences in one parameter.In combination with the reservoir characteristics,the reasonableness of the PCA-RF-VIM model is verified,and the importance ranking of the parameters by this method is more consistent with the reservoir characteristics of the study area.Ultimately,the ten parameters are selected as the controlling factors that have the potential to influence post-fracturing gas production,as the combined importance of these top ten parameters is 91.95%on driving natural gas production.Analyzing and obtaining these ten controlling factors provides engineers with a new insight into the reservoir selection for fracturing stimulation and fracturing parameter optimization to improve fracturing efficiency and productivity.
基金funded by the project of the Major Scientific and Technological Projects of CNOOC in the 14th Five-Year Plan(No.KJGG2022-0701)the CNOOC Research Institute(No.2020PFS-03).
文摘To analyze the differences in the transport and distribution of different types of proppants and to address issues such as the short effective support of proppant and poor placement in hydraulically intersecting fractures,this study considered the combined impact of geological-engineering factors on conductivity.Using reservoir production parameters and the discrete elementmethod,multispherical proppants were constructed.Additionally,a 3D fracture model,based on the specified conditions of the L block,employed coupled(Computational Fluid Dynamics)CFD-DEM(Discrete ElementMethod)for joint simulations to quantitatively analyze the transport and placement patterns of multispherical proppants in intersecting fractures.Results indicate that turbulent kinetic energy is an intrinsic factor affecting proppant transport.Moreover,the efficiency of placement and migration distance of low-sphericity quartz sand constructed by the DEM in the main fracture are significantly reduced compared to spherical ceramic proppants,with a 27.7%decrease in the volume fraction of the fracture surface,subsequently affecting the placement concentration and damaging fracture conductivity.Compared to small-angle fractures,controlling artificial and natural fractures to expand at angles of 45°to 60°increases the effective support length by approximately 20.6%.During hydraulic fracturing of gas wells,ensuring the fracture support area and post-closure conductivity can be achieved by controlling the sphericity of proppants and adjusting the perforation direction to control the direction of artificial fractures.
基金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.
基金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.
文摘The prediction of the fracture plane orientation in fatigue is a scientific topic and remains relevant for every type of material. However, in this work, we compared the orientation of the fracture plane obtained experimentally through tests on specimens under multiaxial loading with that calculated by the variance method. In the statistical approach criteria, several methods have been developed but we have presented only one method, namely the variance method using the equivalent stress. She assumes that the fracture plane orientation is the one on which the variance of the equivalent stress is maximum. Three types of equivalent stress are defined for this method [1]: normal stress, shear stress and combined normal and shear stress. The results obtained were compared with experimental results for multiaxial cyclic stress states, and it emerges that the variance method for the case of combined loading is conservative as it gives a better prediction of the fracture plane.
基金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.
基金the project of the State Key Laboratory of Petroleum Resources and Engineering(No.PRE/open-2307)the CNOOC Research Institute(No.2020PFS-03).
文摘Hydraulic fracturing,an effective method for enhancing coal seam productivity,largely determines coalbed methane(CBM)production,which is significantly influenced by geological and engineering factors.This study focuses on the L block to investigate the mechanisms influencing efficient fracture propagation and enhanced stimulated reservoir volume(SRV)in fracturing.To explore the mechanisms influencing effective fracture propagation and enhanced SRV,the L block was selected as the research object,with a comprehensive consideration of geological background,reservoir properties,and dynamic production data.By combining the discrete lattice method with numer-ical analysis and true triaxial experimental simulation,the fracture morphology of a single cluster and the propagation patterns of multiple clusters of complex fractures were obtained.Additionally,the optimization of temporary plugging timing and the fracture map under multiple factors were innovatively proposed.Results indicate that greater flow rate and viscosity can effectively overcome the stress shadow effect of the outermost fractures(1st and 6th clusters),increasing the fracture pressure of the single cluster and the equilibrium degree of multiple fracture propagation,thus forming a more complex fracture network.Moreover,when viscosity exceeds 45 pressure concentrates at fracture mPa⋅s,tips,promoting discontinuous propagation and reducing flow resistance.Conversely,increased gangue thickness and spacing between horizontal wells increase the vertical propagation pressure,suppressing fracture growth and reducing central flow velocity.This study provides a multi-cluster fracture propagation map for optimizing volumetric fracturing in coal seams and suggests that the optimal temporary plugging time significantly enhances the SRV.
文摘Hydraulic fracturing plays a critical role in enhancing shale gas production in deep shale reservoirs.Conventional hydraulic fracturing simulation methods rely on prefabricated grids,which can be hindered by the challenge of being computationally overpowered.This study proposes an efficient fracturing simulator to analyze fracture morphology during hydraulic fracturing processes in deep shale gas reservoirs.The simulator integrates the boundary element displacement discontinuity method and the finite volume method to model the fluid-solid coupling process by employing a pseudo-3D fracture model to calculate the fracture height.In particular,the Broyden iteration method was introduced to improve the computational efficiency and model robustness;it achieved a 46.6%reduction in computation time compared to the Newton-Raphson method.The influences of horizontal stress differences,natural fracture density,and natural fracture angle on the modified zone of the reservoir were simulated,and the following results were observed.(1)High stress difference reservoirs have smaller stimulated reservoir area than low stress difference reservoirs.(2)A higher natural fracture angle resulted in larger modification zones at low stress differences,while the effect of a natural fracture angle at high stress differences was not significant.(3)High-density and long natural fracture zones played a significant role in enhancing the stimulated reservoir area.These findings are critical for comprehending the impact of geological parameters on deep shale 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.
基金Financial support of this work by the Technology Development program of China(Grant No.2022204B003)National Natural Science Foundation of China(12272083 and 12172078)the Fundamental Research Funds for the Central Universities(DUT24YJ136)is gratefully acknowledged.
文摘This paper presents a novel element differential method for modeling cracks in piezoelectric materials,aiming to simulate fracture behaviors and predict the fracture parameter known as the J-integral accurately.The method leverages an efficient collocation technique to satisfy traction and electric charge equilibrium on the crack surface,aligning internal nodes with piezoelectric governing equations without needing integration or variational principles.It combines the strengths of the strong form collocation and finite element methods.The J-integral is derived analytically using the equivalent domain integral method,employing Green's formula and Gauss's divergence theorem to transform line integrals into area integrals for solving two-dimensional piezoelectric material problems.The accuracy of the method is validated through comparison with three typical examples,and it offers fracture prevention strategies for engineering piezoelectric structures under different electrical loading patterns.
基金supported by National Natural Science Foundation of China(Grant Nos.52074312 and 52211530097)CNPC Science and Technology Innovation Foundation(Grant No.2021DQ02-0505).
文摘Wellbore breakout is one of the critical issues in drilling due to the fact that the related problems result in additional costs and impact the drilling scheme severely.However,the majority of such wellbore breakout analyses were based on continuum mechanics.In addition to failure in intact rocks,wellbore breakouts can also be initiated along natural discontinuities,e.g.weak planes and fractures.Furthermore,the conventional models in wellbore breakouts with uniform distribution fractures could not reflect the real drilling situation.This paper presents a fully coupled hydro-mechanical model of the SB-X well in the Tarim Basin,China for evaluating wellbore breakouts in heavily fractured rocks under anisotropic stress states using the distinct element method(DEM)and the discrete fracture network(DFN).The developed model was validated against caliper log measurement,and its stability study was carried out by stress and displacement analyses.A parametric study was performed to investigate the effects of the characteristics of fracture distribution(orientation and length)on borehole stability by sensitivity studies.Simulation results demonstrate that the increase of the standard deviation of orientation when the fracture direction aligns parallel or perpendicular to the principal stress direction aggravates borehole instability.Moreover,an elevation in the average fracture length causes the borehole failure to change from the direction of the minimum in-situ horizontal principal stress(i.e.the direction of wellbore breakouts)towards alternative directions,ultimately leading to the whole wellbore failure.These findings provide theoretical insights for predicting wellbore breakouts in heavily fractured rocks.
基金supported by the National Natural Science Foundation of China(Grant No.52304003)the Natural Science Foundation of Sichuan Province(Grant No.2024NSFSC0961)the Postdoctoral Fellowship Program of CPSF(Grant No.GZB20230090).
文摘The glutenite reservoir is strongly heterogeneous due to the random distribution of gravels, making it challenging to perform hydraulic fracturing effectively. To solve this issue, it is essential to study interaction behavior between hydraulic fractures(HFs) and gravels. A coupled hydro-mechanical model is proposed for HF propagation in glutenite using a grain-based discrete element method. This paper first investigates the dynamic evolution of HFs in glutenite, then analyzes the influences of various factors such as horizontal stress difference(Δσ), minimum horizontal stress(σh), gravel content(Vg), gravel size(dg), and stiffness ratio of gravel to matrix(Rs) on HF propagation geometries. Results show that penetrating the gravel is the primary HF-gravel interaction behavior, which follows sequential and staggered initiation modes. Bypassing the gravel is the secondary behavior, which obeys the sequential initiation mode and occurs when the orientation of the gravel boundary is inclined to the maximum horizontal stress(σH). An offset along the gravel boundary is usually formed while penetrating gravels, and the offsets may cause fracture widths to decrease by 37.8%-84.4%. Even if stress dominates the direction of HF propagation, HFs still tend to deflect within gravels. The deviation angle from σH decreases with rising Δσand increases with the increase of dgand Rs. Additionally, intra-gravel shear HFs(IGS-HFs) are prone to be generated in coarse-grained glutenite under high Δσ, while more gravel-bypassing shear HFs(GBSHFs) tend to be created in argillaceous glutenite with high Rsthan in sandy glutenite with low Rs. The findings above prompt the emergence of a novel HF propagation pattern in glutenite, which helps to understand the real HF geometries and to provide theoretical guidance for treatments in the field.
基金sponsorship of the National Natural Science Foundation of China(42430809,42274157,42030103,42404132)the Fund of State Key Laboratory of Deep Oil and Gas,China University of Petroleum(East China)(SKLDOG2024-ZYTS-02)+5 种基金the Postdoctoral Fellowship Program of CPSF(GZB20240850)the Postdoctoral Project of Qingdao(QDBSH20240102082)the Fundamental Research Funds for the Central Universities(24CX07004A,24CX06036A)the CNPC Innovation Fund(2024DQ02-0505,2024DQ02-0136)the Innovation fund project for graduate student of China University of Petroleum(East China)the Fundamental Research Funds for the Central Universities(24CX04002A).
文摘Reservoirs with a group of vertical fractures in a vertical transversely isotropic(VTI)background are considered as orthorhombic(ORT)medium.However,fracture detection in ORT medium using seismic inversion methods remains challenging,as it requires the estimation of more than eight parameters.Assuming the reservoir to be a weakly anisotropic ORT medium with small contrasts in the background elastic parameters,a new azimuthal elastic impedance equation was first derived using parameter combinations and mathematical approximations.This equation exhibited almost the same accuracy as the original equation and contained only six model parameters:the compression modulus,anisotropic shear modulus,anisotropic compression modulus,density,normal fracture weakness,and tangential fracture weakness.Subsequently,a stepwise inversion method using second-order derivatives of the elastic impedance was developed to estimate these parameters.Moreover,the Thomsen anisotropy parameter,epsilon,was estimated from the inversion results using the ratio of the anisotropic compression modulus to the compression modulus.Synthetic examples with moderate noise and field data examples confirm the feasibility and effectiveness of the inversion method.The proposed method exhibited accuracy similar to that of previous inversion strategies and could predict richer vertical fracture information.Ultimately,the method was applied to a three-dimensional work area,and the predictions were consistent with logging and geological a priori information,confirming the effectiveness of this method.Summarily,the proposed stepwise inversion method can alleviate the uncertainty of multi-parameter inversion in ORT medium,thereby improving the reliability of fracture detection.
基金funded by the Science and Technology Cooperation Project of the CNPC-SWPU Innovation Alliance,grant numbers“2020CX020202,2020CX030202 and 2020CX010403”.
文摘Horizontal well drilling and multi-stage hydraulic fracturing are key technologies for the development of shale gas reservoirs.Instantaneous acquisition of hydraulic fracture parameters is crucial for evaluating fracturing effectiveness,optimizing processes,and predicting gas productivity.This paper establishes a transient flow model for shale gas wells based on the boundary element method,achieving the characterization of stimulated reservoir volume for a single stage.By integrating pressure monitoring data following the pumping shut-in period of hydraulic fracturing for well testing interpretation,a workflow for inverting fracture parameters of shale gas wells is established.This new method eliminates the need for prolonged production testing and can interpret parameters of individual hydraulic fracture segments,offering significant advantages over the conventional pressure transient analysismethod.The practical application of thismethodology was conducted on 10 shale gaswellswithin the Changning shale gas block of Sichuan,China.The results show a high correlation between the interpreted single-stage total length and surface area of hydraulic fractures and the outcomes of gas production profile tests.Additionally,significant correlations are observed between these parameters and cluster number,horizontal stress difference,and natural fracture density.This demonstrates the effectiveness of the proposed fracture parameter inversion method and the feasibility of field application.The findings of this study aim to provide solutions and references for the inversion of fracture parameters in shale gas wells.
基金the National Natural Science Foundation of China(No.51974042)National Key Research and Development Program of China(No.2023YFC3009005).
文摘Ground hydraulic fracturing plays a crucial role in controlling the far-field hard roof,making it imperative to identify the most suitable target stratum for effective control.Physical experiments are conducted based on engineering properties to simulate the gradual collapse of the roof during longwall top coal caving(LTCC).A numerical model is established using the material point method(MPM)and the strain-softening damage constitutive model according to the structure of the physical model.Numerical simulations are conducted to analyze the LTCC process under different hard roofs for ground hydraulic fracturing.The results show that ground hydraulic fracturing releases the energy and stress of the target stratum,resulting in a substantial lag in the fracturing of the overburden before collapse occurs in the hydraulic fracturing stratum.Ground hydraulic fracturing of a low hard roof reduces the lag effect of hydraulic fractures,dissipates the energy consumed by the fracture of the hard roof,and reduces the abutment stress.Therefore,it is advisable to prioritize the selection of the lower hard roof as the target stratum.
基金the funding support from the National Natural Science Foundation of China(Nos.52274093 and 52004236).
文摘The exploitation of shale gas is promising due to depletion of the conventional energy and intensification of the greenhouse effect.In this paper,we proposed a heat-fluid-solid coupling damage model of supercritical CO_(2)(SC-CO_(2))compound fracturing which is expected to be an efficient and environmentally friendly way to develop shale gas.The coupling model is solved by the finite element method,and the results are in good agreement with the analytical solutions and fracturing experiments.Based on this model,the fracture propagation characteristics at the two stages of compound fracturing are studied and the influence of pressurization rate,in situ stress,bedding angle,and other factors are considered.The results show that at the SC-CO_(2)fracturing stage,a lower pressurization rate is conducive to formation of the branches around main fractures,while a higher pressurization rate inhibits formation of the branches around main fractures and promotes formation of the main fractures.Both bedding and in situ stress play a dominant role in the fracture propagation.When the in situ stress ratio(δ_(x)/δ_(y))is 1,the presence of bedding can reduce the initiation pressure and failure pressure.Nevertheless,it will cause the fracture to propagate along the bedding direction,reducing the fracture complexity.In rocks without bedding,hydraulic fracturing has the lengthening and widening effects for SC-CO_(2)induced fracture.In shale,fractures induced at the hydraulic fracturing stage are more likely to be dominated by in situ stresses and have a shorter reorientation radius.Therefore,fracture branches propagating along the maximum principal stress direction may be generated around the main fractures induced by SC-CO_(2)at the hydraulic fracturing stage.When the branches converge with the main fractures,fracture zones are easily formed,and thus the fracture complexity and damage area can be significantly increased.The results are instructive for the design and application of SC-CO_(2)compound fracturing.
基金supported by the National Natural Science Foundation of China(Grant No.42277165)the Fundamental Research Funds for the Central Universities,China University of Geosciences(Wuhan)(Grant No.CUGCJ1821)the National Overseas Study Fund(Grant No.202106410040).
文摘As a calculation method based on the Galerkin variation,the numerical manifold method(NMM)adopts a double covering system,which can easily deal with discontinuous deformation problems and has a high calculation accuracy.Aiming at the thermo-mechanical(TM)coupling problem of fractured rock masses,this study uses the NMM to simulate the processes of crack initiation and propagation in a rock mass under the influence of temperature field,deduces related system equations,and proposes a penalty function method to deal with boundary conditions.Numerical examples are employed to confirm the effectiveness and high accuracy of this method.By the thermal stress analysis of a thick-walled cylinder(TWC),the simulation of cracking in the TWC under heating and cooling conditions,and the simulation of thermal cracking of the SwedishÄspöPillar Stability Experiment(APSE)rock column,the thermal stress,and TM coupling are obtained.The numerical simulation results are in good agreement with the test data and other numerical results,thus verifying the effectiveness of the NMM in dealing with thermal stress and crack propagation problems of fractured rock masses.
文摘In this study,a fully coupled hydromechanical model within the extended finite element method(XFEM)-based cohesive zone method(CZM)is employed to investigate the simultaneous height growth behavior of multi-cluster hydraulic fractures in layered porous reservoirs with modulus contrast.The coupled hydromechanical model is first verified against an analytical solution and a laboratory experiment.Then,the fracture geometry(e.g.height,aperture,and area)and fluid pressure evolutions of multiple hydraulic fractures placed in a porous reservoir interbedded with alternating stiff and soft layers are investigated using the model.The stress and pore pressure distributions within the layered reservoir during fluid injection are also presented.The simulation results reveal that stress umbrellas are easily to form among multiple hydraulic fractures’tips when propagating in soft layers,which impedes the simultaneous height growth.It is also observed that the impediment effect of soft layer is much more significant in the fractures suppressed by the preferential growth of adjoining fractures.After that,the combined effect of in situ stress ratio and fracturing spacing on the multi-fracture height growth is presented,and the results elucidate the influence of in situ stress ratio on the height growth behavior depending on the fracture spacing.Finally,it is found that the inclusion of soft layers changes the aperture distribution of outmost and interior hydraulic fractures.The results obtained from this study may provide some insights on the understanding of hydraulic fracture height containment observed in filed.
基金funded by National Natural Science Foundation of China(Grant No.52105209)funded by the National Engineering Research Centre for Oil and Gas Drilling and Completion Technology(Grant No.F2023199).
文摘During the implementation of CO_(2) fracturing for oil and gas development,the force transfer effect caused by the unsteady flow of high-pressure CO_(2) fluid can lead to forced vibration of the tubing and ensuing structural fatigue.In this study,a forced vibration analysis of tubing under CO_(2) fracturing conditions is carried out by taking into account the fluid-structure coupling and related interaction forces by means of the method of characteristics(MOC).The results show that for every 1 m^(3)/min increase in pumping displacement,the fluid flow rate increases up to 3.67 m/s.The flow pressure in the pipe tends to be consistent with the pumping pressure at the initial stage and then decreases with an increase in the pump starting time.When the pumping pressure increases by 10 MPa,the additional stress in the tubing increases by 11.8%,and the peak value of the additional stress at the bottom of the well is the largest.The temperature in the tubing grows with well depth,which causes a phase change in CO_(2) due to heat absorption.At this time the pressure in the tubing decreases,the fluid flow rate increases by about 1.12 m/s,and the additional stress grows by about 1.5 MPa.
基金financial support from the National Natural Science Foundation of China(Grant Nos.51839003 and 41827806)Liaoning Revitalization Talents Program of China(Grant No.XLYCYSZX1902)。
文摘To analyze and predict the mechanical behaviors of deep hard rocks,some key issues concerning rock fracturing mechanics for deep hard rock excavations are discussed.First,a series of apparatuses and methods have been developed to test the mechanical properties and fracturing behaviors of hard rocks under high true triaxial stress paths.Evolution mechanisms of stress-induced disasters in deep hard rock excavations,such as spalling,deep cracking,massive roof collapse,large deformation and rockbursts,have been recognized.The analytical theory for the fracturing process of hard rock masses,including the three-dimensional failure criterion,stress-induced mechanical model,fracturing degree index,energy release index and numerical method,has been established.The cracking-restraint method is developed for mitigating or controlling rock spalling,deep cracking and massive collapse of deep hard rocks.An energy-controlled method is also proposed for the prevention of rockbursts.Finally,two typical cases are used to illustrate the application of the proposed methodology in the Baihetan caverns and Bayu tunnels of China.
