同震重力

基于GRACE重力卫星数据研究同震形变
一、研究背景

二、国内外进展
GRACE观测用于大地震同展及展后变形的研究现状

ＧＰＳ观测用于日本地震同震及震后变形的研究现状

三、理论与方法
1.利用GRACE时变重力数据提取同震重力
1)GRACE数据处理
GRACE卫星的数据处理和存储主要由美国德克萨斯大学空间研究中也德国地学中也GFZ、UTCSR及美国宇航局喷气推进实验室JPL合作完成。而GRACE数据产品根据处理的程度分为Leve1-0、Leve1-1A、Leve1-1B和Leve1-2不同类型。Leve1-2数据中的月重力场模型是研究时变重为场的基础，提供时变地球重为场模型。GRACE在项的不足，利用SLR数据解算得到的项替换。GRACE月重力场模型的阶数越大（尤其是在90阶Ｗ后），其误差也越大，因此根据具体问题选取适当阶数的重力场模型非常关键。一般选择60阶的重力场模型。现在已发布了较多的重力场模型，因此选取精度最高的重力场模型。通常利用内符合精度和外符合精度来评定重力场模型的精度。外符合精度和内符合精度评定公式分别如下：

式中、为模型系数与参考模型系数之差，、为模型系数对应阶次的误差方差，为地球赤道半径。比较外符合精度时，采用EGM2008模型。

利用GRACE时变重力场模型得到的重力变化可以表示为：

R为地球的平均半径，r= a(1−fcos2θ)为地面点到地球质心的距离。θ和λ为地理余纬和地理经度，ΔCnm, ΔSnm为真实重力场与正常重力场的完全规格化球谐系数之差。Pnm为规格化勒让德函数。n、m分别为重力场模型的阶数和次数。由球谐系数得到的重力场变化的网格大小为：

GRACE时变重力场模型由于球谐系数自身的相关性及高阶噪声大的叠频问题，获取的时变信号中存在明显的南北向异常条带。为消除这些影响可以进行高斯平滑滤波和去相关滤波。
2)提取同震重力场变化
地震震前孕育、同震释放和震后调整等不同阶段对应着不同的物质运移过程，具有不同的研究价值。震前重力变化的研究能为地展孕震机理研究及地震预报震例实践提供有益的构造性观测资稱；同震重力变化研究有利于验证重力卫星提取构造性重力变化方法的可行性及有效性；震后变化分析能为研究震中区的地下结构等提供观测资料。因此如何有效分离出地震事件在不同发生阶段的重力变化，具有很重要的现实意义。[11]
有效分离同震重力的方法有：差分法、最小二乘法、经验正交函数法。再此介绍最小二乘法。
为有效分离出震前、同震和震后重力效应，可对震源区周边的各点位重力变化时间序列进行分析，建立涉及震前、同展、震后相关参数的分段函数，从而实现不同阶段的重力信号的有效分离。最直观的半年及年度周期性信号，需要通过周期函数给出；考虑到潮汐频率误差是重力时变