CN113065242A - 一种基于kplsr模型的污水处理出水总氮浓度软测量方法 - Google Patents
一种基于kplsr模型的污水处理出水总氮浓度软测量方法 Download PDFInfo
- Publication number
- CN113065242A CN113065242A CN202110300719.7A CN202110300719A CN113065242A CN 113065242 A CN113065242 A CN 113065242A CN 202110300719 A CN202110300719 A CN 202110300719A CN 113065242 A CN113065242 A CN 113065242A
- Authority
- CN
- China
- Prior art keywords
- vector
- effluent
- kplsr
- formula
- column
- 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
- IJGRMHOSHXDMSA-UHFFFAOYSA-N Atomic nitrogen Chemical compound N#N IJGRMHOSHXDMSA-UHFFFAOYSA-N 0.000 title claims abstract description 95
- 239000010865 sewage Substances 0.000 title claims abstract description 72
- 229910052757 nitrogen Inorganic materials 0.000 title claims abstract description 48
- 238000000691 measurement method Methods 0.000 title claims abstract description 11
- 238000005259 measurement Methods 0.000 claims abstract description 34
- 238000000034 method Methods 0.000 claims abstract description 30
- 238000005457 optimization Methods 0.000 claims abstract description 5
- 239000013598 vector Substances 0.000 claims description 124
- 239000011159 matrix material Substances 0.000 claims description 36
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims description 30
- 239000007787 solid Substances 0.000 claims description 17
- VEXZGXHMUGYJMC-UHFFFAOYSA-M Chloride anion Chemical compound [Cl-] VEXZGXHMUGYJMC-UHFFFAOYSA-M 0.000 claims description 8
- 238000005273 aeration Methods 0.000 claims description 6
- 238000004364 calculation method Methods 0.000 claims description 5
- 125000004432 carbon atom Chemical group C* 0.000 claims description 5
- QVGXLLKOCUKJST-UHFFFAOYSA-N atomic oxygen Chemical compound [O] QVGXLLKOCUKJST-UHFFFAOYSA-N 0.000 claims description 4
- 239000001301 oxygen Substances 0.000 claims description 4
- 229910052760 oxygen Inorganic materials 0.000 claims description 4
- 238000004062 sedimentation Methods 0.000 claims description 4
- 239000010802 sludge Substances 0.000 claims description 4
- 238000010606 normalization Methods 0.000 claims description 3
- 238000009827 uniform distribution Methods 0.000 claims 1
- 239000000460 chlorine Substances 0.000 description 10
- -1 nitrogen-containing organic compounds Chemical class 0.000 description 5
- QGZKDVFQNNGYKY-UHFFFAOYSA-N Ammonia Chemical compound N QGZKDVFQNNGYKY-UHFFFAOYSA-N 0.000 description 4
- XKMRRTOUMJRJIA-UHFFFAOYSA-N ammonia nh3 Chemical compound N.N XKMRRTOUMJRJIA-UHFFFAOYSA-N 0.000 description 4
- 229910052801 chlorine Inorganic materials 0.000 description 4
- QJGQUHMNIGDVPM-UHFFFAOYSA-N nitrogen group Chemical group [N] QJGQUHMNIGDVPM-UHFFFAOYSA-N 0.000 description 4
- 239000002351 wastewater Substances 0.000 description 4
- 230000001105 regulatory effect Effects 0.000 description 3
- QGZKDVFQNNGYKY-UHFFFAOYSA-O Ammonium Chemical compound [NH4+] QGZKDVFQNNGYKY-UHFFFAOYSA-O 0.000 description 2
- 229910021529 ammonia Inorganic materials 0.000 description 2
- 239000000126 substance Substances 0.000 description 2
- JVMRPSJZNHXORP-UHFFFAOYSA-N ON=O.ON=O.ON=O.N Chemical compound ON=O.ON=O.ON=O.N JVMRPSJZNHXORP-UHFFFAOYSA-N 0.000 description 1
- 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 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000005842 biochemical reaction Methods 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 238000004939 coking Methods 0.000 description 1
- 238000012937 correction Methods 0.000 description 1
- 238000000354 decomposition reaction Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 230000007613 environmental effect Effects 0.000 description 1
- 239000003344 environmental pollutant Substances 0.000 description 1
- 239000003337 fertilizer Substances 0.000 description 1
- 239000010842 industrial wastewater Substances 0.000 description 1
- 150000002500 ions Chemical class 0.000 description 1
- 238000012423 maintenance Methods 0.000 description 1
- 239000000463 material Substances 0.000 description 1
- 125000001477 organic nitrogen group Chemical group 0.000 description 1
- 230000003647 oxidation Effects 0.000 description 1
- 238000007254 oxidation reaction Methods 0.000 description 1
- 238000010238 partial least squares regression Methods 0.000 description 1
- 231100000719 pollutant Toxicity 0.000 description 1
- 229920005594 polymer fiber Polymers 0.000 description 1
- 238000005070 sampling Methods 0.000 description 1
- 230000001131 transforming effect Effects 0.000 description 1
- 238000004065 wastewater treatment Methods 0.000 description 1
Images
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
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2111/00—Details relating to CAD techniques
- G06F2111/10—Numerical modelling
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Activated Sludge Processes (AREA)
Abstract
本发明公开一种基于KPLSR模型的污水处理出水总氮浓度软测量方法,首先通过差分进化算法优化求解近邻成分分析的目标函数,从而实施不依赖于经验的辅助变量选择;其次,本发明方法通过融合利用不同核函数的KPLSR模型,来避免KPLSR模型核函数的选择问题;最后,利用选择的辅助变量实施基于KPLSR模型的出水总氮浓度实时软测量。本发明方法的最大优势在于:不依赖于污水处理过程的机理知识以及技术人员的主观经验,直接通过数据驱动的角度实现了特征变量的优选。此外,充分利用三类经典的核函数,分别建立三个不同的KPLSR模型,并进一步通过最小二乘回归将三个不同KPLSR模型的输出进行合并,保证了软测量结果的精度。
Description
技术领域
本发明涉及一种污水出水总氮浓度的软测量方法,特别涉及一种基于KPLSR模型的污水处理出水总氮浓度软测量方法。
背景技术
近年来,随着经济的发展,越来越多含氮污染物的任意排放给环境造成了极大的危害。氮在污水中以有机态氮、氨态氮(NH4+-N)、硝态氮(NO3--N)以及亚硝态氮(NO2--N)等多种形式存在,而氨态氮是最主要的存在形式之一,污水中的总氮浓度就是这些含氮有机化合物的浓度总和。污水中的氨氮是指以游离氨和离子铵形式存在的氮,主要来源于生活污水中含氮有机物的分解,焦化、合成氨等工业废水等。氨氮污染源多,并且还会散发出恶臭味,会直接影响人们的工作与生活。随着对环境问题的愈发重视,已经相应的规章制度规定经过污水处理后的出水总氮浓度的上限。
污水处理厂或污水处理系统为了降低污水出水口的总氮含量,通常是在曝气池氧化沟中进行投料,通过生物和化学反应将污水中的含氮成分转化为可利用的其他物质,比如:农作物的肥料。由于污水的来源众多,且总氮浓度波动较大,需要实时测量出水口总氮含量,从而及时调控曝气池中的投料量,从而使排放的污水中总氮浓度符合相应法律法规的要求。然而,现阶段能直接快速测量污水中总氮浓度的仪器仪表购置价格高昂,后期使用与维护成分同样是一笔不小的开支。在当前自动化程度普遍较低的污水处理厂中,通常不会采购高价格的仪器仪表设备。取而代之的是,通过人工取样的手段,进行总氮浓度的测量。虽然这种技术手段精确度较高,但是主要依赖人工操作,实时性较差。
在这种背景下,利用历史数据实施出水总氮浓度的软测量就被提上议程了。出水总氮浓度软测量的基本思想是利用其它测量数据,包括:入水流量、污水温度、悬浮物固体浓度、投料量等等,实时预测出水口的总氮浓度。然而,软测量技术实施成功与否关键在于,如何选择用于预测出水总氮浓度的辅助变量。虽然,污水处理厂可每天需要测量的数据很多,但是只有部分变量在建立软测量模型时能起到关键的作用。考虑到污水处理过程的生化反应复杂性,依靠经验选择辅助变量通常是不可靠。若是经验不可用或经验有误,则会导致出水总氮浓度测量精度达不到要求。
发明内容
本发明所要解决的主要技术问题是:如何选择有益的辅助变量以及如何确定核偏最小二乘回归(Kernel Partial Least Sqaure Regression,缩写:KPLSR)模型的参数,从而提升出水总氮浓度软测量模型的精度。具体来讲,本发明方法首先通过差分进化算法优化求解近邻成分分析的目标函数,从而实施不依赖于经验的辅助变量选择;其次,本发明方法通过融合利用不同核函数的KPLSR模型,来避免KPLSR模型参数的选择问题;最后,利用选择的辅助变量实施基于KPLSR模型的出水总氮浓度实时软测量。
本发明方法解决上述问题所采用的技术方案为:一种基于KPLSR模型的污水处理出水总氮浓度软测量方法,包括以下所示步骤:
步骤(1):确定污水处理厂每天可测量的变量;其中,污水处理厂入水处可测量的变量依次包括:入水流量,入水温度,入水色度,入水氯离子浓度,入水悬浮物固体浓度,入水PH值;曝气池中可测量的变量依次包括:投料量,污水色度,污水悬浮物固体浓度,污水PH值,污水氯离子浓度,生化需氧量;沉淀池中可测量的变量依次是:污水色度,污水悬浮物固体浓度,污水PH值,污水氯离子浓度,污泥量;污水处理厂出水处可测量的变量具体包括:出水色度,出水PH值,出水悬浮物固体浓度,和出水总氮浓度。
步骤(2):连续采集n天的数据,并将每天的数据存储为一个m×1维的数据向量xi;其中,第i天的数据向量xi中的元素是按照步骤(1)中所述变量的先后顺序排列,m等于步骤(1)中可测量的变量总数,i∈{1,2,…,n}。
步骤(3):将n天的数据向量x1,x2,…,xn组建成一个数据矩阵X=[x1,x2,…,xn]T后,对X∈Rn×m中的各个列向量分别实施标准化处理,得到标准化后的数据矩阵并将X中第m列的列向量的均值与标准差分别记录为μ与σ;其中,上标号T表示矩阵或向量的转置,Rn ×m表示n×m维的实数矩阵,R表示实数集。
步骤(4.1):初始化迭代次数g=1并确定差分进化优化算法的参数,具体包括:种群个数N,缩放因子zf,交叉概率cp,最大迭代次数G。
步骤(4.2):随机产生N个(m-1)×1维的种群向量w1,w2,…,wN,每个种群向量中的元素都按照均匀分布随机取值于区间[-1,1]。
步骤(4.3):计算种群向量w1,w2,…,wN分别对应的目标函数值F1,F2,…,FN;其中,计算第c个种群向量wc对应的目标函数值Fc的具体实施过程如步骤(4.3-1)至步骤(4.3-4)所示。
其中,i∈{1,2,…,n},j∈{1,2,…,n},diag{wc}表示将种群向量wc转变为一个对角矩阵,符号||||表示计算向量的长度,c∈{1,2,…,N}。
上式中,exp()表示以自然常数e为底数的指数函数。
步骤(4.3-3):根据如下所示公式计算输出概率误差f1,f2,…,fn:
步骤(4.3-4):计算第c个种群向量wc对应的目标函数值Fc=f1+f2+…+fn。
步骤(4.4):将F1,F2,…,FN中的最大值及其对应的种群向量分别记录为Fbest和wbest后,执行差分进化算法的种群更新操作,得到更新后的N个种群向量w1,w2,…,wN,具体的实施过程如步骤(4.4-1)至步骤(4.4-4)所示。
步骤(4.4-1):根据如下所示公式为每个种群向量产生一个对应的变异向量vc:
vc=wc+zf×(wbest-wc)+zf×(wa-wb) ④
上式中,下标号a与b是从区间[1,N]中随机产生的2个互不相等的整数。
步骤(4.4-2):按照如下所示公式对变异向量vc进行修正:
上式中,vc(k)表示变异向量vc中的第k个元素,k∈{1,2,…,m-1}。
步骤(4.4-3):根据如下所示公式产生N个尝试向量u1,u2,…,uN:
其中,uc(k)与wc(k)分别为uc与wc中的第k个元素,表示0至1之间的随机数。
步骤(4.4-4):依据如下所示公式更新种群向量w1,w2,…,wN:
步骤(4.5):判断是否满足条件g>G;若否,则设置g=g+1后返回步骤(4.3);若是,则得到最优的种群向量wbest。
步骤(4.6):按照如下所示公式确定出选择向量v∈R1×(m-1)中的各个元素:
上式中,v(k)和wbest(k)分别表示v和wbest中的第k个元素,η为阈值,其建议取值范围是0.01≤η≤0.1。
步骤(6):分别根据如下所示的三个核函数G1(ε1,ε2),G2(ε1,ε2),和G3(ε1,ε2),建立三个不同的KPLSR模型后,计算得到输出估计向量具体的实施过程如步骤(6.1)至步骤(6.6)所示。
其中,ε1和ε2表示核函数的两个输入行向量,α和β为核函数的参数。
步骤(6.1):根据公式Kd(i,j)=Gd(zi,zj)计算核矩阵Kd∈Rn×n中的第i行第j列元素Kd(i,j);其中,d∈{1,2,3},zi和zj分别表示输入矩阵Z中的第i行和第j行的行向量。
其中,方阵Θ∈Rn×n中所有元素都等于1。
步骤(6.3):依次设置d=1,2和3,并执行如下所示步骤(6.3-1)至步骤(6.3-4),从而得到KPLSR模型的核回归系数向量θ1,θ2,θ3。
步骤(6.3-3):判断是否满足条件:若是,则将得到的γ个得分向量s1,s2,…,sγ组成核得分矩阵Sd=[s1,s2,…,sγ];若否,则根据公式更新K0,再设置γ=γ+1后,返回步骤(6.3-2);其中,δ表示阈值,其建议的取值范围是:0.0001≤δ≤0.01。
步骤(8):采集污水处理厂新一天的测量数据,并将其存储为一个1×(m-1)维的数据向量xt;其中,数据向量xt中的元素是按照步骤(1)中所述的前m-1个变量的先后顺序依次排列,下标号t表示最新的一天。
步骤(10):调用步骤(4)中辅助变量的选择向量v∈R1×(m-1),并根据v中等于1的元素所在的列,对应的将xt中相同列的元素组成输入向量zt。
步骤(11):根据公式kd(i)=Gd(zt,zi)计算核向量kd∈R1×n中的第i个元素kd(i),重复本步骤直至得到三个核向量k1,k2,k3。
步骤(12):根据如下所示步骤(12.1)至步骤(12.3)计算得到出水总氮浓度的软测量值yt。
其中,φ∈R1×n中所有元素都等于1。
通过以上所述实施步骤,本发明方法的优势介绍如下。
本发明方法的最大优势在于:不依赖于污水处理过程的机理知识以及技术人员的主观经验,直接通过数据驱动的角度实现了特征变量的优选。此外,充分利用三类经典的核函数,分别建立三个不同的KPLSR模型,避免了核函数的选择问题。最后,进一步通过最小二乘回归将三个不同KPLSR模型的输出进行合并,对于预测结果的精度能有较充分的保障。
附图说明
图1为本发明方法的实施流程示意图。
图2为污水处理厂的流程示意图。
图3为出水总氮浓度的软测量与实测值之间的误差对比图。
具体实施方式
下面结合附图和具体实施方式对本发明进行详细说明。
如图1所示,本发明公开了一种基于KPLSR模型的污水处理出水总氮浓度软测量方法,下面结合一个具体应用实例来说明本发明方法的具体实施方式。
本实施案例中的污水处理厂的流程如图2所示。在本实施案例中,共连续采集了212天的数据;其中,将前112天的数据向量用于软测量模型的离线建模阶段,后100天的数据向量分别被当成新一天采样的数据向量,用于验证出水总氮的软测量精度。首先,实施如下所示步骤(1)至步骤(7),从而建立出水总氮的软测量模型。
步骤(1):确定污水处理厂每天可测量的变量:污水处理厂入水处可测量的变量具体包括:入水流量Qin,入水温度Tin,入水色度Cin,入水氯离子浓度Clin,入水悬浮物固体浓度Din,入水PH值Pin;曝气池中可测量的变量具体包括:投料量Qa,污水色度Ca,污水悬浮物固体浓度Da,污水PH值Pa,污水氯离子浓度Cla,生化需氧量Oa;沉淀池中可测量的变量具体包括:污水色度Cs,污水悬浮物固体浓度Ds,污水PH值Ps,污水氯离子浓度Cls,污泥量Vs;污水处理厂出水处可测量的变量具体包括:出水色度Co,出水PH值Po,出水悬浮物固体浓度Do,和出水总氮浓度TN。
步骤(2):连续采集n天的数据,并将每天的数据存储为一个m×1维的数据向量xi;其中,m表示步骤(1)中可测量的变量总数。
步骤(3):将n天的数据向量x1,x2,…,xn组建成一个数据矩阵X=[x1,x2,…,xn]T后,对X∈Rn×m中的各个列向量分别实施标准化处理,得到标准化后的数据矩阵并将X中最后一列向量的均值与标准差分别记录为μ与σ。
步骤(6):分别根据如下所示的三个核函数G1(ε1,ε2),G2(ε1,ε2),和G3(ε1,ε2),建立三个不同的KPLSR模型后,计算得到输出估计向量具体的实施过程如前述步骤(6.1)至步骤(6.6)所示。
完成上述步骤(1)至步骤(7)后,即可实施出水总氮浓度的软测量,即实施如下所示步骤(8)至步骤(12)。
步骤(8):采集污水处理厂新一天的m-1个测量数据,并将其存储为一个1×(m-1)维的数据向量xt;其中,下标号t表示最新的一天,数据向量xt中的元素的先后排列顺序依次为:入水流量Qin,入水温度Tin,入水色度Cin,入水氯离子浓度Clin,入水悬浮物固体浓度Din,入水PH值Pin;曝气池中可测量的变量具体包括:投料量Qa,污水色度Ca,污水悬浮物固体浓度Da,污水PH值Pa,污水氯离子浓度Cla,生化需氧量Oa;沉淀池中可测量的变量具体包括:污水色度Cs,污水悬浮物固体浓度Ds,污水PH值Ps,污水氯离子浓度Cls,污泥量Vs;污水处理厂出水处可测量的变量具体包括:出水色度Co,出水PH值Po,出水悬浮物固体浓度Do。
步骤(10):根据选择向量v∈R1×(m-1)中等于1的元素所在列,对应的将xt中相同列的元素组建成输入向量zt。
步骤(11):根据公式kd(i)=Gd(zt,zi)计算核向量kd∈R1×n中的第i个元素kd(i),重复本步骤直至得到三个核向量k1,k2,k3。
步骤(12):通过步骤(6)中的三个KPLSR模型以及步骤(7)中的回归系数向量θo,计算得到出水总氮浓度的软测量值yt,具体的实施过程如步骤(12.1)至步骤(12.3)所示。
通过对比步骤(12)中的出水总氮浓度的软测量值与实际通过化学手段测量得到的出水总氮浓度测量值进行对比,相应的对比结果展示于图3中。从图3中可以看出,本发明方法具备较强的软测量精度,误差变化在可接受变化范围内。
Claims (4)
1.一种基于KPLSR模型的污水处理出水总氮浓度软测量方法,其特征在于,包括以下步骤:
步骤(1):确定污水处理厂每天可测量的变量;其中,污水处理厂入水处可测量的变量依次包括:入水流量,入水温度,入水色度,入水氯离子浓度,入水悬浮物固体浓度,入水PH值;曝气池中可测量的变量依次包括:投料量,污水色度,污水悬浮物固体浓度,污水PH值,污水氯离子浓度,生化需氧量;沉淀池中可测量的变量依次是:污水色度,污水悬浮物固体浓度,污水PH值,污水氯离子浓度,污泥量;污水处理厂出水处可测量的变量具体包括:出水色度,出水PH值,出水悬浮物固体浓度,和出水总氮浓度;
步骤(2):连续采集n天的数据,并将每天的数据存储为一个m×1维的数据向量xi;其中,第i天的数据向量xi中的元素是按照步骤(1)中所述变量的先后顺序排列,m等于步骤(1)中可测量的变量总数,i∈{1,2,…,n};
步骤(3):将n天的数据向量x1,x2,…,xn组建成一个数据矩阵X=[x1,x2,…,xn]T后,对X∈Rn×m中的各个列向量分别实施标准化处理,得到标准化后的数据矩阵并将X中第m列的列向量的均值与标准差分别记录为μ与σ;其中,上标号T表示矩阵或向量的转置,Rn×m表示n×m维的实数矩阵,R表示实数集;
步骤(4.1):初始化迭代次数g=1,确定差分进化优化算法的参数,具体包括:种群个数N,缩放因子zf,交叉概率cp,最大迭代次数G;
步骤(4.2):随机产生N个(m-1)×1维的种群向量w1,w2,…,wN,每个种群向量中的元素都按照均匀分布随机取值于区间[-1,1];
步骤(4.3):计算种群向量w1,w2,…,wN分别对应的目标函数值F1,F2,…,FN;
步骤(4.4):将F1,F2,…,FN中的最大值及其对应的种群向量分别记录为Fbest和wbest后,执行差分进化算法的种群更新操作,得到更新后的N个种群向量w1,w2,…,wN;
步骤(4.5):判断是否满足条件g>G;若否,则设置g=g+1后返回步骤(4.3);若是,则得到最优的种群向量wbest;
步骤(4.6):按照如下所示公式确定出选择向量v∈R1×(m-1)中的各个元素:
上式中,v(k)和wbest(k)分别表示v和wbest中的第k个元素,η为阈值;
步骤(6):分别根据如下所示的三个核函数G1(ε1,ε2),G2(ε1,ε2),和G3(ε1,ε2),建立三个不同的KPLSR模型后,计算得到输出估计向量具体的实施过程如步骤(6.1)至步骤(6.6)所示;
其中,ε1和ε2表示核函数的两个输入行向量,α和β为核函数的参数;
步骤(6.1):根据公式Kd(i,j)=Gd(zi,zj)计算核矩阵Kd∈Rn×n中的第i行第j列元素Kd(i,j);其中,d∈{1,2,3},zi和zj分别表示输入矩阵Z中的第i行和第j行的行向量;
其中,方阵Θ∈Rn×n中所有元素都等于1;
步骤(6.3):依次设置d=1,2和3,并执行如下所示步骤(6.3-1)至步骤(6.3-4),从而得到KPLSR模型的核回归系数向量θ1,θ2,θ3;
步骤(8):采集污水处理厂新一天的m-1个测量数据,并将其存储为一个1×(m-1)维的数据向量xt;其中,数据向量xt中的元素按照步骤(1)中所述的前m-1个变量的先后顺序依次排列,下标号t表示最新的一天;
步骤(10):利用步骤(4)中辅助变量的选择向量v∈R1×(m-1),根据v中等于1的元素所在的列,对应的将xt中相同列的元素组成输入向量zt;
步骤(11):根据公式kd(i)=Gd(zt,zi)计算核向量kd∈R1×n中的第i个元素kd(i),重复本步骤直至得到三个核向量k1,k2,k3;
步骤(12):根据如下所示步骤(12.1)至步骤(12.3)计算得到出水总氮浓度的软测量值yt;
其中,φ∈R1×n中所有元素都等于1;
2.根据权利要求1所述的一种基于KPLSR模型的污水处理出水总氮浓度软测量方法,其特征在于,所述步骤(4.3)中计算第c个种群向量wc对应的目标函数值Fc的具体实施过程如下所示:
上式中,i∈{1,2,…,n},j∈{1,2,…,n},diag{wc}表示将种群向量wc转变为一个对角矩阵,符号|| ||表示计算向量的长度,c∈{1,2,…,N};
上式中,exp( )表示以自然常数e为底数的指数函数;
步骤(4.3-3):根据如下所示公式计算输出概率误差f1,f2,…,fn:
步骤(4.3-4):计算第c个种群向量wc对应的目标函数值Fc=f1+f2+…+fn。
3.根据权利要求1所述的一种基于KPLSR模型的污水处理出水总氮浓度软测量方法,其特征在于,所述步骤(4.4)中差分进化优化算法的种群更新操作的具体实施过程如下所示:
步骤(4.4-1):根据如下所示公式为每个种群向量产生一个对应的变异向量vc:
vc=wc+zf×(wbest-wc)+zf×(wa-wb) ⑩
上式中,下标号a与b是从区间[1,N]中随机产生的2个互不相等的整数;
步骤(4.4-2):按照如下所示公式对变异向量vc进行修正:
上式中,vc(k)表示变异向量vc中的第k个元素,k∈{1,2,…,m-1};
步骤(4.4-3):根据如下所示公式产生N个尝试向量u1,u2,…,uN:
其中,uc(k)与wc(k)分别为uc与wc中的第k个元素,rk表示0至1之间的随机数;
步骤(4.4-4):依据如下所示公式更新种群向量w1,w2,…,wN:
4.根据权利要求1所述的一种基于KPLSR模型的污水处理出水总氮浓度软测量方法,其特征在于,所述步骤(6.3)中建立KPLSR模型的具体实施过程如下所示:
步骤(6.3-3):判断是否满足条件:若是,则将得到的γ个得分向量s1,s2,…,sγ组成核得分矩阵Sd=[s1,s2,…,sγ];若否,则根据公式更新K0,并设置γ=γ+1后,再返回步骤(6.3-2);其中,δ表示阈值;
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110300719.7A CN113065242A (zh) | 2021-03-13 | 2021-03-13 | 一种基于kplsr模型的污水处理出水总氮浓度软测量方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110300719.7A CN113065242A (zh) | 2021-03-13 | 2021-03-13 | 一种基于kplsr模型的污水处理出水总氮浓度软测量方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN113065242A true CN113065242A (zh) | 2021-07-02 |
Family
ID=76562858
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202110300719.7A Pending CN113065242A (zh) | 2021-03-13 | 2021-03-13 | 一种基于kplsr模型的污水处理出水总氮浓度软测量方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN113065242A (zh) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN115952728A (zh) * | 2022-11-29 | 2023-04-11 | 淮阴工学院 | 一种亚硝态氮浓度软测量预测方法 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104914227A (zh) * | 2015-06-16 | 2015-09-16 | 华南理工大学 | 基于多高斯核自优化相关向量机的污水水质软测量方法 |
CN107085372A (zh) * | 2017-05-10 | 2017-08-22 | 湖南工业大学 | 一种基于改进萤火虫算法和最小二乘支持向量机的污水节能处理优化控制方法 |
CN110033175A (zh) * | 2019-03-12 | 2019-07-19 | 宁波大学 | 一种基于集成多核偏最小二乘回归模型的软测量方法 |
CN110929809A (zh) * | 2019-12-14 | 2020-03-27 | 北京工业大学 | 一种特征自增强的循环神经网络的污水关键水质指标软测量方法 |
-
2021
- 2021-03-13 CN CN202110300719.7A patent/CN113065242A/zh active Pending
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104914227A (zh) * | 2015-06-16 | 2015-09-16 | 华南理工大学 | 基于多高斯核自优化相关向量机的污水水质软测量方法 |
CN107085372A (zh) * | 2017-05-10 | 2017-08-22 | 湖南工业大学 | 一种基于改进萤火虫算法和最小二乘支持向量机的污水节能处理优化控制方法 |
CN110033175A (zh) * | 2019-03-12 | 2019-07-19 | 宁波大学 | 一种基于集成多核偏最小二乘回归模型的软测量方法 |
CN110929809A (zh) * | 2019-12-14 | 2020-03-27 | 北京工业大学 | 一种特征自增强的循环神经网络的污水关键水质指标软测量方法 |
Non-Patent Citations (2)
Title |
---|
刘鸿斌 等: "废水处理过程的KPLS-GPR软测量建模", 江苏大学学报(自然科学版), no. 05, 10 September 2020 (2020-09-10), pages 1 - 6 * |
杨冲: "基于改进PLS算法的造纸废水处理过程故障检测与水质预测", 中国优秀硕士论文全文库工程科技Ⅰ辑, no. 1, 15 January 2021 (2021-01-15), pages 1 - 76 * |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN115952728A (zh) * | 2022-11-29 | 2023-04-11 | 淮阴工学院 | 一种亚硝态氮浓度软测量预测方法 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US10570024B2 (en) | Method for effluent total nitrogen-based on a recurrent self-organizing RBF neural network | |
CN102854296B (zh) | 一种基于集成神经网络的污水处理软测量方法 | |
CN111291937A (zh) | 基于支持向量分类与gru神经网络联合的处理污水水质预测方法 | |
Weijers | Modelling, identification and control of activated sludge plants for nitrogen removal | |
CN111354423A (zh) | 一种基于多元时间序列分析的自组织递归模糊神经网络的出水氨氮浓度预测方法 | |
CN111650834B (zh) | 基于极限学习机的污水处理过程预测控制方法 | |
CN112884056A (zh) | 基于优化的lstm神经网络的污水水质预测方法 | |
CN110320335B (zh) | 一种关于污水处理出水水质指标的多元鲁棒软测量方法 | |
CN110542748B (zh) | 一种基于知识的鲁棒型出水氨氮软测量方法 | |
CN103632032A (zh) | 一种城市污水处理过程出水指标在线软测量预测方法 | |
CN112989704B (zh) | 一种基于de算法的irfm-cmnn出水bod浓度预测方法 | |
CN109657790A (zh) | 一种基于pso的递归rbf神经网络出水bod预测方法 | |
CN108549740A (zh) | 一种基于混合智能算法的厌氧系统出水氨氮软测量方法 | |
CN113065242A (zh) | 一种基于kplsr模型的污水处理出水总氮浓度软测量方法 | |
Baeza et al. | In-line fast OUR (oxygen uptake rate) measurements for monitoring and control of WWTP | |
CN201330211Y (zh) | 污水处理厂运行参数自寻优模拟系统 | |
Lim et al. | A systematic model calibration methodology based on multiple errors minimization method for the optimal parameter estimation of ASM1 | |
CN113222324B (zh) | 一种基于pls-pso-rbf神经网络模型的污水质量监测方法 | |
CN114861543A (zh) | 一种数据驱动的石化污水可生化性智能评价方法 | |
CN113051818B (zh) | 一种基于特征变量极限学习机的排污口磷含量测量方法 | |
CN113221436A (zh) | 一种基于改进型rbf神经网络的污水悬浮物浓度软测量方法 | |
Jeppsson | On the verifiability of the activated sludge system dynamics | |
Krivec et al. | Integrated theoretical and data-driven Gaussian Process NARX Model for the Simulation of Effluent Concentrations in Wastewater Treatment Plant | |
Liu et al. | Neural network identification of wastewater treatment plants | |
Al-Asheh et al. | Forecasting influent-effluent wastewater treatment plant using time series analysis and artificial neural network techniques |
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 |