王 竹，吳 鵬，張銳鋒，龐文濤，盛德仁

?

基于歷史數(shù)據(jù)的汽輪機調(diào)節(jié)閥流量特性優(yōu)化

王 竹1，吳 鵬2，張銳鋒2，龐文濤1，盛德仁1

（1.浙江大學熱工與動力系統(tǒng)研究所，浙江 杭州 310027； 2.貴州電網(wǎng)有限責任公司電力科學研究院，貴州 貴陽 550000）

汽輪機在經(jīng)過長期運行、汽輪機數(shù)字電液控制系統(tǒng)（DEH）或通流部分改造后，汽輪機調(diào)節(jié)閥的實際流量特性曲線將偏離設(shè)計值，從而影響機組一次調(diào)頻和負荷控制能力，而現(xiàn)場流量特性試驗存在試驗條件苛刻、耗時長、精度低等缺點。本文通過對機組歷史運行數(shù)據(jù)進行分析和挖掘，針對不同的情況提出2種方法以辨識汽輪機調(diào)節(jié)閥實際流量特性，有閥后壓力數(shù)據(jù)時采用迭代計算求出各閥門流量特性曲線，無閥后壓力數(shù)據(jù)時采用弗留格爾公式求出閥組流量特性曲線，分別通過仿真和分段線性化對DEH閥門管理曲線進行優(yōu)化修正，確定閥門合理的重疊度。將該方法應(yīng)用于某660 MW機組，結(jié)果表明在順序閥模式下，汽輪機調(diào)節(jié)閥流量特性曲線的線性度有了很大改善，機組自動發(fā)電控制（AGC）變負荷和一次調(diào)頻的能力得到了提高。

汽輪機；調(diào)節(jié)閥；流量特性；DEH；閥門管理曲線；變負荷；調(diào)頻

汽輪機調(diào)節(jié)閥流量特性是指汽輪機進汽調(diào)節(jié)閥開度與通過調(diào)節(jié)閥的蒸汽流量之間的關(guān)系。在生產(chǎn)過程中，汽輪機長期運行或調(diào)節(jié)閥解體檢修后，調(diào)節(jié)閥的流量特性將發(fā)生改變，在機組自動發(fā)電控制（AGC）變負荷和一次調(diào)頻時，可能出現(xiàn)負荷突變而調(diào)節(jié)緩慢，或者導致汽輪機轉(zhuǎn)子系統(tǒng)劇烈振蕩，嚴重影響機組運行的安全性和穩(wěn)定性[1-3]。準確獲取汽輪機調(diào)節(jié)閥的實際流量特性，使汽輪機調(diào)節(jié)閥總閥位與總流量呈線性關(guān)系，從而實現(xiàn)機組的精確控制與優(yōu)化運行，是汽輪機控制的重要環(huán)節(jié)。

目前，電廠多通過現(xiàn)場流量特性試驗來獲取汽輪機調(diào)節(jié)閥流量特性，但現(xiàn)場流量特性試驗普遍存在試驗條件苛刻、耗時長、精度低等缺點[4-6]。

在機組運行過程中，分布式控制系統(tǒng)（DCS）存儲了大量實際運行數(shù)據(jù)，是機組運行特性的全面、客觀反映。尚星宇等[7]采用DCS數(shù)據(jù)擬合建立蒸汽流量、閥門開度和總閥位指令之間的數(shù)學模型，利用流量校正的方法對閥門開度進行優(yōu)化，最后利用BP神經(jīng)網(wǎng)絡(luò)對優(yōu)化后的閥門流量特性進行模擬還原，結(jié)果表明該方法能有效提高優(yōu)化的可靠性。鄒包產(chǎn)等[8]基于試驗測得的實際調(diào)節(jié)閥流量特性，運用最小二乘法對實測數(shù)據(jù)進行辨識建模，確定最優(yōu)調(diào)節(jié)閥流量特性曲線，并利用BP神經(jīng)網(wǎng)絡(luò)模型預測出修正量，從而對總閥位進行校正。陳堅紅等[9]通過理論計算求出汽輪機調(diào)節(jié)級的臨界壓比，得出調(diào)節(jié)級中全開閥后的噴嘴組最先達到臨界狀態(tài)，噴嘴組壓比始終低于動葉壓比；并且級臨界壓比是僅與調(diào)節(jié)級結(jié)構(gòu)尺寸有關(guān)的常數(shù)，與新蒸汽參數(shù)、機組配汽方式無關(guān)。李存文等[10]采用功率值表征機組進汽量，提出一種基于噴嘴流量計算公式的汽輪機調(diào)節(jié)閥組流量特性辨識方法，優(yōu)化了流量管理函數(shù)，提高了順序閥模式下流量曲線的線性度。王剛等[11]等通過對某660 MW機組汽輪機運行歷史數(shù)據(jù)的挖掘和分析，結(jié)合其特征通流面積的特點，提出了一種新的流量特性曲線辨識方法。

本文對機組歷史運行數(shù)據(jù)進行深入挖掘，對運行數(shù)據(jù)進行預處理，剔除非穩(wěn)態(tài)數(shù)據(jù)和噪聲數(shù)據(jù)，針對機組有無調(diào)節(jié)閥后壓力測點的情況提出了不同的調(diào)節(jié)閥流量特性辨識方法。有閥后壓力測點時采用迭代計算求出各閥門的流量特性曲線，通過MATLAB軟件配汽仿真對DEH閥門管理曲線進行優(yōu)化；無閥后壓力測點時通過弗留格爾公式，求出閥組流量特性曲線，通過分段線性化方法對汽輪機數(shù)字電液控制系統(tǒng)（DEH）閥門管理曲線進行優(yōu)化。以某660 MW機組為例，分別使用2種方法，結(jié)果顯示調(diào)節(jié)閥流量特性曲線的線性度有很大改善，機組AGC變負荷和一次調(diào)頻的能力得到提高。

1 數(shù)據(jù)預處理

1.1 流量特征參數(shù)選取

本文主要研究汽輪機各閥門和閥組的流量特性，因此只需選取與調(diào)節(jié)級閥門及汽輪機流量相關(guān)的參數(shù)即可，包括機組負荷、主蒸汽參數(shù)、調(diào)節(jié)級參數(shù)、高壓缸排汽口（高排）參數(shù)、總閥開度、各閥門開度、各閥后壓力等，具體見表1。

表1 機組調(diào)節(jié)閥流量特征參數(shù)

Tab.1 Flow characteristic parameters of the control valve

1.2 穩(wěn)態(tài)數(shù)據(jù)篩選

2 調(diào)節(jié)閥流量特性辨識

2.1 有閥后壓力測點機組

2.1.1 迭代計算

汽輪機調(diào)節(jié)級各閥門的流量可以由通過其噴嘴組的流量確定

2.1.2 仿真優(yōu)化

通過迭代計算求出調(diào)節(jié)級各個閥門全開度范圍下的流量后，對流量結(jié)果進行標幺化處理，得到各個閥門的流量-開度關(guān)系曲線。由于關(guān)系曲線包含的工況數(shù)據(jù)點較多，為后續(xù)仿真的閥門流量特性函數(shù)的確定造成了困難，運用-means聚類算法對流量-開度關(guān)系曲線進行擬合[14]，擬合出各個閥門的流量特性函數(shù)。采用MATLAB軟件中的仿真軟件Simulink搭建閥門管理模塊，對汽輪機DEH閥門管理過程進行模擬仿真。根據(jù)各個閥門的實際流量特性函數(shù)確定閥門間合理的重疊度，優(yōu)化DEH閥門管理曲線，實現(xiàn)汽輪機調(diào)節(jié)閥配汽運行優(yōu)化。仿真結(jié)構(gòu)流程如圖2所示。

圖1 迭代計算流程

圖2 汽輪機DEH閥門管理曲線仿真結(jié)構(gòu)流程

2.2 無閥后壓力測點機組

2.2.1 弗留格爾公式

弗留格爾公式常用來計算汽輪機級組變工況時的流量特性

式中，、、分別為主蒸汽流量、壓力、溫度，下標0表示額定工況參數(shù)，下標1表示級組前參數(shù)，下標2表示級組后參數(shù)。

弗留格爾公式不適用于蒸汽超音速流動以及閥門切換過程中的流量計算，并且其誤差主要來源于理想氣體假設(shè)[15]。因此對于無閥后壓力測點的機組，將研究對象確定為調(diào)節(jié)級后第一壓力級至高壓缸排氣口間的壓力級組成的級組，并用壓力和比容的乘積代替溫度，從而有效減小誤差[16]，改進后的弗留格爾公式為

式中，下標2表示調(diào)節(jié)級后參數(shù)，下標d表示高壓缸排汽口參數(shù)。

通過式(3)將機組的實際工況流量修正到額定工況流量，并對流量結(jié)果進行標幺化處理，得出額定工況下機組總流量-總閥位關(guān)系曲線，同樣，運用-means聚類算法對總流量-總閥位關(guān)系曲線進行擬合，即得調(diào)節(jié)級閥組流量特性曲線。

2.2.2 分段線性化

總閥位是機組DCS根據(jù)當前主蒸汽壓力及負荷指令等運行參數(shù)通過比例積分產(chǎn)生的計算值，汽輪機各閥門根據(jù)總閥位值確定各自開度，即閥門管理曲線[17]。在定壓運行下，機組負荷與總流量近似成正比關(guān)系，因此，總閥位是連接總流量與各閥開度（流量）的中間值。采用分段線性化的方法對閥組流量特性曲線進行優(yōu)化，通過改變流量拐點和非線性段對應(yīng)的總閥位值，達到改變總閥位與各閥位對應(yīng)關(guān)系的目的，即優(yōu)化了DEH閥門管理曲線，實現(xiàn)了汽輪機調(diào)節(jié)閥配汽運行優(yōu)化。

3 660 MW機組算例及結(jié)果

本算例機組采用N660-24.2/566/566型超臨界、凝汽式汽輪機，調(diào)節(jié)閥為常規(guī)的四閥門配置，開啟順序為CV2/4→CV1→CV3從機組DCS數(shù)據(jù)庫中提取指定參數(shù)連續(xù)3天的運行數(shù)據(jù)，采樣周期為1 s，采集工況數(shù)據(jù)點259 200個。對提取的數(shù)據(jù)進行預處理，剔除非穩(wěn)態(tài)數(shù)據(jù)和噪聲數(shù)據(jù)，得到穩(wěn)態(tài)工況數(shù)據(jù)點9 583個，數(shù)據(jù)精簡效率達到96.3%。穩(wěn)態(tài)工況數(shù)據(jù)庫覆蓋機組運行區(qū)間為：機組負荷321.56~668.75 MW，主蒸汽溫度477.12~566.37 ℃，主蒸汽壓力9.95~25.63 MPa。

3.1 迭代計算

通過Visual Studio平臺的MFC模塊，使用C++語言編寫軟件，實現(xiàn)上述迭代計算和-means聚類算法，軟件主界面如圖3所示。

從穩(wěn)態(tài)工況數(shù)據(jù)庫中篩選出單個閥門全開度下的運行數(shù)據(jù)并存入Excel，所需參數(shù)包括主蒸汽壓力、主蒸汽溫度、閥后壓力、調(diào)節(jié)級壓力和閥門開度。在主界面手動輸入調(diào)節(jié)級結(jié)構(gòu)尺寸參數(shù)（噴嘴出口角1、動靜葉出口面積比、調(diào)節(jié)級平均直徑），然后點擊文件讀取Excel的運行數(shù)據(jù)，對各個閥門全開度范圍下的流量進行計算，最后通過聚類分析完成各個閥門流量特性函數(shù)擬合。該660 MW機組的CV1―CV4流量特性曲線以及與原設(shè)計流量特性曲線的對比如圖4所示。

分析可得CV1－CV4實際流量特性與原設(shè)計流量特性的主要偏差有2點：1）存在死區(qū)，即前5%開度段的流量為0；2）滯后性，特定開度下的實際流量小于原設(shè)計流量。因此必須根據(jù)CV1－CV4實際流量特性重新整定閥門管理函數(shù)，確定閥門間合理的重疊度，通過MATLAB軟件仿真優(yōu)化后閥門管理函數(shù)見表2。

表2 優(yōu)化后閥門管理函數(shù)

Tab.2 The optimized valve management function %

3.2 分段線性化

