基于声发射的岩石破裂应力场动态反演

为获取岩石在破裂过程中的动态应力场,首先进行室内物理试验获取声发射数据,在此基础上基于能量耗散原理,利用实测声发射信息建立起岩石细观损伤表征方程,并通过FLAC3D内嵌FISH语言进行二次开发,自动搜索声发射损伤影响范围内岩石单元并对其力学参数进行弱化,计算得到较为真实的岩石破裂动态应力场。结果表明,计算得到的岩石破裂应力场与实测结果有良好的对应关系,验证了岩石细观损伤表征方法的合理性。同时,获得的岩石动态破裂应力场不仅可以从岩石力学的角度解释破坏产生的原因,也对预测岩石下一步破坏具有指导意义。但这种方法仍具有一定的局限性,需要进一步的探索研究。     

Inverse calculation of in situ stress in rock mass using the Surrogate-Model Accelerated Random Search algorithm

A new approach is presented that incorporates the Surrogate Model Accelerated Random Search (SMARS) Algorithm to inversely determine the overall stress state in a rock mass based on sparse stress measurements. The SMARS algorithm relies on a random search component to maintain global search capabilities, while using the surrogate-model method to accelerate convergence to a solution estimate. Two examples sets are carried out in this presentation to display the utility of the proposed SMARS-based inverse characterization process. The first example set compares the SMARS algorithm with two other popular methods, namely a multiple linear regression analysis method and a neural network method, to estimate the in situ stress field for a simple numerically simulated test case. The results of the numerical testing verify that SMARS provides a relatively stable approach and gives rise to a relatively high accuracy and efficiency (i.e., with less computational expense) compared to the other contemporary approaches considered. Finally, an example is shown for utilizing the SMARS approach for in situ stress estimation based on an actual underground mine located in Pennsylvania. The SMARS results are shown to produce a realistic estimate of the distribution of stress within the area investigated, and overall, the approach has potential for practical use in realistic scenarios to efficiently and accurately estimate in situ stresses in rock mass.     

Effects of size and strain rate on the mechanical behaviors of rock specimens under uniaxial compression

Size and strain rate are two key factors that dramatically influence the estimation of rock mechanical behaviors. To better understand the effects of size and strain rate on measured rocks, rock specimens with six different sizes were tested at six different strain rates under uniaxial compression using the MTS 815 Rock Mechanics Test System. Having determined that the size and strain rate significantly affect the peak strain, peak stress, elastic modulus, acoustic emission (AE), and failure pattern of the rock specimens, the relation was established between the strength and the size and strain rate of red sandstone. And the variation was revealed among the size and strain rate, the AE, and the failure pattern. It turned out that the peak stress was negatively correlated with the rock size and was positively correlated with the strain rate. When the length to diameter ratio (L/D) of the rock specimen was less than 2.0, the AE appeared mildly. The AE quantities gradually increased before the peak stress, and then sharply decreased after the peak stress. The failure pattern of the rock specimen was relatively complicated, with a fracture plane appearing along the axial direction. Conical failure type was also presented. When the L/D ratio of the specimen was greater than 2.0, the AE characteristics of red sandstone showed the radical model. There were relatively few AE rings before the peak stress. But the AE rings increased suddenly and dramatically during the peak stress. The rock specimens primarily failed with a single shear plane. Moreover, with an increase in the strain rate, the AE activities were enhanced and the AE quantities increased. When the strain rate of the rock specimen was less than 5.0?×?10^?4/s, the rock specimen failed with a shear or tensile-shear pattern. And when the strain rate was greater than 5.0?×?10^?4/s, the rock specimen tended to fail in a conical pattern.     

A new method of modelling the rock micro-fracturing process in double-torsion experiments using neural networks

Microfracturing of rock is a complicated damage evolution process. Inaccurate prediction of micro-fracturing behaviours suggests a need for the development of a better modelling method. Analysis of acoustic emission (AE) measurements in double-torsion tests indicates that micro-fracturing behaviours during the loading stage have fractal time structures. This fractal behaviour can be described by C(t) ∝ t^D, where D is the correlation exponent, t is the time and C(t) is the correlation integral. Furthermore, by utilizing measured AE data, a new method has been developed to model the AE behaviours of micro-fracturing in rock, in air, and following soaking in water and in a chemical solution of DTAB. The neutral models NN (10,21,2) and NN (10,20,2) were found to describe reasonably well the AE behaviours of micro-fracturing in rock under air and DTAB conditions, and water conditions, respectively. The cumulative AE events and the cumulative AE counts predicted by the neural models agreed well with those measured in experiments. Copyright © 1999 John Wiley & Sons, Ltd.     

The determination ofin situ rock thermal properties and the simulation of climate in an underground mine

Summary This paper describes the results from a continuing research project started in 1984, which is investigating improved methods for predicting the underground climate in hot mines. The Homestake Gold Mine in South Dakota was the site for testing methods, based on radial heat flow theory, to measurein situ values of the thermal conductivity and diffusivity of rock. The methods used proved to be quite within the capability of mine staff to carry out for themselves. The results were acceptable and verified thatin situ thermal rock properties can differ considerably from laboratory determined values, in this case by a factor of approximately two.A mine climate simulation program, CLIMSIM, was applied to case studies from the Homestake mine using the measured thermal rock parameters. The results showed good correlation between the measured and simulated values.     

An anisotropic strength criterion for jointed rock masses and its application in wellbore stability analyses

In this paper, an anisotropic strength criterion is established for jointed rock masses. An orientation distribution function (ODF) of joint connectivity, is introduced to characterize the anisotropic strength of jointed rock masses related to directional distributed joint sets. Coulomb failure condition is formulated for each plane of jointed rock masses by joint connectivity, where the friction coefficient and cohesion of the jointed rock mass are related to those of the intact rock and joint and become orientation dependent. When approximating joint connectivity by its second‐order fabric tensor, an anisotropic strength criterion is derived through an approximate analytical solution to the critical plane problem. To demonstrate the effects of joint distribution on the anisotropic strength of jointed rock masses, the failure envelopes are worked out for different relative orientations of material anisotropy and principal stress axes. The anisotropic strength criterion is also applied to wellbore stability analyses. It is shown that a borehole drilled in the direction of the maximum principal in situ stress is not always the safest due to the anisotropic strength of the jointed rock mass. Copyright © 2007 John Wiley & Sons, Ltd.     

Effect of rock discontinuities on certain rock strength and fracture energy parameters under uniaxial compression

Summary Five series of test blocks of Pendeli marble with artificially created discontinuities of different crack densities (simulating three mutually orthogonal joint sets) were tested in uniaxial compression in order to study the effect of discontinuities on: (a) the compressive strength and the modulus of elasticity, and (b) certain fracture energy parameters expressed by the ratio W _A/W _V, where W _A is the surface energy and W _V the volume elastic strain energy. Mathematical relationships are derived similar to those suggested by other authors relating strength parameters to crack densities. Such relationships clearly show a reduction in strength with increased crack density. The experimental results obtained permit the extension of Persson's relation (which refers to ideal intact rock) to the more realistic case of discontinuous rock mass by introducing the appropriate term that takes into consideration the effect of rock mass discontinuities on the energy ratio W _A/W _V. A comparison between laboratory results and field observations was subsequently carried out assuming the rock mass to behave as a linearly elastic material, obeying the Hoek and Brown failure criterion. This comparison showed that laboratory results can be extended to larger scale. Furthermore, in order to predict the in situ strength and stability of a rock mass in uniaxial compression (which is of major importance in underground excavations) certain concepts are proposed based on laboratory tests, in situ investigations and first principles of linear elastic fracture mechanics.     

A diametral deformation method forin situ stress and rock property measurement

Summary Failure deformation of a borehole boundary in response to hydraulic loading via a plastic cylinder is directly related to both the initial stress state and the material properties of the ground. These relationships are utilized for determination ofin situ stress states and material properties of complex geological formations, including fractured hard rock and soft ductile media, by directly measuring diametral deformation of the borehole in relation to the applied pressure. This new method was first developed by means of finite element simulation models, and later confirmed by laboratory model studies. This method was then field tested in an outcropping of homogeneous sandstone and later applied to three underground geological formations: highly stratified shale, rock salt, and fractured green schist. Some of the results of these field applications are presented here to illustrate the validity and usefulness of this method, especially for complex ground where the conventional methods of overcoring and hydrofracturing may be ineffective. The instrumentation of this method involves a borehole probe which utilizes electrically powered hydraulic pumping and computerized data acquisition for rapid measurement and on-site analysis. The significance of the new method is its ability to obtain input data ofin situ stress states and material properties for finite element analysis of earth structures in complex ground.     

基于三维扫描和三维雕刻技术的岩石结构面原状重构方法及其力学特性

为了客观重复制作与原岩结构面的岩性、表面形态和力学性质完全一致的试样,集成三维扫描和三维雕刻技术,提出了原岩结构面试样重构新方法。采用该方法制作了具有相同表面形态的大理岩结构面试样,开展了不同法向力下的直剪试验,并采用声发射技术对结构面剪切破坏过程进行了监测。结果表明：重构结构面试样与原始结构面试样具有高度相似性;结构面剪应力-剪切位移曲线可分为两类：剪胀滑移型和剪断跌落型;剪切破坏后结构面可明显地划分为两个区域：剪胀脱开区和剪断磨损区;提出了结构面剪切强度新公式,理论公式与试验结果更接近;结构面剪切过程中剪应力、能量率、计数、撞击率随时间的变化规律具有较好的一致性,声发射定位位置与剪切过程中结构面损伤破坏区域基本吻合。     

岩体的声发射特征试验研究

在单轴压缩条件下,讨论了4种岩石的声发射特征;在静态胀裂剂作用下,讨论了岩体的声发射特征。研究表明,在低应力阶段,岩石几乎没有声发射活动,一般当强度达到80 %以上,即临近破坏时,声发射活动才显著增加;岩体破坏的声发射过程分为4个阶段,即初始区(Ⅰ)、剧烈区(Ⅱ)、下降区(Ⅲ)和沉寂区(Ⅳ);有的岩石的凯瑟效应不明显;岩石的声发射主频与岩石的强度有关,强度越高,主频也越高;随着应力的增加,岩石的声发射主频不会平移;应用1 kHz的探头,能满足工程岩体的稳定性分析和冒落预报的要求。     

An approach to large-scale field stress determination

The construction of stable structures in rock masses requires knowledge of the in situ stresses at the scale of excavations. However, the measurements obtained by the conventional overcoring technique are related to a small scale (centimetres). To extrapolate them to the scales of interest to rock mechanics (from meters to kilometres) requires a large number of individual stress measurements, followed by statistical analysis to avoid a considerable scatter of the measured values. In this paper, a method is proposed based on (a) large-scale surface stress and modulus measurements using the cylindrical jack method complemented by a special measuring scheme and then (b) back analysis for a given excavation shape. The method allows the simultaneous reconstruction of the stress components at the scale of excavation. A numerical simulation for a cylindrical excavation in an isotropic rock mass demonstrates the high accuracy and robustness of the method. The presence of a fractured zone surrounding the excavation can hamper the stress reconstruction, hence special measures should be taken to conduct the measurements in competent rock.     

岩石声发射混沌特征分析

岩石等脆性材料在加载过程中,随着荷载的增加,材料内部的微裂纹产生、扩展并伴随着声发射现象的发生。声发射是研究脆性材料破坏的良好工具。对砂岩、细砂岩和高丽山砂岩3类岩石进行了声发射试验,记录了加载及破坏过程中产生的声发射信号,并且采用混沌动力学理论研究了3类岩石的声发射活动规律,计算了岩石的关联维数和最大Lyapunov指数。研究结果表明,岩石加载及破坏过程具有混沌特征,用相空间重构法可以较好地揭示岩石破坏过程的动力学特征,这为混沌理论在岩石、岩体声发射其他领域的研究及应用奠定了基础。     

三维地应力计算模型研究

简要回顾了现有的三维地应力实测方法在实际工程应用中存在的局限性。对钻孔局部壁面应力全解除法进行了系统的研究,从理论上说明了由它确定三维地壳应力张量的可行性。采用线弹性有限元法对钻孔局部壁面应力解除过程做了数值模拟研究,分析了应力解除曲线的变化特征,给出了达到完全应力解除时解除深度和被解除岩芯直径之间的定量关系,还进一步分析了环形解除槽引起的应力集中对其附近局部孔壁应力分布的影响,从而为在钻孔孔壁的局部壁面上如何合理的布置应变计（花）提供了参考依据。最后对本课题的研究前景做了展望。     

基于逾渗模型的盐岩损伤与破坏研究

逾渗是一种分形模型,利用声发射点定位盐岩内部破坏网格来建立逾渗模型,可以分析三轴压缩条件下盐岩逾渗特征以及损伤演变发展。研究发现,各试样逾渗模型团簇(cluster)数与最大团簇占有率关系曲线斜率在70%应力峰值后基本相等,利用最大团簇沿试样轴向的延伸来描述裂缝的扩展规律,并通过最大团簇在轴向的延伸终点得到了各试样逾渗的临界破坏比率,它对于研究渗透陡增点非常重要。提出利用逾渗关联长度(关联长度?代表处于同一团簇中的两个点的平均距离)确定试样损伤起始点,并通过计算逾渗模型破坏比率得到损伤变量,它和基于AE振铃计数以及AE能量计算得到的损伤变量结果相近。研究结果表明,逾渗可以形象地表述岩石内部破裂过程和损伤情况,为研究岩石破坏失效及裂缝衍生发展提供了新的思路。     

持续开挖效应下结构面剪切力学性质与能量特征研究

基于地下岩体工程开挖扰动的实际情况,开展考虑持续开挖效应,充分反映岩体结构面应力调整过程的结构面剪切试验更具理论意义和工程应用价值。采用人工劈裂方法制备岩体结构面,开展了常规应力路径和持续开挖效应下的结构面剪切试验,系统研究了两种条件下结构面剪切力学性质、声发射特征和能量的演化规律。研究结果表明：开挖扰动荷载强度越大,结构面发生剪切破坏时的剪应力降整体上越大,但持续开挖效应下的剪应力降最大值仅为常规应力路径下剪应力降的48.57%;考虑持续开挖效应的结构面剪切过程中声发射活动主要集中于结构面受剪破坏和剪应力降产生时,且声发射活动强度与开挖扰动强度正相关,声发射振铃计数变化率极值明显小于常规应力路径下的相应值;持续开挖效应下试样裂纹发育密集度、结构面磨损区域及破坏程度均随开挖扰动强度增加而增加,但试样裂纹发育的密集程度和结构面破坏程度相对常规应力路径下的情况轻;考虑持续开挖效应的结构面发生剪切破坏时弹性应变能随扰动荷载增大而增大,均低于常规应力路径条件下的弹性应变能值,持续开挖效应降低了结构面剪切破坏的强度,但更易引发结构面发生剪切破坏。     

Three-dimensional in situ stress evaluation using a new under-coring technique: the Tseng-Wen Reservoir Transbasin water tunnel

The Tseng-Wen Reservoir Transbasin water tunnel has broken the historical record of the maximum overburden of a tunnel within Taiwan. The three-dimensional in situ stresses at the tunnel were evaluated using acoustic emission (AE) and deformation rate analysis (DRA) with an under-coring technique. This is the first time that the AE and DRA have been used in evaluating an in situ stress in Taiwan. As cores drilled in different directions are required for the three-dimensional stress evaluation, but the in situ borehole drilling is usually limited to the vertical direction, small samples 2.28 cm in diameter and 4.56 cm in length were under-cored from the borehole cores in six directions to evaluate the in situ stresses. Cyclic uniaxial compressive loadings were artificially applied to the Changchikeng sandstones under-cored from the outcrops. The experimental results validate the definitiveness of using under-cored samples at different pre-stresses with AE and DRA. Then, cores taken from the Tseng-Wen Reservoir Transbasin water tunnel site at a great depth of 416–739 m underground were under-cored. The experimental results indicate that the major principle in situ stress is from 30.8 to 39 MPa and is in the directions of SE–NW and SW–NE at different depths.     

Numerical computations of rock dissolution and geomechanical effects for CO2 geological storage

The paper is motivated by the long‐term safety analysis of the CO₂ geological storage. We present a methodology for the assessment of the geomechanical impact of progressive rock dissolution. The method is based on the use of X‐ray tomography and the numerical dissolution technique. The influence of evolution of the microstructure on the macroscopic properties of the rock is analysed by using periodic homogenization method. The numerical computations show progressive degradation of all components of the stiffness (orthotropic) tensor. Moreover, the evolution of associated mass transfer properties (as tortuosity and conductivity tensors), by using the periodic homogenization method, is also calculated. The correlation between the mechanical parameters and the transfer properties during the dissolution process is presented. The results show that the highest increase of the hydraulic conductivity (in direction Y) is not associated with the highest decrease of Young modulus in this direction. Moreover, the highest decrease of Young modulus (in the direction X) is not associated with percolation in this direction. Finally, an incremental law to calculate settlement, in case of a rock with evolving microstructure, is proposed. The solution of the macroscopic settlement problem under constant stress and drained conditions showed that the geomechanical effects of the rock dissolution are rather limited. Copyright © 2014 John Wiley & Sons, Ltd.     

A new method for identification of protocatagenesis, mesocatagenesis, and apocatagenesis zones in a terrigenous sedimentary cover by applying layer seismic velocities

This method pertains to oil and gas geology and to geology of sedimentary basins prospective for oil and gas. It includes identifying catagenesis zones in drilled areas within the sedimentary cover of the basin based on assay results for drill cores and cuttings using the standard methods. Analysis is primarily made in order to determine rock catagenesis based on the optical characteristics of vitrinite. A correlation between catagenesis zones and layer seismic velocities obtained from regional and exploration seismic data is made for a drilled area. Both the layer seismic velocities and the degree of rock catagenesis increase with depth under the influence of increasing rock density. Correlations between layer velocities and the degree of rock catagenesis have been established. The following ratios have been determined for the Scotian shelf, Canada, and the Barents shelf, Russia: a protocatagenesis zone (the cap) corresponds to layer seismic velocities (V _lay) of 1.5–3.3 km/s, a mesocatagenesis zone (the principal hydrocarbon generation area) corresponds to V _lay of 3.3–5.0 km/s, and an apocatagenesis zone (an area with a very low hydrocarbon potential) corresponds to V _lay of over 5.0 km/s. An advantage of the new method of identification of catagenesis zones is that it can be used prior to drilling. Its conceptual originality and cost efficiency lie precisely in this.     

岩石声速与其损伤及声发射关系研究

首先,建立了单轴压缩过程岩石损伤参量、应变与声速之间的定量关系式,分析了不同均质度对单轴压缩过程岩石声速的影响,结果表明,随着均匀度的增加,单轴压缩过程中声速由平缓变化到急剧变化,这与已有的岩石声发射数值模拟分析结果是一致的;其次,根据建立的声速与应变的关系公式,通过单轴压缩过程岩石声速与应变实测结果的回归分析,得到了具有较高精度的回归方程,从而通过试验验证了所建立的关系式的正确性;最后,从损伤的角度讨论了单轴压缩过程岩石声速与声发射的关系,得出了Kaiser点应位于声速初始下降点附近的结论,为岩石声发射测量地应力试验中Kaiser点的确定提供了新的方法。     

Limiting equilibrium versus BS3D analysis of a tetrahedral rock block

In order to analyze the stability of a 3D rock block using a limiting equilibrium method, some simplifications must be introduced, such as neglecting the effects of deformability of the rock block, fractures, and the surrounding rock mass, the progressive failure and the mobilization of shear strength of fractures, and changes in in situ stresses with block displacement. These limitations may cause some errors when calculating the factor of safety of the block. This paper quantifies the error caused by these simplifications comparing the results against those obtained with a numerical method (BS3D). In addition, the effects of the normal stiffness of the fractures, dilatancy, the tunnel radius, and the block size on stability of the tetrahedron are investigated using the numerical tool to demonstrate the importance of the effects eliminated by the limiting equilibrium approach. Copyright © 2010 John Wiley & Sons, Ltd.