CN110635485B - Load flow model calculation method for multi-terminal direct current system self-adaptive droop control - Google Patents
Load flow model calculation method for multi-terminal direct current system self-adaptive droop control Download PDFInfo
- Publication number
- CN110635485B CN110635485B CN201910956096.1A CN201910956096A CN110635485B CN 110635485 B CN110635485 B CN 110635485B CN 201910956096 A CN201910956096 A CN 201910956096A CN 110635485 B CN110635485 B CN 110635485B
- Authority
- CN
- China
- Prior art keywords
- droop control
- adaptive droop
- direct current
- adaptive
- current
- 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
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
Landscapes
- Engineering & Computer Science (AREA)
- Power Engineering (AREA)
- Supply And Distribution Of Alternating Current (AREA)
Abstract
The method for calculating the power flow of the multi-terminal direct current system with the self-adaptive droop control comprises the following steps of: (1) firstly, assuming an initial value of the node voltage as a voltage reference value, and calculating to obtain a node current through a load flow equation; (2) calculating a self-adaptive droop control parameter according to the node current; (3) performing first-order Taylor expansion on the adaptive droop control equation according to the magnitude relation between the node current and the reference current to obtain a corresponding adaptive droop control linear equation; (4) calculating to obtain updated node current by combining a multi-terminal direct current system power flow formula; (5) judging whether the updated node current and the node current before updating have the same size relationship with the reference current, and if the same size relationship does not exist, repeating the steps (2) - (5) according to the modified node current; if the same magnitude relationship exists, the loop is ended.
Description
Technical Field
The invention relates to a power system load flow calculation method, in particular to a multi-terminal direct current system load flow calculation method with self-adaptive droop control.
Background
So far, many scholars consider that basic control strategies of a multi-terminal direct-current system can be divided into three types, namely master-slave control, voltage margin control, direct-current voltage droop control and the like, wherein the direct-current voltage droop control can reasonably set the proportion of droop control coefficients of each converter station in the multi-terminal direct-current system, so that each converter station simultaneously participates in power regulation to improve the active power distribution of the system.
In recent years, scholars of Liuying culture and the like of the university of North China electric power propose a self-adaptive droop control method by improving V-I droop control.
However, the adaptive droop control equation has strong nonlinear characteristics, and brings great obstruction to load flow calculation and optimization planning of a multi-terminal direct-current power system.
Disclosure of Invention
The technical problem to be solved by the invention is to overcome the defects in the prior art and provide a power flow calculation method for a multi-terminal direct current system with adaptive droop control.
The technical scheme adopted by the invention for solving the technical problem is that the load flow calculation method for the multi-terminal direct current system with the self-adaptive droop control comprises the following steps:
step 1: initial bus voltage of DC bus assuming connection of VSC converter stationThe injection power of the initial node on the DC bus can be calculated through a first direct current flow formula, and then the injection current of the initial node on the DC bus is obtained
Wherein, the reference source of the formula I is Roger Wiget,Andersson,“DC Optimal Power Flow Including HVDC Grids”,2013 IEEE Electrical Power&Energy Conference.
step 2: determining alpha according to a formula II, wherein alpha is a parameter in the adaptive droop control equation
And step 3: according toThe size of the droop control system is represented by a self-adaptive droop control formulaA first order taylor expansion is performed.
Wherein, the reference sources of the formula two and the formula three are as follows: yingpei Liu, La Zhang and Haiping Liung, "DC Voltage Adaptive Droop Control Stratagy for a Hybrid Multi-Terminal HVDC System", Energies, 2019;
when in useAnd when the adaptive droop control linear equation is obtained by performing first-order Taylor expansion on the formula III, the adaptive droop control linear equation is as follows:
when in useAnd when the adaptive droop control linear equation is obtained by performing first-order Taylor expansion on the formula III, the adaptive droop control linear equation is as follows:
and 4, step 4: combining the adaptive droop control linear equation obtained in the step 3 with a direct current system power flow formula I, and calculating to obtain idc;
And 5: judgment of idcAndwhether simultaneously greater or less thanIf yes, then idcFor the final solution, if not, letAnd returning to the step 2.
The method can effectively process the nonlinear characteristics of the power flow model of the multi-terminal direct current system with the self-adaptive droop control.
Drawings
FIG. 1 is a block flow diagram of the method of the present invention;
FIG. 2 is a topology diagram of a multi-terminal DC system used in an embodiment;
FIG. 3 is a comparison graph of DC1 point current values calculated by the method of the present invention versus a method that does not employ droop control equation approximation;
FIG. 4 is a comparison graph of DC2 point current values calculated by the method of the present invention versus a method that does not employ droop control equation approximation;
FIG. 5 is a graph comparing the voltage value at DC1 calculated by the method of the present invention with a method that does not use the droop control equation approximation;
FIG. 6 is a comparison graph of voltage values at DC2 calculated by the method of the present invention compared to a method that does not employ the droop control equation approximation.
Detailed Description
The present invention will be described in further detail with reference to the accompanying drawings and examples.
The specific flow of this embodiment is shown in fig. 1, and the adopted topology structure of the multi-terminal dc system is shown in fig. 2. Ten scenarios were randomly generated from the normal distribution (the loads on the DC3 and DC4 nodes obeyed a normal distribution with expected values of 0.3p.u. and 0.6p.u., respectively, and a standard deviation of 10% of the expected values).
And carrying out load flow calculation by using the calculation method of the invention aiming at the ten scenes.
Bus voltage reference for DC1 and DC2Are all set to 1.0p.u., current reference valueSet to 0.26p.u. and 0.7p.u., respectively, converter station capacity PmaxSet to 0.5p.u. and 1.0p.u., respectively, and the droop control coefficient K is set to 2.
The calculation method comprises the following steps:
step 1: assuming that there is a bus voltage on the DC bus of the VSC converter station connected theretoCalculating to obtain the injection power of the upper node of the DC bus by combining with the first direct current flow formula, and further obtaining the injection current of the upper node of the DC bus
Wherein, the reference source of the formula I is Roger Wiget,Andersson,“DC Optimal Power Flow Including HVDC Grids”,2013 IEEE Electrical Power&Energy Conference.
|
|
|
|
|
|
|
|
|
|
|
DC1 | 0.2281 | 0.2995 | 0.2247 | 0.2484 | 0.2527 | 0.2229 | 0.2335 | 0.2631 | 0.2973 | 0.2886 |
DC2 | 0.6071 | 0.8376 | 0.6511 | 0.6737 | 0.6256 | 0.6256 | 0.6461 | 0.7366 | 0.7946 | 0.7795 |
Step 2: judgment ofAnd the reference value of the droop control device, and accordingly determining alpha, which is a parameter in the adaptive droop control equation.
And step 3: according toThe size of the droop control system is represented by a self-adaptive droop control formulaA first order taylor expansion is performed.
Wherein, the reference sources of the formula two and the formula three are as follows: yingpei Liu, La Zhang and Haiping Liung, "DC Voltage Adaptive Droop Control Stratagy for a Hybrid Multi-Terminal HVDC System", Energies, 2019;
when in useAnd when the adaptive droop control linear equation is obtained by performing first-order Taylor expansion on the formula III, the adaptive droop control linear equation is as follows:
and 4, step 4: combining the adaptive droop control linear equation obtained in the step 3 with a direct current system power flow formula I, and calculating to obtain idcThe results are shown in the following table.
And 5: judgment of idcAndwhether simultaneously greater or less thanWith scenario 5 i on DC1 busdcAnddoes not satisfy the condition, orderAnd returning to the step 2 and entering a loop.
And (3) circulation 1: to obtain idcThe results are shown in the following table.
|
|
|
|
|
|
|
|
|
|
|
DC1 | 0.2237 | 0.3591 | 0.2456 | 0.2523 | 0.2581 | 0.2332 | 0.2422 | 0.2845 | 0.3345 | 0.3220 |
DC2 | 0.6114 | 0.7780 | 0.6301 | 0.6698 | 0.6944 | 0.6153 | 0.6373 | 0.7152 | 0.7574 | 0.7461 |
With scenario 5 i on DC1 busdcAnddoes not satisfy the condition, orderAnd returning to the step 2 and entering a loop.
And (3) circulation 2: to obtain idcThe results are shown in the following table.
|
|
|
|
|
|
|
|
|
|
|
DC1 | 0.2237 | 0.3591 | 0.2456 | 0.2523 | 0.2591 | 0.2332 | 0.2422 | 0.2845 | 0.3345 | 0.3220 |
DC2 | 0.6114 | 0.7780 | 0.6301 | 0.6698 | 0.6933 | 0.6153 | 0.6373 | 0.7152 | 0.7574 | 0.7461 |
All buses under all scenesAll satisfy the condition idcIs the final solution.
Fig. 3 to fig. 6 are schematic diagrams comparing the calculation result of the method with the calculation result of the load flow algorithm which does not adopt the linear approximation of the droop control equation, and the results show that the errors of the two calculation results are very close, which also shows that the calculation method of the invention is feasible and effective.
Various modifications and variations of the embodiments of the present invention may be made by those skilled in the art, and they are still within the scope of the present patent invention, provided they are within the scope of the claims and their equivalents.
What is not described in detail in the specification is prior art that is well known to those skilled in the art.
Claims (1)
1. The method for calculating the power flow of the multi-terminal direct current system with the self-adaptive droop control is characterized by comprising the following steps of:
step 1: DC bus voltage assuming connection of VSC converter stationThe injection power of the upper node of the DC bus can be calculated through a first direct current flow formula, and then the injection current of the upper node of the DC bus is obtained
Step 2: determining alpha according to the formula II, wherein the alpha is a parameter in the adaptive droop control equation
And step 3: according toThe size of the droop control system is represented by a self-adaptive droop control formulaPerforming first-order Taylor expansion;
when in useAnd when the adaptive droop control linear equation is obtained by performing first-order Taylor expansion on the formula III, the adaptive droop control linear equation is as follows:
when in useAnd when the adaptive droop control linear equation is obtained by performing first-order Taylor expansion on the formula III, the adaptive droop control linear equation is as follows:
and 4, step 4: combining the adaptive droop control linear equation obtained in the step 3 with a direct current system power flow formula I, and calculating to obtain idc;
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910956096.1A CN110635485B (en) | 2019-10-09 | 2019-10-09 | Load flow model calculation method for multi-terminal direct current system self-adaptive droop control |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910956096.1A CN110635485B (en) | 2019-10-09 | 2019-10-09 | Load flow model calculation method for multi-terminal direct current system self-adaptive droop control |
Publications (2)
Publication Number | Publication Date |
---|---|
CN110635485A CN110635485A (en) | 2019-12-31 |
CN110635485B true CN110635485B (en) | 2021-06-15 |
Family
ID=68975955
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910956096.1A Active CN110635485B (en) | 2019-10-09 | 2019-10-09 | Load flow model calculation method for multi-terminal direct current system self-adaptive droop control |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN110635485B (en) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111181189A (en) * | 2020-01-07 | 2020-05-19 | 中国电力科学研究院有限公司 | Nonlinear droop control method and system for multi-terminal direct current transmission system |
Family Cites Families (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US7940105B2 (en) * | 2008-08-08 | 2011-05-10 | Beckman Coulter, Inc. | High-resolution parametric signal restoration |
EP2688191A1 (en) * | 2012-07-17 | 2014-01-22 | ABB Research Ltd. | Multi terminal HVDC control |
CN104167756B (en) * | 2014-08-18 | 2017-02-15 | 国家电网公司 | Power flow determination method of alternating current and direct current system containing multi-terminal high voltage direct current transmission |
CN107171354B (en) * | 2017-06-19 | 2020-01-14 | 天津大学 | Method for calculating droop coefficient of converter station of flexible direct current power grid |
CN109494746B (en) * | 2018-11-08 | 2021-11-02 | 国网甘肃省电力公司电力科学研究院 | Island alternating current-direct current series-parallel micro-grid load flow calculation method based on improved adaptive droop control |
CN110148946B (en) * | 2019-05-24 | 2022-03-25 | 浙江大学 | Island microgrid load flow calculation method based on two-step solution of auxiliary factors |
-
2019
- 2019-10-09 CN CN201910956096.1A patent/CN110635485B/en active Active
Also Published As
Publication number | Publication date |
---|---|
CN110635485A (en) | 2019-12-31 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Peyghami et al. | On secondary control approaches for voltage regulation in DC microgrids | |
Tabari et al. | A mathematical model for stability analysis of a DC distribution system for power system integration of plug-in electric vehicles | |
Song et al. | Cost-based adaptive droop control strategy for VSC-MTDC system | |
CN110350597B (en) | Distribution network voltage control method based on distributed photovoltaic active-reactive optimization | |
CN103746583A (en) | MMC (Modular Multilevel Converter) bridge arm equivalent simulation method | |
CN109301870B (en) | Capacity optimization method for power electronic multi-feed-in power system | |
CN109617112B (en) | Improved direct-current voltage control strategy applicable to multi-terminal flexible direct-current system | |
CN110635485B (en) | Load flow model calculation method for multi-terminal direct current system self-adaptive droop control | |
KR20150045223A (en) | Method and Apparatus for Controlling Doubly-fed Induction Generator using Adaptive Backstepping Control Scheme | |
CN109494787B (en) | Grid-connected capacity limit calculation method for photovoltaic multi-feed-in power system | |
Zaery et al. | Distributed economic dispatch for islanded DC microgrids based on finite-time consensus protocol | |
CN109687499B (en) | Voltage balance control method for series connection converter valve of flexible direct current converter station | |
CN112531715B (en) | Droop control multi-terminal direct-current micro-grid power flow calculation method based on virtual resistor | |
CN115579967B (en) | Power adjustment method and system for large power grid AC/DC multi-section | |
CN111697607A (en) | Multi-terminal flexible direct-current transmission receiving-end power grid access method and system | |
Meng et al. | Adaptive control of energy storage systems for power smoothing applications | |
CN107612290B (en) | Optimization control method and system for capacitor voltage sequencing frequency of current converter | |
Dehaghania et al. | A Distributed Cooperative Secondary Control Scheme for Obtaining Power and Voltage References of Distributed Generations in Islanded DC Microgrids | |
CN109560561B (en) | Dynamic simulation method, system and storage medium for three-phase asymmetric operation of active power distribution network | |
CN109449968B (en) | Power electronic transformer and AC/DC source network load multi-current equipment integration method | |
CN111864782A (en) | PI parameter optimization method and system in MMC control system of flexible direct current power distribution network | |
CN108933442B (en) | Uncertain distribution network voltage optimization control method | |
Pathak et al. | Optimal Power Transmission of Offshore Wind Farms with Droop Based Control Strategies of Multi-terminal HVDC Systems | |
Thau et al. | EMT implementation and validation of MPC for VSC-HVDC embedded in AC meshed grid | |
Yi et al. | Automatic Parameter Optimization Based on High-Efficiency Simulation Model of Regional Energy Router |
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 |