[發(fā)明專利]一種地震數(shù)據(jù)快速相位校正處理的方法

有效

• 申請(qǐng)?zhí)枺?/td>: 201210016370.5
• 申請(qǐng)日: 2012-01-18
• 公開（公告）號(hào): CN103217716A
• 公開（公告）日: 2013-07-24
• 發(fā)明（設(shè)計(jì)）人: 高少武;趙波;祝樹云;祝寬海
• 申請(qǐng)（專利權(quán)）人: 中國(guó)石油天然氣集團(tuán)公司;中國(guó)石油集團(tuán)東方地球物理勘探有限責(zé)任公司
• 主分類號(hào): G01V1/36
• 分類號(hào): G01V1/36
• 代理公司: 北京市中實(shí)友知識(shí)產(chǎn)權(quán)代理有限責(zé)任公司 11013
• 代理人: 劉天語(yǔ);任潔

1.一種地震數(shù)據(jù)快速相位校正處理的方法，特點(diǎn)是采用以下步驟：

1)用人工震源激發(fā)和采集地震數(shù)據(jù)并做預(yù)處理；

2)計(jì)算凱澤窗希爾伯特因子序列；

3)計(jì)算希爾伯特地震道序列；

4)確定相位校正因子特征方程的系數(shù)值；

按照以下公式分別計(jì)算相位校正因子特征方程的系數(shù)值A(chǔ)、B、C、D和E：

A＝Pb

B＝2(Pc-2Pa)

C＝3(Pd-Pb)

D＝2(2Pe-Pc)

E＝-Pd???????????(5)

式中，Pa、Pb、Pc、Pd和Pe是五個(gè)常數(shù)值，僅僅與地震數(shù)據(jù)序列x[n]和希爾伯特地震道序列x_H[n]有關(guān)；

按照以下公式計(jì)算Pa、Pb、Pc、Pd和Pe五個(gè)常數(shù)值：

Pa=Σi=1Nx4[i]-3(Σi=1Nx2[i])2]]>

Pb=12Σi=1Nx2[i]Σi=1Nx[i]xH[i]-4Σi=1Nx3[i]xH[i]]]>

Pc=6Σi=1Nx2[i]xH2[i]-6Σi=1Nx2[i]Σi=1NxH2[i]-12(Σi=1Nx[i]xH[i])2]]>

Pd=12Σi=1NxH2[i]Σi=1Nx[i]xH[i]-4Σi=1Nx[i]xH3[i]]]>

Pe=Σi=1NxH4[i]-3(Σi=1NxH2[i])2---(6)]]>

式中，x[i]為地震數(shù)據(jù)序列，x_H[i]為地震數(shù)據(jù)序列對(duì)應(yīng)的希爾伯特地震道序列，由方程(4)計(jì)算確定，i為地震數(shù)據(jù)樣點(diǎn)順序號(hào)，N表示地震數(shù)據(jù)樣點(diǎn)個(gè)數(shù)；

5)計(jì)算相位校正因子特征方程的根；

按照以下公式構(gòu)造相位校正因子特征方程：

A+Bx+Cx²+Dx³+Ex⁴＝0??????(7)

式中，A、B、C、D和E為根據(jù)相位校正因子特征方程的五個(gè)系數(shù)值。x表示相位校正因子特征方程自變量，是相位校正因子的正切函數(shù)值。求解方程(7)，可以得到相位校正因子特征方程的對(duì)應(yīng)四個(gè)根x₁、x₂、x₃和x₄；

按照以下公式計(jì)算相位校正因子φ：

φ＝arctgx，x＝tgφ????????(8)

式中，x為相位校正因子特征方程自變量的根。

根據(jù)相位校正因子特征方程對(duì)應(yīng)的四個(gè)根x₁、x₂、x₃和x₄，由方程(8)可以得到四個(gè)相位校正因子φ₁、φ₂、φ₃和φ₄；

6)計(jì)算相位校正因子峰度函數(shù)值并確定最佳相位校正因子；

按照以下公式計(jì)算峰度函數(shù)值kurt(φ)

kurt(φ)＝E{y⁴[n]}-3(E{y²[n]})²????(9)

式中，E{·}表示期望值運(yùn)算符，y[n]是地震數(shù)據(jù)相位校正序列。已知地震數(shù)據(jù)序列x[n]，相位校正因子φ，由方程(11)計(jì)算地震數(shù)據(jù)相位校正序列。

按照以下公式確定最佳相位校正因子φ_best：

φbest=Maxφbest∈[φ1,φ2,φ3,φ4]{kurt(φ1),kurt(φ2),kurt(φ3),kurt(φ4)}---(10)]]>

式中，kurt(φ₁)、kurt(φ₂)、kurt(φ₃)和kurt(φ₄)分別是相位校正因子特征方程計(jì)算的四個(gè)相位校正因子φ₁、φ₂、φ₃和φ₄對(duì)應(yīng)的峰度函數(shù)值，其中四個(gè)峰度函數(shù)值中最大值對(duì)應(yīng)的相位校正因子就是最佳相位校正因子；

7)實(shí)現(xiàn)相位校正；

y[n]＝x[n]cosφ_best-x_H[n]sinφ_best????(11)

式中，φ_best是最佳相位校正因子，由方程(10)確定，單位是弧度，x[n]是地震數(shù)據(jù)序列，x_H[n]是希爾伯特地震道序列，由方程(4)計(jì)算確定，y[n]是地震數(shù)據(jù)相位校正序列，n是地震數(shù)據(jù)樣點(diǎn)順序號(hào)，n＝1，2，3，...，N，N是地震數(shù)據(jù)樣點(diǎn)個(gè)數(shù)；

由步驟2)----步驟6)確定出最佳相位校正因子φ_best，根據(jù)相位校正因子φ_best，由方程(11)進(jìn)行相位校正處理；

8)實(shí)現(xiàn)串聯(lián)快速相位校正；

重復(fù)步驟2)----步驟7)，就實(shí)現(xiàn)了串聯(lián)快速相位校正因子確定和快速相位校正處理，以實(shí)現(xiàn)消除地震數(shù)據(jù)剩余子波的快速相位校正處理；

9)繪制消除剩余子波后的地震數(shù)據(jù)剖面和存儲(chǔ)消除剩余子波后的地震數(shù)據(jù)。