利用2号点的绝对重力值可以计算各点的近似绝对重力值(单位:mGal): 2:979349.1833 6:979349.0836
3:979349.2265 7:979349.2063
4:979349.2474 8:979349.1625
5:979349.2358
沿着测量路径,两点间高差为观测值, 数学模型——观测方程:
5
L1+v1=X5-X4 L2+v2=X6-X5 L3+v3=X3-X6 L4+v4=X2-X3 L5+v5=X4-X2 L6+v6=X3-X4 L7+v7=X7-X3 L8+v8=X6-X7 L9+v9=X8-X6 L10+v10=X7-X8 L11+v11=X8-X7 L12+v12=X2-X8 L13+v13=X7-X2 L14+v14=X3-X7 L15+v15=X5-X3
注意到2号点是已知点,利用上面算出的初值,可以把观测方程写成:V=Bx-l的形式,其中:
随机模型——权阵P使用单位阵。
经过间接平差计算,结果如下(单位:mGal):
点号 2 3 4 5 6 7 8 平差值 改正数 中误差 979349.183 0 0 979349.186 -0.041 0.015 979349.195 -0.053 0.017 979349.188 -0.048 0.018 979349.182 0.099 0.018 979349.182 -0.024 0.015 979349.184 0.021 0.016 2.重力垂直梯度测量
原始数据和计算结果如下表,重力及重力梯度随高度的变化如下图:
6
7
可以看出随着高度的增加,重力大小逐渐减小,而重力垂直梯度基本不变。
3.重力网格测量
重力网格测量数据处理以及各组测量结果汇总见后表。 对于实测数据,我采用了两种归算方法计算重力异常, 空间重力异常:Δg=g+Δ1g-γ0; 布格重力异常:Δg=g+Δ1g+Δ2g-γ0.
其中空间改正Δ1g=0.3086H-0.72X10-7H2(mGal),
层间改正Δ2g=-0.1116H(mGal)
用程序算出的重力异常减去自己算的重力异常得出差值,
对于空间改正,差值的平均值-27.3896mGal,最大最小值-24.8092mGal和-33.0783mGal,标准差2.0683;
对于布格改正,差值的平均值-22.8114mGal,最大最小值-18.9078mGal和-24.0995mGal,标准差0.9794。
可见程序算出的结果和实测归算所得的结果有较大的偏差,估计EGM96模型的精度不足以反映小范围的重力异常变化。
测区的等高线图和重力异常等值线图见图1-图4.
8
相关推荐: