丁 華，中村高紀(jì)，德宮健男

?

高比轉(zhuǎn)速混流式水輪機水力優(yōu)化設(shè)計方法

丁 華1，中村高紀(jì)2，德宮健男1

（1. 東芝水電設(shè)備（杭州）有限公司，杭州 310016；2. 株式會社東芝，東京230-0034）

本文主要介紹了有實際運用價值的低水頭、大流量的高比轉(zhuǎn)速混流式水輪機的水力開發(fā)方法。概述了該類混流式水輪機的基本特性、設(shè)計理念，以及綜合運用現(xiàn)代CFD流體解析方法和模型試驗進行水力模型開發(fā)的流程。東芝水電設(shè)備（杭州）有限公司（以下簡稱東芝水電），結(jié)合日本（株式會社）東芝公司（以下簡稱東芝）在水力研究方面擁有多年系統(tǒng)開發(fā)經(jīng)驗，運用經(jīng)大量試驗數(shù)據(jù)驗證的CFD流體解析方法進行性能預(yù)測，以及在日本工業(yè)領(lǐng)域普遍采用的實用優(yōu)化方法，保證了開發(fā)過程的快速、可靠和有效。通過與模型試驗的密切配合，無論新制機組方案還是老機組的改造方案都能夠提供有力的支持，為客戶創(chuàng)造更大的價值。

混流式轉(zhuǎn)輪；高比轉(zhuǎn)速；CFD解析優(yōu)化設(shè)計；試驗參數(shù)法；L18直交法

0 前言

低水頭、大流量的混流式水輪機開發(fā)，對于國民經(jīng)濟的發(fā)展，具有特別的意義。一方面，目前梯級電站開發(fā)過程中，有不少低水頭（設(shè)計水頭低于60m）、大流量的水輪機采用混流式結(jié)構(gòu)。此類機組由于流量較大、葉片較高，轉(zhuǎn)輪內(nèi)部空間扭曲大，流動復(fù)雜，效率、空化、脈動特性要同時滿足客戶要求是項艱巨的任務(wù)。另一方面，從世界范圍來看，很多此類混流式水輪機組已經(jīng)運行多年。它們多數(shù)基于當(dāng)時的技術(shù)進行開發(fā)、設(shè)計并制作。隨著近年來水輪機的水力模型開發(fā)技術(shù)得到飛速的提高，采用最新技術(shù)對老機組進行改造，從而在保證機組穩(wěn)定性、安全性的同時，提高效率、出力，或者減少由于空化、非設(shè)計工況等引起的過流部件的損傷，是目前各大設(shè)備制造廠家關(guān)心的課題。

從國內(nèi)現(xiàn)狀來講，新中國水輪機行業(yè)經(jīng)過60年的發(fā)展，創(chuàng)造了令世人矚目的成績。不同時期的開發(fā)、設(shè)計和制造的機組，體現(xiàn)了當(dāng)時的技術(shù)水平[1]。從水電站經(jīng)濟性的觀點來講，國內(nèi)現(xiàn)有的很多大流量高比轉(zhuǎn)速的水輪機組，通過轉(zhuǎn)輪等關(guān)鍵部件的有效改造，能夠更有效地利用水能，擴大穩(wěn)定高效運行區(qū)間，提高機組的抗空蝕能力及穩(wěn)定性，對電站本身的經(jīng)濟效益，電網(wǎng)的安全供電具有重要意義。

圖1為典型的采用高比轉(zhuǎn)速混流式水輪機的水電站，水流自上庫經(jīng)壓力鋼管進入水輪機帶動轉(zhuǎn)輪旋轉(zhuǎn)，從而推動發(fā)電機運轉(zhuǎn)并送電入電網(wǎng)，經(jīng)過轉(zhuǎn)輪的水再經(jīng)過尾水管回收能量后進入下庫。圖2為東芝水電設(shè)計、制造的典型立式高比轉(zhuǎn)速混流式機組（轉(zhuǎn)輪附近）的剖面圖。

圖1 典型水力發(fā)電站示意圖

圖2 立式高比轉(zhuǎn)速混流式水輪機的構(gòu)造

1 混流式水輪機的基本特性

1.1 性能基本方程式

水力發(fā)電站的上、下庫水位形成的有效水頭，扣除自蝸殼進口經(jīng)過固定導(dǎo)葉、活動導(dǎo)葉到轉(zhuǎn)輪進口，以及自轉(zhuǎn)輪出口到尾水管出口為止的靜止流路內(nèi)發(fā)生的水頭損失loss，實際轉(zhuǎn)輪得到的理論水頭th為：th=–loss。

混流式水輪機轉(zhuǎn)輪作為水輪機組的核心部件，其作用可以基于歐拉方程通過轉(zhuǎn)輪進出口的角運動量的變化來評價。圖3為水輪機運行時轉(zhuǎn)輪進出口的速度系數(shù)三角形，分別表述轉(zhuǎn)輪進出口處的流道內(nèi)的平均流動狀態(tài)。設(shè)轉(zhuǎn)輪進出口對應(yīng)的圓周速度為1、2（m/s），水的絕對流速為1、2（m/s），相對流速為1、2（m/s），各自除以(2)0.5得到的無量綱速度系數(shù)分別為1，2，1，2，1，2。絕對流速與圓周速度的夾角為1、2，相對流速與圓周速度的夾角則為1、2。轉(zhuǎn)輪的理論水頭th可按下式：

th= (11cos1–22cos2) /

圖3 高比轉(zhuǎn)速水輪機進出口速度三角形示意圖

由于轉(zhuǎn)輪內(nèi)部本身存在各種能量損失，轉(zhuǎn)輪內(nèi)水流實際獲得的水頭R要比理論水頭th更低，要扣除轉(zhuǎn)輪內(nèi)部的流動損失，包括壁面引起的摩擦損失f、葉片進口的沖擊損失s、葉片內(nèi)部的減速損失d以及流場不均勻引起的混合損失m，有R=th–f–s–d–m

從而水對轉(zhuǎn)輪的作用功率可按R=RR計算。

式中，為水密度（kg/m3），R為通過轉(zhuǎn)輪的水流量（m3/s）。

對水輪機整體而言，令水輪機進口流量為1，對具備1能量的水流，經(jīng)轉(zhuǎn)輪作用在水輪機的輸出功率為，有=R–Ld–Lm

其中Ld表示上冠、下環(huán)內(nèi)表面的圓盤摩擦損失，Lm為導(dǎo)軸承、密封等位置發(fā)生的機械損失。

令Q=R/1為流量效率，E=R/為比能效率，R= (R–Ld) /R為動力效率，m=/ (R–Ld)為機械效率，則水輪機的效率可表示為：

=/1=m′(R–Ld)/1=m′R′RR/1=R′E′Q′m

對于混流式水輪機而言，效率是反映水輪機能量特性的重要指標(biāo)。由上式可知，水輪機的效率損失由轉(zhuǎn)輪內(nèi)部流動損失、上冠下環(huán)表面圓盤摩擦損失、泄漏損失以及機械損失組成，其中前三項可以表述為水力效率h=R′E′Q，在模型水力開發(fā)階段，可以通過模型試驗，或者CFD流體解析方法預(yù)估水力效率，并通過現(xiàn)代化的優(yōu)化方法來優(yōu)化轉(zhuǎn)輪形狀，從而提高模型水輪機的水力性能。

1.2 高比轉(zhuǎn)速混流式水輪機的特征與設(shè)計理念

具有代表性的高比轉(zhuǎn)速混流式水輪機的模型水力特性曲線如圖4所示。圖中橫軸為單位轉(zhuǎn)速11，縱坐標(biāo)為模型效率，單位出力11以及單位流量11，以模型活動導(dǎo)葉開度為參數(shù)。

圖4 典型高比轉(zhuǎn)速混流式水輪機的模型特性曲線

隨著活動導(dǎo)葉開度的增大，11相應(yīng)地增加，在最大開度附近11的增加量逐漸變小，相應(yīng)地，單位出力11增加量也在變小。實際上若開度超過一定限度繼續(xù)增大，則11反過來會變小（圖中未顯示）。隨開度增大，各開度下的最高效率點所對應(yīng)的單位轉(zhuǎn)速11也在變大。

在一定導(dǎo)葉開度下流量11~11曲線略微朝下，隨著導(dǎo)葉接近最大開度，該曲線近似水平，表明流量幾乎不受轉(zhuǎn)速變化的影響。對高比轉(zhuǎn)速混流式水輪機而言，開度很大時，在高n11處仍能保證很高的效率[2]。

最高效率點附近的轉(zhuǎn)輪進出口的速度三角形如圖3所示。高比轉(zhuǎn)速下活動導(dǎo)葉開度較大，轉(zhuǎn)輪入口的相對液流角1也較大，使得入口處的相對速度C1較大影響了轉(zhuǎn)輪進口形狀。對轉(zhuǎn)輪出口而言，由于需要尾水管回收的水的運動比能22/2較大，考慮到尾水管的設(shè)計，為了避免脫流一般要求設(shè)計葉型時留有一定的旋回速度C2。

高比轉(zhuǎn)速混流式機組常被要求在一定水頭范圍下運行。為了兼顧不同水頭下的運行效率及穩(wěn)定性，一方面應(yīng)當(dāng)將葉片頭部背面的進口空化線排除在最高水頭max之外，并排除葉道渦初生限制線，另一方面，應(yīng)將葉片頭部正面脫流線排除在最低水頭min之外。按工程設(shè)計要求[2]，一般規(guī)定應(yīng)滿足(max-min)/0在30%~40%[3]，其中0為設(shè)計水頭。對實際模型開發(fā)而言，通過設(shè)計葉片特別是頭部附近的厚度分布[4]，和導(dǎo)葉形成有效匹配，改善其寬水頭范圍內(nèi)的水力特性。

2 CFD解析技術(shù)在水力開發(fā)上的運用

無論是新水電站項目水輪機模型的開發(fā)，還是針對現(xiàn)有老機組的改造，都需要在充分把握設(shè)計要求、性能指標(biāo)，分析現(xiàn)有問題和運行條件的前提下，制定經(jīng)過充分論證的性能開發(fā)目標(biāo)，并通過先進的CFD解析技術(shù)和最優(yōu)化方法，開發(fā)設(shè)計出一個性能較優(yōu)的水力模型。之后，通過模型試驗驗證該設(shè)計方案，根據(jù)試驗結(jié)果以此方案為基礎(chǔ)進行形狀微調(diào)后再進行試驗，通常需要反復(fù)數(shù)次后才能得到水力性能較為滿意的方案，該最終方案經(jīng)客戶確認(rèn)后才進入原型機設(shè)計、制作。

對于老機組改造項目，在進入優(yōu)化設(shè)計之前需要先進行性能診斷。包括參考原有機組的模型試驗結(jié)果，現(xiàn)場效率試驗結(jié)果，以及對待改造部件進行專門的CFD流體解析分析。也有省去模型試驗步驟，流體解析結(jié)果直接作為模型最終設(shè)計優(yōu)化的依據(jù)由客戶確認(rèn)的。

開發(fā)目標(biāo)的設(shè)定，一般需要根據(jù)電站的運行情況，從效益、穩(wěn)定性、安全性等方面綜合考慮[5,6]。對于改造方案尺寸已經(jīng)有限制的項目，僅討論轉(zhuǎn)輪的基本參數(shù)的變化，如轉(zhuǎn)輪進口/出口直徑、流路高度、葉片角度的分布等的影響，這樣的例子也不少見。

圖5為東芝水電水力機械研究所現(xiàn)有的水輪機模型開發(fā)的流程示意。（a）和（b）分別表示應(yīng)對新機組開發(fā)和舊機組改造項目的開發(fā)流程。模型開發(fā)的關(guān)鍵指標(biāo)就是使原型機的水力性能滿足客戶要求。為了預(yù)測原型機的水力特性，一般需采用模型試驗測量或CFD解析模擬（改造項目）各運行水頭、流量所對應(yīng)的二元變量下主要的水力特性，包括效率（功率）、空化指標(biāo)和水壓脈動指標(biāo)等。為了完全模擬整個運行范圍內(nèi)的模型水力特性，通常需要在覆蓋整個運行范圍內(nèi)的多個開度下的若干工況點上進行解析或測量。上述獲得的模型特性在指定水頭，按照IEC60193、GB/T15613-2008標(biāo)準(zhǔn)或者客戶指定的其他換算方法換算，從而預(yù)測原型機的水力性能。

圖5 混流式水輪機轉(zhuǎn)輪模型開發(fā)流程圖

目前比較主流的運用于模型開發(fā)的CFD解析模型，兼顧了計算精度和開發(fā)時間上的平衡。靜止部件及轉(zhuǎn)輪全部采用結(jié)構(gòu)化網(wǎng)格。湍流模型目前最常用的為RNG-湍流模型或SST-湍流模型，后者對網(wǎng)格的要求更高一些。

代表性的方法如下[7,8]：

①蝸殼、固定導(dǎo)葉及活動導(dǎo)葉、轉(zhuǎn)輪以及尾水管，考慮水的粘性進行三維湍流CFD數(shù)值解析全流道解析，通過計算全壓差估算水力損失，通過翼面壓力分布或歐拉兩相流模型評估葉片進/出口的空化情況，通過非定常計算，預(yù)測水壓脈動的存在情形。

②部分項目考慮到開發(fā)時間、成本上的限制，對于蝸殼形狀的優(yōu)化，可對蝸殼進行單獨的CFD解析。將蝸殼出口的計算結(jié)果作為入口條件，將固定導(dǎo)葉、活動導(dǎo)葉和轉(zhuǎn)輪聯(lián)合解析優(yōu)化，以考察固定部件、旋轉(zhuǎn)部件之間的耦合作用。對尾水管形狀的優(yōu)化，可以由上述轉(zhuǎn)輪出口計算結(jié)果作為入口條件單獨CFD解析而來，高比轉(zhuǎn)速水輪機為了衡量尾水管在回收轉(zhuǎn)輪出口動能方面的效果，也有將轉(zhuǎn)輪和尾水管聯(lián)合模型進行分析。

③對于泄漏損失以及機械損失，在不同的國際、國家標(biāo)準(zhǔn)內(nèi)有相應(yīng)的經(jīng)驗公式。在實際開發(fā)中，如針對不同的密封形式、形狀，也可對泄漏損失等進行單獨的分析計算。

3 高比轉(zhuǎn)速水輪機水力開發(fā)事例

目前國內(nèi)已建、在建或籌建的水電站中，有不少要求大流量、中低水頭的機組。

下文中，針對某低水頭大流量電站的設(shè)計要求開發(fā)的轉(zhuǎn)輪模型，在高比轉(zhuǎn)速基礎(chǔ)轉(zhuǎn)輪的前提下，如圖6所示。性能提高的目標(biāo)主要為以下4點：

（1）通過改善轉(zhuǎn)輪葉片的全體形狀，使得葉片表面的壓力分布更均勻，特別是下環(huán)外圓側(cè)的葉片形狀的調(diào)整，通過改善此處的流動偏離，使得設(shè)計流量偏向原轉(zhuǎn)輪的過負(fù)荷區(qū)，從而改善加權(quán)平均效率的分布。

（2）通過改變進口及出口附近的葉片角度，改善水輪機的空化特性。

（3）通過加長泄水錐等措施，改善水輪機出口的穩(wěn)定性，從而改善尾水管壓力脈動。

（4）高比轉(zhuǎn)速水輪機的靜止部件的優(yōu)化設(shè)計，考慮如何在大流量工況下仍能保證效率及穩(wěn)定性為主要參考因素。

圖6 優(yōu)化前高比轉(zhuǎn)速混流式水輪機轉(zhuǎn)輪示意圖

另外，高比轉(zhuǎn)速水輪機一般流量較大，提高水力效率的代價往往是惡化空化性能，因此優(yōu)化時應(yīng)當(dāng)綜合考慮效率和空化的影響。另一方面，由于水頭較低，水頭變動幅度大，為了保證運行范圍避開葉片進口空化線，現(xiàn)在一般都會對葉片入口進行加厚的圓頭處理。

3.1 轉(zhuǎn)輪優(yōu)化方法

為了便于討論，減少設(shè)計參數(shù)，先假定轉(zhuǎn)輪以外的其他固定部件尺寸不變，考慮到設(shè)計上的制約，在不改變轉(zhuǎn)輪出口直徑（作為基準(zhǔn)尺寸）和上冠、下環(huán)大致形狀的前提下，對葉片形狀優(yōu)化。葉片數(shù)保持不變。

3.1.1 解析方法

水力優(yōu)化通過CFD解析分析實現(xiàn)。各部件采用自主開發(fā)的網(wǎng)格程序自動生成高質(zhì)量的結(jié)構(gòu)化網(wǎng)格。

一般用于全流道計算時，考慮轉(zhuǎn)輪與尾水管的耦合，須模擬全部葉片的轉(zhuǎn)輪網(wǎng)格。本文在轉(zhuǎn)輪優(yōu)化時，考慮計算時間和資源的限制，先進行固定導(dǎo)葉/活動導(dǎo)葉區(qū)域的CFD解析計算，然后再采用轉(zhuǎn)輪單枚葉片單個流道網(wǎng)格進行定常計算。轉(zhuǎn)輪與尾水管的聯(lián)合計算，用于作為評價綜合性能時的驗算。計算域模型如圖7所示。下文計算中，進口條件采用相應(yīng)工況點的固定導(dǎo)葉、活動導(dǎo)葉區(qū)域CFD流體解析結(jié)果，出口采用壓力邊界條件。差分格式采用相當(dāng)于2次精度的自過濾中心差分格式，湍流模型采用RNG-模型。

