基于改進(jìn)粒子群算法的牽引供電系統(tǒng)多目標(biāo)優(yōu)化設(shè)計(jì)

陳民武，許臣友，黃文勛，邱　璐，孫小凱

(1.西南交通大學(xué) 電氣工程學(xué)院，四川 成都　610031；2.中鐵第一勘察設(shè)計(jì)院集團(tuán)有限公司，陜西 西安　710043)

傳統(tǒng)的牽引供電系統(tǒng)設(shè)計(jì)過(guò)于依賴設(shè)計(jì)人員的專業(yè)經(jīng)驗(yàn)，主要供電設(shè)施數(shù)量和位置的確定易受主觀因素影響，并且由于有限方案的比選也缺乏全局尋優(yōu)能力，難以保證設(shè)計(jì)容量的最優(yōu)[1]。近年來(lái)，國(guó)內(nèi)多個(gè)鐵路設(shè)計(jì)部門先后引進(jìn)了國(guó)際公認(rèn)的牽引供電系統(tǒng)設(shè)計(jì)軟件，例如德國(guó)SIGNON公司的ELBAS/WEBANET和IFB公司的OpenPowerNet軟件，利用牽引負(fù)荷過(guò)程仿真手段，動(dòng)態(tài)計(jì)算系統(tǒng)內(nèi)潮流分布，計(jì)算精度和效率均優(yōu)于之前廣泛采用的平均運(yùn)量法和概率統(tǒng)計(jì)法。但上述商業(yè)軟件的核心算法不公開(kāi)，且僅能對(duì)給定的設(shè)計(jì)方案進(jìn)行計(jì)算和分析，同樣不具備方案自動(dòng)尋優(yōu)能力。

隨著國(guó)內(nèi)外高速和重載鐵路的快速發(fā)展，針對(duì)牽引供電系統(tǒng)建模和仿真的研究也較多[2-5]。但這些研究均未涉及實(shí)際供電系統(tǒng)的優(yōu)化問(wèn)題。文獻(xiàn)[6]在已知全線牽引能耗分布的前提下提出了供電系統(tǒng)優(yōu)化設(shè)計(jì)的基本思路，針對(duì)給定優(yōu)化目標(biāo)推導(dǎo)了嚴(yán)格的解析表達(dá)式；但實(shí)際線路的供電節(jié)點(diǎn)多，加之模型的復(fù)雜性給求最優(yōu)解帶來(lái)了一定困難。群智能算法的研究為解決復(fù)雜尋優(yōu)問(wèn)題提供了有效途徑，并已經(jīng)在電力系統(tǒng)的無(wú)功優(yōu)化、負(fù)荷預(yù)測(cè)和經(jīng)濟(jì)運(yùn)行等領(lǐng)域取得了成功應(yīng)用[7-8]。因此，解決實(shí)際工程問(wèn)題的尋優(yōu)方法可以建立在相關(guān)數(shù)學(xué)模型和群智能算法的基礎(chǔ)上，充分利用計(jì)算機(jī)強(qiáng)大的計(jì)算能力，使得結(jié)果不斷地沿著改善方案自動(dòng)進(jìn)行，最終收斂于最優(yōu)方案，其分析比較的范圍也是人工選優(yōu)不可比擬的。

本文將群智能算法和車—網(wǎng)耦合系統(tǒng)交互仿真相結(jié)合，開(kāi)展?fàn)恳╇娤到y(tǒng)多目標(biāo)優(yōu)化設(shè)計(jì)研究，并結(jié)合實(shí)際工程算例驗(yàn)證優(yōu)化模型和算法的正確性及有效性。

1　基于車—網(wǎng)耦合系統(tǒng)交互仿真的牽引供電系統(tǒng)多目標(biāo)優(yōu)化模型

1.1　車—網(wǎng)耦合系統(tǒng)模型以及潮流分布計(jì)算算法

