CN107688701A - 基于wasp模型的水质软测量及水体富营养化评价方法 - Google Patents
基于wasp模型的水质软测量及水体富营养化评价方法 Download PDFInfo
- Publication number
- CN107688701A CN107688701A CN201710729131.7A CN201710729131A CN107688701A CN 107688701 A CN107688701 A CN 107688701A CN 201710729131 A CN201710729131 A CN 201710729131A CN 107688701 A CN107688701 A CN 107688701A
- Authority
- CN
- China
- Prior art keywords
- mrow
- msub
- mtd
- mtr
- water quality
- 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 152
- 238000011156 evaluation Methods 0.000 title claims abstract description 84
- 238000012851 eutrophication Methods 0.000 title claims abstract description 67
- 238000005259 measurement Methods 0.000 title claims abstract description 50
- 238000000034 method Methods 0.000 claims description 51
- 238000000691 measurement method Methods 0.000 claims description 23
- 239000011159 matrix material Substances 0.000 claims description 22
- 230000008569 process Effects 0.000 claims description 14
- 238000013441 quality evaluation Methods 0.000 claims description 11
- IJGRMHOSHXDMSA-UHFFFAOYSA-N Atomic nitrogen Chemical compound N#N IJGRMHOSHXDMSA-UHFFFAOYSA-N 0.000 claims description 7
- QVGXLLKOCUKJST-UHFFFAOYSA-N atomic oxygen Chemical compound [O] QVGXLLKOCUKJST-UHFFFAOYSA-N 0.000 claims description 7
- 229910052760 oxygen Inorganic materials 0.000 claims description 7
- 239000001301 oxygen Substances 0.000 claims description 7
- OAICVXFJPJFONN-UHFFFAOYSA-N Phosphorus Chemical compound [P] OAICVXFJPJFONN-UHFFFAOYSA-N 0.000 claims description 6
- 229910052698 phosphorus Inorganic materials 0.000 claims description 6
- 239000011574 phosphorus Substances 0.000 claims description 6
- 230000015572 biosynthetic process Effects 0.000 claims description 5
- 125000001477 organic nitrogen group Chemical group 0.000 claims description 5
- 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 3
- 230000001186 cumulative effect Effects 0.000 claims description 3
- 229910052757 nitrogen Inorganic materials 0.000 claims description 3
- 238000003786 synthesis reaction Methods 0.000 claims description 3
- OKTJSMMVPCPJKN-UHFFFAOYSA-N Carbon Chemical compound [C] OKTJSMMVPCPJKN-UHFFFAOYSA-N 0.000 claims description 2
- XKMRRTOUMJRJIA-UHFFFAOYSA-N ammonia nh3 Chemical compound N.N XKMRRTOUMJRJIA-UHFFFAOYSA-N 0.000 claims description 2
- 229910052799 carbon Inorganic materials 0.000 claims description 2
- 230000008859 change Effects 0.000 claims description 2
- 238000006243 chemical reaction Methods 0.000 claims description 2
- 229930002868 chlorophyll a Natural products 0.000 claims description 2
- ATNHDLDRLWWWCB-AENOIHSZSA-M chlorophyll a Chemical compound C1([C@@H](C(=O)OC)C(=O)C2=C3C)=C2N2C3=CC(C(CC)=C3C)=[N+]4C3=CC3=C(C=C)C(C)=C5N3[Mg-2]42[N+]2=C1[C@@H](CCC(=O)OC\C=C(/C)CCC[C@H](C)CCC[C@H](C)CCCC(C)C)[C@H](C)C2=C5 ATNHDLDRLWWWCB-AENOIHSZSA-M 0.000 claims description 2
- 238000007619 statistical method Methods 0.000 claims description 2
- 238000005516 engineering process Methods 0.000 abstract description 2
- 238000004088 simulation Methods 0.000 description 11
- 230000004927 fusion Effects 0.000 description 7
- 230000007246 mechanism Effects 0.000 description 6
- 238000013473 artificial intelligence Methods 0.000 description 4
- 238000004422 calculation algorithm Methods 0.000 description 4
- 238000004364 calculation method Methods 0.000 description 4
- 238000012544 monitoring process Methods 0.000 description 4
- 235000016709 nutrition Nutrition 0.000 description 3
- 230000035764 nutrition Effects 0.000 description 3
- 239000002028 Biomass Substances 0.000 description 2
- 238000009825 accumulation Methods 0.000 description 2
- 239000005422 algal bloom Substances 0.000 description 2
- 230000001351 cycling effect Effects 0.000 description 2
- 230000005059 dormancy Effects 0.000 description 2
- 229910000069 nitrogen hydride Inorganic materials 0.000 description 2
- 235000003784 poor nutrition Nutrition 0.000 description 2
- 238000011160 research Methods 0.000 description 2
- 241000192710 Microcystis aeruginosa Species 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000009826 distribution Methods 0.000 description 1
- 238000004134 energy conservation Methods 0.000 description 1
- 230000007613 environmental effect Effects 0.000 description 1
- 238000001914 filtration Methods 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 239000000463 material Substances 0.000 description 1
- 238000013178 mathematical model Methods 0.000 description 1
- 238000005457 optimization Methods 0.000 description 1
- 230000001737 promoting effect Effects 0.000 description 1
- 238000011158 quantitative evaluation Methods 0.000 description 1
- 238000011084 recovery Methods 0.000 description 1
- 238000005070 sampling Methods 0.000 description 1
- 239000010865 sewage Substances 0.000 description 1
- 239000000126 substance Substances 0.000 description 1
- 238000012360 testing method Methods 0.000 description 1
- 238000012546 transfer Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
-
- 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
-
- 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
- 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/186—Water using one or more living organisms, e.g. a fish
-
- 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/188—Determining the state of nitrification
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06N—COMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
- G06N7/00—Computing arrangements based on specific mathematical models
- G06N7/02—Computing arrangements based on specific mathematical models using fuzzy logic
-
- 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)
- Life Sciences & Earth Sciences (AREA)
- Health & Medical Sciences (AREA)
- Chemical & Material Sciences (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- General Health & Medical Sciences (AREA)
- Biochemistry (AREA)
- Food Science & Technology (AREA)
- Medicinal Chemistry (AREA)
- Analytical Chemistry (AREA)
- Immunology (AREA)
- Pathology (AREA)
- Theoretical Computer Science (AREA)
- Molecular Biology (AREA)
- Biomedical Technology (AREA)
- Software Systems (AREA)
- General Engineering & Computer Science (AREA)
- Evolutionary Computation (AREA)
- Fuzzy Systems (AREA)
- Algebra (AREA)
- Artificial Intelligence (AREA)
- Computational Mathematics (AREA)
- Data Mining & Analysis (AREA)
- Automation & Control Theory (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Pure & Applied Mathematics (AREA)
- Computing Systems (AREA)
- Emergency Medicine (AREA)
- Mathematical Physics (AREA)
- Biodiversity & Conservation Biology (AREA)
- Computer Hardware Design (AREA)
- Geometry (AREA)
- Other Investigation Or Analysis Of Materials By Electrical Means (AREA)
- Measuring Or Testing Involving Enzymes Or Micro-Organisms (AREA)
Abstract
本发明公开了一种基于WASP模型的水质软测量及水体富营养化评价方法,属于环境工程技术领域。本发明首先建立基于WASP模型的水质软测量模型,同时结合UKF对水质软测量模型的未知参数进行估计,以提高水质软测量模型的准确性。然后利用改进的模糊综合评价法对水体进行富营养化评价,以进一步验证水质软测量模型的有效性。本发明可解决水质指标不能在线测量的问题,在提高对水质软测量模型的准确性的同时,还可以得到实时的评价结果并提高富营养化评价的准确度。
Description
技术领域
本发明涉及一种基于WASP模型的水质软测量方法及水体富营养评价方法,属于环境工程技术领域。具体地说,是在深入研究WASP机理模型的基础上,对模型中的未知参数利用无迹卡尔曼滤波(UKF)进行估计,从而构建融合水质软测量方法,最后利用改进的模糊综合评价法对水体状态进行富营养化评价。
背景技术
随着现代社会工业科技的飞速发展,工业和生活污水日益增多,湖库富营养化现象越来越严重,其中氮和磷是湖库等慢流水体富营养化现象加剧的主要原因。目前,亚太地区中54%的湖泊存在富营养化现象,因此,深入研究富营养化的生成过程,获取各水质指标的实时浓度值,及时地评价或预测水体的富营养化状态,对促进水环境保护和技术进步具有重要意义。
然而,一些水质指标例如生化需氧量(BOD)和总氮(TN)等不能在线测量,便提出了水质软测量模型来解决此问题。目前,软测量建模有三大类基本方法:机理建模、辨识建模和基于人工智能的建模。机理建模是在分析系统机理的基础上,利用基本的物理和化学定律,如物料或能量守恒关系,得到数学关系表达式的方法;辨识建模方法直接利用输入输出数据所提供的信息建立数学模型;基于人工智能的建模是利用人工智能方法对实际系统或系统的某一部分进行描述和表达的过程。上述三类软测量建模方法各有所长,但对于实际工业过程,机理建模可能代价很高,引入的各种假设条件也会影响模型的精度;而辨识建模和基于人工智能的建模方法均存在精度不高等缺点。
发明内容
本发明为了解决现有的水质指标未能实现实时测量,以及对富营养化评价方法存在不确定性的问题,提出一种基于WASP模型的水质软测量及水体富营养化评价方法,为水体富营养化状态进一步预测或评价提供了理论支持。
本发明是一种基于WASP模型的水质软测量及水体富营养化评价方法,其设计思想是:首先建立基于WASP模型的水质软测量模型,同时结合UKF对水质软测量模型的未知参数进行估计,以提高水质软测量模型的准确性。然后利用改进的模糊综合评价法对水体进行富营养化评价,以进一步验证水质软测量模型的有效性。
本发明提供一种基于WASP模型的水质软测量方法及水体富营养化评价技术方法,该方法具体步骤如下:
步骤一:建立基于WASP模型的水质软测量模型。
WASP模型中的富营养化(EUTRO)模块,可描述水质指标的动态变化,包括溶解氧(DO)、生化需氧量(BOD)、浮游植物碳(Phyt)(或叶绿素a(Chl-a))、硝酸盐氮(NO3-N)、氨氮(NH3-N)、有机氮(ON)和有机磷(OP)等。且他们之间的相互关系可用四个反应系统表示,即浮游植物动力学、磷循环、氮循环和溶解氧平衡。水质软测量模型便由上述七个水质指标方程组成的一组非线性连续微分方程组得到。
步骤二:基于UKF的模型参数估计。
在实际中,UKF是状态估计的一种有效方法。UKF摒弃了非线性函数进行线性化的传统做法,采用Kalman线性滤波框架,对于一步预测方程,使用无迹变换(UT)来处理均值和协方差的非线性传递问题。UKF算法是用一系列确定样本来逼近状态的后验概率密度,而不是对非线性函数进行近似,不需要对Jacobian矩阵进行求导,且UKF没有把高阶项忽略,因此对于非线性分布的统计量有较高的计算精度。故将UKF应用于水质软测量模型中未知参数的优化估计,可进一步提高水质软测量模型的准确性。
步骤三:融合的水质软测量方法的构建。
基于WASP模型和基于UKF的模型参数估计,构建融合的水质软测量方法,结合水质指标观测值,便可得到各水质指标的实时浓度值。
步骤四:基于改进模糊综合评价法的水体富营养化评价。
在传统模糊综合评价法的基础上,对评价指标选取及其权重的计算方法进行了改良。采用累积频率法对评价指标进行选取,是通过计算各水质指标的超标倍数百分比的累积频率,选取对水环境影响较大的水质指标作为评价指标;利用聚类权法确定各评价指标的权重,此法将实测值与标准值相结合,更加客观地反映了各评价指标在所有指标中的相对重要性。然后,利用改进的模糊综合评价法对水体富营养化状态进行评价,结合融合水质软测量方法得到的水质指标模拟值与实测值,并将两者的评价结果进行对比,可进一步验证方法的有效性。
本发明的优点在于:
1、本发明将WASP模型和UKF结合起来,构建融合的水质软测量方法,可解决BOD和TN等水质指标不能在线测量的问题。
2、本发明利用UKF对水质软测量模型中的未知参数进行优化估计,可进一步提高对水质软测量模型的准确性。
3、本发明通过对传统的模糊综合评价法中评价指标的选取及其权重的计算方法进行改进,提出改进的模糊综合评价法进行水体富营养化评价,可得到实时的评价结果并提高富营养化评价的准确度。
附图说明
图1:本发明基于WASP模型的水质软测量方法及水体富营养化评价技术方法的流程图。
图2:本发明基于UKF的模型参数实时估计结果图。
图3:本发明基于融合水质软测量方法的水质指标模拟结果与实测值、基于试算法和非线性最小二乘法的模拟值对比图。
图4:本发明基于改进的模糊综合评价法中水质评价指标的三角线性隶属函数图。
图5:本发明基于融合水质软测量方法的水质指标模拟值与实测值的富营养化等级的隶属度评价结果图。
图6:本发明基于融合水质软测量方法的水质指标模拟值与实测值的实时富营养化评价结果图。
具体实施方式
下面将结合附图和实施例对本发明作进一步的详细说明。
本发明是一种基于WASP模型的水质软测量方法及水体富营养化评价方法,具体实施流程如图1所示,通过如下步骤实现:
步骤一:建立基于WASP模型的水质软测量模型。
水质软测量模型是基于WASP模型中EUTRO模块的机理方程构建的,即由DO、BOD、Phyt、NO3-N、NH3-N、ON和OP七个水质指标方程组成的一组非线性连续微分方程,具体方程如下:
其中,分别表示DO、BOD、Phyt、NO3-N、NH3-N、ON和OP的浓度值,公式(1)中未知模型参数的名称如表1所示。
表1未知模型参数的名称或定义
令θ=(k1,k2,k3,k4,k5,k6,k7,k8,k9,k10)Τ,则过程方程可写为:
又令X=[xT,θT]T,然后加入过程噪声w,可得扩展的过程方程为:
其中,w为过程噪声,且满足w~N(0,Q),Q为过程协方差矩阵。同样地,也可得到测量方程为:
其中,y为观测量,h为观测矩阵,v为测量噪声,且满足v~N(0,R),R为测量协方差矩阵。因此,便可得到水质指标的连续动态模型为:
公式(5)即为水质软测量模型。
步骤二:基于UKF的模型参数估计。
为了得到离散时间模型,利用四阶龙格库塔法进行方程的离散化,便可得到公式(5)的离散形式为:
其中,k表示采样时间点,也是离散点,X(k)表示离散点k的各水质指标的浓度值,v(k)为离散点k的测量噪声,w(k-1)为离散点k-1的过程噪声,X(k)=[xT(k),θT]T。
在此基础上,便可利用UKF进行未知模型参数的优化估计,结合一组观测量Y(k)=[y(1),y(2),…,y(k)]。最终,可得到模型参数估计值
步骤三:融合的水质软测量方法的构建。
通过建立基于WASP模型的水质软测量模型,然后利用UKF对水质软测量模型中的未知模型参数进行估计,从而构建融合的水质软测量方法。结合水质指标观测值,便可得到BOD和TN等水质指标的实时浓度值。
步骤四:基于改进模糊综合评价法的水体富营养化评价。
结合融合的水质软测量方法的输出值和水质指标在线测量值,采用改进的模糊综合评价法对水体富营养化状态进行评价。模糊综合评价法基于模糊数学理论,从对象的复杂性和评价指标的模糊性出发,从而获得定量的评价结果。具体步骤如下:
(4.1)水质评价指标的选取。
在传统的模糊综合评价法中,通常采用经验判断来选取对水环境影响较大的水质指标作为评价指标,此法虽相对简单,但缺乏理论依据,因此选用累积频率法来进行评价指标的选取。
评价集是对评价对象进行评价的标准集合,假设水体富营养化状态为n个等级,即:
V={V1,V2,…,Vj,…,Vn} (7)
其中,Vj为第j个富营养化评价等级。进一步假设,存在M个评价指标,即{C1,C2,…,CM}。则累积频率法具体的公式如下:
其中,i为水质评价指标的标签;j为富营养化评价等级的标签;Ci为第i个指标的浓度值;σij为第i个指标的j类标准值;βi为第i个指标的超标倍数;和为计算中间变量;Ki为前i个指标的累积频率。根据统计分析要求,在选取评价指标时,一般取Ki≥85%。在此基础上,便可得到m个关键的水质评价指标,即C={C1,C2,…,Cm}(m∈M)。
(4.2)隶属函数和模糊关系矩阵的建立。
采用模糊数学理论进行评价研究,最关键的问题是建立隶属函数,三角线性隶属函数在实际中使用广泛,便依此来进行模糊关系矩阵的构建。可得模糊关系矩阵R:C→V如下:
其中,rij为第j个富营养化评价等级Vj中第i个水质评价指标Ci的隶属度,且
(4.3)评价指标权重的确定。
评价指标权重的确定是影响最终评价结果的重要因素之一。利用聚类权法进行指标权重的确定,此法将实测值与标准值相结合,更加客观地反映了各评价指标在所有指标中的相对重要性。具体公式如下所示:
其中,wij是第j类富营养化等级的第i个评价指标权重。因此,便可得到指标权重矩阵W,即
W=(W1,W2,…,Wj,…,Wn) (11)
其中,Wj是第j类富营养化等级的指标权重矩阵,即
Wj=(w1j,w2j,…,wij,…,wmj) (12)
(4.4)模糊合成运算。
结合指标权重W和模糊关系矩阵R,基于全部评价指标的模糊综合评价结果B便可得到。采用“加权平均型”的相乘相加法,来进行评价结果的计算,该方法根据权重的大小对所有指标均衡兼顾,保留了单指标评价的所有信息,能够切实地反映出水体富营养化的状况。具体公式如下所示:
其中,bj是对应于第j类富营养化等级的相对隶属度。根据最大隶属度原则,最终的水体富营养化等级可得出,即
实施例1:
实验数据来源于太湖某一监测点的数据,包括DO、BOD、Ch-a、TN、NH3-N和总磷(TP)等水质指标的实测浓度值,每天一组数据,共包括336组数据。应用本发明提供的方法进行水质软测量和水体富营养化评价,具体步骤如下:
步骤一:建立基于WASP模型的水质软测量模型:
根据太湖的水质指标监测点数据,构建基于WASP模型的太湖水质软测量模型。
步骤二、基于UKF的模型参数估计:
由于监测数据的限制,用TN浓度值代替NO3-N数据,结合已有的DO,BOD,Phyt和NH3-N数据,观测矩阵可设置为:
便可得到十个未知模型参数的实时估计值,具体见图2所示。
藻类水华是富营养化的主要特征,其形成阶段可分为复苏期、生物量增加和积聚期、休眠期三个阶段,可将一年时间具体划分为:1至3月份为复苏期,4至10月中旬为生物量增加和积聚阶段,其余月份为休眠期。依据水华形成的多阶段原理,分别计算每个阶段的模型参数平均估计值,具体结果如表2所示。
表2每个阶段的模型参数平均估计值
步骤三:融合水质软测量方法的构建:
基于WASP模型和UKF的模型参数估计,构建融合水质软测量方法,结合水质指标观测值,便可得到各水质指标的实时浓度模拟值,然后分别与实测值、基于试算法的模拟值和基于非线性最小二乘的模拟值进行比较。图3表示DO、BOD、NH3-N、TN和Chl-a的浓度值的模拟值,其中Chl-a的值可由Phyt值换算得出。
为了定量地表示模型误差,采用均方根误差(RMSE)值计算水质指标模拟值与实测值之间的偏差,具体的计算结果如表3所示。
表3不同方法获得的水质指标RMSE值
由图3和表3均可得出,融合水质软测量方法得出的水质指标模拟值与实测值吻合较好,可验证此方法的有效性。且基于UKF的模型参数估计较试算法与非线性最小二乘法,可进一步提高水质软测量方法的准确性。
步骤四:基于改进模糊综合评价法的水体富营养化评价:
利用改进的模糊综合评价法,分别对基于太湖的水质指标模拟值与实测值进行水体富营养化评价,且将两者的评价结果进行对比。
(1)水质评价指标的选取。
通过查阅相应的文献,湖库水体富营养化状态可分为五个等级:I(极贫营养),II(贫营养),III(中营养),IV(富营养)和V(重富营养),即V={I,II,III,IV,V}。在此基础上,根据式(8)具体选取出水质评价指标,包括透明度(SD),BOD,TN,TP,DO和Chl-a,其中BOD,TN,DO和Chl-a是融合水质软测量方法的在线测量值,SD和TP值是取自太湖实测数据。表4为水质指标各富营养化等级标准值。
表4水质指标各富营养化等级标准值
等级 | SD(m) | BOD(mg/L) | TN(mg/L) | TP(mg/L) | DO(mg/L) | Chl-a(mg/L) |
I | 10.00 | 0.18 | 0.02 | 0.001 | 0.015 | 0.001 |
II | 5.00 | 0.24 | 0.31 | 0.004 | 0.0085 | 0.002 |
III | 1.50 | 1.20 | 1.20 | 0.023 | 0.005 | 0.004 |
IV | 0.55 | 6.00 | 3.60 | 0.110 | 0.002 | 0.01 |
V | 0.17 | 15.00 | 4.70 | 0.660 | 0.001 | 0.065 |
(2)隶属函数和模糊关系矩阵的建立。
水质指标的三角线性隶属函数如图4,根据式(9)和(10),结合实时的观测值,便可得到与时间相关的模糊关系矩阵和指标权重矩阵,即R(t)和W(t)。
(3)模糊合成运算。
根据式(12),与时间相关的模糊综合评价结果B(t)可求出,图5所示为基于实测值与模拟值的各富营养化等级隶属度评价结果。然后依据最大隶属度原则,基于两者的水体富营养化评价结果可得出,结果曲线如图6所示。
由图6可得出,太湖某监测点的水体富营养化状态处于富营养等级,且基于实测值与模拟值的结果基本吻合,以验证改进的模糊综合评价方法的有效性。
Claims (3)
1.基于WASP模型的水质软测量及水体富营养化评价方法,其特征在于:包括如下步骤,
步骤一:建立基于WASP模型的水质软测量模型;
采用WASP模型中的富营养化模块描述水质指标的动态变化,所述的水质指标包括溶解氧、生化需氧量、浮游植物碳或叶绿素a、硝酸盐氮、氨氮、有机氮和有机磷,各水质指标之间的相互关系用四个反应系统表示,即浮游植物动力学、磷循环、氮循环和溶解氧平衡,根据上述水质指标的方程组成的一组非线性连续微分方程组,便得到水质软测量模型;
步骤二:基于UKF的模型参数估计;
步骤三:融合的水质软测量方法的构建;
基于WASP模型和基于UKF的模型参数估计,构建融合的水质软测量方法,结合水质指标观测值,得到各水质指标的实时浓度值;
步骤四:基于改进模糊综合评价法的水体富营养化评价;
结合融合的水质软测量方法的输出值和水质指标在线测量值,采用改进的模糊综合评价法对水体富营养化状态进行评价。
2.根据权利要求1所述的基于WASP模型的水质软测量及水体富营养化评价方法,其特征在于:
水质软测量模型的非线性连续微分方程组如下:
<mrow>
<msub>
<mover>
<mi>C</mi>
<mo>&CenterDot;</mo>
</mover>
<mrow>
<mi>D</mi>
<mi>O</mi>
</mrow>
</msub>
<mo>=</mo>
<msub>
<mi>k</mi>
<mn>1</mn>
</msub>
<mrow>
<mo>(</mo>
<msub>
<mi>C</mi>
<mrow>
<msub>
<mi>NO</mi>
<mn>3</mn>
</msub>
<mo>-</mo>
<mi>N</mi>
</mrow>
</msub>
<mo>-</mo>
<msub>
<mi>C</mi>
<mrow>
<mi>D</mi>
<mi>O</mi>
</mrow>
</msub>
<mo>)</mo>
</mrow>
<mo>-</mo>
<mn>1.047</mn>
<msub>
<mi>k</mi>
<mn>2</mn>
</msub>
<mrow>
<mo>(</mo>
<mfrac>
<msub>
<mi>C</mi>
<mrow>
<mi>D</mi>
<mi>O</mi>
</mrow>
</msub>
<mrow>
<mn>0.5</mn>
<mo>+</mo>
<msub>
<mi>C</mi>
<mrow>
<mi>D</mi>
<mi>O</mi>
</mrow>
</msub>
</mrow>
</mfrac>
<mo>)</mo>
</mrow>
<msub>
<mi>C</mi>
<mrow>
<mi>B</mi>
<mi>O</mi>
<mi>D</mi>
</mrow>
</msub>
<mo>-</mo>
<mn>4.94</mn>
<msub>
<mi>k</mi>
<mn>3</mn>
</msub>
<mo>-</mo>
<mn>5.4</mn>
<msub>
<mi>k</mi>
<mn>4</mn>
</msub>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>14.53</mn>
<msub>
<mi>k</mi>
<mn>5</mn>
</msub>
<mo>+</mo>
<mn>0.35</mn>
<mo>)</mo>
</mrow>
<msub>
<mi>C</mi>
<mrow>
<mi>P</mi>
<mi>h</mi>
<mi>y</mi>
<mi>t</mi>
</mrow>
</msub>
</mrow>
<mrow>
<msub>
<mover>
<mi>C</mi>
<mo>&CenterDot;</mo>
</mover>
<mrow>
<mi>B</mi>
<mi>O</mi>
<mi>D</mi>
</mrow>
</msub>
<mo>=</mo>
<mn>2.67</mn>
<msub>
<mi>k</mi>
<mn>6</mn>
</msub>
<msub>
<mi>C</mi>
<mrow>
<mi>P</mi>
<mi>h</mi>
<mi>y</mi>
<mi>t</mi>
</mrow>
</msub>
<mo>-</mo>
<mn>1.047</mn>
<msub>
<mi>k</mi>
<mn>2</mn>
</msub>
<mrow>
<mo>(</mo>
<mfrac>
<msub>
<mi>C</mi>
<mrow>
<mi>D</mi>
<mi>O</mi>
</mrow>
</msub>
<mrow>
<mn>0.5</mn>
<mo>+</mo>
<msub>
<mi>C</mi>
<mrow>
<mi>D</mi>
<mi>O</mi>
</mrow>
</msub>
</mrow>
</mfrac>
<mo>)</mo>
</mrow>
<msub>
<mi>C</mi>
<mrow>
<mi>B</mi>
<mi>O</mi>
<mi>D</mi>
</mrow>
</msub>
<mo>-</mo>
<mn>2.99</mn>
<msub>
<mi>k</mi>
<mn>7</mn>
</msub>
<mrow>
<mo>(</mo>
<mfrac>
<mn>0.1</mn>
<mrow>
<mn>0.1</mn>
<mo>+</mo>
<msub>
<mi>C</mi>
<mrow>
<mi>D</mi>
<mi>O</mi>
</mrow>
</msub>
</mrow>
</mfrac>
<mo>)</mo>
</mrow>
<msub>
<mi>C</mi>
<mrow>
<msub>
<mi>NO</mi>
<mn>3</mn>
</msub>
<mo>-</mo>
<mi>N</mi>
</mrow>
</msub>
</mrow>
<mrow>
<msub>
<mover>
<mi>C</mi>
<mo>&CenterDot;</mo>
</mover>
<mrow>
<mi>P</mi>
<mi>h</mi>
<mi>y</mi>
<mi>t</mi>
</mrow>
</msub>
<mo>=</mo>
<mo>&lsqb;</mo>
<msub>
<mi>k</mi>
<mn>5</mn>
</msub>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>0.125</mn>
<mo>+</mo>
<msub>
<mi>k</mi>
<mn>6</mn>
</msub>
<mo>)</mo>
</mrow>
<mo>-</mo>
<mn>0.02</mn>
<mo>&rsqb;</mo>
<msub>
<mi>C</mi>
<mrow>
<mi>P</mi>
<mi>h</mi>
<mi>y</mi>
<mi>t</mi>
</mrow>
</msub>
</mrow>
<mrow>
<msub>
<mover>
<mi>C</mi>
<mo>&CenterDot;</mo>
</mover>
<mrow>
<msub>
<mi>NO</mi>
<mn>3</mn>
</msub>
<mo>-</mo>
<mi>N</mi>
</mrow>
</msub>
<mo>=</mo>
<mn>1.08</mn>
<msub>
<mi>k</mi>
<mn>3</mn>
</msub>
<mrow>
<mo>(</mo>
<mfrac>
<msub>
<mi>C</mi>
<mrow>
<mi>D</mi>
<mi>O</mi>
</mrow>
</msub>
<mrow>
<msub>
<mi>k</mi>
<mn>8</mn>
</msub>
<mo>+</mo>
<msub>
<mi>C</mi>
<mrow>
<mi>D</mi>
<mi>O</mi>
</mrow>
</msub>
</mrow>
</mfrac>
<mo>)</mo>
</mrow>
<msub>
<mi>C</mi>
<mrow>
<msub>
<mi>NH</mi>
<mn>3</mn>
</msub>
<mo>-</mo>
<mi>N</mi>
</mrow>
</msub>
<mo>+</mo>
<mn>0.645</mn>
<msub>
<mi>k</mi>
<mn>5</mn>
</msub>
<msub>
<mi>C</mi>
<mrow>
<mi>P</mi>
<mi>h</mi>
<mi>y</mi>
<mi>t</mi>
</mrow>
</msub>
<mo>-</mo>
<mn>1.045</mn>
<msub>
<mi>k</mi>
<mn>7</mn>
</msub>
<mrow>
<mo>(</mo>
<mfrac>
<mn>0.1</mn>
<mrow>
<mn>0.1</mn>
<mo>+</mo>
<msub>
<mi>C</mi>
<mrow>
<msub>
<mi>NH</mi>
<mn>3</mn>
</msub>
<mo>-</mo>
<mi>N</mi>
</mrow>
</msub>
</mrow>
</mfrac>
<mo>)</mo>
</mrow>
<msub>
<mi>C</mi>
<mrow>
<mi>B</mi>
<mi>O</mi>
<mi>D</mi>
</mrow>
</msub>
</mrow>
<mrow>
<msub>
<mover>
<mi>C</mi>
<mo>&CenterDot;</mo>
</mover>
<mrow>
<msub>
<mi>NH</mi>
<mn>3</mn>
</msub>
<mo>-</mo>
<mi>N</mi>
</mrow>
</msub>
<mo>=</mo>
<mn>0.075</mn>
<mrow>
<mo>(</mo>
<mn>0.125</mn>
<mo>+</mo>
<msub>
<mi>k</mi>
<mn>6</mn>
</msub>
<mo>)</mo>
</mrow>
<msub>
<mi>C</mi>
<mrow>
<mi>P</mi>
<mi>h</mi>
<mi>y</mi>
<mi>t</mi>
</mrow>
</msub>
<mo>+</mo>
<mn>1.08</mn>
<msub>
<mi>k</mi>
<mn>9</mn>
</msub>
<mrow>
<mo>(</mo>
<mfrac>
<msub>
<mi>C</mi>
<mrow>
<mi>P</mi>
<mi>h</mi>
<mi>y</mi>
<mi>t</mi>
</mrow>
</msub>
<mrow>
<mn>1</mn>
<mo>+</mo>
<msub>
<mi>C</mi>
<mrow>
<mi>P</mi>
<mi>h</mi>
<mi>y</mi>
<mi>t</mi>
</mrow>
</msub>
</mrow>
</mfrac>
<mo>)</mo>
</mrow>
<msub>
<mi>C</mi>
<mrow>
<mi>O</mi>
<mi>N</mi>
</mrow>
</msub>
<mo>+</mo>
<mn>0.645</mn>
<msub>
<mi>k</mi>
<mn>5</mn>
</msub>
<msub>
<mi>C</mi>
<mrow>
<mi>P</mi>
<mi>h</mi>
<mi>y</mi>
<mi>t</mi>
</mrow>
</msub>
<mo>-</mo>
<mn>1.08</mn>
<msub>
<mi>k</mi>
<mn>3</mn>
</msub>
<mrow>
<mo>(</mo>
<mfrac>
<msub>
<mi>C</mi>
<mrow>
<mi>D</mi>
<mi>O</mi>
</mrow>
</msub>
<mrow>
<msub>
<mi>k</mi>
<mn>8</mn>
</msub>
<mo>+</mo>
<msub>
<mi>C</mi>
<mrow>
<mi>D</mi>
<mi>O</mi>
</mrow>
</msub>
</mrow>
</mfrac>
<mo>)</mo>
</mrow>
<msub>
<mi>C</mi>
<mrow>
<msub>
<mi>NH</mi>
<mn>3</mn>
</msub>
<mo>-</mo>
<mi>N</mi>
</mrow>
</msub>
</mrow>
<mrow>
<msub>
<mover>
<mi>C</mi>
<mo>&CenterDot;</mo>
</mover>
<mrow>
<mi>O</mi>
<mi>N</mi>
</mrow>
</msub>
<mo>=</mo>
<mn>0.075</mn>
<mrow>
<mo>(</mo>
<mn>0.125</mn>
<mo>+</mo>
<msub>
<mi>k</mi>
<mn>6</mn>
</msub>
<mo>)</mo>
</mrow>
<msub>
<mi>C</mi>
<mrow>
<mi>P</mi>
<mi>h</mi>
<mi>y</mi>
<mi>t</mi>
</mrow>
</msub>
<mo>+</mo>
<mn>1.08</mn>
<msub>
<mi>k</mi>
<mn>9</mn>
</msub>
<mrow>
<mo>(</mo>
<mfrac>
<msub>
<mi>C</mi>
<mrow>
<mi>P</mi>
<mi>h</mi>
<mi>y</mi>
<mi>t</mi>
</mrow>
</msub>
<mrow>
<mn>1</mn>
<mo>+</mo>
<msub>
<mi>C</mi>
<mrow>
<mi>P</mi>
<mi>h</mi>
<mi>y</mi>
<mi>t</mi>
</mrow>
</msub>
</mrow>
</mfrac>
<mo>)</mo>
</mrow>
<msub>
<mi>C</mi>
<mrow>
<mi>O</mi>
<mi>P</mi>
</mrow>
</msub>
</mrow>
其中,分别表示DO、BOD、Phyt、NO3-N、NH3-N、ON和OP的浓度值,公式(1)中未知模型参数的名称或定义如下表所示:
令θ=(k1,k2,k3,k4,k5,k6,k7,k8,k9,k10)Τ,则过程方程写为:
<mrow>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<mover>
<mi>x</mi>
<mo>&CenterDot;</mo>
</mover>
</mtd>
</mtr>
<mtr>
<mtd>
<mover>
<mi>&theta;</mi>
<mo>&CenterDot;</mo>
</mover>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>=</mo>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<mrow>
<mi>f</mi>
<mrow>
<mo>(</mo>
<mi>x</mi>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<msub>
<mn>0</mn>
<mrow>
<mn>10</mn>
<mo>&times;</mo>
<mn>1</mn>
</mrow>
</msub>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>2</mn>
<mo>)</mo>
</mrow>
</mrow>
又令X=[xT,θT]T,然后加入过程噪声w,得扩展的过程方程为:
<mrow>
<mover>
<mi>X</mi>
<mo>&CenterDot;</mo>
</mover>
<mo>=</mo>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<mrow>
<mi>f</mi>
<mrow>
<mo>(</mo>
<mi>x</mi>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<msub>
<mn>0</mn>
<mrow>
<mn>10</mn>
<mo>&times;</mo>
<mn>1</mn>
</mrow>
</msub>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>+</mo>
<mi>w</mi>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>3</mn>
<mo>)</mo>
</mrow>
</mrow>
其中,w为过程噪声,且满足w~N(0,Q),Q为过程协方差矩阵;同样地,得到测量方程为:
<mrow>
<mi>y</mi>
<mo>=</mo>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<mi>h</mi>
</mtd>
<mtd>
<mn>0</mn>
</mtd>
</mtr>
</mtable>
</mfenced>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<mi>x</mi>
</mtd>
</mtr>
<mtr>
<mtd>
<mi>&theta;</mi>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>+</mo>
<mi>v</mi>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>4</mn>
<mo>)</mo>
</mrow>
</mrow>
其中,y为观测量,h为观测矩阵,v为测量噪声,且满足v~N(0,R),R为测量协方差矩阵;得到水质指标的连续动态模型为:
<mrow>
<mfenced open = "{" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<mover>
<mi>X</mi>
<mo>&CenterDot;</mo>
</mover>
<mo>=</mo>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<mrow>
<mi>f</mi>
<mrow>
<mo>(</mo>
<mi>x</mi>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<msub>
<mn>0</mn>
<mrow>
<mn>10</mn>
<mo>&times;</mo>
<mn>1</mn>
</mrow>
</msub>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>+</mo>
<mi>w</mi>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mi>y</mi>
<mo>=</mo>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<mi>h</mi>
</mtd>
<mtd>
<mn>0</mn>
</mtd>
</mtr>
</mtable>
</mfenced>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<mi>x</mi>
</mtd>
</mtr>
<mtr>
<mtd>
<mi>&theta;</mi>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>+</mo>
<mi>v</mi>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>5</mn>
<mo>)</mo>
</mrow>
</mrow>
公式(5)即为水质软测量模型。
3.根据权利要求1所述的基于WASP模型的水质软测量及水体富营养化评价方法,其特征在于:
步骤四具体步骤如下:
(4.1)水质评价指标的选取;
假设水体富营养化状态为n个等级,即:
V={V1,V2,…,Vj,…,Vn} (7)
其中,Vj为第j个富营养化评价等级,假设存在M个评价指标,即{C1,C2,…,CM};则累积频率法具体的公式如下:
<mrow>
<mfenced open = "{" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<msub>
<mi>&beta;</mi>
<mi>i</mi>
</msub>
<mo>=</mo>
<mfrac>
<msub>
<mi>C</mi>
<mi>i</mi>
</msub>
<mrow>
<mfrac>
<mn>1</mn>
<mi>n</mi>
</mfrac>
<munderover>
<mi>&Sigma;</mi>
<mrow>
<mi>j</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>n</mi>
</munderover>
<msub>
<mi>&sigma;</mi>
<mrow>
<mi>i</mi>
<mi>j</mi>
</mrow>
</msub>
</mrow>
</mfrac>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mover>
<mi>&beta;</mi>
<mo>~</mo>
</mover>
<mi>i</mi>
</msub>
<mo>=</mo>
<mrow>
<mo>{</mo>
<mrow>
<msub>
<mi>&beta;</mi>
<mi>i</mi>
</msub>
<mo>|</mo>
<msub>
<mi>&beta;</mi>
<mn>1</mn>
</msub>
<mo>&GreaterEqual;</mo>
<msub>
<mi>&beta;</mi>
<mi>i</mi>
</msub>
<mo>&GreaterEqual;</mo>
<msub>
<mi>&beta;</mi>
<mi>M</mi>
</msub>
</mrow>
<mo>}</mo>
</mrow>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>K</mi>
<mi>i</mi>
</msub>
<mo>=</mo>
<mfrac>
<mrow>
<munderover>
<mi>&Sigma;</mi>
<mrow>
<mi>j</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>i</mi>
</munderover>
<msub>
<mover>
<mi>&beta;</mi>
<mo>~</mo>
</mover>
<mi>j</mi>
</msub>
</mrow>
<mrow>
<munderover>
<mi>&Sigma;</mi>
<mrow>
<mi>j</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>M</mi>
</munderover>
<msub>
<mover>
<mi>&beta;</mi>
<mo>~</mo>
</mover>
<mi>j</mi>
</msub>
</mrow>
</mfrac>
<mo>&times;</mo>
<mn>100</mn>
<mi>%</mi>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>,</mo>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
<mo>,</mo>
<mn>2</mn>
<mo>,</mo>
<mn>3</mn>
<mo>,</mo>
<mn>...</mn>
<mo>,</mo>
<mi>M</mi>
<mo>;</mo>
<mi>j</mi>
<mo>=</mo>
<mn>1</mn>
<mo>,</mo>
<mn>2</mn>
<mo>,</mo>
<mn>3</mn>
<mo>,</mo>
<mn>...</mn>
<mo>,</mo>
<mi>n</mi>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>8</mn>
<mo>)</mo>
</mrow>
</mrow>
其中,i为水质评价指标的标签;j为富营养化评价等级的标签;Ci为第i个指标的浓度值;σij为第i个指标的j类标准值;βi为第i个指标的超标倍数;和为计算中间变量;Ki为前i个指标的累积频率;
根据统计分析要求,在选取评价指标时,取Ki≥85%,在此基础上,便得到m个关键的水质评价指标,即C={C1,C2,…,Cm},m∈M;
(4.2)隶属函数和模糊关系矩阵的建立;
构建模糊关系矩阵R:C→V如下:
<mrow>
<mi>R</mi>
<mo>=</mo>
<msub>
<mrow>
<mo>(</mo>
<msub>
<mi>r</mi>
<mrow>
<mi>i</mi>
<mi>j</mi>
</mrow>
</msub>
<mo>)</mo>
</mrow>
<mrow>
<mi>m</mi>
<mo>&times;</mo>
<mi>n</mi>
</mrow>
</msub>
<mo>=</mo>
<mfenced open = "(" close = ")">
<mtable>
<mtr>
<mtd>
<msub>
<mi>r</mi>
<mn>11</mn>
</msub>
</mtd>
<mtd>
<msub>
<mi>r</mi>
<mn>12</mn>
</msub>
</mtd>
<mtd>
<mo>...</mo>
</mtd>
<mtd>
<msub>
<mi>r</mi>
<mrow>
<mn>1</mn>
<mi>n</mi>
</mrow>
</msub>
</mtd>
</mtr>
<mtr>
<mtd>
<msub>
<mi>r</mi>
<mn>21</mn>
</msub>
</mtd>
<mtd>
<msub>
<mi>r</mi>
<mn>22</mn>
</msub>
</mtd>
<mtd>
<mo>...</mo>
</mtd>
<mtd>
<msub>
<mi>r</mi>
<mrow>
<mn>2</mn>
<mi>n</mi>
</mrow>
</msub>
</mtd>
</mtr>
<mtr>
<mtd>
<mo>...</mo>
</mtd>
<mtd>
<mo>...</mo>
</mtd>
<mtd>
<mo>...</mo>
</mtd>
<mtd>
<mo>...</mo>
</mtd>
</mtr>
<mtr>
<mtd>
<msub>
<mi>r</mi>
<mrow>
<mi>m</mi>
<mn>1</mn>
</mrow>
</msub>
</mtd>
<mtd>
<msub>
<mi>r</mi>
<mrow>
<mi>m</mi>
<mn>2</mn>
</mrow>
</msub>
</mtd>
<mtd>
<mo>...</mo>
</mtd>
<mtd>
<msub>
<mi>r</mi>
<mrow>
<mi>m</mi>
<mi>n</mi>
</mrow>
</msub>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>9</mn>
<mo>)</mo>
</mrow>
</mrow>
其中,rij为第j个富营养化评价等级Vj中第i个水质评价指标Ci的隶属度,且
(4.3)评价指标权重的确定;
利用聚类权法进行指标权重的确定,具体公式如下所示:
<mrow>
<mtable>
<mtr>
<mtd>
<mrow>
<msub>
<mi>w</mi>
<mrow>
<mi>i</mi>
<mi>j</mi>
</mrow>
</msub>
<mo>=</mo>
<mfrac>
<mrow>
<msub>
<mi>C</mi>
<mi>i</mi>
</msub>
<mo>/</mo>
<msub>
<mi>&sigma;</mi>
<mrow>
<mi>i</mi>
<mi>j</mi>
</mrow>
</msub>
</mrow>
<mrow>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>m</mi>
</munderover>
<msub>
<mi>C</mi>
<mi>i</mi>
</msub>
<mo>/</mo>
<msub>
<mi>&sigma;</mi>
<mrow>
<mi>i</mi>
<mi>j</mi>
</mrow>
</msub>
</mrow>
</mfrac>
</mrow>
</mtd>
<mtd>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
<mo>,</mo>
<mn>2</mn>
<mo>,</mo>
<mo>...</mo>
<mo>,</mo>
<mi>m</mi>
<mo>;</mo>
<mi>j</mi>
<mo>=</mo>
<mn>1</mn>
<mo>,</mo>
<mn>2</mn>
<mo>,</mo>
<mo>...</mo>
<mo>,</mo>
<mi>n</mi>
</mrow>
</mtd>
</mtr>
</mtable>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>10</mn>
<mo>)</mo>
</mrow>
</mrow>
其中,wij是第j类富营养化等级的第i个评价指标权重,指标权重矩阵W为,
W=(W1,W2,…,Wj,…,Wn) (11)
其中,Wj是第j类富营养化等级的指标权重矩阵,即
Wj=(w1j,w2j,…,wij,…,wmj) (12)
(4.4)模糊合成运算;
结合指标权重W和模糊关系矩阵R,基于全部评价指标的模糊综合评价结果B如下所示:
其中,bj是对应于第j类富营养化等级的相对隶属度;根据最大隶属度原则,最终的水体富营养化等级为,
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710729131.7A CN107688701B (zh) | 2017-08-23 | 2017-08-23 | 基于wasp模型的水质软测量及水体富营养化评价方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710729131.7A CN107688701B (zh) | 2017-08-23 | 2017-08-23 | 基于wasp模型的水质软测量及水体富营养化评价方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN107688701A true CN107688701A (zh) | 2018-02-13 |
CN107688701B CN107688701B (zh) | 2020-10-20 |
Family
ID=61152518
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710729131.7A Active CN107688701B (zh) | 2017-08-23 | 2017-08-23 | 基于wasp模型的水质软测量及水体富营养化评价方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107688701B (zh) |
Cited By (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108615090A (zh) * | 2018-03-29 | 2018-10-02 | 北京工商大学 | 基于蒙特卡洛仿真的河湖水质预测及富营养化风险评估方法 |
CN109033648A (zh) * | 2018-08-01 | 2018-12-18 | 北京工商大学 | 基于果蝇优化算法的水质机理建模及水质预测方法 |
CN109214699A (zh) * | 2018-09-19 | 2019-01-15 | 朱明华 | 一种水体评价方法及装置 |
CN109272231A (zh) * | 2018-09-19 | 2019-01-25 | 朱明华 | 一种污水评估数据处理方法及装置 |
CN110766282A (zh) * | 2019-09-20 | 2020-02-07 | 天津大学 | 一种湿地净化能力评估与提升方法 |
CN110889201A (zh) * | 2019-11-05 | 2020-03-17 | 北京工商大学 | 基于混杂动态系统的水质机理建模及健康风险评估方法 |
CN111398548A (zh) * | 2020-04-09 | 2020-07-10 | 中国水利水电科学研究院 | 一种基于植物作用的表面流湿地水体氮磷浓度预测方法 |
CN112499766A (zh) * | 2020-10-13 | 2021-03-16 | 北京林业大学 | 一种快速鉴定植物净化富营养化水体能力的方法 |
CN114295707A (zh) * | 2021-12-28 | 2022-04-08 | 南京大学 | 一种基于机器学习的污水有机氮生物有效性评价方法 |
CN117009887A (zh) * | 2023-07-28 | 2023-11-07 | 中科景感(厦门)生态环境科技有限公司 | 流域水环境质量精细化估算及分析的方法及系统 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102034214A (zh) * | 2010-12-20 | 2011-04-27 | 上海勘测设计研究院 | 一种富营养化湖泊的健康评价方法 |
US20140012556A1 (en) * | 2012-07-06 | 2014-01-09 | Electronics And Telecommunications Research Institute | Apparatus and method for real time water quality prediction using hydrodynamic model |
CN106485053A (zh) * | 2016-09-16 | 2017-03-08 | 天津大学 | 一种基于模糊标识指数的水功能区水质评价方法 |
-
2017
- 2017-08-23 CN CN201710729131.7A patent/CN107688701B/zh active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102034214A (zh) * | 2010-12-20 | 2011-04-27 | 上海勘测设计研究院 | 一种富营养化湖泊的健康评价方法 |
US20140012556A1 (en) * | 2012-07-06 | 2014-01-09 | Electronics And Telecommunications Research Institute | Apparatus and method for real time water quality prediction using hydrodynamic model |
CN106485053A (zh) * | 2016-09-16 | 2017-03-08 | 天津大学 | 一种基于模糊标识指数的水功能区水质评价方法 |
Non-Patent Citations (2)
Title |
---|
S. Y. CHANG等: "Using Unscented Kalman Filter in Subsurface Contaminant Transport Models", 《JOURNAL OF ENVIRONMENTAL INFORMATICS》 * |
SUBRINA TAHSIN等: "Fast eutrophication assessment for stormwater wet detention", 《ENVIRONMENTAL MONITORING AND ASSESSMENT》 * |
Cited By (14)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108615090B (zh) * | 2018-03-29 | 2020-07-17 | 北京工商大学 | 基于蒙特卡洛仿真的河湖水质预测及富营养化风险评估方法 |
CN108615090A (zh) * | 2018-03-29 | 2018-10-02 | 北京工商大学 | 基于蒙特卡洛仿真的河湖水质预测及富营养化风险评估方法 |
CN109033648A (zh) * | 2018-08-01 | 2018-12-18 | 北京工商大学 | 基于果蝇优化算法的水质机理建模及水质预测方法 |
CN109214699A (zh) * | 2018-09-19 | 2019-01-15 | 朱明华 | 一种水体评价方法及装置 |
CN109272231A (zh) * | 2018-09-19 | 2019-01-25 | 朱明华 | 一种污水评估数据处理方法及装置 |
CN110766282A (zh) * | 2019-09-20 | 2020-02-07 | 天津大学 | 一种湿地净化能力评估与提升方法 |
CN110889201A (zh) * | 2019-11-05 | 2020-03-17 | 北京工商大学 | 基于混杂动态系统的水质机理建模及健康风险评估方法 |
CN111398548A (zh) * | 2020-04-09 | 2020-07-10 | 中国水利水电科学研究院 | 一种基于植物作用的表面流湿地水体氮磷浓度预测方法 |
CN111398548B (zh) * | 2020-04-09 | 2021-05-28 | 中国水利水电科学研究院 | 一种基于植物作用的表面流湿地水体氮磷浓度预测方法 |
CN112499766A (zh) * | 2020-10-13 | 2021-03-16 | 北京林业大学 | 一种快速鉴定植物净化富营养化水体能力的方法 |
CN114295707A (zh) * | 2021-12-28 | 2022-04-08 | 南京大学 | 一种基于机器学习的污水有机氮生物有效性评价方法 |
CN114295707B (zh) * | 2021-12-28 | 2024-02-09 | 南京大学 | 一种基于机器学习的污水有机氮生物有效性评价方法 |
CN117009887A (zh) * | 2023-07-28 | 2023-11-07 | 中科景感(厦门)生态环境科技有限公司 | 流域水环境质量精细化估算及分析的方法及系统 |
CN117009887B (zh) * | 2023-07-28 | 2024-02-13 | 中科景感(厦门)生态环境科技有限公司 | 流域水环境质量精细化估算及分析的方法及系统 |
Also Published As
Publication number | Publication date |
---|---|
CN107688701B (zh) | 2020-10-20 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107688701B (zh) | 基于wasp模型的水质软测量及水体富营养化评价方法 | |
CN102854296B (zh) | 一种基于集成神经网络的污水处理软测量方法 | |
Newhart et al. | Data-driven performance analyses of wastewater treatment plants: A review | |
Bagherzadeh et al. | Prediction of energy consumption and evaluation of affecting factors in a full-scale WWTP using a machine learning approach | |
Wei et al. | Comprehensive evaluation model for water environment carrying capacity based on VPOSRM framework: A case study in Wuhan, China | |
CN101630376B (zh) | 多模型神经网络的生物发酵过程软测量建模方法及软仪表 | |
CN101387632B (zh) | 一种污水处理过程中生化需氧量bod的软测量方法 | |
CN101699477B (zh) | 精确预测大坝变形的神经网络方法 | |
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 | |
CN105510546A (zh) | 一种基于自组织递归rbf神经网络的生化需氧量bod智能检测方法 | |
CN111260117B (zh) | 一种基于气象因子的ca-narx水质预测方法 | |
Sajib et al. | Developing a novel tool for assessing the groundwater incorporating water quality index and machine learning approach | |
CN104462850A (zh) | 基于模糊高斯混合模型的多阶段间歇过程软测量方法 | |
CN105574326A (zh) | 一种基于自组织模糊神经网络的出水氨氮浓度软测量方法 | |
CN101566612A (zh) | 一种污水化学需氧量软测量方法 | |
CN103632032A (zh) | 一种城市污水处理过程出水指标在线软测量预测方法 | |
Wang et al. | A full-view management method based on artificial neural networks for energy and material-savings in wastewater treatment plants | |
CN113343601A (zh) | 一种复杂水系湖泊水位和污染物迁移的动态模拟方法 | |
Singh et al. | Review on Data Mining Techniques for Prediction of Water Quality. | |
CN109408896B (zh) | 一种污水厌氧处理产气量多元智能实时监控方法 | |
CN111204867B (zh) | 膜生物反应器-mbr膜污染智能决策方法 | |
CN103675005B (zh) | 最优模糊网络的工业熔融指数软测量仪表及方法 | |
CN113191689B (zh) | 一种耦合主成分分析与bp神经网络的土地适宜评价方法 | |
Putro et al. | Modeling methane emission of wastewater Anaerobic pond at Palm oil mill using radial basis function neural network | |
Liu et al. | ANN-based prediction of ammonia nitrogen for wastewater discharge indicators under carbon neutral trend |
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 |