CN103591901B - 一种气液两相流气泡直径检测方法 - Google Patents

Abstract

本发明公开了一种气液两相流气泡直径检测方法，本发明的检测装置采用对射式光路，包括双路红外检测对管，与红外检测对管相连的信号调理模块包括红外发光二极管的脉冲恒流源驱动电路、红外探测器的光信号检测电路及红外发光二极管的恒温控制电路。单路红外检测对管实现泡状流条件下气泡大小的检测，检测电压峰值对应气泡直径大小。本发明可实现不同流型条件下的气液两相流气泡直径检测，具有测量精度高、响应时间快等优势。

Description

一种气液两相流气泡直径检测方法

技术领域

本发明涉及一种测量气液两相流气泡直径的方法，特别适用于密闭环境下电解制氧装置中循环水中气泡直径检测。

背景技术

对于载人空间站、潜艇等装备，一般都是采用电解水制氧方法来实现对乘员的供氧需求。SPE(solidpilymerelectrode)电解制氧技术是以固体聚合物为电解质的电解水制氧技术。它目前广泛应用在航天器、核潜艇等装备上，取代传统的碱性电解制氧装置。

在SPE电解制氧装置中，作为反应物的水经过净化，在电解槽、换热器、水箱之间循环。水气分离装置将反应产生的气/水混和物进行分离，分离出的水将重新投入循环。在此过程中，未电解的水需要循环利用，而返回的水通常带有一定量的氢气和氧气，如果气泡大小及体积含量超过一定比例，容易引起“气蚀”现象，导致电解制氧装置循环泵的损坏，这对电解制氧装置乃至整个装备都将造成极大的安全隐患，因此需要实时检测水中的气泡大小及体积含量，判定电解制氧装置是否正常工作，防止安全事故的发生。

目前，国内外检测液态水中气泡大小及体积含量比较常见的方法为医疗器械领域中采用超声波测量法。如专利《一种气泡检测装置》(申请号：CN102335476A)基本只能监测血液中有无气泡，而无法实现气泡大小的准确测量。专利《光气泡检测系统》(申请号：01823139.X)基于光的折射实现气泡大小的测量，其效果只能实现管道内有无气泡，而无法实现不同流型条件下气泡累计体积的计量。专利《在线气泡检测仪器》(授权号：CN2935140Y)，基于激光技术实现涂布流体在制备、处理和输送过程中微气泡的监测，该技术只能实现单一流型条件下的气泡测量，基于单传感器基本无法实现复杂流型条件下的气泡体积准确测量。

因此，需要一种能够在不同流型条件下实现气泡体积的测量装置及方法，以便对电解制氧装置的安全运行提供保障。

发明内容

本发明所要解决的技术问题是，针对现有技术不足，提供一种气液两相流气泡直径检测方法，在不同流型条件下准确测量气泡直径，为电解制氧装置及其他需要检测气泡大小的装置的安全运行提供保障。

为解决上述技术问题，本发明所采用的技术方案是：一种气液两相流气泡直径检测方法，包括气液两相流气泡直径检测装置，所述检测装置包括分离气液两相流的测量通道(1)，所述测量通道(1)的两个出口上各安装有红外检测对管，所述红外检测对管包括分别安装于所述测量通道(1)出口两侧且位置相对的红外发光二极管(2)与红外探测器(3)；所述红外发光二极管(2)和红外探测器(3)均与信号调理模块(8)电连接，所述信号调理模块(8)通过数据采集卡(13)接入计算机(14)，该方法为：

1)数据采集卡采集信号调理模块的输出电压信号以采样频率f，并对采集的电压信号进行平滑、去噪处理；

2)设定分析周期T，将分析周期T内信号调理模块的输出电压信号存入数组a[M]，其中M＝T·f；其中T的取值范围为0.5s～60s；

3)比较数据a[M]中输出电压信号与弹状流电压信号最大阈值V_slug-max的大小，若a[M]≥V_slug-max，则判定为弹状流气泡；比较数组a[M]中电压信号与泡状流电压信号最大阈值V_bubble-max、泡状流电压信号最小阈值V_bubble-mix的大小，若V_bubble-mix≤a[M]≤V_bubble-max，则判定为泡状流气泡；比较数组a[M]中电压信号与泡状流电压信号最大阈值V_bubble-max、弹状流电压信号最大阈值V_slug-max的大小，若V_bubble-max＜a[M]＜V_slug-max，则判定为环状流气泡；

4)弹状流气泡直径D_slug的大小与测量通道内径相同；

5)第i个泡状流气泡直径D_bubble[i]的计算公式为：

