地震屬性參數在煤層厚度預測中的應用對比

馮躍封

（河南省煤田地質局 地質二隊，河南 洛陽 471023）

0 前言

煤礦井下瓦斯抽排、礦注漿堵水和大型、超大型煤礦井建設，對煤田地質勘查工作提出了許多新的更高的要求。煤田高分辨率三維地震勘探，通過地下空間三維數據體數據的采集和處理，不僅可以探明采區內的細微構造，還可提供下賦煤層厚度變化的相關信息［1］。

通常情況下，采區煤層厚度是由鉆孔資料對比、內插的方法獲得的［2］，但在鉆孔稀少或鉆孔煤層資料缺失的地區，僅靠少數鉆孔資料來確定采取煤炭儲量說服力就明顯不足。因此，確定煤礦采區內煤層厚度，就顯得尤為重要。通過對三維地震勘探中采集到的煤層參數的分析，來判斷采區煤層厚度的變化趨勢，這為煤礦礦井設計提供了科學依據［3］。

1 煤層厚度研究現狀

根據幾何地震學原理，利用時距曲線將某個層或段的雙程旅行時間，轉換為厚度數據來判斷大套地層的厚度，長期以來一直在煤田地震勘探中普遍使用［4］，但對厚度小于10m （一般為1.0m～6.0m）的煤層來說，由于在頂底板位置形成的地震反射波表現為復合波，從而為煤層厚度的判斷造成困難。

為解決薄煤層地區通過三維地震采集數據分析來判斷煤層厚度的問題，不少學者從理論上討論了煤層反射波的形成機制，提出了多種煤層厚度的解釋方法，利用地震波的波形、振幅與頻率的的變化來研究判斷煤層厚度的變化，為通過煤層地震屬性參數變化進行厚度定量的預測，提供了理論依據［5］。利用地震屬性參數研究煤層厚度的方法主要包括振幅法、直接反演法和統計分析法三種。特別是隨著計算機技術的普遍應用，使薄層煤厚度預測成為可能，較好地滿足了煤礦生產對煤層厚度精確判斷的要求。

2 煤層厚度地震屬性預測

煤層構造或巖性的變化，導致了地震波在煤層中的傳播時間、振幅、相位、頻率等方面發生變化［6］。在地震剖面上可以清楚地看到：厚度變化較大的煤層，地震波同相軸的走勢、振幅和相位變化明顯，但頻率域的變化卻難以直觀地分析。因此，用常規人工解釋方法來判斷煤層厚度的變化，效果很不理想。通過對煤層厚度變化引起的地震屬性變化規律的研究，建立煤層地震參數變化與煤層厚度之間的統計關系，就可以直接利用地震屬性參數來預測煤層厚度的變化［7］。

由于煤層厚度變化會引起煤層物理力學性質的變化，而這些變化必然引起地震時間剖面上反射波組的振幅、相位和頻率的變化。因此，可以根據煤層這些物性的變化來提取煤層的時間域反射波能量、最大振幅值，頻率域的主頻、頻寬等地震屬性參數，從而對煤層的厚度進行預測。

2.1 振幅譜平方比法

設已知厚度與未知厚度的煤層反射波振幅譜分別為A0（f）、A（f），可得：

式中 A（f）為反射波振幅譜；B（f）為地震子波振幅頻譜；ω為入射波圓頻率；K1、K2為反射系數。

令：

將式（1）代入式（2）得式（3）。

式中 P為反射系數參量，且

根據式（3），應用反余弦函數可得到第一周期內的τ值：

由τ即可求得煤層厚度h＝Vτ／2（因為0≤ωτ≤π，所以式（5）的應用條件是f≤1／2τ，即h≤λ／4）。

2.2 振幅譜比法

根據簡化公式（1）得：

子波的振幅譜可通過反射記錄的自相關求取，根據地震道可得：

令：C＝4πk／（1－K2）υ2。

由式（7）得式（8）。

若某地震道煤層厚度已知（如過井道），由式（8）可得式（9）。

由此可得該地震道的煤層厚度為式（10）。

2.3 譜矩法

設某采區煤層厚度為H、縱波的速度為V、頂底板反射系數絕對值為R，則縱波在煤層中雙程傳播的時間t＝2 H／V。當地震子波b（x，t）入射到煤層時，煤層反射波s（x，t）為式（11）。

式中 h（t）為煤層的脈沖響應。對于薄煤層，在近似一階條件下，式（11）可變為式（12）。

地震子波b（x，t）可通過理論估算、地震道自相關和同態系統估算等方法求取。

對式（12）兩邊取傅氏變換、積分，得式（13）。

通過鉆孔附近地震道數可標定參數C，根據勘探區內第i個鉆孔附近煤層厚度，可求得煤層反射波頻譜積分At與地震子波的一階頻譜Dt。

3 煤層厚度預測方法對比

根據上述理論和方法，我們編制了相應的程序，以方便廣大技術人員使用。該系統采用當今流行的VB開發環境，運行于Windows操作系統下，主要運用了地震子波的自相關法、對數統計法，綜合分析三維地震中獲得的地質信息對煤層厚度進行預測。下面以山東某煤礦三維地震資料為例進行分析。

3.1 不同預測方法

依據鉆孔標定煤層厚度數據，分別采用振幅譜平方比法、振幅譜比法和譜矩法進行對比，結果如表1所示。從表1可以看出，采用譜矩法，預測出的煤層厚度與實際煤層厚度的絕對誤差和相對誤差浮動最小，算法最穩定。利用開發的軟件，求取煤層的厚度曲線如圖1所示。

由以上分析可知，譜矩法和振幅譜比法對煤層厚度的預測結果較為一致，譜矩法預測結果與鉆孔標定煤層厚度更加接近。

3.2 預測結果驗證

從山東某煤礦三維地震資料中提取地震子波，采用譜矩法對礦區煤層厚度進行預測，選擇25Hz～50Hz的頻率區間，逐道計算煤層厚度，對照采區內鉆孔煤層資料，預測結果與揭露煤層厚度基本吻合（見表2）。

表1 煤層厚度預測方法對比Tab.1 The comparison of the prediction method of coal seam thickness

圖1 不同方法預測煤層厚度曲線（121線）Fig.1 Curves of coal seam thickness predicted by different methods

表2 煤層實際厚度與預測厚度對比Tab.2 The comparison of coal seam thickness predicted by different methods

4 結論

（1）譜矩法預測煤層厚度精度：①對厚度在4m～10m的煤層預測效果好；②隨著煤層厚度的增加，預測誤差相對越大；③譜矩法預測煤厚，應多標定常數C值；④振幅譜預測煤層厚度應經過濾波設置；⑤振幅譜平方比法預測煤層厚度，預測范圍較大，但精度有限。

（2）作者開發的煤層厚度預測軟件，利用地震屬性參數，借助計算機對煤層厚度進行測算，精度高，操作簡單。

［1］蔡軍.譜矩法預測煤層厚度的應用［J］.內江科技，2008（11）：124.

［2］董守華，許永忠.地震資料譜矩法反演煤層厚度［J］.遼寧工程技術大學學報，2005，24（1）：38－40.

［3］張恩珂，柏冠軍，鄭中剛.地震資料頻率域預測煤層厚度方法研究及應用［J］.山東科技大學學報：自然科學版，2005（12）：59－62.

［4］徐衍和，李東會.煤田地震資料橫向預測技術［J］.勘察科學技術，2006（6）：35－38.

［5］張玉峰，夏玲燕.3D地震資料煤層厚度變化規律的研究［J］.山西師范大學學報：自然科學版，2008（03）：89－91.

［6］蔣賢東，柏冠軍.地震資料定量預測煤層厚度譜矩法及其應用［J］.能源技術與管理，2005（6）：26－28.

［7］董守華，馬彥良，周明.煤層厚度與振幅、頻率地震屬性的正演模擬［J］.中國礦業大學學報，2004（01）：29－32.