秦永軍，馬　麗，2，薛海軍，2，孫文華

（1.陜西省煤田物探測繪有限公司，西安710005；2.國土資源部煤炭資源勘查與綜合利用重點實驗室，西安710026）

地震波形分類技術(shù)在煤層分叉解釋中的應(yīng)用

秦永軍1，馬麗1，2，薛海軍1，2，孫文華1

（1.陜西省煤田物探測繪有限公司，西安710005；2.國土資源部煤炭資源勘查與綜合利用重點實驗室，西安710026）

地震波形分類技術(shù)是地震相分析技術(shù)的重要延伸。通過對時窗選取、地震屬性體選擇、波形分類數(shù)、模型道選取等主要技術(shù)參數(shù)的分析，認(rèn)為時窗選取和地震屬性體是影響波形分類結(jié)果最大的兩個因素。以河?xùn)|北部煤田為例，利用地震波形分類技術(shù)對其13#煤層的分叉現(xiàn)象進(jìn)行了分析，在此基礎(chǔ)上，對13#煤反射波上下各20 ms的層段，進(jìn)行了三次計算，經(jīng)對比認(rèn)為當(dāng)波形分類數(shù)為6時的地震相可很好的區(qū)分煤層分叉的不同狀態(tài)。對比振幅、頻率及相位等常規(guī)屬性與分類結(jié)果的地震相，發(fā)現(xiàn)常規(guī)屬性平面雖也能反映出煤層分叉趨勢，但屬性平面圖對分叉線表達(dá)欠清晰，而三種屬性的波形分類平面圖可清晰的反映出煤層從西南至東北方向的分叉邊界。

波形分類、煤層分叉、模型道、地震相

0　引言

隨著地震勘探技術(shù)的不斷提高，由地震數(shù)據(jù)獲得的地質(zhì)信息越來越多，但其中許多信息難以得到直觀的結(jié)果，需通過數(shù)學(xué)變換、并在已知信息標(biāo)定下得到較明了的地質(zhì)成果。近年來發(fā)展的地震相分析技術(shù)，是在劃分地震層序的基礎(chǔ)上，利用地震參數(shù)特征差別，將地震層序劃分為不同的地震相區(qū)，然后作出巖相和沉積環(huán)境的推斷。地震相技術(shù)目前已經(jīng)推廣到巖性及儲層的預(yù)測中，常采用的參數(shù)有反射振幅、反射頻率、同相軸連續(xù)性等。

地震波形分類技術(shù)是地震相分析技術(shù)的一種，是通過分析地震道間的振幅、頻率、層速度、連續(xù)性等屬性的異常來分析沉積相的變化。目前該項技術(shù)在油氣藏勘探領(lǐng)域已經(jīng)成功地應(yīng)用于油氣層分布預(yù)測、砂體儲層預(yù)測、碳酸鹽巖礁灘體巖性識別等方面［1-3］。程增慶、劉天放等人指出，煤層厚度變化會引起地震反射波振幅、頻率等屬性的變化［4］，當(dāng)煤層沉積穩(wěn)定時屬性參數(shù)變化平穩(wěn)；而當(dāng)煤層分叉時，其沉積環(huán)境和地層巖性組分也表現(xiàn)出明顯的不同［5］，使地震波形分類技術(shù)對煤層分叉這一地質(zhì)現(xiàn)象的預(yù)測成為可能。本文充分利用豐富的地震信息，將地震波形變化定量刻畫出來，對目的層段波形逐道對比分析，得到不同分類的地震波形橫向變化的地震相平面圖［6］，從面實現(xiàn)預(yù)測煤層分叉的目的。

1　波形分類技術(shù)的原理

波形分類技術(shù)通過提取地震信息在空間上的相似性來描述地質(zhì)信息的空間變化。在特定目的層時窗內(nèi)，根據(jù)地震道波形特征，逐道對比，求同存異，突出各地震道的相似性，刻畫地震波形在橫向上的變化，最后在平面上以不同顏色代表不同地震相類型，以反映沉積相變化，從而認(rèn)識沉積相在平面上的分布規(guī)律。

