CN103592152B - 确定模型水轮机转轮叶片进口部位空化的声学方法 - Google Patents
确定模型水轮机转轮叶片进口部位空化的声学方法 Download PDFInfo
- Publication number
- CN103592152B CN103592152B CN201310585565.6A CN201310585565A CN103592152B CN 103592152 B CN103592152 B CN 103592152B CN 201310585565 A CN201310585565 A CN 201310585565A CN 103592152 B CN103592152 B CN 103592152B
- Authority
- CN
- China
- Prior art keywords
- cavitation
- prime
- underwater sound
- piecewise function
- coefficient
- 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
Landscapes
- Hydraulic Turbines (AREA)
Abstract
确定模型水轮机转轮叶片进口部位空化的声学方法涉及模型水轮机转轮叶片进口部位空化判别领域。本发明通过分析表征水声信号低频能量分布趋势的Δf1值随空化系数变化的趋势来判定模型水轮机转轮叶片进口部位处是否发生空化。在获取模型水轮机运行转轮的声波信号后,通过建立水声信号低频能量分布趋势的数学模型,利用发生模型水轮机转轮叶片进口部位空化时,表征水声信号低频能量分布趋势的一次分段函数的分界点恰好为模型水轮机转轮叶片进口部位空化初始发生点这一特性,通过确定该一次分段函数的分界点的空化系数来确定模型水轮机转轮叶片进口部位空化的发生。
Description
技术领域:本发明涉及一种确定模型水轮机转轮叶片进口部位空化发生的声学方法。
背景技术:随着对水轮机内特性研究的进一步深入及用户对水轮机稳定性要求的逐步提高,保证水轮机特别是巨型水轮机在无空化状态下,特别是转轮叶片进口部位处无空化状态下,安全稳定运行已成为考核机组运行状态的一项重要指标。
鉴于在目前的技术水平下,根本无法在现场确定水轮机转轮叶片进口部位处的空化现象何时发生。只能够在模型水轮机上对转轮叶片进口部位处的空化现象进行研究。传统的研究方法为通过布置在模型水轮机导叶与转轮间的专用设备观察模型水轮机转轮叶片进口部位处是否有连续气泡产生,如果有,则表明模型水轮机转轮叶片进口部位处有空化发生,否则,则未发生转轮叶片进口部位空化。由于该方法只能依靠目测,没有一个统一的可以量化的标准,致使模型水轮机转轮叶片进口部位处空化的确定随着观测者及观测位置的不同而产生很大的差异,从而影响到对转轮叶片进口部位空化的确定。因此迫切需要一种可以完全利用测量数据来判定模型水轮机转轮叶片进口部位处空化何时发生的方法。
发明内容:本发明要解决的技术问题是提供一种能够采用计算机程序确定模型水轮机转轮叶片进口部位处空化的声学方法。本发明的技术方案是:一种确定模型水轮机转轮叶片进口部位处空化的声学方法,该方法包括如下步骤:
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;
X*(k)为X(k)的共轭频谱函数,
6)利用一次分段函数来表示水声功率谱中低频能量与高频能量的分布趋势,具体方法如下:
在水轮机水声信号的功率谱中,能量的变化趋势是连续的,且低频区域呈现出随着水声频率的升高能量值陡降的趋势,而在高频区域能量值随着水声频率的升高而降低的幅度则要小得多,低频区域与高频区域之间存在着一个明显的交界点,将水轮机水声信号的功率谱中高低频区域能量随水声频率的变化趋势采用最小二乘法拟合的一次分段函数来拟合,则该一次分段函数可表示为:
式中:
是一次分段函数;
xm是一次分段函数的交点;
k0、k1是一次分段函数中一次项的系数;
b0、b1是一次分段函数中常数项;
7)利用一次分段函数来表示水声功率谱中低频能量的分布趋势,具体方法如下:
利用项6所述方法,在水声功率谱(0,f0)频率区间内对水轮机水声信号的能量变化趋势采用最小二乘法拟合的一次分段函数来描述,其中,f0为利用项6所述方法计算出的水轮机水声信号功率谱中高低频区域能量的交界点所对应的频率,该一次分段函数的交点所对应的频率记为f1,
8)计算Δf1
按下式计算Δf1:
Δf1=f0-f1
9)采集不同空化系数下的水声信号
逐步降低模型水轮机的空化系数,在不同的空化系数下重复步骤3至步骤8,直至水轮机转轮叶片进口部位处出现明显的空化现象;
10)确定模型水轮机转轮叶片进口部位处空化:
具体判定方法如下:
将上述Δf1随空化系数σ的变化趋势采用最小二乘法拟合的一次分段函数y′来拟合,则偏差平方和最小的一次分段函数的交点x′m′即为模型水轮机转轮叶片进口部位处空化初始发生的位置,该一次分段函数y′可表示为:
式中:
y'是一次分段函数;
x′m是一次分段函数y′的交点;
k′0、k′1是一次分段函数y′中一次项的系数;
b′0、b′1是一次分段函数y′中常数项;
令任一组测量Vi=yi-y′i数据(xi,yi)对一次分段函数y'的偏差为Vi,则有Vi=yi-y′i,即:
假设x′i≤x′m时有n1个测量数据,x′i>x′m时有n2个测量数据,即n1+n2=N,则偏差的平方和Qi为:
令即可计算出Δf1随空化系数的变化趋势偏差的平方和Qi,其中最小Qi值所对应的x′m值即为一次分段函数y'的交界点,而该点所对应的空化系数即为模型水轮机转轮叶片进口部位处空化初始发生的位置;
11)关闭模型水轮机试验系统。
技术效果:本发明在判定模型水轮机转轮叶片进口部位处是否发生空化的研究领域首次提出了利用Δf1值随空化系数变化的趋势来判定模型水轮机转轮叶片进口部位处是否发生空化的方法。如图1所示,模型水轮机转轮叶片进口部位处发生空化前,Δf1随着空化系数的减小呈现出小幅平缓单调下降的趋势;空化发生后,则呈现出随着空化系数的减小,Δf1剧烈大幅波动的规律性。上述两种不同趋势间存在着一个明显的分界点,点1。而点1即为模型水轮机转轮叶片进口部位空化发生的初始点。点1所对应的空化系数即为模型水轮机转轮叶片进口部位处空化初始发生时的空化系数。
附图说明:
图1为Δf1值随着空化系数的变化趋势
具体实施方式如下:
一种确定模型水轮机转轮叶片进口部位处空化的声学方法,其操作流程为:
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;
X*(k)为X(k)的共轭频谱函数,
6)利用一次分段函数来表示水声功率谱中低频能量与高频能量的分布趋势,具体方法如下:
在水轮机水声信号的功率谱中,能量的变化趋势是连续的,且低频区域呈现出随着水声频率的升高能量值陡降的趋势,而在高频区域能量值随着水声频率的升高而降低的幅度则要小得多,低频区域与高频区域之间存在着一个明显的交界点,将水轮机水声信号的功率谱中高低频区域能量随水声频率的变化趋势采用最小二乘法拟合的一次分段函数来拟合,则该一次分段函数可表示为:
式中:
是一次分段函数
xm是一次分段函数的交点
k0、k1是一次分段函数中一次项的系数;
b0、b1是一次分段函数中常数项;
7)利用一次分段函数来表示水声功率谱中低频能量的分布趋势,具体方法如下:
利用项6所述方法,在水声功率谱(0,f0)频率区间内对水轮机水声信号的能量变化趋势采用最小二乘法拟合的一次分段函数来描述,其中,f0为利用项6所述方法计算出的水轮机水声信号功率谱中高低频区域能量的交界点所对应的频率。该一次分段函数的交点所对应的频率记为f1。
8)计算Δf1
按下式计算Δf1:
Δf1=f0-f1
9)采集不同空化系数下的水声信号
逐步降低模型水轮机的空化系数,在不同的空化系数下重复步骤3至步骤8,直至水轮机转轮叶片进口部位处出现明显的空化现象;
10)确定模型水轮机转轮叶片进口部位处空化:
具体判定方法如下:
将上述Δf1随空化系数σ的变化趋势采用最小二乘法拟合的一次分段函数y′来拟合,则偏差平方和最小的一次分段函数的交点x′m′即为模型水轮机转轮叶片进口部位处空化初始发生的位置。该一次分段函数y′可表示为:
式中:
y'是一次分段函数
x′m是一次分段函数y′的交点
k′0、k′1是一次分段函数y′中一次项的系数;
b′0、b′1是一次分段函数y′中常数项;
令任一组测量Vi=yi-y′i数据(xi,yi)对一次分段函数y'的偏差为Vi,则有Vi=yi-y′i,即:
假设x′i≤x′m时有n1个测量数据,x′i>x′m时有n2个测量数据,即n1+n2=N。则偏差的平方和Qi为:
令即可计算出Δf1随空化系数的变化趋势偏差的平方和Qi。其中最小Qi值所对应的x′m值即为一次分段函数y'的交界点,而该点所对应的空化系数即为模型水轮机转轮叶片进口部位处空化初始发生的位置。
11)关闭模型水轮机试验系统。
Claims (1)
1.一种确定模型水轮机转轮叶片进口部位处空化的声学方法,其特征是:该方法包括如下步骤:
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;
X*(k)为X(k)的共轭频谱函数,
6)利用一次分段函数来表示水声功率谱中低频能量与高频能量的分布趋势,具体方法如下:
在水轮机水声信号的功率谱中,能量的变化趋势是连续的,且低频区域呈现出随着水声频率的升高能量值陡降的趋势,而在高频区域能量值随着水声频率的升高而降低的幅度则要小得多,低频区域与高频区域之间存在着一个明显的交界点,将水轮机水声信号的功率谱中高低频区域能量随水声频率的变化趋势采用最小二乘法拟合的一次分段函数来拟合,则该一次分段函数可表示为:
式中:
是一次分段函数;
xm是一次分段函数的交点;
k0、k1是一次分段函数中一次项的系数;
b0、b1是一次分段函数中常数项;
7)利用一次分段函数来表示水声功率谱中低频能量的分布趋势,具体方法如下:
利用项6所述方法,在水声功率谱(0,f0)频率区间内对水轮机水声信号的能量变化趋势采用最小二乘法拟合的一次分段函数来描述,其中,f0为利用项6所述方法计算出的水轮机水声信号功率谱中高低频区域能量的交界点所对应的频率,该一次分段函数的交点所对应的频率记为f1,
8)计算Δf1
按下式计算Δf1:
Δf1=f0-f1
9)采集不同空化系数下的水声信号
逐步降低模型水轮机的空化系数,在不同的空化系数下重复步骤3至步骤8,直至水轮机转轮叶片进口部位处出现明显的空化现象;
10)确定模型水轮机转轮叶片进口部位处空化:
具体判定方法如下:
将上述Δf1随空化系数σ的变化趋势采用最小二乘法拟合的一次分段函数y′来拟合,则偏差平方和最小的一次分段函数的交点x′m′即为模型水轮机转轮叶片进口部位处空化初始发生的位置,该一次分段函数y′可表示为:
式中:
y'是一次分段函数;
x′m是一次分段函数y′的交点;
k′0、k′1是一次分段函数y′中一次项的系数;
b′0、b′1是一次分段函数y′中常数项;
令任一组测量Vi=yi-y′i数据(xi,yi)对一次分段函数y'的偏差为Vi,则有Vi=yi-y′i,即:
假设x′i≤x′m时有n1个测量数据,x′i>x′m时有n2个测量数据,即n1+n2=N,则偏差的平方和Qi为:
令即可计算出Δf1随空化系数的变化趋势偏差的平方和Qi,其中最小Qi值所对应的x′m值即为一次分段函数y'的交界点,而该点所对应的空化系数即为模型水轮机转轮叶片进口部位处空化初始发生的位置;
11)关闭模型水轮机试验系统。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310585565.6A CN103592152B (zh) | 2013-11-20 | 2013-11-20 | 确定模型水轮机转轮叶片进口部位空化的声学方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310585565.6A CN103592152B (zh) | 2013-11-20 | 2013-11-20 | 确定模型水轮机转轮叶片进口部位空化的声学方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN103592152A CN103592152A (zh) | 2014-02-19 |
CN103592152B true CN103592152B (zh) | 2016-09-21 |
Family
ID=50082379
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201310585565.6A Active CN103592152B (zh) | 2013-11-20 | 2013-11-20 | 确定模型水轮机转轮叶片进口部位空化的声学方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN103592152B (zh) |
Families Citing this family (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111189915B (zh) * | 2020-01-13 | 2022-08-19 | 明君 | 一种水力机械空化发生的实时判定方法 |
CN111706455A (zh) * | 2020-05-28 | 2020-09-25 | 哈尔滨泽级源码大数据科技有限公司 | 一种水力机械转轮叶片进口脱流的判定方法 |
CN112065629B (zh) * | 2020-08-06 | 2022-01-07 | 西安理工大学 | 一种贯流式水轮机间隙空化初生的检测方法 |
CN112834142B (zh) * | 2020-12-29 | 2023-02-07 | 哈动国家水力发电设备工程技术研究中心有限公司 | 一种确定轴流式模型水轮机转轮叶片空化初始位置的方法 |
CN112818475B (zh) * | 2020-12-31 | 2023-07-21 | 广西科技大学 | 一种水轮机转轮叶片动水压力分析方法 |
CN113155266B (zh) * | 2021-03-08 | 2022-11-01 | 西安理工大学 | 综合振动测试和压力脉动测试的水轮机空化初生判定方法 |
Family Cites Families (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPH0659864B2 (ja) * | 1985-04-08 | 1994-08-10 | 三井造船株式会社 | 水中音の計測装置 |
CN101813512B (zh) * | 2009-12-07 | 2012-05-30 | 哈尔滨电机厂有限责任公司 | 采用计算机程序确定模型水轮机转轮叶片初生空化的声学方法 |
CN101813568B (zh) * | 2009-12-07 | 2011-07-27 | 哈尔滨电机厂有限责任公司 | 一种确定空化初生的判定方法 |
CN102043908A (zh) * | 2010-12-29 | 2011-05-04 | 哈尔滨电机厂有限责任公司 | 计算机确定模型水轮机转轮叶片气化空泡初生的方法 |
-
2013
- 2013-11-20 CN CN201310585565.6A patent/CN103592152B/zh active Active
Also Published As
Publication number | Publication date |
---|---|
CN103592152A (zh) | 2014-02-19 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103592152B (zh) | 确定模型水轮机转轮叶片进口部位空化的声学方法 | |
Boutilier et al. | Parametric study of separation and transition characteristics over an airfoil at low Reynolds numbers | |
CN101487763B (zh) | 一种测量大噪声环境下振动结构频率响应函数的方法 | |
CN104330471B (zh) | 航空结构损伤的Lamb波时变概率模型监测方法 | |
CN101813512B (zh) | 采用计算机程序确定模型水轮机转轮叶片初生空化的声学方法 | |
CN102043908A (zh) | 计算机确定模型水轮机转轮叶片气化空泡初生的方法 | |
CN103631989A (zh) | 一种离心泵流动诱导噪声数值预测方法 | |
Liu et al. | Numerical and experimental studies of hydraulic noise induced by surface dipole sources in a centrifugal pump | |
CN103149276A (zh) | 确定模型水轮机转轮叶片空蚀起始点的方法 | |
Wang et al. | Applicability of eddy viscosity turbulence models in low specific speed centrifugal pump | |
Zhu | Characteristic correlation between propellers cavitating wake and cavitation noise | |
Lian et al. | ERA modal identification method for hydraulic structures based on order determination and noise reduction of singular entropy | |
CN103336060A (zh) | 确定模型水轮机转轮叶片汽蚀发生的方法 | |
CN106907285A (zh) | 一种水轮机泄水锥打2孔位置的确定方法 | |
CN106934134B (zh) | 一种用于表征离心泵空化非定常强度的方法 | |
CN106202906A (zh) | 一种腐蚀疲劳性能表征和寿命估算方法 | |
CN105841714A (zh) | 一种高速x射线脉冲星脉冲轮廓时延测量方法 | |
Sharma et al. | Anisotropy properties of turbulence in flow over seepage bed | |
CN103557933A (zh) | 确定模型水轮机转轮叶片头部空化的声学方法 | |
CN104090143A (zh) | 变电站母线电压幅值瞬时值测量方法 | |
CN103411666B (zh) | 确定模型水轮机转轮密封空化的声学方法 | |
Hur et al. | Propeller loads of large commercial vessels at crash stop | |
Lyu et al. | Numerical investigation of the cavitation noise around a marine propeller focusing on the influence of ventilation | |
CN106321465A (zh) | 一种离心泵叶轮旋转失速的动力学表征方法 | |
CN103604607A (zh) | 确定模型水轮机转轮间隙汽蚀的声学方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant |