CN114024338B - Large-scale wind power collection grid split-phase power flow optimization method and system - Google Patents

Large-scale wind power collection grid split-phase power flow optimization method and system Download PDF

Info

Publication number
CN114024338B
CN114024338B CN202111414460.5A CN202111414460A CN114024338B CN 114024338 B CN114024338 B CN 114024338B CN 202111414460 A CN202111414460 A CN 202111414460A CN 114024338 B CN114024338 B CN 114024338B
Authority
CN
China
Prior art keywords
phase
node
phi
voltage
power
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
CN202111414460.5A
Other languages
Chinese (zh)
Other versions
CN114024338A (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.)
State Grid Corp of China SGCC
North China Electric Power University
Electric Power Research Institute of State Grid Jibei Electric Power Co Ltd
Original Assignee
State Grid Corp of China SGCC
North China Electric Power University
Electric Power Research Institute of State Grid Jibei Electric Power Co Ltd
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 State Grid Corp of China SGCC, North China Electric Power University, Electric Power Research Institute of State Grid Jibei Electric Power Co Ltd filed Critical State Grid Corp of China SGCC
Priority to CN202111414460.5A priority Critical patent/CN114024338B/en
Publication of CN114024338A publication Critical patent/CN114024338A/en
Application granted granted Critical
Publication of CN114024338B publication Critical patent/CN114024338B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • HELECTRICITY
    • H02GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02JCIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
    • H02J3/00Circuit arrangements for ac mains or ac distribution networks
    • H02J3/38Arrangements for parallely feeding a single network by two or more generators, converters or transformers
    • H02J3/381Dispersed generators
    • HELECTRICITY
    • H02GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02JCIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
    • H02J3/00Circuit arrangements for ac mains or ac distribution networks
    • H02J3/04Circuit arrangements for ac mains or ac distribution networks for connecting networks of the same frequency but supplied from different sources
    • H02J3/06Controlling transfer of power between connected networks; Controlling sharing of load between connected networks
    • HELECTRICITY
    • H02GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02JCIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
    • H02J3/00Circuit arrangements for ac mains or ac distribution networks
    • H02J3/12Circuit arrangements for ac mains or ac distribution networks for adjusting voltage in ac networks by changing a characteristic of the network load
    • H02J3/14Circuit arrangements for ac mains or ac distribution networks for adjusting voltage in ac networks by changing a characteristic of the network load by switching loads on to, or off from, network, e.g. progressively balanced loading
    • HELECTRICITY
    • H02GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02JCIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
    • H02J2300/00Systems for supplying or distributing electric power characterised by decentralized, dispersed, or local generation
    • H02J2300/20The dispersed energy generation being of renewable origin
    • H02J2300/28The renewable source being wind energy
    • YGENERAL 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
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02BCLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO BUILDINGS, e.g. HOUSING, HOUSE APPLIANCES OR RELATED END-USER APPLICATIONS
    • Y02B70/00Technologies for an efficient end-user side electric power management and consumption
    • Y02B70/30Systems integrating technologies related to power network operation and communication or information technologies for improving the carbon footprint of the management of residential or tertiary loads, i.e. smart grids as climate change mitigation technology in the buildings sector, including also the last stages of power distribution and the control, monitoring or operating management systems at local level
    • Y02B70/3225Demand response systems, e.g. load shedding, peak shaving
    • YGENERAL 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
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02EREDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
    • Y02E10/00Energy generation through renewable energy sources
    • Y02E10/70Wind energy
    • Y02E10/76Power conversion electric or electronic aspects
    • YGENERAL 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
    • Y04INFORMATION OR COMMUNICATION TECHNOLOGIES HAVING AN IMPACT ON OTHER TECHNOLOGY AREAS
    • Y04SSYSTEMS INTEGRATING TECHNOLOGIES RELATED TO POWER NETWORK OPERATION, COMMUNICATION OR INFORMATION TECHNOLOGIES FOR IMPROVING THE ELECTRICAL POWER GENERATION, TRANSMISSION, DISTRIBUTION, MANAGEMENT OR USAGE, i.e. SMART GRIDS
    • Y04S20/00Management or operation of end-user stationary applications or the last stages of power distribution; Controlling, monitoring or operating thereof
    • Y04S20/20End-user application control systems
    • Y04S20/222Demand response systems, e.g. load shedding, peak shaving

Landscapes

  • Engineering & Computer Science (AREA)
  • Power Engineering (AREA)
  • Supply And Distribution Of Alternating Current (AREA)

Abstract

The invention relates to a split-phase power flow optimization method and a split-phase power flow optimization system for a large-scale wind power collection grid, which belong to the technical field of wind power generation, wherein a split-phase power flow optimization model for optimizing control from a system level is constructed, a Pareto archive multi-target particle swarm optimization algorithm is utilized for optimizing and solving the split-phase power flow optimization model, and the compensation capacity of each phase of a collection station in the large-scale wind power collection grid and reactive power compensation devices of all wind power plants at each wind speed is obtained, so that the voltage unbalance degree of all optimized nodes in the large-scale wind power collection grid can meet the requirement, and the three-phase voltage unbalance phenomenon of a wind power collection area is effectively inhibited at the system level from the global angle.

Description

Large-scale wind power collection grid split-phase power flow optimization method and system
Technical Field
The invention relates to the technical field of wind power generation, in particular to a split-phase tide optimization method and system for a large-scale wind power collection grid.
Background
Three-phase voltage unbalance phenomenon often occurs in large-scale wind power centralized grid-connected areas, and even when the phenomenon is serious, a fan is disconnected, so that the safe and stable operation of a power system is greatly adversely affected. At present, the research of current optimization of a wind power collection area is concentrated on economic problems such as minimum network loss, minimum total capacity of a compensation device, and system safety problems such as minimum node voltage deviation, but the research of improving the three-phase voltage unbalance problem of the wind power collection area from the system level through current optimization is lacking, so that a split-phase current optimization model for inhibiting the voltage unbalance of a wind power plant is established, and an algorithm capable of accurately solving the optimization model to obtain the compensation capacity of each phase of each wind power plant and the reactive power compensation device of a collection station is introduced, so that the method has important significance and engineering value.
Disclosure of Invention
The invention aims to provide a split-phase tide optimization method and system for a large-scale wind power collection grid, which can effectively inhibit the phenomenon of three-phase voltage unbalance of a wind power collection area from a global angle in a system level.
In order to achieve the above object, the present invention provides the following solutions:
a split-phase power flow optimization method for a large-scale wind power collection grid, the method comprising:
constructing a split-phase tide optimization model of a large-scale wind power collection grid; the split-phase power flow optimization model comprises an objective function and constraint conditions; the objective function comprises a minimum negative sequence voltage unbalance mean objective function and a minimum voltage deviation objective function;
acquiring line parameters, load data and output power data of a wind power plant of a large-scale wind power collection grid at different wind speeds;
and solving the split-phase power flow optimization model by adopting a Pareto archive multi-objective particle swarm optimization algorithm according to line parameters, load data and wind farm output power data of the large-scale wind power collection grid at each wind speed to obtain the compensation capacity of each phase of the reactive power compensation device of the collection station and each wind farm in the large-scale wind power collection grid at each wind speed.
Optionally, the minimum negative sequence voltage unbalance mean value objective function isWherein F is 1 Is the mean value of the negative sequence voltage unbalance degree, U i1 Is the positive sequence voltage effective value of node i, < ->U i2 Is the negative sequence voltage effective value of node i, < -> And->The voltage effective values of a phase, b phase and c phase of the node i are respectively represented, alpha is a twiddle factor, and n is the total number of nodes;
the minimum voltage deviation objective function isWherein F is 2 For voltage deviation>Actual value of node voltage for the phi-th phase of node i, for example>The node voltage reference value for the phi-th phase of node i.
Optionally, the constraint condition includes: equality constraints and inequality constraints.
Optionally, the equality constraint includes: node current equation constraint, branch voltage equation constraint and three-phase power flow equation constraint;
the node current equation is constrained to beWherein (1)>For phi phase current flowing on the branch connected between node i and node j, a, b and c are respectively a phase, b phase and c phase in three phases, < >>Represents the phi-phase active power of node j, < +.>Phi-phase reactive power, which represents node j, +.>Reactive power representing phase compensation of the reactive compensation means phi of node j +.>Phi-phase voltage, which represents node j, +.>The phi phase current flowing on the first branch connected with the node j is represented, and p is the number of branches connected between the node i and the node j;
the branch voltage equation is constrained to beWherein (1)>Phi-phase voltage, which represents node i, +.>Branch resistance representing the branch connected between node i and node jAn anti-matrix;
the three-phase tide equation is constrained as
Wherein P is i φ For phi-phase active power of node i, N is the number of nodes of the whole collecting power grid, ++>For m-phase voltage of node k, G ik,φm For the conductance between the phi and m phases of nodes i, k, theta ik,φm B is the voltage phase angle difference between phi and m phases of the node i and the node k ik,φm For susceptances between phi and m phases of node i, node k, ++>Phi-phase reactive power, which represents node i, +.>And representing the reactive power compensated by the node i reactive power compensation device phi.
Optionally, the inequality constraint includes: the negative sequence voltage unbalance degree constraint, the reactive compensation device compensates the capacity upper limit constraint and the capacity lower limit constraint and the node voltage safety constraint;
the negative sequence voltage unbalance degree constraint is that
The upper limit and the lower limit of the compensation capacity of the reactive compensation device are constrained as followsWherein (1)> Reactive compensation installed separatelyThe reactive lower limit and the reactive upper limit of the phi phase of the device can be compensated;
the voltage safety constraint of the node is thatWherein (1)>The phi-phase voltages at node i, respectively, allow a lower and an upper limit.
A split-phase power flow optimization system for a large-scale wind power collection grid, the system comprising:
the split-phase tide optimization model construction module is used for constructing a split-phase tide optimization model of the large-scale wind power collection grid; the split-phase power flow optimization model comprises an objective function and constraint conditions; the objective function comprises a minimum negative sequence voltage unbalance mean objective function and a minimum voltage deviation objective function;
the input data acquisition module is used for acquiring line parameters, load data and output power data of the wind power plant of the large-scale wind power collection grid at different wind speeds;
the compensation capacity obtaining module is used for solving the split-phase power flow optimization model by adopting a Pareto archive multi-objective particle swarm optimization algorithm according to line parameters, load data and wind power plant output power data of the large-scale wind power collection grid at each wind speed to obtain the compensation capacity of each phase of the reactive power compensation device of the collection station and each wind power plant in the large-scale wind power collection grid at each wind speed.
Optionally, the minimum negative sequence voltage unbalance mean value objective function isWherein F is 1 Is the mean value of the negative sequence voltage unbalance degree, U i1 Is the positive sequence voltage effective value of node i, < ->U i2 Is the negative sequence voltage effective value of node i, < -> And->The voltage effective values of a phase, b phase and c phase of the node i are respectively represented, alpha is a twiddle factor, and n is the total number of nodes;
the minimum voltage deviation objective function isWherein F is 2 For voltage deviation>Actual value of node voltage for the phi-th phase of node i, for example>The node voltage reference value for the phi-th phase of node i.
Optionally, the constraint condition includes: equality constraints and inequality constraints.
Optionally, the equality constraint includes: node current equation constraint, branch voltage equation constraint and three-phase power flow equation constraint;
the node current equation is constrained to beWherein (1)>For phi phase current flowing on the branch connected between node i and node j, a, b and c are respectively a phase, b phase and c phase in three phases, < >>Represents the phi-phase active power of node j, < +.>Phi-phase reactive power, which represents node j, +.>Reactive power representing phase compensation of the reactive compensation means phi of node j +.>Phi-phase voltage, which represents node j, +.>The phi phase current flowing on the first branch connected with the node j is represented, and p is the number of branches connected between the node i and the node j;
the branch voltage equation is constrained to beWherein (1)>Phi-phase voltage, which represents node i, +.>A branch impedance matrix representing a branch connected between node i and node j;
the three-phase tide equation is constrained as
Wherein P is i φ For phi-phase active power of node i, N is the number of nodes of the whole collecting power grid, ++>For m-phase voltage of node k, G ik,φm For the conductance between the phi and m phases of nodes i, k, theta ik,φm B is the voltage phase angle difference between phi and m phases of the node i and the node k ik,φm For susceptances between phi and m phases of node i, node k, ++>Phi-phase reactive power, which represents node i, +.>And representing the reactive power compensated by the node i reactive power compensation device phi.
Optionally, the inequality constraint includes: the negative sequence voltage unbalance degree constraint, the reactive compensation device compensates the capacity upper limit constraint and the capacity lower limit constraint and the node voltage safety constraint;
the negative sequence voltage unbalance degree constraint is that
The upper limit and the lower limit of the compensation capacity of the reactive compensation device are constrained as followsWherein (1)> The lower limit and the upper limit of reactive power which can be compensated by phi phase of the installed reactive power compensation device are respectively;
the voltage safety constraint of the node is thatWherein (1)>The phi-phase voltages at node i, respectively, allow a lower and an upper limit.
According to the specific embodiment provided by the invention, the invention discloses the following technical effects:
the invention discloses a split-phase power flow optimization method and a split-phase power flow optimization system for a large-scale wind power collection grid, which are characterized in that a split-phase power flow optimization model which is optimally controlled from a system level is constructed, the split-phase power flow optimization model is optimally solved by using a Pareto archive multi-target particle swarm optimization algorithm, and the compensation capacity of each phase of a collection station in the large-scale wind power collection grid and reactive power compensation devices of all wind power fields at each wind speed is obtained, so that the voltage unbalance degree of all optimized nodes in the large-scale wind power collection grid can meet the requirement, and the three-phase voltage unbalance phenomenon of a wind power collection area is effectively inhibited at the system level from the global angle.
Drawings
In order to more clearly illustrate the embodiments of the present invention or the technical solutions of the prior art, the drawings that are needed in the embodiments will be briefly described below, it being obvious that the drawings in the following description are only some embodiments of the present invention, and that other drawings may be obtained according to these drawings without inventive effort for a person skilled in the art.
FIG. 1 is a flow chart of a split-phase power flow optimization method of a large-scale wind power collection grid;
FIG. 2 is a partial equivalent circuit diagram of a large-scale wind power collection grid provided by the invention;
FIG. 3 is a schematic diagram of a Pareto archive multi-target particle swarm algorithm according to the present invention;
fig. 4 is a network structure diagram of a wind power collection area according to an embodiment of the present invention.
Detailed Description
The following description of the embodiments of the present invention will be made clearly and completely with reference to the accompanying drawings, in which it is apparent that the embodiments described are only some embodiments of the present invention, but not all embodiments. All other embodiments, which can be made by those skilled in the art based on the embodiments of the invention without making any inventive effort, are intended to be within the scope of the invention.
The invention aims to provide a split-phase tide optimization method and system for a large-scale wind power collection grid, which can effectively inhibit the phenomenon of three-phase voltage unbalance of a wind power collection area from a global angle in a system level.
In order that the above-recited objects, features and advantages of the present invention will become more readily apparent, a more particular description of the invention will be rendered by reference to the appended drawings and appended detailed description.
The invention provides a split-phase power flow optimization method for a large-scale wind power collection grid, which is shown in figure 1 and comprises the following steps:
step 101, constructing a split-phase tide optimization model of a large-scale wind power collection grid; the split-phase power flow optimization model comprises an objective function and constraint conditions.
Split-phase tide optimizing model
(1) Objective function
The split-phase power flow optimization model of the large-scale wind power collection grid established by the invention mainly comprises 2 optimization targets. Firstly, aiming at the problem that the voltage unbalance of a large-scale wind power collection area and the negative sequence voltage unbalance of a strong wind working condition are overlarge, the voltage unbalance of each bus after three-phase tide optimization needs to be as small as possible. Secondly, in order to ensure the safety and stability of the power system and have good electric energy quality, the deviation between the actual voltage value and the voltage reference value of the collecting bus and each wind power plant access point should be as small as possible.
The objective functions include a minimum negative sequence voltage imbalance mean objective function and a minimum voltage deviation objective function.
1) Minimum negative sequence voltage imbalance mean objective function
According to national standard GB/T15543, the expression of the voltage unbalance degree can be obtained as follows:
then, the node negative sequence voltage balance degree mean expression for optimization is:
the minimum negative sequence voltage unbalance mean value objective function is
Wherein F is 1 Is the mean value of the negative sequence voltage unbalance degree, U i1 Is the positive sequence voltage effective value of the node i, U i2 Is the negative sequence voltage effective value of the node i,and->The voltage effective values of a phase, b phase and c phase of the node i are respectively represented, alpha is a twiddle factor, and n is the total number of nodes.
2) The deviation between the actual value of each phase voltage of the optimization node and the voltage reference value is minimum:
the minimum voltage deviation objective function is
Wherein F is 2 In order to be able to vary the voltage,actual value of node voltage for the phi-th phase of node i, for example>The node voltage reference value for the phi-th phase of node i.
(2) Constraint conditions
The constraint conditions include: equality constraints and inequality constraints.
1) Equation constraint
The equality constraints include: node current equation constraints, branch voltage equation constraints, and three-phase power flow equation constraints.
Assuming that a local equivalent circuit diagram of the large-scale wind power collection grid is shown in fig. 2, according to the local equivalent circuit diagram, the equation constraint conditions to be satisfied by the split-phase power flow optimization model can be listed as follows:
node current equation constraint:
wherein,for phi phase current flowing on the branch connected between node i and node j, a, b and c are respectively a phase, b phase and c phase in three phases, < >>Represents the phi-phase active power of node j, < +.>Phi-phase reactive power, which represents node j, +.>Reactive power representing phase compensation of the reactive compensation means phi of node j +.>Phi-phase voltage, which represents node j, +.>Represents the phi phase current flowing on the first branch connected to node j, p is the number of branches connected between node i and node j, +.>Representing the sum of all phi-phase currents flowing on the branch connected to node j.
Branch voltage equation constraint:
in the method, in the process of the invention,
wherein,phi-phase voltage, which represents node i, +.>A branch impedance matrix representing a branch connected between node i and node j; />For the respective phase self-impedance between nodes i, j, < >> Is the transimpedance between the phases of nodes i, j.
Three-phase tide equation constraint:
wherein P is i φ The phi-phase active power of the node i is N is the number of nodes of the whole collecting power grid,for m-phase voltage of node k, G ik,φm For the conductance between the phi and m phases of nodes i, k, theta ik,φm B is the voltage phase angle difference between phi and m phases of the node i and the node k ik,φm For susceptances between phi and m phases of node i, node k, ++>Phi-phase reactive power, which represents node i, +.>And representing the reactive power compensated by the node i reactive power compensation device phi.
In the view of figure 2,phi phase voltage representing the head node, +.>Represents the apparent power absorbed by the phi-phase of node i from the grid,apparent power absorbed from the grid by phi phase representing node j, < >>Represents the phi phase current flowing on the 1 st branch connected to node j,/>Represents the phi phase current flowing on the mth branch connected to node j,/and>representing the phi-phase voltage of node k.
2) Inequality constraint
Inequality constraints include: the negative sequence voltage unbalance degree constraint, the reactive compensation device compensates the capacity upper limit constraint and the capacity lower limit constraint and the voltage safety constraint of the node.
Negative sequence voltage imbalance constraint:
the national standard GB/T15543 indicates that under the normal operation condition of the power system, the voltage unbalance degree of the public connection point cannot exceed 2 percent and cannot exceed 4 percent in short time, so that the constraint on the negative sequence voltage unbalance degree of the optimized node is introduced:
the reactive compensation device compensates the upper limit constraint and the lower limit constraint of capacity:
wherein,the lower limit and the upper limit of reactive power which can be compensated by phi phase of the installed reactive power compensation device are respectively;
voltage safety constraint of node:
wherein,the phi-phase voltages at node i, respectively, allow a lower and an upper limit.
And 102, acquiring line parameters, load data and output power data of the wind power plant of the large-scale wind power collection grid at different wind speeds.
And 103, solving a split-phase power flow optimization model by adopting a Pareto archive multi-objective particle swarm optimization algorithm according to line parameters, load data and wind farm output power data of the large-scale wind power collection grid at each wind speed, and obtaining the compensation capacity of each phase of a collection station and a reactive power compensation device of each wind farm in the large-scale wind power collection grid at each wind speed.
The split-phase power flow optimization model of the large-scale wind power collection grid has the characteristics of multiple targets, nonlinearity and multiple constraints, and a Pareto archive multi-target particle swarm optimization algorithm is introduced to perform optimization solving aiming at the characteristics of the provided optimization model.
In the algorithm, a solution in the optimization model of the required solution is represented by particles with no volume in a D-dimensional search space, and this solution can be represented by the velocity vector v i =[v i1 ,v i2 ,...,v iD ]And a position vector x i =[x i1 ,x i2 ,...,x iD ]Representing, x i Substituting the objective function F in the optimization model 1 (x),F 2 (x) Then x=x can be found i The function value of each objective function, and all particles constitute a population. Each particle in the population updates its own speed and position through its own and population information, and the update formula is:
where t is the current iteration number and ω is the inertial weight coefficient. The invention adopts linear decreasing inertia weight, namelyWherein omega min 、ω max Respectively, is the maximum value and the minimum value of omega, T is the maximum iteration number, c 1 、c 2 Learning factors of individuals and populations, r 1 And r 2 Is in [0,1 ]]Random number distributed on the upper part, pbest i As the individual optimal position of particle i, gbest is the global optimal position.
Referring to fig. 3, the specific calculation steps are as follows:
step 1: and (3) inputting system parameters: inputting line parameters, load data and wind power plant output power data of a large-scale wind power collection grid;
step 2: algorithm parameter input: setting a population scale N and a learning factor c 1 、c 2 Inertia weight coefficient omega, iteration times T and external file scale N a
Step 3: randomly initializing a population: according to the compensation capacity of the reactive compensation device in the optimization model, randomly initializing the initial speed and the position of an individual in the particle population, substituting the initial speed and the position into an objective function as solutions of the optimization model to calculate F 1 (x i ),F 2 (x i ) In the iteration, a solution that all objective function values are not inferior to other particles and at least one objective function value is superior to other particles is used as a non-inferior solution, and the obtained non-inferior solution is added into an external file for storage;
step 4: determining an initial individual optimum position pbest for each particle i An initial global optimum position gbest;
step 4.1: pbest (p best) i Is selected from the group consisting of: all objective function values of the new solution x of the particles are not inferior to the original pbest i Corresponding objective function values and at least one of the objective function values is better than the pbest i At this time, it is called x-dominant pbest i Then take x as the new pbest i Otherwise, pbest i Remain unchanged; if x is equal to the original pbest i Independent of each other, randomly selecting one as a pbest between the two i
Step 4.2: selection of gbest: adopting a self-adaptive grid algorithm, forming a target space by the change range of the objective function value of each objective function, equally dividing the target space into a plurality of grids, and enabling all non-inferior solutions in the step 3 to have a position corresponding to the target space; defining fitness values for each grid containing non-bad solutions; selecting a grid by adopting a roulette method; one particle is randomly selected from the grid as gbest.
Step 5: updating the position and speed of each particle according to the formula (2-1), and selecting a new pbest of the particle i
Step 6: updating, maintaining and selecting a new gbest mixing process for the external files;
step 6.1: if the number of the non-inferior solutions in the external file is smaller than or equal to a specified value, the newly generated non-inferior solutions x are directly entered into the external file, if x dominates a part of individuals of the external file, x replaces the dominated solution, and at the same time, all particles taking the dominated solution as the gbest select x as the new gbest. If the number of particles in the external file exceeds the file size, selecting a grid with the highest density to randomly select one particle and reject the particle on the basis of obtaining each grid fitness value by using a self-adaptive grid method until the file size is within a required range.
Step 7: judging whether the iteration number reaches the maximum iteration number, if so, terminating the algorithm, and if not, turning to the step 5.
The reactive power optimization concept is applied to improving the electric energy quality, and particularly an optimization model for split-phase power flow optimization is constructed according to the actual problem background of three-phase voltage unbalance of a large-scale wind power collection grid, the optimization model is an optimization model for optimizing control from a system level, and a Pareto archive multi-target particle swarm optimization algorithm is introduced to perform optimization solution. The split-phase power flow optimization method can effectively inhibit the three-phase voltage unbalance phenomenon of the wind power collection area from the global angle at the system level, and is beneficial to the safe and stable operation of the power system under the access of high-proportion new energy.
According to the method for optimizing the split-phase power flow of the large-scale wind power collection grid, which is provided by the invention, the matlab simulation software is utilized to perform split-phase power flow optimization calculation on the matlab simulation software at different wind speeds so as to verify the correctness of the method.
Taking a certain wind power collection area with three-phase voltage imbalance as an example, the network structure of the area is shown in fig. 4.
Assuming that the optimized collection station and the subordinate wind farm are provided with static reactive compensators, the compensation capacity is 30Mvar, and the reactive capacity which can be compensated for each phase is +/-10 Mvar.
When split-phase tide optimization is not performed, the voltage unbalance degree of the collecting station and the subordinate wind power plant under different wind speeds is shown in table 1.
TABLE 1 degree of imbalance in voltage at each node before optimization
Therefore, under different wind speeds, three-phase voltage unbalance phenomena exist in the collecting station and the subordinate wind field, and under the working condition of strong wind, namely, when the wind speed is over 12m/s, the negative sequence voltage unbalance degree exceeds the national standard requirement.
The split-phase power flow optimization algorithm provided by the invention can calculate the reactive power required to be compensated for each phase of the optimization node under different wind speeds, and the reactive power required to be compensated for each phase of the reactive power compensation device installed in the collection station is listed when the wind speeds are 9-13m/s respectively by taking the collection station as an example, as shown in the table 2.
Table 2 reactive power per phase compensation (MVA) for each node
After split-phase tide optimization, the voltage unbalance degree of the collecting station and the subordinate wind power plant under different wind speeds is shown in table 3.
TABLE 3 degree of imbalance in voltage at each node after optimization
The method for optimizing the split-phase power flow of the large-scale wind power collection grid can be seen to well reduce the negative sequence voltage unbalance degree at the collection station and the wind field outlet under different wind speeds.
The invention has the following advantages:
(1) The split-phase power flow optimization method for the large-scale wind power collection grid can effectively inhibit the problem of three-phase voltage unbalance of a collection station and a subordinate wind field under the working condition of strong wind, reduces the node voltage unbalance of all optimization nodes to below 2%, and meets the requirement that the voltage unbalance of the current national standard on PCC points (namely the collection station nodes in the invention) of a power system is not more than 2% and is not more than 4% in short time; the method can lower the unbalance degree of the PCC point negative sequence voltage of the large-scale wind power collection grid meeting national standard requirements under the working condition of small wind, improves the electric energy quality, improves the safety and stability of the large-scale wind power collection grid-connected system, and has a certain reference value for the development of large-scale wind power collection grid-connected system.
(2) The three-phase voltage unbalance degree is reduced, the node voltage reference value can be considered, and the deviation from the node voltage reference value can be as small as possible while the node voltage actual value is within a certain safety requirement range.
The reason for the advantages described above is:
(1) The split-phase power flow optimization model provided by the invention coordinates the split-phase dispatching and collecting station and the reactive compensation device of the subordinate wind field uniformly from the global angle, so that the voltage unbalance degree of all the optimization nodes can meet the requirement.
(2) The invention adds the target with minimum node voltage deviation into the split-phase power flow optimization model, so that the aims of inhibiting three-phase voltage unbalance and reducing the deviation between the actual value and the reference value of the node voltage can be simultaneously achieved.
The invention also provides a split-phase power flow optimizing system of the large-scale wind power collection grid, which comprises the following components:
the split-phase tide optimization model construction module is used for constructing a split-phase tide optimization model of the large-scale wind power collection grid; the split-phase power flow optimization model comprises an objective function and constraint conditions; the objective function comprises a minimum negative sequence voltage unbalance mean objective function and a minimum voltage deviation objective function;
the input data acquisition module is used for acquiring line parameters, load data and output power data of the wind power plant of the large-scale wind power collection grid at different wind speeds;
the compensation capacity obtaining module is used for solving the split-phase power flow optimization model by adopting a Pareto archive multi-target particle swarm optimization algorithm according to line parameters, load data and wind power plant output power data of the large-scale wind power collection grid at each wind speed, and obtaining the compensation capacity of each phase of the reactive power compensation device of the collection station and each wind power plant in the large-scale wind power collection grid at each wind speed.
The minimum negative sequence voltage unbalance mean value objective function isWherein F1 is the average value of the negative sequence voltage unbalance degree, U i1 Is the positive sequence voltage effective value of node i, < ->U i2 Is the negative sequence voltage effective value of node i, < -> And->The voltage effective values of a phase, b phase and c phase of the node i are respectively represented, alpha is a twiddle factor, and n is the total number of nodes;
the minimum voltage deviation objective function isWherein F2 is the voltage deviation,actual value of node voltage for the phi-th phase of node i, for example>The node voltage reference value for the phi-th phase of node i.
The constraint conditions include: equality constraints and inequality constraints.
The equality constraints include: node current equation constraint, branch voltage equation constraint and three-phase power flow equation constraint;
the node current equation is constrained to beWherein (1)>For phi phase current flowing on the branch connected between node i and node j, a, b and c are respectively a phase, b phase and c phase in three phases, < >>Represents the phi-phase active power of node j, < +.>Phi-phase reactive power, which represents node j, +.>Reactive power representing phase compensation of the reactive compensation means phi of node j +.>Phi-phase voltage, which represents node j, +.>The phi phase current flowing on the first branch connected with the node j is represented, and p is the number of branches connected between the node i and the node j;
the branch voltage equation is constrained toWherein (1)>Phi-phase voltage, which represents node i, +.>A branch impedance matrix representing a branch connected between node i and node j;
the three-phase tide equation is constrained as
Wherein P is i φ The phi-phase active power of the node i is N, the number of the nodes of the whole collecting power grid is V k m For m-phase voltage of node k, G ik,φm For the conductance between the phi and m phases of nodes i, k, theta ik,φm B is the voltage phase angle difference between phi and m phases of the node i and the node k ik,φm For susceptances between phi and m phases of node i, node k, ++>Phi-phase reactive power, which represents node i, +.>And representing the reactive power compensated by the node i reactive power compensation device phi.
Inequality constraints include: the negative sequence voltage unbalance degree constraint, the reactive compensation device compensates the capacity upper limit constraint and the capacity lower limit constraint and the node voltage safety constraint;
the negative sequence voltage unbalance degree is constrained as
The upper limit and the lower limit of the compensation capacity of the reactive compensation device are constrained as followsWherein (1)>The lower limit and the upper limit of reactive power which can be compensated by phi phase of the installed reactive power compensation device are respectively;
the voltage safety constraint of the node is thatWherein (1)>The phi-phase voltages at node i, respectively, allow a lower and an upper limit.
In the present specification, each embodiment is described in a progressive manner, and each embodiment is mainly described in a different point from other embodiments, and identical and similar parts between the embodiments are all enough to refer to each other. For the system disclosed in the embodiment, since it corresponds to the method disclosed in the embodiment, the description is relatively simple, and the relevant points refer to the description of the method section.
The principles and embodiments of the present invention have been described herein with reference to specific examples, the description of which is intended only to assist in understanding the methods of the present invention and the core ideas thereof; also, it is within the scope of the present invention to be modified by those of ordinary skill in the art in light of the present teachings. In view of the foregoing, this description should not be construed as limiting the invention.

