楊 軍，卿 彪

(國(guó)網(wǎng)四川省電力公司映秀灣水力發(fā)電總廠(chǎng)，成都 611830)

映秀灣水電站位于四川省汶川縣境內(nèi)，坐落于岷江上游左岸，系一低閘隧洞引水式徑流電站。其機(jī)組在經(jīng)過(guò)多年的運(yùn)行之后，機(jī)組效率和水力穩(wěn)定性降低，空化空蝕嚴(yán)重，嚴(yán)重影響了水輪發(fā)電機(jī)組的穩(wěn)定、高效運(yùn)行。因此為了提高機(jī)組性能，有必要對(duì)映秀灣電站機(jī)組轉(zhuǎn)輪進(jìn)行優(yōu)化設(shè)計(jì)。本文先對(duì)轉(zhuǎn)輪進(jìn)行優(yōu)化設(shè)計(jì)，而后采用計(jì)算流體動(dòng)力學(xué)的方法(CFD)對(duì)轉(zhuǎn)輪和活動(dòng)導(dǎo)葉進(jìn)行聯(lián)合計(jì)算，進(jìn)而分析新轉(zhuǎn)輪的指標(biāo)是否符合改造要求。

1 轉(zhuǎn)輪設(shè)計(jì)

1.1 基礎(chǔ)轉(zhuǎn)輪的選擇

現(xiàn)代水輪機(jī)水力設(shè)計(jì)的特點(diǎn)是針對(duì)某一電站的具體參數(shù)條件，首先選擇盡量能夠適應(yīng)或接近于該電站參數(shù)的、性能較好的基礎(chǔ)轉(zhuǎn)輪，在此基礎(chǔ)上應(yīng)用現(xiàn)代水力設(shè)計(jì)手段進(jìn)行轉(zhuǎn)輪的優(yōu)化設(shè)計(jì)[1]，使新設(shè)計(jì)出的轉(zhuǎn)輪的各項(xiàng)性能指標(biāo)全面優(yōu)于原基礎(chǔ)轉(zhuǎn)輪，并且更加適應(yīng)于電站的參數(shù)條件，表1為映秀灣電站原水輪機(jī)的設(shè)計(jì)參數(shù)，原水輪機(jī)轉(zhuǎn)輪為A36。

此次改造的目的是為了提高水輪機(jī)效率并改善機(jī)組的水力穩(wěn)定性。希望在水頭未47～48 m ，水輪機(jī)出力為40 MW左右處于高效率區(qū)運(yùn)行，而改造后水輪機(jī)額定出力不變，同時(shí)提高水輪機(jī)抗磨蝕能力。因此在此次改造過(guò)程中，能更換的只是轉(zhuǎn)輪，導(dǎo)葉或者對(duì)固定導(dǎo)葉進(jìn)行局部的修型處理，同時(shí)還得考慮實(shí)際通道條件的限制。

綜合映秀灣電站的基本參數(shù)以及此次改造的目的，此次轉(zhuǎn)輪的優(yōu)化設(shè)計(jì)選擇的基礎(chǔ)轉(zhuǎn)輪型號(hào)為A982，其能量指標(biāo)，空化指標(biāo)，以及壓力脈動(dòng)等性能指標(biāo)都達(dá)到了很高的水平(見(jiàn)表2)。

1.2 轉(zhuǎn)輪軸面圖的對(duì)比

表1 原水輪機(jī)設(shè)計(jì)參數(shù)Tab1.Design parameters of raw water turbine

表2 A9822轉(zhuǎn)輪基本能量指標(biāo)Tab.2 Basic energy index of A982 runner

A982轉(zhuǎn)輪與A36轉(zhuǎn)輪的軸面圖如圖1所示，前者的導(dǎo)葉高度是0.288，后者的導(dǎo)葉高度是0.25，因此為了滿(mǎn)足映秀灣電站的通道條件，在優(yōu)化過(guò)程當(dāng)中的，將A982轉(zhuǎn)輪的上冠下壓，將導(dǎo)葉高度由0.288下降到0.25，如圖1中粗實(shí)線(xiàn)所示，這就是新轉(zhuǎn)輪A606c轉(zhuǎn)輪的流道圖。新轉(zhuǎn)輪 A606c比 A982的流量減小，基本適合于映秀灣電站參數(shù)要求，A982轉(zhuǎn)輪下環(huán)軸面高度以及轉(zhuǎn)輪出口直徑與A36基本吻合， 故采用 A982做基礎(chǔ)轉(zhuǎn)輪的方案可使映秀灣電站原錐管保持不變。

