自然通風(fēng)冷卻塔預(yù)熱式防凍系統(tǒng)控制研究

金慶輝， 宋傳教， 楊建明， 黃恩和， 姚麗萍

(1. 東南大學(xué) 能源與環(huán)境學(xué)院， 南京 210096； 2. 江蘇省國(guó)信集團(tuán)有限公司， 南京 210005)

自然通風(fēng)冷卻塔是汽輪發(fā)電機(jī)組閉式循環(huán)冷端系統(tǒng)的重要設(shè)備，因其運(yùn)行時(shí)具有對(duì)周邊生態(tài)環(huán)境影響小、運(yùn)行費(fèi)用及故障率低、易于維護(hù)等優(yōu)點(diǎn)，而被廣泛應(yīng)用于電力行業(yè)[1-2]，并且逐漸取代從江河直接取水的開式供水冷卻方式。凝汽器出口的高溫循環(huán)水進(jìn)入冷卻塔的中央豎井后，由配水系統(tǒng)分配至配水槽，經(jīng)若干個(gè)噴淋裝置噴濺成小水滴，均勻?yàn)R射在填料表面形成水膜，與干冷空氣充分接觸?？諝馕昭h(huán)水的熱量后，溫度升高、濕度變大、密度變小，濕熱空氣由冷卻塔出風(fēng)口排向大氣，冷卻后的循環(huán)水匯集于集水池，由循環(huán)水泵返送至凝汽器[3-5]。

寒冷季節(jié)，部分地區(qū)氣溫降至零下十幾攝氏度甚至零下幾十?dāng)z氏度，淋水與冷空氣直接接觸，水溫易降至冰點(diǎn)以下，在進(jìn)風(fēng)口上緣、支撐梁柱、填料或噴淋裝置等位置會(huì)出現(xiàn)結(jié)冰[6-7]，并加劇成長(zhǎng)形成冰凌、冰柱。這不僅會(huì)阻擋進(jìn)風(fēng)，影響冷卻效果，惡化氣水傳熱，更嚴(yán)重的是會(huì)增大填料與支撐構(gòu)件的載荷，甚至造成坍塌事故，威脅冷卻塔的安全運(yùn)行；另外，混凝土的反復(fù)凍融也會(huì)使塔體結(jié)構(gòu)受損，降低使用壽命[8]。

目前，電廠較多采取冷卻塔進(jìn)風(fēng)口加裝擋風(fēng)圍板的防凍方式，依靠降低通風(fēng)量防止淋水結(jié)冰[8]。實(shí)際運(yùn)行時(shí)，由于擋風(fēng)板安裝高度和面積較大，裝卸具有一定難度，開度不可調(diào)，雖然滿足了最低溫度時(shí)的防凍要求，避免了支承結(jié)構(gòu)載荷過(guò)大等安全問(wèn)題；但是，其他工況時(shí)因過(guò)度防凍而導(dǎo)致出塔水溫高于經(jīng)濟(jì)值，冷卻塔喪失主要冷卻功能，甚至造成凝汽器壓力接近夏秋高溫季時(shí)的壓力，導(dǎo)致機(jī)組發(fā)電效率降低，熱耗及煤耗增加，冷端經(jīng)濟(jì)性較差。

潘昌遠(yuǎn)[9]針對(duì)300 MW燃煤機(jī)組研究了擋風(fēng)板設(shè)置層數(shù)對(duì)機(jī)組煤耗的影響，計(jì)算得到冬季因懸掛過(guò)多的擋風(fēng)板而多消耗標(biāo)準(zhǔn)煤3 802.13 t，造成235.24萬(wàn)元經(jīng)濟(jì)損失。李永華等[10]針對(duì)北方寒冷季填料下表面及進(jìn)風(fēng)口易結(jié)冰的問(wèn)題，研究了擋風(fēng)面積對(duì)不同特征面水滴溫度的影響，得出當(dāng)環(huán)境溫度為-10 ℃、-17 ℃、-23 ℃時(shí)，分別加裝1層、3層、5層擋風(fēng)板，內(nèi)部各特征面的水滴平均溫度均已達(dá)到冰點(diǎn)以上，能剛好防止塔內(nèi)結(jié)冰；但是，由于進(jìn)風(fēng)區(qū)域與大氣環(huán)境直接接觸，特征面呈現(xiàn)出中間高、周圍低的溫度分布，塔內(nèi)靠近進(jìn)風(fēng)邊緣仍存在氣溫低于冰點(diǎn)的區(qū)域，甚至接近環(huán)境溫度，無(wú)法完全規(guī)避結(jié)冰風(fēng)險(xiǎn)。溫智慧[11]借助FLUENT軟件研究了不同擋風(fēng)板布置方式的防凍效果，得出縱向間隔懸掛擋風(fēng)板優(yōu)于迎風(fēng)側(cè)懸掛方式，并且前者需要更小的擋風(fēng)面積。多數(shù)研究集中在對(duì)不同擋風(fēng)板布置方式及懸掛層數(shù)的防凍效果研究，并且未解決變工況運(yùn)行下因擋風(fēng)板開度不可調(diào)而導(dǎo)致經(jīng)濟(jì)性差等問(wèn)題。目前，對(duì)冷卻塔的防凍仍處于不可控階段。

筆者對(duì)自然通風(fēng)冷卻塔提出一種預(yù)熱式防凍系統(tǒng)，在冷卻塔底部進(jìn)風(fēng)口設(shè)置表面式氣水空氣預(yù)熱器，利用凝汽器出口高溫循環(huán)水通過(guò)表面換熱加熱冷空氣，在防凍約束條件下，基于氣象條件及機(jī)組負(fù)荷變化控制進(jìn)風(fēng)量和加熱循環(huán)水量，使出塔水溫盡可能低、機(jī)組冷端經(jīng)濟(jì)性最優(yōu)；同時(shí)，基于防凍系統(tǒng)的組成和工藝流程，對(duì)循環(huán)水和冷卻空氣的流動(dòng)與傳熱建立數(shù)學(xué)模型，開發(fā)1 000 MW機(jī)組冷卻塔預(yù)熱式防凍系統(tǒng)的仿真計(jì)算軟件，驗(yàn)證防凍系統(tǒng)的有效性，在防凍約束下優(yōu)化出塔水溫，給出最優(yōu)控制參數(shù)，為電廠寒冷季節(jié)冷卻塔安全運(yùn)行及冷端優(yōu)化提供參考。

1 預(yù)熱式防凍系統(tǒng)

冷卻塔預(yù)熱式防凍系統(tǒng)的基本工藝流程[12]見圖1。塔底局部進(jìn)風(fēng)區(qū)域全周覆蓋擋風(fēng)板，未覆蓋區(qū)域均開設(shè)百葉窗，并設(shè)置利用凝汽器出口高溫循環(huán)水加熱進(jìn)風(fēng)的翅片管式氣水換熱系統(tǒng)。原中央豎井出口的循環(huán)水分為兩路：一路仍沿進(jìn)水母管進(jìn)入冷卻塔配水系統(tǒng)；另一路進(jìn)入預(yù)熱器管側(cè)。

圖1 冷卻塔預(yù)熱式防凍系統(tǒng)的基本工藝流程

自然通風(fēng)冷卻塔最易結(jié)冰的區(qū)域?yàn)檫M(jìn)風(fēng)口邊緣，冷卻塔進(jìn)風(fēng)口冰幕現(xiàn)場(chǎng)圖見圖2。

圖2 冷卻塔進(jìn)風(fēng)口冰幕現(xiàn)場(chǎng)圖

預(yù)熱器的設(shè)置使塔進(jìn)風(fēng)得到充分預(yù)熱，平均溫度升高，降低了低溫運(yùn)行工況下淋水與冷空氣因接觸散熱造成的結(jié)冰風(fēng)險(xiǎn)，同時(shí)可以克服擋風(fēng)板防凍方式下進(jìn)塔風(fēng)溫仍接近實(shí)際環(huán)境溫度的缺陷。由調(diào)節(jié)閥調(diào)節(jié)進(jìn)水量，以及由百葉窗調(diào)節(jié)進(jìn)風(fēng)量，管側(cè)循環(huán)水經(jīng)冷卻后匯聚于集水池，承擔(dān)了冷卻塔部分冷卻負(fù)荷。

防凍系統(tǒng)設(shè)有溫度監(jiān)測(cè)及控制模塊，其運(yùn)行狀態(tài)可調(diào)可控，可依據(jù)氣象條件或機(jī)組負(fù)荷變化調(diào)節(jié)預(yù)熱器的冷熱流體量，嚴(yán)格控制進(jìn)塔風(fēng)溫與出塔水溫。另外，在非寒冷季節(jié)運(yùn)行時(shí)，考慮到預(yù)熱器阻力對(duì)冷卻塔冷卻性能的影響，應(yīng)及時(shí)拆卸防凍裝置。

2 數(shù)學(xué)模型及計(jì)算方法

2.1 預(yù)熱器模型

2.1.1 預(yù)熱器熱力模型

預(yù)熱器由多支翅片管、管箱及構(gòu)架組成，翅片管是核心傳熱元件，管箱作為連接翅片管兩端的箱體，構(gòu)架用以支撐和固定翅片管束。預(yù)熱器管束排列設(shè)置為管中心呈三角形的叉排布置，熱力性能計(jì)算采用的試驗(yàn)關(guān)聯(lián)式[13]為：

(1)

式中：ho為氣側(cè)對(duì)流傳熱系數(shù)，W/(m2·K)；λ為導(dǎo)熱系數(shù)，W/(m·K)；dr為翅根直徑，m；qm,max為最窄截面處的質(zhì)量流量，kg/(m2·s)；Sf、Hf分別為翅片間距和翅片高度，m；μ為動(dòng)力黏度，kg/(m2·s)；cp為比定壓熱容，J/(kg·K)。

預(yù)熱器的氣水流動(dòng)形式是以順流為總趨勢(shì)的多次交叉流，具體連接方式和溫度變化見圖3，其中：t1,in、t1,out分別為預(yù)熱器內(nèi)循環(huán)水的進(jìn)口、出口溫度；t2,in、t2,out分別為預(yù)熱器內(nèi)空氣的進(jìn)口、出口溫度；t1,i、t2,i分別為循環(huán)水和空氣在管程i的溫度。

圖3 預(yù)熱器的流動(dòng)形式及溫度變化

熱力模型求解方法采用效能-傳熱單元數(shù)法[14]，具體計(jì)算公式為：

ε1=(t2,in-t2,out)/(t1,in-t2,in)

(2)

ε2=[1-e-NTU(1-C)φ]/[1-Ce-NTU(1-C)φ]

(3)

式中：ε1、ε2均為效能；NTU為傳熱單元數(shù)；Cmin、Cmax分別為最小、最大質(zhì)量熱容流率，W/K；K為總傳熱系數(shù)，W/(m2·K)；Atotal為總傳熱面積，m2；n為總管程數(shù)；φ、φi分別為總體和管程i的溫差修正(相比于純逆流的修正)系數(shù)；P、Pi分別為總體和管程i的熱效率；R、Ri分別為總體和管程i的熱容比，求解模型時(shí)假設(shè)了Ri=R。

2.1.2 預(yù)熱器阻力模型

預(yù)熱器的氣流阻力計(jì)算公式[15]為：

式中：Δpz2為預(yù)熱器氣流阻力，Pa；ρ為空氣密度，kg/m3；f為阻力系數(shù)；S1、S2分別為翅片管的橫向、縱向間距，m；ρ為濕空氣密度，kg/m3。

2.2 冷卻塔模型

2.2.1 冷卻塔熱力性能模型

冷卻塔熱力性能計(jì)算采用Mekel模型[16]，具體計(jì)算公式為：

N2=mApγn

(6)

根據(jù)濕空氣飽和焓近似假定下的平均焓差法求解計(jì)算，具體計(jì)算公式為：

(7)

2.2.2 冷卻塔阻力及抽力模型

冷卻塔阻力主要集中在塔進(jìn)出口、雨區(qū)、淋水填料、配水系統(tǒng)、收水層，其阻力方程為：

式中：Δpz1為冷卻塔氣流阻力，Pa；vf為填料斷面風(fēng)速，m/s；ξ1、ξ2均為阻力系數(shù)[17]；A0、A1分別為填料、出風(fēng)口的面積，m2。

冷卻塔依靠塔內(nèi)外空氣密度差形成抽力，其計(jì)算公式[18]為：

Δpd=Htowerg(ρin-ρout)

(11)

式中：Δpd為抽力，Pa；Htower為冷卻塔有效高度，m；g為重力加速度，m/s2；ρin、ρout分別為進(jìn)塔、出塔空氣的密度，kg/m3。

2.3 循環(huán)水泵變速調(diào)節(jié)模型

變速調(diào)節(jié)是通過(guò)改變循環(huán)水泵轉(zhuǎn)速來(lái)改變運(yùn)行工況，達(dá)到調(diào)節(jié)水量目的。由相似定律得到變轉(zhuǎn)速比下的特性方程為：

式中：H為實(shí)際揚(yáng)程，m；P為功率，kW；qV為容積流量，m3/s；k為變轉(zhuǎn)速比；a1、a2、b1、b2、b3均為擬合系數(shù)。

2.4 仿真計(jì)算平臺(tái)開發(fā)

基于LabVIEW軟件，開發(fā)1 000 MW機(jī)組自然通風(fēng)冷卻塔預(yù)熱式防凍系統(tǒng)仿真計(jì)算軟件，以驗(yàn)證防凍系統(tǒng)的有效性。圖4為特性計(jì)算程序流程。

圖4 特性計(jì)算程序流程

3 模型驗(yàn)證及運(yùn)行控制

3.1 基于實(shí)測(cè)數(shù)據(jù)的模型驗(yàn)證

以寧夏石嘴山電廠1 000 MW超超臨界機(jī)組自然通風(fēng)冷卻塔[19]為例，塔型設(shè)計(jì)參數(shù)見表1。該塔淋水面積為12 000 m3，并采用雙斜波型淋水填料，高度按1.00 m、1.25 m、1.50 m由塔中心向外呈階梯狀布置，單個(gè)組裝塊尺寸為2 000 mm×500 mm×500 mm，且呈交錯(cuò)排列；噴濺裝置為TP-Ⅱ型，布置間距為1 m；收水器使用改良的聚氯乙烯(PVC)材料，片距為45 mm。

表1 自然通風(fēng)冷卻塔設(shè)計(jì)參數(shù)

根據(jù)該電廠運(yùn)行的實(shí)測(cè)數(shù)據(jù)[19]對(duì)冷卻塔數(shù)學(xué)模型進(jìn)行數(shù)值驗(yàn)證，試驗(yàn)工況分2臺(tái)泵并聯(lián)(簡(jiǎn)稱雙泵并聯(lián))運(yùn)行工況與3臺(tái)泵并聯(lián)(簡(jiǎn)稱三泵并聯(lián))運(yùn)行工況。2種工況均為全塔配水，機(jī)組滿負(fù)荷運(yùn)行。大氣干球溫度、濕球溫度采用電動(dòng)通風(fēng)干濕表測(cè)量，并且在距離冷卻塔20 m處設(shè)置2個(gè)氣象測(cè)點(diǎn)，測(cè)點(diǎn)高度為7 m；出塔水溫測(cè)點(diǎn)位于回水溝，并且使用Pt-100型鉑電阻溫度探頭；塔內(nèi)風(fēng)速測(cè)點(diǎn)共有76個(gè)，均勻布置在4條主水槽，測(cè)量?jī)x表為MSF-1型電子微風(fēng)表，試驗(yàn)結(jié)果取平均值。將測(cè)量得到的干球溫度作為環(huán)境溫度。

圖5為雙泵并聯(lián)運(yùn)行下出塔水溫的試驗(yàn)結(jié)果與仿真結(jié)果。冷卻塔模型出塔水溫的試驗(yàn)結(jié)果與仿真結(jié)果較為接近，且整體變化趨勢(shì)一致。

圖5 雙泵并聯(lián)運(yùn)行下出塔水溫的試驗(yàn)的結(jié)果與仿真的結(jié)果

出塔水溫偏差(仿真得到的出塔水溫和試驗(yàn)得到的出塔水溫的差)分布見圖6。雙泵并聯(lián)運(yùn)行下，出塔水溫偏差為0.35～1.49 K；三泵并聯(lián)運(yùn)行下，出塔水溫偏差為0.07～0.86 K。雙泵并聯(lián)運(yùn)行的總循環(huán)水量少，在相似氣象條件下溫降幅度較大，總體上溫度偏差略大于三泵并聯(lián)運(yùn)行。

圖6 出塔水溫偏差分布

表2為出塔水溫及冷卻數(shù)的誤差。考慮到冷卻塔試驗(yàn)情況下存在的測(cè)量精度、測(cè)點(diǎn)位置選擇等，可認(rèn)為冷卻塔模型計(jì)算結(jié)果合理有效，滿足工程設(shè)計(jì)與計(jì)算的精度要求。

表2 出塔水溫及冷卻數(shù)的誤差

預(yù)熱器的設(shè)計(jì)參考文獻(xiàn)[15]，經(jīng)優(yōu)化設(shè)計(jì)，翅片管總體尺寸型號(hào)為φ12×1/26/1.5/0.5-AL/AL-高頻焊接，翅化比為15.44。預(yù)熱器管束中心呈三角形叉排布置，管間距為36 mm，有60根單排管束，共有3程管程，箱體尺寸為2.16 m×0.11 m×2.00 m；初步設(shè)計(jì)冷卻塔的實(shí)際擋風(fēng)覆蓋率為65%，預(yù)熱器設(shè)置于冷卻塔進(jìn)風(fēng)口最底部(即擋風(fēng)板下方)，并沿周向連續(xù)布置，總覆蓋率為35%。

3.2 溫度控制模塊

防凍系統(tǒng)的溫度控制模塊見圖7，其主要包括溫度監(jiān)測(cè)單元、信號(hào)處理單元和調(diào)控執(zhí)行單元。具體設(shè)備分為百葉窗式調(diào)節(jié)風(fēng)門、進(jìn)水調(diào)節(jié)閥、風(fēng)門執(zhí)行機(jī)構(gòu)、調(diào)節(jié)閥執(zhí)行機(jī)構(gòu)、可編程序邏輯控制器(PLC)及各測(cè)溫元件(T)。

圖7 防凍系統(tǒng)的溫度控制模塊

測(cè)溫元件用于實(shí)時(shí)監(jiān)測(cè)進(jìn)塔風(fēng)溫、管側(cè)出口水溫及出塔水溫，PLC基于各溫度測(cè)點(diǎn)的輸入信號(hào)，通過(guò)控制邏輯判斷并輸出反饋調(diào)節(jié)信號(hào)，使系統(tǒng)執(zhí)行相應(yīng)的調(diào)控指令。溫度控制過(guò)程為閉環(huán)反饋調(diào)節(jié)過(guò)程，通過(guò)調(diào)節(jié)預(yù)熱器的進(jìn)風(fēng)量及循環(huán)水量，對(duì)預(yù)熱器及冷卻塔的運(yùn)行狀態(tài)實(shí)現(xiàn)有效監(jiān)測(cè)與控制。

3.3 運(yùn)行控制策略

