吳士龍

（河南理工大學(xué)能源科學(xué)與工程學(xué)院，河南焦作454000）

通過(guò)氣化清潔高效地使用煤炭被認(rèn)為是確保能源市場(chǎng)對(duì)這種燃料的持續(xù)需求的一種有希望的替代方法[1]。這項(xiàng)技術(shù)的特點(diǎn)是通用性強(qiáng)，可以使用多種燃料以及對(duì)產(chǎn)品開放的各種應(yīng)用。產(chǎn)生的合成氣可用于發(fā)電、化學(xué)原料、氫氣等[2]。煤氣化是一種成熟的生產(chǎn)合成氣（H2+CO）的技術(shù)，建立煤氣化過(guò)程的代理模型對(duì)于理解操作參數(shù)對(duì)氣化性能的影響和優(yōu)化過(guò)程性能及確定操作條件是重要的。Kriging代理模型具有高度擬合的優(yōu)點(diǎn)，已經(jīng)被廣泛使用[3]。另外蒙特卡洛模擬也是進(jìn)行不確定性量化分析的一個(gè)重要方法[4]。在這項(xiàng)工作中，研究的目標(biāo)是：評(píng)估操作參數(shù)的波動(dòng)性對(duì)合成氣產(chǎn)量的影響。

1 數(shù)據(jù)和方法

1.1 煤氣化數(shù)據(jù)

在Kriging代理模型的建立過(guò)程中使用了印度的燃料研究所在文獻(xiàn)中公布的煤氣化數(shù)據(jù)以及煤的性質(zhì)[5]。煤的性質(zhì)如表1所示，煤氣化數(shù)據(jù)如表2所示。

表1 煤性質(zhì)

表2 煤氣化數(shù)據(jù)

1.2 Kriging代理模型

高壓煤氣化技術(shù)的主要優(yōu)點(diǎn)之一是產(chǎn)生的合成氣可用于多種應(yīng)用。為此，該過(guò)程可用于優(yōu)化不同的響應(yīng)變量。因此，在煤氣化過(guò)程中，最重要的變量是合成氣產(chǎn)量。同時(shí)，還必須考慮發(fā)熱量，發(fā)熱量的高低代表了煤的燃燒是否充分，發(fā)熱量越高代表煤燃燒越充分。

Kriging模型是一種插值模型，由Danie G.Krige開發(fā)，最早用于用于預(yù)測(cè)采礦孔[6]。隨后Sacks等[7]將其引入計(jì)算機(jī)實(shí)驗(yàn)。為了檢驗(yàn)代理模型對(duì)實(shí)驗(yàn)數(shù)據(jù)擬合的統(tǒng)計(jì)顯著性，通過(guò)方差分析來(lái)對(duì)其進(jìn)行檢驗(yàn)。

1.3 蒙特卡洛模擬

蒙特卡洛模擬[8，9]也稱計(jì)算機(jī)隨機(jī)模擬，是以隨機(jī)數(shù)集為基礎(chǔ)的一種概率統(tǒng)計(jì)方法。蒙特卡洛模擬作為統(tǒng)計(jì)學(xué)模型工具，可以在使用時(shí)進(jìn)行模擬仿真上萬(wàn)次，這相對(duì)于進(jìn)行試驗(yàn)或者進(jìn)行大型的計(jì)算機(jī)仿真來(lái)說(shuō)非常的便捷。可以在使用時(shí)模擬仿真數(shù)千至上萬(wàn)次，模擬一次相當(dāng)于做了一次煤氣化過(guò)程的實(shí)驗(yàn)。

2 結(jié)果與討論

2.1 代理模型

建立關(guān)于合成氣產(chǎn)量和發(fā)熱量的Kriging代理模型，我們將其代理模型通過(guò)圖1和圖2來(lái)對(duì)其進(jìn)行直觀地表示。通過(guò)計(jì)算得出關(guān)于合成氣產(chǎn)量的方差為0.215%，發(fā)熱量的方差為0.199%，表明所建立的代理模型精度能夠滿足接下來(lái)的研究。

圖1 合成氣產(chǎn)量代理模型三維圖

2.2 全局靈敏度分析

作為不確定性量化評(píng)估的一部分，基本分析之一是全局敏感性分析。它旨在回答以下問(wèn)題:對(duì)于給定的參數(shù)空間，哪些輸入因素對(duì)觀察到的感興趣的量的可變性影響最大，因此可以分配有限的資源來(lái)減少這些輸入因素的不確定性。在傳統(tǒng)的靈敏度分析中，設(shè)計(jì)空間中固定點(diǎn)的梯度(通常為平均值)用于評(píng)估單個(gè)因素的靈敏度，這提供了非常有限的靈敏度視圖。因此，全局靈敏度提供了整個(gè)設(shè)計(jì)空間中靈敏度的整體視圖，從而為設(shè)計(jì)提供了完整的覆蓋范圍。

作為不確定性量化評(píng)估的一部分，基本分析之一是全局敏感性分析。它旨在回答以下問(wèn)題:對(duì)于給定的參數(shù)空間，哪些輸入因素對(duì)觀察到的感興趣的量的可變性影響最大，因此可以分配有限的資源來(lái)減少這些輸入因素的不確定性。在傳統(tǒng)的靈敏度分析中，設(shè)計(jì)空間中固定點(diǎn)的梯度(通常為平均值)用于評(píng)估單個(gè)因素的靈敏度，這提供了非常有限的靈敏度視圖。因此，全局靈敏度提供了整個(gè)設(shè)計(jì)空間中靈敏度的整體視圖，從而為設(shè)計(jì)提供了完整的覆蓋范圍。

根據(jù)上面的數(shù)據(jù)建立了Kriging代理模型，在該代理模型上進(jìn)行了全局靈敏度分析。表3給出了全局靈敏度分析的結(jié)果，結(jié)果表明對(duì)合成氣產(chǎn)量產(chǎn)生的主要影響是操作參數(shù)的變化而不是參數(shù)之間的相互作用，其中空氣流量對(duì)合成氣質(zhì)量的影響最大，占比為63.95%，水蒸氣占比為18.41%，溫度占比為12.96%；此外，參數(shù)之間的相互作用只有空氣與水蒸氣之間的相互作用大于1%，為4.6%。

表3 參數(shù)對(duì)合成氣產(chǎn)量的貢獻(xiàn)

2.3 不確定性傳播

不確定性的正向傳播分析了輸入變量的不確定性對(duì)感興趣參數(shù)數(shù)量的影響。使用的方法為蒙特卡洛，我們此節(jié)是為了說(shuō)明參數(shù)的波動(dòng)性對(duì)合成氣產(chǎn)量的影響。值得注意的是，這些參數(shù)的波動(dòng)性是不容易測(cè)量的，因此我們只能假設(shè)參數(shù)是保持不變的或者是波動(dòng)的并賦予其一個(gè)概率密度函數(shù)(PDF)。表4、表5展示了我們所假設(shè)的幾種參數(shù)分布情況。

表4 不確定性案例

表5 天然氣產(chǎn)量均值與標(biāo)準(zhǔn)差

從表3和圖2可以看出：空氣流量對(duì)于天然氣產(chǎn)量有著重大影響，根據(jù)圖2的不確定性分析結(jié)果，在給定規(guī)定的操作參數(shù)不確定性的情況下，評(píng)估天然氣產(chǎn)量達(dá)到2.9（m3/kg煤）的概率，我們會(huì)得到不同結(jié)果。對(duì)于案例1，空氣流量的概率密度函數(shù)為均勻分布，天然氣產(chǎn)量達(dá)到2.9（m3/kg煤）的概率為0；對(duì)于案例2，空氣流量的概率密度函數(shù)是均值為2.095方差為0.2095的正態(tài)分布，天然氣產(chǎn)量達(dá)到2.9（m3/kg煤）的概率為62.14%；對(duì)于案例3，水蒸氣流量的概率密度函數(shù)是均值為0.385方差為0.0385的正態(tài)分布，天然氣產(chǎn)量達(dá)到2.9（m3/kg煤）的概率為100%；對(duì)于案例4，溫度的概率密度函數(shù)是均勻分布，天然氣產(chǎn)量達(dá)到2.9（m3/kg煤）的概率為66.22%；對(duì)于案例5，空氣流量的概率密度函數(shù)是均值為2.095方差為0.2095的正態(tài)分布，水蒸氣流量的概率密度函數(shù)是均值為0.385方差為0.0385的正態(tài)分布，溫度的概率密度函數(shù)是均勻分布，天然氣產(chǎn)量達(dá)到2.9（m3/kg煤）的概率為69.74%。

圖2 天然氣產(chǎn)量經(jīng)驗(yàn)累積分布圖

3 結(jié)論

不確定性量化在煤氣化中的應(yīng)用是通過(guò)現(xiàn)有的煤氣化數(shù)據(jù)集對(duì)實(shí)驗(yàn)中適用的條件來(lái)進(jìn)行模擬，輸入?yún)?shù)和目標(biāo)函數(shù)是固定的。我們通過(guò)原始數(shù)據(jù)建立了Kriging代理模型，在該模型的基礎(chǔ)上進(jìn)行了全局靈敏度分析，得到了空氣、水蒸氣和溫度以及它們之間的相互作用對(duì)天然氣產(chǎn)量的影響，此外進(jìn)行了參數(shù)的不確定性傳播，并用概率密度函數(shù)表示它們，結(jié)果表明參數(shù)的不同分布情況下天然氣產(chǎn)量達(dá)到某一概率密度的概率不同。