楊 東，李淑英，苗安利，王志濤

（1.海軍駐西安地區(qū)航空軍事代表室，西安 710021；2.哈爾濱工程大學(xué)動力與能源工程學(xué)院，哈爾濱 150001）

0 引言

目前，燃?xì)廨啓C(jī)在世界船舶上大量使用，并因可以提高經(jīng)濟(jì)性而受到越來越多的重視。在能源危機(jī)之后，余熱鍋爐-透平發(fā)電和推進(jìn)系統(tǒng)的余熱回收成為船舶動力裝置最有效的1種節(jié)能措施[1]。因此，對船舶燃?xì)?蒸汽聯(lián)合循環(huán)(COGAS)進(jìn)行研究和性能分析，可以為其在中國船舶中的實際應(yīng)用提供一定的參考[2]。

1 艦船COGAS裝置工作原理

COGAS裝置循環(huán)系統(tǒng)原理如圖1所示。利用燃?xì)廨啓C(jī)排出的燃?xì)饧訜釓U氣鍋爐中的水，在完成余熱回收后排入大氣。在廢氣鍋爐中，水由除氧器經(jīng)給水泵進(jìn)入經(jīng)濟(jì)器，在其中加熱到接近或達(dá)到飽和溫度，進(jìn)入蒸發(fā)器中蒸發(fā)，并在過熱器中過熱到額定參數(shù)，然后進(jìn)入蒸汽輪機(jī)作功；作功后的乏汽排入冷凝器凝結(jié)，再由凝水泵將其排入除氧器加熱、除氧后，經(jīng)給水泵送入廢氣鍋爐。在蒸汽輪機(jī)中有時還設(shè)置1個撤汽點，將部分蒸汽抽出用于船上生活[2]。

2 艦船COGAS的仿真建模

國外開發(fā)出ASPENPLUS一類的專門用于聯(lián)合循環(huán)仿真的平臺，以及GTPro和GTMaster等專門用來進(jìn)行燃?xì)廨啓C(jī)聯(lián)合循環(huán)發(fā)電機(jī)組的方案設(shè)計和性能計算（包括變工況計算）的軟件。國內(nèi)對艦船COGAS研究得很少[3]，原因是沒有通用仿真軟件平臺。為此，本文以Matlab/Simulink作為仿真平臺，采用微分方程的方法，對COGAS進(jìn)行建模和仿真，并深入了解其穩(wěn)態(tài)與動態(tài)特性。模塊化仿真是近年來發(fā)展起來的1種計算機(jī)仿真建模方法。其核心思想是按照一定規(guī)則對大系統(tǒng)進(jìn)行模塊劃分，各模塊相互獨立，實現(xiàn)特定的功能，模塊具有通用性，可重復(fù)使用。COGAS系統(tǒng)模塊劃分如圖2所示。

根據(jù)如圖2所示的COGAS系統(tǒng)模塊的劃分，分別建立各部件的數(shù)學(xué)模型，并以Matlab/Simulink作為仿真平臺，根據(jù)其數(shù)學(xué)模型建立該系統(tǒng)的整體仿真模型，如圖3所示。

3 COGAS的特性分析

影響COGAS機(jī)組變工況性能的因素很多，主要包括機(jī)組負(fù)荷、環(huán)境條件(包括大氣溫度、大氣壓力和大氣濕度)、機(jī)組燃料、余熱鍋爐蒸汽循環(huán)方式、燃?xì)廨啓C(jī)進(jìn)出口壓損等[4]。本文以燃?xì)廨啓C(jī)的負(fù)荷變化為變量研究COGAS的穩(wěn)態(tài)與動態(tài)特性。

3.1 COGAS的穩(wěn)態(tài)性能分析

COGAS總體設(shè)計的首要工作就是進(jìn)行COGAS總體性能分析，得到燃?xì)廨啓C(jī)在設(shè)計工況及變工況下的余熱鍋爐所能產(chǎn)生的蒸汽流量、溫度、壓力以及汽輪機(jī)發(fā)出的功率等。因此，利用所編制COGAS整體仿真軟件，進(jìn)行COGAS穩(wěn)態(tài)特性仿真分析，為COGAS總體設(shè)計提供設(shè)計基礎(chǔ)。

在燃?xì)廨啓C(jī)設(shè)計點即100%負(fù)荷工況下進(jìn)行COGAS性能分析，結(jié)果見表1。表中機(jī)組主要參數(shù)的仿真計算值與額定設(shè)計參數(shù)的相對誤差都不超過2%；若要進(jìn)行其他負(fù)荷工況下的穩(wěn)態(tài)計算，只要將負(fù)載輸入量改成相應(yīng)的負(fù)荷即可。因此，用部件通用仿真模塊搭建的COGAS系統(tǒng)整體模型能夠較好地滿足穩(wěn)態(tài)仿真的要求。

表1 機(jī)組100%負(fù)荷工況下穩(wěn)態(tài)仿真結(jié)果與設(shè)計工況參數(shù)值比較

