安徽黄山市徽州区地质灾害危险性评价研究

刘 乐，杨 智，孙 健，刘 钦，彭 鹏，段俊斌

（1.安徽省地质调查院，安徽 合肥 230001；2.安徽省地质矿产勘查局332 地质队，安徽 黄山 245000）

0 引言

地质灾害的发生是内因和外因共同作用的结果。内因即为区内的地质环境，决定着地质灾害的类型、分布、规模和强度，主要包括地形地貌、岩性构造、工程地质岩组等[1]。外因是触发地质灾害的外部条件，主要包括降雨、地震、工程切坡、植被破坏、矿山开采等。

开展大比例尺高精度的地质灾害评价是当前地质灾害调查亟需的[2]。利用GIS 平台结合高分遥感影像、DEM 数据等，可以快速有效的进行大范围地质灾害危险性评价[3]。即根据研究区特点，通过对已发的地质灾害进行提取，分析其孕灾、成灾的因素，并选用恰当的数学评价模型，综合评价区内的地质灾害危险性程度，并预测区域内地质灾害易发生空间位置。进而建立针对研究区的地质灾害危险性评价模型。

1 研究区概况

研究区位于皖南山区，行政区划隶属黄山市徽州区，面积约398.49 km2（图1）。以中低山地貌为主，总体呈北高南低，海拔标高一般在200～1 400 m。地层岩性以浅变质砂岩、粉砂岩、凝灰质粉砂岩为主。区域变质、变形为板岩、千枚状板岩。岩石易风化，残坡积层较发育。区域构造上，属下扬子陆块、江南造山带、江南古陆隆起带之历口构造带，构造线走向主要为北东向，其次为北西向。区内水系属钱塘江流域、新安江水系，最低侵蚀基准面标高为196 m 左右。区内多年平均降雨量1 708.1 mm，春、夏汛期为地质灾害高发期[4]。

图1 研究区交通位置图Fig.1 Traffic map of the research area

2 信息量模型原理

信息量模型是一种定量分析方法，其物理意义明确，广泛应用于区域地质灾害危险性评价[5-6]，评价过程中可较好的反映致灾因子和地质灾害的关联性[7]。通过模型评价，能够直观的反映出研究区内各致灾因子对于形成地质灾害的敏感度和贡献率。

式中：x i——成灾因素x中的第i区间；

Y——成灾因素x中第i区间地质灾害发生的信息量值；

B——地质灾害事件；

Ni——研究区内包含评价因素xi的单元数且存在地质灾害的单元数；

Si——研究区内包含评价因素xi的单元数；

N——研究区存在地质灾害单元总数；

S——研究区中评价单元的总数。

3 研究区地质灾害危险性评价

3.1 评价因子的选取

研究区所处的皖南山区，地质构造复杂，山体风化强度较高，残坡积发育，丰乐河及其支流广布于研究区内。区内道路、房屋的建设多采用切坡施工。另外，该地区广泛种植茶叶，茶园多开垦于陡峻的山坡坡面，人类活动强度大（图2）。结合研究区实际情况和前人工作经验，本次地质灾害危险性评价选取了高程、坡度、坡向、断裂构造、土地覆盖类型、水系、工程地质岩组、人类活动强度，共8 项主要因素作为评价指标。

通过收集以往地质灾害调查成果及野外查证，研究区内已发生地质灾害点215 处。根据已知灾害点，分别结合选取的8 项主要因素进行分析，求得评价因素中各因子对地质灾害的“贡献度”，将8 项评价因素叠加，进而得到研究区地质灾害危险性评价结果。

3.1.1 高程

研究区地形地貌跨度较大，总体地势西北高南东低，北部为中低山区，山体走向以北西、北东为主，海拔500～1 400 m，相对高差600～800 m，局部地区800～1 000 m；南部为高丘地貌，广泛分布在中低山侧，高程200～500 m，相对高差50～250 m。

根据研究区地形特点，将高程划分为0～200 m、200～300 m、300～400 m、400～500 m、500～700 m、700～1 300 m 共6 个级别，见图2（b）。

通过高程因素信息量分析，300～700 m 高程区间地质灾害发生的可能性较大，而高程低于200 m 或高于700 m，地质灾害发生概率较低（表1）。

表1 高程分级及信息量统计表Table 1 Satistics of elevation classification and information

3.1.2 坡度

坡度直接影响坡面上坡积物的厚度、物质的稳定性和水动力条件，从而影响地质灾害发生的强度和规模[8]。

本次将坡度划分为0°～10°、10°～20°、20°～30°、30°～40°、40°～50°、50°～90°共6 个级别，见图2（c）。

图2 地质灾害点密度与致灾因子分级图Fig.2 Density maps of geological disaster points andclassification diagram of the impact factors of geological disasters

通过坡度对地质灾害提供的信息量图，可见坡度在20°～40°区间，地质灾害发生可能性较大，而坡度小于20°时，灾害发生概率较小，对于坡度大于50°，由于研究区内此坡度区域极少，故地灾信息量较低（表2）。

表2 坡度分级及信息量统计表Table 2 Satistics of slope classification and information

3.1.3 坡向

坡向即坡面的朝向。不同朝向的坡面，坡体受太阳辐射强度，各种物理化学的风化作用程度不同。

将坡向每隔90°划分为四个方位区间，将分别为东坡45°～135°、南坡135°～225°、西坡225°～315°、北坡315°～45°，见图2（d）。

通过坡向对地质灾害提供的信息量图，可见北坡地质灾害发生可能性较小，其他方向地质发生概率较大，特别是南坡（表3）。

表3 坡向分级及信息量统计表Table 3 Satistics of slope classification and information

3.1.4 断裂

研究区内断裂构造发育，且以北东及北西向断裂为主北东东向最为发育（规模大、密集），北西向次之。区内大型断裂切割本区古老地层，同时伴生次级断裂和构造裂隙，使岩体结构松散，岩石破碎及风化，是引发地质灾害的主要因素[9]。

根据遥感解译断裂构造的位置，以与断裂构造的距离为依据划分缓冲区，共划分5 个级别，分别为0～300 m、300～600 m、600～1 000 m、1 000～1 500 m、＞1 500 m，见图2（e）。

通过距离断裂级别对地质灾害提供的信息量图，可见在距离断裂带300 m 内发生地质灾害的可能性最大，1 000 m 以上则地质灾害的发生概率较低（表4）。

表4 断裂分级及信息量统计表Table 4 Satistics of fracture classification and information

3.1.5 水系

河流的侵蚀是导致地质灾害发生的重要因素，主要表现为侵蚀作用对斜坡前缘抗力的削弱和临空面的增加造成斜坡失稳[10]。

研究区内主要河流为丰乐河及其支流，对河流根据其线密度进行统计，密度越大，说明河流沟谷越多，地面越破碎，地面物质稳定性越低，造成的灾害越多。

在ArcGIS 中，将河流进行遥感解译，并形成线矢量文件，将线文件进行密度分析，计算公式为：

式中：Ds——水系密度；

∑L——研究区水系总长度/km；

A——研究区面积/km2。

通过密度分析，研究区水系密度值范围为0～6.264 4，通过与遥感影像图进行比对，将密度范围划分为3 个区间，分别为0～1.405 2、1.405 2～2.869 5、2.869 5～6.264 4，见图2（f）。通过水系密度对地质灾害提供的信息量可见，在水系密度高的区域，地质灾害相对概率也较高（表5）。

表5 水系密度分级及信息量统计表Table 5 Satistics of the river system density classification and information

3.1.6 土地利用类型

研究区为黄山毛峰的主产区，陡坡茶园在区内分布广泛，因此在土地类型的划分上，将茶园进行单独分类，共划分为道路、耕地、城镇用地、水域、林地、茶园6 种类型，分析不同用地类型对引起地质灾害所提供的信息量，见图2（g）。

依据各土地类型对地质灾害提供的信息量，可见道路、茶园、城镇用地3 种土地类型发生地质灾害的可能性较大（表6）。

表6 各土地利用类型信息量统计表Table 6 Statistical table of information quantity of land use types

3.1.7 工程地质岩组

岩土体作为斜坡的基本组成，其控制着地质灾害的形成、分布和规模[11]。研究区按岩石强度划分为①坚硬中厚层砂岩岩组（NH1x）、②较坚硬层状、板状砂岩夹板岩岩组（PT2n）、③坚硬中厚层状变质砂岩岩组（PT2d）、④坚硬块状花岗闪长岩岩组（γδ）、⑤第四系松散岩组（Q4）、⑥残坡积层碎石土（Qdl+el）、⑦强风化层、⑧冲洪积层（Qapl），见图2（h）。

通过分析，8 类工程地质岩组中，残坡积层碎石土（Qdl+el）和强风化层两类对地质灾害贡献的信息量最大（表7）。

表7 各工程地质岩组信息量统计表Table 7 Statistical table of information of each engineering geological rock formation

3.1.8 人类活动强度

在ArcGIS 中，将道路路网和房屋范围进行遥感解译，并形成线矢量文件，将线文件进行密度分析，线网密度越大，说明人类活动越强烈，引发的地质灾害越多。

通过线密度分析，徽州区人类活动密度值范围为0～12.295 5，通过与遥感影像图进行比对，将密度范围划分为4 个区间，分别为0～1.418 5（人类活动强度极低）、1.418 5～4.530 9（人类活动强度一般）、4.530 9～5.882 1（人类活动强度较高）、5.882 1～12.295 5（人类活动强度极高），见图2（i）。

通过人类活动强度对地质灾害提供的信息量图，可见在人类活动强度较高、极高区间，地质灾害发生较频繁，在人类活动强度极低区间地质灾害发生可能性较低（表8）。

表8 人类活动强度分级及信息量统计表Table 8 Satistics of human activity intensity classification and information

3.2 地质灾害危险性评价

3.2.1 评价网格划分

对于基于GIS 栅格运算的地质灾害易发性区段评价中，研究区中各致灾因子图层的评价单元，选用正方形标准栅格作为评价单元，确定栅格单元大小的经验公式为：

式中：Gs——适宜栅格大小；

S——地质灾害评价比例尺的倒数。

评价单元的划分会直接影响评价结果的合理性[12]，利用公式（3）并结合本次地质灾害评估的精度要求，使用13.5 m×13.5 m 的栅格大小作为评价单元，研究区共划分约2 186 478 个评价单元。

3.2.2 地质灾害危险性评价结果

将所有致灾因子信息量求和，得到2 186 478 个均一条件单元的总信息量值，其范围为-4.732 721～3.449 854，数值越大，对地质灾害发生的“贡献率”越大，地质灾害越容易发生（图2）[13-16]。

根据各评价单元的信息量值，采用自然间断点法，取-1.18、-0.03、1.26 为分界点，将研究区按危险度划分为地质灾害不易发区、低易发区、中易发区和高易发区（表9、图3）[17-19]。

表9 地质灾害危险度分区面积统计表Table 9 Statistical table of geological hazard area

图3 研究区地质灾害危险性评价图Fig.3 Geological hazard risk assessment map of the research area

4 成果验证

通过本次地质灾害危险性评价结果与野外调查成果的215 处地质灾害点进行叠合分析，灾点分布情况为：高易发区132 处、中易发区80 处、低易发区3 处、不易发区0 处（表10）。地质灾害点的分布结果与地质灾害危险性评价区划相符，评价结果较合理。

表10 不同危险度分区内地质灾害点数量统计表Table 10 Statistical table of the number of geological disaster points in different risk zones

5 结束语

通过研究区数据分析，得出以下结论：

（1）结合野外实际调查的地质灾害点成果，选取高程、坡度、坡向、断裂、水系、土地利用类型、工程地质岩组、人类活动强度等8 个致灾因子，使用信息量模型法对研究区进行地质灾害危险性评价，其中中易发区和高易发区面积分别为116.35，19.86 km2，占研究区总面积的29.2%和4.98%。

（2）根据评价结果，研究区内地质灾害的分布主要受到断裂带的控制，同时在河流、道路、茶园附近，坡度20°～40°的坡面，第四系松散层和强风化层覆盖区域均为地质灾害高发区，需要在重点区域加强地质灾害的防治工作。