吳奉亮 李智勝 常心坦

摘要：為模擬礦井采空區(qū)漏風(fēng)引起的煤自燃與瓦斯問(wèn)題，建立了采空區(qū)“瓦斯-氧”對(duì)流擴(kuò)散過(guò)程的迎風(fēng)有限元統(tǒng)一求解模型。分析了基于伽遼金法的采空區(qū)流場(chǎng)有限元求解技術(shù)。給出了采空區(qū)“瓦斯-氧”對(duì)流擴(kuò)散過(guò)程的統(tǒng)一模型，針對(duì)標(biāo)準(zhǔn)伽遼金有限元法求解統(tǒng)一模型時(shí)計(jì)算結(jié)果振蕩的問(wèn)題，采用迎風(fēng)有限元方法得到模型的穩(wěn)定解?；谌切螁卧x散的采空區(qū)，給出了以所有節(jié)點(diǎn)的壓力以及氣體濃度為未知數(shù)的控制方程組。使用VC++開發(fā)了模擬軟件，采用典型算例對(duì)軟件的穩(wěn)定性與實(shí)用性進(jìn)行了驗(yàn)證，模擬結(jié)果表明：煤自燃耗氧導(dǎo)致采空區(qū)氧濃度降低，可爆瓦斯?jié)舛确秶?%～15%不再適用，采空區(qū)內(nèi)具有瓦斯爆炸危險(xiǎn)性覆蓋的區(qū)域隨遺煤耗氧與瓦斯涌出情況而動(dòng)態(tài)變化;開發(fā)的軟件可以直接劃定采空區(qū)內(nèi)煤自燃三帶及具有爆炸危險(xiǎn)性的采空區(qū)區(qū)域。

關(guān)鍵詞：礦業(yè)工程;采空區(qū)漏風(fēng);煤自燃;瓦斯;迎風(fēng)有限元

中圖分類號(hào)：TD 75+2.2文獻(xiàn)標(biāo)志碼：A

DOI：10.13800/j.cnki.xakjdxxb.2019.0208文章編號(hào)：1672-9315（2019）02-0234-07

0引言

煤礦采空區(qū)漏風(fēng)容易引發(fā)煤自燃，形成可爆瓦斯區(qū)域[1]甚至引發(fā)瓦斯爆炸。由于對(duì)采空區(qū)進(jìn)行直接觀測(cè)十分困難，數(shù)值模擬成為國(guó)內(nèi)外學(xué)者研究采空區(qū)漏風(fēng)問(wèn)題的有效手段，如判定采空區(qū)內(nèi)的煤自燃三帶[2-3]，模擬煤自燃過(guò)程[4-5]及瓦斯運(yùn)移規(guī)律[6-7]。目前模擬采空區(qū)漏風(fēng)問(wèn)題的途徑主要有采用商業(yè)CFD與自編程2種，相關(guān)的數(shù)值方法包括有限元、有限體積法等。盡管商業(yè)CFD（如Fluent，COMSOL等）已被大量學(xué)者使用，但它對(duì)理論知識(shí)要求高，目前還很難被現(xiàn)場(chǎng)人員所接受。自編程研究采空區(qū)漏風(fēng)問(wèn)題很可能是將CFD技術(shù)推向現(xiàn)場(chǎng)應(yīng)用的有效途徑，因?yàn)橥ㄟ^(guò)自編程可將礦井通風(fēng)問(wèn)題的相關(guān)概念作為軟件的界面元素進(jìn)行設(shè)計(jì)，將CFD算法封裝到后臺(tái)[8]，以此來(lái)降低CFD技術(shù)的使用門檻。有限元法相對(duì)其它數(shù)值方法具有對(duì)網(wǎng)格劃分要求低、易編程的特點(diǎn)，許多偏微分方程都可通過(guò)標(biāo)準(zhǔn)的伽遼金法拆分成代數(shù)方程組，這些特點(diǎn)使其成為自編程研究采空區(qū)漏風(fēng)問(wèn)題的主要數(shù)值方法：章夢(mèng)濤、趙陽(yáng)升等給出了采空區(qū)線性滲流模型的有限元求解方法[9-10];丁廣嚷、李宗翔等研究了采空區(qū)非線性滲流的有限元求解模型[11-12]，基于Matlab開發(fā)了可以模擬采空區(qū)瓦斯涌出與耗氧情況[13]的G3程序。采空區(qū)漏風(fēng)引起的瓦斯運(yùn)移與遺煤耗氧驅(qū)動(dòng)的氧濃度場(chǎng)變化都屬于“對(duì)流-擴(kuò)散”過(guò)程。在對(duì)流起主導(dǎo)作用的采空區(qū)流體輸運(yùn)過(guò)程中，基于標(biāo)準(zhǔn)伽遼金法導(dǎo)出的有限元控制方程給出的是一個(gè)不正確的振蕩解，Zienkiewicz，章本照等的研究表明這個(gè)問(wèn)題可以通過(guò)迎風(fēng)有限元法來(lái)解決[14-15]。盡管采空區(qū)中瓦斯與氧氣的運(yùn)移過(guò)程都是“對(duì)流-擴(kuò)散”問(wèn)題，但在數(shù)學(xué)模型描述中具有不同的源項(xiàng)表現(xiàn)，采用普通對(duì)流擴(kuò)散模型的控制方程進(jìn)行求解迭代量大，并不高效。文中基于迎風(fēng)有限元法思想，導(dǎo)出采空區(qū)“瓦斯-氧”對(duì)流擴(kuò)散的統(tǒng)一模型，通過(guò)自編程將其用于煤自燃三帶與可爆瓦斯?jié)舛雀采w區(qū)域的判定。

