堆芯功率分布展平因子与焓升因子计算书

堆芯功率分布展平因子与焓升因子计算书第1页堆芯功率分布展平因子与焓升因子计算书文件编号：THERM-HYD-2026-CALC-001阶段：施工图设计阶段类型：热工水力计算报告编制：__________日期：2026年10月26日审核：__________日期：2026年__月__日批准：__________日期：2026年__月__日

目录TOC\o"1-2"\h\z\u1概述 12引用标准与依据文件 13术语、定义与符号说明 13.1功率展平因子（功率峰因子） 13.2焓升因子 23.3主要符号说明 24计算原理与数学模型 34.1功率展平因子计算模型 34.2焓升因子计算模型 35堆芯功率分布展平技术措施 55.1径向分区燃料装载 55.2可燃毒物优化布置 55.3反射层应用 55.4控制棒与化学补偿优化 55.5运行功率分布控制 56典型压水堆算例计算 66.1计算输入参数 66.2功率展平因子计算 66.3焓升因子计算 77安全核验与裕量分析 87.1设计限值要求 87.2计算结果核验 87.3不确定性与敏感性分析 88结果分析与工程指导意义 98.1计算结果综合分析 98.2对燃料管理的指导意义 98.3对运行控制的指导意义 108.4对安全分析的输入要求 109结论 10免责声明 12第1页1概述本计算书依据压水堆核电厂热工水力设计相关标准，对压水堆堆芯功率分布展平因子（功率峰因子）与焓升因子进行系统计算与安全核验。计算结果将为堆芯热工水力设计、燃料管理、运行限值制定以及安全分析提供基础依据，确保反应堆在各类工况下满足热工安全准则要求。堆芯功率分布不均匀性是核反应堆设计与运行中的核心问题之一。由于中子通量分布的固有不均匀性、燃料装载方式、控制棒位置、燃耗深度以及氙振荡等因素影响，堆芯内功率密度存在显著的空间差异。功率展平因子和焓升因子是定量表征这种不均匀性的关键参数，直接关系到燃料元件的热工裕量、偏离泡核沸腾比（DNBR）限值以及堆芯的整体安全性能。2引用标准与依据文件本计算书的编制严格遵循以下国家、行业标准及技术规范：（1）NB/T20057.2-2012《压水堆核电厂反应堆系统设计堆芯第2部分：热工水力设计准则》；（2）GB/T4960.1~4960.9《核科学技术术语》；（3）NB/T20035-2011《压水堆核电厂反应堆堆芯核设计准则》；（4）GJB843.18A-2017《潜艇核动力装置设计安全规定第18部分：压水型反应堆核设计要求》；（5）HAF102《核动力厂设计安全规定》；（6）《核反应堆热工水力学》（第3版），西安交通大学出版社；（7）《核反应堆热工分析》，重庆大学出版社。3术语、定义与符号说明3.1功率展平因子（功率峰因子）功率展平因子，又称功率峰因子或功率不均匀系数，是表征堆芯内功率分布不均匀程度的无量纲参数，定义为堆芯内最大功率密度与堆芯平均功率密度的比值。功率峰因子越小，说明堆芯功率分布越平坦，堆芯的热工裕量越大，经济性越好。功率展平因子包含三个维度：径向功率峰因子、轴向功率峰因子以及总功率峰因子。总功率峰因子为径向与轴向功率峰因子的乘积，考虑了三维空间内最极端的功率集中情况。3.2焓升因子焓升因子，又称热通道焓升因子，是表征堆芯内冷却剂焓升不均匀程度的关键参数，定义为热通道内冷却剂焓升与平均通道冷却剂焓升的比值。焓升因子直接决定了堆芯最热通道的出口焓值，是影响偏离泡核沸腾比（DNBR）计算的核心参数。焓升因子由核焓升因子和工程焓升因子两部分组成，分别考虑了核功率分布不均匀和工程因素（如流量分配偏差、制造公差等）造成的焓升不均匀性。3.3主要符号说明符号名称定义与说明F_r径向功率峰因子堆芯径向最大功率密度与平均功率密度之比F_z轴向功率峰因子燃料棒轴向最大功率密度与平均功率密度之比F_Q总功率峰因子三维空间最大功率密度与平均功率密度之比，F_Q=F_r×F_zF_ΔH^N核焓升因子仅考虑核功率分布不均匀的焓升因子F_ΔH^E工程焓升因子考虑工程不利因素的焓升因子F_ΔH总焓升因子热通道总焓升与平均通道焓升之比，F_ΔH=F_ΔH^N×F_ΔH^Eq'(z)局部线功率密度轴向高度z处燃料棒的线功率密度，W/cmDNBR偏离泡核沸腾比临界热流密度与实际热流密度之比

4计算原理与数学模型4.1功率展平因子计算模型堆芯功率分布本质上由中子通量密度分布决定，任意位置的功率密度与该位置的中子通量密度近似成正比。对于均匀裸圆柱形堆芯，热中子通量密度的径向分布为零阶贝塞尔函数分布，轴向分布为余弦分布。（1）径向功率峰因子F_r计算：径向功率峰因子定义为堆芯横截面上最大功率密度与该截面平均功率密度的比值，计算公式为：Fr=max[q'(r,θ)]/q'avg对于分区装载的实际压水堆堆芯，径向功率峰因子需要通过三维中子扩散计算或蒙特卡罗输运计算获得，同时需要考虑燃料组件内局部功率峰（棒功率峰）的影响。（2）轴向功率峰因子F_z计算：轴向功率峰因子定义为单根燃料棒上最大线功率密度与该棒平均线功率密度的比值，计算公式为：Fz=max[q'(z)]/(1/H)∫0Hq'(z)dz其中H为堆芯活性段高度。轴向功率分布受控制棒插入深度、氙振荡、慢化剂密度分布以及燃耗深度等因素显著影响。（3）总功率峰因子F_Q计算：总功率峰因子（三维功率峰因子）为径向功率峰因子与轴向功率峰因子的乘积，并需考虑方位角倾斜等局部功率畸变因素：FQ=Fr×Fz×FU其中F_U为计算不确定性因子，通常取1.03~1.05，用于包络核计算方法、核数据以及制造公差带来的不确定性。4.2焓升因子计算模型焓升因子是热通道因子的重要组成部分，表征了由于功率分布不均匀和工程因素导致的冷却剂焓升不均匀性。与功率峰因子不同，焓升因子是沿整个通道高度的积分效应，而非某一点的局部值。（1）核焓升因子F_ΔH^N计算：核焓升因子定义为最热燃料组件的积分功率与堆芯平均组件功率的比值，反映了径向功率分布不均匀性对整个通道焓升的贡献：FΔHN=(∫0Hq'hot(z)dz)/q'avg×H（2）工程焓升因子F_ΔH^E计算：工程焓升因子用于包络各种工程不利因素对焓升的影响，主要包括：流量分配不均匀、燃料组件制造公差、定位格架阻力偏差、入口流量偏差等。各子因子按乘积法则合成：FΔHE=Fflow×Fmech×Finlet（3）总焓升因子F_ΔH计算：FΔH=FΔHN×FΔHE

5堆芯功率分布展平技术措施为降低功率峰因子、展平堆芯功率分布，现代压水堆核电厂普遍采用以下多项技术措施进行综合优化：5.1径向分区燃料装载将堆芯按径向分为多个富集度区域，在中心区域布置较低富集度的燃料组件，在外围区域布置较高富集度的燃料组件。通过这种“低内高外”的装载方式，补偿因中子泄漏导致的外围功率下降，有效展平径向功率分布。典型百万千瓦级压水堆通常采用三区装载方案，初始富集度分别为1.8%、2.4%、3.1%左右。5.2可燃毒物优化布置在功率较高的中心区域或功率峰位置布置含可燃毒物（如硼、钆、铒等）的燃料棒，通过可燃毒物的中子吸收作用压低局部功率峰值。随着燃耗加深，可燃毒物逐渐消耗，其抑制功率的作用也逐渐减弱，与燃料燃耗反应性损失相匹配，实现全寿期内的功率分布展平。一体化可燃毒物（如Gd₂O₃-UO₂芯块）是目前的主流技术。5.3反射层应用在堆芯外围设置水反射层或径向反射层，利用反射层将泄漏出堆芯的中子散射回堆芯，提高堆芯边缘区域的中子通量水平，从而降低径向功率峰因子。合理的反射层设计可使径向功率峰因子降低约0.1~0.2。5.4控制棒与化学补偿优化采用化学补偿控制（硼溶液）减少控制棒的插入数量，避免控制棒造成的严重功率畸变。采用部分长度控制棒（灰棒）进行轴向功率形状控制，抑制轴向功率峰。优化提棒程序，采用“提棒外推”策略，避免在堆芯中部形成大的功率峰。5.5运行功率分布控制运行期间通过堆芯仪表系统实时监测轴向功率偏移（AO）和象限功率倾斜比（QPTR），将功率偏移控制在目标带（通常±5%）范围内，防止氙振荡导致的功率分布畸变。实施运行图控制，确保在负荷跟踪、功率升降等瞬态过程中功率峰因子不超限值。6典型压水堆算例计算6.1计算输入参数以某典型百万千瓦级三代压水堆为计算对象，堆芯主要设计参数如下表所示：参数名称数值单位反应堆额定热功率3050MW堆芯活性段高度365.8cm燃料组件数177组燃料组件排列17×17-平均线功率密度173.8W/cm堆芯入口温度292.4℃堆芯出口平均温度327.2℃一回路压力15.5MPa堆芯冷却剂总流量68500m³/h6.2功率展平因子计算基于三维中子扩散计算程序，对堆芯寿期初（BOL）、寿期中（MOL）、寿期末（EOL）三个典型燃耗状态进行功率分布计算，得到各功率峰因子如下表：燃耗状态径向功率峰因子F_r轴向功率峰因子F_z总功率峰因子F_Q寿期初（BOL）1.3281.4521.987寿期中（MOL）1.2851.4161.872寿期末（EOL）1.2641.3981.816全寿期最大值1.3281.4521.987计算结果表明：总功率峰因子最大值为1.987，出现在寿期初零氙状态。随着燃耗加深，可燃毒物逐渐消耗，功率分布逐渐展平，总功率峰因子呈下降趋势。采用三区装载加钆可燃毒物的设计方案有效展平了堆芯功率分布。6.3焓升因子计算根据径向功率分布计算各燃料组件的积分功率，结合工程不利因子分析，得到焓升因子计算结果如下表：燃耗状态核焓升因子F_ΔH^N工程焓升因子F_ΔH^E总焓升因子F_ΔH寿期初（BOL）1.4321.0721.535寿期中（MOL）1.3951.0721.495寿期末（EOL）1.3681.0721.466全寿期最大值1.4321.0721.535其中工程焓升因子由以下子因子组成：流量分配因子1.04，机械制造公差因子1.02，入口流量偏差因子1.01，按乘积法则合成得到1.04×1.02×1.01≈1.072。工程焓升因子在整个寿期内取保守恒定值。

7安全核验与裕量分析7.1设计限值要求根据NB/T20057.2-2012及三代压水堆设计要求，功率展平因子与焓升因子需满足以下安全限值：参数设计限值限值依据总功率峰因子F_Q≤2.35燃料线功率密度限值、LOCA分析限值总焓升因子F_ΔH≤1.65DNBR限值、15%热工裕量要求最大线功率密度q'_max≤410W/cm燃料中心熔化温度限值最小DNBR≥1.18W-3临界热流密度关系式7.2计算结果核验将本计算得到的全寿期最大功率展平因子和焓升因子与设计限值进行对比核验，结果如下：（1）总功率峰因子：计算最大值为1.987，小于设计限值2.35，裕量为(2.35-1.987)/2.35×100%=15.4%，满足安全要求。（2）总焓升因子：计算最大值为1.535，小于设计限值1.65，裕量为(1.65-1.535)/1.65×100%=7.0%，满足URD要求的至少15%热工裕量中的焓升因子部分要求。（3）最大线功率密度：q'_max=F_Q×q'_avg=1.987×173.8=345.3W/cm，小于限值410W/cm，裕量为15.8%，确保燃料中心温度低于熔点。（4）DNBR核验：基于子通道分析程序计算，在额定工况下考虑F_ΔH=1.535时，全堆芯最小DNBR为1.37，大于限值1.18，具有充足的偏离泡核沸腾裕量。7.3不确定性与敏感性分析为确保计算结果的保守性，对主要不确定性来源进行敏感性分析：（1）核计算不确定性：采用统计抽样方法对核数据、计算模型不确定性进行分析，95%概率、95%置信度下F_Q和F_ΔH的单侧上限不超过计算值的1.05倍，已在计算中通过不确定性因子予以考虑。（2）燃耗不确定性：考虑燃耗计算偏差±5%时，功率峰因子变化不超过2%，焓升因子变化不超过1.5%，对安全裕量影响较小。（3）控制棒插入敏感性：在最大价值控制棒卡棒极限工况下，F_Q最大值增加至2.12，F_ΔH最大值增加至1.58，仍满足安全限值要求。（4）氙振荡敏感性：在氙振荡导致的最不利功率分布下，轴向功率偏移可达到±10%，此时F_Q最大值约为2.18，仍有7.2%的安全裕量。8结果分析与工程指导意义8.1计算结果综合分析本计算结果表明，采用三区装载加一体化钆可燃毒物的堆芯设计方案，配合优化的控制棒策略和功率分布控制，能够在整个燃料循环寿期内将功率展平因子和焓升因子控制在安全限值以内，并保持充足的热工裕量。从寿期变化趋势来看，功率峰因子和焓升因子均随燃耗加深而逐渐降低，说明可燃毒物的燃耗特性与功率展平需求匹配良好。寿期初为功率分布控制的最不利阶段，这一结论与工程运行经验一致。径向功率峰因子（1.26~1.33）显著低于轴向功率峰因子（1.40~1.45），说明轴向功率分布是制约总功率峰因子的主要因素。因此，运行期间应重点关注轴向功率形状控制，通过合理的控制棒移动和硼浓度调节，将轴向功率偏移维持在目标范围内。8.2对燃料管理的指导意义（1）换料设计优化：本计算结果验证了三区装载方案的功率展平效果，后续换料设计可在此基础上进一步优化高燃耗组件的布置位置，在满足功率峰限值的前提下加深卸料燃耗，提高燃料经济性。（2）可燃毒物配置：钆棒数量和布置位置对径向功率分布影响显著。寿期初功率峰主要受钆棒配置控制，应避免钆棒过度集中导致局部功率压低过多、反而引起外围功率升高的反效果。（3）长循环换料适应性：若将燃料循环长度从18个月延长至24个月，需要重新评估可燃毒物装载量和功率峰因子变化，确保长循环下功率分布仍满足限值。8.3对运行控制的指导意义（1）功率分布运行图：基于本计算得到的功率峰因子限值，制定堆芯运行图（PowerCoastdownChart），明确不同功率水平和燃耗深度下的允许轴向功率偏移范围，为运行人员提供操作指导。（2）负荷跟踪策略：在负荷跟踪运行期间，应采用灰棒控制而非硼稀释来跟踪负荷变化，减少对轴向功率分布的扰动。建议采用机械补偿（MSHIM）运行模式时，严格控制灰棒插入深度，避免局部功率峰过大。（3）氙振荡抑制：当轴向功率偏移超出目标带时，应及时通过控制棒和硼浓度调整进行纠正，防止氙振荡发散导致功率峰因子持续升高。（4）功率倾斜监督：密切监测象限功率倾斜比（QPTR），当QPTR超过2%时启动降功率程序并查找原因，防止局部功率异常升高危及燃料安全。8.4对安全分析的输入要求本计算得到的功率展平因子和焓升因子是后续安全分析的重要输入参数：（1）LOCA分析：大破口失水事故分析中需采用F_Q最大值作为初始条件，确保峰值包壳温度（PCT）低于1204℃的限值。（2）DNB分析：在二类、三类工况瞬态分析中，需采用F_ΔH最大值计算DNBR变化，确保