Multiscale modeling of gas-induced fracturing in anisotropic clayey rocks

In the context of repositories for nuclear waste,understanding the behavior of gas migration through clayey rocks with inherent anisotropy is crucial for assessing the safety of geological disposal facilities.The primary mechanism for gas breakthrough is the opening of micro-fractures due to high gas pressure.This occurs at gas pressures lower than the combined strength of the rock and its minimum principal stress under external loading conditions.To investigate the mechanism of microscale mode-I ruptures,it is essential to incorporate a multiscale approach that includes subcritical microcracks in the modeling framework.In this contribution,we derive the model from microstructures that contain periodically distributed microcracks within a porous material.The damage evolution law is coupled with the macroscopic poroelastic system by employing the asymptotic homogenization method and considering the inherent hydro-mechanical(HM)anisotropy at the microscale.The resulting permeability change induced by fracture opening is implicitly integrated into the gas flow equation.Verification examples are presented to validate the developed model step by step.An analysis of local macroscopic response is undertaken to underscore the influence of factors such as strain rate,initial damage,and applied stress,on the gas migration process.Numerical examples of direct tension tests are used to demonstrate the model’s efficacy in describing localized failure characteristics.Finally,the simulation results for preferential gas flow reveal the robustness of the two-scale model in explicitly depicting gas-induced fracturing in anisotropic clayey rocks.The model successfully captures the common behaviors observed in laboratory experiments,such as a sudden drop in gas injection pressure,rapid build-up of downstream gas pressure,and steady-state gas flow following gas breakthrough.

Impact of effective stress on permeability for carbonate fractured-vuggy rocks

To gain insight into the flow mechanisms and stress sensitivity for fractured-vuggy reservoirs,several core models with different structural characteristics were designed and fabricated to investigate the impact of effective stress on permeability for carbonate fractured-vuggy rocks(CFVR).It shows that the permeability performance curves under different pore and confining pressures(i.e.altered stress conditions)for the fractured core models and the vuggy core models have similar change patterns.The ranges of permeability variation are significantly wider at high pore pressures,indicating that permeability reduction is the most significant during the early stage of development for fractured-vuggy reservoirs.Since each obtained effective stress coefficient for permeability(ESCP)varies with the changes in confining pressure and pore pressure,the effective stresses for permeability of four representative CFVR show obvious nonlinear characteristics,and the variation ranges of ESCP are all between 0 and 1.Meanwhile,a comprehensive ESCP mathematical model considering triple media,including matrix pores,fractures,and dissolved vugs,was proposed.It is proved theoretically that the ESCP of CFVR generally varies between 0 and 1.Additionally,the regression results showed that the power model ranked highest among the four empirical models mainly applied in stress sensitivity characterization,followed by the logarithmic model,exponential model,and binomial model.The concept of“permeability decline rate”was introduced to better evaluate the stress sensitivity performance for CFVR,in which the one-fracture rock is the strongest,followed by the fracture-vug rock and two-horizontalfracture rock;the through-hole rock is the weakest.In general,this study provides a theoretical basis to guide the design of development and adjustment programs for carbonate fractured-vuggy reservoirs.

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.

Development and application of novel high‐efficiency composite ultrafine cement grouts for roadway in fractured surrounding rocks

The fractured surrounding rocks of roadways pose major challenges to safe mining.Grouting has often been used to reinforce the surrounding rocks to mitigate the safety risks associated with fractured rocks.The aim of this study is to develop highly efficient composite ultrafine cement(CUC)grouts to reinforce the roadway in fractured surrounding rocks.The materials used are ultrafine cement(UC),ultrafine fly ash(UF),ultrafine slag(US),and additives(superplasticizer[SUP],aluminate ultrafine expansion agent[AUA],gypsum,and retarder).The fluidity,bleeding,shrinkage,setting time,chemical composition,microstructure,degree of hydration,and mechanical property of grouting materials were evaluated in this study.Also,a suitable and effective CUC grout mixture was used to reinforce the roadway in the fractured surrounding rock.The results have shown that the addition of UF and US reduces the plastic viscosity of CUC,and the best fluidity can be obtained by adding 40%UF and 10%US.Since UC and UF particles are small,the pozzolanic effect of UF promotes the hydration reaction,which is conductive to the stability of CUC grouts.In addition,fine particles of UC,UF,and US can effectively fill the pores,while the volumetric expansion of AUA and gypsum decreases the pores and thus affects the microstructure of the solidified grout.The compressive test results have shown that the addition of specific amounts of UF and US can ameliorate the mechanical properties of CUC grouts.Finally,the CUC22‐8 grout was used to reinforce the No.20322 belt roadway.The results of numerical simulation and field monitoring have indicated that grouting can efficaciously reinforce the surrounding rock of the roadway.In this research,high‐performance CUC grouts were developed for surrounding rock reinforcement of underground engineering by utilizing UC and some additives.

