一种介观膜表面污染模拟方法技术

本发明专利技术提供一种介观膜表面污染模拟方法，包括以下步骤：确定膜分离体系中膜性质、污水性质表征参数以及系统工况参数；设计膜表面反应体系的几何空间构造，确定涉及模块、计算流程与计算精度，构建模型；将所有参数输入所构建模型，运算导出膜表面污染过程数据；代入膜通量计算公式，计算通量变化情况；可视化输出。所述方法通过对膜表面空间网格化，采用动力学蒙特卡罗方法对污染物颗粒相对膜表面运动进行仿真，还原膜污染过程，从而实现对膜分离系统膜污染、膜通量等关键参数的实时模拟。针对不同性质的膜材料、运行工况可以高效的建模计算，为膜分离系统的设计提供理论基础，也可以在系统运行期间提供实时监控的手段与维护管理的依据。

【技术实现步骤摘要】
一种介观膜表面污染模拟方法
本专利技术属于污水处理模拟仿真算法
，具体涉及一种介观膜表面污染模拟方法。
技术介绍
随着膜分离技术的日渐成熟，膜技术在污水处理、饮用水深度净化、资源回收、物质分离等领域有着广阔的前景与发展动力。此外，在新型污染物去除、微藻培养固液分离等方面也日益受到关注，对水中藻类、病原微生物、细菌、消毒副产物、药物及个人护理品(PPCPs)等的去除效果良好。相对于传统污水处理与资源化技术，膜分离技术有着运行管理简便、自动化程度高、占地面积小、维护成本低、出水水质稳定等特点，尤其在微污染物去除方面有着巨大的优势。虽然微滤、超滤等膜分离技术有着诸多优点，但是运行过程中，随着过滤的进行，膜污染问题也会愈发严重。膜污染会提高系统运行成本、降低系统出水水质与膜组件寿命，对系统的稳定运行造成威胁，这也是膜分离技术发展的主要瓶颈问题之一。因此对膜污染的预测模拟也是膜分离技术发展的关键研究，精确地模拟不仅能帮助运维团队制定膜分离系统反冲洗策略，也能实时监测系统运行情况，提高管理精细度。膜污染过程难以观测且复杂，受到压差、通量等运行工况、水质水量、等一系列外界因素的影响，模拟预测难度大。现在主流上采用半经验公式，将膜污染过程简化为完全堵塞模型、标准堵塞模型、中间堵塞模型、滤饼层过滤模型四个机理进行模拟，只能对通量、压强等表观物理量进行预测，缺少对膜表面污染的物理性质与形成过程等重要特征的模拟。膜过滤现在普遍采用实验手段进行设计以及工况优化，导致实验周期长、成本高，且无法获取反应器内实时信息。因此急需开发能够对膜表面污染进行更准确、更全面预测的模型。
技术实现思路
本专利技术针对上述缺陷，提供一种可提供三维视角的膜污染过程并突破性的实现膜表面污染物形态的模拟的介观膜表面污染模拟方法。本专利技术提供如下技术方案：一种介观膜表面污染模拟方法，包括如下步骤：S1：确定膜分离体系中膜性质、污水性质表征参数以及系统工况参数；S2：设计膜表面反应体系的几何空间构造，确定涉及模块、计算流程与计算精度，构建模型；S3：将所有参数输入所构建模型，运算导出膜表面污染过程数据；S4：代入膜通量计算公式，计算通量变化情况；S5：可视化输出。进一步地，所述S1步骤通过实验法、经验法或文献法确定膜分离体系中膜性质、污水性质表征参数以及系统工况参数。进一步地，所述参数包括液体性质参数、膜性质参数和系统工况参数；所述液体性质参数包括液体比阻、液体浓度、液体粘度、污染物粒径；所述膜性质参数为膜自身阻力、膜与污染物颗粒间自由能；所述系统工况参数为过滤压力、过滤时间。进一步地，所述S2步骤的几何空间构造包括以下步骤：1)根据过滤液体中颗粒的粒径分布进行分类，在考虑计算精度的前提下将膜表面空间进行网格划分，设计颗粒在网格空间的运动规则；2)确定模型中颗粒与膜表面碰撞黏附、不同粒径颗粒出现概率；所述碰撞粘附计算公式如下：其中，所述N1为粘附在膜表面颗粒物个数；所述N2为未粘附颗粒个数；所述ΔEb为颗粒粘附能垒；所述ΔEIL为颗粒所受水力拖拽势能；所述β为渗透阻力；所述kb为玻尔兹曼常数；所述J为水通量；所述γ为颗粒粘附概率；所述T为体系温度；所述不同粒径颗粒出现概率计算公式如下：Pd＝f(d)；其中，所述d为颗粒物粒径；所述Pd为粒径为d的颗粒物出现概率；所述f(x)为颗粒物直径分布函数；3)根据计算资源与体系实际运行情况设定计算精度与停止条件。进一步地，所述S3步骤的模型主体采用动力学蒙特卡罗方法，包括以下步骤：S31：确定计算格位置与计算时间点；S32：依次代入模块进行颗粒反应判定，确定该时刻的颗粒状态；S33：依据颗粒变化更新格状态并更新系统时间；S34：若满足终止条件则进入后续运算或输出可视化结果，反之则继续重复步骤(1)至(4)。进一步地，所述S32步骤的颗粒状态为是否碰撞粘附。进一步地，所述S33步骤中的颗粒变化为颗粒空间堆积、压缩情况；所述更新格状态并更新系统时间为：是否循环全部小格，若达到循环全部小格则达到停止条件，若未达到循环全部小格则进入下一格继续重复所述步骤S31-S33。进一步地，所述S4步骤中，根据所述S3步骤计算所得膜表面污染情况，通过其滤饼性质计算滤饼层阻力，并计算膜通量J；所述滤饼性质包括颗粒性质、污染层厚度和孔隙率：Rf＝αf×mf；其中，所述mf为膜污染质量；所述Cb为膜表面附近颗粒物浓度；所述αf为污染比阻；所述ΔP为跨膜压差；所述μ为流体粘度；所述Rf为滤饼层阻力；所述Rp为膜自身阻力；所述mf为滤饼层质量；所述μ为液体粘度；所述T为体系温度。进一步地，所述S5步骤为根据S3、S4计算结果进行可视化输出，输出的数据包括可绘制膜表面污染高程图、泥饼孔隙率随时间变化曲线图和通量变化图。本专利技术的有益效果为：本专利技术提供的方法通过对膜表面空间网格化，采用动力学蒙特卡罗方法对污染物颗粒相对膜表面运动进行仿真，还原膜污染过程，从而实现对膜分离系统膜污染、膜通量等关键参数的实时模拟。针对不同性质的膜材料、运行工况可以高效的建模计算，为膜分离系统的设计提供理论基础，也可以在系统运行期间提供实时监控的手段与维护管理的依据现在对反应器的模拟仿真常采用实验手段，导致传统经验拟合法无法实现实时可视化的仿真，本方法采用机理模拟手段，可提供三维视角的膜污染过程并突破性的实现膜表面污染物形态的模拟。附图说明在下文中将基于实施例并参考附图来对本专利技术进行更详细的描述。其中：图1为本专利技术的设计优化流程图；图2为1分钟时膜表面污染沉积高程图；图3为150分钟时膜表面污染沉积高程图；图4为300分钟时膜表面污染沉积高程图；图5为600分钟时膜表面污染沉积高程图；图6为1-600分钟时膜表面污染沉积高程对比图；图7为膜表面泥饼总孔隙率变化曲线与膜表面不同位置泥饼孔隙率变化曲线；图8为膜通量变化曲线图。具体实施例方式下面将结合本专利技术实施例中的附图，对本专利技术实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本专利技术一部分实施例，而不是全部的实施例。基于本专利技术中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本专利技术保护的范围。针对某振动膜分离藻液设备，膜组件采用0.1μm孔径的平板超滤膜过滤，振幅2cm、频率5Hz。利用本专利所开发的一种介观膜表面污染模拟方法对所述反应器运行过程中膜污染与膜通量情况进行模拟预测，实时监控。包括如下步骤：S1：通过实验法、经验法、文献法等手段，确定膜分离体系中膜性质、污水性质表征参数与系统工况参数，主要参数包括有：液体比阻、液体浓度、液体粘度、膜自身阻力、膜与颗粒间自由能、过滤本文档来自技高网...

【技术保护点】
1.一种介观膜表面污染模拟方法，其特征在于，包括如下步骤：/nS1：确定膜分离体系中膜性质、污水性质表征参数以及系统工况参数；/nS2：设计膜表面反应体系的几何空间构造，确定涉及模块、计算流程与计算精度，构建模型；/nS3：将所有参数输入所构建模型，运算导出膜表面污染过程数据；/nS4：代入膜通量计算公式，计算通量变化情况；/nS5：可视化输出。/n

【技术特征摘要】
1.一种介观膜表面污染模拟方法，其特征在于，包括如下步骤：
S1：确定膜分离体系中膜性质、污水性质表征参数以及系统工况参数；
S2：设计膜表面反应体系的几何空间构造，确定涉及模块、计算流程与计算精度，构建模型；
S3：将所有参数输入所构建模型，运算导出膜表面污染过程数据；
S4：代入膜通量计算公式，计算通量变化情况；
S5：可视化输出。


2.根据权利要求1所述的一种介观膜表面污染模拟方法，其特征在于，所述S1步骤通过实验法、经验法或文献法确定膜分离体系中膜性质、污水性质表征参数以及系统工况参数。


3.根据权利要求1所述的一种介观膜表面污染模拟方法，其特征在于，所述参数包括液体性质参数、膜性质参数和系统工况参数；所述液体性质参数包括液体比阻、液体浓度、液体粘度、污染物粒径；所述膜性质参数为膜自身阻力、膜与污染物颗粒间自由能；所述系统工况参数为过滤压力、过滤时间。


4.根据权利要求1所述的一种介观膜表面污染模拟方法，其特征在于，所述S2步骤的几何空间构造包括以下步骤：
1)根据过滤液体中颗粒的粒径分布进行分类，在考虑计算精度的前提下将膜表面空间进行网格划分，设计颗粒在网格空间的运动规则；
2)确定模型中颗粒与膜表面碰撞黏附、不同粒径颗粒出现概率；
所述碰撞粘附计算公式如下：






其中，所述N1为粘附在膜表面颗粒物个数；所述N2为未粘附颗粒个数；所述ΔEb为颗粒粘附能垒；所述ΔEIL为颗粒所受水力拖拽势能；所述β为渗透阻力；所述kb为玻尔兹曼常数；所述J为水通量；所述γ为颗粒粘附概率；所述T为体系温度；
所述不同粒径颗粒出现概率计算公式如下：
Pd＝f(d)；
其中，所述d为颗粒物粒径；所述Pd为粒径为d的颗粒物出现概率；所述f(x)为颗粒物直径分布函数...

【专利技术属性】
技术研发人员：褚华强，张亚雷，周雪飞，肖绍赜，孙晶晶，
申请(专利权)人：同济大学，
类型：发明
国别省市：上海;31