CN104101687B - 基于超声多普勒与电学多传感器的多相流可视化测试方法 - Google Patents

Abstract

本发明属于流体测量技术领域，涉及一种基于电学与超声双模态传感器阵列的多相流可视化测量方法，该测量方法包括以下步骤：电学传感器阵列采集电信号；超声传感器阵列采集超声信号；利用电学信号与超声信号提取被测对象流动特征实现流型识别；判断连续相导电性；连续相导电选择电导率计算水相含率并重建水相分布；连续相不导电选择介电常数计算水相含率并重建水相分布；利用超声透射信号进行气相分布重建；利用超声多普勒信号计算液相与水相流速；计算各分相流速；融合电学重建结果与超声重建结果获得油气水多相流相分布。本发明的测量装置具有适用范围宽、可视化、成本低、信息丰富等特点。

Description

基于超声多普勒与电学多传感器的多相流可视化测试方法

技术领域

本发明属于流体测量技术领域，涉及一种电学与超声双模态多传感器多相流过程参数测试方法，用于实现多相流分相含率、流速的在线获取。本发明以油/气/水多相流测量为描述对象，但并不仅限于该应用，在其他工业过程和化学反应中存在的多相流动现象中，本发明的测试方法仍适用。

技术背景

在工业生产与日常生活中广泛存在着多相流动现象。多相流中的“相”定义为物质的存在形式，即气态、液态或固态，因此多相流即为具有两种或两种以上“相”物质同时流动的流体。多相流现象广泛存在于能源、动力、石油、化工、冶金、医药等工业过程中，在工业生产与科学研究中有着十分重要的作用，对其流动过程监测和描述、以及对流动过程参数的准确检测也给工程师和科研人员提出挑战。近年来，国际上对多相流的研究兴趣在持续增长，其原因在于多相流不仅在一系列现代工程中得到广泛应用，而且对促进这些工程设备的发展和创新也起到了重要作用。

由于多相流各相间存在界面效应和相对速度，相界面在时间和空间上均呈随机变化，致使多相流的流动特性远比单相流复杂，特征参数也比单相流多。多相流过程参数检测策略随着工况与对象属性变化，可以利用的物理现象与关系很多，因此检测手段也多种多样。检测方式可分为两大类：直接法与间接法。直接法指的是对象参数能通过测量直接得到，而间接法需要在测量值(辅助测量值)与被测参数间建立关系式通过计算得到，即“软测量”方法。在多相流测量领域内，许多直接采用单相流仪表的检测方法属于直接法，而新兴的检测技术则多采用软测量的方法，如电学法、超声法以及射线法等。在科学研究与工业应用中，多相流的检测手段需要对被测流体不产生任何扰动，因此电学法与超声法由于其结构简单、非扰动、造价低而备受关注。

电学测量方法根据其传感器结构、形状、激励方式的不同又分为多种形式，如电导探针、环形电导阵列等，其中旋转场电阻抗方法即电学过程层析成像技术是电学法中的一种重要改进形式，能够提供丰富的物质截面分布信息，且能将不透明管道中的多相介质分布进行可视化重建。此外，该技术具有非扰动与多点测量的特点，可以实现被测多相流体的相含率与流速等过程参数，具有很好的应用前景。电学过程层析成像技术根据测量敏感原理不同又分为电容式、电阻式、阻抗式和电磁式层析成像等。

超声检测是一种应用较为广泛的技术，在医学监测、流体测量中有其独特的优势，超声波在流体中传播时不会破坏流体的流场，没有压力损失，同时若将检测元件置于管道外壁，可以避免与流体直接接触，降低传感器的腐蚀程度。超声层析成像方法可通过安装于管道同一截面处的多个超声收发探头，以非扰动的形式获得被测截面内部不同声阻抗介质的分布信息，并通过图像重建算法实现该分布的可视化。同时，当入射超声波遇到离散相液滴和气泡时会发生反射作用，并在超声接收端产生与反射物体(离散相液滴和气泡)流速有关的声波频率差，即超声多普勒效应。利用该效应可有效获得油气水多相流的流速分布信息。

发明内容

本发明的目的是提供一种精确地、非扰动的多相流过程参数可视化测试方法。本发明的传感器和测量方法，将电学传感器阵列与超声传感器相组合提取多相流流动信息，应用信息融合技术对多相流的多源信息进行处理，不对多相流动产生扰动且无需对测量流体进行预分离或混合，实现多相流过程分相含率与分相流速的准确获取，流型的正确识别以及相分布的可视化重建。

本发明的技术方案如下：

一种基于超声多普勒与电学多传感器的多相流可视化测试方法，所采用的测试装置包括电学传感器阵列、超声传感器阵列、电学信号发生与检测单元、超声信号发生与检测单元、流动参数计算与可视化单元；所述的电学传感器阵列包括两组分布在被测管道相同截面位置的电极阵列；所述的超声传感器阵列包括两组分布在被测管道不同截面位置的探头阵列，每组探头阵列由两个以上的超声探头构成，其中第一组以超声多普勒方式工作，包括高频超声探头与低频超声探头，前者用于发射高频超声信号，以获取被测多相流体内的液滴流速，后者用于发射低频的超声信号，以获取被测多相流体内的气泡流速，第二组以超声透射法工作，通过超声幅值衰减获取被测多相流体相含率；所述的电学传感器阵列与超声传感器阵列组成一套双模态传感器阵列；电学信号发生与检测单元选通电学传感器阵列中的一对电极作为激励电极对，其中一个电极为激励电极，另一个电极为接地电极，在所述激励电极对之间建立电学敏感场，测量其余电极间电势差；对于第一组超声探头，超声信号发生与检测单元选通该组探头阵列中的一个探头发出超声波，该探头自身或另一个探头接收超声波，将超声强度转换为电信号并提取频移信息；对于第二组超声探头，超声信号发生与检测单元选通该组探头阵列中的一个探头发出超声波，其余探头中的一个或多个探头接收超声波，并将超声强度转换为电信号并提取幅值信息；电学信号发生与检测单元与超声信号发生与检测单元所获信号一起，送入流动参数计算与可视化单元，实现流型在线识别、相含率与分相流速的计算，以及相分布的可视化重建与显示，该测试方法包括下列步骤：

1).采集电学传感器阵列信号与超声传感器阵列信号；

2).根据电学传感器阵列测量数据和超声传感器阵列测量数据，结合特征提取与分类算法，实现多相流流型的识别；

3).根据超声透射衰减幅度，利用超声传感器测量数据带入计算气相含率α_g；其中U为超声穿透多相流混合物后的幅值，U₀为超声穿透单相水后的幅值，L为超声传播路径长度，X_g为衰减系数，通过实验标定获得；

4).采用多频率超声激励，超声探头测量端获得的超声频率f'与原始发射频率f间存在关系式其中u为超声波在被测介质中的行进速度，u_s为液滴或气泡的移动速度，当液滴或气泡流向测量端时，u_s前方运算符号为-号，反之则为+号；根据将高频超声波反射频率变化带入上式获取液滴流速u_o＝u_s，将低频超神波反射频率变化带入上式获取气泡速度u_g＝u_s，假设油相与水相间不存在速度差，因此水相流速u_w＝u_o。多频超声高频频率范围典型取值为4MHz以上，低频频率范围典型取值为500KHz-1MHz，以上范围根据被测流体密度与流体内声速变化；

5).根据流型识别结果选取相含率估计模型，并判断连续相是否导电，

若连续相为非导电相，采用电容传感器测量数据，计算模型如下：

a.根据电容电极对采集的电压信号，计算多相流体混合介电常数ε_m；

b.已知油相介电常数ε_o、水相介电常数ε_w及气相介电常数ε_g，根据公式

$\sqrt{{\mathrm{\ϵ}}_m}={\mathrm{\α}}_g\sqrt{{\mathrm{\ϵ}}_g}+(1-{\mathrm{\α}}_g)({\mathrm{\α}}_o\sqrt{{\mathrm{\ϵ}}_o}+{\mathrm{\α}}_w\sqrt{{\mathrm{\ϵ}}_w})$ 以及多相流分相含率和为1的关系式

α_o+α_g+α_w＝1，计算油相含率α_o、气相含率α_g与水相含率α_w。

若连续相导电，采用电导传感器测量数据，计算模型如下：

a.根据电阻电极对采集的电压信号，计算多相流体混合电导率σ_m；

b.已知油相电导率σ_o、气相电导率σ_g与水相电导率σ_w，另α_l＝α_o+α_w为油气水多相流的液相含率，根据公式和公式α_o+α_g+α_w＝1计算出油相含率α_o、水相含率α_w和气相含率α_g，式中σ_m为多相混合流体电导率，σ_l＝α_oσ_o+α_wσ_w为液相混合电导率；

6).利用所述电学传感器阵列与超声传感器阵列测量数据进行油气水多相流相分布重建：

a.根据流型识别结果判断连续相是否导电，若连续相导电，将电学传感器阵列视为电导传感器阵列，设采用电导传感器测量数据结合图像重建算法实现电学传感器截面内电导率分布

σ(x,y)，利用二值化方法获得水相分布f_w(x,y)，其中x与y为像素坐标，若像素(x,y)

为水相，则其值为1，若像素(x,y)为气相或油相，则该值为0；

b.若连续相为非导电相，将电学传感器阵列视为电容传感器阵列，采用电容传感器测量数据结合图像重建算法实现电学传感器截面内的介电常数分布ε(x,y)，利用二值化方法获得水相分布f_w(x,y)，若像素(x,y)为水相，则其值为1，若像素(x,y)为气相或油相，则该值为0；

c.利用超声传感器阵列测量数据结合图像重建算法实现超声传感器被测截面内的声阻抗分布Z(x,y)，利用二值化方法获得气相分布f_g(x,y)，若像素(x,y)为气相，则其值为1，若像素(x,y)为水相或油相，则该值为0；

d.设f_A(x,y)为截面内的相分布，且各像素值均为1，由于油、气、水三相充满管道截面，因此可得油相分布f_o(x,y)＝f_A(x,y)-f_w(x,y)-f_g(x,y)，其中，若像素(x,y)为油相，则其值为1，若像素(x,y)为水相或气相，则该值为0。

基于电学敏感原理的信息检测技术对被测流体的电学性质变化敏感，如电导率或介电常数，但对密度差较大的气液相间边界信息模糊，不易获得。而多相流中，由于气相与液相界面处声阻抗的巨大差异，使得超声在气液界面的反射特性极为明显(可达99％)，因此超声对气液分界面有极好的分辨能力。同时，超声波在液滴或气泡的界面处产生反射波，该反射波的频率变化与反射体的速度有关，即超声多普勒效应，可实现分散相(液滴或气泡)流速的测量。因此，将基于电学敏感原理与基于超声敏感原理的传感器阵列相结合，利用电学传感器阵列获取液相流动信息，利用超声传感器阵列获取气相流动信息，可准确、全面地实现多相流流型在线识别、相含率和相流速在线估计，以及相分布的可视化重建。本发明利用不同粒径的两相流对不同频率超声多普勒的频移调制作用不同的原理获取气相与液相流速，利用超声衰减法获取气相分布计算气相含率，利用电导/电容法获取电导率/介电常数分布(水相分布)与含率，综合利用超声与电学两类传感器的测量信息实现多相分布重建，进行多相流流型的在线识别。本发明的有益效果及优点如下：

1、电学测量法对多相流体的导电性有很好的区分作用，也即对油气水混合物中的水相变化敏感。而超声对多相流体的密度变化敏感，两类测量模态有互补性；

2、超声多普勒法可直接获得被测流体的流速信息；

3、该方法为非侵入的测量手段，不会对流体产生任何的扰动；

4、多传感器阵列可获得对多相流动过程更全面的描述；

5、测量速度快，成本低；

6、可提供被测流体瞬时相分布状态的图像，并对过程参数的测量结果进行可视化显示。

附图说明

以下图描述了本发明所选择的实施例，均为示例性图而非穷举或限制性，其中：

图1本发明测试方法所用测试装置整体结构图；其中，0-来流方向；1-管道；2-电学传感器阵列，包括电容式电极与电导式电极的单独或组合使用，以及同时实现电容与电导检测的电极；3-电学信号发生与采集单元；4-电学传感器阵列测量数据相含率提取单元；5-超声传感器阵列；6-超声信号发生与采集单元；7-超声传感器阵列测量数据多普勒测速与相含率提取单元；8-流动参数计算与可视化单元。

图2本发明测试方法所用测试装置电学传感器阵列结构图，其中图2(a)为传感器阵列结构侧视图，包括被测管道1，安装于不同截面位置的弧形电极阵列2；图2(b)为纵向A-A截面剖视图；图2(c)为横向B-B截面剖视图；

图3本发明测试方法所用测试装置超声传感器阵列结构图，其中图3(a)为传感器阵列结构侧视图，包括被测管道1与超声探头阵列(5a，5b)；图3(b)为纵向B-B截面剖视图；图3(c)为横向A-A截面剖视图；

图4本发明测试方法所用测试装置的超声与电学双模态测量系统结构图；

图5本发明测试方法所用测试装置电学信号发生与检测单元结构图；

图6本发明测试方法所用测试装置超声信号发生与检测单元结构图；

图7本发明测试方法所用测试装置计算中心功能结构图；

图8本发明测试方法计算步骤；

图9本发明测试方法图像重建步骤；

具体实施方式

以下详细描述制造和操作本发明的步骤，旨在作为本发明的实施例描述，并非是可被制造或利用的唯一形式，对其他可实现相同功能的实施例也应包括在本发明的范围内。本发明中的电学测量包括电阻、电容或电磁信息的检测，既可同时使用几种电学测量方式亦可分别使用。

下面结合说明书图详细说明本发明的实施例。

图1描述了本发明测试方法所用测试装置的整体结构图，包括一个被测管段1，一组电学传感器阵列2以及与其相连的电学信号发生于采集单元3，一组超声传感器阵列5以及超声信号发生与采集单元6。电学传感器阵列2包括一组或一组以上由多个分布在被测管道相同截面位置的电极构成；超声传感器阵列5包括两组或两组以上以一定间隔分布在被测管道上的超声探头阵列，每组探头由多个分布在被测管道相同截面位置的探头构成，其中至少一组以超声多普勒方式工作，获取被测多相流体流速，另外至少一组以超声透射法工作，获取被测多相流体相含率；所述的电学传感器阵列2与超声传感器阵列5在被测管道上同时安装(无安装顺序要求，可交叉位置安装)、同时工作，组成一套双模态传感器阵列，且不对多相流动造成任何干扰，因此所述的两类传感器阵列在管道上下游的安装顺序不影响测量效果。

当被测多相流从来流方向0进入被测管段时，电学传感器阵列2可通过电学信号发生与采集单元3获得包含被测流体水相含率波动信息的测量数据，超声传感器阵列5可通过超声信号发生与采集单元6获得被测流体气相含率波动以及流速的测量数据。将以上信息同时送入流动参数计算与可视化单元8进行流型识别，判断连续相是否导电，并依据判别结果选择使用电学传感器阵列2中的电容传感器或电导传感器。之后将电学信号发生与采集单元3获得的电容或电导测量数据送入电学传感器相含率计算单元4中进行水相含率计算，同时将超声信号发生与采集单元6获得测量数据送入超声传感器多普勒测速与相含率计算单元7中进行流速与气相含率的计算。依据流型识别结果选择计算模型，对所述的电学传感器相含率计算单元4与超声传感器多普勒测速与相含率计算单元7的计算结果进行处理，获得混合与分相流速。利用电学传感器阵列测量数据与超声传感器阵列测量数据进行信息融合，实现多相流相分布的重建与可视化显示。

图2为本发明测试方法所用测试装置的电学传感器阵列结构图，传感器阵列为安装在被测管道1的弧形电极阵列2组成，包含安装在被测管道同一截面位置内的两个及以上弧形电极。弧形电极尺寸与数目视应用条件可以改变。工作时可在任意弧形电极上施加激励信号，任意其他电极接地形成激励电极对。在该激励电极对中间形成电学测量敏感场，当多相流体流过该敏感场时，由于场内电学参数随着多相介质的含率与分布变化，导致激励电极对间，以及其他电极之间的电势差发生相应的变化，通过测量该电势差可实现多相流参数的测量。本发明中的电学测量包括电阻、电容或电磁信息的检测，既可同时使用几种电学测量方式亦可分别使用，可在同一截面内通过安装不同尺寸的电极实现同一截面内信息的多样化获取，故本装置的电极结构具有多种形式。本发明装置的电学传感器阵列可使用不同频率的激励信号，通过获得被测多相流体在不同激励信号频率下的响应数据实现多相流信息更全面的提取。