A novel true triaxial test system for microwave-induced fracturing of hard rocks

This study introduces a test system for microwave-induced fracturing of hard rocks under true triaxial stress.The test system comprises a true triaxial stress loading system,an open-ended microwaveinduced fracturing system,a data acquisition system,an acoustic emission(AE)monitoring system,and an auxiliary specimen loading system.Microwave-induced surface and borehole fracturing tests under true triaxial stress were fulfilled for the first time,which overcomes the problem of microwave leakage in the coupling loading of true triaxial stress and microwave.By developing the dynamic monitoring system,the thermal response and fracture evolution were obtained during microwave irradiation.The monitoring system includes the infrared thermometry technique for monitoring rock surface temperature,the distributed optic fiber sensing technique for monitoring temperature in borehole in rock,the AE technique and two-dimensional digital speckle correlation technique for monitoring the evolution of thermal damage and the rock fracturing process.To validate the advantages of the test system and investigate the characteristics of microwave-induced fracturing of hard rocks,the study demonstrates the experimental methods and results for microwave-induced surface and borehole fracturing under true triaxial stress.The results show that thermal cracking presented intermittent characteristics(calm eactiveecalm)during microwave-induced surface and borehole fracturing of basalt.In addition,true triaxial stress can inhibit the development and distribution of thermal cracks during microwave-induced surface fracturing.When microwave-induced borehole fracturing occurs,it promotes the distribution of thermal cracks in rock,but inhibits the width of cracks.The results also prove the reliability of the test system.

Shear mechanical properties and fracturing responses of layered rough jointed rock-like materials

This study aims to investigate mechanical properties and failure mechanisms of layered rock with rough joint surfaces under direct shear loading.Cubic layered samples with dimensions of 100 mm×100 mm×100 mm were casted using rock-like materials,with anisotropic angle(α)and joint roughness coefficient(JRC)ranging from 15°to 75°and 2-20,respectively.The direct shear tests were conducted under the application of initial normal stress(σ_(n)) ranging from 1-4 MPa.The test results indicate significant differences in mechanical properties,acoustic emission(AE)responses,maximum principal strain fields,and ultimate failure modes of layered samples under different test conditions.The peak stress increases with the increasingαand achieves a maximum value atα=60°or 75°.As σ_(n) increases,the peak stress shows an increasing trend,with correlation coefficients R² ranging from 0.918 to 0.995 for the linear least squares fitting.As JRC increases from 2-4 to 18-20,the cohesion increases by 86.32%whenα=15°,while the cohesion decreases by 27.93%whenα=75°.The differences in roughness characteristics of shear failure surface induced byαresult in anisotropic post-peak AE responses,which is characterized by active AE signals whenαis small and quiet AE signals for a largeα.For a given JRC=6-8 andσ_(n)=1 MPa,asαincreases,the accumulative AE counts increase by 224.31%(αincreased from 15°to 60°),and then decrease by 14.68%(αincreased from 60°to 75°).The shear failure surface is formed along the weak interlayer whenα=15°and penetrates the layered matrix whenα=60°.Whenα=15°,as σ_(n) increases,the adjacent weak interlayer induces a change in the direction of tensile cracks propagation,resulting in a stepped pattern of cracks distribution.The increase in JRC intensifies roughness characteristics of shear failure surface for a smallα,however,it is not pronounced for a largeα.The findings will contribute to a better understanding of the mechanical responses and failure mechanisms of the layered rocks subjected to shear loads.

Effect of dynamic loading orientation on fracture properties of surrounding rocks in twin tunnels

For expedited transportation,vehicular tunnels are often designed as two adjacent tunnels,which frequently experience dynamic stress waves from various orientations during blasting excavation.To analyze the impact of dynamic loading orientation on the stability of the twin-tunnel,a split Hopkinson pressure bar(SHPB)apparatus was used to conduct a dynamic test on the twin-tunnel specimens.The two tunnels were rotated around the specimen’s center to consider the effect of dynamic loading orientation.LS-DYNA software was used for numerical simulation to reveal the failure properties and stress wave propagation law of the twin-tunnel specimens.The findings indicate that,for a twin-tunnel exposed to a dynamic load from different orientations,the crack initiation position appears most often at the tunnel corner,tunnel spandrel,and tunnel floor.As the impact direction is created by a certain angle(30°,45°,60°,120°,135°,and 150°),the fractures are produced in the middle of the line between the left tunnel corner and the right tunnel spandrel.As the impact loading angle(a)is 90°,the tunnel sustains minimal damage,and only tensile fractures form in the surrounding rocks.The orientation of the impact load could change the stress distribution in the twin-tunnel,and major fractures are more likely to form in areas where the tensile stress is concentrated.

Coupled hydro-thermo-mechanical modeling of hydraulic fracturing in quasi-brittle rocks using BPM-DEM

This paper presents an improved understanding of coupled hydro-thermo-mechanical(HTM) hydraulic fracturing of quasi-brittle rock using the bonded particle model(BPM) within the discrete element method(DEM). BPM has been recently extended by the authors to account for coupled convective econductive heat flow and transport, and to enable full hydro-thermal fluidesolid coupled modeling.The application of the work is on enhanced geothermal systems(EGSs), and hydraulic fracturing of hot dry rock(HDR) is studied in terms of the impact of temperature difference between rock and a flowing fracturing fluid. Micro-mechanical investigation of temperature and fracturing fluid effects on hydraulic fracturing damage in rocks is presented. It was found that fracture is shorter with pronounced secondary microcracking along the main fracture for the case when the convectiveeconductive thermal heat exchange is considered. First, the convection heat exchange during low-viscosity fluid infiltration in permeable rock around the wellbore causes significant rock cooling, where a finger-like fluid infiltration was observed. Second, fluid infiltration inhibits pressure rise during pumping and delays fracture initiation and propagation. Additionally, thermal damage occurs in the whole area around the wellbore due to rock cooling and cold fluid infiltration. The size of a damaged area around the wellbore increases with decreasing fluid dynamic viscosity. Fluid and rock compressibility ratio was found to have significant effect on the fracture propagation velocity.

Integrated hydraulic fracturing techniques to enhance oil recovery from tight rocks

Two main challenges exist in enhancing oil recovery rate from tight oil reservoirs,namely how to create an effective complicated fracture network and how to enhance the imbibition effect of fracturing fluid.In response to the challenges,through modeling experiment in laboratory and evaluation of field application results,a set of integrated efficient fracturing and enhanced oil recovery(EOR)techniques suitable for tight oil development in China has been proposed.(1)Fracturing with temporary plugging agents to realize stimulation in multiple clusters,to form dense fracture network,and thus maximizing the drainage area;(2)Supporting induced fractures with micro-sized proppants during the prepad fluid fracture-making stage,to generate dense fracture network with high conductivity;(3)Using the liquid nanofluid as a fracturing fluid additive to increase oil-water displacement ratio and take advantage of the massive injected fracturing fluid and maximize the oil production after hydraulic fracturing.

Investigating the Effect of Fatigue on Fracturing Resistance of Rocks

There is no doubt that an understanding of brittle rock fracturing is a key element in the solution of many engineering problems that involve rock structures. Some rock structures such as bridge and dam abutments and foundations, and tunnel walls, undergo both static and cyclic loading caused by drilling and blasting, and vehicle-induced vibrations. This type of loading often causes rock to fail at a lower than its static strength due to the effect of rock fatigue. A series of laboratory diametrical compression tests was performed on Brisbane tuff disc specimens to investigate their mode-I fracture toughness response to static and cyclic loading, as a function of the applied load. Both the static and cyclic loading tests were carried out on CCNBD （cracked chevron notched Brazilian disc） rock specimens. In the tests described herein, the reduction in fracture toughness under dynamic cyclic loading was found to be up to 48% of the static fracture toughness. Contrary to the static tests, the cyclic tests produced much more crushed material in front of the tip of the chevron notched crack.

