Algorithmic approach to discrete fracture network flow modeling in consideration of realistic connections in large-scale fracture networks

Analyzing rock mass seepage using the discrete fracture network(DFN)flow model poses challenges when dealing with complex fracture networks.This paper presents a novel DFN flow model that incorporates the actual connections of large-scale fractures.Notably,this model efficiently manages over 20,000 fractures without necessitating adjustments to the DFN geometry.All geometric analyses,such as identifying connected fractures,dividing the two-dimensional domain into closed loops,triangulating arbitrary loops,and refining triangular elements,are fully automated.The analysis processes are comprehensively introduced,and core algorithms,along with their pseudo-codes,are outlined and explained to assist readers in their programming endeavors.The accuracy of geometric analyses is validated through topological graphs representing the connection relationships between fractures.In practical application,the proposed model is employed to assess the water-sealing effectiveness of an underground storage cavern project.The analysis results indicate that the existing design scheme can effectively prevent the stored oil from leaking in the presence of both dense and sparse fractures.Furthermore,following extensive modification and optimization,the scale and precision of model computation suggest that the proposed model and developed codes can meet the requirements of engineering applications.

Fractured reservoir modeling by discrete fracture network and seismic modeling in the Tarim Basin,China

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.

Development of an improved three-dimensional rough discrete fracture network model:Method and application

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.

Borehole stability in naturally fractured rocks with drilling mud intrusion and associated fracture strength weakening:A coupled DFN-DEM approach

Borehole instability in naturally fractured rocks poses significant challenges to drilling.Drilling mud invades the surrounding formations through natural fractures under the difference between the wellbore pressure(P w)and pore pressure(P p)during drilling,which may cause wellbore instability.However,the weakening of fracture strength due to mud intrusion is not considered in most existing borehole stability analyses,which may yield significant errors and misleading predictions.In addition,only limited factors were analyzed,and the fracture distribution was oversimplified.In this paper,the impacts of mud intrusion and associated fracture strength weakening on borehole stability in fractured rocks under both isotropic and anisotropic stress states are investigated using a coupled DEM(distinct element method)and DFN(discrete fracture network)method.It provides estimates of the effect of fracture strength weakening,wellbore pressure,in situ stresses,and sealing efficiency on borehole stability.The results show that mud intrusion and weakening of fracture strength can damage the borehole.This is demonstrated by the large displacement around the borehole,shear displacement on natural fractures,and the generation of fracture at shear limit.Mud intrusion reduces the shear strength of the fracture surface and leads to shear failure,which explains that the increase in mud weight may worsen borehole stability during overbalanced drilling in fractured formations.A higher in situ stress anisotropy exerts a significant influence on the mechanism of shear failure distribution around the wellbore.Moreover,the effect of sealing natural fractures on maintaining borehole stability is verified in this study,and the increase in sealing efficiency reduces the radial invasion distance of drilling mud.This study provides a directly quantitative prediction method of borehole instability in naturally fractured formations,which can consider the discrete fracture network,mud intrusion,and associated weakening of fracture strength.The information provided by the numerical approach(e.g.displacement around the borehole,shear displacement on fracture,and fracture at shear limit)is helpful for managing wellbore stability and designing wellbore-strengthening operations.

Characterizing the influence of stress-induced microcracks on the laboratory strength and fracture development in brittle rocks using a finite-discrete element method-micro discrete fracture network FDEM-μDFN approach

Heterogeneity is an inherent component of rock and may be present in different forms including mineralheterogeneity, geometrical heterogeneity, weak grain boundaries and micro-defects. Microcracks areusually observed in crystalline rocks in two forms： natural and stress-induced; the amount of stressinducedmicrocracking increases with depth and in-situ stress. Laboratory results indicate that thephysical properties of rocks such as strength, deformability, P-wave velocity and permeability areinfluenced by increase in microcrack intensity. In this study, the finite-discrete element method （FDEM）is used to model microcrack heterogeneity by introducing into a model sample sets of microcracks usingthe proposed micro discrete fracture network （mDFN） approach. The characteristics of the microcracksrequired to create mDFN models are obtained through image analyses of thin sections of Lac du Bonnetgranite adopted from published literature. A suite of two-dimensional laboratory tests including uniaxial,triaxial compression and Brazilian tests is simulated and the results are compared with laboratory data.The FDEM-mDFN models indicate that micro-heterogeneity has a profound influence on both the mechanicalbehavior and resultant fracture pattern. An increase in the microcrack intensity leads to areduction in the strength of the sample and changes the character of the rock strength envelope. Spallingand axial splitting dominate the failure mode at low confinement while shear failure is the dominantfailure mode at high confinement. Numerical results from simulated compression tests show thatmicrocracking reduces the cohesive component of strength alone, and the frictional strength componentremains unaffected. Results from simulated Brazilian tests show that the tensile strength is influenced bythe presence of microcracks, with a reduction in tensile strength as microcrack intensity increases. Theimportance of microcrack heterogeneity in reproducing a bi-linear or S-shape failure envelope and itseffects on the mechanisms leading to spalling damage near an underground opening are also discussed.

A discrete model for prediction of radon flux from fractured rocks

Prediction of radon flux from the fractured zone of a propagating cave mine is basically associated with uncertainty and complexity. For instance, there is restricted access to these zones for field measure- ments, and it is quite difficult to replicate the complex nature of both natural and induced fractures in these zones in laboratory studies. Hence, a technique for predicting radon flux from a fractured rock using a discrete fracture network （DFN） model is developed to address these difficulties. This model quantifies the contribution of fractures to the total radon flux, and estimates the fracture density from a measured radon flux considering the effects of advection, diffusion, as well as radon generation and decay. Radon generation and decay are classified as reaction processes. Therefore, the equation solved is termed as the advection-diffusion-reaction equation （ADRE）. Peclet number （Pe）, a conventional dimensionless parameter that indicates the ratio of mass transport by advection to diffusion, is used to classify the transport regimes. The results show that the proposed model effectively predicts radon flux from a fractured rock. An increase in fracture density for a rock sample with uniformly distributed radon generation rate can elevate radon flux significantly compared with another rock sample with an equivalent increase in radon generation rate. In addition to Pe, two other independent dimensionless parameters （derived for radon transport through fractures） significantly affect radon dimensionless flux. Findings provide insight into radon transport through fractured rocks and can be used to improve radon control measures for proactive mitigation.

Complex fracture propagation model and plugging timing optimization for temporary plugging fracturing in naturally fractured shale

In this paper,a viscoelasticity-plastic damage constitutive equation for naturally fractured shale is deduced,coupling nonlinear tensile-shear mixed fracture mode.Dynamic perforation-erosion on fluid re-distribution among multi-clusters are considered as well.DFN-FEM(discrete fracture network combined with finite element method)was developed to simulate the multi-cluster complex fractures propagation within temporary plugging fracturing(TPF).Numerical results are matched with field injection and micro-seismic monitoring data.Based on geomechanical characteristics of Weiyuan deep shale gas reservoir in Sichuan Basin,SW China,a multi-cluster complex fractures propagation model is built for TPF.To study complex fractures propagation and the permeability-enhanced region evolution,intersecting and competition mechanisms between the fractures before and after TPF treatment are revealed.Simulation results show that:fracture from middle cluster is restricted by the fractures from side-clusters,and side-clusters plugging is benefit for multi fractures propagation in uniformity;optimized TPF timing should be delayed within a higher density or strike of natural fractures;Within a reservoir-featured natural fractures distribution,optimized TPF timing for most clustered method is 2/3 of total fluid injection time as the optimal plugging time under different clustering modes.

Multi-dimensional size effects and representative elements for nonpersistent fractured rock masses: A perspective of geometric parameter distribution

This study takes a fractured rock mass in the Datengxia Hydropower Station,China as an example to analyze the size effects and determine the representative elementary sizes.A novel method considering geometric parameter distributions is proposed in this work.The proposed method can quickly and simply determine the size effects and representative elementary sizes.Specifically,geometric parameter distributions,including fracture frequency,size and orientation,are generated on the basis of the Bernoulli trial and Monte Carlo simulation.The distributions are assessed using the coefficient of variation(CV),and the acceptable variations for CV(5%,10%and 20%)are used to determine representative elementary sizes.Generally,the representative element of rock masses is the representative elementary volume(REV).The present study extends the representative element to other dimensions,i.e.representative elementary length(REL)and representative elementary area(REA)for one and two dimensions,respectively.REL and REA are useful in studying the size effects of one-(1D)and twodimensional(2D)characteristics of rock masses.The relationships among multi-dimensional representative elementary sizes are established.The representative elementary sizes reduce with the increase in the dimensions,and REA and REV can be deduced by REL.Therefore,the proposed method can quickly and simply determine REL and further estimate REA and REV,which considerably improves the efficiency of rock mass analysis.

Multi-domain equivalent method for prediction of elastic modulus of complex fractured rock mass

Reliable estimation of deformation and failure behaviors of fractured rock mass is important for practical engineering design.This study proposes a multi-domain equivalent method for fracture network to estimate the deformation properties of complex fractured rock mass.It comprehends both the advantages of the discrete fracture network model and the equivalent continuum model to capture the features of discontinuities explicitly while reducing computational intensity.The complex fracture network is stochastically split into a number of subfracture networks according to the domain,length or angle.An analytical solution is derived to infer theoretically the relationship between the elastic moduli of the original complex fractured rock mass and the split subfractured rock masses by introducing a correction term based on the deformation superposition principle.Numerical simulations are conducted to determine the elastic moduli of split subfractured rock masses using universal distinct element code(UDEC),while the elastic modulus of the original model is estimated based on the currently proposed analytical relationship.The results show that the estimation accuracy with the current domainbased splitting model is far superior compared to those with the other two splitting models.Thus,the estimation method of elastic modulus of complex fractured rock mass based on domain splitting mode of fracture network is identified as the multi-domain equivalent method proposed in this paper.The reliability of this method is evaluated,and its high computational efficiency is demonstrated through exemplification with regard to different geometric configurations for stochastically artificial discrete fracture network.The proposed multi-domain equivalent method constructs the theoretical framework except for the regression analysis hypothesis compared to the density-reduced model equivalent method.

Development of a DFN-based probabilistic block theory approach for bench face angle design in open pit mining

In open pit mining,uncontrolled block instabilities have serious social,economic and regulatory consequences,such as casualties,disruption of operation and increased regulation difficulties.For this reason,bench face angle,as one of the controlling parameters associated with block instabilities,should be carefully designed for sustainable mining.This study introduces a discrete fracture network(DFN)-based probabilistic block theory approach for the fast design of the bench face angle.A major advantage is the explicit incorporation of discontinuity size and spatial distribution in the procedure of key blocks testing.The proposed approach was applied to a granite mine in China.First,DFN models were generated from a multi-step modeling procedure to simulate the complex structural characteristics of pit slopes.Then,a modified key blocks searching method was applied to the slope faces modeled,and a cumulative probability of failure was obtained for each sector.Finally,a bench face angle was determined commensurate with an acceptable risk level of stability.The simulation results have shown that the number of hazardous traces exposed on the slope face can be significantly reduced when the suggested bench face angle is adopted,indicating an extremely low risk of uncontrolled block instabilities.

A three-dimensional feature extraction-based method for coal cleat characterization using X-ray μCT and its application to a Bowen Basin coal specimen

Cleats are the dominant micro-fracture network controlling the macro-mechanical behavior of coal.Improved understanding of the spatial characteristics of cleat networks is therefore important to the coal mining industry.Discrete fracture networks(DFNs)are increasingly used in engineering analyses to spatially model fractures at various scales.The reliability of coal DFNs largely depends on the confidence in the input cleat statistics.Estimates of these parameters can be made from image-based three-dimensional(3D)characterization of coal cleats using X-ray micro-computed tomography(m CT).One key step in this process,after cleat extraction,is the separation of individual cleats,without which the cleats are a connected network and statistics for different cleat sets cannot be measured.In this paper,a feature extraction-based image processing method is introduced to identify and separate distinct cleat groups from 3D X-ray m CT images.Kernels(filters)representing explicit cleat features of coal are built and cleat separation is successfully achieved by convolutional operations on 3D coal images.The new method is applied to a coal specimen with 80 mm in diameter and 100 mm in length acquired from an Anglo American Steelmaking Coal mine in the Bowen Basin,Queensland,Australia.It is demonstrated that the new method produces reliable cleat separation capable of defining individual cleats and preserving 3D topology after separation.Bedding-parallel fractures are also identified and separated,which has his-torically been challenging to delineate and rarely reported.A variety of cleat/fracture statistics is measured which not only can quantitatively characterize the cleat/fracture system but also can be used for DFN modeling.Finally,variability and heterogeneity with respect to the core axis are investigated.Significant heterogeneity is observed and suggests that the representative elementary volume(REV)of the cleat groups for engineering purposes may be a complex problem requiring careful consideration.

Evaluating the mechanical properties of anisotropic shale containing bedding and natural fractures with discrete element modeling

Natural fracture data from one of the Carboniferous shale masses in the eastern Qaidam Basin were used to establish a stochastic model of a discrete fracture network and to perform discrete element simulation research on the size efect and mechanical parameters of shale.Analytical solutions of fctitious joints in transversely isotropic media were derived,which made it possible for the proposed numerical model to simulate the bedding and natural fractures in shale masses.The results indicate that there are two main factors infuencing the representative elementary volume(REV)size of a shale mass.The frst and most decisive factor is the presence of natural fractures in the block itself.The second is the anisotropy ratio:the greater the anisotropy is,the larger the REV.The bedding angle has little infuence on the REV size,whereas it has a certain infuence on the mechanical parameters of the rock mass.When the bedding angle approaches the average orientation of the natural fractures,the mechanical parameters of the shale blocks decrease greatly.The REV representing the mechanical properties of the Carboniferous shale masses in the eastern Qaidam Basin were comprehensively identifed by considering the infuence of bedding and natural fractures.When the numerical model size is larger than the REV,the fractured rock mass discontinuities can be transformed into equivalent continuities,which provides a method for simulating shale with natural fractures and bedding to analyze the stability of a borehole wall in shale.

Numerical modeling of thermal breakthrough induced by geothermal production in fractured granite

It is well known that the complicated channeling of fluid flow and heat transfer is strongly related with the intricate natural fracture system.However,it is still challenging to set up the fracture network model which is strong heterogeneous.Compared with other methods(e.g.equivalent continuum model(ECM),discrete fracture model(DFM),and ECM-DFM),the fracture flow module in the COMSOL Multiphysics simulator is powerful in definition of fractures as the inner flow boundary existing in the porous media.Thus it is selected to simulate the fluid flow and heat transfer in the geothermal-developed fractured granite of Sanguliu area located at Liaodong Peninsula,Eastern China.The natural faults/fractures based on field investigation combined with the discrete fracture network(DFN)generated by the MATLAB are used to represent the two-dimensional geological model.Numerical results show that early thermal breakthrough occurs at the production well caused by quick flow of cold water along the highly connected fractures.Suitable hydraulic fracturing treatments with proper injection rates,locations,etc.can efficiently hinder the thermal breakthrough time in the natural fracture system.Large well spacing helps the long-term operation of geothermal production,but it is highly dependent on the geometrical morphology of the fracture network.The enhancement of reservoir properties at the near-well regions can also increase the geothermal production efficiency.The results in this study can provide references to achieve a sustainable geothermal exploitation in fractured granitic geothermal reservoirs or hot dry rocks at depth.

Effect of seismic waves on the hydro-mechanical properties of fractured rock masses

The transmission of seismic waves in a particular region may influence the hydraulic properties of a rock mass, including permeability, which is one of the most important. To determine the effect of a seismic wave on the hydraulic behavior of a fractured rock mass, systematic numerical modeling was conducted. A number of discrete fracture network（DFN） models with a size of 20 m × 20 m were used as geometrical bases, and a discrete element method（DEM） was employed as a numerical simulation tool. Three different boundary conditions without（Type Ⅰ） and with static（Type Ⅱ） and dynamic（Type Ⅲ） loading were performed on the models, and then their permeability was calculated. The results showed that permeability in Type Ⅲ models was respectively 62.7% and 44.2% higher than in Type I and Type Ⅱ models. This study indicates that seismic waves can affect deep earth, and, according to the results, seismic waves increase the permeability and change the flow rate patterns in a fractured rock mass.

Anisotropy of strength and deformability of fractured rocks

Anisotropy of the strength and deformation behaviors of fractured rock masses is a crucial issue for design and stability assessments of rock engineering structures, due mainly to the non-uniform and non- regular geometries of the fracture systems. However, no adequate efforts have been made to study this issue due to the current practical impossibility of laboratory tests with samples of large volumes con- taining many fractures, and the difficulty for controlling reliable initial and boundary conditions for large-scale in situ tests. Therefore, a reliable numerical predicting approach for evaluating anisotropy of fractured rock masses is needed. The objective of this study is to systematically investigate anisotropy of strength and deformability of fractured rocks, which has not been conducted in the past, using a nu- merical modeling method. A series of realistic two-dimensional （2D） discrete fracture network （DFN） models were established based on site investigation data, which were then loaded in different directions, using the code UDEC of discrete element method （DEM）, with changing confining pressures. Numerical results show that strength envelopes and elastic deformability parameters of tested numerical models are significantly anisotropic, and vary with changing axial loading and confining pressures. The results indicate that for design and safety assessments of rock engineering projects, the directional variations of strength and deformability of the fractured rock mass concerned must be treated properly with respect to the directions of in situ stresses. Traditional practice for simply positioning axial orientation of tunnels in association with principal stress directions only may not be adequate for safety requirements. Outstanding issues of the present study and su^zestions for future study are also oresented.

A smart productivity evaluation method for shale gas wells based on 3D fractal fracture network model

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.

A Review of the Dynamic Modeling Approaches for Characterizing Fluid Flow in Naturally Fractured Reservoirs

Fluid flow in fractured media has been studied for decades and received considerable attention in the oil and gas industry because of the high productivity of naturally fractured reservoirs.Due to formation complexity and reservoir heterogeneity,characterizing fluid flow with an appropriate reservoir model presents a challenging task that differs relatively from homogeneous conventional reservoirs in many aspects of view,including geological,petrophysical,production,and economics.In most fractured reservoirs,fracture networks create complex pathways that affect hydrocarbon flow,well performance,hence reservoir characterization.A better and comprehensive understanding of the available reservoir modeling approaches is much needed to accurately characterize fluid flow behavior in NFRs.Therefore,in this paper,a perspective review of the available modeling approaches was presented for fluid flow characterization in naturally fractured medium.Modeling methods were evaluated in terms of their description,application,advantages,and disadvantages.This study has also included the applications of these reservoir models in fluid flow characterizing studies and governing equations for fluid flow.Dual continuum models were proved to be better than single continuum models in the presence of large scale fractures.In comparison,discrete models were more appropriate for reservoirs that contain a smaller number of fractures.However,hybrid modeling was the best method to provide accurate and scalable fluid flow modeling.It is our understanding that this paper will bridge the gap between the fundamental understanding and application of NFRs modeling approaches and serve as a useful reference for engineers and researchers for present and future applications.

Estimation of fracture size and azimuth in the universal elliptical disc model based on trace information

The geometric characteristics of fractures within a rock mass can be inferred by the data sampling from boreholes or exposed surfaces.Recently,the universal elliptical disc(UED)model was developed to represent natural fractures,where the fracture is assumed to be an elliptical disc and the fracture orientation,rotation angle,length of the long axis and ratio of short-long axis lengths are considered as variables.This paper aims to estimate the fracture size-and azimuth-related parameters in the UED model based on the trace information from sampling windows.The stereological relationship between the trace length,size-and azimuth-related parameters of the UED model was established,and the formulae of the mean value and standard deviation of trace length were proposed.The proposed formulae were validated via the Monte Carlo simulations with less than 5%of error rate between the calculated and true values.With respect to the estimation of the size-and azimuth-related parameters using the trace length,an optimization method was developed based on the pre-assumed size and azimuth distribution forms.A hypothetical case study was designed to illustrate and verify the parameter estimation method,where three combinations of the sampling windows were used to estimate the parameters,and the results showed that the estimated values could agree well with the true values.Furthermore,a hypothetical three-dimensional(3D)elliptical fracture network was constructed,and the circular disc,non-UED and UED models were used to represent it.The simulated trace information from different models was compared,and the results clearly illustrated the superiority of the proposed UED model over the existing circular disc and non-UED models。

New approaches to quantify progressive damage and associated dynamic rock mass blockiness

In the past decade, numerical modelling has been increasingly used for simulating the mechanical behaviour of naturally fractured rock masses. In this paper, we introduce new algorithms for spatial and temporal analyses of newly generated fractures and blocks using an integrated discrete fracture network (DFN)-finite-discrete element method (FDEM) (DFN-FDEM) modelling approach. A fracture line calculator and analysis technique (i.e. discrete element method (DEM) fracture analysis, DEMFA) calculates the geometrical aspects of induced fractures using a dilation criterion. The resultant two-dimensional (2D) blocks are then identified and characterised using a graph structure. Block tracking trees allow track of newly generated blocks across timesteps and to analyse progressive breakage of these blocks into smaller blocks. Fracture statistics (number and total length of initial and induced fractures) are then related to the block forming processes to investigate damage evolution. The combination of various proposed methodologies together across various stages of modelling processes provides new insights to investigate the dependency of structure's resistance on the initial fracture configuration.

Numerical Modeling of Natural Fracture Distributions in Shale

The production efficiency of shale gas is affected by the interaction between hydraulic and natural fractures.This study presents a simulation of natural fractures in shale reservoirs,based on a discrete fracture network(DFN)method for hydraulic fracturing engineering.Fracture properties of the model are calculated from core fracture data,according to statistical mathematical analysis.The calculation results make full use of the quantitative information of core fracture orientation,density,opening and length,which constitute the direct and extensive data of mining engineering.The reliability and applicability of the model are analyzed with regard to model size and density,a calculation method for dominant size and density being proposed.Then,finite element analysis is applied to a hydraulic fracturing numerical simulation of a shale fractured reservoir in southeastern Chongqing.The hydraulic pressure distribution,fracture propagation,acoustic emission information and in situ stress changes during fracturing are analyzed.The results show the application of fracture statistics in fracture modeling and the influence of fracture distribution on hydraulic fracturing engineering.The present analysis may provide a reference for shale gas exploitation.