項(xiàng)洪一++馮曉露++劉建光++劉道光++于航

DOI：10.16660/j.cnki.1674-098X.2017.14.167

摘 要：隨著環(huán)境的日益惡化，核能以其清潔、高效等優(yōu)點(diǎn)成為我國(guó)能源戰(zhàn)略的重要組成部分之一。目前核能發(fā)電占總發(fā)電量不到3%。到2020年裝機(jī)容量將從當(dāng)前的約23GW增加到58GW，核電站發(fā)展前景廣闊。蒸汽發(fā)生器SG（steam generator）作為核島的主設(shè)備之一，它的液位控制好壞直接影響到核島的安全運(yùn)行，目前核電站SG液位普遍采用PID控制，由于SG是一個(gè)非線性、時(shí)變、非最小相位、存在假水位現(xiàn)象、小穩(wěn)定裕度的復(fù)雜系統(tǒng)，所以在實(shí)際的運(yùn)行中，PID的控制效果比較差。為了解決這些問(wèn)題，本文對(duì)采用模型預(yù)測(cè)控制算法對(duì)SG液位進(jìn)行控制的可行性進(jìn)行了研究：首先根據(jù)現(xiàn)場(chǎng)機(jī)組實(shí)際的運(yùn)行數(shù)據(jù)，應(yīng)用浙江大學(xué)研發(fā)的FRONT-ID多變量閉環(huán)模型辨識(shí)軟件建立SG被控對(duì)象數(shù)學(xué)模型，然后對(duì)現(xiàn)有SG液位控制系統(tǒng)的特點(diǎn)進(jìn)行分析，最后運(yùn)用MATLAB軟件平臺(tái)設(shè)計(jì)和搭建SG的MPC模型預(yù)測(cè)控制系統(tǒng)，并與PID控制系統(tǒng)進(jìn)行仿真比較。仿真結(jié)果表明，MPC的控制效果要明顯好于PID控制，減小了蒸汽發(fā)生器液位的波動(dòng)，使蒸汽發(fā)生器的液位更好地跟蹤設(shè)定值。

關(guān)鍵詞：模型預(yù)測(cè)控制 蒸汽發(fā)生器水位控制 仿真 模型辨識(shí)

中圖分類(lèi)號(hào)：TK 172 文獻(xiàn)標(biāo)識(shí)碼：A 文章編號(hào)：1674-098X（2017）05（b）-0167-04

Model Predictive Control for Nuclear Steam Generator Water Level

Xiang Hongyi Feng Xiaoxia Liu Jianguang Liu Daoguang Yu Hang

（Cnpec adjustable center， Shenzhen Guangdong 518124， China）

Abstract：With the worsening of environment， nuclear power， with its advantages of clean and efficient，become one of the important components of China's energy strategy. nuclear power plant has a very good prospect. Steam generator （SG ） as one of the main equipments of nuclear island， its level control quality directly affects the operation safety of nuclear power unit. At present SG level control in nuclear power plant mainly adopts PID control. SG is a complex system with these characters of nonlinear， time-varying， non-minimum phase， false water level， small stability margin， so in the actual operation， the control effect of PID is relatively poor. In order to solve these problems， the model predictive control method is used to optimize the SG level control in this paper. According to the actual operating data of a nuclear power plant， the FRONT-ID software researched and developed by Zhejiang University is applied to establish the object model of SG. After the analysis on the characteristics and control effect of SG， the contrasting simulation tests between MPC and PID on MATLAB platform are performed in this paper. Simulating results show that the control effect of MPC is better than PID obviously， reducing the level fluctuation of SG and making SG level tracking its set point better.

Key Words：Model predictive control； Steam generator level control； Simulation； Model identification

模型預(yù)測(cè)控制（MPC）可以處理非方系統(tǒng)，它可以根據(jù)優(yōu)先級(jí)[1]控制目標(biāo)排序。1997年Mayuresh， Metler與Morari等[2]應(yīng)用模型預(yù)測(cè)控制（MPC）來(lái)控制蒸汽發(fā)生器的水位。

近年來(lái)國(guó)外很多自動(dòng)化工程公司，已經(jīng)對(duì)模型預(yù)測(cè)控制技術(shù)和策略在生產(chǎn)過(guò)程控制中的應(yīng)用進(jìn)行了廣泛研究，許多預(yù)測(cè)控制軟件已經(jīng)商品化。美國(guó)、加拿大、英國(guó)成立了Setpoint、DMC、Speedup、Simcon、Trieber Controls等公司，取得了巨額利潤(rùn)。

在我國(guó)，預(yù)測(cè)控制在火力發(fā)電站已經(jīng)開(kāi)始應(yīng)用，在主再熱蒸汽溫度控制、燃燒控制、AGC控制和SCR脫硝控制中已經(jīng)取得很好的控制效果，并獲得了很好的經(jīng)濟(jì)效益，目前正處于推廣階段；在世界的壓水堆核電站中，目前應(yīng)用最多的還是傳統(tǒng)的PID控制方法[3]，還沒(méi)有使用預(yù)測(cè)控制算法等其它智能控制算法。但是隨著核電的大力發(fā)展，核電站需要執(zhí)行G模式調(diào)峰運(yùn)行，以保證電網(wǎng)的安全運(yùn)行和電力需求[4-5]。傳統(tǒng)的固定參數(shù)PID控制方法本身不具備優(yōu)化、自適應(yīng)、自學(xué)習(xí)等功能，無(wú)法滿(mǎn)足電網(wǎng)的要求。

該文應(yīng)用MPC控制方法，利用現(xiàn)場(chǎng)機(jī)組實(shí)際運(yùn)行數(shù)據(jù)，通過(guò)浙江大學(xué)研發(fā)的FRONT-ID多變量閉環(huán)辨識(shí)軟件建立SG水位控制的模型函數(shù)，通過(guò)MATLAB平臺(tái)進(jìn)行SG液位控制組態(tài)的搭建，控制組態(tài)有MPC和PID兩種控制方案，同時(shí)進(jìn)行仿真比對(duì)。比對(duì)結(jié)果表明，通過(guò)MPC在SG液位控制中的應(yīng)用，減小了核電站在大瞬態(tài)試驗(yàn)時(shí)蒸汽發(fā)生器水位的波動(dòng)。

1 模型預(yù)測(cè)控制在核電站SG水位控制中的應(yīng)用研究

1.1 模型辨識(shí)

模型辨識(shí)過(guò)程中模型結(jié)構(gòu)的選擇是最重要的環(huán)節(jié)，它決定著模型的質(zhì)量。對(duì)于多變量的模型辨識(shí)，傳統(tǒng)的方案解決模型評(píng)估、階次評(píng)估、閉環(huán)數(shù)據(jù)比較困難。FRONT-ID多變量閉環(huán)辨識(shí)軟件可以解決MPC控制的這些難題。本篇論文的各個(gè)傳遞函數(shù)模型：小閥開(kāi)度對(duì)應(yīng)給水流量；大閥開(kāi)度對(duì)應(yīng)給水流量；給水流量對(duì)應(yīng)液位；蒸汽流量對(duì)應(yīng)液位等傳遞函數(shù)，都是通過(guò)FRONT-ID多變量閉環(huán)辨識(shí)軟件辨識(shí)得到的。

模型辨識(shí)分為大、小功率工況分別進(jìn)行辨識(shí)。大、小功率模型辨識(shí)的分界點(diǎn)為核功率25%Pn（Pn：核功率），當(dāng)核功率大于25%Pn時(shí)為大功率，核功率小于25%Pn時(shí)為小功率。

小功率階段，當(dāng)核功率小于21%Pn時(shí)，給水流量由小閥單獨(dú)控制，大閥保持關(guān)閉；當(dāng)核功率在21%Pn～25%Pn之間時(shí)，小閥逐漸全開(kāi)，大閥由全關(guān)逐漸開(kāi)啟，大小閥共同控制給水流量，所以小功率工況下SG液位控制組態(tài)由大、小閥兩個(gè)副調(diào)節(jié)器共同耦合控制流量。

