CN1554936A

CN1554936A - 基于软测量技术的气液两相流流型辨识方法及系统

Info

Publication number: CN1554936A
Application number: CNA2003101227739A
Authority: CN
Inventors: 冀海峰; 黄志尧; 王保良; 李海青
Original assignee: Zhejiang University ZJU
Current assignee: Zhejiang University ZJU
Priority date: 2003-12-19
Filing date: 2003-12-19
Publication date: 2004-12-15
Anticipated expiration: 2023-12-19
Also published as: CN1246683C

Abstract

本发明公开了一种基于软测量技术的气液两相流流型辨识方法及系统。采用了软测量技术框架，利用易测可靠的差压波动信号作为辅助变量，进行主导变量两相流流型的检测。对采集到的差压信号进行小波分析，并在此基础上提取相应的特征值，根据特征值采用模糊模式识别技术对水平管气液两相流流型进行客观辨识。其中，关键的软测量模型采用了小波分析技术与模糊模式识别技术相结合进行建模。在此基础上建立了一套气液两相流流型辨识系统。本系统由实验管路、取压孔、取压管、差压传感器、放大器、数据采集卡、微型计算机以及打印机构成。本系统将检测到的差压信号通过放大器和数据采集卡送入计算机，在计算机中实现对数据处理，进而实现气液两相流流型的辨识。

Description

基于软测量技术的气液两相流流型辨识方法及系统

技术领域

本发明涉及一种基于软测量技术的气液两相流流型辨识方法及系统。

背景技术

随着科学技术的迅速发展，两相流动体系在国民经济和人类生活中的地位日益重要。在自然界和工业生产中涉及范围十分广泛。在化工、石油、冶金、动力及原子能等工业中，两相流动过程更是普遍存在。

两相流动中各相介质的分布状况称为流型，或流态。流型的不同，不但影响两相流的流动特性和传热传质性能，而且影响对两相流参数的准确测量。例如压力的波动与被输送的两相流混合物相互作用会产生管路的强烈振动现象。因此研究两相流流型的自动识别不仅具有重要的工业应用价值和学术价值，而且它能够为相关工业的安全与自动化生产、管路系统的设计与运行、两相流量计量的开发等提供重要的技术支持。因此，两相流的流型辨识是两相流检测中的一个重要方向。

但气液两相流中两相之间存在着多变的相界面，相界面的形状及其在两相流中的分布情况是随着流动过程随时在变化的，两相流的相界面的这种多变性，致使两相流的流型不仅是多种多样的，而且其变化带有随机性。这给两相流的流型检测带来了很大的挑战性。

传统的两相流流型检测(判别)方法有两大类：一类是采用实验方法做出流型图。早期对流型的研究主要是靠流型的状态图和流型的模型判据；另一类是根据对流型转变机理得到的转变关系式，利用现场的流动参数来确定具体的流型。但传统的流型辨识方法在实际应用中有很大的局限性。

流型的现代测量方法从工作原理上可分为直接测量法和间接测量法。常见的直接测量法有目测法、高速摄影法、射线衰减法和接触探头法等几种；间接测量法主要有概率密度函数分析法(PDF)、功率谱密度函数分析法(PSD)。随着科学技术的不断发展，近些年又有许多新的流型测量方法产生：过程层析成像技术、模糊数学判别法、分形几何分析法、WVD分析法等，它们都从信号处理的角度出发对所得到的信息经过不同的方法进行处理，最终得到与流型的某种关系进行流型的分析和辨识。传统的通常建立在传感器等硬件基础上的各种现有测试手段，目前在两相流领域还未能获得令人十分满意的结果，将软测量技术(如状态估计，过程参数辨识，人工神经网络，模式识别等)引入多相流参数测试领域中来是两相流检测的一个发展趋势。通过软测量方法(是利用较易在线测量的辅助过程变量和离线分析信息提供主要过程参数的在线估计的方法)，解决具有复杂性，不确定性，且很难用数学模型精确描述的两相流系统的测试问题。

发明内容

本发明的目的是提供一种基于软测量技术的气液两相流流型辨识方法及系统。

辨识方法：首先进行水平管气液两相流差压信号的检测，对差压信号进行小波分析处理，再利用模糊模式识别技术建立差压信号辅助变量和两相流流型主导变量之间的数学模型，通过此软测量模型实现两相流流型的自动客观辨识。

系统具有两相流水平管上的第一取压孔、第二取压孔、第三取压孔、第四取压孔、并分别与取压管相连、取压管依次与两个差压传感器、信号调制放大模块、数据采集卡、微型计算机以及打印机相连。

本发明在原有的仪表检测的基础上，依托计算机技术，采用软测量技术进行水平管气液两相流流型的辨识。与传统流型辨识系统相比，采用信号准确可靠，辨识结果准确，能够区分水平管气液两相流的主要流型，且程序功能丰富，显示界面友好。本发明利用所提出的方法对水平管气液两相流流型进行了辨识。试验分别在本实验室试验装置和大庆油田多相流试验基地进行，取得了较好的辨识效果，表1示出了部分试验结果。

表1

	气泡流	塞状流	弹状流	环状流
	气泡流	塞状流	弹状流	环状流	实验数据个数	20	16	20	20
判别正确个数	18	13	17	19	实验数据个数	20	16	20	20
判别正确个数	18	13	17	19	判别正确率	95％	81.25％	85％	90％

附图说明

图1是基于软测量技术的气液两相流流型辨识方法方框图；

图2是不同流型下压力波动信号小波分解示意图；

图3是基于软测量技术的气液两相流流型辨识系统结构示意图；

图4是基于软测量技术的气液两相流流型辨识系统取压距离示意图；

图5是本发明信号调制放大模块电路图；

图6是本发明软件程序流程图。

具体实施方式

如图1所示，介绍了基于软测量技术的气液两相流流型辨识方法方框图。软测量技术是近年来在过程控制和检测领域涌现出来的一种新技术。它采用间接测量思路，利用易于获取的其它测量信息，通过计算来实现被检测量的估计。所谓软测量技术就是利用易测过程变量(常称为辅助变量或二次变量)，依据这些易测过程变量与难以直接测量的待测过程变量(常称为主导变量)之间的数学关系(软测量模型)，通过各种数学计算和估计方法，从而实现待测过程变量的测量。可以看出，软测量技术是依据某种最优化准则，利用由辅助变量构成的可测信息，通过软件计算实现对主导变量的测量，软测量的核心是表征辅助变量和主导变量之间的数学关系的软测量模型，因此软测量的本质就是如何建立软测量模型，即一个数学建模问题。软测量系统构造的核心是如何建立软测量模型。软测量的分类一般都是依据软测量模型的建立方法，可以分为：机理建模、回归分析、状态估计、模式识别、人工神经网络、模糊数学、过程层析成像、相关分析和现代非线性信息处理技术等。图1中描述了本发明所采用软测量方案的技术路线。由于两相流系统的复杂性，流型判别的不确定性，使得很多传统的方法都得不到很好的效果。本发明采用小波分析技术和模糊数学进行数学建模，利用软测量技术对不同管径的气液两相流的流型进行了分析和判别。在两相流系统中，压力或差压波动信号的检测方便可靠，而且包含了丰富的信息，因此将气液两相流的差压波动信号作为软测量系统的易测辅助变量。而主导变量就是气液两相流中重要参数——流型。选择各尺度细节的能量占信号的总能量的百分比作为信号各尺度的特征值，即信号经小波分解后，各尺度细节的能量占信号的总能量的百分比(“尺度能量百分比”)。所谓的尺度能量百分比定义如下：

对于正交小波，可以用小波分解系数来表示信号能量，其总能量可以表示为：

E = \underset{j, k}{Σ} {| | D_{j, k} | |}^{2} + \underset{k}{Σ} {| | C_{J_{2, k}} | |}^{2} - - - - J_{1} - 1 \leq j \leq J_{2}

而频率介于2^-j和2^-(j-1)之间的成分的“尺度能量”可以表示为：

E_{j} = \underset{k &Element; Z}{Σ} {| | D_{j, k} | |}^{2} - - - - J_{1} - 1 \leq j \leq J_{2}

本文以各尺度细节的能量占信号的总能量的百分比作为信号各尺度的特征值，即取能量特征值为“尺度能量百分比”

P_{j} = \frac{E_{j}}{E} .

如图2所示，在不同流型下差压波动信号的小波分解的图形是不同的，其各尺度能量占总能量的百分比也有很大区别，可以作为特征值进行流型的判别。

根据模糊数学理论，流型辨识问题就是一个决策问题，即通过已知的信息来进行流型的分类。

首先，设定论域为

U = {\underset{~}{A_{1}}, \underset{~}{A_{2}}, \underset{~}{A_{3}}, \underset{~}{A_{4}}},

其中

即在实验中观测到的实际流型。

然后确定待辨识流型的特征向量X＝{x₁，x₂，x₃，x₄，x₅，x₆，x₇}，这里，x_j＝P_j＝E_j/E₀，1≤j≤J，x₇＝E₀，即选用信号的尺度能量百分比P_j和信号能量E₀作为特征向量。并确定它们的主要分布区间[min_ij，max_ij]，1≤j≤J，，其中i代表某种具体流型，1≤i≤4。也就是说，将根据不同尺度的能量分量P_j和信号的能量E₀的分布状况来确定流型的判别准则，以此作为流型判别的依据。

