关键词:
地基合成孔径雷达
干涉测量
形变监测
序贯平差
大气补偿
摘要:
地表承载着人类的生活,是人们生产活动的基石。地表稳定性关系着群众的人身安全以及人们生产生活的有序进行。近二十年来,在全球气候变暖的背景下,极端天气的发生频率逐年增多,地质灾害的危害也不断加剧,这也对地表形变监测提出了更高的要求。相对于传统地表形变测量技术,合成孔径雷达(Synthetic Aperture Radar,SAR)因具有全天时、全天候、高分辨率与大范围的优点,成为了地表形变监测的有利手段之一。星载SAR在重访周期以及成像几何等方面有着固有的缺陷,无法满足突发性自然灾害监测的需求。与星载SAR相比,地基合成孔径雷达(Ground-based Synthetic Aperture Radar,GBSAR)因具有高时间分辨率(分钟级至秒级)、高空间分辨率(亚米级)、高精度(亚毫米级)等优点,在精细化实时形变监测中展现出了卓越的性能,成为了星载SAR的一种有效补充。因此,基于地基合成孔径雷达干涉测量(GB-In SAR)技术,积极开展地表形变监测技术研究,对灾害监测预警、抢险救灾应急测绘具有重要的实际意义。目前,地基SAR的数据处理方法大多借鉴了星载In SAR技术。由于地基SAR独特的高时间采样频率特性,使得某些常规的In SAR数据处理技术在地基SAR中无法得到较好应用。此外,目前的地基SAR数据处理方法不能快速应对高时间分辨率带来的大数据量实时处理的情景。鉴于此,本文开展了地基SAR形变监测数据处理技术研究与应用,主要研究内容及创新之处如下:(1)针对利用传统FFT进行距离向成像时出现的频谱泄漏的问题,提出了一种基于全相位FFT的地基SAR距离向成像算法,实现了抑制频谱泄漏的效果,提高了观测目标的辨识度,减弱噪声。该算法将传统的FFT替换为全相位FFT,结合了全相位FFT具有优良的抑制频谱泄漏性能与“初始相位不变”的重要特性,实现最高谱线的初始相位严格等于输入信号的初始相位,保证了最高谱线不受频谱扩散的影响。使用模拟实验与真实实验对该算法进行了验证,与传统FFT相比,全相位FFT的频谱图中,观测目标所对应的主频线更加明显,周围的旁频线得到了有效抑制。在该成像算法下,地基SAR在室内的形变测量精度可达到亚毫米级。(2)针对已有的大气延迟相位补偿方法难以区分地基SAR高频观测信号中的形变与大气延迟变化的问题,提出了一种地基SAR地表形变与大气模型参数联合估计模型,利用形变的时空特征以及大气延迟的空间特征,使用最小二乘实现了形变参数与大气模型参数的联合估计,提高了形变估计的精度。该方法的核心思想是对选取的PS点的地表形变相位和大气延迟相位分别进行建模,再构建干涉相位与地表形变和大气模型参数之间的函数模型与随机模型,并利用最小二乘方法进行参数估计,从而反算求解大气延迟相位。通过设计实验对该算法进行可行性分析,同时为了验证该估计模型的精度优越性,对该联合估计方法与经典的外部气象参数校正法、多项式拟合校正法进行了对比实验,使用均值、标准差以及均方根误差等指标进行精度评估。实验结果证明,本文提出的联合估计方法处理得到的结果标准差和均方根误差相对于其它两种方法都有显著的提高,验证了该联合估计方法的可靠性和精度优越性。此外,将该方法应用在湖南省长沙市实地的大桥振动监测中,成功地将桥梁振动与大气延迟分离。(3)针对地基In SAR后处理效率低,现有的地基In SAR实时处理算法实时性差的问题,提出了一种基于参数可变的序贯平差地基In SAR实时形变函数模型,科学建模了形变参数的动态变化;基于建立的实时形变函数模型及形变参数和大气联合估计方法,构建了一套地基In SAR实时动态处理方法,实现了地基In SAR数据的实时高精度动态处理。使用地基SAR对四川省丹巴县的一处滑坡进行7个小时的连续观测,共计137景SLC影像。使用本文提出的实时动态处理方法对其进行处理,得到了观测期间内滑坡的累积形变结果与平均速率图。从处理结果可以看出,该观测区域有多处滑坡正在发生,累积形变量最大可达20厘米。对实验结果进行分析,参数估计的残余形变的均方根误差90%分布在1mm之内,均值为0.61mm。解缠误差在第一阶段为4%左右;在第二阶段开始,由于观测目标数量突减,严重影响了相位时空三维解缠的精度,解缠误差突增至20%;随着观测目标个数趋于稳定,影像解缠误差稳定在1%左右。随着观测量增多,待求的未知参数也逐渐增多,单个观测耗时也逐渐增长,其中“形变参数解算”环节所耗时间逐渐成为观测总耗时的主要部分,但仍然满足实地观测场景中数据实时处理的需求。(4)针对地基SAR软件可视化程度低,人机交互性弱的需求,研发了地基SAR成像、D-In SAR与PS-In SAR处理软件;开展了地基SAR与高频GNSS桥梁稳定性监测综合实验,成功监测了桥梁的高频振动,验证了本文发展的