陳 上，陳金保，王安林，丘濤基，陳齊燈，肖志懷

（1.武漢大學動力與機械學院，武漢430072；2.福建華電金湖電力有限公司，福建三明353399）

0 引 言

隨著大規(guī)模風能、太陽能等新能源電能的接入，輸電電網(wǎng)架構和區(qū)域電網(wǎng)間互聯(lián)方式日趨復雜，在電力系統(tǒng)的規(guī)劃、調(diào)峰、安全備用、頻率電壓穩(wěn)定等方面帶來更多新的問題。作為電網(wǎng)中起主要調(diào)節(jié)作用的水電站，對其調(diào)節(jié)性能提出了更高要求。由于水輪機內(nèi)部流體運動復雜，各運行變量間的動態(tài)關系復雜、耦合程度高，難以用簡單的數(shù)學模型描述，目前只能通過真機試驗、模型試驗和近似解析方法來建立水輪機模型。在工程上通常采用穩(wěn)態(tài)特性近似代替動態(tài)特性建立水輪機模型［1］，進而對水輪機的穩(wěn)態(tài)特性、動態(tài)特性、控制策略等進行研究。水輪機建模主要包括外特性法［2-6］和內(nèi)特性法［7-9］，內(nèi)特性建模由于過度依賴水輪機幾何參數(shù)和結構參數(shù)，且實際建模時很難獲得具體準確的參數(shù)，故很少采用。外特性建模主要基于模型運轉(zhuǎn)綜合特性曲線（或原型運轉(zhuǎn)綜合特性曲線），位于模型綜合特性曲線內(nèi)的工況通?？奢^好表達，但小開度區(qū)間建模還需要進一步深入研究；此外，也有學者提出其他水輪機建模方法，如適用于電力系統(tǒng)仿真的遞推水輪機線性模型［10］；反映系統(tǒng)能量特性的非線性水輪機哈密爾頓模型［11］。為了分析各類模型之間的建模精度與差異，為水輪機建模提供理論指導，本文結合目前常用的分段線性化模型、理想水輪機模型、神經(jīng)網(wǎng)絡模型和插值法模型，借助matlab 仿真平臺，針對實際機組開展仿真建模研究，對比分析不同建模方法在穩(wěn)態(tài)特性、動態(tài)特性、建模復雜程度、模型仿真時間等方面效果與差異。

1 水輪機建模原理

1.1 分段線性化模型

水輪機動態(tài)特性通常使用水輪機穩(wěn)態(tài)工況下的力矩與流量特性表示，其非線性穩(wěn)態(tài)表達式為式（1）和式（2）：

式中：Mt為水輪機力矩；a為導葉開度；H為水輪機水頭；z為槳葉角度；n為機組轉(zhuǎn)速；Q為水輪機流量。

將這兩個式子展開為泰勒級數(shù)并略去二次及以上高次項得式（3）與（4）。

為方便描述，將其偏導數(shù)替換成傳遞系數(shù)，如式（5）：

該方程描述穩(wěn)態(tài)工況附近微小偏差時力矩和流量的變化。即該方程的應用條件是已知該工況的傳遞系數(shù)，且波動幅度在工程允許范圍之內(nèi)。當工況波動過大時，該方程計算數(shù)據(jù)會與實際結果產(chǎn)生較大的誤差，故當波動較大時，此方法應用受到限制。

分段線性化建模中，關鍵是工況的劃分和傳遞系數(shù)的計算。工況劃分通常以水頭和導葉開度等比例選取，工況劃分越細，模型精度越高。

傳遞系數(shù)計算通常采用兩點法或者曲線擬合法，適合計算有限工況點的傳遞系數(shù)，但計算繁瑣，難以清晰全面認識傳遞系數(shù)隨工況變化規(guī)律，故而出現(xiàn)通過模型綜合特性曲線結合BP神經(jīng)網(wǎng)絡對數(shù)據(jù)處理后直接進行求偏導數(shù)的處理方法，此方法計算既保證了計算精度又減少了計算量［12］。

計算傳遞系數(shù)，得出如圖1 所示的水輪機分段線性化模型方框圖［13］，圖1中Tw為水流慣性時間常數(shù)。

1.2 理想水輪機模型

主要分析調(diào)節(jié)系統(tǒng)其他環(huán)節(jié)時，一般采用理想水輪機模型表示水輪機特性。理想水輪機模型是對1.1 節(jié)分段線性化模型進行簡化，應用時不考慮槳葉角度，通過槳葉角度與導葉開度協(xié)聯(lián)，將轉(zhuǎn)漿式水輪機視為混流式水輪機，考慮到轉(zhuǎn)速對力矩和流量影響較小，理想水輪機模型中傳遞系數(shù)ex和eqx取為0，在額定工況附近，取ey=1.0，eh=1.5，eqy=1.0，eqh=0.5，將傳遞系數(shù)合并得到傳遞函數(shù)如圖2所示。

圖2 理想水輪機模型Fig.2 Ideal hydraulic turbine model

1.3 神經(jīng)網(wǎng)絡模型

隨著人工智能發(fā)展，基于神經(jīng)網(wǎng)絡技術的非線性建模技術得到廣泛應用。針對BP 神經(jīng)網(wǎng)絡，只要有足夠多隱含層和隱含層節(jié)點，就可以逼近任意非線性映射關系，許多學者將神經(jīng)網(wǎng)絡引入水輪機建模之中［14-17］。

實時仿真中，利用神經(jīng)網(wǎng)絡實現(xiàn)水輪機非線性特性計算具有計算簡單、精度可控且導數(shù)連續(xù)等優(yōu)點?？晒┻x擇的神經(jīng)網(wǎng)絡很多，如前饋神經(jīng)網(wǎng)絡、循環(huán)網(wǎng)絡和對稱連接網(wǎng)絡等，水輪機建模中使用前饋神經(jīng)網(wǎng)絡中BP 神經(jīng)網(wǎng)絡較多，其網(wǎng)絡結構如圖3所示。

圖3 BP神經(jīng)網(wǎng)絡結構圖Fig.3 BP neural network structure diagram

該網(wǎng)絡包含三輸入、單輸出和一個隱層，實現(xiàn)軸流轉(zhuǎn)漿式機組流量或力矩模擬，其隱層神經(jīng)元個數(shù)可調(diào)整，輸出層為線性，隱層神經(jīng)元傳遞函數(shù)為log-signoid 函數(shù)。神經(jīng)網(wǎng)絡算法表達式如式（6）：

式中：NNout為神經(jīng)網(wǎng)絡輸出；N11為單位轉(zhuǎn)速，w1i、w2i、w3i為隱層神經(jīng)元權值；bi為隱層神經(jīng)元閾值；w0i為輸出層權值；b0為輸出層閾值；W1、W2為權值向量。

通過神經(jīng)網(wǎng)絡學習可得到網(wǎng)絡相應權值與閾值，從而建立基于神經(jīng)網(wǎng)絡的水輪機模型，實際應用中，只需輸入單位轉(zhuǎn)速、導葉開度和槳葉角度，利用已建立的水輪機神經(jīng)網(wǎng)絡模型，就可計算水輪機輸出流量和力矩。水輪機神經(jīng)網(wǎng)絡模型如圖4所示。

圖4 神經(jīng)網(wǎng)絡模型方框圖Fig.4 Neural network model block diagram

神經(jīng)網(wǎng)絡樣本數(shù)據(jù)獲取方法如下：

神經(jīng)網(wǎng)絡建模的依據(jù)為機組模型綜合特性曲線，引入單位流量、單位轉(zhuǎn)速和單位力矩如式：

式中：M11為單位力矩；D為水輪機直徑；H為水輪機水頭；M為水輪機力矩；Q11為單位流量；Q為水輪機流量；n11為單位轉(zhuǎn)速；n為水輪機轉(zhuǎn)速。

單位轉(zhuǎn)速n11、單位流量Q11、導葉開度α、槳葉角度β以及效率η是神經(jīng)網(wǎng)絡訓練的輸入數(shù)據(jù)。首先以n11和α為輸入，單位流量Q11為輸出得到第一組數(shù)據(jù)，訓練出的神經(jīng)網(wǎng)絡即Qnet。其次以n11和Q11為輸入，效率η為輸出得到第二組數(shù)據(jù)。以一個擬合函數(shù)或者BP神經(jīng)網(wǎng)絡，使用第二組數(shù)據(jù)建立輸入n11、Q11與輸出η的關系式（10）。坐標變換后進行神經(jīng)網(wǎng)絡訓練。利用關系式（10）與關系式（11），計算得出輸入導葉開度α與單位轉(zhuǎn)速n11對應輸出單位力矩M11，即為神經(jīng)網(wǎng)絡Mnet。

式中：η為水輪機效率；α為導葉開度；γ=ρ g，γ為容重。

1.4 插值法水輪機模型

插值法水輪機模型與神經(jīng)網(wǎng)絡類似，也可描述水輪機非線性特性。神經(jīng)網(wǎng)絡模型是通過數(shù)據(jù)處理訓練得到權值與閾值，在仿真時通過歸一化后的單位轉(zhuǎn)速n11與導葉開度α進行矩陣計算，再反歸一化得出輸出單位力矩與單位流量。插值法模型通過模型綜合特性曲線數(shù)據(jù)計算得到插值表，在仿真時輸入單位轉(zhuǎn)速n11與導葉開度α，通過二元三次插值計算輸出力矩與流量。神經(jīng)網(wǎng)絡模型需要前期樣本數(shù)據(jù)多，但實時仿真輸出數(shù)據(jù)誤差相對小，插值法使用的插值表劃分以及不同插值方法都會對最終結果產(chǎn)生影響。

插值法水輪機模型同樣基于模型綜合特性曲線建模，按照水頭和開度劃分工況，計算出對應工況點流量和力矩。仿真時輸入單位轉(zhuǎn)速與開度，采用二元三次插值法進行插值計算出全工況的流量與力矩。

二元三次插值公式為式（12）。

其中k≠i，l≠j，xp＜x＜xp+1，yq＜y＜yq+1。

其原理如圖5所示。具體步驟為：

圖5 二元三次插值示意圖Fig.5 Binary cubic interpolation schematic diagram

設插值點為S，其坐標為S（xs，ys，zs），

（1）在插值點S周圍找出9 個點A、B、C、D、E、F、G、H、I，且滿足以下關系式：

其中：xp＜xs＜xp+1，yq＜ys＜yq+1。

（2）分別對（A、D、G），（B、E、F），（C、F、I）進行一元三點拉格朗日插值，得到點（U、V、W）。

（3）對點（U、V、W）再進行一次一元三點拉格朗日插值，得點S，S點的值zs，s即為該插值點的函數(shù)值。

仿真時插值法模型與神經(jīng)網(wǎng)絡模型相似，都是輸入單位轉(zhuǎn)速與導葉開度，神經(jīng)網(wǎng)絡模型通過s 函數(shù)里編寫的歸一化及矩陣計算代碼計算得出輸出力矩與流量，插值法模型通過simulink 中2-D lookup table 模塊實現(xiàn)插值，輸出力矩與流量。模型方框圖如圖6所示。

圖6 插值法模型方框圖Fig.6 Block diagram of interpolation model

2 仿真及對比實驗

2.1 模型參數(shù)

某機組為混流式機組，該機組各參數(shù)如表1，神經(jīng)網(wǎng)絡權值與閾值如表2，神經(jīng)網(wǎng)絡輸入輸出關系曲面如圖7 所示，插值表輸入輸出關系如圖8所示。圖中X軸為單位轉(zhuǎn)速，Y軸為導葉開度，Z軸為輸出流量或力矩。

圖7 流量與力矩神經(jīng)網(wǎng)絡Fig.7 Flow and torque neural network

圖8 流量與力矩插值表Fig.8 Flow and torque interpolation table

表1 額定工況參數(shù)Tab.1 Rated condition parameter

表2 神經(jīng)網(wǎng)絡權值與閾值Tab.2 Neural network weights and thresholds

調(diào)速器為并聯(lián)型調(diào)速器，PID 參數(shù)為KP=0.20，KI=0.05，KD=0.05。分別以功率上調(diào)1%和10%為階躍輸入得到仿真力矩輸出如圖9（a）所示，流量輸出如圖9（b）所示，水頭輸出如圖9（c）所示。

2.2 仿真結果分析

圖9（a）為模型力矩輸出，如圖所示其中4 種模型穩(wěn)態(tài)力矩相近，只有動態(tài)過程存在一定區(qū)別，其中分段線性化模型動態(tài)過程最快，理想水輪機模型動態(tài)過程最慢。圖9（b）為模型流量輸出，如圖所示4種模型流量波動基本不同，其中理想水輪機模型的波動量最小，而神經(jīng)網(wǎng)絡模型和插值法模型的波動量最大，且穩(wěn)態(tài)值相近。圖9（c）為模型水頭輸出，如圖9所示除插值法外水頭波動穩(wěn)態(tài)后基本相同，插值法在波動開始時水頭產(chǎn)生較大瞬時變化然后趨于平緩。為更便于觀察各參數(shù)變化，引入公式（13），對力矩、流量和水頭變化進行對時間求導后，除以功率變化，處理得圖10和11。

圖9 仿真輸出Fig.9 The simulation output

式中：t為時間；ΔP為機組功率變化；φ1為力矩導數(shù)對功率變化相對值；φ2為流量導數(shù)對功率變化相對值；φ3為水頭導數(shù)對功率變化相對值。

由圖10 可知在功率波動初期，插值法水輪機模型的φ1、φ2與φ3經(jīng)過較大波動然后趨于平緩，其他3 個模型的φ1與φ2參數(shù)相近，只有φ3存在較大差別。由不同功率波動對比可知：功率波動對于φ波形變化影響較小。由圖11可知：分段線性化模型仿真時間約為0.22 s，理想水輪機模型為0.21 s，插值法模型接近0.53 s，約為前兩種方法二點五倍，神經(jīng)網(wǎng)絡模型仿真時間更長，為18 s左右。

圖10 不同功率波動時φ值變化Fig.10 The φ value changes with different powerfluctuations

圖11 實際仿真時間Fig.11 Actual simulation time

3 結 語

本文建立了4種水輪機模型，并進行了仿真對比分析，結果表明：4種不同模型都能有效仿真力矩穩(wěn)態(tài)特性，動態(tài)過程相差不大，而流量仿真時4 種模型輸出相差較大。進行水頭仿真時相比于其他3種方法，插值法模型水頭波動會偏大，水頭波動反映到蝸殼水壓后會造成仿真蝸殼水壓大于實際水壓，進而使調(diào)節(jié)保證計算產(chǎn)生一定誤差。仿真穩(wěn)態(tài)力矩時四種模型精度相近，在動態(tài)過程存在一定差別，在波動瞬間若超出插值表邊界處會產(chǎn)生水頭瞬時較大波動，產(chǎn)生一定誤差。

分段線性化模型在需要仿真時間較短和建模更簡單時是較好的選擇，當更多考慮動態(tài)過程精確度時，選擇神經(jīng)網(wǎng)絡模型會更好。 □