CN112383063B - Transmission and distribution cooperative power flow calculation method based on least square fitting correction - Google Patents

Transmission and distribution cooperative power flow calculation method based on least square fitting correction Download PDF

Info

Publication number
CN112383063B
CN112383063B CN202011214553.9A CN202011214553A CN112383063B CN 112383063 B CN112383063 B CN 112383063B CN 202011214553 A CN202011214553 A CN 202011214553A CN 112383063 B CN112383063 B CN 112383063B
Authority
CN
China
Prior art keywords
power
distribution network
flow calculation
voltage
node
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.)
Active
Application number
CN202011214553.9A
Other languages
Chinese (zh)
Other versions
CN112383063A (en
Inventor
唐坤杰
董树锋
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Zhejiang University ZJU
Original Assignee
Zhejiang University ZJU
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Zhejiang University ZJU filed Critical Zhejiang University ZJU
Priority to CN202011214553.9A priority Critical patent/CN112383063B/en
Publication of CN112383063A publication Critical patent/CN112383063A/en
Application granted granted Critical
Publication of CN112383063B publication Critical patent/CN112383063B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

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 power flow calculation method based on least square fitting correction, which comprises the following seven main steps: initialization, power distribution network load flow calculation, power transmission network load flow calculation, convergence inspection, fitting correction based on a least square method, current iteration number updating and maximum iteration number inspection. According to the method, the least square method is used for fitting and correcting the iteration variables, the efficiency and the robustness of the overall algorithm of the power transmission and distribution cooperative power flow calculation can be remarkably improved, and particularly, when the power distribution network has large disturbance on the power transmission network, namely the sensitivity of the boundary node injection power on the boundary node voltage amplitude is high, the problem of divergence of the traditional master-slave splitting method can be solved.

Description

Transmission and distribution cooperative power flow calculation method based on least square fitting correction
Technical Field
The invention relates to a power grid power flow analysis method, relates to a power transmission network and power distribution network collaborative global power flow analysis method, and particularly relates to a transmission and distribution collaborative power flow calculation method based on least square fitting correction.
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 power flow calculation method based on least square fitting correction so as to improve the overall calculation efficiency and convergence.
Therefore, the invention adopts the following technical scheme:
a transmission and distribution cooperative power flow calculation method based on least square fitting correction comprises the following steps:
1) initialization: setting initial value of voltage amplitude of boundary node of power distribution network
Figure BDA0002759913380000021
Initial value of voltage phase angle
Figure BDA0002759913380000022
Setting convergence accuracy eta; setting a maximum iteration number K; setting a memory capacity m; and 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
Figure BDA0002759913380000023
Phase angle of voltage
Figure BDA0002759913380000024
Calculating the power flow of the power distribution network to obtain the node voltage amplitude of each node of the power distribution network
Figure BDA0002759913380000025
Phase angle of voltage
Figure BDA0002759913380000026
Computing power injection for border nodes
Figure BDA0002759913380000027
And will be
Figure BDA0002759913380000028
And transmitted to the transmission grid.
3) And (3) power transmission network load flow calculation: injecting power according to the current boundary node of the power distribution network
Figure BDA0002759913380000029
Calculating the power flow of the power transmission network to obtain the node voltage amplitude of each node of the power transmission network
Figure BDA00027599133800000210
Phase angle of voltage
Figure BDA00027599133800000211
And new distribution network boundary node voltage amplitude
Figure BDA00027599133800000212
Phase angle of voltage
Figure BDA00027599133800000213
4) And (3) testing the convergence: if it is not
Figure BDA00027599133800000214
The algorithm converges and terminates, otherwise go to step 5).
5) Fitting correction based on a least square method: if k is<m, turning to the step 6); if k is larger than or equal to m, solving the following optimization problem to obtain alpha0，α1，…，αm：
Figure BDA00027599133800000215
Figure BDA00027599133800000216
The iteration variables are then corrected using:
Figure BDA00027599133800000217
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 number of iterations without convergence, terminates the algorithm,otherwise, the current boundary node voltage amplitude value is obtained
Figure BDA0002759913380000031
Phase angle of voltage
Figure BDA0002759913380000032
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:
according to the method, the least square method is used for fitting and correcting the iteration variables, the efficiency and the robustness of the overall algorithm of the power transmission and distribution cooperative power flow calculation can be remarkably improved, and particularly, when the power distribution network has large disturbance on the power transmission network, namely the sensitivity of the boundary node injection power on the boundary node voltage amplitude is high, the problem of divergence of the traditional master-slave splitting method can be solved. The method has strong operability and is easy to realize.
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, fitting correction based on a least square method, current iteration number updating and maximum iteration number inspection.
Step S1: setting initial values of boundary node voltage amplitudes
Figure BDA0002759913380000033
Initial value of voltage phase angle
Figure BDA0002759913380000034
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
Figure BDA0002759913380000035
Phase angle of voltage
Figure BDA0002759913380000036
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
Figure BDA0002759913380000041
Phase angle of voltage
Figure BDA0002759913380000042
Computing power injection for border nodes
Figure BDA0002759913380000043
And will be
Figure BDA0002759913380000044
And transmitted to the transmission grid.
Step S3: based on current boundary node power injection
Figure BDA0002759913380000045
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
Figure BDA0002759913380000046
Phase angle of voltage
Figure BDA0002759913380000047
And new boundary node voltage amplitude
Figure BDA0002759913380000048
Phase angle of voltage
Figure BDA0002759913380000049
Step S4: if it is not
Figure BDA00027599133800000410
The algorithm ends with convergence, otherwise it proceeds to step S5.
Step S5: if k is<m, go to step S6; if k is larger than or equal to m, solving the following optimization problem by using a convex optimization solver to obtain alpha0，α1，…，αm：
Figure BDA00027599133800000411
Figure BDA00027599133800000412
The iteration variables are then corrected using:
Figure BDA00027599133800000413
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
Figure BDA00027599133800000414
Phase angle of voltage
Figure BDA00027599133800000415
Transmitting to the distribution network and turning to step 2).
The method for calculating power distribution network power flow based on least square fitting correction shown in fig. 1 is a specific embodiment of the present invention, and already embodies the substantial features and the improvements of the present invention, and those skilled in the art can calculate power distribution network power flow by using different calculation methods in step S2, calculate power transmission network power flow by using different calculation methods in step S3, solve the optimization problem by using different calculation methods in step S5, or make equivalent modifications to the calculation method in terms of structure and the like according to the practical use requirements.
In order to further embody the substantive characteristics and remarkable 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. The calculation example A 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, the No. 8, No. 15 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. The memory capacity m is set to 3. Initial value of boundary node voltage amplitude
Figure BDA0002759913380000051
Set to 1, initial value of voltage phase angle
Figure BDA0002759913380000052
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, which takes 78.8ms, while the method of the invention only needs 7 iterations of convergence, which takes 31.5ms, which shows that the method of the invention 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, while the method of the invention needs 9 convergence times and takes 40.8ms, 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
Figure BDA0002759913380000061

Claims (1)

1. A transmission and distribution cooperative power flow calculation method based on least square fitting correction is characterized by comprising the following steps:
1) initialization: setting initial value of voltage amplitude of boundary node of power distribution network
Figure FDA0002759913370000011
Initial value of voltage phase angle
Figure FDA0002759913370000012
Setting convergence accuracy eta; setting a maximum iteration number K; setting a memory capacity m; 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
Figure FDA0002759913370000013
Phase angle of voltage
Figure FDA0002759913370000014
Calculating the power flow of the power distribution network to obtain the node voltage amplitude of each node of the power distribution network
Figure FDA0002759913370000015
Phase angle of voltage
Figure FDA0002759913370000016
Computing injected power of boundary nodes
Figure FDA0002759913370000017
And will be
Figure FDA0002759913370000018
Transmitting to the transmission grid;
3) and (3) power transmission network load flow calculation: injecting power according to the current boundary node
Figure FDA0002759913370000019
Calculating the power flow of the power transmission network to obtain the node voltage amplitude of each node of the power transmission network
Figure FDA00027599133700000110
Phase angle of voltage
Figure FDA00027599133700000111
And new boundary node voltage amplitude
Figure FDA00027599133700000112
Phase angle of voltage
Figure FDA00027599133700000113
4) And (3) testing the convergence: if it is not
Figure FDA00027599133700000114
The algorithm is converged, and the algorithm is terminated, otherwise, the step 5) is carried out;
5) fitting correction based on a least square method: if k is<m, turning to the step 6); if k is larger than or equal to m, solving the following optimization problem to obtain alpha0,α1,…,αm：
Figure FDA00027599133700000115
Figure FDA00027599133700000116
The iteration variables are then corrected using:
Figure FDA00027599133700000117
6) updating the current iteration times: let k be k + 1;
7) reach maximum iterationAnd (4) checking the 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
Figure FDA0002759913370000021
Phase angle of voltage
Figure FDA0002759913370000022
Transmitting to the distribution network and turning to step 2).
CN202011214553.9A 2020-11-04 2020-11-04 Transmission and distribution cooperative power flow calculation method based on least square fitting correction Active CN112383063B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202011214553.9A CN112383063B (en) 2020-11-04 2020-11-04 Transmission and distribution cooperative power flow calculation method based on least square fitting correction

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202011214553.9A CN112383063B (en) 2020-11-04 2020-11-04 Transmission and distribution cooperative power flow calculation method based on least square fitting correction

