四川罗家青杠岭崩塌风险的定量评价研究

我国多山,崩塌灾害频繁发生,相应的风险评价也得到了越来越多的关注。由于崩塌发生和运移过程的高度不确定性以及历史数据的不完备,往往很难进行相应的定量风险评价。四川罗家青杠岭的崩塌现场非常典型,而且现场资料较全、历史数据较多并且明确,是开展崩塌风险定量研究的很好实例。通过现场工程地质调查、测绘和统计分析,确定了历史崩塌的物源区、堆积区、最大运移距离、年发生概率以及坡体上的4块典型危岩体A-D。基于历史崩塌堆积区的块石统计特征以及物源区危岩体失稳启动位置的不确定性,利用二维Rockfall模拟软件对所在坡面的恢复系数及摩擦系数进行了反演。在此基础上,对危岩体A-D失稳后的运动特征进行了随机性数值模拟和统计分析,从而确定了崩塌的到达概率。基于崩塌发生概率、到达概率、承灾体时空分布概率和易损性的乘积,作者对罗家青杠岭崩塌进行了定量风险评价。评价结果表明,危岩体A和D的风险值处于不可接受的风险区间,块石B和C的风险值处于警告的风险区间,严重威胁着坡脚附近居民的生命财产安全,有必要采取相应的防灾减灾措施。     

Probabilistic key block analysis of a mine ventilation shaft stability – a case study

In this study, the probabilistic key block analysis was applied to evaluate the stability of a mine ventilation shaft developed in a rock mass of granite. The key blocks were identified based on the block theory. The variations of discontinuity orientations were fitted with the Beta distribution and taken into consideration. The key block forming probabilities were analyzed. For simplification of calculations the first-order second-moment (FOSM) approximation was employed for probability estimation. With the considerations of the rock properties as random variables and applications of several statistical analysis tools, the key block failure probabilities, the probabilistic distribution of safety factors, and the probabilistic distribution of potential maximum key block volumes were analyzed. The analysis indicated that although the safety factor calculated based on the mean values of the variables was above 1.0 for the stability of the most critical key block, the block had a considerable probability of failure with a significant rock volume due to variations in discontinuity orientations and rock properties. Without promptly applying supports to the rock excavation, the shaft had a significant likelihood of failure.     

输电线路杆塔危岩威胁等级评价方法研究

从影响滚石击中塔体时的能量以及塔体被滚石击中的可能性出发,提出现场评价危岩对塔体威胁等级的方法。该方法考虑了边坡形状、塔位与山体的相对位置、坡上危岩的数量和大小以及坡脚滚石的数量和大小等因素。当边坡总体坡度相同时,上缓下陡的折线边坡,滚石落地速度最大;而内凹的曲线边坡,滚石落地速度最小。位于山体内凹一侧的塔位,塔体受危岩威胁最大,而位于山体凸出一侧的塔位,受危岩威胁最小。     

节理粗糙度对岩块体积的影响规律研究

岩体结构直接控制着岩体力学、水力学特性,岩块体积分布规律可直接反映岩体结构特征。目前研究多将节理面简化为平面,忽略了粗糙度对岩块体积的影响。利用Hurst指数H和高度均方根R_q表征节理粗糙度,研究了节理粗糙度对岩块数量和体积分布规律的影响。结果表明：(1)H、R_q控制节理的粗糙起伏形态,粗糙度随着R_q的增大而增大、随着H的增大而减小;(2)块体的数量总体上与节理粗糙度正相关,即随着粗糙度的增大而增加;(3)R_q和H对岩块体积分布的影响,主要通过改变小体积块体所占比例实现,块体的平均体积和中位体积随着粗糙度的增大而减小;(4)当节理正交且间距相近时,块体数量随粗糙度的变化可以分为3个阶段,即稳定阶段、初始增长阶段和快速增长阶段,R_q和H共同控制着各区间的相对分布范围。最后,基于摄影测量采集的数据建立了岩体节理模型,对大连某公路边坡的岩块体积分布进行了研究,验证了上述结论的准确性。     

A fractal fragmentation model for rockfalls

The impact-induced rock mass fragmentation in a rockfall is analyzed by comparing the in situ block size distribution (IBSD) of the rock mass detached from the cliff face and the resultant rockfall block size distribution (RBSD) of the rockfall fragments on the slope. The analysis of several inventoried rockfall events suggests that the volumes of the rockfall fragments can be characterized by a power law distribution. We propose the application of a three-parameter rockfall fractal fragmentation model (RFFM) for the transformation of the IBSD into the RBSD. A discrete fracture network model is used to simulate the discontinuity pattern of the detached rock mass and to generate the IBSD. Each block of the IBSD of the detached rock mass is an initiator. A survival rate is included to express the proportion of the unbroken blocks after the impact on the ground surface. The model was calibrated using the volume distribution of a rockfall event in Vilanova de Banat in the Cadí Sierra, Eastern Pyrenees, Spain. The RBSD was obtained directly in the field, by measuring the rock block fragments deposited on the slope. The IBSD and the RBSD were fitted by exponential and power law functions, respectively. The results show that the proposed fractal model can successfully generate the RBSD from the IBSD and indicate the model parameter values for the case study.     

某电站地下厂房裂隙发育规律及其对围岩稳定性的影响

裂隙岩体内的洞室围岩块体稳定问题是抽水蓄能电站地下厂房建设中要解决的重要工程地质问题。为了在厂房开挖之前较准确地预测不稳定块体的发育位置,本文以某抽水蓄能电站为例,对厂房探洞内的裂隙进行统计分析,结合厂房区构造断层的发育情况,确定了裂隙发育规律;根据裂隙发育规律将厂房区岩体划分为A、B、C、D、E 5个构造分区;并利用Unwedge程序,以各分区裂隙组合为基础,确定可能形成的块体位置,评价了其稳定性状况。研究结果表明,裂隙受多期构造控制,在一定范围内具有方向性,A、E区优势裂隙与构造产状接近,B、C、D区优势裂隙与构造产状斜交;裂隙在厂房上、下游侧边墙各形成6处潜在不稳定块体,顶拱处形成10处潜在不稳定块体,右端墙处形成1处潜在不稳定块体。     

三峡库区望霞危岩体孕灾地质环境条件分析

危岩体是一种典型的灾害地质体,其形成、变形与破坏受控于地质环境条件,地质基础要素在孕灾的过程中起着决定性作用,各基本地质要素与灾害形成的关系是研究危岩体变形规律和破坏机理的关键性问题。通过对三峡库区望霞危岩体区域地质构造作用分析,横石溪背斜构造解析,各类地质结构面发育特征统计和危岩体岩性组合方式及工程地质条件剖析,结合危岩体两次变形破坏过程,分析各因素对危岩体的形成及变形的关系。研究认为:孕育望霞危岩体的横石溪背斜雏形形成于印支期的隆升阶段,至燕山主期形成褶皱主体,后期则以继承性改造为主,危岩体的形成及孕育主要集中在背斜核部区域,枢纽位置尤为发育。控制危岩体变形的5条关键裂缝与顶部发育的"V"字形裂缝关系密切,其是在纵弯褶皱作用下背斜顶部脆性岩层受到强烈的拉裂影响,叠加台阶式层间滑动作用,形成与轴面近于平行的拉张裂缝。在望霞危岩体中发育多组地质结构面,其中原生结构面对各条关键裂缝及关键块体提供基础背景值,控制其规模、范围及变形强度。控制着变形关键部位特别是前缘软弱层变形位置、牵引区段及强烈程度。构造结构面主要控制关键裂缝的延伸方向、距离及样式,控制关键块体的变形方向,裂解范围及规模,掉块及崩裂位置,控制危岩体底部关键部位即压张裂缝的形成位置、张开度及位移方向。次生结构面主要控制后缘变形区,对各条关键裂缝控制其继承性变形,控制关键块体扩张,对岩体裂解破坏至关重要。同时,望霞危岩体在地层岩性上构成了典型的硬—软—硬岩性组合,为危岩体的变形提供了地质结构组合条件,其底部大面积的采空区控制着危岩体整体滑塌变形。     

差异风化型危岩力学模型及破坏机制研究

差异风化型危岩多形成于由砂岩和泥岩等软、硬岩性组成的互层陡倾岩质边坡中,其崩塌破坏属地质灾害中的常见形态。以普洱渡危岩高边坡为例,推导出差异风化型危岩岩腔后壁泥岩压应力随岩腔深度增加而增大的计算公式。据普洱渡危岩高边坡地形地貌和地质构造特征建立了两类4种差异风化型危岩力学模型：第1类为主控结构面贯通率等于1时的泥岩基座压裂破坏型和转动破坏型危岩;第2类为主控结构面贯通率小于1时的坠落型和倾倒型危岩。利用岩石强度理论推导出研究区砂岩的二次抛物线型摩尔强度包络方程,以极限平衡理论和摩尔强度理论推导了4种危岩体在自重、地震力和裂隙水压力共同作用下的崩塌破坏判别表达式,并由此反演出危岩体的临界崩塌边界方程,得出危岩体崩塌边界（岩腔深度、危岩体厚度和高度、主控裂隙深度）之间的关系,为现场预测和判别危岩体的稳定性和崩落时序提供简便可靠的依据。     

Evaluation of rock falls in an urban area: the case of Boğaziçi (Erzincan/Turkey)

Disasters caused by events such as earthquake, flooding, rock falls, landslides are often encountered. However, generally, the reasons for the destructive and devastating effects of these nature events are that settlement locations were chosen without site investigation studies, or that available studies were inadequate. Such inadequacies in the field are related to inappropriate settlement location and the resulting damage caused by rock falls. This study evaluated rockfall risk in a settlement that developed in a similar manner. The study was carried out in Bo?aziçi village of Kemah (Erzincan/Turkey), which is located in a very important tectonic zone. The study site is located on the lower sections of an area with very steep cliffs and 50–75° slopes. This cliff, which is the source of rockfalls, has a slope dip of approximately 90°. The cliff comprises 25–30 m high, fractured and cracked basaltic volcanic mass. To determine block size in the study area, scanline survey measurements and block size measurements were performed on blocks that loosened and fell from the cliff face. It was found that block sizes reached 6 m³. Rockfall analyses were performed along the selected profiles using the Rockfall V.4.0 software. Kinetic energy, bounce height, horizontal location of rock end-points, and velocity of the rocks along each section were evaluated separately for each profile. This data were used to produce distribution maps for each profile and the settlement was evaluated in terms of rockfall risk. The results indicate that the study area was at risk of future rockfalls and that it would be appropriate to relocate one part of the settlement.     

Stochastic Representation of Sedimentary Geology

Discrete fracture network representations of discontinuities in rock masses have been shown to be useful in capturing heterogeneity in rock mass properties. Providing computational efficiency in the resulting simulations and analyses is attained, these fracture representations can be combined with structural modelling and sampling algorithms. Multiple fracture network realisations can be generated and the resulting rock mass properties interrogated. Statistical analyses based on fracture connectivity, block size distribution and slope stability can be performed and provide results defined in terms of confidence intervals. For sedimentary geology consisting of dense bedding, equivalent medium continuum methods have traditionally been used in preference to discrete fracture representations due to the large numbers of structures involved and resulting computational complexity. In this paper, it is shown that stochastic representation of these layers can be employed. An analytical solution to accommodate bedding given an assumed block size distribution has been derived. Using this formulation, polyhedral modelling has been used to investigate the influence of bedding on block formation and block size distributions using field data. It is shown that the analysis is both computationally efficient and can capture truncation of size distribution by such layers without numerical methods.     

Evaluating Kinematic Controls on Planar Translational Slope Failure Mechanisms Using Three-Dimensional Distinct Element Modelling

This paper investigates the importance of kinematic release mechanisms on planar translational slope failure using three-dimensional distinct element codes. The importance of the dip and dip direction of the rear, basal and lateral release surfaces and their influence on failure mechanism, dilation, and the development of step-path failures is illustrated. The three-dimensional block shape and volume of the unstable rock masses simulated with the different discontinuity set geometries are characterized. Two assumed three-dimensional slope models are investigated in order to assess the importance of varying kinematic confinement/release mechanisms. These two assumed boundary conditions are shown to be critical in the development of asymmetrical rock mass deformation patterns. Scale effects due to the block size and discontinuity persistence are shown to control the calculated displacement and failure mechanisms. The numerical modelling results are also demonstrated to be sensitive to the assumed normal and shear stiffness of the discontinuities. The influence of the factors investigated on the failure of a single rock block versus a rock mass are compared and discussed.     

岩体结构稳定分析原理和方法

工程实践表明,在边坡、地基或地下工程岩体中几组地质结构面和临空面的组合往往构成可能失稳的岩石结构体、存在向临空面运动的可能性、影响到工程岩体的稳定性。这时必须针对这种可能失稳的结构体进行分析,采取必要的加固或处理。关于岩石结构体的稳定分析,文献中已有一些报导,但偏重于图解和计算方法、一般是针对特定单一块体的解析,仅应用于边坡工程。     

Numerical Analysis of the 650,000 m<Superscript>2</Superscript> Åknes Rock Slope based on Measured Displacements and Geotechnical Data

The unstable 650,000 m² Åknes rock slope (Western Norway) poses a hazard, as a sudden failure may cause a destructive tsunami in the fjord. In this study the slope was divided into blocks based on displacements measured at the slope surface. Discontinuous deformation analysis (DDA) showed that three or four blocks in the upper half may be considered as potential subareas that may fail catastrophically. The lower half may be divided into two or three blocks, but more limited data introduces more uncertainty into block definition. The Universal Distinct Element code (UDEC) was used for two-dimensional (2D) stability analyses. By varying fracture geometry, fracture friction, and groundwater conditions within reasonable limits based on site-specific data a number of possible models were compared. The conclusions show that models that were unstable to great depths were in closer agreement with shear strength parameters derived from an earlier study than models that were unstable to smaller depths. The length (depth) of the outcropping fracture, along which shear displacements are shown to occur, plays an important role. A (shallow) slide at 30 m, in which displacements have been documented by borehole measurements, will reduce the stability at greater depths. Increased groundwater pressure is demonstrated to be less critical for very deep slope instability. The results of the DDA and UDEC modelling will be useful for planning of future investigations, interpretation of the subsequent results, further development of the early warning system and in the tsunami modelling.     

Rock slopes risk assessment based on advanced geostructural survey techniques

The rock mass structure determines the possible unstable blocks that can induce rock fall phenomena. The stability analyses must therefore be based on an accurate geo-structural survey. In this work, the stability conditions of several steep slopes along a motorway in the Far East have been evaluated through key block analysis based on traditional surveys and on laser scanner acquisitions. Discontinuity orientations and positions on the rock face are derived from the point cloud in order to perform the reconstruction of the rock mass and to identify blocks in the slope. Results obtained from both the traditional and the new method is in good agreement. Stability analyses have been performed for evaluating the kinematic feasibility of different failure mechanisms. The rock block shapes and volumes are computed by performing 2D and 3D analyses whereas the failure mechanisms are examined using the key block method. Parametrical analyses have been carried on to evaluate the influence of slope angle variation. DEM models have also been set up. The relative hazard is determined by statistically evaluating the kinematical feasibility of different failure mechanisms. Hazard mapping has been utilized to identify the best methodology for risk mitigation.     

裂隙岩体块体化程度评价方法的若干问题修正

块体化程度是评价岩体完整性的一种新指标,能从三维角度表征岩体破碎程度,但目前该方法仍存在未充分考虑岩体切割程度及块体规模、未限定基础应用条件、块体化程度等级划分不合理等缺陷。针对上述缺陷,深入分析了其产生原因,并借鉴岩体切割程度、三维块度模数、体积RQD等计算原理,限定了块体体积百分比相关概念,提出了考虑岩体完整性块度尺寸效应的块体体积综合百分比计算方法,确立了块体化程度等级及分级指标取值依据,建立了修正的块体化程度评价方法。通过对比块体化程度评价修正方法、岩土工程规范对岩体完整性的划分结果,分析了修正方法的合理性。分别以广西铜坑矿锌多金属矿体、乌东德水电工程的块体研究数据为基础,开展了修正块体化程度评价方法的实例验证,结果表明：利用修正块体化程度评价方法计算的铜坑矿+255 m中段4#试验区岩体块体体积百分比为11.18%,属轻度块状化岩体;乌东德水电站PD49-1平硐、PD4支硐块体体积百分比分别为12.847%、10.168%,均属于轻度块状化岩体,岩体完整性程度为较完整级别。与现有块体化程度评价方法比较,修正方法计算的块体体积百分比能够更准确地从三维角度表征真实岩体的完整性。研究成果对精确刻画岩体三维完整性具有重要意义。     

Use of backpropagation neural network for landslide monitoring: a case study in the higher Himalaya

Static and dynamic rock slope stability analyses were performed using a numerical discontinuum modelling technique for a 700-m high rock slope in western Norway. The rock slope has been investigated by the Geological Survey of Norway (NGU), which has been carrying out rock slide studies for the county Møre and Romsdal in western Norway. The purpose of numerical modelling was to estimate the volume of the rock mass that could potentially slide under static and dynamic forces. This estimation was required to assess the run-up heights (tsunami) in a fjord that could potentially be caused by the rockslide. Three cases have been simulated for predicting the behaviour of the rock slope. First, an initial static loading is applied in the numerical model to simulate the prevailing rock mass conditions at the site. Second, saturated and weathered joint conditions are modelled by reducing the residual friction angle along the discontinuities of the rock mass. In doing so, the model simulates the effect of degradation of discontinuities in the rock slope. Third, a dynamic loading, based on peak ground accelerations expected in the area, is applied to simulate dynamic earthquake conditions.

These numerical studies have provided some useful insights into the deformation mechanisms in the rock slope. Both sliding and rotation of blocks start to occur once the residual friction angle along the discontinuities is reduced and when the region is shaken by a strong earthquake. The results indicate that, due to variations in the inclination of discontinuities, the entire slope does not become unstable and that down-slope sliding and rotation of blocks occur mainly on the top layers of the slope. Within the range of parameter values considered for this study, it is unlikely that the whole rock slope can be destabilised. The study provides an illustration of how the geo-mechanical properties of a rock mass can be integrated in a discontinuum rock slope model, which is used for predicting the behaviour of the slope under existing environmental and earthquake conditions. This model has helped not only to better understand the dynamics of the rockslide but also to estimate the potential rock volume that can become unstable when subjected to static and dynamic loads.     

基于结构面网络模拟的桥基岩体空间状态分析

分析结构面在岩体中的分布和组合特征是研究岩体稳定性的前提。根据结构面网络随机模拟方法,通过对岩体表面结构面的测量,可以建立岩体结构面空间分布的概率统计模型,从而确立了岩体结构面的网络结构模型。在此基础上分析了局部临空块体的空间几何形态,指出了块体的可能破坏形式,从宏观角度对块体系统整体稳定性进行了预测。     

大型地下厂房块体稳定性简便分析方法

块体是结构控制型岩体中常见的潜在危险源之一。利用极限平衡法及强度折减法两种方法计算了某在建特大型水电站地下厂房开挖揭露的部分块体的安全系数,并根据计算结果提出一种利用块体几何及力学参数判断其稳定性的简便图解方法,经现场监测数据验证计算结果可靠性可满足工程要求。研究表明,对同一块体而言,极限平衡法和强度折减法得到的安全系数以及对其稳定性的总体判别结果并不一致。强度折减法受软件算法及网格尺寸影响,结果偏于保守。简单块体的安全系数计算应以极限平衡法为主,而复杂形态块体的安全系数用强度折减法计算较为方便。利用垂向地应力、块体体积、最大角点深度及结构面等效强度等4个指标并结合块体稳定性判别分区图,可满足快速判断块体稳定性的需要。对于判别为不稳定的块体,应及时支护并考虑加强支护。研究成果可用于类似工程块体稳定性的快速分析。     

Rockfall hazard analysis for an historical Castle in Kastamonu (Turkey)

The Kastamonu Castle located on a sandstone hill with Eocene age is one of the most historical and touristic places in Kastamonu city center. The settlement of the city expanded towards the hill of the Kastamonu Castle and adversely affected by rockfalls in the past. The rockfall problems around the castle could be related to jointing, weathering, freezing-thawing and earthquake effects or a combination. In this study, the rockfall hazard at the castle is evaluated by two-dimensional rockfall analyses along 17 profiles selected in different orientations. Different size of rock blocks and various types of movements are taken into consideration in the analyses. Fall-out distance, bounce height, kinetic energy and velocity of the sandstone blocks are separately evaluated. The obtained data are used to define the possible rockfall hazard zones. Finally, the areas having potential rockfall risks are distinguished. Based on the evaluation of the data, rock bolting after removing of unstable blocks and supporting the area with the protective fences are suggested.     

Preventing and analysis of falling rocks: A case of Sarica village (Gürün,Turkey)

Rock fall is a common natural hazard causing significant damage to infrastructure and loss of life and property. This paper describes selection of prevention method, construction processes and the results of a rockfall hazard assessment for the village of Sarica (Gürün region, Turkey) based on an engineering-geological model. The study area is located on the lower parts of an area with very steep cliffs and 40-50⁰ slopes with fractured sedimentary and volcanic rocks. Detailed studies have been conducted in the region to prevent the risk of any rockfalls. Maps and charts are prepared including the rock blocks that are likely to fall and the related housing area under the risk. It was found that block sizes reached 3 m³. Rockfall analyses were performed along the selected profiles using the Rocfall V.5.0 software. Kinetic energy, bounce height, horizontal location of rock end-points and velocity of the rocks along each section were evaluated separately for each profile. This data was used to select design and construction of rockfall preventing system. In this study professional type rock fall barrier (catchment fence) was used. This method was chosen because it was cheaper, more practical and reliable than others.