CN104933310A - 基于ga-lssvm模型预测膜蒸馏水处理系统中膜污染趋势的方法 - Google Patents

Abstract

本发明公开了基于GA-LSSVM模型预测膜蒸馏水处理系统中膜污染趋势的方法，首先利用LSSVM算法建立膜蒸馏污水处理过程的预测模型；其次分别在拟稳定状态和非稳定状态下，利用GA算法优化预测模型的参数；然后分别在拟稳定状态和非稳定状态下，利用优化的预测模型对膜通量以及膜污染阻力的变化趋势进行预测，并分析膜蒸馏的基本操作参数对膜通量和膜污染阻力的影响；最后对预测结果进行灵敏度分析计算，确定出影响膜通量和膜污染阻力的主导因素。本发明利用GA-LSSVM实时预测膜通量以及膜污染阻力变化情况，将膜蒸馏基本操作参对膜污染的影响明晰化、量化。

Description

基于GA-LSSVM模型预测膜蒸馏水处理系统中膜污染趋势的方法

技术领域

本发明属于膜蒸馏技术、污水废水处理领域，特别涉及了基于GA-LSSVM模型预测膜蒸馏水处理系统中膜污染趋势的方法。

背景技术

近年来，膜分离技术发展迅速，包括纳滤膜、超滤膜、反渗透膜、膜蒸馏，它们的快速发展为城市污水处理和工业废水处理提供新途径。膜蒸馏是基于气液平衡及传热传质原理的热驱动过程，以微孔疏水膜为介质，在膜两侧蒸气压差的作用下，进水中的挥发性组分以蒸汽形式透过膜孔，在膜的低温侧冷凝为液态，而非挥发性组分留在热侧进料液中，从而达到分离并去除污染物的目的。膜蒸馏技术应用于污水处理过程的优点在于：1、膜蒸馏浓缩过程在常压下进行，设备简单，操作方便；2、在非挥发性水溶液的膜蒸馏过程中，仅有水蒸气能透过疏水膜孔，因此出水水质好；3、与反渗透相比，盐浓度以及浓差极化对膜蒸馏影响较小，因此可以处理含有高浓度无机盐的溶液，甚至可以将溶液浓缩到过饱和状态。

然而，利用膜蒸馏技术处理污水最主要的限制来源于膜污染，膜污染会严重导致膜通量的降低，膜蒸馏处理效率的下降。此外，膜污染还将增加额外的膜清洗费用，缩短膜的使用寿命，并且因污水处理过程中需要频繁地更换新膜，而增加了处理费用。膜的运行过程中，其污染程度主要取决于3个因素：1、系统的运行参数；2、进水水质；3、膜本身性能。由于膜分离过程的非线性特征以及膜污染与系统操作参数之间关系的模糊性，导致了膜蒸馏处理过程中，其基本操作参数(比如：温度差、错流速度、进水pH、膜孔径、溶质粒子的粒径大小等)对膜污染的具体量化的影响还不明晰，有待进一步的研究。因此，为探明系统操作参数对膜污染的具体影响并将这种影响进行量化，需对膜蒸馏处理过程进行数学模拟，即建立膜蒸馏污水处理过程的预测模型。

目前，用于预测膜分离过程的数学模型大致分为两类：1、理论模型(参数模型)；2、非参数模型；其中理论模型是基于机理知识的(即基本参数对膜污染的影响机理)，要求对相应的膜污染影响机理非常了解，才可以应用理论模型进行膜污染预测，但运算复杂、计算处理费用昂贵。而非参数模型则不需要对相应机理(基本参数对膜污染的影响机理)了解得很透彻，弥补了理论模型存在的弊端。由于膜蒸馏处理过程的复杂性以及基本操作参数对膜污染影响机理的不明确性，因此需要选取合适的非参数模型对膜蒸馏过程进行预测。目前，应用较多的算法模型有很多，包括GP(遗传规划算法)、模糊系统、ANN(人工神经网络)等，其中当前应用最多的是ANN模型，它被广泛运用于各种膜分离处理过程(例如：微滤MF、超滤UF、反渗透RO、纳滤NF、电渗析ED),其中ANN还成功用于膜蒸馏过程的模拟预测。但是ANN存在以下缺点：1、预测结果仅为局部极小值，而非全局极小值；2、难以确定其网络拓扑结构；3、易出现过度拟合问题。

