陳 浩， 王 健

( 國(guó)網(wǎng)安徽省電力有限公司馬鞍山供電公司，安徽 馬鞍山 243011)

0 引言

近年來(lái)，隨著我國(guó)城市化進(jìn)程的加快，電力設(shè)施建設(shè)規(guī)模不斷擴(kuò)大，迫切需要進(jìn)行電力設(shè)施布局專項(xiàng)規(guī)劃，變電站的選址定容是其中的核心部分，關(guān)系到規(guī)劃區(qū)域內(nèi)未來(lái)電網(wǎng)布局和網(wǎng)架結(jié)構(gòu)。長(zhǎng)期以來(lái)，針對(duì)這一問(wèn)題，國(guó)內(nèi)外科技工作者已有大量研究。文獻(xiàn)[1]提出了布谷鳥(niǎo)算法解決變電站規(guī)劃問(wèn)題，但該方法采用就近原則分配供電負(fù)荷，不能保證變電站負(fù)載率滿足要求；文獻(xiàn)[2]應(yīng)用漁夫捕魚(yú)算法求解變電站選址模型，算法需要已知變電站的規(guī)劃容量，否則變電站的容量可能會(huì)超過(guò)允許上限；文獻(xiàn)[3—4]基于加權(quán)Voronoi圖進(jìn)行站址優(yōu)化，可以確定各變電站的供電范圍和容量，但使用交替定位算法進(jìn)行選址難以得到全局最優(yōu)解；文獻(xiàn)[5—6]引入改進(jìn)粒子群算法進(jìn)行變電站選址定容規(guī)劃，然而迭代參數(shù)的變化只能擴(kuò)大搜索范圍，難以完全解決粒子群算法(particle swarm optimization, PSO)的局部最優(yōu)問(wèn)題。

文中提出基于PSO算法與差分進(jìn)化算法(differential evolution, DE)的混合優(yōu)化算法差分粒子群算法(differential particle swarm optimization, DEPSO)進(jìn)行變電站的選址定容優(yōu)化規(guī)劃。該算法形成2個(gè)子群，以PSO子群為尋優(yōu)主體，通過(guò)在PSO子群和DE子群之間建立信息共享和優(yōu)勝劣汰機(jī)制，提高了種群的多樣性，改善了PSO子群的收斂方向，有效解決了PSO算法的局部最優(yōu)問(wèn)題，且搜索效率較高。同時(shí)運(yùn)用Voronio圖在求解過(guò)程中直接劃分出了每個(gè)變電站的供電范圍，得到變電站規(guī)劃容量，校驗(yàn)變電站實(shí)際負(fù)載率，提高了搜索效率和實(shí)際應(yīng)用能力。算例分析表明DEPSO算法全局搜索能力比PSO算法高，收斂速度比DE算法快，規(guī)劃結(jié)果準(zhǔn)確、可靠。

1 變電站選址的數(shù)學(xué)模型

變電站選址問(wèn)題是一個(gè)多目標(biāo)的非線性約束問(wèn)題[7-10]，在滿足目標(biāo)水平年的負(fù)荷需求前提下，以最小投資和年運(yùn)行費(fèi)用為目標(biāo)函數(shù)，確定新建變電站的位置、數(shù)量、容量，已建變電站的遠(yuǎn)期容量以及所有變電站的供電范圍，具體數(shù)學(xué)模型如下：

(1)

約束條件為：

(2)

2 Voronoi圖的定義與性質(zhì)

Voronoi圖，又稱泰森多邊形或Dirichlet圖，由俄國(guó)數(shù)學(xué)家Voronoi于1908年提出，由連接兩鄰點(diǎn)直線的垂直平分線組成的連續(xù)多邊形組成。N個(gè)在平面上有區(qū)別的點(diǎn)，按照最鄰近原則劃分平面，每個(gè)點(diǎn)與它的最近鄰區(qū)域相關(guān)聯(lián)。

設(shè)有二維歐幾里得平面上離散生長(zhǎng)點(diǎn)的集合為pi與pj間的歐式距離，由V(pi)={pd(p,pi)≤d(p,pj),j≠i,i,j∈In}給出的區(qū)域稱為生長(zhǎng)點(diǎn)pi的Voronoi多邊形，而所有生長(zhǎng)點(diǎn)p1,p2,…pi,pj,…pn的Voronoi多邊形的集合V={V(p1),V(p2),V(p3),V(pi),…,V(pn)}構(gòu)成了P的Voronoi圖，如圖1所示。

圖1 離散生長(zhǎng)點(diǎn)的Voronoi圖Fig.1 Voronoi diagram of discrete growth point

Voronoi圖有許多有趣而驚人的數(shù)學(xué)特性，文中主要應(yīng)用的是勢(shì)力范圍特性[11]：由IR2中一組離散生長(zhǎng)點(diǎn)P={p1,…,pn}(n≥2)所生成的Voronoi多邊形。對(duì)于一個(gè)空間生長(zhǎng)目標(biāo)而言，凡落在其Voronoi多邊形范圍內(nèi)的空間點(diǎn)均距其最近，這一特性能夠約束變電站處于負(fù)荷中心，并得出變電站供電范圍，繼而得到變電站的規(guī)劃容量Si，公式如下：

(3)

式中：e為規(guī)定的變電站標(biāo)準(zhǔn)負(fù)載率。

由于變電站規(guī)劃容量只能選擇有限的離散值，故根據(jù)就近原則確定相應(yīng)變電站規(guī)劃容量后，需重新校驗(yàn)各變電站的實(shí)際負(fù)載率e(Si)，對(duì)超出負(fù)載率允許范圍的粒子，在其適應(yīng)值中加入懲罰因子Pu，文中取較大值10 000。

3 智能算法

3.1 PSO算法

PSO算法[12-15]的基本思想是隨機(jī)初始化一群沒(méi)有體積和質(zhì)量的粒子，所有粒子在可行解空間里運(yùn)動(dòng)，粒子的方向和距離將由一個(gè)速度變量來(lái)決定。在每一次迭代中，粒子將跟隨兩個(gè)極值：一個(gè)是該粒子自身目前為止找到的最優(yōu)解，另一個(gè)是整個(gè)粒子群體目前為止找到的最優(yōu)解。標(biāo)準(zhǔn)PSO算法的公式如下：

(4)

(5)

3.2 DE算法

(6)

(7)

式中：j∈[1,D]；rand(j)為[0,1]內(nèi)均勻分布的隨機(jī)數(shù)；PCR為變異概率，PCR∈[0,1]；rand(i)為隨機(jī)選擇指數(shù)，rand(i)∈[1,2,…,D]。

最后進(jìn)行選擇操作，采用貪婪策略，公式如下：

(8)

式中：φ(x)為適應(yīng)度函數(shù)。

3.3 DEPSO算法

PSO算法參數(shù)設(shè)置簡(jiǎn)單，優(yōu)化前期種群的多樣性較高，搜索能力較強(qiáng)，可以在可行解空間內(nèi)快速收斂，但在搜索后期，隨著群體中的粒子逐漸向種群最優(yōu)粒子靠近，群體的多樣性逐步消失，表現(xiàn)為適應(yīng)度值變化緩慢或停止變化，出現(xiàn)早熟現(xiàn)象。DE算法是一種并行隨機(jī)搜索策略算法[17-19]，其特有的記憶能力可以根據(jù)當(dāng)前的搜索狀態(tài)動(dòng)態(tài)調(diào)整搜索策略，進(jìn)行協(xié)同搜索，但在搜索后期，由于其高隨機(jī)性和種群多樣性的下降，DE算法收斂速度變得緩慢而難以達(dá)到全局最優(yōu)。

針對(duì)PSO算法易陷入局部最優(yōu)而DE算法搜索較慢的缺點(diǎn)，文中將兩種算法結(jié)合，提出了基于雙子群混合優(yōu)化的DEPSO算法。該算法將整個(gè)搜索種群分為PSO和DE 2個(gè)子群，以PSO子群為尋優(yōu)主體，DE子群不斷對(duì)其進(jìn)行調(diào)整和優(yōu)化，其具體步驟如下：

(1) 初始化PSO、DE 2個(gè)子群，設(shè)置參數(shù)；

(2) PSO子群粒子按照PSO搜索機(jī)制進(jìn)行迭代，記錄PSO子群最優(yōu)解Pbest及位置，DE子群粒子按照DE搜索機(jī)制進(jìn)行迭代，記錄DE子群最優(yōu)解Dbest及位置，比較Pbest、Dbest大小，以此更新當(dāng)前全局最優(yōu)解Gbest及其位置Xbest；