Formation of fractures in carbonate rocks by pad acid fracturing with different states of carbon dioxide

Carbonate outcrops were taken from Ma 51 sub-member in the Lower Paleozoic in the Yan’an gas field to conduct true tri-axial hydraulic fracturing experiments with water, liquid CO_(2) and supercritical CO_(2). CT scan was applied to analyze initiation and propagation laws of hydraulic fractures in carbonate rocks. The experiments show that supercritical CO_(2) has low viscosity, strong diffusivity and large filtration during fracturing, which is more liable to increase pore pressure of rocks around wellbore and decrease breakdown pressure of carbonate rocks. However, it would cost much more volume of supercritical CO_(2) than water to fracture rocks since the former increases the wellbore pressure more slowly during fracturing. For carbonate rocks with few natural fractures, tensional fractures are generated by fracturing with water and liquid CO_(2), and these fractures propagate along the maximum horizontal principal stress direction;while fracturing with supercritical CO_(2) can form shear fractures, whose morphology is rarely influenced by horizontal stress difference. Besides, the angle between propagation direction of these shear fractures near the wellbore and the maximum horizontal principal stress is 45°, and the fractures would gradually turn to propagate along the maximum horizontal principal stress when they extend to a certain distance from the wellbore, leading to an increase of fracture tortuosity compared with the former. For carbonate rocks with well-developed natural fractures, fracturing with fresh water is conducive to connect natural fractures with low approaching angle and form stepped fractures with simple morphology. The key to forming complex fractures after fracturing carbonate rocks is to connect the natural fractures with high approaching angle. It is easier for liquid CO_(2) with low viscosity to realize such connection. Multi-directional fractures with relatively complex morphology would be formed after fracturing with liquid CO_(2).

Fracture evolution around pre-existing cylindrical cavities in brittle rocks under uniaxial compression

The development of fracture around pre-existing cylindrical cavities in brittle rocks was examined using physical models and acoustic emission technique. The experimental results indicate that when granite blocks containing one pre-existing cylindrical cavity are loaded in uniaxial compression condition, the profiles of cracks around the cavity can be characterized by tensile cracking （splitting parallel to the axial compression direction） at the roof-floor, compressive crack at two side walls, and remote or secondary cracks at the perimeter of the cavity. Moreover, fracture around cavity is size-dependent. In granite blocks containing pre-existing half-length cylindrical cavities, compressive stress concentration is found to initiate at the two sidewalls and induce shear crack propagation and coalescence. In granite blocks containing multiple parallel cylindrical cavities, the adjacent cylindrical cavities can influence each other and the eventual failure mode is determined by the interaction of tensile, compressive and shear stresses. Experimental results show that both tensile and compressive stresses play an important role in fracture evolution process around cavities in brittle rocks.

Numerical simulation of optimum mining design for high stress hard-rock deposit based on inducing fracturing mechanism

The 3D numerical simulation model of deep hard-rock deposit in Kaiyang Phosphate Mine of Guiyang was established based on the practical engineering using 3DEC numerical simulation software. The distribution characteristics of displacement fields and plastic zones of the orebody were simulated in three different excavation cases, including the case of excavation artificial inducted roadway in the orebody, the case of horizontal or vertical excavation direction and the case of the upward or downward excavation order. The simulation results indicate that the plastic zone and displacement field of surrounding rock around the inducted roadway are continuously increasing with the increase of the exposure time after digging an artificial inducted roadway in the orebody. Thus the raw rock ore becomes easier to be fragmented, which provides advantageous conditions for roadheader to cut high stress hard-rock. It is worthy noting that there is a large difference in effective utilization of deep ground pressure between horizontal and vertical excavation directions. The later can produce larger deformation and fracture zone than the former on the rock mass around the deduced roadway, which means that the later may utilize the high ground pressure more effectively to break hard-rock. And the obtained results also show that upward excavation order is more helpful for ground pressure to break rock than downward excavation order.

A method to interpret fracture aperture of rock slope using adaptive shape and unmanned aerial vehicle multi-angle nap-of-the-object photogrammetry

The aperture of natural rock fractures significantly affects the deformation and strength properties of rock masses,as well as the hydrodynamic properties of fractured rock masses.The conventional measurement methods are inadequate for collecting data on high-steep rock slopes in complex mountainous regions.This study establishes a high-resolution three-dimensional model of a rock slope using unmanned aerial vehicle(UAV)multi-angle nap-of-the-object photogrammetry to obtain edge feature points of fractures.Fracture opening morphology is characterized using coordinate projection and transformation.Fracture central axis is determined using vertical measuring lines,allowing for the interpretation of aperture of adaptive fracture shape.The feasibility and reliability of the new method are verified at a construction site of a railway in southeast Tibet,China.The study shows that the fracture aperture has a significant interval effect and size effect.The optimal sampling length for fractures is approximately 0.5e1 m,and the optimal aperture interpretation results can be achieved when the measuring line spacing is 1%of the sampling length.Tensile fractures in the study area generally have larger apertures than shear fractures,and their tendency to increase with slope height is also greater than that of shear fractures.The aperture of tensile fractures is generally positively correlated with their trace length,while the correlation between the aperture of shear fractures and their trace length appears to be weak.Fractures of different orientations exhibit certain differences in their distribution of aperture,but generally follow the forms of normal,log-normal,and gamma distributions.This study provides essential data support for rock and slope stability evaluation,which is of significant practical importance.

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.

A review of reservoir damage during hydraulic fracturing of deep and ultra-deep reservoirs

Deep and ultra-deep reservoirs have gradually become the primary focus of hydrocarbon exploration as a result of a series of significant discoveries in deep hydrocarbon exploration worldwide.These reservoirs present unique challenges due to their deep burial depth(4500-8882 m),low matrix permeability,complex crustal stress conditions,high temperature and pressure(HTHP,150-200℃,105-155 MPa),coupled with high salinity of formation water.Consequently,the costs associated with their exploitation and development are exceptionally high.In deep and ultra-deep reservoirs,hydraulic fracturing is commonly used to achieve high and stable production.During hydraulic fracturing,a substantial volume of fluid is injected into the reservoir.However,statistical analysis reveals that the flowback rate is typically less than 30%,leaving the majority of the fluid trapped within the reservoir.Therefore,hydraulic fracturing in deep reservoirs not only enhances the reservoir permeability by creating artificial fractures but also damages reservoirs due to the fracturing fluids involved.The challenging“three-high”environment of a deep reservoir,characterized by high temperature,high pressure,and high salinity,exacerbates conventional forms of damage,including water sensitivity,retention of fracturing fluids,rock creep,and proppant breakage.In addition,specific damage mechanisms come into play,such as fracturing fluid decomposition at elevated temperatures and proppant diagenetic reactions at HTHP conditions.Presently,the foremost concern in deep oil and gas development lies in effectively assessing the damage inflicted on these reservoirs by hydraulic fracturing,comprehending the underlying mechanisms,and selecting appropriate solutions.It's noteworthy that the majority of existing studies on reservoir damage primarily focus on conventional reservoirs,with limited attention given to deep reservoirs and a lack of systematic summaries.In light of this,our approach entails initially summarizing the current knowledge pertaining to the types of fracturing fluids employed in deep and ultra-deep reservoirs.Subsequently,we delve into a systematic examination of the damage processes and mechanisms caused by fracturing fluids within the context of hydraulic fracturing in deep reservoirs,taking into account the unique reservoir characteristics of high temperature,high pressure,and high in-situ stress.In addition,we provide an overview of research progress related to high-temperature deep reservoir fracturing fluid and the damage of aqueous fracturing fluids to rock matrix,both artificial and natural fractures,and sand-packed fractures.We conclude by offering a summary of current research advancements and future directions,which hold significant potential for facilitating the efficient development of deep oil and gas reservoirs while effectively mitigating reservoir damage.

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.

Discontinuity development patterns and the challenges for 3D discrete fracture network modeling on complicated exposed rock surfaces

Natural slopes usually display complicated exposed rock surfaces that are characterized by complex and substantial terrain undulation and ubiquitous undesirable phenomena such as vegetation cover and rockfalls.This study presents a systematic outcrop research of fracture pattern variations in a complicated rock slope,and the qualitative and quantitative study of the complex phenomena impact on threedimensional(3D)discrete fracture network(DFN)modeling.As the studies of the outcrop fracture pattern have been so far focused on local variations,thus,we put forward a statistical analysis of global variations.The entire outcrop is partitioned into several subzones,and the subzone-scale variability of fracture geometric properties is analyzed(including the orientation,the density,and the trace length).The results reveal significant variations in fracture characteristics(such as the concentrative degree,the average orientation,the density,and the trace length)among different subzones.Moreover,the density of fracture sets,which is approximately parallel to the slope surface,exhibits a notably higher value compared to other fracture sets across all subzones.To improve the accuracy of the DFN modeling,the effects of three common phenomena resulting from vegetation and rockfalls are qualitatively analyzed and the corresponding quantitative data processing solutions are proposed.Subsequently,the 3D fracture geometric parameters are determined for different areas of the high-steep rock slope in terms of the subzone dimensions.The results show significant variations in the same set of 3D fracture parameters across different regions with density differing by up to tenfold and mean trace length exhibiting differences of 3e4 times.The study results present precise geological structural information,improve modeling accuracy,and provide practical solutions for addressing complex outcrop issues.

Numerical modeling of fracture propagation of supercritical CO_(2)compound fracturing

The exploitation of shale gas is promising due to depletion of the conventional energy and intensification of the greenhouse effect.In this paper,we proposed a heat-fluid-solid coupling damage model of supercritical CO_(2)(SC-CO_(2))compound fracturing which is expected to be an efficient and environmentally friendly way to develop shale gas.The coupling model is solved by the finite element method,and the results are in good agreement with the analytical solutions and fracturing experiments.Based on this model,the fracture propagation characteristics at the two stages of compound fracturing are studied and the influence of pressurization rate,in situ stress,bedding angle,and other factors are considered.The results show that at the SC-CO_(2)fracturing stage,a lower pressurization rate is conducive to formation of the branches around main fractures,while a higher pressurization rate inhibits formation of the branches around main fractures and promotes formation of the main fractures.Both bedding and in situ stress play a dominant role in the fracture propagation.When the in situ stress ratio(δ_(x)/δ_(y))is 1,the presence of bedding can reduce the initiation pressure and failure pressure.Nevertheless,it will cause the fracture to propagate along the bedding direction,reducing the fracture complexity.In rocks without bedding,hydraulic fracturing has the lengthening and widening effects for SC-CO_(2)induced fracture.In shale,fractures induced at the hydraulic fracturing stage are more likely to be dominated by in situ stresses and have a shorter reorientation radius.Therefore,fracture branches propagating along the maximum principal stress direction may be generated around the main fractures induced by SC-CO_(2)at the hydraulic fracturing stage.When the branches converge with the main fractures,fracture zones are easily formed,and thus the fracture complexity and damage area can be significantly increased.The results are instructive for the design and application of SC-CO_(2)compound fracturing.

Predicting the probability distribution of Martian rocks mechanical property based on microscale rock mechanical experiments and accurate grain-based modeling

The exploration of Mars would heavily rely on Martian rocks mechanics and engineering technology.As the mechanical property of Martian rocks is uncertain,it is of utmost importance to predict the probability distribution of Martian rocks mechanical property for the success of Mars exploration.In this paper,a fast and accurate probability distribution method for predicting the macroscale elastic modulus of Martian rocks was proposed by integrating the microscale rock mechanical experiments(micro-RME),accurate grain-based modeling(AGBM)and upscaling methods based on reliability principles.Firstly,the microstructure of NWA12564 Martian sample and elastic modulus of each mineral were obtained by micro-RME with TESCAN integrated mineral analyzer(TIMA)and nanoindentation.The best probability distribution function of the minerals was determined by Kolmogorov-Smirnov(K-S)test.Secondly,based on best distribution function of each mineral,the Monte Carlo simulations(MCS)and upscaling methods were implemented to obtain the probability distribution of upscaled elastic modulus.Thirdly,the correlation between the upscaled elastic modulus and macroscale elastic modulus obtained by AGBM was established.The accurate probability distribution of the macroscale elastic modulus was obtained by this correlation relationship.The proposed method can predict the probability distribution of Martian rocks mechanical property with any size and shape samples.