1. 引言

主蒸汽管道大破口是“一次侧热输出增加”这类事故中最严重的反应堆冷却剂降温瞬态，而发生在安全壳外的主蒸汽管道破裂事故(MSLB)将旁通安全壳这一屏障，放射性物质可以通过破口直接释放到环境中。MSLB对核电厂安全具有重要影响，在环境影响评价及安全分析报告中均需进行相应评价。

针对最极限的安全壳外主蒸汽管道双端断裂事故，本文结合国内某AP1000滨海厂址具体特征开展了事故放射性后果评价，全面评价了该事故对主控制室工作人员的剂量影响。同时，本文深入研究确定了对剂量起主导作用的放射源、剂量主要产生时段，并且识别出了剂量影响关键参数，为主控制室设计优化指明了方向。

2. 事故描述

MSLB事故将引起二回路蒸汽从破口释放，起初在短时间内蒸汽流量快速增加，而后，随着蒸汽压力下降，蒸汽流量逐渐减小。由于一次侧排出热量过多导致反应堆冷却剂系统降温降压，反应堆紧急停堆。事故过程中，破损主蒸汽管道及其启动给水将被隔离，反应堆由完好回路或者非能动余热排出热交换器冷却。事故放射性后果分析中保守的假设破损回路的蒸汽发生器(SG)被蒸干，从一回路泄漏到该破损回路SG的全部放射性通过破口直接排放入环境。完好回路SG不会裸露，大部分放射性会滞留在二次侧。为了简化分析，也保守假设其初始放射性全部释放。

对于主蒸汽管道破裂事故剂量影响评价，通常考虑两种独立工况：

· 工况1为MSLB叠加事故前碘尖峰释放；

· 工况2为MSLB引起事故并发碘尖峰释放。

放射性产生途径可归纳为如下三方面：

· 通过SG传热管泄漏到二次侧的RCS初始碘、碱金属、惰性气体；

· 二回路初始存在的碘、碱金属；

· 事故前碘尖峰释放(适用工况1)或事故并发碘尖峰释放(适用工况2)。

3. 源项分析模型

3.1. RCS初始放射性释放

事故发生后，RCS中放射性核素将伴随事故进程而有所损耗，主要包括核素自身衰变以及向二次侧泄漏损耗，则其活度为：

(1)

式中，C_i(t)为事故发生后t时刻RCS中核素i的活度浓度，Bq/g；C_0,i为事故发生初始时刻RCS中核素i的活度浓度，Bq/g；λ_i= λ_d,i+ 2F_l/m₁，λ_d,i为核素i的衰变常数，h⁻¹；F_l为每台SG一次侧向二次侧泄漏率，g/h；m₁为RCS冷却剂质量，g。

通过蒸汽发生器由一次侧泄漏到二次侧，并进一步释放到环境中的核素活度微分方程如下：

(2)

求解上述微分方程可以计算得到下式：

(3)

式中，Q_1,i(t)为RCS中初始核素i在t时刻释放到环境的活度，Bq；P_c,i为完好回路SG汽水分配系数(破损回路保守忽略稀释及分配效应)，只有碘及碱金属需要考虑该汽水分配效应，惰性气体假设不滞留而直接释放；其余参数物理意义同上。

3.2. 二回路初始放射性释放

二回路冷却剂中主要包括初始存在的碘及碱金属，伴随事故进程释放到环境的源项如下式：

(4)

式中，Q_2,i(t)为二回路中初始核素i在t时刻释放到环境的活度，Bq；B_0,i为事故发生初始时刻二回路中核素i的活度浓度，Bq/g；m₂为每台SG二回路的冷却剂质量，g；t₁为破损回路SG喷放持续时间，h；t₂为完好回路SG释放持续时间，h。

3.3. 碘尖峰释放

1) 事故前碘尖峰释放

事故发生前的功率、压力等瞬态作用，可以使存在破损的燃料棒向反应堆冷却剂的碘释放速率短时间内快速增加，而产生一种尖峰释放现象，即事故前碘尖峰释放。通常，假设事故前碘尖峰导致一回路冷却剂碘活度浓度达到公式(1) C_0,i的60倍 [1] 。

2) 事故并发碘尖峰释放

与事故前碘尖峰有所不同，事故并发碘尖峰是由于事故发生过程的瞬态作用而产生的尖峰释放现象。通过并发碘尖峰释放到一回路冷却剂内核素活度微分方程为：

(5)

式中，A_i(t)为t时刻通过并发碘尖峰释放到RCS内的碘核素i活度，Bq；R_i为并发碘尖峰向RCS释放率，Bq/h，其余参数物理意义同上。

求解上述微分方程得到下式：

(6)

假设泄漏到破损回路SG的碘不滞留全部直接释放，而泄漏到完好回路SG的碘需要考虑汽水分配效应，则释放到环境碘活度微分方程为：

(7)

求解上述微分方程得到下式：

(8)

