多能互補條件下轉輪優化對水輪機低負荷區穩定性能的影響

趙亞萍，鄭小波，張 歡，郭鵬程

多能互補條件下轉輪優化對水輪機低負荷區穩定性能的影響

趙亞萍1，鄭小波1※，張 歡2，郭鵬程1

（1. 西安理工大學水利水電學院，西安 710048；2. 浙江富春江水電設備有限公司，杭州 311121）

多能互補系統中新能源發電的不穩定性使得作為調能機組的水電機組頻繁在水力效率低、振動劇烈的低負荷區運行，嚴重影響機組的壽命。該研究以多能互補系統中的混流式水輪機為研究對象，在前期考慮工況權重系數的轉輪多工況優化設計結果基礎上，對比分析了優化前后轉輪葉片的幾何參數變化，不同負荷區的水輪機內部流動狀態及壓力脈動特征差異。研究結果表明：優化后葉片包角、安放角以及葉片長度均有所增加，葉片表面壓力分布及轉輪進出水邊速度矩分布更加均勻，有助于改善水輪機低負荷區的空化性能、提高能量轉換能力。轉輪進出口安放角的增加很好地抑制了轉輪進口背面脫流渦及出水邊的脫流渦區，改善了尾水管的入流條件，使得尾水管渦帶的強度和影響范圍明顯減小。葉片優化后，轉輪內各頻率的壓力脈動幅值均有不同程度的降低，尾水管內壓力脈動改善明顯。尾水管內0.2f（f為轉頻）和14f壓力脈動在低負荷工況（OP1）幅值降幅分別為45%和40%，額定工況（OP4）尾水管內0.2f壓力脈動基本消除，14f壓力脈動幅值降幅為31%。本文所得研究結果對多能互補系統中水輪機轉輪的設計優化及運行具有參考意義。

水力發電；混流式水輪機；壓力脈動；多能互補；低負荷工況；渦帶

0 引 言

碳達峰碳中和背景下，能源結構的變化使得新能源發電飛速發展[1-2]，風、光、水多能互補已成為中國能源轉型的主流方向之一。由于風、光等新能源的波動性、間歇性和不穩定性，使得運行靈活、啟動迅速的水力發電在多能互補系統中常需要承擔調節負荷的任務[3]。

常規水電機組轉為多能互補發電系統中的調能機組后，水輪機被迫長時間在水流流動紊亂、水力效率低、振動劇烈的低負荷區運行，嚴重影響機組的安全運行[4-5]。所以，改善水輪機過流部件的優化設計方法、抑制水力不穩定現象、拓寬水輪機的高效穩定運行范圍成為多能互補系統中水電機組研究的迫切需求[6-7]。

張軍智等[8]針對多能互補下水電站的穩定運行進行了較早的研究，總結了李家峽水電站水輪機的轉輪改造關鍵點，提出了風光水多能互補條件下混流式水輪機轉輪全負荷運行下轉輪改造關鍵，為風光水多能互補下水電站水輪機轉輪設計提供了思路。ZHAO等[9-11]對水輪機穩定工況及過渡過程的振動信號和壓力脈動特性進行系統研究，揭示了誘發不穩定高振幅壓力脈動的主要原因，討論了在不穩定情況發生前檢測不穩定的可能性。劉靜[12]對風-水互補發電系統的供電質量和水力機組的運行穩定性進行了大量研究，分析了該系統中不同參數（如不同風水容量配比）對其頻率、電壓等方面的影響。馬騰宇等[13-16]針對混流式水輪機在水光互補系統中非設計工況下水輪機的穩定性進行了研究，發現轉輪進口水流角與安放角的不匹配所引起的主頻為導葉通過頻率的壓力脈動以及尾水管渦帶所引起的壓力脈動，是水輪機振動劇烈的主要因素，會引起葉片出水邊與上冠和下環連接部分的較大形變，甚至斷裂。由此可見，風、光、水多能互補系統中，考慮水電機組的調節功能時，水輪機內水力振動是無法忽視的問題。

對水輪機性能的常規改善方法主要是結合數值模擬和優化算法對轉輪進行優化設計[17]。MARCOS等[18-21]在實現轉輪葉片的批量建模的基礎上，建立了包含計算流體動力學、試驗設計法、響應模型和多目標遺傳算法的多目標優化設計系統，以葉片效率、空化、應力等參數為目標函數，有效的提高水輪機的效率和穩定性，改善了轉輪的空化及受力情況，利用可視化系統較好的展示了優化前后轉輪的形狀與葉片附近的速度流線和應變。ZBIGNIEW等[22]基于速度矢量理論，結合算流體動力學結果對高比轉速的混流式水輪機葉片進行逆向分析，調整流道和葉片形狀，經過模型試驗驗證了優化后葉片效率升高，表明了該優化方法的可行性。

目前，通過轉輪優化設計改善水輪機穩定性的研究對象主要以常規水力發電機組為主，而對多能互補系統中調能水輪機穩定性的研究主要集中在控制策略和系統評價上[23-27]。多能互補機組由于頻繁運行在低負荷工況，且需頻繁變換工況，機組運行的穩定差，因此對轉輪水力設計也提出了更高的要求，文獻[8]在總結李家峽水電站作為風光水多能互補系統中的調能機組進行改造的具體經驗時，指出不同工況權重因子對于水輪機的整體性能具有重要影響。可見多能互補系統中工況權重因子是水輪機水力設計過程中需著重考慮的重要因素之一。

