CN107665288A - 一种化学需氧量的水质软测量预测方法 - Google Patents
一种化学需氧量的水质软测量预测方法 Download PDFInfo
- Publication number
- CN107665288A CN107665288A CN201610619777.5A CN201610619777A CN107665288A CN 107665288 A CN107665288 A CN 107665288A CN 201610619777 A CN201610619777 A CN 201610619777A CN 107665288 A CN107665288 A CN 107665288A
- Authority
- CN
- China
- Prior art keywords
- cod
- water quality
- model
- quality monitoring
- soft
- 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.)
- Pending
Links
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 title claims abstract description 128
- 238000005259 measurement Methods 0.000 title claims abstract description 49
- 238000000034 method Methods 0.000 title claims abstract description 31
- 238000012544 monitoring process Methods 0.000 claims abstract description 65
- 238000012360 testing method Methods 0.000 claims abstract description 50
- 238000012549 training Methods 0.000 claims abstract description 17
- 238000013528 artificial neural network Methods 0.000 claims abstract description 11
- 238000004422 calculation algorithm Methods 0.000 claims abstract description 8
- IJGRMHOSHXDMSA-UHFFFAOYSA-N Atomic nitrogen Chemical compound N#N IJGRMHOSHXDMSA-UHFFFAOYSA-N 0.000 claims description 14
- XKMRRTOUMJRJIA-UHFFFAOYSA-N ammonia nh3 Chemical compound N.N XKMRRTOUMJRJIA-UHFFFAOYSA-N 0.000 claims description 10
- 238000010219 correlation analysis Methods 0.000 claims description 10
- OAICVXFJPJFONN-UHFFFAOYSA-N Phosphorus Chemical compound [P] OAICVXFJPJFONN-UHFFFAOYSA-N 0.000 claims description 8
- 229910052698 phosphorus Inorganic materials 0.000 claims description 8
- 239000011574 phosphorus Substances 0.000 claims description 8
- 238000011478 gradient descent method Methods 0.000 claims description 7
- 229910052757 nitrogen Inorganic materials 0.000 claims description 6
- 238000004458 analytical method Methods 0.000 claims description 5
- 238000005457 optimization Methods 0.000 claims description 5
- 230000007613 environmental effect Effects 0.000 abstract description 2
- 239000010865 sewage Substances 0.000 description 13
- 238000005516 engineering process Methods 0.000 description 7
- 230000008569 process Effects 0.000 description 7
- 230000000694 effects Effects 0.000 description 6
- 239000003153 chemical reaction reagent Substances 0.000 description 5
- 230000006872 improvement Effects 0.000 description 5
- 230000006870 function Effects 0.000 description 4
- 230000008901 benefit Effects 0.000 description 3
- 238000004364 calculation method Methods 0.000 description 3
- 238000001514 detection method Methods 0.000 description 3
- 238000012423 maintenance Methods 0.000 description 3
- 239000000047 product Substances 0.000 description 3
- 241000208340 Araliaceae Species 0.000 description 2
- 235000005035 Panax pseudoginseng ssp. pseudoginseng Nutrition 0.000 description 2
- 235000003140 Panax quinquefolius Nutrition 0.000 description 2
- 238000005273 aeration Methods 0.000 description 2
- 230000008859 change Effects 0.000 description 2
- 238000013461 design Methods 0.000 description 2
- 238000011161 development Methods 0.000 description 2
- 235000008434 ginseng Nutrition 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 230000007935 neutral effect Effects 0.000 description 2
- 238000012545 processing Methods 0.000 description 2
- 239000007787 solid Substances 0.000 description 2
- 238000004065 wastewater treatment Methods 0.000 description 2
- 238000003911 water pollution Methods 0.000 description 2
- ZAMOUSCENKQFHK-UHFFFAOYSA-N Chlorine atom Chemical compound [Cl] ZAMOUSCENKQFHK-UHFFFAOYSA-N 0.000 description 1
- VYZAMTAEIAYCRO-UHFFFAOYSA-N Chromium Chemical compound [Cr] VYZAMTAEIAYCRO-UHFFFAOYSA-N 0.000 description 1
- QVGXLLKOCUKJST-UHFFFAOYSA-N atomic oxygen Chemical compound [O] QVGXLLKOCUKJST-UHFFFAOYSA-N 0.000 description 1
- 239000003795 chemical substances by application Substances 0.000 description 1
- 239000000460 chlorine Substances 0.000 description 1
- 229910052801 chlorine Inorganic materials 0.000 description 1
- 229910052804 chromium Inorganic materials 0.000 description 1
- 239000011651 chromium Substances 0.000 description 1
- 238000011109 contamination Methods 0.000 description 1
- 238000012937 correction Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 201000010099 disease Diseases 0.000 description 1
- 208000037265 diseases, disorders, signs and symptoms Diseases 0.000 description 1
- 235000013399 edible fruits Nutrition 0.000 description 1
- 238000012851 eutrophication Methods 0.000 description 1
- 238000013507 mapping Methods 0.000 description 1
- QSHDDOUJBYECFT-UHFFFAOYSA-N mercury Chemical compound [Hg] QSHDDOUJBYECFT-UHFFFAOYSA-N 0.000 description 1
- 229910052753 mercury Inorganic materials 0.000 description 1
- 239000000203 mixture Substances 0.000 description 1
- 210000005036 nerve Anatomy 0.000 description 1
- 230000003647 oxidation Effects 0.000 description 1
- 238000007254 oxidation reaction Methods 0.000 description 1
- 230000033116 oxidation-reduction process Effects 0.000 description 1
- 239000001301 oxygen Substances 0.000 description 1
- 229910052760 oxygen Inorganic materials 0.000 description 1
- 238000005375 photometry Methods 0.000 description 1
- 238000002360 preparation method Methods 0.000 description 1
- 238000000746 purification Methods 0.000 description 1
- 230000008439 repair process Effects 0.000 description 1
- 238000012827 research and development Methods 0.000 description 1
- 238000005070 sampling Methods 0.000 description 1
- 230000001568 sexual effect Effects 0.000 description 1
- 238000004088 simulation Methods 0.000 description 1
- 239000010802 sludge Substances 0.000 description 1
- 238000005728 strengthening Methods 0.000 description 1
- 239000000126 substance Substances 0.000 description 1
- 239000013589 supplement Substances 0.000 description 1
- 239000000725 suspension Substances 0.000 description 1
- 238000010998 test method Methods 0.000 description 1
- 238000012956 testing procedure Methods 0.000 description 1
- 238000004448 titration Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N33/00—Investigating or analysing materials by specific methods not covered by groups G01N1/00 - G01N31/00
- G01N33/18—Water
- G01N33/1806—Biological oxygen demand [BOD] or chemical oxygen demand [COD]
-
- G—PHYSICS
- G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
- G16C—COMPUTATIONAL CHEMISTRY; CHEMOINFORMATICS; COMPUTATIONAL MATERIALS SCIENCE
- G16C20/00—Chemoinformatics, i.e. ICT specially adapted for the handling of physicochemical or structural data of chemical particles, elements, compounds or mixtures
- G16C20/70—Machine learning, data mining or chemometrics
Landscapes
- Life Sciences & Earth Sciences (AREA)
- Engineering & Computer Science (AREA)
- Health & Medical Sciences (AREA)
- Chemical & Material Sciences (AREA)
- General Health & Medical Sciences (AREA)
- General Physics & Mathematics (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Emergency Medicine (AREA)
- Food Science & Technology (AREA)
- Medicinal Chemistry (AREA)
- Physics & Mathematics (AREA)
- Analytical Chemistry (AREA)
- Biochemistry (AREA)
- Biomedical Technology (AREA)
- Biodiversity & Conservation Biology (AREA)
- Immunology (AREA)
- Pathology (AREA)
- Artificial Intelligence (AREA)
- Molecular Biology (AREA)
- Data Mining & Analysis (AREA)
- Databases & Information Systems (AREA)
- Evolutionary Computation (AREA)
- Medical Informatics (AREA)
- Software Systems (AREA)
- Crystallography & Structural Chemistry (AREA)
- Bioinformatics & Cheminformatics (AREA)
- Bioinformatics & Computational Biology (AREA)
- Computing Systems (AREA)
- Theoretical Computer Science (AREA)
- Sampling And Sample Adjustment (AREA)
Abstract
本发明实施方式涉及环境监测领域,公开了一种化学需氧量的水质软测量预测方法。包括:获取待测水环境中多个水质监测指标的数据值,每个水质监测指标被获取到的数据值有多个,多个水质监测指标中至少有一个为COD;从除COD之外的各水质监测指标中选取若干个水质监测指标,将所选的水质监测指标的数据值形成数据样本集,将数据样本集划分为训练集及测试集,利用模糊神经网络算法对训练集进行训练,获得COD软测量模型;利用测试集,对COD软测量模型进行测试,获得测试结果;重复执行模型的建立和模型测试,直至获得的测试结果符合预设条件,将符合预设条件的COD软测量模型作为最终的软测量模型。所得模型预测效果精确,降低COD的水质软测量预测成本。
Description
技术领域
本发明涉及环境监测领域,特别涉及化学需氧量的水质软测量预测技术。
背景技术
众所周知,我国水资源十分稀缺,人均占有量仅占世界平均水平的1/4,同时严重的水环境的污染问题,如水污染、水体富营养化、城市黑臭水体、地下水污染等。使水资源问题愈加突出。在农村,我国将近96%的、约80~90亿t/a的农村污水没有治理,严重农村环境污染和河海流域污染,在城市,黑臭水体已然成为一种城市病,所有城市几乎无一幸免。国内大部分河流均受到不同程度的污染。缓解水资源紧张、加大污水净化、减轻水环境污染刻不容缓。
随着水环境保护及污染治理的紧迫,水质监测提供的水质信息尤为重要,水质监测是水环境保护及污染治理的重要基础,其重要性关乎掌握水资源质量状况、以及对用水、污水处理、排污等有效监控。一方面,受污染水体的水质情况都需要通过化学需氧量(COD)、氨氮、总氮(TN)、总磷(TP)等关键水质指标的监测提供准确及时的水质动态信息,另一方面,及时有效的污水处理参数监测,对污水处理系统有重大意义,污水的排放必须符合国家污水排放标准中的相关规定,这就要求我们必须检测处理出水中TN、TP、COD、氨氮等关键参数。因此加强水环境保护及污染治理需要水质监测工作超前发展。
COD、氨氮、TN、TP等关键水质指标的监测,对于水质处理的控制、过程的优化及诊断起着重要的作用,然而这类水质指标难以测量或不易在线测量,目前主要有人工实验室检测及检测水质的自动分析仪表(水质在线检测仪)。
水质在线自动分析仪,如COD在线监测仪、TN在线检测仪等,在我国发展比较晚,国内产品的测量精确度和可靠性与国外存在一定的差距,存在品种单一、测量精准度低,测量周期长等缺陷,没有完全达到满足测量需求,污水生物处理过程的水质参数无法准确地进行测量。但使用国外的产品存在检修维护困难、价格昂贵等问题,限制了其在水环境保护及污染治理领域的应用。
以COD在线监测仪举例,由于水样中部分有机物很难被氧化剂氧化,有的甚至根本不能氧化。因此COD在线监测仪难以应用于高氯污水、强碱污水、浓度大幅变动污水及地表水的自动监测,仅能满足部分污染源在线自动监测的需要。另外,采用消解-氧化还原滴定法、消解-光度法的COD在线监测仪的分析周期长,需要60分钟左右。同时COD在线监测仪价格昂贵,一次性投入成本高,日常检测需要定期消耗试剂,试剂费用大;还可能存在铬、汞的二次污染问题;试剂的更换以及泵管的更换维护方面非常烦琐,试剂费用和维护工作量都很大。
目前国内污水处理厂的出水水质关键参数如化学需氧量(COD)、氨氮、总氮(TN)、总磷(TP)等水指标大多依靠实验室人工化验检测来确知,显然,实验室人工化验存在滞后性,从采样到输出测定结果,有时间差值,时间滞后;精确度受人为影响,较难判断;无法实时监测,不能及时依据水质变化指导,难以控制污水处理情况。
目前COD监测解决思路有以下两种:
(1)改进直接测量仪表:按照传统的检测技术发展思路,以硬件形式研制新型的过程测量仪表。但由于污水处理系统是一个复杂的生物化学处理系统,是典型的非线性、多变量、不稳定、大时滞系统,干扰因素众多且不确定,直接研发、改进在线监测仪表,难度大,且监测成本(一次性投入及试剂费用)较高,难以降低。
(2)间接测量:采用间接测量的思路,利用容易获取并且与被测变量相关的测量信息,通过计算估计被测变量。软测量即是这一思想的集中体现,将软测量应用于水环境保护及污染治理中,可实现投资少、实时监测的目的。
发明内容
本发明实施方式的目的在于提供一种化学需氧量的水质软测量预测方法,在保证COD的水质软测量预测准确度的同时,降低COD的水质软测量预测成本。
为解决上述技术问题,本发明的实施方式提供了一种化学需氧量的水质软测量预测方法,包括:
数据的获取及分析:获取待测水环境中多个水质监测指标的数据值,每个水处理监测指标被获取到的数据值有多个,多个所述水质监测指标中至少有一个为化学需氧量COD;
模型的建立:从除所述COD之外的各水质监测指标中选取若干个水质监测指标,将所选的水质监测指标的数据值形成数据样本集,将所述数据样本集划分为训练集及测试集,利用模糊神经网络算法对所述训练集进行训练,获得COD软测量模型;
模型测试:利用所述测试集,通过多次仿真实验的测试,对所述COD软测量模型进行测试,获得测试结果;
重复执行所述模型的建立和所述模型测试,直至获得的测试结果符合预设条件,将符合所述预设条件的COD软测量模型作为最终的软测量模型。
本发明实施方式相对于现有技术而言,主要区别及其效果在于:在COD的软测量预测中,运用多次建模并测试的方法,从众多水质监测指标中,选取到与待测水环境更为合适的指标,使得最终获得的模型可以更为准确稳定地预测出待测水环境中COD的值。解决了现有技术中,利用同一组水处理监测指标建模,所获得的模型不能准确符合待测水环境的COD预测,拓展本发明实施方式的应用场景。另外,进一步限定获取到的数据被分为训练集和测试集,两者的数据值可以为不同的数据值,使得测试时,可以利用集外测试,使COD软测量模型的效果说服力大大提高。
作为进一步改进,所述建模步骤之前,还包括:相关性分析:利用各个水质监测指标的数据值,进行相关性分析,获得除所述COD之外的各个水质监测指标和所述COD之间的相关性;相应的,所述选取若干个水质监测指标中,利用所述相关性的高低,选取若干个水质监测指标。
利用相关性分析,使得选取指标时可以根据分析结果选取,使得选取具有了较佳的依据,提高获得符合要求的软测量模型的效率。
作为进一步改进,所述模型测试后,和所述模型测试前,还包括:通过梯度下降法对所获得的COD软测量模型进行优化;所述对COD软测量模型进行测试中,对优化后的COD软测量模型进行测试。
利用梯度下降法进行模型优化,使得模型的误差得到大幅度降低,大大提高模型精确度,通过上述模糊神经网络算法和梯度下降法,使获得的最终模型更为精确稳定。
作为进一步改进,所述根据各测试结果选定最终的软测量模型中,所述模型测试中获得的所述测试结果包括:所述COD软测量模型的时间复杂度、准确率和稳定性。进一步限定测试结果的内容,提高选定的最终的软测量模型的有效性。
作为进一步改进,所述获取步骤中,如果获取到的数据值属于:总氮、总磷或氨氮,则舍弃。由于总氮、总磷或氨氮这三个指标的直接测试成本过高,不采用这三个指标作为输入,进一步降低COD的水质软测量预测成本。
附图说明
图1是根据本发明第一实施方式中的COD的软测量预测方法流程图;
图2是根据本发明第一实施方式中的实际出水COD值与软测量所得COD值的对比及相对误差的示意图;
图3是根据本发明第二实施方式中的COD的软测量预测方法流程图。
具体实施方式
为使本发明实施方式的目的、技术方案和优点更加清楚,下面将结合附图对本发明的各实施方式进行详细的阐述。然而,本领域的普通技术人员可以理解,在本发明各实施方式中,为了使读者更好地理解本申请而提出了许多技术细节。但是,即使没有这些技术细节和基于以下各实施方式的种种变化和修改,也可以实现本申请所要求保护的技术方案。
本发明的第一实施方式涉及一种COD的水质软测量预测方法。其流程如图1所示,具体如下:
步骤101,获取待测水环境中多个水质监测指标的数据值。
具体的说,本实施方式中的水质监测指标可以是以下指标:COD、进水流量、出水流量、进水PH、出水PH、出水SS(即固体悬浮物)、曝气池MLSS(即混合液悬浮固体浓度)、回流MLSS、MLVSS(即混合液挥发性悬浮固体浓度)、曝气SVI(即污泥体积指数)、回流SVI、1#池DO(即溶解氧)、2#池DO、ORP(即氧化还原电位)、温度、HRT(即水力停留时间)。当然,在实际应用中,还可以选取其他水质监测指标,不限于上述指标。更具体的说,获取的水质监测指标中至少一个为化学需氧量COD。
更具体的说,上述水质监测指标将划分为两类,分别是输入指标及输出指标,输出指标指水质监测指标COD,COD的数据值作为输出数据样本;输入指标包括除COD外的其他水质监测指标,这些指标的数据值作为输入数据样本。
需要说明的是,每个水质监测指标被获取到的数据值有多个。一般会选择近阶段的历史数据作为各指标的数据值。
值得一提的是,本步骤中还可以进一步筛选数据,比如:如果获取到的数据值属于:总氮、总磷或氨氮,则舍弃。这样就可以不采用总氮、总磷或氨氮这三个指标,由于总氮、总磷或氨氮这三个指标的直接测试成本过高,不采用这三个指标作为输入,进一步降低COD的水质软测量预测成本。
步骤102,利用各个水质监测指标的数据值,进行相关性分析,获得除COD之外的各个水质监测指标和COD之间的相关性。
具体的说,相关性分析公式为下式(1):
其中,相关系数r表示两个变量x和y之间的关系密切程度。
值得一提的是,在实际应用中,也可以不进行相关分析,直接利用经验选取,但本实施方式还是增加了相关性分析的步骤,可以使得选取指标时可以根据分析结果选取,使得选取具有了较佳的依据,加快获得符合要求的软测量模型的速度。
步骤103,从除COD之外的各水质监测指标中选取若干个水质监测指标组合,形成数据样本集。
具体的说,本实施方式中的选取依据可以是步骤102的相关性分析结果,如选取的是相关性超过30%的若干个指标,并将所选的水质监测指标的数据值形成数据样本集。在实际应用中,还可以根据经验,设定一些相关性较大的必选参数,和没有相关性的不选参数。
具体的说,相关性结果以相关系数r表示,r的值介于–1与+1之间,即–1≤r≤+1。其性质如下:
当r>0时,表示两变量正相关,r<0时,两变量为负相关;
当|r|=1时,表示两变量为完全线性相关,即为函数关系;
当r=0时,表示两变量间无线性相关关系;
当0<|r|<1时,表示两变量存在一定程度的线性相关。且|r|越接近1,两变量间线性关系越密切;|r|越接近于0,表示两变量的线性相关越弱。一般可按三级划分:|r|<0.4为低度线性相关;0.4≤|r|<0.7为显著性相关;0.7≤|r|<1为高度线性相关。
更具体的说,选取的指标数量也不做限定,可以选取任意数量的指标。
步骤104,将数据样本集划分为两类:一类作为训练集,一类作为测试集。
具体的说,本实施方式划分后的训练集中包含的数据值可以多于测试集中包含的数据值。其中,划分时可以采用各种原则。
值得一提的是,后续建模步骤中使用的数据样本集为训练集;测试步骤中所使用的数据样本集为测试集。进一步限定训练集和测试集的数据值为不同的数据值,使得测试时,可以利用集外测试,使COD软测量模型的效果说服力大大提高。
步骤105,利用模糊神经网络算法对训练集进行训练,获得COD的软测量模型。
具体的说,尽管模糊逻辑系统和神经网络系统在结构上截然不同,但是在两者之间可以实现优势互补。神经网络具有非线性映射能力、学习能力、并行处理能力和容错能力。模糊逻辑具有处理不确定性的能力、该发明将两者的优势有机结合起来,组成在功能上更加完善和强大的系统应用到污水水质监测中。即采用模糊神经网络(FNN)建立COD的软测量模型。
更具体的说,本实施方式中建立的模糊神经网络结构分为4层,分别为输入层、模糊化层、模糊规则层、输出层。各层之间依据模糊逻辑系统的语言变量、模糊if-then规则、模糊推理方法、反模糊函数相联系。在规则为Ri的情况下,模糊推理如下:
其中,为模糊系统的模糊集;为模糊系统参数;yi为根据模糊规则得到的输出,输入部分(即if部分)是模糊的,输出部分(即then部分)是确定的,该模糊推理表示输出为输入的线性组合。
假设对于输入量x=[x1,x2,…,xk],首先根据模糊规则计算各输入变量xj的隶属度。
上式(2)中,分别为隶属度函数的中心和宽度;k为输入参数个数;n为模糊子集数。
下式(3)中,将各隶属度进行模糊计算,采用模糊算子为连乘算子。
下式(4)中,根据模糊计算结果计算模糊模型的输出值yi。
模糊神经网络分为输入层、模糊化层、模糊规则计算层和输出层等四层。输入层与输入向量xi连接,节点数与输入的维数相同。模糊化层采用隶属度函数(1)对输入值进行模糊化得到模糊隶属度值u。模糊规则计算层采用模糊连乘(2)计算得到w。输出层采用个公式(3)计算模糊神经网络的输出。
步骤106,利用梯度下降法对获得COD软测量模型进行优化,获得优化后的COD水质软测量模型。
具体的说,梯度下降法优化计算包括
(1)误差计算
式(5)中,yd是网络期望输出;yc是网络实际输出;e为期望输出和实际输出的误差。
(2)系数修正
式(6)和(7)中,为神经网络系数;α为网络学习率;xj为网络输入参数;wi为输入参数隶属度连乘积。
(3)参数修正
式(8)和(9)中,分别为隶属度函数的中心和宽度。
步骤107,利用测试集对所获得优化后的COD的软测量模型进行测试,获得测试结果。
具体的说,本实施方式中测试结果包括:软测量模型的时间复杂度、准确率和稳定性。为了获得更佳的预测模型,所以在这些方面都需要进行考量,才能获得运行时间短、准确率高、稳定性好的预测模型。当然,在实际应用中,测试结果还可能包含更多项目,比如:模型中涉及的指标数量。
步骤108,判断获得的测试结果是否符合预设条件;若是,则执行步骤109;若否,则返回步骤103。
具体的说,可以本实施方式中的判断依据是步骤107的测试结果是否符合预设的条件。比如,准确率是否高于90%,由于准确率高于90%时,预测模型可以被认为是比较准确的,当然准确率能高于95%,对应的预测模型就相当成功了。又比如,涉及的指标数量是否合适,由于实际预测时,如果所涉及的指标数量越少,影响模型的稳定性,反之,则影响该模型的时间复杂度。
当然,如果判定测试结果还不能符合预设条件,可能所获得的模型还不够准确,那么就重复执行建模步骤,不断改进所涉及指标的组合。具体的说,每次执行步骤103至106后,对应获得一个COD的软测量模型。
比如,实际应用中,预定需要找到准确度在95%以上的COD软测量模型。第一次建模,选取了14个水质指标,在对第一次获得的COD软测量模型测试后,获得的准确度是90%;进行第二次建模,重新选取了7个水质指标,在对第二次获得的COD软测量模型测试后,获得的准确度是93%;再进行第三次建模,重新选取了3个水质指标,在对第三次获得的COD软测量模型测试后,获得的准确度是97%达到了准确度在95%以上的要求。
还需要说明的是,在多次执行建模步骤中,任意两次训练中所选取的水质监测指标组合不重复。
步骤109,将测试结果符合预设条件的COD的软测量模型作为最终的软测量模型。
具体的说,本发明人对最终的软测量模型的测试如图2所示,其中需要说明的是,图中线1表示COD的实际测量值,线2表示利用最终的软测量模型进行的COD预测值,线3表示两者间的相对误差。
本实施方式相对于现有技术而言,在COD的软测量预测中,运用多次建模并测试的方法,从众多水质监测指标中,选取到与待测水环境更为合适的指标,使得最终获得的模型可以更为准确稳定地预测出待测水环境中COD的值。解决了现有技术中,利用同一组水处理监测指标建模,所获得的模型不能准确符合各种不同条件的待测水环境,拓展本发明实施方式的应用场景。另外,由于采用梯度下降法,减少获得较佳预测模型的时间。
上面各种方法的步骤划分,只是为了描述清楚,实现时可以合并为一个步骤或者对某些步骤进行拆分,分解为多个步骤,只要包含相同的逻辑关系,都在本专利的保护范围内;对算法中或者流程中添加无关紧要的修改或者引入无关紧要的设计,但不改变其算法和流程的核心设计都在该专利的保护范围内。
本发明的第二实施方式涉及一种COD的水质软测量预测方法。其流程如图3所示,具体如下:
步骤201,数据的获取及分析。
具体的说,获取待测水环境中多个水质监测指标的数据值,每个水质监测指标被获取到的数据值有多个,多个水质监测指标中至少有一个为化学需氧量COD。
步骤202,指标选取。
具体的说,从除COD之外的各水质监测指标中选取若干个水质监测指标。
步骤203,数据样本集的准备。
具体的说,将所选的水质监测指标的数据值形成数据样本集,同时,将数据样本集划分为训练集及测试集。
步骤204,模型的建立。
具体的说,利用模糊神经网络算法对训练集进行训练,获得COD软测量模型。
步骤205,模型测试。
具体的说,利用测试集,对COD软测量模型进行测试,获得测试结果。
步骤206,判断获得的测试结果是否符合预设条件;若是,则执行步骤207;若否,则返回步骤202。
步骤207,将符合预设条件的COD软测量模型作为最终的软测量模型。
本实施方式相对于现有技术而言,主要区别及其效果在于:在COD的软测量预测中,运用多次建模并测试的方法,从众多水质监测指标中,选取到与待测水环境更为合适的指标,使得最终获得的模型可以更为准确稳定地预测出待测水环境中COD的值。解决了现有技术中,利用同一组水处理监测指标建模,所获得的模型不能准确符合待测水环境的COD预测,拓展本发明实施方式的应用场景。另外,进一步限定获取到的数据被分为训练集和测试集,两者的数据值可以为不同的数据值,使得测试时,可以利用集外测试,使COD软测量模型的效果说服力大大提高。
本领域的普通技术人员可以理解,上述各实施方式是实现本发明的具体实施例,而在实际应用中,可以在形式上和细节上对其作各种改变,而不偏离本发明的精神和范围。
Claims (6)
1.一种化学需氧量的水质软测量预测方法,其特征在于,包括:
数据的获取及分析:获取待测水环境中多个水质监测指标的数据值,每个水质监测指标被获取到的数据值有多个,多个所述水质监测指标中至少有一个为化学需氧量COD;
模型的建立:从除所述COD之外的各水质监测指标中选取若干个水质监测指标,将所选的水质监测指标的数据值形成数据样本集,将所述数据样本集划分为训练集及测试集,利用模糊神经网络算法对所述训练集进行训练,获得COD软测量模型;
模型测试:利用所述测试集,对所述COD软测量模型进行测试,获得测试结果;
重复执行所述模型的建立和所述模型测试,直至获得的测试结果符合预设条件,将符合所述预设条件的COD软测量模型作为最终的软测量模型。
2.根据权利要求1所述的化学需氧量的水质软测量预测方法,其特征在于,所述模型测试前,和所述模型测试后,还包括:通过梯度下降法对所获得的COD软测量模型进行优化;
所述对COD软测量模型进行测试中,对优化后的COD软测量模型进行测试。
3.根据权利要求1所述的化学需氧量的水质软测量预测方法,其特征在于,所述模型的建立之前,还包括:
相关性分析:利用各个水质监测指标的数据值,进行相关性分析,获得除所述COD之外的各个水质监测指标和所述COD之间的相关性;
相应的,所述选取若干个水质监测指标中,利用所述相关性的高低,选取若干个水质监测指标。
4.根据权利要求1所述的化学需氧量的水质软测量预测方法,其特征在于,所述模型测试中获得的所述测试结果包括:所述COD软测量模型的时间复杂度、准确率和稳定性。
5.根据权利要求1所述的化学需氧量的水质软测量预测方法,其特征在于,所述数据的获取及分析中,如果获取到的数据值属于:总氮、总磷或氨氮,则舍弃。
6.根据权利要求1至5中任意一项所述的化学需氧量的水质软测量预测方法,其特征在于,所述训练集中包含的数据值多于所述测试集中包含的数据值。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610619777.5A CN107665288A (zh) | 2016-07-30 | 2016-07-30 | 一种化学需氧量的水质软测量预测方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610619777.5A CN107665288A (zh) | 2016-07-30 | 2016-07-30 | 一种化学需氧量的水质软测量预测方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN107665288A true CN107665288A (zh) | 2018-02-06 |
Family
ID=61122326
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201610619777.5A Pending CN107665288A (zh) | 2016-07-30 | 2016-07-30 | 一种化学需氧量的水质软测量预测方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107665288A (zh) |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109975366A (zh) * | 2019-03-25 | 2019-07-05 | 浙江大学 | 农村生活污水a2o处理终端出水cod浓度软测量方法及装置 |
CN113159130A (zh) * | 2021-03-25 | 2021-07-23 | 中电建电力检修工程有限公司 | 一种施工污水处理的方法 |
CN115385445A (zh) * | 2022-09-21 | 2022-11-25 | 南京晓庄学院 | 一种基于大数据的污水生化处理用软测量管理系统 |
CN115824993A (zh) * | 2023-02-14 | 2023-03-21 | 北京英视睿达科技股份有限公司 | 水体化学需氧量的确定方法、装置、计算机设备及介质 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101334395A (zh) * | 2008-08-07 | 2008-12-31 | 北京工业大学 | Cod软测量的方法 |
CN102262147A (zh) * | 2011-07-15 | 2011-11-30 | 华南理工大学 | 一种废水处理系统出水cod的软测量方法和系统 |
CN104965971A (zh) * | 2015-05-24 | 2015-10-07 | 北京工业大学 | 一种基于模糊神经网络的氨氮浓度软测量方法 |
-
2016
- 2016-07-30 CN CN201610619777.5A patent/CN107665288A/zh active Pending
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101334395A (zh) * | 2008-08-07 | 2008-12-31 | 北京工业大学 | Cod软测量的方法 |
CN102262147A (zh) * | 2011-07-15 | 2011-11-30 | 华南理工大学 | 一种废水处理系统出水cod的软测量方法和系统 |
CN104965971A (zh) * | 2015-05-24 | 2015-10-07 | 北京工业大学 | 一种基于模糊神经网络的氨氮浓度软测量方法 |
Non-Patent Citations (3)
Title |
---|
朱立琴等: "《新建水库蓄水初期富营养化发生机理研究》", 31 January 2013 * |
谢福民: "基于支持向量机的软测量方法及其应用研究", 《中国优秀硕士学位论文全文数据库 信息科技辑》 * |
郭民等: "基于模糊神经网络的出水总磷软测量方法研究", 《计算机与应用化学》 * |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109975366A (zh) * | 2019-03-25 | 2019-07-05 | 浙江大学 | 农村生活污水a2o处理终端出水cod浓度软测量方法及装置 |
CN113159130A (zh) * | 2021-03-25 | 2021-07-23 | 中电建电力检修工程有限公司 | 一种施工污水处理的方法 |
CN115385445A (zh) * | 2022-09-21 | 2022-11-25 | 南京晓庄学院 | 一种基于大数据的污水生化处理用软测量管理系统 |
CN115824993A (zh) * | 2023-02-14 | 2023-03-21 | 北京英视睿达科技股份有限公司 | 水体化学需氧量的确定方法、装置、计算机设备及介质 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN104376380B (zh) | 一种基于递归自组织神经网络的氨氮浓度预测方法 | |
CN107664682A (zh) | 一种氨氮的水质软测量预测方法 | |
Heddam et al. | Predicting effluent biochemical oxygen demand in a wastewater treatment plant using generalized regression neural network based approach: a comparative study | |
Feng et al. | A risk assessment model of water shortage based on information diffusion technology and its application in analyzing carrying capacity of water resources | |
Semiromi et al. | Water quality index development using fuzzy logic: A case study of the Karoon River of Iran | |
Han et al. | Data-driven intelligent monitoring system for key variables in wastewater treatment process | |
CN111291937A (zh) | 基于支持向量分类与gru神经网络联合的处理污水水质预测方法 | |
CN108564229A (zh) | 一种基于arima模型的工业污水进水量预测的方法 | |
CN107358021A (zh) | 一种基于优化bp神经网络的do预测模型建立方法 | |
CN104182794B (zh) | 一种基于神经网络的污水处理过程出水总磷tp软测量方法 | |
CN100514060C (zh) | 水体化学需氧量在线虚拟监测方法 | |
CN109344971B (zh) | 一种基于自适应递归模糊神经网络的出水氨氮浓度预测方法 | |
CN104360035A (zh) | 一种基于自组织粒子群-径向基神经网络的污水总磷tp软测量方法 | |
CN107402586A (zh) | 基于深度神经网络的溶解氧浓度控制方法及系统 | |
CN103632032A (zh) | 一种城市污水处理过程出水指标在线软测量预测方法 | |
CN107665288A (zh) | 一种化学需氧量的水质软测量预测方法 | |
CN109657790A (zh) | 一种基于pso的递归rbf神经网络出水bod预测方法 | |
CN107664683A (zh) | 一种总氮的水质软测量预测方法 | |
CN115392617A (zh) | 一种基于环境异质性的水库库尾河段水环境安全评估方法 | |
CN105160422B (zh) | 基于自组织级联神经网络的污水处理出水总磷预测方法 | |
Pan et al. | A new approach to estimating oxygen off-gas fraction and dynamic alpha factor in aeration systems using hybrid machine learning and mechanistic models | |
CN107665363A (zh) | 一种总磷的水质软测量预测方法 | |
CN109019717A (zh) | 一种火电厂脱硫废水智能处理方法及系统 | |
CN109359741A (zh) | 一种废水处理进水水质时序变化智能预测方法 | |
Fan et al. | Neural-network-based water quality monitoring for wastewater treatment processes |
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 | ||
RJ01 | Rejection of invention patent application after publication |
Application publication date: 20180206 |
|
RJ01 | Rejection of invention patent application after publication |