通常地震波形分類技術(shù)主要是應(yīng)用Manhattan距離公式。Manhattan距離是一種求兩點之間的相似度和非相似度的統(tǒng)計方法，表征了波形樣點之間的采樣差異值的總和［7-8］。Manhattan距離形式如下：

M=總值/標(biāo)準(zhǔn)測量值

M為Manhattan距離，A為參考子波，B為目標(biāo)子波，i為子波的采樣數(shù)。

當(dāng)兩個子波完全相同時將導(dǎo)致Manhattan距離為零，不完全相同的子波將產(chǎn)生正的不同的Manhat?tan距離。如果參考子波和目標(biāo)子波的M值是0.25或更小，那么那些子波被分在同一個組或類。根據(jù)這一原則，就可以用參考子波確定分類，然后將周圍的目標(biāo)子波分到不同的參考子波組。如果一個目標(biāo)子波屬于更多的參考子波組，將目標(biāo)子波分配給M值最小的一組。這樣就可將地震波形進(jìn)行分類。

本次波形分類技術(shù)主要依靠的波形包括地震波的振幅、頻率、相位等屬性的綜合特征。根據(jù)波形分類結(jié)果，各個地震道形成離散的“地震相”。對這些離散的“地震相”進(jìn)行平面歸類，得到平面地震相圖。通過已知的地質(zhì)信息進(jìn)行標(biāo)定，對波形分類的結(jié)果進(jìn)行綜合地質(zhì)解釋，最終得出與地震相分布相對應(yīng)的沉積相圖，即由地震相轉(zhuǎn)化為沉積相。

2　波形分類主要技術(shù)參數(shù)

波形分類的質(zhì)量取決于時窗的選取、地震屬性體的選擇、波形分類數(shù)，波形計算中的迭代次數(shù)等主要技術(shù)參數(shù)的選擇。其中，時窗的選取和地震屬性體是影響波形分類結(jié)果最大的兩個因素。

2.1時窗的選取

地震時窗是在兩個層位之間或某個層位上、下限范圍的地震數(shù)據(jù)的集合［9］。時窗的選取是波形對比最為敏感的參數(shù)，合理的時窗才可以保證分類結(jié)果能真實反映目標(biāo)層位的特征。時窗的選擇有兩種：一種是固定時窗；另一種是變時窗。本次是以固定時窗來進(jìn)行研究，選取的時窗不僅保證在選定的時窗格架內(nèi)，避免出現(xiàn)穿時現(xiàn)象，也要完整地涵蓋一個地質(zhì)單元，目標(biāo)地質(zhì)特征層段在選定時窗格架內(nèi)，才能獲得一個較好的效果。

2.2地震屬性體

波形分類技術(shù)采用的數(shù)據(jù)體并不局限于原始地震數(shù)據(jù)，還可以使用轉(zhuǎn)換過的地震屬性體，即針對不同地震屬性體的波形分類對比［10］?；诘卣饘傩泽w的波形分析，其本質(zhì)是通過不同的數(shù)學(xué)變換突出地震信號的某些特征，其目的是為了突出地震信息中的有利信息，削弱干擾信息。如：轉(zhuǎn)化為振幅類屬性體，使其在地震振幅上的差異變得更明顯；轉(zhuǎn)換為頻率類的屬性，使其在頻率域上的區(qū)分更加突出。盡管地震屬性體綜合了頻率、相位、振幅等綜合信息，但通過相關(guān)的數(shù)學(xué)變換可使振幅、頻率等波形特征在波形對比中處于優(yōu)勢地位。針對同一地質(zhì)體，采用不同的地震屬性體進(jìn)行波形分類可以得到不同的結(jié)果。在地震信噪比低或斷裂較發(fā)育地區(qū)，可以通過適當(dāng)?shù)膶傩泽w轉(zhuǎn)換來進(jìn)行波形分類，弱化無關(guān)的噪聲和雜亂地震信息對波形分類所帶來的影響，強(qiáng)調(diào)不同巖相主要波形特征在波形對比中的權(quán)重。

