CN108507924B - 一种超滤膜完整性识别方法及装置 - Google Patents
一种超滤膜完整性识别方法及装置 Download PDFInfo
- Publication number
- CN108507924B CN108507924B CN201810274603.9A CN201810274603A CN108507924B CN 108507924 B CN108507924 B CN 108507924B CN 201810274603 A CN201810274603 A CN 201810274603A CN 108507924 B CN108507924 B CN 108507924B
- Authority
- CN
- China
- Prior art keywords
- pressure
- raw water
- air chamber
- water end
- ultrafiltration membrane
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
- 239000012528 membrane Substances 0.000 title claims abstract description 165
- 238000000108 ultra-filtration Methods 0.000 title claims abstract description 141
- 238000000034 method Methods 0.000 title claims abstract description 56
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims abstract description 180
- 239000007788 liquid Substances 0.000 claims abstract description 42
- 230000008859 change Effects 0.000 claims abstract description 28
- 230000008569 process Effects 0.000 claims abstract description 27
- 230000015572 biosynthetic process Effects 0.000 claims abstract description 13
- 238000012360 testing method Methods 0.000 claims abstract description 13
- 238000005259 measurement Methods 0.000 claims abstract description 10
- 230000036962 time dependent Effects 0.000 claims abstract description 7
- 238000004364 calculation method Methods 0.000 claims description 34
- 238000012937 correction Methods 0.000 claims description 12
- 238000012935 Averaging Methods 0.000 claims description 3
- 238000005070 sampling Methods 0.000 claims description 3
- 239000007789 gas Substances 0.000 description 41
- 238000001514 detection method Methods 0.000 description 9
- IJGRMHOSHXDMSA-UHFFFAOYSA-N nitrogen Substances N#N IJGRMHOSHXDMSA-UHFFFAOYSA-N 0.000 description 8
- 238000005516 engineering process Methods 0.000 description 6
- 238000002474 experimental method Methods 0.000 description 4
- 239000012530 fluid Substances 0.000 description 4
- 229910052757 nitrogen Inorganic materials 0.000 description 4
- 238000010924 continuous production Methods 0.000 description 2
- 230000007547 defect Effects 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 239000003651 drinking water Substances 0.000 description 2
- 235000020188 drinking water Nutrition 0.000 description 2
- 239000000835 fiber Substances 0.000 description 2
- 230000006872 improvement Effects 0.000 description 2
- 238000000926 separation method Methods 0.000 description 2
- 239000010865 sewage Substances 0.000 description 2
- 241000894006 Bacteria Species 0.000 description 1
- 230000001133 acceleration Effects 0.000 description 1
- 238000001467 acupuncture Methods 0.000 description 1
- 238000005273 aeration Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000004140 cleaning Methods 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 239000000356 contaminant Substances 0.000 description 1
- 230000007812 deficiency Effects 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 229910001873 dinitrogen Inorganic materials 0.000 description 1
- 238000011066 ex-situ storage Methods 0.000 description 1
- 239000000706 filtrate Substances 0.000 description 1
- 230000036541 health Effects 0.000 description 1
- 230000002706 hydrostatic effect Effects 0.000 description 1
- 238000011065 in-situ storage Methods 0.000 description 1
- 238000009434 installation Methods 0.000 description 1
- 238000012423 maintenance Methods 0.000 description 1
- 244000005700 microbiome Species 0.000 description 1
- 238000012821 model calculation Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- QJGQUHMNIGDVPM-UHFFFAOYSA-N nitrogen group Chemical group [N] QJGQUHMNIGDVPM-UHFFFAOYSA-N 0.000 description 1
- 238000000053 physical method Methods 0.000 description 1
- 230000009467 reduction Effects 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 238000004088 simulation Methods 0.000 description 1
- 238000010998 test method Methods 0.000 description 1
- 238000009736 wetting Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N15/00—Investigating characteristics of particles; Investigating permeability, pore-volume or surface-area of porous materials
- G01N15/08—Investigating permeability, pore-volume, or surface area of porous materials
- G01N15/088—Investigating volume, surface area, size or distribution of pores; Porosimetry
-
- B—PERFORMING OPERATIONS; TRANSPORTING
- B01—PHYSICAL OR CHEMICAL PROCESSES OR APPARATUS IN GENERAL
- B01D—SEPARATION
- B01D65/00—Accessories or auxiliary operations, in general, for separation processes or apparatus using semi-permeable membranes
- B01D65/10—Testing of membranes or membrane apparatus; Detecting or repairing leaks
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
- Y02A20/00—Water conservation; Efficient water supply; Efficient water use
- Y02A20/124—Water desalination
- Y02A20/131—Reverse-osmosis
Landscapes
- Chemical & Material Sciences (AREA)
- Chemical Kinetics & Catalysis (AREA)
- Dispersion Chemistry (AREA)
- Physics & Mathematics (AREA)
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Analytical Chemistry (AREA)
- Biochemistry (AREA)
- General Health & Medical Sciences (AREA)
- General Physics & Mathematics (AREA)
- Immunology (AREA)
- Pathology (AREA)
- Separation Using Semi-Permeable Membranes (AREA)
Abstract
本发明公开了一种超滤膜完整性识别方法,其包括如下步骤:(1)将超滤膜原水端的液体排干,获得原水端气室;(2)对原水端气室进行充气,使得所述原水端气室的压力大于初始条件下气泡内的压力;(3)对充气后的原水端气室进行保压,测试所述原水端气室的压力,获得原水端气室压力随时间的变化数据;(4)建立超滤膜气泡模型,超滤膜气泡模型用于模拟超滤膜不同破损孔径对应的破损孔口处气泡形成过程中原水端气室压力随时间的变化情况;(5)将超滤膜特征参数输入超滤膜气泡模型,并将步骤(3)所述原水端气室压力随时间的变化实测数据导入超滤膜气泡模型进行计算,获得超滤膜破损尺寸,从而识别所述超滤膜的完整性。
Description
技术领域
本发明属于检测超滤膜完整性技术领域,更具体地,涉及一种超滤膜完整性识别方法及装置。
背景技术
随着生活水平的日益提升,以膜技术为核心的新型水处理工艺为污水和饮用水的治理提供了技术保障。膜分离技术对微生物、细菌和颗粒物的有效截留能力已经在水处理领域中得到了普遍的承认,其操作压力较低,出水水质相对稳定,并且可以与多种技术原位、异位耦合,可置于传统工艺的任何段位,已经具有取代传统污水和饮用水处理工艺的潜力。然而,由于膜组件在安装和维护方面的不足,剧烈的膜曝气清洗引起的过度应力变化,以及原水中锋利的污染物质等容易造成膜纤维破损。膜纤维一旦出现破损,将严重影响出水水质,对人的健康造成危害。
现有的超滤膜完整性识别的方法主要为间接测试法,比如通过检测滤出液的水质情况来间接推测超滤膜是否完整,其存在的主要缺陷为精度较低,检测不及时,操作复杂。
发明内容
针对现有技术的以上缺陷或改进需求,本发明提供了一种超滤膜完整性识别方法,其充分结合超滤膜破损时原水端气室一定压力下,在超滤膜破损孔口处产生、生长气泡的特点,提出一种超滤膜气泡模型描述超滤膜破损孔口孔径与原水端气室压力以及超滤膜本身特征参数之间的关系,通过将待测超滤膜原水端气室的实测压力数据代入该模型,计算超滤膜实际破损孔径,获知超滤膜是否完整,其为一种精度高、能够获知超滤膜破损绝对孔径大小的超滤膜完整性的直接识别方法,由此解决现有技术的超滤膜间接识别方法存在的精度较低、检测不及时或操作复杂等的技术问题。
为实现上述目的,按照本发明的一个方面,提供了一种超滤膜完整性识别方法,包括如下步骤:
(1)将超滤膜原水端的液体排干,获得原水端气室;
(2)对步骤(1)所述原水端气室进行充气,使得所述原水端气室的压力大于初始条件下超滤膜出水端破损孔口处的压力;
(3)对充气后的原水端气室进行保压,采用压力测试装置测试所述原水端气室的压力,获得所述原水端气室压力随时间的变化数据;如果该变化数据与相同条件下全新的超滤膜压降变化数据相吻合,表面该超滤膜完整;否则执行步骤(4);
(4)建立超滤膜气泡模型,所述超滤膜气泡模型用于模拟所述超滤膜不同破损孔径对应的破损孔口处气泡形成过程中原水端气室压力随时间的变化情况;
(5)将超滤膜特征参数输入所述超滤膜气泡模型,并将步骤(3)所述原水端气室压力随时间的变化数据导入所述超滤膜气泡模型进行计算,获得超滤膜破损尺寸。
优选地,步骤(2)所述原水端气室的压力为135~220kPa。
优选地,步骤(3)采用压力测试装置测试所述原水端气室的压力,每秒采样至少100个点。
优选地,步骤(4)所述超滤膜气泡模型的建立包括如下子步骤:
(4-1)建立单个气泡模型,所述单个气泡模型为:
t=0:Pg|0>Pb|0,u0=0
其中,t为时间,Pg|0为在初始时刻的气室压力,Pb|0为在初始时刻的气泡压力,Patm为大气压力,ρl为液体密度,H为膜上液体高度,σ为液体表面张力,Pe为膜弹性压力,Vb为气泡体积,u0表示通过超滤膜破损孔口时气体的流速;m为原水端气室内含有的气体的总质量,M为气室内气体的摩尔质量,R为理想气体常数,T为温度,V为体积,P为压力,Vg为气室体积,Pg为气室内的压力,Pb为气泡内的压力,ρg为气体密度,ΔP为气室内的压力Pg与气泡内的压力Pb之差,μ表示孔口系数,dOR为破损孔径,χ为温度修正系数;
(4-2)建立连续气泡模型:
将初始超滤膜特征参数输入所述单个气泡模型,所述初始超滤膜特征参数包括初始条件的原水端气室压力、膜弹性压力、液体表面张力、膜上液体高度、液体密度以及假定破损孔径,通过计算,获得原水端气室压力Pg与形成的气泡内的压力Pb随时间变化的值,比较第一个时间步长下原水端气室压力Pg与形成的气泡内的压力Pb的大小;
如果原水端气室压力Pg与气泡内的压力Pb的差大于10Pa,就认为单个气泡处于继续生长状态,没有脱离超滤膜表面,将此时刻计算所得的原水端气室压力、气泡体积、气泡内的压力和原水端气室体积的计算结果分别取代初始条件的原水端气室压力、气泡体积、气泡内的压力和原水端气室体积的值进行计算,得到同一个气泡形成过程中下一个时间步长的原水端气室压力、气泡体积、气泡内的压力和原水端气室体积的计算结果,重复该过程,不断记录循环过程中的各项数据,直到原水端气室压力Pg与气泡内的压力Pb的差小于或等于10Pa;
如果原水端气室压力Pg与气泡内的压力Pb的差小于或等于10Pa,表明气泡生长完毕,气泡脱离超滤膜表面;提取此时计算所得的原水端气室压力Pg数值,将此数值替代初始条件的原水端气室压力,其他参数仍用初始条件的各参数值,并将此数组作为初始条件,重复单个气泡连续生长状态的计算过程,直到计算所得原水端气室内压力小于或等于初始条件下气泡内部压力;
当原水端气室压力小于或等于初始条件下气泡内部压力时,气泡难以形成,则此时气泡连续形成过程结束,记录计算过程中不同时刻对应的气室压力Pg,获得连续气泡模型。
优选地,步骤(5)具体包括如下子步骤:
(5-1)将根据超滤膜气泡模型得到的某一时刻下的原水端气室压力计算数据与步骤(3)获得的同一时刻下的原水端气室压力实测数据视为一组数据,比较从气泡开始形成到气泡停止形成整个时间段中每组数据中计算数据与实测数据的比值,得到若干个比值,对该若干个比值求平均数,将该平均数作为初始假定破损孔径的修正系数;这里计算数据指的是原水端气室压力Pg计算数据,实测数据指的是原水端气室压力Pg的实测数据。
(5-2)将初始假定破损孔径与步骤(5-1)所述初始假定破损孔径的修正系数相乘后获得修正后的破损孔径,将该修正后的破损孔径作为新的初始假定破损孔径再次代入超滤膜气泡模型进行计算,其他初始条件参数值不变,得到新的初始假定破损孔径条件下的不同时刻原水端气室压力计算数据;
(5-3)返回执行步骤(5-1)和步骤(5-2),直至不同时刻下原水端气室压力的计算数据与实测数据的比值的平均值与1的差值的绝对值小于0.01,则视为计算数据与实测数据基本一致,此时计算所得的修正后的破损孔径就是该方法下返算的等效破损孔径。
按照本发明的另一个方面,提供了一种用于识别超滤膜完整性的压力检测装置,包括气体提供单元、超滤膜单元、压力传感器、数据采集单元和数据记录单元;其中,
所述超滤膜单元包括超滤膜、原水端气室和出水端液体;所述原水端气室和出水端液体分别位于所述超滤膜的两侧,所述原水端气室为将原水端液体排干后获得的气室;所述超滤膜破损后所述原水端气室的气体通过破损孔口外泄,在所述出水端液体内形成气泡;
所述气体提供单元和所述原水端气室相连通,用以向所述原水端气室提供气体以确保所述原水端气室足够的压力;
所述压力传感器与所述原水端气室相连通,用于测量所述原水端气室的压力;并将测量得到的压力信号转变为电信号;
所述数据采集单元与所述压力传感器相连接,用于采集所述压力传感器检测到的压力电信号;
所述数据记录单元与所述数据采集单元相连接,所述数据记录单元用于记录保存所述数据采集单元采集到的压力电信号数据。
优选地,所述气体提供单元为氮气瓶或空气瓶,即氮气高压瓶或空气高压瓶。
优选地,原水端气室内的压力大于初始条件下超滤膜出水端破损孔口处的压力。
优选地,所述原水端气室的压力为135~220kPa
优选地,所述数据采集单元每秒采样至少100个点。
总体而言,通过本发明所构思的以上技术方案与现有技术相比,能够取得下列有益效果:
(1)本发明提出了一种超滤膜完整性识别的方法,首先在超滤膜原水端气室进行加压充气,采用压力检测装置收集气室内部压力变化的数据,然后将该压力变化数据代入超滤膜气泡模型,根据该模型识别该超滤膜是否破损,且能够获得破损孔径大小。
(2)本发明基于连续流单个气泡形成的原理,提出了封闭气室条件下的气泡形成过程,借助于计算方法,完成封闭气室连续气泡形成的建模过程。本发明通过数学建模的方法模拟破损超滤膜在封闭气室的条件下气室内气压变化情况,为研究封闭气室内压力变化以及超滤膜完整性识别提供理论依据。
(3)本发明通过建立模型的方法监测超滤膜原水端气室压力的情况,属于非破坏性直接识别方法,检测精度高,准确度高,检测及时,且操作简便。通过超滤膜气泡模型计算方法,可以定量获得超滤膜破损孔径大小,且通过增大原水端气室初始压力可以提高检测识别精度。
附图说明
图1是本发明超滤膜完整性识别方法流程示意图;
图2是本发明用于识别超滤膜完整性的压力检测装置;
图3是本发明超滤膜单元示意图;
图4是本发明连续气泡形成流程图。
在所有附图中,相同的附图标记用来表示相同的元件或结构,其中:
1-气体提供单元;2-超滤膜单元;21-超滤膜;22-原水端气室;23-出水端液体;3-压力传感器;4-数据采集单元;5-数据记录单元,6-气泡模型。
具体实施方式
为了使本发明的目的、技术方案及优点更加清楚明白,以下结合附图及实施例,对本发明进行进一步详细说明。应当理解,此处所描述的具体实施例仅仅用以解释本发明,并不用于限定本发明。此外,下面所描述的本发明各个实施方式中所涉及到的技术特征只要彼此之间未构成冲突就可以相互组合。
本发明提供了一种超滤膜完整性识别方法,如图1所示,包括如下步骤:
(1)将超滤膜原水端的液体排干,获得原水端气室,该原水端即超滤膜工作时的进水端,在进行本发明的超滤膜完整性识别方法中,原水端液体排干后,获得封闭的原水端气室。超滤膜两侧分别为出水端液体和原水端气室;
(2)对步骤(1)所述原水端气室进行充气,使得所述原水端气室的压力大于初始条件下超滤膜出水端破损孔口处的压力。即假定超滤膜包括一个破损孔,未对原水端气室进行充气时该假定破损孔口处的压力即为初始条件下超滤膜出水端破损孔口处的压力。原水端气室压力大于超滤膜出水端破损孔口处的压力,方能在超滤膜破损孔口处向出水端方向形成气泡。但是较高的原水端气室压力可以保证破损孔径的识别精度和准确度。如果要求超滤膜破损孔口尺寸识别精度达到3微米,根据本发明的超滤膜完整性识别方法,原水端气室压力至少为135kPa,优选压力为135~220kPa,即优选压力为1.35~2.2个标准大气压。初始条件下超滤膜出水端破损孔口处的压力也可看做初始条件下气泡内的压力。
(3)对原水端气室进行保压实验,采用压力测试装置测试原水端气室压力随时间的变化,每秒采样至少100个点;获得原水端气室压力随时间的压力变化数据。如果该变化数据与相同条件下全新的超滤膜压降变化数据相吻合,表明该超滤膜完整;否则执行步骤(4)。
本发明提供了一种用于识别超滤膜完整性的压力检测装置,如图2所示,包括气体提供单元、超滤膜单元、压力传感器、数据采集单元和数据记录单元;其中,
所述超滤膜单元包括超滤膜、原水端气室和出水端液体;所述原水端气室和出水端液体分别位于所述超滤膜的两侧,所述原水端气室为将原水端液体排干后获得的气室;所述超滤膜破损后所述原水端气室的气体通过破损孔口外泄,在所述出水端液体内形成气泡。
气体提供单元和所述原水端气室相连通,用以向原水端气室提供气体以确保原水端气室足够的压力;本实施例中,气体提供单元为一个高压氮气瓶。实际应用时可以为任意提供高压氮气或空气的设备或装置。对原水端气室进行充压,使得原水端气室内的压力大于初始条件下超滤膜出水端破损孔口处的压力。优选压力为135~220kPa。
压力传感器与原水端气室相连通,用于测量原水端气室的压力;并将测量得到的压力信号转变为电信号。
数据采集单元与压力传感器相连接,用于采集压力传感器检测到的压力电信号;数据采集单元为一个数据采集卡。
数据记录单元与数据采集单元相连接,数据记录单元用于记录保存数据采集单元采集到的压力电信号数据。数据记录单元在本实施例中为一台电脑,该电脑上安装有与数据采集卡配套的数据记录软件,可以将采集卡上的数据记录保存在该软件中。数据采集卡采样率最大可为100kSa/s,即每秒采集10万个数据,采样时钟精度为20ns,分辨率12-bit,表示能分辨的数字最小变化是量程的1/4096,可以反映压力传感器所测量的数值。
(4)建立超滤膜气泡模型,所述超滤膜气泡模型用于模拟所述超滤膜破损情况下,不同破损孔径对应的破损孔口处气泡形成过程中原水端气室压力随时间的变化情况。
超滤膜气泡模型包括单个气泡模型以及连续气泡模型,其中单个气泡产生、生长模型的构建步骤具体为:
图3为超滤膜单元示意图,其包括超滤膜分割了原水端气室和出水端液体,假定超滤膜中间有一破损,破损孔径为dOR,超滤膜下端的原水端气室体积为Vg,原水端气室内气体压力为Pg,气体密度为ρg,原水端气室内含有的气体的总质量为mg。在超滤膜上方存在着滞止的液体,本发明中称之为出水端液体,液体密度ρl,液体液位高度为H,液面的压力为标准大气压Patm。气泡产生及生长过程中,气泡内的压力不断变化,记为Pb,同时气泡不断生长,气泡体积Vb不断发生增大,直到气泡脱离滤膜。气体通过超滤膜破损微孔形成气泡,过孔流速为u0。
将超滤膜下方的原水端气室、产生的气泡以及产生气泡的微小孔口(超滤膜破损微孔)视作一个控制体整体,把该整体作为研究对象,根据质量守恒定律,在任意一段时间内,进入该控制体整体内部的气体的质量,等于在该时间段内,气泡内的气体的质量变化和气室内的气体的质量变化之和(忽略气体在微小孔口及膜内部毛细结构中存留的质量变化)。根据理想状态下的气体方程PV=mRT/M,可得,
因此,气体质量的增量为原水端气室内的质量增量为由于当气室内压力达到设计压力后,关闭氮气瓶气阀,所以该控制体的气体体积增量为0,因而可得到该控制体如式(2)的质量守恒方程,
当气室内的压力Pg较气泡内部压力Pb大时,气室内的气体将通过孔口流向气泡内部,气体在滤膜破损的微小孔口内的流速为u0,单位m/s,过孔口的气体为有压流体,则气泡形成过程可视为一个薄壁小孔口的淹没出流,根据水力学公式中的薄壁小孔口淹没出流公式(3),其中Q为经过薄壁小孔口的流体出流流量,μ表示薄壁小孔口的孔口系数,Δh为薄壁两侧流体液面高度差,A为孔口面积,薄壁两侧流体压差ΔP计算公式如式(4)所示,ρ为流体密度。
ΔP=ρgΔh (4)
将超滤膜破损孔口视为水力学公式中的薄壁小孔口,超滤膜破损孔口处两侧压差即为原水端气室内的压力与气泡内的压力之差,相应地,气泡形成过程中从原水端气室进入气泡内的气体经过破损孔口的流速为
其中式(5)中μ表示超滤膜破损微孔的孔口出流系数,借鉴水力学公式中的取值,一般取0.60~0.62,本实施例中取值0.62。
气室内气体通过破损孔口的质量即为气泡的质量增量,故有,
气室内部存在的质量守恒方程如下:
其中Q入、Q出表示进、出气室的流量,在该模型中Q入为0。又因为
则有
本实验中气泡形成是一个连续不断的过程,其气室内的压力衰减也是一个连续的过程,所以认为有充足的时间来平衡内外的温度,所以该过程可以视为是一个等温的变化,因而温度修正系数χ取值为1。
气泡形成的实际情况非常复杂,实际的初始条件往往是未知的,初始气泡的形状选择是一种综合考虑后的结果。在气泡形成之初,所有物理量的速度和加速度的最初值取为0。初始条件的公式如下,其中气室的初始压力则是孔口静水压力、表面张力以及膜弹性压力Pe之和,如式(11)所示:
式中Pe为膜弹性压力,可根据现有技术测定方法对超滤膜进行实验测定,t=0表示初始时刻,Pg|0为在初始时刻的气室压力,Pb|0为在初始时刻的气泡压力,Patm为大气压力,ρl为液体密度,H为膜上液体高度,σ为湿润液体表面张力,Vb为气泡体积,u0表示通过超滤膜破损孔口时气体的流速。
当气室内部压力开始下降时,由于膜两侧压力差,气室内气流通过膜的破损处往膜另一侧移动,气泡开始形成,气泡内部的压力随着气室内部压力的降低而增大,当膜两侧达到压力平衡时,此时气泡处于脱离膜表面的阶段,气泡形成颈部收缩,然后老的气泡脱离膜表面,新的气泡进入膨胀阶段。如此周而复始,气室内压力的不断衰减将会使通过孔口的气泡形成过程不断循环往复。
最后得到的超滤膜单个气泡模型为:
根据该模型,基于四阶龙格库塔法在matlab软件里编写程序的m文件,完成单个气泡形成模型的建立。基于建立的单个气泡形成模型,用if以及for语句使之进行循环运算,达到气泡连续形成生长的状态,程序流程图如图4所示。
将初始超滤膜特征参数输入所述单个气泡模型,所述初始超滤膜特征参数包括初始条件的原水端气室压力、膜弹性压力、液体表面张力、膜上液体高度、液体密度以及假定破损孔径,通过计算,获得原水端气室压力Pg与形成的气泡内的压力Pb随时间变化的值,比较第一个时间步长下原水端气室压力Pg与形成的气泡内的压力Pb的大小;
如果原水端气室压力Pg与气泡内的压力Pb的差大于10Pa,就认为单个气泡处于继续生长状态,没有脱离超滤膜表面,将此时刻计算所得的原水端气室压力、气泡体积、气泡内的压力和原水端气室体积的计算结果分别取代初始条件的原水端气室压力、气泡体积、气泡内的压力和原水端气室体积的值进行计算,得到同一个气泡形成过程中下一个时间步长的原水端气室压力、气泡体积、气泡内的压力和原水端气室体积的计算结果,重复该过程,不断记录循环过程中的各项数据,直到原水端气室压力Pg与气泡内的压力Pb的差小于10Pa。
如果原水端气室压力Pg与气泡内的压力Pb的差小于等于10Pa,表明气泡生长完毕,气泡脱离超滤膜表面;提取此时计算所得的原水端气室压力Pg数值,将此数值替代初始条件的原水端气室压力,其他参数仍用初始条件的各参数值,并将此数组作为初始条件,重复单个气泡连续生长状态的计算过程,直到计算所得原水端气室内压力小于或等于初始条件下气泡内部压力。
当原水端气室压力小于或等于初始条件下气泡内部压力时,气泡难以形成,则此时气泡连续形成过程结束,记录计算过程中不同时刻对应的气室压力Pg,获得连续气泡模型。
(5)将超滤膜特征参数输入所述超滤膜气泡模型,并将步骤(3)所述原水端气室压力随时间的变化数据导入所述超滤膜气泡模型进行计算,获得超滤膜破损尺寸,从而识别所述超滤膜的完整性。步骤(5)具体包括如下子步骤:
(5-1)将根据超滤膜气泡模型得到的某一时刻下的原水端气室压力计算数据与步骤(3)获得的同一时刻下的原水端气室压力实测数据视为一组数据,比较从气泡开始形成到气泡停止形成整个时间段中每组数据中计算数据与实测数据的比值,得到若干个比值,对该若干个比值求平均数,将该平均数作为初始假定破损孔径的修正系数;这里计算数据指的是原水端气室压力Pg计算数据,实测数据指的是原水端气室压力Pg的实测数据。
(5-2)将初始假定破损孔径与步骤(5-1)所述初始假定破损孔径的修正系数相乘后获得修正后的破损孔径,将该修正后的破损孔径作为新的初始假定破损孔径再次代入超滤膜气泡模型进行计算,其他初始条件参数值不变,得到新的初始假定破损孔径条件下的不同时刻原水端气室压力计算数据;
(5-3)返回执行步骤(5-1)和步骤(5-2),直至不同时刻下原水端气室压力的计算数据与实测数据的比值的平均值与1的差值的绝对值小于0.01,则视为计算数据与实测数据基本一致,此时计算所得的修正后的破损孔径就是该方法下返算的等效破损孔径。
该方法能够有效识别超滤膜破损的问题,不需要对膜进行处理,操作简单,精度高。
实施例1
初始条件的参数物理量取值为:ρg=1.29kg/m3,ρl=997kg/m3,χ=1,μ=0.62,σ=72.8mN/m以及Patm=1.01×105Pa,Pg=180kPa,Pe=20mbar,H=30mm,假设破损尺寸0.24mm,Vg=365.59mL。
根据之前所建立的模型以及实验方法,获得模拟及实测Pg值,并将数值进行比较,第一次修正系数为0.708,修正后孔径0.170mm,第二次修正系数1.039,修正后孔径0.177,如此反复计算,最后修正的修正系数为1.00009,修正后孔径为0.178,实验所用针灸针直径0.180mm,经显微镜及图形处理软件处理后可知实际孔径直径为0.183mm。相对误差2.7%。
本领域的技术人员容易理解,以上所述仅为本发明的较佳实施例而已,并不用以限制本发明,凡在本发明的精神和原则之内所作的任何修改、等同替换和改进等,均应包含在本发明的保护范围之内。
Claims (4)
1.一种超滤膜完整性识别方法,其特征在于,包括如下步骤:
(1)将超滤膜原水端的液体排干,获得原水端气室;
(2)对步骤(1)所述原水端气室进行充气,使得所述原水端气室的压力大于初始条件下超滤膜出水端破损孔口处的压力;
(3)对充气后的原水端气室进行保压,采用压力测试装置测试所述原水端气室的压力,获得所述原水端气室压力随时间的变化数据;如果该变化数据与相同条件下全新的超滤膜压降变化数据相吻合,表面该超滤膜完整;否则执行步骤(4);
(4)建立超滤膜气泡模型,所述超滤膜气泡模型用于模拟所述超滤膜不同破损孔径对应的破损孔口处气泡形成过程中原水端气室压力随时间的变化情况;
所述超滤膜气泡模型的建立包括如下子步骤:
(4-1)建立单个气泡模型,所述单个气泡模型为:
t=0:Pg|0>Pb|0,u0=0
其中,t为时间,Pg|0为在初始时刻的气室压力,Pb|0为在初始时刻的气泡压力,Patm为大气压力,ρl为液体密度,H为膜上液体高度,σ为液体表面张力,Pe为膜弹性压力,Vb为气泡体积,u0表示通过超滤膜破损孔口时气体的流速;m为原水端气室内含有的气体的总质量,M为气室内气体的摩尔质量,R为理想气体常数,T为温度,V为体积,P为压力,Vg为气室体积,Pg为气室内的压力,Pb为气泡内的压力,ρg为气体密度,ΔP为气室内的压力Pg与气泡内的压力Pb之差,μ表示孔口系数,dOR为破损孔径,χ为温度修正系数;
(4-2)建立连续气泡模型:
将初始超滤膜特征参数输入所述单个气泡模型,所述初始超滤膜特征参数包括初始条件的原水端气室压力、膜弹性压力、液体表面张力、膜上液体高度、液体密度以及假定破损孔径,通过计算,获得原水端气室压力Pg与形成的气泡内的压力Pb随时间变化的值,比较第一个时间步长下原水端气室压力Pg与形成的气泡内的压力Pb的大小;
如果原水端气室压力Pg与气泡内的压力Pb的差大于10Pa,就认为单个气泡处于继续生长状态,没有脱离超滤膜表面,将此时刻计算所得的原水端气室压力、气泡体积、气泡内的压力和原水端气室体积的计算结果分别取代初始条件的原水端气室压力、气泡体积、气泡内的压力和原水端气室体积的值进行计算,得到同一个气泡形成过程中下一个时间步长的原水端气室压力、气泡体积、气泡内的压力和原水端气室体积的计算结果,重复该过程,不断记录循环过程中的各项数据,直到原水端气室压力Pg与气泡内的压力Pb的差小于或等于10Pa;
如果原水端气室压力Pg与气泡内的压力Pb的差小于或等于10Pa,表明气泡生长完毕,气泡脱离超滤膜表面;提取此时计算所得的原水端气室压力Pg数值,将此数值替代初始条件的原水端气室压力,其他参数仍用初始条件的各参数值,并将此数组作为初始条件,重复单个气泡连续生长状态的计算过程,直到计算所得原水端气室内压力小于或等于初始条件下气泡内部压力;
当原水端气室压力小于或等于初始条件下气泡内部压力时,气泡难以形成,则此时气泡连续形成过程结束,记录计算过程中不同时刻对应的气室压力Pg,获得连续气泡模型;
(5)将超滤膜特征参数输入所述超滤膜气泡模型,所述特征参数为运行过程中原水端气室压力、膜弹性压力、液体表面张力、膜上液体高度、液体密度以及假定破损孔径;并将步骤(3)所述原水端气室压力随时间的变化数据导入所述超滤膜气泡模型进行计算,获得超滤膜破损尺寸。
2.如权利要求1所述的识别方法,其特征在于,步骤(2)所述原水端气室的压力为135~220kPa。
3.如权利要求1所述的识别方法,其特征在于,步骤(3)采用压力测试装置测试所述原水端气室的压力,每秒采样至少100个点。
4.如权利要求1所述的识别方法,其特征在于,步骤(5)具体包括如下子步骤:
(5-1)将根据超滤膜气泡模型得到的某一时刻下的原水端气室压力计算数据与步骤(3)获得的同一时刻下的原水端气室压力实测数据视为一组数据,比较从气泡开始形成到气泡停止形成整个时间段中每组数据中计算数据与实测数据的比值,得到若干个比值,对该若干个比值求平均数,将该平均数作为初始假定破损孔径的修正系数;这里计算数据指的是原水端气室压力Pg计算数据,实测数据指的是原水端气室压力Pg的实测数据;
(5-2)将初始假定破损孔径与步骤(5-1)所述初始假定破损孔径的修正系数相乘后获得修正后的破损孔径,将该修正后的破损孔径作为新的初始假定破损孔径再次代入超滤膜气泡模型进行计算,其他初始条件参数值不变,得到新的初始假定破损孔径条件下的不同时刻原水端气室压力计算数据;
(5-3)返回执行步骤(5-1)和步骤(5-2),直至不同时刻下原水端气室压力的计算数据与实测数据的比值的平均值与1的差值的绝对值小于0.01,则视为计算数据与实测数据基本一致,此时计算所得的修正后的破损孔径就是该方法下返算的等效破损孔径。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810274603.9A CN108507924B (zh) | 2018-03-30 | 2018-03-30 | 一种超滤膜完整性识别方法及装置 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810274603.9A CN108507924B (zh) | 2018-03-30 | 2018-03-30 | 一种超滤膜完整性识别方法及装置 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN108507924A CN108507924A (zh) | 2018-09-07 |
CN108507924B true CN108507924B (zh) | 2023-08-04 |
Family
ID=63379569
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201810274603.9A Active CN108507924B (zh) | 2018-03-30 | 2018-03-30 | 一种超滤膜完整性识别方法及装置 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN108507924B (zh) |
Families Citing this family (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108939927A (zh) * | 2018-09-14 | 2018-12-07 | 楚格(武汉)水处理有限公司 | 超滤膜修复组件及方法 |
CN114405283A (zh) * | 2021-12-30 | 2022-04-29 | 武汉艾科滤膜技术有限公司 | 用于超滤膜完整性诊断的水膜发生装置 |
CN114771897B (zh) * | 2022-04-21 | 2024-07-19 | 深圳赛动智造科技有限公司 | 一种智能放射性药物稀释分装方法 |
Citations (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2004212230A (ja) * | 2002-12-27 | 2004-07-29 | Uerushii:Kk | 膜濾過装置に於ける膜損傷探知システム |
JP2006289309A (ja) * | 2005-04-14 | 2006-10-26 | Fuji Electric Systems Co Ltd | 膜ろ過装置の膜破損診断方法および装置 |
CN1971241A (zh) * | 2006-12-06 | 2007-05-30 | 东华大学 | 安全气囊织物动态透气性测试装置及其方法 |
CN101762561A (zh) * | 2009-11-10 | 2010-06-30 | 浙江天元生物药业股份有限公司 | 一种采用紫外监测法在线检测超滤膜完整性的方法 |
CN102087103A (zh) * | 2010-12-08 | 2011-06-08 | 上海斯纳普膜分离科技有限公司 | 超滤膜泡点测定仪及超滤膜平均孔径测定方法 |
CN102430340A (zh) * | 2011-11-15 | 2012-05-02 | 天津瑞普高科生物药业有限公司 | 一种超滤膜膜包的完整性检测方法 |
CN102989324A (zh) * | 2011-09-07 | 2013-03-27 | 克朗斯股份公司 | 超滤系统中的卫生的完整性测试 |
CN103278432A (zh) * | 2013-06-08 | 2013-09-04 | 南京拓鉒医药科技有限公司 | 一种滤膜完整性的检测装置及其检测方法 |
CN207937324U (zh) * | 2018-03-30 | 2018-10-02 | 华中科技大学 | 一种用于识别超滤膜完整性的压力检测装置 |
-
2018
- 2018-03-30 CN CN201810274603.9A patent/CN108507924B/zh active Active
Patent Citations (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2004212230A (ja) * | 2002-12-27 | 2004-07-29 | Uerushii:Kk | 膜濾過装置に於ける膜損傷探知システム |
JP2006289309A (ja) * | 2005-04-14 | 2006-10-26 | Fuji Electric Systems Co Ltd | 膜ろ過装置の膜破損診断方法および装置 |
CN1971241A (zh) * | 2006-12-06 | 2007-05-30 | 东华大学 | 安全气囊织物动态透气性测试装置及其方法 |
CN101762561A (zh) * | 2009-11-10 | 2010-06-30 | 浙江天元生物药业股份有限公司 | 一种采用紫外监测法在线检测超滤膜完整性的方法 |
CN102087103A (zh) * | 2010-12-08 | 2011-06-08 | 上海斯纳普膜分离科技有限公司 | 超滤膜泡点测定仪及超滤膜平均孔径测定方法 |
CN102989324A (zh) * | 2011-09-07 | 2013-03-27 | 克朗斯股份公司 | 超滤系统中的卫生的完整性测试 |
CN102430340A (zh) * | 2011-11-15 | 2012-05-02 | 天津瑞普高科生物药业有限公司 | 一种超滤膜膜包的完整性检测方法 |
CN103278432A (zh) * | 2013-06-08 | 2013-09-04 | 南京拓鉒医药科技有限公司 | 一种滤膜完整性的检测装置及其检测方法 |
CN207937324U (zh) * | 2018-03-30 | 2018-10-02 | 华中科技大学 | 一种用于识别超滤膜完整性的压力检测装置 |
Also Published As
Publication number | Publication date |
---|---|
CN108507924A (zh) | 2018-09-07 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN207937324U (zh) | 一种用于识别超滤膜完整性的压力检测装置 | |
CN108507924B (zh) | 一种超滤膜完整性识别方法及装置 | |
CN106814018B (zh) | 一种致密岩石气相相对渗透率测量装置及方法 | |
CN110160929B (zh) | 实现两种溶液连续渗透的柔性壁渗透仪及使用方法 | |
CN103967806B (zh) | 一种水泵汽蚀试验方法和装置 | |
CN201348601Y (zh) | 有压渗透仪 | |
CN111411930B (zh) | 一种致密气藏压裂液可视化动态滤失及排采模拟装置及模拟方法 | |
CN104697913A (zh) | 致密砂岩气体渗透性测试装置及测试方法 | |
CN110608953A (zh) | 三轴试验土单元中心位置孔隙水压力测试系统 | |
CN111239023A (zh) | 一种高温高压条件下岩石压缩系数的测试系统及测试方法 | |
CN101975735B (zh) | 多孔材料多场耦合渗透率测量装置及其测量方法 | |
WO2014087374A1 (pt) | Método e sonda de monitorização da fermentação alcoólica com espectroscopia uv-vis-swnir | |
CN105158512A (zh) | 声学多普勒剖面流速仪的检测方法 | |
CN112560237B (zh) | 一种燃气胶管的寿命快速测试方法及其测试系统 | |
CN201262547Y (zh) | 滤膜孔径测定仪 | |
CN210571877U (zh) | 低应力条件下岩土体渗透变形测试装置 | |
CN207610849U (zh) | 一种固定拱顶罐顶板稳定性试验检测装置 | |
CN113639844B (zh) | 基于光栅测温的污水处理厂曝气设备健康情况监测方法 | |
CN105181027B (zh) | 管内气液两相逆向流动检测装置 | |
CN211652637U (zh) | 一种空气-蒸汽给热系数测定装置 | |
CN204086066U (zh) | 平板膜扩渗性测试仪 | |
CN206248512U (zh) | 一种应力相关的粗粒土渗透变形特性测试仪 | |
CN212301548U (zh) | 一种岩石自由膨胀率测量仪 | |
CN206639473U (zh) | 一种浮力实验装置 | |
CN112808017B (zh) | 一种外压式中空纤维超滤膜沿程过滤特性与膜污染分布特征的检测装置与方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |