Since its introduction,discontinuous deformation analysis(DDA)has been widely used in different areas of rock mechanics.By dividing large blocks into subblocks and introducing artificial joints,DDA can be applied to r...Since its introduction,discontinuous deformation analysis(DDA)has been widely used in different areas of rock mechanics.By dividing large blocks into subblocks and introducing artificial joints,DDA can be applied to rock fracture simulation.However,parameter calibration,a fundamental issue in discontinuum methods,has not received enough attention in DDA.In this study,the parameter calibration of DDA for intact rock is carefully studied.To this end,a subblock DDA with Voronoi tessellation is presented first.Then,a modified contact constitutive law is introduced,in which the tensile and shear meso-strengths are modified to be independent of the bond lengths.This improvement can prevent the unjustified preferential failure of short edges.A method for imposing confining pressure is also introduced.Thereafter,sensitivity analysis is performed to investigate the influence of the calculated parameters and meso-parameters on the mechanical properties of modeled rock.Based on the sensitivity analysis,a unified calibration procedure is suggested for both cases with and without confining pressure.Finally,the calibration procedure is applied to two examples,including a biaxial compression test.The results show that the proposed Voronoi-based DDA can simulate rock fracture with and without confining pressure very well after careful parameter calibration.展开更多
Deep-seated rock fractures(referred to as DSRF hereafter)in valley slopes are uncommon geological phenomena that challenge our previous understanding of slope unloading processes.These fractures weaken the strength an...Deep-seated rock fractures(referred to as DSRF hereafter)in valley slopes are uncommon geological phenomena that challenge our previous understanding of slope unloading processes.These fractures weaken the strength and integrity of the rock mass,potentially forming unstable block boundaries with significant volume,thereby affecting the stability of slopes,chambers,and dam abutments.DSRF has emerged as a critical environmental and engineering geological issue that hinders large-scale projects in deep canyon areas.Despite the attention and practical treatment given to DSRF in engineering practice,theoretical research on this topic still lags behind the demands of engineering applications.To garner widespread attention and promote the resolution of DSRF-related problems,this review aims to redefine DSRF through comprehensive data collection and analysis,engineering geological analogies,and field investigations,and provide a summary and analysis of the research progress on DSRF,along with future research directions.The study defines DSRF as the intermittent tension cracks or relaxation zones within a slightly weathered or fresh,and intact or relatively intact rock mass distributed below the surface unloading zones of a deep canyon slope,and should be distinguished from"loose rock mass"and"deep-seated gravitational slope deformations".The article provides an overview of the development and distribution,rupture characteristics,and genesis mechanism of DSRF.It proposes that DSRF is formed based on the fluvial deviation-undercutting evolution mode,wherein the energy accumulated in the rock mass is violently released when the river further down cuts the slope after the rock mass has undergone cyclical loadingunloading.However,further research is necessary to establish a comprehensive database for DSRF,refine exploration techniques,understand evolutionary processes,develop engineering evaluation methods,and predict the distribution of DSRF.展开更多
Underground space creation and energy extraction, which induce unloading on rock fractures, commonly occur in various rock engineering projects, and rock engineering projects are subjected to high temperatures with in...Underground space creation and energy extraction, which induce unloading on rock fractures, commonly occur in various rock engineering projects, and rock engineering projects are subjected to high temperatures with increasing depth. Fluid flow behavior of rock fractures is a critical issue in many subsurface rock engineering projects. Previous studies have extensively considered permeability evolution in rock fractures under loading phase, whereas changes in fracture permeability under unloading phase have not been fully understood. To examine the unloading-induced changes in fracture permeability under different temperatures, we performed water flow-through tests on fractured rock samples subjected to decreasing confining pressures and different temperatures. The experimental results show that the permeability of fracture increases with unloading of confining pressure but decreases with loading-unloading cycles. Temperature may affect fracture permeability when it is higher than a certain threshold. An empirical model of fracture hydraulic aperture including two material parameters of initial normal stiffness and maximum normal closure can well describe the permeability changes in rough rock fracture subjected to loading-unloading cycles and heating. A coupled thermo-mechanical model considering asperity damage is finally used to understand the influences of stress paths and temperatures on fracture permeability.展开更多
A fractal model governing saw-tooth fractures was first introduced to replicate sandstone samples containing an inclined 3D penetrating rough fracture surface with various joint roughness coefficients(JRC).In conventi...A fractal model governing saw-tooth fractures was first introduced to replicate sandstone samples containing an inclined 3D penetrating rough fracture surface with various joint roughness coefficients(JRC).In conventional triaxial compression,the peak strength for fractured samples increased with both confining pressure and JRC.During the unloading confining pressure process,the normal stress of fractures declined but the shear stress increased,resulting in shear sliding of fractures.The shear displacement of fractures exponentially increased,and the positive normal displacement decreased gradually to negative values under coupling effects of shear contraction caused by normal stress and shear dilation due to climbing effects of fractures.Transition from quasi-static to dynamic sliding of the fractures was identified.The sliding resistance duration increased with confining pressure but decreased with JRC.After prepeak unloading,the fracture surfaces presented a more significant surface wear response and JRC values decreased by 1.70%–59.20%due to more remarkable asperity degradation compared with those after conventional triaxial compression.The theoretical model for shear strength of fractures was established through improving the Ladanyi&Archambault model by introducing the relations between normal stress and surface wear ratios of fractures,which agreed well with the experimental results.展开更多
This study experimentally analyzes the nonlinear flow characteristics and channelization of fluid through rough-walled fractures during the shear process using a shear-flow-visualization apparatus.A series of fluid fl...This study experimentally analyzes the nonlinear flow characteristics and channelization of fluid through rough-walled fractures during the shear process using a shear-flow-visualization apparatus.A series of fluid flow and visualization tests is performed on four transparent fracture specimens with various shear displacements of 1 mm,3 mm,5 mm,7 mm and 10 mm under a normal stress of 0.5 MPa.Four granite fractures with different roughnesses are selected and quantified using variogram fractal dimensions.The obtained results show that the critical Reynolds number tends to increase with increasing shear displacement but decrease with increasing roughness of fracture surface.The flow paths are more tortuous at the beginning of shear because of the wide distribution of small contact spots.As the shear displacement continues to increase,preferential flow paths are more distinctly observed due to the decrease in the number of contact spots caused by shear dilation;yet the area of single contacts in-creases.Based on the experimental results,an empirical mathematical equation is proposed to quantify the critical Reynolds number using the contact area ratio and fractal dimension.展开更多
Three-dimensional rock fracture induced by blasting is a highly complex problem and has received considerable attention in geotechnical engineering.The material point method is firstly applied to treat this challengin...Three-dimensional rock fracture induced by blasting is a highly complex problem and has received considerable attention in geotechnical engineering.The material point method is firstly applied to treat this challenging task.Some inherent weaknesses can be overcome by coupling the generalized interpolation material point(GIMP)and the convected particle domain interpolation technique(CPDI).For the media in the borehole,unchanged GIMP-type particles are used to guarantee a homogenous blast pressure.CPDITetrahedron type particles are employed to avoid the fake numerical fracture near the borehole for the rock material.A blasting experiment using three-dimensional single-borehole rock was simulated to examine the applicability of the coupled model under realistic loading and boundary conditions.A good agreement was achieved between the simulation and experimental results.Moreover,the mechanism of three-dimensional rock fracture was analyzed.It was concluded that rock particle size and material parameters play an important role in rock damage.The reflected tensile waves cause severe damage in the lower part of the model.Rayleigh waves occur on the top face of the rock model to induce a hoop failure band.展开更多
A convenient approach was proposed by which to evaluate and monitor the permeability of a rock fracture by verifying the quantitative correlation between the electrical resistivity and permeability at laboratory scale...A convenient approach was proposed by which to evaluate and monitor the permeability of a rock fracture by verifying the quantitative correlation between the electrical resistivity and permeability at laboratory scale.For this purpose,an electrical resistivity measurement system was applied to the laboratory experiments using artificial cells with the shape of a single rock fracture.Sixty experiments were conducted using rock fractures according to the geometry,aperture sizes,wavelengths,and roughness amplitudes.The overall negative relationship between the normalized electrical resistivity values and the aperture sizes directly linked with the permeability,was well fitted by the power-law function with a large determination coefficient(≈0.86).The effects of wavelength and roughness amplitude of the rock fracture on the electrical resistivity were also analyzed.Results showed that the electrical resistivity was slightly increased with decreasing wavelength and increasing roughness amplitude.An empirical model for evaluating the permeability of a rock fracture was proposed based on the experimental data.In the field,if the electrical resistivity of pore groundwater could be measured in advance,this empirical model could be applied effectively for simple,quick monitoring of the fracture permeability.Although uncertainty may be associated with the permeability estimation due to the limited control parameters considered in this research,this electrical resistivity approach could be helpful to monitor the rock permeability in deep underground facilities such as those used for radioactive waste repositories or forms of energy storage.展开更多
Accurate determination of rock mass properties is a critical part of open-pit mine planning activities to enable more prescriptive blast designs to achieve improved loading and hauling and downstream process efficienc...Accurate determination of rock mass properties is a critical part of open-pit mine planning activities to enable more prescriptive blast designs to achieve improved loading and hauling and downstream process efficiency. Better and more accurate blast practice that delivers enhanced outcomes(better fragmentation, improved diggability, less dilution, etc.) is a critical and fundamental element of being able to achieve an effective Mine-Mill approach at a mining operation. Based on previous work, it has been demonstrated that an accurate representation of the rock mass properties can be obtained from the analysis of variations in blasthole drill performance as derived from measurement while drilling(MWD) systems when using tricone bits. This paper further investigates how monitored rate of penetration,pulldown force, rotary torque, rotary speed and bailing air pressure responses can be used to determine the presence of open and partially open fractures having varying dip angles. Based on a correlation of geophysically measured fracture logs and monitored drill performance variables in the same blastholes, the results show that the latter responses can accurately determine open versus closed fractures. The results also identified that variations in rate of penetration and rotary torque show the most sensitivity in the presence of open fractures that intersect a vertical blasthole at near orthogonal angles.展开更多
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.展开更多
In order to improve our understanding of rock fracture and fault instability driven by high-pressure fluid sources, the authors carried out rock fracture tests using granite under a confining pressure of 80 MPa with f...In order to improve our understanding of rock fracture and fault instability driven by high-pressure fluid sources, the authors carried out rock fracture tests using granite under a confining pressure of 80 MPa with fluid injection in the laboratory. Furthermore, we tested a number of numerical models using the FLAC^(3D) modeling software to find the best model to represent the experimental results. The high-speed multichannel acoustic emission(AE) waveform recording system used in this study made it possible to examine the total fracture process through detailed monitoring of AE hypocenters and seismic velocity.The experimental results show that injecting high-pressure oil into the rock sample can induce AE activity at very low stress levels and can dramatically reduce the strength of the rock. The results of the numerical simulations show that major experimental results, including the strength, the temporal and spatial patterns of the AE events, and the role of the fluid can be represented fairly well by a model involving(1) randomly distributed defect elements to model pre-existing cracks,(2) random modification of rock properties to represent inhomogeneity introduced by different mineral grains, and(3)macroscopic inhomogeneity. Our study, which incorporates laboratory experiments and numerical simulations, indicates that such an approach is helpful in finding a better model not only for simulating experimental results but also for upscaling purposes.展开更多
To identify refl ector fractures near borehole by using dipole-source refl ected-shearwave logging, we need to understand the relation between the amplitude of the refl ected shear wave and the source radiation, boreh...To identify refl ector fractures near borehole by using dipole-source refl ected-shearwave logging, we need to understand the relation between the amplitude of the refl ected shear wave and the source radiation, borehole conditions, and attenuation owing to the surrounding formations. To assess the effect of these factors on the amplitude of the refl ected waves, we first studied the radiation performance and radiation direction of the dipole source in fast, medium, and slow formations by using the asymptotic solution in the far fi eld of the borehole. Then, the relation between the fracture parameters, and the refl ected-shear-wave amplitude as well as the ratio of the refl ected-shear-wave amplitude to the direct-wave amplitude (relative amplitude, RA) was evaluated by the three-dimensional fi nite-difference (3D FDTD) method. Finally, the fracture detection capability of the dipole reflected-shear-wave logging tool in different formations was analyzed by using the RA. The results suggest that the radiation amplitude of the SH-wave in the slow formation is weaker than those in the fast and medium formations, and the amplitude of the refl ected shear wave is lower. However, the RA in the slow formation is close to or even greater than in the fast and medium formations, which means that dipole-source shear-wave logging has the same or even better fracture detection capability in the slow formation as in the fast and medium formations. In addition, when RA is small, there is a good correlation between the RA and the various fracture parameters in the different types of formation, which can be used in determining the lower limit of the fracture parameters identifi ed by refl ection logging.展开更多
This paper presents a novel integrated method for interactive characterization of fracture spacing in rock tunnel sections.The main procedure includes four steps:(1)Automatic extraction of fracture traces,(2)digitizat...This paper presents a novel integrated method for interactive characterization of fracture spacing in rock tunnel sections.The main procedure includes four steps:(1)Automatic extraction of fracture traces,(2)digitization of trace maps,(3)disconnection and grouping of traces,and(4)interactive measurement of fracture set spacing,total spacing,and surface rock quality designation(S-RQD)value.To evaluate the performance of the proposed method,sample images were obtained by employing a photogrammetrybased scheme in tunnel faces.Experiments were then conducted to determine the optimal parameter values(i.e.distance threshold,angle threshold,and number of fracture trace grouping)for characterizing rock fracture spacing.By applying the identified optimal parameters involved in the model,the proposed method could lead to excellent qualitative results to a new tunnel face.To perform a quantitative analysis,three methods(i.e.field,straightening,and the proposed method)were employed in the same study and comparisons were made.The proposed method agrees well with the field measurement in terms of the maximum and average values of measured spacing distribution.Overall,the proposed method has reasonably good accuracy and interactive advantage for estimating the ultimate fracture spacing and S-RQD.It can be a possible extension of existing methods for fracture spacing characterization for two-dimensional(2D)rock tunnel faces.展开更多
Numerous intersected rock fractures constitute the fracture network in enhanced geothermal systems.The complicated convective heat transfer behavior in intersected fractures is critical to the heat recovery in fractur...Numerous intersected rock fractures constitute the fracture network in enhanced geothermal systems.The complicated convective heat transfer behavior in intersected fractures is critical to the heat recovery in fractured geothermal reservoirs.A series of three-dimensional intersected fracture models is constructed to perform the flow-through heat transfer simulations.The geometry effects of dead-end fractures(DEFs)on the heat transfer are evaluated in terms of intersected angles,apertures,lengths,and the connectivity.The results indicate that annular streamlines appear in the rough DEF and cause an ellipse distribution of the cold front.Compared to plate DEFs,the fluid flow in the rough DEF enhances the heat transfer.Both the increment of outlet water temperatureΔToutand the ratio of heat production Qrpresent the largest at the intersected angle of 90°while decline with the decrease of the intersected angle between the main flow fracture(MFF)and the DEFs.The extension of the length of intersected DEFs is beneficial to heat production while enhancing its aperture is not needed.Solely increasing the number of intersected DEFs induces a little increase of heat extraction,and more significant heat production can be obtained through connecting these DEFs with the MFF forming the flow network.展开更多
Anti-plane punch-through shear test and anti-planefour-point bending test are used to study the crack initiation and propagation under anti-plane shear (Mode Ⅲ) loading. The tensile and shear stresses at the crack ti...Anti-plane punch-through shear test and anti-planefour-point bending test are used to study the crack initiation and propagation under anti-plane shear (Mode Ⅲ) loading. The tensile and shear stresses at the crack tip are calculated by finite element method. The results show that under Mode Ⅲ loading the maximum principal stress σ1 at crack tip is smaller or a little larger than the maximum shear stress τmax. Since the tensile strength of brittle rock is much lower than its shear strength, σ1 is easy to reach its critical value before τmax reaches its critical value and thus results in Mode I fracture. The fracture trajectory is helicoid and the normal direction of tangential plane with the fractured helicoid is along the predicted direction of the maximum principal stress at the notch tip. It is further proved that Mode Ⅰ instead of Mode Ⅲ fracture occurs in brittle rock under Mode Ⅲ loading. The fracture mode depending on the fracture mechanism must be distinguished from the loading form.展开更多
The joint roughness coefficient (JRC) is an important geometric parameter in rock mechanics, seismo-geology, engineering geology, and structural geology. In recent years, the fractol geometry has been widely applied i...The joint roughness coefficient (JRC) is an important geometric parameter in rock mechanics, seismo-geology, engineering geology, and structural geology. In recent years, the fractol geometry has been widely applied in researching JRC In this paper, a simpler fractal model for estimating JRC was developed, whick can be applied in simulating rock discontinuity profiles, aud the rougher the surface, the greater the fractal dimension D The epirical equation of JRC and D was established according to the model’ The JRC values of six joint profiles were estimated according to Eq (10),which are ratker approximate to the typical JRC values; so the method may be applied to estimate the roughness of the fractured surfaces of natural rocks For a give profile the estimates of D and JRC can be obtained from Eqs (5) and (l0) with only one parameter e’ (the average angle of the profile fluctuation) being needed.展开更多
The application of the non-explosive expansion material (NEEM) is widely used as the controlled fracture method in quarry mining,especially in hard rocks.The pressure of NEEM is an important parameter in causing rock ...The application of the non-explosive expansion material (NEEM) is widely used as the controlled fracture method in quarry mining,especially in hard rocks.The pressure of NEEM is an important parameter in causing rock fracture.An empirical model based on hole spacing was developed to determine the pressure of NEEM in the rock fracture process.Primarily,the empirical model was developed by the mathematical method,utilizing dimensional analysis.Then,the Phase2 code,which is based on the finite element method,was utilized to predict crack growth in rocks.The results of numerical analysis show slight deviations from the empirical model.Hence,the polynomial regression analysis was used to modify the model.Finally,the modified model shows a good agreement with the results gained from numerical modeling.展开更多
Through seepage tests under different loading and unloading confining pressures and different hydraulic gradients,the authors studied the effects of stress states and stress history on fracture permeability evolution ...Through seepage tests under different loading and unloading confining pressures and different hydraulic gradients,the authors studied the effects of stress states and stress history on fracture permeability evolution for single granite fracture and sandstone fracture. The results show that there exists a linear relationship between the seepage discharge and osmotic pressure in sandstone fissure under each level of confining pressure. With the increasing in the confining pressure,the permeability of the fracture decreases,but the decreasing rate is changeing. During the unloading process,the fracture seepage velocity cannot be fully recovered to the size of the loading process. Therefore,in the unloading process of the confining pressure,the recovery of fracture permeability shows obvious hysteresis effects. The flow rate of the fracture remains unchanged during five cycles of loading and unloading processes of the confining pressure. In each cycle,the evolution character of the flow rate with the confining pressure remains unchanged. These experiments show that the seepage characteristics of sandstone and granite fractures are not the same under the same stress state.展开更多
It is observed that the parameter of seismic inhomogeneous degree (GL value) calculated from the earthquake catalog shows obvious abnormal changes prior to strong earthquakes, indicating the state change of local seis...It is observed that the parameter of seismic inhomogeneous degree (GL value) calculated from the earthquake catalog shows obvious abnormal changes prior to strong earthquakes, indicating the state change of local seismic activity. This paper focuses on the mechanism for the abnormal changes of the GL values based on the sequences of acoustic emission for three types of rock samples containing macro-asperity fracture; compressional en-echelon fracture and model-III shear fracture. The results show that for the three types of rock samples, there are continuous abnormal changes of GL value (>1) just before the non-elastic deformation occurs or during the process of nucleation prior to the instability. Based on the experimental results, it seems that the process of creep sliding and resistance-uniformization along fault zone is the possible mechanism for the abnormal changes of GL value before rock fractures.展开更多
Fracture and ground vibration of rock subjected to different decoupling decked charges are investigated based on the numerical simulation.The dynamic pressure value is studied,which demonstrates that simulation of fra...Fracture and ground vibration of rock subjected to different decoupling decked charges are investigated based on the numerical simulation.The dynamic pressure value is studied,which demonstrates that simulation of fracture zone is feasible.Attenuation index of dynamic pressure is 2.06,2.05 and 1.93 for air,water and sand intervals respectively.The small attenuation of sand interval results in the large ground vibration.The predicted vertical vibration waveform and peak particle velocities(PPV) in far-field are in agreement with the monitoring results.The results show that the air and water decked charges can improve the effect of rock fracture in near-field and reduce ground vibration in far-field.展开更多
A true triaxial apparatus which is composed of three units was presented.The apparatus allows for investigations on deformation and seepage behaviors of a single rock fracture subjected to lateral stress and normal st...A true triaxial apparatus which is composed of three units was presented.The apparatus allows for investigations on deformation and seepage behaviors of a single rock fracture subjected to lateral stress and normal stress.The first unit has three jacks which can apply loads independently in three orthogonal directions.The second unit is used to supply water inflow,control seepage pressure and measure flow velocity in real time.The third unit is for measuring the normal deformation of rock fractures.Some tests for investigating the normal deformation and seepage behaviors of rock fractures subjected to normal and lateral loads on hard granite specimens with an artificial persistent fracture,were introduced.The results show that both the normal deformation and the hydraulic conductivity are influenced not only by the normal stress but also by the lateral stress.It is also shown that the aperture and the hydraulic conductivity decrease with the increasing normal stress but increase with the increasing lateral stress and both the aperture and the hydraulic conductivity obey exponential relationships with the normal stress and the lateral stress.展开更多
基金The authors would like to thank the National Natural Science Foundation of China(Grant Nos.51879184 and 52079091)for funding this work.
文摘Since its introduction,discontinuous deformation analysis(DDA)has been widely used in different areas of rock mechanics.By dividing large blocks into subblocks and introducing artificial joints,DDA can be applied to rock fracture simulation.However,parameter calibration,a fundamental issue in discontinuum methods,has not received enough attention in DDA.In this study,the parameter calibration of DDA for intact rock is carefully studied.To this end,a subblock DDA with Voronoi tessellation is presented first.Then,a modified contact constitutive law is introduced,in which the tensile and shear meso-strengths are modified to be independent of the bond lengths.This improvement can prevent the unjustified preferential failure of short edges.A method for imposing confining pressure is also introduced.Thereafter,sensitivity analysis is performed to investigate the influence of the calculated parameters and meso-parameters on the mechanical properties of modeled rock.Based on the sensitivity analysis,a unified calibration procedure is suggested for both cases with and without confining pressure.Finally,the calibration procedure is applied to two examples,including a biaxial compression test.The results show that the proposed Voronoi-based DDA can simulate rock fracture with and without confining pressure very well after careful parameter calibration.
基金supported by the National Natural Science Foundation of China(Grant Nos.41272333)the National Key Research and Development Program of China(Grant Nos.2011CB013501)。
文摘Deep-seated rock fractures(referred to as DSRF hereafter)in valley slopes are uncommon geological phenomena that challenge our previous understanding of slope unloading processes.These fractures weaken the strength and integrity of the rock mass,potentially forming unstable block boundaries with significant volume,thereby affecting the stability of slopes,chambers,and dam abutments.DSRF has emerged as a critical environmental and engineering geological issue that hinders large-scale projects in deep canyon areas.Despite the attention and practical treatment given to DSRF in engineering practice,theoretical research on this topic still lags behind the demands of engineering applications.To garner widespread attention and promote the resolution of DSRF-related problems,this review aims to redefine DSRF through comprehensive data collection and analysis,engineering geological analogies,and field investigations,and provide a summary and analysis of the research progress on DSRF,along with future research directions.The study defines DSRF as the intermittent tension cracks or relaxation zones within a slightly weathered or fresh,and intact or relatively intact rock mass distributed below the surface unloading zones of a deep canyon slope,and should be distinguished from"loose rock mass"and"deep-seated gravitational slope deformations".The article provides an overview of the development and distribution,rupture characteristics,and genesis mechanism of DSRF.It proposes that DSRF is formed based on the fluvial deviation-undercutting evolution mode,wherein the energy accumulated in the rock mass is violently released when the river further down cuts the slope after the rock mass has undergone cyclical loadingunloading.However,further research is necessary to establish a comprehensive database for DSRF,refine exploration techniques,understand evolutionary processes,develop engineering evaluation methods,and predict the distribution of DSRF.
基金supported by the National Natural Science Foun-dation of China(Grant Nos.U2067203 and 42277140)Tsinghua University Initiative Scientific Research Program(Grant No.2022Z11QYJ006).
文摘Underground space creation and energy extraction, which induce unloading on rock fractures, commonly occur in various rock engineering projects, and rock engineering projects are subjected to high temperatures with increasing depth. Fluid flow behavior of rock fractures is a critical issue in many subsurface rock engineering projects. Previous studies have extensively considered permeability evolution in rock fractures under loading phase, whereas changes in fracture permeability under unloading phase have not been fully understood. To examine the unloading-induced changes in fracture permeability under different temperatures, we performed water flow-through tests on fractured rock samples subjected to decreasing confining pressures and different temperatures. The experimental results show that the permeability of fracture increases with unloading of confining pressure but decreases with loading-unloading cycles. Temperature may affect fracture permeability when it is higher than a certain threshold. An empirical model of fracture hydraulic aperture including two material parameters of initial normal stiffness and maximum normal closure can well describe the permeability changes in rough rock fracture subjected to loading-unloading cycles and heating. A coupled thermo-mechanical model considering asperity damage is finally used to understand the influences of stress paths and temperatures on fracture permeability.
基金The financial support from the National Natural Science Foundation of China(Nos.52174092,51904290,52004272,and 52274145)Natural Science Foundation of Jiangsu Province,China(Nos.BK20220157 and BK20200660)+2 种基金the Fundamental Research Funds for the Central Universities,China(No.2022YCPY0202)Xuzhou Science and Technology Project,China(Nos.KC21033 and KC22005)Yunlong Lake Laboratory of Deep Underground Science and Engineering Project,China(No.104023002)。
文摘A fractal model governing saw-tooth fractures was first introduced to replicate sandstone samples containing an inclined 3D penetrating rough fracture surface with various joint roughness coefficients(JRC).In conventional triaxial compression,the peak strength for fractured samples increased with both confining pressure and JRC.During the unloading confining pressure process,the normal stress of fractures declined but the shear stress increased,resulting in shear sliding of fractures.The shear displacement of fractures exponentially increased,and the positive normal displacement decreased gradually to negative values under coupling effects of shear contraction caused by normal stress and shear dilation due to climbing effects of fractures.Transition from quasi-static to dynamic sliding of the fractures was identified.The sliding resistance duration increased with confining pressure but decreased with JRC.After prepeak unloading,the fracture surfaces presented a more significant surface wear response and JRC values decreased by 1.70%–59.20%due to more remarkable asperity degradation compared with those after conventional triaxial compression.The theoretical model for shear strength of fractures was established through improving the Ladanyi&Archambault model by introducing the relations between normal stress and surface wear ratios of fractures,which agreed well with the experimental results.
基金This study has been partially funded by National Key Research and Development Program of China(Grant No.2020YFA0711800)the National Natural Science Foundation of China(Grant No.51979272)the Natural Science Foundation of Shandong Province,China(Grant No.ZR2021QE069).
文摘This study experimentally analyzes the nonlinear flow characteristics and channelization of fluid through rough-walled fractures during the shear process using a shear-flow-visualization apparatus.A series of fluid flow and visualization tests is performed on four transparent fracture specimens with various shear displacements of 1 mm,3 mm,5 mm,7 mm and 10 mm under a normal stress of 0.5 MPa.Four granite fractures with different roughnesses are selected and quantified using variogram fractal dimensions.The obtained results show that the critical Reynolds number tends to increase with increasing shear displacement but decrease with increasing roughness of fracture surface.The flow paths are more tortuous at the beginning of shear because of the wide distribution of small contact spots.As the shear displacement continues to increase,preferential flow paths are more distinctly observed due to the decrease in the number of contact spots caused by shear dilation;yet the area of single contacts in-creases.Based on the experimental results,an empirical mathematical equation is proposed to quantify the critical Reynolds number using the contact area ratio and fractal dimension.
基金This research was funded by the Natural Science Foundation of Sichuan,China(No.2022NSFSC1915)the National Natural Science Foundation of China(No.U19A2098)+1 种基金State Key Laboratory of Precision Blasting and Hubei Key Laboratory of Blasting Engineering,Jianghan University(No.PBSKL2022B06)the Fundamental Research Funds for the Central Universities。
文摘Three-dimensional rock fracture induced by blasting is a highly complex problem and has received considerable attention in geotechnical engineering.The material point method is firstly applied to treat this challenging task.Some inherent weaknesses can be overcome by coupling the generalized interpolation material point(GIMP)and the convected particle domain interpolation technique(CPDI).For the media in the borehole,unchanged GIMP-type particles are used to guarantee a homogenous blast pressure.CPDITetrahedron type particles are employed to avoid the fake numerical fracture near the borehole for the rock material.A blasting experiment using three-dimensional single-borehole rock was simulated to examine the applicability of the coupled model under realistic loading and boundary conditions.A good agreement was achieved between the simulation and experimental results.Moreover,the mechanism of three-dimensional rock fracture was analyzed.It was concluded that rock particle size and material parameters play an important role in rock damage.The reflected tensile waves cause severe damage in the lower part of the model.Rayleigh waves occur on the top face of the rock model to induce a hoop failure band.
基金supported by the National Research Foundation of Korea(KRF)grant funded by the Korea government(MSIT)(No.NRF-2019R1G1A1100517)the Basic Research and Development Project of the Korea Institute of Geoscience and Mineral Resources(KIGAM),which was funded by the Ministry of Science and ICT,Korea。
文摘A convenient approach was proposed by which to evaluate and monitor the permeability of a rock fracture by verifying the quantitative correlation between the electrical resistivity and permeability at laboratory scale.For this purpose,an electrical resistivity measurement system was applied to the laboratory experiments using artificial cells with the shape of a single rock fracture.Sixty experiments were conducted using rock fractures according to the geometry,aperture sizes,wavelengths,and roughness amplitudes.The overall negative relationship between the normalized electrical resistivity values and the aperture sizes directly linked with the permeability,was well fitted by the power-law function with a large determination coefficient(≈0.86).The effects of wavelength and roughness amplitude of the rock fracture on the electrical resistivity were also analyzed.Results showed that the electrical resistivity was slightly increased with decreasing wavelength and increasing roughness amplitude.An empirical model for evaluating the permeability of a rock fracture was proposed based on the experimental data.In the field,if the electrical resistivity of pore groundwater could be measured in advance,this empirical model could be applied effectively for simple,quick monitoring of the fracture permeability.Although uncertainty may be associated with the permeability estimation due to the limited control parameters considered in this research,this electrical resistivity approach could be helpful to monitor the rock permeability in deep underground facilities such as those used for radioactive waste repositories or forms of energy storage.
基金Peck Tech Consulting Ltd.for its support and help
文摘Accurate determination of rock mass properties is a critical part of open-pit mine planning activities to enable more prescriptive blast designs to achieve improved loading and hauling and downstream process efficiency. Better and more accurate blast practice that delivers enhanced outcomes(better fragmentation, improved diggability, less dilution, etc.) is a critical and fundamental element of being able to achieve an effective Mine-Mill approach at a mining operation. Based on previous work, it has been demonstrated that an accurate representation of the rock mass properties can be obtained from the analysis of variations in blasthole drill performance as derived from measurement while drilling(MWD) systems when using tricone bits. This paper further investigates how monitored rate of penetration,pulldown force, rotary torque, rotary speed and bailing air pressure responses can be used to determine the presence of open and partially open fractures having varying dip angles. Based on a correlation of geophysically measured fracture logs and monitored drill performance variables in the same blastholes, the results show that the latter responses can accurately determine open versus closed fractures. The results also identified that variations in rate of penetration and rotary torque show the most sensitivity in the presence of open fractures that intersect a vertical blasthole at near orthogonal angles.
基金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.
基金supported by State Key Laboratory of Earthquake Dynamics,China
文摘In order to improve our understanding of rock fracture and fault instability driven by high-pressure fluid sources, the authors carried out rock fracture tests using granite under a confining pressure of 80 MPa with fluid injection in the laboratory. Furthermore, we tested a number of numerical models using the FLAC^(3D) modeling software to find the best model to represent the experimental results. The high-speed multichannel acoustic emission(AE) waveform recording system used in this study made it possible to examine the total fracture process through detailed monitoring of AE hypocenters and seismic velocity.The experimental results show that injecting high-pressure oil into the rock sample can induce AE activity at very low stress levels and can dramatically reduce the strength of the rock. The results of the numerical simulations show that major experimental results, including the strength, the temporal and spatial patterns of the AE events, and the role of the fluid can be represented fairly well by a model involving(1) randomly distributed defect elements to model pre-existing cracks,(2) random modification of rock properties to represent inhomogeneity introduced by different mineral grains, and(3)macroscopic inhomogeneity. Our study, which incorporates laboratory experiments and numerical simulations, indicates that such an approach is helpful in finding a better model not only for simulating experimental results but also for upscaling purposes.
基金supported by the National Petroleum Major Projects(No.2017ZX05019-005)CNPC Fundamental Research Projects(No.2016A-3605)
文摘To identify refl ector fractures near borehole by using dipole-source refl ected-shearwave logging, we need to understand the relation between the amplitude of the refl ected shear wave and the source radiation, borehole conditions, and attenuation owing to the surrounding formations. To assess the effect of these factors on the amplitude of the refl ected waves, we first studied the radiation performance and radiation direction of the dipole source in fast, medium, and slow formations by using the asymptotic solution in the far fi eld of the borehole. Then, the relation between the fracture parameters, and the refl ected-shear-wave amplitude as well as the ratio of the refl ected-shear-wave amplitude to the direct-wave amplitude (relative amplitude, RA) was evaluated by the three-dimensional fi nite-difference (3D FDTD) method. Finally, the fracture detection capability of the dipole reflected-shear-wave logging tool in different formations was analyzed by using the RA. The results suggest that the radiation amplitude of the SH-wave in the slow formation is weaker than those in the fast and medium formations, and the amplitude of the refl ected shear wave is lower. However, the RA in the slow formation is close to or even greater than in the fast and medium formations, which means that dipole-source shear-wave logging has the same or even better fracture detection capability in the slow formation as in the fast and medium formations. In addition, when RA is small, there is a good correlation between the RA and the various fracture parameters in the different types of formation, which can be used in determining the lower limit of the fracture parameters identifi ed by refl ection logging.
基金supported by Key Innovation Team Program of Innovation Talents Promotion Plan by Ministry of Science and Technology(MOST)of China(Grant No.2016RA4059)Science and Technology Project of Yunnan Provincial Transportation Department(Grant No.25 of 2018)Shanghai Science and Technology Committee Program(Grant No.20dz1202200).
文摘This paper presents a novel integrated method for interactive characterization of fracture spacing in rock tunnel sections.The main procedure includes four steps:(1)Automatic extraction of fracture traces,(2)digitization of trace maps,(3)disconnection and grouping of traces,and(4)interactive measurement of fracture set spacing,total spacing,and surface rock quality designation(S-RQD)value.To evaluate the performance of the proposed method,sample images were obtained by employing a photogrammetrybased scheme in tunnel faces.Experiments were then conducted to determine the optimal parameter values(i.e.distance threshold,angle threshold,and number of fracture trace grouping)for characterizing rock fracture spacing.By applying the identified optimal parameters involved in the model,the proposed method could lead to excellent qualitative results to a new tunnel face.To perform a quantitative analysis,three methods(i.e.field,straightening,and the proposed method)were employed in the same study and comparisons were made.The proposed method agrees well with the field measurement in terms of the maximum and average values of measured spacing distribution.Overall,the proposed method has reasonably good accuracy and interactive advantage for estimating the ultimate fracture spacing and S-RQD.It can be a possible extension of existing methods for fracture spacing characterization for two-dimensional(2D)rock tunnel faces.
基金financially supported by the National Key R&D Program of China(Grant No.2019YFB1504103)the China Postdoctoral Science Foundation(Grant Nos.2019TQ0174)。
文摘Numerous intersected rock fractures constitute the fracture network in enhanced geothermal systems.The complicated convective heat transfer behavior in intersected fractures is critical to the heat recovery in fractured geothermal reservoirs.A series of three-dimensional intersected fracture models is constructed to perform the flow-through heat transfer simulations.The geometry effects of dead-end fractures(DEFs)on the heat transfer are evaluated in terms of intersected angles,apertures,lengths,and the connectivity.The results indicate that annular streamlines appear in the rough DEF and cause an ellipse distribution of the cold front.Compared to plate DEFs,the fluid flow in the rough DEF enhances the heat transfer.Both the increment of outlet water temperatureΔToutand the ratio of heat production Qrpresent the largest at the intersected angle of 90°while decline with the decrease of the intersected angle between the main flow fracture(MFF)and the DEFs.The extension of the length of intersected DEFs is beneficial to heat production while enhancing its aperture is not needed.Solely increasing the number of intersected DEFs induces a little increase of heat extraction,and more significant heat production can be obtained through connecting these DEFs with the MFF forming the flow network.
基金Project (50374073)supported by the National Natural Science Foundation of China project (2002032256)supported bythe Postdoctor Science Foundation of China
文摘Anti-plane punch-through shear test and anti-planefour-point bending test are used to study the crack initiation and propagation under anti-plane shear (Mode Ⅲ) loading. The tensile and shear stresses at the crack tip are calculated by finite element method. The results show that under Mode Ⅲ loading the maximum principal stress σ1 at crack tip is smaller or a little larger than the maximum shear stress τmax. Since the tensile strength of brittle rock is much lower than its shear strength, σ1 is easy to reach its critical value before τmax reaches its critical value and thus results in Mode I fracture. The fracture trajectory is helicoid and the normal direction of tangential plane with the fractured helicoid is along the predicted direction of the maximum principal stress at the notch tip. It is further proved that Mode Ⅰ instead of Mode Ⅲ fracture occurs in brittle rock under Mode Ⅲ loading. The fracture mode depending on the fracture mechanism must be distinguished from the loading form.
文摘The joint roughness coefficient (JRC) is an important geometric parameter in rock mechanics, seismo-geology, engineering geology, and structural geology. In recent years, the fractol geometry has been widely applied in researching JRC In this paper, a simpler fractal model for estimating JRC was developed, whick can be applied in simulating rock discontinuity profiles, aud the rougher the surface, the greater the fractal dimension D The epirical equation of JRC and D was established according to the model’ The JRC values of six joint profiles were estimated according to Eq (10),which are ratker approximate to the typical JRC values; so the method may be applied to estimate the roughness of the fractured surfaces of natural rocks For a give profile the estimates of D and JRC can be obtained from Eqs (5) and (l0) with only one parameter e’ (the average angle of the profile fluctuation) being needed.
文摘The application of the non-explosive expansion material (NEEM) is widely used as the controlled fracture method in quarry mining,especially in hard rocks.The pressure of NEEM is an important parameter in causing rock fracture.An empirical model based on hole spacing was developed to determine the pressure of NEEM in the rock fracture process.Primarily,the empirical model was developed by the mathematical method,utilizing dimensional analysis.Then,the Phase2 code,which is based on the finite element method,was utilized to predict crack growth in rocks.The results of numerical analysis show slight deviations from the empirical model.Hence,the polynomial regression analysis was used to modify the model.Finally,the modified model shows a good agreement with the results gained from numerical modeling.
基金National Natural Science Foundation of China(No.41372239)Graduate Innovation Fund of Jilin University(No.2015096)
文摘Through seepage tests under different loading and unloading confining pressures and different hydraulic gradients,the authors studied the effects of stress states and stress history on fracture permeability evolution for single granite fracture and sandstone fracture. The results show that there exists a linear relationship between the seepage discharge and osmotic pressure in sandstone fissure under each level of confining pressure. With the increasing in the confining pressure,the permeability of the fracture decreases,but the decreasing rate is changeing. During the unloading process,the fracture seepage velocity cannot be fully recovered to the size of the loading process. Therefore,in the unloading process of the confining pressure,the recovery of fracture permeability shows obvious hysteresis effects. The flow rate of the fracture remains unchanged during five cycles of loading and unloading processes of the confining pressure. In each cycle,the evolution character of the flow rate with the confining pressure remains unchanged. These experiments show that the seepage characteristics of sandstone and granite fractures are not the same under the same stress state.
文摘It is observed that the parameter of seismic inhomogeneous degree (GL value) calculated from the earthquake catalog shows obvious abnormal changes prior to strong earthquakes, indicating the state change of local seismic activity. This paper focuses on the mechanism for the abnormal changes of the GL values based on the sequences of acoustic emission for three types of rock samples containing macro-asperity fracture; compressional en-echelon fracture and model-III shear fracture. The results show that for the three types of rock samples, there are continuous abnormal changes of GL value (>1) just before the non-elastic deformation occurs or during the process of nucleation prior to the instability. Based on the experimental results, it seems that the process of creep sliding and resistance-uniformization along fault zone is the possible mechanism for the abnormal changes of GL value before rock fractures.
文摘Fracture and ground vibration of rock subjected to different decoupling decked charges are investigated based on the numerical simulation.The dynamic pressure value is studied,which demonstrates that simulation of fracture zone is feasible.Attenuation index of dynamic pressure is 2.06,2.05 and 1.93 for air,water and sand intervals respectively.The small attenuation of sand interval results in the large ground vibration.The predicted vertical vibration waveform and peak particle velocities(PPV) in far-field are in agreement with the monitoring results.The results show that the air and water decked charges can improve the effect of rock fracture in near-field and reduce ground vibration in far-field.
基金Projects(50809069,10872207) supported by the National Natural Science Foundation of ChinaProject(2009CAD017) supported by Outstanding Youth of Hubei Province,China
文摘A true triaxial apparatus which is composed of three units was presented.The apparatus allows for investigations on deformation and seepage behaviors of a single rock fracture subjected to lateral stress and normal stress.The first unit has three jacks which can apply loads independently in three orthogonal directions.The second unit is used to supply water inflow,control seepage pressure and measure flow velocity in real time.The third unit is for measuring the normal deformation of rock fractures.Some tests for investigating the normal deformation and seepage behaviors of rock fractures subjected to normal and lateral loads on hard granite specimens with an artificial persistent fracture,were introduced.The results show that both the normal deformation and the hydraulic conductivity are influenced not only by the normal stress but also by the lateral stress.It is also shown that the aperture and the hydraulic conductivity decrease with the increasing normal stress but increase with the increasing lateral stress and both the aperture and the hydraulic conductivity obey exponential relationships with the normal stress and the lateral stress.
