（煤基清潔能源國家重點(diǎn)實(shí)驗(yàn)室(中國華能集團(tuán)清潔能源技術(shù)研究院有限公司））

0 引言

2012 年11 月6 日，我國首座、世界第6 座大型整體煤氣化聯(lián)合循環(huán)（IGCC）電站——天津IGCC 示范電站投產(chǎn)發(fā)電[1]。與傳統(tǒng)燃煤電廠相比，IGCC 作為多種設(shè)備、多項(xiàng)技術(shù)集成（包括空分裝置、煤氣化及凈化系統(tǒng)、燃?xì)廨啓C(jī)、余熱鍋爐以及蒸汽輪機(jī)等）的一個(gè)復(fù)雜系統(tǒng)[2-3]，對(duì)于運(yùn)行人員技術(shù)知識(shí)、操作水平、事故分析和處理能力等方面的要求都很高，也對(duì)運(yùn)行人員的技術(shù)培訓(xùn)提出了更高的要求。這就需要研究開發(fā)針對(duì)IGCC的仿真系統(tǒng)，提供完善的實(shí)時(shí)培訓(xùn)環(huán)境，以便能夠可靠、有效的檢驗(yàn)學(xué)員的操作過程和操作表現(xiàn)，加快培訓(xùn)速度并節(jié)約大量的培訓(xùn)費(fèi)用，還可用于研究IGCC設(shè)計(jì)、優(yōu)化及控制技術(shù)。

天津IGCC 示范電站所用燃?xì)廨啓C(jī)為SGT5-2000E(LC)型低熱值合成氣燃?xì)廨啓C(jī)，該燃?xì)廨啓C(jī)由常規(guī)燃用天然氣的SGT5-2000E 型燃?xì)廨啓C(jī)改造而來，其額定負(fù)荷為171MW，采用單缸單軸、軸向排氣、冷端驅(qū)動(dòng)的結(jié)構(gòu)，有兩個(gè)筒形燃燒室，每個(gè)燃燒室安裝8支燃燒器，啟動(dòng)燃料為柴油，正常運(yùn)行燃料為氣化爐所產(chǎn)合成氣，燃燒器燃燒方式為雙通道擴(kuò)散燃燒[4-5]。燃?xì)廨啓C(jī)由壓氣機(jī)、燃燒室、透平等三大部件組成，其結(jié)構(gòu)示意圖可參見圖1。目前已有大量關(guān)于天然氣燃?xì)廨啓C(jī)的研究，采用數(shù)值計(jì)算[6-8]、計(jì)算流體力學(xué)模擬[9]到大渦模擬[10]，對(duì)于燃?xì)廨啓C(jī)性能預(yù)測(cè)、問題分析及改進(jìn)均達(dá)到了一定的深度，并取得了一定的成效。雖然已有采用商業(yè)軟件建立壓縮機(jī)的仿真模型[11-12]和采用神經(jīng)網(wǎng)絡(luò)算法建立軸流通風(fēng)機(jī)仿真模型[13]的研究，但對(duì)于低熱值合成氣燃?xì)廨啓C(jī)的仿真研究尚不多見。本文將針對(duì)低熱值合成氣燃?xì)廨啓C(jī)開展仿真系統(tǒng)設(shè)計(jì)及驗(yàn)證。

圖1 低熱值合成氣燃?xì)廨啓C(jī)結(jié)構(gòu)示意圖Fig.1 Structure low heating value syngas gas turbine

1 模型搭建

1.1 仿真系統(tǒng)平臺(tái)

本項(xiàng)目采用成熟的商業(yè)軟件平臺(tái)DYNSIM 來開發(fā)低熱值合成氣燃?xì)廨啓C(jī)仿真機(jī)，以便將精力集中于低熱值燃?xì)廨啓C(jī)的特有過程特點(diǎn)分析和機(jī)理建模。DYNSIM是SimSci公司開發(fā)的基于嚴(yán)格計(jì)算的、成熟的動(dòng)態(tài)過程模擬系統(tǒng)，是Pro/II 的動(dòng)態(tài)版本，采用基于機(jī)理的技術(shù)和嚴(yán)格的熱力學(xué)數(shù)據(jù)，提供準(zhǔn)確可靠的計(jì)算結(jié)果，用于解決從工程分析、控制系統(tǒng)校核到操作員培訓(xùn)系統(tǒng)等工作中遇到的動(dòng)態(tài)模擬問題。DYNSIM 目前已廣泛應(yīng)用于化工過程的動(dòng)態(tài)模擬[15-17]和人員培訓(xùn)[18]。

1.2 低熱值燃?xì)廨啓C(jī)主要模型搭建

1.2.1 壓氣機(jī)

圖2 低熱值合成氣燃?xì)廨啓C(jī)仿真模型Fig.2 Simulation model of gas turbine for low heating value syngas

在DYNSIM中利用自帶的壓縮機(jī)（Compressor）模型作為低熱值合成氣燃?xì)廨啓C(jī)壓氣機(jī)仿真模型，見圖2。壓縮后的空氣一部分與燃料（柴油或合成氣）混合燃燒，另一部分冷卻燃燒室和透平，不過最終都將進(jìn)入燃?xì)廨啓C(jī)透平而做功，所以本模塊不再單獨(dú)設(shè)置冷卻通道，將壓氣機(jī)和透平做單通道連接。為保證燃?xì)廨啓C(jī)總體質(zhì)量和能量平衡，假設(shè)壓氣機(jī)壓縮后的空氣全部進(jìn)入燃燒室燃燒，燃燒后的高溫廢氣驅(qū)動(dòng)透平做功。

1.2.2 燃燒室

為減少仿真系統(tǒng)復(fù)雜性，同時(shí)能夠準(zhǔn)確反映低熱值合成氣燃?xì)廨啓C(jī)燃燒室的運(yùn)行特性，對(duì)于燃燒過程的仿真做了大幅簡化。在DYNSIM 中利用自帶的燃燒器（Combustor）模型作為低熱值合成氣燃?xì)廨啓C(jī)燃燒室仿真模型，見圖2。在進(jìn)入燃燒室仿真模型之前，就將壓縮后的空氣和燃料進(jìn)行混合，然后進(jìn)入燃燒室仿真模型，不考慮真實(shí)燃?xì)廨啓C(jī)燃燒室中的擴(kuò)散燃燒。同時(shí)假設(shè)燃料的燃燒效率為100%，即所有可燃物完全燃燒。

1.2.3 透平

在DYNSIM 中利用自帶的膨脹機(jī)（Expander）模型作為低熱值合成氣燃?xì)廨啓C(jī)透平仿真模型，見圖2。轉(zhuǎn)動(dòng)軸仿真模型將壓氣機(jī)仿真模型及透平仿真模型連接。膨脹機(jī)的通流能力可以用流導(dǎo)（flow conductance）J來表征，當(dāng)膨脹機(jī)處于臨界狀態(tài)（即壅塞流）時(shí)，膨脹機(jī)的通流能力將不隨膨脹機(jī)出口壓力的降低而增大，流導(dǎo)可以采用下式計(jì)算[19]。

式中，F(xiàn)為摩爾流量，kmol/s，Mi為進(jìn)口氣體分子質(zhì)量，kg/kmol，Ri為進(jìn)口氣體摩爾密度kmol/m3，ΔP為膨脹機(jī)壓差，kPa。

由式(2)可見，ΔP為膨脹機(jī)進(jìn)口壓力Pi與出口壓力Po的差值加上膨脹機(jī)出進(jìn)口靜壓差dPhead與膨脹機(jī)發(fā)生壅塞時(shí)的壓差比率KCPR乘以膨脹機(jī)進(jìn)口壓力Pi之間的最小值。在額定工況下，計(jì)算所得透平仿真模型流導(dǎo)為

1.2.4 輔助設(shè)備

