基于權(quán)重自適應(yīng)調(diào)整的混沌量子粒子群算法的城市電動(dòng)汽車充電站優(yōu)化布局

于擎，李菁華，趙前扶，邢春陽(yáng)

（東北電力大學(xué)電氣工程學(xué)院，吉林吉林132012）

0 引 言

隨著環(huán)境的惡化以及能源的短缺，越來(lái)越多的政府、企業(yè)開始關(guān)注城市電動(dòng)汽車。而建立一個(gè)可靠的和多樣化的充電設(shè)施以滿足不同用戶的需要成為了首要解決的問(wèn)題。文獻(xiàn)［1］中預(yù)計(jì)到2020年，中國(guó)絕大部分的乘用車的能量來(lái)源將多是以電力和清潔能源為主，城市中主要的交通工具將以電動(dòng)汽車為主，電力與交通運(yùn)輸?shù)挠行ЫY(jié)合將成為我國(guó)能源發(fā)展的新方向。電動(dòng)汽車不但有著減少碳排放、噪聲低等優(yōu)點(diǎn)，還可以與可再生能源結(jié)合，協(xié)調(diào)調(diào)度，以此增加可再生能源的利用效率，平抑系統(tǒng)負(fù)荷波動(dòng)，減小負(fù)荷曲線的峰谷差，有利于電力系統(tǒng)的安全、穩(wěn)定運(yùn)行。

隨著電動(dòng)汽車的快速發(fā)展，作為配套設(shè)施的電動(dòng)汽車充電站成為學(xué)者關(guān)注的重點(diǎn)，若無(wú)法解決合理布局問(wèn)題將嚴(yán)重影響電動(dòng)汽車的全面推廣。如何對(duì)城市電動(dòng)汽車充電站進(jìn)行合理布局，使其既能減少消費(fèi)者的費(fèi)用又能降低建設(shè)成本，使得總體的效益最大化，成為國(guó)家、企業(yè)重點(diǎn)關(guān)注的問(wèn)題之一。針對(duì)此問(wèn)題，文獻(xiàn)［2］不僅考慮了電動(dòng)汽車的使用率還考慮了城市建設(shè)等因素，建立多目標(biāo)規(guī)劃模型；文獻(xiàn)［3］提出兩階段優(yōu)化布局算法，即先預(yù)測(cè)用戶的使用需求，再針對(duì)用戶的需求建立相應(yīng)的數(shù)學(xué)模型實(shí)現(xiàn)最優(yōu)化；文獻(xiàn)［4］針對(duì)實(shí)際問(wèn)題可變因素較多的問(wèn)題，建立將一些可變因素考慮在內(nèi)的數(shù)學(xué)模型，并用粒子群算法進(jìn)行求解；文獻(xiàn)［5］根據(jù)排隊(duì)論，建立一個(gè)以建設(shè)成本最小為目標(biāo)的模型，將建站所用的費(fèi)用，充電者充電途中的費(fèi)用、排隊(duì)等待的費(fèi)用考慮在內(nèi)，利用差分進(jìn)化混合算法進(jìn)行研究，最后通過(guò)南方某城市為實(shí)例，采用差分進(jìn)化混合粒子群算法對(duì)該模型求解；文獻(xiàn)［6］分別對(duì)國(guó)內(nèi)外電動(dòng)汽車充電站建設(shè)的情況進(jìn)行了深入的說(shuō)明，并對(duì)國(guó)內(nèi)充電站建設(shè)的不足之處給出建議；文獻(xiàn)［7］提出一種新的配電網(wǎng)輻射狀約束，經(jīng)過(guò)算例驗(yàn)證可有效降低配電網(wǎng)的網(wǎng)損和運(yùn)行成本。

國(guó)內(nèi)外對(duì)充電站優(yōu)化問(wèn)題的研究還沒(méi)有做到十分深入，在建模過(guò)程中，考慮的因素越全面，模型的準(zhǔn)確性越高、越合理，因此，在現(xiàn)有模型基礎(chǔ)上應(yīng)考慮更多的影響因素，充電站的規(guī)劃應(yīng)盡量靠近負(fù)荷中心，減少線路損耗，增加電網(wǎng)經(jīng)濟(jì)性；應(yīng)與城市建設(shè)相結(jié)合，不影響交通；應(yīng)盡量選擇在用戶方便充電的位置，如住宅區(qū)或辦公樓宇附近。

根據(jù)我國(guó)的有關(guān)政策和現(xiàn)狀，經(jīng)濟(jì)性是首要關(guān)注的問(wèn)題，因此，本文記及道路地理信息、交通流量、土地成本以及相鄰充電站間距離等因素，對(duì)整個(gè)區(qū)域進(jìn)行統(tǒng)一規(guī)劃，建立了以經(jīng)濟(jì)性最好為目標(biāo)的數(shù)學(xué)模型，包含建站的經(jīng)濟(jì)性和用戶的經(jīng)濟(jì)性。約束條件記及了充電站容量、充電半徑以及相鄰兩充電站間的距離。針對(duì)該模型，采用權(quán)重自適應(yīng)調(diào)整的混沌量子粒子群優(yōu)化算法來(lái)進(jìn)行求解，該算法在迭代中會(huì)動(dòng)態(tài)改變權(quán)重，以加強(qiáng)粒子的搜索能力。最后將結(jié)果按照相鄰充電站間的距離約束進(jìn)行篩選，得出最終建站位置及容量。經(jīng)算例驗(yàn)證，模型合理、算法的精度高并且迭代快。

1 電動(dòng)汽車充電站建設(shè)的數(shù)學(xué)模型

為避免電動(dòng)汽車出現(xiàn)充電不便的情況，文中建立了兩充電站間的距離約束，用于對(duì)最后結(jié)果的篩選；電動(dòng)汽車建站的位置會(huì)直接影響用戶的行駛距離，因此，模型也考慮了用戶的費(fèi)用。本文以文獻(xiàn)［8］中模型為基礎(chǔ)，實(shí)現(xiàn)對(duì)整個(gè)區(qū)域的規(guī)劃。以規(guī)劃年限內(nèi)充電站年均綜合費(fèi)用最小為目標(biāo)，以充電站的服務(wù)能力，服務(wù)半徑和相鄰兩充電站間的距離為約束條件，建立最優(yōu)化數(shù)學(xué)模型：

其中：式（1）為規(guī)劃年限內(nèi)建站的年均總費(fèi)用；Cc為折算后的年均建設(shè)成本；Cr為經(jīng)過(guò)折算后的年均運(yùn)行、維護(hù)成本；Cd為表示用戶的費(fèi)用；式（2）中，r0表示投資回報(bào)率，本文取0.07；m表示計(jì)劃運(yùn)行年限，本文取20年；AL和CL分別表示征地面積及其單價(jià)，Nt和Ct分別表示配電變壓器個(gè)數(shù)及其單價(jià)，NC和CC分別表示充電機(jī)個(gè)數(shù)及其單價(jià)，CG表示固定成本即其他輔助設(shè)備的投資；式（3）中Ncar為所有進(jìn)行充電的電動(dòng)汽車的輛數(shù)，εp為每輛車的平均充電費(fèi)用，α為充電損耗率，β為充電站運(yùn)行成本的系數(shù)，是將維護(hù)、檢修成本和員工工資折算后得到的；式（4）中，將行駛到最近充電站的距離折算為用戶費(fèi)用，θ為先將直線距離折算為曲線距離再折算為費(fèi)用的系數(shù)，r為道路編號(hào)，j為充電站編號(hào)；式（5）為計(jì)算電動(dòng)汽車數(shù)量的公式，Tr表示道路r上的交通流量，ζ表示電動(dòng)汽車的占有率，Ntr表示每臺(tái)電動(dòng)汽車平均產(chǎn)生的交通流量；式（6）～式（8）為約束條件。式（6）表示充電站的充電能力大于最大充電要求。式（7）表示充電站的服務(wù)半徑要大于充電汽車到充電站的距離。式（8）表示任意的相鄰兩充電站之間的距離不能過(guò)遠(yuǎn)，以保證電動(dòng)汽車能夠及時(shí)充電，否則將影響電池的使用壽命甚至影響電動(dòng)汽車的正常使用。γmax為最大同時(shí)充電率；NC為附近充電機(jī)個(gè)數(shù)；ru為用戶到充電站的距離；rc為充電站服務(wù)半徑；d為兩相鄰充電站間的距離。

2 權(quán)重自適應(yīng)調(diào)整的混沌量子優(yōu)化算法

2.1 量子粒子群優(yōu)化算法

在標(biāo)準(zhǔn)粒子群算法（PSO）的基礎(chǔ)上，Sun等結(jié)合量子力學(xué)提出了量子粒子群優(yōu)化算法［9］（QPSO），該算法在搜索過(guò)程中沒(méi)有固定的路線，能夠搜索整個(gè)可行域。因此，QPSO全局搜索能力比PSO好，但存在著容易早熟收斂的缺陷。

2.2 自適應(yīng)混沌量子粒子群優(yōu)化算法

混沌系統(tǒng)由非線性系統(tǒng)演變而來(lái)，它對(duì)初始條件的細(xì)微變化十分敏感并且具有遍歷性，已經(jīng)成為一種有效的優(yōu)化工具。以常用的Logistic映射為例：

當(dāng)x∈（0，1）并且3.56≤μ≤4時(shí)，系統(tǒng)由非線性變?yōu)榛煦纾疚睦没煦绲谋闅v性，對(duì)可行域內(nèi)的粒子進(jìn)行混沌搜索，大大增加搜索多樣性，避免早熟收斂，從而可以跳出局部最優(yōu)解，尋得全局最優(yōu)解。

在QPSO中，慣性權(quán)重ω影響著粒子的運(yùn)動(dòng)性能，進(jìn)而影響粒子的搜索能力。當(dāng)ω值大時(shí)，粒子的全局搜索能力強(qiáng)，收斂的速度較快，但難以獲得精確的解；當(dāng)ω值較小時(shí)，粒子的局部搜索能力強(qiáng)，容易獲得更精確的解，但收斂的速度慢。在搜索過(guò)程中，合理地調(diào)整ω值，很大程度上能夠提高算法的性能，因此，本文采用權(quán)重自適應(yīng)調(diào)整的混沌量子粒子群優(yōu)化算法［10］（ACQPSO），即在粒子進(jìn)化過(guò)程中動(dòng)態(tài)改變權(quán)值且加入混沌算子，提高搜索精度，獲得精確解。

對(duì)粒子群早熟程度評(píng)價(jià)如下：第i個(gè)粒子的適應(yīng)值為fi；fg為群體的最優(yōu)適應(yīng)值；favg為種群的平均適應(yīng)值；f′avg為所有比f(wàn)avg適應(yīng)值高的粒子的平均適應(yīng)值；規(guī)定 Δ＝｜fg-f′avg｜，來(lái)判斷粒子群的早熟程度，Δ值較小時(shí)，粒子群趨于早熟。ACQPSO的自適應(yīng)調(diào)整策略是：將種群分為三個(gè)子群，根據(jù)種群中個(gè)體適應(yīng)值的不同，采用相應(yīng)的自適應(yīng)調(diào)整方式，慣性權(quán)重較小的粒子適合局部尋優(yōu)；慣性權(quán)重較大的粒子適合在早期用于全局尋優(yōu)，在后期跳出局部最優(yōu)，慣性權(quán)重的調(diào)整方式［10］如下：

（1）當(dāng)fi優(yōu)于f′avg時(shí)：此時(shí)的粒子較好，只需賦給粒子較小權(quán)重即可：

（2）當(dāng)fi優(yōu)于favg次于f′avg時(shí)：此時(shí)的粒子情況一般。

式中iter表示當(dāng)前迭代次數(shù)，Maxgen表示最大迭代次數(shù)。

（3）當(dāng)fi次于favg時(shí)：此時(shí)的粒子較差，需要利用參數(shù)k1和k2進(jìn)行調(diào)整：

3 基于ACQPSO算法的電動(dòng)汽車充電站優(yōu)化布局

綜合以上分析，利用ACQPSO算法在具體的電動(dòng)汽車充電站優(yōu)化問(wèn)題中的步驟如下：

（1）初始化各參數(shù)。根據(jù)需要同時(shí)充電的電動(dòng)汽車數(shù)量和各等級(jí)充電機(jī)個(gè)數(shù)，計(jì)算出充電站個(gè)數(shù)n的變化范圍，設(shè)定n值、最大迭代次數(shù)、種群規(guī)模以及終止條件等，輸入充電站和規(guī)劃區(qū)相關(guān)的數(shù)據(jù)；

