CN109459397A - 一种基于光谱高差特征的水质参数预测方法 - Google Patents
一种基于光谱高差特征的水质参数预测方法 Download PDFInfo
- Publication number
- CN109459397A CN109459397A CN201811596625.3A CN201811596625A CN109459397A CN 109459397 A CN109459397 A CN 109459397A CN 201811596625 A CN201811596625 A CN 201811596625A CN 109459397 A CN109459397 A CN 109459397A
- Authority
- CN
- China
- Prior art keywords
- area
- cod
- water
- turbidity
- water sample
- 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.)
- Granted
Links
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 title claims abstract description 109
- 238000000034 method Methods 0.000 title claims abstract description 25
- 238000001228 spectrum Methods 0.000 title claims abstract description 14
- 238000002835 absorbance Methods 0.000 claims abstract description 27
- 238000005259 measurement Methods 0.000 claims abstract description 8
- 238000001514 detection method Methods 0.000 claims abstract description 5
- 238000004611 spectroscopical analysis Methods 0.000 claims description 26
- 239000002131 composite material Substances 0.000 claims description 14
- 230000003287 optical effect Effects 0.000 claims description 4
- 238000012544 monitoring process Methods 0.000 description 7
- 238000010586 diagram Methods 0.000 description 3
- 239000003651 drinking water Substances 0.000 description 3
- 235000020188 drinking water Nutrition 0.000 description 3
- 239000000203 mixture Substances 0.000 description 2
- 239000003153 chemical reaction reagent Substances 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 238000011109 contamination Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- JEGUKCSWCFPDGT-UHFFFAOYSA-N h2o hydrate Chemical compound O.O JEGUKCSWCFPDGT-UHFFFAOYSA-N 0.000 description 1
- 239000000126 substance Substances 0.000 description 1
- 239000002352 surface water Substances 0.000 description 1
Classifications
-
- 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
Landscapes
- Physics & Mathematics (AREA)
- Spectroscopy & Molecular Physics (AREA)
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Chemical & Material Sciences (AREA)
- Analytical Chemistry (AREA)
- Biochemistry (AREA)
- General Health & Medical Sciences (AREA)
- General Physics & Mathematics (AREA)
- Immunology (AREA)
- Pathology (AREA)
- Investigating Or Analysing Materials By Optical Means (AREA)
Abstract
本发明公开了一种基于光谱高差特征的水质参数预测方法，包括使用无浊度水样数据在特征波段建立COD参数值与COD特征面积的相关性模型FC，使用普通水样数据在特征波段建立浊度参数值与浊度特征面积的相关性模型FT，同时使用梯形形状描述浊度特征，以及建立梯形自由底边与特定可见光段的吸光度高差值H相关性模型Fh，在预测时，提取待测水样在特征波段的面积特征和高差特征，通过相关性模型计算浊度和COD参数。本发明基于实际水样数据与其特征，能够最大程度避免浊度因素对于水质COD参数测算的影响，从而精确测算出水质浊度与COD参数值，为水质检测提供便利。
Description
技术领域
本发明涉及水质监测技术领域，尤其涉及一种基于光谱高差特征的水质参数预测方法。
背景技术
饮水安全问题直接关系到广大人民群众的健康，积极探索采用先进可行的监测技术与方法，实时掌握重要饮用水水源地水质变化情况，加强突发水污染事件水质预警和应急处理能力是十分必要和迫切的。
中国国家水资源监控能力建设2012-2014年项目(以下简称一期项目)已于2015年完成全部建设内容。一期项目完成后，在填补中国水资源监控手段缺乏、改善水资源管理基础设施薄弱状况、提高水资源源管理信息化永平等方面发挥了重要作用。通过一期项目建设基本实现了列入《全国重要饮用水水源地名录》(水资源函[2011]109号)的重要地表水水源地水质在线监测全覆盖。
近年来，中国质相关指标仍然使用化学法，虽然测量较为准确，但是耗时长，所需化学试剂多，操作复杂，稳定性差，二次污染严重，因此光谱法应运而生，其监测周期短，速度快，是在线监测的必备手段。
发明内容
发明目的：本发明针对现有技术存在的问题，提供一种基于光谱高差特征的水质参数预测方法，能够有效区分不同水质的浊度，测算出较为精准的水质浊度和COD参数。
技术方案：本发明所述的基于光谱高差特征的水质参数预测方法包括：
(1)根据无浊度水样的光谱数据和COD参数值，在特征波段建立COD参数值与COD特征面积的相关性模型FC；
(2)根据普通水样的光谱数据以及COD参数值和浊度参数值，在特征波段建立浊度参数值与浊度特征面积的相关性模型FT，以及被描述为梯形的浊度特征面积的自由底边与可见光波段的吸光度差值H的相关性模型Fh；
(3)对于待测水样，计算其光谱数据在特征波段的混合特征面积S'和在可见光波段的吸光度差值H'；
(4)根据可见光波段的吸光度高差值H'和模型Fh计算得到浊度特征面积，并将混合特征面积S'和浊度特征面积的差作为COD特征面积；
(5)将待测水样的浊度特征面积代入模型FT计算得到待测水样的浊度参数值，将待测水样的COD特征面积代入模型FC计算得到待测水样的COD参数值。
进一步的，步骤(1)具体包括：
(1.1)获取多条无浊度水样的光谱数据和COD参数值；
(1.2)计算每条无浊度水样在特征波段范围内的包络面积，作为COD特征面积SC；
(1.3)采用多条无浊度水样的COD参数值和COD特征面积SC数据进行拟合，得到COD参数值与COD特征面积SC的相关性模型FC。
进一步的，步骤(2)具体包括：
(2.1)获取多条普通水样的光谱数据以及其COD参数值和浊度参数值；
(2.2)计算每条普通水样在特征波段范围内的包络面积，作为包含浊度特征和COD特征的混合特征面积STC；
(2.3)将每条普通水样的COD参数值带入模型FC，从而得到对应COD特征面积SC；
(2.4)将每条普通水样的混合特征面积STC减去其COD特征面积SC，得到对应的浊度特征面积ST；
(2.5)采用多条普通水样的浊度参数值和浊度特征面积ST数据进行拟合，得到浊度参数值与浊度特征面积ST的相关性模型FT；
(2.6)计算每条普通水样的光谱数据在可见光波段的吸光度差值H；
(2.7)计算每条普通水样的浊度特征面积ST的自由底边h；其中，浊度特征面积ST为梯形；
(2.8)采用多条普通水样的吸光度差值H和自由底边h数据进行拟合，得到自由底边h与吸光度差值H的相关性模型Fh。
进一步的，步骤(3)具体包括：
(3.1)计算待测水样光谱数据在特征波段范围内的包络面积，作为混合特征面积S'；
(3.2)计算待测水样光谱数据在可见光波段的吸光度差值，记为H'。
进一步的，步骤(4)具体包括：
(4.1)将可见光波段的吸光度高差值H'代入模型Fh，得到浊度特征面积的自由底边Fh(H')；
(4.2)获取待测水样光谱数据在特征波段终点波长处的吸光度值h'；
(4.3)将h'和Fh(H')作为梯形的浊度特征面积的两个底边，按照下式计算得到浊度特征面积S'T：
式中，λ1、λ2分别为特征波段的起点波长值和终点波长值；
(4.4)计算混合特征面积S'和浊度特征面积的差S'-S'T，作为COD特征面积S'C。
可选的，所述特征波段具体为[250nm,310nm]，所述可见光波段具体为[310nm,680nm].
进一步的，所述无浊度水样是浊度测定值低于国标测定法检出限的水样，普通水样是参数测定值均在国标测定法可检测范围内的水样。
有益效果：本发明与现有技术相比，其显著优点是：本发明基于实际水样数据与其特征，能够最大程度避免浊度因素对于水质COD参数测算的影响，从而精确测算出水质浊度与COD参数值，为水质检测提供便利。
附图说明
图1是本发明的一个实施例的流程示意图；
图2为使用无浊度水样数据建模示意图；
图3为使用普通水样数据建模示意图。
具体实施方式
本实施例提供了一种基于光谱高差特征的水质参数预测方法，如图1所示，包括：
(1)根据无浊度水样的光谱数据和COD参数值，在特征波段建立COD参数值与COD特征面积的相关性模型FC。具体包括：
(1.1)获取多条无浊度水样的光谱数据和COD参数值；无浊度水样是浊度测定值低于国标测定法检出限的水样；
(1.2)计算每条无浊度水样在250nm-310nm范围内的包络面积，作为COD特征面积SC；如图2所示；
(1.3)采用多条无浊度水样的COD参数值和COD特征面积SC数据进行拟合，得到COD参数值与COD特征面积SC的相关性模型FC。
(2)根据普通水样的光谱数据以及COD参数值和浊度参数值，在特征波段建立浊度参数值与浊度特征面积的相关性模型FT，以及被描述为梯形的浊度特征面积的自由底边与可见光波段的吸光度差值H的相关性模型Fh。具体包括：
(2.1)获取多条普通水样的光谱数据以及其COD参数值和浊度参数值；所述普通水样是参数测定值均在国标测定法可检测范围内的水样；
(2.2)计算每条普通水样在250nm-310nm范围内的包络面积，作为包含浊度特征和COD特征的混合特征面积STC；
(2.3)将每条普通水样的COD参数值带入模型FC，从而得到对应COD特征面积SC；
(2.4)将每条普通水样的混合特征面积STC减去其COD特征面积SC，得到对应的浊度特征面积ST；
(2.5)采用多条普通水样的浊度参数值和浊度特征面积ST数据进行拟合，得到浊度参数值与浊度特征面积ST的相关性模型FT；
(2.6)计算每条普通水样的光谱数据从310nm至680nm的吸光度差值H，如图3所示；
(2.7)计算每条普通水样的浊度特征面积ST的自由底边h；其中，浊度特征面积ST为梯形，如图3所示，底边分别为h0和h，h为自由底边，h0为固定边；
(2.8)采用多条普通水样的吸光度差值H和自由底边h数据进行拟合，得到自由底边h与吸光度差值H的相关性模型Fh。
(3)对于待测水样，计算其光谱数据在250nm-310nm范围内的包络面积，作为混合特征面积S'，并计算光谱数据从310nm至680nm的吸光度差值，记为H'。
(4)根据可见光波段的吸光度高差值H'和模型Fh计算得到浊度特征面积，并将混合特征面积S'和浊度特征面积的差作为COD特征面积。具体包括：
(4.1)将可见光波段的吸光度高差值H'代入模型Fh，得到浊度特征面积的自由底边Fh(H')；
(4.2)获取待测水样光谱数据在310nm的吸光度值h'；
(4.3)将h'和Fh(H')作为梯形的浊度特征面积的两个底边，按照下式计算得到浊度特征面积S'T：
(4.4)计算混合特征面积S'和浊度特征面积的差S'-S'T，作为COD特征面积S'C。
(5)将待测水样的浊度特征面积S'T代入模型FT计算得到待测水样的浊度参数值TURB，将待测水样的COD特征面积S'C代入模型FC计算得到待测水样的COD参数值COD。
对公式进行整理得：
以上所揭露的仅为本发明一种较佳实施例而已，不能以此来限定本发明之权利范围，因此依本发明权利要求所作的等同变化，仍属本发明所涵盖的范围。
Claims (7)
1.一种基于光谱高差特征的水质参数预测方法，其特征在于该方法包括：
(1)根据无浊度水样的光谱数据和COD参数值，在特征波段建立COD参数值与COD特征面积的相关性模型FC；
(2)根据普通水样的光谱数据以及COD参数值和浊度参数值，在特征波段建立浊度参数值与浊度特征面积的相关性模型FT，以及被描述为梯形的浊度特征面积的自由底边与可见光波段的吸光度差值H的相关性模型Fh；
(3)对于待测水样，计算其光谱数据在特征波段的混合特征面积S'和在可见光波段的吸光度差值H'；
(4)根据可见光波段的吸光度高差值H'和模型Fh计算得到浊度特征面积，并将混合特征面积S'和浊度特征面积的差作为COD特征面积；
(5)将待测水样的浊度特征面积代入模型FT计算得到待测水样的浊度参数值，将待测水样的COD特征面积代入模型FC计算得到待测水样的COD参数值。
2.根据权利要求1所述的基于光谱高差特征的水质参数预测方法，其特征在于：步骤(1)具体包括：
(1.1)获取多条无浊度水样的光谱数据和COD参数值；
(1.2)计算每条无浊度水样在特征波段范围内的包络面积，作为COD特征面积SC；
(1.3)采用多条无浊度水样的COD参数值和COD特征面积SC数据进行拟合，得到COD参数值与COD特征面积SC的相关性模型FC。
3.根据权利要求1所述的基于光谱高差特征的水质参数预测方法，其特征在于：步骤(2)具体包括：
(2.1)获取多条普通水样的光谱数据以及其COD参数值和浊度参数值；
(2.2)计算每条普通水样在特征波段范围内的包络面积，作为包含浊度特征和COD特征的混合特征面积STC；
(2.3)将每条普通水样的COD参数值带入模型FC，从而得到对应COD特征面积SC；
(2.4)将每条普通水样的混合特征面积STC减去其COD特征面积SC，得到对应的浊度特征面积ST；
(2.5)采用多条普通水样的浊度参数值和浊度特征面积ST数据进行拟合，得到浊度参数值与浊度特征面积ST的相关性模型FT；
(2.6)计算每条普通水样的光谱数据在可见光波段的吸光度差值H；
(2.7)计算每条普通水样的浊度特征面积ST的自由底边h；其中，浊度特征面积ST为梯形；
(2.8)采用多条普通水样的吸光度差值H和自由底边h数据进行拟合，得到自由底边h与吸光度差值H的相关性模型Fh。
4.根据权利要求1所述的基于光谱高差特征的水质参数预测方法，其特征在于：步骤(3)具体包括：
(3.1)计算待测水样光谱数据在特征波段范围内的包络面积，作为混合特征面积S'；
(3.2)计算待测水样光谱数据在可见光波段的吸光度差值，记为H'。
5.根据权利要求1所述的基于光谱高差特征的水质参数预测方法，其特征在于：步骤(4)具体包括：
(4.1)将可见光波段的吸光度高差值H'代入模型Fh，得到浊度特征面积的自由底边Fh(H')；
(4.2)获取待测水样光谱数据在特征波段终点波长处的吸光度值h'；
(4.3)将h'和Fh(H')作为梯形的浊度特征面积的两个底边，按照下式计算得到浊度特征面积S'T：
式中，λ1、λ2分别为特征波段的起点波长值和终点波长值；
(4.4)计算混合特征面积S'和浊度特征面积的差S'-S'T，作为COD特征面积S'C。
6.根据权利要求1-5中任意一项所述的基于光谱高差特征的水质参数预测方法，其特征在于：所述特征波段具体为[250nm,310nm]，所述可见光波段具体为[310nm,680nm]。
7.根据权利要求1-5中任意一项所述的基于光谱高差特征的水质参数预测方法，其特征在于：所述无浊度水样是浊度测定值低于国标测定法检出限的水样，普通水样是参数测定值均在国标测定法可检测范围内的水样。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201811596625.3A CN109459397B (zh) | 2018-12-26 | 2018-12-26 | 一种基于光谱高差特征的水质参数预测方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201811596625.3A CN109459397B (zh) | 2018-12-26 | 2018-12-26 | 一种基于光谱高差特征的水质参数预测方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN109459397A true CN109459397A (zh) | 2019-03-12 |
| CN109459397B CN109459397B (zh) | 2021-02-23 |
Family
ID=65614787
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201811596625.3A Active CN109459397B (zh) | 2018-12-26 | 2018-12-26 | 一种基于光谱高差特征的水质参数预测方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN109459397B (zh) |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN113049512A (zh) * | 2021-03-12 | 2021-06-29 | 厦门斯坦道科学仪器股份有限公司 | 基于全波长紫外可见吸收光谱的水质在线监测方法 |
| CN116202975A (zh) * | 2022-12-15 | 2023-06-02 | 中国科学院西安光学精密机械研究所 | 一种水体参数预测方法、存储介质及终端设备 |
Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JPS5720648A (en) * | 1980-07-14 | 1982-02-03 | Tadashi Enatsu | Method for presuming bod and ss concentration applied with absorbance of visible light and constant turbidity cod concentration |
| CN101329251A (zh) * | 2007-06-18 | 2008-12-24 | 北京安控科技股份有限公司 | 化学需氧量、生物需氧量检测装置 |
| CN101871882A (zh) * | 2009-04-22 | 2010-10-27 | 恩德莱斯和豪瑟尔测量及调节技术分析仪表两合公司 | 用于检查液体的方法和设备 |
| CN105784617A (zh) * | 2016-04-25 | 2016-07-20 | 河海大学 | 一种基于光谱分析的水质检测方法 |
| CN107148567A (zh) * | 2014-10-29 | 2017-09-08 | 堀场仪器株式会社 | 基于吸光度和荧光确定水处理参数 |
-
2018
- 2018-12-26 CN CN201811596625.3A patent/CN109459397B/zh active Active
Patent Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JPS5720648A (en) * | 1980-07-14 | 1982-02-03 | Tadashi Enatsu | Method for presuming bod and ss concentration applied with absorbance of visible light and constant turbidity cod concentration |
| CN101329251A (zh) * | 2007-06-18 | 2008-12-24 | 北京安控科技股份有限公司 | 化学需氧量、生物需氧量检测装置 |
| CN101871882A (zh) * | 2009-04-22 | 2010-10-27 | 恩德莱斯和豪瑟尔测量及调节技术分析仪表两合公司 | 用于检查液体的方法和设备 |
| CN107148567A (zh) * | 2014-10-29 | 2017-09-08 | 堀场仪器株式会社 | 基于吸光度和荧光确定水处理参数 |
| CN105784617A (zh) * | 2016-04-25 | 2016-07-20 | 河海大学 | 一种基于光谱分析的水质检测方法 |
Non-Patent Citations (3)
| Title |
|---|
| SHOSHANA FOGELMAN MICHAEL BLUMENSTEIN等: "Estimation of chemical oxygen demand by ultraviolet spectroscopic", 《OPTIK - INTERNATIONAL JOURNAL FOR LIGHT AND ELECTRON OPTICS》 * |
| XUSONG QIN等: "Wastewater quality monitoring system using sensor fusion", 《WATER RESEARCH》 * |
| 汤斌 等: "一种紫外－可见光谱法检测水质COD的浊度影响实验研究", 《光谱学与光谱分析》 * |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN113049512A (zh) * | 2021-03-12 | 2021-06-29 | 厦门斯坦道科学仪器股份有限公司 | 基于全波长紫外可见吸收光谱的水质在线监测方法 |
| CN116202975A (zh) * | 2022-12-15 | 2023-06-02 | 中国科学院西安光学精密机械研究所 | 一种水体参数预测方法、存储介质及终端设备 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN109459397B (zh) | 2021-02-23 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN103512829B (zh) | 一种通过浑水密度测量泥沙浓度的方法 | |
| CN104655572B (zh) | 溶液显色反应定量分析检测装置 | |
| CN101852722B (zh) | 一种实现水体叶绿素a遥感反演精度评估方法 | |
| CN105486655B (zh) | 基于红外光谱智能鉴定模型的土壤有机质快速检测方法 | |
| CN103729650A (zh) | 近红外光谱建模样本的选择方法 | |
| CN104298883A (zh) | 一种油气资源评价中的烃源岩产烃率图版的建立方法 | |
| US20200271590A1 (en) | Rapid joint detection device and detection method for deterioration degree of fried oil | |
| CN112034480B (zh) | 一种双波长海洋激光雷达探测的波长择优方法 | |
| CN101216419A (zh) | 一种快速检测黄酒品质指标的方法 | |
| CN116202975A (zh) | 一种水体参数预测方法、存储介质及终端设备 | |
| CN109459397A (zh) | 一种基于光谱高差特征的水质参数预测方法 | |
| CN103592234A (zh) | 一种评价景观水体表观质量的方法 | |
| CN101477036A (zh) | 基于分段协同模型的内陆水体叶绿素a浓度遥感监测方法 | |
| CN117491301A (zh) | 一种高坝水库水环境垂向监测方法、系统及设备 | |
| CN107505283A (zh) | 近红外光谱测定次生盐渍土壤中硝酸根离子含量的方法 | |
| CN106814035A (zh) | 漂浮海表的大型藻类覆盖度估算方法 | |
| CN103293175B (zh) | 测定液体水玻璃化学成分的方法 | |
| CN105548027A (zh) | 基于近红外光谱测定调和油中茶油含量的分析模型及方法 | |
| CN105548038A (zh) | 葡萄酒颜色的测定方法 | |
| CN118961604B (zh) | 探测海洋漂浮垃圾的卫星遥感方法、装置 | |
| CN104849234A (zh) | 基于近红外光谱分析吡虫啉原药主成分含量的测定方法 | |
| WO2013076512A3 (en) | Colour strength measurement and its use in production processes | |
| CN110609017A (zh) | 一种基于反射率光谱的溢油风化时间分析方法及应用 | |
| CN105891124B (zh) | 近岸河口水色参数反演装置及方法 | |
| CN105987879B (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 | ||
| TA01 | Transfer of patent application right | ||
| TA01 | Transfer of patent application right |
Effective date of registration: 20200508 Address after: Room 712, building 1, Qilin artificial intelligence Industrial Park, 266 Chuangyan Road, Jiangning District, Nanjing City, Jiangsu Province Applicant after: NANJING BOSITU INTELLIGENT TECHNOLOGY Co.,Ltd. Address before: 211135 2, B unit 300, Zhihui Road, Qilin science and Technology Innovation Park, Nanjing, Jiangsu. Applicant before: NANJING BOSITU INTELLIGENT TECHNOLOGY Co.,Ltd. Applicant before: NANJING BOSITU ELECTRONIC SCIENCE TECHNOLOGY Co.,Ltd. |
|
| GR01 | Patent grant | ||
| GR01 | Patent grant |