核电站冷却剂平均温度的DMC预测控制方法

潘瑾宜, 杨 婷,2, 钱 虹,2, 罗 非

(1.上海电力大学, 上海 200090; 2.上海市电站自动化技术重点实验室, 上海 200090;3.大连海事大学, 辽宁 大连 116026)

核电站压水堆堆芯冷却剂平均温度控制系统的性能会直接影响二回路产生的蒸汽质量以及功率调节性能,因此对压水堆堆芯冷却剂平均温度控制的要求越来越高[1-2]。

文献[3]设计了一种冷却剂平均温度模糊控制器,并与传统的棒速程序控制的控制效果进行对比,结果表明,模糊控制器在一定程度上减少了冷却剂平均温度跟踪时间,但隶属函数和模糊推理规则等参数的确定依赖人的主观经验,控制精度较低,动态品质较差。文献[4]提出将对角神经网络(Diagonal Recurrent Neural Networks,DRNN)应用于冷却剂平均温度控制,并证明了DRNN的鲁棒性。文献[5]利用冷却剂平均温度的脉冲响应建立非参数模型,设计了内膜控制器,系统的响应速度得到了提升,但出现了较大的超调。

预测控制是近年来被提出的一种适用于过程控制的先进控制算法,可适用于核电站冷却剂平均温度这类大惯性及较强非线性的过程。预测控制包含多种算法,例如:动态矩阵控制(Dynamic Matrix Control,DMC)、模型算法控制(Model Algorithm Control,MAC)、广义预测控制(Generalized Predictive Control,GPC)等。基于冷却剂温度控制难题,文献[6]利用MAC算法,将非线性堆芯方程在固定工况点附近线性化,建立了基于状态空间模型的预测控制器,具有良好的控制效果,但这种线性化的状态空间模型不能完整反映非线性系统的特性。文献[7-8]通过读取仿真平台实时数据的方式,方便快捷地建立了R棒棒位和冷却剂平均温度的非参数模型,然后搭建了一种无需辨识被控对象的DMC控制器,但没有对棒速与冷却剂平均温度的关系进行直接建模。

棒速信号输入下的冷却剂平均温度的变化是一种非自衡过程,然而DMC算法无法直接控制非自衡过程,会产生截断误差。为了改进这一点,许多文献提出了解决方法。文献[9-10]根据带积分环节的被控对象其阶跃响应值线性递增这一特点,修改了DMC算法的预测模型部分,使得控制结果不产生截断误差。但这一控制方法的品质优劣对截断点的选择很敏感。

本文根据压水堆冷却剂平均温度被控系统的非自衡特性及响应需求,设计了一种带反馈环节的DMC控制器。首先利用MATLAB仿真平台建立了冷却剂平均温度非自衡被控对象的机理模型,然后通过加入负反馈环节的手段使其变为自衡被控对象,进行DMC控制器的设计,求解出实际送入执行机构的棒速控制信号,最后将改进的DMC控制效果与PI控制的效果进行对比,以验证其控制的优越性。

1 自衡化冷却剂平均温度模型的建立

1.1 堆芯冷却剂平均温度模型的建立

联合中子物理模型、堆芯热力学模型和温度反馈模型,得到堆芯冷却剂平均温度模型。代入相关参数,利用MATLAB仿真平台求解该方程组,以实现棒速Rs输入下的堆芯冷却剂平均温度Tm的建模。

温度反馈模型为

[ρfuel(t-1)+αf×ΔTf]+

ρm(t-1)+αm×ΔTm

(1)

式中:ρ——反应性;

ρrod——控制棒反应性;

ρI——R棒组的微分价值;

ρfuel——燃料温度反馈反应性;

αf——燃料温度反馈系数;

ΔTf——燃料平均温度变化量;

ρm——冷却剂温度反馈反应性;

αm——冷却剂温度反应性系数;

ΔTm——冷却剂温度变化量。

核反应堆的反馈作用极为复杂,本文主要考虑燃料温度和慢化剂温度变化引起的反应性变化,即ρfuel和ρm。根据文献[11],常规压水堆αf的经验范围为-4×10-5～-1×10-5,αm的经验范围为-50×10-5～-8×10-5,这里取αf=-1×10-5,αm=-8×10-5。实际工程应用时,ρI随着插入的深度而改变,ρI取均一化近似值3.5 pcm[2]。

中子物理模型为

(2)

式中:n——堆芯热功率;

l——中子代时间;

βi——第i组缓发中子有效份额;

λi——第i组缓发中子衰变常数;

Ci——先驱核衰变功率。

表1 6组缓发中子常数

堆芯热力学模型为

(3)

式中:c1——单位长度燃料芯块比热容,c1=πrfρfcf;

Tf——燃料平均温度;

Tcd——包壳平均温度;

R1——芯块与包壳之间的等效热阻,若考虑气隙导热,则R1=1/8πλf+1/2πrfhg;

rf——燃料芯块半径;

qv——燃料的体积释热率,qv=Fn/V;

c2——单位长度包壳比热容,c2=2πrcdρcdCcdδcd;

Tm——冷却剂平均温度;

R2——包壳与冷却剂之间的等效热阻,R2=1/2πrcsh;

cpm——冷却剂比热容;

M——堆芯冷却剂装量;

Tout,Tin——冷却剂出入口温度;

F——燃料元件释热比;

h——对流换热系数;

Acs——总传热面积;

将反应堆内核燃料产生的热量输出到堆外,一般要经历3个过程:燃料及原件包壳的热传导;原件壁面与冷却剂之间的对流传热;冷却剂将热量传到堆外的输热。式(3)分别对应这3个传热过程。其中,冷却剂平均温度为冷却剂进口温度和冷却剂出口温度的平均值,即Tm=(Tin+Tout)/2,冷却剂入口温度Tin为292 ℃[13]。

在代入数据求解微分方程组之前,需要确定Tf,Tcd,Tout,Tm,Ci的初始值。在堆芯处于满功率稳定状态下,Tf,Tcd,Tout,Tm处于平衡状态,Ci也处于平衡状态,即:

(4)

将式(4)代入式(3),即可得Tf,Tcd,Tout,Tm,Ci的初值。

本文以与秦山二期核电站相关的参考文献[13-17],整理得到式(1)～式(4)中涉及的堆芯热工水力参数如表2所示。根据常规压水堆特性,对流换热系数与流动区间、流型、流量、热流密度等诸多因素有关,气隙等效导热系数与燃耗有关,燃料元件释热比也与燃料元件结构有关。为了简化计算,以上参数本文均取工程经验值。

表2 堆芯热工水力参数

MATLAB求解微分方程组有多种不同的算法,如ode45,ode23,ode15s等。其中,ode15s 是基于1～5阶数值微分公式的可变步长、可变阶次求解器,常用于解微分代数方程 (Differential Algebraic Equations,DAE)[18]。由于中子物理模型方程组、堆芯传热方程组是刚性方程组,故本文选择求解刚性方程组的ode15s算法,以实现棒速输入条件下的堆芯冷却剂平均温度被控对象的建模。

1.2 冷却剂平均温度模型的自衡化

设仿真总时间为60 s,初始时间为0 s,仿真步数为1 s。当引入棒速信号为1步/s的阶跃信号时,启动MATLAB仿真,冷却剂平均温度实时值与初始值相比的增量变化曲线如图1所示。

由图1可知,当棒速信号uc为阶跃信号时,ΔTm呈线性增长趋势,故棒速控制下的冷却剂平均温度系统是非自衡系统。若在开环系统加入适当的反馈值K时,则这种带有反馈环节的被控对象是稳定的。本文取反馈系数K=1,u(k)=1步/s。

图1 被控对象自衡化效果

2 DMC预测控制器设计

DMC算法是一种基于被控对象阶跃响应的预测控制算法,适用于渐进稳定的系统。只有被控对象的阶跃响应值趋于稳定的情况下,才能采集稳定区域内的数据点,构建模型向量[19]。本文利用上述构建自衡被控对象的方法,加入反馈通道,将冷却剂平均温度非自衡模型转化为自衡模型。基于DMC预测算法的冷却剂平均温度控制系统结构图如图2所示。

图2 控制系统整体结构示意

首先,完成自衡化被控对象的建立;然后,将预测控制理论应用到自衡化被控系统中,求得预测控制量u;最后,由u计算出实际送入控制棒驱动机构的棒速信号uc,完成非自衡系统预测控制器的设计。

2.1 自衡化被控对象DMC预测控制器的设计

2.1.1 非参数模型的建立

