【正文】 條線路只是輸送發(fā)電機(jī)組出力的一部分，它的系數(shù)就不會(huì)大于1，也不可能大于1。（二）、阻塞費(fèi)用計(jì)算規(guī)則的設(shè)計(jì)：對(duì)阻塞費(fèi)用規(guī)則的設(shè)計(jì)，既要考慮電力市場(chǎng)交易規(guī)則，又要公平對(duì)待序內(nèi)容量不能出力的部分和報(bào)價(jià)高于清算價(jià)的序外容量出力的部分。 序內(nèi)容量：通過競(jìng)價(jià)取得發(fā)電權(quán)的發(fā)電容量；u序內(nèi)容量和序外容量最小損失的計(jì)算根據(jù)市場(chǎng)交易規(guī)則，全部機(jī)組的所有出力都是按清算價(jià)結(jié)算，因此，調(diào)整后序內(nèi)容量和序外容量對(duì)應(yīng)的機(jī)組都會(huì)有經(jīng)濟(jì)損失。由于報(bào)價(jià)高于清算價(jià)，機(jī)組在以清算價(jià)賣給電網(wǎng)公司的出力中就有經(jīng)濟(jì)損失。調(diào)整后，序外容量不能出力的部分，電網(wǎng)公司按清算價(jià)結(jié)算，并將資金付給出力的機(jī)組，因此該部分的損失只在于沒能將其按高于報(bào)價(jià)的清算價(jià)出力；而序內(nèi)容量不出力的部分，電網(wǎng)公司將不再以清算價(jià)結(jié)算，該部分的損失只在于以清算價(jià)賣出時(shí)能賺的最小利潤(rùn)。在此基礎(chǔ)上，建立起各機(jī)組在下一時(shí)段的出力分配預(yù)案模型。于是得到下列模型：（）其中：C為給定值，和為在表3和4中給出的定值。此值就是機(jī)組i在下一時(shí)段出力值的下限，不可能更小。（d）、最后選取的段的報(bào)價(jià)就是該時(shí)段所有機(jī)組出力的清算價(jià)。u u利用問題(1)得到的各線路上有功潮流的計(jì)算表達(dá)式，可以求出各線路上的有功潮流值，然后與潮流限值比較，若大于限值，說明這條線路上會(huì)發(fā)生輸電阻塞。 調(diào)整各機(jī)組出力分配方案使得輸電阻塞消除。 如果（1）做不到，還可以使用線路的安全裕度輸電，以避免拉閘限電（強(qiáng)制減少負(fù)荷需求），但要使每條線路上潮流的絕對(duì)值超過限值的百分比盡量小。 如果無論怎樣分配機(jī)組出力都無法使每條線路上的潮流絕對(duì)值超過限值的百分比小于相對(duì)安全裕度，則必須在用電側(cè)拉閘限電。若a=()那么調(diào)整后機(jī)組i的第k段有a的發(fā)電容量,在分配預(yù)案中取得了發(fā)電權(quán)但不能出力，即，為序內(nèi)容量。又因?yàn)?，所以0。如果對(duì)所給的負(fù)荷要求。輸電的安全和經(jīng)濟(jì)兩個(gè)原則不可能同時(shí)達(dá)到最優(yōu)。安全原則要求每條線路上潮流的絕對(duì)值超過限值的百分比盡量小，首先應(yīng)該小于線路的安全預(yù)度，即： 然后可以認(rèn)為，如果超過限值的百分比最大的線路安全問題讓人滿意，那么其余線路超出限值得百分比均可滿意，即要求每條線路上潮流的絕對(duì)值超過限值的百分比都小于某個(gè)數(shù)(0),表示為：目標(biāo)函數(shù)即是求的最小值： min 而且要去除各線路潮流值小于限值的約束。因?yàn)檩旊娋€路在超負(fù)荷下使用已經(jīng)在所難免，優(yōu)化目標(biāo)就變成使限制的電量最小，因?yàn)榫W(wǎng)方應(yīng)該盡量滿足用電方的用電要求，體現(xiàn)用戶至上的原則。計(jì)算得到對(duì)照表：下面根據(jù)建立的機(jī)組出力分配方案的調(diào)整模型對(duì)兩個(gè)分配預(yù)案進(jìn)行調(diào)整。用Lingo程序（見附件[2]）求解，因?yàn)槎稳萘繛?的對(duì)各機(jī)組出力的調(diào)整毫無影響，所以在求解時(shí)可以去掉對(duì)應(yīng)的變量，以減少變量個(gè)數(shù)。（2）表六 調(diào)整后各線路潮流絕對(duì)值超過限值百分比線路123456調(diào)整后%000%%對(duì)求解結(jié)果的改進(jìn)可以發(fā)現(xiàn)其中線路1已經(jīng)達(dá)到了超出限值的百分比的最小值，且線路5超出限值的百分比相當(dāng)小。在求解時(shí)考慮到有正有負(fù)，選取松弛變量進(jìn)行代換，求和是只求所有超出限值的百分比，最后得到調(diào)整后的各線路潮流絕對(duì)值超過限值的百分比，見表八：表九 調(diào)整后的分配方案機(jī)組12345678總出力值（）228152155117八、 結(jié)果分析，雖然可以使得輸電阻塞消除，但是阻塞費(fèi)用很大。最后求得6條線路中僅有1條超出潮流值的限值，這非常符合輸電線路的安全原則，且阻塞費(fèi)用不是很大。九、 模型評(píng)價(jià)首先建立了各機(jī)組出力分配預(yù)案模型，嚴(yán)格遵循電力市場(chǎng)交易規(guī)則，考慮爬坡速度的限制，先滿足各機(jī)組出力的下限值，然后按段價(jià)從低到高選取未選滿的段容量。依次考慮輸電阻塞管理原則中的三種情況建立了三個(gè)模型來求解不同的負(fù)荷要求。這樣求解可以更接近最優(yōu)值。還可以考慮，先分析各線路上潮流值的關(guān)于分配方案的靈敏度，然后以靈敏度大的潮流值的絕對(duì)值超出限值的百分比為限值條件，阻塞費(fèi)用最小為目標(biāo)函數(shù)，進(jìn)行求解。 參考文獻(xiàn)[1].Matlab程序：給定預(yù)報(bào)負(fù)荷需求時(shí)，求各機(jī)組出力的分配預(yù)案[2].lingo程序：，求解調(diào)整后的方案[3].lingo程序：，用原模型求解調(diào)整后的方案[4].lingo程序：，用改后模型求解調(diào)整后的方案附件附件[1]：給定預(yù)報(bào)負(fù)荷需求時(shí)，求各機(jī)組出力的分配預(yù)案的程序load load load load lowlimit=currentspeed.*15。segment=zeros(1,length(current))。)。 while next(i)lowlimit(i) j=j+1。 nextprice(i)=price(i,j)。endminprice=0。 next(j)=next(j)+max。 sum=destination。 [minprice,j]=min(nextprice)。segment=segment+ones(1,length(segment))。 currentprice=nextprice(j)。 next(j)=uplimit(j)。 next(j)=uplimit(j)。 else nextprice(j)=inf。505 0 124 168 210 252 312 330 363 489560 0 182 203 245 300 320 360 410 495610 0 152 189 233 258 308 356 415 500500 150 170 200 255 302 325 380 435 800590 0 116 146 188 215 250 310 396 510607 0 159 173 205 252 305 380 405 520500 120 180 251 260 306 315 335 348 548800 153 183 233 253 283 303 318 400 800x11040。x2515。x31110。x3840。x4310。x4715。x5315。x5910。x6420。x61010。x7510。x793。x8520。x8105。*(x11+x13+x16+x110)*(x21+x23+x24+x25+x26+x27+x210)+*(x31+x33+x35+x37+x38+x310)*(x51+x52+x53+x55+x56+x58+x59+x510)+*(x81+x83+x85+x87+x88+x89+x810)。(x11+x13+x16+x110)+(x21+x23+x24+x25+x26+x27+x210)+(x31+x33+x35+x37+x38+x310)+(x41+x42+x43+x44+x45+x46+x47+x410)+(x51+x52+x53+x55+x56+x58+x59+x510)+(x61+x63+x64+x66+x67+x68+x610)+(x71+x72+x73+x74+x75+x76+x77+x78+x79+x710)+(x81+x83+x85+x87+x88+x89+x810)=。x21+x23+x24+x25+x26+x27+x21058。x41+x42+x43+x44+x45+x46+x47+x410。x61+x63+x64+x66+x67+x68+x61095。x81+x83+x85+x87+x88+x89+x81063。(x11+x13+x16+x110)+(x21+x23+x24+x25+x26+x27+x210)+(x31+x33+x35+x37+x38+x310)+(x41+x42+x43+x44+x45+x46+x47+x410)+(x51+x52+x53+x55+x56+x58+x59+x510)+(x61+x63+x64+x66+x67+x68+x610)+(x71+x72+x73+x74+x75+x76+x77+x78+x79+x710)+(x81+x83+x85+x87+x88+x89+x810)=。x21+x23+x24+x25+x26+x27+x21088。x41+x42+x43+x44+x45+x46+x47+x410。x61+x63+x64+x66+x67+x68+x610155。x81+x83+x85+x87+x88+x89+x810117。x1630。x248。x2108。x3720。x425。x4610。x525。x5810。x6310。x6820。x7415。x7810。x8320。x8915。p1p10=(*(x11+x13+x16+x110)+*(x21+x23+x24+x25+x26+x27+x210)+*(x31+x33+x35+x37+x38+x310)+*(x41+x42+x43+x44+x45+x46+x47+x410)*(x51+x52+x53+x55+x56+x58+x59+x510)+*(x61+x63+x64+x66+x67+x68+x610)+*(x71+x72+x73+x74+x75+x76+x77+x78+x79+x710)+)/165。p5p50=(*(x21+x23+x24+x25+x26+x27+x210)*(x31+x33+x35+x37+x38+x310)*(x41+x42+x43+x44+x45+x46+x47+x410)*(x51+x52+x53+x55+x56+x58+x59+x510)+*(x61+x63+x64+x66+x67+x68+x610)+)/132。pp3pp30=((p3p30)+(pp2pp20)+abs((p3p30)(pp2pp20)))/2。p4p400。附件[4]：，用改后模型求解調(diào)整后方案的程序model:min=p1+p2+p3+p4+p5+p6。x21+x23+x24+x25+x26+x27+x21058。x41+x42+x43+x44+x45+x46+x47+x410。x61+x63+x64+x66+x67+x68+x61095。x81+x83+x85+x87+x88+x89+x81063。x1350。x2320。x272。x3530。x4155。x4510。x5175。x5615。x6195。x6710。x7x735。x775。x8170。x8810。p4p40=(*(x11+x13+x16+x110)*(x21+x23+x24+x25+x26+x27+x210)+*(x31+x33+x35+x37+x38+x310)*(x41+x42+x43+x44+x45+x46+x47+x410)*(x51+x52+x53+x55+x56+x58+x59+x510)+*(x71+x72+x73+x74+x75+x76+x77+x78+x79+x710)+*(x81+x83+x85+x87+x88+x89+x810)+)/155。pp2pp20=((p1p10)+(p2p20)+abs((p1p10)(p2p20)))/2。pp6pp60=((p6p60)+(pp5pp50)+abs((p6p60)(pp5pp50)))