[發明專利]一種地震數據子波相位轉換處理的方法

有效

• 申請號: 201010197068.5
• 申請日: 2010-06-02
• 公開（公告）號: CN102269824A
• 公開（公告）日: 2011-12-07
• 發明（設計）人: 高少武;趙海珍
• 申請（專利權）人: 中國石油天然氣集團公司;中國石油集團東方地球物理勘探有限責任公司
• 主分類號: G01V1/32
• 分類號: G01V1/32
• 代理公司: 北京市中實友知識產權代理有限責任公司 11013
• 代理人: 劉天語;石筱巖

1.一種地震數據子波相位轉換處理的方法，其特征是包括以下步驟：

1)用地震人工震源激發和采集地震數據并做預處理；

2)確定地震數據子波參數和計算地震數據零相位子波；

3)對子波進行Z變換；

4)對子波Z變換按下式分解成N個相乘的因式：

B(z)=Πi=1N(1-b^iz-1)=(1-b^1z-1)(1-b^2z-1)Λ(1-b^Nz-1)---(3)]]>

且每個因式的特征因子按照絕對值從小到大依次排列，

|b^1|≤|b^2|≤Λ≤|b^N|---(4)]]>

5)將特征因子分成單位圓內和單位圓外兩部分，按特征因子絕對值的大小，將特征因子分成兩部分：

第一部分是特征因子的絕對值小于或者等于1，共有K(K≥1)個，其中K是Z平面上單位圓內和單位圓上特征因子的總個數；

第二部分是特征因子的絕對值大于1，共有M(M≥1)個，其中M是Z平面上單位圓外特征因子的總個數，且M+K＝N；

6)根據單位圓外的特征因子，按下式構成將零相位子波轉換為最小相位子波的純相位轉換因子H₁(z)：

H1(z)=Πi=1M(1-b^K+i*z-1)(z-1-1b^K+i)---(6)]]>

根據特征方程單位圓外的特征因子，按下式構成將最小相位子波轉換為零相位子波的純相位轉換因子H₂(z)：

H2(z)=Πi=1M(z-1-1b^K+i)(1-1b^K+i*z-1)---(7)]]>

式中：“*”表示取共軛；

7)根據最小相位子波的純相位轉換因子H₁(z)，按下式將零相位子波轉換為最小相位子波：

W_N(z)＝H₁(z)B_N(z)????????????(8)

根據零相位子波的純相位轉換因子H₂(z)，按下式將最小相位子波轉換為零相位子波：

B_N(z)＝H₂(z)W_N(z)????????????(9)

式中：W_N(z)表示最小相位子波的Z變換，B_N(z)表示零相位子波的Z變換；

8)根據純相位轉換因子H₁(z)，按下式將零相位子波地震數據轉換為最小相位子波地震數據：

Y_W(z)＝H₁(z)X_B(z)????????????(10)

式中：Y_W(z)表示最小相位子波地震數據的Z變換，X_B(z)表示零相位子波地震數據的Z變換；

根據零相位子波的純相位轉換因子H₂(z)，按下式將最小相位子波地震數據轉換為零相位子波地震數據：

Y_B(z)＝H₂(z)X_W(z)????????????(11)

式中：Y_B(z)表示零相位子波地震數據的Z變換，X_W(z)表示最小相位子波地震數據的Z變換；

9)根據純相位轉換因子實現地震子波和地震數據的純相位轉換，繪制相位轉換后的地震數據剖面和存儲相位轉換后的地震數據。