肇毓鋒
(遼寧省水利事務(wù)服務(wù)中心,遼寧 沈陽 110003)
水能資源是清潔能源,但由于部分電站過分追求經(jīng)濟效益,不合理開發(fā)利用水能資源,由此導(dǎo)致了許多生態(tài)問題,因此科學(xué)、合理、高效利用水能資源,既能滿足電站經(jīng)濟效益,同時能夠保證生態(tài)安全,意義十分重大[1]。水能資源開發(fā)利用受利用水量、經(jīng)濟效益和生態(tài)環(huán)境等動態(tài)因素影響,因此其優(yōu)化問題是一個典型的多目標問題,必須統(tǒng)籌考慮各種制約因素的相互關(guān)系、相互影響,建立符合實際情況的多目標、群決策模型。
石橋水電站位于遼寧省愛河流域,愛河流域水能資源蘊藏量十分豐富,水系發(fā)育良好。流域整體開發(fā)程度較低,沒有大型控制性水利工程,但隨著社會的發(fā)展對能源需求量的不斷加大以及對清潔能源的開發(fā)程度不斷提高,加之早期水能資源開發(fā)未能充分重視生態(tài)環(huán)境因素,已開發(fā)水電站多以引水式為主,導(dǎo)致區(qū)域內(nèi)的水生態(tài)環(huán)境受到了影響。由此可見,采取原始單目標調(diào)度模式,單純以追求電站經(jīng)濟效益為目標,勢必會對流域水生態(tài)環(huán)境產(chǎn)生不良影響。為統(tǒng)籌考慮流域水生態(tài)環(huán)境保護與電站經(jīng)濟效益,解決單一目標調(diào)度無法同時兼顧水生態(tài)保護與經(jīng)濟效益的問題,本文建立了石橋水電站水能資源優(yōu)化開發(fā)利用模型,有效解決了水能資源利用過程中追求經(jīng)濟效益目標與生態(tài)目標存在相互矛盾、相互制約的問題,找到最優(yōu)交集使2個目標同時達到要求,對在滿足生態(tài)目標要求的前提下充分開發(fā)利用水能資源意義重大[2-5]。
石橋水電站為日調(diào)節(jié)水電站,電站工程壩址位于愛河干流的鳳城市東湯鎮(zhèn)石橋村的唐家堡,是愛河流域水能梯級開發(fā)規(guī)劃中,草河河口以下河段水電開發(fā)的第一級水電站工程,在流域水能資源開發(fā)序列中處于比較重要的地位。石橋水電站工程屬于小型水電站工程,以發(fā)電為主,其壩址以上控制流域面積約4 000多km2,占總流域面積的80%以上。
電站廠房內(nèi)布置3臺水輪發(fā)電機組,總裝機容量為9 600 kW。其校核洪水位以下的水庫總?cè)萘考s為3 664萬m3,水庫正常蓄水運行時,總水量為465.5萬m3。本文以石橋電站為研究對象,應(yīng)用水能資源優(yōu)化開發(fā)利用模型對其進行優(yōu)化運行調(diào)度的計算并提出優(yōu)化運行調(diào)度方案[6-9]。通過優(yōu)化運行調(diào)度方案的選擇,可以在滿足石橋水電站發(fā)電效益的基礎(chǔ)上,最大程度地滿足河道生態(tài)流量要求,保證區(qū)域內(nèi)水生態(tài)環(huán)境的健康發(fā)展,平衡經(jīng)濟目標與生態(tài)目標之間的關(guān)系,使得二者在矛盾中得到統(tǒng)一,并對愛河流域內(nèi)其他區(qū)域水能開發(fā)提供參考依據(jù)[10-17]。
在滿足河流生態(tài)需求、保護河流天然形態(tài)為前提下充分開發(fā)利用水能資源過程中,受到水庫的庫容(水位)、發(fā)電流量約束、水電站水頭限制、水輪機組的出力限制等諸多要素影響,而各要素之間還存在矛盾關(guān)系,因此必須依托建立水電站發(fā)電多目標運行實施,即在發(fā)揮水電站社會經(jīng)濟效益的同時,能夠保證水電站下游河段生態(tài)流量的基本要求。水能資源優(yōu)化開發(fā)是水電站多目標運行在生態(tài)層面的集中表現(xiàn),也是水電站多目標運行的有機組成部分。水能資源開發(fā)利用優(yōu)化調(diào)度運行多目標模型目標函數(shù)如下。
a)發(fā)電量最大目標。追求總發(fā)電量最大為目標,充分開發(fā)利用水能資源,同時滿足最大經(jīng)濟效益。
(1)
式中E——河段水流理論發(fā)電能量,kW·h;N——河段水流理論出力,kW;Q——發(fā)電流量,m3/s;H——河段上、下游兩斷面水面高程差,m;Δt——發(fā)電時段長度,h;K——出力系數(shù),取7.5;T——時段,以月為單位進行計算時段的選取,即T取值12,下同。
b)保證出力最大。追求水電站保證出力最大為目標,本文以枯水期水電站最小出力近似為水電站的保證出力,在保證水能資源利用率最高的同時也保證了發(fā)電下泄流量的合理性。
F2=Max(minNt)t=1,2,… ,T
(2)
式中Nt——第t時段水電站的總出力,以t時段最小出力近似為各時段的保證出力。
a)下泄流量約束。為滿足河道必要的生態(tài)環(huán)境用水,保護流域生態(tài)、水文系統(tǒng)的天然狀態(tài),水電站下泄流量不得小于該段河道適宜生態(tài)流量,由于汛期總來水量較大,可以滿足下游河道的生態(tài)流量要求,所以主要考慮枯水時期的水電站下泄流量能否滿足下游河道生態(tài)流量的要求。因此對下泄流量進行如下約束:
Qt≥Qst=1,2,… ,T
(3)
式中Qs——電站下游河道各月適宜生態(tài)流量;Qt——水電站第t月發(fā)電后下游河道流量。
b)水輪機可承受過水量約束。為確保電站水輪機正常運轉(zhuǎn),電站水輪機最大過流量不得超過其設(shè)計值。因此對電站水輪機過水量進行如下約束:
Qt≤QTt=1,2,…,T
(4)
式中QT——水電站達到最大設(shè)計水頭時,水輪機允許發(fā)電時可承受的最大過流量。
c)水庫最低蓄水位限制約束。正常情況下,水庫水位不會低于其死水位,因此對水庫最低蓄水位進行如下約束:
Ht≥Hmint=1,2,…,T
(5)
式中Hmin——t時段水庫水位達到其設(shè)計允許消落的最低水位,一般情況下特指水庫達到死庫容時的水位。
d)防洪限制約束。在汛期過程中,主要針對汛期水庫水量的限制,水庫最高水位即為水庫水位達到防洪限制水位;非汛期過程中,以正常蓄水位進行約束,即水庫在正常運行情況下對應(yīng)的水庫水位。因此對水庫防洪限制進行如下約束:
Ht≤Hmaxt=1,2,…,T
(6)
式中Hmax——第t時段水庫蓄水量的上限水平。
e)水輪機出力限制約束。水電站在其水輪機設(shè)計最小出力及最大出力區(qū)間能夠正常運轉(zhuǎn),因此對水輪機出力進行如下約束:
NN≤Nt≤NTt=1,2,…,T
(7)
式中Nt——第t時段水輪機機組的實際出力值;NN——在滿足水輪機設(shè)計最小出力情況下的發(fā)電出力;NT——水輪機在區(qū)域現(xiàn)有水能資源條件下的最大發(fā)電出力。
f)水量平衡約束。水庫運用調(diào)度過程中,其蓄水量、來水量、泄流量保持動態(tài)平衡關(guān)系,因此對水庫水量平衡進行如下約束:
Vt+1=Vt+(St-Qt)Δtt=1,2,…,T
(8)
式中Vt+1——t+1時段水庫的總水量;Vt——第t時段初水庫的蓄水總量;St、Qt——t時段水電站水庫的來水量及下泄流量;Δt——時段長度。
g)非負約束。所有變量均為非負值。
綜合各目標及約束條件,得出了石橋水電站優(yōu)化運行調(diào)度多目標函數(shù)。
本文所研究的水能資源優(yōu)化開發(fā)是指在考慮水生態(tài)環(huán)境前提下,提高水能資源利用率的最優(yōu)決策。在保證開發(fā)利用水能資源過程中不對水生態(tài)環(huán)境造成嚴重破壞的情況下,使其獲得的經(jīng)濟效益最大化。根據(jù)目前水能資源開發(fā)方面多目標理論研究狀況,有多種方法可以用于模型求解,包括最優(yōu)化方法、差分進化算法、數(shù)學(xué)模擬法等。水能資源優(yōu)化開發(fā)模型具有多個約束條件,是一個復(fù)雜的多目標優(yōu)化模型,對模型的求解較為困難,同時愛河流域水能資源利用具有多目標性和復(fù)雜性,所以本文采用數(shù)學(xué)模擬方法進行研究,即利用計算機仿真模擬技術(shù)來對石橋水電站水能資源優(yōu)化開發(fā)利用模型進行求解。
數(shù)學(xué)模擬方法是將實際案例通過數(shù)學(xué)的方式進行轉(zhuǎn)化,得出可以貼切描述實際情況的數(shù)學(xué)方程式或方程組,通過計算機編程進行計算的一種高度仿真算法。在實際問題較為復(fù)雜且存在多個變量及約束條件時,常規(guī)方法難以進行計算,采用數(shù)學(xué)模擬法可以在短時間內(nèi)進行大量數(shù)據(jù)的求解,有著明顯的優(yōu)勢。同時,數(shù)學(xué)模擬法具有直觀性較強,其原理容易被大多數(shù)人接受及實際應(yīng)用過程可以省去很多的復(fù)雜程序等特點,因此成為近年來應(yīng)用頻率較高的算法之一。
采用計算機仿真模擬技術(shù)來對石橋水電站水能資源優(yōu)化開發(fā)利用模型進行求解的步驟如下。
步驟一確定石橋水電站的初始調(diào)度狀態(tài)。根據(jù)石橋水電站壩高及水電站防洪限制水位,在典型年來水條件下,確定石橋水電站各月入庫流量及經(jīng)計算后的下游河道生態(tài)流量。根據(jù)石橋水電站H-W(水位-庫容)曲線,可得水庫的水位與庫容之間關(guān)系。聯(lián)系典型年來水條件下水庫庫容與來水量及泄流量關(guān)系,下泄流量初始值設(shè)置為0,單位m3/s,迭代步長為0.1,進行試算。
步驟二約束條件判定。模型共有7個約束條件,在保證所有變量值非負的前提下,分別從生態(tài)流量約束、水庫約束、電站出力約束三方面進行判定。其中生態(tài)流量約束主要對下泄流量進行判定;水庫約束主要對水庫水位最大值、最低值及水量平衡進行判定;電站出力約束主要對水輪機最大過流量、水輪機出力上下限進行判定。具體過程如下:①生態(tài)流量判定。根據(jù)水庫的下泄流量,確定石橋水電站下游斷面流量并判斷各時段斷面流量是否滿足生態(tài)流量要求。如相應(yīng)時段內(nèi)能夠滿足下游河道生態(tài)流量要求,則繼續(xù)進行判定,否則返回步驟一重新確定下泄流量。②水庫約束判定。在滿足生態(tài)流量辨識前提下,在不考慮水量傳播損失的條件下,計算滿足條件的水量平衡方程。同時,判斷相應(yīng)流量下石橋水電站水庫水位是否滿足約束條件,從石橋水電站上游到下游尾水河道依次對水庫水位進行辨識。如果水庫水位超過最大值或低于最小值,則返回步驟一重新確定下泄流量。③電站出力約束判定。采用同時滿足生態(tài)流量判定、水庫約束判定情況下的下泄流量及相應(yīng)水位進行計算,得出水電站出力值。對水電站出力值進行判定,如出力值符合約束條件則對該組數(shù)據(jù)輸出,如出力值不符合約束判定,則返回初始狀態(tài)重新試算。
步驟三確定調(diào)度備選方案。將符合各約束條件的水電站保證出力及發(fā)電流量作為調(diào)度方案,并以求得的水電站調(diào)度方案來進行仿真模擬,重復(fù)以上步驟,輸出滿足約束條件的30組最優(yōu)備選方案。
基于以上的求解步驟,本研究在Visual Studio平臺下采用C#編寫模擬程序求解,并通過txt文本文件輸出備選方案。
當(dāng)枯水期來水量較小時,如電站以追求經(jīng)濟效益為主要目標,在這一時期將控制水電站的發(fā)電流量,減少低水能利用率情況下發(fā)電的時段。同時水庫進行蓄水,抬高發(fā)電水頭,直到水庫水位達到滿足其經(jīng)濟運行水位時進行發(fā)電。這種調(diào)度方案的優(yōu)點是可以提高總發(fā)電量,將水電站經(jīng)濟效益達到最大化。但是采用這種控制流量提高發(fā)電效益的做法會導(dǎo)致下游河道的流量減小。如果電站在追求經(jīng)濟目標的同時兼顧生態(tài)目標,則電站在枯水時期將增大下泄流量,保證下游河道的生態(tài)需要,采用這種方案整體發(fā)電量較上一種方案略低,但是對保護下游河道水生態(tài)環(huán)境有著積極的作用。
采用石橋水電站在典型年來水情況下枯水期最小出力為電站的保證出力,經(jīng)計算,輸出了30組在典型年來水條件下,石橋水電站優(yōu)化運行調(diào)度方案。將輸出數(shù)據(jù)整理后見表1。
表1 石橋水電站典型年來水條件下優(yōu)化運行調(diào)度方案
表1共列出了滿足各項約束條件的30組優(yōu)化運行調(diào)度方案,由于追求的主要目標不同,針對石橋水電站的優(yōu)化調(diào)度方案仍有差別。綜合表1中30組優(yōu)化運行調(diào)度方案,分別選取1號生態(tài)效益最大化方案、15號兼顧經(jīng)濟效益及生態(tài)效益方案、30號經(jīng)濟效益最大化調(diào)度方案作為石橋水電站優(yōu)化后運行調(diào)度的典型方案進行對比分析。石橋水電站典型年來水情況下滿足各約束條件的出力過程見圖1,石橋水電站典型年來水情況下滿足各約束條件的發(fā)電量見圖2,石橋水電站典型年來水情況下各月發(fā)電流量見圖3。
圖2 石橋水電站典型年來水條件下各月發(fā)電量
圖3 石橋水電站典型年來水情況下各月發(fā)電流量
從圖1、2、3可以看出,各運行方案在7—12月份水電站的發(fā)電出力情況、發(fā)電量及電站的發(fā)電流量基本相同,差別不大。各項指標的不同之處主要體現(xiàn)在1—6月份這段時間,即整個來水量相對較小的時期。而水電站的最小出力一般也多出現(xiàn)在枯水期。
圖1 石橋水電站出力過程
第1號方案,石橋水電站在實際運行過程中年發(fā)電量為1 996.8萬kW·h,保證出力為1 358 kW。電站全年保證出力略大于其他2組方案,并且電站全年發(fā)電過程出力波動幅度較小,能夠保證區(qū)域電網(wǎng)安全。但是發(fā)電量相對最小,雖然該種方案能夠保證生態(tài)流量最大,但是未能兼顧電站經(jīng)濟效益。
第30號方案,石橋水電站年發(fā)電量為2 044.8萬kW·h,保證出力為1 305.2 kW。全年發(fā)電量略大于其他2組方案,發(fā)電效益最大,但是該種方案下水電站保證出力小于其他2組方案,年利用小時數(shù)較少,且全年發(fā)電過程出力波動幅度較大,不利于區(qū)域電網(wǎng)的安全性,未能兼顧生態(tài)效益及電網(wǎng)安全需要。
第15號方案,石橋水電站年發(fā)電量為2 023.2萬kW·h,保證出力為1 336 kW。在電站出力、發(fā)電量等方面,指標均處于另2種方案之間,同時全年發(fā)電過程出力平穩(wěn),符合區(qū)域電網(wǎng)安全需要。兼顧了經(jīng)濟目標與生態(tài)目標,二者協(xié)調(diào)統(tǒng)一。
由以上3種典型優(yōu)化方案的對比分析可以看出,在滿足石橋水電站下游河道生態(tài)流量及其他約束條件的情況下,可以有針對性地進行調(diào)度方案的選擇。如果地區(qū)水生態(tài)環(huán)境基礎(chǔ)較為薄弱且電網(wǎng)公司對電站出力平穩(wěn)程度有特殊要求情況下,可以采用以滿足電站保證出力為主的調(diào)度方案,在能夠保證較大生態(tài)流量和利用小時數(shù)的情況下同時能夠滿足發(fā)電出力的平穩(wěn)性,雖然該種方案年發(fā)電量稍微有所降低,但是其保證出力有所提高,對下游河道的生態(tài)環(huán)境保護作用最好,如果考慮生態(tài)效益比重較大情況下可選擇第1號方案;當(dāng)?shù)貐^(qū)水生態(tài)環(huán)境基礎(chǔ)情況處于中等水平,且地區(qū)電網(wǎng)公司對電站出力過程沒有明確要求的情況下,綜合考慮發(fā)電效益及水生態(tài)環(huán)境保護可以優(yōu)先考慮第15號方案;如地區(qū)水生態(tài)環(huán)境基礎(chǔ)較好,電站可以采用以經(jīng)濟效益為主的調(diào)度方案,即以發(fā)電量為主要目標,可以考慮選擇第30號方案進行發(fā)電。
由于愛河流域水能資源豐富,地表水資源量較大,水生態(tài)環(huán)境基礎(chǔ)較好,且遼東山區(qū)供電網(wǎng)絡(luò)安全等級有限,所以綜合研究區(qū)的具體情況,經(jīng)分析推薦第15號方案為石橋水電站在典型年來水情況下合理運行調(diào)度方案。
a)建立的石橋水電站水能資源優(yōu)化開發(fā)模型,以水電站發(fā)電量最大及保證出力最大為目標建立了多目標函數(shù),以水電站下游生態(tài)流量為主要約束,指標獲取相對方便,直觀性、可操作性較強,與以往單目標調(diào)度運行方式相比,實用性、針對性、綜合性更加突出,對保護水電站下游區(qū)水生態(tài)環(huán)境并充分利用水能資源有著重要意義。
b)采用數(shù)學(xué)模擬法對模型進行求解,可以在短時間內(nèi)處理大量數(shù)據(jù),同時,數(shù)學(xué)模擬法具有直觀性較強的特點,求解模型過程可以省去很多的復(fù)雜的程序,實際應(yīng)用過程高效、準確。
c)石橋水電站工程是草河河口以下河段水電開發(fā)的第一級水電站工程,在流域水能資源開發(fā)序列中處于比較重要的地位,采用第15號方案,在保證發(fā)電效益的同時,滿足豐水、枯水不同時期的河道生態(tài)需水,避免了由于不合理調(diào)度導(dǎo)致的下游河段出現(xiàn)減脫水情況,有力保障了河段水生態(tài)環(huán)境,符合生態(tài)發(fā)展理念要求。