一种用于预测矿井水在反渗透膜组件内部结垢的分析方法技术

本发明专利技术公开了一种用于预测矿井水在反渗透膜组件内部结垢的分析方法，以反渗透膜组件进水中两种以上的盐组份的浓度、进口流速、操作压力为输入条件，通过MATLAB进行多元回归分析，自动拟建出二次元回归方程，取其中一种工况下的一个输入条件作为变量、其他输入条件作为定值，从而可以算出根据变量的数值得到反渗透膜组件内部结垢位置距离进口的距离。本发明专利技术不仅可以预测反渗透膜组件内部的结垢位置，还可以预测进水中某一个输入条件改变对结垢位置的影响，为反渗透系统运行和化学清洗提供参考依据。

【技术实现步骤摘要】
一种用于预测矿井水在反渗透膜组件内部结垢的分析方法
本专利技术涉及反渗透膜组件领域，尤其涉及一种用于预测矿井水在反渗透膜组件内部结垢的分析方法。
技术介绍
矿井水作为煤炭开采过程中的伴生水资源，不加处理直接外排是巨大的污染源，处理后加以利用就是优质的水资源。近年来，随着生态、环境保护要求的不断提高，矿井水反渗透处理越来越广泛。矿井水由于长期埋藏于深部地层，水中溶解性盐分种类复杂、含量极端，与地表水相比性质完全不同。目前矿井水处理主要采用预处理+多级反渗透处理工艺，通常回收率可以达到95％以上，最终的浓盐水含盐量通常在80000mg/L以上，水中的盐类接近或者超过了溶解度极限；由于矿井水水质的特殊性及浓缩后极高的含盐量，极易造成反渗透膜结垢。然而针对矿井水水质的多级反渗透膜工艺目前尚无设计规范，加上矿井水水质的多变性以及无机盐成分的复杂性，比如硬度、氟以及硫化物浓度较高等，使得在实际运行中虽然出水水质基本可以得到保证，但膜元件结垢、堵塞严重，更换成本较高。因此急需研究不同水质成分与膜结垢的关系，对不同矿井水水质进行多级反渗透系统过程设计的结垢预测与分析。目前对矿井水反渗透系统设计多依赖于膜厂商提供的过程设计软件(WAVE、ROSA、IMSD、Q+等等)，其操作系统成熟，对于反渗透膜广泛应用的海水/苦咸水淡化领域进水水质数据库较完整，可以提供较完善的过程设计指导。但由于矿井水中无机物成分复杂而且水质多变，过程设计软件对于矿井水的数据库尚不完善，且均为“黑箱”设计，得不到RO反渗透膜内部的详细信息(内部盐度分布及结垢趋势分布等)，不能提供较全面的设计方案。
技术实现思路
本专利技术的目的是提供一种用于预测矿井水在反渗透膜组件内部结垢的分析方法。本专利技术的创新点在于本专利技术可以预测反渗透膜组件内部的结垢位置，还可以通过计算预测进水中某一个输入条件改变对结垢位置的影响，为反渗透系统运行和化学清洗提供参考依据。为实现上述专利技术目的，本专利技术的技术方案是：一种用于预测矿井水在反渗透膜组件内部结垢的分析方法，包括以下步骤：(1)以反渗透膜组件进水中两种以上的盐组份的浓度、进口流速、操作压力为输入条件，利用MATLAB中的Centralcompositedesigns进行试验设计，模拟不同工况条件下的反渗透运行状态，得到若干种工况组合，对每种工况组合进行CFD分析计算，得到各工况下的反渗透膜组件的内部结垢数据，通过MATLAB进行多元回归分析，自动拟建出二次元回归方程，得到反渗透膜组件内部结垢位置距离进口的距离与各工况下输入条件的二次元回归方程；(2)取其中一种工况下的一个输入条件作为变量、其他输入条件作为定值，得到反渗透膜组件内部结垢位置距离进口的距离与变量之间的方程式，画出反渗透膜组件内部结垢位置距离进口的距离随变量的变化图；从而可以算出根据变量的数值得到反渗透膜组件内部结垢位置距离进口的距离。进一步地，还包括以下步骤：以ANSYS建立反渗透膜组件的三维模型，输入反渗透膜组件进水中两种以上的盐组份的浓度、进口流速、操作压力，并对其进行CFD分析，得到反渗透膜组件内部所输入的盐组份的分布图。利用内部盐组份分布图，可以对多级反渗透膜工艺进行优化设计，通过改变某些输入条件，使反渗透内部结垢位置后移，从而减少结垢，提高反渗透系统的运行效率。本专利技术的有益效果是：1、本专利技术可以预测反渗透膜组件内部的结垢位置，还可以通过计算预测进水某一个变量改变对结垢位置的影响，为反渗透系统运行和化学清洗提供参考依据。2、本专利技术中还可以获得反渗透膜内部盐组份的分布图，利用其可以对多级反渗透膜工艺进行优化设计。附图说明图1为本专利技术预测膜元件内CaSO4发生结垢的位置与反渗透进水口之间的距离随进口CaSO4浓度的变化图。具体实施方式下面将结合附图对本专利技术实施例中的技术方案进行清楚、完整地描述。实施例1：一种用于预测矿井水在反渗透膜组件内部结垢的分析方法,以反渗透膜组件进水中两种以上的盐组份的浓度、进口流速、操作压力为输入条件，利用MATLAB中的Centralcompositedesigns进行试验设计，模拟不同工况条件下的反渗透运行状态，得到若干种工况组合，对每种工况组合进行CFD分析计算，得到各工况下的反渗透膜组件的内部结垢数据，通过MATLAB进行多元回归分析，自动拟建出二次元回归方程，得到反渗透膜组件内部结垢位置距离进口的距离与各工况下输入条件的二次元回归方程；取其中一种工况下的一个输入条件作为变量、其他输入条件作为定值，得到反渗透膜组件内部结垢位置距离进口的距离与变量之间的方程式，画出反渗透膜组件内部结垢位置距离进口的距离随变量的变化图；从而可以算出根据变量的数值得到反渗透膜组件内部结垢位置距离进口的距离。以ANSYS建立反渗透膜组件的三维模型，输入反渗透膜组件进水中两种以上的盐组份的浓度、进口流速、操作压力，并对其进行CFD分析，得到反渗透膜组件内部所输入的盐组份的分布图。实施例2：以进水水质为CaSO4、Na2SO4浓度，进口流速v、操作压力P为输入条件作为分析案例，利用MATLAB中的Centralcompositedesigns进行试验设计，具体试验和试验结果见表一：表一：DOE试验设计及结果从表一可以看出一共得到36个工况组合，因此，对每一种工况进行CFD分析计算，得到各工况下的反渗透膜组件的内部结垢数据，结合表一的试验数据通过MATLAB进行多元回归分析，构建二次元回归方程，如下所示：Y＝-1.1489+0.001446cCASO4+0.0000299cNa2SO4+0.09891P+9.7162v-0.00000003753cCASO4*cTDS-0.00002876cNa2SO4*P+0.004627cCASO4*v+0.000000448cNa2SO4*P-0.0000842cNa2SO4*v-0.09016P*v-0.000001325cCASO42-0.0000000002315cNa2SO42-0.002157P2-25.5567v2其中，Y为RO反渗透膜组件内部结垢位置距离进口的距离。以进口微溶盐浓度CaSO4为浓度变量，取进口Na2SO4浓度为19000mg/L，压力P＝19.48Bar，流速v＝0.15m/s为定值。分析初始CaSO4浓度对膜结垢产生的的影响。进口CaSO4浓度与临界结垢位置的拟合结果为：Y＝-0.000001325cCaSO42+0.0008667352cCaSO4+0.9886781372图1为膜元件内CaSO4发生结垢位置与反渗透进水口的之间的距离随进口CaSO4浓度变化图，由于整个膜元件的长度为0.97m，当膜元件内CaSO4发生结垢的位置大于0.97m时，说明在整个膜元件内CaSO4不发生结垢；而当膜元件内CaSO4发生结垢坐标等于0m时，说明当溶液刚进入膜元件时，CaSO4就发生了结垢。在进口Ca本文档来自技高网...

