王 良

(同煤集團(tuán)朔州煤電公司，山西 朔州 038300)

1 研究背景

本文提出了利用火星溫度預(yù)測(cè)模型預(yù)測(cè)煤地下氣化過(guò)程中最關(guān)鍵的參數(shù)溫度，即解決在活動(dòng)礦山地下開采過(guò)程對(duì)溫度的控制。因?yàn)樘烊粴庖夯且粋€(gè)強(qiáng)放熱過(guò)程，提高合成氣的質(zhì)量，需要較高的溫度。Hamanaka等[1]也表明,在煤炭地下氣化效率或能量回收率過(guò)程中可以通過(guò)增加反應(yīng)溫度和擴(kuò)大氣化面積增加輸出質(zhì)量。因此,保持較高的溫度且沒(méi)有超過(guò)安全限度避免火災(zāi)的風(fēng)險(xiǎn)是一項(xiàng)至關(guān)重要的技術(shù)。

2 數(shù)據(jù)來(lái)源

針對(duì)5.4 m厚煤層(位于大約460 m深)進(jìn)行試驗(yàn),氣化試驗(yàn)持續(xù)了60 d,理論最大的煤炭氣化能力600 kg/h,相當(dāng)于3 MW氣體焓值。表1給出了風(fēng)干煤煤層特點(diǎn)。

為了測(cè)試幾種替代方法并探索如何提高生產(chǎn)率，實(shí)驗(yàn)采用了不同的氣化劑:空氣、空氣和氧氣、空氣和CO、空氣和氮?dú)饧翱諝夂退?/p>

試驗(yàn)的主要階段為冷卻前的1 343 h，其主要階段為：

a) 第一階段:用煙火藥點(diǎn)火后，煤通過(guò)富氧空氣點(diǎn)火。這個(gè)階段持續(xù)了193 h。

b) 第二階段:氣化劑為空氣。在這一階段，由于一個(gè)分離器的冷凝水平過(guò)高，出現(xiàn)了技術(shù)故障。通過(guò)增加空氣量來(lái)重新開始操作。這個(gè)階段持續(xù)了695 h。

表1 風(fēng)干煤煤層特點(diǎn)

c) 第三階段:氣化劑為空氣和二氧化碳，測(cè)試二氧化碳的行為。這個(gè)階段持續(xù)了120 h。

d) 第四階段:氣化劑為空氣，在二氧化碳階段后對(duì)反應(yīng)器進(jìn)行再加熱。這個(gè)階段持續(xù)了173 h。

e) 第五階段:氣化劑為空氣和氮?dú)?，通過(guò)逐漸增加氮?dú)獗壤龁?dòng)實(shí)驗(yàn)滅火。這個(gè)階段持續(xù)了162 h。

f) 第六階段:氣化劑為氮?dú)?，用于滅火和冷卻試驗(yàn)。該階段持續(xù)時(shí)間超過(guò)500 h，不被認(rèn)為是實(shí)驗(yàn)的特定部分。

由于實(shí)驗(yàn)是在一個(gè)活躍的礦井中進(jìn)行的，所以它的設(shè)計(jì)要滿足所有的采礦安全要求。此外，國(guó)家也對(duì)地質(zhì)反應(yīng)器出口氣體的最高溫度進(jìn)行規(guī)定[2]。為了防止發(fā)生火災(zāi)，溫度被設(shè)定在550 ℃。此外，當(dāng)溫度在600 ℃以上時(shí)，如果工藝氣體泄漏到通道中，也有自燃的可能。

3 模型來(lái)源

火星模型最初是由Ref[3]引入的。這類模型屬于一組多元非參數(shù)回歸技術(shù)，這些技術(shù)根據(jù)具體數(shù)據(jù)調(diào)整其參數(shù)。

采用分段多邊形方法建立了火星模型。

通過(guò)將基函數(shù)擬合到X的不同區(qū)間。多項(xiàng)式的連接點(diǎn)用t表示，稱為節(jié)點(diǎn)、擊穿點(diǎn)或節(jié)。

為了確定應(yīng)該參與模型的基本函數(shù)或它們的估計(jì)重要性，MARS使用了通用交叉驗(yàn)證(GCV)。通過(guò)計(jì)算每個(gè)變量的模型子集n的數(shù)目，計(jì)算剩余平方和(RSS)的均值，除以一個(gè)懲罰，這個(gè)懲罰取決于模型的復(fù)雜性[4-5]:

用于開發(fā)火星模型來(lái)預(yù)測(cè)本次實(shí)驗(yàn)溫度的工具是R語(yǔ)言，使用EARTH工具包，它是由MDA衍生而來(lái)。為了開發(fā)圖形，使用了GGPLOT2工具包。

4 實(shí)驗(yàn)結(jié)果

預(yù)測(cè)模型對(duì)溫度值的魯棒性是通過(guò)對(duì)模型進(jìn)行5次獨(dú)立的再訓(xùn)練來(lái)確定的，每次都使用數(shù)據(jù)集的80%作為訓(xùn)練子集，并對(duì)剩下的20%進(jìn)行驗(yàn)證。

隨后，在實(shí)驗(yàn)的各個(gè)階段建立火星模型，分析影響預(yù)測(cè)精度的因素和原因。

盡管火星模型建立了僅使用的數(shù)據(jù)從第一階段氣化劑使用富氧空氣,該模型能夠預(yù)測(cè)溫度的值提前1 h的差異小于11 ℃在95%的情況(見表2)。這個(gè)階段持續(xù)了193 h。

表2 實(shí)驗(yàn)各階段實(shí)際值與預(yù)測(cè)值之間的絕對(duì)差的百分比

表中95%百分比值，即為分析中使用的值，比用整個(gè)實(shí)驗(yàn)的數(shù)據(jù)所作的預(yù)測(cè)更好。圖1為實(shí)驗(yàn)第一階段所用氣化劑的溫度和用量。

然而第二階段時(shí),氣化劑使用的是空氣,能夠獲得的模型計(jì)算溫度提前一小時(shí)只有的差異小于19 ℃在95%的情況下(見表2)。這個(gè)階段持續(xù)了695 h,因此在實(shí)驗(yàn)中是最長(zhǎng)的階段。

這種較低的精度可能與此階段發(fā)生的技術(shù)故障有關(guān)，這是由于一個(gè)分離器的冷凝水平過(guò)高，導(dǎo)致實(shí)驗(yàn)溫度范圍更大。

圖1 實(shí)驗(yàn)第一階段的實(shí)際溫度和氣化劑用量

然而，與實(shí)驗(yàn)其他階段的模型相比，第二階段模型的準(zhǔn)確性是最好的。圖2為同期真實(shí)值與預(yù)測(cè)值的溫度散點(diǎn)圖。

圖2 同一時(shí)間對(duì)應(yīng)的實(shí)驗(yàn)第二階段的實(shí)際(Temp_out)與預(yù)測(cè)值(Temp_outP1)的溫度散點(diǎn)圖

第三階段，持續(xù)120 h，使用空氣和二氧化碳作為氣化劑。在這種情況下，得到的模型能夠在95%的情況下提前1 h計(jì)算出溫度，且差異小于5 ℃(見表2)。圖3為實(shí)驗(yàn)第三階段所用氣化劑的溫度和用量。第66頁(yè)圖4給出了與圖2相同尺度下，這一階段的溫度真實(shí)值散點(diǎn)圖與同期的預(yù)測(cè)值，以解釋相對(duì)于第二階段的精度差異。

圖3 實(shí)際溫度和氣化劑在實(shí)驗(yàn)三階段使用的數(shù)量

第四階段，在二氧化碳階段結(jié)束后，以空氣為氣化劑對(duì)反應(yīng)器進(jìn)行再加熱，持續(xù)173 h，所得到的模型能夠在95%的情況下提前1 h計(jì)算出溫差小于8 ℃的溫度(見表2)。在這種情況下，預(yù)測(cè)比使用整個(gè)實(shí)驗(yàn)數(shù)據(jù)時(shí)要好，但差異很小。這可能是由于，使用空氣作為唯一的氣化劑，類似于第二階段，但沒(méi)有使用其他氣化劑的特殊情況下，沒(méi)有技術(shù)停機(jī)時(shí)間。

圖4 同一時(shí)間對(duì)應(yīng)的實(shí)驗(yàn)第三階段的真實(shí)值(Temp_out)與預(yù)測(cè)值(Temp_outP1)的溫度散點(diǎn)圖

第五階段通過(guò)逐漸增加氮?dú)獗壤龁?dòng)實(shí)驗(yàn)滅火。在此情況下，得到的模型能夠在95%的情況下提前1 h計(jì)算出溫差小于11 ℃(見表2)。

此外，在這種情況下，由于參數(shù)之間的關(guān)系更穩(wěn)定，只有空氣和N2，所以比使用整個(gè)實(shí)驗(yàn)的數(shù)據(jù)得到了更好的預(yù)測(cè)使用。

第六階段僅使用氣化劑Nas對(duì)實(shí)驗(yàn)進(jìn)行了熄滅和冷卻。在這種情況下，所得到的模型不能用所有的實(shí)驗(yàn)數(shù)據(jù)來(lái)改善溫度預(yù)測(cè)。在95%的情況下，提前1 h估計(jì)了溫度，誤差小于13 ℃(見表2)。

5 結(jié)論

本研究的重點(diǎn)是使用多變量自適應(yīng)回歸樣條(MARS)方法，煤礦試點(diǎn)開發(fā)的實(shí)驗(yàn)數(shù)據(jù)，利用UCG過(guò)程中的溫度預(yù)測(cè)模型。實(shí)現(xiàn)了活動(dòng)礦井防止火災(zāi)風(fēng)險(xiǎn)，同時(shí)在煤層氣開采過(guò)程中獲得最佳的合成氣質(zhì)量。

為了描述這個(gè)問(wèn)題，輸入變量以小時(shí)為單位進(jìn)行測(cè)量:不同氣化劑的用量、合成氣的總流量和組成、溫度和熱值。輸出變量是一小時(shí)前的合成氣溫度。

在95%的情況下，預(yù)測(cè)的合成氣溫度差小于15 ℃。該模型在經(jīng)過(guò)5次再訓(xùn)練后證明了一致性，使用80%的數(shù)據(jù)集作為訓(xùn)練子集，其余20%用于驗(yàn)證。

在實(shí)驗(yàn)的不同階段在較低的溫度范圍和較穩(wěn)定的氣化劑使用情況下，可以獲得較好的預(yù)測(cè)結(jié)果?？偠灾?，火星模型能夠高精度地預(yù)測(cè)氣化過(guò)程參數(shù)之間的相互作用，即使只有少量的觀測(cè)，并且在關(guān)于氣化劑的數(shù)量和組成不斷變化的環(huán)境中也能做到。在預(yù)期的一小時(shí)內(nèi)預(yù)測(cè)合成氣的溫度，將在活動(dòng)礦井內(nèi)進(jìn)行地下天然氣的開采，在獲得最佳合成氣質(zhì)量的同時(shí)，充分防止火災(zāi)風(fēng)險(xiǎn)。