[發明專利]基于測量的重力異常數據提取局部重力異常的方法

在審

• 申請號: 202111055317.1
• 申請日: 2021-09-09
• 公開（公告）號: CN113761457A
• 公開（公告）日: 2021-12-07
• 發明（設計）人: 孟慶奎;舒晴;高維;周堅鑫;徐光晶;張文志
• 申請（專利權）人: 中國自然資源航空物探遙感中心
• 主分類號: G06F17/10
• 分類號: G06F17/10;G06F17/13;G01V7/00
• 代理公司: 北京云嘉湃富知識產權代理有限公司 11678
• 代理人: 程凌軍

一種基于測量的重力異常數據提取局部重力異常的方法，包括以下步驟：根據重力異常數據Z(x，y)和切割半徑r計算與點(x，y)相距r的某4點重力異常的平均值A(x，y)；計算加權參數e；比較區域重力異常與局部重力異常的曲線彎曲方向；對于彎曲方向相反的情形，通過切割區域重力異常修正量進行修正，確定切割半徑r區域重力異常；計算并提取局部重力異常L(x，y)＝Z(x，y)?R(x，y)。該方法可以改善和提高局部重力異常的準確度。還提供了一種基于該局部重力異常的油氣勘探方法。

技術領域

本發明屬于重力異常數據處理以及油氣勘探的技術領域，涉及一種基于測量的重力異常數據提取局部重力異常的方法。

背景技術

一般獲得的重力數據為布格重力異常數據，它是當代重力研究中獲取最為簡單且研究最為廣泛的重力異常數據。布格重力異常是由區域重力異常和局部重力異常疊加組合而成。其中，由分布范圍較大且影響較深的地質因素(例如：區域地質構造)所引起的重力異常稱為區域重力異常，由地質因素(例如：礦產，油氣)所引起的小范圍的重力異常稱為局部重力異常。由此可見，局部重力異常的提取不僅是資源勘探的重要研究內容，同時也有助于地質構造的研究。在此基礎下，如何準確地提取出局部重力異常也就成為重要問題。

插值切割法提取局部重力異常是現有技術中常用的方法之一，其以當前計算點場值與周圍多點平均值的插值運算為切割算子，通過連續切割，得到切割區域場(區域重力異常)，進而將其從位場異常中減去，便得到切割局部場。根據插值切割算子的不同，依次形成了基于四點圓周插值切割算子、八點窗口插值切割算子和動態改進型插值切割算子的插值切割方法。

然而，現有方法的準確性有待進一步提升。

因此，有必要研究一種基于測量的重力異常數據提取局部重力異常的方法來解決上述的一個或多個技術問題。

發明內容

為解決上述至少一個技術問題，申請人通過研究發現，現有的插值切割法忽視了產生異常的不同地質屬性，得到的計算結果經常會產生一定的誤差。對此，申請人通過從實際地質屬性出發，結合重力位場局部異常(局部重力異常)和區域異常(區域重力異常)的四種疊合關系，對插值切割法加以改進，進而提出一種基于異常約束條件的插值切割法，提高了獲取的局部重力異常的準確度。

根據本發明一方面，提供了一種基于測量的重力異常數據提取局部重力異常的方法，其特征在于包括以下步驟：

通過航空、航天、船載或地面平臺測量重力異常數據Z(x,y)，該重力異常數據Z(x，y)由區域重力異常R(x，y)和局部重力異常L(x，y)組成；

確定切割半徑r；

根據重力異常數據Z(x，y)和切割半徑r計算與點(x，y)相距r的某4點重力異常的平均值A(x，y)；

計算加權參數e；

比較區域重力異常與局部重力異常的曲線彎曲方向；

對于彎曲方向相反的情形，區域重力異常R_n(x，y)確定為：

ΔZ(x，y)為切割區域重力異常修正量，第1次切割的區域重力異常用 R₁(x，y)表示，令R₁(x，y)作為下次迭代的重力異常數據Z(x，y)，通過該迭代的重力異常數據Z(x，y)確定第2次切割的區域重力異常R₂(x，y)；依次迭代下去，最終有

于是切割半徑r區域重力異常

計算并提取局部重力異常L(x，y)＝Z(x，y)-R(x，y)。

根據本發明又一方面，計算加權參數e具體為；

e＝c+d(0≤e≤2)