BHDB-M701F4-Q1 余熱鍋爐動(dòng)態(tài)建模與仿真

嚴(yán)夢(mèng)月 康英偉

（上海電力大學(xué) 自動(dòng)化工程學(xué)院，上海 200090）

在燃?xì)?蒸汽聯(lián)合循環(huán)機(jī)組中，余熱鍋爐對(duì)煙氣能量的利用嚴(yán)重影響著整個(gè)系統(tǒng)的效率。目前，國(guó)內(nèi)外諸多學(xué)者針對(duì)余熱鍋爐數(shù)學(xué)模型的建立及其動(dòng)態(tài)特性的仿真分析進(jìn)行大量研究[1-4]。丁艷軍等采用序貫?zāi)K法，建立了各單元之間連接關(guān)系的系統(tǒng)結(jié)構(gòu)模型[5]。高雷等以面向?qū)ο蟮乃枷?，分析并歸納了以傳遞函數(shù)為研究對(duì)象的序貫?zāi)K法的建模過(guò)程[6]。陳維春等采用模塊化建模方法，建立了雙壓余熱鍋爐的動(dòng)態(tài)數(shù)學(xué)模型，并通過(guò)數(shù)值計(jì)算分析了其動(dòng)態(tài)特性[7]。馬博等介紹了余熱鍋爐的熱力特性以及系統(tǒng)劃分、建模方法、參數(shù)選取、二次建模等，采用集總參數(shù)法詳細(xì)闡述設(shè)計(jì)、動(dòng)態(tài)建模及計(jì)算方面的問(wèn)題[8]。

文章針對(duì)燃?xì)?蒸汽聯(lián)合循環(huán)機(jī)組中使用的BHDBM701F4-Q1 型號(hào)的三壓再熱自然循環(huán)余熱鍋爐，采取機(jī)理建模和動(dòng)態(tài)仿真的方法分析余熱鍋爐的動(dòng)態(tài)特性。將余熱鍋爐整體系統(tǒng)劃分成單相介質(zhì)換熱器、汽包、下降管和上升管等模塊，并分別構(gòu)建相應(yīng)的動(dòng)態(tài)數(shù)學(xué)模型。使用MATLAB 軟件，建立余熱鍋爐系統(tǒng)的仿真模型，根據(jù)仿真結(jié)果研究其動(dòng)態(tài)特性。

1 余熱鍋爐的構(gòu)成及原理

余熱鍋爐主要由省煤器、汽包、蒸發(fā)器、過(guò)熱器等換熱器組成。文章研究的余熱鍋爐采用臥式布置，水平方向流動(dòng)的煙氣與垂直布置的受熱面進(jìn)行熱量交換，產(chǎn)生過(guò)熱蒸汽。

凝結(jié)水泵將給水送至低壓省煤器，水在省煤器內(nèi)預(yù)熱至近飽和溫度時(shí)進(jìn)入低壓汽包。低壓汽包中的一部分水經(jīng)過(guò)低壓蒸發(fā)器和過(guò)熱器，產(chǎn)生低壓過(guò)熱蒸汽進(jìn)入低壓汽缸；另一部分水作為中壓汽包與高壓汽包的給水，各自流至其給水系統(tǒng)。中壓給水通過(guò)增壓泵流入中壓省煤器，再進(jìn)入中壓汽包，流入中壓蒸發(fā)器和過(guò)熱器，生成中壓過(guò)熱蒸汽，與冷段再熱蒸汽匯合產(chǎn)生再熱蒸汽，送往中壓汽缸。高壓系統(tǒng)的水經(jīng)高壓給水泵增壓后流經(jīng)高壓省煤器，再進(jìn)入高壓汽包，流入高壓蒸發(fā)器與高壓過(guò)熱器，產(chǎn)生高壓過(guò)熱蒸汽，送進(jìn)高壓汽缸。

2 余熱鍋爐各模塊數(shù)學(xué)模型

根據(jù)管內(nèi)工質(zhì)狀態(tài)將余熱鍋爐分為單相換熱器和雙相換熱器。余熱鍋爐數(shù)學(xué)模型主要由單相介質(zhì)換熱模型和蒸發(fā)系統(tǒng)模型構(gòu)成。省煤器、過(guò)熱器和再熱器為單相換熱器，只需構(gòu)建一個(gè)基本模塊。蒸發(fā)系統(tǒng)由汽包、下降管和上升管構(gòu)成，主要完成相變換熱的過(guò)程，為雙相換熱器。

2.1 省煤器和過(guò)熱器模型

建立數(shù)學(xué)模型時(shí)，把單相受熱面的模型簡(jiǎn)化為一根無(wú)限長(zhǎng)的單管，并進(jìn)行如下假設(shè)：第一，用一根受熱管等效所有并聯(lián)管，其長(zhǎng)度與單管相同，通流面積為各管通流面積之和；第二，將煙氣與管壁、管壁與管內(nèi)工質(zhì)之間視為只有徑向換熱，且換熱強(qiáng)度一致；第三，管壁徑向無(wú)溫差；第四，沿管長(zhǎng)方向無(wú)導(dǎo)熱；第五，煙氣進(jìn)出口流量相同；第六，煙氣為理想氣體。

工質(zhì)側(cè)質(zhì)量平衡方程為

式中：Dw,in、Dw,out分別為工質(zhì)入口流量和出口流量；ρw,out為工質(zhì)出口密度；Vw為工質(zhì)體積；t為時(shí)間。

工質(zhì)側(cè)能量平衡方程為

式中：Q2為金屬對(duì)工質(zhì)的放熱量；hw,in、hw,out分別為進(jìn)口工質(zhì)焓值和出口工質(zhì)焓值。

動(dòng)量方程為

式中：Pw,in、Pw,out分別為入口工質(zhì)壓力和出口工質(zhì)壓力；Kd為壓損系數(shù)；ρw,in為工質(zhì)進(jìn)口密度。

管壁的吸熱方程為

式中：Q1為煙氣向管壁金屬的傳熱量；Mj為金屬質(zhì)量；Cj為金屬比熱容；Tj為金屬溫度。

金屬向介質(zhì)放熱的放熱方程為

式中：A為換熱面積；K2為換熱系數(shù)；Tw,in、Tw,out分別為工質(zhì)入口溫度和出口溫度。

煙氣向金屬放熱的放熱方程為

式中：K1為換熱系數(shù)；Dg為煙氣流量；Tg,in、Tg,out分別為進(jìn)口煙氣溫度和出口煙氣溫度。

煙氣側(cè)能量平衡方程為

式中：Vg為煙氣體積；ρg為煙氣密度；Cg為煙氣比熱容。

2.2 余熱鍋爐蒸發(fā)系統(tǒng)模型

為建立更精確的動(dòng)態(tài)數(shù)學(xué)模型，分別建立汽包、下降管和上升管這3 個(gè)部分的模型。3 個(gè)部分形成循環(huán)回路，因管內(nèi)工質(zhì)的密度不同而產(chǎn)生的壓力差推動(dòng)工質(zhì)循環(huán)流動(dòng)[9]。余熱鍋爐自然循環(huán)的蒸發(fā)系統(tǒng)如圖1 所示。

建立蒸發(fā)系統(tǒng)的數(shù)學(xué)模型時(shí)，對(duì)蒸發(fā)系統(tǒng)進(jìn)行如下假設(shè)：第一，全部給水直接進(jìn)入下降管；第二，汽包內(nèi)工質(zhì)壓力均勻分布；第三，蒸發(fā)區(qū)內(nèi)水的密度按飽和水密度計(jì)算；第四，汽水完全分離下降管不帶汽；第五，用出口參數(shù)代表每個(gè)環(huán)節(jié)集總參數(shù)；第六，煙氣進(jìn)出口流量相同；第七，煙氣為理想氣體。

2.2.1 汽包模型

在汽水系統(tǒng)中，汽包是汽和水的銜接部分，是產(chǎn)汽過(guò)程中不可或缺的裝置，因此對(duì)其模型的要求很高。

汽包水側(cè)質(zhì)量平衡方程為

式中：Dsm,out省煤器出口工質(zhì)流量；Dss,out為上升管出口工質(zhì)流量；Xqs為含汽率；Ddz、Dlp分別為汽包內(nèi)動(dòng)態(tài)蒸發(fā)量和連續(xù)排污量；Dxj,in為下降管入口流量；Mw為汽包水質(zhì)量。

汽包內(nèi)動(dòng)態(tài)蒸發(fā)量的計(jì)算公式為[10]

式中：kzf為動(dòng)態(tài)蒸發(fā)系數(shù)；hw為汽包水焓；hbs為飽和水焓。

汽包水位L的計(jì)算公式為

式中：vw為水比容；Vqb為汽包總?cè)莘e；A為水位截面積。

水側(cè)能量平衡方程為

式中：hsm,out為省煤器出口水焓；Qb為汽包水對(duì)底部的傳熱量。

汽側(cè)質(zhì)量平衡方程為

式中：Ds為主蒸汽流量；Dfq為汽包放汽量；Mv為蒸汽質(zhì)量。

主蒸汽流量和汽包放汽量的計(jì)算公式分別為