大多數(shù)時間內(nèi)，機組實際運行在40%~100%額定負荷區(qū)間[18]，較低負荷段的運行數(shù)據(jù)為非穩(wěn)態(tài)數(shù)據(jù)。在無閥后壓力測點的情況下，對穩(wěn)態(tài)工況數(shù)據(jù)庫中的主蒸汽流量應(yīng)用式(3)進行修正，同樣，對結(jié)果進行標幺化處理，并運用-means聚類算法擬合出閥組流量特性曲線，結(jié)果如圖5所示。

由圖5可知：擬合曲線能夠完全表征實際的閥組流量特性，無失真；通過與理想流量特性對比發(fā)現(xiàn)，在總閥位開度為75%~95%時，流量特性曲線出現(xiàn)下凹現(xiàn)象，流量變化較小。該660 MW機組總閥位開度為75%~95%時對應(yīng)CV1、CV3開啟段，因此需要重新設(shè)置CV1和CV3開啟段的開度曲線以及合理的重疊度。采用分段線性化方法修正總閥位指令，優(yōu)化后的閥門管理函數(shù)對比見表3中加黑部分。對比表2和表3可知，2種方法優(yōu)化后的閥門管理函數(shù)基本相同，也驗證了MATLAB軟件仿真優(yōu)化的準確性。

表3 優(yōu)化前后閥門管理函數(shù)對比

Tab.3 The original and optimized valve management functions %

該660 MW機組采用優(yōu)化后的閥門管理函數(shù)實際運行后，提取一段時間內(nèi)的運行數(shù)據(jù)，同理，辨識出汽輪機閥組流量特性，優(yōu)化后的閥組流量特性以及閥門管理曲線如圖6所示。由圖6可見，總閥位開度為75%~95%時的流量下凹現(xiàn)象得到解決，流量特性曲線的線性度得到很大改善。

4 結(jié) 語

本文對機組歷史運行數(shù)據(jù)進行深入挖掘，針對機組有無調(diào)節(jié)閥后壓力測點的情況提出了不同的調(diào)節(jié)閥流量特性辨識與優(yōu)化方法。通過編寫軟件實現(xiàn)上述功能，較好地解決了流量特性試驗復雜、耗時長的問題，可用于指導機組實際運行優(yōu)化。以某660 MW機組為例，2種方法所得結(jié)果基本一致，汽輪機調(diào)節(jié)閥流量特性得到了極大的改善，機組變負荷、一次調(diào)頻能力得到了提高。

［1］ 馬琳, 胥建群, 曹祖慶. 基于量綱分析的汽輪機調(diào)節(jié)級變工況特性研究[J]. 中國電機工程學報, 2014, 34(5): 778-783.

MA Lin, XU Jianqun, CAO Zuqing. Research on varying condition characteristic of steam turbine governing stage based on dimensional analysis[J]. Proceedings of the CSEE, 2014, 34(5): 778-783.

［2］ 李勁柏, 劉復平. 汽輪機閥門流量特性函數(shù)優(yōu)化和對機組安全性經(jīng)濟性的影響[J]. 中國電力, 2008, 41(12): 50-53.

LI Jinbai, LIU Fuping. Optimization for the flow characteristic function of turbine valves and the influences on safe and economic operation[J]. Electric Power, 2008, 41(12): 50-53.

［3］ 李虹銳, 王智微. 基于歷史數(shù)據(jù)計算母管制機組各臺汽輪機進汽流量的研究[J]. 中國電機工程學報, 2016, 36(增刊1): 154-160.

LI Hongrui, WANG Zhiwei. Research on the calculation of the steam rate of turbines based on the historical data in the units with common steam header system[J]. Proceedings of the CSEE, 2016, 36(Suppl.1): 154-160.

［4］ 孫長生, 尹淞, 蘇燁, 等. 電力行業(yè)燃煤機組熱控調(diào)節(jié)系統(tǒng)運行情況調(diào)研[J]. 中國電力, 2013, 46(4): 16-20.

SUN Changsheng, YIN Song, SU Ye, et al. Survey and analysis on thermal control regulation system for coal-fired units in electric power industry[J]. Electric Power, 2013, 46(4): 16-20.

