作者简介: 马超(1983-),男,工程师,东华理工学院本科毕业,主要从事物探方面的工作。
死时间是FD-3019型 γ测井仪的一个重要参数,笔者主要介绍死时间的定义、测量、计算和修正的方法,通过对多个钻孔测井数据死时间修正前后进行对比,并分别计算出钻孔内铀矿段含量,求出相对误差,进而说明死时间修正的必要性。
Dead time is a very important parameter of FD-3019 γ logging instrument.This paper described the definition of dead time,introduced measurement,calculation and correction methods of dead time,explained the necessity of dead time correction by comparing revised data and calculated the errors of uranium content in the drill holes.
在放射性测井中, 通常使用上海申核电子仪器有限公司生产的FD-3019型γ 测井仪来确定铀矿层的厚度、品位以及储量等。按照γ 测井规范要求, 仪器每年都需经过核工业放射性勘查计量站的检定, 在检定证书中, 有一个不能忽视的重要参数— — 死时间τ 。
由于FD-3019型γ 测井仪采用的是脉冲式γ 总量测井的方法, 其测量的能量范围为0.4~3 MeV。对于FD-3019型γ 测井仪来说, 死时间可以被定义为能区分两个脉冲之间的最小时间间隔, 其大小由探测器件的分辨率和配套的电子学电路共同决定, 一般取决于两者之中最大的一个。死时间也可以理解为γ 测井仪对辐射的响应恢复时间。
死时间又分“ 扩展型” (可延续型或累积型)和“ 非扩展型” (可延续型)两种。它们的区别在于对死时间内所到达脉冲的响应不同, 扩展型是在死时间内, 由最后脉冲到达时刻起再把死时间延长; 非扩展型是对在死时间内到达的脉冲完全没有响应。而FD-3019型γ 测井仪的死时间是非扩展型死时间, 因为在死时间内到达的脉冲没有引起计数, 但也不会进一步引起死时间的增加。
辐射具有随机性, 在死时间内, 存在着一定的辐射几率, 这个几率与辐射强度成正比, 由于γ 测井仪的死时间是非扩展型, 所以探头在死时间内将不能区分下一个脉冲的到来, 造成两个或多个脉冲被堆叠到了一起。计数率越高, 堆叠在一起的概率就越大, 对总的脉冲计数影响也就越大(图1)。
按照γ 测井仪检定规程要求, 应选择模型体源或固体镭标准源来进行测量。前者至少应选择两个含量合适的铀模型体源, 后者至少应选择两个合适的空气比释动能率值, 且要求两个模型体源的铀含量或两个空气比释动能率值应相差5~10倍。
为了评价不同铀含量或不同空气比释动能率值的死时间变化, 可采用几对模型体源或几对空气比释动能率值。要求每个点或每个模型体源的读数应不少于10次。
核工业放射性勘查计量站采用模型体源的测量方式, 优点是放射源强度稳定, 测井仪处于密闭环境, 不受其他条件影响。
放射性强度公式为
其中:X为放射性强度, 单位Bq; A为镭标准源常数, 单位(nGy· m2)/h; r为测量点到放射源的距离, 单位m。
本次实验采用的是固体镭源测量方式。先给定若干个强度值, 按照放射性强度公式(1), 计算出各给定强度所对应的距离。把γ 测井仪固定在检定架上, 将固体镭源中心与γ 测井仪测量点标线处于同一水平高度, 且距离地面高度约2 m, 然后用钢卷尺准确测量距离, 并按照各距离逐点测量。测量方式采用连续测量, 每个测点的读数都在500次以上, 取平均值参与计算, 结果见表1。
![]() | 表1 各强度下的测量距离和计数率 |
按上述方法测量, γ 测井仪的死时间τ 应按公式(2)计算
其中:τ 为γ 测井仪的死时间, 单位s; NL为在低含量模型体源(或固体镭标准源低空气比释动能率值)上的测量结果, 单位为s-1; R为低模型体源的含量(或固体镭标准源低空气比释动能率值)与高模型体源的含量(或固体镭标准源高空气比释动能率值)的比值; NH为在高含量模型体源(或固体镭标准源高空气比释动能率值)上的测量结果, 单位s-1。
按照两个空气比释动能率值应相差5~10倍的要求, 选出13组数据, 计算结果如表2所示。
![]() | 表2 实验测得的计数率及死时间 |
计算出的结果平均值为3.99 μ s, 和放射性勘查计量站测算出的3.90 μ s相对误差为2.3%, 但个别测点的误差较大, 超过了5%的要求, 其原因是由于实验所使用的固体5号镭源活度小, 在测得较大计数率时, γ 测井仪和镭源距离过近(最小距离不足半米), 稍有误差, 就会产生较大影响; 另外由于在户外进行测量, 微风等气候因素也容易造成镭源摆动, 影响结果的准确性。
γ 测井仪的死时间修正应按公式
计算, 其中:Nτ 为经过死时间修正后的测量结果, 单位s-1; N0为未经过修正的原始测量值, 单位s-1。
依据式(3), 对K548孔的伽玛测井原始曲线进行死时间修正, 绘制对比曲线(图2)。从图中可以看出, 计数率越高, 误差越大。
选取了近年来几个品位较高钻孔(砂岩型铀矿)的矿段, 按式(3)对死时间进行修正, 并计算出修正前后伽玛计数率的相对误差。再应用核工业地质局内部研发的测井资料处理解释系统中γ 测井解释模块, 按式(4)计算出该计数率峰值所处矿段的平米铀含量。由于砂岩型铀矿的工业指标是以平米铀含量作为参考, 故文中以平米铀含量进行对比, 并计算出修正前后的平米铀含量相对误差(表3)。计算公式为
其中:Y为平米铀含量, 单位kg/m2; H为由边界品位确定的矿层厚度, m; C为通过五点式反褶积分层解释法求得的品位, 单位eU0.01%; ρ 为岩石密度, 单位g/cm3。
![]() | 表3 死时间修正前后误差统计 |
从表3中可以看出, 当计数率达到50 000时, 修正前后误差超过20%, 单层矿段的平米铀含量误差近10%, 也就是说矿层的品位越高, 影响就越大。
综上所述, γ 测井数据的死时间修正非常重要, 尤其是在高品位铀矿储量的计算上, 品位越高, 影响越大, 如果用在整个矿床的储量统计上, 影响是相当大的。随着找矿技术和方法的不断进步, 或将发现更大规模储量的矿床, 因此在γ 测井数据的处理解释上, 必须进行死时间修正。
The authors have declared that no competing interests exist.
[1] |
|
[2] |
|
[3] |
|
[4] |
|
[5] |
|
[6] |
|
[7] |
|
[8] |
|
[9] |
|
[10] |
|
[11] |
|
[12] |
|
[13] |
|
[14] |
|
[15] |
|