[發(fā)明專利]一種基于方向小波分析的重磁異常處理方法

有效

• 申請?zhí)枺?/td>: 201210449200.6
• 申請日: 2012-11-12
• 公開（公告）號: CN102944905A
• 公開（公告）日: 2013-02-27
• 發(fā)明（設(shè)計）人: 徐亞;郝天珧;黃松;胡衛(wèi)劍;李志偉
• 申請（專利權(quán)）人: 中國科學(xué)院地質(zhì)與地球物理研究所
• 主分類號: G01V11/00
• 分類號: G01V11/00
• 代理公司: 暫無信息
• 代理人: 暫無信息

1.一種基于方向小波分析的重磁異常處理方法，其適用于對地球物理勘探中的重磁數(shù)據(jù)進(jìn)行處理，該方法包括以下步驟：

第一步、原始數(shù)據(jù)獲取：首先利用重力儀或磁力儀在野外采集重磁資料數(shù)據(jù)，按規(guī)則網(wǎng)格間距網(wǎng)格化，形成網(wǎng)格化數(shù)據(jù)文件作為輸入數(shù)據(jù)；

第二步、對數(shù)據(jù)進(jìn)行處理：

(1)重磁異常頻譜F(k)分析，利用二維傅里葉變換計算輸入數(shù)據(jù)頻譜特征；

(2)根據(jù)重磁異常頻譜特征及異常提取目的，確定所需方向小波基及方向小波的參數(shù)；

其中重磁異常頻譜F(k)為波數(shù)k的函數(shù)，將重磁異常頻譜F(k)變換到極坐標(biāo)系可表述為波數(shù)和方位角(k，φ)的函數(shù)，其中tanφ＝k_y/k_x，k_x為x方向的波數(shù)、k_y為y方向的波數(shù)。；在極坐標(biāo)系下按照k對頻譜異常進(jìn)行統(tǒng)計，獲取重磁異常的振幅譜隨波數(shù)的變換規(guī)律，按照能量集中的區(qū)域選取小波基的中心頻率及相對應(yīng)的帶寬范圍；在極坐標(biāo)系下按照角度φ對頻譜異常進(jìn)行統(tǒng)計，獲取重磁異常的振幅譜隨方位角的變換規(guī)律，通過選取能量集中的優(yōu)勢方位角或與優(yōu)勢方位角呈正交的方向，確定方向小波基的方位角參數(shù)；進(jìn)而提取異常主要分量、提取異常優(yōu)勢變化方向分量，并進(jìn)而分析提取異常中隱含的地質(zhì)體走向、斷裂走向的地質(zhì)信息。

如方向小波基為空間域表達(dá)式，則轉(zhuǎn)化為頻率域表達(dá)式進(jìn)行計算；

(3)重磁異常頻譜與方向小波基相乘獲取各階次小波變換信息：

WTf(s,b)=s2π∫F(k)Ψ(sk)ejkbdk]]>

其中，F(xiàn)(k)為輸入數(shù)據(jù)重磁異常頻譜，Ψ(k)為基本小波或母小波函數(shù)(x)的頻譜，k為波數(shù)，s＞0是尺度因子，b表示位移，ψ_s(x)的頻譜為

將相乘后結(jié)果(即：WT_f(s，b))通過逆傅里葉變換，轉(zhuǎn)化為空間域，實(shí)現(xiàn)重磁異常數(shù)據(jù)的方向小波分析；

(4)執(zhí)行步驟(2)-步驟(3)，對多尺度參數(shù)及多個小波基參數(shù)對重磁數(shù)據(jù)進(jìn)行處理；

第三步：對所獲取的所有處理結(jié)果(即不同尺度s和位移b所獲取的小波變換結(jié)果WT_f(s，b))進(jìn)行異常分析及提取，獲取斷裂、構(gòu)造邊界、特殊地質(zhì)體邊界的各類目標(biāo)體信息。

2.依據(jù)如權(quán)利要求1所述的基于方向小波分析的重磁異常處理方法，其中，篩選的重磁資料數(shù)據(jù)在空間域或頻率域具有特定方向性。

3.依據(jù)如權(quán)利要求2所述的基于方向小波分析的重磁異常處理方法，其中，篩選的重磁資料數(shù)據(jù)在空間域或頻率域在某一方向或多個方向?yàn)橹鲗?dǎo)形態(tài)。

4.依據(jù)如權(quán)利要求1所述的基于方向小波分析的重磁異常處理方法，其中，在第二步的步驟(3)中，使用公式(1)實(shí)現(xiàn)重磁異常數(shù)據(jù)的方向小波分析：

WTf(s,b)=s2π∫F(k)Ψ(sk)ejkbdk---(1)]]>

其中，F(xiàn)(k)為重磁異常頻譜，Ψ(k)為基本小波或母小波函數(shù)ψ(x)的頻譜，k為波數(shù)，s＞0是尺度因子，b反映位移，ψ_s(x)的頻譜為