发明内容

为了解决上述背景技术提出的技术问题，本发明旨在提供基于GA-LSSVM模型预测膜蒸馏水处理系统中膜污染趋势的方法，利用GA-LSSVM实时预测膜通量以及膜污染阻力变化情况，将膜蒸馏基本操作参对膜污染的影响明晰化、量化。

为了实现上述技术目的，本发明的技术方案为：

基于GA-LSSVM模型预测膜蒸馏水处理系统中膜污染趋势的方法，包括以下步骤：

(1)利用LSSVM算法建立膜蒸馏污水处理过程的预测模型；

(2)分别在拟稳定状态和非稳定状态下，利用GA算法优化步骤(1)建立的预测模型的参数；

(3)分别在拟稳定状态和非稳定状态下，利用步骤(2)优化的预测模型对膜通量以及膜污染阻力的变化趋势进行预测，并分析膜蒸馏的基本操作参数对膜通量和膜污染阻力的影响；

(4)对步骤(3)得到的预测结果进行灵敏度分析计算，确定出影响膜通量和膜污染阻力的主导因素。

步骤(1)的具体过程如下：

首先，输入数据记为X，输出数据记为Y，这些数据分为训练数据和测试数据，所述训练数据的集合记为A＝{(x₁,y₁),…,(x_i,y_i),…,(x_N,y_N)}，其中，x_i∈X，y_i∈Y，1≤i≤N，N为训练数据集合的大小，利用非线性映射函数建立如下式(1)的回归模型，将输入数据X映射到高维特征空间，

式(1)中，w为权值，b为偏置项，为非线性映射函数；

然后，建立式(2)所示的LSSVM优化问题，

$\min J(w,b,\mathrm{\ξ})=\frac{1}{2}{w}^{T}w+\frac{1}{2}\mathrm{\γ}\underset{i=1}{\overset{N}{\Σ}}{\mathrm{\ξ}}_i^2---\left(2\right)$

式(2)中，ξ_i为松弛因子，ξ_i＞0，γ为正则常数；

最后，构造式(3)所示的拉格朗日函数，求解LSSVM最优化问题，

式(3)中，a_i为拉格朗日乘数；分别对w,b,ξ_i,a_i求偏导数，并令偏导数为0，消去w和ξ_i，得式(4)，

$\left[\begin{array}{cc}0& 1_N^T\\ I_N& \mathrm{\Ω}+{\mathrm{\γ}}^{-1}{I}_N\end{array}\right]\left[\begin{array}{c}b\\ a\end{array}\right]=\left[\begin{array}{c}0\\ Y\end{array}\right]---\left(4\right)$

式(4)中，Y＝[y₁,...,y_N]^T，1_N＝[1,...,1]^T，a＝[a₁,...,a_N]^T，I_N为N×N的单位矩阵，Ω为N×N的矩阵，它的第i行第j列的元素κ(x_i,x_j)是内核函数。

所述内核函数采用径向基函数：

$\mathrm{\κ}(x_i,x_j)=\exp \left(\frac{-\left|\right|x_i-x_j^2\left|\right|}{{\mathrm{\σ}}^2}\right)---\left(5\right)$

式(5)中，σ²为内核参数。

步骤(2)中所述利用GA算法优化的参数包括正则常数γ和内核参数σ²。

步骤(3)的具体过程为，分别在拟稳定状态和非稳定状态下进行如下操作：

将温度差、错流速度、进水pH、膜孔径和膜蒸馏处理时间作为预测模型的输入数据，将膜通量和膜污染阻力作为预测模型的输出数据，每次仅改变温度差、错流速度、进水pH、膜孔径和膜蒸馏处理时间中的一个参数，保证其他参数不变，运用GA算法优化LSSVM模型的参数，并将优化后的参数代入LSSVM模型进行预测，再利用均方根误差RMSE和回归系数R²这两个评价参数对预测性能进行评价，当RMSE越接近于0且R²越接近于1时，表明预测性能越好。

在步骤(4)中，分别采用斯皮尔曼相关性、伽马相关性、肯德尔相关性和皮尔逊相关性对预测结果进行灵敏度分析计算。

