壓水堆核電機(jī)組火用成本分析的矩陣算法

李永華，蒲亮，周思遠(yuǎn)

(1. 電站設(shè)備狀態(tài)檢測(cè)與控制教育部重點(diǎn)實(shí)驗(yàn)室(華北電力大學(xué))，河北省保定市 071003；2.國核工程有限公司，上海市 200233)

壓水堆核電機(jī)組火用成本分析的矩陣算法

李永華1，蒲亮1，周思遠(yuǎn)2

(1. 電站設(shè)備狀態(tài)檢測(cè)與控制教育部重點(diǎn)實(shí)驗(yàn)室(華北電力大學(xué))，河北省保定市 071003；2.國核工程有限公司，上海市 200233)

為進(jìn)行復(fù)雜能量系統(tǒng)的熱經(jīng)濟(jì)學(xué)分析和評(píng)價(jià)，找出提高核電機(jī)組熱經(jīng)濟(jì)性的新途徑，基于火用經(jīng)濟(jì)學(xué)理論和矩陣算法，建立了壓水堆核電機(jī)組熱力系統(tǒng)火用成本分析通用模型。以某900 MW壓水堆核電機(jī)組熱力系統(tǒng)為例，進(jìn)行了39股火用流的單位火用成本計(jì)算及分析，并對(duì)具體設(shè)備提出了明確的優(yōu)化改進(jìn)措施。結(jié)果表明:沿?zé)崃ρh(huán)方向，各股火用流的單位火用成本先增加后減?。荒鬏敵龌鹩昧鞯膯挝换鹩贸杀咀罡?，為7.1754，其次是凝結(jié)水泵和1號(hào)低壓加熱器，分別為5.249和4.1911。矩陣算法較常規(guī)火用成本分析法具有構(gòu)造簡單、矩陣元素填寫法則簡便、物理意義明確和規(guī)律性強(qiáng)等優(yōu)點(diǎn)；利用矩陣算法的易編程性，便于開發(fā)出相應(yīng)的計(jì)算及分析軟件，進(jìn)而為核電機(jī)組熱力系統(tǒng)優(yōu)化及故障診斷奠定理論基礎(chǔ)。

核電機(jī)組；壓水堆；熱力系統(tǒng)；單位火用成本；矩陣方程

0 引 言

核電作為緩解能源緊張和減排溫室氣體的重要手段，具有廣泛的發(fā)展前景[1]。隨著我國核電機(jī)組裝機(jī)容量的增加，核燃料(不可再生資源)消耗巨大，且乏燃料處理難度大、成本高，對(duì)環(huán)境危害大。在保證核安全的基礎(chǔ)下，提高核電機(jī)組的熱經(jīng)濟(jì)性對(duì)節(jié)約核燃料和降低發(fā)電成本具有重要的現(xiàn)實(shí)意義。

先進(jìn)的能量系統(tǒng)分析方法為評(píng)價(jià)熱功轉(zhuǎn)換過程的完善性及提高核電廠的熱經(jīng)濟(jì)性提供新的途徑。Dunbar[2-4]基于熱力學(xué)第二定律的火用分析法對(duì)運(yùn)行中的核電廠進(jìn)行了詳細(xì)的火用分析，得出火用損失最大設(shè)備為反應(yīng)堆，其次為汽輪機(jī)，并給出了具體的改進(jìn)措施。Sayyaadi[5]應(yīng)用熱力學(xué)和火用經(jīng)濟(jì)學(xué)理論對(duì)典型的1 000 MW級(jí)壓水堆核電機(jī)組進(jìn)行了優(yōu)化分析研究，對(duì)比分析3個(gè)不同層次的優(yōu)化方法，最終得出火用經(jīng)濟(jì)學(xué)多目標(biāo)優(yōu)化為最佳優(yōu)化方案?；鹩媒?jīng)濟(jì)學(xué)理論的發(fā)展，使傳統(tǒng)熱力學(xué)分析得到進(jìn)一步完善，Valero[6-10]等多位學(xué)者已將此理論應(yīng)用于火電機(jī)組的成本分析、系統(tǒng)優(yōu)化及故障診斷，但國內(nèi)對(duì)核電機(jī)組進(jìn)行火用成本分析的研究成果尚不多見。

