CN111189915B - 一种水力机械空化发生的实时判定方法 - Google Patents
一种水力机械空化发生的实时判定方法 Download PDFInfo
- Publication number
- CN111189915B CN111189915B CN202010032574.2A CN202010032574A CN111189915B CN 111189915 B CN111189915 B CN 111189915B CN 202010032574 A CN202010032574 A CN 202010032574A CN 111189915 B CN111189915 B CN 111189915B
- Authority
- CN
- China
- Prior art keywords
- frequency domain
- cavitation
- underwater acoustic
- standard deviation
- hydraulic
- 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.)
- Active
Links
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N29/00—Investigating or analysing materials by the use of ultrasonic, sonic or infrasonic waves; Visualisation of the interior of objects by transmitting ultrasonic or sonic waves through the object
- G01N29/04—Analysing solids
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N29/00—Investigating or analysing materials by the use of ultrasonic, sonic or infrasonic waves; Visualisation of the interior of objects by transmitting ultrasonic or sonic waves through the object
- G01N29/44—Processing the detected response signal, e.g. electronic circuits specially adapted therefor
- G01N29/449—Statistical methods not provided for in G01N29/4409, e.g. averaging, smoothing and interpolation
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N29/00—Investigating or analysing materials by the use of ultrasonic, sonic or infrasonic waves; Visualisation of the interior of objects by transmitting ultrasonic or sonic waves through the object
- G01N29/44—Processing the detected response signal, e.g. electronic circuits specially adapted therefor
- G01N29/46—Processing the detected response signal, e.g. electronic circuits specially adapted therefor by spectral analysis, e.g. Fourier analysis or wavelet analysis
Abstract
本发明提出一种水力机械空化发生的实时判定方法,所述方法采用小波分析的方法对水力机械水声信号进行小波分析,获得基于时间序列的水声信号频域能量分布情况并据此对水力机械是否发生空化进行判定。在本发明中采用标准差对水力机械水声信号频域内能量分布离散程度和整个时间序列的各时刻频域主要能量的分布范围的离散程度进行评价。最终依靠比较该标准差与水力机械空化发生门槛值大小的方法来判定水力机械是否发生空化。所述方法准确率高,能够有效的判定出当前时刻水力机械是否发生了空化现象。
Description
技术领域
本发明属于水力机械运行技术领域,特别是涉及一种水力机械空化发生的实时判定方法。
背景技术
随着用户对水力机械稳定性要求的逐步提高,保证水力机械特别是巨型水力机械在无空化状态下安全稳定运行已成为考核机组运行状态的一项重要指标。
空化是液体内局部压力降低时,液体内部或液固交界面上蒸气或气体的空穴(空泡)的形成、发展和溃灭的过程,是液体与气体的相变过程,包括空泡的发生、发展和溃灭三个阶段。空泡发生阶段会伴有噪声辐射的大幅度上升;空泡溃灭阶段能够释放出巨大的压能和热能,是造成液体与固体交界面材料破坏的主要过程;在空泡初生与溃灭之间的是空化发展阶段,它主要是干扰空化物体的运动、使作用力发生变化,譬如推力下降、部件振动摆度增大等。
鉴于在目前的技术水平,仅仅当空化现象严重到足以影响水力机械外部特性,即振动摆度增大、噪声增大时,才能够通过外部设备信号信息再依靠技术人员的经验按照概率估计出可能是发生了空化现象。这种方法的准确率不高,而且在空化现象的初始阶段现有的所有检测设备均无法检测。因此,有必要开发一种能够实时准确判定水力机械是否发生空化的方法。
发明内容
本发明目的是为了解决现有技术中的问题,提出了一种水力机械空化发生的实时判定方法。
本发明是通过以下技术方案实现的,本发明提出一种水力机械空化发生的实时判定方法,所述方法包括以下步骤:
步骤1、启动计算机系统,使得水力机械处于非过渡工况稳定运行状态;
步骤2、采集水力机械水声信号;
步骤3、对水力机械水声信号进行小波分析,从而获得基于时间序列ti={t1,t2,...,tm}的水声信号频域能量分布Ef:
式中:ti表示第i时刻,i=1,2,...,m;m为自然数;Ef为水声信号频域能量分布;Fi表示第i时刻水声信号对应不同频率的能量分布,i=1,2,...,mfij表示水声信号第i时刻第j点频率,j=1,2,...,n;n为自然数;vij表示对应第i时刻第j点频率的水声信号频域幅值;
步骤4、对第i时刻水声信号频域幅值系列vij={vi1,vi2,…,vin}进行相对值计算:
式中:vijN表示频域幅值vij的相对值;vimax表示第i时刻水声信号频域幅值系列中的最大值;
步骤5、计算第i时刻频域幅值的相对值系列vijN={vi1N,vi2N,…,vinN}的标准差:
步骤6、计算基于时间序列ti={t1,t2,…,tm}的各时刻频域幅值的相对值的标准差SiN结果的标准差:
式中:表示基于时间序列ti={t1,t2,…,tm}的各时刻频域幅值的相对值的标准差SiN结果的标准差;表示基于时间序列ti={t1,t2,...,tm}的各时刻频域幅值的相对值的标准差SiN结果的平均值;
步骤7、利用空化发生指征系数δ判定水力机械空化的发生,所述空化发生指征系数δ按下式计算:
当空化发生指征系数δ大于零,即δ>0时,水力机械发生空化;
当空化发生指征系数δ不大于零,即δ≤0时,水力机械未发生空化。
本发明提出的一种水力机械空化发生的实时判定方法采用小波分析的方法对水力机械水声信号进行小波分析,获得基于时间序列的水声信号频域能量分布情况并据此对水力机械是否发生空化进行判定。在本发明中采用标准差对水力机械水声信号频域内能量分布离散程度和整个时间序列的各时刻频域主要能量的分布范围的离散程度进行评价。所述标准差就是样本平均数方差的开平方,是表征样本数据离散程度的一种数学方法。当样本数据越聚集时,计算出的标准差就越小;当样本数据越离散,则计算出的标准差就越大。鉴于不同安装位置也就是距离空化发生区域不同的情况下水力机械水声信号相应频率对应的频域幅值(能量)是不同的:距离空化发生区域越近,则水力机械水声信号相应频率对应的频域幅值越大;距离空化发生区域越远,则水力机械水声信号相应频率对应的频域幅值越小。为保证能够采用统一的水力机械空化发生门槛值判定水力机械空化发生与否,采用按照水声信号频域幅值相对值,即该时刻水声信号频域幅值的最大值分别除以该时刻水声信号频域幅值系列中的各个频域幅值,用标准差来计算水力机械水声信号频域内能量分布离散程度和整个时间序列的各时刻频域主要能量的分布范围的离散程度进行评价。根据水力机械在发生空化后水声信号在整个频域内能量分布离散水力机械水声信号频域内能量分布离散程度和整个时间序列的各时刻频域主要能量的分布范围的离散程度进行评价均急剧增大的特点,通过计算标准差的方法评价采用小波分析获得的水力机械水声信号在不同时刻下整个频域内的能量分布离散程度和整个时间序列的各时刻频域主要能量的分布范围的离散程度变化的总体趋势,最终依靠比较以此得到的标准差与水力机械空化发生门槛值的大小,即空化发生指征系数δ是否大于零来判定水力机械是否发生空化:当空化发生指征系数δ大于零,即δ>0时,则表明水力机械发生空化;当空化发生指征系数δ不大于零,即δ≤0时,则说明水力机械未发生空化。
附图说明
图1为水力机械发生空化前后水声信号频域的能量分布示意图;
图2为未发生空化时的水声信号时频特性示意图;
图3为发生空化时的水声信号时频特性示意图。
具体实施方式
下面将结合本发明实施例中的附图对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
本发明提出一种水力机械空化发生的实时判定方法,所述方法包括以下步骤:
步骤1、启动计算机系统,使得水力机械处于非过渡工况稳定运行状态;
步骤2、采集水力机械水声信号;
步骤3、对水力机械水声信号进行小波分析,从而获得基于时间序列ti={t1,t2,…,tm}的水声信号频域能量分布Ef:
式中:ti表示第i时刻,i=1,2,...,m;m为自然数;Ef为水声信号频域能量分布;Fi表示第i时刻水声信号对应不同频率的能量分布,i=1,2,...,mfij表示水声信号第i时刻第j点频率,j=1,2,...,n;n为自然数;vij表示对应第i时刻第j点频率的水声信号频域幅值;
步骤4、对第i时刻水声信号频域幅值系列vij={vi1,vi2,...,vin}进行相对值计算:
式中:vijN表示频域幅值vij的相对值;vimax表示第i时刻水声信号频域幅值系列中的最大值;
步骤5、计算第i时刻频域幅值的相对值系列vijN={vi1N,vi2N,...,vinN}的标准差:
步骤6、计算基于时间序列ti={t1,t2,...,tm}的各时刻频域幅值的相对值的标准差SiN结果的标准差:
式中:表示基于时间序列ti={t1,t2,...,tm}的各时刻频域幅值的相对值的标准差SiN结果的标准差;表示基于时间序列ti={t1,t2,...,tm}的各时刻频域幅值的相对值的标准差SiN结果的平均值;
步骤7、利用空化发生指征系数δ判定水力机械空化的发生,所述空化发生指征系数δ按下式计算:
当空化发生指征系数δ大于零,即δ>0时,水力机械发生空化;
当空化发生指征系数δ不大于零,即δ≤0时,水力机械未发生空化。
就某一时刻而言,相对于未发生空化的状态,水力机械在发生空化现象后,水声信号的特性发生很大的变化,如图1所示,在整个频域范围内,未发生空化的水力机械水声信号能量分布系列1的能量比较低且分布比较均匀、变化较小,相应的水声信号在整个频域内的能量分布离散度较小;而空化发生后的水力机械水声信号能量分布系列2的能量比较高且分布发生剧烈变化,各频段能量均呈现增大的趋势且频率越低能量增加的绝对值越大,水声信号在整个频域内的能量分布离散度急剧增大。
比较图2和图3所示发生空化前后的水力机械水声信号的时频特性可以发现,未发生空化时,由于水力机械水声信号能量分布系列1的能量比较低且分布比较均匀、变化较小,故而水力机械水声信号在整个时间序列的各时刻频域主要能量的分布范围比较小,且各时刻频域主要能量的分布范围变化也比较小;而一旦发生空化,由于水力机械水声信号能量分布系列2的能量比较高且分布发生剧烈变化,则水力机械水声信号在整个时间序列的各时刻频域主要能量的分布范围大幅度增大,且各时刻频域主要能量的分布范围也随之发生剧烈地变化。
综上所述,通过分析发生空化前后的水力机械水声信号在整个时间序列的各时刻频域主要能量的分布范围的离散程度即可区分出空化现象发生与否,即水力机械水声信号在整个时间序列的各时刻频域主要能量的分布范围的离散程度小,则水力机械未发生空化;反之,空化现象就发生了。
以上对本发明所提出的一种水力机械空化发生的实时判定方法,进行了详细介绍,本文中应用了具体个例对本发明的原理及实施方式进行了阐述,以上实施例的说明只是用于帮助理解本发明的方法及其核心思想;同时,对于本领域的一般技术人员,依据本发明的思想,在具体实施方式及应用范围上均会有改变之处,综上所述,本说明书内容不应理解为对本发明的限制。
Claims (1)
1.一种水力机械空化发生的实时判定方法,其特征在于:所述方法包括以下步骤:
步骤1、启动计算机系统,使得水力机械处于非过渡工况稳定运行状态;
步骤2、采集水力机械水声信号;
步骤3、对水力机械水声信号进行小波分析,从而获得基于时间序列ti={t1,t2,...,tm}的水声信号频域能量分布Ef:
式中:ti表示第i时刻,i=1,2,...,m;m为自然数;Ef为水声信号频域能量分布;Fi表示第i时刻水声信号对应不同频率的能量分布,i=1,2,...,mfij表示水声信号第i时刻第j点频率,j=1,2,...,n;n为自然数;vij表示对应第i时刻第j点频率的水声信号频域幅值;
步骤4、对第i时刻水声信号频域幅值序列vij={vi1,vi2,...,vin}进行相对值计算:
式中:vijN表示频域幅值vij的相对值;vimax表示第i时刻水声信号频域幅值序列中的最大值;
步骤5、计算第i时刻频域幅值的相对值序列vijN={vi1N,vi2N,...,vinN}的标准差:
步骤6、计算基于时间序列ti={t1,t2,...,tm}的各时刻频域幅值的相对值的标准差SiN结果的标准差:
式中:表示基于时间序列ti={t1,t2,...,tm}的各时刻频域幅值的相对值的标准差SiN结果的标准差;表示基于时间序列ti={t1,t2,...,tm}的各时刻频域幅值的相对值的标准差SiN结果的平均值;
步骤7、利用空化发生指征系数δ判定水力机械空化的发生,所述空化发生指征系数δ按下式计算:
当空化发生指征系数δ大于零,即δ>0时,水力机械发生空化;
当空化发生指征系数δ不大于零,即δ≤0时,水力机械未发生空化。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010032574.2A CN111189915B (zh) | 2020-01-13 | 2020-01-13 | 一种水力机械空化发生的实时判定方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010032574.2A CN111189915B (zh) | 2020-01-13 | 2020-01-13 | 一种水力机械空化发生的实时判定方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN111189915A CN111189915A (zh) | 2020-05-22 |
CN111189915B true CN111189915B (zh) | 2022-08-19 |
Family
ID=70706228
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202010032574.2A Active CN111189915B (zh) | 2020-01-13 | 2020-01-13 | 一种水力机械空化发生的实时判定方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN111189915B (zh) |
Families Citing this family (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111695212B (zh) * | 2020-05-25 | 2022-09-09 | 武汉大学 | 精确判别水力机械空化区域的方法 |
CN111706455A (zh) * | 2020-05-28 | 2020-09-25 | 哈尔滨泽级源码大数据科技有限公司 | 一种水力机械转轮叶片进口脱流的判定方法 |
Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPH09146586A (ja) * | 1995-11-20 | 1997-06-06 | Sony Corp | 音声認識装置 |
JP2007303889A (ja) * | 2006-05-09 | 2007-11-22 | Ho Jinyama | キャビテーション検知方法、評価方法、コンピュータプログラム、及び、キャビテーション検知装置 |
CN102546061A (zh) * | 2012-01-09 | 2012-07-04 | 西安交通大学 | 基于小波变换的自适应时频空穴检测方法 |
CN103162815A (zh) * | 2012-12-11 | 2013-06-19 | 哈尔滨电机厂有限责任公司 | 确定水力机械振动信号功率谱高低能量区域的方法 |
CN103592152A (zh) * | 2013-11-20 | 2014-02-19 | 哈尔滨电机厂有限责任公司 | 确定模型水轮机转轮叶片进口部位空化的声学方法 |
CN104535275A (zh) * | 2014-12-11 | 2015-04-22 | 天津大学 | 基于气泡声学的水下气体泄漏量的检测方法和检测装置 |
EP3039329A1 (en) * | 2013-08-27 | 2016-07-06 | Fisher Controls International Llc | Method of cavitation/flashing detection in or near process control valve |
DE102015201203A1 (de) * | 2015-01-26 | 2016-07-28 | Voith Patent Gmbh | Verfahren zur Detektion von Kavitation |
Family Cites Families (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
GB0916635D0 (en) * | 2009-09-22 | 2009-11-04 | Isis Innovation | Ultrasound systems |
CL2014003603A1 (es) * | 2014-12-30 | 2015-07-10 | Univ Chile | Dispositivo para clasificación de tamaño de burbujas en un medio líquido, que comprende un transductor eléctrico emisor de señales ultrasónicas, un transductor eléctrico receptor de señales ultrasónicas, circuitos transmisores y receptores de señales ultrasónicas, un convertidor analógico-digital, y un procesador para señales digitalizadas; método asociado |
-
2020
- 2020-01-13 CN CN202010032574.2A patent/CN111189915B/zh active Active
Patent Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPH09146586A (ja) * | 1995-11-20 | 1997-06-06 | Sony Corp | 音声認識装置 |
JP2007303889A (ja) * | 2006-05-09 | 2007-11-22 | Ho Jinyama | キャビテーション検知方法、評価方法、コンピュータプログラム、及び、キャビテーション検知装置 |
CN102546061A (zh) * | 2012-01-09 | 2012-07-04 | 西安交通大学 | 基于小波变换的自适应时频空穴检测方法 |
CN103162815A (zh) * | 2012-12-11 | 2013-06-19 | 哈尔滨电机厂有限责任公司 | 确定水力机械振动信号功率谱高低能量区域的方法 |
EP3039329A1 (en) * | 2013-08-27 | 2016-07-06 | Fisher Controls International Llc | Method of cavitation/flashing detection in or near process control valve |
CN103592152A (zh) * | 2013-11-20 | 2014-02-19 | 哈尔滨电机厂有限责任公司 | 确定模型水轮机转轮叶片进口部位空化的声学方法 |
CN104535275A (zh) * | 2014-12-11 | 2015-04-22 | 天津大学 | 基于气泡声学的水下气体泄漏量的检测方法和检测装置 |
DE102015201203A1 (de) * | 2015-01-26 | 2016-07-28 | Voith Patent Gmbh | Verfahren zur Detektion von Kavitation |
Non-Patent Citations (5)
Title |
---|
Statistical Analysis of Acoustic Signal for Cavitation Detection;Salem. A. Al-Hashmi;《International Journal of Emerging Technology and Advanced Engineering》;20131031(第3期);全文 * |
基于小波包频带分析技术的水轮机空化特性研究;史会轩等;《大电机技术》;20090115(第01期);全文 * |
基于深度学习的离心泵空化状态识别;曹玉良等;《西安交通大学学报》;20171031;全文 * |
离心泵空化故障识别的时域特征分析方法研究;苏永生等;《水泵技术》;20100831(第04期);全文 * |
离心泵空化现象的模拟试验研究;夏远志等;《机电工程》;20151130(第11期);全文 * |
Also Published As
Publication number | Publication date |
---|---|
CN111189915A (zh) | 2020-05-22 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN111189915B (zh) | 一种水力机械空化发生的实时判定方法 | |
CN102842305B (zh) | 一种基音检测的方法和装置 | |
CN107908863A (zh) | 一种基于emd理论与hht变换的水轮机运转状态判定方法 | |
CN102982347B (zh) | 一种基于kl距离的电能质量扰动分类方法 | |
CN103235953B (zh) | 一种光纤分布式扰动传感器模式识别的方法 | |
JP5760899B2 (ja) | 電磁干渉検出装置、電磁干渉検出方法および電磁干渉検出プログラム | |
Cencîc et al. | Study of erosive cavitation detection in pump mode of pump–storage hydropower plant prototype | |
CN114893360B (zh) | 风电机组塔筒异常振动识别和运行状态监测方法及系统 | |
CN111695461A (zh) | 一种基于图频特性的雷达杂波智能分类方法 | |
Wang et al. | Effective component extraction for hydraulic pump pressure signal based on fast empirical mode decomposition and relative entropy | |
CN113688502A (zh) | 一种基于随机共振系统的变频信号去噪方法 | |
CN111706455A (zh) | 一种水力机械转轮叶片进口脱流的判定方法 | |
Jeans et al. | Improved current profile criteria for deepwater riser design | |
Gao et al. | Leakage detection of pipeline based on weighted-permutation entropy | |
CN115932773A (zh) | 基于谱形特征的目标角度检测方法、装置、设备及介质 | |
Faghidi et al. | Detection of bearing fault detection from heavily contaminated signals: a higher-order analytic energy operator method | |
CN112632726B (zh) | 一种面向叶轮机械叶片气动弹性模拟的流场重构方法 | |
CN117005975A (zh) | 一种水轮机空化程度表示方法 | |
CN114708885A (zh) | 一种基于声音信号的风机故障预警方法 | |
Liu et al. | Fault detection of rotating machinery based on marine predator algorithm optimized resonance-based sparse signal decomposition and refined composite multiscale fluctuation dispersion entropy | |
CN113611284A (zh) | 语音库构建方法、识别方法、构建系统和识别系统 | |
Zhu et al. | A novel extraction method for useful component of vibration signal combining variational mode decomposition and relative entropy | |
CN113218683A (zh) | 一种基于振动信号的石油井下电动套管切割器故障识别方法 | |
Tromans et al. | A spectral response surface method for calculating crest elevation statistics | |
CN111524036A (zh) | 一种暂态电能质量扰动分类方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |