TWI421705B - Distribution power flow calculating system and method - Google Patents
Distribution power flow calculating system and method Download PDFInfo
- Publication number
- TWI421705B TWI421705B TW98103895A TW98103895A TWI421705B TW I421705 B TWI421705 B TW I421705B TW 98103895 A TW98103895 A TW 98103895A TW 98103895 A TW98103895 A TW 98103895A TW I421705 B TWI421705 B TW I421705B
- Authority
- TW
- Taiwan
- Prior art keywords
- matrix
- row
- column
- relationship matrix
- relationship
- Prior art date
Links
- 238000000034 method Methods 0.000 title claims description 26
- 239000011159 matrix material Substances 0.000 claims description 205
- 238000005206 flow analysis Methods 0.000 claims description 33
- 238000002347 injection Methods 0.000 claims description 21
- 239000007924 injection Substances 0.000 claims description 21
- 239000006185 dispersion Substances 0.000 claims 1
- 238000010586 diagram Methods 0.000 description 14
- 238000007796 conventional method Methods 0.000 description 9
- 238000004364 calculation method Methods 0.000 description 4
- 238000002940 Newton-Raphson method Methods 0.000 description 1
- 241001604129 Polydactylus Species 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 230000006870 function Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000000750 progressive effect Effects 0.000 description 1
- 238000004088 simulation Methods 0.000 description 1
Classifications
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02E—REDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
- Y02E40/00—Technologies for an efficient electrical power generation, transmission or distribution
- Y02E40/30—Reactive power compensation
Landscapes
- Supply And Distribution Of Alternating Current (AREA)
Description
本發明係關於一種分散式電力潮流分析系統及方法。The present invention relates to a decentralized power flow analysis system and method.
參考圖1，其顯示習知匯流排之三相阻抗模型示意圖。其為匯流排0 至匯流排k 的三相阻抗模型，以Z 矩陣表示為一4×4矩陣如(1)式。Referring to Figure 1, there is shown a schematic diagram of a three-phase impedance model of a conventional bus. It is a three-phase impedance model of bus bar 0 to bus bar k , represented by a Z matrix as a 4×4 matrix such as (1).
經由卡農降階法(Kron's Reduction)可降成3×3矩陣如(2)式。It can be reduced to a 3 × 3 matrix such as (2) by Kron's Reduction.
因此匯流排0 至匯流排k 的電壓與電流的關係式可透過Z 矩陣表示成(3)式Therefore, the relationship between the voltage and current of bus bar 0 to bus bar k can be expressed as (3) through the Z matrix.
上式以電壓差量ΔV可表示成一般式，如(4)式The above formula can be expressed as a general formula with a voltage difference ΔV, such as (4)
[ΔV abc ]=[Z abc ][I abc ] (4)[Δ V abc ]=[ Z abc ][ I abc ] (4)
習知分散式電力潮流分析系統及方法必須利用複雜的計算，且於新增匯流排或阻抗時，阻抗(Z)矩陣之變動相當大，計算亦相當繁複，造成執行速度降低及需要大記憶體空間，且強健性亦不佳。Conventional distributed power flow analysis systems and methods must utilize complex calculations, and when new busbars or impedances are added, the impedance (Z) matrix varies considerably, and the calculation is quite complicated, resulting in reduced execution speed and large memory. Space, and the robustness is not good.
因此，有必要提供一種創新且具進步性的分散式電力潮流分析系統及方法，以解決上述問題。Therefore, it is necessary to provide an innovative and progressive distributed power flow analysis system and method to solve the above problems.
本發明提供一種分散式電力潮流分析系統，包括：一第一關係矩陣建立裝置、一第二關係矩陣建立裝置及一分散式電力潮流分析裝置。該第一關係矩陣建立裝置用以建立一第一關係矩陣，其為匯流排注入電流矩陣與支路電流矩陣之關係，其中匯流排注入電流矩陣係為複數個匯流排之注入電流，支路電流矩陣係為複數個匯流排間之電流。該第二關係矩陣建立裝置用以建立一第二關係矩陣，其為匯流排電壓差矩陣與支路電流矩陣之關係，其中匯流排電壓差矩陣係為一參考匯流排與其他匯流排之電壓差。該分散式電力潮流分析裝置用以依據該第一關係矩陣及該第二關係矩陣，分析分散式電力潮流。The invention provides a distributed power flow analysis system, comprising: a first relationship matrix establishing device, a second relationship matrix establishing device and a distributed power flow analyzing device. The first relationship matrix establishing device is configured to establish a first relationship matrix, which is a relationship between a bus matrix injection current matrix and a branch current matrix, wherein the bus bar injection current matrix is an injection current of a plurality of bus bars, and a branch current The matrix is the current between a plurality of bus bars. The second relationship matrix establishing device is configured to establish a second relationship matrix, which is a relationship between a bus voltage difference matrix and a branch current matrix, wherein the busbar voltage difference matrix is a voltage difference between a reference bus and other busbars. . The distributed power flow analysis device is configured to analyze the distributed power flow according to the first relationship matrix and the second relationship matrix.
本發明另提供一種分散式電力潮流分析方法，包括以下步驟：(a)建立一第一關係矩陣，其為匯流排注入電流矩陣與支路電流矩陣之關係，其中匯流排注入電流矩陣係為複數個匯流排之注入電流，支路電流矩陣係為複數個匯流排間之電流；(b)建立一第二關係矩陣，其為匯流排電壓差矩陣與支路電流矩陣之關係，其中匯流排電壓差矩陣係為一參考匯流排與其他匯流排之電壓差；(c)依據該第一關係矩陣及該第二關係矩陣，分析分散式電力潮流。The invention further provides a distributed power flow analysis method, comprising the following steps: (a) establishing a first relationship matrix, which is a relationship between a bus current injection matrix and a branch current matrix, wherein the bus bar injection current matrix is a plurality The injection current of the busbars, the branch current matrix is the current between the plurality of busbars; (b) establishing a second relationship matrix, which is the relationship between the busbar voltage difference matrix and the branch current matrix, wherein the busbar voltage The difference matrix is the voltage difference between a reference bus and other bus bars; (c) analyzing the distributed power flow according to the first relationship matrix and the second relationship matrix.
本發明之系統及方法利用該第一關係矩陣及該第二關係矩陣，可分析分散式電力潮流，且可應用於新增匯流排或阻抗之情形，與其他習知方法比較，本發明之系統及方法在針對分散式電力系統之電力潮流計算具備良好的強健性及執行速度且所需之記憶體空間較少。The system and method of the present invention can analyze the distributed power flow by using the first relationship matrix and the second relationship matrix, and can be applied to the situation of adding a bus bar or impedance, compared with other conventional methods, the system of the present invention The method has good robustness and execution speed for the power flow calculation for the distributed power system and requires less memory space.
參考圖2，其顯示本發明實施例具5匯流排之分散式電力系統之示意圖。參考圖7，其顯示本發明分散式電力潮流分析方法之流程示意圖。參考圖11，其顯示本發明分散式電力潮流分析系統之電路示意圖。配合參考圖2、圖7及圖11，說明本發明之系統及方法。本發明之分散式電力潮流分析系統11包括：一第一關係矩陣建立裝置110、一第二關係矩陣建立裝置120及一分散式電力潮流分析裝置150。配合參考步驟S71，該第一關係矩陣建立裝置110用以建立一第一關係矩陣，其為匯流排注入電流矩陣與支路電流矩陣之關係，其中匯流排注入電流矩陣係為複數個匯流排之注入電流，支路電流矩陣係為複數個匯流排間之電流。Referring to Figure 2, there is shown a schematic diagram of a decentralized power system having five bus bars in accordance with an embodiment of the present invention. Referring to Figure 7, there is shown a flow chart of the distributed power flow analysis method of the present invention. Referring to Figure 11, there is shown a circuit diagram of a distributed power flow analysis system of the present invention. The system and method of the present invention will be described with reference to Figs. 2, 7 and 11. The distributed power flow analysis system 11 of the present invention includes a first relationship matrix establishing device 110, a second relationship matrix establishing device 120, and a distributed power flow analyzing device 150. Referring to the reference step S71, the first relationship matrix establishing means 110 is configured to establish a first relationship matrix, which is a relationship between the bus matrix injection current matrix and the branch current matrix, wherein the bus bar injection current matrix is a plurality of bus bars. Injecting current, the branch current matrix is the current between a plurality of busbars.
利用等效電流注入法，匯流排i 的電流可以其第k次疊代的注入功率及電壓來表示其等效電流值，如式(5)所示With the equivalent current injection method, the current of the bus i can express its equivalent current value by the injection power and voltage of the kth iteration, as shown in equation (5).
其中，V i k 及I i k 為匯流排i 第k次疊代的電壓及電流值。Where V i k and I i k are the voltage and current values of the kth iteration of the bus i .
本發明實施例圖2之匯流排注入電流I 與支路電流B 之關係如式(6)所示。The relationship between the bus current injection current I and the branch current B of FIG. 2 in the embodiment of the present invention is as shown in the formula (6).
因此，支路電流B 與注入電流I 的關係式可透過第一關係矩陣B I 表示如(7a)式。Therefore, the relationship between the branch current B and the injection current I can be expressed as the equation (7a) through the first relation matrix B I .
(7a)式可改表示成一般式，如(7b)式(7a) can be expressed as a general formula, such as (7b)
[B ]=[B I ][I ] (7b)[ B ]=[ B I ][ I ] (7b)
其中，匯流排注入電流矩陣[I]係為複數個匯流排之注入電流，支路電流矩陣[B]係為複數個匯流排間之電流。因此，[B I ]為匯流排注入電流矩陣[I]與支路電流矩陣[B]的第一關係矩陣且為一僅包含0與1的上三角矩陣。The bus current injection current matrix [I] is an injection current of a plurality of bus bars, and the branch current matrix [B] is a current between a plurality of bus bars. Therefore, [ B I ] is the first relation matrix in which the bus matrix injection current matrix [I] and the branch current matrix [B] are injected and is an upper triangular matrix containing only 0 and 1.
配合參考步驟S72，該第二關係矩陣建立裝置120用以建立一第二關係矩陣，其為匯流排電壓差矩陣與支路電流矩陣之關係，其中匯流排電壓差矩陣係為一參考匯流排與其他匯流排之電壓差。Referring to the reference step S72, the second relationship matrix establishing means 120 is configured to establish a second relationship matrix, which is a relationship between the bus voltage difference matrix and the branch current matrix, wherein the busbar voltage difference matrix is a reference bus and The voltage difference between other busbars.
再參考圖2，匯流排1、2、3之電壓可改寫如(8)式。Referring again to Figure 2, the voltages of the busbars 1, 2, 3 can be rewritten as in equation (8).
其中V i 為匯流排i 的電壓，Z ij 為匯流排i,j 間的阻抗。上式(8)，V 3 可改寫如(9)式。 I wherein V i is the voltage of the bus, Z ij is the impedance between the busbar i, j. In the above formula (8), V 3 can be rewritten as in the formula (9).
V 3 =V 0 -B 1 Z 1 -B 2 Z 12 -B 3 Z 23 (9) V 3 = V 0 - B 1 Z 1 - B 2 Z 12 - B 3 Z 23 (9)
由(9)式，匯流排的電壓差可表示成支路電流、支路阻抗的函數。因此匯流排電壓差矩陣可表示如(10a)式。From equation (9), the voltage difference of the busbar can be expressed as a function of the branch current and the branch impedance. Therefore, the bus voltage difference matrix can be expressed as in the equation (10a).
(10a)式之電壓差量ΔV可表示成一般式，如(10b)式。The voltage difference ΔV of the formula (10a) can be expressed as a general formula such as the formula (10b).
[ΔV ]=[Z V-BC ][B] (10b)[Δ V ]=[ Z V-BC ][B] (10b)
其中，匯流排電壓差矩陣[ΔV]係為一參考匯流排電壓V 0 與其他匯流排之電壓差，[Z V-BC ]為匯流排電壓差矩陣[ΔV]與支路電流矩陣[B]的第二關係矩陣。該第二關係矩陣係為匯流排間阻抗，且係為一下三角矩陣。Wherein, the bus voltage difference matrix [ΔV] is a voltage difference between a reference bus voltage V 0 and other bus bars, [ Z V-BC ] is a bus bar voltage difference matrix [ΔV] and a branch current matrix [B] The second relational matrix. The second relational matrix is the busbar impedance and is a lower triangular matrix.
配合參考步驟S75，該分散式電力潮流分析裝置150用以依據該第一關係矩陣[B I ]及該第二關係矩陣[Z V-BC ]，分析分散式電力潮流。因此，將(7b)式代入(10b)式可得(11)式。Referring to step S75, the distributed power flow analysis device 150 is configured to analyze the distributed power flow according to the first relationship matrix [ B I ] and the second relationship matrix [ Z V-BC ]. Therefore, substituting the formula (7b) into the formula (10b) gives the formula (11).
[ΔV ]=[Z V-BC ][B I ][I ]=[ZDPF ][I ] (11)[Δ V ]=[ Z V-BC ][ B I ][ I ]=[Z DPF ][ I ] (11)
[ΔV k+l ]=[ZDPF ][I k ] (12)[Δ V k+l ]=[Z DPF ][ I k ] (12)
求解(5)式及(12)式，可分析分散式電力潮流。Solve the (5) and (12) equations and analyze the distributed power flow.
參考圖3，其顯示本發明第一新增匯流排實施例之分散式電力系統之示意圖。配合參考圖2及3，本發明第一新增匯流排實施例係增加一新匯流排bus 6連接至該參考匯流排bus 0，且一新阻抗Z new 於該新匯流排bus 6與該參考匯流排bus 0之間。Referring to Figure 3, there is shown a schematic diagram of a decentralized power system of the first new busbar embodiment of the present invention. Referring to Figures 2 and 3, the first new busbar embodiment of the present invention adds a new bus bus 6 to the reference bus bus 0, and a new impedance Z new to the new bus bus 6 and the reference. Bus between bus 0.
配合參考圖3、圖7及圖11，參考步驟S73及S74，判斷是否新增匯流排或阻抗，若有則修正第一關係矩陣及第二關係矩陣。本發明之分散式電力潮流分析系統11另包括一第一關係矩陣修正裝置130及一第二關係矩陣修正裝置140，用以於新增匯流排或阻抗時，修正該第一關係矩陣及該第二關係矩陣。依據上述第一新增匯流排實施例，該第一關係矩陣修正裝置130將該第一關係矩陣新增一行及一列，在該新增行列之對角線位置為1，其餘為0。表示如(13a)式。Referring to FIG. 3, FIG. 7, and FIG. 11, referring to steps S73 and S74, it is determined whether or not a bus bar or impedance is added, and if so, the first relationship matrix and the second relationship matrix are corrected. The distributed power flow analysis system 11 of the present invention further includes a first relationship matrix correction device 130 and a second relationship matrix correction device 140 for modifying the first relationship matrix and the first when adding a bus bar or impedance Two relational matrices. According to the first new bus bar embodiment, the first relationship matrix correcting device 130 adds a row and a column to the first relationship matrix, and the diagonal position of the newly added row is 1 and the rest is 0. Expressed as in (13a).
該第二關係矩陣修正裝置140將該第二關係矩陣新增一行及一列，在該新增行列之對角線位置為該新阻抗，其餘為0。表示如(13b)式。The second relationship matrix correcting device 140 adds a row and a column to the second relationship matrix, and the diagonal position of the newly added row and column is the new impedance, and the rest is 0. Expressed as in (13b).
(13b)式可表示成一般式，如(13c)式。The formula (13b) can be expressed as a general formula such as the formula (13c).
將(13a)式代入(13c)式可得(13d)式。Substituting the formula (13a) into the formula (13c) gives the formula (13d).
參考圖4，其顯示本發明第二新增匯流排實施例之分散式電力系統之示意圖。配合參考圖2及4，本發明第二新增匯流排實施例係增加一新匯流排bus 6連接至一第k個匯流排(在本實施例為第5個匯流排bus 5)，且一新阻抗Z new 於該新匯流排bus 6與該第5個匯流排bus 5之間。Referring to Figure 4, there is shown a schematic diagram of a distributed power system of a second new busbar embodiment of the present invention. With reference to FIGS. 2 and 4, the second new bus bar embodiment of the present invention adds a new bus bar 6 to a kth bus (in the present embodiment, the fifth bus bus 5), and one A new impedance Z new is between the new bus bus 6 and the fifth bus bus 5.
配合參考圖4及圖11，該第一關係矩陣修正裝置130將該第一關係矩陣新增一行及一列，並複製在該第一關係矩陣之第k行(在本實施例為第5行)之值至該新增行，該新增行列之對角線位置為1，其餘為0。表示如(14a)式。Referring to FIG. 4 and FIG. 11, the first relationship matrix correcting device 130 adds a row and a column to the first relationship matrix, and copies the kth row of the first relationship matrix (in the fifth row of the embodiment) The value is up to the new row. The new row has a diagonal position of 1 and the rest is 0. Expressed as in (14a).
其中，col.(k) 為原第一關係矩陣之第k行(在本實施例為第5行)之值。Where col.(k) is the value of the kth row (the fifth row in this embodiment ) of the original first relationship matrix.
該第二關係矩陣修正裝置140將該第二關係矩陣新增一行及一列，並複製在該第二關係矩陣之第k列(在本實施例為第5列)之值至該新增列，在該新增行列之對角線位置為該新阻抗，其餘為0。表示如(14b)式。The second relationship matrix correcting means 140 adds a row and a column to the second relationship matrix, and copies the value of the kth column (the fifth column in this embodiment) to the newly added column. The diagonal position of the newly added row is the new impedance, and the rest is 0. Expressed as in (14b).
(14b)式可表示成一般式，如(14c)式。The formula (14b) can be expressed as a general formula such as the formula (14c).
其中，row.(k) 為原第二關係矩陣之第k列(在本實施例為第5列)之值。將(14a)式代入(14c)式可得(14d)式。Where row.(k) is the value of the kth column (the fifth column in this embodiment ) of the original second relationship matrix. Substituting the formula (14a) into the formula (14c) gives the formula (14d).
參考圖5，其顯示本發明第三新增阻抗實施例之分散式電力系統之示意圖。配合參考圖2及5，本發明第三新增阻抗實施例係增加一新阻抗Znew 於一第i個匯流排(在本實施例為第4個匯流排bus 4)與一第j個匯流排(在本實施例為第5個匯流排bus 5)之間。Referring to Figure 5, there is shown a schematic diagram of a distributed power system of a third new impedance embodiment of the present invention. Referring to Figures 2 and 5, the third new impedance embodiment of the present invention adds a new impedance Z new to an i-th bus bar (in the present embodiment, the fourth bus bar 4) and a j-th confluence. The row (in this embodiment is the fifth bus bus 5).
配合參考圖6及圖11，該第一關係矩陣修正裝置130將該第一關係矩陣新增一行及一列，並將該第一關係矩陣之第i行(在本實施例為第4行)之值減第j行(在本實施例為第5行)之值置於該新增行，該新增行列之對角線位置為1，其餘為0。在本實施例中係增加一新阻抗Znew 於第4個匯流排bus4與第5個匯流排bus 5之間，在第4個匯流排bus 4與第5個匯流排bus 5之注入電流表示如(15a)式。Referring to FIG. 6 and FIG. 11, the first relationship matrix correcting device 130 adds a row and a column to the first relationship matrix, and the ith row of the first relationship matrix (the fourth row in this embodiment) The value minus the jth row (in the fifth row of this embodiment) is placed in the newly added row, the diagonal position of the newly added row and column is 1, and the rest is 0. In this embodiment, a new impedance Z new is added between the fourth bus bar bus 4 and the fifth bus bar 5, and the injection currents in the fourth bus bar 4 and the fifth bus bus 5 are expressed. Such as (15a).
第一關係矩陣修正如(15a)式。The first relationship matrix is modified as in (15a).
將B new 移項，則(15b)式可表示如(15c)式。If B new is shifted, then (15b) can be expressed as (15c).
(15c)式可表示成一般式，如(15d)式。The formula (15c) can be expressed as a general formula such as (15d).
其中，col.(i-j) 為原第一關係矩陣之第i行(在本實施例為第4行)之值減第j行(在本實施例為第5行)之值。Where col.(ij) is the value of the i-th row (the fourth row in the present embodiment ) of the original first relationship matrix minus the value of the jth row (the fifth row in this embodiment).
該第二關係矩陣修正裝置140將該第二關係矩陣新增一行及一列，並將該第一關係矩陣之第i列之值減第j列之值置於該新增列，在該新增行列之對角線位置為該新阻抗，其餘為0。應用KVL於該新增回路，表示如(16a)式。The second relationship matrix correcting device 140 adds a row and a column to the second relationship matrix, and places the value of the i-th column of the first relationship matrix minus the value of the j-th column in the newly added column. The diagonal position of the row and column is the new impedance, and the rest is 0. Apply KVL to the new loop, which is expressed as (16a).
Z 23 B 3 +Z 34 B 4 +Z new B new -Z 25 B 5 =0 (16a) Z 23 B 3 + Z 34 B 4 + Z new B new - Z 25 B 5 =0 (16a)
合併(16a)式及(10a)式，匯流排電壓差矩陣可表示如(16b)式。Combining the equations (16a) and (10a), the bus voltage difference matrix can be expressed as (16b).
(16b)式可表示成一般式，如(16c)式。The formula (16b) can be expressed as a general formula such as the formula (16c).
其中，row.(i-j) 為原第二關係矩陣之第i列(在本實施例為第4列)之值減第j列(在本實施例為第5列)之值。將(15c)式代入(16c)式可得(16d)式。Wherein, row.(ij) is the value of the i-th column (the fourth column in the present embodiment ) of the original second relationship matrix minus the value of the j-th column (the fifth column in the present embodiment). Substituting the formula (15c) into the formula (16c) gives the formula (16d).
參考圖6，其顯示本發明第四新增匯流排及阻抗實施例之分散式電力系統之示意圖。配合參考圖2及6，本發明第四新增匯流排及阻抗實施例係增加一新匯流排bus 6於一第i個匯流排(在本實施例為第2個匯流排bus 2)與一第j個匯流排之間(在本實施例為第5個匯流排bus 5)，且一新阻抗Znew 於該新匯流排bus 6與該第i個匯流排(在本實施例為第2個匯流排bus 2)之間。Referring to Figure 6, there is shown a schematic diagram of a fourth new busbar of the present invention and a distributed power system of an impedance embodiment. With reference to FIGS. 2 and 6, the fourth new bus bar and impedance embodiment of the present invention adds a new bus bar 6 to an i-th bus bar (in this embodiment, a second bus bar 2) and a Between the jth busbars (in the present embodiment, the fifth busbar bus 5), and a new impedance Z new in the new busbar bus 6 and the i-th busbar (in this embodiment, the second Bus between bus 2).
配合參考圖5及圖11，該第一關係矩陣修正裝置130將該第一關係矩陣新增二行及二列，並複製在該第一關係矩陣之第i行(在本實施例為第2行)之值至該新增第一行，且將該第一關係矩陣之新增第一行之值減第j行(在本實施例為第5行)之值置於該新增第二行，該新增行列之對角線位置為1，其餘為0，表示如(17a)式。Referring to FIG. 5 and FIG. 11, the first relationship matrix correcting device 130 adds two rows and two columns to the first relationship matrix, and copies the i-th row of the first relationship matrix (in this embodiment, the second row) The value of the row is added to the first row, and the value of the first row of the first relationship matrix is reduced by the value of the jth row (the fifth row in this embodiment). Row, the new row and column has a diagonal position of 1, and the rest is 0, indicating that it is of (17a).
(17a)式可表示成一般式，如(17b)式。The formula (17a) can be expressed as a general formula such as the formula (17b).
其中，col.(i) 為原第一關係矩陣之第i行(在本實施例為第2行)之值；col.(k-j) 為第一關係矩陣之新增第一行之值減第j行(在本實施例為第5行)之值。Where col.(i) is the value of the ith row of the original first relationship matrix (the second row in this embodiment); col.(kj) is the value of the first row added by the first relationship matrix minus The value of line j (in the fifth line of this embodiment).
該第二關係矩陣修正裝置140將該第二關係矩陣新增二行及二列，並複製在該第二關係矩陣之第i列(在本實施例為第2列)之值至該新增第一列，且將該第二關係矩陣之新增第一列之值減第j列(在本實施例為第5列)之值置於該新增第二列，在該新增第一行及新增第一列之對角線位置為該新阻抗，其餘為0。應用KVL於該新增回路，表示如(18a)式。The second relationship matrix correcting means 140 adds two rows and two columns to the second relationship matrix, and copies the value of the i-th column (the second column in this embodiment) of the second relationship matrix to the new one. The first column, and the value of the first column added by the second relationship matrix minus the value of the jth column (the fifth column in this embodiment) is placed in the newly added second column, and the first one is added. The diagonal position of the row and the first column added is the new impedance, and the rest is 0. Apply KVL to the new loop, which is expressed as (18a).
Z new B new -B 5 Z 25 =0 (18a) Z new B new - B 5 Z 25 =0 (18a)
合併(18a)式及(10a)式，匯流排電壓差矩陣可表示如(18b)式。Combining equations (18a) and (10a), the busbar voltage difference matrix can be expressed as in equation (18b).
(18b)式可表示成一般式，如(18c)式。The formula (18b) can be expressed as a general formula such as the formula (18c).
其中，row.(i) 為原第二關係矩陣之第i列(在本實施例為第2列)之值；row.(k-j) 為第二關係矩陣之新增第一列之值之值減第j列(在本實施例為第5列)之值。將(17b)式代入(18c)式可得(18d)式。Where row.(i) is the value of the i-th column of the original second relationship matrix (in the second column in this embodiment); row.(kj) is the value of the value of the first column added to the second relationship matrix. Subtract the value of column j (in the fifth column of this embodiment). Substituting the formula (17b) into the formula (18c) gives the formula (18d).
(18d)式可表示成如(19)式。The formula (18d) can be expressed as in the formula (19).
利用卡農降階法(Kron's Reduction)，匯流排電壓差矩陣可表示如(20)式。With Kron's Reduction, the busbar voltage difference matrix can be expressed as equation (20).
ΔV =[A -M T N -1 M ][I ] (20)Δ V =[ A - M T N -1 M ][ I ] (20)
參考圖8，其顯示本發明模擬之分散式電力系統示意圖。該模擬之分散式電力系統具有8個匯流排，本發明係以習知之高斯Z矩陣法(Gauss implicit Z-matrix)、習知牛頓拉福生方法(Newton-Raphson)及本發明之方法進行模擬比較。Referring to Figure 8, there is shown a schematic diagram of a simulated distributed power system of the present invention. The simulated distributed power system has eight bus bars, and the present invention performs simulation comparison by the conventional Gaussian Z-matrix, the Newton-Raphson method and the method of the present invention. .
參考圖9，其顯示本發明與習知方法之執行時間比較圖。其中，第一習知方法為習知之高斯Z矩陣法，第二習知方法為牛頓拉福生方法。由圖9可知，本發明之執行時間(normalized execution time,NET)遠小於習知方法，且當匯流排數越多時，本發明之方法更有效率且更快。Referring to Figure 9, there is shown a comparison of execution times of the present invention with conventional methods. Among them, the first conventional method is the conventional Gaussian Z matrix method, and the second conventional method is the Newton Rafson method. As can be seen from Fig. 9, the normalized execution time (NET) of the present invention is much smaller than the conventional method, and the method of the present invention is more efficient and faster when the number of bus bars is larger.
參考圖10，其顯示本發明與習知方法之疊代次數比較圖。在不同之KR 比值下，本發明之疊代次數較習知快速解偶法(Fast Decoupled Load Flow Method)為低。Referring to Figure 10, there is shown a comparison of the number of iterations of the present invention with conventional methods. At different K R ratios, the number of iterations of the present invention is lower than the Fast Decoupled Load Flow Method.
本發明之系統及方法利用該第一關係矩陣及該第二關係矩陣，可分析分散式電力潮流，且可應用於新增匯流排或阻抗之情形，與其他習知方法比較，本發明之系統及方法在針對分散式電力系統之電力潮流計算具備良好的強健性及執行速度且所需之記憶體空間較少。The system and method of the present invention can analyze the distributed power flow by using the first relationship matrix and the second relationship matrix, and can be applied to the situation of adding a bus bar or impedance, compared with other conventional methods, the system of the present invention The method has good robustness and execution speed for the power flow calculation for the distributed power system and requires less memory space.
惟上述實施例僅為說明本發明之原理及其功效，而非限制本發明。因此，習於此技術之人士對上述實施例進行修改及變化仍不脫本發明之精神。本發明之權利範圍應如後述之申請專利範圍所列。However, the above embodiments are merely illustrative of the principles of the invention and its effects, and are not intended to limit the invention. Therefore, those skilled in the art can make modifications and changes to the above embodiments without departing from the spirit of the invention. The scope of the invention should be as set forth in the appended claims.
11．．．本發明之分散式電力潮流分析系統11. . . Decentralized power flow analysis system of the present invention
110．．．第一關係矩陣建立裝置110. . . First relationship matrix establishing device
120．．．第二關係矩陣建立裝置120. . . Second relationship matrix establishing device
130．．．第一關係矩陣修正裝置130. . . First relationship matrix correction device
140．．．第二關係矩陣修正裝置140. . . Second relationship matrix correction device
150．．．分散式電力潮流分析裝置150. . . Decentralized power flow analysis device
圖1顯示習知匯流排之三相阻抗模型示意圖；Figure 1 shows a schematic diagram of a three-phase impedance model of a conventional bus bar;
圖2顯示本發明實施例具5匯流排之分散式電力系統之示意圖；2 is a schematic view showing a distributed power system with 5 bus bars according to an embodiment of the present invention;
圖3顯示本發明第一新增匯流排實施例之分散式電力系統之示意圖；3 is a schematic diagram showing a distributed power system of the first new bus bar embodiment of the present invention;
圖4顯示本發明第二新增匯流排實施例之分散式電力系統之示意圖；4 is a schematic diagram showing a distributed power system of a second new bus bar embodiment of the present invention;
圖5顯示本發明第三新增阻抗實施例之分散式電力系統之示意圖；5 is a schematic diagram showing a distributed power system of a third new impedance embodiment of the present invention;
圖6顯示本發明第四新增匯流排及阻抗實施例之分散式電力系統之示意圖；6 is a schematic diagram showing a fourth power bus of the present invention and a distributed power system of an impedance embodiment;
圖7顯示本發明分散式電力潮流分析方法之流程示意圖；7 is a schematic flow chart showing a method for analyzing a distributed power flow according to the present invention;
圖8顯示本發明模擬之分散式電力系統示意圖；Figure 8 is a schematic view showing the simulated distributed power system of the present invention;
圖9顯示本發明與習知方法之執行時間比較圖；Figure 9 is a graph showing the comparison of the execution time of the present invention with the conventional method;
圖10顯示本發明與習知方法之疊代次數比較圖；及Figure 10 is a graph showing the comparison of the number of iterations of the present invention with the conventional method;
圖11顯示本發明分散式電力潮流分析系統之電路示意圖。Figure 11 is a circuit diagram showing the distributed power flow analysis system of the present invention.
11．．．本發明之分散式電力潮流分析系統11. . . Decentralized power flow analysis system of the present invention
110．．．第一關係矩陣建立裝置110. . . First relationship matrix establishing device
120．．．第二關係矩陣建立裝置120. . . Second relationship matrix establishing device
130．．．第一關係矩陣修正裝置130. . . First relationship matrix correction device
140．．．第二關係矩陣修正裝置140. . . Second relationship matrix correction device
150．．．分散式電力潮流分析裝置150. . . Decentralized power flow analysis device
Claims (16)
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| TW98103895A TWI421705B (en) | 2009-02-06 | 2009-02-06 | Distribution power flow calculating system and method |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| TW98103895A TWI421705B (en) | 2009-02-06 | 2009-02-06 | Distribution power flow calculating system and method |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| TW201030534A TW201030534A (en) | 2010-08-16 |
| TWI421705B true TWI421705B (en) | 2014-01-01 |
Family
ID=44854285
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| TW98103895A TWI421705B (en) | 2009-02-06 | 2009-02-06 | Distribution power flow calculating system and method |
Country Status (1)
| Country | Link |
|---|---|
| TW (1) | TWI421705B (en) |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US5734586A (en) * | 1995-05-05 | 1998-03-31 | Cornell Research Foundation, Inc. | System for achieving optimal steady state in power distribution networks |
| TWI264864B (en) * | 2005-04-08 | 2006-10-21 | Univ Chang Gung | Power flow calculation method of power grid with unified power flow controller |
| TWI269137B (en) * | 2005-03-18 | 2006-12-21 | Univ Chang Gung | Power flow models of interline power flow controllers (IPFC) for large-scale power systems |
| TWI270762B (en) * | 2005-03-25 | 2007-01-11 | Univ Chang Gung | Control method for enhancing dynamic model of power flow |
-
2009
- 2009-02-06 TW TW98103895A patent/TWI421705B/en not_active IP Right Cessation
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US5734586A (en) * | 1995-05-05 | 1998-03-31 | Cornell Research Foundation, Inc. | System for achieving optimal steady state in power distribution networks |
| TWI269137B (en) * | 2005-03-18 | 2006-12-21 | Univ Chang Gung | Power flow models of interline power flow controllers (IPFC) for large-scale power systems |
| TWI270762B (en) * | 2005-03-25 | 2007-01-11 | Univ Chang Gung | Control method for enhancing dynamic model of power flow |
| TWI264864B (en) * | 2005-04-08 | 2006-10-21 | Univ Chang Gung | Power flow calculation method of power grid with unified power flow controller |
Also Published As
| Publication number | Publication date |
|---|---|
| TW201030534A (en) | 2010-08-16 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| JP5426716B2 (en) | Distribution network power flow analysis system and method | |
| Elsaiah et al. | Analytical approach for placement and sizing of distributed generation on distribution systems | |
| Gözel et al. | An analytical method for the sizing and siting of distributed generators in radial systems | |
| Dufour et al. | A combined state-space nodal method for the simulation of power system transients | |
| KR102066894B1 (en) | Calculation method and program of electric potential, electric current, and peripheral electromagnetic field in electric circuit | |
| CN102156760A (en) | Saber-based circuit failure simulation analyzing method | |
| Wu et al. | A study of the impacts of flow direction and electrical constraints on vulnerability assessment of power grid using electrical betweenness measures | |
| Yang et al. | Decomposed Newton algorithm-based three-phase power-flow for unbalanced radial distribution networks with distributed energy resources and electric vehicle demands | |
| Lisboa et al. | A fast power flow method for radial networks with linear storage and no matrix inversions | |
| CN101291062B (en) | A Method for Analyzing DC Power Flow in Dynamic Topological Power Grid | |
| CN104734148B (en) | Three-phrase power-distributing network continuation power flow analysis of distributed power supply | |
| CN115982914B (en) | Method and system for evaluating grounding grid corrosion based on elastic grid theory | |
| Yao et al. | Contingency analysis based on partitioned and parallel holomorphic embedding | |
| CN111682546B (en) | DC power flow improvement algorithm based on sensitivity analysis | |
| Moura et al. | IMICV fault analysis method with multiple PV grid-connected inverters for distribution systems | |
| CN104881566B (en) | A kind of Three-Phase Transformer modeling method for power network asymmetry Load flow calculation | |
| TWI421705B (en) | Distribution power flow calculating system and method | |
| CN104142419B (en) | Power grid short-circuit current acquiring method with influences of loads considered | |
| CN118278323A (en) | Electric arc ground fault simulation method, system and computer equipment for transmission and distribution network | |
| Alsaadi et al. | An effective approach for distribution system power flow solution | |
| Liu et al. | Development of phase domain frequency-dependent transmission line model on FPGA for real-time digital simulator | |
| WO2018188393A1 (en) | Linearized determination method and system for power flow of flexible direct-current power grid, and computer storage medium | |
| CN114441835B (en) | Dual-conductor current measurement method and device based on magnetic field sensing chip | |
| Wang et al. | A waveform relaxation technique for steady state initialization of circuits with nonlinear elements and ideal diodes | |
| Li et al. | Whole‐line fault analysis method for unbalanced distribution networks with inverter interfaced distributed generators |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| MM4A | Annulment or lapse of patent due to non-payment of fees |