Claims (8)

1. The split-phase power flow optimization method for the large-scale wind power collection grid is characterized by comprising the following steps of:
constructing a split-phase tide optimization model of a large-scale wind power collection grid; the split-phase power flow optimization model comprises an objective function and constraint conditions; the objective function comprises a minimum negative sequence voltage unbalance mean objective function and a minimum voltage deviation objective function;
acquiring line parameters, load data and output power data of a wind power plant of a large-scale wind power collection grid at different wind speeds;
according to line parameters, load data and output power data of the wind power plant of the large-scale wind power collection grid at each wind speed, solving the split-phase power flow optimization model by adopting a Pareto archive multi-objective particle swarm optimization algorithm to obtain compensation capacity of each phase of a collection station and a reactive power compensation device of each wind power plant in the large-scale wind power collection grid at each wind speed;
the minimum negative sequence voltage unbalance mean value objective function isWherein F is 1 Is the mean value of the negative sequence voltage unbalance degree, U i1 Is the positive of node iSequence voltage valid value, < >>U i2 Is the negative sequence voltage effective value of node i, < -> And->The voltage effective values of a phase, b phase and c phase of the node i are respectively represented, alpha is a twiddle factor, and n is the total number of nodes;
the minimum voltage deviation objective function isWherein F is 2 In order to be able to vary the voltage,actual value of node voltage for the phi-th phase of node i, for example>The node voltage reference value for the phi-th phase of node i.
2. The method for optimizing split-phase power flow of a large-scale wind power collection grid according to claim 1, wherein the constraint conditions comprise: equality constraints and inequality constraints.
3. The method for split-phase power flow optimization of a large-scale wind power collection grid according to claim 2, wherein the equality constraint comprises: node current equation constraint, branch voltage equation constraint and three-phase power flow equation constraint;
the node current equation is constrained to beWherein (1)>For phi phase current flowing on the branch connected between node i and node j, a, b and c are respectively a phase, b phase and c phase in three phases, < >>Represents the phi-phase active power of node j, < +.>Phi-phase reactive power, which represents node j, +.>Reactive power representing phase compensation of the reactive compensation means phi of node j +.>Phi-phase voltage, which represents node j, +.>The phi phase current flowing on the first branch connected with the node j is represented, and p is the number of branches connected between the node i and the node j;
the branch voltage equation is constrained to beWherein (1)>Phi-phase voltage, which represents node i, +.>A branch impedance matrix representing a branch connected between node i and node j;
the three-phase tide equation is constrained as
Wherein (1)>For phi-phase active power of node i, N is the number of nodes of the whole collecting power grid, ++>For m-phase voltage of node k, G ik,φm For the conductance between the phi and m phases of nodes i, k, theta ik,φm B is the voltage phase angle difference between phi and m phases of the node i and the node k ik,φm For susceptances between phi and m phases of node i, node k, ++>Phi-phase reactive power, which represents node i, +.>And representing the reactive power compensated by the node i reactive power compensation device phi.
4. A method of split-phase power flow optimization of a large-scale wind power collection grid according to claim 3, wherein the inequality constraints include: the negative sequence voltage unbalance degree constraint, the reactive compensation device compensates the capacity upper limit constraint and the capacity lower limit constraint and the node voltage safety constraint;
the negative sequence voltage unbalance degree constraint is that
The upper limit and the lower limit of the compensation capacity of the reactive compensation device are constrained as followsWherein (1)> The lower limit and the upper limit of reactive power which can be compensated by phi phase of the installed reactive power compensation device are respectively;
the voltage safety constraint of the node is thatWherein (1)>The phi-phase voltages at node i, respectively, allow a lower and an upper limit.
5. A split-phase power flow optimization system for a large-scale wind power collection grid, the system comprising:
the split-phase tide optimization model construction module is used for constructing a split-phase tide optimization model of the large-scale wind power collection grid; the split-phase power flow optimization model comprises an objective function and constraint conditions; the objective function comprises a minimum negative sequence voltage unbalance mean objective function and a minimum voltage deviation objective function;
the input data acquisition module is used for acquiring line parameters, load data and output power data of the wind power plant of the large-scale wind power collection grid at different wind speeds;
the compensation capacity obtaining module is used for solving the split-phase power flow optimization model by adopting a Pareto archive multi-objective particle swarm optimization algorithm according to line parameters, load data and wind power plant output power data of the large-scale wind power collection grid at each wind speed to obtain the compensation capacity of each phase of a collection station and reactive power compensation devices of each wind power plant in the large-scale wind power collection grid at each wind speed;
the minimum negative sequence voltage unbalance mean value objective function isWherein F is 1 Is the mean value of the negative sequence voltage unbalance degree, U i1 Is the positive sequence voltage effective value of node i, < ->U i2 Is the negative sequence voltage effective value of node i, < -> And->The voltage effective values of a phase, b phase and c phase of the node i are respectively represented, alpha is a twiddle factor, and n is the total number of nodes;
the minimum voltage deviation objective function isWherein F is 2 In order to be able to vary the voltage,actual value of node voltage for the phi-th phase of node i, for example>The node voltage reference value for the phi-th phase of node i.
6. The split-phase power flow optimization system of a large-scale wind power collection grid according to claim 5, wherein the constraint condition comprises: equality constraints and inequality constraints.
7. The large-scale wind power collection grid split-phase power flow optimization system of claim 6, wherein the equality constraint comprises: node current equation constraint, branch voltage equation constraint and three-phase power flow equation constraint;
the node current equation is constrained to beWherein (1)>For phi phase current flowing on the branch connected between node i and node j, a, b and c are respectively a phase, b phase and c phase in three phases, < >>Represents the phi-phase active power of node j, < +.>Phi-phase reactive power, which represents node j, +.>Reactive power representing phase compensation of the reactive compensation means phi of node j +.>Phi-phase voltage, which represents node j, +.>The phi phase current flowing on the first branch connected with the node j is represented, and p is the number of branches connected between the node i and the node j;
the branch voltage equation is constrained to beWherein (1)>Phi-phase voltage, which represents node i, +.>A branch impedance matrix representing a branch connected between node i and node j;
the three-phase tide equation is constrained as
Wherein (1)>For phi-phase active power of node i, N is the number of nodes of the whole collecting power grid, ++>For m-phase voltage of node k, G ik,φm For the conductance between the phi and m phases of nodes i, k, theta ik,φm B is the voltage phase angle difference between phi and m phases of the node i and the node k ik,φm For susceptances between phi and m phases of node i, node k, ++>Phi-phase reactive power, which represents node i, +.>And representing the reactive power compensated by the node i reactive power compensation device phi.
8. The large scale wind power collection grid split-phase power flow optimization system of claim 7, wherein the inequality constraint comprises: the negative sequence voltage unbalance degree constraint, the reactive compensation device compensates the capacity upper limit constraint and the capacity lower limit constraint and the node voltage safety constraint;
the negative sequence voltage unbalance degree constraint is that
The reactive powerThe upper limit and the lower limit of the compensation capacity of the compensation device are constrained asWherein (1)> The lower limit and the upper limit of reactive power which can be compensated by phi phase of the installed reactive power compensation device are respectively;
the voltage safety constraint of the node is thatWherein (1)>The phi-phase voltages at node i, respectively, allow a lower and an upper limit.
CN202111414460.5A 2021-11-25 2021-11-25 Large-scale wind power collection grid split-phase power flow optimization method and system Active CN114024338B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202111414460.5A CN114024338B (en) 2021-11-25 2021-11-25 Large-scale wind power collection grid split-phase power flow optimization method and system

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202111414460.5A CN114024338B (en) 2021-11-25 2021-11-25 Large-scale wind power collection grid split-phase power flow optimization method and system

Publications (2)

Publication Number Publication Date
CN114024338A CN114024338A (en) 2022-02-08
CN114024338B true CN114024338B (en) 2024-01-19

Family

ID=80066309

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202111414460.5A Active CN114024338B (en) 2021-11-25 2021-11-25 Large-scale wind power collection grid split-phase power flow optimization method and system

Country Status (1)

Country Link
CN (1) CN114024338B (en)

Families Citing this family (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN115425695B (en) * 2022-11-03 2023-02-03 国网四川省电力公司电力科学研究院 Power distribution network joint planning method suitable for distributed photovoltaic and energy storage
CN117096962B (en) * 2023-08-08 2024-04-09 国网浙江省电力有限公司宁海县供电公司 Photovoltaic-considered power grid dynamic reactive power compensation optimization method and system

Citations (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
GB8330819D0 (en) * 1983-11-18 1983-12-29 Ass Elect Ind Thyristor-controlled-reactor arrangement
WO2014173081A1 (en) * 2013-04-27 2014-10-30 国家电网公司 Optimization control method for reactive voltage of wind farm cluster
CN105870968A (en) * 2016-05-13 2016-08-17 河海大学 Three-phase imbalance reactive voltage control method metering system negative sequence voltage
CN106356838A (en) * 2016-08-31 2017-01-25 湘潭大学 Negative-sequence real-time optimized compensation method of RPC (railway power conditioner)
CN107508298A (en) * 2017-09-28 2017-12-22 合肥工业大学 A kind of micro-capacitance sensor unbalance voltage hierarchy optimization control method
CN107834539A (en) * 2017-10-13 2018-03-23 上海电力学院 A kind of excess load compensation method based on railway electric energy quality comprehensive treatment device
CN108539766A (en) * 2018-04-28 2018-09-14 广东电网有限责任公司 Three-phase imbalance virtual resistance optimization method based on coevolution
CN111900732A (en) * 2020-08-04 2020-11-06 四川大学 Wind power cluster access mode planning method considering wind power uncertainty and reactive power flow and voltage constraints
CN112564138A (en) * 2020-11-18 2021-03-26 国网山西省电力公司晋城供电公司 Three-phase unbalanced reactive power optimization method and system thereof

Patent Citations (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
GB8330819D0 (en) * 1983-11-18 1983-12-29 Ass Elect Ind Thyristor-controlled-reactor arrangement
WO2014173081A1 (en) * 2013-04-27 2014-10-30 国家电网公司 Optimization control method for reactive voltage of wind farm cluster
CN105870968A (en) * 2016-05-13 2016-08-17 河海大学 Three-phase imbalance reactive voltage control method metering system negative sequence voltage
CN106356838A (en) * 2016-08-31 2017-01-25 湘潭大学 Negative-sequence real-time optimized compensation method of RPC (railway power conditioner)
CN107508298A (en) * 2017-09-28 2017-12-22 合肥工业大学 A kind of micro-capacitance sensor unbalance voltage hierarchy optimization control method
CN107834539A (en) * 2017-10-13 2018-03-23 上海电力学院 A kind of excess load compensation method based on railway electric energy quality comprehensive treatment device
CN108539766A (en) * 2018-04-28 2018-09-14 广东电网有限责任公司 Three-phase imbalance virtual resistance optimization method based on coevolution
CN111900732A (en) * 2020-08-04 2020-11-06 四川大学 Wind power cluster access mode planning method considering wind power uncertainty and reactive power flow and voltage constraints
CN112564138A (en) * 2020-11-18 2021-03-26 国网山西省电力公司晋城供电公司 Three-phase unbalanced reactive power optimization method and system thereof

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
"考虑分布式电源随机特性的配电网重构研究";陈澎;《中国优秀硕士学位论文全文数据库工程科技Ⅱ辑》;第2019年卷(第05期);C042-907 *
"风电汇集地区电压不平衡特性分析及抑制策略";刘辉等;《电力***自动化》;第45卷(第14期);第132-139页 *

Also Published As

Publication number Publication date
CN114024338A (en) 2022-02-08

Similar Documents

Publication Publication Date Title
CN114024338B (en) Large-scale wind power collection grid split-phase power flow optimization method and system
CN110265991B (en) Distributed coordination control method for direct-current micro-grid
CN111446721B (en) Power distribution network voltage regulation control method based on transient voltage sensitivity
WO2021253368A1 (en) Coordinated control system and method of wind turbine and statcom for suppressing unbalanced voltage in dispersed wind farm
CN112467748A (en) Double-time-scale distributed voltage control method and system for three-phase unbalanced active power distribution network
Xiao et al. Optimal sizing and siting of soft open point for improving the three phase unbalance of the distribution network
Aydin et al. Comparative analysis of multi-criteria decision making methods for the assessment of optimal SVC location
CN109888817B (en) Method for carrying out position deployment and capacity planning on photovoltaic power station and data center
CN105896547B (en) A kind of bulk power grid hierarchical voltage control method under wind power integration
CN111082470B (en) Multi-target dynamic robust reconstruction method for power distribution network containing low wind speed distributed wind power
Wang et al. Pareto optimization of power system reconstruction using NSGA-II algorithm
CN111799846A (en) Distributed photovoltaic three-phase access capacity distribution calculation method
CN117035335A (en) Multi-stage energy storage and power transmission network collaborative planning method and system
CN107134783B (en) Bus voltage optimization adjustment method based on sensitivity rapid screening
Sun et al. Optimal voltage/var control of pv inverters in distribution network based on genetic algorithm
Rudion et al. MaWind-tool for the aggregation of wind farm models
CN115483715A (en) Virtual synchronous generator self-adaptive control method and system for centralized photovoltaic power station
CN110930263B (en) Medium-voltage distribution network short-circuit current calculation method containing photovoltaic power supply and induction motor based on black hole particle swarm algorithm
CN110829440B (en) Three-phase medium-low voltage integrated power distribution network voltage control method and system
CN104143826A (en) Reactive compensation method based on improved differential evolutionary algorithm and applied to electric power system with wind power plants
Maya et al. A generalised three phase power flow algorithm incorporating the uncertainty of Photo Voltaic (PV) source for unbalanced distribution network
CN109995088B (en) Safety control method for large-scale renewable energy source access to power grid
Fu et al. Control of PV systems for distribution network voltage regulation with communication delays
CN113742907A (en) Photovoltaic power station short-circuit current unified calculation method
CN114204613A (en) Reactive compensation method and system for offshore wind farm access power system

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