管式膜过滤过程中的流体力学研究

1、第32卷第6期2013年12月天津工业大学学报JOURNALOFTIANJINPOLYTECHNICUNIVERSITYVol.32No.6December2013管式膜过滤过程中的流体力学研究刘恩华，杨丽（天津工业大学中空纤维膜材料与膜过程省部共建国家重点实验室培育基地，天津300387）摘要：通过实验方法，研究了管式膜在过滤过程中进口流量、出口压力、压力降、膜通量等工艺参数之间的相互关系.实验结果显示，出口压力在0.10.15MPa时，管式膜透过液单位能耗较低；进口流量在10001200L/h时，管式膜单位能耗增加较慢，是较优的运行工艺参数.通过对实验结果进行计算和推导，得到膜管压力与膜管

2、长度的关系和通量与压力的关系.关键词：管式膜；过滤；流体力学；工艺优化中图分类号：TQ0288；X7031文献标志码：A文章编号：1671024X（2013）06000108StudyonhydromechanicsinprocessoftubularmembranefiltrationLIUEn-hua，YANGLi（StateKeyLaboratoryofHollowFiberMembraneMaterialandMembraneProcess，TianjinPolytechnicUniversity，Tianjin300387，China）Abstract：Therelationamon

3、ginletflux，outletpressure，pressuredropandpermeatefluxintheprocessoftubularmembranefiltrationisstudiedbyexperiment.Theresultsshowthatpermeateenergyconsumptionislowerwhenoutletpressureis0.1-0.15MPa，andpermeateenergyconsumptionincreasingtrendbecomesslowlywheninletfluxisupto1000-1200L/h.Bycalculationand

4、derivation，theformulaofpressureandlengthandtheformulaofpermeatefluxandpressureareobtained.Keywords：tubularmembrane；filtration；hydromechanics；processoptimization膜技术是一门应用技术，主要包括膜和膜材料的制备技术和膜应用技术1-6.膜和膜材料的制备技术主要是膜材料的选择和制备、膜的制备以及膜性能的表征技术，膜应用技术主要是膜在使用过程中的工艺参数、条件及方法等技术.前者能制备出高性能的分离膜，后者保证高性能的分离膜能很好的应用到实践中，只

5、有二者有机结合，才能实现膜技术的成功应用.膜技术有分离效率高、节能、可连续分离等优点，但也有易污染、受浓差极化影响、寿命较短等缺点，因此确定优化的膜应用工艺，对发挥膜技术优点、弥补缺点有着重要意义，有时甚至关系到膜技术应用的成功与否.管式膜的特点是流道宽、流速快、膜不易污染；强度高、运行压力大、寿命长；流动阻力小，可多支组件串联使用，但料液流速快，循环流量就大，运行压力高，也会增大管式膜使用过程中的能耗，增加运行成本.因此，如何确定管式膜的运行工艺参数，在保持管收稿日期：2013-07-05式膜的耐污染特性的基础上，尽可能减小流速，降低压力，减少能耗，对管式膜的实际应用至关重要.本文通过实验方

6、法，考察了管式膜进口流量、运行压力、压力降、膜通量等工艺参数之间的相互关系，并对其相互之间的影响进行了分析，在实验的基础上，对膜管内的流体动力学进行初步理论研究，为管式膜的工程应用设计提供指导.1实验部分1.1实验装置实验装置流程如图1所示.实验装置由水箱、泵、流量计、压力表、管式膜组件、换热器、温度计及阀门管路组成.其中管式膜由8支1m或3m长单芯膜组件通过弯头串联组成，由下至上分别编号为1、2、8.实验中使用的管式膜组件为单芯，膜材料为PVDF，基金项目：国家高技术研究发展计划项目（863计划）（2012AA03A602）第一作者：刘恩华（1973），男，博士，副研究员.E-mail：45

7、02215082膜管内径12mm，由天津工业大学制备.天津工业大学学报第32卷0.060.05压力降/MPa0.040.030.020.010.00压力表温度计纯水水槽换热器泵阀门阀门流量计压力表阀门流量计40080012001600-12000进口流量（/Lh）图4Fig.4管式膜组件单位长度膜管压力降与进口流量关系Relationbetweeninletfluxandunitmembranetubepressuredrop0.0140.012图1实验流程Fig.1Processofexperiment压力降/MPa1.2实验方法实验中，水从水箱中通过泵送入管式膜组件，自下而上流经若干根膜管

8、，每根膜管上设置透过液出0.0100.0080.0060.0040.00240080012001600-1口，透过液集合后流回水箱.通过调节进水和出水阀门来改变进口流量和出口压力；水箱中设有温度计和换热器，用于监测和控制流体温度；并测定膜管进出口压力及每根膜管的透过液量，如图2所示.其中压力表1、1的压力差为弯头压力降.11流量计流量计2000进口流量（/Lh）图5Fig.5弯头压力降与进口流量的关系Relationbetweenelbowpressuredropandinletflux16弯头压力降比率/%141210864400800120016002000图2压力测试点示意图Fig.2T

9、estpoint2.1结果与讨论无透过液系统的压力变化进口流量（/Lh-1）图6弯头压力降比率与进口流量的关系Fig.6Relationbetweenelbowpressuredropratioandinletflux采用4支3m长单芯管式膜组件串联，在没有透过液的情况下测试各点压力，实验结果如图3图6所示.0.60.5压力降/MPa0.40.30.20.10.00.000.050.100.150.20由图4和图5可知，不同进口流量下，膜管和弯头的压力降不同，进口流量越大，压力降越大，且随着进口流量的增加，压力降增加的幅度越大，在进口流量超过1200L/h后，压力降增加幅度明显加大.根据膜管内

10、流体的机械能守恒方程式7z1g+1u12+p1=z2g+1u22+p2+hf（1）式中：z1=z2，u1=u2，可知膜管直管段压力降p=p1-p2=hf=（l/d）u2/2，说明膜管压力降与膜管内液体流速平方成正比关系，即与膜管内流量平方成正比关系.图5和图6所示为连接弯头的压力降与进口流量的关系.由图中可以看出，随着进口流量的增大，弯600L/h1000L/h1800L/h出口压力/MPa图3系统压力降与出口压力关系Fig.3Relationbetweenpressuredropandoutletpressure第6期刘恩华，等：管式膜过滤过程中的流体力学研究3透过液单位能耗（/kWht-1

11、）头压力降增大，其占总压力降的比例减小.对于膜系统来讲，弯头的压力降是无效压力降，膜管的压力降为有效压降.进口流量在10002000L/h之间时，弯头的压力降占总压力降的比例变化不大，均在6.5%8.5%之间.2.2出口压力对膜过滤性能的影响分别采用8支1m和8支3m长单芯管式膜组件，使用纯水作为测试液体，在有透过液的情况下测试各点压力，分析出口压力对膜过滤性能的影响，结果如图7图9所示.0.60.5总压力降/MPa0.40.30.20.10.00.000.050.100.150.200.25出口压力/MPa1000L/h，1m1m2000L/h，1000L/h，3m3m2000L/h，1.5

12、1.20.90.60.30.00.000.050.100.150.200.251m，1000L/h2000L/h1m，3m，1000L/h3m，2000L/h出口压力/MPa图9Fig.9出口压力和透过液能耗关系Relationbetweenpermeateenergyconsumptionandouletpressure大，系统的透过液量增加，膜管内流量减小，压力降减小.比较1m和3m组件系统，3m组件系统的变化幅度更大，这主要是因为3m膜组件膜面积大，相同条件下膜透过液量较大，膜管内流量减小幅度更大，因此膜系统压力降随压力变化更显著.由图9可以看出，透过液单位能耗一开始随出口压力增加而减小

13、，之后随出口压力增加而增加，出口压力为0.100.15MPa时能耗最低.比较1m和3m组件系统，3m膜组件系统透过液单位能耗较低.2.3进口流量对膜过滤性能的影响分别采用8支1m和8支3m长单芯管式膜组件，使用纯水作为测试液体，在有透过液的情况下考察进口流量对膜过滤性能的影响，结果如图10图13所示.0.4总压力降/MPa0.0MPa0.1MPa图7出口压力与压力降的关系Fig.7Relationbetweenpressuredropandouletpressure500透过液量（/Lh-1）40030020010000.000.050.100.150.200.251000L/h2000L/h

14、0.30.20.10.0400800120016002000进口流量（/Lh-1）出口压力/MPa（a）1m组件24002100透过液量（/Lh-1）18001500120090060030000.000.050.10出口压力/MPa0.150.201000L/h2000L/h（a）1m组件0.50.4总压力降/MPa0.30.20.10.040080012001600-10.0MPa0.1MPa（b）3m组件图8出口压力与膜透过液量的关系Fig.8Relationbetweenmembranepermeatefluxandouletpressure20002400由图7、图8可以看出，管式膜

15、系统压力降随着出口压力的增大而减小，膜通量随出口压力增大而增大.出口压力越大，系统内总的压力水平升高，膜通量增进口流量（/Lh）（b）3m组件图10进口流量与总压力降关系Fig.10Relationbetweentotalpressuredropandinletflux4天津工业大学学报第32卷透过液单位能耗/（kWht-1）1.51.20.90.60.3040080012001600-11600透过液流量/（Lh-1）120080040004008001200-10.0MPa0.1MPa0.0MPa1m，1m，0.1MPa3m，0.0MPa3m，0.1MPa1600200020002400进

16、口流量（/Lh）进口流量（/Lh）（a）1m组件2000透过液流量/（Lh-1）16001200800400040080012001600-1图13进口流量对透过液能耗影响0.0MPa0.1MPaFig.13Influenceofinletfluxtopermeateenergyconsumption由图11可见，在保持出口压力一定的情况下，进口流量越大，进口压力越大，膜系统的透过液量越大，系统透过液量基本与进口流量成线性关系.以系统的进口流量和出口流量的平均值计为系统的平均流量.由图12可知，在相同出口压力下，系20002400进口流量（/Lh）（b）3m组件统平均流量均小于进口流量，与进口

17、流量基本呈线性关系，1m系统的平均流量较3m系统要高.由图13可知，透过液单位能耗随进口流量的增加而增加，1m膜组件系统较3m膜组件系统能耗大.2.4压力降变化情况分析分别采用8支1m和8支3m长单芯管式膜组件，使用纯水作为测试液体，在有透过液的情况下测试各点压力和透过液流量，分析压力降变化情况，结果如图14图17所示.0.5图11进口流量与总透过液量的关系Fig.11Relationbetweenpermeatefluxandinletflux2400平均流量/（Lh-1）20001600120080040040080012001600200024000.0MPa0.1MPa0.4压力降/M

18、Pa0.30.20.10.00.000.050.10进口流量（/Lh-1）系统膜管弯头（a）1m组件2400平均流量/（Lh-1）2000160012008004004008001200160020002400进口流量（/Lh-1）0.0MPa0.1MPa0.15出口压力/MPa（a）1m组件0.3系统膜管弯头压力降/MPa（b）3m组件0.2图12Fig.12进口流量与系统平均流量的关系Relationbetweenevenfluxintubeandinletflux0.1由图10可见，随着进口流量增加，系统总的压力降增大，增加幅度也随进口流量的增大而增大，变化趋势与无透过液系统一致.比较1

19、m膜组件系统和3m膜组件系统总压力降，1m系统总压力降随进口流量增加幅度较大.0.00.000.050.100.150.20出口压力/MPa（b）3m组件图14Fig.14出口压力与压力降的关系Relationbetweenpressuredropandoutletpressure第6期刘恩华，等：管式膜过滤过程中的流体力学研究510080%/率60膜管，1000L/h比弯头，1000L/h降膜管，2000L/h力40压弯头，2000L/h2000.000.050.100.150.200.250.30出口压力/MPa（a）1m组件10080%/率60膜管，2000L/h比弯头，2000L/h降

20、膜管，1000L/h力40弯头，1000L/h压2000.000.050.100.150.20出口压力/MPa（b）3m组件图15出口压力与压力降比率关系Fig.15Relationbetweenpressuredropratioandoutletpressure0.60.0MPa，系统0.50.0MPa，膜管0.0MPa，弯头aP0.1MPa，系统M0.40.1MPa，膜管/降0.30.1MPa，弯头力压0.20.10.04008001200160020002400进口流量（/Lh-1）（a）1m组件0.60.0MPa，系统0.50.0MPa，膜管0.0MPa，弯头aP0.1MPa，系统M0

21、.40.1MPa，膜管/降0.30.1MPa，弯头力压0.20.10.04008001200160020002400进口流量（/Lh-1）（b）3m组件图16进口流量对压力降的影响Fig.16Influenceofinletfluxtopressuredrop2.4.1出口压力的影响图14所示为不同出口压力对压力降的影响.可以看出，随着出口压力的增加，系统的总压力降、膜管1009080%/70率600.0MPa，膜管比0.0MPa，弯头降500.1MPa，膜管力400.1MPa，弯头压30201004008001200160020002400进口流量（/Lh-1）（a）1m组件10080%/率

22、600.0MPa，膜管比0.0MPa，弯头降0.1MPa，膜管力400.1MPa，弯头压2004008001200160020002400进口流量（/Lh-1）（b）3m组件图17进口流量对压力降比率的影响Fig.17Influenceofinletfluxtoratioofpressuredrop压力降、弯头压力降均下降；相比之下，3m膜组件系统下降趋势更显著.在系统中，膜管有透过液产生，而弯头仅起工艺连接作用，因此，可将膜管产生的压力降看作有效压力降，弯头产生的压力降看作无效压力降.图15所示为膜管压力降、弯头压力降占系统总压力降比率与出口压力的关系图.可以看出，膜管压力降占总压力降的比率