D_bubble[i]＝a₂₁·U[i]²+a₂₂·U[i]+a₂₃，

其中，a₂₁，a₂₂，a₂₃为二次拟合曲线系数，U[i]为分析周期T内数组a[M]中电压信号在V_bubble-mix～V_bubble-max之间的第i个脉冲信号对应的峰值；

6)环状流气泡直径D_annular的计算公式为：

D_annular＝a₃₁·V_annular-avr ²+a₃₂·V_annular-avr+a₃₃，

其中，a₃₁，a₃₂，a₃₃为二次拟合曲线系数，V_annular-avr为数组a[M]中大于泡状流电压信号最大阈值V_bubble-max，小于弹状流电压信号最大阈值V_slug-max的大小的所有电压数据的平均值。

与现有技术相比，本发明所具有的有益效果为：本发明可在不同流型条件下，实现气液两相流中气泡直径的测量，与超声波测量技术相比，具有气泡直径测量精度高的优势。

附图说明

图1为本发明一实施例检测装置结构示意图；

图2为本发明一实施例信号调理模块结构框图；

图3为本发明一实施例弹状流信号采集结果示意图；

图4为本发明一实施例泡状流信号采集结果示意图；

图5为本发明一实施例环状流信号采集结果示意图。

具体实施方式

如图1所示，本发明一实施例检测装置包括分离气液两相流的测量通道1，所述测量通道1的两个出口上各安装有红外检测对管，所述红外检测对管包括分别安装于所述测量通道1出口两侧且位置相对的红外发光二极管2与红外探测器3；所述红外发光二极管2和红外探测器3均与信号调理模块8电连接，所述信号调理模块8通过数据采集卡12接入计算机14。

如图2所示，信号调理模块8包括与所述红外探测器3电连接的光信号检测电路10、与所述红外发光二极管2电连接的恒温控制电路11、以及与所述光信号检测电路10、恒温控制电路11连接的脉冲恒流源驱动电路9。

本实施例中，测量通道为蓝宝石玻璃管，其内径D为8mm。

测量通道两个出口上的红外检测对管平行，且位置相对，即两个红外发光二极管的位置相对，两个红外探测器的位置相对；测量通道两个出口上的红外检测对管中心轴之间的距离L(即两个红外发光二极管中心轴或两个红外探测器中心轴之间的距离)的范围为20mm。

本实施例的气泡直径检测方法如下：

(1)数据采集及预处理：数据采集卡13对信号调理模块8的输出电压信号以采样频率f≥1KHz进行采集，并对采集的电压信号进行数据平滑、去噪处理。

(2)设定分析周期T，通过对单个分析周期内的气泡体积进行累加求和，最终实现气泡体积计量。将分析周期T内的采集数据存在数组a[M]，M＝T·f。比较数据a[M]中电压信号与弹状流电压信号阈值V_slug-max，若a[M]≥V_slug-max，则判别为弹状流。

对数组a[M]中大于弹状流电压信号阈值V_slug-max的所有数据求平均值，得到弹状流信号平均值V_Slug-avr，其与弹状流气柱直径大小D_slug的关系见式(1)。

D_Slug＝a₁₁·V_Slug-avr ²+a₁₂·V_Slug-avr+a₁₃(1)

其中：a₁₁，a₁₂，a₁₃为二次拟合曲线系数。累计统计信号大于V_slug-max的采集信号个数P。

弹状流气泡直径大小D_slug也可以根据测量通道内径确定，本发明中，D_slug与测量通道内径相同。

(3)比较数组a[M]中电压信号与泡状流电压信号最大阈值V_bubble-max及最小阈值V_bubble-mix。若V_bubble-mix≤a[M]≤V_bubble-max，则判定为泡状流。

通过峰值检测，统计电压信号在V_bubble-mix～V_bubble-max之间的脉冲个数及对应的脉冲信号峰值，存入数组b[N]，N为脉冲个数。气泡直径大小D_bubble与脉冲信号峰值电压U之间的关系见式(2)。

D_bubble＝a₂₁·U²+a₂₂·U+a₂₃(2)

其中：a₂₁，a₂₂，a₂₃为二次拟合曲线系数。

(4))比较数组a[M]中电压信号与泡状流电压信号最大阈值V_bubble-max及环状流阈值V_slug-max。若V_bubble-max＜a[M]＜V_slug-max，则判定为环状流。并累计相应的电压脉冲个数Q。

对数组a[M]中大于弹状流电压信号阈值V_slug-max，小于弹状流电压信号阈值V_slug-max的所有数据求平均值，得到环状流采集信号平均值V_annular-avr。环状流气柱直径大小D_annular与采集信号平均值V_annular-avr的关系见式(3)：

D_annular＝a₃₁·V_annular-avr ²+a₃₂·V_annular-avr+a₃₃(3)

气液两相流气泡体积检测装置测量的气泡直径可在计算机14上实时显示出来。

本实施例中，信号调理模块输出电压信号范围为0.5V～4.5V。0.5V对应管道内全水状态时红外探测器的信号输出值，4.5V对应管道内全气状态时红外探测器的信号输出值。

图3、图4、图5对应装置对环状流、弹状流、泡状流流型条件下的信号采集波形。本发明中，红外检测对管的输出电压范围为0.5～4.5V，弹状流气泡对应的电压大小≥4.4V，泡状流气泡对应的电压大小为0.5～3.6V，环状流气泡对应的电压大小为3.6～4.4V。其中，弹状流电压信号阈值V_slug-max为4.4V，泡状流电压信号最大阈值V_bubble-max为3.6V，最小阈值V_bubble-mix为0.6V。

其中：

弹状流气泡直径为8mm；

在标准气泡信号源发生装置上，选取三个或三个以上气泡直径标定点(选取的标定点越多，拟合误差越小)，本实施例中选取了直径为2mm、4mm、6mm的泡状流气泡，其脉冲峰值电压信号分别为0.8V，1.5V，3.6V。采用二阶函数拟合得到如下关系：

D_bubble[i]＝-0.80·U[i]²+5.03·U[i]-1.83

其中：a₁₁＝-0.8，a₁₂＝5.03，a₁₃＝-1.83。

环状流气泡直径与峰值电压的拟合曲线计算过程如下：

在标准气泡信号源发生装置上，选取三个或三个以上标定点(选取的标定点越多，拟合误差越小)，本实施例中选取了平均直径为6.5mm，7mm，7.5mm的环状流气泡，其信号平均值输出为3.81V，4.12V，4.38V。采用二阶函数拟合得到如下关系：

D_annular＝-0.62·V_annular-avr ²+6.82·V_annular-avr-10.48

其中：a₂₁＝-0.62，a₁₂＝6.82，a₁₃＝-10.48。

对采集信号按照本发明的方法流程进行数据处理，得到气液两相流气泡直径的测量结果。

Claims (1)

1.一种气液两相流气泡直径检测方法，包括气液两相流气泡直径检测装置，所述检测装置包括测量气液两相流气泡直径的测量通道(1)，所述测量通道(1)的两个出口上各安装有红外检测对管，所述红外检测对管包括分别安装于所述测量通道(1)出口两侧且位置相对的红外发光二极管(2)与红外探测器(3)；所述红外发光二极管(2)和红外探测器(3)均与信号调理模块(8)电连接，所述信号调理模块(8)通过数据采集卡(13)接入计算机(14)，其特征在于，该方法为：

1)数据采集卡采集信号调理模块的输出电压信号以采样频率f，并对采集的电压信号进行平滑、去噪处理；

2)设定分析周期T，将分析周期T内信号调理模块的输出电压信号存入数组a[M]，其中M＝T·f；其中T的取值范围为0.5s～60s；

3)比较数据a[M]中输出电压信号与弹状流电压信号最大阈值V_slug-max的大小，若a[M]≥V_slug-max，则判定为弹状流气泡；比较数组a[M]中电压信号与泡状流电压信号最大阈值V_bubble-max、泡状流电压信号最小阈值V_bubble-mix的大小，若V_bubble-mix≤a[M]≤V_bubble-max，则判定为泡状流气泡；比较数组a[M]中电压信号与泡状流电压信号最大阈值V_bubble-max、弹状流电压信号最大阈值V_slug-max的大小，若V_bubble-max＜a[M]＜V_slug-max，则判定为环状流气泡；

4)弹状流气泡直径D_slug的大小与测量通道内径相同；

5)第i个泡状流气泡直径D_bubble[i]的计算公式为：

D_bubble[i]＝a₂₁·U[i]²+a₂₂·U[i]+a₂₃，

其中，a₂₁，a₂₂，a₂₃为二次拟合曲线系数，U[i]为分析周期T内数组a[M]中电压信号在V_bubble-mix～V_bubble-max之间的第i个脉冲信号对应的峰值；

6)环状流气泡直径D_annular的计算公式为：

D_annular＝a₃₁·V_annular-avr ²+a₃₂·V_annular-avr+a₃₃，

其中，a₃₁，a₃₂，a₃₃为二次拟合曲线系数，V_annular-avr为数组a[M]中大于泡状流电压信号最大阈值V_bubble-max，小于弹状流电压信号最大阈值V_slug-max的大小的所有电压数据的平均值。