核素裂变截面分析.docx
- 文档编号:5681609
- 上传时间:2022-12-31
- 格式:DOCX
- 页数:14
- 大小:286.66KB
核素裂变截面分析.docx
《核素裂变截面分析.docx》由会员分享,可在线阅读,更多相关《核素裂变截面分析.docx(14页珍藏版)》请在冰豆网上搜索。
核素裂变截面分析
核素裂变截面分析
6.1核反应截面测量原理
6.1.1裂变截面的计算公式推导
假定样品质量为M,母核X的原子量为A,丰度为η,测量时间为T3,反应道X(n,b)Y的截面值为σ,样品辐照时间为T1,辐照期间的平均中子注量率为φ(由27Al(n,α)24Na监测反应测得)冷却时间为T2,如图6.1.1所示。
t0t1t2t3
图6.1.1时间示意图
根据人工放射性随时间的生长规律[20],在样品辐照时间T1内的任一时刻t,单位时间内产生放射性子核Y的净数目为:
(6.1.1)
式中,NA为阿伏伽德罗常量,λ为生成子核Y的衰变常量。
由初始条件t=0时,N=0解式(6.1.1)可求得,在辐照时间T1内的任一时刻t,样品中生成子核Y的数目为:
(6.1.2)
照射结束(t=T1)时,样品中生成子核Y的数目为:
(6.1.3)
式中,
称为饱和因子。
开始测量(t=T1+T2)时,样品中生成子核Y的数目为:
(6.1.4)
设生成的放射性核Y放出的特征γ射线的分支比为Iγ,探测器对其全能峰探测效率为
,在测量时间T3内,测到特征γ射线的全能峰计数为C,计数的总校正因子F=Fs*Fc*Fg(随后给出计算公式的详细推导),中子注量率随时间波动的校正因子K,则有效的特征γ射线的全能峰计数应为:
(6.1.5)
式中,
称为测量收集因子。
得,
(6.1.6)
上式为绝对法测量核反应截面的一般计算公式,在实际实验中的中子注量率φ的数值很难准确测定,故而使用相对法测量更容易得出截面值。
顾名思义,相对法就是在实验计算中使用某种核素的核反应已知准确截面值。
它的子核与待测核反应的子核半衰期相差不大,用其做成两块监督片来测量中子注量率,在实验中将待测样品夹在这两块监督片中间,辐照完成后,把已知监督片的准确核反应截面值带入绝对法测量核反应截面计算公式中,再使用实验测得的相关数据就可得到准确的平均中子注量率φ的数值。
利用得到的φ就可以求解待测样品的核反应截面数值。
使用下标0来表示监督片相应的量,下标x来表示待测核反应的相应的量,则根据上式,可得
(6.1.7)
(6.1.8)
通过对上面两个公式的联立(两个公式中中子注量率是相等的),得到
(6.1.9)
上式即为相对法测量核反应截面中的常用公式。
6.1.2β-衰变引起的计数扣除
以计算过程中遇到的138Xeβ-衰变到138Cs为例:
先忽略中子注量率波动时,在中子辐照期间T1内138Xe和138Cs(仅由138Xe的β-衰变产生)的计数N1(t)和N2(t)为:
(6.1.10)
(6.1.11)
(6.1.12)
(6.1.13)
其中Xit代表138Xe发生β-衰变的概率,下标1和2分别对应138Xe和138Cs的物理量。
同理,在冷却期间T2内138Xe和138Cs(仅由138Xe衰变产生的)的计数N1(t')和N2(t')为:
(6.1.14)
(6.1.15)
在测量时间T3内,138Cs(仅由138Xe衰变产生的)的计数N2(t")为:
(6.1.16)
从而,138Cs(仅由138Xe衰变产生的)的全能峰计数为:
(6.1.17)
当考虑中子通量注量率修正因子K时,上式变为:
(6.1.18)
(6.1.19)
其中,
(6.1.20)
因此,232Th(n,f)138Cs裂变反应的全能峰计数C'为:
(6.1.21)
6.1.3校正因子F的计算
(1)特征γ射线在样品中的自吸收校正因子Fs
γ射线与物质相互作用,穿过该物质,射线强度由I0变为I,它们有关系
(式中,μm为样品对某种特征γ射线的质量吸收系数,可以在NationalInstituteofStandardsandTechnology网站查出;dm为质量厚度,可计算得到),对于低能γ射线与原子序数高的样品,发生光电效应占优势,此时γ射线在样品中的自吸收变得尤为显著,所以对自吸收的校正非常重要。
样品的自吸收校正因子Fs为:
(6.1.22)
(2)样品几何校正因子Fg
样品几何校正因子Fg的计算公式为:
(6.1.23)
式中,d为样品的厚度,h为样品距HPGe晶体表面的距离。
(3)特征γ射线级联符合校正因子Fc[21-23]
核反应中伴随的各条γ射线的寿命一般较短,它们大多数是有级联的,由于探测器的分辨时间的限制,可能不能分辨级联辐射中的各条γ,发生级联符合效应,这就会使测得的特征γ射线的全能峰计数增加(或减少),级联符合的程度与测量距离和级联链的个数有关。
造成测量的准确性,因此需要对测得的全能峰计数进行修正。
由于级联符合校正较为复杂,计算难度高。
我们仅对简单的三个能级的衰变级联符合修正进行讨论,衰变纲图如下,
图6.1.2简单的三个能级衰变纲图
对于γ20、γ10和γ21的级联符合校正因子F20、F10和F21分别为:
(6.1.24)
(6.1.25)
(6.1.26)
上述几个公式中,因子
;
;
;
Pij为γij的发射率,αij为γij的内转换系数,xij为与γij对应的分支比;
表示母核衰变到子核i能级的分支比;
、
分别表示探测器对γij的全能峰探测效率和全谱探测效率。
(4)中子注量率校正因子K
在实验进行中,中子源出射中子并不是绝对稳定的,中子注量率是会随时间波动的。
中子注量率波动校正因子K的计算公式为:
(6.1.26)式中,
为照射时间内的平均中子注量率。
6.2232Th(n,f)X截面计算
6.2.1监督片Al相关参数计算
对截面计算6.1.9式
观察,可令式中与监督片Al有关量相关为
,则截面计算公式变为
,则对于每一个样品来说,B0是一个常数,所以首先确定B0的值是截面计算的重点。
测得监督片的27Al(n,α)24Na反应道的γ能谱只有一个峰,全能峰的值为1368.63MeV,由NNDC数据库(http:
//www.nndc.bnl.gov/)给出的截面可以查出,活化中子能量为14.1MeV,14.5MeV和14.7MeV时,σ0的值分别为121.2055mb,115.9875mb和113.3785mb。
表6.2.1监督片测量相关读数值
样品
照射时间T1
冷却时间T2
测量时间T3
Al(计数C)
Al(I_gamma)
1
6420
6108
1492.06
5887
0.999936
2
7520
5955
1219.12
5369
0.999936
3
7080
6240
1230.30
6413
0.999936
表6.2.2监督片相关因子计算值
样品
Al(S)
Al(D)
Al(ε)
K(Al)
B0
1
0.0791
0.0175
0.0024256
0.9993
1.44257E-05
2
0.0920
0.0144
0.0024256
1.0320
1.68499E-05
3
0.0869
0.0145
0.0024256
1.0003
1.24162E-05
6.2.2样品截面分析
简化后的样品截面计算公式为
,其中B0因子已经计算得出。
其余样品的相关因子的具体计算过程数据在正文不再一一列出,计算方法参照6.1中已详细给出的相关公式。
在对应每个核素截面分析时,根据其半衰期选择相应测量时间合适的能谱进行峰计数的读取,一般选择4~5个半衰期为测量时间。
从前面对探测器效率刻度知道,所使用的标准源的能量最低为121.7824keV,而理论的刻度曲线应该是先有一段上升的曲线,但刻度结果只有下降部分的,所以在截面分析中能量低于121.7824keV的峰应该舍去。
利用相关核数据库NEDF(如:
半衰期、伽玛能量和分支比等),得到各232Th(n,f)X裂变通道截面数据,如下表6.2.3~6.2.6所示,。
(1)232Th(n,f)138Xe截面分析
由核素衰变纲图可以知道,138Xe的生成不只一种途径,除了样品的裂变产生,还可以由:
138Iβ-138Xe产生,但由于138I的半衰期只有6.49s,存在时间极短,实验使用的活化法只能甄别出中等半衰期的放射性核素,所以138I的存在与否,不能确定,如果它存在,则对本反应道的截面计算有一些影响,但由于实验无法甄别出如此短半衰期的核素,在认为138Xe的截面不存在β-衰变的先驱核的影响下,计算结果如下。
表6.2.3232Th(n,f)138Xe通道在不同中子能量点的裂变截面数据
裂变
核素
伽玛能量
(keV)
Iγ
截面σx(mb)
14.1±0.3
(MeV)
14.5±0.3
(MeV)
14.7±0.3
(MeV)
138Xe
153.9(E1)
0.0595
13.39±0.83
18.72±1.14
16.33±0.97
258.4(E2)
0.3150
13.95±0.77
16.89±0.94
15.68±0.84
396.5(E3)
0.0630
14.09±0.90
15.39±1.00
17.27±1.06
434.56(E4)
0.2032
42.81±2.35
38.47±2.13
41.49±2.22
如表6.2.3所示,通过对四条特征伽玛进行甄别和分析,我们得到了232Th(n,f)138Xe裂变通道在中子能量为14.1MeV,14.5MeV和14.7MeV时的裂变截面数据。
然而,对同一中子能量点不同特征伽玛,其结果有所差异,这主要是由于峰计数的统计误差造成的。
仔细分析数据可以看出,E4能量处的结果比E1、E2和E3处的截面值都要大。
原因在于,在E4处,由于混合了134Te特征γ射线能量435.06keV处的计数,而谱仪不能将两者甄别,造成在该能量峰的叠加。
故不能作为232Th(n,f)138Xe反应的截面值。
考虑到一些系统误差和统计误差的存在,并经过相关物理量的修正,计算所得到的138Xe不同特征γ峰能量处的截面值确定了232Th(n,f)138Xe反应的存在。
对同一中子能量点的三个数据取平均,得到中子能量为14.1MeV,14.5MeV和14.7MeV时,截面分别是13.81±0.83mb,17.00±1.03mb,16.43±0.96mb。
(2)232Th(n,f)130Sb截面分析
从核素衰变纲图可知,130Sb也可以由130Snβ-衰变而来,130Sn的半衰期3.72m,在特征γ能谱中没有发现130Sn存在,故没有其对该反应道的影响,可见130Sb的截面值较为准确。
表6.2.4232Th(n,f)130Sb通道在不同中子能量点的裂变截面数据
裂变
核素
伽玛能量
(keV)
Iγ
截面σx(mb)
14.1±0.3
(MeV)
14.5±0.3
(MeV)
14.7±0.3
(MeV)
130Sb
182.3(E1)
0.65
1.73±0.097
1.56±0.089
1.76±0.099
330.9(E2)
0.78
1.81±0.102
1.89±0.108
2.10±0.118
468.0(E3)
0.18
2.74±0.170
3.74±0.230
3.06±0.190
793.4(E4)
1
1.35±0.079
1.67±0.099
1.83±0.106
839.5(E5)
1
1.43±0.084
1.95±0.115
2.16±0.124
如表6.2.4所示,通过对五条特征伽玛进行甄别和分析,我们得到了232Th(n,f)130Sb裂变通道在中子能量为14.1MeV,14.5MeV和14.7MeV时的裂变截面数据,可以明确裂变碎片的存在。
然而,对同一中子能量点不同特征伽玛,其结果有所差异,这主要是由于峰计数的统计误差造成的,也可以看出能量为E3这组的截面值明显比其他组的值大,可能是因为与其他核素峰的重叠造成,所以对截面值分析时,舍弃这个峰。
对同一中子能量点的四个数据取平均,得到中子能量为14.1MeV,14.5MeV和14.7MeV时,截面分别是1.58±0.09mb,1.77±0.10mb,1.97±0.11mb。
(3)232Th(n,f)141Ba截面分析
由衰变纲图有,141Cs有到141Ba的β-途径,要准确计算141Ba的裂变截面,就要先扣除其先驱核141Cs通过β-衰变对计数的贡献。
但由于141Cs的半衰期只有24.94s,在裂变生成141Cs后,它在很短的时间内消失,实验使用的活化法只能甄别出中等半衰期的放射性核素,所以141Cs的存在与否,不能确定,如果它存在,则对本反应道的截面计算有一些影响,在认为141Ba不受141Cs的影响的情况下进行截面值分析。
表6.2.5232Th(n,f)141Ba通道在不同中子能量点的裂变截面数据
裂变
核素
伽玛能量
(keV)
Iγ
截面σx(mb)
14.1±0.3
(MeV)
14.5±0.3
(MeV)
14.7±0.3
(MeV)
141Ba
190.30(E1)
0.46
16.60±0.90
17.13±0.94
19.03±1.03
276.95(E2)
0.23414
55.48±3.00
59.76±3.25
61.64±3.34
304.20(E3)
0.25438
18.86±1.04
20.13±1.12
21.28±1.17
343.70(E4)
0.14444
16.61±0.94
18.65±1.06
16.49±0.94
462.30(E5)
0.04876
34.20±2.00
29.41±1.80
37.12±2.18
如表6.2.5所示,通过对五条特征伽玛进行甄别和分析,我们得到了232Th(n,f)141Ba裂变通道在中子能量为14.1MeV,14.5MeV和14.7MeV时的裂变截面数据。
然而,对同一中子能量点不同特征伽玛,其结果有所差异,这主要是由于峰计数的统计误差造成的,仔细观察数据,可以发现全能峰能量为E2和E5得到的这两组截面值较其它几个峰的值大很多,认真分析发现这两个峰的位置处有较大的本底峰,造成计数的叠加,故舍弃这两个峰。
从其余三组截面值看,数据符合较好,证实裂变碎片141Ba的存在。
对同一中子能量点的三个数据取平均,得到中子能量为14.1MeV,14.5MeV和14.7MeV时,截面分别是17.36±0.96mb,18.64±1.04mb,18.94±1.05mb。
(4)232Th(n,f)134Te截面分析
与上面的情况一样,134Te的先驱核134Sb,它的半衰期远小于1分钟,而实验中样品照射后的冷却时间是5分钟左右,根本无法测量得到134Sb的活化能谱。
表6.2.6232Th(n,f)134Te通道在不同中子能量点的裂变截面数据
裂变
核素
伽玛能量
(keV)
Iγ
截面σx(mb)
14.1±0.3
(MeV)
14.5±0.3
(MeV)
14.7±0.3
(MeV)
134Te
180.9(E1)
0.1829
14.77±0.806
15.09±0.831
16.24±0.864
210.4(E2)
0.22715
18.78±1.018
19.76±1.079
20.45±1.080
435.1(E3)
0.1888
21.57±1.178
18.93±1.047
21.03±1.122
767.2(E4)
0.295
16.84±0.926
17.00±0.945
17.65±0.947
如表6.2.6所示,通过对四条特征伽玛进行甄别和分析,我们得到了232Th(n,f)134Te裂变通道在中子能量为14.1MeV,14.5MeV和14.7MeV时的裂变截面数据。
然而,对同一中子能量点不同特征伽玛,其结果有所差异,这主要是由于峰计数的统计误差造成的。
对同一中子能量点的四个数据取平均,得到中子能量为14.1MeV,14.5MeV和14.7MeV时,截面分别是17.99±0.98mb,17.70±0.98mb,18.84±1.00mb,可见232Th(n,f)134Te裂变通道是确实存在的。
(5)232Th(n,f)134I和232Th(n,f)138Cs截面分析
甄别得到的碎片138Cs和134I,首先在不考虑被其它因素的影响下,按照截面计算公式得到了它们的截面值。
但这两个反应道的截面较其它几个反应道不一样,见如下两个核素图。
图6.2.1138Cs的核素图
图6.2.2134I的核素图
由上面两幅图可以看出,138Cs和134I的截面计算,会受到先驱核β-衰变和同质异能态退激的影响,且138Xe和134Te两个核素是存在的,基于已经计算得到138Xe和134Te的截面值,再利用相关公式[24]进行扣除。
在对两个碎片截面的详细分析中,首先运用6.1.2中给出的理论计算公式扣除核素138Cs和134I的截面中138Xe的β-衰变以及134Te的β-衰变对截面值的影响,结果计数很小,截面值极小(中间计算数据不在文中一一列举),此时就不再需要考虑同质异能态退激的影响。
说明,裂变本身产生极少量的138Cs和134I,它们主要来自于138Xe和134Te的β-衰变。
6.3比较分析
在核素甄别后对130Sb,138Xe,141Ba和134Te几个核素的裂变生成截面进行了分析,使得对232Th的裂变有了更深的了解,有的核素虽然可以在能谱上得到,但不能说明是完全由裂变反应产生它,有可能是有别的核素转变而来,比如β-衰变,以及同质异能态的退激,由上面对138Cs和134I的截面计算结果可以证明。
其中130Sb的截面值最为准确,因为它没有先驱核的影响,而其它几个核素是可能会有β-影响的,所以它们的截面值没有太大的意义。
在数据的获取过程中,γ谱上的许多峰存在重叠,交叉,本底的影响,使得在读取峰计数是存在较大的误差,所以这种峰不能采用,导致截面的计算上可用的峰减少,致使核素的裂变生成截面会有一些误差。
实验中使用了不同能量的种子照射样品,使用了14.1MeV,14.5MeV和14.7MeV三个能量,从计算结果可以看出能量高一点的中子所导致的裂变截面略微大一点,说明在14MeV能区中能量较高的中子对232Th活化程度深一些,但由于仅只有这几个通道的截面值,且截面值可能会有一些误差,所以这样的规律仅只是一个片面的推测,在今后大量的工作后就可以得到14MeV能区中子诱发裂变截面值随能量变化的关系。
通过对上面几组截面值的比较可以发现,130Sb的截面较其它几个碎片截面值小很多,这种情况可以由下图解释,在裂变产额计算时可以发现,产额随着质量数的变化类似于一个马鞍峰[25],在质量数为130时,裂变产额很小,所以130Sb的裂变生成截面很小。
而其它几个核素可以从图中看出,它们的裂变产额较大,通过计算得出的截面也较大。
图6.3.1裂变产额随质量数的变化
把上述的计算结果与235U的中子活化相对应的碎片裂变截面值相对比,可以发现,235U的裂变截面大得多,说明14MeV中子对235U的活化程度更高,也可以说明232Th比235U较为稳定。
从这些反应道截面分析中,可以发现,其中被甄别出的全能峰能量都很低,而那些能量较高的峰因为其分支比小,所以计数较低,即便峰存在也大多被本底湮没。
γ探测器是实验中的关键部分,如何提高探测器的探测效率和降低本底是截面精准计算的前提。
实验的一些主要误差来源有:
特征γ射线计数统计误差,特征γ射线探测效率误差,标准截面误差,样品称重误差,γ射线自吸收校正误差,测量几何校正误差,级联符合校正误差和中子注量率波动因子误差。
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 核素 裂变 截面 分析
![提示](https://static.bdocx.com/images/bang_tan.gif)