采用上述技术方案带来的有益效果：

本发明采用最小二乘支持向量机LSSVM来建立预测模型，克服了常用模型的缺点，其适应性和预测性能较好，预测的速度也更快；本发明选用遗传算法GA对LSSVM模型参数进行优化，具有既可处理离散变量又可处理连续变量、针对非线性目标以及不需要梯度信息来约束函数等优点，优化精度较高；本发明还对预测结果进行灵敏度分析计算，筛选出影响膜通量以及膜污染阻力的主导因素。

附图说明

图1是本发明的基本流程图；

图2是本发明GA算法流程图。

具体实施方式

以下将结合附图，对本发明的技术方案进行详细说明。

如图1所示本发明的基本流程图，基于GA-LSSVM模型预测膜蒸馏水处理系统中膜污染趋势的方法，包括以下步骤：

第一步：利用LSSVM算法建立膜蒸馏污水处理过程的预测模型；其具体过程如下：

首先，输入数据记为X，输出数据记为Y，这些数据分为训练数据和测试数据，所述训练数据的集合记为A＝{(x₁,y₁),…,(x_i,y_i),…,(x_N,y_N)}，其中，x_i∈X，y_i∈Y，1≤i≤N，N为训练数据集合的大小，利用非线性映射函数建立如下式(1)的回归模型，将输入数据X映射到高维特征空间，

式(1)中，w为权值，b为偏置项，为非线性映射函数；

然后，建立式(2)所示的LSSVM优化问题，

$\min J(w,b,\mathrm{\ξ})=\frac{1}{2}{w}^{T}w+\frac{1}{2}\mathrm{\γ}\underset{i=1}{\overset{N}{\Σ}}{\mathrm{\ξ}}_i^2---\left(2\right)$

式(2)中，ξ_i为松弛因子，ξ_i＞0，γ为正则常数，用于调整模型的复杂性和经验风险的折中；

最后，构造式(3)所示的拉格朗日函数，求解LSSVM最优化问题，

式(3)中，a_i为拉格朗日乘数；分别对w,b,ξ_i,a_i求偏导数，并令偏导数为0，消去w和ξ_i，得式(4)，

$\left[\begin{array}{cc}0& 1_N^T\\ I_N& \mathrm{\Ω}+{\mathrm{\γ}}^{-1}{I}_N\end{array}\right]\left[\begin{array}{c}b\\ a\end{array}\right]=\left[\begin{array}{c}0\\ Y\end{array}\right]---\left(4\right)$

式(4)中，Y＝[y₁,...,y_N]^T，1_N＝[1,...,1]^T，a＝[a₁,...,a_N]^T，I_N为N×N的单位矩阵，Ω为N×N的矩阵，它的第i行第j列的元素κ(x_i，x_j)是内核函数。

因径向基函数能够较好实现输入样本到高维特征空间的非线性映射，且需要的调节参数较少，因此本发明选用径向基函数作为内核函数：

$\mathrm{\κ}(x_i,x_j)=\exp \left(\frac{-\left|\right|x_i-x_j^2\left|\right|}{{\mathrm{\σ}}^2}\right)---\left(5\right)$

式(5)中，σ²为内核参数。

第二步：分别在拟稳定状态和非稳定状态下，利用GA算法优化第一步建立的预测模型的参数。在本发明中，需要对上述正则常数γ和内核参数σ²进行优化，优化的具体步骤如图2所示：

1、初始化：设置进化代数计算器t＝0，设置最大进化代数T，随机生成M个个体作为初始群体P(0)；

2、个体评价：计算群体P(t)中各个个体的适应度，在本发明中，利用平均值绝对相关偏差(AARD)这一适应度函数计算适应度；

3、选择运算：将选择算子作用于群体，其目的在于把优化的个体直接遗传到下一代或通过配对交叉产生新的个体再遗传到下一代。选择操作是建立在群体中个体的适应度评估基础上的；

4、交叉运算：将交叉算子作用于群体，遗传算法中起核心作用的就是交叉算子；

5、变异运算：将变异算子作用于群体，即对群体中的个体串的某些基因座上的基因值作变动。群体P(t)经过选择、交叉、变异运算之后得到下一代群体P(t+1)；

6、终止条件判断：若t＝T，则以进化过程中所得到的具有最大适应度个体作为最优解输出，终止计算。

第三步：分别在拟稳定状态和非稳定状态下，利用第二步优化的预测模型对膜通量以及膜污染阻力的变化趋势进行预测，并分析膜蒸馏的基本操作参数对膜通量和膜污染阻力的影响；其具体过程如下：

分别在拟稳定状态和非稳定状态下进行如下操作：

将温度差、错流速度、进水pH、膜孔径和膜蒸馏处理时间作为预测模型的输入数据，将膜通量和膜污染阻力作为预测模型的输出数据，每次仅改变温度差、错流速度、进水pH、膜孔径和膜蒸馏处理时间中的一个参数，保证其他参数不变，运用GA算法优化LSSVM模型的参数，并将优化后的参数代入LSSVM模型进行预测，再利用均方根误差RMSE和回归系数R²这两个评价参数对预测性能进行评价，当RMSE越接近于0且R²越接近于1时，表明预测性能越好。

第四步：对步骤(3)得到的预测结果进行灵敏度分析计算，确定出影响膜通量和膜污染阻力的主导因素。

本发明中分别利用4种灵敏度分析计算方法——斯皮尔曼相关性、伽马相关性、肯德尔相关性、皮尔逊相关性分别计算：1、各个输入变量(温度差ΔT、错流速度CFV、进水pH、膜孔径D)与膜通量之间的相关系数，从而确定出影响膜蒸馏污水处理过程中膜通量的主导因素；2、各个输入变量(温度差ΔT、错流速度CFV、进水pH、膜孔径D)与膜污染阻力之间的相关系数，从而确定出影响膜蒸馏污水处理过程中膜污染的主导因素；从而使得膜蒸馏基本操作参数(温度差ΔT、错流速度CFV、进水pH、膜孔径D)对膜污染、膜通量的影响程度得到量化跟明晰化。

在本实施例中，自配生活污水进行模拟试验：

膜蒸馏反应器的进水CODcr浓度为157.9mg/L，氨氮(NH₄ ⁺-N)浓度为31.8mg/L，电导率为0.978ms/cm，多糖浓度为1.96mg/L，蛋白质浓度为17.5mg/L。

试验运行中，膜蒸馏反应器热侧水温梯度设为40℃、60℃、80℃，进水蠕动泵转速梯度设为25r/s、50r/s、75r/s；冷侧水温为11℃，出水蠕动泵转速梯度设为25r/s、50r/s、75r/s。反应器持续运行16d，热侧污水初始pH分别为4、7、10,16d运行后其pH稳定在10左右。本模拟试验采用的是聚四氟乙烯PTFE材质的疏水膜，膜孔径D分别为0.22μm、1.20μm、2.0μm。

结果显示，在利用不同孔径的PTFE输水膜条件下，污水中COD的去除率均可高达95％以上，氨氮的去除率可高达75％以上，多糖的去除率可高达73％以上，蛋白质的去除率可高达98％以上。

模拟试验的具体结果如下：

(1)根据回归系数(R²)的大小判别GA-LSSVM模型的预测性能好坏。模拟实验结果显示，在非稳态情况下，当输入数据属于实时训化数据时，该GA-LSSVM模型的预测值与实验值的拟合度最高，拟合度R²＝1，此时GA-LSSVM模型的预测性能最优；当输入数据属于实时测试数据时，该GA-LSSVM模型的预测值与实验值的拟合度也比较高，拟合度R²＝0.9997。

(2)非稳态情况下，当利用GA-LSSVM模型预测膜蒸馏过程中膜通量变化趋势时，该模型预测的膜通量值与膜通量的实际测量值的匹配程度较高，预测的百分比偏差小于1％。

以上实施例仅为说明本发明的技术思想，不能以此限定本发明的保护范围，凡是按照本发明提出的技术思想，在技术方案基础上所做的任何改动，均落入本发明保护范围之内。

Claims (6)

1.基于GA-LSSVM模型预测膜蒸馏水处理系统中膜污染趋势的方法，其特征在于，包括以下步骤：

(1)利用LSSVM算法建立膜蒸馏污水处理过程的预测模型；

