In the development of coalbed methane(CBM)reservoirs using multistage fractured horizontal wells,there often exist areas that are either repeatedly stimulated or completely unstimulated between fracturing stages,leadi...In the development of coalbed methane(CBM)reservoirs using multistage fractured horizontal wells,there often exist areas that are either repeatedly stimulated or completely unstimulated between fracturing stages,leading to suboptimal reservoir performance.Currently,there is no well-established method for accurately evaluating the effectiveness of such stimulation.This study introduces,for the first time,the concept of the Fracture Network Bridging Coefficient(FNBC)as a novel metric to assess stimulation performance.By quantitatively coupling the proportions of unstimulated and overstimulated volumes,the FNBC effectively characterizes the connectivity and efficiency of the fracture network.A background grid calibration method is developed to quantify the stage-controlled volume,effectively stimulated volume,unstimulated volume,and repeatedly stimulated volume among different stages of horizontal wells.Furthermore,an optimization model is constructed by taking the FNBC as the objective function and the fracturing injection rate and fluid volume as optimization variables.The Simultaneous Perturbation Stochastic Approximation(SPSA)algorithm is employed to iteratively perturb and optimize these variables,progressively improving the FNBC until the optimal displacement rate and fluid volume corresponding to the maximum FNBC are obtained.Field application in a typical CBM multistage fractured horizontal well in China demonstrates that the FNBC increased from 0.358 to 0.539(a 50.6% improvement),with the injection rate rising from 16 m^(3)/min to 24 m^(3)/min and the average fluid volume per stage increasing from 2490 m^(3) to 3192 m^(3),significantly enhancing the stimulation effectiveness.This research provides theoretical support for designing high-efficiency stimulation strategies in unconventional reservoirs under dynamic limits.展开更多
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.展开更多
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.展开更多
This review summarizes the development of particle-based numerical manifold method(PNMM)and its applications to rock dynamics.The fundamental principle of numerical manifold method(NMM)is first briefly introduced.Then...This review summarizes the development of particle-based numerical manifold method(PNMM)and its applications to rock dynamics.The fundamental principle of numerical manifold method(NMM)is first briefly introduced.Then,the history of the newly developed PNMM is given.Basic idea of PNMM and its simulation procedure are presented.Considering that PNMM could be regarded as an NMM-based model,a comparison of PNMM and NMM is discussed from several points of view in this paper.Besides,accomplished applications of PNMM to the dynamic rock fracturing are also reviewed.Finally,some recommendations are provided for the future work of PNMM.展开更多
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.展开更多
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.展开更多
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.展开更多
Hydraulic fracturing （HF） technique has been extensively used for the exploitation of unconventional oiland gas reservoirs. HF enhances the connectivity of less permeable oil and gas-bearing rock formationsby fluid ...Hydraulic fracturing （HF） technique has been extensively used for the exploitation of unconventional oiland gas reservoirs. HF enhances the connectivity of less permeable oil and gas-bearing rock formationsby fluid injection, which creates an interconnected fracture network and increases the hydrocarbonproduction. Meanwhile, microseismic （MS） monitoring is one of the most effective approaches to evaluatesuch stimulation process. In this paper, the combined finite-discrete element method （FDEM） isadopted to numerically simulate HF and associated MS. Several post-processing tools, includingfrequency-magnitude distribution （b-value）, fractal dimension （D-value）, and seismic events clustering,are utilized to interpret numerical results. A non-parametric clustering algorithm designed specificallyfor FDEM is used to reduce the mesh dependency and extract more realistic seismic information.Simulation results indicated that at the local scale, the HF process tends to propagate following the rockmass discontinuities; while at the reservoir scale, it tends to develop in the direction parallel to themaximum in-situ stress. 2014 Institute of Rock and Soil Mechanics, Chinese Academy of Sciences. Production and hosting byElsevier B.V. All rights reserved.展开更多
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.展开更多
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.展开更多
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.展开更多
Compared with the middle-shallow shale gas reservoirs,deep shale formation often reveals lower fracture complexity,smaller effective stimulated volumes,lower conductivity,and faster declined when the same injection mo...Compared with the middle-shallow shale gas reservoirs,deep shale formation often reveals lower fracture complexity,smaller effective stimulated volumes,lower conductivity,and faster declined when the same injection model and technical parameters were used.Therefore,it is necessary to research a special injection model for the deep shale formation.In this paper,a series of laboratory stimulated experiments were carried to analyze the injection method on hydraulic fracturing initiation and propagation in the deep shale formation.The outcrop shale of Longmaxi formation in Sichuan Basin was selected and three different injection methods(constant,stepped-rise and cyclic progressive uplift pump rate)were used in these experiments.The test results showed that the stepped-rise pump rate is conductive to opening the pre-existing natural fracture and bedding planes,the fracture complexity of this injection method is the highest.Meanwhile,the fracturing pressure and total injection volume are the lowest.In contrast,the cyclic-uplift injection method can only form a transverse fracture,and the fracturing pressure is the highest.The research results have an important theoretical guiding value for deep shale hydraulic fracturing operation.展开更多
Injection pressure is one of the key parameters used in the design of shale gas reservoir stimulation.Higher injection pressure helps to initiate fractures,inject fracturing fluid and proppant smoothly and maximize th...Injection pressure is one of the key parameters used in the design of shale gas reservoir stimulation.Higher injection pressure helps to initiate fractures,inject fracturing fluid and proppant smoothly and maximize the stimulated reservoir volume as soon as possible.If the injection pressure is too high,however,the risk of casing deformation is increased significantly.In this paper,the numerical calculation format of the injection pressure safety window while ensuring casing integrity(i.e.,the maximum safety injection pressure)was proposed based on an example of an actual engineering project.Then,it was verified based on the actual situation of one shale gas well in Weiyuan.The numerical calculation format is as follows.First,a 3D finite element model of initial fine geostress field is established in the scale of block.Second,a primary submodel for introducing the asymmetric characteristics of reservoir stiffness to simulate the asymmetry of fracture distribution.Third,a secondary submodel containing the attributes of casing,cement sheath and reservoir material properties.Fourth,submodels are used to calculate the casing deformation generated by different injection loads and estimate the maximum allowable fracturing injection pressure(pj)while the safety requirement of casing deformation is satisfied.The calculation result of the case well shows that when the cementing quality is poor and fractures are distributed asymmetrically,the lateral and vertical displacements at the maximum displacement point of the casing under the injection pressure of 80 MPa are obviously lower than those under 90 MPa.According to the yield limit criterion of P-110 casing steel,the casing deformation in this case is elastic strain,so the injection pressure of 80 MPa is safe.In conclusion,this method has rational precision and accuracy,for its numerical result is consistent with the actual engineering phenomenon.展开更多
基金the financial support from National Natural Science Foundation of China(No.52474029)Strategic and Applied Scientific Research Project of PetroChina Company Limited(2023ZZ18,2023ZZ18YJ04).
文摘In the development of coalbed methane(CBM)reservoirs using multistage fractured horizontal wells,there often exist areas that are either repeatedly stimulated or completely unstimulated between fracturing stages,leading to suboptimal reservoir performance.Currently,there is no well-established method for accurately evaluating the effectiveness of such stimulation.This study introduces,for the first time,the concept of the Fracture Network Bridging Coefficient(FNBC)as a novel metric to assess stimulation performance.By quantitatively coupling the proportions of unstimulated and overstimulated volumes,the FNBC effectively characterizes the connectivity and efficiency of the fracture network.A background grid calibration method is developed to quantify the stage-controlled volume,effectively stimulated volume,unstimulated volume,and repeatedly stimulated volume among different stages of horizontal wells.Furthermore,an optimization model is constructed by taking the FNBC as the objective function and the fracturing injection rate and fluid volume as optimization variables.The Simultaneous Perturbation Stochastic Approximation(SPSA)algorithm is employed to iteratively perturb and optimize these variables,progressively improving the FNBC until the optimal displacement rate and fluid volume corresponding to the maximum FNBC are obtained.Field application in a typical CBM multistage fractured horizontal well in China demonstrates that the FNBC increased from 0.358 to 0.539(a 50.6% improvement),with the injection rate rising from 16 m^(3)/min to 24 m^(3)/min and the average fluid volume per stage increasing from 2490 m^(3) to 3192 m^(3),significantly enhancing the stimulation effectiveness.This research provides theoretical support for designing high-efficiency stimulation strategies in unconventional reservoirs under dynamic limits.
基金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.
基金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.
基金the financial support to the development of PNMM from the Laboratory of Rock Mechanics at école Polytechnique Fédérale de Lausanne (LMR-EPFL), Monash Universitythe National Natural Science Foundation of China (Grant No. 11802058)
文摘This review summarizes the development of particle-based numerical manifold method(PNMM)and its applications to rock dynamics.The fundamental principle of numerical manifold method(NMM)is first briefly introduced.Then,the history of the newly developed PNMM is given.Basic idea of PNMM and its simulation procedure are presented.Considering that PNMM could be regarded as an NMM-based model,a comparison of PNMM and NMM is discussed from several points of view in this paper.Besides,accomplished applications of PNMM to the dynamic rock fracturing are also reviewed.Finally,some recommendations are provided for the future work of PNMM.
基金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.
基金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.
基金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 Natural Sciences and Engineering Research Council of Canada through Discovery Grant 341275 (G. Grasselli) and Engage EGP 461019-13
文摘Hydraulic fracturing （HF） technique has been extensively used for the exploitation of unconventional oiland gas reservoirs. HF enhances the connectivity of less permeable oil and gas-bearing rock formationsby fluid injection, which creates an interconnected fracture network and increases the hydrocarbonproduction. Meanwhile, microseismic （MS） monitoring is one of the most effective approaches to evaluatesuch stimulation process. In this paper, the combined finite-discrete element method （FDEM） isadopted to numerically simulate HF and associated MS. Several post-processing tools, includingfrequency-magnitude distribution （b-value）, fractal dimension （D-value）, and seismic events clustering,are utilized to interpret numerical results. A non-parametric clustering algorithm designed specificallyfor FDEM is used to reduce the mesh dependency and extract more realistic seismic information.Simulation results indicated that at the local scale, the HF process tends to propagate following the rockmass discontinuities; while at the reservoir scale, it tends to develop in the direction parallel to themaximum in-situ stress. 2014 Institute of Rock and Soil Mechanics, Chinese Academy of Sciences. Production and hosting byElsevier B.V. All rights reserved.
基金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 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 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 National Science and Technology Major Project of China,2016ZX05060-004,2017ZX05036-003the National Natural Science Foundation of China,51574218Strategic Priority Research Program of the Chinese Academy of Sciences,grant number XDB10040200。
文摘Compared with the middle-shallow shale gas reservoirs,deep shale formation often reveals lower fracture complexity,smaller effective stimulated volumes,lower conductivity,and faster declined when the same injection model and technical parameters were used.Therefore,it is necessary to research a special injection model for the deep shale formation.In this paper,a series of laboratory stimulated experiments were carried to analyze the injection method on hydraulic fracturing initiation and propagation in the deep shale formation.The outcrop shale of Longmaxi formation in Sichuan Basin was selected and three different injection methods(constant,stepped-rise and cyclic progressive uplift pump rate)were used in these experiments.The test results showed that the stepped-rise pump rate is conductive to opening the pre-existing natural fracture and bedding planes,the fracture complexity of this injection method is the highest.Meanwhile,the fracturing pressure and total injection volume are the lowest.In contrast,the cyclic-uplift injection method can only form a transverse fracture,and the fracturing pressure is the highest.The research results have an important theoretical guiding value for deep shale hydraulic fracturing operation.
基金supported by the General Program of National Natural Science Foundation of China“Theoretical and Experimental Study on Key Mechanical Problems in Unconventional Natural Gas Exploitation”(No.:11272216).
文摘Injection pressure is one of the key parameters used in the design of shale gas reservoir stimulation.Higher injection pressure helps to initiate fractures,inject fracturing fluid and proppant smoothly and maximize the stimulated reservoir volume as soon as possible.If the injection pressure is too high,however,the risk of casing deformation is increased significantly.In this paper,the numerical calculation format of the injection pressure safety window while ensuring casing integrity(i.e.,the maximum safety injection pressure)was proposed based on an example of an actual engineering project.Then,it was verified based on the actual situation of one shale gas well in Weiyuan.The numerical calculation format is as follows.First,a 3D finite element model of initial fine geostress field is established in the scale of block.Second,a primary submodel for introducing the asymmetric characteristics of reservoir stiffness to simulate the asymmetry of fracture distribution.Third,a secondary submodel containing the attributes of casing,cement sheath and reservoir material properties.Fourth,submodels are used to calculate the casing deformation generated by different injection loads and estimate the maximum allowable fracturing injection pressure(pj)while the safety requirement of casing deformation is satisfied.The calculation result of the case well shows that when the cementing quality is poor and fractures are distributed asymmetrically,the lateral and vertical displacements at the maximum displacement point of the casing under the injection pressure of 80 MPa are obviously lower than those under 90 MPa.According to the yield limit criterion of P-110 casing steel,the casing deformation in this case is elastic strain,so the injection pressure of 80 MPa is safe.In conclusion,this method has rational precision and accuracy,for its numerical result is consistent with the actual engineering phenomenon.
