地震資料解釋巖漿巖侵入煤層范圍的探索研究

湯紅偉

(中煤科工集團西安研究院有限公司，陜西省西安市，710077)

地震資料解釋巖漿巖侵入煤層范圍的探索研究

湯紅偉

(中煤科工集團西安研究院有限公司，陜西省西安市，710077)

巖漿侵入煤層對煤礦安全高效生產有較大影響。在有鉆孔揭露巖漿巖區域，應用波阻抗反演同時輔以地震屬性分析等三維地震勘探技術預測巖漿巖侵入煤層現象能取得較好的效果。對于無井區域來說，反演方法不能運用，單純的屬性分析技術存在多解性。正演模擬表明，巖漿巖侵入煤層后會使得反射波主相位與輔助相位存在較大的時間差及振幅差，此差異為利用巖漿巖地震響應預測巖漿巖侵入煤層范圍提供了可能。以陳四樓煤礦八采區二2煤層為例，利用地震時間剖面、屬性分析相結合的方法，圈定了巖漿侵入煤層范圍。通過巷探實際驗證，解釋結果吻合較好。

巖漿巖侵入 時間剖面 地震屬性 地震資料 三維地震勘探

AbstractMagma intruding into coal seam had a greater impact on coal mine efficient and safety production. In the magmatic rock area exposed by drillhole, 3D seismic exploration techniques such as wave impedance inversion and seismic attribute analysis, had preferable effects on the prediction of magma intruding into coal seam. For the non-well area, the retrieval method couldn't be applied well, and simple attribute analysis had multiple solutions. The forward modelling showed that after magmatic rocks intruding into coal seam, there were big time difference and amplitude difference between master phase and assisted phase of reflected wave, and the differences made it possible for magmatic rock seismic response to predict the range of magmatic rock intrusion. Taking two-2 coal seam of No. 8 panel of Chensilou Coal Mine as an example, the range of magmatic rock intrusion in coal seam was determined by using combined method of seismic time profile and attribute analysis. The actual roadway protection results was identical with the interpretation results.

Keywordsmagmatic rock intrusion, time profile, seismic attribute, seismic data, 3D seismic exploration

巖漿巖的侵入對煤層的影響很大，一方面，巖漿巖侵入煤層后使得煤層的完整性受到了破壞，而且由于巖漿的烘烤，使與巖漿巖接觸部分的煤質變差甚至變成天然焦，影響了煤的工業價值，降低了可采煤炭資源儲量。另一方面，由于巖漿巖的硬度較大，掘進速度受到影響，而且使采掘設備的使用壽命減少，從而增加了開采成本。因此，準確圈定巖漿巖侵入煤層的范圍對煤礦安全高效生產具有重大的現實意義。

在地質勘查階段，多應用鉆孔資料圈定巖漿巖侵入煤層的范圍，但受鉆孔數量及鉆孔密度的限制，圈定邊界的精度較差。相比之下，三維地震勘探資料具有大面積密集采集及橫向分辨率高的優勢，因此，波阻抗反演、神經網絡及灰色關聯分析、地震相分析及頻譜分解、地震屬性技術等三維地震勘探技術預測巖漿巖侵入煤層的方法在不同地區得到應用。波阻抗反演是利用反射波阻抗值的大小來區分巖漿巖侵入煤層范圍的，巖漿巖一般表現為高阻抗值、煤層為低阻抗值；屬性分析技術主要依據反射波連續性、反射波主頻的高低來識別巖漿巖侵入區，在時間剖面上，正常煤層區域反射波能量較強、連續性好，天然焦區域的反射波有連續性變差或者反射波消失、能量變弱、頻率增高的地震響應特點。這些方法在鉆孔揭露巖漿巖侵入煤層的信息較豐富、測井資料齊全且巖漿巖侵入的煤層為主要可采煤層的區域都取得了較好的效果。

永城礦區為巖漿巖活動頻繁區域，根據收集的陳四樓井田的237個鉆孔資料統計來看，有59%的鉆孔穿見巖體，侵入的主要為局部可采的三煤組。揭露巖漿巖侵入二2主采煤層的鉆孔只有5個，且鉆孔施工年代久遠，測井資料缺失嚴重，在采掘過程中，經常能發現小范圍的巖漿巖侵入煤層現象，為此，必須采用新的技術方法來預測巖漿巖侵入范圍。

