林建東，魏名地，孟凡彬

(1.中國(guó)煤炭地質(zhì)總局地球物理勘探研究院，河北 涿州 072750；2.中國(guó)礦業(yè)大學(xué)資源與地球科學(xué)學(xué)院，江蘇 徐州 221008)

高密度三維地震勘探譜矩法反演煤厚技術(shù)

林建東1，魏名地2，孟凡彬1

(1.中國(guó)煤炭地質(zhì)總局地球物理勘探研究院，河北 涿州 072750；2.中國(guó)礦業(yè)大學(xué)資源與地球科學(xué)學(xué)院，江蘇 徐州 221008)

結(jié)合研究區(qū)內(nèi)地質(zhì)資料，建立4個(gè)譜矩法反演煤層厚度的測(cè)試模型，模型參數(shù)分別為子波頻率、煤層頂?shù)装鍘r性、煤層夾矸及積分區(qū)間。不同模型的反演結(jié)果表明，地震子波主頻大小、煤層頂?shù)装鍘r性的差異、以及煤層中夾矸的變化基本上不影響譜矩法反演煤層厚度，但不同的積分區(qū)間對(duì)譜矩法反演煤厚影響較大。以子波主頻50 Hz、積分區(qū)25～50Hz對(duì)研究區(qū)內(nèi)8煤層的厚度進(jìn)行譜矩法反演預(yù)測(cè)，結(jié)果顯示，研究區(qū)西南部煤層厚度變化較大，最大厚度達(dá)14m；中西部煤層較厚，變化較緩；西北部煤層厚度較薄，一般小于2m。對(duì)比25個(gè)鉆孔數(shù)據(jù)，三維地震數(shù)據(jù)譜矩法反演的8煤層厚度，有17個(gè)鉆孔相對(duì)誤差小于18%。

譜矩法；煤層厚度；正演模型；高密度三維地震

0 引言

幾十年來(lái)，煤炭地質(zhì)科技工作者致力于對(duì)煤層厚度變化的研究，但其研究成果受限于所獲地質(zhì)資料，僅能反映出煤層厚度變化的宏觀趨勢(shì)，其精度往往難以滿足礦區(qū)安全高效生產(chǎn)的要求[1-3]。究其根本原因在于煤層賦存狀態(tài)極為復(fù)雜，難于利用有限的地質(zhì)資料精細(xì)分析采區(qū)中復(fù)雜的煤層厚度變化。常規(guī)采區(qū)三維地震受地震勘探方法的局限性，包括來(lái)自煤層反射波信號(hào)的頻帶寬度、分辨率及信噪比等因素影響，僅僅利用煤層反射波信息的運(yùn)動(dòng)學(xué)特征，很難對(duì)落差大于(1/4～1/8)λ的斷裂構(gòu)造進(jìn)行解釋。地震反演技術(shù)是綜合運(yùn)用地震、測(cè)井、地質(zhì)等資料以揭示地下目標(biāo)層(如煤層等)的空間幾何形態(tài)(包括目標(biāo)層厚度、頂?shù)讟?gòu)造形態(tài)、延伸方向、延伸范圍、尖滅位置等)和目標(biāo)層微觀特征。煤礦勘探區(qū)分布數(shù)量較多的鉆孔與巷道資料，可以將縱向高分辨率的鉆孔資料與橫向穩(wěn)定的高分辨率地震資料結(jié)合，預(yù)測(cè)煤層空間變化及分布[4-12]。

據(jù)此，本文以研究煤層厚度變化導(dǎo)致煤層反射波運(yùn)動(dòng)學(xué)、動(dòng)力學(xué)特征變化的地震響應(yīng)為切入點(diǎn)，利用譜矩法反演研究煤層厚度變化。根據(jù)多種地質(zhì)模型試算，研究、分析譜矩法反演煤層厚度的主要影響因素。利用高密度三維地震勘探資料，基于鉆井約束下定量預(yù)測(cè)研究區(qū)煤層厚度,變化規(guī)律，探索出一種利用三維地震信息精細(xì)判別煤層厚度變化的方法。

1 井約束譜矩法定量預(yù)測(cè)煤厚原理

煤層的構(gòu)造或巖性變化主要反映在密度、速度及其它彈性參量的差異上。這些差異導(dǎo)致了地震波在傳播時(shí)間、振幅、相位、頻率等方面存在異常。當(dāng)煤層厚度變化較大時(shí)，會(huì)引起介質(zhì)的彈性物理力學(xué)性質(zhì)的變化，這些變化表現(xiàn)在時(shí)間剖面上為反射與反射波組的振幅、相位和頻率發(fā)生較大變化。通過建立地震屬性與煤層厚度之間的統(tǒng)計(jì)關(guān)系，就可以直接利用地震屬性參數(shù)來(lái)預(yù)測(cè)煤層厚度的變化。

煤層是典型的低速薄層，煤層反射波主要是煤層頂?shù)装宸瓷浏B加的復(fù)合波，其波形具有入射地震子波的微分形式。在一定的條件下，煤層厚度近似正比于煤層反射波振幅譜積分與地震子波振幅譜一階矩的比值。根據(jù)一定數(shù)量的鉆井煤厚資料，標(biāo)定其比例系數(shù)，便可直接從地震反射資料逐道定量估算煤層厚度。

在無(wú)限介質(zhì)中有一低速煤層，如圖1，設(shè)煤層厚度為H、縱波速度為ν、煤層的頂?shù)装宸瓷湎禂?shù)絕對(duì)值為R，τ為煤層時(shí)間厚度，那么縱波在煤層中雙程傳播時(shí)間為τ=2H/v。

圖1 薄煤層模型示意圖Figure 1 A schematic diagram of thin coal seam model

當(dāng)震源脈沖bx,t(地震子波)入射到煤層時(shí)，則煤層反射波sx,t為

s(x,t)=h(t)*b(x,t)

(1)

式中ht—煤層的脈沖響應(yīng)。

對(duì)于薄煤層來(lái)說(shuō)，在其一階近似條件下，式(1)變?yōu)?/p>

(2)

對(duì)式(2)兩邊取傅氏變換、積分，得

(3)

那么式(3)變?yōu)?/p>

(4)

式(4)為譜矩法地震資料反演煤層厚度的基本公式。由于D為地震子波的譜一階矩，故稱此法為譜矩法。

地震子波bx,t的求取可依據(jù)理論估算、地震道自相關(guān)、同態(tài)系統(tǒng)估算等方法，系數(shù)C由鉆孔附近地震道數(shù)據(jù)來(lái)標(biāo)定。假設(shè)勘探區(qū)內(nèi)有N個(gè)鉆孔，第i個(gè)鉆孔處煤層厚度為Hi，求該鉆孔處煤層反射波頻譜積分Ai，地震子波的一階頻譜Di，那么根據(jù)式(4)求得第i個(gè)鉆孔處的標(biāo)定系數(shù)Ci為

對(duì)于離鉆孔較遠(yuǎn)處的標(biāo)定系數(shù)由鉆孔處的Ci(i=1,2,…,N)通過內(nèi)插方法求得。

2 井約束譜矩法參數(shù)選擇

從地震資料中預(yù)測(cè)煤層厚度是根據(jù)一個(gè)無(wú)限均勻高速介質(zhì)中間夾有一煤層的簡(jiǎn)單模型而建立的。根據(jù)鉆孔資料與地震解釋的結(jié)果，依據(jù)研究區(qū)8煤層實(shí)際地質(zhì)情況，給定模型參數(shù)如表1所示。其中，表中W1、W2分別代表頂、底板巖性；C代表煤層；b代表夾矸，本研究區(qū)夾矸主要是炭質(zhì)泥巖。煤層設(shè)為構(gòu)造煤，煤層厚度為0～20m，煤層中含夾矸時(shí)，煤層厚度不變，為20m，夾矸厚度為0～0.8m。

表1 模型物性參數(shù)

建立譜矩法參數(shù)測(cè)試模型如圖2所示。該地質(zhì)模型考慮多種因素，如子波頻率不同、煤層頂?shù)装鍘r性不對(duì)稱、煤層中含有夾矸等參數(shù)對(duì)譜矩法反演的影響。

模型1(圖2a)主要研究地震子波主頻對(duì)與譜矩法反演的影響，分別選用主頻為50、60、70、80Hz及90Hz的地震子波輸入；模型2(圖2b)主要研究積分區(qū)間對(duì)譜矩法反演煤厚的影響，選擇的積分區(qū)間為25～50、25～75、25～100Hz以及50～100Hz；模型3(圖2c)與模型1、模型2以相同的參數(shù)進(jìn)行譜矩法反演，研究頂、底板巖性對(duì)譜矩法反演的影響；模型4(圖2d)主要研究煤層中含夾矸時(shí)對(duì)譜矩法反演煤層厚度是否有影響。

圖2 譜矩法反演地質(zhì)模型Figure 2 Spectral moment inversion geological model

正演模擬采用聲波波動(dòng)方程法進(jìn)行計(jì)算。采樣間隔為0.5ms，體積元剖分大小為0.5m×0.5m，接收道距為5m，共201道。

2.1 子波主頻對(duì)煤厚反演的影響

針對(duì)模型1使用相同的積分區(qū)間(25～50Hz)，地震子波的主頻分別選取50、60、70、80Hz及90Hz，分析不同主頻對(duì)譜矩法反演的影響，如圖3所示。當(dāng)子波主頻不同時(shí)，譜矩法反演的煤厚曲線結(jié)果基本重合，子波主頻的大小基本不影響譜矩法反演結(jié)

果。結(jié)合研究區(qū)地震主頻大小， 本課題選用主頻為50Hz的雷克子波進(jìn)行譜矩法反演煤層厚度。

2.2 積分區(qū)間對(duì)煤厚反演的影響

以主頻為50Hz模型2作為研究對(duì)象，分別選取25～50、25～75、25～100、50-100Hz的4個(gè)不同積分區(qū)間進(jìn)行反演，結(jié)果如圖4所示。由于積分區(qū)間不同，反演的結(jié)果存在明顯差異。當(dāng)煤厚較小(小于5m)時(shí)，不同積分區(qū)間相對(duì)誤差較小。煤層厚度越大，積分區(qū)間的影響越明顯，積分區(qū)間相對(duì)誤差隨著煤厚的增大而增大。積分區(qū)間向高頻移動(dòng)，差異也越大；隨著煤厚增加，反演偏離線性的誤差越大；當(dāng)積分區(qū)間過大(25～100Hz、50～100Hz) 時(shí)，積分頻率區(qū)間影響曲線規(guī)律出現(xiàn)偏差，不適合用來(lái)預(yù)測(cè)煤厚值。經(jīng)驗(yàn)證，積分區(qū)間25～50Hz與25～75Hz均可用于譜矩法反演，而積分區(qū)間25～50Hz比積分區(qū)間25～75Hz的變化范圍大，更適合來(lái)預(yù)測(cè)煤厚。因而本課題選用積分區(qū)間 25～50Hz進(jìn)行譜矩法反演煤層厚度分析與研究。

圖3 不同頻率子波譜矩法反演圖Figure 3 Spectral moment inversion for different frequency wavelets

圖4 不同積分頻率區(qū)間譜矩法反演圖Figure 4 Spectral moment inversion for different integral intervals

2.3 頂、底板巖性對(duì)煤厚識(shí)別的影響

分別對(duì)模型1、2、3進(jìn)行正演模擬，主頻選用50Hz的Ricker子波，積分區(qū)間選用25～50Hz，其結(jié)果如圖5所示?？梢姰?dāng)頂?shù)装鍘r石物性不對(duì)稱時(shí)，反演的煤厚與對(duì)稱的情況基本一致，尤其是模型1即頂?shù)装鍖?duì)稱( 均為W1) 和模型1即頂( W1) 底板( W2) 不對(duì)稱時(shí)，兩者很接近；當(dāng)煤厚大于9m時(shí)，前兩種情況與模型3即頂( W2) 底板( W1) 不對(duì)稱(高速圍巖做頂板) 有較明顯差距。

2.4 夾矸對(duì)煤厚識(shí)別的影響

對(duì)于煤層來(lái)說(shuō)，由于沉積過程中環(huán)境的變化，在煤層中存在夾矸，采用譜矩法對(duì)于含有夾矸煤層來(lái)說(shuō)，對(duì)解釋它的厚度的結(jié)果是否有影響需要研究。模型1、2、3即圖2a、圖2b、圖2c夾矸為零，煤層由薄到厚逐漸增加，模型4(圖2d)中夾矸由薄變厚，

煤層厚度不變。利用譜矩法反演煤厚進(jìn)行預(yù)測(cè)對(duì)比分析，結(jié)果表明：夾矸的存在基本不影響煤層厚度變化的預(yù)測(cè)。

多種地質(zhì)模型試算的結(jié)果表明，地震子波主頻、煤層頂?shù)装鍘r性的對(duì)稱性、以及煤層中夾矸的存在基本上不影響譜矩法反演煤層厚度；積分區(qū)間不同，對(duì)譜矩法反演煤厚的結(jié)果差異較大。經(jīng)過驗(yàn)證，本課題選用積分頻率區(qū)間25～50Hz進(jìn)行研究區(qū)煤層厚度譜矩法反演預(yù)測(cè)。

3 研究區(qū)8煤層地質(zhì)概況

依據(jù)研究區(qū)內(nèi)的地質(zhì)、測(cè)井等資料可知，受區(qū)域構(gòu)造控制，區(qū)內(nèi)煤層揉皺變形產(chǎn)生層滑或斷裂式層滑現(xiàn)象極為嚴(yán)重，因而導(dǎo)致煤層厚度發(fā)育極不穩(wěn)定，并在小區(qū)域內(nèi)煤層厚度產(chǎn)生劇烈變化。

8煤組含煤1～2層，一般含煤1層(即8-1煤層與8-2煤層合并)，8-2煤層為本區(qū)的主要可采煤層，區(qū)內(nèi)發(fā)育良好，煤層厚度較大。8煤組頂部為一厚層砂巖，中、下部圍巖多以泥巖或粉砂巖為主。

圖5 不同頂、底板巖性譜矩法反演圖Figure 5 Spectral moment inversion for different roof and floor lithology

①主采煤層8煤層厚度一般較大，變化也較大，個(gè)別點(diǎn)突然變薄，可能受小構(gòu)造影響。

②區(qū)內(nèi)中淺部形成厚煤帶，煤厚為5～8m，局部有>10m的特厚點(diǎn)，但向深部煤層厚度逐漸變薄。

③在非合并區(qū)煤厚相對(duì)較薄，為3～5m，厚度變化相對(duì)也較小。

④煤層結(jié)構(gòu)簡(jiǎn)單～較簡(jiǎn)單，部分含1層夾矸，個(gè)別含2～3層夾矸，偶見4層夾矸；夾矸多為泥巖或炭質(zhì)泥巖。

⑤受巖漿巖侵蝕影響，部分孔可見侵蝕變質(zhì)天然焦，形成一小范圍的不可采區(qū)。煤類除受巖漿巖侵蝕影響外，變化不大，為全區(qū)可采的較穩(wěn)定煤層。

研究區(qū)典型地震剖面如圖6所示。從圖中可以看出，區(qū)內(nèi)淺層和深層地震地質(zhì)條件較好。

圖6 青東礦研究區(qū)典型地震剖面Figure 6 Qingdong coalmine study area typical seismic section

通過對(duì)研究區(qū)地質(zhì)資料的搜集、分析與研究，結(jié)合實(shí)際巷道揭露，確定本區(qū)煤層厚度變化的主要因素為后生變化，層滑構(gòu)造是煤層厚度后生變化的主要影響因素之一。因此本研究以地震煤厚反演為基礎(chǔ)，根據(jù)延伸小、走向變化快的斷層及構(gòu)造煤的分布輔助驗(yàn)證煤層厚的變化。

4 8煤層厚度變化與精度分析

4.1 高密度三維地震資料反演成果分析

本次課題選取積分頻率區(qū)間25～50Hz進(jìn)行譜矩法反演煤層厚度應(yīng)用。首先利用高密度三維地震數(shù)據(jù)對(duì)研究區(qū)進(jìn)行譜矩法反演8煤層厚度。區(qū)內(nèi)共有鉆孔25口，利用其附近地震道數(shù)據(jù)得到標(biāo)定系數(shù)，進(jìn)而在研究區(qū)內(nèi)插值，聯(lián)合反射波頻譜積分與地震子波頻譜一階矩反演研究區(qū)煤層厚度(圖7)。

從圖7中可以看出，研究區(qū)內(nèi)煤層厚度整體變化較大且變化復(fù)雜，與本研究區(qū)8煤層受層滑構(gòu)造影響煤厚變化情況吻合。研究區(qū)區(qū)西南部8煤層厚度變化較大，中西部煤層變化較緩，煤厚較穩(wěn)定，西北部煤層厚度較薄。在研究區(qū)西北部煤層厚度一般小于2m，鉆孔7-8-8處8煤層厚度僅為1.04m；西南部、南部部分煤層厚度較大，最大厚度為14m，鉆孔8-d-1處8煤層厚度為11.91m。

4.2 譜矩法反演煤厚精度分析與對(duì)比

研究區(qū)內(nèi)833機(jī)巷為已揭露巷道，根據(jù)對(duì)833機(jī)巷的實(shí)際資料分析，833機(jī)巷中、右部層滑構(gòu)造較為發(fā)育，具體表現(xiàn)為煤層厚度變薄，平均煤層厚度小于2m(圖7)?？梢钥闯?，833機(jī)巷穿過煤層厚度變薄區(qū)，煤厚變薄區(qū)的平均煤厚小于2m。

833機(jī)巷穿過井7-x-2，井7-x-2處8煤層厚度為4.17m，利用高密度三維地震數(shù)據(jù)譜矩法反演得到的井7-x-2處8煤層厚度為3.91m，絕對(duì)誤差僅為-0.26m，相對(duì)誤差為6.24%。

利用研究區(qū)內(nèi)鉆孔信息，驗(yàn)證高密度三維地震數(shù)據(jù)與常規(guī)數(shù)據(jù)譜矩法反演煤層厚度的精度。提取研究區(qū)25個(gè)鉆孔的8煤厚度與相應(yīng)位置三維地震數(shù)據(jù)譜矩法反演的8煤層厚度進(jìn)行對(duì)比分析(圖8)。從圖中可以看出： 基于三維地震數(shù)據(jù)譜矩法反演的煤層厚度精度較高，基本滿足煤層厚度預(yù)測(cè)精度。其中，25個(gè)鉆孔中有17個(gè)鉆孔， 其高密度三維地震數(shù)據(jù)譜矩法反演的煤層厚度預(yù)測(cè)精度較高，表明該區(qū)域基于高密度三維地震數(shù)據(jù)譜矩法反演煤層厚度的精度較高。

