歸平軍，吳永鋒，李　新

(中石化華北油氣分公司 勘探開發(fā)研究院，河南 鄭州 450006)

?

西藏羌塘盆地二維地震資料處理探討

歸平軍，吳永鋒，李新

(中石化華北油氣分公司 勘探開發(fā)研究院，河南 鄭州 450006)

西藏羌塘盆地地區(qū)在地層的演變形成過程中，由于地質(zhì)條件復雜，地下碳酸鹽巖目的層波阻抗差異小、界面反射信號弱。針對這些特點，探討了本地區(qū)地震資料處理方法，分析原始記錄及地質(zhì)特征，認識各種干擾波的特征與地表條件的關(guān)系，并采取相應的措施對地震數(shù)據(jù)進行靜校正、去噪、疊加等處理，取得了良好的效果，可為今后在此類地區(qū)進行地震資料現(xiàn)場處理提供借鑒。

靜校正；低信噪比；處理參數(shù)；西藏羌塘盆地

1　引　言

西藏羌塘盆地內(nèi)表層沉積結(jié)構(gòu)復雜多變，地形起伏劇烈，低降速帶的厚度和速度極不穩(wěn)定，同時存在強能量的面波和規(guī)則線性干擾波，資料的信噪比低。此時，若按照一般常規(guī)方法對地震資料處理，不能取得好的地質(zhì)效果。因此如何解決好此類地區(qū)的靜校正問題，有效地壓制各種干擾提高中深層的信噪比，是復雜地表條件下地震資料處理的一個關(guān)鍵。需要認真做好方法和參數(shù)試驗，優(yōu)化處理流程，精選處理參數(shù)，建立適合于本區(qū)資料的處理流程[1]。

針對以上特點，本文采用交互迭代折射靜校正方法進行靜校正，將多種去噪方法相結(jié)合，最大限度地壓制干擾波，迭后剩余靜校正，消除高頻剩余時差，在處理中取得了良好的應用效果。

2　原始記錄的特點

工區(qū)地貌以丘陵為主，個別地區(qū)為山地，地形起伏較大，區(qū)內(nèi)河灘發(fā)育，沼澤較多，個別地方有基巖出露地表；一般海拔高程在4 970～5 600 m左右；本區(qū)地表地質(zhì)條件復雜，低降速帶厚度變化不大，厚度在2～12 m左右，低速層速度在800～2 900 m/s，高速層速度在2 500～4 700 m/s[2]。

2.1幾種地表類型地震記錄

工區(qū)地表主要分為3類[3]：

1)山地區(qū)：地表高差變化大，地表基巖出露或覆蓋一層薄的第四紀覆蓋物，高速層為基巖。記錄品質(zhì)較差，信噪比較低，斷面反射及側(cè)反射互相干涉，嚴重影響有效波組的反射特征(圖1)。

2)平原區(qū)：地表較為平坦，為第四紀覆蓋物，高速層為基巖。記錄干擾以面波、折射多次波為主，3 s以上可見有效波，信噪比較高(圖2)。

3)河灘礫石區(qū)：主要分布河道或古河道，地表為小礫石，高速層為含水沙石。資料品質(zhì)較好，但面波、折射多次干擾嚴重(圖3)。

2.2干擾波分析

圖4為工區(qū)主要干擾波的分布圖。由圖4可知，本區(qū)干擾以面波、折射與折射多次波、高頻隨機噪音為主。其中，面波能量強，頻率高，在本區(qū)分布最為廣泛。面波干擾在記錄上分為兩組，第一組視速度約為1 200 m/s，頻率最高可達17 Hz，第二組視速度在700 m/s左右，頻率較低，約在10 Hz以下。

圖1　山地原始記錄品質(zhì)　Fig.1　Original record quality in the mountain

圖2　平原原始記錄品質(zhì)　Fig.2　Original record quality in the plain

圖3　河灘礫石區(qū)原始記錄品質(zhì)Fig.3　The original recording quality in beach gravel area

圖4　工區(qū)主要干擾波　Fig.4　The main interference wave in areas

折射與折射多次波以線性干擾形式出現(xiàn)在記錄的初至部分，速度較高，嚴重干擾剖面淺層及記錄中、遠偏移距的有效反射。環(huán)境噪音以高頻隨機噪音的形式出現(xiàn)在記錄上，其共有特征為干擾能量強，頻率高，隨機性強，嚴重影響了原始資料的品質(zhì)。

3　處理難點分析及對策

3.1影響資料品質(zhì)及處理工作的主要因素

區(qū)域相對高差大,地形起伏劇烈，溝壑縱橫;表層結(jié)構(gòu)復雜多變，是本工區(qū)表層的主要特征。地層起伏劇烈，褶皺嚴重，深部地層傾角大，產(chǎn)狀多變，逆掩斷層發(fā)育；地層切割性強，反射界面物性差,是本工區(qū)深層的主要特征。這些特征不僅是影響資料品質(zhì)的主要因素[4]，而且是各種干擾、次聲干擾產(chǎn)生的根源。

資料處理中多種問題交織在一起，相互制約、相互影響，極大地增加了處理的難度。歸納起來，主要有以下幾點問題：

1)復雜多變的表層，導致有效波組同相軸扭曲、錯斷(圖5)。

2)斷塊區(qū)內(nèi)繞射、斷面波、側(cè)反射、面波互相干涉，嚴重影響記錄的信噪比(圖6)。

3)低、降速帶速度、厚度的變化頻率遠遠大于折射點的分布頻率，因而工區(qū)的低降速帶資料不能滿足資料處理的要求，采用野外一次靜校正方法處理，其成果難以達到構(gòu)造解釋的要求(圖7)。

4)表層變化劇烈，深層干擾種類多、能量強、物性差，加劇了地震波、各種干擾能量在時、空間分布的復雜性，增加了處理工作中去噪、地震資料真振幅恢復的難度。

5)導致地震子波地表一致性變差，其振幅、頻率、相位發(fā)生畸變。

綜上分析，地震資料的信噪比、分辨率與工區(qū)的表層地震地質(zhì)條件息息相關(guān)，消除表層地震地質(zhì)條件變化劇烈所帶來的一系列問題是搞好處理工作的基礎(chǔ)。

3.2處理難點分析

根據(jù)對工區(qū)地震地質(zhì)條件的了解以及對原始資料深入細致的分析，結(jié)合以往對低信噪比地區(qū)資料的處理經(jīng)驗，筆者認為完成本次處理任務的關(guān)鍵在于解決好以下幾個重點環(huán)節(jié)[5]。

1)靜校正問題

特殊的沼澤—丘陵—山地地形，速度、厚度空間變化劇烈的表層，是本區(qū)的一大特征。為了保證良好的疊加效果，克服靜校不準的問題，需要選擇正確、合理的靜校方法。

2)“去噪凈化”問題

本區(qū)各類干擾發(fā)育，能量強，分布范圍廣，搞好本區(qū)地震資料的“去噪凈化”，改善資料的信噪比，是提高速度分析精度，保證疊加品質(zhì)的基礎(chǔ)。

圖5　淺層地震地質(zhì)條件對記錄同相軸的影響　Fig.5　The effect of shallow seismic geological conditions to the record of events

圖6　各種干擾對記錄同相軸的影響　Fig.6　The interference effect to the record of events

圖7　一次靜校正方法的缺陷Fig.7　The effect of a static correction method to the field

所以根據(jù)資料中的具體干擾特征選取最佳的去噪方法、加強去噪效果的監(jiān)控與分析是本次處理工作的一個重點。

3)地表一致性及分辨率問題

地震資料的保幅處理是后續(xù)巖性處理的基礎(chǔ)，消除非地質(zhì)因素對地震資料的影響，保護好地震反射波組的巖性特征是后續(xù)解釋處理成功的前提。隨著對表層因素影響的深入研究，人們越來越注意到，表層條件的變化不僅影響到記錄的時移，而且會使地震波的振幅、頻率、相位在傳播過程中均受到極大的改造，導致原始資料分辨率下降，地震資料受表層因素影響產(chǎn)生了一系列地表一致性問題。因此，消除表層因素的影響，恢復資料在空間上的真振幅特征，保證地震子波在空間上的一致性，進而達到壓縮子波，提高資料分辨率的目的，是處理中的關(guān)鍵環(huán)節(jié)。

3.3針對性的處理措施

3.3.1本區(qū)靜校正問題

針對該區(qū)地形起伏，低降速帶速度、厚度劇烈變化所引發(fā)的靜校正難題，此次處理項目擬采用層析靜校正方法求取長、中波長的靜校正量；地表一致性靜校正解決短波長的剩余靜校正量問題。處理中通過認真細致地初至拾取、建模及嚴密的質(zhì)量監(jiān)控手段，力爭在本區(qū)靜校正工作上取得良好的效果。

3.3.2多系統(tǒng)、多方法的聯(lián)合去噪凈化工作

選用不同處理系統(tǒng)的最佳去噪手段組合應用，對疊前地震資料進行充分去噪靜化，可有效地提高疊前數(shù)據(jù)信噪比，為資料的后續(xù)處理提供有力的保證。具體措施如下：

①充分發(fā)揮FOCUS系統(tǒng)良好的交互處理功能進行資料的壞炮、干擾道編輯及其處理效果、處理質(zhì)量的監(jiān)控工作；②根據(jù)頻譜分析和濾波掃描選取合理的濾波參數(shù)來壓制高能量的高、低頻隨機干擾；③針對尖脈沖干擾，采用小時窗、大門檻值的強脈沖能量自動壓制手段給予充分的壓制；④根據(jù)面波的頻率、頻散特性，通過時頻分析，采用自適應濾波方法，壓制記錄中不同頻率成分的面波能量；⑤采用多道識別、單道壓制的線性濾波方法消除折射、折射多次波等線性干擾，避免多道處理的混波、假頻等副作用；⑥對于部分原始資料中的聲波干擾，采用交互切除的方法，在反褶積后進行時窗切除，以避免截斷處理噪音的影響。

3.3.3地表一致性、分辨率問題

為校正本區(qū)表層地震地質(zhì)條件的劇烈變化的導致的地震波振幅、頻率、相位在時、空上的畸變、衰減，必須進行地表一致性處理。首先，擬應用時頻、頻時分析的手段，詳細掌握資料的振幅、頻率在時、空方向的變化規(guī)律，并根據(jù)不同深度地層在速度上的變化信息，在疊前采用能夠保留波組之間相對強弱特征的變速球面擴散補償來調(diào)整疊前數(shù)據(jù)能量；根據(jù)資料的振幅統(tǒng)計特征進行地表一致性振幅補償、地層反射、透射振幅補償，力求準確、合理地恢復記錄在時、空間上的能量損失，并通過保幅的整道均衡，使淺、中、深層波組間能量處在同一能量級上；疊后在試驗的基礎(chǔ)上，通過增益曲線分析、時頻分析等進行剩余振幅分析與補償，使淺、中、深層能量均一，波組特征清楚、強弱特征明顯。然后，采用多道地表一致性反褶積對子波整形，消除表層因素對子波的影響，改善資料的頻率和相位特征。在兼顧信噪比的原則下，壓縮子波，適度提高剖面的分辨率，并在此基礎(chǔ)上采用最佳信噪比濾波進一步改善資料的信噪比。最后，采用地表一致性剩余靜校，消除受表層因素影響產(chǎn)生的剩余時差，改善疊加效果。

4　流程試驗與效果

4.1靜校正試驗與效果

