Processing math: 100%
  • 中文核心期刊
  • 中国科技核心期刊
  • 中国科学引文数据库核心期刊

区域地质灾害易发性评价的证据权法原理与实践

冯振, 陈亮, 王立朝, 侯圣山, 田怡帆, 刘明学

冯振, 陈亮, 王立朝, 侯圣山, 田怡帆, 刘明学. 2024: 区域地质灾害易发性评价的证据权法原理与实践. 地质通报, 43(7): 1255-1265. DOI: 10.12097/gbc.2023.02.034
引用本文: 冯振, 陈亮, 王立朝, 侯圣山, 田怡帆, 刘明学. 2024: 区域地质灾害易发性评价的证据权法原理与实践. 地质通报, 43(7): 1255-1265. DOI: 10.12097/gbc.2023.02.034
Feng Z, Chen L, Wang L C, Hou S S, Tian Y F, Liu M X. Principle and application of the weight-of-evidence method in regional landslide susceptibility assessment. Geological Bulletin of China, 2024, 43(7): 1255−1265. DOI: 10.12097/gbc.2023.02.034
Citation: Feng Z, Chen L, Wang L C, Hou S S, Tian Y F, Liu M X. Principle and application of the weight-of-evidence method in regional landslide susceptibility assessment. Geological Bulletin of China, 2024, 43(7): 1255−1265. DOI: 10.12097/gbc.2023.02.034

区域地质灾害易发性评价的证据权法原理与实践

基金项目: 中国地质调查局项目《重点地区地质灾害风险调查与评价》(编号:DD20221748)、《洮河流域中游灾害地质调查》(编号:DD20190645)和贵州电科院科技项目《山区电网地质灾害危险性快速与动态评价技术研究》(编号:0002200000089405)
详细信息
    作者简介:

    冯振(1985− ),男,博士,正高级工程师,从事工程地质、地质灾害等研究。E−mail: fengzhencgs@126.com

    通讯作者:

    陈亮(1983− ),男,硕士,正高级工程师,从事地质灾害调查与风险评价研究。E−mail: 109209183@qq.com

  • 中图分类号: P628;P642.2

Principle and application of the weight-of-evidence method in regional landslide susceptibility assessment

  • 摘要:

    在阐述证据权法概念与原理的基础上,以西北黄土高原滑坡、崩塌地质灾害高易发的兰州市区为评价区,详解基于证据权法的区域地质灾害易发性评价的技术流程,利用GIS工具,开展评价指标选取与相关性分析、证据因子权重与后验概率分析、评价模型检验、易发程度评价。分析结果显示,坡度、坡向、地层岩性、断层、绿化灌溉、土地利用类型等地质灾害影响因素之间呈无相关或弱相关性,符合相互独立的要求,组成的模型ROC曲线(接受者操作特征曲线)AUC值(曲线下面积)为0.85,具有较高的预测精度。理论与实践表明,与其他统计学方法相比,证据权法同时考虑了评价指标中存在地质灾害与不存在地质灾害2种情况,避免了当地质灾害不存在时权重偏大导致结果失真的问题。证据权法结果为地质灾害发生的绝对概率,可以作为不同区域、地区的地质灾害易发程度对比的依据。

    Abstract:

    Based on elaboration of the concept and principle of WoE, Lanzhou urban area, which is highly prone to slides and falls in the Loess Plateau of Northwest China, is taken as example to detail technical process of regional geological hazards susceptibility assessment based on the WoE. GIS tools are used to carry out evaluation index selection and correlation analysis, weight and posterior probability analysis of evidence factors, assessment model verification, and susceptibility degree classification and assessment. The analysis results show that there is no correlation or weak correlation between the factors affecting geological hazards such as slope, slope direction, stratum lithology, faults, waterfalls, green irrigation and land use type, which meets the requirements of mutual independence. AUC (Area Under Curve) of the ROC curve (Receiver Operating Characteristic Curve) of the model is 0.85, which means high prediction accuracy. Theoretical analysis and empirical application have demonstrated that, in comparison with other statistical methods, the WoE takes into account both probability of the presence and absence of geological hazard within the evaluation factors. Thereby, it prevents distortion caused by overestimation of weights when geological hazards are absence. The result of the WoE is the absolute probability of the occurrence of geological hazards, which can serve as a basis for comparing the susceptibility of geological hazards from different areas.

  • 吉塘变质杂岩系分布于澜沧江西南,沿他念他翁山链呈北西—南东向长条状展布,西起巴青县老巴青乡,经类乌齐、吉塘,东至左贡县(图 1),构成北澜沧江北侧他念他翁山脉主体,东部以前称为北澜沧江变质岩系,交通要道位于吉塘镇,也称为吉塘变质岩系、吉塘群、吉塘变质核杂岩,在类乌齐一带又称类乌齐群[1]。在时代归属上有多种划分方案,20世纪50年代初李璞等在三江地区开展地质调查,将山露于吉塘—类乌齐地区的变质岩命名为“吉塘变质岩、“泄巴变质岩和“类乌齐变质岩”,认为其时代为前寒武纪[2-3]。四川地质局第三区域地质测量大队将其划分为古生界变质岩系,西藏地勘局第一地质大队将其划分为变质岩系酉西杂岩群,形成时代为古生代。周云生等[4]将其称为北澜沧江变质岩带,划分为下部酉西杂岩群,上部卡贡群,时代为晚古生代。《“三江”区域地质图》称为吉塘群,其时代划归前石炭纪[5]。艾长兴等[6]称其为吉塘群,时代为泥盆系;雍永源等[7]称其为北澜沧江变质岩带,划分为上部酉西群,下部吉塘群,酉西群为早古生代,吉塘群归属前寒武纪。王根厚[1]称为吉塘变质杂岩系,划分为酉西岩组和吉塘岩组构造片岩,时代置于前石炭纪。

    图  1  研究区大地构造位置图(a)和藏东吉塘变质杂岩空间分布图(b)
    Figure  1.  Tectonic location map of study area(a)and spatial distribution map of Jitang metamorphic complex, eastern Tibet(b)

    本文沿用雍永源等[7]的北澜沧江变质岩带,划分为上部酉西群、下部吉塘群。西藏自治区地热地质大队对吉塘群中具片麻状、条纹(条痕)状、眼球状、混合岩化现象的岩石进行锆石U-Pb定年研究,表明其形成于印支期,原岩为(碎裂岩化、糜棱岩化)中—酸性花岗岩,与吉塘复式花岗岩形成时代一致,将其归于吉塘复式花岗岩。本文所指吉塘群岩性组合为片麻岩、斜长角闪岩、片岩等。

    近年随着锆石U-Pb同位素定年技术广泛应用,众多研究者从吉塘群片麻岩中不断得出众多年龄。李才等[8]在类乌齐地区吉塘群花岗片麻岩得到锆石U-Pb年龄为254~227 Ma。沙绍礼等[9]利用锆石U-Pb和白云母Ar-Ar法分别得到240±12 Ma和226±2 Ma的年龄,认为吉塘变质杂岩体形成于中生代,很可能不存在变质基底。陶琰等[10]对吉塘镇吉塘群出露的花岗岩采用U-Pb定年,获得其年龄为249~220 Ma,花岗岩岩浆源岩为具有变质杂砂岩成分的古元古代吉塘群片麻岩,由此认为吉塘变质杂岩体中生代的锆石年龄很可能代表古特提斯洋俯冲岩浆事件的年龄,并不代表变质基底的形成时代。

    目前关于吉塘群的形成时代认识不统一,报道的形成时代从元古宙—中生代,跨度很大。本文以察雅县吉塘镇多穷沟出露的吉塘群为研究对象,对黑云二长片麻岩进行LA-ICP-MS锆石年龄研究,同时与类乌齐地区吉塘群对比,在此基础上限定吉塘群形成时代的下限,讨论成岩后期变质事件,研究结果有助于解析藏东吉塘群变质结晶基底的形成时代及后期的变质演化过程。

    研究区地处羌北-昌都地块、双湖-查吾拉-澜沧江结合带及羌塘-左贡地块交汇部位,大地构造位置独特(图 1-a)。地层分区为羌塘-三江地层大区之北澜沧江构造地层分区(图 2-a), 构造运动十分发育,吉塘镇吉塘群位于桑多-吉塘逆冲推覆断裂之间,地层呈北西—南东向展布,与区域构造迹线一致(图 2-b)。吉塘群主要出露岩性为灰色、浅灰色、灰白色黑云二长片麻岩片麻岩、黑云斜长片麻岩、斜长角闪岩、二云片岩,以残留体或捕虏体形式产于印支期吉塘复式花岗岩中,受印支期岩浆改造程度较高。

    图  2  藏东吉塘地区地层分区简图(a)和吉塘镇吉塘岩群地质图(b)
    T3—上三叠统沉积岩地层;C1kg—下石炭统卡贡组;Pt3Jbmg—吉塘岩群黑云二长片麻岩;Pt3Jabl—吉塘岩群斜长角闪岩残留体;Pt3Jplg—吉塘岩群黑云斜长片麻岩残留体;Pt3Jbmg+mis—吉塘岩群黑云二长片麻岩、二云片岩残留体;cs—酉西岩群变质砾岩;PzYqs+mis—酉西岩群石英片岩、二云片岩;sls—酉西岩群钠长片岩;hos—酉西岩群角闪片岩;ggT3—上三叠统片麻状花岗岩;ηγT3—晚三叠世二长花岗岩;γσοT3—晚三叠世英云闪长岩;ξγT3—晚三叠世正长花岗岩;γσT3—晚三叠世花岗闪长岩;γσπT3—晚三叠世花岗闪长斑岩;1—逆断层;2—性质不明断层;3—侵入接触界线;4—混合岩化作用;5—样品采集位置及编号
    Figure  2.  Simplified map showing stratigraphic division in Jitang area(a)and simplified geological map of Jitang Group in Jitang town(b), eastern Tibet

    锆石分选、微量元素分析及U-Pb同位素定年均在武汉上谱分析科技有限责任公司完成。原岩样品经过粉碎、淘洗后去除轻矿物部分,重砂部分经电磁选后得到含有少量杂质的锆石样品,最后在双目镜下挑选出锆石晶体。选择晶形较好、无裂隙的锆石颗粒粘贴在环氧树脂表面,制成锆石样品靶,打磨样品靶,使锆石的中心部位暴露出来进行抛光。对锆石进行反射光、透射光显微照相和阴极发光(CL)图像分析,根据反射光、透射光及锆石CL图像,选择代表性的锆石颗粒和区域进行U-Pb测年。

    锆石微量元素及U-Pb同位素定年利用LA-ICP-MS分析完成,激光剥蚀系统为GeolasPro,ICP-MS为Agilent 7700e,激光剥蚀过程中采用氦气作载气、氩气为补偿气以调节灵敏度,分析精度优于2%,在等离子体中心气流中加入少量氮气,以提高仪器灵敏度、降低检出限和改善分析精密度[11]。另外,激光剥蚀系统配置了信号平滑装置,即使激光脉冲频率低达1 Hz,采用该装置后能获得光滑分析信号[12];时间分辨分析数据包括20~30 s空白信号和50 s样品信号;对分析数据离线处理采用ICP-MS DataCal[13-14]软件完成;详细的仪器操作条件和数据处理方法见Liu等[14-15]

    锆石微量元素含量利用USGS参考玻璃Nist610作为外标,Si为内标的方法进行定量计算;U-Pb同位素定年中采用标准锆石91500为外标进行同位素分馏校正,每分析5个样品点分析测试1次91500;对于与分析时间有关的U-Th-Pb同位素比值漂移,利用91500的变化采用线性内插的方式进行校正,标准锆石91500的U-Th-Pb同位素比值推荐值据Wiedenbeck等[16]。样品锆石的U-Pb年龄谐和图绘制和年龄加权平均值计算均采用Isoplot/Ex_ver3[17]完成。

    从黑云二长片麻岩样品D7769中分选出较多的锆石颗粒,多呈圆状、次圆状、碎片状、短柱状等,外部形态显示磨蚀程度不等。大多数锆石为40~180 μm,长宽比介于1.2~3.3之间。阴极发光图像按锆石形态可分为2类,一类锆石具有明显的核-边结构,核部较宽,有明显的岩浆结晶环带,是残留核(图 3,点5、54、11等),边部较窄,隐约可见环带结构,为变质边(图 3,点3、4、7等);另一类为无核-边结构,隐约可见岩浆环带结构的板状锆石,未受到变质作用影响(图 3,点39、50)。本次研究从中随机挑选105颗锆石进行U-Pb年龄测定,测试时避开了裂隙和包体。

    图  3  吉塘群黑云二长片麻岩代表性锆石阴极发光(CL)图像
    Figure  3.  Representative zircon CL images of the biotite monzogneiss from the Jitang Group

    样品LA-ICP-MS原位U-Pb年龄测试结果见表 1,各测点Th和U含量及Th/U值差异较大,U含量为10×10-6~1752×10-6,Th含量为4×10-6~2135×10-6,Th/U值为0.02~2.87。锆石原位U-Pb测年结果显示,年龄跨度较大,排除谐和度较差及受热液混杂影响的18、20、56、74、77、85、86、102、103测点后,获得96个有效分析点,谐和度较高(图 4-a),锆石谐和年龄变化范围为2818~203 Ma,年龄结果大于1000 Ma时采用207Pb/206Pb年龄,其余采用206Pb/238U年龄。

    表  1  吉塘岩群黑云二长片麻岩LA-ICP-MS锆石U-Th-Pb分析结果
    Table  1.  LA-ICP-MS zircon U-Th-Pb dating results of the biotite monzogneiss in the Jitang Group
    测点号 元素含量10-6 Th/U 同位素比值 年龄/Ma 谐和
    度/%
    采用年
    龄/Ma
    Th U 207Pb/206Pb 207Pb/235U 206Pb/238U 207Pb/206Pb 207Pb/235U 206Pb/238U
    D7769-01 237 429 0.55 0.0553 0.0017 0.5254 0.0160 0.0683 0.0006 433 70.4 429 10.7 426 3.9 99 426
    D7769-02 78.5 562 0.14 0.0541 0.0011 0.3668 0.0126 0.0489 0.0013 372 48.1 317 9.3 308 7.7 96 308
    D7769-03 18.1 265 0.07 0.0509 0.0024 0.2264 0.0109 0.0320 0.0004 239 109.2 207 9.0 203 2.5 98 203
    D7769-04 10.8 489 0.02 0.0526 0.0015 0.2918 0.0111 0.0399 0.0008 309 58.3 260 8.8 252 5.1 96 252
    D7769-05 72.3 143 0.51 0.0720 0.0017 1.3352 0.0325 0.1343 0.0011 985 17.6 861 14.1 812 6.5 94 812
    D7769-06 69.3 521 0.13 0.0728 0.0013 1.5615 0.0353 0.1554 0.0023 1007 37.8 955 14.0 931 13.0 97 931
    D7769-07 163 968 0.17 0.0498 0.0014 0.2473 0.0073 0.0357 0.0004 187 66.7 224 5.9 226 2.4 99 226
    D7769-08 106 220 0.48 0.0975 0.0016 3.6266 0.0632 0.2692 0.0022 1577 31.5 1555 13.9 1537 11.4 98 1577
    D7769-09 262 1388 0.19 0.0495 0.0013 0.2270 0.0060 0.0330 0.0003 172 63.0 208 4.9 210 1.6 99 210
    D7769-10 7.85 426 0.02 0.0554 0.0014 0.3327 0.0103 0.0432 0.0007 428 53.7 292 7.9 273 4.3 93 273
    D7769-11 156 388 0.40 0.0742 0.0013 1.7650 0.0356 0.1720 0.0021 1048 35.2 1033 13.1 1023 11.8 99 1048
    D7769-12 69.6 628 0.11 0.0573 0.0017 0.4574 0.0171 0.0573 0.0012 502 63.0 382 11.9 359 7.2 93 359
    D7769-13 325 873 0.37 0.0547 0.0015 0.4543 0.0144 0.0602 0.0011 398 58.3 380 10.1 377 6.9 99 377
    D7769-14 449 673 0.67 0.0601 0.0011 0.7638 0.0192 0.0917 0.0016 609 43.5 576 11.1 566 9.5 98 566
    D7769-15 211 685 0.31 0.0511 0.0016 0.2506 0.0080 0.0355 0.0004 256 74.1 227 6.5 225 2.5 99 225
    D7769-16 39.3 209 0.19 0.0638 0.0015 0.9704 0.0269 0.1100 0.0017 744 50.0 689 13.9 673 9.9 97 673
    D7769-17 286 486 0.59 0.0711 0.0011 1.5282 0.0252 0.1554 0.0012 961 33.3 942 10.1 931 6.5 98 931
    D7769-18 69.4 494 0.14 0.0516 0.0021 0.2854 0.0169 0.0393 0.0014 265 99.1 255 13.4 248 8.7 97 248
    D7769-19 17.2 48.4 0.35 0.0711 0.0024 1.3701 0.0450 0.1399 0.0014 961 63.9 876 19.3 844 8.2 96 844
    D7769-20 287 510 0.56 0.0702 0.0012 0.9926 0.0187 0.1022 0.0008 1000 36.0 700 9.6 627 4.6 89 627
    D7769-21 320 1752 0.18 0.0783 0.0010 1.7517 0.0243 0.1616 0.0012 1155 26.2 1028 9.0 966 6.5 93 966
    D7769-22 156 391 0.40 0.0729 0.0011 1.5791 0.0254 0.1565 0.0011 1011 26.9 962 10.0 937 6.3 97 937
    D7769-23 113 412 0.27 0.0693 0.0012 1.1428 0.0403 0.1182 0.0035 909 37.0 774 19.1 720 20.4 92 720
    D7769-24 166 240 0.69 0.0726 0.0014 1.7374 0.0341 0.1727 0.0013 1011 40.7 1022 12.7 1027 7.2 99 1011
    D7769-25 69.4 167 0.42 0.0701 0.0017 1.4609 0.0497 0.1487 0.0033 931 54.6 914 20.5 894 18.7 97 894
    D7769-26 209 220 0.95 0.0709 0.0014 1.5404 0.0303 0.1571 0.0011 954 36.1 947 12.1 941 5.9 99 941
    D7769-27 27.8 35.6 0.78 0.0814 0.0029 1.8069 0.0631 0.1610 0.0016 1231 68.5 1048 22.8 962 9.1 91 962
    D7769-28 17.0 306 0.06 0.0549 0.0016 0.3046 0.0089 0.0403 0.0005 409 64.8 270 6.9 255 3.2 94 255
    D7769-29 43.6 992 0.04 0.0528 0.0011 0.2955 0.0079 0.0404 0.0006 320 48.1 263 6.2 255 4.0 96 255
    D7769-30 101 217 0.47 0.0811 0.0016 1.7554 0.0382 0.1567 0.0020 1233 37.5 1029 14.1 938 11.3 90 938
    D7769-31 146 323 0.45 0.0716 0.0014 1.3836 0.0273 0.1400 0.0011 974 6.5 882 11.7 844 6.1 95 844
    D7769-32 123 295 0.42 0.1990 0.0028 12.4358 0.2954 0.4512 0.0084 2818 24.2 2638 22.4 2401 37.5 90 2818
    D7769-33 329 986 0.33 0.0604 0.0009 0.7903 0.0125 0.0948 0.0006 617 33.3 591 7.1 584 3.5 98 584
    D7769-34 302 638 0.47 0.0694 0.0011 1.2213 0.0236 0.1275 0.0016 911 33.0 810 10.8 773 8.9 95 773
    D7769-35 3.50 10.4 0.34 0.0708 0.0056 1.3547 0.0974 0.1418 0.0027 950 157.3 870 42.0 855 15.2 98 855
    D7769-36 69.1 485 0.14 0.0714 0.0013 1.5198 0.0337 0.1539 0.0016 969 38.9 938 13.6 923 9.0 98 923
    D7769-37 245 469 0.52 0.0679 0.0014 1.0308 0.0499 0.1071 0.0041 866 43.7 719 24.9 656 23.8 90 656
    D7769-38 902 881 1.02 0.0685 0.0011 1.1177 0.0189 0.1182 0.0009 883 33.3 762 9.1 720 4.9 94 720
    D7769-39 551 400 1.38 0.1323 0.0019 5.9745 0.1035 0.3263 0.0035 2129 23.9 1972 15.1 1821 17.1 92 2129
    D7769-40 110 637 0.17 0.0770 0.0011 1.8225 0.0285 0.1713 0.0014 1120 25.0 1054 10.3 1019 7.8 96 1120
    D7769-41 108 133 0.82 0.1169 0.0020 5.3808 0.0927 0.3329 0.0026 1910 30.4 1882 14.8 1852 12.7 98 1910
    D7769-42 201 497 0.40 0.0765 0.0013 1.8877 0.0336 0.1782 0.0014 1109 34.4 1077 11.8 1057 7.9 98 1109
    D7769-43 188 322 0.58 0.0687 0.0013 1.4397 0.0321 0.1512 0.0020 900 40.7 906 13.4 908 11.3 99 908
    D7769-44 73.2 442 0.17 0.0709 0.0012 1.4831 0.0241 0.1512 0.0009 954 33.3 924 9.9 907 5.1 98 907
    D7769-45 99.0 172 0.58 0.0759 0.0013 1.8348 0.0337 0.1748 0.0015 1100 34.1 1058 12.1 1039 8.2 98 1100
    D7769-46 235 285 0.82 0.0701 0.0012 1.4392 0.0262 0.1485 0.0013 931 37.0 905 10.9 893 7.3 98 893
    D7769-47 349 498 0.70 0.0764 0.0011 1.8723 0.0289 0.1771 0.0013 1106 34.3 1071 10.2 1051 7.2 98 1106
    D7769-48 388 545 0.71 0.0716 0.0012 1.4402 0.0254 0.1454 0.0012 974 35.2 906 10.6 875 7.0 96 875
    D7769-49 68.2 203 0.34 0.0770 0.0015 1.9403 0.0392 0.1820 0.0015 1122 38.9 1095 13.6 1078 8.3 98 1122
    D7769-50 173 300 0.58 0.0713 0.0013 1.2264 0.0242 0.1244 0.0012 965 41.7 813 11.0 756 6.7 92 756
    D7769-51 38.0 33.7 1.13 0.0605 0.0032 0.8973 0.0443 0.1097 0.0015 620 114.8 650 23.7 671 9.0 96 671
    D7769-52 169 521 0.32 0.0712 0.0012 1.3493 0.0238 0.1372 0.0012 963 35.2 867 10.3 829 6.8 95 829
    D7769-53 525 813 0.65 0.0615 0.0011 0.8507 0.0154 0.1001 0.0007 657 38.9 625 8.4 615 3.9 98 615
    D7769-54 144 138 1.05 0.0627 0.0018 0.9182 0.0262 0.1063 0.0010 698 65.7 661 13.9 651 6.0 98 651
    D7769-55 77.4 179 0.43 0.0679 0.0018 1.2718 0.0353 0.1361 0.0019 866 55.6 833 15.8 823 10.5 98 823
    D7769-56 51.7 456 0.11 0.0658 0.0021 0.6411 0.0388 0.0668 0.0028 800 66.7 503 24.0 417 16.8 81 417
    D7769-57 378 761 0.50 0.0740 0.0016 1.8320 0.0421 0.1789 0.0014 1043 44.1 1057 15.1 1061 7.9 99 1043
    D7769-58 117 235 0.50 0.0683 0.0019 1.3782 0.0376 0.1460 0.0009 880 89.8 880 16.0 878 5.3 99 878
    D7769-59 745 872 0.86 0.0679 0.0019 1.0327 0.0313 0.1096 0.0013 878 57.4 720 15.7 670 7.4 92 670
    D7769-60 28.9 135 0.21 0.0642 0.0025 1.0702 0.0401 0.1213 0.0011 750 83.3 739 19.7 738 6.1 99 738
    D7769-61 231 598 0.39 0.0712 0.0012 1.5191 0.0257 0.1542 0.0011 965 35.2 938 10.4 924 6.2 98 924
    D7769-62 92.4 719 0.13 0.0689 0.0011 0.9865 0.0188 0.1034 0.0012 898 33.3 697 9.6 634 7.1 90 634
    D7769-63 55.2 561 0.10 0.0570 0.0013 0.3676 0.0087 0.0467 0.0005 500 48.1 318 6.5 294 3.2 92 294
    D7769-64 117 179 0.65 0.0736 0.0015 1.7893 0.0377 0.1757 0.0018 1031 40.7 1042 13.7 1044 10.1 99 1031
    D7769-65 96.7 1179 0.08 0.0641 0.0009 0.8916 0.0145 0.1003 0.0009 746 29.6 647 7.8 616 5.5 95 616
    D7769-66 261 285 0.92 0.1034 0.0016 4.4268 0.0750 0.3090 0.0029 1687 28.2 1717 14.1 1736 14.4 98 1687
    D7769-67 166 358 0.46 0.0661 0.0012 1.3818 0.0260 0.1508 0.0016 809 41.7 881 11.1 906 8.7 97 906
    D7769-68 103 212 0.49 0.0746 0.0015 1.3879 0.0336 0.1336 0.0016 1057 40.7 884 14.3 808 8.9 91 808
    D7769-69 101 140 0.72 0.0714 0.0015 1.4178 0.0333 0.1436 0.0020 969 44.4 896 14.0 865 11.5 96 865
    D7769-70 85.4 207 0.41 0.1693 0.0022 11.2806 0.1743 0.4808 0.0048 2551 21.9 2547 14.5 2531 21.1 99 2551
    D7769-71 30.9 48.5 0.64 0.0688 0.0023 1.1984 0.0380 0.1272 0.0018 894 63.9 800 17.6 772 10.0 96 772
    D7769-72 164 330 0.50 0.0656 0.0013 1.1804 0.0256 0.1298 0.0014 794 40.7 792 11.9 787 7.9 99 787
    D7769-73 137 224 0.61 0.0723 0.0014 1.4871 0.0333 0.1484 0.0019 995 40.0 925 13.6 892 10.8 96 892
    D7769-74 172 679 0.25 0.2058 0.0028 11.5290 0.2069 0.4040 0.0051 2873 21.9 2567 16.9 2188 23.6 84 2873
    D7769-75 251 1393 0.18 0.0720 0.0009 1.3489 0.0244 0.1352 0.0018 987 25.9 867 10.6 818 10.5 94 818
    D7769-76 156 194 0.80 0.0774 0.0014 2.0020 0.0384 0.1872 0.0020 1132 37.0 1116 13.0 1106 11.0 99 1132
    D7769-77 71.4 220 0.32 0.0993 0.0026 2.3484 0.1437 0.1585 0.0079 1611 50.0 1227 43.6 949 43.8 74 949
    D7769-78 119 667 0.18 0.0727 0.0013 1.2233 0.0223 0.1215 0.0010 1006 36.0 811 10.2 739 5.8 90 739
    D7769-79 109 130 0.84 0.0770 0.0016 1.9506 0.0419 0.1833 0.0020 1120 41.2 1099 14.5 1085 11.2 98 1120
    D7769-80 237 216 1.10 0.0781 0.0016 1.9294 0.0394 0.1785 0.0016 1150 40.7 1091 13.7 1059 8.5 96 1150
    D7769-81 210 324 0.65 0.0778 0.0014 1.6667 0.0335 0.1545 0.0016 1143 41.7 996 12.8 926 8.7 92 926
    D7769-82 118 186 0.63 0.0724 0.0015 1.6279 0.0339 0.1623 0.0013 996 42.6 981 13.1 969 7.0 98 969
    D7769-83 165 243 0.68 0.0756 0.0014 1.4400 0.0271 0.1374 0.0011 1085 34.3 906 11.3 830 6.1 91 830
    D7769-84 417 1000 0.42 0.0749 0.0012 1.4251 0.0226 0.1372 0.0010 1066 32.6 899 9.5 829 5.5 91 829
    D7769-85 79.4 408 0.19 0.0643 0.0015 0.6266 0.0167 0.0701 0.0009 750 48.1 494 10.4 437 5.5 87 437
    D7769-86 148 543 0.27 0.0776 0.0013 1.3782 0.0251 0.1280 0.0014 1137 33.0 880 10.7 777 7.7 87 777
    D7769-87 83.1 312 0.27 0.0760 0.0012 1.7573 0.0300 0.1668 0.0016 1094 32.9 1030 11.1 995 9.1 96 995
    D7769-88 49.0 98.8 0.50 0.0716 0.0017 1.4495 0.0373 0.1457 0.0018 976 53.9 910 15.5 877 10.0 96 877
    D7769-89 153 151 1.01 0.0621 0.0017 0.8529 0.0223 0.0995 0.0010 676 63.9 626 12.2 611 5.7 97 611
    D7769-90 448 522 0.86 0.0907 0.0015 2.5088 0.0470 0.1993 0.0022 1440 32.6 1275 13.6 1172 12.0 91 1440
    D7769-91 223 409 0.54 0.0801 0.0015 1.9418 0.0413 0.1743 0.0019 1198 37.0 1096 14.3 1036 10.3 94 1198
    D7769-92 108 299 0.36 0.0765 0.0013 1.4680 0.0250 0.1385 0.0011 1109 34.4 917 10.3 836 6.5 90 836
    D7769-93 146 373 0.39 0.0781 0.0012 1.9649 0.0344 0.1819 0.0022 1150 36.1 1104 11.8 1077 12.0 97 1150
    D7769-94 294 513 0.57 0.0711 0.0015 1.1359 0.0326 0.1158 0.0026 961 45.4 771 15.5 707 14.9 91 707
    D7769-95 132 277 0.47 0.0746 0.0015 1.6725 0.0341 0.1622 0.0012 1058 40.7 998 13.0 969 6.6 97 969
    D7769-96 124 170 0.73 0.0731 0.0015 1.6443 0.0351 0.1630 0.0015 1017 43.1 987 13.5 973 8.2 98 973
    D7769-98 2135 743 2.87 0.2801 0.0042 24.4270 0.4092 0.6307 0.0056 3365 23.5 3286 16.5 3152 22.3 95 3365
    D7769-100 112 231 0.49 0.0683 0.0019 1.2870 0.0438 0.1361 0.0031 877 89.8 840 19.5 823 17.8 97 823
    D7769-101 254 814 0.31 0.0688 0.0013 1.4036 0.0305 0.1470 0.0015 894 39.4 890 12.9 884 8.5 99 884
    D7769-102 185 522 0.36 0.0694 0.0016 1.4180 0.0399 0.1477 0.0026 911 52.8 897 16.8 888 14.9 99 888
    D7769-103 354 1144 0.31 0.0699 0.0014 1.4330 0.0280 0.1477 0.0012 928 41.8 903 11.7 888 6.7 98 888
    D7769-104 142 511 0.28 0.0744 0.0016 1.7229 0.0379 0.1673 0.0013 1051 44.4 1017 14.2 997 7.1 98 997
    D7769-105 401 974 0.41 0.0845 0.0015 2.5123 0.0458 0.2142 0.0017 1303 30.6 1276 13.3 1251 8.9 98 1303
    D7769-106 165 218 0.76 0.0876 0.0020 2.6309 0.0580 0.2163 0.0019 1373 43.7 1309 16.2 1262 10.3 96 1373
    D7769-107 105 188 0.56 0.1078 0.0025 4.4692 0.1049 0.2985 0.0031 1763 42.6 1725 19.5 1684 15.4 97 1763
    下载: 导出CSV 
    | 显示表格
    图  4  吉塘群黑云二长片麻岩锆石U-Pb年龄谐和图
    Figure  4.  Zircon U-Pb concordia diagrams for the biotite monzogneiss from the Jitang Gruop

    对96个年龄值进行统计分析,年龄值频率分布图上形成多个组/峰(图 5),U-Pb年龄值可分为2组:第一组437~203 Ma,锆石有效分析点13个,具有明显的核-边结构,年龄值均测试于锆石边部(图 3),为锆石变质增生边,隐约见环带结构,Th含量为7.75×10-6~325×10-6,U含量为265×10-6~1388×10-6,Th/U值除1号点为0.55外均小于0.1,其年龄代表原岩成岩后期遭受变质作用的时间,13个变质年龄集中分布在5个组/峰值(图 5),分别为226~ 203 Ma(n=4),255~252 Ma(n=3),308~273 Ma(n=3),377~359 Ma(n=2),427 Ma(n=1);第二组2818~566 Ma(图 4-c),锆石有效分析点83个,具核-边或板状结构,且具核-边结构的锆石核部具有环带结构,Th含量为3.5×10-6~2134×10-6,U含量为10×10-6~1752×10-6,Th/U值为0.1~2.87,大多数大于0.4,显示出岩浆成因的锆石特征[18],年龄测试点位于锆石核部;板状结构锆石数量少,无变质边,无环带,外形浑圆,年龄测试点位于锆石中心部位(图 3);该段年龄值代表锆石的形成时间,83个年龄值在相对频率分布图上形成多个密集连续的组/峰(图 5),分别为584~566 Ma(n=2),615~611 Ma(n=3),673~634 Ma(n=6),739~707 Ma(n=5),787~756 Ma(n=4),855~808 Ma(n=12),894~865 Ma(n=8),941~906 Ma(n=11),973~962 Ma(n=5),1048~995 Ma(n=6),1150~ 1100 Ma(n=9),>1198 Ma(n=12),最年轻年龄值为566± 9.51 Ma。

    图  5  黑云二长片麻岩锆石U-Pb年龄相对频数分布图
    Figure  5.  Relative frequency distribution diagram of zircon U-Pb age in the biotite monogenesis

    锆石作为副矿物广泛存在于各类岩石中,具有非常高的封闭温度和稳定性,对锆石进行结构分析和原位微区定年可以提供岩石形成、变质时代等重要信息[19-20]。碎屑沉积岩锆石U-Pb定年的基本原理为:随机选取足够量的碎屑锆石进行单颗粒锆石定年,然后对结果进行统计分析,利用最年轻的一组年龄限定最大沉积时间。利用这种方法确定地层时代的准确度有限,难以获得精准的绝对年龄值,但用于限定无化石或火山岩夹层的地层时代范围非常有效;对于副变质岩,碎屑锆石没有受到后期变质事件影响,同位素体系保持封闭,可以采用这种方法限定其原始沉积时代。

    黑云二长片麻岩样品第二组锆石U-Pb年龄值为2818~566 Ma,测点锆石具有核-边结构或板状结构,且测点均位于锆石的核部或中心。具核-边结构的锆石核部形态显示,它们是不同磨蚀程度的碎屑锆石,CL图像揭示核部为岩浆成因,后期有变质作用生成的变质边;板状结构锆石外部形态表明,其为磨蚀程度不同的碎屑锆石,成岩后期未遭受变质作用。锆石年龄谐和度高(图 4),碎屑锆石U-Pb同位素体系保持封闭,锆石核部或中心部位未受到后期变质事件的影响,可用第二组年龄限定寄主岩石的原始沉积时代。黑云二长片麻岩的碎屑锆石年龄分布在566~2818 Ma较宽的范围,表明其物源丰富多样,既有太古宙古老物质,也有新元古代物质。这些年龄形成多个群组,年龄分布连续,最年轻锆石年龄为566±9.5 Ma,表明吉塘群最大沉积时限不早于566±9.5 Ma,形成时代为新元古代晚期;该年龄值与陈奇等[21]间接推测的多穷沟吉塘群黑云斜长片麻岩年龄560 Ma相近。

    雍永源等[7]认为,吉塘群原岩具有碎屑岩或杂砂岩的特征;陶琰等[10]认为,吉塘群原岩建造中包含火山弧成因的英安岩或英安质火山碎屑岩;陈奇等[21]认为,斜长片麻岩原岩以长英质粘土岩和砂砾岩为主,沉积岩碎屑物质主要来自花岗质岩石、长英质火山岩或类似TTG的源岩;野外地质调查发现,黑云二长片麻岩虽然变质程度较高,且由于受到部分混合岩化和构造改造强烈,原岩组构特征残留很少,但在采样点附近的黑云二长片麻岩中保留了“成分层”(图版Ⅰ-d)。这种“成分层”体现在主要矿物长石、石英、黑云母含量上的不同,部分“成分层”长石矿物减少变成云母片岩。这种“成分层”差异是沉积岩特有的现象。锆石CL图像显示,锆石多呈圆状、次圆状,具有较明显的搬运磨蚀再沉积特征,说明黑云二长片麻岩原岩为碎屑岩;第二段2818~566 Ma的锆石核部具有环带结构,Th/U值均大于0.1,平均0.55,这是岩浆锆石典型的特征,结合岩石矿物成分说明,碎屑物多来源于中—酸性花岗质岩石、长英质火山岩等。

      图版Ⅰ 
    a.黑云二长片麻岩残留于印支期二长花岗岩;b.黑云二长片麻岩手标本照片;c.二长片麻岩弱混合岩化现象;d.黑云二长片麻岩夹石英片岩组成的“成分层”差异;e.黑云二长片麻岩镜下照片;f.石英、长石和云母组成的S-C组构;Pl—斜长石;Qz—石英;Ms—白云母;Bit—黑云母;Kf—钾长石
      图版Ⅰ. 

    陈奇等[21]对吉塘群黑云斜长片麻岩的锆石U-Pb年代学研究认为,吉塘群变质结晶基底源岩由泛非期和格林威尔造山期构造-岩浆热事件形成,在560 Ma的沉积时代之前,源区可能发生了4次岩浆事件;结合本次第二段2818~566 Ma出现多组或峰,锆石为岩浆锆石,可能代表在566±9.51 Ma沉积时代之前物源区中—酸性岩浆活动频繁,多组/峰年龄值分别对应频繁的岩浆活动时间,表明频繁的岩浆活动为后期吉塘群黑云二长片麻岩原岩沉积提供了丰富的物源。

    西藏地质调查院对类乌齐地区吉塘群(报告称类乌齐变质岩套,本文称类乌齐地区吉塘群)在以往基础上进行了解体;西藏地勘局地热地质大队对吉塘地区吉塘群在以往基础上进行了解体,二者的对比结果见表 2

    表  2  藏东吉塘地区吉塘群与类乌齐地区片麻岩套对比
    Table  2.  Comparison of the Jitang Group in Jitang area and Leiwuqi gneiss suite, eastern Tibet
    单元名称 单元代号 岩性组合 原岩恢复 特征变质矿物 变质相带 构造环境 形成时代
    类乌齐吉塘群 宁达片麻岩单元 Nd<gn>T3 浅灰白色花岗质片麻岩、花岗质片糜岩、二长片麻岩、黑云斜长片麻岩、花岗质粗糜棱岩、碎裂花岗岩等 正片麻岩 矽线石/红柱石铁铝榴石/角闪石/黑云母/斜长石 高绿片岩—低角闪岩相 火山弧及同碰撞花岗岩 227±0.6 Ma[8]
    加隆片麻岩单元 JL<gn>T1 眼球状、条纹条带状黑云母(角闪)斜长片麻岩、黑云母花岗质片麻岩、黑云斜长(二长)片麻岩、花岗质片糜岩、碎斑岩、碎粒岩等 正片麻岩 同碰撞花岗岩 246.3±1.7 Ma、248.3±1.5 Ma、247.3±1.3 Ma
    加压玛片麻岩单元 JY<gn>P3 浅灰色、灰白色角闪花岗质片糜岩、黑云斜长片麻岩和浅灰色片麻状黑云二长花岗岩等 正片麻岩 同碰撞花岗岩 253±9 Ma(片麻状花
    岗岩)254±8 Ma[8](片
    麻状花岗岩)、246.3 ±
    0.8[23](片麻状花岗岩)
    都拉扎片麻岩单元 DL<gn>P1 片麻状碎裂细粒二长花岗岩 正片麻岩 同碰撞花岗岩 281.8±1.9 Ma
    查秋赛组 AnOcq. 浅灰色-灰色石榴子石黑云斜长片麻岩,灰色-深灰色条纹条带状钙铝榴石二云斜长片麻岩、浅灰绿色(角闪)斜长片麻岩、黑云斜长片麻岩、黑云二长片麻岩,局部夹少量黑云斜长片岩和矽线石云英岩等 碎屑岩夹火山岩建造 主体区间年龄为438.2±3.6~700 Ma,时代置于前奥陶纪
    吉塘地区吉塘群 纽多岩体 ηγT2/γο 灰白色碎裂岩化细粒二长花岗岩/灰白色碎裂岩化黑云母 低—高角闪岩相带 同碰撞花岗岩 不早于566±9.5 Ma(黑云
    二长片麻岩)
    吉塘复式花岗岩 γσοT3/ξγT3/γοT3/ηγT3 眼球状、条纹条带状、片麻状黑云母二长花岗岩、浅灰色片麻状黑云二长花岗岩、灰白色二长花岗岩/条纹条痕条带状混合岩化英云闪长岩、浅灰色片麻状英云闪长岩、灰白色英云闪长岩/灰色糜棱岩化花岗岩闪长岩、灰白色碎裂花岗岩闪长岩/灰白色条纹条痕弱混合岩化正长花岗岩、灰白色正长花岗岩等 碎屑岩夹火山岩建造 绿泥石/绢云母/绿帘石/角闪石/铁铝榴石/云母斜长石/碱性长石 被动陆源
    吉塘群 Pt3Jabl 斜长片麻岩、二长片麻岩、黑云母(二云)片岩(非正式填图单元) 岛弧玄武岩
    Pt3Jabl 斜长角闪岩 基性火山岩
    下载: 导出CSV 
    | 显示表格

    表 2分析,类乌齐地区吉塘群经解体后可分为两部分,第一部分为形成于印支期的宁达片麻岩单元(227 Ma)、加隆片麻岩单元(246.3 Ma)、加压片麻岩单元(253 Ma)、都拉扎片麻岩单元(281.8 Ma),为同碰撞正片麻岩;吉塘地区吉塘群围岩或附近同样分布印支期花岗岩体,吉塘岩体(225.7~222.8 Ma)、纽多岩体(262~243.6 Ma)形成的大地构造环境为同碰撞—后碰撞。第二部分为查秋赛组,岩石组合为二长(斜长)片麻岩、片岩、云英岩,原岩为碎屑岩夹火山岩,这部分与解体后的吉塘地区吉塘群岩性一致。

    吉塘镇多穷沟地质路线调查发现,印支期岩浆岩具有片麻状构造、条纹(条痕)状构造、眼球状构造及碎裂岩化、糜棱岩化混合岩化(图版Ⅰ-c)等现象,且多发育在岩体与围岩接触带、韧性剪切带及脆-韧性断裂附近,说明吉塘地区在印支期中—酸性岩体形成后,经历了构造-岩浆的改造作用,这与类乌齐地区吉塘群经解体后的正片麻岩相似。这部分不属于吉塘群的变质基底物质,属于华力西期—印支期岩浆活动的产物,应将类乌齐吉塘群中的正片麻岩部分归属于华力西期—印支期岩浆岩。

    吉塘群最大沉积时限不早于566±9.51 Ma,其后年龄记录了变质事件发生的时间;锆石U-Pb变质年龄主要为437~203 Ma,测点位于具核-边结构锆石的变质增生边上,在相对频率分布图上出现多个峰值/组(图 5),代表了不同的变质事件。

    陶琰等[10]认为,吉塘复式花岗岩岩浆源岩为杂砂岩成分的吉塘群副片麻岩。樊炳良等[22]认为,吉塘复式花岗岩源岩可能为吉塘群。西藏地勘局地热地质大队在吉塘复式花岗岩中测得各类中—酸性花岗岩形成时间在225.7~213 Ma范围,为印支期岩浆活动的产物,该年龄与变质组/峰226~ 203 Ma(n=4)在误差范围内一致,为吉塘复式花岗岩岩浆活动引起吉塘群变质事件的时间。吉塘群主要以残留体或捕虏体形式分布于吉塘复式花岗岩中,说明印支期岩浆活动对吉塘群改造程度大,在吉塘复式花岗岩锆石U-Pb年龄值中出现大于吉塘群黑云二长片麻岩沉积下限的年龄566±9.51 Ma,可能与本次岩浆活动捕虏吉塘群锆石关系密切,导致吉塘群重融。

    李才等[8]于恩达附近的片麻状花岗岩中获得锆石U-Pb年龄为254±8 Ma。王新雨等[24]在纽多岩体南部花岗闪长岩中获得锆石U-Pb年龄为262±1.3 Ma。周新等[25]在纽多细粒花岗岩中获得锆石U-Pb年龄为258.1±1.9 Ma。樊炳良等[26]测得纽多二长花岗岩年龄为243.6±1.4 Ma。这些年龄值与变质年龄组/峰255~252 Ma接近,说明印支期岩浆活动引起吉塘群的又一次变质事件的发生。

    变质年龄组/峰273~427 Ma的变质事件目前在吉塘地区尚未报道,但是在类乌齐地区吉塘群中正片麻岩部分相关年龄报道较多。据区域地质调查报告,在吉塘群博日松多岩组二云片麻岩中获取单颗粒锆石U-Pb年龄300±14 Ma;在吉塘群都拉扎片麻岩单元的片麻状碎裂细粒二长花岗岩中获得锆石U-Pb年龄281.8±1.9 Ma。据此认为,308~ 273 Ma(n=3)可能为华力西期岩浆活动对吉塘群产生的变质事件的记录。

    雍永源等[7]在多穷沟酉西群片岩中获得Rb-Sr等时线年龄371.5±0.5 Ma,认为是变质年龄。本次变质年龄组/峰377~359 Ma(n=2),所测锆石(图 3,点12、13)具明显的核-边结构,边部具明显的环带,为岩浆活动引起的锆石生长产生的增生边,可能为岩浆活动引起的吉塘群在华力西期的又一次变质事件的记录。

    河南省地勘局区域地质调查队在侵入吉塘群觉拉片麻岩的片麻状细粒黑云二长花岗岩中获得Rb-Sr等时线年龄438.2±3.6 Ma。该年龄与本次变质年龄组/峰427 Ma(n=1)在误差范围内一致,可能代表吉塘群在加里东期岩浆活动中同样受到了变质作用。

    综上,可分析得出吉塘群原岩成岩后经历了5期岩浆变质事件,变质时期从加里东期不断持续到印支期。其中印支期遭受到2期岩浆变质事件,变质事件年龄峰/组年龄为255~252 Ma、226~203 Ma;华力西期遭受到2期岩浆变质事件,变质事件年龄峰/组年龄为377~359 Ma、308~273 Ma;加里东期遭受到1期岩浆变质事件,变质事件年龄峰/组年龄为427 Ma。其中,印支期岩浆活动对吉塘群改造程度大,可能导致吉塘群重融。

    (1) 吉塘群黑云二长片麻岩原岩为副变质岩,碎屑物质多来源于太古宙、新元古代的中—酸性花岗质岩石、长英质火山岩等,吉塘群最大沉积时限不早于566±9.51 Ma,形成时代为新元古代,在此时间之前岩浆活动频繁,为吉塘群沉积提供了丰富的物质来源。

    (2) 吉塘地区吉塘群与类乌齐吉塘群在物质组成上具有相似性;类乌齐地区吉塘群中正片麻岩单元形成于华力西期—印支期,不属于吉塘群结晶基底物质,应该从吉塘群中分解,归属于华力西期—印支期岩浆活动产生的岩浆岩部分。

    (3) 吉塘群原岩成岩后经历了5期岩浆变质事件,吉塘群沉积后从加里东期—印支期不断遭受到岩浆活动引起的变质,变质事件年龄峰/组为427 Ma、377~359 Ma、308~273 Ma、255~252 Ma、226~203 Ma;其中226~203 Ma岩浆活动对吉塘群改造程度大,可能导致吉塘群重融。

  • 图  1   评价区地质灾害分布图

    Figure  1.   Landslide distribution in study area

    图  2   评价指标证据权重(证据因子,权重值)

    a—坡度;b—坡向;c—断层距离;d—地层岩性;e—土地利用类型;f—绿化灌溉

    Figure  2.   The WoE of indexes

    图  3   综合证据权

    Figure  3.   Integrated WoE

    图  4   后验概率分布图

    Figure  4.   Posterior probability distribution

    图  5   ROC曲线

    Figure  5.   The ROC curve

    图  6   兰州地质灾害易发性评价图

    Figure  6.   Landslide susceptibility map of Lanzhou

    表  1   评价指标及证据因子划分

    Table  1   Assessment indexes and weight factors

    评价指标 证据因子
    序号i 名称 序号j 类别
    1坡度10°~10°
    210°~20°
    320°~30°
    430°~40°
    540°~50°
    650°~60°
    2坡向1342.5°~22.5°(北)
    222.5°~67.5°(东北)
    367.5°~112.5°(东)
    4112.5°~157.5°(东南)
    5157.5°~202.5°(南)
    6202.5°~247.5°(西南)
    7247.5°~292.5°(西)
    8292.5°~342.5°(西北)
    3断层距离10~1 km
    21~2 km
    32~3 km
    43~5 km
    55~10 km
    4地层岩性1前寒武系皋兰群(AnЄgl)黑云片岩、
    片麻岩
    2始新统西柳沟组(E1-2x)砂砾岩、粉砂岩
    3渐新统野狐城组(E3y)粉砂质泥岩、细砂岩、砾岩
    4白垩系河口群(K1hk4)砾岩、泥质砂岩
    5中新统咸水河组(N1x)粉砂质泥岩、砂岩、砂砾岩
    6下更新统(Q1)午城黄土、砂砾石、亚粘土等
    7中更新统(Q3)离石黄土、砂砾石、亚粘土等
    8上更新统(Q2)马兰黄土、砂砾石、亚粘土等
    9全新统(Q4)风积黄土、砂砾石、亚粘土等
    10水域
    11侵入岩
    5土地利用类型1城乡、工矿及居民建设用地
    2草地
    3林业用地与未利用地
    4水利水电、交通建设用地
    5农业用地
    6水域
    6绿化灌溉1未灌溉
    21999年起喷灌区
    32000年起喷灌区
    41999年起漫灌区
    52000年起漫灌区
    下载: 导出CSV

    表  2   皮尔逊相关系数

    Table  2   Pearson correlation coefficient

    相关性
    极强相关 0.8~1.0 −1.0~−0.8
    强相关 0.6~0.8 −0.8~−0.6
    中等程度相关 0.4~0.6 −0.6~−0.4
    弱相关 0.2~0.4 −0.4~−0.2
    极弱相关或无相关 0.0~0.2 −0.2~−0.0
    下载: 导出CSV

    表  3   评价指标相关系数矩阵

    Table  3   The matrix of correlation coefficients for evaluation indicators

    相关
    系数
    坡度 坡向 地层
    岩性
    断层
    距离
    绿化
    灌溉
    土地利用
    类型
    坡度 1.00 0.07 −0.17 0.23 0.09 −0.05
    坡向 0.07 1.00 −0.02 0.01 0.02 0.00
    地层岩性 −0.17 −0.02 1.00 −0.02 −0.20 0.04
    断层距离 0.23 0.01 −0.02 1.00 −0.06 0.04
    绿化灌溉 0.09 0.02 −0.20 −0.06 1.00 0.04
    土地利用类型 −0.05 0.00 0.04 0.04 0.04 1.00
    下载: 导出CSV

    表  4   坡度评价指标证据权计算结果

    Table  4   The WoE calculation process of the slope index

    证据因子序号(j) 1 2 3 4 5 6 计算公式
    证据因子类别 0°~10° 10°~20° 20°~30° 30°~40° 40°~50° 50°~60° 已知
    证据因子栅格数(Aj 231104 146950 172670 95956 16030 330 已知
    证据因子中灾害栅格数(AjAL 101 341 1360 2010 972 25 已知
    评价区栅格数(At 663040 jAj
    评价区灾害栅格数(AL 4809 j(AjAL)
    评价区的非灾害栅格数(¯AL 658231 AtAL
    先验几率(O(D)) 0.00731 AL/(AtAL)
    证据因子中非灾害栅格数(Aj¯AL 231003 146609 171310 93946 15058 305 Aj−(AjAL
    非证据因子中灾害栅格数(¯AjAL 4708 4468 3449 2799 3837 4784 AL−(AjAL
    非证据因子中非灾害栅格数(¯Aj¯AL 427228 511622 486921 564285 643173 657926 AtAj−(AjAL
    正权(W+ −2.81600 −1.14458 0.08308 1.07448 2.17876 2.41763 In(AjALALAj¯AL¯AL)
    负权(W 0.41101 0.17842 −0.03095 −0.38723 −0.20266 −0.00475 In(¯AjALAL¯Aj¯AL¯AL)
    相对系数(C −3.22701 −1.32300 0.11403 1.46171 2.38141 2.42238 W+jWj
    证据权重(W −3.26316 −1.35915 0.07788 1.42556 2.34526 2.38623 Cj+jWj
    下载: 导出CSV

    表  5   地质灾害发生概率与易发性程度分级

    Table  5   Classification of the probability and landslide susceptibility

    概率可能性发生频率易发程度面积比例灾害面积比例频率比
    >0.1非常可能常常会发生高易发3.61%49.59%13.75
    0.01~0.1很可能某些情况下会发生中易发11.89%24.0%2.02
    0.001~0.01有可能极少情况下发生低易发25.07%18.41%0.73
    ≤0.001不可能一般情况下不会发生非易发59.43%8.0%0.13
    下载: 导出CSV
  • Agterberg F P, Bonham−Carter G F, Cheng Q, et al. 1993. Weights of evidence modeling and weighted logistic regression for mineral potential mapping[C]//Davis J C, Herzfeld U C. Computers in geology-25 years of progress. New York: Oxford University Press: 13−32.

    Armaş I. 2012. Weights of evidence method for landslide susceptibility mapping. Prahova Subcarpathians, Romania[J]. Natural Hazards, 60: 937−950. doi: 10.1007/s11069-011-9879-4

    Cheng Q, Agterberg F P. 1999. Fuzzy Weights of Evidence Method and Its Application in Mineral Potential Mapping[J]. Natural Resources Research, 8(1): 27−35. doi: 10.1023/A:1021677510649

    Dahal R K, Hasegawa S, Nonomura A, et al. 2008. GIS−based weights−of−evidence modelling of rainfall−induced landslides in small catchments for landslide susceptibility mapping[J]. Environmental Geology, 54(2): 311−324. doi: 10.1007/s00254-007-0818-3

    Goyes−Penafiel P. 2021. Landslide susceptibility index based on the integration of logistic regression and weights of evidence: A case study in Popayan, Colombia[J]. Engineering Geology, 280: 105958. doi: 10.1016/j.enggeo.2020.105958

    Ilia I, Tsangaratos P. 2016. Applying weight of evidence method and sensitivity analysis to produce a landslide susceptibility map[J]. Landslides, 13(2): 379−397. doi: 10.1007/s10346-015-0576-3

    Kobayashi H, Sanga T. 2005. Application of spatial information technology to landslides Landslide vulnerability mapping using GIS and spatial statistical analysis[J]. Journal of the Japan Landslide Society, 42(4): 281−292. doi: 10.3313/jls.42.4_281

    Lee S, Choi J. 2004. Landslide susceptibility mapping using GIS and the weight−of−evidence model[J]. International Journal of Geographical Information Science, 18(8): 789−814. doi: 10.1080/13658810410001702003

    Neuhaeuser B, Terhorst B. 2007. Landslide susceptibility assessment using "weights−of−evidence" applied to a study area at the Jurassic escarpment (SW−Germany)[J]. Geomorphology, 86(1/2): 12−24.

    Velmurugan J, Venkatesan M. 2018. Soft Bayesian model for landslide risk analysis[J]. International Journal of Soft Computing, 13(2): 31−37.

    Zêzere J L, Ferreira A B, Vieira G, et al. 2000. The use of Bayesian Probability for landslide susceptibility evaluation: a case study in the area north of Lisbon (Portugal) [C]//Bromhead E, Dixon N, Ibsen M L. Landslides in research, theory and practice: Proceedings of the 8th International Symposium on Landslides, Cardiff. London: Thomas Telford: 1635−1640.

    白光顺, 杨雪梅, 朱杰勇, 等. 2022. 基于证据权法的昆明五华区地质灾害易发性评价[J]. 中国地质灾害与防治学报, 33(5): 128−138.
    邓勇, 邱瑞山, 罗鑫, 等. 2007. 基于证据权重法的成矿预测——以广东省钨锡矿的成矿预测为例[J]. 地质通报, 26(9): 1228−1234.
    范海明, 王翔, 茹湘兰. 2017. 基于证据权−奇异性−灰色理论的山西五台地区金矿成矿研究[J]. 地质通报, 36(8): 1462−1466.
    方然可, 刘艳辉, 苏永超, 等. 2021. 基于逻辑回归的四川青川县区域滑坡灾害预警模型[J]. 水文地质工程地质, 48(1): 181−187.
    郭富赟. 2019. 兰州市地质灾害特征与风险管控对策[J]. 城市与减灾, (3): 59−63.
    黎志恒, 张永军, 梁收运. 2014. 兰州城市地质灾害与人类工程活动[J]. 兰州大学学报: 自然科学版, 50(5): 588−593.
    李媛, 孟晖, 董颖, 等. 2004. 中国地质灾害类型及其特征——基于全国县市地质灾害调查成果分析[J]. 中国地质灾害与防治学报, (2): 32−37.
    刘帅, 王涛, 曹佳文, 等. 2024. 基于优化随机森林横型的降雨群发滑坡易发性评价——以西秦岭极端降雨事件为例[J]. 地质通报, 43(6): 958−970.
    饶品增, 曹冉, 蒋卫国. 2017. 基于地理加权回归模型的云南省地质灾害易发性评价[J]. 自然灾害学报, 26(2): 134−143.
    尚瑾瑜, 张永军, 李松. 2017. 兰州市区地质灾害发育特征及类型研究[J]. 甘肃地质, 26(3): 70−77.
    孙岩, 王训练, 陈建平, 等. 2010. 基于证据权重法的沱沱河地区铅锌银矿成矿预测[J]. 地质通报, 29(4): 556−564.
    谭玉敏, 郭栋, 白冰心, 等. 2015. 基于信息量模型的涪陵区地质灾害易发性评价[J]. 地球信息科学学报, 17(12): 1554−1562.
    王立朝, 陈亮, 侯圣山, 等. 2020. 洮河流域下游地质灾害调查成果报告[R]. 中国地质调查局.
    吴明堂, 薛正海, 崔振华, 等. 2022. 基于斜坡单元和证据权−Logistic回归的滑坡易发性评价[J]. 人民长江, 53(10): 87−94.
    许冲, 戴福初, 徐锡伟. 2011. 基于GIS平台与证据权的地震滑坡易发性评价[J]. 地球科学: 中国地质大学学报, 36(6): 1155−1164.
    杨华阳, 许向宁, 杨鸿发. 2020. 基于证据权法的九寨沟地震滑坡危险性评价[J]. 中国地质灾害与防治学报, 31(3): 20−29.
    殷跃平. 2022. 地质灾害风险调查评价方法与应用实践[J]. 中国地质灾害与防治学报, 33(4): 5−6.
    俞布, 潘文卓, 宋健, 等 2012. 杭州市滑坡地质灾害危险性区划与评价[J]. 岩土力学, 33 (S1): 193−199, 216.
    张艳玲. 2009. 汶川地震区地质灾害易发性快速区划与制图[D]. 中国地质大学(北京)硕士学位论文: 19−23.
    赵志芳, 党伟, 王瑞雪, 等. 2010. 基于证据权重法的区域地质环境稳定性遥感评价研究[J]. 国土资源遥感, (S1): 9−13.
    周粤, 王超生, 赵逊, 等. 2022. 怒江支流迪麻洛河流域泥石流易发性评价[J]. 地质通报, 41(4): 702−712.
    自然资源部. 2021. 地质灾害风险调查技术要求(1∶50000)[S].
    自然资源部. 2022. 全国地质灾害防治“十四五”规划[R].
图(6)  /  表(5)
计量
  • 文章访问数:  7829
  • HTML全文浏览量:  6903
  • PDF下载量:  1334
  • 被引次数: 0
出版历程
  • 收稿日期:  2023-02-16
  • 修回日期:  2023-08-18
  • 刊出日期:  2024-07-29

目录

/

返回文章
返回