基于GNSS网主成分分析的断层形变检测和反演.pdfVIP

中国地球科学联合学术年会 2014 ·1571 · 基于 GNSS 网主成分分析的断层形变检测与反演 * 徐克科 伍吉仓 同济大学测绘与地理信息学院 上海 200092 摘 要：基于OKADA 弹性位错理论，利用GNSS 网位移时空数据主成分析，研究了一套集断层形变检 测与滑移时空分布反演于一体的理论方法。针对主成分分解时因 GNSS 缺失数据并由此导致矩阵中心 化时产生的偏差问题，采用了附加偏差参数加权低秩分解的方法。根据断层形变高空间相关性的特点， 通过主成分时空响应分析，揭示了断层活动方式及滑移特征。并以此作为约束条件构建了主成分反演模 型，实现了断层滑移时空分布的快速反演。以2005 年苏门答腊 8.6 级大地震震后滑移为例，通过地表 GNSS 站的主成分分析，检测了断层形变特征，精细地反演了断层滑移时空分布。 关键词： GNSS 网；断层形变；滑移分布；主成分分析 1．引言 GNSS 位移时空序列低信噪比制约了活动断层微动态形变信息的提取，且从庞大的GNSS 观测网络 大量台站中检测断层形变十分耗时与不易。目前，基于连续 GNSS 观测数据进行断层形变的研究主要 集中于两方面。一是通过GNSS 数据时空相关噪声处理，提高板块整体速度场的解算精度。如Dong 等 （2008 ）提出PCA 和KLE 相结合，提取了南加州GNSS 基准站2000-2004 年的区域共性误差[1-3] 。二 是利用地球物理反演的方法反演断层滑移时空分布，如 Segall、McGuire （2003 ）和Kositsky （2010 ） [4] [5] 基于OKADA 断层位错模型，分别通过卡尔曼滤波反演 和主成分反演 的方法研究了断层滑移分布。 但这种反演必须事先提供较为准确的断层几何参数和滑移特征作为先验信息。近年来，Kang Hyeun Ji(2013)等利用滤波和PCA 方法检测到了赫克托耳矿山地震和地下水引起的地壳瞬态形变[6]，但没有进 行反演工作。本文以活动断层带为研究区域，研究了一套基于 GNSS 网时空数据主成分分析的断层形 变检测与滑移时空分布反演于一体的理论方法。首先对地表 GNSS 站位移时空序列去除长期趋势项、 年/半年周期项、阶跃等，再进行PCA 和KLE 联合滤波，剔除共模误差的影响。根据断层形变高空间 相关性的特点，利用主成分时空响应分析得出了不同断层活动方式和演化过程对地表测站位移造成的不 同时空影响规律，揭示了断层活动方式及滑移特征。以此作为约束条件构建了主成分反演模型，实现了 断层滑移时空分布的快速反演。 2 ．主成分时空响应分析方法 对于一个共有 站，观测了 天的区域网GNSS台站位移时空序列，组成矩阵为X 。其中，每 m n m×n 一行代表一个给定站点北（或东或高程）方向所有历元的值；每一列表示在一个给定的历元所有站点的 某个分量位移值。为保证所有站点位移都是以零为基点，按（1）式对矩阵实施中心化， i, j 分别表示 矩阵的行列号。设： m ∑X 0 (i ,k ) k 1