CN107665363A - 一种总磷的水质软测量预测方法 - Google Patents
一种总磷的水质软测量预测方法 Download PDFInfo
- Publication number
- CN107665363A CN107665363A CN201610615700.0A CN201610615700A CN107665363A CN 107665363 A CN107665363 A CN 107665363A CN 201610615700 A CN201610615700 A CN 201610615700A CN 107665363 A CN107665363 A CN 107665363A
- Authority
- CN
- China
- Prior art keywords
- water quality
- model
- quality monitoring
- soft measurement
- measurement model
- 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 149
- 238000005259 measurement Methods 0.000 title claims abstract description 109
- 238000000034 method Methods 0.000 title claims abstract description 39
- OAICVXFJPJFONN-UHFFFAOYSA-N Phosphorus Chemical compound [P] OAICVXFJPJFONN-UHFFFAOYSA-N 0.000 title claims abstract description 19
- 229910052698 phosphorus Inorganic materials 0.000 title claims abstract description 19
- 239000011574 phosphorus Substances 0.000 title claims abstract description 19
- 238000012360 testing method Methods 0.000 claims abstract description 94
- 238000012544 monitoring process Methods 0.000 claims abstract description 90
- 238000012549 training Methods 0.000 claims abstract description 31
- 238000013528 artificial neural network Methods 0.000 claims abstract description 16
- 238000004422 calculation algorithm Methods 0.000 claims abstract description 11
- XKMRRTOUMJRJIA-UHFFFAOYSA-N ammonia nh3 Chemical compound N.N XKMRRTOUMJRJIA-UHFFFAOYSA-N 0.000 claims description 41
- IJGRMHOSHXDMSA-UHFFFAOYSA-N Atomic nitrogen Chemical compound N#N IJGRMHOSHXDMSA-UHFFFAOYSA-N 0.000 claims description 24
- 238000010219 correlation analysis Methods 0.000 claims description 11
- 229910052757 nitrogen Inorganic materials 0.000 claims description 10
- 238000004458 analytical method Methods 0.000 claims description 9
- QVGXLLKOCUKJST-UHFFFAOYSA-N atomic oxygen Chemical compound [O] QVGXLLKOCUKJST-UHFFFAOYSA-N 0.000 claims description 9
- 229910052760 oxygen Inorganic materials 0.000 claims description 9
- 239000001301 oxygen Substances 0.000 claims description 9
- 239000000126 substance Substances 0.000 claims description 8
- 238000011478 gradient descent method Methods 0.000 claims description 7
- 230000007613 environmental effect Effects 0.000 abstract description 2
- 239000010865 sewage Substances 0.000 description 12
- 238000004364 calculation method Methods 0.000 description 7
- 230000006872 improvement Effects 0.000 description 6
- 230000000694 effects Effects 0.000 description 5
- 230000006870 function Effects 0.000 description 5
- 230000008569 process Effects 0.000 description 5
- 230000008901 benefit Effects 0.000 description 4
- 239000007787 solid Substances 0.000 description 3
- 238000005273 aeration Methods 0.000 description 2
- 238000012937 correction Methods 0.000 description 2
- 230000002596 correlated effect Effects 0.000 description 2
- 238000013461 design Methods 0.000 description 2
- 238000001514 detection method Methods 0.000 description 2
- 238000011161 development Methods 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 238000005457 optimization Methods 0.000 description 2
- 238000012545 processing Methods 0.000 description 2
- 238000010992 reflux Methods 0.000 description 2
- 238000003911 water pollution Methods 0.000 description 2
- 230000033228 biological regulation Effects 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 239000003153 chemical reaction reagent Substances 0.000 description 1
- 230000000295 complement effect Effects 0.000 description 1
- 230000001276 controlling effect Effects 0.000 description 1
- 230000000875 corresponding effect Effects 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
- 238000003912 environmental pollution Methods 0.000 description 1
- 238000012851 eutrophication Methods 0.000 description 1
- 238000012423 maintenance Methods 0.000 description 1
- 230000014759 maintenance of location Effects 0.000 description 1
- 238000013507 mapping Methods 0.000 description 1
- 230000033116 oxidation-reduction process Effects 0.000 description 1
- 238000002360 preparation method Methods 0.000 description 1
- 238000003825 pressing Methods 0.000 description 1
- 238000000746 purification Methods 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 238000005070 sampling Methods 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
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06N—COMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
- G06N3/00—Computing arrangements based on biological models
- G06N3/02—Neural networks
- G06N3/04—Architecture, e.g. interconnection topology
- G06N3/043—Architecture, e.g. interconnection topology based on fuzzy logic, fuzzy membership or fuzzy inference, e.g. adaptive neuro-fuzzy inference systems [ANFIS]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06N—COMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
- G06N3/00—Computing arrangements based on biological models
- G06N3/02—Neural networks
- G06N3/08—Learning methods
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q10/00—Administration; Management
- G06Q10/04—Forecasting or optimisation specially adapted for administrative or management purposes, e.g. linear programming or "cutting stock problem"
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q50/00—Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
- G06Q50/10—Services
- G06Q50/26—Government or public services
-
- G—PHYSICS
- G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
- G16Z—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS, NOT OTHERWISE PROVIDED FOR
- G16Z99/00—Subject matter not provided for in other main groups of this subclass
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2219/00—Indexing scheme relating to application aspects of data processing equipment or methods
- G06F2219/10—Environmental application, e.g. waste reduction, pollution control, compliance with environmental legislation
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
- Y02A20/00—Water conservation; Efficient water supply; Efficient water use
- Y02A20/152—Water filtration
Landscapes
- Engineering & Computer Science (AREA)
- Business, Economics & Management (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- General Health & Medical Sciences (AREA)
- Tourism & Hospitality (AREA)
- Strategic Management (AREA)
- Human Resources & Organizations (AREA)
- Software Systems (AREA)
- Health & Medical Sciences (AREA)
- Economics (AREA)
- Mathematical Physics (AREA)
- Development Economics (AREA)
- Biophysics (AREA)
- Computational Linguistics (AREA)
- Data Mining & Analysis (AREA)
- Evolutionary Computation (AREA)
- General Business, Economics & Management (AREA)
- Molecular Biology (AREA)
- Computing Systems (AREA)
- General Engineering & Computer Science (AREA)
- Marketing (AREA)
- Life Sciences & Earth Sciences (AREA)
- Artificial Intelligence (AREA)
- Biomedical Technology (AREA)
- Mathematical Optimization (AREA)
- Automation & Control Theory (AREA)
- Computational Mathematics (AREA)
- Pure & Applied Mathematics (AREA)
- Mathematical Analysis (AREA)
- Primary Health Care (AREA)
- Educational Administration (AREA)
- Fuzzy Systems (AREA)
- Game Theory and Decision Science (AREA)
- Entrepreneurship & Innovation (AREA)
- Operations Research (AREA)
- Quality & Reliability (AREA)
- Sampling And Sample Adjustment (AREA)
Abstract
本发明实施方式涉及环境监测领域,公开了一种总磷的水质软测量预测方法。包括:获取待测水环境中多个水质监测指标的数据值,每个水质监测指标被获取到的数据值有多个,多个水质监测指标中至少有一个为TP;从除TP之外的各水质监测指标中选取若干个水质监测指标,将所选的水质监测指标的数据值形成数据样本集,将数据样本集划分为训练集及测试集,利用模糊神经网络算法对训练集进行训练,获得TP软测量模型;利用测试集,对TP软测量模型进行测试,获得测试结果;重复执行模型的建立和模型测试,直至获得的测试结果符合预设条件,将符合预设条件的TP软测量模型作为最终的软测量模型。所得模型预测效果精确,降低TP的水质软测量预测成本。
Description
技术领域
本发明涉及环境监测领域,特别涉及总磷的水质软测量预测技术。
背景技术
众所周知,我国水资源十分稀缺,人均占有量仅占世界平均水平的1/4,同时严重的水环境的污染问题,如水污染、水体富营养化、城市黑臭水体、地下水污染等。使水资源问题愈加突出。在农村,我国将近96%的、约80~90亿t/a的农村污水没有治理,严重农村环境污染和河海流域污染,在城市,黑臭水体已然成为一种城市病,所有城市几乎无一幸免。国内大部分河流均受到不同程度的污染。缓解水资源紧张、加大污水净化、减轻水环境污染刻不容缓。
随着水环境保护及污染治理的紧迫,水质监测提供的水质信息尤为重要,水质监测是水环境保护及污染治理的重要基础,其重要性关乎掌握水资源质量状况、以及对用水、污水处理、排污等有效监控。一方面,受污染水体的水质情况都需要通过化学需氧量(COD)、氨氮、总氮(TN)、总磷(TP)等关键水质指标的监测提供准确及时的水质动态信息,另一方面,及时有效的污水处理参数监测,对污水处理系统有重大意义,污水的排放必须符合国家污水排放标准中的相关规定,这就要求我们必须检测处理出水中TN、TP、COD、氨氮等关键参数。因此加强水环境保护及污染治理需要水质监测工作超前发展。
COD、氨氮、TN、TP等关键水质指标的监测,对于水质处理的控制、过程的优化及诊断起着重要的作用,然而这类水质指标难以测量或不易在线测量,目前主要有人工实验室检测及检测水质的自动分析仪表(水质在线检测仪)。
水质在线自动分析仪,如COD在线监测仪、TP在线检测仪等,在我国发展比较晚,国内产品的测量精确度和可靠性与国外存在一定的差距,存在品种单一、测量精准度低,测量周期长等缺陷,没有完全达到满足测量需求,污水生物处理过程的水质参数无法准确地进行测量。但使用国外的产品存在检修维护困难、价格昂贵等问题,限制了其在水环境保护及污染治理领域的应用。
目前国内污水处理厂的出水水质关键参数如化学需氧量(COD)、氨氮、总氮(TN)、总磷(TP)等水指标大多依靠实验室人工化验检测来确知,显然,实验室人工化验存在滞后性,从采样到输出测定结果,有时间差值,时间滞后;精确度受人为影响,较难判断;无法实时监测,不能及时依据水质变化指导,难以控制污水处理情况。
目前TP监测解决思路有以下两种:
(1)改进直接测量仪表:按照传统的检测技术发展思路,以硬件形式研制新型的过程测量仪表。但由于污水处理系统是一个复杂的生物化学处理系统,是典型的非线性、多变量、不稳定、大时滞系统,干扰因素众多且不确定,直接研发、改进在线监测仪表,难度大,且监测成本(一次性投入及试剂费用)较高,难以降低。
(2)间接测量:采用间接测量的思路,利用容易获取并且与被测变量相关的测量信息,通过计算估计被测变量。软测量即是这一思想的集中体现,将软测量应用于水环境保护及污染治理中,可实现投资少、实时监测的目的。
发明内容
本发明实施方式的目的在于提供一种总磷的水质软测量预测方法,在保证TP的水质软测量预测准确度的同时,降低TP的水质软测量预测成本。
为解决上述技术问题,本发明的实施方式提供了一种总磷的水质软测量预测方法,包括:
数据的获取及分析:获取待测水环境中多个水质监测指标的数据值,每个水处理监测指标被获取到的数据值有多个,多个所述水质监测指标中至少有一个为总磷TP;
模型的建立:从除所述TP之外的各水质监测指标中选取若干个水质监测指标,将所选的水质监测指标的数据值形成数据样本集,将所述数据样本集划分为训练集及测试集,利用模糊神经网络算法对所述训练集进行训练,获得TP软测量模型;
模型测试:利用所述测试集,通过多次仿真实验的测试,对所述TP软测量模型进行测试,获得测试结果;
重复执行所述模型的建立和所述模型测试,直至获得的测试结果符合预设条件,将符合所述预设条件的TP软测量模型作为最终的软测量模型。
本发明实施方式相对于现有技术而言,主要区别及其效果在于:在TP的软测量预测中,运用多次建模并测试的方法,从众多水质监测指标中,选取到与待测水环境更为合适的指标,使得最终获得的模型可以更为准确稳定地预测出待测水环境中TP的值。解决了现有技术中,利用同一组水处理监测指标建模,所获得的模型不能准确符合待测水环境的TP预测,拓展本发明实施方式的应用场景。另外,进一步限定获取到的数据被分为训练集和测试集,两者的数据值可以为不同的数据值,使得测试时,可以利用集外测试,使TP软测量模型的效果说服力大大提高。
作为进一步改进,所述建模步骤之前,还包括:相关性分析:利用各个水质监测指标的数据值,进行相关性分析,获得除所述TP之外的各个水质监测指标和所述TP之间的相关性;相应的,所述选取若干个水质监测指标中,利用所述相关性的高低,选取若干个水质监测指标。
利用相关性分析,使得选取指标时可以根据分析结果选取,使得选取具有了较佳的依据,提高获得符合要求的软测量模型的效率。
作为进一步改进,所述模型测试后,和所述模型测试前,还包括:通过梯度下降法对所获得的TP软测量模型进行优化;所述对TP软测量模型进行测试中,对优化后的TP软测量模型进行测试。
利用梯度下降法进行模型优化,使得模型的误差得到大幅度降低,大大提高模型精确度,通过上述模糊神经网络算法和梯度下降法,使获得的最终模型更为精确稳定。
作为进一步改进,所述根据各测试结果选定最终的软测量模型中,所述模型测试中获得的所述测试结果包括:所述TP软测量模型的时间复杂度、准确率和稳定性。进一步限定测试结果的内容,提高选定的最终的软测量模型的有效性。
作为进一步改进,所述获取步骤中,如果获取到的数据值属于:总磷、总磷或氨氮,则舍弃。由于总磷、总磷或氨氮这三个指标的直接测试成本过高,不采用这三个指标作为输入,进一步降低TP的水质软测量预测成本。
作为进一步改进,所述数据的获取中,包括:利用预设的化学需氧量COD的软测量模型所预测的所述COD的预测值,和利用预设的氨氮的软测量模型所预测的所述氨氮的预测值,以及利用预设的总氮TN的软测量模型所预测的所述TN的预测值;所述COD的预测值、所述氨氮和所述TN的预测值分别有多个;所形成的数据样本集中,还包括:所述COD的预测值、所述氨氮的预测值和所述TN的预测值。避免使用仪器直接测量COD、TN和氨氮,改用这三个指标的预测值参与模型的生成,使得在保证准确度时,模型的生成成本更低。
附图说明
图1是根据本发明第一实施方式中的TP的水质软测量预测方法流程图;
图2是根据本发明第一实施方式中的实际出水TP值与软测量所得TP值的对比及相对误差的示意图;
图3是根据本发明第二实施方式中的TP的水质软测量预测方法流程图;
图4是根据本发明第三实施方式中的TP的软测量预测方法流程图。
具体实施方式
为使本发明实施方式的目的、技术方案和优点更加清楚,下面将结合附图对本发明的各实施方式进行详细的阐述。然而,本领域的普通技术人员可以理解,在本发明各实施方式中,为了使读者更好地理解本申请而提出了许多技术细节。但是,即使没有这些技术细节和基于以下各实施方式的种种变化和修改,也可以实现本申请所要求保护的技术方案。
本发明的第一实施方式涉及一种TP的水质软测量预测方法。其流程如图1所示,具体如下:
步骤101,获取待测水环境中多个水质监测指标的数据值。
具体的说,本实施方式中的水质监测指标可以是以下指标:TP、进水流量、出水流量、进水PH、出水PH、出水SS(即固体悬浮物)、曝气池MLSS(即混合液悬浮固体浓度)、回流MLSS、MLVSS(即混合液挥发性悬浮固体浓度)、曝气SVI(即污泥体积指数)、回流SVI、1#池DO(即溶解氧)、2#池DO、ORP(即氧化还原电位)、温度、HRT(即水力停留时间)。当然,在实际应用中,还可以选取其他水质监测指标,不限于上述指标。更具体的说,获取的水质监测指标中至少一个为总磷TP。
更具体的说,上述水质监测指标将划分为两类,分别是输入指标及输出指标,输出指标指水质监测指标TP,TP的数据值作为输出数据样本;输入指标包括除TP外的其他水质监测指标,这些指标的数据值作为输入数据样本。
需要说明的是,每个水质监测指标被获取到的数据值有多个。一般会选择近阶段的历史数据作为各指标的数据值。
值得一提的是,本步骤中还可以进一步筛选数据,比如:如果获取到的数据值属于:总氮、COD或氨氮,则舍弃。这样就可以不采用总氮、COD或氨氮这三个指标,由于总氮、COD或氨氮这三个指标,在利用在线测量仪表获取时,测试成本过高,不采用这三个指标作为输入,进一步降低TP的水质软测量预测成本。
步骤102,利用各个水质监测指标的数据值,进行相关性分析,获得除TP之外的各个水质监测指标和TP之间的相关性。
具体的说,相关性分析公式为下式(1):
其中,相关系数r表示两个变量x和y之间的关系密切程度。
值得一提的是,在实际应用中,也可以不进行相关分析,直接利用经验选取,但本实施方式还是增加了相关性分析的步骤,可以使得选取指标时可以根据分析结果选取,使得选取具有了较佳的依据,加快获得符合要求的软测量模型的速度。
步骤103,从除TP之外的各水质监测指标中选取若干个水质监测指标组合,形成数据样本集。
具体的说,本实施方式中的选取依据可以是步骤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,将数据样本集划分为两类:一类作为训练集,一类作为测试集。
具体的说,本实施方式划分后的训练集中包含的数据值可以多于测试集中包含的数据值。其中,划分时可以采用各种原则。
值得一提的是,后续建模步骤中使用的数据样本集为训练集;测试步骤中所使用的数据样本集为测试集。进一步限定训练集和测试集的数据值为不同的数据值,使得测试时,可以利用集外测试,使TP软测量模型的效果说服力大大提高。
步骤105,利用模糊神经网络算法对训练集进行训练,获得TP的软测量模型。
具体的说,尽管模糊逻辑系统和神经网络系统在结构上截然不同,但是在两者之间可以实现优势互补。神经网络具有非线性映射能力、学习能力、并行处理能力和容错能力。模糊逻辑具有处理不确定性的能力、该发明将两者的优势有机结合起来,组成在功能上更加完善和强大的系统应用到污水水质监测中。即采用模糊神经网络(FNN)建立TP的软测量模型。
更具体的说,本实施方式中建立的模糊神经网络结构分为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,利用梯度下降法对获得TP软测量模型进行优化,获得优化后的TP水质软测量模型。
具体的说,梯度下降法优化计算包括
(1)误差计算
式(5)中,yd是网络期望输出;yc是网络实际输出;e为期望输出和实际输出的误差。
(2)系数修正
式(6)和(7)中,为神经网络系数;α为网络学习率;xj为网络输入参数;wi为输入参数隶属度连乘积。
(3)参数修正
式(8)和(9)中,分别为隶属度函数的中心和宽度。
步骤107,利用测试集对所获得优化后的TP的软测量模型进行测试,获得测试结果。
具体的说,本实施方式中测试结果包括:软测量模型的时间复杂度、准确率和稳定性。为了获得更佳的预测模型,所以在这些方面都需要进行考量,才能获得运行时间短、准确率高、稳定性好的预测模型。当然,在实际应用中,测试结果还可能包含更多项目,比如:模型中涉及的指标数量。
步骤108,判断获得的测试结果是否符合预设条件;若是,则执行步骤109;若否,则返回步骤103。
具体的说,可以本实施方式中的判断依据是步骤107的测试结果是否符合预设的条件。比如,准确率是否高于90%,由于准确率高于90%时,预测模型可以被认为是比较准确的,当然准确率能高于95%,对应的预测模型就相当成功了。又比如,涉及的指标数量是否合适,由于实际预测时,如果所涉及的指标数量越少,影响模型的稳定性,反之,则影响该模型的时间复杂度。
当然,如果判定测试结果还不能符合预设条件,可能所获得的模型还不够准确,那么就重复执行建模步骤,不断改进所涉及指标的组合。具体的说,每次执行步骤103至106后,对应获得一个TP的软测量模型。
比如,实际应用中,预定需要找到准确度在95%以上的TP软测量模型。第一次建模,选取了14个水质指标,在对第一次获得的TP软测量模型测试后,获得的准确度是90%;进行第二次建模,重新选取了7个水质指标,在对第二次获得的TP软测量模型测试后,获得的准确度是93%;再进行第三次建模,重新选取了3个水质指标,在对第三次获得的TP软测量模型测试后,获得的准确度是97%达到了准确度在95%以上的要求。
还需要说明的是,在多次执行建模步骤中,任意两次训练中所选取的水质监测指标组合不重复。
步骤109,将测试结果符合预设条件的TP的软测量模型作为最终的软测量模型。
具体的说,本发明人对最终的软测量模型的测试如图2所示,其中需要说明的是,图中线1表示TP的实际测量值,线2表示利用最终的软测量模型进行测试所获得的TP预测值,线3表示两者间的相对误差。
本实施方式相对于现有技术而言,在TP的软测量预测中,运用多次建模并测试的方法,从众多水质监测指标中,选取到与待测水环境更为合适的指标,使得最终获得的模型可以更为准确稳定地预测出待测水环境中TP的值。解决了现有技术中,利用同一组水处理监测指标建模,所获得的模型不能准确符合各种不同条件的待测水环境,拓展本发明实施方式的应用场景。另外,由于采用梯度下降法,减少获得较佳预测模型的时间。
上面各种方法的步骤划分,只是为了描述清楚,实现时可以合并为一个步骤或者对某些步骤进行拆分,分解为多个步骤,只要包含相同的逻辑关系,都在本专利的保护范围内;对算法中或者流程中添加无关紧要的修改或者引入无关紧要的设计,但不改变其算法和流程的核心设计都在该专利的保护范围内。
本发明的第二实施方式涉及一种TP的水质软测量预测方法。其流程如图3所示,具体如下:
步骤201,数据的获取及分析。
具体的说,获取待测水环境中多个水质监测指标的数据值,每个水质监测指标被获取到的数据值有多个,多个水质监测指标中至少有一个为总磷TP。
步骤202,指标选取。
具体的说,从除TP之外的各水质监测指标中选取若干个水质监测指标。
步骤203,数据样本集的准备。
具体的说,将所选的水质监测指标的数据值形成数据样本集,同时,将数据样本集划分为训练集及测试集。
步骤204,模型的建立。
具体的说,利用模糊神经网络算法对训练集进行训练,获得TP软测量模型。
步骤205,模型测试。
具体的说,利用测试集,对TP软测量模型进行测试,获得测试结果。
步骤206,判断获得的测试结果是否符合预设条件;若是,则执行步骤207;若否,则返回步骤202。
步骤207,将符合预设条件的TP软测量模型作为最终的软测量模型。
本实施方式相对于现有技术而言,主要区别及其效果在于:在TP的软测量预测中,运用多次建模并测试的方法,从众多水质监测指标中,选取到与待测水环境更为合适的指标,使得最终获得的模型可以更为准确稳定地预测出待测水环境中TP的值。解决了现有技术中,利用同一组水处理监测指标建模,所获得的模型不能准确符合待测水环境的TP预测,拓展本发明实施方式的应用场景。另外,进一步限定获取到的数据被分为训练集和测试集,两者的数据值可以为不同的数据值,使得测试时,可以利用集外测试,使TP软测量模型的效果说服力大大提高。
本发明的第三实施方式还涉及一种总磷的水质软测量预测方法。本实施方式是在第一实施方式上做的进一步改进,主要改进之处在于,在软测量模型生成时,参与生成的原始数据除了直接获取到的各个水质监测指标的数据值,还包括了COD的预测值、总氮和氨氮的预测值。本发明发明人发现,TN的软测量模型中如果包括COD、总氮和氨氮,模型的准确度更高,但为了控制成本,避免使用仪器直接测量COD、总氮和氨氮,改用这三个指标的预测值参与模型的生成,使得在保证准确度时,模型的生成成本更低。
本实施方式的流程如图4所示,具体如下:
步骤401和第一实施方式中的步骤101相类似,在此不再赘述。
步骤402,分别获取COD、总氮和氨氮的预测值。
具体的说,利用预设的化学需氧量COD的软测量模型所预测的COD的预测值,和利用预设的氨氮的软测量模型所预测的氨氮的预测值。更具体的说,COD的预测值和氨氮的预测值分别有多个。
值得一提的是,COD的软测量模型利用以下方法获得:
数据的获取及分析:获取待测水环境中多个水质监测指标的数据值,每个水质监测指标被获取到的数据值有多个,多个水质监测指标中至少有一个为化学需氧量COD;
模型的建立:从除COD之外的各水质监测指标中选取若干个水质监测指标,将所选的水质监测指标的数据值形成数据样本集,将数据样本集划分为训练集及测试集,利用模糊神经网络算法对训练集进行训练,获得COD软测量模型;
模型测试:利用测试集,对COD软测量模型进行测试,获得测试结果;
重复执行模型的建立和模型测试,直至获得的测试结果符合预设条件,将符合预设条件的COD软测量模型作为最终的软测量模型。
还需说明的是,总氮的软测量模型利用以下方法获得:
数据的获取及分析:获取待测水环境中多个水质监测指标的数据值,每个水质监测指标被获取到的数据值有多个,多个水质监测指标中至少有一个为总氮TN;
模型的建立:从除TN之外的各水质监测指标中选取若干个水质监测指标,将所选的水质监测指标的数据值形成数据样本集,将数据样本集划分为训练集及测试集,利用模糊神经网络算法对训练集进行训练,获得TN软测量模型;
模型测试:利用测试集,对TN软测量模型进行测试,获得测试结果;
重复执行模型的建立和模型测试,直至获得的测试结果符合预设条件,将符合预设条件的TN软测量模型作为最终的软测量模型。
还需一提的是,氨氮的软测量模型利用以下方法获得:
数据的获取及分析:获取待测水环境中多个水质监测指标的数据值,每个水质监测指标被获取到的数据值有多个,多个水质监测指标中至少有一个为氨氮;
模型的建立:从除氨氮之外的各水质监测指标中选取若干个水质监测指标,将所选的水质监测指标的数据值形成数据样本集,将数据样本集划分为训练集及测试集,利用模糊神经网络算法对训练集进行训练,获得氨氮软测量模型;
模型测试:利用测试集,对氨氮软测量模型进行测试,获得测试结果;
重复执行模型的建立和模型测试,直至获得的测试结果符合预设条件,将符合预设条件的氨氮软测量模型作为最终的软测量模型。
可见,利用上述COD的软测量模型、TN的软测量模型和氨氮的软测量模型对待测水环境进行预测,就可以获取到COD的预测值和氨氮的预测值。
步骤403,利用各个水处理监测指标的数据值、COD预测值、TN预测值和氨氮预测值进行相关性分析。
步骤404,选取若干个指标、COD的预测值、TN的预测值和氨氮的预测值,形成数据样本集。
具体的说,从除TP之外的各水质监测指标中选取若干个水质监测指标,将所选的水质监测指标的数据值、COD的预测值、TN的预测值和氨氮的预测值形成数据样本集。也就是说,所形成的数据样本集中,还包括:COD的预测值、TN的预测值和氨氮的预测值。
步骤405,将数据样本集划分为两类。
具体的说,是将包括了COD的预测值、TN的预测值和氨氮的预测值的数据样本集划分为两类。
步骤406至步骤410和第一实施方式中的步骤105至步骤109相类似,在此不再赘述。
本实施方式和现有技术相比,利用软测量模型预测了COD、TN和氨氮的值,使得COD、TN和氨氮得以参与TP的软测量模型的生成,提高了获得的TP软测量模型的准确性,也由于避免使用了在线测量仪器获取的数据,控制了TP软测量模型的生成成本。
本领域的普通技术人员可以理解,上述各实施方式是实现本发明的具体实施例,而在实际应用中,可以在形式上和细节上对其作各种改变,而不偏离本发明的精神和范围。
Claims (7)
1.一种总磷的水质软测量预测方法,其特征在于,包括:
数据的获取及分析:获取待测水环境中多个水质监测指标的数据值,每个水质监测指标被获取到的数据值有多个,多个所述水质监测指标中至少有一个为总磷TP;
模型的建立:从除所述TP之外的各水质监测指标中选取若干个水质监测指标,将所选的水质监测指标的数据值形成数据样本集,将所述数据样本集划分为训练集及测试集,利用模糊神经网络算法对所述训练集进行训练,获得TP软测量模型;
模型测试:利用所述测试集,对所述TP软测量模型进行测试,获得测试结果;
重复执行所述模型的建立和所述模型测试,直至获得的测试结果符合预设条件,将符合所述预设条件的TP软测量模型作为最终的软测量模型。
2.根据权利要求1所述的总磷的水质软测量预测方法,其特征在于,所述模型测试前,和所述模型测试后,还包括:通过梯度下降法对所获得的TP软测量模型进行优化;
所述对TP软测量模型进行测试中,对优化后的TP软测量模型进行测试。
3.根据权利要求1所述的总磷的水质软测量预测方法,其特征在于,所述模型的建立之前,还包括:
相关性分析:利用各个水质监测指标的数据值,进行相关性分析,获得除所述TP之外的各个水质监测指标和所述TP之间的相关性;
相应的,所述选取若干个水质监测指标中,利用所述相关性的高低,选取若干个水质监测指标。
4.根据权利要求1所述的总磷的水质软测量预测方法,其特征在于,所述模型测试中获得的所述测试结果包括:所述TP软测量模型的时间复杂度、准确率和稳定性。
5.根据权利要求1所述的总磷的水质软测量预测方法,其特征在于,所述数据的获取及分析中,利用在线测量仪表获取各所述水质监测指标的数据值;如果获取到的数据值属于:化学需氧量、总氮或氨氮,则舍弃。
6.根据权利要求1所述的总磷的水质软测量预测方法,其特征在于,所述数据的获取中,包括:利用预设的化学需氧量COD的软测量模型所预测的所述COD的预测值,和利用预设的氨氮的软测量模型所预测的所述氨氮的预测值,以及利用预设的总氮TN的软测量模型所预测的所述TN的预测值;所述COD的预测值、所述氨氮和所述TN的预测值分别有多个;
所形成的数据样本集中,还包括:所述COD的预测值、所述氨氮的预测值和所述TN的预测值。
7.根据权利要求1至6中任意一项所述的总磷的水质软测量预测方法,其特征在于,所述训练集中包含的数据值多于所述测试集中包含的数据值。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610615700.0A CN107665363A (zh) | 2016-07-30 | 2016-07-30 | 一种总磷的水质软测量预测方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610615700.0A CN107665363A (zh) | 2016-07-30 | 2016-07-30 | 一种总磷的水质软测量预测方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN107665363A true CN107665363A (zh) | 2018-02-06 |
Family
ID=61114604
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201610615700.0A Pending CN107665363A (zh) | 2016-07-30 | 2016-07-30 | 一种总磷的水质软测量预测方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107665363A (zh) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109534486A (zh) * | 2018-12-21 | 2019-03-29 | 哈尔滨工程大学 | 一种基于随机递归小波神经网络的船舶污水处理控制预测系统及预测方法 |
CN109934419A (zh) * | 2019-04-11 | 2019-06-25 | 苏州热工研究院有限公司 | 一种核电厂取水口海洋生物量变化预测方法 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103235096A (zh) * | 2013-04-16 | 2013-08-07 | 广州铁路职业技术学院 | 污水水质检测方法与装置 |
CN103886369A (zh) * | 2014-03-27 | 2014-06-25 | 北京工业大学 | 一种基于模糊神经网络的出水总磷tp预测方法 |
CN104182794A (zh) * | 2014-09-05 | 2014-12-03 | 北京工业大学 | 一种基于神经网络的污水处理过程出水总磷tp软测量方法 |
CN104965971A (zh) * | 2015-05-24 | 2015-10-07 | 北京工业大学 | 一种基于模糊神经网络的氨氮浓度软测量方法 |
CN105574326A (zh) * | 2015-12-12 | 2016-05-11 | 北京工业大学 | 一种基于自组织模糊神经网络的出水氨氮浓度软测量方法 |
-
2016
- 2016-07-30 CN CN201610615700.0A patent/CN107665363A/zh active Pending
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103235096A (zh) * | 2013-04-16 | 2013-08-07 | 广州铁路职业技术学院 | 污水水质检测方法与装置 |
CN103886369A (zh) * | 2014-03-27 | 2014-06-25 | 北京工业大学 | 一种基于模糊神经网络的出水总磷tp预测方法 |
CN104182794A (zh) * | 2014-09-05 | 2014-12-03 | 北京工业大学 | 一种基于神经网络的污水处理过程出水总磷tp软测量方法 |
CN104965971A (zh) * | 2015-05-24 | 2015-10-07 | 北京工业大学 | 一种基于模糊神经网络的氨氮浓度软测量方法 |
CN105574326A (zh) * | 2015-12-12 | 2016-05-11 | 北京工业大学 | 一种基于自组织模糊神经网络的出水氨氮浓度软测量方法 |
Non-Patent Citations (1)
Title |
---|
郭民: "基于模糊神经网络的出水总磷软测量方法研究", 《计算机与应用化学》 * |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109534486A (zh) * | 2018-12-21 | 2019-03-29 | 哈尔滨工程大学 | 一种基于随机递归小波神经网络的船舶污水处理控制预测系统及预测方法 |
CN109934419A (zh) * | 2019-04-11 | 2019-06-25 | 苏州热工研究院有限公司 | 一种核电厂取水口海洋生物量变化预测方法 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107664682A (zh) | 一种氨氮的水质软测量预测方法 | |
Han et al. | Data-driven intelligent monitoring system for key variables in wastewater treatment process | |
CN107665288A (zh) | 一种化学需氧量的水质软测量预测方法 | |
CN110378533B (zh) | 一种基于大数据分析的智能曝气管理方法 | |
Tian et al. | An optimization of artificial neural network model for predicting chlorophyll dynamics | |
CN103886369B (zh) | 一种基于模糊神经网络的出水总磷tp预测方法 | |
Honggui et al. | A fuzzy neural network approach for online fault detection in waste water treatment process | |
Xu et al. | Data-driven prediction of neutralizer pH and valve position towards precise control of chemical dosage in a wastewater treatment plant | |
Han et al. | Hierarchical neural network modeling approach to predict sludge volume index of wastewater treatment process | |
CN103606006B (zh) | 基于自组织t‑s模糊神经网络的污泥沉降指数软测量方法 | |
CN102262147A (zh) | 一种废水处理系统出水cod的软测量方法和系统 | |
CN110981021B (zh) | 一种基于模糊bp神经网络的废水高级氧化处理智能加药系统及方法 | |
CN103632032A (zh) | 一种城市污水处理过程出水指标在线软测量预测方法 | |
CN203772781U (zh) | 一种基于特征变量的污水总磷测量装置 | |
CN111398539A (zh) | 一种基于大数据和分子生物技术的水质微生物指示方法 | |
CN107664683A (zh) | 一种总氮的水质软测量预测方法 | |
Wang et al. | Prediction of effluent ammonia nitrogen in wastewater treatment plant based on self-organizing hybrid neural network | |
CN107665363A (zh) | 一种总磷的水质软测量预测方法 | |
CN105160422B (zh) | 基于自组织级联神经网络的污水处理出水总磷预测方法 | |
CN113433086B (zh) | 一种模糊神经网络结合分光光度法预测水质cod的方法 | |
Al-Musawi et al. | Prediction and assessment of water quality index using neural network model and GIS case study: Tigris River in Baghdad city | |
Razif | BOD, COD, and TSS Predictions from DO measurement results for the Surabaya River, Indonesia | |
Hore et al. | Application of an artificial neural network in wastewater quality monitoring: prediction of water quality index | |
CN114330951A (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 | ||
RJ01 | Rejection of invention patent application after publication |
Application publication date: 20180206 |