電力系統(tǒng)的潮流計算

電力系統(tǒng)的潮流計算第1頁，課件共153頁，創(chuàng)作于2023年2月主要內(nèi)容:基本概念數(shù)學(xué)模型-潮流方程潮流方程的求解相關(guān)技術(shù)電力系統(tǒng)潮流計算程序設(shè)計電力系統(tǒng)狀態(tài)估計電力系統(tǒng)最優(yōu)潮流2第2頁，課件共153頁，創(chuàng)作于2023年2月一.基本概念電力系統(tǒng)潮流計算是對復(fù)雜電力系統(tǒng)正常和故障條件下穩(wěn)態(tài)運行狀態(tài)的計算。其目的是求取電力系統(tǒng)在給定運行方式下的節(jié)點電壓和功率分布，用以檢查系統(tǒng)各元件是否過負荷、各點電壓是否滿足要求、功率分布和分配是否合理以及功率損耗等。潮流計算的計算機算法是以電網(wǎng)絡(luò)理論為基礎(chǔ)的，應(yīng)用數(shù)值計算方法求解一組描述電力系統(tǒng)穩(wěn)態(tài)特性的方程。3第3頁，課件共153頁，創(chuàng)作于2023年2月一.基本概念

潮流計算方法的要求：計算速度快內(nèi)存需要小計算結(jié)果有良好的可靠性和可信性適應(yīng)性好，即能處理變壓器變比調(diào)整、系統(tǒng)元件的不同描述和與其他程序配合的能力強簡單4第4頁，課件共153頁，創(chuàng)作于2023年2月一.基本概念潮流計算方法的步驟：建立潮流的數(shù)學(xué)模型確定適宜的計算方法制定計算流程圖編制計算機程序?qū)τ嬎憬Y(jié)果進行分析和確定，檢查程序的正確性5第5頁，課件共153頁，創(chuàng)作于2023年2月二.數(shù)學(xué)模型-潮流方程節(jié)點導(dǎo)納矩陣節(jié)點阻抗矩陣功率方程節(jié)點分類6第6頁，課件共153頁，創(chuàng)作于2023年2月(一).節(jié)點導(dǎo)納矩陣y12y23y13y20y10y30I20I107第7頁，課件共153頁，創(chuàng)作于2023年2月(一).節(jié)點導(dǎo)納矩陣運用基爾霍夫電流定律可以得到:8第8頁，課件共153頁，創(chuàng)作于2023年2月(一).節(jié)點導(dǎo)納矩陣整理:9第9頁，課件共153頁，創(chuàng)作于2023年2月(一).節(jié)點導(dǎo)納矩陣整理:

節(jié)點電壓法10第10頁，課件共153頁，創(chuàng)作于2023年2月(一).節(jié)點導(dǎo)納矩陣節(jié)點電壓方程

IB：為節(jié)點注入電流的列向量，可理解為各節(jié)點電源電流與負荷電流之和，并規(guī)定電源流向網(wǎng)絡(luò)的注入電流為正；

UB：為節(jié)點電壓的列向量；

YB：為節(jié)點導(dǎo)納矩陣。11第11頁，課件共153頁，創(chuàng)作于2023年2月(一).節(jié)點導(dǎo)納矩陣互導(dǎo)納自導(dǎo)納12第12頁，課件共153頁，創(chuàng)作于2023年2月(一).節(jié)點導(dǎo)納矩陣節(jié)點導(dǎo)納矩陣元素的物理意義對角元Yii稱為自導(dǎo)納：節(jié)點i注入電流與該節(jié)點i電壓之比，數(shù)值上等于該節(jié)點直接連接的所有支路導(dǎo)納的總和；非對角元Yij稱為互導(dǎo)納：節(jié)點i注入電流與該節(jié)點j電壓之比，數(shù)值上等于連接節(jié)點i，j支路導(dǎo)納的負值。13第13頁，課件共153頁，創(chuàng)作于2023年2月(一).節(jié)點導(dǎo)納矩陣節(jié)點導(dǎo)納矩陣的特點N×N階方陣對稱復(fù)數(shù)矩陣高度稀疏矩陣14第14頁，課件共153頁，創(chuàng)作于2023年2月(一).節(jié)點導(dǎo)納矩陣節(jié)點導(dǎo)納矩陣的修改原網(wǎng)絡(luò)節(jié)點增加一接地支路(設(shè)在節(jié)點i增加一接地支路X)

15第15頁，課件共153頁，創(chuàng)作于2023年2月節(jié)點導(dǎo)納矩陣的修改原網(wǎng)絡(luò)節(jié)點i，j增加一條支路X16第16頁，課件共153頁，創(chuàng)作于2023年2月節(jié)點導(dǎo)納矩陣的修改設(shè)原網(wǎng)絡(luò)有n個節(jié)點，從節(jié)點i(i≤n)引出一條支路X及新增一節(jié)點j。17第17頁，課件共153頁，創(chuàng)作于2023年2月節(jié)點導(dǎo)納矩陣的修改刪除網(wǎng)絡(luò)中的一條支路

與增加相反，可理解為增加了一條負支路修改原網(wǎng)絡(luò)中的支路參數(shù)

可理解為先將被修改支路刪除，然后增加一條參數(shù)為修改后導(dǎo)納值的支路。因此，修改原網(wǎng)絡(luò)中的支路參數(shù)可通過給原網(wǎng)絡(luò)并聯(lián)一條支路來實現(xiàn)。18第18頁，課件共153頁，創(chuàng)作于2023年2月節(jié)點導(dǎo)納矩陣的修改增加一臺變壓器(阻抗+理想變壓器)

19第19頁，課件共153頁，創(chuàng)作于2023年2月節(jié)點導(dǎo)納矩陣的修改增加一臺變壓器(理想變壓器+阻抗)

20第20頁，課件共153頁，創(chuàng)作于2023年2月節(jié)點導(dǎo)納矩陣的修改將節(jié)點i、j之間變壓器的變比由k改為k’21第21頁，課件共153頁，創(chuàng)作于2023年2月節(jié)點導(dǎo)納矩陣的形成根據(jù)節(jié)點導(dǎo)納矩陣的特點形成自導(dǎo)納Yii為所有連接節(jié)點i的支路的導(dǎo)納之和?；?dǎo)納Yij為連接節(jié)點i，j支路導(dǎo)納的負值。根據(jù)節(jié)點導(dǎo)納矩陣的修改形成無支路支路添加22第22頁，課件共153頁，創(chuàng)作于2023年2月節(jié)點導(dǎo)納矩陣的形成例題23第23頁，課件共153頁，創(chuàng)作于2023年2月例題24第24頁，課件共153頁，創(chuàng)作于2023年2月參數(shù)化為導(dǎo)納值25第25頁，課件共153頁，創(chuàng)作于2023年2月所求的節(jié)點導(dǎo)納矩陣26第26頁，課件共153頁，創(chuàng)作于2023年2月（二）.節(jié)點阻抗矩陣以地為參考節(jié)點的節(jié)點導(dǎo)納矩陣Y是N×N階稀疏矩陣；如果網(wǎng)絡(luò)中存在接地支路，Y是非奇異的，其逆矩陣是節(jié)點阻抗矩陣：用節(jié)點阻抗矩陣Z表示的網(wǎng)絡(luò)方程是：

27第27頁，課件共153頁，創(chuàng)作于2023年2月節(jié)點阻抗矩陣的特點是對稱矩陣。對于連通的電力系統(tǒng)網(wǎng)絡(luò)，當網(wǎng)絡(luò)中有接地支路時，Z是非奇異滿矩陣。對純電阻性或電感性支路組成的電網(wǎng)，節(jié)點對的自阻抗不為零。

28第28頁，課件共153頁，創(chuàng)作于2023年2月節(jié)點阻抗矩陣的形成支路追加法

實質(zhì)上是與根據(jù)定義直接求節(jié)點導(dǎo)納矩陣的方法相對應(yīng)。根據(jù)自阻抗和互阻抗的定義直接節(jié)點阻抗矩陣的方法。29第29頁，課件共153頁，創(chuàng)作于2023年2月節(jié)點阻抗矩陣的形成利用節(jié)點導(dǎo)納矩陣逐列形成節(jié)點阻抗矩陣的方法

