CN113065095B - 一种基于紫外光谱的水中氮含量检测算法 - Google Patents
一种基于紫外光谱的水中氮含量检测算法 Download PDFInfo
- Publication number
- CN113065095B CN113065095B CN202110352070.3A CN202110352070A CN113065095B CN 113065095 B CN113065095 B CN 113065095B CN 202110352070 A CN202110352070 A CN 202110352070A CN 113065095 B CN113065095 B CN 113065095B
- Authority
- CN
- China
- Prior art keywords
- matrix
- nitrogen content
- concentration
- absorbance
- value
- 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
- IJGRMHOSHXDMSA-UHFFFAOYSA-N Atomic nitrogen Chemical compound N#N IJGRMHOSHXDMSA-UHFFFAOYSA-N 0.000 title claims abstract description 84
- 229910052757 nitrogen Inorganic materials 0.000 title claims abstract description 42
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 title claims abstract description 35
- 238000001514 detection method Methods 0.000 title claims abstract description 33
- 238000002211 ultraviolet spectrum Methods 0.000 title claims abstract description 13
- 238000000034 method Methods 0.000 claims abstract description 37
- 238000002835 absorbance Methods 0.000 claims abstract description 29
- 239000000126 substance Substances 0.000 claims abstract description 16
- 238000001228 spectrum Methods 0.000 claims abstract description 9
- 239000011159 matrix material Substances 0.000 claims description 48
- 238000005259 measurement Methods 0.000 claims description 12
- 238000012795 verification Methods 0.000 claims description 12
- 238000012937 correction Methods 0.000 claims description 11
- 229910002651 NO3 Inorganic materials 0.000 claims description 10
- NHNBFGGVMKEFGY-UHFFFAOYSA-N Nitrate Chemical compound [O-][N+]([O-])=O NHNBFGGVMKEFGY-UHFFFAOYSA-N 0.000 claims description 10
- MMDJDBSEMBIJBB-UHFFFAOYSA-N [O-][N+]([O-])=O.[O-][N+]([O-])=O.[O-][N+]([O-])=O.[NH6+3] Chemical compound [O-][N+]([O-])=O.[O-][N+]([O-])=O.[O-][N+]([O-])=O.[NH6+3] MMDJDBSEMBIJBB-UHFFFAOYSA-N 0.000 claims description 9
- 238000004364 calculation method Methods 0.000 claims description 7
- 238000007430 reference method Methods 0.000 claims description 5
- 230000003595 spectral effect Effects 0.000 claims description 5
- 230000002452 interceptive effect Effects 0.000 claims description 4
- 238000004458 analytical method Methods 0.000 claims description 3
- 238000000975 co-precipitation Methods 0.000 claims description 3
- 238000011156 evaluation Methods 0.000 claims description 3
- 238000005189 flocculation Methods 0.000 claims description 3
- 230000016615 flocculation Effects 0.000 claims description 3
- 238000012417 linear regression Methods 0.000 claims description 3
- 230000007935 neutral effect Effects 0.000 claims description 3
- 239000011347 resin Substances 0.000 claims description 3
- 229920005989 resin Polymers 0.000 claims description 3
- 238000001179 sorption measurement Methods 0.000 claims description 3
- 239000002245 particle Substances 0.000 claims description 2
- 238000007781 pre-processing Methods 0.000 claims 1
- 238000000870 ultraviolet spectroscopy Methods 0.000 claims 1
- 238000002798 spectrophotometry method Methods 0.000 description 5
- XTEGARKTQYYJKE-UHFFFAOYSA-M Chlorate Chemical compound [O-]Cl(=O)=O XTEGARKTQYYJKE-UHFFFAOYSA-M 0.000 description 4
- 229910052793 cadmium Inorganic materials 0.000 description 4
- BDOSMKKIYDKNTQ-UHFFFAOYSA-N cadmium atom Chemical compound [Cd] BDOSMKKIYDKNTQ-UHFFFAOYSA-N 0.000 description 4
- 238000004255 ion exchange chromatography Methods 0.000 description 4
- USHAGKDGDHPEEY-UHFFFAOYSA-L potassium persulfate Chemical compound [K+].[K+].[O-]S(=O)(=O)OOS([O-])(=O)=O USHAGKDGDHPEEY-UHFFFAOYSA-L 0.000 description 4
- 238000010521 absorption reaction Methods 0.000 description 3
- 238000005206 flow analysis Methods 0.000 description 3
- IOVCWXUNBOPUCH-UHFFFAOYSA-M Nitrite anion Chemical compound [O-]N=O IOVCWXUNBOPUCH-UHFFFAOYSA-M 0.000 description 2
- JVMRPSJZNHXORP-UHFFFAOYSA-N ON=O.ON=O.ON=O.N Chemical compound ON=O.ON=O.ON=O.N JVMRPSJZNHXORP-UHFFFAOYSA-N 0.000 description 2
- 238000002038 chemiluminescence detection Methods 0.000 description 2
- JOPOVCBBYLSVDA-UHFFFAOYSA-N chromium(6+) Chemical compound [Cr+6] JOPOVCBBYLSVDA-UHFFFAOYSA-N 0.000 description 2
- 230000029087 digestion Effects 0.000 description 2
- 238000012544 monitoring process Methods 0.000 description 2
- QJGQUHMNIGDVPM-UHFFFAOYSA-N nitrogen group Chemical group [N] QJGQUHMNIGDVPM-UHFFFAOYSA-N 0.000 description 2
- -1 nitrogen-containing compound Chemical class 0.000 description 2
- 238000004611 spectroscopical analysis Methods 0.000 description 2
- BVKZGUZCCUSVTD-UHFFFAOYSA-M Bicarbonate Chemical class OC([O-])=O BVKZGUZCCUSVTD-UHFFFAOYSA-M 0.000 description 1
- VEXZGXHMUGYJMC-UHFFFAOYSA-M Chloride anion Chemical compound [Cl-] VEXZGXHMUGYJMC-UHFFFAOYSA-M 0.000 description 1
- 238000003556 assay Methods 0.000 description 1
- 150000001649 bromium compounds Chemical class 0.000 description 1
- 150000004649 carbonic acid derivatives Chemical class 0.000 description 1
- 239000012159 carrier gas Substances 0.000 description 1
- 239000003054 catalyst Substances 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 239000003153 chemical reaction reagent Substances 0.000 description 1
- 229940079593 drug Drugs 0.000 description 1
- 239000003814 drug Substances 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 230000005281 excited state Effects 0.000 description 1
- 230000005283 ground state Effects 0.000 description 1
- 150000004767 nitrides Chemical class 0.000 description 1
- 229910017464 nitrogen compound Inorganic materials 0.000 description 1
- 150000002830 nitrogen compounds Chemical class 0.000 description 1
- 125000001477 organic nitrogen group Chemical group 0.000 description 1
- 239000007800 oxidant agent Substances 0.000 description 1
- 230000001590 oxidative effect Effects 0.000 description 1
- 230000000737 periodic effect Effects 0.000 description 1
- 239000004094 surface-active agent Substances 0.000 description 1
- 230000004083 survival effect Effects 0.000 description 1
- 238000012360 testing method Methods 0.000 description 1
- 230000007704 transition Effects 0.000 description 1
- 238000001429 visible spectrum Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/16—Matrix or vector computation, e.g. matrix-matrix or matrix-vector multiplication, matrix factorization
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N21/00—Investigating or analysing materials by the use of optical means, i.e. using sub-millimetre waves, infrared, visible or ultraviolet light
- G01N21/17—Systems in which incident light is modified in accordance with the properties of the material investigated
- G01N21/25—Colour; Spectral properties, i.e. comparison of effect of material on the light at two or more different wavelengths or wavelength bands
- G01N21/31—Investigating relative effect of material at wavelengths characteristic of specific elements or molecules, e.g. atomic absorption spectrometry
- G01N21/33—Investigating relative effect of material at wavelengths characteristic of specific elements or molecules, e.g. atomic absorption spectrometry using ultraviolet light
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/18—Complex mathematical operations for evaluating statistical data, e.g. average values, frequency distributions, probability functions, regression analysis
Abstract
一种基于紫外光谱的水中氮含量检测算法，它包括以下步骤：步骤一：对溶液进行预处理并根据溶液的吸光度进行光谱区的选择；步骤二：构建基于偏最小二乘法的从光谱到氮浓度的检测模型；步骤三：获取步骤二中所构建的检测模型中的各主因子；步骤四：采用所获得的检测模型对样本进行氮含量的检测。本发明的目的是为了精确的对水中氮物质进行氮含量的检测，且实现对水中氮含量无添加、无接触、无二次污染的检测，而提供的一种基于紫外光谱的水中氮含量检测算法。
Description
技术领域
本发明属于环境检测技术领域，具体涉及一种基于紫外光谱的水中氮含量检测算法。
背景技术
水资源作为人类赖以生存的重要资源，其品质问题始终是人们关注的重点问题。而通过行之有效的水质监测，可改善水环境污染现状并可有效保障水环境稳定与安全。随着我国不断加大生态环境保护的力度，水环境保护作为污染防治的三大攻坚战之一，备受关注，而总氮作为水环境质量的重要衡量指标，受到了越来越多的重视，且精准测定水体中总氮含量是水环境监测的基础。总氮是水体中有机氮和各种无机氮化物的总称。近年来，水中总氮的测定方法有紫外分光光度法、高温氧化-化学发光检测法、离子色谱法、连续流动分析法等。
紫外分光光度法中，碱性过硫酸钾消解紫外分光光度法是目前实验室应用最广泛的水中总氮测定方法。该方法的原理是在高温的碱性介质中，过硫酸钾可充分分解产生具有强氧化性的硫酸根自由基，该自由基能将水中的含氮化合物氧化为硝酸盐，而后用紫外分光光度法测定硝酸盐的吸光度，进而计算出总氮含量。但紫外分光光度法存在干扰物多、选择性差的问题；高温氧化-化学发光检测法是基于总氮浓度与化学发光强度呈良好的线性关系，通过检测化学发光强度，可测得样品中总氮浓度。样品由载气带入高温炉中，经高温完全气化、氧化裂解，含氮类化合物被定量转为NO，进而被O3氧化成激发态的NO2*,NO2*向基态跃迁发射光子，通过对光子浓度的检测，进而确定总氮的浓度。但该方法需定期更换催化剂；采用离子色谱法测定水中总氮，通常是先将样品中的含氮类化合物氧化成硝酸盐，而后用离子色谱分析检测出硝酸盐浓度，进而计算得到总氮浓度。离子色谱法虽然可以避免有机物的干扰，但在碱性过硫酸钾消解的过程中过硫酸钾会把体系中的氯离子氧化为氯酸盐，氯酸盐会干扰硝酸盐的测定；我国于2013年颁布了关于连续流动分析法测定水中总氮的标准(HJ667—2013)，其化学反应原理是首先利用氧化剂将氮类化合物氧化成硝酸盐，后经镉柱还原为亚硝酸盐，再采用亚硝酸盐氮的测定方法(GB7493—1987)进行测定，从而测得水中总氮含量。但连续流动分析法中，将硝酸盐还原为亚硝酸盐的镉柱需要周期性地进行效能检定，并且镉柱中镉会溶出，会对环境造成污染。
发明内容
本发明的目的是为了精确的对水中氮物质进行氮含量的检测，且实现对水中氮含量无添加、无接触、无二次污染的检测，而提供的一种基于紫外光谱的水中氮含量检测算法。
一种基于紫外光谱的水中氮含量检测算法，它包括以下步骤：
步骤一：对溶液进行预处理并根据溶液的吸光度进行光谱区的选择；
步骤二：构建基于偏最小二乘法的从光谱到氮浓度的检测模型；
步骤三：获取步骤二中所构建的检测模型中的各主因子；
步骤四：采用所获得的检测模型对样本进行氮含量的检测。
在步骤一中，通过测定硝酸根粒子在指定波长点A的吸光度测定硝酸盐氮的含量，在指定波长点B再进行吸光度的测量来校正硝酸盐氮值，所采用的硝酸盐氮含量的计算公式为：
A校＝A220-2A275 (1.1)。
在步骤一中，在进行预处理时，采用絮凝共沉淀和/或大孔中性吸附树脂，对水样中具有干扰性的有机物进行滤除。
在步骤二中，还包括分析模型评价指标步骤，所评价的指标包括校正均方根误差、预测均方根误差和RMSEP/RMSEC比值；
(1)校正均方根误差为：
yi,actual为第i个样本参考方法测定值，yi,predicted为用所建立的模型对校正集中第i个样本的预测值，n为校正集样本数；
(2)预测均方根误差为：
yi,actual为第i个样本参考方法测定值，yi,predicted为验证集验证过程中第i个样本的光谱方法预测值，m为验证集样本数；
(3)RMSEP/RMSEC比值
一般来说RMSEP越小结果越准确，但是RMSEP一般不能大于1.2倍RMSEC，不能小于0.8倍的RMSEC，也就是0.8≤RMSEP/RMSEC≤1.2，符合这种条件的模型就会具有较好的稳定性与准确性，按照概率统计，通过RMSEP可以估计预测值与化学方法得到的浓度值之间的偏差，若光谱方法得到的浓度预测值为ypre，则化学方法测得的实际值落在[ypre±RMSEP]范围的概率为67％左右，落在[ypre±2×RMSEP]范围为95％左右。
在步骤二中，所述检测模型为线性模型，其中，吸光度矩阵与浓度矩阵之间的关系为：
Y＝XB+E；(3.1)
其中吸光度矩阵为X，物质浓度矩阵为Y，B为X与Y之间的回归系数；
(1)首先对X和Y进行分解，
X＝TPT+EX； (3.2)
Y＝UQT+EY； (3.3)
T和U分别为X和Y矩阵的得分矩阵，P和Q分别为X和Y矩阵的载荷矩阵，EX和EY分别为X和Y的PLS拟合残差矩阵；
(2)对T和U做线性回归，b是回归系数矩阵，
U＝Tb； (3.4)
b＝(TTT)-1TTU； (3.5)
(3)对未知样品进行浓度预测，
Yunknow＝TunknowbQ；(3.6)
Tunknow由P和未知样本的光谱矩阵Xunknow求出。
在步骤四中，获取步骤二中所构建的检测模型中的各主因子时，采用的步骤如下：
对于校正过程，忽略残差矩阵E，主成分数取1有：
对X＝tpT，如下式所示：
t＝Xp/(pTp)；(4.1)
对Y＝uqT，如下式所示：
u＝Y/qT；(4.2)
(1)求吸光度矩阵X的权重向量w，取浓度矩阵Y的某一列作为u的起始迭代值，如下式所示用(4.4)式计算w并对其归一化：
wT＝uTX/(uTu)； (4.3)
wT＝wT/||wT||； (4.4)
(2)求吸光度矩阵X的因子得分t，由归一化后的w计算t，如下式所示：
t＝Xw/(wTw)；(4.5)
(3)求浓度矩阵Y的载荷q值，以t代替u计算q并对其归一化，如下式所示：
qT＝tTY/(tTt)； (4.6)
qT＝qT/||qT||； (4.7)
(4)求浓度矩阵Y的因子得分u，由qT计算u，如下式所示：
u＝Yq/(qTq)；(4.8)
(5)再以u代替t返回步骤(1)计算wT，由wT计算tnew，如此反复迭代，如果t已经满足收敛条件：||tnew-told||≤10-6||tnew||，转入(6)计算，否则返回(1)；
(6)由收敛后的t求吸光度矩阵X的载荷向量p并对其归一化，如下式所示：
pT＝tTY/(tTt)； (4.9)
pT＝pT/||pT||； (4.10)
(7)标准化X的得分因子t与权重向量w，
t＝t||p||； (4.11)
w＝w||p||； (4.12)
(8)计算t与u之间的内在关系b，
b＝uTt/(tTt)；(4.13)
(9)计算残差矩阵E，
EX＝X-tpT； (4.14)
EY＝Y-btqT； (4.15)
(10)以EX代替X，EY代替Y，返回(1)，以此类推，求出X、Y的各主因子(主成分)w、t、p、u、q、b，用交互验证的方法来得到最佳主因子数(最佳主成分数)f。
在步骤四中，在对样本进行氮含量的检测时，对于未知浓度的样本，已知吸光度矩阵为xun，对其物质浓度矩阵yun的预测，具体步骤如下：
(1)设置初始值h＝0，yun＝0；
(2)设h＝h+1，进行以下运算：
如果h＜f，则继续进行步骤(2)，否则停止运算，最终得到的yun就是预测值。
采用上述技术方案，能带来以下技术效果：
1、本发明提供了一种利用氮物质对紫外光谱的吸收这一原理进行氮浓度检测的方案，检测中无需添加检测试剂等药物，可避免二次污染；
2、本发明实现了无接触检测，在检测中避免了外部因素对检测精度的干扰；
3、本发明检测过程自动化、实时化，提高了检测的时效性和鲁棒性。
附图说明
下面结合附图和实施例对本发明作进一步说明：
图1是本发明中模型参数算法流程图。
具体实施方式
如图1所示一种基于紫外光谱的水中氮含量检测算法，本发明包括以下步骤：
步骤一：光谱区域选择。
紫外-可见光吸光度光谱法是根据被测物质分子对紫外-可见波段范围(100nm-800nm)光的吸收性质来进行定量和定性分析的一种方法。其中，100nm-200nm为远紫外区，200nm-400nm为紫外光谱区，400nm-800nm为可见光谱区。利用硝酸根离子在220nm波长点的吸光度测定硝酸盐氮的含量。由于溶解的有机物在220nm处也会有吸收，而硝酸根离子在275nm没有吸收，所以在275nm波长点做另一次测量，来校正硝酸盐氮值。测定下限为0.32mg/L，测定上限为4mg/L。
硝酸盐氮含量计算公式：
A校＝A220-2A275 (1.1)
溶解的有机物、表面活性剂、亚硝酸氮、六价铬、溴化物、碳酸氢盐和碳酸盐等这些物质会干扰测定，需要通过适当的预处理将这些干扰去掉。采用絮凝共沉淀和大孔中性吸附树脂，可排除水样中大部分常见有机物、浊度和Fe3+、Cr6+对测定的干扰。
步骤二：性能指标设定。
分析模型评价指标有校正均方根误差、预测均方根误差和RMSEP/RMSEC比值。
(1)校正均方根误差为：
yi,actual为第i个样本参考方法测定值，yi,predicted为用所建立的模型对校正集中第i个样本的预测值，n为校正集样本数。
(2)预测均方根误差为：
yi,actual为第i个样本参考方法测定值，yi,predicted为验证集验证过程中第i个样本的光谱方法预测值，m为验证集样本数。
(3)RMSEP/RMSEC比值
一般来说RMSEP越小结果越准确，但是RMSEP一般不能大于1.2倍RMSEC，不能小于0.8倍的RMSEC，也就是0.8≤RMSEP/RMSEC≤1.2。符合这种条件的模型就会具有较好的稳定性与准确性。按照概率统计，通过RMSEP可以估计预测值与化学方法得到的浓度值之间的偏差。若光谱方法得到的浓度预测值为ypre，则化学方法测得的实际值落在[ypre±RMSEP]范围的概率为67％左右，落在[ypre±2×RMSEP]范围为95％左右。
步骤三：构建基于偏最小二乘法(PLS)的从紫外光谱到氮浓度的检测模型。
结合实地采集的三峡库区及长江、清江、运河等水域的样本和人工检测与本项目模型检测的氮浓度之综合误差最小来获得模型参数。
一般情况下水体中所含的物质不多或者相互影响较小时浓度检测模型为线性模型。对于吸光度与浓度值呈线性的溶液可以采用PLS有效预测物质浓度。吸光度矩阵与浓度矩阵一般满足以下关系：
Y＝XB+E；(3.1)
其中吸光度矩阵为X，物质浓度矩阵为Y，B为X与Y之间的回归系数。
(3)首先对X和Y进行分解，
X＝TPT+EX； (3.2)
Y＝UQT+EY； (3.3)
T和U分别为X和Y矩阵的得分矩阵，P和Q分别为X和Y矩阵的载荷矩阵，EX和EY分别为X和Y的PLS拟合残差矩阵。
(4)对T和U做线性回归，b是回归系数矩阵，
U＝Tb； (3.4)
b＝(TTT)-1TTU； (3.5)(3)对未知样品进行浓度预测，
Yunknow＝TunknowbQ；(3.6)
Tunknow由P和未知样本的光谱矩阵Xunknow求出。
步骤四：算法实现
偏最小二乘法具体的计算由H Wold提出的非线性迭代偏最小二乘法(NIPALS)来完成。其具体的计算过程如下：
对于校正过程，忽略残差矩阵E，主成分数取1有：
对X＝tpT，如下式所示：
t＝Xp/(pTp)；(4.1)
对Y＝uqT，如下式所示：
u＝Y/qT；(4.2)
(1)求吸光度矩阵X的权重向量w。取浓度矩阵Y的某一列作为u的起始迭代值，如下式所示用(4.4)式计算w并对其归一化：
wT＝uTX/(uTu)； (4.3)
wT＝wT/||wT||； (4.4)
(2)求吸光度矩阵X的因子得分t。由归一化后的w计算t，如下式所示：
t＝Xw/(wTw)；(4.5)
(3)求浓度矩阵Y的载荷q值。以t代替u计算q并对其归一化，如下式所示：
qT＝tTY/(tTt)； (4.6)
qT＝qT/||qT||； (4.7)
(4)求浓度矩阵Y的因子得分u。由qT计算u，如下式所示：
u＝Yq/(qTq)；(4.8)
(5)再以u代替t返回步骤(1)计算wT，由wT计算tnew，如此反复迭代。如果t已经满足收敛条件：||tnew-told||≤10-6||tnew||，转入(6)计算，否则返回(1)。
(6)由收敛后的t求吸光度矩阵X的载荷向量p并对其归一化，如下式所示：
pT＝tTY/(tTt)； (4.9)
pT＝pT/||pT|| (4.10)
pT＝pT/||pT||；
(7)标准化X的得分因子t与权重向量w，
t＝t||p||； (4.11)
w＝w||p||； (4.12)
(8)计算t与u之间的内在关系b，
b＝uTt/(tTt)；(4.13)
(9)计算残差矩阵E，
EX＝X-tpT； (4.14)
EY＝Y-btqT； (4.15)
(10)以EX代替X，EY代替Y，返回(1)，以此类推，求出X、Y的各主因子(主成分)w、t、p、u、q、b。用交互验证的方法来得到最佳主因子数(最佳主成分数)f。
步骤五：对于未知浓度的样本，已知吸光度矩阵为xun，对其物质浓度矩阵yun的预测。
(1)设置初始值h＝0，yun＝0；
(2)设h＝h+1，进行以下运算：
如果h＜f，则继续进行步骤(2)，否则停止运算，最终得到的yun就是预测值。
Claims (5)
1.一种基于紫外光谱的水中氮含量检测算法，其特征在于，它包括以下步骤：
步骤一：对溶液进行预处理并根据溶液的吸光度进行光谱区的选择；
步骤二：构建基于偏最小二乘法的从光谱到氮浓度的检测模型；
步骤三：获取步骤二中所构建的检测模型中的各主因子；
步骤四：采用所获得的检测模型对样本进行氮含量的检测；
在步骤二中，所述检测模型为线性模型，其中，吸光度矩阵与浓度矩阵之间的关系为：
Y＝XB+E； (3.1)
其中吸光度矩阵为X，物质浓度矩阵为Y，B为X与Y之间的回归系数；
(1)首先对X和Y进行分解，
X＝TPT+EX； (3.2)
Y＝UQT+EY； (3.3)
T和U分别为X和Y矩阵的得分矩阵，P和Q分别为X和Y矩阵的载荷矩阵，EX和EY分别为X和Y的PLS拟合残差矩阵；
(2)对T和U做线性回归，b是回归系数矩阵，
U＝Tb； (3.4)
b＝(TTT)-1TTU； (3.5)
(3)对未知样品进行浓度预测，
Yunknow＝TunknowbQ； (3.6)
Tunknow由P和未知样本的光谱矩阵Xunknow求出；
在步骤四中，获取步骤二中所构建的检测模型中的各主因子时，采用的步骤如下：
对于校正过程，忽略残差矩阵E，主成分数取1有：
对X＝tpT，如下式所示：
t＝Xp/(pTp)； (4.1)
对Y＝uqT，如下式所示：
u＝Y/qT； (4.2)
(1)求吸光度矩阵X的权重向量w，取浓度矩阵Y的某一列作为u的起始迭代值，如下式所示用(4.4)式计算w并对其归一化：
wT＝uTX/(uTu)； (4.3)
wT＝wT/||wT||； (4.4)
(2)求吸光度矩阵X的因子得分t，由归一化后的w计算t，如下式所示：
t＝Xw/(wTw)； (4.5)
(3)求浓度矩阵Y的载荷q值，以t代替u计算q并对其归一化，如下式所示：
qT＝tTY/(tTt)； (4.6)
qT＝qT/||qT||； (4.7)
(4)求浓度矩阵Y的因子得分u，由qT计算u，如下式所示：
u＝Yq/(qTq)； (4.8)
(5)再以u代替t返回步骤(1)计算wT，由wT计算tnew，如此反复迭代，如果t已经满足收敛条件：||tnew-told||≤10-6||tnew||，转入(6)计算，否则返回(1)；
(6)由收敛后的t求吸光度矩阵X的载荷向量p并对其归一化，如下式所示：
pT＝tTY/(tTt)； (4.9)
pT＝pT/||pT||； (4.10)
(7)标准化X的得分因子t与权重向量w，
t＝t||p||； (4.11)
w＝w||p||； (4.12)
(8)计算t与u之间的内在关系b，
b＝uTt/(tTt)； (4.13)
(9)计算残差矩阵E，
EX＝X-tpT； (4.14)
EY＝Y-btqT； (4.15)
(10)以EX代替X，EY代替Y，返回(1)，以此类推，求出X、Y的各主因子w、t、p、u、q、b，用交互验证的方法来得到最佳主因子数f。
2.根据权利要求1所述的基于紫外光谱的水中氮含量检测算法，其特征在于，在步骤一中，通过测定硝酸根粒子在指定波长点A的吸光度测定硝酸盐氮的含量，在指定波长点B再进行吸光度的测量来校正硝酸盐氮值，所采用的硝酸盐氮含量的计算公式为：
A校＝A220-2A275 (1.1)。
3.根据权利要求1所述的基于紫外光谱的水中氮含量检测算法，其特征在于，在步骤一中，在进行预处理时，采用絮凝共沉淀和/或大孔中性吸附树脂，对水样中具有干扰性的有机物进行滤除。
4.根据权利要求1所述的基于紫外光谱的水中氮含量检测算法，其特征在于，在步骤二中，还包括分析模型评价指标步骤，所评价的指标包括校正均方根误差、预测均方根误差和RMSEP/RMSEC比值；
(1)校正均方根误差为：
yi,actual为第i个样本参考方法测定值，yi,predicted为用所建立的模型对校正集中第i个样本的预测值，n为校正集样本数；
(2)预测均方根误差为：
yi,actual为第i个样本参考方法测定值，yi,predicted为验证集验证过程中第i个样本的光谱方法预测值，m为验证集样本数；
(3)RMSEP/RMSEC比值
RMSEP越小结果越准确，但是RMSEP不能大于1.2倍RMSEC，不能小于0.8倍的RMSEC，也就是0.8≤RMSEP/RMSEC≤1.2，符合这种条件的模型就会具有好的稳定性与准确性，按照概率统计，通过RMSEP可以估计预测值与化学方法得到的浓度值之间的偏差。
5.根据权利要求1所述的基于紫外光谱的水中氮含量检测算法，其特征在于，在步骤四中，在对样本进行氮含量的检测时，对于未知浓度的样本，已知吸光度矩阵为xun，对其物质浓度矩阵yun的预测，具体步骤如下：
(1)设置初始值h＝0，yun＝0；
(2)设h＝h+1，进行以下运算：
如果h＜f，则继续进行步骤(2)，否则停止运算，最终得到的yun就是预测值。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202110352070.3A CN113065095B (zh) | 2021-03-31 | 2021-03-31 | 一种基于紫外光谱的水中氮含量检测算法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202110352070.3A CN113065095B (zh) | 2021-03-31 | 2021-03-31 | 一种基于紫外光谱的水中氮含量检测算法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN113065095A CN113065095A (zh) | 2021-07-02 |
| CN113065095B true CN113065095B (zh) | 2023-09-19 |
Family
ID=76565219
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202110352070.3A Active CN113065095B (zh) | 2021-03-31 | 2021-03-31 | 一种基于紫外光谱的水中氮含量检测算法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN113065095B (zh) |
Families Citing this family (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN113552085B (zh) * | 2021-07-15 | 2022-06-28 | 大连理工大学 | 一种氯离子浓度的紫外光谱检测方法 |
| CN114486786A (zh) * | 2022-03-03 | 2022-05-13 | 上海园林绿化建设有限公司 | 土壤有机质测定方法及测定系统 |
Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2000338038A (ja) * | 1999-05-28 | 2000-12-08 | Jasco Corp | スペクトルデータ処理方法 |
| CN105319179A (zh) * | 2014-06-30 | 2016-02-10 | 中国石油化工股份有限公司 | 一种利用中红外光谱预测脱硫胺液中硫化氢含量的方法 |
| CN106248617A (zh) * | 2016-07-12 | 2016-12-21 | 上海创和亿电子科技发展有限公司 | 基于近红外的烟叶焦油检测方法 |
| CN111239071A (zh) * | 2020-02-19 | 2020-06-05 | 中国科学院烟台海岸带研究所 | 一种光谱法海水硝酸盐浓度检测方法 |
| CN111351762A (zh) * | 2020-04-22 | 2020-06-30 | 青岛理工大学 | 一种紫外-可见光全波长扫描污水水质在线快速检测方法和应用 |
Family Cites Families (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| EP2430500B1 (en) * | 2009-05-14 | 2021-07-14 | Pioneer Hi-Bred International, Inc. | Inverse modeling for characteristic prediction from multi-spectral and hyper-spectral remote sensed datasets |
-
2021
- 2021-03-31 CN CN202110352070.3A patent/CN113065095B/zh active Active
Patent Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2000338038A (ja) * | 1999-05-28 | 2000-12-08 | Jasco Corp | スペクトルデータ処理方法 |
| CN105319179A (zh) * | 2014-06-30 | 2016-02-10 | 中国石油化工股份有限公司 | 一种利用中红外光谱预测脱硫胺液中硫化氢含量的方法 |
| CN106248617A (zh) * | 2016-07-12 | 2016-12-21 | 上海创和亿电子科技发展有限公司 | 基于近红外的烟叶焦油检测方法 |
| CN111239071A (zh) * | 2020-02-19 | 2020-06-05 | 中国科学院烟台海岸带研究所 | 一种光谱法海水硝酸盐浓度检测方法 |
| CN111351762A (zh) * | 2020-04-22 | 2020-06-30 | 青岛理工大学 | 一种紫外-可见光全波长扫描污水水质在线快速检测方法和应用 |
Non-Patent Citations (2)
| Title |
|---|
| 基于K-S算法的水质硝酸盐含量光谱检测方法研究;潘国锋;;光谱实验室(05);第556-560页 * |
| 紫外可见光谱检测水体COD算法研究;鲍灵利;;信息通信(02);第11-13页 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN113065095A (zh) | 2021-07-02 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN113065095B (zh) | 一种基于紫外光谱的水中氮含量检测算法 | |
| CN109709057B (zh) | 水质指标预测模型构建方法及水质指标监测方法 | |
| CN100541171C (zh) | 紫外光协同臭氧消解光度法测量水体总氮总磷的方法 | |
| CN104034684A (zh) | 一种基于紫外-可见吸收光谱的水质多指标检测方法 | |
| Florescu et al. | Validation procedure for assessing the total organic carbon in water samples | |
| CN113049512A (zh) | 基于全波长紫外可见吸收光谱的水质在线监测方法 | |
| CN203101266U (zh) | 一种水质化学需氧量快速在线监测仪器 | |
| CN106226257A (zh) | 一种水中cod在线监测装置及其监测方法 | |
| CN115810403B (zh) | 一种基于环境特征信息评估水质污染的方法 | |
| Li et al. | Global calibration model of UV-Vis spectroscopy for COD estimation in the effluent of rural sewage treatment facilities | |
| CN100543455C (zh) | 臭氧氧化紫外扫描光谱积分方式测量总有机碳的方法 | |
| CN114894725A (zh) | 一种水质多参数的光谱数据Stacking融合模型及水质多参数测量方法 | |
| Jinjun et al. | Chemiluminescence detection of permanganate index (CODMn) bya luminol-KMnO4 based reaction | |
| TWI477760B (zh) | 一種量測水中成分及其濃度之變動光徑量測裝置及其量測方法 | |
| CN101793820B (zh) | 双模联动水质在线监测方法 | |
| CN112014344B (zh) | 一种污水在线监测方法 | |
| CN113283072A (zh) | 一种适用于多场景条件的水体cod检测方法 | |
| CN111948303A (zh) | 一种利用探针化合物检测羟基自由基浓度的方法 | |
| Xie et al. | Waveband Selection with Equivalent Prediction Performance for FTIR/ATR Spectroscopic Analysis of COD in Sugar Refinery Waste Water. | |
| CN108680716B (zh) | 一种含氯废水化学需氧量初筛方法 | |
| CN1632525A (zh) | 在线光催化总磷总氮消解分析仪 | |
| JP2019113357A (ja) | 分析方法および分析装置 | |
| Liu et al. | A review on optical measurement method of chemical oxygen demand in water bodies | |
| Kong et al. | A rapid determination method of chemical oxygen demand in printing and dyeing wastewater using ultraviolet spectroscopy | |
| CN106596433A (zh) | 一种带有加热系统的水质检测系统 |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | 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 |