[发明专利]一种自动识别和消除地震勘探工业电干扰的方法

摘要

本发明涉及油田的勘探、开发技术中自动识别和消除地震勘探工业电干扰的方法，确定工业电干扰的初始频率，计算工业电干扰余弦函数和正弦函数构成余弦函数矩阵和正弦函数矩阵，计算确定工业电干扰振幅向量，计算余弦函数自适应减和正余弦函数自适应减确定工业电干扰，得到地震有效信号。本发明不需要确定工业电干扰频率，有效提高计算效率。本发明既可以消除地震数据中由高压输电线产生的工业电干扰，也可以消除地震数据中由周期性震动(如发电机等)产生的干扰。

主权项

一种自动识别和消除地震勘探工业电干扰的方法，其特征是包括以下步骤：1)用地震震源激发和采集地震数据并做预处理；2)确定工业电干扰的初始频率；3)采用下式计算工业电干扰余弦函数和正弦函数： cos 2 πft = R xx ( t ) R xx ( 0 ) - - - ( 1 ) 式中，f 工业电干扰频率；t 工业电干扰时间采样，也是地震数据自相关函数时间采样；Rxx(t) 地震数据自相关函数；Rxx(0) 地震数据自相关函数零延迟值；4)采用下式计算工业电干扰正弦函数： sin 2 πft = P xx ( 0 ) R xx ( t ) - R xx ( 0 ) P xx ( t ) R xx ( 0 ) R xx 2 ( 0 ) - P xx 2 ( 0 ) - - - ( 2 ) 式中，Pxx(t) 地震数据褶积函数；Pxx(0) 地震数据褶积函数零延迟值；5)根据工业电干扰余弦函数和正弦函数构成余弦函数矩阵和正弦函数矩阵；6)按照如下公式计算确定工业电干扰振幅向量：a＝(CCT) 1Cx (5)对于工业电余弦函数振幅向量，C是工业电干扰余弦函数矩阵c，a表示工业电干扰余弦函数系数向量；对于工业电正余弦函数振幅向量，C是由工业电干扰余弦函数矩阵c和工业电干扰余弦函数矩阵s按照下式构成， C = c s - - - ( 6 ) a是由工业电干扰正弦函数系数向量A和余弦函数系数向量B构成， a = A B ; - - - ( 7 ) 7)采用下式计算余弦函数自适应减和正余弦函数自适应减确定工业电干扰：余弦函数计算：y＝CTa (8)式中，y 工业电干扰向量；C 工业电干扰余弦函数矩阵；a 工业电干扰余弦函数系数向量；T 向量或者矩阵转置。且yT＝(y1，y2，Λ，yN)aT＝(A L，A L+1，A L+2，Λ，A 1，A0，A1，A，AL 2，AL 1，AL)正余弦函数计算：y＝sTA +cTB (9)式中，y 工业电干扰向量；s 工业电干扰正弦函数矩阵c 工业电干扰余弦函数矩阵A 工业电干扰正弦函数系数向量；B 工业电干扰余弦函数系数向量。yT＝(y1，y2，Λ，yN)AT＝(A L，A L+1，A L+2，Λ，A 1，A0，A1，Λ，AL 2，AL 1，AL)BT＝(B L，B L+1，B L+2，Λ，B 1，B0，B1，Λ，BL 2，BL 1，BL)；8)采用下述公式计算地震有效信号：Si＝xi yi (10)式中：原始地震数据xi，由野外数据采集得到；估算的工业电干扰yi，由工业电干扰计算公式计算得到；Si是消除工业电干扰后的地震有效信号。