CN106645396A - 一种气液两相流含气率测量方法 - Google Patents
一种气液两相流含气率测量方法 Download PDFInfo
- Publication number
- CN106645396A CN106645396A CN201610858540.2A CN201610858540A CN106645396A CN 106645396 A CN106645396 A CN 106645396A CN 201610858540 A CN201610858540 A CN 201610858540A CN 106645396 A CN106645396 A CN 106645396A
- Authority
- CN
- China
- Prior art keywords
- exp
- projection
- revised
- sample
- void fraction
- 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.)
- Pending
Links
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N29/00—Investigating or analysing materials by the use of ultrasonic, sonic or infrasonic waves; Visualisation of the interior of objects by transmitting ultrasonic or sonic waves through the object
- G01N29/02—Analysing fluids
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N29/00—Investigating or analysing materials by the use of ultrasonic, sonic or infrasonic waves; Visualisation of the interior of objects by transmitting ultrasonic or sonic waves through the object
- G01N29/44—Processing the detected response signal, e.g. electronic circuits specially adapted therefor
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N29/00—Investigating or analysing materials by the use of ultrasonic, sonic or infrasonic waves; Visualisation of the interior of objects by transmitting ultrasonic or sonic waves through the object
- G01N29/44—Processing the detected response signal, e.g. electronic circuits specially adapted therefor
- G01N29/4463—Signal correction, e.g. distance amplitude correction [DAC], distance gain size [DGS], noise filtering
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N2291/00—Indexing codes associated with group G01N29/00
- G01N2291/02—Indexing codes associated with the analysed material
- G01N2291/024—Mixtures
- G01N2291/02433—Gases in liquids, e.g. bubbles, foams
Landscapes
- Physics & Mathematics (AREA)
- Analytical Chemistry (AREA)
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Chemical & Material Sciences (AREA)
- Biochemistry (AREA)
- General Health & Medical Sciences (AREA)
- General Physics & Mathematics (AREA)
- Immunology (AREA)
- Pathology (AREA)
- Signal Processing (AREA)
- Engineering & Computer Science (AREA)
- Acoustics & Sound (AREA)
- Investigating Or Analyzing Materials By The Use Of Ultrasonic Waves (AREA)
Abstract
本发明涉及一种气液两相流含气率测量方法。本发明通过扇形束扫描方式,去除了伪影,提高了重建精度和分辨率,从而提高含气率的计算精度。在相同数量投影值的进行图像重建的结果,扇形束扫描方式明显优于平行束扫描方式,而且扇形束较平行束扫描方式效率高,操作较为简便和快捷。
Description
技术领域
本发明属于检测技术领域,具体涉及一种气液两相流含气率测量方法。
背景技术
超声层析成像技术是以Radon变换作为理论基础的一种无损检测技术,其定义为“利用声源对待测物体进行扫描,再利用相应的探测器获得与物体物理参数相关的特征值,然后利用反演算法,将特征值作为投影数据,进行图像反演,重建物体的形状、尺寸、坐标位置、数目和温度等参数”。
超声层析成像技术可分为两步:1)通过声源对待测区域进行扫描,使用相应位置的超声探测器获取接收信号,并在接收信号中获取有效的特征值;即Radon变换的正变换的过程;2)将提取的特征值作为投影数据,进行反演成像,即Radon逆变换。
根据平行束的扫描机制,以线阵形式,利用时延控制线阵上的每个阵元发射超声波扫描待测区域,在相应一端也以线阵形式的探测器接收传播过来的超声波。这种扫描方式表面上较为简单,其实存在很多的不足之处,首先以时延控制线阵上每个阵元的发射声波,其实验必须大于声波在扫描区域的传播时间,且当扫描区域更换时,需调节时延;其次,在某一方向扫描完毕后,需利用电机同步转动两个线阵至另一方向进行扫描,操作上较为繁琐。综上所述,以平行束扫描机制所涉及的线阵获取投影值,其效率较低,操作复杂;时延的设置不确定,鲁棒性差。
发明内容
本发明针对现有技术的不足,提供了一种气液两相流含气率测量方法
为解决以上问题,本发明采用了以下技术手段:
步骤1)选择超声波作为激励信号,使用超声环形阵获取走时投影数据,使用时延去控制环形阵中每个阵元的发射时间顺序,分别获取扫描区域为无气泡水域和向含气泡水域的走时投影值texp-water和texp-sample。
步骤2)利用MATLAB仿真出扫描区域为无气泡水域时的仿真走时投影值tsimulated。
步骤3)利用校准公式计算走时校准值trevised。
步骤4)将步骤3)所计算得到的走时校准值trevised对texp-sample进行校准得到经校准的走时投影值texp-sample-revised。
步骤5)将步骤4)中经校准修正的texp-sample-revised进行反演成像,得到目标的几何结构、坐标位置、数目和尺寸。
步骤6)结合几何结构、数目和尺寸,对含气率进行计算。
本发明的有益效果在于:本发明通过扇形束扫描方式,去除了伪影,提高了重建精度和分辨率,从而提高含气率的计算精度。在相同数量投影值的进行图像重建的结果,扇形束扫描方式明显优于平行束扫描方式,而且扇形束较平行束扫描方式效率高,操作较为简便和快捷。
附图说明
图1为超声环形阵示意图。
图2为本发明方法流程图。
图3为扇形束扫描坐标位置处于(-80,0)处的气泡重建结果图。
图4为平行束扫描坐标位置处于(-80,0)处的气泡重建结果图。
图5为扇形束扫描坐标位置处于(20,40)、(-20,-40)、(0,80)、(-80,0)和(0,0)处的五个气泡重建结果图。
图6为平行束扫描坐标位置处于(20,40)、(-20,-40)、(0,80)、(-80,0)和(0,0)处的五个气泡重建结果图。
具体实施方式
下面结合附图对本发明作进一步的说明。
针对平行束扫描方式所存在缺陷,采用扇形束扫描方式,设计了一阵超声环形阵,以获取扇形束投影值,弥补上述方案设计的不足,以期获得推广应用。
本发明提供了一种超声环形阵,以实现扇形束扫描,如图1所示,其直径210mm,在环形等弧长分布16个换能器,可作为发射端和接收端。
以气泡作为研究对象,使用超声环形阵,以合适频率的超声波作为激励信号,对目标进行扇形束扫描。具体操作步骤如下:选择16个阵元中的任意一个作为扫描源,发射超声波,经过扫描区域,由相应的阵元接收信号,依顺序重复执行此操作,直至16个阵元发射完毕截止。最后,在相应扫描方向的接收信号中提取走时作为投影数据。将投影数据记作Ti,j,其中i表示发射端序号,j表示接收端序号,共获取i×j投影值。将所获得的投影数据进行输入数据进行反演成像,对目标进行图像重建。
利用MATLAB对超声环形阵获取走时投影值进行仿真,表1所示为部分仿真结果和实验结果。
表1某一方向扇形束投影值
从表1可以看到,仿真结果和实验结果相差很大,这主要是在实验中由于水下环境复杂和环形阵本身存在系统性误差,导致实验结果准确性较差,但是对比两组走时数据,不难看到其变化趋势呈抛物状,因此,因此只要对实验结果进行校准,即可获得相对准确的实验走时投影数据,以供反演成像。其校准公式为:trevised=tsimulated-texp-water,其中trevised为走时校准值,tsimulated为仿真走时,texp-water为在无气泡水域中利用超声环形阵所获得的走时投影值。则校准过的实验走时投影值为texp-sample-revised=trevised+texp-sample,其中texp-sample-revised为经校准后的含气泡水域所获取的走时投影值,texp-revised为校准后的实验走时,texp-sample为校准前的含气泡水域所获取的走时投影值。
综上所述,该发明的实施步骤如下:1)选择合适频率的超声波作为激励信号,使用超声环形阵获取走时投影数据,使用时延去控制环形阵中每个阵元的发射时间顺序,分别获取扫描区域为无气泡水域和含气泡水域的走时投影值texp-water和texp-sample;2)利用MATLAB仿真出扫描区域为无气泡水域时的仿真走时投影值tsimulated;3)利用上述的校准公式计算走时校准值trevised;4)将3)所计算得到的走时校准值trevised对texp-sample进行校准得到经校准的走时投影值texp-sample-revised;5)将4)中经校准修正的texp-sample-revised进行反演成像,得到目标的几何结构、尺寸、坐标位置和数目,流程图如图2所示。
如图3、4、5和6所示,分别为扇形束和平行束走时投影利用上述校准公式经过校准后的投影值通过插值后,经过反演成像,重建气泡形状,计算出尺寸、数目和坐标位置。图3和图4分别为在扇形束和平行束扫描方式的重建结果,气泡位于坐标(-80,0),半径为20mm;通过对两种扫描方式成像结果进行对比,很明显地,扇形束扫描得到的成像结果在成像精度和分辨率上高于平行束扫描方式。通过对图5和图6进行相应计算,得到气泡的坐标位置分别为(20,40)、(-20,-40)、(0,80)、(-80,0)和(0,0),半径分别为5mm、10mm、15mm、20mm和25mm,结合气泡几何形状、数目、尺寸和实验时环形阵的尺寸参数等即可对含气率进行计算。然而,通过对成像结果进行分析,采用平行束扫描存方式的成像结果存在大量的伪影,使得重建精度和重建分辨率降低,而通过扇形束扫描方式,大大地克服了这一现象,大大地提高了重建精度和分辨率,从而提高气液两相流中含气率的精度至两个数量级。因此,在相同数量投影值进行图像重建的结果,扇形束扫描方式明显优于平行束扫描方式,而且扇形束较平行束扫描方式效率高,实验过程中操作较为简便和快捷。
Claims (5)
1.一种气液两相流含气率测量方法,其特征在于该方法包括以下步骤:
步骤1)选择超声波作为激励信号,使用超声环形阵获取走时投影数据,使用时延去控制环形阵中每个阵元的发射时间顺序,分别获取扫描区域为无气泡水域和向含气泡水域的走时投影值texp-water和texp-sample;
步骤2)利用MATLAB仿真出扫描区域为无气泡水域时的仿真走时投影值tsimulated;
步骤3)利用校准公式计算走时校准值trevised;
步骤4)将步骤3)所计算得到的走时校准值trevised对texp-sample进行校准得到经校准的走时投影值texp-sample-revised;
步骤5)将步骤4)中经校准修正的texp-sample-revised进行反演成像,得到目标的几何结构、坐标位置、树木和尺寸;
步骤6)结合几何结构、数目和尺寸,对含气率进行计算。
2.根据权利要求1所述的一种气液两相流含气率测量方法,其特征在于:所述的超声环形阵上的阵元等距离分布。
3.根据权利要求2所述的一种气液两相流含气率测量方法,其特征在于:超声环形阵获取走时投影数据的过程具体是:选择阵元中的任意一个作为扫描源,发射超声波,经过扫描区域,由相应的阵元接收信号,依顺序重复执行此操作,直至所有阵元发射完毕截止;最后,在相应扫描方向的接收信号中提取走时作为投影数据;将投影数据记作Ti,j,其中i表示发射端序号,j表示接收端序号,共获取i×j投影值。
4.根据权利要求1所述的一种气液两相流含气率测量方法,其特征在于:所述的校准公式为:trevised=tsimulated-texp-water。
5.根据权利要求4所述的一种气液两相流含气率测量方法,其特征在于:经校准的走时投影值texp-sample-revised=trevised+texp-sample。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610858540.2A CN106645396A (zh) | 2016-09-28 | 2016-09-28 | 一种气液两相流含气率测量方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610858540.2A CN106645396A (zh) | 2016-09-28 | 2016-09-28 | 一种气液两相流含气率测量方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN106645396A true CN106645396A (zh) | 2017-05-10 |
Family
ID=58854770
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201610858540.2A Pending CN106645396A (zh) | 2016-09-28 | 2016-09-28 | 一种气液两相流含气率测量方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106645396A (zh) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108490068A (zh) * | 2018-01-19 | 2018-09-04 | 天津大学 | 超声平面波扫描式多相流可视化测量装置 |
CN108536952A (zh) * | 2018-04-03 | 2018-09-14 | 东北大学 | 一种确定铁水包内气液两相流气含率的计算方法 |
Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1484821A (zh) * | 2001-01-05 | 2004-03-24 | 比约恩・A・J・安杰尔森 | 环形阵列 |
CN1718999A (zh) * | 2005-07-19 | 2006-01-11 | 大庆石油学院 | 一种气液两相流中含气率的检测方法 |
CN102305828A (zh) * | 2011-05-13 | 2012-01-04 | 中北大学 | 基于环绕式阵列的超声层析成像检测系统及其方法 |
CN103336145A (zh) * | 2013-06-18 | 2013-10-02 | 清华大学 | 一种针对管道内流体轴向流场的超声成像方法及设备 |
WO2014016110A1 (de) * | 2012-07-25 | 2014-01-30 | Siemens Aktiengesellschaft | Verfahren zur ermittlung einer blasengrössenverteilung und messeinrichtung |
CN104764804A (zh) * | 2015-03-16 | 2015-07-08 | 西安交通大学 | 超声Lamb波局部循环扫描概率重构层析成像方法 |
CN105044205A (zh) * | 2015-04-17 | 2015-11-11 | 北京理工大学 | 一种基于概率算法的气液两相流超声柔性阵列检测方法 |
CN105181793A (zh) * | 2015-05-14 | 2015-12-23 | 中国人民解放军国防科学技术大学 | 基于超声波双频信号测量两相流含气率的方法 |
-
2016
- 2016-09-28 CN CN201610858540.2A patent/CN106645396A/zh active Pending
Patent Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1484821A (zh) * | 2001-01-05 | 2004-03-24 | 比约恩・A・J・安杰尔森 | 环形阵列 |
CN1718999A (zh) * | 2005-07-19 | 2006-01-11 | 大庆石油学院 | 一种气液两相流中含气率的检测方法 |
CN102305828A (zh) * | 2011-05-13 | 2012-01-04 | 中北大学 | 基于环绕式阵列的超声层析成像检测系统及其方法 |
WO2014016110A1 (de) * | 2012-07-25 | 2014-01-30 | Siemens Aktiengesellschaft | Verfahren zur ermittlung einer blasengrössenverteilung und messeinrichtung |
CN103336145A (zh) * | 2013-06-18 | 2013-10-02 | 清华大学 | 一种针对管道内流体轴向流场的超声成像方法及设备 |
CN104764804A (zh) * | 2015-03-16 | 2015-07-08 | 西安交通大学 | 超声Lamb波局部循环扫描概率重构层析成像方法 |
CN105044205A (zh) * | 2015-04-17 | 2015-11-11 | 北京理工大学 | 一种基于概率算法的气液两相流超声柔性阵列检测方法 |
CN105181793A (zh) * | 2015-05-14 | 2015-12-23 | 中国人民解放军国防科学技术大学 | 基于超声波双频信号测量两相流含气率的方法 |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108490068A (zh) * | 2018-01-19 | 2018-09-04 | 天津大学 | 超声平面波扫描式多相流可视化测量装置 |
CN108536952A (zh) * | 2018-04-03 | 2018-09-14 | 东北大学 | 一种确定铁水包内气液两相流气含率的计算方法 |
CN108536952B (zh) * | 2018-04-03 | 2020-09-29 | 东北大学 | 一种确定铁水包内气液两相流气含率的计算方法 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN104501938B (zh) | 一种在置于空气中的矩形混响水槽内测量水下声源低频辐射声功率的方法 | |
CN101907704A (zh) | 多模式合成孔径雷达仿真成像评估方法 | |
CN102608205B (zh) | 基于变波速相位迁移的多层物体无损检测超声成像方法 | |
CN105866740B (zh) | 一种基于压缩感知的水声匹配场定位方法 | |
CN103454695B (zh) | 一种gps电离层tec层析方法 | |
CN104007444A (zh) | 一种基于中心投影的地面激光雷达反射强度图像的生成方法 | |
CN102305828A (zh) | 基于环绕式阵列的超声层析成像检测系统及其方法 | |
CN110261485A (zh) | 一种超声波测量材料内部各处弹性模量及泊松比的方法 | |
CN102590804A (zh) | 一种多普勒声纳陆上测试系统及其测试方法 | |
CN106645396A (zh) | 一种气液两相流含气率测量方法 | |
RU2507531C1 (ru) | Способ пассивного определения координат источников излучения | |
Nakahata et al. | Three-dimensional numerical modeling of ultrasonic wave propagation in concrete and its experimental validation | |
CN106525976A (zh) | 基于声发射层析成像定量分析混凝土结构损伤部位的方法 | |
CN103018333A (zh) | 分层物体的合成孔径聚焦超声成像方法 | |
CN108627841A (zh) | 一种浅滩掩埋管道超声探测中的圆心定位方法 | |
CN105157597A (zh) | 一种激光测量效应靶变形的方法 | |
CN106841382B (zh) | 基于三波耦合互作用非均匀混合介质非线性系数测量方法 | |
JP6621656B2 (ja) | 磁界のモデリング | |
CN110244305A (zh) | 一种水下目标信号散射的仿真方法 | |
CN113777654B (zh) | 一种基于伴随状态法初至波走时层析的海水速度建模方法 | |
CN109655525A (zh) | 一种用于超声层析成像系统的相对激励方法 | |
CN113805233B (zh) | 一种点扩散函数的计算方法 | |
Langener et al. | Simulation of full-angle ultrasound process tomography with two-phase media using a ray-tracing technique | |
CN109188378B (zh) | 一种超声多普勒回波信号仿真方法 | |
CN109858149B (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 | ||
WD01 | Invention patent application deemed withdrawn after publication | ||
WD01 | Invention patent application deemed withdrawn after publication |
Application publication date: 20170510 |