Numerical manifold method for thermo-mechanical coupling simulation of fractured rock mass

As a calculation method based on the Galerkin variation,the numerical manifold method(NMM)adopts a double covering system,which can easily deal with discontinuous deformation problems and has a high calculation accuracy.Aiming at the thermo-mechanical(TM)coupling problem of fractured rock masses,this study uses the NMM to simulate the processes of crack initiation and propagation in a rock mass under the influence of temperature field,deduces related system equations,and proposes a penalty function method to deal with boundary conditions.Numerical examples are employed to confirm the effectiveness and high accuracy of this method.By the thermal stress analysis of a thick-walled cylinder(TWC),the simulation of cracking in the TWC under heating and cooling conditions,and the simulation of thermal cracking of the SwedishÄspöPillar Stability Experiment(APSE)rock column,the thermal stress,and TM coupling are obtained.The numerical simulation results are in good agreement with the test data and other numerical results,thus verifying the effectiveness of the NMM in dealing with thermal stress and crack propagation problems of fractured rock masses.

Effect mechanism of seepage force on the hydraulic fracture propagation

The flow of fluid through the porous matrix of a reservoir rock applies a seepage force to the solid rock matrix.Although the seepage force exerted by fluid flow through the porous matrix of a reservoir rock has a notable influence on rock deformation and failure,its effect on hydraulic fracture(HF)propagation remains ambiguous.Therefore,in this study,we improved a traditional fluid–solid coupling method by incorporating the role of seepage force during the fracturing fluid seepage,using the discrete element method.First,we validated the simulation results of the improved method by comparing them with an analytical solution of the seepage force and published experimental results.Next,we conducted numerical simulations in both homogeneous and heterogeneous sandstone formations to investigate the influence of seepage force on HF propagation.Our results indicate that fluid viscosity has a greater impact on the magnitude and extent of seepage force compared to injection rate,and that lower viscosity and injection rate correspond to shorter hydraulic fracture lengths.Furthermore,seepage force influences the direction of HF propagation,causing HFs to deflect towards the side of the reservoir with weaker cementation and higher permeability.

Study on the Impact of Massive Refracturing on the Fracture Network in Tight Oil Reservoir Horizontal Wells

Class III tight oil reservoirs have low porosity and permeability,which are often responsible for low production rates and limited recovery.Extensive repeated fracturing is a well-known technique to fix some of these issues.With such methods,existing fractures are refractured,and/or new fractures are created to facilitate communication with natural fractures.This study explored how different refracturing methods affect horizontal well fracture networks,with a special focus on morphology and related fluid flow changes.In particular,the study relied on the unconventional fracture model(UFM).The evolution of fracture morphology and flow field after the initial fracturing were analyzed accordingly.The simulation results indicated that increased formation energy and reduced reservoir stress differences can promote fracture expansion.It was shown that the length of the fracture network,the width of the fracture network,and the complexity of the fracture can be improved,the oil drainage area can be increased,the distance of oil and gas seepage can be reduced,and the production of a single well can be significantly increased.

Effect of fracture fluid flowback on shale microfractures using CT scanning

The field data of shale fracturing demonstrate that the flowback performance of fracturing fluid is different from that of conventional reservoirs,where the flowback rate of shale fracturing fluid is lower than that of conventional reservoirs.At the early stage of flowback,there is no single-phase flow of the liquid phase in shale,but rather a gas-water two-phase flow,such that the single-phase flow model for tight oil and gas reservoirs is not applicable.In this study,pores and microfractures are extracted based on the experimental results of computed tomography(CT)scanning,and a spatial model of microfractures is established.Then,the influence of rough microfracture surfaces on the flow is corrected using the modified cubic law,which was modified by introducing the average deviation of the microfracture height as a roughness factor to consider the influence of microfracture surface roughness.The flow in the fracture network is simulated using the modified cubic law and the lattice Boltzmann method(LBM).The results obtained demonstrate that most of the fracturing fluid is retained in the shale microfractures,which explains the low fracturing fluid flowback rate in shale hydraulic fracturing.

Numerical investigation of geostress influence on the grouting reinforcement effectiveness of tunnel surrounding rock mass in fault fracture zones

Grouting is a widely used approach to reinforce broken surrounding rock mass during the construction of underground tunnels in fault fracture zones,and its reinforcement effectiveness is highly affected by geostress.In this study,a numerical manifold method(NMM)based simulator has been developed to examine the impact of geostress conditions on grouting reinforcement during tunnel excavation.To develop this simulator,a detection technique for identifying slurry migration channels and an improved fluid-solid coupling(FeS)framework,which considers the influence of fracture properties and geostress states,is developed and incorporated into a zero-thickness cohesive element(ZE)based NMM(Co-NMM)for simulating tunnel excavation.Additionally,to simulate coagulation of injected slurry,a bonding repair algorithm is further proposed based on the ZE model.To verify the accuracy of the proposed simulator,a series of simulations about slurry migration in single fractures and fracture networks are numerically reproduced,and the results align well with analytical and laboratory test results.Furthermore,these numerical results show that neglecting the influence of geostress condition can lead to a serious over-estimation of slurry migration range and reinforcement effectiveness.After validations,a series of simulations about tunnel grouting reinforcement and tunnel excavation in fault fracture zones with varying fracture densities under different geostress conditions are conducted.Based on these simula-tions,the influence of geostress conditions and the optimization of grouting schemes are discussed.

Efficient production optimization for naturally fractured reservoir using EDFM

Naturally fractured reservoirs make important contributions to global oil and gas reserves and production.The modeling and simulation of naturally fractured reservoirs are different from conventional reservoirs as the existence of natural fractures.To address the development optimization problem of naturally fractured reservoirs,we propose an optimization workflow by coupling the optimization methods with the embedded discrete fracture model(EDFM).Firstly,the effective and superior performance of the workflow is verified based on the conceptual model.The stochastic simplex approximate gradient(StoSAG)algorithm,the ensemble optimization(EnOpt)algorithm,and the particle swarm optimization(PSO)algorithm are implemented for the production optimization of naturally fractured reservoirs based on the improved versions of the Egg model and the PUNQ-S3 model.The results of the two cases demonstrate the effectiveness of this optimization workflow by finding the optimal well controls which yield the maximum net present value(NPV).Compared to the initial well control guess,the final NPV obtained from the production optimization of fractured reservoirs based on all three optimization algorithms is significantly enhanced.Compared with the optimization results of the PSO algorithm,StoSAG and EnOpt have significant advantages in terms of final NPV and computational efficiency.The results also show that fractures have a significant impact on reservoir production.The economic efficiency of fractured reservoir development can be significantly improved by the optimization workflow.

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.

A new approach for flow simulation in complex hydraulic fracture morphology and its application:Fracture connection element method

Efficient flow simulation and optimization methods of hydraulic fracture morphology in unconventional reservoirs are effective ways to enhance oil/gas recovery.Based on the connection element method(CEM)and distribution of stimulated reservoir volume,the complex hydraulic fracture morphology was accurately described using heterogeneous node connection system.Then a new fracture connection element method(FCEM)for fluid flow in stimulated unconventional reservoirs with complex hydraulic fracture morphology was proposed.In the proposed FCEM,the arrangement of dense nodes in the stimulated area and sparse nodes in the unstimulated area ensures the calculation accuracy and efficiency.The key parameter,transmissibility,was also modified according to the strong heterogeneity of stimulated reservoirs.The finite difference and semi-analytical tracking were used to accurately solve the pressure and saturation distribution between nodes.The FCEM is validated by comparing with traditional numerical simulation method,and the results show that the bottom hole pressure simulated by the FCEM is consistent with the results from traditional numerical simulation method,and the matching rate is larger than 95%.The proposed FCEM was also used in the optimization of fracturing parameters by coupling the hydraulic fracture propagation method and intelligent optimization algorithm.The integrated intelligent optimization approach for multi-parameters,such as perforation number,perforation location,and displacement in hydraulic fracturing is proposed.The proposed approach was applied in a shale gas reservoir,and the result shows that the optimized perforation location and morphology distribution are related to the distribution of porosity/permeability.When the perforation location and displacement are optimized with the same fracture number,NPV increases by 70.58%,which greatly improves the economic benefits of unconventional reservoirs.This work provides a new way for flow simulation and optimization of hydraulic fracture morphology of multi-fractured horizontal wells in unconventional reservoirs.

A method for predicting the water-flowing fractured zone height based on an improved key stratum theory

In the process of using the original key stratum theory to predict the height of a water-flowing fractured zone(WFZ),the influence of rock strata outside the calculation range on the rock strata within the calculation range as well as the fact that the shape of the overburden deformation area will change with the excavation length are ignored.In this paper,an improved key stratum theory(IKS theory)was proposed by fixing these two shortcomings.Then,a WFZ height prediction method based on IKS theory was established and applied.First,the range of overburden involved in the analysis was determined according to the tensile stress distribution range above the goaf.Second,the key stratum in the overburden involved in the analysis was identified through IKS theory.Finally,the tendency of the WFZ to develop upward was determined by judging whether or not the identified key stratum will break.The proposed method was applied and verified in a mining case study,and the reasons for the differences in the development patterns between the WFZs in coalfields in Northwest and East China were also fully explained by this method.

An efficient 3D cell-based discrete fracture-matrix flow model for digitally captured fracture networks

Complex hydraulic fracture networks are critical for enhancing permeability in unconventional reservoirs and mining indus-tries.However,accurately simulating the fluid flow in realistic fracture networks(compared to the statistical fracture net-works)is still challenging due to the fracture complexity and computational burden.This work proposes a simple yet efficient numerical framework for the flow simulation in fractured porous media obtained by 3D high-resolution images,aiming at both computational accuracy and efficiency.The fractured rock with complex fracture geometries is numerically constructed with a cell-based discrete fracture-matrix model(DFM)having implicit fracture apertures.The flow in the complex fractured porous media(including matrix flow,fracture flow,as well as exchange flow)is simulated with a pipe-based cell-centered finite volume method.The performance of this model is validated against analytical/numerical solutions.Then a lab-scale true triaxial hydraulically fractured shale sample is reconstructed,and the fluid flow in this realistic fracture network is simu-lated.Results suggest that the proposed method achieves a good balance between computational efficiency and accuracy.The complex fracture networks control the fluid flow process,and the opened natural fractures behave as primary fuid pathways.Heterogeneous and anisotropic features of fluid flow are well captured with the present model.

A modified single edge V-notched beam method for evaluating surface fracture toughness of thermal barrier coatings

The surface fracture toughness is an important mechanical parameter for studying the failure behavior of air plasma sprayed(APS)thermal barrier coatings(TBCs).As APS TBCs are typical multilayer porous ceramic materials,the direct applications of the traditional single edge notched beam(SENB)method that ignores those typical structural characters may cause errors.To measure the surface fracture toughness more accurately,the effects of multilayer and porous characters on the fracture toughness of APS TBCs should be considered.In this paper,a modified single edge V-notched beam(MSEVNB)method with typical structural characters is developed.According to the finite element analysis(FEA),the geometry factor of the multilayer structure is recalculated.Owing to the narrower V-notches,a more accurate critical fracture stress is obtained.Based on the Griffith energy balance,the reduction of the crack surface caused by micro-defects is corrected.The MSEVNB method can measure the surface fracture toughness more accurately than the SENB method.

Blount Method in the Management of Supracondylar Fractures of the Humerus in Children

Supracondylar fractures of the humerus are the most common traumatic elbow injury in children. Several techniques are described in their support. Blount’s method is one of the widely used orthopedic means. The objective of our study was to evaluate the results of supracondylar fractures of the humerus treated by this method in our department. This was a descriptive, longitudinal study with the prospective collection over a period of 12 months from May 2018 to April 2019. It focused on children aged 0 to 15 years treated in the Orthopedic Traumatology Department of the CHU of Kati for a supracondylar fracture of the humerus in extension. A minimum follow-up period of 6 months was necessary for inclusion. 12 patients were collected. The sociodemographic, anatomopathologic, therapeutic and evolutionary aspects were studied. The results were evaluated according to Flynn’s criteria. The male sex was predominant with a ratio of 2. The average age was 5.5 years. The average admission time was 12.5 hours. Domestic accidents were the main etiology (66.66%) and Type II fracture was the most frequent with 58.33%. The duration of immobilization was 42 days for 11 patients. Complications were 1 case of secondary displacement and 1 case of 12&#730 flexion deficit. The average Baumann angle was 72&#730. The result was satisfactory for all our patients. Supracondylar fractures of the humerus are common in children, especially boys. The method of Blount when it is well carried out by respecting the indications gives good results.

Fracture propagation and evolution law of indirect fracturing in the roof of broken soft coal seams

Indirect fracturing in the roof of broken soft coal seams has been demonstrated to be a feasible technology.In this work,the No.5 coal seam in the Hancheng block was taken as the research object.Based on the findings of true triaxial hydraulic fracturing experiments and field pilot under this technology and the cohesive element method,a 3D numerical model of indirect fracturing in the roof of broken soft coal seams was established,the fracture morphology propagation and evolution law under different conditions was investigated,and analysis of main controlling factors of fracture parameters was conducted with the combination weight method,which was based on grey incidence,analytic hierarchy process and entropy weight method.The results show that“士”-shaped fractures,T-shaped fractures,cross fractures,H-shaped fractures,and“干”-shaped fractures dominated by horizontal fractures were formed.Different parameter combinations can form different fracture morphologies.When the coal seam permeability is lower and the minimum horizontal principal stress difference between layers and fracturing fluid injection rate are both larger,it tends to form“士”-shaped fractures.When the coal seam permeability and minimum horizontal principal stress between layers and perforation position are moderate,cross fractures are easily generated.Different fracture parameters have different main controlling factors.Engineering factors of perforation location,fracturing fluid injection rate and viscosity are the dominant factors of hydraulic fracture shape parameters.This study can provide a reference for the design of indirect fracturing in the roof of broken soft coal seams.

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.

Modeling the Interaction between Vacancies and Grain Boundaries during Ductile Fracture

The experimental results in previous studies have indicated that during the ductile fracture of pure metals,vacancies aggregate and form voids at grain boundaries.However,the physical mechanism underlying this phenomenon remains not fully understood.This study derives the equilibrium distribution of vacancies analytically by following thermodynamics and the micromechanics of crystal defects.This derivation suggests that vacancies cluster in regions under hydrostatic compression to minimize the elastic strain energy.Subsequently,a finite element model is developed for examining more general scenarios of interaction between vacancies and grain boundaries.This model is first verified and validated through comparison with some available analytical solutions,demonstrating consistency between finite element simulation results and analytical solutions within a specified numerical accuracy.A systematic numerical study is then conducted to investigate the mechanism that might govern the micromechanical interaction between grain boundaries and the profuse vacancies typically generated during plastic deformation.The simulation results indicate that the reduction in total elastic strain energy can indeed drive vacancies toward grain boundaries,potentially facilitating void nucleation in ductile fracture.

A phase-field model for simulating the propagation behavior of mixed-mode cracks during the hydraulic fracturing process in fractured reservoirs

A novel phase-field model for the propagation of mixed-mode hydraulic fractures,characterized by the formation of mixed-mode fractures due to the interactions between fluids and solids,is proposed.In this model,the driving force for the phase field consists of both tensile and shear components,with the fluid contribution primarily manifesting in the tension driving force.The displacement and pressure are solved simultaneously by an implicit method.The numerical solution's iterative format is established by the finite element discretization and Newton-Raphson(NR)iterative methods.The correctness of the model is verified through the uniaxial compression physical experiments on fluid-pressurized rocks,and the limitations of the hydraulic fracture expansion phase-field model,which only considers mode I fractures,are revealed.In addition,the influence of matrix mode II fracture toughness value,natural fracture mode II toughness value,and fracturing fluid injection rate on the hydraulic fracture propagation in porous media with natural fractures is studied.

The Boundary Element Method for Ordinary State-Based Peridynamics

The peridynamics(PD),as a promising nonlocal continuum mechanics theory,shines in solving discontinuous problems.Up to now,various numerical methods,such as the peridynamic mesh-free particlemethod(PD-MPM),peridynamic finite element method(PD-FEM),and peridynamic boundary element method(PD-BEM),have been proposed.PD-BEM,in particular,outperforms other methods by eliminating spurious boundary softening,efficiently handling infinite problems,and ensuring high computational accuracy.However,the existing PD-BEM is constructed exclusively for bond-based peridynamics(BBPD)with fixed Poisson’s ratio,limiting its applicability to crack propagation problems and scenarios involving infinite or semi-infinite problems.In this paper,we address these limitations by introducing the boundary element method(BEM)for ordinary state-based peridynamics(OSPD-BEM).Additionally,we present a crack propagationmodel embeddedwithin the framework ofOSPD-BEM to simulate crack propagations.To validate the effectiveness of OSPD-BEM,we conduct four numerical examples:deformation under uniaxial loading,crack initiation in a double-notched specimen,wedge-splitting test,and threepoint bending test.The results demonstrate the accuracy and efficiency of OSPD-BEM,highlighting its capability to successfully eliminate spurious boundary softening phenomena under varying Poisson’s ratios.Moreover,OSPDBEMsignificantly reduces computational time and exhibits greater consistencywith experimental results compared to PD-MPM.

Micromechanical modeling of asphalt concrete fracture using a user-defined three-dimensional discrete element method

A user-defined micromechanical model was developed to investigate the fracture mechanism of asphalt concrete(AC)using the discrete element method(DEM).A three-dimensional(3D)AC beam was built using the "Fish" language provided by PFC3D and was employed to simulate the three-point bending beam test at two temperature levels:-10°C and 15°C.The AC beam was modeled with the consideration of the microstructural features of asphalt mixtures.Uniaxial complex modulus test and indirect tensile strength test were conducted to obtain material input parameters for numerical modeling.The 3D predictions were validated using laboratory experimental measurements of AC beams prepared by the same mixture design.Effects of mastic stiffness,cohesive and adhesive strength on AC fracture behavior were investigated using the DEM model.The results show that the 3D DEM fracture model can accurately predict the fracture patterns of asphalt concrete.The ratio of stress at interfaces to the stress in mastics increases as the mastic stiffness decreases;however,the increase in the cohesive strength or adhesive strength shows no significant influence on the tensile strength.

Mode Ⅲ fracture analysis by Trefftz boundary element method

This paper presents a hybrid Trefftz （HT） boundary element method （BEM） by using two indirect techniques for mode III fracture problems. Two Trefftz complete functions of Laplace equation for normal elements and a special purpose Trefftz function for crack elements are proposed in deriving the Galerkin and the collocation techniques of HT BEM. Then two auxiliary functions are introduced to improve the accuracy of the displacement field near the crack tips, and stress intensity factor （SIF） is evaluated by local crack elements as well. Furthermore, numerical examples are given, including comparisons of the present results with the analytical solution and the other numerical methods, to demonstrate the efficiency for different boundary conditions and to illustrate the convergence influenced by several parameters. It shows that HT BEM by usingthe Galerkin and the collocation techniques is effective for mode III fracture problems.

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.