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 ...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.展开更多
The geometric properties of fracture surfaces significantly influence shear-seepage in rock fractures,introducing complexities to fracture modelling.The present study focuses on the hydro-mechanical behaviours of roug...The geometric properties of fracture surfaces significantly influence shear-seepage in rock fractures,introducing complexities to fracture modelling.The present study focuses on the hydro-mechanical behaviours of rough rock fractures during shear-seepage processes to reveal how dilatancy and fracture asperities affect these phenomena.To achieve this,an improved shear-flow model(SFM)is proposed with the incorporation of dilatancy effect and asperities.In particular,shear dilatancy is accounted for in both the elastic and plastic stages,in contrast to some existing models that only consider it in the elastic stage.Depending on the computation approaches for the peak dilatancy angle,three different versions of the SFM are derived based on Mohr-Coulomb,joint roughness coefficient-joint compressive strength(JRC-JCS),and Grasselli’s theories.Notably,this is a new attempt that utilizes Grasselli’s model in shearseepage analysis.An advanced parameter optimization method is introduced to accurately determine model parameters,addressing the issue of local optima inherent in some conventional methods.Then,model performance is evaluated against existing experimental results.The findings demonstrate that the SFM effectively reproduces the shear-seepage characteristics of rock fracture across a wide range of stress levels.Further sensitivity analysis reveals how dilatancy and asperity affect hydraulic properties.The relation between hydro-mechanical properties(dilatancy displacement and hydraulic conductivity)and asperity parameters is analysed.Several profound understandings of the shear-seepage process are obtained by exploring the phenomenon under various conditions.展开更多
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 fracture...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.展开更多
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 accura...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.展开更多
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 st...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.展开更多
To achieve the loading of the stress path of hard rock,the spherical discrete element model(DEM)and the new flexible membrane technology were utilized to realize the transient loading of three principal stresses with ...To achieve the loading of the stress path of hard rock,the spherical discrete element model(DEM)and the new flexible membrane technology were utilized to realize the transient loading of three principal stresses with arbitrary magnitudes and orientations.Furthermore,based on the deep tunnel of China Jinping Underground Laboratory II(CJPL-II),the deformation and fracture evolution characteristics of deep hard rock induced by excavation stress path were analyzed,and the mechanisms of transient loading-unloading and stress rotation-induced fractures were revealed from a mesoscopic perspective.The results indicated that the stressestrain curve exhibits different trends and degrees of sudden changes when subjected to transient changes in principal stress,accompanied by sudden changes in strain rate.Stress rotation induces spatially directional deformation,resulting in fractures of different degrees and orientations,and increasing the degree of deformation anisotropy.The correlation between the degree of induced fracture and the unloading magnitude of minimum principal stress,as well as its initial level is significant and positive.The process of mechanical response during transient unloading exhibits clear nonlinearity and directivity.After transient unloading,both the minimum principal stress and minimum principal strain rate decrease sharply and then tend to stabilize.This occurs from the edge to the interior and from the direction of the minimum principal stress to the direction of the maximum principal stress on theε1-ε3 plane.Transient unloading will induce a tensile stress wave.The ability to induce fractures due to changes in principal stress magnitude,orientation and rotation paths gradually increases.The analysis indicates a positive correlation between the abrupt change amplitude of strain rate and the maximum unloading magnitude,which is determined by the magnitude and rotation of principal stress.A high tensile strain rate is more likely to induce fractures under low minimum principal stress.展开更多
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 ef...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.展开更多
Accurately picking P-and S-wave arrivals of microseismic(MS)signals in real-time directly influences the early warning of rock mass failure.A common contradiction between accuracy and computation exists in the current...Accurately picking P-and S-wave arrivals of microseismic(MS)signals in real-time directly influences the early warning of rock mass failure.A common contradiction between accuracy and computation exists in the current arrival picking methods.Thus,a real-time arrival picking method of MS signals is constructed based on a convolutional-recurrent neural network(CRNN).This method fully utilizes the advantages of convolutional layers and gated recurrent units(GRU)in extracting short-and long-term features,in order to create a precise and lightweight arrival picking structure.Then,the synthetic signals with field noises are used to evaluate the hyperparameters of the CRNN model and obtain an optimal CRNN model.The actual operation on various devices indicates that compared with the U-Net method,the CRNN method achieves faster arrival picking with less performance consumption.An application of large underground caverns in the Yebatan hydropower station(YBT)project shows that compared with the short-term average/long-term average(STA/LTA),Akaike information criterion(AIC)and U-Net methods,the CRNN method has the highest accuracy within four sampling points,which is 87.44%for P-wave and 91.29%for S-wave,respectively.The sum of mean absolute errors(MAESUM)of the CRNN method is 4.22 sampling points,which is lower than that of the other methods.Among the four methods,the MS sources location calculated based on the CRNN method shows the best consistency with the actual failure,which occurs at the junction of the shaft and the second gallery.Thus,the proposed method can pick up P-and S-arrival accurately and rapidly,providing a reference for rock failure analysis and evaluation in engineering applications.展开更多
Rock slope with horizontal-layered fractured structure(HLFS)has high stability in its natural state.However,a strong earthquake can induce rock fissure expansion,ultimately leading to slope failure.In this study,the d...Rock slope with horizontal-layered fractured structure(HLFS)has high stability in its natural state.However,a strong earthquake can induce rock fissure expansion,ultimately leading to slope failure.In this study,the dynamic response,failure mode,and spectral characteristics of rock slope with HLFS under strong earthquake conditions were investigated based on the large-scale shaking table model test.On this basis,multiple sets of numerical calculation models were further established by UDEC discrete element program.Five influencing factors were considered in the parametric study of numerical simulations,including slope height,slope angle,bedding-plane spacing and secondary joint spacing as well as bedrock dip angle.The results showed that the failure process of rock slope with HLFS under earthquake action is mainly divided into four phases,i.e.,the tensile crack of the slope shoulder joints and shear dislocation at the top bedding plane,the extension of vertical joint cracks and increase of shear displacement,the formation of step-through sliding surfaces and the instability,and finally collapse of fractured rock mass.The acceleration response of slopes exhibits elevation amplification effect and surface effect.Numerical simulations indicate that the seismic stability of slopes with HLFS exhibits a negative correlation with slope height and angle,but a positive correlation with bedding-plane spacing,joint spacing,and bedrock dip angle.The results of this study can provide a reference for seismic stability evaluation of weathered rock slopes.展开更多
Due to the presence of ice and unfrozen water in pores of frozen rock,the rock fracture behaviors are susceptible to temperature.In this study,the potential thawing-induced softening effects on the fracture behaviors ...Due to the presence of ice and unfrozen water in pores of frozen rock,the rock fracture behaviors are susceptible to temperature.In this study,the potential thawing-induced softening effects on the fracture behaviors of frozen rock is evaluated by testing the tension fracture toughness(KIC)of frozen rock at different temperatures(i.e.-20℃,-15℃,-12℃,-10℃,-8℃,-6℃,-4℃,-2℃,and 0℃).Acoustic emission(AE)and digital image correlation(DIC)methods are utilized to analyze the microcrack propagation during fracturing.The melting of pore ice is measured using nuclear magnetic resonance(NMR)method.The results indicate that:(1)The KIC of frozen rock decreases moderately between-20℃ and-4℃,and rapidly between-4℃ and 0℃.(2)At-20℃ to-4℃,the fracturing process,deduced from the DIC results at the notch tip,exhibits three stages:elastic deformation,microcrack propagation and microcrack coalescence.However,at-4℃e0℃,only the latter two stages are observed.(3)At-4℃e0℃,the AE activities during fracturing are less than that at-20℃ to-4℃,while more small events are reported.(4)The NMR results demonstrate a reverse variation trend in pore ice content with increasing temperature,that is,a moderate decrease is followed by a sharp decrease and-4℃ is exactly the critical temperature.Next,we interpret the thawing-induced softening effect by linking the evolution in microscopic structure of frozen rock with its macroscopic fracture behaviors as follow:from-20℃ to-4℃,the thickening of the unfrozen water film diminishes the cementation strength between ice and rock skeleton,leading to the decrease in fracture parameters.From-4℃ to 0℃,the cementation effect of ice almost vanishes,and the filling effect of pore ice is reduced significantly,which facilitates microcrack propagation and thus the easier fracture of frozen rocks.展开更多
Identifying the real fracture of rock hidden in acoustic emission(AE)source clusters(AE-depicted microcrack zone)remains challenging and crucial.Here we revealed the AE energy(representing dissipated energy)distributi...Identifying the real fracture of rock hidden in acoustic emission(AE)source clusters(AE-depicted microcrack zone)remains challenging and crucial.Here we revealed the AE energy(representing dissipated energy)distribution rule in the rock microcrack zone and proposed an AE-energy-based method for identifying the real fracture.(1)A set of fracture experiments were performed on granite using wedgeloading,and the fracture process was detected and recorded by AE.The microcrack zone associated with the energy dissipation was characterized by AE sources and energy distribution,utilizing our selfdeveloped AE analysis program(RockAE).(2)The accumulated AE energy,an index representing energy dissipation,across the AE-depicted microcrack zone followed the normal distribution model(the mean and variance relate to the real fracture path and the microcrack zone width).This result implies that the nucleation and coalescence of massive cracks(i.e.,real fracture generation process)are supposed to follow a normal distribution.(3)Then,we obtained the real fracture extension path by joining the peak positions of the AE energy normal distribution curve at different cross-sections of the microcrack zone.Consequently,we distinguished between the microcrack zone and the concealed real fracture within it.The deviation was validated as slight as 1–3 mm.展开更多
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...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.展开更多
Photogrammetry,reconstructing three-dimensional(3D)models from overlapping two-dimensional(2D)photos,finds application in rock mechanics and rock engineering to extract geometrical details of reconstructed objects,for...Photogrammetry,reconstructing three-dimensional(3D)models from overlapping two-dimensional(2D)photos,finds application in rock mechanics and rock engineering to extract geometrical details of reconstructed objects,for example rock fractures.Fracture properties are important for determining the mechanical stability,permeability,strength,and shear behavior of the rock mass.Photogrammetry can be used to reconstruct detailed 3D models of two separated rock fracture surfaces to characterize fracture roughness and physical aperture,which controls the fluid flow,hydromechanical and shear behavior of the rock mass.This research aimed to determine the optimal number of scale bars required to produce high-precision 3D models of a fracture surface.A workflow has been developed to define the physical aperture of a fracture using photogrammetry.Three blocks of Kuru granite(25 cm×25 cm×10 cm)with an artificially induced fracture,were investigated.For scaling 3D models,321 markers were used as ground control points(GCPs)with predefined distances on each block.When the samples were wellmatched in their original positions,the entire block was photographed.Coordinate data of the GCPs were extracted from the 3D model of the blocks.Each half was surveyed separately and georeferenced by GCPs and merged into the same coordinate system.Two fracture surfaces were extracted from the 3D models and the vertical distance between the two surfaces was digitally calculated as physical aperture.Accuracy assessment of the photogrammetric reconstruction showed a 20-30 mm digital control distance accuracy when compared to known distances defined between markers.To attain this accuracy,the study found that at least 200 scale bars were required.Furthermore,photogrammetry was employed to measure changes in aperture under normal stresses.The results obtained from this approach were found to be in good agreement with those obtained using linear variable displacement transducers(LVDTs),with differences ranging from 1 mm to 8μm.展开更多
The surrounding rock of roadways exhibits intricate characteristics of discontinuity and heterogeneity.To address these complexities,this study employs non-local Peridynamics(PD)theory and reconstructs the kernel func...The surrounding rock of roadways exhibits intricate characteristics of discontinuity and heterogeneity.To address these complexities,this study employs non-local Peridynamics(PD)theory and reconstructs the kernel function to represent accurately the spatial decline of long-range force.Additionally,modifications to the traditional bondbased PD model are made.By considering the micro-structure of coal-rock materials within a uniform discrete model,heterogeneity characterized by bond random pre-breaking is introduced.This approach facilitates the proposal of a novel model capable of handling the random distribution characteristics of material heterogeneity,rendering the PD model suitable for analyzing the deformation and failure of heterogeneous layered coal-rock mass structures.The established numerical model and simulation method,termed the sub-homogeneous PD model,not only incorporates the support effect but also captures accurately the random heterogeneous micro-structure of roadway surrounding rock.The simulation results obtained using this model show good agreement with field measurements from the Fucun coal mine,effectively validating the model’s capability in accurately reproducing the deformation and failure mode of surrounding rock under bolt-supported(anchor cable).The proposed subhomogeneous PD model presents a valuable and effective simulation tool for studying the deformation and failure of roadway surrounding rock in coal mines,offering new insights and potential advancements.展开更多
Complexities in mechanical behaviours of rock masses mainly stem from inherent discontinuities,which calls for advanced bolt-grouting techniques for stability enhancement.Understanding the mechanical properties of bol...Complexities in mechanical behaviours of rock masses mainly stem from inherent discontinuities,which calls for advanced bolt-grouting techniques for stability enhancement.Understanding the mechanical properties of bolt-grouted fractured rock mass(BGFR)and developing accurate prediction methods are crucial to optimize the BGFR support strategies.This paper establishes a new elastoplastic(E-P)model based on the orthotropic and the Mohr-Coulomb(M-C)plastic-yielding criteria.The elastic parameters of the model were derived through a meso-mechanical analysis of composite materials mechanics(CMM).Laboratory BGFR specimens were prepared and uniaxial compression test and variable-angle shear test considering different bolt arrangements were carried out to obtain the mechanical parameters of the specimens.Results showed that the anisotropy of BGFR mainly depends on the relative volume content of each component material in a certain direction.Moreover,the mechanical parameters deduced from the theory of composite materials which consider the short fibre effect are shown to be in good agreement with those determined by laboratory experiments,and the variation rules maintained good consistency.Last,a case study of a real tunnel project is provided to highlight the effectiveness,validity and robustness of the developed E-P model in prediction of stresses and deformations.展开更多
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 aquifer...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.展开更多
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 ...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.展开更多
Rock fracture mechanics and accurate characterization of rock fracture are crucial for understanding a variety of phenomena interested in geological engineering and geoscience.These phenomena range from very large-sca...Rock fracture mechanics and accurate characterization of rock fracture are crucial for understanding a variety of phenomena interested in geological engineering and geoscience.These phenomena range from very large-scale asymmetrical fault structures to the scale of engineering projects and laboratory-scale rock fracture tests.Comprehensive study can involve mechanical modeling,site or post-mortem investigations,and inspection on the point cloud of the source locations in the form of earthquake,microseismicity,or acoustic emission.This study presents a comprehensive data analysis on characterizing the forming of the asymmetrical damage zone around a laboratory mixed-mode rock fracture.We substantiate the presence of asymmetrical damage through qualitative analysis and demonstrate that measurement uncertainties cannot solely explain the observed asymmetry.The implications of this demonstration can be manifold.On a larger scale,it solidifies a mechanical model used for explaining the contribution of aseismic mechanisms to asymmetrical fault structures.On a laboratory scale,it exemplifies an alternative approach to understanding the observational difference between the source location and the in situ or post-mortem inspection on the rock fracture path.The mechanical model and the data analysis can be informative to the interpretations of other engineering practices as well,but may face different types of challenges.展开更多
In practical engineering applications,rock mass are often found to be subjected to a triaxial stress state.Concurrently,defects like joints and fractures have a notable impact on the mechanical behavior of rock mass.S...In practical engineering applications,rock mass are often found to be subjected to a triaxial stress state.Concurrently,defects like joints and fractures have a notable impact on the mechanical behavior of rock mass.Such defects are identified as crucial contributors to the failure and instability of the surrounding rock,subsequently impacting the engineering stability.The study aimed to investigate the impact of fracture geometry and confining pressure on the deformation,failure characteristics,and strength of specimens using sand powder 3D printing technology and conventional triaxial compression tests.The results indicate that the number of fractures present considerably influences the peak strength,axial peak strain and elastic modulus of the specimens.Confining pressure is an important factor affecting the failure pattern of the specimen,under which the specimen is more prone to shear failure,but the initiation,expansion and penetration processes of secondary cracks in different fracture specimens are different.This study confirmed the feasibility of using sand powder 3D printing specimens as soft rock analogs for triaxial compression research.The insights from this research are deemed essential for a deeper understanding of the mechanical behavior of fractured surrounding rocks when under triaxial stress state.展开更多
This study employs similar simulation testing and discrete element simulation coupling to analyze the failure and deformation processes of a model coal seam's roof.The caving area of the overburden rock is divided...This study employs similar simulation testing and discrete element simulation coupling to analyze the failure and deformation processes of a model coal seam's roof.The caving area of the overburden rock is divided into three zones:the delamination fracture zone,broken fracture zone,and compaction zone.The caving and fracture zones'heights are approximately 110 m above the coal seam,with a maximum subsidence of 11 m.The delamination fracture zone's porosity range is between 0.2 and 0.3,while the remainder of the roof predominantly exhibits a porosity of less than 0.1.In addition,the numerical model's stress analysis revealed that the overburden rock's displacement zone forms an'arch-beam'structure starting from 160 m,with the maximum and minimum stress values decreasing as the distance of advancement increases.In the stress beam interval of the overburden rock,the maximum value changes periodically as the advancement distance increases.Based on a comparative analysis between observable data from on-site work and numerical simulation results,the stress data from the numerical simulation are essentially consistent with the actual results detected on-site,indicating the validity of the numerical simulation results.展开更多
基金This work was supported by the National Nature Science Foundation of China(Grant Nos.42177139 and 41941017)the Natural Science Foundation Project of Jilin Province,China(Grant No.20230101088JC).The authors would like to thank the anonymous reviewers for their comments and suggestions.
文摘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.
基金support from the National Natural Science Foundation of China(Grant Nos.51991392 and 42293355).
文摘The geometric properties of fracture surfaces significantly influence shear-seepage in rock fractures,introducing complexities to fracture modelling.The present study focuses on the hydro-mechanical behaviours of rough rock fractures during shear-seepage processes to reveal how dilatancy and fracture asperities affect these phenomena.To achieve this,an improved shear-flow model(SFM)is proposed with the incorporation of dilatancy effect and asperities.In particular,shear dilatancy is accounted for in both the elastic and plastic stages,in contrast to some existing models that only consider it in the elastic stage.Depending on the computation approaches for the peak dilatancy angle,three different versions of the SFM are derived based on Mohr-Coulomb,joint roughness coefficient-joint compressive strength(JRC-JCS),and Grasselli’s theories.Notably,this is a new attempt that utilizes Grasselli’s model in shearseepage analysis.An advanced parameter optimization method is introduced to accurately determine model parameters,addressing the issue of local optima inherent in some conventional methods.Then,model performance is evaluated against existing experimental results.The findings demonstrate that the SFM effectively reproduces the shear-seepage characteristics of rock fracture across a wide range of stress levels.Further sensitivity analysis reveals how dilatancy and asperity affect hydraulic properties.The relation between hydro-mechanical properties(dilatancy displacement and hydraulic conductivity)and asperity parameters is analysed.Several profound understandings of the shear-seepage process are obtained by exploring the phenomenon under various conditions.
基金supported by the National Natural Science Foundation of China(No.42307258)the technological research projects in Sichuan Province(No.2022YFSY0007)the China Atomic Energy Authority(CAEA)through the Geological Disposal Program.
文摘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.
基金supported by the National Natural Science Foundation of China(Grant No.42277165)the Fundamental Research Funds for the Central Universities,China University of Geosciences(Wuhan)(Grant No.CUGCJ1821)the National Overseas Study Fund(Grant No.202106410040).
文摘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.
基金supported by the National Key Research and Development Program of China(Grant No.2022YFC3080200)the National Natural Science Foundation of China(Grant No.42022053)the China Postdoctoral Science Foundation(Grant No.2023M731264).
文摘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.
基金the financial support from the National Natural Science Foundation of China(Grant No.51839003)Liaoning Revitalization Talents Program(Grant No.XLYCYSZX 1902)Hubei Key Laboratory for Efficient Utilization and Agglomeration of Metallurgic Mineral Resources(Grant No.2023zy002).
文摘To achieve the loading of the stress path of hard rock,the spherical discrete element model(DEM)and the new flexible membrane technology were utilized to realize the transient loading of three principal stresses with arbitrary magnitudes and orientations.Furthermore,based on the deep tunnel of China Jinping Underground Laboratory II(CJPL-II),the deformation and fracture evolution characteristics of deep hard rock induced by excavation stress path were analyzed,and the mechanisms of transient loading-unloading and stress rotation-induced fractures were revealed from a mesoscopic perspective.The results indicated that the stressestrain curve exhibits different trends and degrees of sudden changes when subjected to transient changes in principal stress,accompanied by sudden changes in strain rate.Stress rotation induces spatially directional deformation,resulting in fractures of different degrees and orientations,and increasing the degree of deformation anisotropy.The correlation between the degree of induced fracture and the unloading magnitude of minimum principal stress,as well as its initial level is significant and positive.The process of mechanical response during transient unloading exhibits clear nonlinearity and directivity.After transient unloading,both the minimum principal stress and minimum principal strain rate decrease sharply and then tend to stabilize.This occurs from the edge to the interior and from the direction of the minimum principal stress to the direction of the maximum principal stress on theε1-ε3 plane.Transient unloading will induce a tensile stress wave.The ability to induce fractures due to changes in principal stress magnitude,orientation and rotation paths gradually increases.The analysis indicates a positive correlation between the abrupt change amplitude of strain rate and the maximum unloading magnitude,which is determined by the magnitude and rotation of principal stress.A high tensile strain rate is more likely to induce fractures under low minimum principal stress.
基金This work was supported by the Joint Fund of NSFC for Enterprise Innovation and Development(Grant No.U19B6003-02-06)the National Natural Science Foundation of China(Grant No.51974331)+1 种基金the Natural Science Foundation of Jiangsu Province(Grant No.BK20200525)The authors would like to sincerely acknowledge these funding programs for their financial support.Particularly,the support provided by the China Scholarship Council(CSC)during a visit of Ke Sun(File No.202106440065)to the University of Alberta is also sincerely acknowledged.
文摘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.
基金We acknowledge the funding support from National Natural Science Foundation of China(Grant No.42077263).
文摘Accurately picking P-and S-wave arrivals of microseismic(MS)signals in real-time directly influences the early warning of rock mass failure.A common contradiction between accuracy and computation exists in the current arrival picking methods.Thus,a real-time arrival picking method of MS signals is constructed based on a convolutional-recurrent neural network(CRNN).This method fully utilizes the advantages of convolutional layers and gated recurrent units(GRU)in extracting short-and long-term features,in order to create a precise and lightweight arrival picking structure.Then,the synthetic signals with field noises are used to evaluate the hyperparameters of the CRNN model and obtain an optimal CRNN model.The actual operation on various devices indicates that compared with the U-Net method,the CRNN method achieves faster arrival picking with less performance consumption.An application of large underground caverns in the Yebatan hydropower station(YBT)project shows that compared with the short-term average/long-term average(STA/LTA),Akaike information criterion(AIC)and U-Net methods,the CRNN method has the highest accuracy within four sampling points,which is 87.44%for P-wave and 91.29%for S-wave,respectively.The sum of mean absolute errors(MAESUM)of the CRNN method is 4.22 sampling points,which is lower than that of the other methods.Among the four methods,the MS sources location calculated based on the CRNN method shows the best consistency with the actual failure,which occurs at the junction of the shaft and the second gallery.Thus,the proposed method can pick up P-and S-arrival accurately and rapidly,providing a reference for rock failure analysis and evaluation in engineering applications.
基金supported by Central Guiding Local Science and Technology Development Special Fund Project(No.ZYYD2023B02)the National Natural Science Foundation of China(Nos.52078432 and 52168066)the Scientific Research Project of China Railway First Survey and Design Institute Group Co.(No.20-06).
文摘Rock slope with horizontal-layered fractured structure(HLFS)has high stability in its natural state.However,a strong earthquake can induce rock fissure expansion,ultimately leading to slope failure.In this study,the dynamic response,failure mode,and spectral characteristics of rock slope with HLFS under strong earthquake conditions were investigated based on the large-scale shaking table model test.On this basis,multiple sets of numerical calculation models were further established by UDEC discrete element program.Five influencing factors were considered in the parametric study of numerical simulations,including slope height,slope angle,bedding-plane spacing and secondary joint spacing as well as bedrock dip angle.The results showed that the failure process of rock slope with HLFS under earthquake action is mainly divided into four phases,i.e.,the tensile crack of the slope shoulder joints and shear dislocation at the top bedding plane,the extension of vertical joint cracks and increase of shear displacement,the formation of step-through sliding surfaces and the instability,and finally collapse of fractured rock mass.The acceleration response of slopes exhibits elevation amplification effect and surface effect.Numerical simulations indicate that the seismic stability of slopes with HLFS exhibits a negative correlation with slope height and angle,but a positive correlation with bedding-plane spacing,joint spacing,and bedrock dip angle.The results of this study can provide a reference for seismic stability evaluation of weathered rock slopes.
基金We acknowledge the funding support from the National Natural Science Foundation of China(Grant No.42271148).
文摘Due to the presence of ice and unfrozen water in pores of frozen rock,the rock fracture behaviors are susceptible to temperature.In this study,the potential thawing-induced softening effects on the fracture behaviors of frozen rock is evaluated by testing the tension fracture toughness(KIC)of frozen rock at different temperatures(i.e.-20℃,-15℃,-12℃,-10℃,-8℃,-6℃,-4℃,-2℃,and 0℃).Acoustic emission(AE)and digital image correlation(DIC)methods are utilized to analyze the microcrack propagation during fracturing.The melting of pore ice is measured using nuclear magnetic resonance(NMR)method.The results indicate that:(1)The KIC of frozen rock decreases moderately between-20℃ and-4℃,and rapidly between-4℃ and 0℃.(2)At-20℃ to-4℃,the fracturing process,deduced from the DIC results at the notch tip,exhibits three stages:elastic deformation,microcrack propagation and microcrack coalescence.However,at-4℃e0℃,only the latter two stages are observed.(3)At-4℃e0℃,the AE activities during fracturing are less than that at-20℃ to-4℃,while more small events are reported.(4)The NMR results demonstrate a reverse variation trend in pore ice content with increasing temperature,that is,a moderate decrease is followed by a sharp decrease and-4℃ is exactly the critical temperature.Next,we interpret the thawing-induced softening effect by linking the evolution in microscopic structure of frozen rock with its macroscopic fracture behaviors as follow:from-20℃ to-4℃,the thickening of the unfrozen water film diminishes the cementation strength between ice and rock skeleton,leading to the decrease in fracture parameters.From-4℃ to 0℃,the cementation effect of ice almost vanishes,and the filling effect of pore ice is reduced significantly,which facilitates microcrack propagation and thus the easier fracture of frozen rocks.
基金supported by the National Natural Science Foundation of China(No.52274013)the Fundamental Research Funds for the Central Universities(No.2024ZDPYYQ1005)+1 种基金the National Key Research and Development Program of China(No.2021YFC2902103)the Independent Research Project of State Key Laboratory for Fine Exploration and Intelligent Development of Coal Resources,CUMT(No.SKLCRSM23X002).
文摘Identifying the real fracture of rock hidden in acoustic emission(AE)source clusters(AE-depicted microcrack zone)remains challenging and crucial.Here we revealed the AE energy(representing dissipated energy)distribution rule in the rock microcrack zone and proposed an AE-energy-based method for identifying the real fracture.(1)A set of fracture experiments were performed on granite using wedgeloading,and the fracture process was detected and recorded by AE.The microcrack zone associated with the energy dissipation was characterized by AE sources and energy distribution,utilizing our selfdeveloped AE analysis program(RockAE).(2)The accumulated AE energy,an index representing energy dissipation,across the AE-depicted microcrack zone followed the normal distribution model(the mean and variance relate to the real fracture path and the microcrack zone width).This result implies that the nucleation and coalescence of massive cracks(i.e.,real fracture generation process)are supposed to follow a normal distribution.(3)Then,we obtained the real fracture extension path by joining the peak positions of the AE energy normal distribution curve at different cross-sections of the microcrack zone.Consequently,we distinguished between the microcrack zone and the concealed real fracture within it.The deviation was validated as slight as 1–3 mm.
基金financially supported by National Natural Science Foundation of China(Grant Nos.52074312 and 52211530097)CNPC Science and Technology Innovation Foundation(Grant No.2021DQ02-0505).
文摘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.
基金funding provided by the State Nuclear Waste Management Fund(VYR)and the support of the Ministry of Economic Affairs and Employment of Finland on the Finnish Research Program on Nuclear Waste Management KYT2018 and KYT2022 of the Nuclear Energy Act(990/1987)in the research projects Fluid flow in fractured hard rock mass(RAKKA),funding numbers KYT 1/2021 and KYT 1/2022Additional support was received from the National Nuclear Safety and Waste Management Research Program SAFER2028,funding numbers SAFER 25/2023(MIRKA)and SAFER 42/2023(CORF).
文摘Photogrammetry,reconstructing three-dimensional(3D)models from overlapping two-dimensional(2D)photos,finds application in rock mechanics and rock engineering to extract geometrical details of reconstructed objects,for example rock fractures.Fracture properties are important for determining the mechanical stability,permeability,strength,and shear behavior of the rock mass.Photogrammetry can be used to reconstruct detailed 3D models of two separated rock fracture surfaces to characterize fracture roughness and physical aperture,which controls the fluid flow,hydromechanical and shear behavior of the rock mass.This research aimed to determine the optimal number of scale bars required to produce high-precision 3D models of a fracture surface.A workflow has been developed to define the physical aperture of a fracture using photogrammetry.Three blocks of Kuru granite(25 cm×25 cm×10 cm)with an artificially induced fracture,were investigated.For scaling 3D models,321 markers were used as ground control points(GCPs)with predefined distances on each block.When the samples were wellmatched in their original positions,the entire block was photographed.Coordinate data of the GCPs were extracted from the 3D model of the blocks.Each half was surveyed separately and georeferenced by GCPs and merged into the same coordinate system.Two fracture surfaces were extracted from the 3D models and the vertical distance between the two surfaces was digitally calculated as physical aperture.Accuracy assessment of the photogrammetric reconstruction showed a 20-30 mm digital control distance accuracy when compared to known distances defined between markers.To attain this accuracy,the study found that at least 200 scale bars were required.Furthermore,photogrammetry was employed to measure changes in aperture under normal stresses.The results obtained from this approach were found to be in good agreement with those obtained using linear variable displacement transducers(LVDTs),with differences ranging from 1 mm to 8μm.
基金supported by the National Natural Science Foundation of China(Nos.12302264,52104004,12072170,and 12202225)the Natural Science Foundation of Shandong Province(No.ZR2021QA042)Special Fund for Taishan Scholar Project(No.Tsqn202211180).
文摘The surrounding rock of roadways exhibits intricate characteristics of discontinuity and heterogeneity.To address these complexities,this study employs non-local Peridynamics(PD)theory and reconstructs the kernel function to represent accurately the spatial decline of long-range force.Additionally,modifications to the traditional bondbased PD model are made.By considering the micro-structure of coal-rock materials within a uniform discrete model,heterogeneity characterized by bond random pre-breaking is introduced.This approach facilitates the proposal of a novel model capable of handling the random distribution characteristics of material heterogeneity,rendering the PD model suitable for analyzing the deformation and failure of heterogeneous layered coal-rock mass structures.The established numerical model and simulation method,termed the sub-homogeneous PD model,not only incorporates the support effect but also captures accurately the random heterogeneous micro-structure of roadway surrounding rock.The simulation results obtained using this model show good agreement with field measurements from the Fucun coal mine,effectively validating the model’s capability in accurately reproducing the deformation and failure mode of surrounding rock under bolt-supported(anchor cable).The proposed subhomogeneous PD model presents a valuable and effective simulation tool for studying the deformation and failure of roadway surrounding rock in coal mines,offering new insights and potential advancements.
基金funded by the National Key Research and Development Plan(No.2022YFC3203200)Department of Science and Technology of Guangdong Province(No.2021ZT09G087)the National Natural Science Foundation Project of China(No.42167025).
文摘Complexities in mechanical behaviours of rock masses mainly stem from inherent discontinuities,which calls for advanced bolt-grouting techniques for stability enhancement.Understanding the mechanical properties of bolt-grouted fractured rock mass(BGFR)and developing accurate prediction methods are crucial to optimize the BGFR support strategies.This paper establishes a new elastoplastic(E-P)model based on the orthotropic and the Mohr-Coulomb(M-C)plastic-yielding criteria.The elastic parameters of the model were derived through a meso-mechanical analysis of composite materials mechanics(CMM).Laboratory BGFR specimens were prepared and uniaxial compression test and variable-angle shear test considering different bolt arrangements were carried out to obtain the mechanical parameters of the specimens.Results showed that the anisotropy of BGFR mainly depends on the relative volume content of each component material in a certain direction.Moreover,the mechanical parameters deduced from the theory of composite materials which consider the short fibre effect are shown to be in good agreement with those determined by laboratory experiments,and the variation rules maintained good consistency.Last,a case study of a real tunnel project is provided to highlight the effectiveness,validity and robustness of the developed E-P model in prediction of stresses and deformations.
基金the Water and Wastewater Company of East Azarbaijan Province for providing the funding for this research
文摘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.
基金supported by the National Natural Science Foundation of China(NSFC)(grant No.52074169,No.51704280)the China Postdoctoral Science Foundation(No.2023M732109)the Opening Foundation of Shandong Key Laboratory of Civil Engineering Disaster Prevention and Mitigation(No.CDPM2021FK02).
文摘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.
文摘Rock fracture mechanics and accurate characterization of rock fracture are crucial for understanding a variety of phenomena interested in geological engineering and geoscience.These phenomena range from very large-scale asymmetrical fault structures to the scale of engineering projects and laboratory-scale rock fracture tests.Comprehensive study can involve mechanical modeling,site or post-mortem investigations,and inspection on the point cloud of the source locations in the form of earthquake,microseismicity,or acoustic emission.This study presents a comprehensive data analysis on characterizing the forming of the asymmetrical damage zone around a laboratory mixed-mode rock fracture.We substantiate the presence of asymmetrical damage through qualitative analysis and demonstrate that measurement uncertainties cannot solely explain the observed asymmetry.The implications of this demonstration can be manifold.On a larger scale,it solidifies a mechanical model used for explaining the contribution of aseismic mechanisms to asymmetrical fault structures.On a laboratory scale,it exemplifies an alternative approach to understanding the observational difference between the source location and the in situ or post-mortem inspection on the rock fracture path.The mechanical model and the data analysis can be informative to the interpretations of other engineering practices as well,but may face different types of challenges.
基金Project(2021YFC2900600)supported by the Young Scientist Project of National Key Research and Development Program of ChinaProject(52074166)supported by the National Natural Science Foundation of China+1 种基金Projects(ZR2021YQ38,ZR2020QE121)supported by the Natural Science Foundation of Shandong Province,ChinaProject(2022KJ101)supported by the Science and Technology Support Plan for Youth Innovation of Colleges and Universities in Shandong Province,China。
文摘In practical engineering applications,rock mass are often found to be subjected to a triaxial stress state.Concurrently,defects like joints and fractures have a notable impact on the mechanical behavior of rock mass.Such defects are identified as crucial contributors to the failure and instability of the surrounding rock,subsequently impacting the engineering stability.The study aimed to investigate the impact of fracture geometry and confining pressure on the deformation,failure characteristics,and strength of specimens using sand powder 3D printing technology and conventional triaxial compression tests.The results indicate that the number of fractures present considerably influences the peak strength,axial peak strain and elastic modulus of the specimens.Confining pressure is an important factor affecting the failure pattern of the specimen,under which the specimen is more prone to shear failure,but the initiation,expansion and penetration processes of secondary cracks in different fracture specimens are different.This study confirmed the feasibility of using sand powder 3D printing specimens as soft rock analogs for triaxial compression research.The insights from this research are deemed essential for a deeper understanding of the mechanical behavior of fractured surrounding rocks when under triaxial stress state.
基金National Key R&D Program of China(2023YFC3009100,2023YFC3009102)National Natural Science Foundation of China(52304198)Open Fund of the National and Local Joint Engineering Research Center for Safe and Accurate Coal Mining(EC2021016).
文摘This study employs similar simulation testing and discrete element simulation coupling to analyze the failure and deformation processes of a model coal seam's roof.The caving area of the overburden rock is divided into three zones:the delamination fracture zone,broken fracture zone,and compaction zone.The caving and fracture zones'heights are approximately 110 m above the coal seam,with a maximum subsidence of 11 m.The delamination fracture zone's porosity range is between 0.2 and 0.3,while the remainder of the roof predominantly exhibits a porosity of less than 0.1.In addition,the numerical model's stress analysis revealed that the overburden rock's displacement zone forms an'arch-beam'structure starting from 160 m,with the maximum and minimum stress values decreasing as the distance of advancement increases.In the stress beam interval of the overburden rock,the maximum value changes periodically as the advancement distance increases.Based on a comparative analysis between observable data from on-site work and numerical simulation results,the stress data from the numerical simulation are essentially consistent with the actual results detected on-site,indicating the validity of the numerical simulation results.