電力市場(chǎng)的輸電阻塞模型 帶matlab程序

1、電力市場(chǎng)的輸電阻塞模型 摘要: 本文針對(duì)第一問(wèn)，建立了多元線性均值回歸模型，用于描述各線路上的有功潮流與各發(fā)電機(jī)組出力的關(guān)系。模型中利用均值定理，運(yùn)用matlab軟件求得解，此模型的誤差較小，較符合實(shí)際情況。 對(duì)于第二問(wèn)，我們引入基本補(bǔ)償價(jià)概念，利用限上和限下的方法設(shè)計(jì)阻塞費(fèi)用計(jì)算規(guī)則，在阻塞費(fèi)用計(jì)算規(guī)則里面我們認(rèn)為序外容量在低于對(duì)應(yīng)報(bào)價(jià)的清算價(jià)上出力的時(shí)候，網(wǎng)方已經(jīng)按清算價(jià)給予了發(fā)電商一部分的費(fèi)用，這部分費(fèi)用不作為經(jīng)濟(jì)補(bǔ)償，即其是不包含在阻塞費(fèi)用中，對(duì)于序外容量的經(jīng)濟(jì)補(bǔ)償我們就用對(duì)應(yīng)段容量的報(bào)價(jià)和清算價(jià)之差作為補(bǔ)償價(jià)格來(lái)計(jì)算阻塞費(fèi)用。 對(duì)于第三問(wèn)，在考慮爬坡速率和負(fù)荷需求的前提下，根據(jù)電力市

2、場(chǎng)規(guī)則，考慮最小購(gòu)電成本為目標(biāo)函數(shù)，利用最優(yōu)化的方法，我們得到了各機(jī)組的出力分配預(yù)案（如下表所示）。機(jī)組12345678出力1507918099.512514095113.9此時(shí)各機(jī)組所對(duì)應(yīng)的潮流值（如下表所示）線路123456潮流值173.315796.8095-150.5051101.8151136.8380168.5319對(duì)于第四問(wèn)，利用第一問(wèn)求得的結(jié)果，檢驗(yàn)第三問(wèn)各機(jī)組的出力分配預(yù)案是否會(huì)引起輸電阻塞，如果引起輸電阻塞，根據(jù)輸電阻塞管理原則并考慮經(jīng)濟(jì)原則，用計(jì)算機(jī)編程搜索調(diào)整預(yù)案，得出阻塞費(fèi)用為405.8004元，新的方案見(jiàn)下表機(jī)組12345678出力1538822890.177515

3、29560.1116.1225此時(shí)各線路所對(duì)應(yīng)的潮流值（如下表所示）線路123456潮流值165100.6408-154.8553105.5524131.5169159.5725 對(duì)于第五問(wèn)，重復(fù)三，四問(wèn)的工作，我們得出方案（如下表所示），阻塞費(fèi)用為94.0377元。機(jī)組12345678出力1538822899.515215560.3117此時(shí)各線路對(duì)應(yīng)的潮流值（如下表所示）線路123456潮流值173.461994.1143-154.9968104.9631135.3452160.4779問(wèn)題的重述：電網(wǎng)公司在組織交易，調(diào)度和配送時(shí),必須遵循電網(wǎng)“安全第一”的原則，同時(shí)制定一個(gè)電力市場(chǎng)交易規(guī)

4、則，按照購(gòu)電費(fèi)用最小的經(jīng)濟(jì)目標(biāo)來(lái)運(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ì)劃的過(guò)程中，還需實(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ì)值可以超過(guò)限值的百分比的上限）。如果各機(jī)組出力分配方案使某條線路上的有功潮流的絕對(duì)值超出限值，稱為輸電阻塞。當(dāng)發(fā)生輸電阻塞時(shí)，需要研

