CN108871211B - 一种圆管内气液两相环状流液膜厚度的测量及修正方法 - Google Patents
一种圆管内气液两相环状流液膜厚度的测量及修正方法 Download PDFInfo
- Publication number
- CN108871211B CN108871211B CN201810327886.9A CN201810327886A CN108871211B CN 108871211 B CN108871211 B CN 108871211B CN 201810327886 A CN201810327886 A CN 201810327886A CN 108871211 B CN108871211 B CN 108871211B
- Authority
- CN
- China
- Prior art keywords
- liquid film
- circular tube
- thickness
- liquid
- wall
- 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
Images
Classifications
-
- 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/02—Measuring arrangements characterised by the use of optical techniques for measuring length, width or thickness
- G01B11/06—Measuring arrangements characterised by the use of optical techniques for measuring length, width or thickness for measuring thickness ; e.g. of sheet material
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Length Measuring Devices By Optical Means (AREA)
Abstract
本发明公开了一种圆管内气液两相环状流液膜厚度的测量及修正方法,包括以下步骤:1)通过摄影装置1结合光源3对圆管内气液两相环状流的相界面进行拍摄,得包含液膜和气芯的灰度图;2)对灰度图进行数字图像处理,得包含液膜和气芯的二值化图像;3)提取二值化图像中气液相界面信息,得液膜厚度的成像值;4)建立预测圆管壁厚及液膜厚度成像值与真实值之间定量关系的光学成像理论模型;5)利用步骤4)构建的液膜光学成像理论模型对待测量液膜厚度的成像值h成像进行修正,得待测量液膜厚度的真实值h真实,该方法能够实现圆管内气液两相环状流液膜厚度的可视化检测,具有较高的检测精度、安全性及可靠性,并且检测成本低。
Description
技术领域
本发明属于气液两相环状流参数测量技术领域,具体涉及一种圆管内气液两相环状流液膜厚度的测量及修正方法。
背景技术
气液两相环状流是一种重要的、并在工程领域应用广泛的两相流动,管内环状流主要由管壁液膜和中心雾状高速气芯组成。环状流中管壁液膜具有很强的热质交换能力,液膜厚度的准确测量,对深入研究环状流的流动及传递特征具有重要作用。根据测量原理不同,液膜厚度的测量方法主要包括声学法、电导法、光学法和射线法等。声学法基于超声波在气液界面的衰减与反射传播的时间来确定液膜厚度,但超声波波长的不确定性限制了声学法在液膜极薄情况下的应用。电导法简单可靠、成本较低,但它只适用导电液膜的测量,电导探针侵入流体对流场也会产生一定的干扰,导致测量偏差。射线法是一种成熟且应用较好的测量方法,但实际应用中涉及射线防护、放射源存储及设备维护等安全问题。光学法是目前应用较多的方法,但光学测量设备普遍价格昂贵,操作过程复杂,且对被测介质和应用环境的清洁度有严格要求。
发明内容
本发明的目的在于克服上述现有技术的缺点,提供了一种圆管内气液两相环状流液膜厚度的测量及修正方法,该方法能够精确测量气液两相环状流液膜的厚度,并且测量成本低,安全可靠。
为达到上述目的,本发明所述的圆管内气液两相环状流液膜厚度的测量及修正方法包括以下步骤:
1)在圆管的测量段外布置矩形透明水箱,在矩形透明水箱的正面及背面分别设置摄影装置及光源,开启光源,通过摄影装置拍摄圆管内气液两相环状流相界面的图像,得包含管壁、液膜及气芯的灰度图;
2)对步骤1)得到的灰度图进行数字图像处理,得包含液膜及气芯的二值化图像;
3)提取步骤2)得到的二值化图像中的气液相界面信息,得液膜厚度的成像值h成像;
4)建立预测圆管壁厚及液膜厚度成像值h成像与真实值之间定量关系的光学成像理论模型,得液膜厚度缩放因子K液膜;
5)利用步骤4)得到的液膜厚度缩放因子K液膜对待测量液膜厚度的成像值h成像进行修正,得待测量液膜厚度的真实值h真实。
步骤2)的具体操作过程为:
a)从灰度图上分割出包含气液相界面的液膜有效区域;
b)利用对比度调整函数提高灰度图的对比度;
c)利用填充函数对灰度图的空洞区域进行填充,消除液膜中夹带的气泡;
d)将灰度图转化为二值图像;
e)通过填充函数及求反函数消除气芯中夹带的液滴,得包含液膜及气芯的二值化图像。
设圆管内径及外径分别为r及R,圆管管壁厚度的真实值为R-r,光源入射光线进入圆管后在圆管外管壁处的入射角为θ2,光线穿过圆管管壁在矩形透明水箱水侧的折射角为θ1,其中,折射光线经矩形透明水箱的壁面射出后平行进入摄影装置中,圆管管壁厚度的成像值为R–r–L,其中,变量L为射出圆管的光线与圆管外壁之间的垂直距离;
折射角θ1满足的几何关系为:
根据斯涅尔定律,入射角θ2与折射角θ1满足:
其中,n水为水的折射率,n管壁为圆管的折射率;
联立式(1)、式(2)及式(3),得:
由式(4)可知,当已知圆管内径及外径、光线在水中的折射率n水及圆管管壁中的折射率n管壁,即可求得变量L的大小。
来自光源的入射光线AB从液膜与圆管内壁边界折射后以光线BC进入圆管管壁,光线BC随后在圆管外壁与矩形透明水箱水侧发生折射后平行射入摄影装置中,液膜与圆管内管壁处的入射角及折射角分别为θ4及θ3,光线在圆管外管壁与矩形透明水箱水侧处的入射角及折射角分别为θ2及θ1。管壁液膜厚度的真实值为h,设y1为射出圆管的光线与圆管轴线之间的垂直距离,则折射角θ1满足的几何关系为:
圆管外壁处的入射角θ2与矩形透明水箱2水侧的折射角θ1满足:
光线BC的轨迹方程为:
y=y1-tan(θ1-θ2)(Rcosθ1-x) (7)
圆OB的轨迹方程为:
x2+y2=r2 (8)
由式(7)及式(8),得B点的坐标位置(xB,yB);
则线段OB的斜率为:
光线在液膜处的入射角θ4与圆管内壁处的折射角θ3满足:
经过液膜的光线AB的方程为:
y=yB+xBtan[θ4-θ3-(θ1-θ2)] (11)
光线OA的方程为:
y=r-h (12)
由式(11)及式(12)得A点的坐标(xA,yA),则A点处气液相界面与圆心的距离rA为:
液膜的真实厚度h为:
h=r-rA (14)
液膜厚度的成像值h成像为:
h成像=r+L-y1 (15)
其中,变量L通过式(4)得到;
则液膜厚度的缩放因子K液膜为:
液膜厚度的真实值h真实为:
本发明具有以下有益效果:
本发明所述的圆管内气液两相环状流液膜厚度的测量及修正方法在具体操作时,以摄影装置作为监测仪器,通过摄影装置结合光源对圆管内气液两相环状流的相界面进行拍摄,得包含管壁、液膜及气芯的灰度图;然后对灰度图进行数字图像处理得二值化图像,提取二值化图像中气液相界面信息,得液膜厚度的成像值;最后构建光学成像理论模型,通过该光学成像理论模型对待测量液膜厚度的成像值进行修正,得待测量液膜厚度的真实值,操作简单,易于实现,并且整个测量过程只需通过非侵入方式对圆管内气液两相环状流相界面进行拍摄,整个检测过程中对圆管内液体的流动无干扰,检测精度较高,并且安全、可靠,成本低。
附图说明
图1a为本发明中摄影装置1、矩形透明水箱2及光源3的位置关系图;
图1b为本发明中圆管内气液两相环状流相界面的示意图;
图2a为本发明中包含液膜和气芯有效区域的灰度图像;
图2b为本发明中增强对比度后的图像;
图2c为本发明中填洞消去液膜中夹带气泡后的图像
图2d为本发明中二值化图像
图2e为本发明中填洞消去气芯中夹带液滴后的图像;
图3a为本发明中圆管管壁成像的光路图;
图3b为本发明中管壁液膜成像的光路图;
图4a为本发明中液膜厚度成像值h成像的分布图;
图4b为本发明中液膜厚度真实值h真实的分布图。
其中,1为摄影装置、2为矩形透明水箱、3为光源。
具体实施方式
下面结合附图对本发明做进一步详细描述:
参考图1a至图4b,本发明所述的圆管内气液两相环状流液膜厚度的测量及修正方法包括以下步骤:
1)在圆管的测量段外布置矩形透明水箱2,在矩形透明水箱2的正面及背面分别设置摄影装置1及光源3,开启光源3,通过摄影装置1拍摄圆管内气液两相环状流相界面的图像,得包含管壁、液膜及气芯的灰度图;
2)对步骤1)得到的灰度图进行数字图像处理,得包含液膜及气芯的二值化图像;
3)提取步骤2)得到的二值化图像中的气液相界面信息,得液膜厚度的成像值h成像;
4)建立预测圆管壁厚及液膜厚度成像值与真实值之间定量关系的光学成像理论模型,得液膜厚度缩放因子K液膜;
5)利用步骤4)得到的液膜厚度缩放因子K液膜对待测量液膜厚度的成像值h成像进行修正,得待测量液膜厚度的真实值h真实。
步骤2)的具体操作过程为:
a)从灰度图上分割出包含气液相界面的液膜有效区域;
b)利用对比度调整函数提高灰度图的对比度;
c)利用填充函数对灰度图的空洞区域进行填充,消除液膜中夹带的气泡;
d)将灰度图转化为二值图像;
e)通过填充函数及求反函数消除气芯中夹带的液滴,得包含液膜及气芯的二值化图像。
步骤4)的具体操作为:
设圆管内径及外径分别为r及R,圆管管壁厚度的真实值为R-r,光源3入射光线进入圆管后在圆管外管壁处的入射角为θ2,光线穿过圆管管壁在矩形透明水箱2水侧的折射角为θ1,其中,折射光线经矩形透明水箱2的壁面射出后平行进入摄影装置1中,圆管管壁厚度的成像值为R–r–L,其中,变量L为射出圆管的光线与圆管外壁之间的垂直距离;
折射角θ1满足的几何关系为:
根据斯涅尔定律,入射角θ2与折射角θ1满足:
其中,n水为水的折射率,n管壁为圆管的折射率;
联立式(1)、式(2)及式(3),得:
由式(4)可知,当已知圆管内径及外径、光线在水中的折射率n水及圆管管壁中的折射率n管壁,即可求得变量L的大小。
来自光源3的入射光线AB从液膜与圆管内壁边界折射后以光线BC进入圆管管壁,光线BC随后在圆管外壁与矩形透明水箱2水侧发生折射后平行射入摄影装置1中,液膜与圆管内管壁处的入射角及折射角分别为θ4及θ3,光线在圆管外管壁与矩形透明水箱2水侧处的入射角及折射角分别为θ2及θ1,设y1为射出圆管的光线与圆管轴线之间的垂直距离,则折射角θ1满足几何关系:
圆管外壁处的入射角θ2与矩形透明水箱2水侧的折射角θ1满足:
光线BC的轨迹方程为:
y=y1-tan(θ1-θ2)(Rcosθ1-x) (7)
圆OB的轨迹方程为:
x2+y2=r2 (8)
由式(7)及式(8),得B点的坐标位置(xB,yB);
则线段OB的斜率为:
光线在液膜处的入射角θ4与圆管内壁处的折射角θ3满足:
经过液膜的光线AB的方程为:
y=yB+xBtan[θ4-θ3-(θ1-θ2)] (11)
光线OA的方程为:
y=r-h (12)
由式(11)及式(12)得A点的坐标(xA,yA),则A点处气液相界面与圆心的距离rA为:
液膜的真实厚度h为:
h=r-rA (14)
液膜厚度的成像值h成像为:
h成像=r+L-y1 (15)
其中,变量L通过式(4)得到;
则液膜厚度的缩放因子K液膜为:
步骤5)中液膜厚度的真实值h真实为:
本发明在具体实施时,只需采用摄影装置1对圆管内气液两相环状流相界面进行拍摄,通过数字图像处理获得液膜厚度成像值,并利用构建的液膜厚度光学成像理论模型对液膜厚度成像值进行修正,即可实现对液膜真实厚度的非侵入、高精度测量,对工业实践中的环状流检测具有重要的指导意义。
Claims (2)
1.一种圆管内气液两相环状流液膜厚度的测量及修正方法,其特征在于,包括以下步骤:
1)在圆管的测量段外布置矩形透明水箱(2),在矩形透明水箱(2)的正面及背面分别设置摄影装置(1)及光源(3),开启光源(3),通过摄影装置(1)拍摄圆管内气液两相环状流相界面的图像,得包含管壁、液膜及气芯的灰度图;
2)对步骤1)得到的灰度图进行数字图像处理,得包含液膜及气芯的二值化图像;
3)提取步骤2)得到的二值化图像中的气液相界面信息,得液膜厚度的成像值h成像;
4)建立预测圆管壁厚及液膜厚度成像值h成像与真实值之间定量关系的光学成像理论模型,得液膜厚度缩放因子K液膜;
5)利用步骤4)得到的液膜厚度缩放因子K液膜对待测量液膜厚度的成像值进行修正,得待测量液膜厚度的真实值h真实;
步骤2)的具体操作过程为:
a)从灰度图上分割出包含气液相界面的液膜有效区域;
b)利用对比度调整函数提高灰度图的对比度;
c)利用填充函数对灰度图的空洞区域进行填充,消除液膜中夹带的气泡;
d)将灰度图转化为二值图像;
e)通过填充函数及求反函数消除气芯中夹带的液滴,得包含液膜及气芯的二值化图像;
设圆管内径及外径分别为r及R,圆管管壁厚度的真实值为R-r,光源(3)入射光线进入圆管后在圆管外管壁处的入射角为θ2,光线穿过圆管管壁在矩形透明水箱(2)水侧的折射角为θ1,其中,折射光线经矩形透明水箱(2)的壁面射出后平行进入摄影装置(1)中,圆管管壁厚度的成像值为R–r-L,其中,变量L为射出圆管的光线与圆管外壁之间的垂直距离;
折射角θ1满足的几何关系为:
根据斯涅尔定律,入射角θ2与折射角θ1满足:
其中,n水为水的折射率,n管壁为圆管的折射率;
联立式(1)、式(2)及式(3),得:
由式(4)可知,当已知圆管内径及外径、光线在水中的折射率n水及圆管管壁中的折射率n管壁,即可求得变量L的大小;
来自光源(3)的入射光线AB从液膜与圆管内壁边界折射后以光线BC进入圆管管壁,光线BC随后在圆管外壁与矩形透明水箱(2)水侧发生折射后平行进入摄影装置(1)中,液膜与圆管内管壁处的入射角及折射角分别为θ4及θ3,光线在圆管外管壁与矩形透明水箱(2)水侧处的入射角及折射角分别为θ2及θ1,设y1为射出圆管的光线与圆管轴线之间的垂直距离,则折射角θ1满足几何关系:
圆管外壁处的入射角θ2与矩形透明水箱(2)水侧的折射角θ1满足:
光线BC的轨迹方程为:
y=y1-tan(θ1-θ2)(Rcosθ1-x) (7)
圆OB的轨迹方程为:
x2+y2=r2 (8)
由式(7)及式(8),得B点的坐标位置(xB,yB);
则线段OB的斜率为:
光线在液膜处的入射角θ4与圆管内壁处的折射角θ3满足:
经过液膜的光线AB的方程为:
y=yB+xBtan[θ4-θ3-(θ1-θ2)] (11)
光线OA的方程为:
y=r-h (12)
由式(11)及式(12)得A点的坐标(xA,yA),则A点处气液相界面与圆心的距离rA为:
液膜的真实厚度h为:
h=r-rA (14)
液膜厚度的成像值h成像为:
h成像=r+L-y1 (15)
其中,变量L通过式(4)得到;
则液膜厚度的缩放因子K液膜为:
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810327886.9A CN108871211B (zh) | 2018-04-12 | 2018-04-12 | 一种圆管内气液两相环状流液膜厚度的测量及修正方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810327886.9A CN108871211B (zh) | 2018-04-12 | 2018-04-12 | 一种圆管内气液两相环状流液膜厚度的测量及修正方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN108871211A CN108871211A (zh) | 2018-11-23 |
CN108871211B true CN108871211B (zh) | 2020-05-22 |
Family
ID=64326325
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201810327886.9A Active CN108871211B (zh) | 2018-04-12 | 2018-04-12 | 一种圆管内气液两相环状流液膜厚度的测量及修正方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN108871211B (zh) |
Families Citing this family (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110907456B (zh) * | 2019-12-09 | 2024-07-23 | 湖南科技大学 | 一种管道环状流气液界面波波动特征参数测量方法 |
CN114518647B (zh) * | 2020-11-20 | 2024-06-25 | 中国石油天然气股份有限公司 | 微米级孔道光学信息确定方法、装置及光源装置 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2007139606A (ja) * | 2005-11-18 | 2007-06-07 | Fuji Xerox Co Ltd | 膜厚測定装置及び膜厚測定方法 |
CN103424406A (zh) * | 2013-09-03 | 2013-12-04 | 上海理工大学 | 管道内气液两相流的图像法测量装置及方法 |
CN105701842A (zh) * | 2016-01-08 | 2016-06-22 | 天津大学 | 基于链码轮廓特征的液膜波动速度测量方法 |
CN106770115A (zh) * | 2016-12-26 | 2017-05-31 | 天津大学 | 管道内周向液膜特征光学畸变校正与测量方法 |
CN107270980A (zh) * | 2017-07-21 | 2017-10-20 | 西安交通大学 | 一种垂直管内气液两相环状流液膜厚度及流量的测量方法 |
-
2018
- 2018-04-12 CN CN201810327886.9A patent/CN108871211B/zh active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2007139606A (ja) * | 2005-11-18 | 2007-06-07 | Fuji Xerox Co Ltd | 膜厚測定装置及び膜厚測定方法 |
CN103424406A (zh) * | 2013-09-03 | 2013-12-04 | 上海理工大学 | 管道内气液两相流的图像法测量装置及方法 |
CN105701842A (zh) * | 2016-01-08 | 2016-06-22 | 天津大学 | 基于链码轮廓特征的液膜波动速度测量方法 |
CN106770115A (zh) * | 2016-12-26 | 2017-05-31 | 天津大学 | 管道内周向液膜特征光学畸变校正与测量方法 |
CN107270980A (zh) * | 2017-07-21 | 2017-10-20 | 西安交通大学 | 一种垂直管内气液两相环状流液膜厚度及流量的测量方法 |
Non-Patent Citations (1)
Title |
---|
数字图像技术在液膜厚度测量中的应用;盛伟等;《东北电力技术》;20071231(第2期);第45页第2.1节及图2 * |
Also Published As
Publication number | Publication date |
---|---|
CN108871211A (zh) | 2018-11-23 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN111006610B (zh) | 一种基于结构光三维测量的水下三维测量数据校正方法 | |
CN107449756B (zh) | 一种icf靶丸内冰层折射率与厚度的测量方法及装置 | |
CN108871211B (zh) | 一种圆管内气液两相环状流液膜厚度的测量及修正方法 | |
CN103017684A (zh) | 同轴光检测柱孔圆度和直线度的装置及方法 | |
CN104316443A (zh) | 一种基于ccd后向散射的pm 2.5浓度监测方法 | |
CN108303570B (zh) | 一种多普勒海流计声波散射区域的标定装置及方法 | |
CN108844959B (zh) | 一种圆管内气液两相环状流截面相含率的测量及修正方法 | |
CN204831220U (zh) | 氟化钙平晶两面平行度高精度测试装置 | |
CN101922932B (zh) | 一种角锥棱镜坐标测量误差的补偿装置 | |
CN107121712A (zh) | 差分光柱像运动与闪烁激光雷达实时测量湍流廓线的方法 | |
CN107727744B (zh) | 用于岩石力学三轴试验的声发射源定位方法及系统 | |
CN117006958A (zh) | 一种高深径比小孔内表面几何特征精密测量方法 | |
CN103868455B (zh) | 一种视觉重建水槽内目标点空间坐标的方法 | |
CN103453850B (zh) | 基于数字散斑相关技术的透明液面微形貌测量方法和系统 | |
CN117091660A (zh) | 一种基于机器视觉的液体微流量测量方法 | |
CN105806780A (zh) | 基于多视觉的小通道气液两相流相含率测量系统和方法 | |
CN113916139B (zh) | 一种用于测量泡沫尺寸分布的装置及方法 | |
CN115096768A (zh) | 可同时测量颗粒粒径与体积浓度的背光成像系统和方法 | |
CN113252603B (zh) | 一种多层透明球床的最优折射率测量方法 | |
CN105931250B (zh) | 一种非球形气泡的等效粒径转换方法 | |
CN112964209B (zh) | 一种基于接触式测量的离轴检测法 | |
CN108801377A (zh) | 一种针对特殊流体流速和流量测定的光学装置 | |
CN110907456B (zh) | 一种管道环状流气液界面波波动特征参数测量方法 | |
CN105180820B (zh) | 一种透明毛细管内表面面型测试系统及测试方法 | |
CN112285059A (zh) | 基于ccd方法测量液体折射率装置 |
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 |