(3) 在以Xbest為中心，R為半徑的鄰域內(nèi)隨機(jī)生成N個(gè)粒子X(jué)i代替DE子群的所有粒子進(jìn)行迭代，其公式如下，Yi各維為隨機(jī)分量，Yi∈[-1,1]：

Xi=Xbest+RYii=1,2,3…N

(9)

(4) 計(jì)算DE子群的適應(yīng)值，記錄其全局最優(yōu)解Dbest及位置，將整個(gè)種群所有粒子的適應(yīng)值進(jìn)行排序，用適應(yīng)值較好的DE粒子代替相同數(shù)量適應(yīng)值較差的PSO粒子，比較Pbest、Dbest大小，以此更新Gbest及其位置Xbest；

(5) 如果適應(yīng)值滿足精度要求或者達(dá)到最大迭代次數(shù)則算法終止，否則轉(zhuǎn)至步驟(2)。

DEPSO算法的核心在于用整個(gè)種群的當(dāng)前全局最優(yōu)解Gbest代替PSO子群的最優(yōu)解Pbest進(jìn)行參數(shù)調(diào)整，并利用種群當(dāng)前的全局最優(yōu)位置Xbest對(duì)DE子群進(jìn)行更換繼而更新部分較差的PSO粒子。這樣PSO子群可以參照整個(gè)種群的全局最優(yōu)解改變?cè)械牡较?，并通過(guò)DE子群的不斷替換保持種群的多樣性，最終引導(dǎo)PSO子群逃離局部最優(yōu)點(diǎn)。

根據(jù)上述DEPSO算法，變電站選址定容規(guī)劃的流程如圖2所示。

圖2 基于DEPSO算法的變電站選址定容流程Fig.2 The flow chart of substation site and volume selection based on DEPSO algorithm

4 算例分析

依據(jù)DEPSO算法該方法對(duì)某市配電網(wǎng)進(jìn)行遠(yuǎn)期規(guī)劃，該市遠(yuǎn)期建成區(qū)面積165 km2，人口139萬(wàn)；飽和負(fù)荷2 644.7 MW，負(fù)荷點(diǎn)588個(gè)；220 kV變電站6座，容量4760 MW；110 kV變電站33座，其中需新增26座；110 kV變電站單臺(tái)主變?nèi)萘窟x擇為40 MW或50 MW，主變臺(tái)數(shù)最多為3臺(tái)，標(biāo)準(zhǔn)負(fù)載率0.57，允許范圍0.55～0.59，規(guī)劃區(qū)標(biāo)準(zhǔn)容載比1.8。電力平衡如表1所示。

表1 電力平衡表Tab.1 Power balance sheet MW

設(shè)種群規(guī)模為40，PSO算法慣性權(quán)重初始值為0.9，末值為0.4；DE算法的變異概率設(shè)為0.8，縮放因子最大值為0.2，最小值為0.1，鄰域半徑R取200，最大迭代次數(shù)為500 。圖3所示為該市遠(yuǎn)期配電網(wǎng)變電站的規(guī)劃結(jié)果示意，圖中紅色圓點(diǎn)為已有變電站，藍(lán)色圓點(diǎn)為新建變電站。

圖3 變電站遠(yuǎn)期規(guī)劃示意Fig.3 Substation long term planning diagram

本算法所得變電站負(fù)載率在0.55～0.59范圍內(nèi)，遠(yuǎn)期變電站的容量總和為4810 MW，規(guī)劃區(qū)容載比為1.82，滿足規(guī)劃要求。規(guī)劃結(jié)果如表2所示。

表2 規(guī)劃結(jié)果Tab.2 Planning results

為比較算法性能，分別采用文中所提DEPSO算法、PSO算法和DE算法進(jìn)行變電站選址定容規(guī)劃，其尋優(yōu)過(guò)程如圖4所示。

圖4 幾種算法的尋優(yōu)過(guò)程對(duì)比Fig.4 Comparison of optimization algorithms

從圖中可以看出，DE算法在450次仍在搜索，搜索速度慢； PSO算法雖然收斂較快，但迭代不到150次即停止搜索，陷入局部最優(yōu)；而DESPO算法在迭代100次左右接近全局最優(yōu)解，搜索速度較快，有較好的全局尋優(yōu)能力。原因?yàn)镈E算法搜索過(guò)程隨機(jī)性較大，導(dǎo)致后期收斂較慢且難以找到全局最優(yōu)解；PSO算法搜索導(dǎo)向性較強(qiáng)，表現(xiàn)為向局部最優(yōu)解快速收斂，難以逃離局部最優(yōu)； DEPSO算法有效平衡了搜索中的導(dǎo)向性和隨機(jī)性，引入一種新的信息交流機(jī)制，以隨機(jī)性的概率轉(zhuǎn)換代替序貫的確定性機(jī)理轉(zhuǎn)換，使信息能夠在2個(gè)子群中傳遞，調(diào)節(jié)PSO子群的收斂方向，并通過(guò)淘汰機(jī)制，剔除可能造成局部最優(yōu)的個(gè)體，保持種群多樣性，有利于個(gè)體避免因錯(cuò)誤的信息判斷而陷入局部最優(yōu)。

5 結(jié)語(yǔ)

PSO算法將全局最優(yōu)粒子信息共享給其他粒子，使得子群中其他粒子向其單向流動(dòng)，在求解多目標(biāo)非線性復(fù)雜函數(shù)問(wèn)題時(shí)，容易陷入局部極值點(diǎn)，而單一DE算法搜索速度較慢，難以全局尋優(yōu)。文中提出的DEPSO算法將二者混合優(yōu)化，利用DE算法隨機(jī)性的概率轉(zhuǎn)換機(jī)制改變PSO子群的流動(dòng)方向，保持群體的多樣性，引導(dǎo)PSO子群找到全局最優(yōu)點(diǎn)。同時(shí)將變電站選址模型與Voronoi圖相結(jié)合，使變電站處于負(fù)荷中心，確定變電站規(guī)劃容量，滿足負(fù)載率和規(guī)劃區(qū)容載比要求，減少了計(jì)算量。算例結(jié)果表明，文中提出的變電站選址模型和算法可以有效的解決城市配電網(wǎng)變電站規(guī)劃問(wèn)題，具有較高的運(yùn)算速度和精度。

參考文獻(xiàn)：

[1] 邵 萌，李郁俠，譚雅嵐，等. 基于布谷鳥(niǎo)搜索算法的變電站選址方法研究[J]. 電網(wǎng)與清潔能源，2016，32(11)：51-55．

SHAO Meng，LI Yuxia，TAN Yalan，et al. Study on substation locating method based on cuckoo search algorithm [J].Power System and Clean Energy，2016，32(11)：51-55．

[2] 王澤黎. 基于小生境漁夫捕魚(yú)算法的變電站規(guī)劃[J]. 電力系統(tǒng)保護(hù)與控制，2014，42(16)：84-88．

WANG Zeli. Optimal substation palnning based on niche fisher fishing algorithm [J]. Power System Protection and Control，2014，42(16)：84-88．

[3] 葛少云，李 慧，劉 洪. 基于加權(quán)Voronoi圖的變電站優(yōu)化規(guī)劃[J]. 電力系統(tǒng)自動(dòng)化，2007，31(3)：29-33.

GE Shaoyun，LI Hui，LIU Hong. Substation optimization based on weighted Voronoi diagram[J]. Automation of Electric Power Systems，2007，31(3)：29-33.

[4] 曹 昉，孟琦斌，苗培青，等. 基于改進(jìn)加權(quán)Voronoi圖和遺傳算法的變電站規(guī)劃[J]. 電網(wǎng)技術(shù)，2015，39(2):511-516．

CAO Fang，MENG Qibin，MIAO Peiqing，et al. Substation planning based on improved weighted Voronoi diagram and genetic algorithm[J].Power System Technology，2015，39(2):511-516．

[5] 嚴(yán)司瑋． 基于改進(jìn)粒子群算法的變電站兩階段優(yōu)化選址[D]. 長(zhǎng)沙：湖南大學(xué)，2009：1-12．

YAN Siwei．Two stage optimal location of substation based on improved particle swarm optimization[D]. Changsha：Hunan University，2009：1-12．

[6] 劉自發(fā)，張建華. 基于改進(jìn)多組織粒子群體優(yōu)化算法的配電網(wǎng)變電站選址定容[J]. 中國(guó)電機(jī)工程學(xué)報(bào)，2007，38(19)：235-239．

LIU Zifa，ZHANG Jianhua. Optimal planning of substation locating and sizing based on refined muti-team PSO algorithm[J]. Proceedings of the CSEE. 2007，38(19): 235-239．

[7] 王成山，魏海洋. 變電站選址定容兩階段優(yōu)化規(guī)劃 [J]． 電力系統(tǒng)自動(dòng)化，2005，29(3)：41-46．

WANG Chengshan，WEI Haiyang. Two-stage optimal planning of substation locating and sizing [J]. Automation of Electric Power Systems，2005，29(3)：41-46．

[8] 劉友強(qiáng)，李欣然. 云理論在配電網(wǎng)絡(luò)變電站選址定容中的應(yīng)用[J]. 中國(guó)電機(jī)工程學(xué)報(bào)，2014(4):48-52．

LIU Youqiang，LI Xinran． Application of cloud theory to optimal planning of substation locating and sizing [J]. Proceedings of the CSEE，2014(4) :48-52．

[9] HAFHIFAM M R，SHAOHABI M． Optimal location and sizing of HV/MV substations in uncertainty load environment using genetic algorithm[J]． Electric Power Systems Research，2012，63(11)：37-50．

[10] 李鑫濱，朱慶軍. 變電站選址定容新模型及其遺傳算法優(yōu)化[J].電力系統(tǒng)及其自動(dòng)化學(xué)報(bào)，2009，21(3)：32-35．

LI Xinbin，ZHU Qingjun．New model optimized by genetic algorithm for distribution substation locating and sizing[J]. Proceedings of the CSU-EPSA，2009，21(3) ：32-35．

[11] 陳 軍. Voronoi動(dòng)態(tài)空間數(shù)據(jù)模型 [M]. 北京：測(cè)繪出版社，2002：101-120．

CHEN Jun. Voronoi dynamic spatial data model [M]. Beijing: Mapping Press，2002：101-120．

[12] 周建華，孫 蓉，陳久林，等. 基于改進(jìn)PSO算法的含風(fēng)電場(chǎng)電力系統(tǒng)網(wǎng)架擴(kuò)展規(guī)劃研究[J]. 電力工程技術(shù)，2014，33(5)： 28-31．

ZHOU Jianhua，SUN Rong，CHEN Jiulin，et al． Study on expansion planning of power system network containing wind farm based on improved PSO algorithm [J]. Electric Power Engineering Technology，2014，33(5)： 28-31．

[13] 李 麗，牛 奔. 粒子群優(yōu)化算法 [M]. 北京：冶金工業(yè)出版社，2009：70-99．

LI Li，LIU Ben. Particle swarm optimization [M]. Beijing: Metallurgical Industry Press, 2009: 70-99．

[14] 劉文霞，張林峰，張建華． 基于GIS和PSO的智能變電站選址定容[C]∥中國(guó)國(guó)際供電會(huì)，2006．

LIU Wenxia，ZHANG linfeng. Intelligent substation locating and sizing based on GIS and PSO [C]∥CICED，2006．

[15] L Z S，HOU Z R. Particle swarm optimization with adaptive mutation[J]. ActaElectronica Sinica，2009，32(3):416-420.

[16] BECERRA R L，COELLO C A. Culturizing differential evolution for constrained optimization[J] .IEEE Service Center，Mexico，2004，70(2):300-312.

[17] 黃 映. 基于改進(jìn)混合差分進(jìn)化算法的輸電網(wǎng)規(guī)劃研究[J]. 電力工程技術(shù)，2015，34(4)： 36-43．

HUANG Ying. Substation optimization based on weighted Voronoi diagram[J]. Electric Power Engineering Technology, 2015，34(4)： 36-43．

[18] 賴文海，陳賢陽(yáng)，明國(guó)鋒，等. 基于雙種群小生境差分進(jìn)化算法的動(dòng)態(tài)經(jīng)濟(jì)調(diào)度[J]. 廣東電力，2016，29(7)：83-87,92.

LAI Wenhai, CHEN Xianyang, MING Guofeng, et al. Dynamic economic dispatching based on dual population Niche differential evolution algorithm[J]. Guangdong Electric Power，2016，29(7)：83-87,92.

[19] 林藝城, 李錦焙, 孟安波，等. 基于改進(jìn)差分縱橫交叉算法的多區(qū)域動(dòng)態(tài)經(jīng)濟(jì)調(diào)度 [J]. 廣東電力，2017，30(5)：69-76.

LIN Yicheng, LI Jinpei, MENG Anbo, et al. Multi-area dynamic economic dispatching based on improved differential evolution with crisscrossed optimization algorithm[J]. Guangdong Electric Power，2017，30(5)：69-76.