压水堆一回路pH控制策略对积垢燃料包壳完整性的影响研究

蒙舒祺，毛玉龙，胡艺嵩，阮天鸣，胡友森

压水堆一回路pH控制策略对积垢燃料包壳完整性的影响研究

蒙舒祺，毛玉龙，胡艺嵩，阮天鸣，胡友森*

（中广核研究院有限公司，广东 深圳 518000）

压水堆（Pressurized Water Reactor，PWR）通常采用不同的硼酸和氢氧化锂浓度配比维持一回路pH。一回路各种材料在高温高压水环境下长期服役，其腐蚀产物在燃料包壳表面沉积形成污垢（Chalk Rivers Unidentified Deposits，CRUD），污垢内部的灯芯沸腾会造成硼酸和氢氧化锂局部浓缩，改变污垢内部冷却剂的热物性，进而影响污垢和燃料包壳交界面处的温度和pH。由于温度和pH对燃料包壳局部腐蚀及裂纹萌生拓展影响显著，因此需要耦合热力学和化学动力学评估一回路pH控制策略对积垢燃料包壳完整性的影响。研究结果表明，采用天然硼和氢氧化锂控制策略的积垢燃料包壳完整性失效风险最高，采用无硼运行控制策略的积垢燃料包壳完整性失效风险最低。

压水堆；pH控制策略；污垢；燃料包壳；完整性

目前国内外较成熟的PWR堆型大多使用天然硼和氢氧化锂作为控制一回路pH的药剂[1]，而新研发的堆型一回路pH控制策略有所不同，如台山CEPR堆型使用富集硼和氢氧化锂作为药剂[2]，韩国设计的一体化小型PWR采用无硼策略控制一回路pH[3]。长期在一回路高温高压水环境下服役的金属材料释放的腐蚀产物会在燃料表面形成CRUD[4,5]，CRUD多为疏松多孔结构，从主流体进入CRUD内部的冷却剂在灯芯沸腾效应作用下会造成积垢燃料表面硼锂浓度富集，影响局部冷却剂热物性，进而导致燃料包壳表面温度和pH发生显著变化[6-8]。国内外已有研究表明，表面温度升高、pH降低会加速燃料包壳局部腐蚀和裂纹萌生拓展[9-11]。不同一回路pH控制策略对积垢燃料表面温度和pH的影响存在差异，因此需要耦合热力学和化学动力学评估一回路pH控制策略对积垢燃料包壳完整性的影响。

1　CRUD内部的热力学和化学动力学

1.1　溶液热物性参数计算模型

国际上对PWR的检测结果表明CRUD大多疏松多孔且零星分布着蒸汽通道[12,13]，因此通常将CRUD形貌抽象成由均匀分布的蒸汽通道和固体区域组成的结构[6,14]，蒸汽通道和固体区域分别以蒸发换热和导热形式传递能量。含硼酸和氢氧化锂的冷却剂从固体区域流入CRUD内部，在靠近燃料包壳表面处发生灯芯沸腾，大部分纯水以气态形式从蒸汽通道回到主流体中重新冷凝，而硼酸和氢氧化锂在局部富集[6,15]。基于能量守恒可得到沿CRUD厚度方向的温度分布关系式[6]，解析解形式为：

其中：——冷却剂温度，K；

s——冷却剂饱和温度，K；

——CRUD固体区域占比；

c——导热系数，W/（m·K）；

c——蒸汽通道平均半径，mm；

c——蒸汽通道分布密度，个/mm2；

e——沸腾换热系数，W/（m2·K）。

在PWR主流体中硼酸和氢氧化锂浓度相对较低，可近似用纯水表征冷却剂热物性。当硼酸和氢氧化锂在CRUD内部富集后，根据非理想稀溶液定义[16]，用纯水表征冷却剂热物性所带来的偏差较大，因此需要考虑溶液活度变化对冷却剂饱和温度的影响。Deshon给出了冷却剂饱和温度与溶液活度之间的经验关系式[17]：

其中：——冷却剂活度；

CRUD内部溶液活度的变化还会影响汽化潜热和沸腾换热系数，对应的经验关系式如下：

1.2　硼锂浓度富集计算模型

假设硼锂在CRUD内部的流动和富集瞬间达到稳态，同时忽略灯芯沸腾造成的硼锂蒸汽挥发，可得到表征沿CRUD厚度方向硼锂浓度变化的微分方程：

其中：——厚度为处的硼锂浓度，10-6；

bulk——主流体中硼锂浓度，10-6；

0——富集的硼锂浓度，10-6；

——硼锂在CRUD内部的流速，cm/s；

——硼锂扩散系数[18]，cm2/s。

基于能量守恒，燃料产生的热量等于灯芯沸腾带走的热量，由此可得：

其中：v——蒸汽换热面积，cm2；

——燃料包壳总换热面积，cm2。

对给定厚度的CRUD，根据孔隙率的定义[19]：

联立式（8）～式（12）得到硼锂富集浓度与主流体浓度之间的关系式：

1.3　CRUD内部pH计算模型

纯水、硼酸和氢氧化锂在溶液中的电离，会影响沿CRUD厚度方向的pH分布，相关的化学方程式如下：

式（14）～式（18）对应的化学电离平衡常数表达式为[20-22]：

根据酸碱电离平衡和电荷守恒，可得到用于迭代计算pH的表达式：

2　pH控制策略对燃料包壳完整性的影响

2.1　评估方法及结果

在编程软件中对CRUD内部的热力学和化学动力学耦合模型进行了功能实现。首先根据输入的热工水力参数和硼锂浓度，完成温度和硼锂浓度迭代计算；然后根据温度和硼锂浓度，完成pH迭代计算。沿CRUD厚度方向的溶液温度和pH迭代计算过程如图1所示。

硼酸中的同位素10B能够控制PWR反应性，由于富集硼中10B的丰度更高，在控制反应性需求不变的前提下，相对于天然硼，使用富集硼后主流体中硼酸浓度更低。在给定的热工水力条件和CRUD总体参数[6,23,24]基础上，分别构造3组算例模拟天然硼+氢氧化锂（1 500×10-6Boron+2×10-6LiOH）、富集硼+氢氧化锂（800×10-6Boron+2×10-6LiOH）和无硼控制（0 Boron+0.5×10-6LiOH）策略下燃料包壳表面温度和pH的变化情况。评估使用的输入参数如表1所示，这些参数取值与PWR一回路真实运行工况相近，具有一定的代表性。

图1　CRUD内部温度和pH迭代计算的逻辑图

表1　用于评估pH控制策略的输入参数

不同pH控制策略下沿CRUD厚度方向的温度和pH分布计算结果，如图2、图3所示。可以看出：

（1）采用目前应用最为广泛的pH控制策略（天然硼+氢氧化锂），积垢燃料包壳表面的温度最高、pH值最低，燃料包壳完整性失效风险最高；

（2）采用富集硼+氢氧化锂的pH控制策略，由于主流体硼浓度较低，在CRUD内部中硼酸的富集也相对较少，对溶液热物性的影响程度下降，因此积垢燃料包壳完整性失效的风险也降低；

（3）采用无硼运行的pH控制策略，氢氧化锂对溶液热物性的影响程度较小，此时积垢燃料包壳表面温度和pH与主流体相近，燃料包壳完整性失效风险最低。

图2　不同pH控制策略下CRUD厚度方向温度分布

图3　不同pH控制策略下CRUD厚度方向pH分布

2.2　对评估结果的讨论

当前的评估结果可初步反映一回路pH控制策略对积垢燃料包壳完整性的影响，但针对实际运行的PWR开展评估工作，还需要在以下方面开展更深入的研究：

（1）目前已有实验结果表明，富集硼可抑制PWR一回路金属材料的腐蚀行为[25]，从而降低CRUD厚度。针对采用富集硼和氢氧化锂作为一回路pH控制策略的PWR，需要在当前耦合的热力学和化学动力学模型基础上，考虑一回路pH控制策略对材料腐蚀和CRUD增长速度的影响；

（2） PWR燃料表面沉积的CRUD主要成分为NiO和NiFe2O4组成的混合尖晶石[17,26]，CRUD内部温度和pH值的变化会影响混合尖晶石氧化物的溶解度，改变CRUD沉淀溶解平衡状态及固相[27,28]，进而影响CRUD厚度及导热系数。因此，针对实际运行的PWR，还需要在当前耦合的热力学和化学动力学模型基础上，考虑CRUD增长速度和固相变化对其内部温度和pH影响的反馈效应；

（3）对实际运行的PWR，硼在CRUD中的富集不仅会造成局部pH值和温度分布异常，还可能导致燃料组件轴向功率偏移和偏离泡核沸腾起始点改变，进一步加剧燃料包壳完整性失效的风险。后续将在当前工作的基础上，耦合堆芯物理与热工水力模型，进一步提高对PWR真实工况的模拟精度。