图3为本发明测试方法所用测试装置的超声传感器阵列，包括多套安装在被测管道不同截面位置的超声探头阵列(5a，5b)组成，每套探头阵列包含安装在被测管道1同一截面位置内的两个及两个以上超声探头，其中一组以超声多普勒方式工作，获取被测多相流体流速，另外一组以超声透射法工作，获取被测多相流体相含率。超声探头尺寸与数目视应用条件可以改变。工作时可在同一截面内的任意探头上施加激励信号，任一其他探头接收超声信号，成激励电极对。在该激励电极对中间形成超声调制测量敏感场，当多相流体流过该敏感场时，由于场内介质密度随着多相介质的含率与分布变化，导致超声调制测量敏感场内超声强度发生相应的变化，通过声电转换装置(如压电陶瓷等)获取该超声信息，进而实现超声测量通路内相含率的测量，以及气相分布的重建。

由于不同尺寸的颗粒对不同频率的超声入射波所产生的多普勒频移调制特性不同，且油气水多相流动过程中，液相中的油滴(或水滴)的尺寸要明显小于气泡，且超声在油水界面与超声在气液界面处的反射率不同，因此本发明装置的超声传感器阵列可使用不同频率的激励信号，通过采集不同频率超声在不同强度和时间尺度上的频移实现多相流中不同尺寸液滴与气泡流速的测量，进而结合分相含率计算结果实现各分相流速的测量。

图4描述了本发明测试方法所用测试装置的超声与电学双模态测量系统结构，包括电学传感器阵列激励与检测单元，超声传感器阵列激励与检测单元，计算机总线以及计算中心，在同一总线平台下实现双模态传感器的测量信息获取。电学传感器阵列激励与检测单元用于在被测管道1中建立电学敏感场，进而获得多相流含水率信息，可实现电学与电导不同模态传感器信息的获取，具体结构如图5所示；超声传感器阵列激励与检测单元用于向被测管道1中的多相流体发射超声波，进而获得多相流含气率与流速信息，具体结构如图6所示。电学传感器阵列激励与检测单元，以及超声传感器阵列激励与检测单元所获得的测量数据通过计算机总线送入计算中心实施流型识别、相含率和流速的计算以及流动相分布的重建。

图5描述了电学信号发生与检测单元结构。可将电容与电导两种电学模式传感器的激励与信号采集在同一测量系统中实现。系统控制与设置信息由计算机通过计算机总线传送至逻辑控制单元，并通过逻辑控制单元对系统整体工作逻辑和参数进行控制和设定，根据系统设定要求在激励信号发生模块中实现激励信号的产生以及信号驱动能力的提升，即采用电压控电压源(VCVS)或电压控电流源(VCCS)将原始激励信号转换为幅值和相位可调的恒定电压或电流信号，并通过逻辑控制按一定规律选通电容/电导电极阵列的对应电极，使激励信号施加至测量空间中形成电学敏感场。当多相流体流过所述电学敏感场时，由于流体的电学参数变化导致电场强度分布发生变化，进而在电极阵列中会获得不同的电势差，利用信号解调模块与参考信号将该电势差进行解调处理，提取出包含多相流分相含率变化的模拟信号，并送入模拟-数字转换(A/D转换)阵列实现模拟信号的数字化，再通过计算机总线将测量数据送至计算中心实现流动参数的计算与电学参数分布的可视化重建。

图6描述了超声信号发生与检测单元结构。系统控制与设置信息由计算机通过计算机总线传送至逻辑控制单元，并通过逻辑控制单元对系统整体工作逻辑和参数进行控制和设定，根据系统设定要求在激励信号发生模块中实现激励信号的产生。通过逻辑控制单元按一定规律选通超声传感器阵列的对应探头，使激励信号通过电声转换功能产生超声波，当多相流体流过所述超声传播路径时，由于气液两相密度差会在气液界面处将一部分超声反射，并有一部分未反射超声波被接收探头接收，透射部分超声幅值的衰减幅度与超声传播路径中的气相尺寸有直接关系，反射波的频率变化与流速有直接关系，并通过差分放大与滤波电路对信号进行整形后通过A/D转换阵列转为数字信号，通过计算机总线送入计算中心，利用图像重建算法与流动参数计算模型实现气相分布图像重建与流速参数提取。

图7描述了计算中心的功能结构，包括电学传感器阵列与超声传感器阵列测量数据的相含率提取功能，超声多普勒测速功能，以及流动过程相分布重建与参数可视化功能。电学传感器阵列与超声传感器阵列的测量数据分别送入各自的相含率计算单元实现相含率信息的获取，超声传感器阵列利用所测反射波的频率变化实现流速计算，同时对电学传感器阵列与超声传感器阵列的测量数据进行信息融合与相分布重建。利用电学传感器阵列与超声传感器阵列所获得的流速与相含率信息，与重建的相分布结果一起送入参数计算与可视化单元实现最终结果的汇总与可视化输出。

图8描述了本发明测试方法的计算步骤，具体计算过程如下：

(一).采集电学传感器阵列信号与超声传感器阵列信号；

(二).根据电学传感器阵列测量数据和超声传感器阵列测量数据，结合特征提取与分类算法，实现多相流流型的识别。

(三).流型识别的一般步骤为：首先从测量数据中提取能够反映流型变化的特征值，如统计方法、时频分析方法、非线性分析方法等，然后利用分类算法将特征值归类识别，比如支持向量机、人工神经网络等。举例实现过程可参考专利[1]，不同之处在于本发明分别从电学传感器阵列与超声传感器阵列中的每一传感器输出中提取一定时间内的采集数据组成时间序列，从每一时间序列中提取特征组成识别特征向量，对电学传感器阵列与超声传感器阵列采集到的数据进行分类，最终可得到连续相是否导电以及典型流型识别结果。

(四).当超声波在介质中传输时发生衰减，其衰减的幅度与介质的声阻抗有直接关系。由于超声气液之间声阻抗的巨大差异，使得超声在气液界面的反射特性极为明显(可达99％)，采用穿透法测量时，由于气相的反射作用，使得接收到的信号幅值可以表征流动过程的气相含率。

将超声传感器阵列采集到的测量数据带入计算气相含率α_g。其中U为超声穿透多相流混合物后的幅值，U₀为超声穿透单相水后的幅值，L为超声传播路径长度，X_g为衰减系数，由实验标定获得。

(五).利用超声多普勒传感器阵列测量数据，实现多相流速度测量u_g。

当油气水多相流中的液滴或气泡流经超声多普勒测速的敏感空间中时，固定频率的超声波作用在所述液滴或气泡上形成反射，在接收端接收到的信号频率变化与流速有直接关系。根据多普勒效应，超声波在液滴或气泡与测量方向相对移动时接收的反射声波频率变高，而在与测量方向相反移动时接收的反射声波频率变低。即固定位置测量端获得的超声频率f'与原始发射频率f间存在以下关系式

其中u为超声波在被测介质中的行进速度，u_s为液滴或气泡的移动速度，当液滴或气泡流向测量端时，u_s前方运算符号为-号，反之则为+号。

由于油气水多相流中的液滴尺寸小于气泡尺寸，且不同波长的超声会对不同尺度的液滴或气泡有不同的响应特性，因此采用多频率超声激励，分别获取液滴流速u_o或气泡速度u_g，假设油相与水相间不存在速度差，因此水相流速u_w＝u_o。

(六).根据流型识别结果选取相含率估计模型，并判断连续相是否导电。

若连续相为非导电介质，如油相与气相为连续相时，采用电容传感器测量数据。此时多相流体的混合介电常数率ε_m可由公式计算，其中ε_o为油相介电常数、ε_w为水相介电常数，ε_g为气相介电常数,α_o、α_g与α_w分别为油相、气相与水相含率。将该公式与多相流分相含率和为1的关系式α_o+α_g+α_w＝1联立，带入利用超声传感器阵列测得的气相含率α_g，计算油相含率α_o与水相含率α_w。

若连续相为导电介质，如矿化水时，采用电导传感器测量数据。此时多相流体的混合电导率σ_m可由公式计算，其中σ_l＝α_oσ_o+α_wσ_w为液相混合电导率，σ_g与σ_w分别为油相、气相与水相电导率。将该公式多相流分相含率和为1的关系式α_o+α_g+α_w＝1联立，带入利用超声传感器阵列测得的气相含率α_g，计算油相含率α_o与水相含率α_w。

