确定模型水轮机转轮止漏环汽蚀的声学方法.pdf

(19)中华人民共和国国家知识产权局(12)发明专利申请(10)申请公布号(10)申请公布号CNCN103616440103616440A(43)申请公布日2014.03.05(21)申请号201310585673.3(22)申请日2013.11.20(71)申请人哈尔滨电机厂有限责任公司地址150040黑龙江省哈尔滨市香坊区三大动力路99号哈尔滨电机厂有限责任公司技术管理部申请人哈尔滨大电机研究所(72)发明人赵越(51)Int.Cl.G01N29/14(2006.01)权权利要求书3页利要求书3页说明书7页说明书7页附图1页附图1页(54)发明名称确定模型水轮机转轮止漏环汽蚀的声学方法(57)摘要确定模型水轮机转轮止漏环汽蚀的声学方法涉及水轮机汽蚀判别领域。本发明利用表征水声信号功率谱中能量分布规律的三次分段函数与实际测量数据偏差的平方和随汽蚀系数变化的趋势,实现了准确地判定模型水轮机转轮止漏环汽蚀的发生。在获取水轮机转轮的声波信号后，通过分析水声信号功率谱中能量分布规律的三次分段函数与实际测量数据间偏差的平方和随汽蚀系数的变化趋势，采用一次分段函数拟合的方式来确定模型水轮机转轮止漏环汽蚀发生的位置。CN103616440ACN10364ACN103616440A权利要求书1/3页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是三次分段函数的交点；2CN103616440A权利要求书2/3页q0、q1是三次分段函数中三次项的系数；f0、f1是三次分段函数中二次项的系数；k0、k1是三次分段函数中一次项的系数；b0、b1是三次分段函数中常数项；7）采集不同空化系数下的水声信号：逐步降低模型水轮机的空化系数，在不同的空化系数下重复步骤3至步骤6，直至模型水轮机转轮止漏环出现明显的汽蚀现象；8）计算测量数据与三次分段函数偏差的平方和：令任一组测量数据(xi，yi)对三次分段函数的偏差为Vi，则有即：假设xi≤xm时有n1个测量数据，xi＞xm时有n2个测量数据，即n1+n2＝N，则偏差的平方和Qi为：令即可计算出上述三次分段函数与测量数据偏差的平方和Qi；9）确定模型水轮机转轮止漏环汽蚀的发生位置：发生模型水轮机转轮止漏环汽蚀前后，三次分段函数与测量数据偏差的平方和Qi随着汽蚀系数的减小均呈现出下降的趋势：模型水轮机转轮止漏环汽蚀发生前，三次分段函数与测量数据偏差的平方和Qi随着汽蚀系数的减小而降低的幅度较大,而模型水轮机转轮止漏环汽蚀发生后,三次分段函数与测量数据偏差的平方和Qi随着汽蚀系数的减小而降低的幅度大幅减缓，几乎呈现出一种不随汽蚀系数变化的趋势，这两种变化趋势之间存在着一个明显的交界点，该交界点即为模型水轮机转轮止漏环汽蚀的发生位置；将上述三次分段函数与测量数据偏差的平方和Qi随着汽蚀系数的变化趋势采用最小二乘法拟合的一次分段函数来拟合，则偏差平方和最小的一次分段函数的交点即为模型水轮机转轮止漏环汽蚀发生的位置，该一次分段函数可表示为：式中：3CN103616440A权利要求书3/3页是一次分段函数；是一次分段函数的交点；k′0、k′1是一次分段函数中一次项的系数；b′0、b′1是一次分段函数中常数项；具体步骤为：令每一(x′i,y′i)对拟合曲线的偏差为V′i，则有即：假设时有n′1个数据，时有n′2个数据，即n′1+n′2＝N′。则偏差的平方和Q′i为：令即可确定上述一次分段函数与Qi偏差的平方和Q′i，在区间(x′0,x′n-1)内，以x′m＝x′m+ih