徐義洋，褚廷湘，馮學(xué)文，代文程

(1. 北京中煤礦山工程有限公司，北京 100013；2. 華北科技學(xué)院 安全監(jiān)管學(xué)院，河北 三河 065201；3. 中國礦業(yè)大學(xué)(北京) 機電與信息工程學(xué)院，北京 100083； 4. 安徽省皖北煤電集團(tuán)臨汾天煜恒昇煤業(yè)有限責(zé)任公司，山西 臨汾 041204)

井工煤礦其煤自燃主要發(fā)生在采空區(qū)、破裂煤柱等地點，其中以采空區(qū)的發(fā)火次數(shù)最多、造成的損失也最嚴(yán)重。因此采空區(qū)自然發(fā)火的研究，一直都是礦井火災(zāi)防治領(lǐng)域所關(guān)注的重點方向。在煤自燃發(fā)生及發(fā)展過程中，涉及到流固介質(zhì)之間的傳質(zhì)及傳熱過程，而在這一過程中，空隙率對煤自燃的發(fā)生及發(fā)展進(jìn)程具有內(nèi)在的作用。王繼仁等應(yīng)用分形理論，分析了氧氣在散體煤中的擴散問題，建立了氧氣在散體煤中的分形擴散方程和散體煤顆粒分形表面氧氣吸附速率方程。高建良和王海生根據(jù)多孔介質(zhì)滲流理論，分別模擬了采空區(qū)滲透率為均勻、分段均勻和連續(xù)性分布情況下的采空區(qū)漏風(fēng)流場，發(fā)現(xiàn)漏風(fēng)量、漏風(fēng)風(fēng)速分布和自燃三帶的位置和寬度有很大差別。李樹剛和錢鳴高通過理論分析和現(xiàn)場觀測，分析了綜放采空區(qū)垮落巖體碎脹特性與采空區(qū)流態(tài)的影響。楊永良等利用數(shù)值模擬及實測的方法，分析了煤巖垮落規(guī)律與采空區(qū)自燃規(guī)律之間的關(guān)系。自分形理論提出以來，在礦業(yè)領(lǐng)域此理論也得到了廣泛使用。謝和平等運用分形理論研究了巖石的斷裂和破碎問題，及綜放工作面巷道裂隙發(fā)育情況。李宏艷等通過采動條件下裂隙場演化實驗室相似模擬實驗，利用分形理論定量描述了裂隙時空的演化規(guī)律。鄭克洪等根據(jù)分形理論建立了煤矸破碎粒度分布的分形表達(dá)式。姜秀民等根據(jù)實驗數(shù)據(jù)研究了煤粉顆粒的分形特征。KATZ等把分形幾何理論用于分析多孔介質(zhì)內(nèi)部的孔隙結(jié)構(gòu)。他們的研究表明：多孔介質(zhì)的孔隙空間和孔隙界面都具有分形特征，可以用分形維數(shù)來預(yù)測多孔介質(zhì)的孔隙度。以上研究可以看出不論是從宏觀巖石破斷、覆巖裂隙、煤矸破碎，還是在微觀尺度分形理論描述空隙率也是可行的。

上述研究對描述采空區(qū)破碎煤巖空隙率提供了途徑，但是由于采空區(qū)垮落的不規(guī)則性及采空區(qū)內(nèi)部的不可接近性，造成對采空區(qū)自燃環(huán)境的識別比較局限。而遺煤自燃的發(fā)生卻與其所在的堆積狀態(tài)及環(huán)境條件息息相關(guān)，在自燃發(fā)生及發(fā)展過程中滲流場、氧濃度場、溫度場之間的相互作用、相互聯(lián)系，這其中最關(guān)鍵的問題在于采空區(qū)破碎遺煤的堆積結(jié)構(gòu)的復(fù)雜性，此外進(jìn)行全尺度的研究的設(shè)備手段還不能實現(xiàn)。因此，筆者基于煤自燃角度及采空區(qū)遺煤垮落特征，剔除矸石堆積部分，只關(guān)注垮落帶內(nèi)破碎遺煤承壓下空隙率演變，結(jié)合分形理論，建立了顆粒煤堆積體承壓下的空隙率模型，并對加載應(yīng)力條件下顆粒煤的空隙率演化進(jìn)行嘗試。

1 空隙率對煤自燃進(jìn)程影響分析