【技术保护点】
1.一种用于预测矿井水在反渗透膜组件内部结垢的分析方法，其特征在于，包括以下步骤：/n（1）以反渗透膜组件进水中两种以上的盐组份的浓度、进口流速、操作压力为输入条件，利用MATLAB中的Central composite designs进行试验设计，模拟不同工况条件下的反渗透运行状态，得到若干种工况组合，对每种工况组合进行CFD分析计算，得到各工况下的反渗透膜组件的内部结垢数据，通过MATLAB进行多元回归分析，自动拟建出二次元回归方程，得到反渗透膜组件内部结垢位置距离进口的距离与各工况下输入条件的二次元回归方程；/n（2）取其中一种工况下的一个输入条件作为变量、其他输入条件作为定值，得到反渗透膜组件内部结垢位置距离进口的距离与变量之间的方程式，画出反渗透膜组件内部结垢位置距离进口的距离随变量的变化图；从而可以算出根据变量的数值得到反渗透膜组件内部结垢位置距离进口的距离。/n

【技术特征摘要】
1.一种用于预测矿井水在反渗透膜组件内部结垢的分析方法，其特征在于，包括以下步骤：
（1）以反渗透膜组件进水中两种以上的盐组份的浓度、进口流速、操作压力为输入条件，利用MATLAB中的Centralcompositedesigns进行试验设计，模拟不同工况条件下的反渗透运行状态，得到若干种工况组合，对每种工况组合进行CFD分析计算，得到各工况下的反渗透膜组件的内部结垢数据，通过MATLAB进行多元回归分析，自动拟建出二次元回归方程，得到反渗透膜组件内部结垢位置距离进口的距离与各工况下输入条件的二次元回归方程；
（2）取其中...

【专利技术属性】
技术研发人员：郭中权，邹湘，毛维东，孙邃，肖艳，王远，马赛，刘雪菲，
申请(专利权)人：中煤科工集团杭州研究院有限公司，江苏新宜中澳环境技术有限公司，
类型：发明
国别省市：浙江;33