摘 要: 针对地质灾害易发性评价指标识别与构建的问题，以准确获取吉林省泥石流孕灾因子权重系数为目标，采用“3S”技术获取大量数据为基础，选取对吉林省泥石流易发性有显著影响的地貌、坡度、地势高程等 7 个孕灾因子，首次将偏好比率法与粗糙集理论组合赋权并优化，分析吉林省泥石流孕灾因子权重，并将分析结果与以往所采用的方法所取得的权重结果进行对比，结果表明: 坡度、植被覆盖率、地势高程、人口密集度因素对吉林省泥石流发育影响较显著，与以往方法相比，植被覆盖率对吉林省泥石流易发性的影响显著提高，岩组的影响显著降低，在 ARCGIS 平台下开展评价，证明该方法取得结果与实地灾情较好吻合。研究成果对地质灾害风险评价及防灾减灾具有重要参考价值。
关键词: 泥石流; 孕灾因子; 组合赋权; ARCGIS
泥石流是一种复杂的地质灾害，爆发突然，来势凶猛，危害巨大。长期以来，国内外对泥石流的研究经历了从定性到定量分析的发展过程，其中权重分析对其易发性评价结果起着重要作用，而层次分析、粗糙集理论、灰色关联度模型、神经网络、回归分析模型等在泥石流灾害权重分析中均得到了广泛的应用[1]。如刘涛等[2]，采用层次分析法对北京密云县历史活动及潜在泥石流指标量化，进而分析危险度，温智虹等[3]采用层次分析法计算长白山火山泥石流分布影响因素的权重值，李炫等[4]采用层次分析法计算西藏 30 条冰川泥石流沟影响因素权重系数，将其引入到 WKFCM 聚类算法判定泥石流危险性。层次分析法虽可使问题简化、层次化，但仅采用此客观分析法可因原始资料缺乏，选取因子较少影响评价精度。粗糙集是解决不确定和模糊问题的一种工具，从给定问题的描述集合入手，分析其内在规律，当描述集合数据较少时，其精度也随之降低。王禹萌等[5]基于粗糙集理论针对长白山火山泥石流选取 5 个因子进行重要性分析，但因样本数据相对较少，而评价过程中客观因素所占比例较大，降低了其准确性; 向灵芝等[6]利用灾害熵计算汶川震后道路泥石流各个因子的权重，分析研究区泥石流敏感性 。现今，随着 GIS 手段的不断发展，越来越多学者将 GIS 应用到泥石流灾害分析评价中，且取得了较好的效果。如苏凤环等[7]在利用层次分析法获取汶川县泥石流影响因子权重基础上，利用 GIS 分析其敏感性; 潘赟等[8]利用 GIS 对辽宁省岫岩县泥石流灾害危险性进行区划; 张具琴等[9]利用“3S”技术结合层次分析法确定汶川县泥石流各因子权重系数。通过现今对泥石流孕灾因子权重分析方法及所取得的成果可以看出，这些方法往往是基于样本数据的客观分析，而忽略了基于野外调查的主观权重判断的重要性，其判定结果与实地灾情存在一定误差。因此分析不同区域泥石流成灾特点，准确判定孕灾因子，在此基础上选取合理、适宜的权重判定方法确定权重大小是泥石流危险 性 评 价、机理分析及风险评价的重要前提。
相关期刊推荐：《中国地质灾害与防治学报》创刊于1991年。是经原国家科委正式批准，由中华人民共和国国土资源部主管，中国地质灾害研究会主办，中国地质调查局、国土资源部地质环境司、中国地质环境监测院、中国科学院力学研究所等9单位协办，反映地质灾害学科，并向国内外公开发行的学术性刊物。有投稿需求的作者，可以咨询期刊天空在线编辑。
吉林省东南部地区以山区为主，每年地质灾害发生频率较高，尤其是泥石流灾害，给当地的经济造成了较大的损失，制约了其经济发展。吉林省根据其泥石流成灾特点，开展了 1: 5 万县 ( 市) 地质灾害调查与区划工作，本文就是在此基础上，对吉林省泥石流地质灾害的地质环境进行研究，利用“3S”技术获取数据，提取泥石流影响因子，利用 TM 影像图、DEM 数字等高线图、地质图等数据将其在 ARCGIS 中矢量化，分析及筛选影响泥石流灾害的主要孕灾因子。首次提出将主观判断的偏好比率法与客观分析的粗糙集理论进行优化组合，避免了以往客观分析法中样本数据不足而影响结果精度，并应用到吉林省泥石流孕灾因子的权重分析中，将结果与其他方法所取得的权重值进行对比，在 ARCGIS 平台下归纳灾害点、植被覆盖率、岩组图层中的单元集，分析吉林省泥石流灾害点发育个数与植被覆盖率及岩组的相互关系，进而分析差异性，其研究成果对完善吉林省地区泥石流风险区划有着重要参考价值，同时对长白山火山崩滑流灾害风险评价、其它地区地质灾害孕灾因子权重获取、风险区划、防灾减灾工作有着重要意义。
1 研究区地质条件概况
吉林省位于中国东北地区中部，处于日 本、俄罗斯、韩国、朝鲜与中国东北部组成的东北亚腹心地带。东南高，西北低，地貌类型主要为东部低山丘陵、长白山中山低山等 15 个单元。其中西部主要为中高山地区，坡度可达 30° ～ 70°，而东部主要为中低山地区。气候特点是: 春季干旱、夏季炎热多雨，6 - 9 月降水量较多，对吉林省泥石流出现的汛期分析发现，25% 的泥石流均发生在暴雨当天，且为持续集中的降雨型。岩组类型较多，土体及软质岩石主要分布在平原区，极硬岩及硬质岩石主要为中厚层火山喷发而沉积形成的碎屑岩组及中厚至薄层的砂砾岩组，山区分布较多，地形坡度相对陡峭。植被覆盖率较高，林业面积居全国第 12 位，其中东部区主要为森林区域。而近代由于人类活动加强、植被破坏、风化加强，对地质环境破坏较显著，其中对于一些不合理的开挖边坡，矿山建设，废弃土体的不当堆砌等均是诱发崩塌、滑坡的主要原因，为泥石流提供丰富的物源，植被的破坏也加大了水源的汇集。
吉林省泥石流共发育有 1 407 条，如表 1 所列，主要为暴雨诱发，较多为沟谷型泥石流，规模以小型群发为主，其具有重度大、流速快、破坏力强的特点。主要分布在松花江、鸭绿江等流域。其中典型的泥石流如大砬子沟泥石流其坡度 25°，坡降 80‰，沟长约 2 000 m，流域面积 4. 5 km2 ，堆积物厚度可达 2 m 左右，由块石、碎石组成，是沟谷型泥石流; 长青村泥石流坡度在 25° ～ 50°左右，坡降大于 200‰，流域面积 1 km2 左右，相对高差 50 ～ 100 m，冲沟横断面形态呈“V”型，为沟谷型稀性泥石流，2010 年突然爆发导致 2 人死亡; 瀑布东泥石流其坡度 40°，坡降 300‰，沟长约 500 m，流域面积 0. 1 km2 ，堆积物为 2 000 m3 ，由块石、碎石组成，是坡面型泥石流，严重威胁下游聚龙泉旅游点。
2 基于组合赋权法的孕灾因子重要度分析
泥石流易发性分区评价对泥石流灾害危险性分析、风险分析及防灾减灾工作有着重要意义，其孕灾因子选取及权重分析结果直接影响着泥石流易发性评价的精度，科学、合理的分析结果可有效进行风险分区与评价，进而避免各种工程与非工程减灾措施实施过程中的盲目性，为防患于未然和减轻泥石流灾害对当地的影响提供可靠的依据[10]。
通过对吉林省泥石流点与各孕灾因子关系的统计分析如图 1 所示( 其它图略) ，选取与吉林省泥石流发育有较显著影响的地貌、坡度、地势高程、岩组、年平均降水量、人口密集度、植被覆盖率 7 个孕灾因子进行权重分析。
3 结果分析
对吉林省泥石流易发性有显著影响的孕灾因子采用组合赋权理论所取得的权重系数结果与以往客观分析方法取得的结果进行对比分析( 坡度 > 地势高程 > 人口密度 > 植被覆盖率 > 岩组 > 年平均降水量> 地貌) [15]，植被覆盖率对吉林省泥石流易发性的影响显著提高，岩组的影响显著性降低。文中在 Mpagis 中将植被覆盖率、岩组地质图纸矢量化，在图形处理模块中利用转换功能将矢量化的数据转为 shape 文件后导入到 ARCGIS 平台，并转换为1 km × 1 km 栅格数据。在 ARCGIS 平台下归纳、分析灾害点、植被覆盖率、岩组图层中的单元集，分析吉林省泥石流灾害点发育个数与植被覆盖率及岩组的相互关系( 图 2、图 3) ，进而分析其原因。
植被覆盖率较少，岩体相对裸露，有利于风化形成丰富松散的物源条件，并且有利于水源的汇集而为泥石流发生提供动力条件。吉林省植被覆盖率较高，林业面积居全国第 12 位，其中东部区主要为森林区域。而近代由于人类活动加强、植被破坏、风化加强，进而为泥石流的发育提供了物源条件。吉林省植被覆盖率在20% ～ 50%区域泥石流最发育，共 481 条。植被覆盖率大于 50% 区域泥石流发育程度降低，共 325 条，植被覆盖率小于 20% 区域因地处平原区，因此泥石流相对较少。由此可见，植被覆盖率对泥石流发育有着较为重要影响，其影响显著性大于地势高程及人口密度。
岩组是泥石流的物质来源之一，其松散碎裂及岩石本身的抗风化性能影响着物源的丰富程度[6]。吉林省岩组类型较多，通过分析 ARCGIS 所得到的吉林省泥石流点及岩组相互关系可见，土体及软质岩石主要分布在平原区，地形地貌条件不利于泥石流的发育，因此泥石流相对较少。吉林省极硬岩及硬质岩石主要为中厚层火山喷发而沉积形成的碎屑岩组及中厚至薄层的砂砾岩组，其因地形坡度相对陡峭，泥石流分布相对较多，岩组与泥石流发育点数量相关性不显著，因此岩组对吉林省泥石流的易发性影响小于年平均降水量及地貌因子。
上述结论与实地灾情更加吻合，避免了层次分析法、模糊数学法、粗糙集理论等方法在地质灾害权重分析中过多采用客观评价而忽略了主观判断重要性的弊端，结果更加科学、合理。
4 结论
文中在吉林省 1∶ 5 万地质灾害调查的基础上，根据吉林省泥石流的发育特点，选取了影响吉林省泥石流易发性的坡度、地势高程、年平均降水量等 7 个孕灾因子，采用偏好比率—粗糙集理论组合优化，分析各孕灾因子权重的大小，并将所取得的权重系数与其它方法取得的结果对比分析，得出以下结论:
( 1) 吉林省泥石流孕灾因子影响程度由大到小为坡度 > 植被覆盖率 > 地势高程 > 人口密度 > 年平均降水量 > 地貌 > 岩组。
( 2) 与其它方法取得权重系数结果主要区别为: 植被覆盖率对吉林省泥石流易发性的影响显著提高，岩组的影响显著性降低。
( 3) 在 ARCGIS 平台下分析植被覆盖率与岩组因子权重大小变化的原因后，认为该组合赋权方法在泥石流孕灾因子权重分析中所得出的结论与实地灾情更加吻合，该方法即考虑了主观因素定性分析重要性，又结合了客观分析方法的定量评价，使结果更加科学、合理。
文中采用组合赋权法得到的泥石流孕灾因子权重值将进一步完善吉林省泥石流风险区划，所采用的组合赋权法可结合承灾体暴露性、脆弱性及防灾减灾能力因子应用于 长白山地区火山崩滑流灾害的风险评价、预警中，且可拓展应用到其它地区地质灾害易发性分区评价中，对地质灾害风险分析、减灾防灾的准确评价及保证公共安全具有重要意义。
