CN103411666A - 确定模型水轮机转轮密封空化的声学方法 - Google Patents
确定模型水轮机转轮密封空化的声学方法 Download PDFInfo
- Publication number
- CN103411666A CN103411666A CN2013103846510A CN201310384651A CN103411666A CN 103411666 A CN103411666 A CN 103411666A CN 2013103846510 A CN2013103846510 A CN 2013103846510A CN 201310384651 A CN201310384651 A CN 201310384651A CN 103411666 A CN103411666 A CN 103411666A
- Authority
- CN
- China
- Prior art keywords
- prime
- cavitation
- coefficient
- function
- model
- 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.)
- Granted
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)利用二次分段函数来表示水声功率谱中低频能量与高频能量的分布趋势,具体方法如下:
在模型水轮机水声信号的功率谱中,能量的变化趋势是连续的,且低频区域呈现出随着水声频率的升高能量值陡降的趋势,而在高频区域能量值随着水声频率的升高而降低的幅度则要小得多,低频区域与高频区域之间存在着一个明显的交界点,将水轮机水声信号的功率谱中高低频区域能量随水声频率的变化趋势采用最小二乘法拟合的二次分段函数来拟合,则该二次分段函数可表示为:
式中:
xm是二次分段函数的交点
b0、b1是二次分段函数中常数项;
7)采集不同空化系数下的水声信号
逐步降低模型水轮机的空化系数,在不同的空化系数下重复步骤3至步骤6,直至模型水轮机转轮密封出现明显的空化现象;
8)计算测量数据与二次分段函数偏差的平方和
假设xi≤xm时有n1个测量数据,xi>xm时有n2个测量数据,即n1+n2=N。则偏差的平方和Qi为:
9)确定模型水轮机转轮密封空化的发生位置
发生模型水轮机转轮密封空化前后,二次分段函数与测量数据偏差的平方和Qi随着空化系数的减小均呈现出下降的趋势:模型水轮机转轮密封空化发生前,二次分段函数与测量数据偏差的平方和Qi随着空化系数的减小而降低的幅度较大,而模型水轮机转轮密封空化发生后,二次分段函数与测量数据偏差的平方和Qi随着空化系数的减小而降低的幅度大幅减缓,几乎呈现出一种不随空化系数变化的趋势。这两种变化趋势之间存在着一个明显的交界点,该交界点即为模型水轮机转轮密封空化的发生位置。
将上述二次分段函数与测量数据偏差的平方和Qi随着空化系数的变化趋势采用最小二乘法拟合的一次分段函数来拟合,则偏差平方和最小的一次分段函数的交点x′m′即为模型水轮机转轮密封空化发生的位置。该一次分段函数可表示为:
式中:
具体步骤为:令每一(x′j,y′j)对拟合曲线的偏差为V′i,则有 即:
假设x′i≤x′m′时有n′1个数据,x′i>x′m′时有n′2个数据,即n′1+n′2=N′。则偏差的平方和Q′i为:
令 即可确定上述一次分段函数与Qi偏差的平方和Q′i。在区间(x′0,x′n-1)内,以x′m=x′m+ih′递增方式计算不同x′m时的Q′i值(i=1,2,…,N′-1。h′为频率分辨率),其中Q′i的最小值所对应的x′m值即为一次分段函数的交界点,而该点所对应的空化系数即为模型水轮机转轮密封空化发生时的空化系数。具体判定方法如下:
假定改变空化系数m′次,空化系数及其对应的一次分段函数与Qi偏差的平方和Q′i按空化系数由小到大的顺序分别记为(σ1,Q′1),(σ2,Q′2),…,(σm′,Q′m′),其中σ为空化系数,对相邻空化系数所对应的Q′i值进行求差运算:
d′i=Q′i+1-Q′i,i=1,2,…,m′-1
式中:
d′i是相邻空化系数所对应的Q′i值之差;
当d′i>0时,表明没有发生模型水轮机转轮密封空化;
当首次出现d′i<0时,则点(σi,Q′i)即为发生模型水轮机转轮密封空化的位置;
10)关闭模型水轮机试验系统。
如图1所示,发生模型水轮机转轮密封空化前后,二次分段函数与测量数据偏差的平方和Qi随着空化系数的减小均呈现出下降的趋势:模型水轮机转轮密封空化发生前,二次分段函数与测量数据偏差的平方和Qi随着空化系数的减小而降低的幅度较大,而模型水轮机转轮密封空化发生后,二次分段函数与测量数据偏差的平方和Qi随着空化系数的减小而降低的幅度大幅减缓,几乎呈现出一种不随空化系数变化的趋势。这两种变化趋势之间存在着一个明显的交界点,即图1中的点1,该交界点即为模型水轮机转轮密封空化的发生位置。将上述二次分段函数与测量数据偏差的平方和Qi随着空化系数的变化趋势采用最小二乘法拟合的一次分段函数来拟合,则拟合偏差平方和最小的一次分段函数的交点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;
6)利用二次分段函数来表示水声功率谱中低频能量与高频能量的分布趋势,具体方法如下:
在模型水轮机水声信号的功率谱中,能量的变化趋势是连续的,且低频区域呈现出随着水声频率的升高能量值陡降的趋势,而在高频区域能量值随着水声频率的升高而降低的幅度则要小得多,低频区域与高频区域之间存在着一个明显的交界点,将水轮机水声信号的功率谱中高低频区域能量随水声频率的变化趋势采用最小二乘法拟合的二次分段函数来拟合,则该二次分段函数可表示为:
式中:
xm是二次分段函数的交点
k0、k1是二次分段函数中一次项的系数;
7)采集不同空化系数下的水声信号
逐步降低模型水轮机的空化系数,在不同的空化系数下重复步骤3至步骤6,直至模型水轮机转轮密封出现明显的空化现象;
8)计算测量数据与二次分段函数偏差的平方和
令任一组测量数据(xi,yi)对二次分段函数的偏差为Vi,则有 即:
假设xi≤xm时有n1个测量数据,xi>xm时有n2个测量数据,即n1+n2=N。则偏差的平方和Qi为:
9)确定模型水轮机转轮密封空化的发生位置
发生模型水轮机转轮密封空化前后,二次分段函数与测量数据偏差的平方和Qi随着空化系数的减小均呈现出下降的趋势:模型水轮机转轮密封空化发生前,二次分段函数与测量数据偏差的平方和Qi随着空化系数的减小而降低的幅度较大,而模型水轮机转轮密封空化发生后,二次分段函数与测量数据偏差的平方和Qi随着空化系数的减小而降低的幅度大幅减缓,几乎呈现出一种不随空化系数变化的趋势。这两种变化趋势之间存在着一个明显的交界点,该交界点即为模型水轮机转轮密封空化的发生位置。
将上述二次分段函数与测量数据偏差的平方和Qi随着空化系数的变化趋势采用最小二乘法拟合的一次分段函数来拟合,则偏差平方和最小的一次分段函数的交点x′m′即为模型水轮机转轮密封空化发生的位置。该一次分段函数可表示为:
式中:
具体步骤为:令每一(x′i,y′i)对拟合曲线的偏差为V′i,则有 即:
假设x′i≤x′m′时有n′1个数据,x′i>x′m′时有n′2个数据,即n′1+n′2=N′。则偏差的平方和Q′i为:
令 即可确定上述一次分段函数与Qi偏差的平方和Q′i。在区间(x′0,x′n-1)内,以x′m=x′m+ih′递增方式计算不同x′m时的Q′i值(i=1,2,…,N′-1。h′为频率分辨率),其中Q′i的最小值所对应的x′m值即为一次分段函数的交界点,而该点所对应的空化系数即为模型水轮机转轮密封空化发生时的空化系数。具体判定方法如下:
假定改变空化系数m′次,空化系数及其对应的一次分段函数与Qi偏差的平方和Q′i按空化系数由小到大的顺序分别记为(σ1,Q′i),(σ2,Q′2),…,(σm′,Q′m′),其中σ为空化系数,对相邻空化系数所对应的Q′i值进行求差运算:
d′i=Q′i+1-Q′i,i=1,2,…,m′-1
式中:
d′i是相邻空化系数所对应的Q′i值之差;
当d′i>0时,表明没有发生模型水轮机转轮密封空化;
当首次出现d′i<0时,则点(σi,Q′i)即为发生模型水轮机转轮密封空化的位置;
10)关闭模型水轮机试验系统。
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)计算测量数据与二次分段函数偏差的平方和
假设xi≤xm时有n1个测量数据,xi>xm时有n2个测量数据,即n1+n2=N,则偏差的平方和Qi为:
9)确定模型水轮机转轮密封空化的发生位置
发生模型水轮机转轮密封空化前后,二次分段函数与测量数据偏差的平方和Qi随着空化系数的减小均呈现出下降的趋势:模型水轮机转轮密封空化发生前,二次分段函数与测量数据偏差的平方和Qi随着空化系数的减小而降低的幅度较大,而模型水轮机转轮密封空化发生后,二次分段函数与测量数据偏差的平方和Qi随着空化系数的减小而降低的幅度大幅减缓,几乎呈现出一种不随空化系数变化的趋势,这两种变化趋势之间存在着一个明显的交界点,该交界点即为模型水轮机转轮密封空化的发生位置;
将上述二次分段函数与测量数据偏差的平方和Qi随着空化系数的变化趋势采用最小二乘法拟合的一次分段函数来拟合,则偏差平方和最小的一次分段函数的交点x′m′即为模型水轮机转轮密封空化发生的位置,该一次分段函数可表示为:
式中:
具体步骤为:令每一(x′i,y′i)对拟合曲线的偏差为V′i,则有 即:
假设x′i≤x′m′时有n′1个数据,x′i>x′m′时有n′2个数据,即n′1+n′2=N′。则偏差的平方和Q′i为:
令 即可确定上述一次分段函数与Qi偏差的平方和Q′i,在区间(x′0,x′n-1)内,以x′m=x′m+ih′递增方式计算不同x′m时的Q′i值,其中,i=1,2,…,N′-1,h′为频率分辨率,其中Q′i的最小值所对应的x′m值即为一次分段函数的交界点,而该点所对应的空化系数即为模型水轮机转轮密封空化发生时的空化系数,具体判定方法如下:
假定改变空化系数m′次,空化系数及其对应的一次分段函数与Qi偏差的平方和Q′i按空化系数由小到大的顺序分别记为(σ1,Q′1),(σ2,Q′2),…,(σm′,Q′m′),其中σ为空化系数,对相邻空化系数所对应的Q′i值进行求差运算:
d′i=Q′i+1-Q′i,i=1,2,…,m′-1
式中:
d′i是相邻空化系数所对应的Q′i值之差;
当d′i>0时,表明没有发生模型水轮机转轮密封空化;
当首次出现d′i<0时,则点(σi,Q′i)即为发生模型水轮机转轮密封空化的位置;
10)关闭模型水轮机试验系统。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310384651.0A CN103411666B (zh) | 2013-08-29 | 2013-08-29 | 确定模型水轮机转轮密封空化的声学方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310384651.0A CN103411666B (zh) | 2013-08-29 | 2013-08-29 | 确定模型水轮机转轮密封空化的声学方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN103411666A true CN103411666A (zh) | 2013-11-27 |
CN103411666B CN103411666B (zh) | 2015-06-03 |
Family
ID=49604693
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201310384651.0A Active CN103411666B (zh) | 2013-08-29 | 2013-08-29 | 确定模型水轮机转轮密封空化的声学方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN103411666B (zh) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103616440A (zh) * | 2013-11-20 | 2014-03-05 | 哈尔滨电机厂有限责任公司 | 确定模型水轮机转轮止漏环汽蚀的声学方法 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101813512A (zh) * | 2009-12-07 | 2010-08-25 | 哈尔滨电机厂有限责任公司 | 采用计算机程序确定模型水轮机转轮叶片初生空化的声学方法 |
CN102043908A (zh) * | 2010-12-29 | 2011-05-04 | 哈尔滨电机厂有限责任公司 | 计算机确定模型水轮机转轮叶片气化空泡初生的方法 |
CN103149276A (zh) * | 2013-03-04 | 2013-06-12 | 哈尔滨电机厂有限责任公司 | 确定模型水轮机转轮叶片空蚀起始点的方法 |
-
2013
- 2013-08-29 CN CN201310384651.0A patent/CN103411666B/zh active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101813512A (zh) * | 2009-12-07 | 2010-08-25 | 哈尔滨电机厂有限责任公司 | 采用计算机程序确定模型水轮机转轮叶片初生空化的声学方法 |
CN102043908A (zh) * | 2010-12-29 | 2011-05-04 | 哈尔滨电机厂有限责任公司 | 计算机确定模型水轮机转轮叶片气化空泡初生的方法 |
CN103149276A (zh) * | 2013-03-04 | 2013-06-12 | 哈尔滨电机厂有限责任公司 | 确定模型水轮机转轮叶片空蚀起始点的方法 |
Non-Patent Citations (1)
Title |
---|
赵越等: "确定模型水泵水轮机水泵工况叶片初生空化的声学方法研究", 《大电机技术》 * |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103616440A (zh) * | 2013-11-20 | 2014-03-05 | 哈尔滨电机厂有限责任公司 | 确定模型水轮机转轮止漏环汽蚀的声学方法 |
Also Published As
Publication number | Publication date |
---|---|
CN103411666B (zh) | 2015-06-03 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN102043908A (zh) | 计算机确定模型水轮机转轮叶片气化空泡初生的方法 | |
CN101813512B (zh) | 采用计算机程序确定模型水轮机转轮叶片初生空化的声学方法 | |
CN103592152B (zh) | 确定模型水轮机转轮叶片进口部位空化的声学方法 | |
CN110410336A (zh) | 一种泵汽蚀状态自动识别方法 | |
CN103149276A (zh) | 确定模型水轮机转轮叶片空蚀起始点的方法 | |
CN103454537A (zh) | 基于小波分析的风力发电低电压穿越检测设备及方法 | |
CN107478280A (zh) | 一种基于励磁频率高次谐波分析的水煤浆电磁流量计信号处理方法 | |
CN105403915A (zh) | 基于谱模拟提取地层瞬时吸收衰减参数的方法 | |
CN103644460A (zh) | 一种输气管线泄漏音波信号的滤波优选方法 | |
CN103220241A (zh) | 一种低信噪比条件下信号盒维数特征提取方法 | |
CN102022348B (zh) | 一种水泵汽蚀测量方法 | |
CN104089186A (zh) | 一种基于组合滤波和动态阈值的管道压力异常诊断方法 | |
CN103336060A (zh) | 确定模型水轮机转轮叶片汽蚀发生的方法 | |
CN112065629B (zh) | 一种贯流式水轮机间隙空化初生的检测方法 | |
Fleming et al. | In-situ orifice calibration for reversing oscillating flow and improved performance prediction for oscillating water column model test experiments | |
CN103207931A (zh) | 基于mm和arma算法的次同步振荡模态辨识方法 | |
CN103411666A (zh) | 确定模型水轮机转轮密封空化的声学方法 | |
CN103557933A (zh) | 确定模型水轮机转轮叶片头部空化的声学方法 | |
CN103604607A (zh) | 确定模型水轮机转轮间隙汽蚀的声学方法 | |
CN103424260A (zh) | 确定模型水轮机转轮间隙空化的声学方法 | |
CN103616440A (zh) | 确定模型水轮机转轮止漏环汽蚀的声学方法 | |
CN103778282A (zh) | 一种基于最大横流累积频率的航道有效宽度设计方法 | |
CN103411665A (zh) | 确定模型水轮机转轮止漏环空化的声学方法 | |
CN104239614A (zh) | 压气机气动失稳信号的模拟方法 | |
McNaughton et al. | Experimental testing of the performance and interference effects of a cross-stream array of tidal turbines |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant |