牟曉春, 王佳佳, 孫守鑫, 張超, 吳雪瑞
(1.國電南瑞科技股份有限公司,江蘇 南京 210000;2. 河海大學 能源與電氣學院,江蘇 南京 210000;3.國網(wǎng)濰坊供電公司,山東 濰坊 310000; 4. 河南科技大學 電氣工程學院,河南 洛陽 471023)
以高效、經(jīng)濟和環(huán)保的可再生能源為主的分布式發(fā)電技術(shù)得到越來越多的應(yīng)用[1-3]。作為微電網(wǎng)的核心元件——并網(wǎng)變流器,其幾乎沒有轉(zhuǎn)動慣量,當分布式電源的滲透率越來越高時,導(dǎo)致電力系統(tǒng)中的轉(zhuǎn)動慣量和旋轉(zhuǎn)備用容量的相對減少,分布式電源和負載的功率波動,極易引起電網(wǎng)系統(tǒng)失去穩(wěn)定。由于缺少轉(zhuǎn)動慣量和阻尼的作用,在微電網(wǎng)模式轉(zhuǎn)換時,其輸出的電壓和頻率常為階躍式,對電網(wǎng)系統(tǒng)造成極大的沖擊,更有可能引起并網(wǎng)失敗。為此,將并網(wǎng)逆變器和同步發(fā)電機的特性結(jié)合起來,同時在分布式發(fā)電電源交流側(cè)配置儲能系統(tǒng),并配以基于同步發(fā)電機模型的并網(wǎng)逆變器控制算法,這就是虛擬同步發(fā)電機(Virtual Synchronous Generator,VSG)技術(shù)[4-6]。
文獻[7-8]介紹了逆變器在孤島模式下的下垂控制策略,在孤島運行時提供電壓和頻率的參考值。文獻[9-14]提出微電網(wǎng)在并網(wǎng)時采用PQ控制,孤島時采用V/f控制,并給出了在不同工況下的仿真波形,但是當微電網(wǎng)在兩種運行模式切換時,需要改變控制測略。
本文詳細地講述了虛擬同步機的算法、組成結(jié)構(gòu)以及理論基礎(chǔ),并用MATLAB軟件分本體模型、有功功率-頻率控制器以及勵磁控制器三部分建立了VSG仿真模型?;赩SG的風儲一體發(fā)電系統(tǒng),進行了并網(wǎng)、孤島、并網(wǎng)孤島模式轉(zhuǎn)換及孤島下切增負荷的仿真驗證。仿真結(jié)果表明,所設(shè)計的系統(tǒng)能夠根據(jù)電網(wǎng)功率、電壓和頻率的變化,進行實時有效調(diào)節(jié),整個過程無需改變控制策略,保證了系統(tǒng)的穩(wěn)定性和可靠性。
電能在實際中的主要儲存方式有三大類型,分別是:物理儲能、電磁儲能和電化學儲能。
物理儲能指利用動能或勢能等方式對電能進行存儲,主要形式有飛輪儲能、抽水蓄能和壓縮空氣儲能等。電磁儲能是指通過電磁能實現(xiàn)儲能,基本方式有超導(dǎo)儲能和超級電容儲能等。電化學儲能關(guān)鍵在于電池,它可以在很長時間內(nèi)穩(wěn)定保存化學能,也可在隨時釋放電能。電池分為很多種類,不同的電池因為成分不同所產(chǎn)生的反應(yīng)也不一樣。下面介紹鋰電池的常見等效模型。
圖1 等效電路模型
1)理想模型
圖1(a)所示為理想等效電路模型。該鋰電池模型電路中包含一個電壓源和一個電阻,并且認為電阻的阻值及電壓源的電壓在工作過程中保持不變,此模型只能應(yīng)用在假定電池無限的仿真中。
2)通用等效電路模型
圖1(b)所示為通用等效電路模型。該電路中包含一個受控電壓源Eb和電池內(nèi)阻Rb。
忽略電池的自放電與記憶特性以及溫度的影響,默認電池的容量和內(nèi)阻為定值,受控源Eb的表達式為:
其中:
(1)
式中:Q為電池容量;Eb,E0分別為電池的空載電壓和電池恒定電壓;K為極化電壓;Efull,Eexp,Enom分別為完全充電電壓,末端點的電壓和標稱點的電壓;Qexp,Qnom分別為末端點的容量和標稱點的容量;A、B、K均由電池的電壓—容量放電曲線中相關(guān)的點計算得到。
為了避開復(fù)雜的磁場效應(yīng)以及體現(xiàn)同步發(fā)電機的基本特性,采用同步發(fā)電機的二階模型,其中包括同步發(fā)電機定子電壓方程以及轉(zhuǎn)子的機械方程。
(2)
(3)
式中:E0為勵磁電動勢;U、I、Ra、Xs分別為電樞端電壓、電樞電流、電樞電阻和電樞電感;J為轉(zhuǎn)動慣量;Ω為機械角速度;MT為機械轉(zhuǎn)矩;Me為電磁轉(zhuǎn)矩。
對式(2)、式(3)進行推導(dǎo)與變換,得到以電角度形式所表示的轉(zhuǎn)子機械運動方程:
(4)
式中:Ω為電角速度;ωN為同步電角速度;PT為機械功率;Pe為電磁功率;θ為電角度。其算法實現(xiàn)結(jié)構(gòu)如圖2所示,控制結(jié)構(gòu)流程如圖3所示。
圖2 虛擬同步機算法實現(xiàn)結(jié)構(gòu)圖
圖3 虛擬同步機算法控制流程圖
圖3的具體推導(dǎo)過程如下:
(1)先由轉(zhuǎn)子的機械方程計算出系統(tǒng)的電角速度ω,然后積分得到相位角。假定相角初始值為0,則A相相位初始值為ωt,代入正弦sin函數(shù),乘以由勵磁控制器得到的勵磁電動勢幅值Emag,則得到A相的勵磁電動勢向量。分別在相角上減去120°和240°就可以得到B相與C相的勵磁電動勢向量。三相疊加就得到三相合成勵磁電動勢。
(2)定子壓降的生成過程:設(shè)虛擬電阻和同步電阻抗值為Ra+jXs,電樞電流為儲能電池逆變器輸出電流Ia的實際采樣值,采用公式Ia=Ia∠Ψi將電流分解為幅值和相角;將電樞電阻作為實部,同步電抗作為虛部,合成向量,再轉(zhuǎn)換為幅值與相角形式;最后將幅值相乘相角相加,合成定子線路的電壓降。
(3)儲能單元輸出端電壓的生成:將之前的運算結(jié)果代入式(2)即可得到儲能單元輸出端電壓向量U。將之與實際端電壓值進行比較,經(jīng)PI控制器輸出到SVPWM調(diào)制生成控制變流器的脈沖觸發(fā)信號,完成整個本體模型算法。
2.2.1功頻控制器設(shè)計
發(fā)電機組發(fā)出的有功功率P1與電網(wǎng)頻率f1的相互關(guān)聯(lián)性能夠用功頻曲線的傾斜度來展現(xiàn),即:
(5)
式中:R為發(fā)電機組的調(diào)差系數(shù),反應(yīng)功率與頻率相反的變化趨勢;PN為額定有功功率;fN為額定頻率;ΔP為發(fā)電機組輸出有功功率的變化量;Δf為對應(yīng)頻率的變化量;ΔP*為有功功率變化量標么值;Δf*為頻率變化量的標么值;R*為發(fā)電機組頻率調(diào)差系數(shù)的標么值。
為了使在實際運行中的同步發(fā)電機組能恰當?shù)胤謸β首兓梢宰尣⒃谝黄鸬耐桨l(fā)電機的R*相等,一起參與一次調(diào)頻,從而維持電網(wǎng)頻率穩(wěn)定。例如,兩臺發(fā)電機并聯(lián)運行時,可得:
(6)
式中:ΔP1*和ΔP2*為兩臺發(fā)電機有功功率增加量的標么值;R1*和R2*為兩臺發(fā)電機組頻率調(diào)差系數(shù)的標么值。發(fā)電機組間的頻率調(diào)差系數(shù)與機組的有功功率分配成反比。
2.2.2勵磁控制器設(shè)計
定義電壓調(diào)差系數(shù)δ:
(7)
式中:QN為發(fā)電機組的額定無功功率;VN為額定電壓;Q1為發(fā)電機組實際的輸出無功功率;V1為實際電壓;ΔQ為發(fā)電機組輸出無功功率的變化量;ΔV為對應(yīng)電壓的變化量;ΔQ*為無功功率變化量的標么值;ΔV*為電壓變化量的標么值;δ*為電壓調(diào)差系數(shù)的標幺值。
當系統(tǒng)電壓下降時,發(fā)電機無功電流增大,具有下垂外特性,并聯(lián)的發(fā)電機合理分擔無功功率負荷。例如,兩臺并聯(lián)的發(fā)電機,當發(fā)電機的端電壓因為無功負荷增加而下降時,兩臺機承擔無功變化量為:
(8)
式中:ΔQ1*和ΔQ2*為兩臺發(fā)電機無功功率增加量的標么值;δ1*和δ2*為兩臺發(fā)電機組電壓調(diào)差系數(shù)的標么值。因此,在無功功率波動時,兩臺并聯(lián)發(fā)電機的無功分配與電壓調(diào)差系數(shù)大小成反比。假設(shè)并聯(lián)運行的每臺發(fā)電機的電壓調(diào)差系數(shù)相等,各發(fā)電機就會依容量大小自動分擔負荷功率的波動值。在真實的電網(wǎng)環(huán)境中,電壓調(diào)差系數(shù)通常取3%~5%。
圖4為風儲一體發(fā)電系統(tǒng)并網(wǎng)模型圖。
圖4 風儲一體發(fā)電系統(tǒng)并網(wǎng)模型圖
模型的仿真運行模式如表1所示。
表1 系統(tǒng)運行模式與負荷狀態(tài)表
3.2.1鋰電池性能
電池荷電狀態(tài)特性說明:電池荷電狀態(tài)(State of Charge,SOC)定義為:在規(guī)定放電倍率的情況下,電池中剩余電量與在相同條件下電池的額定容量的比值,即:
(9)
圖5 鋰電池電壓曲線圖
圖6 鋰電池電流曲線圖
圖7 鋰電池SOC曲線圖
式中:Qrem為電池的剩余容量;Qrat為電池的額定容量;當電池充滿電時SOC=100%;電池放完電時SOC=0%。
由圖5~圖7可知:2.5~4 s、6~7 s、12~14 s和17~20 s,鋰電池吸收功率,SOC值平穩(wěn)上升;10~11.5 s、14~17 s,鋰電池發(fā)出功率,SOC值平穩(wěn)下降。由波形可見:鋰電池放電曲線較平穩(wěn),且能保持相對穩(wěn)定的電壓水平,與鋰電池的充放電特性相符。
3.2.2風力發(fā)電功率平抑
圖8、圖9中,P2、P1和PT分別為風力發(fā)電輸出有功功率、儲能電池輸出(吸收)有功功率和風儲一體發(fā)電系統(tǒng)輸出功率。由圖可知,初始狀態(tài)下,風力發(fā)電系統(tǒng)的電能全部向負荷和電網(wǎng)輸送,儲能系統(tǒng)處于備用。
由圖8可知,儲能電池發(fā)出(吸收)功率的變化根據(jù)風力發(fā)電單元輸出功率的變化而變化。當風力發(fā)電單元輸出有功大于1.2 MW時,儲能電池會實時吸收多余的有功功率;當風力發(fā)電單元輸出有功小于1.2 MW時,儲能電池實時發(fā)出差額的有功功率以彌補功率的不足,從而平抑了由于風速的改變而造成的輸出功率波動。
圖8 儲能與風力發(fā)電單元的輸出功率
圖9 風儲一體發(fā)電系統(tǒng)輸出功率
3.2.3有功功率與頻率的仿真結(jié)果
圖10和圖11中:P1、P2、Pgrid分別為負荷1、2吸收的有功功率和無窮大電網(wǎng)提供的有功功率。由圖中可知:0~12 s,風儲系統(tǒng)與無窮大電網(wǎng)共同承擔系統(tǒng)中負荷1的有功功率;12 s時運行模式由并網(wǎng)變?yōu)楣聧u,由風儲系統(tǒng)承單獨擔負荷的全部有功功率;12~14 s,儲能電池吸收有功功率與0~12 s無窮大電網(wǎng)吸收的有功功率基本相等;14~17 s,儲能系統(tǒng)發(fā)出0.5 MW,以滿足14 s時增投的負荷2所導(dǎo)致的有功功率需求的增加;17 s,由于負荷2被切除,儲能系統(tǒng)迅速由發(fā)出0.5 MW轉(zhuǎn)變?yōu)槲?.075 MW。
圖10 負荷1、2上消耗的有功功率
圖11 無窮大電網(wǎng)提供的有功功率
圖12 系統(tǒng)頻率
如圖12所示,f為系統(tǒng)頻率。由于VSG能平抑功率,所以f在0~12 s內(nèi)始終穩(wěn)定于50 Hz。在12 s時由于運行狀態(tài)的改變以及一次調(diào)頻的作用,頻率升高并穩(wěn)定在50.06 Hz左右。在14 s時,在運行模式不變的情況下,增投負荷2,頻率相應(yīng)下降并且穩(wěn)定在49.6 Hz左右。在17 s時,在孤島的情況下,切除增投的負荷2,頻率經(jīng)過短時振蕩即回到并且穩(wěn)定在50.06 Hz左右。
3.2.4無功功率與電壓的仿真結(jié)果
圖13中,Q1、Q2分別為負荷1、2吸收的無功功率;圖14中Qgrid為電網(wǎng)提供的無功功率;圖15為并網(wǎng)點電壓波形。
圖13 負荷1、2上消耗的無功功率
圖14 無窮大電網(wǎng)提供的無功功率
圖15 風儲一體系統(tǒng)的并網(wǎng)相電壓幅值
由圖13可知:0~12 s,系統(tǒng)所帶的負荷一直恒定。風儲一體系統(tǒng)與無窮大電網(wǎng)共同向負荷1提供無功功率。在12 s時由于運行狀態(tài)的改變,風儲系統(tǒng)單獨承擔系統(tǒng)中負荷1的無功功率,儲能系統(tǒng)發(fā)出約80 kvar的無功功率給負荷1,使系統(tǒng)的端電壓和負荷上的電壓基本穩(wěn)定。
在14 s時,增投負荷,從Q1波形圖可以看出:儲能單元立刻增發(fā)約20 kvar的無功功率以滿足因負荷2引發(fā)的無功功率需求的增加。系統(tǒng)中的無功功率依靠增發(fā)的無功功率仍然維持平衡,如圖14所示。即使輸電線上無功功率的增多使線路壓降的增大,導(dǎo)致負荷上的電壓有了輕微的下降,但由電壓波形圖可知,無論是風儲一體系統(tǒng)的端電壓還是負荷1、負荷2上的電壓均能夠維持穩(wěn)定在額定值附近且波動不超過1%。
在17 s時,負荷2切除,儲能單元發(fā)出的無功功率由100 kvar重新回到80 kvar,系統(tǒng)電壓依舊保持穩(wěn)定。
本文研究了傳統(tǒng)同步發(fā)電機的特性,提出將VSG作為儲能系統(tǒng)的控制方法,掌握了VSG的算法和模型之后,最終搭建了以鋰電池儲能為核心的風儲一體發(fā)電系統(tǒng),并在運行模式變化和負荷狀態(tài)改變等情況下對其進行了仿真驗證。最終,通過仿真結(jié)果證明,VSG技術(shù)可以實現(xiàn)儲能輸出功率根據(jù)電網(wǎng)有功功率、無功功率的波動以及系統(tǒng)電壓、頻率的變化進行實時地調(diào)節(jié),并讓系統(tǒng)的電壓、頻率保持在標準允許的范圍內(nèi)。