張艷偉，屠世浩，王 沉，張 村，魏帥鋒，宋 啟

(1.中國礦業(yè)大學(xué) 礦業(yè)工程學(xué)院，江蘇 徐州 221116；2.中國礦業(yè)大學(xué) 煤炭資源與安全開采國家重點實驗室，江蘇 徐州 221116)

?

基于趨勢面分析的煤層厚度變化規(guī)律研究

張艷偉1,2，屠世浩1,2，王 沉1,2，張 村1,2，魏帥鋒1,2，宋 啟1,2

(1.中國礦業(yè)大學(xué) 礦業(yè)工程學(xué)院，江蘇 徐州 221116；2.中國礦業(yè)大學(xué) 煤炭資源與安全開采國家重點實驗室，江蘇 徐州 221116)

以金莊煤礦8203工作面為工程背景，通過統(tǒng)計工作面煤層鉆孔資料，運用趨勢面分析方法預(yù)測5(3～5)#煤層厚度變化趨勢，得到煤層厚度趨勢面方程，采用F法檢驗煤層厚度趨勢面擬合度，最后采用工程驗證方法證明了趨勢面分析方法預(yù)測煤層工作面煤厚變化規(guī)律的正確性和可靠性。結(jié)果表明：一次趨勢面方程擬合度41.31%，四次趨勢面方程擬合度74.79%，工作面開采期間鉆探鉆孔數(shù)據(jù)驗證顯示和四次趨勢面方程預(yù)測數(shù)據(jù)最大正負誤差率為-9.5%和9.1%，最小正負誤差率為-1.6%和0.8%，且與實際鉆孔測量煤層厚度吻合；四次趨勢面方程計算得到工作面平均采出率84.14%，頂煤回采率81.71%；實際工作面回采率為85.24%，頂煤回采率80.69%，兩者誤差率為1.29%和1.26%.趨勢面分析方法預(yù)測工作面煤層厚度變化規(guī)律是正確可靠的。

特厚煤層；趨勢面分析；煤層厚度變化；F分布

0 引 言

煤層厚度是煤礦設(shè)計與開采必不可少的數(shù)據(jù)，煤層厚度的變化直接影響采掘部署、造成采掘失調(diào)、降低回采率、增加煤厚的研究與探測工作。準(zhǔn)確地預(yù)測煤層厚度，不僅能給煤礦提供有力的地質(zhì)保障，而且還能給煤礦帶來巨大的經(jīng)濟效益。了解并掌握煤層厚度及其變化規(guī)律問題是合理確定開采方法，正確估算煤炭儲量的重要理論基礎(chǔ)[1-3]。

目前國內(nèi)外學(xué)者對煤層厚度變化規(guī)律研究的研究較多，如王峰[4]對邯鄲地區(qū)云陶擴大勘探區(qū)煤厚進行了相關(guān)研究得到了煤層厚度變化等值線；韓艷玲[5]對核桃峪勘探區(qū)煤層賦存進行了研究并進行了3次趨勢面擬合并得到了煤厚變化的影響因素；郭鵬程[6]對邢臺礦區(qū)煤層厚度變化規(guī)律進行了趨勢面四次擬合研究并得到了煤厚與古地形的關(guān)系。以上學(xué)者從理論分析、現(xiàn)場實測和數(shù)值計算等角度均取得了可喜成就。以上研究成果對金莊煤礦的煤厚變化規(guī)律研究奠定了基礎(chǔ)，但是以上方法中并未就趨勢面方程進行現(xiàn)場實測數(shù)據(jù)驗證，文中就趨勢面方程進行開采期間現(xiàn)場鉆孔和采出率對比2方面進行了實測驗證，證明趨勢面分析預(yù)測煤層厚度變化的正確可靠性。國內(nèi)外針對煤層厚度變化的分析方法主要有地質(zhì)分析法、儀器探測法和數(shù)學(xué)地質(zhì)方法[7-14]。由于小范圍煤層厚度突變，其確切變化位置通過地質(zhì)分析方法研究往往難以查明，儀器探測方法由于野外施工困難，探測費用高而難以應(yīng)用，而通過數(shù)學(xué)地質(zhì)方法能夠全面的分析煤層局部和整體的變化規(guī)律，進而達到研究地質(zhì)特征的空間分布和變化規(guī)律的目的[15-19]。