2.3提取模型道

在提取模型道時，初始劃分幾種已知的典型形狀，然后對每一實際道賦予一個基于相似性的典型形狀。在目標(biāo)層段內(nèi)與實際地震數(shù)據(jù)進(jìn)行對比，通過自適應(yīng)和誤差分析，最后使模型道和實際地震道之間存在更好的相關(guān)性。

2.4波形分類數(shù)的選取

波形分類數(shù)的選擇是指在整個目標(biāo)層段內(nèi)所遇到的地震道的種類，是影響波形分類結(jié)果較小的一個因素。最為理想的分類數(shù)是不容易定義的。至少經(jīng)過三次的計算分析。根據(jù)經(jīng)驗估算方法：①把目標(biāo)層段的厚度除以6作為第一次計算的分類數(shù)；②把上次計算的分類數(shù)的50%作為第二次計算的分類數(shù)；③把第一次計算分類數(shù)150%作為第三次計算的分類數(shù)［11］。最佳分類數(shù)應(yīng)根據(jù)研究目標(biāo)與地震數(shù)據(jù)的實際情況確定，分類數(shù)越大，結(jié)果過于詳細(xì)復(fù)雜，分類數(shù)越小，結(jié)果過于粗糙無實際地質(zhì)意義，一般情況波形的分類數(shù)在6～12。

3　應(yīng)用效果

3.1研究區(qū)概況

研究區(qū)位于河?xùn)|煤田北部，整體為走向近南北，傾向西的單斜構(gòu)造，區(qū)內(nèi)大面積被第四系黃土層覆蓋，僅在中部溝谷出露新近系保德組，區(qū)內(nèi)賦存地層有奧陶系上馬家溝組，石炭系本溪組、太原組，二疊系山西組、下石盒子組，新近系及第四系。含煤地層為二疊系山西組及石炭系太原組。主采煤層有山西組3號煤層，太原組11號、13號煤層，其中13號煤層在研究區(qū)存在煤層分叉現(xiàn)象。

13號煤層煤厚度為6.6～14.45 m，平均煤厚10.52 m。在研究區(qū)西部13號煤層逐漸分為兩層。整體煤厚呈現(xiàn)由西到東逐漸變厚的趨勢，煤層結(jié)構(gòu)復(fù)雜，含夾石1～4層，夾石厚度為0.30～0.50 m，巖性為泥巖或炭質(zhì)泥巖，頂板底板均為泥巖或砂質(zhì)泥巖。煤層可采性指數(shù)Km=1.00，煤厚變異系數(shù)γ= 20.20%，屬全區(qū)穩(wěn)定可采的特厚煤層。

3.2主要技術(shù)參數(shù)

通過分析對比，以標(biāo)定的13煤反射波上下各20 ms為目標(biāo)層段（上下均不超過相鄰煤層反射波）。經(jīng)過三次計算，發(fā)現(xiàn)當(dāng)波形分類數(shù)為6時，地震相較為合理，能很好地區(qū)分煤層分叉的不同狀態(tài)。最終在該時窗范圍內(nèi)，針對不同典型分叉狀態(tài)選擇提取6類不同的模型道（圖1），相應(yīng)的地震波形特征如圖2。從模型（圖1）波形特征（圖2）可以看出1-2類模型道及波形完整，能量強(qiáng)，代表煤層未分叉前的煤層；3-4類模型道及波形出現(xiàn)復(fù)波，頻率較高，能量變?nèi)?，代表分叉過渡段的煤層；5-6類模型道及波形頻率降低，能量較強(qiáng)，代表分叉后的煤層。

3.3煤層分叉解釋與效果分析

3.3.1煤層分叉解釋

圖1　6類模型道Figure1　Six kinds of model trace

圖2　T13波±20ms層段不同類模型道的波形特征Figure2　Waveform features of different model traces in an interval above and below 20ms each from T13 wave

圖3　常規(guī)地震屬性與波形分類結(jié)果Figure3　Conventional seismic attributes and waveform classified results

按照地震不同屬性（振幅、頻率及相位）的常規(guī)屬性與分類結(jié)果的地震相進(jìn)行對比，見圖3，從圖中可以看出，常規(guī)屬性平面與波形分類均能表達(dá)出煤層分叉趨勢為走向NE的條帶狀，但是屬性平面圖對分叉線的表達(dá)不清晰，尤其是振幅屬性僅表現(xiàn)為平緩地過渡；相位屬性表達(dá)效果略好，有部分地段表達(dá)也不清晰；頻率屬性表達(dá)效果最好，但是三種屬性表達(dá)的分叉線形態(tài)與平面擺布差異大，不能用于解釋分叉位置。比較而言，三種屬性的波形分類平面圖對煤層分叉的表達(dá)效果很直觀，三者具有一定的統(tǒng)一性，都能清晰地指示出煤層從西南至東北方向的分叉邊界走向，受地震資料品質(zhì)影響，在東北部煤層分叉連續(xù)性變差，但趨勢清晰，說明波形分類對煤層分叉解釋較好。綜合不同屬性，形成地震分類平面圖（圖4），其中色標(biāo)1～1.5為煤層分叉區(qū)域，圖中灰線為常規(guī)地震屬性解釋的煤層分叉線，黑線為波形分類解釋的煤層分叉線，可以看到黑線向前推進(jìn)5～60 m，相對于常規(guī)地震解釋的分辨率限制，波形分類解釋結(jié)果更符合實際地質(zhì)規(guī)律，相對更接近真實分叉線，控制也更為精細(xì)。

3.3.2分類解釋煤層分叉區(qū)效果分析

圖4　波形分類解釋成果示意圖Figure4　A schematic diagram of waveform classification interpreted results

通過波形分類后可以看出區(qū)內(nèi)分為三種地震相：分叉相、過渡相和合并相。從目前施工的鉆孔來看合并相內(nèi)的SD-11孔及SD-12孔煤層厚度分別為12.6 m、13.5 m，均含煤一層，且厚度較大；南部巷道揭露13煤均為一層煤，厚度為13 m左右。過渡相內(nèi)的SD-07孔煤層結(jié)構(gòu)為兩層，煤厚分別為6.9 m、4.7 m，煤間距5.2 m，煤層開始分叉；分叉相內(nèi)的SD-03鉆孔含煤兩層，煤厚分別為6.6 m、6.8 m，煤層間距33.5 m，煤層完全分開。地震波形分類解釋成果得到了鉆孔的驗證。

4　結(jié)論

通過對研究區(qū)煤層分叉的實踐，證實了地震波形分類技術(shù)可用于煤層分叉解釋，解釋結(jié)果比常規(guī)地震解釋成果更直觀、更準(zhǔn)確，為地震波形分類技術(shù)在煤田勘探中煤層宏觀結(jié)構(gòu)的變化研究提供了一定參考。在地震波形分類應(yīng)用過程中，應(yīng)注意以下幾點：

（1）在煤層分叉解釋中，時窗的選擇不宜過大，薄-厚煤層一般選擇一個完整的波形時窗即可；對于特厚煤層（厚度大于8 m），可以根據(jù)情況選擇地震波形的一至幾個周期，應(yīng)盡量避免穿時現(xiàn)象發(fā)生。

（2）典型波形分類的選擇應(yīng)以地質(zhì)鉆井的測井曲線為約束，確保波形分類的正確性。

（3）波形分類技術(shù)與地震資料品質(zhì)有密切聯(lián)系，適于地震資料品質(zhì)好的地區(qū)。本次研究區(qū)東北部的地震數(shù)據(jù)品質(zhì)較差，導(dǎo)致這一區(qū)域內(nèi)波形分類結(jié)果跨度較大、連續(xù)性較差，這樣目標(biāo)波形可能會匹配不同的模型道，盡管最終還是分配給距離最小的一組，但仍對波形分類結(jié)果有一定影響。

［1］任芳祥，孫巖，樸永紅，等.地震技術(shù)在油田開發(fā)中的應(yīng)用問題探討［J］.地球物理學(xué)進(jìn)展，2010，25（1）：4-8.

［2］王玉學(xué)，叢玉梅，黃見，等.地震波表分類技術(shù)在河道預(yù)測中的應(yīng)用［J］.資源與產(chǎn)業(yè)，2006，8（2）：71-75.

［3］李雷濤，肖秋紅.地震波形分類技術(shù)在地震相分析中的應(yīng)用——以清溪場地區(qū)為例［J］.天然氣技術(shù)，2010，4（5）：4-7.

［4］董守華，馬彥良，周明.煤層厚度與振幅、頻率地震屬性的正演模擬［J］.中國礦業(yè)大學(xué)學(xué)報，2004，33（1）：29-32.

［5］張玉峰.朔縣平原分叉煤層對比方法的探討［J］.中國煤炭地質(zhì)，2008，20（9）：23-24.

［6］鄭文召.波形分類技術(shù)在王146地區(qū)地震相分析中的應(yīng)用［J］.中國石油大學(xué)勝利學(xué)院學(xué)報，2014，28（2）：5-7.

［7］劉海燕，劉財，劉洋，等.基于Manhattan距離的相干體技術(shù)的應(yīng)用［J］.世界地質(zhì)，2013，32（1）：144-146.

［8］范洪軍，范廷恩，王暉，等.地震波形分類技術(shù)在河流相儲層研究中的應(yīng)用［J］.CT理論與應(yīng)用研究，2014，23（1）：71-80.

［9］江青春，王海，李丹，等.地震波形分類技術(shù)應(yīng)用條件及其在葡北地區(qū)沉積微相研究中的應(yīng)用［J］.石油與天然氣，2012，33（1）：135-140.［10］洪忠，張猛剛，蘇明君，等.應(yīng)用地震波形分類技術(shù)識別巖相的適用性和局限性［J］.物探與化探，2013，37（5）：904-910.

［11］鄧傳偉，李莉華，金銀姬，等.波形分類技術(shù)在儲層沉積微相預(yù)測中的應(yīng)用［J］.石油物探，2008，47（3）：262-265.

Application of Seismic Waveform Classification Technology in Coal Seam Bifurcation Interpretation

Qin Yongjun1，Ma Li1，2，Xue Haijun1，2and Sun Wenhua1
（1.Shaanxi Provincial Coal Geophysical Prospecting，Surveying and Mapping Co.Ltd.，Xi'an Shaanxi 710005；2.Key Laboratory of Coal Resource Exploration and Comprehensive Utilization，Ministry of Land and Resources，Xi'an，Shaanxi 710026）

The seismic waveform classification technology is a major development of seismic facies analysis technology.Through analy?ses of major technical parameters including time window designation，seismic attribute cube selection，number of waveform classified results，model trace extraction，considered that the former two are the major factors impacting waveform classified results.Taking the Hedong north coalfield as example，through seismic waveform classification technology，analyzed the bifurcation phenomenon of the coal No.13.To an interval above and below 20ms each of reflection wave，after 3 computations，when classified waveform number is 6，seismic facies can discriminate different states of coal bifurcation well.Comparison of conventional attributes of amplitude，frequency and phase with classified result seismic facies have found that although conventional attributes plan can reflect coal bifurcation trend，but expression of attribute plan is not clear，while the waveform classification plan of three attributes can reflect coal bifurcation bound?ary from SW to NE clearly.

waveform classification；coal seam bifurcation；model trace；seismic facies

P631.4

A

10.3969/j.issn.1674-1803.2016.10.17

1674-1803（2016）10-0076-05

陜西省科技統(tǒng)籌項目資助（編號：2014KTZB01-03-03；2016FWRT-16）

秦永軍（1987—），男，工程師，畢業(yè)于西安石油大學(xué)，從事煤田地質(zhì)勘探工作。

2016-08-12

責(zé)任編輯：孫常長