[發明專利]一種基于積分迭代算法的航空重力向下延拓方法

有效

• 申請號: 201310322111.X
• 申請日: 2013-07-29
• 公開（公告）號: CN103399350A
• 公開（公告）日: 2013-11-20
• 發明（設計）人: 吳美平;趙磊;張開東
• 申請（專利權）人: 中國人民解放軍國防科學技術大學
• 主分類號: G01V7/00
• 分類號: G01V7/00
• 代理公司: 湖南兆弘專利事務所 43008
• 代理人: 周長清

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)].]]>