張焱，王艷，王磊

（1.東南大學(xué)電氣工程學(xué)院，江蘇南京210096；2.常州供電公司，江蘇常州213003）

動(dòng)態(tài)優(yōu)化調(diào)度是研究在連續(xù)多個(gè)時(shí)間段內(nèi)、機(jī)組運(yùn)行狀態(tài)已經(jīng)確定的情況下，如何有效地調(diào)整各個(gè)機(jī)組在不同時(shí)段的有功功率輸出，使得在滿(mǎn)足各種靜態(tài)約束和動(dòng)態(tài)約束下，整個(gè)調(diào)度周期內(nèi)的某些目標(biāo)函數(shù)達(dá)到最優(yōu)。由于動(dòng)態(tài)優(yōu)化調(diào)度問(wèn)題的規(guī)模隨著系統(tǒng)規(guī)模以及時(shí)段數(shù)的增加而迅速增長(zhǎng)，因此提高算法的效率是迫切需要解決的問(wèn)題［1］。

多年來(lái)，學(xué)者們通過(guò)不斷探索，提出了各種求解優(yōu)化調(diào)度問(wèn)題的計(jì)算方法。其中，簡(jiǎn)化梯度法具有簡(jiǎn)單、容易實(shí)現(xiàn)等優(yōu)點(diǎn)，但是收斂性較差，尤其是在接近最優(yōu)點(diǎn)附近時(shí)收斂很慢［2］；牛頓法和簡(jiǎn)化梯度法相比，收斂速度較快，具有二階收斂性，但是也存在著數(shù)值不穩(wěn)定、不能總是有效地確定迭代過(guò)程中起作用的約束集等缺點(diǎn)；人工智能算法對(duì)函數(shù)的性態(tài)要求比較低，近年來(lái)在最優(yōu)潮流領(lǐng)域中得到了廣泛的運(yùn)用［3］，但是存在著收斂速度比較慢的缺點(diǎn)。相比于以上算法，原對(duì)偶內(nèi)點(diǎn)法具有迭代次數(shù)與系統(tǒng)規(guī)模關(guān)系不大、收斂性好、計(jì)算速度快等優(yōu)點(diǎn)。文獻(xiàn)［4］將原對(duì)偶內(nèi)點(diǎn)法應(yīng)用于非線(xiàn)性規(guī)模的無(wú)功優(yōu)化調(diào)度問(wèn)題，充分顯示了內(nèi)點(diǎn)法的收斂特性；文獻(xiàn)［5］則采用線(xiàn)性約束網(wǎng)絡(luò)流來(lái)描述電力系統(tǒng)動(dòng)態(tài)經(jīng)濟(jì)調(diào)度問(wèn)題，并用內(nèi)點(diǎn)法進(jìn)行求解，克服了傳統(tǒng)方法難以精確處理不等式約束的弊端，使計(jì)算結(jié)果更精確；文獻(xiàn)［6］則采用引入離散懲罰的非線(xiàn)性原對(duì)偶內(nèi)點(diǎn)法求解動(dòng)態(tài)無(wú)功優(yōu)化模型，具有較快的計(jì)算速度。

隨著電力行業(yè)節(jié)能減排政策的實(shí)施，在動(dòng)態(tài)優(yōu)化調(diào)度中同時(shí)將煤耗和污染物排放一起優(yōu)化更具現(xiàn)實(shí)意義［7］。本文首先建立了同時(shí)考慮調(diào)度周期內(nèi)煤耗量、污染物氣體排放量最小的多目標(biāo)動(dòng)態(tài)優(yōu)化調(diào)度模型，并利用模糊集理論使之轉(zhuǎn)化為求取滿(mǎn)意度最高的單目標(biāo)問(wèn)題；接著推導(dǎo)出了求解該模型的原對(duì)偶解耦內(nèi)點(diǎn)算法；最后對(duì)算例結(jié)果進(jìn)行了比較分析。

1 基于節(jié)能減排的多目標(biāo)動(dòng)態(tài)優(yōu)化調(diào)度模型

1.1 多目標(biāo)動(dòng)態(tài)優(yōu)化調(diào)度模型

1.1.1 目標(biāo)函數(shù)

本文以常規(guī)燃煤機(jī)組的煤耗和污染物氣體排放最小作為目標(biāo)函數(shù)：

式(1，2)中：T為調(diào)度周期；Ng表示發(fā)電機(jī)臺(tái)數(shù)；PGi(t)為第i臺(tái)發(fā)電機(jī)組t時(shí)段的有功出力；ai，bi，ci為煤耗系數(shù)；αi，βi，γi為污染物排放系數(shù)。

1.1.2 靜態(tài)等式約束

動(dòng)態(tài)優(yōu)化調(diào)度問(wèn)題的靜態(tài)等式約束即為各個(gè)時(shí)段的節(jié)點(diǎn)潮流平衡約束：

式中：省略下標(biāo)t；Vi，θi為節(jié)點(diǎn)電壓與相角；θij=θi-θj；Gij，Bij分別為系統(tǒng)的導(dǎo)納矩陣的實(shí)部與虛部；PDi，QDi為節(jié)點(diǎn)i的有功負(fù)荷與無(wú)功負(fù)荷；PGi，QGi為節(jié)點(diǎn)i的有功注入與無(wú)功注入。

1.1.3 靜態(tài)不等式約束

式中不等式約束依次為發(fā)電機(jī)有功出力約束、發(fā)電機(jī)無(wú)功出力約束、節(jié)點(diǎn)電壓約束、線(xiàn)路有功潮流約束。各個(gè)變量的上下限且分別用上劃線(xiàn)和下劃線(xiàn)表示。

1.1.4 動(dòng)態(tài)約束

本文考慮發(fā)電機(jī)爬坡約束：

式中：RampGi為發(fā)電機(jī)在相鄰時(shí)段間能增加或減小的最大功率。

1.2 多目標(biāo)動(dòng)態(tài)優(yōu)化調(diào)度的模糊模型

上述模型有2個(gè)相互競(jìng)爭(zhēng)的目標(biāo)函數(shù)，是多目標(biāo)規(guī)劃模型，有學(xué)者提出用模糊集理論來(lái)解決這個(gè)問(wèn)題，并取得較好效果。文獻(xiàn)［8］用模糊集理論將多目標(biāo)函數(shù)和部分可伸縮的約束條件模糊化，使最優(yōu)潮流問(wèn)題在更加符合實(shí)際情況的模型上實(shí)現(xiàn)優(yōu)化；文獻(xiàn)［9］將模糊集理論應(yīng)用于來(lái)水量不確定性分析，提出了一種水火電系統(tǒng)短期經(jīng)濟(jì)調(diào)度的新方法，較好地提高了系統(tǒng)的經(jīng)濟(jì)性；文獻(xiàn)［10］用多目標(biāo)模糊優(yōu)化理論和動(dòng)態(tài)規(guī)劃法解決火電廠(chǎng)多目標(biāo)負(fù)荷優(yōu)化問(wèn)題，得到科學(xué)、合理的負(fù)荷分配方案。因此，本文采用模糊集理論使多目標(biāo)動(dòng)態(tài)優(yōu)化調(diào)度轉(zhuǎn)化為單目標(biāo)問(wèn)題。在滿(mǎn)足所有約束條件的前提下，發(fā)電煤耗和污染物氣體排放總量越小越好，有上限而無(wú)下限，因此選用降半直線(xiàn)作為它們的隸屬度函數(shù)［11］，如圖1、圖2所示。

圖1 煤耗隸屬度函數(shù)

圖2 污染物氣體排放隸屬度函數(shù)

2個(gè)目標(biāo)隸屬度函數(shù)可分別用式（6，7）表示：式中：f1(x)，f2(x)分別為調(diào)度周期內(nèi)的煤耗、污染物氣體排放總量；c01，c02分別為單獨(dú)以煤耗、污染物氣體排放為優(yōu)化目標(biāo)時(shí)的結(jié)果；c01+δ01，c02+δ02分別為煤耗、污染物排放的最大可接受值。取λ為所有隸屬函數(shù)中最小的隸屬變量，可稱(chēng)之為滿(mǎn)意度。

根據(jù)最大隸屬度原理，可以將原來(lái)的多目標(biāo)動(dòng)態(tài)優(yōu)化調(diào)度問(wèn)題轉(zhuǎn)化為滿(mǎn)足所有約束條件的滿(mǎn)意度最大化的問(wèn)題。用以下非線(xiàn)性規(guī)劃模型表示。

