氡子体测量的基本公式

1、 第三章 氡子体测量的基本公式 收集氡子体的方法很多，但用于常规测量中最可行的还是滤膜取样，即用纤维过滤器滤取一定体积的空气，收集氡子体。用贝特曼方程描述滤膜上氡子体的变化规律，邮此导出一系列测量氡子体产物的公式。本章介绍基本公式的导出过程和应用举例。 第一节 取样中的氡子体积累一、 基本假设1 符号命名 i种子体的原子浓度，原子； i种子体的放射性浓度，Bq； i种子体的衰变常数，； 采样t时刻滤膜上i种子体原子数； T衰变时间，min，采样结束时T=0； 衰变到T时刻i种子体的原子数；取样流速，L；E仪器的计数效率；R本底计数率，cpm；n样品的净计数率，cpm；N样品的积分计数；衰变率，

2、dmp；衰变率，dmp；g与取样有关的时间因子；与计数时间有关的时间因子；与计数时间间隔有关的时间因子；滤膜的过滤效率；滤膜对粒子的自吸收修正系数；滤膜对粒子的反散射系数。2 三点假定（1） 取样过程中空气中各氡子体浓度不变；（2） 取样过程中取样流速不变；（3） 滤膜对各种氡子体均有相同的过滤效率。3简化衰变链 氡及其子体对简化衰变链是推导公式的理论基础（见图1-15）。二、 基本方程式 取样过程中，由于氡是气体，滤膜无法扣留，而氡子体因是固体粒子却不断地被截留在滤膜上。与此同时，被截留在滤膜上的氡子体又在不停地衰变而减少。浓集和衰变同时进行着，这种子体的富集过程可用下面方程式来描述： （3

3、-1）式中，表示取样过程中i种子体在滤膜上的直接沉积数。表示滤膜上已沉积的前一种子体衰变产生i种子体，而使i种子体增加的数目。是i种子体自身衰变而减少的数目。三者的代数和，便是滤膜上i种子体的净增率，也就是总的变化率。因为0，故有 （3-3） 三、 方程的定解这三个方程属于一元一阶非奇次微分方程，其一般形式为 其通解为 （3-5）由初始条件定出常数C，即可得到方程的定解。在我们的具体条件下，p(x)= ,x=t,f(x)为每个微分方程式右边的系数项。 1解方程（3-2） 按照式（3-5）形式，得其通解为 由初始条件t=0, =0定常数C 0 代入上式得 (3-6)2.解方程(3-3) 把式(3

4、-6)代入式(3-3)，整理后有 用同样的积分方法，有 由初始条件t=0, =0，定常数C， 代入上式得 整理合并后，最后得 （3-7）3解方程（3-4） 把式（3-7）代入式（3-4），用同样的方法得其定解为 （3-8）四、 简化表达式 令，即可把，的表达式化简，得到描述取样过程的最终表达式，即 （3-9）第二节 取样后的衰变一、 基本方程式 取样一停止，氡子体在滤膜上的直接沉积现像已不存在，即式（3-1）中的右边第一项=0，只存在氡子体的衰变过程。于是式（3-1）变为如下形式： （3-10）同样具有下列方程组： （3-12）这三个方程同样也是一元一阶微分方程，式（3-11）为奇次方程，式（

5、3-12），（3-13）为非奇次方程。其初始条件是T=0，=。 二、 方程的定解1解式（3-11） 式（3-11）的通解为 由初始条件T=0，定常数C C=于是 = （3-14） 2解式（3-12） 把式（3-14）代入式（3-12），有 其通解为 由初始条件T=0，=，定常数C C=于是得定解 （3-15） 3解式（3-13） 把式（3-15）代入式（3-13），并运用初始条件T=0，得其定解为 （3-16）三、 简化表达式 令，即把，化简成如下形式： （3-17） 描述氡子体（也适用于Tn子体）在滤膜积累和衰变的两个方程给(3-1)，(3-10)，就是通常所说的贝特曼方程由贝特曼（Bate

6、man H.）所创立。在不同的文章中最后表达式的形式有所不同，但所描述的物理过程和最终结果是一致的。第三节 基本公式的建立取样过程中氡子体在滤膜上的积累衰变和取样停止后的氡子体衰变是两个基本过程。有了描述这两个过程的数学表达式，就可以导出测量各自氡子体浓度的基本公式。一、 完整表达式把式（3-9）中值分别代入式（3-17）中，得到 （3-18）式（3-18）就是描述两个基本过程的完整表达式，揭示了，原子数的变化规律。，是两个时间因子。在任意时间T时刻的原子数取决于：空气中自身的浓度，取样时间t，衰变时间T。在T时刻，原子数的多少不仅取决于上述三个因素，还与前一种子体的贡献有关。衰变常数是固定值

7、，不能说明变化规律，但定量地说明了数值的大小。 三、基本公式的导出放出粒子，放出粒子，而是和的混合发射体。由，贡献的总放射性活度为 （3-19）由和贡献的总放射性活度为 （3-20）原子浓度与放射性浓度C的关系是 （3-21）把关系式（3-21）分别代入式（3-19）和（3-20）中，有 （3-22） （3-23）衰变率D和计数率之间的关系为 n=DE放射性浓度C以37Bq为单位，并考虑到滤膜对粒子的自吸收，从式（3-12）导出T时刻的计数率n()为 (3-24)在实际测量中,往往是测量一段时间间隔内的积分计数，因为这样容易操作，并且会减少误差，提高方法的灵敏度。这只需要把时间因子从积分到即可

8、。 （3-25）式中 （3-26）也是时间因子。只须把式（3-24）中的时间因子换成，即可得到总积分计数的表达式。 （3-27）以=0.227261，=0.025864，=0.035185代入式(3-24)和(3-27)中，得 （3-28） （3-29）同理可得到计数率和积分计数的表达式，即 （3-30） （3-31） 式（3-28），（3-29），（3-30）（3-31）就是最后导出的四个基本公式，在氡子体产物测量中是至关重要的。，是三个未知数，为待求量。，N是可测定的量，称为待测量。，E，为已知量，只要测量装置和取样滤膜选定，这几个参数就可以定出来。变化量是时间因子g，。只与取样时间t，计

9、数时间T，计数间隔有关，而这些恰恰是人为选定的。确定了t，T，后，即可求出g，。将所有已知量代入公式后通过解联立方程即可求出三个未知量，进而求出氡子体潜能的浓度。所有测量氡子体浓度和子体潜能的方法都是从这两组基本公式推导出来的。 第四节 应用举例在实际测量中，多数是测量样品的活度，因为测量仪器和滤材的本底比较低，可得到较高的灵敏度，而且粒子的能量高，产生的脉冲幅度大，电子学测量线路较容易做到，但滤膜对粒子的自吸收修正是一个相当复杂的问题。活度的测量受本底的干扰较严重，对测量设备的要求较高，其优点是可免去自吸收修正，因为粒子的穿透能力强。现已推导出的测量方法中，多数是测粒子的，但也有测粒子的。

10、一、测量活度1推导三点法 所谓三点法就是测量取样结束后三个时刻的计数率，从而求出，三种子体的放射性浓度，。（1）确定取样时间 规定取样时间t=5min，求出的时间因子g，为=1=0.678997，=1=0.121307，=0.161319。 （2）确定计数时间 规定在取样结束后5，15，30min时刻测样品的计数率，得到的时间因子为 T515 300.3210030.0330770.0010940.8786930.6784390.4602800.8386810.5899160.348001 把，值代入式（3-28），即得到下列三个方程式： 式中：n(5)，n(15)，n(30)为取样结束后5，

11、15，30min实测的样品的计数率，变换一下上式的形式就得到了含有三个未知数的一个方程组，即 （3-32）（3）解方程组 式（3-32）属于线性方程组，求解的方法中较简单的是矩阵法。现用矩阵法解之。矩阵元为=/ ，37同样可以求出， ，37（3-33）这就是三点法测量子体浓度公式。若以为单位，则公式变为， （3-34）本公式推导过程中没有考虑本底的影响，另外，不同文章中给出的公式中系数尾数有所不同，主要是由于计算中取的有效数字不同所造成的。2推导三段法三段法就是通过测量三段时间间隔内样品的积分计数，进而求出各自氡子体浓度，。其推导过程如下：确定取样时间 规定取样时间为5min，由前面所知时间因

12、子为 0.678997 0.121307 0.161319确定计数时间间隔 选定计数时间间隔为取样后的25，620，2130min，计算时间因子。256221300.3137480.2451290.0073650.0708940.2601180.1206390.0933680.3149340.129645列联立方程 把,值代入式（3-29）中得到联立方程为 (3-35)解方程组 同样用矩阵法解此联立方程,得到,的表达式。式中的，是在取样结束后25，620，2130min三段时间间隔内测得的样品的积分计数。如果以为单位，则公式为，（3-36） 3推导马尔柯夫法 马尔柯夫法是测量氡子体潜能浓度的一

13、种快速方法，应用广泛。其导出过程如下： 确定取样时间 规定取样时间t为5min，时间因子为 0.678997 0.121307 0.161391 确定计数时间间隔 选定在取样结束后的710min内测量样品的积分计数，求出值 =710min 0.100714 0.062295 0.078306列方程式 把,值代入式(3-29),得 因为是取样结束后710min的积分计数，故写成N（7,10），由上式得 （3-37） 子体的潜能 在第一章中表1-13已列出氡及其子体的潜能值。同样是把的粒子归结到，并且只计算粒子的潜能。 由子体的放射性浓度，再利用表1-13中的数值可计算出子体的总潜能浓度，即 （3

14、-28）式中是以为单位的。 计算k值 由式(3-27)得 (3-29)把式(3-29)代入式(3-28)，有 而 其意义是每一个计数所具有的潜能值。不同平衡比下的值列于表3-1，平均值为40.3，于是有 （3-40） 表3.1 不同平衡比下的k值子体平衡比:k子体平衡比: k1 : 1 : 136.021 : 0.5 : 0.243.881 : 0.9 : 0.639.931 : 0.4 : 0.146.521 : 0.8 : 0.638.731 : 0.3 : 0.239.181 : 0.7 : 0.539.141 : 0.2 : 0.140.511 : 0.7 : 0.343.911 :

15、0.1 : 0.0539.231 : 0.6 : 0.342.301 : 0 : 036.591 : 0.5 : 0.438.00k40.3±3.06 (3-41)这两个公式中没有进行本底影响的修正,这一点将在以后的章节中进行讨论。、 二、测量活度通过测量样品的活度，同样可以导出一系列测量氡子体浓度和氡子体潜能浓度的公式，在活度的测量中，不能忽略本底的影响。1推导三段法确定取样时间 取样时间为5min，时间因子为 =0.678997, =0.121307, =0.161319确定计数时间间隔 计数时间间隔同样取采样结束后的25，620，2130min，得到的时间因子分别为256202

16、1300.3137480.2451290.0073650.0708940.2601180.1206390.0933680.3149340.129645列联立方程 把,值找入式(3-31)中,得到联立方程 (3-42)解方程 用矩阵法解方程组(3-42),得出,的表达式为 , (3-43)消除本底影响 净计数N和本底计数R的关系为 式中:G是样品加本底的总计数,将此关系代入式(3-43),得 , （3-45）2.推导马尔柯夫法计算时间因子 取样时间t=5min,计数时间间隔为710min,总的时间因子值为 =0.068385 =0.007557 =0.012632列方程式 把值代入式(3-31)中,得 (3-45)把式(3-45)代入式(3-28),有 计算k值 假定不同的子体平衡比,求得一系列k值列于表3-2,平均值为23.78,取其近似值24。此值代入上式，即得 （3-45）考虑到本底的影响，上式作如下修正： （3-47） 表3-2 子体不同