劉云生，馮永明，陳華清，劉祥源，宋漢江

(1.92537部隊(duì)，北京100161;2.哈爾濱工程大學(xué)能源與動(dòng)力工程學(xué)院，黑龍江哈爾濱15001)

0 引言

增壓鍋爐裝置由增壓鍋爐和為其提供增壓助燃空氣的煙氣渦輪增壓機(jī)組組成。與常規(guī)鍋爐裝置相比，增壓鍋爐裝置在單機(jī)功率、體積、重量、經(jīng)濟(jì)性、機(jī)動(dòng)性等方面占有較大優(yōu)勢，適用于大型水面艦船和偏遠(yuǎn)地區(qū)的大型工廠電站，目前已經(jīng)在俄羅斯海軍和大型工廠的電站系統(tǒng)中得到廣泛應(yīng)用。作為增壓鍋爐裝置的重要組成部分，增壓鍋爐和渦輪增壓機(jī)組的熱力耦合關(guān)系極其復(fù)雜:一方面，增壓鍋爐向渦輪增壓機(jī)組煙氣渦輪和輔助汽輪機(jī)提供的煙氣和蒸汽的熱力性能參數(shù)取決于鍋爐負(fù)荷狀態(tài)和性能;另一方面，鍋爐性能又受到煙氣渦輪和輔助汽輪機(jī)共同驅(qū)動(dòng)的壓氣機(jī)的影響。煙氣渦輪增壓機(jī)組和增壓鍋爐相互聯(lián)系又相互影響，是一個(gè)有機(jī)整體。國內(nèi)在增壓鍋爐裝置研究方面研究較少，積累不多［1－6］，更缺乏從系統(tǒng)層面對(duì)增壓鍋爐和渦輪增壓機(jī)組匹配特性的深入分析。本文以船用增壓鍋爐裝置為研究對(duì)象，建立增壓鍋爐、渦輪增壓機(jī)組及相應(yīng)控制系統(tǒng)、輔助系統(tǒng)的仿真模型，對(duì)增壓鍋爐和渦輪增壓機(jī)組穩(wěn)動(dòng)態(tài)匹配問題進(jìn)行深入系統(tǒng)的仿真分析。

1 仿真模型及校驗(yàn)

1.1 仿真模型及關(guān)鍵技術(shù)

渦輪增壓機(jī)組主要由壓氣機(jī)、煙氣渦輪和輔助汽輪機(jī)組成，增壓鍋爐主要由爐膛、過熱器、經(jīng)濟(jì)器等組成。上述設(shè)備或部件的數(shù)學(xué)模型都已比較成熟，限于篇幅這里不再一一列出。下面重點(diǎn)對(duì)部分仿真模型中需要解決的關(guān)鍵技術(shù)進(jìn)行分析。本文在通用熱能動(dòng)力系統(tǒng)仿真平臺(tái)SimuWorks上開發(fā)的渦輪增壓機(jī)組與增壓鍋爐協(xié)配性仿真程序組態(tài)圖如圖1所示(給水、燃油等輔助系統(tǒng)的組態(tài)圖略去)，其中鍋爐各部件、輔助汽輪機(jī)采用軟件模型庫中的模型搭建，部分鍋爐部件模型根據(jù)增壓鍋爐的特點(diǎn)進(jìn)行了修正;壓氣機(jī)、煙氣渦輪、容積模塊等仿真模塊是在仿真平臺(tái)的基礎(chǔ)上利用Fortan語言2次開發(fā)得到。

圖1 渦輪增壓機(jī)組與增壓鍋爐協(xié)配性仿真程序組態(tài)圖Fig.1 Simulation program configuration figure of turbin charger and superboiler

在蒸汽動(dòng)力裝置建模中，需要解決的難點(diǎn)之一是解算由蒸汽、水、油、煙氣等流體工質(zhì)所組成的流體網(wǎng)絡(luò)［7］。由于SimuWorks軟件模型庫中自帶的仿真模型已經(jīng)包含流體網(wǎng)絡(luò)計(jì)算功能，這就需要解決自編模塊的流體網(wǎng)絡(luò)計(jì)算問題。由于反映壓氣機(jī)容積慣性影響的容積模塊建模相對(duì)簡單，所以本文重點(diǎn)解決了煙氣渦輪模塊和壓氣機(jī)模塊的流體網(wǎng)絡(luò)計(jì)算問題。

