CN108871211A - 一种圆管内气液两相环状流液膜厚度的测量及修正方法 - Google Patents
一种圆管内气液两相环状流液膜厚度的测量及修正方法 Download PDFInfo
- Publication number
- CN108871211A CN108871211A CN201810327886.9A CN201810327886A CN108871211A CN 108871211 A CN108871211 A CN 108871211A CN 201810327886 A CN201810327886 A CN 201810327886A CN 108871211 A CN108871211 A CN 108871211A
- Authority
- CN
- China
- Prior art keywords
- liquid film
- thickness
- round tube
- liquid
- gas
- 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.)
- Granted
Links
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 (5)
1.一种圆管内气液两相环状流液膜厚度的测量及修正方法,其特征在于,包括以下步骤:
1)在圆管的测量段外布置矩形透明水箱(2),在矩形透明水箱(2)的正面及背面分别设置摄影装置(1)及光源(3),开启光源(3),通过摄影装置(1)拍摄圆管内气液两相环状流相界面的图像,得包含管壁、液膜及气芯的灰度图;
2)对步骤1)得到的灰度图进行数字图像处理,得包含液膜及气芯的二值化图像;
3)提取步骤2)得到的二值化图像中的气液相界面信息,得液膜厚度的成像值h成像;
4)建立预测圆管壁厚及液膜厚度成像值h成像与真实值之间定量关系的光学成像理论模型,得液膜厚度缩放因子K液膜;
5)利用步骤4)得到的液膜厚度缩放因子K液膜对待测量液膜厚度的成像值进行修正,得待测量液膜厚度的真实值h真实。
2.根据权利要求1所述的圆管内气液两相环状流液膜厚度的测量与修正方法,其特征在于,步骤2)的具体操作过程为:
a)从灰度图上分割出包含气液相界面的液膜有效区域;
b)利用对比度调整函数提高灰度图的对比度;
c)利用填充函数对灰度图的空洞区域进行填充,消除液膜中夹带的气泡;
d)将灰度图转化为二值图像;
e)通过填充函数及求反函数消除气芯中夹带的液滴,得包含液膜及气芯的二值化图像。
3.根据权利要求2所述的圆管内气液两相环状流液膜厚度的测量及修正方法,其特征在于,
设圆管内径及外径分别为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的大小。
4.根据权利要求3所述的圆管内气液两相环状流液膜厚度的测量及修正方法,其特征在于,
来自光源(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.根据权利要求4所述的圆管内气液两相环状流液膜厚度的测量及修正方法,其特征在于,液膜厚度的真实值h真实为:
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 true CN108871211A (zh) | 2018-11-23 |
CN108871211B 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) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110907456A (zh) * | 2019-12-09 | 2020-03-24 | 湖南科技大学 | 一种管道环状流气液界面波波动特征参数测量方法 |
CN114518647A (zh) * | 2020-11-20 | 2022-05-20 | 中国石油天然气股份有限公司 | 微米级孔道光学信息确定方法、装置及光源装置 |
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 |
---|
盛伟等: "数字图像技术在液膜厚度测量中的应用", 《东北电力技术》 * |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110907456A (zh) * | 2019-12-09 | 2020-03-24 | 湖南科技大学 | 一种管道环状流气液界面波波动特征参数测量方法 |
CN114518647A (zh) * | 2020-11-20 | 2022-05-20 | 中国石油天然气股份有限公司 | 微米级孔道光学信息确定方法、装置及光源装置 |
Also Published As
Publication number | Publication date |
---|---|
CN108871211B (zh) | 2020-05-22 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Walsh et al. | Photographic measurements of bubble populations from breaking wind waves at sea | |
CN102305598B (zh) | 半刚性自回弹反射器型面精度的水下摄影测量方法 | |
CN111006610A (zh) | 一种基于结构光三维测量的水下三维测量数据校正方法 | |
CN108871211A (zh) | 一种圆管内气液两相环状流液膜厚度的测量及修正方法 | |
US20130057675A1 (en) | Method and arrangement for measuring flow rate of optically non-homogeneous material | |
CN108844959B (zh) | 一种圆管内气液两相环状流截面相含率的测量及修正方法 | |
CN111458108A (zh) | 一种发射接收光轴平行度测量装置及测量方法 | |
CN109001217A (zh) | 一种手机外壳缺陷的检测方法及检测装置 | |
CN114235809A (zh) | 一种光学非接触式的测量管道内部大颗粒的方法 | |
CN103868455B (zh) | 一种视觉重建水槽内目标点空间坐标的方法 | |
JP5035606B2 (ja) | 河川流量算出装置、河川流量算出方法およびコンピュータプログラム | |
Liu et al. | Development of a three-axis accelerometer and large-scale particle image velocimetry (LSPIV) to enhance surface velocity measurements in rivers | |
Paolillo et al. | Perspective camera model with refraction correction for optical velocimetry measurements in complex geometries | |
CN110320143A (zh) | 一种煤岩裂隙腔内三维粒子图像测速试验系统 | |
CN107505240B (zh) | 低气压环境下水体中气核自然分布的观测装置及方法 | |
CN105806780A (zh) | 基于多视觉的小通道气液两相流相含率测量系统和方法 | |
Monica et al. | Application of photogrammetric 3D-PTV technique to track particles in porous media | |
Machicoane et al. | Recent developments in particle tracking diagnostics for turbulence research | |
JPS58109875A (ja) | レ−ザ−光を利用した気泡検出方法 | |
CN108801377A (zh) | 一种针对特殊流体流速和流量测定的光学装置 | |
Padmanaban | Film thickness measurements in falling annular films | |
CN113358191A (zh) | 基于条纹投影结构光的全域洪水水位实时监测方法 | |
CN110907456A (zh) | 一种管道环状流气液界面波波动特征参数测量方法 | |
CN104535462A (zh) | 一种原位实时测量悬移质浓度和级配的装置及方法 | |
JPH03289504A (ja) | 気泡計測装置 |
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 |