堆芯预计发生DNB燃料棒数量计算的统计学方法论主要分为以下几个步骤：

1) 选取需要分析的特定工况，确定堆芯各个功率水平的燃料棒数量；

2) 使用堆芯热工水力子通道程序计算该工况下堆芯各处燃料棒表面的最小DNBR；

3) 使用统计学方法计算95%置信度水平下每根燃料棒表面最小DNBR处发生DNB的概率；

4) 计算堆芯预计发生DNB燃料棒数量的期望值。

此统计学方法论的核心步骤在于计算95%置信度水平下最小DNBR处发生DNB概率。

如果将真实DNBR记为z，那么z可以写为以下形式：

z = CHF M Φ (1)

其中：是试验测得的临界热流密度(CHF, Critical Heat Flux)；是实际热流密度。

z是一个随机变量，包含了以下几个方面的不确定度：

1) CHF关系式的不确定度；

2) 热工参数的不确定度；

3) 瞬态的不确定度；

4) 程序的不确定度。

其中热工参数的不确定度主要考虑以下热工参数：

1) 堆芯功率；

2) 一回路压力；

3) 反应堆压力容器(RPV, Reactor Pressure Vessel)入口冷却剂温度；

4) 核焓升热管因子；

5) 工程焓升热管因子。

对于给定的不确定度 U i ，其对应的标准差 σ i 由不确定度的分布类型决定：

-正态分布：

σ i = U i 1.645 (2)

-均匀分布：

σ i = U i 3 (3)

CHF关系式的不确定度与热工参数的不确定度一并计算。CHF关系式与参数的不确定度、瞬态的不确定度和程序的不确定度彼此独立。因此，总的不确定度所对应的标准差 σ z 可以写为以下形式：

σ z = σ 参 数 - 关 系 式 2 + σ 程 序 2 + σ 瞬 态 2 (4)

其中： σ 参 数 - 关 系 式 是CHF关系式和参数的标准差； σ 程 序 是程序的标准差； σ 瞬 态 是瞬态的标准差。

对于选定的工况， σ 参 数 - 关 系 式 通过以下统计方法计算：

1) 以选定工况的热工参数初始值为平均值，按照热工参数的分布类型和标准差，采用拉丁超立方的方法随机生成n (例如1000)个热工参数组(新的工况)；

2) 使用堆芯热工水力子通道程序计算选定工况和第1步随机生成的n个工况的最小DNBR；

3) 为了考虑CHF关系的不确定度，需要对第2步随机生成的n个工况的最小DNBR按以下算法进行修正：

DNBR 修 正 = DNBR × [ μ c + σ 关 系 式 × C ] (5)

其中：DNBR为随机生成工况的最小DNBR计算值； μ c 为随机变量 CHF M CHF P 的平均值， CHF P 是依据CHF关系式预测得到的CHF； σ 关 系 式 为CHF关系式的标准差；C为符合标准正态分布的随机变量； DNBR 修 正

为修正之后的DNBR值。

1) 基于n个修正之后的DNBR值，计算出随机变量 DNBR 修 正 样本的标准差 σ a l l ；

2) 基于随机变量 DNBR 修 正 样本的标准差 σ a l l ，保守计算出CHF关系式与参数的标准差：

σ 参 数 - 关 系 式 = ( ( n − 1 ) / χ n − 1 , 95 2 ) × σ a l l (6)

其中： χ n − 1 , 95 2 是95%概率和n-1自由度对应的卡方分布值 [ 4 ] 。

z符合平均值为 μ z ，标准差为 σ z 的正态分布：

z ~ N ( μ z , σ z 2 ) (7)

基于堆芯热工水力子通道程序和CHF关系计算得到的DNBR的定义如下：

DNBR = CHF P Φ (8)

在堆芯热工水力子通道程序和CHF关系式确定之后，DNBR的计算结果只取决于上述热工参数的取值。依据中心极限定理可知，无论上述热工参数符合何种分布，DNBR都符合正态分布。而DNBR分布的均值 DNBR ¯ 就等于当各个热工参数都取其最佳估算值时计算得到的DNBR。