大功率階段ARE小閥已經(jīng)全開(kāi)，對(duì)流量的影響不變，小閥全開(kāi)對(duì)給水流量的貢獻(xiàn)為350 t/h，在仿真組態(tài)時(shí)將小閥開(kāi)度對(duì)應(yīng)給水流量，簡(jiǎn)化為常數(shù)350 t/h。

1.2 MPC和PID控制組態(tài)搭建

完整工況下的MPC和PID的控制組態(tài)見(jiàn)圖1。上半部分是MPC控制組態(tài)，下半部分是PID控制組態(tài)。

MPC控制回路：主調(diào)節(jié)器（MPC Controller2） 調(diào)節(jié)液位，副調(diào)節(jié)器（PID1）調(diào)節(jié)給水流量。對(duì)于主調(diào)節(jié)器（MPC Controller2），SG的液位設(shè)定值經(jīng)過(guò)給水溫度修正后，作為設(shè)定值連接到MPC控制器的ref端口。液位測(cè)量值（ayeweiceliang是工程現(xiàn)場(chǎng)實(shí)際的液位值）經(jīng)給水溫度修正后，作為測(cè)量值連接到MPC控制器的mo端口。副調(diào)節(jié)器（PID1） 用于控制大閥開(kāi)度，從而控制給水流量。

線性擬合函數(shù)模塊（dgl-xgl MPC linear fitting）和轉(zhuǎn)換開(kāi)關(guān)（Switch3）的用途是實(shí)現(xiàn)核功率大于25%時(shí)，選擇函數(shù)模型為大功率工況下的調(diào)節(jié)器（MPC Controller2）進(jìn)行液位控制，當(dāng)功率小于25%時(shí)，選擇函數(shù)模型為小功率工況下的調(diào)節(jié)器（MPC Controller2經(jīng)過(guò)線性擬合函數(shù)模塊dgl-xgl MPC linear fitting處理后，表征小功率工況下的MPC調(diào)節(jié)器）進(jìn)行液位控制。線性擬合函數(shù)模塊（dgl-xgl MPC linear fitting）的系數(shù)，是通過(guò)將大功率工況下MPC Controller2的控制輸出與小功率工況下的MPC Controller2的輸出經(jīng)過(guò)線性擬合得出的（即另外搭建大、小功率工況下兩頁(yè)組態(tài)，分別仿真出MPC的輸出值，將輸出值進(jìn)行線性擬合，y代表小功率工況下的MPC輸出值，x代表大功率工況下的MPC輸出值，y=0.1774x +0.0862）。

給水流量減去蒸汽流量（蒸汽發(fā)生器凈進(jìn)水量）作為測(cè)量值鏈接到副調(diào)節(jié)器（PID1）。副調(diào)節(jié)器（PID1）的輸出經(jīng)過(guò)大閥開(kāi)度對(duì)應(yīng)給水流量的傳遞函數(shù)（dfkd-gsll ch）后，轉(zhuǎn)換為給水流量，連接到加法器。

小閥控制回路的控制策略是，根據(jù)蒸汽負(fù)荷對(duì)應(yīng)小閥開(kāi)度，從而實(shí)現(xiàn)給水流量控制，但這樣是開(kāi)環(huán)控制，所以將蒸汽負(fù)荷減去了主調(diào)的輸出，從而實(shí)現(xiàn)閉環(huán)控制。具體方案如下：蒸汽流量（azhengqiliuliang from workspace20），通過(guò)函數(shù)（zhengqifuhe）轉(zhuǎn)換為蒸汽負(fù)荷，蒸汽負(fù)荷減去主調(diào)節(jié)器的輸出，經(jīng)過(guò)蒸汽負(fù)荷對(duì)應(yīng)小閥開(kāi)度的線性函數(shù)發(fā)生器（zhengqifuhe-xfkd），開(kāi)環(huán)控制小閥開(kāi)度（蒸汽負(fù)荷0%～25%對(duì)應(yīng)小閥開(kāi)度為0%～100%），經(jīng)過(guò)小閥開(kāi)度對(duì)應(yīng)給水流量的傳遞函數(shù)（xfkd-gsll ch），轉(zhuǎn)換為小閥控制的給水流量，當(dāng)核功率小于25%時(shí)，該流量經(jīng)過(guò)Switch2直接輸出，當(dāng)核功率大于25%時(shí)Switch2選擇350t/h輸出。

