如何解析電力市場(chǎng)的輸電阻塞管理

PAGEPAGE21電力市場(chǎng)的輸電阻塞管理摘要本模型根據(jù)下一時(shí)段的負(fù)荷預(yù)報(bào)、每臺(tái)機(jī)組的報(bào)價(jià)、當(dāng)前出力以及交易規(guī)則建立了一個(gè)分配預(yù)案模型。當(dāng)執(zhí)行分配預(yù)案時(shí)，如果有線路出現(xiàn)輸電阻塞，遵循電網(wǎng)“安全第一”的原則和輸電阻塞管理原則建立了輸電阻塞管理的優(yōu)化模型。就題目所給的五個(gè)問(wèn)題本文建立了四個(gè)模型。模型1通過(guò)一個(gè)物理學(xué)的公式，用差值的計(jì)算方法求出了有功潮流關(guān)于各機(jī)組出力的近似表達(dá)式，巧妙的避開(kāi)了因改變量相對(duì)較小，誤差不可忽略時(shí)而使擬和精度很低的情況，并通過(guò)一個(gè)數(shù)學(xué)公式給出了阻塞費(fèi)用的計(jì)算規(guī)則。在模型3中引入了一個(gè)定理,使得只要計(jì)算較少的負(fù)荷差的優(yōu)化選取即可求得最優(yōu)方案，同時(shí)根據(jù)爬坡速率的限制給出了可以選取的范圍，大大減少了運(yùn)算量。對(duì)于模型4建立了一個(gè)優(yōu)化模型，求解阻塞費(fèi)用時(shí)，根據(jù)安全且經(jīng)濟(jì)的原則，建立了以阻塞費(fèi)用、電網(wǎng)系統(tǒng)的穩(wěn)定性和安全性為目標(biāo)的多目標(biāo)函數(shù)，但通過(guò)將電網(wǎng)系統(tǒng)的穩(wěn)定性和安全性的約束轉(zhuǎn)換為“損失的錢(qián)數(shù)”的約束，將問(wèn)題最終轉(zhuǎn)化為求各項(xiàng)“錢(qián)數(shù)”之和最少的單目標(biāo)優(yōu)化問(wèn)題，有效的降低了問(wèn)題求解的難度，并且模型具有很高的推廣性。針對(duì)982.4的負(fù)荷需求，通過(guò)調(diào)整各機(jī)組的出力分配方案可以使輸電阻塞消除，通過(guò)matlab編程求出最小的阻塞費(fèi)用為3957.5元，但對(duì)于1052.8的負(fù)荷需求，無(wú)論怎樣調(diào)整各機(jī)組的出力分配方案都不能消除輸電阻塞，此時(shí)要使用線路的安全裕度輸電，同樣用matlab編程求出“錢(qián)數(shù)”之和最小為4643.5元。一、問(wèn)題重述電力從生產(chǎn)到使用的四大環(huán)節(jié)——發(fā)電、輸電、配電和用電是瞬間完成的。我國(guó)電力市場(chǎng)初期是發(fā)電側(cè)電力市場(chǎng)，采取交易與調(diào)度一體化的模式。電網(wǎng)公司在組織交易、調(diào)度和配送時(shí)，必須遵循電網(wǎng)“安全第一”的原則，同時(shí)要制訂一個(gè)電力市場(chǎng)交易規(guī)則，按照購(gòu)電費(fèi)用最小的經(jīng)濟(jì)目標(biāo)來(lái)運(yùn)作。設(shè)某電網(wǎng)有若干臺(tái)發(fā)電機(jī)組和若干條主要線路，每條線路上的有功潮流（輸電功率和方向）取決于電網(wǎng)結(jié)構(gòu)和各發(fā)電機(jī)組的出力。電網(wǎng)每條線路上的有功潮流的絕對(duì)值有一安全限值，限值還具有一定的相對(duì)安全裕度。如果各機(jī)組出力分配方案使某條線路上的有功潮流的絕對(duì)值超出限值，稱為輸電阻塞。當(dāng)發(fā)生輸電阻塞時(shí)，需要研究如何制訂既安全又經(jīng)濟(jì)的調(diào)度計(jì)劃。我們要做的工作有：（電力市場(chǎng)交易規(guī)則和輸電阻塞管理原則見(jiàn)原題。）某電網(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ì)一種簡(jiǎn)明、合理的阻塞費(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.設(shè)各機(jī)組在運(yùn)行期間能正常運(yùn)行，不會(huì)出現(xiàn)因故障等情況而停機(jī)且模型只考慮網(wǎng)方的費(fèi)用支出，不考慮電力廠的利潤(rùn)。2.設(shè)在一個(gè)時(shí)段內(nèi)用電量不會(huì)相差太大，設(shè)用電量為一定值。3.設(shè)對(duì)于任一時(shí)段，機(jī)組的出力只要在該時(shí)段的最后時(shí)刻滿足方案即認(rèn)為滿足要求，而不考慮中間變化是否滿足。4.設(shè)機(jī)組的報(bào)價(jià)正常，沒(méi)有因企業(yè)追逐利潤(rùn)而出現(xiàn)的策略報(bào)價(jià)等情況。5.設(shè)調(diào)度中心當(dāng)前時(shí)段采用的調(diào)度方案是最優(yōu)的。三、符號(hào)說(shuō)明 下一時(shí)段預(yù)報(bào)的總負(fù)荷第臺(tái)機(jī)組在當(dāng)前時(shí)段的實(shí)際出力第臺(tái)機(jī)組在下一時(shí)段的預(yù)測(cè)出力第臺(tái)機(jī)組在下一時(shí)段的實(shí)際出力第條線路的預(yù)測(cè)潮流值第條線路的實(shí)際潮流值第條線路上有功潮流絕對(duì)值的安全限值第條線路上安全限值的相對(duì)安全裕度第臺(tái)機(jī)組的爬坡速率第條線路的有功潮流的絕對(duì)值超過(guò)限值的百分比四、模型分析本問(wèn)題是一個(gè)電力市場(chǎng)的輸電阻塞管理問(wèn)題，即在電力市場(chǎng)輸電阻塞管理中，電力交易調(diào)度中心如何給出下個(gè)時(shí)段的各機(jī)組的出力分配預(yù)案，以及發(fā)生輸電阻塞時(shí)如何實(shí)施阻塞管理。由于下個(gè)時(shí)段各機(jī)組的出力分配預(yù)案和下個(gè)時(shí)段的負(fù)荷預(yù)報(bào)、每臺(tái)機(jī)組的報(bào)價(jià)、當(dāng)前出力、以及出力改變速率有關(guān)。我們可以跟據(jù)下個(gè)時(shí)段的負(fù)荷預(yù)報(bào)、每臺(tái)機(jī)組的報(bào)價(jià)、當(dāng)前出力、可用出力以及出力改變速率，依據(jù)電力市場(chǎng)交易規(guī)則構(gòu)造一個(gè)數(shù)學(xué)模型來(lái)給出下個(gè)時(shí)段各個(gè)機(jī)組的出力分配預(yù)案。發(fā)生輸電阻塞時(shí)就需要遵循“安全第一”的電網(wǎng)原則以及輸電阻力管理規(guī)則來(lái)調(diào)整原來(lái)的各機(jī)組出力分配預(yù)案。而在調(diào)整根據(jù)電力市場(chǎng)交易規(guī)則得到的各機(jī)組出力分配預(yù)案時(shí)，就會(huì)造成網(wǎng)方和發(fā)電方的利益沖突。一方面是由于序內(nèi)容量不能出力而對(duì)發(fā)電方造成的損失。網(wǎng)方要給予一定的經(jīng)濟(jì)賠償。另一方面由于段價(jià)高于清算價(jià)而對(duì)發(fā)電方造成的損失。網(wǎng)方也要給予經(jīng)濟(jì)賠償。這兩方面費(fèi)用稱為阻塞費(fèi)用。因此在調(diào)整電力分配預(yù)案時(shí)就要考慮在安全運(yùn)行的前提下盡量減少阻塞費(fèi)用。我們可以以阻塞費(fèi)用最小為目標(biāo)，在考慮安全第一、有功潮流和爬坡速率等約束下。建一個(gè)優(yōu)化模型來(lái)求最佳的阻塞調(diào)度方案。五、模型建立及求解模型1為了便于處理和分析，我們首先引入以下的定理定理：在電網(wǎng)結(jié)構(gòu)一定的情況下，即一定的情況下，如果不變，則。并稱為變換因子。證明（見(jiàn)附錄4）：對(duì)于本問(wèn)題題目給出了各機(jī)組出力和各線路上對(duì)應(yīng)的有功潮流值的一些實(shí)驗(yàn)數(shù)據(jù)，要求用這些實(shí)驗(yàn)數(shù)據(jù)確定各線路上有功潮流關(guān)于各發(fā)電機(jī)組出力的近似表達(dá)式。根據(jù)引入的定理很容易就能得出如下的關(guān)系式：式中為機(jī)組對(duì)線路的變換因子。當(dāng)電機(jī)組的出力有變化，而其他機(jī)組的出力不變時(shí)，線路有變化由表達(dá)式解得：對(duì)于組不同的變化，從而求出個(gè)不同的變換因子。加和求平均就可得出機(jī)組對(duì)線路的變換因子：，對(duì)于本題，由題中表1可知，當(dāng)?shù)谂_(tái)機(jī)組出力變化而其它機(jī)組的出力不變都給出了四組數(shù)據(jù)，所以取。下面給出第一臺(tái)機(jī)組對(duì)第一條線路的權(quán)重的計(jì)算過(guò)程。取0方案為基準(zhǔn)數(shù)據(jù)，具體計(jì)算如下：；；；；可求得：判斷有無(wú)異常數(shù)據(jù)：根據(jù)物理中的拉依達(dá)準(zhǔn)則[2]：由于，所以沒(méi)有異常數(shù)據(jù)需要剔除。其它的具體計(jì)算與此相同不在重復(fù)。（實(shí)際計(jì)算我們是編程用計(jì)算機(jī)處理的，具體程序見(jiàn)附錄1）。具體的計(jì)算結(jié)果如下：表1變換因子0.0789-0.0583-0.0681-0.0373-0.00470.23600.03070.11810.0698-0.10950.2352-0.07290.0489-0.0018-0.15340.2020-0.0661-0.08140.11980.0332-0.0118-0.0202-0.04120.0961-0.02870.08320.1271-0.0135-0.07000.04340.1183-0.11540.00560.00210.0647-0.00450.1262-0.0146-0.00740.14930.00170.1738-0.00710.0937-0.19960.0743-0.0133-0.0043由此得到有潮功率關(guān)于機(jī)組的出力的近似表達(dá)式：模型2由題目可知，阻塞管理的關(guān)鍵是給出一個(gè)合理、簡(jiǎn)明的阻塞費(fèi)用計(jì)算規(guī)則。而電網(wǎng)要付的阻塞費(fèi)用可分成兩部分：一部分是由于電網(wǎng)阻塞而調(diào)整各機(jī)組出力，使得部分機(jī)組的出力少于預(yù)算分配方案，從而對(duì)機(jī)組造成一定的經(jīng)濟(jì)損失。網(wǎng)方要給予一定的經(jīng)濟(jì)賠償。也就是對(duì)序內(nèi)容量不出力部分的賠償；記為；另一部分是由于調(diào)整出力方案，使的個(gè)別機(jī)組出力比預(yù)算方案多，當(dāng)額外增加的段容量所對(duì)應(yīng)的段價(jià)高于清算價(jià)時(shí)，由于統(tǒng)一按清算價(jià)處理而對(duì)機(jī)組造成一定的經(jīng)濟(jì)損失。網(wǎng)方要給予一定的補(bǔ)償。這就是對(duì)序外容量出力部分的賠償。記為；要解決的問(wèn)題就是給出、的計(jì)算規(guī)則。只有當(dāng)調(diào)整后的機(jī)組出力低于分配預(yù)案出力時(shí)，電網(wǎng)公司才會(huì)支付序內(nèi)容量不能出力部分的補(bǔ)償，否則不予補(bǔ)償。因此定義一個(gè)0—1變量來(lái)描述是否需要給予補(bǔ)償。其中表示第臺(tái)機(jī)組的預(yù)測(cè)出力，表示第臺(tái)機(jī)組的實(shí)際出力。由上面的分析可知對(duì)序內(nèi)容量不出力部分的賠償可用公式表示為：其中是對(duì)序內(nèi)容量不能出力部分的補(bǔ)償系數(shù)，由于序內(nèi)容量不能出力的部分仍由電機(jī)廠所保留，所以不能完全按清算價(jià)給予補(bǔ)償，我們?nèi)?duì)序外容量出力部分的賠償可用公式表示為：表示只有當(dāng)電力公司的實(shí)際出力大于分配預(yù)案出力且超出部分的段價(jià)高于清算價(jià)時(shí)電力公司才支付序外容量出力部分賠償；是對(duì)序外容量的補(bǔ)償系數(shù)。由于對(duì)于一個(gè)序外容量，其未被入選的原因就是該段報(bào)價(jià)高，所以在阻塞調(diào)整時(shí)，對(duì)其進(jìn)行補(bǔ)償時(shí)不能取段價(jià)與清算價(jià)的差。因?yàn)檫@樣就可能會(huì)導(dǎo)致電力公司采用策略性報(bào)價(jià)的情況。所以我們應(yīng)使其低于報(bào)價(jià)與清算價(jià)的差用公式可表為：式中表示序外容量在高于清算價(jià)出力部分的報(bào)價(jià)。由上可知阻塞費(fèi)用的計(jì)算規(guī)則為：由此公式，只要知道相關(guān)信息即可方便的求出阻塞費(fèi)用。模型3本模型考慮的網(wǎng)方的購(gòu)電費(fèi)用，即以最小為目標(biāo)給出一個(gè)出力分配預(yù)案。由電力市場(chǎng)交易規(guī)則，對(duì)某一時(shí)段的全部機(jī)組出力都按清算價(jià)結(jié)算。所以有：式中代表清算價(jià)。對(duì)某一時(shí)段來(lái)講為一定值，其大小就等于該時(shí)段內(nèi)所選中段容量中最高的段價(jià)?？杀硎緸椋浩渲斜硎緳C(jī)組在段的報(bào)價(jià)表示被選入段的段價(jià)。由于負(fù)荷的要求對(duì)于機(jī)組的出力應(yīng)滿足：對(duì)于某一固定時(shí)段，清算價(jià)為定值。所以有：這樣以為最小目標(biāo)的求分配預(yù)案可轉(zhuǎn)化為以為最小目標(biāo)求分配預(yù)案。同時(shí)由于受到爬坡速率的限制，由假設(shè)（3）對(duì)于任一時(shí)段，機(jī)組的出力只要在該時(shí)段的最后時(shí)刻滿足方案即認(rèn)為滿足要求，而不考慮中間變化是否滿足。可知要使機(jī)組的出力滿足方案要求，需滿足：綜上可得問(wèn)題的優(yōu)化模型：這是一個(gè)非線性規(guī)劃的離散模型，可由計(jì)算機(jī)編程求解（見(jiàn)附錄2）。具體的求解流程如下：由市場(chǎng)交易規(guī)則可知，市場(chǎng)交易調(diào)度中心在調(diào)配各機(jī)組出力時(shí)，按段價(jià)從低到高選取各機(jī)組的段容量或其部分，直到它們之和等于預(yù)報(bào)的負(fù)荷?？紤]到下個(gè)時(shí)段的方案是以當(dāng)前時(shí)段的方案為基礎(chǔ)變化而來(lái)，同時(shí)還要受到式的約束而使機(jī)組的出力變動(dòng)在一定范圍之內(nèi)。機(jī)組出力的變動(dòng)范圍為：由于選取段容量時(shí)按段價(jià)由低到高進(jìn)行選取?？蓪⒏鞫味蝺r(jià)進(jìn)行一個(gè)變換，將變動(dòng)范圍之外段容量所對(duì)應(yīng)的段價(jià)取為無(wú)窮，即讓該段不可能被取到。為了求解的方便，我們引入以下定理。定理：如果當(dāng)前時(shí)段的方案是最優(yōu)的，且下個(gè)時(shí)段的負(fù)荷大于當(dāng)前時(shí)段。那么當(dāng)前時(shí)段方案所包含的部分一定包含在下個(gè)時(shí)段的最優(yōu)方案中。證明（反證法）：假設(shè)下個(gè)時(shí)段的最優(yōu)方案不包括當(dāng)前時(shí)段的方案，則下個(gè)時(shí)段在當(dāng)前時(shí)段的負(fù)荷點(diǎn)，也會(huì)有一個(gè)最優(yōu)的方案。而當(dāng)前時(shí)段也是最優(yōu)的。這時(shí)對(duì)于同一個(gè)負(fù)荷就有兩個(gè)最優(yōu)方案，而這于事實(shí)不符，假設(shè)不成立。即當(dāng)前時(shí)段所包含的部分一定包含在下個(gè)時(shí)段的最優(yōu)方案中。由假設(shè)（5）和定理可知，對(duì)于本問(wèn)題當(dāng)前時(shí)段的方案一定包含于下時(shí)段的最優(yōu)方案中。所以我們只考慮下時(shí)段的負(fù)荷高于當(dāng)前時(shí)段的部分即只要使在余下的段容量中獲得最優(yōu)分配，就能使所有的分配在當(dāng)前時(shí)段獲得最優(yōu)分配。同上也可以將已被選取的段容量的段價(jià)置為無(wú)窮。即以后不再考慮這些段的選取。同時(shí)將段容量為0的段價(jià)置為無(wú)窮。即不考慮這些段的選取，因?yàn)槠涠稳萘繛?，即使選取也毫無(wú)意義。這樣經(jīng)過(guò)變換的段價(jià)表如表2；表2：各機(jī)組變換后的段價(jià)（單位：元/兆瓦小時(shí)，記作元/MWh）機(jī)組\段123456789101InfInfInfInfInf252InfInfInf4892InfInfInfInfInf300320InfInf4953InfInfInfInfInfInf308356InfInf4InfInfInfInf255302InfInfInfinf5InfInfInfInfInfInfInf3103965106InfInfInfInfInf252305380InfInf7InfInfInf251260306InfInfInfInf8InfInfInfInf253Inf303InfInfInf同樣也可將段容量表進(jìn)行變換，使各機(jī)組的段容量轉(zhuǎn)換為還可被選取的段容量。即將已選取完的段容量和選取范圍之外的段容量置為0，而部分能被選取的段容量則只保留能被選取的部分，變換后可得表2表3：各機(jī)組變換后的段容量（單位：MW）機(jī)組\段1234567891010000030000320000062007300000020280040000109.500005000000010107600000151050070003.9107.10000800002007000轉(zhuǎn)換后再計(jì)算時(shí)只需從段價(jià)表中找出最小的段價(jià)，取出它所對(duì)應(yīng)的段容量，看能否滿足負(fù)荷增量，如果不能滿足則取出其段容量或其部分，令其段價(jià)為無(wú)窮，再重復(fù)上述過(guò)程直到取出的段容量之和滿足為之。此時(shí)求得機(jī)組出力的變動(dòng)量加上其當(dāng)前的出力就可求的下一時(shí)段的出力值。綜上求得機(jī)組出力的改變?yōu)椋?30，6，0，19.5，0，15，13.9，23.9)；機(jī)組的分配預(yù)案為：(150，79，180，99.5，125，140，95，113.9)；此時(shí)清算價(jià)q最小且值為：=303網(wǎng)方的最小購(gòu)電費(fèi)用：=297670模型4由題意可知，是否會(huì)引起輸電阻塞就是線路上的有功潮流的絕對(duì)值是否超過(guò)限值。即只要有成立，就會(huì)發(fā)生輸電阻塞。發(fā)生輸電阻塞時(shí)，由模型分析可知，可分做兩種情形來(lái)處理。情形在使各線路的有功潮流值不超過(guò)限值的前提下進(jìn)行調(diào)整，盡量消除阻塞。情形在安全裕度范圍內(nèi)，在阻塞費(fèi)用最小和每條線路上潮流的絕對(duì)值超過(guò)限值的百分比盡量小的前提下進(jìn)行優(yōu)化。對(duì)于阻塞費(fèi)用，進(jìn)行調(diào)整時(shí)應(yīng)使阻塞費(fèi)用盡可能的小。由模型2可知同時(shí)考慮阻塞調(diào)整方案要求線路上潮流的絕對(duì)值超過(guò)限值的百分比盡量小。可考慮目標(biāo)函數(shù)：考慮到爬坡速率的約束，根據(jù)假設(shè)可知，要使機(jī)組出力在時(shí)段的最后時(shí)刻滿足條件。有約束條件：考慮到電網(wǎng)安全第一的原則，線路上的有功潮流的絕對(duì)值超過(guò)限值的百分比不能超過(guò)相對(duì)安全裕度，有約束條件：根據(jù)模型1所求出的線路上的有功潮流和機(jī)組出力的關(guān)系，有約束條件：線路上的有功潮流的和要滿足負(fù)荷要求，有約束條件對(duì)于情形由于要求線路上的有功潮流小于各線路潮流限值。只需考慮阻塞費(fèi)用最小和系統(tǒng)的穩(wěn)定性即可。綜上可得情形的目標(biāo)函數(shù)為：當(dāng)情形不滿足時(shí)，我們考慮情形，即允許超過(guò)線路潮流限值，但要在線路的安全裕度內(nèi)輸電。而超過(guò)線路的潮流限值進(jìn)行輸電必然會(huì)造成額外的損失。當(dāng)超過(guò)的比例比較小時(shí)，隨著超過(guò)比例的增加，損失會(huì)增加很快。但當(dāng)比例達(dá)到一定程度時(shí)，損失將增長(zhǎng)緩慢最后將趨于一定值。引入Michaelis-Menten模型[3]表示條線路上因超過(guò)限值而造成的損失。對(duì)于易知當(dāng)取最大時(shí)，將超過(guò)最大的安全裕度，會(huì)帶來(lái)電網(wǎng)的損壞，損失將趨于無(wú)窮，所以我們?nèi)《鴮?duì)于，為使所算的與實(shí)際相符，我們?nèi)?，用?jì)算機(jī)進(jìn)行處理，得出合理值為=20000。用代替優(yōu)化目標(biāo)綜上可得情形的優(yōu)化模型：解這個(gè)單目標(biāo)非線性優(yōu)化模型即可給出調(diào)度后的各機(jī)組的出力以及相應(yīng)的阻塞費(fèi)用。模型4的求解根據(jù)模型(1)求得的線路上有功潮流關(guān)于機(jī)組出力的近似表達(dá)式以及機(jī)組預(yù)測(cè)出力，通過(guò)計(jì)算機(jī)的求解（具體程序見(jiàn)附錄3）。求得各線路上的潮流的值如下表：表4：根據(jù)預(yù)測(cè)方案求出的線路上的有功潮流值123456173.0263140.7825-150.8936120.8026136.5831168.4521表5：各線路的潮流限值（單位：MW）和相對(duì)安全裕度線路123456限值165150160155132162安全裕度13%18%9%11%15%14%最大安全裕度時(shí)的值186．45177174．4172．05151．8184．68比較線路上的有功潮流值和潮流限值可知線路1和線路5和線路6上發(fā)生輸電阻塞。根據(jù)模型分兩種情形進(jìn)行阻塞調(diào)整根據(jù)情形(1)進(jìn)行阻塞調(diào)整，盡量消除阻塞。這是一個(gè)多目標(biāo)的非線性離散模型，用計(jì)算機(jī)編程進(jìn)行求解（具體程序見(jiàn)附錄3）具體的求解流程如下調(diào)整后的方案如下表：表612345678149.98780.976225.84073.485152.000123.01360.100117.000.調(diào)整后各線路上的有功潮流的值如下表：表7123456165.0000145.0772-133.9965124.9914132.0000157.2429阻塞費(fèi)用的值為：=3957.5由計(jì)算結(jié)果可知，經(jīng)過(guò)調(diào)整可以消除阻塞。對(duì)于問(wèn)題5按模型(3)計(jì)算可得分配預(yù)案如下表表81234567815081218.299.5135150102.1117線路上的有功潮流計(jì)算如下表表9123456176.726140.815-155.958129.476134.446166.807和線路的潮流限值比較可知，線路1、線路5、線路6上產(chǎn)生輸電阻塞。按模型(4)進(jìn)行計(jì)算，對(duì)于情形(1)，經(jīng)計(jì)算發(fā)現(xiàn)無(wú)論怎樣調(diào)整出力方案都無(wú)法消除阻塞。在情形(2)的情況下進(jìn)行求解。計(jì)算結(jié)果如下調(diào)整后方案如下表表1012345678150.00082.175228.00099.500152.000140.00084.125117.000調(diào)整后各線路上的有功潮流如下表表11123456173.3016143.7667-133.3582128.3932132.2072163.582阻塞費(fèi)用的值為：=4643.5七、模型評(píng)價(jià)優(yōu)點(diǎn)：1求解問(wèn)題1的近似表達(dá)式時(shí)，我們引入了一個(gè)模型中證明的定理，有效地消除了傳輸過(guò)程中的電網(wǎng)損耗對(duì)計(jì)算結(jié)果的影響，提高了計(jì)算結(jié)果的精度，這就提高了以后模型最終結(jié)果的準(zhǔn)確性。有很強(qiáng)的操作性，在給出了一組初始出力初值，輸入程序即可得到輸電阻塞的解決方案。有很強(qiáng)的推廣性，可推廣到機(jī)組數(shù)任意多和傳輸線路任意多的情形。4求解阻塞費(fèi)用時(shí)，根據(jù)安全且經(jīng)濟(jì)的原則，建立了以阻塞費(fèi)用、電網(wǎng)系統(tǒng)的穩(wěn)定性和安全性為目標(biāo)的多目標(biāo)函數(shù)，但通過(guò)將電網(wǎng)系統(tǒng)的穩(wěn)定性和安全性的約束轉(zhuǎn)換為“損失的錢(qián)數(shù)”的約束，將問(wèn)題最終轉(zhuǎn)化為求各項(xiàng)“錢(qián)數(shù)”之和最少的單目標(biāo)優(yōu)化問(wèn)題，有效的降低了問(wèn)題求解的難度。缺點(diǎn)：模型的結(jié)果是計(jì)算機(jī)迭代得出的，與真實(shí)解有一定的差距。參考文獻(xiàn)：[1]蘇金明,張蓮花,劉波.MATLAB工具箱應(yīng)用[M].北京:電子工業(yè)出版社,2004[2]黃建偉.普通物理實(shí)驗(yàn)[M].哈爾濱:黑龍江教育出版社,2002[3]姜啟源,謝金星,葉俊.數(shù)學(xué)模型[M].北京:高等教育出版社,2004[4]張永平等:電力市場(chǎng)管理綜述http:///，04/9/17附錄一%此程序用于計(jì)算模型(1)中的機(jī)組1到各線路的變換因子%機(jī)組2--8到各線路的變換因子可用同樣的算法求出jizu=[120133.02129.63158.77145.32];xianlu=[164.78 140.87 -144.25 119.09 135.44 157.69;...165.81 140.13 -145.14 118.63 135.37 160.76;...165.51 140.25 -144.92 118.7 135.33 159.98;...167.93 138.71 -146.91 117.72 135.41 166.81;...166.79 139.45 -145.92 118.13 135.41 163.64];p=[];fori=2:size(xianlu,1)p(i-1,:)=(xianlu(i,:)-xianlu(1,:))/((jizu(i)-jizu(1)));endfori=1:6tt(i)=sum(p(:,i))/4;endtt附錄二%此算法用于計(jì)算模型(3)中的各機(jī)組的初始分配方案（初始交易值）drl=[70 0 50 0 0 30 0 0 0 40;...30 0 20 8 15 6 2 0 0 8;...110 0 40 0 30 0 20 40 0 40;...55 5 10 10 10 10 15 0 0 1;...75 5 15 0 15 15 0 10 10 10;...95 0 10 20 0 15 10 20 0 10;...50 15 5 15 10 10 5 10 3 2;...70 0 20 0 20 0 20 10 15 5];dj=[-505 0 124 168 210 252 312 330 363 489;...-560 0 182 203 245 300 320 360 410 495;...-610 0 152 189 233 258 308 356 415 500;...-500 150 170 200 255 302 325 380 435 800;...-590 0 116 146 188 215 250 310 396 510;...-607 0 159 173 205 252 305 380 405 520;...-500 120 180 251 260 306 315 335 348 548;...-800 153 183 233 253 283 303 318 400 800];cl0=[120 73 180 80 125 125 81.1 90];fhl=sum(cl0);form=1:size(drl,1)forn=1:size(drl,2)cl(m,n)=sum(drl(m,1:n));endendforu=1:size(cl,1)forv=1:size(cl,2)ifcl(u,v)>=cl0(u)%DJ(u)=dj(u,v);dj(u,1:v-1)=inf*ones(1,v-1);drl(u,v)=cl(u,v)-cl0(u);breakendendendv=[2.2 1 3.2 1.3 1.8 2 1.4 1.8];club=cl0+15*v;foru=1:size(cl,1)forv=1:size(cl,2)ifcl(u,v)>=club(u)%DJ(u)=dj(u,v);dj(u,v+1:10)=inf*ones(1,10-v);drl(u,v+1:10)=zeros(1,10-v);drl(u,v)=drl(u,v)-cl(u,v)+club(u)cl(u,v)=club(u);breakendendendfh0=1052.8-fhl;fh=0;djwhilefh<fh0[m,n,min]=min2(dj);dj(m,n)=inf;if(drl(m,n)~=0)fh=fh+drl(m,n);fdj(m,:)=[m,n,min];cl0(m)=cl(m,n);endendiffh>fh0fdrl=drl(m,n)-(fh-fh0);endfdjfdrlcl0%算法中用到的函數(shù)min2的代碼如下：function[m,n,min]=min2(x)m=1;n=1;min=x(1,1);fora=1:size(x,1)forb=1:size(x,2)ifx(a,b)<minmin=x(a,b);m=a;n=b;endendend附錄三%此算法用于求模型(4)中的阻塞費(fèi)用B=[165 150 160 155 132 162];%線路的潮流限值s=[0.13 0.180.09 0.11 0.15 0.14];%安全裕度v=[2.2 1 3.2 1.3 1.8 2 1.4 1.8];%爬坡速率QXJ=315;%清算價(jià)%x0=[1507918099.512514095113.9];%初始出力1%x0=[15081218.299.5135150102.1117];%初始出力2x0=[120 73 180 80 125 125 81.1 90];x0=x0';x1=[120 73 180 80 125 125 81.1 90];%基準(zhǔn)出力x1=x1';p1=[164.78 140.87 144.25 119.09 135.44 157.69]';%基準(zhǔn)線功率%Aeq=[ones(1,8),zeros(1,16)];%beq=0;Aeq=[];beq=[];bh=[0.0789-0.0583-0.0681-0.0373-0.00470.2360;...0.03070.11810.0698-0.10950.2352-0.0729;...0.0489-0.0018-0.15340.2020-0.0661-0.0814;...0.11980.0332-0.0118-0.0202-0.04120.0961;...-0.02870.08320.1271-0.0135-0.07000.0434;...0.1183-0.11540.00560.00210.0647-0.0045;...0.1262-0.0146-0.00740.14930.00170.1738;...-0.00710.0937-0.19960.0743-0.0133-0.0043];bh=bh';p0=bh*(x0-x1)+p1;p0=p0';X0=zeros(1,24);%[zeros(1,8),zeros(1,16)];%初始迭代值A(chǔ)=[bh,zeros(6,16)];b=B.*(1+s)-p1';xL=[-15*v,zeros(1,16)];xU=[15*v,15*v,15*v];%ones(1,16)*inf];globalnumber;number=0;options=optimset('MaxFunEvals',inf,'MaxIter',1000);%100[y,fval]=fmincon('zsfy',X0,A,b,Aeq,beq,xL,xU,'',options)numberp0pw=bh*y(1:8)'+p1%目標(biāo)函數(shù)zsfy的代碼如下：functionG=zsfy(x)x0=[120 73 180 80 125 125 81.1 90];%初始出力p1=[164.78 140.87 144.25 119.09 135.44 157.69]';bh=[0.0789-0.0583-0.0681-0.0373-0.00470.2360;...0.03070.11810.0698-0.10950.2352-0.0729;...0.0489-0.0018-0.15340.2020-0.0661-0.0814;...0.11980.0332-0.0118-0.0202-0.04120.0961;...-0.02870.08320.1271-0.0135-0.07000.0434;...0.1183-0.11540.00560.00210.0647-0.0045;...0.1262-0.0146-0.00740.14930.00170.1738;...-0.00710.0937-0.19960.0743-0.0133-0.0043];bh=bh';p0=bh*x(1:8)'+p1;p0=p0';globalnumber;number=number+1;B=[165 150 160 155 132 162];%線路的潮流限值drl=[70 0 50 0 0 30 0 0 0 40;...30 0 20 8 15 6 2 0 0 8;...110 0 40 0 30 0 20 40 0 40;...55 5 10 10 10 10 15 0 0 1;...75 5 15 0 15 15 0 10 10 10;...95 0 10 20 0 15 10 20 0 10;...50 15 5 15 10 10 5 10 3 2;...70 0 20 0 20 0 20 10 15 5];form=1:size(drl,1)forn=1:size(drl,2)cl(m,n)=sum(drl(m,1:n));endenddj=[-505 0 124 168 210 252 312 330 363 489;...-560 0 182 203 245 300 320 360 410 495;...-610 0 152 189 233 258 308 356 415 500;...-500 150 170 200 255 302 325 380 435 800;...-590 0 116 146 188 215 250 310 396 510;...-607 0 159 173 205 252 305 380 405 520;...-500 120 180 251 260 306 315 335 348 548;...-800 153 183 233 253 283 303 318 400 800];QXJ=303;DJ=zeros(1,8);chx=zeros(1,6);end%x0=[1507918099.512514095113.9];%初始出力1%x0=[15081218.299.5135150102.1117];%初始出力2x00=[1507918099.512514095113.9];fork=1:8ifx(k)+x0(k)-x00(k)>0x(8+k)=0;elsex(8+k)=-(x(k)+x0(k)-x00(k));endifx(k)+x0(k)-x00(k)<=0x(k+16)=0;elseifcl2dj(cl,dj,k,x(k)+x0(k))>QXJx(k+16)=x(k)+x0(k)-x00(k);DJ(k)=cl2dj(cl,dj,k,x(k)+x0(k));