采空區(qū)遺煤自燃與覆巖垮落關(guān)系密切，煤巖垮落狀態(tài)決定了采空區(qū)破碎遺煤的堆積狀態(tài)，而遺煤堆積空間結(jié)構(gòu)及變化，直接影響煤的蓄熱環(huán)境和采空區(qū)的漏風(fēng)狀態(tài)及氧氣分布。

1.1 覆巖垮落與堆積

據(jù)巖層的移動特征及礦壓理論將煤層的上覆巖層在垂直方向上分為垮落帶、裂隙帶、彎曲下沉帶。工作面后方的采空區(qū)，初次來壓之前只有采空卸荷區(qū)，初次來壓之后開始逐漸形成應(yīng)力升高區(qū)和應(yīng)力恢復(fù)區(qū)，但初次來壓后的工作面與應(yīng)力升高區(qū)之間仍存在采空卸荷區(qū)。采空區(qū)覆巖垮落與應(yīng)力分布如圖1所示，其中，為工作面走向卸壓區(qū)寬度，為工作面走向應(yīng)力集中區(qū)寬度，為采空區(qū)走向應(yīng)力逐漸回復(fù)長度，m。在緊靠工作面后方，受支架及上部未完全破斷的巖梁掩護(hù)，未受到基本頂來壓的應(yīng)力影響，浮煤及矸石處于自然堆積狀態(tài)，稱之為自然堆積區(qū)；隨著采空區(qū)應(yīng)力逐漸恢復(fù)，浮煤及矸石不再成自然狀態(tài)堆積，在壓應(yīng)力作用下浮煤及矸石進(jìn)一步破碎，密實度增加，稱為破碎堆積區(qū)；當(dāng)隨著采空區(qū)深度的增加，采空區(qū)底板應(yīng)力恢復(fù)到原巖應(yīng)力水平而重新壓實區(qū)，密實度最大，稱為壓實區(qū)。

圖1 采空區(qū)覆巖垮落及應(yīng)力分布示意Fig.1 Diagram of the strata caving and stress distribution in goaf

采空區(qū)遺煤先后經(jīng)歷卸壓應(yīng)力區(qū)、應(yīng)力恢復(fù)區(qū)及原巖應(yīng)力區(qū)3個階段，在每個階段遺煤所承載的應(yīng)力不同，致使遺煤堆積結(jié)構(gòu)發(fā)生變化，進(jìn)而影響遺煤堆積區(qū)域的空隙率。從煤自燃角度分析，進(jìn)而影響到采空區(qū)的漏風(fēng)狀態(tài)、漏風(fēng)強度及溫度場環(huán)境，即影響采空區(qū)遺煤的自然發(fā)火的進(jìn)程。文獻(xiàn)[15]在實驗室內(nèi)通過對散體煤進(jìn)行加載實驗，得到了軸向應(yīng)力與碎脹系數(shù)的關(guān)系，發(fā)現(xiàn)隨著軸向應(yīng)力增加，碎脹系數(shù)降低，逐漸趨于壓實，進(jìn)而改變破碎煤體的空隙率。所以如何確定采空區(qū)承壓破碎煤體在應(yīng)力變化條件下的空隙率響應(yīng)特征，對研究采空區(qū)煤自燃的發(fā)生顯得尤為重要。

采空區(qū)覆巖應(yīng)力對空隙率的影響是間接的，是通過影響采空區(qū)煤巖的碎脹系數(shù)來影響空隙率的，空隙率與碎脹系數(shù)存在以下關(guān)系：

=1-1

(1)

由式(1)可以看出，空隙率與碎脹系數(shù)成正比關(guān)系：碎脹系數(shù)越大，空隙率亦越大。

1.2 空隙率與煤自燃關(guān)系

對于一個特定的綜放工作面，其通風(fēng)方式、供風(fēng)量、風(fēng)壓分布等基本確定，影響漏風(fēng)強度的主要因素是松散煤體的空隙率。漏風(fēng)強度是風(fēng)壓梯度、空隙率的函數(shù)。

(2)

由式(2)可以看出，當(dāng)漏風(fēng)通道兩端壓差一定情況下，空隙率越大，漏風(fēng)量越大，空隙率影響漏風(fēng)區(qū)域的漏風(fēng)強度。在采空區(qū)漏風(fēng)強度一定的情況下，采空區(qū)遺煤堆積環(huán)境的氧氣體積分?jǐn)?shù)的分布與空隙率分布有關(guān)。

(3)

由式(3)可以看出，采空區(qū)某點的氧氣體積分?jǐn)?shù)與空隙率的分布有關(guān)。

根據(jù)能量守恒定律，考慮氣固兩相時，采空區(qū)遺煤的溫度方程為

(4)

由式(4)可以看出，采空區(qū)遺煤堆積空間的空隙率對溫度場的演化具有重要影響。

綜上，可見空隙率對采空區(qū)漏風(fēng)強度、氧氣體積分?jǐn)?shù)的分布、溫度場的分布影響較大，因此基于采空區(qū)應(yīng)力恢復(fù)特征，如何定量描述其在應(yīng)力擾動下的空隙率計算模型，對煤自燃研究來說具有一定的意義。

2 承壓破碎遺煤的分型模型

采空區(qū)遺煤堆積空間四周受力，軸向覆巖應(yīng)力作用為主動應(yīng)力，同時由于遺煤散體的變形、滑動受兩側(cè)煤壁、老開切眼及支架的束縛限制，宏觀表現(xiàn)為以軸向壓實為主，四周相對位移固定的受力狀態(tài)。在承壓破碎散體滲透性方面劉衛(wèi)群、李順才等研究了不同破碎程度下壓實的滲透系數(shù)變化規(guī)律，指出破碎巖石的滲透性不僅和外載水平及初始孔隙度有關(guān)，還與碎石的粒徑、顆粒的排列方式以等有關(guān)。PAPPAS和MARK通過20組單軸壓縮實驗，詳細(xì)研究了采空區(qū)垮落材料的力學(xué)性質(zhì)，通過應(yīng)力-應(yīng)變關(guān)系，獲得了有關(guān)巖石的力學(xué)參數(shù)及散體破斷特征。梁軍等采用大型壓縮儀進(jìn)行了堆石蠕變試驗和顆粒破碎測試，指出堆石破碎有主壓縮破碎和蠕變破碎2類，其中蠕變破碎是堆石料骨架應(yīng)力的重新分布導(dǎo)致粗大顆粒棱角或者軟弱顆粒發(fā)生破碎、細(xì)化，顆粒排列進(jìn)一步調(diào)整，如圖2所示。

圖2 顆粒破碎類型Fig.2 Particle breakage types

采空區(qū)遺煤依次經(jīng)歷應(yīng)力三區(qū)時，會再次發(fā)生破碎，遺煤顆粒破碎在覆巖荷載作用下產(chǎn)生結(jié)構(gòu)的破裂或破損，進(jìn)而影響其堆積結(jié)構(gòu)與狀態(tài)。

2.1 遺煤堆積空隙率分形模型

TYLER等提出了質(zhì)量和孔徑關(guān)系的標(biāo)準(zhǔn)化方程，且通過假設(shè)不同土壤粒級具有相同密度的條件，建立了土壤粒徑分布分形維數(shù)的質(zhì)量分布模型。由于采空區(qū)遺煤粒徑、形狀不一，根據(jù)TYLER的研究結(jié)果，可以通過不同直徑的篩子篩分試樣，通過篩下試樣的質(zhì)量和總試樣的質(zhì)量建立函數(shù)關(guān)系，分析遺煤粒度的分布。

設(shè)()為篩下試樣質(zhì)量，()為篩上試樣質(zhì)量，為試樣總質(zhì)量，定義關(guān)聯(lián)函數(shù)()為

(5)

在三維空間中，大于某一粒度的粗粒顆粒所構(gòu)成的體積為，則有

=[1-()3-]

(6)

式中，為破碎分形維數(shù)；，均為常數(shù)，與遺煤顆粒的大小、形狀有關(guān)。

假定遺煤顆粒具有相同的密度，則大于某一粒徑的遺煤質(zhì)量為

()=[1-()3-]

(7)

若令=0，可得到試樣的總質(zhì)量；若令等于試樣中的最大粒徑，則有()=0，代入式(7)可得=，則可得到遺煤顆粒質(zhì)量與粒徑之間的分形關(guān)系:

(8)

對式(8)兩邊取對數(shù)，則有

ln[()]=(3-)ln()

(9)

通過物理篩分、稱重統(tǒng)計，可得到()的統(tǒng)計值，并依據(jù)線性回歸求得直線部分的斜率，即(3-)，進(jìn)而得到遺煤破碎分形維數(shù)。

采空區(qū)垮落遺煤空隙率為破碎遺煤間的空隙體積與整個垮落堆積體體積之比，設(shè)垮落遺煤的顆粒尺寸為，根據(jù)密度定義，則有

(10)

其中，()為顆粒特征尺寸小于的遺煤分形體積，m；為垮落巖體總體積，m，為破碎煤體密度，kg/m。則對于粒徑小于部分遺煤堆積體積為

(11)

則，粒徑處于區(qū)間(，+d)內(nèi)的堆積體積為

(12)

則，垮落遺煤的分形總體積為

(13)

則有

(14)

式中，為垮落遺煤自然碎脹空隙率；，分別為最小及最大顆粒粒徑；為煤的質(zhì)量；為自然碎脹遺煤密度；為煤層的密度。

2.2 承壓破碎煤體空隙率演變模型

采空區(qū)遺煤堆以軸向覆巖應(yīng)力作用為主動應(yīng)力，同時由于遺煤散體的變形、滑動受兩側(cè)煤壁、開切眼及支架的束縛限制，宏觀表現(xiàn)為以軸向壓實為主，因此可以通過側(cè)限壓實實驗，輔助應(yīng)力-應(yīng)變關(guān)系，得到承壓破碎遺煤的視密度與真密度之間的關(guān)系。

軸向應(yīng)變?yōu)?/p>

=Δ

(15)

式中，Δ為壓縮變形的煤體高度；為壓縮變形后煤體高度。

假設(shè)在側(cè)限壓實實驗過程中完整煤樣初始質(zhì)量為，初始碎脹密度為，自然碎脹堆積裝置內(nèi)煤樣高度為；承載破碎煤體的密度為，承壓過程中應(yīng)變?yōu)?。則承壓后煤樣在不同應(yīng)力狀態(tài)下的高度=(1-)，破碎煤體承壓過程中，不斷被壓實，空隙體積隨壓實程度不斷減小，由質(zhì)量守恒得

=1(1-)

(16)

根據(jù)空隙率定義，則承壓破碎煤體的空隙率，可表達(dá)為

(17)

式中，為顆粒堆積體的空隙體積，m；為顆粒堆積體表觀體積，m。

式(17)有效使用，需滿足：≤，理論上，數(shù)值上當(dāng)=時，煤樣完全被壓實，無任何空隙。

在壓實過程中應(yīng)力-應(yīng)變的關(guān)系可表達(dá)為

=e

(18)

式中，為軸向壓應(yīng)力；，為壓實回歸系數(shù)。

承壓破碎遺煤的空隙率分型模型為

(19)

3 承壓破碎煤體空隙率實驗

3.1 實驗煤樣及方法

所選煤樣來自山西恒昇煤業(yè)9106工作面，在該工作面傾向支架后部3個不同位置取樣，選取15 mm以下粒徑的混合煤樣作為試驗煤樣，經(jīng)自然風(fēng)干后，進(jìn)行篩分，確定了自然碎脹區(qū)遺煤15 mm以下遺煤堆積體的級配的構(gòu)成，見表1,其中為顆粒粒徑。

表1 初始各粒徑組質(zhì)量

其中相同位置煤樣各準(zhǔn)備2份分別標(biāo)記為A，B，C，A，B，C，保持級配相同，一份用于應(yīng)力應(yīng)變測試，一份用于特定軸向加載應(yīng)力下的分形維數(shù)確定。

利用實驗室內(nèi)自主設(shè)計的承壓破碎遺煤滲透率演化及自燃測試裝置，如圖3所示。通過此裝置，可對采空區(qū)篩選的煤樣開展軸向加載實驗，通過設(shè)定不同的軸向加載強度，通過重新篩分，可確定加載后煤樣的級配構(gòu)成及承壓破碎煤樣的空隙率。

圖3 承壓破碎煤體空隙率測試實驗裝置Fig.3 Porosity test device of the compacted broken coal

3.2 原始煤樣分形維數(shù)確定

結(jié)合表1，對煤樣粒徑與對應(yīng)篩下質(zhì)量比，通過對數(shù)取值，然后線性回歸，根據(jù)擬合曲線斜率，通過式(9)可回歸得到原始煤樣的破碎分形維數(shù)，如圖4所示，其中，為粒徑以下的煤樣質(zhì)量。通過對原始煤樣的級配篩分及分形維數(shù)回歸，利用空隙率分形模型式(14)，得到了自然堆積狀態(tài)下的A，B，C煤樣的空隙率(表2)。發(fā)現(xiàn)利用分形模型計算得到的空隙率與空隙率定義計算出的空隙率誤差在0.045～0.086，空隙率分形模型可用于采空區(qū)遺煤孔隙率計算。

