[發明專利]一種空地井立體重磁數據協同目標位置反演方法

有效

• 申請號: 201810185427.1
• 申請日: 2018-03-07
• 公開（公告）號: CN108508479B
• 公開（公告）日: 2020-09-11
• 發明（設計）人: 馬國慶;李麗麗;杜曉娟
• 申請（專利權）人: 吉林大學
• 主分類號: G01V1/28
• 分類號: G01V1/28;G01V1/40
• 代理公司: 暫無信息
• 代理人: 暫無信息

1.一種空地井立體重磁數據協同目標位置反演方法，其特征在于，包括如下步驟：

S1、首先在需要勘測地點的地平面計算地面重磁數據，并利用泰勒展開式方程來計算數據精度ε₂；

S2、然后計算航空重磁數據和下井測試井中重磁數據，并分別利用泰勒展開式方程來計算數據精度ε₁和ε₃；

S3、再對單一空間重磁數據進行地質體位置反演，并且根據不同的高度和深度重復步驟S2和地質體位置反演；

S4、根據不同數據的精度分配不同數據權重W，并根據多空間重磁數據按照相應公式進行反演；

S5、最后將單層反演窗口縮小而利用其它空間數據形成補充，重復上述操作，根據不同數據的精度和根據精度設定相應的權重系數，進而組合為矩陣進行反演，最后得出地質體空間位置信息。

2.根據權利要求1所述的一種空地井立體重磁數據協同目標位置反演方法，其特征在于：所述S1和S2中重磁數據的泰勒展開式方程為：

其中，T₀(x′，y′，0)表示觀測高度為0的平面上的異常；其中x′，y′分別表示采樣點的水平坐標；T(x，y，h)表示觀測高度為h的平面上的異常；x，y分別表示計算點的水平坐標；當h＞0時表示向上延拓，當h＜0時表示向下延拓，對公式進行Fourier變換可得到

其中，表示異常T₀(x′，y′，0)的Fourier譜；k_x，k_y分別表示x和y方向上的波數；稱為延拓因子。

3.根據權利要求1所述的一種空地井立體重磁數據協同目標位置反演方法，其特征在于：所述S3中對單一空間重磁數據進行地質體位置反演的原理是重磁數據及其導數滿足歐拉齊次方程：

(x-x₀)·f_x+(y-y₀)·f_y+(z-z₀)·f_z＝-N(f(x-x₀，y-y₀，z-z₀)+A)

其中，f為原始重磁異常；f_x，f_y，f_z分別為異常f在不同方向的導數；x，y，z為觀測點坐標；A為背景異常；N為構造指數；x₀，y₀，z₀為需反演的地質體位置坐標。

4.根據權利要求1所述的一種空地井立體重磁數據協同目標位置反演方法，其特征在于：所述S4中根據多空間重磁數據進行反演時的公式為：

其中，f為原始重磁異常；f_x，f_y，f_z分別為異常f在不同方向的導數；x，y，z為觀測點坐標；N為構造指數；x₀，y₀，z₀為需反演的地質體位置坐標；h和-h均代表觀測面高度。

5.根據權利要求1所述的一種空地井立體重磁數據協同目標位置反演方法，其特征在于：所述S5中根據不同數據的精度和根據精度設定相應的權重系數，進而組合為矩陣進行反演的公式為：

其中，W₁、W₂、W₃分別代表航空、地面、井中三個觀測高度異常反演時的權重系數；f為原始重磁異常；f_x，f_y，f_z分別為異常f在不同方向的導數；x，y，z為觀測點坐標；N為構造指數；x₀，y₀，z₀為需反演的地質體位置坐標；h和-h均代表觀測面高度。