蒸汽、空氣、煙氣等都是可壓縮流體，可壓縮流體支路(所謂支路就是流體流經(jīng)的一條通道)的線性化模型可由下式［7］表示:

因?yàn)?個(gè)蒸汽動(dòng)力系統(tǒng)中有很多條流體支路，而不同的流體支路之間有交匯點(diǎn)(又稱節(jié)點(diǎn))，通過交匯點(diǎn)的流量要符合質(zhì)量守恒方程，所以，流體網(wǎng)絡(luò)求解的本質(zhì)就是求解系統(tǒng)中所有節(jié)點(diǎn)按照質(zhì)量守恒方程建立起來的多個(gè)方程組，而由式(1)可見，質(zhì)量流量又取決于Rb1和Rb2。在本仿真平臺(tái)中，流體網(wǎng)絡(luò)方程組的解算由仿真平臺(tái)完成，支路特性參數(shù)需要各模塊計(jì)算。所以這里需要重點(diǎn)解決的就是煙氣渦輪和壓氣機(jī)的支路特性計(jì)算問題。

1)煙氣渦輪模塊

通過分析簡化，渦輪模塊的流網(wǎng)特性參數(shù)的計(jì)算式如下:

其中:GTp為渦輪流量特性曲線中的折合流量;G0為流量參考值。

根據(jù)式(2)，通過煙氣渦輪進(jìn)口的壓力溫度值可以計(jì)算得到渦輪的膨脹比，根據(jù)渦輪特性圖可以通過插值獲得相應(yīng)的折合流量，然后就可求得相應(yīng)的Rb1值。

2)壓氣機(jī)模塊

壓氣機(jī)流體網(wǎng)絡(luò)處理的方法相對(duì)簡單，由于壓氣機(jī)的進(jìn)口為外界大氣，壓力已知，因此壓氣機(jī)不參與流網(wǎng)計(jì)算，只作為一個(gè)邊界模塊(即其壓力已知)，計(jì)算壓氣機(jī)后的空氣壓力，提供給流體網(wǎng)絡(luò)，由流網(wǎng)計(jì)算出壓氣機(jī)所需要的空氣流量，然后再把這個(gè)流量直接賦給壓氣機(jī)進(jìn)口。

仿真模塊運(yùn)行時(shí)，各個(gè)部分間的運(yùn)行關(guān)系如圖2所示，由控制系統(tǒng)控制輔助汽輪機(jī)閥門前的開度，以此來調(diào)節(jié)輔助汽輪機(jī)發(fā)出的功率，輔助汽輪機(jī)發(fā)出的功率傳遞給轉(zhuǎn)動(dòng)慣量模塊，在轉(zhuǎn)動(dòng)慣量模塊中，根據(jù)輔助汽輪機(jī)和煙氣渦輪機(jī)發(fā)出的功率以及壓氣機(jī)消耗的功率等關(guān)系，得到壓氣機(jī)的轉(zhuǎn)速和輔助汽輪機(jī)的轉(zhuǎn)速;壓氣機(jī)模塊根據(jù)得到的轉(zhuǎn)速和流網(wǎng)算出來的流量，通過插值得到壓氣機(jī)的壓比和效率，進(jìn)一步計(jì)算出壓氣機(jī)消耗的功率，排氣溫度，最后將得到的計(jì)算值往下一級(jí)傳遞(物質(zhì)端口接容積慣性模塊，信息端口接轉(zhuǎn)動(dòng)慣量模塊);容積模塊根據(jù)壓氣機(jī)的容積慣性，計(jì)算壓氣機(jī)實(shí)際出口的空氣值，然后經(jīng)壓氣機(jī)排氣管路進(jìn)入鍋爐爐膛;完成鍋爐相關(guān)熱力計(jì)算后，鍋爐爐膛經(jīng)濟(jì)器出口接煙氣渦輪機(jī)，煙氣渦輪仿真模塊根據(jù)爐膛經(jīng)濟(jì)器出口的煙氣焓值、壓力、流量等，計(jì)算渦輪效率和渦輪發(fā)出的功率等，最后煙氣排入大氣，渦輪產(chǎn)生的功率信號(hào)傳遞給轉(zhuǎn)動(dòng)慣量模塊。到此，整個(gè)渦輪增壓機(jī)組和增壓鍋爐包括輔助汽輪機(jī)在內(nèi)，完成整個(gè)的循環(huán)計(jì)算。

