CN106323589B - 一种密相气力输送两相流流型的归类方法 - Google Patents
一种密相气力输送两相流流型的归类方法 Download PDFInfo
- Publication number
- CN106323589B CN106323589B CN201610659189.4A CN201610659189A CN106323589B CN 106323589 B CN106323589 B CN 106323589B CN 201610659189 A CN201610659189 A CN 201610659189A CN 106323589 B CN106323589 B CN 106323589B
- Authority
- CN
- China
- Prior art keywords
- flow pattern
- classification
- scale
- electrostatic signal
- phase
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Expired - Fee Related
Links
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01M—TESTING STATIC OR DYNAMIC BALANCE OF MACHINES OR STRUCTURES; TESTING OF STRUCTURES OR APPARATUS, NOT OTHERWISE PROVIDED FOR
- G01M10/00—Hydrodynamic testing; Arrangements in or on ship-testing tanks or water tunnels
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01M—TESTING STATIC OR DYNAMIC BALANCE OF MACHINES OR STRUCTURES; TESTING OF STRUCTURES OR APPARATUS, NOT OTHERWISE PROVIDED FOR
- G01M9/00—Aerodynamic testing; Arrangements in or on wind tunnels
- G01M9/06—Measuring arrangements specially adapted for aerodynamic testing
Landscapes
- Physics & Mathematics (AREA)
- Fluid Mechanics (AREA)
- General Physics & Mathematics (AREA)
- Information Retrieval, Db Structures And Fs Structures Therefor (AREA)
Abstract
本发明公开了一种密相气力输送两相流流型的归类方法,获取密相气力输送两相流流型的静电信号;对该静电信号进行多尺度分解;计算静电信号的不同尺度能量比重大小;将不同尺度能量比重及大小关系作为必要参数,对该流型所属的类别进行判断,从而实现密相气力输送两相流流型的归类。该归类方法从气固两相流流型的动力学本质出发,具有较高的科学性及合理性。相对于传统按流动形态划分流型的方法,该方法划分的流型种类减少且彼此间有区别明确的动力学机理。归类方法实用范围广,对密相气力输送两相流中的输送介质的物性没有要求。
Description
技术领域
本发明涉及气固两相流流动检测技术领域,具体地讲,涉及一种密相气力输送两相流流型的归类方法。
背景技术
密相气力输送两相流广泛存在于能源、化工、冶金及医药食品加工等领域。气固两相流的流动状态,即流型,极大的影响着流动参数的准确测量以及流动系统的优化设计和运行稳定性。流型的研究及其在线测量技术的研究对两相流的流动特性,传热传质性能以及其他问题的分析研究,具有重要的科学意义和工业应用价值。
气固两相在流动过程中,气固两相界面及固相的存在状态均发生变化,流动形态多样,流型形成机理复杂。因此,在流型划分方面,不同的研究者给出流型的不同定义和划分。于遵宏等对其进行了总结,目前常用的流型分类方法是根据流动形态,相邻典型流型之间还存在过渡流型。流型种类繁多,典型流型就有7种之多。由于流动条件变化的多样性以及研究角度的多样性,基于流动形态上的流型分类依然有诸多问题,如流型定义模糊、没有定量的划分标准、流型种类较多,相邻流型之间的动力学特征区别不明显等。由此衍生出的问题极大的影响了流型的研究及流型识别技术在气固两相流中的发展和应用,亟待解决。因此,为解决上述问题,需要从气固两相流流动本质出发,对传统流型进行归类,减少流型种类,进而提高流型识别效率,推动流型识别技术在气力输送中的应用。
发明内容
本发明要解决的技术问题是提供一种密相气力输送两相流流型的归类方法,根据两相流波动信号的多尺度能量比重,实现对密相气力输送两相流流型的归类。
本发明采用如下技术方案实现发明目的:
一种密相气力输送两相流流型的归类方法,其特征是:包括如下步骤:
(1)获取气力输送两相流某一流型的静电信号数据;
(2)对获取的每种流型的静电信号进行多尺度分解,具体步骤如下:
首先对静电信号进行经验模态分解处理,得到10个本征模态分量IMF,然后对每个IMF分量进行R/S分形分析得到关系曲线ln[R(τ)/S(τ)]~lnτ,τ表示R/S分形方法中的一个时间延迟变量,最后对关系曲线的直线段部分做线性拟合,得到拟合直线的斜率,即Hurst指数H;
(3)根据分形特征及Hurst指数H的大小,将静电信号划分为尺度1、尺度2和尺度3,具体划分方法如下:
将Hurst指数H小于0.5的IMF分量划入尺度1;将在小的时间延迟τ下Hurst指数H大于0.5小于1而大的时间延迟τ下的Hurst指数H小于0.5的IMF分量划入尺度2;将Hurst指数H大于0.5小于1的IMF划入尺度3;
(4)计算静电信号的尺度1、尺度2及尺度3的能量比重,能量比重的计算方法如下:
静电信号被分解成10个IMF分量,其中尺度1占有M个IMF分量,分别是IMF1,IMF2…,IMFM,一个IMF分量为一组时间序列{x1,x2,…xH},其中,能量EIMF1的计算公式为:
EIMF1=x1 2+x2 2+…+xH 2
其余的能量EIMF2、EIMF3…EIMFM的计算公式与EIMF1的计算公式雷同,利用上述公式分别计算每个IMF分量的能量,则静电信号尺度1的能量比重R1:
按照计算静电信号尺度1的能量比重R1的方法,计算静电信号尺度2能量比重R2以及尺度3能量比重R3;
(5)以步骤(4)计算得到的R1,R2,R3作为必要参数,对该流型所属的流型类别进行判断。
作为对本技术方案的进一步限定,所述尺度1比重R1、尺度2比重R2和尺度3的比重R3之和为1,即:
R1+R2+R3=1。
作为对本技术方案的进一步限定,所述步骤(5)的流型类别共有5种,分别定义为类别1、类别2、类别3、类别4和类别5,类别1、类别2、类别3、类别4和类别5分别代表微尺度流型、微介均衡流型、介尺度流型、介宏均衡流型和宏尺度流型,具体判断方法如下:
当满足以下关系式时,可归为“类别1”,
当满足以下关系式时,可归为“类别2”,
当满足以下关系式时,可归为“类别3”,
当满足以下关系式时,可归为“类别4”,
当满足以下关系式时,可归为“类别5”,
现有技术相比,本发明的优点和积极效果是:该归类方法从气固两相流流型的动力学本质出发,具有较高的科学性及合理性。相对于传统按流动形态划分流型的方法,该方法划分的流型种类减少且彼此间有区别明确的动力学机理。归类方法实用范围广,对气力输送两相流中的输送介质的物性没有要求。
附图说明
图1是密相气力输送两相流的一组静电信号。
图2是静电信号的10个IMF分量。
图3是静电信号的10个IMF分量的关系曲线ln[R(τ)/S(τ)]~lnτ。
图4是5种流型类别的分布图。
具体实施方式
下面结合附图,对本发明的一个具体实施方式进行详细描述,但应当理解本发明的保护范围并不受具体实施方式的限制。
如图1-图4所示,1、获取密相气力输送两相流流型的静电信号数据,图1所示为一组静电信号数据。
2、静电信号进行多尺度分解,具体步骤如下:
首先对静电信号进行经验模态分解处理(此为现有技术,在此不再赘述),得到10个本征模态分量IMFs,然后对每个IMF分量进行R/S分形分析得到关系曲线ln[R(τ)/S(τ)]~lnτ,最后对关系曲线的直线段部分做线性拟合,得到拟合直线的斜率,即Hurst指数H;图2所示为一组静电信号的10个IMF分量,每个分量的关系曲线ln[R(τ)/S(τ)]~lnτ,如图3所示,其中,对每个IMF分量进行R/S分形处理时所选的时间延迟τ的范围为10~100。
根据分形特征及Hurst指数H的大小,将静电信号划分为尺度1、尺度2和尺度3,具体划分方法如下:
将Hurst指数H小于0.5的IMF分量划入尺度1;将在小的时间延迟τ下Hurst指数H大于0.5小于1而大的时间延迟τ下的Hurst指数H小于0.5的IMF分量划入尺度2;将Hurst指数H恒大于0.5小于1的IMF划入尺度3;
4、计算静电信号的尺度1、尺度2及尺度3的能量比重,能量比重的计算方法如下:
静电信号被分解成10个IMF分量,其中尺度1占有M个IMF分量,分别是IMF1,IMF2…,IMFM。一个IMF分量为一组时间序列{x1,x2,…xH},其能量EIMF1的计算公式为:
EIMF1=x1 2+x2 2+…+xH 2
其余的能量EIMF2、EIMF3…EIMFM的计算公式与EIMF1的计算公式雷同,IMF分量的表述方式也雷同,利用上述公式分别计算每个IMF分量的能量,则静电信号尺度1的能量比重R1:
按照计算静电信号尺度1的能量比重R1的方法,计算静电信号尺度2能量比重R2以及尺度3能量比重R3。并且,R1+R2+R3=1
5、根据步骤4计算得到的静电信号3个尺度的能量比重,对该流型所属的流型类别进行判断。流型类别共有5种,分别定义为“类别1”,“类别2”,“类别3”,“类别4”,和“类别5”,类别1、类别2、类别3、类别4和类别5分别代表微尺度流型、微介均衡流型、介尺度流型、介宏均衡流型和宏尺度流型,具体判断方法如下:
当满足以下关系式时,可归为“类别1”,
当满足以下关系式时,可归为“类别2”,
当满足以下关系式时,可归为“类别3”,
当满足以下关系式时,可归为“类别4”,
当满足以下关系式时,可归为“类别5”,
分别以R1和R2为横纵坐标,绘制5种流型类别的分布图,如图4所示。
以上公开的仅为本发明的一个具体实施例,但是,本发明并非局限于此,任何本领域的技术人员能思之的变化都应落入本发明的保护范围。
Claims (3)
1.一种密相气力输送两相流流型的归类方法,其特征是:包括如下步骤:
(1)获取气力输送两相流某一流型的静电信号数据;
(2)对获取的每种流型的静电信号进行多尺度分解,具体步骤如下:
首先对静电信号进行经验模态分解处理,得到10个本征模态分量IMF,然后对每个IMF分量进行R/S分形分析得到关系曲线ln[R(τ)/S(τ)]~lnτ,最后对关系曲线的直线段部分做线性拟合,得到拟合直线的斜率,即Hurst指数H;
(3)根据分形特征及Hurst指数H的大小,将静电信号划分为尺度1、尺度2和尺度3,具体划分方法如下:
将Hurst指数H小于0.5的IMF分量划入尺度1;将在小的时间延迟τ下Hurst指数H大于0.5小于1而大的时间延迟τ下的Hurst指数H小于0.5的IMF分量划入尺度2;将Hurst指数H大于0.5小于1的IMF划入尺度3;
(4)计算静电信号的尺度1、尺度2及尺度3的能量比重,能量比重的计算方法如下:
静电信号被分解成10个IMF分量,其中尺度1占有M个IMF分量,分别是IMF1,IMF2…,IMFM,一个IMF分量为一组时间序列{x1,x2,…xH},x1,x2,…xH是按时间顺序排列的一组数据,xi表示一个数据,H是数据的个数或者是时间序列的长度,其中,能量EIMF1的计算公式为:
EIMF1=x1 2+x2 2+…+xH 2
其余的能量EIMF2、EIMF3…EIMFM的计算公式与EIMF1的计算公式雷同,利用上述公式分别计算每个IMF分量的能量,则静电信号尺度1的能量比重R1:
按照计算静电信号尺度1的能量比重R1的方法,计算静电信号尺度2能量比重R2以及尺度3能量比重R3;
(5)以步骤(4)计算得到的R1,R2,R3作为必要参数,对该流型所属的流型类别进行判断。
2.根据权利要求1所述的密相气力输送两相流流型的归类方法,其特征是:所述尺度1比重R1、尺度2比重R2和尺度3的比重R3之和为1,即:R1+R2+R3=1。
3.根据权利要求1所述的密相气力输送两相流流型的归类方法,其特征是:所述步骤(5)的流型类别共有5种,分别定义为类别1、类别2、类别3、类别4和类别5,类别1、类别2、类别3、类别4和类别5分别代表微尺度流型、微介均衡流型、介尺度流型、介宏均衡流型和宏尺度流型,具体判断方法如下:
当满足以下关系式时,可归为“类别1”,
当满足以下关系式时,可归为“类别2”,
当满足以下关系式时,可归为“类别3”,
当满足以下关系式时,可归为“类别4”,
当满足以下关系式时,可归为“类别5”,
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610659189.4A CN106323589B (zh) | 2016-08-11 | 2016-08-11 | 一种密相气力输送两相流流型的归类方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610659189.4A CN106323589B (zh) | 2016-08-11 | 2016-08-11 | 一种密相气力输送两相流流型的归类方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN106323589A CN106323589A (zh) | 2017-01-11 |
CN106323589B true CN106323589B (zh) | 2018-07-06 |
Family
ID=57740739
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201610659189.4A Expired - Fee Related CN106323589B (zh) | 2016-08-11 | 2016-08-11 | 一种密相气力输送两相流流型的归类方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106323589B (zh) |
Families Citing this family (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107764889B (zh) * | 2017-09-27 | 2020-04-21 | 西安理工大学 | 一种基于信号能量法的气固两相流粉尘含碳量测量方法 |
CN108573236B (zh) * | 2018-04-22 | 2021-06-25 | 西安电子科技大学 | 基于离散分数布朗随机场的云天背景下红外弱小目标检测方法 |
CN117519067A (zh) * | 2023-10-20 | 2024-02-06 | 东北大学 | 一种连轧过程多机架控制性能评价方法 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1987485A (zh) * | 2006-11-29 | 2007-06-27 | 东南大学 | 气固两相管流颗粒速度的静电感应空间滤波测量方法 |
CN101900743A (zh) * | 2010-05-28 | 2010-12-01 | 东南大学 | 颗粒速度的线性静电传感器阵列测量方法及装置 |
CN103499516A (zh) * | 2013-10-22 | 2014-01-08 | 东南大学 | 一种高压密相气力输送煤粉流动状态的检测方法及检测装置 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP2399115B1 (en) * | 2009-02-18 | 2016-01-06 | Battelle Memorial Institute | Small area electrostatic aerosol collector |
-
2016
- 2016-08-11 CN CN201610659189.4A patent/CN106323589B/zh not_active Expired - Fee Related
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1987485A (zh) * | 2006-11-29 | 2007-06-27 | 东南大学 | 气固两相管流颗粒速度的静电感应空间滤波测量方法 |
CN101900743A (zh) * | 2010-05-28 | 2010-12-01 | 东南大学 | 颗粒速度的线性静电传感器阵列测量方法及装置 |
CN103499516A (zh) * | 2013-10-22 | 2014-01-08 | 东南大学 | 一种高压密相气力输送煤粉流动状态的检测方法及检测装置 |
Non-Patent Citations (4)
Title |
---|
密相气力输送中气固两相流动特性多源信息分析;付飞飞 等;《化工学报》;20121031;第63卷(第10期);3070-3079 * |
密相气力输送颗粒静电波动信号多尺度分析;许传龙 等;《工程热物理学报》;20110228;第32卷(第2期);239-242 * |
气固两相流颗粒荷电及流动参数检测方法研究;许传龙;《中国优秀博硕士学位论文全文数据库 (博士) 工程科技Ⅰ辑》;20070415(第04期);B015-2 * |
稠密气固两相流静电与压力信号多尺度分析;付飞飞 等;《中国电机工程学报》;20120915;第32卷(第26期);72-78 * |
Also Published As
Publication number | Publication date |
---|---|
CN106323589A (zh) | 2017-01-11 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106323589B (zh) | 一种密相气力输送两相流流型的归类方法 | |
CN106845717B (zh) | 一种基于多模型融合策略的能源效率评价方法 | |
CN104317838B (zh) | 一种基于耦合鉴别性字典的跨媒体哈希索引方法 | |
CN106121622B (zh) | 一种基于示功图的有杆泵抽油井的多故障诊断方法 | |
CN103063345B (zh) | 一种检测轴销所承受径向力的方法 | |
CN100412869C (zh) | 一种改进的基于文档结构的文档相似性度量方法 | |
CN104331591A (zh) | 基于支持向量回归的粮仓储粮数量检测方法 | |
CN109545290B (zh) | 一种基于Voronoi分形技术的非晶合金自由体积检测方法 | |
Brown et al. | ON THE OFFSET OF BARRED GALAXIES FROM THE BLACK HOLE MBH–σ RELATIONSHIP | |
Wang et al. | Cross-correlation focus method with an electrostatic sensor array for local particle velocity measurement in dilute gas–solid two-phase flow | |
Liu et al. | Estimating inter-regional trade flows in China: A sector-specific statistical model | |
CN104216993A (zh) | 一种标签共现的标签聚类方法 | |
CN105574265B (zh) | 面向模型检索的装配体模型定量描述方法 | |
CN114510974B (zh) | 一种多孔介质内气液两相流流型智能识别方法 | |
CN105447521A (zh) | 一种K-means聚类的初值选择方法 | |
CN102147383A (zh) | 多截面阻抗式长腰内锥传感器及多相流测量装置 | |
CN106407601A (zh) | 一种基于数据挖掘技术的气动特性数据处理方法 | |
CN109614074A (zh) | 基于概率转移矩阵模型的近似加法器可靠度计算方法 | |
CN104778363A (zh) | 基于多元时间序列的河流混沌特性识别方法 | |
Han et al. | Numerical methods for analyzing the aerodynamic characteristics of cross parachute with permeability | |
CN106250917A (zh) | 一种基于加速近端梯度pca的时序评分异常值检测方法 | |
CN106295569A (zh) | 一种密相气力输送两相流流型的识别方法 | |
Zhang et al. | Numerical and experimental research on pentagonal cross-section of the averaging Pitot tube | |
CN104897331B (zh) | 面向家电产品的宽量程空气性能检测系统及其方法 | |
Gu et al. | Application of the new pattern recognition system in the new e-nose to detecting Chinese spirits |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20180706 Termination date: 20190811 |