防凍系統(tǒng)的溫度控制方式為基于百葉窗開度或調(diào)節(jié)閥開度的連續(xù)調(diào)節(jié)，以實(shí)現(xiàn)對(duì)預(yù)熱氣溫(冷卻塔進(jìn)風(fēng)溫度)的閉環(huán)控制。當(dāng)進(jìn)塔風(fēng)溫在-5 ℃左右，未加擋風(fēng)板的冷卻塔的結(jié)冰風(fēng)險(xiǎn)極低[11]，因此控制系統(tǒng)的防凍約束條件設(shè)定為預(yù)熱氣溫≥-5 ℃。設(shè)計(jì)工況時(shí)各參數(shù)如下：環(huán)境溫度為-15 ℃，相對(duì)濕度為40%，機(jī)組負(fù)荷為熱耗率驗(yàn)收(THA)工況，背壓為4.8 kPa。為比較不同控制方式對(duì)系統(tǒng)運(yùn)行特性的影響，在設(shè)計(jì)工況下，對(duì)不同調(diào)節(jié)方式進(jìn)行仿真模擬計(jì)算，得到預(yù)熱器及冷卻塔的熱力性能參數(shù)見圖8(基于調(diào)節(jié)閥開度)及圖9(基于百葉窗開度)，并且計(jì)算中考慮了循環(huán)水泵的功耗。

圖8 基于調(diào)節(jié)閥開度的調(diào)節(jié)特性

圖9 基于百葉窗開度的調(diào)節(jié)特性

由圖8及圖9可得：調(diào)節(jié)調(diào)節(jié)閥或百葉窗均可以有效控制預(yù)熱氣溫，并且隨著調(diào)節(jié)閥開度的增大或百葉窗開度的減小，預(yù)熱氣溫升高，降低了冷卻塔結(jié)冰的概率，但2種調(diào)節(jié)方式冷卻循環(huán)水的效果相差較多。

表3為預(yù)熱氣溫變化1 K時(shí)，不同調(diào)節(jié)方式下的經(jīng)濟(jì)性對(duì)比。當(dāng)預(yù)熱氣溫變化1 K時(shí)，調(diào)節(jié)閥開度平均變化9.6百分點(diǎn)，百葉窗開度平均變化為3.9百分點(diǎn)。與調(diào)節(jié)百葉窗開度相比，調(diào)節(jié)調(diào)節(jié)閥開度單獨(dú)控制下出塔水溫升高幅度較小，表明維持凝汽器壓力所需的循環(huán)水量增量及水泵功耗增量更小，冷端經(jīng)濟(jì)性更優(yōu)。

表3 不同調(diào)節(jié)方式下的經(jīng)濟(jì)性對(duì)比

基于上述分析，優(yōu)化控制策略為：優(yōu)先控制百葉窗全開，改變調(diào)節(jié)閥開度，維持預(yù)熱氣溫高于防凍閾值，并在防凍約束下對(duì)調(diào)節(jié)閥開度進(jìn)行尋優(yōu)，降低出塔水溫，減少水泵功耗；當(dāng)調(diào)節(jié)閥開度達(dá)到水量調(diào)節(jié)的上限，減小百葉窗開度，降低冷卻塔的進(jìn)風(fēng)量。

4 仿真結(jié)果及分析

4.1 基于氣象條件變化的最優(yōu)控制

寒冷季節(jié)時(shí)，冷卻塔內(nèi)循環(huán)水的降溫方式以接觸散熱為主，蒸發(fā)散熱強(qiáng)度較弱，空氣濕度變化對(duì)冷卻性能的影響并不大。圖10為環(huán)境溫度變化下防凍系統(tǒng)的最優(yōu)控制參數(shù)，其中仿真計(jì)算中運(yùn)行工況的負(fù)荷條件設(shè)置同設(shè)計(jì)工況。

圖10 環(huán)境溫度變化時(shí)的最優(yōu)控制參數(shù)

由圖10可得：當(dāng)環(huán)境溫度從-10 ℃降低至-28 ℃，最優(yōu)控制下的預(yù)熱氣溫始終穩(wěn)定于防凍閾值范圍，驗(yàn)證了防凍系統(tǒng)的有效性。當(dāng)環(huán)境溫度高于-22 ℃，調(diào)節(jié)方式以調(diào)節(jié)循環(huán)水量為主，控制百葉窗全開，最優(yōu)調(diào)節(jié)閥開度隨氣溫的降低而增大，出塔水溫呈下降趨勢(shì)，且存在最低溫度(8.28 ℃)。當(dāng)氣溫降至-22 ℃，水量調(diào)節(jié)已達(dá)上限，調(diào)節(jié)方式轉(zhuǎn)變?yōu)閷?duì)進(jìn)風(fēng)量的調(diào)節(jié)，減小百葉窗開度，以確保預(yù)熱氣溫仍滿足防凍要求，出塔水溫因進(jìn)風(fēng)量的降低而逐漸升高。

4.2 基于機(jī)組負(fù)荷變化的最優(yōu)控制

保持氣象條件不變(環(huán)境溫度為-15 ℃，相對(duì)濕度為40%)，改變機(jī)組的運(yùn)行負(fù)荷，得到負(fù)荷變化時(shí)的最優(yōu)控制參數(shù)見圖11。循環(huán)水量達(dá)到調(diào)節(jié)臨界值，機(jī)組負(fù)荷在(40%～50%)THA，預(yù)熱器出口水溫、出塔水溫變化的總體趨勢(shì)均隨負(fù)荷的減小而降低，但循環(huán)水量調(diào)節(jié)方式對(duì)管側(cè)水溫影響較大，出塔水溫主要受風(fēng)量調(diào)節(jié)的影響。另外，負(fù)荷率每降低10百分點(diǎn)，管側(cè)水溫(預(yù)熱器出口水溫)平均降低1.72 K，出塔水溫平均降低0.94 K。

圖11 負(fù)荷變化時(shí)的最優(yōu)控制參數(shù)

4.3 優(yōu)化控制實(shí)例

對(duì)自然通風(fēng)冷卻塔預(yù)熱式防凍系統(tǒng)按冷端運(yùn)行最優(yōu)參數(shù)進(jìn)行控制，不僅提高了進(jìn)風(fēng)塔溫度，滿足了防凍要求，而且降低了出塔水溫，提高了冷端系統(tǒng)運(yùn)行經(jīng)濟(jì)性。優(yōu)化控制的目標(biāo)函數(shù)可表示為冷端系統(tǒng)的最大凈收益，即當(dāng)減少的發(fā)電煤耗帶來(lái)的收益與循環(huán)水泵功耗增加的費(fèi)用之差達(dá)到最大值時(shí)，此時(shí)對(duì)應(yīng)的循環(huán)水量為最優(yōu)循環(huán)水量，凝汽器壓力為經(jīng)濟(jì)背壓。具體計(jì)算公式為：

maxE=max[-(L×ΔB+R×ΔP)]

(14)

式中：E為凈收益，元/h；L為標(biāo)準(zhǔn)煤價(jià)格，元/t；ΔB為標(biāo)準(zhǔn)煤質(zhì)量流量增量，t/h；R為電價(jià)，元/(kW·h)；ΔP為所有循環(huán)水泵功耗的增量，kW。

