CN102411997A - 高温气冷堆功率控制方法及系统 - Google Patents
高温气冷堆功率控制方法及系统 Download PDFInfo
- Publication number
- CN102411997A CN102411997A CN2011103282716A CN201110328271A CN102411997A CN 102411997 A CN102411997 A CN 102411997A CN 2011103282716 A CN2011103282716 A CN 2011103282716A CN 201110328271 A CN201110328271 A CN 201110328271A CN 102411997 A CN102411997 A CN 102411997A
- Authority
- CN
- China
- Prior art keywords
- core
- hot helium
- flux density
- value
- relative
- 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.)
- Granted
Links
- 238000000034 method Methods 0.000 title claims abstract description 68
- 239000001307 helium Substances 0.000 claims abstract description 101
- 229910052734 helium Inorganic materials 0.000 claims abstract description 101
- SWQJXJOGLNCZEY-UHFFFAOYSA-N helium atom Chemical compound [He] SWQJXJOGLNCZEY-UHFFFAOYSA-N 0.000 claims abstract description 101
- 238000005070 sampling Methods 0.000 claims description 28
- 238000012937 correction Methods 0.000 claims description 12
- 230000003111 delayed effect Effects 0.000 claims description 6
- 238000012986 modification Methods 0.000 abstract description 2
- 230000004048 modification Effects 0.000 abstract description 2
- 230000004907 flux Effects 0.000 abstract 3
- 239000007789 gas Substances 0.000 abstract 3
- 230000001052 transient effect Effects 0.000 description 13
- 238000013461 design Methods 0.000 description 10
- 230000001276 controlling effect Effects 0.000 description 9
- 238000005516 engineering process Methods 0.000 description 6
- 238000011161 development Methods 0.000 description 4
- 238000010586 diagram Methods 0.000 description 4
- 238000005259 measurement Methods 0.000 description 4
- 239000002826 coolant Substances 0.000 description 2
- 239000000463 material Substances 0.000 description 2
- 230000001105 regulatory effect Effects 0.000 description 2
- 238000004088 simulation Methods 0.000 description 2
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 2
- OKTJSMMVPCPJKN-UHFFFAOYSA-N Carbon Chemical compound [C] OKTJSMMVPCPJKN-UHFFFAOYSA-N 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 239000000919 ceramic Substances 0.000 description 1
- 238000012938 design process Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 239000000446 fuel Substances 0.000 description 1
- 229910002804 graphite Inorganic materials 0.000 description 1
- 239000010439 graphite Substances 0.000 description 1
- 230000010354 integration Effects 0.000 description 1
- 238000011160 research Methods 0.000 description 1
Images
Classifications
-
- 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
- Monitoring And Testing Of Nuclear Reactors (AREA)
Abstract
本发明公开了一种高温气冷堆功率控制方法及系统,涉及核能科学与工程技术领域,该方法包括:S1:设置核功率控制器参数和堆芯出口热氦温度控制器参数;S2:所述堆芯出口热氦温度控制器根据堆芯出口热氦温度测量值与所述堆芯出口热氦温度参考值之间的偏差,通过计算获得相对核功率参考值的修正;所述核功率控制器根据由堆芯相对中子通量密度测量值和所述相对堆芯中子通量密度参考值之间的偏差通过计算给出驱动控制棒运动的棒速信号。本发明通过实时获取高温气冷堆的堆芯出口热氦温度测量值和堆芯相对中子通量密度测量值,并对测量值进行相应的计算,获得驱动控制棒运动的棒速信号,以实现高温气冷堆功率控制。
Description
技术领域
本发明涉及核能科学与工程技术领域,特别涉及一种高温气冷堆功率控制方法及系统。
背景技术
由于具有固有安全性和潜在经济竞争力,高温气冷堆已被世界核能界公认为具有第四代核能系统特征的首选堆型之一,发展高温气冷堆技术对于我国可持续发展具有重要意义。高温气冷堆采用陶瓷包覆型燃料元件,以氦气作为冷却剂,以石墨作为慢化剂和结构材料,其固有安全性由低功率密度、瘦长型堆芯和强温度负反馈效应来保证。安全、稳定和高效的运行是高温气冷堆最为关键的性能指标,而功率控制是保证反应堆动态性能并进而保证其稳定高效运行的最重要的技术手段之一。所谓反应堆功率控制就是通过调节控制棒棒位使得堆芯中子通量密度(核功率)和堆芯出口热氦温度满足工程需求。功率控制的本质是增强闭环系统的稳定性,性能良好的功率控制器对于发展高温气冷堆技术具有十分重要的意义。由于高温气冷堆的慢化剂、冷却剂和反射层材料均不同于压水堆,这就导致二者在动态特性上具有显著的差别,从而适用于压水堆的功率控制方法并不适用于高温气冷堆。虽然高温气冷堆是一种极具发展前景的堆型,但目前在高温气冷堆功率控制方法的设计及其参数整定方面的报道却很少。《核动力工程》于2001年公开了“10MW高温气冷堆控制方案研究”,给出了通过数值仿真并凭借经验整定的PID型高温气冷堆功率控制方法,由于该控制方法的参数是凭借经验整定的,因而仅能保证局部稳定性且在性能上具有保守性。《中国智能自动化学术会议论文集》于1999年公开了“HTR-10核功率调节系统模糊控制及其性能的仿真研究”,给出了一种高温气冷堆的模糊PID功率控制方法,且该控制方法由模糊规则和在不同功率水平下凭借经验整定的一组PID控制器组成,它虽能保证闭环系统大范围的稳定性,但这种大范围稳定性是靠增加PID控制器的数量来保证并必然导致其实现过于复杂,从而不适用于实际工程。由此可见,目前的高温气冷堆功率控制方法大都是凭借经验整定的,且在形式上要么是简单的但仅能保证工作点附近闭环稳定性的单个PID控制器,要么是保证大范围闭环稳定性但实现复杂的一组PID控制器及其对应切换规则,即存在形式和性能上的矛盾。
发明内容
(一)要解决的技术问题
本发明要解决的技术问题是:如何解决高温气冷堆功率控制方法存在形式和性能上矛盾的问题。
(二)技术方案
为解决上述技术问题,本发明提供了一种高温气冷堆功率控制方法,包括以下步骤:
S1:设置核功率控制器参数和堆芯出口热氦温度控制器参数;
S2:所述堆芯出口热氦温度控制器根据堆芯出口热氦温度测量值与所述堆芯出口热氦温度参考值之间的偏差,通过计算获得相对核功率参考值的修正;所述核功率控制器根据由堆芯相对中子通量密度测量值和所述堆芯相对中子通量密度参考值之间的偏差通过计算给出驱动控制棒运动的棒速信号,以实现对反应堆功率的控制,所述相对堆芯中子通量密度的参考值为所述相对核功率的参考值和所述修正值相加。
优选地,所述核功率控制器的离散化形式为:
其中,en=ls-nr,KP,1=ΓP,1+ΓP,2,KD,1=ΓD,1+ΓD,2,
σ1、σ2、F1以及F2均为预设的正常数,Ts,n为核功率控制回路的采样周期,Gr为控制棒微分价值,ls为相对堆芯中子通量密度参考值,nr为相对堆芯中子通量密度测量值,en为相对堆芯中子通量密度测量值和所述相对堆芯中子通量密度参考值之间的偏差,en(k)为当前采样时刻所对应的堆芯相对中子通量密度测量值和所述堆芯相对中子通量密度参考值之间的偏差,u1(k)为当前所述核功率控制器的输出量;
所述堆芯出口热氦温度控制器的离散化形式为:
其中,eT=Tout,r-Tout,KP,2=θΓP,2,KD,2=θΓD,2,
Ts,T为堆芯出口热氦温度控制回路的采样周期,Tout为堆芯出口热氦温度测量值,Tout,r为堆芯出口热氦温度参考值,u2(k)为当前相对核功率参考值的修正值,eT为堆芯出口热氦温度测量值与所述堆芯出口热氦温度参考值之间的偏差,eT(k)为当前采样时刻所对应的堆芯出口热氦温度测量值与所述堆芯出口热氦温度参考值之间的偏差。
优选地,所述核功率控制器的离散化形式为:
其中,en=ls-nr,KP,1=ΓP,1+ΓP,2,KD,1=ΓD,1+ΓD,2,
σ1、σ2、F1以及F2均为预设的正常数,κ为预设的正常数、且取值范围为0.001≤κ≤0.01,Ts,n为核功率控制回路的采样周期,Gr为控制棒微分价值,ls为堆芯相对中子通量密度参考值,nr为堆芯相对中子通量密度测量值,en为堆芯相对中子通量密度测量值和所述堆芯相对中子通量密度参考值之间的偏差,en(k)为当前采样时刻所对应的堆芯相对中子通量密度测量值和所述堆芯相对中子通量密度参考值之间的偏差,u1(k)为当前所述核功率控制器的输出量;
所述堆芯出口热氦温度控制器的离散化形式为:
其中,eT=Tout,r-Tout,KP,2=θΓP,2,KD,2=θΓD,2,
θ为预设的正常数,Ts,T为堆芯出口热氦温度控制回路的采样周期,Tout为堆芯出口热氦温度测量值,Tout,r为堆芯出口热氦温度参考值,u2(k)为当前相对核功率参考值的修正,eT为堆芯出口热氦温度测量值与所述堆芯出口热氦温度参考值之间的偏差,eT(k)为当前采样时刻所对应的堆芯出口热氦温度测量值与所述堆芯出口热氦温度参考值之间的偏差。
优选地,控制器参数σ1、σ2、F1、以及F2的取值范围分别为:
其中,Λ为中子每代时间,β为缓发中子份额,vm为控制棒运动的最大速率,ε1、ε2和ε3均为预设的常数且取值范围为大于0小于1。
本发明还公开了一种高温气冷堆功率控制系统,包括:
参数设置模块,用于设置核功率控制器参数和堆芯出口热氦温度控制器参数;
反应堆功率控制模块,用所述堆芯出口热氦温度控制器根据堆芯出口热氦温度测量值与所述堆芯出口热氦温度参考值之间的偏差,通过计算获得相对核功率参考值的修正;所述核功率控制器根据由堆芯相对中子通量密度测量值和所述堆芯相对中子通量密度参考值之间的偏差通过计算给出驱动控制棒运动的棒速信号,以实现对反应堆功率的控制,所述堆芯相对中子通量密度参考值为所述相对核功率参考值和所述修正值相加。
(三)有益效果
本发明通过实时获取高温气冷堆的堆芯出口热氦温度测量值和堆芯相对中子通量密度测量值,并对测量值进行相应的计算,获得驱动控制棒运动的棒速信号,以实现对高温气冷堆热功率的控制,解决了高温气冷堆功率控制方法存在形式和性能上矛盾的问题。
附图说明
图1是高温气冷堆功率控制的回路结构图;
图2是按照本发明一种实施方式的高温气冷堆功率控制方法的流程图;
图3是图2所示的高温气冷堆功率控制方法的参数整定流程图。
具体实施方式
下面结合附图和实施例,对本发明的具体实施方式作进一步详细描述。以下实施例用于说明本发明,但不用来限制本发明的范围。
图1是高温气冷堆功率控制的回路结构图,参照图1,其中,u1(·)为核功率控制器,u2(·)为堆芯出口热氦温度控制器。核功率控制器u1(·)根据由测量给出的堆芯相对中子通量密度nr与其参考值ls之间偏差en,通过计算给出驱动控制棒运动的棒速信号,从而完成对堆芯中子通量密度的调节。热氦温度控制器u2(·)根据堆芯出口热氦温度的测量值Tout与其参考值Tout,r之间的偏差eT,通过计算给出对相对核功率参考值l1的修正l2,从而完成对堆芯出口热氦温度的控制。这里控制棒运动的最大速率是vm,即核功率控制器u1(·)给出的棒速信号需要经过图1中的饱和环节后再用来驱动控制棒运动。
造成高温气冷堆功率控制方法存在形式和性能上矛盾的根源在于无论是控制方法形式还是控制方法参数都采用了基于经验的设计方法,本发明技术方案的总体思路是:通过控制高温气冷堆的堆芯出口热氦温度和堆芯中子通量密度来调节高温气冷堆的热功率,在选择描述过渡过程品质指标函数的基础上,通过保证指标函数在任意初值下收敛于零,即闭环系统全局渐近稳定来给出高温气冷堆功率控制方法的形式和参数取值范围,并进而给出便于工程实现的离散化形式;
其中,过程指标函数是指刻画系统过渡过程品质的函数,本发明采用指标函数
和
来描述过渡过程品质,其中
en=ls-nr (3)
eT=Tout,r-Tout (4)
t>0表示时间,且σ1、σ2和θ皆为给定正常数。这里指标函数H1的物理意义是核功率测量值对其参考值的趋近程度,若H1随时间的推移而越来越小,则核功率的测量值渐近收敛于其参考值。指标函数H2描述了热氦温度测量值相对其参考值的趋近程度,若H2随时间的推移而越来越小,则堆芯出口氦气温度的测量值渐近收敛于其参考值。表达式(1)和(2)中的正常数σ1、σ2和θ也具有明显的物理意义,其中σ1表示核功率的偏差积分在指标函数H1中的比重越大,θ表示热氦温度偏差在指标函数H2中的比重,而σ2越大则表明误差项en+θeT的积分在指标函数H2中的比重越大。
本发明中可以通过设定参数σ1、σ2和θ的值来获取需要的过渡过程,通过增大σ1和σ2的值可以降低稳态误差,而通过增大θ的值可以提升堆芯出口氦气温度的控制性能。
在选定过渡过程指标函数H1的基础上,可设计得到图1中的核功率控制器u1(·)和热氦温度控制器u2(·)的表达式,且这两个控制器在频域内的表达式分别为
和
其中
KP,1=ΓP,1+ΓP,2 (7)
KD,1=ΓD,1+ΓD,2 (8)
KP,2=θΓP,2 (9)
KD,2=θΓD,2 (10)
Fi(i=1,2)为预设的正常数,且Gr为控制棒微分价值。参数Fi称为反馈系数其物理意义是表明了指标函数Hi的收敛速度,即若Fi越大则指标函数Hi收敛至0的速度也就越快,这里i=1,2。
由式(7)至(12)可知,核功率控制器(5)和堆芯出口热氦温度控制器(6)的参数与选定的过渡过程指标函数(1)和(2)之间具有明确的关系,即具备鲜明的物理意义。由于式(5)中包含纯微分环节,在具有中高频测量噪声的情况下不方便实现,因而也可应用如下的带有一阶惯性环节的核功率控制器,
其中κ为给定正常数,且其推荐取值范围是
0.001≤κ≤0.01 (14)
控制方法的形式设计给出了核功率控制器和热氦温度控制器的表达式及其参数与过渡过程指标函数之间的关系,经过计算只有当控制器参数在特定范围内时,闭环系统才是全局渐近稳定的。下面给出计算得到的控制器参数所须满足的取值范围。
核功率控制器的参数F1和σ1需满足的取值范围是
其中β为缓发中子份额,Λ是中子每代时间,ε1和ε2为给定常数,且推荐范围是
0<εi<1,i=1,2 (16)
堆芯出口热氦温度控制器的参数F2和σ2需满足的取值范围是
其中ε2为给定常数,且推荐范围是
0<ε3<1. (18)
基于式(7)-(12)以及(15)和(17)可得Γp,i和Γd,i(i=1,2)的取值范围是
在工程实际中,由核功率控制器(5)或(13)以及堆芯出口热氦温度控制器(6)构成的高温气冷堆功率控制方法通常实现为功率控制系统中的程序代码,这就需要给出上述连续控制器的离散化形式。
核功率控制器(5)的离散化形式为
其中Ts,n为核功率控制回路的采样周期。
堆芯出口热氦温度控制器(6)的离散化形式为
其中Ts,T为热氦温度控制的回路采样周期。
带有一阶惯性环节的核功率控制器(13)的离散化形式为
图2是按照本发明一种实施方式的高温气冷堆功率控制方法设计的流程图,本实施方式的方法包括以下步骤:
S1:设置核功率控制器参数和堆芯出口热氦温度控制器参数;
S2:所述堆芯出口热氦温度控制器根据堆芯出口热氦温度测量值与所述堆芯出口热氦温度参考值之间的偏差,通过计算获得相对核功率参考值的修正;所述核功率控制器根据由堆芯相对中子通量密度测量值和所述堆芯相对中子通量密度参考值之间的偏差通过计算给出驱动控制棒运动的棒速信号,以实现对反应堆功率的控制,所述堆芯相对中子通量密度参考值为所述相对核功率参考值和所述修正值相加口。
优选地,所述核功率控制器的离散化形式为:
其中,en=ls-nr,KP,1=ΓP,1+ΓP,2,KD,1=ΓD,1+ΓD,2,
σ1、σ2、θ、F1、以及F2均为预设的正常数,Ts,n为核功率控制回路的采样周期,Gr为控制棒微分价值,ls为相对堆芯中子通量密度参考值,nr为相对堆芯中子通量密度测量值,en为相对堆芯中子通量密度测量值和所述相对堆芯中子通量密度参考值之间的偏差,en(k)为当前采样时刻所对应的相对堆芯中子通量密度测量值和所述相对堆芯中子通量密度参考值之间的偏差,u1(k)为当前所述核功率控制函数的输出量;
所述堆芯出口热氦温度控制器的离散形式为:
其中,eT=Tout,r-Tout,KP,2=θΓP,2,KD,2=θΓD,2,
Ts,T为堆芯出口热氦温度控制回路的采样周期,Tout为堆芯出口热氦温度测量值,Tout,r为堆芯出口热氦温度参考值,u2(k)为当前相对核功率参考值的修正,eT为堆芯出口热氦温度测量值与所述堆芯出口热氦温度参考值之间的偏差,eT(k)为当前采样时刻所对应的堆芯出口热氦温度测量值与所述堆芯出口热氦温度参考值之间的偏差。
由于式(5)中包含纯微分环节,在具有中高频测量噪声的情况下不方便实现,因而也可应用如下的带有一阶惯性环节的核功率控制器,优选地,所述核功率控制器的离散化形式为:
其中,en=ls-nr,KP,1=ΓP,1+ΓP,2,KD,1=ΓD,1+ΓD,2,
σ1、σ2、F1以及F2均为预设的正常数,κ为预设的正常数、且取值范围为0.001≤κ≤0.01,Ts,n为核功率控制回路的采样周期,Gr为控制棒微分价值,ls为堆芯相对中子通量密度参考值,nr为堆芯相对中子通量密度测量值,en为堆芯相对中子通量密度测量值和所述堆芯相对中子通量密度参考值之间的偏差,en(k)为当前采样时刻所对应的堆芯相对中子通量密度测量值和所述堆芯相对中子通量密度参考值之间的偏差,u1(k)为当前所述核功率控制函数的输出量;
所述堆芯出口热氦温度控制器的离散化形式为:
其中,eT=Tout,r-Tout,KP,2=θΓP,2,KD,2=θΓD,2,
θ为预设的正常数,Ts,T为堆芯出口热氦温度控制回路的采样周期,Tout为堆芯出口热氦温度测量值,Tout,r为堆芯出口热氦温度参考值,u2(k)为当前相对核功率参考值的修正,eT为堆芯出口热氦温度测量值与所述堆芯出口热氦温度参考值之间的偏差,eT(k)为当前采样时刻所对应的堆芯出口热氦温度测量值与所述堆芯出口热氦温度参考值之间的偏差。
优选地,控制器参数σ1、σ2、F1以及F2的取值范围分别为:
其中,Λ为中子每代时间,β为缓发中子份额,vm为控制棒运动的最大速率,ε1、ε2和ε3均为预设的常数且取值范围为大于0小于1。
本发明提出的高温气冷堆功率控制方法解决了目前高温气冷堆功率控制方法所存在的形式和性能上的矛盾,其根本出发点就是通过严格的计算给出组成功率控制方法的核功率控制器和热氦温度控制器的具体形式和参数取值范围,该功率控制方法的主要特点为:
1、控制方法形式简洁,工程实现方便
由式(5)、(6)和(13)可知,组成功率控制方法的核功率控制器以及热氦温度控制器具有简洁的形式。此外,由式(20)、(21)和(22)可知,控制器的离散化形式也不复杂。因而,该功率控制方法便于工程实现。
2、控制方法参数直接与过渡过程指标函数对应,物理意义明确
通过选定指标函数(1)和(2)可以刻画所需的过渡过程,进而由于指标函数又决定了控制器(5)和(6)的形式和参数,因而控制方法能够保证过渡过程具有良好的品质。此外,由于控制器参数与指标函数参数是严格对应的,从而控制器参数的物理意义明确。
3、当控制方法参数在给定范围内时,控制方法能够保证闭环系统全局渐近稳定
目前存在的高温气冷堆功率控制方法仅能在给定工作点附近保证系统闭环的稳定性,而本发明提出的功率控制方法的参数在(15)和(17)给出的范围内时,就可以保证闭环系统的全局渐近稳定性,这将对高温气冷堆安全、稳定和高效的运行具有极为重要的意义。
实施例1
本发明所提出的高温气冷堆功率控制方法中的参数整定的具体实施方式可参照图3给出的流程图,具体步骤为:第一步,依据过渡过程指标函数的物理意义选择其参数σ1、σ2和θ的取值;第二步,给出反应堆物理参数,即缓发中子份额β、中子每代时间Λ和控制棒微分价值Gr的取值;第三步,根据式(15)和(17)中的第一个不等式计算指标函数参数σ1和σ2的取值范围;第四步,判断所选择的σ1、σ2和θ是否满足指标函数参数的取值范围,如果不满足则要重新选取σ1、σ2和θ,反之则继续设计流程;第五步,选择反馈系数F1和F2的取值;第六,依据式(7)至(12)计算控制器参数ΓP,1、ΓP,2、ΓD,1和ΓD,2的值;第七步,按照(16)和(18)给定正常数εi(i=1,2,3)和最大棒速率vm;第八步,根据(19)计算控制函数的参数ΓP,1、ΓP,2、ΓD,1和ΓD,2应满足的取值范围;第九步,判断控制器参数计算值是否在取值范围内,如果在则继续设计,反之则返回第五步;第十步,计算控制器参数KP,1、KP,2、KD,1和KD,2的值;第十一步,分别选定采样周期Ts,n和Ts,T,并给出功率控制方法的离散化形式,从而完成功率控制方法的设计和参数选择。所设计的核功率控制器的控制函数和热氦温度控制器最终固化为分布式控制系统中的程序代码。
下面给出一个高温气冷堆功率控制方法参数整定的具体实例:
第一步、选择过渡过程指标函数的参数满足σ1=0.2,σ2=0.1且θ=0.1。
第二步、给定某高温气冷堆物理参数:缓发中子份额β=0.005、中子每代时间Λ=0.001秒以及控制棒微分价值Gr=0.001。
第三步、根据式(15)和(17)中的第一个不等式可σ1和σ2的取值范围是
0<σi<10,i=1,2
第四步、根据第一步和第三步的结果可判断出过渡过程指标函数的参数值选择符合要求,设计过程继续。
第五步、选择设反馈系数F1=0.004,F2=0.002。
第六步、计算控制器参数ΓP,1、ΓP,2、ΓD,1和ΓD,2的取值为
第七步、给定正常数ε1=0.05,ε2=0.02,ε3=0.01和最大棒速率vm=0.05m/s。
第八步、根据(19)计算控制器参数ΓP,1、ΓP,2、ΓD,1和ΓD,2应满足的取值范围是
第九步、由第六步和第八步的结果可知ΓP,1、ΓP,2、ΓD,1和ΓD,2满足给定范围要求,设计继续进行。
第十步、依据式(7)至(12)计算控制器参数KP,1、KP,2、KD,1和KD,2的值分别为
第十一步、选择核功率和热氦温度回路的采样周期分别为Ts,n=0.1秒和Ts,T=0.5秒,则依据式(20)和(21)可计算离散化核功率控制器和热氦温度控制器的表达式分别为:
u1(k)=u1(k-1)+61en(k)-121en(k-1)+60en(k-2)
和
u2(k)=1.923u2(k-1)-0.923u2(k-2)+0.0321eT(k)
-0.063eT(k-1)+0.031eT(k-2).
这就完成了一个高温气冷堆功率控制方法的形式设计和参数整定流程,开发人员可根据最后得到的离散化控制方法迭代公式编写功率控制系统相关控制设备中的可执行代码。若采用带有惯性环节的核功率控制器(13),则也可按照上述步骤完成设计和参数整定。
本发明还公开了一种高温气冷堆功率控制系统,包括:
参数设置模块,用于设置核功率控制器的参数和堆芯出口热氦温度控制器的参数;
反应堆功率控制模块,用于所述堆芯出口热氦温度控制器根据堆芯出口热氦温度测量值与所述堆芯出口热氦温度参考值之间的偏差,通过计算获得相对核功率参考值的修正;所述核功率控制器根据由堆芯相对中子通量密度测量值和所述堆芯相对中子通量密度参考值之间的偏差通过计算给出驱动控制棒运动的棒速信号,以实现对反应堆功率的控制,所述堆芯相对中子通量密度参考值为所述相对核功率参考值和所述修正值相加。
以上实施方式仅用于说明本发明,而并非对本发明的限制,有关技术领域的普通技术人员,在不脱离本发明的精神和范围的情况下,还可以做出各种变化和变型,因此所有等同的技术方案也属于本发明的范畴,本发明的专利保护范围应由权利要求限定。
Claims (5)
1.一种高温气冷堆功率控制方法,其特征在于,包括以下步骤:
S1:设置核功率控制器参数和堆芯出口热氦温度控制器参数;
S2:所述堆芯出口热氦温度控制器根据堆芯出口热氦温度测量值与所述堆芯出口热氦温度参考值之间的偏差,通过计算获得相对核功率参考值的修正;所述核功率控制器根据由堆芯相对中子通量密度测量值和所述堆芯相对中子通量密度参考值之间的偏差通过计算给出驱动控制棒运动的棒速信号,以实现对反应堆功率的控制,所述堆芯相对中子通量密度参考值为所述相对核功率参考值和所述修正值相力。
2.如权利要求1所述的方法,其特征在于,所述核功率控制器的离散化形式为:
其中,en=ls-nr,KP,1=ΓP,1+ΓP,2,KD,1=ΓD,1+ΓD,2,
σ1、σ2、θ、F1以及F2均为预设的正常数,Ts,n为核功率控制回路的采样周期,Gr为控制棒微分价值,ls为堆芯相对中子通量密度的参考值,nr为堆芯相对中子通量密度测量值,en为相对堆芯中子通量密度测量值和所述堆芯相对中子通量密度参考值之间的偏差,en(k)为当前采样时刻所应对的堆芯相对中子通量密度测量值和所述堆芯相对中子通量密度参考值之间的偏差,u1(k)为当前所述核功率控制器的输出量;
所述堆芯出口热氦温度控制器的离散化形式为:
其中,eT=Tout,r-Tout,KP,2=θΓP,2,KD,2=θΓD,2,
Ts,T为堆芯出口热氦温度控制回路的采样周期,Tout为堆芯出口热氦温度测量值,Tout,r为堆芯出口热氦温度参考值,u2(k)为当前相对核功率参考值的修正,eT为堆芯出口热氦温度测量值与所述堆芯出口热氦温度参考值之间的偏差,eT(k)为当前采样时刻所对应的堆芯出口热氦温度测量值与所述堆芯出口热氦温度参考值之间的偏差。
3.如权利要求1所述的方法,其特征在于,所述核功率控制器的离散化形式为:
其中,en=ls-nr,KP,1=ΓP,1+ΓP,2,KD,1=ΓD,1+ΓD,2,
σ1、σ2、F1以及F2均为预设的正常数,κ为预设的正常数、且取值范围为0.001≤κ≤0.01,Ts,n为核功率控制回路的采样周期,Gr为控制棒微分价值,ls为堆芯相对中子通量密度参考值,nr为堆芯相对中子通量密度测量值,en为堆芯相对中子通量密度测量值和所述堆芯相对中子通量密度参考值之间的偏差,en(k)为当前采样时刻所对应的堆芯相对中子通量密度测量值和所述堆芯相对中子通量密度参考值之间的偏差,u1(k)为当前所述核功率控制器的输出量;
所述堆芯出口热氦温度控制器的离散化形式为:
其中,eT=Tout,r-Tout,KP,2=θΓP,2,KD,2=θΓD,2,
θ为预设的正常数,Ts,T为堆芯出口热氦温度控制回路的采样周期,Tout为堆芯出口热氦温度测量值,Tout,r为堆芯出口热氦温度参考值,u2(k)为当前相对核功率参考值的修正,eT为堆芯出口热氦温度测量值与所述堆芯出口热氦温度参考值之间的偏差,eT(k)为当前采样时刻所对应的堆芯出口热氦温度测量值与所述堆芯出口热氦温度参考值之间的偏差。
4.如权利要求2或3所述的方法,其特征在于,控制器参数σ1、σ2、F1、以及F2的取值范围分别为:
其中,Λ为中子每代时间,β为缓发中子份额,vm为控制棒运动的最大速率,ε1,ε2和ε3均为预设的常数、且取值范围为大于0小于1。
5.一种高温气冷堆功率控制系统,其特征在于,包括:
参数设置模块,用于设置核功率控制器参数和堆芯出口热氦温度控制器参数;
功率控制模块,用于所述堆芯出口热氦温度控制器根据堆芯出口热氦温度的测量值与所述堆芯出口热氦温度的参考值之间的偏差,通过计算获得相对核功率参考值的修正;所述核功率控制器根据由堆芯相对中子通量密度的测量值和所述堆芯相对中子通量密度的参考值之间的偏差通过计算给出驱动控制棒运动的棒速信号,以实现对反应堆功率的控制,所述堆芯相对中子通量密度参考值为所述相对核功率参考值和所述修正值相加。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201110328271.6A CN102411997B (zh) | 2011-10-25 | 2011-10-25 | 高温气冷堆功率控制方法及系统 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201110328271.6A CN102411997B (zh) | 2011-10-25 | 2011-10-25 | 高温气冷堆功率控制方法及系统 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN102411997A true CN102411997A (zh) | 2012-04-11 |
CN102411997B CN102411997B (zh) | 2014-04-02 |
Family
ID=45914032
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201110328271.6A Active CN102411997B (zh) | 2011-10-25 | 2011-10-25 | 高温气冷堆功率控制方法及系统 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN102411997B (zh) |
Cited By (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103218512A (zh) * | 2013-03-14 | 2013-07-24 | 西安交通大学 | 一种获取核燃料组件中子角通量密度的方法 |
CN103400621A (zh) * | 2013-07-29 | 2013-11-20 | 国核电力规划设计研究院 | 一种高温气冷堆核电站的启停堆系统的控制方法和装置 |
CN106297925A (zh) * | 2016-07-28 | 2017-01-04 | 董哲 | 一种多模块式高温气冷堆核电站的给水方案、动态模型与调控方法 |
CN111192696A (zh) * | 2020-01-06 | 2020-05-22 | 中国原子能科学研究院 | 用于核反应堆的功率控制装置及功率控制方法 |
CN111564226A (zh) * | 2020-04-03 | 2020-08-21 | 中国原子能科学研究院 | 一种基于模糊控制的研究堆功率自动调节方法 |
CN112382418A (zh) * | 2020-11-20 | 2021-02-19 | 西安热工研究院有限公司 | 带有增量式调节功能的高温气冷堆氦气流量控制系统及方法 |
CN114496319A (zh) * | 2021-12-08 | 2022-05-13 | 华能核能技术研究院有限公司 | 高温气冷堆机组的运行状态的确定方法及装置 |
CN116994787A (zh) * | 2023-07-28 | 2023-11-03 | 华能核能技术研究院有限公司 | 一种高温气冷堆核电厂核功率控制的方法及系统 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5158738A (en) * | 1988-11-14 | 1992-10-27 | Framatome | Method of controlling a pressurized water nuclear reactor |
US20080069288A1 (en) * | 2006-05-22 | 2008-03-20 | Areva Np | Method for regulation of operational parameters of the core of a pressurised water nuclear reactor |
CN101669176A (zh) * | 2007-03-19 | 2010-03-10 | 阿海珐核能公司 | 确定核反应堆堆芯三维功率分布的方法 |
CN101755308A (zh) * | 2007-06-08 | 2010-06-23 | 阿海珐核能公司 | 确定关于核反应堆堆芯中功率分布的不确定性分量的方法 |
US20110002432A1 (en) * | 2009-07-01 | 2011-01-06 | Westinghouse Electric Company Llc | Incore instrument core performance verification method |
-
2011
- 2011-10-25 CN CN201110328271.6A patent/CN102411997B/zh active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5158738A (en) * | 1988-11-14 | 1992-10-27 | Framatome | Method of controlling a pressurized water nuclear reactor |
US20080069288A1 (en) * | 2006-05-22 | 2008-03-20 | Areva Np | Method for regulation of operational parameters of the core of a pressurised water nuclear reactor |
CN101669176A (zh) * | 2007-03-19 | 2010-03-10 | 阿海珐核能公司 | 确定核反应堆堆芯三维功率分布的方法 |
CN101755308A (zh) * | 2007-06-08 | 2010-06-23 | 阿海珐核能公司 | 确定关于核反应堆堆芯中功率分布的不确定性分量的方法 |
US20110002432A1 (en) * | 2009-07-01 | 2011-01-06 | Westinghouse Electric Company Llc | Incore instrument core performance verification method |
Non-Patent Citations (3)
Title |
---|
冯俊婷等: "核反应堆功率调节系统控制特性研究", 《原子能科学技术》 * |
曾海: "中国先进研究堆功率调节系统仿真分析", 《原子能科学技术》 * |
李海鹏等: "模块式高温气冷堆核电站的运行特性研究", 《原子能科学技术》 * |
Cited By (13)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103218512A (zh) * | 2013-03-14 | 2013-07-24 | 西安交通大学 | 一种获取核燃料组件中子角通量密度的方法 |
CN103400621A (zh) * | 2013-07-29 | 2013-11-20 | 国核电力规划设计研究院 | 一种高温气冷堆核电站的启停堆系统的控制方法和装置 |
CN103400621B (zh) * | 2013-07-29 | 2016-02-24 | 国核电力规划设计研究院 | 一种高温气冷堆核电站的启停堆系统的控制方法和装置 |
CN106297925A (zh) * | 2016-07-28 | 2017-01-04 | 董哲 | 一种多模块式高温气冷堆核电站的给水方案、动态模型与调控方法 |
CN111192696A (zh) * | 2020-01-06 | 2020-05-22 | 中国原子能科学研究院 | 用于核反应堆的功率控制装置及功率控制方法 |
CN111564226B (zh) * | 2020-04-03 | 2022-05-13 | 中国原子能科学研究院 | 一种基于模糊控制的研究堆功率自动调节方法 |
CN111564226A (zh) * | 2020-04-03 | 2020-08-21 | 中国原子能科学研究院 | 一种基于模糊控制的研究堆功率自动调节方法 |
CN112382418A (zh) * | 2020-11-20 | 2021-02-19 | 西安热工研究院有限公司 | 带有增量式调节功能的高温气冷堆氦气流量控制系统及方法 |
CN112382418B (zh) * | 2020-11-20 | 2021-08-31 | 西安热工研究院有限公司 | 带有增量式调节功能的高温气冷堆氦气流量控制系统及方法 |
WO2022105357A1 (zh) * | 2020-11-20 | 2022-05-27 | 西安热工研究院有限公司 | 带有增量式调节功能的高温气冷堆氦气流量控制系统及方法 |
CN114496319A (zh) * | 2021-12-08 | 2022-05-13 | 华能核能技术研究院有限公司 | 高温气冷堆机组的运行状态的确定方法及装置 |
CN116994787A (zh) * | 2023-07-28 | 2023-11-03 | 华能核能技术研究院有限公司 | 一种高温气冷堆核电厂核功率控制的方法及系统 |
CN116994787B (zh) * | 2023-07-28 | 2024-05-24 | 华能核能技术研究院有限公司 | 一种高温气冷堆核电厂核功率控制的方法及系统 |
Also Published As
Publication number | Publication date |
---|---|
CN102411997B (zh) | 2014-04-02 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN102411997B (zh) | 高温气冷堆功率控制方法及系统 | |
Deng et al. | Model predictive control of central chiller plant with thermal energy storage via dynamic programming and mixed-integer linear programming | |
Dong et al. | Multi-layer perception based model predictive control for the thermal power of nuclear superheated-steam supply systems | |
CN104503228B (zh) | 一种水轮机调速器功率模式下一次调频稳定域的确定方法 | |
CN109212974A (zh) | 区间时变时滞系统的鲁棒模糊预测容错控制方法 | |
Banerjee et al. | An interval approach for robust control of a large PHWR with PID controllers | |
CN104950670B (zh) | 一种连续搅拌釜式反应器的一体化多模型控制方法 | |
Wang et al. | Modification and analysis of load follow control without boron adjustment for CPR1000 | |
Sun et al. | Application of model free active disturbance rejection controller in nuclear reactor power control | |
CN105353811A (zh) | 一种用于液氮低温装置的智能温度控制系统及其控制方法 | |
CN105387449A (zh) | 一种在锅炉蒸汽温度控制中使用二阶微分的控制方法 | |
CN206075835U (zh) | 一种电加热模拟核释热的系统 | |
Dong | Nonlinear coordinated control for MHTGR-based nuclear steam supply systems | |
WO2023184899A1 (zh) | 核反应堆的自主控制方法、装置和计算机设备 | |
Dong et al. | Power-level control of nuclear reactors based on feedback dissipation and backstepping | |
CN113653589B (zh) | 一种变参双幂趋近律的水轮机组微分预测滑模控制方法 | |
CN205210740U (zh) | 一种用于液氮低温装置的智能温度控制系统 | |
CN112992292B (zh) | 一种压水堆一回路冷却剂中硼10丰度的实时预测方法 | |
CN104035336A (zh) | Mimo非最小相位cstr反应器的非线性控制器设计方法 | |
Lee et al. | Robust controller design of nuclear power reactor by parametric method | |
Cai et al. | Experimental study on gas temperature control for a high-speed heat-airflow wind tunnel | |
Ashaari et al. | State space modeling of reactor core in a pressurized water reactor | |
CN112186752B (zh) | 一种单目标多时段精准调节方法 | |
Zhang et al. | Switching controller for enthalpy difference laboratory temperature based on fuzzy control and PID control | |
Dong | Power-level control for MHTGRs with time-delay in helium temperature measurement |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
TR01 | Transfer of patent right | ||
TR01 | Transfer of patent right |
Effective date of registration: 20210120 Address after: 100193 building 26, Zhongguancun Software Park, 8 Dongbeiwang West Road, Haidian District, Beijing Patentee after: CHINERGY Co.,Ltd. Address before: 100084 mailbox, 100084-82 Tsinghua Yuan, Beijing, Haidian District, Beijing Patentee before: TSINGHUA University |