圖2 仿真計(jì)算流程圖Fig.2 Simulation calculation process

1.2 仿真模型校驗(yàn)

為了準(zhǔn)確評(píng)價(jià)仿真模型的準(zhǔn)確性，本文將仿真值和測試數(shù)據(jù)進(jìn)行了對(duì)比分析，根據(jù)對(duì)比分析結(jié)果可知，額定工況下各主要參數(shù)的仿真值與測試值相差不大;非額定工況下，燃油耗量、空氣耗量仿真值比設(shè)計(jì)值略大一些，煙氣渦輪功率增大，會(huì)導(dǎo)致輔助汽輪機(jī)需要功率減少，計(jì)算得到的輔助汽輪機(jī)蒸汽流量偏小?？傮w而言，所建立的仿真模型基本可信。

2 渦輪增壓機(jī)組與增壓鍋爐協(xié)配性分析

本文利用上述仿真模型，對(duì)渦輪增壓機(jī)組和增壓鍋爐進(jìn)行穩(wěn)、動(dòng)態(tài)仿真分析。圖3是30℃環(huán)境溫度(夏季工況)和10℃環(huán)境溫度(冬季工況)2種穩(wěn)態(tài)工況情況下，隨著鍋爐負(fù)荷變化渦輪增壓機(jī)組主要性能參數(shù)的變化規(guī)律。從圖中可以看出，無論是夏季工況還是冬季工況，輔助汽輪機(jī)始終需投入運(yùn)行，即渦輪增壓機(jī)組不存在功率平衡點(diǎn);其次，夏季工況，隨著鍋爐負(fù)荷增加即鍋爐空氣流量增大，壓氣機(jī)消耗的功率逐漸增大，盡管煙氣渦輪發(fā)出的功率也逐漸增大，但二者的差值逐漸增大，從而導(dǎo)致輔助汽輪機(jī)需補(bǔ)充的功率逐漸加大，在機(jī)組全工況運(yùn)行時(shí)，輔助汽輪機(jī)需補(bǔ)充的功率達(dá)到最大值;而在冬季工況，隨著鍋爐空氣流量增大，壓氣機(jī)消耗的功率逐漸增大，但由于煙氣渦輪發(fā)出的功率逐漸增大的更快，導(dǎo)致二者的差值先增大后減小，即較高負(fù)荷時(shí)，輔助汽輪機(jī)需補(bǔ)充的功率達(dá)到最大值。實(shí)際上根據(jù)壓氣機(jī)折合流量計(jì)算公式

圖4給出了機(jī)組升工況加速過程渦輪增壓機(jī)組主要性能參數(shù)和工作參數(shù)隨鍋爐負(fù)荷變化的瞬態(tài)變化規(guī)律。由圖4～圖7可以看出，在機(jī)組加速過程中，為滿足鍋爐空氣流量的需要，壓氣機(jī)的空氣流量增大，由于系統(tǒng)超調(diào)，其間參數(shù)波動(dòng)劇烈:在鍋爐參數(shù)切換過程中，壓氣機(jī)空氣流量首先急劇增大，隨后再減小并又增大，產(chǎn)生這種復(fù)雜波動(dòng)的原因?yàn)榧铀僦噶钍馆o助汽輪機(jī)功率增大，壓氣機(jī)流量增大鍋爐負(fù)荷升高，渦輪做功能力增強(qiáng)，自動(dòng)調(diào)節(jié)系統(tǒng)關(guān)小輔助汽輪機(jī)進(jìn)氣閥門以避免機(jī)組超速，鍋爐達(dá)到系統(tǒng)滑參數(shù)運(yùn)行的汽包壓力目標(biāo)，主汽輪機(jī)調(diào)節(jié)閥門開大其消耗的蒸汽流量相對(duì)瞬時(shí)增加，輔助汽輪機(jī)的蒸汽流量減小，其做功能力下降，機(jī)組轉(zhuǎn)速下降壓氣機(jī)空氣流量減小，煙氣渦輪的做功能力下降，當(dāng)機(jī)組提供的空氣不能滿足鍋爐運(yùn)行的需要時(shí)，自動(dòng)調(diào)節(jié)系統(tǒng)控制輔助汽輪機(jī)進(jìn)氣閥門開大，機(jī)組轉(zhuǎn)速升高壓氣機(jī)空氣流量再次增大，最終機(jī)組穩(wěn)定運(yùn)行，且提供滿足鍋爐需要的空氣。此外，煙氣渦輪的進(jìn)口壓力、溫度、流量的變化趨勢與壓氣機(jī)的也基本一致，而且在整個(gè)鍋爐負(fù)荷渦輪增壓機(jī)組加速過程，輔助汽輪機(jī)需補(bǔ)充的功率遠(yuǎn)大于穩(wěn)態(tài)工況的功率平衡計(jì)算值，而這主要是由于加速過程中煙氣渦輪發(fā)出功率滯后，機(jī)組轉(zhuǎn)子加速過程也要消耗功率造成的。這表明僅僅依靠穩(wěn)態(tài)功率平衡計(jì)算是不能準(zhǔn)確確定渦輪增壓機(jī)組功率平衡即輔助汽輪機(jī)需補(bǔ)充的功率及其蒸汽消耗量的，足夠精度的動(dòng)態(tài)匹配計(jì)算是渦輪增壓機(jī)組與增壓鍋爐匹配設(shè)計(jì)所必需的。