23、要远高于弯头压力降所占比率；膜管压力降比率与出口压力大小的关系不大；1m膜组件系统膜管压力降比率为76.5%81.5%，3m膜组件系统膜管压力降比率为88.4%89.7%，3m膜组件有效压力降比率较1m组件高.2.4.2进口流量的影响图16所示为不同进口流量对压力降的影响.可以看出，随着进口流量的增加，系统的总压力降、膜管压力降、弯头压力降均增加；相比之下，1m膜组件系统增加趋势和幅度稍高，主要是因为3m膜组件系统在进口流量增加时透过液量增加较大，系统总的流动阻力较小.图17所示为膜管压力降、弯头压力降比率与进口流量的关系图.可以看出，膜管压力降占总压力降的比率要远高于弯头压力降所占比率；进口

24、流量小于800L/h时，膜管压力降比率随进口流量增加而增大的幅度较大；当进口流量大于800L/h后，膜管压力降比率随进口流量增加而增大的幅度减小，只是略有增加.6天津工业大学学报第32卷1m膜组件系统膜管压力降比率为67.4%81.5%，3m膜组件系统膜管压力降比率为86%91.3%，3m膜组件有效压力降比率较高.2.5管式膜过滤过程压力及通量方程推导管式膜在应用过程中，需要按一定的工艺流程和要求进行设计、组装，制成成套设备才能完成其分离功能.前文已经通过实验方法研究了管式膜使用过程中主要参数与膜性能的关系，但如此大量的实验数据在设计中很难直接运用，因此需要研究各工艺参数之间的关系，建立工艺参

25、数间的函数关系，揭示其中的内在联系，为管式膜的应用设计提供理论和技术依据.在管式膜的应用过程中，主要的设计参数为操作压力和膜通量.本实验测定了不同条件下的压力和膜管的透过液量，通过对实验数据分析，推导压力与膜管长度的关系方程及通量与压力的关系方程.2.5.1压力与膜管长度的关系方程实验装置为多根自下而上串联连接的膜管，实验过程中料液顺序流过每根膜管，取第1根膜管进口定为长度原点，设距离原点任意长度l处的压力随长度变化的关系式为：P（l）=ab+l（2）式中：P（l）为压力对膜管长度的函数；l为膜管长度；a、b为常数.代入第1根膜管的进出口压力，求得a、b，即得压力与膜管长度的关系式.随后再代入

26、其余7根膜管的进出口压力，验证方程的准确性.2.5.2通量与压力的关系方程假定膜管外侧压力为0，则膜通量的变化只受膜管内压力影响，即为实验实测压力.以膜管进出口压力平均值为膜管的平均压力，同时将实验测得的透过液量换算为膜通量.取各膜管的平均压力和通量作图，作线性拟合，得通量与压力的关系方程为：F=AP+B（3）式中：F为通量（L/（m2h）；P为平均压力（MPa）；A为膜通量压力系数（L/（m2hMPa）；B为校正常数（L/m2h）.则膜管的透过液方程为m+3Q=0.012乙mA（P-P）+Bdl（4）式中：m为该膜管距离第1根膜管进口端的长度；P为各膜管之间的压力损失.并利用式（4）验证通量

27、方程的准确性.2.5.3压力与膜管长度关系方程的验证采用实验数据验证式（2）的准确性，选在温度40、出口压力0.1MPa、进口流量800L/h时，各膜管进出口压力的实验值与按式（2）计算的值进行比较，结果如表1所示.表1各膜管进出口压力的实验与计算值比较Tab.1Comparetestedvaluewithcaculatedvalueoftubepressure压力测试点长度/m实验值/MPa计算值/MPa误差/%10.000.2110.209-0.9523.000.1960.1970.5133.060.1960.193-1.5346.060.1780.1811.6956.160.1760.1

28、75-0.5769.160.1650.164-0.6179.260.1630.159-2.45812.260.1480.1480.00912.340.1470.144-2.041015.340.1370.135-1.461115.420.1360.131-3.681218.420.1210.1220.831318.500.1200.118-1.671421.500.1110.110-0.901521.610.1090.105-3.671624.610.0970.0970.00由表1可看出，推导出的公式计算值与实验值很接近，误差在-3.68%1.69%之间，表明推导出的关系式能够较准确的表征压力

29、与膜管长度的关系.选取温度40、出口压力0.1MPa时不同进口流量的实验数据，按式（2）计算得到a、b值，其随进口流量变化的情况如图18所示.0.9900.9850.980a0.9750.9700.9650.9606008001000120014001600进口流量（/Lh-1）（a）a值0.500.450.40b0.35a0.300.250.200.156008001000120014001600进口流量（/Lh-1）（b）ab值图18a、ab值随进口流量变化Fig.18Influenceofinletfluxtovalueofaandab由图18可以看出，a值均小于1，说明压力随着膜管的长

30、度增加而减小；进口流量越大，a值越小，压力（第6期刘恩华，等：管式膜过滤过程中的流体力学研究7随着膜管长度增加而下降的速率越快；ab值则表示l为0时的压力值，也就是膜管的进口压力，ab值随进口流量的增加而变大，即进口流量越大，膜管进口压力越大，与实验结果相吻合.选取在温度40、进口流量1000L/h时的不同出口压力实验数据，按式（2）计算得到a、b值，其随出口压力变化的情况如图19所示.0.990.98a0.970.960.950.000.050.100.150.200.25出口压力/MPa（a）a值0.450.400.35ba0.300.250.200.150.000.050.100.150

31、.200.25出口压力/MPa（b）ab值图19a、ab值随出口压力变化Fig.19Influenceofoutletpresssuretovalueofaandab由图19可以看出，在相同温度和进口流量下，随着出口压力的增大，a值变大，即出口压力越大，压力随着膜管长度增加而下降的速率越慢，在出口压力达到0.1MPa后，a值变化渐趋平缓；而ab值随出口压力的增加而变大，即出口压力越大，膜管进口压力越大，与实验结果相吻合.2.5.4通量与压力关系方程的验证通量与压力关系如式（3）所示，利用各膜管平均压力与通量进行线性拟合，得到拟合结果.取温度40、出口压力0.1MPa、进口流量800L/h和10

32、00L/h时，各膜管通量的实验值与计算值的比较结果如表2所示.由表2可以看出，计算值与实验值较为接近，误差很小，均小于5%，表明推导出的公式准确性较高，基本能表征通量与压力的关系.取温度40、出口压力0.1MPa时，不同进口流量的实验数据进行拟合，计算A、B值，结果如图20所示.表2膜管通量的实验与计算值比较Tab.2Comparetestedvaluewithcaculatedvalueoftubepermeate膜管进口流量/膜管通量（/Lm-2h-1）编号（Lm-2h-1）实验值计算值误差/372800381379-0.533800366351-4.104800

33、327324-0.925800311301-3.226800287279-2.7978002492594.0288002332382.15110004674813.0021000451440-2.4431000405400-1.2341000373359-3.7551000335326-2.6961000311296-4.8271000265264-0.38810002332403.0020001800A1600140012006008001000120014001600进口流量（/Lh-1）（a）A值与进口流量关系1009080B7060506008001000120014001600进口流

34、量（/Lh-1）（b）B值与进口流量关系图20A、B值与进口流量关系Fig.20RelationbetweenAandBvalueandinletflux由图20可知，在相同温度和出口压力下，A值随进口流量的增加而减小，即进口流量越大，膜通量随着压力增加而增大的变化速率越慢；而B值随进口流量的增加而增大，与实验结果相吻合.取温度40，进口流量1000L/h时，不同出口压8天津工业大学学报第32卷力下的实验数据进行拟合，计算A、B值，结果如图21所示.20001800160014001200A能耗最低的出口压力为0.100.15MPa；3m膜组加，件系统渗透液单位能耗低于1m膜组件系统.（2）进

35、口流量增加，系统总的压力降增大，增加幅度也随进口流量的增大而增大，1m膜组件系统总压力降随进口流量增加的幅度较3m膜组件系统更大；进口流量越大，进口压力越大，膜系统的透过液量越大，系统透过液量基本与进口流量成线性关系；渗透液单位能耗随进口流量的增加而增加，在进口流量为10001200L/h时，单位能耗增加较慢.（3）管式膜系统中弯头压力降随出口压力增加而减小，随进口流量增加而增大，有效压力降比率随出口压力变化不大，随进口流量增加稍有增加，3m膜组件系统的有效压力降比率高于1m膜组件系统.（4）通过计算和推导，得到压力与膜管长度关系式P（l）=ab+l，通量与压力的关系式F=AP+B，经过验证，误差在5%以内.参考文献：0.000.050.100.150.20出口压力/MPa（a）A值与出口压力关60300-300.000.050.100.150.20出口压力/MPa12赵宁，王启山，李思思.膜技术研究进展J.科技资讯，2010（10）：6-9