2011214049管式连续微滤膜法处理热电厂循环冷却排污水过滤模型

1、管式连续微滤膜法处理热电厂循环冷却排污水过滤模型研究史昱骁 20112140491.试验背景在技术及经济高速发展的现代社会，部分地区由于环保没跟上工业的快速发展的步伐， 环境污染严重，水体水质日趋恶化。 当人们的文化及物质生活水平不断提高，人们最为关心的环境问题之一将是如何保障人们用上卫生、安全、优质的饮用水。美国自来水协会的调查研究表明：膜法水处理技术将是未来最重要的水处理方法，因为到目前为止好没有一种单独的水处理工艺能与膜法水处理过程相媲美，用单一的物理方法同时达到完全消除致病微生物，显著降低有机和无机污染物以及供水系统的协调、紧凑、美观的水处理目的。在1996年于日本召开的题为 “21世

2、纪的水处理技术”专题会议上，国际水协 （IWA ）把膜法水处理技术列为 21世纪最重要的替代水处理技术。2管式连续微滤膜污染机理及模型开发2.1管式连续微滤膜污染机理分析膜污染是指料液中的悬浮颗粒、胶体粒子或溶解性大分子有机物在膜表面和膜孔内吸附沉积，造成膜孔径减小或堵塞，使膜通量下降的现象。根据发生的位置可以分为表膜堵塞（污染物吸附沉积在膜的表面，增加了水流过膜的阻力） 和膜孔堵塞（污染物在膜孔内吸附沉积，减小了膜孔径，从而降低了膜的通量）。根据造成污染的物质不同可分为无机污染、有机污染和生物污染，其中以有机污染和生物污染最为普遍，对膜通量的影响最大1。确定膜过滤过程的污染机理，对于了解过滤

3、过程膜通量的下降原因、预测膜通量、控制 和消除膜污染都具有一定的指导意义。对于大分子和胶体悬浮液膜过滤过程的污染机理，Marshall等2和Belfort等已经做过详细的综述，但对于颗粒悬浮液、尤其是无机颗粒悬浮 液通过微滤膜过程的污染机理，文献报导很少，普遍认为颗粒在膜表面的沉积为主4-5。本文选用的阻力模型是以Yeh and Cheng6以及Bellevile的研究为基础的， 将滤饼过滤模型改进以用于澄清过滤。2.2膜污染阻力的计算方法以达西（Darcy）定律为基础，可得出下列透水通量表达式:, PPJ（ 1）R叫Rm Rf）其中：J膜通量，L/（m 2 s）;P膜两侧压力差，kPa;卩一

4、一料液粘度，Pas;Rt 过滤总阻力，它可以分解为两个部分，即Rm膜本身的阻力和 Rf污染阻力，m-1。对过滤总阻力（Rt ）进行不同方式的分解，可以确定不同过程的污染物阻力分布810,本文提到的管式微滤膜过滤属于澄清过滤，不考虑污泥层和滤渣的阻力影响， 同时，我们没有办法测量膜孔堵塞造成的不可恢复的阻力影响，在试验阶段将其忽略（这部分阻力表现为膜的不可恢复衰老，根据美国陶氏化学给出的膜通量衰减值一一年衰减约为7%,故在短期的试验过程中，这部分变化是可以忽略的），故，只简单地将阻力分解为纯膜阻力（Rm）和污染阻力（Rf）两大类，而膜的污染阻力（Rf）包括了表膜堵塞阻力（Rc）和膜孔堵塞阻力（R

5、p）,前者指悬浮或吸附在膜表面的微滤颗粒物，可以通过反冲洗的方法冲刷掉污染物加 以恢复，后者则需要通过加强洗或化学等方法加以去除，当然也包括物理和化学方法都不能去除的永久性污染。2.2.1 次纯水通量和纯膜阻力：P膜的一次纯水通量（Jo）,是用干净的膜过滤清水测得的通量，用于计算纯膜阻力（Rm ）°J0Rm-0 Jo(3)其中：J。一一膜的一次纯水通量L/（m2s）;P纯水粘度，Pa s°2.2.2料液通量和污染阻力 在同样的操作条件下，对料液进行过滤，用稳定时的料液膜通量（Rf）°Ji）,计算污染阻力由（4）和（5）式可得，膜的污染阻力为:叫RmRc Rp)=(

6、Rm Rc Rp)Rf 二 RcRp 二 Rt - Ri其中：Ji料液透水通量，L/（m 2 s）,其它符号同前。2.2.3二次纯水通量和永久堵塞阻力膜的二次纯水通量（J2 ）,是用受污染并经过反冲洗清洗后（认为表膜堵塞阻力已经完 的膜过滤清水测得的通量，全恢复）用于计算膜孔堵塞阻力（Rp) °(4)(5)(6)(3)和（7）式得:(7)J2 二 %(Rm Rp)(8)(3)和（8）式得:Rc - Rt - Rm - RpP . P Jio J 2(9)其中：通过上述方法，可以计算各部分阻力的大小，从而对膜的污染机理有所认识。3实验结果与讨论通过在一次纯水实验、 的不同变化值，得到如

