劉義磊 黃素果

摘要：針對采空區(qū)瓦斯運移包括滲流和非均質(zhì)氣體的擴散兩個方面的問題，利用滲流理論、氣體擴散理論為基礎(chǔ)建立起采空區(qū)風流流動數(shù)學模型，以及采空區(qū)瓦斯運移的數(shù)學模型，并采用計算機數(shù)值求解，將數(shù)值模擬結(jié)果與現(xiàn)場實測數(shù)據(jù)結(jié)果進行對比分析，提出了采空區(qū)瓦斯?jié)舛确植继卣骷巴咚惯\移規(guī)律。

關(guān)鍵詞：綜放工作面 采空區(qū) 瓦斯 氣象參數(shù)

中圖分類號：TD712 文獻標識碼：A 文章編號：1672-3791（2012）10（b）-0048-02

采空區(qū)瓦斯?jié)舛确植疾坏芰鲌鏊俣确植加绊?，同時流場的速度分布又反過來影響濃度分布，這個兩個問題的方程需要耦合求解。本文利用滲流理論、氣體擴散理論為基礎(chǔ)建立風流流動及瓦斯運移的數(shù)學模型，采用計算機數(shù)值求解，分析研究采場氣體流動和瓦斯分布規(guī)律。

1 采空區(qū)風流流動數(shù)學模型

在礦井通風系統(tǒng)中，風流的流動可分為管道流動和采空區(qū)冒落區(qū)域內(nèi)的流動。在采場內(nèi)的管道流動時風壓穩(wěn)定，各處的壓力差值相對很小，不考慮溫度的變化影響，此時風流流動可以看作不可壓縮流體的穩(wěn)定流動，符合連續(xù)方程。而采空區(qū)的冒落區(qū)域被破碎的巖石填滿，巖石之間形成多處縫隙，風流在這些縫隙中流動的速度很低，可看作是滲流，在忽略結(jié)構(gòu)變形的影響時，符合達西定律：。式中：V為滲流速度；為滲流場壓力；k為滲流系數(shù)。因此可以得到風流在采空區(qū)范圍內(nèi)穩(wěn)定流動的微分方程：

由于在實際的流場分布中風速測量起來比較簡單方便，因此選擇第二類邊界條件作為求解條件，此時該方程存在定解的充要條件是：

式中：為工作面或沿空巷與采空區(qū)交界面上采空區(qū)側(cè)任意兩點的壓頭差，m；

為工作面或沿空巷內(nèi)和相對應(yīng)的靜壓差，Pa；

為工作面或沿空巷內(nèi)和相對應(yīng)的標高差，m。

式（1）、（2）聯(lián)立即可對采空區(qū)風流流動的數(shù)學問題進行求解。

2 采空區(qū)瓦斯運移數(shù)學模型

采空區(qū)的氣體狀態(tài)可以視為一個由瓦斯和空氣這兩種元素構(gòu)成的系統(tǒng)，將瓦斯用字母來表示，將空氣用字母來表示。那么它們流速和密度分別用、、、來表示，因此可以定義該二元系統(tǒng)的平均流速為：

在不考慮熱擴散影響的情況下，該二元系統(tǒng)中描述其中一種物質(zhì)擴散的定律為：

式中：為密度，g/m3；為的擴散通量，g/（s·m2）；為擴散系數(shù)，m2/s。

在二元系統(tǒng)中取出一種物質(zhì)作為單元體，由質(zhì)量守恒定律可以得到二元系統(tǒng)中的連續(xù)方程式為：，式中：為頂?shù)装鍐挝幻娣e上向采空區(qū)單位高度上的瓦斯涌出速率，g/（s·m2）。將式4代入其中可以得到：

式（5）適用于二元系統(tǒng)中一個單元體的質(zhì)量守恒問題，當牽涉到瓦斯在縫隙結(jié)構(gòu)中的遷移問題時，只適用上面的方程式無法計算出結(jié)果的。由于縫隙結(jié)構(gòu)的復(fù)雜性，巖石顆粒表面的形狀和尺寸是幾乎無法具體給出的。因此，要將式（5）用于縫隙結(jié)構(gòu)中，就必須把方程中的各個物理量都適當?shù)剞D(zhuǎn)化取其平均值，從而導(dǎo)出縫隙結(jié)構(gòu)中的瓦斯運移方程。由于通風系統(tǒng)中瓦斯的濃度不會太高，因此可近似的看作常數(shù)，可以得出：

將式（6）代入式（5）可以得到：

為了表達上的方便，往往用字母來代替，定義其為瓦斯?jié)舛?；，稱為平均流速，這時式（7）的形式為：

其穩(wěn)定型為：

式（9）即為瓦斯在采空區(qū)的縫隙結(jié)構(gòu)中的擴散方程，從該方程的求解過程中可以看出：瓦斯在采空區(qū)的運移可以分解成由瓦斯?jié)舛纫鸬南鄬τ谄骄\移的分子擴散和對流造成的瓦斯運移。

3 采空區(qū)瓦斯流場的數(shù)值計算

根據(jù)上述數(shù)學模型和計算方法，采用計算機數(shù)值求解，即可求出采空區(qū)瓦斯的流動狀態(tài)及濃度的分布規(guī)律。根據(jù)礦井生產(chǎn)實際條件，可將采空區(qū)劃分為生產(chǎn)區(qū)和封閉區(qū)。其中生產(chǎn)采空區(qū)又可劃分為獨立生產(chǎn)采空區(qū)、鄰空生產(chǎn)采空區(qū)。從本文所研究的某礦5404和6404綜放面的實際情況看，5404采空區(qū)為獨立采空區(qū)，6404采空區(qū)為進風鄰空采空區(qū)。工作面長度為220m，采空區(qū)模擬深度為410m，采空區(qū)瓦斯?jié)舛葹?5%。在正常生產(chǎn)條件下，5404面的配風量1100m3/min，6404工作面的配風量1250m3/min，數(shù)值模擬結(jié)果見表1。

4 結(jié)論

將數(shù)值模擬的結(jié)果和現(xiàn)場實測數(shù)據(jù)進行對比分析可以得總結(jié)出以下結(jié)論。

（1）采空區(qū)瓦斯的分布和壓力分布是相互影響的。在靠近工作面的地方由于風量和風速都比較高，瓦斯迅速被帶走，其瓦斯?jié)舛认鄬^低；沿工作面方向，從進風側(cè)到回風側(cè)，瓦斯?jié)舛瓤偟内厔菔侵饾u增高的，到上隅角處增到最大。

（2）當受到相鄰采空區(qū)的漏風影響時，其進風區(qū)域瓦斯含量會有所升高。

（3）采空區(qū)內(nèi)部有高濃度瓦斯存在，從計算模擬結(jié)果可以看出高濃度瓦斯大部分處在采動影響區(qū)的后部，只有在回風的部分區(qū)域離工作面較近，應(yīng)加強防治。

（4）風量增加，雖然使得采空區(qū)漏風范圍和強度增大，但是并沒有明顯增加采空區(qū)瓦斯涌出的平均濃度。

參考文獻

[1] 程遠平，俞其香.中國煤礦區(qū)域性瓦斯治理技術(shù)的發(fā)展[J].采礦與安全工程學報，2007（12）.

[2] 兗礦集團有限公司.兗州礦區(qū)礦井通風安全技術(shù)[M].煤炭工業(yè)出版社，2001.

[3] 鮮學福，李曉紅，姜德義，等.煤與瓦斯突出的激發(fā)和發(fā)生條件[J].巖土力學，2009（3）.

[4] 譚允禎.礦井通風系統(tǒng)管理技術(shù)理論[M].煤炭工業(yè)出版社，1998.

[5] 朱委申，何滿潮.復(fù)雜條件下圍巖穩(wěn)定性與巖體動態(tài)施工力學[M].北京：科學出版社，1996.