楊海生，王旻昕

（1.河北省電力研究院，石家莊市，050021；2.西柏坡發(fā)電有限責(zé)任公司，石家莊市，050400）

0 引言

給水流量作為大型機組運行中的主要測量流量，直接參與汽包水位自動調(diào)節(jié)。隨著對機組經(jīng)濟性的日異重視，許多機組的性能監(jiān)測系統(tǒng)采用了分散控制系統(tǒng)（distribution control system，DCS）中的給水流量對機組的性能進行在線計算，這對DCS中給水流量的準(zhǔn)確計算提出了很高的要求。

目前，大多數(shù)機組DCS中給水流量的計算公式為

式中：qmax、Δpmax、ρfmax分別表示流量噴嘴的刻度流量、差壓上限及設(shè)計工作參數(shù)下（對應(yīng)于刻度流量）的密度；qc、Δp、ρf分別表示實際的流量、差壓及密度。

而在密度修正計算中，通常并不直接進行密度計算，認(rèn)為液體的密度只與溫度有關(guān)。目前多數(shù)機組DCS中采用式（2）進行修正：

式中：tmax及t分別表示流量噴嘴的設(shè)計工作溫度及實際工作溫度。

實際運行表明，這種計算方法在tmax及t存在較大偏差時會引起較大的計算誤差。

如何在DCS中實現(xiàn)給水密度的相對準(zhǔn)確及簡捷的計算，對于準(zhǔn)確計算給水流量具有重要意義，本文將介紹一種簡化計算方法。

1 水密度的常用計算方法

1.1 水和水蒸氣熱力學(xué)性質(zhì)IAPW IFC97與IFC67計算公式

水和水蒸氣作為一種常規(guī)工質(zhì)，其熱力學(xué)性質(zhì)計算中主要采用IAPWS IFC計算公式。目前IFC97為最新標(biāo)準(zhǔn)，IFC97與IFC67公式在使用范圍、區(qū)域劃分、計算精度上有一定區(qū)別，但在工業(yè)蒸汽流量測量常用范圍內(nèi)，2公式計算結(jié)果偏差極小。

IFC計算公式由于涉及大量數(shù)值計算，在DCS中采用的難度較大。

1.2 查表法

這種方法需要把水蒸氣密度表數(shù)據(jù)裝入計算機中，并根據(jù)實際運行工況的溫度、壓力，從表中查出相應(yīng)的密度值。

這種方法的缺點是需要存儲大量的密度數(shù)據(jù)，也不便于在DCS中得到應(yīng)用。

2 新的水密度簡化計算公式

2.1 計算方法的基本原理

由于水物理性質(zhì)的復(fù)雜性，給水密度除與給水溫度有關(guān)外，還與給水壓力存在一定關(guān)系。壓力、溫度對水密度的影響中，溫度的作用是主要的。

對于某一類型的機組（如亞臨界機組），鍋爐給水介質(zhì)主要運行參數(shù)如壓力、溫度只能在某一特定范圍內(nèi)變化。以亞臨界機組為例，正常運行中給水壓力變化范圍為5～24 MPa，給水溫度變化范圍為100～280℃。

密度公式擬合中，可選取某一基準(zhǔn)運行工況，并將壓力、溫度對密度的影響分別進行考慮?；鶞?zhǔn)工況下的運行參數(shù)為常用機組運行參數(shù)。例如，對于亞臨界機組，基準(zhǔn)工況選取參數(shù)為壓力16 MPa，溫度240℃。

擬合的密度公式可表示為

式中：ρ及ρb分別表示實際密度與基準(zhǔn)工況的密度，kg/m3；f（p）與f（t）分別表示實際運行壓力與實際運行溫度對密度的修正系數(shù)。

2.2 溫度對密度修正系數(shù)的擬合

根據(jù)基準(zhǔn)工況的工作壓力，計算不同給水溫度下的密度變化，并轉(zhuǎn)化為修正系數(shù)，最終得到溫度與密度修正系數(shù)對應(yīng)表（或擬合函數(shù)f（t））。