（4）大量PWR运行经验表明，燃料组件轴向的温度、功率密度、CRUD厚度及固相分布差异显著[5,29]，针对实际运行的PWR，需要系统模拟燃料组件轴向CRUD内部温度和pH分布，才能更准确地评估积垢燃料包壳发生完整性失效风险最大的位置。

3　结论

（1）PWR一回路腐蚀产物形成的CRUD在热力学和化学动力学耦合作用下，会造成积垢燃料包壳表面温度和pH发生显著改变，进而影响燃料包壳完整性；

（2）对积垢燃料包壳，在不考虑硼酸对燃料组件轴向功率分布和偏离泡核沸腾起始点影响的前提下，采用天然硼和氢氧化锂控制策略的完整性失效风险最高，采用富集硼和氢氧化锂风险次之，采用无硼运行的风险最低；

（3）对实际运行的PWR，还需要考虑pH控制策略对材料腐蚀、CRUD生长速度、燃料组件轴向CRUD分布的影响，才能准确评估积垢燃料包壳发生完整性失效风险最大的位置。

［1］ Millett．P．PWR Primary Water Chemistry Guidelines，Volume 1，Revision 4［R］．EPRI Technical Report，TR-105714-V1R4，1999．

［2］ 王永刚，向文欣，姚波，等．台山核电厂CEPR堆芯设计及燃料管理［J］．核动力工程，2015，36（S1）：51-53．

［3］ 张国旭，解衡，谢菲．小型模块式压水堆设计综述［J］．原子能科学技术，2015，49（S）：40-47．

［4］ Joe．H．J，Jones．B．G，RAGHEB．M．M．Axial Offset Anomaly and Radiolysis Product Transport in PWRs Fuel Assembly Crud Layer［C］．International Conference on Nuclear Technology：Achieving Global Economic Growth While Safeguarding the Environment，November 16-20，New Orleans，USA，2003．

［5］ Deshon．J．PWR Axial Offset Anomaly（AOA）Guidelines［R］．EPRI Technical Report，1008102，2004．

［6］ Henshaw．J，McGURK．J．C，Sims．H．E，et al．A Model of Chemistry and Thermal Hydraulics in PWR Fuel Crud Deposits［J］．Journal of Nuclear Materials，2008，353（1-2），1-11．

［7］ Walter．D．J．A High Fidelity Multiphysics Framework for Modeling CRUD Deposition On PWR Fuel Rods［D］．University of Michigan，USA，2016．

［8］ Short．M．P．The particulate Nature of the Crud Source Term in Light Water Reactors［J］．Journal of Nuclear Materials，2018，509：478-481．

［9］ Lange．T．L．Methodology for an Advanced Risk Assessment of Crud Induced Power Shift Using Coupled Multiphysics Simulations and a Monte Carlo Scenario Analysis of the Potential Financial Benefits［D］．University of Tennessee，USA，2017．

［10］ Wang．G，Young．M．Y，Byers．W．A，et al．Methods to Reduce CIPS/CILC Risk for The Zero Fuel Failure by 2010 Initiative［C］．Proceedings of the 18th International Conference on Nuclear Engineering，Xi’an，China，2010．

［11］刘新福．压水堆一回路水化学对燃料包壳完整性的影响［D］．上海交通大学，2007．

［12］ Deshon．J．Evaluation of Fuel Cladding Corrosion and Corrosion Product Deposits from Callaway Cycle 13，Results of Poolside Examinations Following One Cycle of Zinc Addition［R］．EPRI Technical Report，1011088，2005．

［13］ IAEA．Water Chemistry and Corrosion Control of Cladding and Primary Circuit Components［R］．IAEA Technical Document，IAEA-TECDOC-1128，1999．

［14］ Haq．I．U，Cinosi．N，Bluck．M，et al．Modelling Heat Transfer and Dissolved Species Concentrations within PWR Crud［J］．Nuclear Engineering and Design，2011，241：155-162．

［15］ Frattini．P．L，Blok．J，Chauffriat．S，et al．Axial Offset Anomaly：Coupling PWR Primary Chemistry with Core Design［J］．Nuclear Energy，2001，40（2）：123-135．

［16］傅献彩，沈文霞，姚天扬，等．物理化学［M］．5版．北京：高等教育出版社，2005：233-252．

［17］ Deshon．J．Modeling PWR Fuel Corrosion Product Deposition and Growth Processes［R］．EPRI Technical Report，1009734，2004．

［18］ Manacorda．A，Schehr．G，Zamponi．F．Numerical Solution of the Dynamical Mean Field Theory of Infinite- Dimensional Equilibrium Liquids［J］．The Journal of Chemical Physics，2020，152（16）：164506．

［19］ Zou．L，Zhang．H．B，Gehin．J，et al．Coupled Thermal-Hydraulic/Neutronics-Crud Framework in Prediction of Crud Induced Power Shift Phenomenon［J］．Nuclear Technology，2013，183（3）：535-542．

［20］ Marshall．W．L，Frank．E．Ion Product of Water Substance，0～100 ℃，1-10，000 Bars，New International Formulation and its Background［J］．The Journal of Chemical and Physics，1981，10（2），295-304．

［21］ Harvey．A．H，Friend．D．G．Aqueous Systems at Elevated Temperatures and Pressures：Physical Chemistry in Water，Steam and Hydrothermal Solutions，Chapter 1，Physical Properties of Water［M］．USA：Elsevier Ltd，2003：478．

［22］ Park．B．G，Seo．S，Kim．S．J，et al．Meso-Scale Multi-Physics Full Coupling within Porous CRUD Deposits on Nuclear Fuel［J］．Journal of Nuclear Materials，2018，512：100-117．

［23］ Yeo．D．Y，NO．H．C．Modeling Heat Transfer Through Chimney-Structrued Porous Deposit Formed in Pressurized Water Reactors［J］．International Journal of Heat and Mass Transfer，2017，108：868-879．

［24］ Wang．G，Byers．W．A，Young．M．Y，et al．Thermal Conductivity Measurements for Simulated PWR Crud［C］．Proceedings of the 2013 21st International Conference on Nuclear Engineering，Chengdu，China，2013．

［25］温菊花，姜峨，刘金华，等．大型压水堆一回路冷却剂富集10B的硼锂水化学应用研究［J］．腐蚀与防护，2015，36（S2）：208-212．

［26］ Betova．I，Bojinov．M，Saario．T．Start-Up and Shut-Down Water Chemistries in Pressurized Water Reactors［R］．VTT Research Report，VTT-R-00699-12，2012．

［27］李璐．水冷反应堆主回路腐蚀产物活化及迁移模型的研究［D］．华北电力大学，2017．

［28］高希龙，景福庭，刘嘉嘉．富集硼酸对压水堆一回路腐蚀产物沉积的影响研究［J］．中国核电，2020，13（4），428-432．

［29］ Riess．R．Chemistry Experience in The Primary Heat Transfer Circuit of Kraftwerk Union Pressurized Water Reactors［M］．Nuclear Technology，Taylor & Francis Group，2017：153-159．

The Effect of PWR Primary pH Controlling Strategy on Fuel Cladding Integrity with CRUD

MENG Shuqi，MAO Yulong，HU Yisong，RUAN Tianming，HU Yousen*

（China Nuclear Power Technology Research Institute，Shenzhen of Guangdong Prov. 518000，China）

The PWR primary pH is normally maintained by compounding boron acid and lithium hydroxide concentrations. Corrosion products from multiple materials served in high temperature pressurized water may deposit on fuel surface and form CRUD. Wick boiling inside CRUD causes local enrichment of boron acid and lithium hydroxide, and the thermal properties of coolant inside CRUD correspondingly change. Consequently, the pH values and temperature of CRUD-fuel cladding interface will be affected. Since the temperature and pH values are key factors on localized corrosion and crack initiation and growth of fuel cladding, it is necessary to assess the effect of pH controlling strategy on PWR fuel cladding integrity with CRUD by coupling thermal and chemical dynamics. Research results show that the highest fuel cladding integrity disable risks is corresponding to natural boric acid and lithium hydroxide strategy, and the lowest risks is corresponding to non-lithium hydroxide strategy.

PWR; pH controlling strategy; CRUD; Fuel cladding; Integrity

TL349

A

0258-0918（2022）06-1242-06

2022-03-12

国家自然科学基金（U20B0211，针对堆芯氧化腐蚀产物材料-热工-中子行为的多物理耦合机理）；国家自然科学基金（52171085，模拟压水堆一回路冷却剂中燃料包壳管表面污垢沉积行为与机理研究）

蒙舒祺（1992—），男，广东惠州人，工程师，硕士，现主要从事反应堆一回路热工水力和水化学研究

胡友森，E-mail：huyousen@cgnpc.com.cn