1. 引言
核电厂受堆芯特性的限制，反应堆功率不能总是随着汽轮机负荷快速地变化。汽轮机旁路排放系统(GCTc)可作为一个“人工负荷”，通过向冷凝器排放蒸汽的方式解决核电厂一、二回路间功率不平衡的问题。
CPR1000型核电厂在高功率运行期间，主要使用温度控制模式(T模式)应对甩负荷到厂用电、汽轮机跳闸等大幅、快速变化的瞬态 [1]。在国内的在运核电厂中，华龙一号、AP1000和CPR1000类似，GCTc在高功率下使用T模式 [2] [3] [4]。根据CPR1000运行经验反馈，T模式在寿期初应对大幅甩负荷的能力较强，但是在寿期末由于堆芯中子学参数发生变化，其能力有所下降，甚至有触发反应堆紧急停堆的风险。此外，GCTc控制系统性能提升相关的研究相对较少 [5] [6]，故从提高运行安全可靠性角度出发，有必要设计一套应对甩负荷工况的控制器，以提高GCTc的控制性能。
2. 系统建模
使用由中广核研究院有限公司自主研发的非破口类(Non LOCA)热工水力系统分析程序GINKGO对CPR1000的核蒸汽供应系统(NSSS)以及控制系统进行了建模，GINKGO程序具有完全自主知识产权，并经过模型测试和验证，采用点堆模型、单相流模型以及两相流平衡态和/或非平衡态模型、控制系统模型等模拟压水堆的正常运行瞬态过程，具有准确模拟NSSS、堆芯保护与反应堆控制系统等热工水力现象和控制特性的能力 [7]。
2.1. 反应堆热工水力模型
NSSS模型主要包括：反应堆堆芯、压力容器、主泵、稳压器、一回路管道(包括上充和下泄)、蒸汽发生器(SG)、主蒸汽管道、汽轮机旁排阀(GCTc配置了总共4组阀门，排放容量为额定蒸汽流量的85%左右，控制阀通常具有调制和快速两种开启方式 [1])等。
一回路冷却剂从堆芯带出热量，流出压力容器经热管段进入SG，充分换热后回到冷管段由主泵带入堆芯；给水在SG中被加热为饱和蒸汽汇集于主蒸汽母管，再流入汽轮机或由GCTc进入冷凝器，其中汽轮机和冷凝器设置为边界条件。具体的建模节点如图1所示。
2.2. 反应堆控制系统模型
反应堆控制系统模型主要包括：堆芯功率与平均温度控制、稳压器压力与水位控制、蒸汽发生器(SG)水位控制以及GCTc控制系统，其能够使核电厂应对基本的运行瞬态，而不发生紧急停堆事件。
为了研究与分析GCTc的控制特性，需要筛选包络事件工况进行研究(如额定满功率工况下甩负荷至厂用电)，同时为了模拟变负荷过程中复杂的热工水力现象和控制过程特性，需要使用上述控制系统模型。由于本文研究重点为GCTc控制系统，因此仅对GCTc控制系统进行控制原理分析，其余控制系统模型(堆芯功率与平均温度控制、稳压器压力与水位控制、SG水位控制)详见文献8 [8]。
基于问题研究和分析需要，根据GINKGO程序建模方法，在GINKGO程序中建立了上述各个控制系统的模型。以下对GCTc甩负荷控制器的设计和控制性能进行了研究。
3. GCTc控制原理分析
3.1. T模式控制
T模式应对甩负荷工况的控制原理是：由于汽轮机负荷突降进汽阀门关小，使二回路传热恶化，导致一回路冷却剂平均温度升高。冷却剂实测平均温度与通过汽轮机负荷表征的冷却剂参考平均温度存在偏差(相对于反应堆功率，汽轮机负荷变化速度更快)，其反映了一、二回路功率的失配程度，即T模式利用实测与参考平均温度的偏差进行控制。可以发现：在甩负荷瞬态中，冷却剂平均温度是一个应变量，其对一、二回路功率失配的响应有物理参数变化过程的延迟。
T模式控制过程如图2所示。
3.2. 甩负荷控制器
甩负荷控制器在堆芯功率与平均温度控制系统输出的功率整定值大于零时生效(辅助识别甩负荷工况)，使用汽轮机第一级压力额定值与实测值的偏差作为控制信号，通过函数发生器直接生成阀门开度。可以发现：在发生甩负荷工况时，汽轮机第一级压力紧随汽轮机进汽阀动作，由于压力是一个快速响应参数，预计其对一、二回路功率失配的响应比T模式延迟低，可在一定程度上优化控制系统响应。
甩负荷控制器控制过程见图3所示。
Figure 3. Diagram of load rejection controller
图3. 甩负荷控制器控制示意图
4. 控制性能研究与分析
使用需要GCTc应对的最恶劣的甩负荷工况，满功率下甩负荷到厂用电，对T模式和加入甩负荷控制器的控制性能进行了数值模拟和结果分析，评价标准主要是功率量程中子高变化率停堆(DFDT)裕量、超温\超功率停堆(OT/OPΔT)裕量，其次是GCTc调节阀的响应时间。
4.1. 甩负荷至厂用电瞬态
首先，NSSS系统在满功率额定工况运行50 s，50 s后触发甩负荷瞬态过程，寿期初(BOL) T模式下相关参数的数值模拟曲线见图4~7，对数据结果进行了归一化处理。
机组运行在额定满功率下，电网断路器打开导致汽轮发电机立即丧失电磁转矩，从而汽轮机转速上升。在此过程中，由于主泵切换为厂用电供电，将导致一回路冷却剂流量短时上升。接着汽轮机控制系统通过调节汽轮机入口阀门来调节汽轮机转速，汽轮机入口阀有一个从关闭到再次打开的过程。在此之后，汽轮机维持一定的负荷向厂内辅助系统供电。
汽轮机入口阀门的突然关闭导致蒸汽旁排阀快速打开，并通过控制棒下插使核功率降低。在汽轮机入口阀门关闭后的短时间内，二次侧蒸汽排放能力不足导致反应堆冷却剂系统升温。
Figure 4. Diagram of nuclear power transient trend
图4. 核功率瞬态变化趋势
Figure 5. Diagram of coolant average temperature transient trend
图5. 一回路冷却剂平均温度瞬态趋势
Figure 6. Diagram of SG outlet pressure
图6. SG出口压力瞬态趋势
Figure 7. Diagram of GCTc opening transient trend
图7. GCTc调节阀开度瞬态趋势
在甩负荷至厂用电过程中，反应堆堆芯产生的大部分热量由GCT排放至冷凝器，随着反应堆功率的降低，GCT调节阀开度逐渐减小，直至稳定在新的状态，与NSSS系统热功率相匹配。T模式基本可以达到满功率甩负荷至厂用电瞬态的控制目标，但在冷却剂平均温度极值和GCT调节阀响应速度控制上有一定程度的优化空间。
4.2. 甩负荷控制器性能研究
基于3.2节的甩负荷控制器方案，在计算程序中建立了甩负荷控制器的控制逻辑。使用与4.1节相同的瞬态对其控制性能进行了模拟分析，计算结果主要参数的变化曲线与4.1节相似，不再单独给出，瞬态中关注的主要裕量等参数见表1。
从结果可以发现：在瞬态开始时刻，由于堆芯功率调节的作用堆芯功率首先开始下降，在GCTc调节阀完全发挥功能之前，冷却剂平均温度依然上升叠加较大的慢化剂温度系数，使得寿期末的功率量程中子注量率负变化率高的停堆裕量较低。
在T模式的基础上增加甩负荷控制器，使其直接对汽轮机负荷变化进行响应，在一定程度上加快GCTc调节阀开启速率，可以获得相对较高的功率量程中子注量率负变化率高的停堆裕量。此外，也能改善OT/OPΔT停堆裕量。
Table 1. Comparison of main parameters
表1. 主要结果对比
5. 结论
本文结合CPR1000型核电厂甩负荷瞬态中参数变化的物理特性，研究了基于T模式的甩负荷控制器的性能，结果表明：基于T模式并增加甩负荷控制器后，机组在寿期末应对满功率甩负荷至厂用电瞬态时，GCTc调节阀响应时间加快，功率量程中子注量率负变化率高的停堆裕量增大50%，OT/OPΔT裕量也得到一定程度提升。
本研究成果有助于降低CPR1000型机组在寿期末由于甩负荷瞬态触发反应堆紧急停堆的风险，可为CPR1000型机组后续的GCTc控制优化提供理论支持。