[發(fā)明專利]一種自動識別和消除地震勘探工業(yè)電干擾的方法

有效

• 申請?zhí)枺?/td>: 201010205215.9
• 申請日: 2010-06-11
• 公開（公告）號: CN101907726A
• 公開（公告）日: 2010-12-08
• 發(fā)明（設計）人: 高少武;馬玉寧;趙海珍
• 申請（專利權(quán)）人: 中國石油天然氣集團公司;中國石油集團東方地球物理勘探有限責任公司
• 主分類號: G01V1/36
• 分類號: G01V1/36
• 代理公司: 北京市中實友知識產(chǎn)權(quán)代理有限責任公司 11013
• 代理人: 劉天語

1.一種自動識別和消除地震勘探工業(yè)電干擾的方法，其特征是包括以下步驟：

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

2)確定工業(yè)電干擾的初始頻率；

3)采用下式計算工業(yè)電干擾余弦函數(shù)和正弦函數(shù)：

cos2πft=Rxx(t)Rxx(0)---(1)]]>

式中，

f----工業(yè)電干擾頻率；

t----工業(yè)電干擾時間采樣，也是地震數(shù)據(jù)自相關(guān)函數(shù)時間采樣；

R_xx(t)----地震數(shù)據(jù)自相關(guān)函數(shù)；

R_xx(0)----地震數(shù)據(jù)自相關(guān)函數(shù)零延遲值；

4)采用下式計算工業(yè)電干擾正弦函數(shù)：

sin2πft=Pxx(0)Rxx(t)-Rxx(0)Pxx(t)Rxx(0)Rxx2(0)-Pxx2(0)---(2)]]>

式中，

P_xx(t)----地震數(shù)據(jù)褶積函數(shù)；

P_xx(0)----地震數(shù)據(jù)褶積函數(shù)零延遲值；

5)根據(jù)工業(yè)電干擾余弦函數(shù)和正弦函數(shù)構(gòu)成余弦函數(shù)矩陣和正弦函數(shù)矩陣；

6)按照如下公式計算確定工業(yè)電干擾振幅向量：

a＝(CC^T)^-1Cx????????????????(5)

對于工業(yè)電余弦函數(shù)振幅向量，C是工業(yè)電干擾余弦函數(shù)矩陣c，a表示工業(yè)電干擾余弦函數(shù)系數(shù)向量；

對于工業(yè)電正余弦函數(shù)振幅向量，C是由工業(yè)電干擾余弦函數(shù)矩陣c和工業(yè)電干擾余弦函數(shù)矩陣s按照下式構(gòu)成，

C=cs---(6)]]>

a是由工業(yè)電干擾正弦函數(shù)系數(shù)向量A和余弦函數(shù)系數(shù)向量B構(gòu)成，

a=AB;---(7)]]>

7)采用下式計算余弦函數(shù)自適應減和正余弦函數(shù)自適應減確定工業(yè)電干擾：

余弦函數(shù)計算：

y＝C^Ta????????????????????????????(8)

式中，

y----工業(yè)電干擾向量；

C----工業(yè)電干擾余弦函數(shù)矩陣；

a----工業(yè)電干擾余弦函數(shù)系數(shù)向量；

T----向量或者矩陣轉(zhuǎn)置。

且

y^T＝(y₁，y₂，Λ，y_N)

a^T＝(A_-L，A_-L+1，A_-L+2，Λ，A_-1，A₀，A₁，A，A_L-2，A_L-1，A_L)

正余弦函數(shù)計算：

y＝s^TA?+c^TB????????????????????????(9)

式中，

y-----工業(yè)電干擾向量；

s----工業(yè)電干擾正弦函數(shù)矩陣

c----工業(yè)電干擾余弦函數(shù)矩陣

A----工業(yè)電干擾正弦函數(shù)系數(shù)向量；

B----工業(yè)電干擾余弦函數(shù)系數(shù)向量。

y^T＝(y₁，y₂，Λ，y_N)

A^T＝(A_-L，A_-L+1，A_-L+2，Λ，A_-1，A₀，A₁，Λ，A_L-2，A_L-1，A_L)

B^T＝(B_-L，B_-L+1，B_-L+2，Λ，B_-1，B₀，B₁，Λ，B_L-2，B_L-1，B_L)；

8)采用下述公式計算地震有效信號：

S_i＝x_i-y_i????????????????????(10)

式中：原始地震數(shù)據(jù)x_i，由野外數(shù)據(jù)采集得到；估算的工業(yè)電干擾y_i，由工業(yè)電干擾計算公式計算得到；S_i是消除工業(yè)電干擾后的地震有效信號。