李冬玉，婁 柯，尹 杰，鞏冠華

（安徽工程大學 電氣工程學院,安徽 蕪湖 241000）

近年來,化石燃料過度使用,導致環(huán)境與能源問題愈發(fā)嚴峻［1］。在此背景下,可再生能源發(fā)電如風力發(fā)電和光伏發(fā)電等備受青睞,包含多種分布式能源設(shè)備的微電網(wǎng)系統(tǒng)得到了快速發(fā)展［2］。隨著環(huán)境惡化與能源科技的進步,人們對微電網(wǎng)的經(jīng)濟環(huán)保性有了更高的要求。因此,研究一種算法來提高微電網(wǎng)的經(jīng)濟環(huán)保性是很有必要的。

微電網(wǎng)優(yōu)化調(diào)度問題,一直是國內(nèi)外學者研究的熱點。目前,微電網(wǎng)優(yōu)化的主要算法有粒子群優(yōu)化（PSO）算法［3-5］、遺傳優(yōu)化（GA）算法［6］和模擬退火（SA）算法［7］。PSO 算法易于實現(xiàn)且收斂速度快,但此算法受初始值影響大,可能會導致算法陷入局部最優(yōu)解,無法找到全局最優(yōu)解。GA算法具有很強的全局尋優(yōu)能力,可以有效避免算法陷入局部最優(yōu)解,但算法收斂速度慢,計算成本高。SA 算法具有很好的魯棒性,可以處理更高維度的尋優(yōu)問題,但其過于依賴超參數(shù)的選擇,容易陷入局部最優(yōu)解。

綜上所述,針對目前主流微電網(wǎng)優(yōu)化算法的種種缺陷。本文采用魯棒性強的原子搜索優(yōu)化算法作為理論基礎(chǔ),并引入自適應思想,加強了算法的收斂性與全局尋優(yōu)能力。最后,通過在孤島和并網(wǎng)運行模式下的對比算例仿真,證明了該算法在微電網(wǎng)優(yōu)化調(diào)度上的有效性和優(yōu)越性。

1 微電網(wǎng)系統(tǒng)

本文的微電網(wǎng)系統(tǒng)由光伏發(fā)電（PV）、風力發(fā)電（WT）、燃料電池（FC）、微型燃氣輪（MT）和蓄電池模塊（BT）構(gòu)成,如圖1所示。

圖1 微電網(wǎng)結(jié)構(gòu)圖

1.1 光伏發(fā)電

光伏發(fā)電的功率主要取決于太陽輻射和環(huán)境溫度,其發(fā)電功率預測公式PPV和維護成本公式CPV如下所示［8］：

式中,PSTC、TSTC、ISTC分別為標準測試下的光伏發(fā)電功率、電池板溫度和光照強度；IT、T分別為在t時刻下的光照強度和電池板溫度；λt為溫度系數(shù)；KPV,i為第i臺光伏微電源的維護系數(shù)。

1.2 風力發(fā)電機

風速的大小直接決定風力發(fā)電設(shè)備的發(fā)電功率,其發(fā)電功率PWT和維護成本CWT如下所示［9］：

式中,Prwt為t時刻風力發(fā)電的額定功率,v、vci、vco、vr分別為風力發(fā)電設(shè)備的實時風速、切入風速、切出風速和額定風速；KWT,i為風力微電源的維護系數(shù)。

1.3 微型燃氣輪機

微型燃氣輪機通過燃燒燃料來輸出電能,其發(fā)電成本包括維護成本和燃料費用。微型燃氣輪機的發(fā)電效率ηMT和發(fā)電成本CMT公式［10］為：

式中,PMT為t 時刻下微型燃氣輪的發(fā)電功率；LHV 為天然氣低熱值；Ci為第i 個微型燃氣輪設(shè)備的燃料單價；KMT,i為微型燃氣輪設(shè)備的維護系數(shù)。

1.4 燃料電池

燃料電池在微電網(wǎng)中的發(fā)電成本CFC同樣包括維護成本和燃料費用［11］：

式中,PFC、ηFC分別為t時刻下燃料電池的輸出效率和發(fā)電效率；KFC,i為燃料電池維護系數(shù)。

1.5 蓄電池

蓄電池在微電網(wǎng)中應用最為廣泛,其荷電狀態(tài)SOC和維護成本CBT［12］為：

式中,Pch、Pdis、ηch、ηdis分別為t 時刻下的充電功率、放電功率、充電效率和放電效率；KBT,i為蓄電池的維護系數(shù)。

2 微電網(wǎng)多目標優(yōu)化模型

2.1 目標函數(shù)

微電網(wǎng)優(yōu)化策略的主要思想是在滿足用電負荷需求的前提下,合理分配微電網(wǎng)的分布式能源,以實現(xiàn)經(jīng)濟效益最大化、環(huán)境污染最小化。本文的優(yōu)化調(diào)度策略的目標函數(shù)F由兩部分組成,分別是經(jīng)濟成本函數(shù)F1和環(huán)保成本函數(shù)F2：

2.1.1 經(jīng)濟性

為了實現(xiàn)微電網(wǎng)運行經(jīng)濟最優(yōu)化,本研究考慮的經(jīng)濟成本主要包括發(fā)電成本（維護成本和燃料費用）、設(shè)備折舊成本以及與外部電網(wǎng)的交互成本,如式（11）所示。

式中,CMC、CDE、Cgrid分別為微電網(wǎng)的發(fā)電成本、折舊成本以及與外部電網(wǎng)的交互成本。

1）發(fā)電成本

2）折舊成本

式中,α為折現(xiàn)率；n 為投資償還期；Ci為第i個微電源的安裝成本；fi為第i 臺微電源的容量因數(shù)。

3）與外部電網(wǎng)交互成本

式中,Pb、Ps、xb、xs分別為t時刻的購電效率、售電效率、購電單價和售電單價；Cb、Cs分別是t 時刻的購電和售電成本。

2.1.2 環(huán)保性

微型燃氣輪機和燃料電池作為微電網(wǎng)的發(fā)電單元,其在發(fā)電過程中不可避免地會產(chǎn)生污染氣體,主要成分為一氧化碳（CO）、二氧化碳（CO2）、二氧化硫（SO2）和氮氧化合物（NOx）。因此,微電網(wǎng)的環(huán)境保護成本主要為治理微電網(wǎng)系統(tǒng)產(chǎn)生污染氣體的成本,其公式為：

式中,K1、K2、K3、K4分別為CO、CO2、SO2和NOx的治理費用；E1、E2、E3、E4分別為t時刻CO、CO2、SO2和NOx的排放量。

2.2 約束條件

1）功率平衡約束

式中,PL、Pgrid分別為t 時刻微電網(wǎng)用電負荷和與外部電網(wǎng)的交互功率。

2）分布式電源輸出功率約束

式中,Pi、分別為各分布式電源的輸出功率、最大和最小輸出功率。

3）與外部電網(wǎng)的交互功率約束

3 自適應原子搜索算法

3.1 原子搜索算法

原子搜索優(yōu)化算法由Zhao［13］于2018年提出的一種新穎智能算法,通過原子之間的相互作用力和約束力,促使原子探索系統(tǒng)空間解決優(yōu)化問題。ASO 算法易實現(xiàn)、尋優(yōu)能力強、魯棒性好,因此被廣泛使用。本文以ASO 算法作為理論基礎(chǔ),引入自適應思想來加強算法的全局尋優(yōu)能力,以及收斂速度和精度。

ASO算法具體如下：

根據(jù)牛頓第二定律,計算第n 個原子的加速度：

式中,Fn、Gn、mn分別為原子n 的作用力、約束力和原子質(zhì)量；T、t分別為最大迭代次數(shù)和當前迭代次數(shù)；randi為［0,1］的隨機數(shù)；λ,ξ為深度權(quán)重和系數(shù)因子；Kbest為對原子n產(chǎn)生作用力的所有原子；hni為原子n與原子i之間的距離。