30第30頁，課件共153頁，創(chuàng)作于2023年2月作業(yè)1：已知一簡單系統(tǒng)的等值電路圖及元件參數(shù)，其中阻抗支路的參數(shù)以阻抗標注，導(dǎo)納支路的參數(shù)以導(dǎo)納標注（均為標么值），試求：1）該等值電路的節(jié)點導(dǎo)納矩陣；2）若支路34斷開，節(jié)點導(dǎo)納矩陣有何變化？3）若節(jié)點2接地，節(jié)點導(dǎo)納矩陣有何變化？4）若變壓器變比變?yōu)?：1.1，節(jié)點導(dǎo)納矩陣有何變化？31第31頁，課件共153頁，創(chuàng)作于2023年2月作業(yè)132第32頁，課件共153頁，創(chuàng)作于2023年2月（三）.功率方程在實際電力系統(tǒng)中，已知的運行條件往往不是節(jié)點的注入電流而是負荷和發(fā)電機的功率，而且這些功率一般不隨節(jié)點電壓的變化而變化，因此在節(jié)點功率不變的情況下，節(jié)點的注入電流隨節(jié)點電壓的變化而變化。在已知節(jié)點導(dǎo)納矩陣的情況下，必須用已知的節(jié)點功率來代替未知的節(jié)點注入電流，才能求出節(jié)點電壓。33第33頁，課件共153頁，創(chuàng)作于2023年2月（三）.功率方程34第34頁，課件共153頁，創(chuàng)作于2023年2月（三）.功率方程每節(jié)點的注入功率方程式為：極坐標直角坐標35第35頁，課件共153頁，創(chuàng)作于2023年2月（三）.功率方程功率方程36第36頁，課件共153頁，創(chuàng)作于2023年2月（三）.功率方程對于N個節(jié)點的電力網(wǎng)絡(luò)，可以列出2N個功率方程。每個節(jié)點具有四個變量，N個節(jié)點有4N個變量，但只有2N個關(guān)系方程式。37第37頁，課件共153頁，創(chuàng)作于2023年2月定解條件：在具有N個節(jié)點的系統(tǒng)中，給定(N-1)對控制變量PGi、QGi，余下一對控制變量待定PGs、QGs，其將使系統(tǒng)功率，包括電源功率、負荷功率和損耗功率保持平衡。給定一對狀態(tài)變量δs、Us，要求確定(n-1)對狀態(tài)變量δi、Ui，δs給定的通常為0，Us一般取標幺值為1，以使系統(tǒng)中各節(jié)點的電壓水平在額定值附近。除此之外，還應(yīng)滿足一些約束條件：U的約束條件：Umin<Ui<Umaxδ的約束條件：|δi-δj|<|δi-δj|max38第38頁，課件共153頁，創(chuàng)作于2023年2月（四）節(jié)點分類PQ節(jié)點注入有功和無功功率是給定的。相應(yīng)于實際電力系統(tǒng)中的一個負荷節(jié)點，或有功和無功功率給定的發(fā)電機母線。39第39頁，課件共153頁，創(chuàng)作于2023年2月（四）節(jié)點分類PV節(jié)點（電壓控制母線）注入有功功率Pi為給定值，電壓Ui也保持在給定數(shù)值。這種類型節(jié)點相當于發(fā)電機母線節(jié)點，其注入的有功功率由汽輪機調(diào)速器設(shè)定，而電壓則大小由裝在發(fā)電機上的勵磁調(diào)節(jié)器控制；或者相應(yīng)于一個裝有調(diào)相機或靜止補償器的變電所母線，其電壓由可調(diào)無功功率的控制器設(shè)定。要求有連續(xù)可調(diào)的無功設(shè)備，調(diào)無功來調(diào)電壓值。40第40頁，課件共153頁，創(chuàng)作于2023年2月（四）節(jié)點分類平衡節(jié)點平衡節(jié)點的電壓和相位大小是給定的，通常以它的相角為參考量，即取其電壓相角為0。這種節(jié)點用來平衡全電網(wǎng)的功率，一般選用一容量足夠大的發(fā)電廠（通常是承擔(dān)系統(tǒng)調(diào)頻任務(wù)的發(fā)電廠）來擔(dān)任。一個獨立的電力網(wǎng)絡(luò)只設(shè)一個平衡節(jié)點。41第41頁，課件共153頁，創(chuàng)作于2023年2月注意：三類節(jié)點的劃分并不是絕對不變的。PV節(jié)點之所以能控制其節(jié)點的電壓為某一設(shè)定值，重要原因在于它具有可調(diào)節(jié)的無功功率出力。一旦它的無功功率出力達到可調(diào)節(jié)的上限或下限，就不能使電壓保持在設(shè)定值，PV節(jié)點將轉(zhuǎn)化成PQ節(jié)點。42第42頁，課件共153頁，創(chuàng)作于2023年2月43第43頁，課件共153頁，創(chuàng)作于2023年2月三.潮流方程的求解高斯迭代法牛頓-拉夫遜法（N-R法）快速解耦法（P-Q分解法）44第44頁，課件共153頁，創(chuàng)作于2023年2月（一）.高斯-塞德爾法潮流計算迭代法隱式方程反復(fù)迭代確定數(shù)列{xk}有極限高斯迭代猜測值45第45頁，課件共153頁，創(chuàng)作于2023年2月迭代法是一種逐次逼近法基本思想是將隱式方程歸結(jié)為一組顯式的計算公式實質(zhì)上是一個逐步顯式化的過程46第46頁，課件共153頁，創(chuàng)作于2023年2月高斯迭代法潮流計算假設(shè)有n個節(jié)點的電力系統(tǒng)，沒有PV節(jié)點，平衡節(jié)點編號為s。47第47頁，課件共153頁，創(chuàng)作于2023年2月高斯-塞德爾法潮流計算最新計算出來的第k+1次近似分量加以利用

本次迭代新值48第48頁，課件共153頁，創(chuàng)作于2023年2月高斯-塞德爾法潮流計算PV節(jié)點做第k+1次迭代前如果所求得PV節(jié)點的無功功率越限，則無功功率在限，該PV節(jié)點轉(zhuǎn)化為PQ節(jié)點。設(shè)定的Up0對電壓進行修正，但其相角仍保持上式所求得的值49第49頁，課件共153頁，創(chuàng)作于2023年2月（二）.牛頓-拉夫遜法潮流計算牛頓-拉夫遜法 牛頓-拉夫遜法是求解非線性代數(shù)方程有效的迭代計算方法。在牛頓-拉夫遜法的每一次迭代過程中，非線性問題通過線性化逐步近似。50第50頁，課件共153頁，創(chuàng)作于2023年2月牛頓-拉夫遜法設(shè)解的初值為x0，與真解的誤差為Δx0泰勒展開非線性函數(shù)雅可比因子51第51頁，課件共153頁，創(chuàng)作于2023年2月牛頓-拉夫遜法幾何解釋方程f(x)=0的根x*可解釋為曲線y=f(x)與x軸的交點的橫坐標。設(shè)xk是根x*

的某個近似值，過曲線y=f(x)上橫坐標為xk的點Pk引切線，并將該切線與x軸的交點的橫坐標xk+1做為x*的新的近似值。52第52頁，課件共153頁，創(chuàng)作于2023年2月牛頓-拉夫遜法將單變量問題推廣到具有n個未知變量的X的n階非線性聯(lián)立代數(shù)方程組F（X）其中：J為函數(shù)向量F（X）對變量X的一階偏導(dǎo)數(shù)的雅可比矩陣，是n階方陣。每次迭代的修正量為：53第53頁，課件共153頁，創(chuàng)作于2023年2月牛頓-拉夫遜法收斂條件：54第54頁，課件共153頁，創(chuàng)作于2023年2月注意：55第55頁，課件共153頁，創(chuàng)作于2023年2月牛頓-拉夫遜法計算潮流F(X)是節(jié)點功率方程：X是節(jié)點的電壓復(fù)值和相角：56第56頁，課件共153頁，創(chuàng)作于2023年2月節(jié)點電壓以極坐標形式表示功率方程：57第57頁，課件共153頁，創(chuàng)作于2023年2月節(jié)點電壓以極坐標形式表示對功率方程求導(dǎo)，得到修正方程為：其中雅可比矩陣的各元素分別為：58第58頁，課件共153頁，創(chuàng)作于2023年2月雅可比矩陣元素59第59頁，課件共153頁，創(chuàng)作于2023年2月修正方程中對各類節(jié)點的處理：PQ節(jié)點：每個PQ節(jié)點有兩個變量待求，都要參加聯(lián)立求解；PV節(jié)點：節(jié)點電壓給定，為零，只有一個變量。因此，該類節(jié)點只有有功部分參加聯(lián)立求解，而雅可比矩陣中該類節(jié)點無功部分則除去相應(yīng)的行和列，但每次迭代完成需計算該節(jié)點的無功功率，以校驗是否越限；60第60頁，課件共153頁，創(chuàng)作于2023年2月修正方程中對各類節(jié)點的處理：平衡節(jié)點：因其電壓大小、相位均為已知，所以不需要參加聯(lián)立求解，一般處理為，在雅可比矩陣中對應(yīng)該節(jié)點的對角元素為一大數(shù)，其他部分為0，當?shù)Y(jié)束后再求該節(jié)點的有功功率和無功功率。61第61頁，課件共153頁，創(chuàng)作于2023年2月方程組維數(shù)的確定：100個節(jié)點，有10個PV節(jié)點，雅可比矩陣的維數(shù)？

2*（100-1）-10=188n個節(jié)點，有m個PV節(jié)點，雅可比矩陣的維數(shù)：

2*（n-1）-m62第62頁，課件共153頁，創(chuàng)作于2023年2月修正方程：63第63頁，課件共153頁，創(chuàng)作于2023年2月節(jié)點電壓以直角坐標形式表示功率方程64第64頁，課件共153頁，創(chuàng)作于2023年2月節(jié)點電壓以直角坐標形式表示對功率方程求導(dǎo)，得到修正方程為：其中雅可比矩陣的各元素分別為：65第65頁，課件共153頁，創(chuàng)作于2023年2月修正方程中對各類節(jié)點的處理：PQ節(jié)點：每個PQ節(jié)點有兩個變量待求，都要參加聯(lián)立求解；PV節(jié)點：節(jié)點電壓有效值給定，它們之間的關(guān)系為：，用這個關(guān)系式來代替該節(jié)點無功功率表達式，并改變雅可比矩陣中對應(yīng)該節(jié)點相應(yīng)的部分；平衡節(jié)點：因其電壓大小、相位均為已知，所以不需要參加聯(lián)立求解。66第66頁，課件共153頁，創(chuàng)作于2023年2月方程組維數(shù)的確定：100個節(jié)點，有10個PV節(jié)點，雅可比矩陣的維數(shù)？

2*（100-1）=198n個節(jié)點，有m個PV節(jié)點，雅可比矩陣的維數(shù)：

2*（n-1）67第67頁，課件共153頁，創(chuàng)作于2023年2月節(jié)點電壓以完全極坐標形式表示功率方程（節(jié)點電壓和節(jié)點導(dǎo)納矩陣都以極坐標形式表示）：68第68頁，課件共153頁，創(chuàng)作于2023年2月雅可比矩陣元素69第69頁，課件共153頁，創(chuàng)作于2023年2月雅可比矩陣的特點雅可比矩陣為一非奇異方陣。傳統(tǒng)的，當節(jié)點電壓以極坐標表示時，該矩陣為2(n-1)-m階方陣（m為PV節(jié)點數(shù)）；當節(jié)點電壓以直角坐標表示時，該矩陣為2(n-1)階方陣。現(xiàn)在，為了便于編程，一般為經(jīng)過處理的2n階。矩陣元素與節(jié)點電壓有關(guān)，故每次迭代時都要重新計算。與導(dǎo)納矩陣具有相似的結(jié)構(gòu)，當Yij=0，Hij、Nij、Jij、Lij均為0，因此也是高度稀疏的矩陣。具有結(jié)構(gòu)對稱性，但數(shù)值不對稱。70第70頁，課件共153頁，創(chuàng)作于2023年2月輸入數(shù)據(jù)形成節(jié)點導(dǎo)納矩陣求修正方程式的常數(shù)項向量求出雅可比矩陣元素求解修正方程式，求出變量的修正向量求出節(jié)點電壓的新值檢查是否收斂計算支路功率分布，PV節(jié)點無功功率和平衡節(jié)點注入功率YN71第71頁，課件共153頁，創(chuàng)作于2023年2月牛頓-拉夫遜法的收斂特性牛頓-拉夫遜法具有平方收斂特性，高斯-塞德爾法為一階收斂特性。牛頓-拉夫遜法對初值設(shè)定很敏感。因此，在實際應(yīng)用當中，常常在牛頓-拉夫遜法計算潮流以前先用對初值不敏感的高斯-塞德爾法（迭代1-2次）計算電壓的初值。72第72頁，課件共153頁，創(chuàng)作于2023年2月輸入數(shù)據(jù)（完全極坐標形式）73第73頁，課件共153頁，創(chuàng)作于2023年2月形成節(jié)點導(dǎo)納矩陣74第74頁，課件共153頁，創(chuàng)作于2023年2月形成節(jié)點導(dǎo)納矩陣幅值相角75第75頁，課件共153頁，創(chuàng)作于2023年2月設(shè)初值（平啟動）76第76頁，課件共153頁，創(chuàng)作于2023年2月解修正方程常數(shù)向量77第77頁，課件共153頁，創(chuàng)作于2023年2月解修正方程常數(shù)向量78第78頁，課件共153頁，創(chuàng)作于2023年2月解修正方程常數(shù)向量第一次迭代79第79頁，課件共153頁，創(chuàng)作于2023年2月雅可比矩陣

80第80頁，課件共153頁，創(chuàng)作于2023年2月雅可比矩陣81第81頁，課件共153頁，創(chuàng)作于2023年2月求解修正方程式82第82頁，課件共153頁，創(chuàng)作于2023年2月求解修正方程式83第83頁，課件共153頁，創(chuàng)作于2023年2月求出變量的修正向量δ的單位是度，計算時是弧度84第84頁，課件共153頁，創(chuàng)作于2023年2月85第85頁，課件共153頁，創(chuàng)作于2023年2月收斂條件86第86頁，課件共153頁，創(chuàng)作于2023年2月節(jié)點功率不平衡量變化迭代次數(shù)0-0.2773-0.05097-0.5260.01960.51-4.0E-05-0.0438-0.02047-0.024540.0045120.0001-0.00045-0.00042-0.000320.0000831.0E-091.00E-081.00E-091.00E-081.00E-0987第87頁，課件共153頁，創(chuàng)作于2023年2月牛頓-拉夫遜法的收斂特性88第88頁，課件共153頁，創(chuàng)作于2023年2月牛頓-拉夫遜法的收斂特性89第89頁，課件共153頁，創(chuàng)作于2023年2月節(jié)點電壓變化迭代次數(shù)1-0.00881-0.00654-0.10781-0.023680.115132-8.9E-05-0.00879-0.00467-0.011660.0023231.1E-05-0.0001-8.9E-05-0.000160.0023243.3E-09-1.7E-08-2.2E-08-3.3E-088.9E-0890第90頁，課件共153頁，創(chuàng)作于2023年2月牛頓-拉夫遜法的收斂特性91第91頁，課件共153頁，創(chuàng)作于2023年2月高斯-塞德爾法收斂特性92第92頁，課件共153頁，創(chuàng)作于2023年2月高斯-塞德爾法收斂特性93第93頁，課件共153頁，創(chuàng)作于2023年2月（三）.P-Q分解法牛頓-拉夫遜法修正方程：牛頓-拉夫遜法的核心就是反復(fù)形成并求解修正方程式94第94頁，課件共153頁，創(chuàng)作于2023年2月（三）.P-Q分解法牛頓-拉夫遜法的缺點牛頓-拉夫遜法的雅可比矩陣在每一次迭代過程中都有變化，需要重新形成和求解，這占據(jù)了計算的大部分時間，成為牛頓-拉夫遜法計算速度不能提高的主要原因。牛頓-拉夫遜法的修正方程維數(shù)高，為2N，內(nèi)存需量大。95第95頁，課件共153頁，創(chuàng)作于2023年2月（三）.P-Q分解法P-Q分解法（快速解耦法）利用了電力系統(tǒng)的一些特有的運行特性，對牛頓-拉夫遜法做了簡化，以改進和提高計算速度。雅可比矩陣為常數(shù)修正方程解耦，降低維數(shù)P-Q分解法是牛頓-拉夫遜法潮流計算的一種簡化方法。96第96頁，課件共153頁，創(chuàng)作于2023年2月牛拉法簡化形成P-Q法的過程牛頓-拉夫遜法修正方程展開為：97第97頁，課件共153頁，創(chuàng)作于2023年2月牛拉法簡化形成P-Q法的過程考慮到電力系統(tǒng)中有功功率變化主要受節(jié)點電壓相角變化的影響，無功功率變化主要受節(jié)點電壓幅值變化的影響，所以可以近似的忽略電壓幅值變化對有功功率和電壓相位變化對無功功率分布的影響：

