何 滔，汪廣明，汪文元，熊 璽

(國(guó)家能源大渡河沙坪發(fā)電有限公司，四川 樂山 614300)

充足的調(diào)峰能力和快速的響應(yīng)能力是保持電網(wǎng)穩(wěn)定對(duì)電力系統(tǒng)電源端提出的十分重要的要求之一。梯級(jí)水電短期內(nèi)分級(jí)調(diào)能效果受到水庫(kù)耦合性、水流時(shí)滯性、機(jī)組爬坡率、水文氣象、電站調(diào)度等多種原因的影響，是一類典型的大規(guī)模、強(qiáng)耦合、多約束、動(dòng)態(tài)、離散的復(fù)雜非線性優(yōu)化問題,也是水電能源優(yōu)化運(yùn)行領(lǐng)域的研究熱點(diǎn)之一[1-6]。在梯級(jí)水電聯(lián)合調(diào)度過程中，往往易忽略控制性龍頭水庫(kù)的負(fù)荷波動(dòng)情況，而下游電站的不規(guī)律調(diào)節(jié)和被動(dòng)跟隨與其相關(guān)性顯著。采用斯皮爾曼系數(shù)的加權(quán)樸素貝葉斯算法定量描述條件屬性與決策類之間的相關(guān)性，從而為目標(biāo)電站分解協(xié)調(diào)算法調(diào)控中各條件屬性賦予相應(yīng)權(quán)值，部分消除了耦合條件獨(dú)立假設(shè)對(duì)影響因子分類效果的影響。

1 水電站的不規(guī)律調(diào)節(jié)及其影響因子關(guān)聯(lián)性分析方法

采用斯皮爾曼系數(shù)的加權(quán)樸素貝葉斯分類算法來確定電站不規(guī)律調(diào)節(jié)與各影響因子之間的關(guān)聯(lián)性及關(guān)聯(lián)強(qiáng)度。其分類式可表達(dá)為：

式中：Cj為電站不規(guī)律調(diào)節(jié)影響因子數(shù)據(jù)集中m個(gè)類，j=1,2,3,…,m；P(Cj)表示Cj類發(fā)生的先驗(yàn)概率；wk為電站負(fù)荷區(qū)間、流量變化區(qū)間、庫(kù)區(qū)水位變化區(qū)間等條件屬性與上游電站負(fù)荷波動(dòng)耦合、站內(nèi)條件要求、自身庫(kù)容約束、控制水位要求、調(diào)度要求等決策類之間的斯皮爾曼系數(shù)ρk的權(quán)值。

2 水電集群聯(lián)合調(diào)能的必要性

大渡河流域上游往下游選取鄰近幾個(gè)梯級(jí)電站：瀑布溝、深溪溝、枕頭壩、沙南、龔嘴、銅街子，其中電站瀑布溝、深溪溝、枕頭壩已經(jīng)實(shí)現(xiàn)了三站聯(lián)調(diào)且瀑布溝是下游幾個(gè)電站的控制性水庫(kù)，龔嘴、銅街子也已實(shí)現(xiàn)了兩站聯(lián)調(diào)，而沙南電站由于2019年5月進(jìn)入全投，機(jī)組尚未經(jīng)歷完整檢修期的設(shè)備穩(wěn)定性檢驗(yàn)，設(shè)備定值、缺陷也尚需調(diào)整和治理，水庫(kù)調(diào)度規(guī)律仍需在實(shí)際投產(chǎn)發(fā)電的完整周期年中完善探索，因此沙南電站尚未納入聯(lián)調(diào)隊(duì)列。在各種因素的影響下，未納入集群調(diào)能電站局限性表現(xiàn)的相對(duì)明顯，具體表現(xiàn)如表1所示。

表1 大渡河流域梯級(jí)電站對(duì)應(yīng)月泄流設(shè)施動(dòng)作統(tǒng)計(jì)表

3 水電集群調(diào)能關(guān)鍵耦合因素分析

