[發明專利]基于廣義瞬時相位及P范數負模的地震弱信號處理方法

無效

• 申請號: 201310692354.2
• 申請日: 2013-12-17
• 公開（公告）號: CN103645504A
• 公開（公告）日: 2014-03-19
• 發明（設計）人: 李緒宣;曹思遠;胡光義;范廷恩;董建華;王宗俊;井涌泉;張顯文;樊鵬軍;高云峰
• 申請（專利權）人: 中國海洋石油總公司;中海油研究總院;中國石油大學（北京）
• 主分類號: G01V1/28
• 分類號: G01V1/28;G01V1/36
• 代理公司: 北京紀凱知識產權代理有限公司 11245
• 代理人: 徐寧;關暢

1.一種基于廣義瞬時相位及P范數負模的地震弱信號處理方法，其包括以下內容：

1）對地震信號進行識別，具體過程為：

1.1）對輸入的地震信號x(n)進行P階廣義希爾伯特變換并得到廣義瞬時相位φ(t)；

1.2）求取廣義瞬時相位P范數負模γ(t)；

1.3）根據廣義瞬時相位P范數負模γ(t)提取地震弱信號的相關屬性參數，通過提取的屬性參數對地震有效信號進行識別；

2）對地震信號的噪聲進行去除，得到信噪比比較高的地震信號，具體過程為：

2.1）對輸入的地震信號x(n)進行P階廣義希爾伯特變換并得到廣義瞬時相位φ(t)；

2.2）對計算得到的廣義瞬時相位φ(t)進行平滑濾波處理，得到濾波后的廣義瞬時相位信號

2.3）采用逐步迭代的方法求取地震信號的瞬時振幅a(n)；

2.4）根據求取的廣義瞬時相位信號和瞬時振幅a(n)重構地震信號

2.5）將重構地震信號與測井資料得到的數據進行比對，對重構地震信號的合理性進行判斷，如果合理則輸出重構信號如果不合理，則返回步驟2.2）。

2.如權利要求1所述的基于廣義瞬時相位及P范數負模的地震弱信號處理方法，其特征在于：廣義瞬時相位式中，hr(t)和hi(t)是復數道h(t)的實部和虛部。

3.如權利要求1所述的基于廣義瞬時相位及P范數負模的地震弱信號處理方法，其特征在于：廣義瞬時相位P范數負模γ(t)：

γ(t)=-|arctanhi(t)hr(t)|]]>

式中，hr(t)和hi(t)是復數道h(t)的實部和虛部。

4.如權利要求2所述的基于廣義瞬時相位及P范數負模的地震弱信號處理方法，其特征在于：廣義瞬時相位P范數負模γ(t)：

γ(t)=-|arctanhi(t)hr(t)|]]>

式中，hr(t)和hi(t)是復數道h(t)的實部和虛部。

5.如權利要求1或2或3或4所述的基于廣義瞬時相位及P范數負模的地震弱信號處理方法，其特征在于：所述步驟2.3）采用逐步迭代的方法求取地震信號的瞬時振幅a(n)，包括以下步驟：

①將地震信號表示為廣義瞬時相位信號的形式，即：

②采用一假設的正數c代替①中振幅的初始值迭代值|a₀(n)|，令

③計算的N點DFT得到A₁(k)，由于解析信號的傅立葉變換是的單邊函數，在處取值為0，所以有A1′(k)=A1(k)0≤N2-10N2≤k≤N-1,]]>式中N都取2^M形式；

④求A′₁(k)的傅立葉逆變換IDFT得到其中|a′₁(n)|是a(n)的新估算值；

⑤判斷|a′₁(n)|是否滿足精度要求，如果不滿足精度要求，則返回步驟③繼續迭代；假設迭代次數為p，則迭代時采用代替迭代結果如果滿足精度要求，則輸出得到結果a(n)=|a'_p(n)|。