高曉光，漆智鵬通信作者，黃斌斌

(1.華電電力科學研究院有限公司，浙江 杭州 310030；2.貴州烏江水電開發(fā)有限責任公司思林發(fā)電廠，貴州 銅仁 565109)

0 引言

過渡過程是指一個動力系統(tǒng)在其自身特性變化或在外力干擾的作用下，由一個穩(wěn)定狀態(tài)到達另一個狀態(tài)的過程[1]。計算機技術(shù)的飛速發(fā)展為水力過渡過程研究提供了強有力的技術(shù)支持。結(jié)合數(shù)值分析的研究成果[2]，計算機技術(shù)使得對具有復雜邊界的非線性系統(tǒng)進行水力過渡過程計算和分析成為了可能，從而更好地進行水電站引水發(fā)電系統(tǒng)優(yōu)化設(shè)計和電站運行控制仿真研究[3]，也為水力過渡過程的工程應用展現(xiàn)了前所未有的廣闊前景。

本文建立了水電站水力過渡過程完整的數(shù)學模型以及邊界方程，包括水庫、引水隧洞、壓力鋼管、岔管、調(diào)壓室、調(diào)速器和水輪發(fā)電機組等[4-5]。將模型中的各個部件實現(xiàn)模塊化，在Windows平臺下，采用Visual C++和Fortan語言混合編程的方式實現(xiàn)程序主體代碼，并通過數(shù)據(jù)庫存儲過的特定的數(shù)據(jù)結(jié)構(gòu)使其能夠自由組合，構(gòu)成水電站完整的模型數(shù)據(jù)并進行水力過渡過程計算。

1 數(shù)學模型及數(shù)據(jù)結(jié)構(gòu)

在Windows平臺上開發(fā)水力過渡過程通用計算軟件需要掌握瞬變流的基本方程和計算方法，建立水電站完整的數(shù)學模型以及上下游邊界、調(diào)壓室邊界、水輪機邊界等典型邊界的基本解法。本文從理論上了分析水電站系統(tǒng)各部分組成，以及水電站過渡過程計算軟件開發(fā)的數(shù)據(jù)結(jié)構(gòu)。

1.1 基本方程

有壓非恒定流的基本方程如下。

1)運動方程：

(1)

2)連續(xù)方程：

(2)

式中：H為壓頭；g為重力加速度；a為水擊波速；V為水流速；f為沿程損失系數(shù)；D為管徑；t為時間。

由運動方程和連續(xù)方程可以得到特征線方程：

C+:HP=CP-B×QP

(3)

C-:HP=CM+B×QP

(4)

CP=HA+B×QA-R×|QA|×QA

(5)

CM=HB-B×QB+R×|QB|×QB

(6)

(7)

(8)

式中：HA、HB分別為t-Δt時刻管段第i-1和第i+1節(jié)點處的壓頭；QA、QB分別為t-Δt時刻管段第i-1和第i+1節(jié)點處的流量；CP、CM分別與t-Δt時刻的壓頭和流量有關(guān)，t時刻是已知量；Δx為相鄰節(jié)點間的距離；R為阻力系數(shù)。

1.2 邊界方程

1.2.1 上游邊界

上游邊界方程(考慮進水口的局部水頭損失)：

(9)

C-:HP=CM1+B1×QP

(10)

式中：k1為進口損失系數(shù)；下標1為1號管的參量。

1.2.2 調(diào)壓室

調(diào)壓室邊界方程主要有：

C+:HP1=CP1-B1×QP1

(11)

C-:HP2=CM2+B2×QP2

(12)

HP1=HP2=HPS

(13)

QP1-QP2=QS

(14)

(15)

(16)

式中：下標1為調(diào)壓室前隧洞的參量；下標2為調(diào)壓室后壓力鋼管的參量；Qs為流進調(diào)壓室的流量；k為流進調(diào)壓室的阻力系數(shù)；Z為調(diào)壓室水位。

1.2.3 水輪機邊界

水輪機邊界條件主要方程為：

HPU=CP1-B1×QP

(17)

HPD=CM2+B2×QP

(18)

(19)

(20)

(21)

(22)

ω=2πn

(23)

1.2.4 下游河道

下游河道的邊界方程(考慮尾水管出口的局部水頭損失)：

(24)

C+:HP=CPS-BS×QP

(25)

1.3 數(shù)據(jù)結(jié)構(gòu)

本文采用Visual C++與Fortran語言混合編程，同時以Access數(shù)據(jù)庫作為存儲媒介。軟件中引水發(fā)電系統(tǒng)各個組件由元件來反映，并采用樹來存儲各個元件的屬性及協(xié)聯(lián)關(guān)系。

1.3.1 元件

每個元件都含有一個水庫對象、管道對象、調(diào)壓室對象、水輪機對象以及尾水對象。各元件之間的關(guān)系由元件指針指出，用于訪問前一元件或者后面的元件。其中，父節(jié)點指向唯一的前一元件，3個子節(jié)點指向后面的元件。具體元件的定義如下。

class CUnit

{

public:

CUnit()；

virtual～CUnit()；

CPipe pipe； //管道對象

CFm fm； //閥門對象

CSur sur； //調(diào)壓室對象

CTur tur； //水輪機對象

CSk sk； //水庫對象

CWs ws； //尾水對象

CString type； //元件類型

CString nname； //下一元件名稱

CString name； //當前元件名稱

CUnit *f； //父節(jié)點

CUnit *Lchild； //子節(jié)點

CUnit *Mchild； //子節(jié)點

CUnit *Rchild； //子節(jié)點

}；

1.3.2 樹

樹是由一個集合以及在該集合上定義的一種層次關(guān)系所構(gòu)成的一個具有層次結(jié)構(gòu)的集合。集合中的元素是樹的節(jié)點，節(jié)點間的關(guān)系為父子關(guān)系。創(chuàng)建樹的過程是通過元件間的關(guān)系，由樹根(即水庫)開始，讀取元件的下一節(jié)點名稱并開辟新內(nèi)存，給內(nèi)存賦入相對應的值。引水發(fā)電系統(tǒng)的水力過渡過程計算實際上就是以水庫為根節(jié)點的樹的計算。以樹結(jié)構(gòu)進行水電站引水發(fā)電系統(tǒng)模型的存儲，在數(shù)據(jù)結(jié)構(gòu)和算法上都是比較合理的，節(jié)約了程序?qū)ぶ窌r間，很大程度上提高了程序的運行效率。

圖1 樹的結(jié)構(gòu)圖

2 實例計算

使用上述研究開發(fā)的水電站水力過渡過程通用軟件對某水電站進行甩負荷大波動過渡過程計算，驗證軟件的有效性。某水電站設(shè)計水頭86.2 m，引水隧洞長15.2 km，直徑8.0 m，隧洞末端布置調(diào)壓井，調(diào)壓井底部隧洞底高程2 278.0 m，調(diào)壓井底板高程2 286.4 m，頂部高程2 345.0 m。壓力鋼管φ 6.4 m，分叉引水至三臺混流式水輪發(fā)電機組，支管φ 3.8 m，水庫正常蓄水位2 315.6 m，尾水位2 215.0 m。

將水電站引水系統(tǒng)和水輪機模型相關(guān)數(shù)據(jù)輸入至水電站水力過渡過程通用軟件中(見圖2)，選擇相應的工況進行計算分析。本次選擇大波動過渡過程工況進行計算，設(shè)置計算時間為1 800 s，具體計算結(jié)果如圖3～4所示。

圖2 水力過渡過程通用計算軟件界面

圖3 調(diào)壓室水位變化過程

圖4 轉(zhuǎn)速變化過程

從圖3和圖4可以看出，甩負荷瞬間，調(diào)壓室水位開始上升，223 s左右達到最大涌浪水位2 333.8 m，隨后水位波幅不斷衰減。機組最大轉(zhuǎn)速上升發(fā)生在甩最負荷時﹐最大轉(zhuǎn)速為426 r/min。計算結(jié)果能完整地反映機組甩負荷大波動過渡過程，驗證了軟件的有效性。

3 結(jié)語

本文結(jié)合有壓非恒定流的基本方程和特征線方程，通過基于Windows平臺的混合編程，開發(fā)具有通用建模功能的水力過渡過程計算軟件。計算結(jié)果表明：軟件對于常規(guī)引水發(fā)電系統(tǒng)水電站工程具有較強的通用性，可以高效應用于水電站過渡過程仿真計算，指導工程實踐。