主权项 |
1.统计分布式γ或X射线能谱解谱方法:数字化谱仪或多道能谱仪向计算机输出的能谱数据组的长度为m,m为正整数,能谱数据组由道址j,j=1,2,Λ,m,射线能量E<sub>j</sub>,射线计数率C<sub>j</sub>构成,E<sub>j</sub>=a+b·j,a,b为常数,称为刻度系数,解谱步骤如下:(1)寻峰:从能谱数据组得到全能峰的峰位射线能量E<sub>k</sub>和道址j<sub>k</sub>,小于或大于j<sub>k</sub>的道址的射线计数率都小于道址j<sub>k</sub>的射线计数率C<sub>k</sub>,全能峰峰位道址为1-n个,其道址不同,k为1-n内的正整数,(2)确定左右边界:以道址j<sub>k</sub>为中心,得到小于j<sub>k</sub>的左边界道址L和大于j<sub>k</sub>的右边界道址H以及对应的边界道址射线计数率C<sub>L</sub>、C<sub>H</sub>,C<sub>L</sub>满足以下条件:C<sub>L</sub><C<sub>L-1</sub>且C<sub>L</sub><C<sub>i</sub>,i=L+1,L+2,Λ,j<sub>k</sub>C<sub>H</sub>满足以下条件:C<sub>H</sub><C<sub>H+1</sub>且C<sub>H</sub><C<sub>y</sub>,y=H-1,H-2,Λ,j<sub>k</sub>(3)确定全能峰的道址分布范围:如果C<sub>L</sub>>C<sub>H</sub>,则n=j<sub>k</sub>-L,如果C<sub>L</sub><C<sub>H</sub>,则n=H-j<sub>k</sub>,这样,全能峰的道址分布范围为:[j<sub>k</sub>-n,j<sub>k</sub>+n],(4)计算全能峰道址分布范围内各道址对应射线能量E<sub>j</sub>的分布律p<sub>j</sub>:<img file="FSA00000678099000011.GIF" wi="291" he="193" />j=j<sub>k</sub>-n,j<sub>k</sub>-n+1,Λ,j<sub>k</sub>+n,其中C<sub>j</sub>为第j道的射线计数率,(5)计算能量标准差σ<sub>k</sub>:全能峰分布范围内能量标准差,<maths num="0001"><![CDATA[<math><mrow><msub><mi>σ</mi><mi>k</mi></msub><mo>=</mo><msqrt><munderover><mi>Σ</mi><mrow><mi>j</mi><mo>=</mo><msub><mi>j</mi><mi>k</mi></msub><mo>-</mo><mi>n</mi></mrow><mrow><msub><mi>j</mi><mi>k</mi></msub><mo>+</mo><mi>n</mi></mrow></munderover><msup><mrow><mo>(</mo><msub><mi>E</mi><mi>j</mi></msub><mo>-</mo><msub><mi>E</mi><mi>k</mi></msub><mo>)</mo></mrow><mn>2</mn></msup><msub><mi>p</mi><mi>j</mi></msub></msqrt><mo>,</mo></mrow></math>]]></maths>j=j<sub>k</sub>-n,j<sub>k</sub>-n+1,Λ,j<sub>k</sub>+n,(6)建立服从N(E<sub>k</sub>,σ<sub>k</sub><sup>2</sup>)正态分布的全能峰能量概率密度函数f(E<sub>j</sub>):<maths num="0002"><![CDATA[<math><mrow><mi>f</mi><mrow><mo>(</mo><msub><mi>E</mi><mi>j</mi></msub><mo>)</mo></mrow><mo>=</mo><mfrac><mn>1</mn><mrow><msqrt><mn>2</mn><mi>π</mi></msqrt><msub><mi>σ</mi><mi>k</mi></msub></mrow></mfrac><mi>exp</mi><mrow><mo>(</mo><mfrac><msup><mrow><mo>(</mo><msub><mi>E</mi><mi>j</mi></msub><mo>-</mo><msub><mi>E</mi><mi>k</mi></msub><mo>)</mo></mrow><mn>2</mn></msup><msup><msub><mrow><mn>2</mn><mi>σ</mi></mrow><mi>k</mi></msub><mn>2</mn></msup></mfrac><mo>)</mo></mrow><mo>,</mo></mrow></math>]]></maths>j=j<sub>k</sub>-n,j<sub>k</sub>-n+1,Λ,j<sub>k</sub>+n(7)计算峰位射线能量E<sub>k</sub>处全能峰分布概率积S<sub>k</sub>,以及射线能量E<sub>j</sub>的概率比q<sub>j</sub>;<maths num="0003"><![CDATA[<math><mrow><msub><mi>S</mi><mi>k</mi></msub><mo>=</mo><mfrac><mn>1</mn><mn>0.9973</mn></mfrac><msubsup><mo>∫</mo><mrow><msub><mi>E</mi><mi>k</mi></msub><mo>-</mo><mn>3</mn><msub><mi>σ</mi><mi>k</mi></msub></mrow><mrow><msub><mi>E</mi><mi>k</mi></msub><mo>+</mo><msub><mrow><mn>3</mn><mi>σ</mi></mrow><mi>k</mi></msub></mrow></msubsup><mi>f</mi><mrow><mo>(</mo><msub><mi>E</mi><mi>j</mi></msub><mo>)</mo></mrow><msub><mi>dE</mi><mi>j</mi></msub></mrow></math>]]></maths>或<maths num="0004"><![CDATA[<math><mrow><msub><mi>S</mi><mi>k</mi></msub><mo>=</mo><mfrac><mn>1</mn><mn>0.9973</mn></mfrac><munderover><mi>Σ</mi><mrow><msub><mi>E</mi><mi>k</mi></msub><mo>-</mo><mn>3</mn><msub><mi>σ</mi><mi>k</mi></msub></mrow><mrow><msub><mi>E</mi><mi>k</mi></msub><mo>+</mo><mn>3</mn><msub><mi>σ</mi><mi>k</mi></msub></mrow></munderover><mi>f</mi><mrow><mo>(</mo><msub><mi>E</mi><mi>j</mi></msub><mo>)</mo></mrow></mrow></math>]]></maths><maths num="0005"><![CDATA[<math><mrow><msub><mi>q</mi><mi>j</mi></msub><mo>=</mo><mfrac><mrow><mi>f</mi><mrow><mo>(</mo><msub><mi>E</mi><mi>j</mi></msub><mo>)</mo></mrow></mrow><msub><mi>S</mi><mi>k</mi></msub></mfrac><mo>,</mo></mrow></math>]]></maths>j=j<sub>k</sub>-n,j<sub>k</sub>-n+1,Λ,j<sub>k</sub>+n于是峰位射线能量E<sub>k</sub>的概率比:<img file="FSA00000678099000024.GIF" wi="233" he="122" />(8)计算全能峰面积W<sub>k</sub>和峰位射线能量为E<sub>k</sub>的全能峰拟合函数F(E<sub>j</sub>):<maths num="0006"><![CDATA[<math><mrow><msub><mi>W</mi><mi>k</mi></msub><mo>=</mo><mfrac><msub><mi>C</mi><mi>k</mi></msub><msub><mi>q</mi><mi>k</mi></msub></mfrac></mrow></math>]]></maths>F(E<sub>j</sub>)=W<sub>k</sub>·q<sub>j</sub>,即为解谱函数如果只有一个全能峰,解谱结束,(9)假设有两个全能峰,对应峰位射线能量分别为E<sub>k1</sub>和E<sub>k2</sub>,且E<sub>k1</sub><E<sub>k2</sub>,分以下三种情况:a.无重叠情况:即E<sub>k1</sub><E<sub>k2</sub>-3σ<sub>k2</sub>且E<sub>k2</sub>>E<sub>k1</sub>+3σ<sub>k1</sub>W<sub>k1</sub>和W<sub>k2</sub>值均直接按照(1)-(8)的求解过程获得,b.单边重叠情况即E<sub>k1</sub>≥E<sub>k2</sub>-3σ<sub>k2</sub>且E<sub>k2</sub>>E<sub>k1</sub>+3σ<sub>k1</sub><maths num="0007"><![CDATA[<math><mrow><msub><mi>W</mi><mrow><mi>k</mi><mn>1</mn></mrow></msub><mo>=</mo><mfrac><mrow><msub><mi>C</mi><mrow><mi>k</mi><mn>1</mn></mrow></msub><mo>-</mo><msub><mi>W</mi><mrow><mi>k</mi><mn>2</mn></mrow></msub><msub><mrow><mo>(</mo><msub><mi>q</mi><mrow><mi>k</mi><mn>2</mn></mrow></msub><mo>)</mo></mrow><msub><mi>E</mi><mrow><mi>k</mi><mn>1</mn></mrow></msub></msub></mrow><msub><mi>q</mi><mrow><mi>k</mi><mn>1</mn></mrow></msub></mfrac></mrow></math>]]></maths>式中,<img file="FSA00000678099000027.GIF" wi="156" he="79" />表示E<sub>k2</sub>全能峰在射线能量E<sub>k1</sub>处的概率比,W<sub>k2</sub>的求解方式不变,或E<sub>k2</sub>≤E<sub>k1</sub>+3σ<sub>k1</sub>且E<sub>k1</sub><E<sub>k2</sub>-3σ<sub>k2</sub><maths num="0008"><![CDATA[<math><mrow><msub><mi>W</mi><mrow><mi>k</mi><mn>2</mn></mrow></msub><mo>=</mo><mfrac><mrow><msub><mi>C</mi><mrow><mi>k</mi><mn>2</mn></mrow></msub><mo>-</mo><msub><mi>W</mi><mrow><mi>k</mi><mn>1</mn></mrow></msub><msub><mrow><mo>(</mo><msub><mi>q</mi><mrow><mi>k</mi><mn>1</mn></mrow></msub><mo>)</mo></mrow><msub><mi>E</mi><mrow><mi>k</mi><mn>2</mn></mrow></msub></msub></mrow><msub><mi>q</mi><mrow><mi>k</mi><mn>2</mn></mrow></msub></mfrac></mrow></math>]]></maths>式中,<img file="FSA00000678099000029.GIF" wi="152" he="78" />表示E<sub>k1</sub>全能峰在射线能量E<sub>k2</sub>的概率比,W<sub>k1</sub>求解方式不变,c.双重叠或多重叠情况即E<sub>k1</sub>≥E<sub>k2</sub>-3σ<sub>k2</sub>且E<sub>k2</sub>≤E<sub>k1</sub>+3σ<sub>k1</sub>采用两个全能峰范围内未发生重叠的能量下相应的C<sub>j</sub>和q<sub>j</sub>值来分别求解W<sub>k1</sub>和W<sub>k2</sub>,或者利用γ或X的不同能量射线对应的谱线分数SF来开展重叠峰的分解和解谱工作。 |