李江俊

(吕梁学院汾阳师范分校，山西 汾阳 032200)

时至今日，对重力学的研究已有400年的历史[1]。有学者在16世纪初就发现地球存在重力。有学者在17世纪发现，重力在不同的位置存在一定的差别。在上述研究理论基础下，有学者提出了地球重力变化产生的原因：地球两极扁平，促使重力由两极向赤道增大，从而导致不同位置出现不同程度的重力，并提出水准面概念[2]。20世纪初，众多科学家投入重力研究中，并不断完善重力网，采用地球卫星采集地球重力数据，健全地球重力模型[3]。随着地球重力模型的完善，有学者针对地球重力异常部分进行了研究，并采用动力学对重力变化的形变和运动进行推断[4]。由于已知地球地质、板块运行等地球地质信息，只需要根据地球重力场数据，通过统一归算改正，就可以得到重力改正、重力异常等重力场信息，为此提出重力场深部结构与动力学特征分析这一课题。

1 分析重力场深部结构

此次分析重力场深部结构，将采用反演的方式，分析重力场深部结构，形成重力场深部莫霍面结构图。

1.1 计算重力场场源边界

计算重力场边界，需要计算重力场局部梯度，确定重力场中心和边缘位置，从而确定重力场场源边界。因此，假设重力场数据获取点为x、y、z三点，其重力场场源点则为x0、y0、z0三点，则重力场水平梯度的最大值max{T(x，y)，θ}为

(1)

式中：∂x和∂y为重力场数据获取点x、y两个方向上的导数；θ表示重力场水平梯度方向[5]。此时，根据式(1)，所提取出的重力场水平梯度的最大值，即可确定重力场场源边界及梯度变化情况，为重力场深部结构的莫霍面反演奠定基础。……