确定模型水轮机转轮叶片汽蚀发生的方法技术

一种确定模型水轮机转轮叶片汽蚀发生的方法涉及水轮机转轮叶片汽蚀发生的判别领域。通过分析水声信号高低频能量在分界点处拟合曲线夹角随空化系数的变化趋势来判断模型水轮机初生空化。在获取水轮机转轮的声波信号后，通过一次分段函数拟合确定水声功率谱中低频能量与高频能量交界处的频率值，并在此基础上计算交界点处分段函数夹角随空化系数的变化趋势，当该夹角值出现局部最大值时，水轮机初生汽蚀系数随之确定。该发明专利技术完成了水轮机汽蚀发生从人工甄别到智能识别的跨越。

【技术实现步骤摘要】

【技术保护点】
一种确定模型水轮机转轮叶片汽蚀发生的方法，其特征是：1）启动计算机系统；2）使模型水轮机处于未空化状态；3）保持模型水轮机运行工况稳定，对水声信号进行采集；4）计算水轮机水声信号的功率谱，具体方法如下：首先，对水声信号进行采样量化与编码，形成采样后的时间序列，采用汉宁（Hanning）窗函数W(n)对采样序列s(n)进行截取，那么加窗后的时间序列x(n)为：x(n)=s(n)W(n)其中x(n)进行截取后的主值序列s(n)进行采样后的时间序列W(n)汉宁（Hanning）窗函数经截取处理的时间序列再经过傅立叶变换得到频谱：X(k)=Σn=0N-1x(n)WNkn式中WNkn=e-j2πkn/NX(k)对x(n)进行傅立叶变换后的频谱函数其共轭为X*(k)=Σn=0N-1x(n)ej2πkn/NX*(k)X(k)的共轭频谱函数则时间序列s(n)的功率谱离散值为Sxx=2ΔtN|X*(k)X(k)|Sxx功率谱的离散值由于功率谱计算具有对称性的特点，输出功率谱线正负谐波关于奈奎斯特（Nyquist）采样频率对称。因此，在计算完功率谱之后，采用单边输出的方式，去掉负谐波；5）确定水声功率谱中高低频能量交界处的频率值，具体方法如下：在水轮机水声信号的功率谱中，幅值的变化趋势是连续的，且低频区域呈现出随着水声频率的升高幅值陡降的趋势，而在高频区域幅值随着水声频率的升高而降低的幅度则要小得多，低频区域与高频区域之间存在着一个明显的交界点，将水轮机水声信号的功率谱中高低频区域幅值随水声频率的变化趋势用采用最小二乘法拟合的一次分段函数来拟合，则偏差平方和最小的一次分段函数的交点xm即为水声功率谱中一次分段函数的交界点；一次分段函数可表示为：y‾=k0x+b0,x≤xmk1x+b1,x>xm式中：一次分段函数xm一次分段函数的交点k0、k1一次分段函数中一次项的系数；b0、b1一次分段函数中常数项；令每一测量数据(xi，yi)对拟合曲线的偏差为Vi，则有即：Vi=yi-k0xi-b0,xi≤xmyi-k1xi-b1,xi>xm假设xi≤xm时有n1个测量数据，xi>xm时有n2个测量数据，即n1+n2=N。则偏差的平方和Qi为：Qi=ΣVi2=Σn1Vi2+Σn2Vi2=Σi=0n1-1[yi-(k0xi+b0)]2+Σi=n1N-1[yi-(k1xi+b1)]2令∂Qi∂k0=0,∂Qi∂k1=0,∂Qi∂b0=0,∂Qi∂b1=0,即可确定一次分段函数和偏差的平方和Qi；在区间(x0,xn?1)内，以xm=x0+ih递增方式计算不同xm时的Qi值（i=1，2，…，N?1。h为频率分辨率），其中最小Qi值所对应的xm值即为水声功率谱中一次分段函数的交界点频率。6）确定水声功率谱中低频区域和高频区域能量梯度及及其比例关系，具体方法如下：a1k0+a2b0=c1a2k0+a3b0=c2a1/k1+a2/b1=c1/a2/k1+a3/b1=c2/水声功率谱中低频能量梯度值k0=Σn1yi×Σn1xi-Σn1yixi×Σn1(Σn1xi)2-Σn1xi2×Σn1水声功率谱中高频能量梯度值k1=Σn2yi×Σn2xi-Σn2yixi×Σn2(Σn2xi)2-Σn2xi2×Σn2则水声功率谱中低频能量与高频能量梯度比例:k0/k1=(Σn1yi×Σn1xi-Σn1yixi×Σn1(Σn1xi)2-Σn1xi2×Σn1)/(Σn2yi×Σn2xi-Σn2yixi×Σn2(Σn2xi)2-Σn2xi2×Σn2)7）确定模型水轮机转轮叶片发生初生空化时的空化系数，在发生空化前,随着空化系数的减小，水声功率谱高低频能量梯...

【技术特征摘要】

【专利技术属性】
技术研发人员：赵越，刘智良，郭全宝，陈金霞，刘登峰，郭彦峰，
申请(专利权)人：哈尔滨电机厂有限责任公司，
类型：发明
国别省市：