汽輪機變工況模型的簡便建立方法及應用

余興剛，李 旭，蔣北華，李 明，謝誕梅

?

汽輪機變工況模型的簡便建立方法及應用

余興剛1，李 旭2，蔣北華2，李 明2，謝誕梅3

（1.高效清潔火力發(fā)電技術湖南省重點實驗室，湖南 長沙 410007；2.湖南省湘電試驗研究院有限公司，湖南 長沙 410004；3.武漢大學動力與機械學院，湖北 武漢 430072）

鑒于汽輪機詳細變工況計算方法在實際計算過程中很難開展，本文提出了一種根據(jù)機組熱平衡圖構造變工況計算模型的簡便方法，并采用該方法建立了汽輪機定流量和定功率變工況計算模型。以某超臨界660 MW機組為例，運用所建立的定功率變工況模型計算得到機組典型工況下的各參數(shù)值，計算結果與設計值相比誤差較?。贿\用所建立的定流量變工況模型計算得到了四閥全開工況下主蒸汽壓力、主蒸汽溫度、再熱蒸汽溫度、再熱壓損和排汽壓力變化對熱耗率和功率的修正曲線，與制造廠提供的修正曲線相比誤差較小。這表明所建立的變工況計算模型精度較高，滿足工程計算的要求。

汽輪機；變工況；熱平衡圖；定流量；定功率；修正曲線

由于環(huán)境溫度、負荷、設備狀態(tài)等的變化，汽輪機通常處于變工況運行狀態(tài)，導致實際運行初終參數(shù)偏離設計值，這將對汽輪機運行的經(jīng)濟性產(chǎn)生一定的影響[1-3]?，F(xiàn)有研究汽輪機參數(shù)偏離對機組經(jīng)濟性影響的方法主要有特性試驗法、特性曲線法、偏導算法、變工況算法等[4]。其中，變工況算法應用較為廣泛，然而汽輪機詳細變工況計算方法考慮的因素較多，精度雖較高，但需要汽輪機及其附屬設備的詳細結構尺寸數(shù)據(jù)[5]，使汽輪機詳細變工況計算很難開展[6-7]。鑒于此，本文提出了一種僅根據(jù)機組熱平衡圖構造汽輪機變工況計算模型的簡便方法，采用該方法建立了汽輪機定流量和定功率變工況計算模型，并將此模型應用至某超臨界660 MW機組，計算得到了典型工況下的各參數(shù)值和四閥全開工況下機組初終參數(shù)對熱耗率和功率的修正曲線。

1 汽輪機變工況計算模型

根據(jù)各級組的變工況特性，將汽輪機各級組分為調節(jié)級、中間級組和末級組3類考慮[8]。

1.1 調節(jié)級變工況計算

參照文獻[9]，根據(jù)熱平衡圖上典型工況下的調節(jié)級后壓力和焓值，繪制出調節(jié)級焓降與調節(jié)級壓比的關系曲線，并將該曲線擬合成相應的公式作為調節(jié)級變工況計算的依據(jù)。當主蒸汽參數(shù)偏離額定值時，調節(jié)級的實際壓比應采用式(1)進行修正[10]：

式中：為蒸汽壓力，MPa；下標0表示主蒸汽參數(shù)，g表示調節(jié)級參數(shù)，1表示實際工況參數(shù)；為蒸汽比熱容，m3/kg；為過熱蒸汽等熵指數(shù)。

考慮到調節(jié)級后為高溫高壓蒸汽，蒸汽性質偏離理想氣體特性較遠，采用改進后的弗留格爾公式（式(2)）計算工況變化后的調節(jié)級后蒸汽壓力[11]：

式中：表示調節(jié)級后蒸汽流量，t/h；1為一段抽汽壓力，MPa

1.2 中間級組變工況計算

現(xiàn)有文獻計算變工況下中間各級組抽汽口處蒸汽壓力時多采用弗留格爾公式，然而由于低壓缸各級組壓比通常比高中壓缸各級組小，當工況偏離較大時，會出現(xiàn)弗留格爾公式根號內為負值的情況，導致計算無法進行。鑒于此，本文參照ASME PTC6—2004[12]推薦的方法計算變工況下各抽汽口處和各汽缸入口處蒸汽壓力，根據(jù)熱耗率驗收工況（THA工況）的熱平衡圖，采用式(4)計算各抽汽口處和汽缸入口處流量系數(shù)：

