一种气液两相流气泡体积检测装置及检测方法

摘要

本发明公开了一种气液两相流气泡体积检测装置及检测方法，本发明的检测装置采用对射式光路，包括双路红外检测对管，与红外检测对管相连的信号调理模块包括红外发光二极管的脉冲恒流源驱动电路、红外探测器的光信号检测电路及红外发光二极管的恒温控制电路。单路红外检测对管实现泡状流条件下气泡大小的检测，通过对累计脉冲个数统计实现气泡体积计量；计算同一气泡流经固定距离内的二对红外发光二极管所需时间，实现弹状流、环状流的气泡流速测量，并结合气泡大小，实现气泡体积的测量。本发明可实现不同流型条件下的气液两相流气泡体积检测，具有测量精度高、响应时间快及体积小等优势。

说明书

技术领域

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

背景技术

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

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

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

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

发明内容

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

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

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

所述测量通道两个出口上的两个红外发光二极管的波长范围分别为 1.4μm±0.1μm、1.9μm±0.1μm。

所述测量通道两个出口上的红外检测对管平行，且所述两个出口上的红 外检测对管中心轴之间的距离L的范围为6mm～50mm。

所述测量通道为蓝宝石玻璃管，所述蓝宝石玻璃管的内径D的范围为 1mm～20mm。

本发明利用所述检测装置检测气液两相流气泡体积的方法为：

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）弹状流气泡体积V_slug的计算公式为：

V_slug=π·（D_slug/2）²·v_slug·P·（1/f），

其中，D_slug为弹状流气泡直径，D_slug大小与测量通道内径相同，P为数组a[M] 中大于弹状流电压信号阈值V_slug-max的电压脉冲的个数；v_slug为弹状流气泡的 流速，v_slug=L/t，L为测量通道两个出口上的红外检测对管中心轴之间的距离， t为测量通道两个出口上红外检测对管之间同一弹状流的电压信号序列的峰值 间隔时间；

5）泡状流气泡体积V_bubble的计算公式为：

其中，D_bubble[i]为第i个泡状流气泡直径，N为一个分析周期T内累计的泡状 流气泡个数，气泡个数与分析周期T内V_bubble-mix≤a[M]≤V_bubble-max的a[M]中 的电压脉冲信号个数相同；

6）环状流气泡体积V_annular的计算公式为：

V_annular=π·（D_annular/2）²·Q·（1/f），

其中，D_annular为环状流气泡直径，Q为数组a[M]中，满足V_bubble-max＜a[M]＜ V_slug-max的电压脉冲的个数；

7）计算整个采用周期内的气液两相流气泡体积V：

V=∑（V_bubble+V_slug+V_annular）。

第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为本发明一实施例环状流信号采集结果示意图；

图6为本发明一实施例方法流程图。

具体实施方式

如图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）

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

其中：a₁₁，a₁₂，a₁₃为二次拟合曲线系数。累计统计信号大于V_slug-max的 采集信号个数P。比较脉冲信号峰值，计算出测量通道两个出口上的红外检测 对管之间同一弹状流的电压信号序列的峰值间隔时间t，并通过两个出口上红 外检测对管中心轴之间的固定距离L，计算出弹状流气泡的流速v_slug=L/t，弹 状流气泡体积V_slug可根据式（2）计算：

V_slug=π·（D_slug/2）²·v_slug·P·（1/f）   （2）

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

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

D_bubble[i]=a₂₁·U[i]²+a₂₂·U[i]+a₂₃   （3）

其中：a₂₁，a₂₂，a₂₃为二次拟合曲线系数。并根据式（4），计算泡状流气 泡体积V_bubble。

（4））比较数组a[M]中电压信号与泡状流电压信号最大阈值V_bubble-max及 环状流阈值V_slug-max。若V_bubble-max＜a[M]＜V_slug-max，则判定为环状流。并累计 满足U_bubble-max＜a[M]＜U_slug-max的a[M]中的电压脉冲数量Q。

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

D_annular=a₃₁·V_annular-avr²+a₃₂·V_annular-avr+a₃₃   （5）

其中：a₃₁，a₃₂，a₃₃为二次拟合曲线系数。环状流气泡体积V_annular可根据 式（6）计算：

V_annular=π·（D_annular/2）²·Q·（1/f）   （6）

（5）整个采样周期内，气泡体积V可通过对三种流型条件下的气泡体 积进行累加求和。

V=∑（V_bubble+V_slug+V_annular）   （7）

气液两相流气泡体积检测装置测量的气泡体积可在计算机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。

对采集信号按照图6所示流程图进行数据处理，得到气液两相流气泡大 小及体积的测量结果。