Creep constitutive model considering nonlinear creep degradation of fractured rock

Stability analysis of underground constructions requires a model study of rock masses’ long-term performance. Creep tests under different stress conditions was conducted on intact granite and granite samples fractured at 30° and 45° angles. The experimental results indicate that the steady creep strain rates of intact and fractured rock present an exponential increase trend with the increase of stress level. A nonlinear creep model is developed based on the experimental results, in which the initial damage caused by fracture together with the damage caused by constant load have been taken into consideration. The fitting analysis results indicated that the model proposed is more accurate at identifying the full creep regions in fractured granite, especially the accelerated stage of creep deformation. The least-square fit error of the proposed creep model is significantly lower than that of Nishihara model by almost an order of magnitude. An analysis of the effects of elastic modulus, viscosity coefficient, and damage factors on fractured rock strain rate and creep strain is conducted. If no consideration is given to the effects of the damage, the proposed nonlinear creep model can degenerate into to the classical Nishihara model.

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.

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.

A new elastoplastic model for bolt-grouted fractured rock

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.

Identifying flow regime in the aquifer of fractured rock system in Germi Chai Basin,Iran

Considering the importance of fractured rock aquifers in the hydrogeologic process,this research aimed to analyze the flow regime,internal degree of karstification,and estimate storage volume in fractured rock aquifers of the Germi Chai Basin in northwest Iran,which is attributed to its active tectonics,erosion,and the lithological diversity.Given the geological setting,the hypothesis is that this basin is characterized by a high degree of karstification and diffuse or intermediate flow regime leading to variation in discharge flow rate.The hydrodynamic and hadrochemical analysis was conducted on 9 well distributed springs across the basin from 2019 to 2020.The maximum flow rate in most of the springs appeared in the early wet season despite their different levels of fluctuations on the monthly discharge time series.Analyzing the spring recession curve form revealed an aquifer containing multiple micro-regimes withαrecession coefficients and a degree of karstification ranging between 0.001 to 0.06 and 0.55 to 2.61,respectively.These findings indicated a dominant diffuse and intermediate flow system resulting from the development of a high density of fractures in this area.The electrical conductivity of the spring changes inversely proportional to the change in flow discharge,indicating the reasonable hydrological response of the aquifer to rainfall events.Hydrograph analysis revealed that the delay time of spring discharge after rainfall events mostly varies between 10 to 30 days.The total dynamic storage volume of the spring for a given period(2019-2020)was estimated to be approximately 1324 million cubic meters reflecting the long-term drainage potential and high perdurability of dynamic storage.Estimating the maximum and minimum ratio revealed that the springs recharging system in Germi Chai Basin comes under the slow aquifers category.This finding provides valuable insight into the hydrogeological properties of fractured rock aquifers contributing to effective water management strategy.

Experimental Investigation on Condensate Revaporization During Gas Injection Development in Fractured Gas Condensate Reservoirs

The gas field in the Bohai Bay Basin is a fractured metamorphic buried-hill reservoir with dual-media characteristics. The retrograde vaporization mechanism observed in this type of gas condensate reservoir differs significantly from that observed in sand gas condensate reservoirs. However, studies on improving the recovery of fractured gas condensate reservoirs are limited;thus, the impact of retrograde vaporization on condensate within fractured metamorphic buried-hill reservoirs remains unclear. To address this gap, a series of gas injection experiments are conducted in pressure-volume-temperature(PVT) cells and long-cores to investigate the retrograde vaporization effect of condensate using different gas injection media in fractured gas condensate reservoirs. We analyze the variation in condensate volume, gas-to-oil ratio, and condensate recovery during gas injection and examine the influence of various gas injection media(CO_(2), N_(2), and dry gas) under different reservoir properties and varying gas injection times. The results demonstrate that the exchange of components between injected gas and condensate significantly influences condensate retrograde vaporization in the formation. Compared with dry gas injection and N_(2) injection,CO_(2) injection exhibits a superior retrograde vaporization effect. At a CO_(2) injection volume of 1 PV, the percentage shrinkage volume of condensate is 13.82%. Additionally, at the maximum retrograde condensation pressure, CO_(2) injection can increase the recovery of condensate by 22.4%. However, the condensate recovery is notably lower in fractured gas condensate reservoirs than in homogeneous reservoirs, owing to the creation of dominant gas channeling by fractures, which leads to decreased condensate recovery. Regarding gas injection timing, the effect of gas injection at reservoir pressure on improving condensate recovery is superior to that of gas injection at the maximum retrograde condensation pressure. This research provides valuable guidance for designing gas injection development plans and dynamic tracking adjustments for fractured gas condensate reservoirs.

Evaluation of Well Spacing for Primary Development of Fractured Horizontal Wells in Tight Sandstone Gas Reservoirs

Methods for horizontal well spacing calculation in tight gas reservoirs are still adversely affected by the complexity of related control factors,such as strong reservoir heterogeneity and seepage mechanisms.In this study,the stress sensitivity and threshold pressure gradient of various types of reservoirs are quantitatively evaluated through reservoir seepage experiments.On the basis of these experiments,a numerical simulation model(based on the special seepage mechanism)and an inverse dynamic reserve algorithm(with different equivalent drainage areas)were developed.The well spacing ranges of Classes I,II,and III wells in the Q gas field are determined to be 802–1,000,600–662,and 285–400 m,respectively,with their average ranges as 901,631,and 342.5 m,respectively.By considering both the pairs of parallel well groups and series well groups as examples,the reliability of the calculation results is verified.It is shown that the combination of the two models can reduce errors and provide accurate results.

The Impact of Deformable Surrounding Rocks on Naturally-Fractured Reservoir Simulation

For the case of a fractured reservoir surrounded by deformable rocks, the appropriateness and applicability of the two common methods of coupling of flow and deformation, explicit (coupled) and implicit (uncoupled) methods are investigated. The explicit formulation is capable of modelling surrounding media;while the implicit coupling is unable to do so as deformation vector does not appear as a primary variable in the formulation. The governing differential equations and the finite element approximation of the governing equations for each of the methods are presented. Spatial discretization is achieved using the Galerkin method, and temporal discretisation using the finite difference technique. In the explicit model, coupling between flow and deformation is captured through volumetric strain compatibility amongst the phases within the system. In the implicit model, this is achieved by defining the pore space storativity as a function of the formation compressibility and the compressibility of the fluid phases within the pore space. The impact of rock deformability on early, intermediate and late time responses of fractured reservoir is investigated through several numerical examples. Salient features of each formulation are discussed and highlighted. It is shown that the implicit model is unable to capture the constraining effects of a non-yielding, surrounding rock, leading to incorrect projections of reservoir production irrespective of the history matching strategy adopted.

Mechanisms of fracture propagation from multi-cluster using a phase field based HMD coupling model in fractured reservoir

Natural fractures(NFs)are common in shale and tight reservoirs,where staged multi-cluster fracturing of horizontal wells is a prevalent technique for reservoir stimulation.While NFs and stress interference are recognized as significant factors affecting hydraulic fracture(HF)propagation,the combined influence of these factors remains poorly understood.To address this knowledge gap,a novel coupled hydromechanical-damage(HMD)model based on the phase field method is developed to investigate the propagation of multi-cluster HFs in fractured reservoirs.The comprehensive energy functional and control functions are established,while incorporating dynamic fluid distribution between multiple perforation clusters and refined changes in rock mechanical parameters during hydraulic fracturing.The HMD coupled multi-cluster HF propagation model investigates various scenarios,including single HF and single NF,reservoir heterogeneity,single HF and NF clusters,and multi-cluster HFs with NF clusters.The results show that the HMD coupling model can accurately capture the impact of approach angle(θ),stress difference and cementation strength on the interaction of HF and NF.The criterion of the open and cross zones is not fixed.The NF angle(a)is not a decisive parameter to discriminate the interaction.According to the relationship between approach angle(θ)and NF angle(a),the contact relationship of HF can be divided into three categories(θ=a,θ<a,andθ>a).The connected NF can increase the complexity of HF by inducing it to form branch fracture,resulting in a fractal dimension of HF as high as2.1280 at angles of±45°.Inter-fracture interference from the heel to the toe of HF shows the phenomenon of no,strong and weak interference.Interestingly,under the influence of NFs,distant HFs from the injection can become dominant fractures.However,as a gradually increases,inter-fracture stress interference becomes the primary factor influencing HF propagation,gradually superseding the dominance of NF induced fractures.

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.

Pressure transient characteristics of non-uniform conductivity fractured wells in viscoelasticity polymer flooding based on oil-water two-phase flow

Polymer flooding in fractured wells has been extensively applied in oilfields to enhance oil recovery.In contrast to water,polymer solution exhibits non-Newtonian and nonlinear behavior such as effects of shear thinning and shear thickening,polymer convection,diffusion,adsorption retention,inaccessible pore volume and reduced effective permeability.Meanwhile,the flux density and fracture conductivity along the hydraulic fracture are generally non-uniform due to the effects of pressure distribution,formation damage,and proppant breakage.In this paper,we present an oil-water two-phase flow model that captures these complex non-Newtonian and nonlinear behavior,and non-uniform fracture characteristics in fractured polymer flooding.The hydraulic fracture is firstly divided into two parts:high-conductivity fracture near the wellbore and low-conductivity fracture in the far-wellbore section.A hybrid grid system,including perpendicular bisection(PEBI)and Cartesian grid,is applied to discrete the partial differential flow equations,and the local grid refinement method is applied in the near-wellbore region to accurately calculate the pressure distribution and shear rate of polymer solution.The combination of polymer behavior characterizations and numerical flow simulations are applied,resulting in the calculation for the distribution of water saturation,polymer concentration and reservoir pressure.Compared with the polymer flooding well with uniform fracture conductivity,this non-uniform fracture conductivity model exhibits the larger pressure difference,and the shorter bilinear flow period due to the decrease of fracture flow ability in the far-wellbore section.The field case of the fall-off test demonstrates that the proposed method characterizes fracture characteristics more accurately,and yields fracture half-lengths that better match engineering reality,enabling a quantitative segmented characterization of the near-wellbore section with high fracture conductivity and the far-wellbore section with low fracture conductivity.The novelty of this paper is the analysis of pressure performances caused by the fracture dynamics and polymer rheology,as well as an analysis method that derives formation and fracture parameters based on the pressure and its derivative curves.

Exploring pore-scale production characteristics of oil shale after CO_(2) huff‘n’puff in fractured shale with varied permeability

Recent studies have indicated that the injection of carbon dioxide(CO_(2))can lead to increased oil recovery in fractured shale reservoirs following natural depletion.Despite advancements in understanding mass exchange processes in subsurface formations,there remains a knowledge gap concerning the disparities in these processes between the matrix and fractures at the pore scale in formations with varying permeability.This study aims to experimentally investigate the CO_(2) diffusion behaviors and in situ oil recovery through a CO_(2) huff‘n’puff process in the Jimsar shale oil reservoir.To achieve this,we designed three matrix-fracture models with different permeabilities(0.074 mD,0.170 mD,and 0.466 mD)and experimented at 30 MPa and 91℃.The oil concentration in both the matrix and fracture was monitored using a low-field nuclear magnetic resonance(LF-NMR)technique to quantify in situ oil recovery and elucidate mass-exchange behaviors.The results showed that after three cycles of CO_(2) huff‘n’puff,the total recovery degree increased from 30.28%to 34.95%as the matrix permeability of the core samples increased from 0.074 to 0.466 mD,indicating a positive correlation between CO_(2) extraction efficiency and matrix permeability.Under similar fracture conditions,the increase in matrix permeability further promoted CO_(2) extraction efficiency during CO_(2) huff‘n’puff.Specifically,the increase in matrix permeability of the core had the greatest effect on the extraction of the first-cycle injection in large pores,which increased from 16.42%to 36.64%.The findings from our research provide valuable insights into the CO_(2) huff‘n’puff effects in different pore sizes following fracturing under varying permeability conditions,shedding light on the mechanisms of CO_(2)-enhanced oil recovery in fractured shale oil reservoirs.

A multi-mechanism numerical simulation model for CO_(2)-EOR and storage in fractured shale oil reservoirs

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.

Cemented vertebra and adjacent vertebra refractured in a chronic kidney disease-mineral and bone disorder patient: A case report

BACKGROUND Although percutaneous vertebral augmentation(PVA)is a commonly used procedure for treating vertebral compression fracture(VCF),the risk of vertebral refracture should be considered.Chronic kidney disease-mineral and bone disorder(CKD-MBD)is a systemic disease of mineral and bone metabolism.It is associated with an increased risk of fracture.Few studies have reported the use of PVA in patients with CKD-MBD.We herein report a rare case wherein the cemented vertebra and the adjacent vertebra refractured simultaneously in a CKD-MBD patient after PVA.CASE SUMMARY A 74-year-old man suffered from low back pain after taking a fall about 3 wk ago.According to physical examination,imaging and laboratory findings,diagnoses of T12 VCF,CKD-MBD,and chronic kidney disease stage 5 were established.He then received percutaneous vertebroplasty at T12 vertebra.Fourteen weeks later,he presented with T12 and L1 vertebral refractures caused by lumbar sprain.Once again,he was given PVA which was optimized for the refractured vertebrae.Although the short-term postoperative effect was satisfactory,he reported chronic low back pain again at the 3-month follow-up.CONCLUSION It is necessary that patients with CKD-MBD who have received PVA are aware of the adverse effects of CKD-MBD.It may increase the risk of vertebral refracture.Furthermore,the PVA surgical technique needs to be optimized according to the condition of the patient.The medium-and long-term effects of PVA remain uncertain in patients with CKD-MBD.

Groundwater flow through fractured rocks and seepage control in geotechnical engineering: Theories and practices

Groundwater flow through fractured rocks has been recognized as an important issue in many geotechnical engineering practices.Several key aspects of fundamental mechanisms,numerical modeling and engineering applications of flow in fractured rocks are discussed.First,the microscopic mechanisms of fluid flow in fractured rocks,especially under the complex conditions of non-Darcian flow,multiphase flow,rock dissolution,and particle transport,have been revealed through a com-bined effort of visualized experiments and theoretical analysis.Then,laboratory and field methods of characterizing hydraulic properties(e.g.intrinsic permeability,inertial permeability,and unsaturated flow parameters)of fractured rocks in different flow regimes have been proposed.Subsequently,high-performance numerical simulation approaches for large-scale modeling of groundwater flow in frac-tured rocks and aquifers have been developed.Numerical procedures for optimization design of seepage control systems in various settings have also been proposed.Mechanisms of coupled hydro-mechanical processes and control of flow-induced deformation have been discussed.Finally,three case studies are presented to illustrate the applications of the improved theoretical understanding,characterization methods,modeling approaches,and seepage and deformation control strategies to geotechnical engi-neering projects.

Laboratory insights into the effects of methane hydrate on the anisotropic joint elastic-electrical properties in fractured sandstones

Fractured hydrate-bearing reservoirs show significantly anisotropic geophysical properties. The joint application of seismic and electromagnetic explorations is expected to accurately assess hydrate resources in the fractured reservoirs. However, the anisotropic joint elastic-electrical properties in such reservoirs that are the key to the successful application of the joint explorations, remain poorly understood. To obtain such knowledge, we designed and implemented dedicated laboratory experiments to study the anisotropic joint elastic-electrical properties in fractured artificial silica sandstones (with fracture density of about 6.2%, porosity of approximately 25.7%, and mean grainsize of 0.089 mm) with evolving methane hydrate. The experimental results showed that the anisotropic compressional wave velocities respectively increased and decreased with the forming and dissociating hydrate, and the variation in the increasing trend and the decreasing extent of the velocity perpendicular to the fractures were more significant than that parallel to the fractures, respectively. The experimental results also showed that the overall decreasing trend of the electrical conductivity parallel to the fractures was steeper than that perpendicular to the fractures during hydrate formation, and the general variations of the two conductivities with complex trend were similar during hydrate dissociation. The variations in the elastic and electrical anisotropic parameters with forming and dissociating hydrate were also found to be distinct. Interpretation of the experimental results suggested that the hydrate binding to the grains evolved to bridge the surfaces of fractures when saturation exceeded 10% during hydrate formation, and the bridging hydrate gradually evolved to floating in fractures during dissociation. The experimental results further showed that the anisotropic velocities and electrical conductivities were correlated with approximately consistent trends of different slopes during hydrate formation, and the joint elastic-electrical anisotropic parameters exhibited a sharp peak at the hydrate saturation of about 10%. The results suggested that the anisotropic joint properties can be employed not only to accurately estimate hydrate saturation but also possibly to identify hydrate distribution in the fractures.

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.

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.

Grouting theories and technologies for the reinforcement of fractured rocks surrounding deep roadways

Grouting is an effective method to improve the integrity and stability of fractured rocks that surround deep roadways.After years of research and practice,various theories and a complete set of grouting technologies for deep roadways with fractured rocks have been developed and are widely applied in Chinese coal mining production.This paper systematically summarizes and analyzes the research results concerning the theory,design,materials,processes,and equipment for the grouting and reinforcement of fractured rocks surrounding deep roadways.Specifically,in terms of grouting methods,pregrouting,groutingwhile-excavation,and postgrouting methods are explored;in terms of grouting theory,backfill grouting,compaction grouting,infiltration grouting,and fracture grouting theories are studied.In addition,this paper also studies grouting borehole arrangement,water-cement ratio,grouting pressure,grouting volume,grout diffusion radius,and other grouting parameters and their determination methods.On this basis,this paper explores the physical and mechanical properties of organic and organic-inorganic composite grouting materials,and assess grouting reinforcement quality testing methods and instruments.Taken as the field cases,the application of pregrouting in front of heading faces,groutingwhile-excavation,and postgrouting in the Kouzidong coal mine are then introduced,and the effects of the grouting reinforcements are evaluated.This paper proposes a development direction for grouting technology based on problems existing in the grouting reinforcement of fractured rocks surrounding deep roadways.

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.