(七).图9描述了本发明测试方法的相分布重建步骤。基于单一模态电学传感器阵列测量数据或超声传感器阵列测量数据的相分布重建算法已有很多，如线性反投影等，因此在本发明测量方法中可采用任意图像重建算法实现单一测量模态的相分布图像重建，利用电阻数据进行图像重建的举例实现过程可参考文献[2]，利用电容数据进行图像重建的举例实现过程可参考文献[3]，利用超声数据进行图像重建的举例实现过程可参考文献[4]。

a.根据流型识别结果判断连续相是否导电，若连续相导电，将电学传感器阵列视为电导传感器阵列，设采用电导传感器测量数据结合图像重建算法实现电学传感器截面内电导率分布σ(x,y)，利用二值化方法获得水相分布f_w(x,y)，其中x与y为像素坐标，若像素(x,y)为水相，则其值为1，若像素(x,y)为气相或油相，则该值为0；

b.若连续相为非导电相，将电学传感器阵列视为电容传感器阵列，采用电容传感器测量数据结合图像重建算法实现电学传感器截面内的介电常数分布ε(x,y)，利用二值化方法获得水相分布f_w(x,y)，若像素(x,y)为水相，则其值为1，若像素(x,y)为气相或油相，则该值为0；

c.利用超声传感器阵列测量数据结合图像重建算法实现超声传感器被测截面内的声阻抗分布Z(x,y)，利用二值化方法获得气相分布f_g(x,y)，若像素(x,y)为气相，则其值为1，若像素(x,y)为水相或油相，则该值为0；

d.设f_A(x,y)为截面内的相分布，且各像素值均为1，由于油、气、水三相充满管道截面，因此可得油相分布f_o(x,y)＝f_A(x,y)-f_w(x,y)-f_g(x,y)，其中，若像素(x,y)为油相，则其值为1，若像素(x,y)为水相或气相，则该值为0。

参考文献

[1]董峰，魏灿，谭超，基于多截面阻抗式长腰内锥传感器及相关测速的多相流测量方法，发明专利，专利号：ZL201110048191.5，申请日期：2011.02.28，授权日期：2012.12.5

[2]Tan Chao,Xu Yaoyuan,Dong Feng.Determining the boundary of inclusions with known conductivitiesusing a Levenberg-Marquardt algorithm by electrical resistance tomography.Measurement Science andTechnology,2011,22(10):104005

[3]Li Yi,Yang Wuqiang,Xie Cheng-Gang,et al.Gas/oil/water flow measurement by electrical capacitancetomography.Measurement Science and Technology,2013,24(7):74001

[4]Xu L.J.,Xu L.A.Gas/liquid two-phase flow regime identification by ultrasonic tomography.FlowMeasurement and Instrumentation,1998,8(3-4):145～155

Claims (1)

1.一种基于超声多普勒与电学多传感器的多相流可视化测试方法，所采用的测试装置包括电学传感器阵列、超声传感器阵列、电学信号发生与检测单元、超声信号发生与检测单元、流动参数计算与可视化单元；所述的电学传感器阵列包括两组分布在被测管道相同截面位置的电极阵列；所述的超声传感器阵列包括两组分布在被测管道不同截面位置的探头阵列，每组探头阵列由两个以上的超声探头构成，其中第一组以超声多普勒方式工作，包括高频超声探头与低频超声探头，前者用于发射高频超声信号，以获取被测多相流体内的液滴流速，后者用于发射低频的超声信号，以获取被测多相流体内的气泡流速，第二组以超声透射法工作，通过超声幅值衰减获取被测多相流体相含率；所述的电学传感器阵列与超声传感器阵列组成一套双模态传感器阵列；电学信号发生与检测单元选通电学传感器阵列中的一对电极作为激励电极对，其中一个电极为激励电极，另一个电极为接地电极，在所述激励电极对之间建立电学敏感场，测量其余电极间电势差；对于第一组超声探头，超声信号发生与检测单元选通该组探头阵列中的一个探头发出超声波，该探头自身或另一个探头接收超声波，将超声强度转换为电信号并提取频移信息；对于第二组超声探头，超声信号发生与检测单元选通该组探头阵列中的一个探头发出超声波，其余探头中的一个或多个探头接收超声波，并将超声强度转换为电信号并提取幅值信息；电学信号发生与检测单元与超声信号发生与检测单元所获信号一起，送入流动参数计算与可视化单元，实现流型在线识别、相含率与分相流速的计算，以及相分布的可视化重建与显示；该测试方法包括下列步骤：

1).采集电学传感器阵列信号与超声传感器阵列信号；

2).根据电学传感器阵列测量数据和超声传感器阵列测量数据，结合特征提取与分类算法，实现多相流流型的识别；

3).根据超声透射衰减幅度，利用超声传感器测量数据带入计算气相含率α_g；其中U为超声穿透多相流混合物后的幅值，U₀为超声穿透单相水后的幅值，L为超声传播路径长度，X_g为衰减系数，通过实验标定获得；

4).采用多频率超声激励，超声探头测量端获得的超声频率f'与原始发射频率f间存在关系式其中u为超声波在被测介质中的行进速度，u_s为液滴或气泡的移动速度，当液滴或气泡流向测量端时，u_s前方运算符号为-号，反之则为+号；根据将高频超声波反射频率变化带入上式获取液滴流速u_o＝u_s，将低频超声波反射频率变化带入上式获取气泡速度u_g＝u_s，假设油相与水相间不存在速度差，因此水相流速u_w＝u_o；多频超声高频频率范围典型取值为4MHz以上，低频频率范围典型取值为500KHz-1MHz，以上范围根据被测流体密度与流体内声速变化；

5).根据流型识别结果选取相含率估计模型，并判断连续相是否导电，

若连续相为非导电相，采用电容传感器测量数据，计算模型如下：

a.根据电容电极对采集的电压信号，计算多相流体混合介电常数ε_m；

b.已知油相介电常数ε_o、水相介电常数ε_w及气相介电常数ε_g，根据公式 $\sqrt{{\mathrm{\ϵ}}_m}={\mathrm{\α}}_g\sqrt{{\mathrm{\ϵ}}_g}+(1-{\mathrm{\α}}_g)({\mathrm{\α}}_o\sqrt{{\mathrm{\ϵ}}_o}+{\mathrm{\α}}_w\sqrt{{\mathrm{\ϵ}}_w})$ 以及多相流分相含率和为1的关系式α_o+α_g+α_w＝1，计算油相含率α_o、气相含率α_g与水相含率α_w；

若连续相导电，采用电导传感器测量数据，计算模型如下：

a.根据电阻电极对采集的电压信号，计算多相流体混合电导率σ_m；

b.已知油相电导率σ_o、气相电导率σ_g与水相电导率σ_w，另α_l＝α_o+α_w为油气水多相流的液相含率，根据公式和公式α_o+α_g+α_w＝1计算出油相含率α_o、水相含率α_w和气相含率α_g，式中σ_m为多相混合流体电导率，σ_l＝α_oσ_o+α_wσ_w为液相混合电导率；

6).利用所述电学传感器阵列与超声传感器阵列测量数据进行油气水多相流相分布重建：

a.根据流型识别结果判断连续相是否导电，若连续相导电，将电学传感器阵列视为电导传感器阵列，设采用电导传感器测量数据结合图像重建算法实现电学传感器截面内电导率分布σ(x,y)，利用二值化方法获得水相分布f_w(x,y)，其中x与y为像素坐标，若像素(x,y)为水相，则其值为1，若像素(x,y)为气相或油相，则该值为0；

b.若连续相为非导电相，将电学传感器阵列视为电容传感器阵列，采用电容传感器测量数据结合图像重建算法实现电学传感器截面内的介电常数分布ε(x,y)，利用二值化方法获得水相分布f_w(x,y)，若像素(x,y)为水相，则其值为1，若像素(x,y)为气相或油相，则该值为0；

c.利用超声传感器阵列测量数据结合图像重建算法实现超声传感器被测截面内的声阻抗分布Z(x,y)，利用二值化方法获得气相分布f_g(x,y)，若像素(x,y)为气相，则其值为1，若像素(x,y)为水相或油相，则该值为0；

d.设f_A(x,y)为截面内的相分布，且各像素值均为1，由于油、气、水三相充满管道截面，因此可得油相分布f_o(x,y)＝f_A(x,y)-f_w(x,y)-f_g(x,y)，其中，若像素(x,y)为油相，则其值为1，若像素(x,y)为水相或气相，则该值为0。