CN112398134A - Transmission and distribution cooperative rapid robust load flow calculation method - Google Patents
Transmission and distribution cooperative rapid robust load flow calculation method Download PDFInfo
- Publication number
- CN112398134A CN112398134A CN202011215316.4A CN202011215316A CN112398134A CN 112398134 A CN112398134 A CN 112398134A CN 202011215316 A CN202011215316 A CN 202011215316A CN 112398134 A CN112398134 A CN 112398134A
- Authority
- CN
- China
- Prior art keywords
- power
- transmission
- flow calculation
- voltage
- phase angle
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Pending
Links
Images
Classifications
-
- H—ELECTRICITY
- H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
- H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
- H02J3/00—Circuit arrangements for ac mains or ac distribution networks
- H02J3/04—Circuit arrangements for ac mains or ac distribution networks for connecting networks of the same frequency but supplied from different sources
- H02J3/06—Controlling transfer of power between connected networks; Controlling sharing of load between connected networks
-
- H—ELECTRICITY
- H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
- H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
- H02J2203/00—Indexing scheme relating to details of circuit arrangements for AC mains or AC distribution networks
- H02J2203/10—Power transmission or distribution systems management focussing at grid-level, e.g. load flow analysis, node profile computation, meshed network optimisation, active network management or spinning reserve management
-
- H—ELECTRICITY
- H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
- H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
- H02J2203/00—Indexing scheme relating to details of circuit arrangements for AC mains or AC distribution networks
- H02J2203/20—Simulating, e g planning, reliability check, modelling or computer assisted design [CAD]
Landscapes
- Engineering & Computer Science (AREA)
- Power Engineering (AREA)
- Supply And Distribution Of Alternating Current (AREA)
Abstract
The invention discloses a transmission and distribution cooperative rapid robust power flow calculation method which mainly comprises seven substeps of initialization, power distribution network power flow calculation, power transmission network power flow calculation, convergence inspection, iteration variable correction, current iteration time updating and maximum iteration time inspection. The method of the invention corrects the iteration variable by using the iteration variable values of two adjacent iteration steps, can obviously improve the efficiency and robustness of the overall algorithm of the transmission and distribution cooperative power flow calculation, and can solve the problem of divergence of the traditional master-slave splitting method particularly when the power distribution network has large disturbance on the transmission network, namely the sensitivity of the boundary node injection power on the boundary node voltage amplitude is high. Meanwhile, the method is high in operability and easy to implement.
Description
Technical Field
The invention relates to a power grid load flow analysis method, relates to a global load flow analysis method for power transmission network and power distribution network cooperation, and particularly relates to a rapid robust load flow calculation method for power transmission and distribution cooperation.
Background
In a traditional dispatching mode, the transmission and distribution network is managed by different dispatching departments, and analysis and calculation are also independently performed. However, with the massive access of distributed power supplies, the coupling of the transmission network and the distribution network is greatly enhanced, and the traditional independent analysis and calculation mode lacks synchronism and cannot meet the requirement of calculation accuracy. Therefore, the transmission network and the distribution network need to calculate an accurate global power flow cooperatively. The traditional transmission and distribution cooperative power flow calculation based on the master-slave splitting method can realize accurate calculation of the global power flow when the power distribution network has small disturbance on the power transmission network. However, with the development of an active power distribution network, the disturbance of the power distribution network to the transmission network is large, that is, the sensitivity of the injected power of a root node (boundary node) of the power distribution network with respect to the voltage amplitude of the root node is high, and the convergence performance cannot be ensured by the traditional master-slave splitting method.
Disclosure of Invention
The invention aims to solve the technical problem of perfecting and improving the prior technical scheme and provides a transmission and distribution cooperative fast robust load flow calculation method so as to improve the overall calculation efficiency and convergence.
Therefore, the invention adopts the following technical scheme:
a transmission and distribution cooperative rapid robust power flow calculation method comprises the following steps:
1) initialization: setting initial values of boundary node voltage amplitudes
Initial value of voltage phase angle
Setting convergence accuracy eta; setting a maximum iteration number K; setting the current iteration number k to be 0;
2) and (3) power distribution network load flow calculation: according to the current boundary node voltage amplitude
Phase angle of voltage
Calculating the power flow of the power distribution network to obtain the node voltage amplitude of each node of the power distribution network
Phase angle of voltage
Computing injected power of boundary nodes
And will be
Transmitting to the transmission grid;
3) and (3) power transmission network load flow calculation: injecting power according to the current boundary node
Calculating the power flow of the power transmission network to obtain the node voltage amplitude of each node of the power transmission network
Phase angle of voltage
And new boundary node voltage amplitude
Phase angle of voltage
4) And (3) testing the convergence: if it is not
The algorithm is converged, and the algorithm is terminated, otherwise, the step 5) is carried out;
5) iterative variable correction: if the current iteration times k are odd, directly jumping to the step 6); if the current iteration number k is an even number, updating the value of the iteration variable according to the following formula:
6) updating the current iteration times: let k be k + 1;
7) and (5) checking the maximum iteration times: if K is equal to K, then,the algorithm reaches the maximum iteration times without convergence, and the algorithm is terminated, otherwise the current boundary node voltage amplitude value is obtained
Phase angle of voltage
Transmitting to the distribution network and turning to step 2).
In the invention, the boundary node is a node at the intersection of a power distribution network and a power transmission network; each node of the power distribution network and each node of the power transmission network do not comprise boundary nodes.
The invention has the beneficial effects that:
the method of the invention corrects the iteration variable by using the iteration variable values of two adjacent iteration steps, can obviously improve the efficiency and robustness of the overall algorithm of the transmission and distribution cooperative power flow calculation, and can solve the problem of divergence of the traditional master-slave splitting method particularly when the power distribution network has large disturbance on the transmission network, namely the sensitivity of the boundary node injection power on the boundary node voltage amplitude is high. Meanwhile, the method is high in operability and easy to implement.
Drawings
FIG. 1 is a flow chart of the method of the present invention.
Detailed Description
The technical scheme of the invention is further explained in detail by combining the drawings in the specification.
As shown in fig. 1, the method of the present invention comprises seven main steps: initialization, power distribution network load flow calculation, power transmission network load flow calculation, convergence inspection, iteration variable correction, current iteration times and maximum iteration time inspection. The method comprises the following specific steps:
step S1: setting initial values of boundary node voltage amplitudes
Initial value of voltage phase angle
Setting convergence accuracy eta; setting a maximum iteration number K; and setting the current iteration number k to be 0.
Step S2: according to the current boundary node voltage amplitude
Phase angle of voltage
Calculating the power flow of the power distribution network by using a forward-backward deduction method to obtain the node voltage amplitude of each node of the power distribution network
Phase angle of voltage
Computing injected power of boundary nodes
And will be
And transmitted to the transmission grid.
Step S3: injecting power according to the current boundary node
Calculating the power flow of the power transmission network by using a Newton-Raphson method to obtain the node voltage amplitude of each node of the power transmission network
Phase angle of voltage
And new boundary node voltage amplitude
Phase angle of voltage
Step S4: if it is not
The algorithm ends with convergence, otherwise it proceeds to step S5.
Step S5: if the current iteration times k are odd, directly jumping to the step 6); if the current iteration number k is an even number, updating the value of the iteration variable according to the following formula:
step S6: updating the current iteration times: let k be k + 1.
Step S7: if K is equal to K, the algorithm reaches the maximum iteration times without convergence, and the algorithm is terminated, otherwise, the current boundary node voltage amplitude value is used
Phase angle of voltage
Transmitting to the distribution network and turning to step 2).
The above fast robust power flow calculation method for power transmission and distribution coordination shown in fig. 1 is a specific embodiment of the present invention, and already embodies the substantial features and advances of the present invention, and those skilled in the art can calculate the power distribution network power flow by using different calculation methods in step S2, calculate the power transmission network power flow by using different calculation methods in step S3, or make equivalent modifications to the calculation methods in terms of structure, etc., according to the practical needs of use, and the modifications are within the scope of protection of the present solution.
In order to further embody the substantive characteristics and progress of the invention, the two transmission and distribution calculation examples are utilized to respectively carry out transmission and distribution cooperative load flow calculation under the traditional master-slave split method and the calculation method of the invention, and the convergence times and the calculation time are compared. Example A is formed by splicing an IEEE14 node transmission network and an IEEE69 node distribution network, wherein the distribution network is connected to the transmission network through an ideal transformerNo. 14 nodes, No. 8 nodes, No. 15 nodes and No. 20 nodes of the power distribution network are respectively connected to a PV-type distributed power supply, and the active power is 0.5 MW. The calculation example B is formed by splicing an IEEE14 node power transmission network and an IEEE69 node power distribution network, the power distribution network is connected to a No. 14 node of the power transmission network through an ideal transformer, No. 45 and No. 61 nodes of the power distribution network are respectively connected to a PV type distributed power supply, and the active power is 0.5 MW. The convergence accuracy eta is set to 10-6. The maximum number of iterations K is set to 100. Initial value of boundary node voltage amplitude
Set to 1, initial value of voltage phase angle
Is set to 0.
Table 1 records the number of iterations and the time spent in the calculations for both methods. For example A, the traditional master-slave splitting method needs 20 iterations of convergence and takes 78.8ms, while the method only needs 6 iterations of convergence and takes 27.9ms, which shows that the method can improve the efficiency of the transmission and distribution cooperative power flow calculation. For example B, the traditional master-slave splitting method does not converge, but the method of the invention needs 8 convergence times and takes 37.4ms, which shows that the method of the invention can improve the convergence performance of the transmission and distribution cooperative power flow calculation and has better robustness.
TABLE 1
Claims (1)
1. A transmission and distribution cooperative rapid robust power flow calculation method is characterized by comprising the following steps:
1) initialization: setting initial values of boundary node voltage amplitudes
Initial value of voltage phase angle
Setting convergence accuracy eta; setting a maximum iteration number K; setting the current iteration number k to be 0;
2) and (3) power distribution network load flow calculation: according to the current boundary node voltage amplitude
Phase angle of voltage
Calculating the power flow of the power distribution network to obtain the node voltage amplitude of each node of the power distribution network
Phase angle of voltage
Computing injected power of boundary nodes
And will be
Transmitting to the transmission grid;
3) and (3) power transmission network load flow calculation: injecting power according to the current boundary node
Calculating the power flow of the power transmission network to obtain the node voltage amplitude of each node of the power transmission network
Phase angle of voltage
And new boundary node voltage amplitude
Phase angle of voltage
4) And (3) testing the convergence: if it is not
The algorithm is converged, and the algorithm is terminated, otherwise, the step 5) is carried out;
5) iterative variable correction: if the current iteration times k are odd, directly jumping to the step 6); if the current iteration number k is an even number, updating the value of the iteration variable according to the following formula:
6) updating the current iteration times: let k be k + 1;
7) and (5) checking the maximum iteration times: if K is equal to K, the algorithm reaches the maximum iteration times without convergence, and the algorithm is terminated, otherwise, the current boundary node voltage amplitude value is used
Phase angle of voltage
Transmitting to the distribution network and turning to step 2).
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202011215316.4A CN112398134A (en) | 2020-11-04 | 2020-11-04 | Transmission and distribution cooperative rapid robust load flow calculation method |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202011215316.4A CN112398134A (en) | 2020-11-04 | 2020-11-04 | Transmission and distribution cooperative rapid robust load flow calculation method |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN112398134A true CN112398134A (en) | 2021-02-23 |
Family
ID=74598726
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202011215316.4A Pending CN112398134A (en) | 2020-11-04 | 2020-11-04 | Transmission and distribution cooperative rapid robust load flow calculation method |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN112398134A (en) |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN105470955A (en) * | 2015-12-25 | 2016-04-06 | 国家电网公司 | Newton-Raphson and forward pushing coordinated flow algorithm for transmission and distribution network comprising distributed power supply |
| CN110148946A (en) * | 2019-05-24 | 2019-08-20 | 浙江大学 | A kind of isolated island micro-capacitance sensor tidal current computing method solved based on two step of confactor |
| CN111181164A (en) * | 2020-02-11 | 2020-05-19 | 山东大学 | Improved master-slave split transmission and distribution cooperative power flow calculation method and system |
| CN111463794A (en) * | 2020-04-29 | 2020-07-28 | 山东大学 | Transmission and distribution network cooperative power flow control method and system with additional earth impedance |
-
2020
- 2020-11-04 CN CN202011215316.4A patent/CN112398134A/en active Pending
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN105470955A (en) * | 2015-12-25 | 2016-04-06 | 国家电网公司 | Newton-Raphson and forward pushing coordinated flow algorithm for transmission and distribution network comprising distributed power supply |
| CN110148946A (en) * | 2019-05-24 | 2019-08-20 | 浙江大学 | A kind of isolated island micro-capacitance sensor tidal current computing method solved based on two step of confactor |
| CN111181164A (en) * | 2020-02-11 | 2020-05-19 | 山东大学 | Improved master-slave split transmission and distribution cooperative power flow calculation method and system |
| CN111463794A (en) * | 2020-04-29 | 2020-07-28 | 山东大学 | Transmission and distribution network cooperative power flow control method and system with additional earth impedance |
Non-Patent Citations (1)
| Title |
|---|
| KUNJIE TANG: "Successive-Intersection-Approximation-Based Power Flow Method for Integrated Transmission and Distribution Networks", 《IEEE》 * |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| US10222404B2 (en) | General load flow calculation method for power systems with unified power flow controller | |
| US11050248B2 (en) | Optimization model for quick track of SVSR boundary of power system | |
| Eltamaly et al. | Load flow analysis by gauss-seidel method; a survey | |
| CN105514971A (en) | Flow calculation method suitable for microgrids in various operation modes | |
| CN109617080B (en) | Rectangular coordinate Newton method load flow calculation method based on improved Jacobian matrix | |
| CN102420427A (en) | Method for calculating voltage stability margin of regional power grid by considering outer network equivalence | |
| CN107529644A (en) | A kind of power system steady state voltage stability domain border linear approximation method | |
| CN103632046A (en) | Power grid load flow calculation method | |
| CN112636358A (en) | Power system load flow calculation method based on multivariable difference-of-quotient method | |
| CN103956735A (en) | Harmonic power flow analysis method of distributed power generation system | |
| CN114204564B (en) | Short-circuit current calculation method and device for power grid containing inversion type new energy | |
| CN117498433A (en) | Transient stability power limit quantization method for hybrid parallel system | |
| CN106532712A (en) | Rectangular coordinate Newton method load flow calculation method for small-impedance-branch-containing power grid based on compensation method | |
| CN109494748A (en) | Newton load flow calculation method based on node type and modified Jacobian matrix | |
| CN111049144B (en) | Active power distribution network load flow calculation PV node processing method based on constant change rate matrix | |
| CN114336635B (en) | Full-pure embedded power flow calculation method and device based on constant term value and priori node | |
| CN112383063B (en) | Transmission and distribution cooperative power flow calculation method based on least square fitting correction | |
| CN106410811B (en) | Iteration small impedance branches endpoint changes the tidal current computing method of Jacobian matrix for the first time | |
| CN112398134A (en) | Transmission and distribution cooperative rapid robust load flow calculation method | |
| CN111835029A (en) | Method and system for determining resonance stability of direct-current transmission sending system | |
| CN112383064B (en) | Transmission and distribution cooperative power flow calculation method based on iteration variable historical value updating | |
| CN114188945B (en) | Method and device for calculating short-circuit current of power distribution network containing photovoltaic power supply | |
| CN111884226A (en) | Power grid probabilistic power flow analysis method based on generalized semi-invariant and maximum entropy method | |
| CN103199524A (en) | Load flow calculation method adapting to multiple distributed power access | |
| CN107465195B (en) | Optimal power flow double-layer iteration method based on micro-grid combined power flow calculation |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| WD01 | Invention patent application deemed withdrawn after publication | ||
| WD01 | Invention patent application deemed withdrawn after publication |
Application publication date: 20210223 |