Multistage fracturing technology has been used to enhance tight hydrocarbon resource recovery.Determining the proper well spacing and fracturing strategy is crucial for generating a complex fracture network that facil...Multistage fracturing technology has been used to enhance tight hydrocarbon resource recovery.Determining the proper well spacing and fracturing strategy is crucial for generating a complex fracture network that facilitates oil and gas flow in reservoirs.The stress-shadow effect that occurs between multiple wells significantly affects the development of fracture networks in reservoirs.However,the quantification of the stress-shadow effect and its influence on fracture networks has not been satisfactorily resolved because of the difficulties in detecting and identifying fracture propagation and reorientation in reservoirs.In this study,based on the geological information from the Shengli oilfield,we applied a hybrid finite element-discrete element method to analyze engineering-scale three-dimensional fracture propagation and reorientation by altering well spacings and fracturing strategies.The results indicate that the fracturing area generated by the synchronous fracturing scheme is much smaller than those generated by the sequential and alternative schemes.An alternative hydrofracturing scheme is optimal with respect to fracturing area.The stress-blind area was defined to quantify the mechanical disturbance between adjacent wells.Our study improves the understanding of the effect of fracturing schemes on fracture networks and the impact of independent factors contributing to stress-shadow effects.展开更多
Complexities in mechanical behaviours of rock masses mainly stem from inherent discontinuities,which calls for advanced bolt-grouting techniques for stability enhancement.Understanding the mechanical properties of bol...Complexities in mechanical behaviours of rock masses mainly stem from inherent discontinuities,which calls for advanced bolt-grouting techniques for stability enhancement.Understanding the mechanical properties of bolt-grouted fractured rock mass(BGFR)and developing accurate prediction methods are crucial to optimize the BGFR support strategies.This paper establishes a new elastoplastic(E-P)model based on the orthotropic and the Mohr-Coulomb(M-C)plastic-yielding criteria.The elastic parameters of the model were derived through a meso-mechanical analysis of composite materials mechanics(CMM).Laboratory BGFR specimens were prepared and uniaxial compression test and variable-angle shear test considering different bolt arrangements were carried out to obtain the mechanical parameters of the specimens.Results showed that the anisotropy of BGFR mainly depends on the relative volume content of each component material in a certain direction.Moreover,the mechanical parameters deduced from the theory of composite materials which consider the short fibre effect are shown to be in good agreement with those determined by laboratory experiments,and the variation rules maintained good consistency.Last,a case study of a real tunnel project is provided to highlight the effectiveness,validity and robustness of the developed E-P model in prediction of stresses and deformations.展开更多
Under the policy background and advocacy of carbon capture,utilization,and storage(CCUS),CO_(2)-EOR has become a promising direction in the shale oil reservoir industry.The multi-scale pore structure distribution and ...Under the policy background and advocacy of carbon capture,utilization,and storage(CCUS),CO_(2)-EOR has become a promising direction in the shale oil reservoir industry.The multi-scale pore structure distribution and fracture structure lead to complex multiphase flow,comprehensively considering multiple mechanisms is crucial for development and CO_(2) storage in fractured shale reservoirs.In this paper,a multi-mechanism coupled model is developed by MATLAB.Compared to the traditional Eclipse300 and MATLAB Reservoir Simulation Toolbox(MRST),this model considers the impact of pore structure on fluid phase behavior by the modified Peng—Robinson equation of state(PR-EOS),and the effect simultaneously radiate to Maxwell—Stefan(M—S)diffusion,stress sensitivity,the nano-confinement(NC)effect.Moreover,a modified embedded discrete fracture model(EDFM)is used to model the complex fractures,which optimizes connection types and half-transmissibility calculation approaches between non-neighboring connections(NNCs).The full implicit equation adopts the finite volume method(FVM)and Newton—Raphson iteration for discretization and solution.The model verification with the Eclipse300 and MRST is satisfactory.The results show that the interaction between the mechanisms significantly affects the production performance and storage characteristics.The effect of molecular diffusion may be overestimated in oil-dominated(liquid-dominated)shale reservoirs.The well spacing and injection gas rate are the most crucial factors affecting the production by sensitivity analysis.Moreover,the potential gas invasion risk is mentioned.This model provides a reliable theoretical basis for CO_(2)-EOR and sequestration in shale oil reservoirs.展开更多
Due to inherent limits of data acquisition and geophysical data resolution, there are large uncertainties in the characterization of subsurface fractures. However, outcrop analogies can provide qualitative and quantit...Due to inherent limits of data acquisition and geophysical data resolution, there are large uncertainties in the characterization of subsurface fractures. However, outcrop analogies can provide qualitative and quantitative information on a large number of fractures, based on which the accuracy of subsurface fracture characterization can be improved. Here we take the tectonic fracture modeling of an ultra-low permeability sandstone reservoir based on an outcrop analogy, a case study of the Chang6t~ Formation of the Upper Triassic Yanchang Group of the Wangyao Oilfield in the Ordos Basin of China. An outcrop at the edge of the basin is a suitable analog for the reservoir, but the prerequisite is that they must have equivalent previous stress fields, similar final structural characteristics, relative timing and an identical depositional environment and diagenesis. The relationship among fracture density, rock type and bed thickness based on the outcrop is one of the most important fracture distribution models, and can be used to interpret fracture density in individual wells quantitatively. Fracture orientation, dip, geometry and scale, also should be described and measured in the outcrop, and can be used together with structure restoration and single well fracture density interpretation to guide fracture intensity prediction on bed surfaces and to constrain the construction of the 3D fracture geometry model of the subsurface reservoir. The application of the above principles shows the outcrop-based tectonic fracture models of the target ultra-low permeability sandstone reservoir are consistent with fractures inferred from microseismic interpretation and tracer tests. This illustrated that the fracture modeling based on the outcrop analogy is reliable and can reduce the uncertainty in stochastic fracture modeling.展开更多
Fractured reservoirs are an important target for oil and gas exploration in the Tarim Basin and the prediction of this type of reservoir is challenging.Due to the complicated fracture system in the Tarim Basin,the con...Fractured reservoirs are an important target for oil and gas exploration in the Tarim Basin and the prediction of this type of reservoir is challenging.Due to the complicated fracture system in the Tarim Basin,the conventional AVO inversion method based on HTI theory to predict fracture development will result in some errors.Thus,an integrated research concept for fractured reservoir prediction is put forward in this paper.Seismic modeling plays a bridging role in this concept,and the establishment of an anisotropic fracture model by Discrete Fracture Network （DFN） is the key part.Because the fracture system in the Tarim Basin shows complex anisotropic characteristics,it is vital to build an effective anisotropic model.Based on geological,well logging and seismic data,an effective anisotropic model of complex fracture systems can be set up with the DFN method.The effective elastic coefficients,and the input data for seismic modeling can be calculated.Then seismic modeling based on this model is performed,and the seismic response characteristics are analyzed.The modeling results can be used in the following AVO inversion for fracture detection.展开更多
To investigate the height growth of multi-cluster fractures during variable fluid-viscosity fracturing in a layered shale oil reservoir,a two-dimensional finite element method(FEM)-discrete fracture network(DFN)model ...To investigate the height growth of multi-cluster fractures during variable fluid-viscosity fracturing in a layered shale oil reservoir,a two-dimensional finite element method(FEM)-discrete fracture network(DFN)model coupled with flow,stress and damage is proposed.A traction-separation law is used to describe the mixed-mode response of the damaged adhesive fractures,and the cubic law is used to describe the fluid flow within the fractures.The rock deformation is controlled by the in-situ stress,fracture cohesion and fluid pressure on the hydraulic fracture surface.The coupled finite element equations are solved by the explicit time difference method.The effects of the fracturing treatment parameters including fluid viscosity,pumping rate and cluster spacing on the geometries of multifractures are investigated.The results show that variable fluid-viscosity injection can improve the complexity of the fracture network and height of the main fractures simultaneously.The pumping rate of15 m^(3)/min,variable fluid-viscosity of 3-9-21-36-45 mPa s with a cluster spacing of 7.5 m is the ideal treatment strategy.The field application shows that the peak daily production of the application well with the optimized injection procedu re of variable fluid-viscosity fracturing is 171 tons(about 2.85 times that of the adjacent well),which is the highest daily production record of a single shale oil well in China,marking a strategic breakthrough of commercial shale oil production in the Jiyang Depression,Shengli Oilfield.The variable fluid-viscosity fracturing technique is proved to be very effective for improving shale oil production.展开更多
Since natural fractures are often non-equidimensional,the circular disc model still has great limitations.By contrast,the elliptical disc model is more applicable to representing natural fractures,especially for slend...Since natural fractures are often non-equidimensional,the circular disc model still has great limitations.By contrast,the elliptical disc model is more applicable to representing natural fractures,especially for slender ones.This paper developed a universal elliptical disc(UED)model by incorporating the center point,size,and azimuth of fractures as variables.Specifically,with respect to the azimuth of elliptical fractures in three-dimensional(3D)space,we proposed a paradigm to construct its probability density function(PDF)by coupling the orientation and rotation angle of long axis based on three coordinate transformations.To illustrate the construction process of the PDF of the fracture azimuth,we took the orientation following the Fisher distribution and the rotation angle following Von Mises distribution as an example.A rock slope is used to show the use of the developed UED model,and the 3D DFNs for the slope rock mass are generated by Monte Carlo simulation.In addition,the DFNs for the rock mass are also generated based on the existing circular disc model and non-universal elliptical disc model.The comparison results from the three models clearly illustrate the superiority of the UED model over the existing circular and non-universal elliptical disc models.展开更多
Based on break characteristics of roofs in fully mechanized top-coal mining of thick shallow coal seams, a fracture mechanics model was built, and the criterion of crack propagation in the main roof was derived using ...Based on break characteristics of roofs in fully mechanized top-coal mining of thick shallow coal seams, a fracture mechanics model was built, and the criterion of crack propagation in the main roof was derived using the fracture mechanics theory. The relationships between the fracture length of the roof and the working resistance of the supports were discovered, and the correlations between the load on the overlying strata and the ratio of the crack's length to the thickness of the roof were obtained. Using a working face of Jindi Coal Mine, Xing county Shanxi province as an example, the relationships between the fracture length of the roof and the working resistance of the supports were analysed in detail. The results give a design basis in hydraulic top coal caving supports, which could provide useful references in the practical application. On-site experiment proves that the periodic weighting step interval of the caving face is 15–16 m, which is basically consistent with the theoretical analysis results, and indicates that the mechanized caving hydraulic support is capable of meeting the support requirements in the mining of a super-thick but shallowly buried coal seam.展开更多
The main objective is to optimize the development of shale gas-rich areas by predicting seismic sweet spot parameters in shale reservoirs. We systematically assessed the fracture development, fracture gas content, and...The main objective is to optimize the development of shale gas-rich areas by predicting seismic sweet spot parameters in shale reservoirs. We systematically assessed the fracture development, fracture gas content, and rock brittleness in fractured gas-bearing shale reservoirs. To better characterize gas-bearing shale reservoirs with tilted fractures, we optimized the petrophysical modeling based on the equivalent medium theory. Based on the advantages of shale petrophysical modeling, we not only considered the brittle mineral fraction but also the combined effect of shale porosity, gas saturation, and total organic carbon(TOC) when optimizing the brittleness index. Due to fractures generally functioning as essential channels for fluid storage and movement, fracture density and fracture fluid identification factors are critical geophysical parameters for fractured reservoir prediction. We defined a new fracture gas indication factor(GFI) to detect fracture-effective gas content. A new linear PP-wave reflection coefficient equation for a tilted transversely isotropic(TTI) medium was rederived, realizing the direct prediction of anisotropic fracture parameters and the isotropic elasticity parameters from offset vector tile(OVT)-domain seismic data. Synthetic seismic data experiments demonstrated that the inversion algorithm based on the L_P quasinorm sparsity constraint and the split-component inversion strategy exhibits high stability and noise resistance. Finally, we applied our new prediction method to evaluate fractured gas-bearing shale reservoirs in the Sichuan Basin of China, demonstrating its effectiveness.展开更多
Heterotopic ossification(HO)is a consequence of traumatic bone and tissue damage,which occurs in 65%of military casualties with blast-associated amputations.However,the mechanisms behind blast-induced HO remain unclea...Heterotopic ossification(HO)is a consequence of traumatic bone and tissue damage,which occurs in 65%of military casualties with blast-associated amputations.However,the mechanisms behind blast-induced HO remain unclear.Animal models are used to study blast-induced HO,but developing such models is challenging,particularly in how to use a pure blast wave(primary blast)to induce limb fracture that then requires an amputation.Several studies,including our recent study,have developed platforms to induce limb fractures in rats with blast loading or a mixture of blast and impact loading.However,these models are limited by the survivability of the animal and repeatability of the model.In this study,we developed an improved platform,aiming to improve the animal's survivability and injury repeatability as well as focusing on primary blast only.The platform exposed only one limb of the rat to a blast wave while providing proper protection to the rest of the rat's body.We obtained very consistent fracture outcome in the tibia(location and pattern)in cadaveric rats with a large range of size and weight.Importantly,the rats did not obviously move during the test,where movement is a potential cause of uncontrolled injury.We further conducted parametric studies by varying the features of the design of the platform.These factors,such as how the limb is fixed and how the cavity through which the limb is placed is sealed,significantly affect the resulting injury.This platform and test setups enable well-controlled limb fracture induced directly by pure blast wave,which is the fundamental step towards a complete in vivo animal model for blast-induced HO induced by primary blast alone,excluding secondary and tertiary blast injury.In addition,the platform design and the findings presented here,particularly regarding the proper protection of the animal,have implications for future studies investigating localized blast injuries,such as blast induced brain and lung injuries.展开更多
The continental shale reservoirs of Jurassic Lianggaoshan Formation in Sichuan Basin contain thin lamina,which is characterized by strong plasticity and developed longitudinal shell limestone interlayer.To improve the...The continental shale reservoirs of Jurassic Lianggaoshan Formation in Sichuan Basin contain thin lamina,which is characterized by strong plasticity and developed longitudinal shell limestone interlayer.To improve the production efficiency of reservoirs by multi-cluster fracturing,it is necessary to consider the unbalanced propagation of hydraulic fractures and the penetration effect of fractures.This paper constructed a numerical model of multi-fracture propagation and penetration based on the finite element coupling cohesive zone method;considering the construction cluster spacing,pump rate,lamina strength and other parameters studied the influencing factors of multi-cluster fracture interaction propagation;combined with AE energy data and fracture mode reconstruction method,quantitatively characterized the comprehensive impact of the strength of thin interlayer rock interfaces on the initiation and propagation of fractures that penetrate layers,and accurately predicted the propagation pattern of hydraulic fractures through laminated shale oil reservoirs.Simulation results revealed that in the process of multi-cluster fracturing,the proportion of shear damage is low,and mainly occurs in bedding fractures activated by outer fractures.Reducing the cluster spacing enhances the fracture system's penetration ability,though it lowers the activation efficiency of lamina.The high plasticity of the limestone interlayer may impact the vertical propagation distance of the main fracture.Improving the interface strength is beneficial to the reconstruction of the fracture height,but the interface communication effect is limited.Reasonable selection of layers with moderate lamina strength for fracturing stimulation,increasing the pump rate during fracturing and setting the cluster spacing reasonably are beneficial to improve the effect of reservoir stimulation.展开更多
The relationship between hydrogen concentration and crack tip stress, strain field, hydrogen diffusion and internal pressure respectively in the crack tip process zone was investigated, and the length of the crack tip...The relationship between hydrogen concentration and crack tip stress, strain field, hydrogen diffusion and internal pressure respectively in the crack tip process zone was investigated, and the length of the crack tip process zone of hydrogen-induced cracking （HIC） was determined. Based on the mechanism of fracture of micro-crack nucleation, a dislocation model was presented for the fracture criteria of HIC. The influence factors of pipeline tube fracture ductile KISCC in the presence of hydrogen was analyzed, and the critical pressure bearing capability of a pipeline with hydrogeninduced cracking and the critical J-integrity （JISCC） were calculated, which is very important for pipeline safety.展开更多
The paper presents a novel hydraulic fracturing model for the characterization and simulation of the complex fracture network in shale gas reservoirs. We go beyond the existing method that uses planar or orthogonal co...The paper presents a novel hydraulic fracturing model for the characterization and simulation of the complex fracture network in shale gas reservoirs. We go beyond the existing method that uses planar or orthogonal conjugate fractures for representing the ＇＇complexity＇＇ of the network. Bifurcation of fractures is performed utilizing the Lindenmayer system based on fractal geometry to describe the fracture propagation pattern, density and network connectivity. Four controlling parameters are proposed to describe the details of complex fractures and stimulated reservoir volume（SRV）. The results show that due to the multilevel feature of fractal fractures, the model could provide a simple method for contributing reservoir volume calibration. The primary-and second-stage fracture networks across the overall SRV are the main contributions to the production, while the induced fracture network just contributes another 20% in the late producing period. We also conduct simulation with respect to different refracturing cases and find that increasing the complexity of the fracture network provides better performance than only enhancing the fracture conductivity.展开更多
Accurate fluid flow simulation in geologically complex reservoirs is of particular importance in construction of reservoir simulators.General approaches in naturally fractured reservoir simulation involve use of unstr...Accurate fluid flow simulation in geologically complex reservoirs is of particular importance in construction of reservoir simulators.General approaches in naturally fractured reservoir simulation involve use of unstructured grids or a structured grid coupled with locally unstructured grids and discrete fracture models.These methods suffer from drawbacks such as lack of flexibility and of ease of updating.In this study,I combined fracture modeling by elastic gridding which improves flexibility,especially in complex reservoirs.The proposed model revises conventional modeling fractures by hard rigid planes that do not change through production.This is a dubious assumption,especially in reservoirs with a high production rate in the beginning.The proposed elastic fracture modeling considers changes in fracture properties,shape and aperture through the simulation.This strategy is only reliable for naturally fractured reservoirs with high fracture permeability and less permeable matrix and parallel fractures with less cross-connections.Comparison of elastic fracture modeling results with conventional modeling showed that these assumptions will cause production pressure to enlarge fracture apertures and change fracture shapes,which consequently results in lower production compared with what was previously assumed.It is concluded that an elastic gridded model could better simulate reservoir performance.展开更多
Structure plane is one of the important factors affecting the stability and failure mode of rock mass engineering.Rock mass structure characterization is the basic work of rock mechanics research and the important con...Structure plane is one of the important factors affecting the stability and failure mode of rock mass engineering.Rock mass structure characterization is the basic work of rock mechanics research and the important content of numerical simulation.A new 3-dimensional rough discrete fracture network(RDFN3D)model and its modeling method based on the Weierstrass-Mandelbrot(W-M)function were presented in this paper.The RDFN3D model,which improves and unifies the modelling methods for the complex structural planes,has been realized.The influence of fractal dimension,amplitude,and surface precision on the modeling parameters of RDFN3D was discussed.The reasonable W-M parameters suitable for the roughness coefficient of JRC were proposed,and the relationship between the mathematical model and the joint characterization was established.The RDFN3D together with the smooth 3-dimensional discrete fracture network(DFN3D)models were successfully exported to the drawing exchange format,which will provide a wide application in numerous numerical simulation codes including both the continuous and discontinuous methods.The numerical models were discussed using the COMSOL Multiphysics code and the 3-dimensional particle flow code,respectively.The reliability of the RDFN3D model was preliminarily discussed and analyzed.The roughness and spatial connectivity of the fracture networks have a dominant effect on the fluid flow patterns.The research results can provide a new geological model and analysis model for numerical simulation and engineering analysis of jointed rock mass.展开更多
In practical development of unconventional reservoirs,fracture networks are a highly conductive transport media for subsurface fluid flow.Therefore,it is crucial to clearly determine the fracture properties used in pr...In practical development of unconventional reservoirs,fracture networks are a highly conductive transport media for subsurface fluid flow.Therefore,it is crucial to clearly determine the fracture properties used in production forecast.However,it is different to calibrate the properties of fracture networks because it is an inverse problem with multi-patterns and highcomplexity of fracture distribution and inherent defect of multiplicity of solution.In this paper,in order to solve the problem,the complex fracture model is divided into two sub-systems,namely"Pattern A"and"Pattern B."In addition,the generation method is grouped into two categories.Firstly,we construct each sub-system based on the probability density function of the fracture properties.Secondly,we recombine the sub-systems into an integral complex fracture system.Based on the generation mechanism,the estimation of the complex fracture from dynamic performance and observation data can be solved as an inverse problem.In this study,the Bayesian formulation is used to quantify the uncertainty of fracture properties.To minimize observation data misfit immediately as it occurs,we optimize the updated properties by a simultaneous perturbation stochastic algorithm which requires only two measurements of the loss function.In numerical experiments,we firstly visualize that small-scale fractures significantly contribute to the flow simulation.Then,we demonstrate the suitability and effectiveness of the Bayesian formulation for calibrating the complex fracture model in the following simulation.展开更多
Understanding the impact of mining disturbances and creep deformation on the macroscopic deformation and the microscopic pore and fracture structures(MPFS)of coal is paramount for ensuring the secure extraction of coa...Understanding the impact of mining disturbances and creep deformation on the macroscopic deformation and the microscopic pore and fracture structures(MPFS)of coal is paramount for ensuring the secure extraction of coal resources.This study conducts cyclic loading-unloading and creep experiments on coal using a low-field nuclear magnetic resonance(NMR)experimental apparatus which is equipped with mechanical loading units,enabling real-time monitoring the T2spectrum.The experiments indicated that cyclic loading-unloading stress paths initiate internal damage within coal samples.Under identical creep stress conditions,coal samples with more initial damages had more substantial instantaneous deformation and creep deformation during the creep process.After undergoing nearly 35 h of staged creep,the total strains for coal samples CC01,CC02,and CC03 reach 2.160%,2.261%,and 2.282%,respectively.In the creep stage,the peak area ratio of seepage pores and microfractures(SPM)gradually diminishes.A higher degree of initial damage leads to a more pronounced compaction trend in the SPM of coal samples.Considering the porosity evolution of SPM during the creep process,this study proposes a novel fractional derivative model for the porosity evolution of SPM.The efficacy of the proposed model in predicting porosity evolution of SPM is substantiated through experimental validation.Furthermore,an analysis of the impact mechanisms on key parameters in the model was carried out.展开更多
A new numerical model is presented to simulate fracture initiation and propagation in geological structures. This model is based on the recent amalgamation of established failure and fracture mechanics theory, which h...A new numerical model is presented to simulate fracture initiation and propagation in geological structures. This model is based on the recent amalgamation of established failure and fracture mechanics theory, which has been implemented to the finite difference FLAC code as a constitutive FISH userdefined-model. Validation of the model has been studied on the basis of comparing the transitional failure modes in rock. It is shown that the model is capable of accurately simulating fracture distributions over entire brittle to ductile rock phases. The application of the model during longwall retreat simulation highlighted several caving characteristics relevant to varying geological condition. The distribution and behaviour of modelled fractures were both realistic and shown to provide an enhanced post failure analysis to geological structures in FLAC. Moreover, the model introduces new potential insight towards the failure analysis of more complicated problems. This is best suited towards improving safety and efficiency in mines through the prediction of various key fractures and caving characteristics of geological structures.展开更多
Investigation of paper cutting process is vital for the design of cutting tools,but the fracture mechanism of paper cutting is still unclear.Here,we focus on the cutting process of paper,including the key parameters o...Investigation of paper cutting process is vital for the design of cutting tools,but the fracture mechanism of paper cutting is still unclear.Here,we focus on the cutting process of paper,including the key parameters of cohesive zone model(CZM)for the orthotropic paper,to simulate the shear fracture process.Firstly,the material constants of the orthotropic paper are determined by longitudinal and transverse tensile test.Secondly,based on the tensile stressstrain curves,combined with damage theory and numerical simulations,the key parameters of the CZM for the orthotropic paper are obtained.Finally,a model III fracture is simulated to verify the accuracy of the model.Results show that the load-displacement curves obtained by the simulation is consistent with the test results.展开更多
The generation method of three-dimensional fractal discrete fracture network(FDFN)based on multiplicative cascade process was developed.The complex multi-scale fracture system in shale after fracturing was characteriz...The generation method of three-dimensional fractal discrete fracture network(FDFN)based on multiplicative cascade process was developed.The complex multi-scale fracture system in shale after fracturing was characterized by coupling the artificial fracture model and the natural fracture model.Based on an assisted history matching(AHM)using multiple-proxy-based Markov chain Monte Carlo algorithm(MCMC),an embedded discrete fracture modeling(EDFM)incorporated with reservoir simulator was used to predict productivity of shale gas well.When using the natural fracture generation method,the distribution of natural fracture network can be controlled by fractal parameters,and the natural fracture network generated coupling with artificial fractures can characterize the complex system of different-scale fractures in shale after fracturing.The EDFM,with fewer grids and less computation time consumption,can characterize the attributes of natural fractures and artificial fractures flexibly,and simulate the details of mass transfer between matrix cells and fractures while reducing computation significantly.The combination of AMH and EDFM can lower the uncertainty of reservoir and fracture parameters,and realize effective inversion of key reservoir and fracture parameters and the productivity forecast of shale gas wells.Application demonstrates the results from the proposed productivity prediction model integrating FDFN,EDFM and AHM have high credibility.展开更多
基金supported in part by the National Key Research and Development Project of China(No.2022YFC3004602)in part by the National Natural Science Foundation of China(Nos.52121003 and 52342403).
文摘Multistage fracturing technology has been used to enhance tight hydrocarbon resource recovery.Determining the proper well spacing and fracturing strategy is crucial for generating a complex fracture network that facilitates oil and gas flow in reservoirs.The stress-shadow effect that occurs between multiple wells significantly affects the development of fracture networks in reservoirs.However,the quantification of the stress-shadow effect and its influence on fracture networks has not been satisfactorily resolved because of the difficulties in detecting and identifying fracture propagation and reorientation in reservoirs.In this study,based on the geological information from the Shengli oilfield,we applied a hybrid finite element-discrete element method to analyze engineering-scale three-dimensional fracture propagation and reorientation by altering well spacings and fracturing strategies.The results indicate that the fracturing area generated by the synchronous fracturing scheme is much smaller than those generated by the sequential and alternative schemes.An alternative hydrofracturing scheme is optimal with respect to fracturing area.The stress-blind area was defined to quantify the mechanical disturbance between adjacent wells.Our study improves the understanding of the effect of fracturing schemes on fracture networks and the impact of independent factors contributing to stress-shadow effects.
基金funded by the National Key Research and Development Plan(No.2022YFC3203200)Department of Science and Technology of Guangdong Province(No.2021ZT09G087)the National Natural Science Foundation Project of China(No.42167025).
文摘Complexities in mechanical behaviours of rock masses mainly stem from inherent discontinuities,which calls for advanced bolt-grouting techniques for stability enhancement.Understanding the mechanical properties of bolt-grouted fractured rock mass(BGFR)and developing accurate prediction methods are crucial to optimize the BGFR support strategies.This paper establishes a new elastoplastic(E-P)model based on the orthotropic and the Mohr-Coulomb(M-C)plastic-yielding criteria.The elastic parameters of the model were derived through a meso-mechanical analysis of composite materials mechanics(CMM).Laboratory BGFR specimens were prepared and uniaxial compression test and variable-angle shear test considering different bolt arrangements were carried out to obtain the mechanical parameters of the specimens.Results showed that the anisotropy of BGFR mainly depends on the relative volume content of each component material in a certain direction.Moreover,the mechanical parameters deduced from the theory of composite materials which consider the short fibre effect are shown to be in good agreement with those determined by laboratory experiments,and the variation rules maintained good consistency.Last,a case study of a real tunnel project is provided to highlight the effectiveness,validity and robustness of the developed E-P model in prediction of stresses and deformations.
基金supported by the National Natural Science Foundation of China(No.52174038 and No.52004307)China Petroleum Science and Technology Project-Major Project-Research on Tight Oil-Shale Oil Reservoir Engineering Methods and Key Technologies in Ordos Basin(No.ZLZX2020-02-04)Science Foundation of China University of Petroleum,Beijing(No.2462018YJRC015)。
文摘Under the policy background and advocacy of carbon capture,utilization,and storage(CCUS),CO_(2)-EOR has become a promising direction in the shale oil reservoir industry.The multi-scale pore structure distribution and fracture structure lead to complex multiphase flow,comprehensively considering multiple mechanisms is crucial for development and CO_(2) storage in fractured shale reservoirs.In this paper,a multi-mechanism coupled model is developed by MATLAB.Compared to the traditional Eclipse300 and MATLAB Reservoir Simulation Toolbox(MRST),this model considers the impact of pore structure on fluid phase behavior by the modified Peng—Robinson equation of state(PR-EOS),and the effect simultaneously radiate to Maxwell—Stefan(M—S)diffusion,stress sensitivity,the nano-confinement(NC)effect.Moreover,a modified embedded discrete fracture model(EDFM)is used to model the complex fractures,which optimizes connection types and half-transmissibility calculation approaches between non-neighboring connections(NNCs).The full implicit equation adopts the finite volume method(FVM)and Newton—Raphson iteration for discretization and solution.The model verification with the Eclipse300 and MRST is satisfactory.The results show that the interaction between the mechanisms significantly affects the production performance and storage characteristics.The effect of molecular diffusion may be overestimated in oil-dominated(liquid-dominated)shale reservoirs.The well spacing and injection gas rate are the most crucial factors affecting the production by sensitivity analysis.Moreover,the potential gas invasion risk is mentioned.This model provides a reliable theoretical basis for CO_(2)-EOR and sequestration in shale oil reservoirs.
基金supported by Open Fund (PLC201203) of State Key Laboratory of Oil and Gas Reservoir Geology and Exploitation (Chengdu University of Technology)Major Project of Education Department in Sichuan Province (13ZA0177)
文摘Due to inherent limits of data acquisition and geophysical data resolution, there are large uncertainties in the characterization of subsurface fractures. However, outcrop analogies can provide qualitative and quantitative information on a large number of fractures, based on which the accuracy of subsurface fracture characterization can be improved. Here we take the tectonic fracture modeling of an ultra-low permeability sandstone reservoir based on an outcrop analogy, a case study of the Chang6t~ Formation of the Upper Triassic Yanchang Group of the Wangyao Oilfield in the Ordos Basin of China. An outcrop at the edge of the basin is a suitable analog for the reservoir, but the prerequisite is that they must have equivalent previous stress fields, similar final structural characteristics, relative timing and an identical depositional environment and diagenesis. The relationship among fracture density, rock type and bed thickness based on the outcrop is one of the most important fracture distribution models, and can be used to interpret fracture density in individual wells quantitatively. Fracture orientation, dip, geometry and scale, also should be described and measured in the outcrop, and can be used together with structure restoration and single well fracture density interpretation to guide fracture intensity prediction on bed surfaces and to constrain the construction of the 3D fracture geometry model of the subsurface reservoir. The application of the above principles shows the outcrop-based tectonic fracture models of the target ultra-low permeability sandstone reservoir are consistent with fractures inferred from microseismic interpretation and tracer tests. This illustrated that the fracture modeling based on the outcrop analogy is reliable and can reduce the uncertainty in stochastic fracture modeling.
基金co-supported by the National Basic Research Program of China(Grant No.2011CB201103)the National Science and Technology Major Project(GrantNo.2011ZX05004003)
文摘Fractured reservoirs are an important target for oil and gas exploration in the Tarim Basin and the prediction of this type of reservoir is challenging.Due to the complicated fracture system in the Tarim Basin,the conventional AVO inversion method based on HTI theory to predict fracture development will result in some errors.Thus,an integrated research concept for fractured reservoir prediction is put forward in this paper.Seismic modeling plays a bridging role in this concept,and the establishment of an anisotropic fracture model by Discrete Fracture Network （DFN） is the key part.Because the fracture system in the Tarim Basin shows complex anisotropic characteristics,it is vital to build an effective anisotropic model.Based on geological,well logging and seismic data,an effective anisotropic model of complex fracture systems can be set up with the DFN method.The effective elastic coefficients,and the input data for seismic modeling can be calculated.Then seismic modeling based on this model is performed,and the seismic response characteristics are analyzed.The modeling results can be used in the following AVO inversion for fracture detection.
基金funded by the National Natural Science Foundation of China(Nos.52192622,51874253,U20A202)
文摘To investigate the height growth of multi-cluster fractures during variable fluid-viscosity fracturing in a layered shale oil reservoir,a two-dimensional finite element method(FEM)-discrete fracture network(DFN)model coupled with flow,stress and damage is proposed.A traction-separation law is used to describe the mixed-mode response of the damaged adhesive fractures,and the cubic law is used to describe the fluid flow within the fractures.The rock deformation is controlled by the in-situ stress,fracture cohesion and fluid pressure on the hydraulic fracture surface.The coupled finite element equations are solved by the explicit time difference method.The effects of the fracturing treatment parameters including fluid viscosity,pumping rate and cluster spacing on the geometries of multifractures are investigated.The results show that variable fluid-viscosity injection can improve the complexity of the fracture network and height of the main fractures simultaneously.The pumping rate of15 m^(3)/min,variable fluid-viscosity of 3-9-21-36-45 mPa s with a cluster spacing of 7.5 m is the ideal treatment strategy.The field application shows that the peak daily production of the application well with the optimized injection procedu re of variable fluid-viscosity fracturing is 171 tons(about 2.85 times that of the adjacent well),which is the highest daily production record of a single shale oil well in China,marking a strategic breakthrough of commercial shale oil production in the Jiyang Depression,Shengli Oilfield.The variable fluid-viscosity fracturing technique is proved to be very effective for improving shale oil production.
基金funded by the National Natural Science s of China(No.41972264)the Key R&D Project of Zhejiang Province(No.2021C03159)the Field Scientific Observation&Research Station of Geological Hazard in Zhejiang,Ministry of Natural Resources,China(No.ZJDZGCZ-2021)。
文摘Since natural fractures are often non-equidimensional,the circular disc model still has great limitations.By contrast,the elliptical disc model is more applicable to representing natural fractures,especially for slender ones.This paper developed a universal elliptical disc(UED)model by incorporating the center point,size,and azimuth of fractures as variables.Specifically,with respect to the azimuth of elliptical fractures in three-dimensional(3D)space,we proposed a paradigm to construct its probability density function(PDF)by coupling the orientation and rotation angle of long axis based on three coordinate transformations.To illustrate the construction process of the PDF of the fracture azimuth,we took the orientation following the Fisher distribution and the rotation angle following Von Mises distribution as an example.A rock slope is used to show the use of the developed UED model,and the 3D DFNs for the slope rock mass are generated by Monte Carlo simulation.In addition,the DFNs for the rock mass are also generated based on the existing circular disc model and non-universal elliptical disc model.The comparison results from the three models clearly illustrate the superiority of the UED model over the existing circular and non-universal elliptical disc models.
文摘Based on break characteristics of roofs in fully mechanized top-coal mining of thick shallow coal seams, a fracture mechanics model was built, and the criterion of crack propagation in the main roof was derived using the fracture mechanics theory. The relationships between the fracture length of the roof and the working resistance of the supports were discovered, and the correlations between the load on the overlying strata and the ratio of the crack's length to the thickness of the roof were obtained. Using a working face of Jindi Coal Mine, Xing county Shanxi province as an example, the relationships between the fracture length of the roof and the working resistance of the supports were analysed in detail. The results give a design basis in hydraulic top coal caving supports, which could provide useful references in the practical application. On-site experiment proves that the periodic weighting step interval of the caving face is 15–16 m, which is basically consistent with the theoretical analysis results, and indicates that the mechanized caving hydraulic support is capable of meeting the support requirements in the mining of a super-thick but shallowly buried coal seam.
基金financially supported by the Sichuan Science and Technology Program (Grant No. 2023ZYD0158)the National Natural Science Foundation of China (Grant Nos. 42304147 and 42304076)。
文摘The main objective is to optimize the development of shale gas-rich areas by predicting seismic sweet spot parameters in shale reservoirs. We systematically assessed the fracture development, fracture gas content, and rock brittleness in fractured gas-bearing shale reservoirs. To better characterize gas-bearing shale reservoirs with tilted fractures, we optimized the petrophysical modeling based on the equivalent medium theory. Based on the advantages of shale petrophysical modeling, we not only considered the brittle mineral fraction but also the combined effect of shale porosity, gas saturation, and total organic carbon(TOC) when optimizing the brittleness index. Due to fractures generally functioning as essential channels for fluid storage and movement, fracture density and fracture fluid identification factors are critical geophysical parameters for fractured reservoir prediction. We defined a new fracture gas indication factor(GFI) to detect fracture-effective gas content. A new linear PP-wave reflection coefficient equation for a tilted transversely isotropic(TTI) medium was rederived, realizing the direct prediction of anisotropic fracture parameters and the isotropic elasticity parameters from offset vector tile(OVT)-domain seismic data. Synthetic seismic data experiments demonstrated that the inversion algorithm based on the L_P quasinorm sparsity constraint and the split-component inversion strategy exhibits high stability and noise resistance. Finally, we applied our new prediction method to evaluate fractured gas-bearing shale reservoirs in the Sichuan Basin of China, demonstrating its effectiveness.
基金the auspices of the Royal British Legion Centre for Blast Injury Studies at Imperial College Londonthe financial support of the Royal British Legion。
文摘Heterotopic ossification(HO)is a consequence of traumatic bone and tissue damage,which occurs in 65%of military casualties with blast-associated amputations.However,the mechanisms behind blast-induced HO remain unclear.Animal models are used to study blast-induced HO,but developing such models is challenging,particularly in how to use a pure blast wave(primary blast)to induce limb fracture that then requires an amputation.Several studies,including our recent study,have developed platforms to induce limb fractures in rats with blast loading or a mixture of blast and impact loading.However,these models are limited by the survivability of the animal and repeatability of the model.In this study,we developed an improved platform,aiming to improve the animal's survivability and injury repeatability as well as focusing on primary blast only.The platform exposed only one limb of the rat to a blast wave while providing proper protection to the rest of the rat's body.We obtained very consistent fracture outcome in the tibia(location and pattern)in cadaveric rats with a large range of size and weight.Importantly,the rats did not obviously move during the test,where movement is a potential cause of uncontrolled injury.We further conducted parametric studies by varying the features of the design of the platform.These factors,such as how the limb is fixed and how the cavity through which the limb is placed is sealed,significantly affect the resulting injury.This platform and test setups enable well-controlled limb fracture induced directly by pure blast wave,which is the fundamental step towards a complete in vivo animal model for blast-induced HO induced by primary blast alone,excluding secondary and tertiary blast injury.In addition,the platform design and the findings presented here,particularly regarding the proper protection of the animal,have implications for future studies investigating localized blast injuries,such as blast induced brain and lung injuries.
基金financial support by the National Key Research and Development Program of China (No.2022YFE0129800)the National Natural Science Foundation of China (No.52074311)。
文摘The continental shale reservoirs of Jurassic Lianggaoshan Formation in Sichuan Basin contain thin lamina,which is characterized by strong plasticity and developed longitudinal shell limestone interlayer.To improve the production efficiency of reservoirs by multi-cluster fracturing,it is necessary to consider the unbalanced propagation of hydraulic fractures and the penetration effect of fractures.This paper constructed a numerical model of multi-fracture propagation and penetration based on the finite element coupling cohesive zone method;considering the construction cluster spacing,pump rate,lamina strength and other parameters studied the influencing factors of multi-cluster fracture interaction propagation;combined with AE energy data and fracture mode reconstruction method,quantitatively characterized the comprehensive impact of the strength of thin interlayer rock interfaces on the initiation and propagation of fractures that penetrate layers,and accurately predicted the propagation pattern of hydraulic fractures through laminated shale oil reservoirs.Simulation results revealed that in the process of multi-cluster fracturing,the proportion of shear damage is low,and mainly occurs in bedding fractures activated by outer fractures.Reducing the cluster spacing enhances the fracture system's penetration ability,though it lowers the activation efficiency of lamina.The high plasticity of the limestone interlayer may impact the vertical propagation distance of the main fracture.Improving the interface strength is beneficial to the reconstruction of the fracture height,but the interface communication effect is limited.Reasonable selection of layers with moderate lamina strength for fracturing stimulation,increasing the pump rate during fracturing and setting the cluster spacing reasonably are beneficial to improve the effect of reservoir stimulation.
文摘The relationship between hydrogen concentration and crack tip stress, strain field, hydrogen diffusion and internal pressure respectively in the crack tip process zone was investigated, and the length of the crack tip process zone of hydrogen-induced cracking （HIC） was determined. Based on the mechanism of fracture of micro-crack nucleation, a dislocation model was presented for the fracture criteria of HIC. The influence factors of pipeline tube fracture ductile KISCC in the presence of hydrogen was analyzed, and the critical pressure bearing capability of a pipeline with hydrogeninduced cracking and the critical J-integrity （JISCC） were calculated, which is very important for pipeline safety.
基金supported by National Natural Science Foundation of China(No.51674279)China Postdoctoral Science Foundation(No.2016M602227)a grant from National Science and Technology Major Project(No.2017ZX05049-006)
文摘The paper presents a novel hydraulic fracturing model for the characterization and simulation of the complex fracture network in shale gas reservoirs. We go beyond the existing method that uses planar or orthogonal conjugate fractures for representing the ＇＇complexity＇＇ of the network. Bifurcation of fractures is performed utilizing the Lindenmayer system based on fractal geometry to describe the fracture propagation pattern, density and network connectivity. Four controlling parameters are proposed to describe the details of complex fractures and stimulated reservoir volume（SRV）. The results show that due to the multilevel feature of fractal fractures, the model could provide a simple method for contributing reservoir volume calibration. The primary-and second-stage fracture networks across the overall SRV are the main contributions to the production, while the induced fracture network just contributes another 20% in the late producing period. We also conduct simulation with respect to different refracturing cases and find that increasing the complexity of the fracture network provides better performance than only enhancing the fracture conductivity.
文摘Accurate fluid flow simulation in geologically complex reservoirs is of particular importance in construction of reservoir simulators.General approaches in naturally fractured reservoir simulation involve use of unstructured grids or a structured grid coupled with locally unstructured grids and discrete fracture models.These methods suffer from drawbacks such as lack of flexibility and of ease of updating.In this study,I combined fracture modeling by elastic gridding which improves flexibility,especially in complex reservoirs.The proposed model revises conventional modeling fractures by hard rigid planes that do not change through production.This is a dubious assumption,especially in reservoirs with a high production rate in the beginning.The proposed elastic fracture modeling considers changes in fracture properties,shape and aperture through the simulation.This strategy is only reliable for naturally fractured reservoirs with high fracture permeability and less permeable matrix and parallel fractures with less cross-connections.Comparison of elastic fracture modeling results with conventional modeling showed that these assumptions will cause production pressure to enlarge fracture apertures and change fracture shapes,which consequently results in lower production compared with what was previously assumed.It is concluded that an elastic gridded model could better simulate reservoir performance.
基金This work was financially supported by the National Key R&D Program of China(No.2021YFC2900500)the National Natural Science Foundation of China(Nos.52074020 and 42202306)+2 种基金the Open Fund of State Key Laboratory of Water Resource Protection and Utilization in Coal Mining(No.WPUKFJJ2019-06)the Interdisciplinary Research Project for Young Teachers of USTB(Fundamental Research Funds for the Central Universities)(No.FRF-IDRY-21001)the Natural Science Foundation of Jiangsu Province,China(No.BK20200993).
文摘Structure plane is one of the important factors affecting the stability and failure mode of rock mass engineering.Rock mass structure characterization is the basic work of rock mechanics research and the important content of numerical simulation.A new 3-dimensional rough discrete fracture network(RDFN3D)model and its modeling method based on the Weierstrass-Mandelbrot(W-M)function were presented in this paper.The RDFN3D model,which improves and unifies the modelling methods for the complex structural planes,has been realized.The influence of fractal dimension,amplitude,and surface precision on the modeling parameters of RDFN3D was discussed.The reasonable W-M parameters suitable for the roughness coefficient of JRC were proposed,and the relationship between the mathematical model and the joint characterization was established.The RDFN3D together with the smooth 3-dimensional discrete fracture network(DFN3D)models were successfully exported to the drawing exchange format,which will provide a wide application in numerous numerical simulation codes including both the continuous and discontinuous methods.The numerical models were discussed using the COMSOL Multiphysics code and the 3-dimensional particle flow code,respectively.The reliability of the RDFN3D model was preliminarily discussed and analyzed.The roughness and spatial connectivity of the fracture networks have a dominant effect on the fluid flow patterns.The research results can provide a new geological model and analysis model for numerical simulation and engineering analysis of jointed rock mass.
基金supported by the National Natural Science Foundation of China(Grant Nos.51722406,61573018 and 51874335)the Shandong Provincial Natural Science Foundation(Grant JQ201808)+1 种基金the Fundamental Research Funds for the Central Universities(Grant 18CX02097A)the National Science and Technology Major Project of China(Grant 2016ZX05025001-006)
文摘In practical development of unconventional reservoirs,fracture networks are a highly conductive transport media for subsurface fluid flow.Therefore,it is crucial to clearly determine the fracture properties used in production forecast.However,it is different to calibrate the properties of fracture networks because it is an inverse problem with multi-patterns and highcomplexity of fracture distribution and inherent defect of multiplicity of solution.In this paper,in order to solve the problem,the complex fracture model is divided into two sub-systems,namely"Pattern A"and"Pattern B."In addition,the generation method is grouped into two categories.Firstly,we construct each sub-system based on the probability density function of the fracture properties.Secondly,we recombine the sub-systems into an integral complex fracture system.Based on the generation mechanism,the estimation of the complex fracture from dynamic performance and observation data can be solved as an inverse problem.In this study,the Bayesian formulation is used to quantify the uncertainty of fracture properties.To minimize observation data misfit immediately as it occurs,we optimize the updated properties by a simultaneous perturbation stochastic algorithm which requires only two measurements of the loss function.In numerical experiments,we firstly visualize that small-scale fractures significantly contribute to the flow simulation.Then,we demonstrate the suitability and effectiveness of the Bayesian formulation for calibrating the complex fracture model in the following simulation.
基金the National Science Fund for Distinguished Young Scholars(No.52225403)the Natural Science Foundation of Shanxi Province(No.202303021212073)the National Natural Science Foundation of China(No.52104210)。
文摘Understanding the impact of mining disturbances and creep deformation on the macroscopic deformation and the microscopic pore and fracture structures(MPFS)of coal is paramount for ensuring the secure extraction of coal resources.This study conducts cyclic loading-unloading and creep experiments on coal using a low-field nuclear magnetic resonance(NMR)experimental apparatus which is equipped with mechanical loading units,enabling real-time monitoring the T2spectrum.The experiments indicated that cyclic loading-unloading stress paths initiate internal damage within coal samples.Under identical creep stress conditions,coal samples with more initial damages had more substantial instantaneous deformation and creep deformation during the creep process.After undergoing nearly 35 h of staged creep,the total strains for coal samples CC01,CC02,and CC03 reach 2.160%,2.261%,and 2.282%,respectively.In the creep stage,the peak area ratio of seepage pores and microfractures(SPM)gradually diminishes.A higher degree of initial damage leads to a more pronounced compaction trend in the SPM of coal samples.Considering the porosity evolution of SPM during the creep process,this study proposes a novel fractional derivative model for the porosity evolution of SPM.The efficacy of the proposed model in predicting porosity evolution of SPM is substantiated through experimental validation.Furthermore,an analysis of the impact mechanisms on key parameters in the model was carried out.
文摘A new numerical model is presented to simulate fracture initiation and propagation in geological structures. This model is based on the recent amalgamation of established failure and fracture mechanics theory, which has been implemented to the finite difference FLAC code as a constitutive FISH userdefined-model. Validation of the model has been studied on the basis of comparing the transitional failure modes in rock. It is shown that the model is capable of accurately simulating fracture distributions over entire brittle to ductile rock phases. The application of the model during longwall retreat simulation highlighted several caving characteristics relevant to varying geological condition. The distribution and behaviour of modelled fractures were both realistic and shown to provide an enhanced post failure analysis to geological structures in FLAC. Moreover, the model introduces new potential insight towards the failure analysis of more complicated problems. This is best suited towards improving safety and efficiency in mines through the prediction of various key fractures and caving characteristics of geological structures.
基金supported by the National Natural Science Foundation of China(No.11702147)。
文摘Investigation of paper cutting process is vital for the design of cutting tools,but the fracture mechanism of paper cutting is still unclear.Here,we focus on the cutting process of paper,including the key parameters of cohesive zone model(CZM)for the orthotropic paper,to simulate the shear fracture process.Firstly,the material constants of the orthotropic paper are determined by longitudinal and transverse tensile test.Secondly,based on the tensile stressstrain curves,combined with damage theory and numerical simulations,the key parameters of the CZM for the orthotropic paper are obtained.Finally,a model III fracture is simulated to verify the accuracy of the model.Results show that the load-displacement curves obtained by the simulation is consistent with the test results.
基金Supported by the National Science and Technology Major Project(2017ZX05063-005)Science and Technology Development Project of PetroChina Research Institute of Petroleum Exploration and Development(YGJ2019-12-04)。
文摘The generation method of three-dimensional fractal discrete fracture network(FDFN)based on multiplicative cascade process was developed.The complex multi-scale fracture system in shale after fracturing was characterized by coupling the artificial fracture model and the natural fracture model.Based on an assisted history matching(AHM)using multiple-proxy-based Markov chain Monte Carlo algorithm(MCMC),an embedded discrete fracture modeling(EDFM)incorporated with reservoir simulator was used to predict productivity of shale gas well.When using the natural fracture generation method,the distribution of natural fracture network can be controlled by fractal parameters,and the natural fracture network generated coupling with artificial fractures can characterize the complex system of different-scale fractures in shale after fracturing.The EDFM,with fewer grids and less computation time consumption,can characterize the attributes of natural fractures and artificial fractures flexibly,and simulate the details of mass transfer between matrix cells and fractures while reducing computation significantly.The combination of AMH and EDFM can lower the uncertainty of reservoir and fracture parameters,and realize effective inversion of key reservoir and fracture parameters and the productivity forecast of shale gas wells.Application demonstrates the results from the proposed productivity prediction model integrating FDFN,EDFM and AHM have high credibility.