3 結(jié)語

本文重點(diǎn)解決了壓氣機(jī)、渦輪模塊的流體網(wǎng)絡(luò)計(jì)算與仿真平臺(tái)SimuWorks流體網(wǎng)絡(luò)計(jì)算方法的匹配問題，在仿真平臺(tái)上完成了渦輪增壓機(jī)組與增壓鍋爐二者聯(lián)合工作仿真系統(tǒng)的開發(fā)和集成，通過實(shí)驗(yàn)數(shù)據(jù)驗(yàn)證了仿真模型的正確性，并在此基礎(chǔ)上研究分析了環(huán)境溫度和鍋爐負(fù)荷變化等對(duì)增壓鍋爐與渦輪增壓機(jī)組的穩(wěn)動(dòng)態(tài)匹配特性的影響，得出的結(jié)論可用于指導(dǎo)增壓鍋爐渦輪增壓機(jī)組的研制和試驗(yàn)。

［1］李章，等.艦用增壓鍋爐裝置［M］.北京:海潮出版社，2009.LI Zhang，et al.Naval supercharged boiler installation.Beijing:Haichao Press，2009.

［2］張利劍，等.增壓鍋爐經(jīng)濟(jì)器結(jié)構(gòu)與熱力計(jì)算研究［J］.哈爾濱工程大學(xué)學(xué)報(bào)，2004，25(2):229－231.ZHANG Li-jian，et al.Research on the structure and thermal calculation method of economizer in supercharged boiler［J］.Journal of Harbin Engineering University，2004，25(2):229－231.

［3］姜任秋，等.艦用增壓鍋爐煙風(fēng)阻力計(jì)算方法研究［J］.哈爾濱工程大學(xué)學(xué)報(bào)，2004，25(5):566－568.JIANG Ren-qiu，et al.Calculation methods of gas resistance of supercharged boiler［J］.Journal of Harbin Engineering University，2004，25(5):566 －568.

［4］姜任秋，等.船用主增壓鍋爐煙氣工作特性分析與研究［J］.哈爾濱工程大學(xué)學(xué)報(bào)，2005，26(6):722 －725.JIANG Ren-qiu，et al. Analysis of flue gas working characteristic in marine primary supercharged boiler［J］.Journal of Harbin Engineering University，2005，26(6):722－725.

［5］李彥軍，等.船用增壓鍋爐熱力計(jì)算方法有關(guān)問題分析［J］.熱能動(dòng)力工程，2005，26(1):93 －96.LI Yan-jun，et al.A study of some issues related to the thermodynamic calculation of a super charged marine boiler［J］.Journal of Engineering for Thermal Energy＆ Power，2005，26(1):93 －96.

［6］韓靜，等.船用增壓鍋爐的熱平衡計(jì)算方法研究［J］.工程熱物理學(xué)報(bào)，2004，25(s):173 －175.HAN Jing，et al.Research of heat calculation method of the marine supercharged boiler［J］.Journal of Engineering Thermophysics，2004，25(s):173 －175.

［7］呂崇德，任挺進(jìn)，姜學(xué)智，等.大型火電機(jī)組系統(tǒng)仿真與建模［M］.清華大學(xué)出版社，2002.LV Chong-de，REN Ting-jin，JIANG Xue-zhi，etal.Simulation and modeling for power plant［M］.Beijing:Tsinghua University Press，2002.