上述模型包含T個(gè)時(shí)段，每個(gè)時(shí)段有n個(gè)變量、m個(gè)靜態(tài)等式約束、r個(gè)靜態(tài)不等式約束，各時(shí)段間有q個(gè)動(dòng)態(tài)不等式約束，p個(gè)含有λ的不等式。多目標(biāo)動(dòng)態(tài)優(yōu)化調(diào)度的模糊模型比傳統(tǒng)動(dòng)態(tài)優(yōu)化調(diào)度模型增加了p個(gè)含有λ的不等式，λ不屬于靜態(tài)變量，因此求解該模型時(shí)更為復(fù)雜。

2 原對(duì)偶解耦內(nèi)點(diǎn)算法

由于動(dòng)態(tài)優(yōu)化調(diào)度問(wèn)題系統(tǒng)規(guī)模很大，而原對(duì)偶內(nèi)點(diǎn)法具有迭代次數(shù)與系統(tǒng)規(guī)模關(guān)系不大、收斂性好、計(jì)算速度快等優(yōu)點(diǎn)，因此本文采用原對(duì)偶內(nèi)點(diǎn)法求解多目標(biāo)動(dòng)態(tài)優(yōu)化調(diào)度模型。根據(jù)多目標(biāo)動(dòng)態(tài)優(yōu)化調(diào)度模型的特點(diǎn)，對(duì)原對(duì)偶內(nèi)點(diǎn)算法作出相應(yīng)改進(jìn)，通過(guò)對(duì)算法中修正方程的降階與解耦，將動(dòng)態(tài)優(yōu)化調(diào)度的大規(guī)模線(xiàn)性修正方程的求解，等值解耦轉(zhuǎn)換成小規(guī)模的動(dòng)態(tài)變量求解和各時(shí)段靜態(tài)變量求解，從而提高算法的求解效率。

2.1 引入滿(mǎn)意度后的修正方程

式中：yt為各時(shí)段靜態(tài)等式約束對(duì)應(yīng)的拉格朗日乘子，各時(shí)段靜態(tài)不等式約束對(duì)應(yīng)的松弛變量為lt，ut＞0，拉格朗日乘子為zt＞0，wt＞0，對(duì)數(shù)壁壘參數(shù)μt＞0；動(dòng)態(tài)約束對(duì)應(yīng)的松弛變量為sl，su＞0，拉格朗日乘子為yl，yu＞0，對(duì)數(shù)壁壘參數(shù)μd＞0；含有λ的不等式的松弛變量為ll，uu＞0，拉格朗日乘子為zz＞0，ww＞0，對(duì)數(shù)壁壘參數(shù)。

為了表述方便，將各時(shí)段的優(yōu)化變量以及與靜態(tài)約束相關(guān)的拉格朗日乘子和松弛變量構(gòu)成靜態(tài)變量ρt=[xt，yt，zt，wt，lt，ut]T，每個(gè)時(shí)段的靜態(tài)變量一共構(gòu)成T組靜態(tài)變量；定義與動(dòng)態(tài)約束相關(guān)的拉格朗日乘子和松弛變量為動(dòng)態(tài)變量ρd=[λ，yu，yl，su，sl，ww，zz，uu，ll]T。

原對(duì)偶內(nèi)點(diǎn)法通過(guò)式(10)的最優(yōu)性KKT條件形成一組非線(xiàn)性方程，然后采用牛頓法迭代求解。其中，線(xiàn)性化牛頓修正方程可以表示如下：

式(11)中：Wt、Lt推導(dǎo)過(guò)程參照文獻(xiàn)［13］，修正方程的維數(shù)高達(dá)(4r+m+n)T+4(q+p)維，需進(jìn)行降階。

2.2 修正方程的降階

通過(guò)線(xiàn)性變換消去方程(11)中的Δzt，Δwt，Δlt，Δut，Δsu，Δsl，Δuu，Δll后得到：

式中：ΔρR1=[Δxt，Δyt]T，ΔρRd=[Δλ，Δyu，Δyl，Δww，Δzz]T，各子矩陣WRT、LRT推導(dǎo)過(guò)程也可參照文獻(xiàn)[13]。求出ΔρRT，ΔρRd后，可以通過(guò)線(xiàn)性變換求出其他變量。降階后的修正方程維數(shù)為(m+n)T+2(q+p)，大大減少了方程規(guī)模和計(jì)算量。

2.3 修正方程的解耦

由多目標(biāo)動(dòng)態(tài)優(yōu)化調(diào)度的模型可知，各時(shí)段的靜態(tài)約束和目標(biāo)函數(shù)相互獨(dú)立，相應(yīng)地靜態(tài)變量ρt相互獨(dú)立，彼此不相關(guān)，從而使得式(11)中常數(shù)項(xiàng)Lt只與t時(shí)段的靜態(tài)變量有關(guān)，與動(dòng)態(tài)變量無(wú)關(guān)，且Δρt之間的關(guān)聯(lián)系數(shù)矩陣為零矩陣；另一方面，動(dòng)態(tài)約束使部分靜態(tài)變量（各時(shí)段的發(fā)電機(jī)有功出力）與動(dòng)態(tài)變量相互關(guān)聯(lián)，因此式(11)中Δρt與Δρd之間的關(guān)聯(lián)系數(shù)矩陣Et非0。由此可以得出，多目標(biāo)動(dòng)態(tài)優(yōu)化調(diào)度問(wèn)題的特殊耦合關(guān)系決定了式(11)的系數(shù)矩陣具有分塊對(duì)角帶邊結(jié)構(gòu)，降階KKT方程(12)的系數(shù)矩陣也同樣具有分塊對(duì)角帶邊結(jié)構(gòu)。針對(duì)這種特性，可以通過(guò)線(xiàn)性變換將式(12)等值簡(jiǎn)化成解耦方程［14］：

首先根據(jù)式(13)解出ΔρRd，然后分別求出各個(gè)時(shí)段的ΔρRt。解耦后，需要求解一個(gè)2(q+p)維方程和T個(gè)(m+n)維方程，計(jì)算量比直接求解一個(gè)(m+n)T+2(q+p)維方程要小得多。

2.4 相關(guān)參數(shù)設(shè)置

所有的對(duì)數(shù)壁壘參數(shù)都由式(15)求得，靜態(tài)變量中原變量、對(duì)偶變量對(duì)應(yīng)的步長(zhǎng)根據(jù)式(16)計(jì)算，各動(dòng)態(tài)變量對(duì)應(yīng)步長(zhǎng)αpd，αdd，αpλ，αdλ的計(jì)算方法同式(16)。

2.5 算法流程

算法流程如圖3所示。

3 算例分析

根據(jù)上述模型和算法，本文對(duì)IEEE-30節(jié)點(diǎn)測(cè)試系統(tǒng)進(jìn)行仿真計(jì)算與分析。

3.1 初始條件設(shè)置

調(diào)度周期取為一天，劃分成24個(gè)時(shí)段，每個(gè)時(shí)段為1h，圖4為負(fù)荷波動(dòng)系數(shù)曲線(xiàn)，假設(shè)系統(tǒng)有功功率日負(fù)荷波動(dòng)曲線(xiàn)和無(wú)功功率日負(fù)荷波動(dòng)曲線(xiàn)是相同的，并且各節(jié)點(diǎn)的負(fù)荷值在一天中以相同的負(fù)荷系數(shù)波動(dòng)。各常規(guī)發(fā)電機(jī)組的煤耗系數(shù)ai，bi，ci和污染物氣體排放系數(shù)αi，βi，γi如表1所示。本文以煤耗費(fèi)用（元）來(lái)表征煤耗量，以排放氣體的重量（t）來(lái)表征污染物氣體的排放量。

圖3 原對(duì)偶解耦內(nèi)點(diǎn)法的具體流程圖

表1 各常規(guī)機(jī)組煤耗和污染物氣體排放系數(shù)

圖4 系統(tǒng)日負(fù)荷波動(dòng)系數(shù)曲線(xiàn)

以發(fā)電機(jī)有功功率上限的百分比來(lái)表示爬坡約束，本文取為5%。

單獨(dú)優(yōu)化煤耗的結(jié)果為10502元，單獨(dú)優(yōu)化污染物氣體排放時(shí)結(jié)果為3660.2t，取各自的最大可接受值為單獨(dú)優(yōu)化結(jié)果的1.4倍。

3.2 動(dòng)態(tài)與靜態(tài)優(yōu)化調(diào)度結(jié)果對(duì)比分析

本文首先選取調(diào)度周期內(nèi)煤耗量最小為單目標(biāo)進(jìn)行優(yōu)化調(diào)度，圖5比較了當(dāng)爬坡約束為5%和未考慮爬坡約束時(shí)的27號(hào)發(fā)電機(jī)有功出力，即動(dòng)態(tài)優(yōu)化調(diào)度和靜態(tài)優(yōu)化調(diào)度時(shí)27號(hào)發(fā)電機(jī)有功出力變化的不同。

圖5 27號(hào)發(fā)電機(jī)有功出力變化曲線(xiàn)