Kbest、hni的定義為：

式中,N 為原子總數(shù),φ為原子的長度,rni為原子n與原子i之間的歐幾里得距離。

原子質(zhì)量mn的定義如下：

式中：fn（t）為原子n 的適應值,fmax（t）、fmin（t）分別為原子的最大適應值和最小適應值。

每次迭代都會通過更新加速度來對原子的速度和位置進行更新,如下所示：

ASO算法的偽代碼如圖2所示：

圖2 ASO算法偽代碼

3.2 自適應調(diào)節(jié)

ASO 算法在每次迭代中,最優(yōu)原子會施加給相鄰原子一個共價鍵約束力Gn,表示當前原子和最優(yōu)原子之間的位置差值,這種約束增強了最優(yōu)信息在原子之間的傳播,從而提高了算法全局尋優(yōu)的性能。為進一步加強最優(yōu)信息在原子種群中的影響,本文引入歷史原子最優(yōu)信息Pn,來更新原子的加速度,以此提高算法的優(yōu)化性能。

歷史原子最優(yōu)信息Pn和更新之后的原子加速度an,如下所示：

式中,xp為原子n的歷史最優(yōu)距離,μ為歷史系數(shù)因子。

式中的μ和ξ為算法中的超參數(shù),分別代表當前最優(yōu)原子信息和歷史原子最優(yōu)信息的置信度。如果μ值遠大于ξ,則表示算法更相信當前最優(yōu)原子信息,而忽略歷史最優(yōu)原子信息,會降低算法的收斂速度；如果ξ值遠大于μ,則表示算法更相信歷史最優(yōu)原子信息,而忽略當前最優(yōu)原子信息,會降低算法的尋優(yōu)能力。因此,超參數(shù)的選取很大程度上影響著算法的收斂速度和精度,為了避免兩種超參數(shù)在迭代過程中出現(xiàn)上述現(xiàn)象,而導致算法的尋優(yōu)能力下降。本文引入自適應的改進思想,避免兩種超參數(shù)在更新中產(chǎn)生突變的情況,以此平衡歷史最優(yōu)解和當前最優(yōu)解的置信度問題,兩種超參數(shù)的更新公式如下：

式中,t為當前迭代次數(shù)；T為最大迭代次數(shù)。

4 算例結(jié)果分析

4.1 算例參數(shù)

為了驗證本文提出的自適應ASO 算法,在微電網(wǎng)優(yōu)化調(diào)度中的可行性,以圖1所述的微電網(wǎng)作為優(yōu)化算法的研究對象,在時長為24 h 的條件下進行研究。微電網(wǎng)中各分布式發(fā)電單元的參數(shù)如表1 所示,表2 為污染物的治理成本和排放系數(shù),圖3 為某地24 h 的光伏出力、風機出力和電負荷功率曲線。

表1 分布式發(fā)電單元參數(shù)

表2 污染物治理成本和排放系數(shù)

圖3 電負荷、光伏、風機功率曲線

4.2 算例結(jié)果分析

本文設(shè)定原子數(shù)目N=100、深度權(quán)重λ=50,來避免AASO 算法陷入局部最優(yōu)解；系數(shù)因子ξ=0.45、歷史系數(shù)因子μ=0.55,來平衡當前最優(yōu)解和歷史最優(yōu)解的置信度問題,以此提高算法收斂速度和精度,并在最大迭代次數(shù)T=160的條件下進行算例仿真。為了證明AASO 算法在微電網(wǎng)優(yōu)化上的可行性,文中分別在孤島和并網(wǎng)模式下的微電網(wǎng)系統(tǒng)進行算例求解,并與經(jīng)典PSO 算法進行對比分析。

4.2.1 孤島模式下算例結(jié)果分析

孤島運行模式指微電網(wǎng)系統(tǒng)與外部電網(wǎng)系統(tǒng)分離,自主獨立運行的一種模式。圖4為AASO算法優(yōu)化的各發(fā)電單元功率曲線,圖5為兩種算法的運行成本對比。

圖4 孤島模式下各發(fā)電單元功率曲線

圖5 孤島模式下微電網(wǎng)運行成本對比曲線

在孤島運行模式下,微電網(wǎng)系統(tǒng)需要根據(jù)自身的負載需求和可用資源來管理電力的生產(chǎn)、消耗和儲存,以確保能夠滿足其內(nèi)部負載需求并保持穩(wěn)定運行。由圖4 可知孤島模式下的各分布式發(fā)電單元的發(fā)電功率,會根據(jù)用電需求的變化不斷地做出調(diào)整。圖4表明光伏、風機的發(fā)電功率不能達到用電需求時,燃料電池、微型燃氣輪機和蓄電池將進行發(fā)電；當發(fā)電功率滿足用電需求時,多余的電量就會向蓄電池進行充電,以備不時之需。由圖5 可知AASO 算法迭代16 次就找到了最優(yōu)值,而PSO 算法需要迭代38 次,AASO 算法優(yōu)化下的微電網(wǎng)運行成本低于PSO算法10.22%。

4.2.2 并網(wǎng)模式下算例結(jié)果分析

并網(wǎng)模式下的微電網(wǎng)系統(tǒng)可以與外部電網(wǎng)連接,以便在需要時購買或出售額外電力,來保持用電負荷的平衡。圖4、圖5 分別為并網(wǎng)運行模式下AASO算法優(yōu)化的各發(fā)電單元功率曲線,以及兩種算法的運行成本對比曲線。

根據(jù)圖6中的算法優(yōu)化結(jié)果可知,在0～7時段用電負荷較少,此時段電價便宜,主要由外部電網(wǎng)以及燃料電池提供電量,微型燃氣輪機發(fā)電成本較高不進行發(fā)電,多余的電量向蓄電池充電；在7～17時段用電負荷上升,此時段電價較高,微型燃氣輪機、燃料電池和蓄電池開始工作,并把多余的電量出售給外部電網(wǎng)。在17～24時段,此時段微型燃氣輪機發(fā)電成本大于購電成本,主要供電方式為燃料電池和外部電網(wǎng),并根據(jù)用電需求來調(diào)整購電策略。由圖7 可知,AASO 算法和PSO 算法的迭代次數(shù)和運行成本分別為22 次、43 次和1 882.74元、2 103.36元,AASO算法優(yōu)化下的微電網(wǎng)運行成本低于PSO算法10.49%。

圖6 并網(wǎng)模式下各發(fā)電單元功率曲線

圖7 并網(wǎng)模式下微電網(wǎng)運行成本對比曲線

通過兩種模式下的算例仿真結(jié)果可知,本文提出的自適應原子搜索優(yōu)化（AASO）算法的收斂速度和優(yōu)化能力,均強于經(jīng)典PSO算法,可以有效降低微電網(wǎng)的運行成本,提高微電網(wǎng)運行的經(jīng)濟性和環(huán)保性。

5 結(jié)論

本文針對微電網(wǎng)優(yōu)化調(diào)度的問題,提出了一種新穎的自適應原子搜索優(yōu)化算法。首先根據(jù)各分布式電源的工作特性,建立了兼顧微電網(wǎng)運行經(jīng)濟性和環(huán)保性的優(yōu)化調(diào)度模型。其次,以原子搜索優(yōu)化算法為理論基礎(chǔ),引入自適應的思想來加強算法的尋優(yōu)收斂能力。最后,通過在孤島和并網(wǎng)模式下的算例仿真,驗證了本文提出的AASO算法在微電網(wǎng)優(yōu)化調(diào)度上的可行性與優(yōu)越性,有效地提高了微電網(wǎng)運行的經(jīng)濟環(huán)保性。