5、究如何制訂既安全又經(jīng)濟(jì)的調(diào)度計(jì)劃。 電力市場(chǎng)交易規(guī)則：1. 以15分鐘為一個(gè)時(shí)段組織交易，每臺(tái)機(jī)組在當(dāng)前時(shí)段開(kāi)始時(shí)刻前給出下一個(gè)時(shí)段的報(bào)價(jià)。各機(jī)組將可用出力由低到高分成至多10段報(bào)價(jià)，每個(gè)段的長(zhǎng)度稱為段容量，每個(gè)段容量報(bào)一個(gè)價(jià)（稱為段價(jià)），段價(jià)按段序數(shù)單調(diào)不減。在最低技術(shù)出力以下的報(bào)價(jià)一般為負(fù)值，表示愿意付費(fèi)維持發(fā)電以避免停機(jī)帶來(lái)更大的損失。2. 在當(dāng)前時(shí)段內(nèi)，市場(chǎng)交易-調(diào)度中心根據(jù)下一個(gè)時(shí)段的負(fù)荷預(yù)報(bào),每臺(tái)機(jī)組的報(bào)價(jià)、當(dāng)前出力和出力改變速率，按段價(jià)從低到高選取各機(jī)組的段容量或其部分（見(jiàn)下面注釋），直到它們之和等于預(yù)報(bào)的負(fù)荷，這時(shí)每個(gè)機(jī)組被選入的段容量或其部分之和形成該時(shí)段該機(jī)組的出力分配預(yù)案

6、（初始交易結(jié)果）。最后一個(gè)被選入的段價(jià)（最高段價(jià)）稱為該時(shí)段的清算價(jià)，該時(shí)段全部機(jī)組的所有出力均按清算價(jià)結(jié)算。注釋：（a） 每個(gè)時(shí)段的負(fù)荷預(yù)報(bào)和機(jī)組出力分配計(jì)劃的參照時(shí)刻均為該時(shí)段結(jié)束時(shí)刻。（b） 機(jī)組當(dāng)前出力是對(duì)機(jī)組在當(dāng)前時(shí)段結(jié)束時(shí)刻實(shí)際出力的預(yù)測(cè)值。（c） 假設(shè)每臺(tái)機(jī)組單位時(shí)間內(nèi)能增加或減少的出力相同，該出力值稱為該機(jī)組的爬坡速率。由于機(jī)組爬坡速率的約束，可能導(dǎo)致選取它的某個(gè)段容量的部分。（d） 為了使得各機(jī)組計(jì)劃出力之和等于預(yù)報(bào)的負(fù)荷需求，清算價(jià)對(duì)應(yīng)的段容量可能只選取部分。市場(chǎng)交易-調(diào)度中心在當(dāng)前時(shí)段內(nèi)要完成的具體操作過(guò)程如下：1、 監(jiān)控當(dāng)前時(shí)段各機(jī)組出力分配方案的執(zhí)行，調(diào)度AGC輔助服

7、務(wù)，在此基礎(chǔ)上給出各機(jī)組的當(dāng)前出力值。2、 作出下一個(gè)時(shí)段的負(fù)荷需求預(yù)報(bào)。3、 根據(jù)電力市場(chǎng)交易規(guī)則得到下一個(gè)時(shí)段各機(jī)組出力分配預(yù)案。4、 計(jì)算當(dāng)執(zhí)行各機(jī)組出力分配預(yù)案時(shí)電網(wǎng)各主要線路上的有功潮流，判斷是否會(huì)出現(xiàn)輸電阻塞。如果不出現(xiàn)，接受各機(jī)組出力分配預(yù)案；否則，按照如下原則實(shí)施阻塞管理：輸電阻塞管理原則：（1） 調(diào)整各機(jī)組出力分配方案使得輸電阻塞消除。（2） 如果（1）做不到，還可以使用線路的安全裕度輸電，以避免拉閘限電（強(qiáng)制減少負(fù)荷需求），但要使每條線路上潮流的絕對(duì)值超過(guò)限值的百分比盡量小。（3） 如果無(wú)論怎樣分配機(jī)組出力都無(wú)法使每條線路上的潮流絕對(duì)值超過(guò)限值的百分比小于相對(duì)安全裕度，則必

8、須在用電側(cè)拉閘限電。（4） 當(dāng)改變根據(jù)電力市場(chǎng)交易規(guī)則得到的各機(jī)組出力分配預(yù)案時(shí)，一些通過(guò)競(jìng)價(jià)取得發(fā)電權(quán)的發(fā)電容量（稱序內(nèi)容量）不能出力；而一些在競(jìng)價(jià)中未取得發(fā)電權(quán)的發(fā)電容量（稱序外容量）要在低于對(duì)應(yīng)報(bào)價(jià)的清算價(jià)上出力。因此，發(fā)電商和網(wǎng)方將產(chǎn)生經(jīng)濟(jì)利益沖突。網(wǎng)方應(yīng)該為因輸電阻塞而不能執(zhí)行初始交易結(jié)果付出代價(jià)，網(wǎng)方在結(jié)算時(shí)應(yīng)該適當(dāng)?shù)亟o發(fā)電商以經(jīng)濟(jì)補(bǔ)償，由此引起的費(fèi)用稱之為阻塞費(fèi)用。網(wǎng)方在電網(wǎng)安全運(yùn)行的保證下應(yīng)當(dāng)同時(shí)考慮盡量減少阻塞費(fèi)用?，F(xiàn)要解決的問(wèn)題如下：（表格見(jiàn)附錄1）1. 某電網(wǎng)有8臺(tái)發(fā)電機(jī)組，6條主要線路，表1和表2中的方案0給出了各機(jī)組的當(dāng)前出力和各線路上對(duì)應(yīng)的有功潮流值，方案132給出

9、了圍繞方案0的一些實(shí)驗(yàn)數(shù)據(jù)，試用這些數(shù)據(jù)確定各線路上有功潮流關(guān)于各發(fā)電機(jī)組出力的近似表達(dá)式。2. 設(shè)計(jì)一種簡(jiǎn)明、合理的阻塞費(fèi)用計(jì)算規(guī)則，除考慮上述電力市場(chǎng)規(guī)則外，還需注意：在輸電阻塞發(fā)生時(shí)公平地對(duì)待序內(nèi)容量不能出力的部分和報(bào)價(jià)高于清算價(jià)的序外容量出力的部分。3. 假設(shè)下一個(gè)時(shí)段預(yù)報(bào)的負(fù)荷需求是982.4MW，表3、表4和表5分別給出了各機(jī)組的段容量、段價(jià)和爬坡速率的數(shù)據(jù)，試按照電力市場(chǎng)規(guī)則給出下一個(gè)時(shí)段各機(jī)組的出力分配預(yù)案。4. 按照表6給出的潮流限值，檢查得到的出力分配預(yù)案是否會(huì)引起輸電阻塞，并在發(fā)生輸電阻塞時(shí)，根據(jù)安全且經(jīng)濟(jì)的原則，調(diào)整各機(jī)組出力分配方案，并給出與該方案相應(yīng)的阻塞費(fèi)用。5.

10、 假設(shè)下一個(gè)時(shí)段預(yù)報(bào)的負(fù)荷需求是1052.8MW，重復(fù)34的工作。問(wèn)題的分析：1 初步分析題意，根據(jù)表1和表2所給的數(shù)據(jù)，我們利用數(shù)理統(tǒng)計(jì)知識(shí)和線性回歸分析方法，擬建立多元線性均值回歸模型，推算出各線路上有功潮流關(guān)于各發(fā)電機(jī)組出力的近似表達(dá)式；2 設(shè)計(jì)簡(jiǎn)明，合理的阻塞費(fèi)用的計(jì)算規(guī)則，關(guān)鍵在于體現(xiàn)網(wǎng)方是怎樣充分補(bǔ)償發(fā)電廠的損失的，在這里，我們引入基本補(bǔ)償價(jià)的概念，利用限上和限下的方法得出阻塞費(fèi)用的計(jì)算規(guī)則；3 僅僅考慮電力市場(chǎng)交易規(guī)則和各機(jī)組的爬坡速率的影響，要得出下一時(shí)段各機(jī)組的出力分配預(yù)案，就是在各機(jī)組的總段容量不超出爬坡速率所限制的極限值的前提條件下，選取價(jià)格最低的全部機(jī)組的段容量總和，使

11、得這個(gè)總和等于負(fù)荷需求； 4 用第一問(wèn)所得出來(lái)的有功潮流的表達(dá)式檢查預(yù)案是否會(huì)引起輸電阻塞，要是會(huì)發(fā)生輸電阻塞，則根據(jù)輸電阻塞管理原則并考慮電力市場(chǎng)規(guī)則的情況下調(diào)整預(yù)案得出一個(gè)新的方案，并根據(jù)設(shè)計(jì)的阻塞費(fèi)用計(jì)算規(guī)則求出新方案對(duì)于預(yù)案的阻塞費(fèi)用。模型的假設(shè)：1 各線路上的有功潮流只受到這個(gè)電網(wǎng)里面八個(gè)發(fā)電機(jī)組所發(fā)的電量的影響，與外界的電網(wǎng)無(wú)關(guān)，即電網(wǎng)結(jié)構(gòu)不會(huì)改變；2 發(fā)電機(jī)組出力是按照爬坡速率一直增加或減少到符合出力分配方案的要求后，出力就是一直恒定不變的；3 發(fā)電機(jī)出力的增量是由低價(jià)格的容量增加到高價(jià)格的容量的，而出力的減量就是由高價(jià)格的容量減少到低價(jià)格的容量；符號(hào)的說(shuō)明： ： 第i條線路的潮

12、流值（i1，2，6）； ： 序外容量要出力的段容量或其部分； ： 序內(nèi)容量不能出力的段容量或其部分； ： 第i個(gè)機(jī)組的出力（i1，2，8）； ： 當(dāng)前方案0中第i個(gè)機(jī)組的出力（i1，2，8）； ： 第i個(gè)機(jī)組在單位時(shí)間內(nèi)的功率增量（i1，2，8）； ： 第i個(gè)機(jī)組的爬坡速率（i1，2，8）； ： 一個(gè)時(shí)間段（15分鐘）； ： 對(duì)應(yīng)段容量/對(duì)應(yīng)段容量所在機(jī)組的爬坡速率； ： 對(duì)應(yīng)段容量的前面容量總和變化到方案0的出力所需的時(shí)間；模型的建立：?jiǎn)栴}1 多元線性均值回歸模型 根據(jù)題中所給統(tǒng)計(jì)數(shù)據(jù)，我們用回歸分析方法，考慮線性分析模型，討論得到：運(yùn)用全回歸分析雖然計(jì)算簡(jiǎn)單，但是最終的計(jì)算結(jié)果誤差比較大。

13、所以我們建立多元線性均值回歸模型。則有： （1）對(duì)（1）式的算法如下： （2） （3） 由表1和表2的132方案中知，為稀疏矩陣，元素只有一項(xiàng)i不為零（此項(xiàng)機(jī)組表現(xiàn)為相對(duì)0方案的變動(dòng)），且同類的有4組，剛好8個(gè)機(jī)組的單獨(dú)變動(dòng)都有4組數(shù)，很明顯我們可以均值就可求出i組的k的列向量，運(yùn)用循環(huán)即可全部求得8組k的列向量，合并得k值。（程序附錄2和3）系數(shù)矩陣k觀察以上求出的系數(shù)k可以發(fā)現(xiàn)：k值比較小，而且相鄰k值的變化也比較小。 求c時(shí)（程序見(jiàn)附錄4），利用以上求得k值，有： （4） 與對(duì)應(yīng)于132方案中的某一方案，由此運(yùn)用循環(huán)當(dāng)求36組，求均值可得。常數(shù)c所以寫出y和p的近似表達(dá)式為：對(duì)上述方法進(jìn)

14、行誤差分析，由求出的表達(dá)式返回32組方案與原方案對(duì)比，得出y的相對(duì)誤差圖（如圖一所示）： 如圖一可知，此模型的誤差較小。問(wèn)題2對(duì)于第二問(wèn)，考慮的是當(dāng)阻塞發(fā)生的時(shí)候網(wǎng)方所要補(bǔ)償給發(fā)電商的費(fèi)用，這個(gè)費(fèi)用稱為阻塞費(fèi)用S，一共包括兩個(gè)部分：一部分是序內(nèi)容量不能出力的部分，決定這部分阻塞費(fèi)用的主要因素是不能出力的機(jī)組停發(fā)或少發(fā)電的度數(shù)以及少發(fā)一度電所需的補(bǔ)償價(jià)格；另一部分是報(bào)價(jià)高于清算價(jià)的序外容量出力的部分，決定這部分阻塞費(fèi)用的主要因素是序外容量機(jī)組多發(fā)電的度數(shù)以及多發(fā)一度電所需的補(bǔ)償價(jià)格。模型一：限上和限下的方法根據(jù)上面的思想，我們應(yīng)用限上和限下的方法1來(lái)計(jì)算阻塞費(fèi)用，規(guī)則如下：1 當(dāng)出現(xiàn)輸電阻塞時(shí)，

15、在競(jìng)價(jià)中未取得發(fā)電權(quán)的發(fā)電容量（序外容量）按報(bào)價(jià)由低到高的順序安排機(jī)組發(fā)電(限上),在競(jìng)價(jià)中本已經(jīng)取得發(fā)電權(quán)的發(fā)電容量（序內(nèi)容量）按報(bào)價(jià)由高到低的順序安排機(jī)組停發(fā)或少發(fā)電(限下)。2 在競(jìng)價(jià)中未取得發(fā)電權(quán)的發(fā)電容量（序外容量）是在低于其對(duì)應(yīng)報(bào)價(jià)的清算 價(jià)上出力的，這樣就會(huì)使得發(fā)電商利益虧損而網(wǎng)方從中盈利，網(wǎng)方應(yīng)該對(duì)這 些發(fā)電商給與經(jīng)濟(jì)補(bǔ)償，這里補(bǔ)償我們就用對(duì)應(yīng)段容量的報(bào)價(jià)和清算價(jià)之差作為補(bǔ)償價(jià)格；而那些在競(jìng)價(jià)中本已經(jīng)取得發(fā)電權(quán)的發(fā)電容量（序內(nèi)容量）卻不能出力的發(fā)電商，由于不能出力或部分不能出力，使得這些發(fā)電商由于網(wǎng)方線路原因造成經(jīng)濟(jì)利益的損失，這時(shí)網(wǎng)方也應(yīng)該給與這些發(fā)電商經(jīng)濟(jì)補(bǔ)償，這里我們就以

16、清算價(jià)和其不能出力的對(duì)應(yīng)的段容量的報(bào)價(jià)之差作為補(bǔ)償價(jià)格。根據(jù)計(jì)算規(guī)則，我們可以列出阻塞費(fèi)用S的式子： （5）其中 清算價(jià)對(duì)應(yīng)段容量的報(bào)價(jià)對(duì)應(yīng)段容量的報(bào)價(jià)清算價(jià)模型二：模型二是對(duì)于模型一所進(jìn)行的改進(jìn)，在模型一中可以發(fā)現(xiàn)：對(duì)于序內(nèi)容量，當(dāng)某個(gè)的段容量的報(bào)價(jià)等于它的清算價(jià)的時(shí)候，此時(shí)0，這表明這個(gè)段容量是得不到經(jīng)濟(jì)補(bǔ)償?shù)?。由于每一個(gè)發(fā)電廠發(fā)的每度電都存在一個(gè)利潤(rùn)空間，這個(gè)利潤(rùn)是包含在報(bào)價(jià)里面的，按照模型一的結(jié)論表明發(fā)電廠按清算價(jià)作為報(bào)價(jià)的段容量的利潤(rùn)得不到經(jīng)濟(jì)補(bǔ)償，顯然這是不合理的。為了解決這個(gè)問(wèn)題，我們?cè)谶@里引入一個(gè)概念：基本補(bǔ)償價(jià)B，根據(jù)文獻(xiàn)資料 2 ，我們定 （6），于是我們就以清算價(jià)和其不能

17、出力的對(duì)應(yīng)的段容量的報(bào)價(jià)之差再加上基本補(bǔ)償價(jià)作為序內(nèi)容量不能出力的發(fā)電商的經(jīng)濟(jì)補(bǔ)償。于是式子（5）變?yōu)椋?（7）問(wèn)題3由題意可知: （8）在下一個(gè)時(shí)間段的結(jié)束時(shí)刻各機(jī)組的出力大小為 （9）根據(jù)式子（8）（9），結(jié)合表1中方案0各機(jī)組的數(shù)據(jù)和表5中各機(jī)組的爬坡速率，相對(duì)于方案0，我們可以算出在下一個(gè)時(shí)間段的結(jié)束時(shí)刻各機(jī)組所允許能選取的出力范圍（如下表31所示）表31 各機(jī)組所允許能選取的出力范圍機(jī)組12345678最大出力1538822899.5152155102.1117最少出力875813260.5989560.163根據(jù)電力市場(chǎng)規(guī)則，我們對(duì)附錄1中的表3，表4進(jìn)行處理。（按照段價(jià)由低到高排

18、列如下表格群所示）表格群32 價(jià)格-800-610-607-590-560-505-500-500機(jī)組83651174段容量70110957570705055價(jià)格00000116120124機(jī)組12356571段容量00050151550價(jià)格146150152153159168170173機(jī)組54386146段容量054001001020價(jià)格180182183188189200203205機(jī)組72853426段容量520201501080價(jià)格210215233233245250251252機(jī)組15382571段容量0153001501530價(jià)格252253255258260283300302機(jī)

19、組68437824段容量1520100100610價(jià)格303305306308310312315318機(jī)組86735178段容量20101020100510價(jià)格320325330335348356360363機(jī)組24177321段容量21501034000價(jià)格380380396400405410415435機(jī)組64586234段容量20010150000價(jià)格489495500510520548800800機(jī)組12356784段容量408401010251按照以上表格，僅僅考慮電力市場(chǎng)規(guī)則，而不考慮各機(jī)組所允許的出力范圍的時(shí)候，選取到第7個(gè)表格中的303價(jià)格位的段容量就能夠滿足負(fù)荷需求（這個(gè)價(jià)格

20、位之前的段容量全部選?。?，這時(shí)候我們得到各個(gè)機(jī)組的出力分別為：表33機(jī)組12345678出力1507918010012514095113.4由于受到爬坡速率的限制，第四機(jī)組的出力最大值為99.5，不能達(dá)到100，將第四組不能達(dá)到的0.5段容量按照電力市場(chǎng)規(guī)則選取還未選取的且對(duì)應(yīng)段容量?jī)r(jià)格最低的，由此，我們得出下一個(gè)時(shí)段各機(jī)組的出力分配預(yù)案為（如下表34所示）：表34 下一時(shí)段各機(jī)組的出力分配預(yù)案（負(fù)荷需求為982.4MW時(shí)）機(jī)組12345678出力1507918099.512514095113.9問(wèn)題4 將上題所得的各機(jī)組的出力分配預(yù)案代入第一問(wèn)的模型二中求出的y和p的關(guān)系式中，得到各線路的潮

21、流值如下表41所示：表41 各線路的潮流值線路123456潮流值173.315796.8095-150.5051101.8151136.8380168.5319上表41與附錄1中的表6相比較可以發(fā)現(xiàn)：線路1，5，6會(huì)引起輸電阻塞。 圖 二根據(jù)算法（如圖二所示）編寫程序（見(jiàn)附錄6），由計(jì)算機(jī)搜索得到：改變后的出力分配方案如下表42所示：表42 各機(jī)組的出力分配方案機(jī)組12345678出力15388228901116.1225此時(shí)各線路的潮流值如下表43所示：表43 改變方案后各線路的潮流值線路123456潮流值165100.6408-154.8553105.552413

22、1.5169159.5725并得到對(duì)應(yīng)的阻塞費(fèi)用為405.8004元。問(wèn)題5 利用第三問(wèn)和第四問(wèn)同樣的方法，得出在下一個(gè)時(shí)段預(yù)報(bào)的負(fù)荷需求是1052.8MW時(shí)，各機(jī)組的出力分配預(yù)案如下（如下表51所示）表51 下一時(shí)段各機(jī)組的出力分配預(yù)案（負(fù)荷需求為1052.8MW時(shí)）機(jī)組12345678出力15081218.299.5135150102.1117此時(shí)各線路的潮流值如下表52所示表52 各線路的潮流值線路123456潮流值177.260894.6307-155.7527110.0612134.8521167.0870調(diào)整之后各機(jī)組的出力分配方案為：（如表53所示） 表53 調(diào)整后的各機(jī)組出力分

23、配方案機(jī)組12345678出力1538822899.515215560.3117各線路的潮流值為：（如表54所示） 表54 調(diào)整后的各線路的潮流值線路123456潮流值173.461994.1143-154.9968104.9631135.3452160.4779從表54可以知道，第1，5條線路都是超過(guò)限值的，但是都在安全裕度內(nèi)，第一條線路超出限值的百分比為5.128，第五條線路超出限值的百分比為2.5342。并得出相應(yīng)的阻塞費(fèi)用為：94.0377元。模型的優(yōu)缺點(diǎn)和改進(jìn)方向：優(yōu)缺點(diǎn)：本論文中的平均值分析法能在一定范圍內(nèi)很好地吻合數(shù)據(jù)，經(jīng)誤差分析亦可得相對(duì)誤差0.0001左右，而且機(jī)組系數(shù)比較小

24、，能很好地說(shuō)明此電網(wǎng)有比較好的穩(wěn)定性，不會(huì)因一個(gè)機(jī)組的小變化而有大的改動(dòng)，但是，此模型只適應(yīng)機(jī)組出力都較大的情況下，在所有機(jī)組都比較小時(shí)，此模型有較大誤差，但此范圍并不在此題的計(jì)算范圍內(nèi)，此模型已經(jīng)能非常很好地解決一般的電網(wǎng)調(diào)動(dòng)。另外，本文中的計(jì)算較多應(yīng)用matlab編程，對(duì)大量的數(shù)據(jù)計(jì)算提供很好的橋梁，針對(duì)每一個(gè)問(wèn)題都有一個(gè)程序來(lái)解決，為以后的利用提供方便。橋梁改進(jìn)方向：由于給出數(shù)據(jù)中有突變數(shù)據(jù)，所以在計(jì)算得出的相對(duì)誤差中有6個(gè)是大于0.05的，當(dāng)踢除這些值代以估算的均值，能更完美地完成數(shù)據(jù)的擬合。對(duì)于賠償算法中有發(fā)電商的隱藏利潤(rùn)，此利潤(rùn)應(yīng)該賠償過(guò)加入更加合理準(zhǔn)確，對(duì)此中的專業(yè)估計(jì)，是改進(jìn)的