由圖5可以看出：計(jì)及爬坡約束時(shí)，為了可以順利到達(dá)12點(diǎn)時(shí)的負(fù)荷高峰，27號(hào)發(fā)電機(jī)組在第5時(shí)段作出了相應(yīng)的調(diào)整，第5時(shí)段系統(tǒng)的負(fù)荷水平降低，但是此時(shí)煤耗較小的27號(hào)發(fā)電機(jī)就開(kāi)始增加出力，從第6時(shí)段開(kāi)始則全力爬坡，以保證在12點(diǎn)的高峰期能夠承擔(dān)更多的負(fù)荷；同樣，為了應(yīng)付22-24時(shí)段負(fù)荷的陡降，即使是在21-22時(shí)段負(fù)荷上升的情況下，27號(hào)發(fā)電機(jī)仍然減少了出力，并在23、24時(shí)段全速減小出力，保證系統(tǒng)能夠順利到達(dá)負(fù)荷的低點(diǎn)。

而如果未考慮爬坡約束，在12點(diǎn)負(fù)荷高峰前的多個(gè)爬坡時(shí)段和負(fù)荷陡降的22-24時(shí)段，27號(hào)發(fā)電機(jī)有功出力調(diào)節(jié)速率均超出了允許范圍。由此可見(jiàn)，靜態(tài)優(yōu)化調(diào)度的結(jié)果顯然是不可行的，只有動(dòng)態(tài)優(yōu)化調(diào)度的結(jié)果才是正確有效的。

3.3 雙目標(biāo)優(yōu)化與單目標(biāo)優(yōu)化結(jié)果對(duì)比分析

圖6顯示了雙目標(biāo)優(yōu)化時(shí)，各時(shí)段的煤耗量與污染物氣體排放量都比單目標(biāo)時(shí)的優(yōu)化結(jié)果大。這是因?yàn)榭紤]多目標(biāo)時(shí)，不僅要使得系統(tǒng)煤耗量最小，同時(shí)還要使污染物氣體排放量最小，但這2個(gè)目標(biāo)是相互矛盾的，污染物氣體排放滿(mǎn)意度的增加必然要以煤耗滿(mǎn)意度的下降來(lái)?yè)Q取，反之亦然，因此只能使得兩者的綜合滿(mǎn)意度達(dá)到最高。

表3比較了雙目標(biāo)與單目標(biāo)時(shí)整個(gè)調(diào)度周期內(nèi)的總煤耗量和污染物氣體總排放量。

圖6 雙目標(biāo)與單目標(biāo)各時(shí)段優(yōu)化結(jié)果對(duì)比

表3 雙目標(biāo)優(yōu)化和單目標(biāo)優(yōu)化目標(biāo)值比較

從表3可以看到，雙目標(biāo)優(yōu)化之后，煤耗量比單獨(dú)優(yōu)化煤耗時(shí)增加了256元，相對(duì)增長(zhǎng)2.44%，但比單獨(dú)優(yōu)化排放時(shí)的煤耗量減少了498元，相對(duì)減少4.42%；同樣，雙目標(biāo)優(yōu)化之后，污染物氣體排放量比單獨(dú)優(yōu)化排放時(shí)增加了89.4 t，也僅相對(duì)增長(zhǎng)2.44%，但比單獨(dú)優(yōu)化煤耗時(shí)的排放量減小了461.3 t，相對(duì)減小10.95%。

由此可見(jiàn)，雖然調(diào)度周期內(nèi)的總煤耗量和排放量都沒(méi)有單目標(biāo)模型時(shí)的理想，但從綜合效益的角度看，基于節(jié)能減排的多目標(biāo)動(dòng)態(tài)優(yōu)化調(diào)度模型具有明顯的優(yōu)勢(shì)，它能夠更加合理的協(xié)調(diào)各個(gè)目標(biāo)之間的關(guān)系。

4 結(jié)束語(yǔ)

本文以節(jié)能、減排作為動(dòng)態(tài)優(yōu)化調(diào)度的雙重目標(biāo)，建立了多目標(biāo)動(dòng)態(tài)優(yōu)化調(diào)度模型，并根據(jù)模糊集理論的最大隸屬度原則使之轉(zhuǎn)化為求解滿(mǎn)意度最大的單目標(biāo)問(wèn)題，提出了適合求解多目標(biāo)動(dòng)態(tài)優(yōu)化問(wèn)題的原對(duì)偶解耦內(nèi)點(diǎn)算法。IEEE-30節(jié)點(diǎn)算例結(jié)果表明：

（1）本文通過(guò)對(duì)原對(duì)偶內(nèi)點(diǎn)法中修正方程的降階與解耦，大幅度提高了算法的求解效率，收斂特性較好，數(shù)值魯棒性高。

（2）本文提出的基于節(jié)能減排的多目標(biāo)動(dòng)態(tài)優(yōu)化調(diào)度模型能夠更好的協(xié)調(diào)節(jié)約能源、減少環(huán)境污染之間的關(guān)系，并能根據(jù)機(jī)組的煤耗特性、排放特性以及爬坡約束等更加合理的安排發(fā)電機(jī)組的出力，在整個(gè)調(diào)度周期內(nèi)實(shí)現(xiàn)各個(gè)目標(biāo)函數(shù)綜合滿(mǎn)意度的最大化。

（3）本文所提的模型與算法不僅適合于雙目標(biāo)動(dòng)態(tài)優(yōu)化調(diào)度問(wèn)題，還能根據(jù)實(shí)際情況靈活擴(kuò)展到更多的目標(biāo)函數(shù)。

［1］ 覃智君，陽(yáng)育德，吳杰康.矢量化動(dòng)態(tài)最優(yōu)潮流計(jì)算的步長(zhǎng)控制內(nèi)點(diǎn)法實(shí)現(xiàn)[J].中國(guó)電機(jī)工程學(xué)報(bào)，2009，29(7)：52-58.

［2］ 張伯明，陳壽孫，嚴(yán)正.高等電力網(wǎng)絡(luò)分析（第2版）[M].北京：清華大學(xué)出版社，2007.

［3］ IMMANUEL SELVAKUMAR A，THANUSHKODI K.A New Particle Swarm Optimization Solution to Nonconvex Economic Dispatch Problems[J].IEEE Transactions on Power Systems,2007,22(1)：42-51.

［4］ GRANVILLE S.Optimal Reactive Dispatch Through Interior Point Method[J].IEEE Transactions on Power Systems,1994,9(1)：136-146.

［5］ 李秀卿,羅忠游,郭芳義.基于網(wǎng)絡(luò)流和內(nèi)點(diǎn)法的電力系統(tǒng)動(dòng)態(tài)經(jīng)濟(jì)調(diào)度[J].繼電器，2006，34(10)：37-40.

［6］ 賴(lài)永生,劉明波.電力系統(tǒng)動(dòng)態(tài)無(wú)功優(yōu)化問(wèn)題的快速解耦算法[J].中國(guó)電機(jī)工程學(xué)報(bào)，2008，28(7)：32-39.

［7］ 馮士剛，艾芊.火電站多目標(biāo)負(fù)荷調(diào)度及其算法的研究[J].動(dòng)力工程，2008，28(3)：404-407.

［8］ 劉明波,段曉軍,趙艷.多目標(biāo)最優(yōu)潮流問(wèn)題的模糊建模及內(nèi)點(diǎn)解法[J].電力系統(tǒng)自動(dòng)化，1999，23(14)：37-40.

［9］ 艾琳，華棟.基于模糊優(yōu)化和內(nèi)點(diǎn)法的水火電系統(tǒng)短期經(jīng)濟(jì)調(diào)度[J].電力自動(dòng)化設(shè)備，2008，28(8)：46-50.

［10］ 馮培燕.基于模糊多目標(biāo)優(yōu)化的火電廠(chǎng)負(fù)荷優(yōu)化分配研究[D].北京:華北電力大學(xué)，2007.

［11］ 馬瑞，穆大慶，李欣然.電力市場(chǎng)中日有功負(fù)荷多目標(biāo)分配模糊決策的研究[J].電網(wǎng)技術(shù)，2001，25(2)：25-29.

［12］ XIE K,SONG Y H.Optimal Power Flow with Time-related Constrains by a Nonlinear Interior Point Method[A].Proceedings of 2000 IEEE Power Engineering Society Winter Meeting[C].Singapore，2000.

［13］ 王錫凡，方萬(wàn)良，杜正春.現(xiàn)代電力系統(tǒng)分析[M].北京：科學(xué)出版社，2003.

［14］ XIE K，SONG Y H.Dynamic Optimal Power Flow by Interior Point Methods[J].IEE Proceedings：Generation,Transmission and Distribution，2001，148(1)：76-84.