水電站控制的主要戰(zhàn)略對(duì)象是發(fā)電量與水資源，其控制主要是調(diào)節(jié)發(fā)電量與水資源的耦合動(dòng)態(tài)平衡關(guān)系。而影響此類耦合動(dòng)態(tài)平衡關(guān)系的幾個(gè)主要影響因素有：電站的固有屬性、耦合電站負(fù)荷波動(dòng)、計(jì)劃負(fù)荷與實(shí)時(shí)發(fā)電之間關(guān)系、自身內(nèi)部控制等。

3.1 電站固有屬性分析

選取流域幾個(gè)鄰近電站進(jìn)行固有屬性分析，如表2所示。

表2 大渡河流域梯級(jí)電站固有屬性表

從固有屬性可知，沙南電站在集群電站中具有最大耗水率、最小可調(diào)節(jié)庫(kù)容兩大典型特征，電站自身庫(kù)容小影響了其調(diào)節(jié)能力。選擇與沙南電站庫(kù)容相近的深溪溝，深溪溝調(diào)節(jié)庫(kù)容(848萬m3)約為沙南的1.44倍，單從調(diào)節(jié)時(shí)間裕度來說，深溪溝、沙南兩站差別不大，但從最終泄流設(shè)施動(dòng)作情況來看，深溪溝明顯少于沙南同期。

在研究對(duì)象樣本容量中隨機(jī)抽樣選取連續(xù)時(shí)間段瀑布溝、深溪溝、枕頭壩、沙南四站有功數(shù)據(jù)進(jìn)行聯(lián)合分析，得出集群調(diào)能電站瀑布溝、深溪溝、枕頭壩三站全站實(shí)時(shí)出力曲線跟隨性較好，三站負(fù)荷變化頻率、步調(diào)基本一致，三站期間泄流設(shè)施動(dòng)作總次數(shù)0次，充分體現(xiàn)了水電集群聯(lián)合調(diào)能的優(yōu)越性；期間沙南負(fù)荷波動(dòng)量較少，波動(dòng)幅度與波動(dòng)頻率低，閘門動(dòng)作次數(shù)21次，負(fù)荷調(diào)節(jié)更多是“跟趨勢(shì)”，步調(diào)難與上游電站保持一致。由于沙南電站現(xiàn)階段僅能依靠人工觀測(cè)上游三站負(fù)荷增減情況來向調(diào)度中心申請(qǐng)負(fù)荷，導(dǎo)致負(fù)荷難與上游三站步調(diào)一致，最終上游庫(kù)水位的控制主要依靠泄流設(shè)施進(jìn)行控制。因此調(diào)節(jié)庫(kù)容小不是泄流設(shè)施動(dòng)作頻繁的最主要影響因素。

3.2 電站負(fù)荷波動(dòng)耦合關(guān)系分析

為盡量減少沙南電站閘門動(dòng)作次數(shù)，降低閘門動(dòng)作頻率，認(rèn)為沙南電站閘門每次必要?jiǎng)幼鳛? m，則對(duì)應(yīng)期望流量變化范圍在198.11～216.05 m3/s。為便于計(jì)算，取變化流量為200 m3/s。

經(jīng)統(tǒng)計(jì)，選擇2018年1月至10月瀑布溝、深溪溝、枕頭壩三站對(duì)應(yīng)整點(diǎn)負(fù)荷及沙南泄流設(shè)施動(dòng)作為樣本進(jìn)行分析，可以發(fā)現(xiàn)沙南臨近上游三站在此時(shí)段中負(fù)荷波動(dòng)情況及波動(dòng)概率較大，如表3所示。

表3 沙南電站上游三電站對(duì)應(yīng)負(fù)荷波動(dòng)及波動(dòng)頻率表

從負(fù)荷波動(dòng)及其頻率中可以看出，上游三站負(fù)荷平均波動(dòng)概率為0.286 7，總波動(dòng)概率為0.284 7，偏差為0.002，在95%置信區(qū)間范圍內(nèi)。