1 概況

陳四樓井田地層由老至新為寒武系下統、中統、上統，奧陶系下統、中統，石炭系中統、上統，二疊系下統、上統，新近系漸新—始新統、中新統、上新統，第四更新統、全新統，含煤地層為石炭系和二疊系。上石炭統太原組含煤1～6層，稱一煤組，煤層總厚2.06 m，含煤系數為1.4%。下二疊統山西組含煤3層，稱二煤組，編號為二1～二3，其中二2煤層為礦井主采煤層，賦存穩定，結構單一，平均厚度2.45 m；下石盒子組含煤7層，稱三煤組，編號為三1～三7。三1煤層厚度為1.30 m，為不穩定煤層；三2煤層厚度為1.50 m，為較穩定煤層；三4煤層厚度為1.60 m，為不穩定煤層；三5煤層平均厚度為0.96 m，為不穩定的局部可采煤層。二煤組與三煤組平均間距為80 m左右。

井田內巖漿活動頻繁，統計井田內237個見煤鉆孔，有59%的鉆孔穿見巖漿巖體。其中巖漿巖侵入三煤組的鉆孔有136個，巖漿巖侵入二2煤層的鉆孔有5個，鉆孔資料揭露巖漿巖最大厚度為31.57 m，平均厚度為2.59 m。巖漿活動有2個期次，其中，輝綠巖偏老，為華力西晚期的產物；花崗斑巖較新，為燕山早—晚期的產物。井田內巖漿巖一般沿三煤組煤層及其底板侵入，產狀多為巖床或巖墻。二2煤層巖漿巖侵入區主要在井田中西部及北部，其中0716孔、6406孔、0908孔、1006孔和1308孔附近發育大片巖漿巖侵入區及天然焦，四采區2408工作面、2409工作面、五采區2513工作面和十二采區21203工作面均受到不同程度的影響。另外，二、七、十一采區局部發育條帶狀巖漿巖墻或天然焦，十采區有巖漿巖順斷層面侵入煤層現象。

2 正演模擬

將原生煤層、天然焦、天然焦與巖漿巖互層、巖漿巖連續存在的現象設計在同一模型中，采用專業正演模擬軟件進行模擬，分析正演模型中反射波的地震響應特征。正演模型如圖1所示，模型中的物理參數見表1，經過處理后得到的正演模擬時間剖面如圖2所示，地震時間剖面中巖漿巖侵入區段和天然焦區段交界會產生明顯的煤層反射同相軸錯斷，同時反射波主相位與輔助相位之間的時間差存在差異(巖漿巖侵入區的主輔相位間的時間差大于天然焦區段)，這是由于這些區段的巖石參數的改變引起橫向上縱波速度變化顯著。焦化區段與正常煤層交界區段反射波同相軸有輕微的振幅變化和波形變化，并未引起目的層反射波同相軸的錯斷和扭曲，為了更好地區分正常煤層與天然焦分布范圍，利用多個不同屬性進行分析。瞬時振幅剖面圖見圖3，從圖中可以看出正常煤層區、天然焦區、巖漿巖侵入區能夠較好地區分。

圖1 放大的正演模型

表1 正演模型物理參數

圖2 60 Hz子波疊前偏移疊加剖面

圖3 60 Hz子波瞬時振幅剖面

3 八采區巖漿巖侵入煤層后反射波屬性特性

3.1 三煤組被巖漿巖侵入后反射波特征

八采區三煤組被巖漿巖嚴重侵蝕，多個鉆孔揭露了巖漿巖體。勘探區聯井剖面圖見圖4，方框內為巖漿巖侵入區域。由圖4可以看出，1004鉆孔被巖漿巖侵蝕，1019鉆孔揭露的煤層為正常煤層，巖漿巖侵入部分的煤層反射波變弱、頻率降低。在地震屬性響應中，巖漿巖侵入部分與原生煤層直接接觸時，能量差異較大，從屬性平面圖中可以較好地區分巖漿巖侵入煤層范圍。

3.2 二2煤層被巖漿巖侵入后的反射波特征

八采區沒有鉆孔揭露巖漿巖侵入二2煤層的現象，但在2802工作面運輸巷的掘進過程中出現了巖漿巖侵入煤層現象。沿2802工作面運輸巷方向的地震時間剖面如圖5所示，圖中yj5處為揭露巖漿巖的位置，yj5右側部分為正常煤層區域。

由圖5可以看出，二2煤層反射波表現為雙相位，上相位較連續、反射波能量較一般；下相位連續性好、能量較強。層位標定中，下相位為二2煤層構造發育的主要反射波，在主相位上，被巖漿巖侵入部分與正常煤層部分的反射波波形特征基本一致。在巖漿巖侵入煤層部位，煤層反射波的上相位頻率變高、上下相位的時間差較大。

圖4 勘探區聯井剖面圖

圖5 2802工作面運輸巷時間剖面圖

二2煤層反射波屬性平面圖見圖6。

圖6 反射波屬性平面圖

反射波屬性平面圖如圖6所示。從圖6(a)中可以看出，巖漿巖侵入地段與正常煤層區域有一定的差異，其振幅能量減弱，能量值大小基本與斷層帶、資料品質較差區域的能量值一致，但斷層帶呈線條狀，而巖漿巖侵入帶與資料品質較差區域的形狀不規則。資料品質較差區域形成的主要原因是地表障礙物的影響。通過炮檢位置分布圖可以進一步區分巖漿巖侵入體與資料品質變差區域。

由圖6(b)可以看出，巖漿巖侵入地段與正常煤層區域差異較小。

4 地質成果

利用煤層反射波均方根振幅值、炮檢位置關系、野外資料品質分布圖、時間剖面上兩個相位的相對時差值大小綜合圈定了八采區二2煤層中的巖漿巖侵入煤層帶2處，侵入帶2在地震時間剖面中的反應如圖7所示。其中，通過巷探揭露確定的巖漿巖侵入邊界與地震資料解釋邊界吻合，另外一處有待進一步驗證。

圖7 巖漿巖侵入帶2在地震時間剖面上的表現特征

5 結論

在無鉆孔資料揭露巖漿巖侵入煤層范圍的區域，不能有效應用反演方法來確定巖漿巖侵入煤層邊界，但可以結合地震時間剖面上反射波(組)的特征變化、煤層反射波均方根振幅值、煤層反射波能量值大小、單炮記錄品質分布圖來確定巖漿巖侵入煤層范圍，實際揭露驗證效果較好。

[1] 朱書階，陳青峰，何文欣.利用波阻抗反演技術確定巖漿巖侵入煤層范圍[J].煤田地質與勘探，2003(1)

[2] 李仁海，崔若飛，毛欣榮等.利用巖性解釋方法確定巖漿巖侵入煤層范圍[J].地球物理學進展，2008(1)

[3] 錢進，崔若飛，陳同俊等.波阻抗反演預測煤層巖漿巖侵入范圍[J].中國煤炭地質，2010(1)

[4] 劉俊杰,王彥春. 煤田數字高密度地震勘探的應用效果探析[J]. 中國煤炭,2013(1)

[5] 唐漢平.復雜地震地質條件下煤礦采空區三維地震勘探技術[J] .中國煤炭，2013(12)

[6] 介偉.地震屬性技術在圈定巖漿侵入煤層范圍的應用[J].中國煤炭地質，2014(9)

[7] 劉萬金，崔凡，梁平.煤田三維地震勘探采集設計中的偏移孔徑探討[J].中國煤炭，2011(4)

[8] 湯紅偉. 基于探采對比分析的三維地震資料精細解釋方法探討[J]. 中國煤炭，2017(4)

[9] 王遠，崔若飛，秦軻.利用地震相分析方法圈定煤層變焦區[J].中國煤炭地質，2010(12)

(責任編輯 郭東芝)

Exploratoryresearchoninterpretationofmagmaticrockintrusionrangeincoalseambyusingseismicdata

Tang Hongwei

(Xi'an Research Institute, China Coal Technology & Engineering Group, Xi'an, Shaanxi 710077, China)

P631.4

A

湯紅偉. 地震資料解釋巖漿巖侵入煤層范圍的探索研究[J]. 中國煤炭，2017，43(9)：35-38. Tang Hongwei.Exploratory research on interpretation of magmatic rock intrusion range in coal seam by using seismic data[J]. China Coal, 2017, 43(9):35-38.

湯紅偉(1980-)，女，河南信陽人，博士，副研究員，畢業于煤炭科學研究總院西安研究院地球探測與信息技術專業，現在中國煤炭科工集團西安研究院從事地震勘探技術研究工作。