第四章重力资料处理与分析.PPTVIP

第四章重力资料处理与分析 按最简单的假设δ与地球平均密度ρ、地球半径R、地表重力加速度g0及地球刚性系数μ存在如下关系： h=5f/(2f+1) （2-2-2） k=3f/(2f+1) （2-2-3） 其中f=ρRg0/15μ。 地球模型的不同，f值也不同。故借助δ的异常变化，再结合倾斜潮汐因子和应变潮汐因子的变化，可反算勒夫数h和k，研究地壳构造运动，为地震预报和地球动力学研究服务。 维尼迪科夫调和分析计算步骤 1、计算每个48段的中央时刻T的儒略世纪数,并计算六个天文元素； 2、根据天文天文元素和给定的波群分类数，计算T时刻波群内各波的理论相位和理论振幅； 3、对48小时（两整天的观测数据）读数组使用了一组滤波器，计算各相应波振幅与相位，即观测值；（数字滤波与最小二乘法相结合） 4、进行误差分析； 5、求解潮汐因子与相位差。 现在台站使用常规重力观测仪器主要为阿斯卡尼亚GS仪器和我自己研制的DZW微伽重力仪，其它仪器主用于实验观测。 非潮汐观测 非潮汐观测结果主要指在重力固体潮观测结果中扣除潮汐效应后的部分。通过适当的方法排除其中的仪器漂移影响后，可以获得测点位置的重力长期变化结果。重力长期变化包括极移、大气、地下水、物质迁移、点位移动、地震孕育等信息，认真加以研究和鉴别，可以提取与地震预报有关的前兆信息。 固体潮观测曲线包括的成份: 1、固体潮汐 2、外界干扰 3、仪器零漂 4、地球内部物质、密度等变化（地震） 外界干扰: 温度变化干扰及校正 气压改正（对仪器影响、对重力场影响） 地下水干扰的校正（地下水类型） 人为因素造成的重力效应(大型工程、抽水造成地面下降等 其它 固体潮观测资料的零漂的计算 日均值法 25小时均值法 Φ0=（y0+ y±1+ y±2+ y±3+…+ y±12)/25 彼尔采夫法 Φ0=（y0+ y±2+ y±3+ y±5 +y±8 + y±10 + y±13+ y±18)/15 重力场非潮汐变化 一、引力常数G变化可能引起重力非潮汐变化 二、地球的核幔边界变化可能引起的非潮汐变化 三、区域性的短期非潮汐变化 非潮汐观测结果主要指在重力固体潮观测结果中扣除潮汐效应后的部分。通过适当的方法排除其中的仪器漂移影响后，可以获得测点位置的重力长期变化结果。重力长期变化包括极移、大气、地下水、物质迁移、点位移动、地震孕育等信息，认真加以研究和鉴别，可以提取与地震预报有关的前兆信息。 课堂讨论1、讨论大仓点重力变化原因？能给我们什么启示？2、讨论重力非潮汐变化的原因？说明研究重力非潮汐变化意义与困难之处。3、台站重力观测曲线包含内容？资料处理程序，怎样提取有用信息？ 形变观测数据流程 “十五”处理软件功能 “十五”观测数据产出流程 数据预处理 年初归零：每年1月1日0时起数据归零。 突跳处理：对严重影响动态曲线的单点突跳数据做补插处理。 缺数处理：将因停电、仪器标定等造成的缺记、非正常数据标记为“缺数”。 台阶处理：对因标定、仪器调零等造成的观测数据序列突变掉格进行处理以形成连续观测值。 预处理记录：对所有预处理过程记录、标注其操作内容和特征参数。 彼尔采夫组合计算仪器零漂 Φ=1/15(y0+y±2+y±3+y±5+y±8+y±10+y±13+y±18) 重力台站观测成果的质量评价 一、资料的连续性 二、观测资料内符合精度 1、nakai检验均方差 2、潮汐因子的稳定性 台站观测数据质量评价——Nakai方差检验 原理：概率统计原理（大概率事件与小概率事件） 1当观测值g（t)落于[Y（t) -2δ、Y（t) +2 δ]区间的概率为95% 2当残差︱ν （t) ︱≥ 2δ的概率为5%为小概率事件 Nakai方差检验过程 1、输入台站经度纬度及观测时间 2、计算48整点时刻的固体潮理论值及其导数，组合系数矩阵A，M 3、输入重力观测值g（t)，形成列向量B； 4、解方程X= 5、计算各时刻拟合值、偏差值及各段数据的方均差 值δ 6、判别：若有某个观测值残差满足≥ 2δ，则用Y（t)代替 g（t)，再重复4-6步运算，直至所有的观测值都满足︱ν （t) ︱≤ 2δ为止。 讨论：台站观测数据的质量 1、方均差 值δ大小问题 结论：Nakai处理方法只能恢复观测数据的原有精度，而不能提高数据的精度。数据的质量是由仪器的观测精度和观测工作水平所决定的。 2、数据质量评价： 当η1和η2均在90%以上者为优 3、潮汐因子及潮汐因子的中误差 计算每月和年度重力潮汐因子δ和潮汐因子的中误差mδ。按月统计M2、O