基于多導(dǎo)體傳輸線理論構(gòu)建適用于各種牽引供電方式的牽引網(wǎng)分布參數(shù)等效模型，將在線列車、牽引變電所(AT供電方式下含AT所和分區(qū)所)以及其他與牽引網(wǎng)橫向連接的電氣設(shè)備等效為分割切面模型。利用牽引變電所三相—兩相電氣量通用變換關(guān)系，構(gòu)建牽引變電所數(shù)學(xué)模型[1]。當(dāng)前高速和重載鐵路推廣采用自耦變壓器(AT)供電方式，該供電方式下供電臂鏈?zhǔn)杰嚒W(wǎng)耦合系統(tǒng)模型如圖1所示。

圖1　AT供電方式下車—網(wǎng)耦合系統(tǒng)供電臂鏈?zhǔn)骄W(wǎng)絡(luò)模型

結(jié)合AT供電方式下車—網(wǎng)耦合系統(tǒng)供電臂鏈?zhǔn)骄W(wǎng)絡(luò)模型，可以得到統(tǒng)一的節(jié)點(diǎn)導(dǎo)納矩陣Y[2-3]。與電力系統(tǒng)其他電氣設(shè)備的負(fù)荷特性不同，電力機(jī)車(含動(dòng)車組)是時(shí)間—空間移動(dòng)負(fù)荷，隨著它在不同時(shí)刻的分布位置不同，牽引供電系統(tǒng)的拓?fù)浣Y(jié)構(gòu)也會(huì)發(fā)生變化，上述系統(tǒng)的鏈?zhǔn)骄W(wǎng)絡(luò)結(jié)構(gòu)和導(dǎo)納矩陣也是時(shí)變的。在設(shè)計(jì)階段，根據(jù)運(yùn)行組織計(jì)劃安排，可以得到任意仿真時(shí)間步長(zhǎng)T對(duì)應(yīng)全線列車位置。結(jié)合線路條件和機(jī)車特性，通過(guò)牽引計(jì)算得到電力機(jī)車在當(dāng)前位置的取用功率，假設(shè)有功功率為P，無(wú)功功率為Q，則機(jī)車初始電流I(0)為

(1)

式中：U(0)為列車對(duì)應(yīng)切面上接觸網(wǎng)與鋼軌之間的初始電壓；j為虛數(shù)單位；*表示求復(fù)數(shù)的共軛。

在建立車—網(wǎng)耦合系統(tǒng)仿真模型的基礎(chǔ)上，對(duì)應(yīng)已知牽引供電系統(tǒng)設(shè)施數(shù)目和位置，通過(guò)交互迭代計(jì)算，可以求解牽引負(fù)荷過(guò)程的潮流分布，從而得到第i個(gè)牽引變電所的計(jì)算容量Si和全線牽引變電所計(jì)算容量之和，即全線總?cè)萘縎sub。負(fù)荷潮流分布計(jì)算時(shí)的迭代算法流程如圖2所示。

按照?qǐng)D2所示的迭代算法流程，前期研究已經(jīng)完成了上述牽引計(jì)算和供電計(jì)算仿真平臺(tái)，并在國(guó)內(nèi)多條高速鐵路和重載鐵路牽引供電系統(tǒng)設(shè)計(jì)中加以應(yīng)用，與ELBAS/WEBANET軟件對(duì)比，計(jì)算結(jié)果基本一致，達(dá)到了工程設(shè)計(jì)要求[3,9]。

1.2　牽引供電系統(tǒng)多目標(biāo)優(yōu)化模型

AT供電方式下?tīng)恳冸娝毡椴捎?臺(tái)中抽式單相變壓器構(gòu)造Vx接線[10]。假設(shè)某電氣化復(fù)線區(qū)段牽引供電系統(tǒng)初始設(shè)計(jì)方案中包含NTS個(gè)牽引變電所，以全線總?cè)萘縎sub和平均有功功率損耗Ploss均最小為優(yōu)化目標(biāo)，以牽引變電所、AT所和分區(qū)所的數(shù)量NTS,NAT,NSP及其位置LTS,LAT,LSP為優(yōu)化變量(其中LTS,LAT,LSP分別為長(zhǎng)度是NTS,NAT,NSP的一維數(shù)組)，通過(guò)調(diào)整上述變量的取值，反復(fù)進(jìn)行負(fù)荷潮流分布計(jì)算，進(jìn)而得到各個(gè)牽引變電所饋線電流的大小，最終決定各個(gè)牽引變電所輸出的功率以及整個(gè)系統(tǒng)的功率損失。

圖2　迭代算法流程圖

以牽引供電系統(tǒng)主要設(shè)計(jì)原則為約束條件，以Ssub和Ploss最小為優(yōu)化目標(biāo)[11]，則牽引供電系統(tǒng)的多目標(biāo)優(yōu)化模型可以表示為

(2)

f2=minPloss

(3)

s.t.

ULmin≤ULk≤ULmax

NTS,NAT,NSP≠0

LTS∈A

LAT∈B

LSP∈C

式中：NS為采樣點(diǎn)數(shù)；Pih為仿真時(shí)刻h下?tīng)恳冸娝鵬輸出的有功功率；K為在線列車數(shù)目；PLkh為仿真時(shí)刻h下電力機(jī)車k的有功功率；Tsum為仿真時(shí)間；ULk為電力機(jī)車k端的接觸網(wǎng)對(duì)地電壓；A，B和C分別為牽引變電所、AT所和分區(qū)所所址位置的可選值域，且根據(jù)實(shí)際設(shè)計(jì)需要，還可以補(bǔ)充其他約束條件；ULmin和ULmax分別為機(jī)車牽引特性對(duì)ULk的限值要求。

在車—網(wǎng)耦合系統(tǒng)模型中，調(diào)整優(yōu)化變量NTS,NAT,NSP和LTS,LAT,LSP的取值，在改變牽引供電設(shè)施布置和供電臂長(zhǎng)度的同時(shí)，也改變了牽引供電系統(tǒng)負(fù)荷潮流分布和機(jī)車位置的接觸網(wǎng)對(duì)地電壓ULk。在優(yōu)化過(guò)程中如果ULk超出了上述限值要求，則將對(duì)應(yīng)的牽引供電系統(tǒng)設(shè)計(jì)方案排除。

2　基于Pareto熵的混沌多目標(biāo)粒子群優(yōu)化算法

2.1　多目標(biāo)粒子群優(yōu)化算法

在1.2節(jié)優(yōu)化模型中，D維優(yōu)化變量x=(NTS,NAT,NSP,LTS,LAT,LSP)T，其中維數(shù)D=NTS+NAT+NSP+3，目標(biāo)向量F=(f1,f2)T=(minSsub, minPloss)T。將上述最小多目標(biāo)優(yōu)化問(wèn)題描述為

F(x)=(f1(x),f2(x))T

(4)

多目標(biāo)優(yōu)化問(wèn)題中的各個(gè)目標(biāo)往往是相互對(duì)抗的，大多數(shù)情況下各個(gè)目標(biāo)同時(shí)達(dá)到最優(yōu)值的可能性比較小，只能是這些目標(biāo)函數(shù)之間的折中解，這就需要引入Pareto占優(yōu)的概念。設(shè)優(yōu)化變量集合中任意2個(gè)D維優(yōu)化變量x1和x2，當(dāng)滿足l=1,2，fl(x1)

(5)

式中：t為迭代次數(shù)；w為慣性權(quán)重系數(shù)，w≥0；b1和b2為在[0,1]上均勻分布的隨機(jī)數(shù)；pBestz表示第z個(gè)粒子在第t次迭代時(shí)的個(gè)體最優(yōu)位置；c1,c2為學(xué)習(xí)因子(或稱加速因子)。

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

2.2.1Pareto熵及進(jìn)化狀態(tài)檢測(cè)

熵是微觀狀態(tài)均勻性和多樣性的一種度量，表示系統(tǒng)的無(wú)序程度。儲(chǔ)存在外部檔案中的Pareto解的分布均勻性可以由Pareto熵來(lái)反映，Pareto熵間接地表征了Pareto前沿解的多樣性狀態(tài)，前后2次迭代熵的變化(差熵)反映了Pareto前沿重新分布的情況，從而通過(guò)差熵可以推斷種群的進(jìn)化狀態(tài)，如收斂狀態(tài)、多樣化狀態(tài)和停滯狀態(tài)，通過(guò)狀態(tài)反饋信息動(dòng)態(tài)的調(diào)整進(jìn)化策略，平衡種群的開(kāi)采和開(kāi)發(fā)能力。Pareto熵的求解見(jiàn)文獻(xiàn)[14]。