風光水多能互補條件下，對作為調能機組運行的水輪機的基本要求是在保證額定負荷區運行性能的前提下，提高低負荷區甚至超低負荷區的水力效率，改善和抑制水輪機振動，拓寬水輪機的高效穩定運行范圍。前期的研究中以多能互補系統中混流式水輪機為例，建立了基于超傳遞近似法確定最優的運行工況權重系數的多能互補條件下的水輪機轉輪多工況多目標優化方法[28]。優化得到的轉輪葉片能夠在保證額定工況的能量特性的同時提高水輪機低負荷工況的運行性能。

振動劇烈是限制混流式水輪機低負荷運行的主要原因。基于此，本文在轉輪葉片優化的基礎上[28]，對優化前后的水輪機進行非定常數值計算，詳細對比優化前后轉輪葉片幾何形狀差異，研究優化前后水輪機轉輪內不穩定流動特征、渦結構的分布特征及產生機理，系統分析水輪機內不同頻率壓力脈動的產生機理，探討葉片幾何參數對水輪機內部流動的影響，為考慮運行工況權重因子的多工況多目標優化系統的可行性提供理論支持。

1 數值計算方法

1.1 計算域及工況點選取

本文以某多能互補條件下的混流式水輪機為例，選取水輪機整機為計算域（圖1），包括蝸殼、固定導葉、活動導葉、轉輪以及尾水管。水輪機轉輪直徑1=2.46 m，額定轉速=250 r/min，活動導葉和固定導葉數為24，轉輪葉片數為14。

圖1 計算域及網格劃分

由于多能互補系統中水輪機頻繁在低負荷區運行，因此當混流式水輪機轉為調能機組時，希望其具有較寬的運行范圍，尤其是低負荷區運行性能良好。因此本文根據水輪機模型特性曲線，考慮到多能互補條件下水輪機的運行范圍有可能會涉及全工況范圍，因此在同時兼顧額定工況和偏工況下水輪機的穩定運行，又能提高水輪機在低負荷區的運行性能的條件下，最終分別選取40%、60%、70%及100%出力的4個工況點，分別記為：OP1、OP2、OP3和OP4。其中OP3和OP4分別為最優工況和額定工況。各個工況的特征參數如表1所示。

表1 工況點參數

1.2 計算域網格劃分

本文采用六面體網格對計算域進行離散。在優化設計過程中涉及轉輪批量劃分網格的過程，利用Turbogrid軟件在旋轉機械領域能短時間畫出高質量網格的優勢，采用Turbogrid對參數化后的轉輪及優化設計中生成的樣本空間進行網格劃分。其余部件均采用ICEM CFD軟件進行網格劃分，為保證計算數據準確性，在劃分各流體域網格時，針對蝸殼隔舌處、導葉及轉輪葉片進出口、尾水管壁面等幾何突變位置進行了網格加密處理。為了減小網格數對數值計算結果的影響，盡量提高計算效率、縮短優化時長，本文以水頭為變量，選取額定工況點進行網格無關性驗證，如圖2所示。可見當整個流體域的網格總數為4.75×106時，水頭增長的趨勢逐漸趨于穩定。因此最終選取流體域網格數量為475萬，其中蝸殼網格數為20.5萬、固定導葉網格數為90.1萬、活動導葉網格數為126.7萬、轉輪網格數為162.5萬、尾水管網格數為72.4萬，各個部件的網格分布如圖1中的網格局部圖。為了更好的捕捉近壁區的流動特征，網格劃分過程對轉輪區域的近壁面網格進行控制，最終轉輪內除下環出口局部高湍流區外，大部分區域的Yplus小于10，滿足湍流模型。

圖2 網格無關性驗證

1.3 數值計算模型

水輪機內的流動是以水為介質的三維非穩態流動，一般認為水是不可壓縮流體，熱交換量很小，可以不考慮能量守恒，因此水輪機內復雜的三維黏性不可壓縮流動，可由連續性方程和動量方法來描述[29]：

1）連續性方程

2）動量方程

式中為流體的密度，kg/m3；τ為流體的切應力，N/m2；為壓強，Pa；為流體速度矢量，m/s；；x與x為空間坐標分量，m；S為附加源項，N/m3。

通過N-S方程描述湍流運動時具有不封閉性，因此引入湍流模型來封閉方程組，本文采用考慮了湍流剪切應力運輸、不會對渦流黏度造成過度預測[30]的-SST湍流模型來求解水輪機內部復雜的三維不可壓縮流動。

數值計算過程中，以蝸殼進口作為計算域的進口，尾水管出口為計算域出口。具體的邊界條件為：進口給定質量流量，并假設速度方向垂直于蝸殼進口面；出口給定相對壓力；假設固壁面無滑移，靠近壁面區域采用標準壁面函數法進行處理；導葉域與轉輪域以及轉輪域與尾水管域之間的靜止和轉動部分采用凍結轉子法。