大、小閥控制的給水流量求和，計(jì)算出總的給水流量?？偟慕o水流量經(jīng)過(guò)并行的兩個(gè)傳遞函數(shù)，轉(zhuǎn)換為給水流量影響的液位值：大功率工況下給水流量對(duì)應(yīng)液位的傳遞函數(shù)（dgl gs-yw ch）；小功率工況下給水流量對(duì)應(yīng)液位的傳遞函數(shù)（xgl gs-yw ch），兩者通過(guò)Switch1進(jìn)行選擇，當(dāng)功率大于25%時(shí)取大功率一路信號(hào)，當(dāng)功率小于25%時(shí)取小功率一路信號(hào)。同時(shí)蒸汽流量也經(jīng)過(guò)并行的兩個(gè)傳遞函數(shù)，轉(zhuǎn)換為蒸汽流量影響的液位值：大功率工況下蒸汽流量對(duì)應(yīng)液位的傳遞函數(shù)（dgl zq-yw ch）；小功率工況下蒸汽流量對(duì)應(yīng)液位的傳遞函數(shù)（xgl zq-yw ch），兩者通過(guò)Switch5進(jìn)行選擇，當(dāng)功率大于25%時(shí)取大功率一路信號(hào)，當(dāng)功率小于25%時(shí)取小功率一路信號(hào)。兩者求和后代表了現(xiàn)場(chǎng)SG液位的測(cè)量值。

PID控制回路：PID3控制SG液位，PID4控制給水流量。對(duì)于主調(diào)節(jié)器，液位設(shè)定值減去液位測(cè)量值，其差值經(jīng)過(guò)給水溫度的修正后，作為主調(diào)節(jié)器的輸入。副調(diào)節(jié)器負(fù)責(zé)調(diào)節(jié)給水流量。副調(diào)節(jié)器計(jì)算出的大閥開(kāi)度經(jīng)過(guò)傳遞函數(shù)（xfkd-gsll ch1）轉(zhuǎn)換為大閥控制的給水流量值。

小閥控制回路根據(jù)蒸汽流量（azhengqiliuliang from workspace23），通過(guò)函數(shù)（zhengqifuhe1）轉(zhuǎn)換為蒸汽負(fù)荷，蒸汽負(fù)荷減去主調(diào)節(jié)器的輸出，經(jīng)過(guò)蒸汽負(fù)荷對(duì)應(yīng)小閥開(kāi)度的線性函數(shù)發(fā)生器（zhengqifuhe-xfkd1），開(kāi)環(huán)控制小閥開(kāi)度（蒸汽負(fù)荷0%～25%對(duì)應(yīng)小閥開(kāi)度為0%～100%），經(jīng)過(guò)小閥開(kāi)度對(duì)應(yīng)給水流量的傳遞函數(shù)（xfkd-gsll ch），轉(zhuǎn)換為小閥控制的給水流量，當(dāng)核功率小于25%時(shí)，該流量經(jīng)過(guò)Switch4直接輸出，當(dāng)核功率大于25%時(shí)Switch4選擇350 t/h輸出。大、小閥控制的給水流量求和，計(jì)算出總的給水流量?？偟慕o水流量經(jīng)過(guò)并行的兩個(gè)傳遞函數(shù)，轉(zhuǎn)換為給水流量影響的液位值：大功率工況下給水流量對(duì)應(yīng)液位的傳遞函數(shù)（dgl gs-yw ch1）；小功率工況下給水流量對(duì)應(yīng)液位的傳遞函數(shù)（xgl gs-yw ch1），兩者通過(guò)Switch6進(jìn)行選擇，當(dāng)功率大于25%時(shí)取大功率一路信號(hào)，當(dāng)功率小于25%時(shí)取小功率一路信號(hào)。同時(shí)蒸汽流量也經(jīng)過(guò)并行的兩個(gè)傳遞函數(shù)，轉(zhuǎn)換為蒸汽流量影響的液位值：大功率工況下蒸汽流量對(duì)應(yīng)液位的傳遞函數(shù)（dgl zq-yw ch1）；小功率工況下蒸汽流量對(duì)應(yīng)液位的傳遞函數(shù)（xgl zq-yw ch1），兩者通過(guò)Switch7進(jìn)行選擇，當(dāng)功率大于25%時(shí)取大功率一路信號(hào)，當(dāng)功率小于25%時(shí)取小功率一路信號(hào)。兩者求和后代表了現(xiàn)場(chǎng)SG液位的測(cè)量值。

對(duì)于液位的設(shè)定值，當(dāng)蒸汽負(fù)荷大于20%時(shí)，液位的設(shè)定值為常數(shù)0由（setpoint of level）產(chǎn)生。

1.3 仿真結(jié)果

當(dāng)液位設(shè)定值不變機(jī)組功率發(fā)生變化時(shí)的具體控制曲線詳見(jiàn)下面兩圖：

由圖2、3可知，當(dāng)核功率在17%～100%～17%之間變化時(shí)，MPC控制液位的范圍是（-0.401～0.025 m），PID控制液位的范圍是（-0.591～0.172m），在整個(gè)仿真時(shí)間內(nèi)MPC控制的波動(dòng)范圍更小，更穩(wěn)定，控制效果更好。

3 結(jié)語(yǔ)

該篇論文將模型預(yù)測(cè)控制方法應(yīng)用到SG的水位控制中。并通過(guò)仿真試驗(yàn)，仿真了核功率在17%～100%～17%之間變化時(shí)，進(jìn)行負(fù)荷線性、階躍變化的工況下，SG液位分別在MPC和PID控制下的控制效果，并將MPC與PID控制效果進(jìn)行同步比對(duì)：PID控制回路的液位波動(dòng)量范圍是0.763m，MPC控制回路的液位波動(dòng)量是0.426m，MPC控制回路的液位波動(dòng)比PID控制回路的液位波動(dòng)減小了44.17%。比對(duì)結(jié)果表明：MPC的控制效果要明顯好于PID的控制效果，使蒸汽發(fā)生器的液位更好的跟蹤設(shè)定值。

我國(guó)現(xiàn)在大力發(fā)展核電，具有自主品牌的華龍一號(hào)已經(jīng)開(kāi)工建設(shè)，三代核電已經(jīng)開(kāi)工建設(shè)多臺(tái)機(jī)組，4代核電技術(shù)也會(huì)有很大的飛躍，核電領(lǐng)域?qū)φ羝l(fā)生器水位控制的要求也會(huì)越來(lái)越高，希望本篇論文的研究成果會(huì)對(duì)MPC引入到核電控制領(lǐng)域具有指導(dǎo)意義。

參考文獻(xiàn)

[1] 錢(qián)積新，趙均，徐祖華.預(yù)測(cè)控制[M].北京：化學(xué)工業(yè)出版社，2007（9）：7.

[2] Kothare MV，Mettler B，Morari M， et al. Linear parameter varying model predictive control for steam generator level control [J]. Computers & chemical engineering，1997.

[3] 周剛，彭威，張大發(fā).核動(dòng)力蒸汽發(fā)生器水位控制方法分析[J].原子能科學(xué)技術(shù)，2004（38）：19-23.

[4] Na M G， Jung D W. A model predictive controller for Load-following operation of PWR reactor[J].IEEE Transactions on Nuclear Science， 2005，52（4）：1009-1020.

[5] Na M G，Upadhyaya B R. Model predictive control of an SP-100 space reactor using support vector regression and genetic optimization[J].IEEE Transactions on Nuclear Science， 2006，53（4）：2318-2327.