利用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
相关推荐：
loading