我國核電發(fā)展路線以壓水堆技術(shù)為主，壓水堆核電機(jī)組與火電機(jī)組熱力系統(tǒng)有著明顯的不同，主要表現(xiàn)為核電站以核反應(yīng)堆替代火電站的鍋爐提供熱能，汽輪機(jī)主要工作在濕蒸汽區(qū)，高壓缸排汽進(jìn)入再熱器之前需要進(jìn)行汽水分離等。因此，火電機(jī)組的火用成本分析方法很難應(yīng)用于現(xiàn)行的核電機(jī)組，且常規(guī)方法求取能流的單位火用成本時(shí)需要對(duì)每個(gè)研究對(duì)象列出火用成本平衡方程，不易編制通用計(jì)算程序，變工況或在線計(jì)算工作量大。

為此，本文針對(duì)壓水堆核電機(jī)組熱力系統(tǒng)的特點(diǎn)，基于火用經(jīng)濟(jì)學(xué)理論和矩陣算法，研究適用于壓水堆核電機(jī)組的火用成本分析模型，推導(dǎo)熱力系統(tǒng)火用成本分布的通用矩陣方程，通過實(shí)例計(jì)算驗(yàn)證矩陣方程的正確性，分析熱力系統(tǒng)單位火用成本的分布規(guī)律，并對(duì)具體設(shè)備提出明確的優(yōu)化改進(jìn)措施。

1 火用成本分析模型

1.1 物理模型

根據(jù)常規(guī)火電機(jī)組熱力系統(tǒng)劃分原則，結(jié)合壓水堆核電機(jī)組熱力系統(tǒng)的特點(diǎn)，將其劃分為6個(gè)局部熱力子系統(tǒng)：回?zé)嶙酉到y(tǒng)、原子鍋爐子系統(tǒng)、汽輪發(fā)電機(jī)子系統(tǒng)、汽水分離再熱器子系統(tǒng)、給水泵子系統(tǒng)、凝汽器及凝結(jié)水泵子系統(tǒng)，如圖1所示。本文中各種輔助汽水的分類可參考文獻(xiàn)[11]，主要參數(shù)含義如下：e為比火用 ；D為流量；W為功率；k*為單位火用成本，下標(biāo)表示流的屬性和位置。

圖1 壓水堆核電機(jī)組熱力系統(tǒng)

1.2 數(shù)學(xué)模型

在劃分好子系統(tǒng)的基礎(chǔ)上，根據(jù)質(zhì)量平衡方程和火用成本平衡方程，推導(dǎo)出局部熱力子系統(tǒng)的火用成本分析方程。原子鍋爐子系統(tǒng)、汽水分離再熱器子系統(tǒng)的方程具體推導(dǎo)過程如下。

(1)原子鍋爐子系統(tǒng)火用成本分析的通用物理模型見圖2。根據(jù)質(zhì)量平衡和火用成本平衡，建立原子鍋爐子系統(tǒng)火用成本平衡方程，見下式。

(1)

(2)

(3)

(4)

圖2 原子鍋爐子系統(tǒng)火用成本分析通用物理模型

(2)汽水分離再熱器子系統(tǒng)火用成本分析的通用物理模型見圖3。根據(jù)質(zhì)量平衡和火用成本平衡，建立汽水分離再熱器子系統(tǒng)火用成本平衡方程，見下式。

圖3 汽水分離再熱器子系統(tǒng)火用成本分析通用物理模型

(5)

(6)

(7)

(8)

回?zé)嶙酉到y(tǒng)、汽輪發(fā)電機(jī)子系統(tǒng)、給水泵子系統(tǒng)和凝結(jié)水及凝結(jié)水泵子系統(tǒng)與常規(guī)火電機(jī)組熱力子系統(tǒng)類似，故方程推導(dǎo)及符號(hào)說明詳見參考文獻(xiàn)[12]，這里不再贅述。

將各子系統(tǒng)的火用成本分析方程歸并到一個(gè)統(tǒng)一的矩陣中去，得到整個(gè)熱力系統(tǒng)的火用成本分析矩陣方程

(9)

式中：

D1、D2……Dn為各級(jí)回?zé)峒訜崞鞒槠髁?；DH、DM、DL為汽輪機(jī)高、中、低壓缸的進(jìn)汽流量；Dc為低壓缸排入凝汽器的蒸汽流量。

∑Diki為進(jìn)出第i級(jí)加熱器的第一類輔助蒸汽流量和，i取1，2…n；∑DtlH、∑DtlM、∑DtlL為進(jìn)出汽輪機(jī)高、中、低壓缸的輔助蒸汽流量和，如軸封漏氣等；∑Dc,kj為進(jìn)入凝汽器的第一類輔助汽水流量和。

Dw1、Dw2…Dwn為各級(jí)回?zé)峒訜崞鹘o水(凝結(jié)水)流量；Dgs為給水泵中的給水流量。

2 算例分析

2.1 機(jī)組簡介

本文以某900 MW級(jí)壓水堆核電機(jī)組為研究對(duì)象，圖4為其原則性熱力系統(tǒng)圖，該機(jī)組由1個(gè)雙流道高壓缸和3個(gè)雙流道、雙排汽低壓缸組成，回?zé)嵯到y(tǒng)為“兩高、四低、一除氧”。反應(yīng)堆額定熱功率為2 895 MW，汽輪發(fā)電機(jī)組額定輸出功率為983.8 MW，其他參數(shù)見文獻(xiàn)[13]。

2.2 補(bǔ)充方程的建立原則

(1)從外部輸入系統(tǒng)火用流的單位火用成本均規(guī)定為1。(2)子系統(tǒng)的非產(chǎn)品性輸出火用的單位火用成本與輸入火用的單位火用成本相等，所有產(chǎn)品火用流的單位火用成本相等。(3)若子系統(tǒng)的“燃料”為雙線流，則構(gòu)成雙線流的2股火用流的單位火用成本相等[14]。

2.3 單位火用成本計(jì)算及分析

在熱經(jīng)濟(jì)學(xué)分析中，將系統(tǒng)中的火用流賦予能量意義上的價(jià)值——火用成本，即獲得該股火用流所需要消耗的外部火用，生產(chǎn)單位火用流所需要消耗的外部火用稱為單位火用成本。該熱力系統(tǒng)有39股火用流(如圖4上火用流編號(hào)為1～39)，這就需要計(jì)算39個(gè)單位火用成本。根據(jù)式(1)～(4)可以求解15個(gè)不同的單位火用成本，因此矩陣方程需要補(bǔ)充24個(gè)單位火用成本相等的關(guān)系式進(jìn)行簡化。由補(bǔ)充方程的建立原則可得，外部輸入系統(tǒng)的核能火用的單位火用成本為1；第2、3、4、5、6、7、8、29、30、35、36、37股火用流的單位火用成本相等；第9、10、11、12、13、14、15、31、32、33、34股火用流的單位火用成本相等；第18、38、39股火用流的單位火用成本相等。綜上所述，可得24個(gè)單位火用成本相等的關(guān)系式使矩陣方程封閉，有唯一解。通過Matlab軟件編程計(jì)算，結(jié)果見表1。

圖4 某900 MW壓水堆核電機(jī)組原則性熱力系統(tǒng)火用流圖

表1 某900 MW壓水堆核電機(jī)組單位火用成本計(jì)算結(jié)果

由表1可看出：

(1)本文所述方法與傳統(tǒng)算法(符號(hào)模式熱經(jīng)濟(jì)學(xué))計(jì)算的結(jié)果完全一致，驗(yàn)證了矩陣方程的正確性和準(zhǔn)確性。

(2)沿?zé)崃ρh(huán)方向，各股火用流的單位火用成本先增加后減小?？拷鼰崃ο到y(tǒng)前端的單位火用成本較低，如主蒸汽(火用流編號(hào)為2)和主給水(火用流編號(hào)為28)的單位火用成本分別為2.2761、2.5776；遠(yuǎn)離熱力系統(tǒng)前端的單位火用成本較大，如凝汽器熱阱出口凝結(jié)水(火用流編號(hào)為19)的單位火用成本為7.1754。給水回?zé)嵯到y(tǒng)中，單位火用成本從低加到高加逐級(jí)減小，其中1號(hào)低加出口凝結(jié)水(火用流編號(hào)為21)的單位火用成本最大，為4.1911；7號(hào)高加出口給水(火用流編號(hào)為28)的單位火用成本最小，為2.5776。

(3)輸出火用流的單位火用成本較高的設(shè)備有凝汽器、凝結(jié)水泵和1號(hào)低壓加熱器，相比其他設(shè)備，若減少相同的火用損，其節(jié)能潛力更大。因此，應(yīng)設(shè)法減少凝汽器和1號(hào)低加內(nèi)不可逆性引起的火用損失，如采取及時(shí)清理凝汽器和1號(hào)低加內(nèi)的污垢，排除不凝結(jié)氣體，確保1號(hào)低加疏水暢通等措施[15]。

(4)發(fā)電機(jī)輸出電能(火用流編號(hào)為18)的單位火用成本為2.826 3，若引入核燃料的價(jià)格，可計(jì)算出發(fā)電機(jī)輸出電能火用的熱經(jīng)濟(jì)學(xué)成本，再考慮設(shè)備折舊、維修及勞務(wù)工資等非能量費(fèi)用，可得出發(fā)電的實(shí)時(shí)成本，進(jìn)而為上網(wǎng)電價(jià)的制定提供參考。

3 結(jié) 論

(1)基于火用經(jīng)濟(jì)學(xué)理論和矩陣算法，推導(dǎo)出壓水堆核電機(jī)組熱力系統(tǒng)的火用成本分布矩陣方程。該矩陣方程填寫簡單，各項(xiàng)物理意義明確，通用性強(qiáng)，且矩陣算法的易編程性使得單位火用成本的計(jì)算更加清晰和簡便。

(2)應(yīng)用該方程將全面揭示核電機(jī)組熱力系統(tǒng)中各股能流的單位火用成本的分布規(guī)律，進(jìn)而得出局部火用損對(duì)系統(tǒng)整體能量消耗的影響，為節(jié)能潛力挖掘指明了方向。

(3)在計(jì)算出單位火用成本的基礎(chǔ)上，引入核燃料的價(jià)格及非能量費(fèi)用，可對(duì)核電機(jī)組熱力系統(tǒng)進(jìn)行進(jìn)一步的熱經(jīng)濟(jì)學(xué)分析。從而為確定合理的發(fā)電成本提供依據(jù)，為核電機(jī)組熱力系統(tǒng)優(yōu)化及故障診斷奠定理論基礎(chǔ)。

[1]祁恩蘭.我國核電的建設(shè)形勢(shì)及思考[J].電力建設(shè)，2009，30(5)：1-4.

[2]Dunbar W R, Moody S D, Lior N. Energy analysis of an operating boiling-water-reactor nuclear power station [J]. Energy conversion and management, 1995, 36(3): 149-159.

[3]Durmayaz A, Yavuz H. Energy analysis of a pressurized-water reactor nuclear-power plant [J]. Applied Energy, 2001, 69(1): 39-57.

[4]崔雨.基于火用分析的900 MW級(jí)PWR核電機(jī)組熱經(jīng)濟(jì)性研究[D].北京：華北電力大學(xué)，2011.

[5]Sayyaadi H, Sabzaligol T. Energoeconomic optimization of a 1000 MW light water reactor power generation system[J]. International Journal of Energy Research,2009,33(4):378-395.

[6]Valero A, Lozano M A, Bartolomé J L. On-line monitoring of power-plant performance, using exergetic cost techniques[J]. Applied thermal engineering, 1996, 16(12): 933-948.