由于隶属度函数的确定并没有一定的标准，多数是根据经验知识加以确定。经过分析，对每个特征变量选用梯形函数作为隶属度函数：

选取每一特征变量隶属度值的加权均值作为每一的隶属度函数

即：

μ_{\underset{~}{A_{i}}} (X) = \underset{j}{Σ} w_{ij} \underset{~}{μ_{i}} (x_{j}),

其中w_ij为加权系数，且

\underset{j}{Σ} w_{ij} = 1,

加权系数是在实验中确立的。由此可以计算出某一特征向量对于每一模糊子集的隶属度值，按最大隶属度原则就可以判断待辨识流型。

在利用样本数据确立了特征变量的隶属度函数后，整个辨识过程的方法步骤如下：

(1)首先对待辨识信号进行小波分解，提取特征向量值：能量E₀和尺度能量百分比P_j；

(2)然后计算各特征值的隶属度值，从而得到此信号对于每一流型(模糊子集)的隶属度值；

(3)根据最大隶属度原则，得到辨识的结果，如果有两个相同的隶属度值，则认为流型处于两种流型的过渡阶段；

(4)最后，根据实际流型与判别得到的流型进行验证，得到辨识成功率。

本发明还实现了采样频率计算机自动确定。首先进行压力波动信号的预采集，在这一步中，采用较高的采集频率，一般为500Hz的采样频率。将采集到的压力波动信号进行FFT变换(快速Fourier变换)进而求得压力波动信号的频谱值，根据谱图，确定压力波动信号的主要集中的频率范围。得到压力波动信号所处的频率范围后，根据采样定理，将采样频率设置到压力波动信号最高频率的两倍，在实际试验中，对得到的采样频率乘以一个系数作为最后确定的采样频率。这样通过软件本发明实现了采样频率的自动设定。

如图3所示，描述基于软测量技术的气液两相流流型辨识系统结构示意图：采用不同取压距离(第一取压孔(1)，第二取压孔(2)之间和第三取压孔(3)，第四取压孔(4)之间)获取水平管中气液两相流的差压信号，差压信号可通过两种方式进行外传，其一，当差压传感器(6)可输出标准信号时，可直接接入数据采集卡(7)中；其二，若考虑成本所用的差压传感器(6)不能输出标准信号，通过自制的信号调制放大模块进行信号的放大和远传。标准信号或放大后的信号通过工业上已成熟应用的模拟量采集板卡送入计算机进行相应的处理。计算机中可以进行采样频率、采样通道和数据长度等的选择。

如图4所示，描述了基于软测量技术的气液两相流流型辨识系统取压距离示意图。在管道上进行差压信号的获取，采集的差压信号通道多，获取的信息就完整，但由于管路的限制，再则考虑采集的信号通道多，成本提高，但相应的效果并不会明显改善，因此本实用新型中选取不同取压距离的两路差压信号进行测量分析。取压距离的不同，采集到的信号的侧重点是不同的，而且水平管道管径的不同也对取压距离提出一定的要求，综合进行考虑，选取的取压距离分别为L/D＝10和L/D＝2，其中L为选用的取压距离，D为水平管管径。需要说明的一点是这里的选取原则并不是一成不变的，可以根据实际情况进行一定的调整。

如图5所示，描述了本发明信号调制放大模块电路图。包括差压传感器DPS、电压源芯片P1、电压源芯片P2、仪表放大器U1、运算放大器U2、运算放大器U3、三极管T1、第一电阻R1、第二电阻R2、第三电阻R3、第四电阻R4、第五电阻R5、调零电阻R_Z、增益调节电阻R_G、第一电容C1、第二电容C2和第三电容C3。其中，+15V电源接到P1的输入管脚，P1的输出接到DPS的输入管脚1、3。DPS的输出管脚2、4接到U1的输入端。-15V电源接到P2的输入管脚，P2的-5V输出通过R_Z接到P1的+5V输出上，R_Z的可调端接到U2的输入正极，U2的输入负极接到U2的输出端，组成一电压跟随器，用来调节零点。第二电容C2接到P1的+5V电压输出和地，第三电容C3接到P2的-5V输出和地。U2的输出接到U1的参考电源端。U1的放大倍数通过接到U1的R_G进行调节。U1的输出接第一电阻R1和第一电容C1作为信号调制放大模块的电压输出V₀。为进行远传，V。通过第二电阻R2和第一电容C1组成的低通滤波器接到U3的输入正极，U3的负极通过第三电阻R3接到三极管T1的发射极，三极管T1的发射极通过第五电阻R5接到地。U3的输出通过第四电阻R4接到三极管T1基极，+15V电源通过输出负载接到三极管T1的集电极。输出负载上的电流I_OUT即为信号调制放大模块电流输出。仪表放大器U1的输出V₀既此信号调制放大模块的电压输出。电流输出利用运算放大器U3和三极管T1组成的典型电压-电流转换电路来实现，I_OUT为输出电流。