定義相鄰2次迭代的差熵為

ΔE(t)=E(t)-E(t-1)

(6)

式中：E為Pareto解的熵；ΔE為差熵。

2.2.2引入個(gè)性化混沌因子的自適應(yīng)進(jìn)化策略

粒子群優(yōu)化算法可以通過(guò)進(jìn)化狀態(tài)反饋信息，動(dòng)態(tài)調(diào)整慣性權(quán)重系數(shù)w和學(xué)習(xí)因子c1，c2以平衡算法的全局和局部搜索能力。利用混沌序列的遍歷性提出自適應(yīng)調(diào)整進(jìn)化策略。

首先，對(duì)每個(gè)粒子引入個(gè)性化混沌因子γ(t)=μ(z)γ(t-1)(1-γ(t-1))，γ(0)∈rand(0,1)，μ(z)∈rand(3.75, 4]，通過(guò)混沌序列可以提高算法的個(gè)性化遍歷能力[15]。利用種群收斂狀態(tài)和差熵對(duì)其所有運(yùn)動(dòng)參數(shù)進(jìn)行混沌擾動(dòng)，建立如下調(diào)整關(guān)系。

1)個(gè)性化混沌局部學(xué)習(xí)因子

c1(t)=

(7)

2)個(gè)性化混沌全局學(xué)習(xí)因子

c2(t)=

(8)

3)個(gè)性化混沌慣性系數(shù)

w(t)=

(9)

式中：Stepw，Stepc1和Stepc2分別為w，c1和c2的調(diào)節(jié)步長(zhǎng)，分別等于對(duì)應(yīng)參數(shù)區(qū)間長(zhǎng)度除以最大迭代次數(shù)。所有參數(shù)限制在預(yù)先確定的區(qū)間內(nèi)，若超出邊界，則取邊界值。

2.2.3自適應(yīng)變異操作

(10)

變異概率計(jì)算公式為

(11)

式中：Tmax為最大迭代次數(shù)。

可見(jiàn)，當(dāng)Pareto熵E(t)減小(即群多樣性下降)時(shí)變異概率ρm增大，需要擴(kuò)大搜索空間，以避免陷入局部極值。

3　基于改進(jìn)粒子群算法的牽引供電系統(tǒng)多目標(biāo)優(yōu)化過(guò)程

將基于Pareto熵的混沌多目標(biāo)粒子群優(yōu)化算法與牽引供電系統(tǒng)優(yōu)化模型相結(jié)合，確定優(yōu)化變量為全線牽引變電所、分區(qū)所和AT所的數(shù)目和位置，進(jìn)而由潮流計(jì)算得到對(duì)應(yīng)的供電容量需求。首先，在優(yōu)化變量可行域范圍內(nèi)初始化，生成1個(gè)粒子群，粒子初始速度是在速度范圍內(nèi)的1組隨機(jī)數(shù)。然后，代入車—網(wǎng)耦合系統(tǒng)交互仿真平臺(tái)，計(jì)算潮流分布和優(yōu)化目標(biāo)函數(shù)值。在迭代計(jì)算過(guò)程中，采用近似Pareto分布熵及差熵評(píng)估種群的進(jìn)化狀態(tài)，并以此作為反饋信息，動(dòng)態(tài)跟蹤和調(diào)節(jié)進(jìn)化策略和變異算子，利用混沌擾動(dòng)對(duì)變量進(jìn)行調(diào)整。通過(guò)協(xié)調(diào)多個(gè)目標(biāo)函數(shù)之間的關(guān)系，計(jì)算滿足約束條件的Pareto解集。具體優(yōu)化過(guò)程如圖3所示。

圖3　牽引供電系統(tǒng)優(yōu)化流程圖

為了進(jìn)一步確定最優(yōu)方案，應(yīng)用模糊隸屬度函數(shù)評(píng)價(jià)每個(gè)Pareto解中各個(gè)目標(biāo)函數(shù)對(duì)應(yīng)的隸屬度，定義模糊隸屬度函數(shù)如式(12)所示。

(12)

對(duì)于Pareto解集中每個(gè)解，應(yīng)用式(13)求其對(duì)應(yīng)的滿意度值，將滿意度值最大的解所對(duì)應(yīng)的設(shè)計(jì)方案記為最優(yōu)方案。

(13)

式中：Ψd為第d個(gè)Pareto解的滿意度值；G為外部檔案最大容量。

4　算例分析

以某高速鐵路牽引供電系統(tǒng)的優(yōu)化設(shè)計(jì)為例，主要設(shè)計(jì)原則和技術(shù)條件如下：①外部電源電壓等級(jí)220 kV；②復(fù)線區(qū)段，公里標(biāo)DK 331.2—DK460.6，全長(zhǎng)129.4 km；③供電臂采用全并聯(lián)供電方式設(shè)計(jì)，在供電臂末端分區(qū)所、供電臂中部AT所實(shí)現(xiàn)并聯(lián)；④牽引變壓器采用Vx接線方式；⑤重聯(lián)機(jī)車類型為CRH3型動(dòng)車組；⑥牽引網(wǎng)持續(xù)最高電壓為27.5 kV、短時(shí)最高電壓為29 kV、最低電壓正常供電時(shí)不低于20 kV；⑦運(yùn)行組織近期按照4 min追蹤運(yùn)行。根據(jù)線路條件和機(jī)車運(yùn)行特性，得到單車牽引負(fù)荷曲線計(jì)算結(jié)果如圖4所示。

圖4　單車牽引負(fù)荷曲線

此外，以牽引變電所、AT所和分區(qū)所的位置作為多目標(biāo)粒子群優(yōu)化算法中的優(yōu)化變量，全線牽引變電所、分區(qū)所和AT所的選址受多種實(shí)際因素影響。為了避免設(shè)計(jì)方案無(wú)法實(shí)施，比如受線路平縱面、外部電源方案或者設(shè)備運(yùn)輸及生活便利等條件限制，某些位置不能或者不宜作為所址等，可以通過(guò)修改優(yōu)化模型中所址值域A，B和C，在全局尋優(yōu)過(guò)程中自動(dòng)避開(kāi)。

結(jié)合本文算例，基于Pareto熵的混沌多目標(biāo)粒子群優(yōu)化算法主要參數(shù)取值如下：最大迭代次數(shù)Tmax=50，種群規(guī)模N=20，外部檔案最大容量G=25，學(xué)習(xí)因子c1max=c2max=2.5，c1min=c2min=0.5，學(xué)習(xí)因子c1和c2的初值均取1.5，wmax=0.9，wmin=0.4，慣性權(quán)重系數(shù)w初值取0.9。在迭代過(guò)程中，基于Pareto熵的混沌自適應(yīng)參數(shù)調(diào)整策略如圖5所示。由圖5可見(jiàn)，在迭代初期取相對(duì)較大的w，c1和相對(duì)較小的c2，使得粒子較均勻地分布到搜索空間中去，以增加粒子群的多樣性和較強(qiáng)的開(kāi)發(fā)能力，避免早熟；隨著迭代次數(shù)的增加w，c1逐漸減小，c2逐漸增大，粒子群處于收斂狀態(tài)，不斷地跳出局部極值向全局極值靠攏；迭代后期的w，c1和c2趨于平穩(wěn)，迭代處于多樣化狀態(tài)，種群具有較大開(kāi)采能力，局部搜索能力增強(qiáng)。

圖5　自適應(yīng)參數(shù)調(diào)整策略曲線

最終得到Pareto解集在目標(biāo)空間中的分布情況，如圖6所示。

圖6　Pareto前沿

由圖6可見(jiàn)，Pareto前沿具有相當(dāng)直觀的均勻性，體現(xiàn)了算法既兼顧收斂性又不失多樣性的特征，說(shuō)明了基于Pareto熵的混沌自適應(yīng)參數(shù)調(diào)整策略是有效的。通過(guò)模糊隸屬度函數(shù)評(píng)價(jià)每個(gè)Pareto解，求出每個(gè)解的滿意度值，選出滿意度值最大的解作為最優(yōu)方案，見(jiàn)表1。表中牽引變電所、分區(qū)所和AT所的最佳位置以線路公里標(biāo)表示。

表1　滿意優(yōu)化方案中供電設(shè)施布局

為了驗(yàn)證上述優(yōu)化方案的正確性和有效性，由牽引供電專業(yè)設(shè)計(jì)人員采用ELBAS/WEBANET軟件(人工方案1)和傳統(tǒng)的供電計(jì)算方法(人工方案2)，分別獨(dú)立給出牽引供電設(shè)計(jì)方案，3種方法的計(jì)算結(jié)果對(duì)比情況見(jiàn)表2。由于牽引變電所均采用中抽式單相變壓器構(gòu)成Vx接線方式，因此每個(gè)供電臂計(jì)算容量對(duì)應(yīng)單相變壓器原邊繞組容量。表中Smax代表負(fù)荷過(guò)程中最大容量需求，Srms代表負(fù)荷過(guò)程對(duì)應(yīng)的均方根值容量。由表2可見(jiàn)，利用多目標(biāo)尋優(yōu)模型，調(diào)整牽引變電所、分區(qū)所和AT所的布局，使得各供電臂的牽引負(fù)荷均衡性更好。以牽引變電所s1的供電區(qū)間為例，方案差異較大的原因分析如下。

(1)該區(qū)間位于山區(qū)長(zhǎng)大坡道范圍內(nèi)，機(jī)車牽引特性決定了其持續(xù)滿功率運(yùn)行，如圖4中牽引負(fù)荷曲線所示。為了保證接觸網(wǎng)對(duì)地電壓的約束條件，需要牽引供電系統(tǒng)提供足夠的供電能力。因此，該牽引變電所內(nèi)單相變壓器的計(jì)算容量顯著偏大。

(2)該區(qū)段橋隧比例較大，限制了各類型供電設(shè)施的設(shè)所條件。設(shè)計(jì)人員僅能依賴專業(yè)經(jīng)驗(yàn)和主觀判斷，手動(dòng)調(diào)整供電設(shè)施布置方案，缺乏可行域內(nèi)全局搜索能力。優(yōu)化算法則以目標(biāo)函數(shù)滿意度最大為收斂條件，通過(guò)不斷調(diào)整優(yōu)化變量的取值，沿著改善供電方案目標(biāo)函數(shù)的趨勢(shì)自動(dòng)迭代。

不同設(shè)計(jì)方案在主要目標(biāo)函數(shù)上也存在一定差異，對(duì)比結(jié)果見(jiàn)表3。由表3可見(jiàn)，滿意優(yōu)化方案中全線牽引變電所總計(jì)算容量最小，基于仿真計(jì)算軟件得到的結(jié)果居中，傳統(tǒng)手工供電計(jì)算方法得到的結(jié)果偏大。此外，滿意優(yōu)化方案中牽引供電系統(tǒng)的平均有功損耗也較小。造成上述差異的原因分析如下。

(1)既有牽引供電系統(tǒng)仿真計(jì)算軟件僅能開(kāi)展?fàn)恳?fù)荷過(guò)程模擬，不具備設(shè)計(jì)尋優(yōu)能力。本文通過(guò)將供電仿真與智能優(yōu)化相結(jié)合，充分利用計(jì)算機(jī)強(qiáng)大的計(jì)算和高效的搜索能力，實(shí)現(xiàn)了多目標(biāo)函數(shù)下的滿意優(yōu)化設(shè)計(jì)。

(2)傳統(tǒng)供電計(jì)算方法一般取區(qū)間平均帶電電流描述單車牽引負(fù)荷特性，在給定設(shè)所數(shù)目及所址的前提下計(jì)算供電容量，通過(guò)對(duì)有限供電方案的反復(fù)校核和比選，最終確定牽引供電方案。雖然簡(jiǎn)化了計(jì)算過(guò)程，但降低了計(jì)算精度，從而在一定程度上放大了供電容量需求和平均功耗。