式中：Bs、Bfq為線性化導(dǎo)納矩陣[11]；Pqb為汽包壓強(qiáng)；Ps為主蒸汽壓強(qiáng)；Pa為大氣壓強(qiáng)。

將式（13）和式（14）代入式（12），可得汽包蒸汽質(zhì)量微分方程，表達(dá)式為

假設(shè)汽包內(nèi)蒸汽容積Vv為常數(shù)，可得

式中：ρv為蒸汽密度。

2.2.2 下降管模型

能量平衡方程為

式中：Vxj為下降管容積；ρxj,out、hxj,out、Dxj,out分別為下降管工質(zhì)出口水密度、焓和流量；Qxj為下降管工質(zhì)的放熱量。

工質(zhì)對(duì)管壁的放熱量計(jì)算公式為

式中：kxj為下降管中工質(zhì)對(duì)管壁的放熱系數(shù)；Txj,in、Txj,out分別為下降管工質(zhì)入口溫度和出口溫度。

下降管金屬能量平衡表達(dá)式為

式中：Qa為下降管管壁對(duì)大氣的傳熱量，計(jì)算公式為

式中：ka為管壁對(duì)大氣的傳熱系數(shù)；Ta為大氣溫度。

下降管出口工質(zhì)流量計(jì)算公式為

下降管工質(zhì)出口壓強(qiáng)Pxj,out的計(jì)算公式為

式中：Pqb為汽包壓強(qiáng)；lxj為下降管的長(zhǎng)度；ζxj為下降管的阻力系數(shù)。

2.2.3 上升管模型

上升管通過(guò)吸收煙氣熱量加熱下降管中的水，將其中一部分水變?yōu)轱柡驼羝?。上升管模型除上升管外，還包括引入管、引出管和汽水分離器。

工質(zhì)能量平衡方程為

式中：ρss,out為上升管出口工質(zhì)密度；Vss為上升管容積；hss,out為上升管出口工質(zhì)焓值；Q2為上升管壁對(duì)工質(zhì)的傳熱量。

管壁吸熱方程為

煙氣向金屬的對(duì)流放熱方程為

金屬對(duì)介質(zhì)的放熱方程為

式中：Tss,in、Tss,out分別為上升管入口工質(zhì)溫度和出口工質(zhì)溫度。

煙氣向金屬的輻射放熱方程為

上升管出口工質(zhì)流量的計(jì)算公式為

式中：Dpz為膨脹流量；Dxl為泄漏流量。

熱水段高度lrs的計(jì)算公式為

式中：C1為修正系數(shù)；lssg為上升管總高度。

上升管系統(tǒng)總壓頭ΔPss的計(jì)算公式為[12]

式中：ΔPyr為引入管壓頭；ΔPrs為熱水段壓頭；ΔP2為汽液雙相和汽水引出管壓頭；ΔPyc為分離器壓頭。

3 余熱鍋爐仿真模型建立

利用MATLAB 仿真平臺(tái)，基于各模塊的動(dòng)態(tài)數(shù)學(xué)模型分別建立算法模塊，根據(jù)各模塊間的關(guān)系將其連接起來(lái)得到仿真模型，如圖2 所示。各壓力級(jí)的仿真模塊相同，只需將建立的低、中、高3 個(gè)壓力級(jí)的模型組合成整體的仿真模型即可。

圖2 余熱鍋爐單個(gè)壓力級(jí)仿真模型連接圖

4 余熱鍋爐的動(dòng)態(tài)仿真分析

4.1 給水流量擾動(dòng)下的動(dòng)態(tài)仿真分析

煙氣流量保持不變，對(duì)給水流量加入階躍擾動(dòng)，流量降低10%，即省煤器入口流量從116.860 kg·s-1階躍下降到105.14 kg·s-1。在該擾動(dòng)下余熱鍋爐高壓模型的動(dòng)態(tài)響應(yīng)曲線如圖3 所示。

圖3 給水流量擾動(dòng)下的高壓模型動(dòng)態(tài)響應(yīng)曲線

圖3（a）為汽包壓強(qiáng)的動(dòng)態(tài)響應(yīng)曲線。當(dāng)給水流量減少時(shí)，隨著時(shí)間的增加汽包壓強(qiáng)先變大然后逐漸趨向穩(wěn)定。流入汽包的水的過(guò)冷度會(huì)隨著流量的減少而下降，使飽和蒸汽產(chǎn)量增加，壓強(qiáng)升高，但高壓蒸發(fā)器入口處的煙溫沒(méi)有大幅度變化，因此汽包壓強(qiáng)增加幅度較小。圖3（b）為高壓過(guò)熱蒸汽溫度變化曲線。當(dāng)給水流量降低時(shí)，過(guò)熱蒸汽溫度隨著時(shí)間的增加而下降，最后逐漸穩(wěn)定。單位質(zhì)量的蒸汽吸熱量隨著給水流量的減少而增加，吸熱量增加使蒸汽流量增多。受蒸汽流量的影響，過(guò)熱蒸汽溫度先下降，然后趨于穩(wěn)定。圖3（c）為汽包水位變化曲線。當(dāng)給水流量減少時(shí)，高壓汽包的水位隨著時(shí)間的增加而不斷降低，主要原因是給水流量不足導(dǎo)致汽包壓強(qiáng)增大，使得汽包蒸發(fā)量增大，當(dāng)給水量少于汽包蒸發(fā)量時(shí)，汽包水位下降。圖3（d）為過(guò)熱器出口蒸汽流量變化曲線，高壓過(guò)熱蒸汽流量隨著時(shí)間的增加而增大。汽包的水欠焓隨著給水流量的減少而變小，部分飽和水變成飽和汽，使得汽包壓強(qiáng)瞬間增大，蒸汽流量隨之上升，而后逐漸趨于穩(wěn)定。

4.2 煙氣流量擾動(dòng)下的動(dòng)態(tài)仿真分析

給水流量保持不變，對(duì)煙氣流量加入階躍擾動(dòng)，流量降低10%，即入口煙氣流量從623.0 kg·s-1階躍下降到560.7 kg·s-1。該擾動(dòng)下余熱鍋爐高壓模型的動(dòng)態(tài)響應(yīng)曲線如圖4 所示。

圖4 煙氣流量擾動(dòng)下的高壓模型動(dòng)態(tài)響應(yīng)曲線

圖4（a）為汽包壓強(qiáng)的動(dòng)態(tài)響應(yīng)曲線。當(dāng)煙氣流量降低時(shí)，汽包壓強(qiáng)隨著時(shí)間的增加先減小然后逐漸趨向穩(wěn)定。煙氣與蒸發(fā)器間的換熱量降低，使汽包產(chǎn)汽量減少，當(dāng)汽包產(chǎn)汽量小于出口蒸汽量時(shí)，汽包壓強(qiáng)降低。隨著壓強(qiáng)的下降，汽包出口蒸汽量降低，當(dāng)汽包產(chǎn)汽量與汽包出口蒸汽量相近時(shí)，壓強(qiáng)再次趨于穩(wěn)定。圖4（b）為過(guò)熱蒸汽出口溫度的變化曲線。由于吸熱量降低，蒸汽溫度急劇下降。在蒸汽溫度下降的同時(shí)，汽包蒸發(fā)量也在降低，導(dǎo)致過(guò)熱器入口處的蒸汽量下降。而蒸汽流量下降的影響遠(yuǎn)大于吸熱量下降的影響，使得蒸汽溫度緩慢上升，最后達(dá)到新的平衡點(diǎn)。圖4（c）為汽包水位的動(dòng)態(tài)響應(yīng)曲線。當(dāng)排煙流量減少時(shí)，汽包壓強(qiáng)降低，汽包蒸發(fā)量降低，當(dāng)給水量大于蒸發(fā)量時(shí)，汽包水位上升。因給水量保持不變，汽包水位持續(xù)上升。圖4（d）為過(guò)熱器出口蒸汽流量的變化曲線。蒸汽流量隨著時(shí)間的增加先降低然后逐漸趨于穩(wěn)定，其變化趨勢(shì)與汽包壓強(qiáng)的變化趨勢(shì)基本一致。換熱面的吸熱量會(huì)隨著煙氣流量的減少而降低，導(dǎo)致蒸汽流量下降。

5 結(jié)語(yǔ)

文章給出BHDB-M701F4-Q1 型號(hào)的三壓再熱自然循環(huán)余熱鍋爐的模型，描述了在給水流量擾動(dòng)與煙氣流量擾動(dòng)下不同參數(shù)的動(dòng)態(tài)過(guò)程，并對(duì)其動(dòng)態(tài)特性進(jìn)行分析，為余熱鍋爐運(yùn)行特性和控制系統(tǒng)的研究提供了可靠的理論依據(jù)和數(shù)學(xué)模型。仿真結(jié)果表明，所建立的模型能夠正確反映對(duì)象的動(dòng)態(tài)特性，尤其是擾動(dòng)下各參數(shù)的動(dòng)態(tài)變化，具有一定的工程實(shí)用價(jià)值，可以用于該余熱鍋爐的控制系統(tǒng)的設(shè)計(jì)、優(yōu)化等。