一种燃气管道重点区域地质灾害SAR动态监测方法及装置与流程

本发明属于地质灾害识别与监测技术领域，尤其涉及一种燃气管道重点区域地质灾害sar动态监测方法及装置。

背景技术：

随着我国的经济建设发展，燃气管网已成为国家重要的基础设施和民生工程，是现代能源体系和现代综合交通运输体系的重要组成部分，在国民经济中占有十分重要的战略地位。已建成管网和规划管网不可避免地要穿越高山、丘陵及平原地区，管网所经之地可能大量发育着各种类型的地质灾害。燃气管网一旦被线路两侧的滑坡、泥石流等地质灾害损坏，就会造成巨大的经济损失和人员伤亡。由此，开展燃气管道区域地质灾害识别与监测工作显得尤为重要。

燃气管道区域地质灾害主要表现为滑坡、崩塌、泥石流和地面塌陷等，传统的地质灾害调查以人工线路巡查为主，周期长，成本高，调查成果质量也难以保证。燃气管道覆盖区域较大，单点地质灾害调查难以全面准确的反应区域地质灾害整体情况，在区域地质灾害动态监测方面有一定的应用局限性。

近年来，以insar为主的遥感新技术凭借其全天时、全天候、远距离获取地表微小形变的能力，已经广泛应用于地质灾害的前期大范围识别与后期高精度监测的研究与应用中。在利用insar对滑坡等灾害识别的过程中，常遇到相干目标较少、相干性较差等不能有效获取形变信息的情况，现有的技术中虽提出了多源遥感调查方法、光学遥感+insar融合调查方法，但适用性存在一定的限制，例如没有考虑区域地质灾害形成条件与影响因素，光学遥感+insar方法对于植被覆盖密集的地区，无法提取有效的形变信息；多源遥感调查方法需要在无人机设备、光学遥感数据的支撑下实现，对硬件和数据要求较高；多源遥感调查方法中堆栈insar技术解算结果精度不够高。

技术实现要素：

本发明旨在解决上述问题，提供一种燃气管道重点区域地质灾害sar动态监测方法及装置。

第一方面，本发明提供一种燃气管道重点区域地质灾害sar动态监测方法，包括：从燃气管道重点区域获取k幅原始格式的sar数据，生成单视复数影像；将单视复数影像中的高分辨率影像sar数据采用dinsar两轨法提取影像覆盖时间范围内形变信息；

将单视复数影像中的中等分辨率影像sar数据采用psinsar和sbas技术提取影像覆盖时间范围内形变信息；

将中等分辨率和高分辨率影像sar数据提取的形变信息进行形变结果解译，进行地质灾害形成条件分析；

在地质灾害形成条件分析基础上再选择坡度、坡高、坡形、岩土类型、植被指数、水系影响和人类工程活动作为影响因子，基于gis的加权信息量模型进行地质灾害易发性评价。

进一步，本发明所述燃气管道重点区域地质灾害sar动态监测方法，所述采用dinsar两轨法提取影像覆盖时间范围内形变信息，包括步骤：影像配准-干涉组合生成干涉图-dem编码-差分干涉图生成-自适应滤波和趋势性相位去除-相位解缠-垂向形变量解算-地理编码，从而提取影像覆盖时间范围内垂向形变量空间分布信息；

所述采用psinsar和sbas技术提取影像覆盖影像时间范围内形变信息，包括步骤：dem编码-影像配准-小基线原则生成干涉图-差分干涉图生成-自适应滤波和去除趋势相位-提取相干目标上的差分干涉相位-回归分析提取时序形变信息，从而获得影像时间范围内垂向形变速率空间分布。

进一步，本发明所述燃气管道重点区域地质灾害sar动态监测方法，所述形变结果解译包括点位时空变化、燃气管道沿线沉降剖面图、统计图表信息分析形变时空分布特征。

进一步，本发明所述燃气管道重点区域地质灾害sar动态监测方法，所述点位时空变化、燃气管道沿线沉降剖面图、统计图表信息分析形变时空分布特征，具体包括：

利用地理编码表将某些点位的地理坐标转换为雷达坐标，获取点位的行列数，进行回归分析得到每个点位的时间序列变化；

将沉降量空间分布图进行空间差值，提取燃气管道沿线的沉降值，生成剖面图；将沉降量和沉降速率结果按照像元或相干目标的值进行形变等级划分，并统计不同等级的占比，以及空间分布占比；

结合点位时间序列变化、沉降剖面图和统计图表，总结燃气管道沿线的形变时空分布特征。

进一步，本发明所述燃气管道重点区域地质灾害sar动态监测方法，所述基于gis的加权信息量模型进行地质灾害易发性评价包括：采用栅格单元进行易发程度划分；在地质灾害形成条件分析基础上再选择坡度、坡高、坡形、岩土类型、植被指数、水系影响和人类工程活动作为影响因子，在arcgis平台下，基于gis的信息量分析模型实现地质灾害易发性评价，圈定地质灾害易发程度，进行燃气管道重点区域地质灾害早期识别和动态监测。

进一步，本发明所述燃气管道重点区域地质灾害sar动态监测方法，所述基于gis的加权信息量模型进行地质灾害易发性评价包括以下步骤：1)计算单因素(指标)xi提供斜坡失稳(a)的信息量i(xi,a)：

其中：s-已知样本总单元数；n-已知样本中变形破坏的单元总数；si-有xi的单元个数；ni-有指标xi的变形破坏单元个数；

2)计算某一单元p种因素组合情况下，提供斜坡变形破坏的信息量ii，即：

3)根据ii的大小，给单元确定稳定性等级；

其中ii＜0表示该单元变形破坏的可能性小于区域平均变形破坏的可能性；

ii＝0表示该单元变形破坏的可能性等于区域平均变形破坏的可能性；

ii＞0表示该单元变形破坏的可能性大于区域平均变形破坏的可能性；

即单元信息量值越大越有利于斜坡变形破坏；

4)经统计分析找出突变点作为分界点，将区域分成极高易发区、高易发区、中易发区和低易发区四个不同等级。

第二方面，本发明提供一种燃气管道重点区域地质灾害sar动态监测装置，包括：采集模块，用于从燃气管道重点区域获取k幅原始格式的sar数据，生成单视复数影像；

形变提取模块，用于将单视复数影像中的高分辨率影像sar数据采用dinsar两轨法提取影像覆盖时间范围内形变信息；将单视复数影像中的中等分辨率影像sar数据采用psinsar和sbas技术提取影像覆盖时间范围内形变信息；

解译模块，用于将中等分辨率和高分辨率影像sar数据提取的形变信息进行形变结果解译，进行地质灾害形成条件分析；

评价模块，用于在地质灾害形成条件分析基础上再选择坡度、坡高、坡形、岩土类型、植被指数、水系影响和人类工程活动作为影响因子，基于gis的加权信息量模型进行地质灾害易发性评价。

第三方面，本发明提供一种燃气管道重点区域地质灾害sar动态监测电子装置，包括存储器和处理器；所述存储器，用于存储计算机程序；所述处理器，用于当执行所述计算机程序时，实现如第一方面所述的燃气管道重点区域地质灾害sar动态监测方法。

第四方面，本发明提供一种计算机可读存储介质，所述存储介质上存储有计算机程序，当所述计算机程序被处理器执行时，实现如第一方面所述的燃气管道重点区域地质灾害sar动态监测方法。

本发明所述燃气管道重点区域地质灾害sar动态监测方法及装置，实现了insar技术在燃气管道重点区域动态形变信息的提取，弥补了原有调查方法对数据和设备要求严苛的不足；针对不同类型不同数据量的sar数据，采取不同类型的insar技术方法，针对不同的研究区设置合理的模型参数，提高了insar监测的精度；将insar时序技术得到的平均形变速率和坡度、坡高、坡形、岩土类型、植被指数、水系影响和人类工程活动等共同构成地质灾害易发性评价的影响因子，对燃气管道重点区域的地质灾害程度的划定更加合理准确。

附图说明

图1为本发明实施例所述燃气管道重点区域地质灾害sar动态监测方法流程示意图；

图2为本发明实施例所述采用dinsar两轨法提取影像覆盖时间范围内形变信息流程示意图；

图3为本发明实施例所述燃气管道重点区域地质灾害sar动态监测装置结构示意图。

具体实施方式

下面通过附图及实施例对本发明所述燃气管道重点区域地质灾害sar动态监测方法及装置进行详细说明。

实施例说明本发明的实施方式，熟悉此技术的人士可由本说明书所揭露的内容轻易地了解本发明的其他优点及功效，显然，所描述的实施例是本发明一部分实施例，而不是全部的实施例。

实施例一

在本公开实施例中，所公开的燃气管道重点区域地质灾害sar动态监测方法，如图1、图2所示，包括以下步骤：

步骤一：从中等分辨率和高分辨率数据融合的insar形变信息提取，结合dinsar、psinsar和sbas技术，获取形变信息。

从燃气管道重点区域获取k幅原始格式的sar数据后，生成单视复数影像；对于数据量较少的高分辨率数据采用dinsar两轨法，基于外部dem，经过影像配准-干涉组合生成干涉图-dem编码-差分干涉图生成-自适应滤波和趋势性相位去除-相位解缠-垂向形变量解算-地理编码，从而提取影像覆盖时间范围内垂向形变量空间分布信息；针对数据量较大的中等分辨率sar影像，融合psinsar和sbas关键技术，包含：影像配准-dem编码-小基线原则组合干涉图-差分干涉图生成-滤波和去除趋势相位-提取相干目标-回归分析提取时序形变信息；中等分辨率和高分辨率影像提取的形变结果在形变趋势上具有一致性，进行下一步的形变结果解译。

步骤二：形变结果解译：点位时空变化、燃气管道沿线沉降剖面图、统计图表信息分析形变时空分布特征。

利用地理编码表将点位的地理坐标转换为雷达坐标，获取点位的行列数，进行回归分析得到每个点位的时间序列变化；此外，将沉降量空间分布图进行空间差值，提取燃气管道沿线的沉降值，生成剖面图；根据沉降最大最小值划分几个等级，将沉降量和沉降速率结果按照像元或相干目标的值进行形变等级划分，并统计不同等级的占比，以及空间分布占比；然后，结合点位时间序列变化、沉降剖面图和统计图表，总结燃气管道沿线的形变时空分布特征，包括燃气管道沿线有无沉降快速变化区、沉降较大区、时间域有无沉降加速的情况等。

步骤三：使用基于gis的加权信息量模型对重点区域进行地质灾害易发性分析，圈定燃气管道重点区域地质灾害易发程度。

首先采用栅格单元进行易发程度划分；采用栅格单元的优点是可利用gis实现单元的快速剖分，同时栅格数据为矩阵形式，可借助计算机快速完成运算；在重点区域地质灾害形成条件分析的基础上，选择坡度、坡高、坡形、岩土类型、植被指数、水系影响和人类工程活动，再结合利用本实施例步骤二中时序insar技术得到的沉降速率共同组成地质灾害易发性分析的影响因子，并对影响因子进行重分类；在arcgis平台下，利用基于gis的信息量分析模型实现地质灾害易发性评价，圈定地质灾害易发程度，实现燃气管道重点区域地质灾害早期识别和动态监测，具体包括以下步骤：

1)计算单因素(指标)xi提供斜坡失稳(a)的信息量i(xi,a)：

其中：s—已知样本总单元数；n—已知样本中变形破坏的单元总数；si—有xi的单元个数；ni—有指标xi的变形破坏单元个数。

2)计算某一单元p种因素组合情况下，提供斜坡变形破坏的信息量ii，即：

3)根据ii的大小，给单元确定稳定性等级。ii＜0表示该单元变形破坏的可能性小于区域平均变形破坏的可能性；ii＝0表示该单元变形破坏的可能性等于区域平均变形破坏的可能性；ii＞0表示该单元变形破坏的可能性大于区域平均变形破坏的可能性；即单元信息量值越大越有利于斜坡变形破坏。

4)经统计分析找出突变点作为分界点，将区域分成极高易发区、高易发区、中易发区和低易发区等四个不同等级。

实施例二

本公开另一具体实施例所述燃气管道重点区域地质灾害sar动态监测装置，如图3所示，包括：采集模块，用于从燃气管道重点区域获取k幅原始格式的sar数据，生成单视复数影像。

形变提取模块，用于将单视复数影像中的高分辨率影像sar数据采用dinsar两轨法提取影像覆盖时间范围内形变信息；将单视复数影像中的中等分辨率影像sar数据采用psinsar和sbas技术提取影像覆盖时间范围内形变信息。

解译模块，用于将中等分辨率和高分辨率影像sar数据提取的形变信息进行形变结果解译，进行地质灾害形成条件分析。

评价模块，用于在地质灾害形成条件分析基础上再选择坡度、坡高、坡形、岩土类型、植被指数、水系影响和人类工程活动作为影响因子，基于gis的加权信息量模型进行地质灾害易发性评价。

实施例三

本公开另一具体实施例所述燃气管道重点区域地质灾害sar动态监测电子装置，包括存储器和处理器；所述存储器，用于存储计算机程序；所述处理器，用于当执行所述计算机程序时，使得计算机执行上述实施例一中的燃气管道重点区域地质灾害sar动态监测方法，具体提取步骤与前述实施例一相同，不再赘述。

实施例四

本公开另一具体实施例计算机可读介质，所述计算机可读介质存储有程序代码，当所述计算机程序代码在计算机上运行时，使计算机执行上述实施例一中的燃气管道重点区域地质灾害sar动态监测方法，具体提取步骤与前述实施例一相同，不再赘述。

实施例所述计算机可以是通用计算机、专用计算机、计算机网络、或者其他可编程装置。所述计算机指令可以存储在计算机可读存储介质中，或者从一个计算机可读存储介质向另一个计算机可读存储介质传输。所述计算机可读存储介质可以是计算机能够读取的任何可用介质或者是包含一个或多个可用介质集成的服务器、数据中心等数据存储设备。所述可用介质可以是磁性介质，(例如，软盘、硬盘、磁带)、光介质(例如，数字通用光盘(dvd))或者半导体介质(例如，固态硬盘(ssd))等。计算机存储代码所形成软件可以位于随机存储器，闪存、只读存储器，可编程只读存储器或者电可擦写可编程存储器、寄存器等本领域成熟的存储介质中。

在本申请各实施例中的各功能单元可以集成在一个处理单元中，也可以是各个单元单独物理存在，也可以两个或两个以上单元集成在一个单元中。在上述实施例中，可以全部或部分地通过软件、硬件、固件或者其任意组合来实现。当使用软件实现时，可以全部或部分地以计算机程序产品的形式实现。所述计算机程序产品包括一个或多个计算机指令。在计算机上加载和执行所述计算机程序指令时，全部或部分地产生按照本申请实施例所述的流程或功能。

以上所述，仅为本申请的具体实施方式，但本申请的保护范围并不局限于此，任何熟悉本技术领域的技术人员在本申请揭露的技术范围内，可轻易想到变化或替换，都应涵盖在本申请的保护范围之内。因此，本申请的保护范围应以所述权利要求的保护范围为准。