式中為抽氣口處蒸汽濕度。

當汽輪機工況發(fā)生變化后，根據(jù)THA工況下各抽汽口處和汽缸入口處流量系數(shù)便可計算得到對應位置的蒸汽壓力。筆者計算時發(fā)現(xiàn)，如果所有抽汽口和汽缸入口處蒸汽壓力均用式(4)計算，迭代計算將無法收斂。經(jīng)嘗試發(fā)現(xiàn)，當一段抽汽壓力改用弗留格爾公式或改進后的弗留格爾公式[11]計算時，計算便可收斂。鑒于此，當計算一段抽汽壓力時采用改進后的弗留格爾公式，而其他抽汽口處和汽缸入口處蒸汽壓力則采用式(4)計算。

對于汽輪機中間各級組，當工況發(fā)生變化后級組的前后壓比變化較小，級組的效率變化也較小。因此，在本文的變工況計算模型中，汽輪機中間各級組的效率取熱平衡圖上各工況下的平均值[13]。

1.3 末級組變工況計算

汽輪機末級組變工況計算的目的主要為確定機組的排汽焓。由于缺乏末級組的詳細結構尺寸數(shù)據(jù)，考慮到影響末級組效率的因素主要為末級組壓比，本文參照文獻[14]，根據(jù)熱平衡圖上典型工況下的數(shù)據(jù)，繪制出末級組效率隨壓比變化的曲線，并將該曲線擬合成相應的公式供編程計算時使用。

1.4 其他系統(tǒng)變工況計算

本文在計算回熱系統(tǒng)各汽水參數(shù)時，未考慮加熱器的變工況，并假定各加熱器的上、下端差和各抽汽管道壓損均為設計值。當計算給水泵汽輪機進汽量時，假定給水泵汽輪機和給水泵效率不變，工況變化后的給水泵出口壓力根據(jù)式(5)計算[15]：

式中：gs為給水壓力，MPa；gs為給水流量，t/h；zq為主蒸汽壓力，MPa。

在統(tǒng)計熱平衡圖中各工況下的汽封漏汽量時發(fā)現(xiàn)，過橋汽封漏汽總量隨主蒸汽流量的變化、過橋汽封至中壓缸漏汽量隨中壓缸進汽流量的變化、高壓缸軸封至四段抽汽漏汽量和至汽封母管漏汽量隨高壓缸排汽流量的變化以及中壓缸軸封至汽封母管漏汽量隨中壓缸排汽流量的變化基本呈線性。將這些關系曲線擬合成相應的公式用于計算漏汽量。熱平衡圖中高壓缸軸封至汽封回汽母管漏汽量、中壓缸軸封至汽封回汽母管漏汽量、低壓缸軸封供汽量和漏汽量都不隨工況變化而變化，因此這些參數(shù)均取為設計值。

2 變工況計算方法

本文所有計算均通過MATLAB軟件編程實現(xiàn)。當計算某負荷下機組的參數(shù)時，采用定功率模型，即已知機組功率求取機組主蒸汽流量及其他參數(shù)；而當研究某一工況下汽輪機初終參數(shù)變化對熱耗率和功率的影響時，則采用定流量模型，即假定閥門開度不變，機組主蒸汽流量根據(jù)式(3)計算，求取初終參數(shù)變化后對應的機組熱耗率和功率。

定功率變工況模型的計算步驟如下：1）已知機組功率，假定機組主蒸汽流量；2）根據(jù)調節(jié)級焓降與調節(jié)級壓比的關系曲線和式(3)迭代求解調節(jié)級后壓力和焓值，根據(jù)式(2)求解一段抽汽壓力和焓值；3）假定各段抽汽口處蒸汽壓力，求取各段抽汽口處蒸汽焓值，運用加熱器熱平衡，計算各段抽汽流量；4）根據(jù)式(4)求解新的各段抽汽口處和汽缸入口處蒸汽壓力，并與假定值進行對比，如誤差大于10-5，則用新的各抽汽口處和汽缸入口處蒸汽壓力代替假定值重新進行計算，直至前后2次計算誤差小于10-5；5）根據(jù)末級組效率和壓比的關系曲線，求解排汽焓值，然后計算機組新的功率；6）將機組新的功率與已知機組功率比較，如誤差大于10-5，則需重新假定主蒸汽流量，直至計算得到的機組功率和已知功率誤差小于10-5。

定流量變工況模型由于主蒸汽流量已知，與定功率變工況模型的計算步驟相比，無需步驟1）和步驟6），其余步驟相同。

3 變工況模型計算精度驗證

本文以某超臨界660 MW、一次中間再熱、三缸四排汽、單軸、凝汽式CLN660-24.2/566/566 型汽輪機為研究對象。該機組共有8級非調整抽汽，分別供給3臺高壓加熱器、1臺除氧器和4臺低壓加熱器，各加熱器的疏水采用逐級自流的方式。

為驗證改進后的弗留格爾公式的計算精度，分別采用弗留格爾公式和改進后的弗留格爾公式計算了THA工況、75%THA工況和50%THA工況下的調節(jié)級壓力和一段抽汽壓力，計算結果和設計值的對比見表1。

表2為采用定功率變工況模型計算得到的3種工況下主蒸汽流量、各段抽汽壓力和熱耗率與設計值的比較。

表1 弗留格爾公式和改進后的弗留格爾公式計算結果比較

Tab.1 Comparison of calculated results of the Flugel formula and the improved Flugel formula MPa

表2 變工計算結果與設計值對比

Tab.2 Comparison between the calculated results and designed values

從表1可以看出：對于THA工況，由于工況未變化，2種弗留格爾公式的計算結果精度都較高；而對于75%THA工況和50%THA工況，由于工況變化較大，弗留格爾公式的計算結果與設計值的偏差較大，而改進后的弗留格爾公式的計算結果與設計值相比偏差較小。

從表2可以看出，各參數(shù)的最大誤差為2.367%。考慮到熱平衡圖中水和水蒸氣的熱力性質采用IFC67進行計算，而本文則采用IAPWS-IF97，本身就存在一定誤差。由此可見，本文所建立的定功率變工況計算模型可滿足工程計算精度的要求。

4 初終參數(shù)修正曲線的計算結果

以該超臨界660 MW機組汽輪機為例，采用所建立的定流量變工況計算模型，計算得到四閥全開工況下主蒸汽壓力、主蒸汽溫度、再熱蒸汽溫度、再熱壓損、排汽壓力對機組熱耗率和功率的修正曲線，并與制造廠提供的修正曲線進行比較，結果如圖1—圖5所示。

圖2 主蒸汽溫度修正曲線

圖3 再熱蒸汽溫度修正曲線

圖4 再熱壓損修正曲線

圖5 排汽壓力修正曲線

從圖1和圖2可以看出，計算得到的主蒸汽壓力和主蒸汽溫度修正曲線與制造廠提供的修正曲線的變化趨勢一致，且修正率的數(shù)值也基本相同。

從圖3和圖4可以看出，計算得到的再熱蒸汽溫度和再熱壓損的修正曲線與制造廠提供的修正曲線的變化趨勢一致，但再熱蒸汽溫度對熱耗率和再熱壓損對功率和熱耗率的修正率與制造廠提供的修正率存在一定誤差。由于當再熱蒸汽溫度和再熱壓損變化后，機組的高壓缸排汽壓力將發(fā)生變化，進而導致一段抽汽口至高壓缸排汽口之間級組的壓比發(fā)生變化，這可能使這一級組的效率發(fā)生變化。在本文所建立的變工況計算模型中，由于缺乏汽輪機的詳細尺寸，未考慮這一級組效率的變化，可能導致出現(xiàn)一定的計算誤差。

從圖5可以看出，計算得到排汽壓力修正曲線與制造廠提供的修正曲線的變化趨勢一致，但計算得到的修正率絕對值與制造廠提供的修正率絕對值相比偏低。

3 結 論

1）提出了一種根據(jù)機組熱平衡圖構造汽輪機變工況計算模型的簡便方法，采用該方法建立了汽輪機定流量和定功率變工況計算模型。

