|
|
|
 |
|
|
| 一種地球物理勘探中的重磁數據處理方法，其包括步驟：在波數域將重磁資料頻譜分別與任意方向的ｎ階導數或積分頻率響應函數、延拓回返響應函數及差分法下延算子頻率響應函數相乘，得到下延后的頻譜；再經過快速傅氏反變換，獲得向下延拓后的三維重磁數據體。利用向下延拓之后的三維重磁數據體結合地震、地質、鉆井資料進行對比、分析并提交下一步勘探的有利區帶及目標。該方法不僅有更大的向下延拓深度并且有較好的抗噪能力，從而獲得穩定的下半空間三維重磁資料分布圖像。 |
|
|
|
|
|
|
|
|
一種地球物理勘探中的重磁數據處理方法
一種地球物理勘探中的重磁數據處理方法，該方法包括以下步驟： �。ǎ保┦紫壤弥亓x或磁力儀在野外采集重磁資料，或者在室內對原有的重磁異常圖進行數字化而得到重磁資料，然后對野外采集或室內數字化的重磁資料進行各種常規校正及改正； （２）利用Ｓｕｒｆｅｒ中的Ｋｒｉｇｉｎｇ插值方法將重磁資料網格化，形成規則網重磁資料ｆ↓［ａ］（ｘ，ｙ），其中ｘ、ｙ分別為重磁資料的平面網格點坐標； （３）進行鑲邊處理：將實際重磁資料ｆ↓［ａ］（ｘ，ｙ）的線數和點數擴大為２的整數冪的重磁資料ｆ↓［ｂ］（ｘ，ｙ）；利用低精度的重磁資料進行補充鑲邊； �。ǎ矗├每焖俑凳险儞Q（ＦＦＴ）將ｆ↓［ｂ］（ｘ，ｙ）變換到波數域，獲得實測重磁資料的頻譜Ｆ↓［ｂ］（ｕ，ｖ），Ｆ↓［ｂ］（ｕ，ｖ）＝ＦＦＴ［ｆ↓［ｂ］（ｘ，ｙ）］，其中ｕ、ｖ分別為ｘ和ｙ方向的波數； �。ǎ担┰诓〝涤颍瑢㈣傔吅蟮闹卮刨Y料頻譜Ｆ↓［ｂ］（ｕ，ｖ）與重磁場任意方向的ｎ階導數或積分頻率響應函數ｑ↑［ｎ］（ｕ，ｖ）＝［２πｉ（αｕ＋βｖ）＋γ＊］↑［ｎ］相乘，增強局部或區域重磁異常，即： �。啤郏悖荩ǎ�，ｖ）＝Ｆ↓［ｂ］（ｕ，ｖ）×ｑ↑［ｎ］（ｕ，ｖ） 其中ｎ為有理數，ｉ為虛數；α、β、γ分別為所求導數或積分方向的三個方向余弦；ｎ為導數或積分的階數，即分辨率參數； （６）在波數域，根據向下延拓深度ｈ，將增強后的局部或區域重磁異常Ｆ↓［ｃ］（ｕ，ｖ）與延拓回返頻域響應函數 φ（ｕ，ｖ）＝［ｅｘｐ（－２πｈ＊）］↑［ｐ］（６－ｅｘｐ（－２πｈ＊－２ｃｏｓ（２πｈｕ）－２ｃｏｓ（２πｈｖ）］↑［ｐ］相乘，抑制高頻成分，使其滿足向下延拓的條件，即： �。啤郏洌荩ǎ�，ｖ）＝Ｆ↓［ｃ］（ｕ，ｖ）×φ（ｕ，ｖ） 其中ｈ為下延深度，ｐ為延拓回返次數；ｈ，ｐ為正有理數； �。ǎ罚┰诓〝涤颍谝种聘哳l成分的基礎上，與差分法向下延拓頻域響應函數ｄ（ｕ，ｖ）＝（６－ｅｘｐ（－２πｈ＊－２ｃｏｓ（２πｈｕ）－２ｃｏｓ（２πｈｖ）））相乘，得到下延后的頻譜Ｆ↓［ｅ］（ｕ，ｖ） �。啤郏澹荩ǎ�，ｖ）＝Ｆ↓［ｄ］（ｕ，ｖ）×ｄ（ｕ，ｖ）； �。ǎ福┰俳涍^快速傅氏反變換（ＩＦＦＴ），獲得向下延拓后的重磁資料ｆ↓［ｅ］（ｘ，ｙ），ｆ↓［ｅ］（ｘ，ｙ）＝ＩＦＦＴ［Ｆ↓［ｅ］（ｕ，ｖ）］； �。ǎ梗┻M行反鑲邊處理，即從ｆ↓［ｅ］（ｘ，ｙ）中提取原數據大小的向下延拓深度ｈ后的重磁數據ｆ↓［ｆ］（ｘ，ｙ）； �。ǎ保埃┲貜筒襟E（６），從而獲得下半空間三維重磁數據體ｆ↓［ｇ］（ｘ，ｙ，ｚ），其中，ｘ，ｙ，ｚ為重磁異常的空間網格點坐標； �。ǎ保保├孟蛳卵油刂蟮娜S重磁數據體結合地震、地質、鉆井資料進行對比、分析并提交下一步勘探的有利區帶及目標。
|
|
|
|
|
|
|
| 專利號: |
200810110500 |
| 申請日: |
2008年6月13日 |
| 公開/公告日: |
2008年10月15日 |
| 授權公告日: |
|
| 申請人/專利權人: |
楊輝 |
| 國家/省市: |
北京(11) |
| 郵編: |
100083 |
| 發明/設計人: |
楊輝 |
| 代理人: |
|
| 專利代理機構: |
() |
| 專利代理機構地址: |
() |
| 專利類型: |
發明 |
| 公開號: |
101285896 |
| 公告日: |
|
| 授權日: |
|
| 公告號: |
000000000 |
| 優先權: |
|
| 審批歷史: |
|
| 附圖數: |
3 |
| 頁數: |
11 |
| 權利要求項數: |
3 |
| |
| |
|
|