［5］ 鄭國, 張甘泉, 胡晏銘. 300 MW汽輪機主汽調(diào)節(jié)閥組流動特性分析及優(yōu)化改造[J]. 中國電力, 2016, 49(4): 119-123.

ZHENG Guo, ZHANG Ganquan, HU Yanmin. Flow characteristics analysis and optimization of the main steam control valve in 300 MW steam turbine units[J]. Electric Power, 2016, 49(4): 119-123.

［6］ 張寶, 顧正皓, 樊印龍, 等. 通過試驗計算汽輪機的流量特性[J]. 汽輪機技術(shù), 2013, 55(3): 215-218.

ZHANG Bao, GU Zhenghao, FAN Yinlong, et al. Calculating the steam turbine flow characteristics through testing[J]. Turbine Technology, 2013, 55(3): 215-218.

［7］ 尚星宇, 何永君, 王瑞, 等. 基于數(shù)據(jù)擬合的汽輪機調(diào)節(jié)閥流量特性優(yōu)化[J]. 熱力發(fā)電, 2017, 46(3): 121-125.

SHANG Xingyu, HE Yongjun, WANG Rui, et al. Flow characteristics optimization for steam turbine valves based on data fitting[J]. Thermal Power Generation, 2017, 46(3): 121-125.

［8］ 鄒包產(chǎn), 趙宇, 李云, 等. 基于BP神經(jīng)網(wǎng)絡(luò)的汽輪機調(diào)閥流量特性校正[J]. 電力科學與工程, 2017, 33(5): 60-64.

ZOU Baochan, ZHAO Yu, LI Yun, et al. Correction of flow characteristic of steam turbine governing valve based on BP neural network[J]. Electric Power Science and Engineering, 2017, 33(5): 60-64.

［9］陳堅紅, 喬慶, 張寶, 等. 汽輪機調(diào)節(jié)級臨界壓比特性[J]. 浙江大學學報(工學版), 2014, 48(11): 2072-2079.

CHEN Jianhong, QIAO Qing, ZHANG Bao, et al. Critical pressure ratio characteristics of steam turbine governing stage[J]. Journal of Zhejiang University (Engineering Science), 2014, 48(11): 2072-2079.

［10］李存文, 王剛, 梁正玉, 等. 基于噴嘴流量計算公式的汽輪機調(diào)節(jié)閥組流量特性曲線辨識方法研究[C]//2015年中國電機工程學會年會論文集, 武漢: 中國電機工程學會, 2015: 440.

LI Cunwen, WANG Gang, LIANG Zhengyu, et al. Research on the identification and optimization method to steam turbine regulating valve flow characteristic based by nozzle flow calculation[C]//Annual Meeting of China Society of Electrical Engineering in 2015, Wuhan: China Society of Electrical Engineering, 2015: 440.

［11］ 王剛, 李存文, 梁正玉, 等. 基于特征通流面積的汽輪機流量特性辨識方法[J]. 熱力發(fā)電, 2016, 45(6): 33-39.

WANG Gang, LI Cunwen, LIANG Zhengyu, et al. Identification method for flow characteristics of steam turbine based on characteristic flow area[J]. Thermal Power Generation, 2016, 45(6): 33-39.

［12］ 萬祥, 胡念蘇, 韓鵬飛, 等. 大數(shù)據(jù)挖掘技術(shù)應(yīng)用于汽輪機組運行性能優(yōu)化的研究[J]. 中國電機工程學報, 2016, 36(2): 459-467.

WAN Xiang, HU Niansu, HAN Pengfei, et al. Research on application of big data mining technology in performance optimization of steam turbines[J]. Proceedings of the CSEE, 2016, 36(2): 459-467.

［13］ 徐婧, 顧煜炯, 王仲, 等. 基于數(shù)據(jù)挖掘的煤電機組能效特征指標及其基準值的研究[J]. 中國電機工程學報, 2017, 37(7): 2009-2015.

XU Jing, GU Yujiong, WANG Zhong, et al. Research on indexes of energy efficiency and its reference-value for coal-fired power units based on data-mining[J]. Proceedings of the CSEE, 2017, 37(7): 2009-2015.

［14］ 趙歡, 王培紅, 錢瑾, 等. 基于模糊c-均值聚類的鍋爐監(jiān)控參數(shù)基準值建模[J]. 中國電機工程學報, 2011, 31(32): 16-22.

ZHAO Huan, WANG Peihong, QIAN Jin, et al. Modeling for target-value of boiler monitoring parameters based on fuzzy c-means clustering algorithm[J]. Proceedings of the CSEE, 2011, 31(32): 16-22.

［15］ 徐大懋, 鄧德兵, 王世勇, 等. 汽輪機的特征通流面積及弗留格爾公式改進[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]. Power Engineering, 2010, 30(7): 473-477.

［16］ 汪軍, 周建新. 汽輪機主蒸汽流量在線計算方法及應(yīng)用[J]. 熱力發(fā)電, 2010, 39(10): 42-45.

WANG Jun, ZHOU Jianxin. On-line calculation method for the main steam flow rate of steam turbine and its application[J]. Thermal Power Generation, 2010, 39(10): 42-45.

［17］ 高亞杰, 闞玉英. 基于歷史數(shù)據(jù)分析的汽機調(diào)門曲線優(yōu)化[J]. 華北電力技術(shù), 2012(10): 35-39.

GAO Yajie, KAN Yuying. Optimization of turbine governor valve curve based on historical data[J]. North China Electric Power, 2012(10): 35-39.

［18］王剛, 梁正玉, 李存文, 等. 大數(shù)據(jù)分析在汽輪機調(diào)節(jié)門特性參數(shù)辨識及優(yōu)化中的應(yīng)用[J]. 中國電力, 2016, 49(6): 15-19.

WANG Gang, LIANG Zhengyu, LI Cunwen, et al. Application of large data analysis on characteristics parameter identification and optimization of steam turbine regulating valves[J]. Electric Power, 2016, 49(6): 15-19.

Research on flow characteristics of steam turbine regulating valve based on historical data

WANG Zhu1, WU Peng2, ZHANG Ruifeng2, PANG Wentao1, SHENG Deren1

(1. Institute of Thermal Science and Power System, Zhejiang University, Hangzhou 310027, China; 2. Electric Power Research Institute of Guizhou Power Grid Co., Ltd., Guiyang 550000, China)

After long-term operation and retrofit of digital electro-hydraulic (DEH) control system or flow path of steam turbine, the actual flow characteristic curve of the turbine’s regulating valve will deviate from the design values, thus affecting the unit’s primary frequency and load control ability. However, the field flow characteristics test has common shortcomings of severe test conditions, long duration and low precision. In this article, through analysis and excavation of the unit’s historical operation data, two methods are proposed for different situations to identify the actual flow characteristics of the steam turbine’s regulating valve. The flow characteristic curves of each valve are obtained by iteration calculation when the downstream pressure data can be achieved, or by the Freugel’s formula when the downstream pressure data can not be achieved. The DEH valve management curve is optimized by simulation and piece-wise linearization to determine the reasonable valve overlap degree. The application of this method in a 660 MW unit shows that, in sequence valve operation mode, the linearity of the flow characteristics curve has been improved greatly, and the capability of the load variation and primary frequency regulation have been enhanced.

steam turbine, regulating valve, flow characteristics, DEH, valve management curve, variable load, frequency modulation

Science and Technology Project of Guizhou Power Grid Co., Ltd. (GZ2015-2-0091)

王竹(1995—)，男，碩士，主要研究方向為火電機組運行優(yōu)化，50439923@qq.com。

TK37

B

10.19666/j.rlfd.201804079

王竹, 吳鵬, 張銳鋒, 等. 基于歷史數(shù)據(jù)的汽輪機調(diào)節(jié)閥流量特性優(yōu)化[J]. 熱力發(fā)電, 2019, 48(2): 39-44. WANG Zhu, WU Peng, ZHANG Ruifeng, et al. Research on flow characteristics of steam turbine regulating valve based on historical data[J]. Thermal Power Generation, 2019, 48(2): 39-44.

2018-04-19

貴州電網(wǎng)有限責任公司科技項目(GZ2015-2-0091)

盛德仁(1960—)，男，教授，shengdr@zju.edu.cn。

（責任編輯 杜亞勤）