CN109543222B - 一种不同功率水平下压水堆核电机组甩负荷保护仿真方法 - Google Patents
一种不同功率水平下压水堆核电机组甩负荷保护仿真方法 Download PDFInfo
- Publication number
- CN109543222B CN109543222B CN201811223308.7A CN201811223308A CN109543222B CN 109543222 B CN109543222 B CN 109543222B CN 201811223308 A CN201811223308 A CN 201811223308A CN 109543222 B CN109543222 B CN 109543222B
- Authority
- CN
- China
- Prior art keywords
- nuclear power
- load shedding
- valve
- unit
- steam
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
- 238000004088 simulation Methods 0.000 title claims abstract description 55
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 title claims abstract description 42
- 238000000034 method Methods 0.000 title claims abstract description 24
- 238000012360 testing method Methods 0.000 claims abstract description 28
- 239000002826 coolant Substances 0.000 claims description 25
- 230000009471 action Effects 0.000 claims description 18
- 238000012546 transfer Methods 0.000 claims description 15
- 230000002829 reductive effect Effects 0.000 claims description 13
- 230000008859 change Effects 0.000 claims description 9
- 238000005259 measurement Methods 0.000 claims description 7
- 230000001105 regulatory effect Effects 0.000 claims description 6
- 230000008569 process Effects 0.000 claims description 5
- 230000001276 controlling effect Effects 0.000 claims description 4
- 230000003321 amplification Effects 0.000 claims description 3
- 239000000446 fuel Substances 0.000 claims description 3
- 238000012544 monitoring process Methods 0.000 claims description 3
- 238000003199 nucleic acid amplification method Methods 0.000 claims description 3
- 238000012937 correction Methods 0.000 claims description 2
- 230000006870 function Effects 0.000 description 11
- 238000011160 research Methods 0.000 description 10
- 238000010586 diagram Methods 0.000 description 9
- 230000000694 effects Effects 0.000 description 5
- 238000004458 analytical method Methods 0.000 description 4
- 230000000670 limiting effect Effects 0.000 description 3
- 230000001960 triggered effect Effects 0.000 description 3
- 238000013459 approach Methods 0.000 description 2
- 230000033228 biological regulation Effects 0.000 description 2
- 230000006378 damage Effects 0.000 description 2
- 230000007423 decrease Effects 0.000 description 2
- 230000002401 inhibitory effect Effects 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000008901 benefit Effects 0.000 description 1
- 238000009835 boiling Methods 0.000 description 1
- 238000005094 computer simulation Methods 0.000 description 1
- 230000003750 conditioning effect Effects 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 238000001816 cooling Methods 0.000 description 1
- 230000003247 decreasing effect Effects 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 230000007123 defense Effects 0.000 description 1
- 230000004907 flux Effects 0.000 description 1
- 230000010354 integration Effects 0.000 description 1
- 238000013178 mathematical model Methods 0.000 description 1
- 230000000116 mitigating effect Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 239000003758 nuclear fuel Substances 0.000 description 1
- 238000002360 preparation method Methods 0.000 description 1
- 230000001681 protective effect Effects 0.000 description 1
- 230000009257 reactivity Effects 0.000 description 1
- 230000009467 reduction Effects 0.000 description 1
- 230000004044 response Effects 0.000 description 1
- 230000035945 sensitivity Effects 0.000 description 1
- 238000010998 test method Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q10/00—Administration; Management
- G06Q10/20—Administration of product repair or maintenance
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q50/00—Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
- G06Q50/06—Energy or water supply
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02E—REDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
- Y02E30/00—Energy generation of nuclear origin
Landscapes
- Engineering & Computer Science (AREA)
- Business, Economics & Management (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Economics (AREA)
- General Physics & Mathematics (AREA)
- Human Resources & Organizations (AREA)
- Health & Medical Sciences (AREA)
- Marketing (AREA)
- General Business, Economics & Management (AREA)
- Strategic Management (AREA)
- Tourism & Hospitality (AREA)
- General Health & Medical Sciences (AREA)
- Primary Health Care (AREA)
- Water Supply & Treatment (AREA)
- Public Health (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- Entrepreneurship & Innovation (AREA)
- Operations Research (AREA)
- Quality & Reliability (AREA)
- Monitoring And Testing Of Nuclear Reactors (AREA)
- Control Of Turbines (AREA)
Abstract
本发明公开了一种不同功率水平下压水堆核电机组甩负荷保护仿真方法,在通用的压水堆核电机组仿真模型下,加入调速器和汽轮机,再根据实际情况引入旁路蒸汽控制系统和汽门快速关闭系统两种保护装置,在此基础上对多个功率水平下运行的核电机组进行甩负荷仿真试验,作为实际中核电机组甩负荷试验的参考;通过仿真试验,验证了本发明所搭建的压水堆核电机组仿真模型的准确性和实用性,同时验证了旁路蒸汽控制系统和汽门快速系统能够对有效降低甩负荷对机组的冲击,显著提高核电机组的甩负荷性能,有利于核电机组安全稳定运行,并提出了机组甩负荷保护措施与策略。
Description
技术领域
本发明属于甩负荷试验方法技术领域,尤其涉及一种不同功率水平下压水堆核电机组甩负荷保护仿真方法。
背景技术
由于核电机组往往带基荷运行,具有单机容量大、参数低、易超速的特点。一旦发生甩负荷孤岛运行这种极端工况,机组安全稳定将面临较大的威胁。甩负荷特性是衡量发电机组安全稳定运行能力的一项重要指标,较为常见的有机组从100%额定功率甩负荷至厂用电或至空载。而在现实状况下,核电机组并不都是以100%额定功率满发运行,除此之外还存在不同功率水平下运行的工况。例如广西的防城港核电1号机在投运前经历了30%功率水平跳机不跳堆、50%甩空载、50%甩负荷至厂用电等试验;中广核宁德核电2号机组也完成了30%、50%和75%功率水平试验;还有部分核电机组由于各种原因被迫降功率运行或是长期低功率运行。因此,研究核电机组在不同功率水平下的甩负荷过程及其保护措施具有重要的现实意义。
目前已有部分涉及机组甩负荷及其保护的研究。廖诗武、曾凯文等提出的火电快速甩负荷机组动态仿真建模,王立地、姚金环提出的FCB功能的成功应用与一种新的实现方案,以及冯伟忠提出的1000MW级火电机组旁路蒸汽控制系统作用及配置都是针对快速甩负荷机组提出了一系列的保护措施并进行仿真试验,但这些文献仅仅局限于火电机组;吴文超、王子奇等提出的AP1000汽轮发电机组甩负荷试验研究,周修铭、干福麟等提出的适用于分析大扰动下机网相互影响的压水堆核电站数学模型研究,王梦、丁剑等提出的考虑核电接入的大电网严重故障下孤网高频问题及协调控制措施,高慧敏、孙蓉等提出的田湾核电机组动态模型研究及其对江苏电网系统扰动的响应分析,以及王宝生、王冬青等提出的压水堆核电厂蒸汽排放控制系统实时仿真研究,这些文献研究了核电机组满发甩负荷至厂用电、甩负荷至空载两种工况,而没有对其他功率水平下运行的核电机组进行试验,也缺少对核电机组内部变量的监测;朱军.提出的某核电厂甩负荷孤岛运行试验及问题分析,以及胡彦令提出的M310核电机组调试前准备过程探究,这两篇文献提出了核电机组多个功率水平甩负荷试验条件和方法,但未进行具体试验研究。目前尚无核电机组在多功率水平甩负荷的模型及相关仿真数据,也缺少对核电机组甩负荷保护的针对性研究。
基于此,本发明提出一种不同功率水平下压水堆核电机组甩负荷保护仿真方法,针对1000MW等级大型压水堆核电机组,基于电力系统分析综合程序PSASP的用户自定义模型建立了核反应堆热力学模型,在此基础上加入调速器、汽轮机,并根据实际情况引入旁路蒸汽模型、快速关闭汽门两种典型的保护装置,对多个功率水平下运行的核电机组进行甩负荷仿真研究。
发明内容
针对现有技术的不足,本发明提供一种不同功率水平下压水堆核电机组甩负荷保护仿真方法。
本发明是通过如下的技术方案来解决上述技术问题的:一种不同功率水平下压水堆核电机组甩负荷保护仿真方法,包括以下几个步骤:
步骤1:压水堆核电机组仿真模型的建立;
将实际核电机组的动力学、热力学过程量化后,经过一系列的简化得到通用的压水堆核电机组仿真模型,具体的压水堆核电机组仿真模型包括堆芯燃料和冷却剂热量传递模型、热线温度模型、冷线温度模型、蒸汽发生器模型、中子动态模型、反应堆功率控制系统、一回路平均温度模型、调速器模型以及汽轮机模型;
步骤2:保护装置的引入;
在所述步骤1的压水堆核电机组仿真模型基础上,引入旁路蒸汽控制系统和汽门快速关闭系统两种保护装置,避免旁通阀频繁开闭以及触发汽门快速关闭系统动作;
步骤3:不同功率水平下甩负荷保护的仿真试验;
在不同功率水平下,采用所述步骤2的压水堆核电机组仿真模型模拟甩负荷保护仿真试验;选取具有代表性的机组参数进行监测,通过压水堆核电机组仿真模型获得其波形,观察各具有代表性机组参数在不同功率水平下甩负荷时的变化情况。
进一步的,所述步骤1中,调速器模型是由调节系统、电液伺服系统和原动机组成。
进一步的,所述调节系统的传递函数表达式为:
其中,K表示转速放大倍数,KP表示PID控制比例环节系数,KD表示PID控制微分环节系数,KI表示PID控制积分环节系数,T1表示转速测量时间常数,TR表示电磁功率测量时间常数,K2表示负荷控制前馈系数,△ω表示机组转速偏差值,Pref表示发电机组的功率整给定值,PE表示电磁功率,PCV表示汽门开度指令值。
进一步的,所述电液伺服系统的传递函数表达式为:
其中,TC表示油动机关闭时间常数,TO表示油动机开启时间常数,T2表示反馈环节时间常数,KP表示PID控制比例环节系数,KD表示PID控制微分环节系数,KI表示PID控制积分环节系数,PCV表示汽门开度指令值,PGV表示电液伺服系统的汽门开度输出值。
进一步的,所述步骤1中,汽轮机模型的传递函数表达式为:
其中,TCH表示蒸汽容积时间常数,TRH表示再热器时间常数,TCO表示交叉管时间常数,FHP表示高压缸占比,FIP表示中压缸占比,FLP表示低压缸占比,PS表示一回路输出的主蒸汽压力,λ表示高压缸功率自然过调系数,μ表示汽门开度,PM表示汽轮机输出的机械功率。
进一步的,所述步骤2中,旁路蒸汽控制系统安装在汽轮机的进汽阀前,能够把一部分蒸汽通过旁路阀门直接导入冷凝器,在一定程度上缓和主蒸汽压力激增的同时降低进汽量,从而减小汽轮机输出的机械功率。
进一步的,所述旁路蒸汽控制系统包括蒸汽发生器、冷凝器、再热器、发电机以及高压缸、中压缸和低压缸;所述蒸汽发生器的蒸汽排出口经主汽门和高压汽门后与所述高压缸连接,所述高压缸与再热器连接;所述再热器经过中压汽门与所述中压缸连接,所述中压缸经连接管道与低压缸连接;所述冷凝器还通过旁通阀与主蒸汽管连接;所述发电机与高压缸、中压缸以及低压缸连接。
进一步的,所述步骤2中,触发汽门快速关闭系统动作的条件为:甩负荷信号、功率偏差越限或者转速越限,这三个触发信号中任何一项均能触发汽门快速关闭系统动作;汽门快速关闭系统通过对汽轮机内部高压汽门和中压汽门的快速关闭,减小机组输出的机械功率从而限制机组转速过快上升;
在机组发生甩负荷、解列或是孤岛运行时,汽门快速关闭系统在功率偏差超过30%阈值的情况下触发动作;功率偏差恢复之后,由转速偏差信号承担对汽门快速关闭系统的控制直到机组转速回到103%额定转速以下,最终调速器重新承担起转速控制的功能。
进一步的,所述步骤3中,具有代表性的机组参数包括冷却剂平均温度、主蒸汽压力、机组转速以及反应堆功率变化量。
与现有技术相比,本发明所提供的不同功率水平下压水堆核电机组甩负荷保护仿真方法,在通用的压水堆核电机组仿真模型下,加入调速器和汽轮机,再根据实际情况引入旁路蒸汽控制系统和汽门快速关闭系统两种保护装置,在此基础上对多个功率水平下运行的核电机组进行甩负荷仿真试验,作为实际中核电机组甩负荷试验的参考;通过仿真试验,验证了本发明所搭建的压水堆核电机组仿真模型的准确性和实用性,同时验证了旁路蒸汽控制系统和汽门快速系统能够对有效降低甩负荷对机组的冲击,显著提高核电机组的甩负荷性能,有利于核电机组安全稳定运行,并提出了机组甩负荷保护措施与策略。
附图说明
为了更清楚地说明本发明的技术方案,下面将对实施例描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本发明的一个实施例,对于本领域普通技术人员来说,在不付出创造性劳动的前提下,还可以根据这些附图获得其他的附图。
图1是本发明压水堆核电机组仿真模型的结构框图;
图2是本发明调节系统的传递函数结构框图;
图3是本发明电液伺服系统的传递函数结构框图;
图4是本发明汽轮机模型的传递函数结构框图;
图5是本发明旁路蒸汽控制系统的结构示意图;
图6是本发明旁路蒸汽控制系统的具体控制模型图;
图7是本发明汽门快速关闭系统的控制逻辑图;
图8是本发明30%PN功率水平负荷阶跃试验,其中,图8(a)为负荷阶跃输入信号PG/MW,图8(b)为汽轮机输出功率PM/MW,图8(c)为二回路主蒸汽压力PS/MPa;
图9是本发明30%、50%、100%PN功率水平甩负荷仿真,其中,图9(a)为机组输出电磁功率PG,图9(b)为汽轮机输出机械功率PM,图9(c)为机组转速ω,图9(d)为主蒸汽压力PS,图9(e)为反应堆堆芯功率变化△N,图9(f)为冷却剂平均温度Tav;
图10是本发明旁路蒸汽控制系统、汽门快速关闭系统作用下的机组甩负荷,其中,图10(a)为冷却剂平均温度Tav,图10(b)为主蒸汽压力PS,图10(c)为机组转速ω,图10(d)为反应堆功率变化量△N;
图11是本发明旁路蒸汽控制系统、汽门快速关闭系统作用下的机组参数波动峰值,其中,图11(a)为转速波动峰值ωmax,图11(b)为冷却剂平均温度波动峰值Tav max,图11(c)为主蒸汽压力波动峰值PS max。
具体实施方式
下面结合本发明实施例中的附图,对本发明中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有作出创造性劳动的前提下所获得的所有其他实施例,都属于本发明保护的范围。
本发明所提供的一种不同功率水平下压水堆核电机组甩负荷保护仿真方法,包括以下几个步骤:
步骤1:压水堆核电机组仿真模型的建立;
传统的压水堆核电站由反应堆堆芯、一回路、二回路以及一些辅助回路构成。一回路包括冷却剂循环和蒸汽循环的回路,也就是核电机组的热力系统;其主要功能在于由通过冷却剂的流动将反应堆输出的能量以热能的形式带出,冷却剂处于高压高沸点的状态下,由冷却机泵带动循环流动。吸收了反应堆热量的冷却剂在蒸汽发生器处释放热能产生蒸汽,经过热交换之后的冷却剂再回流至反应堆芯处,继续吸收反应堆热量以达到密闭循环将核反应热量不断送出的目的。吸收冷却剂热量的高压蒸汽进入汽轮机做功,带动发电机输出功率,做功结束后产生的乏蒸汽经过冷凝器,加热器之后回到蒸汽发生器,实现了完整的一回路循环。核电机组二回路包括汽轮机发电机组、冷凝器、凝结水泵等装置,与传统火电机组相似,不再赘述。
将实际核电机组的动力学、热力学过程量化后,经过一系列的简化得到通用的压水堆核电机组仿真模型,具体的压水堆核电机组仿真模型包括堆芯燃料和冷却剂热量传递模型、热线温度模型、冷线温度模型、蒸汽发生器模型、中子动态模型、反应堆功率控制系统、一回路平均温度模型、调速器模型以及汽轮机模型,如图1所示。
在图1中,△THL和△TCL分别表示热线和冷线温度偏差值;△Tav表示冷却剂平均温度偏差值;ΔTe和△Tθ2分别表示与堆芯热交换前和热交换后的冷却剂温度偏差值;△ρext表示控制棒组引入的反应性偏差值;△N表示中子通量偏差值;△TP表示蒸汽发生器处热交换管的温度偏差值;μ和△μ表示汽门开度实际值和偏差值;PM和△PM分别表示机组输出机械功率的实际值和偏差值。
根据某核电站的实际情况,本发明中核电机组的调速器模型是由调节系统、电液伺服系统和原动机组成,其中,如图2所示,调节系统的传递函数表达式为:
其中,K表示转速放大倍数,KP表示PID控制比例环节系数,KD表示PID控制微分环节系数,KI表示PID控制积分环节系数,T1表示转速测量时间常数,TR表示电磁功率测量时间常数,K2表示负荷控制前馈系数,△ω表示机组转速偏差值,Pref表示发电机组的功率整给定值,PE表示电磁功率,PCV表示汽门开度指令值。
如图3所示,电液伺服系统的传递函数表达式为:
其中,TC表示油动机关闭时间常数,TO表示油动机开启时间常数,T2表示反馈环节时间常数,KP表示PID控制比例环节系数,KD表示PID控制微分环节系数,KI表示PID控制积分环节系数,PCV表示汽门开度指令值,PGV表示电液伺服系统的汽门开度输出值
在本实施例中,调速系统的参数如表1所示:
表1调速系统参数
其中,ε表示调节死区,PIDMAX表示PID环节输出最大值,PIDMIN表示PID环节输出最小值,μmax表示汽门开度最大值,μmin表示汽门开度最小值。
如图4所示,汽轮机模型的传递函数表达式为:
其中,TCH表示蒸汽容积时间常数,TRH表示再热器时间常数,TCO表示交叉管时间常数,FHP表示高压缸占比,FIP表示中压缸占比,FLP表示低压缸占比,PS表示一回路输出的主蒸汽压力,λ表示高压缸功率自然过调系数,μ表示汽门开度,PM表示汽轮机输出的机械功率。
在本实施例中,汽轮机的参数如表2所示:
表2汽轮机参数
步骤2:保护装置的引入;
在步骤1的压水堆核电机组仿真模型基础上,引入旁路蒸汽控制系统和汽门快速关闭系统两种保护装置,避免旁通阀频繁开闭以及触发汽门快速关闭系统动作。
旁路蒸汽控制系统安装在汽轮机的进汽阀处,能够把一部分蒸汽通过旁路阀门直接导入冷凝器,在一定程度上缓和主蒸汽压力激增的同时降低进汽量,从而减小汽轮机输出的机械功率。旁路蒸汽控制系统的旁通阀正常情况下处于关闭状态,主要在以下两种情况下动作:(1)机组负荷出现超过10%的阶跃下降;(2)汽轮机高功率运行时,突然脱扣或是发生甩负荷。
如图5所示,旁路蒸汽控制系统包括蒸汽发生器、冷凝器、再热器、发电机以及高压缸、中压缸和低压缸;蒸汽发生器的蒸汽排出口经主汽门和高压汽门后与高压缸连接,高压缸与再热器连接;再热器经过中压汽门与中压缸连接,中压缸经连接管道与低压缸连接;冷凝器还通过旁通阀与主蒸汽管连接;发电机与高压缸、中压缸以及低压缸连接。
应用于核电机组的旁路蒸汽控制系统主要通过温度控制。冷却剂平均温度Tav经过超前滞后校正环节后与平均温度的定值比较,由产生的偏差信号来控制旁通阀的开度,当偏差值减小到死区以内,旁通阀处于关闭状态,避免了旁通阀频繁开闭。其中,平均温度的定值由汽轮机输出功率决定,两者成线性关系;旁路蒸汽控制系统的具体模型如图6所示。
汽门快速关闭系统通过对汽轮机内部高压汽门和中压汽门的快速关闭,减小机组输出的机械功率从而限制机组转速过快上升。
在机组发生甩负荷、解列或是孤岛运行时,汽门快速关闭装置在功率偏差超过30%阈值的情况下触发动作。功率偏差恢复之后,由转速偏差信号承担对汽门快速关闭系统的控制直到机组转速回到103%额定转速以下,最终调速器重新承担起转速控制的功能。触发汽门快速关闭动作的条件按照动作速度从快到慢依次有甩负荷信号、功率偏差越限、转速越限,三个触发信号有一项就能够触发汽门快速关闭保护动作,其控制逻辑如图7所示。
步骤3:不同功率水平下甩负荷保护的仿真试验;
在不同功率水平下,采用所述步骤2的压水堆核电机组仿真模型模拟甩负荷保护仿真试验;选取具有代表性的机组参数进行监测,通过压水堆核电机组仿真模型获得其波形,观察各具有代表性机组参数在不同功率水平下甩负荷时的变化情况。
30%PN功率水平负荷阶跃试验:
用压水堆核电机组仿真模型模拟以下工况:带有汽门快速关闭系统和旁路蒸汽控制系统保护的机组以30%PN运行,在220s降负荷至15%PN,然后在750s时刻升负荷至27%PN功率水平,升降负荷均以阶跃的形式,得到结果如图8(a)、8(b)和8(c)所示,通过与实测录波的对比,机组功率和主蒸汽压力的变化情况吻合,验证了本发明搭建的压水堆核电机组仿真模型的准确性与实用性。
30%、50%、100%PN功率水平甩负荷仿真分析:
甩负荷对于核电机组的内部参数会造成较大的扰动,过大冲击会危害机组运行,或是触发机组保护导致紧急停堆,如冷却剂平均温度Tav过高会导致冷却剂接近沸腾,所要求的一回路压力过高导致超出稳压器限度;主蒸汽压力PS过高会增加蒸汽发生器和汽轮机的负担;反应堆功率N变化过快会引起反应控制棒的过快动作导致磨损或是意外损坏,减少使用寿命,同时引起过快的堆芯温度变化会危害核燃料的完整性;机组转速ω飞升过快会危害汽轮机转子。
国际原子能机构对这些机组参数有相关的规范要求,用于保证核电机组的安全稳定运行,如表3所示。
表3国际原子能机构规范要求
因此选取上述几个有代表性的机组参数——冷却剂平均温度Tav、主蒸汽压力PS、机组转速ω、反应堆功率变化量△N进行监测,通过仿真模型获得其波形,观察各参数在不同功率水平下甩负荷时的变化情况,如图9所示。
由图9可知,在甩负荷的工况下,各功率水平运行的核电机组输出功率PG瞬间从初始值跌落到5%PN的厂用电,如图9(a)所示;通过调速系统的调节作用,机组输出的机械功率PM缓慢下降,经过约50s到达5%PN,与PG一致,如图9(b)所示;甩负荷初期PM与PG有较大的差值,引起机组的转速ω飞升,如图9(c)所示;机组甩负荷引起主汽门迅速关闭,蒸汽总量变化速率相对缓慢,引起了甩负荷初期主蒸汽压力PS的迅速上升,后随着核反应堆反应程度的下降,蒸汽量减少,主蒸汽压力回落,如图9(d)所示;反应堆堆芯功率在甩负荷初期较快地下降一部分,最后以非常缓慢的速率继续下降至与目标负荷5%PN的厂用电相适应,如图9(e);由于机组输出功率瞬间减小较多,反应堆功率短时间内仍处在原来的水平,一回路通过热交换送出功率减少导致冷却剂平均温度Tav迅速上升,后随着反应堆功率下降,Tav也开始缓慢下降,如图9(f)所示。
总体上看,随着核电机组所处功率水平的上升,其甩负荷至厂用电带来的核电机组各项参数扰动也增大。在核电机组投运前的调试过程中,从小扰动到大扰动、从低功率水平到高功率水平的调试顺序有其合理性。
结合图9和表3,核电机组在30%额定功率甩负荷试验中,各项参数能够保持在规范要求范围内,允许不加入保护措施,但PS为7.6MPa(1.07pu),二回路压力仍然较大;在50%额定功率甩负荷试验中,PS的峰值达到了7.75MPa,超过108%的阈值;在100%功率水平甩负荷情况下,机组转速ω达到1.040pu,达到了触发超速保护的阈值1.03pu,Tav的峰值达到324℃,接近紧急停堆的触发值325℃,PS达到了8.4MPa(1.18pu),严重越限。因此,为满足各种功率水平下核电机组甩负荷工况的安全性要求,需要提出合适的保护策略并投入相应的保护措施。
带旁路、带汽门快速关闭系统核电机组甩负荷仿真:
在压水堆核电机组模型中添加旁路蒸汽排放装置和汽门快速关闭系统,以100%功率水平甩负荷工况为例进行仿真试验,仿真结果如图10所示。
由图10可知,在旁路蒸汽控制系统和汽门快速关闭系统的共同作用下,机组转速ω的峰值从1.04pu降至1.032pu;冷却剂平均温度Tav的上升得到一定程度的抑制,其峰值从324℃下降至314℃;主蒸汽压力PS的峰值从1.18pu抑制到1.105pu;反应堆功率△N的下降速率显著减小。虽然在保护装置的作用下机组参数得到了较大程度的抑制,但是PS的峰值仍然处在停堆保护阈值以上,因此在实际运行中应当尽量避免出现满发机组瞬时甩负荷至厂用电。
通过带保护核电机组100%功率水平甩负荷仿真试验,初步验证了旁路蒸汽控制系统和汽门快关系统对机组的保护作用。
为了更进一步研究旁路蒸汽控制系统和汽门快速系统对核电机组的作用特点,对不同功率水平、不同保护装置投入方案下的核电机组甩负荷特性分别进行试验。选择记录试验过程中机组转速ω、冷却剂平均温度Tav、主蒸汽压力PS的波动峰值;经过统计,结果如图11所示。
其中横坐标A表示不装设机组保护;B表示装设旁路蒸汽控制系统而不装设汽门快速关闭系统;C表示装设汽门快速关闭装置而不装设旁路蒸汽控制系统;D表示同时装设旁路蒸汽控制系统和汽门快速关闭系统。
根据图11可知,旁路蒸汽控制系统和汽门快速关闭系统对核电机组的保护作用具有一定针对性。旁路蒸汽控制系统对机组一回路参数Tav、PS波动有显著的抑制作用,而其对二回路参数汽轮机转速ω的限制作用不明显;汽门快速关闭系统对ω飞升有明显的抑制作用,而对Tav、PS的限制作用不明显。从原理上看,由于旁路蒸汽控制系统通过打开排气阀直接作用于主蒸汽压力PS,又因为PS对循环蒸汽的温度有重要影响,是核电机组热力学部分的重要反馈量,所以旁路蒸汽控制系统对机组一回路有较显著的影响;汽门快速关闭系统直接作用于减小汽门开度μ从而减小汽轮机出力PM,因此能够显著减小功率差值进而抑制机组转速飞升从而保护机组二回路。
基于上述特点,核电机组的甩负荷保护可以采取以下策略:当机组一回路的冷却剂温度过高、堆芯功率和堆芯温度变化过快转速飞升不明显时,优先选择增大旁路阀门的开度或是提升旁路蒸汽控制系统动作的灵敏度;当汽轮机超速严重而一回路参数扰动较小时,优先选择提升汽门关闭的速度,缩短汽门关闭时间。
根据本发明所述的甩负荷保护仿真方法可知:
(1)核电机组在不同功率水平下的甩负荷特性不同,随着所处功率水平的上升,机组内部各重要参数受到的扰动增大。如果在较高功率水平下甩负荷,未装设保护装置的核电机组内部会出现参数越限甚至引发紧急停堆。实际运行中应当尽量避免满功率运行压水堆核电机组甩负荷至厂用电。
(2)旁路蒸汽控制系统和汽门快速关闭系统能够对有效降低甩负荷对机组的冲击,显著提高核电机组的甩负荷性能,有利于核电机组安全稳定运行。
(3)根据试验结果,旁路蒸汽控制系统对核电机组一回路参数的保护作用较为明显,汽门快速关闭系统对机组二回路保护作用较为明显。实际运行中可利用这一特点有针对性地制定合理有效的保护策略。
以上所揭露的仅为本发明的具体实施方式,但本发明的保护范围并不局限于此,任何熟悉本技术领域的技术人员在本发明揭露的技术范围内,可轻易想到变化或变型,都应涵盖在本发明的保护范围之内。
Claims (6)
1.一种不同功率水平下压水堆核电机组甩负荷保护仿真方法,其特征在于,包括以下几个步骤:
步骤1:压水堆核电机组仿真模型的建立;
将实际核电机组的动力学、热力学过程量化后,经过一系列的简化得到通用的压水堆核电机组仿真模型,具体的压水堆核电机组仿真模型包括堆芯燃料和冷却剂热量传递模型、热线温度模型、冷线温度模型、蒸汽发生器模型、中子动态模型、反应堆功率控制系统、一回路平均温度模型、调速器模型以及汽轮机模型;
步骤2:保护装置的引入;
在所述步骤1的压水堆核电机组仿真模型基础上,引入旁路蒸汽控制系统和汽门快速关闭系统两种保护装置;旁路蒸汽控制系统安装在汽轮机的进汽阀前;
所述旁路蒸汽控制系统包括蒸汽发生器、冷凝器、再热器、发电机以及高压缸、中压缸和低压缸;所述蒸汽发生器的蒸汽排出口经主汽门和高压汽门后与所述高压缸连接,所述高压缸与再热器连接;所述再热器经过中压汽门与所述中压缸连接,所述中压缸经连接管道与低压缸连接;所述冷凝器还通过旁通阀与主蒸汽管连接;所述发电机与高压缸、中压缸以及低压缸连接;旁路蒸汽控制系统主要通过温度控制,冷却剂平均温度Tav经过超前滞后校正环节后与平均温度的定值比较,由产生的偏差信号来控制旁通阀的开度,当偏差值减小到死区以内,旁通阀处于关闭状态,避免了旁通阀频繁开闭; 其中,平均温度的定值由汽轮机输出功率决定,两者成线性关系;
触发汽门快速关闭系统动作的条件为:甩负荷信号、功率偏差越限或者转速越限;在机组发生甩负荷、解列或是孤岛运行时,汽门快速关闭装置在功率偏差超过30%阈值的情况下触发动作; 功率偏差恢复之后,由转速偏差信号承担对汽门快速关闭系统的控制直到机组转速回到103%额定转速以下,最终调速器重新承担起转速控制的功能; 触发汽门快速关闭动作的条件按照动作速度从快到慢依次有甩负荷信号、功率偏差越限、转速越限,三个触发信号有一项就能够触发汽门快速关闭保护动作;
步骤3:不同功率水平下甩负荷保护的仿真试验;
在不同功率水平下,采用所述步骤2的压水堆核电机组仿真模型模拟甩负荷保护仿真试验;选取具有代表性的机组参数进行监测,通过压水堆核电机组仿真模型获得其波形,观察各具有代表性机组参数在不同功率水平下甩负荷时的变化情况。
2.如权利要求1所述的一种不同功率水平下压水堆核电机组甩负荷保护仿真方法,其特征在于,所述步骤1中,调速器模型是由调节系统、电液伺服系统和原动机组成。
6.如权利要求1所述的一种不同功率水平下压水堆核电机组甩负荷保护仿真方法,其特征在于,所述步骤3中,具有代表性的机组参数包括冷却剂平均温度、主蒸汽压力、机组转速以及反应堆功率变化量。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811223308.7A CN109543222B (zh) | 2018-10-19 | 2018-10-19 | 一种不同功率水平下压水堆核电机组甩负荷保护仿真方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811223308.7A CN109543222B (zh) | 2018-10-19 | 2018-10-19 | 一种不同功率水平下压水堆核电机组甩负荷保护仿真方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN109543222A CN109543222A (zh) | 2019-03-29 |
CN109543222B true CN109543222B (zh) | 2023-04-18 |
Family
ID=65844295
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201811223308.7A Active CN109543222B (zh) | 2018-10-19 | 2018-10-19 | 一种不同功率水平下压水堆核电机组甩负荷保护仿真方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109543222B (zh) |
Families Citing this family (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110365016B (zh) * | 2019-08-05 | 2021-06-22 | 国网福建省电力有限公司 | 一种基于混合仿真的发电机组一次调频特性的评估方法 |
CN110532713A (zh) * | 2019-09-03 | 2019-12-03 | 贵州电网有限责任公司 | 一种水轮机调速器电液伺服仿真系统 |
CN111192696A (zh) * | 2020-01-06 | 2020-05-22 | 中国原子能科学研究院 | 用于核反应堆的功率控制装置及功率控制方法 |
CN111305914B (zh) * | 2020-02-11 | 2022-08-02 | 西安热工研究院有限公司 | 一种基于能量平衡的核电汽轮机高压缸效率测试方法 |
CN111522253A (zh) * | 2020-04-17 | 2020-08-11 | 广西电网有限责任公司电力科学研究院 | 一种锅炉-汽轮机系统的喷水减温系统仿真方法及装置 |
CN112394639B (zh) * | 2020-11-20 | 2021-08-24 | 西安热工研究院有限公司 | 带有增量式调节功能的核电机组控制棒调节方法及系统 |
CN113985775B (zh) * | 2021-10-27 | 2024-01-02 | 中广核工程有限公司 | 一种考虑快关汽门动作的核电机组控制方法及系统 |
CN114784805A (zh) * | 2022-06-20 | 2022-07-22 | 中国科学院合肥物质科学研究院 | 一种电源系统甩负荷过电压抑制方法 |
CN117270379A (zh) * | 2023-04-11 | 2023-12-22 | 四川大学 | 一种反应堆堆芯功率分数阶PIλDμ控制器控制方法 |
Citations (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4977529A (en) * | 1973-02-23 | 1990-12-11 | Westinghouse Electric Corp. | Training simulator for a nuclear power plant |
CN103089342A (zh) * | 2011-11-04 | 2013-05-08 | 华北电力科学研究院有限责任公司 | 汽轮机及其调速系统模型参数修正方法 |
CN104100307A (zh) * | 2014-07-02 | 2014-10-15 | 中广核工程有限公司 | 核电站汽轮机调节系统调试方法及系统 |
CN104865830A (zh) * | 2015-04-29 | 2015-08-26 | 华北电力大学(保定) | 一种机组负荷双重智能优化控制方法 |
CN105653782A (zh) * | 2015-12-29 | 2016-06-08 | 中国电力科学研究院 | 基于实际调速控制柜的数字水轮发电机甩负荷仿真方法 |
CN106326529A (zh) * | 2016-08-09 | 2017-01-11 | 广东电网有限责任公司电力科学研究院 | 一种水电机组一次调频贡献电量的考核系统及方法 |
CN106528978A (zh) * | 2016-11-02 | 2017-03-22 | 中国电力科学研究院 | 一种先进压水堆核电机组调速系统的动态仿真方法和系统 |
CN106527174A (zh) * | 2015-09-11 | 2017-03-22 | 中国电力科学研究院 | 一种风电机组电压源仿真系统 |
CN106680620A (zh) * | 2016-11-23 | 2017-05-17 | 云南电网有限责任公司电力科学研究院 | 一种新型模拟直流停甩负荷试验方法 |
CN107506297A (zh) * | 2017-07-27 | 2017-12-22 | 北京广利核系统工程有限公司 | 基于仿真平台的反应堆保护系统功能测试方法和装置 |
CN107784168A (zh) * | 2017-10-11 | 2018-03-09 | 上海交通大学 | 核电站半速汽轮机系统的全工况仿真系统 |
Family Cites Families (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US7742555B1 (en) * | 2007-07-13 | 2010-06-22 | Westinghouse Electric Company Llc | Enhanced steam dump (bypass) control system |
JP5725817B2 (ja) * | 2010-11-30 | 2015-05-27 | 三菱日立パワーシステムズ株式会社 | ガスタービン制御装置および発電システム |
CN102279901B (zh) * | 2011-05-17 | 2013-01-09 | 湖北省电力公司电力科学研究院 | 一种针对第三代压水堆核电机组的建模方法 |
CN103559924B (zh) * | 2013-11-11 | 2017-03-15 | 中广核工程有限公司 | 核电站汽轮机控制系统的调试方法、装置及系统 |
US9574511B2 (en) * | 2014-07-24 | 2017-02-21 | Basler Electric Company | System and method for a load anticipation feature and its tuning method for a generating set |
CN104484528A (zh) * | 2014-12-18 | 2015-04-01 | 国家电网公司 | 一种基于pscad的核电机组动力、电气混合仿真平台 |
CN106526519B (zh) * | 2016-10-28 | 2019-12-24 | 黑龙江省电力科学研究院 | 一种基于adpss仿真系统的有功功率变送器暂态特性指标评价方法 |
CN106640229B (zh) * | 2016-12-16 | 2018-09-07 | 中广核工程有限公司 | 对核电厂半速汽轮机保护系统进行调试的系统及方法 |
CN107103167B (zh) * | 2017-06-22 | 2020-04-03 | 山东中实易通集团有限公司 | 针对甩负荷工况的deh调速功能诊断方法及系统 |
-
2018
- 2018-10-19 CN CN201811223308.7A patent/CN109543222B/zh active Active
Patent Citations (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4977529A (en) * | 1973-02-23 | 1990-12-11 | Westinghouse Electric Corp. | Training simulator for a nuclear power plant |
CN103089342A (zh) * | 2011-11-04 | 2013-05-08 | 华北电力科学研究院有限责任公司 | 汽轮机及其调速系统模型参数修正方法 |
CN104100307A (zh) * | 2014-07-02 | 2014-10-15 | 中广核工程有限公司 | 核电站汽轮机调节系统调试方法及系统 |
CN104865830A (zh) * | 2015-04-29 | 2015-08-26 | 华北电力大学(保定) | 一种机组负荷双重智能优化控制方法 |
CN106527174A (zh) * | 2015-09-11 | 2017-03-22 | 中国电力科学研究院 | 一种风电机组电压源仿真系统 |
CN105653782A (zh) * | 2015-12-29 | 2016-06-08 | 中国电力科学研究院 | 基于实际调速控制柜的数字水轮发电机甩负荷仿真方法 |
CN106326529A (zh) * | 2016-08-09 | 2017-01-11 | 广东电网有限责任公司电力科学研究院 | 一种水电机组一次调频贡献电量的考核系统及方法 |
CN106528978A (zh) * | 2016-11-02 | 2017-03-22 | 中国电力科学研究院 | 一种先进压水堆核电机组调速系统的动态仿真方法和系统 |
CN106680620A (zh) * | 2016-11-23 | 2017-05-17 | 云南电网有限责任公司电力科学研究院 | 一种新型模拟直流停甩负荷试验方法 |
CN107506297A (zh) * | 2017-07-27 | 2017-12-22 | 北京广利核系统工程有限公司 | 基于仿真平台的反应堆保护系统功能测试方法和装置 |
CN107784168A (zh) * | 2017-10-11 | 2018-03-09 | 上海交通大学 | 核电站半速汽轮机系统的全工况仿真系统 |
Non-Patent Citations (2)
Title |
---|
基于PSO-SVM的快速甩负荷工况下过热蒸汽温度预测建模;孙灵芳等;《化工自动化及仪表》(第02期);第51-55页 * |
火电快速甩负荷机组动态仿真建模;廖诗武等;《电工技术学报》(第04期);第218-226页 * |
Also Published As
Publication number | Publication date |
---|---|
CN109543222A (zh) | 2019-03-29 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109543222B (zh) | 一种不同功率水平下压水堆核电机组甩负荷保护仿真方法 | |
Ichikawa et al. | Light water reactor plant modeling for power system dynamics simulation | |
CN107084011A (zh) | 一种核电汽轮机快速冷却方法 | |
Baek et al. | Transient analyses of the S-CO2 Cycle coupled to PWR for nuclear marine propulsion | |
CN109098788B (zh) | 压水堆核电机组超速保护控制优化系统及优化控制方法 | |
CN214660397U (zh) | 双低压缸汽轮机低负荷运行时提升供热抽汽参数的装置 | |
CN111262280A (zh) | 一种用于压水堆核电机组一次调频的建模及分析方法 | |
Tanim et al. | Analysis of various types of possible fault and consequences in VVER-1200 using PCTRAN | |
Zhang et al. | Assessments of water ingress accidents in a modular high-temperature gas-cooled reactor | |
Sheng et al. | Research on load rejection protection of PWR power plants at different power levels | |
Volkanovski et al. | Electric power frequency and nuclear safety-Subsynchronous resonance case study | |
CN109543223A (zh) | 一种核电机组外特性分析方法 | |
Jang et al. | Applications of integrated safety analysis methodology to reload safety evaluation | |
Khan et al. | Modeling of a Turbine System for Digital Instrumentation and Control System Design | |
WO2024073994A1 (zh) | 反应堆事故模拟试验装置及方法 | |
Tanaka et al. | A 10% cold leg break test at ROSA-IV large scale test facility | |
Qu et al. | Optimization of rapid power reduction (RPR) design for advanced pressurized water reactor | |
Borysenko et al. | Study of accelerated unit unloading mode initiated by turbine feed pump trip with tvsa fuel assemblies operation in WWER-1000 | |
Wen | Modeling technology of large-scale half speed unit and its thermal actuator participating in grid frequency modulation | |
Tang et al. | Modeling and Simulation of Nuclear Power Units for Primary Frequency Regulation Research | |
Wang et al. | ICONE23-2027 ANALYSIS OF NSSS CONTROL SYSTEM RESPONSES DURING LOAD REJECTION TRANSIENT FOR AP1000 | |
Baea et al. | TASS/SMR-S Analysis on a Single-Phase Natural Circulation Phenomenon in SMART-ITL during 3-Different Operating Conditions | |
Hu et al. | Study on the influence of SG structure on the natural circulation characteristics of PWR units | |
Muzafarov et al. | Analysis of the connection of an idle MCP with three operating without a preliminary reduction in the power of the VVER-1000 reactor | |
Hu et al. | Safety Evaluation of Isolated Island Operation on 1000MW Steam Turbine in NPP |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |