一种基于安全壳压升监测的安全壳泄漏快速判断方法

1.本发明涉及核电厂安全技术领域，具体涉及一种基于安全壳压升监测的安全壳泄漏快速判断方法。

背景技术：

2.安全壳是核反应堆的第三道安全屏障，外形为带圆穹顶的圆柱形筒体，自由容积达5万立方米，属于大型容器设备，其结构上贯穿件数量众多，加之混凝土气孔和裂缝的客观存在，当安全壳内气压高于大气压时，安全壳内空气会通过这些缝隙向外泄漏，对环境造成影响。
3.为了控制安全壳内气体泄漏，核电厂《运行技术规范》中规定了安全壳在正常运行状况下的泄漏率限值和控制准则，即在60mbar.g下反应堆厂房的正向泄漏率大于5n
·
m3/h以及小于10n
·
m3/h时，应查并消除泄漏的原因，否则14天内机组开始向mcs模式后撤。目前，我国改进型百万千瓦级(1000mw)压水堆核电技术(cpr1000)，安全壳泄漏监测系统(epp)需要经过筛选后至少五天的安全壳内气体监测数据，才可以计算安全壳的泄漏率。在机组启动期间，数据量不足导致无法计算安全壳泄漏率，无法评估安全壳的密封状况，作为安全壳密封性、机组运行安全性的重要指标之一，有必要掌握安全壳在机组启动阶段的泄漏率，以了解安全壳密封性的变化，确保泄漏率在限值内。
4.在核电机组启动或运行期间，引起安全壳内气体压力变化的主要因素包括五方面：安全壳内气体的温度、湿度变化、控制安全壳内气动调节阀所用的核岛仪表压缩空气系统(sar)的压缩空气因调节阀的动作所注入安全壳的空气、开启气闸门时安全壳与气闸间内空气的交换和通过安全壳的气体泄漏。在机组启动的短时间内安全壳内压力变化明显，目前尚缺乏机组启动阶段，基于压力增长特性来快速判定泄漏情况的方法。

技术实现要素：

5.针对以上问题，本发明要解决的技术问题在于，提供一种基于安全壳压升监测的安全壳泄漏快速判断方法。
6.本发明的技术方案如下：
7.一种基于安全壳压升监测的安全壳泄漏快速判断方法，其特征在于，可分为定性判断法和定量判断法。安全壳内湿空气组分满足质量平衡关系、道尔顿分压定律和理想气体状态方程。安全壳内从t1时刻到t2时刻，采用以下干空气质量平衡关系：
8.m
t2
，＝m
t1
，+m
sar
+δm，+m
leak
，
ꢀꢀꢀ
(1)
9.式中，m
t2，空
为t2时刻安全壳内干空气质量；m
t1，空
为t1时刻安全壳内干空气质量；m
sar
为从t1时刻到t2时刻δt时间段内sar系统注入安全壳的干空气质量(可由sar系统检测数据获取)；δm
闸，空
为δt时间段内安全壳与气闸间交换干空气质量，从壳内到气闸间的数值为负；m
leak，空
是δt时间段内安全壳的干空气泄漏量，从壳内泄漏到壳外时数值为负。
10.当进行定性判断时，此时式(1)中m
t2，空
为未知数，先设定安全壳标准泄漏率阈值，
再从标准泄漏率阈值计算出从t1时刻到t2时刻δt时间段内的泄漏量m
leak，空
；然后计算δt时间段内安全壳内湿空气压力上升值δp
′
，并与安全壳监测压升δp进行对比；当δp
′
大于δp，说明安全壳内湿空气泄漏率超过限值；当δp
′
小于δp，说明安全壳标准泄漏率未超限；
11.δp
′
＝p
′
t2-p
t1
ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
(2)
12.δp＝p
t2-p
t1
ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
(3)
13.式中，p
t2
为t2时刻安全壳内监测的湿空气压力；p
t1
为t1时刻安全壳内监测的湿空气压力；
14.p
′
t2
＝p
′
t2，空
+p
t2，水
ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
(4)
15.式中，p
′
t2，空
为t2时刻安全壳内干空气分压力；p
t2，水
为t2时刻安全壳内水蒸气的分压力；
16.当进行定量判断时，设定δp
′
与δp相等，即t2时刻安全壳内湿空气的计算压力p
′
t2
＝p
t2
，此时式(1)中泄漏量m
leak，空
为未知数；根据基于组分叠加的安全壳内气体的质量平衡关系和理想气体状态方程，求取δt小时时间段内安全壳内干空气的泄漏量m
leak，空
，并根据m
leak，空
计算δt小时时间段内安全壳泄漏率；
17.将m
leak，空
进行湿度修正和标准体积转换，设定安全壳内气压与大气压差值为60mbar时为标准，求取标准泄漏率q
l60
；
[0018][0019]
式中：rg为干空气气体常数；tn为标准温度；pn为标准压力；p
con
为δt小时时间段内安全壳内平均压力；p
atm
为δt小时时间段内平均大气压。
[0020]
当5nm3/h＜q
l60
＜10nm3/h时，应查并消除泄漏的原因，否则14天内机组开始向mcs模式后撤。当10nm3/h＜q
l60
时，3天内机组开始向mcs模式后撤，如果这个泄漏可以通过自动隔离或通过一个手动阀门隔离，那么机组开始后撤的时间延长为14天。
[0021]
所述的δt小时时间段δt可以是任意时长，一般取6至48小时，指机组启动阶段气闸门双关前的时间；设定安全壳标准泄漏率为5n
·
m3/h。
[0022]
所述的基于组分叠加的安全壳内气体的质量平衡关系是指任意时刻安全壳和气闸间内湿空气质量等于干空气质量和水蒸气质量之和。
[0023]
所述的δt小时时间段内安全壳与气闸间交换干空气质量的计算方法是基于理想气体状态方程，通过内气闸门开关前后的气闸间内的监测数据计算安全壳与气闸间的气体交换质量；
[0024][0025][0026]
式中：m
开门前，空，i
和m
开门后，空，i
分别为第i次开门前后气闸间内干空气的质量；p
开门前，空，i
和p
开门后，空，i
分别为第i次开门前后气闸间内干空气的压力，通过开门前后的相对湿度和温度计算获得；t
开门前，i
和t
开门后，i
分别为第i次开门前后气闸气闸间内湿空气的温度。
[0027]
所述的从标准泄漏率阈值计算出从t1时刻到t2时刻δt时间段内的泄漏量m
leak，空
，
具体计算方法是当q
l60
已知时，根据式(5)获得m
leak，空
的表达式：
[0028][0029]
根据质量守恒关系，任意时刻安全壳和气闸间内湿空气质量m，等于安全壳内干空气质量m
空
和水蒸气质量m
水
之和，表达式如下：
[0030]
m＝m
空
+m
水
ꢀꢀꢀꢀꢀꢀ
(9)
[0031]
根据理想气体状态方程，任意时刻安全壳内干空气质量m
t1，空
可表示为：
[0032][0033]
式中：p
空
为任意时刻安全壳内干空气分压力；v
con
为安全壳有效容积；rg为干空气气体常数；t为任意时刻安全壳内空气平均温度。
[0034]
根据道尔顿分压定律可知，任意时刻安全壳和气闸间内湿空气监测压力p等于安全壳内干空气分压力p
空
和水蒸气分压力p
水
之和，表达式如下：
[0035]
p＝p
空
+p
水
ꢀꢀꢀꢀ
(11)
[0036]
根据理想气体状态方程，所述安全壳内干空气分压力在t2时刻的计算公式为：
[0037][0038]
式中：t
t2
为t2时刻安全壳内的平均空气温度；v
con
为安全壳的有效容积。
[0039]
进一步的，任意时刻安全壳和气闸间内水蒸气分压力p
水
可表示为当前温度下水蒸气的饱和压力p
sat
(t)与当前相对湿度的乘积：
[0040]
p
水
＝p
sat
(t)
×
rh
ꢀꢀꢀꢀꢀꢀꢀ
(13)
[0041]
进一步的，安全壳内水蒸气的分压力p
t2，水
可表示为t2时刻温度下水蒸气的饱和压力p
sat
(t
t2
)与t2时刻安全壳内湿空气的相对湿度的乘积：
[0042]
p
t2，水
＝p
sat
(t
t2
)
×
rh
t2
ꢀꢀꢀꢀꢀꢀꢀ
(14)
[0043]
某一时刻的监测压力，则是对各测点的测量值取平均，作为某一时刻的监测压力。由于sar系统动作而注入的干空气质量m
sar
对安全壳内气体质量变化影响较大，所以需2秒采集一次。所述的安全壳内某一刻平均温度、相对湿度均为n个测点在某一时刻所获取的体积加权平均值。
[0044]
平均温度t
avg
的计算方法为：
[0045][0046]
式中：v
l
为安全壳的自由容积；vi为测点的有效测量体积；ti为某一时刻第i个测点监测的温度值。
[0047]
平均相对湿度h
avg
的计算方法为：
[0048][0049]
式中：rhi为某一时刻第i个测点监测的相对湿度值。
[0050]
本发明的有益效果是：在核电机组启动阶段的短时间内，根据安全壳和气闸间的短期(6-48小时)监测数据，一方面，通过假设安全壳标准泄漏率阈值为5n
·
m3/h时，可以在较短时间内实现对判断安全壳泄漏情况的判断；另一方面，假设计算压升等于监测压升，可
以根据定性判断法反推出短期内安全壳的泄漏量，进而定量计算短期安全壳的标准泄漏率。本发明可以通过安全壳和气闸间的短期监测数据，定性和定量的判断安全壳内空气泄漏情况，为判断安全壳的密封性提供支撑，为机组安全运行提供保障。
附图说明
[0051]
图1为安全壳内气体质量平衡关系示意图。
[0052]
图2为基于安全壳压升监测的安全壳泄漏快速判断方法流程图。
[0053]
图3为基于安全壳压升监测的安全壳标准泄漏率计算流程图。
具体实施方式
[0054]
为了对本发明的技术特征、目的和效果有更加清楚的理解，现以某cpr1000机组为例，根据其气闸门双关前24小时启机阶段监测数据，对一种基于安全壳压升监测的安全壳泄漏快速判断方法的判断过程进行展示，详细说明本发明的具体实施方式。
[0055]
一方面，提供一种基于安全壳压升监测的安全壳泄漏定性判断方法；另一方面，提供一种基于安全壳压升监测的安全壳泄漏率定量计算方法。
[0056]
一方面，根据安全壳和气闸间24小时监测数据，基于组分叠加的安全壳内气体的质量平衡关系、道尔顿分压定律的压力平衡关系和理想气体状态方程，计算0h到24h时假设安全壳标准泄漏率为5n
·
m3/h时的安全壳内的计算压升，并与安全壳监测压升进行对比，若计算压升大于监测压升，则说明安全壳泄漏率超过5n
·
m3/h的限值；若计算压升小于监测压升，则说明安全壳泄漏率处于正常水平。
[0057]
根据质量守恒关系，任意时刻安全壳内湿空气质量m，等于安全壳内干空气质量m
空
和水蒸气质量m
水
之和，表达式如下；
[0058]
m＝m
空
+m
水
[0059]
根据道尔顿分压定律可知，任意时刻安全壳内湿空气监测压力p等于安全壳内干空气分压力p
空
和水蒸气分压力p
水
之和，表达式如下；
[0060]
p＝p
空
+p
水
[0061]
进一步的，安全壳内干空气分压力p
0h，水
可表示为初始温度下水蒸气的饱和压力p
sat
(t
0h
)与初始时刻安全壳内湿空气的相对湿度的乘积；
[0062]
p
0h，水
＝p
sat
(t
0h
)
×
rh
0h
＝3456.97
×
46.53％＝1609pa
[0063]
进一步的，安全壳内干空气分压力p
0h，空
可表示为；
[0064]
p
0h，空
＝p
0h-p
0h，水
＝102015-1609＝100406pa
[0065]
式中：v
con
为安全壳有效容积；rg为干空气气体常数。
[0066]
根据理想气体状态方程m
0h，空
表示为：
[0067][0068]
根据理想气体状态方程，可获得24小时安全壳与气闸间交换干空气质量δm
闸，空
的表达式：
[0069][0070][0071]
根据气闸间监测数据，24小时内气闸间内门共开启18次，即i＝18，求得δm
闸，空
＝-30157.49g。
[0072]
当进行定性判断时，先设定安全壳标准泄漏率q
l60
＝5n
·
m3/h，由式(8)可知，当q
160
＝5n
·
m3/h时，对应的24小时安全壳泄漏的干空气质量m
leak，空
＝79670.62g，根据sar监测数据，24小时sar系统注入干空气质量m
sar
＝212164.74g。
[0073]
根据安全壳内气体质量平衡关系可知，24小时后安全壳内计算干空气质量m
24h，空
等于初始时刻安全壳内干空气m
0h，空
、m
sar
、δm
闸，空
和m
leak，空
之和。
[0074]m24h，空
＝m
0h，空
+m
sar
+δm
闸，空
+m
leak，空
[0075]
＝57765511.03+212164.74-30158.49-79670.62
[0076]
＝57867847.67g
[0077]
根据质量守恒关系，初始时刻安全壳内湿空气质量m
0h
等于初始时刻安全壳内干空气质量m
0h，空
和水蒸气质量m
0h，水
之和，表达式如下：
[0078]m0h
＝m
0h，空
+m
0h，水
[0079]
进一步的，安全壳内水蒸气分压力p
24h，水
可表示为初始温度下水蒸气的饱和压力p
sat
(t
24h
)与初始时刻安全壳内湿空气的相对湿度的乘积：
[0080]
p
24h，水
＝p
sat
(t
24h
)
×
rh
24h
＝4886.98
×
36.06％＝1762.4pa
[0081]
进一步的，根据理想气体状态方程，24小时后安全壳内计算干空气分压力p
′
24h，空
可表示为：
[0082][0083]
根据道尔顿分压定律可知，24小时后安全壳内湿空气计算压力p
24h
等于24小时后安全壳内干空气分压力p
′
24h，空
和水蒸气分压力p
24h，水
之和，表达式如下：
[0084]
p
′
24h
＝p
′
24h，空
+p
24h，水
＝102598.24+1762.4＝104360.64pa
[0085]
24小时后安全壳内湿空气计算压升可表示为：
[0086]
δp
′
＝p
′
24h-p
0h
＝104360.64-102015＝2345.64pa
[0087]
24小时后安全壳内湿空气监测压升可表示为：
[0088]
δp＝p
24h-p
0h
＝104468-102015＝2453pa
[0089]
式中：p
24h
为24小时后安全壳内湿空气的监测的空气压力。
[0090]
从压升的计算结果可以看出，计算压升δp
′
小于监测压升δp，则说明安全壳内湿空气泄漏率未超过5n
·
m3/h的限值，安全壳泄漏率处于正常水平。
[0091]
第二方面，提供一种基于安全壳压升监测的安全壳泄漏率定量计算方法，其特征在于，假设δp
′
与δp相等，即p
′
24h
＝p
24h
，根据基于组分叠加的安全壳内气体的质量平衡关系和理想气体状态方程，求取24小时安全壳干空气的泄漏量m
leak，空
，并根据m
leak，空
计算24小
时安全壳标准泄漏率。
[0092]
根据安全壳内气体质量平衡关系可知，24小时安全壳泄漏干空气质量m
leak，空
等于24小时后安全壳内干空气质量m
24h，空
分别减去初始时刻安全壳内干空气m
0h，空
、24小时sar系统注入干空气质量m
sar
、24小时安全壳与气闸间交换干空气质量δm
闸，空
，表达式如下：
[0093]mleak，空
＝m
24h，空-m
0h，空-m
sar-δm
闸，空
[0094]
其中，根据理想气体状态方程可获得m
24h，空
表达式：
[0095][0096]
进一步的，24小时安全壳泄漏干空气质量m
leak，空
可表示为：
[0097]mleak，空
＝57928400.5g-57765511.03-212164.74+30158.49
[0098]
＝-18637.13g
[0099]
根据24小时后安全壳内湿空气的分压力p
24h，水
对m
leak，空
进行湿度修正：
[0100][0101]
根据理想气体状态方程将m
leak
进行标注体积转换：
[0102][0103]
式中：tn为标准温度；pn为标准压力。
[0104]
以安全壳内气压与大气压差值为60mbar时为标准，求取标准泄漏率q
160
：
[0105][0106]
式中：p
con
为24小时安全壳内平均空气压力；p
atm
为24小时平均大气压力。
[0107]
当5nm3/h＜q
l60
＜10nm3/h时，应查并消除泄漏的原因，否则14天内机组开始向mcs模式后撤。当10nm3/h＜q
l60
时，3天内机组开始向mcs模式后撤，如果这个泄漏可以通过自动隔离或通过一个手动阀门隔离，那么机组开始后撤的时间延长为14天。
[0108]
所述的某一时刻的监测压力，则是对各测点的测量值取平均，作为某一时刻的监测压力。由于sar系统动作而注入的干空气质量m
sar
对安全壳内气体质量变化影响较大，所以需2秒采集一次。所述的安全壳内某一刻平均温度、相对湿度均为n个测点在某一时刻所获取的体积加权平均值。
[0109]
平均温度t
avg
的计算方法为；
[0110][0111]
式中：v
l
为安全壳的自由容积；vi为测点的有效测量体积；ti为某一时刻第i个测点监测的温度值。
[0112]
平均相对湿度h
avg
的计算方法为；
[0113][0114]
式中：rhi为某一时刻第i个测点监测的相对湿度值。
[0115]
综上所述，本发明实施例通过核电机组安全壳和气闸间24小时监测数据，一方面，通过假设安全壳标准泄漏率为5n
·
m3/h时，基于组分叠加的安全壳内气体的质量平衡关系、道尔顿分压定律的压力平衡关系和理想气体状态方程，定性判断了安全壳泄漏率是否超限；另一方面，假设计算压升等于监测压升，根据基于组分叠加的安全壳内气体的质量平衡关系以及理想气体状态方程，定量计算除安全壳的标准泄漏率。

技术特征：

1.一种基于安全壳压升监测的安全壳泄漏快速判断方法，其特征在于，安全壳泄漏快速判断方法分为定性判断法和定量判断法；安全壳内湿空气组分满足质量平衡关系、道尔顿分压定律和理想气体状态方程；安全壳内从t1时刻到t2时刻，采用以下干空气质量平衡关系：m
t2，空
＝m
t1，空
+m
sar
+δm
闸，空
+m
leak，空
ꢀꢀꢀꢀ
(1)式中，m
t2，空
为t2时刻安全壳内干空气质量；m
t1，空
为t1时刻安全壳内干空气质量；m
sar
为从t1时刻到t2时刻δt时间段内sar系统注入安全壳的干空气质量；δm
闸，空
为δt时间段内安全壳与气闸间交换干空气质量，从壳内到气闸间的质量为负；m
leak，空
是δt时间段内安全壳的干空气泄漏量，从壳内泄漏到壳外的质量为负；当进行定性判断时，此时式(1)中m
t2，空
为未知数，先设定安全壳标准泄漏率阈值，再从标准泄漏率阈值计算出从t1时刻到t2时刻δt时间段内的泄漏量m
leak，空
；然后计算δt时间段内安全壳内湿空气压力上升值δp
′
，并与安全壳监测压升δp进行对比；当δp
′
大于δp，说明安全壳内湿空气泄漏率超过限值；当δp
′
小于δp，说明安全壳标准泄漏率未超限；且δp
′
＝p
′
t2-p
t1
ꢀꢀꢀꢀ
(2)δp＝p
t2-p
t1
ꢀꢀꢀꢀ
(3)式中，p
′
t2
为t2时刻安全壳内的湿空气压力计算值，p
t2
为t2时刻安全壳内监测的湿空气压力；p
t1
为t1时刻安全壳内监测的湿空气压力；且p
′
t2
＝p
′
t2，空
+p
t2，水
ꢀꢀꢀꢀ
(4)式中，p
′
t2，空
为t2时刻安全壳内干空气分压力；p
t2，水
为t2时刻安全壳内水蒸气的分压力；当进行定量判断时，设定δp
′
与δp相等，即t2时刻安全壳内湿空气的计算压力p
′
t2
＝p
t2
，此时式(1)中泄漏量m
leak，空
为未知数；根据基于组分叠加的安全壳内气体的质量平衡关系和理想气体状态方程，求取δt小时时间段内安全壳内干空气的泄漏量m
leak，空
，并根据m
leak，空
计算δt小时时间段内安全壳泄漏率；将m
leak，空
进行湿度修正和标准体积转换，设定安全壳内气压与大气压差值为60mbar时为标准，求取标准泄漏率q
l60
；式中：r
g
为干空气气体常数；t
n
为标准温度；p
n
为标准压力；p
con
为δt小时时间段内安全壳内平均压力；p
atm
为δt小时时间段内平均大气压。2.根据权利要求1中所述的基于安全壳压升监测的安全壳泄漏快速判断方法，其特征在于，所述的δt小时时间段取6至48小时，指机组启动阶段气闸门双关前的时间。3.根据权利要求1中所述的基于安全壳压升监测的安全壳泄漏快速判断方法，其特征在于，所述的基于组分叠加的安全壳内气体的质量平衡关系是指任意时刻安全壳和气闸间内湿空气质量等于干空气质量和水蒸气质量之和。4.根据权利要求3中所述的基于安全壳压升监测的安全壳泄漏快速判断方法，其特征在于，基于组分叠加的安全壳内气体的质量平衡关系，δt时间段内安全壳与气闸间交换干空气质量的计算方法是基于理想气体状态方程，通过内气闸门开关前后的气闸间内的监测数据计算安全壳与气闸间的气体交换质量；
式中：m
开门前，空，i
和m
开门后，空，i
分别为第i次开门前后气闸间内干空气的质量；p
开门前，空，i
和p
开门后，空，i
分别为第i次开门前后气闸间内干空气的压力，通过开门前后的相对湿度和温度计算获得；t
开门前，i
和t
开门后，i
分别为第i次开门前后气闸间内湿空气的温度。5.根据权利要求1中所述的基于安全壳压升监测的安全壳泄漏快速判断方法，其特征在于，所述的从标准泄漏率阈值计算出从t1时刻到t2时刻δt时间段内的泄漏量m
leak，空
，具体计算方法是当q
l60
已知时，根据式(5)获得m
leak，空
的表达式：6.根据权利要求1中所述的基于安全壳压升监测的安全壳泄漏快速判断方法，其特征在于，所述的基于道尔顿分压定律的压力平衡关系，指任意时刻安全壳和气闸间内湿空气压力均表示为干空气分压力和水蒸气分压力之和。7.根据权利要求1中所述的基于安全壳压升监测的安全壳泄漏快速判断方法，其特征在于，所述安全壳内干空气分压力在t2时刻的计算公式为：式中：t
t2
为t2时刻安全壳内的平均空气温度；v
con
为安全壳的有效容积。

技术总结

本发明属于核电厂安全技术领域，提出一种基于安全壳压升监测的安全壳泄漏快速判断方法。当进行泄漏率是否超标的定性判断时，通过假设安全壳标准泄漏率，计算安全壳内湿空气压力上升值获得计算压升，并与安全壳监测压升进行对比。若计算压升大于监测压升，则说明安全壳内湿空气泄漏率超过限值；若计算压升小于监测压升，则说明安全壳标准泄漏率未超限。当需要定量计算泄漏率时，则假设计算压升等于监测压升，然后求取对应的安全壳泄漏率。本发明可以通过安全壳和气闸间内短期监测数据，定性和定量的判断安全壳内空气泄漏情况，为机组启动期间判断安全壳的密封性提供支撑，为机组安全运行提供保障。运行提供保障。运行提供保障。