A New Method for In-situ Non-contact Roughness Measurement of Large Rock Fracture Surfaces

Summary This paper presents a new method for in-situ non-contact measurements of fracture roughness by using a total station (TS). The TS is a non-reflector geodetic instrument usually used for measuring control points in surveying and mapping. By using a special-developed program, the TS can be used as a point-sensor laser scanner to scan a defined area of the fracture surface automatically, in field or in laboratory, at a distance away from the target surface. A large fracture surface can be automatically scanned with a constant interval of the sampling points, both within a defined area or along a cross-section of the exposed rock face. To quantify fracture roughness at different scales and obtain different densities of the scanned points, the point interval can be selected with the minimum interval of 1 mm. A local Cartesian co-ordinate system needs to be established first by the TS in front of the target rock face to define the true North or link the measurements to a known spatial co-ordinate system for both quantitative and spatial analysis of fracture roughness. To validate the method, fracture roughness data recorded with a non-reflector TS was compared with the data captured by a high-accuracy 3D-laser scanner. Results of this study revealed that both primary roughness and waviness of fracture surfaces can be quantified by the TS in the same accuracy level as that of the high accuracy laser scanner. Roughness of a natural fracture surface can be sampled without physical contact in a maximum distance of tens of meters from the rock faces. Received May 24, 2001; accepted July 24, 2002; Published online November 19, 2002     

断层相关褶皱概念模型中的裂缝域

断层相关褶皱对裂缝的控制作用主要体现在岩层的几何形态对裂缝方位的影响和变形对岩石破裂的控制程度。本文在断层相关褶皱概念模型的基础上提出了裂缝域,并以Suppe的平行膝折概念模型为对象,利用几何学方法对模型中区域性裂缝的裂缝域进行分区描述,并提出了断层相关褶皱的曲线概念模型中区域性裂缝的分布规律。以新疆卡普沙良河一处断层转折褶皱剖面上裂缝的实测数据为依据,对模型进行验证,发现测量结果同模型基本一致。该方法为前陆冲断带地区断层相关褶皱对裂缝控制作用的定量预测提供了一种全新的思路。     

Mapping hydraulic fractures from tiltmeter data using the ensemble Kalman filter

Hydraulic fracturing involves the initiation and propagation of fractures in rock formations by the injection of pressurized fluid. The largest use of hydraulic fracturing is in enhancing oil and gas production. Tiltmeters are sometimes used in the process to monitor the generated fracture geometry by measuring the fracture‐induced deformations. Fracture growth parameters obtained from tiltmeter mapping can be used to study the effectiveness of such stimulations. In this work, we present a novel scheme that uses the ensemble Kalman Filter (EnKF) to assimilate tiltmeter data using a simple process model to describe the evolution of fracture growth parameters, and an observation model that maps the fracture geometry with the observed tilt. The forward observation model is based on the analytical solution for computing the displacements and tilts due to a point source displacement discontinuity in an elastic half‐space developed by Okada 1 . The displacement and tilts for any given fracture geometry are then obtained by numerical integration of this solution, by considering multiple point sources to be located at the quadrature points. The proposed method is validated using synthetic data sets generated from polygon and elliptical shaped fracture geometries. Finally, real data from a field site, where asymmetry was measured from the intersections of the hydraulic fracture with offset boreholes, have been analyzed. Preliminary results show that, in addition to extracting the fracture dip, orientation, and volume, the procedure is able to satisfactorily predict fracture growth parameters when the fracture is relatively close to the tiltmeter array and provides some insight into the development of asymmetry when the measurements are relatively far from the fracture plane. Copyright © 2015 John Wiley & Sons, Ltd.     

Advanced Geostructural Survey Methods Applied to Rock Mass Characterization

The location and orientation of rock discontinuities, which are traditionally obtained from geological surveys with obvious drawbacks (safety, rock face accessibility, etc.), may also be derived from a detailed and accurate photogrammetric or laser scanning survey. Selecting from the point cloud determined on the rock face a set of points distributed on a particular discontinuity, location, dip, and dip direction can be computed from the least-squares estimate of the plane interpolating the set of points. Likewise, the normal vector to the surface may be computed from an interpolation or approximation of the surface by appropriate functions. To become a real alternative (both in terms of productivity as well as accuracy) to a traditional survey, interactive or automated software tools are necessary, to allow the efficient selection of the point sets on the discontinuities or the interpretation of the normal vector pattern. After introducing the two best technologies available today for data acquisition and their performance, the paper presents an approach, based on the random sample consensus (RANSAC) procedure, to the segmentation of the point cloud into subsets, each made of points measured on a discontinuity plane of the rock face. For each subset, the plane’s equations coefficients are first determined by robust estimation and then refined by least-squares estimation after outlier removal. The segmentation algorithm has been implemented in RockScan, a software tool developed to facilitate the interaction with the point cloud in the identification of the discontinuities; rather than using the three-dimensional (3D) data, selection of regions of interest is performed on oriented images of the rock face. Finally, application of RockScan to four different test sites is discussed and results presented. The sites differ in size (from tens to hundreds of meters), rock surface characteristics, and the technology used to produce the point cloud (in three cases photogrammetry, in the fourth laser scanning), giving the opportunity to test the methodology in different contexts. In the first and in the fourth site an extensive traditional survey has been performed, providing reference data to validate the RockScan results.     

Strike Standard Deviation for Shallow-dipping Rock Fracture Sets

Summary The well-known problem that the observed strike of a shallow-dipping rock fracture can vary widely, even if one re-measures the strike in the same location, is dealt with in this paper. The overall purpose of the paper is to place the statistical estimation of the strike confidence interval for shallow-dipping rock fracture sets on a sound foundation. We have made series of strike measurements at each of various dips on a wooden laboratory “fracture” to allow estimation of strike standard deviation at all shallow dips. This is accurate for smooth rock and serves as a lower bound estimate for rock with a rougher surface. These estimates are used, along with two independent iterative statistical procedures, to re-estimate the dip direction confidence interval of an actual set of nineteen fractures. The dip direction confidence interval is more than three times as wide as the corresponding dip confidence interval. This will be made visible in our transformed confidence region plot for the mean pole.     

回采工作面围岩采动应变场-渗透系数场耦合

采动岩体变形与渗透特性是工作面突水防治研究的基本问题。采用理论分析和数值模拟计算方法,分析了岩石变形-渗透特征及其三维定量关系,研究了煤炭开采中回采工作面围岩应变场、渗透系数场分布及其控制因素。研究结果表明,采动岩体渗透性变化主要取决于应变状态及应变增量,且随着垂直于裂隙的张应变的增加而增加;工作面后方垮落带和煤壁边缘的剪碎带产生剧烈的采动拉伸变形,渗透系数较采前显著增大,而支承压力区和整体移动带岩层产生较大的采动压缩变形,渗透系数较采前明显减小;工作面围岩垂向渗透系数较水平渗透系数增加的幅度及增加区的范围小,但对水体下采煤工作面涌水起主导作用。加快工作面推进速度、减小工作面斜长和采高能降低采场围岩渗透系数增加幅度,并将渗透系数场变化范围局限在采区附近,能有效减小采动对原始煤岩层渗透性的影响。      

Permeability tensor and representative elementary volume of fractured rock masses

Based on a simulation of three-dimensional fracture networks and a superposition principle of liquid dissipation energy for fractured rock masses, a model of the fracture permeability tensor is proposed. An elastic constitutive model of rock fractures, considering fracture closure and dilation during shearing, is also proposed, based on the dilation angle of the fracture. Algorithms of flow-path searching and calculation of the effective flow coefficients for fracture networks are presented, together with a discussion on the influence of geometric parameters of the fractures (trace length, spacing, aperture, orientation and the number of fracture sets) on magnitude, anisotropy of hydraulic permeability and the size of a representative elementary volume (REV). The anisotropy of hydraulic permeability of fractured rock masses is mainly affected by orientation and the number of fracture sets, and the REV size is mainly influenced by trace length, spacing and the number of fracture sets. The results of studies on REV size and the influence of in-situ stress on hydraulic conductivity of the rock mass on the slope of Jinping-I hydropower station, China, are presented using the developed models and methods. The simulation results agreed well with the results obtained from field water-pressure measurements, with an error of less than 10 %.     

逆断层控制构造裂缝发育的力学机制模拟

在探讨逆断层对构造裂缝控制作用的基础上, 应用Comsol有限元软件, 模拟不同水平作用力、不同岩性、不同断层倾角和距断层面远近等因素影响下断层及其周边区域的应力、应变情况, 并分析了不同控制因素下构造裂缝发育的规律。模拟结果表明, 构造裂缝的发育程度随施加的水平应力作用增大而线性增大; 岩石破裂前产生的应变量可用来描述岩石的脆性, 抗压强度对岩石破裂起主导作用, 而岩石的抗剪强度与岩石破裂发育有一定的相关性; 逆断层倾角存在一个临界角度, 使构造裂缝发育最为强烈; 逆断层控制裂缝发育, 距离断裂面越远裂缝发育程度越低, 并存在一个发育程度骤减的范围, 称为"断裂控制带", 该断裂控制带的形成和分布应该与断裂的性质、规模、断距及岩石力学参数等有关。      

基于连续-非连续方法的裂隙破坏与相互作用研究

岩体内裂隙等非连续结构面对岩体的强度及变形等力学特性有着显著的影响,研究岩体裂隙起裂、扩展、相互作用和贯通机制,对工程岩体力学行为的表征和工程性能的评价十分重要。本文基于连续介质力学模型的离散元方法,通过考虑裂隙分布、模型加载条件及其与裂隙产状的关系,建立了一系列裂隙力学计算模型,研究了不同模型裂隙扩展演化特征和岩体破裂机制,分析了岩体裂隙扩展规律及其对岩体破坏路径和强度的影响,研究结果表明:（1）裂隙岩体模型加载条件下的破坏起裂点、最终贯通破坏特征及损伤分布受控于裂隙的产状及其与最大主压应力取向角度大小及围压大小。（2）裂隙弱面走向与最大主压应力取向斜交时,裂隙弱面在加载条件下其端部裂隙扩展、贯通破坏表现比较明显,反之,当裂隙弱面走向与最大主压应力取向一致时,裂隙弱面被动影响裂隙模型内新生裂隙的萌生、扩展和贯通模式,自身未出现新的扩展破坏。（3）裂隙数目的增多和围压的增大会显著增加模型内部剪切裂缝的数量和模型破坏后的破碎程度,模型内部的损伤区域主要围绕破裂面呈滑移线型交叉分布,非破裂面区域损伤呈条带状X型分布。（4）裂隙弱面走向与最大主压应力取向斜交时,裂隙对岩体模型强度的弱化程度高于裂隙弱面走向与最大主压应力取向一致的情况,而裂隙模型破坏后的残余强度则正好相反。     

Fracturing around a preconditioned deep level gold mine stope

The mining faces of the highly stressed tabular stopes in South Africa's deep level gold mines are prone to a type of rock burst termed a face-burst. As a means of ameliorating these face-burst conditions, a destressing technique termed preconditioning has been employed. Part of the studies into the quantification of the effects and mechanisms of preconditioning was a detailed investigation of the fracture pattern around unpreconditioned and preconditioned stopes. Techniques included mapping of the fractures exposed in the mined areas of the stope, measurement of profiles of the stope hangingwall and the use of Ground Penetrating Radar (GPR). These studies indicate that preconditioning does not produce new sets of fractures and stress redistribution ahead of the face occurs by re-activation of specific pre-existing fractures, thereby reducing the potential for face-bursting.     

天然地应力作用下不同产状岩体裂隙渗流特性研究

为了研究天然地应力作用下裂隙产状等因素对深部岩体裂隙渗流特性的影响,基于单裂隙面渗透性服从负指数变化规律,建立了三维应力作用下不同产状裂隙的渗透系数计算公式,利用Lagrange乘子法分析裂隙面产状变化对其渗透性的影响,并分析了岩体裂隙有、无充填物对其渗透性的影响及敏感性;然后,以我国大陆地区地应力统计规律为例,分析了地表以下5 000 m范围内在天然地应力作用下裂隙渗透性随深度、产状的变化规律。结果表明：裂隙产状的变化对其渗流特性有明显影响,对于浅层岩体,在大主应力大致呈水平方向分布时,随着裂隙面倾角的增加,裂隙渗透系数逐渐降低;但随着深度的增加,在裂隙深度超过约200 m和裂隙面走向与大或中主应力方向大致一致时,裂隙渗透性反而会随着裂隙面倾角的增加逐渐增加,在裂隙面走向与小主应力方向垂直时增加最为明显;对于深部岩体,裂隙的渗透性很小,裂隙面产状的变化对其渗透性影响很弱;对于有充填物裂隙,岩块与充填物的弹模比和充填物泊松比的变化对裂隙渗透性的影响很小。研究结果可为深入研究我国深部岩体渗透特性变化规律提供借鉴意义。     

乌东德水电站右岸地下厂房随机块体特征研究

随机块体的几何特征和稳定性是水利水电工程中地下厂房支护设计的重要依据。对乌东德水电站右岸地下厂房随机块体的特征进行了深入研究。对研究区域内实测的裂隙参数进行了分组统计,确定了各组裂隙产状、迹长等参数的分布形式和大小。利用逆建模方法建立了三维裂隙网络模型,获取了各组裂隙的直径和三维密度。采用一般块体方法对地下厂房进行随机块体的识别和稳定性分析,利用GeneralBlock软件进行了10次随机实现,对10次计算结果进行了统计分析和讨论。研究结果表明,结构面与开挖面形成的随机块体集中在地下厂房的顶拱部位,10次随机实现中,地下厂房全长范围内平均每次形成的随机块体为414个;随机块体的平均体积为2.9 m³,最大块体的体积达152.0 m³;可移动块体中,大部分为稳定块体,不稳定块体均以坠落形式破坏;构成可移动块体的结构面多为3~4条,最多可达12条,其中倾W向的中等倾角裂隙是构成块体并可能造成块体失稳的最危险结构面;随机块体的平均深度为1.2 m,最大深度为8.8 m。建议支护锚杆应尽量穿透倾W向中等倾角的长大裂隙,且锚杆支护长度大于8.8 m。     

Origin of sheet structure, 1. Morphology and boundary conditions

Sheet structure, or large-scale exfoliation, is the division of a rock mass into lenses, plates or “sheets” approximately parallel to the earth's surface. Sheet fractures, which separate the plates, are characterized by surface markings resembling those formed during the brittle fracture of metals, glass and ceramics and those on joints in rock. Lineations, or hackle marks, parallel to the direction of fracture propagation, are common. Branching occurs during propagation so that the fractures appear as many echelon segments in profile. The opposing surfaces of sheet fractures observed in quarry walls and natural exposures are typically in contact. Damage, if any, to rock adjacent to sheet fractures is generally limited to a zone less than one centimeter wide. Sheets tend to parallel preferred orientations of microcracks in many rock masses. Most sheets that are exposed today are of prehistoric origin, but in some rock masses, such as Chelmsford granite in Massachusetts, new sheets continue to form.

Sheet structure forms in environments of high differential stress, dominated by large-scale compression parallel to an exposed rock surface. In parts of the Chelmsford granite the magnitude of this compression is greater than 30 MPa in the zone of sheeting. High differential stresses near the ground surface can result from several natural agents, including contemporary tectonic forces, vertical unloading of a rock mass that formed at depth under high triaxial compression, and suppression of expansion that would otherwise result from temperature increases or chemical alteration of the rock. Exfoliation of boulders during forest fires and spheroidal weathering of boulders appears to occur under similar states of stress, although the agents responsible for the stresses differ.     

Influence of three-dimensional roughness of rock fracture on seepage characteristics based on the digital image technology

Digital image processing technology can objectively reflect the surface roughness of coal and rock mass fractures and turn physical properties such as fluctuation height to physical data of the rock fracture surface to be used for numerical analyses; thus, it can effectively be applied to the seepage flow analysis of rock fractures. First, clear digital images of the fracture surface height are obtained under the same camera conditions, and then characteristic function values specific to various structure are normalized. Subsequently, the fluctuation height distribution of the rock surface is restored according to the measured fluctuation heights of local points. Finally, according to the surface fluctuation of fractures, invasive restructuring is carried out for three-dimensional fractures in the same coordinate system. Obtained physical parameters on the three-dimensional fractures of the rock are inputted into the COMSOL Multiphysics software to obtain the characteristic three-dimensional model of the rock fracture surface. The numerical analysis results are compared with experimental data on fracture seepage obtained from a seepage coupling true triaxial test system, and the relevance between the simulated result and the physical experiment is higher than 90%, which confirms that the integration of the digital image technology and the numerical analysis method can effectively simulate seepage in rough fractures. The hydraulic gradient of rough fractures and seepage velocity are also consistent with Forchheimer flow characteristics.     

隧洞控制测量中免棱镜全站仪测量方向与结构面夹角的限差分析

准确获取隧洞岩体结构面的几何信息对于隧洞的稳定性评价至关重要。研究区位于长春市净月潭吉林大学岩体力学隧洞实习基地。应用免棱镜全站仪进行岩体隧洞的控制测量,按照测量工作中"先控制后碎部"的基本原则,进行全站仪棱镜模式下的导线测量及在各导线点上极坐标法免棱镜模式的岩体结构面控制点测量;通过对同测站及异测站相同结构面控制点坐标数据的对比分析,探讨免棱镜全站仪测量方向与岩体结构面夹角的限差。研究表明:在隧洞控制测量中,免棱镜全站仪测量方向与结构面间夹角大于30°是准确获取岩体目标控制点坐标的必要条件之一。     

反倾岩坡倾倒变形结构面影响效应研究

为探究斜坡内赋存不同角度裂隙对反倾岩质斜坡倾倒变形影响效应,设计9组底摩擦试验,对比无裂隙、含一组陡裂隙和一陡一缓两组裂隙的3种类型试验,并研究不同裂隙角度对破裂面影响效应。研究发现:陡倾裂隙倾角的变化对岩坡变形及主破裂面形态有明显规律性影响,陡倾裂隙倾角越陡,发生倾倒破坏的初始破坏部位逐渐偏浅,破坏面积相对减小,破裂面从近直线型逐渐转变成近弧形;缓倾裂隙的赋存,使陡倾裂隙顺尖端开裂增长并且相互沟通形成贯通的破裂面过程更加快捷;主破裂面上覆层状岩体在倾倒变形时,会在其中部产生反向弯曲折断;通过定量分析,发现随着陡倾裂隙倾角增大,岩坡倾倒变形破坏幅度降低,而主破裂面无论是长度还是反映迹线复杂度的分维值亦随之降低。     

瀑布沟水电站库首右岸深部裂缝成因分析

水电工程常常建在高山峡谷地带,其天然岸坡通常由坡面向内有一个强卸荷带和弱卸荷带以及相应的强风化和弱风化带,内侧则为完整新鲜的岩石。对涉及的工程岸坡在正常卸荷带以内发育的一系列张性破裂或破裂带,称之为深拉裂缝。瀑布沟水电站库首右岸存在两个拉裂变形体,通过对其岸坡深部拉裂缝空间发育分布、变形特征的考察,综合分析造成深部裂缝发育规律与变形特征的因素。在此基础上提出,库首右岸深部拉裂缝是岸坡快速卸荷条件下浅表生改造的产物,其形成时期相当于河谷由宽谷深切为峡谷这一转换时期。     

改进遗传算法和支持向量机的岩体结构面聚类分析

岩体结构面控制着岩质边坡和地下洞室等岩体工程的稳定性,在岩体力学及水力学分析中起到关键作用。为对岩体结构面进行合理分组,精确地模拟岩体结构面网络的分布,提出一种融合改进遗传算法和支持向量机的聚类方法。首先,根据岩体结构面产状信息建立结构面分组的数学模型,采用改进的遗传算法计算结构面样本的全局最优聚类中心,再以聚类中心为训练样本,利用支持向量机方法将结构面样本进行完全划分。通过随机产生的结构面数据以及实际工程的运用表明,遗传-支持向量机聚类算法对岩体结构面的分组合理,获得的优势结构面结果可靠。     

基于RJM断续节理岩体强度与破坏特征的数值模拟研究

裂隙结构的存在对于工程岩体的强度和稳定性具有重要影响,岩石宏观裂隙的产生源自于微破裂的积累。针对岩体裂隙的粗糙特性,通过Matlab建立考虑粗糙度的节理模型（Roughness Joint Model）,采用简化的正弦曲线来表示粗糙节理,并将其导入到颗粒流试验模型中进行单轴压缩试验。对比完整岩体、直线型裂隙岩体、RJM岩体三者破坏的应力-应变曲线,改变裂隙倾角（与水平方向夹角）α,岩桥倾角β,裂隙密度γ,建立不同裂隙分布的断续节理岩体数值试样,开展一系列数值模拟试验。研究结果发现:（1）裂隙的存在明显降低了岩体的抗压强度,RJM模型峰值强度和峰值应变均高于直线型裂隙岩体;（2）岩体抗压强度总体上随裂隙倾角增大而增加,随裂隙密度增加而减小,但随岩桥倾角的改变呈非线性变化,岩桥倾角45°时峰值强度最低,峰值应变最小;（3）裂隙分布会影响岩体的破裂模式,微裂隙的扩展反映了岩体力学性质的各向异性;（4）不同倾角下增加裂隙密度,岩体强度下降程度不同,倾角75°时密度对强度影响最小,30°和60°时影响最大。