低熱值合成氣燃?xì)廨啓C(jī)其他輔助設(shè)備及管道，如各種泵及管路等均采用DYNSIM 中相關(guān)模型搭建。低熱值合成氣燃?xì)廨啓C(jī)的燃料合成氣來自于IGCC系統(tǒng)的氣化系統(tǒng)，在進(jìn)入燃?xì)廨啓C(jī)燃燒室之前需要在合成氣中注入水蒸汽以調(diào)節(jié)合成氣熱值，降低燃?xì)廨啓C(jī)燃燒過程中的NOx排放。在燃?xì)廨啓C(jī)燃燒室前合成氣管路設(shè)有合成氣混合加熱系統(tǒng)，通過注入水蒸汽調(diào)節(jié)合成氣熱值和溫度。注入蒸汽后，合成氣中水的體積分?jǐn)?shù)控制在30%～33%之間。在實(shí)際運(yùn)行過程以及仿真模型中，進(jìn)入燃?xì)廨啓C(jī)燃燒室的合成氣熱值及壓力均易受到前序工藝即合成氣注蒸汽混合加熱工序的影響。如圖3 所示，即為合成氣混合加熱系統(tǒng)仿真模型，由換熱器模型和壓力調(diào)節(jié)閥模型組成，可以實(shí)現(xiàn)合成氣加熱和壓力調(diào)節(jié)。

圖3 合成氣混合加熱系統(tǒng)仿真模型Fig.3 Simulation model of syngas mixing and heating system

1.3 控制邏輯

1.3.1 主要控制系統(tǒng)

低熱值合成氣燃?xì)廨啓C(jī)仿真模型控制的難點(diǎn)在于對(duì)空氣量、燃料量的精確控制，通過分別編寫六個(gè)主要約束因素的控制邏輯，根據(jù)現(xiàn)場數(shù)據(jù)，不斷優(yōu)化調(diào)整，最終以最小值選擇器的形式實(shí)現(xiàn)燃料量精確控制，如圖4所示。

圖4 低熱值合成氣燃?xì)廨啓C(jī)仿真模型主要控制邏輯Fig.4 Control system of low heating value syngas gas turbine simulation model

1.3.2 燃料切換邏輯

燃?xì)廨啓C(jī)在啟動(dòng)過程、停機(jī)過程、及低負(fù)荷運(yùn)行時(shí)（低于50%額定負(fù)荷），需要運(yùn)行柴油燃料系統(tǒng)，當(dāng)燃?xì)廨啓C(jī)在正常工況下工作時(shí)，運(yùn)行合成氣燃料系統(tǒng)。燃料切換的時(shí)機(jī)和穩(wěn)定混合燃燒的燃料比例選擇不是任意的，應(yīng)遵循下列原則，燃料的選擇必須在啟動(dòng)之前或者在25%額定負(fù)荷以上才能切換。切換的過渡過程為150s。燃料切換功能在DYNSIM 中可用速率限制器（Rate limiter）控制器實(shí)現(xiàn)。切換過程可見圖5。

圖5 燃料切換過程Fig.5 Fuel switching process

2 低熱值合成氣燃?xì)廨啓C(jī)仿真系統(tǒng)驗(yàn)證與應(yīng)用

當(dāng)?shù)蜔嶂岛铣蓺馊細(xì)廨啓C(jī)仿真系統(tǒng)搭建完畢，并根據(jù)現(xiàn)場運(yùn)行數(shù)據(jù)進(jìn)行調(diào)試之后，對(duì)低熱值合成氣燃?xì)廨啓C(jī)仿真系統(tǒng)分別進(jìn)行了不同負(fù)荷下穩(wěn)定運(yùn)行工況，低熱值合成氣燃?xì)廨啓C(jī)仿真模型與現(xiàn)場運(yùn)行數(shù)據(jù)對(duì)比驗(yàn)證，以及燃油切換過程的動(dòng)態(tài)測(cè)試。

2.1 穩(wěn)定運(yùn)行工況靜態(tài)測(cè)試