針對(duì)1 000 MW超超臨界機(jī)組經(jīng)濟(jì)性進(jìn)行分析，通過(guò)調(diào)節(jié)循環(huán)水泵轉(zhuǎn)速改變循環(huán)水量，得到不同環(huán)境溫度下優(yōu)化控制結(jié)果見表4。在不同環(huán)境溫度下，防凍系統(tǒng)均存在最優(yōu)調(diào)節(jié)閥開度，經(jīng)濟(jì)背壓低于設(shè)計(jì)背壓(4.8 kPa)，凝汽器真空度增大，提高了汽輪發(fā)電機(jī)組的實(shí)際循環(huán)熱效率。經(jīng)計(jì)算，環(huán)境溫度為-20～-10 ℃時(shí)，環(huán)境溫度每降低1 K，發(fā)電熱耗平均降低60.70 kJ/(kW·h)。同時(shí)，循環(huán)水泵功耗也因出塔水溫的優(yōu)化而降低，冷端經(jīng)濟(jì)性變好，且環(huán)境溫度每降低1 K，經(jīng)濟(jì)背壓平均降低0.09 kPa，凈收益平均增加101元/h。

表4 不同環(huán)境溫度下優(yōu)化控制結(jié)果

圖12為寧夏石嘴山2020年12月12日中各個(gè)時(shí)刻的環(huán)境溫度下經(jīng)優(yōu)化控制后的凈收益變化。

圖12 優(yōu)化控制后的凈收益變化

凈收益隨環(huán)境溫度變化的差異較大，若以標(biāo)準(zhǔn)煤價(jià)格為700元/t、廠用電價(jià)為0.4元/(kW·h)計(jì)算，經(jīng)優(yōu)化控制后發(fā)電機(jī)組一天可增加凈收益20 974元，整個(gè)冬季凈收益為188.76萬(wàn)元，具有工程應(yīng)用及推廣價(jià)值。

5 結(jié)語(yǔ)

針對(duì)自然通風(fēng)濕式冷卻塔寒冷季節(jié)結(jié)冰問(wèn)題，提出了一種預(yù)熱式防凍系統(tǒng)，通過(guò)設(shè)置氣水表面式預(yù)熱器，使冷卻塔的進(jìn)風(fēng)充分得到預(yù)熱，既提高了進(jìn)風(fēng)溫度避免結(jié)冰，又冷卻了循環(huán)水。同時(shí)，根據(jù)系統(tǒng)組成和工藝流程，建立了流動(dòng)與傳熱數(shù)學(xué)模型，并在防凍約束下建立了運(yùn)行控制模型，得到結(jié)論為：

(1) 基于LabVIEW軟件開發(fā)了可用于預(yù)熱器及冷卻塔特性計(jì)算的仿真軟件，并依據(jù)實(shí)測(cè)數(shù)據(jù)對(duì)冷卻塔模型進(jìn)行了數(shù)值驗(yàn)證，得到仿真出塔水溫變化的總體趨勢(shì)與試驗(yàn)數(shù)據(jù)一致。雙泵并聯(lián)運(yùn)行下出塔水溫偏差較大，最大平均溫度偏差為0.83 K，最大平均相對(duì)誤差為8.33%，滿足工程計(jì)算精度要求。

(2) 研究了基于調(diào)節(jié)閥或百葉窗調(diào)節(jié)的特性，得到調(diào)節(jié)調(diào)節(jié)閥開度控制預(yù)熱氣溫方式下的冷端經(jīng)濟(jì)性更優(yōu)，調(diào)節(jié)百葉窗開度對(duì)出塔水溫的影響較大。

(3) 給出了最優(yōu)控制參數(shù)，當(dāng)環(huán)境溫度或機(jī)組負(fù)荷降低時(shí)，優(yōu)先增大調(diào)節(jié)閥開度，預(yù)熱氣溫高于防凍閾值，驗(yàn)證了防凍系統(tǒng)的有效性，并且當(dāng)環(huán)境溫度為-22 ℃左右或機(jī)組負(fù)荷為50%THA左右時(shí)，循環(huán)水量達(dá)到調(diào)節(jié)臨界值，出塔水溫可維持在3～12 ℃。

(4) 對(duì)某1 000 MW機(jī)組經(jīng)濟(jì)性進(jìn)行分析和優(yōu)化，優(yōu)化控制后發(fā)電熱耗平均降低60.70 kJ/(kW·h)，整個(gè)冬季凈收益為188.76萬(wàn)元，經(jīng)濟(jì)效益明顯。