[发明专利]一种滑坡作用下长输油气管道易损性评价方法

权利要求书

1.一种滑坡作用下长输油气管道易损性评价方法，其特征在于，包括以下步骤：

S1、对研究区目标管道进行评价单元划分：将目标管道线图层与研究区斜坡单元面图层叠加，切分管段，并导出叠加后的管段图层数据，使每个管段都处于独立的斜坡单元内，然后在导出的管段图层上加载预处理后的高分一号遥感影像图层；

S2、建立管道滑坡易损性评价的指标体系：依次建立一级评价指标体系和二级评价指标体系，获取各管段初始指标量化值，求得管道滑坡易损性评价初始矩阵B，将矩阵B和矩阵B的转置矩阵做各指标相关性分析，得到初始指标两两之间的相关系数ρ；

所述一级评价指标体系包括管道本体类指标和管道与滑坡空间关系类指标；

所述管道本体类指标包括：管道缺陷密度、管道壁厚、内压、管材、管径、输送介质；管道与滑坡空间关系类指标包括：管道埋深、管道位置、管道与滑坡夹角；

所述管道缺陷密度的缺陷点坐标的计算方法为：

假设内检测缺陷点个数为n，管道中线测量坐标点个数为m，以内检测缺陷点P序号i为循环结构外层循环变量，以管道中线测量坐标点Q序号j为循环结构内层循环变量；

当i＝1时，比较缺陷点P₁到起算点Q₁的距离L_1,1与中线测量坐标点Q₂到起算点Q₁的距离l_1,2；

若L_1,1＜l_1,2，则P₁位于Q₁和Q₂之间，根据Q₁和Q₂的坐标值、L_1,1和l_1,2比值，计算P₁坐标；

若L_1,1＞l_1,2，则P₁位于Q₂与Q_j之间，j＞2，计算缺陷点P₁到中线测量坐标点Q₂的距离L_1,2＝L_1,1-l_1,2，比较L_1,2与中线测量坐标点Q₃到起算点Q₂的距离l_2,3的大小；

若L_1,2＜l_2,3，则P₁位于Q₂和Q₃之间，根据Q₂和Q₃的坐标值、L_1,2和l_2,3比值，计算P₁坐标；

若L_1,2＞l_2,3，则P₁位于Q₃与Q_j之间，此时j＞3，计算缺陷点P₁到中线测量坐标点Q₃的距离L_1,3＝L_1,2-l_2,3，比较L_1,3与中线测量坐标点Q₄到起算点Q₃的距离l_3,4的大小，以此类推，直到找到P₁位于Q_j与Q_j+1之间，然后根据Q_j和Q_j+1的坐标值、L_1,j和l_j,j+1比值，计算P₁坐标；

以此类推，直到i＝n时，所有缺陷点坐标计算完毕；

所述管道壁厚的指标值计算公式为：

Ψ_i′＝max(Ψ)+min(Ψ)-Ψ_i 3-1

式中，Ψ_i′—管道易损性评价矩阵中第i条管段的管道壁厚指标值；

Ψ_i—第i条管段的管道壁厚实际值；

max(Ψ)—所有管段实际壁厚值的最大值；

min(Ψ)—所有管段实际壁厚值的最小值；

S3、基于熵权法的管道滑坡易损性评价的熵权计算：将原始数据进行归一化处理，获得第j个评价单元在第i个评价指标上的归一化值r_ij,r_ij∈[0,1]，根据公式，定义第个i评价指标的熵e_i，计算第个i评价指标的熵权w_i,0w_i 1；

S4、管道滑坡易损性评价：将管道滑坡易损性评价的评价原始矩阵B进行归一化处理得到计算管道滑坡易损性值，将各条管段的易损性值植入到管段图层，得到研究区管道滑坡易损性分区结果。