由(1)式和(8)式可知，z可以写为以下形式：

z = CHF M CHF P × DNBR (9)

根据CHF试验能够得到 CHF M CHF P 的数据库。 CHF M CHF P 数据库的平均值 m c 和标准差 s c 直接计算即可。然而随机变量 CHF M CHF P 真实的平均值 μ c 和标准差 σ c 则是基于 m c 和 s c 计算得到。在95%置信水平下， μ c 和 σ c

可以按照以下公式保守计算得到：

μ c = m c − ( | t n − 1 , 95 | / n ) s c (10)

其中： t n − 1.95 是95%概率和n-1自由度对应的Student分布值 [ 5 ] 。

σ c = ( ( n − 1 ) / χ n − 1 , 95 2 ) s c (11)

当一根燃料棒表面的最小真实DNBR为z时，这根燃料棒表面发生DNB的概率 P ( DNB ) 有如下形式：

P ( DNB ) = P ( z < I ) (12)

其中：I是燃料棒弯曲的设计惩罚因子，I大于1。

由(7)式可知随机变量 z − μ z σ z 符合标准正态分布：

z − μ z σ z ~ N ( 0 , 1 ) (13)

由于 σ z 是非0正数，所以(12)式可以转化为如下形式：

P ( DNB ) = P ( z − μ z σ z < I − μ z σ z ) (14)

标准正态分布的概率密度函数如下：

f ( t ) = 1 2 π e − t 2 2 (15)

由(14)式和(15)式可知：

P ( DNB ) = 1 2 π ∫ − ∞ I − μ z σ z e − t 2 2 d t (16)

由(9)式可知：

μ z = μ c × DNBR ¯ (17)

将(17)式代入(16)式可得：

P ( DNB ) = 1 2 π ∫ − ∞ I − μ c × D N B R ¯ σ z e − t 2 2 d t (18)

根据以上推导，当一根燃料棒表面的最小DNBR值(最佳估算值)确定时，将其带入(18)式即可计算对应的此根燃料棒表面发生DNB的概率。

当堆芯内每一根燃料棒表面的最小DNBR值确定时，堆芯预计发生DNB燃料棒数量的期望值 E ( DNB ) 表达式如下：

E ( DNB ) = ∑ i P i ( DNB ) (19)

其中i为燃料棒的编号， P i ( DNB ) 为编号为i的燃料棒发生DNB的概率。

选取的典型工况为I、II类工况中非常恶劣的工况。该工况的运行参数如下表1所示。

表1. 典型工况运行参数

各种热工参数的分布类型和不确定度如下表2所示。

表2. 热工参数的分布及不确定度

堆芯轴向功率分布采用峰值为1.55的截尾余弦分布。

堆芯径向功率分布采用一种包络的参考堆芯径向功率分布，热组件在堆芯中心，热组件及其相邻组件的组件F_ΔH为1.52，子通道的峰值F_ΔH为1.59。

燃料棒统计曲线反映了堆芯内高于某一种棒功率水平的燃料棒份额。与典型工况分析采用的堆芯径向功率分布对应的燃料棒统计曲线如图1所示。图中横坐标表示F_ΔH数值，纵坐标表示堆芯中功率高于该F_ΔH数值的燃料棒份额。

图1. 燃料棒统计曲线

按照第2章介绍的统计学方法结合相关输入，计算了典型工况对应的预计发生DNB燃料棒的数量，如图2所示。

图2. 预计发生DNB燃料棒的数量

由图2可知，对于此典型工况，当F_ΔH为1.519时，预计发生DNB燃料棒数量的期望值最大，为0.205根。除此之外，其他F_ΔH值对应的预计发生DNB燃料棒数量的期望值均小于0.1根。图2的计算结果取决于各个F_ΔH值对应的燃料棒表面发生DNB的概率和燃料棒的数量。整个堆芯预计发生DNB的燃料棒数量的期望值仅为0.389根。