林诚格、易湘红
一、随着核电机组大量增加，核风险随之增大。这个非定量的叙述可以为大众理解。所以为了保持和降低总的核风险，就要提高单堆机组的安全水平。
二、应用概率论与数理统计知识，可以对群堆的安全水平进行定量的分析。
1、方法一[1]：
假设两个独立的事件a和b，他们各自发生的概率分别为p（a）和p（b）。我们可以计算出发生a事件或b事件的概率为：
p(a或b)=p(a) p(b)-p(a和b)
=p(a) p(b)-p(a)p(b) …公式(1)
如果是a、b、c三个独立事件,则
p(a或b或c)=p(a) p(b) p(c)-p(a和b)-p(a和c）
-p(b和c) p(a和b和c） …公式(2)
对于n个独立事件e1，e2，……en，则
p(e1或e2或……或en)
= - …
…-(-1)n p(e1和e2和……和en) …公式(3)
如果p(e1)=p(e2)= ……=p(en)=p
当概率p很小，且np<<1时，则
p（e1或e2或……或en）≈np …公式(4)
定义是不发生ei的状态,即p（）=1-p（ei）, 则有
p（至少有一个ei发生）=1-p（无ei发生）=1-p（和…）=1-p()p() …p()=1-{[1-p(e1)][1-p(e2)] …[1-p(en)]}
如果p(e1)=p(e2)= ……=p(en)=p, 则
p(至少有一个ei发生)=1-(1-p)n …公式(5)
2、方法二[2]：
定义 设将试验独立重复进行n次，每次试验中，事件a发生的概率均为p，则称这n次试验为n重贝努里试验。若以x表示n重贝努里试验事件a发生的次数，则称x服从参数为n,p的二项分布。记作x~b（n,p)分布律为：
至少发生一次的概率等于1减去一次也不发生的概率（也就是发生零次的概率），即p{x31}＝1－ p{x＝0}=1-(1-p)n，与公式(5)完全一致。二项分布的数学期望为np.
3、应用概率论与数理统计知识对核电机组安全性的定量分析。
(1) 计算堆芯熔化概率为10-4，100个堆一年发生事故的概率[3]
由于n=100,p=10-4,np<<1, 可直接应用公式(4), 得出p=np=0.01/年。
如应用公式（5）, 同样p{x31}＝1-(1-p)n =0.01/年
(2) 计算堆芯熔化概率为10-4，1000个堆一年发生事故的概率
当n=1000时, p{x31}＝1-(1-p)n =0.095，仍然可以用np来进行近似计算。
(3) 计算堆芯熔化概率为10-4， 1万堆年左右发生事故的概率。
当n=7000时, p{x31}＝1－ p{x＝0}=1-(1-p) n =0.5034
当n=10000时, p{x31}＝1－ p{x＝0}=1-(1-p) n =0.6321
当n=20000时, p{x31}＝1－ p{x＝0}=1-(1-p) n =0.8647
也就是接近一万堆年时，发生严重事故的可能性就比较大了。
不同的n和p值的事故概率表
发生事故的概率 | n=10 | n=100 | n=1000 | n=7000 | n=10000 | n=20000 |
二代 p=10-4 | 0.001 | 0.01 | 0.095 | 0.5034 | 0.6321 | 0.8647 |
ap1000 p=5×10-7 | 0.000005 | 0.00005 | 0.0005 | 0.0035 | 0.005 | 0.01 |
三、分析和结论
1、当np<<1时，有p（e1或e2或……或en）≈np，因为堆芯熔化概率为<10-4，几百个堆每年发生事故的概率满足np<<1的条件，所以群堆每年发生事故的概率可以近似于群堆数量乘以单堆熔化概率，与反应堆数量呈正比；
2、当np接近1时，发生事故的概率就比较大；
从上表中不难看出，当累计运行2万堆年时，对ap1000而言，发生事故的概率仍然可以视为小概率事件（0.01）；而二代机组发生事故的可能性已经比较大了(0.8647)。
这就是为什么要在推进核电规模发展的情况下要强调提高单堆机组安全水平的缘故。
参考资料:
[1] nureg-0492, fault tree handbook, usnrc, january,1981.
[2]盛骤等,概率论与数理统计，第三版，浙江大学, 2001年12月.
[3] 林诚格, 国际核电安全的发展历史对我国发展核电的影响，核电，2007年2月.