【正文】 條線路只是輸送發(fā)電機組出力的一部分，它的系數(shù)就不會大于1，也不可能大于1。（二）、阻塞費用計算規(guī)則的設計：對阻塞費用規(guī)則的設計，既要考慮電力市場交易規(guī)則，又要公平對待序內(nèi)容量不能出力的部分和報價高于清算價的序外容量出力的部分。 序內(nèi)容量：通過競價取得發(fā)電權的發(fā)電容量；u序內(nèi)容量和序外容量最小損失的計算根據(jù)市場交易規(guī)則，全部機組的所有出力都是按清算價結算，因此，調(diào)整后序內(nèi)容量和序外容量對應的機組都會有經(jīng)濟損失。由于報價高于清算價，機組在以清算價賣給電網(wǎng)公司的出力中就有經(jīng)濟損失。調(diào)整后，序外容量不能出力的部分，電網(wǎng)公司按清算價結算，并將資金付給出力的機組，因此該部分的損失只在于沒能將其按高于報價的清算價出力；而序內(nèi)容量不出力的部分，電網(wǎng)公司將不再以清算價結算，該部分的損失只在于以清算價賣出時能賺的最小利潤。在此基礎上，建立起各機組在下一時段的出力分配預案模型。于是得到下列模型：（）其中：C為給定值，和為在表3和4中給出的定值。此值就是機組i在下一時段出力值的下限，不可能更小。（d）、最后選取的段的報價就是該時段所有機組出力的清算價。u u利用問題(1)得到的各線路上有功潮流的計算表達式，可以求出各線路上的有功潮流值，然后與潮流限值比較，若大于限值，說明這條線路上會發(fā)生輸電阻塞。 調(diào)整各機組出力分配方案使得輸電阻塞消除。 如果（1）做不到，還可以使用線路的安全裕度輸電，以避免拉閘限電（強制減少負荷需求），但要使每條線路上潮流的絕對值超過限值的百分比盡量小。 如果無論怎樣分配機組出力都無法使每條線路上的潮流絕對值超過限值的百分比小于相對安全裕度，則必須在用電側拉閘限電。若a=()那么調(diào)整后機組i的第k段有a的發(fā)電容量,在分配預案中取得了發(fā)電權但不能出力，即，為序內(nèi)容量。又因為0，所以0。如果對所給的負荷要求。輸電的安全和經(jīng)濟兩個原則不可能同時達到最優(yōu)。安全原則要求每條線路上潮流的絕對值超過限值的百分比盡量小，首先應該小于線路的安全預度，即： 然后可以認為，如果超過限值的百分比最大的線路安全問題讓人滿意，那么其余線路超出限值得百分比均可滿意，即要求每條線路上潮流的絕對值超過限值的百分比都小于某個數(shù)(0),表示為：目標函數(shù)即是求的最小值： min 而且要去除各線路潮流值小于限值的約束。因為輸電線路在超負荷下使用已經(jīng)在所難免，優(yōu)化目標就變成使限制的電量最小，因為網(wǎng)方應該盡量滿足用電方的用電要求，體現(xiàn)用戶至上的原則。計算得到對照表：下面根據(jù)建立的機組出力分配方案的調(diào)整模型對兩個分配預案進行調(diào)整。用Lingo程序（見附件[2]）求解，因為段容量為0的對各機組出力的調(diào)整毫無影響，所以在求解時可以去掉對應的變量，以減少變量個數(shù)。（2）表六 調(diào)整后各線路潮流絕對值超過限值百分比線路123456調(diào)整后%000%%對求解結果的改進可以發(fā)現(xiàn)其中線路1已經(jīng)達到了超出限值的百分比的最小值，且線路5超出限值的百分比相當小。在求解時考慮到有正有負，選取松弛變量進行代換，求和是只求所有超出限值的百分比，最后得到調(diào)整后的各線路潮流絕對值超過限值的百分比，見表八：表九 調(diào)整后的分配方案機組12345678總出力值（）228152155117八、 結果分析，雖然可以使得輸電阻塞消除，但是阻塞費用很大。最后求得6條線路中僅有1條超出潮流值的限值，這非常符合輸電線路的安全原則，且阻塞費用不是很大。九、 模型評價首先建立了各機組出力分配預案模型，嚴格遵循電力市場交易規(guī)則，考慮爬坡速度的限制，先滿足各機組出力的下限值，然后按段價從低到高選取未選滿的段容量。依次考慮輸電阻塞管理原則中的三種情況建立了三個模型來求解不同的負荷要求。這樣求解可以更接近最優(yōu)值。還可以考慮，先分析各線路上潮流值的關于分配方案的靈敏度，然后以靈敏度大的潮流值的絕對值超出限值的百分比為限值條件，阻塞費用最小為目標函數(shù)，進行求解。 參考文獻[1].Matlab程序：給定預報負荷需求時，求各機組出力的分配預案[2].lingo程序：，求解調(diào)整后的方案[3].lingo程序：，用原模型求解調(diào)整后的方案[4].lingo程序：，用改后模型求解調(diào)整后的方案附件附件[1]：給定預報負荷需求時，求各機組出力的分配預案的程序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)))