98第98頁，課件共153頁，創(chuàng)作于2023年2月牛拉法簡化形成P-Q法的過程根據(jù)電力系統(tǒng)的正常運行條件還可作下列假設(shè)：電力系統(tǒng)正常運行時線路兩端的電壓相位角一般變化不大（不超過10~20度）；電力系統(tǒng)中一般架空線路的電抗遠大于電阻；節(jié)點無功功率相應(yīng)的導(dǎo)納Q/U*U遠小于該節(jié)點的自導(dǎo)納的虛部。99第99頁，課件共153頁，創(chuàng)作于2023年2月牛拉法簡化形成P-Q法的過程用算式表示如下：+=1=0100第100頁，課件共153頁，創(chuàng)作于2023年2月牛拉法簡化形成P-Q法的過程U為節(jié)點電壓有效值的對角矩陣B為電納矩陣（由節(jié)點導(dǎo)納矩陣中各元素的虛部構(gòu)成）101第101頁，課件共153頁，創(chuàng)作于2023年2月牛拉法簡化形成P-Q法的過程根據(jù)不同的節(jié)點還要做一些改變：在有功功率部分，要除去與有功功率和電壓相位關(guān)系較小的因素，如不包含各輸電線路和變壓器支路等值Π型電路的對地電納。在無功功率部分，PV節(jié)點要做相應(yīng)的處理。一般，對稱的常數(shù)矩陣B’和B’’是不相同的102第102頁，課件共153頁，創(chuàng)作于2023年2月收斂條件103第103頁，課件共153頁，創(chuàng)作于2023年2月P-Q分解法的特點：以一個n-1階和一個n-m-1階線性方程組代替原有的2n-m-1階線性方程組；修正方程的系數(shù)矩陣B’和B”為對稱常數(shù)矩陣，且在迭代過程中保持不變；P-Q分解法具有線性收斂特性，與牛頓-拉夫遜法相比，當收斂到同樣的精度時需要的迭代次數(shù)較多；P-Q分解法一般只適用于110KV及以上電網(wǎng)的計算。因為35KV及以下電壓等級的線路r/x比值很大，不滿足上述簡化條件，可能出現(xiàn)迭代計算不收斂的情況。104第104頁，課件共153頁，創(chuàng)作于2023年2月105第105頁，課件共153頁，創(chuàng)作于2023年2月收斂特性比較106第106頁，課件共153頁，創(chuàng)作于2023年2月作業(yè)2簡單電力系統(tǒng)如圖所示，已知各段線路阻抗和節(jié)點功率為：。節(jié)點1為平衡節(jié)點，，試用牛頓-拉夫遜法計算潮流。

1）形成節(jié)點導(dǎo)納矩陣。

2）求第一次迭代用的雅克比矩陣。

3）求解第一次的修正方程。132G110kVSLD3SLD2107第107頁，課件共153頁，創(chuàng)作于2023年2月（五）相關(guān)技術(shù)初值設(shè)定解方程——高斯消去法——LU分解節(jié)點優(yōu)化編號稀疏技術(shù)結(jié)果分析不收斂分析潮流的發(fā)展潮流軟件的介紹108第108頁，課件共153頁，創(chuàng)作于2023年2月初值設(shè)定平啟動

用高斯-塞德爾迭代法處理109第109頁，課件共153頁，創(chuàng)作于2023年2月解方程110第110頁，課件共153頁，創(chuàng)作于2023年2月解方程-高斯消去法111第111頁，課件共153頁，創(chuàng)作于2023年2月解方程-高斯消去法矩陣A元素表達式（消元過程）112第112頁，課件共153頁，創(chuàng)作于2023年2月解方程-高斯消去法求解三角形方程組，得到求解公式（回代過程）：113第113頁，課件共153頁，創(chuàng)作于2023年2月解方程-LU分解114第114頁，課件共153頁，創(chuàng)作于2023年2月解方程-LU分解

對于r=2,3,…n計算：計算U的第r行元素：計算L的第r列元素：

115第115頁，課件共153頁，創(chuàng)作于2023年2月節(jié)點優(yōu)化編號節(jié)點編號順序與稀疏度的關(guān)系41321432經(jīng)過3次消去經(jīng)過1次消去經(jīng)過3次消去經(jīng)過2次消去116第116頁，課件共153頁，創(chuàng)作于2023年2月節(jié)點優(yōu)化編號節(jié)點的編號順序?qū)τ谟嬎阈ЯΦ挠绊懼陵P(guān)重要，特別是采用了稀疏技術(shù)后，它直接影響到矩陣A的因子表矩陣的稀疏度。最優(yōu)編號是一個組合優(yōu)化問題，求其最優(yōu)解是困難的，但在實際工程中，有許多實用的次優(yōu)的編號方法得到了廣泛的應(yīng)用。

117第117頁，課件共153頁，創(chuàng)作于2023年2月Tinney-1編號方法

也稱靜態(tài)節(jié)點優(yōu)化編號法。這種方法在有向圖上統(tǒng)計每個節(jié)點的出線度，即該節(jié)點和其他節(jié)點相連接的支路樹，然后按節(jié)點出線度由大到小按順序進行編號。方法簡單，但編號效果較差。

118第118頁，課件共153頁，創(chuàng)作于2023年2月Tinny-2編號方法

這種方法也稱最小度算法，或半動態(tài)節(jié)點優(yōu)化編號方法。首先統(tǒng)計所有節(jié)點的出線度，然后選擇出線度最小的節(jié)點進行編號。編號過程中，按圖上因子分解的方法消去該節(jié)點，只進行網(wǎng)絡(luò)結(jié)構(gòu)變化的處理，而不進行邊權(quán)計算。然后消去已編號的節(jié)點和其相關(guān)支路，在剩下的子圖上重復(fù)上述編號過程。119第119頁，課件共153頁，創(chuàng)作于2023年2月Tinney-2編號方法

這種方法也比較簡單，圖上因子分解產(chǎn)生新支路以及處理過的支路這些變化可用在原來的圖上修正來實現(xiàn)。這種編號方法可使有向因子圖上新增加的支路數(shù)大大減少，而程序復(fù)雜性和計算量又增加不多，是一種使用十分廣泛的編號方法。120第120頁，課件共153頁，創(chuàng)作于2023年2月Tinney-3編號方法這種方法也稱動態(tài)節(jié)點優(yōu)化編號方法。它和上面的Tinney-2編號方法的不同之處是對所有待編號的節(jié)點，統(tǒng)計消去該節(jié)點時產(chǎn)生的新支路的數(shù)目，并以該數(shù)目最小為優(yōu)先編號的準則。某一節(jié)點編號完成之后，要立即修改因子圖。

121第121頁，課件共153頁，創(chuàng)作于2023年2月稀疏技術(shù)稀疏矩陣的存儲特點是排零存儲，即只存儲其中的非零元素和有關(guān)的檢索信息。存儲的目的是為了在計算中能方便地訪問使用。要求所采用的存儲格式既節(jié)省內(nèi)存，又能夠方便地檢索和存取，同時還要考慮網(wǎng)絡(luò)矩陣結(jié)構(gòu)變化時能方便地對存儲的信息加以修改。

122第122頁，課件共153頁，創(chuàng)作于2023年2月稀疏技術(shù)線性表單鏈表雙鏈表位結(jié)構(gòu)123第123頁，課件共153頁，創(chuàng)作于2023年2月結(jié)果分析124第124頁，課件共153頁，創(chuàng)作于2023年2月結(jié)果分析線路功率損耗125第125頁，課件共153頁，創(chuàng)作于2023年2月不收斂分析可能出現(xiàn)的問題（系統(tǒng)問題）網(wǎng)絡(luò)參數(shù)多島或孤島標幺值特別大（?。┑闹?/p>

單線變壓器支路過長運行參數(shù)（調(diào)度方式）有功不平衡無功不平衡126第126頁，課件共153頁，創(chuàng)作于2023年2月不收斂分析相應(yīng)的調(diào)試方法（系統(tǒng)問題）放開收斂判據(jù)記錄分析不平衡量，根據(jù)不平衡量出現(xiàn)的節(jié)點號查找相應(yīng)的網(wǎng)絡(luò)參數(shù)和運行參數(shù)直流潮流，分區(qū)檢查缺額

127第127頁，課件共153頁，創(chuàng)作于2023年2月不收斂分析可能出現(xiàn)的問題（初始條件）計算的初始條件設(shè)置不合理PV點的設(shè)置平衡機的選點128第128頁，課件共153頁，創(chuàng)作于2023年2月不收斂分析相應(yīng)的調(diào)試方法（初始條件）應(yīng)盡量選擇無功調(diào)節(jié)能力強的發(fā)電機作為PV節(jié)點PV節(jié)點的電壓應(yīng)設(shè)在對應(yīng)的發(fā)電機額定電壓小的鄰域內(nèi),保證全網(wǎng)各處的電壓水平相當PV節(jié)點應(yīng)在電網(wǎng)中均勻分布,防止大量無功在系統(tǒng)中長距離流動造成大的電壓降129第129頁，課件共153頁，創(chuàng)作于2023年2月不收斂分析相應(yīng)的調(diào)試方法（初始條件）在1個廠站內(nèi)最好只選1個PV節(jié)點，應(yīng)在臨近廠站上設(shè)多個PV節(jié)點,因為這要求兩節(jié)點間的潮流滿足規(guī)定的電壓差,很可能是不合理的平衡機盡可能選在電網(wǎng)重心附近，即到全網(wǎng)各母線的電氣距離盡量短130第130頁，課件共153頁，創(chuàng)作于2023年2月6.不收斂分析3）保證潮流收斂的人工干預(yù)手段要重點考察的兩個指標:系統(tǒng)最低電壓和平衡機功率必要的電壓調(diào)節(jié):提高發(fā)電機PQ節(jié)點的無功出力降低負荷的無功消耗投退無功設(shè)備協(xié)調(diào)平衡機功率:如果平衡機上有大的功率吞吐,應(yīng)將此不平衡功率預(yù)先分配到協(xié)調(diào)機組上去,減輕平衡機的輸送壓力,防止出現(xiàn)大功率長距離輸送的情況。131第131頁，課件共153頁，創(chuàng)作于2023年2月潮流的發(fā)展發(fā)展的方向（算法）改善收斂提高速度減少內(nèi)存132第132頁，課件共153頁，創(chuàng)作于2023年2月潮流的發(fā)展研究的方向（問題）動態(tài)潮流諧波潮流狀態(tài)估計潮流概率潮流最優(yōu)潮流133第133頁，課件共153頁，創(chuàng)作于2023年2月8.潮流計算軟件介紹國際上幾種電力系統(tǒng)分析計算軟件包134第134頁，課件共153頁，創(chuàng)作于2023年2月BPA潮流計算程序美國幫涅維爾電力局（BPA，BonnevillePowerAdministr-ation）開發(fā)；被中國電力科學(xué)院引進吸收，從1984年開始在中國得到推廣應(yīng)用。程序提供兩種潮流計算方法：P_Q分解法和牛頓法。135第135頁，課件共153頁，創(chuàng)作于2023年2月PSASP潮流計算程序簡介：中國電力科學(xué)院開發(fā)。程序提供五種潮流計算方法：P_Q分解法、牛頓法(功率式)、最佳乘子法、牛頓法（電流式）、P_Q分解法轉(zhuǎn)牛頓法(電流式)136第136頁，課件共153頁，創(chuàng)作于2023年2月（六）電力系統(tǒng)潮流計算程序設(shè)計以牛頓-拉夫遜法為例137第137頁，課件共153頁，創(chuàng)作于2023年2月

開始讀入網(wǎng)絡(luò)參數(shù)及系統(tǒng)運行參數(shù)數(shù)據(jù)文件，形成節(jié)點導(dǎo)納矩陣初始化，設(shè)定電壓初值V（0）=1，δ（0）=0，Kmax=50，k=0

計算節(jié)點功率不平衡量ΔP和ΔQ，max|ΔPi，ΔQi|=e，i=1,2,…n

判斷k>Kmax或e<ε

形成雅可比矩陣，修正方程，求出V，δ

計算節(jié)點功率不平衡量ΔP和ΔQ，max|ΔPi，ΔQi|=e，i=1,2,…n

k>Kmax計算支路功率，輸出潮流結(jié)果

結(jié)束YNYN第138頁，課件共153頁，創(chuàng)作于2023年2月數(shù)據(jù)文件線路參數(shù) 在電力系統(tǒng)程序設(shè)計中，線路參數(shù)一般采用線路的Π型數(shù)學(xué)模型，即線路用節(jié)點間的阻抗和節(jié)點對地容性電納來表示，由于線路的對地電導(dǎo)很小，一般可忽略不計。對于線路參數(shù)的數(shù)據(jù)文件格式一般可寫為：

線路參數(shù)（序號，節(jié)點i，節(jié)點j，r，x，b/2）139第139頁，課件共153頁，創(chuàng)作于2023年2月數(shù)據(jù)文件變壓器參數(shù) 在電力系統(tǒng)程序設(shè)計中，變壓器參數(shù)一般采用Π型等值變壓器模型，這是一種可等值地體現(xiàn)變壓器電壓變換功能的模型。在多電壓級網(wǎng)絡(luò)計算中采用這種變壓器模型后，就可不必進行參數(shù)和變量的歸算。 對于變壓器參數(shù)的數(shù)據(jù)文件格式一般可寫為：

變壓器參數(shù)（序號，節(jié)點i，節(jié)點j，r，x，k0） 其中，k0表示變壓器變比。140第140頁，課件共153頁，創(chuàng)作于2023年2月數(shù)據(jù)文件對地支路參數(shù)對地支路參數(shù)一般以導(dǎo)納形式表示，其等價回路如下：對地支路參數(shù)的數(shù)據(jù)文件格式一般可寫為：

接地支路參數(shù)（序號，節(jié)點i，gi，bi）

141第141頁，課件共153頁，創(chuàng)作于2023年2月數(shù)據(jù)文件節(jié)點功率參數(shù)節(jié)點功率參數(shù)的數(shù)據(jù)文件格式一般可寫為：

節(jié)點功率數(shù)據(jù)（序號，節(jié)點i，PGi，QGi，PLi，QLi）142第142頁，課件共153頁，創(chuàng)作于2023年2月數(shù)據(jù)文件三類節(jié)點參數(shù)平衡節(jié)點：給出節(jié)點編號，節(jié)點電壓。PQ節(jié)點：在節(jié)點功率參數(shù)中就可表示。PV節(jié)點：需單列，其數(shù)據(jù)文件格式一般可寫為：

PV節(jié)點數(shù)據(jù)（序號，節(jié)點i，電壓Vi，無功功率下限，無功功率上限）。143第143頁，課件共153頁，創(chuàng)作于2023年2月一般潮流數(shù)據(jù)文件格式

節(jié)點數(shù)，平衡節(jié)點，平衡節(jié)點電壓，計算精度

線路參數(shù)（序號，節(jié)點i，節(jié)點j，r，x，b/2）

變壓器參數(shù)（序號，節(jié)點i，節(jié)點j，r，x，k0）

接地支路參數(shù)（序號，節(jié)點i，gi，bi）

節(jié)點功率數(shù)據(jù)（序號，節(jié)點i，