CN104484832A - Method for evaluating total supplying capability of 220KV Lashou net - Google Patents
Method for evaluating total supplying capability of 220KV Lashou net Download PDFInfo
- Publication number
- CN104484832A CN104484832A CN201410706099.7A CN201410706099A CN104484832A CN 104484832 A CN104484832 A CN 104484832A CN 201410706099 A CN201410706099 A CN 201410706099A CN 104484832 A CN104484832 A CN 104484832A
- Authority
- CN
- China
- Prior art keywords
- node
- net
- population
- load
- handle
- 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.)
- Granted
Links
- 238000000034 method Methods 0.000 title claims abstract description 81
- 238000012545 processing Methods 0.000 claims abstract description 5
- 230000003044 adaptive effect Effects 0.000 claims description 28
- 239000013598 vector Substances 0.000 claims description 28
- 239000011159 matrix material Substances 0.000 claims description 10
- 230000005540 biological transmission Effects 0.000 claims description 9
- 238000013178 mathematical model Methods 0.000 claims description 7
- 230000015572 biosynthetic process Effects 0.000 claims description 6
- 238000005755 formation reaction Methods 0.000 claims description 6
- 230000035772 mutation Effects 0.000 claims description 6
- 230000001131 transforming effect Effects 0.000 claims description 6
- 230000005283 ground state Effects 0.000 claims description 4
- 238000004364 calculation method Methods 0.000 claims description 3
- 230000008649 adaptation response Effects 0.000 claims description 2
- 238000012937 correction Methods 0.000 claims description 2
- 230000009191 jumping Effects 0.000 claims description 2
- 230000008929 regeneration Effects 0.000 claims description 2
- 238000011069 regeneration method Methods 0.000 claims description 2
- 238000005316 response function Methods 0.000 claims description 2
- 238000012795 verification Methods 0.000 claims description 2
- 230000000694 effects Effects 0.000 description 6
- 238000010586 diagram Methods 0.000 description 2
- 230000005611 electricity Effects 0.000 description 2
- 238000011160 research Methods 0.000 description 1
- 238000012546 transfer Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q10/00—Administration; Management
- G06Q10/06—Resources, workflows, human or project management; Enterprise or organisation planning; Enterprise or organisation modelling
- G06Q10/063—Operations research, analysis or management
- G06Q10/0637—Strategic management or analysis, e.g. setting a goal or target of an organisation; Planning actions based on goals; Analysis or evaluation of effectiveness of goals
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q50/00—Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
- G06Q50/06—Energy or water supply
Landscapes
- Business, Economics & Management (AREA)
- Human Resources & Organizations (AREA)
- Engineering & Computer Science (AREA)
- Economics (AREA)
- Strategic Management (AREA)
- General Physics & Mathematics (AREA)
- General Business, Economics & Management (AREA)
- Health & Medical Sciences (AREA)
- Theoretical Computer Science (AREA)
- Marketing (AREA)
- Entrepreneurship & Innovation (AREA)
- Educational Administration (AREA)
- Tourism & Hospitality (AREA)
- Physics & Mathematics (AREA)
- Operations Research (AREA)
- Quality & Reliability (AREA)
- Development Economics (AREA)
- Game Theory and Decision Science (AREA)
- Public Health (AREA)
- Water Supply & Treatment (AREA)
- General Health & Medical Sciences (AREA)
- Primary Health Care (AREA)
- Supply And Distribution Of Alternating Current (AREA)
Abstract
The invention provides a method for evaluating the total supplying capability of a 220KV Lashou net. The method comprises the following steps of obtaining a space truss structure of the 220KV Lashou net, determining the network node number and the branch number of the space truss structure, and numbering the nodes and the branches of the space truss structure; reading the parameter and the load parameter of the space truss structure as well as the voltage amplitude and the phase angle of a balance node for the numbered space truss structure of the 220KV Lashou net; according to the constraint of the preset node voltage Vi and the phase angle difference deltaij as well as a PQnode and a balance node, obtaining an initialized Newton-Ralfson's method; obtaining an initialized self-adaptive differential evolutionary algorithm; according to the Newton-Ralfson's method, solving the tidal current of the 220KV Lashou net, and processing the tidal current by the self-adaptive differential evolutionary algorithm to obtain the total supplying capability of 220KV Lashou net under a N-1 constraint. According to the method, a district can be taken as a minimum unit to evaluate the total supplying capability, and the total supplying capability of the 220KV Lashou net can be quickly and accurately obtained when the N-1 condition is satisfied.
Description
Technical field
The present invention relates to mains supply enabling technology field, particularly relate to a kind of method assessing 220KV handle net net capability.
Background technology
Along with China's expanding economy, living standards of the people improve constantly, and the demand for electric energy also constantly increases thereupon.The growth of electric load is obviously accelerated, and causes the quality of power supply, power supply capacity and power supply reliability to can not meet the electricity needs of user, defines bottleneck of much powering.And urban distribution network has substantially built up at present, if the underpass wanting the site and new feeder line obtaining new transformer station from the planning and reconstruction of system is very difficult.Therefore, when not building the underpass in new power station and feeder line, the net capability of research electrical network becomes very important.
Electrical network net capability (Total Supplying Capacity, TSC) refers to that in certain power supply area, electrical network meets N-1 safety criterion, and considers the peak load deliverability under network practical operation situation.Solve the linear law of planning of common method of electrical network net capability, interior point method, hit-and-miss method and peak load method of multiplicity etc.Compared with net capability, available transmission capacity (Available Transfer Capacity, ATC) refers on existing transmission of electricity Contract basis, remaining in actual physics power transmission network, to can be used for business use transmission capacity.TSC it is emphasised that electrical network to meet under certain constraint can with peak load, what ATC paid close attention to is deduct basic trend and adequate allowance on the basis of maximum transmitted ability in power transmission network after, the maximum power also can transmitted.
At present, the linear law of planning of the common method of electrical network net capability, hit-and-miss method and peak load method of multiplicity etc. are determined.
Linear programming technique, based on DC flow model, does not consider impact that is idle and voltage, so can affect accuracy and the validity of result; Hit-and-miss method solution procedure is quite consuming time, and the accuracy of result is difficult to ensure; Peak load method of multiplicity solving speed is fast, but it supposes that each node load all increases with same ratio, and accuracy can decrease.
Summary of the invention
Based on this, the invention provides a kind of method assessing 220KV handle net net capability, can be that least unit is assessed net capability with section, meet under N-1 condition, determine the net capability of 220KV handle net fast and accurately.
Assess a method for 220KV handle net net capability, comprise the steps:
Obtain the grid structure of 220KV handle net, determine number of network node and the circuitry number of described grid structure, and the node of described grid structure and branch road are numbered; Wherein, the balance node of 220kV handle net is the 220kV side of 500kV transforming plant main transformer, and PQ node is the load bus in 220kV handle net;
To the grid structure of the 220KV handle net after numbering, read the voltage magnitude V of grid structure parameter, load parameter and balance node
nand phase angle theta
n;
According to the node voltage V preset
iwith phase angle difference δ
ijconstraint, and described PQ node and balance node, obtain the Newton-Laphson method after initialization;
According to the initialization of each node load, the initial population preset carried out to population X, default differential evolution algorithm maximum iteration time Gm, self adaptive pantographic factor minimum value F0 that group expanding process obtains
minwith maximal value F0
max, adaptive crossover mutation factor minimum value CR
minwith maximal value CR
max, obtain the adaptive differential evolution algorithm after initialization; Wherein, according to following formula, the mean value initialization of node load being got to each load upper lower limit value generates initial population:
in formula
s i be respectively maximal value and the minimum value of the nominal load of node i;
Solve the trend of described 220KV handle net according to described Newton-Laphson method, by described adaptive differential evolution algorithm, described trend is processed, obtain the net capability that the described 220KV handle under N-1 constraint is netted.
The method of above-mentioned assessment 220KV handle net net capability, embedded Newton-Laphson method in adaptive differential evolvement method, the effect of embedded Newton-Laphson method is the trend of Exact Solution handle net, the effect of adaptive differential evolution algorithm processes the trend that Newton-Laphson method is obtained fast, to reach under N-1 constraint, obtain the net capability of handle net fast and accurately; Can be that least unit calculates net capability with section, the method has the features such as good stability, precision is high, ability of searching optimum is strong.
Accompanying drawing explanation
Fig. 1 is the method schematic flow sheet in one embodiment that the present invention assesses 220KV handle net net capability.
Fig. 2 is a kind of typical handle net schematic diagram.
Embodiment
Below in conjunction with embodiment and accompanying drawing, the present invention is described in further detail, but embodiments of the present invention are not limited thereto.
As shown in Figure 1, be a kind of method assessing 220KV handle net net capability of the present invention, comprise the steps:
S11, obtain the grid structure of 220KV handle net, determine number of network node and the circuitry number of described grid structure, and the node of described grid structure and branch road are numbered; Wherein, the balance node of 220kV handle net is the 220kV side of 500kV transforming plant main transformer, and PQ node is the load bus in 220kV handle net;
S12, the grid structure of netting the 220KV handle after numbering, read the voltage magnitude V of grid structure parameter, load parameter and balance node
nand phase angle theta
n;
The node voltage V that S13, basis are preset
iwith phase angle difference δ
ijconstraint, and described PQ node and balance node, obtain the Newton-Laphson method after initialization;
S14, according to the initialization of each node load, the initial population preset carried out to population X, default differential evolution algorithm maximum iteration time Gm, self adaptive pantographic factor minimum value F0 that group expanding process obtains
minwith maximal value F0
max, adaptive crossover mutation factor minimum value CR
minwith maximal value CR
max, obtain the adaptive differential evolution algorithm after initialization; Wherein, according to following formula, the mean value initialization of node load being got to each load upper lower limit value generates initial population:
in formula
s i be respectively maximal value and the minimum value of the nominal load of node i;
S15, to solve described 220KV according to described Newton-Laphson method and to shake hands the trend of net, by described adaptive differential evolution algorithm, described trend is processed, obtain the net capability that the lower described 220KV handle of N-1 constraint is netted;
Net capability, refers to that in certain power supply area, electrical network meets N-1 safety criterion, and considers the peak load deliverability under network practical operation situation.The calculating 220KV handle net net capability method of the present embodiment, it is the adaptive differential evolvement method of embedded Newton-Laphson method, the effect of embedded Newton-Laphson method is the trend of Exact Solution handle net, the effect of adaptive differential evolvement method processes the trend that Newton-Laphson method is obtained fast, to reach under N-1 constraint, obtain the net capability of handle net fast and accurately.
Concrete:
One, 220kV typical case handle web frame is refined.In 220kV handle net, its higher level's mains supply point is generally two, is generally the 220kV side of 500kV transforming plant main transformer, is generally regarded as balance node; Other load buses generally regard PQ node as.Determine number of network node and circuitry number in rack, and from left to right number successively;
Wherein, described handle net, refers to a kind of topological structure of the electrical network formed in hand-in-hand mode, as shown in Figure 2, shows a kind of typical handle net schematic diagram.
Two, the rack data of input simplification.Read 220kV handle grid structure parameter (i.e. each branch resistance R
i, reactance X
i, susceptance value B over the ground
iand maximum carrying capacity I
max), load parameter (the peak load value of each load point
and minimal negative charge values
sd i ) and the given voltage magnitude V of balance node
nand phase angle theta
n.
Three, Newton-Laphson method initialization.Key step comprises node voltage V
iwith phase angle difference δ
ijthe setting of constraint and the identification etc. to PQ node and balance node.
At handle net owing to there are two balance node, simultaneously in order to ensure the accuracy calculated, Newton-Laphson method is adopted to calculate trend.In computation process, the load iteration of 3 220KV just can restrain for 10 times, and the error calculated is almost nil.Use the step of Newton-Laphson method calculating trend as follows:
Step 21, input raw data: as inputted voltage Vn and phase angle theta n, the parameter of rack and the constraint of voltage constraint and phase angle difference of the 220kV side of two 500KV transforming plant main transformers;
Step 22, formation bus admittance matrix, and according to N-1 fault criteria, bus admittance matrix is modified;
Step 23, calculate each node power amount of unbalance
(in formula, Δ P, Δ Q distinguish active power deviation and the reactive power deviation of dactylus point), judges whether strength of current deviation meets the condition of convergence; As met, then jumping to step 26, if do not met, then carry out step 24; Wherein, the computing formula of trend deviation is as follows:
Step 24, generate Load flow calculation Jacobi matrix J by the variable inputted and existing bus admittance matrix:
Wherein, H is (n-1) rank square formations, and its element is
n is (n-1) × m rank matrixes, and its element is
k is m × (n-1) rank matrix, and its element is
l is m rank square formations, and its element is
Step 25, solve linear revise system of equations
Wherein
obtain the correction amount θ of each node voltage amplitude and phase angle, Δ V, upgrade node voltage, jump to step 23;
Step 26, be calculated as follows all branch powers:
Wherein, i is the first node of branch road, and j is branch road end-node, and tilde represents the conjugate getting plural number.
Four, the initialization of adaptive differential algorithm (DE).Main step comprises the initialization of each node load and carries out group expanding process to initial population and obtain population X, the setting of differential evolution algorithm maximum iteration time Gm, self adaptive pantographic factor minimum value F0
minwith maximal value F0
maxsetting, adaptive crossover mutation factor minimum value CR
minwith maximal value CR
maxisoparametric.Wherein, the mean value generation initial population of each load upper lower limit value is got in the initialization of node load, and concrete formula is:
in formula
s i be respectively maximal value and the minimum value of the nominal load of node i;
To the population at individual X preset
i={ S
i1, S
i2..., S
imcarry out group expanding operation, obtain initialization population Ω={ X
1, X
2..., X
np; Wherein np is default group expanding parameter, each individual X in population
iin the generation rule of each variable as follows:
In formula,
with
be respectively X
ithe minimum value of a middle jth component and maximal value, rand (0,1) is the uniform random number between (0,1);
Wherein, adaptive differential evolution algorithm committed step comprises:
1, cross processing
The individual vector of Stochastic choice two generates difference vector, by the difference vector of generation and another addition of vectors of Stochastic choice, generates variation vector.Concrete formula is:
x in formula
r1, x
r2, x
r3to represent in population 3 different individualities.T represents current state, and t+1 represents state of future generation.
Wherein, the zoom factor of mutation operation adopts adaptive strategy, and concrete formula is:
f0 in formula
l, F0
ube respectively the bound of F0, f
t1, f
t2, f
t3be respectively
fitness;
2, variation process
To be made a variation vector
intersect with object vector
generate intersection vector
concrete formula is:
In formula, rand (j) ∈ [0,1] is equally distributed random function, and CR is the crossover probability factor.
Wherein, the crossover probability factor adopts adaptive strategy, and concrete formula is:
in formula, CRmin, CRmax are respectively the maximum crossover probability factor of minimum crossover probability Summing Factor, and T is maximum iteration time;
3, process is selected
If
fitness
be better than
fitness
then use
replace
and be selected as the next generation; Otherwise
as the next generation.Adopt the greedy search strategy preset, with maximum load level for objective function carries out selection operation.
Five, realization is optimized.
Step 31, judge whether current iterations reaches maximal value, if do not reach, carry out step 32, if the precision of the iterations iterations or calculating that reach setting reaches requirement, jump procedure 39;
Step 32, the population X after group expanding carried out to variation process, obtain variation population Xvar: the individual vector of Stochastic choice two generates difference vector, by difference vector and another addition of vectors of Stochastic choice of generating, generate variation vectorial; Wherein, the zoom factor of mutation operation adopts adaptive strategy, and concrete formula is:
f0 in formula
min, F0
maxbe respectively the bound of default F0, f
t1, f
t2, f
t3be respectively
fitness;
Step 33, the population Xvar of variation is carried out cross processing, obtain cross-species Xcros: the vector that will make a variation intersects with object vector, generates intersection vectorial; Wherein, the crossover probability factor adopts adaptive strategy, and concrete formula is:
in formula, CRmin, CRmax are respectively the maximum crossover probability factor of minimum crossover probability Summing Factor, and T is maximum iteration time;
Step 34, judge whether the traversal of population completes, if do not complete, carry out step 35, otherwise 1 process will be added to iterations, jump to step 31;
Step 35, initial population X and cross-species Xcros to be calculated simultaneously, be substituting in the inferior method of described newton's pressgang by population X and population Xcros respectively according to dimension (representative of every one dimension be the load value of each node), Load flow calculation when ground state and N-1 fault carried out to it; Wherein, the load value of every each node of one-dimensional representation of described dimension;
Step 36, trend when ground state and N-1 fault to be verified, if all verifications by; carry out step 37, otherwise adds 1 process, jump procedure 34 by the variable number of traversal;
Step 37, with preset net capability mathematical model be adaptive response function, computing formula is:
wherein
the active load of representation node i, the i-th number in the one dimension namely in population X or Xcros; The numerical value that fitness is larger will be saved, and this is exactly likely the net capability of trying to achieve;
Step 38, Population Regeneration X: dimension value fitness being greater than preset value is updated in the same one dimension of population X;
Step 39, exit iteration, export the numerical value of net capability.
Wherein, described net capability mathematical model can be described as: take net capability as objective function, according to the definition of net capability, take into account N-1 safety criterion, and consider the actual conditions of the network operation, comprise the constraints such as main transformer capacity, network topology structure, circuit overload ability, the net capability obtained.Concrete mathematical model is as follows:
State variable: the applied power S comprising each load point
diand power-factor angle (angle of impedance)
the voltage magnitude V of node
i, phase angle difference δ
ij.Consider the ruuning situation of 220kV electrical network reality, herein by the power factor of all load point
all be set to 0.98.
Objective function: with maximum power supply capacity for objective function, namely with the active power sum of each load bus for objective function
Constraint condition comprises:
1) load constraint
2) line transmission power constraint
3) node voltage bound constraint
4) phase angle retrains up and down
|δ
i-δ
j|<|δ
i-δ
j|
max(5)
In formula,
s di,
be respectively lower limit and the higher limit of the applied power of node i;
s ij,
be respectively from node i to the lower limit of the line transmission power of node j branch road and higher limit;
v i,
be respectively lower limit and the higher limit of the voltage of node i.
The present invention assesses the method for 220KV handle net net capability, embedded Newton-Laphson method in adaptive differential evolvement method, the effect of embedded Newton-Laphson method is the trend of Exact Solution handle net, the effect of adaptive differential evolvement method processes the trend that Newton-Laphson method is obtained fast, to reach under N-1 constraint, obtain the net capability of handle net fast and accurately; Can be that least unit calculates net capability with section, the method has the features such as good stability, precision is high, ability of searching optimum is strong.
The above embodiment only have expressed several embodiment of the present invention, and it describes comparatively concrete and detailed, but therefore can not be interpreted as the restriction to the scope of the claims of the present invention.It should be pointed out that for the person of ordinary skill of the art, without departing from the inventive concept of the premise, can also make some distortion and improvement, these all belong to protection scope of the present invention.Therefore, the protection domain of patent of the present invention should be as the criterion with claims.
Claims (6)
1. assess a method for 220KV handle net net capability, it is characterized in that, comprise the steps:
Obtain the grid structure of 220KV handle net, determine number of network node and the circuitry number of described grid structure, and the node of described grid structure and branch road are numbered; Wherein, the balance node of 220kV handle net is the 220kV side of 500kV transforming plant main transformer, and PQ node is the load bus in 220kV handle net;
To the grid structure of the 220KV handle net after numbering, read the voltage magnitude V of grid structure parameter, load parameter and balance node
nand phase angle theta
n;
According to the node voltage V preset
iwith phase angle difference δ
ijconstraint, and described PQ node and balance node, obtain the Newton-Laphson method after initialization;
According to the initialization of each node load, the initial population preset carried out to population X, default differential evolution algorithm maximum iteration time Gm, self adaptive pantographic factor minimum value F0 that group expanding process obtains
minwith maximal value F0
max, adaptive crossover mutation factor minimum value CR
minwith maximal value CR
max, obtain the adaptive differential evolution algorithm after initialization; Wherein, according to following formula, the mean value initialization of node load being got to each load upper lower limit value generates initial population:
in formula
be respectively maximal value and the minimum value of the nominal load of node i;
Solve the trend of described 220KV handle net according to described Newton-Laphson method, by described adaptive differential evolution algorithm, described trend is processed, obtain the net capability that the described 220KV handle under N-1 constraint is netted.
2. assess the method for 220KV handle net net capability according to claim 1, it is characterized in that, the described step solving the trend of described 220KV handle net according to described Newton-Laphson method comprises:
Step 21, input raw data: the constraint of the parameter of the voltage Vn of the 220kV side of two 500KV transforming plant main transformers and phase angle theta n, rack and voltage constraint and phase angle difference;
Step 22, formation bus admittance matrix, and according to N-1 fault, bus admittance matrix is modified;
Step 23, calculate each node power amount of unbalance
Judge whether strength of current deviation meets the condition of convergence; As met, then jumping to step 26, if do not met, then carry out step 24; Wherein, the computing formula of trend deviation is as follows:
Wherein, Δ P, Δ Q are respectively active power deviation and the reactive power deviation of node, P
is, Q
isbe the active power and reactive power that i-th node is given; V
i, V
jbe respectively the voltage of i-th node and a jth node; G
ij, B
ijbe respectively from node i to the conductance of node j branch road and susceptance; θ
ijfor the phase angle difference of node i and node j;
Step 24, generate trend by the variable inputted and existing bus admittance matrix, calculate Jacobi matrix J:
Wherein, H is (n-1) rank square formations, and its element is
n is (n-1) × m rank matrixes, and its element is
k is m × (n-1) rank matrix, and its element is
l is m rank square formations, and its element is
Step 25, solve linear revise system of equations
Wherein
obtain the correction amount θ of each node voltage amplitude and phase angle, Δ V, upgrade node voltage, jump to step 23;
Step 26, be calculated as follows all branch powers:
Wherein, i is the first node of branch road, and j is branch road end-node, and tilde represents the conjugate getting plural number.
3. calculate the method for 220KV handle net net capability according to claim 2, it is characterized in that, describedly processed described trend by described adaptive differential evolution algorithm, the step obtaining the net capability that the lower described 220KV handle of N-1 constraint is netted comprises:
Step 31, judge whether current iterations reaches maximal value, if do not reach, carry out step 32, if the precision of the iterations iterations or calculating that reach setting reaches requirement, jump procedure 39;
Step 32, the population X after group expanding carried out to variation process, obtain variation population Xvar: the individual vector of Stochastic choice two generates difference vector, by difference vector and another addition of vectors of Stochastic choice of generating, generate variation vectorial;
Step 33, the population Xvar of variation is carried out cross processing, obtain cross-species Xcros: the vector that will make a variation intersects with object vector, generates intersection vectorial;
Step 34, judge whether the traversal of population completes, if do not complete, carry out step 35, otherwise 1 process will be added to iterations, jump to step 31;
Step 35, initial population X and cross-species Xcros to be calculated simultaneously, be substituting to respectively in the inferior method of described newton's pressgang according to dimension by population X and population Xcros, carry out Load flow calculation when ground state and N-1 fault; Wherein, the load value of every each node of one-dimensional representation of described dimension;
Step 36, trend when ground state and N-1 fault to be verified, if all verifications by; carry out step 37, otherwise adds 1 process, jump procedure 34 by the variable number of traversal;
Step 37, with preset net capability mathematical model be adaptive response function, computing formula is:
wherein
the active load of representation node i, the i-th number in the one dimension namely in population X or Xcros;
Step 38, Population Regeneration X: the dimension value being greater than preset value from appropriateness is updated in the same one dimension of population X;
The numerical value of step 39, output net capability.
4. calculate the method for 220KV handle net net capability according to claim 3, it is characterized in that, the state variable of described net capability mathematical model comprises: the applied power S of each load point
di, power-factor angle
the voltage magnitude V of node
iand phase angle difference δ
ij;
The objective function of described net capability mathematical model is maximum power supply capacity, calculates the active power sum of each load bus as described objective function by following formula:
The constraint condition of described net capability mathematical model is: load retrains:
line transmission power constraint:
node voltage bound retrains:
and phase angle retrains up and down: | δ
i-δ
j| < | δ
i-δ
j|
max;
Wherein,
be respectively lower limit and the higher limit of the applied power of node i;
be respectively from node i to the lower limit of the line transmission power of node j branch road and higher limit;
be respectively lower limit and the higher limit of the voltage of node i.
5. assess the method for 220KV handle net net capability according to claim 3, it is characterized in that, described cross processing is: generate difference vector according to the individual vector of following formula Stochastic choice two, by the difference vector of generation and another addition of vectors of Stochastic choice, generates variation vector:
x in formula
r1, x
r2, x
r3to represent in population 3 different individualities.
6. assess the method for 220KV handle net net capability according to claim 3, it is characterized in that, described variation process comprises:
To make a variation vector according to following formula
intersect with object vector
generate intersection vector
The described crossover probability factor adopts adaptive strategy:
in formula, CRmin, CRmax are respectively the maximum crossover probability factor of minimum crossover probability Summing Factor, and T is maximum iteration time.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410706099.7A CN104484832B (en) | 2014-11-26 | 2014-11-26 | The method of assessment 220KV handle net net capability |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410706099.7A CN104484832B (en) | 2014-11-26 | 2014-11-26 | The method of assessment 220KV handle net net capability |
Publications (2)
Publication Number | Publication Date |
---|---|
CN104484832A true CN104484832A (en) | 2015-04-01 |
CN104484832B CN104484832B (en) | 2015-08-19 |
Family
ID=52759372
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201410706099.7A Active CN104484832B (en) | 2014-11-26 | 2014-11-26 | The method of assessment 220KV handle net net capability |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN104484832B (en) |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105490288A (en) * | 2016-01-06 | 2016-04-13 | 华南理工大学 | Reactive compensation optimization configuration method for 220kV power network |
CN105552906A (en) * | 2016-02-14 | 2016-05-04 | 华南理工大学 | Regional power grid load margin analysis method based on prime-dual interior point method |
CN105552907A (en) * | 2016-02-14 | 2016-05-04 | 华南理工大学 | Method for evaluating maximum power supply capacity of 220kV double-loop hand-holding grid |
CN105576653A (en) * | 2016-01-06 | 2016-05-11 | 华南理工大学 | 220kV district power grid power supply capacity optimization method |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20070027642A1 (en) * | 2005-07-15 | 2007-02-01 | Chang Gung University | Method for Calculating Power Flow Solution of a Power Transmission Network that Includes Interline Power Flow Controller (IPFC) |
CN103199544A (en) * | 2013-03-26 | 2013-07-10 | 上海理工大学 | Reactive power optimization method of electrical power system |
CN103701117A (en) * | 2013-12-19 | 2014-04-02 | 浙江大学 | Active distribution network dynamic topology reconstitution method based on mixed artificial intelligence |
CN104037764A (en) * | 2014-07-03 | 2014-09-10 | 大连海事大学 | Rectangular coordinate Newton method load flow calculation method with changeable Jacobian matrix |
-
2014
- 2014-11-26 CN CN201410706099.7A patent/CN104484832B/en active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20070027642A1 (en) * | 2005-07-15 | 2007-02-01 | Chang Gung University | Method for Calculating Power Flow Solution of a Power Transmission Network that Includes Interline Power Flow Controller (IPFC) |
CN103199544A (en) * | 2013-03-26 | 2013-07-10 | 上海理工大学 | Reactive power optimization method of electrical power system |
CN103701117A (en) * | 2013-12-19 | 2014-04-02 | 浙江大学 | Active distribution network dynamic topology reconstitution method based on mixed artificial intelligence |
CN104037764A (en) * | 2014-07-03 | 2014-09-10 | 大连海事大学 | Rectangular coordinate Newton method load flow calculation method with changeable Jacobian matrix |
Non-Patent Citations (2)
Title |
---|
杨旭英等: "含有分布式发电的配电网潮流计算", 《中国高等学校电力***及其自动化专业第二十四届学术年会论文集》 * |
樊宇璐,等: "基于拟牛顿法的电力***潮流计算", 《中国高等学校电力***及其自动化专业第二十四届学术年会论文集》 * |
Cited By (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105490288A (en) * | 2016-01-06 | 2016-04-13 | 华南理工大学 | Reactive compensation optimization configuration method for 220kV power network |
CN105576653A (en) * | 2016-01-06 | 2016-05-11 | 华南理工大学 | 220kV district power grid power supply capacity optimization method |
CN105490288B (en) * | 2016-01-06 | 2018-02-27 | 华南理工大学 | A kind of 220kV electric network reactive compensations Optimal Configuration Method |
CN105576653B (en) * | 2016-01-06 | 2018-02-27 | 华南理工大学 | A kind of 220kV sections power network power supply capacity optimization method |
CN105552906A (en) * | 2016-02-14 | 2016-05-04 | 华南理工大学 | Regional power grid load margin analysis method based on prime-dual interior point method |
CN105552907A (en) * | 2016-02-14 | 2016-05-04 | 华南理工大学 | Method for evaluating maximum power supply capacity of 220kV double-loop hand-holding grid |
CN105552906B (en) * | 2016-02-14 | 2018-04-13 | 华南理工大学 | A kind of area power grid load nargin analysis method based on prim al- dual interior point m ethod |
Also Published As
Publication number | Publication date |
---|---|
CN104484832B (en) | 2015-08-19 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN104463701B (en) | A kind of distribution system and the coordinated planning method of charging electric vehicle network | |
Ghadimi et al. | A PSO-based fuzzy long-term multi-objective optimization approach for placement and parameter setting of UPFC | |
CN106487005A (en) | A kind of Electric power network planning method considering T-D tariff | |
CN105321003B (en) | A kind of ac and dc systems multiple target tide optimization method containing VSC-HVDC | |
CN107979092A (en) | It is a kind of to consider distributed generation resource and the power distribution network dynamic reconfiguration method of Sofe Switch access | |
CN106253335A (en) | A kind of distributed power source capacity and on-position uncertain distribution network planning method | |
CN105225163A (en) | The reconstructing method of active distribution network and device | |
CN107947192A (en) | A kind of optimal reactive power allocation method of droop control type isolated island micro-capacitance sensor | |
CN103077480B (en) | Safety checking method for power system | |
CN106374513B (en) | A kind of more microgrid dominant eigenvalues optimization methods based on leader-followers games | |
CN106972504A (en) | Interval idle work optimization method based on genetic algorithm | |
CN103746374A (en) | Closed loop control method comprising multi-microgrid power distribution network | |
CN106786493A (en) | A kind of practical calculation method of multi-infeed HVDC interaction factor | |
CN106532710A (en) | Microgrid power flow optimization method considering voltage stability constraint | |
CN103208797A (en) | Estimation method for new-energy-containing power distribution network state based on intelligent optimization technology | |
CN105529703B (en) | A kind of urban network reconstruction planing method based on power supply capacity bottleneck analysis | |
CN105243432A (en) | NSGA-III-based UPFC (unified power flow controller) location and capacity multi-target configuration method | |
CN104484832A (en) | Method for evaluating total supplying capability of 220KV Lashou net | |
CN104484555B (en) | The method of assessment 220kV self-healing looped network net capability | |
CN104767207B (en) | Voltage stabilization prevention and control method based on generator reactive deposit sensitivity | |
CN105490288B (en) | A kind of 220kV electric network reactive compensations Optimal Configuration Method | |
CN105633955B (en) | A kind of multifunctional solid current limiter Optimal Configuration Method | |
CN104410065B (en) | A kind of Multiobjective Decision Making Method of receiving end electrical network limiting short-circuit current | |
CN103887823A (en) | Micro-grid connection position selection method based on fuzzy hierarchical analysis | |
CN105449672A (en) | Method for estimating total supply capability of 220KV loop ring network |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant |