郭学飞，焦润成，曹 颖，闫 驰，王晟宇，赵丹凝

(北京市地质研究所，北京 100120)

0 引言

北京地区西部为太行山余脉，北部属燕山山脉，山区面积广泛，地质构造复杂，人类工程活动频繁，崩塌地质灾害较为发育[1]。为了有效防止崩塌灾害的发生，开展崩塌隐患的勘查和评估工作是十分重要的。传统的勘查方法属于接触式、单点调查[2]，但这仅限于工作人员能够到达的隐患点，对于很多所处位置高、地形复杂的崩塌隐患点，难以开展实地测量。这时倾斜摄影技术可以以其非接触、全视角、高精度测量的优势，在崩塌地质灾害调查监测中发挥重要的作用。

倾斜摄影是利用无人机搭载多台航摄仪，同时从多个侧面和垂直面获取目标区影像和位置信息的测量技术[3]。通过后期数据处理，可获取目标区高精度、可测量的三维实景模型，更加真实地反映出地物的原本面貌和形态特征。本文主要介绍利用倾斜摄影测量技术，对北京市一处典型崩塌隐患——中山村崩塌隐患进行了特征分析和稳定性评价，并进行了典型危岩体崩塌后落石运动轨迹的模拟分析。

1 研究区倾斜摄影测量

倾斜摄影测量主要包括外业无人机倾斜摄影获取数据和内业建立三维实景模型两部分[4-5]。本研究测量工作外业采用的飞行平台为旋翼电动无人机，数码相机选用五拼倾斜相机系统；为满足后期制作三维的需要，航向重叠及旁向重叠均设定为80%。内业采用大场景快速三维建模技术和模型精细化处理技术，通过多视立体空中三角测量、生成密集点云、构建TIN三角网、3D纹理映射等过程[6]，建立了研究区1∶500比例尺三维实景模型(图1)。该模型可实现陡坡危岩体的具体坐标、高程、几何尺寸等要素的量测。通过测量空间点坐标参数，采用空间三点坐标法计算可得到结构面产状要素[7]，从而得到其准确的几何特征、空间特征以及地质特征。据此，可进一步判定危岩体的稳定性、崩落岩石的运动轨迹和破坏能量，为崩塌地质灾害的防治提供参考依据。

图1 中山村崩塌隐患及其周边倾斜摄影三维实景模型Fig.1 Real-time 3D model of the potential collapse site

2 崩塌隐患基本特征

2.1 研究区地质背景

中山村崩塌隐患位于房山区中山村北侧的山体上。对倾斜摄影三维实景模型进行地形地貌分析和构造岩性推测，可得出如下推论：研究区在地貌上分为上陡下缓两个部分。其中，北侧山体高耸，坡体陡峭，覆盖层较薄，植被以灌木为主，且覆盖程度相对较低，说明其岩性性质相对强硬，风化剥蚀程度较低；同时，在该侧坡脚位置可见一系列明显的断层三角面，预示该处有断裂构造通过。而南侧坡体地势低矮，坡度较缓，覆盖层厚度大，植被覆盖程度高，说明其岩性性质相对软弱，受风化剥蚀程度高。

结合区域地质背景的先验知识，研究区发育有霞云岭逆冲推覆构造，展布方向与三维实景模型显示方位一致；断裂带南侧为石炭-二叠纪砂、板岩地层，北侧为寒武纪灰岩地层；受逆冲推覆作用，寒武纪灰岩地层叠覆于石炭-二叠纪砂、板岩地层之上(图2)。

图2 中山村崩塌隐患地质背景示意图Fig.2 Geological background sketch map of the potential collapse site

2.2 崩塌隐患基本特征