Publications (2)

Publication Number Publication Date
CN112383063A CN112383063A (en) 2021-02-19
CN112383063B true CN112383063B (en) 2022-04-08

Family

ID=74578292

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202011214553.9A Active CN112383063B (en) 2020-11-04 2020-11-04 Transmission and distribution cooperative power flow calculation method based on least square fitting correction

Country Status (1)

Country Link
CN (1) CN112383063B (en)

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105356456A (en) * 2015-11-18 2016-02-24 中国电力科学研究院 Boundary-node-equivalence-based combined power flow calculation method of power transmission network and distribution network
CN109274091A (en) * 2018-10-15 2019-01-25 同济大学 A kind of transmission & distribution integration parallel state estimation method
CN110687352A (en) * 2019-10-21 2020-01-14 国网上海市电力公司 PMU-based power distribution network T-type line parameter identification method and system
CN110829411A (en) * 2019-10-25 2020-02-21 国网天津市电力公司电力科学研究院 Flexible power distribution network operation domain observation method considering alternating current power flow
CN111463794A (en) * 2020-04-29 2020-07-28 山东大学 Transmission and distribution network cooperative power flow control method and system with additional earth impedance

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US9619848B2 (en) * 2013-05-06 2017-04-11 Sas Institute Inc. Techniques to determine settings for an electrical distribution network

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105356456A (en) * 2015-11-18 2016-02-24 中国电力科学研究院 Boundary-node-equivalence-based combined power flow calculation method of power transmission network and distribution network
CN109274091A (en) * 2018-10-15 2019-01-25 同济大学 A kind of transmission & distribution integration parallel state estimation method
CN110687352A (en) * 2019-10-21 2020-01-14 国网上海市电力公司 PMU-based power distribution network T-type line parameter identification method and system
CN110829411A (en) * 2019-10-25 2020-02-21 国网天津市电力公司电力科学研究院 Flexible power distribution network operation domain observation method considering alternating current power flow
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 (2)

* Cited by examiner, † Cited by third party
Title
Successive-Intersection-Approximation-Based Power Flow Method for Integrated Transmission and Distribution Networks;Kunjie Tang et al.;《IEEE》;20200514;第35卷(第6期);第4836-4846页 *
计及输入变量强相关性的概率潮流计算模型;别佩 等;《湖北工业大学学报》;20140228;第29卷(第1期);第49-52页 *

Also Published As

Publication number Publication date
CN112383063A (en) 2021-02-19

Similar Documents

Publication Publication Date Title
WO2018103317A1 (en) Universal power flow calculation method for power system comprising upfc
Eltamaly et al. Load flow analysis by gauss-seidel method; a survey
CN107332240A (en) The method of power system steady state voltage stability domain boundary search based on Optimized model
CN109167366B (en) Method for calculating static voltage stability critical point of power system
CN109617080B (en) Rectangular coordinate Newton method load flow calculation method based on improved Jacobian matrix
CN114336635B (en) Full-pure embedded power flow calculation method and device based on constant term value and priori node
CN107529644A (en) A kind of power system steady state voltage stability domain border linear approximation method
CN103632046A (en) Power grid load flow calculation method
CN109494748B (en) Newton method load flow calculation method based on node type and modified Jacobian matrix
CN117498433A (en) Transient stability power limit quantization method for hybrid parallel system
CN106856327B (en) A kind of compensation of line series containing small impedance branches algorithm quicksort tidal current computing method
CN108808681A (en) Grid-connected tidal current computing method based on mixed injection model
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
CN111049144B (en) Active power distribution network load flow calculation PV node processing method based on constant change rate matrix
Sun et al. Distributed power flow calculation for whole networks including transmission and distribution
CN111884226A (en) Power grid probabilistic power flow analysis method based on generalized semi-invariant and maximum entropy method
CN111835029A (en) Method and system for determining resonance stability of direct-current transmission sending system
CN107465195B (en) Optimal power flow double-layer iteration method based on micro-grid combined power flow calculation
CN112383064B (en) Transmission and distribution cooperative power flow calculation method based on iteration variable historical value updating
CN106712029B (en) The Newton load flow calculation method of small impedance branches PQ endpoint change Jacobian matrix
CN112398134A (en) Transmission and distribution cooperative rapid robust load flow calculation method
Wang et al. Load flow analysis
CN114188945B (en) Method and device for calculating short-circuit current of power distribution network containing photovoltaic power supply
CN110046450B (en) Initial value selection method suitable for Newton method load flow calculation of superconducting cable-containing power grid

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
GR01 Patent grant
GR01 Patent grant