[發明專利]基于時頻譜模擬的時變混合相位地震子波提取方法

在審

• 申請號: 201410212287.4
• 申請日: 2014-05-19
• 公開（公告）號: CN103995289A
• 公開（公告）日: 2014-08-20
• 發明（設計）人: 戴永壽;張漫漫;張亞南;王蓉蓉;張鵬
• 申請（專利權）人: 中國石油大學（華東）
• 主分類號: G01V1/28
• 分類號: G01V1/28
• 代理公司: 暫無信息
• 代理人: 暫無信息

1.基于時頻譜模擬的時變混合相位地震子波提取方法的思想，其特征一在于，該方法依次含有以下步驟：

步驟(1)對非平穩地震記錄進行改進的廣義S變換時頻分析：

NGST(τ,f)=∫-∞∞s(t){12πq|f|pexp(-(τ-t)22q2f2p)×exp(-i2πft)}dt,f≠0]]>

其中s(t)為非平穩地震記錄，為改進的廣義S變換的窗函數；

步驟(2)引入時頻濾波器對地震記錄的時頻譜進行濾波處理：

S(τ,f)=∫-∞∞NGST(t,f)HT(t-tr,f)dtr,f≠0]]>

其中H(t,f)為時頻濾波器，T為時頻濾波器的時窗寬度，t_r為反射系數出現的時間點。時頻濾波器表示如下：

Hτ(t,ω)=[1+(t-τ)2ω2v]v+12×[1+2πarctan(λ|(t-τ)ω|p)]]]>

該濾波器是一個三參數(λ,p,v)的時間和頻率的光滑函數；

步驟(3)固定時間τ的取值得到某一時間點處的地震記錄的頻譜S(T,f)，采用譜模擬的方法擬合出該時刻的子波振幅譜W(T,f)，根據反射層出現的位置改變時間τ的取值，得到相應反射點處的子波振幅譜；

步驟(4)利用時變維納濾波法從地震記錄時頻譜中去除時變子波振幅譜的影響，則處理后的地震記錄中殘留著子波相位；

步驟(5)假設子波相位是時不變的，利用高階累積量的雙譜法估計子波的混合相位譜φ(f)；

步驟(6)將上述時變子波的振幅譜和相位譜結合，得到完整的時變混合相位子波，而將步驟(5)得到的混合相位用于步驟(4)中進一步除去子波相位的影響，便得到反射系數序列。

2.基于時頻譜模擬的時變混合相位地震子波提取方法，其特征二在于該方法充分保留了地震記錄中的非平穩成分，且不需要對反射系數序列進行白噪假設，解決了常規分段提取時變子波方法在相鄰層段子波屬性存在較大變化的情況下，子波估計不準的問題，實現更精細的時變地震子波估計；采用高階雙譜法提取子波的混合相位，拓寬了常規子波提取方法中子波零相位或最小相位的假設，提取的子波相位屬性更符合實際；在提高地震資料分辨率方面，該方法不需要研究地震記錄的動態衰減機制，直接對地震資料進行處理，避免了誤差累積；時頻濾波器的引入，大大提高了時間分辨率，消除了時頻分析時由于頻散現象相鄰反射層之間存在的影響，在此基礎上更準確地進行后續的處理工作。