針對該區(qū)地形起伏，低降速帶速度、厚度劇烈變化所引發(fā)的靜校正難題，此次處理項目擬采用層析靜校正方法求取長、中波長的靜校正量。通過認真細致的初至拾取、建模及嚴密的質(zhì)量監(jiān)控手段(圖8)，在本區(qū)靜校正工作上取得良好的效果[7]。

圖9是本區(qū)原始炮記錄與折射靜校正的效果對比圖。其中圖9(a)為原始炮記錄，圖9(b)為折射靜校正的效果。由圖9可見，折射靜校正在本區(qū)應用有明顯的效果， 折射靜校正記錄有效波同相軸光滑、連續(xù)性好。

圖8　初至模型反演效果Fig.8　The first break model inversion effect

圖9　靜校正方法試驗Fig.9　The static method test

4.2“去噪凈化”試驗與效果

針對干擾波不同的特征，采用不同的方法進行疊前去噪工作。

圖10為傾角濾波去噪效果圖。由圖10可知，針對地震記錄中折射多次波嚴重影響資料中、遠道的信噪比的問題，采用傾角濾波壓制手段取得了較好的效果。

圖11為面波自適應壓制效果對比圖。其中，圖11(a)為原始記錄，圖11(b)為面波自適應壓制效果。從圖11可以看到，高通濾波壓制面波易傷及低頻有效波，且近道效果差；自適應面波壓制根據(jù)面波與有效波在頻率、速度等方面的差異，自動識別壓制，壓制效果好，無混波作用。記錄信噪比進一步得到提高，有效波能量更加突出。

圖10　傾角濾波去噪效果　Fig.10　The dip filtering denoising effect

圖11　面波自適應壓制效果　Fig.11　The surface wave adaptive suppression effect

4.3振幅補償試驗與效果

根據(jù)地震波能量在傳播中損失的各個原因，有針對性地進行了球面擴散補償和炮域、檢波域與共偏移距域的地表一致性振幅補償，用于補償?shù)卣鸩▊鞑r在時、空間的能量衰減[8]。

圖12圖為振幅補償效果圖。其中，圖12(a)是補償前的炮記錄，它在時空間上的能量衰減是很明顯的，經(jīng)過球面擴散補償以后，地震波深部能量得到了有效補償(圖12b)。圖12(c)是炮域、檢波域與共偏移距域的地表一致性振幅補償效果，記錄能量通過兩種有效補償在時、空間分布趨于真實的能量反射。

4.4反褶積處理試驗與效果

在充分去噪的基礎(chǔ)上，應用地表一致性預測

反褶積在兼顧信噪比的原則下，通過子波整形，壓縮子波，消除表層因素對子波的影響，適度提高資料的分辨率。圖13為反褶積前、后子波壓縮效果對比。在此基礎(chǔ)上采用最佳信噪比濾波，通過壓制反褶積引起的甚高、低頻的噪音，加強了優(yōu)勢頻帶能量，改善了記錄的信噪比。圖14從左到右依次為去噪后記錄、地表一致性反褶積及最佳信噪比濾波處理效果。

4.5剩余靜校正處理試驗與效果

在靜校正處理的基礎(chǔ)上，為了進一步提高剖面同相軸的連續(xù)性，改善疊加效果， 剩余靜校正處理是不可或缺的一環(huán)，為此，進行了常規(guī)的兩次地表一致性剩余靜校正處理試驗， 提高了反射波同相軸的連續(xù)性，改善了剖面的信噪比(圖15)。

圖12　振幅補償效果Fig.12　The amplitude compensation effect

圖13　反褶積前、后炮記錄子波整形、壓縮效果Fig.13　The shot record of wavelet shaping and compression effect before and after deconvolution

圖14　處理效果分析Fig.14　The processing effect analysis

圖15　自動剩余靜校正對比試驗效果Fig.15　The automatic residual statics contrast test

4.6處理效果分析