2）指出了弗留格爾公式應用至汽輪機變工況計算中存在誤差較大的問題，并提出采用改進后的弗留格爾替代弗留格爾公式。

3）以某超臨界660 MW機組為例，采用定功率模型計算得到了典型工況下主蒸汽流量、各段抽汽壓力和熱耗率，與設計值相比，各參數(shù)的最大誤差為2.367%。

4）采用定流量變工況模型計算得到了四閥全開工況下主蒸汽壓力、主蒸汽溫度、再熱蒸汽溫度、再熱壓損和排汽壓力變化對熱耗率和功率的修正曲線，與制造廠提供的修正曲線相比誤差較小。

［1］ 王惠杰, 張春發(fā), 宋之平. 火電機組運行參數(shù)能耗敏感性分析[J]. 中國電機工程學報, 2008, 28(29): 2-10.

WANG Huijie, ZHANG Chunfa, SONG Zhiping. Sensitive analysis of energy consumption of operating parameters for coal-fired unit[J]. Proceedings of the CSEE, 2008, 28(29): 2-10.

［2］ 徐曙. 用Baumann公式估算汽輪機排汽壓力對機組經(jīng)濟性的影響[J]. 汽輪機技術, 2016, 58(1): 72-74.

XU Shu. Estimate exhaust pressure of steam turbine on the economic impact with Baumann formula[J]. Turbine Technology, 2016, 58(1): 72-74.

［3］ 徐大懋, 柯嚴, 王世勇. 汽輪機功率背壓特性的通用計算方法及應用[J]. 熱能動力工程, 2010, 25(6): 605-608.

XU Damao, KE Yan, WANG Shiyong. Universal method for calculating the power/back-pressure characteristics of steam turbine and its application[J]. Journal of Enginee- ring for Thermal Energy Power, 2010, 25(6): 605-608.

［4］ 劉振剛, 李太興, 王俊有, 等. 凝汽機組主蒸汽參數(shù)耗差分析方法的研究及其應用[J]. 熱力透平, 2008, 37(1): 59-62.

LIU Zhengang, LI Taixing, WANG Junyou, et al. Research and application of energy loss analysis of main steam parameter for condensing steam turbine[J]. Thermal Turbine, 2008, 37(1): 59-62.

［5］ 盛德仁, 任浩仁, 陳堅紅, 等. 汽輪機調峰運行時蒸汽初壓的優(yōu)化確定[J]. 動力工程, 2000, 20(5): 809-812.

SHENG Deren, REN Haoren, CHEN Jianhong, et al. Optimization of main steam pressure of steam turbine unit under cycling loading[J]. Power Engineering, 2000, 20(5): 809-812.

［6］ 陳小慧, 高丹, 胡三高. 凝汽機組主蒸汽參數(shù)與熱耗修正曲線的變工況簡便算法[J]. 汽輪機技術, 2014, 56(1): 53-55.

CHEN Xiaohui, GAO Dan, HU Sangao. A simple variable condition calculation method of main Steam parameters and the heat rate crrection curves of the condensate unit[J]. Turbine Technology, 2014, 56(1): 53-55.

［7］ 李慧君, 劉學敏. 基于新型變工況算法的汽輪機初壓優(yōu)化研究[J]. 熱力發(fā)電, 2013, 42(6): 1-5.

LI Huijun, LIU Xuemin. A new variable condition algorithm based study on initial steam pressure optimization for steam turbines[J]. Thermal Power Generation, 2013, 42(6): 1-5.

［8］ 周蘭欣, 華敏, 王為, 等. 機組參數(shù)與熱耗修正曲線的變工況計算法[J]. 動力工程學報, 2011, 31(5): 387-390.

ZHOU Lanxin, HUA Min, WANG Wei, et al. A variable condition calculation method for main steam parameters and the heat rate correction curves[J]. Journal of Chinese Society of Power Engineering, 2011, 31(5): 387-390.

［9］ 方超. 漏汽率試驗時調節(jié)級無溫度測點的解決方法[J]. 熱力透平, 2012, 41(4): 276-279.

FANG Chao. Solution to regulating stage with no temperature measure point during steam leakage rate examination[J]. Thermal Turbine, 2012, 41(4): 276-279.

［10］ 江寧, 曹祖慶. 調節(jié)級通用曲線的計算及應用[J]. 汽輪機技術, 2003, 45(1): 17-19.

JIANG Ning, CAO Zuqing. The calculation and application of turbines’s velocity stage master curve[J]. Turbine Technology, 2003, 45(1): 17-19.

［11］ 徐大懋, 鄧德兵, 王世勇, 等. 汽輪機的特征通流面積及弗留格爾公式改進[J]. 動力工程學報, 2010, 30(7): 473-477.

XU Damao, DENG Debing, WANG Shiyong, et al. Application of characteristic flow area of steam turbines and improvement on Flügel formula[J]. Journal of Chinese Society of Power Engineering, 2010, 30(7): 473-477.

［12］ Performance test codes on steam turbines: ASME PTC6-2004[S].

［13］ 曹麗華, 周云龍, 李勇. 汽輪機通流部分故障診斷基準值的研究[J]. 汽輪機技術, 2009, 51(2): 132-134.

CAO Lihua, ZHOU Yunlong, LI Yong. Study on the norm for fault diagnosis of steam turbine flow passage[J]. Turbine Technology, 2009, 51(2): 132-134.

［14］ 張樂樂, 張燕平, 黃樹紅, 等. 汽輪機濕蒸汽區(qū)蒸汽焓值的計算模型研究[J]. 動力工程學報, 2013, 33(4): 272-277.

ZHANG Lele, ZHANG Yanping, HUANG Shuhong, et al. Steam enthalpy calculation model for wet steam region of steam turbines[J]. Journal of Chinese Society of Power Engineering, 2013, 33(4): 272-277.

［15］ 吳仲. 給水泵汽輪機的運行特性及熱經(jīng)濟性分析[D]. 保定: 華北電力大學, 2008: 13-14.

WU Zhong. Operational characteristic and heat economy analysis of boiler feed water pump turbine[D]. Baoding: North China Electric Power University, 2008: 13-14.

A simple method to construct variable condition model for steam turbine and its application

YU Xinggang1, LI Xu2, JIANG Beihua2, LI Ming2, XIE Danmei3

(1. Hunan Province Key Laboratory of Efficient & Clean Thermal Power Generation Technologies, Changsha 410007, China; 2. Hunan Xiangdian Test & Research Institute Co., Ltd., Changsha 410004, China; 3. School of Power and Mechanical Engineering, Wuhan University, Wuhan 430072, China)

Since the detailed variable condition calculating method of steam turbine is very difficult to be used in the actual calculation process, a simple method to construct variable condition model according to heat balance diagram of the unit was proposed. By employing this method, calculation model of constant flow and constant power variable condition for steam turbines were developed. Moreover, by taking a supercritical 660 MW unit as an example, this model was applied to calculate the values of typical parameters under different typical conditions, and comparison between the calculated values and designed valves indicates the deviation was small. Furthermore, the correcting curves of main steam pressure, main steam temperature, reheat steam temperature, reheater pressure drop and exhaust steam pressure on heat rate and power under VWO condition were also obtained. The difference between the calculated correcting curves and that provided by the manufacturer was small. It demonstrates that the accuracy of the above variable condition model is high and can meet requirements of engineering calculations.

steam turbine, variable condition, heat balance diagram, constant flow, constant power, correcting curve

National Natural Science Foundation of China (51376140); Scientific and Technological Innovation Platform and Talent Plan of Hunan Province (2016TP1027)

TK262

A

10.19666/j.rlfd.201806123

格余興剛, 李旭, 蔣北華, 等. 汽輪機變工況模型的簡便建立方法及應用[J]. 熱力發(fā)電, 2019, 48(2): 96-100. YU Xinggang, LI Xu, JIANG Beihua, et al. A simple method to construct variable condition model for steam turbine and its application[J]. Thermal Power Generation, 2019, 48(2): 96-100.

2018-06-25

國家自然科學基金項目(51376140)；湖南省科技創(chuàng)新平臺與人才計劃(2016TP1027)

余興剛(1989—)，男，博士，工程師，主要研究方向為電廠熱力系統(tǒng)節(jié)能技術，849370151@qq.com。

（責任編輯 劉永強）