王凱潤, 馬素霞, 霍 茂

(太原理工大學(xué) 電氣與動力工程學(xué)院, 太原 030024)

火電廠運(yùn)行性能在線監(jiān)測是節(jié)能提效降碳[1]的基本保障。汽輪機(jī)性能指標(biāo)如熱耗率、相對熱效率等的計(jì)算都需要確定汽輪機(jī)的排汽狀態(tài),但因排汽處于濕蒸汽區(qū),僅通過監(jiān)測其壓力或溫度無法得到排汽焓值[2],使得汽輪機(jī)的運(yùn)行性能精準(zhǔn)在線計(jì)算難以實(shí)現(xiàn)。

近年來,很多學(xué)者進(jìn)行了低壓缸排汽焓在線計(jì)算方法的研究。任浩仁等[3]使用能量平衡法計(jì)算汽輪機(jī)排汽焓,并通過曲線擬合確定排汽焓迭代初值。韓中合等[4]將包括汽輪機(jī)及凝汽器在內(nèi)的回?zé)嵯到y(tǒng)視為一閉口系統(tǒng),使用能量平衡法計(jì)算汽輪機(jī)排汽焓。張赟等[5]根據(jù)濕氣損失與排汽焓的關(guān)系提出一種計(jì)算方法,但計(jì)算過程中涉及到的一些初始參數(shù)一般情況下難以獲得。齊繼鵬等[6]提出利用偏最小二乘法和多種群遺傳算法優(yōu)化極限學(xué)習(xí)機(jī)的汽輪機(jī)排汽焓預(yù)測模型。楊斌等[7]建立了基于主成分分析與神經(jīng)網(wǎng)絡(luò)復(fù)合模型的汽輪機(jī)排汽焓計(jì)算模型,并采用汽輪機(jī)實(shí)際運(yùn)行數(shù)據(jù)進(jìn)行訓(xùn)練,形成了便于在線監(jiān)測實(shí)時調(diào)用的模型。

綜上所述,目前求取低壓缸排汽焓的方法主要有能量平衡法、曲線外推法、級內(nèi)損失法、數(shù)據(jù)分析法等。其中,能量平衡法和曲線外推法屬于熱力學(xué)方法,其特點(diǎn)在于易于計(jì)算,普遍性強(qiáng),但計(jì)算精度不高。級內(nèi)損失法則需已知大量汽輪機(jī)結(jié)構(gòu)數(shù)據(jù),這些數(shù)據(jù)通常難以獲取,但該類方法針對性強(qiáng),計(jì)算精度較高。數(shù)據(jù)分析法主要分為黑箱子法與機(jī)理法,其計(jì)算精度取決于數(shù)據(jù)集的準(zhǔn)確性。黑箱子法未分析輸入與輸出變量之間的關(guān)系,僅將其作為數(shù)學(xué)問題求解,且訓(xùn)練計(jì)算模型時往往需要大量的數(shù)據(jù),通常情況下,只有少部分?jǐn)?shù)據(jù)為試驗(yàn)數(shù)據(jù),大部分?jǐn)?shù)據(jù)為根據(jù)能量平衡求出的排汽焓,因此難以擺脫能量平衡法誤差累積的影響。

筆者基于汽輪機(jī)變工況原理,以反映排汽焓與做功能量轉(zhuǎn)換能效關(guān)系的多變指數(shù)為中間參數(shù),將不能直接測量的排汽焓同能夠直接測量的溫度、壓力等數(shù)據(jù)相關(guān)聯(lián),提出了一種在線計(jì)算汽輪機(jī)排汽焓的新方法——多變指數(shù)法。通過對額定功率為300～1 000 MW的8臺汽輪機(jī)的設(shè)計(jì)數(shù)據(jù)和運(yùn)行數(shù)據(jù)進(jìn)行分析,得出了不同容量汽輪機(jī)低壓缸多變指數(shù)與低壓缸進(jìn)汽參數(shù)和排汽參數(shù)的關(guān)聯(lián)式,為排汽焓在線計(jì)算提供技術(shù)方法。多變指數(shù)法本質(zhì)上是對汽輪機(jī)熱平衡系統(tǒng)的充分分析、提取及外推,所需初始參數(shù)簡單、易得,計(jì)算過程不涉及迭代,在線性能好,且計(jì)算精度能夠滿足工程需要。

1 低壓缸多變指數(shù)與汽輪機(jī)排汽焓關(guān)系

1.1 多變指數(shù)

蒸汽在噴嘴和動葉中的膨脹過程是熵增的多變過程,多變指數(shù)n隨級內(nèi)損失的增大而減小。以噴嘴為例,噴嘴膨脹多變指數(shù)n隨噴嘴速度系數(shù)φ的減小而減小。n隨φ2的變化規(guī)律如圖1所示。汽輪機(jī)低壓缸末幾級工作在濕蒸汽區(qū),當(dāng)濕汽損失增大,級的多變指數(shù)n減小。總之,工況偏離設(shè)計(jì)工況越遠(yuǎn),汽輪機(jī)各級內(nèi)的流動損失越大,葉柵速度系數(shù)越小,多變指數(shù)越小,故可以用多變指數(shù)來反映汽輪機(jī)的運(yùn)行工況。

圖1 n隨φ2的變化規(guī)律Fig.1 Variation of n with φ2

蒸汽在汽輪機(jī)低壓缸內(nèi)的做功可認(rèn)為是一個多變過程[8],則

(1)

(2)

式中:p1為低壓缸進(jìn)汽壓力,MPa;v1為進(jìn)汽比容,m3/kg;p2為低壓缸排汽壓力,MPa;v2為排汽比容,m3/kg。

式(1)中共有5個變量,其中:壓力p1、p2可通過測量得到;v1可根據(jù)入口過熱蒸汽壓力和溫度計(jì)算得到;由于低壓缸排汽常為濕蒸汽,排汽狀態(tài)點(diǎn)難以確定,故v2為未知量;多變指數(shù)n既與進(jìn)汽參數(shù)有關(guān)又受排汽參數(shù)的影響,因此不能直接求出。

汽輪機(jī)設(shè)計(jì)熱平衡圖中含有不同工況下的低壓缸排汽焓值,該值是汽輪機(jī)設(shè)計(jì)廠家以熱力試驗(yàn)為基礎(chǔ),結(jié)合汽輪機(jī)的結(jié)構(gòu)參數(shù)通過級內(nèi)損失計(jì)算給出,準(zhǔn)確可靠。以這些值為樣本數(shù)據(jù),并考慮到影響排汽焓值的主要因素,依據(jù)多元回歸方法[9]分析多變指數(shù)的變化規(guī)律,根據(jù)低壓缸不同工況下的多變指數(shù)計(jì)算排汽焓是實(shí)現(xiàn)其在線計(jì)算的有效途徑之一。

1.2 多變指數(shù)與排汽焓關(guān)系

首先,通過數(shù)據(jù)分析建立多變指數(shù)n隨汽輪機(jī)組低壓缸進(jìn)汽量G、進(jìn)汽溫度T、進(jìn)汽壓力p1及排汽壓力p2的關(guān)聯(lián)式:

n=f(G,T,p1,p2)+ε

(3)

式中:ε為隨機(jī)誤差。

然后,通過在線測得的G、p1、p2、T實(shí)時計(jì)算出口排汽比容v2,即

(4)

最后,由排汽比容v2及排汽壓力p2計(jì)算排汽比焓h2。

h2=f(p2,v2)

(5)

2 多變指數(shù)建模

2.1 變量選取及數(shù)據(jù)分析

以某電廠一臺CZK350/295-24.2/0.4/566/566型汽輪機(jī)為研究對象,對其全工況熱平衡圖中相關(guān)數(shù)據(jù)進(jìn)行分析計(jì)算,得到不同工況下的排汽焓和低壓缸多變指數(shù)。低壓缸的主要熱力參數(shù)有低壓缸進(jìn)、出口蒸汽參數(shù)以及低壓缸各級回?zé)岢槠麉?shù)。低壓缸回?zé)岢槠麨榉钦{(diào)整式抽汽,抽汽量與回?zé)岢槠麎毫σ约暗蛪焊走M(jìn)汽量具有很強(qiáng)的正相關(guān)性,即各抽汽量與低壓缸進(jìn)汽量成比例變化[10]。

綜上,回?zé)岢槠麉?shù)與低壓缸進(jìn)汽參數(shù)具有多重共線性,在數(shù)據(jù)分析的過程中,具有多重共線性的多個變量選用其中一個變量進(jìn)行分析即可。因此,選取低壓缸多變指數(shù)n與低壓缸進(jìn)汽量G、低壓缸進(jìn)汽溫度T、低壓缸進(jìn)汽壓力p1及背壓pb等參數(shù)進(jìn)行分析。多變指數(shù)n隨低壓缸進(jìn)、出口參數(shù)的變化結(jié)果如圖2所示。

(a)

(b)

(c)

(d)圖2 CZK350機(jī)組低壓缸多變指數(shù)n隨低壓缸進(jìn)、出口參數(shù)的變化Fig.2 Fig.2 Variation of polytropic exponent n with the inlet and exhaust parameters of LP turbine for CZK350 unit

由圖2可知,低壓缸進(jìn)汽量G對多變指數(shù)n具有顯著的負(fù)向影響;低壓缸進(jìn)汽溫度T對n具有顯著的正向影響;背壓pb對n具有顯著的正向影響。G對n呈對數(shù)函數(shù)關(guān)系變化,這一點(diǎn)在設(shè)計(jì)背壓工況及供熱工況時有明顯的體現(xiàn)。設(shè)計(jì)背壓工況下,機(jī)組運(yùn)行背壓均為11 kPa,不同工況下的低壓缸進(jìn)汽溫度略有差異,導(dǎo)致各點(diǎn)沒有嚴(yán)格地落在一條對數(shù)函數(shù)曲線上。相同進(jìn)汽量下,供熱工況下的n相比于設(shè)計(jì)背壓工況下的n有所減小。造成該現(xiàn)象的原因有兩方面:第一,供熱工況下平均低壓缸進(jìn)汽溫度為218 ℃,相比設(shè)計(jì)背壓工況下的降低了27 K;第二,供熱工況下運(yùn)行背壓為9 kPa,相比設(shè)計(jì)背壓工況下降低了2 kPa。變背壓工況下,機(jī)組低壓缸進(jìn)汽量與進(jìn)汽溫度均變動不大,n隨pb升高而逐漸增大,呈對數(shù)函數(shù)關(guān)系。變背壓工況包含滿發(fā)背壓工況、銘牌工況(TRL工況)及阻塞背壓工況,背壓變動范圍為5.8～45.0 kPa。

2.2 精確模型構(gòu)建

采用多元回歸分析方法,建立低壓缸多變指數(shù)n與低壓缸進(jìn)、出口參數(shù)的定量關(guān)系。依據(jù)表1中的15組數(shù)據(jù),采用交叉驗(yàn)證法中的留一法[11]劃分出建立模型所需的訓(xùn)練集與測試集,同時對各變量進(jìn)行無量綱化與線性化變量替換[12]。多變指數(shù)的多元線性回歸模型為:

(6)

式中:G0為熱耗率驗(yàn)收工況(THA工況)下低壓缸進(jìn)汽量,t/h;pb0為THA工況下的運(yùn)行背壓,MPa;T0為THA工況下低壓缸進(jìn)汽溫度,℃。

表1 CZK350機(jī)組低壓缸進(jìn)、排汽參數(shù)及多變指數(shù)

2.3 模型檢驗(yàn)

求出回歸模型后,采用擬合優(yōu)度檢驗(yàn)、方程總體的顯著性檢驗(yàn)判定模型的可靠度。模型的擬合優(yōu)度由可決系數(shù)R2判定,R2越接近1,模型的擬合優(yōu)度越高;使用F檢驗(yàn)驗(yàn)證模型的整體顯著性水平,當(dāng)F統(tǒng)計(jì)量的概率值P<0.001時判定模型通過檢驗(yàn)[13]。

對式(6)所示的回歸模型進(jìn)行統(tǒng)計(jì)檢驗(yàn),檢驗(yàn)結(jié)果見表2。由表2可得,F統(tǒng)計(jì)量的概率值為5.43×10-14,小于0.001,可認(rèn)為自變量和因變量之間存在顯著的線性關(guān)系。從擬合優(yōu)度角度進(jìn)行分析,該回歸模型的R2=0.998,表明方程擬合優(yōu)度較好。

由模型計(jì)算得到的不同工況下低壓缸多變指數(shù)和排汽焓,與多變指數(shù)和排汽焓設(shè)計(jì)值的對比結(jié)果分別如圖3、圖4所示。

表2 多元線性回歸模型擬合優(yōu)度及顯著性檢驗(yàn)

圖3 CZK350機(jī)組低壓缸多變指數(shù)模型計(jì)算誤差Fig.3 Calculation error of polytropic exponential model of LP turbine for CZK350 unit

圖4 CZK350機(jī)組低壓缸排汽焓模型計(jì)算誤差Fig.4 Calculation error of exhaust enthalpy model of LP turbine for CZK350 unit

由圖3可知,多變指數(shù)模型計(jì)算值的平均相對誤差為0.078%,最大相對誤差為0.36%。由圖4可知,排汽焓計(jì)算值與設(shè)計(jì)值的平均相對誤差為0.165%,最大相對誤差為0.738%,具有較高的準(zhǔn)確度。

3 汽輪機(jī)排汽焓通用計(jì)算模型

3.1 通用模型

式(6)為以CZK350機(jī)組汽輪機(jī)熱平衡數(shù)據(jù)建立的排汽焓模型,該模型對CZK350機(jī)組的排汽焓計(jì)算精度較高,對其余機(jī)組的排汽焓計(jì)算精度較差。因此對8臺不同容量機(jī)組的熱平衡系統(tǒng)進(jìn)行分析,擬構(gòu)建排汽焓通用計(jì)算模型。8臺機(jī)組的編號及型號見表3。

表3 8臺機(jī)組的編號及型號

8臺機(jī)組的建模數(shù)據(jù)共計(jì)129組,采用十折交叉驗(yàn)證法[14]劃分出建立模型所需的訓(xùn)練集與測試集。各工況下低壓缸多變指數(shù)n與低壓缸相對進(jìn)汽量、進(jìn)汽溫度和排汽壓力的關(guān)系如圖5所示,圖5中氣泡的大小用于區(qū)別低壓缸多變指數(shù)的大小。其多元線性回歸通用模型如式(7)所示。

(7)

圖5 機(jī)組A～H低壓缸多變指數(shù)n隨低壓缸進(jìn)出口參數(shù)的變化Fig.5 Variation of A-H units' LP polytropic exponent n with inlet and outlet parameters

3.2 通用模型檢驗(yàn)及誤差分析

對式(7)的通用模型進(jìn)行統(tǒng)計(jì)檢驗(yàn),該通用模型的R2=0.413,表明方程擬合優(yōu)度較差。

這是由于利用式(7)計(jì)算機(jī)組的排汽焓時采用了該機(jī)組THA工況下的設(shè)計(jì)數(shù)據(jù)作為公式中的pb0和T0。不同機(jī)組THA工況下低壓缸進(jìn)汽溫度T0與設(shè)計(jì)背壓pb0不盡相同,將T與pb同各機(jī)組設(shè)計(jì)值相除進(jìn)行無量綱化不合理。因此,引入基準(zhǔn)背壓,以4 kPa作為改進(jìn)通用模型的基準(zhǔn)背壓。類似地,以200 ℃作為改進(jìn)通用模型的低壓缸基準(zhǔn)進(jìn)汽溫度。

3.3 改進(jìn)通用模型

采用基準(zhǔn)背壓和基準(zhǔn)進(jìn)汽溫度對自變量進(jìn)行無量綱化,以建立改進(jìn)通用模型。改進(jìn)后的通用模型為:

(8)

式中:ps為基準(zhǔn)背壓,kPa;Ts為低壓缸基準(zhǔn)進(jìn)汽溫度,℃。

對式(8)進(jìn)行統(tǒng)計(jì)檢驗(yàn),結(jié)果見表4。

表4 改進(jìn)通用模型擬合優(yōu)度及顯著性檢驗(yàn)

改進(jìn)通用模型的R2=0.942,方程擬合優(yōu)度較好。圖6為通用模型與改進(jìn)通用模型計(jì)算誤差的對比。

圖6 通用模型與改進(jìn)通用模型計(jì)算誤差對比Fig.6 Comparison of calculation errors between general model and improved general model

由圖6可知,通用模型得到的排汽焓計(jì)算值與設(shè)計(jì)值的平均相對誤差為3.42%,最大相對誤差為8.57%;經(jīng)改進(jìn)后平均相對誤差降為0.93%,最大相對誤差為3.52%,可滿足工程計(jì)算精度要求。

3.4 低負(fù)荷及關(guān)閉回?zé)岢槠r下改進(jìn)通用模型誤差分析

2.1節(jié)進(jìn)行變量選取時忽略了低壓缸內(nèi)回?zé)岢槠挠绊懸韵鳒p變量的個數(shù),現(xiàn)對此項(xiàng)簡化所帶來的影響進(jìn)行分析。在所分析的129個工況中,共有11個工況關(guān)閉了低壓缸回?zé)岢槠?。由改進(jìn)通用模型得到的該11個工況下排汽焓計(jì)算值與設(shè)計(jì)值的平均相對誤差為0.46%,最大相對誤差為1.21%,如圖7所示。此平均誤差相較于整體誤差下降了0.47%,因此即使因?yàn)楣r偏離太遠(yuǎn)導(dǎo)致低壓缸不能抽汽,改進(jìn)通用模型的計(jì)算精度仍能夠滿足工程要求。

圖7 關(guān)閉回?zé)岢槠r下改進(jìn)通用模型排汽焓計(jì)算誤差Fig.7 Calculation error of exhaust steam enthalpy under heaters cut condition by improved general model

在129個工況中共有16個工況汽輪機(jī)處于40%負(fù)荷及以下的低負(fù)荷運(yùn)行狀態(tài),此時的排汽焓計(jì)算誤差如圖8所示。

圖8 低負(fù)荷工況下改進(jìn)通用模型排汽焓計(jì)算誤差Fig.8 Calculation error of exhaust steam enthalpy under low load condition by improved general model

低負(fù)荷下,由改進(jìn)通用模型得到的計(jì)算值與設(shè)計(jì)值的排汽焓平均相對誤差為1.01%,最大相對誤差為3.52%。低負(fù)荷工況下排汽焓計(jì)算值與設(shè)計(jì)值的平均誤差相比于總體的平均誤差稍大,且包含了相對誤差最大的工況點(diǎn),但是,5%以下的相對誤差仍在工程誤差允許范圍之內(nèi)。

4 結(jié) 論

(1) 汽輪機(jī)排汽焓可通過實(shí)時確定低壓缸做功多變指數(shù)實(shí)現(xiàn)在線計(jì)算。汽輪機(jī)低壓缸多變指數(shù)與低壓缸進(jìn)汽量、進(jìn)汽溫度、進(jìn)汽壓力、排汽壓力等密切相關(guān)。

(2) 使用多變指數(shù)法建立了針對CZK350機(jī)組運(yùn)行特性的排汽焓在線計(jì)算模型。由模型得到的排汽焓計(jì)算值與設(shè)計(jì)值的平均相對誤差為0.165%,最大相對誤差為0.738%。

(3) 建立了適用于不同機(jī)組排汽焓計(jì)算的改進(jìn)通用模型,其排汽焓計(jì)算值與設(shè)計(jì)值的平均相對誤差為0.93%,最大相對誤差為3.52%,可滿足工程計(jì)算精度要求。