表1 對(duì)比了合成氣工況100%負(fù)荷（環(huán)境溫度平均約12.6°C，大氣壓力約為101.6kPa）、71.4%負(fù)荷（環(huán)境溫度平均約1.9°C，大氣壓力約為102.0kPa）下，低熱值合成氣燃?xì)廨啓C(jī)仿真模型計(jì)算主要參數(shù)與現(xiàn)場運(yùn)行數(shù)據(jù)（未折合到標(biāo)準(zhǔn)狀況下的實(shí)際值，下同）的對(duì)比。如圖6所示100%負(fù)荷，無論是燃料流量和壓力，燃?xì)廨啓C(jī)轉(zhuǎn)速、負(fù)荷、排煙溫度等與現(xiàn)場運(yùn)行數(shù)據(jù)誤差均小于4%。在71.4%負(fù)荷下，除燃料流量和排煙溫度與現(xiàn)場運(yùn)行數(shù)據(jù)誤差在5%～7%之間，其它參數(shù)與現(xiàn)場運(yùn)行數(shù)據(jù)誤差均不高于4%。表2 對(duì)比了柴油工況50%負(fù)荷下（環(huán)境溫度平均約0.8°C，大氣壓力為102.7kPa），低熱值合成氣燃?xì)廨啓C(jī)仿真模型計(jì)算主要參數(shù)與現(xiàn)場運(yùn)行數(shù)據(jù)的對(duì)比。

表1 合成氣工況，低熱值合成氣燃?xì)廨啓C(jī)仿真模型數(shù)據(jù)與現(xiàn)場運(yùn)行數(shù)據(jù)對(duì)比Tab.1 Comparison of main simulation parameters and on-site experiments parameters under syngas as fuel

圖6 合成氣工況，低熱值合成氣燃?xì)廨啓C(jī)仿真模型數(shù)據(jù)與現(xiàn)場運(yùn)行數(shù)據(jù)誤差Fig.6 The error between simulation and on-site experiments under syngas as fuel

表2 燃油工況，低熱值合成氣燃?xì)廨啓C(jī)仿真模型數(shù)據(jù)與現(xiàn)場運(yùn)行數(shù)據(jù)對(duì)比Tab.2 Comparison of main simulation parameters by and on-site experiments under fuel oil as fuel

由圖7 可見，50%負(fù)荷下燃料流量和壓力，燃?xì)廨啓C(jī)轉(zhuǎn)速、負(fù)荷和與現(xiàn)場運(yùn)行數(shù)據(jù)誤差均小于2%，僅排煙溫度誤差約為9%。在不同燃料工況以及不同負(fù)荷穩(wěn)定運(yùn)行工況下的靜態(tài)參數(shù)對(duì)比，表明搭建的低熱值合成氣燃?xì)廨啓C(jī)仿真系統(tǒng)能夠?qū)崿F(xiàn)對(duì)于現(xiàn)場實(shí)際燃?xì)廨啓C(jī)穩(wěn)定運(yùn)行工況的準(zhǔn)確仿真。

2.2 燃料切換動(dòng)態(tài)測(cè)試

圖7 柴油工況，低熱值合成氣燃?xì)廨啓C(jī)仿真模型數(shù)據(jù)與現(xiàn)場運(yùn)行數(shù)據(jù)誤差Fig.7 The error between simulation and on-site experiments under diesel as fuel

