单站高频GNSS求解同震位移的新方法及联合强震仪的地震预警应用研究

【摘要】：地震是地球上最为常见的自然灾害之一,其中极具破坏力的大地震带来的原生灾害、次生灾害和衍生灾害给人类造成了巨大的威胁。虽然无法控制和避免地震的发生,但是人类也在积极地寻求应对之策。目前,在准确的地震短期预测还无法实现的情况下,地震预警便是一种行之有效的防震减灾措施。历经三十年的探索和努力,很多国家和地区建立了以地震仪/强震仪为基础的地震预警系统,并且取得了众多成功预警的实例。但是,这样的地震预警系统在应对大地震时(M7)还存在问题,尤其在地震震级估算上容易“失灵”。这是由于地震震级估算与震时地表位移密切相关,受基线偏差和量程限幅的影响,地/强震仪很难准确把其测量值转化为与地震过程相吻合的地面位移值。因此,如何快速、准确获取震时地表位移成为了地震预警系统可靠运行必须解决的关键问题。近年来,随着高频GNSS技术(高采样率数据,存储技术等)的快速发展,其已能用于监测短期瞬时的地表运动状态。由于GNSS不会受到类似于地/强震仪器记录限幅和倾斜的影响,因而即使在大地震中,也可以较好地获取地震近场同震位移,为地震预警提供了一种新的、可靠的观测手段。虽然高频GNSS在获取震时地表位移上具有优势,但是相对于地/强震仪而言,它在实时地震监测和预警应用中还存在一些限制因素和挑战：(1)GNSS的直接观测量是伪距、载波相位和多普勒观测值,准确的测站位置和速度信息需要通过高精度的数据处理算法才能获得,而且这些算法还需具备实时性；(2)GNSS的测量噪声较大,对微弱地表震动不敏感,不易拾取到初至地震波；(3)GNSS受限于接收机的频带宽度,无法获得与地/强震仪相当的采样率(200 Hz甚至400 Hz)。因此,高频GNSS和地/强震仪具有明显的优势互补性。联合两者将有助于提高地震预警的准确性和时效性,这也是当前地震学的研究热点和前沿课题。本论文旨在将高频GNSS技术引入到现行地震预警系统中,以增强其在面对大地震时的预警能力。围绕这一核心目标,本文主要从数据处理的角度出发,深入研究并提出了单站GNSS单历元求解测站同震位移的新方法,详细推导并建立了GNSS和强震仪组合模型,实现了两者联合进行地震预警应用的目的。论文的主要工作和贡献如下：(1)系统研究了常用的两种利用单站GNSS求解同震位移的方法：精密单点定位(PPP)和单站测速方法。首先介绍了各方法的基本原理,然后建立了相应的求解同震位移模型,最后结合具体地震事件分析了各方法的应用效果。深入分析了接收机钟跳对单站测速法的不利影响,提出了一种顾及接收机钟跳影响的多普勒观测值构造方法,并用动态和静态算例验证了所提方法的可行性和有效性。(2)以单站GNSS实时准确获取同震位移为目标,并结合地震应用的实际需求(获取一段时间的同震位移,并且震前测站坐标已知),从GNSS基本观测模型出发,提出了一种直接求解测站同震位移的新方法——(附有已知坐标的)历元间位置差法。根据数据处理策略的不同,历元间位置差法又可细分成TPP法和R-variometric法。从函数模型的角度出发,证明了PPP法,TPP法和R-variometric法的等价性。详细分析了该方法中所包含的各种主要误差对位移的影响,并给出相应的解决策略。在顾及各项误差改正后,采用历元间位置差法模型能够获得与事后PPP法精度相当的同震位移。实例结果表明：在20分钟的时间里,该法所获得同震位移的精度在水平方向上约为3厘米,在垂直方向上约为5厘米。(3)针对地震监测和预警需要高密度的站网需求以及双频GNSS接收机成本高昂等问题,提出了一种电离层空域建模模型——MSEID模型,将历元间位置差法的适用范围从双频接收机拓展至单频接收机,从而使单频接收机能够高精度的获取同震位移并用于大范围加密现有的双频GNSS观测台网。详细分析和讨论了MSEID模型恢复电离层延迟的有效性,所提模型获得位移的准确性以及适用性。实例结果表明：MSEID模型能够足够精确的恢复出电离层延迟变化情况,并且比SEID模型的计算效率更高；所提模型存在一定的适用范围,对于网内加密,双频测站间距最好小于80公里,对于网外加密,单频测站距离双频测站的平均间距最好在50公里范围内；以双频观测数据计算的同震位移作为参照,所提模型能够获得准确的同震位移,精度在平面方向为2厘米,高程方向上为4厘米。(4)详细对比和分析了强震仪和高频GNSS两种观测手段在获取同震位移时的优缺点。为了克服单一仪器的不足,提出了顾及基线偏差改正的GNSS和强震仪组合模型。实例验证表明：通过Kalman滤波算法能够将两种观测仪器的优势充分结合起来,充分发挥其各自优势,获得了更为准确的宽频带同震位移。在滤波的早期阶段具有更多强震仪结果的特性,可用于探测测站或地面微弱地运动；而在滤波的中后期,则受到更多GNSS位移结果的约束,使得组合位移不会发散,因而能够准确地记录地震波的低频信息。(5)开展了联合高频GNSS和强震仪在地震预警中的应用研究工作,详细探讨并分析了基于组合位移结果进行地震波到时拾取、地震震中位置确定以及地震震级估算的方法及其效果。以东日本大地震为例,联合使用STA/LTA方法和AIC准则进行地震波P波到时拾取,所得结果与理论到时平均时延约0.8秒；采用Geiger法和网格搜索法进行地震震中确定,能在震后18.52秒,分别获得17公里和4公里偏差的震中位置,相应的地震发震时刻偏差在0.5s范围内；根据峰值位移(Pd与PGD)与地震震级的经验关系进行震级估算和预警信息播报,在震后26.35秒可进行第一次预警播报,相应震级为Mw 7.4,在震后102.19秒进行最后一次播报,可获得Mw 8.9的预警最终震级。上述结果表明基于组合的方式,能够获得足够精确的位移结果,可用于P波探测,能够在较短时间内(102.19秒)获得与真实震级十分接近的最终震级而未发生震级饱和现象,这对于增强现有地震预警系统面对大地震时的能力具有极其重要的作用。