（a） 轉(zhuǎn)輪葉片流道網(wǎng)格

（b） 轉(zhuǎn)輪與尾水管的聯(lián)合網(wǎng)格

3.1.2 實驗設(shè)計方法(DOE)

本研究采用實驗設(shè)計方法（DOE）法來對轉(zhuǎn)輪葉片的主要形狀參量進行優(yōu)化。

實驗設(shè)計方法以概率論和數(shù)理統(tǒng)計為理論基礎(chǔ)[9]，是經(jīng)濟、科學(xué)地安排試驗的一項技術(shù)。先確定試驗指標(biāo)和影響因子，再通過直交表來選擇試驗設(shè)計組合，最后可以通過一定的統(tǒng)計分析方法科學(xué)地分析試驗結(jié)果，得到所需的優(yōu)化參數(shù)組合。在本研究中，考慮采用8~10個設(shè)計變量。這些變量包括子午面參數(shù)、葉片進出口流入（出）角及厚度等。運用L18直交表安排這些變量的組合，將可能存在的2′37=4374個待解析方案簡化到18個方案。

3.1.3 轉(zhuǎn)輪性能評價

根據(jù)上節(jié)得到的18個參數(shù)組合，采用三維轉(zhuǎn)輪設(shè)計程序生成相應(yīng)的轉(zhuǎn)輪，根據(jù)第3.1.1節(jié)論述的CFD流體解析方法，即可以對各自的轉(zhuǎn)輪進行性能評價，并預(yù)估靜止部件（蝸殼、固定導(dǎo)葉、活動導(dǎo)葉及尾水管）的水力損失，從而可以對水輪機全體的效率進行預(yù)測。

CFD流體解析的計算工況，一般取包括最優(yōu)效率點導(dǎo)葉開度在內(nèi)的若干開度，在大致覆蓋主要水頭范圍的若干單位轉(zhuǎn)速進行計算。

3.1.4 穩(wěn)健設(shè)計法與反應(yīng)曲面法

根據(jù)3.1.3節(jié)得到的18個設(shè)計方案及相應(yīng)的性能評價結(jié)果（這里，即為水輪機效率Y，為各自水準(zhǔn)對應(yīng)的數(shù)據(jù)個數(shù)），可以分析評價各設(shè)計參數(shù)對性能的相關(guān)影響。

根據(jù)穩(wěn)健設(shè)計法[10]的理念，采用源自通信工程的信噪比SN比作為指標(biāo)，由于效率值希望盡可能地大，作為望大特性的SN比可定義為：

SN比越大，根據(jù)穩(wěn)健設(shè)計的理論，表示該處性能受到周圍環(huán)境影響的敏感度越小，也即性能越穩(wěn)定。通過計算不同參數(shù)各自水準(zhǔn)的SN比，可以清楚地看到不同參數(shù)對最終設(shè)計目標(biāo)（水輪機效率）的影響效果，如圖8所示，橫軸為8個設(shè)計變量各自相應(yīng)的3個水準(zhǔn)，縱軸表示計算得到的SN比，SN比大的變量水準(zhǔn)意味著設(shè)計上更優(yōu)的點。

圖8 設(shè)計參數(shù)影響因子（SN比）圖

根據(jù)反應(yīng)曲面法（RSM）的兩階段優(yōu)化理論[11]，以上述得到的最優(yōu)SN比的參數(shù)為基礎(chǔ)，可以進入第二階段的優(yōu)化，也即在保證SN比不下降的情況下，通過調(diào)整參數(shù)使得敏感度提高進行優(yōu)化設(shè)計。為了加快開發(fā)進程可針對敏感度較高的參數(shù)進行調(diào)整，結(jié)合設(shè)計經(jīng)驗可更快達到設(shè)計結(jié)果。

在本例中，經(jīng)過優(yōu)化后的轉(zhuǎn)輪子午面形狀及葉片形狀如圖9所示。

圖9 優(yōu)化后的轉(zhuǎn)輪葉片形狀

3.1.5 再現(xiàn)性驗證

