CN103426487A - 一种堆芯水位和过冷裕度通道互校方法 - Google Patents

一种堆芯水位和过冷裕度通道互校方法 Download PDF

Info

Publication number
CN103426487A
CN103426487A CN2012101637982A CN201210163798A CN103426487A CN 103426487 A CN103426487 A CN 103426487A CN 2012101637982 A CN2012101637982 A CN 2012101637982A CN 201210163798 A CN201210163798 A CN 201210163798A CN 103426487 A CN103426487 A CN 103426487A
Authority
CN
China
Prior art keywords
sat
vsl
avg
nargin
cold
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
Application number
CN2012101637982A
Other languages
English (en)
Other versions
CN103426487B (zh
Inventor
李小芬
苟拓
何正熙
余俊辉
李白
王远兵
王华金
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Nuclear Power Institute of China
Original Assignee
Nuclear Power Institute of China
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Nuclear Power Institute of China filed Critical Nuclear Power Institute of China
Priority to CN201210163798.2A priority Critical patent/CN103426487B/zh
Publication of CN103426487A publication Critical patent/CN103426487A/zh
Application granted granted Critical
Publication of CN103426487B publication Critical patent/CN103426487B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

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
    • Y02E30/30—Nuclear fission reactors

Landscapes

  • Monitoring And Testing Of Nuclear Reactors (AREA)

Abstract

本发明涉及一种堆芯水位和过冷裕度通道互校方法。本发明基于堆芯最低过冷裕度和反应堆压力容器水位测量通道不确定度和冗余系列数据组合与物理可能区域的交叉校核，通过获得A、B系列的测量数据的步骤、互校启动条件判断的步骤、信号有效性状态统计和初步计算的步骤、理论一致性判断的步骤、以及数据校核计算的步骤，将A、B列测量数据最终处理成唯一一组最可信的堆芯最低过冷裕度和压力容器水位数据，提高了数据可信度，简化了操纵员接口数量，同时避免了操纵员在遇到A、B列数据不一致时无法处理的情况。

Description

一种堆芯水位和过冷裕度通道互校方法
技术领域
本发明涉及一种堆芯水位和过冷裕度通道互校方法，特别是涉及一种核电厂事故后监测系统数据处理领域的堆芯水位和过冷裕度通道互校方法。
背景技术
根据文献检索，在《核动力工程》2008Vol.29No.1第5~9页中发表的《秦山核电二期扩建工程堆芯冷却监测系统设计》与本发明较为接近。
在传统的基于事件导向规程（EOP）的电厂中，反应堆堆芯最低过冷裕度和反应堆压力容器水位测量系统均是根据三里岛事故后经验反馈进行设计。三里岛事故后，为了在事故条件下知道堆芯最低过冷裕度和水位变化趋势，各国核安全当局要求新建核电厂必须增加直接监测堆芯最低过冷裕度和水位的手段，并需要满足RG1.97 Category 1（1级）参数要求。根据RG1.97要求，该系统设置了冗余的A、B两个系列，且两个系列间完全独立运行。但是随着状态导向规程（SOP）的逐渐成熟和大量应用，操纵员越来越依靠堆芯最低过冷裕度和水位提供的堆芯冷却状态信息进行规程判断，这对堆芯最低过冷裕度和水位监测数据的准确性要求越来越高，且要避免冗余系列出现数据不一致的情况，否则对操纵员规程选取带来障碍，影响电厂安全。因此亟需提供一种新型的堆芯水位和过冷裕度通道互校方法。
发明内容
本发明要解决的技术问题是提供一种满足状态导向规程的要求，进一步提高堆芯最低过冷裕度和水位监测数据的准确性，且尽可能地提供一组唯一的堆芯最低过冷裕度和水位数据供操纵员在选取规程时使用的堆芯水位和过冷裕度通道互校方法。
为解决上述技术问题，本发明一种堆芯水位和过冷裕度通道互校方法，依次包括以下步骤：
步骤一、获得A、B系列的测量数据；
通过A系列测量得到下列参数：
-目前正处于正常运行的主泵数量NRCP-A；
-堆芯最低过冷裕度ΔTSAT-A；
-堆芯最低过冷裕度有效性状态SΔTSAT-A；1表示有效，0表示无效；
-堆芯最低过冷裕度不确定度ΣΔTSAT-A；
-反应堆压力容器水位LVSL-A；
-反应堆压力容器水位有效性状态SLVSL-A；
通过B系列测量得到下列参数；
-目前正处于正常运行的主泵数量NRCP-B；
-堆芯最低过冷裕度ΔTSAT-B；
-堆芯最低过冷裕度有效性状态SΔTSAT-B；1表示有效，0表示无效；
-堆芯最低过冷裕度不确定度ΣΔTSAT-B；
-反应堆压力容器水位LVSL-B；
-反应堆压力容器水位有效性状态SLVSL-B；
A系列系统将上述测量到的数据通过点对点通信发送到B系列，同时B系列将上述测量到的数据通过点对点通信发送到A系列；
根据标准JJF1059-1999，得到反应堆压力容器水位测量通道的不确定度ΣLVSL；
步骤二、互校启动条件判断；
1)对于A系列系统，判断B系列到A系列的数据通信状态是否正常，如果通信状态不正常，则不进行数据互校处理，最终显示本系列原始测量数据；
2）对于B系列系统，判断A系列到B系列的数据通信状态是否正常，如果通信状态不正常，则不进行数据互校处理，最终显示本系列原始测量数据；
3）如果两个系列间数据通信正常，则判断NRCP-A是否与NRCP-B相等；如果不相等，A、B系列均不进行数据互校，最终显示本系列原始测量数据；如果相等，则进行步骤三；
步骤三、信号有效性状态统计和初步计算；
1)得到堆芯最低过冷裕度有效信号的个数VT，得到反应堆压力容器水位有效信号个数VL；
VT=SΔTSAT-A+SΔTSAT-B；
VL=SLVSL-A+SLVSL-B；
2)得到平均堆芯最低过冷裕度ΔTSAT-avg和平均堆芯最低过冷裕度不确定度ΣΔTSAT-avg；
Figure BDA00001679561100031
当SΔTSAT-A=1，SΔTSAT-B=0时，
ΔTSAT-avg=ΔTSAT-A，ΣΔTSAT-avg=ΣΔTSAT-A；
Figure BDA00001679561100032
当SΔTSAT-A=0，SΔTSAT-B=1时，
ΔTSAT-avg=ΔTSAT-B，ΣΔTSAT-avg=ΣΔTSAT-B；
Figure BDA00001679561100033
当SΔTSAT-A=1，SΔTSAT-B=1；
■当|ΔTSAT-A-ΔTSAT-B|≤ΣΔTSAT-A+ΣΔTSAT-B时，
ΔTSAT-avg=(ΔTSAT-A+ΔTSAT-B)/2，
ΣΔTSAT-avg=(ΣΔTSAT-A+ΣΔTSAT-B)/2；
■当|ΔTSAT-A-ΔTSAT-B|>ΣΔTSAT-A+ΣΔTSAT-B时，
◆当ΔTSAT-A≤ΔTSAT-B时，ΔTSAT-avg=ΔTSAT-A，ΣΔTSAT-avg=ΣΔTSAT-A；
◆当ΔTSAT-A>ΔTSAT-B时，ΔTSAT-avg=ΔTSAT-B，ΣΔTSAT-avg=ΣΔTSAT-B；
3)得到平均水位LVSL-avg；
当SLVSL-A=1，SLVSL-B=0时，LVSL-avg=LVSL-A；
Figure BDA00001679561100042
当SLVSL-A=0，SLVSL-B=1时，LVSL-avg=LVSL-B；
Figure BDA00001679561100043
当SLVSL-A=1，SLVSL-B=1时，
■当|LVSL-A-LVSL-B|≤2×ΣLVSL；LVSL-avg=(LVSL-A+LVSL-B)/2
■当|LVSL-A-LVSL-B|>2×ΣLVSL；LVSL-avg=min(LVSL-A，LVSL-B)
Figure BDA00001679561100044
当SLVSL-A=0，SLVSL-B=0时，
(1)当没有主泵运行时；
√若ΔTSAT-avg≥ΣΔTSAT-avg+24，LVSL-F=100%；
√若ΣΔTSAT-avg≤ΔTSAT-avg<ΣΔTSAT-avg+24，LVSL-F=64%；
√若-ΣΔTSAT-avg≤ΔTSAT-avg<ΣΔTSAT-avg，LVSL-F=55%；
√若-200℃≤ΔTSAT-avg<-ΣΔTSAT-avg，LVSL-F=40%；
√若ΔTSAT-avg<-200℃，LVSL-F=0%；
(2)当至少要有一台主泵运行时；
√若ΔTSAT-avg≥ΣΔTSAT-avg+24，LVSL-F=100%；
√若ΣΔTSAT-avg≤ΔTSAT-avg<ΣΔTSAT-avg+24，LVSL-F=76%；
√若-ΣΔTSAT-avg≤ΔTSAT-avg<ΣΔTSAT-avg，LVSL-F=60%；
√若ΔTSAT-avg<-ΣΔTSAT-avg，LVSL-F=0%；
步骤四：理论一致性判断；
根据步骤三中得到的平均堆芯最低过冷裕度不确定度ΣΔTSAT-avg在“堆芯最低过冷裕度─反应堆压力容器水位”坐标系中建立“物理不可能区域”，“物理不可能区域”定义如下：
（1）如果至少有一台主泵运行，物理不可能区域定义为：
Figure BDA00001679561100045
堆芯最低过冷裕度小于-ΣΔTSAT-avg且反应堆压力容器水位高于堆芯底部；
堆芯最低过冷裕度大于ΣΔTSAT-avg+24且反应堆压力容器水位低于堆芯底部；
（2）如果所有主泵都停止运行，物理不可能区域定义为：
Figure BDA00001679561100052
堆芯最低过冷裕度小于-ΣΔTSAT-avg且反应堆压力容器水位高于热管段顶部；
Figure BDA00001679561100053
堆芯最低过冷裕度大于ΣΔTSAT-avg且反应堆压力容器水位低于堆芯顶部；
根据A、B列测量数据，形成如下四个数组：（ΔTSAT-A，LVSL-A），（ΔTSAT-A，LVSL-B），（ΔTSAT-B，LVSL-A），（ΔTSAT-B，LVSL-B）；这四个数组中，如果任何一个数组含有无效数据分量，则删除该数组；将剩余有效的数组放入“堆芯最低过冷裕度─反应堆压力容器水位”坐标系中，得到落入“物理不可能区域”的组合数Nbad；
步骤五：数据校核计算，得到最可信的唯一一组（ΔTSAT-V，LVSL-V）；
Figure BDA00001679561100054
当VT=1，VL=0时，ΔTSAT-V=ΔTSAT-avg，LVSL-V=LVSL-F；
Figure BDA00001679561100055
当VT=2，VL=0时，ΔTSAT-V=ΔTSAT-avg，LVSL-V=LVSL-F；
Figure BDA00001679561100056
当VT=1，VL=1时，ΔTSAT-V=ΔTSAT-avg，LVSL-V=LVSL-avg；
Figure BDA00001679561100057
当VT=2，VL=1时；
■如果Nbad=0，ΔTSAT-V=ΔTSAT-avg，LVSL-V=LVSL-avg；
■如果Nbad=1，将落入“物理不可能区域”的组合中的堆芯最低过冷裕度强制为无效，重新执行步骤三的计算，更新ΔTSAT-avg和ΔLVSL-avg，且
ΔTSAT-V=ΔTSAT-avg，LVSL-V=LVSL-avg；
■如果Nbad=2，ΔTSAT-V=min(ΔTSAT-A，ΔTSAT-A)，LVSL-V=LVSL-avg；
Figure BDA00001679561100058
当VT=1，VL=2时；
■如果Nbad=0，ΔTSAT-V=ΔTSAT-avg，LVSL-V=LVSL-avg；
■如果Nbad=1，将落入“物理不可能区域”的组合中的反应堆压力容器水位强制为无效，重新执行步骤三的计算，更新ΔTSAT-avg和LVSL-avg，且ΔTSAT-V=ΔTSAT-avg，LVSL-V=LVSL-avg；
■如果Nbad=2，ΔTSAT-V=ΔTSAT-avg，LVSL-V=min(LVSL-A,LVSL-B)；
Figure BDA00001679561100061
当VT=2，VL=2时；
■如果Nbad=0，ΔTSAT-V=ΔTSAT-avg，LVSL-V=LVSL-avg；
■如果Nbad=1，将落入“物理不可能区域”的组合中的堆芯最低过冷裕度和反应堆压力容器水位强制为无效，重新执行步骤三的计算，更新ΔTSAT-avg和LVSL-avg，且ΔTSAT-V=ΔTSAT-avg，LVSL-V=LVSL-avg；
■如果Nbad=2，
◆如果落入“物理不可能区域”的两个数组中有相同的反应堆压力容器水位，则该反应堆压力容器水位强制为无效，重新执行步骤三的计算，更新ΔTSAT-avg和LVSL-avg，且ΔTSAT-V=ΔTSAT-avg，LVSL-V=LVSL-avg；
◆如果落入“物理不可能区域”的两个数组中没有相同的堆芯最低过冷裕度，也没有相同的反应堆压力容器水位，则ΔTSAT-V=min(ΔTSAT-A，ΔTSAT-A)，LVSL-V=min(LVSL-A,LVSL-B)；
■如果Nbad=4，ΔTSAT-V=min(ΔTSAT-A，ΔTSAT-A)，LVSL-V=min(LVSL-A,LVSL-B)。
本发明通过基于堆芯最低过冷裕度和反应堆压力容器水位测量通道不确定度和冗余系列数据组合与物理可能区域的交叉校核，将A、B系列测量数据最终处理成唯一一组最可信的堆芯最低过冷裕度和压力容器水位数据，提高了数据可信度，简化了操纵员接口数量，同时避免了操纵员在遇到A、B列数据不一致时无法处理的情况。本发明能够及时将错误数据和漂移过大的数据剔除，并且提供唯一一组最可信的测量数据给操作员，简化了运行人员判断过程，避免了在冗余系列测量数据不一致时运行人员处理困难的情况，提高了电厂的安全性和经济性。
具体实施方式
本发明依次包括以下步骤：
步骤一、获得A、B系列的测量数据；
通过A系列测量得到下列参数：
-目前正处于正常运行的主泵数量NRCP-A；
-堆芯最低过冷裕度ΔTSAT-A；
-堆芯最低过冷裕度有效性状态SΔTSAT-A；1表示有效，0表示无效；
-堆芯最低过冷裕度不确定度ΣΔTSAT-A；
-反应堆压力容器水位LVSL-A；
-反应堆压力容器水位有效性状态SLVSL-A；
通过B系列测量得到下列参数；
-目前正处于正常运行的主泵数量NRCP-B；
-堆芯最低过冷裕度ΔTSAT-B；
-堆芯最低过冷裕度有效性状态SΔTSAT-B；1表示有效，0表示无效；
-堆芯最低过冷裕度不确定度ΣΔTSAT-B；
-反应堆压力容器水位LVSL-B；
-反应堆压力容器水位有效性状态SLVSL-B；
A系列系统将上述测量到的数据通过点对点通信发送到B系列，同时B系列将上述测量到的数据通过点对点通信发送到A系列；
A、B系列之间的通信满足电气隔离准则；
根据标准JJF1059-1999，得到反应堆压力容器水位测量通道的不确定度ΣLVSL；
步骤二、互校启动条件判断；
1)对于A系列系统，判断B系列到A系列的数据通信状态是否正常，如果通信状态不正常，则不进行数据互校处理，最终送主控室显示的是本系列原始测量数据（非互校后数据），并提示操纵员；
2）对于B系列系统，判断A系列到B系列的数据通信状态是否正常，如果通信状态不正常，则不进行数据互校处理，最终送主控室显示的是本系列原始测量数据（非互校后数据），并提示操纵员；
3）如果两个系列间数据通信正常，则判断NRCP-A是否与NRCP-B相等；如果不相等，A、B系列均不进行数据互校，分别显示本系列原始测量数据（非互校后数据），并提示操纵员；如果相等，则进行步骤三；
步骤三、信号有效性状态统计和初步计算；
1)得到堆芯最低过冷裕度有效信号的个数VT，得到反应堆压力容器水位有效信号个数VL；
VT=SΔTSAT-A+SΔTSAT-B；
VL=SLVSL-A+SLVSL-B；
2)得到平均堆芯最低过冷裕度ΔTSAT-avg和平均堆芯最低过冷裕度不确定度ΣΔTSAT-avg；
Figure BDA00001679561100081
当SΔTSAT-A=1，SΔTSAT-B=0时，ΔTSAT-avg=ΔTSAT-A，ΣΔTSAT-avg=ΣΔTSAT-A；
Figure BDA00001679561100082
当SΔTSAT-A=0，SΔTSAT-B=1时，ΔTSAT-avg=ΔTSAT-B，ΣΔTSAT-avg=ΣΔTSAT-B；
当SΔTSAT-A=1，SΔTSAT-B=1；
■当|ΔTSAT-A-ΔTSAT-B|≤ΣΔTSAT-A+ΣΔTSAT-B时，
ΔTSAT-avg=(ΔTSAT-A+ΔTSAT-B)/2，
ΣΔTSAT-avg=(ΣΔTSAT-A+ΣΔTSAT-B)/2；
■当|ΔTSAT-A-ΔTSAT-B|>ΣΔTSAT-A+ΣΔTSAT-B时，
◆当ΔTSAT-A≤ΔTSAT-B时，ΔTSAT-avg=ΔTSAT-A，ΣΔTSAT-avg=ΣΔTSAT-A；
◆当ΔTSAT-A>ΔTSAT-B时，ΔTSAT-avg=ΔTSAT-B，ΣΔTSAT-avg=ΣΔTSAT-B；
注：不考虑SΔTSAT-A=0同时SΔTSAT-B=0的情况；
3)得到平均水位LVSL-avg；
Figure BDA00001679561100091
当SLVSL-A=1，SLVSL-B=0时，LVSL-avg＝LVSL-A；
Figure BDA00001679561100092
当SLVSL-A=0，SLVSL-B=1时，LVSL-avg=LVSL-B；
Figure BDA00001679561100093
当SLVSL-A=1，SLVSL-B=1时，
■当|LVSL-A-LVSL-B|≤2×ΣLVSL；LVSL-avg=(LVSL-A+LVSL-B)/2
■当|LVSL-A-LVSL-B|>2×ΣLVSL；LVSL-avg=min(LVSL-A，LVSL-B)
Figure BDA00001679561100094
当SLVSL-A=0，SLVSL-B=0时，
LVSL-A和LVSL-B都无效，则需要根据ΔTSAT-avg计算缺省值水位LVSL-F；
(1)当没有主泵运行时；
√若ΔTSAT-avg≥ΣΔTSAT-avg+24，LVSL-F=100%；
√若ΣΔTSAT-avg≤ΔTSAT-avg<ΣΔTSAT-avg+24，LVSL-F=64%；
√若-ΣΔTSAT-avg≤ΔTSAT-avg<ΣΔTSAT-avg，LVSL-F=55%；
√若-200℃≤ΔTSAT-avg<-ΣΔTSAT-avg，LVSL-F=40%；
√若ΔTSAT-avg<-200℃，LVSL-F=0%；
(2)当至少要有一台主泵运行时；
√若ΔTSAT-avg≥ΣΔTSAT-avg+24，LVSL-F=100%；
√若ΣΔTSAT-avg≤ΔTSAT-avg<ΣΔTSAT-avg+24，LVSL-F=76%；
√若-ΣΔTSAT-avg≤ΔTSAT-avg<ΣΔTSAT-avg，LVSL-F=60%；
√若ΔTSAT-avg<-ΣΔTSAT-avg，LVSL-F=0%；
步骤四：理论一致性判断；
根据步骤三中得到的平均堆芯最低过冷裕度不确定度ΣΔTSAT-avg在“堆芯最低过冷裕度─反应堆压力容器水位”坐标系中建立“物理不可能区域”，“物理不可能区域”定义如下：
（1）如果至少有一台主泵运行，物理不可能区域定义为：
Figure BDA00001679561100101
堆芯最低过冷裕度小于-ΣΔTSAT-avg且反应堆压力容器水位高于堆芯底部；
堆芯最低过冷裕度大于ΣΔTSAT-avg+24且反应堆压力容器水位低于堆芯底部；
（2）如果所有主泵都停止运行，物理不可能区域定义为：
堆芯最低过冷裕度小于-ΣΔTSAT-avg且反应堆压力容器水位高于热管段顶部；
堆芯最低过冷裕度大于ΣΔTSAT-avg且反应堆压力容器水位低于堆芯顶部；
根据A、B列测量数据，形成如下四个数组：（ΔTSAT-A，LVSL-A），（ΔTSAT-A，LVSL-B），（ΔTSAT-B，LVSL-A），（ΔTSAT-B，LVSL-B）；这四个数组中，如果任何一个数组含有无效数据分量，则删除该数组；例如，如果B系列测量到的水位LVSL-B无效，即SLVSL-B=0，则（ΔTSAT-A，LVSL-B）和（ΔTSAT-B，LVSL-B）这两个数组就无效，应被删除；将剩余有效的数组放入堆芯最低过冷裕度─反应堆压力容器水位坐标系中，得到落入“物理不可能区域”的组合数Nbad；
步骤五：数据校核计算，得到最可信的唯一一组（ΔTSAT-V，LVSL-V）；
Figure BDA00001679561100105
当VT=1，VL=0时，ΔTSAT-V=ΔTSAT-avg，LVSL-V=LVSL-F；
当VT=2，VL=0时，ΔTSAT-V=ΔTSAT-avg，LVSL-V=LVSL-F；
Figure BDA00001679561100107
当VT=1，VL=1时，ΔTSAT-V=ΔTSAT-avg，LVSL-V=LVSL-avg；
Figure BDA00001679561100108
当VT=2，VL=1时；
■如果Nbad=0，ΔTSAT-V=ΔTSAT-avg，LVSL-V=LVSL-avg；
■如果Nbad=1，将落入“物理不可能区域”的组合中的堆芯最低过冷裕度强制为无效，重新执行步骤三的计算，更新ΔTSAT-avg和ΔLVSL-avg，且ΔTSAT-V=ΔTSAT-avg，LVSL-V=LVSL-avg；
■如果Nbad=2，ΔTSAT-V=min(ΔTSAT-A，ΔTSAT-A)，LVSL-V=LVSL-avg；
Figure BDA00001679561100111
当VT=1，VL=2时；
■如果Nbad=0，ΔTSAT-V=ΔTSAT-avg，LVSL-V=LVSL-avg；
■如果Nbad=1，将落入“物理不可能区域”的组合中的反应堆压力容器水位强制为无效，重新执行步骤三的计算，更新ΔTSAT-avg和LVSL-avg，且ΔTSAT-V=ΔTSAT-avg，LVSL-V=LVSL-avg；
■如果Nbad=2，ΔTSAT-V=ΔTSAT-avg，LVSL-V=min(LVSL-A,LVSL-B)；
Figure BDA00001679561100112
当VT=2，VL=2时；
■如果Nbad=0，ΔTSAT-V=ΔTSAT-avg，LVSL-V=LVSL-avg；
■如果Nbad=1，将落入“物理不可能区域”的组合中的堆芯最低过冷裕度和反应堆压力容器水位强制为无效，重新执行步骤三的计算，更新ΔTSAT-avg和LVSL-avg，且ΔTSAT-V=ΔTSAT-avg，LVSL-V=LVSL-avg；
■如果Nbad=2
◆如果落入“物理不可能区域”的两个数组中有相同的反应堆压力容器水位，则该反应堆压力容器水位强制为无效，重新执行步骤三的计算，更新ΔTSAT-avg和LVSL-avg，且ΔTSAT-V=ΔTSAT-avg，LVSL-V=LVSL-avg；
◆如果落入“物理不可能区域”的两个数组中没有相同的堆芯最低过冷裕度，也没有相同的反应堆压力容器水位，则ΔTSAT-V=min(ΔTSAT-A，ΔTSAT-A)，LVSL-V=min(LVSL-A,LVSL-B)；
■如果Nbad=4，ΔTSAT-V=min(ΔTSAT-A，ΔTSAT-A)，LVSL-V=min(LVSL-A,LVSL-B)；
注：Nbad不可能等于3。

Claims (1)

1.一种堆芯水位和过冷裕度通道互校方法，依次包括以下步骤：
步骤一、获得A、B系列的测量数据；
通过A系列测量得到下列参数：
-目前正处于正常运行的主泵数量NRCP-A；
-堆芯最低过冷裕度ΔTSAT-A；
-堆芯最低过冷裕度有效性状态SΔTSAT-A；1表示有效，0表示无效；
-堆芯最低过冷裕度不确定度ΣΔTSAT-A；
-反应堆压力容器水位LVSL-A；
-反应堆压力容器水位有效性状态SLVSL-A；
通过B系列测量得到下列参数；
-目前正处于正常运行的主泵数量NRCP-B；
-堆芯最低过冷裕度ΔTSAT-B；
-堆芯最低过冷裕度有效性状态SΔTSAT-B；1表示有效，0表示无效；
-堆芯最低过冷裕度不确定度ΣΔTSAT-B；
-反应堆压力容器水位LVSL-B；
-反应堆压力容器水位有效性状态SLVSL-B；
A系列系统将上述测量到的数据通过点对点通信发送到B系列，同时B系列将上述测量到的数据通过点对点通信发送到A系列；
根据标准JJF1059-1999，得到反应堆压力容器水位测量通道的不确定度ΣLVSL；
步骤二、互校启动条件判断；
1)对于A系列系统，判断B系列到A系列的数据通信状态是否正常，如果通信状态不正常，则不进行数据互校处理，最终显示本系列原始测量数据；
2）对于B系列系统，判断A系列到B系列的数据通信状态是否正常，如果通信状态不正常，则不进行数据互校处理，最终显示本系列原始测量数据；
3）如果两个系列间数据通信正常，则判断NRCP-A是否与NRCP-B相等；如果不相等，A、B系列均不进行数据互校，最终显示本系列原始测量数据；如果相等，则进行步骤三；
步骤三、信号有效性状态统计和初步计算；
1)得到堆芯最低过冷裕度有效信号的个数VT，得到反应堆压力容器水位有效信号个数VL；
VT=SΔTSAT-A+SΔTSAT-B；
VL=SLVSL-A+SLVSL-B；
2)得到平均堆芯最低过冷裕度ΔTSAT-avg和平均堆芯最低过冷裕度不确定度ΣΔTSAT-avg；
Figure FDA00001679561000021
当SΔTSAT-A=1，SΔTSAT-B=0时，
ΔTSAT-avg=ΔTSAT-A，ΣΔTSAT-avg=ΣΔTSAT-A；
Figure FDA00001679561000022
当SΔTSAT-A=0，SΔTSAT-B=1时，
ΔTSAT-avg=ΔTSAT-B，ΣΔTSAT-avg=ΣΔTSAT-B；
当SΔTSAT-A=1，SΔTSAT-B=1；
■当|ΔTSAT-A-ΔTSAT-B|≤ΣΔTSAT-A+ΣΔTSAT-B时，
ΔTSAT-avg=(ΔTSAT-A+ΔTSAT-B)/2，
ΣΔTSAT-avg=(ΣΔTSAT-A+ΣΔTSAT-B)/2；
■当|ΔTSAT-A-ΔTSAT-B|>ΣΔTSAT-A+ΣΔTSAT-B时，
◆当ΔTSAT-A≤ΔTSAT-B时，ΔTSAT-avg=ΔTSAT-A，ΣΔTSAT-avg=ΣΔTSAT-A；
◆当ΔTSAT-A>ΔTSAT-B时，ΔTSAT-avg=ΔTSAT-B，ΣΔTSAT-avg=ΣΔTSAT-B；
3)得到平均水位LVSL-avg；
Figure FDA00001679561000031
当SLVSL-A=1，SLVSL-B=0时，LVSL-avg=LVSL-A；
当SLVSL-A=0，SLVSL-B=1时，LVSL-avg=LVSL-B；
Figure FDA00001679561000033
当SLVSL-A=1，SLVSL-B=1时，
■当|LVSL-A-LVSL-B|≤2×ΣLVSL；LVSL-avg=(LVSL-A+LVSL-B)/2
■当|LVSL-A-LVSL-B|>2×ΣLVSL；LVSL-avg=min(LVSL-A，LVSL-B)
Figure FDA00001679561000034
当SLVSL-A=0，SLVSL-B=0时，
(1)当没有主泵运行时；
√若ΔTSAT-avg≥ΣΔTSAT-avg+24，LVSL-F=100%；
√若ΣΔTSAT-avg≤ΔTSAT-avg<ΣΔTSAT-avg+24，LVSL-F=64%；
√若-ΣΔTSAT-avg≤ΔTSAT-avg<ΣΔTSAT-avg，LVSL-F=55%；
√若-200℃≤ΔTSAT-avg<-ΣΔTSAT-avg，LVSL-F=40%；
√若ΔTSAT-avg<-200℃，LVSL-F=0%；
(2)当至少要有一台主泵运行时；
√若ΔTSAT-avg≥ΣΔTSAT-avg+24，LVSL-F=100%；
√若ΣΔTSAT-avg≤ΔTSAT-avg<ΣΔTSAT-avg+24，LVSL-F=76%；
√若-ΣΔTSAT-avg≤ΔTSAT-avg<ΣΔTSAT-avg，LVSL-F=60%；
√若ΔTSAT-avg<-ΣΔTSAT-avg，LVSL-F=0%；
步骤四：理论一致性判断；
根据步骤三中得到的平均堆芯最低过冷裕度不确定度ΣΔTSAT-avg在“堆芯最低过冷裕度─反应堆压力容器水位”坐标系中建立“物理不可能区域”，“物理不可能区域”定义如下：
（1）如果至少有一台主泵运行，物理不可能区域定义为：
堆芯最低过冷裕度小于-ΣΔTSAT-avg且反应堆压力容器水位高于堆芯底部；
Figure FDA00001679561000042
堆芯最低过冷裕度大于ΣΔTSAT-avg+24且反应堆压力容器水位低于堆芯底部；
（2）如果所有主泵都停止运行，物理不可能区域定义为：
Figure FDA00001679561000043
堆芯最低过冷裕度小于-ΣΔTSAT-avg且反应堆压力容器水位高于热管段顶部；
Figure FDA00001679561000044
堆芯最低过冷裕度大于ΣΔTSAT-avg且反应堆压力容器水位低于堆芯顶部；
根据A、B列测量数据，形成如下四个数组：（ΔTSAT-A，LVSL-A），（ΔTSAT-A，LVSL-B），（ΔTSAT-B，LVSL-A），（ΔTSAT-B，LVSL-B）；这四个数组中，如果任何一个数组含有无效数据分量，则删除该数组；将剩余有效的数组放入“堆芯最低过冷裕度─反应堆压力容器水位”坐标系中，得到落入“物理不可能区域”的组合数Nbad；
步骤五：数据校核计算，得到最可信的唯一一组（ΔTSAT-V，LVSL-V）；
Figure FDA00001679561000051
当VT=1，VL=0时，ΔTSAT-V=ΔTSAT-avg，LVSL-V=LVSL-F；
Figure FDA00001679561000052
当VT=2，VL=0时，ΔTSAT-V=ΔTSAT-avg，LVSL-V=LVSL-F；
Figure FDA00001679561000053
当VT=1，VL=1时，ΔTSAT-V=ΔTSAT-avg，LVSL-V=LVSL-avg；
Figure FDA00001679561000054
当VT=2，VL=1时；
■如果Nbad=0，ΔTSAT-V=ΔTSAT-avg，LVSL-V=LVSL-avg；
■如果Nbad=1，将落入“物理不可能区域”的组合中的堆芯最低过冷裕度强制为无效，重新执行步骤三的计算，更新ΔTSAT-avg和ΔLVSL-avg，且ΔTSAT-V=ΔTSAT-avg，LVSL-V=LVSL-avg；
■如果Nbad=2，ΔTSAT-V=min(ΔTSAT-A，ΔTSAT-A)，LVSL-V=LVSL-avg；
Figure FDA00001679561000055
当VT=1，VL=2时；
■如果Nbad=0，ΔTSAT-V=ΔTSAT-avg，LVSL-V=LVSL-avg；
■如果Nbad=1，将落入“物理不可能区域”的组合中的反应堆压力容器水位强制为无效，重新执行步骤三的计算，更新ΔTSAT-avg和LVSL-avg，且ΔTSAT-V=ΔTSAT-avg，LVSL-V=LVSL-avg；
■如果Nbad=2，ΔTSAT-V=ΔTSAT-avg，LVSL-V=min(LVSL-A,LVSL-B)；
Figure FDA00001679561000056
当VT=2，VL=2时；
■如果Nbad=0，ΔTSAT-V=ΔTSAT-avg，LVSL-V=LVSL-avg；
■如果Nbad=1，将落入“物理不可能区域”的组合中的堆芯最低过冷裕度和反应堆压力容器水位强制为无效，重新执行步骤三的计算，更新ΔTSAT-avg和LVSL-avg，且ΔTSAT-V=ΔTSAT-avg，LVSL-V=LVSL-avg；
■如果Nbad=2，
◆如果落入“物理不可能区域”的两个数组中有相同的反应堆压力容器水位，则该反应堆压力容器水位强制为无效，重新执行步骤三的计算，更新ΔTSAT-avg和LVSL-avg，且ΔTSAT-V=ΔTSAT-avg，LVSL-V=LVSL-avg；
◆如果落入“物理不可能区域”的两个数组中没有相同的堆芯最低过冷裕度，也没有相同的反应堆压力容器水位，则ΔTSAT-V=min(ΔTSAT-A，ΔTSAT-A)，LVSL-V=min(LVSL-A,LVSL-B)；
■如果Nbad=4，ΔTSAT-V=min(ΔTSAT-A，ΔTSAT-A)，LVSL-V=min(LVSL-A,LVSL-B)。
CN201210163798.2A 2012-05-24 2012-05-24 一种堆芯水位和过冷裕度通道互校方法 Active CN103426487B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201210163798.2A CN103426487B (zh) 2012-05-24 2012-05-24 一种堆芯水位和过冷裕度通道互校方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201210163798.2A CN103426487B (zh) 2012-05-24 2012-05-24 一种堆芯水位和过冷裕度通道互校方法

Publications (2)

Publication Number Publication Date
CN103426487A true CN103426487A (zh) 2013-12-04
CN103426487B CN103426487B (zh) 2016-01-27

Family

ID=49651101

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201210163798.2A Active CN103426487B (zh) 2012-05-24 2012-05-24 一种堆芯水位和过冷裕度通道互校方法

Country Status (1)

Country Link
CN (1) CN103426487B (zh)

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106205754A (zh) * 2016-07-14 2016-12-07 中广核工程有限公司 一种用于检测核电站反应堆压力容器水位的装置和方法
CN106558351A (zh) * 2015-09-28 2017-04-05 苏州热工研究院有限公司 核电厂堆芯水位测量方法
CN106992029A (zh) * 2017-03-30 2017-07-28 中广核工程有限公司 一种核电站堆芯饱和裕度计算系统及方法

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5024802A (en) * 1990-10-22 1991-06-18 Westinghouse Electric Corp. Method for steam generator water level measurement
US5323430A (en) * 1991-05-17 1994-06-21 N.V. Tot Keuring Van Elektrotechnische Materialen (Kema) Method and system for fine control of thermal power of a boiling water reactor with natural circulation
US5365555A (en) * 1993-06-23 1994-11-15 General Electric Company Water level measurement system
US5426677A (en) * 1994-01-03 1995-06-20 Westinghouse Electric Corporation Device and method for measuring temperature of a liquid contained in a pressurizer vessel

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5024802A (en) * 1990-10-22 1991-06-18 Westinghouse Electric Corp. Method for steam generator water level measurement
US5323430A (en) * 1991-05-17 1994-06-21 N.V. Tot Keuring Van Elektrotechnische Materialen (Kema) Method and system for fine control of thermal power of a boiling water reactor with natural circulation
US5365555A (en) * 1993-06-23 1994-11-15 General Electric Company Water level measurement system
US5426677A (en) * 1994-01-03 1995-06-20 Westinghouse Electric Corporation Device and method for measuring temperature of a liquid contained in a pressurizer vessel

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
何正熙等: "秦山核电二期扩建工程堆芯冷却监测系统设计", 《核动力工程》, vol. 29, no. 1, 28 February 2008 (2008-02-28), pages 5 - 9 *

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106558351A (zh) * 2015-09-28 2017-04-05 苏州热工研究院有限公司 核电厂堆芯水位测量方法
CN106558351B (zh) * 2015-09-28 2018-06-29 苏州热工研究院有限公司 核电厂堆芯水位测量方法
CN106205754A (zh) * 2016-07-14 2016-12-07 中广核工程有限公司 一种用于检测核电站反应堆压力容器水位的装置和方法
CN106992029A (zh) * 2017-03-30 2017-07-28 中广核工程有限公司 一种核电站堆芯饱和裕度计算系统及方法

Also Published As

Publication number Publication date
CN103426487B (zh) 2016-01-27

Similar Documents

Publication Publication Date Title
CN102682862B (zh) 压水堆核电站堆芯热电偶有效性状态在线诊断方法
CN103426487A (zh) 一种堆芯水位和过冷裕度通道互校方法
CN106340332A (zh) 核电站数字化保护控制系统
CN104240781B (zh) 核电厂数字化仪控系统的信号分配方法及系统
CN106448764A (zh) 核电站乏燃料水池监测系统、方法
CN110415848A (zh) 一种应对排热减少叠加swccf事故的保护系统
CN110380514A (zh) 一种智能变电站继电保护二次回路故障诊断方法
CN106782719B (zh) 压水堆机组汽轮机已跳闸表征信号产生方法
CN109359414B (zh) 机械零部件关键部位测点优化布局下的应变实时监测方法
CN104240779B (zh) 核电站堆芯功率象限倾斜的监督方法和装置
CN103426488B (zh) 一种基于平均密度的反应堆压力容器水位测量方法
CN104330990A (zh) 一种偶校验电缆型号识别方法及装置
CN107887038A (zh) 核电站冷超压防护方法及装置
CN108172314A (zh) 数字化压水堆核电厂主控室不可用故障处理策略分析方法
CN106325131B (zh) 一种海上小型堆严重事故仪控系统及其控制方法
CN103983879A (zh) 一种油浸式变压器在线故障诊断方法及模拟装置
CN110415849A (zh) 一种应对排热增加叠加swccf事故的保护系统
CN106816191A (zh) 中子温度测量通道与电缆连接正确性的检查方法
CN103871512A (zh) 一种用于反应堆堆芯饱和温度计算的主回路压力校核方法
CN109712729B (zh) 一种核电站报警方法、饱和裕度的计算方法及系统
CN104332193B (zh) 一种基于数字化技术的蒸汽发生器水位测量方法
CN110706834B (zh) 一种堆芯冷却监视信号有效性处理方法及装置
CN104181420B (zh) 核电站atws系统的响应时间测试方法和系统
CN106531268A (zh) 一种用于核电厂数字化事故程序的链接布置方法
CN105160171B (zh) 一种基于多状态转换推理的电子系统热可靠性分析及预测方法

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