基于魚類友好理念的水輪機(jī)改型及特性分析

趙文龍，周大慶，林奇峰

(河海大學(xué)能源與電氣學(xué)院， 南京 211100)

隨著人們對(duì)于生態(tài)環(huán)境的日益重視和關(guān)注，設(shè)計(jì)建立魚類友好型水電設(shè)施已經(jīng)成為一種趨勢(shì)和研究方向。為解決這一問題，除了修建專門的魚梯、魚道，或進(jìn)行人工增殖放流外[1]，對(duì)于水輪機(jī)組本身的改造使其低傷害或無傷害過魚的能力提升也是極其重要的方式。趙青山等[2]概括了魚類友好型水輪機(jī)設(shè)計(jì)的一般過程并提出魚類友好型水輪機(jī)僅能部分解決下行魚過壩的問題；D．烏帕德亞伊等[3]對(duì)魚類在通過水輪機(jī)時(shí)的損傷情況進(jìn)行分析并引入STRIKER模型來比較改型前后水輪機(jī)對(duì)魚類的傷害概率。但目前對(duì)于魚類友好型水輪機(jī)的研究及其工程實(shí)例仍不多見。故有必要根據(jù)現(xiàn)有的水輪機(jī)組進(jìn)行理論分析和CFD數(shù)值模擬來為今后的工程實(shí)踐奠定基礎(chǔ)。

魚類進(jìn)入水輪機(jī)，損傷和死亡的3個(gè)主要因素是機(jī)械損傷、剪切速度引起的高梯度力以及壓力波動(dòng)[4]，以此得出魚類損傷的具體表現(xiàn)為[5]：葉片撞擊；固定機(jī)械部件擦傷；間隙擦傷；壓降損傷；水力剪切和湍流；空化現(xiàn)象。本文針對(duì)上述葉片撞擊、固定機(jī)械部件擦傷及間隙擦傷3個(gè)表現(xiàn)形式確立以下改型設(shè)計(jì)理念[6]：減少活動(dòng)導(dǎo)葉數(shù)量，合理布置導(dǎo)葉相對(duì)位置；減少和消除轉(zhuǎn)輪室縫隙；轉(zhuǎn)輪入口及出口消除間隙。

本文采用上述理念對(duì)水輪機(jī)組原始模型改型，每一處幾何改型均采取單一變量的方式進(jìn)行改型前后的對(duì)比，并在最后整合所有改型方式，設(shè)計(jì)出最終改型后的水輪機(jī)并利用CFD數(shù)值模擬得出結(jié)論。

1 計(jì)算模型及方法

1.1 模型基本參數(shù)

本文采用軸流式水輪機(jī)模型，水輪機(jī)轉(zhuǎn)輪直徑D=460 mm，設(shè)計(jì)水頭Hr=0.732 m，轉(zhuǎn)速為270 r/min，轉(zhuǎn)輪葉片數(shù)Z1=5，活動(dòng)導(dǎo)葉數(shù)Z2=32。模型試驗(yàn)結(jié)果表明在設(shè)計(jì)水頭下，機(jī)組流量Q=281.97 L/s，出力P=1.63 kW，效率為80.53%。模型包括進(jìn)水流道、蝸殼、固定導(dǎo)葉、活動(dòng)導(dǎo)葉、轉(zhuǎn)輪室以及尾水管，整個(gè)模型流道如圖1所示。

圖1 軸流式水輪機(jī)的模型三維圖Fig.1 Three dimensional model of Kaplan turbine

1.2 計(jì)算域及網(wǎng)格劃分

模型由進(jìn)水流道、蝸殼、導(dǎo)葉、導(dǎo)流裝置、轉(zhuǎn)輪、尾水管幾部分構(gòu)成計(jì)算域。考慮到模型結(jié)構(gòu)復(fù)雜，故采用gambit對(duì)模型劃分適應(yīng)性較強(qiáng)的非結(jié)構(gòu)網(wǎng)格。經(jīng)過網(wǎng)格無關(guān)性驗(yàn)證,初始模型數(shù)值計(jì)算得到的流量為282.22 L/s，出力為1.64 kW，效率為80.87%，與模型試驗(yàn)數(shù)據(jù)十分接近，故確定以上過流部件網(wǎng)格單元數(shù)分別為：27.1、12.2、54.2、15.9、88.7、18 萬個(gè)，共計(jì)216.1 萬個(gè)。

1.3 控制方程

軸流式水輪機(jī)數(shù)值模擬介質(zhì)為水，為不可壓縮的黏性流體。計(jì)算時(shí)采用不可壓縮的連續(xù)方程和N-S方程[7]，方程如下：

▽u=0

(1)

(u▽)u=f-▽p/ρ+υ▽2u

(2)

式中：u為流體速度矢量；f為單位質(zhì)量力矢量；ρ為流體密度；▽為哈密頓算符；p為流體微元上的壓強(qiáng)；υ為流體運(yùn)動(dòng)黏度。

計(jì)算采用Spalart-Allmaras湍流模型，因其網(wǎng)格質(zhì)量要求相對(duì)較低，收斂較快，故在流體機(jī)械的數(shù)值模擬中得以廣泛運(yùn)用。計(jì)算中采用二階迎風(fēng)格式來保證計(jì)算精度，對(duì)湍流流場(chǎng)中的速度和壓力方程用SIMPLEC算法進(jìn)行耦合[8]，計(jì)算時(shí)為確保收斂性，監(jiān)測(cè)殘差值設(shè)為10-5。

1.4 邊界條件

由于已知水輪機(jī)的設(shè)計(jì)水頭，進(jìn)口及尾水管出口分別定義為壓力進(jìn)口和壓力出口邊界條件，轉(zhuǎn)輪設(shè)置為轉(zhuǎn)動(dòng)邊界。壁面設(shè)置為無滑移邊界條件，近壁區(qū)域使用標(biāo)準(zhǔn)壁面函數(shù)的方法。

2 改型及數(shù)值模擬結(jié)果

2.1 導(dǎo)葉改型及數(shù)值模擬結(jié)果

導(dǎo)葉改型旨在減少導(dǎo)葉數(shù)量以及改變固定導(dǎo)葉和活動(dòng)導(dǎo)葉的相對(duì)位置，使其位置、方向趨向一致，減少或避免擾流和機(jī)械傷害[9]，以降低魚類通過導(dǎo)葉時(shí)撞擊卡殼幾率，提高水輪機(jī)組低傷或無傷過魚的能力。導(dǎo)葉改型前后模型及中間截面速度分布如圖2所示。

圖2 導(dǎo)葉改型前后模型及速度分布圖(單位：m/s)Fig.2 The guide vane model and velocity distribution before and after retrofitting

在按照魚類友好理念改型導(dǎo)葉前其中間截面平均流速多處于1.3～1.7 m/s間，改型后其流速多低于1.3 m/s，且流態(tài)穩(wěn)定繞流減少，魚類碰撞和損傷幾率減小。另外，在相同工況下改型后水輪機(jī)流量提升至291.81 L/s，出力為1.66 kW，比改型前提升0.02 kW，同時(shí)其改型后效率為79.11%，降低了1.76%。由此表明這種改型方式在損失少部分效率時(shí)既能保證水輪機(jī)所需出力，又能夠起到保護(hù)魚類的作用。

2.2 輪轂改型及數(shù)值模擬結(jié)果

輪轂改型主要是將傳統(tǒng)的柱形-球形-錐形改為球形輪轂，減少和消除了葉片與輪轂間縫隙并使流道平滑，降低魚類在葉片輪轂間卡殼撞擊的概率從而提高水輪機(jī)組保護(hù)魚類的能力，減少魚類的受傷和死亡幾率。輪轂改型前后的局部示意圖如圖3所示。

圖3 輪轂改型前后轉(zhuǎn)輪室局部正視圖Fig.3 Local graph of runner chamber with different hub shapes

改型后水輪機(jī)流量下降為279.05 L/s，出力下降為1.62 kW，均為小幅度下降，而機(jī)組效率仍維持在80.8%左右。改型前后輪轂及葉片局部壓力分布如圖4所示，從內(nèi)特性上來看并未新增如空化點(diǎn)、壓力明顯變化點(diǎn)等影響轉(zhuǎn)輪運(yùn)行的情況，且輪轂處壓力梯度線分布更為均勻。所以可得出結(jié)論，輪轂由傳統(tǒng)的柱形-球形-錐形改為球形可以適用于水輪機(jī)轉(zhuǎn)輪基于魚類友好理念的改型。

圖4 輪轂改型前后局部壓力分布圖(單位：kPa)Fig.4 Local pressure distribution map with different hub shapes

2.3 輪緣改型及數(shù)值模擬結(jié)果

輪緣改型旨在減少或消除葉片與轉(zhuǎn)輪室間的間隙。將轉(zhuǎn)輪室外壁由傳統(tǒng)的柱形-球形-錐形改為球形并將葉片延伸至轉(zhuǎn)輪室外壁，目的是為了使流道平滑并減少葉片與轉(zhuǎn)輪室壁面的縫隙，從而降低魚類因縫隙擠壓受傷和死亡的幾率。輪緣改型前后二維模型如圖5所示。

圖5 輪緣改型前后轉(zhuǎn)輪室二維示意圖Fig.5 2D runner chamber model with different shroud shapes

改型后在相同工況下機(jī)組流量為287.32 L/s，機(jī)組效率由于減少和消除縫隙后降低了水力損失而有明顯提高，由原始模型的80.87%上升至83.05%，出力也在原始模型的基礎(chǔ)上提升了0.07 kW變?yōu)?.71 kW。改型前后葉片壓力分布如圖6所示，轉(zhuǎn)輪葉片壓力面5 kPa以上高壓區(qū)面積增大，吸力面壓力分布層次更為均勻也可證明其轉(zhuǎn)輪做功能力提高使效率提升?？傻贸鲞@樣的改型方式既滿足了魚類友好理念的設(shè)計(jì)需求，還能在一定程度上提升水輪機(jī)組的水力性能。

圖6 輪緣改型前后葉片壓力分布圖(單位：kPa)Fig.6 Pressure distribution diagram of blade with different shroud shapes

2.4 轉(zhuǎn)輪室進(jìn)出口改型及數(shù)值模擬結(jié)果

轉(zhuǎn)輪室進(jìn)出口改型主要包括修改出水環(huán)形狀以及修改轉(zhuǎn)輪室入口形狀兩方面。轉(zhuǎn)輪室出口改型是將原始的錐形改為從葉片出水端先經(jīng)過球形流道再漸變?yōu)殄F形以減少葉片出口與轉(zhuǎn)輪室壁面間隙。轉(zhuǎn)輪室入口即導(dǎo)葉出口改型主要指將該部位傳統(tǒng)的柱形改為更加平穩(wěn)圓滑的球形和橢球形，以減少導(dǎo)葉懸垂距離從而減少導(dǎo)葉和轉(zhuǎn)輪入口間隙并且有效避免水流泄漏和剪切力對(duì)魚類的傷害。兩種改型可降低魚類進(jìn)出轉(zhuǎn)輪室時(shí)的碰撞卡殼損傷概率。改型前后模型如圖7所示。

圖7 轉(zhuǎn)輪室進(jìn)出口改型前后局部正視圖Fig.7 The inlet and outlet profile of the runner chamber

改型前后模型水力特性如表1所示，在相同工況下，出水環(huán)和轉(zhuǎn)輪入口即導(dǎo)葉出口處經(jīng)過改型后，均不會(huì)對(duì)原始機(jī)組的 水力性能造成過多的影響，所以可按此方式進(jìn)行改型以達(dá)到保護(hù)魚類的目的。

表1 改型前后水輪機(jī)水力特性Tab.1 Turbine’s hydraulic characteristics before and after the reform

2.5 整合改型及數(shù)值模擬結(jié)果

經(jīng)過上述幾種幾何形狀的單一改型，最后整合所有改型于一體。用相同的工況對(duì)其進(jìn)行CFD數(shù)值模擬并與原始模型對(duì)比分析。改型后水輪機(jī)流量小幅提升至284.18 L/s，而效率變?yōu)?7.76%，降低了3.11%，出力為1.59 kW較改型前下降0.05 kW。由此可見，所有單一改型方式的整合并不是完美的搭配，基于魚類友好理念的改型使得原始機(jī)組效率及水力性能有所下降在所難免，但數(shù)值實(shí)驗(yàn)中性能的下降仍然處于可接受范圍內(nèi)，且在改型后并無過度影響原始機(jī)組的性能。

3 結(jié) 語

本文利用CFD數(shù)值模擬對(duì)該軸流式水輪機(jī)進(jìn)行建模并基于魚類友好理念對(duì)機(jī)組進(jìn)行單一改型和整合改型的計(jì)算，得出以下結(jié)論。

(1)導(dǎo)葉改型后出力提升，效率下降1.76%；轉(zhuǎn)輪輪轂和轉(zhuǎn)輪室進(jìn)出口改型后效率、出力和流量基本不變；轉(zhuǎn)輪輪緣改型后出力提高0.07 kW，效率提升2.18%；整合所有單一改型設(shè)計(jì)最終模型后，效率降低3.11%，出力下降0.05 kW。

(2)CFD數(shù)值模擬表明所有改型對(duì)水力特性的影響均在可接受范圍內(nèi)浮動(dòng)，對(duì)魚類友好理念水輪機(jī)的改型具有指導(dǎo)意義。今后還需進(jìn)一步針對(duì)不同魚類的生物特性，研究更多魚類損傷因素并通過更多試驗(yàn)來驗(yàn)證水輪機(jī)組基于魚類友好理念的改型設(shè)計(jì)措施的有效性。

□

[1] 李海鋒. 環(huán)境友好水輪機(jī)流動(dòng)機(jī)理研究[D]. 北京：清華大學(xué), 2002.

[2] 趙青山, 王 顥, 杜 楓,等. 魚類友好型水輪機(jī)的設(shè)計(jì)研究[J]. 電站系統(tǒng)工程, 2014,(2):67-68.

[3] D 烏帕德亞伊, 王振華, 趙建達(dá). 英國(guó)觀察:魚類友好的水輪機(jī)設(shè)計(jì)的最新進(jìn)展[J]. 小水電, 2007,(5):5-7.

[4] Richmond M C, Serkowski J A, Ebner L L, et al. Quantifying barotrauma risk to juvenile fish during hydro-turbine passage[J]. Fisheries Research, 2014,154(154):152-164.

[5] Deng, Zhiqun, Carlson, et al. Fish passage assessment of an advanced hydropower turbine and conventional turbine using Blade-Strike modeling[J]. Energies, 2011,4(1):57-67.

[6] 姜莉萍, 馮順新, 廖文根,等. 環(huán)境友好型水輪機(jī)技術(shù)最新進(jìn)展[C]∥ 水利水電工程新技術(shù)推廣研討會(huì)暨中國(guó)水利學(xué)會(huì)水工結(jié)構(gòu)專業(yè)委員會(huì)第九次年會(huì)，2010.

[7] 尚 勇, 劉小兵, 曾永忠. 葉片數(shù)對(duì)離心泵性能影響的數(shù)值模擬及試驗(yàn)對(duì)比[J]. 中國(guó)農(nóng)村水利水電, 2014,(5):124-127.

[8] 朱 李, 賴喜德. 長(zhǎng)短葉片混流式水輪機(jī)流動(dòng)特性分析[J]. 中國(guó)農(nóng)村水利水電, 2015,(6):162-165.

[9] Robb D. Hydro’s fish-friendly turbines[J]. Renewable Energy Focus, 2011,12(2):16-17.