表2　多種設(shè)計(jì)方案各變電所供電容量對(duì)比

表3　多種設(shè)計(jì)方案目標(biāo)函數(shù)值對(duì)比

通過(guò)上述工程算例的分析，驗(yàn)證了本文提出的牽引供電系統(tǒng)多目標(biāo)優(yōu)化方法的正確性和有效性, 增強(qiáng)了最終方案的技術(shù)經(jīng)濟(jì)性。

5　結(jié)　語(yǔ)

以牽引供電系統(tǒng)主要設(shè)計(jì)原則作為約束條件，以全線總?cè)萘孔钚『推骄泄β蕮p耗最小為目標(biāo)函數(shù)，建立了牽引供電系統(tǒng)多目標(biāo)優(yōu)化模型。提出了基于Pareto熵的混沌多目標(biāo)粒子群優(yōu)化算法，結(jié)合模糊隸屬度函數(shù)計(jì)算，使得目標(biāo)函數(shù)盡可能達(dá)到最優(yōu)，得到滿意的優(yōu)化設(shè)計(jì)方案。結(jié)合實(shí)際工程算例，對(duì)比既有設(shè)計(jì)方法，驗(yàn)證了上述模型和算法的正確性及有效性。

實(shí)際工程設(shè)計(jì)中除了技術(shù)和經(jīng)濟(jì)因素外，還需考慮社會(huì)因素和環(huán)境因素。在本文研究的基礎(chǔ)上可以調(diào)整優(yōu)化變量的可行域及約束條件，進(jìn)一步完善多目標(biāo)優(yōu)化數(shù)學(xué)模型，提高尋優(yōu)算法的效率，更好地滿足電氣化鐵路牽引供電系統(tǒng)精細(xì)化設(shè)計(jì)的要求。

[1]陳民武，李群湛，智慧，等.牽引供電系統(tǒng)設(shè)計(jì)方案的綜合評(píng)判[J].高電壓技術(shù)，2010，36(2)： 525-530.

(CHEN Minwu，LI Qunzhan，ZHI Hui，et al.Comprehensive Evaluation for Design Scheme of Traction Power Supply System[J].High Voltage Engineering，2010，36(2)：525-530. in Chinese)

[2]BATTISTELLI L，PAGANO M，PROTO D.2×25-kV 50Hz High-Speed Traction Power System：Short-Circuit Modeling[J].IEEE Transactions on Power Delivery，2011，26(3)：1459-1466.

[3]吳命利.牽引供電系統(tǒng)電氣參數(shù)與數(shù)學(xué)模型研究[D].北京：北京交通大學(xué)，2006.

[4]陳宏偉，耿光超，江全元.電氣化鐵路牽引供電系統(tǒng)車網(wǎng)耦合的潮流計(jì)算方法[J].電力系統(tǒng)自動(dòng)化，2012，36(3)：76-80，110.

(CHEN Hongwei, GENG Guangchao, JIANG Quanyuan. Power Flow Algorithm for Traction Power Supply System of Electric Railway Based on Locomotive and Network Coupling[J]. Automation of Electric Power Systems, 2012, 36(3): 76-80，110.in Chinese)

[5]胡海濤，何正友，王江峰，等.基于車網(wǎng)耦合的高速鐵路牽引網(wǎng)潮流計(jì)算[J].中國(guó)電機(jī)工程學(xué)報(bào)，2012， 32(19)：101-108.

(HU Haitao，HE Zhengyou，WANG Jiangfeng，et al.Power Flow Calculation of High-Speed Railway Traction Network Based on Train-Network Coupling Systems[J].Proceedings of the CSEE，2012，32 (19): 101-108. in Chinese)

[6]宮衍圣.牽引供電系統(tǒng)優(yōu)化設(shè)計(jì)技術(shù)的應(yīng)用研究[D].成都：西南交通大學(xué)，2002.

