駱 釗, 高培淇, 和 婧, 劉可真, 梁俊宇, 董 敏, 劉 果
(1.昆明理工大學(xué) 電力工程學(xué)院,云南 昆明 650500;2.云南電網(wǎng)有限責(zé)任公司電力科學(xué)研究院,云南 昆明 650000)
抽水蓄能電站(PPS)作為成熟的儲能技術(shù),具有容量大、經(jīng)濟性好、環(huán)保清潔等優(yōu)勢,對于我國“十四五”規(guī)劃及未來電力系統(tǒng)的發(fā)展至關(guān)重要。通過抽水蓄能電站,可充分發(fā)揮傳統(tǒng)機組的調(diào)頻潛力,提高常規(guī)機組的調(diào)頻性能[1-2],有效解決新能源規(guī)?;⒕W(wǎng)帶來的電網(wǎng)頻率波動問題[3]。然而PPS在夏季頻繁投切會導(dǎo)致電網(wǎng)間聯(lián)絡(luò)線波動增大,影響線路的輸電能力。
在此背景下,可引入電化學(xué)儲能(電池儲能系統(tǒng),BESS)參與調(diào)峰、調(diào)頻以平滑可再生能源出力的任務(wù)。BESS憑借快速響應(yīng)和精確的跟蹤特性來改善調(diào)頻效果,可有效彌補PSS在調(diào)頻上的缺陷[4]。目前BESS在提高電力系統(tǒng)的穩(wěn)定性以及促進可再生能源消納等方面有著廣闊的應(yīng)用潛力,將成為新一代智能電網(wǎng)的重要組成部分[5-6]。
目前投運的抽水蓄能電站與電化學(xué)儲能系統(tǒng)協(xié)同運行時存在功率備用上限、荷電狀態(tài)(SOC)限制及其他種種約束條件。傳統(tǒng)PID控制處理約束條件效能較差,難以同時應(yīng)對這些非線性約束[7]。自抗擾控制(ADRC)是一種非線性控制算法,能夠在不依賴于系統(tǒng)精確模型的基礎(chǔ)上自動補償對象模型的內(nèi)外干擾,實現(xiàn)對不確定、強耦合系統(tǒng)的有效控制,解決非線性約束問題[8-9]。
然而,現(xiàn)有研究僅僅考慮抽水蓄能或儲能系統(tǒng)單獨參與調(diào)頻,并未將兩者整合起來,發(fā)揮其協(xié)同優(yōu)勢。盡管PPS與BESS在功率等級、響應(yīng)時間等方面存在著一定差異[7],但BESS聯(lián)合PPS參與負(fù)荷頻率控制(LFC)的優(yōu)勢有以下幾點:(1)響應(yīng)速度快,可減少功率儲備裕度;(2)控制精準(zhǔn),可在任何功率點穩(wěn)定輸出;(3)雙向調(diào)節(jié);(4)減少PSS機組旋轉(zhuǎn)備用容量。
基于以上分析,本文結(jié)合BESS與PSS調(diào)頻性能,提出一種基于ADRC技術(shù)的抽/儲聯(lián)合LFC策略,并引入需求響應(yīng)(DR)參與調(diào)頻。與傳統(tǒng)PID控制器和分?jǐn)?shù)階PID控制器進行對比,設(shè)置不同情景,對調(diào)節(jié)時間、頻率偏移、振蕩幅值三個方面進行仿真。仿真結(jié)果驗證了本文所提控制方法是正確且有效的。
機組模型包括以汽輪機和水輪機為主的原動機及其調(diào)速器[10]。本文汽輪機選用再熱式汽輪機,模型為
(1)
式中:Tt為蒸汽室時間常數(shù);Kr為高壓缸段蒸汽與總功率之比;Tr為汽輪機時間常數(shù);s為Laplace算子。
汽輪機調(diào)速器模型為
(2)
式中:Tg為調(diào)速器慣性環(huán)節(jié)的時間常數(shù)。
水輪機模型為
(3)
式中:Tw為水起動時間。
水輪機調(diào)速器模型為
(4)
式中:R為水輪機機組調(diào)差系數(shù);f為系統(tǒng)額定頻率;Kp、Ki、Kd分別為比例環(huán)節(jié)、積分環(huán)節(jié)和微分環(huán)節(jié)的增益。
電網(wǎng)的二次調(diào)頻方式一般采用自動發(fā)電控制(AGC),根據(jù)受到擾動時機組的區(qū)域控制誤差(ACE)信號變化,調(diào)整該區(qū)域的發(fā)電機組輸出。按照不同的ACE策略,AGC的調(diào)控方式分為3種[11]。本文采用聯(lián)絡(luò)線偏差控制(TBC)模式,其表達式為
eACE=ΔPtie+β×Δf
(5)
式中:eACE為區(qū)域控制誤差;Δf為頻率偏移量;ΔPtie為聯(lián)絡(luò)線交換功率偏移量。
系統(tǒng)的頻率偏差系數(shù)β定義為
βi=Di+1/Ri
(6)
為了使Δf和ACE為零,系統(tǒng)控制信號以ACE為輸入,Ri為調(diào)節(jié)系數(shù),負(fù)載阻尼系數(shù)Di取值為8.33×10-3MW/Hz。
(1) 調(diào)速器死區(qū)。由于調(diào)速器的遲滯性和非線性,其死區(qū)一般用間隙特性環(huán)節(jié)來描述。假設(shè)輸入輸出關(guān)系為[12]
(7)
式中:X為非線性元件的輸入;Y為非線性元件的輸出。
當(dāng)函數(shù)Y呈正弦關(guān)系時,有:
X=Asinω0t
(8)
式中:A為正弦輸入信號幅值;ω0為正弦輸入信號頻率。
(9)
其中:各項系數(shù)為
(10)
將式(10)代入式(9)并經(jīng)Laplace變換,結(jié)合式(2),可得到線性化后的計及死區(qū)的汽輪機調(diào)速器傳遞函數(shù)為
(11)
(2) 發(fā)電速度約束(GRC)定義為爬坡速率限制,指機組出力受到某個最大變化率的約束。其限位裝置數(shù)值設(shè)置為1.7×10-3MW/s[9]。
忽略網(wǎng)絡(luò)線損,兩區(qū)域間的聯(lián)絡(luò)線功率流動表示為[13]
(12)
兩區(qū)域聯(lián)絡(luò)線功率微增量模型為
ΔPtieij=
(13)
式中:ΔPtieij為區(qū)域i、j流動功率微增量;Δfi和Δfj分別為i區(qū)域和j區(qū)域的頻率偏差。ΔPtieij經(jīng)Laplace變換后有:
(14)
式中:Tij為聯(lián)絡(luò)線同步系數(shù),可表示為
(15)
忽略區(qū)域互聯(lián)電網(wǎng)中各機組之間的差異性,對電網(wǎng)進行精簡建模,用聯(lián)絡(luò)線將兩區(qū)域連接起來,每個區(qū)域仍以再熱式汽輪機組為基礎(chǔ),再考慮非線性因素,得到兩區(qū)域互聯(lián)的LFC模型,如圖1所示。
圖1中,ACE為輸入,Δf為系統(tǒng)輸出;ΔPDi,ΔPDj為不同區(qū)域施加的負(fù)荷擾動;B為頻率偏差系數(shù);aij為同步功率系數(shù),有:
(16)
Tij=0.545,aij=-1,B=0.425
(17)
電網(wǎng)中發(fā)電機與系統(tǒng)連接的傳遞函數(shù)關(guān)系為
(18)
式中:Kpi1/D;Tri2H/fD;ΔPg為發(fā)電機組的輸出功率;D為負(fù)荷阻尼系數(shù);ΔPd為負(fù)荷狀態(tài)下功率變化量;H為再熱式汽輪機組慣性常數(shù)。
圖1 抽/儲聯(lián)合系統(tǒng)LFC模型
本文提出的抽/儲聯(lián)合LFC模型以兩區(qū)域再熱式汽輪機組LFC模型為基礎(chǔ)。負(fù)荷低谷期時,電化學(xué)儲能充足,PPS依賴盈余電力從下游水庫向上抽水,處于抽水狀態(tài),通過勢能的形式將能量儲存起來。負(fù)荷高峰期時,電化學(xué)儲能總量不足,蓄能電站為系統(tǒng)提供輸出,處于發(fā)電狀態(tài),上游水庫的水通過水輪機流向下游水庫發(fā)電。本文將PPS的發(fā)電工況模型綜合簡化為一臺水輪機機組,模擬作為一個水電站對外供電。在抽水工況下,相當(dāng)于一個負(fù)荷。在此基礎(chǔ)上建立了兩段式PPS傳遞函數(shù)模型。其仿真結(jié)構(gòu)如圖1虛線部分發(fā)電工況和抽水工況所示。
兩段式PPS傳遞函數(shù)模型如下:
(19)
式中:ΔPpd為抽水工況停止運行時輸送到電網(wǎng)的功率;Gt(s)為水輪機的傳遞函數(shù)模型,如式(3)所示;Gd(s)為發(fā)電工況下水輪機調(diào)速器的傳遞函數(shù)模型,如式(4)所示。圖1中η為PPS在抽水工況下機組工作的數(shù)量,η={1,2,…,n}。
按照基爾霍夫定律簡化得到如圖2所示的電化學(xué)儲能等效電源模型[4]。圖2中,Rc表示能量轉(zhuǎn)換系統(tǒng)與儲能電池單元模型的連接阻抗;n為串聯(lián)電池單體數(shù)目;m為儲能電池單元模型中由電池單體串聯(lián)形成的子系統(tǒng)并聯(lián)數(shù)目;k為所有單元模型數(shù)目;Rseries為電池內(nèi)阻;Ct、Rt分別為過電容和電阻;Voc為開路電壓。
圖2 電化學(xué)儲能電源等效模型
根據(jù)戴維寧等效電路模型(II型),將圖2的等效模型圖轉(zhuǎn)化為傳遞函數(shù)的形式[14],得到對應(yīng)的傳遞函數(shù)數(shù)學(xué)模型表達式:
(20)
(21)
(22)
(23)
(24)
ΔPb(s)=Δib(s)Vb(s)=
ib(s)[Voc(s)+ΔVt(s)+
ΔVseries(s)+ΔVc(s)]
(25)
式中:control(s)為控制信號,電池內(nèi)部以電流作為響應(yīng)變量;Tb為一階慣性環(huán)節(jié)時間常數(shù);開路電壓Voc為SOC的函數(shù),常用百分?jǐn)?shù)表示;C0為電池初始容量;Cp為電池額定容量;ΔVc(s)為連接阻抗電壓增量;ΔVb(s)為儲能電池前端電壓值;ΔVseries(s)為電池內(nèi)阻電壓增量;ΔVt(s)為暫態(tài)電壓增量;ΔPb(s)為實際輸出功率偏差。電池SOC表達式如下[15]:
(26)
式中:S(t)和S(t0)分別為電池在t時刻和初始時刻t0的SOC;ηbat為電池充放電效率,代表充放電循環(huán)中內(nèi)部耗散的電量。
電化學(xué)儲能電源傳遞函數(shù)模型如圖3所示。
圖3 電化學(xué)儲能電源傳遞函數(shù)模型
該模型具體實現(xiàn)步驟如下:在接收外部控制信號后,電池內(nèi)部以電流作為響應(yīng)變量,電流經(jīng)過積分環(huán)節(jié)后與初始電池電量作差并除以電池額定容量得到儲能電池的SOC值,然后根據(jù)SOC值與儲能電池開路電壓之間的函數(shù)關(guān)系求得開路電壓Voc值,此外,電流經(jīng)過電池的內(nèi)部電阻與過電壓電阻時會產(chǎn)生相應(yīng)的電壓增量,電壓增量與開路電壓Voc進行疊加后得到儲能電池當(dāng)前端電壓值Vb,Vb同Δib的乘積即為當(dāng)前儲能電池的功率輸出量。
在本文所描述的短時調(diào)頻場景中,不考慮ACE死區(qū)限額、最大負(fù)荷調(diào)整限幅和機組容量限額,抽/儲聯(lián)合參與電網(wǎng)二次調(diào)頻方法如圖4所示。
圖4 抽/儲聯(lián)合系統(tǒng)二次調(diào)頻方法
當(dāng)擾動發(fā)生時,負(fù)荷需求增大,負(fù)荷頻率特性曲線由L1(Δf)平移至L2(Δf),系統(tǒng)穩(wěn)定運行點A隨之移動至非穩(wěn)定點B,產(chǎn)生Δf1的頻率偏差。系統(tǒng)二次調(diào)頻傳統(tǒng)機組出力為ΔPG,使功率曲線G1(Δf)平移至G2(Δf),系統(tǒng)運行點移至點C,頻率偏差恢復(fù)至Δf2。此時PSS與BESS工作,出力分別為PE和PPE,以系統(tǒng)ACE為控制變量進行功率分配及補償,參與系統(tǒng)整體調(diào)節(jié),保持電網(wǎng)頻率穩(wěn)定。