提取不同三維地震數(shù)據(jù)譜矩法反演煤層厚度的成果， 建立不同三維地震數(shù)據(jù)譜矩法反演煤層厚度的相對(duì)誤差對(duì)比圖，如圖9，研究不同三維地震數(shù)據(jù)對(duì)譜矩法反演煤層厚度精度的影響。從圖中可以看出：不同三維地震數(shù)據(jù)利用譜矩法反演煤層厚度的精度差異較大，研究區(qū)內(nèi)共有鉆孔25個(gè)，其中17個(gè)鉆孔處高密度三維地震數(shù)據(jù)譜矩法反演煤層厚度精度較高，相對(duì)誤差小于18%。存在個(gè)別鉆孔反演煤厚相對(duì)誤差較大，分析發(fā)現(xiàn)，可能是由于該井處于采區(qū)邊界。由此看出，基于不同三維地震數(shù)據(jù)利用譜矩法反演煤層厚度的精度不同，高密度三維地震數(shù)據(jù)對(duì)譜矩法反演煤層厚度的精度更高，因此，在研究區(qū)內(nèi)首選高密度三維地震數(shù)據(jù)進(jìn)行反演煤層厚度的計(jì)算。

圖7 研究區(qū)巷道驗(yàn)證高密度數(shù)據(jù)譜矩法反演8煤層厚度精度圖Figure 7 Precision chart of coal No.8 thickness high density 3D seismic data spectral moment inversion identified by study area roadways

圖8 研究區(qū)鉆孔處8煤厚度與譜矩法反演煤層厚度對(duì)比圖Figure 8 Comparison of study area coal No.8 thicknesses from boreholes and high density 3D seismic data spectral moment inversion

圖9 三維地震數(shù)據(jù)體譜矩法反演煤厚相對(duì)誤差對(duì)比圖Figure 9 Comparison of coal thickness fractional errors from 3D seismic data volume and spectral moment inversion

綜上所述，針對(duì)本研究區(qū)而言，不同三維地震數(shù)據(jù)對(duì)譜矩法反演煤層厚度的精度有影響，基于高密度三維地震數(shù)據(jù)利用譜矩法反演煤層厚度的精度高于基于常規(guī)三維地震數(shù)據(jù)利用譜矩法反演煤層厚度的精度，因此針對(duì)本研究區(qū)，本課題推薦使用高密度三維地震數(shù)據(jù)進(jìn)行反演、預(yù)測(cè)煤層厚度。

5 結(jié)論

在充分搜集、整理研究區(qū)的地質(zhì)、測(cè)井資料的基礎(chǔ)上，對(duì)研究區(qū)煤層厚度變化規(guī)律進(jìn)行了精細(xì)的地質(zhì)分析與研究。結(jié)合研究區(qū)地質(zhì)資料建立正演模型，研究、分析影響譜矩法反演煤層厚度的主要因素，利用高密度三維地震勘探資料，基于鉆井約束下譜矩法編寫應(yīng)用程序，定量預(yù)測(cè)研究區(qū)煤層厚度變化。豐富了復(fù)雜煤層厚度變化反演方法，其適用性更為廣泛，預(yù)測(cè)精度較以往鉆孔對(duì)比、插值等煤厚反演方法有明顯提高，為煤礦安全生產(chǎn)提供了保障。

[1] 劉程, 李向東, 楊守國(guó).地質(zhì)構(gòu)造對(duì)煤層厚度的影響研究[J].煤礦安全, 2008, 39(5):14-16.

[2] 張德亮.煤層厚度變化的地質(zhì)成因分析[J].內(nèi)蒙古煤炭經(jīng)濟(jì), 2015(12):211-212.

[3] 李剛.基于BP算法的煤層厚度預(yù)測(cè)技術(shù)應(yīng)用研究[J].中國(guó)煤炭地質(zhì), 2011, 23(5):45-48.

[4] 吳秀英, 胡海平.基于小波變換和人眼視覺系統(tǒng)特性的圖像降噪方法[J].應(yīng)用數(shù)學(xué)與計(jì)算數(shù)學(xué)學(xué)報(bào), 2011, 25(1):74-79.

[5] 黃卉.基于小波變換的圖像融合方法研究[D].合肥：合肥工業(yè)大學(xué), 2005.

[6] 陳曉東, 朱俊杰, 郭華東,等.基于小波變換和局部相關(guān)系數(shù)改進(jìn)IHS變換的圖像融合方法[J].地理與地理信息科學(xué), 2005, 21(6):22-24.

[7] 王建青, 楊曉東, 田忠斌,等.基于地震頻譜屬性的煤層厚度預(yù)測(cè)方法[J].物探與化探, 2017,41(2): 364-368.

[8] 索重輝,常鎖亮,彭仕宓,等.地震屬性在煤層厚度預(yù)測(cè)上的分析與應(yīng)用[J]. 科學(xué)技術(shù)與工程,2011,11(34): 8429-8433.

[9] 潘冬明,李娟娟,胡明順,等.多煤層的反射地震響應(yīng)特征研究[J].地球物理學(xué)進(jìn)展,2010,25(5):1769-1744.

[10] 王建軍.優(yōu)選多地震屬性線性回歸方法在預(yù)測(cè)煤層厚度中的應(yīng)用[J].工程地球物理學(xué)報(bào), 2014, 11(5):609-613.

[11]馮李治欣, 鄧國(guó)成.地震多屬性技術(shù)解釋煤礦薄煤區(qū)的應(yīng)用[J].科技創(chuàng)新與應(yīng)用, 2016(19):27-28.

[12] 陳優(yōu)闊.基于核方法的煤層厚度變化預(yù)測(cè)模型及應(yīng)用研究[D].江蘇徐州：中國(guó)礦業(yè)大學(xué), 2016.

CoalThicknessHighDensity3DSeismicProspectingSpectralMomentInversionTechnology

Lin Jiandong1, Wei Mingdi2， Meng Fanbin1

(1.Geophysical Prospecting and Research Institute, CNACG, Zhuozhou, Hebei 072750；2.College of Resource and Earth Science, China University of Mining and Technology, Jiangshu，Xuzhou 221008)

Based on geological data from study area, four testing models of coal thickness from spectral moment inversion have been modeled. Model parameters have wavelet frequency, lithology of coal roof and floor, coal parting and integral interval respectively. The inversion results from different models have shown that the seismic wavelet master frequency range, difference in coal roof and floor lithology and parting variation in coal seam are basically no impacts on coal thickness spectral moment inversion, but impacts from integral interval are large. Taking the wavelet master frequency 50Hz, integral interval 25～50Hz carried out coal seam No.8 thickness spectral moment inversion prediction. The result has shown that coal seams in the southwestern study area have larger variations, maximum thickness can be 14m as high; in middle western part rather thick and gently varied; in northwestern part rather thin, generally less than 2m. Comparison of data from 25 boreholes has found that the thickness fractional errors of 17 boreholes are less than 18% from coal No.8 3D seismic data spectral moment inversion.

spectral moment method; coal thickness; forward model; high density 3D seismic prospecting

國(guó)家“十三五”科技重大專項(xiàng)經(jīng)費(fèi)資助項(xiàng)目(2016ZX05041002)、(2016ZX05043001),淮北礦業(yè)(集團(tuán))有限責(zé)任公司科研及技術(shù)創(chuàng)新計(jì)劃項(xiàng)目(青東礦8煤厚度變化規(guī)律研究及應(yīng)用)。

林建東(1967—)，男,高級(jí)工程師,主要從事煤田地震勘探生產(chǎn)與科研。

2017-08-01

10.3969/j.issn.1674-1803.2017.09.13

1674-1803(2017)09-0067-07

A

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