确定模型水轮机转轮叶片进口部位空化的声学方法

摘要

确定模型水轮机转轮叶片进口部位空化的声学方法涉及模型水轮机转轮叶片进口部位空化判别领域。本发明通过分析表征水声信号低频能量分布趋势的Δf₁值随空化系数变化的趋势来判定模型水轮机转轮叶片进口部位处是否发生空化。在获取模型水轮机运行转轮的声波信号后，通过建立水声信号低频能量分布趋势的数学模型，利用发生模型水轮机转轮叶片进口部位空化时，表征水声信号低频能量分布趋势的一次分段函数的分界点恰好为模型水轮机转轮叶片进口部位空化初始发生点这一特性，通过确定该一次分段函数的分界点的空化系数来确定模型水轮机转轮叶片进口部位空化的发生。

主权项

一种确定模型水轮机转轮叶片进口部位处空化的声学方法，其特征是：该方法包括如下步骤：1）启动模型水轮机试验系统；2）调整模型水轮机运行工况，使模型水轮机转轮叶片进口部位处于未空化状态；3）保持模型水轮机运行工况稳定，对水声信号进行采集；4）建立模型水轮机水声信号截取后的时间序列：x(n)＝s(n)W(n)式中：x(n)为模型水轮机水声信号截取后的时间序列；s(n)为模型水轮机水声信号采样后的时间序列；W(n)为窗函数；5）计算模型水轮机水声信号的频谱：$S_{\mathrm{xx}}=\frac{2\mathrm{\Delta t}}{N}\left|X^{*}\left(k\right)X\left(k\right)\right|$式中：S_xx为模型水轮机水声信号的功率谱；Δt为模型水轮机水声信号的采样时间间隔；N为模型水轮机水声信号的采样数；X(k)为对x(n)进行傅立叶变换后的频谱函数，$W_N^{\mathrm{kn}}=e^{-j2\mathrm{\pi kn}/N};$k＝0,1,2,3,4,…，N；X^*(k)为X(k)的共轭频谱函数，6）利用一次分段函数来表示水声功率谱中低频能量与高频能量的分布趋势，具体方法如下：在水轮机水声信号的功率谱中，能量的变化趋势是连续的，且低频区域呈现出随着水声频率的升高能量值陡降的趋势，而在高频区域能量值随着水声频率的升高而降低的幅度则要小得多，低频区域与高频区域之间存在着一个明显的交界点，将水轮机水声信号的功率谱中高低频区域能量随水声频率的变化趋势采用最小二乘法拟合的一次分段函数来拟合，则该一次分段函数可表示为：$\bar{y}=\left\{\begin{array}{c}{k}_{0}x+b_0,x\le x_m\\ k_{1}x+b_1,x>x_m\end{array}\right.$式中：是一次分段函数；x_m是一次分段函数的交点；k₀、k₁是一次分段函数中一次项的系数；b₀、b₁是一次分段函数中常数项；7）利用一次分段函数来表示水声功率谱中低频能量的分布趋势，具体方法如下：利用项6所述方法，在水声功率谱(0，f₀)频率区间内对水轮机水声信号的能量变化趋势采用最小二乘法拟合的一次分段函数来描述，其中，f₀为利用项6所述方法计算出的水轮机水声信号功率谱中高低频区域能量的交界点所对应的频率，该一次分段函数的交点所对应的频率记为f₁，8）计算Δf₁按下式计算Δf₁：Δf₁＝f₀‑f₁9）采集不同空化系数下的水声信号逐步降低模型水轮机的空化系数，在不同的空化系数下重复步骤3至步骤8，直至水轮机转轮叶片进口部位处出现明显的空化现象；10）确定模型水轮机转轮叶片进口部位处空化：具体判定方法如下:将上述Δf₁随空化系数σ的变化趋势采用最小二乘法拟合的一次分段函数y′来拟合，则偏差平方和最小的一次分段函数的交点x′_m′即为模型水轮机转轮叶片进口部位处空化初始发生的位置，该一次分段函数y′可表示为：$y^{\prime }=\left\{\begin{array}{c}{k}_0^{\prime }{x}^{\prime }+b_0^{\prime },x^{\prime }\le x_m^{\prime }\\ k_1^{\prime }{x}^{\prime }+b_1^{\prime },x^{\prime }>x_m^{\prime }\end{array}\right.$式中：y'是一次分段函数；x′_m是一次分段函数y′的交点；k′₀、k′₁是一次分段函数y′中一次项的系数；b′₀、b′₁是一次分段函数y′中常数项；令任一组测量V_i＝y_i‑y′_i数据(x_i，y_i)对一次分段函数y'的偏差为V_i，则有V_i＝y_i‑y′_i，即：$V_i=\left\{\begin{array}{c}{y}_i-k_0^{\prime }{x}_i^{\prime }-b_0^{\prime },x_i^{\prime }\le x_m^{\prime }\\ y_i-k_1^{\prime }{x}_i^{\prime }-b_1^{\prime },x_i^{\prime }>x_m^{\prime }\end{array}\right.$假设x′_i≤x′_m时有n₁个测量数据，x′_i＞x′_m时有n₂个测量数据，即n₁+n₂＝N，则偏差的平方和Q_i为：令即可计算出Δf₁随空化系数的变化趋势偏差的平方和Q_i，其中最小Q_i值所对应的x′_m值即为一次分段函数y'的交界点，而该点所对应的空化系数即为模型水轮机转轮叶片进口部位处空化初始发生的位置；11）关闭模型水轮机试验系统。