一种黄土滑坡灾害综合监测方法及系统

技术特征：
1.一种黄土滑坡灾害综合监测方法，其特征在于，包括以下步骤：获取待监测区域的各项监测数据；根据各项监测数据，利用gnss相对定位技术计算待监测区域的地表三维形变时间序列；利用gnss折射遥感技术计算待监测区域的大气水汽含量时间序列；利用gnss反射遥感技术计算待监测区域周围的地表土壤含水率时间序列；获取的待监测区域的降雨量数据，分析其与大气水汽含量时间序列之间的响应关系，得到大气水汽含量与降雨量之间的关联性；分析三维地表形变时间序列与地表土壤含水率时间序列之间的响应关系，得到待监测区域的滑坡变形与地表土壤含水率之间的关联性；根据待监测区域的大气水汽含量与降雨量之间的关联性，以及解算的待监测区域大气水汽含量，对待监测区域进行降雨预警；根据待监测区域的滑坡变形与地表土壤含水率之间的关联性，以及解译的待监测区域的地表土壤含水率，对待监测区域进行滑坡综合监测。2.根据权利要求1所述的一种黄土滑坡灾害综合监测方法，其特征在于：所述获取待监测区域的各项监测数据，包括：在所述待监测区域设置gnss监测基准站、流动站，并获取所述待监测区域的gnss测站观测数据；获取所述待监测区域的广播星历数据；获取所述待监测区域的精密星历数据；获取所述待监测区域的气象数据。3.根据权利要求2所述的一种黄土滑坡灾害综合监测方法，其特征在于：所述根据各项监测数据，利用gnss相对定位技术计算待监测区域的地表三维形变时间序列，包括以下步骤：根据所述gnss测站观测数据，计算卫星的位置、速度、钟差；根据所述gnss测站观测数据，利用伪距单点定位计算流动站、基准站的概略坐标；根据所述gnss测站观测数据，计算所述基准站与流动站的非差残差项；对所述gnss测站观测数据进行潮汐改正、对流层改正、天线相位中心改正；对所述基准站与流动站的非差残差项进行卡尔曼滤波；利用周跳探测中的m-w组合观测算法，并根据基准站与流动站的概略坐标，以及经过卡尔曼滤波处理后的非差残差项，确定整周模糊度；提取流动站相对基准站的站心坐标neu；根据流动站相对基准站的站心坐标neu，得到待监测区域的地表三维形变时间序列。4.根据权利要求3所述的一种黄土滑坡灾害综合监测方法，其特征在于：所述根据各项监测数据，利用gnss折射遥感技术计算待监测区域的大气水汽含量时间序列，包括：计算gnss测站观测数据中的天顶方向总延迟ztd；从所述气象数据中获取待监测区域的气压、温度，并计算大气加权平均温度；计算gnss测站观测数据中的天顶干延迟zhd；计算湿延迟zwd，其计算公式为：zwd＝ztd-zhd
ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
(1)根据湿延迟zwd，以及大气加权平均温度，计算不同时刻的大气水汽含量pwv，得到待监
测区域的大气水汽含量时间序列。5.根据权利要求1所述的一种黄土滑坡灾害综合监测方法，其特征在于：所述根据各项监测数据，利用gnss反射遥感技术计算待监测区域周围的地表土壤含水率时间序列，包括：提取所述gnss测站观测数据中的信噪比snr数据；选取观测卫星的高度角，并设置方位角，截取所述信噪比snr数据中发生高频振荡的信噪比snr数据；利用低阶多项式将发生高频振荡的信噪比snr数据中的直、反射信号分离；对信噪比snr数据中分离出的反射信号进行重采样；对重采样后的所述反射信号进行ls频谱分析，得到gnss测站天线相位中心到地面的反射高h；将所述反射高h、反射信号的波长λ、信噪比snr数据以及卫星高度角θ代入公式(2)，并利用非线性最小二乘法拟合得到相位参数ψ和幅值a；根据所获取的相位参数ψ，建立线性反演模型获取待监测区域周围的地表土壤含水率时间序列。6.根据权利要求1所述的一种黄土滑坡灾害综合监测方法，其特征在于：所述分析大气水汽含量时间序列与降雨量数据之间的响应关系，得到大气水汽含量与降雨量之间的关联性，包括：绘制所述大气水汽含量时间序列和降雨量数据之间的双坐标轴关系图，从图中分析大气水汽含量时间序列和降雨量数据之间的关系，得到大气水汽含量与降雨量之间的关联性。7.根据权利要求1所述的一种黄土滑坡灾害综合监测方法，其特征在于：所述分析三维地表形变时间序列与地表土壤含水率时间序列之间的响应关系，得到待监测区域的滑坡变形与地表土壤含水率之间的关联性，包括：根据所述三维地表形变时间序列，得到滑坡累积位移；绘制所述滑坡累计位移与地表土壤含水率之间的双坐标轴关系图，得到所述滑坡累计位移与地表土壤含水率之间的关联性。8.根据权利要求7所述的一种黄土滑坡灾害综合监测方法，其特征在于：在所述绘制滑坡累计位移与地表土壤含水率之间的双坐标轴关系图时，通过时滞互相关分析方法定量分析滑坡位移速率对地表土壤含水率变化响应的滞后时间；依据所述滑坡位移速率对地表土壤含水率变化响应的滞后时间，绘制所述滑坡位移速率与地表土壤含水率之间的双坐标轴关系图，从图中分析所述滑坡位移速率与地表土壤含水率之间的相关性。9.一种黄土滑坡灾害综合监测系统，其特征在于：包括：数据获取模块，用于获取待监测区域的各项监测数据；数据解算模块，用于根据各项监测数据，利用gnss相对定位技术计算待监测区域的地表三维形变时间序列；利用gnss折射遥感技术计算待监测区域的大气水汽含量时间序列；利用gnss反射遥感技术计算待监测区域周围的地表土壤含水率时间序列；第一响应关系获取模块，用于获取的待监测区域的降雨量数据，并分析其与大气水汽
含量时间序列之间的响应关系，得到大气水汽含量与降雨量之间的关联性；第二响应关系获取模块，用于分析三维地表形变时间序列与地表土壤含水率时间序列之间的响应关系，得到待监测区域的滑坡变形与地表土壤含水率之间的关联性；综合监测模块，用于根据待监测区域的大气水汽含量与降雨量之间的关联性，以及解算的待监测区域的大气水汽含量，对待监测区域进行降雨预警；根据待监测区域的滑坡变形与地表土壤含水率之间的关联性，以及解译的待监测区域的地表土壤含水率，对待监测区域进行滑坡综合监测。

技术总结
本发明公开了一种黄土滑坡灾害综合监测方法及系统，其黄土滑坡灾害综合监测方法包括以下步骤：获取待监测区域的各项监测数据；利用GNSS相对定位技术及GNSS遥感技术，获取待监测区域的地表三维形变时间序列、大气水汽含量时间序列、地表土壤含水率时间序列；分析大气水汽含量时间序列与待监测区域的降雨量数据之间的关联性；分析三维地表形变时间序列与地表土壤含水率时间序列之间的关联性；对待监测区域进行降雨预警，及对待监测区域进行滑坡综合监测。本发明利用GNSS技术获取待监测区域的形变信息、大气水汽含量和土壤含水率的变化等环境信息，并分析其相互之间的响应关系，实现将GNSS技术更好地服务于滑坡灾害的综合监测。将GNSS技术更好地服务于滑坡灾害的综合监测。将GNSS技术更好地服务于滑坡灾害的综合监测。


技术研发人员：张双成 周昕 田静 刘奇 樊茜佑 马中民
受保护的技术使用者：长安大学
技术研发日：2022.08.15
技术公布日：2022/11/15