COGAS參數(shù)在不同工況下的變化如圖4所示。從圖中可見，隨著燃?xì)廨啓C(jī)負(fù)荷的減小，蒸汽流量、蒸汽溫度、聯(lián)合效率、汽輪機(jī)功率和機(jī)組總功率都減小。這是由于燃?xì)廨啓C(jī)負(fù)荷減小，使燃?xì)饬髁亢腿細(xì)鉁囟榷冀档?，?dǎo)致燃?xì)鈧?cè)放熱量減少，可回收的余熱量減少，使蒸汽側(cè)吸熱量也減少，從而使蒸汽產(chǎn)量降低。雖然燃?xì)鉁囟冉档筒粫拐羝隹跍囟冉档停侨細(xì)饬髁拷档蜁?dǎo)致燃?xì)鈱α鞣艧嵯禂?shù)減小，燃?xì)庠谶^熱區(qū)放熱量減少。因此，在過熱區(qū)，出口蒸汽溫度降低。由于蒸汽出口的流量和溫度都降低，則汽輪機(jī)的輸出功率也降低。而機(jī)組功率是汽輪機(jī)功率和燃?xì)廨啓C(jī)功率的總和，所以機(jī)組功率也降低。從圖4（c）、（d）中可見，采用COGAS后，在額定工況下，機(jī)組功率增加31%，效率增加到47%；在50%工況下，機(jī)組功率增加29%，效率增加到39%。由此可知，無論是在高工況還是低工況下，燃?xì)廨啓C(jī)采用COGAS后，其效率都會顯著提高。

3.2 COGAS的動態(tài)仿真

為了更深入研究COGAS系統(tǒng)的運(yùn)行特性，分析燃?xì)廨啓C(jī)負(fù)荷變化時，COGAS各參數(shù)的動態(tài)特性。燃?xì)廨啓C(jī)由100%負(fù)荷突降到80%負(fù)荷時,各主要參數(shù)變化規(guī)律如圖5所示。

從圖5中可見，由于燃?xì)廨啓C(jī)的負(fù)荷突然減小，使燃?xì)廨啓C(jī)的輸出功率和噴油量減小，并且經(jīng)過一段時間后達(dá)到新的平衡。這是由于燃?xì)廨啓C(jī)負(fù)荷變化引起其轉(zhuǎn)速變化，控制模塊通過改變噴油量使轉(zhuǎn)速回到設(shè)定值，噴油量的改變使燃?xì)廨啓C(jī)的輸出功率與其負(fù)荷相匹配，最后在該工況下穩(wěn)定運(yùn)行。

隨著燃?xì)廨啓C(jī)負(fù)荷的突降，使蒸汽流量、壓力和溫度都降低，并在一段時間后達(dá)到新的穩(wěn)定。這是由于燃?xì)廨啓C(jī)負(fù)荷的減小，使燃?xì)饬髁亢腿細(xì)鉁囟榷冀档?，?dǎo)致燃?xì)鈧?cè)放熱量減少，可回收的余熱量減少，所以蒸汽側(cè)吸熱量也減少，從而使蒸汽產(chǎn)量和壓力降低。雖然燃?xì)鉁囟鹊慕档筒粫拐羝隹跍囟冉档?，但是燃?xì)饬髁拷档蜁?dǎo)致燃?xì)鈱α鞣艧嵯禂?shù)減小，燃?xì)庠谶^熱區(qū)放熱量減少，因此在過熱區(qū)，出口蒸汽溫度降低。由于蒸汽出口的流量和溫度都降低，則汽輪機(jī)的輸出功率減小。而當(dāng)負(fù)荷大于50%時，蒸汽系統(tǒng)采用滑壓運(yùn)行，所以主汽門開度為1，并保持不變。當(dāng)燃?xì)廨啓C(jī)達(dá)到新的穩(wěn)定運(yùn)行以后，蒸汽系統(tǒng)的各參數(shù)也達(dá)到新的穩(wěn)定工作狀態(tài)。

燃?xì)廨啓C(jī)負(fù)荷由80%突增到100%時，各參數(shù)的變化規(guī)律如圖6所示。這些參數(shù)都是先達(dá)到1個極限值，而后達(dá)到新的平衡。這是因為燃?xì)廨啓C(jī)負(fù)荷的增大，使燃?xì)饬髁吭龃蠛蜏囟壬?，COGAS的幾個參數(shù)都增大，并且隨著燃?xì)廨啓C(jī)負(fù)荷的增大達(dá)到新的穩(wěn)定工作狀態(tài)。

4 結(jié)束語

采用模塊化建模的方式對COGAS系統(tǒng)進(jìn)行了穩(wěn)態(tài)與動態(tài)特性分析，得到了在燃?xì)廨啓C(jī)變工況下的COGAS穩(wěn)態(tài)與動態(tài)特性。由分析結(jié)果可知，無論是在高工況下，還是在其低工況下，采用COGAS后，燃?xì)廨啓C(jī)效率都顯著提高。因此，COGAS在中國未來船舶推進(jìn)動力裝置的應(yīng)用中具有廣闊的前景。

[1]徐利軍，曹渝白.船舶余熱動力回收系統(tǒng)熱力學(xué)參數(shù)及余熱鍋爐結(jié)構(gòu)參數(shù)優(yōu)化[J].熱能動力工程，1992，7（2）：79-82.

[2]蔡振飛，朱秀章.艦船燃?xì)?蒸汽聯(lián)合循環(huán)(COGAS)裝置的探討[J].船舶工程，1978（1）:34-41.

[3]張俊禮.S109FA型燃?xì)猓羝?lián)合循環(huán)機(jī)組的全工況仿真[D].南京：東南大學(xué)，2005.

[4]常靜華.燃?xì)?蒸汽聯(lián)合循環(huán)機(jī)組變工況建模及其特性分析[D].云南：昆明理工大學(xué)，2008.

[5]Surgeor B W.A Compact model of power house boiler [J].Transactions of the Society for Computer Simulation,1987,4（4）：271- 297.

[6]Akber Pasha, Sanjeev Jolly.Combined cycle heat recovery steam generators optimum capabilities and selection criteria [J]. Heat Recovery Systems＆CHP, 1995,15（2）：147- 154.