25、重點(diǎn)。參考文獻(xiàn)：1 柯進(jìn),管霖 電力市場(chǎng)下的輸電阻塞管理技術(shù) 電力系統(tǒng)自動(dòng)化 第26卷 第14期 第四頁(yè) 2002年2 西部水利水電投資項(xiàng)目介紹29/xbshj/xbbshd/xbbshd_indx.htm 訪問(wèn)時(shí)間：2004-9-19 韓忠誠(chéng) 呂伯雄 金昌市招標(biāo)項(xiàng)目 訪問(wèn)時(shí)間：2004-9-19附錄1：表1 各機(jī)組出力方案 （單位：兆瓦，記作MW）方案機(jī)組123 4 5 6780120731808012512581.1901133.02731808012512581.1902129.63731808012512581.1903158.7773180801

26、2512581.1904145.32731808012512581.190512078.5961808012512581.190612075.451808012512581.190712090.4871808012512581.190812083.8481808012512581.190912073231.398012512581.1901012073198.488012512581.1901112073212.648012512581.1901212073190.558012512581.190131207318075.85712512581.190141207318065.95812512

27、581.190151207318087.25812512581.190161207318097.82412512581.190171207318080150.7112581.190181207318080141.5812581.190191207318080132.3712581.190201207318080156.9312581.190211207318080125138.8881.190221207318080125131.2181.190231207318080125141.7181.190241207318080125149.2981.19025120731808012512560.

28、5829026120731808012512570.9629027120731808012512564.8549028120731808012512575.5299029120731808012512581.1104.8430120731808012512581.1111.2231120731808012512581.198.09232120731808012512581.1120.44表2 各線路的潮流值（各方案與表1相對(duì)應(yīng)，單位：MW）方案線路1234560164.78140.87-144.25119.09135.44157.691165.81140.13-145.14118.63135.

29、37160.762165.51140.25-144.92118.7135.33159.983167.93138.71-146.91117.72135.41166.814166.79139.45-145.92118.13135.41163.645164.94141.5-143.84118.43136.72157.226164.8141.13-144.07118.82136.02157.57165.59143.03-143.16117.24139.66156.598165.21142.28-143.49117.96137.98156.969167.43140.82-152.26129.58132.

30、04153.610165.71140.82-147.08122.85134.21156.2311166.45140.82-149.33125.75133.28155.0912165.23140.85-145.82121.16134.75156.7713164.23140.73-144.18119.12135.57157.214163.04140.34-144.03119.31135.97156.3115165.54141.1-144.32118.84135.06158.2616166.88141.4-144.34118.67134.67159.2817164.07143.03-140.9711

31、8.75133.75158.8318164.27142.29-142.15118.85134.27158.3719164.57141.44-143.3119134.88158.0120163.89143.61-140.25118.64133.28159.1221166.35139.29-144.2119.1136.33157.5922165.54140.14-144.19119.09135.81157.6723166.75138.95-144.17119.15136.55157.5924167.69138.07-144.14119.19137.11157.6525162.21141.21-14

32、4.13116.03135.5154.2626163.54141-144.16117.56135.44155.9327162.7141.14-144.21116.74135.4154.8828164.06140.94-144.18118.24135.4156.6829164.66142.27-147.2120.21135.28157.6530164.7142.94-148.45120.68135.16157.6331164.67141.56-145.88119.68135.29157.6132164.69143.84-150.34121.34135.12157.64表3 各機(jī)組的段容量 （單位

33、：MW）機(jī)組段1234567891017005000300004023002081562008311004003002040040455510101010150015755150151501010106950102001510200107501551510105103287002002002010155表4 各機(jī)組的段價(jià)（單位：元/兆瓦小時(shí)，記作元/MWh）機(jī)組段123456789101-50501241682102523123303634892-56001822032453003203604104953-61001521892332583083564155004-50015017020025

34、53023253804358005-59001161461882152503103965106-60701591732052523053804055207-5001201802512603063153353485488-800153183233253283303318400800表5 各機(jī)組的爬坡速率 （單位：MW/分鐘）機(jī)組12345678速率1.821.41.8表6 各線路的潮流限值（單位：MW）和相對(duì)安全裕度線路123456限值165150160155132162安全裕度13%18%9%11%15%14%附錄2function k=zhen(fa8,r6) %求k初值

35、 p0=164.78 140.87 -144.25 119.09 135.44 157.69;pj=120 73 180 80 125 125 81.1 90;j=0;for i=1:8 j=fa8(i)-pj(i)+j;endpi=r6-p0;k=pi./j;附錄3function ok=zhen1(fa8,r6) %求k系數(shù)clcfor i=1:32 a=fa8(i,:);b=r6(i,:); k(i,:)=zhen(a,b);endoptions=;for i=1:8 for j=1:6 ok(i,j)=mean(k(4*i,j),k(4*i-1,j),k(4*i-2,j),k(4*i-

36、3,j); end end附錄4function c=zhen4(ok,fa8,r6) %求c值for k=1:32 if k=11 c(k,:)=c(k-1,:); end for j=1:6 c(k,j)=r6(k,j)-fa8(k,:)*ok(:,j); endendc=mean(c);附錄5function pout=zhen5(pin) %求潮流值ok=0.08284 0.04828 0.05297 0.11993 -0.02544 0.12201 0.12158 -0.00123;. -0.05456 0.12785 -0.00003 0.03328 -0.08685 -0.112

37、44 -0.01893 -0.09873;. -0.06954 0.06105 -0.15662 -0.00992 0.12449 0.00212 0.00251 -0.20139;. -0.03446 -0.10241 0.20516 -0.02083 -0.01183 -0.00595 0.14492 -0.07655;. 0.00053 0.24329 -0.06455 -0.04113 -0.06522 0.07034 -0.00426 -0.00891;. 0.23781 -0.06017 -0.07787 0.09298 0.04690 0.00008 0.16593 0.0006

38、9;c=110.2965 131.2289 -108.873 77.4817 132.97447 120.6633;pout=ok*pin'+c'附錄6function f=zhenfun2(x)x0=120 73 180 80 125 125 81.1 90;x1=150 81 218.2 99.5 135 150 102.1 117;v=2.2 1 3.2 1.3 1.8 2 1.4 1.8;duan=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

39、,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; price=-505 0 124 168 210 252 312 330 363 489;. -560 0 182 203 245 300 320 360 410 495;. -610 0

40、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; n0=6 6 5 6 6 6 5 10;b=0 0 0 0 0 0 0 0 ;p=0 0 0 0 0 0 0 0

41、 ;for i=1:8 if x(i)-x0(i)<0 xi=x(i)-x0(i); t=0; if xi<0 for n=n0(i):-1:1 b(i)=duan(i,n)+b(i); if b(i)+xi>=0 b(i)=-duan(i,n)+b(i); break; end end for ni=n0(i):-1:n if ni=n if duan(i,ni)=0 p(i)=-duan(i,ni)2/v(i)/2*(303-price(i,ni)+p(i); t=t+duan(i,ni)/v(i); end else if xi<=x0(i)-x1(i) p(i)

42、=-(xi+b(i)2/v(i)/2*(303-price(i,ni)+(0.25-t-(x1(i)-x0(i)/v(i)*(xi+b(i)*(303-price(i,ni)+p(i); else p(i)=-(xi+b(i)2/v(i)/2*(303-price(i,ni)+(0.25-t-(x1(i)+xi-x0(i)/v(i)*(xi+b(i)*(303-price(i,ni)+p(i); end end end else for n=n0(i):10 b(i)=duan(i,n)+b(i); if b(i)-xi>=0 b(i)=-duan(i,n)+b(i); break; end end for ni=n0(i):n if ni=n if duan(i,ni)=0 p(i)=duan(i,ni)2/v(i)/2*(price(i,ni)-303)+p(i); t=t+duan(i,ni)/v(i); end else p(i)=(xi-b(i)2/v(i)/2*(price(i,ni)-303)+(0.25-t)*(xi-b(i)*(303-price(i,ni)+p(i); end endend elseif x(i)-x1(i)<0 xi=x(i)-x1(i); t=0; if xi<0 for n=n0(i):-1:1 b(