发明名称 统计分布式γ或X射线能谱解谱方法
摘要 本发明为统计分布式γ或X射线能谱解谱方法,解决已有方法通用性差、工作量大、精度差的问题。本发明针对γ或X射线全能峰中的高斯部分,建立了多个根据γ或X射线探测特点得到的具有各自物理意义的参数,并阐明了这些参数的求解过程。其技术特点是在γ或X射线探测仪器稳定或模拟条件不变的情况下,只需计算一次参数,便可长期利用这些参数进行能谱拟合,使解谱实现过程简单化、通用化,避免了传统能谱拟合函数或模拟能谱展宽函数解谱方法中参数值多变和意义缺乏的不足,同时提高了γ或X射线能谱分析的分析速度和精度。
申请公布号 CN102608649A 申请公布日期 2012.07.25
申请号 CN201210052742.X 申请日期 2012.03.02
申请人 成都理工大学 发明人 李哲;庹先国;杨剑波;刘明哲;成毅;王磊;石睿
分类号 G01T1/36(2006.01)I 主分类号 G01T1/36(2006.01)I
代理机构 成都立信专利事务所有限公司 51100 代理人 冯忠亮
主权项 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>&sigma;</mi><mi>k</mi></msub><mo>=</mo><msqrt><munderover><mi>&Sigma;</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>&pi;</mi></msqrt><msub><mi>&sigma;</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>&sigma;</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>&Integral;</mo><mrow><msub><mi>E</mi><mi>k</mi></msub><mo>-</mo><mn>3</mn><msub><mi>&sigma;</mi><mi>k</mi></msub></mrow><mrow><msub><mi>E</mi><mi>k</mi></msub><mo>+</mo><msub><mrow><mn>3</mn><mi>&sigma;</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>&Sigma;</mi><mrow><msub><mi>E</mi><mi>k</mi></msub><mo>-</mo><mn>3</mn><msub><mi>&sigma;</mi><mi>k</mi></msub></mrow><mrow><msub><mi>E</mi><mi>k</mi></msub><mo>+</mo><mn>3</mn><msub><mi>&sigma;</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来开展重叠峰的分解和解谱工作。
地址 610059 四川省成都市二仙桥东三路1号