由图1中曲线2取采样周期为1 s的离散值,并取预测时域为180 s,得到自衡化冷却剂平均温度被控系统的阶跃响应。离散化模型记为A1=[a1,a2,a3,…,aN]T,即为非参数模型的值,如表3所示。

表3 自衡化被控对象非参数模型值

2.1.2 预测值计算

2.1.3 滚动优化

(5)

式中:qi,rj——误差加权系数和控制加权系数,分别表示对跟踪误差和控制作用变化的抑制。

通过约束条件J(k)的极值导数为零这一条件,可求得最优控制增量Δu(k)。

2.1.4 反馈校正

2.2 棒速信号求解

由于冷却剂平均温度为非自衡被控对象,上述控制量u仅为作用于自衡化模型的控制信号,因此还需求解出直接作用于控制棒驱动程序的实际棒速控制信号uc。由被控对象自衡化过程可得如下结构关系

uc(k)=u(k-1)+Δu(k)-KΔy(k)=

u(k-1)+Δu(k)-Δy(k)

(6)

其中,K=1。

由式(6)可以直接解得uc。

3 仿真分析

本文基于MATLAB仿真平台,搭建冷却剂平均温度DMC预测算法控制回路,并通过与PI控制进行对比,以验证本文所提出的非自衡系统预测控制方法对被控对象控制性能的改善。设冷却剂平均温度Tm的初始值为306.9 ℃,在t=0 s时给定冷却剂平均温度设定值为309.9 ℃。由于核反应堆二回路负荷变化会对冷却剂平均温度产生扰动,故在180～200 s加入+0.3 K的阶跃信号,用以检测系统的抗扰性能。

设预测算法采样周期Ts为1 s,建模时域N为180 s,预测时域P为10 s,控制时域U为1 s。利用DMC预测控制算法离线计算出控制量u,再解得R棒棒速控制信号uc。

PI控制是工业应用中的一种基本控制方法,控制原理是根据设定值309.9 ℃与实际输出值y构成控制偏差e,将偏差e的比例和积分通过线性组合构成控制量uc,达到消除系统的稳态误差的目的,从而实现对被控对象的控制。本文通过人工调节PI参数,得到控制效果最佳的一组控制参数为kp=0.1,ki=0.2。

启动MATLAB仿真后,两种控制算法下,R棒棒速控制信号uc的变化曲线对比如图3所示。燃料平均温度Tf的对比如图4所示,包壳平均温度Tcd的对比如图5所示,冷却剂平均温度Tm的对比如图6所示。

图3 棒速控制信号uc变化对比示意

图4 燃料平均温度Tf对比示意

图5 包壳平均温度Tcd对比示意

图6 冷却剂平均温度Tm对比示意

由上述仿真曲线可知:预测控制下冷却剂平均温度在60 s左右到达稳定,而PI控制下冷却剂平均温度在100 s左右到达稳定;在+0.3 K的阶跃输出扰动测试下,预测控制比PI控制能更快地使系统冷却剂平均温度恢复设定值。仿真结果表明,将DMC算法应用于棒速控制下的冷却剂平均温度系统,能够较为快速、稳定地使系统达到设定值,有效缩短了系统的响应时间,并使系统具有更好的抗扰能力,在一定程度上改善了系统的控制性能。

4 结 语

针对压水堆冷却剂平均温度系统具有较强的刚性、开环不稳定性、复杂非线性等特点,本文提出了一种带反馈环节的DMC预测控制算法,对冷却剂平均温度进行了控制。首先,利用MATLAB的ode15s算法,求解压水堆堆芯热工水力刚性方程组,建立冷却剂平均温度被控对象模型。然后,为了弥补DMC算法的控制局限性,加入反馈环节将非自衡被控对象自衡化,进行DMC预测控制器设计,最终求解出作用于R棒的实际棒速控制信号。此外,本文还将经典PI控制器的控制结果与DMC预测控制器的控制结果进行了对比,在300 s的仿真时间内,控制效果对比表明,加入DMC预测控制器的冷却剂平均温度被控系统的响应速度更快,对冷却剂平均温度设定值的跟踪能力更强,具有良好的控制性能。在后续的研究中还可以考虑压水堆核电站二回路对冷却剂平均温度的影响,以及这种预测控制器对小型核反应堆等其他堆型的控制效果。