CN105928847B - 一种多相体系中颗粒浓度和粒径的在线测量方法 - Google Patents
一种多相体系中颗粒浓度和粒径的在线测量方法 Download PDFInfo
- Publication number
- CN105928847B CN105928847B CN201610245330.6A CN201610245330A CN105928847B CN 105928847 B CN105928847 B CN 105928847B CN 201610245330 A CN201610245330 A CN 201610245330A CN 105928847 B CN105928847 B CN 105928847B
- Authority
- CN
- China
- Prior art keywords
- particle
- method described
- image
- heterogeneous
- measuring instrument
- 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
- 239000002245 particle Substances 0.000 title claims abstract description 120
- 238000000034 method Methods 0.000 title claims abstract description 65
- 239000007787 solid Substances 0.000 claims abstract description 30
- 239000008187 granular material Substances 0.000 claims abstract description 10
- 238000003384 imaging method Methods 0.000 claims description 30
- 238000007654 immersion Methods 0.000 claims description 23
- 238000012545 processing Methods 0.000 claims description 12
- 229910001220 stainless steel Inorganic materials 0.000 claims description 8
- 239000010935 stainless steel Substances 0.000 claims description 8
- 239000011521 glass Substances 0.000 claims description 5
- 239000000463 material Substances 0.000 claims description 3
- 238000005538 encapsulation Methods 0.000 claims 1
- 230000010354 integration Effects 0.000 claims 1
- 238000002386 leaching Methods 0.000 claims 1
- 238000005259 measurement Methods 0.000 description 20
- 239000000523 sample Substances 0.000 description 9
- 239000007788 liquid Substances 0.000 description 5
- 238000002474 experimental method Methods 0.000 description 4
- 238000005286 illumination Methods 0.000 description 4
- 239000013307 optical fiber Substances 0.000 description 4
- 230000001360 synchronised effect Effects 0.000 description 4
- 230000003287 optical effect Effects 0.000 description 3
- 210000001367 artery Anatomy 0.000 description 2
- 239000012141 concentrate Substances 0.000 description 2
- 239000006185 dispersion Substances 0.000 description 2
- 239000000835 fiber Substances 0.000 description 2
- 238000005070 sampling Methods 0.000 description 2
- 238000012360 testing method Methods 0.000 description 2
- 230000005514 two-phase flow Effects 0.000 description 2
- 210000003462 vein Anatomy 0.000 description 2
- 208000035126 Facies Diseases 0.000 description 1
- 206010070834 Sensitisation Diseases 0.000 description 1
- 238000005273 aeration Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000005540 biological transmission Effects 0.000 description 1
- 238000001816 cooling Methods 0.000 description 1
- 230000007812 deficiency Effects 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 239000006196 drop Substances 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 230000007613 environmental effect Effects 0.000 description 1
- 239000012530 fluid Substances 0.000 description 1
- 230000006870 function Effects 0.000 description 1
- 238000001033 granulometry Methods 0.000 description 1
- 238000005272 metallurgy Methods 0.000 description 1
- 239000003595 mist Substances 0.000 description 1
- 239000011236 particulate material Substances 0.000 description 1
- 238000002360 preparation method Methods 0.000 description 1
- 238000003672 processing method Methods 0.000 description 1
- 229910052594 sapphire Inorganic materials 0.000 description 1
- 239000010980 sapphire Substances 0.000 description 1
- 238000007789 sealing Methods 0.000 description 1
- 238000004062 sedimentation Methods 0.000 description 1
- 230000008313 sensitization Effects 0.000 description 1
- 239000000126 substance Substances 0.000 description 1
- 238000012546 transfer Methods 0.000 description 1
- 230000001052 transient effect Effects 0.000 description 1
- 238000002834 transmittance Methods 0.000 description 1
- 230000003519 ventilatory effect Effects 0.000 description 1
- 230000000007 visual effect Effects 0.000 description 1
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01B—MEASURING LENGTH, THICKNESS OR SIMILAR LINEAR DIMENSIONS; MEASURING ANGLES; MEASURING AREAS; MEASURING IRREGULARITIES OF SURFACES OR CONTOURS
- G01B21/00—Measuring arrangements or details thereof, where the measuring technique is not covered by the other groups of this subclass, unspecified or not relevant
- G01B21/10—Measuring arrangements or details thereof, where the measuring technique is not covered by the other groups of this subclass, unspecified or not relevant for measuring diameters
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01B—MEASURING LENGTH, THICKNESS OR SIMILAR LINEAR DIMENSIONS; MEASURING ANGLES; MEASURING AREAS; MEASURING IRREGULARITIES OF SURFACES OR CONTOURS
- G01B11/00—Measuring arrangements characterised by the use of optical techniques
- G01B11/08—Measuring arrangements characterised by the use of optical techniques for measuring diameters
- G01B11/10—Measuring arrangements characterised by the use of optical techniques for measuring diameters of objects while moving
- G01B11/105—Measuring arrangements characterised by the use of optical techniques for measuring diameters of objects while moving using photoelectric detection means
-
- 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/02—Investigating particle size or size distribution
- G01N15/0205—Investigating particle size or size distribution by optical means, e.g. by light scattering, diffraction, holography or imaging
- G01N15/0227—Investigating particle size or size distribution by optical means, e.g. by light scattering, diffraction, holography or imaging using imaging, e.g. a projected image of suspension; using holography
-
- 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/06—Investigating concentration of particle suspensions
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/60—Analysis of geometric attributes
- G06T7/62—Analysis of geometric attributes of area, perimeter, diameter or volume
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04N—PICTORIAL COMMUNICATION, e.g. TELEVISION
- H04N23/00—Cameras or camera modules comprising electronic image sensors; Control thereof
- H04N23/50—Constructional details
- H04N23/51—Housings
-
- G01N15/075—
-
- 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
- G01N2015/0042—Investigating dispersion of solids
- G01N2015/0053—Investigating dispersion of solids in liquids, e.g. trouble
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/10—Image acquisition modality
- G06T2207/10024—Color image
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/10—Image acquisition modality
- G06T2207/10068—Endoscopic image
Abstract
本发明提供了一种多相体系中颗粒浓度和粒径的在线测量方法,所述方法基于一种在线多相测量仪,所述测量方法包括如下步骤:(1)将在线多相测量仪置于多相体系中,得到多相体系中颗粒的图像;(2)确定有效颗粒:将颗粒边界处灰度梯度Grad(Φ)≥Grad(Φl/2)的颗粒表示为有效颗粒;(3)确定单个像素的实际长度,并测得有效颗粒所占的像素个数,则颗粒粒径di为有效颗粒所占的像素数与单个像素实际长度的乘积;根据计算有效颗粒的浓度。该方法能够实时在线获得多相反应器内气泡、液滴或固体等颗粒浓度和粒径分布,并且测量的准确率较高。
Description
技术领域
本发明属于物料测量技术领域,涉及一种多相体系中颗粒浓度和粒径的在线测量方法,尤其涉及一种两相及三相反应器内分散相颗粒浓度和粒径分布的测量方法。
背景技术
在化学工业、冶金及环境等领域常常存在一些复杂多相流过程,涉及一个、两个甚至两个以上的分散相,该类反应器内分散相颗粒浓度和粒径分布的测量是迫切需要解决的难题。
目前广泛应用于两相流中测量的光纤探针法和电导探针法,当存在多个分散相时,不同分散相的产生同向的电信号,相互之间的干扰导致无法获得准确的测量结果。例如,多相管流中相含率和相界面的单丝电容探针测量系统(中国发明专利,授权号ZL200610042792.4)和基于双头电容探针的两相流参数测量方法和装置(中国发明专利,授权号ZL 101413911A)。取样法可以同时测量多相反应器内的分散相浓度和粒径分布,原理简单,操作方便,然而实际操作中很难达到准确取样需要的等动量条件,同时因其是离线测量不适用于非稳态过程。粒度仪是用物理的方法测试固体颗粒的大小和分布的一种仪器,根据测试原理的不同分为沉降式粒度仪、沉降天平、激光粒度仪、光学颗粒计数器、电阻式颗粒计数器及颗粒图像分析仪等。各种粒度仪在使用过程中往往需要先进行繁琐的制样,故粒度仪测量方法也不是真正的在线测量方法。
光学拍照法是一种最直观的测量方法,非侵入式照相法虽然具有不干扰多相流场的优点,但其要求被测装置透明或者安装有视窗,同时被测气泡或液滴浓度要求较稀,只能测量有限厚度反应器或壁面附近的流动。侵入式照相法虽然会对反应器内流场产生一定的干扰,但是该方法对反应器没有特殊要求。其中比较典型的有利用光纤传相原理的光纤内窥镜法,但该方法受限于光纤直径导致分辨率较低,同时由于光纤束排丝困难导致其有效通光直径往往较小,为了进行有效测量前端需要安装角度较大的广角物镜,造成了较大的图像畸变。
发明内容
针对现有技术存在的不足,本发明的目的在于提供一种多相体系中颗粒浓度和粒径的在线测量方法,能够实时在线获得多相反应器内气泡、液滴或固体等颗粒浓度和粒径分布。
为达此目的,本发明采用以下技术方案:
一种多相体系中颗粒浓度和粒径的在线测量方法,所述方法基于一种浸入式在线多相测量仪,所述浸入式在线多相测量仪包括:
封装管;
视窗,密封连接于封装管一端;
用于照明多相流的照明系统,包括LED灯和与LED灯相连的可调制光源,可调制光源包括电源、信号发生器和示波器;
用于拍照的照相系统,包括远心镜头和图像传感器;
与信号发生器和图像传感器相连的控制器;
与图像传感器相连的信号处理与输出系统;
与信号处理与输出系统相连的显示系统;
其中,LED灯、远心镜头和图像传感器位于封装管中,控制器控制图像传感器的曝光周期小于信号发生器的脉冲周期;
所述测量方法包括如下步骤:
(1)将浸入式在线多相测量仪置于多相体系中,得到多相体系中颗粒的图像;
(2)确定有效颗粒:首先,确定远心镜头的焦平面位置;之后,将被测物置于封装管前方、焦平面前后l/2处,其中,l为远心镜头的景深,单位为mm,采用在线多相测量仪对被测物进行拍照,得到被测物图像,识别读取图像中被测物边界灰度梯度值Grad(Φl/2),Φl/2为焦平面前后两侧l/2平面处灰度值,将颗粒边界处灰度梯度Grad(Φ)≥Grad(Φl/2)的颗粒表示为有效颗粒;
(3)确定单个像素的实际长度,并测得有效颗粒的像素个数,则有效颗粒粒径di为有效颗粒的像素数与单个像素实际长度的乘积;根据计算得到有效颗粒的浓度,其中,S为图像传感器感光部件的有效面积,单位为mm2;l为远心镜头的景深,单位为mm,di为有效颗粒粒径,单位为mm;Vc为有效颗粒的总体积,V为测量区域的体积;n为有效颗粒的个数。
S值由所选图像传感器型号确定,l值则由远心镜头型号确定。
本发明提供的测量方法基于远心照相的浸入式在线多相测量仪,所述在线多相测量仪具有独特的平行光路,其获得的图像几乎不畸变,具有超长物距能够伸入反应器内各个位置进行拍照测量。通过设计适当的操作程序结合图像处理方法,能够实现高分散相浓度的两相甚至三相反应器中分散相浓度和粒径分布的实时在线测量。
步骤(1)所述多相体系中的颗粒为浅色不透明固体颗粒时,光线在其表面的反射和周围介质存在巨大差异,图像拍摄时只使用多相测量仪同向照明、控制信号发生器的脉冲周期和图像传感器的曝光周期同步即可得到高锐度的图像。
所述常规的信号发生器的脉冲周期和图像传感器的曝光周期同步控制步骤主要包括:
(a)根据所测体系和周围环境情况,调节在线多相测量仪的可调制光源的电源的输出电流,通过信号发生器和示波器对脉冲信号的波长和周期进行显示;
(b)通过控制器设置图像拍摄的曝光时间、光平衡、帧频和增益;
(c)通过控制器设置照明信号的脉冲周期大于图像传感器的曝光时间,实现脉冲光和图像拍摄的同步。
由于光线在深色固体颗粒、气泡或液滴表面的反射和周围介质差距较小,只通过浅色颗粒常规的调节难以实现高锐度图像的拍摄。因此,步骤(2)所述多相体系中的颗粒为深色固体颗粒、气泡或液滴时,除进行常规的信号发生器的脉冲周期和图像传感器的曝光周期同步控制外,在线多相测量仪拍照时需加装反光板,被测颗粒置于反光板和在线多相测量仪前端面之间。
优选地,所述反光板与在线多相测量仪前端面的距离大于远心镜头的景深,如是远心镜头景深的1.2倍、1.4倍、1.6倍、1.8倍、2.0倍或2.5倍等。
优选地,所述定向反光板与多相测量仪前端面的距离为1.5倍的远心镜头景深。
优选地,所述反光板的尺寸与图像传感器的有效感光区域相同。
优选地,所述反光板为中心凸起的方形板。
优选地,所述反光板表面为亮白色或银白色。
步骤(2)确定焦平面位置的方法为:将一个精确度至少为0.1mm的刻度尺(如精确度为0.05mm、0.02mm或0.01mm等)置于与所测多相体系相同的介质中,采用在线多相测量仪对所述介质拍照,调整刻度尺到浸入式在线测量仪前端的距离,观察照片中刻度尺的清晰程度,刻度尺最清晰的位置为焦平面位置。
步骤(3)所述确定单个像素的实际长度具体为:将一个精确度至少为0.1mm的刻度尺置于在线多相测量仪前方,采用图像传感器驱动控制软件拍摄图像,通过图像处理软件拾取刻度尺上10mm距离所对应的像素个数N10,从而确定单个像素的实际长度。所述在线多相测量仪前方是指视窗前方。所述图像处理软件也可拾取刻度尺上20mm、30mm、40mm、45mm或50mm距离所对应的像素数。
优选地,所述图像处理软件为Image-Pro Plus;
优选地,步骤(3)所述有效颗粒的粒径为di=10*ni/N10,其中,di为颗粒尺寸,mm;ni为图像中所测有效颗粒对应的像素个数;N10为刻度尺上10mm长度对应的像素数,“*”表示相乘的意思。
为了测量时探头能深入到多相反应器的各个位置,所述远心镜头的工作距离为250~550mm,如260mm、300mm、350mm、380mm、420mm、470mm或520mm等,景深为1~3.7mm,如1.2mm、1.5mm、1.8mm、2.0mm、2.2mm、2.5mm、2.8mm、3.0mm或3.5mm等。
为了尽可能小的干扰流场,可适当放弃镜头的放大倍率。优选地,所述远心镜头的放大倍率为0.5~1倍,如0.6倍、0.7倍、0.8倍或0.9倍等。
优选地,所述远心镜头的外径为19~25mm,如20mm、21mm、22mm、23mm或24mm等。
所述图像传感器为CCD相机或CMOS相机。
优选地,所述CCD相机或CMOS相机的曝光时间≤1ms,如0.1ms,0.5ms,1ms等,分辨率为5~15μm,如6μm、7μm、8μm、9μm、10μm、12μm或14μm等,长和宽方向像素数至少为800(H)×600(V),如2560(H)×1920(V)、2048(H)×1536(V)、1600(H)×1200(V)、1280(H)×1024(V)或800(H)×600(V)等,帧频至少为60fps,如60fps,100fps,150fps,200fps或1000fps等。
所述LED灯位于封装管中,所述可调制光源位于封装管外。通过可调制光源能够调制不同多相体系照相所需要的不同波长的脉冲光。为了照明更均匀,所述LED灯的个数至少为12个,如12个、16个、20个或24个等。
优选地,所述LED灯在封装管中组成环形均匀排列。在保证亮度的前提下环形LED灯的内径应尽可能小。
优选地,所述LED灯与可调制光源之间通过导线连接。
所述封装管由直径不同的前段管和后段管组成;
优选地,所述前段管的外径为25~30mm,如25mm、26mm或28mm等,长度为300~600mm,如320mm、350mm、400mm、450mm、500mm、550mm或580mm等。本领域技术人员可根据所选远心镜头的具体参数决定封装管前段管的具体尺寸。
优选地,所述后段管的外径为50mm,长度为50mm,本领域技术人员可根据所选图像传感器的尺寸决定后段封装管的具体尺寸。
优选地,所述封装管的材质为不锈钢。
所述视窗、LED灯和远心镜头封装于前段管中,所述视窗置于前段管远离后段管的一端,其后依次设置LED灯和远心镜头,所述图像传感器(CCD或CMOS相机)封装于后段管中。
所述视窗为内侧镀增透膜的圆形玻璃,以使光线透射率超过95%,视窗与封装管之间做防水机械密封。
所述图像传感器与控制器通过高速数据线连接,实现图像的高速传输。
所述显示系统为显示屏。所述显示系统用于显示从信号处理与输出系统接收到的信号。
优选地,所述信号处理与输出系统、控制器及显示系统集成为计算机。所述计算机可实现信号处理与输出系统、控制器及显示系统的功能。
与现有技术相比,本发明的有益效果为:
本发明提供的一种多相体系中颗粒浓度和粒径的在线测量方法,利用一台基于远心照相的多相测量仪所拍摄多相反应器内的多相流动图像,能够实时在线获得多相反应器内气泡、液滴或固体等颗粒浓度和粒径分布,并且测量的准确率较高,粒度测量误差<3%,浓度误差<10%。
附图说明
图1是本发明一种实施方式提供的利用浸入式在线多相测量仪测量多相体系中颗粒粒径和浓度的实验装置。
其中:1,搅拌槽,2,搅拌轴,3,挡板,4,搅拌桨,5,固体颗粒,6,气泡,7,照相探头,8,脉冲光源,9,图像采集计算机。
图2是本发明一种实施方式提供的照明闪光和CCD拍照同步的控制模式。
图3是实施例1提供的采用浸入式在线多相测量仪得到的液-固体系中颗粒的图像,其中,(a)是测量点瞬时图像;(b)是固体颗粒粒度测量结果。
图4是实施例2提供的采用浸入式在线多相测量仪得到的气-液体系中颗粒的图像,其中,(a)是测量点瞬时图像;(b)是气泡粒度测量结果。
图5是实施例3提供的采用浸入式在线多相测量仪得到的气-液-固体系中颗粒的图像,其中,(a)是测量点瞬时图像;(b)固体颗粒粒度测量结果;(c)气泡粒度测量结果。
具体实施方式
下面结合附图并通过具体实施方式来进一步说明本发明的技术方案。
如图1所示,为浸入式在线多相测量仪测量多相体系中颗粒粒径和浓度的实验装置。
所述多相测量仪包括:
不锈钢封装管;
视窗,密封连接于不锈钢封装管一端;
用于照明多相流的照明系统,包括LED灯和与LED灯相连的可调制光源,可调制光源(即脉冲光源8)包括电源、信号发生器和示波器;
用于拍照的照相系统,包括远心镜头和图像传感器,所述图像传感器为微型高速CMOS相机;
与信号发生器和图像传感器相连的控制器;
与图像传感器相连的信号处理与输出系统;
与信号处理与输出系统相连的显示系统;
其中,LED灯、远心镜头和图像传感器位于不锈钢封装管中形成照相探头7,可调制光源、控制器、信号处理与输出系统及显示系统位于不锈钢封装管外;控制器控制图像传感器的曝光周期小于信号发生器的脉冲周期。
信号处理与输出系统、控制器及显示系统集成为图像采集计算机9。
具体地,不锈钢封装管内最前端为一视窗,视窗为一圆形内侧镀增透膜蓝宝石玻璃;视窗后均布20个高亮度LED灯,20个LED灯组成环形,并均匀分布;在LED灯后端安装远心镜头,远心镜头相关参数为:倍率1,物方及像方视野均为φ8mm(φ为直径),工作距离250mm±3%,远心度<0.1°,景深2.1mm,分辨率14.3μm,光学畸变<0.12%,为清晰成像,远心镜头的前端距离视窗1外侧表面距离为远心镜头的工作距。远心镜头通过标准C口与微型高速CMOS相机连接,CMOS相机参数为:分辨率1280×1024,颜色黑白,帧速150fps,接口USB3.0。视窗1、LED灯2、远心镜头4和微型高速CMOS相机5封装在不锈钢封装管内。所述在线多相测量仪外配置一个可调制光源,可调制光源通过导线与LED灯相连。远心镜头通过USB3.0数据传输线与安装有高速图像采集卡的图像采样计算机9相连。
为了获得清晰的图像,采用如图2所示的控制方式实现照明闪光和CCD拍照的同步,所述控制方式为开启光源开关,通过光源驱动器设置脉冲光的强度和周期,通过计算机上的控制器设置图像拍摄的曝光时间、光平衡、帧频和增益,使照明信号的脉冲周期与图像传感器的曝光时间相匹配(图像传感器的曝光周期小于信号发生器的脉冲周期),以实现脉冲光和图像拍摄的同步。
所述实验装置为搅拌槽1,所述搅拌槽中放置有搅拌轴2和搅拌桨3,用于将多相体系混合均匀。
实施例1:利用浸入式在线多相测量仪对液-固体系内固体颗粒的粒径分布进行测量
实验在如图1所示的搅拌槽中进行,搅拌转速为480rpm。所用固体颗粒为粒径大小在1mm左右的白色塑料珠,整体平均固含率(体积比)为0.01,测量点位置为r=0.07m,z=0.055m。
所述测量方法包括如下步骤:
(1)将在线多相测量仪置于所述液-固体系中,得到液-固体系中固体颗粒的图像,如图3-(a)所示;
(2)确定有效颗粒:
首先,确定焦平面位置:将一个精确度至少为0.1mm的刻度尺置于与所测多相体系相同的介质中,采用在线多相测量仪对其拍照,调整刻度尺到封装管前端的距离,观察照片中刻度尺的清晰程度,刻度尺最清晰的位置为焦平面位置;
之后,将被测物置于封装管前方、焦平面前后l/2处,其中,l为远心镜头的景深,单位为mm,采用在线多相测量仪对被测物进行拍照,得到被测物图像,识别读取图像中被测物边界灰度梯度值Grad(Φl/2),Φl/2为焦平面前后两侧l/2平面处灰度值,将颗粒边界处灰度梯度Grad(Φ)≥Grad(Φl/2)的颗粒表示为有效颗粒;
(3)将一个精确度至少为0.1mm的刻度尺置于在线多相测量仪前方,采用图像传感器驱动控制软件拍摄图像,通过图像处理软件拾取刻度尺上10mm距离所对应的像素个数N10,从而确定单个像素的实际长度,并测得有效颗粒所占的像素个数,所述有效颗粒的粒径为di=10*ni/N10,其中,di为有效颗粒的粒径,单位为mm;ni为图像中所测有效颗粒对应的像素数;N10为刻度尺上10mm长度对应的像素数;
根据计算有效颗粒的浓度,其中,S为图像传感器感光部件的有效面积,单位为mm2;l为远心镜头的景深,单位为mm;di为有效颗粒的粒径,单位为mm;Vc为有效颗粒的总体积,V为测量区域的体积;n为有效颗粒的个数。
通过对近4000个颗粒进行分析,得到了测量点处的固体颗粒的粒径分布,结果如图3-(b)所示,测量结果与实际固体颗粒尺寸很接近。统计得到该点处的固体颗粒浓度为0.0112,与光纤颗粒浓度测量仪得到的0.0109相近,误差为2.75%。
实施例2:利用浸入式在线多相测量仪对气-液体系内气泡的粒径分布进行测量
实验在如图1所示的搅拌槽中进行,搅拌转速为450rpm,通过环形气体分布器进气,通气量为800L/h。测量点位置设置与实施例1中相同。
所述测量方法包括如下步骤:
(1)将在线多相测量仪置于所述气-液体系中,得到气-液体系中气泡的图像,如图4-(a)所示;
(2)确定有效颗粒:
首先,确定焦平面位置:将一个精确度至少为0.1mm的刻度尺置于与所测多相体系相同的介质中,采用在线多相测量仪对其拍照,调整刻度尺到封装管前端的距离,观察照片中刻度尺的清晰程度,刻度尺最清晰的位置为焦平面位置;
之后,将被测物置于封装管前方、焦平面前后l/2处,其中,l为远心镜头的景深,单位为mm,采用在线多相测量仪对被测物进行拍照,得到被测物图像,识别读取图像中被测物边界灰度梯度值Grad(Φl/2),Φl/2为焦平面前后两侧l/2平面处灰度值,将颗粒边界处灰度梯度Grad(Φ)≥Grad(Φl/2)的颗粒表示为有效颗粒;
(3)将一个精确度至少为0.1mm的刻度尺置于在线多相测量仪前方,采用图像传感器控制软件拍摄图像,通过图像处理软件拾取刻度尺上10mm距离所对应的像素个数N10,从而确定单个像素的实际长度,并测得有效颗粒所占的像素个数,所述有效颗粒的粒径为di=10*ni/N10,其中,di为有效颗粒粒径,单位为mm;ni为图像中所测有效颗粒对应的像素数;N10为刻度尺上10mm长度对应的像素数;
根据计算得到有效颗粒的浓度,其中,S为图像传感器感光部件的有效面积,单位为mm2;l为远心镜头的景深,单位为mm;di为有效颗粒粒径,单位为mm;Vc为有效颗粒的总体积,V为测量区域的体积;n为有效颗粒的个数。
与固体颗粒相比,气泡尺寸较大易发生变形。通过对近4000个气泡进行分析,气泡尺寸计算结果如图4-(b)所示,反应器内气泡尺寸主要集中在1~3mm范围内。计算得到该点处气含率为0.022。
实施例3:利用浸入式在线多相测量仪测量气-液-固三相搅拌槽内固体颗粒以及气泡的粒径分布
实验在如图1所示的搅拌槽中进行,使用的固体颗粒材料及体积与实施例1相同,通气条件与实施例2相同,搅拌转速480rpm。测量点位置设置与实施例1中相同
所述测量方法包括如下步骤:
(1)将在线多相测量仪置于所述气-液-固体系中,得到图像,如图5-(a)所示;
(2)确定有效颗粒:
首先,确定焦平面位置:将一个精确度至少为0.1mm的刻度尺置于与所测多相体系相同的介质中,采用在线多相测量仪对其拍照,调整刻度尺到封装管前端的距离,观察照片中刻度尺的清晰程度,刻度尺最清晰的位置为焦平面位置;
之后,将被测物置于封装管前方、焦平面前后l/2处,其中,l为远心镜头的景深,单位为mm,采用在线多相测量仪对被测物进行拍照,得到被测物图像,识别读取图像中被测物边界灰度梯度值Grad(Φl/2),Φl/2为焦平面前后两侧l/2平面处灰度值,将颗粒边界处灰度梯度Grad(Φ)≥Grad(Φl/2)的颗粒表示为有效颗粒;
(3)将一个精确度至少为0.1mm的刻度尺置于在线多相测量仪前方,采用图像传感器控制软件拍摄图像,通过图像处理软件拾取刻度尺上10mm距离所对应的像素个数N10,从而单个像素的实际长度,并测得有效颗粒所占的像素个数,所述有效颗粒的粒径为di=10*ni/N10,其中,di为有效颗粒粒径,单位为mm;ni为图像中所测有效颗粒对应的像素数;N10为刻度尺上10mm长度对应的像素数;
根据计算得到有效颗粒的浓度,其中,S为图像传感器感光部件的有效面积,单位为mm2;l为远心镜头的景深,单位为mm,di为有效颗粒粒径,单位为mm;Vc为有效颗粒的总体积,V为测量区域的体积;n为有效颗粒的个数。
通过对图像进行处理,得到了测量点处的固体颗粒和气泡尺寸分布,结果如图5-(b)和5-(c)所示,固体颗粒的粒径集中在0.9-1.1mm,气泡的粒径集中在0.5-2.8mm。计算得到该点处气含率和固含率分别为0.016和0.0115。
申请人声明,以上所述仅为本发明的具体实施方式,但本发明的保护范围并不局限于此,所属技术领域的技术人员应该明了,任何属于本技术领域的技术人员在本发明揭露的技术范围内,可轻易想到的变化或替换,均落在本发明的保护范围和公开范围之内。
Claims (28)
1.一种多相体系中气泡或颗粒的浓度和粒径的在线测量方法,其特征在于,所述方法基于一种浸入式在线多相测量仪,所述浸入式在线多相测量仪包括:
封装管;
视窗,密封连接于封装管一端;
用于照明多相流的照明系统,包括LED灯和与LED灯相连的可调制光源,可调制光源包括电源、信号发生器和示波器;
用于拍照的照相系统,包括远心镜头和图像传感器;
与信号发生器和图像传感器相连的控制器;
与图像传感器相连的信号处理与输出系统;
与信号处理与输出系统相连的显示系统;
其中,LED灯、远心镜头和图像传感器位于封装管中,控制器控制图像传感器的曝光周期小于信号发生器的脉冲周期;
所述测量方法包括如下步骤:
(1)将浸入式在线多相测量仪置于多相体系中,得到多相体系中气泡或颗粒的图像;
(2)确定有效气泡或颗粒:首先,确定远心镜头的焦平面位置;之后,将被测物置于封装管前方、焦平面前后l/2处,其中,l为远心镜头的景深,单位为mm,采用浸入式在线多相测量仪对被测物进行拍照,得到被测物图像,识别读取图像中被测物边界灰度梯度值Grad(Φl/2),Φl/2为焦平面前后两侧l/2平面处灰度值,将气泡或颗粒边界处灰度梯度Grad(Φ)≥Grad(Φl/2)的气泡或颗粒表示为有效气泡或颗粒;
(3)确定单个像素的实际长度,并测得有效气泡或颗粒所占的像素个数,则有效气泡或颗粒粒径di为有效气泡或颗粒所占的像素个数与单个像素实际长度的乘积;
根据计算得到有效气泡或颗粒的浓度,其中,S为图像传感器感光部件的有效面积,单位为mm2;l为远心镜头的景深,单位为mm,di为有效气泡或颗粒粒径,单位为mm;Vc为有效气泡或颗粒的总体积,V为测量区域的体积;n为有效气泡或颗粒的个数。
2.根据权利要求1所述的方法,其特征在于,步骤(1)所述多相体系中的颗粒为深色固体颗粒或液滴时,浸入式在线多相测量仪拍照时需加装反光板,被测颗粒置于反光板和浸入式在线多相测量仪前端面之间。
3.根据权利要求2所述的方法,其特征在于,所述反光板与浸入式在线多相测量仪前端面的距离大于远心镜头的景深。
4.根据权利要求2所述的方法,其特征在于,所述反光板与浸入式在线多相测量仪前端面的距离为1.5倍的远心镜头景深。
5.根据权利要求2所述的方法,其特征在于,所述反光板的尺寸与图像传感器的有效感光区域相同。
6.根据权利要求2所述的方法,其特征在于,所述反光板为中心凸起的方形板。
7.根据权利要求2所述的方法,其特征在于,所述反光板表面为亮白色或银白色。
8.根据权利要求1所述的方法,其特征在于,步骤(2)确定焦平面位置的方法为:将一个精确度至少为0.1mm的刻度尺置于与所测多相体系相同的介质中,采用浸入式在线多相测量仪对所述刻度尺拍照,调整刻度尺到浸入式多相测量仪前端的距离,观察照片中刻度尺的清晰程度,刻度尺最清晰的位置为焦平面位置。
9.根据权利要求1所述的方法,其特征在于,步骤(3)所述确定单个像素的实际长度具体为:将一个精确度至少为0.1mm的刻度尺置于浸入式在线多相测量仪前方,采用图像传感器驱动控制软件拍摄图像,通过图像处理软件拾取刻度尺上10mm距离所对应的像素个数N10,从而确定单个像素的实际长度。
10.根据权利要求9所述的方法,其特征在于,所述图像处理软件为Image-Pro Plus。
11.根据权利要求1所述的方法,其特征在于,步骤(3)所述有效气泡或颗粒的粒径为di=10*ni/N10,其中,di为气泡或颗粒尺寸,mm;ni为图像中所测有效气泡或颗粒所占的像素个数;N10为刻度尺上10mm长度对应的像素个数。
12.根据权利要求1所述的方法,其特征在于,所述远心镜头的工作距离为250~550mm,景深为1~3.7mm。
13.根据权利要求1所述的方法,其特征在于,所述远心镜头的放大倍率为0.5~1倍。
14.根据权利要求1所述的方法,其特征在于,所述远心镜头的外径为19~25mm。
15.根据权利要求1所述的方法,其特征在于,所述图像传感器为CCD相机或CMOS相机。
16.根据权利要求15所述的方法,其特征在于,所述CCD相机或CMOS相机的曝光时间≤1ms,分辨率为5~15μm,长和宽方向像素数至少为800×600,帧频至少为60fps。
17.根据权利要求1所述的方法,其特征在于,所述LED灯的个数至少为12个。
18.根据权利要求1所述的方法,其特征在于,所述LED灯在封装管中组成环形均匀排列。
19.根据权利要求1所述的方法,其特征在于,所述LED灯与可调制光源之间通过导线连接。
20.根据权利要求1所述的方法,其特征在于,所述封装管由直径不同的前段管和后段管组成。
21.根据权利要求20所述的方法,其特征在于,所述前段管的外径为25~30mm,长度为300~600mm。
22.根据权利要求20所述的方法,其特征在于,所述后段管的外径为50mm,长度为50mm。
23.根据权利要求1所述的方法,其特征在于,所述封装管的材质为不锈钢。
24.根据权利要求20所述的方法,其特征在于,所述视窗、LED灯和远心镜头封装于前段管中,所述视窗置于前段管远离后段管的一端,其后依次设置LED灯和远心镜头,所述图像传感器封装于后段管中。
25.根据权利要求1所述的方法,其特征在于,所述视窗为内侧镀增透膜的圆形玻璃。
26.根据权利要求1所述的方法,其特征在于,所述信号发生器和图像传感器通过高速数据线与控制器相连。
27.根据权利要求1所述的方法,其特征在于,所述显示系统为显示屏。
28.根据权利要求1所述的方法,其特征在于,所述信号处理与输出系统、控制器及显示系统集成为计算机。
Priority Applications (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610245330.6A CN105928847B (zh) | 2016-04-19 | 2016-04-19 | 一种多相体系中颗粒浓度和粒径的在线测量方法 |
US15/368,838 US10209064B2 (en) | 2016-04-19 | 2016-12-05 | Online measuring method of concentration and diameter of particles in multiphase system |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610245330.6A CN105928847B (zh) | 2016-04-19 | 2016-04-19 | 一种多相体系中颗粒浓度和粒径的在线测量方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN105928847A CN105928847A (zh) | 2016-09-07 |
CN105928847B true CN105928847B (zh) | 2018-06-26 |
Family
ID=56839322
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201610245330.6A Active CN105928847B (zh) | 2016-04-19 | 2016-04-19 | 一种多相体系中颗粒浓度和粒径的在线测量方法 |
Country Status (2)
Country | Link |
---|---|
US (1) | US10209064B2 (zh) |
CN (1) | CN105928847B (zh) |
Families Citing this family (18)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106373519A (zh) * | 2016-09-23 | 2017-02-01 | 京东方科技集团股份有限公司 | 显示装置及其控制方法 |
CN106872526B (zh) * | 2017-01-09 | 2019-05-21 | 浙江大学 | 一种非接触式气液固三相鼓泡床相含率测量装置及方法 |
CN107389509B (zh) * | 2017-07-07 | 2019-06-04 | 中国科学院过程工程研究所 | 一种提高多相体系中低亮度分散相颗粒辨识度的方法 |
US10959411B2 (en) | 2018-01-04 | 2021-03-30 | Atlantic Sapphire IP, L.L.C. | Bottom grading apparatuses for aquaculture systems |
CN108693078B (zh) * | 2018-04-09 | 2020-10-16 | 青岛科技大学 | 一种利用界面现象确定co2/油相体系一次混相压力的新方法 |
CN110333127B (zh) | 2019-06-24 | 2020-06-12 | 北京大学 | 一种气相半挥发性有机物在线测量系统、方法和应用 |
CN111986175A (zh) * | 2020-08-19 | 2020-11-24 | 北京科技大学 | 工业喷嘴所喷液滴的粒径测量方法 |
CN112485166B (zh) * | 2020-11-18 | 2023-03-31 | 青岛科技大学 | 一种在线测量结晶形状和粒度分布的装置 |
CN112577859B (zh) * | 2020-12-02 | 2023-05-30 | 苏州海狸生物医学工程有限公司 | 一种测量磁性微球基本物理参数的实验装置及测量方法 |
CN112730174A (zh) * | 2020-12-17 | 2021-04-30 | 华南理工大学 | 一种实时监测结晶过程中类球形颗粒浓度的方法 |
US11662291B1 (en) * | 2021-02-03 | 2023-05-30 | Atlantic Sapphire Ip, Llc | System and method for feed validation measurement |
CN113642609B (zh) * | 2021-07-15 | 2024-03-26 | 东华大学 | 基于图像识别技术的共混聚合物中分散相形貌的表征方法 |
CN113739847A (zh) * | 2021-08-30 | 2021-12-03 | 中国大唐集团科学技术研究院有限公司西北电力试验研究院 | 一种测量锅炉煤粉管道煤粉浓度及速度偏差的系统及方法 |
CN114046748A (zh) * | 2021-09-27 | 2022-02-15 | 广州超音速自动化科技股份有限公司 | 一种锂电池隔膜的颗粒检测方法、装置及系统 |
CN115096768A (zh) * | 2022-08-24 | 2022-09-23 | 浙江省水利水电勘测设计院有限责任公司 | 可同时测量颗粒粒径与体积浓度的背光成像系统和方法 |
CN115624947B (zh) * | 2022-12-21 | 2023-03-14 | 广州市凯闻食品发展有限公司 | 一种高粘度固液混合分散的反应方法 |
CN117129413B (zh) * | 2023-04-07 | 2024-04-05 | 南京工业大学 | 一种颗粒黏附力和表面能的测量装置及使用方法 |
CN116359086B (zh) * | 2023-06-01 | 2023-09-08 | 南昌科晨电力试验研究有限公司 | 一种基于测量装置的煤粉粒径与质量浓度测量方法 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2014138939A1 (en) * | 2013-03-13 | 2014-09-18 | Queen's University At Kingston | Methods and systems for characterizing laser machining properties by measuring keyhole dynamics using interferometry |
DE102013226932A1 (de) * | 2013-12-20 | 2015-06-25 | Carl Zeiss Ag | Verfahren zur Ermittlung der Phasenverteilung |
CN204831972U (zh) * | 2015-06-29 | 2015-12-02 | 广州机械科学研究院有限公司 | 一种在线镶嵌式成像励磁传感器 |
CN105136803A (zh) * | 2015-04-27 | 2015-12-09 | 内蒙古科技大学 | 用于气力输送中检测固相物质浓度/颗粒度大小的装置 |
CN103868831B (zh) * | 2014-02-26 | 2016-01-20 | 天津大学 | 云粒子谱分布测量方法及测量系统 |
Family Cites Families (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6879708B2 (en) * | 2001-05-24 | 2005-04-12 | Case Western Reserve University | Planar particle/droplet size measurement technique using digital particle image velocimetry image data |
US20060116531A1 (en) * | 2004-11-29 | 2006-06-01 | Wonders Alan G | Modeling of liquid-phase oxidation |
US9234852B2 (en) * | 2005-07-29 | 2016-01-12 | Mitutoyo Corporation | Systems and methods for controlling strobe illumination |
FI122997B (fi) * | 2011-09-06 | 2012-09-28 | Janesko Oy | Menetelmä ja sovitelma optisesti epähomogeenisen aineen virtausnopeuden mittaamiseksi |
JP6276158B2 (ja) * | 2014-09-30 | 2018-02-07 | 富士フイルム株式会社 | 細胞撮像装置および方法 |
-
2016
- 2016-04-19 CN CN201610245330.6A patent/CN105928847B/zh active Active
- 2016-12-05 US US15/368,838 patent/US10209064B2/en active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2014138939A1 (en) * | 2013-03-13 | 2014-09-18 | Queen's University At Kingston | Methods and systems for characterizing laser machining properties by measuring keyhole dynamics using interferometry |
DE102013226932A1 (de) * | 2013-12-20 | 2015-06-25 | Carl Zeiss Ag | Verfahren zur Ermittlung der Phasenverteilung |
CN103868831B (zh) * | 2014-02-26 | 2016-01-20 | 天津大学 | 云粒子谱分布测量方法及测量系统 |
CN105136803A (zh) * | 2015-04-27 | 2015-12-09 | 内蒙古科技大学 | 用于气力输送中检测固相物质浓度/颗粒度大小的装置 |
CN204831972U (zh) * | 2015-06-29 | 2015-12-02 | 广州机械科学研究院有限公司 | 一种在线镶嵌式成像励磁传感器 |
Also Published As
Publication number | Publication date |
---|---|
US20170299382A1 (en) | 2017-10-19 |
CN105928847A (zh) | 2016-09-07 |
US10209064B2 (en) | 2019-02-19 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN105928847B (zh) | 一种多相体系中颗粒浓度和粒径的在线测量方法 | |
US10247750B2 (en) | Online measuring method of particle velocity in multiphase system | |
CN104303026B (zh) | 双光束非接触式位移传感器 | |
CN105928841B (zh) | 一种浸入式在线多相测量仪及测量方法 | |
CN106404623B (zh) | 悬浮泥沙浓度监测系统及监测方法 | |
US9612144B2 (en) | Method and arrangement for measuring flow rate of optically non-homogenous material | |
CN107063943A (zh) | 液体中颗粒物自动计数、图像采集及形态测量装置 | |
CN109253992A (zh) | 一种等离子体荧光寿命测量装置及方法 | |
JPS61159135A (ja) | 粒子解析装置 | |
CN109187373A (zh) | 一种基于机器视觉的番茄品质检测系统及检测方法 | |
WO2020119600A1 (zh) | 液体中颗粒物的图像采集设备及检测装置 | |
CN208071854U (zh) | 用于监控通过电纺丝或静电纺丝来生产纳米纤维的过程的装置 | |
EP0507746A2 (en) | Method and device for measuring the particle size distribution of a flowing solid particulate substance | |
US10481062B2 (en) | Method for improving identification degree of low-luminosity dispersed-phase particles in multiphase system | |
CN109297554B (zh) | 用于测量t型管内液相流量的方法 | |
CN207395666U (zh) | 基于线阵ccd的鱼苗长度测量仪结构 | |
CN111141746B (zh) | 一种笔芯尾油长度自动检测方法与系统 | |
CN110095316A (zh) | 一种粒子采集装置、检测系统及检测方法 | |
CN203101029U (zh) | 一种用于检测中继镜的装置 | |
WO2010108449A1 (en) | Improved paint mixing device | |
CN206223209U (zh) | 光学机构和光学测量仪 | |
CN210572335U (zh) | 一种高灵敏的特定蛋白分析仪 | |
JPS6314768B2 (zh) | ||
CN208118692U (zh) | 一种面膜排错系统 | |
CN209727759U (zh) | 一种光学玻璃折射率的快速测量仪器 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |