（河南省煤炭地質勘察研究總院，河南鄭州市，450009） 郝 珺

（河南省地質礦產開發(fā)局第一地質勘察院，河南鄭州市，450009） 劉登鋒

（河南工程學院，河南鄭州市，451191） 黃 波

近年來，隨著煤炭資源的持續(xù)穩(wěn)定開采，接替采區(qū)已有向深部或者外圍拓展的趨勢，三維地震（測井）被廣泛應用于煤層預測與評價[1-4]。通過三維地震在礦井生產和補充勘探中的應用，可以較為準確的解釋礦井新生界厚度、構造、陷落柱的發(fā)育、煤層底板標高和煤層賦存情況等[5-9]。而煤厚變化(包括薄煤帶圈定)解釋精度對于礦井開拓尤為重要[10-11]。

煤層厚度變化預測主要有構造分析，鉆孔對比、內插法等[12]，受地質條件及鉆孔數量限制，這些方法要么定性預測，要么預測的煤厚誤差較大。因地質條件的復雜性和三維地震的多解性，常規(guī)地震解釋成果往往不理想[10-11]。采用三維地震方法可提供高精度的含煤地層地質信息，基于提取的地質信息，可較為準確的圈定煤層埋深、厚度變化、甚至薄煤帶等[12]。以趙樓煤礦八、十采區(qū)補充勘探為例，對趙樓煤礦三號煤層埋深、煤層厚度以及薄煤帶進行解釋。

1 方法原理

根據由已知到未知、由簡單到復雜、由點到面、由大到小的原則，充分利用解釋工作站的靈活、高效、直觀的特性，綜合已知地質、測井資料、屬性技術進行綜合解釋，在分析勘探區(qū)地質、地震資料特點的基礎上，制定適合勘探區(qū)的三維資料解釋流程（圖1）。解釋過程中，利用縱向、橫向和任意方向時間剖面相結合，時間剖面和水平切片、順層切片相結合等方法，全方位的反復對比、反復檢查、反復修改確認，確保解釋結果的正確可靠[12]。

圖1 資料解釋流程圖

1.1 多屬性提取

地震屬性是由疊前或疊后地震數據經過數學變換而導出的有關地震波幾何形態(tài)、運動學、動力學和統(tǒng)計特征的可定量化表征參數[13]。根據Widess的薄層理論，當薄層厚度小于λ/4 時（λ 反射波主波長），薄層反射波波峰與波谷視時差近似為一個常數，而反射波振幅隨薄層厚度呈準線性變化，即薄層的厚度信息包含在反射波振幅之中[14]。

利用目前功能強大的的Geoframe 及VVA（VisualVaxAt）解釋軟件，提取了常用的12種地震屬性，即最大振幅；最大數量；正極性振幅和；瞬時相位；瞬時頻率；平均數量；平均波峰數；均方根振幅；弧長；最小振幅。

經優(yōu)化選擇了對薄煤區(qū)反映較好的地震屬性參與了解釋，最終通過對優(yōu)選的屬性進行多元回歸計算，實現對煤層厚度預測及薄煤帶范圍圈定[13]。

1.2 波阻抗反演

常規(guī)的地震波阻抗反演就是指利用地震資料反演地層波阻抗的特殊地震處理技術。該方法有效地將測井資料垂向分辨率高的特點和地震資料橫向分辨率高的特點結合起來，達到識別巖性和煤層的目的。

波阻抗反演主要流程為：①對初始模型的波阻抗進行方形濾波；②用塊狀化的波阻抗和子波進行褶積合成地震記錄；③將合成的地震記錄和真實地震記錄進行比較；④根據步驟3 中的殘差改變塊狀波阻抗提高其與真實波阻抗的擬合程度；⑤重復若干次步驟3、4 最終得到反演結果(一般要求殘差小于5%）；⑥反演地質參數，預測煤層分布。

2 應用實例

2.1 補充勘探區(qū)概況

趙樓煤礦位于山東省鄆城縣城東南約22km，巨野縣城西13km。行政區(qū)劃歸鄆城縣、巨野縣管轄?？碧絽^(qū)位于趙樓煤礦八、十采區(qū)。

勘探區(qū)煤系地層以奧陶系為基底，沉積了石炭系本溪組、石炭-二疊系太原組、二疊系山西組及石盒子組，其上被新近系和第四系所覆蓋。主要含煤地層為太原組和山西組，兩組地層平均總厚242.46m，含煤25 層。其中山西組含煤3 層（1、2、3（3上、3下））；太原組含煤22層?？碧絽^(qū)主采煤層為3煤層，分上、下分層，屬穩(wěn)定煤層。頂板多為泥巖、細砂巖，局部見中砂巖；底板多為泥巖、粉砂巖。煤層與頂底巖石板密度差異較大，波阻抗分異明顯，形成了一個良好的波阻抗界面，地震條件良好。

2.2 補充勘探區(qū)3煤層圈定

3煤厚度穩(wěn)定，與上下頂底板圍巖間有明顯的波阻抗差異，在時間剖面上,所形成的反射波T3波能量強，連續(xù)性好,全區(qū)可連續(xù)追蹤。本次除了采用常規(guī)的時間剖面解釋外，還綜合利用地震資料多種屬性并結合區(qū)內鉆孔揭示的煤層厚度數據，以及采用Strata 軟件進行了疊后反演（縱波波阻抗），對勘探區(qū)煤厚變化趨勢進行了預測。繪制了煤層厚度變化趨勢平面圖（如圖2）。

對于薄煤層，在時間剖面上主要表現為煤層同相軸反射波變弱、消失(圖3)。波阻抗疊后反演切片圖上，排除巖漿巖影響后，薄煤帶解釋依據為波阻抗值一般大于7 000(g/cm3)·(m/s)，以黃色為主（如圖4）。屬性分析以分析以均方根振幅屬性為例，薄煤帶與均方根振幅呈負相關變化，屬性值越小為薄煤帶的可能性越大，以紅色、黃色為主（如圖5）。經綜合分析，勘探區(qū)內3 煤厚度在0.5～9m 之間,大部分在4～8m 之間。最厚處位于勘探區(qū)中部，約9m。在勘探區(qū)內解釋了9 處薄煤區(qū)，總面積94 183m2。（如圖2）

圖2 3號煤煤厚及薄煤帶分布圖

圖3 薄煤區(qū)在時間剖面上的顯示

圖4 3號煤層波阻抗切片圖

圖5 3號煤層均方根振幅屬性圖

3 結論

論述了常規(guī)的時間剖面、波阻抗反演技術、屬性提取技術的處理流程，并以趙樓煤層八、十采區(qū)為實例進行了應用研究，解釋了煤層厚度并圈定了薄煤帶。

依據解釋成果，勘探區(qū)內3 煤厚度在0.5～9m之間,大部分在4～8m 之間。最厚處位于勘探區(qū)中部，約9m。在勘探區(qū)內解釋了9 處薄煤區(qū)，總面積94 183m2。

三維地震綜合解釋過程中，應獲取較高質量的原始野外數據，對屬性提取技術應做好保幅處理，結合已有鉆孔數據和生產揭露煤厚數據，對預測數據進行驗證和修訂。