4應(yīng)用實(shí)例

4.1模型實(shí)現(xiàn)與算例概況

文中基于VC++與ObjectARX[20]技術(shù)將以上模型及前后處理過(guò)程集成于AutoCAD中，方程組的求解采用PARDISO[21-22]并行求解器。應(yīng)用實(shí)例為某礦U形工作面通風(fēng)系統(tǒng)采空區(qū)，寬180 m，煤層開采厚度6 m，工作面設(shè)計(jì)風(fēng)量為25 m3/s，入風(fēng)隅角與回風(fēng)流隅角的相對(duì)壓力分別為308.1和268.2 Pa;工作面推進(jìn)度為3 m/d，圖2為回采90 d后形成的270 m長(zhǎng)采空區(qū)（按滲透系數(shù)渲染）。模擬中涉及的其它參數(shù)取值見(jiàn)表2.

4.2計(jì)算結(jié)果

圖3采用流網(wǎng)顯示采空區(qū)流場(chǎng)的模擬結(jié)果，其中流線的流函數(shù)值間距為1.9×10-3，等壓線間距為0.8 Pa，流網(wǎng)顏色按流速大小渲染，綠、黃、藍(lán)分別表示按流速劃定的采空區(qū)自燃三帶：散熱帶（流速大于0.004 m/s），氧化升溫帶（流速位于0.001 7至0.004 m/s）和窒息帶（流速小于0.001 7 m/s）。圖4（a）為模擬得到的采空區(qū)的氧濃度，采用氧濃度劃分三帶時(shí)，通常將濃度介于8%至18%的區(qū)域劃定了氧化升溫帶，如圖中黃色區(qū)域。由于風(fēng)速過(guò)大將使熱量無(wú)法積聚，因此更合理的方法是用風(fēng)速與氧濃度雙指標(biāo)進(jìn)行劃分，即氧化升溫帶應(yīng)滿足

模擬結(jié)果中三帶的分布及自燃帶的形狀與文獻(xiàn)[23]采用FLUENT模擬的結(jié)果相似：采空區(qū)自燃帶位于工作面后一定位置，其寬度沿進(jìn)風(fēng)側(cè)向回風(fēng)側(cè)逐漸減小。

圖5為模擬得到的采空區(qū)瓦斯?jié)舛葓?chǎng)。藍(lán)色顯示采空區(qū)深部是瓦斯富積的區(qū)域，濃度大于15%，綠色是濃度低于5%的區(qū)域，黃色區(qū)域瓦斯?jié)舛冉橛?%～15%、在正常空氣氧濃度下具有爆炸危險(xiǎn)。從圖3，圖5可知漏風(fēng)主要從2個(gè)隅角流入采空區(qū)與流回工作面，工作面與采空區(qū)邊界內(nèi)側(cè)、近上隅角的49.5 m長(zhǎng)范圍內(nèi)（回風(fēng)隅角）瓦斯?jié)舛却笥?%（5%～23%）。由于遺煤耗氧，使得采空區(qū)中氧濃度降低，不一定在任何地方都滿足瓦斯爆炸的需要。圖6左圖的瓦斯爆炸三角形給出了瓦斯與氧氣形成可爆氣體的濃度范圍[24]，當(dāng)氧濃度低于12%時(shí)，將不滿足瓦斯爆炸的需要。據(jù)此將圖4（a）與圖5疊加，得到圖6右圖紅色區(qū)域所示的具有瓦斯爆炸危階性區(qū)域，其中綠色、黃色分別表示瓦斯含量過(guò)低和過(guò)高的區(qū)域。因此在圖6右圖所示的紅色區(qū)域內(nèi)，任何點(diǎn)火源（包括煤自燃）均可引發(fā)瓦斯爆炸。因此針對(duì)模擬結(jié)果，對(duì)圖中標(biāo)定的危險(xiǎn)區(qū)域進(jìn)行惰化，以及保證必要的推進(jìn)度來(lái)防止煤自燃對(duì)避免采空區(qū)火災(zāi)或爆炸事故的發(fā)生是十分必要的。

1）基于伽遼金法建立了采空區(qū)流場(chǎng)的求解模型。給出了工作面動(dòng)態(tài)回采下采空區(qū)滲流參數(shù)的計(jì)算方法。針對(duì)采空區(qū)流場(chǎng)的定解模型，基于三角形單元離散的采空區(qū)，采用標(biāo)準(zhǔn)伽遼金法導(dǎo)出了以所有節(jié)點(diǎn)壓力為未知數(shù)的有限元控制方程組，為自編程研究采空區(qū)漏風(fēng)問(wèn)題提供了基礎(chǔ);

2）基于迎風(fēng)有限元法建立了采空區(qū)“瓦斯-氧”輸運(yùn)過(guò)程的統(tǒng)一模型。分析了采空區(qū)瓦斯涌出與遺煤耗氧的特點(diǎn)，建立了采空區(qū)中“瓦斯-氧”輸運(yùn)過(guò)程的對(duì)流擴(kuò)散模型。針對(duì)標(biāo)準(zhǔn)伽遼金法求解對(duì)流擴(kuò)散模型時(shí)解的振蕩問(wèn)題，給出了迎風(fēng)有限元法求解“瓦斯-氧”輸運(yùn)過(guò)程的控制方程組;

3）開發(fā)了采空區(qū)漏風(fēng)問(wèn)題模擬軟件。基于以上模型，以及VC++，ObjectARX等技術(shù)開發(fā)了可視化模擬軟件，應(yīng)用實(shí)例表明使用此軟件可方便地進(jìn)行采空區(qū)模型的前后處理，并根據(jù)模擬結(jié)果給出采空區(qū)煤自燃三帶及具有瓦斯爆炸危險(xiǎn)性的區(qū)域。

參考文獻(xiàn)（References）：

[1]Brune J F.Methaneair explosion hazard within coal mine gobs[J].Transactions of the Society for Mining，Metallurgy，and Exploration，2013，334：376-390.