(2)分别在拟稳定状态和非稳定状态下，利用GA算法优化步骤(1)建立的预测模型的参数；

(3)分别在拟稳定状态和非稳定状态下，利用步骤(2)优化的预测模型对膜通量以及膜污染阻力的变化趋势进行预测，并分析膜蒸馏的基本操作参数对膜通量和膜污染阻力的影响；

(4)对步骤(3)得到的预测结果进行灵敏度分析计算，确定出影响膜通量和膜污染阻力的主导因素。

2.根据权利要求1所述基于GA-LSSVM模型预测膜蒸馏水处理系统中膜污染趋势的方法，其特征在于，步骤(1)的具体过程如下：

首先，输入数据记为X，输出数据记为Y，这些数据分为训练数据和测试数据，所述训练数据的集合记为A＝{(x₁,y₁),…,(x_i,y_i),…,(x_N,y_N)}，其中，x_i∈X，y_i∈Y，1≤i≤N，N为训练数据集合的大小，利用非线性映射函数建立如下式(1)的回归模型，将输入数据X映射到高维特征空间，

式(1)中，w为权值，b为偏置项，为非线性映射函数；

然后，建立式(2)所示的LSSVM优化问题，

$\min J(w,b,\mathrm{\ξ})=\frac{1}{2}{w}^{T}w+\frac{1}{2}\mathrm{\γ}\underset{i=1}{\overset{N}{\Σ}}{\mathrm{\ξ}}_i^2---\left(2\right)$

式(2)中，ξ_i为松弛因子，ξ_i＞0，γ为正则常数；

最后，构造式(3)所示的拉格朗日函数，求解LSSVM最优化问题，

式(3)中，a_i为拉格朗日乘数；分别对w,b,ξ_i,a_i求偏导数，并令偏导数为0，消去w和ξ_i，得式(4)，

$\left[\begin{array}{cc}0& 1_N^T\\ 1_N& \mathrm{\Ω}+{\mathrm{\γ}}^{-1}{I}_N\end{array}\right]\left[\begin{array}{c}b\\ a\end{array}\right]=\left[\begin{array}{c}0\\ Y\end{array}\right]---\left(4\right)$

式(4)中，Y＝[y₁,...,y_N]^T，1_N＝[1,...,1]^T，a＝[a₁,...,a_N]^T，I_N为N×N的单位矩阵，Ω为N×N的矩阵，它的第i行第j列的元素κ(x_i,x_j)是内核函数。

3.根据权利要求2所述基于GA-LSSVM模型预测膜蒸馏水处理系统中膜污染趋势的方法，其特征在于：所述内核函数为径向基函数：

$\mathrm{\κ}(x_i,x_j)=\exp \left(\frac{-||x_i-x_j^2||}{{\mathrm{\σ}}^2}\right)---\left(5\right)$

式(5)中，σ²为内核参数。

4.根据权利要求3所述基于GA-LSSVM模型预测膜蒸馏水处理系统中膜污染趋势的方法，其特征在于：步骤(2)中所述利用GA算法优化的参数包括正则常数γ和内核参数σ²。

5.根据权利要求2所述基于GA-LSSVM模型预测膜蒸馏水处理系统中膜污染趋势的方法，其特征在于，步骤(3)的具体过程为，分别在拟稳定状态和非稳定状态下进行如下操作：

将温度差、错流速度、进水pH、膜孔径和膜蒸馏处理时间作为预测模型的输入数据，将膜通量和膜污染阻力作为预测模型的输出数据，每次仅改变温度差、错流速度、进水pH、膜孔径和膜蒸馏处理时间中的一个参数，保证其他参数不变，运用GA算法优化LSSVM模型的参数，并将优化后的参数代入LSSVM模型进行预测，再利用均方根误差RMSE和回归系数R²这两个评价参数对预测性能进行评价，当RMSE越接近于0且R²越接近于1时，表明预测性能越好。

6.根据权利要求1所述基于GA-LSSVM模型预测膜蒸馏水处理系统中膜污染趋势的方法，其特征在于：在步骤(4)中，分别采用斯皮尔曼相关性、伽马相关性、肯德尔相关性和皮尔逊相关性对预测结果进行灵敏度分析计算。