為了驗證數值模擬結果的可靠性，本文以水輪機效率為評價值，將數值模擬結果與原模型試驗結果進行對比，如圖3所示。可見，數值模擬值與試驗值具有很好的一致性，表明采用該計算模型預估水輪機轉輪在不同運行工況下的水力性能是可行的。

圖3 試驗及數值模擬結果對比

1.4 壓力測點布置

本文重點對轉輪與尾水管內的水力不穩定因素及其特征進行分析，并在轉輪葉片壓力面、吸力面及尾水管內布置壓力脈動測點。圖4a、圖4b為數值計算過程中轉輪葉片壓力測點布置，測點命名規律如下：葉片壓力面記為RNPS，在靠近葉片上冠附近，沿流線方向分別布置3個監測點。其中RN表示轉輪，PS表示壓力面，PS后第一位數字表示處于葉片上冠附近，最后一位數字表示監測點編號。在葉片中部布置監測點（RNPS21、RNPS22、RNPS23）、葉片下環附近布置監測點（RNPS31、RNPS32、RNPS33）。葉片吸力面記為RNSS，吸力面上的監測點布置情況與壓力面類似，轉輪葉片兩側監測點總計18個。尾水管監測點布置如圖4c、圖4d所示，從尾水管錐管段至彎肘段共選取4個監控截面，每個截面布置4個壓力監控點，監控點布置規則為：每個沿截面圓周順時針方向等間距對稱布置4個監控點，編號為DT01至DT04。其中DT表示尾水管，DT后第一位數字表示該點所在平面，后兩位數字表示監測點編號，尾水管中監測點總計16個。

圖4 水輪機內監測位置示意圖

為探究轉輪優化前后葉片幾何參數的變化對水輪機內部水力振動特征的影響，對各個監測點所得的壓力脈動統一壓力系數C進行無量綱化如式（3）。通過快速傅里葉變換獲得水輪機內不同部位的壓力脈動頻譜特性。

2 優化前后葉片幾何參數對比

基于文獻[28]建立的多能互補系統中水輪機轉輪的優化設計方法得到的優化前后的轉輪葉片。圖5為水輪機優化前后的幾何形狀及葉片表面均勻分布的5個三維截面對比，其中Span0為葉片上冠處截面；Span1.00為葉片下環處截面；Span0.25、Span0.50、Span0.75為葉片中間均勻布置的3個截面。為了更為詳細的描述優化前后葉片幾何差異，圖6分別給出了優化前后葉片不同截面的安放角，包角及厚度分布。

從圖6a可以看出，葉片包角從頭部到尾部均有增加，且越靠近葉片出水邊（流向=1），包角的增加越明顯，尤其葉片下環處翼型包角由45°增加至51°，增幅約為13%，使得葉片長度增加10%，葉片的表面積也相應增加。因此單位面積葉片上承受的壓差減小，有助于改善葉片的空化性能。圖6b中葉片優化前后安放角差異最大的位置位于葉片靠上冠處，其余截面的安放角從頭部到尾部逐漸增加。葉片進口靠上冠處安放角較原始葉片增大約7%，葉片進口靠下環處安放角較原始葉片增大約5%，靠近葉片頭部10%～30%處安放角增加并不明顯。葉片出口從上冠至下環，葉片安放角增加明顯，增加幅度約為12%。葉片進出口安放角的增加必然會對尾水管內的流動狀態產生影響。圖6c中，優化前后葉片厚度差異主要集中在Span 0.75和Span 1.00處，其余截面厚度變化較小，究其原因主要是由于葉片包角的增加使得葉片在軸面上的投影厚度增加。

注：Span表示葉高位置。

圖6 優化前后葉片幾何參數對比

3 結果與分析

3.1 水輪機內部流動特性分析

圖7為工況OP1和OP4轉輪葉片表面壓力分布，從圖可以看出，優化后葉片包角增大長度增加使得葉片表面低壓區面積減小，壓力分布更加均勻。在OP1工況40%出力時，優化后葉片正面出水邊近上冠處以及葉片背面近下環處的低壓區明顯減小。葉片進口處壓力增加了35.39 kPa，葉片尾部最低壓力提高了34.51 kPa。在100%出力工況時，優化后葉片進口處壓力提升較為明顯，葉片從頭部到尾部的壓力梯度變化較大，葉片出水邊低壓區沿葉展方向的分布明顯減小。由圖8的轉輪軸面速度矩分布可以看出，額定工況（OP4），優化前后轉輪軸截面上的速度矩分布差別較小，而低負荷工況（OP1），優化后的轉輪進出水邊速度矩分布更加均勻，因此葉片沿葉展方向的做功能力更加均衡，有助于改善低負荷區葉片表面的受力狀態。可見，優化后的轉輪葉片在保證100%出力工況水輪機性能的情況下，能夠很好的改善低負荷區葉片頭部撞擊空化以及葉片背面翼型空化，從而降低水力損失，提高水輪機水力效率。

圖7 葉片表面壓力分布

圖8 轉輪軸面速度矩分布

水輪機內部的復雜不穩定渦是造成低負荷區運行時性能惡劣的主要原因，圖9給出了40%出力工況時（OP1）轉輪內部的渦分布。從圖中可以看出，優化前轉輪內的渦主要由葉道渦（渦1）、葉片靠上冠處過大正沖角引起的葉片背面脫流渦（渦2）、葉片靠下環處頭部過大負沖角引起的葉片正面脫流渦（渦3）、以及葉片背面出水邊的脫流渦區（渦4）組成。葉片靠上冠處，水流到達葉片頭部位置時以較大的正沖擊擾流葉片在葉片頭部背面形成脫流，當脫流遠離葉片背面時便產生回流，形成渦1；受渦1影響，使得葉道內水流產生由葉片正面向相鄰葉片背面運動的趨勢，并與轉輪的旋轉方向恰恰相反，因此形成渦2。由于優化前葉片中部靠近下環處為負沖角，過大的負沖角使得葉片頭部正面產生脫流，隨著轉輪的旋轉，當葉片背面的水流運動至葉片中后方時產生大面積的脫流區，形成渦4；而靠近下環處，受下環壁面的約束水流無法充分發展，因此在下環處形成從葉片頭部發展至尾部的旋轉渦，形成渦3。通過對轉輪葉片進行優化，葉片靠上冠處從頭部到尾部安放角增加改善了此處水流的繞流狀態，有效的減小了葉道渦（渦1），破壞了渦2的形成條件；同時在葉片包角和安放角的增大使得葉片長度增加，延長了水流繞流葉片的路徑，使得水流能夠得到充分發展，有效的抑制了葉片背面出水邊的脫流渦4。轉輪內渦流的改善將有助于于降低水力損失和渦流引起的壓力脈動。

圖10為40%出力工況下水輪機內部湍動能分布，從圖可以看出，轉輪內的高湍動能區域也主要分布在產生渦的位置。葉片靠上冠處高湍動能區主要位于葉片背面，其位置與渦1、渦2和渦4一致，如圖10c；葉片靠下環處的高湍動能區主要存在于葉片頭部正面及葉片背面出水邊，其位置與渦3和渦4一致。通過對轉輪葉片進行優化，葉片進口安放角增加，有效的減小了葉片頭部的正撞擊，抑制了葉片頭部渦2的形成，使得渦1的強度得以削弱；葉片出水邊安放角減加，改善了葉片背面的翼型脫流，很好的控制了葉片背面靠出水邊的脫流渦，使得此處的湍動能減小。

轉輪幾何參數的變化不僅會引起轉輪內部流動特性的變化，同時會對尾水管內流態產生影響。圖11為轉輪優化前后各個工況下尾水管渦帶及圖4c所示截面的壓力分布情況。可以看出，在OP1和OP2工況下，轉輪優化前，尾水管錐管段形成一段明顯螺旋狀偏心渦帶，轉輪優化后的尾水管渦帶強度明顯減小，渦帶偏心度也較小，各個截面的低壓區面積減小，壓力分布趨于均勻。OP3工況為最優工況點，水輪機整體流態較好，在轉輪出口與尾水管交界面中心位置形成柱狀渦帶，且低壓區一直處于尾水管中心位置，轉輪優化后，由于尾水管最低壓力值升高，尾水管中心柱狀渦的體積明顯減小，在截面2處柱狀渦帶已經消失。在OP4工況下，尾水管內形成體積較大中心柱狀渦，尾水管中心壓力值更低，壓力分布具有明顯的對稱性，轉輪優化后的尾水管內的柱狀渦帶消失，形成半徑增加，但長度顯著減小的柱狀渦帶，很好的減小了尾水管渦帶的影響范圍。

3.2 水輪機壓力脈動特征分析

3.2.1 轉輪內壓力脈動分析

圖12和圖13分別為OP1工況和OP4工況轉輪優化前后葉片兩側各個監測點的時域與頻域圖。從圖可以看出，優化前后葉片上各個監測點的壓力脈動系數隨時間變化規律相似，葉片進口以及葉片中部處的壓力脈動系數隨時間變化的趨勢較出口處比較平緩，出口處的壓力脈動系數幅值相對較高。轉輪內頻率為轉頻f的壓力脈動幅值很小。此外，葉片進口邊靠下環位置（監測點RNPS31、RNSS31）與導葉的距離最近，轉輪旋轉過程動靜干涉現象在該位置最為突出，因此出現了與導葉數相關的24f頻率的壓力脈動。葉片優化前，轉輪內以0.20f和0.80f頻率的壓力脈動為主，葉片靠近上冠及中間位置的壓力脈動幅值沿進口至出口逐漸增加，葉片下環處壓力脈動幅值變化均勻。優化后轉輪內壓力脈動幅值下降明顯，尤其是0.20f的低頻壓力脈動基本消失。0.80f頻率的壓力脈動在轉輪優化前后均存在，且吸力面的幅值明顯大于壓力面，壓力脈動幅值最大的位置分別為位于葉片出水邊的RNSS/PS13、RNSS/PS23這4個點，可見，壓力脈動幅值最大位置存在于葉片出水邊。轉輪優化后靠葉片進水邊及中間位置的幅值下降明顯，該現象與葉片背面脫流渦區4分布一致，而在葉片出水邊靠下環處（RNSS/PS33處），優化后的壓力脈動幅值有微小的上升，主要是由于優化后渦4的體積明顯減小，且緊貼葉片與下環形成的三角區流出葉道，導致下環位置流動紊亂，并引起壓力脈動幅值的增加，可見該頻率的壓力脈動主要由于葉片出水邊脫流渦4引起。轉輪進口脫流渦（渦2、渦3）以及葉道渦（渦1）雖然并未產生明顯的壓力脈動，因此對低負荷情況下轉輪的穩定性差的貢獻較小。從OP4工況（圖12）葉片表面壓力脈動可知，優化前后葉片表面的壓力脈動差異較小，主要以轉頻f為主，在葉片進口邊靠下環處由于與導葉的動靜干涉作用存在主頻為24f的壓力脈動。可見，優化后的轉輪能夠在保證額定工況的性能的情況下使得很好的改善低負荷區的壓力脈動特征，使得各頻率的壓力脈動幅值整體降低。

圖10 轉輪內部湍動能分布

圖11 尾水管渦帶及截面壓力分布

注：fn為轉輪的旋轉頻率，Hz；Cp為壓力系數。

圖13 OP4工況轉輪壓力脈動時域圖與頻域圖

3.2.2 尾水管內壓力脈動分析

圖14和圖15為OP1和OP4兩個工況下尾水管4個平面內各個監測點的壓力脈動頻譜特性。尾水管上各個平面以主頻為0.20f的壓力脈動為主，且其幅值遠大于轉輪內存在的0.20f頻率的壓力脈動，并且在離轉輪較近位置（截面1）存在14f（葉片通過頻率）的壓力脈動。

圖14 OP1工況尾水管壓力脈動時域圖與頻域圖

圖15 OP4工況尾水管壓力脈動時域圖與頻域圖（100%出力）

從圖14可以看出，在OP1工況，14f的壓力脈動幅值較低，對尾水管內部整體流態影響不大。在截面1、2上0.20f的壓力脈動最大值出現在監測點3、4位置，隨著流態的發展，受逆時針轉動的尾水管渦帶的影響，壓力脈動的最大位置沿逆時針方向轉動轉移至截面3、4的測點2、3位置。轉輪優化后，尾水管內的壓力脈動仍然為0.20f和14f兩個頻率，但壓力脈動幅值有顯著的降低，其中一階主頻（0.2f）壓力脈動最大幅值降低約45%，二階主頻（14f）壓力脈動最大幅值降低約40%。從圖15可以看出，在額定工況下，尾水管錐管段的壓力脈動較小，水流進入彎肘段（截面4），過流通道的劇烈變化引起水流流動狀態的變異，使得0.2f的低頻壓力脈動幅值增加。同時位于截面1的14f頻率壓力脈動幅值與OP1的差異較小，并未向下游傳遞。轉輪優化前，尾水管內0.2f的低頻壓力脈動的壓力脈動最大值的壓力系數為0.121、而優化后幅值僅為0.0015，同時14f頻率壓力脈動幅值也由0.0055降為0.0038，降幅約31%。可見，轉輪的通過頻率對水輪機尾水管的穩定性影響較小，通過對轉輪進行優化，能夠很好地改善不同負荷區尾水管內的低頻壓力脈動，進而提高水輪機的穩定性。

由上面分析可知，轉輪和尾水管內均存在0.2f的低頻壓力脈動，且該頻率的壓力脈動在尾水管內幅值遠大于轉輪內，因此推測該壓力脈動來源于尾水管錐管段，與尾水管渦帶有關。為了進一步證明該推測，截取如圖16a所示的監測點DT105的時域信息圖（圖16b）的各個時刻尾水管中渦帶的形態（圖16c）進行對比分析。從尾水管渦帶形態可以看出，時刻A與C尾水管渦帶形態相似，相位相差360°，時刻B與D的尾水管渦帶形態似度極高。時刻在A，尾水管渦帶距監測點較遠，壓力脈動系數為最小值。時刻B，螺旋狀的尾水管渦轉至監測點附近，此時對應壓力脈動的波峰。可見在從時刻B到D（=3.22～4.51 s）尾水管渦帶運動一個周期，頻率約為0.2f，其后尾水管渦帶進行重復性運動，引發錐管段產生周期性壓力脈動。由此可見，尾水管中0.20f是由尾水管渦帶引起。

注：A、B、C、D為波峰或波谷時刻。

通過對轉輪葉片優化前后的水輪機性能進行分析可知，優化后的轉輪能夠很好的改善水輪機在低負荷區的空化性能及穩定性，更適用于風光水多能互補系統。

4 結 論

本文以多能互補運行條件下的混流式水輪機轉輪為研究對象，對其進行考慮工況權重系數的轉輪多工況優化設計，獲得適用于多能互補條件下的水輪機轉輪。并對優化前后水輪機不同負荷區的運行性能進行詳細分析，得出以下結論：

1）基于前期的優化結果，分析優化前后葉片幾何參數可知，葉片包角及安放角從頭部到尾部均逐漸增加，優化后葉片包角、安放角以及葉片長度的最大增幅分別為13%、7%和10%，且優化后葉片表面壓力分布及轉輪進出水邊速度矩分布更加均勻，有助于很好的改善低負荷區空化特性及能量轉換能力。

2）優化后轉輪進出口安放角的增加很好的抑制了轉輪進口背面脫流渦及出水邊的脫流渦區，改善了尾水管的入流條件，使得低負荷區尾水管最低壓力值升高，尾水管渦帶強度和影響范圍明顯減小。

3）低負荷區域，引起水輪機內明顯壓力脈動的因素有靠下環處的轉輪與導葉的動靜干涉引起的頻率為24f壓力脈動、葉片背面靠出水邊的不穩定脫流渦引起的頻率為0.8f低頻壓力脈動、轉輪出口動靜干涉引起的頻率為14f壓力脈動以及尾水管渦帶引起的頻率為0.2f低頻壓力脈動，且尾水管渦帶引起的壓力脈動會向上游傳遞，使得轉輪內的壓力脈動狀況更加復雜。通過對轉輪進行優化，雖并未消除某個頻率的壓力脈動，但使得壓力脈動的幅值明顯減小，尤其尾水管內壓力脈動改善明顯。低負荷工況（OP1）尾水管內0.2f和14f壓力脈動在幅值降幅可達45%和40%，額定工況（OP4）尾水管內0.2f壓力脈動基本消除，14f壓力脈動幅值降幅為31%，很好的改善水輪機在低負荷區的運行穩定性。

[1] 莊貴陽. 我國實現“雙碳”目標面臨的挑戰及對策[J]. 人民論壇，2021(18)：50-53.

[2] 曹軍文，鄭云，張文強，等. 能源互聯網推動下的氫能發展[J]. 清華大學學報（自然科學版），2021，61(4)：302-311.

CAO Junwen, ZHENG Yun, ZHANG Wenqiang, et al. Hydrogen energy development driven by the energy internet[J]. Journal of Tsinghua University (Science and Technology), 2021, 61(4): 302-311. (in Chinese with English abstract)

[3] 羅興锜，朱國俊，馮建軍. 水輪機技術進展與發展趨勢[J]. 水力發電學報，2020，39(8)：1-18.

LUO Xingqi, ZHU Guojun, FENG Jianjun. Progress and development trends in hydraulic turbine technology[J]. Journal of Hydroelectric Engineering, 2020, 39(8): 1-18. (in Chinese with English abstract)

[4] 劉攀，陳學力，汪泉，等. 高水頭混流式水輪機的動靜干涉與振動問題研究[J]. 水力發電學報，2016，35(3)：91-98.

LIU Pan, CHEN Xueli, WANG Quan, et al. Analysis of rotor-stator interaction and vibration in high-head Francis turbines[J]. Journal of Hydroelectric Engineering, 2016, 35(3): 91-98. (in Chinese with English abstract)

[5] 馮金海，凌祖光，趙子文，等. 偏負荷工況混流式水輪機轉輪結構穩定性研究[J]. 水力發電學報，2021，40(1)：107-114.

FENG Jinhai, LING Zuguang, ZHAO Ziwen, et al. Study on structural stability of Francis turbine impellers under partial load condition[J]. Journal of Hydroelectric Engineering, 2021, 40(1): 107-114. (in Chinese with English abstract)

[6] YU G Y. Analysis of the development status and prospect of multi-energy complementary technology[C]. Jinan, 6th International Conference on Advances in Energy, Environment and Chemical Engineering (AEECE 2020),2020

[7] 鮑海艷，龍麗婷，付亮，等. 水輪機調速器功率調節模式下負荷調節過渡過程穩定性研究[J]. 農業工程學報，2019，35(17)：50-57.

Bao Haiyan, Long Liting, Fu Liang, et al. Stability of load regulation transition process of hydro turbine governor in powerregulation mode[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2019, 35(17): 50-57. (in Chinese with English abstract)

[8] 張軍智，何振忠，王子瑞，等. 風光水多能互補混流式水輪機全負荷范圍穩定經濟運行設計[J].西北水電，2021(6)：15-21.

ZHANG Junzhi, HE Zhenzhong, WANG Zirui, et al. Design for stable and economical operation of francis turbine of Full-Load range in Wind-PV-Hydro complementary power system[J]. Northwest Hydropower, 2021(6): 15-21. (in Chinese with English abstract)

[9] ZHAO W Q, PRESAS A, EGUSQUIZA M, et al. Increasing the operating range and energy production in Francis turbines by an early detection of the overload instability[J]. Measurement, 2021, 181: 109580.

[10] 陳會向，劉漢中，王胤淞，等. 抽水蓄能機組低水頭起動過渡過程壓力脈動分析[J]. 農業工程學報，2023，39(6)：63-72.

CHEN Huixiang, LIU Hanzhong, WANG Yinsong, et al. Pressure pulsation during low head start-up transient in a pumpedstorage hydropower unit[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2023, 39(6): 63-72. (in Chinese with English abstract)

[11] 孫龍剛，郭鵬程，鄭小波，等. 混流式水輪機葉道空化渦誘發高振幅壓力脈動特性[J]. 農業工程學報，2021，37(21)：62-70.

SUN Longgang, GUO Pengcheng, ZHENG Xiaobo, et al. Characteristics of high-amplitude pressure fluctuation induced byinter-blade cavitation vortex in Francis turbine[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions ofthe CSAE), 2021, 37(21): 62-70. (in Chinese with English abstract)

[12] 劉靜. 風水互補發電系統的電能質量與穩定性研究[D]. 楊凌：西北農林科技大學，2020.

LIU Jing. Research on the Power Quality and Stability of Wind-Hydro Complementary Power System[D]. Yangling: Northwest A&F University, 2020. (in Chinese with English abstract)

[13] 馬騰宇，唐雯，徐連琛，等. 水光互補系統中混流式水輪機轉輪結構穩定性研究[J]. 水利水電快報，2022，43(6)：78-84.

MA Tengyu, TANG Wen, XU Lianchen, et al. Research on structural stability of runner of Francis turbinebased on water-light complementarity[J]. Express Water Resources & Hydropower Information, 2022, 43(6): 78-84. (in Chinese with English abstract)

[14] ZUO Z G, LIU S H, LIU D M, et al. Numerical predictions and stability analysis of cavitating draft tube vortices at high head in a model Francis turbine[J]. Science China Technological Sciences, 2014, 57(11): 2106-2114.

[15] 龐嘉揚，劉小兵，宋罕，等. 基于水光蓄互補聯合發電系統混流式水輪機在超低出力區工況穩定性研究[J]. 水電與抽水蓄能，2021，7(5)：39-46.

PANG Jiayang, LIU Xiaobing, SONG Han. Study on stability of francis turbine in ultra-low output region based on hydro-photovoltaic-storage combined power generation system[J]. Hydropower and Pumped Storage, 2021, 7(5): 39-46. (in Chinese with English abstract)

[16] 楊鈞翔，劉小兵，龐嘉揚，等. 多能互補發電系統中楊家灣水電站水輪機壓力脈動特性的數值研究[J]. 水電與抽水蓄能，2021，7(5)：59-63，110.

YANG Junxiang, LIU Xiaobing, PANG Jiayang, et al. Numerical study on pressure fluctuation characteristics ofturbine in yangjiawan hydropower station in multi-energy complementary power generation system[J]. Hydropower and Pumped Storage, 2021, 7(5): 59-63, 110. (in Chinese with English abstract)

[17] ZHOU Y K, ZHENG S Q. Multi-level uncertainty optimisation on phase change materials integrated renewable systems with hybrid ventilations and active cooling[J]. Energy, 2020, 202: 117747.

[18] TSUZUKI M, VATANABE S L, CASTRO E G, et al. Development of a complete methodology to reconstruct, optimize, analyze and visualize francis turbine runners[C]// Ottawa, 15th IFAC Symposium on Information Control Problems in Manufacturing (INCOM 2015), 2015.

[19] ZHU B S, WANG X H, TAN L, et al. Optimization design of a reversible pump-turbine runner with high efficiency and stability[J]. Renewable Energy, 2015, 81: 366-376.

[20] LEONEL A T, FRANCISCO J L, Sara A R. Performance improvement of a 500 kW Francis turbine based on CFD[J]. Renewable Energy, 2016, 96: 977-992.

[21] LIU L H, ZHU B S, BAI L, et al. Parametric design of an ultrahigh-head pump-turbine runner based on multiobjective optimization[J]. Energies, 2017, 10: 01169.

[22] ZBIGNIEW K, JANUSZ S. High specific speed Francis turbine for small hydro purposes- Design methodology based on solving the inverse problem in fluid mechanics and the cavitation test experience[J]. Renewable Energy, 2021, 169: 1210-1228.

[23] XIONG H L, EGUSQUIZA M, ALBERG ? P, et al. Multi-objective optimization of a hydro-wind-photovoltaic power complementary plant with a vibration avoidance strategy[J]. Applied Energy, 2021, 301: 117459.

[24] SEYED M M, MOHAMMAD R S, HAMID M. A Pareto optimal multi-objective optimization for a horizontal axis wind turbine blade airfoil sections utilizing exergy analysis and neural networks[J]. Journal of Wind Engineering and Industrial Aerodynamics, 2015, 136: 62-72.

[25] 李劍鋒，郝曉光，曾四鳴，等. 多能互補的綜合能源微網多目標優化運行方法研究[J/OL]. 華北電力大學學報（自然科學版）：1-10[2022-09-09]. http://kns.cnki.net/kcms/detail/ 13.1212.TM.20211122.1515.002.html

LI Jianfeng, HAO Xiaoguang, ZENG Siming, et al. Research on multi-energy complementary integrated energy micro-grid multi-objective optimal operation method[J/OL]. Journal of North China Electric Power University (Natural Science Edition): 1-10[2022-09-09]. http://kns.cnki.net/kcms/detail/ 13.1212.TM.20211122.1515.002.html

[26] 趙鳳展，李奇，張啟承，等. 多能互補能源系統多維度綜合評價方法[J]. 農業工程學報，2021，37(17)：204-210.

ZHAO Fengzhan, LI Qi, ZHANG Qicheng, et al. Multi-dimensional comprehensive evaluation of multi-energy complementary energy system[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2021, 37(17): 204-210. (in Chinese with English abstract)

[27] 王玨，廖溢文，韓文福，等. 碳達峰背景下抽水蓄能-風電聯合系統建模及有功功率控制特性研究[J]. 水利水電技術，2021，52(9)：172-181.

WANG Jue, LIAO Yiwen, HAN Wenfu, et al. Modeling and active power control characteristics of pumped storage-wind hybrid power system in the context of peak carbon dioxide emission[J]. Water Resources and Hydropower Engineering, 2021, 52(9):172-181. (in Chinese with English abstract)

[28] Zheng X B, Zhao Y P, Zhang H, et al. Optimization and performance analysis of Francis turbine runner based on super-transfer approximate method under multi-energy complementary conditions[J]. Sustainability, 2022, 14(16): 10331

[29] 陶文銓. 數值傳熱學[M]. 西安：西安交通大學出版社，2001.

[30] 周宇，錢煒祺，鄧有奇，等.-SST兩方程湍流模型中參數影響的初步分析[J]. 空氣動力學學報，2010，28(2)：213-217.

ZHOU Yu, QIAN Weiqi, DENG Youqi, et al. Introductory analysis of the influence of Menter’s-SST turbulence model’s parameters[J]. Acta Aerodynamica Sinica, 2010, 28(2): 213-217. (in Chinese with English abstract)

Influences of runner optimization on the stability performance of hydraulic turbine in the low-load range under the condition of multi-energy complementary

ZHAO Yaping1, ZHENG Xiaobo1※, ZHANG Huan2, GUO Pengcheng1

(1.-,,710048,; 2.,.,311121,)

Hydropower is often required to adjust the load in the multi-energy complementary system, due to the strong volatility, intermittency, and instability of new energy power generation. Therefore, the turbine is forced to operate in a low-load area with the low efficiency and severe vibration over a long time. The operating conditions vary frequently to threaten the stability and operating life of the unit. It is necessary to optimize the overcurrent components of the turbine for the hydroelectric unit in the multi-energy complementary system. The hydraulically unstable flow can be suppressed to broaden the high-efficiency operation range of the turbine. In this study, a multi-operating optimization of the runner was implemented to consider the weight coefficient in the operation of the turbine under the condition of multi-energy complementation of wind, solar and water. The runner of turbine was obtained suitable for the multi-energy complementary condition. The unsteady numerical analysis and comparison were also carried out on the turbines before and after optimization. The research results show that: An appropriate increase in the weight coefficient of the operating conditions in the low-load area was effectively improved the cavitation performance and the efficiency of the turbine after the multi-condition optimization, particularly with the operational performance of the turbine in the high-load area. There was an increase in the inlet and outlet placement angles of the optimized runner blades, leading to effectively reduce the attack angle of the heading edge of the blade and the flow angle of the tailing edge of the blade under low-load conditions. The vortex was better restrained from the heading edge to the back of the blade. There was an increase in the flow separation area at the tailing edge of the blade, and the inflow conditions of the draft tube. The strength of the vortex band in the draft tube was significantly reduced. Under the low load conditions, the pressure pulsation in the draft tube was mainly the 0.2flow-frequency pressure pulsation that caused by the vortex, and the low-amplitude pressure pulsation with the blade passing frequency at the inlet of the draft tube that caused by the rotation of the runner. The pressure pulsation in the runner was mainly the 0.8fn pressure pulsation that caused by the flow separation area at the tailing edge of the blade, while the 24fhigh-frequency pressure pulsation was caused by the dynamic and static interference between the guide vane and the runner near the band. The 0.2flow-frequency pressure pulsation was passed up from the draft tube. The pressure pulsation amplitudes of different frequencies in the runner and the draft tube were effectively reduced after optimization of the runner blade, indicating the particularly outstanding improvement of pressure pulsation in the draft tube. The amplitudes of pressure pulsation were reduced by 45% and 40%, respectively, in the draft tube with the frequency of 0.2fand 14funder the low load condition (OP1). There was no pressure pulsation with the frequency of 0.2fin the draft tube under rated condition (OP4). The amplitude of pressure pulsation of 14fn was reduced by 31%. The operation stability of the turbine was better improved in the low load area. The finding can provide a strong reference to optimize the operation of the turbine runner in the multi-energy complementary system.

hydro power; Francis turbine; pressure pulsation; multi-energy complementation; low load condition; vortex

2022-09-10

2023-01-13

國家自然科學基金項目（52009105；51839010）

趙亞萍，博士，講師，研究方向為水力機械優化設計及流動分析。Email：zyp0168@xaut.edu.cn.

鄭小波，博士，副教授，研究方向為流體機械振動與穩定性分析。Email：zxb@xaut.edu.cn

10.11975/j.issn.1002-6819.202209081

S277.9+4; TK733+.3

A

1002-6819(2023)-07-0067-10

趙亞萍，鄭小波，張歡，等. 多能互補條件下轉輪優化對水輪機低負荷區穩定性能的影響[J]. 農業工程學報，2023，39(7)：67-76. doi：10.11975/j.issn.1002-6819.202209081 http://www.tcsae.org

ZHAO Yaping, ZHENG Xiaobo, ZHANG Huan, et al. Influences of runner optimization on the stability performance of hydraulic turbine in the low-load range under the condition of multi-energy complementary[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2023, 39(7): 67-76. (in Chinese with English abstract) doi：10.11975/j.issn.1002-6819.202209081 http://www.tcsae.org