《1 前言》

1 前言

动物实验和临床试验都证明硼中子俘获治疗 (boron neutron capture therapy,BNCT)过程中血液中10B 的含量和肿瘤组织以及正常组织中10B 的含量成一定比例[1] ,所以通过测量血液中硼的含量可以推算出肿瘤组织和正常组织中硼的含量,且10B 产生的剂量是患者接受总剂量的主要组成部分,因而准确测量10B 的含量是准确计算患者所接受剂量的前提。

瞬发 中子活化分析(prompt gamma neutron activation analysis,PGNAA)方法利用中子束轰击血液样品,样品中10B 俘获中子后生成激发态 7Li *,其在 10-14 s内退激,释放出能量为 0.478 MeV 的瞬发 射线,通过高纯锗(HPGe)探测器探测瞬发 射线,根据特征峰的能量和强度(峰面积)对10B 浓度定量分析。该方法具有测量精度较高、速度快的特点,可以实现 BNCT 治疗过程中在线快速测量血硼浓度。国际上日本 KUR[2] ,美国 MIT[3,4] ,瑞典的 SINQ[5~7] ,韩国的 SNU-KAERI[8]等都已建立了 PGNAA 系统用以测量10B 及其他元素含量。

医院中子照射器是专门设计用于 BNCT 的中子源装置,其共有 3 条中子束流:热中子束、超热中子束和分析束,其中热中子束和超热中子束主要用于患者治疗照射,分析束是专门为用 PGNAA 方法测量硼浓度而设计建造的。笔者等根据医院中子照射器测量束实际参数,初步设计医院中子照射器 PGNAA 系统,并进行了硼浓度测量实验。

《2 医院中子照射器 PGNAA 系统初步设计》

2 医院中子照射器 PGNAA 系统初步设计

《2.1 系统组成》

2.1 系统组成

医院中子照射器 PGNAA 系统主要分为两部分:HPGe 探测器及其电子学系统和屏蔽系统。笔者等利用现有的 HPGe 谱仪及电子学系统,主要进行了屏蔽系统的设计,其中包括中子束屏蔽、探测器屏蔽和中子捕集装置 3 部分。

医院中子照射器 PGNAA 系统主要用于测量微量硼 -10 元素含量,需要低本底的测量环境。测量时本底主要来自 3 个方面:a. 从中子引出孔道周边缝隙的泄露和从地板和墙壁散射的中子及二次射线形成的环境本底;b. 热中子束上的中子和射线通过在样品上反应及分析室中散射产生的中子及射线进入探测器形成本底;c. 中子与血液中的杂质元素作用发生(n , )反应形成本底,如 23Na(n , 24Na 反应放出能量为 0.475 MeV 的瞬发射线。屏蔽系统设计主要考虑减少前两种本底。

因此医院中子照射器 PGNAA 系统屏蔽系统三部分的主要作用是:a. 中子束屏蔽装置:准直中子束,保证到达样品处的中子注量率尽量大,同时减少散射中子;减少反应堆开启状态下分析室漏射中子及二次射线进入测量系统。b. 探测器屏蔽装置:减少散射中子对 HPGe 探测器的损伤,同时屏蔽散射射线。c. 中子捕集器:吸收束上中子,减小中子散射,同时降低测量室内剂量水平,减小操作人员受照剂量。

《2.2 材料选择》

2.2 材料选择

 由前面分析可知:PGNAA 屏蔽装置设计主要考虑吸收散射中子和二次 射线。

医院中子照射器分析束为热中子束,束孔直径 2 cm,热中子通量密度设计值 ( < 0.4 eV ) = 3.08 × 106 cm-2 · s-1 ,方向性  =96.8 %。 因此设计时以屏蔽热中子为主。

常用热中子屏蔽材料一般含有硼元素、锂元素或氢元素。虽然10B 的热中子截面(3 838 b)比6Li(942 b)大,天然丰度也大于6Li,但是10B 为待测元素,会直接影响测量结果;由于测量时会使用硼/氢光峰比,所以屏蔽材料中不能含有 H 元素。因此选用含锂材料作为该系统热中子吸收材料。

常用的含锂材料有:含锂聚乙烯,LiF 粉末和 Li2CO3 粉末。含锂聚乙烯易机械加工,但因其含有H元素,不能使用在本系统中;Li2CO粉末可以溶于水,且 Li 元素质量百分比低于 LiF,因此选用 LiF。此外由于含 6Li 材料的购买受到严格控制且造价高 ,最终选择天然 LiF 作为热中子吸收材料。

对于二次射线的屏蔽,测量室内射线主要来自于墙面重混凝土屏蔽层漏射,混凝土漏射射线主要来自于混凝土中的氢俘获中子产生 2. 224 MeV 的射线以及重混凝土中的铁俘获中子产生 2~3 MeV 的射线。采用低本底测量通用的 10 cm 厚铅屏蔽层屏蔽此能量下的射线。

《2.3 系统设计》

2.3 系统设计

首先根据已有探测器情况和测量室实际情况设计样品位置,探测器摆放方向等。HPGe 探测器冷指为水平方向,因此设计探测器水平中心高度与束流中心高度一致,且垂直于中子束流。样品摆放于位于中子束流上,探测器正前方,应尽量靠近中子束流孔道出口。

中子吸收材料为天然 LiF 粉末,其实际封装密度 ρ =1.0 g/cm3。对于中子束和探测器屏蔽装置,采用厚度为 3 cm 的 LiF 粉末,热中子通量下降倍数为 1.59E+02;对于中子捕集装置,采用厚度为 20 cm 的 LiF 粉末,热中子通量的下降倍数为 4.67E+14 。 所有 LiF 粉末均封装于厚度为 3 mm 的聚四氟乙烯盒内。

图 1 为医院中子照射器 PGNAA 系统初步设计示意图。样品摆放位置距束流出口 45 cm , 距探测器铍窗表面 6 cm,水平高度 110 cm 。中子束屏蔽系统中心为直径 3 cm,壁厚 3 mm 的圆柱形空心聚四氟乙烯管,长度为 60 cm 。管外为最小厚度为 3 cm 的 LiF 粉末 , 封装于长方体聚四氟乙烯盒内 , 横截面为边长是 12 cm 的正方形。中子捕集装置总长 20 cm , 也由聚四氟乙烯盒组成 , 横截面也为边长是 12 cm 的正方形。最外层为 10 cm Pb 屏蔽层。探测器屏蔽系统结构与束屏蔽系统相同,也为双层屏蔽结构,但探测器前端增加厚度为 2 mm 的 6Li2CO3 粉末,粉末封装于 Al 盒内。图 2 为屏蔽系统搭建完成后的照片。 

《图1》

图1 医院中子照射器 PGNAA 系统初步设计图

Fig.1 Diagram of preliminary designed PGNAA system of IHNI

《图2》

图2 按初步设计方案搭建的 PGNAA 屏蔽系统照片

Fig.2 The photo of PGNAA system as the design

《3 硼浓度测量实验》

3 硼浓度测量实验

《3.1 实验方法》

3.1 实验方法

当热中子束轰击样品时,样品中的10B 发生10B (n , α)7Li * 反应,7Li * 在 10-14 s 内退激。一般血液 样品中10B 和 H 分布均匀 ,中子轰击样品后,10B (n , α)7Li 反应和1H(n , 2H反应所产生的瞬发 射线计数之比与中子束的参数(中子注量率、中子能谱和中子束均整性等)无关。所以10B(n , α)7Li * 反应产生瞬发射线的计数与1H(n , 2H 反应所产生的瞬发射线的计数之比正比于样品中10B 和 H 的浓度之比。待测样品的10B 可由标准样品的10B/H 光峰比推出。该方法属于相对测量方法,用该方法测量未知样品中的血硼浓度时结果与中子束参数无关,与 HPGe 探测器的探测效率无关。对样品的量以及样品容器的几何尺寸无严格要求,但样品的量会影响到测量结果的统计误差。文章使用样品由标准硼酸溶液配制。

《3.2 实验装置》

3.2 实验装置

实验使用 HPGe 射线探测器型号为 Ortec - GMX35P4,高纯锗晶体尺寸为 58. 8 × 80. 0 mm3,Be 窗厚 0.5 mm 。实测 60Co 1. 332 MeV 特征峰半高宽 FWHM = 1. 82 MeV,峰康比 65∶1。图 3 为 HPGe 探测器电子学线路框图。屏蔽系统如上所述。标准样品由硼标准溶液加去离子水配置而成。样品容器为圆柱形聚四氟乙烯瓶。由于束流出口直径为 2 cm,因此样品溶液高度应小于 2 cm , 瓶直径小于 2 cm , 因此样品容器直径为 16 mm , 高 30 mm 。 同时考虑尽量增加计数率因此实验选用 2 mL 样本进行测量。由血液样本中添加不同质量硼酸得到含硼血液。

《图3》

图3 HPGe 探测器电子学线路框图

Fig.3 Electronic circuit diagram of HPGe detector

《3.3 实验结果》

3.3 实验结果

测量样品时,反应堆功率 30 kW,分析室 HPGe 探测器一侧环境 剂量率 <7 μSv/h,中当量剂量率 <3 μSv/h。图 4 为满功率情况下测量 50 ppm 10B, 测量时间 600 s,射线全谱。

《图4》

图4 满功率情况下测量 50 ppm 10B 2 mL 样品600 s得到射线全谱图

Fig.4 -ray full spectra of 50 ppm 10B 2 ml sample with measuring time 600 s for the reactor at full power

实验主要测量了不同10B 浓度样本。图 5 为测量硼 -10 浓度为10 ppm,50 ppm标准样品,水和空瓶时 450~550 keV 能量段能谱。10B 多普勒展宽峰能段计数由多到少依次为 50 ppm 标样、10 ppm 标样、水样品和空样品瓶。

《图5》

图5 满功率情况下测量不同10B 浓度 2 mL 样品 300 s 得到的能谱图

Fig.5 -ray spectra of 2 mL sample of different 10B concentrations with measuring time of 300 s

该能谱表明使用该系统可以分辨10B 浓度为 10 ppm,体积为 2 mL 的样品。但同时也发现在只放置空瓶和样品为水时依然可以探测到较高的 10B 0.478 MeV 多普勒展宽峰 ( 0.470~0.485 MeV)。此外反应堆满功率运行情况下还测量了 2 mL 血液样品,如图 6 所示。含 50 ppm 10B 血液样本与 50 ppm 标样测量结果基本一致。血液中含有 Na,会发生 23Na(n,24Na反应,产生能量为 475 keV 的光子,但其热中子截面较小,该能量光峰湮没于本底中。单纯血液样本计数明显少于 50 ppm 10B 血液样本与 50 ppm 标样。

《图6》

图6 满功率情况下测量 2 mL 血液样品及标样 300 s 得到的能谱图

Fig.6 -ray spectra of 2 mL blood sample and standard sample with measuring time of 300 s

反应堆满功率 30 kW 运行,50 ppm 10B 样品净计数率为 41. 578 ± 0. 658 cps;2mL,10 ppm 10B 样品净计数率为 8. 223 ± 0. 38 6cps。10B 元素测量灵敏度为 0. 822 cps/ppm。

《4 讨论》

4 讨论

根据测量结果可知,10B 光峰能段本底较高,该本底可能产生的原因有:a. 系统所用材料中杂质 10B;b. 除10B 外其他元素产生激发态 7Li * 其退激产生 0. 478 MeV多普勒展宽峰。

如图 7 所示,将系统中所用到材料做成样品放入系统中测量,由图可知 LiF 样品产生计数最多,铝样品和空样品瓶在本能段计数基本相同且均少于 LiF 样品。 LiF 屏蔽系统中子吸收体正为此种材料。根据7Li 各反应道反应截面 , 中子能量大于 0. 5MeV 时会发生7Li(n , n′)7Li反应,反应截面小于 1 b , 该反应产生的激发态 7Li退激同样产生 0. 478 MeV 的射线。

《图7》

图7 满功率情况下测量不同样品 600 s 得到的能谱图

Fig.7 -ray spectra of samples of different materials with measuring time of 300 s

虽然 Li ( n , n′ ) 7Li反应截面较 10B ( n , α ) 7Li截面小 3~4 个量级,但由于分析束根据计算值能量大于 0.5 MeV 中子的通量密度约 6.1 × 10cm-2 · s-1 ,且 LiF 用量很大,样品室体积较小,因此该反应造成影响不可忽略。

如果在吸收中子的同时又要避免使用硼和锂作为吸收体,而镉等中子吸收材料在吸收中子的同时会产生大量的光子,使得探测器的死时间过大,不利于测量。因此考虑在 LiF 中子吸收层内侧添加屏蔽层,该材料应具备以下 3 个特点 : a. 中子截面小;b. 射线屏蔽能力强;c. 与中子反应不产生 0. 478 MeV 能区附近的二次射线。此外,应考虑增大样品室体积,减少散射中子及屏蔽材料产生二次射线进入探测器的几率。