1、 1 中华人民共和国国家环境保护标准 HJ 1129-2020 就地高纯锗谱仪测量土壤中 核素技术规范 Technical Specifications for Determination of -emitting Radionuclides in Soil by In-situ HPGe Spectrometry （发布稿） 2020-06-03 发布 2020-06-30 实施 生态环境部 发 布 2 目 次 前 言 . 3 1 适用范围 . 4 2 规范性引用文件 . 4 3 术语和定义 . 4 4 方法原理 . 6 5 仪器和设备 . 6 6 测量分析步骤 . 6 7 探测效率计算 .
2、 7 8 结果计算与表示 . 9 9 质量保证和质量控制 . 10 10 注意事项 . 11 附录 A （资料性附录） 就地高纯锗谱仪的视域 Rs . 12 附录 B （资料性附录） 质量张弛深度 的确定方法和经验数据 . 14 附录 C （资料性附录） 放射性核素的注量活度比 F . 16 附录 D （资料性附录） 标准 点源校准就地高纯锗谱仪有效前面积 及角度修正因子 W . 36 附录 E （资料性附录） 蒙特卡洛模拟方法计算就地高纯锗谱仪的有效前面积 和角响应修 正因子 W . 38 附录 F （资料性附录） 就地高纯锗谱仪有效前面积 及角度修正因子 W 经验值 . 41 附录 G （
3、资料性附录） 就地高纯锗谱仪 测量土壤中 核素活度浓度的影响因素 . 43 3 前 言 为贯彻中华人民共和国环境保护法中华人民共和国放射性污染防治法中华人 民共和国核安全法，保护生态环境，保护公众健康， 规范就地高纯锗谱仪测量 土壤中 核素 放射性活度浓度和放射性沉降灰沉积通量的 方法，特制定本标准。 本标准的附录 A附录 G为资料性附录。 本标准为首次发布。 本标准由生态环境部核设施安全监管司、法规与标准司组织制订。 本标准起草单位：生态环境部核与辐射安全中心、生态环境部辐射环境监测技术 中心、 华北电力大学、安徽省辐射环境监督站。 本标准生态环境部 2020年 06月 03日批准。 本标准
4、自 2020年 06月 30日起实施。 本标准由生态环境部解释。 4 就地高纯锗谱仪测量土壤中 核素技术规范 1 适用范围 本标准规定了 用于测量土壤中 核素放射性活度浓度和放射性沉降灰沉积通量的就地 高纯锗谱仪测量方法 。 本标准适用于核事故应急监测时地表土壤人工 放射性 核素 的测定， 辐射环境质量监测 可参照执行。 对于新沉降的放射性沉降灰，本标准测定的 137Cs 活度浓度探测下限为 100 Bq/m2。 2 规范性引用文件 本标准引用了下 列文件或其中的条款。凡是不注日期的引用文件，其有效版本适用于本 标准。 GB/T 27418 测量不确定度评定和表示 ISO 18589-7 环境
5、放射性的测量 土壤 第 7 部分： 核素的就地测量（ Measurement of radioactivity in the environment - Soil - Part 7: In situ measurement of gamma-emitting radionuclides） 3 术语和定义 下列术语和定义适用于本标准。 3.1 相对探测效率 relative detection efficiency 在点源至探测器前表面距离为 25cm 时，高纯锗探测器和标准圆柱形 NaI(Tl)闪烁体探测 器（直径 7.62 cm，高度 7.62 cm）测量 60Co 源 1332.5 keV
6、 射线的全能吸收峰探测效率之 比。 注：标准圆柱形 NaI(Tl)闪烁体晶体对 1332.5 keV 射线的全能吸收峰探测效率理论值为 1.210-3。 3.2 能量分辨力 energy resolution 探测器分辨能量不同却又非常相近的入射 射线的能力。能量分辨力与入射 射线能量 有关，对于指定能量的单能 射线 ，常用该能量的全能吸收峰的半高宽来表示， keV。 3.3 半高宽 full width at half maximum 在脉冲高度谱中，单峰峰值一半处两点的横坐标之间的距离， keV。半高宽又记作 FWHM。 注：如果能谱中的峰由几个峰叠加组成，则每个峰都有一个半高宽。 5 3
7、.4 探测器的视阈 field-of-view of a detector 表征探测器测量地表土壤范围的量，即探测器测量到 90%的非散射 射线来源于该范围 内的土壤。一般用测量土壤范围的半径表示， m。 3.5 探测下限 lower limit of detection （ LLD） 在给定的置信度下，就地高纯锗谱仪可探测到最低的活度浓度值或最小沉积通量。 3.6 就地测量 in-situ measurement 不对待测对象做任何物理或化学性质处理或改变，在其所处的现场进行的测量。 3.7 探测器高度 detection height 高纯锗晶体几何中心距地表的垂直高度， m。 3.8 活
8、度浓度 activity concentration 单位质量待测物质中放射性核素的活度， Bqg-1。 3.9 沉积通量 sedimentation flux 放射性沉降灰在单位面积上沉降的特定放射性核素的活度， Bqm-2。 3.10 就地高纯锗谱仪的探测效率 detection efficiency of in-suit HPGe Spectrometry 探测器高度为 1m 时，就地高纯锗谱仪测量能量为 E 的特征 射线计数率与土壤中发射 该射线的放射性核素的沉积通量或活度浓度之比， s-1(Bq.m-2)-1 或 s-1(Bqg-1)-1 。 3.11 有效前面积 effective
9、 front area 表征探测器全能吸收峰探测效率的参数，即能量为 E 的 射线束沿探测器轴向平行入射 的情况下，单位注量率所产生的全能吸收峰计数率， m2。 3.12 角 度修正因子 angle correction factor 对于能量为 E 的 射线，就地高纯锗谱仪在野外地表实际测量时与测量沿探测器轴向平 行入射时，单位 射线注量率所得全能吸收峰计数率之比，量纲一的量。 3.13 注量活度比 fluence activity ratio 土壤中单位活度浓度或单位沉积通量的放射性核素发射能量为 E 的射线，在地面上 方 h 高度处的初级射线注量率， (m-2s-1) (Bqg-1)-1
10、 或 (m-2s-1) (Bqm-2) -1。 3.14 张弛深度 relaxation depth 描述放 射性核素随土壤深度变化的数学参数，表示活度浓度减少至表面活度浓度的 e-1 （ 37%）时，距离土壤表面的深度， cm。 6 3.15 质量张弛深度 relaxation mass per unit area 地表至待测放射性核素的张弛深度范围内，单位面积土壤的质量， gcm-2。 3.16 等效表面沉积通量 equivalent surface sedimentation flux 对于放射性沉降灰中的核素，假设质量张弛深度为 0 gcm-1 时 ，计算得到的沉积通量即 为等效表面沉
11、积通量， Bqm-2。 4 方法原理 地表射线的注量率与土 壤中放射性核素的活度浓度（沉积通量）、核素在土壤中的深 度分布，以及土壤的组分、含水率等因素有关。探测器高度为 1 m 时，就地高纯锗谱仪测量 的全能吸收峰计数率与土壤中放射性核素的活度浓度成正比。 fNA （ 1） 式中： A 待测放射性核素的活度浓度或沉积通量， Bqg-1 或 Bqm-2； Nf 待测放射性核素的全能吸收峰净计数率， s-1； 就地高纯锗谱仪的 探测效率 ， s-1(Bqg-1)-1 或 s-1(Bqm-2)-1。 5 仪器和设备 5.1 高纯锗谱仪：由高纯锗探测器、高压电源、主放大器、多道脉冲幅度分析器、制冷单
12、 元（液氮或电制冷）、数据处理系统等 6 个主要部分构成。 5.2 支架：能够固定探测器，使探测器高度距地面 1 m 高的支架，应使用低原子序数、低 密度的材料制造（如：铝、塑料、木材）。 5.3 卷尺 ：卷尺的校准长度应大于 1 m，且校准分度不小于 1 mm。 5.4 其它辅助设备：应配备地理定位系统，记录测试时的地理坐标。 6 测量分析步骤 6.1 布点原则 一般情况，选择一个开阔、平坦的区域，而且探测器的视阈 Rs 内应无干扰物体（如建 筑物、树林 等），探测器的视域 Rs 见附录 A。测量放射性核素的沉积通量时，应选择自放 射核素沉积以来，放射性核素在土壤中分布未被破坏的区域（如：没
13、有农业活动）。 6.2 测量条件记录 记录测量点位地理坐标，以及 气象条件、 地形地貌的描述等 。 7 6.3 仪器状态的检查 测量前，应检查高纯锗探测器的晶体温度是否达到工作温度。处于工作温度后，再打开 高纯锗谱仪的偏压电源，检查仪器的能量分辨力。 6.4 就地高纯锗谱仪的安装 在选定的测量区域中间位置，将高纯锗探测器稳定地安装在测量支架上，探测器的端面 朝下。调整探测器高度，使之距离地表 1 m。连接高纯锗谱仪电子 学仪器的所有电缆。 6.5 测量设置 根据现场测量的待测核素特征峰的全能吸收峰计数率，设置测量时间，使得感兴趣区的 计数统计涨落在 10%以内。 6.6 能谱的采集 在预定的时
14、间内采集 能谱。停止测量后，存储采集到的 能谱。 6.7 能谱分析 用 能谱分析软件查看采集的 能谱，若能谱中存在可明显辨别的峰，根据峰形特征设 置适合的感兴趣区域，读取全能吸收峰的净峰面积；若能谱中无可明显辨别的峰，则根据待 测核素的特征峰的能量和探测器的能量分辨力，设置感兴趣区，读取全能吸收峰的积分面积。 7 探测效率计算 7.1 原理 对于就地高纯锗谱仪而言，很难有已 知不同深度分布的大面积 射线标准场用于探测效 率校准。因此，需根据土壤、空气的物理参数，采用放射性核素在土壤中分布的理论模型， 并结合探测器的特征进行探测效率计算。就地高纯锗谱仪的探测效率取决于放射性核素在土 壤中的深度分
15、布，探测器高度，探测器的性质， 射线的能量以及在土壤和空气的吸收性质 等因素，可利用如下 Beck 公式进行计算： WF （ 2） 式中： F注量活度比， (m-2s-1)(Bqg-1)-1 或 (m-2s-1)(Bqm-2)-1。该物理量由放射性核素在土壤 中分布、土壤组分与密度等因素决定，与探测器无关。 W角度修正因子，由探测器的角响应因子 R()和地表 射线注量率角分布 ()共同决 定的，量纲一的量。这是因为探测器非标准球体，对不同角度入射的 射线探测效率并不相 同。而探测器所在的位置， 射线是从土壤表面发射而来的（ 2 立体角），而非沿探测器轴 向平行入射的 射线，故需要进行角度修正。
16、该物理量与探测器几何结构、 射线能量、放 射性核素在土壤的分布有关。 8 有效前面积， m2。该物理量只与所用的探测器 有关，由探测器晶体尺寸、材料、几 何结构等因素决定，与土壤中放射性核素分布无关。 7.2 探测 效率 计算 的方法 7.2.1 方法概述 由 7.1 可知，由注量活度比 F、角度修正因子 W 和有效前面积 这三个不相关的物理 量乘积可得就地高纯锗谱仪的探测效率 。其中，注量活度比 F 和角度修正因子 W 与放射 性核素的质量张弛深度 相关，确定质量张弛深度 是开展探测效率计算的前提条件。在已 知质量张弛深度 条件下，按照以下步骤开展探测效率计算：（ 1）根据理论模型计算土壤
17、中放射性核素的注量活度比 F。 （ 2）利用标准源实验测量或蒙特卡洛模拟计算探测器的有 效前面 积 。（ 3）理论模型计算探测器高度处初级 注量率的角分布 ()；利用标准源实 验测量或蒙特卡洛计算探测器的角度响应因子 R()；联合初级 注量率的角分布 ()和探测 器的角度响应因子 R()计算得到角度修正因子 W。 7.2.2 质量张弛深度 的确定 一般来说，新近沉降的放射性沉降灰，放射性核素均匀分布在土壤表面， = 0 gcm-2； 天然放射性核素在土壤中呈均匀分布的， = gcm-2。放射性沉降灰中 核素随着时间推移 不断向土壤深度渗透，活度浓度随着土壤深度 z 的增加而呈指数下降， 0 g
18、cm-2 gcm-2。 若质量张弛深度 未知， 可先假设 = 0 gcm-2，计算等效表面沉积通量。之后，根据现场 分层取样或经验值得到实际的质量张弛深度 ，质量张弛深度 的确定方法和经验数据见附 录 B。 7.2.3 注量活度比 F 的计算 待测放射性核素的质量张弛深度 确定后，可由以下三种方法得到注量活度比 F：（ 1） 对于沉降在土壤的人工放射性核素（ 0 gcm-2 gcm-2）查阅附录 C 中表 C2 得到注量 活度比 F；对于天然放射性核素（ = gcm-2）查阅附录 C 中 C3 得到注量活度比 F。（ 2） 对于沉降在 土壤的放射性核素（ 0 gcm-2 gcm-2）， 若质量
19、张弛深度 在附录 C 表 C2 中无对应值时，根据质量张弛深度 的数据 插值计算得到注量活度比 F。（ 3）按照附录 C 的理论模型，计算得到注量活度比 F。 7.2.4 有效前面积 的校准 就地高纯锗谱仪的有效前面积 可通过以下三种方法得到：（ 1） 标准 点源校准方法， 具体的 方法及步骤 见附录 D。（ 2）蒙特卡洛模拟计算方法，具体的 方法及步骤 见附录 E。 （ 3）经验公式快速校准方法，根据探测器的基本性质利用经验公式快速校准的方法见附录 F。 7.2.5 角度修正因子 W 的校准 就地高纯锗谱仪 的角度修正因子 W 可通过以下三种方法得到：（ 1） 标准 点源校准 方 法，具体的
20、方法及步骤 见附录 D。（ 2）蒙特卡洛模拟计算 方法，具体的方法及步骤 见附录 9 E。（ 3）对于新沉降的人工放射性核素（ = 0 g.cm-2）和天然放射性核素（ = g.cm-2）， 根据探测器的基本性质，参考典型高纯锗谱仪的角度修正因子进行校准，见附录 F。 8 结果计算与表示 8.1 结果计算 土壤中放射性核素的 沉积通量或 活度 浓度可用 式（ 3）计算 ： tNA net （ 3） 式中： A待测 放射性核素的沉积通量或活度浓度， Bqm-2 或 Bqg-1； Nnet 能量为 E 的 射线全能吸收峰净计数，量纲一的量； t有效 测量时间， s； 对于 能量为 E 的 射线，
21、就地高纯锗谱仪 的 探测效率 ， s-1(Bqm-2)-1 或 s-1(Bq.g-1)-1。 8.2 测量不确定度评估 根据 GB/T 27418 推荐的 不确定度传递传播公式，对公式（ 2）和（ 3）各不确定度分量 的来源进行计算。公式（ 3）中不确定度最大来源为 和 Nnet，其它量的不确定度可以忽略 不计，因此合成标准不确定度用下式计算 ： 2 2 2 2 2 2 2n e tn e t 2 )()()()()( F FuW WuuN NuAAu （ 4） 式中： u (A)待测 核素活度浓度或沉降通量 A 测量结果的不确定度； u(Nnet)能量为 E 的 射线全能吸收峰净计数 Nne
22、t 的不确定度； u()能量为 E 的 射线有效前面积 的不确定度； u(F)能量为 E 的 射线注量活度比 F 的不确定度； u(W)能量为 E 的 射线角度修正因子 W 的不确定度。 8.3 探测下限 在 95%置信度下， 就地高纯锗谱仪测量土壤中 核素的活度浓度 或沉积通量的 探测下限 LLD 用公式 （ 5）计算。 g4.6 5 NLLD t （ 5） 式中： Ng能量为 E 的 射线全能吸收峰的总计数，量纲一的量。 10 8.4 结果表示 根据报告要求，对于大于探测下限的核素，测量结果应以活度浓度或沉积通量及测量不 确定度 A u(A)来表述；对于小于探测下限的核素，测量结果应以“小
23、于 LLD”表述，并给 出探测下限 LLD 值。测量结果的小数点后位数的保留与方法探测下限一致，最多保留三位 有效数字。 9 质量保证和质量控制 9.1 影响变量 附录 G 中讨论了影响测量结果的变量，主要包括测量点位的环境特征、地表的平整度、 地表植被特征、土壤组分、探测器高度、有效前面积和角 度修正因子 的校准。测量时，需要 特别注意并尽可能地降低这些参数的影响。 9.2 仪器验证 按照实验室质量控制体系的要求，应定期对 高纯锗谱仪 的关键技术指标（如探测效率、 能量分辨力、 测量 本底）进行期间检查 。 9.3 质量控制计划 9.3.1 总则 就地高纯锗谱仪应建立质量保证计划，规定 质量
24、控制 程序和 质量 措施以保证测量结果的 质量。 9.3.2 期间核查 9.3.2.1 性能测试 将就地高纯锗谱仪放在固定位置，分析环境中天然放射性核素（如 40K， 214Pb， 214Bi， 226Ra） 特征峰的能量和峰形参数。定期检查高纯锗谱仪性能指标，以及测量结果是否可复 现。 9.3.2.2 能量分辨力和探测效率 应定期对就地高纯锗谱仪的能量分辨力和探测效率进行期间核查。一般情况，每年至少 开展一次期间核查。若高纯锗探测器经常回温，应缩小期间核查的时间间隔。用固定检验源 测量高纯锗谱仪的能量分辨力和全能吸收峰探测效率，绘制质控图。检验源在低能、中能和 高能区应至少各有 1 条 射线
25、。若期间核查结果与之前相比存在显著差异，需要进行重新校 准或维修。 9.3.2.3 本底 在固定测量条件下，应定期进行高纯锗谱仪的测量 本底进行期间核查。一般情况下，每 11 年至少开展一次期间核查，若怀疑有污染可能性的时候，也应尽快开展本底期间核查，绘制 质控图。若期间核查结果与之前相比存在显著差异，需要进行去污处理。 9.4 文件文档 就地高纯锗谱仪的测量过程应有详细的原始记录，以及现场测量的 能谱文件，并规定 适当的保存期限。一般而言，应急监测、辐射环境质量监测和重点污染源监测的原始记录应 永久保存，期间核查报告等质量保证记录应至少保存 6 年。 10 注意事项 10.1 就地高纯锗谱仪
26、必须 注意 防潮和防水。 10.2 就地高纯锗谱仪的工作环境温度在 -20 +50 范围内，否则可能使得高纯锗探测 器的真空变差。 10.3 制冷达到工作温度 后， 高纯锗晶体易碎， 在 运输过程中应保证机械冲击不会导致探测 器损坏。 10.4 高纯锗 谱仪 在工作温度才能工作，测量前需要预估探测器的运输和测量时间，保证有 足够的液氮制冷或电池电量。 注 :对于液氮制冷系统，高纯锗探测器在运输过程中，液氮消耗量会变大。 10.5 若放射性核素随雨水沉降时，因降水的流动，可能导致测量区域放射性核素分布并非 均匀，可能产生较大的偏差 ，不宜进行测量。 10.6 空气中 存在人工 放射性核素 时，须
27、 防止高纯锗谱仪被污 染，并注意监测人员的防护。 10.7 超出视阈范围的山脉、山脊也可能影响最终测量结果，建议尽量避免在周边开展测量。 10.8 测量放射性沉降灰时，应尽可能避开周围存在树林的测量环境。 12 附录 A （资料性附录） 就地高纯锗谱仪的视域 Rs 在 开阔、平坦的区域，且 探测器高度为 1 m 时，就地高纯锗谱仪的视域 Rs 见图 A1。 就地高纯锗谱仪的视域 Rs 与质量张弛深度 ， 射线的能量 E 有关。质量张弛深度 越小， 射线的能量 E 越高，就地高纯锗谱仪的视域 Rs 越大。 100 1000 10000 1 10 100 1000 1 10 100 1000 视域
28、 R s （ m ） E ( k e V) 0 g c m -2 0 .3 g c m -2 0 .5 g c m -2 1 .0 g c m -2 3 .0 g c m -2 1 0 .0 g c m -2 g c m -2 图 A1 不同质量 张弛深度 ， 就地高纯锗谱仪的视域 Rs 0 10 20 30 40 50 60 70 80 90 100 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1.0 比值 半径 ( m ) 121 k e V 662 k e V 898 k e V 1332 k e V 1461 k e V 图 A2 不同半径范围内土壤中放射性核素在地表上方
29、1m 高处的初级射线注量率与总注量率 13 的比值（ = 0 gcm-2） 0 10 20 30 40 50 60 70 80 90 100 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1.0 距离 比值 半径 ( m ) 121 k e V 662 k e V 898 k e V 1332 k e V 1461 k e V 图 A3 不同半径范围内土壤中放射性核素在地表上方 1m 高处的初级射线注量率与总注量率 的比值（ = gcm-2） 14 附录 B （资料性附录） 质量张弛深度 的确定方法和经验数据 B.1 分层取样法确定放射性核素质量张弛深度 a） 在就地
30、高纯锗谱仪的视域 Rs 内，采 用简单随机布点，确定采集土壤的点位。 b） 在每个测量点，按照等间距或随土壤深度增加间距变大（表层间距一般小于 1 cm），分 若干层采集土壤样品。 c） 将各个采集点相同土壤深度 Hi 的样品混合在一起，作为一个分层土壤样品，张贴样品标 识。 d） 在实验室内，将 土壤样品中杂草、碎石等异物剔除， 测量样品的密度 s（ gcm-3）。 e） 利用实验室高纯锗谱仪 分别测量各层 (序号为 i)土壤样品中 各放射性 核素的活度浓度 Ai （ Bqg-1）； f） 建立坐标系， 横坐标为 土壤深度 H(cm)，纵坐标为 各层土壤深度的活度浓度 Ai（ Bq.g-1）
31、 与表层土壤活度 浓度 Am,0（ Bqg-1）之比的自然对数 ln(Ai/Am,0)（量纲为一的量）。以各层 土 壤深度 Hi 及相应的 ln(Ai/Am,0)在建立的坐标系中作图。 g） 用最小二乘法 拟合直线 ，求出 ln(Ai/Am,0)和 Hi 的直线斜率 k（ cm-1）。 h） 放射性核素的质量张弛深度 （ gcm-2）用式 B-1 计算可得。 k- (B1) 式中， 放射性核素的质量张弛深度 ， gcm-2； 土壤的密度 ， gcm-3； k拟合的直线斜率 ， cm-1； B.2 137Cs 质量张弛深度 的经验数据 表 B.1 137Cs 质量张弛深度 的经验数据 a 序号
32、沉降时间 来源 测量地点 （ g.cm-2） 1 34 周 切尔诺贝利核事故 俄罗斯西部 01.2.0b 2 56 周 切尔诺贝利核事故 德国 0.51.0b 3 1 年 切尔诺贝利核事故 瑞典 2.20.2 4 1 年 切尔诺贝利核事故 德国 1.4 5 15 年 大气核试验 美国西部 144c 15 表 B.1 137Cs 质量张弛深度 的经验数据 （续 ) 序号 沉降时间 来源 测量地点 （ g.cm-2） 16 15 年 大气核试验 美国南部 2.91.6d 17 15 年 大气核试验 美国东部 27（森林） 89（平原） 18 30 年 大气核试验、切尔诺贝 利 中国吉林 11.2e
33、 a 序号 117 的数据来源于 ICRU Report 53 环境测量中的射线谱仪法（ Gamma-Ray Spectrometry in the Environment） b 湿沉降 c 干旱区 ，灌溉的草原 d 干旱区，自然保护区 e 长白 山自然保护区，高山苔原带 16 附录 C （资料性附录） 放射性核素的注量活度比 F 一般来说，新沉降的放射性核素在地表呈 面 分布， = 0 g.cm-2；天然放射性核素在土 壤中均匀分布的， = g.cm-2。随着时间推移 ， 放射性沉降灰中 的 放射性核素不断向土壤 深度渗透，活度浓度随着土壤深度 z 的增加而呈指数下降， 0 g.cm-2 1
34、MeV 的 射线，晶体的有效中心近似在晶体的几何中心； 2） 对 E0.1MeV 的 射线，晶体的有效中心近似在晶体的表面； 3） 对能量处在上面两个值之间的 射线，必须根据晶体的吸收系数作出平均穿透 厚度的估算： 01 1 (1 )=. 1- d dedr d xe （ D2） 能量为 E 的 射线在高纯锗晶体中的减弱系数， cm-1； d 高纯锗晶体的厚度， cm； d0 探测器前表面到晶体的距离， cm。 b） 采集并储存 能谱，进行解谱，确定校准源特征峰的全能吸收峰的计数率 ng。 c） 采集不放源时的本底谱，进行谱分析，获取本底计数率 nb。 d） 由修正后的计数率 (ng-nb)除
35、以注量率 ，得到 。 = gbnn （ D3） e） 对于其它校准源 ，重复 a） d）步，得到探测器不同能量的有效前面积 。 37 D.2 标准 点源校准 就地高纯锗谱仪角度修正因子 D.2.1 用标准 点源校准 探测器角响应因子 R() 校准高纯锗探测器角响应因子 R()的标准 点源应与 D.1 相同，一般包括 241Am、 133Ba、 137Cs、 60Co、 152Eu 等。 a）尽可能远离探测器（ 10 探测器直径）的固定长度放置标准点源，在入射角 = 0(垂 直于探测器的正表面 )和 = 90之间，每隔 10间隔测量净峰计数率 Nf()。 b） 用 = 0位置得到的计数率 N0
36、去归一各角度的净计数率 Nf()，得到相对响应 Nf()/ N0 与角度 的关系，然后用多项式进行拟合。 0( ) cos n i iiRk （ D4） 式中： n多项式的阶数， 量纲一的量，一般情况下采用四阶多项式，即 n = 4； ki拟合系数 ，量纲一的量。 D.2.2 角度修正因子 W 计算 对于能量为 E 的 射线，高纯锗探测器的 角度修正因子 W 用公式（ D5）进行计算。 dRW )()(2/ 0 （ D5） 式中： () 射线注量率的角分布，可以根据放射性核素在土壤中的分布计算出来，详见 附录 C； 弧度， rad。 能量不同的的 射线 重复步骤 D.2.1 和 D.2.2 ，
37、得到不同能量的角度修正因子 W。 38 附录 E （资料性附录） 蒙特卡洛模拟方法 计算就地高纯锗谱仪的 有效前面积 和角 度修正因子 W 蒙特卡洛方法能够准确地模拟粒子在物质中的输运过程，在处理复杂几何条件具有很大 的优势。本附录提供了蒙特卡洛模拟计算高纯锗探测器的有效前面积和角响应的方法。 蒙特卡洛模拟计算的准确度取决于仿真计算模型与实际情况的一致性。因此，在构建几 何模型前，用户需要获取探测器的准确材料组成和几何尺寸。一般情况下，用户虽然无法获 得确切的探测器组件材料和几何尺寸，但能够根据所有已知的探测器组件信 息和实验测量数 据，构建一个合适的几何模型。 E.1 高纯锗探测器 建模的准
38、确度评估 以下步骤过程用于验证高纯锗探测器建模的准确度。 E.1.1 实验测量 a）仪器与设备 1）高纯锗谱仪：见 5.1.1 所述的高纯锗谱仪。 2）标准 点源：实验所用的标准 点源都应可溯源至国家基准。标准点源的特征 射 线应至少涵盖低能、中能和高能区间，一般推荐使用 241Am 和 152Eu 标准点源。 3）支架：用于将放射源固定在固定位置的支架。该支架可在 0、 45和 90的方位角上 固定，而且放射源与探测器之间的距离可准确获得。 b）实验 测量 1） 将标准 点源固定在高纯锗探测器的支架上，分别在 0、 45和 90的方位角位置测 量 能谱。当感兴趣区全能吸收峰净计数的超过 10
39、 000 时，停止测量（即计数统计涨落小 于 1%）。 2） 在 90的方位角位置，与第一次测量位置相差 120和 240的极角的情况下，各重复 测量下测量 能谱一次。 3） 比较三次方位角 90的测量值。如果有任何一个测量值大于平均值的 5%，就应计算 该探测器的径向偏移量，或者拒绝使用该探测器。若计算探测器的径向偏移量，则在后续的 的计算和测量中都应用偏移量的数值。 4） 解谱分析标准源 射线的特征峰探测 效率。 E.1.2 蒙特卡洛建模 a） 根据仪器制造商或其它方式获取的高纯锗探测器几何结构和材料信息，构建蒙特卡 洛计算模拟模型，确保高纯锗晶体、探测器外壳、冷指孔、死层和其它外壳内的组
40、件都得到 恰当的建模； 注：除厂家提供尺寸外，还可以通过 X 射线透射成像、 电子计算机断层扫描 (CT)或射线 39 透射成像等获得晶体尺寸。 b） 根据实验测量的几何条件，建立放射源、探测器、空气、支架及其它影响探测器效 率的组建模型。 c） 利用蒙特卡洛模型程序计算标准点源各特征峰的全能吸收峰探测效率。尽量增加计 算的粒子数或多次运行计算，保证各能量点的标准 偏差小于 1%。 d） 比较测量和模拟中得到的探测效率，调整探测器的相关参数。保证最终模拟计算结 果，每个能量在各位置的相对偏差都在 5%以内。 结合实验测量结果，得到最终的高纯锗探测器模型。 E.2 蒙特卡洛模拟计算有效前面积 下
41、面的步骤用于计算高纯锗探测器的有效前面积。只有经过 E.1 的步骤检验后，下面的 步骤才能够进行。 a） 建立平面源模型：源平面与探测器的端面相平行；平面源的形状和面积与探测器相 同；平面源位于探测器端面正前方； 射线是平行光束，入射至高 纯锗探测器端面。 b） 计算不同能量的 射线探测效率，能量范围涵盖低能、中能、高能区间。 c） 根据计算的探测效率与高纯锗探测器的灵敏面积相乘，计算得到高纯锗探测器的有 效前面积 。 E.3 蒙特卡洛模拟计算角响应因子 下面的步骤用于计算高纯锗探测器的角响应因子 R()。只有经过 E.1 的步骤检验后，下 面的步骤才能够进行。 a） 建立平面源模型：平面源在
42、探测器端面前方，距灵敏体积的距离大于 22rh ； 平面源所在平面的法向量与探测器端面法向量呈 度角（ 角 在 18090之间）； 射线 是平行光束，从平面源平行入射至高纯锗探测器；平面源在任意 角位置探测器晶体上的投 晶体 端面弧 外部死层 外壳 冷指孔 低温恒温器 真空 图 E.1 典型同轴 HPGe 结构示意图 40 影应大于探测器的面积。 b） 计算不同能量的 射线探测效率，能量范围涵盖低能、中能、高能区间。 c） 根据计算的探测效率与平面源的面积，得到高纯锗探测器在 角时平行 射线单位 注量率下的计数率 (Nf)。 d） 依次改变平面源所在平面的法向量与探测器端面法向量的角度，在 1
43、8090之间 角度间隔不大于 10，重新计算注量率 (Nf)。 e） 根据计算得到的注量率 (Nf)，计算不同能量在各个角度注量率与 180 角之比，得 到 (Nf)/(Nf)180。利用附录 D 中式 D4，计算得到 角响应因子 R()。 f） 对于能量为 E 的 射线，高纯锗探测器的 角度修正因子 W 用附录 D 中式 D5 计算 得到。 能量不同的的 射线 重复步骤 a) f)，得到不同能量的角度修正因子 W。 41 附录 F （资料性附录） 就地高纯锗谱仪有效前面积 及角度修正因子 W 经验值 F.1 典型就地高纯锗谱仪的有效前面积校准 对于 P 型高纯锗谱仪，有效前面积 利用下式进行
44、校准 )ln()ln( Eba (F1) 式中： 就地高纯锗谱仪的有效前面积， cm2； E射线的能量， MeV； a,b 与探测器有关的拟合参数，可用式 F2， F3 计算。 2( l n 0 . 0 9 6 9ln0 . 4 9 9 6 + 2 . 6 8 9 ）a (F2) 2 0 . 0 0 0 1 2 + 0 . 0 2 0 4 4-1 . 3 1 5 = b (F3) 式中： 高纯锗谱仪的相对探测效率，量纲一的量。 F.2 典型就地高纯锗谱仪的角度修正因子 W 高纯锗探测器晶体长度为 L cm，直径为 D cm， 可以根据 长度与直径的 L/D，参考图 F1 和图 F2 校准 P
45、型高纯锗谱仪的角度修正因子 W。 0 500 1000 1500 2000 2500 0.8 0.9 1.0 1.1 1.2 1.3 1.4 角度修正 因子 W E (k e V ) L /D= 0. 5 L /D= 0. 6 L /D= 0. 7 L /D= 0. 8 L /D= 0. 9 L /D= 1 L /D= 1. 1 L /D= 1. 2 L /D= 1. 3 图 F1 质量张弛深度 =0 gcm-2 时，就地高纯锗谱仪的角度修正因子 W 42 0 500 1000 1500 2000 2500 0.8 0.9 1.0 1.1 1.2 1.3 0.8 0.9 1.0 1.1 1.2 1.3 角度修正 因子 W E (k e V ) L /D=
