多相流研究背景.ppt

气液两相流小波能量熵分析 宋西姊 1014203040 吴 昊 1014203046 多相流现象广泛存在于自然界和生产生活中。其中，多相管流在石油工业、化学工业、核工业、航空航天、能源工业等生产过程中都得到了广泛的应用。 多相流研究背景 多相管流中，流型的准确预测不但对多相流的工业设备设计和正常运行具有重要意义，通过结合测量模型还能够提高其他流动参数的测量精度。 水平管气液两相流流型 气液两相流流型 多相流参数检测目前的研究方向 传统的单相仪表与流动模型结合的方法： 差压法、容积法、科里奥利力流量计等。 利用现代信号处理技术的软测量方法： 模糊理论、神经网络 、小波分析、复杂网 络等。 应用新兴的两相流测量技术： 示踪粒子检测、太赫兹射线技术等。 两相流研究领域中，除流型外的主要流动参数包括：相含率、流量、速度、滑速比、密度和压力降等。 多相流参数检测 环形电极阵列，简称电导环 以非侵入的方式获得被测流体电导率的空间平均值来反映介质分布变化。 方波电流激励： 2mAP-P、20kHz，电压测量输出频率1kHz（采样频率200kHz）。 特点：平均含率、平均流量；连续相须为导电相 电导环介绍 Mallat 小波分解 任意信号经过式Mallat小波分解后，各尺度单支系数（重构信号）被分布在不同频带上，Mallat小波分解的示意图如图所示， Mallat小波分解示意图 Mallat 小波分解 若信号的采样频率为Fs，经 j 尺度分解后，信号分布在各频带上的信号长度以及频率范围 与各子频带的频率范围的关系如表2-1所示。 信号在各层上的信号长度及子频带的频率分布范围 信号 数据长度 频率范围 S(0) N [0，Fs/21] S(1) N/2 [Fs/22，Fs/21] S(2) N/4 [Fs/23，Fs/22] … … … S(j) N/2j [0，Fs/2j] 小波相对能量 选用正交db4小波对长度为 的两相流的电导波动信号在 空间上进行层正交分解， ，得到的小波系数为 ， ，因此被分解到每层上的信号能量 为： 根据能量守恒定理，经小波变：换后两相流的电导信号的总能量可表示为： 定义相对能量 为： 因此，可以把 视作系统在时频域上的能量概率分布，而且 。 Shannon小波能量熵 将Mallat小波变换与Shannon熵相结合，研究基于小波相对能量的Shannon小波相对能量熵，为了方便表达，在此将Shannon小波相对能量熵称为Shannon小波能量熵。 为了更好理解两相流在时频域上的时频谱分布情况及频率结构随时间变化的情况，在两相流电导波动信号的小波系数 上定义滑动数据窗，设窗宽的取值范围为 ，滑动步长的取值范围 。则滑动数据窗可表示为： 式中， ，其中 ，将滑动窗口划分为 个区间，且有： 式中， 且 。 Shannon小波能量熵 设 表示小波系数 对应的能量落在区间 中的概率，则在第 尺度下，时间段 的Shannon小波能量熵（Shannon Wavelet Energy Entropy，SWEE）可表示为： 经过数据窗在小波分解系数上的滑动，计算每个窗的小波能量熵，作为衡量该窗口数据的复杂程度，滑动步长为小波能量熵的时间分辨率，结合信号的小波分解系数所对应的频率范围，这样就能从时频多尺度上考察两相流动的复杂性。 Shannon小波能量熵的仿真分析 将Shannon小波熵用于由两个含有单频成分的正弦信号与白噪声共同激励下的非线性系统的多尺度复杂性，下式为相应的仿真信号表达式。在仿真实验中，采样频率均为1kHz，两个激励正弦信号的采集时间均为10s。 式中， 为白噪声，其长度与仿真信号中采集10s的正弦信号长度保持一致。在