由于合成氣與柴油密度、溫度、壓力等特性的差別，導(dǎo)致燃料切換過程中負(fù)荷波動(dòng)大以及燃燒室內(nèi)壓力波動(dòng)大。通過合成氣混熱加熱仿真模型精確調(diào)節(jié)燃料切換過程中合成氣特性參數(shù)的穩(wěn)定，以及嚴(yán)格控制燃料量基準(zhǔn)及燃料閥門的前后壓差，使得燃燒室內(nèi)燃燒熱值穩(wěn)定。實(shí)現(xiàn)燃?xì)廨啓C(jī)燃料切換過程平穩(wěn)。圖8顯示了低熱值合成氣燃?xì)廨啓C(jī)燃料切換過程中燃料流量、合成氣熱值和壓力的變化情況?？梢姡?50s 內(nèi)，柴油流量從最大約8.5kg/s 降至0，而合成氣流量從0 升至最高50.3kg/s，同時(shí)合成氣熱值、壓力波動(dòng)均小于1%。由圖9可見，燃?xì)廨啓C(jī)燃料切換過程前后，負(fù)荷均穩(wěn)定在115MW，燃料切換時(shí)，低熱值合成氣燃?xì)廨啓C(jī)仿真模型負(fù)荷波動(dòng)控制在5MW 之內(nèi)，仿真模型負(fù)荷與現(xiàn)場實(shí)際運(yùn)行負(fù)荷（環(huán)境溫度平均約-0.7°C，大氣壓力為102.8kPa）平均誤差小于1%，最大誤差小于10%，基本符合現(xiàn)場燃料切換過程。

圖8 燃料切換過程中燃料流量、合成氣熱值和壓力的變化Fig.8 Flow rates of fuel,syngas heating value and pressure during fuel switchover process

圖9 燃料切換過程燃?xì)廨啓C(jī)仿真模型和現(xiàn)場實(shí)際負(fù)荷對(duì)比Fig.9 Comparison of load between gas turbine simulation model and on-site experiments during fuel switchover process

2.3 應(yīng)用

搭建的低熱值合成氣燃?xì)廨啓C(jī)仿真系統(tǒng)已多次應(yīng)用于天津IGCC 操作人員培訓(xùn)，幫助現(xiàn)場運(yùn)行人員熟悉掌握燃?xì)廨啓C(jī)開停車、升降負(fù)荷以及燃料切換操作，減少人為失誤，提升運(yùn)行效率和穩(wěn)定性。利用該仿真系統(tǒng)還可以對(duì)燃?xì)廨啓C(jī)進(jìn)行的技術(shù)改造預(yù)先進(jìn)行仿真驗(yàn)證，以減少現(xiàn)場試驗(yàn)的風(fēng)險(xiǎn)和成本。

3 結(jié)論

以國內(nèi)首臺(tái)天津IGCC關(guān)鍵設(shè)備低熱值合成氣燃?xì)廨啓C(jī)為研究對(duì)象，通過機(jī)理分析等手段，基于商業(yè)動(dòng)態(tài)過程模擬系統(tǒng)搭建了針對(duì)低熱值合成器燃?xì)廨啓C(jī)的仿真系統(tǒng)。合成氣工況100%負(fù)荷、71.4%負(fù)荷，柴油50%負(fù)荷穩(wěn)定運(yùn)行工況下，低熱值合成氣燃?xì)廨啓C(jī)仿真系統(tǒng)運(yùn)行數(shù)據(jù)與現(xiàn)場實(shí)際運(yùn)行數(shù)據(jù)相比，絕大部分主要運(yùn)行參數(shù)誤差小于5%。在115MW 進(jìn)行的燃料切換過程動(dòng)態(tài)測(cè)試，低熱值合成氣燃?xì)廨啓C(jī)仿真系統(tǒng)負(fù)荷波動(dòng)控制在5MW 之內(nèi)，仿真系統(tǒng)負(fù)荷與現(xiàn)場實(shí)際運(yùn)行負(fù)荷平均誤差小于1%，最大誤差小于10%，基本符合現(xiàn)場燃料切換過程。這表明，搭建的低熱值合成氣燃?xì)廨啓C(jī)仿真系統(tǒng)能夠?qū)崿F(xiàn)對(duì)于現(xiàn)場實(shí)際燃?xì)廨啓C(jī)穩(wěn)定運(yùn)行工況以及動(dòng)態(tài)過程的準(zhǔn)確仿真。該仿真系統(tǒng)已應(yīng)用于現(xiàn)場操作人員的培訓(xùn)，以提高運(yùn)行效率和穩(wěn)定性。