CN116859710A - 高温气冷堆氦气流量控制系统仿真方法及装置 - Google Patents
高温气冷堆氦气流量控制系统仿真方法及装置 Download PDFInfo
- Publication number
- CN116859710A CN116859710A CN202310600268.8A CN202310600268A CN116859710A CN 116859710 A CN116859710 A CN 116859710A CN 202310600268 A CN202310600268 A CN 202310600268A CN 116859710 A CN116859710 A CN 116859710A
- Authority
- CN
- China
- Prior art keywords
- helium
- proportional
- mass flow
- action
- control system
- 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.)
- Pending
Links
- 239000001307 helium Substances 0.000 title claims abstract description 433
- 229910052734 helium Inorganic materials 0.000 title claims abstract description 433
- SWQJXJOGLNCZEY-UHFFFAOYSA-N helium atom Chemical compound [He] SWQJXJOGLNCZEY-UHFFFAOYSA-N 0.000 title claims abstract description 433
- 238000000034 method Methods 0.000 title claims abstract description 115
- 238000004088 simulation Methods 0.000 title claims abstract description 59
- 239000007789 gas Substances 0.000 title claims abstract description 42
- 230000009471 action Effects 0.000 claims abstract description 156
- 230000004044 response Effects 0.000 claims abstract description 95
- 230000008569 process Effects 0.000 claims abstract description 58
- 230000033228 biological regulation Effects 0.000 claims abstract description 13
- 230000010354 integration Effects 0.000 claims description 61
- 238000005259 measurement Methods 0.000 claims description 40
- 230000010355 oscillation Effects 0.000 claims description 11
- 230000000694 effects Effects 0.000 claims description 9
- 230000000087 stabilizing effect Effects 0.000 claims description 4
- 239000003607 modifier Substances 0.000 claims 1
- 230000009467 reduction Effects 0.000 description 17
- 238000012360 testing method Methods 0.000 description 15
- 238000004458 analytical method Methods 0.000 description 8
- 238000012937 correction Methods 0.000 description 7
- 238000012795 verification Methods 0.000 description 5
- 238000007405 data analysis Methods 0.000 description 4
- 230000007423 decrease Effects 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 230000006872 improvement Effects 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 230000009257 reactivity Effects 0.000 description 2
- 230000008901 benefit Effects 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 238000011217 control strategy Methods 0.000 description 1
- 230000008878 coupling Effects 0.000 description 1
- 238000010168 coupling process Methods 0.000 description 1
- 238000005859 coupling reaction Methods 0.000 description 1
- 230000003247 decreasing effect Effects 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 230000007246 mechanism Effects 0.000 description 1
- 230000001105 regulatory effect Effects 0.000 description 1
- 239000000126 substance Substances 0.000 description 1
Classifications
-
- G—PHYSICS
- G05—CONTROLLING; REGULATING
- G05B—CONTROL OR REGULATING SYSTEMS IN GENERAL; FUNCTIONAL ELEMENTS OF SUCH SYSTEMS; MONITORING OR TESTING ARRANGEMENTS FOR SUCH SYSTEMS OR ELEMENTS
- G05B11/00—Automatic controllers
- G05B11/01—Automatic controllers electric
- G05B11/36—Automatic controllers electric with provision for obtaining particular characteristics, e.g. proportional, integral, differential
- G05B11/42—Automatic controllers electric with provision for obtaining particular characteristics, e.g. proportional, integral, differential for obtaining a characteristic which is both proportional and time-dependent, e.g. P. I., P. I. D.
-
- 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
- Y02E30/30—Nuclear fission reactors
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Automation & Control Theory (AREA)
- Feedback Control In General (AREA)
Abstract
本公开实施例提供一种高温气冷堆氦气流量控制系统仿真方法及装置,在对氦气流量控制器的仿真调试过程中,先对氦气流量控制器进行比例作用,然后再对氦气流量控制器分别进行比例作用和积分作用,根据氦气流量控制系统各过程状态变量的闭环动态响应曲线得到目标比例增益和目标积分时间常数,通过目标比例增益和目标积分时间常数,在实际运行时实现对氦气流量的精准控制。本方法可以提前仿真验证全厂功率自动控制与调节系统控制氦气流量控制系统的调节特性,验证了高温气冷堆示范工程氦气流量控制系统调试方法的合理性和可靠性,可掌握关键系统设备的调试方法及其运行特性,避免工程实际应用调试过程中频繁引入扰动。
Description
技术领域
本公开实施例属于高温气冷堆自动控制技术领域,具体涉及一种高温气冷堆氦气流量控制系统仿真方法及装置。
背景技术
氦气流量控制系统为国内核电站首次应用,探索制定完善的调试方法和策略,采用可行的仿真验证试验,验证高温气冷堆示范工程氦气流量控制系统调试方法的合理性和可靠性,掌握关键系统设备的调试方法及其运行特性,圆满完成高温堆示范工程调试任务,保障高温堆示范工程安全可靠运行,助力高温堆的商业化推广。
由于高温气冷堆示范工程是一个多输入多输出的复杂大系统,各控制量与被调量之间都存在紧密的耦合关系,控制策略至执行机构都是首次应用于工程实践,氦气流量控制系统的调试工作将是全新的挑战。
针对上述问题,有必要提出一种设计合理且有效解决上述问题的高温气冷堆氦气流量控制系统仿真方法及装置。
发明内容
本公开实施例旨在至少解决现有技术中存在的技术问题之一,提供一种高温气冷堆氦气流量控制系统仿真方法及装置。
本公开实施例的一方面提供一种高温气冷堆氦气流量控制系统仿真方法,所述方法包括:
将氦气流量控制器设置于自动调节状态,并对所述氦气流量控制器进行比例作用;
保持氦气质量流量测量值不变,将氦气质量流量测量值反馈回路断开,并对所述氦气流量控制器施加第一氦气质量流量设定值扰动信号;
根据所述第一氦气质量流量设定值扰动信号,获得在不同比例增益下所述氦气流量控制系统各过程状态变量在比例作用下的开环响应曲线;
将所述氦气流量控制器的所述氦气质量流量测量值反馈回路闭合,并对所述氦气流量控制器施加第二氦气质量流量设定值扰动信号;
根据所述第二氦气质量流量设定值扰动信号,获得在不同比例增益下所述氦气流量控制系统各过程状态变量在比例作用下的闭环响应曲线;
根据在不同比例增益下所述在比例作用下的闭环响应曲线,选择目标比例增益;
对所述氦气流量控制器分别进行比例作用和积分作用,并对所述氦气流量控制器施加第三氦气质量流量设定值扰动信号,其中,所述氦气流量控制器的比例增益为所述目标比例增益;
根据所述第三氦气质量流量设定值扰动信号,获得在不同积分时间常数下所述氦气流量控制系统各过程状态变量在比例作用和积分作用下的闭环动态响应曲线;
根据在不同积分时间常数下所述在比例作用和积分作用下的闭环动态响应曲线,选择目标积分时间常数。
可选的,所述根据所述第二氦气质量流量设定值扰动信号,获得在不同比例增益下所述氦气流量控制系统各过程状态变量在比例作用下的闭环响应曲线,包括:
根据所述第二氦气质量流量设定值扰动信号,调整所述氦气流量控制器的比例增益,获得在不同比例增益下所述氦气流量控制系统各过程状态变量从初始到稳态的闭环动态响应曲线,以得到所述氦气流量控制系统在比例作用下的闭环动态响应曲线;其中,
所述第二氦气质量流量设定值扰动信号为幅度范围为0.5kg/s~1.5kg/s氦气质量流量设定值负/正阶跃扰动信号。
可选的,所述调整所述氦气流量控制器的比例增益,包括:
先将所述比例增益设定为一个较小值,并逐渐增大所述比例增益;其中,所述比例增益的调整范围为1~50。
可选的,所述根据在不同比例增益下所述在比例作用下的闭环响应曲线,选择目标比例增益,包括:
根据所述在比例作用下的闭环动态响应曲线的振荡大小、稳定时间短长短和偏差大小中的至少一者,在各不同比例增益中选择所述目标比例增益。
可选的,所述根据所述第三氦气质量流量设定值扰动信号,获得在不同积分时间常数下所述氦气流量控制系统各过程状态变量在比例作用和积分作用下的闭环动态响应曲线,包括:
根据所述第三氦气质量流量设定值扰动信号,调整所述氦气流量控制器的积分时间常数,获得在不同积分时间常数下所述氦气流量控制系统各过程状态变量从初始到稳态的闭环动态响应曲线,以得到所述氦气流量控制系统在比例作用和积分作用下的闭环动态响应曲线;其中,
所述第三氦气质量流量设定值扰动信号为幅度范围为0.5kg/s~1.5kg/s的氦气质量流量设定值正/负阶跃扰动信号。
可选的,所述调整所述氦气流量控制器的积分时间常数,包括:
先将所述积分时间常数设定为一个较大值,并逐渐减小所述积分时间常数;其中,所述积分时间常数的调整范围为10~0.1。
可选的,所述根据在不同积分时间常数下所述在比例作用和积分作用下的闭环动态响应曲线,选择目标积分时间常数,包括:
根据在在比例作用和积分作用下的闭环动态响应曲线的振荡大小、稳定时间短长短和偏差大小中的至少一者,在各不同积分时间常数中选择所述目标积分时间常数。
可选的,所述对所述氦气流量控制器施加第一氦气质量流量设定值扰动信号,包括:
对所述氦气流量控制器施加幅值范围为0.5kg/s~1.5kg/s的氦气质量流量设定值正/负阶跃扰动信号。
可选的,所述氦气流量控制系统各过程状态变量包括主氦风机转速设定值及其修正量、主氦风机转速、一回路氦气质量流量、核功率、反应堆出口热氦温度、蒸发器出口蒸汽温度和输出热功率中的至少一者。
本公开实施例的另一方面提供一种高温气冷堆氦气流量控制系统的仿真装置,所述装置包括:
第一设定模块,用于将氦气流量控制器设置于自动调节状态,并对所述氦气流量控制器进行比例作用;
第二设定模块,用于保持氦气质量流量测量值不变,将氦气质量流量测量值反馈回路断开,并对所述氦气流量控制器施加第一氦气质量流量设定值扰动信号;
第一响应曲线获取模块,用于根据所述第一氦气质量流量设定值扰动信号,获得在不同比例增益下所述氦气流量控制系统各过程状态变量在比例作用下的开环响应曲线;
第二设定模块,还用于将所述氦气流量控制器的所述氦气质量流量测量值反馈回路闭合,并对所述氦气流量控制器施加第二氦气质量流量设定值扰动信号;
第二响应曲线获取模块,用于根据所述第二氦气质量流量设定值扰动信号,获得在不同比例增益下所述氦气流量控制系统各过程状态变量在比例作用下的闭环响应曲线;
第一选择模块,用于根据在不同比例增益下所述在比例作用下的闭环响应曲线,选择目标比例增益;
第三设定模块,用于对所述氦气流量控制器分别进行比例作用和积分作用,并对所述氦气流量控制器施加第三氦气质量流量设定值扰动信号,其中,所述氦气流量控制器的比例增益为所述目标比例增益;
第三响应曲线获取模块,用于根据所述第三氦气质量流量设定值扰动信号,获得在不同积分时间常数下所述氦气流量控制系统各过程状态变量在比例作用和积分作用下的闭环动态响应曲线;
第二选择模块,用于根据在不同积分时间常数下所述在比例作用和积分作用下的闭环动态响应曲线,选择目标积分时间常数。
本公开实施例的高温气冷堆氦气流量控制系统的仿真方法及装置,在对氦气流量控制器的仿真调试过程中,先对氦气流量控制器进行比例作用,然后再对氦气流量控制器分别进行比例作用和积分作用,根据氦气流量控制系统各过程状态变量的闭环动态响应曲线得到目标比例增益和目标积分时间常数,通过目标比例增益和目标积分时间常数,在实际运行时实现对氦气流量的精准控制。本公开实施例可以提前仿真验证全厂功率自动控制与调节系统控制氦气流量控制系统的调节特性,验证了高温气冷堆示范工程氦气流量控制系统调试方法的合理性和可靠性,可掌握关键系统设备的调试方法及其运行特性,避免工程实际应用调试过程中频繁引入扰动,使高温气冷堆可以根据调试结果进行安全可靠的运行。
附图说明
图1为本公开实施例中一实施例的一种高温气冷堆核氦气流量控制系统的氦气流量控制器的工作流程示意图;
图2为本公开实施例中另一实施例的一种高温气冷堆氦气流量控制系统的仿真方法的流程示意图;
图3为本公开实施例中另一实施例中一种高温气冷堆氦气流量控制系统的仿真方法的氦气流量控制系统闭环试验比例作用下氦气质量流量仿真曲线;
图4为本公开实施例中另一实施例中一种高温气冷堆氦气流量控制系统的仿真方法的仿真方法的氦气流量控制系统闭环试验比例作用下核功率仿真曲线;
图5为本公开实施例中另一实施例中一种高温气冷堆氦气流量控制系统的仿真方法的氦气流量控制系统闭环试验比例作用下反应堆出口热氦温度仿真曲线;
图6为本公开实施例中另一实施例中一种高温气冷堆氦气流量控制系统的仿真方法的氦气流量控制系统闭环试验比例作用下蒸发器出口蒸汽温度仿真曲线;
图7为本公开实施例中另一实施例中一种高温气冷堆氦气流量控制系统的仿真方法的氦气流量控制系统闭环试验比例作用下输出热功率仿真曲线;
图8为本公开实施例中另一实施例中一种高温气冷堆氦气流量控制系统的仿真方法的氦气流量控制系统闭环试验比例作用和积分作用下氦气质量流量仿真曲线;
图9为本公开实施例中另一实施例中一种高温气冷堆氦气流量控制系统的仿真方法的氦气流量控制系统闭环试验比例作用和积分作用下核功率仿真曲线;
图10为本公开实施例中另一实施例中一种高温气冷堆氦气流量控制系统的仿真方法的氦气流量控制系统闭环试验比例作用和积分作用下反应堆出口热氦温度仿真曲线;
图11为本公开实施例中另一实施例中一种高温气冷堆氦气流量控制系统的仿真方法的氦气流量控制系统闭环试验比例作用和积分作用下蒸发器出口蒸汽温度仿真曲线;
图12为本公开实施例中另一实施例中一种高温气冷堆氦气流量控制系统的仿真方法的氦气流量控制系统闭环试验比例作用和积分作用下输出热功率仿真曲线;
图13为本公开实施例中另一实施例中一种高温气冷堆氦气流量控制系统的仿真装置的结构示意图。
具体实施方式
为使本领域技术人员更好地理解本公开实施例的技术方案,下面结合附图和具体实施方式对本公开实施例作进一步详细描述。
如图1所示,示例性的说明了氦气流量控制器的工作流程,也就是氦气流量控制器自动控制的工作原理。具体地,前馈回路为氦气质量流量设定值对应氦风机转速,基本可以确定氦气流量的需求量与氦风机转速是一一对应的;PID控制回路将氦气流量测量值与氦气质量流量设定值进行比较,根据两者的偏差和偏差的变化信号经过预定的自动调节算法产生氦风机转速的指令信号,当氦气流量测量值高于氦气质量流量设定值时,调节器会根据偏差计算出降低氦风机转速的指令;当氦气流量测量值低于氦气质量流量设定值时,调节器会根据偏差计算出增加氦风机转速的指令,最终通过产生氦风机变频器的设定频率信号,通过升高或降低氦风机的转速,来抑制氦气流量的扰动。
如图2所示,本公开实施例的一个方面提供一种高温气冷堆氦气流量控制系统仿真方法S100,所述方法S100包括:
S110、将氦气流量控制器设置于自动调节状态,并对所述氦气流量控制器进行比例作用。
具体地,将氦气流量控制器置于自动调节状态,通过将积分时间常数P设置成无穷大,关闭积分作用,使氦气流量控制器只进行比例作用。
S120、保持氦气质量流量测量值不变,将氦气质量流量测量值反馈回路断开,并对所述氦气流量控制器施加第一氦气质量流量设定值扰动信号。
具体地,对氦气流量控制器进行比例作用后,保持氦气质量流量测量值不变,将氦气质量流量测量值反馈回路断开,并对氦气流量控制器施加幅值范围为0.5kg/s~1.5kg/s的氦气质量流量设定值正/负阶跃扰动信号,氦气质量流量设定值正阶跃扰动信号可以使主氦风机升速,氦气质量流量设定值负阶跃扰动信号可以使主氦风机降速。
优选的,在本实施例中,第一氦气质量流量设定值扰动信号为幅值为1.0kg/s的氦气质量流量设定值正/负阶跃扰动信号。也就是说,在比例作用下,对氦气流量控制器施加1.0kg/s或者-1.0kg/s的氦气质量流量设定值扰动信号。具体地,在本实施例中,在比例作用下,对氦气流量控制器施加-1.0kg/s的氦气质量流量设定值扰动信号。
S130、根据所述第一氦气质量流量设定值扰动信号,获得在不同比例增益下所述氦气流量控制系统各过程状态变量在比例作用下的开环响应曲线。
在本实施例中,根据施加的幅值为1.0kg/s的氦气质量流量设定值正/负阶跃扰动信号,获得在不同比例增益下氦气流量控制系统各过程状态变量在比例作用下的开环响应曲线。
需要说明的是,在本实施例中,氦气流量控制系统各过程状态变量包括主氦风机转速设定值及其修正量、主氦风机转速、一回路氦气质量流量、核功率、反应堆出口热氦温度、蒸发器出口蒸汽温度和输出热功率。
也就是说,在本实施例中,根据施加的幅值为1.0kg/s的氦气质量流量设定值正/负阶跃扰动信号,获得在不同比例增益下氦气流量控制系统中主氦风机转速设定值及其修正量、主氦风机转速、一回路氦气质量流量、核功率、反应堆出口热氦温度、蒸发器出口蒸汽温度和输出热功率等重要过程状态变量从初始到稳态的开环响应曲线。
本实施例中,根据第一氦气质量流量设定值扰动信号,获得在不同比例增益下氦气流量控制系统各过程状态变量在比例作用下的开环响应曲线,通过测量氦气流量控制器的开环特性,验证主氦风机转速和氦气流量的对应关系,可以更准确的得出比例增益参数。
S140、将所述氦气流量控制器的所述氦气质量流量测量值反馈回路闭合,并对所述氦气流量控制器施加第二氦气质量流量设定值扰动信号。
具体地,将氦气流量控制器的氦气质量流量测量值反馈回路闭合,并对氦气流量控制器施加幅度范围为0.5kg/s~1.5kg/s氦气质量流量设定值负/正阶跃扰动信号。优选的,在本实施例中,第二氦气质量流量设定值扰动信号为幅值为1.0kg/s的氦气质量流量设定值负/正阶跃扰动信号。
S150、根据所述第二氦气质量流量设定值扰动信号,获得在不同比例增益下所述氦气流量控制系统各过程状态变量在比例作用下的闭环响应曲线。
具体地,根据对氦气流量控制器施加幅度范围为0.5kg/s~1.5kg/s氦气质量流量设定值负/正阶跃扰动信号,调整氦气流量控制器的比例增益,获得在不同比例增益下氦气流量控制系统各过程状态变量从初始到稳态的闭环动态响应曲线,以得到氦气流量控制系统在比例作用下的闭环动态响应曲线。
需要说明的是,在本实施例中,氦气流量控制系统各过程状态变量包括主氦风机转速设定值及其修正量、主氦风机转速、一回路氦气质量流量、核功率、反应堆出口热氦温度、蒸发器出口蒸汽温度和输出热功率。
示例性的,调整氦气流量控制器的比例增益,具体包括:
先将比例增益P设定为一个较小值,并逐渐增大比例增益P,这样可以消除偏差。其中,比例增益P的调整范围为1~50。具体地,在本实施例中,比例增益P分别为P=1、P=5、P=10和P=50。
进一步具体地,根据对氦气流量控制器施加幅度为1.0kg/s氦气质量流量设定值负/正阶跃扰动信号,调整氦气流量控制器的比例增益,获得比例增益P分别为P=1、P=5、P=10和P=50时,氦气流量控制系统中主氦风机转速设定值及其修正量、主氦风机转速、一回路氦气质量流量、核功率、反应堆出口热氦温度、蒸发器出口蒸汽温度和输出热功率等从初始到稳态的闭环动态响应曲线,以得到氦气流量控制系统在比例作用下的闭环动态响应曲线。
S160、根据在不同比例增益下所述在比例作用下的闭环响应曲线,选择目标比例增益。
具体地,在本步骤中,比较和分析在不同比例增益下的在比例作用下各过程状态变量的闭环动态响应曲线,可使用人眼观察或参数对照的方法比较和分析响应曲线,从而根据比较和分析的结果,找到性能最好的响应曲线,将该响应曲线对应的比例增益作为目标比例增益。
需要说明的是,在比较和分析过程中,可根据单一的过程状态变量的响应曲线进行比较和分析,也可以结合多个过程状态变量的响应曲线进行比较和分析。
示例性的,根据在不同比例增益下在比例作用下的闭环响应曲线,选择目标比例增益,具体包括:
根据在比例作用下的闭环动态响应曲线的振荡大小、稳定时间短长短和偏差大小中的至少一者,在各不同比例增益中选择目标比例增益。进一步具体地,在比较和分析在不同比例增益下的在比例作用下的闭环动态响应曲线时,选择振荡小、稳定时间短、偏差小的响应曲线。也就是说,选择一个既使过渡过程超调量和稳定时间尽可能小,又使偏差最小的比例增益。也就是说,选择一个最优比例增益。
在本实施例中,施加-1kg/s氦气质量流量设定值阶跃扰动信号,积分时间常数I为0,比例增益P分别为P=1、P=5、P=10和P=50,依次对氦气质量流量、核功率、反应堆出口热氦温度、蒸发器出口蒸汽温度、输出热功率进行仿真,综合数据分析,选取比例增益P=5时为目标比例增益。
通过目标比例增益,可以验证高温气冷堆示范工程氦气流量控制系统调试方法的合理性和可靠性,进而掌握关键系统设备的调试方法及其运行特性,在氦气流量控制器实际运行时实现对氦气流量的精准控制,避免了工程实际应用时为机组频繁引入扰动。
S170、对所述氦气流量控制器分别进行比例作用和积分作用,并对所述氦气流量控制器施加第三氦气质量流量设定值扰动信号,其中,所述氦气流量控制器的比例增益为所述目标比例增益。
首先,对氦气流量控制器分别进行比例作用和积分作用。具体地,将目标比例增益作为氦气流量控制器进行比例作用的比例增益,也就是说,比例增益保持为得到的目标比例增益。在本实施例中,氦气流量控制器进行比例作用的比例增益P=5。将积分时间常数I设置成一个较大值,也就是说,分别对氦气流量控制器进行比例作用和积分作用。
然后,对氦气流量控制器施加幅度范围为0.5kg/s~1.5kg/s的氦气质量流量设定值正/负阶跃扰动信号。
优选的,在本实施例中,氦气流量控制器施加幅度为1.0kg/s的氦气质量流量设定值正/负阶跃扰动信号,测量回路的响应曲线,并逐渐减小积分时间常数I。也就是说,在比例作用和积分作用下,对氦气流量控制器施加1.0kg/s或者-1.0kg/s的氦气质量流量设定值阶跃扰动信号。具体地,在本实施例中,在比例作用和积分作用下,对氦气流量控制器施加-1.0kg/s的氦气质量流量设定值阶跃扰动信号。
S180、根据所述第三氦气质量流量设定值扰动信号,获得在不同积分时间常数下所述氦气流量控制系统各过程状态变量在比例作用和积分作用下的闭环动态响应曲线。
具体地,根据施加的幅度范围为0.5kg/s~1.5kg/s的氦气质量流量设定值正/负阶跃扰动信号,调整氦气流量控制器的积分时间常数,获得在不同积分时间常数下所述氦气流量控制系统各过程状态变量从初始到稳态的闭环动态响应曲线,以得到氦气流量控制系统在比例作用和积分作用下的闭环动态响应曲线。
需要说明的是,在本实施例中,氦气流量控制系统各过程状态变量包括主氦风机转速设定值及其修正量、主氦风机转速、一回路氦气质量流量、核功率、反应堆出口热氦温度、蒸发器出口蒸汽温度和输出热功率。
示例性的,所述调整所述氦气流量控制器的积分时间常数,具体包括:
先将积分时间常数I设定为一个较大值,并逐渐减小积分时间常数I。这样一方面可以以减小测量值收敛到设定值的时间,另一方面可以避免积分时间过短所引发的振荡和不稳定。其中,积分时间常数I的调整范围为10~0.1。在本实施例中,积分时间常数I先后调整为I=10、I=1和I=0.1。
也就是说,根据施加的幅度为1.0kg/s的氦气质量流量设定值正/负阶跃扰动信号,调整氦气流量控制器的积分时间常数I,获得在积分时间常数为I=10、I=1和I=0.1时,氦气流量控制系统中主氦风机转速设定值及其修正量、主氦风机转速、一回路氦气质量流量、核功率、反应堆出口热氦温度、蒸发器出口蒸汽温度和输出热功率等重要过程状态变量从初始到稳态的闭环动态响应曲线,以得到氦气流量控制系统在比例作用和积分作用下的闭环动态响应曲线。
S190、根据在不同积分时间常数下所述在比例作用和积分作用下的闭环动态响应曲线,选择目标积分时间常数。
具体地,在本步骤中,比较和分析在不同积分时间常数下核功率控制系统在比例作用和积分作用下的闭环动态响应曲线时,可使用人眼观察或参数对照的方法比较和分析响应曲线,从而根据比较和分析的结果,找到性能最好的响应曲线,将该响应曲线对应的积分时间常数作为目标积分时间常数。
需要说明的是,在比较和分析过程中,可根据单一的过程状态变量的响应曲线进行比较和分析,也可以结合多个过程状态变量的响应曲线进行比较和分析。
示例性的,根据在不同积分时间常数下所述在比例作用和积分作用下的闭环动态响应曲线,选择目标积分时间常数,具体包括:
根据在比例作用和积分作用下的闭环动态响应曲线的振荡大小、稳定时间短长短和偏差大小中的至少一者,在各不同积分时间常数中选择目标积分时间常数。进一步具体地,在比较和分析在不同积分时间常数下的在在比例作用和积分作用下的闭环动态响应曲线时,选择振荡小、稳定时间短、偏差小的响应曲线。也就是说,选择一个既使过渡过程超调量和稳定时间尽可能小,又使偏差最小的积分时间常数。也就是说,选择一个最优积分时间常数。
在本实施例中,施加幅值为1kg/s氦气质量流量设定值阶跃扰动信号;比例增益为P=5,积分时间常数分别为I=0.1;I=1;I=10,依次对氦气质量流量、核功率、反应堆出口热氦温度、蒸发器出口蒸汽温度、输出热功率进行仿真,综合数据分析,选取P=5,I=1时,较合适,曲线整体最平稳。
通过目标积分时间常数,可以验证高温气冷堆示范工程氦气流量控制系统调试方法的合理性和可靠性,进而掌握关键系统设备的调试方法及其运行特性,在氦气流量控制器实际运行时实现对氦气流量的精准控制,避免了工程实际应用时为机组频繁引入扰动。
示例性的,在本实施例中,图3至图7为一种高温气冷堆氦气流量控制系统仿真方法的氦气流量控制系统在比例作用下的仿真曲线图,图3至图7中,参数的设置包括:
以100%RFP功率水平试验为例,1#NSSS模块氦气流量闭环控制回路试验:
扰动为施加-1kg/s氦气质量流量设定值阶跃扰动信号;
积分时间常数I为0,比例作用分别为P=1、P=5、P=10和P=50;
依次截取氦气质量流量、核功率、反应堆出口热氦温度、蒸发器出口蒸汽温度、输出热功率仿真结果:
如图3所示,施加-1kg/s氦气质量流量设定值阶跃扰动信号,设定值低于测量值,比例作用下,氦气流量测量值会下降,从曲线可以看出,比例作用强,氦气流量测量值下降幅度越大,波形幅度也会越大。当P=10时,氦气质量流量测量值的阶跃变化基本为-1kg/s,后有小幅度回调,可验证在该偏差扰动下,P=10的比例系数基本合适。
如图4所示,施加-1kg/s氦气质量流量设定值阶跃扰动信号,设定值低于测量值,比例作用下,氦气流量测量值会下降,同时,核功率会随着氦气流量测量值的下降而下降,后由于反应堆负反应性,核功率会有回调。从曲线可以看出,比例作用强,核功率波动大,振荡幅度。当比例相对较小时,氦气流量的波动对核功率的影响最小,对系统稳定性越好。
如图5所示,施加-1kg/s氦气质量流量设定值阶跃扰动信号,设定值低于测量值,比例作用下,氦气流量测量值会下降,从曲线可以看出,氦气流量下降,核功率下降,氦气温度也随着下降,随着反应堆负反应性作用,氦气温度出现升高,后随着核功率最终下降,氦气温度也趋势稳定,稳定值低于阶跃扰动前的温度数值。
如图6所示,施加-1kg/s氦气质量流量设定值阶跃扰动信号,设定值低于测量值,比例作用下,氦气流量测量值会下降,同时,核功率会随着氦气流量测量值的下降而下降,氦气温度也随着下降,蒸发器出口蒸汽温度会随着核功率的下降而下降。从曲线可以看出,比例作用强,核功率下降作用明显,会导致蒸发器出口蒸汽温度降幅大。从曲线看出,比例越小,氦气流量的波动对蒸发器出口蒸汽温度的影响最小,对系统稳定性越好。
如图7所示:施加-1kg/s氦气质量流量设定值阶跃扰动信号,设定值低于测量值,比例作用下,氦气流量测量值会下降,同时,核功率会随着氦气流量测量值的下降而下降,输出热功率会随着核功率的下降而下降。从曲线可以看出,比例作用强,输出热功率波动大,振荡幅度大。当比例小时,氦气流量的波动对输出热功率的影响最小,对系统稳定性越好。
综上所述,扰动为施加-1kg/s氦气质量流量设定值阶跃扰动信号;积分时间常数I为0,比例增益P分别为P=1、P=5、P=10和P=50,依次对氦气质量流量、核功率、反应堆出口热氦温度、蒸发器出口蒸汽温度、输出热功率进行仿真,综合数据分析,选取比例增益P=5时,作为后续仿真的暂定参数。
需要说明的是,本实施例中的仿真验证仅为举例说明,也就是说,实验参数的选择和分析的过程均为举例说明,并不限制仿真验证仅能使用上述方法。
示例性的,在本实施例中,图8至图12为一种高温气冷堆氦气流量控制系统的仿真方法的核功率控制系统在比例作用和积分作用下的仿真曲线图,图8至图12中,参数的设置包括:
以100%RFP功率水平试验为例,1#NSSS模块氦气流量闭环控制回路试验:
扰动为施加-1kg/s氦气质量流量设定值阶跃扰动信号;
比例增益P=5,积分时间常数I分别为I=0.1、I=1和I=10;
依次截取氦气质量流量、核功率、反应堆出口热氦温度、蒸发器出口蒸汽温度、输出热功率仿真结果:
如图8所示,施加-1kg/s氦气质量流量设定值阶跃扰动信号,设定值低于测量值,比例增益P=5作用下,氦气流量测量值会下降,从曲线可以看出,积分作用合适,氦气流量测量值会快速趋于稳定,积分时间常数I=1时的积分参数基本合适。
如图9所示,施加-1kg/s氦气质量流量设定值阶跃信号,设定值低于测量值,比例作用下,氦气流量测量值会下降,同时,核功率会随着氦气流量测量值的下降而下降,后由于反应堆负反应性,核功率会有回调。从曲线可以看出,积分作用强,核功率下降大,当积分作用相对较弱时,氦气流量的波动对核功率的影响最小,对系统稳定性越好。
如图10所示,施加-1kg/s氦气质量流量设定值阶跃扰动信号,设定值低于测量值,比例作用下,氦气流量测量值会下降,从曲线可以看出,氦气流量下降,核功率下降,氦气温度也随着下降,随着反应堆负反应性作用,氦气温度出现升高,后随着核功率最终下降,氦气温度也趋势稳定,稳定值低于阶跃扰动前的温度数值。当积分作用相对较弱时,氦气流量的波动对氦气温度的影响最小,对系统稳定性越好。
如图11所示,施加-1kg/s氦气质量流量设定值阶跃扰动信号,设定值低于测量值,比例作用下,氦气流量测量值会下降,同时,核功率会随着氦气流量测量值的下降而下降,氦气温度也随着下降,蒸发器出口蒸汽温度会随着核功率的下降而下降。从曲线可以看出,比例作用下,蒸发器出口蒸汽温度会阶跃降低,后续稳定过程由积分作用完成,仿真3个积分参数下,积分作用越强,蒸发器出口蒸汽温度趋于稳定时间越短。
如图12所示,施加-1kg/s氦气质量流量设定值阶跃扰动信号,设定值低于测量值,比例作用下,氦气流量测量值会下降,同时,核功率会随着氦气流量测量值的下降而下降,输出热功率会随着核功率的下降而下降。从曲线可以看出,比例作用下,输出热功率会阶跃降低,后续稳定过程由积分作用完成,仿真3个积分参数下,积分作用越强,输出热功率趋于稳定时间越短。
综上所述,扰动为施加-1kg/s氦气质量流量设定值阶跃信号;比例增益P=5,积分时间常数I分别为I=0.1;I=1;I=10,依次对氦气质量流量、核功率、反应堆出口热氦温度、蒸发器出口蒸汽温度、输出热功率进行仿真,综合数据分析,选取比例增益P=5,积分时间常数I=1时,作为后续仿真的暂定参数。
需要说明的是,本实施例中的仿真验证仅为举例说明,也就是说,实验参数的选择和分析的过程均为举例说明,并不限制仿真验证仅能使用上述方法。
本公开实施例的高温气冷堆氦气流量控制系统的仿真方法,在对氦气流量控制器的仿真调试过程中,先对氦气流量控制器进行比例作用,然后再对氦气流量控制器分别进行比例作用和积分作用,根据氦气流量控制系统各过程状态变量的闭环动态响应曲线得到目标比例增益和目标积分时间常数,通过目标比例增益和目标积分时间常数,在实际运行时实现对氦气流量的精准控制。本公开实施例可以提前仿真验证全厂功率自动控制与调节系统控制氦气流量控制系统的调节特性,验证了高温气冷堆示范工程氦气流量控制系统调试方法的合理性和可靠性,可掌握关键系统设备的调试方法及其运行特性,避免工程实际应用调试过程中频繁引入扰动,使高温气冷堆可以根据调试结果进行安全可靠的运行。
如图13所示,本公开实施例的另一方面提供一种高温气冷堆氦气流量控制系统的仿真装置100,所述装置100包括:
第一设定模块110,用于将氦气流量控制器设置于自动调节状态,并对所述氦气流量控制器进行比例作用。
第二设定模块120,用于保持氦气质量流量测量值不变,将氦气质量流量测量值反馈回路断开,并对所述氦气流量控制器施加第一氦气质量流量设定值扰动信号,
第一响应曲线获取模块130,用于根据所述第一氦气质量流量设定值扰动信号,获得在不同比例增益下所述氦气流量控制系统各过程状态变量在比例作用下的开环响应曲线。
第二设定模块120,还用于将所述氦气流量控制器的所述氦气质量流量测量值反馈回路闭合,并对所述氦气流量控制器施加第二氦气质量流量设定值扰动信号。
第二响应曲线获取模块140,用于根据所述第二氦气质量流量设定值扰动信号,获得在不同比例增益下所述氦气流量控制系统各过程状态变量在比例作用下的闭环响应曲线。
第一选择模块150,用于根据在不同比例增益下所述在比例作用下的闭环响应曲线,选择目标比例增益。
第三设定模块160,用于对所述氦气流量控制器分别进行比例作用和积分作用,并对所述氦气流量控制器施加第三氦气质量流量设定值扰动信号,其中,所述氦气流量控制器的比例增益为所述目标比例增益。
第三响应曲线获取模块170,用于根据所述第三氦气质量流量设定值扰动信号,获得在不同积分时间常数下所述氦气流量控制系统各过程状态变量在比例作用和积分作用下的闭环动态响应曲线。
第二选择模块180,用于根据在不同积分时间常数下所述在比例作用和积分作用下的闭环动态响应曲线,选择目标积分时间常数。
本公开实施例的高温气冷堆氦气流量控制系统的仿真装置,在对氦气流量控制器的仿真调试过程中,先对氦气流量控制器进行比例作用,然后再对氦气流量控制器分别进行比例作用和积分作用,根据氦气流量控制系统各过程状态变量的闭环动态响应曲线得到目标比例增益和目标积分时间常数,通过目标比例增益和目标积分时间常数,在实际运行时实现对氦气流量的精准控制。本装置可以提前仿真验证全厂功率自动控制与调节系统控制氦气流量控制系统的调节特性,验证了高温气冷堆示范工程氦气流量控制系统调试方法的合理性和可靠性,可掌握关键系统设备的调试方法及其运行特性,避免工程实际应用调试过程中频繁引入扰动,使高温气冷堆可以根据调试结果进行安全可靠的运行。
可以理解的是,以上实施方式仅仅是为了说明本公开实施例的原理而采用的示例性实施方式,然而本公开实施例并不局限于此。对于本领域内的普通技术人员而言,在不脱离本公开实施例的精神和实质的情况下,可以做出各种变型和改进,这些变型和改进也视为本公开实施例的保护范围。
Claims (10)
1.一种高温气冷堆氦气流量控制系统仿真方法,其特征在于,所述方法包括:
将氦气流量控制器设置于自动调节状态,并对所述氦气流量控制器进行比例作用;
保持氦气质量流量测量值不变,将氦气质量流量测量值反馈回路断开,并对所述氦气流量控制器施加第一氦气质量流量设定值扰动信号;
根据所述第一氦气质量流量设定值扰动信号,获得在不同比例增益下所述氦气流量控制系统各过程状态变量在比例作用下的开环响应曲线;
将所述氦气流量控制器的所述氦气质量流量测量值反馈回路闭合,并对所述氦气流量控制器施加第二氦气质量流量设定值扰动信号;
根据所述第二氦气质量流量设定值扰动信号,获得在不同比例增益下所述氦气流量控制系统各过程状态变量在比例作用下的闭环响应曲线;
根据在不同比例增益下所述在比例作用下的闭环响应曲线,选择目标比例增益;
对所述氦气流量控制器分别进行比例作用和积分作用,并对所述氦气流量控制器施加第三氦气质量流量设定值扰动信号,其中,所述氦气流量控制器的比例增益为所述目标比例增益;
根据所述第三氦气质量流量设定值扰动信号,获得在不同积分时间常数下所述氦气流量控制系统各过程状态变量在比例作用和积分作用下的闭环动态响应曲线;
根据在不同积分时间常数下所述在比例作用和积分作用下的闭环动态响应曲线,选择目标积分时间常数。
2.根据权利要求1所述的方法,其特征在于,所述根据所述第二氦气质量流量设定值扰动信号,获得在不同比例增益下所述氦气流量控制系统各过程状态变量在比例作用下的闭环响应曲线,包括:
根据所述第二氦气质量流量设定值扰动信号,调整所述氦气流量控制器的比例增益,获得在不同比例增益下所述氦气流量控制系统各过程状态变量从初始到稳态的闭环动态响应曲线,以得到所述氦气流量控制系统在比例作用下的闭环动态响应曲线;其中,
所述第二氦气质量流量设定值扰动信号为幅度范围为0.5kg/s~1.5kg/s氦气质量流量设定值负/正阶跃扰动信号。
3.根据权利要求2所述的方法,其特征在于,所述调整所述氦气流量控制器的比例增益,包括:
先将所述比例增益设定为一个较小值,并逐渐增大所述比例增益;其中,所述比例增益的调整范围为1~50。
4.根据权利要求1至3任一项所述的方法,其特征在于,所述根据在不同比例增益下所述在比例作用下的闭环响应曲线,选择目标比例增益,包括:
根据所述在比例作用下的闭环动态响应曲线的振荡大小、稳定时间短长短和偏差大小中的至少一者,在各不同比例增益中选择所述目标比例增益。
5.根据权利要求1至3任一项所述的方法,其特征在于,所述根据所述第三氦气质量流量设定值扰动信号,获得在不同积分时间常数下所述氦气流量控制系统各过程状态变量在比例作用和积分作用下的闭环动态响应曲线,包括:
根据所述第三氦气质量流量设定值扰动信号,调整所述氦气流量控制器的积分时间常数,获得在不同积分时间常数下所述氦气流量控制系统各过程状态变量从初始到稳态的闭环动态响应曲线,以得到所述氦气流量控制系统在比例作用和积分作用下的闭环动态响应曲线;其中,
所述第三氦气质量流量设定值扰动信号为幅度范围为0.5kg/s~1.5kg/s的氦气质量流量设定值正/负阶跃扰动信号。
6.根据权利要求5所述的方法,其特征在于,所述调整所述氦气流量控制器的积分时间常数,包括:
先将所述积分时间常数设定为一个较大值,并逐渐减小所述积分时间常数;其中,所述积分时间常数的调整范围为10~0.1。
7.根据权利要求1至3任一项所述的方法,其特征在于,所述根据在不同积分时间常数下所述在比例作用和积分作用下的闭环动态响应曲线,选择目标积分时间常数,包括:
根据在在比例作用和积分作用下的闭环动态响应曲线的振荡大小、稳定时间短长短和偏差大小中的至少一者,在各不同积分时间常数中选择所述目标积分时间常数。
8.根据权利要求1至3任一项所述的方法,其特征在于,所述对所述氦气流量控制器施加第一氦气质量流量设定值扰动信号,包括:
对所述氦气流量控制器施加幅值范围为0.5kg/s~1.5kg/s的氦气质量流量设定值正/负阶跃扰动信号。
9.根据权利要求1至3任一项所述的方法,其特征在于,所述氦气流量控制系统各过程状态变量包括主氦风机转速设定值及其修正量、主氦风机转速、一回路氦气质量流量、核功率、反应堆出口热氦温度、蒸发器出口蒸汽温度和输出热功率中的至少一者。
10.一种高温气冷堆氦气流量控制系统的仿真装置,其特征在于,所述装置包括:
第一设定模块,用于将氦气流量控制器设置于自动调节状态,并对所述氦气流量控制器进行比例作用;
第二设定模块,用于保持氦气质量流量测量值不变,将氦气质量流量测量值反馈回路断开,并对所述氦气流量控制器施加第一氦气质量流量设定值扰动信号;
第一响应曲线获取模块,用于根据所述第一氦气质量流量设定值扰动信号,获得在不同比例增益下所述氦气流量控制系统各过程状态变量在比例作用下的开环响应曲线;
第二设定模块,还用于将所述氦气流量控制器的所述氦气质量流量测量值反馈回路闭合,并对所述氦气流量控制器施加第二氦气质量流量设定值扰动信号;
第二响应曲线获取模块,用于根据所述第二氦气质量流量设定值扰动信号,获得在不同比例增益下所述氦气流量控制系统各过程状态变量在比例作用下的闭环响应曲线;
第一选择模块,用于根据在不同比例增益下所述在比例作用下的闭环响应曲线,选择目标比例增益;
第三设定模块,用于对所述氦气流量控制器分别进行比例作用和积分作用,并对所述氦气流量控制器施加第三氦气质量流量设定值扰动信号,其中,所述氦气流量控制器的比例增益为所述目标比例增益;
第三响应曲线获取模块,用于根据所述第三氦气质量流量设定值扰动信号,获得在不同积分时间常数下所述氦气流量控制系统各过程状态变量在比例作用和积分作用下的闭环动态响应曲线;
第二选择模块,用于根据在不同积分时间常数下所述在比例作用和积分作用下的闭环动态响应曲线,选择目标积分时间常数。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202310600268.8A CN116859710A (zh) | 2023-05-25 | 2023-05-25 | 高温气冷堆氦气流量控制系统仿真方法及装置 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202310600268.8A CN116859710A (zh) | 2023-05-25 | 2023-05-25 | 高温气冷堆氦气流量控制系统仿真方法及装置 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN116859710A true CN116859710A (zh) | 2023-10-10 |
Family
ID=88229330
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202310600268.8A Pending CN116859710A (zh) | 2023-05-25 | 2023-05-25 | 高温气冷堆氦气流量控制系统仿真方法及装置 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN116859710A (zh) |
-
2023
- 2023-05-25 CN CN202310600268.8A patent/CN116859710A/zh active Pending
Similar Documents
Publication | Publication Date | Title |
---|---|---|
KR0149864B1 (ko) | 터빈 제어기의 컴퓨터 보조운전방법 | |
US5461559A (en) | Hierarchical control system for molecular beam epitaxy | |
CN112382427B (zh) | 核电厂蒸发器的液位控制方法及系统 | |
CN107061132A (zh) | 一种水电站调速器功率闭环调节提高功率调节速度和调节精度的调节系统 | |
CN111637444B (zh) | 一种基于q学习的核电蒸汽发生器水位控制方法 | |
US4912732A (en) | Automatic steam generator control at low power | |
Huang et al. | Auto-tuning of PID controllers for second order unstable process having dead time | |
CN114999687B (zh) | 一种核反应堆热工水力瞬态试验电功率调节方法及系统 | |
CN116859710A (zh) | 高温气冷堆氦气流量控制系统仿真方法及装置 | |
Jiang et al. | Nonlinear steam valve adaptive controller design for the power systems | |
CN103032869B (zh) | 超临界机组汽温观测优化控制方法 | |
Hambali et al. | Reformulated tangent method of various PID controller tuning for air pressure control | |
CN116466570A (zh) | 用于高温气冷堆蒸汽温度控制系统的仿真方法 | |
CN116466569A (zh) | 用于高温气冷堆输出热功率控制系统的仿真方法 | |
CN116627025A (zh) | 高温气冷堆核功率控制系统的仿真方法及装置 | |
CN110994642A (zh) | 励磁增益对凸极发电机同步转矩影响的量化方法及装置 | |
CN113791587A (zh) | 流化床发电机组的协调控制方法及系统 | |
CN103021489B (zh) | 核电站化学与容积控制系统下泄流温度的控制装置及方法 | |
CN114838347A (zh) | 再热汽温的控制方法、装置、存储介质及电子设备 | |
CN115691850A (zh) | 全厂功率协调控制系统氦气温度仿真方法和装置 | |
Trojnar et al. | Simulation of the fuel gas valve PID controller in closed loop system | |
CN113885311B (zh) | 一种基于广义频率法的燃气轮机燃烧室模型闭环辨识方法 | |
CN111694348B (zh) | 一种火电机组调速系统功率反馈性能的测试方法及装置 | |
Gilkey | Prompt critical control of the ACRR using a linear quadratic regulator design | |
CN117039923A (zh) | 一种用于送端电网高频控制的方法及系统 |
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 |