重力异常的数据处理.ppt

从位场理论可知，一个未知的异常体在观测平面上所引起的重力异常若已知时，则可将这个观测面展布成一个无限大而面密度不均匀的等效物质面，使这个面上各点的面密度?(ξ，?，0)满足下式这时,在其外部空间任意点引起的重力异常及其各阶导数都将与原来场源在该点产生的异常各阶导数是等效的。由引力位的定义可知，一个密度分布不均匀的无限大物质面，在其上部空间任意点A的引力位为将(9-18)代入式(9-19)可得将式(9-20)对z求偏导数,则可得该点的重力异常表达式同理，还可计算出该点所在平面上的Vxzvzz及vzzz及其它各阶导数。若令式(9-21)式中z=0时，便又可计算出原来重力异常所在观测面上的Vxz，Vzz和Vzzz了。为了便于应用，可将式(9-21)改用柱坐标来计算，并把计算点选在坐标原点的正上方、高度为h的P点，见图9-24,图9-24位场转换计算时的坐标选择对于二度体而言，式（9-21）可变为这时式(9-21)则变为二、重力异常的导数换算重力异常的导数在不同形状地质体上有不同的特征，有助于对异常的解释和分类;重力异常的导数可以突出浅而小的地质体的异常特征而压制区域性深部地质因素的影响，在一定程度上可以划分不同深度和大小异常源产生的叠加异常，且导数的次数越高，这种分辨能力就越强(见图补);重力高阶导数可以将几个互相靠近、埋深相差不大的相邻地质因素引起的叠加异划分开来。这些功能主要是因为导数阶次越高，则异常随中心埋深加大而衰减越快，从水平方向看，基于同样道理，阶次越高的异常范围越小，因而无论从垂向看或从水平方向看，高阶导数异常的分辨能力都提高了。（一）vxz的计算01.由Vxz的物理意义可知，实际是vz在其x方向的变化率，其表达式可为02.式中(Vxz)0表示两点相距为2?x时，中间点0的平均变化率，这只能近似代表水平梯度.01若想提高计算精度,可用最小二乘法对实测值进行函数拟合,然后对x求导数即可.下面具体介绍其原理。02010203设坐标原点定在计算点上，并令对上式求x的导数返回返回根据（9-24），在x=0点的?g对x的导数是vxz(0)=a1，当(1)m=2时，n=5时计算的公式如下（9-26）上述两式中的?x为取数点距。当m=2,n=7时计算的公式如下（9-27）

（二）vzz（gz）的换算

式中?个（0，0，0）为计算点上的重力异常。Ki为各环系数。具体计算时，借助同心圆取数量板，若量板的半径为ri和ri+1，取数半径为并以km为单位，则各环相应系数见表9-2.010203将表中的ki依次代入(9-33)中,得具体计算vzz计算公式（9-34）。01对于二度异常，可用类似方法求出其近似公式为(9-35)02vzzz的换算已知在场源外部,引力位是空间坐标的调和函数,满足拉普拉斯方程对于由于所以在场源外部空间有01其中代入上式解得02（9-36）至今导出的计算公式很多,基本原理相似,常用公式为艾勒金斯第II公式、第I公式和第III公式,介绍如下:01若用符号表示以坐标原点0为圆心,R为半径的一个圆周上重力异常的平均值,则02式中g(R,?)为圆周上某一点的重力值,由于它是坐标位置的调和函数,因此,当R不大时,g(R,?)可以写成对坐标原点的台劳展开式由于x=Rcos?,y=Rsin?,代入(9-37)式,直接积分后得到01其中R的奇次项积分后代入上下限均已消去.a0a1a2的表达式为02这样计算原点0处的重力垂向二次导数的问题,就变成了确定(9-39)中的系数a1了,求得a1再乘上（-4）,就得到计算点的gzz，即1（9-40）2当采用不同方法确定系数a1时,就可以得到不同的计算公式.1哈克公式2在公式(9-38)中略去高次项则有3即由此导出哈克公式为艾勒金斯公式同理,当略去(9-38)中的高次项,得并分别取半径为R,的圆周上的值代入上式,得到图9-8铅垂台阶异常与单斜区域异常的叠加区域异常和局部异常区域异常是叠加异常中的一部分，主要是由分布较广的中、深部地质因素所引起的重力异常。这种异常特征是异常幅值较大，异常范围也较大，但异常梯度小。局部异常是叠加异常中的一部分，主要是指相对区域因素而言范围有限的研究对象引起的范围和幅度较小的异常，但异常梯度相对较大。由于局部异常是布格异常中去掉区域异常后的剩余部分，局部异常也称为剩余异常。区域异常和局部异常是相对而言的，绝对的划分标准，应视研究的问题而言。由图9-9可知，相对异常A而言，异常B都可以看成区域异常;而