式中，Q_3,i(t)为事故后t时刻通过并发碘尖峰释放到环境的碘核素i活度，Bq；其余参数物理意义同上。

4. 主控制室剂量评价

4.1. 通风模式

主控制室剂量计算时，分别考虑两种通风模式：应急可居留系统(VES)通风模式和非放射性通风系统(VBS)新风过滤模式。

VES模式分析中，假设通风系统初始工作在VBS正常运行模式。当惰性气体探测器发出“高”报信号后，保守假设非安全有关的VBS新风过滤模式未被触发投入。主控制室VBS正常运行模式将持续运行，直至放射性粒子探测器或碘探测器发出“高”报信号后，安全有关的VES模式被触发投入使用。本文分析了该VES模式情况下主控制室剂量。

同时，按照纵深防御的原则，当惰性气体探测器发出“高”报后，非安全有关的VBS新风过滤模式也有可能被触发投入，以保护主控制室人员。因此，本文也评估了从VBS正常运行模式切换到新风过滤模式时主控制室的剂量(此分析不需要切换到VES模式)。

4.2. 评价模型

1) MCR内核素活度

VBS新风过滤模式情况下，主控制室内放射性核素活度可由下式计算：

(9)

式中，A_c,i为事故后各时刻主控制室内核素i活度，Bq；Q_i为单位时间内向环境释放的核素i活度，Bq/s；(χ/Q)_k为k时段大气弥散因子，s/m³；F₁为新风流量，m³/h；F₂为内循环流量，m³/h；F₃为非过滤渗入率，m³/h；F₄为主控制室向环境排风流量，m³/h；η_i为过滤器对核素i过滤效率；λ_d,i为核素i衰变常数，h⁻¹；V为主控制室容积，m³。

上式中，放射性通过新风及内渗漏进入主控制室，在主控制室内考虑衰减、内循环过滤、排风等衰减效应。通过求解该微分方程可以计算主控制室在事故后各时刻放射性性核素活度，进一步对时间进行积分可得到累积活度，除以主控制室容积即可得到时间积分浓度。

对于主控制室VES模式，已经切断外界送风，故需将公式(9)中的新风流量项去掉。

2) MCR内剂量

主控制室内工作人员有效剂量包括内照射有效剂量及外照射有效剂量两部分总和。主控制室是一个有限大的空间，相比于半无限大烟羽外照射模型，主控制室工作人员外照射有效剂量需要进行几何修正。同时，剂量评价需要考虑事故后30天内接受剂量，时间跨度比较长，需要引入居留时间份额，即居留因子。

有效剂量计算公式如下：

(10)

式中，D_E为有效剂量，Sv；D_EDE为外照射有效剂量，Sv；D_CEDE为内照射有效剂量，Sv；GF为主控制室几何修正因子 [1] ，GF = 351.6/V^0.338；DCF_1,i为核素i外照射有效剂量转换因子 [2] ，Sv·m³/Bq·s；DCF_2,i为核素i内照射有效剂量转换因子 [3] ，Sv/Bq；(IAR)_i,j为在时间段j内主控制室累积的核素i活度，Bq·s/m³；BR_j为时间段j内工作人员的呼吸率 [1] ，m³/s；O_j为时间段j内工作人员在主控制室内的居留因子 [1] 。

4.3. 剂量结果

表1计算给出了两种工况下主控制室工作人员接受的剂量结果，工况1发生30天后VES模式及VBS新风过滤模式下有效剂量分别为15.38 mSv、4.72 mSv，甲状腺剂量分别为288.49 mSv、135.11 mSv；工况2发生30天后VES模式及VBS新风过滤模式下有效剂量分别为16.39 mSv、14.62 mSv，甲状腺剂量分别为322.22 mSv、453.94 mSv，剂量结果均满足HAD 002/01-2010 [4] 规定的剂量限值要求(有效剂量限值50 mSv、甲状腺剂量限值500 mSv)。

5. 敏感性分析

鉴于AP1000核电厂主控制室VBS通风模式为非安全级，仅作为纵深防御手段，因此本文主要以VES模式为基础开展关键参数敏感性分析。

5.1. 主导放射源

MSLB事故后，一回路、二回路冷却剂中初始存在的放射性以及碘尖峰释放的放射性将通过大气释放阀或SG破口释放到环境。为了定量评估上述放射源的相对贡献、确定主导放射源，本文针对MSLB引起并发碘尖峰工况，定量的评估了各类放射源对主控制室有效剂量的贡献(针对VES模式情况)，具体剂量结果见图1。结果表明，二回路初始存在的碘及碱金属剂量贡献最大(分别为8.17 mSv及6.86 mSv)，并发碘尖峰释放的碘也具有一定贡献(1.34 mSv)，一回路冷却剂初始存在放射性贡献最小(小于0.01 mSv)。

5.2. 各时间段接受剂量

事故发生后初始0.052 h，主控制室工作在VBS正常通风模式下；之后启动VES模式，由灌装压缩空气提供送风，并且由非能动空气过滤管线提供内循环过滤，整个过程持续到72 h；空气耗光后打开房门由辅助风机进行送风，一直持续到厂外供气恢复(事故后168 h)。为了定量的评估事故后主控制室内剂量变化情况，图2给出了各时段所获得剂量。结果表明：事故后主控制室内剂量主要来自于0.052~4 h时间段，主要是VBS正常通风引入的放射性产生的剂量；72~168 h，虽然主控制室敞开门，但事故后放射性释放已经结束，该时段剂量贡献非常小。

表1. 主控制室剂量

5.3. 关键参数敏感性分析

敏感性分析是指从定量分析的角度研究有关因素发生某种变化对某一个或一组关键指标影响程度的一种不确定分析技术。敏感性因素一般可选择主要参数进行分析。若某参数的小幅度变化能导致关键指标的较大变化，则称此参数为敏感性因素，反之则称其为非敏感性因素。

本文用敏感性系数来评价主要参数对主控制室有效剂量的影响程度。敏感性系数计算公式如下：

式中，S_DF为敏感性系数；为有效剂量结果的变化率；为不确定因素的变化率。

当S_DF> 0时，表示剂量结果与不确定性因素同方向变化；当S_DF< 0时，表示剂量结果与不确定性因素反方向变化。越大，表明剂量结果D_E对于不确定性因素F越敏感；反之，则不敏感。

针对MSLB事故破损SG喷放时间、VES模式投用时间、MCR内渗透风量、内循环风量、辅助风机风量等，本文计算给出了主控制室有效剂量对这些参数的敏感性系数，具体见表2。上述参数基准值依次为：600 s、188 s、25.5 m³/h、1020 m³/h、2890 m³/h；敏感性系数针对变化±33%、±67%、±100%情况进行评估。结果表明：

l 喷放时间直接决定二回路初始放射性释放速率，喷放时间越短，主控制室剂量越大，并且越敏感；

l VES模式投用时间决定进入主控制室内放射性量，投用时间越晚，主控制室剂量越大，并且敏感性系数几乎为常值；

l 内渗透率对主控制室剂量影响不大，为不敏感参数；

l 内循环风量直接决定对进入主控制室内放射性的循环去除能力，内循环风量越小，主控制室剂量越大，并且越敏感；

l MSLB事故在72 h内放射性释放结束，剂量主要是在辅助风机投用之前产生，辅助风机风量对主控制室剂量几乎无影响。

按照传统观点，主控制室VES模式启动情况下，内渗透是放射性进入主控制室的关键途径，其取值应对主控制室剂量产生较大影响，这与内渗透率敏感性分析结论存在一定差异。为了解释该现象，本文进一步研究了不同内渗透率情况下主控制室内DE.I-131放射性水平随时间的动态变化，具体如图3所示。

表2. 关键参数敏感性系数

结果表明，事故后初始4 h，内渗透率对主控室内放射性水平无明显影响。这是因为VES模式启用前送入MCR的放射性远大于内渗透漏入的放射性，即使VBS正常通风模式被隔离，不再有放射性送入，但滞留的放射性也占据主导。但随着VES内循环过滤器对MCR内放射性的循环过滤，这部分滞留的放射性呈逐渐减少趋势。直到4 h后，内渗透进入MCR的放射性逐渐超过了原来滞留的放射性，则此时刻之后内渗透率越大MCR内放射性越大。而图2分析结果表明，VES模式情况下MCR剂量基本上在前4 h内产生，而内渗透率对MCR剂量产生影响是在事故发生4 h后，故最终MCR剂量对其不敏感。

上述分析表明内循环风量对MCR有效剂量影响较大，为了进一步反映整个事故持续过程中的具体影响，图4给出了不同内循环风量情况下MCR内放射性动态变化。结果表明，内循环机制启用后(0.052 h)，将显著降低MCR内放射性水平；同时内循环风量越大，放射性过滤效果越明显。

6. 结论

基于AP1000滨海厂址，本文开展了MSLB事故对主控制室剂量影响研究，并进行了全面的敏感性分析(针对叠加事故并发碘尖峰工况、主控制室处于VES通风模式)，得到结论如下：

1) MSLB事故后主控制室开启任何一种通风过滤模式，均能够对工作人员提供一定的防护，有效剂量均可以满足HAD 002/01-2010规定的50 mSv剂量限值要求。

2) 对于MSLB事故，源项包括来自于一回路、二回路中初始存在的放射性，以及碘尖峰释放的放射性，主控制室剂量主要取决于二回路放射性的释放。

3) 对于MSLB事故，主控制室剂量主要是在前4 h时段产生，72 h后门敞开阶段将无显著剂量贡献。

4) 由敏感性分析可知，主控制剂量对二回路冷却剂喷放时间、内循环风量、VES模式投用时间比较敏感，而内渗透率及辅助风机风量为不敏感参数。

参考文献