[7]湯慶峰，劉念，張建華，等.基于EMD-KELM-EKF與參數(shù)優(yōu)選的用戶側(cè)微電網(wǎng)短期負(fù)荷預(yù)測(cè)方法[J].電網(wǎng)技術(shù)，2014，38(10)：2691-2699.

(TANG Qingfeng，LIU Nian，ZHANG Jianhua，et al.A Short-Term Load Forecasting Method for Micro-Grid Based on EMD-KELM-EKF and Parameter Optimization[J].Power System Technology，2014，38(10): 2691-2699. in Chinese)

[8]沈玉明，胡博，謝開(kāi)貴，等.計(jì)及儲(chǔ)能壽命損耗的孤立微電網(wǎng)最優(yōu)經(jīng)濟(jì)運(yùn)行[J].電網(wǎng)技術(shù)，2014，38(9)：2371-2378.

(SHEN Yuming，HU Bo，XIE Kaigui，et al.Optimal Economic Operation of Isolated Microgrid Considering Battery Life Loss[J].Power System Technology, 2014, 38(9): 2371-2378. in Chinese)

[9]陳民武.電氣化鐵路電能質(zhì)量預(yù)測(cè)與對(duì)策分析研究[R].西安：中鐵第一勘察設(shè)計(jì)院集團(tuán)有限公司，2012.

[10]陳民武, 尚國(guó)旭, 智慧, 等. 高速鐵路牽引變壓器容量與配置方案優(yōu)化研究[J].中國(guó)鐵道科學(xué), 2013, 34(5):70-75.

(CHEN Minwu, SHANG Guoxu, ZHI Hui, et al. Optimization on the Capacity and Configuration Scheme for Traction Transformer of High Speed Railway [J]. China Railway Science, 2013, 34(5):70-75. in Chinese)

[11]中華人民共和國(guó)鐵道部.TB 10009—2005 鐵路電力牽引供電設(shè)計(jì)規(guī)范[S].北京：中國(guó)鐵道出版社，2005.

[12]劉文穎，謝昶，文晶，等.基于小生境多目標(biāo)粒子群算法的輸電網(wǎng)檢修計(jì)劃優(yōu)化[J].中國(guó)電機(jī)工程學(xué)報(bào)，2013,33(4)：141-148.

(LIU Wenying，XIE Chang，WEN Jing，et al.Optimization of Transmission Network Maintenance Scheduling Based on Niche Multi-Objective Particle Swarm Algorithm[J].Proceedings of the CSEE，2013，33(4):141-148. in Chinese)

[13]趙晶晶，李新，彭怡，等.基于粒子群優(yōu)化算法的配電網(wǎng)重構(gòu)和分布式電源注入功率綜合優(yōu)化算法[J].電網(wǎng)技術(shù)，2009，33(17)：162-166.

(ZHAO Jingjing，LI Xin，PENG Yi，et al.A Comprehensive Optimization Algorithm for Injection Power of Distributed Generation and Distribution Network Reconfiguration Based on Particle Swarm Optimization[J].Power System Technology，2009，33(17): 162-166. in Chinese)

[14]胡旺，GARY G YEN，張?chǎng)?基于Pareto熵的多目標(biāo)粒子群優(yōu)化算法[J].軟件學(xué)報(bào)， 2014，25(5)：1025-1050.

(HU Wang，GARY G YEN，ZHANG Xin.Multiobjective Particle Swarm Optimization Based on Pareto Entropy[J].Journal of Software，2014，25(5): 1025-1050. in Chinese)

[15]周孝法，陳陳，楊帆，等.基于自適應(yīng)混沌粒子群優(yōu)化算法的多饋入直流輸電系統(tǒng)優(yōu)化協(xié)調(diào)直流調(diào)制[J].電工技術(shù)學(xué)報(bào)，2009，24(4): 193-201.

(ZHOU Xiaofa，CHEN Chen，YANG Fan，et al.Optimal Coordinated HVDC Modulation Based on Adaptive Chaos Particle Swarm Optimization Algorithm in Multi-Infeed HVDC Transmission System[J].Transactions of China Electrotechnical Society，2009，24(4): 193-201. in Chinese)