以金莊煤礦8203工作面為工程背景，采用趨勢面分析方法研究8203首采工作面煤層厚度變化規(guī)律。金莊煤礦主采煤層為5(3～5)#煤層，煤層傾角3°～4°。工作面開采前期現(xiàn)場鉆探100個煤層鉆孔，煤層厚度分布頻度如圖1所示。通過多元線性回歸的方法，用一個多項式曲面來代表地質(zhì)體的空間形態(tài)，通過開采前期鉆孔煤層厚度建立趨勢面方程，并驗證分析趨勢面方程煤厚情況。開采期間鉆孔煤厚和工作面采出率驗證趨勢面方程預(yù)測工作面煤厚變化的正確性和可靠性。

圖1 煤層厚度分布頻度

1 趨勢面分析

1.1 趨勢面方程的建立

1.2 原始數(shù)據(jù)預(yù)處理

由于原始數(shù)據(jù)的坐標(biāo)采用“平面直角坐標(biāo)系”，坐標(biāo)值與煤層厚度等值線值之間會產(chǎn)生巨大的級差，造成正規(guī)方程組病態(tài)，又會使計算失去實際意義，甚至還會在計算時產(chǎn)生溢出，必須對坐標(biāo)值進行量綱統(tǒng)一，先對其最小化，然后再均值化，以達到趨勢分析的要求。

對坐標(biāo)進行最小化變換

(1)

進行均勻化變化

(2)

分別對ui′和vi′乘以系數(shù)C，使坐標(biāo)數(shù)據(jù)滿足趨勢面分析處理要求，計算完成后，再將數(shù)據(jù)轉(zhuǎn)化為原始數(shù)據(jù)格式輸出。

2 趨勢面方程的確立、檢驗及分析

2.1 趨勢面方程確立

根據(jù)工作面開采前期100組鉆孔數(shù)據(jù)進行趨勢面分析，建立模型，運用mathematics軟件對煤層厚度和x,y坐標(biāo)值進行1～4次回歸分析，得出趨勢面方程。各次趨勢面分析的趨勢面方程見表1.

表1 煤層厚度趨勢面方程Tab.1 Trend surface equation of coal seam thickness

2.2 趨勢面方程檢驗

為了判斷得到的煤層厚度趨勢面方程能否準(zhǔn)確反映5(3～5)#煤層整體變化趨勢和局部異常情況，用擬合度C來衡量趨勢面方程的擬合程度。計算公式為

(3)

對所建立的趨勢面方程是否顯著，采用F檢驗法，其統(tǒng)計量F的計算式為

(4)

根據(jù)數(shù)理統(tǒng)計理論，式(4)服從自由度為(p，n-p-1)的F分布，給定置信水平的情況下，查F分布表，得到臨界值F，若當(dāng)F>Fα(p，n-p-1)時，認為自變量與因變量相關(guān)，稱趨勢面方程在置信水平α下顯著；若當(dāng)F≤Fα(p，n-p-1)時，認為自變量與因變量不相關(guān)，稱趨勢面方程在置信水平α下不顯著。趨勢面擬合度及F檢驗見表2.

表2 趨勢面擬合度及F檢驗表Tab.2 Trend surface fitting degree and F test

由表2可知，在煤層厚度變化趨勢面分析中，各次趨勢面方程的F>F0.01且F>F0.05，說明所建立的趨勢面方程均顯著，即所建立的趨勢面方程均可靠。趨勢面方程擬合度隨其次數(shù)的增加而增加，也增加了其反映地質(zhì)信息的可信程度，根據(jù)趨勢面擬合度及F檢驗表，選擇四次趨勢面分析的結(jié)果作為研究依據(jù)。

2.3 煤層厚度變化整體趨勢分析

根據(jù)工作面開采前期100組鉆孔數(shù)據(jù)作出的煤層厚度等值線如圖2所示。鉆孔數(shù)據(jù)中8203工作面煤層厚度變化區(qū)間在12.97～24.82m，研究區(qū)域最小厚度12.97m(k90鉆孔)，最大厚度24.82m(k10鉆孔)。等值線圖中煤層厚度分布在13.0～24.0m之間，煤層總體分布較均勻，局部出現(xiàn)煤層厚度變大，或者變小。

