本文采 用 MELCOR 程序,以某船用压水堆为研究对象,分析了冷却剂主 管道冷管段中破口失水事故,将所得结果与该船用堆的安全分析报告进行了比对,验证了 MELCOR 在计算设计基准事故时的准确性。并针对安注方式、停泵时间对事故 后果的影响进行了分析,对预测事故进程、缓解事故后果具有重要意义。<br/> In this paper, MELCOR code was applied to analyze the break loss of coolant accident occurring at the cold leg in the main pipe, taking the ship reactor for the object of the research. And the results were compared with the safety analysis report for the ship reactor, which proved the accuracy of the calculation for design basic accident (DBA) by MELCOR. The consequences of the accident was analyzed according to the safety injection manner and the pump stopping time, which is important for predicting the accident progress and alleviating the consequence of the accident.
船用核动力装置是以核裂变能作为产生推进动力的能源,由于其在航行和停靠岸过程中,需要经常改变功率[
本文以某船用压水堆为对象,采用MELCOR1.8.5程序对设计基准失水事故(DBLOCA)进行研究[
MELCOR程序是由桑地亚国家实验室(SNL)为美国核管会开发的严重事故分析程序,能模拟轻水反应堆严重事故进程的主要现象,计算放射性核素的释放及其后果。MELCOR程序是模拟严重事故全过程的一体化程序包,经过许多试验数据的验证和严重事故分析程序的相互验证,对事故进程的模拟具有很高的可信度,也是IAEA向世界推荐使用的严重事故分析程序之一,在国际上有数十个国家和地区的用户在使用这个程序。
MELCOR程序不仅能计算严重事故,也能进行设计基准事故的分析。为验证本程序对船用堆计算结果的准确性,首先我们进行了稳态运行计算,并与实际运行参数进行了比对,验证了稳态计算的准确性;然后针对设计基准事故,将相同事故序列下的MELCOR的计算结果与该船用堆的安全分析报告进行了比对,进一步验证本程序的对设计基准事故计算的准确性。在此基础上所进行的严重事故分析,才能确保合理、可信。
反应堆冷却剂系统由主冷却剂管道、主泵、蒸汽发生器一次侧、稳压器、安注箱、安注泵和堆芯等部件组成[
MELCOR采用COR程序包计算堆芯、下腔室内部结构的热工、水力响应,为了将堆芯的传热分布细化,以便于研究轴向和径向的热量产生及传递,将堆芯活性区及下腔室进行区域节点划分,堆芯及下腔室节点的划分如图1所示,为准确描述设计基准事故后果,将船用堆堆芯径向划分为7个同心圆环,轴向划分为14层,其中堆芯活性区部分分为12段,下腔室分为2段,包括下管板和下腔室。
DBLOCA事故发生时,安注水若不能及时弥补泄漏的水量,将会导致堆芯裸露,包壳温度迅速上升,锆合金包壳与过热水蒸汽发生显著的锆–水反应,如果化学反应过度,会导致包壳脆化破裂。失水事故分析时,包壳破损的限制准则[4,5]如下:
1) 峰值包壳表面温度不超过限值(1204℃);
2) 包壳总氧化率不超过包壳总厚度的17%;
3) 包壳与水或水蒸汽发生化学反应后的产氢量不超过全部产氢量的1%;
4) 堆芯几何形状不改变。
为验证计算结果的准确性,首先进行稳态计算,并将各参数与实际运行参数进行比对,完全吻合,证明程序稳态计算正确性。然后,在事故序列完全一致的情况下,将MELCOR的计算结果与安全分析报告
图1. 堆芯节块划分
和RELAP5的结果进行比对,符合较好。表明程序对设计基准事故计算的准确性。限于篇幅,将主要事故序列的数据进行对比,列于表1中。
计算中的主要输入参数及假设见表2。
1) 对破口位置在冷管段,破口当量直径为29.4%RM的中破口失水事故,采用冷端安注、热端安注方式,将两种方式的后果进行对比。图2~图4给出了破口流量M、堆芯水位H、包壳最高温度TH随时间的变化。
由图2~图4可以看出,冷端安注相对于热端安注,其破口流量较大,堆芯水位较低,包壳的最高温度要高出很多,其原因是:冷端安注时,部分安注水会通过破口直接流出,不能有效注入堆芯;其次,由于堆芯冷却剂开始汽化,导致冷却剂流量减少、滞留甚至倒流,从堆芯排出的蒸汽与下腔室继续蒸发产生的蒸汽一起,通过下降段向上流动,与冷端安注水形成逆向流,阻止冷端的安注水进入;另外,冷端安注水还有很大一部分被蒸汽夹带到破口,并不通过下降段而直接被带到破口流出。因此,堆芯中上部仍裸露在外面,未得到有效冷却,导致包壳温度一直上升;热端安注时,安注水由上至下流经堆芯,不会出现堆芯中上部裸露的情况,从而避免了其温度的继续上升[
通过计算分析得,冷管段29.4%RM破口采用冷端安注时,元件包壳出现破损;采用热瑞安注时,元件包壳未出现破损。热端安注的效果比冷端安注好。
2) 破口位置在冷管段,破口当量直径为14.9%RM的中破口失水事故,采用不同停泵时间,将其结果进行分析。图5~图7给出了破口流量M、堆芯水位H、包壳最高温度TH随时间的变化。
由图7可以看出,当主泵停止运行的整定值为50.7%P0 (P0为额定运行压力)时,所得包壳最高的温度要高于整定值为21.7%P0的情况。由于前者主泵过早停转,强迫循环终止,导致堆芯产生的热量难以及时排除,包壳温度上升。
计算结果表明,若主冷却剂压力降至50.7%P0时停止转动,会导致包壳破损,若压力降至21.7%P0主
表1. 主要事故序列
表2. 事故主要输入参数及假设
图2. 破口流量
图3. 堆芯水位
图4. 元件包壳温度变化
图5. 破口流量
图6. 堆芯水位
图7. 元件包壳温度变化
泵停止转动,包壳则不会破损。
1) 本文用MELCOR程序分析设计失水事故,所得结果与RELAP5的结果一致,验证了MELCOR程序在计算设计基准事故时的准确性。
2) 对于冷管段29.4%RM破口的中破口失水事故,采用热瑞安注时,元件包壳未出现破损,采用冷端安注包壳破损。冷管段破口热端安主可缓解事故后果。
3) 对于14.9%RM中破口失水事故,主泵在主回路压力降至21.7%P0时停转,包壳破损;主泵在主回路压力降至50.7%P0时停转,包壳不破损。事故后主泵保持低速运行可缓解事故后果。
[
[
[
[
[
[