辛康康，鐘建偉*，戴小劍，田家俊，龍玉雪，黃謀甫

(1.湖北民族大學(xué) 信息工程學(xué)院，湖北 恩施 445000；2.國網(wǎng)湖北省電力有限公司恩施供電公司，湖北 恩施 445000)

現(xiàn)代配電網(wǎng)(distributed network,DN)負(fù)荷供給需要不斷填補(bǔ)，導(dǎo)致DN運(yùn)行在負(fù)荷限度的邊緣[1-3].由此人們開始考量分布式電源(distributed generation,DG)的開辟層面，但當(dāng)其接到DN時(shí)會引發(fā)像電壓畸變等很多不利變化[4-7]，因此采用DN重構(gòu)來研究這一問題.該優(yōu)化工具有助于改變網(wǎng)絡(luò)配置，并利用網(wǎng)絡(luò)中的最優(yōu)潮流來改進(jìn)其功能.

DN重構(gòu)[8]是改變DN中某些分支的開關(guān)狀態(tài)，優(yōu)化某些DN參數(shù)，從而提高系統(tǒng)性能.本文DN重構(gòu)中網(wǎng)損的計(jì)算關(guān)聯(lián)到潮流軟件MATPOWER[9]，該軟件操作簡便、計(jì)算結(jié)果精確直觀，非常適用于本文模型.

1 DG潮流模型[10-13]

1.1 PQ型

此類型DG的P、Q不變，可看作負(fù)的負(fù)荷，體現(xiàn)在公式中即為將正號改變?yōu)樨?fù)號，即：

(1)

式中：Pi、Qi為DG等效處理后的功率；PDG、QDG為DG實(shí)際接入的功率.

1.2 PV型

此類型DG的P、V不變，即：

(2)

式中：Pi、Vi為DG等效處理后的有功功率及電壓；PDG、VDG為DG實(shí)際接入的有功功率及電壓.但由于計(jì)算只能用PQ型，解決要領(lǐng)為：

(3)

式中：t為更新數(shù)，ΔQ為無功變化量，Qt為DG等效處理后的無功功率.

1.3 PI型

此類型DG的P、I不變，即：

(4)

式中：Pi、Ii為DG等效處理后的有功功率及電流；PDG、IDG為DG實(shí)際接入的有功功率及電壓.但由于計(jì)算只能用PQ型，解決要領(lǐng)為:采用電壓迭代式.

(5)

1.4 P維系一致而Q因U而不同型

此類型DG即：

Im≤Immax.

(6)

式中：t為更新數(shù)，Ut-1為第t-1次更新電壓的最大絕對值.

2 MATPOWER闡釋

MATPOWER是以MATLAB為底層基礎(chǔ)的，聚焦于最佳潮流[14].根據(jù)功率損耗與潮流計(jì)算的關(guān)系，利用MATPOWER包來計(jì)算潮流.文獻(xiàn)[15]通過MATPOWER計(jì)算某地區(qū)實(shí)際電網(wǎng)潮流，得出MATPOWER操作簡便、計(jì)算精確、簡明直觀.MATPOWER的過程包括三個(gè)步驟.

步驟1：載入數(shù)值，將電力系統(tǒng)數(shù)據(jù)如母線數(shù)據(jù)、支路數(shù)據(jù)、發(fā)電機(jī)燃油成本函數(shù)等進(jìn)行讀取，再到下一個(gè)步驟.

步驟2：MATPOWER的依據(jù)是牛拉法.

步驟3：顯示結(jié)論.

在步驟1中，輸進(jìn)數(shù)據(jù)，它在一組矩陣包中確定為“MATPOWER case”結(jié)構(gòu)，這種結(jié)構(gòu)通常在示例中定義.MATPOWER適用于低壓配電網(wǎng)理論線損的計(jì)算.

3 模型構(gòu)造

3.1 目標(biāo)函數(shù)

網(wǎng)絡(luò)耗損最低，即：

(7)

式中，Pm、Qm是在DG連入之后記值的支路功率，Nm為支路個(gè)數(shù)，Rm為支路阻值，Um為節(jié)點(diǎn)電壓.

3.2 約束條件

1) 潮流約束：

(8)

式中，Pm、Qm為節(jié)點(diǎn)m處原始功率，PDGm、QDGm為DG給節(jié)點(diǎn)m處提供的功率，PLm、QLm為節(jié)點(diǎn)m處負(fù)載消耗的功率，Um、Un為節(jié)點(diǎn)m、n的電壓，Gmn、Bmn、δmn為節(jié)點(diǎn)m、n的電導(dǎo)值、電納值及相角差.

2) 電壓約束：

Ummin≤Um≤Ummax.

(9)

式中，Um為m處電壓值，Ummin、Ummax為m處電壓能承載的最小、最大值.

3) 電流約束：

Im≤Immax.

(10)

式中，Im為節(jié)點(diǎn)m處電流值，Immax為節(jié)點(diǎn)m處能承受電流的最大值.

圖1 算法流程圖Fig.1 Algorithm flow chart

圖2 節(jié)點(diǎn)連接圖Fig.2 Node connection diagram

3.3 模擬植物生長算法[16-18]流程

第1步：設(shè)置DG未連前的全部節(jié)點(diǎn)為根X0，并將其設(shè)置為最優(yōu)解Ymin和基點(diǎn)Y0.

第2步：導(dǎo)入原始支路和節(jié)點(diǎn)數(shù)據(jù)，利用MATPOWER軟件計(jì)算潮流，得到未接入DG前的網(wǎng)絡(luò)損耗及節(jié)點(diǎn)電壓分布情況，求出X0的目標(biāo)函數(shù)值Z0，并將其設(shè)置為最優(yōu)解Zmin.

第3步：迭代數(shù)設(shè)置為iter=50，步長s=5，q=s依次遞減1.

第4步：從Y0開始，在滿足約束的情況下查找新生長點(diǎn)Ynew，去除不滿足的點(diǎn)，保存滿足的點(diǎn)到Y(jié)中.

第5步：求取Y中生長點(diǎn)的Z并與Z0比照，將低于Z0的生長點(diǎn)設(shè)置為新的Ymin和Zmin，將等于Z0的生長點(diǎn)納入生長點(diǎn)聚合中，去除Y中相鄰次數(shù)生長過后Z值相同的生長點(diǎn).

第6步：求Y中各點(diǎn)的濃度值cm，由系統(tǒng)相關(guān)函數(shù)生成1以內(nèi)的任意數(shù)beta，由此確定Y中哪個(gè)生長點(diǎn)先生長.

第7步：如果目標(biāo)函數(shù)值Zmin反復(fù)出現(xiàn)，就終止更迭顯示結(jié)論，如果沒有的話，則再進(jìn)入第4步.

如圖1為該算法更迭流程圖.

4 算例仿真

如圖2所示，采用的是IEEE33模型，1為電源點(diǎn)，其他節(jié)點(diǎn)均為PQ節(jié)點(diǎn)，包含32條分段開關(guān)和32支路，虛線部分為聯(lián)絡(luò)開關(guān)，其序號為34～38.

本文擬在該模型中某點(diǎn)處接入一臺DG，使得線路的網(wǎng)損降到最小.采用MATPOWER軟件進(jìn)行該網(wǎng)絡(luò)的潮流計(jì)算，得到初始網(wǎng)損為211.270 3 kW，再采用基于模擬植物生長算法的配電網(wǎng)重構(gòu)模型對該系統(tǒng)進(jìn)行仿真，結(jié)果如圖3所示.

由圖可知，該模型在迭代2～3次后就完成了配電網(wǎng)重構(gòu)，重構(gòu)后網(wǎng)損降低到了75.278 6 kW，降低了64%，接入DG的位置為6節(jié)點(diǎn)處，其容量為2.526 7 kW，表明DG接入DN的合適地方和大小會降低整個(gè)網(wǎng)絡(luò)的損失大小.圖4為加和不加DG各點(diǎn)對照圖.

由圖4可知，經(jīng)過DN重構(gòu)后該系統(tǒng)的所有節(jié)點(diǎn)的電壓都得到了一定的提升，表明含DG的DN重構(gòu)能有效地改善網(wǎng)絡(luò)電壓質(zhì)量，抬高網(wǎng)絡(luò)最低電壓值;同時(shí)6處的電壓也有一定的提升，說明DG的接入會對其連接點(diǎn)的電壓起到抬高功能.

圖3 運(yùn)作數(shù)與丟損數(shù)圖 圖4 加/不加分布式電源各點(diǎn)對照圖Fig.3 Operation number and lost number graph Fig.4 Voltage profile with & without DG

5 結(jié)語

中MATPOWER被用作相關(guān)計(jì)算，使用簡單快速且結(jié)果精準(zhǔn)直觀，并構(gòu)建了采用PGSA的DN重構(gòu)構(gòu)型做仿真.最后得出：DG接入DN的合適地方和大小會降低整個(gè)網(wǎng)絡(luò)的網(wǎng)損大小，改善網(wǎng)絡(luò)電壓的質(zhì)量，同時(shí)DG的接入會對其連接點(diǎn)的電壓起到提升作用.