[發明專利]一種雙磁性層磁異常的剝離方法

有效

• 申請號: 201210212070.4
• 申請日: 2012-06-21
• 公開（公告）號: CN103513284A
• 公開（公告）日: 2014-01-15
• 發明（設計）人: 劉云祥
• 申請（專利權）人: 中國石油天然氣集團公司;中國石油集團東方地球物理勘探有限責任公司
• 主分類號: G01V3/38
• 分類號: G01V3/38
• 代理公司: 北京市中實友知識產權代理有限責任公司 11013
• 代理人: 劉天語;周慧麟

技術領域

本發明屬于地球物理勘探的磁力勘探數據處理領域，具體是磁力勘探資料處理中剝離雙層火成巖磁性層的一種雙磁性層磁異常的剝離方法。

背景技術

磁異常分離是地球物理勘探中磁法（或稱磁力）勘探的數據處理的重要內容；濾波法分離磁異常是常用的方法之一，它的目前處理步驟是：1.磁異常網格化，2.選取濾波方法（如向上延拓、低通濾波、滑動平均濾波等），3.按第二步選定的方法對磁異常網格數據進行濾波，獲得區域磁異常，4.磁異常減去區域磁異常，獲得剩余磁異常，即將磁異常分離為區域磁異常和剩余磁異常。

火山巖是油氣勘探中的一類重要的儲層類型，上述使用的濾波法分離磁異常的方法是定性分離方法，不能準確分離含有上下兩層火成巖的磁異常，且經常與實際資料不符或存在矛盾，效果較差。特別是對于雙層磁性層的磁異常剝離或分離，國內外尚無很好的辦法，尤其對于沉積盆地內的雙層火成巖的磁異常，更是磁力勘探中的難題。地質出版社2005.8，P190-196公開了《地磁場與磁力勘探》的一種方法，均是基于濾波的定性處理，這些方法存在兩個問題，一是分離出的深、淺層磁異常經常與已有地質資料存在不符或矛盾的問題，二是實際處理時存在不同人處理出的結果差異大的問題。因而，現實需要研發一種能較好分離（剝離）含有上下兩層火成巖磁異常的磁異常分離方法。

發明內容

本發明的目的是提供一種利用磁力勘探資料能分離出上、下兩層火成巖的磁異常的雙磁性層磁異常的剝離方法。

本發明通過以下步驟實現：

1）實測工區總場磁異常ΔT_總數據將數據進行網格化處理，再進行化極處理；

2）用工區已知地震勘探得到的構造圖，構置上磁性層頂面深度圖和底界深度圖，得到上磁性層頂面埋深平面網格數據和上磁性層底界埋深平面網格數據；

步驟2）所述的構置是對地震構造圖進行數字化，將其數據加上地面海拔，變為埋深數據。

3）用工區已知鉆井得到的上層火成巖巖芯磁化率值k，構置上磁性層磁化率平面分布，形成上磁性層磁化率網格數據；

步驟3）所述的火成巖巖芯磁化率值，在上磁性層不存在的區域，上磁性層的磁化率網格值取0。

步驟3）所述的上磁性層磁化率平面分布的構置方法是采用對工區已知鉆井得到的上層火成巖巖芯磁化率值k進行網格化。

4）按照以下公式計算節點Ｐ的上磁性層的磁異常ΔT_上：

磁性體V(ξ,η,ζ)在P(x,y,z)點產生的磁異常ΔT_上為：

其中：

磁化強度j＝k.T₀+j_r＝(1+Q).k.T₀；Q為剩感比，k為磁化率，T₀為地磁場強度，μ₀是真空磁常數，取1.0，j是磁化強度，j_r是剩磁強度；R是V(ξ，η，ζ)到P(x,y,z)點的距離，P(x,y,z)是計算點及坐標（x,y,z）位置，V(ξ，η，ζ)是磁性體體元及坐標(ξ，η，ζ)；體積分中V的空間區域為上磁性層頂面和底面包圍的空間；

步驟4）所述的Q剩感比根據測區物性資料而定，為0.1；磁化率k是步驟3）中的網格數據各節點上的磁化率值。

5）利用步驟2）得到的上磁性層頂面網格數據和上磁性層底面網格數據及步驟3）得到的磁化率k的網格數據，按步驟4）的公式計算得到全工區所有網格節點的上磁性層磁性體的磁異常ΔT_上的網格數據。

6）用步驟1)得到的實測工區總場化極磁異常ΔT_總網格數據減去步驟5）得到的上磁性層的磁異常ΔT_上，其差值即下磁性層的磁異常ΔT_下網格數據，ΔT_下＝ΔT_總－ΔT_上，從實測總場化極磁異常ΔT_總網格數據中剝離出上、下磁性層的磁異常ΔT_上、ΔT_下。

本發明對實測數據處理實驗，獲得了上、下兩套火成巖的磁異常，實現了雙層磁異常的定量剝離，滿足了為研究上、下兩層火成巖提供對應的磁異常的需要，處理結果是唯一，且與鉆井吻合，實際資料表明效果明顯。

附圖說明

圖1為工區實測總場化極磁異常ΔT_總圖；

圖2為工區上磁性層頂面埋深圖；

圖3為工區上磁性層底界埋深圖；