如图6所示，描述了本发明软件程序流程图。程序首先进行硬件的初始化，检查硬件的驱动是否装好，硬件是否都连接，如果成功则进行下面的操作，如果不成功则等待进行检查。若硬件设备初始化成功则进行软件的配置(设置)，可以设置相应的采样频率、采样时间、所选用的通道编号等等。这些按要求配置好后就可进入到程序的主要部分进行数据采集、处理、分析，最后进行水平管气液两相流流型的辨识。

Claims

1.一种基于软测量技术的气液两相流流型辨识方法，其特征在于首先进行水平管气液两相流差压信号的检测，对差压信号进行小波分析处理，再利用模糊模式识别技术建立差压信号辅助变量和两相流流型主导变量之间的数学模型，通过此软测量模型实现两相流流型的自动客观辨识。

2.根据权利要求1所述的一种基于软测量技术的气液两相流流型辨识方法，其特征在于方法的步骤为：

1)数据采集分析：选取的取压距离分别为L/D＝10和L/D＝2，其中L为选用的取压距离，D为水平管管径，对两种不同取压距离的差压信号进行采集，并进行采样频率的自动设定；对差压波动信号进行小波分析，并提取相应的特征值：尺度能量和尺度能量百分比；

2)流型的模糊模式识别：针对水平管气液两相流的气泡流、弹状流、塞状流和环状流主要流型建立论域，并根据小波分析得到的特征值，采用最大隶属度原则进行流型的判别；

3)软测量模型的建立过程：软测量模型的建立是通过对样本数据的训练进行的，寻求不同流型下由小波分析技术得到的特征值的不同特点，确定隶属度函数进而建立软测量模型，即选用信号的尺度能量百分比P_j和信号能量E₀作为特征向量，并确定它们的主要分布区间[min_ij，max_ij]，1≤j≤J，其中i代表某种具体流型，1≤i≤4。根据特征向量的主要分布区间，对每个特征变量选用梯形函数作为隶属度函数，选取每一特征变量隶属度值的加权均值作为每一A_i的隶属度函数

即：

μ_{\underset{&OverBar;}{A_{i}}} (X) = \underset{j}{Σ} w_{ij} \underset{&OverBar;}{μ_{i}} (x_{j}),

其中w_ij为加权系数，且

\underset{j}{Σ} w_{ij} = 1,

最后采用最大隶属度原则进行流型辨识，从而建立了软测量模型。

3.一种基于软测量技术的气液两相流流型辨识系统，其特征在于：它具有两相流水平管上的多个取压孔、并分别与取压管(5)相连、取压管(5)依次与多个差压传感器(6)、信号调制放大模块(7)、数据采集卡(8)、微型计算机(9)以及打印机(10)相连。

4.根据权利要求3所述的一种基于软测量技术的气液两相流流型辨识系统，其特征在于所说多个取压孔为第一取压孔(1)、第二取压孔(2)、第三取压孔(3)、第四取压孔(4)。

5.根据权利要求3所述的一种基于软测量技术的气液两相流流型辨识系统，其特征在于所说信号调制放大模块(7)的电路为：它具有差压传感器DPS、电压源芯片P1、电压源芯片P2、仪表放大器U1、运算放大器U2、运算放大器U3、三极管T1、第一电阻R1、、第二电阻R2、第三电阻R3、第四电阻R4、第五电阻R5、调零电阻R_Z、增益调节电阻R_G、第一电容C1、第二电容C2和第三电容C3。其中，+15V电源接到P1的输入管脚，P1的输出接到DPS的输入管脚1、3。DPS的输出管脚2、4接到U1的输入端。-15V电源接到P2的输入管脚，P2的-5V输出通过R_Z接到P1的+5V输出上，R_Z的可调端接到U2的输入正极，U2的输入负极接到U2的输出端，组成一电压跟随器，用来调节零点。第二电容C2接到P1的+5V电压输出和地，第三电容C3接到P2的-5V输出和地。U2的输出接到U1的参考电源端。U1的放大倍数通过接到U1的R_G进行调节。U1的输出接第一电阻R1和第一电容C1作为信号调制放大模块的电压输出V₀。为进行远传，V_o通过第二电阻R2和第一电容C1组成的低通滤波器接到U3的输入正极，U3的负极通过第三电阻R3接到三极管T1的发射极，三极管T1的发射极通过第五电阻R5接到地。U3的输出通过第四电阻R4接到三极管T1基极，+15V电源通过输出负载接到三极管T1的集电极。输出负载上的电流I_OUT即为信号调制放大模块电流输出。