首先，計算溫度與密度修正的函數(shù)?；?6 MPa設(shè)計工作壓力的數(shù)據(jù)如表1所示。

表1 基準(zhǔn)壓力、不同溫度的密度及修正系數(shù)Tab.1 Density and its correction coefficients under datum pressure and various temperatures

由表1擬合出的多項式公式可表示為

2.3 壓力對密度修正系數(shù)的擬合

根據(jù)基準(zhǔn)工況的工作溫度，計算不同給水壓力下的密度變化，并轉(zhuǎn)化為修正系數(shù)，最終得到壓力與密度修正系數(shù)對應(yīng)表（或擬合函數(shù)f（p））。

基于240℃的基準(zhǔn)工作溫度的數(shù)據(jù)如表2所示。

表2 基準(zhǔn)溫度、不同壓力的密度及修正系數(shù)Tab.2 Density and its correction coefficients under datum temperature and various pressures

由表2擬合出的多項式公式可表示為

2.4 密度擬合計算公式的誤差分析

對6～24 MPa壓力范圍，100～280 ℃溫度范圍的給水密度，采用不同方法進行了計算比較，結(jié)果如表3所示。

從計算結(jié)果來看，在12～20 MPa給水壓力范圍內(nèi)，2種計算方法的結(jié)果偏差在0.3%之內(nèi)。采用擬合公式計算給水流量，引起的誤差在0.15%以內(nèi)。

表3 不同計算方法的給水密度比較Tab.3 Comparison of feedwater densities calculated by different methods

當(dāng)給水壓力在6～8 MPa范圍內(nèi)時，雖然最大偏差達到了0.73%，但考慮到正常運行中高壓加熱器投運時的鍋爐最終給水溫度一般均高于160℃，因此實際運行中最大偏差為0.574%，由此引起給水流量計算誤差在0.28%以內(nèi)。

綜合以上分析，對于亞臨界機組，在鍋爐給水正常參數(shù)變化范圍內(nèi)，采用IFC97計算方法得到的密度與采用擬合公式得到的結(jié)果最大偏差為0.574%，可以滿足現(xiàn)場給水流量密度修正的精度要求。

3 結(jié)論

（1）提出了一種考慮實際給水壓力與給水溫度影響的密度計算方法。在密度計算中，選取某一基準(zhǔn)運行工況，并將壓力、溫度對密度的影響分別進行考慮。

（2）實例計算表明，在任意實際鍋爐給水運行參數(shù)下，采用簡化計算方法可以得到與IFC97理論計算公式偏差很小的計算結(jié)果，既能夠滿足實際工程精度需要，又便于在DCS中實施。

[1]王彥生.給水流量測量精度的提高[J].安徽電氣工程職業(yè)技術(shù)學(xué)院學(xué)報，2006（3）：74-76.

[2]鄭金龍.關(guān)于水蒸汽流量測量中的密度計算[J].計量與測試技術(shù)，2007（6）：21-21.

[3]趙慶煒.汽水流量測量的密度校正與專用校正裝置和研制[J].熱力發(fā)電，1983（1）：43-45.

[4]林萬超.火電廠熱系統(tǒng)節(jié)能理論[M].西安：西安交通大學(xué)出版社，1994.

[5]鄭體寬.熱力發(fā)電廠[M].北京：中國電力出版社，2001.

[6]沈維道.工程熱力學(xué)[M].北京：高等教育出版社，2001.

[7]趙常興.汽輪機組技術(shù)手冊[M].北京：中國電力出版社，2007.

[8]IAPWS IFC97,IAPWS Industrial Formulation 1997 for the Thermodynamic Properties of Water and Steam[S].

[9]IAPWS IFC67,IAPWS Industrial Formulation 1967 for the Thermodynamic Properties of Water and Steam[S].

[10]American Society of Mechanical Engineers.ASME Steam Properties for Industrial Use[M].New York:ASME Press,1998.