CN103424260A - 确定模型水轮机转轮间隙空化的声学方法 - Google Patents
确定模型水轮机转轮间隙空化的声学方法 Download PDFInfo
- Publication number
- CN103424260A CN103424260A CN201310384491XA CN201310384491A CN103424260A CN 103424260 A CN103424260 A CN 103424260A CN 201310384491X A CN201310384491X A CN 201310384491XA CN 201310384491 A CN201310384491 A CN 201310384491A CN 103424260 A CN103424260 A CN 103424260A
- Authority
- CN
- China
- Prior art keywords
- cavitation
- coefficient
- underwater sound
- model
- turbine
- 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.)
- Pending
Links
Images
Landscapes
- Hydraulic Turbines (AREA)
Abstract
Description
技术领域
本发明涉及一种确定模型水轮机转轮间隙空化发生的声学方法。
背景技术
随着对水轮机内特性研究的进一步深入及用户对水轮机稳定性要求的逐步提高,保证水轮机特别是巨型水轮机在无空化状态下,特别是转轮间隙无空化状态下,安全稳定运行已成为考核机组运行状态的一项重要指标。
鉴于在目前的技术水平下,根本无法在现场确定水轮机转轮间隙空化现象何时发生。只能够在模型水轮机上对转轮间隙空化现象进行研究。传统的研究方法为通过透明锥管观察模型水轮机转轮间隙处是否有连续气泡产生,如果有,则表明已发生模型水轮机转轮间隙空化,否则,则未发生间隙空化。由于该方法只能依靠目测,没有一个统一的可以量化的标准,致使模型水轮机转轮间隙空化的确定随着观测者及观测位置的不同而产生很大的差异,从而影响到对间隙空化的确定。因此迫切需要一种可以完全利用测量数据来判定模型水轮机转轮间隙空化何时发生的方法。
发明内容
本发明要解决的技术问题是提供一种能够采用计算机程序确定模型水轮机转轮间隙空化的声学方法。本发明的技术方案是:确定模型水轮机转轮间隙空化的声学方法,其操作流程为:
1)启动模型水轮机试验系统;
2)调整模型水轮机运行工况,使模型水轮机转轮间隙处于未空化状态;
3)保持模型水轮机运行工况稳定,对水声信号进行采集;
4)建立模型水轮机水声信号截取后的时间序列:
x(n)=s(n)W(n)
式中:
x(n)为模型水轮机水声信号截取后的时间序列;
s(n)为模型水轮机水声信号采样后的时间序列;
W(n)为窗函数;
5)计算模型水轮机水声信号的频谱:
式中:
Sxx为模型水轮机水声信号的功率谱;
Δt为模型水轮机水声信号的采样时间间隔;
N为模型水轮机水声信号的采样数;
k=0,1,2,3,4,…,N;
6)利用一次分段函数来表示水声功率谱中低频能量与高频能量的分布趋势,具体方法如下:
在水轮机水声信号的功率谱中,能量的变化趋势是连续的,且低频区域呈现出随着水声频率的升高能量值陡降的趋势,而在高频区域能量值随着水声频率的升高而降低的幅度则要小得多,低频区域与高频区域之间存在着一个明显的交界点,将水轮机水声信号的功率谱中高低频区域能量随水声频率的变化趋势采用最小二乘法拟合的一次分段函数来拟合,则该一次分段函数可表示为:
式中:
b0、b1是一次分段函数中常数项;
7)采集不同空化系数下的水声信号
逐步降低模型水轮机的空化系数,在不同的空化系数下重复步骤3至步骤6,直至水轮机转轮间隙出现明显的空化现象;
8)确定模型水轮机转轮间隙处发生空化:
模型水轮机转轮间隙处发生空化前,一次分段函数中的b1随着空化系数的减小呈现出单调上升的趋势;空化发生后,则呈现出随着空化系数的减小,b1单调下降的规律性。而一次分段函数中的b1最大值处所对应的空化系数即为模型水轮机转轮间隙处空化初始发生时的空化系数。具体判定方法如下:
假定改变空化系数m次,空化系数及其对应的一次分段函数中的b1按空化系数由小到大的顺序分别记为(σ1,(b1)1),(σ2,(b1)2),…,(σm,(b1)m),其中σ为空化系数,对相邻空化系数所对应的b1值进行求差运算:
di=(b1)i+1-(b1)i,i=2,3,…,m
式中:
di是相邻空化系数所对应的b1值之差
当di<0时,表明没有发生模型水轮机转轮间隙空化;
当首次出现di>0时,则点(σi,(b1)i)即为发生模型水轮机转轮间隙空化的位置;
9)关闭模型水轮机试验系统。
本发明在判定模型水轮机转轮间隙是否发生空化的研究领域首次提出了利用表征水声信号能量分布趋势的一次分段函数中的b1值随空化系数变化的趋势来判定模型水轮机转轮间隙是否发生空化的方法,即水轮机在有无空化状态下运行,其水声功率谱均为连续的带频,且低频区域呈现出随着水声频率的升高能量值陡降的趋势,而在高频区域能量值随着水声频率的升高而降低的幅度则要小得多,低频区域与高频区域之间存在着一个明显的交界点。将水轮机水声信号的功率谱中高低频区域能量随水声频率的变化趋势采用最小二乘法拟合的一次分段函数来拟合,则偏差平方和最小的一次分段函数的交点xm即为水声功率谱中低频能量与高频能量的交界点,该一次分段函数可表示为:
式中:
具体步骤为:令每一测量数据(xi,yi)对拟合曲线的偏差为Vi,则有 即:
假设xi≤xm时有n1个测量数据,xi>xm时有n2个测量数据,即n1+n2=N。则偏差的平方和Qi为:
在区间(x0,xn-1)内,以xm=xm+ih递增方式计算不同xm时的Qi值(i=1,2,…,N-1。h为频率分辨率),其中最小Qi值所对应的xm值即为水声功率谱中低频能量与高频能量的交界点频率。
附图说明
具体实施方式:
一种确定模型水轮机转轮间隙空化的声学方法,其操作流程为:
1、启动模型水轮机试验系统;
2、调整模型水轮机运行工况,使模型水轮机转轮间隙处于未空化状态;
3、保持模型水轮机运行工况稳定,对水声信号进行采集;
4、建立模型水轮机水声信号截取后的时间序列:
x(n)=s(n)W(n)
式中:
x(n)为模型水轮机水声信号截取后的时间序列;
s(n)为模型水轮机水声信号采样后的时间序列;
W(n)为窗函数;
5、计算模型水轮机水声信号的频谱:
式中:
Sxx为模型水轮机水声信号的功率谱;
Δt为模型水轮机水声信号的采样时间间隔;
N为模型水轮机水声信号的采样数;
X(k)为对x(n)进行傅立叶变换后的频谱函数,
k=0,1,2,3,4,…,N;
6、利用一次分段函数来表示水声功率谱中低频能量与高频能量的分布趋势,具体方法如下:
在水轮机水声信号的功率谱中,能量的变化趋势是连续的,且低频区域呈现出随着水声频率的升高能量值陡降的趋势,而在高频区域能量值随着水声频率的升高而降低的幅度则要小得多,低频区域与高频区域之间存在着一个明显的交界点,将水轮机水声信号的功率谱中高低频区域能量随水声频率的变化趋势采用最小二乘法拟合的一次分段函数来拟合,则该一次分段函数可表示为:
式中:
7、采集不同空化系数下的水声信号
逐步降低模型水轮机的空化系数,在不同的空化系数下重复步骤3至步骤6,直至模型水轮机转轮间隙出现明显的空化现象;
8、确定模型水轮机转轮间隙空化:
如图1所示,未发生模型水轮机转轮间隙空化时,一次分段函数中的b1值随着空化系数的降低而呈不断升高的趋势;而当发生模型水轮机转轮间隙空化时,b1值随着空化系数的不断减小呈不断降低的趋势。而此过程中一次分段函数中的b1值的最大值,图1中点1的位置,所对应的空化系数即为发生模型水轮机转轮间隙空化时的空化系数,一次分段函数中的b1值的最大值所对应的位置即为模型水轮机转轮间隙空化发生的位置。具体判定方法如下:
假定改变空化系数m次,空化系数及其对应的一次分段函数中的b1按空化系数由小到大的顺序分别记为(σ1,(b1)1),(σ2,(b1)2),…,(σm,(b1)m),其中σ为空化系数,对相邻空化系数所对应的b1值进行求差运算:
di=(b1)i+1-(b1)i,i=2,3,…,m
式中:
di是相邻空化系数所对应的b1值之差
当di<0时,表明没有发生模型水轮机转轮间隙空化;
当首次出现di>0时,则点(σi,(b1)i)即为发生模型水轮机转轮间隙空化的位置。
9、关闭模型水轮机试验系统。
Claims (1)
1.一种确定模型水轮机转轮间隙空化的声学方法,其特征是:
1)启动模型水轮机试验系统;
2)调整模型水轮机运行工况,使模型水轮机转轮间隙处于未空化状态;
3)保持模型水轮机运行工况稳定,对水声信号进行采集;
4)建立模型水轮机水声信号截取后的时间序列:
x(n)=s(n)W(n)
式中:
x(n)为模型水轮机水声信号截取后的时间序列;
s(n)为模型水轮机水声信号采样后的时间序列;
W(n)为窗函数;
5)计算模型水轮机水声信号的频谱:
式中:
Sxx为模型水轮机水声信号的功率谱;
Δt为模型水轮机水声信号的采样时间间隔;
N为模型水轮机水声信号的采样数;
k=0,1,2,3,4,…,N;
6)利用一次分段函数来表示水声功率谱中低频能量与高频能量的分布趋势,具体方法如下:
在水轮机水声信号的功率谱中,能量的变化趋势是连续的,且低频区域呈现出随着水声频率的升高能量值陡降的趋势,而在高频区域能量值随着水声频率的升高而降低的幅度则要小得多,低频区域与高频区域之间存在着一个明显的交界点,将水轮机水声信号的功率谱中高低频区域能量随水声频率的变化趋势采用最小二乘法拟合的一次分段函数来拟合,则该一次分段函数可表示为:
式中:
7)采集不同空化系数下的水声信号;
逐步降低模型水轮机的空化系数,在不同的空化系数下重复步骤3)至步骤6),直至模型水轮机转轮间隙出现明显的空化现象;
8)确定模型水轮机转轮间隙空化:
模型水轮机转轮间隙发生空化前,一次分段函数中的b1随着空化系数的减小呈现出单调上升的趋势;空化发生后,则呈现出随着空化系数的减小,b1单调下降的规律性;而一次分段函数中的b1最大值处所对应的空化系数即为模型水轮机转轮间隙空化初始发生时的空化系数,具体判定方法如下:
假定改变空化系数m次,空化系数及其对应的一次分段函数中的b1按空化系数由小到大的顺序分别记为(σ1,(b1)1),(σ2,(b1)2),…,(σm,(b1)m),其中σ为空化系数,对相邻空化系数所对应的b1值进行求差运算:
di=(b1)i+1-(b1)i,i=2,3,…,m
式中:
di是相邻空化系数所对应的b1值之差;
当di<0时,表明没有发生模型水轮机转轮间隙空化;
当首次出现di>0时,则点(σi,(b1)i)即为发生模型水轮机转轮间隙空化的位置;
9)关闭模型水轮机试验系统。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310384491XA CN103424260A (zh) | 2013-08-29 | 2013-08-29 | 确定模型水轮机转轮间隙空化的声学方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310384491XA CN103424260A (zh) | 2013-08-29 | 2013-08-29 | 确定模型水轮机转轮间隙空化的声学方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN103424260A true CN103424260A (zh) | 2013-12-04 |
Family
ID=49649343
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201310384491XA Pending CN103424260A (zh) | 2013-08-29 | 2013-08-29 | 确定模型水轮机转轮间隙空化的声学方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN103424260A (zh) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112065629A (zh) * | 2020-08-06 | 2020-12-11 | 西安理工大学 | 一种贯流式水轮机间隙空化初生的检测方法 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102043908A (zh) * | 2010-12-29 | 2011-05-04 | 哈尔滨电机厂有限责任公司 | 计算机确定模型水轮机转轮叶片气化空泡初生的方法 |
CN103149276A (zh) * | 2013-03-04 | 2013-06-12 | 哈尔滨电机厂有限责任公司 | 确定模型水轮机转轮叶片空蚀起始点的方法 |
-
2013
- 2013-08-29 CN CN201310384491XA patent/CN103424260A/zh active Pending
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102043908A (zh) * | 2010-12-29 | 2011-05-04 | 哈尔滨电机厂有限责任公司 | 计算机确定模型水轮机转轮叶片气化空泡初生的方法 |
CN103149276A (zh) * | 2013-03-04 | 2013-06-12 | 哈尔滨电机厂有限责任公司 | 确定模型水轮机转轮叶片空蚀起始点的方法 |
Non-Patent Citations (3)
Title |
---|
L. ALFAYEZ 等: "The application of acoustic emission for detecting incipient cavitation and the best efficiency point of a 60 kW centrifugal pump: case study", 《NDT&E INTERNATIONAL》 * |
王辉斌 等: "水轮机故障声学诊断技术研究及其应用探讨", 《大电机技术》 * |
赵越 等: "确定模型水泵水轮机水泵工况叶片初生空化的声学方法研究", 《大电机技术》 * |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112065629A (zh) * | 2020-08-06 | 2020-12-11 | 西安理工大学 | 一种贯流式水轮机间隙空化初生的检测方法 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN102043908A (zh) | 计算机确定模型水轮机转轮叶片气化空泡初生的方法 | |
CN101813512B (zh) | 采用计算机程序确定模型水轮机转轮叶片初生空化的声学方法 | |
CN103592152A (zh) | 确定模型水轮机转轮叶片进口部位空化的声学方法 | |
Van Gent et al. | Stability of rock slopes with shallow foreshores | |
CN102252748B (zh) | 基于经验模态的空化噪声调制特征提取方法 | |
CN103149276A (zh) | 确定模型水轮机转轮叶片空蚀起始点的方法 | |
CN103389341A (zh) | 风力机叶片裂纹检测方法 | |
CN103454537A (zh) | 基于小波分析的风力发电低电压穿越检测设备及方法 | |
CN110410336A (zh) | 一种泵汽蚀状态自动识别方法 | |
CN103336060A (zh) | 确定模型水轮机转轮叶片汽蚀发生的方法 | |
CN114091368B (zh) | 一种轴流式水轮机空化状态的识别方法 | |
CN105067025A (zh) | 一种利用单稳系统随机共振效应检测微弱信号的方法 | |
CN103424183A (zh) | 机械振动信号检测异常干扰消除方法 | |
CN103674235B (zh) | 基于短时傅里叶变换的单一频率报警声音特征检测方法 | |
CN106769054A (zh) | 一种基于声发射信号的水轮机组空蚀空化状态诊断方法 | |
CN103424260A (zh) | 确定模型水轮机转轮间隙空化的声学方法 | |
CN103674550B (zh) | 一种滚动轴承静电监测信号实时混合去噪方法 | |
CN103207931A (zh) | 基于mm和arma算法的次同步振荡模态辨识方法 | |
CN104989633A (zh) | 一种基于仿生小波变换的飞机液压泵故障诊断方法 | |
CN103557933A (zh) | 确定模型水轮机转轮叶片头部空化的声学方法 | |
CN103411666A (zh) | 确定模型水轮机转轮密封空化的声学方法 | |
CN103604607A (zh) | 确定模型水轮机转轮间隙汽蚀的声学方法 | |
CN103411665A (zh) | 确定模型水轮机转轮止漏环空化的声学方法 | |
Lan et al. | Wavelet and improved Hilbert–Huang transform method are used to study the spectrum distribution and energy of turbine pressure pulsation | |
CN103616440A (zh) | 确定模型水轮机转轮止漏环汽蚀的声学方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C02 | Deemed withdrawal of patent application after publication (patent law 2001) | ||
WD01 | Invention patent application deemed withdrawn after publication |
Application publication date: 20131204 |