The deep earth,deep sea,and deep space are the main parts of the national“three deep”strategy,which is in the forefront of the strategic deployment clearly defined in China’s 14th Five-Year Plan(2021-2025)and the L...The deep earth,deep sea,and deep space are the main parts of the national“three deep”strategy,which is in the forefront of the strategic deployment clearly defined in China’s 14th Five-Year Plan(2021-2025)and the Long-Range Objectives Through the Year 2035.It is important to reveal the evolutionary process and mechanism of deep tectonics to understand the earth’s past,present and future.The academic con-notation of Geology in Time has been given for the first time,which refers to the multi-field evolution response process of geological bodies at different time and spatial scales caused by geological processes inside and outside the Earth.Based on the deep in situ detection space and the unique geological envi-ronment of China Jinping Underground Laboratory,the scientific issue of the correlation mechanism and law between deep internal time-varying and shallow geological response is given attention.Innovative research and frontier exploration on deep underground in situ geo-information detection experiments for Geology in Time are designed to be carried out,which will have the potential to explore the driving force of Geology in Time,reveal essential laws of deep earth science,and explore innovative technologies in deep underground engineering.展开更多
The geometric characteristics of fractures within a rock mass can be inferred by the data sampling from boreholes or exposed surfaces.Recently,the universal elliptical disc(UED)model was developed to represent natural...The geometric characteristics of fractures within a rock mass can be inferred by the data sampling from boreholes or exposed surfaces.Recently,the universal elliptical disc(UED)model was developed to represent natural fractures,where the fracture is assumed to be an elliptical disc and the fracture orientation,rotation angle,length of the long axis and ratio of short-long axis lengths are considered as variables.This paper aims to estimate the fracture size-and azimuth-related parameters in the UED model based on the trace information from sampling windows.The stereological relationship between the trace length,size-and azimuth-related parameters of the UED model was established,and the formulae of the mean value and standard deviation of trace length were proposed.The proposed formulae were validated via the Monte Carlo simulations with less than 5%of error rate between the calculated and true values.With respect to the estimation of the size-and azimuth-related parameters using the trace length,an optimization method was developed based on the pre-assumed size and azimuth distribution forms.A hypothetical case study was designed to illustrate and verify the parameter estimation method,where three combinations of the sampling windows were used to estimate the parameters,and the results showed that the estimated values could agree well with the true values.Furthermore,a hypothetical three-dimensional(3D)elliptical fracture network was constructed,and the circular disc,non-UED and UED models were used to represent it.The simulated trace information from different models was compared,and the results clearly illustrated the superiority of the proposed UED model over the existing circular disc and non-UED models。展开更多
Rock joints are one of the vital discontinuities in a natural rock mass.How to accurately and conveniently determine joint normal stiffness is therefore significant in rock mechanics.Here,first,seven existing methods ...Rock joints are one of the vital discontinuities in a natural rock mass.How to accurately and conveniently determine joint normal stiffness is therefore significant in rock mechanics.Here,first,seven existing methods for determining joint normal stiffness were introduced and reviewed,among which MethodⅠ(the indirect measurement method),MethodⅡ(the direct determination method),MethodⅢ(the across-joint strain gauge measurement method)and MethodⅣ(the deformation measuring ring method)are via destructive uniaxial compression testing,while MethodⅤ(the best fitting method),MethodⅥ(the rapid evaluation method)and MethodⅦ(the effective modulus method)are through wave propagation principles and nondestructive ultrasonic testing.Subsequently,laboratory tests of intact and jointed sandstone specimens were conducted following the testing requirements and pro-cedures of those seven methods.A comparison among those methods was then performed.The results show that Method I,i.e.the benchmark method,is reliable and stable.MethodⅡhas a conceptual drawback,and its accuracy is acceptable at only very low stress levels.Relative errors in the results from MethodⅢare very large.With MethodⅣ,the testing results are sufficiently accurate despite the strict testing environment and complicated testing procedures.The results from MethodⅤare greatly unstable and significantly dependent on the natural frequency of the transducers.The joint normal stiffness determined with MethodⅥis stable and accurate,although data processing is complex.MethodⅦcould be adopted to determine the joint normal stiffness corresponding to the rock elastic deformation phase only.Consequently,it is suggested that MethodsⅠ,ⅣandⅥshould be adopted for the mea-surement of joint normal stiffness.The findings could be helpful in selecting an appropriate method to determine joint normal stiffness and,hence,to better solve discontinuous rock mass problems.展开更多
Since natural fractures are often non-equidimensional,the circular disc model still has great limitations.By contrast,the elliptical disc model is more applicable to representing natural fractures,especially for slend...Since natural fractures are often non-equidimensional,the circular disc model still has great limitations.By contrast,the elliptical disc model is more applicable to representing natural fractures,especially for slender ones.This paper developed a universal elliptical disc(UED)model by incorporating the center point,size,and azimuth of fractures as variables.Specifically,with respect to the azimuth of elliptical fractures in three-dimensional(3D)space,we proposed a paradigm to construct its probability density function(PDF)by coupling the orientation and rotation angle of long axis based on three coordinate transformations.To illustrate the construction process of the PDF of the fracture azimuth,we took the orientation following the Fisher distribution and the rotation angle following Von Mises distribution as an example.A rock slope is used to show the use of the developed UED model,and the 3D DFNs for the slope rock mass are generated by Monte Carlo simulation.In addition,the DFNs for the rock mass are also generated based on the existing circular disc model and non-universal elliptical disc model.The comparison results from the three models clearly illustrate the superiority of the UED model over the existing circular and non-universal elliptical disc models.展开更多
The paper proposes a new multiple-factor clustering method(NMFCM)with consideration of both box fractal dimension(BFD)and orientation of joints.This method assumes that the BFDs of different clusters were uneven,and c...The paper proposes a new multiple-factor clustering method(NMFCM)with consideration of both box fractal dimension(BFD)and orientation of joints.This method assumes that the BFDs of different clusters were uneven,and clustering was performed by redistributing the joints near the boundaries of clusters on a polar map to maximize an index for estimating the difference of the BFD(DBFD).Three main aspects were studied to develop the NMFCM.First,procedures of the NMFCM and reasonableness of assumptions were illustrated.Second,main factors affecting the NMFCM were investigated by numerical simulations with disk joint models.Finally,two different sections of a rock slope were studied to verify the practicability of the NMFCM.The results demonstrated that:(1)The NMFCM was practical and could effectively alleviate the problem of hard boundary during clustering;(2)The DBFD tended to increase after the improvement of clustering accuracy;(3)The improvement degree of the NMFCM clustering accuracy was mainly influenced by three parameters,namely,the number of clusters,number of redistributed joints,and total number of joints;and(4)The accuracy rate of clustering could be effectively improved by the NMFCM.展开更多
The deformation modulus of a rock mass is an important parameter to describe its mechanical behavior.In this study,an analytical method is developed to determine the deformation modulus of jointed rock masses,which co...The deformation modulus of a rock mass is an important parameter to describe its mechanical behavior.In this study,an analytical method is developed to determine the deformation modulus of jointed rock masses,which considers the mechanical properties of intact rocks and joints based on the superposition principle.Due to incorporating the variations in the orientations and sizes of joint sets,the proposed method is applicable to the rock mass with persistent and parallel joints as well as that with nonpersistent and nonparallel joints.In addition,an anisotropy index AIdmfor the deformation modulus is defined to quantitatively describe the anisotropy of rock masses.The range of AIdmis from 0 to 1,and the more anisotropic the rock mass is,the larger the value of AIdmwill be.To evaluate the proposed method,20 groups of numerical experiments are conducted with the universal distinct element code(UDEC).For each experimental group,the deformation modulus in 24 directions are obtained by UDEC(numerical value)and the proposed method(predicted value),and then the mean error rates are calculated.Note that the mean error rate is the mean value of the error rates of the deformation modulus in 24 directions,where for each direction,the error rate is equal to the ratio of numerical value minus predicted value to the numerical value.The results show that(i)for different experimental groups,the mean error rates vary between 5.06%and 22.03%;(ii)the error rates for the discrete fracture networks(DFNs)with two sets of joints are at the same level as those with one set of joints;and(iii)therefore,the proposed method for estimating the deformation modulus of jointed rock masses is valid.展开更多
Flexural toppling occurs when a series of layered rock masses bend towards their free face.It is important to evaluate the maximum bending degree and the requirement of supports of flexural toppling rock mass to preve...Flexural toppling occurs when a series of layered rock masses bend towards their free face.It is important to evaluate the maximum bending degree and the requirement of supports of flexural toppling rock mass to prevent rock mass cracking and even failure leading to a landslide.Based on the rock tensile strain-softening model,this study proposes a method for calculating the maximum curvature(C_(ppmax))of flexural toppling rock masses.By applying this method to calculate Cppmax of 9 types of rock masses with different hardness and rock layer thickness,some conclusions are drawn:(1)the internal key factors affecting C_(ppmax)are E^(⋆)(E^(⋆)=E_(ss)/E_(0),where E_(0)and E_(ss)are the mean deformation moduli of the rock before and after reaching its peak tensile strength,respectively),the strainεt corresponding to the tensile strength of rock,and the thickness(h)of rock layers;(2)hard rock layers are more likely to develop into block toppling than soft rock layers;and(3)thin rock layers are more likely to remain in flexural toppling state than thick rock layers.In addition,it is found that C_(ppmax)for flexural toppling rock masses composed of bedded rocks such as gneiss is related to the tensile direction.展开更多
This work used a modified direct shear apparatus, created newly by the authors, to explore effects of the gap between shear box halves and specimen size on the shear resistance of coarse-grained soil. The shear boxes ...This work used a modified direct shear apparatus, created newly by the authors, to explore effects of the gap between shear box halves and specimen size on the shear resistance of coarse-grained soil. The shear boxes of this apparatus were assembled from a series of steel structures capable of superimposition and nesting. Such characteristics facilitated variation of specimen size in both diameter and height. The new device can also maintain a constant gap during shearing. We performed a series of gap-effect and size-effect tests for two uniformly graded, coarse^grained soil samples. The test results showed that both the gap space and specimen size bad significant influences on shear resistance of the coarse-grained soil. Further, analysis of variations in shear strength indices led to a reasonable gap dimension and specimen size of the two soil samples.展开更多
This paper studied the rock dynamic fracture propagation under impact loads elaborately with a determination method proposed to calculate crack propagation dynamic stress intensity factor(DSIF).By utilizing the split-...This paper studied the rock dynamic fracture propagation under impact loads elaborately with a determination method proposed to calculate crack propagation dynamic stress intensity factor(DSIF).By utilizing the split-Hopkinson pressure bar,the impact experiments with an improved single cleavage semi-circle(ISCSC)specimen were conducted to illuminate the dynamic crack propagation behaviour.Meanwhile,the fracture characteristics and crack propagation velocity were obtained by the crack propagation gauges.Coordinating experiments with a numerical approach,the crack propagation dynamic stress intensity factors were calculated by an experimental—numerical method with fractal theory.Then,a finite difference model was developed based on the tensile fracture softening damage criterion.With the analysis of numerical and experimental results,the crack propagation behaviour and mechanism of crack arrest were discussed sophisticatedly.The results demonstrate that the novel ISCSC specimen shows a definite advantage in determining crack propagation and arrest DSIF.Additionally,the crack arrest DSIF is larger than the average propagation DSIF with a sharp increase.Meanwhile,the numerical simulation results which agree well with the actual crack propagation illustrate that the crack arrest should be dominated by the compressive stress perpendicular to the crack path,and there were several arrest pauses existing in the transitory crack arrest process.展开更多
基金supported by the National Natural Science Foundation of China(Nos.52125402 and 52174084)the Natural Science Foundation of Sichuan Province of China(No.2022NSFSC0005).
文摘The deep earth,deep sea,and deep space are the main parts of the national“three deep”strategy,which is in the forefront of the strategic deployment clearly defined in China’s 14th Five-Year Plan(2021-2025)and the Long-Range Objectives Through the Year 2035.It is important to reveal the evolutionary process and mechanism of deep tectonics to understand the earth’s past,present and future.The academic con-notation of Geology in Time has been given for the first time,which refers to the multi-field evolution response process of geological bodies at different time and spatial scales caused by geological processes inside and outside the Earth.Based on the deep in situ detection space and the unique geological envi-ronment of China Jinping Underground Laboratory,the scientific issue of the correlation mechanism and law between deep internal time-varying and shallow geological response is given attention.Innovative research and frontier exploration on deep underground in situ geo-information detection experiments for Geology in Time are designed to be carried out,which will have the potential to explore the driving force of Geology in Time,reveal essential laws of deep earth science,and explore innovative technologies in deep underground engineering.
基金funded by National Natural Science Foundation of China(Grant No.41972264)Zhejiang Provincial Natural Science Foundation of China(Grant No.LR22E080002)the Observation and Research Station of Geohazards in Zhejiang,Ministry of Natural Resources,China(Grant No.ZJDZGCZ-2021).
文摘The geometric characteristics of fractures within a rock mass can be inferred by the data sampling from boreholes or exposed surfaces.Recently,the universal elliptical disc(UED)model was developed to represent natural fractures,where the fracture is assumed to be an elliptical disc and the fracture orientation,rotation angle,length of the long axis and ratio of short-long axis lengths are considered as variables.This paper aims to estimate the fracture size-and azimuth-related parameters in the UED model based on the trace information from sampling windows.The stereological relationship between the trace length,size-and azimuth-related parameters of the UED model was established,and the formulae of the mean value and standard deviation of trace length were proposed.The proposed formulae were validated via the Monte Carlo simulations with less than 5%of error rate between the calculated and true values.With respect to the estimation of the size-and azimuth-related parameters using the trace length,an optimization method was developed based on the pre-assumed size and azimuth distribution forms.A hypothetical case study was designed to illustrate and verify the parameter estimation method,where three combinations of the sampling windows were used to estimate the parameters,and the results showed that the estimated values could agree well with the true values.Furthermore,a hypothetical three-dimensional(3D)elliptical fracture network was constructed,and the circular disc,non-UED and UED models were used to represent it.The simulated trace information from different models was compared,and the results clearly illustrated the superiority of the proposed UED model over the existing circular disc and non-UED models。
基金supported by the Shenzhen Fundamental Research Program(Grant No.JCYJ20220818095605012)the National Natural Science Foundation of China(Grant No.51909026)the Fund of Guangdong Provincial Key Laboratory of Deep Earth Sciences and Geothermal Energy Exploitation and Utilization(Grant No.2020-08).
文摘Rock joints are one of the vital discontinuities in a natural rock mass.How to accurately and conveniently determine joint normal stiffness is therefore significant in rock mechanics.Here,first,seven existing methods for determining joint normal stiffness were introduced and reviewed,among which MethodⅠ(the indirect measurement method),MethodⅡ(the direct determination method),MethodⅢ(the across-joint strain gauge measurement method)and MethodⅣ(the deformation measuring ring method)are via destructive uniaxial compression testing,while MethodⅤ(the best fitting method),MethodⅥ(the rapid evaluation method)and MethodⅦ(the effective modulus method)are through wave propagation principles and nondestructive ultrasonic testing.Subsequently,laboratory tests of intact and jointed sandstone specimens were conducted following the testing requirements and pro-cedures of those seven methods.A comparison among those methods was then performed.The results show that Method I,i.e.the benchmark method,is reliable and stable.MethodⅡhas a conceptual drawback,and its accuracy is acceptable at only very low stress levels.Relative errors in the results from MethodⅢare very large.With MethodⅣ,the testing results are sufficiently accurate despite the strict testing environment and complicated testing procedures.The results from MethodⅤare greatly unstable and significantly dependent on the natural frequency of the transducers.The joint normal stiffness determined with MethodⅥis stable and accurate,although data processing is complex.MethodⅦcould be adopted to determine the joint normal stiffness corresponding to the rock elastic deformation phase only.Consequently,it is suggested that MethodsⅠ,ⅣandⅥshould be adopted for the mea-surement of joint normal stiffness.The findings could be helpful in selecting an appropriate method to determine joint normal stiffness and,hence,to better solve discontinuous rock mass problems.
基金funded by the National Natural Science s of China(No.41972264)the Key R&D Project of Zhejiang Province(No.2021C03159)the Field Scientific Observation&Research Station of Geological Hazard in Zhejiang,Ministry of Natural Resources,China(No.ZJDZGCZ-2021)。
文摘Since natural fractures are often non-equidimensional,the circular disc model still has great limitations.By contrast,the elliptical disc model is more applicable to representing natural fractures,especially for slender ones.This paper developed a universal elliptical disc(UED)model by incorporating the center point,size,and azimuth of fractures as variables.Specifically,with respect to the azimuth of elliptical fractures in three-dimensional(3D)space,we proposed a paradigm to construct its probability density function(PDF)by coupling the orientation and rotation angle of long axis based on three coordinate transformations.To illustrate the construction process of the PDF of the fracture azimuth,we took the orientation following the Fisher distribution and the rotation angle following Von Mises distribution as an example.A rock slope is used to show the use of the developed UED model,and the 3D DFNs for the slope rock mass are generated by Monte Carlo simulation.In addition,the DFNs for the rock mass are also generated based on the existing circular disc model and non-universal elliptical disc model.The comparison results from the three models clearly illustrate the superiority of the UED model over the existing circular and non-universal elliptical disc models.
基金funded by the National Natural Science Foundation of China(Grant Nos.41972264 and 52078093)Liaoning Revitalization Talents Program,China(Grant No.XLYC1905015)。
文摘The paper proposes a new multiple-factor clustering method(NMFCM)with consideration of both box fractal dimension(BFD)and orientation of joints.This method assumes that the BFDs of different clusters were uneven,and clustering was performed by redistributing the joints near the boundaries of clusters on a polar map to maximize an index for estimating the difference of the BFD(DBFD).Three main aspects were studied to develop the NMFCM.First,procedures of the NMFCM and reasonableness of assumptions were illustrated.Second,main factors affecting the NMFCM were investigated by numerical simulations with disk joint models.Finally,two different sections of a rock slope were studied to verify the practicability of the NMFCM.The results demonstrated that:(1)The NMFCM was practical and could effectively alleviate the problem of hard boundary during clustering;(2)The DBFD tended to increase after the improvement of clustering accuracy;(3)The improvement degree of the NMFCM clustering accuracy was mainly influenced by three parameters,namely,the number of clusters,number of redistributed joints,and total number of joints;and(4)The accuracy rate of clustering could be effectively improved by the NMFCM.
基金funded by the National Key R&D Program of China(Grant Nos.2017YFE0119500 and 2018YFC1505005)the National Natural Science Foundation of China(Grant No.41972264)。
文摘The deformation modulus of a rock mass is an important parameter to describe its mechanical behavior.In this study,an analytical method is developed to determine the deformation modulus of jointed rock masses,which considers the mechanical properties of intact rocks and joints based on the superposition principle.Due to incorporating the variations in the orientations and sizes of joint sets,the proposed method is applicable to the rock mass with persistent and parallel joints as well as that with nonpersistent and nonparallel joints.In addition,an anisotropy index AIdmfor the deformation modulus is defined to quantitatively describe the anisotropy of rock masses.The range of AIdmis from 0 to 1,and the more anisotropic the rock mass is,the larger the value of AIdmwill be.To evaluate the proposed method,20 groups of numerical experiments are conducted with the universal distinct element code(UDEC).For each experimental group,the deformation modulus in 24 directions are obtained by UDEC(numerical value)and the proposed method(predicted value),and then the mean error rates are calculated.Note that the mean error rate is the mean value of the error rates of the deformation modulus in 24 directions,where for each direction,the error rate is equal to the ratio of numerical value minus predicted value to the numerical value.The results show that(i)for different experimental groups,the mean error rates vary between 5.06%and 22.03%;(ii)the error rates for the discrete fracture networks(DFNs)with two sets of joints are at the same level as those with one set of joints;and(iii)therefore,the proposed method for estimating the deformation modulus of jointed rock masses is valid.
基金funded by the National Natural Science Foundation of China(No.41972264)Zhejiang Provincial Natural Science Foundation of China(No.LR22E080002)the Key R&D Project of Zhejiang Province(No.2021C03159).
文摘Flexural toppling occurs when a series of layered rock masses bend towards their free face.It is important to evaluate the maximum bending degree and the requirement of supports of flexural toppling rock mass to prevent rock mass cracking and even failure leading to a landslide.Based on the rock tensile strain-softening model,this study proposes a method for calculating the maximum curvature(C_(ppmax))of flexural toppling rock masses.By applying this method to calculate Cppmax of 9 types of rock masses with different hardness and rock layer thickness,some conclusions are drawn:(1)the internal key factors affecting C_(ppmax)are E^(⋆)(E^(⋆)=E_(ss)/E_(0),where E_(0)and E_(ss)are the mean deformation moduli of the rock before and after reaching its peak tensile strength,respectively),the strainεt corresponding to the tensile strength of rock,and the thickness(h)of rock layers;(2)hard rock layers are more likely to develop into block toppling than soft rock layers;and(3)thin rock layers are more likely to remain in flexural toppling state than thick rock layers.In addition,it is found that C_(ppmax)for flexural toppling rock masses composed of bedded rocks such as gneiss is related to the tensile direction.
文摘This work used a modified direct shear apparatus, created newly by the authors, to explore effects of the gap between shear box halves and specimen size on the shear resistance of coarse-grained soil. The shear boxes of this apparatus were assembled from a series of steel structures capable of superimposition and nesting. Such characteristics facilitated variation of specimen size in both diameter and height. The new device can also maintain a constant gap during shearing. We performed a series of gap-effect and size-effect tests for two uniformly graded, coarse^grained soil samples. The test results showed that both the gap space and specimen size bad significant influences on shear resistance of the coarse-grained soil. Further, analysis of variations in shear strength indices led to a reasonable gap dimension and specimen size of the two soil samples.
基金the National Natural Science Foundation of China(Grant Nurmbers 11702181,11672194)the Sichuan Science and Technology Program(Grant Numbers 2019YFG0047).
文摘This paper studied the rock dynamic fracture propagation under impact loads elaborately with a determination method proposed to calculate crack propagation dynamic stress intensity factor(DSIF).By utilizing the split-Hopkinson pressure bar,the impact experiments with an improved single cleavage semi-circle(ISCSC)specimen were conducted to illuminate the dynamic crack propagation behaviour.Meanwhile,the fracture characteristics and crack propagation velocity were obtained by the crack propagation gauges.Coordinating experiments with a numerical approach,the crack propagation dynamic stress intensity factors were calculated by an experimental—numerical method with fractal theory.Then,a finite difference model was developed based on the tensile fracture softening damage criterion.With the analysis of numerical and experimental results,the crack propagation behaviour and mechanism of crack arrest were discussed sophisticatedly.The results demonstrate that the novel ISCSC specimen shows a definite advantage in determining crack propagation and arrest DSIF.Additionally,the crack arrest DSIF is larger than the average propagation DSIF with a sharp increase.Meanwhile,the numerical simulation results which agree well with the actual crack propagation illustrate that the crack arrest should be dominated by the compressive stress perpendicular to the crack path,and there were several arrest pauses existing in the transitory crack arrest process.