0 引言

抽水蓄能机组利用电力负荷低谷电能抽水至上水库, 储蓄位能, 在电力负荷高峰时输水至水库。在提高系统中火电站和核电站运行效率、延长机组寿命、降低维护费、减少风电场并网运行对电网的冲击、提高风电场和电网运行的协调性等方面发挥着调频调相、旋转备用、黑启动等多重功能^[1,2,3,4,5]。机组在运行过程会经历水泵与水轮机相互转换的过渡过程, 该过程是决定电站稳定运行的关键因素^[6,7,8]。飞逸偏工况运行时, 机组转速增加, 这些动载荷也增大, 同时在无叶区形成高速旋转的高速水环, 加之转轮不能满足进口的无撞击入流条件和法向出流条件, 导致内部出现回环流 (静态旋流) 、流道间的动态旋涡以及旋转失速等复杂的流动现象, 造成流道堵塞^{[9,10,11,12,13]}。其中, 飞逸工况高速旋转的流体充斥整个无叶区形成高速水环, 阻塞流道使其转轮入口无法满足无冲击进口条件, 进而影响机组稳定运行。目前, 尚缺乏对机组出现飞逸过程高速水环的研究, 即对该情况下的水轮机特性了解不够透彻, 妨碍了对机组进行全面原型的实验。

为此, 本文通过数值模拟与实验相结合的方法, 针对实验表现出的外特性不稳定现象, 对其内流场-无叶区高速水环进行数值计算研究, 通过提取计算高速水环区域速度分布、反作用度以及叶片进口冲角, 探究水泵水轮机飞逸工况下无叶区高速水环形成机理, 以及对机组稳定运行的影响。

1 模型建立与网格划分

以某大型抽水蓄能电站单极立轴混流式水泵水轮机模型为研究对象, 如图1所示, 基本参数如表1所示。

图1 模型水泵水轮机计算区域

Fig.1 Computational domains of model pump-turbine

1.蜗壳2.固定导叶3.活动导叶4.转轮5.尾水管

利用ANSYS-ICEM软件进行网格划分。由于偏离最优工况点, 转轮进口的来流不是对称的, 因此需对全流道进行模拟^{[14,15,16,17]}。考虑计算资源等实际条件下, 保证较小的计算量以及较高的计算精度, 故采用六面体结构网格对模型划分, 图2为水泵水轮机各部分的网格划分示意图, 具体各部分网格划分情况如表2所示。

表1 模型水泵水轮机几何参数 下载原表

图2 局部网格示意图

Fig.2 Sketch of local mesh

表2 各部件网格划分 

2 数值模拟方法

2.1 湍流模型与边界条件

飞逸工况在S形特性区, 水泵水轮机在此工况区运行时内部流动为典型的三维、非定常、湍流、强旋流流动, 因此数值计算需选取合适的湍流模型^[18,19,20]。Realizable k-ε湍流模型能够有效地模拟旋转均匀剪切流, 包含射流、混合流的自由流动和分离流动等^[21], Realizable k-ε模型可以更好地模拟表面涡和附壁涡, 对强旋流有一定适应性, 对水泵水轮机湍流压力脉动定性研究具有较高准确性^[22]。

本文主要研究无叶区压力脉动以及脉动和水力激振相互作用对无叶区的影响, 因此选择Realizable k-ε湍流模型进行数值模拟。根据工况所处的不同位置设置具体边界条件, 如表3所示。

2.2 时间步长

本文主要针对转轮叶片与活动导叶之间无叶区高速水环的研究, 由于转轮叶片与活动导叶组成双列叶栅发生动静干涉, 以及偏离最优工况导致无叶区流动复杂, 导致该区域频率分布相对较丰富, 频率分布跨度较大。为了保证压力信号频率的分辨率以及广度, 以0.000 1 s为时间步长, 并以此时间步长连续计算5 s, 对压力信号进行采样, 以保证压力频率分布精度。

表3 边界条件 下载原表

2.3 水力计算结果

以设计水头下模型水泵水轮机为研究对象, 通过实验数据与模拟数据的对比来验证计算模型及计算方法的可靠性。本次选取了飞逸线上7个不同开度的工况点, 如表4所示, 通过已知数据算出其转速n以及进口质量流量q_m, 进行设置计算。

表4 定常数值计算结果

计算出水头H (单位:m) , 可得出单位流量Q₁₁和单位转速n₁₁, 公式为

式中D₂———转轮名义直径, m

Q———计算流量, m³/s

与实验所得飞逸曲线进行对比, 结果如图3所示, 实验结果与数值模拟结果吻合良好, 说明数值模拟方法能较好地反映转轮内部真实流动, 模拟方法准确可靠。

图3 实验与模拟结果对比 下载原图

Fig.3 Comparison of calculated and experimental results

3 计算结果分析

3.1 水泵水轮机模型S特性实验

可逆式水轮机实验平台如图4所示, 对可逆式水泵水轮机模型进行实验研究。此实验台由哈尔滨大电机研究所设计, 已开展了一定的相关实验研究, 针对本课题做一些相应的改进。实验台的循环系统是封闭式的可以双向运行的系统。实验台的验收实验遵循IEC60193、IEC609等相关的规定实验内容。具体实验台参数如表5所示。

图4 可逆式水轮机实验平台 下载原图

Fig.4 Test platform of pump-turbine

本次模型实验选用30 m水头对原型可逆式水泵水轮机各工况进行测试, 获得S区特性曲线如图5所示。

由图5可以发现, 各开度下随着单位转速升高, 流量先增大随后就进入S区。对比不同开度下S特性曲线, 开度越小其S区越明显。这是由于随着单位转速从零开始升高, 水流相对速度方向将会更加接近转轮的进口安放角方向, 故有利于降低水流阻力。因此在这一区间内, 单位流量将随着单位转速的提升而增加。随着转速的进一步增加, 水流相对速度将与骨线之间形成冲角, 有效过流流量减小, 导致单位流量降低。水流对转轮的阻挡作用, 使得其在流量减小的同时也使得转速略有下降, 使得d Q_ed/dn_ed>0 (Q_ed为相对单位流量, n_ed为相对单位转速) , 此后随着单位转速的升高, 水流在转轮进口产生较强的离心力作用, 产生复杂的流动现象, 即水泵水轮机S特性 (图6) 。而且在低流量工况下, 转轮进口与导叶之间发生显著的脱流, 形成高速水环, 对流道形成阻塞作用 (图6a) , 造成流量进一步下降, 因此其S区较明显。由图可以发现小开度飞逸工况下无叶区高速水环均匀充斥着整个无叶区, 更容易对流道形成堵塞。

表5 实验台参数 下载原表

3.2 无叶区高速水环速度分布

3.2.1 反作用度

静压能在总能量中所占比重, 是叶片式机械一个重要的参数, 为此定义

图5 水泵水轮机模型初步实验S区特性曲线 下载原图

Fig.5 Characteristic curves of preliminary test in S zone for pump turbine model

图6 无叶区速度流线图 下载原图

Fig.6 Streamline maps of vaneless region

式中h_p———静压能, Jh_th———总能量头, J

H_p———势扬程, mH_th———扬程, m

p_p———静压升, Pap_th———总压, Pa

Ω为反作用度, 也称为反击系数、反应度或反动度。Ω值表示静压能在叶轮 (转轮) 的总能量中所占的比例。一般来说叶片机械进出口截面速度相差不大, 介质能量变化主要表现在压力或焓的变化。介质与叶轮所交换的动能部分, 最终仍需依赖静止部件转换为静压能变化。所以, 反作用度表现了转轮和静叶内静压能变化的比例。假定叶轮 (转轮) 进出口处周面速度相等, 于是有

式中h_d———动能, Jg———重力加速度, m/s²

H_d———动扬程, m p_d———动压, Pa

ρ———流体密度, kg/m³

c_p———进口绝对速度, m/s

c_s———出口绝对速度, m/s

c_mp———叶轮进口处周面速度, m/s

c_ms———叶轮出口处周面速度, m/s

c_us———出口绝对速度圆周方向分量, m/s

c_up———进口绝对速度圆周方向分量, m/s

式中u_p———进口圆周速度, m/s

u_s———出口圆周速度, m/s

3.2.2 圆周方向速度、反作用度分布

水泵水轮机飞逸工况下, 作为反S区的界限极易发生一系列复杂的水力现象, 在低流量工况下, 转轮进口与导叶之间发生显著的脱流, 形成高速水环, 对流道形成阻塞作用, 造成流量进一步下降, 因此其S区较明显。为了研究无叶区高速水环对机组稳定运行的影响以及对S特性的影响, 通过数值计算, 获得无叶区速度分布, 并沿转轮径向方向提取无叶区速度, 根据速度三角形 (图7) 计算其在圆周方向分量c_u以及其反作用度Ω, 计算结果如图8所示。

图7 径向速度提取图以及速度三角形 下载原图

Fig.7 Radial velocity feature map and velocity triangle

由图8中速度在圆周方向的分量c_u分布图可以发现由点A到点B无叶区流体径向进口方向, 随着A、B间相对距离的增加, 速度在圆周方向的分量c_u呈增加趋势。这说明水泵水轮机飞逸工况下转轮叶片与活动导叶之间无叶区存在高速水环。以转轮与活动导叶动静交接面为分界线, 靠近活动导叶一侧, c_u变化明显呈递增趋势;靠近转轮叶片进口一侧则c_u变化平缓呈现减小趋势。开度a₀为11 mm时, c_u呈现先增加后减小趋势, 随着开度增加减小趋势越加不明显。一方面来流经过活动导叶静压能转化为动能, 所以随着A、B之间距离的增加, c_u增大;另外飞逸工况下水泵水轮机转轮飞速旋转对无叶区扰流促使此区域c_u增加, 但随着A、B之间距离的继续增加, 受到转轮叶片的扰流影响减弱, 逐渐跨过无叶区高速水环, 流动变平缓, 所以c_u在靠近转轮叶片一侧出现先增加后减小趋势。为了进一步探讨转轮叶片扰流形成高速水环的作用机理, 对无叶区反作用度进行计算分析。

由图8中反作用度Ω分布可发现, 在小开度a₀为11 mm时, 靠近导叶一侧, 反作用度Ω随着A、B间距离增大而增大, 即越靠近转轮该区域静压占总能量头的比率越大;通常情况下在导叶中, 静压能转化为动能, 而使得反作用度Ω降低。这说明飞逸工况下无叶区发生强烈的回旋流即存在高速水环, 高速回转的流体使得流道阻塞。流体因此而堆积致使此区域反作用度升高, 靠近转轮一侧Ω变化与开度为21、41 mm变化趋势相同, 均在靠近转轮一侧出现极小值点, 即先减小后增加并且随着开度的增加极值点向B点方向移动, 这是因为流体经过活动导叶受无叶区高速水环的影响, 静压能转化为动能使得此区域反作用度Ω存在一个减小的过程。随着点A、B距离增加, 跨过高速水环区域, 进入叶片流道区域内流动变缓, 因此Ω回升。并且随着流量的增加, 回升点即极值点向叶片流道内移动。

图8 c_u、Ω分布 下载原图

Fig.8 Distributions of c_uandΩ

3.3 无叶区高速水环对叶片来流冲角的影响

常规混流式水轮机和水泵水轮机随着转速从零升高, 水流相对速度方向 (图9) 将会更加接近转轮进口安放角, 故有利于降低水流阻力。因此, 在这一区间内, 单位流量将随着单位转速的升高而增加。随着单位转速的进一步增加, 水流相对速度与叶片进口安放角之间形成冲角, 使得有效过流量减小, 导致单位流量降低。而水泵水轮机, 随着单位转速的升高, 不仅水流与叶片安放角之间形成冲角, 而且水流将在转轮进口产生较强的离心力作用, 产生复杂的物理现象。非设计飞逸工况下, 叶片与来流方向存在冲角, 使得在转轮叶片进口端产生涡, 如图10所示。由叶片进口端涡结构演化图可以发现, 在叶片进口端叶片工作面一侧分布着因来流与叶片型骨线存在冲角而引起的涡结构。而叶片进口端来流冲角不仅与导叶开度有关还受无叶区高速水环的影响。

图9 叶片进口速度三角形 下载原图

Fig.9 Speed triangle at runner channel inlet

图1 0 叶片进口段位置流线图 下载原图

Fig.10 Streamline of inlet of blade

为了研究无叶区高速水环对来流冲角及转轮区域流态的影响, 对叶片进口端来流冲角进行计算, 根据叶片轮廓绘制其翼形骨线, 以叶片型线切线垂直方向建立速度提取线, 根据流体流动速度三角形计算来流冲角, 如图11所示。为了研究来流位置以及非定常对来流冲角的影响, 选取3个转轮叶片在该位置进行计算分析。各位置冲角α随时间变化如图12所示。

由图12可以发现, 水泵水轮机飞逸工况下来流与翼型骨线存在较大冲角, 这也是转轮叶片进口端工作面一侧产生较大尺度涡结构的原因。同时对冲角α变化规律进行分析发现:冲角α随时间变化浮动较大, 3条提取线冲角波动变化并没有在特定的转轮位置进入波峰或者波谷, 同样时间点上也没有发现相同时刻出现增减相同、共同进入峰谷值, 因此冲角α随着流动变化浮动较大, 与来流位置并无明显关系, 存在很强的非定常性和随机性。这是因为转轮飞速旋转与高速水环之间存在速度差, 动静干涉使得在转轮叶片进口端来流方向随机发生改变, 因此导致冲角随时间波动较大。

图11 叶片进口冲角示意图 下载原图

Fig.11 Schematic of blade inlet angle

图1 2 冲角随时间变化图 下载原图

Fig.12 Changing of attack angle with time

综上所述, 来流与叶片骨线存在着较大冲角是转轮叶片进口端工作面一侧产生涡结构的主要原因, 无叶区高速水环与转轮相互作用使得来流冲角随时间波动较大, 波动规律与单一时间点以及来流位置并无明显关系。

4 结论

(1) 水泵水轮机飞逸工况下, 无叶区存在高速旋转水环, 随着开度变大流量增加, 高速水环逐渐被破坏变弱。

(2) 高速水环区周向速度分布以及反作用度分布存在极值点, 其中周向速度在高速水环靠近转轮叶片一侧存在极大值, 反作用度则在靠近转轮叶片一侧存在极小值, 高速水环是引起无叶区水力损失的主要原因。

(3) 水泵水轮机S特性与飞逸点稳定性存在内在的联系, 而飞逸工况下无叶区高速水环影响飞逸的稳定性;飞逸工况下, 随着开度的增加, 高速水环越发不明显, 小开度飞逸工况下, 由于来流方向与叶片骨线存在较大冲角, 致使转轮叶片进口端产生规律性的旋涡结构, 而冲角随着流动变化浮动较大, 存在较大随机性, 高速水环是导致随机性波动的原因。