圖16為疊加效果分析圖。由圖16可見，處理剖面淺、中、深層反射波組齊全，同相軸連續(xù)性好；主要目的層波組特征清晰，信噪比、分辨率較高，易于追蹤對比。處理剖面斷塊間斷點清楚，繞射波連續(xù)性好；各斷塊反射波組特征清晰，信噪比較高，易于對比解釋。而采用常規(guī)處理方法很難解決靜校正問題及在低信噪比地區(qū)的成像問題，采取本方法可以有效地解決這些問題，達到預期的效果。從剖面總體品質(zhì)看，整條剖面構(gòu)造格局明顯、特征清楚，較好地完成了揭示測區(qū)盆地基底及其上覆沉積層構(gòu)造格架的地質(zhì)任務。

圖16　測線最終疊加剖面效果Fig.16　The final stacked section effect of the survey line

5　結(jié)　語

由于本區(qū)特殊的地震地質(zhì)條件，原始資料信噪比低，各種干擾能量強，分布廣，如何進一步提高原始資料的品質(zhì)是進行地震資料處理的前提。低信噪比地區(qū)資料的處理必須要有疊加次數(shù)作為依托，以保證淺、中層能夠有效的疊加。合適疊加次數(shù)對壓制噪音、提高信噪比有很大的作用。本文通過方法試驗，采用合理的處理流程，在該區(qū)獲得了一個反射波組齊全、主要目的層信息豐富、地質(zhì)構(gòu)造特征清晰的地震剖面，為今后在此類地區(qū)進行地震資料現(xiàn)場處理提供了借鑒。

[1]劉廣信，李新民.CYBER170_720寬線剖面資料處理及效果[J].石油物探，1984,23(4):39-45.

[2]盧占武，高銳, 匡朝陽,等.青藏高原羌塘盆地二維地震數(shù)據(jù)采集方法試驗研究[J].地學前緣，2006，13(5):382-389.

[3]盧占武,高銳,李秋生,等.藏北羌塘盆地反射地震剖面與疊加速度研究[J].中國地質(zhì);2009,36(3):677-681.

[4]翟桐立，劉次源，祝文亮,等.南方山地寬線地震采集方法與效果[J].天然氣工業(yè)，2007,27(增刊A):70-71.

[5]朱漢東.沖斷山區(qū)的地震資料采集[M]//石油物探專題情報成果集.保定:石油部物探局情報所,1997.

[6]張金崗，劉澤源，胡杰,等.柴達木盆地復雜山地實用地震技術(shù)研究[J].青海油田，2006,24(3)：1-8.

[7]張春賀，喬德武，李世臻,等.復雜地區(qū)油氣地球物理勘探技術(shù)集成[J].地球物理學報,2011,54(2):374-387.

[8]陸基孟.地震勘探原理[M].北京:石油工業(yè)出版社,1993.

The Discussion of 2D Seismic Data Processing in Qiangtang Basin of Tibet

Gui Pingjun, Wu Yongfeng, Li Xin

(PetroleumExplorationandProductionResearchInstitute,theNorthChinaOilandGasCompanyofSinopec,ZhengzhouHenan450006,China)

In the process of stratigraphic formation in Qiangtang Basin region, and due to the complex geological condition, underground carbonate target layer wave impedance difference is small, and the interface reflection signal is weak. Based on these features, the regional seismic data processing method is studied, the original data and geological characteristics are described, the characteristics of various interference waves and the relationship of the surface conditions are analyzed. Besides, some measures are taken like the seismic data static correction, de-noising, stacking and other processing methods, which has achieved good results.

static; low signal-to-noise ratio; processing parameter; Qiangtang basin

1672—7940(2016)04—0513—08

10.3969/j.issn.1672-7940.2016.04.019

歸平軍(1986-)，男，工程師，主要從事地震處理及解釋的工作。E-mail:pingjun2004@126.com

P631.4

A

2015-11-26