[發(fā)明專利]一種基于支路電阻偏差的牛頓-拉夫遜法直流電網(wǎng)潮流計算方法

有效

• 申請?zhí)枺?/td>: 201811306817.6
• 申請日: 2018-11-05
• 公開（公告）號: CN109586268B
• 公開（公告）日: 2022-02-11
• 發(fā)明（設(shè)計）人: 王淳;陳宇杰;劉寬;伍惠鋮;王濤;李成豪
• 申請（專利權(quán)）人: 南昌大學(xué)
• 主分類號: H02J1/00
• 分類號: H02J1/00
• 代理公司: 南昌新天下專利商標(biāo)代理有限公司 36115
• 代理人: 施秀瑾

1.一種基于支路電阻偏差的牛頓-拉夫遜法直流電網(wǎng)潮流計算方法，其特征是包括以下步驟：

步驟1、讀取網(wǎng)絡(luò)數(shù)據(jù)，包括節(jié)點和支路數(shù)據(jù)，設(shè)置節(jié)點電壓和回路電流初值，迭代計數(shù)k為0；

步驟2、根據(jù)支路的首末段節(jié)點編號形成支路-節(jié)點關(guān)聯(lián)矩陣A，并根據(jù)式(2)和式(3)得到支路-道路關(guān)聯(lián)矩陣T、支路-回路關(guān)聯(lián)矩陣B；

式中，A_t為樹枝的支路-節(jié)點關(guān)聯(lián)子矩陣，A_l為連枝的支路-節(jié)點關(guān)聯(lián)子矩陣；T_t為樹枝的支路-道路關(guān)聯(lián)子矩陣，T_l為連枝的支路-道路關(guān)聯(lián)子矩陣；B_t為樹枝的支路-回路關(guān)聯(lián)子矩陣，B_l為連枝的支路-回路關(guān)聯(lián)子矩陣；0表示零矩陣；E表示單位矩陣；

步驟3、根據(jù)式(4)、式(5)形成雅克比矩陣；

H′＝A (4)

G′＝diag(R)B (5)

式中，diag()表示將列向量轉(zhuǎn)化成對角矩陣；R表示給定的支路電阻列向量；

步驟4、根據(jù)式(6)、式(7)和式(8)分別求出U_b^(k)、I_b^(k)和ΔR^(k)；

U_b＝A(U-eU_s) (6)

I_b＝TI_n+BI_l＝T(P./U)+BI_l (7)

式中，U_b為支路電壓列向量，U為節(jié)點電壓列向量，U_s為參考節(jié)點電壓；e為單位列向量；I_n為節(jié)點注入電流列向量、P為節(jié)點注入功率列向量、I_b為支路電流列向量，I_l為回路電流列向量；R_i表示第i條支路的支路電阻給定值，R_ci表示第i條支路的支路電阻計算值，ΔR_i表示支路電阻不平衡量；A_i、T_i、B_i分別表示矩陣A、T、B的第i行分量；./表示兩矩陣對應(yīng)元素相除；

步驟5、根據(jù)式(9)求出節(jié)點電壓修正量ΔU^(k)和回路電流修正量ΔI_l^(k)，并按照式(10)對節(jié)點電壓和回路電流進(jìn)行修正：

式中，J'表示修正后的雅克比矩陣，.*表示兩矩陣對應(yīng)元素相乘；

步驟6、按照式(11)判斷支路阻抗的不平衡量ΔR^(k)是否達(dá)到要求的收斂精度ε，如果式(11)成立，迭代到此結(jié)束；反之，繼續(xù)進(jìn)行迭代：

max{|ΔR^(k)|}ε (11)

步驟7、迭代計數(shù)k加1，返回步驟4繼續(xù)迭代過程。