CN105449672A - Method for estimating total supply capability of 220KV loop ring network - Google Patents
Method for estimating total supply capability of 220KV loop ring network Download PDFInfo
- Publication number
- CN105449672A CN105449672A CN201510946257.0A CN201510946257A CN105449672A CN 105449672 A CN105449672 A CN 105449672A CN 201510946257 A CN201510946257 A CN 201510946257A CN 105449672 A CN105449672 A CN 105449672A
- Authority
- CN
- China
- Prior art keywords
- node
- population
- load
- power supply
- delta
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Pending
Links
- 238000000034 method Methods 0.000 title claims abstract description 32
- 238000004364 calculation method Methods 0.000 claims abstract description 38
- 239000013598 vector Substances 0.000 claims description 48
- 239000011159 matrix material Substances 0.000 claims description 30
- 238000002940 Newton-Raphson method Methods 0.000 claims description 16
- 238000013178 mathematical model Methods 0.000 claims description 11
- 230000009191 jumping Effects 0.000 claims description 10
- 230000005540 biological transmission Effects 0.000 claims description 8
- 230000035772 mutation Effects 0.000 claims description 7
- 238000012937 correction Methods 0.000 claims description 6
- 238000012545 processing Methods 0.000 claims description 6
- 230000003044 adaptive effect Effects 0.000 claims description 5
- 241000854350 Enicospilus group Species 0.000 claims description 3
- 241000894007 species Species 0.000 claims description 2
- 238000012795 verification Methods 0.000 abstract 1
- 230000008901 benefit Effects 0.000 description 4
- 230000003068 static effect Effects 0.000 description 4
- 230000007547 defect Effects 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 238000011156 evaluation Methods 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 238000006467 substitution reaction Methods 0.000 description 2
- 101100460704 Aspergillus sp. (strain MF297-2) notI gene Proteins 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000007796 conventional method Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 239000012634 fragment Substances 0.000 description 1
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
-
- G—PHYSICS
- G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
- G16Z—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS, NOT OTHERWISE PROVIDED FOR
- G16Z99/00—Subject matter not provided for in other main groups of this subclass
-
- 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 method for estimating total supply capability of a 220kV looped ring network. The method comprises the following steps of extracting a grid structure of the 220kV looped ring network, determining network node numbers and branch numbers of the grid structure, and numbering nodes and branches of the grid structure; reading a grid structure parameter, a load parameter and a voltage amplitude and a phase angle of a balance node; carrying out group expansion operation on a predetermined load population to obtain an initial load population, and carrying out variation and cross evolution operations on the load population; carrying out power flow calculation and verification on the 220kV looped ring network in a normal running mode and a (N-1) running of each branch; and updating the population according to self-adaptive degree, and outputting a load scheme of the total supply capability. According to the method, the total supply capability is estimated by taking a district as a minimum unit, and the total supply capability of the 220kV looped ring network is rapidly and accurately figured out when the (N-1) condition is met.
Description
Technical Field
The invention relates to the field of power supply capability evaluation of a power grid, in particular to a method for evaluating the maximum power supply capability of a 220kV loop ring network.
Background
Along with the economic development of China, the increase of power load is obviously accelerated, so that the power quality, the power supply capacity and the power supply reliability cannot meet the power requirements of users, and a plurality of power supply bottlenecks are formed. However, the urban power grid is basically built at present, and it is very difficult to obtain the site of a new substation and the underground passage of a new feeder from the planning and reconstruction of the system. Therefore, it becomes important to study the maximum power supply capacity of the grid without constructing underground tunnels for new power stations and feeders.
The maximum power supply capacity (TSC) of a power grid means that the power grid in a certain power supply area meets the N-1 safety criterion, and the maximum load supply capacity under the actual operation condition of the network is considered. The TSC reflects the safety margin of the load supply of the power grid to a certain extent, provides certain reference for operation and planning personnel of the power grid, and particularly after the introduction of the power market, the power grid personnel urgently need to utilize the existing power transmission network to transmit more electric energy, so that the cost is reduced to the maximum extent, and the benefit of the power grid is improved.
At present, the maximum power supply capacity is solved mainly through three stages: the first stage is mainly an evaluation stage of the maximum power supply capacity, a lot of assumptions are made in the first stage for modeling and calculation convenience, load flow calculation is directly ignored, accurate calculation cannot be performed, and the method is generally used for roughly evaluating the TSC of the power grid, such as a capacity-load ratio method. And in the second stage, the maximum power supply capacity of the power grid is solved by adopting a direct current power flow model or a linear programming method, and the defect that the influence of voltage and network loss on the TSC is ignored, so that a large error is generated. The third stage is a nonlinear model based on the power flow, which can accurately solve the power flow of the network and can consider the influence of voltage and network loss on the TSC, so that the calculation result is more reliable.
In addition, when the maximum power supply capacity of a power grid is solved by the conventional method, a complete power grid model or a simplified model is generally adopted for calculation, and if the complete power grid model and an accurate calculation method are adopted and N-1 static safety constraint calculation is considered, the calculation is slow, and even the problem of numerical calculation can occur; the simplified model causes the error of the result to be overlarge, and the application value is greatly reduced.
Disclosure of Invention
The invention aims to overcome the defects of the prior art, provides a method for evaluating the maximum power supply capacity of a 220kV loop network, establishes a nonlinear model based on alternating current power flow, considers the influence of voltage drop and network loss on the maximum power supply capacity of the 220kV loop network, and considers the static safety constraint of an electric power system N-1; the method has the advantages that the fragment area can be used as the minimum unit, the calculation scale is simplified, the calculation dimensionality is reduced, and the calculation precision is guaranteed while calculation is accelerated. In addition, the invention adopts the self-adaptive differential evolution algorithm embedded with the Newton Raphson method to solve the 220kV ring network, and has the advantages of high convergence speed and strong robustness.
The purpose of the invention is realized by the following technical scheme:
a method for evaluating the maximum power supply capacity of a 220kV loop network comprises the following steps:
s1, extracting a grid structure of the 220kV loop ring network, determining the number of network nodes and the number of branches of the grid structure, and numbering the nodes and the branches of the grid structure; the balance node of the 220kV ring-sleeved ring network is a node at the 220kV side of a main transformer of a 500kV transformer substation, and the PQ node is a load node in the 220kV ring-sleeved ring network;
s2, numberedGrid structure of 220kV ring-sleeved ring network, reading grid structure parameters, load parameters and voltage amplitude V of balance nodenAnd phase angle thetan;
S3, carrying out group expanding operation on the preset load population to obtain an initialized load population X, and carrying out variation and cross evolution operation on the initialized load population X through a self-adaptive differential evolution algorithm to respectively obtain variation populations XmutaAnd cross population Xcros;
S4, performing load flow calculation and check on the 220kV ring-sleeved ring network in a normal operation mode and in an N-1 operation mode of each branch by a Newton Raphson method;
s5, according to the self-adaptive degree FitUpdating the population X and outputting a maximum power supply capacity load scheme Sdi。
The step S3 specifically includes the following steps:
s301, judging whether the current iteration time T reaches the maximum value, if not, executing the step S302, and if the iteration time T reaches the set maximum iteration time T or the calculation precision meets the requirement, jumping to the step S309;
s302, carrying out group expanding operation on the population to obtain an initialized load species group, carrying out variation processing on the population X after the group expanding operation to obtain a variation population Xmuta: randomly selecting two individual vectors to generate a difference vector, adding the generated difference vector to another randomly selected vector to generate a variation vector Xmuta;
S303, carrying out mutation on the population XmutaPerforming cross treatment to obtain a cross population Xcros: the variation vector XmutaCrossing with the target vector X to generate a crossing vector Xcros;
S304, judging whether the traversal of the population is finished, if not, performing the step S305, otherwise, performing the processing of adding 1 to the iteration times, and jumping to the step S301;
S305. for individuals in the initial population X and the cross population XcrosThe individuals in the population X and the population X are calculated simultaneouslycrosRespectively substituting the individuals into the Newton Raphson method to carry out load flow calculation under normal operation and under each N-1 operation mode;
s306, checking the power flows in the normal operation mode and each N-1 operation mode, if all the checks are passed, performing the step S307, otherwise, adding 1 to the number of variables traversed by the population, and skipping to the step S304;
s307, taking a preset maximum power supply capacity mathematical model as a self-adaptive function, wherein the calculation formula is as follows:whereinRepresenting the active load of node i, i.e. group X or XcrosThe ith number in one dimension of (1);
s308, updating the population X: updating the individual with high self-adaptive degree into the population X;
and S309, outputting the maximum power supply capacity and the load scheme thereof.
The step S4 specifically includes the following steps:
s401, inputting original data: voltage V of 220kV side for inputting main transformer of 500kV transformer substationnAnd phase angle thetanParameters of the grid frame, voltage constraint and phase angle difference constraint;
s402, forming a node admittance matrix: forming a node admittance matrix under the normal operation mode and modifying the node admittance matrix under the normal operation mode according to each N-1 fault to form a node admittance matrix under each N-1 fault;
s403, calculating the power unbalance amount of each node Wherein, the delta P and the delta Q respectively refer to the deviation of the active power and the reactive power of the node, and whether the maximum power flow deviation meets the convergence condition is judged; if yes, jumping to step S406, and if not, performing step S404; the calculation formula of the power flow deviation is as follows:
wherein, Pis、QisGiving power to the ith node; vi、VjThe voltages of the ith node and the jth node respectively; gij、BijRespectively the conductance and susceptance of the branch from the node i to the node j; thetaijIs the phase angle difference between the node i and the node j;
s404, generating a power flow by the input variables and the existing node admittance matrix, and calculating a Jacobian matrix J:
wherein H is an n-1 order square matrix having elements ofN is an (N-1) × mth order matrix having elements ofK is an m × (n-1) order matrix with elements ofL is an m-th order square matrix having elements of
S405, solving a linear correction equation set WhereinObtaining correction quantities delta theta and delta V of the voltage amplitude and the phase angle of each node, obtaining the voltage and the phase angle of each node, and jumping to the step S403;
s406, calculating the power of all branches according to the following formula:
wherein i is the first node of the branch, j is the last node of the branch, the wave number represents the conjugate value of the complex number,is the conjugate value of the admittance value to ground to node i,is the conjugate of the admittance values between line node i to node j.
In step S307, the state variables of the mathematical model of the maximum power supply capacity include: apparent power S of each load pointdiAngle of power factorVoltage amplitude V of upper nodenAnd phase angle differencen;
The objective function of the mathematical model of the maximum power supply capacity is the maximum power supply capacity, and the sum of the active powers of the load nodes is calculated as the objective function according to the following formula:
the constraint conditions of the maximum power supply capacity mathematical model are as follows: and (3) load restraint:constraint of line transmission power:and (3) limiting the upper and lower limits of the node voltage:and about the phase angleBundling: non-viable cellsi-j|<|i-j|maxWhereinRespectively, a lower limit value and an upper limit value of the apparent power of the node i;respectively is the lower limit value and the upper limit value of the transmission power of the line from the node i to the node j branch;respectively, the lower and upper limit values of the voltage at node i.
In step S303, the interleaving process includes: randomly selecting two individual vectors to generate a difference vector, and adding the generated difference vector to another vector randomly selected to generate a variation vector according to the following formula:
wherein x isr1、xr2、xr3In a presentation population3 different individuals, said F0The scaling factor employs an adaptive strategy:F0u、F0lare respectively F0Upper and lower limits of ft1、ft2、ft3Are respectively asThe fitness of (2).
In step S302, the mutation process includes:
the variation vector is calculated according to the following formulaCrossing the target vectorGenerating cross vectors
the cross probability factor adopts a self-adaptive strategy:in the formula CRmin、CRmaxThe minimum cross probability factor and the maximum cross probability factor are respectively used, T is the current iteration frequency, and T is the maximum iteration frequency.
Compared with the prior art, the invention has the following advantages and beneficial effects:
the method for evaluating the maximum power supply capacity of the 220kV loop ring network establishes a nonlinear model based on alternating current power flow, and considers the influence of voltage drop and network loss on the maximum power supply capacity of the 220kV loop ring network; the maximum power supply capacity of the 220kV ring network is calculated by taking the zone as the minimum unit, a simplified wiring model of the maximum power supply capacity of the 220kV ring network is established, a large amount of redundant data is effectively simplified, the dimensionality of solution is reduced, the calculation speed is accelerated, and therefore the possibility is provided for more precise modeling and solution accuracy; the static safety constraint of the power system N-1 is considered, so that the maximum power supply capacity of the 220kV loop ring network can better meet the safety and reliability of the actual power grid in operation, and the calculation has practical significance; the self-adaptive strategy of differential evolution effectively accelerates the convergence of calculation, ensures the diversity of population and the robustness of calculation, and the Newton-Raphson method solution based on the alternating current model can accurately solve the load flow of the network, thereby ensuring the accuracy of the calculation result.
Drawings
Fig. 1 is a flowchart of a method for evaluating the maximum power supply capability of a 220kV loop ring network according to the present invention.
Fig. 2 is a schematic diagram of a typical looped network.
Detailed Description
The present invention will be described in further detail with reference to examples and drawings, but the present invention is not limited thereto.
As shown in fig. 1, a method for evaluating the maximum power supply capacity of a 220kV loop ring network includes the following steps:
s11, extracting the grid structure of the 220kV ring-sleeved ring network, determining the number of network nodes and the number of branches of the grid structure, and numbering the nodes and the branches of the grid structure.
S12, reading grid structure parameters, load parameters and voltage amplitude V of balance nodenAnd phase angle thetan。
S13, presetting population scale NpThe maximum iteration number T of the preset differential evolution algorithm and the minimum value F of the self-adaptive scaling factor0lAnd maximum value F0uMinimum value of adaptive cross probability factor CRminAnd maximum value CRmax(ii) a And carrying out group expanding operation on a preset load population to obtain an initialized load population X, and carrying out variation and cross evolution operation on the load population through the self-adaptive differential evolution algorithm to obtain a variation load population XmutaAnd cross population Xcros。
And S14, performing load flow calculation and check on the 220kV ring-sleeved ring network in a normal operation mode and in an N-1 operation mode of each branch by the Newton Raphson method.
And S15, updating the population according to the self-adaptive degree, and outputting a maximum power supply capacity load scheme.
The maximum power supply capacity in this embodiment means that a power grid in a certain power supply area meets the N-1 safety criterion, and the maximum load supply capacity under the actual operation condition of the network is considered. The maximum power supply capacity of the 220kV ring-sleeved ring network is calculated by adopting a self-adaptive differential evolution method embedded with a Newton Raphson method, the embedded Newton Raphson method is used for accurately solving the alternating current power flow of the ring-sleeved ring network, the self-adaptive differential evolution method is used for providing a large number of solutions for the Newton Raphson method, and the solutions solved by the Newton Raphson method are selected, so that the maximum power supply capacity of the ring-sleeved ring network can be quickly and accurately solved under the N-1 constraint.
Specifically, the method comprises the following steps:
firstly, a 220kV loop ring network structure is extracted. In the 220kV ring-sleeved ring network, a power supply point of an upper-level power grid is a node on the 220kV side of a main transformer of a 500kV transformer substation, and the node is regarded as a balance node in the 220kV ring-sleeved ring network; the other load nodes are generally referred to as PQ nodes. Determining the number of network nodes and the number of branches in the net rack, and numbering from left to right in sequence;
the 220kV ring-sleeved ring network is a first ring network composed of a 500kV substation and a plurality of 220kV substations, and a second ring network composed of some 220kV substations, or even more ring networks, as shown in fig. 2, a typical ring-sleeved ring network schematic diagram is shown.
And secondly, inputting simplified net rack data. Reading structural parameters of 220kV loop ring network frame (namely, each branch resistance RliReactance XliAnd a ground susceptance value BliAnd maximum ampacity Imax,li) Load parameter (maximum load value at each load point)And minimum load valueSd i ) And a voltage amplitude V given by the balancing nodenAnd phase angle thetan。
And thirdly, initializing the self-adaptive differential algorithm. The main steps include setting of parameters, e.g. population size NpMaximum iteration number T of differential evolution algorithm and minimum value F of self-adaptive scaling factor0lAnd maximum value F0uMinimum value of adaptive cross probability factor CRminAnd maximum value CRmaxSetting of (2); and carrying out group expanding operation on a preset load population to obtain an initialized load population X, and carrying out variation and cross evolution operation on the load population through the self-adaptive differential evolution algorithm to obtain a variation load population XmutaAnd cross population Xcros. The initialization of the node load randomly generates an initial population X according to the upper and lower limit values of the load, and the specific formula is as follows:in the formulaRespectively the maximum value and the minimum value of the rated load of the node j in the ith dimension;
the key steps of the self-adaptive differential evolution algorithm comprise:
1. cross-handling
Two individual vectors are randomly selected to generate a difference vector, and the generated difference vector is added to another vector randomly selected to generate a variation vector. The concrete formula is as follows:in the formula xr1、xr2、xr3Representing 3 different individuals in the population.
Wherein, the scaling factor of the mutation operation adopts a self-adaptive strategy, and the specific formula is as follows:F0u、F0lare respectively F0Upper and lower limits of ft1、ft2、ft3Are respectively asThe fitness of (2);
2. mutation treatment
The variation vectorCrossing the target vectorGenerating cross vectorsThe concrete formula is as follows: wherein rand (j) ∈ [0,1 ]]Random function of CRIs a cross probability factor.
The cross probability factor adopts a self-adaptive strategy, and the specific formula is as follows:
in the formula CRmin、CRmaxRespectively a minimum cross probability factor and a maximum cross probability factor, wherein T is the current iteration frequency, and T is the maximum iteration frequency;
3. selection process
If it is notIs adapted toIs superior toIs adapted toThen useInstead of the formerAnd is selected as the next generation; otherwiseAs the next generation. And adopting a greedy search strategy, and carrying out selection operation by taking the maximum load level as an objective function.
And fourthly, initializing the Newton Raphson method. The method mainly comprises the steps of setting the maximum iteration number T and the convergence progress of the Newton-Raphson method, and setting the node voltage ViPhase angle differenceijAnd the identification of PQ nodes and balancing nodes, etc.
The Newton Raphson method comprises the following steps of:
step 31, inputting original data: voltage V of 220kV side for inputting main transformer of 500kV transformer substationnAnd phase angle thetanParameters of the grid frame, voltage constraint and phase angle difference constraint;
step 32, forming a node admittance matrix: forming a node admittance matrix under the normal operation mode and modifying the node admittance matrix under the normal operation mode according to each N-1 fault to form a node admittance matrix under each N-1 fault;
step 33, calculating the power unbalance amount of each node (in the formula, delta P and delta Q respectively refer to the deviation of the active power and the reactive power of the node), and whether the maximum power flow deviation meets the convergence condition is judged; if yes, jumping to step 36, if not, proceeding to step 34; the calculation formula of the power flow deviation is as follows:
wherein, Pis、QisGiving power to the ith node; vi、VjThe voltages of the ith node and the jth node respectively; gij、BijRespectively the conductance and susceptance of the branch from the node i to the node j; thetaijIs the phase angle difference between node i and node j)
Step 34, generating a power flow by the input variables and the existing node admittance matrix, and calculating a Jacobian matrix J:
wherein H is an n-1 order square matrix having elements ofN is an (N-1) × mth order matrix having elements ofK is an m × (n-1) order matrix with elements ofL is an m-th order square matrix having elements of
Step 35, solving the linear correction equation set WhereinObtaining the correction quantities delta theta and delta V of the voltage amplitude and the phase angle of each node, obtaining the voltage and the phase angle of each node, and skipping to the step 33;
step 36, calculating the power of all branches according to the following formula:
wherein i is the first node of the branch, j is the last node of the branch, the wave number represents the conjugate value of the complex number,is the conjugate value of the admittance value to ground to node i,is the conjugate of the admittance values between line node i to node j.
And fifthly, optimizing and realizing.
Step 21, judging whether the current iteration time T reaches the maximum value, if not, performing step 22, and if the iteration time T reaches the set maximum iteration time T or the calculation precision meets the requirement, skipping to step 29;
step 22, performing population expansion operation on the population to obtain an initialized load population, and performing variation processing on the population X after population expansion to obtain a variation population Xmuta: randomly selecting two individual vectors to generate a difference vector, adding the generated difference vector to another randomly selected vector to generate a variation vector XmutaThe concrete formula isIn the formula xr1、xr2、xr3Representing 3 different individuals in the population, F0The scaling factor of the mutation operation adopts an adaptive strategy, and the specific formula is as follows:F0u、F0lare respectively F0Upper and lower limits of ft1、ft2、ft3Are respectively asThe fitness of (2);
step 23, changing the population XmutaPerforming cross treatment to obtain a cross population Xcros: the variation vector XmutaCrossing with the target vector X to generate a crossing vector XcrosThe concrete formula is Wherein rand (j) ∈ [0,1 ]]Random function of CRThe method adopts a self-adaptive strategy as a cross probability factor, and the specific formula is as follows:in the formula CRmin、CRmaxRespectively a minimum cross probability factor and a maximum cross probability factor, wherein T is the current iteration frequency, and T is the maximum iteration frequency;
step 24, judging whether the traversal of the population is finished, if not, performing step 25, otherwise, performing 1 addition processing on the iteration times, and jumping to the step 21;
step 25, individual in the initial population X and the cross population XcrosThe individuals in the population X and the population X are calculated simultaneouslycrosRespectively substituting the individuals into the Newton Raphson method to carry out load flow calculation under normal operation and under each N-1 operation mode;
step 26, checking the power flow in the normal operation mode and each N-1 operation mode, if all checks pass, such as voltage constraint, phase angle constraint, constraint of line transmission power and the like, performing step 27, otherwise, adding 1 to the number of variables traversed by the population, and skipping to step 24;
step 27, taking a preset maximum power supply capacity mathematical model as a self-adaptive function, and calculating a formula as follows:whereinRepresenting the active load of node i, i.e. group X or XcrosThe ith number in one dimension of (1);
step 28, updating the population X: updating the individual with high self-adaptive degree into the population X;
and step 29, outputting the maximum power supply capacity and the load scheme thereof.
Wherein the maximum power supply capability mathematical model may be described as: and taking the maximum power supply capacity as a target function, taking the N-1 safety criterion into consideration according to the definition of the maximum power supply capacity, and taking the actual conditions of network operation into consideration, wherein the actual conditions comprise the constraints of main transformer capacity, voltage and phase angle constraints, line overload capacity and the like, so as to obtain the maximum power supply capacity. The specific mathematical model is as follows:
the state variables are as follows: including apparent power S of each load pointdi(ii) a Voltage amplitude V of upper nodenPhase angle differencen。
An objective function: taking the maximum power supply capacity as an objective function, namely taking the sum of the active power of each load node as the objective function
Considering the actual operation condition of a 220kV power grid, the power factors of all load points are consideredAre all set to 0.98.
The constraint conditions include:
A. load restraint
In the formula,S dithe value is 30 percent of the sum of the rated capacities of all main transformers in the transformer substation,the value is that when one main transformer with the maximum capacity in the transformer substation quits operation, the other main transformers are overloaded by 20 percent.
B. Line transmission power constraint
Wherein, when k is 0, it represents normal operation mode; when k is>When 0, the fault of the k line N-1 is shown;indicating the maximum transmission capacity of the node i to node j line in the k-th state.
C. Upper and lower limit constraints of node voltage
In the formula,the minimum voltage of the node i in the kth state is generally 0.95;the maximum voltage of the node i in the kth state is generally 1.05.
D. Upper and lower constraint of phase angle
In the formula,which represents the maximum value of the phase angle difference between the node i and the node j in the k-th state.
The invention relates to a method for evaluating the maximum power supply capacity of a 220kV loop ring network, which establishes a nonlinear model based on alternating current power flow and considers the influence of voltage drop and network loss on the maximum power supply capacity of the 220kV loop ring network; the maximum power supply capacity of the 220kV ring network is calculated by taking the zone as the minimum unit, a simplified wiring model of the maximum power supply capacity of the 220kV ring network is established, a large amount of redundant data is effectively simplified, the dimensionality of solution is reduced, the calculation speed is accelerated, and therefore the possibility is provided for more precise modeling and solution accuracy; the static safety constraint of the power system N-1 is considered, so that the maximum power supply capacity of the 220kV loop ring network can better meet the safety and reliability of the actual power grid in operation, and the calculation has practical significance; the self-adaptive strategy of differential evolution effectively accelerates the convergence of calculation, ensures the diversity of population and the robustness of calculation, and the Newton-Raphson method solution based on the alternating current model can accurately solve the load flow of the network, thereby ensuring the accuracy of the calculation result.
The above embodiments are preferred embodiments of the present invention, but the present invention is not limited to the above embodiments, and any other changes, modifications, substitutions, combinations, and simplifications which do not depart from the spirit and principle of the present invention should be construed as equivalents thereof, and all such changes, modifications, substitutions, combinations, and simplifications are intended to be included in the scope of the present invention.
Claims (6)
1. A method for evaluating the maximum power supply capacity of a 220kV loop network is characterized by comprising the following steps:
s1, extracting a grid structure of the 220kV loop ring network, determining the number of network nodes and the number of branches of the grid structure, and numbering the nodes and the branches of the grid structure; the balance node of the 220kV ring-sleeved ring network is a node at the 220kV side of a main transformer of a 500kV transformer substation, and the PQ node is a load node in the 220kV ring-sleeved ring network;
s2, reading parameters of the net rack structure of the numbered 220kV ring sleeve ring networkLoad parameter and voltage amplitude V of balance nodenAnd phase angle thetan;
S3, carrying out group expanding operation on the preset load population to obtain an initialized load population X, and carrying out variation and cross evolution operation on the initialized load population X through a self-adaptive differential evolution algorithm to respectively obtain variation populations XmutaAnd cross population Xcros;
S4, performing load flow calculation and check on the 220kV ring-sleeved ring network in a normal operation mode and in an N-1 operation mode of each branch by a Newton Raphson method;
s5, according to the self-adaptive degree FitUpdating the population X and outputting a maximum power supply capacity load scheme Sdi。
2. The method for evaluating the maximum power supply capability of the 220kV loop network according to claim 1, wherein the step S3 specifically comprises the following steps:
s301, judging whether the current iteration time T reaches the maximum value, if not, executing the step S302, and if the iteration time T reaches the set maximum iteration time T or the calculation precision meets the requirement, jumping to the step S309;
s302, carrying out group expanding operation on the population to obtain an initialized load species group, carrying out variation processing on the population X after the group expanding operation to obtain a variation population Xmuta: randomly selecting two individual vectors to generate a difference vector, adding the generated difference vector to another randomly selected vector to generate a variation vector Xmuta;
S303, carrying out mutation on the population XmutaPerforming cross treatment to obtain a cross population Xcros: the variation vector XmutaCrossing with the target vector X to generate a crossing vector Xcros;
S304, judging whether the traversal of the population is finished, if not, performing the step S305, otherwise, performing the processing of adding 1 to the iteration times, and jumping to the step S301;
s305, carrying out comparison on individuals in the initial population X and the cross population XcrosI.e. the individuals in the population X are calculated simultaneouslyBody and population XcrosRespectively substituting the individuals into the Newton Raphson method to carry out load flow calculation under normal operation and under each N-1 operation mode;
s306, checking the power flows in the normal operation mode and each N-1 operation mode, if all the checks are passed, performing the step S307, otherwise, adding 1 to the number of variables traversed by the population, and skipping to the step S304;
s307, taking a preset maximum power supply capacity mathematical model as a self-adaptive function, wherein the calculation formula is as follows:whereinRepresenting the active load of node i, i.e. group X or XcrosThe ith number in one dimension of (1);
s308, updating the population X: updating the individual with high self-adaptive degree into the population X;
and S309, outputting the maximum power supply capacity and the load scheme thereof.
3. The method for evaluating the maximum power supply capability of the 220kV loop network according to claim 2, wherein the step S4 specifically comprises the following steps:
s401, inputting original data: voltage V of 220kV side for inputting main transformer of 500kV transformer substationnAnd phase angle thetanParameters of the grid frame, voltage constraint and phase angle difference constraint;
s402, forming a node admittance matrix: forming a node admittance matrix under the normal operation mode and modifying the node admittance matrix under the normal operation mode according to each N-1 fault to form a node admittance matrix under each N-1 fault;
s403, calculating the power unbalance amount of each node Wherein, the delta P and the delta Q respectively refer to the deviation of the active power and the reactive power of the node, and whether the maximum power flow deviation meets the convergence condition is judged; if yes, jumping to step S406, and if not, performing step S404; the calculation formula of the power flow deviation is as follows:
wherein, Pis、QisGiving power to the ith node; vi、VjThe voltages of the ith node and the jth node respectively; gij、BijRespectively the conductance and susceptance of the branch from the node i to the node j; thetaijIs the phase angle difference between the node i and the node j;
s404, generating a power flow by the input variables and the existing node admittance matrix, and calculating a Jacobian matrix J:
whereinH is an n-1 order square matrix with elements ofN is an (N-1) × mth order matrix having elements ofK is an m × (n-1) order matrix with elements ofL is an m-th order square matrix having elements of
S405, solving a linear correction equation set WhereinObtaining correction quantities delta theta and delta V of the voltage amplitude and the phase angle of each node, obtaining the voltage and the phase angle of each node, and jumping to the step S403;
s406, calculating the power of all branches according to the following formula:
wherein,i is the first node of the branch, j is the last node of the branch, the wave number represents the conjugate value of the complex number,is the conjugate value of the admittance value to ground to node i,is the conjugate of the admittance values between line node i to node j.
4. The method for calculating the maximum power supply capacity of the 220kV loop network according to claim 3, wherein in step S307, the state variables of the mathematical model of the maximum power supply capacity include: apparent power S of each load pointdiAngle of power factorVoltage amplitude V of upper nodenAnd phase angle differencen;
The objective function of the mathematical model of the maximum power supply capacity is the maximum power supply capacity, and the sum of the active powers of the load nodes is calculated as the objective function according to the following formula:
the constraint conditions of the maximum power supply capacity mathematical model are as follows: and (3) load restraint:constraint of line transmission power:and (3) limiting the upper and lower limits of the node voltage:and phase angle upper and lower constraints: non-viable cellsi-j|<|i-j|maxWhereinS di、Respectively, a lower limit value and an upper limit value of the apparent power of the node i;S ij、respectively is the lower limit value and the upper limit value of the transmission power of the line from the node i to the node j branch; vi、ViRespectively, the lower and upper limit values of the voltage at node i.
5. The method for detecting the maximum power supply capacity of the 220KV loop ring network according to claim 2, wherein in step S303, the crossing process is: randomly selecting two individual vectors to generate a difference vector, and adding the generated difference vector to another vector randomly selected to generate a variation vector according to the following formula:
wherein x isr1、xr2、xr3In a presentation population3 different individuals, said F0The scaling factor employs an adaptive strategy:F0u、F0lare respectively F0Upper and lower limits of ft1、ft2、ft3Are respectively asThe fitness of (2).
6. The method for evaluating the maximum power supply capability of the 220kV loop network according to claim 2, wherein in step S302, the mutation process comprises:
the variation vector is calculated according to the following formulaCrossing the target vectorGenerating cross vectors
the cross probability factor adopts a self-adaptive strategy:in the formula CRmin、CRmaxThe minimum cross probability factor and the maximum cross probability factor are respectively used, T is the current iteration frequency, and T is the maximum iteration frequency.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510946257.0A CN105449672A (en) | 2015-12-16 | 2015-12-16 | Method for estimating total supply capability of 220KV loop ring network |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510946257.0A CN105449672A (en) | 2015-12-16 | 2015-12-16 | Method for estimating total supply capability of 220KV loop ring network |
Publications (1)
Publication Number | Publication Date |
---|---|
CN105449672A true CN105449672A (en) | 2016-03-30 |
Family
ID=55559590
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510946257.0A Pending CN105449672A (en) | 2015-12-16 | 2015-12-16 | Method for estimating total supply capability of 220KV loop ring network |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105449672A (en) |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105762826A (en) * | 2016-04-15 | 2016-07-13 | 华北电力大学 | AC-DC system containing VSC-HVDC total transfer capability calculation method |
CN107528312A (en) * | 2017-07-14 | 2017-12-29 | 同济大学 | A kind of power system state estimation method |
CN108899896A (en) * | 2018-06-28 | 2018-11-27 | 华南理工大学 | A kind of evaluation of power supply capability method based on improvement benders decomposition method |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2012161378A1 (en) * | 2011-05-25 | 2012-11-29 | (주)파워이십일 | Method for estimating the status of a power distribution system |
CN104484555A (en) * | 2014-11-26 | 2015-04-01 | 广州电力设计院 | Method for evaluating maximum power supply capability of 220kV self-healing looped network |
-
2015
- 2015-12-16 CN CN201510946257.0A patent/CN105449672A/en active Pending
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2012161378A1 (en) * | 2011-05-25 | 2012-11-29 | (주)파워이십일 | Method for estimating the status of a power distribution system |
CN104484555A (en) * | 2014-11-26 | 2015-04-01 | 广州电力设计院 | Method for evaluating maximum power supply capability of 220kV self-healing looped network |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105762826A (en) * | 2016-04-15 | 2016-07-13 | 华北电力大学 | AC-DC system containing VSC-HVDC total transfer capability calculation method |
CN107528312A (en) * | 2017-07-14 | 2017-12-29 | 同济大学 | A kind of power system state estimation method |
CN107528312B (en) * | 2017-07-14 | 2020-04-28 | 同济大学 | Power system state estimation method |
CN108899896A (en) * | 2018-06-28 | 2018-11-27 | 华南理工大学 | A kind of evaluation of power supply capability method based on improvement benders decomposition method |
CN108899896B (en) * | 2018-06-28 | 2020-07-28 | 华南理工大学 | Power supply capacity evaluation method based on improved benders decomposition method |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
EP3780307B1 (en) | Online dynamic decision-making method and system for unit restoration | |
CN106972504B (en) | Interval reactive power optimization method based on genetic algorithm | |
US20150112498A1 (en) | Methods of Patel Loadflow Computation for Electrical Power System | |
CN110034581B (en) | Interval electrical medium number vulnerability assessment method for power system under wind power grid-connected condition | |
CN104992266A (en) | Method of determining power grid node importance degree and system thereof | |
Ahmadi et al. | Two approaches to transient stability-constrained optimal power flow | |
Hao et al. | Optimising location of unified power flow controllers by means of improved evolutionary programming | |
US20210049228A1 (en) | Methods of Patel Loadflow Computation for Electrical Power System | |
Sun et al. | A network reconfiguration approach for power system restoration based on preference-based multiobjective optimization | |
CN104484555B (en) | The method of assessment 220kV self-healing looped network net capability | |
CN110417011A (en) | A kind of online dynamic secure estimation method based on mutual information Yu iteration random forest | |
CN104901306A (en) | Power grid operation safety margin calculating method taking cascading faults into consideration | |
CN106655227A (en) | SOP-based active power distribution network feeder load balancing method | |
CN106295911A (en) | A kind of grid branch parameter evaluation method based on chromatographic assays | |
CN105449672A (en) | Method for estimating total supply capability of 220KV loop ring network | |
CN104037760B (en) | Anticipated accident selection method for electric power system with uncertain injection power | |
CN105490288B (en) | A kind of 220kV electric network reactive compensations Optimal Configuration Method | |
Choi et al. | Cost-effective retrofits of power grids based on critical cascading failure scenarios identified by multi-group non-dominated sorting genetic algorithm | |
CN104767207B (en) | Voltage stabilization prevention and control method based on generator reactive deposit sensitivity | |
CN104484832B (en) | The method of assessment 220KV handle net net capability | |
CN105529703A (en) | Urban power net reconstruction planning method based on power supply capability bottleneck analysis | |
CN105811429A (en) | Locating and sizing method of static synchronous series compensator | |
CN104410065B (en) | A kind of Multiobjective Decision Making Method of receiving end electrical network limiting short-circuit current | |
CN105576653B (en) | A kind of 220kV sections power network power supply capacity optimization method | |
CN104578050A (en) | Transient stability strongly-correlated power transmission section identification method for power grid |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
WD01 | Invention patent application deemed withdrawn after publication |
Application publication date: 20160330 |
|
WD01 | Invention patent application deemed withdrawn after publication |