映秀湾电站水轮机转轮沙水流动数值模拟

第47卷第3期
2021年3月
水力发电
映秀湾电站水轮机转轮沙水流动数值模拟
吴晓冬1,马元江1,王忠全1,王凯1,赵啸怡#,卢加兴#,刘小兵#
(1.国网四川省电力公司映秀湾水力发电总厂，四川汶川623000；
2.西华大学流体及动力机械教育部重点实验室，四川成都610039)
摘要：为了研究映秀湾水电站水轮机内部沙水流动特性，对水轮机进行沙水流动数值模拟，分析水轮机转轮部分
沙水流动情况#研究结果表明，转轮内叶片进水边及出水边靠下环处是泥沙分布的主要区域#研究结果对水轮机泥
沙磨损的研究具有指导意义#
关键词：水轮机；蜗壳；导水机构；转轮；数值模拟；沙水流动
Numericai Simulation of Sand-Water Flow in Turbine Runner of Yingxinwan Hydropower Station WU Xiaodong1/MA Yuanjiang1/WANG Zhongquan1/WANG Kai1/ZHAO Xiaoyi2,LU Jiaxing2/LII Xiaobing2
(1.Yingxiuwan Hydropower General Plant of State Grid Sichuan Electric Powes Company,Wenchuan623000,
Sichuan,China; 2.Key Laboratoro of Fluid and Powes Machinero(Xihua University)of Ministio of Education,
Chengdu610039,Sichuan,China)
Abstract:In ordes to study the Oow characteristics of sand-water inside the turbine of Yingxiuwan Hydropowes Station,the
sand-water Oow in the turbine is numerical i y simulated,and the Oow of sand-water in the runner of turbine is analyzed.The results show that the sediment is mainly distributed at the inlet and outlet edges of inner blades near the runner band,which haeeguidingsignioicanceooEthestudyoosedimenteEosion ootuEbine.
Key Word::turbine;spirai case;guide-wne apparatus;runner;numerical simulation;sand-water Oow
中图分类号：TK734文献标识码：A文章编号：0559-9342(2021)03-0098-04
0引言
我国是一个幅员辽阔，水资源丰富且多泥沙河流的国家，泥沙磨损问题对水轮机机组运行稳定和工作效率有很大的影响。

因此，国内外学者对泥沙磨损做了大量研究，以期减少泥沙磨损对水轮机机组的危害。

田长安等1在设计工况下对水轮机沙水流动进行数值模拟，通过分析3个不同流面的压力、速度和泥沙体积分数，探索泥沙对水轮机转轮叶片的磨损规律。

黄先北等121对单吸泵进行清水和沙水数值模拟，设定7组泥沙条件方案探究颗粒运动轨迹和磨损规律。

李远余（3）基于N-S方程对高比转速混流式水轮机进行全流道清水和沙水的数值模拟，结果表明数值模拟方法可以很好的模拟水轮机内部流动。

Biraj Singh Thapa等⑷通过CFD的方法模拟了水轮机的侵蚀速率。

田文文等（5）通过数值模拟研究了小流量工况下沙水流动对导叶端面间隙磨损的规律#朱乔琦等何采用标准湍流模型⑺对长短叶片混流式水轮机进行泥沙磨损的研究，得出水轮机流场内压力分布规律和泥沙浓度分布规律。

刘小兵等（8）通过数值模拟和试验相结合的方法，对水轮机进行沙水流动数值模拟以及转轮叶片泥沙磨损试验，研究水轮机内部泥沙磨损具体情况。

在以上研究内容基础上，本文以固液两流体多相流模型为基础，对四川省汶川县映秀湾电站型号为HLA982a-LJ-418的水轮机进行沙水流动数值模拟。

1工程概况
映秀电站位于四川省阿坝藏族羌族自治州汶川县境内的岷江上游左岸，河流多年平均含沙量为
收稿日期：2020-11-10
基金项目：国家重点研发专项项目资助（2018YFB0905200）；国网四川省电力公司科技项目资助（521901180004）
作者简介：吴晓冬（1973―），男，四川仁寿人，硕士研究生，主要从事水电站生产技术管理.
473
第47卷第3期吴晓冬,等:映秀湾电站水轮机转轮沙水流动数值模拟
0. 72 kg/m 3 ,多年平均过机含沙量为0. 33〜0. 37 kg/
m 3 #在2008年汶川大地震(9-0：后，河床结构变化巨 大，电站下游的水位抬升，此外，库区泥沙淤积加 重，水库库容减小，泥沙含量增大，从而造成水轮
机过流部件泥沙磨损严重，使得机组运行效率减少
等一系列问题，对大坝安全直接构成威胁。

2数学模型与方法
2. 1 固液两相流基本运动方程
从固液两相流动瞬时变量的基本运动方程出发，可以写出如下4组方程：
液相连续方程
+—++(—V % 0
(1%
固相连续方程
+#
+C
+— ++ (— V
P
%
(2%
液相动量方程
+#
+C
+ + 1 + P 云"0 ++C $ )VV A )= -〒丘广
(3 %
1 + P A - +
P f
+C
++C (
)VV )=
VI )(分+斜]-冷)(…)+R (
4 %
固 相 动 量方程
+ + 1 + P 云()V )++C ()VV = -册 +
L +C [ ^
(箸 +
箸)]
-*))( V -V ) +)R 0
(5 %
其中
4 = 18 ( 1+40)p L //d ：，) + ) = 1
( 6 %
式中，V 为速度；Q 为材料密度；L 为材质黏性系数； )为体积分数；P 为压强；R 为重力加速度；C 为坐标
分量；4为相间作用系数；d 为颗粒直径。

4项的引
入是为了考虑除s akes 线性阻力作用外的其他作用因 素，一般情况下，4不为常数，它与颗粒雷诺数等流 动参数有关，这里暂假设为常数；液相和固相分别用
下角标/和卩表示；下角标0, J ，7为张量坐标。

2.2湍流模型
运用标准的7y 湍流模型(11-2::
模化湍动能7方程
和)7)++C ()7 v
))(勺+$ 鸽］+
(7%
其中G 7
L # — — +) A — -----------------\
+ _ —
帀-V 2ll )
—+C ()V )
丄$匝—c 7 $ (
v 吨+v
亟)
P $ + C 0 a f \ fi +C 力 +C 0 丿
N _ f +「L 叼)+V + )vi -2 +2 7
N
込［$C (帀 +兀)］一-矿 +
+—
( +
L
+%7 \ + C 0 Vt
4
-
G 4 = ——(——( LL 一 LL )+
Pf
+
(9 %
L # + ) A 「
()-))(V - V )$丽］+
模化能耗率+方程
知)
+) ++C ()+V )=
制)(
"/+$■
)佥
+C 1 +( G 7 + )C 3+］)¥ - )C 2+7
(10%
(11%
式中，C 7 / 0.1 ； C 1 +/ 1.44 ； -+/ 1.92 ； C 3 +/ 1.2 ； $/ 1.0 ； $/ 1.3 ； G 7为由平均速度梯度产生的湍
流动能；G 4为由浮力产生的湍流动能；N 为由于过
渡的扩散产生耗散率的贡献；V ，P ，-为瞬时量；
L ，P ，2表示脉动量；带符号“-”的变量为平均值； 7定义为湍动能；+定义为湍动能耗散率。

2.3边界条件
水轮机设计水头54 m ，最大水头66 m ，最小水 头47 m ，额定流量94. 1m 3/s ，额定功率46. 7 MW ，
转轮进口直径4.18 m ，转速785 rad/min ，共15个 叶片，额定效率94.71%。

按照1：1的尺寸对水体域 进行三维水体建模，如图1所示。

利用ANSYS-
ICEM 软件，对水轮机蜗壳、导叶区、转轮、尾水管 进行非结构化网格划分，网格总数约为2 847.21万 个，网格质量满足计算要求。

图1三维模型示意
参考压力为一个标准大气压101 325 Pa ，进口
边界条件采用速度进口，确定压力出口为出口边界
条件，方向垂直于出口面。

通过计算得：出口边界 所确定的出口压力为60 247. 1 Pa 。

47
3
水力发电2021年3月
3计算结果与分析
压力/Pa
本文研究混流式水轮机在出力42.189 MW ,流
量Q =97.07 m 3/s 大流量工况下，对泥沙含量1.27 Wg/m 3 （体积分数为0. 000 479 %进行沙水数值模拟计 算，进一步分析得出沙水的速度分布，压力分布和 泥沙的体积分布。

(+05 ■(+05
3. 1转轮内部流场计算分析
转轮内部流场计算压力云图如图2 ~4所示。

根 据转轮压力云图可以分析得出转轮吸力面的压力小 于工作面的压力，压力沿着向圆心方向均匀降低,
呈现出梯度式变化，有较为平稳的压力过渡。

转轮 4.168e
3.634e 3.100€ 2.566s 2.032e 1.498e
9.641e
4.3026 -1.037 -6.377
图3转轮下环面压力分布示意
叶片头部区域是最高压力区域，转轮叶片尾部是最 低压区域。

压力/Pa
050505」4.168e+3.634e+
3.100e^05050504
2.566e+2.032e+
1.498e+9.641 e+4.302e+
•1.03沧 -6.3778H
压力/Pa
3.480e+05
2.977e+0S
1.972e+051.469e+059.669e+044.643e+04
a 工作面
b 背面
图2转轮上冠面压力分布示意
图4叶片压力分布云图
转轮近壁面泥沙绕流速度分布如图5所示。

采 用3个转轮叶高特征流面（20%，50%和80% %作为
为研究对象，对近壁面速度作具体分析。

从速度分
布图可以了解到，在各叶高流面转轮的速度变化的
泥沙绕流速度
19.66917.48415.298
■ - 13.113
10.9278.7426.5564.3712.185"0.000
a 20%转轮叶高截面
b 50%转轮叶高截面
c 80%转轮叶高截面
泥沙浓度
0.00070
0.000620.00055
0.00047
0.00039
0.000310.00023
0.000160.000080.00000
图5 转轮近壁面泥沙绕流速度分析
a 20%转轮叶高截面
b 50%转轮叶高截面
c 80%转轮叶高截面
图6转轮叶高流面泥沙浓度
Watrr Powrr Vo.. 47 No.
4
第47卷第3期吴晓冬，等：映秀湾电站水轮机转轮沙水流动数值模拟
分布规律基本一致，在工作面头部到尾部区域内!速度总体变化规律为先减小后增大，叶片背面速度趋势与工作面大体一致，整体速度高于工作面，速度最大的区域分布在80%叶高接近下环尾部区域。

转轮叶高流面泥沙浓度分布见图6,转轮叶片泥沙浓度分布云图见图7。

在叶片进水边缘和叶片的出水边缘靠下环处泥沙体积分数较大，由图可以看出越靠近出水边缘泥沙体积分数越大，这是因为水流经过活动导叶绕流后，向高速旋转的转轮叶片冲击，转轮叶片头部遭到水流强烈的撞击，致使水流里的泥沙在叶片进水边浓度较大，受到重力作用的影响泥沙向下运动，使叶片靠近下环处有泥沙堆积，泥沙体积分数变大，在出水边缘体积分数达到最大。

由此可以推测泥沙体积分数在叶片进水边和出水边处较大，泥沙磨损程度也相对比较严重。

泥沙浓度
0.00087
0.00081
0.00074
0.00068
0.00061
0.00054
0.00048
0.00041
0.00035
a工作面
B Q.00028
图7转轮叶片泥沙浓度分布云图
4结论
本文通过对水轮机进行沙水流动数值模拟，得出如下结论：
(1%模拟分析结果较好的揭示了映秀湾电站型号HLA982vWJW18的水轮机固液两相流内部的沙水流动规律，使电站人员进行此型号水轮机的防护与维修可以起到指导性作用。

(2%泥沙体积分数分布并不均匀，转轮部分的泥沙分布主要是在叶片进水边和叶片的出水边靠下环处。

参考文献：
:1)田长安，袁帅，王俊雄，等•多泥沙河流水电站水轮机转轮内部流动数值模拟:J)•水电能源科学，2020，38(9%:163-166.
:2)黄先北，杨硕，刘竹青，等•基于颗粒轨道模型的离心泵叶轮泥沙磨损数值预测:J)•农业机械学报，2016，47(8%:35-41.
:3)李远余.高比速混流式水轮机固液两相数值模拟:D).哈尔滨:哈尔滨工业大学，2014.
(4)THAPA B S,THAPA B,DAHLHAUG O G.Empirical modelling of
sediment erosion in Franciu turbinee(J).Eneray,2012，41(1%: 386-391.
(5)田文文，刘小兵，卢加兴，等•高水头水电站混流式水轮机导叶
端面空化的数值模拟:J)•中国农村水利水电，2019(8%:211-216，220.
:6)朱乔琦，袁帅，邓聪，等•含沙水中长短叶片混流式水轮机内部流动特性研究:J)•中国农村水利水电，2019(9%:220-223.
:7)唐宁•混流式水轮机涡带脉动的数值模拟研究与试验[D).杭州：浙江大学，2018.
:8)刘小兵，曾永忠，华红，等•渔子溪电站水轮机转轮叶片泥沙磨损研究:J)•西华大学学报(自然科学版％，2020，39(5%:67-73. :9)李宁，唐川，卜祥航，等•“5.12”地震后汶川县泥石流特征与演化分析[J/OL].工程地质学报，(2020-03-18%:2020-10-01).ht-tps：*/kntbrief/default_result.aspx .
:10)韩用顺，王晶，吴淼，等•汶川县震后潜在泥石流危险性评价研究:J).工程科学与技术，2018，50(3%:158-168.
(11)LIU Xiaobing,CHENG Liangjun.A k-+two-equation turbulence
model for solid-liquid two-phase/ows[J).Applied Mathematico and Mechanice，1996(6%:523-531.
(12)LI Xiaobing.Studiee of solid-liquid two-phase turbulent/ow and
wear in hydraulic machinery(J).Journal of Hydrodynamico，1996 (2%:72-76.
(责任编辑高瑜％
)))))))))))))))))))))))))))))))))))))))))))))))))))
(上接第65页)很小，坝体本身形变比大坝周边地区形变小，大坝本身基本没有形变。

需要注意的是，进行调查的影像面积很大，多视数也比较大，对于发现区域内的形变存在一些困难。

三峡大坝周边以丘陵地区为主，受到了几何畸变的影响，对于不同轨道的影响敏感性也不同。

研究区域覆盖范围不够，从形变序列图可以看出，研究区域覆盖范围的边界已经出现形变，说明区域外侧也很有可能存在形变，下一步应该扩大研究区域#同时，研究的重点也更应放在大坝周边的地区，因为坝体相对而言形变较小。

参考文献：
:1)许才军，何平，温扬茂，等•InSAR技术及应用研究进展:J)•测绘地理信息，2015，40(2%:1-9.
:2)刘国祥•InSAR基本原理:J).四川测绘，2004(4%:187-190.
:3)薛廉，吴森，唐侨，等•PS-nSAR技术在成都市地面形变中的监测应用:J)•四川地质学报，2019，39(3%:512-517.
(责任编辑杨健%
473。