（2）判斷是否滿足最大迭代次數(shù)或收斂條件，若滿足，直接輸出建站的位置、等級(jí)以及最小費(fèi)用，算法結(jié)束，否則繼續(xù)下一步；

（3）計(jì)算所有粒子的適應(yīng)度，找出并更新個(gè)體最優(yōu)位置pBest和mBest全局平均最優(yōu)值；

（4）更新粒子的位置和速度；

（5）將mBest代入式（9）產(chǎn)生混沌序列。利用混沌序列找到全局最優(yōu)解；

（6）用步驟（5）中的最優(yōu)解代替粒子群中任一粒子；

（7）根據(jù)公式（10～12）調(diào)整慣性權(quán)重，再到第（3）步進(jìn)行計(jì)算，直到滿足終止條件。

4 算例分析

以吉林市某區(qū)實(shí)際道路情況為例：規(guī)劃區(qū)占地面積 40.8 km2，東西距離 7.5 km，南北距離 8.8 km；主干道共2條，次干道4條，支路10條，路口節(jié)點(diǎn)310個(gè)；根據(jù)道路實(shí)時(shí)監(jiān)測(cè)經(jīng)驗(yàn)數(shù)據(jù)可知，該區(qū)日均車流總量6萬(wàn)輛。該規(guī)劃區(qū)內(nèi)有住宅區(qū)，商業(yè)區(qū)和工業(yè)區(qū)。該算例借鑒北京市2014年6月27日出臺(tái)的標(biāo)準(zhǔn)文件，《北京市電動(dòng)汽車推廣應(yīng)用行動(dòng)計(jì)劃（2014年～2017年）》中規(guī)定的4個(gè)等級(jí)充電站，各等級(jí)充電站的相關(guān)信息如表1所示。

根據(jù)吉林市國(guó)土資源局發(fā)布的《吉林市人民政府關(guān)于調(diào)整更新吉林市城區(qū)基準(zhǔn)地價(jià)等土地價(jià)格的通知》，各土地類型的價(jià)格如表2所示，實(shí)際的道路情況及實(shí)測(cè)出的交通流量如表3所示；現(xiàn)實(shí)情況中，由于道路分布狀況、車主行駛習(xí)慣等不同，折算系數(shù)θ為一變量，為了簡(jiǎn)化問(wèn)題，文中的折算系數(shù)θ算作固定值，同理，雖然路況信息時(shí)刻在變化，在本文中也算作固定值。模型中各參數(shù)取值如表4所示。

根據(jù)表3、表4中的數(shù)據(jù)以及式（6），可計(jì)算出規(guī)劃區(qū)域內(nèi)所有電動(dòng)汽車同時(shí)充電的最大值為120輛，利用該值與表1中的數(shù)據(jù)可以計(jì)算出若全部按等級(jí)1的充電站建設(shè)，至少需要建設(shè)3座充電站；若全部按等級(jí)4的充電站建設(shè)，最多需要建設(shè)15座充電站。因此，充電站的個(gè)數(shù)n分別取3～15中的值，針對(duì)每個(gè)n的取值，計(jì)算相應(yīng)的目標(biāo)函數(shù)值。

表1 充電站的等級(jí)及相關(guān)信息Tab.1 Station levels and corresponding information

表2 各類型土地價(jià)格Tab.2 Costs of different land

表3 道路交通流量Tab.3 Traffic flow data

表4 模型中的參數(shù)取值Tab.4 Parameters values of model

參數(shù)取為：種群規(guī)模取50；最大迭代次數(shù)取150；通過(guò)對(duì)算法的實(shí)現(xiàn)，得出建設(shè)3～15座充電站的年均最小綜合費(fèi)用分別為：301.61、297.45、334.04、348.88、347.68、364.34、401.29、400.08、452.48、451.28、485.81、484.60、519.13萬(wàn)元。由此可知，建設(shè)4座充電站費(fèi)用最小。建站總成本隨充電站個(gè)數(shù)的變化情況如圖1所示。

圖1 總成本隨充電站個(gè)數(shù)的變化Fig.1 Relationship between total cost and charging station number

從計(jì)算結(jié)果可以看出：建設(shè)費(fèi)用整體上隨充電站的增加而逐漸增大，但充電站個(gè)數(shù)為4時(shí)，是圖中折線的最低點(diǎn)，即建站總成本最低。以規(guī)劃區(qū)域的最低點(diǎn)為坐標(biāo)原點(diǎn)，建設(shè)4座充電站的具體情況如表5所示。