分析上游三站與沙南電站泄流設(shè)施動(dòng)作相關(guān)性，記沙南泄流設(shè)施分為動(dòng)作與不動(dòng)作兩種情況。泄流設(shè)施動(dòng)作大于或等于1均記為1，所有泄流設(shè)施均不動(dòng)作記為0?？紤]到水流時(shí)間滯后性，分別取瀑布溝、深溪溝、枕頭壩、沙南(記為A、B、C、D)響應(yīng)時(shí)間間隔ΔtA,D、ΔtB,D、ΔtC,D為4、3、2 h，若上游對(duì)應(yīng)電站負(fù)荷波動(dòng)且沙南泄流設(shè)施動(dòng)作則記為1，否則記0，以沙南上游三站負(fù)荷波動(dòng)為自變量，沙南泄洪閘是否動(dòng)作為因變量，構(gòu)建基于斯皮爾曼系數(shù)的樸素貝葉斯分析模型，得到協(xié)同動(dòng)作樣本與協(xié)同動(dòng)作概率，結(jié)果如表4所示。

表4 沙南電站泄流設(shè)施動(dòng)作霍斯默-萊梅肖檢驗(yàn)列聯(lián)表 次

基于斯皮爾曼系數(shù)的加權(quán)樸素貝葉斯算法可定量描述條件屬性與決策類之間的相關(guān)性，消除條件獨(dú)立假設(shè)對(duì)耦合關(guān)聯(lián)效果的影響。它是衡量變量間的依賴性的非參數(shù)指標(biāo)，利用單調(diào)方程評(píng)價(jià)統(tǒng)計(jì)變量間的相關(guān)性[7-9]。

對(duì)該模型進(jìn)行擬合優(yōu)度檢驗(yàn)，檢驗(yàn)預(yù)測(cè)值和實(shí)際值之間差別的大小，進(jìn)行多次數(shù)據(jù)擬合試驗(yàn)，以判斷得出結(jié)論的可靠性。H-L檢驗(yàn)的隨機(jī)性表是比較觀測(cè)值與期望值，經(jīng)觀測(cè)，表中期望值與實(shí)測(cè)值誤差極小，表明該模型擬合度較高[10-12]。

通過以上數(shù)據(jù)分析，發(fā)現(xiàn)控制性瀑布溝負(fù)荷波動(dòng)對(duì)未納入集群聯(lián)合調(diào)能隊(duì)列的沙南泄流設(shè)施動(dòng)作總次數(shù)有決定性影響，沙南泄流設(shè)施動(dòng)作與緊鄰上游電站枕頭壩負(fù)荷波動(dòng)規(guī)律跟隨性更好，枕頭壩負(fù)荷波動(dòng)在短時(shí)間內(nèi)可能很大程度直接影響沙南泄流設(shè)施動(dòng)作情況，但在長(zhǎng)期規(guī)律來看沙南泄流設(shè)施動(dòng)作總次數(shù)更大程度受制于瀑布溝負(fù)荷波動(dòng)。如圖1所示。

圖1 沙南上游三電站負(fù)荷波動(dòng)與沙南閘門動(dòng)作表征圖

斯皮爾曼相關(guān)系數(shù)表明X(獨(dú)立變量)和Y(依賴變量)的相關(guān)方向。如果當(dāng)X增加時(shí)，Y趨向于增加，斯皮爾曼相關(guān)系數(shù)則為正。如果當(dāng)X增加時(shí)，Y趨向于減少，斯皮爾曼相關(guān)系數(shù)則為負(fù)[13-14]。斯皮爾曼相關(guān)系數(shù)為0表明當(dāng)X增加時(shí)Y沒有任何趨向性。當(dāng)X和Y越來越接近完全的單調(diào)相關(guān)時(shí)，斯皮爾曼相關(guān)系數(shù)會(huì)在絕對(duì)值上增加。當(dāng)X和Y完全單調(diào)相關(guān)時(shí)，斯皮爾曼相關(guān)系數(shù)的絕對(duì)值為1[15]。

對(duì)沙南泄流設(shè)施動(dòng)作次數(shù)與上游三站負(fù)荷波動(dòng)頻率進(jìn)行斯皮爾曼系數(shù)的加權(quán)樸素貝葉斯算法檢驗(yàn)并做相關(guān)性分析，結(jié)果如表5所示。

表5 四個(gè)電站斯皮爾曼相關(guān)系數(shù)及顯著性表

通過相關(guān)性可看出，從總量角度來看，瀑布溝負(fù)荷波動(dòng)與沙南泄流設(shè)施動(dòng)作相關(guān)系數(shù)為0.912，深溪溝負(fù)荷波動(dòng)與沙南泄流設(shè)施動(dòng)作相關(guān)系數(shù)為0.03，枕頭壩負(fù)荷波動(dòng)與沙南泄流設(shè)施動(dòng)作相關(guān)系數(shù)為0.042。

以上結(jié)果可以印證未納入集群聯(lián)調(diào)隊(duì)列電站泄流設(shè)施動(dòng)作頻率極大相關(guān)于上游控制性水庫(kù)泄流，周期響應(yīng)頻率相關(guān)于上游最鄰近電站泄流量的結(jié)論。

3.3 考慮上游控制性水庫(kù)負(fù)荷波動(dòng)頻率耦合因子的水電集群分解協(xié)調(diào)調(diào)控

水電集群調(diào)控分解協(xié)調(diào)思想是將水電集群大系統(tǒng)暫時(shí)割裂為若干獨(dú)立子系統(tǒng)，并在系統(tǒng)分解層逐步分別尋求局部最優(yōu)解；在各子系統(tǒng)最優(yōu)解或已有解的基礎(chǔ)上，恢復(fù)各子系統(tǒng)之間耦合關(guān)系，并在系統(tǒng)協(xié)調(diào)層通過耦合變量進(jìn)行更高級(jí)協(xié)調(diào)，達(dá)到全系統(tǒng)最優(yōu)。

模型以優(yōu)化調(diào)度模型為切入點(diǎn)，考慮到子系統(tǒng)內(nèi)瀑布溝、深溪溝、枕頭壩、沙南之間的聯(lián)系與制約，分解層分別考慮四站的負(fù)荷調(diào)整時(shí)間與調(diào)整量、泄流設(shè)施調(diào)整時(shí)間與調(diào)整量及四站機(jī)組和泄流設(shè)施的運(yùn)行狀態(tài)，系統(tǒng)協(xié)調(diào)層則考慮到發(fā)電量最大、泄流設(shè)施動(dòng)作最小及水位約束、生態(tài)流量約束、機(jī)組避開振動(dòng)區(qū)的目標(biāo)，應(yīng)用目標(biāo)協(xié)調(diào)法構(gòu)建Lagrange函數(shù)，將水量平衡約束和各站間流量耦合約束納入?yún)f(xié)調(diào)目標(biāo)。

整個(gè)過程中，充分考慮上游控制性水庫(kù)電站(分析模型中的瀑布溝電站)負(fù)荷波動(dòng)頻率耦合因子，與原僅關(guān)注上游物理位置臨近電站負(fù)荷變化(分析模型中的枕頭壩電站)的調(diào)度方案對(duì)比，沙南電站年發(fā)電量可增加1.385億kW·h，減少棄水量34.348億m3，泄流設(shè)施動(dòng)作次數(shù)減少62.21%，水庫(kù)耦合補(bǔ)償效益顯著。

4 結(jié) 語

1)控制性水庫(kù)下游未納入流域統(tǒng)調(diào)的電站負(fù)荷及水調(diào)跟隨性較差，豐水期導(dǎo)致自身泄流設(shè)施頻繁動(dòng)作，增加泄流設(shè)施及下游防汛潛在安全隱患。

2)在實(shí)際調(diào)控過程中，充分考慮上游控制性水庫(kù)電站負(fù)荷波動(dòng)頻幅，能一定程度上提升大渡河流域瀑布溝以下水電集群安全穩(wěn)定性與沙南電站負(fù)荷承擔(dān)能力，可對(duì)其它流域水電調(diào)控提供一定的參考。