陳民武,許臣友,黃文勛,邱 璐,孫小凱
(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)化模型和算法的正確性及有效性。
基于多導(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]。
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ì)方案排除。
在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.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ò)大搜索空間,以避免陷入局部極值。 將基于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為外部檔案最大容量。 以某高速鐵路牽引供電系統(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ì)性。 以牽引供電系統(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)2.2 自適應(yīng)混沌粒子群優(yōu)化算法
3 基于改進(jìn)粒子群算法的牽引供電系統(tǒng)多目標(biāo)優(yōu)化過(guò)程
4 算例分析
5 結(jié) 語(yǔ)