7、表J2膜的二次纯水通量，L/（ m 2 s）,其它符号同前。料液实验和二次纯水实验三种情况下，测定跨膜压差和透水通量1的实验结果：表1不同情况下的通量Rp0J 2一次纯水实验料液实验二次纯水实验跨膜压差AP (kPa)一次纯水通量2J0 ( L/(m 2 s)跨膜压差AP (kPa)料液透水通量2J1 ( L/(m 2 s)跨膜压差AP (kPa)二次纯水通量2J2 ( L/(m 2 s)5.000.420151.50.63324.70.4395.250.445153.80.63425.20.4645.500.468156.30.64625.50.4415.750.485157.70.6562

8、6.70.4796.000.503158.70.65728.30.499158.70.66428.50.5123.1求纯膜阻力由于保定热电厂连续微滤膜法处理回用冷却循环排污水工程已经在运行，故只能取进行彻底化学清洗后的膜过滤清水加以测定其纯水通量(此处数据取自项目运行初期的试运行阶段)，然后计算纯膜阻力(Rm)，曲线拟合情况图 1所示：.52屮量通水纯次一XJ/S 丈.L 48464442A P8455.25.45.65.866.2跨膜压差A P (kPa)图1压差与一次纯水通量曲线图由图1可知，新膜过滤清水时，压差与通量的关系为Jo = 0.0844 A,由公式(2)可知,1 3曲线斜率为：

9、 =0.0844，而温度为20C时，清水粘度 w为1.005 >10- Pa s (查化学%Rm手册即可得到不同温度下水的粘度值)，故在20C时，纯膜阻力为：Rm10.08440= 1.185 104(m4)(10)3.2求总阻力和污染阻力求出纯膜阻力之后，接着就要求出膜的污染阻力(Rf),这就需要对料液过滤的压差和通量进行曲线拟合，此数据取用化学清洗进行之前稳定过滤过程中的压差和通量值，曲线拟合情况如图2:压差与料液通量曲线图J/S 2s丈J量通液滤656463.62Ji = 0.0042 A PR = 0.907550152154156158160跨膜压差A P (kPa)图2压差与

10、料液通量曲线图由图2可以看出，过滤料液时，膜两侧压差与通量关系为：J1= 0.0042AP，可得：10.0042，温度为20C时，料液（热电厂循环冷却水）的粘度与纯水的粘度 由相Rt差不大，可近似看作为纯水，故总过滤阻力为：Rt =10.0042 J= 2.369 105(m)(11)由（10）和（11）知，膜的污染阻力为:Rf 二 Rc Rp 二 R _Rm = 2.251 105(m)(12)3.3表膜堵塞阻力和膜孔堵塞阻力知道膜的污染阻力之后，接着就要确定污染阻力（Rf）中表膜堵塞阻力（Rc）和膜孔堵塞阻力（Rp）的值。可以将过滤料液后的膜管经过反冲洗之后，将可去除的污染物质冲洗掉，然后

11、再过滤清水，即可以测得膜孔堵塞阻力（Rp）的值。此数据取自经过反冲洗之后，过滤初期的压差和通量值，曲线拟合情况如图3:由图3可得，压差与二次纯水通量的曲线关系为:J2 = 0.01781A可得：两“0178,故膜孔堵塞阻力为:Rp10.0178%= 5.590 104(m )(13)由（12）和（13）知，膜的表膜堵塞阻力为：Rc =Rf -Rp =1.692 105（m J（14）3.4结果与讨论经过理论推导并计算，将膜的各种阻力值及其在总阻力中所占的比例汇总如表2所示:表2膜组件的各种阻力值及所占比例阻力系数类型阻力值（x 104m-1）所占比例（%）纯膜阻力Rm1.1854.98表膜堵塞

12、阻力Rc16.9271.42膜孔堵塞阻力Rp5.59023.60污染阻力Rf22.5195.02总阻力Rt23.69100可以看出，该管式膜的纯膜阻力很小，只占总阻力的4.98%,主要阻力为污染阻力，它是纯膜阻力19倍多。在污染阻力里，表膜堵塞阻力占主要，它是膜孔堵塞阻力的 3倍多。封莉等11测定了运行过程中膜过滤压差的变化，考察了膜污染的发展情况，并探讨了减缓膜过滤压差上升和膜污染的几种措施。在温度为24 C时，通过清水和生活污水实验得出膜通量J （ L/（m2s）与压力差 AP（0.615.5kPa）呈线性相关。为此，可以在压差AP为0.615.5kPa之间测定纯水通量Jo和压差AP的值，

13、然后进行曲线拟合来计算。对滤饼过滤 而言，膜的总阻力集中在污染阻力上，膜的污染阻力是纯膜阻力的3.6倍，是影响膜通量的主要因素。这与计算结果是相一致的。以上数据表明：1）膜过滤过程中，主要阻力集中在污染阻力上，纯膜阻力影响较小；2）在污染阻力中，表膜阻力占主要，膜孔堵塞阻力也有较大的影响。所以，如何在过滤过程中解决表膜堵塞，降低膜孔堵塞，延长膜的使用寿命（即降低膜孔堵塞程度）是很重要的。带着这个问题，下面就管式连续微滤膜的过滤性能做更加深入的研究。4管式连续微滤膜过滤模型膜污染是膜在应用过程中最大的阻力之一，如何解决膜污染，发挥膜的最大过滤性能，是研究的难点之一。由于膜过滤过程的影响因素较多，

14、目前预测或关联型的模型实用性不理想12，对于以机械过滤为主要过滤机理的微滤膜来说，膜孔径大小及分布和膜孔密度（即单位膜面积的膜孔数）是考核膜分离功能的两个重要参数13。然而，目前所开发的膜过滤模型与膜参数联系较少，这就需要我们深入过滤的细部，对污染颗粒进行分析， 从而建立更加合理的模型，对微滤膜的过滤性能进行优化分析。4.1污染微粒受力分析通过对料液（澄清过滤时，料液中的污染物基本上都是悬浮颗粒）中的颗粒进行受力分析，可以确定一定操作条件下，悬浮颗粒的受力情况，从而确定膜污染的实际过程，再建立合理的模型进行优化分析。该管式连续微滤膜的过滤方式为全量过滤，膜组件竖向放置，在过滤过程中，料液中的悬

15、浮颗粒沿横向（垂直于膜组件的表面） 靠近膜表面，从而对膜造成污染和堵塞。由于压力的作用，料液在膜组件内基本上是没有紊动效应的，可以认为是沿横向的层流状态流动，所以对悬浮颗粒的运动分析可局限于层流内层，随液体流动的颗粒主要受以下几个力14:颗粒在料液中所受到的垂向浮力 Fl:关于浮力的计算公式， Saffmann 15-16，Vasseur and Cox17, Mclaughin 18等都对其产生的影响进行了研究， 结果相似，此处采用基于Almann and Ripperge严的计算方法，即：Fl =0.761 鳥 X “5/（ 15）其中：Tw 剪切应力，N ；x悬浮颗粒的粒径，m;.OP液

16、体密度，1000kg.m-;n料液粘度，Pas。颗粒受到的重力 Fg：重力的计算公式比较简单，即：Fg 二云（亠 - ：?） g x3（ 16）其中：Ps颗粒密度，kg.m-3;g重力加速度，m/s2。渗透流动产生的横向拉力（也可以理解为有压流产生的压力）Fx：由于颗粒的受力分析局限于层流内层，雷诺数很小，所以可由Stokes方程计算，即：(17)其中：v透水通量，m /m s或L/m s,其匕符号冋前。如果颗粒沉积到膜的表面或膜孔内，还要受另外两个力，即弹力和吸附力。颗粒沉积到膜表面之后，必然会受到膜给予的弹力Fn :此弹力的大小，可以认为与料液的压力相等。吸附力Fa :吸附力是由颗粒间的范

17、德华力和静电相互作用力引起的，其与颗粒形状、粗糙度，接触点数目等多个参数有关，此力是相当小的，现只考虑颗粒间的范德华引力,即：其中：, a方程参数。Fah x32 二 a2从图4中颗粒的受力分析可以看出,颗粒在没有沉降到膜表面之前,（18）在竖向上受到重力和浮力的作用，可能因为受力不平衡而上下移动，大部分是处于悬浮状态， 是不在竖向上产生运动的；在横向上只受到垂直于膜表面的拉力，受力是不可能平衡的，所以颗粒在Fx的作用下，迅速向膜表面靠近，然后粘附在膜表面或钻到膜孔内，使膜的污染阻力增加，影响膜的过滤性能。沉积到膜表面或堵塞在膜孔内部的颗粒，在过滤过程中进一步被压实或钻到膜孔的更深处，在压力过

18、大时，可能造成膜孔的不可恢复堵塞。究膜的过滤性能。连续微 谑 出 水方 向亠F：图4颗粒受力分析图4.2数学模型的建立根据以上受力分析，为了得到比较合适的过滤模型，需要有以下基本假设:1膜表面的薄层膜孔不是不规则的弯曲的，而是一组垂直于膜表面的直的筒状细管，其孔径的大小可用膜的平均孔径来表示，液体在该孔内的流动为层流；2在过滤过程中，膜平均孔径d和膜孔密度 N均随透水体积 V的增加而减少,且这种减少是呈线性变化的。根据第一条假设和液体力学的公式可知，水透过孔径为d2 P3xd的孔的流速可表示为：（19）其中：AP 跨膜压差，Pa; 卩一一料液粘度，Pas;Ax 膜的有效厚度，m； d膜孔平均孔

19、径，m。二d412 x(20)如果膜过滤料液，则可得膜的透水通量（J）为:其中：J膜的透水通量，L/（m2 s）;N膜孔密度，其它符号同前。 如果过滤纯水，则可得初始膜通量（J。），即:Jo 二 Uo 二号 Nod。2 “P32% x2 No皿 No2 o 128 x odo(21)其中：do初始膜平均孔径，m； No初始膜孔密度；3清水的粘度，Pa - s,其它符号同前。根据第二条假设，某时间膜的平均孔径（d）和膜孔密度（N）随累积透水体积（V）的变化而变化，且符合线性关系，则可得：其中：a与累积透水通量相关的膜孔孔径变化系数, 3与累积透水通量相关的膜孔密度变化系数, V 累积透水体积，m

20、3,其它符号同前。m-3，需实验确定;m-3，需实验确定。JJo二：P d0(1 一 ： V)412&iN°(1 - 'V)4: P doNo128%：x(1 一 : V)4 (1 - -V)将（22）代入（20）整理后，再除以（21）可得：考虑到本工程的实际过程，热电厂循环冷却排污水的粘度与纯水的粘度相差不多，即e则可以得到更加简化的模型：(24)(1 _： V)4 (1 - -：V) = (1 - : V)4 (1 - -V)4.3小结本节从管式连续微滤膜的污染机理出发，给出了管式连续微滤膜污染阻力的分类和计算方法，从计算数据看出，管式连续微滤膜的运行总阻力主要集

21、中在污染阻力上（占95.02%）,它是纯膜阻力的19倍多。在污染阻力中，表膜堵塞阻力是膜孔堵塞阻力的3倍多，占主要部分。我们可以通过反冲洗、加强反冲洗、化学清洗等方法来清除堵塞在膜表面及膜孔隙中的污染物质，使膜的污染阻力得到恢复。在膜污染机理的基础上， 讨论了管式连续微滤膜的过滤模型的建立过程，从污染颗粒的受力分析得出：该管式膜属于全量过滤， 所以污染颗粒受力是不平衡的， 不能通过受力分析 得出过滤模型。于是从液体力学的理论出发，建立了比较合理的过滤模型。5模型在工程运行中应用保定热电厂循环冷却排污水处理项目的实际运行数据，用上述模型指导运行参数的调 节，以优化微滤机组的过滤性能。通过matl

22、ab数学工具箱的拟合，对不同形式的模型进行比较和选择，最终确定形如：4； y =（ - ' x） 4 （1 -7 “X）的模型模拟效果最好。从模型处理的过程中可以发现，一个反冲洗周期60min是比较长的，数据拟合在中段和后段偏差较大，为了使过滤过程一直保持比较好的过滤性能，现将过滤周期调整为40min，下面再用模型对实际运行的数据进行模拟。保定热电厂在2007年9月底到10月初进行彻底的化学清洗之后，只对 2#和4#机组进行了现场数据的记录，所以下面只能对2#和4#机组的部分运行周期的数据进行整理并优化处理。5.12#机组运行数据的模拟由于机组在进行化学清洗之后，膜的性能得到比较大的提

23、高，但随着过滤过程的进行， 在开始一段时间内下降比较快，之后便趋于稳定，故取其化学清洗之后第24和第58个周期的数据进行模拟（该次清洗运行62个周期后才进行二次清洗）。表3运行数据的模拟（2 #机组第24周期）时间间流量累积透计算值模拟值隔tQ通量J水体积J/J0 值V残差(min)(m3/h)L/(m s)(m3)4 J/J。(1V)4 (1_ -V)0106.00.05021.00000.00001.00001.00000.0000时间间隔t（min）流量J/J0 值累积透 水体积V（m3）计算值模拟值残差Q(m3/h)通量 JL/(m s)4 J/J。(1 ：V)4(1 一'V)

24、5105.00.04970.99068.83330.99760.9911-0.006510105.00.04970.990617.58330.99760.9824-0.015215104.00.04920.981126.33330.99520.9736-0.021620105.00.04970.990635.00000.99760.9649-0.032725104.00.04920.981143.75000.99520.9562-0.039130104.00.04920.981152.41670.99520.9475-0.047835102.00.04830.962361.08330.9904

25、0.9388-0.051640100.20.04740.945369.58330.98600.9303-0.0557注:模拟值参数为 a =1.00189E-03, 3 =1.00000E-08 ;残差平方和为:0.0114。表4运行数据的模拟（2井机组第58周期）累积透时间间流量通量J水体积计算值模拟值隔tQJ/J0 值V残差（min）(m3/h)2L/(m s)(m3)4 J/J0(1 -： V)4 (1 一 -V)0122.00.05781.00000.00001.00001.00000.00005120.00.05680.983610.16670.99590.9898-0.006110

26、106.00.05020.868920.16670.96550.97980.014315102.00.04830.836129.00000.95620.97090.014720102.00.04830.836137.50000.95620.96240.006225104.00.04920.852546.00000.96090.9539-0.007030101.00.04780.827954.66670.95390.9452-0.00863598.00.04640.803363.08330.94670.9368-0.00994091.00.04310.745971.25000.92930.928

27、6-0.0007注：模拟值参数为 a = 1.00189E-03 , 3 = 1.00000E-08 ;残差平方和为：0.0007。上述数据表还不能直观地显示出第24周期和第58周期数据之间的差异，可以通过图形更加直观地显现出来，如图5所示。由图5可以看出，第 24周期的数据拟合效果明显没有第58周期的数据拟合效果好。4(20406080累积透水体积V (m 3)图5 2#机组计算数据与模拟数据曲线图5.24#机组运行数据的模拟对4#机组，随机抽取运行数据中的第23和第59两个周期进行模拟，结果如下。表5运行数据的模拟(4#机组第23周期)时间间隔t(min)流里Q(m3/h)诵量J累积透 水

28、体积V(3、(m )计算值模拟值(1 -： V)4 (1 - :V)残差2 、L/(m s)J/J0 值4 J/J。0113.000.05351.00000.00001.00001.00000.00005112.600.05330.99659.38330.99910.9883-0.0108时间间隔t(min)流量J/J0 值累积透 水体积V(3、(m )计算值模拟值残差Q(m3/h)通量JL/(m s)4 J/J°(1 ：V)4(1 一'V)10112.800.05340.998218.78330.99960.9790-0.020615107.300.05080.949627

29、.72500.98710.9702-0.016920100.100.04740.885836.06670.97020.9621-0.008125108.700.05150.961945.12500.99030.9529-0.037530102.600.04860.908053.67500.97620.9445-0.031735105.000.04970.929262.42500.98180.9358-0.04614097.000.04590.858470.50830.96260.9279-0.0347注：模拟值参数为a =：9.96371E-04, B = 9.96335E-04;残差平方和为

30、：0.0066。表6运行数据的模拟(4#机组第59周期)时间间隔t(min)流里Q(m3/h)诵量J累积透模拟值(1-： V)4 (1 一 -V)残差水体积V(m3)计算值2L/(m s)J/J0 值4 J/J00118.000.05591.00000.00001.00001.00000.00005117.000.05540.99159.83330.99790.9878-0.010110110.000.05210.932219.58330.98260.9781-0.004515105.000.04970.889828.75000.97120.9691-0.002120101.000.04780

31、.855937.50000.96190.9605-0.00132596.000.04550.813645.91670.94970.95220.002530100.000.04730.847553.91670.95950.9444-0.01513595.000.04500.805162.25000.94720.9360-0.01124093.000.04400.788170.16670.94220.9282-0.0140注：模拟值参数为 a=9.96371E-04 , B =9.96335E-04;残差平方和为:0.0007。从残差平方和可以看出，4#机组的模拟情况比2#机组要好，将上述两表的数

32、据绘制图形如下所示。从残差平方和可以看出，4#机组的模拟情况比 2#机组要好，将上述两表的数据绘制图形如下所示。从图6可以看出，运行初期的拟合情况还是没有运行后期拟合得好，但拟合的残差平方和比较小，还是可以接受的。综上可以看出，尽管将机组的运行周期由60min到40min后，运行数据与模拟值还是存在一定的偏差，但已经远远低于每周期60min时的值，模型的应用是很成功的。第23周期运行数据与模拟值曲线图3累积透水体积V （m ）第59周期运行数据与模拟值曲线图累积透水体积V（m 3）图64#机组计算数据与模拟数据曲线图5.3小结模拟的曲线在同一个周期的数据中， 部分模拟数据有比较大的偏差， 为此

33、，需要对微滤 机组的参数进行优化，即：将微滤机组的过滤周期从60min调整为40min，使膜表面的污染堵塞能及时清除，提高膜的过滤性能。经过优化处理后，再用模型对微滤机组的运行数据进行模拟，可以得到比较高的模拟精度（残差平方和在 0.00070.0114之间），只是运行初期的运行数据与模拟值之间还是有比 较明显的差距。6结论1管式连续微滤膜在过滤过程中的总阻力可分为纯膜阻力、表膜堵塞阻力和膜孔堵塞阻力三个部分，表膜堵塞阻力和膜孔堵塞阻力之和为污染阻力。污染阻力占总阻力的95%以 上，纯膜阻力只占有极小的部分；在污染阻力中，表膜堵塞阻力是膜孔堵塞阻力的 3 倍多， 占主要部分， 在运行过程中，

34、可以通过反冲洗、 加强反洗和化学清洗等方法恢复膜的污染阻 力。2 全量过滤方式下， 污染颗粒的受力是不平衡的， 不能通过解污染颗粒的受力平衡方程 式的方法来分析管式连续微滤膜的过滤情况。3 根据模型的模拟情况，对微滤机组的过滤周期进行调整(由原来的 60min 调整为40min )后，模型的拟合精度较高(残差平方和在0.00070.0114之间)，效果比较好。参考文献:1:杨宗政，顾平膜污染的数学模型研究进展J.化学通报，2005, 68 (W039):1-6.2 Marshall A D, Munro P A, and Tragardh G . The effect of protein f

35、ouling in microfiltration and ultrafiltration on permeate flux, protein retention and selectivty: A literature revieW J. Desalination . 1993, 91:65-108.3 Belfort G, Mikulasek P, Pimbley J M, Chung K Y. Diagnosis of membrane fouling using a rotation annular filter. 2. Dilute particle suspensions of k

36、noWn particle size J. Membrane Sci., 1993, 77:23.4 Baker R J, Fane A G, Fell C J D and Yoo B H. Factors affection flux in crossfloW filtration J. Desalinatilon . 1985, 53:59.5 Benkahla Y K, Ould-DrisA, Jaffrin M Y and Si-Hassen D. Cake groWth mechanism in cross-floW microfiltration of mineral suspen

37、sions J. Menbrane Sci., 1995, 98:107.6 Yeh H M and Cheng T W. Resietance-in-series for membrane ultrafiltration in holloW fibers of tube-andshell arrangement J. Sep. Sci. & Tech., 1993, 28:1341.7 Bellevile M P, Brillouet J M, Fuente B T and Moutounet M. Cross-floW microfiltration of a red Wine on a alumina membrane: inverstigations on fouling colloids. Proc. 2nd Intl. Conf. on Inorganic Membranes J, Montpellier, France .1991,477.8 Ousman M, Bennasar M. Determination of various hydraulic resistances during cross-floW filtration of starch grain suspens