表5 最優(yōu)方案Tab.5 Optimal solution

對(duì)ACQPSO算法和QPSO算法的比較如表6所示。為了驗(yàn)證的準(zhǔn)確性，針對(duì)本文模型，尋優(yōu)算法都運(yùn)行50次；取所有解的平均值作為算法研究的指標(biāo)；為了研究算法的可擴(kuò)展性，對(duì)模型的決策變量維數(shù)進(jìn)行設(shè)置，分別為10、20。對(duì)比結(jié)果如表6所示。

表6 算法比較Tab.6 Comparison of algorithms

從表中可以看出：ACQPSO算法的可擴(kuò)展性較強(qiáng)，在維數(shù)高時(shí)，算法的性能仍然較好。在維數(shù)達(dá)到20時(shí)，達(dá)優(yōu)率為100%，即找到全局最優(yōu)值。由于ACQPSO算法中引入了混沌搜索，即ACQPSO算法的單步時(shí)間較QPSO算法的單步時(shí)間稍長(zhǎng)，但由于迭代次數(shù)少，ACQPSO算法在總的搜索時(shí)間上仍是較優(yōu)的。由本例可知，當(dāng)電動(dòng)汽車的持有率逐漸增加以及規(guī)劃區(qū)域逐漸增大時(shí)，即問(wèn)題的維數(shù)增加，該算法仍然適用。

圖2 兩種算法的最佳適應(yīng)度曲線比較Fig.2 Best fitness curves comparison of two algorithms

從圖中可以看出：ACQPSO算法比QPSO算法大約少迭代20次，收斂速度快，并且ACQPSO算法的目標(biāo)值比QPSO算法的目標(biāo)值小。通過(guò)對(duì)本文模型的實(shí)現(xiàn)，將結(jié)果按照相鄰兩充電站間的距離不能大于兩充電站的服務(wù)半徑的原則進(jìn)行篩選，滿足該條件的作為最終結(jié)果，最終建站位置如圖3所示。

圖3 建站位置Fig.3 Sketch map of charging stations

從圖3中，可以清楚的看到在某區(qū)充電站的建站情況，以及可以在同一時(shí)間進(jìn)行充電的電動(dòng)汽車的情況。例如，充電站1的等級(jí)是1，在所有等級(jí)中，其服務(wù)能力最強(qiáng)、服務(wù)半徑最大、充電機(jī)的個(gè)數(shù)最多，所以，正如圖中充電站1的連線最多，它能夠同時(shí)容納較多的電動(dòng)汽車充電，并且可以保證每輛電動(dòng)汽車行駛到充電站的路程最短，任意兩充電站間的距離沒(méi)有過(guò)遠(yuǎn)，符合模型的約束條件。由此可見模型的合理性。

5 結(jié)束語(yǔ)

針對(duì)充電站的規(guī)劃問(wèn)題，采用量子粒子群算法進(jìn)行求解，經(jīng)對(duì)吉林市某區(qū)的算例驗(yàn)證，文中的模型及方法均合理、可行。

文中的數(shù)學(xué)模型是以規(guī)劃充電站的經(jīng)濟(jì)性為目標(biāo)，對(duì)規(guī)劃區(qū)域進(jìn)行統(tǒng)一規(guī)劃。并按照兩充電站間的距離約束來(lái)篩選結(jié)果，從而得出最小費(fèi)用、最優(yōu)建站位置以及建站等級(jí)。

文中采用的ACQPSO算法是利用QPSO的快速收斂性以及混沌算子的遍歷性，搜索過(guò)程中增加了搜索的多樣性，通過(guò)對(duì)權(quán)重系數(shù)的動(dòng)態(tài)調(diào)節(jié)，增強(qiáng)粒子的搜索能力。經(jīng)驗(yàn)證，該算法性能較好。

綜上所述，從建立的模型來(lái)看，考慮較全面，從采用的算法來(lái)看，要優(yōu)于之前采用的算法，從結(jié)果來(lái)看，符合實(shí)際，且相比于之前的結(jié)果更優(yōu)。本文中的模型和算法都可以進(jìn)一步推廣到更大面積，更高維數(shù)的充電站規(guī)劃問(wèn)題，可供有關(guān)部門參考。