[7]張超.復(fù)雜能量系統(tǒng)的熱經(jīng)濟(jì)學(xué)分析與優(yōu)化[D].武漢：華中科技大學(xué)，2006.

[8]程偉良，王清照，王加璇.300 MW 凝汽機(jī)組的熱經(jīng)濟(jì)學(xué)成本診斷[J].中國電機(jī)工程學(xué)報(bào)，2005，25(8)：126-129.

[9]Zhang C, Wang Y, Zheng C, et al. Energy cost analysis of a coal fired power plant based on structural theory of thermoeconomics[J]. Energy Conversion and Management, 2006,47(7):817-843.

[10]肖衛(wèi)杰.結(jié)構(gòu)理論研究及在熱力系統(tǒng)故障診斷中的應(yīng)用[D].北京：華北電力大學(xué)，2002.

[11]閆順林，張春發(fā)，李永華，等.火電機(jī)組熱力系統(tǒng)汽水分布通用矩陣方程[J].中國電機(jī)工程學(xué)報(bào)，2000，20(8)：69-73.

[12]李永華，劉偉庭.火電機(jī)組熱力系統(tǒng)火用成本分布通用矩陣方程[J].動(dòng)力工程學(xué)報(bào)，2012，32(3)：243-248.

[13]廣東核電培訓(xùn)中心.900 MW 壓水堆核電站系統(tǒng)與設(shè)備[M].北京：原子能出版社，2005：319-372.

[14]王加璇，張恒良.動(dòng)力工程熱經(jīng)濟(jì)學(xué)[M].北京：水利水電出版社，1995：82-83.

[15]廖軍林,黨小建.220 MW機(jī)組低壓加熱器疏水不暢分析及改造[J].電力建設(shè)，2012，33(11)：46-48.

(編輯：蔣毅恒)

MatrixAlgorithmofExergyCostAnalysisonPWRNuclearPowerUnits

LI Yonghua1, PU Liang1, ZHOU Siyuan2

(1. Key Laboratory of Ministry of Education of Condition Monitoring and Control for Power Plant Equipment, North China Electric Power University, Baoding 071003, Hebei Province, China;2. State Nuclear Power Engineering Company, Shanghai 200233, China)

In order to implement thermoeconomics analysis and evaluation for complex energy system and find new ways to improve the thermal economy of nuclear power unit, a general model of exergy cost analysis on thermal system of pressurized water reactor (PWR)nuclear power units was established based on exergoeconomics theory and matrix algorithm. Taking the thermal system of a 900 MW PWR nuclear power unit as an example, the unit exergy cost of 39 exergy flows were calculated and analyzed, and explicit optimization improvements were proposed for the specific equipment. The results show that the unit exergy cost of each exergy flow increases firstly and then decreases along the direction of thermodynamic cycle; the maximum unit exergy cost is that of condenser output exergy flow, which is 7.1754, and followed by that of condensate pump and 1 low pressure heater, 5.249 and 4.1911 respectively. Compared with the traditional exergy cost analysis, the matrix algorithm has advantages in simple structure, convenient filling for each element of matrix, explicit physical meanings and obvious regularity. According to the programmability of matrix algorithm, a corresponding calculation and analysis software can be conveniently developed, which can provide a theoretical basis for the thermal system optimization and fault diagnosis of nuclear power units.

nuclear power units; pressurized water reactor; thermal system; unit exergy cost; matrix equation

TK 123

： A

： 1000-7229(2014)06-0137-05

10.3969/j.issn.1000-7229.2014.06.026

2013-12-16

：2014-02-21

李永華(1961)，女，學(xué)士，教授，主要從事火電機(jī)組節(jié)能理論及節(jié)能技術(shù)方面的研究；

蒲亮(1989)，男，碩士研究生，研究方向?yàn)殡娬驹O(shè)備狀態(tài)監(jiān)測(cè)、控制與信息技術(shù)，E-mail：ncepupuliang@163.com；

周思遠(yuǎn)(1987)，男，碩士，助理工程師，研究方向?yàn)殡姀S運(yùn)行經(jīng)濟(jì)性。