孔昭年，許慶進，周同旭，張曉峰，尤永陶，熊曉蕾，程廣蕾，王柏柏，郭珺瑤

（1.中國水利水電科學(xué)研究院，北京 100044；2.天津電氣科學(xué)研究院有限公司，天津 300186；3.中國電力工程有限公司，北京 100044）

基于Simulink的水輪機調(diào)節(jié)系統(tǒng)通用仿真程序的研發(fā)

孔昭年1，許慶進2，周同旭2，張曉峰3，尤永陶3，熊曉蕾2，程廣蕾2，王柏柏2，郭珺瑤2

（1.中國水利水電科學(xué)研究院，北京 100044；2.天津電氣科學(xué)研究院有限公司，天津 300186；3.中國電力工程有限公司，北京 100044）

基于Simulin開發(fā)了水輪機調(diào)節(jié)系統(tǒng)通用仿真程序，該程序用特征矩陣描敘水輪機非線性特性；提出了各種類型的水擊數(shù)學(xué)模型；配合水輪機PID調(diào)速器開展水輪機自動調(diào)節(jié)系統(tǒng)動態(tài)特性仿真計算研究：自動開機、空載擾動、增負荷調(diào)節(jié)及甩負荷控制等。

水輪機；調(diào)節(jié)系統(tǒng)；仿真程序；非線性特性；水擊數(shù)學(xué)模型

在國家能源局“替代調(diào)壓井的新型調(diào)壓閥及其控制系統(tǒng)研究與電站示范應(yīng)用”科技項目的安排下，中國水利水電科學(xué)研究院與天津電氣科學(xué)研究院有限公司合作開發(fā)了基于Simulink的水輪機調(diào)節(jié)系統(tǒng)通用仿真程序，以滿足工程建設(shè)對水輪機調(diào)節(jié)系統(tǒng)仿真計算的需求。Simulink是MATLAB中的一種可視化仿真工具，是實現(xiàn)動態(tài)系統(tǒng)建模、仿真和分析的一個軟件包，被廣泛應(yīng)用于線性系統(tǒng)、非線性系統(tǒng)、數(shù)字控制及數(shù)字信號處理的建模和仿真中，本研究的總體目標(biāo)是針對調(diào)壓閥代替調(diào)壓井的核心問題展開。水輪機調(diào)節(jié)系統(tǒng)通用框圖如圖1所示。

圖1 混流式水輪機調(diào)節(jié)系統(tǒng)原理圖

由圖1可發(fā)現(xiàn)該圖的核心是非線性特性的水輪機特性，文獻[1]提出了解決此問題的新方法，據(jù)此方法我們開發(fā)了配套程序。技術(shù)核心之一就是研發(fā)能真實反映水輪機動態(tài)特性、調(diào)用方便、能包含水輪機型譜及其他所用常用水輪機特性、數(shù)學(xué)方法先進、在計算中總體耗時小、利于推廣的數(shù)據(jù)文件。在此文件中包含有納入水輪機型譜的所有水輪機特性。

水輪機流量和力矩的特性矩陣應(yīng)在仿真數(shù)據(jù)準備階段根據(jù)有關(guān)數(shù)據(jù)表求得，對于軸流轉(zhuǎn)槳式水輪機對應(yīng)每個定槳特性相類似地求取水輪機流量和力矩的特性矩陣；在實時仿真的主程序段，只要已知某一時刻的x11t和at就可很快計算出該時間的單位流量和單位力矩：

軟件結(jié)構(gòu)設(shè)計時將水輪機流量、力矩特征矩陣的運算求取應(yīng)放在數(shù)據(jù)準備階段進行。根據(jù)我們的經(jīng)驗，在確定的實時仿真計算步長時間內(nèi)，取n=5可獲得滿意的計算結(jié)果，由于式（3）式（4）是參數(shù)平面上的連續(xù)函數(shù)，據(jù)此進行的水擊計算更可保證不會發(fā)生不收斂現(xiàn)象。以一個具體的混流式水輪機展示它的應(yīng)用。

根據(jù)水輪機綜合特性曲線，采用相對值法可有圖2所示混流式水輪機相對單位力矩及相對單位流量特性，據(jù)水輪機單位流量特性參數(shù)表和水輪機單位力矩特性參數(shù)表；根據(jù)此兩表及式（1）（2）可計算出該水輪機的流量、力矩的特征矩陣，從而可求水輪機調(diào)節(jié)系統(tǒng)非線性微分方程式組(1)。

觀察圖1可以看出水輪機調(diào)節(jié)系統(tǒng)的另一個重要環(huán)節(jié)是水電站引水系統(tǒng)的動力學(xué)特性，文獻[2，3]已將引水系統(tǒng)數(shù)學(xué)模型的標(biāo)準化形式加以匯總：

圖2 混流式水輪機相對單位力矩及相對單位流量特性

進一步觀察圖1，該圖為水輪機調(diào)速器部分，圖中取緩沖型調(diào)速器，在調(diào)速器部分共有5個非線性環(huán)節(jié)：環(huán)節(jié)2、5表征調(diào)節(jié)器、接力器輸出控制在0～1.0；環(huán)節(jié)3表征緩沖回路輸出限幅，通常值為-0. 2～0.2；環(huán)節(jié)1、4是本文提出的水輪機調(diào)速器引導(dǎo)閥、主配壓閥的非線性，它是典型的飽和型非線性，可進一步分析如下。

通常導(dǎo)葉開啟時間與關(guān)閉時間相等Tf=Tg，其值通常為4～10 s；定義Ty，水輪機接力器時間常數(shù)，其值通常為0.1～0.3 s。定義Ke=Ty/Tf，配壓閥的飽和系數(shù)，如前述其值通常為0.01～0.05，可見水輪機調(diào)速器配壓閥的飽和特性是非常嚴重的。在應(yīng)用Simulink進行計算分析這種非線性環(huán)節(jié)已毫無困難。

利用Simulink十分便于計算分析水輪機過渡過程,如：可靈活地增加控制回路即可計算分析轉(zhuǎn)槳式水輪機的過渡過程特性；在圖1的基礎(chǔ)上水輪機非線性環(huán)節(jié)增加調(diào)壓閥回路即可用來計算分析帶有調(diào)壓閥的水輪機調(diào)節(jié)系統(tǒng)的動態(tài)特性。

計算實例 完成開發(fā)后項目組承擔(dān)了實際工程

表1 某型號水輪機特性矩陣

（1）剛性水錘

適用條件：Tf/Tr＞6及Tw/Tr＞2；

（2）近似彈性水錘

適用條件：3＜Tf/Tr＜6，0.5＜Tw/Tr＜2；

（3）彈性水錘

其中hw=Tw/Tr；水輪機管道特性系數(shù)；

適用條件：Tf/Tr＜3，Tw/Tr＜0.5；

（4）帶調(diào)壓井的引水系統(tǒng)數(shù)學(xué)模型

式中：的計算分析工作，工程參數(shù)如下：

圖3 混流式水輪發(fā)電機組增負荷自動調(diào)節(jié)過程的仿真結(jié)果

圖4 混流式水輪機空載擾動自動調(diào)節(jié)過程的仿真結(jié)果

電站水系統(tǒng)具有調(diào)壓井Tw1=3.9s；Tw2=1.1s；Te=345.8s；a1=0.046s；a2=0.056s；Ta=4.5s；

分別表征隧洞、管道和調(diào)壓井的時間常數(shù)，a1；a2表征隧洞、管道相對水頭損失；水輪機流量和力矩特性載于圖2。

圖3載有混流式水輪發(fā)電機組增負荷自動調(diào)節(jié)過程的仿真結(jié)果;在圖1Δm（t）處施加負荷擾動，穩(wěn)態(tài)負荷M0=0.5，遞增負荷0.1。調(diào)節(jié)時間約15 s；系統(tǒng)穩(wěn)定，動態(tài)性正常。

圖4載有混流式水輪機空載擾動自動調(diào)節(jié)過程的仿真結(jié)果；進行了大量空載擾動仿真計算，限于篇幅僅發(fā)布一幅上擾仿真計算結(jié)果。在圖1 Cf處施加擾動指令信號，擾動量0.04。

圖5載有混流式水輪發(fā)電機組甩全負荷調(diào)節(jié)過程的仿真結(jié)果；由計算結(jié)果可見當(dāng)Tf=4.5時機組最大速率上升0.43；最大水壓上升0.31。水壓最大值出現(xiàn)在3.1 s；而最大轉(zhuǎn)速上升出現(xiàn)在3.4 s。

圖5 混流式水輪發(fā)電機組甩全負荷調(diào)節(jié)過程的仿真結(jié)果

圖6 混流式水輪機啟動過渡過程的仿真結(jié)果

在圖6上示出混流式水輪機機組啟動過程的仿真結(jié)果，機組啟動的物理過程：接到啟動命令導(dǎo)葉開至啟動開度20%，機組轉(zhuǎn)速陸續(xù)增大；當(dāng)增大至96%時調(diào)速器投入，系統(tǒng)處于自動調(diào)節(jié)狀態(tài)；由于轉(zhuǎn)速尚有4%的偏差，導(dǎo)葉略自動開啟并單調(diào)地迅速調(diào)節(jié)到穩(wěn)定值1.0，歷時40 s啟動過程結(jié)束。

[1]孔昭年.水輪機控制系統(tǒng)的設(shè)計與計算[M].武漢：長江出版社,2012.

[2]DL/T 1120-2009水輪機調(diào)節(jié)系統(tǒng)自動測試及實時仿真裝置技術(shù)條件[S].

[3]DL/T 1548-2016水輪機調(diào)節(jié)系統(tǒng)設(shè)計與應(yīng)用導(dǎo)則[S].

TV734

A

1672-5387（2017）06-0001-03

10.13599/j.cnki.11-5130.2017.06.001

2017-04-27

孔昭年（1941-），男，教授級高級工程師，從事水輪機控制技術(shù)研究工作。