CN103592287B - 通过元素的测量强度和浓度计算被干扰元素浓度的方法 - Google Patents
通过元素的测量强度和浓度计算被干扰元素浓度的方法 Download PDFInfo
- Publication number
- CN103592287B CN103592287B CN201210285349.5A CN201210285349A CN103592287B CN 103592287 B CN103592287 B CN 103592287B CN 201210285349 A CN201210285349 A CN 201210285349A CN 103592287 B CN103592287 B CN 103592287B
- Authority
- CN
- China
- Prior art keywords
- intensity
- concentration
- interference
- disturbed
- coefficient
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
Landscapes
- Investigating Or Analysing Materials By Optical Means (AREA)
- Analysing Materials By The Use Of Radiation (AREA)
- Investigating, Analyzing Materials By Fluorescence Or Luminescence (AREA)
Abstract
本发明公开了一种通过元素的测量强度和浓度计算被干扰元素浓度的方法，该方法包括以下两种校正模型：(1)谱线重叠校正模型： (2)谱线重叠‑自吸校正模型： 本发明的有益效果是，分析精度高，稳定性好。
Description
技术领域
本发明涉及光谱分析技术应用领域，特别是一种通过元素的测量强度和浓度计算被干扰元素浓度的方法。
背景技术
分析物质在外能(电能)作用下，样品中元素的电子吸收能量后发生跃迁，由低能级到高能级，处于高能态的电子处于不稳定状态，由高能态到低能态的过程中要释放能量，这种能量将以光的形式来表现，形成元素的光谱。
即：δE＝E2-E1＝hν＝hc/λ
光谱的能量的高低对应这分析样品中元素浓度的高低，通过光学部件光栅将复合光分解成单元素的光谱。由于光栅的分辨率不足以将不同元素波长接近的元素分辨开，形成光谱分析领域的核心技术问题：光线干扰。
光谱干扰指的是由于伴随物的存在引起分析结果系统误差的效应。伴随物是指除了分析元素之外的存于样品中的所有其它成分，引起干扰的伴随物质成为干扰物。在原子发射光谱分析中，干扰效应是最复杂的问题之一，可分为光谱干扰和非光谱干扰。光谱干扰具有加和性，非光谱干扰不具有加和性。
1、光谱干扰
光谱干扰即分析元素信号(分析信号)与干扰物的辐射信号(干扰信号)分辨不开所引起的干扰，干扰信号的大小取决于干扰物自身的行为，与分析元素是否同时存在无关；分为背景干扰和谱线重叠干扰。背景干扰是指由连续发射形成的带状光谱叠加在分析线上形成的干扰，其主要来源有电子一离子复合辐射、分子辐射、杂散光等。基体含量越高，背景干扰越严重。谱线重叠主要包括直接重叠和分析线变宽引起的线翼重叠。
2、非光谱干扰
非光谱干扰是指试样基体的一个或者几个伴随物使已分辨开的分析信号(谱线强度)增强或减弱的效应，这种干扰物引起的分析信号的增减，只有在分析物和干扰物同时存在时才能显示出来。在发射光谱分析中，常常是引起发射光谱准确度较差的一个原因.
元素谱线的干扰有2种类型：一是由于谱线宽度和狭缝宽度变宽引起的谱线重叠(重叠校正)；另一种是由于谱线的自吸而产生的干扰。谱线干扰里重要的是光谱干扰。
基于传统的技术背景，干扰校正采用干扰系数校正法，它是通过曲线拟合干扰元素对分析元素干扰程度的干扰系数(k系数)，进而对分析元素进行光谱干扰校正的校正技术。干扰系数校正法可校正各类干扰信号，操作简便。干扰系数校正法须事先清楚干扰物的浓度和干扰体系的组成，并要求仪器性能比较稳定。
传统校正方法我们采用的是干扰系数校正法。其校正模型有以下几种：
1、加法浓度校正模型(AC)
浓度加法模型用于解决谱线重叠引起的干扰，其模型为：
C＝AI3+BI2+CI+D+Σ(k1Ci 2+k2Ci)
其中,C被干扰元素的浓度，I为被干扰元素的强度，CiIi为干扰
元素的浓度。
2、乘法浓度校正模型
乘法浓度模型用于解决谱线自吸引起的干扰，其模型为：
C＝(AI3+BI2+CI+D)(1+Σ(κ1Ci 2+κ2Ci))
其中,C被干扰元素的浓度,I为被干扰元素的强度,CiIi为干扰元素的浓度。
3、谱线重叠校正模型
一个元素被其它的元素干扰，拟合时采取净强度，因此Ic＝Ii+
Σκj*Cij，其模型为：
C＝A[Ii+Σκj*Cij]3+B[Ii+Σκj*Cij]2+C[Ii+Σκj*Cij]+D
其中,C被干扰元素的浓度，Cij干扰元素的浓度；Kj为干扰系数，
Ii为被干扰元素的强度。
4、谱线自吸校正模型
一个元素被其它的元素干扰，拟合时采取净强度，因此Ic＝Ii (1+Σ(Kj*Iij)，其模型为：
C＝A[Ii(1+Σκj*Cij)]3+B[Ii(1+Σκj*Cij)]2+C[Ii(1+ΣκΣ*Cij)]+D
其中,C被干扰元素的浓度，Cij干扰元素的浓度，Kj为干扰系数， Ii为被干扰元素的强度。
然而上述传统模型中存在着缺点，由于模型中的强度变量本身不是净信号，本身带有系统误差，因此在实际分析时，这样在模型中的计算可能会引起进一步的误差，甚至有的时候会引起失控。
发明内容
本发明的目的是为了解决上述问题，设计了一种通过元素的测量强度和浓度计算被干扰元素浓度的方法。
实现上述目的本发明的技术方案为，一种通过元素的测量强度和浓度计算被干扰元素浓度的方法，该方法包括以下两种校正模型：
(1)谱线重叠校正模型：其中，Ii为被干扰元素的测量强度，Ic为被干扰元素没有干扰时的净强度，Kj为重叠干扰系数,Cij为第I块标准样品第j个具有叠加干扰性质的干扰元素的浓度，m为同一块标准样品中可能的干扰元素的个数；
(2)谱线重叠-自吸校正模型：其中，Ii为被干扰元素的测量强度，Ic为被干扰元素没有干扰时的净强度，Cij为第I块标准样品第j个具有自吸干扰性质的干扰元素的浓度,Kj＝1、、m为自吸干扰系数；m为同一块标准样品中可能的叠加干扰元素的个数；Cil为第i个标准样品中第l个具有叠加干扰性质的干扰元素的浓度，Kl＝1..k为叠加干扰系数。
所述谱线重叠校正模型的求解方法为：假设标准样品有n块，因此我们共有n组数据，根据加权最小二乘法，有:
设Ic＝Ci*K0,其中Ci为第i块标准样品被干扰元素的浓度， K0为净信号时的系数，Cij为第I块标准样品第j个干扰元素的浓度，Kj为重叠干扰校正系数,m为重叠干扰元素的个数，n为标准样品个数,则上式可写为
由于
得
可整理得到如下方程组
由上述式子分别求解出重叠干扰校正系数K0K1K2…Km。根据经典模型C＝A*I3+B*I2+C*I+D进行曲线拟合，求得曲线的拟合系数A、B、C、D, 然后计算被干扰元素浓度和干扰元素的计算浓度。然后将计算出的干扰元素的浓度和被干扰元素的浓度代入谱线重叠校正模型 求得每块样品被干扰元素的净信号强度Ic，再通过分析元素浓度的经典模型C＝A*Ic3+B*Ic2+C*Ic+D进行曲线拟合，求得曲线的拟合系数A、B、C、D进而求得未知样品的被干扰元素浓度C。
所述谱线重叠-自吸校正模型的求解方法为:假设标准样品有n块，因此我们共有n组数据，根据加权最小二乘法，有
设Ic＝Ci*K0其中Ci第i块标准样品被干扰元素的浓度，K0为没有干扰时的重叠-自吸校正系数，则上式可写为
根据
得
可整理得如下方程组：
通过递归算法，求得分别求解出系数K0、Kj＝1..m、Kl＝1..l。根据经典模型C＝A*I3+B*I2+C*I+D进行曲线拟合，求得曲线的拟合系数A、B、C、 D,然后计算被干扰元素浓度和干扰元素的计算浓度。然后将计算出的干扰元素的浓度和被干扰元素的浓度代入谱线重叠-自吸校正模型 然后计算每块样品的被干扰元素的净信号强度Ic，再通过经典模型C＝A*Ic3+B*Ii(Ic2+C*Ic+D进行曲线拟合，求得曲线的拟合系数A、B、C、D进而求的未知样品被干扰元素的浓度含量C。
本发明的技术方案所述的通过元素的测量强度和浓度计算被干扰元素浓度的方法，采取新的强度干扰校正模型，参与计算的参数都是净信号，这样就能大大的能够提高元素的准确性和分析精度，可使我们国产设备的对元素分析数据的稳定性、精确度、准确度达到国外设备的先进水平；而且相应设备的制作成本低，减轻了使用单位的生产负担，特别是生产铬、镍、铜、铝的生产企业。
附图说明
图1是本发明的所述通过元素的测量强度和浓度计算被干扰元素浓度的方法的分析流程图；
具体实施方式
下面结合附图对本发明进行具体描述，如图1是本发明所述通过元素的测量强度和浓度计算被干扰元素浓度的方法的结构示意图，如图所示，一种通过元素的测量强度和浓度计算被干扰元素浓度的方法，该方法包括以下两种校正模型：
(1)谱线重叠校正模型：其中，Ii为被干扰元素的测量强度，Ic为被干扰元素没有干扰时的净强度，Kj为重叠干扰系数,Cij为第I块标准样品第j个具有叠加干扰性质的干扰元素的浓度，m为同一块标准样品中可能的干扰元素的个数；
(2)谱线重叠-自吸校正模型：其中，Ii为被干扰元素的测量强度，Ic为被干扰元素没有干扰时的净强度，Cij为第I块标准样品第j个具有自吸干扰性质的干扰元素的浓度,Kj＝1、、m为自吸干扰系数；m为同一块标准样品中可能的叠加干扰元素的个数；Cil为第i个标准样品中第l个具有叠加干扰性质的干扰元素的浓度，Kl＝1..k为叠加干扰系数。
所述谱线重叠校正模型的求解方法为：假设标准样品有n块，因此我们共有n组数据，根据加权最小二乘法，有:
设Ic＝Ci*K0,其中Ci为第i块标准样品被干扰元素的浓度， K0为净信号时的系数，Cij为第I块标准样品第j个干扰元素的浓度，Kj为重叠干扰校正系数,m为重叠干扰元素的个数，n为标准样品个数,则上式可写为
由于
得
可整理得到如下方程组
由上述式子分别求解出重叠干扰校正系数K0K1K2…Km。根据经典模型C＝A*I3+B*I2+C*I+D进行曲线拟合，求得曲线的拟合系数A、B、C、D, 然后计算被干扰元素浓度和干扰元素的计算浓度。然后将计算出的干扰元素的浓度和被干扰元素的浓度代入谱线重叠校正模型 求得每块样品被干扰元素的净信号强度Ic，再通过分析元素浓度的经典模型C＝A*Ic3+B*Ic2+C*Ic+D进行曲线拟合，求得曲线的拟合系数A、B、C、D进而求得未知样品的被干扰元素浓度C。
所述谱线重叠-自吸校正模型的求解方法为:假设标准样品有n 块，因此我们共有n组数据，根据加权最小二乘法，有
设Ic＝Ci*K0其中Ci第i块标准样品被干扰元素的浓度，K0为没有干扰时的重叠-自吸校正系数，则上式可写为
根据
得
可整理得如下方程组：
通过递归算法，求得分别求解出系数K0、Kj＝1..m、Kl＝1..l。根据经典模型C＝A*I3+B*I2+C*I+D进行曲线拟合，求得曲线的拟合系数A、B、 C、D,然后计算被干扰元素浓度和干扰元素的计算浓度。然后将计算出的干扰元素的浓度和被干扰元素的浓度代入谱线重叠-自吸校正模型 然后计算每块样品的被干扰元素的净信号强度Ic，再通过经典模型C＝A*Ic3+B*Ii(Ic2+C*Ic+D进行曲线拟合，求得曲线的拟合系数A、B、C、D进而求的未知样品被干扰元素的浓度含量C。
本模型的特点是：通过计算干扰元素的浓度来得到样品中杂质元素的含量。本模型对于计算合金中的非金属元素如C、P、S有很好的效果。
上述技术方案仅体现了本发明技术方案的优选技术方案，本技术领域的技术人员对其中某些部分所可能做出的一些变动均体现了本发明的原理，属于本发明的保护范围之内。
Claims (1)
1.一种通过元素的测量强度和浓度计算被干扰元素浓度的方法，其特征在于，该方法包括以下两种校正模型：
(1)谱线重叠校正模型：其中，Ii为被干扰元素的测量强度，Ic为被干扰元素没有干扰时的净强度，kj为重叠干扰系数，Cij为第I块标准样品第j个具有叠加干扰性质的干扰元素的浓度，m为同一块标准样品中可能的干扰元素的个数；
(2)谱线重叠-自吸校正模型：其中，Ii为被干扰元素的测量强度，Ic为被干扰元素没有干扰时的净强度，Cij为第I块标准样品第j个具有自吸干扰性质的干扰元素的浓度，kj＝1..m为自吸干扰系数；m为同一块标准样品中可能的叠加干扰元素的个数；Cil为第i个标准样品中第l个具有叠加干扰性质的干扰元素的浓度，pl＝1..u为叠加干扰系数；
所述谱线重叠校正模型的求解方法为：假设标准样品有n块，因此我们共有n组数据，根据加权最小二乘法，有：
设Ic＝Ci*k0，其中Ci为第i块标准样品被干扰元素的浓度，k0为净信号时的系数，Cij为第I块标准样品第j个干扰元素的浓度，kj为重叠干扰校正系数，m为重叠干扰元素的个数，n为标准样品个数，则上式可写为：
由于
得
可整理得到如下方程组：
由上述式子分别求解出重叠干扰校正系数k0、k1、k2..km，根据经典模型Conc＝A*I3+B*I2+C*I+D进行曲线拟合，求得曲线的拟合系数A、B、C、D，然后计算被干扰元素浓度和干扰元素的计算浓度，然后
将计算出的干扰元素的浓度和被干扰元素的浓度代入谱线重叠校正模型求得每块样品被干扰元素的净信号强度Ic，再通过分析元素浓度的经典模型Conc＝A*I3+B*I2+C*I+D进行曲线拟 合，求得曲线的拟合系数A、B、C、D进而求得未知样品的被干扰元素浓度Conc；
所述谱线重叠-自吸校正模型的求解方法为：假设标准样品有n块，因此我们共有n组数据，根据加权最小二乘法，有
设Ic＝Ci*k0其中Ci是第i块标准样品被干扰元素的浓度，k0为没有干扰时的重叠-自吸校正系数，则上式可写为
根据
得
可整理得如下方程组：
通过递归算法，分别求解出系数k0、k1..m、p1..u，根据经典模型Conc＝A*I3+B*I2+C*I+D进行曲线拟合，求得曲线的拟合系数A、B、C、D，然后计算被干扰元素浓度和干扰元素的计算浓度， 然后将计算出的干扰元素的浓度和被干扰元素的浓度代入谱线重叠-自吸校正模型然后计算每块样品的被干扰元素的净信号强度Ic，再通过经典模型Conc＝A*I3+B*I2+C*I+D进行曲线拟合，求得曲线的拟合系数A、B、C、D进而求得未知样品被干扰元素的浓度含量Conc。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201210285349.5A CN103592287B (zh) | 2012-08-13 | 2012-08-13 | 通过元素的测量强度和浓度计算被干扰元素浓度的方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201210285349.5A CN103592287B (zh) | 2012-08-13 | 2012-08-13 | 通过元素的测量强度和浓度计算被干扰元素浓度的方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN103592287A CN103592287A (zh) | 2014-02-19 |
| CN103592287B true CN103592287B (zh) | 2019-09-17 |
Family
ID=50082512
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201210285349.5A Active CN103592287B (zh) | 2012-08-13 | 2012-08-13 | 通过元素的测量强度和浓度计算被干扰元素浓度的方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN103592287B (zh) |
Families Citing this family (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN104297216A (zh) * | 2014-10-27 | 2015-01-21 | 天津速伦科技有限公司 | 基于共存元素强度分析样品中相互干扰元素浓度的方法 |
| CN105004707B (zh) * | 2015-07-06 | 2017-07-14 | 浙江大学 | Ccd阵列在线拉曼光谱仪谱图标准化方法 |
| CN106383135B (zh) * | 2016-08-18 | 2018-10-23 | 广州市怡文环境科技股份有限公司 | 基于全反射x射线荧光光谱的元素定量分析方法和系统 |
Family Cites Families (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| DE102006000805B4 (de) * | 2006-01-03 | 2016-09-01 | Spectro Analytical Instruments Gmbh & Co. Kg | Verfahren zur Korrektur von spektralen Störungen in der ICP Emissionsspektroskopie (OES) |
| ITRM20090617A1 (it) * | 2009-11-25 | 2011-05-26 | Consiglio Nazionale Ricerche | Metodo ed apparato per misure di radiazione luminosa isotropica ottenuta da tecniche di spettroscopia laser, in particolare per misure di particolato submicronico. |
| JP5540952B2 (ja) * | 2010-07-09 | 2014-07-02 | ソニー株式会社 | 蛍光強度補正方法及び蛍光強度算出装置 |
| CN102135448B (zh) * | 2010-12-31 | 2012-08-22 | 天津大学 | 一种防止环境光干扰的表面反射率光谱测量方法 |
| CN102033048B (zh) * | 2010-12-31 | 2012-05-23 | 天津大学 | 一种消除环境杂散光干扰的吸收光谱测量方法 |
-
2012
- 2012-08-13 CN CN201210285349.5A patent/CN103592287B/zh active Active
Non-Patent Citations (3)
| Title |
|---|
| 《波长色散X射线荧光光谱法测定古陶瓷胎釉中37个主次痕量元素》;段鸿莺等;《矿盐测试》;20110630;第30卷(第3期);第339页 |
| 《石墨样品中杂质元素的X-射线荧光光谱法直接测定》;包生祥;《分析测试通报》;19891231;第8卷(第3期);第74页 |
| 《钢中痕量硼的光谱分析》;王娟;《一重技术》;20090228(第1期);第43页 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN103592287A (zh) | 2014-02-19 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN102262076B (zh) | 基于谱线组合的激光诱导击穿光谱元素浓度测量方法 | |
| Ruuskanen et al. | Eddy covariance VOC emission and deposition fluxes above grassland using PTR-TOF | |
| Saarikoski et al. | Chemical characterization of springtime submicrometer aerosol in Po Valley, Italy | |
| CN103954567B (zh) | 基于连续统去除法的土壤盐分测定方法 | |
| CN101819168B (zh) | 一种光谱分析方法 | |
| Zheng et al. | Estimation of organic matter content in coastal soil using reflectance spectroscopy | |
| Accomando et al. | Photon-initiated production of a dilepton final state at the LHC: Cross section versus forward-backward asymmetry studies | |
| CN103592287B (zh) | 通过元素的测量强度和浓度计算被干扰元素浓度的方法 | |
| CN103543132B (zh) | 一种基于小波变换的煤质特性测量方法 | |
| Faberi et al. | Fatty acid composition and δ13C of bulk and individual fatty acids as marker for authenticating Italian PDO/PGI extra virgin olive oils by means of isotopic ratio mass spectrometry | |
| CN105510427B (zh) | 一种用于多种元素同位素双稀释剂技术的数值处理方法 | |
| CN104198416A (zh) | 一种光谱仪波长漂移引起的测量误差实时补偿方法 | |
| CN106126879B (zh) | 一种基于稀疏表示技术的土壤近红外光谱分析预测方法 | |
| CN104076003B (zh) | 一种矿物光谱吸收特征参数的提取方法 | |
| CN103592286A (zh) | 通过元素的计算浓度来计算被干扰元素浓度的方法 | |
| JP6725154B2 (ja) | クロロフィル蛍光測定装置 | |
| CN103076387B (zh) | 一种基于dart-tof技术的书写时间鉴别方法 | |
| CN105787518B (zh) | 一种基于零空间投影的近红外光谱预处理方法 | |
| Lin-Sheng et al. | Continuous wavelet analysis for diagnosing stress characteristics of leaf powdery mildew | |
| CN103592226B (zh) | 通过元素的测量强度计算被干扰元素浓度的方法 | |
| JP6619323B2 (ja) | コメ産地推定方法 | |
| Yu et al. | A weighted ensemble method based on wavelength selection for near-infrared spectroscopic calibration | |
| CN111610073A (zh) | 一种含重金属的云南三七标准物质制备方法 | |
| Fang et al. | A inversion model for remote sensing of leaf water content based on the leaf optical property | |
| CN105865625B (zh) | 多维信息一体化获取相机的光谱定标方法 |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |