1 引言

在黄土高原地区，黄土滑坡是一种常见的地质灾害，其数量多、分布广、规模小、监测难度大、危害严重，每年都会造成一定程度的生命财产损失（张茂省等，2006；张树轩，2017）。滑坡危险性评价是减灾防灾的重要措施之一，近年来已成为滑坡研究的热点（王云龙等，2011；胡瑞林等，2013）。随着无人机技术、GIS技术、数理统计分析方法、滑坡稳定性建模技术的发展和推广应用，滑坡危险性评价在技术方法上已取得了很大的突破（Guzzetti et al., 2006；石菊松等，2009；Pradhan et al., 2010；Jia et al., 2012；彭令等，2013；Park et al., 2013；Martha et al., 2013；Grant et al., 2016；赵良军等，2017）。近年来，随着山区城镇建设规模的扩张，造成重大人员伤亡和经济财产损失的滑坡灾害频发。特别是黄土高原山区城镇，由于人类开挖坡脚建窑建房等活动，导致城镇周边斜坡变得高陡，并使斜坡原有的应力平衡状态被打破，在降雨条件下极易发生黄土滑坡灾害，对城镇周边居民构成了严重的威胁，山区城镇已成为滑坡防治工作的重点地区（薛强等，2019）。以往的滑坡危险性评价大多是在县域（几千平方千米）尺度上开展，精度较低，难以满足新时代城镇化建设对滑坡灾害防治工作提出的新要求。如何进行县城（几十平方千米）尺度上的大比例尺滑坡危险性评价，是山区城镇滑坡灾害防治的关键（张茂省等，2019）。

斜坡是滑坡灾害发生的基本地形地貌单元，即滑坡发生的控制因素，以斜坡为单元开展山区城镇大比例尺滑坡危险性评价，能够较准确地预测滑坡发生的具体空间位置（薛强等，2015）。斜坡单元主要是根据河流沟谷的发育情况和地形地貌的完整性划分，主要划分方法有基于GIS的斜坡单元划分法和野外调查划分法（谷天峰等，2013；颜阁等，2017）。同时，滑坡的发生还受降雨等诱发因素的影响（黄润秋，2007；张茂省等，2011；殷跃平等，2017），降雨诱发滑坡发生主要是通过斜坡土体含水率的改变而影响，即降雨过后斜坡土体含水率升高，土体强度指标（黏聚力c和内摩擦角φ）降低，导致黄土斜坡的变形和失稳（常金源等，2015；薛强等, 2014, 2018）。当斜坡的坡度、坡高、坡型不同时，其发生黄土滑坡的概率也不同，而降雨诱发的黄土滑坡则多发生在高陡斜坡地带（庄建琦等，2013）。因此，在一定的坡度、坡高、坡型情况下，当土体含水率达到什么程度时斜坡会变形失稳是进行滑坡危险性评价的关键。笔者在野外调查的基础上，结合黄土斜坡土体含水率监测结果，应用斜坡稳定性评价和GIS技术，以斜坡为评价单元，开展基于土体含水率的陕北黄土高原地区典型山区城镇（清涧县城区）黄土滑坡危险性评价，为城镇地质灾害防治和土地利用规划提供技术支撑。

2 研究区概况

清涧县位于陕西省北部，榆林市东南缘，黄河中游、无定河下游，地处陕北黄土高原腹地。城区为黄土梁峁沟壑地貌，地形切割强烈，清涧河自北向南从清涧县城区流过。就地貌类型划分，清涧县城区及周边可分为河谷阶地区和黄土梁峁区2种地貌类型，河谷阶地区主要分布在县城境内的清涧河沿岸，黄土梁峁区在清涧县城区及周边大部分地区均有分布，为区内最主要的地貌类型。城区出露的地层岩性主要为晚更新世（Qp³）和中更新世（Qp²）黄土，下伏三叠系上统胡家村组（T₃h）砂泥岩，在河谷地带分布有全新世冲积层（Qh^al）。清涧县多年平均降水量为494.06 mm，年内降水量分配不均，在7—9月份多以雷阵雨或暴雨形式降落，占全年总降水量的61.5%。

受地质环境条件的影响，清涧县城区滑坡等地质灾害多发。城区范围内滑坡灾害主要发育在清涧河河谷两侧边坡上和黄土梁峁区的陡坡地段，稳定性差或较差，均为黄土滑坡，且以小型浅层滑坡为主。滑坡形态特征较明显，易识别，后壁平面形态多呈典型的圈椅状，坡度在20°~65°。降雨、人工削坡取土等是滑坡的主要诱发因素。通过清涧县城区地质灾害调查，清涧县城区周边斜坡地带共发育黄土滑坡138处（图 1），以浅层滑坡为主，广泛分布于黄土梁峁侧壁的沟壑地带。

3 城区斜坡单元划分与地形参数提取 3.1 斜坡单元划分与滑坡密度计算

斜坡是滑坡发生的基本地形地貌单元，因此斜坡单元是进行滑坡危险性评价的理想单元，也是跟现场实际情况较为符合的（薛强等，2015）。本研究的基础数据来源于清涧县城区地质灾害野外调查，调查范围为清涧县城及周边区域的斜坡地带，面积31.4 km²。按照河流沟谷的发育情况和地形地貌的完整性，将清涧县城及周边区域的斜坡地带共划分为925个斜坡单元，其中单元面积最小为0.48×10⁴ m²，最大为14.59×10⁴ m²，平均为3.40×10⁴ m²（图 2）。

在滑坡野外调查及填图的基础上，采用滑坡体的面积指标，利用ArcGIS平台，计算清涧县城及周边区域斜坡单元内已发生滑坡的面密度，即滑坡体面积占所处斜坡单元面积的比例，单位：m²/m²（图 3）。

3.2 斜坡单元地形参数提取

清涧县城区及周边地形破碎，可形成滑坡的斜坡单元发育，所有斜坡单元均覆盖有较厚的黄土层，黄土滑坡即发生在黄土层内。黄土斜坡在不同坡度、坡高、坡型情况下，发生滑坡的概率也有所不同，因此需首先提取斜坡单元的基本形态参数。

坡度提取：利用ArcGIS从清涧县城区及周边1:1万DEM数据中提取栅格坡度信息，然后在斜坡单元内求取栅格坡度平均值作为该斜坡体坡度指标（图 4）。

坡高提取：斜坡高度（坡高）可利用地形起伏度来间接获取，地形起伏度是指在一个特定的区域内，最高点高程与最低点高程的差值。从清涧县城区及周边1:1万DEM数据中按5×5单元格（25 m×25 m）间距计算地形起伏度，然后在斜坡单元内求取地形起伏度平均值作为该斜坡体坡高指标（图 5）。

坡型提取：坡型是指地表坡面的形态，可分为直型坡、凸型坡和凹型坡3种基本类型，在GIS中坡型可以利用地表的曲率P进行描述和量化。当地表曲率P>0时，坡型为凸型坡；当地表曲率P=0时，坡型为直型坡；当地表曲率P < 0时，坡型为凹型坡。在实际地形中，地表曲率P=0的平直坡极少，为了更加直观地反映研究区坡型情况，经过计算和现场调查核实，当-0.5 < P < 0.5时，斜坡总体为直线型，且该取值范围内的坡型对地质灾害的贡献一致。因此做如下定义：当P≥0.5时，坡型为凸型坡；当-0.5 < P < 0.5时，坡型为直型坡；当P≤-0.5时，坡型为凹型坡。按照该取值范围定义的坡型与研究区实际情况一致（图 6）。

4 城区斜坡危险性评价

本文所指危险性是指黄土斜坡在一定降雨入渗引发土体含水率增大的情况下，发生黄土滑坡的概率。由于黄土滑坡发生在黄土层内部，其影响因素主要有两类：一类为内在控制因素，主要包括坡度、坡高、坡型等基本形态参数，以及土体物理力学性质；另一类为外在诱发因素，主要为降雨引发土体含水率增大，进而导致抗剪强度降低，最终促使滑坡发生的力学行为。因此，在一定的坡度、坡高、坡型情况下，当土体含水率达到什么程度时斜坡会发生黄土滑坡显得尤为重要（薛强等，2018）。

4.1 斜坡稳定性分析

当斜坡坡度、坡高、坡型不同时，其发生黄土滑坡的概率也不同。因此，将研究区范围内的斜坡单元按照不同的坡度、坡高和坡型分别进行稳定性计算。本文采用Mohr-Coulomb准则作为本构模型，选用Morgenstern- Price作为分析模型，运用GeoStudio软件中的SLOP/W模块分别进行不同含水率工况下不同斜坡类型的斜坡稳定性计算，稳定性计算所选参数如表 1。本文将城区周边斜坡单元作为一个整体进行稳定性分析，当斜坡单元内包含滑坡体时，则稳定性分析剖面切过滑坡体。

图 7为坡高50 m、坡度35°的黄土斜坡在不同坡型、不同含水率工况下的稳定性计算模型及结果，图中所列为最小稳定性系数。在清涧县城区，根据斜坡工程地质测绘，晚更新世黄土一般披覆在中更新世黄土之上，晚更新世黄土披覆厚度不等，一般5~15 m，平均约10 m。因此，在稳定性计算模型中，为简便起见，将晚更新世黄土厚度统一设定为10 m。在同一工况条件下，凸型坡最易发生滑坡，凹型坡最不易发生滑坡。其他工况条件类似。本次将最小稳定性系数>1.0的判定为稳定，最小稳定性系数≤1.0的判定为不稳定，图 8为清涧县城区斜坡单元在含水率分别为5%、10%、15%、20%、25%、30%不同工况下的稳定性计算结果。

清涧县城区及周边区域斜坡在不同含水率工况下稳定性计算结果如表 2，由表 2可知，随着土体含水率的逐渐增加，城区内稳定斜坡的面积逐渐减少，不稳定斜坡的面积逐渐增大（图 9）。当斜坡土体含水率为5%时，稳定斜坡和不稳定斜坡面积占所有斜坡面积的比例分别为98.06%和1.94%；当斜坡土体含水率增加到30%时，稳定斜坡和不稳定斜坡面积占所有斜坡面积的比例分别变为18.21%和81.79%。

同时，随着斜坡土体含水率的增加，不稳定斜坡单元内包含的已发生滑坡的数量亦逐渐增加，说明已发生滑坡的斜坡较原始的自然斜坡更容易失稳（图 10）。当斜坡土体含水率从5%增加到30%时，不稳定斜坡单元内已发生滑坡所占比例从18.84%增加到99.28%，说明有99.28%的已发生滑坡处于不稳定状态。

4.2 斜坡失稳危险性评价

降雨在没有优势入渗通道的情况下，其在黄土中的入渗深度有限（薛强等，2014）。对陕北地区黄土斜坡浅层土体含水率进行了3个水文年的连续监测（图 11），监测时间：2013年10月30日至2017年3月7日。

对监测结果进行统计分析，发现土体含水率w≤0.1的概率为0.132，土体含水率0.1＜w≤0.15的概率为0.490，土体含水率0.15＜w≤0.2的概率为0.296，土体含水率0.2＜w≤0.25的概率为0.082，土体含水率w＞0.25的概率为0（图 12）。根据黄土斜坡在不同含水率工况下的稳定性计算结果，按照土体含水率出现的概率进行斜坡失稳危险性评价，可将斜坡失稳危险性分为危险性很高、危险性高、危险性中等和危险性低4个级别（表 3）。

土体含水率w≤0.15出现的概率为0.622，即出现该含水率工况的概率很高，因此将此含水率下计算的不稳定（失稳概率很高）斜坡单元定为危险性很高区。土体含水率0.15＜w≤0.2出现的概率为0.296，即认为出现该含水率工况的概率高，因此将此含水率下计算的不稳定斜坡单元（除去危险性很高区）为危险性高区。土体含水率0.2＜w≤0.25出现的概率为0.082，即认为出现该含水率工况的概率为中等，因此将此含水率下计算的不稳定斜坡单元（除去危险性很高和危险性高区）为危险性中区。土体含水率w＞0.25出现的概率为0，因此将危险性很高、危险性高和危险性中区之外的斜坡单元均定为危险性低区（图 13）。

危险性很高区：主要分布在清涧县城清涧河一带，以及东沟东门湾一带，其他地段也有零星分布。危险性很高区包含斜坡单元共112个，面积共3.27 km²，共包含已发生滑坡点92个。

危险性高区：主要分布在清涧县城附近西砭沟、赤土沟一带，与危险性很高区交叉分布，其他地段也有零星分布。危险性高区包含斜坡单元共128个，面积共4.19 km²，共包含已发生滑坡点36个。

危险性中区：在清涧河及其一级支流沿岸斜坡地带均有分布。危险性中区包含斜坡单元共251个，面积共8.75 km²，共包含已发生滑坡点6个。

危险性低区：在清涧河及其一级支流沿岸斜坡地带均有分布。危险性低区包含斜坡单元共434个，面积共15.20 km²，共包含已发生滑坡点4个。

5 结论

（1）通过清涧县城区地质灾害调查，清涧县城区周边斜坡地带共发育滑坡138处，以浅层滑坡为主，广泛分布于黄土梁峁侧壁的沟壑地带。按照河流沟谷的发育情况和地形地貌的完整性，将清涧县城及周边区域的斜坡地带共划分为925个斜坡单元，单元面积平均为3.40×10⁴ m²。

（2）当斜坡坡度、坡高、坡型不同时，发生黄土滑坡的概率也不同。将研究区范围内的斜坡单元按照不同的坡度、坡高和坡型分别进行不同含水率工况下的斜坡稳定性计算，并结合土体含水率出现概率确定斜坡危险性。

（3）土体含水率w≤0.15出现的概率很高，在此含水率工况下计算所得的不稳定斜坡单元为危险性很高区，面积共3.27 km²，包含斜坡单元112个，包含已发生滑坡点92个；土体含水率0.15＜w≤0.2出现的概率高，在此含水率工况下计算所得的不稳定斜坡单元（除去危险性很高区）为危险性高区，面积共4.19 km²，包含斜坡单元128个，包含已发生滑坡点36个；土体含水率0.2＜w≤0.25出现的概率中，在此含水率工况下计算所得的不稳定斜坡单元（除去危险性很高和危险性高区）为危险性中区，面积共8.75 km²，包含斜坡单元251个，包含已发生滑坡点6个；土体含水率w＞0.25出现的概率低，在此含水率工况下计算所得的不稳定和稳定斜坡单元（除去危险性很高、危险性高和危险性中区）为危险性低区，面积共15.20 km²，包含斜坡单元434个，包含已发生滑坡点4个。

致谢: 衷心感谢项目组成员王佳运、贾俊、徐继维、张亮在野外调查工作中付出的辛苦劳动。