電力市場(chǎng)的輸電阻塞管理論文

.z.論文2電力市場(chǎng)的輸電阻塞管理參賽隊(duì)員：、、摘要：本文根據(jù)電力市場(chǎng)的交易規(guī)則,就目前我國電力系統(tǒng)中各個(gè)發(fā)電機(jī)組的出力分配預(yù)案和各線路的有功潮流問題進(jìn)行了深入分析,并對(duì)產(chǎn)生輸電阻塞的分配預(yù)案進(jìn)行了調(diào)整,得到了較好的出力分配方案。1.根據(jù)各機(jī)組出力和各線路潮流的關(guān)系建立了一個(gè)多元線性回歸模型(見模型一),利用所給實(shí)驗(yàn)數(shù)據(jù)采用最小二乘法回歸，得到每條線路上的潮流值關(guān)于各發(fā)電機(jī)組出力的的近似表達(dá)式，并對(duì)每一個(gè)表達(dá)式進(jìn)行了誤差分析，得出各表達(dá)式的復(fù)相關(guān)系數(shù),可以看出我們的回歸效果顯著，說明我們的模型是可靠、合理的。2.我們采用pool模式下的輸電阻塞費(fèi)用計(jì)算方法，公平對(duì)待序內(nèi)序外兩種情況,設(shè)計(jì)出了一種簡明、合理的阻塞費(fèi)用計(jì)算規(guī)則:第一、采用序外多發(fā)電量按照發(fā)電報(bào)價(jià)計(jì)算；第二、序內(nèi)少發(fā)電量按清算價(jià)與發(fā)電報(bào)價(jià)之差價(jià)結(jié)算。并建立了一個(gè)合理的計(jì)算阻塞費(fèi)用模型。3．在下一時(shí)段預(yù)報(bào)負(fù)荷需求為982.4MW的條件下,根據(jù)市場(chǎng)規(guī)則，以最小購電費(fèi)用為目標(biāo)、以機(jī)組的段容量,爬坡速率作為約束條件,采用動(dòng)態(tài)規(guī)劃算法建立了一個(gè)單目標(biāo)規(guī)劃模型，通過數(shù)學(xué)軟件MATLAB編程給出各機(jī)組的出力分配預(yù)案，各臺(tái)機(jī)組的出力分別為（MW）：150、79、180、99.5、125、140、95、113.9?？傎M(fèi)用為：=74416.8元。清算價(jià)為:元/MWh4.通過對(duì)預(yù)案分析計(jì)算可得,第一、五、六線路出現(xiàn)輸電阻塞現(xiàn)象,根據(jù)安全且經(jīng)濟(jì)的原則，利用排序算法進(jìn)行了調(diào)整，得到了消除輸電阻塞的分配方案，分別是：152.5947,88.0,228.0,81.7434,152.0,98.0,65.0592,117.0。其清算價(jià)格=303元/MWh阻塞費(fèi)用：12026.0元；總費(fèi)用為：=86140.8元5.同理對(duì)下一時(shí)段預(yù)報(bào)負(fù)荷需求為1052.8MW的條件下，重復(fù)步驟3、4的工作，得到分配預(yù)案為(MW)：152.0,81.0,218.2,99.5,135.0,150.0,102.1,117.0.總費(fèi)用為：=93699.2元。清算價(jià)為:元/MWh;通過調(diào)整預(yù)案不能消除阻塞，然后采用輸電阻塞管理原則第二條，得到新的方案：149.8,88.0,228.2,99.5,152.0,128.0,90.3,117.0。阻塞費(fèi)用為：元;總費(fèi)用為：=94924.8元.最后，對(duì)所得結(jié)果進(jìn)行了詳細(xì)的分析、評(píng)價(jià)和推廣。本文所建模型正確，且建模的思路清晰，簡單易懂，能最大限度地實(shí)現(xiàn)資源優(yōu)化配置的獨(dú)特優(yōu)點(diǎn)，具有較強(qiáng)的實(shí)用使用價(jià)值，對(duì)網(wǎng)方和市場(chǎng)交易－調(diào)度中心的實(shí)時(shí)調(diào)度具有一定的指導(dǎo)意義。-.z.一、問題重述我國電力市場(chǎng)初期是發(fā)電側(cè)電力市場(chǎng)，采取交易與調(diào)度一體化的模式。電網(wǎng)公司在組織交易、調(diào)度和配送時(shí)，必須遵循電網(wǎng)“安全第一”的原則，同時(shí)要制訂一個(gè)電力市場(chǎng)交易規(guī)則，按照購電費(fèi)用最小的經(jīng)濟(jì)目標(biāo)來運(yùn)作。市場(chǎng)交易-調(diào)度中心根據(jù)負(fù)荷預(yù)報(bào)和交易規(guī)則制訂滿足電網(wǎng)安全運(yùn)行的調(diào)度計(jì)劃――各發(fā)電機(jī)組的出力（發(fā)電功率）分配方案；在執(zhí)行調(diào)度計(jì)劃的過程中，還需實(shí)時(shí)調(diào)度承擔(dān)AGC（自動(dòng)發(fā)電控制）輔助服務(wù)的機(jī)組出力，以跟蹤電網(wǎng)中實(shí)時(shí)變化的負(fù)荷。設(shè)*電網(wǎng)有若干臺(tái)發(fā)電機(jī)組和若干條主要線路，每條線路上的有功潮流（輸電功率和方向）取決于電網(wǎng)結(jié)構(gòu)和各發(fā)電機(jī)組的出力。電網(wǎng)每條線路上的有功潮流的絕對(duì)值有一安全限值，限值還具有一定的相對(duì)安全裕度（即在應(yīng)急情況下潮流絕對(duì)值可以超過限值的百分比的上限）。如果各機(jī)組出力分配方案使*條線路上的有功潮流的絕對(duì)值超出限值，稱為輸電阻塞。當(dāng)發(fā)生輸電阻塞時(shí)，需要研究如何制訂既安全又經(jīng)濟(jì)的調(diào)度計(jì)劃。根據(jù)市場(chǎng)規(guī)則（見附錄1），考慮以下問題（表（一）～表（六）見附錄2）：*電網(wǎng)有8臺(tái)發(fā)電機(jī)組，6條主要線路，表1和表2中的方案0給出了各機(jī)組的當(dāng)前出力和各線路上對(duì)應(yīng)的有功潮流值，方案1~32給出了圍繞方案0的一些實(shí)驗(yàn)數(shù)據(jù)，試用這些數(shù)據(jù)確定各線路上有功潮流關(guān)于各發(fā)電機(jī)組出力的近似表達(dá)式。設(shè)計(jì)一種簡明、合理的阻塞費(fèi)用計(jì)算規(guī)則，除考慮上述電力市場(chǎng)規(guī)則外，還需注意：在輸電阻塞發(fā)生時(shí)公平地對(duì)待序內(nèi)容量不能出力的部分和報(bào)價(jià)高于清算價(jià)的序外容量出力的部分。假設(shè)下一個(gè)時(shí)段預(yù)報(bào)的負(fù)荷需求是982.4MW，表3、表4和表5分別給出了各機(jī)組的段容量、段價(jià)和爬坡速率的數(shù)據(jù)，試按照電力市場(chǎng)規(guī)則給出下一個(gè)時(shí)段各機(jī)組的出力分配預(yù)案。按照表6給出的潮流限值，檢查得到的出力分配預(yù)案是否會(huì)引起輸電阻塞，并在發(fā)生輸電阻塞時(shí)，根據(jù)安全且經(jīng)濟(jì)的原則，調(diào)整各機(jī)組出力分配方案，并給出與該方案相應(yīng)的阻塞費(fèi)用。假設(shè)下一個(gè)時(shí)段預(yù)報(bào)的負(fù)荷需求是1052.8MW，重復(fù)3~4的工作。二、基本假設(shè)1.電力從生產(chǎn)到使用的四大環(huán)節(jié)——發(fā)電、輸電、配電和用電是瞬間完成的；2.每臺(tái)機(jī)組均不停機(jī)；3.電網(wǎng)公司在組織交易、調(diào)度和配送時(shí)，必須遵循電網(wǎng)“安全第一”的原則，制訂電力市場(chǎng)交易規(guī)則時(shí)，按照購電費(fèi)用最小的經(jīng)濟(jì)目標(biāo)來運(yùn)作；4.沒有任何一方擁有支配市場(chǎng)運(yùn)作的能力，市場(chǎng)處于自由競爭中，參與市場(chǎng)發(fā)電的部門都是理性的；5.電網(wǎng)公司在組織交易、調(diào)度和配送電時(shí)，嚴(yán)格按照市場(chǎng)交易規(guī)則來運(yùn)作；6.在最低技術(shù)出力以下的報(bào)價(jià)為負(fù)值，表示愿意付費(fèi)維持發(fā)電以避免停機(jī)帶來更大的損失；7.當(dāng)潮流值為負(fù)值時(shí)，電流的通過網(wǎng)方向發(fā)電方向流動(dòng)；8.每條線路上的有功潮流的絕對(duì)值超出其限值時(shí)，就會(huì)發(fā)生輸電阻塞現(xiàn)象；9.當(dāng)線路發(fā)生阻塞時(shí)，按照輸電阻塞管理原則進(jìn)行調(diào)整輸電分配方案；10.電網(wǎng)結(jié)構(gòu)在短時(shí)間內(nèi)不會(huì)改變；11.以15分鐘作為一個(gè)時(shí)間段；12.當(dāng)潮流值剛好為限值時(shí)認(rèn)為沒有產(chǎn)生輸電阻塞；三、模型建立3.１名詞解釋輸電阻塞：當(dāng)線路上的有功潮流值超出其限值時(shí)，稱為輸電阻塞；爬坡速率：機(jī)組在單位時(shí)間內(nèi)所能夠增加或減少出力值稱為該機(jī)組的爬坡速率；限上：當(dāng)輸電發(fā)生阻塞時(shí)，有的機(jī)組會(huì)逆序開機(jī)或多發(fā)，稱為限上；限下：當(dāng)輸電發(fā)生阻塞時(shí)，有的機(jī)組會(huì)逆序關(guān)機(jī)或少發(fā)，稱為限下；3.3模型的建立發(fā)、輸、配電長期垂直壟斷我國電力產(chǎn)業(yè)，使得電力產(chǎn)業(yè)缺乏競爭、效率低下，不利于實(shí)現(xiàn)資源的優(yōu)化配置和降低成本，電力服務(wù)與價(jià)格均難以滿足我國經(jīng)濟(jì)快速發(fā)展的需求。自20世紀(jì)90年代以來，我國進(jìn)行了以提高效率、降低電價(jià)、引入競爭為目標(biāo)的電力工業(yè)改革，由于競爭機(jī)制的引入，就要求保證市場(chǎng)的自由度，而電力系統(tǒng)的基本特征——輸電阻塞和網(wǎng)損極大的限制了這種自由度，本問題就如何協(xié)調(diào)這之間的矛盾提出一些見解。給出計(jì)算框圖：計(jì)算線路潮計(jì)算線路潮流結(jié)束價(jià)格是否最優(yōu)是否滿足開始判斷是是結(jié)束價(jià)格是否最優(yōu)是否滿足開始否調(diào)整調(diào)整調(diào)整調(diào)整問題一各發(fā)電機(jī)組所發(fā)的電通過主要線路向網(wǎng)方供電，電力系統(tǒng)的特點(diǎn)是產(chǎn)、供、銷、用電同時(shí)完成，產(chǎn)品不能儲(chǔ)藏，電力的生產(chǎn)和消費(fèi)要保持平衡，也就是說發(fā)電量的多少是根據(jù)負(fù)荷需求來決定的。這就需要調(diào)度每臺(tái)機(jī)組出力的多少，可以通過一些實(shí)驗(yàn)值，進(jìn)行分析，找出其中的規(guī)律，根據(jù)這一規(guī)律可以分配出下一時(shí)段各機(jī)組的出力多少，因此我們以各機(jī)組出力多少為自變量、各線路的潮流值為因變量建立多元線性規(guī)劃模型，通過數(shù)學(xué)軟件MATLAB來實(shí)現(xiàn)模型的求解。經(jīng)過綜合分析，本問題采用多元線性回歸方法，假若*電網(wǎng)有臺(tái)發(fā)電機(jī)組和條主要線路，每臺(tái)發(fā)電機(jī)組的出力記為（），單位為MW，每條線路上的潮流值設(shè)為（），單位為MW，假定個(gè)自變量與因變量存在*種相關(guān)關(guān)系，可以設(shè)為：（）（MW） (1)其中（，），都為常數(shù),為各發(fā)電機(jī)組的系數(shù)。則式(1)可以進(jìn)一步可寫為：（MW） (2)式（2）即為第（）條線路關(guān)于各發(fā)電機(jī)組出力的多元線性規(guī)劃模型。問題二輸電阻塞是電力系統(tǒng)運(yùn)行中常見的問題，根據(jù)假設(shè)8，當(dāng)線路發(fā)生阻塞時(shí)，按照輸電阻塞管理原則進(jìn)行調(diào)整輸電分配方案，則有部分序內(nèi)容量不能出力，而序外容量要在低于對(duì)應(yīng)報(bào)價(jià)的清算價(jià)上出力，由于電力市場(chǎng)環(huán)境下各成員之間的關(guān)系是經(jīng)濟(jì)利益的關(guān)系，這對(duì)于序內(nèi)容量來說，有力而不出就會(huì)產(chǎn)生經(jīng)濟(jì)損失，而對(duì)于序外容量來說出了力反而還要虧本，這就使得發(fā)電商和網(wǎng)方產(chǎn)生利益沖突，在結(jié)算時(shí)適當(dāng)?shù)慕o發(fā)電方一定的補(bǔ)償，也就是所謂的阻塞費(fèi)用，產(chǎn)生的這一阻塞費(fèi)用如何來計(jì)算是解決本問題的關(guān)鍵，要解決這一問題，需要找到序外與序內(nèi)容量的一個(gè)平衡點(diǎn)，以使雙方均滿意，且公平合理。當(dāng)線路發(fā)生阻塞時(shí)，就必須對(duì)原來的分配方案進(jìn)行調(diào)整，設(shè)調(diào)整前第臺(tái)機(jī)組的出力為（單位：MW），，調(diào)整后第臺(tái)機(jī)組的出力為，調(diào)整前網(wǎng)方的總的購電費(fèi)用為（單位：元），則有： (3)其中（單位：元/MWh）表示清算價(jià)。若所采用的分配預(yù)案會(huì)產(chǎn)生輸電阻塞時(shí)，考慮電網(wǎng)約束，根據(jù)使總的購電費(fèi)用最小的原則，重新調(diào)整各機(jī)組的出力分配方案，調(diào)整后，根據(jù)電力市場(chǎng)規(guī)則，公平的對(duì)待序內(nèi)容量不能出力的部分和電價(jià)高于清算價(jià)的序外容量出力部分，擬采用序外多發(fā)電量按照發(fā)電報(bào)價(jià)計(jì)算，序內(nèi)少發(fā)電量按發(fā)電清算價(jià)與報(bào)價(jià)之差結(jié)算；則支付給多發(fā)電量的發(fā)電商的費(fèi)用（單位：元）為： (4)其中為第臺(tái)機(jī)組的第段容量的段價(jià)，單位為MW。支付少發(fā)電量的發(fā)電商費(fèi)用（單位：元）為： (5)其中表示調(diào)整后的第臺(tái)發(fā)電機(jī)組的出力。所以阻塞總費(fèi)用（單位：元）的模型為： (6)問題三根據(jù)假設(shè)9，當(dāng)每條線路上的有功潮流的絕對(duì)值超出其限值時(shí)，就會(huì)發(fā)生輸電阻塞現(xiàn)象，產(chǎn)生阻塞現(xiàn)象的主要原因是由于負(fù)荷過大，而成本相對(duì)較低的發(fā)電機(jī)組的產(chǎn)量就會(huì)加大，使得線路中的潮流值超出其安全限值，在一定量的需求下，為了盡量不產(chǎn)生輸電阻塞，就需要進(jìn)行合理的對(duì)各個(gè)機(jī)組分配出力預(yù)案，并且還要使得購電費(fèi)用最小，據(jù)此可建立一個(gè)規(guī)劃模型，模型以最小費(fèi)用為目標(biāo)函數(shù)，并且還要受到機(jī)組的爬坡速率、各線路的限值等條件的約束，從而求得最佳的分配預(yù)案。根據(jù)題中的要求，在按照市場(chǎng)規(guī)則進(jìn)行分配各機(jī)組出力時(shí)，給出的分配預(yù)案不僅要求費(fèi)用低，而且還要受負(fù)荷需求、各機(jī)組的爬坡速率、段容量等因素的影響。我們可以考慮成為一個(gè)優(yōu)化模型，以購電最小費(fèi)用為目標(biāo)函數(shù)，以負(fù)荷需求、爬坡速率為約束條件建立優(yōu)化模型，其目標(biāo)函數(shù)為：min (7)其中（單位：MW）表示負(fù)荷需求量，由于最后一個(gè)被選入的段價(jià)作為清算價(jià)，則清算價(jià)約束條件可寫為： (8)又因?yàn)榘l(fā)電與需求要保持平衡，所以負(fù)荷的約束條件可以寫為： (9)各機(jī)組的爬坡速率對(duì)段容量的約束可以寫為： (10)其中為初始時(shí)刻，即方案0的第臺(tái)機(jī)組的出力，為一個(gè)交易時(shí)段。由式(7)～式(10)可得機(jī)組分配預(yù)案的單目標(biāo)優(yōu)化數(shù)學(xué)模型為：minst. (11)問題四所給出力分配預(yù)案是否會(huì)引起輸電阻塞，就需計(jì)算出各線路上的潮流值，看是否超過其限值，這就要借助問題一中所得到的各線路中潮流值關(guān)于各機(jī)組出力的近似表達(dá)式，通過計(jì)算如果線路中的潮流值超過了限值，則將會(huì)發(fā)生阻塞，反之電網(wǎng)運(yùn)行是安全的，當(dāng)線路發(fā)生阻塞時(shí)，就需要對(duì)分配預(yù)案進(jìn)行調(diào)整，使得輸電阻塞消除，如果仍然不能夠消除輸電阻塞，就可以考慮在安全裕度內(nèi)輸電，盡量不限電拉閘，使損失盡量減少。在出現(xiàn)輸電阻塞的情況下，我們要對(duì)分配預(yù)案進(jìn)行調(diào)整。通過上面的分析，出現(xiàn)阻塞的原因是線路上的有功潮流的絕對(duì)值超過了限值。必須對(duì)原來的方案進(jìn)行調(diào)整,由于調(diào)整后會(huì)產(chǎn)生阻塞費(fèi)用,考慮安全且經(jīng)濟(jì)的原則,所以仍然要建立一個(gè)優(yōu)化模型,我們以阻塞費(fèi)用最小為目標(biāo),以各線路的潮流限值、各發(fā)電機(jī)組出力變化為約束條件建立一個(gè)優(yōu)化模型,由式(6)可知,其最小阻塞費(fèi)用為: (12)設(shè)第條線路的潮流限值為,每條線路上的潮流值不應(yīng)該超過其限值,則各線路潮流約束條件可寫為: (13)由于總的負(fù)荷需求在一時(shí)段內(nèi)是不會(huì)改變的,在調(diào)整時(shí)有的機(jī)組的出力將會(huì)增加,有的機(jī)組出力將會(huì)減少,其調(diào)整值為,可取正值也可取負(fù)值,所有發(fā)電機(jī)組調(diào)整值的波動(dòng)總和應(yīng)該為0,則有: (14)又加上各機(jī)組受爬坡速率的影響,所調(diào)整的出力值不能夠超過其爬坡的*圍,即在分配預(yù)案的調(diào)整應(yīng)該在爬坡速率的約束*圍之內(nèi)，故約束條件可寫為: (15)其中,所以模型的具體形式為:minst.問題五當(dāng)下一個(gè)時(shí)段預(yù)報(bào)的負(fù)荷需求經(jīng)過問題三,按照電力市場(chǎng)規(guī)則,給出下一個(gè)時(shí)段各機(jī)組的出力分配預(yù)案,再到問題四,檢查分配預(yù)案是否會(huì)引起輸電阻塞,如果不會(huì),則接受各機(jī)組出力分配方案;否則調(diào)整各機(jī)組出力分配方案來使得輸電阻塞消除。若無論怎樣調(diào)整都不能消除輸電阻塞，則執(zhí)行“輸電阻塞管理原則”的第二條原則，可以使用線路的安全裕度輸電，以避免拉閘限電，但必須保證“安全第一”且較經(jīng)濟(jì)的原則，再來重新調(diào)整各機(jī)組出力分配方案。由以上分析可知，當(dāng)無論怎樣調(diào)整都不能消除輸電阻塞時(shí)，為調(diào)整出安全且經(jīng)濟(jì)的各機(jī)組出力分配方案，則可以各條線路上最大的潮流的絕對(duì)值超過限值的百分比為最小為和購電費(fèi)最小為函數(shù)：為第條線路上最大的潮流的絕對(duì)值超過限值的百分比。除了問題四的全部約束外，同時(shí)要求各條線路上潮流的絕對(duì)值超過限值的百分比均不超過安全裕度其中為第條線路上的安全裕度，則調(diào)整分配方案的模型為： st.四、模型求解4.1問題一求解本問題中共8臺(tái)發(fā)電機(jī)，6條主要線路，即；則式(2)具體形式為： (16)對(duì)于每條線路中的潮流值都各自對(duì)應(yīng)一多元線性回歸模型，為了求得…的最小二乘估計(jì),令: (17)其中表示有32組數(shù)據(jù)實(shí)驗(yàn)值，表示第條線路上的第組潮流實(shí)驗(yàn)值，對(duì)式（4）求導(dǎo)，由得到正則方程組: (18)以上方程組可以近似的表示為： (19)表示各發(fā)電機(jī)組出力的實(shí)驗(yàn)值矩陣，為的轉(zhuǎn)值矩陣，即：由此可得到的最小二乘估計(jì)值為： (20)從而得到經(jīng)驗(yàn)回歸函數(shù)為: (21)將附錄2中表（一）表（二）的數(shù)據(jù)代入上式中，通過MATLAB編程求解可得如下結(jié)果：4.2問題三求解此模型為一單目標(biāo)線形規(guī)劃，可用常規(guī)的優(yōu)化算法來求解，但是，由于模型中各機(jī)組的每個(gè)段容量相應(yīng)的對(duì)應(yīng)一個(gè)段價(jià)，于是此模型的目標(biāo)函數(shù)為分段函數(shù)，用一般的優(yōu)化算法不能求解。所以我們根據(jù)市場(chǎng)交易規(guī)則，按段價(jià)從低到高選取各機(jī)組滿足爬坡的約束*圍內(nèi)的段容量或其部分，直到他們之和等于預(yù)報(bào)的負(fù)荷需求。每臺(tái)機(jī)組在當(dāng)前時(shí)段內(nèi)就要根據(jù)下一個(gè)時(shí)段的負(fù)荷預(yù)報(bào)，即負(fù)荷需求，每臺(tái)機(jī)組的段價(jià)、出力的爬坡速率，按從低到高選取各機(jī)組的段容量或者其部分，直到它們的和等于負(fù)荷需求。此時(shí)每個(gè)機(jī)組被選入的段容量或其部分之和就為該時(shí)段該機(jī)組的出力分配預(yù)案，即求得模型的最優(yōu)解。根據(jù)實(shí)際應(yīng)用，下一時(shí)段的負(fù)荷需求MW，當(dāng)，，即8臺(tái)機(jī)組6條線路。其段價(jià)、出力的爬坡速率和段容量的具體值分別見附錄2表(三)、表(四)、表(五)，通過數(shù)學(xué)軟件MATLAB編程（程序流程圖見附錄3）求解，可以得到:該時(shí)段的分配預(yù)案的清算價(jià)格為:=303元/MWh其最小的費(fèi)用為：=74416.8元各機(jī)組的分配預(yù)案如表(一)：表(一)負(fù)荷需求為982.4MW各機(jī)組的分配預(yù)案機(jī)組12345678出力值150.079.0180.099.5125.0140.095.0113.94.3問題四求解根據(jù)問題三所得出的分配預(yù)案，代入模型一所得的表達(dá)式中進(jìn)行檢驗(yàn)，通過計(jì)算各線路的潮流值見下表表（二）各線路的潮流值線路123456潮流值173.3084141.0167150.9190120.9034136.8083168.5149限值165150160155132162從上表可以看出，第一、第五、第六條線路上的潮流值均超過了其限值，超出值的百分比分別為：5.04%、3.64%、4.02%，說明分配預(yù)案會(huì)使得線路會(huì)產(chǎn)生阻塞，根據(jù)輸電阻塞管理原則，應(yīng)該調(diào)整各機(jī)組分配方案，即對(duì)模型四進(jìn)行求解。由于模型四是以阻塞費(fèi)用最小為目標(biāo)函數(shù)，根據(jù)我們的模型二的規(guī)則，模型四的目標(biāo)函數(shù)為一分段函數(shù)，約束條件都是線形的，故我們可以用軟件MATLAB 中優(yōu)化工具箱的函數(shù)fmincon進(jìn)行求解，其參數(shù)中的目標(biāo)函數(shù)為我們模型二關(guān)于阻塞費(fèi)用的計(jì)算規(guī)則，詳細(xì)的處理方法見附錄，以下是調(diào)整方案后各機(jī)組的出力情況及各線路的潮流值見表（三）和表（四）表（三）各機(jī)組出力方案機(jī)組12345678出力值152.597488.0228.081.7434152.098.065.0592117.0表（四）各線路的潮流值機(jī)組123456潮流值164.9582149.3564155.2325125.4252131.9959159.4773限值165150160155132162其清算價(jià)格取=303元/MWh阻塞費(fèi)用：12026.0元總費(fèi)用為：=86140.8元4.4問題五求解當(dāng)負(fù)荷需求為MW時(shí),按照解決問題三的方法進(jìn)行求解，得到的分配預(yù)案如表（五）:表(五)負(fù)荷需求為1052.8MW時(shí)各機(jī)組的分配預(yù)案機(jī)組12345678出力值152.081.0218.299.5135.0150.0102.1117.0清算價(jià)格為:=356元/MWh,總費(fèi)用為:=93699元.將根據(jù)線路上有功潮流關(guān)于各發(fā)電機(jī)組的出力近似表達(dá)式,得到各線路的潮流值如表（六）：表（六）負(fù)荷需求為1052.8MW時(shí)各線路的潮流值線路123456潮流值177.2415141.1811156.1459126.7333134.8112167.0558限值165150160155132162從上表可以看出第一、第五、第六條線路上的潮流值均超過了其限值，超出值的百分比分別為：7.42%、2.13%、3.12%。說明分配預(yù)案會(huì)產(chǎn)生輸電阻塞,根據(jù)輸電阻塞管理原則，需要調(diào)整各機(jī)組分配方案，調(diào)整后的方案為:表(七)負(fù)荷需求為1052.8MW時(shí)各機(jī)組的調(diào)整分配方案（MW）機(jī)組12345678出力值149.888.0228.299.5152.0128.090.3117.0阻塞費(fèi)用為：元總費(fèi)用為：=94924.8元五、結(jié)果分析與檢驗(yàn)5.1問題一結(jié)果分析由各機(jī)組的出力方案和與之相對(duì)應(yīng)的各線路的潮流值，我們通過多元線性回歸得到其各線路上對(duì)應(yīng)的有功潮流和發(fā)電機(jī)組出力的近似關(guān)系式為線性的。然而在回歸過程中，每一條線路的有功潮流都有一個(gè)方案與之明顯的不成線性關(guān)系，但是對(duì)整個(gè)關(guān)系式?jīng)]有多大的影響，于是我們便忽略了其對(duì)潮流和出力的關(guān)系的影響。我們用的回歸算法是基于LIPSOL,算法的主要原理:首先把問題化為線形規(guī)劃的標(biāo)準(zhǔn)形式,再引入松散變量,用牛頓迭代法,采用預(yù)報(bào)較正算法來搜索.我們利用MATLAB軟件回歸得到了六條線路的近似表達(dá)式.并給出了檢驗(yàn)了回歸模型(程序見附錄)表（八）第一條線路近似表達(dá)式的檢驗(yàn)結(jié)果如下:常數(shù)項(xiàng)機(jī)組1機(jī)組2機(jī)組3機(jī)組4機(jī)組5機(jī)組6機(jī)組7機(jī)組8回歸系數(shù):110.47750.08260.04780.05280.1199-0.02570.12160.1220-0.0015系數(shù)估計(jì)區(qū)間下限109.54210.08080.04370.05140.1166-0.02770.11900.1189-0.0037系數(shù)估計(jì)區(qū)間上限111.41290.08440.05180.05420.1231-0.02370.12430.12510.0007表（九）回歸各點(diǎn)的殘差方案012345678910殘差0.0615-0.01590.0040-0.00880.02010.04580.0356-0.03620.02670.0016-0.0158置信區(qū)間-0.0092-0.0591-0.0801-0.0386-0.0875-0.1185-0.1103-0.0061-0.0941-0.0485-0.05890.13210.09100.07200.05630.04720.02690.03920.07850.04070.04530.0905方案1112131415161718192021殘差-0.00830.0455-0.0080-0.00450.0485-0.0251-0.01230.0223-0.04090.00780.0570置信區(qū)間-0.0600-0.1192-0.0657-0.0534-0.1176-0.0236-0.0524-0.0950-0.0334-0.0629-0.12450.07650.02820.08170.06230.02070.07390.07710.05030.11520.04730.0105方案2223242526272829303132殘差-0.06600.0013-0.01660.0055-0.05820.0366-0.02110.0360-0.01370.0363-0.0177置信區(qū)間-0.0047-0.0698-0.0363-0.0602-0.0111-0.1000-0.0546-0.1078-0.0538-0.1108-0.03380.13680.06720.06950.04920.12750.02670.09670.03580.08110.03830.0691回歸相關(guān)系數(shù)，，對(duì)應(yīng)的概率為，相關(guān)系數(shù)越接近于值越大就說明回歸顯著，從上表中的數(shù)據(jù)也可以看出，回歸各點(diǎn)的殘差都是很小的?；貧w效果顯著，說明我們的模型具有很高的精度，能夠很好的反應(yīng)出各線路的潮流值與各機(jī)組出力的關(guān)系。但是有于在回歸方程中有一個(gè)常數(shù)項(xiàng)，當(dāng)各機(jī)組的出力為而線路中的有功潮流值不為，這在現(xiàn)實(shí)中就不能解釋了。我們初步考慮是這只是一個(gè)近似的表達(dá)式，在一定的*圍內(nèi)滿足就行了。5.2問題三結(jié)果分析根據(jù)表(1)所得結(jié)果，可以算出各機(jī)組的出力總和剛好滿足負(fù)荷需求982.4MW，并且各機(jī)組的段容量是在其爬坡*圍內(nèi)，見下表：表（十）各機(jī)組出力值（MW）機(jī)組12345678上限1538822899.5132155102.1117下限975813260.5989561.163出力150.079.0180.099.5125.0140.095.0113.9可以看出各機(jī)組的出力都在爬坡速率的限制之內(nèi)，符合預(yù)期的要求。在我們的出力分配方案中最后一個(gè)被選入的段價(jià)是第八機(jī)組中的第七個(gè)段容量，清算價(jià)為元，仔細(xì)分析這個(gè)結(jié)果，可以看出其他被選入的段容量所對(duì)應(yīng)的段價(jià)都小于該清算價(jià)。結(jié)合題中的數(shù)據(jù)我們還發(fā)現(xiàn)只有最后一個(gè)被選入的可以再繼續(xù)增大，比如我們這里他還能增大，其它的都不能再增大了。可以看出的這種現(xiàn)象是符合我們選入的規(guī)則，這就說明我們所建立的模型是合理、可靠的；5.3問題四結(jié)果分析我們把原方案與調(diào)整方案進(jìn)行了對(duì)比，見表（十一）表（十一）原預(yù)案與調(diào)整方案的比較機(jī)組12345678預(yù)案150.079.0180