[发明专利]一种交通线沿线区域的滑坡风险定量评估方法

权利要求书

1.一种交通线沿线区域的滑坡风险定量评估方法，其特征在于：包括以下步骤：

S1：依据目标区域内的历史记录资料和实地观测资料编制天然边坡的滑坡清查图，需要包括滑坡位置、发生时间、滑坡形态特征及其参数、滑坡体积的信息，通过ArcGIS将这些滑坡绘制在比例尺为1：10000的数字地形图中；

S2：获取滑坡规模概率，依据滑坡清查图统计滑坡体积在10000立方米以上的滑坡发生的频率，把该频率作为滑坡规模概率PM；

S3：从交通线沿线有关部门收集降雨量数据，结合滑坡清查图，将引发滑坡的前5天内的最小降雨量作为降雨量临界值，把首次在天然边坡上引发山体滑坡的降雨量临界值的超越概率作为滑坡发生的时间概率PT,PT用下式计算：

P_T＝P{(R＞R_T)∩L}＝P{R＞R_T}·P{L|R＞R_T}

其中，R为降雨量，RT为降雨量临界值，

P{R＞R_T}为泊松模型确定的降雨量临界值的超越概率，P{L|R＞R_T}为降雨量超过临界值后滑坡发生的概率，将时间段划分为1年、3年、5年、15年、25年和50年；

S4：获取交通线区域的滑坡空间概率PS；

S401：依据滑坡清查图把滑坡坡向、边坡坡度、区域土地使用情况和区域风化层厚度作为Logistic回归分析的敏感性模型的四个准备因子；使用ArcGIS对描述交通线区域的数字地面模型(DTM)进行插值，每间隔5°将滑坡坡度进一步分类；使用ArcGIS从DTM中提取滑坡坡向，每间隔30°对滑坡坡向进一步分类；按照森林、灌木、茶园、荒地等不同使用情况对区域分类；采用反距离加权的点内插技术，在交通线区域内布置一定数量的测点揭露地层的风化层厚度，使用地质统计学软件R进行球面拟合，得到区域内的风化层厚度分布图；

S402：在ArcGIS中选择基于像素的制图单元，并选择10m×10m的像素作为敏感性模型的基本映射单元；

S403：把四类准备因子和滑坡清查图作为敏感性模型的输入参数且作为一个训练数据集，通过10m×10m的像素对所有输入数据栅格化，转换为等面积的网格，对敏感性模型进行迭代计算；计算完成后由接收器操作特征(ROC)曲线判断训练数据集的性能，用模型输出的估计参数作为每个像素的空间概率值；

将区域划分为若干个以排水线、山脊线和坡度断裂为边界的边坡单元，把边坡单元内像素的最大分布空间概率作为边坡单元的空间概率，得到交通线区域的空间概率分布；

S5：认为滑坡规模概率PM、滑坡时间概率PT、滑坡空间概率PS相互独立，滑坡的危险性评估为这三个概率的联合概率，即滑坡危险性H_L＝P_M×P_T×P_S，即可得到未来1年、3年、5年、15年、25年和50年时间段内交通线区域内的滑坡危险性定量评估情况。

2.根据权利要求1所述的一种交通线沿线区域的滑坡风险定量评估方法，其特征在于：将滑坡危险性定量评估分解为三个独立的基本定量参数，即滑坡规模概率PM、滑坡时间概率PT、滑坡空间概率PS。

3.根据权利要求1所述的一种交通线沿线区域的滑坡风险定量评估方法，其特征在于：计算交通线区域的滑坡空间概率PS时运用了Logistic回归分析。

4.根据权利要求1所述的一种交通线沿线区域的滑坡风险定量评估方法，其特征在于：滑坡危险性H_L＝P_M×P_T×P_S，可得到未来1年、3年、5年、15年、25年和50年时间段内交通线区域内的滑坡危险性定量评估情况。