杜 君，黃定波，鄭程遙

(1.五大連池市山口湖開發(fā)建設(shè)中心 ，黑龍江 五大連池 164500;2.廣州市恩萊吉能源科技有限公司，廣東 廣州 510000;3.廣東水利水電職業(yè)技術(shù)學(xué)院，廣東 廣州 510925)

1 概況

水輪發(fā)電機(jī)組加速時(shí)間常數(shù)Ta較小或者水輪機(jī)Hs較小時(shí)，在甩負(fù)荷過渡過程中，常伴有抬機(jī)現(xiàn)象，損壞機(jī)組零部件，影響水電站經(jīng)濟(jì)效益[1]。山口水電站位于黑龍江省五大連池市嫩江一級(jí)支流訥謨爾河上游，原裝機(jī)為2×13 MW，于2016年增效擴(kuò)容改造至15 MW。1#機(jī)組改造后在甩100%負(fù)荷試驗(yàn)時(shí)，出現(xiàn)了較嚴(yán)重的抬機(jī)現(xiàn)象。事后查明，轉(zhuǎn)動(dòng)部分上抬32 mm，部分零件受損。為防止電站抬機(jī)的發(fā)生，本研究基于水輪機(jī)水泵工況抬機(jī)機(jī)理，提出抬機(jī)的抑制方法及導(dǎo)葉關(guān)閉策略。通過內(nèi)特性法模擬仿真，選擇了最優(yōu)關(guān)閉策略，使之既滿足調(diào)保計(jì)算要求，又滿足不抬機(jī)的要求。電站引水系統(tǒng)及機(jī)組參數(shù)見表1，水輪發(fā)電機(jī)組如圖1所示。

表1 引水系統(tǒng)及機(jī)組參數(shù)

圖1 水輪發(fā)電機(jī)組示意

2 水輪機(jī)水泵工況抬機(jī)機(jī)理及抑制方法

目前，水力機(jī)組抬機(jī)的機(jī)理已基本清楚，一是水輪機(jī)在甩負(fù)荷過程中，將經(jīng)歷制動(dòng)工況及水泵工況，當(dāng)導(dǎo)葉開度關(guān)至0時(shí)，引起最大的反向水推力[2]；二是反水錘。當(dāng)導(dǎo)葉快速關(guān)閉，由于尾水管內(nèi)水流的慣性，轉(zhuǎn)輪區(qū)會(huì)出現(xiàn)較大的真空，水流脫離轉(zhuǎn)輪而“斷流”，這樣在尾水慣性過程完結(jié)后，在下游水壓作用下，反向的尾水水流被高速旋轉(zhuǎn)的轉(zhuǎn)輪制動(dòng)，產(chǎn)生較大的反向水錘力。反向水推力與反向水錘力都會(huì)形成抬機(jī)力。

反水錘抬機(jī)力的計(jì)算，相對簡單，本文僅討論反向水推力引起的抬機(jī)力。

在水輪機(jī)組甩負(fù)荷過程中，水輪機(jī)各參數(shù)如機(jī)組轉(zhuǎn)速上升值、蝸殼壓力上升值、尾水管真空值及軸向水推力等，皆為時(shí)間t的函數(shù)，故將水輪機(jī)基本方程表達(dá)為式(1)，將水輪機(jī)軸向水推力表達(dá)為式(2)[3]：

(1)

(2)

式中V1u(t)和V2u(t)分別為t時(shí)轉(zhuǎn)輪葉片進(jìn)、出口水流質(zhì)點(diǎn)的絕對速度在圓周方向的分量；r1與r2分別為轉(zhuǎn)輪進(jìn)出口計(jì)算半徑；η(t)和H(t)分別為t時(shí)轉(zhuǎn)輪的效率與進(jìn)出水邊的水壓差；ω(t)為t時(shí)水輪機(jī)葉片旋轉(zhuǎn)的角速度。式(2)中，K(t)為軸向水推力系數(shù)，D1為水輪機(jī)直徑，Pt為軸向水推力。

水輪機(jī)甩負(fù)荷，導(dǎo)葉關(guān)閉，從運(yùn)行開度至0開度，此過程中，由于水流慣性形成波動(dòng)過程，出口水流圓周速度分量的減少，總是滯后于進(jìn)口圓周速度分量的減少，故V1u(t)r1-V2u(t)r2是遞減的。甩負(fù)荷時(shí)，ω(t)是連續(xù)的，故必有一最大值，此即為最大轉(zhuǎn)速升高值。過此值后，V1u(t)r1-V2u(t)r2繼續(xù)減少，w(t)下降，但只要V1u(t)r1-V2u(t)r2>0，機(jī)組則仍處于水輪機(jī)工況。V1u(t)r1-V2u(t)r2=0時(shí)，水壓力為0，水輪機(jī)軸向水推力為0。當(dāng)V1u(t)r1-V2u(t)r2<0，水壓力反向，機(jī)組即處于制動(dòng)工況及水泵工況，當(dāng)V1u(t)r1=0，即導(dǎo)葉全關(guān)時(shí)，理論和實(shí)測都表明，H(t)<0，且反向水推力最大。此反向水推力就構(gòu)成了較大的抬機(jī)力。

水輪機(jī)工況—制動(dòng)工況—水泵工況的過渡過程及其抬機(jī)機(jī)理見圖2。

從圖2可見，水輪機(jī)室大波動(dòng)時(shí)抑制抬機(jī)力的辦法，主要在于延緩導(dǎo)葉過零的時(shí)間降低導(dǎo)葉過零時(shí)的轉(zhuǎn)速，因?yàn)閺乃美碚摽芍?，最大反向水推力，與導(dǎo)葉過0時(shí)的轉(zhuǎn)速平方成正比。內(nèi)特性法可模擬關(guān)閉過程以選擇壓力上升值和轉(zhuǎn)速上升值符合規(guī)范要求，又可保證不發(fā)生抬機(jī)的關(guān)閉規(guī)律。

圖2 水輪機(jī)大波動(dòng)過程水泵工況抬機(jī)機(jī)理

3 內(nèi)特性法的數(shù)學(xué)模型

文獻(xiàn)[3]給出了軸流式水輪機(jī)甩負(fù)荷過度過程計(jì)算的數(shù)學(xué)模型和計(jì)算方法。

(3)

(4)

(5)

(6)

(7)

(8)

式(3)～(8)中Mc為水輪機(jī)的靜態(tài)力矩；MH為水輪機(jī)的動(dòng)態(tài)力矩；Qc為水輪機(jī)靜態(tài)流量；qH為動(dòng)態(tài)流量；ωc為靜態(tài)角速度；ωH為動(dòng)態(tài)角速度；D為轉(zhuǎn)輪直徑；ΩJ為水流旋轉(zhuǎn)慣性常數(shù)；ΩM為水流的流動(dòng)慣性常數(shù)；ζp為水輪機(jī)裝置水頭相對升高值；r為轉(zhuǎn)輪中間流面半徑；F為轉(zhuǎn)輪過水?dāng)嗝婷娣e；b0為導(dǎo)葉高度；α為導(dǎo)葉出口水流角；β0為翼柵的水流零向角；σp為水輪機(jī)裝置引排水系統(tǒng)管道特性系數(shù)；J為機(jī)組轉(zhuǎn)動(dòng)慣量；τ為相對時(shí)間；ρ為水體密度；ηc為水輪機(jī)靜態(tài)效率；qH為水輪機(jī)相對流量；Kp為翼柵的透明系數(shù)；Hz0為水輪機(jī)裝置靜態(tài)初始水頭。

模型的計(jì)算方法在文獻(xiàn)[3]中有詳細(xì)論述，不再贅述。

4 過渡過程計(jì)算

山口水電站1#機(jī)組甩100%負(fù)荷試驗(yàn)時(shí)，出現(xiàn)了較嚴(yán)重的抬機(jī)現(xiàn)象。實(shí)測的甩負(fù)荷波形如圖3所示。從圖3中可看出，導(dǎo)葉開度迅速關(guān)至0時(shí)，機(jī)組頻率約為52.7 Hz，機(jī)組轉(zhuǎn)速為225.9 rpm。如前論述，導(dǎo)葉全關(guān)時(shí)將引起最大的反向水推力，引起抬機(jī)。

圖3 1#機(jī)甩負(fù)荷實(shí)際波形示意

為防止抬機(jī)的發(fā)生，在制定導(dǎo)葉關(guān)閉規(guī)律時(shí)，不僅要滿足壓力上升值和轉(zhuǎn)速上升值符合規(guī)范要求，同時(shí)要延緩導(dǎo)葉過0的時(shí)間和降低導(dǎo)葉過0時(shí)的轉(zhuǎn)速。根據(jù)運(yùn)行要求，甩負(fù)荷后，應(yīng)將機(jī)組限制在空載或空轉(zhuǎn)運(yùn)行，正確的選擇是讓導(dǎo)葉停留在空載開度，避免負(fù)水推力而抬機(jī)。因此，對山口電站2#機(jī)組甩負(fù)荷試驗(yàn)，在制定甩負(fù)荷關(guān)閉規(guī)律時(shí)，采取導(dǎo)葉關(guān)至空載開度后轉(zhuǎn)入PID控制的關(guān)閉策略以避免導(dǎo)葉快速過0。下面根據(jù)內(nèi)特性法對表2工況采用表3關(guān)閉策略進(jìn)行計(jì)算仿真，仿真結(jié)果見表4，仿真波形如圖4所示。

表2 甩負(fù)荷前工況參數(shù)

表3 甩負(fù)荷關(guān)閉策略

表4 仿真結(jié)果

根據(jù)規(guī)范，設(shè)計(jì)水頭小于40 m的電站，甩負(fù)荷時(shí)，最大壓力上升允許值為50%～70%Hmax(電站最大水頭);最大轉(zhuǎn)速上升率不超過額定轉(zhuǎn)速的60%；尾水管最大真空度不超過8m[4]。從表4可知，轉(zhuǎn)速、壓力、尾水管真空度指標(biāo)都滿足規(guī)范要求。但在策略2、3下，導(dǎo)葉過0時(shí)轉(zhuǎn)速偏高，分別為221.3 rpm、232.1 rpm，根據(jù)1#機(jī)組的試驗(yàn)結(jié)果，反向水推力較大，存在抬機(jī)風(fēng)險(xiǎn)；策略1延緩了導(dǎo)葉過0的時(shí)間，有效降低了導(dǎo)葉過零時(shí)的轉(zhuǎn)速，導(dǎo)葉過0時(shí)轉(zhuǎn)速接近額定轉(zhuǎn)速214.3 rpm，抬機(jī)力較小，抬機(jī)風(fēng)險(xiǎn)較小。而策略4采用策略1的兩段關(guān)閉規(guī)律，僅在導(dǎo)葉關(guān)至空載開度時(shí)自動(dòng)轉(zhuǎn)入PID控制，這樣既保證轉(zhuǎn)速、壓力、尾水管真空度指標(biāo)滿足規(guī)范要求，又避開了導(dǎo)葉過0工況，達(dá)到防抬機(jī)的目的。

從仿真結(jié)果可見，策略圖4d導(dǎo)葉關(guān)閉時(shí)間略長，且在空載開度時(shí)，投入PID條件，轉(zhuǎn)速緩慢降低至額定轉(zhuǎn)速，近似于空載運(yùn)行，不可能抬機(jī)。同時(shí)轉(zhuǎn)速上升值，壓力上升值，尾水管真空度都符合要求，因此，選擇策略圖4d。

a 策略1

圖5為2#機(jī)按策略4的關(guān)閉規(guī)律進(jìn)行甩負(fù)荷試驗(yàn)實(shí)測圖，圖5中可看出，理論值與實(shí)測值基本吻合，而且確實(shí)沒有發(fā)生抬機(jī)。

圖5 2#機(jī)甩負(fù)荷波形對比示意

5 運(yùn)行效果分析

一般情況下，過渡過程計(jì)算采用外特性法，其須用到水輪機(jī)全特性曲線，將遇到2個(gè)困難，一是該類曲線難以獲取，二是計(jì)算條件容易失真和計(jì)算過程復(fù)雜。

內(nèi)特性法根據(jù)初始工況，將甩負(fù)荷過渡過程的動(dòng)態(tài)參數(shù)與水輪機(jī)流道內(nèi)各元件的物理結(jié)構(gòu)和動(dòng)態(tài)位置相結(jié)合，建立相應(yīng)的非線性微分方程組，利用已知的邊界條件，用數(shù)值方法求得方程組的解，從而得到各動(dòng)態(tài)工況參數(shù)的瞬時(shí)值和演化軌跡。

該方法與傳統(tǒng)的外特性法的區(qū)別如圖6與圖7所示。從圖6和圖7對比分析可知，外特性和內(nèi)特性法都是基于水輪發(fā)電機(jī)組(包括流道)的物理結(jié)構(gòu)和位置參數(shù)，前者通過水輪機(jī)模型試驗(yàn)，得到水輪機(jī)全特性曲線(比一般的模型綜合特性曲線復(fù)雜，且不易獲取)，從而間接取得相關(guān)的微分方程組，而后者則直接通過機(jī)組的物理結(jié)構(gòu)和位置參數(shù)直接求得相關(guān)的微分方程組。因此，無論從簡單性、可操作性，靈活性、全面性考慮，后者明顯優(yōu)于前者。

圖6 外特性法框示意

圖7 內(nèi)特性法框示意

甩負(fù)荷工況引起的反向水推力即水泵工況引起的抬機(jī)力，可通過延緩導(dǎo)葉過0時(shí)間和降低導(dǎo)葉過0時(shí)的轉(zhuǎn)速，或者通過導(dǎo)葉關(guān)至空載開度后轉(zhuǎn)入PID控制的關(guān)閉策略，避免導(dǎo)葉快速過0等關(guān)閉規(guī)律有效抑制?；诖朔椒芍贫ǚ捞C(jī)的關(guān)閉策略，并用內(nèi)特性法進(jìn)行數(shù)值模擬，最后選擇最優(yōu)的關(guān)閉策略，既滿足調(diào)保計(jì)算要求，又滿足不抬機(jī)的要求。從山口電站2#機(jī)組的甩負(fù)荷試驗(yàn)可知，理論值與實(shí)測值基本吻合，且未發(fā)生抬機(jī)。因此，內(nèi)特性法是過渡過程仿真計(jì)算的一種有效方法，對水電站調(diào)保計(jì)算和甩負(fù)荷試驗(yàn)具有一定的指導(dǎo)意義。

6 結(jié)語

1) 甩負(fù)荷工況引起的反向水推力和反水錘力等抬機(jī)力都可以通過選擇調(diào)節(jié)系統(tǒng)的關(guān)閉規(guī)律(α(t)，β(t))得到有效的抑制。

2) 水輪機(jī)甩負(fù)荷時(shí)抑制抬機(jī)力的辦法，主要在于延緩導(dǎo)葉過零的時(shí)間和降低導(dǎo)葉過0時(shí)的轉(zhuǎn)速，也可通過導(dǎo)葉關(guān)至空載開度后轉(zhuǎn)入PID控制的關(guān)閉策略避免導(dǎo)葉過零，降低抬機(jī)風(fēng)險(xiǎn)。

3) 內(nèi)特性法是雙調(diào)節(jié)水輪機(jī)過度過程仿真計(jì)算的一種有效方法，對水電站調(diào)保計(jì)算和甩負(fù)荷試驗(yàn)具有一定的指導(dǎo)意義。