楊 爽， 王彥春， 何俊強

(中國地質(zhì)大學(北京) 地球物理與信息技術學院，北京 100083)

0 引言

在海洋地震數(shù)據(jù)采集過程中，過往船只可以看成是一種震源，當其與采集拖纜達到一定的距離時，可在地震記錄上形成干擾波[1]。干擾波的形態(tài)取決于過往船只與采集拖纜的相對位置，當干擾船在采集拖纜的延長線上時，干擾波為直線型，具有直達波的性質(zhì)；當干擾船在采集拖纜的垂直線上時，干擾波為雙曲線型，具有反射波的性質(zhì)[2]。

直線型船干擾可用高精度拉東變換等線性干擾噪聲壓制方法進行壓制[3]，而雙曲型干擾不同于面波、散射等規(guī)則干擾波，且其頻帶與有效波重疊，用濾波等方法很難壓制。目前壓制雙曲型干擾主要有兩種方法：①將雙曲型干擾在局部時窗內(nèi)當成直線型干擾[4]，并應用噪線性的方法進行壓制，該方法對雙曲型干擾的兩翼有較好的壓制，而對于反射噪音中能量最強的雙曲線頂部卻無能為力；②根據(jù)雙曲型時距曲線和正常時差來得到干擾源的位置，再以實際地震采集觀測系統(tǒng)模擬來自干擾源的干擾波，以此校正炮集記錄，拉平雙曲型干擾[5-6]，該方法對于海上平臺等位置固定的干擾源產(chǎn)生的雙曲型干擾波有很好地壓制效果，但是海上船的位置多變，如果無法對其運動特征進行有效描述，則需要對每一炮地震數(shù)據(jù)進行干擾船位置的掃描預測，這將大大增加計算量。

由于海上船相對于地震采集拖纜的位置一般都是隨機變化的，其產(chǎn)生噪聲干擾在共偏移距道集內(nèi)為離散的不規(guī)則強能量點[7-8]，而異常振幅衰減法[9-13]對于道集中的強能量點有很好的壓制效果，筆者本文提出將炮集數(shù)據(jù)轉(zhuǎn)換至共偏移距道集，再應用異常振幅衰減法分時、頻地壓制海上船噪聲干擾的方法，并對南美某工區(qū)的海上地震資料進行了試驗分析。

1 方法原理

1.1 異常振幅衰減法

異常振幅噪聲和有效信號的主要差別在于振幅和頻率。在某一給定的時窗內(nèi)，如果地震數(shù)據(jù)在某一頻段的平均振幅超出其計算門檻值，則認為該頻段地震數(shù)據(jù)含有噪聲，即有

Eftk≥Tftk

(1)

其中：Eftk、Tftk分別為第k道、時窗t、頻段f的平均振幅和計算門檻值。某一頻段內(nèi)地震數(shù)據(jù)的平均振幅為

(2)

其中：Ai為第k道、時窗t、頻段f的第i個樣點振幅；Nftk為第k道、時窗t、頻段f的樣點數(shù)。相應頻段的地震數(shù)據(jù)的計算門檻值為

Tftk=Mftk·Ct

(3)

其中：Mftk為第k道、時窗t、頻段f的中值振幅；Ct為給定的、隨時間變化的門檻值系數(shù)，用于限定有效信號與噪聲的幅值差異，由于地震信號存在縱向衰減，為了能準確識別噪聲，門檻值系數(shù)在縱向上逐漸減小。中值振幅Mftk由給定時窗內(nèi)所有道的時間方向平均振幅Eftk做中值濾波求取，即

Mftk=Lk*Eftk

(4)

其中：Lk表示中值濾波算子，算子長度為給定時窗內(nèi)的地震道數(shù)。

在識別出異常振幅噪聲之后，給定衰減因子，對含噪聲頻帶的振幅進行衰減，可將噪聲頻率成分充分壓制。

1.2 海上船噪聲干擾壓制

船噪聲在炮集內(nèi)為規(guī)則干擾，在共偏移距道集中為離散的不規(guī)則干擾，因此可將炮集轉(zhuǎn)換成共偏移距道集，然后用異常振幅衰減法進行壓制，其流程如圖1所示。

圖1 海上船噪聲干擾壓制流程Fig.1 The procedures of offshore vessel noise interference suppression

2 實際資料分析

圖2(a)為南美海上某工區(qū)的單炮記錄，從圖2(a)中可以看到，3.5 s～6.0 s船噪聲干擾為多條低速的雙曲線，嚴重覆蓋了深層的有效信號，如果不加以壓制，勢必會影響疊加剖面的效果，為此，對該地區(qū)的船噪聲干擾進行異常振幅衰減壓制。將炮集記錄轉(zhuǎn)換到共偏移距域后，分別設置低頻段和高頻段的時窗、壓制門檻值、濾波道數(shù)、衰減比例的參數(shù)如表1、表2所示。

表1 低頻段異常振幅衰減參數(shù)設置

表2 高頻段異常振幅衰減參數(shù)設置

圖2分別顯示了船噪聲干擾壓制前、后的炮集記錄以及共炮域內(nèi)被壓制的船噪聲干擾，圖3分別顯示了共偏移距域剖面、壓制船噪聲干擾后共偏移距域剖面和共偏移距域內(nèi)船噪聲干擾。從圖2、圖3中可以看出，船噪聲干擾在共偏移距域內(nèi)為異常振幅點，利用異常振幅衰減技術有效地壓制了異常振幅，在轉(zhuǎn)換回的共炮域中，船噪聲干擾基本被壓制，信噪比明顯提高。

圖4是未壓制船噪聲干擾的疊加剖面，在全區(qū)存在斜干擾，在左下角3.5 s～5.0 s區(qū)域還存在弧形干擾，使用上述方法處理后得到圖5，對比剖面可以明顯地看出，船噪聲干擾已經(jīng)被壓制，尤其是深部的位置，基本恢復了反射同相軸的原本面貌。

圖2 共炮域海上船噪聲干擾壓制Fig.2 Offshore vessel noise interference suppression in common shot domain(a)原始炮集；(b)船噪聲干擾壓制后炮集；(c)船噪聲干擾

圖3 共偏移距域海上船噪聲干擾壓制Fig.3 Offshore vessel noise interference suppression in common offset domain(a) 原始共偏移距域剖面；(b)船噪聲干擾壓制后剖面；(c)船噪聲干擾

3 結論

海洋地震資料中通常包含船噪聲干擾，該噪聲多以雙曲形態(tài)出現(xiàn)，常規(guī)濾波法難以將其壓制，由于船的移動性，干擾源預測重構減去技術也很難運用。筆者基于移動船干擾在共偏移距域內(nèi)呈不規(guī)則異常振幅點的特性，將地震炮集轉(zhuǎn)換至共偏移距道集，然后采用異常振幅衰減法分時、分頻地壓制船噪聲干擾，該方法簡單快速，結合實際的海洋地震資料也取得了良好的應用效果，為海洋地震資料的外源干擾波及側(cè)面繞射波的壓制提供了一定的參考。

參考文獻：

[1] 楊振武. 海洋石油地震勘探—資料采集與處理[M]. 北京:石油工業(yè)出版社，2012.

YANG Z W. Offshore oil seismic exploration-data acquisition and processing[M]. Beijing: Petroleum Industry Press, 2012.(In Chinese)

圖4 原始疊加剖面Fig.4 Original stack profile

圖5 海上船噪聲干擾壓制后的疊加剖面Fig.5 Stack profile after suppressing offshore vessel noise interference

[2] 張雅勤，張美華，寧俊瑞. 海上船噪聲干擾的壓制技術[J]. 物探與化探, 2001,25(4):290-293.

ZHANG Y Q, ZHANG M H, NING J R. The technique for suppressing ship noise interference[J]. Geophysical & Geochemical Exploration,2001,25(4):290-293.(In Chinese)

[3] 王興宇. 海洋地震資料處理外源干擾波壓制方法研究[D].北京：中國地質(zhì)大學，2014.

WANG X Y. The study of exogenous interference suppression method on marine seismic data processing[D]. Beijing: China University of Geosciences, 2014.(In Chinese)

[4] 童思友，徐磊磊，崔樹果,等. 南黃海地震資料疊前有源干擾的自動識別與壓制[J]. 海洋地質(zhì)動態(tài), 2008,24(8):14-19

TONG S Y, XU L L, CUI S G, et al. Automatic identification and suppression of exogenous interference wave in the pre-stack seismic data of south Yellow sea[J]. Marine Geology Letters,2008,24(8):14-19.(In Chinese)

[5] 施劍，蔣龍聰. 南黃海地震資料鄰船干擾特征及壓制方法[J]. 海洋地質(zhì)前沿, 2014,30(7):66-70.

SHI J, JIANG L C. Characteristics of vessel noise in seismic survey and the method for noise suppression: a case from the South Yellow sea[J]. Marine Geology Frontiers, 2014,30(7):66-70.(In Chinese)

[6] 郭建卿，張保慶，趙貽水,等. 墨西哥灣海洋地震數(shù)據(jù)處理案例分析[J]. 石油地球物理勘探, 2007,42:132-135

GUO J Q, ZHANG B Q, ZHAO Y S, et al. Case analysis of marine seismic data processing in Mexico Bay[J]. Oil Geophysical Prospecting,2007,42:132-135.(In Chinese)

[7] 萬歡，但志偉，馮全雄. 海上地震勘探外源干擾快速壓制方法[J]. 工程地球物理學報, 2010,7(1):11-14.

WAN H, DAN Z W, FENG Q X. Outsource perturbation suppression in offshore seismic exploration[J]. Chinese Journal of Engineering Geophysics,2010,7(1):11-14.(In Chinese)

[8] 李景葉，陳小宏，張金淼,等. 海上點源干擾噪聲壓制方法[J]. 石油物探, 2013,52(1):23-29.

LI J Y, CHEN X H, ZHANG J M, et al. Point source diffracted noise attenuation in marine data processing[J]. Geophysical Prospecting for Petroleum, 2013,52(1):23-29.(In Chinese)

[9] 李新祥. 三維地震共炮檢距向量道集[J]. 石油物探，2007,46(6):545-549.

Li X X. 3D seismic common offset vector gather[J]. Geophysical Prospecting for Petroleum, 2007,46(6):545-549.(In chinese)

[10] 王在明，朱江，李磊,等. 振幅衰減法分頻自動識別和壓制異常能量道[J]. 新疆石油地質(zhì), 2011,32(1):80-82.

WANG Z M, ZHU J, LI L, et al. Automatic recognition and suppression of abnormal energy seismic Trace by method for amplitude attenuation and frequency division[J]. Xinjiang Petroleum Geology, 2001,25(4):290-293.(In Chinese)

[11] 陳志德，鄭錫娟，關昕，等. 長垣油田地震異常振幅噪聲分頻壓制技術[J]. 大慶石油地質(zhì)與開發(fā), 2010,29(5):164-167.

CHEN Z D, ZHENG X J, GUAN X, et al. Seperate frequency technique for anomalous amplitude noise of the seismic data in Daqing Placanticline[J]. Petroleum Geology and Oil Field Development in Daqing, 2010,29(5): 164-167.(In Chinese)

[12] 王君，周興元，惲春華. 地表一致性約束下異常振幅衰減技術的實現(xiàn)[J]. 石油物探, 2011,50(5):493-498.

WANG J, ZHOU X Y, YUN C H. Realization of surface consistent constrained anomaly amplitude attenuation technology[J]. Geophysical Prospecting for Petroleum, 2011,50(5):493-498.(In Chinese)

[13] 王世青，孟小紅，于春玲，等. 地震異常振幅自適應壓制[J]. 大慶石油地質(zhì)與開發(fā)，2006,25(5):112-115.

WANG S Q, MENG X H, YU C L，et al. Seismic abnormal amplitude self-adapting inhibition[J]. Petroleum Geology and Oil Field Development in Daqing，2006,25(5):112-115.(In Chinese)