[發(fā)明專利]一種在極坐標(biāo)系下減小區(qū)間潮流保守性的方法

有效

• 申請?zhí)枺?/td>: 201410069397.X
• 申請日: 2014-02-27
• 公開（公告）號: CN103779863A
• 公開（公告）日: 2014-05-07
• 發(fā)明（設(shè)計(jì)）人: 張衛(wèi)東;劉靜偉;何星;李智軍;蔡云澤
• 申請（專利權(quán)）人: 上海交通大學(xué)
• 主分類號: H02J3/00
• 分類號: H02J3/00;G06F19/00
• 代理公司: 上海科盛知識產(chǎn)權(quán)代理有限公司 31225
• 代理人: 趙志遠(yuǎn)

1.一種在極坐標(biāo)系下減小區(qū)間潮流保守性的方法，其特征在于，包括以下步驟：

(1)獲取網(wǎng)絡(luò)元件、線路參量的區(qū)間值并用仿射區(qū)間量表示出來，獲取網(wǎng)絡(luò)拓?fù)浣Y(jié)構(gòu)；

(2)在網(wǎng)絡(luò)拓?fù)浣Y(jié)構(gòu)的基礎(chǔ)上建立導(dǎo)納矩陣；

(3)形成雅可比矩陣，并用三角函數(shù)法計(jì)算仿射區(qū)間量；

(4)通過迭代得到減小了保守性的潮流結(jié)果。

2.根據(jù)權(quán)利要求1所述的一種在極坐標(biāo)系下減小區(qū)間潮流保守性的方法，其特征在于，該方法具體過程如下：

步驟一：根據(jù)電網(wǎng)拓?fù)浣Y(jié)構(gòu)和元件參數(shù)求出導(dǎo)納矩陣Y：

Y=[G‾,G‾]+i[B‾,B‾]]]>

其中G是電阻，B是電納，和*分別表示區(qū)間數(shù)的上限和下限，將導(dǎo)納矩陣Y轉(zhuǎn)化為放射區(qū)間的形式為：

IF(Y)=(G‾+G‾)/2+(G‾-G‾)/2+ξ[(B‾+B‾)/2+(B‾-B‾)/2]]]>

其中ξ是噪聲元

步驟二：設(shè)置自變量區(qū)間初始值仿射形式為，：

IF(X(k))=(U‾(K)+U‾(k))/2(θ‾(k)+θ‾(k))/2+(U‾(k)-U‾(k))/2(θ‾(k)+θ‾(k))/2ξ,]]>

其中k=0，U^(k)和θ^(k)是第K次迭代的電壓和相角；

步驟三：求區(qū)間雅可比矩陣

F′(J(k))=[H‾(k),H‾(k)][N‾(k),N‾(k)][K‾(k),K‾(k)][L‾(k),L‾(k)]]]>

其中H^(k)，N^(k)，K^(k)，L^(k)分別是第k次迭代雅可比矩陣的元素矩陣，由于矩陣元素之間的相關(guān)性，此處結(jié)合三角函數(shù)的特性，減小由忽略元素之間的相關(guān)性引起的保守性，以矩陣H為例：

Hij=∂ΔPi∂θj=-UiUj(Gijsinθij-Bijcosθij)]]>

Hij=∂ΔPi∂θ=UiΣj=1,i≠jnUj(Gijsinθij-Bijcosθij)]]>

其中，H_ij是矩陣H中的第i行第j列元素，U是電壓幅值，θ是電壓相角，采用三角函數(shù)分析方法，將故雅可比矩陣中元素H可按照以下方法進(jìn)行計(jì)算：

Hij=∂ΔPi∂θj=-UiUjGij2+Bij2sin(θij+δij)]]>

Hij=∂ΔPi∂θ=-UiΣj=1,i≠jnUjGij2+Bij2sin(θij+δij)]]>

相應(yīng)的N，K，L也可以用這種方法計(jì)算得到，分別如下：

Nij=∂ΔPi∂UθjUj=-UiUjGij2+Bij2cos(θij-δij)]]>

Nij=∂ΔPi∂UjUi=-U12Gij-Pi]]>

Kij=∂ΔQi∂θj=UiUjGij2+Bij2cos(θij-δij)]]>

Kij=∂ΔQi∂θj=Ui2Gii-Pi]]>

Lij=∂ΔQi∂UjUj=-UiUjGij2+Bij2sin(θij-δij)]]>

Lij=∂ΔPi∂UjUi=Ui2Bij-Qi]]>

將區(qū)間值J^(k)的區(qū)間中值表示為J^(k)′：

J(k)′=m(J(k))-1]]>

步驟四：用仿射運(yùn)算計(jì)算得到AF(J^(k)′F′(J^(k)))；

步驟五：得到仿射形式AF(I-J^(k)′F′(J^(k)))和AF(X^(k)-m(X^(k)))＝(rξ)^(k)，求得AF((I-Y^(k)F′(X^(k)))(rξ)^(k))，并轉(zhuǎn)化成區(qū)間形式，其中I是區(qū)間數(shù)的單位，ξ是引入仿射區(qū)間后噪聲元單位，r是噪聲元的系數(shù)，代表噪聲元的大小，也表示不同噪聲元之間的聯(lián)系的強(qiáng)弱，X是要求解的未知區(qū)間數(shù)。

步驟六：按迭代算子：

K(X^(k))=m(X^(k))-Y^(k)f(m(X^(k)))+(I-Y^(k)F′(X^(k)))(X^(k)-m(X^(k)))，計(jì)算出K(X^(k))得到區(qū)間形式的X^(k+1)=X^(k)∩K(X^(k))；

步驟七：判斷是否滿足終止條件，即

步驟八：輸出結(jié)果。