Analysis of the distribution characteristics and influencing factors of geological hazards in Zhejiang Province based on fractal theory
-
摘要: 在浙江省地质灾害空间分布特征研究中引入了多重分形方法,基于分形理论研究地质灾害分布的影响因素。研究结果表明:浙江省地质灾害空间分布具有明显的分形特征,西北和东南地区发育较强烈;浙江省地质灾害空间分布与水系分布、降水量、地形坡度、坡高、地层岩性、道路分布及地形地貌有关,其中在降雨量为2232~2565 mm、坡度为25°~35°范围内、坡高为500~1000 m范围内和侏罗系地层范围内,崩塌、滑坡、泥石流灾害发育明显,在山地地形范围内地质灾害易发生。浙江省的道路分布中,省道分布对地质灾害发育的影响比较大;浙江省钱塘江水系和瓯江水系比较复杂,且对地质灾害发育的影响比较大。研究成果进一步明确了浙江省地质灾害分布特征及主要影响因素影响规律,对灾害防治开展具有指导参考意义。Abstract: The multifractal method is introduced in the study of the spatial distribution characteristics of geological hazards in Zhejiang Province, and the influencing factors of geological hazards distribution are studied based on the fractal theory. The research results show that the spatial distribution of geological hazards in Zhejiang Province has obvious fractal characteristics, and the northwest and southeast regions are more developed; The terrain is related to the topography. Among them, within the range of rainfall of 2232−2565 mm, the slope of 25°−35°, the slope of 500−1000 m, and Jurassic strata, collapse, landslide, and debris flow are obviously developed. Internal geological hazards are prone to occur. In the distribution of roads in Zhejiang Province, the distribution of provincial roads has a greater impact on the development of geological hazards; the Qiantang River and Oujiang River systems in Zhejiang Province are more complex and have a greater impact on the development of geological hazards. The research results further clarified the distribution characteristics of geological hazards in Zhejiang Province and the influence laws of main influencing factors, which have guiding and reference significance for disaster prevention and control.
-
0. 引言
滑坡监测预警是有效减少人员伤亡和财产损失的重要防灾措施之一。近年来随着滑坡监测技术的不断发展,InSAR[1 − 2]、三维激光[3 − 4]、无人机摄影测量[5 − 6]等非接触型边坡监测技术应运而生,可以获取滑坡任意区域的变形数据,但受制于无法进行实时数据分析,一般用于滑坡隐患识别,而对于有实时监测预警需求的滑坡灾害,则主要依靠直接接触的变形监测技术,如裂缝计[7 − 8]、GNSS[9 − 10]等,通过直观可靠的地表变形监测,分析滑坡变形发展趋势进行预警[11 − 12],如近年来广泛应用的切线角预警模型构建的滑坡变形预警技术,成为了实现滑坡实时预警的可靠手段[13]。这类监测技术应用过程中,GNSS设备因其良好的全天候、多方向监测能力得到了广泛的应用,成为滑坡变形监测最主要的设备之一[14]。
GNSS设备应用于滑坡变形监测时,监测精度可达毫米级,常见误差在±5 mm左右[15],但在受天气或设备本身定位解算异常影响时GNSS变形数据也可能出现异常波动,导致根据变形量计算的滑坡变形速率产生明显震荡,引起边坡变形监测的误报。因此,针对监测数据误差问题,一些学者提出了各种方法对变形量监测数据进行过滤,以还原滑坡实际变形速率,包括采用Kalman滤波及其改进[16 − 18]、神经网络法[19 − 22]、小波降噪法及其改进[23 − 25]等,这些方法都能在事后较好地平滑变形数据,并且过滤后的数据趋势也与滑坡实际变形趋势基本相符,但这类方法需要大量的前期数据进行分析,无法真正应用到滑坡的实时预警中。由此,部分学者开始应用具有一定实时数据处理能力的解析式过滤方法,如各类最小二乘和回归拟合[26 − 27]、移动平均法[28]等,以实现变形数据的实时过滤。其中,最小二乘法由于同时考虑了时间序列的影响,相对于移动平均能更真实地还原出滑坡变形发展趋势,得到了更好的应用。
在基于GNSS监测的滑坡变形数据实时过滤过程中,由于GNSS偶然误差的影响,会导致常规最小二乘法过滤时无法有效剔除异常值,导致预警误报的发生,而目前单一的过滤方法都无法实现较好的过滤,因此,解决变形数据的实时过滤问题,并实现异常误差值的实时剔除,是还原真实滑坡变形数据的重要支撑,也是实现滑坡变形真正实时有效预警的关键。
1. GNSS变形数据误差特点
GNSS获取的滑坡变形量误差一般可分为两类,分别是精度误差和偶然误差(图1)。其中,由于设备定位精度影响导致获得的变形量在一定范围内反复波动的误差为精度误差,具有一定的规律和波动幅度;而偶然误差则是由于各类原因导致GNSS设备定位准确度产生的较大偏差,使变形量监测中出现偶发的异常大跳点。在这两类误差影响下,GNSS设备获得的滑坡变形量数据呈现明显的波动,而根据变形量直接计算的滑坡变形速率也呈显著的规律波动和异常抖动,难以进行可靠的变形预警。
2. 变形数据误差过滤方法
2.1 变形数据过滤尺度的确定
对于GNSS设备的变形误差过滤,都是选取一定数量的最新数据,采用各种过滤方法对获得的变形量数据进行实时过滤。而在滑坡变形监测过程中,影响数据过滤效果的最重要因素是数据的选取数量。不同的过滤方法在过滤数据数量不同时,过滤效果也不同。数据选取量越少,平滑后的偏差波动越大;数据选取量越多,平滑效果越好,但对于真正加速后产生的变形量增大则不容易及时识别。目前对于变形量数据的选取数量并没有统一的标准,主要根据经验选取一定数量的变形数据进行平滑,但两者之间应当有一个最佳平衡点,即在确保较好的数据偏差过滤效果的前提下,选取的数据数量最少以减小延迟。同时,还应当考虑不同GNSS设备或区域的差异,构建具有一定普适性的变形数据选取数量,以实现针对这类设备的通用变形数据量选用。
2.2 精度误差的过滤
对于GNSS的变形量精度误差的过滤,目前有很多方法都可以较好地实现数据过滤,如移动平均法、最小二乘法等,其中最小二乘过滤方法简单,拟合后直线斜率为滑坡变形速率,且拟合数据同时考虑了变形量时间序列,在监测频率变化时不会出现局部大幅度波动[29]。因此,本文以常见的最小二乘法过滤为基础,通过搜集国内不同区域、不同型号的GNSS设备获取的边坡不同监测频率的变形量数据进行分析(图2)。通过选取相同长度时间内边坡未产生明显变形的GNSS合位移监测数据,并比较不同数据量选用下的变形量过滤偏差大小的规律,具体方法为:采用不同数量的数据进行最小二乘拟合,获得拟合后的直线位移,并计算其与实际位移量之间的最大偏差值,然后按时间依次向前推进,继续进行最小二乘拟合,并不断计算拟合位移量与实际位移量之间的最大偏差值,最终获得整个时间段上的最大偏移量值。最小二乘拟合的数据量需大于3组,因此分别采用3组~80组数据作为数据选择总量各自进行拟合,最终建立不同数据总量与拟合后最大数据偏差量之间的关系(图3)。
通过对比了甘肃省黑方台陈家3#滑坡(监测频率30 min)、江西省弋阳县三县岭滑坡(监测频率1 h)、四川省理县薛城镇滑坡(监测频率1 h)和云南省丽江市华丽高速边坡(监测频率5 min)4个区域不同监测频率的GNSS监测数据过滤效果的偏差关系发现(图4),整体上GNSS监测数据随着拟合采用的数据总数量增多,拟合后的最大偏差值呈逐渐减小的趋势。不同类型的GNSS监测数据都具有大体相同的减小规律,即整体上随着拟合采用的数据总数量增多,GNSS的最大偏差降低程度可以分为三个阶段:快速降低(数据量3~15)、震荡降低(数据量16~40)、缓慢降低(数据量>40)。在快速降低阶段,随着拟合数据量的不断增多,拟合后的GNSS变形数据最大偏差值呈大幅度下降的特征,数据量越多,偏差降低越明显;在震荡降低阶段,拟合后的GNSS变形数据最大偏差值呈波动下降的特征,整体上随着数据量的增多偏差呈较大幅度降低,同时中间也出现明显的震荡区域,即随着数据量的增多,最大偏差值可能出现小幅度的增大,随后再继续降低;在缓慢降低阶段,随着拟合数据量的继续增多,拟合后的GNSS变形数据最大偏差值呈缓慢下降的特征,如拟合数据总数量增加一倍,从40组数据增加至80组数据时,所有GNSS的最大偏移量降低幅度仅约20%,偏差值的过滤效果并不好。
由此可见,对于常见的GNSS监测数据的过滤,可以以快速降低和震荡降低两阶段过滤效果为参考,选用缓慢降低阶段起点的40个数据量进行数据过滤。为了进一步明确该选取值是否有较好的适用性,通过现有甘肃省黑方台陈家3#滑坡和云南省丽江市华丽高速边坡的监测数据按照每小时1组提取进行修正,采用相同过滤方法对比不同监测频率下相同边坡的数据过滤规律(图5)。通过对比发现,GNSS变形监测数据表现出的过滤阶段与频率无显著关系,不同监测频率下仍然可以采用相同的三阶段进行划分,在大约40个数据量时也处于缓慢降低阶段的起点。而不同频率对GNSS变形监测数据的过滤效果差异主要在于偏差的降低幅度,这是由于监测频率过高使得数据量过于集中,在局部时段内GNSS误差呈整体偏大或偏小时,导致整个数据的偏差过滤都出现一定的偏差。而当监测频率大幅降低时,数据量间的持续时间变长,会克服局部时间段区域内数据整体偏大或偏小的问题,因此总体偏差降低幅度会增加,如图5中云南省丽江市华丽高速边坡GNSS监测频率从5 min降低至1 h时,同样的80组数据对应的时间段也从400 min变成80 h,对应了偏差降低幅度也从60%增加至85%。
根据分析发现,对于GNSS变形监测数据的过滤,可以以缓慢降低阶段起点对应的40个数据量作为过滤数据量尺度,对应变形数据的偏差降低幅度均超过50%,能达到在尽量少的数据量前提下取得较好的过滤效果。
2.3 偶然误差的过滤
由于GNSS变形监测数据通过卫星获取并解算定位信息,在极端恶劣环境或偶然故障情况下,可能出现变形量监测值的大幅度偏差,即偶然误差(图1)。这类误差呈偶发性,通过多个省份多个类型的GNSS数据进行统计分析发现,GNSS均有一定的偶然异常值,且异常数据总量均小于全部数据总量的10%,呈普遍性规律,这部分数据偏移量极大,可能是正偏移值(即出现极大的数据值),也可能是负偏移值(即出现极小的数据值),属于错误数据,本身没有意义。由于偶然误差的出现会导致计算出的变形速率出现显著的震荡,使基于变形速率的预警产生误报,同时常规的过滤方法难以直接过滤这类误差,因此,最好的处理方式是及时剔除偶然误差,但这一过程需要实时条件下进行才能保证滑坡监测预警的实时性。
为了实现监测数据的实时过滤,同时及时判识并剔除偶然误差,可以采用设置一定缓冲过滤区的误差剔除方法(图6)。考虑到GNSS变形监测数据过滤量在40个时,可以较好地实现常见精度误差的过滤,对应缓冲过滤区的数据量也设置为40个,具体过滤方法为:当GNSS变形监测数据获取到最新的变形量数据后,取最新的40组变形量监测数据作为数据缓冲区,由时间正序去掉10%的最大变形量(即从过去到最新时间顺序选取4组最大值),由时间倒序去掉10%的最小变形量数据(即从最新时间到过去顺序选取4组最小值)。采用该方法处理后,数据中的偶然误差极大值和极小值都被剔除,还原出仅有精度波动的有效数据;而剩余32组数据作为有效变形量数据,再采用最小二乘过滤方法进行过滤以获得滑坡变形速率信息。随着GNSS不断获取新的变形量数据,对应更新缓冲过滤区实现偶然误差的实时剔除和精度误差的实时过滤。
该方法的应用中,缓冲过滤区需要40个监测数据作为基数,按照常规GNSS变形监测数据1 h/次的监测频率,对应缓冲过滤区的判断需要连续40个小时的监测数据。需要说明的是,对于具有突发性特点的滑坡该方法难以及时判断其变形发展趋势。对此可以引入具有监测频率动态调节的自适应GNSS监测设备,在滑坡加速变形时通过自适应变频实现监测频率的自动增加。具体为:通过设置GNSS变形监测设备阈值为S(S大于设备监测精度),采用5 min/次的动态监测对比方式进行判断,当实际监测数据和最新记录的GNSS数据之差大于阈值S时,认为变形可能存在加速情况,此时记录最新的监测数据;当实际监测数据和最新记录的GNSS数据之差小于阈值S时,认为变形尚不明显,仍然采用1 h/次的常规监测频率获取数据。以此不断进行比对和记录,实现滑坡变形在进入加速变形期间能自动将监测频率调整到最高5 min/次,对应40组变形监测数据的最小时间间隔仅为不足4小时,对于不具有强烈突发性的滑坡,可以及时有效地识别出加速且不会导致误报。
3. 过滤方法应用效果探讨
四村滑坡位于四川省茂县黑虎乡,为降雨诱发的堆积体变形滑坡,该区域近年来已实施了自动化监测,并布设了1套GNSS监测站。通过前几年的监测获取了一段时间的GNSS变形监测数据(图7a),该滑坡整体处于基本稳定状态。但在2018年4月12日到13日出现了一次较明显的加速过程(图7b),GNSS监测站全程获取了该次加速变形的位移变化数据。通过提取该加速变形区间的GNSS合位移,并分别采用最小二乘法和剔除误差后的最小二乘法拟合该区域的变形速率进行对比发现(图7c),变形速率的变化特征在剔除误差前后基本一致,能有效反应出GNSS合位移先增大后减缓时表现出的滑坡变形速率增大和减小,两者之间一致性较好,并未出现明显的延迟或显著差异。
同时,以最新的彭州市某矿山边坡GNSS变形监测数据为例,进一步验证该方法的过滤效果。该矿山边坡为矿山开采后形成的人工边坡,目前较为稳定,未产生明显变形。GNSS监测站获取了边坡2024年1月份的变形数据(图8a),可见虽然边坡无明显变形,但获取的GNSS原始累计合位移数据有明显的异常跳跃点(图8b),导致了根据累计位移计算的变形速率呈上下波动状态(图8c),最大变形速率达到283.81 mm/d,显然不符合实际情况。采用数据过滤方法处理,剔除10%最大数据和10%最小数据后采用最小二乘拟合,得到的变形速率较稳定地在0附近(图8c),实时反映出了边坡的真实变化状态。
可见,结合了数据过滤和偶然误差剔除的GNSS变形数据过滤方法,可以在实时监测的情况下去掉大偏差数据,且不影响监测数据的实际发展趋势判定,两者的结合为滑坡变形的及时预警提供可靠的数据过滤算法。
4. 结论
针对滑坡地表变形常用的GNSS监测技术在预警过程中出现的精度误差和偶然误差,本文分析了常见的GNSS设备监测获得的滑坡变形数据误差特征,提出了对应的数据过滤方法,以实现滑坡变形数据的实时过滤,为有效的提前预警提供了科学数据过滤方法。本文主要得到以下结论:
(1)滑坡GNSS变形数据存在精度误差和偶然误差,其中精度误差可以通过数据解析式过滤的方法实时过滤,偶然误差则应通过剔除的方式进行过滤。
(2)以最小二乘法进行解析式过滤时发现,随着参与过滤的数据总量增多,GNSS过滤后的变形量数据偏差大小呈快速降低、震荡降低、缓慢降低三个阶段,可以选取缓慢降低阶段起点对应的数据量值40组数据确定为过滤尺度,在确保过滤效果的前提下减少数据延迟。
(3)通过构建数据缓冲区剔除GNSS变形数据可能出现的偶然误差,同时采用自适应变频技术的方式,使构建数据缓冲区导致的数据分析滞后时间减小到最短不足4小时,为滑坡的实时预警预报提供及时可靠的数据支持,满足非突发性滑坡基于变形速率的预警需求。
(4)通过对典型滑坡GNSS变形数据的应用验证,该方法可以实现变形速率误差波动的有效过滤,同时偶然误差的剔除也对滑坡实际变形速率不产生明显延迟,可以较好地还原滑坡实际变形发展趋势,具有一定的普适性和适用性。
-
表 1 各类地质灾害点分布的拟合函数及分维值
Table 1 Fitting function and fractal dimension value of the spatial distribution of various disaster types
灾害类型 拟合函数 分维值 相关系数 总地质灾害 y = −1.6924x + 10.484 1.6924 0.9996 滑坡 y = −1.646x + 10.248 1.6460 0.9994 崩塌 y = −1.4984x + 9.4427 1.4984 0.9973 泥石流 y = −1.4132x + 9.0765 1.4132 0.9948 表 2 2010—2018年浙江省各市地质灾害发生次数
Table 2 The number of geological hazards in cities in Zhejiang Province from 2010 to 2018
灾害类型
地区滑坡/次 崩塌/次 泥石流/次 总计/次 杭州市 231 280 95 606 宁波市 239 192 62 493 温州市 567 114 74 755 嘉兴市 49 46 12 107 湖州市 34 117 10 161 绍兴市 203 151 29 383 金华市 69 47 9 125 衢州市 367 114 45 529 舟山市 21 22 0 43 台州市 24 19 4 47 丽水市 375 241 76 692 表 3 不同水系地质灾害分布分维值
Table 3 Fractal dimension value of the distribution of geological hazards in different water systems
分维值
八大水系水系 地质灾害 分维值 相关系数 分维值 相关系数 苕溪水系 1.3818 0.9437 1.2089 0.9496 运河水系 1.1684 0.9368 0 0 甬江水系 1.2089 0.9923 0.8919 0.9225 钱塘江水系 1.4686 0.9911 1.4204 0.9896 椒江水系 1.0689 0.9722 0.9489 0.9827 瓯江水系 1.4012 0.9941 1.3569 0.9942 飞云江水系 0.9832 0.9394 0.9610 0.9428 鳌江水系 0.7807 0.9347 0.7621 0.9389 表 4 不同降雨量范围内地质灾害分布分维值
Table 4 Fractal dimension value of geological hazard distribution in different rainfall ranges
地质灾害
降雨量/mm分维值 崩塌 滑坡 泥石流 1231~1565 0 0.6229 0.4937 1565~1898 1.0411 1.2517 0.92 1898~2232 1.284 1.449 1.2211 2232~2565 1.3402 1.5138 1.3203 表 5 不同坡度范围内地质灾害分布分维值
Table 5 Fractal dimension value of geological hazards distribution in different slope ranges
地质灾害
坡度/(°)分维值 崩塌 滑坡 泥石流 ≤8 0.6281 0.9113 0.5231 8~15 1.0259 1.2558 1.0099 15~25 1.1265 1.3646 1.0983 25~35 1.1817 1.4418 1.1548 >35 0.9045 1.0422 1.0401 表 6 不同海拔高度范围内地质灾害分布分维值
Table 6 Fractal dimension value of geological hazards distribution in different altitude ranges
地质灾害
坡高/m分维值 崩塌 滑坡 泥石流 0~300 0.7585 1.0925 0.7176 300~500 1.1460 1.3961 1.1182 500~1000 1.2716 1.4482 1.1844 1000~1500 1.0448 1.2646 1.0630 表 7 不同地层岩性分布分维值
Table 7 Fractal dimension value of different strata lithology distribution
分维值
地层分维值 相关系数 中新生代 第四系 1.5007 0.9995 古近系+新近系 1.1900 0.9806 白垩系 1.4483 0.9977 侏罗系 1.6684 0.9999 泥盆-三叠系 三叠系 1.4223 0.9968 二叠系 1.3683 0.9875 石炭系 1.1159 0.9751 泥盆系 1.1590 0.9964 震旦志留系 志留系 0.9884 0.9899 奥陶系 1.0771 0.9935 寒武系 1.2494 0.9963 震旦系 0.8793 0.9900 表 8 不同道路分布分维值
Table 8 The distribution fractal dimension value of the different road
类型 拟合函数 分维值 相关系数 高速 y = −1.4894x + 9.4108 1.4894 0.9993 铁路 y= −1.2661x + 8.2485 1.2661 0.9978 国道 y = −1.4065x + 8.9952 1.4065 0.9947 省道 y = −1.5106x + 9.5783 1.5106 0.9967 全部道路 y = −1.6433x + 10.232 1.6433 0.9993 表 9 不同地形地貌范围内地质灾害分布分维值
Table 9 Fractal dimension value of geological hazards distribution in different topographic and landform ranges
地质灾害
地形地貌分维值 崩塌 滑坡 泥石流 总地质灾害 山地 1.4001 1.6284 1.3874 1.6465 丘陵 0.5138 0.9154 0.4837 0.9306 台地 0.2211 0.7585 0.117 0.6496 平原 0.9666 1.2423 0.7925 1.2922 洼地 − − − − -
[1] 潘懋, 李铁锋. 灾害地质学[M]. 北京: 北京大学出版社, 2002 PAN Mao, LI Tiefeng. Disaster geology[M]. Beijing: Peking University Press, 2002. (in Chinese)
[2] 段晨玉,张以晨,张继权,等. 基于GIS的地质灾害风险预警系统设计[J]. 灾害学,2022,37(1):199 − 204. [DUAN Chenyu,ZHANG Yichen,ZHANG Jiquan,et al. Design of geological disaster risk early warning system based on GIS[J]. Journal of Catastrophology,2022,37(1):199 − 204. (in Chinese with English abstract) DOI: 10.3969/j.issn.1000-811X.2022.01.033 [3] 江南. 盒维数概念的历史演变[J]. 自然辩证法研究,2020,36(10):84 − 90. [JIANG Nan. The historical evolution of the concept of box dimension[J]. Studies in Dialectics of Nature,2020,36(10):84 − 90. (in Chinese with English abstract) DOI: 10.19484/j.cnki.1000-8934.2020.10.015 [4] 周琰嵩. 云南省嵩明县地质灾害发育规律与致灾因素分析[D]. 北京: 中国地质大学(北京), 2018 ZHOU Yansong. Analysis on the development law and disaster factors of geological disasters in Songming County, Yunnan Province[D]. Beijing: China University of Geosciences, 2018. (in Chinese with English abstract)
[5] 菊春燕,贾永刚,潘玉英,等. 基于分形理论的旅游景区地质灾害危险性评估—以青岛崂山为例[J]. 自然灾害学报,2013,22(6):85 − 95. [JU Chunyan,JIA Yonggang,PAN Yuying,et al. Geohazard assessment of tourist attractions based on fractal theory:A case study of Qingdao Laoshan[J]. Journal of Natural Disasters,2013,22(6):85 − 95. (in Chinese with English abstract) [6] 桂蕾,殷坤龙. 基于多重分形理论的滑坡地表监测位移分析[J]. 中南大学学报(自然科学版),2014,45(11):3908 − 3914. [GUI Lei,YIN Kunlong. Analysis of landslide surface monitoring displacement based on multi-fractal theory[J]. Journal of Central South University (Science and Technology),2014,45(11):3908 − 3914. (in Chinese with English abstract) [7] 姜恩三,任光明,陈锦涛. 四川省普格县地质灾害发育的分形特征[J]. 中国地质灾害与防治学报,2016,27(3):122 − 126. [JIANG Ensan,REN Guangming,CHEN Jintao. The fractal characteristics of the geological hazards development of Pugh County Sichuan Province[J]. The Chinese Journal of Geological Hazard and Control,2016,27(3):122 − 126. (in Chinese with English abstract) DOI: 10.16031/j.cnki.issn.1003-8035.2016.03.18 [8] 王森,许强,罗博宇,等. 基于分形理论的南江县滑坡敏感性分析与易发性评价[J]. 水文地质工程地质,2017,44(3):119 − 126. [WANG Sen,XU Qiang,LUO Boyu,et al. Vulnerability analysis and susceptibility evaluation of landslides based on fractal theory in Nanjiang County[J]. Hydrogeology & Engineering Geology,2017,44(3):119 − 126. (in Chinese with English abstract) [9] 王潇,刘婷,姚昆,等. 分形理论在达州市达川区滑坡灾害研究中的应用[J]. 人民长江,2019,50(2):144 − 150. [WANG Xiao,LIU Ting,YAO Kun,et al. Application of fractal theory to landslide disaster research in Dachuan District,Dazhou City[J]. Yangtze River,2019,50(2):144 − 150. (in Chinese with English abstract) DOI: 10.16232/j.cnki.1001-4179.2019.02.025 [10] 乔朋腾. 吕梁山区滑坡崩塌发育规律及易发性研究[D]. 西安: 长安大学, 2019 QIAO Pengteng. Study on development regularity and susceptibility of landslide and collapse in Lv-Liang Mountains[D]. Xi’an: Changan University, 2019. (in Chinese with English abstract)
[11] 陶洁. 基于分形理论和显著性分析的滑坡提取研究[D]. 北京: 中国地质大学(北京), 2020 TAO Jie. Study on extraction of landslide based on fractal theory and saliency analysis[D]. Beijing: China University of Geosciences, 2020. (in Chinese with English abstract)
[12] 都来瀚. 基于分形理论的岩质滑坡失稳判据与形变预测研究[D]. 北京: 中国地质大学(北京), 2020 DU Laihan. Study about instability criterion and deformation prediction of rock slope based on fractal theory[D]. Beijing: China University of Geosciences, 2020. (in Chinese with English abstract)
[13] 史兴旺,管新邦,边筠. 云南省滑坡灾害分形特征研究与分形评价模型探讨[J]. 自然灾害学报,2021,30(3):209 − 216. [SHI Xingwang,GUAN Xinbang,BIAN Yun. Study on fractal characteristics and fractal evaluation model of landslide hazards in Yunnan Province[J]. Journal of Natural Disasters,2021,30(3):209 − 216. (in Chinese with English abstract) [14] 孔庆鸽. 分形插值函数的盒维数及应用[D]. 杭州: 浙江大学, 2018 KONG Qingge. Box dimension and applications of fractal interpolation functions[D]. Hangzhou: Zhejiang University, 2018. (in Chinese with English abstract)
[15] 曹俊贤,宋春林. 分形维数的计算及改进[J]. 信息技术与信息化,2017(10):19 − 23. [CAO Junxian,SONG Chunlin. Calculation and improvement of the fractal dimension[J]. Information Technology and Informatization,2017(10):19 − 23. (in Chinese with English abstract) DOI: 10.3969/j.issn.1672-9528.2017.10.001 [16] 张磊,陈张建,黄丽. 浙江省地质灾害防治管理平台设计与实现[J]. 中国地质灾害与防治学报,2020,31(2):102 − 110. [ZHANG Lei,CHEN Zhangjian,HUANG Li. Design and implementation of Zhejiang Province geological disaster prevention management platform[J]. The Chinese Journal of Geological Hazard and Control,2020,31(2):102 − 110. (in Chinese with English abstract) DOI: 10.16031/j.cnki.issn.1003-8035.2020.02.14 [17] 浙江省自然资源与环境统计年鉴(2019)[M]. 北京: 中国统计出版社, 2019 Zhejiang Province natural resources and environment statistical yearbook (2019)[M]. Beijing: China Statistics Press, 2019. (in Chinese)
[18] 祝杨菲,梁勤欧,林德根. 浙江省滑坡地质灾害风险地图制图综合研究[J]. 国土与自然资源研究,2019(3):34 − 41. [ZHU Yangfei,LIANG Qinou,LIN Degen. Research on mapping of landslide geological hazard risk map in Zhejiang region,China[J]. Territory & Natural Resources Study,2019(3):34 − 41. (in Chinese with English abstract) DOI: 10.3969/j.issn.1003-7853.2019.03.010 [19] 冯杭建,周爱国,唐小明,等. 浙江省泥石流灾害发育分布规律及区域预报[J]. 地球科学,2016,41(12):2088 − 2099. [FENG Hangjian,ZHOU Aiguo,TANG Xiaoming,et al. Development and distribution characteristics of debris flow in Zhejiang Province and its regional forecast[J]. Earth Science,2016,41(12):2088 − 2099. (in Chinese with English abstract) [20] 王若谷. 露天煤矿排土场地貌多重分形特征及其与岩土侵蚀关系[D]. 北京: 中国地质大学(北京), 2019 WANG Ruogu. The relationship between soil erosion and multi-fractal parameters of landform in an opencast coal-mine dump in a loess area of China[D]. Beijing: China University of Geosciences, 2019. (in Chinese with English abstract)