定义“工业透水率VMD”为1000 m 2单位膜面积在1天单位时间和1 kPa单位膜压差下的m 3透水体积。直接采用西安思源学院再生水厂3520天A 2/O-MBR系统每日运转的日产水量和膜压差大数据依序计算VMD。剔除五类异常点后得到71%~78%有效VMD。人为地将每25个有效VMD划分为一个组单元,计算组单元的VMD算术平均值。拟合组VMD算术平均值与组单元得直线型工业透水率衰减方程。因为表示VMD算术平均值与组单元比较发散,所以按“年”为划分界限对组单元进行第二次降维得“年单元”,并计算年单元的VMD算术平均值。拟合年VMD算术平均值与年单元得幂型工业透水率衰减方程。用前期的衰减方程预测后一年的工业透水率,用当期的衰减方程实测后一年的工业透水率。定义误差比为预测值与实测值的相对大小。计算年工业透水率衰减是直接、简单且方便。所得结果既可以预测膜寿命,又可以自我学习不断优化,还得出预测时段不应小于三年,最好选择四到五年为宜。 Defining “industrial permeability VMD” as the m 3 permeable volume of 1000 m 2 membrane area at 1 day and 1 kPa cross-membrane pressure difference. The VMD is directly calculated by using the daily water yield and cross-membrane pressure difference of 3520 days of Xi’an Siyuan A 2/O-MBR system. After excluding five anomalies, 71% to 78% of effective industrial permeability was obtained. Each 25 effective industrial permeability was artificially divided into a group unit by the membrane pool, calculating the VMD arithmetic average of the unit. A linear regressive equation can be got through the VMD arithmetic mean and the group unit. Because it indicates that the arithmetic mean value of VMD is more divergent from the group unit, the group unit is reduced into the annual unit. A power equation is fitted the annual VMD arithmetic mean and the annual unit. The industrial permeability in the following year is predicted by the previous attenuation equation, and the industrial permeability in the following year is measured by the current decay equation. The error ratio is defined as the relative size of the predicted and measured values. Calculating the annual industrial permeability attenuation is straightforward, simple, and convenient. The results can not only predict the membrane life, but also continuously optimize self-learning. It also shows that the prediction period should not be less than three years, so it is best to choose four to five years.
定义“工业透水率VMD”为1000 m2单位膜面积在1天单位时间和1 kPa单位膜压差下的m3透水体积。直接采用西安思源学院再生水厂3520天A2/O-MBR系统每日运转的日产水量和膜压差大数据依序计算VMD。剔除五类异常点后得到71%~78%有效VMD。人为地将每25个有效VMD划分为一个组单元,计算组单元的VMD算术平均值。拟合组VMD算术平均值与组单元得直线型工业透水率衰减方程。因为表示VMD算术平均值与组单元比较发散,所以按“年”为划分界限对组单元进行第二次降维得“年单元”,并计算年单元的VMD算术平均值。拟合年VMD算术平均值与年单元得幂型工业透水率衰减方程。用前期的衰减方程预测后一年的工业透水率,用当期的衰减方程实测后一年的工业透水率。定义误差比为预测值与实测值的相对大小。计算年工业透水率衰减是直接、简单且方便。所得结果既可以预测膜寿命,又可以自我学习不断优化,还得出预测时段不应小于三年,最好选择四到五年为宜。
工业透水率,有效工业透水率,组单元,年单元,膜寿命预测
—Taking Xi’an Siyuan University’s More than Ten Years Operation Data as an Example
Xuemei Zhang1, Dong Li1, Qinghua Ma1, Bo Fu2, Baolan Qin2, Jingyuan Hao1,3
1Energy and Chemical Engineering Research Center, Xi’an Siyuan University, Xi’an Shaanxi
2Regenerative Water Plant, Xi’an Siyuan University, Xi’an Shaanxi
3School of Chemical Engineering and Technology, Xi’an Jiaotong University, Xi’an Shaanxi
Received: Jun. 8th, 2022; accepted: Jul. 8th, 2022; published: Jul. 19th, 2022
Defining “industrial permeability VMD” as the m3 permeable volume of 1000 m2 membrane area at 1 day and 1 kPa cross-membrane pressure difference. The VMD is directly calculated by using the daily water yield and cross-membrane pressure difference of 3520 days of Xi’an Siyuan A2/O-MBR system. After excluding five anomalies, 71% to 78% of effective industrial permeability was obtained. Each 25 effective industrial permeability was artificially divided into a group unit by the membrane pool, calculating the VMD arithmetic average of the unit. A linear regressive equation can be got through the VMD arithmetic mean and the group unit. Because it indicates that the arithmetic mean value of VMD is more divergent from the group unit, the group unit is reduced into the annual unit. A power equation is fitted the annual VMD arithmetic mean and the annual unit. The industrial permeability in the following year is predicted by the previous attenuation equation, and the industrial permeability in the following year is measured by the current decay equation. The error ratio is defined as the relative size of the predicted and measured values. Calculating the annual industrial permeability attenuation is straightforward, simple, and convenient. The results can not only predict the membrane life, but also continuously optimize self-learning. It also shows that the prediction period should not be less than three years, so it is best to choose four to five years.
Keywords:Industrial Permeability, Effective Industrial Permeability, Group Unit, Annual Unit, Membrane Life Prediction
Copyright © 2022 by author(s) and beplay安卓登录
This work is licensed under the Creative Commons Attribution International License (CC BY 4.0).
http://creativecommons.org/licenses/by/4.0/
MBR膜的寿命预测和预防性维护是A2/O-MBR (Anaerobic-Anoxic-Aerobic Membrane Bio-Reactor)系统设备全生命周期管理中的重要内容 [
思源采用日本旭化成公司的偏氟乙烯(PVDF)中空纤维柱式膜组件 [
图1. 思源再生水厂的A2/O-MBR现工艺流程
为了膜寿命研究,选择第一期膜池从2011年10月1日开始到2021年10月31日截止的共计3651天运转数据。扣除从2020年1月8日到2020年5月19日因为疫情停运121天,实际运转天数为3520天,9.644年。
思源再生水厂第一期膜运转的9.644年中进行过一些设备改造或清洗工艺优化,其中有:
① 增加气浮池。2020年11月,思源再生水厂在调节池和转鼓式格栅之间,加入两台气浮池。气浮池内加入PAM聚丙烯酰胺和PAC聚合氯化铝。现仍在使用。
② 优化维护性清洗(EFM)工艺。
③ 优化膜组件进行化学在线清洗(CIP)的顺序。
④ 人工物理维护清理毛发。
目前对膜寿命的进行定性描述较多,如膜寿命是“在正常的使用条件下,膜或膜元件维持预定性能的时间” [
① 化学在线清洗前与后的透水率数据来预测膜寿命 [
② 膜供应商通常会用膜的最大累积氯接触值来衡量在膜使用期限内膜能够承受的累积氯接触量。Memcor公司就B30R膜组件给出的最大累积氯接触值为1百万ppm∙h。实际操作中,膜清洗只是影响膜寿命的一个重要因素,但却不是唯一的因素。用膜最大累积氯接触值来成功预测膜寿命的例子几乎没有。
③ 膜性能检测是对膜所处状态进行量化的手段,能够最直接表征膜的优劣,进而根据膜性能判断其是否能够继续使用。膜性能检测主要受到样品检测量(几根到几十根)和取样代表性的影响以及膜性能(纯水通量、机械强度、表面形态)并没有具体的标准规定其在应用中的限值,所以仅能作为评估膜寿命的一种辅助手段。这类方法的不足主要是取样的数量很少且是离线测量而非在线测量。
④ 国家标准也规定了中空纤维使用寿命评价办法 [
综上所述,运转数据的量、数据来源的合理性、以及实用性决定了膜寿命预测的成败。
作者首先借鉴国家标准“比通量(specific flux)” [
① 适合大工业规模的计算;每个思源膜池里的膜面积是1800 m2。所以思源每个膜组的VMD就是将该膜组的每天产水量除以膜压差后再除以1.8可得。VMD比透水率还有“无需温度校准”的巨大优势。
② 可以从控制室的计算机上直接获取全过程(9.644年)与实际操作有密切关联性的两个重要记录量(产水量、跨膜压差)的数据。
③ 实际操作的有密切关联性。在膜运行中,有两个外界的因素是在一年中因时因地会发生变化。第一个是室外温度,第二个是降雨量。如果年与年地比较膜运转状态,这两个外界的因素的影响会因为系统误差而抵消。
综上所述,采用VMD就为智能诊断与决策数据的量、合理性、和实用性提供了保证。
因为各类系统运行维护操作、仪器的失灵、修理以及记录数据的时间固定等客观原因,所以全过程工业透水率出现五类异常点:① 产水量与膜压差都等于零;② 产水量等于零;③ 膜压差等于零;④ 因为膜压差特别小造成工业透水率异常大;⑤ 因为产水量特别小造成工业透水率异常小。一一剔除所有的异常点后就产生了各膜池“有效工业透水率”。在剔除操作中,依序是强制性要求。第一期的三个膜池的有效VMD的个数列于表1。根据思源再生水厂第一期3个膜池的实际数据看,有效VMD的总天数约占VMD的总天数的71%~78%。
项目 | 1#膜池 | 2#膜池 | 3#膜池 |
---|---|---|---|
运转天数/天 | 3520 | 3520 | 3520 |
有效VMD/天 | 2474 | 2725 | 2652 |
组单元数/个 | 98 | 109 | 106 |
年单元数/个 | 9.644 | 9.644 | 9.644 |
表1. 第一期三个膜池的各项数据
现在进行大数据计算的分类降维,即“人为地”确定25个有效VMD划为一组,并称为“组单元”。第一期的三个膜池的组单元列于表1。按统计数学的方法,计算每个组单元内的VMD算术平均值。
V M D 算 术 平 均 值 = ∑ i = 1 25 x i 25 (1)
式中:小写英文字母xi——第i个有效工业透水率。
显而易见,每个膜池的组单元数与运转年数有特定的比例关系。如表1的数据显示1#膜池共有98个组单元,而98个组单元所代表的运转时间是9.644年。所以对1#膜池而言,98/9.644 = 10.16个组单元等于1年。而对于2#膜池,109/9.644 = 11.30个组单元等于1年。而对于3#膜池,106/9.644 = 10.99个组单元等于1年。
第二次就将组单元按“年”降维至10个“年单元”,并计算年单元算术平均值。
V M D 年 单 元 算 术 平 均 值 = ∑ i = 1 n x i n (2)
式中:小写英文字母xi——为年单元中第i个组单元的算术平均值。在方程(2)中,n代表组单元数。
1#膜池不同运转年数的VMD年算术平均值和组单元数列于表2。
运转年数 | 年算术平均值 | n/个 | 运转年数 | 年算术平均值 | n/个 |
---|---|---|---|---|---|
1 | 19.071 | 10 | 6 | 15.573 | 61 |
2 | 17.279 | 20 | 7 | 15.421 | 72 |
3 | 16.278 | 30 | 8 | 15.056 | 82 |
4 | 16.058 | 41 | 9 | 14.574 | 92 |
5 | 15.731 | 51 | 9.644 | 14.240 | 98 |
表2. 1#膜池不同运转年数的VMD年算术平均值和组单元数
根据表2的数据可以得到VMD年算术平均值与运转全过程的关系,如图2所示。
图2. VMD算术平均值与运转年数的关系
根据图2中VMD算术平均值与运转年数的拟合数学方程为:
V M D 算 术 平 均 值 = q ∗ 运 转 年 数 r (3)
式中:
q为第一年VMD算术平均值;
r为幂,从图1中的拟合方程可以看出幂为负值,也就是说,VMD年算术平均值随着运转年数的增加而减少。所以可以定义为工业透水率年衰减方程。
方程(3)在数学上要求至少有3个点。换句话说,要想预测在现实条件下MBR膜的寿命,系统最少需要运转三年。
现以1#膜池为例对如何进行MBR膜寿命的预测和检验加以说明。图3~5分别显示1#膜池三年到五年的实际运转VMD算术平均值与运转年数的衰减方程。
图3. VMD算术平均值与运转3年的关系
因为西安思源学院再生水厂的MBR系统已经安全有效地运转了10年,它所积累的10年数据都是实际数据,所以可以根据方程(3)利用前期运转数据预测后期运转数据,也可以利用当期运转数据检验前期运转数据的预测。现以1#膜池的数据为例,采用问与答的方式,说明如何计算预测值?如何计算实测值?和如何比较这些值的相对大小?
图4. VMD算术平均值与运转4年的关系
图5. VMD算术平均值与运转5年的关系
问题1:已经运转了3年,如何预测第4年的运转情况?
答:采用图3中已运转3年的衰减方程y = 19.077x−0.144,将x = 4代入即得y = 15.625。这是第4年的预测值。
问题2:第4年的实际运转情况如何?
答:采用图4中已运转4年的衰减方程y = 18.981x−0.129,将x = 4代入即得y = 15.873。这是第4年的实测值。
问题3:第4年的实测值与预测值是否接近?
答:现在定义一个比值,误差比,来衡量对于某个特定时间的预测值与实测值接近程度。
误 差 比 = | 预 测 值 − 实 测 值 实 测 值 | × 100 % (4)
通过计算得第4年的实测值与预测值的误差比为1.56%。误差比越小的时段,就代表预测值越接近实测值。可以用误差比选择合适的预测时段。
问题4:已经运转了4年,如何预测第5年的运转情况?
答:采用图4中已运转4年的衰减方程y = 18.981x−0.129,将x = 5代入即得y = 15.422。这是第5年的预测值。
问题5:第5年的实际运转情况如何?
答:采用图5中已运转5年的衰减方程y = 18.911x−0.121,将x = 5代入即得y = 15.565。这是第5年的实测值。
问题6:第5年的实测值与预测值是否接近?
答:通过计算得第5年的实测值与预测值的误差比为0.92%。
但如果已经运转了4年,要预测第4.5年的VMD,则采用类似图4中运转4年的衰减方程,将x = 4.5代入即可。表3列出根据表2数据计算出VMD算术平均值与运转年数的幂方程的相关数据。
运转年数 | q | r | 年衰减率/% | 预测下一年VMD | 实测下一年VMD | 误差比/% |
---|---|---|---|---|---|---|
3 | 19.077 | −0.144 | 4.88 | 15.625 | 15.873 | 1.56 |
4 | 18.981 | −0.129 | 4.09 | 15.422 | 15.565 | 0.92 |
5 | 18.911 | −0.121 | 3.54 | 15.225 | 15.354 | 0.84 |
6 | 18.834 | −0.114 | 3.08 |
表3. 不同运转年数的幂衰减方程和相关计算
同样可以得出预测时段不应小于三年,应该选择四到五年为宜的结论。
在计划更换新膜前,经营者最为关心的问题通常有,在确定的工作环境、水源条件、运行管理和膜产品性能等现实条件下,① 膜组件开始时的产水量?② 衰减到某一产水量所需时间(年)?和③ 衰减到某一产水量时剩余膜池运行性能如何?等。因为这些问题涉及到经济层面的资金运作以及性价比决策,而所得到膜寿命预测是以VMD的方式出现。这两者之间的换算是需要解释清楚。
从西安思源学院再生水厂的A2/O-MBR系统安全运转10年以上的数据看,尽管系统的产水量的确会逐年降低,但出水的水质从未出现过未达标的现象。满足质量要求是唯一前提,如果质量不达到回用水标准,就不能回用。
思源再生水厂膜池和别的厂一样,也采用恒通量运行模式。这种模式是通过自动调整抽吸泵的频率来保证恒通量运行。其中强制性控制是:当膜池跨膜压差达到50 kPa时,抽吸泵停止,强制该膜池进行必要清洗。所以当VMD = 15和VMD = 19时的区别就是膜压差的变化,而膜池的产水量的变化却不明显。换句话说,当VMD下降,导致跨膜压差上升,需要抽吸泵的频率上升,造成单位产水量的能耗上升,使得单位产水量的效益下降。
1) 因为膜寿命的预测和预防性维护等智能诊断、决策是再生水厂设备全生命周期管理和追求效益最重要的途径,直接采用3520天A2/O-MBR系统每日运转的日产水量和膜压差大数据依序计算VMD,从而保证回归数据的量、数据来源的合理性以及实用性。
2) 剔除五类异常点后得到有效VMD天数,1#膜池2474天,2#膜池2725天,3#膜池2652天。有效VMD的总天数约占VMD的总天数71%~78%。按膜池人为地将每25个有效VMD划分为一个组单元,计算单元的VMD算术平均值。回归VMD算术平均值与组单元得直线型工业透水率衰减方程。其截距表示刚刚开始启动时的VMD,负斜率意味着VMD是随着操作时间的延长而不断衰减。
3) 因为实测的VMD有近十年,所以可以用前期的衰减方程预测后一年的VMD,用当期的衰减方程实测后一年的VMD。定义误差比为预测值与实测值的相对大小,误差比越小的时段,就代表预测值越接近实测值,即选择预测时段越合适。并得出预测时段不应小于三年,应该选择四到五年为宜的结论。
4) 因为VMD算术平均值与组单元的散点图比较发散,所以按“年”为划分界限进行第二次降维得10个“年单元”,并计算年单元的VMD算术平均值。回归VMD算术平均值与年单元得幂型VMD衰减方程。同样可以用幂型VMD衰减方程计算预测值和实测值,如何预测值与实测值的相对大小以决定合适的预测时段。也得出与直线衰减方程相同的结论。
5) 计算年工业透水率衰减是直接、简单且方便。所得结果既可以预测膜寿命,又可以自我学习不断优化,还得出预测时段不应小于三年,最好选择四到五年为宜。
陕西省水利厅科技计划项目(2021slkj-14)。
张学梅,李 东,马青华,付 波,秦宝兰,郝静远. 用年工业透水率衰减预测MBR膜寿命——以西安思源学院超十年运转为例Prediction of MBR Membrane Lifetime by Annual Decay of Industrial Permeability—Taking Xi’an Siyuan University’s More than Ten Years Operation Data as an Example[J]. 水污染及处理, 2022, 10(03): 119-128. https://doi.org/10.12677/WPT.2022.103018