上面得到的最優(yōu)參數(shù)組合，再經(jīng)過CFD流體解析的驗證，確認(rèn)主要的水力性能得到優(yōu)化，或者再考慮其他因素進行微調(diào)后，即可以作為水力優(yōu)化的模型，進入下一階段的模型試驗階段。

3.2 優(yōu)化后的轉(zhuǎn)輪性能

圖10為水輪機在最優(yōu)工況點附近的轉(zhuǎn)輪葉片表面壓力分布等高線。

圖10 最優(yōu)工況附近轉(zhuǎn)輪的壓力等高線分布

可以清楚地看到，優(yōu)化前葉片入口近下環(huán)位置壓力等高線的疏密變化不均勻，出口位置也出現(xiàn)了局部的壓力梯度過高的區(qū)域，相形之下，優(yōu)化后的轉(zhuǎn)輪該不一致現(xiàn)象得到了明顯的改善，顯示得到了更好的水力形狀。

圖11的最優(yōu)點工況附近的葉片表面流線分布亦顯示，優(yōu)化后流動偏向下環(huán)側(cè)，導(dǎo)致負(fù)荷分布不均勻的問題得到了相當(dāng)程度的改善。

圖11 最優(yōu)工況附近轉(zhuǎn)輪葉片的表面流線分布

在設(shè)計水頭對應(yīng)的部分負(fù)荷（約70%額定出力）下，轉(zhuǎn)輪出口軸方向速度分布的比較如圖12所示。可以看到，即便在部分負(fù)荷下，靠近下環(huán)的流動偏離趨勢仍然得到了一定程度的緩解，流速分布整體也更平緩，由此可以推斷該轉(zhuǎn)輪在部分負(fù)荷下的效率也將略有提高。

圖12 設(shè)計水頭，約70%額定出力下轉(zhuǎn)輪出口無量綱軸向速度分布

轉(zhuǎn)輪葉片表面的壓力系數(shù)分布，如圖13所示，考察了最大出力工況下環(huán)側(cè)附近的沿翼面長度分布的壓力系數(shù)?？梢钥吹剑傮w來講優(yōu)化后的葉片其翼面單位負(fù)荷更小，出口附近完全滿足電站空化系數(shù)p條件不發(fā)生空化。葉片進口在負(fù)壓面附近出現(xiàn)了局部的壓力下降，但仍在允許范圍之內(nèi)，而壓力面?zhèn)鹊膲毫Ψ植嘉蓙y及局部壓力梯度則有明顯改善，可以看到，新設(shè)計轉(zhuǎn)輪的空化性能是可以接受的。由于葉片進口仍然保持圓頭設(shè)計使得水流能沿著流入方向角度緩慢地變化，所以入口的脫流、空化能夠滿足寬水頭范圍內(nèi)運行的要求。

圖13 設(shè)計水頭、最大出力點，下環(huán)附近葉片表面壓力系數(shù)沿流線的分布

優(yōu)化前的模型有完整的模型試驗結(jié)果，將模型試驗結(jié)果、優(yōu)化前的CFD流體解析結(jié)果，以及優(yōu)化后的CFD流體解析結(jié)果，按照IEC60193的相應(yīng)換算方法，在設(shè)計水頭下?lián)Q算成真機運行參數(shù)，如圖14所示。

圖14 優(yōu)化前后基于CFD解析的水輪機性能預(yù)測結(jié)果及與優(yōu)化前水力模型試驗結(jié)果的比較

通過與試驗結(jié)果的對比顯示，CFD解析能夠較為準(zhǔn)確地描述流場的總體特征。優(yōu)化后，水輪機在過負(fù)荷附近的流動得到了改善，表明水輪機在較大流量下的流態(tài)得以優(yōu)化。根據(jù)實際電站設(shè)計需求，該優(yōu)化方案優(yōu)化了其加權(quán)效率系數(shù)較高的工作區(qū)域，從而提高了水輪機年發(fā)電量，將為電站增加可觀的經(jīng)濟效益。

部分負(fù)荷處的效率雖然略有下降，但是根據(jù)轉(zhuǎn)輪出口的流態(tài)分布，與優(yōu)化前差別并不大，優(yōu)化前轉(zhuǎn)輪的模型試驗結(jié)果顯示尾水管內(nèi)由于渦帶引起的振動和噪音都滿足IEC60193標(biāo)準(zhǔn)的要求，轉(zhuǎn)輪的優(yōu)化足以保證部分負(fù)荷區(qū)域?qū)λΨ€(wěn)定的要求。

當(dāng)然，對于部分負(fù)荷區(qū)，定常計算對水輪機性能的預(yù)測有一定的局限性。目前，各家機構(gòu)對水輪機全體非定常計算做了大量研究。特別是對描述空化發(fā)展、壓力脈動引起的穩(wěn)定性分析等[12,13]領(lǐng)域與模型試驗結(jié)果做了大量對比驗證。東芝公司長期致力于發(fā)展高精度的“虛擬模型試驗”CFD解析技術(shù)[14,15]，其精度和可靠性得到了公認(rèn)。非定常計算和模型試驗相結(jié)合，可以為轉(zhuǎn)輪流道做進一步優(yōu)化。

3.3 固定導(dǎo)葉與活動導(dǎo)葉的水力性能優(yōu)化

高比轉(zhuǎn)速的混流式水輪機的固定導(dǎo)葉及活動導(dǎo)葉的優(yōu)化，主要考慮兩點：（1）導(dǎo)葉頭部沖擊損失減?。唬?）是否在全部運行范圍內(nèi)保證自關(guān)閉趨勢。一方面，固定導(dǎo)葉或活動導(dǎo)葉出口的尾流，對水壓脈動等穩(wěn)定特性起關(guān)鍵影響。

圖15中（a）為目前常用的固定導(dǎo)葉和活動導(dǎo)葉設(shè)計，（b）為考慮水力矩和尾流影響后優(yōu)化的設(shè)計結(jié)果。該優(yōu)化設(shè)計不僅考慮了水力損失的減小，也考慮了導(dǎo)葉自關(guān)閉特性和水力穩(wěn)定性的確保，從而確保了原型機運行的可靠性。

（a） 常規(guī)形狀

（b） 優(yōu)化形狀

3.4 尾水管的水力性能優(yōu)化

尾水管的設(shè)計，本身受到土建設(shè)計的制約，同時在部分負(fù)荷下運行通常為了避免振動還要進行補氣，所以設(shè)計時通常要與轉(zhuǎn)輪的設(shè)計一并考慮。圖16中，最優(yōu)工況點附近，優(yōu)化后尾水管比常規(guī)的高窄型尾水管深度約低15%，在具備支墩增加水力損失的前提下，水力性能基本保持不變，從而有助于客戶減少開挖量，提高經(jīng)濟效益。

（a） 高窄型常規(guī)尾水管

（b） 新設(shè)計尾水管

4 結(jié)論

本文闡述了高比轉(zhuǎn)速混流式水輪機水力優(yōu)化設(shè)計的主要方法。根據(jù)目前先進的CFD流體解析技術(shù)和優(yōu)化方法，得到所需的主要部件的水力模型，為下一步進入模型試驗驗證和原型機制造打好堅實的基礎(chǔ)。東芝水電采用的模型開發(fā)方法經(jīng)過日本東芝水力研究所多年驗證，是保證向客戶提供穩(wěn)定、可靠、高效的水輪機的重要一環(huán)，作為設(shè)計工作者的重要工具在新機組設(shè)計、老機組改造等項目中將發(fā)揮越來越大的作用。

[1] 王春暖. 優(yōu)秀混流式水輪機轉(zhuǎn)輪的應(yīng)用綜述[J].水力發(fā)電, 2009, 35(7).

[2] 日本ターボ機械協(xié)會. ハイドロタービン[M]. 日本工業(yè)出版社, 2007.

[3] 黃源芳, 等. 原型水輪機運行研究[M]. 中國電力出版社, 2010.

[4] 劉勝柱, 等. 三峽右岸水輪機水力性能優(yōu)化設(shè)計[J]. 大電機技術(shù), 2004,(1).

[5] Abe Y et al.: Technologies for Rehabilitation of Hydraulic Runners[J].Toshiba Review, 2007, 62(6).

[6] Miyakawa K, et al. Updated re-powering technology for old hydraulicturbine by exchanging runner[J].Mitsubishi Heavy Industries Technical Review, 2004, 41(3).

[7] Enomoto Y, et al. Design Optimization of a High Specific Speed Francis Turbine Runner using Multi-Objective Genetic Algorithm 23rd IAHR Symp[C]. Yokohama, Japan, 2006.

[8] P Drtina, et al. Hydraulic turbines – basic principles and state-of-the-art computational fluid dynamics applications, ProcInstnMechEngrs Vol 213 Part C.

[9] Douglas C, Montgomery, John Wiley & Sons. Inc. Design and Analysis of Experiments 4th edition[M].1996.

[10] 田口玄一.第3版実験計畫法(上?下)[M].丸善, 1976.

[11] Yoshitomo I, et al. The optimization of the diaphragm shape based on the quality engineering[J]. Pioneer R&D ,Vol.12(3): 39-46.

[12] Shi Q H. Experimental Investigation of Upper part Load pressure Pulsations for Three Gorges Model Turbine Proc.24th IAHR Symp[C]. Parana, Brazil, 2008.

[13] Wen-Tao Su, et al. On the flow Instabilities and Turbulent Kinetic Energy of Large-Scale Francis Hydro turbine Model at Low Flow Rate Conditions[J]. Advances in Mechanical Engineering, 2014.

[14] Kurosawa S, et al. Numerical Prediction of Critical Cavitation Performance in Hydraulic Turbines, ASME/JSME 2003 4thJoint Fluids Summer Engineering Conference[C], Vol. 2.

[15] Kurosawa S, et a. Virtual model test for a Francis turbine, 25thIAHR Symposium on Hydraulic Machinery and Systems[C]. 2010 IOP Conf. Ser.: Earth Environ. Sci. 12 012063.

The Hydraulic Design and Optimization Method of a High-Specific-Speed Francis Hydro Turbine

DING Hua1, NAKAMURA Takanori2, TOKUMIYA Takeo1

(1. Toshiba Hydro Power (Hangzhou) Co., Ltd., Hangzhou 310016, China; 2. Toshiba Corporation, Tokyo, 230-0034, Japan)

The present paper discusses the basic developing methods of the hydraulic design of high-specific-speed Francis turbines with a relatively low water head and large discharge. The design points and concepts of such type of turbine are summarized, with emphasis on use of computational fluid dynamics (CFD) as a tool to help guide model tests and gain insight into the complex flowing phenomena inside the turbine. The Toshiba Hydro Power (Hangzhou) Co., Ltd. (THPC) takes full advantage of TOSHIBA Corporation’s strong expertise in hydraulic designing and model developing, using a method to integrate high-accuracy CFD techniques validated by experimental results and optimization solutions that are widely used in industrial applications in Japan. These efforts prove to be reliable and effective, which can be a powerful tool to help the design engineer better understand the brand-new designs or technical renovations. The paper finally gives an example on the detailed discussion of the design and optimization of a typical high-specific-speed Francis turbine, which highlights the state-of-art methods in our hydraulic laboratory.

Francis runner; high specific speed; CFD analysis optimization; DOE; L18 orthogonal array

TK730.2

A

1000-3983(2016)05-0030-08

2015-10-27

丁華（1984-），2009年畢業(yè)于清華大學(xué)熱能工程系流體機械研究所，碩士，工程師，2009年起進入東芝水電設(shè)備（杭州）有限公司工作，主要從事水輪機水力設(shè)計，模型開發(fā)及CFD解析研究?，F(xiàn)任水力機械研究所水力開發(fā)課主查。

中村高紀(jì)（1966-），男，1989年畢業(yè)于山口大學(xué)機械工程系,1991年在九州大學(xué)獲得熱能系統(tǒng)工程碩士學(xué)位，2012年在山口大學(xué)科學(xué)和技術(shù)學(xué)院獲得工科博士學(xué)位。1991年起在東芝公司工作，從事水輪機和水泵水輪機開發(fā)工作?，F(xiàn)任東芝能源和工業(yè)系統(tǒng)研究和開發(fā)中心旋轉(zhuǎn)機械部主干（首席工程師）。他是日本機械工程師學(xué)會（JSME）和日本旋轉(zhuǎn)機械學(xué)會（TSJ）成員。

德宮健男（1964-），男，1992年畢業(yè)于東京大學(xué)，株式會社東芝公司水輪機專家?，F(xiàn)任東芝水電設(shè)備(杭州)有限公司常務(wù)副總經(jīng)理。

丁華（1984-），男，上海人，2009年畢業(yè)于清華大學(xué)熱能工程系流體機械研究所，碩士研究生，工程師，主要從事水輪機水力設(shè)計，模型開發(fā)及CFD解析研究。杭州，東芝水電設(shè)備（杭州）有限公司水力機械研究所，310016。Email: ding.hua@toshiba-thpc.com

丁華郵箱：ding.hua@toshiba-thpc.com

聯(lián)系方式：0571-56809073

feng.minjuan@toshiba-thpc.com

審稿人：李任飛