專利名稱：大規模電力系統網絡分割潮流并行計算方法
技術領域：
本發明涉及電力系統運行及在線安全分析技術領域，特別是一種大規模電力系統潮流并行計算方法。
背景技術：
潮流計算是電力系統運行及在線安全分析中最基本的計算任務之一。潮流計算是根據給定電力網絡結構及運行條件計算整個網絡的潮流分布，計算結果中最基本的電氣量是系統各母線的復電壓(稱作系統狀態)的穩態值，其內容屬于電力系統穩態(或靜態)分析。其常規的數學模型和求解方法簡要介紹如下電力系統潮流計算的數學模型通常由式(1)所示非線性代數方程組表示，其中向量X表示非線性代數方程組的解。在數學模型大規模電力系統上潮流計算屬于超大規模非線性代數方程求解問題。
G(X)＝0(1)給定變量X的預測值X0(又稱作初值)后，可通過牛頓迭代法線性化后的牛頓迭代方程式(2)和式(3)求解X的數值解，其中上標k表示迭代次數，J(Xk)為第k次迭代中的雅可比矩陣(Jacobian)矩陣。
J(Xk)ΔXk＝-G(Xk) (2)Xk+1＝Xk+ΔXk(3)當修正量ΔXk元素中的最大絕對值小于某一門檻值ε，例如ε＜10-5，(3)式所得Xk+1作為X的數值解。
上述電力系統潮流計算需要在計算機上進行，其計算量和電力系統的規模相關，電力系統規模越大，潮流計算中的計算量越大，需要的計算時間越多。為提高電力系統潮流計算的速度，在技術層面上可采用(在并行或集群計算機上進行的)分布式計算或并行處理技術；在算法層面上則需要研究計算量少、適于并行處理的新算法。
式(2)和式(3)所示牛頓法迭代法是一種常用的電力系統潮流計算方法，該方法具有收斂快速和不引入交接誤差的優點。但是，常規的聯立求解法將得到的所有非線性代數方程集中在一起進行整體處理，建立一個整體雅可比矩陣方程并加以求解，由于沒有考慮電力系統的網絡結構特點，每個局部網絡的變化都會影響到整體雅可比矩陣方程修改，因此可擴展性和靈活性不足。而且現有的各種聯立求解法在每次迭代中都需要對其整體雅可比矩陣進行代數操作，計算量較大，在計算效率上不能很好地滿足大規模互聯電力系統實時計算的要求。此外，常規的聯立求解法不易實現并行化計算。