研究区崩塌隐患处于断裂构造形成的断层三角面上(图3)，坡体岩性为寒武纪灰岩，山坡相对高差为115 m，平均坡度为80°，坡体陡峭，具有发生崩塌落石的地貌条件。岩层产状为42°∠28°，倾向与坡体方向反向相交。受本身岩性和断裂构造作用的影响，山体发育以节理裂隙为主的多处结构面，采用空间三点坐标法计算，主要结构面产状如下：J1：165°∠82°；J2：262°∠87°；J3：33°∠62°。其中结构面J1倾斜方向与坡向基本一致，为近直立陡倾结构面，将岩体分成竖向的岩块；在其他结构面的共同作用下，岩体被切割为较破碎的块状。岩石下部经长时间的风化剥蚀后出现临空面，形成多处危岩体；崩塌体的类型主要为滑移式崩塌(1#、2#、3#危岩体)和拉裂式崩塌(4#、5#危岩体)，威胁山坡下方道路和中山村居民。各危岩体特征如表1所示。

图3 中山村崩塌隐患遥感判识结果示意图Fig.3 Interpretation results of the potential collapse site on remote sensing map

表1 危岩体特征统计Table 1 Statistics of dangerous rock mass characteristics

3 崩塌隐患稳定性分析

3.1 稳定性宏观分析

中山村崩塌隐患主要发育灰岩，岩体较完整，结构面结合一般，受节理裂隙影响发育多处危岩体。根据边坡岩体的完整性及坚硬程度，该边坡岩体类型属Ⅱ类，风化程度属微风化，危岩带平均高差约76 m。根据坡率法判定崩塌边坡稳定性[8]，中山村崩塌点容许坡率为1∶0.25～1∶0.35，即70.4°～75.6°，而坡体实际坡度平均约80°，超过容许坡率值，稳定性较差，容易发生随机块体崩落。

3.2 典型危岩体稳定性计算

1#危岩体位于边坡中下部，属于典型的滑移式崩塌体。在三维模型上测量，危岩体规模(长宽高)为16 m×4 m×4 m，呈不规则块状，体积约为260 m3，是研究区体积最大的一处危岩体。该危岩体下方具有临空面，且存在明显的后缘裂隙(图4)，裂隙最宽处达45 cm，破裂面的倾向与坡向基本一致，在暴雨或受到扰动时，该危岩体容易发生滑动。

图4 1#危岩体实地照片Fig.4 Picture of 1# dangerous rock mass

危岩体的稳定性受多种因素制约，其稳定性计算涉及较多参数和复杂的边界条件。本文参照《地质灾害防治工程勘察规范》(DB50/43—2018)中的滑移式崩塌体的稳定性计算公式和稳定性划分标准，对1#危岩体的稳定性进行评估[9]，计算公式如式(1)～(4)所示，稳定性划分标准如表2所示。

(1)

Q=ζeW

(2)

(3)

(4)

式中：FS——稳定性系数；

W——危岩体自重与建筑等地面荷载之和/(kN·m-1)；

α——滑面倾角/(°)。

Q——作用于滑坡体或其某条块的地震力/(kN·m-1)；

V——后缘裂隙水压力/(kN·m-1)；

U——滑面水压力/(kN·m-1)；

φ——滑面内摩擦角标准值/(°)；

c——滑面黏聚力标准值/kPa；

l——滑面长度/m；

ζe——地震水平系数；

γw——水的重度/(kN·m-3)；

hw——裂隙充水高度/m。

表2 稳定状态划分标准Table 2 Standard of division for rock stability

注：Fst为滑坡稳定性安全系数。

对上述危岩体分现状、暴雨和地震三种工况进行稳定性计算(计算参数取值如表3所示)，得到三种工况下危岩体的稳定性系数：现状工况下，稳定性系数为1.05，属于欠稳定状态；暴雨工况下，稳定性系数为0.88，属于不稳定状态；地震工况下，稳定性系数为0.73，属于不稳定状态。

根据上述定性评价和定量计算结果，综合分析可知：中山村崩塌隐患点在暴雨和地震条件下表现为不稳定，容易发生随机崩塌落石。

4 典型危岩体崩落轨迹模拟

危岩崩塌是陡峭边坡上的岩块在重力及其他外力作用下突然脱离母岩，并以坠落、滑移、滚动、弹跳或其组合方式向坡下猛烈运动，最终堆积在坡脚的地质灾害；具有突发性强，影响因素多，随机性高，频率大等特点[10-12]。危岩体一旦由于某种原因而崩落，将直接威胁其下方的道路、居民点等建筑物，危害行人和居民的生命财产安全。因此，对危岩体运动破坏的路径进行分析有重要的意义。

表3 1#危岩体稳定性计算参数取值一览表Table 3 Parameter values of stability calculation for 1# dangerous rock mass

4.1 模拟原理介绍

本文利用RocFall软件对1#典型危岩体进行落石运动轨迹模拟分析，其模拟原理为：边坡上的危岩体相对坡体下方具有较大的势能，当危岩体失稳崩落时，危岩体势能转换为动能，开始向下坠落，当落石接触到坡面时开始反弹，由于不同的地表性质具有不同的切向和法向阻尼系数，落石的弹跳高度也不同；同时，接触坡面对落石的动能有消减作用，直至落石的动能为零，整个过程遵循能量转化和守恒定律。据此原理，RocFall软件将落石视为均匀且不会破碎的质点，坡面为各向同性的弹性塑体，忽略空气阻力，通过对危岩体崩落过程所接触的不同的坡面性质设置相应的阻尼系数，可模拟落石在崩落过程中的能量、弹跳高度与运动轨迹等信息。

4.2 危岩体崩落轨迹模拟

首先明确1#危岩体崩落时可能经过的地形剖面，利用高精度三维实景模型，可获取精细的1∶500比例尺地形剖面线，结合实地调查验证，对地形剖面各坡段的地表性质进行判断，以确定各坡段相应的阻尼系数；然后根据RocFall软件的要求，将地形剖面线比例尺转换为1∶1 000，并导入软件中，设定好落石位置和相应参数后，即可进行落石轨迹模拟。

本次模拟对象为1#危岩体，危岩体体积为260 m3，初始速度为0；从坡顶到坡底，坡面的岩土体性质依次为：裸露的基岩，少量的灌木覆盖的风化岩面，少量乔木覆盖的密实边坡，碎石土为主的坝阶地，混凝土路面，碎石土、块石为主的坝阶地，灌木覆盖的土质边坡；各坡段切向阻尼系数Rt和法向阻尼系数Rn取值如表4所示。对1#危岩体进行50次崩落路径模拟的结果如图5所示。

模拟结果显示：危岩体失稳掉落后，兼具滑动和弹跳两种运动模式。落石在贴坡面滑动水平距离约30 m时，能量迅速增加，开始多次弹跳；落石经翻滚、跳跃、碰撞能量逐渐被地面消减，在坡脚处能量达到最低，落石逐渐堆积。落石最大水平运动距离可达143 m，超越了坡下道路和房屋安全距离，将对行人和居民的人身财产安全造成危害；落石最大跳跃高度可达6.4 m，出现在中山村上方的道路附近，相应的落石能量在该处也第一次达到最高值，因此在考虑使用拦挡方式进行灾害防治时，宜将拦挡工程建在上游。

表4 坡段特征与切向/法向阻尼系数取值一览表Table 4 Slope features and tangential/normal damping coefficient values

图5 中山村崩塌隐患1#典型危岩体崩落模拟结果Fig.5 Simulated movement track of 1# dangerous rock mass

5 结论

(1)大比例尺倾斜摄影能够获取工作区高精度的三维实景模型，可实现地质灾害发育背景、灾害体基本特征、威胁对象等关键信息的识别，在崩塌灾害调查评价工作中能够发挥重要的作用。

(2)倾斜摄影测量技术能够获取崩塌危岩体的位置、高程、几何形态等信息，进而可以判断危岩体的稳定性，模拟崩塌落石的运动轨迹，为崩塌灾害的防治提供有效技术支撑。

(3)倾斜摄影测量技术具有非接触、全视角、高精度的特点，能够提高崩塌灾害隐患的调查效率，保障调查人员安全，在高位崩塌灾害隐患调查中具有独特的优势。