[发明专利]一种基于积分迭代算法的航空重力向下延拓方法

权利要求书

1.一种基于积分迭代算法的航空重力向下延拓方法，其特征在于，在对测区航空重力测量数据进行向下延拓时，将航空重力测量数据作为积分迭代的初始值进行向上延拓；将向上延拓得到的向上延拓结果与航空重力测量数据求差后作为向下延拓结果的反馈量；将反馈量与初始值叠加得到新的迭代初始值进行下一次迭代。

2.根据权利要求1所述的基于积分迭代算法的航空重力向下延拓方法，其特征在于，所述向下延拓深度不超过观测数据点距的20倍。

3.根据权利要求1或2所述的基于积分迭代算法的航空重力向下延拓方法，其特征在于，具体步骤为：

（1）将δg_r(x,y)作为待求面的初始值，即

（2）根据下式，将向上延拓至航空重力测量数据高度，得到

δgr(x,y)=F-12{K~(u,v)·F2{δgR(x,y)}}]]>

其中，²F与²F^-1分别表示二维傅立叶变换与傅立叶反变换；K~(u,v)=H2π∫-∞∞∫-∞∞e-i(ux+vy)(x2+y2+H2)3/2dxdy=e-Hu2+v2,]]>积分核K(x,y)=12πH(x2+y2+H2)3/2,]]>H为距离地面的高度；

（3）对向上延拓得到的重力异常数据与航空重力测量数据δg_r(x,y)进行比较，并判断是否满足积分迭代过程的结束条件，即：

当时，结束迭代过程，跳至步骤（6）；

当||δgrk(x,y)-δgr(x,y)||2≥ϵ]]>时，继续下一步；

（4）选择合适的迭代步长s，并按照下列积分迭代公式对待求面的初始值进行更新，得到新的初始值

δgRk+1(x,y)=δgRk(x,y)+s[δgrk(x,y)-δgr(x,y)]]]>

其中k是迭代次数，s是迭代步长；表示经过k次迭代得到的向下延拓结果，是向上延拓至相应高度的结果；

（5）令k=k+1，重复步骤（2）～（4）；

（6）延拓结果为：δgR(x,y)=δgRk(x,y)+s[δgrk(x,y)-δgr(x,y)].]]>