發明內容本發明的目的是解決現有聯立求解算法存在的計算量較大、不能很好地滿足大規模互聯電力系統實時計算要求等問題，提供一種大規模電力系統網絡分割潮流并行計算方法。
本發明是一種大規模電力系統潮流并行計算技術，在計算原理上屬于經典的基于牛頓法的聯立求解法，但在技術實現上，它通過網絡劃分，將整體電力系統高維方程組分解為低維的多個方程組來處理，有效地減少計算量。本發明方法包括以下步驟步驟101將大規模或超大規模電力網絡按照選定的聯絡線集合分割為適度規模的一組子網。考慮大規模電網具有分層互聯的特點(例如我國電網由網級電網互連而成，網級電網由省級電網互連而成，省級電網由地級電網互連而成)，在選擇網絡分割聯絡線時應考慮電網分層互聯特點。
步驟102根據網絡分割結果，構建電網潮流計算樹。潮流計算樹由最上層的一個“根結點”、“中間根結點”和多個最下層的“葉結點”(末梢結點)構成。每一葉結點對應一個電力子網絡，根結點和中間根結點對應一種網絡分割。
步驟103對葉結點對應的電力子網絡中的潮流分析中的變量，即待求母線電壓變量進行分類。依據其母線和參與其網絡分割的聯絡線關系，將其中所有待求變量劃分為三類邊界輸出母線電壓變量、邊界輸入母線電壓變量和內部母線電壓變量。
步驟104給定系統母線電壓(稱作狀態變量)的初始值。例如，所有母線(包括潮流分析中的參考母線)電壓模值初始值為1(標幺值單位)，相角初始值為0(弧度單位)。
步驟105列寫每個葉結點對應網絡的潮流代數方程(即PQ母線的功率平衡方程和PV結點的有功和電壓平衡方程)，對其進行泰勒(Taylor)級數展開，略去Taylor級數中二階及以上的高階項，得到子系統數學上的雅可比矩陣方程。雅可比矩陣方程由雅可比矩陣、母線電壓的修正量向量和不平衡向量組成，在每次迭代中，雅可比矩陣和不平衡向量都是確定的實數矩陣和向量，需要求解的是母線電壓的修正量向量。
步驟106計算母線電壓的修正量向量。具體步驟為106-1.上行替換對于葉結點該上行替換從所有葉結點開始。對于葉結點，通過對雅可比矩陣方程進行初等代數變換(如LU分解和回代計算)，將葉結點中的邊界輸出和內部電壓變量的修正量都寫成由邊界輸入電壓變量的修正量表示的線性表達式。
對于中間根結點(如果有的話)中間根結點要作的上行替換是接受下層結點的邊界母線電壓修正量線性表達式，消掉該中間根結點對應的網絡分割的邊界母線電壓修正量，獲得上一層根結點需要的邊界母線電壓的修正量線性表達式。
對于根結點，合并其所有下層葉結點的邊界輸出、輸入電壓修正量表示的線性方程組。并求解該方程組，解得邊界輸出、輸入電壓修正量。
106-2.下行計算該操作由根節點至潮流計算樹所有葉結點進行。將解得的邊界電壓修正量帶入上行替換所得內部(或下一層邊界輸入輸出電壓修正量)和邊界輸入電壓修正量的線性表達式并最后求得各子網絡內部電壓修正向量。
步驟107進行收斂判斷(判斷變量的修正量向量中的最大值或者不平衡向量中的最大值是否小于給定的容許閾值，若是肯定的，判定計算已收斂)。如果不收斂，使用求出的變量修正量對變量進行修正，并轉步驟105迭代計算；步驟108輸出收斂后的電力系統潮流計算結果。
步驟109結束潮流計算。
本發明的優點和積極效果本發明方法具有如下優點，a、大幅降低了使用牛頓法對大規模電力系統潮流方程進行聯立求解的計算量。該算法利用網絡的分層互聯特點將全網的高維稀疏方程組的求解工作分配到各個網絡計算結點中，將高維電力系統方程組的整體求解問題分解為多個低維方程組的計算問題。阻斷了高維稀疏方程組求解不斷產生新填元的問題，提高了計算效率。
b、由于采用樹形計算結構，潮流計算中同層不同結點上的電壓修正量計算過程互不依賴，便于采用分布集成式計算機并行進行。
c、該算法對每個電力系統子網絡單獨建立數學模型，單獨列寫雅可比方程，因此比常規的聯立求解算法具有更強的模塊性、更好的擴展性、更便于編程實現。
d、該發明提出的網絡分割并行計算技術，可擴展應用于大規模電力系統仿真中求解差分化后的電力系統非線性方程組，實現大規模電力系統的在線超時實仿真。

圖1是本發明的電力系統潮流計算基本算法流程示意圖；圖2是本發明中目標電力網絡的一種分割方式；圖3是本發明的電力網絡潮流計算樹構成示意圖；圖4是典型新英格蘭電力系統結構圖；圖5是典型新英格蘭電力系統第一次網絡分割結果；圖6是典型新英格蘭電力系統第一次網絡分割后得到的潮流計算樹；圖7是典型新英格蘭電力系統第二次網絡分割后結果；圖8是典型新英格蘭電力系統第二次網絡分割后得到的潮流計算樹。
具體實施方式實施例1電力系統由生產、輸送、分配和消費電能的各種電氣設備組成。電力網絡由變壓器、不同電壓等級的輸電線路及其他輸電裝置組成。電力網絡具有較強的地域分區特性，一個電力系統可以視為由多個子網由聯絡線互聯組成。
下面結合附圖對本發明作進一步詳細描述，然后通過集群機系統介紹其并行實施方式。附圖1是本發明的電力系統潮流計算基本算法流程示意圖。
本發明的特征是包括以下步驟步驟101對電力網絡進行劃分，基于大規模電網為地區電網經聯絡線互聯而成的特點，本發明提出通過聯絡線和將電力網絡劃分為適度規模的多個子網的方法。電網劃分結果應滿足以下條件(1)劃分后的電力子網按其網絡聯結方式可組織為多層樹形結構，最上層只有一個根結點，最下層葉結點對應網絡分割后所得的所有子網絡。根結點和中間根結點對應網絡劃分的聯絡線和與其關聯的母線，一般其網絡圖形是不連通的。
(2)每個子網通過選定的聯絡線與一個或多個子網相連接。
(3)與選定的聯絡線相關聯的母線稱為該子網的邊界母線。
圖2給出了一種電力網絡分割方式。如圖2所示，目標電網被分割成共6個電力子網。標號B1~B11表示原始網絡中與6條選定的聯絡線關聯的11個邊界母線。該6條聯絡線(包括其母線)構成潮流計算樹的根結點。其余N1~N5五個子網絡構成五個葉結點。在應用本發明時，在電力網絡沒有發生改變的情況下，網絡分區只需要進行一次即可。分區結果在潮流計算中重復使用。
步驟102根據網絡劃分的結果，構建用潮流計算樹。對照圖2所示分割建立的潮流計算樹如圖3所示。若N1~N5五個子網絡依然龐大，各子網可按圖1所示方法，選定聯絡線進行子網絡劃分，建立子網絡潮流計算樹，子網絡潮流計算樹的根結點為中間根結點。
步驟103對各個結點對應子網絡中的待求電壓變量進行分類。對照圖3所示系統潮流計算樹，我們將B1~B11母線電壓選作邊界輸入輸出變量，其余母線電壓為各子網絡內部變量。
步驟104給定所有電壓變量的初始值。即，所有的待求母線電壓模值初始值為1.0(標幺值單位)，相角初始值為0.0(弧度單位)。
步驟105建立各子網的數學模型。采用直角或極坐標形式(以下敘述按直角坐標形式敘述)，每個子網潮流方程可表示為GNi(Vin，Vbi，Vbo)＝0(4)其中Vin表示子網絡內部母線電壓實部和虛部構成的向量(以下簡稱母線電壓向量)，Vbi為網絡內部邊界母線電壓實部和虛部構成的向量(稱作子網絡潮流計算輸出變量)，Vbo為網絡外部邊界母線電壓實部和虛部構成的向量(稱作子網絡潮流計算輸入變量)。令向量Vin，Vbi和Vbo的維數分別為nin，nbi和nbo，則方程(4)的維數為nin+nbi。對迭代中的值Vink，Vbik和Vbok其線性化后的牛頓迭代增量方程如式(5)所示，其中k表示迭代次數，J(Vink，Vbik，Vbok)為第k次迭代中子網絡潮流方程的雅可比矩陣。
J(Vink,Vbik,Vbok)ΔVinkΔVbikΔVbok=-G(Vink,Vbik,Vbok)···(5)]]>步驟106計算修正量。
106-1.上行替換該操作針對潮流計算樹的所有葉結點進行。對照式(5)，通過對雅可比矩陣方程進行初等代數變換(如LU分解和回代計算)，將葉結點中的邊界輸出和內部變量的修正量都寫成由邊界輸入修正量表示的如下標準形式。
ΔVink+PkΔVbok=Qk···(6)]]>ΔVbik+RkΔVbok=Sk···(7)]]>式(7)稱作上行替換方程，式(6)稱作下行回代方程。在每次迭代中Pk，Rk，Qk和Sk都是確定的實數矩陣和向量，需要求解的是電壓變量的修正量向量。
106-2.邊界輸出、輸入修正量計算對于根結點和中間根結點，合并其所有下層葉結點的標準形式方程(7)得到計算所有邊界輸出、輸入修正量的線性方程組，求解該方程組，可解得所有邊界輸入、輸出修正量ΔVbik和ΔVbok。
106-3.下行計算該操作針對潮流計算樹的所有葉結點進行。將解得的邊界輸入修正量ΔVbok帶入上行替換后所得方程(6)得子網絡內部修正向量ΔVink。
為了加快大規模電力系統計算的速度，人們往往采用由多個計算機組成的集群系統作為計算的硬件設施。近年來集群系統已成為實現并行計算的一種新主流技術。集群系統由多個計算處理機(包括一個主機和若干子機)通過高速網絡互聯而成，各個計算處理機之間采用數據通信機制互相聯絡。下面介紹本發明在集群計算機系統上的并行實施方式。
可以看出，本發明計算對象的數學模型具有樹形的結構，求解各葉結點上行替換的表達式(6)和(7)的過程和下行用式(6)計算子網絡內部修正量ΔVink計算過程互不依賴，可并行進行。
分布式并行潮流計算描述對照圖3所示潮流計算樹，若集群計算機系統的每一臺處理機(包括一個主機和子機)承擔N1至N5各子網絡上行替換和下行計算任務的速度(考慮數據通訊時間)已能滿足電力系統在線安全監控需要，則并行計算可用一個主機和5個子機完成。主機和子機之間通訊上行替換傳送的是Rk(nbi×nbo維矩陣)，和Sk(nbo維向量)，一般維數nbo很小，下行計算傳送的是邊界母線電壓修正量ΔVbik(nbi維向量)和ΔVbok(nbo維向量)，總體數據傳送量十分有限。另一特點是根結點和中間根結點聯立求解各子網絡邊界母線電壓修正量ΔVbik和ΔVbok(見式(7))的計算量十分有限。
若集群計算機系統的某一臺處理機(包括一個主機和子機)承擔某一子網絡上行替換和下行計算任務的速度(考慮數據通訊時間)達不到電力系統在線安全監控需要，則要將該子網絡按本發明的分割方法進一步分割。此時原來子網絡結點變作中間根結點，分割后的子網絡直接作為該中間根結點的葉結點。上述子網絡的分割過程直至集群計算機系統的處理機在潮流計算中達到在線安全分析速度為止。
分布式并行潮流計算的子網絡描述(拓撲結構和元件參數)在每個處理機上存儲。
實施例2以下以圖4所示典型新英格蘭系統說明本發明的網絡分割和構造潮流計算樹的方法。將聯絡線“16-17”，“15-16”，“1-39”和“4-14”選作網絡分割聯絡線，所得的網絡分割結果如圖5所示，對應的潮流計算樹如圖6所示。
采用分布式并行潮流計算，3個子計算機分別承擔N1，N2和N3子網絡上行替換和下行計算任務。對N1~N3子網，建立潮流方程，第k次迭代對雅可比矩陣方程進行LU分解初等代數變換后所得邊界輸出和內部變量標準修正方程為N1下行回代方程ΔVin,N1k+PN1kΔVbo,N1k=QN1k;]]>上行替換方程ΔVbi,N1k+RN1kΔVbo,N1k=SN1k]]>(三個內部邊界母線6個方程)N2下行回代方程ΔVin,N2k+PN2kΔVbo,N2k=QN2k]]>上行替換方程ΔVbi,N2k+RN2kΔVbo,N2k=SN2k]]>(一個內部邊界母線2個方程)N3下行回代方程ΔVin,N3k+PN3kΔVbo,N3k=QN3k]]>上行替換方程ΔVbi,N3k+RN3kΔVbo,N3k=SN3k]]>(三個內部邊界母線6個方程)3個上行替換方程的系數矩陣及右端向量都是實數矩陣，通過數據通訊傳送到主機。3個上行替換方程的個數是14個，待求邊界母線電壓修正量也是14個，求解三個上行替換方程得到14個邊界母線電壓修正量。分別將1，3和17母線電壓修正量(6個)通訊返回一號子計算機；16母線電壓修正量(2個)通訊返回二號子計算機；39，4，15母線電壓修正量(6個)通訊返回三號子計算機；邊界母線電壓修正量帶入三個子計算機就計算出N1，N2和N3三個子網絡所有內部母線電壓修正量。如圖1所示重復上述過程判斷收斂后，潮流計算結束。
假定N3子網規模較大，考慮實時安全分析速度的需要繼續進行分割，選擇線路6-11和4-14作為分割聯絡線，所得的網絡分割結果如圖7所示，對應的兩種潮流計算樹如圖8所示。一種以N1，N2，N5和N6作為葉結點，合并N0和N4結點作為根節點；另一種以N1，N2，N5和N6作為葉結點，NO結點作為根節點，N4結點作為中間根結點。
采用分布式并行潮流計算，由4個子計算機分別承擔N1，N2，N5和N6子網絡上行替換和下行計算任務。對N1，N2，N5和N6子網絡建立潮流方程，第k次迭代對雅可比矩陣方程進行LU分解初等代數變換所得邊界輸出和內部變量標準修正方程為N1下行回代方程ΔVin,N1k+PN1kΔbo,N1k=QN1k;]]>上行替換方程ΔVbi,N1k+RN1kΔVbo,N1k=SN1k]]>(三個內部邊界母線6個方程)N2下行回代方程ΔVin,N2k+PN2kΔVbo,N2k=QN2k;]]>上行替換方程ΔVbi,N2k+RN2kΔVbo,N2k=SN2k]]>(一個內部邊界母線2個方程)N5下行回代方程ΔVin,N5k+PN5kΔVbo,N5k=QN5k;]]>上行替換方程ΔVbi,N5k+RN5kΔVbo,N5k=SN5k]]>(三個內部邊界母線6個方程)N6下行回代方程ΔVin,N6k+PN6kΔVbo,N6k=QN6k;]]>上行替換方程ΔVbi,N6k+RN6kΔVbo,N6k=SN6k]]>(三個內部邊界母線6個方程)
N1，N2，N5和N6子網絡的4個上行替換方程的系數矩陣及右端向量都是實數矩陣，通過數據通訊傳送到主機。4個上行替換方程的個數為20個，待求邊界母線電壓修正量由原來的14個增加了來自6，11，和14母線的6個電壓修正量，待求邊界母線電壓修正量總數也是20個。
對第一種合并N0和N4結點作為根節點的網絡分割，直接求解4個上行替換方程得到20個邊界母線電壓修正量。上述計算在主機進行。
對第二種N0結點作為根節點，N4結點作為中間根結點的網絡分割，計算過程描述如下首先聯立N5和N6上行替換方程，該方程涉及的邊界母線電壓修正量共18個，分別為ΔVbi，N3k(6個變量向量)，ΔVbo，N3k(6個變量向量)和ΔVN5-N6k(二次分割產生(6，11，和14母線)的6個電壓修正量)。整理N5和N6上行替換聯立方程并進行線性變換可得N4結點的上行替換方程(該方程即N3子網的上行替換方程)ΔVbi,N3k+RN3kΔVbo,N3k=SN3k···(7)]]>和下行計算方程ΔVN5+N6k+RN3′kΔVbo,N3k=SN3′k···(8)]]>式(7)和N1和N2的上行替換方程聯立求得第一次分割所有邊界電壓修正量，并將其中的ΔVbo，N3k帶入式(8)計算ΔVN5+N6k同樣可得20個邊界母線電壓修正量。上述計算20個邊界母線電壓修正量的工作皆在主機進行。
分別將1，3和17母線電壓修正量(6個)通訊返回1號子計算機；16母線電壓修正量(2個)通訊返回2號子計算機；15，4，11母線電壓修正量(6個)通訊返回3號子計算機；39，4，6母線電壓修正量(6個)通訊返回4號子計算機；即將邊界母線電壓修正量帶入4個子計算機就計算出N1，N2，N5和N6四個子網絡所有內部母線電壓修正量。如圖1所示重復上述過程判斷收斂后，潮流計算結束。
權利要求
1.一種大規模電力系統網絡分割潮流并行計算方法，其特征是該方法包括1)、將大規模或超大規模電力網絡按照選定的聯絡線集合分割為一組子網；2)、根據網絡分割結果，構建電網潮流計算樹；潮流計算樹由最上層的結點的“根結點”、“中間根結點”和最下層的“葉結點”構成；3)、對分割后的電力子網絡潮流分析中的變量進行分類；4)、列寫各電力子網絡潮流方程；5)、采用牛頓迭代法求解步驟4)得到的各個子網絡的潮流方程，求解過程中變量修正量的計算通過“上行替換”和“下行計算”操作完成；6)、上述變量的修正計算直至牛頓迭代收斂為止，輸出計算結果，結束潮流計算；7)、上述大規模電力系統潮流計算方法在集群計算機系統上并行計算。
2.根據權利要求1所述的潮流并行計算方法，其特征在于步驟1)中，大規模電力系統網絡分割的特點是a)、按照大規模電網分層互聯的特點，分層選擇網絡分割聯絡線將電網分割為一組子網絡；b)、子網絡與子網絡間的電氣聯系通過分層選擇的聯絡線集合描述。
3.根據權利要求1所述的潮流并行計算方法，其特征在于步驟2)中，構建出的電網潮流計算樹具有以下特點a)、每一個“葉結點”對應一個電力子網絡；b)、“根結點”對應進行最第一層網絡分割的聯絡線集合；c)、第一層以后的分割選定的聯絡線對應潮流計算樹的“中間根結點”；d)、電網潮流計算樹描述了潮流分析的邏輯關系，同層結點之間，除了與其上層結點的計算相關外，彼此相互獨立。
4.根據權利要求1所述的潮流并行計算方法，其特征在于步驟3)中，子網絡潮流分析中的待求變量分類具有如下特點a)、子網絡潮流分析中使用的待求變量是其母線電壓變量；b)、依據其母線和參與其網絡分割的聯絡線關系，將其中所有待求變量劃分為三類邊界輸出母線電壓變量、邊界輸入母線電壓變量和內部母線電壓變量；c)、內部母線是和網絡分割聯絡線不關聯的子網絡母線；d)、邊界母線是和網絡分割關聯的聯絡線母線，其屬于該子網絡的聯絡線母線電壓稱作邊界輸出母線電壓，其不屬于該子網絡的聯絡線母線電壓稱作邊界輸入母線電壓。
5.根據權利要求1所述的潮流并行計算方法，其特征在于步驟4)中，電力子網絡潮流方程具有如下特點a)、列寫每個葉結點對應網絡的潮流代數方程，即PQ母線的功率平衡方程和PV結點的有功和電壓平衡方程；b)、所有子網絡潮流方程聯立即是全網潮流的潮流方程；c)、對各子網絡潮流方程進行泰勒(Taylor)級數展開，略去泰勒級數中二階及以上的高階項，得到子系統的雅可比矩陣方程；雅可比矩陣方程由雅可比矩陣、母線電壓的修正量向量和不平衡向量組成，在每次迭代中，雅可比矩陣和不平衡向量都是確定的實數矩陣和向量，需要求解的是母線電壓的修正量向量。
6.根據權利要求1所述的潮流并行計算方法，其特征在于步驟5)中，雅可比矩陣方程變量修正量的計算如下全網母線電壓的修正量向量的計算通過電網潮流計算樹的“上行替換”和“下行計算”操作完成，其中，a)、上行替換操作該上行替換操作從所有葉結點開始，對于葉結點，通過對雅可比矩陣方程進行初等代數變換，將葉結點中的邊界輸出和內部電壓變量的修正量都寫成由邊界輸入電壓變量的修正量表示的線性表達式；對于中間根結點中間根結點要作的上行替換是接受下層結點的邊界母線電壓修正量線性表達式，消去該中間根結點對應的網絡分割的邊界母線電壓修正量，獲得上一層根結點需要的邊界母線電壓的修正量線性表達；對于根結點，聯立其所有下層結點的邊界輸出、輸入電壓修正量表示的線性方程組，并求解該方程組，所得即邊界輸出、輸入電壓修正量；b)、下行計算操作下行計算操作由根節點至潮流計算樹所有葉結點進行，將解得的邊界電壓修正量帶入上行替換所得內部或下一層邊界輸入、輸出電壓修正量和邊界輸入電壓修正量的線性表達式，最后求得各子網絡內部母線電壓修正向量。
7.根據權利要求1所述的潮流并行計算方法，其特征在于步驟6)中，牛頓迭代收斂的判據是所有子網絡母線電壓修正量的絕對值皆小于某一門檻值ε，ε＜10-5。
8.根據權利要求1所述的潮流并行計算方法，其特征在于步驟7)中，大規模電力系統潮流計算方法在集群計算機系統上并行計算的方法是分別在集群計算機系統的子計算機上實現葉結點的上行替代和下行計算，在主機上實現各層根結點的上行替代和下行計算；子計算機上的計算并行同步進行。
全文摘要
本發明提供一種大規模電力系統潮流計算方法，其步驟包括將電力網絡按照選定聯絡線分割為適度規模的多個子網，子網按其網絡聯結方式組成潮流計算樹；根據網絡劃分的結果，電力網絡每一子網絡對應潮流計算樹的一個葉結點；第一層網絡分割選定的聯絡線對應潮流計算樹的根結點；第一層以后的分割選定的聯絡線對應潮流計算樹的中間根結點。對各子網絡分別建立潮流方程，并采用牛頓迭代法求解上述非線性代數方程，求解過程中變量分別為內部、內部邊界和外部邊界的母線電壓，計算通過“上行替換”和“下行計算”操作完成。該計算方法適用于大規模電力系統并行計算實現，計算量小，編程靈活，為大規模電力系統的在線潮流計算提供了新的實現方法。
文檔編號G06F17/00GK1929234SQ200610015718
公開日2007年3月14日 申請日期2006年9月22日 優先權日2006年9月22日
發明者房大中 申請人:天津大學