[2]文虎，張澤，趙慶偉，等.煤層分層前后采空區(qū)自燃“三帶”的數(shù)值模擬[J].煤礦安全，2017，48（3）：178-191.

WEN Hu，ZHANG Ze，ZHAO Qingwei，et al..Numerical simulation on spontaneous combustion “Three Zones” of goaf before and after coal seam layering[J].Safety in Coal Mines，2017，48（3）：178-191.

[3]時(shí)國(guó)慶，胡方坤，王德明，等.采空區(qū)自燃“三帶”分布規(guī)律的四維動(dòng)態(tài)模擬[J].中國(guó)礦業(yè)大學(xué)學(xué)報(bào)，2014，43（2）：191-194.

SHI Guoqing，HU Fangkun，WANG Deming，et al.Unsteady simulation on distribution of three zones for spontaneous combustion in goaf areas[J].Journal of China University of Mining & Technology，2014，43（2）：191-194.

[4]YUAN Liming，Smith Alex C.Numerical study on effects of coal properties on spontaneous heating in longwall gob areas[J].Fuel，2008，87（15-16）：3409-3419.

[5]朱紅青，劉星魁.尾巷瓦斯抽采下采空區(qū)煤自燃升溫的數(shù)值模擬[J].西安科技大學(xué)學(xué)報(bào)，2012，32（1）：1-7.

ZHU Hongqing，LIU Xingkui.Numerical simulation of heating up of coal spontaneous combustion in gob area under tail roadway methane drainage pattern[J].Journal of Xi’an University of Science and Technology，2012，32（1）：1-7.

[6]XIA Tongqiang，ZHOU Fubao，WANG Xinxin，et al.Controlling factors of symbiotic disaster between coal gas and spontaneous combustion in longwall mining gobs[J].Fuel，2016，182：886-896.

[7]XIA Tongqiang，XU Mingjing，WANG Yonglong.Simulation investigation on flow behavior of gob gas by applying a newly developed FE software[J].Environmental Earth Sciences，2017，76（14）：485.

[8]高科，劉劍，鄧立軍.集于AutoCAD的采空區(qū)風(fēng)火瓦斯耦合軟件開發(fā)研究[J].中國(guó)安全生產(chǎn)科學(xué)技術(shù)，2015，11（5）：97-103.

GAO Ke，LIU Jian，DENG Lijun，et al.Development of airfiregas coupling software for goaf based on AutoCAD[J].Journal of Safety Science and Technology，2015，11（5）：97-103.

[9]章夢(mèng)濤，王景琰，梁棟.采場(chǎng)大氣中沼氣運(yùn)移過(guò)程的數(shù)值模擬[J].煤炭學(xué)報(bào)，1987，11（3）：23-30.

ZHANG Mengtao，WANG Jingyan，LIANG Dong.Numerical modelling of methane migration in air in underground workings[J].Journal of China Coal Society，1987，11（3）：23-30.

[10]趙陽(yáng)升.有限元法及其在采礦工程中的應(yīng)用[M].北京：煤炭工業(yè)出版社，1994.

[11]丁廣嚷，柏發(fā)松.采空區(qū)混合氣體運(yùn)動(dòng)基本方程及其有限元解法[J].中國(guó)礦業(yè)大學(xué)學(xué)報(bào)，1996，25（3）：21-26.

DING Guangxiang，BAI Fasong.Basic equations of mixedgas motion in gob and their finite element solution[J].Jounral of China University of Mining & Technology，1996，25（3）：21-26.

[12]李宗翔，劉玉洲，吳強(qiáng).采空區(qū)流場(chǎng)非線性滲流的改進(jìn)迭代算法[J].重慶大學(xué)學(xué)報(bào)，2008，31（2）：186-190.

LI Zongxiang，LIU Yuzhou，WU Qiang.Improved iteration algorithm for nonlinear seepage in flow field of caving goaf[J].Journal of Chongqing University，2008，31（2）：186-190.

[13]李宗翔，王曉冬，王波.采空區(qū)場(chǎng)流數(shù)值模擬程序（G3）實(shí)現(xiàn)與應(yīng)用[J].湖南科技大學(xué)學(xué)報(bào)，2005，20（3）：16-20.

LI Zongxiang，WANG Xiaodong，WANG Bo.Realization and application of numerical simulating program（G3）for field flow of goaf[J].Journal of Hunan University of Science & Technology，2005，20（3）：16-20.

[14]Zienkiewicz O C，Taylor R L，Nithiarasu P.The finite element method for fluid dynamics（7th edition）[M].Waltham：ButterworthHeinemann，2014.

[15]章本照.流體力學(xué)中的有限元方法[M].北京：機(jī)械工業(yè)出版社，1986.

[16]李宗翔.高瓦斯易自燃采空區(qū)瓦斯與自燃耦合研究[D].阜新：遼寧工程技術(shù)大學(xué)，2007.

LI Zongxiang.Study of coupling of gas and spontaneous combustion in highly gassy and spontaneous combustion goafs[D].Fuxin：Liaoning Technical University，2007.

[17]Bear J.Dynamics of fluids in porous media[M].New York，USA：America Elsevier Publishing Company Inc.，1972.

[18]秦躍平，劉偉，楊小彬.基于非達(dá)西滲流的采空區(qū)自然發(fā)火數(shù)值模擬[J].煤炭學(xué)報(bào)，2012，37（7）：1177-1183.

QIN Yueping，LIU Wei，YANG Xiaobin.Numerical simulation of impact of nonDarcy seepage on spontaneous combustion in goaf[J].Journal of China Coal Society，2012，37（7）：1177-1183.

[19]Heinrich J C，Huyakorn P S，Zienkiewicz O C，et al.An upwind finite element scheme for twodimensional convective transport equation[J].International Journal for Numerical Methods in Engineering，1977，11（1）：131-143.

[20]吳奉亮，周澎，李暉，等.基于智能對(duì)象的通風(fēng)CAD模型研究[J].煤炭科學(xué)技術(shù)，2009，37（5）：54-57.

WU Fengliang，ZHOU Peng，LI Hui，et al.Research on the model of mine ventilation CAD based on the intelligent object[J].Coal Science and Technology，2009，37（5）：54-57.

[21]Schenk，Grtner K.Solving unsymmetric sparse systems of linear equations with PARDISO[J].Journal of Future Generation Computer Systems，2004，20（3）：475-487.

[22]曹大志，強(qiáng)洪夫，任革學(xué).基于OpenMP和Pardiso的柔性多體系統(tǒng)動(dòng)力學(xué)并行計(jì)算[J].清華大學(xué)學(xué)報(bào)（自然科學(xué)版），2012，52（11）：1643-1649.

CAO Dazhi，QIANG Hongfu，REN Gexue.Parallel computing studies of flexible multibody system dynamics using OpenMP and Pardiso[J].Journal of Tsinghua University（Sci. & Tech.），2012，52（11）：1643-1649.

[23]張辛亥，席光，徐精彩，等.基于流場(chǎng)模擬的綜放面自燃危險(xiǎn)區(qū)域劃分及預(yù)測(cè)[J].北京科技大學(xué)學(xué)報(bào)，2005，27（6）：641-644.

ZHANG Xinhai，XI Guang，XU Jingcai，et al.Partition of spontaneous combustion dangerous zone and prediction of spontaneous combustion based on numerically modeling the flow field at fully mechanized caving face[J].Journal of University of Science and Technology Beijing，2005，27（6）：641-644.

[24]Coward H F，Jones G W.Limits of flammability of gases and vapors[J].Journal of the tranklin Institute，1928，203（1）：161.

[25]CHENG Jianwei，YI Luo.Modeling atmosphere composition and determining explosibility in a sealed coal mine volume[J].Achieves of Mining Science，2014，59（1）：25-40.