堆芯熔化过程中容器内熔池的热聚集效应.docxVIP

堆芯熔化过程中容器内熔池的热聚集效应 1 堆芯熔融物塌至下封头形成熔池的过程 在严重的堆芯硬化后，通过外部冷却将融化物停留在容器中（ivr）是一种重要的核电站稀释措施。目前在役的核电站和先进堆型都采用了IVR缓解措施并对其进行了有效性评价。 与堆芯熔化及熔融物坍塌至下封头形成熔池的过程相关的现象对下腔室熔池的成分和结构有着十分重要的影响;而熔池的特性又影响容器壁的热传递和下封头的完整性。因此,在核电站IVR措施有效性评价中,有必要研究堆芯熔化及其塌落至下封头形成熔池的具体过程。本文主要针对岭澳核电站的CPR1000堆型,采用一体化严重事故分析程序对选定的严重事故序列进行计算,以研究严重事故下堆芯的熔化、塌落并在下封头形成熔池的整个过程;并通过有、无容器外冷却模型的对比计算评价CPR1000堆型IVR措施的有效性。 2 熔池的结构和特点 2.1 层熔池模型 由于熔池在达到最终稳定分层结构之前所要经历的中间过程仍存在较大的不确定性,不少专家提出了不同的熔池分层结构。本文采用目前最具有代表性的2层熔池模型(图1),这种结构具有一定的保守性。 熔池底层为密度较大的氧化池,包含UO2、Zr O2等氧化物和内热源。此外,氧化池与温度较低的下封头和上部金属层接触,会冷凝成一层薄壳。由于薄壳顶部覆盖的金属层温度与下封头容器壁温度不同,故氧化池上下薄壳厚度并不一致。 熔池上层由密度相对较小的金属层构成,包含熔化的不锈钢和未被氧化的锆等金属物。由于金属层质量与氧化池质量相比往往很小,故其厚度较薄。 2.2 热聚集效应的形成 根据文献,无论采用何种熔池结构,当堆芯熔化并塌落至下封头形成分层熔池后,下封头壁面受热的热流密度沿其轴向是非均匀分布的;熔池上层的金属层由于质量小、厚度薄,造成其热流密度很大,形成热聚集效应。若此处的热流密度超过临界热流密度(CHF)值,会对压力容器的完整性构成巨大威胁。 3 计算程序模型 3.1 下封头节点划分 本文的研究对象为岭澳核电站CPR1000堆型。反应堆冷却剂系统被整合模拟为2个环路,其中“破损环路”模拟带稳压器的环路,“完整环路”综合模拟CPR1000的另外2个环路(图2)。 反应堆堆芯沿径向被划分为7个环;沿轴向被划分为13个节点,其中10个为活性区节点,1个为堆芯上部非活性区节点,另外2个为堆芯下部非活性区节点。堆芯初始功率为2895 MW,停堆后堆芯衰变功率采用美国国家标准学会(ANSI)标准公式计算。 下封头的节点划分如图3所示。在对称的半壳上,沿轴向上平均划分为5个节点,在其径向方向上也等分地划分为5个节点。整个下封头中分布58个贯穿件,并在下封头的可能失效模式中考虑贯穿件失效。 3.2 锆包壳温度和厚度 本文重点关注锆的氧化,因其大小直接影响熔池金属层质量(经计算发现,不锈钢的氧化可以忽略)。锆的氧化主要是通过锆-水反应,其反应公式为: 式中,QOX为每单位锆发生反应后产生的热量。 不同温度下锆的氧化速率由Cathcart和Baker-Just关系式分别给出。 当400KT≤1850K时,氧化层厚度变化率Δxo为: 当T1875K时,Δxo为: 当1850KT≤1875K时,Δxo由式(2)、式(3)插值给出: 上述各式中,R为气体常数;ρZr为锆包壳密度;T为锆包壳温度。 当堆芯温度超过2000 K时,UO2将会溶入锆中而形成U-Zr-O混合物。 3.3 熔池的形成过程 本研究所采用的计算程序提供了一种具有中间分层过程的分层结构,认为熔池在达到最终包络状态之前存在一个过渡结构,这种过渡熔池结构由从下到上的氧化池、金属层和粒子床构成。 堆芯熔融物逐步塌落至下封头中,形成熔池。先期坍塌的熔融物在下封头中经历一段时间后形成具有氧化池和金属层的2层结构,随后熔化的熔融物经历一段时间后分解并溶入相应各层。在这段时间内,熔池是一个顶层覆盖了粒子床的3层结构。随着时间的推移,顶部粒子床(计算发现其主要为氧化物)分解并溶入相应各层,最终粒子床消失,熔池稳定于图1所示的2层结构。 4 计算事故序列的成本序列 4.1 堆型筛选和有效参数确定 综合国外核电站IVR措施有效性评价中的事故序列选取结果[1~3],并根据CPR1000堆型的一级概率安全分析结果,本文选取以下3种事故序列进行详细计算:小破口失水事故(SB-LOCA,具有最大的堆芯损坏频率)、大破口失水事故(LB-LOCA,具有最快的事故进程)和全场断电事故(SBO)。 本文针对以上3种事故序列,对有容器外冷却模型(假定向保温层注水的流量为100 kg/s,称为基准事故)和无此模型的情况进行了对比计算,并根据计算结果初步分析CPR1000堆型的IVR措施的有效性。 根据严重事故管理导则以及对AP系列核电厂的IVR措施有效性评价结果,只要保证反应堆冷却剂系统卸压和