圖1 軸面比較圖Fig.1 Comparison chart of axial plane

2 數(shù)值模擬

2.1 控制方程

為了使得轉(zhuǎn)輪的計(jì)算更加合理，本次計(jì)算采用轉(zhuǎn)輪與活動(dòng)導(dǎo)葉聯(lián)合計(jì)算的方式，其計(jì)算過(guò)程中涉及的控制方程如下[2-4]：

連續(xù)方程：

(1)

動(dòng)量方程：

(2)

1986年Yakhot和Orszag提出RNGk-ε模型。Renormalization Group(譯為重正化群)的縮寫(xiě)為RNG。修正后的k方程和ε方程形式分別如下：

(4)

(5)

式中：μt=ρCμκ2/ε為湍流黏度；μeff=μ+μt為有效黏性系數(shù)。

μt為分析了旋轉(zhuǎn)葉輪式流體機(jī)械流場(chǎng)效果；Boussines，設(shè)：

2.2 模型建立

根據(jù)改型設(shè)計(jì)目標(biāo)參數(shù)的要求，經(jīng)過(guò)多次計(jì)算，比較改型前后轉(zhuǎn)輪葉片區(qū)域工作面及背面壓力分布、速度分布、流線(xiàn)分布的合理性，計(jì)算比較改型前后轉(zhuǎn)輪數(shù)值效率的高低及轉(zhuǎn)輪出口環(huán)量、最低壓力以及流速變化，優(yōu)化篩選出一個(gè)較優(yōu)的設(shè)計(jì)方案。圖2和圖3分別是轉(zhuǎn)輪聯(lián)合計(jì)算的模型及網(wǎng)格劃分，此次網(wǎng)格劃分采用結(jié)構(gòu)化網(wǎng)格，單周期網(wǎng)格節(jié)點(diǎn)數(shù)約為8萬(wàn)個(gè)。根據(jù)映秀灣電站水輪機(jī)運(yùn)行負(fù)荷加權(quán)因子表的具體情況，在 CFD 計(jì)算時(shí)著重分析了最優(yōu)工況、額定水頭額定出力工況、最低水頭額定出力工況和最大水頭額定出力工況(見(jiàn)表3)。

圖2 計(jì)算區(qū)域Fig.2 Calculation area

圖3 網(wǎng)格劃分Fig.3 Mesh generation

計(jì)算工況單位流量Q11/(m3·s-1)單位轉(zhuǎn)速n11/(r·min-1)最優(yōu)工況0．62069．00額定水頭，額定出力0．73071．10最小水頭，額定出力0．72576．20最大水頭0．55064．30

2.3 模擬結(jié)果及分析

為了充分對(duì)改型轉(zhuǎn)輪與原型轉(zhuǎn)輪進(jìn)行比較，分別對(duì)電站原轉(zhuǎn)輪A36和改型后轉(zhuǎn)輪 A606c 進(jìn)行 CFD 分析計(jì)算。以下給出了電站原轉(zhuǎn)輪 A36 和改型后轉(zhuǎn)輪 A606c 在以上各個(gè)工況的流場(chǎng)分布情況，現(xiàn)僅列舉最優(yōu)工況下的模擬結(jié)果。

圖4 A36和A606c葉片正面壓力分布Fig.4 A36 and A606c blade positive pressure distribution

圖5 A36和A606c葉片背面壓力分布Fig.5 A36 and A606c blade negative pressure distribution

圖6 A36和A606c葉片區(qū)域流線(xiàn)圖壓力分布Fig.6 Pressure distribution of flow line in A36 and A606c blade area

圖7 A36和A606c葉片上冠附近速度矢量Fig.7 Velocity vector near the crown of A36 and A606c blade

圖8 A36和A606c葉片中間流面速度矢量Fig.8 Intermediate velocity vector of A36 and A606c blade

圖9 A36和A606c葉片下環(huán)附近速度矢量Fig.9 Velocity vector near the lower ring of A36 and A606c blade

從圖4～圖9的液體流態(tài)分析可以得出，新轉(zhuǎn)輪 A606c葉葉片工作面和背面壓力分布均勻，特別是背面壓力分布明顯好于 A36 轉(zhuǎn)輪。新轉(zhuǎn)輪葉片頭部正面速度矢量分布均勻，沒(méi)有出現(xiàn)葉片進(jìn)口頭部水流撞 擊產(chǎn)生的回流現(xiàn)象，整個(gè)液道間流線(xiàn)分布合理。新轉(zhuǎn)輪葉片頭部背面速度矢量分布均勻，同樣沒(méi)有出現(xiàn)葉片進(jìn)口頭部水流撞擊產(chǎn)生 的回流現(xiàn)象。通過(guò)多工況 CFD 分析結(jié)果分析可知，葉片正、背面脫流現(xiàn)象在映秀灣水輪機(jī)運(yùn)行范圍內(nèi)基本不會(huì)出現(xiàn)。流態(tài)分布的均勻性將會(huì)對(duì)整個(gè)機(jī)組的穩(wěn)定性非常有利。

圖10的數(shù)值效率比較圖計(jì)算了在額定水頭下，不同流量工況點(diǎn)時(shí)數(shù)值效率。從圖10可以看出，A606c轉(zhuǎn)輪的數(shù)值效率始終高于A36轉(zhuǎn)輪，在數(shù)值上提高了2.5%以上。從圖11可以看出，在額定水頭下，不同流量工況點(diǎn)，A606c的最低壓力均大于相同條件下A36最低壓力，這說(shuō)明新轉(zhuǎn)輪A606c的空化性能要優(yōu)于舊轉(zhuǎn)輪A36[5-7]。

圖10 數(shù)值效率比較圖(n11=71.1 r/min)Fig.10 Numerical efficiency comparison chart(n11=71.1 r/min)

圖11 最低壓力比較圖(n11=71.1 r/min)Fig.11 Minimum pressure comparison chart(n11=71.1 r/min)

3 結(jié) 語(yǔ)

本文先對(duì)映秀灣電站機(jī)組的轉(zhuǎn)輪進(jìn)行了優(yōu)化設(shè)計(jì)，在選擇好合適的基礎(chǔ)轉(zhuǎn)輪之后，根據(jù)映秀灣電站的實(shí)際情況進(jìn)行了修型處理得到了新轉(zhuǎn)輪。隨后利用CFD對(duì)映秀灣水電站新轉(zhuǎn)輪進(jìn)行了計(jì)算，從以上分析可以得出以下結(jié)論。

(1)新轉(zhuǎn)輪A606c各項(xiàng)指標(biāo)(能量指標(biāo)，空化指標(biāo)，穩(wěn)定性等)都優(yōu)于舊轉(zhuǎn)輪A36。

(2)通過(guò)選擇的基礎(chǔ)轉(zhuǎn)輪A982改造后的新轉(zhuǎn)輪A606c滿(mǎn)足映秀灣電站的實(shí)際條件。

(3)此次映秀灣電站的轉(zhuǎn)輪改造達(dá)到了當(dāng)初預(yù)計(jì)的改造要求。

□

[1] 鄭 源，鞠小明，程云山.水輪機(jī)[M].北京:中國(guó)水利水電出版社，2007.

[2] 王福軍.計(jì)算流體動(dòng)力學(xué)分析-CFD軟件原理與應(yīng)用[M].北京：清華大學(xué)出版社,2004.

[3] 謝潔飛,李香桂,楊 輝.基于CFD的離心泵內(nèi)部流場(chǎng)數(shù)值模擬與性能預(yù)測(cè)[J]. 中南林業(yè)科技大學(xué)學(xué)報(bào),2010,30(3):129-132.

[4] 楊敏官,顧海飛,劉 棟,等.離心泵葉輪內(nèi)部湍流流動(dòng)的數(shù)值計(jì)算及試驗(yàn)[J].機(jī)械工程學(xué)報(bào)，2006,42(12):181-185.

[5] 卿 彪，余 波，蘭 崴,等.離心泵流動(dòng)空化噪聲測(cè)試[J].排灌機(jī)械工程學(xué)報(bào)，2016，34(3):198-203.

[6] 卿 彪，任建軍，宋曉飛.基于EMD的離心泵空化初生的研究[J].水力發(fā)電學(xué)報(bào)，2016，35(5):102-109.

[7] 王福霞，趙 麗，付顯威,等.空化初生研究綜述[J].中國(guó)石油大學(xué)勝利學(xué)院學(xué)報(bào),2006,20(3):3-6.