圖2 煤層厚度等值線圖

圖3 工作面開采前期/趨勢面方程鉆孔煤層厚度分布頻度

工作面開采前期鉆孔煤層厚度和趨勢面方程煤層厚度分布頻度(圖3)比較可知：煤層厚度主要分布在14.0～20.0m間，工作面開采前期煤厚在此范圍內(nèi)占76%，趨勢面方程煤厚占96%；煤層厚度分布頻度最大在16.0～17.0m間，兩者占的比例分別為19%和22%；工作面開采前期和趨勢面方程鉆孔煤層厚度分布頻度最小值在小于14.0m(包括14.0m)和大于21.0m(包括21.0m)區(qū)間，兩者占的比例為9%，1%和7%，1%，只有局部區(qū)域煤層在此范圍內(nèi)。

2.4 煤層厚度變化區(qū)域性分析

四次趨勢面分析(C=74.79%)如圖4所示，其結(jié)果顯示，煤層厚度變化趨勢總體不均勻，在工作面推進方向上，煤層厚度先降低，再增高，再降低，再增高，最后降低，工作面推進方向中部煤層厚度最大，厚度為20.0m，大巷附近煤層厚度最小，厚度為14.0m，切眼處煤層厚度為17.0～18.0m，煤層的相對富煤區(qū)總體呈現(xiàn)“Z”形，依據(jù)“Z”形逐步降低，局部呈現(xiàn)下凹；工作面中線四次趨勢面煤層厚度見圖5，由圖可知，煤層厚度整體在15.0～17.0m變化，局部區(qū)域煤層厚度增大，在距離切眼長度為700～1 000m間，煤層厚度急劇減小后急劇變大，煤層厚度起伏變化巨大。

圖4 煤層四次趨勢面圖

圖5 工作面中線四次趨勢面煤層厚度

3 工程驗證

3.1 工作面開采期間鉆孔數(shù)據(jù)驗證

根據(jù)煤層厚度四次趨勢面方程，結(jié)合金莊煤礦8203工作面開采期間鉆孔探測(圖6)數(shù)據(jù)，驗證煤層厚度趨勢面方程的可靠性。

由圖7可知，在工作面開采期間28個鉆孔實際煤厚和四次趨勢面方程煤厚誤差率在-10%～10%，其中最大正負誤差率為-9.5%和9.1%，最小正負誤差率為-1.6%和0.8%；實際鉆孔數(shù)據(jù)煤層厚度最大值為18.03m，在距離切眼長度為1 400m處，最小值為13.76m，在距離切眼長度為100m處；四次趨勢面方程計算煤層厚度最大值為18.37m，在距離切眼長度為1 400m處，最小值為14.27m，在距離切眼長度為100m處。實際鉆孔煤層厚度最大值和最小值與四次趨勢面方程計算煤層厚度值所在位置相同，且誤差率分別為1.9%和3.7%.四次趨勢面方程計算煤層厚度和實際鉆孔測量煤層厚度基本吻合，因此四次趨勢面方程能夠準(zhǔn)確預(yù)測首采面煤層分布情況和規(guī)律。

3.2 工作面采出率驗證

通過統(tǒng)計工作面采出煤量，結(jié)合煤層厚度趨勢面方程計算工作面煤量和頂煤量，得到工作面采出率和頂煤回采率，計算得到的工作面采出率和頂煤回采率與金莊煤礦實際采出率和回采率比較進而驗證趨勢面方程預(yù)測煤層厚度變化的正確性。

圖6 工作面開采期間鉆孔布置示意圖

圖7 現(xiàn)場/趨勢面方程鉆孔煤厚和誤差率統(tǒng)計

3.2.1 工作面采出煤量統(tǒng)計

2013年2月到9月的實際生產(chǎn)情況統(tǒng)計見表3，統(tǒng)計當(dāng)月份的進刀數(shù)、頭進尺、尾進尺、月進尺及該月實際煤炭產(chǎn)量。

表3 工作面實際生產(chǎn)情況統(tǒng)計表Tab.3 Statistics of actual production condition in working face

3.2.2 工作面采出率和頂煤回收率計算

8203工作面采用綜采放頂煤開采方法，2013年2月到8月總進尺541.1 m(9月未統(tǒng)計進尺數(shù))，工作面地質(zhì)儲量計算為419.22萬t，工作面凈出煤量352.74萬t，工作面平均采出率84.14%，工作面理論頂煤量304.25萬t，工作面實際頂煤量248.60萬t，頂煤回采率81.71%，工作面頂煤采出情況統(tǒng)計結(jié)果見表4.

3.2.3 工作面采出率和頂煤回采率比較

金莊煤礦實際工作面回采率為85.24%，頂煤回采率80.69%.趨勢面方程計算得到工作面回采率為84.14%，頂煤回收率為81.71%.趨勢面方程工作面回采率誤差率1.29%，頂煤回采率誤差率為1.26%.金莊煤礦工作面實際采出率與頂煤回采率和趨勢面方程比較可得趨勢面方程能夠準(zhǔn)確預(yù)測工作面煤層厚度變化，對首采工作面煤厚變化具有指導(dǎo)意義。

表4 工作面采出率統(tǒng)計Tab.4 Statistics of the top coal extraction ratio in working face

4 結(jié) 論

1)根據(jù)工作面開采前期鉆孔數(shù)據(jù)建立趨勢面方程并驗證，得到一次趨勢面方程擬合度41.31%，四次趨勢面方程擬合度74.79%，且擬合度隨著次數(shù)的增加而升高；趨勢面方程F>F0.01和F>F0.05均滿足，建立的趨勢面方程均顯著；

2)趨勢面方程分析得到煤層厚度主要分布在14.0～20.0 m之間，煤層厚度變化趨勢總體不均勻，在工作面推進方向上，煤層厚度先降低，再增高，再降低，再增高，最后降低，工作面推進方向中部煤層厚度最大，大巷附近煤層厚度最小，煤層的相對富煤區(qū)總體呈現(xiàn)“Z”形，依據(jù)“Z”形逐步降低，局部呈現(xiàn)下凹，趨勢面分析結(jié)果與煤層等厚線圖基本吻合；工作面中線四次趨勢面煤層厚度整體在15.0～17.0 m變化，局部區(qū)域煤層厚度增大，在距離切眼長度為700～1 000 m間，煤層厚度急劇減小后急劇變大，煤層厚度起伏變化很大；

3)工作面開采期間鉆孔數(shù)據(jù)顯示和四次趨勢面方程預(yù)測數(shù)據(jù)最大正負誤差率為-9.5%和9.1%，最小正負誤差率為-1.6%和0.8%.開采期間鉆孔數(shù)據(jù)和四次趨勢面方程計算煤層厚度最大值和最小值所在位置相同，分別距離切眼位置為1 400 m和100 m，且誤差率分別為1.9%和3.7%，能夠準(zhǔn)確的預(yù)測首采面煤層分布情況和規(guī)律；

4)四次趨勢面方程計算得到工作面平均采出率84.14%，頂煤回采率81.71%；實際工作面回采率為85.24%，頂煤回采率80.69%，兩者誤差率為1.29%和1.26%，證明趨勢面方程預(yù)測煤層厚度變化的方法是正確有效的。

References

[1]杜文鳳，彭蘇萍．利用地質(zhì)統(tǒng)計學(xué)預(yù)測煤層厚度[J]．巖石力學(xué)與工程學(xué)報，2010，29(增1)：2 762-2 767．

DU Wen-feng，PENG Su-ping．Coalseam thickness prediction with geostatistics[J]．Chinese Journal of Rock Mechanics and Engineering，2010，29(Supp 1)：2 762-2 767．

[2]李增學(xué)，魏久傳，劉 瑩．煤田地質(zhì)學(xué)[M]．北京：地質(zhì)出版社，2005．

LI Zeng-xue，WEI Jiu-chuan，LIU Ying．Coal geology[M]．Beijing：Geological Press，2005．

[3]賀 清，董康乾，楊 軍．神府煤田沙溝岔煤礦綜采適應(yīng)性評價[J]．西安科技大學(xué)學(xué)報，2013，33(6)：662-668．

HE Qing，DONG Kang-qian，YANG Jun．Evaluation on coal mining conditions and fully mechanized mining suitability of Shagoucha coal field[J]．Journal of Xi’an University of Science and Technology，2013，33(6)：662-668．

[4]王 峰．邯鄲地區(qū)云陶擴大勘探區(qū)煤賦存規(guī)律研究[D]．邯鄲：河北工程大學(xué)，2011．

WANG Feng．Yun tao expansion exploration of the Handan areas coal deposited law research[D]．Handan：Hebei University of Engineering，2011．

[5]韓艷玲．核桃峪勘探區(qū)煤層賦存及煤質(zhì)變化規(guī)律研究[D]．邯鄲：河北工程大學(xué)，2011．

HAN Yan-ling．The discipline of the coal seam occurrence and coal quality variation in the exploration area of Hetaoyu[D]．Handan：Hebei University of Engineering，2011．

[6]郭鵬程．邢臺礦區(qū)南部構(gòu)造及煤層厚度變化規(guī)律研究[D]．西安：西安科技大學(xué)，2010．

GUO Peng-cheng．The study on variation of structure and coal-bed thickness in southern of Xingtai mining area[D]．Xi’an：Xi’an University of Science and Technology，2010．

[7]康永尚，沈金松，諶卓恒．現(xiàn)代數(shù)學(xué)地質(zhì)[M]．北京：石油工業(yè)出版社，2005．

KANG Yong-shang，SHEN Jin-song，CHEN Zhuo-heng．Modern machematical geology[M]．Beijing：Petroleum Industry Press，2005．

[8]徐振邦，婁元仁．?dāng)?shù)學(xué)地質(zhì)基礎(chǔ)[M]．北京：北京大學(xué)出版社，1994．

XU Zhen-bang，LOU Yuan-ren．Mathematical geology basis[M]．Beijing：Peking University Press，1994．

[9]李 翔，王金安，張少杰．復(fù)雜地質(zhì)體三維數(shù)值建模方法研究[J]．西安科技大學(xué)學(xué)報，2012，32(6)：676-681．

LI Xiang，WANG Jin-an，ZHANG Shao-jie．3D modeling method of complicated geological body[J]．Journal of Xi’an University of Science and Technology，2012，32(6)：676-681．

[10]商鐵林．綜放工作面頂煤采出率研究[D]．西安：西安科技大學(xué)，2009．

SHANG Tie-lin．Study on coal recovery ratio in fully mechanized top-coal caving face[D]．Xi’an：Xi’an University of Science and Technology，2009．

[11]李增學(xué)，馬興祥．實用礦井地質(zhì)研究方法與進展[M]．徐州：中國礦業(yè)大學(xué)出版社，1993．

LI Zeng-xue，MA Xing-xiang．The utility of mine geological research methods and progress[M]．Xuzhou：China University of Mining and Technology Press，1993．

[12]王繼承．井田構(gòu)造變化規(guī)律及其對煤層開采的影響[J]．煤礦開采，2002，52(增)：5-6．

WANG Ji-cheng．Variation of Mine field structure and its influence on Mining[J]．Coal Mining，2002，52(Suppl.)：5-6．

[13]路 鵬，周榮福，王 玫．基于趨勢面分析的煤層空間分布特征研究[J]．中州煤炭，2009，161(5)：9-10．

LU Peng，ZHOU Rong-fu，WANG Mei．Research on coal seams spatial distribution features based on trend surface analysis[J]．Zhongzhou Coal，2009，161(5)：9-11．

[14]李小明，劉德民，張品剛．趨勢面分析在東歡坨井田構(gòu)造研究中的應(yīng)用[J]．中國煤炭地質(zhì)，2010，22(4)：6-9．

LI Xiao-ming，LIU De-min，ZHANG Pin-gang．Application of trend surface analysis in Donghuantuo mine field structural study[J]．Coal Geology of China，2010，22(4)：6-9．

[15]李學(xué)錫，周榮福，李小麗．趨勢面分析在東灘礦十四采區(qū)北部3下煤層的應(yīng)用[J]．西部探礦工程，2008(11)：135-138．

LI Xue-xi，ZHOU Rong-fu，LI Xiao-li．Application of trend surface analysis in north of fourteen mining area 3 lower coal seam of Dongtan mine [J]．West China Exploration Engineering，2008(11)：135-138．

[16]曹代勇，劉欽甫，彭蘇萍，等．超化井田煤厚變化的構(gòu)造控制因素研究[J]．中國礦業(yè)大學(xué)學(xué)報，1998，27(4)：393-397．

CAO Dai-yong，LIU Qin-fu，PENG Su-ping，et，all．Teconic control over thickness variation of coal seam in Chaohua mine field[J]．Journal of China University of Mining & Technology，1998，27(4)：393-397．

[17]馬田生，張林山．趨勢面分析在山西朔州王坪井田構(gòu)造研究中的應(yīng)用[J]．新疆地質(zhì)，2004，22(1)：108-111．

MA Tian-sheng，ZHANG Lin-shan．Application of trend surface analysis in the study structure in of Shanxi Shuozhou Wangping coal mine[J]．Xinjiang Geology，2004，22(1)：108-111．

[18]孫學(xué)陽，楊 忠．陜北神木北部礦區(qū)構(gòu)造發(fā)育規(guī)律[J]．黑龍江科技學(xué)院學(xué)報，2008，18(2)：87-91．

SUN Xue-yang，YANG Zhong．Law goverming geological structure development of Shenmu coal mine area in northern Shaanxi[J]．Journal of Heilongjiang Institute of Science & Technology，2008，18(2)：87-91．

[19]劉衡秋，劉欽甫，付 正，等．淮南潘三礦區(qū)13-1#煤層空間分布特征及穩(wěn)定性研究[J]．煤田地質(zhì)與勘探，2004，32(1)：21-23．

LIU Heng-qiu，LIU Qin-fu，F(xiàn)U Zheng，et，all．Study of the space distribution trait and stability of No.13-1 coal seam in the Pansan coal mine,Huainan[J]．Coal Geology & Exploration，2004，32(1)：21-23．

Research on coal seam thickness changing law based on trend surface analysis

ZHANG Yan-wei1,2,TU Shi-hao1,2,WANG Chen1,2,ZHANG Cun1,2,WEI Shuai-feng1,2,SONG Qi1,2

(1.SchoolofMines,ChinaUniversityofMining&Technology,Xuzhou221116,China;2.StateKeyLaboratoryofCoalResourceandMineSafety,ChinaUniversityofMining&Technology,Xuzhou221116,China)

Based on engineering practice in 8203 working face of Jinzhuang coal mine，through the statistics of working face coal seam borehole data，predicting thickness change trend of 5(3～5)#coal seam by using the trend surface analysis method，trend surface equation of coal seam thickness has been obtained.The coal seam thickness trend surface fitting degree has been checked by extraction ratio F.Finally，using the engineering verification method to prove validity and reliability of predicting coal seam thickness changing law by trend surface analysis.Four conclusions were drawn： First trend surface of coal seam fitting degree is 41.31% and fourth is 74.79%.During working face mining，the maximum plus or minus error rate is 9.1% and -9.5%，the smallest positive and negative error rate is 0.8% and -1.6%，and consistent with the actual drilling measuring of coal seam thickness.Fourth trend surface equation calculation indicates that the average recovery coefficient is 84.14% and the top coal recovery rate is 81.71% in working face.The actual recovery coefficient is 85.24% and the top coal recovery rate is 80.69% in working face，both error rate is 1.29% and 1.26%.It is correct and reliable to predicte coal seam thickness changing law by trend surface analysis.

super thick coal seam；trend surface analysis;variation in thickness of coal seam； F distribution

10.13800/j.cnki.xakjdxxb.2015.0104

1672-9315(2015)01-0021-07

2014-10-20 責(zé)任編輯：劉 潔

國家高技術(shù)研究發(fā)展計劃(863)項目(2012AA062101)；江蘇省高校優(yōu)勢學(xué)科建設(shè)項目(PAPD)；國家自然科學(xué)

基金項目(51374200)

張艷偉(1990-)，男，山西呂梁人，碩士研究生，E-mail：zhangyanwei712@163.com

A