圖4 原始煤樣ln(di/dmax)與ln(mi/mt)回歸Fig.4 Regression curves of the initial coal sample between ln(di/dmax)andln(mi/mt)

表2 原始煤樣分形維數(shù)與空隙率

3.3 承壓破碎煤體的分形維數(shù)

分別對煤樣A，B，C進(jìn)行了承壓實驗，通過回歸分析確定了3個煤樣的應(yīng)力-應(yīng)變關(guān)系，見表3。

由表3可知破碎煤體應(yīng)力-應(yīng)變符合指數(shù)關(guān)系。隨著應(yīng)力增加，破碎煤樣的應(yīng)變逐漸增加，即煤樣逐漸被壓實，結(jié)合圖1可知在采空區(qū)橫三區(qū)內(nèi)，隨著采空區(qū)應(yīng)力逐漸恢復(fù)，采空區(qū)遺煤堆積體逐漸被壓實。

表3 煤樣承壓應(yīng)力應(yīng)變擬合關(guān)系

對不同加載應(yīng)力狀態(tài)下的破碎煤樣再次篩分，通過分形理論可得到破碎煤樣的承壓破碎分形維數(shù)。

分別測試了煤樣A，B，C在加載8，12及15 MPa后煤樣承載破碎后的級配情況，并對其回歸，得到加載后A，B，C的破碎分形維數(shù),如圖5所示。發(fā)現(xiàn)承壓破碎煤體分形維數(shù)增大則說明受載煤樣更加破碎，在應(yīng)力加載狀態(tài)下其空隙率逐漸減小。

圖5 承壓破碎后煤樣ln(di/dmax)與ln(mi/mt)回歸Fig.5 Regression curves after compacted between ln(di/dmax)andln(mi/mt)

3.4 模型與實驗結(jié)果對比

通過對煤樣A，B，C的軸向加載實驗，根據(jù)其煤樣質(zhì)量、煤樣真密度、壓實位移，可得到不同應(yīng)力狀態(tài)下的實驗空隙率；通過表2，4的實驗數(shù)據(jù)，根據(jù)式(19)，可確定其分形理論空隙率，見表5。

表4 軸向加載后級配及分形維數(shù)

表5 承壓煤樣分形空隙率計算結(jié)果

經(jīng)實驗測試及理論計算，得到了煤樣A，B，C分別在加載8，12，15 MPa下的空隙率，通過2者比較，3個煤樣的誤差分別為0.106，0.028，0.062，說明通過承壓破碎遺煤的空隙率分形模型是可以用于預(yù)測不同應(yīng)力狀態(tài)下的空隙率變化。

3.5 承壓破碎煤體空隙率預(yù)測方法

以B煤樣的應(yīng)力應(yīng)變曲線為例，結(jié)合承壓破碎遺煤空隙率分形模型，可對不同應(yīng)力狀態(tài)下的破碎煤體空隙率演化進(jìn)行預(yù)測，(式(20))。

(20)

承壓破碎煤體空隙率預(yù)測方法可簡述為：首先獲得其應(yīng)力-應(yīng)變曲線，然后根據(jù)煤樣的級配構(gòu)成確定其破碎分形維數(shù)后，根據(jù)所建立的承壓破碎煤體空隙率分形模型可對其預(yù)測。

4 結(jié) 論

(1)以空隙率對煤自燃影響為切入點，分析了采空區(qū)空隙率對煤自燃發(fā)生中的漏風(fēng)強度、氧氣體積分?jǐn)?shù)分布、蓄熱環(huán)境的影響，明確了空隙率變化對采空區(qū)煤自燃的影響作用。

(2)結(jié)合采空區(qū)應(yīng)力恢復(fù)及承壓破碎煤體壓實特征，通過引入分形理論，建立了承壓破碎煤體空隙率分形模型。

(3)通過承壓破碎煤體壓實實驗，明確了壓實前后分形維數(shù)的變化，基于實驗數(shù)據(jù)與推到模型進(jìn)行了驗證，實驗結(jié)果與理論計算結(jié)果誤差不大，所建立的模型具有一定的可靠性。通過該研究，為后續(xù)研究恒昇煤業(yè)采空區(qū)漏風(fēng)、煤自燃氧化區(qū)域及瓦斯運移提供了理論指導(dǎo)。

(4)因受限于采空區(qū)破碎煤巖全尺度樣本索取困難，本文僅以小樣本尺度煤樣為對象進(jìn)行了研究，這對于研究采空應(yīng)力恢復(fù)狀態(tài)下的空隙率演化提供了一種新思路。