[發(fā)明專利]一種電力系統(tǒng)潮流計算的導(dǎo)納矩陣計算方法

有效

• 申請?zhí)枺?/td>: 201710942153.1
• 申請日: 2017-10-11
• 公開（公告）號: CN109659943B
• 公開（公告）日: 2022-03-29
• 發(fā)明（設(shè)計）人: 姚玉斌;趙偉;吳志良
• 申請（專利權(quán)）人: 大連海事大學(xué)
• 主分類號: H02J3/06
• 分類號: H02J3/06;G06F17/16
• 代理公司: 大連東方專利代理有限責(zé)任公司 21212
• 代理人: 李洪福

1.一種電力系統(tǒng)潮流計算的導(dǎo)納矩陣計算方法，其特征在于：包括以下步驟：

B1、讀支路首節(jié)點號數(shù)組I、末節(jié)點號數(shù)組J、電阻數(shù)組R、電抗數(shù)組X、對地電納數(shù)組B、變壓器變比數(shù)組K；

所述的首節(jié)點號數(shù)組I、末節(jié)點號數(shù)組J、電阻數(shù)組R、電抗數(shù)組X、對地電納數(shù)組B、變壓器變比數(shù)組K分別按順序存放所有支路的首節(jié)點號i_m、末節(jié)點號j_m、電阻r_m、電抗x_m、對地電納b_m、變壓器變比k_m，其中變壓器支路對地電納b_m為0，變壓器非標準變比k_m側(cè)的節(jié)點號加個負號，輸電線路變比k_m設(shè)為1，其中下標m為支路序號；

B2、定義AS和AE為n×l階稀疏矩陣，矩陣AS為有向支路首節(jié)點與支路的關(guān)聯(lián)矩陣，矩陣AE為有向支路末節(jié)點與支路的關(guān)聯(lián)矩陣，n為節(jié)點數(shù)，l為支路數(shù)；；

B3、設(shè)置支路計數(shù)m＝1；

B4、判斷支路m的首節(jié)點號I_m是否小于0，如果不滿足轉(zhuǎn)至步驟B6；

B5、令i＝|I_m|、j＝|J_m|，轉(zhuǎn)至步驟B7；

B6、令i＝|J_m|、j＝|I_m|；

B7、令A(yù)S_im＝1、AE_jm＝1；

B8、令m＝m+1；

B9、判斷m是否大于支路數(shù)l，如果m不大于l，則返回到步驟B4；

B10、計算支路π形等值電路串聯(lián)支路導(dǎo)納數(shù)組；

支路π形等值電路的串聯(lián)支路導(dǎo)納數(shù)組Y_B為：

Y_B＝1·/(R+jX)·/K (1)

式中，“·/”表示兩數(shù)組對應(yīng)元素相除；

B11、計算節(jié)點支路關(guān)聯(lián)矩陣；

節(jié)點支路關(guān)聯(lián)矩陣A為：

A＝AS-AE (2)

B12、由串聯(lián)支路形成初始導(dǎo)納矩陣；

由支路π形等值電路串聯(lián)支路形成的初始導(dǎo)納矩陣Y1為：

Y1＝AY_BDA^T (3)

式中，Y_BD為串聯(lián)支路導(dǎo)納數(shù)組Y_B形成的對角矩陣，上標T表示矩陣的轉(zhuǎn)置；

B13、計算有向支路首節(jié)點和末節(jié)點的對地導(dǎo)納數(shù)組；

有向支路首節(jié)點和末節(jié)點的對地導(dǎo)納數(shù)組YS和YE分別為：

YS＝Y(jié)_B·*(1-K)·/K+jB/2 (4)

YE＝Y(jié)_B·*(K-1)+jB/2 (5)

式中，“·*”表示兩數(shù)組對應(yīng)元素相乘；

B14、形成各節(jié)點對地導(dǎo)納數(shù)組；

各節(jié)點對地導(dǎo)納數(shù)組Y0為：

Y0＝AS·YS+AE·YE (6)

B15、設(shè)置節(jié)點計數(shù)m＝1；

B16、令Y1_mm＝Y(jié)1_mm+Y0_m；

B17、令m＝m+1；

B18、判斷m是否大于節(jié)點數(shù)n，如果m不大于n，則返回到步驟B16；否則，結(jié)束。