【正文】 了 . 17 3. PQ 分解法計(jì)算潮流的程序框圖 開(kāi) 始 輸入原始數(shù)據(jù) 形成矩陣 B ′及 B ″ 形成 B ′， B ″因子表 設(shè) P Q 節(jié)點(diǎn)電壓初值，各節(jié)點(diǎn)電壓相角初值 置迭代次數(shù) K = 0 計(jì)算 Δ Pi/ Vi ( N 1 ) 個(gè) Δ Qi/ Vi ( m ) 個(gè) 迭代收斂了嗎？ 利用 B ′因子表對(duì)右端項(xiàng)變換 回代求解 修正電壓角 利用 B ″對(duì)右端項(xiàng)變換 回代求解 修正電壓模 k = k + 1 計(jì)算支路功率 計(jì)算支路功率損耗 輸出結(jié)果 結(jié)束 18 PQ 分解法源程序 include include include main () { int a=1,k,i,j,d=1,m,n,ls,nl,kk。 DD VBVH ??? ⑥22 DD VBVL ??? 將⑤⑥分別帶入③④ 得： ??????11 39。 2. H、 L 元素都是節(jié)點(diǎn)電壓幅值和相角差函數(shù)，其數(shù)值在迭代過(guò)程中是不斷變化的，從而，最關(guān)鍵的一步是把系數(shù)矩陣 H、 L 簡(jiǎn)化為常數(shù)矩陣。 對(duì) 于 每 一 個(gè) PQ 節(jié) 點(diǎn) 或 每 一 個(gè) PV 節(jié) 點(diǎn) 都 可 以 列 寫(xiě) 一 個(gè) 有 功 功 率 不 平 衡 方 程 ：) 1n , 2 , 1i ( 0)s i nc os(1????????? ??nj ijijijijjiisiisiBGVVPPPP ?? 對(duì)于每一個(gè) PQ節(jié)點(diǎn)可以列寫(xiě)一個(gè)無(wú)功功率不平衡方程式： ) m , 2 , 1i ( 0)c oss i n(1????????? ??nj ijijijijjiisiisiBGVV ?? 所以可以寫(xiě)出修正方程式 ②????????????????????????? V V 1D 2?LKNHQP 其中 ② ??????????????????121 nPPPP? ????????????????? 21m? ??????????????????121 n????? ?????????????????mVVVV 21? ?????????????3212 VVVV D? 16 四、 PQ 分解法 1．在交流高壓電網(wǎng)中，輸電線路的電抗要比電阻大得多，系統(tǒng)中母線有功功率的變化主要受電壓相位的影響，無(wú)功功率的變化規(guī)則主要受母線電壓幅值變化的影響，在修正方程式的系數(shù)矩陣中，偏導(dǎo)數(shù)VP????和????P的數(shù)值相當(dāng)于偏導(dǎo)數(shù)????? P和VQ????是相當(dāng)小的，作為簡(jiǎn)化第一步，可將方程②中 N、 k 略去不計(jì)，既認(rèn)為它們的元素都等于零，這樣，便可分解為 n1 階和 m 階的兩個(gè)方程： ③ ?????? HP ④ VLVQ D ????? ?1 所謂 PQ分解法，節(jié)點(diǎn)有功功率不平衡量又用于修正電壓的相位。 PV節(jié)點(diǎn)的電壓幅值 Vm+1~Vn1也是給定的。方程組②的求解屬于回代過(guò)程，只順用到上三角因子矩陣 以及經(jīng)過(guò)消元變換的右端常數(shù)向量，方程組可以自下而上地逐步算出待求量，其計(jì)算通式為： 三、節(jié)點(diǎn)電壓的表示： 1. 極坐標(biāo)表示法： )( s i n iiiii VVV ?? ??? 節(jié)點(diǎn)功率方程表示為： ???? ni ijijijijjiiBGVVP1)s i nc os( ?? ① ???? nj ijjijijjiiBiGVVQ1)c oss i n( ?? jiij ??? ?? （兩節(jié)點(diǎn)電壓的相位角） 方程式①把節(jié)點(diǎn)功率表示為節(jié)點(diǎn)電壓的幅值和相角的函數(shù) ) 1 , 1n ,n ( )()1(1)1()1(???? ????? ?irxjrfx iijnijiijiii 15 在有 n 個(gè)節(jié)點(diǎn)的系統(tǒng)中，假定第 1 ~ m 號(hào)節(jié)點(diǎn)為 PQ節(jié)點(diǎn)，第 m+1~n1 號(hào)節(jié)點(diǎn)為平衡節(jié)點(diǎn)。 （ 4） 對(duì)于變壓器支路，采用下圖所示等值電路，即非標(biāo)準(zhǔn)變比在 J 側(cè)，變壓器阻抗在 I 側(cè)， K 處填寫(xiě)實(shí)際變化 2． 框圖 12 開(kāi) 始 導(dǎo) 納 矩 陣 清 零 D0 k = 1 , L 取支路 K 數(shù)據(jù) I 、 J 、 R 、X 、 K （ B / 2 ） 是接地支 路嗎？ 是變壓器支路嗎？ 線路容納追加到兩節(jié)點(diǎn)自導(dǎo)納中 P = 1 ， Q = 1 支路導(dǎo)納追加到 I 節(jié)點(diǎn)的自導(dǎo)納中 I 、 J 之間互導(dǎo)納中 將支路導(dǎo)納追加到 J 節(jié)點(diǎn)的自導(dǎo)納中 C O N T I N U E 輸 出 結(jié) 果 結(jié) 束 P = 1 P = k2 、 Q = k 注： B1 —— 不記充電電容和非標(biāo)準(zhǔn)變比時(shí)，導(dǎo)納矩陣的虛部 B2 —— 導(dǎo)納矩陣的虛部 二、形成因子表的方法及框圖 13 1．用行消去過(guò)程形成因子表的程序框圖 開(kāi) 始 輸入矩陣階數(shù) N ，定義數(shù)組 A （ N . N ) 讀入矩陣元素 D 0 I = 2 , N D 0 k = 1 , I 1 A ( I , k ) = A ( I , k ) / A ( k , k ) D 0 J = k + 1 , N A ( I , J ) = A ( I , J ) + A ( I , k ) * A ( k , J ) C O N T I N U E C O N T I N U E C O N T I N U E 輸 出 結(jié) 果 結(jié)束 2. 因子表的形成（三角分解法） 消去法求非線性方程組的一種常用算法是對(duì)方程式的系數(shù)矩陣 A進(jìn)行三角分解，在本次電力系統(tǒng)潮流計(jì)算中采用的三角分解是將非奇方陣 A分解為單位下三角矩陣 L 和上三角矩陣 R 的乘積。 （ 2） 對(duì)于純阻抗支路， K 處置零。 一、 形成導(dǎo)納矩陣的方法及框圖 11 1． 形成導(dǎo)納矩陣的程序框圖 為了形成導(dǎo)納矩陣，必須個(gè)計(jì)算機(jī)輸入必要的原始數(shù)據(jù)，一條支路 的原始數(shù)據(jù)應(yīng)包括兩端的節(jié)點(diǎn)號(hào)和支路阻抗，對(duì)變壓器支路的原始數(shù)據(jù)應(yīng)包括他的變比，對(duì)于電力線路還應(yīng)包括它的容納，即： I ， J ， R ， X ， K（ B/2）。 3． 平衡節(jié)點(diǎn)：網(wǎng)絡(luò)中至少有一個(gè)節(jié)點(diǎn)的有功功率 P 不能給定，此節(jié)點(diǎn)承擔(dān)系統(tǒng)的有功功率平衡。 標(biāo)么值 =實(shí)際有名值（任意單位） /基準(zhǔn)值（與有名值同單位） 對(duì)于直接電氣聯(lián)系的網(wǎng)絡(luò) ,在制訂標(biāo)么值的等值電路時(shí) ,各元件的參數(shù)必須按統(tǒng)一的基準(zhǔn)值進(jìn)行歸算 .由于各元件的額定值可能不同 ,因此 ,必須把不同基準(zhǔn)值的標(biāo)么阻抗換算成統(tǒng)一基準(zhǔn)值的標(biāo)么值 .現(xiàn)統(tǒng)一選定的基準(zhǔn)電壓和基準(zhǔn)功率分別為 V和 S,則電抗的實(shí)際有名值換算成標(biāo)么值 ,即 在工程計(jì)算中規(guī)定 ,各個(gè)電壓等級(jí)都以其平均額定電壓 V 作為基準(zhǔn)電壓 .根據(jù)我國(guó)現(xiàn)行的電壓等級(jí) ,各級(jí)平均額定電壓規(guī)定為 , , , , 37, 115, 230, 345, 525KV 五 .具體計(jì)算及結(jié)果 7 根據(jù)原始數(shù)據(jù)及資料顯示 ,通過(guò)公式計(jì)算可得 :(選擇基準(zhǔn)值 :SB=100MVA,VB=Vav) 雙繞組變壓器 :B1,B2,B3（型號(hào) SFP7150000/220 ） （參數(shù) ΔP0=140KW,ΔPS=450KW,VS%=14%,I0%=,YN—d11) RT=ΔPS UN2 103/SN2=450 2202 103/1500002= RT*=RTS B/VB2= 100/2302= XT=VX%VN2103/100/SN=142202103/100/150000= XT*=XTS B/VB2=100/2302= BT=I0%SN103/100/VN2=150000103/100/2302=103S BT*=BTV B2/SB=1032302/100=