松散煤體內(nèi)液N2/CO2注入過(guò)程的降溫特性對(duì)比實(shí)驗(yàn)

于志金，張志鵬，楊 淞，谷 雨

(西安科技大學(xué) 安全科學(xué)與工程學(xué)院，陜西 西安 710054)

0 引 言

近年來(lái)，隨著礦井開(kāi)采深度逐漸向深部資源延伸以及綜放開(kāi)采技術(shù)的推廣應(yīng)用，煤自燃發(fā)火危險(xiǎn)性不斷提升的同時(shí)，其防治難度也在大幅增加[1]。

針對(duì)煤自燃火災(zāi)具有火源位置隱蔽、高溫范圍大、易復(fù)燃等特點(diǎn)[2]，其防治的基礎(chǔ)和關(guān)鍵在于針對(duì)性降低煤溫和大面積阻斷氧氣與煤體接觸[3]。作為目前有效的大范圍惰化降溫技術(shù)，壓注液氮(liquid nitrogen，LN2)與液態(tài)二氧化碳(liquid carbon dioxide，LCO2)被廣泛應(yīng)用于采空區(qū)煤自燃防治中。SHI和ZHOU等開(kāi)發(fā)了用于地下煤火的LN2控制技術(shù)和設(shè)備，成功控制了采空區(qū)煤的自燃[4-6]。QIN等介紹了LN2在大面積采空區(qū)火災(zāi)中的應(yīng)用，評(píng)估了LN2注入對(duì)滅火的影響[7]。張春華等提出了一種LCO2滅火系統(tǒng)，將LCO2通過(guò)管道直接注入高溫區(qū)域，實(shí)現(xiàn)了對(duì)火區(qū)的快速降溫[8]。LANG等通過(guò)實(shí)驗(yàn)測(cè)量了多孔介質(zhì)內(nèi)LN2注入過(guò)程的溫度變化，研究了LN2的傳熱和傳質(zhì)過(guò)程[9]。LIU等指出，松散煤體內(nèi)LCO2的汽化和干冰的升華均對(duì)溫度有重要影響[10]。YU和SHI等通過(guò)實(shí)驗(yàn)研究了松散煤體內(nèi)LCO2以及LN2的降溫特征與有關(guān)參數(shù)之間的關(guān)系[11-12]。更多研究則集中在數(shù)值模擬方面，王繼仁、郝朝瑜等采用數(shù)值模擬的方法結(jié)合具體工作面采空區(qū)的冒落情況，研究了不同注入條件時(shí)，采空區(qū)氧化帶范圍和溫度場(chǎng)分布，討論了低溫CO2注入位置的合理范圍[13-14]。鄧軍等探討了楊村礦采用注氮方式預(yù)防煤自燃的可行性，并采用數(shù)值模擬確定了楊村礦采空區(qū)注氮參數(shù)以及注氮后采空區(qū)自燃危險(xiǎn)區(qū)域[15]。王國(guó)旗等則對(duì)綜放工作面CO2最佳釋放位置進(jìn)行了研究，并提出了注CO2最小流量的計(jì)算方法[16]。邵昊等通過(guò)Fluent軟件研究了二氧化碳和氮?dú)庠诓煽諈^(qū)中的流動(dòng)特點(diǎn)，認(rèn)為向采空區(qū)注二氧化碳的防滅火效果比注氮?dú)庖肹17]。這些研究表明，松散煤體內(nèi)壓注LN2和LCO2均表現(xiàn)出了極強(qiáng)的冷卻和惰化作用[18-19]，但大多數(shù)研究的側(cè)重點(diǎn)均在于對(duì)關(guān)鍵參數(shù)的確定上，關(guān)于LN2和LCO2在煤體內(nèi)的傳熱特征以及其冷卻能力之間的差異仍少見(jiàn)探究。此外，大多數(shù)研究基于CO2對(duì)煤的吸附性以及阻燃性更強(qiáng)，認(rèn)為與N2相比，CO2具有更好的冷卻和惰化效果[20-23]。但這些結(jié)論主要是基于氣態(tài)介質(zhì)的研究結(jié)果，少有通過(guò)實(shí)驗(yàn)對(duì)比分析LN2和LCO2注入下的煤體降溫特征。

鑒于此，筆者通過(guò)自建實(shí)驗(yàn)裝置對(duì)煤體在注入LN2/LCO2期間和之后的溫度變化進(jìn)行了研究，分析兩者傳熱特性和冷卻機(jī)理，為現(xiàn)場(chǎng)科學(xué)設(shè)計(jì)工藝參數(shù)和效果預(yù)測(cè)提供依據(jù)與指導(dǎo)。

1 實(shí)驗(yàn)系統(tǒng)

1.1 實(shí)驗(yàn)裝置

模擬煤體內(nèi)壓注LCO2/LN2的條件和環(huán)境，自制了松散煤體LCO2/LN2壓注實(shí)驗(yàn)裝置平臺(tái)，如圖1所示。該平臺(tái)由杜瓦瓶、電子平臺(tái)秤、緩沖罐、測(cè)試容器以及管路閥門(mén)與數(shù)據(jù)采集設(shè)備等組成。

其中緩沖罐與內(nèi)徑為1 cm的壓注管路連接，管路另一端位于測(cè)試容器內(nèi)部，末端2 cm處設(shè)有截止閥和溫度、壓力傳感器，用于監(jiān)測(cè)注入介質(zhì)時(shí)的狀態(tài)和控制壓注流量，電子平臺(tái)秤用于計(jì)算平均質(zhì)量流量。測(cè)試容器材質(zhì)為不銹鋼，其內(nèi)徑60 cm，高100 cm，在豎直平面上設(shè)置了5個(gè)測(cè)試平面，分別位于距底部10(Z1)，20(Z2)，30(Z3)，40(Z4)，70 cm(Z5)，平面上共設(shè)32個(gè)Pt100溫度傳感器動(dòng)態(tài)測(cè)定容器內(nèi)的煤體瞬態(tài)溫度，管道末端的注入口與Z2平面平行，出口設(shè)置于容器頂部。

1.2 實(shí)驗(yàn)條件及方案

1)準(zhǔn)備階段，將粒徑小于10 mm的煤樣裝入測(cè)試容器，利用加熱棒對(duì)煤樣品進(jìn)行預(yù)熱之后，開(kāi)始進(jìn)行實(shí)驗(yàn)，表1和表2列出了煤樣條件及實(shí)驗(yàn)初始條件。

表1 煤樣條件Table 1 Conditions of coal

表2 實(shí)驗(yàn)條件Table 2 Conditions of test coolants

2)關(guān)閉放空閥和緩沖罐底部的閥門(mén)，打開(kāi)杜瓦瓶上的氣相閥，將緩沖罐的壓力調(diào)節(jié)至大于0.5 MPa，防止液態(tài)CO2在進(jìn)入緩沖罐后因相變產(chǎn)生干冰。

3)關(guān)閉杜瓦瓶上的氣相閥并打開(kāi)液相閥，將低溫介質(zhì)從杜瓦瓶注入到緩沖罐中。同時(shí)，通過(guò)電子平臺(tái)秤讀取杜瓦瓶的質(zhì)量變化，待緩沖罐內(nèi)壓力升至實(shí)驗(yàn)指定值后，關(guān)閉杜瓦瓶上的液相出口。

4)打開(kāi)緩沖罐底部閥門(mén)，同時(shí)保持壓注口位置旋擰閥關(guān)閉為壓注管路充壓，防止壓注初期因CO2相變過(guò)于劇烈而產(chǎn)生干冰。

5)待壓力達(dá)到實(shí)驗(yàn)預(yù)定壓力并穩(wěn)定后，打開(kāi)壓注口閥門(mén)，開(kāi)始實(shí)驗(yàn)，同時(shí)秒表記錄實(shí)驗(yàn)時(shí)間，采集各測(cè)點(diǎn)溫度數(shù)據(jù)并記錄，到達(dá)預(yù)定時(shí)間后，關(guān)閉釋放口位置旋擰閥，結(jié)束壓注。持續(xù)采集溫度數(shù)據(jù)20 h，以進(jìn)一步分析壓注結(jié)束后的傳熱特征。

2 理論分析

LN2/LCO2注入煤體的傳熱主要涉及2個(gè)方面：顯熱傳遞和潛熱傳遞，分別由等式(1)和等式(2)確定[24]。

Φs=McCp(T2-T1)

(1)

Φl=McΔh

(2)

式中Φs為顯熱熱容，J；Mc為煤體質(zhì)量，kg；Cp為比熱容，J·(kg·K)-1；T1和T2分別為測(cè)試之前和之后的介質(zhì)溫度，K；Φl為潛熱熱容，J；Δh為汽化潛熱，J·kg-1。通過(guò)上述方程即可定義壓注過(guò)程中的總熱容Φ。

Φ=ΦS+Φl

(3)

由于LN2/LCO2離開(kāi)壓注口后會(huì)進(jìn)行連續(xù)的相變，導(dǎo)致Cp和Δh的參數(shù)隨時(shí)變化，無(wú)法僅通過(guò)溫度數(shù)據(jù)進(jìn)行計(jì)算?；诖耍跍囟葓?chǎng)中定義溫度小于0 ℃的區(qū)域?yàn)槔鋮s區(qū)域，通過(guò)該區(qū)域面積的擴(kuò)大率ve來(lái)分析壓注過(guò)程中的冷卻效果，定義為

ve=ΔL/Δt=(Lt2-Lt1)/t2-t1

(4)

式中L為在對(duì)應(yīng)壓注時(shí)間下的最大冷卻距離，cm；t為注入時(shí)間，s。

停注階段，受溫度梯度影響的傳熱過(guò)程可以根據(jù)一維傅立葉傳導(dǎo)定律[25]計(jì)算，表示為

(5)

式中ke通過(guò)等式(6)[26]計(jì)算得出

ke=(1-ε)ks+εkg

(6)

式中Rh為傳熱率，W·m-2；Tx為選定位置的溫度，K；x1-x2為對(duì)應(yīng)溫度傳感器之間的距離，m；ε為煤樣空隙率，實(shí)驗(yàn)條件已給出；ks為煤的導(dǎo)熱系數(shù)，0.12 W·(m·K)-1;kg為混合流體在平均溫度下的導(dǎo)熱系數(shù)。

3 實(shí)驗(yàn)結(jié)果對(duì)比與分析

3.1 壓注過(guò)程

為了直觀(guān)地研究壓注過(guò)程中煤體的溫度分布，將得到的溫度數(shù)據(jù)導(dǎo)入Origin軟件，繪制得到溫度分布云圖。

如圖2所示，LN2/LCO2注入煤體后，在壓差和射流動(dòng)能的作用下高速滲流擴(kuò)散，在節(jié)流膨脹和閃蒸作用下迅速釋放大量的相變潛熱，形成了一個(gè)具有強(qiáng)制冷效果的冷卻區(qū)域，最低溫度達(dá)到大氣壓下的沸騰溫度(LN2為-195.8℃，LCO2為-78.4 ℃)。隨著壓注時(shí)間的增加，2種冷卻劑形成的冷卻區(qū)域均隨之?dāng)U大蔓延，其冷卻能力隨著與壓注口的距離增加而逐漸降低。對(duì)比圖2(a)和(b)可知，LN2注入下，所形成的低溫區(qū)域沿流體滲流路徑緩慢擴(kuò)大，而LCO2的注入產(chǎn)生了一個(gè)明顯的核心冷卻區(qū)域，并且在垂直和水平方向上呈現(xiàn)均勻擴(kuò)張的特點(diǎn)。

為了研究?jī)烧呃鋮s能力的差異，在圖3和圖4中列出了水平和垂直平面中選定位置的溫度和冷卻速率與時(shí)間的關(guān)系。如圖4所示，冷卻劑為L(zhǎng)N2時(shí)，除4#迅速降至沸點(diǎn)溫度以外，其它觀(guān)測(cè)點(diǎn)的溫度都呈線(xiàn)性下降，并且垂直平面上不同測(cè)點(diǎn)溫度變化趨勢(shì)幾乎一致，結(jié)合圖2中LN2的低溫區(qū)域沿滲流路徑緩慢擴(kuò)大的現(xiàn)象，表明LN2對(duì)煤體的冷卻能力主要取決于自身極低的溫度形成的劇烈溫差所造成的對(duì)流換熱。對(duì)于LCO2，溫度曲線(xiàn)則更為復(fù)雜，由于大量相變潛熱的釋放導(dǎo)致了干冰的形成，結(jié)合圖2中的核心冷卻區(qū)域可以發(fā)現(xiàn)，LCO2的冷卻機(jī)理主要由2個(gè)過(guò)程控制：相態(tài)轉(zhuǎn)變(焦耳-湯姆遜效應(yīng))而產(chǎn)生的冷卻效應(yīng)以及對(duì)流傳導(dǎo)引起的熱傳遞，相變潛熱所形成的核心冷卻區(qū)域溫降幅度比對(duì)流傳導(dǎo)區(qū)域更為顯著。

圖2 壓注過(guò)程中不同時(shí)刻煤體溫度分布Fig.2 Temperature profiles of coal at different times during the injection

圖3 壓注過(guò)程Z2平面不同測(cè)點(diǎn)的溫度以及冷卻速率與時(shí)間關(guān)系Fig.3 Temperature and cooling rate versus time for selected locations in the Z2 plane during the injection

圖4 壓注過(guò)程垂直平面不同測(cè)點(diǎn)的溫度以及冷卻速率與時(shí)間關(guān)系Fig.4 Temperature and average cooling rate versus time for selected locations in the vertical plane during the injection

為表征注入過(guò)程中LCO2/LN2對(duì)煤體的冷卻效果，以不同時(shí)刻的0 ℃等溫線(xiàn)定義冷卻區(qū)域的變化。

如圖5所示，與LN2相比，LCO2的冷卻區(qū)域在480 s內(nèi)明顯更大，兩者ve的最大偏差出現(xiàn)在前20 s，這說(shuō)明在壓注初始階段，LCO2的冷卻速率比LN2高，相變潛熱釋放更快，但隨著注入時(shí)間的增加，LN2的ve逐漸超過(guò)LCO2，這意味著與LCO2幾乎不變的核心冷卻區(qū)域相比，LN2的冷卻區(qū)域會(huì)以較小的速率持續(xù)膨脹。

圖5 壓注過(guò)程冷卻區(qū)域擴(kuò)大率Fig.5 Evolution of cooling area(T<0 ℃)during the injection

3.2 停注階段

圖6展示了停注階段20 h內(nèi)的溫度場(chǎng)變化。停注1 h后，壓注LN2/LCO2引起的主要冷卻區(qū)域趨于消失，但壓注過(guò)程中溫度沒(méi)有明顯變化的測(cè)點(diǎn)在該階段有明顯降溫趨勢(shì)，這意味著停注后冷卻能力的差異主要?dú)w因于初始溫度梯度。

圖6 停注后不同時(shí)刻煤體溫度分布Fig.6 Temperature profiles of coal at different times after the injection

圖7和圖8分別列出了水平和垂直平面中選定位置的溫度以及溫度回升率的隨時(shí)間的關(guān)系。停注后，LN2所產(chǎn)生的冷卻能量以較快的速度損耗，煤體溫度回升幅度劇烈，6 h后大部分區(qū)域溫度已回升至常溫水準(zhǔn)；而LCO2由于較小的熱導(dǎo)率和較高的汽化熱，溫度以較為穩(wěn)定的趨勢(shì)回升，特別是核心區(qū)域的溫度在停注6小時(shí)后仍維持在0 ℃以下，因此在持續(xù)冷卻方面LCO2更具競(jìng)爭(zhēng)力。為研究溫度梯度對(duì)溫度回升過(guò)程的影響，將注入方向相鄰觀(guān)測(cè)點(diǎn)之間的傳熱率作為比較指標(biāo)，表征換熱強(qiáng)度隨時(shí)間的變化規(guī)律。

如圖9所示，對(duì)于LN2，R4-5在停注1 h內(nèi)劇烈下降而后保持恒定，R5-6和R6-7則隨時(shí)間持續(xù)不斷降低，這意味著停注后的傳熱率主要由較低的沸騰溫度和相變所引起的溫度梯度所決定，因此LN2停注后的冷卻區(qū)域?qū)⒁宰⑷脒^(guò)程中產(chǎn)生的核心冷卻區(qū)域?yàn)榛A(chǔ)向各個(gè)方向發(fā)展。而LCO2的傳熱率變化則更具爭(zhēng)議性，盡管在R5-6的變化趨勢(shì)上LCO2和LN2呈現(xiàn)一致，但在R4-5和R5-6，R5-6與R6-7之間LCO2的傳熱方向相反。結(jié)合圖7中LCO2的溫度變化規(guī)律可知，這歸因于LCO2注入過(guò)程中產(chǎn)生的冷卻能量大多集中在4#，5#觀(guān)測(cè)點(diǎn)的區(qū)域，另外由于壓注過(guò)程中的6#，7#溫度沒(méi)有明顯變化，導(dǎo)致R6-7的傳熱方向主要由初始溫度決定。

圖7 停注后Z2平面不同測(cè)點(diǎn)的溫度以及溫度回升率與時(shí)間關(guān)系Fig.7 Temperature and temperature recovery rate versus time for selected locations in the Z2 plane after the injection

圖8 停注后垂直平面不同測(cè)點(diǎn)的溫度以及冷卻速率與時(shí)間關(guān)系Fig.8 Temperature and temperature recovery rate versus time for selected locations in the vertical plane after the injection

圖9 停注后不同位置傳熱率與時(shí)間關(guān)系Fig.9 Heat transfer rate versus the time for different locations after injection

4 結(jié) 論

1)隨著注入時(shí)間的增加，LN2的核心冷卻區(qū)域范圍仍以一定速率提高；而對(duì)于LCO2，一旦相變完成，核心冷卻區(qū)域幾乎保持不變。

2)LN2的冷卻機(jī)理主要取決于LN2與煤體間的劇烈溫差導(dǎo)致的熱傳導(dǎo)，相比之下，相變潛熱和對(duì)流傳熱作用是影響LCO2的冷卻效果的關(guān)鍵。

3)停注后，LN2和LCO2的持續(xù)冷卻能力的差異主要?dú)w因于初始溫度梯度的不同。由于較小的熱導(dǎo)率和較高的汽化熱，停注后液態(tài)CO2的持續(xù)冷卻能力優(yōu)于液態(tài)N2。