CN105321588B - 一种核电厂超温超功率回路保护通道定期试验验证准则的整定方法 - Google Patents
一种核电厂超温超功率回路保护通道定期试验验证准则的整定方法 Download PDFInfo
- Publication number
- CN105321588B CN105321588B CN201410380428.3A CN201410380428A CN105321588B CN 105321588 B CN105321588 B CN 105321588B CN 201410380428 A CN201410380428 A CN 201410380428A CN 105321588 B CN105321588 B CN 105321588B
- Authority
- CN
- China
- Prior art keywords
- time
- function
- mrow
- voltage
- domain
- 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
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
Landscapes
- Monitoring And Testing Of Nuclear Reactors (AREA)
Abstract
本发明提供了一种核电厂超温超功率回路保护通道定期试验验证准则的整定方法,可根据事故工况,模拟从现场采集到的事故工况信号,经过计算得到在所述事故工况下阈值模块的输入电压在时间域的函数,当所述输入电压达到阈值电压时,将触发所述事故工况的阈值模块动作,再根据整个保护通道的系统偏差,得到触发所述事故工况动作的时间范围:Tmax,Tmin和电压范围Vmax,Vmin,为核电厂超温超功率回路保护通道定期试验提供验证准则。
Description
技术领域
本发明属于核电站安全检查技术领域,特别涉及一种核电站超温超功率回路保护通道定期试验验证准则的整定方法。
背景技术
仪表控制系统是核电站的神经系统,它直接关系到核电站安全、可靠和稳定的运行,KRG系统是集中控制模拟量系统,是核电站仪表控制系统的关键系统,它的功能定义了从现场传感器输出,到反应堆保护系统(RPR系统)的输入信号,信号在KRG系统内进行处理(包括开方、线性化、滤波、微分、超前-滞后等等)之后,执行阈值比较,当信号超过设定的限值后,就会触发保护信号并发送到RPR系统执行保护逻辑。
在核电厂运行期间,要保证RPR系统的可用性以及有效性,并根据安全法规和国家标准的要求,必须有有效的手段在核电站运行期间可以对KRG系统的各个保护通道进行定期的验证,以确保可以及时、准确地发现系统内部潜在的故障,因此核电站有一套专门的保护系统定期试验装置。
对于核电厂保护系统定期试验装置项目,超温超功率保护回路是保护系统中最为复杂的回路,其涉及到冷管段温度、热管段温度、主泵转速、稳压器压力、中子通量等五个状态参数以及多个动态模块的计算。因此,其相应的定期试验的时间和电压验证准则整定算法也是最为复杂的,其定期试验基本原理为:向保护系统注入模拟事故工况信号,当注入的信号变化到阈值模块的设定值后,就会触发阈值模块动作,而此时,采集动作信号,如果动作电压和时间都在试验准则(Vmax和Vmin,Tmax和Tmin)的范围内,则试验成功,表示该保护通道工作正常。而如果得到的电压和时间落在准则范围外,则试验失败,维护人员需要检查整个保护通道的板件来确定问题所在,进而调整或更换故障板件直到试验可以顺利通过。
当采集冷热管段温度的探头寿命到期或电厂大修被更换后,其属性的微小改变都会使整个超温超功率保护回路的特性发生变化,进而需要操作员手动调整实现冷热管段的平均温度和温度偏差的两个板件的参数。当设备参数有变化的时候,需要重新整定试验准则,即重新确定Vmax、Vmin、Tmax和Tmin的值。
目前保护通道定期试验验证准则的整定方法只能应用于简单通道的保护逻辑,而不适用于超温超功率这样复杂的回路计算。
发明内容
保护通道的超温超功率回路是保护系统中最为复杂的回路,其功能为一回路的三个环路的冷热管段温差超过设定限值后,将触发停堆保护信号,该保护通道主要涉及到冷管段温度、热管段温度、主泵转速、稳压器压力、中子通量五个状态参数的处理和计算,当现场板卡的属性发生改变后,例如:当现场冷热管段温度转换板卡CT更换或转换曲线调整,也就是说对于相同温度转换后的电压是不同的,进而冷\热管段的平均温度和温度偏差的两个模块的参数也需要随之进行调整,因此,当保护通道的参数发生改变后,其阈值触发模块的时间和电压验证准则(Vmax和Vmin,Tmax和Tmin)需要重新进行计算。
本发明所述的超温超功率回路的定期试验是将该保护功能分为超温,超功率以及平均温度三个功能来分别进行验证,而在超温超功率的三个功能验证中又分别模拟了四种不同的事故工况进行验证,这四种事故工况分别为冷热管段温度升高和稳压器压力升高,中子通量升高,主泵转速降低以及冷热管段温度升高和中子通量降低。每一种工况分别注入不同的试验信号,每组信号代表着不同的事故瞬态。因此,每种工况都有着不同时间和电压动作准则。
本发明提供一种定期试验验证准则整定算法,可以解决超温超功率这样复杂回路定期试验的时间和电压验证准则的计算。
本发明提供的定期试验验证准则整定的主要原理为将输入信号进行拉普拉斯变换后,在S域下,根据信号传输路径进行逐步地计算,得出阈值模块输入电压的S域函数,再进行拉普拉斯反变换,得到输入电压与时间t的对应关系,最后带入阈值模块设定的电压动作值,就可以解出动作在这一瞬态的标准时间,为了使得超温超功率试验装置能更加合理的模拟事故工况,再将整个保护通道的系统偏差带入,以此计算出试验项的时间准则范围Tmax和Tmin以及电压准则范围Vmax和Vmin。
本发明通过详细介绍冷热管段温度升高和稳压器压力升高的事故工况下试验电压准则和时间准则的整定方法,其它三种事故工况中子通量升高,主泵转速降低以及冷热管段温度升高和中子通量降低的试验验证准则的推导方法也与上述方法一致,只是初始端的注入信号不同,因此,简化了定期试验验证准则的整定方法,提高了试验验证准则的推导效率。
本发明解决其技术问题所采用的技术方案是一种核电厂超温超功率回路保护通道定期试验验证准则的整定方法,将模拟事故工况的信号,通过计算得到一个输入电压V随时间t变化的函数,当所述输入电压V为阈值电压V0时,将触发所述事故工况,再根据整个保护通道的系统偏差,得到触发所述事故工况动作的验证准则,所述模拟事故工况的信号为输入电压在时间域t上的函数;所述的验证准则为时间验证准则:Tmax、Tmin和电压验证准则:Vmax,Vmin。
其中,所述的模拟事故工况的信号为:
冷管段温度在时间域的函数Vc=0.111429t+2.565714,
热管段温度在时间域的函数Vh=0.037143t+3.605714,
稳压器压力在时间域的函数Vp=0.228571t+3.514286,
主泵转速在时间域的函数Vw,为常数4.3,
中子通量在时间域的函数Vf1为常数1;
所述输入电压V随时间t变化的函数为:
所述计算过程为:
(1)首先,所述冷管段温度在时间域的函数Vc=0.111429t+2.565714经过一阶惯性环节415FI的滤波:
(1.1)由于所述一阶惯性环节415FI的滤波对定值无影响,因此假设冷管段注入信号的初始值为0,即不考虑2.565714,将Vc=0.111429t+2.565714经过拉氏变换后,得到:
(1.2)所述一阶惯性环节415FI的传导函数G415FI:
输入信号为经过拉式变换的冷管段注入信号
由一阶惯性环节415FI的输出信号X1:X1=Vc*G415FI,
得到:然后经过拉式反变换,并将所述初始值2.565714带入方程,得到:X1=0.111429(t+e-t-1)+2.565714;
(2)然后,所述热管段温度在时间域的函数Vh=0.037143t+3.605714经过一阶惯性环节416FI的滤波:
(2.1)由于所述一阶惯性环节416FI对定值无影响,因此假设热管段注入信号的初始值为0,即不考虑3.605714,将Vh=0.037143t+3.605714经过拉氏变换后,得到:
(2.2)所述一阶惯性环节416FI的传导函数G416FI:
输入信号为经过拉式变换的热管段注入信号
由一阶惯性环节416FI的输出信号X2:X2=Vh*G416FI,
得到:然后经过拉式反变换,并将所述初始值3.605714带入方程,得到:X2=0.037143(t+e-t-1)+3.605714;
(3)再将经过所述一阶惯性环节415FI和所述一阶惯性环节416FI滤波后的所述冷管段温度在时间域的函数和所述热管段温度在时间域的函数带入平均温度402ZO的传导公式:G402ZO=a*(X1+X2)+b,并经过计算,得到平均温度
V402ZO=a*0.148572*(t+e-t-1)+a*6.171428+b,其中,0.57≤a≤0.59、-0.517≤b≤-0.465;
(4)由404MT超前滞后校正得到温度响应处理后的平均温度函数,
(4.1)所述的404MT超前滞后的传导函数G404MT:其中T1为超前参数,取值为30,T2为滞后参数,取值为4,因此得到:
(4.2)由于所述超前滞后环节404MT对定值无影响,因此假设平均温度初始值为0,即不考虑-a*0.148572+a*6.171428+b,经过拉氏变换后,得到
(4.3)由于所述超前滞后环节404MT的输入为:
V402ZO=a*0.148572*(t+e-t-1)+a*6.171428+b,
传导函数为:因此可得出所述超前滞后404MT的输出值为:
(4.4)在将式进行拉式反变换,并带入初始值-a*0.148572+a*6.171428+b后,得出:
(5)将经过所述一阶惯性环节415FI和所述一阶惯性环节416FI滤波后的所述冷管段温度在时间域的函数和所述热管段温度在时间域的函数带入温度偏差403ZO的传导公式:G403ZO=c*(X1-X2)+d,并经过计算,得到平均温度
V403ZO=-c*0.074286*(t+e-t-1)+c*1.04+d,其中,-1.556≤c≤1.556、1.711≤d≤1.889;
(6)再将V404MT,稳压器信号Vp=0.228571t+3.514286,主泵转速信号Vw=4.3,中子通量Vf1=1,带入超温计算式:
V404ZO=6.924-1.672*V404MT+0.810*VP+0.267*Vw-Vf1中,得到
(7)最后根据所述工况的输入电压:V=V404ZO-V403ZO,将所述步骤(5)和所述步骤(6)得到V404ZO、V403ZO带入,并经过计算,得到所述输入电压V随时间t变化的函数:
这里,所述事故工况为:冷/热管段温度升高和稳压器压力升高而导致超温的事故工况;
所述阈值为超温阈值444XU1,所述阈值电压V0=0V;
所述系统偏差为:所述超温超功率回路保护通道允许漂移边界MA,MA=1.667%,所述超温超功率回路保护通道中电压值范围为:1V~5V,优选4V,由此得到系统偏差VMA:VMA=4V*MA=4V*1.667%=0.067V,
根据所述输入电压V随时间t变化的函数:
得到:
由此可以得到式中所述时间t的取值范围动作时间上限T+和动作时间下限T-。
本发明所述的时间验证准则还包括触发所述事故工况阈值模块动作的响应时间Trelay和所述信号的采样时间Tsample因此,所述时间验证准则中Tmax:Tmax=T++Trelay,和Tmin:Tmin=T-+Trelay+Tsample;优选的,所述事故工况阈值模块动作响应时间Trelay=27.75ms,所述信号的采样时间Tsample=3.6ms。
本发明所述的电压验证准则,还包括所述阈值444XU1允许的漂移电压所述超温超功率回路保护通道的电压范围,优选1V~5V,再优选4V,所述漂移电压为:Vxu为所述超温超功率回路保护通道允许的漂移量,Vxu=0.41%,因此,得到:
所述的电压验证准则,还包括电压变化率D在时间域上的函数:
所述的电压变化率D在时间域上的函数,还包括触发所述事故工况后的动作时的响应时间Trelay和所述信号的采样时间Tsample;
因此所述的电压验证准则为:
优选的,本发明还包括与所述超温阈值模块互为备用的超温阈值模块444XU2,所述超温阈值模块444XU2与所述超温阈值模块444XU1实现相同的功能。
本发明所述事故工况还包括:
1)中子通量升高而导致平均温度高的事故工况,所述阈值为443XU1,动作电压为0V,所述模拟事故工况的信号为:
冷管段温度在时间域的函数Vc=0.005714t+2.565714,
热管段温度在时间域的函数Vh=0.022857t+4.005714,
稳压器压力在时间域的函数Vp=0.028571t+3.514286,
主泵转速在时间域的函数Vw,为常数4.3,
中子通量在时间域的函数为Vf1=0.025t+0.15;
2)主泵转速降低而导致超功率的事故工况,所述阈值为445XU1,动作电压为0V,所述模拟事故工况的信号为:
冷管段温度在时间域的函数为常数,Vc=2.44,
热管段温度在时间域的函数为常数,Vh=3.477714,
稳压器压力在时间域的函数为常数Vp=3.514286,
主泵转速在时间域的函数,Vw=-0.222222t+4.3,
中子通量在时间域的函数Vf1为常数3.35。
优选的,本发明还包括与所述平均温度阈值模块443XU1互为备用的平均温度阈值模块443XU2;与所述超功率阈值模块445XU1互为备用的超功率阈值模块445XU2。
本发明的有益效果是:
1)在核电站领域,提出了对最为复杂的超温超功率回路保护通道定期试验验证准则的整定方法;
2)将数个电压在时间域上的函数,合并成一个输入电压在时间域的函数,简化了确定事故工况的阈值范围的整定,同时根据模拟工况下,输入的电压的阈值得到触发事故工况时的标准时间;
3)由于在整定过程中加入了超温超功率回路保护通道中各个板件的偏差值,为触发事故工况的动作值设定了一个范围,使得超温超功率试验装置能更加合理的模拟事故工况;
4)同时,考虑了超温阈值动作的响应时间和定期试验装置的采样时间,使得触发所述事故工况动作的时间范围:Tmax,Tmin和电压范围Vmax,Vmin更加精确。
5)由于为超温超功率回路保护通道定义了整定方法,当系统更换板件,人工调整参数后,可以简单、快速的计算出整个超温超功率回路定期试验准则,为整个回路的试验提供精确的数据支持。
附图说明
下面结合附图对本发明所述的一种核电厂超温超功率回路保护通道定期试验验证准则的整定方法进行具体说明。
图1是本发明第一环路的超温超功率保护通道的模拟图;
图2是本发明一回路冷热管段温度升高和稳压器压力升高导致超温,从而触发阈值模块444XU1动作的流程图;
附图说明:033MT为冷管段温度的现场热电阻信号;030MT为热管段温度的现场热电阻信号;419CT和418CT模块是将热电阻信号转换为电压信号;415FI,416FI和401FI为一阶惯性环节模块;402ZO,403ZO,404ZO,405ZO以及406ZO为加法器运算模块;404MT为一阶超前滞后环节模块;401DR为微分环节模块;402MS和416AM为中子通量采集处理以及计算模块;403GD和404GD为函数发生模块;443XU1/2为平均温度阈值模块,444XU1/2超温阈值模块,445XU1/2为超功率阈值模块。
具体实施方式
岭澳I期核电厂中整个一回路是由三个环路组成,而这三个环路的超温超功率保护功能是完全相同的,如图1所示,其中,033MT和030MT分别表示冷管段和热管段温度的现场热电阻信号;419CT和418CT模块是将热电阻信号转换为电压信号;415FI,416FI和401FI为一阶惯性环节模块;402ZO,403ZO,404ZO,405ZO以及406ZO为加法器运算模块;404MT为一阶超前滞后环节模块;401DR为微分环节模块;402MS和416AM为中子通量采集处理以及计算模块;403GD和404GD为函数发生模块;而443XU1/2,444XU1/2以及445XU1/2分别为平均温度阈值模块,超温阈值模块,超功率阈值模块。KRG系统的每个保护通道都是上述的这些模拟量处理以及预算模块组。
1试验工况
以一回路冷热管段温度升高和稳压器压力升高导致超温,从而触发阈值模块(444XU1和444XU2)动作为例,推导出在该事故工况下试验验证准则的整定方法。
下面给出模拟该事故工况时试验的注入信号(时间域):
热管段温度值电压Vh在时间域的函数:Vh=0.037143t+3.605714;
冷管段温度值电压Vc在时间域的函数:Vc=0.111429t+2.565714;
稳压器压力值电压Vp在时间域的函数:Vp=0.228571t+3.514286;
主泵转速值电压Vw在时间域的函数:Vw=4.3;
中子通量值电压Vf1在时间域的函数:Vf1=1。
2通道模块参数
在超温超功率保护通道中,除了平均温度和温度偏差两个模块402ZO和403ZO的增益和偏置两个参数可调外,其它模块的参数设定都是固定的。具体介绍如下(拉氏变换后):
415FI模块的传递函数G415FI:
416FI模块的传递函数G416FI:
402ZO模块的传递函数G402ZO:G402ZO=a*(X1+X2)+b,a和b为402ZO模块设定值,a和b分别在0.57到0.59和-0.517到-0.465之间可调;
403ZO模块的传递函数G403ZO:G403ZO=c*(X1-X2)+d,c和d为403ZO模块设定值,c和d分别在-1.556到1.556和1.711到1.889之间可调;
其中,X1=Vc*G415FI,X2=Vh*G416FI。
404MT模块的传递函数G404MT:其中超前参数为30,滞后参数为4;
404ZO模块的传递函数G404ZO:G404ZO=6.924-1.672*V404MT+0.810*VP+0.267*VW-Vf1。
3试验验证准则算法
3.1时间准则算法
冷管段注入信号为:Vc=0.111429t+2.565724,
由于FI滤波模块(一阶惯性环节)对定值信号无作用,因此,假设冷管段注入信号的初始值为“0”,经过拉氏变换后,得到:
经过415FI模块后,
将上式进行分解,可以得出
将上式拉式反变换,并带入初始值后,得到:X1=0.111429(t+e-t-1)+2.565714;
而热管段注入信号为:Vh=0.037143t+3.605714,
同理,可以得出:X2=Vh*G416FI=0.037143(t+e-t-1)+3.605714;
将X1和X2带入平均温度和温度偏差的方程中,可得出:
402ZO模块的输出值V402ZO=a*0.148572*(t+e-t-1)+a*6.171428+b,
403ZO模块的输出值V403ZO=-c*0.074286*(t+e-t-1)+c*1.04+d;
而超前滞后模块404MT的输入值为V402ZO=a*0.148572*(t+e-t-1)+a*6.171428+b,超前滞后模块也对定值信号无作用,因此,可以继续假设其初始值为“0”,经过拉氏变换后,
并且以此可得出:404MT模块的输出值为V404MT,
将上式进行拉式反变换,并带入初始值后,得出:
而V404ZO=6.924-1.672*V404MT+0.810*VP+0.267*Vw-Vf1,将V404MT,Vp,Vw以及Vf1的输入值带入后,可以得出:
整个超温超功率保护通道所有板件的允许漂移边界MA为1.667%,由此可以得出,VMA=4V*MA=4V*1.667%=0.067V,其中通道计算处理的物理值为1~5V;
根据现场的设定值,444XU1阈值模块的动作电压设定值为0V,因此,可以得出:V404ZO-V403ZO=±0.067V,带入之前推导出的V404ZO和V403ZO时域公式得出下面的时域方程,其中a,b,c,d都是已知的设定值,可以计算出T+和T-;
试验准则除了受保护通道所有板件的允许漂移边界MA影响外,还需考虑XU阈值模块本身的动作响应时间(Trelay)以及定期试验装置的采样时间(Tsample)等因素对动作时间采集准确性的影响,就可以得出此保护功能定期试验的时间验证准则Tmax和Tmin,其中Tmax=T++Trelay,Tmin=T-+Trelay+Tsample,对于KRG系统的Trelay为27.75ms,对于试验装置的Tsample为3.6ms。
而444XU2的时间验证准则也可以通过相同算法计算得出,此处不再重复描述。
以上所推导出来的是超温超功率保护通道在一回路冷热管段温度升高和稳压器压力升高,从而触发停堆保护逻辑的工况下执行试验的时间准则整定方法。
3.2电压准则算法
根据上面时间准则所推导出的VMA以及T+和T-,可以计算出整个通道允许漂移边界的电压变化率D(V/T),
而KRG系统中XU阈值模块的允许漂移VXU=0.41%,因此,可以得出允许漂移的电压值
根据现场的设定值,444XU1阈值模块的动作电压设定值为0V,因此,其允许的动作电压范围为-0.016V~+0.016V;
再考虑XU阈值模块本身的动作响应时间,以及定期试验装置的采样时间等因素对电压采集准确性的影响,得出此保护功能定期试验的电压验证准则Vmax和Vmin,其中而Trelay和Tsample同上,分别为27.75ms和3.6ms。
而444XU2的电压验证准则也可以通过相同算法计算得出,此处也不再重复描述。
以上所推导出来的是超温超功率保护通道在一回路冷热管段温度升高和稳压器压力升高,从而触发停堆保护逻辑的工况下执行试验的电压准则整定方法。
4实施例数据
在定期试验装置研发应用的项目中,通过Matlab编制上面描述的整定方法,然后将其嵌入到Labview开发的就地参数计算软件模块中,可以实现当重新设定CT模块属性以及平均温度和温度偏差模块的参数后,自动计算出超温超功率保护回路定期试验时间和电压准则的功能并更新到数据库中,满足试验要求。
例如,当CT的属性变化后R286=207.06Ω(现场286度时对应的热电阻值),R323=220.2Ω(现场323度时对应的热电阻值),而402ZO的a设为-0.57,b设为-0.5,403ZO的c设为1.276,d设为1.729,则整定后444XU1的Vmax和Vmin分别为0.0165V和-0.0154V,Tmax和Tmin分别为7.0886s和6.3386s,而444XU2的Vmax和Vmin分别为0.0878V和0.0553V,Tmax和Tmin分别为6.9556s和5.7669s;当CT的属性变化后R286=207.06Ω(现场286度时对应的热电阻值),R323=220.2Ω(现场323度时对应的热电阻值),而402ZO的a设为-0.583,b设为-0.5,403ZO的c设为1.276,d设为1.729,则整定后444XU1的Vmax和Vmin分别为0.0108V和-0.0221V,Tmax和Tmin分别为5.7757s和5.0590s,而444XU2的Vmax和Vmin分别为0.0850V和0.0519V,Tmax和Tmin分别为5.3461s和4.7288s;当CT的属性变化后R286=207.06Ω(现场286度时对应的热电阻值),R323=220.2Ω(现场323度时对应的热电阻值),而402ZO的a设为-0.59,b设为-0.5,403ZO的c设为1.276,d设为1.729,则整定后444XU1的Vmax和Vmin分别为0.0095V和-0.0236V,Tmax和Tmin分别为5.1780s和4.6085s,而444XU2的Vmax和Vmin分别为0.0838V和0.0505V,Tmax和Tmin分别为4.8420s和4.3306s。
以上所述仅是本发明的较佳实施例而已,并非对本发明作任何形式上的限制,虽然本发明已以较佳实施例揭露如上,然而并非用以限定本发明,任何熟悉本专利的技术人员在不脱离本发明技术方案范围内,当可利用上述提示的技术内容作出些许更动或修饰为等同变化的等效实施例,但凡是未脱离本发明技术方案的内容,依据本发明的技术实质对以上实施例所作的任何简单修改、等同变化与修饰,均仍属于本发明方案的范围内。
Claims (11)
1.一种核电厂超温超功率回路保护通道定期试验验证准则的整定方法,其特征在于,将模拟的事故工况信号,通过计算得到一个输入电压V随时间t变化的函数,当所述输入电压V达到阈值电压V0时,将触发所述事故工况,再根据整个保护通道的系统偏差,得到触发所述事故工况动作的验证准则,所述影响事故工况的信号为输入电压在时间域t上的函数;所述的验证准则为时间验证准则:Tmax、Tmin和电压验证准则:Vmax,Vmin;
所述事故工况为:冷/热管段温度升高和稳压器压力升高而导致超温的事故工况;
所述阈值为超温阈值444XU1,所述阈值电压V0=0V;
所述系统偏差为:所述超温超功率回路保护通道允许漂移边界MA,MA=1.667%,所述超温超功率回路保护通道中电压值范围为:1V~5V,由此得到系统偏差VMA:VMA=4V*MA=4V*1.667%=0.067V,
根据所述输入电压V随时间t变化的函数:
得到:
由此可以得到式中所述时间t的取值范围动作时间上限T+和动作时间下限T-。
2.根据权利要求1所述的定期试验验证准则的整定方法,其特征在于,其特征在于,将模拟的事故工况信号,通过计算得到一个输入电压V随时间t变化的函数,当所述输入电压V达到阈值电压V0时,将触发所述事故工况,再根据整个保护通道的系统偏差,得到触发所述事故工况动作的验证准则,所述影响事故工况的信号为输入电压在时间域t上的函数;所述的验证准则为时间验证准则:Tmax、Tmin和电压验证准则:Vmax,Vmin;
所述事故工况为:冷/热管段温度升高和稳压器压力升高而导致超温的事故工况;
所述阈值为超温阈值444XU1,所述阈值电压V0=0V;
所述系统偏差为:所述超温超功率回路保护通道允许漂移边界MA,
MA=1.667%,所述超温超功率回路保护通道中电压值范围为:4V,由此得到系统偏差VMA:VMA=4V*MA=4V*1.667%=0.067V,
根据所述输入电压V随时间t变化的函数:
得到:
由此可以得到式中所述时间t的取值范围动作时间上限T+和动作时间下限T-。
3.根据权利要求1所述的定期试验验证准则的整定方法,其特征在于,
所述的模拟事故工况的信号为:
冷管段温度在时间域的函数Vc=0.111429t+2.565714,
热管段温度在时间域的函数Vh=0.037143t+3.605714,
稳压器压力在时间域的函数Vp=0.228571t+3.514286,
主泵转速在时间域的函数Vw,为常数4.3,
中子通量在时间域的函数Vf1为常数1;
所述输入电压V随时间t变化的函数为:
<mfenced open = "" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<mn>10.103815</mn>
<mo>-</mo>
<mi>a</mi>
<mo>*</mo>
<mn>16.700889</mn>
<mo>-</mo>
<mi>b</mi>
<mo>*</mo>
<mn>1.672</mn>
<mo>-</mo>
<mi>a</mi>
<mo>*</mo>
<mn>0.248411</mn>
<mo>-</mo>
<mi>a</mi>
<mo>*</mo>
<mn>2.401308</mn>
<mo>*</mo>
<msup>
<mi>e</mi>
<mrow>
<mo>-</mo>
<mi>t</mi>
</mrow>
</msup>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mo>+</mo>
<mi>a</mi>
<mo>*</mo>
<mn>8.611582</mn>
<mo>*</mo>
<msup>
<mi>e</mi>
<mrow>
<mo>-</mo>
<mfrac>
<mn>1</mn>
<mn>4</mn>
</mfrac>
<mi>t</mi>
</mrow>
</msup>
<mo>+</mo>
<mi>c</mi>
<mo>*</mo>
<mn>0.074286</mn>
<mo>*</mo>
<mrow>
<mo>(</mo>
<mrow>
<mi>t</mi>
<mo>+</mo>
<msup>
<mi>e</mi>
<mrow>
<mo>-</mo>
<mi>t</mi>
</mrow>
</msup>
<mo>-</mo>
<mn>1</mn>
</mrow>
<mo>)</mo>
</mrow>
<mo>-</mo>
<mi>c</mi>
<mo>*</mo>
<mn>1.04</mn>
<mo>-</mo>
<mi>d</mi>
<mo>=</mo>
<mi>V</mi>
<mo>;</mo>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
所述计算过程为:
(1)首先,所述冷管段温度在时间域的函数Vc=0.111429t+2.565714经过一阶惯性环节415FI的滤波:
(1.1)由于所述一阶惯性环节415FI的滤波对定值无影响,因此假设冷管段注入信号的初始值为0,即不考虑2.565714,将Vc=0.111429t+2.565714经过拉氏变换后,得到:
(1.2)所述一阶惯性环节415FI的传导函数G415FI:
输入信号为经过拉式变换的冷管段注入信号
由一阶惯性环节415FI的输出信号X1:X1=Vc*G415FI,
得到:然后经过拉式反变换,并将所述初始值2.565714带入方程,得到:X1=0.111429(t+e-t-1)+2.565714;
(2)然后,所述热管段温度在时间域的函数Vh=0.037143t+3.605714经过一阶惯性环节416FI的滤波:
(2.1)由于所述一阶惯性环节416FI对定值无影响,因此假设热管段注入信号的初始值为0,即不考虑3.605714,将Vh=0.037143t+3.605714经过拉氏变换后,得到:
(2.2)所述一阶惯性环节416FI的传导函数G416FI:
输入信号为经过拉式变换的热管段注入信号
由一阶惯性环节416FI的输出信号X2:X2=Vh*G416FI,
得到:然后经过拉式反变换,并将所述初始值3.605714带入方程,得到:X2=0.037143(t+e-t-1)+3.605714;
(3)再将经过所述一阶惯性环节415FI和所述一阶惯性环节416FI滤波后的所述冷管段温度在时间域的函数和所述热管段温度在时间域的函数带入平均温度402ZO的传导公式:G402ZO=a*(X1+X2)+b,并经过计算,得到平均温度
V402ZO=a*0.148572*(t+e-t-1)+a*6.171428+b,其中,0.57≤a≤0.59、-0.517≤b≤-0.465;
(4)由404MT超前滞后校正得到温度响应处理后的平均温度函数,
(4.1)所述的404MT超前滞后的传导函数G404MT:其中T1为超前参数,取值为30,T2为滞后参数,取值为4,因此得到:
(4.2)由于所述超前滞后环节404MT对定值无影响,因此假设平均温度初始值为0,即不考虑-a*0.148572+a*6.171428+b,经过拉氏变换后,得到
(4.3)由于所述超前滞后环节404MT的输入为:
V402ZO=a*0.148572*(t+e-t-1)+a*6.171428+b,
传导函数为:因此可得出所述超前滞后404MT的输出值为:
(4.4)在将式进行拉式反变换,并带入初始值-a*0.148572+a*6.171428+b后,得出:
<mrow>
<msub>
<mi>V</mi>
<mrow>
<mn>404</mn>
<mi>M</mi>
<mi>T</mi>
</mrow>
</msub>
<mo>=</mo>
<mi>a</mi>
<mo>*</mo>
<mn>9.885699</mn>
<mo>+</mo>
<mi>b</mi>
<mo>+</mo>
<mi>a</mi>
<mo>*</mo>
<mn>0.148571</mn>
<mo>*</mo>
<mi>t</mi>
<mo>+</mo>
<mi>a</mi>
<mo>*</mo>
<mn>1.436189</mn>
<mo>*</mo>
<msup>
<mi>e</mi>
<mrow>
<mo>-</mo>
<mi>t</mi>
</mrow>
</msup>
<mo>-</mo>
<mi>a</mi>
<mo>*</mo>
<mn>5.150468</mn>
<mo>*</mo>
<msup>
<mi>e</mi>
<mrow>
<mo>-</mo>
<mfrac>
<mn>1</mn>
<mn>4</mn>
</mfrac>
<mi>t</mi>
</mrow>
</msup>
<mo>;</mo>
</mrow>
(5)将经过所述一阶惯性环节415FI和所述一阶惯性环节416FI滤波后的所述冷管段温度在时间域的函数和所述热管段温度在时间域的函数带入温度偏差403ZO的传导公式:G403ZO=c*(X1-X2)+d,并经过计算,得到平均温度
V403ZO=-c*0.074286*(t+e-t-1)+c*1.04+d,其中,-1.556≤c≤1.556、1.711≤d≤1.889;
(6)再将V404MT,稳压器信号Vp=0.228571t+3.514286,主泵转速信号Vw=4.3,中子通量Vf1=1,带入式:
V404ZO=6.924-1.672*V404MT+0.810*VP+0.267*Vw-Vf1中,得到
<mrow>
<msub>
<mi>V</mi>
<mrow>
<mn>404</mn>
<mi>Z</mi>
<mi>O</mi>
</mrow>
</msub>
<mo>=</mo>
<mn>9.918672</mn>
<mo>-</mo>
<mi>a</mi>
<mo>*</mo>
<mn>16.700889</mn>
<mo>-</mo>
<mi>b</mi>
<mo>*</mo>
<mn>1.672</mn>
<mo>+</mo>
<mn>0.185143</mn>
<mo>-</mo>
<mi>a</mi>
<mo>*</mo>
<mn>0.248411</mn>
</mrow>
<mrow>
<mo>-</mo>
<mi>a</mi>
<mo>*</mo>
<mn>2.401308</mn>
<mo>*</mo>
<msup>
<mi>e</mi>
<mrow>
<mo>-</mo>
<mi>t</mi>
</mrow>
</msup>
<mo>+</mo>
<mi>a</mi>
<mo>*</mo>
<mn>8.611582</mn>
<mo>*</mo>
<msup>
<mi>e</mi>
<mrow>
<mo>-</mo>
<mfrac>
<mn>1</mn>
<mn>4</mn>
</mfrac>
<mi>t</mi>
</mrow>
</msup>
</mrow>
(7)最后根据所述工况的输入电压:V=V404ZO-V403ZO,将所述步骤(5)和所述步骤(6)得到V404ZO、V403ZO带入,并经过计算,得到所述输入电压V随时间t变化的函数为:
<mfenced open = "" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<mn>10.103815</mn>
<mo>-</mo>
<mi>a</mi>
<mo>*</mo>
<mn>16.700889</mn>
<mo>-</mo>
<mi>b</mi>
<mo>*</mo>
<mn>1.672</mn>
<mo>-</mo>
<mi>a</mi>
<mo>*</mo>
<mn>0.248411</mn>
<mo>-</mo>
<mi>a</mi>
<mo>*</mo>
<mn>2.401308</mn>
<mo>*</mo>
<msup>
<mi>e</mi>
<mrow>
<mo>-</mo>
<mi>t</mi>
</mrow>
</msup>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mo>+</mo>
<mi>a</mi>
<mo>*</mo>
<mn>8.611582</mn>
<mo>*</mo>
<msup>
<mi>e</mi>
<mrow>
<mo>-</mo>
<mfrac>
<mn>1</mn>
<mn>4</mn>
</mfrac>
<mi>t</mi>
</mrow>
</msup>
<mo>+</mo>
<mi>c</mi>
<mo>*</mo>
<mn>0.074286</mn>
<mo>*</mo>
<mrow>
<mo>(</mo>
<mrow>
<mi>t</mi>
<mo>+</mo>
<msup>
<mi>e</mi>
<mrow>
<mo>-</mo>
<mi>t</mi>
</mrow>
</msup>
<mo>-</mo>
<mn>1</mn>
</mrow>
<mo>)</mo>
</mrow>
<mo>-</mo>
<mi>c</mi>
<mo>*</mo>
<mn>1.04</mn>
<mo>-</mo>
<mi>d</mi>
<mo>=</mo>
<mi>V</mi>
<mo>.</mo>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
4.根据权利要求1所述的定期试验验证准则的整定方法,其特征在于,所述的时间验证准则还包括触发所述事故工况阈值模块动作的响应时间Trelay和所述信号的采样时间Tsample因此,所述时间验证准则中Tmax:Tmax=T++Trelay,和Tmin:Tmin=T-+Trelay+Tsample;优选的,所述事故工况阈值模块动作的响应时间Trelay=27.75ms,所述信号的采样时间Tsample=3.6ms。
5.根据权利要求1-4所述的任一定期试验验证准则的整定方法,其特征在于,所述的电压验证准则,还包括所述阈值444XU1允许的漂移电压所述超温超功率回路保护通道的电压范围,为1V~5V,所述漂移电压为:Vxu为所述超温超功率回路保护通道允许的漂移量,Vxu=0.41%,因此,得到:
所述的电压验证准则,还包括电压变化率D在时间域上的函数:
所述的电压变化率D在时间域上的函数,还包括触发所述事故工况后的动作时的响应时间Trelay和所述信号的采样时间Tsample;
因此所述的电压验证准则为:
6.根据权利要求5所述的任一定期试验验证准则的整定方法,其特征在于,所述的电压验证准则,还包括所述阈值444XU1允许的漂移电压所述超温超功率回路保护通道的电压范围,为4V,所述漂移电压为:Vxu为所述超温超功率回路保护通道允许的漂移量,Vxu=0.41%,因此,得到:
所述的电压验证准则,还包括电压变化率D在时间域上的函数:
所述的电压变化率D在时间域上的函数,还包括触发所述事故工况后的动作时的响应时间Trelay和所述信号的采样时间Tsample;
因此所述的电压验证准则为:
7.根据权利要求1所述的定期试验验证准则的整定方法,其特征在于,还包括与所述超温阈值模块互为备用的超温阈值模块444XU2,所述超温阈值模块444XU2与所述超温阈值模块444XU1实现相同的功能。
8.根据权利要求1-4所述的任一定期试验验证准则的整定方法,其特征在于,所述事故工况还包括:
1)中子通量升高而导致平均温度高的事故工况,所述阈值为443XU1,动作电压为0V,所述模拟事故工况的信号为:
冷管段温度在时间域的函数Vc=0.005714t+2.565714,
热管段温度在时间域的函数Vh=0.022857t+4.005714,
稳压器压力在时间域的函数Vp=0.028571t+3.514286,
主泵转速在时间域的函数Vw,为常数4.3,
中子通量在时间域的函数为Vf1=0.025t+0.15;
2)主泵转速降低而导致超功率的事故工况,所述阈值为445XU1,动作电压为0V,所述模拟事故工况的信号为:
冷管段温度在时间域的函数为常数,Vc=2.44,
热管段温度在时间域的函数为常数,Vh=3.477714,
稳压器压力在时间域的函数为常数Vp=3.514286,
主泵转速在时间域的函数,Vw=-0.222222t+4.3,
中子通量在时间域的函数Vf1为常数3.35。
9.根据权利要求5所述的定期试验验证准则的整定方法,其特征在于,所述事故工况还包括:
1)中子通量升高而导致平均温度高的事故工况,所述阈值为443XU1,动作电压为0V,所述模拟事故工况的信号为:
冷管段温度在时间域的函数Vc=0.005714t+2.565714,
热管段温度在时间域的函数Vh=0.022857t+4.005714,
稳压器压力在时间域的函数Vp=0.028571t+3.514286,
主泵转速在时间域的函数Vw,为常数4.3,
中子通量在时间域的函数为Vf1=0.025t+0.15;
2)主泵转速降低而导致超功率的事故工况,所述阈值为445XU1,动作电压为0V,所述模拟事故工况的信号为:
冷管段温度在时间域的函数为常数,Vc=2.44,
热管段温度在时间域的函数为常数,Vh=3.477714,
稳压器压力在时间域的函数为常数Vp=3.514286,
主泵转速在时间域的函数,Vw=-0.222222t+4.3,
中子通量在时间域的函数Vf1为常数3.35。
10.根据权利要求8所述的定期试验验证准则的整定方法,其特征在于,所述阈值模块还包括:与所述平均温度阈值模块443XU1互为备用的平均温度阈值模块443XU2;与所述超功率阈值模块445XU1互为备用的超功率阈值模块445XU2。
11.根据权利要求9所述的定期试验验证准则的整定方法,其特征在于,所述阈值模块还包括:与所述平均温度阈值模块443XU1互为备用的平均温度阈值模块443XU2;与所述超功率阈值模块445XU1互为备用的超功率阈值模块445XU2。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410380428.3A CN105321588B (zh) | 2014-08-05 | 2014-08-05 | 一种核电厂超温超功率回路保护通道定期试验验证准则的整定方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410380428.3A CN105321588B (zh) | 2014-08-05 | 2014-08-05 | 一种核电厂超温超功率回路保护通道定期试验验证准则的整定方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN105321588A CN105321588A (zh) | 2016-02-10 |
CN105321588B true CN105321588B (zh) | 2017-09-01 |
Family
ID=55248798
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201410380428.3A Active CN105321588B (zh) | 2014-08-05 | 2014-08-05 | 一种核电厂超温超功率回路保护通道定期试验验证准则的整定方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105321588B (zh) |
Families Citing this family (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106448775B (zh) * | 2016-09-23 | 2018-09-11 | 中广核工程有限公司 | 一种核电站反应堆保护系统安全参数监控装置与方法 |
CN108172311A (zh) * | 2017-11-28 | 2018-06-15 | 广东核电合营有限公司 | 一种反应堆的超温超功保护系统 |
CN111964931B (zh) * | 2020-08-12 | 2022-08-09 | 哈工智慧(武汉)科技有限公司 | 一种基于互联网的井盖监测方法、系统、存储介质及设备 |
CN113359895B (zh) * | 2021-06-23 | 2022-07-26 | 南方电网电力科技股份有限公司 | 一种壁管壁超温保护系统及其控制方法 |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2003205385A (ja) * | 2002-01-09 | 2003-07-22 | Hitachi Ltd | 構造部材補修方法および構造部材補修装置 |
CN101404409A (zh) * | 2008-11-17 | 2009-04-08 | 国电南瑞科技股份有限公司 | 电气化铁路at供电接触网断线接地故障识别方法 |
CN101533266A (zh) * | 2009-04-02 | 2009-09-16 | 余瑶 | 自控调节算法和控制仪表 |
CN102621883A (zh) * | 2012-04-01 | 2012-08-01 | 广东电网公司电力科学研究院 | Pid参数整定方法及pid参数整定系统 |
CN103064447A (zh) * | 2013-01-18 | 2013-04-24 | 西南科技大学 | 基于窄域论的激光器pid温控参数近似三维整定方法 |
CN103474949A (zh) * | 2013-08-07 | 2013-12-25 | 南京因泰莱配电自动化设备有限公司 | 一种自适应三段保护的整定方法以及装置 |
-
2014
- 2014-08-05 CN CN201410380428.3A patent/CN105321588B/zh active Active
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2003205385A (ja) * | 2002-01-09 | 2003-07-22 | Hitachi Ltd | 構造部材補修方法および構造部材補修装置 |
CN101404409A (zh) * | 2008-11-17 | 2009-04-08 | 国电南瑞科技股份有限公司 | 电气化铁路at供电接触网断线接地故障识别方法 |
CN101533266A (zh) * | 2009-04-02 | 2009-09-16 | 余瑶 | 自控调节算法和控制仪表 |
CN102621883A (zh) * | 2012-04-01 | 2012-08-01 | 广东电网公司电力科学研究院 | Pid参数整定方法及pid参数整定系统 |
CN103064447A (zh) * | 2013-01-18 | 2013-04-24 | 西南科技大学 | 基于窄域论的激光器pid温控参数近似三维整定方法 |
CN103474949A (zh) * | 2013-08-07 | 2013-12-25 | 南京因泰莱配电自动化设备有限公司 | 一种自适应三段保护的整定方法以及装置 |
Also Published As
Publication number | Publication date |
---|---|
CN105321588A (zh) | 2016-02-10 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN105321588B (zh) | 一种核电厂超温超功率回路保护通道定期试验验证准则的整定方法 | |
CN103617816B (zh) | 反应堆堆芯功率分布的测量方法 | |
CN105045233B (zh) | 火电厂热工系统中基于时间量度的pid控制器的优化设计方法 | |
CN105115628B (zh) | 一种热电阻动态响应测试系统及其测试方法 | |
Karamali et al. | A distributed solar collector field temperature profile control and estimation using inlet oil temperature and radiation estimates based on iterative extended kalman filter | |
CN103246279B (zh) | 一种存在执行阀粘滞特性的化工过程的控制性能评估方法 | |
CN107272412A (zh) | 一种暂冲式风洞流场控制的辩识方法 | |
Feng et al. | STRUCTure-based URANS simulations of thermal mixing in T-junctions | |
CN103198184A (zh) | 一种电力系统低频振荡特征类噪声辨识方法 | |
Sugumar et al. | Testing the effectiveness of attack detection mechanisms in industrial control systems | |
CN103869846B (zh) | 一种主回路冷热段平均温度信号处理方法 | |
CN104503260A (zh) | 调速器参数设定方法及装置 | |
CN104535257A (zh) | 一种硅压阻温度补偿评估方法 | |
CN105387449A (zh) | 一种在锅炉蒸汽温度控制中使用二阶微分的控制方法 | |
CN103871518B (zh) | 一种核电厂超温超功率保护温差信号处理方法 | |
CN106599492B (zh) | 一种基于逻辑回归的飞行器颤振分析及其qmu评估方法 | |
CN105422086A (zh) | 吸汽剖面获取方法及装置 | |
CN204270138U (zh) | 汽机调速系统伺服卡控制参数测试系统 | |
CN102938022A (zh) | 一种基于摄动法的频率相关网络等值的无源性校正方法 | |
CN204101980U (zh) | 燃气轮机调速器伺服系统控制参数测试系统 | |
Zhao et al. | Extended forward sensitivity analysis for uncertainty quantification | |
Shih et al. | A Strategy for Very Large Eddy Simulation of Complex Turbulent Flows | |
Zhou et al. | The preliminary research of pressure control system danymic simulation for Ap1000 pressurizer based on parameter adaptive fuzzy pid control algorithm | |
CN103199545B (zh) | 动态无功补偿装置最优二次高斯控制器及其设计方法 | |
van Paridon et al. | On the Real-Time Estimation of Disk Temperature Spatial Distributions in Aeroengines |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |