CN105356481A - Multi-infeed-short-circuit-ratio-based dynamic reactive compensation point selection method - Google Patents
Multi-infeed-short-circuit-ratio-based dynamic reactive compensation point selection method Download PDFInfo
- Publication number
- CN105356481A CN105356481A CN201510794734.6A CN201510794734A CN105356481A CN 105356481 A CN105356481 A CN 105356481A CN 201510794734 A CN201510794734 A CN 201510794734A CN 105356481 A CN105356481 A CN 105356481A
- Authority
- CN
- China
- Prior art keywords
- centerdot
- represent
- theta
- prime
- bus
- 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
Classifications
-
- Y—GENERAL 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
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02E—REDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
- Y02E40/00—Technologies for an efficient electrical power generation, transmission or distribution
- Y02E40/30—Reactive power compensation
Landscapes
- Supply And Distribution Of Alternating Current (AREA)
Abstract
The invention provides a multi-infeed-short-circuit-ratio-based dynamic reactive compensation point selection method. The method comprises the following steps: establishing a jacobian matrix and solving a voltage impact factor; determining a multi-infeed short circuit ratio according to the voltage impact factor; establishing effect evaluation objective functions of dynamic reactive compensation apparatus installation at all stations and corresponding constraint conditions; and selecting an optimal dynamic reactive compensation apparatus installation point by using a genetic algorithm. According to the invention, with consideration of the interactive influence between multi-loop direct-current lines, an effective method and an effective technical support are provided for selection of a dynamic reactive compensation apparatus installation point at a direct-current multi-infeed area of the country; and because the influences on the system stability by response characteristics of all elements during the dynamic process are taken into consideration, the point selection method is practical and has the engineering application value. Dynamic point selection optimization is carried out by using the genetic algorithm, so that the calculation efficiency is high.
Description
Technical field
The invention belongs to technical field of power systems, be specifically related to a kind of dynamic passive compensation reconnaissance method based on many feed-ins short circuit ratio.
Background technology
The load center areas such as China's East China Yangtze River Delta are by the restriction of local normal power supplies development, and will present many direct currents future concentrates feed-in structure, and the system failure may cause multiple-circuit line chain-react, cause Voltage-stabilizing Problems.Domestic and international operation of power networks experience shows, every REACTIVE POWER/VOLTAGE CONTROL problem of relating to, dynamic reactive compensation device is more satisfactory solution, dynamic reactive compensation device control strategy is flexible, fast response time, the Voltage-stabilizing Problems of multi-infeed DC receiving end electrical network can be considered to adopt dynamic reactive compensation device, improves the voltage support ability of extra-high voltage grid, particularly after system malfunctions disturbance, help the fast quick-recovery realizing system weak spot voltage.
Consider that receiving end Net Frame of Electric Network is intensive, land resource scarcity, the restriction that dynamic reactive compensation device is easily subject to site is installed in transformer station, so, its dynamic passive compensation allocation plan can when reaching effect same, the scheme that capacity is minimum or configuration site is minimum should being selected, needing to determine by optimizing.
Many feed-ins short circuit ratio is more common for the index of DC support ability for measurement interchange in multi-infeed DC electrical network, be widely used, but it is defined as static index, only can reflects network topology structure, cannot take into account dynamic element model.And stability of power system focuses on the system performance in dynamic process under each element responds, various dynamic element is as larger on the impact of voltage stability in generator model, load model, DC control model, dynamic passive compensation model etc.At present, domestic and international existing dynamic reactive compensation device collocation method generally carries out reconnaissance based on the weakness zone of static voltage stability, not yet considers multi-infeed DC system performance and dynamic element model characteristics.
Summary of the invention
In order to overcome above-mentioned the deficiencies in the prior art, the invention provides a kind of dynamic passive compensation reconnaissance method based on many feed-ins short circuit ratio, reconnaissance method more gears to actual circumstances, and possesses engineer applied and is worth.
In order to realize foregoing invention object, the present invention takes following technical scheme:
The invention provides a kind of dynamic passive compensation reconnaissance method based on many feed-ins short circuit ratio, said method comprising the steps of:
Step 1: set up Jacobian matrix and solve the voltage influence factor;
Step 2: determine many feed-ins short circuit ratio according to the voltage influence factor;
Step 3: set up each website and install dynamic reactive compensation device effect assessment target function and corresponding constraints;
Step 4: adopt the preferred dynamic reactive compensation device mounting points of genetic algorithm.
Described step 1 comprises the following steps:
Step 1-1: consider dynamic element model, set up power balance equation;
Step 1-2: set up Jacobian matrix, and solve the voltage influence factor.
In described step 1-1, dynamic element model comprises generator model, load model, DC control model and dynamic reactive compensation device model.
In described step 1-1, set up following power balance equation:
Wherein, Δ P
i, Δ Q
irepresent the active power variable quantity that node i is injected and reactive power variable quantity respectively, P
gi, Q
girepresent that generator injects node i meritorious and exerts oneself and idlely to exert oneself, P respectively
li, Q
lirepresent burden with power and the load or burden without work of node i respectively, P
direpresent the direct current power of node i, Q
direpresent that DC filter capacitor injects the reactive power of node i, U
i, U
jrepresent the voltage of node i, j respectively, Q
sirepresent that dynamic reactive compensation device injects the idle of node i and exerts oneself, G
ij, B
ijrepresent the conductance between node i, j and susceptance respectively, θ
ijrepresent the phase difference of voltage between node i, j, i=1,2 ..., n, j=1,2 ..., n, n are node total number.
In described step 1-2, set up following Jacobian matrix equation:
Wherein, MIIF
lkrepresent the voltage influence factor of bus l opposing busbars k, and MIIF
lk=Δ U
l/ Δ U
k, Δ U
lrepresent the voltage variety of bus l, Δ U
krepresent the voltage variety of bus k;
Jacobian matrix element H
ii, N
ii, M
ii, L
iicalculate according to the following formula:
Wherein, G
iirepresent the conductance of node i, B
iirepresent the susceptance of node i;
be expressed as:
Wherein, E
i" represent generator electromotive force, θ
δ irepresent E
i" and U
iphase angle difference, X "
direpresent that generator d axle surpasses transient state reactance; If node i is constant current load bus,
be expressed as:
Wherein, I
pi, I
qirepresent active current and the reactive current of constant current load bus respectively;
If node i is constant impedance load bus,
be expressed as:
Wherein, G
i, B
irepresent conductance and the susceptance of constant impedance load bus respectively;
If node i is DC line access node,
be expressed as:
Wherein, I
drepresent direct current, n
trepresent six pulse conversion device series connection numbers, k
trepresent converter transformer no-load voltage ratio, k
γrepresent converter transformer equivalence no-load voltage ratio, θ
drepresent the direct current angle of overlap of rectification side or the extinguish angle of inverter side, X
crepresent equivalent commutating reactance;
represent Equivalent Power Factor angle, and
If node i is dynamic reactive compensation device install node,
be expressed as:
Wherein, B
irepresent that dynamic passive compensation installs the susceptance of node, and B
i=-K Δ U
i=-K (U
i-U
i0), Δ U
irepresent the voltage deviation installed before and after dynamic reactive compensation device, U
i0represent that dynamic reactive compensation device installs the initial voltage of node, K represents proportionality coefficient.
In described step 2, the drop point that p returns DC line is bus l, and the drop point that q returns DC line is bus k, determines many feed-ins short circuit ratio, have according to the voltage influence factor:
Wherein, MISCR '
prepresent that p returns many feed-ins short circuit ratio of DC line, S
lrepresent the system short circuit capacity of bus l, P
prepresent that p returns the power of DC line, P
qrepresent that q returns the power of DC line, MIIF
lkrepresent the voltage influence factor of bus l opposing busbars k, and MIIF
lk=Δ U
l/ Δ U
k, Δ U
lrepresent the voltage variety of bus l, Δ U
krepresent the voltage variety of bus k, m represents that in direct current many feed-ins electrical network, returning of DC line counts.
In described step 3, set up each website and dynamic reactive compensation device effect assessment target function be installed, have:
Wherein, f represents that each website installs dynamic reactive compensation device effect assessment target function,
represent that p returns the equivalent weight coefficient of DC line, have:
Wherein, ω
prepresent that p returns the weight coefficient of DC line, ω
qrepresent that q returns the weight coefficient of DC line, ω
p, ω
qreflect that p, q return the impact of DC line on other direct currents in direct current many feed-ins electrical network respectively, have:
Wherein, Z
pqrepresent that p returns DC line change of current bus and q and returns equiva lent impedance between DC line change of current bus, Z
pprepresent that p returns the equiva lent impedance of DC line change of current bus, Z
qqrepresent that q returns the equiva lent impedance of DC line change of current bus;
The constraints that each website installs dynamic reactive compensation device effect assessment target function additional corresponding is as follows:
Wherein, P
l, Q
lrepresent active power and the reactive power of bus l respectively, U
l, U
krepresent the voltage of bus l, k respectively,
with
represent upper voltage limit and the lower limit of bus k respectively, G
lk, B
lkrepresent the conductance between bus l, k and susceptance respectively, θ
lkrepresent the phase difference of voltage between bus l, k, I
krepresent the short circuit current of bus k,
represent the short circuit current upper limit of bus k, S
rrepresent the power of branch road r,
represent the power upper limit of branch road r.
In described step 4, adopt the preferred dynamic reactive compensation device mounting points of genetic algorithm, comprising:
(1) arranging evolutionary generation t is 0, and arranges maximum evolutionary generation T and individual in population number M;
(2) in n node, select s node as dynamic reactive compensation device mounting points, then have
the individual dynamic reactive compensation device reconnaissance strategy that may exist, and gene code is carried out to reactive power compensator reconnaissance strategy, gene code length L meets
the value of each bit of gene code length is 0 or 1;
(3) stochastic generation M individuality is as initial population P
0, each dynamic reactive compensation device reconnaissance strategy is as body one by one;
(4) calculate the fitness of each individuality, namely each website installs dynamic reactive compensation device effect assessment target function value;
(5) for each parent individuality distributes a random number, and according to corresponding random number, parent individuality is sorted according to order from big to small, adjacent two parent individualities are hybridized, certain bit in Stochastic choice gene order, 0 of bit or 1 exchange by two parent individualities, produce offspring individual, calculate the fitness of offspring individual, and offspring individual is joined in parent individuality composition parent colony;
(6) parent colony is screened, M individuality before retaining according to ideal adaptation degree size;
(7) certain individuality of Stochastic choice, and a bit in this genes of individuals sequence of Stochastic choice, overturn, obtain colony of future generation;
(8) if t=T, then the individuality in evolutionary process with maximum adaptation degree exports as optimal solution, namely completes the preferred of dynamic reactive compensation device mounting points; If t<T, then repeat (4) ~ (8).
Compared with immediate prior art, technical scheme provided by the invention has following beneficial effect:
1. the present invention utilizes many feed-ins short circuit ratio as the measurement index of idle input, carry out preferably for dynamic reactive compensation device mounting points in direct current many feed-ins electrical network, consider the reciprocal effect between multiple-circuit line circuit, the choosing of dynamic reactive compensation device mounting points for China's direct current many feed-ins area provides effective ways and technical support;
2. the present invention considers multi-infeed DC system performance and dynamic element model characteristics, the dynamic element models such as generator, load, direct current, dynamic reactive compensation device are brought in the measurement index of many feed-ins short circuit ratio, to consider in dynamic process each element responds characteristic to the impact of the stability of a system, reconnaissance method more gears to actual circumstances, and possesses engineer applied and is worth;
3. the present invention utilizes genetic algorithm to carry out dynamic reconnaissance optimization, and not based on time-domain-simulation, computational efficiency is high.
Accompanying drawing explanation
Fig. 1 is the dynamic passive compensation reconnaissance method flow diagram based on many feed-ins short circuit ratio in the embodiment of the present invention;
Fig. 2 is the geographical wiring schematic diagram of the year two thousand twenty East China direct current many feed-ins electrical network in the embodiment of the present invention 2.
Embodiment
Below in conjunction with accompanying drawing, the present invention is described in further detail.
The invention provides a kind of dynamic passive compensation reconnaissance method based on many feed-ins short circuit ratio, as Fig. 1, said method comprising the steps of:
Step 1: set up Jacobian matrix and solve the voltage influence factor;
Step 2: determine many feed-ins short circuit ratio according to the voltage influence factor;
Step 3: set up each website and install dynamic reactive compensation device effect assessment target function and corresponding constraints;
Step 4: adopt the preferred dynamic reactive compensation device mounting points of genetic algorithm.
In described step 1, dynamic element model comprises generator model, load model, DC control model and dynamic reactive compensation device model.
Described step 1 comprises the following steps:
Step 1-1: consider dynamic element model, set up power balance equation;
Step 1-2: set up Jacobian matrix, and solve the voltage influence factor.
In described step 1-1, set up following power balance equation:
Wherein, Δ P
i, Δ Q
irepresent the active power variable quantity that node i is injected and reactive power variable quantity respectively, P
gi, Q
girepresent that generator injects node i meritorious and exerts oneself and idlely to exert oneself, P respectively
li, Q
lirepresent burden with power and the load or burden without work of node i respectively, P
direpresent the direct current power of node i, Q
direpresent that DC filter capacitor injects the reactive power of node i, U
i, U
jrepresent the voltage of node i, j respectively, Q
sirepresent that dynamic reactive compensation device injects the idle of node i and exerts oneself, G
ij, B
ijrepresent the conductance between node i, j and susceptance respectively, θ
ijrepresent the phase difference of voltage between node i, j, i=1,2 ..., n, j=1,2 ..., n, n are node total number.
In described step 1-2, set up following Jacobian matrix equation:
Wherein, MIIF
lkrepresent the voltage influence factor of bus l opposing busbars k, and MIIF
lk=Δ U
l/ Δ U
k, Δ U
lrepresent the voltage variety of bus l, Δ U
krepresent the voltage variety of bus k;
Jacobian matrix element H
ii, N
ii, M
ii, L
iicalculate according to the following formula:
Wherein, G
iirepresent the conductance of node i, B
iirepresent the susceptance of node i;
be expressed as:
Wherein, E
i" represent generator electromotive force, θ
δ irepresent E
i" and U
iphase angle difference, X "
direpresent that generator d axle surpasses transient state reactance;
If node i is constant current load bus,
be expressed as:
Wherein, I
pi, I
qirepresent active current and the reactive current of constant current load bus respectively;
If node i is constant impedance load bus,
be expressed as:
Wherein, G
i, B
irepresent conductance and the susceptance of constant impedance load bus respectively;
If node i is DC line power node,
be expressed as:
Wherein, I
drepresent direct current, n
trepresent six pulse conversion device series connection numbers, k
trepresent converter transformer no-load voltage ratio, k
γrepresent converter transformer equivalence no-load voltage ratio, θ
drepresent the direct current angle of overlap of rectification side or the extinguish angle of inverter side, X
crepresent equivalent commutating reactance;
represent Equivalent Power Factor angle, and
If node i is dynamic reactive compensation device install node,
be expressed as:
Wherein, B
irepresent that dynamic passive compensation installs the susceptance of node, and B
i=-K Δ U
i=-K (U
i-U
i0), Δ U
irepresent the voltage deviation installed before and after dynamic reactive compensation device, U
i0represent that dynamic reactive compensation device installs the initial voltage of node, K represents proportionality coefficient.
In described step 2, the drop point that p returns DC line is bus l, and the drop point that q returns DC line is bus k, determines many feed-ins short circuit ratio, have according to the voltage influence factor:
Wherein, MISCR '
prepresent that p returns many feed-ins short circuit ratio of DC line, S
lrepresent the system short circuit capacity of bus l, P
prepresent that p returns the power of DC line, P
qrepresent that q returns the power of DC line, MIIF
lkrepresent the voltage influence factor of bus l opposing busbars k, and MIIF
lk=Δ U
l/ Δ U
k, Δ U
lrepresent the voltage variety of bus l, Δ U
krepresent the voltage variety of bus k, m represents that in direct current many feed-ins electrical network, returning of DC line counts.
In described step 3, set up each website and dynamic reactive compensation device effect assessment target function be installed, have:
Wherein, f represents that each website installs dynamic reactive compensation device effect assessment target function,
represent that p returns the equivalent weight coefficient of DC line, have:
Wherein, ω
prepresent that p returns the weight coefficient of DC line, ω
qrepresent that q returns the weight coefficient of DC line, ω
p, ω
qreflect that p, q return the impact of DC line on other direct currents in direct current many feed-ins electrical network respectively, have:
Wherein, Z
pqrepresent that p returns DC line change of current bus and q and returns equiva lent impedance between DC line change of current bus, Z
pprepresent that p returns the equiva lent impedance of DC line change of current bus, Z
qqrepresent that q returns the equiva lent impedance of DC line change of current bus;
The constraints that each website installs dynamic reactive compensation device effect assessment target function additional corresponding is as follows:
Wherein, P
l, Q
lrepresent active power and the reactive power of bus l respectively, U
l, U
krepresent the voltage of bus l, k respectively,
with
represent upper voltage limit and the lower limit of bus k respectively, G
lk, B
lkrepresent the conductance between bus l, k and susceptance respectively, θ
lkrepresent the phase difference of voltage between bus l, k, I
krepresent the short circuit current of bus k,
represent the short circuit current upper limit of bus k, S
rrepresent the power of branch road r,
represent the power upper limit of branch road r.
In described step 4, adopt the preferred dynamic reactive compensation device mounting points of genetic algorithm, comprising:
(1) arranging evolutionary generation t is 0, and arranges maximum evolutionary generation T and individual in population number M;
(2) in n node, select s node as dynamic reactive compensation device mounting points, then have
the individual dynamic reactive compensation device reconnaissance strategy that may exist, and gene code is carried out to reactive power compensator reconnaissance strategy, gene code length L meets
the value of each bit of gene code length is 0 or 1;
(3) stochastic generation M individuality is as initial population P
0, each dynamic reactive compensation device reconnaissance strategy is as body one by one;
(4) calculate the fitness of each individuality, namely each website installs dynamic reactive compensation device effect assessment target function value;
(5) for each parent individuality distributes a random number, and according to corresponding random number, parent individuality is sorted according to order from big to small, adjacent two parent individualities are hybridized, certain bit in Stochastic choice gene order, 0 of bit or 1 exchange by two parent individualities, produce offspring individual, calculate the fitness of offspring individual, and offspring individual is joined in parent individuality composition parent colony;
(6) parent colony is screened, M individuality before retaining according to ideal adaptation degree size;
(7) certain individuality of Stochastic choice, and a bit in this genes of individuals sequence of Stochastic choice, overturn, obtain colony of future generation;
(8) if t=T, then the individuality in evolutionary process with maximum adaptation degree exports as optimal solution, namely completes the preferred of dynamic reactive compensation device mounting points; If t<T, then repeat (4) ~ (8).
Embodiment 1
Reconnaissance method provided by the invention is applied to present East China Power Grid, carries out many feed-ins short circuit ratio and calculate.If be MISCR with traditional reconnaissance method acquired results
p, reconnaissance method acquired results provided by the invention is MISCR
p'.East China Power Grid feed-in direct current 9 times, shown in East China Power Grid direct current many feed-ins short circuit ratio result of calculations 1 in 2016:
Table 1
In table 1, MISCR' does not consider hvdc control mode, and load model is " 40% constant impedance+60% firm power ".
Result as can be seen from table, two class methods result of calculations have different, but substantially do not change according to the sequence of many feed-ins short circuit ratio size.
To same 9 times direct currents, the many feed-ins short circuit ratio calculated under different DC control strategy calculates, and under the different control strategies of East China Power Grid direct current in 2016, many feed-ins short circuit ratio result of calculation is as table 2:
Table 2
Table 2 result of calculation shows: under most cases, bigger when the many feed-ins short circuit ratio ratio under Given current controller mode ignores DC control; Constant dc power control mode many feed-ins short circuit ratio is little.This be due to, the reason that constant dc power control mode makes reactive voltage support situation variation is: when on inverter side change of current bus, voltage drop is low, converter makes extinguish angle remain unchanged by regulating Trigger Angle, and inverter side direct voltage reduces, thus causes direct current to raise; Due to the impact of commutating reactance, inverter side direct voltage reduces further, causes inverter side converter power factor to reduce, and reactive requirement raises.For rectification side, because direct current raises, in order to keep power invariability, must direct voltage be reduced, thus also make reactive requirement raise.
Illustrated by above-mentioned result of calculation, the many feed-ins short circuit ratio result of calculation adopting the inventive method to derive and former short circuit ratio result of calculation are consistent sexual intercourse, reflect, the model of all kinds of dynamic element is different, will produce certain influence to the result of calculation of index and stability thereof simultaneously.
Embodiment 2
Reconnaissance method provided by the invention is applied to the year two thousand twenty Jiangsu planning electrical network, as shown in Figure 2, total political affairs are put down, with inner, Liyang, Taizhou, Nanjing, Changshu, 7, Changzhou direct current drop point, are formed typical multi-infeed DC electrical network, the powered ratio 40% of Jiangsu Power Grid.Because feed-in direct current scale is excessive, when not installing dynamic reactive compensation device, Jiangsu has the 14 Flow Line generation three-phase permanent short faults that backcross to cause Voltage Instability.
Choose 25, southern area of Jiangsu Province 500kV transformer station as alternative point, as shown in table 3, select 5 websites that the SVC of 2 × 240Mvar capacity is installed by optimizer;
Table 3
Set different evolutionary generation T, colony individual amount M, and be optimized calculating, southern Jiangsu dynamic passive compensation reconnaissance scenario outcomes is as shown in table 4, from result of calculation, adopt different evolutionary generations and individual amount on result of calculation impact not quite, result of calculation substantially can be stabilized in one and compare the scheme determined.Work as T=30, during M=30, its result and all the other schemes slightly difference, but from target function value, both differences are very little.
Table 4
From result in table, adopt genetic algorithm to solve the problem, good convergence can be reached.Carry out analysis safety and stability to above-mentioned two different reactive power compensation allocation plans to check, result shows, and under two class schemes, N-1 causes the circuit of system unstability to be 7 times, to the lifting successful of the stability of a system.
As can be seen here, the dynamic reactive compensation device optimization reconnaissance method computational efficiency in the direct current many feed-ins area carried herein is high, and effect of optimization is good, has very strong engineer applied and is worth.
Finally should be noted that: above embodiment is only in order to illustrate that technical scheme of the present invention is not intended to limit; those of ordinary skill in the field still can modify to the specific embodiment of the present invention with reference to above-described embodiment or equivalent replacement; these do not depart from any amendment of spirit and scope of the invention or equivalent replacement, are all applying within the claims of the present invention awaited the reply.
Claims (8)
1., based on a dynamic passive compensation reconnaissance method for many feed-ins short circuit ratio, it is characterized in that: said method comprising the steps of:
Step 1: set up Jacobian matrix and solve the voltage influence factor;
Step 2: determine many feed-ins short circuit ratio according to the voltage influence factor;
Step 3: set up each website and install dynamic reactive compensation device effect assessment target function and corresponding constraints;
Step 4: adopt the preferred dynamic reactive compensation device mounting points of genetic algorithm.
2. the dynamic passive compensation reconnaissance method based on many feed-ins short circuit ratio according to claim 1, is characterized in that: described step 1 comprises the following steps:
Step 1-1: consider dynamic element model, set up power balance equation;
Step 1-2: set up Jacobian matrix, and solve the voltage influence factor.
3. the dynamic passive compensation reconnaissance method based on many feed-ins short circuit ratio according to claim 2, it is characterized in that: in described step 1-1, dynamic element model comprises generator model, load model, DC control model and dynamic reactive compensation device model.
4. the dynamic passive compensation reconnaissance method based on many feed-ins short circuit ratio according to claim 2, is characterized in that: in described step 1-1, sets up following power balance equation:
Wherein, Δ P
i, Δ Q
irepresent the active power variable quantity that node i is injected and reactive power variable quantity respectively, P
gi, Q
girepresent that generator injects node i meritorious and exerts oneself and idlely to exert oneself, P respectively
li, Q
lirepresent burden with power and the load or burden without work of node i respectively, P
direpresent the direct current power of node i, Q
direpresent that DC filter capacitor injects the reactive power of node i, U
i, U
jrepresent the voltage of node i, j respectively, Q
sirepresent that dynamic reactive compensation device injects the idle of node i and exerts oneself, G
ij, B
ijrepresent the conductance between node i, j and susceptance respectively, θ
ijrepresent the phase difference of voltage between node i, j, i=1,2 ..., n, j=1,2 ..., n, n are node total number.
5. the dynamic passive compensation reconnaissance method based on many feed-ins short circuit ratio according to claim 4, is characterized in that: in described step 1-2, sets up following Jacobian matrix equation:
Wherein, MIIF
lkrepresent the voltage influence factor of bus l opposing busbars k, and MIIF
lk=Δ U
l/ Δ U
k, Δ U
lrepresent the voltage variety of bus l, Δ U
krepresent the voltage variety of bus k;
Jacobian matrix element H
ii, N
ii, M
ii, L
iicalculate according to the following formula:
Wherein, G
iirepresent the conductance of node i, B
iirepresent the susceptance of node i;
be expressed as:
Wherein, E "
irepresent generator electromotive force, θ
δ irepresent E "
iwith U
iphase angle difference, X "
direpresent that generator d axle surpasses transient state reactance;
If node i is constant current load bus,
be expressed as:
Wherein, I
pi, I
qirepresent active current and the reactive current of constant current load bus respectively;
If node i is constant impedance load bus,
be expressed as:
Wherein, G
i, B
irepresent conductance and the susceptance of constant impedance load bus respectively;
If node i is DC line access node,
be expressed as:
Wherein, I
drepresent direct current, n
trepresent six pulse conversion device series connection numbers, k
trepresent converter transformer no-load voltage ratio, k
γrepresent converter transformer equivalence no-load voltage ratio, θ
drepresent the direct current angle of overlap of rectification side or the extinguish angle of inverter side, X
crepresent equivalent commutating reactance;
represent Equivalent Power Factor angle, and
If node i is dynamic reactive compensation device install node,
be expressed as:
Wherein, B
irepresent that dynamic passive compensation installs the susceptance of node, and B
i=-K Δ U
i=-K (U
i-U
i0), Δ U
irepresent the voltage deviation installed before and after dynamic reactive compensation device, U
i0represent that dynamic reactive compensation device installs the initial voltage of node, K represents proportionality coefficient.
6. the dynamic passive compensation reconnaissance method based on many feed-ins short circuit ratio according to claim 5, it is characterized in that: in described step 2, the drop point that p returns DC line is bus l, and the drop point that q returns DC line is bus k, determine many feed-ins short circuit ratio according to the voltage influence factor, have:
Wherein, MISCR '
prepresent that p returns many feed-ins short circuit ratio of DC line, S
lrepresent the system short circuit capacity of bus l, P
prepresent that p returns the power of DC line, P
qrepresent that q returns the power of DC line, MIIF
lkrepresent the voltage influence factor of bus l opposing busbars k, and MIIF
lk=Δ U
l/ Δ U
k, Δ U
lrepresent the voltage variety of bus l, Δ U
krepresent the voltage variety of bus k, m represents that in direct current many feed-ins electrical network, returning of DC line counts.
7. the dynamic passive compensation reconnaissance method based on many feed-ins short circuit ratio according to claim 6, is characterized in that: in described step 3, sets up each website and installs dynamic reactive compensation device effect assessment target function, have:
Wherein, f represents that each website installs dynamic reactive compensation device effect assessment target function,
represent that p returns the equivalent weight coefficient of DC line, have:
Wherein, ω
prepresent that p returns the weight coefficient of DC line, ω
qrepresent that q returns the weight coefficient of DC line, ω
p, ω
qreflect that p, q return the impact of DC line on other direct currents in direct current many feed-ins electrical network respectively, have:
Wherein, Z
pqrepresent that p returns DC line change of current bus and q and returns equiva lent impedance between DC line change of current bus, Z
pprepresent that p returns the equiva lent impedance of DC line change of current bus, Z
qqrepresent that q returns the equiva lent impedance of DC line change of current bus;
The constraints that each website installs dynamic reactive compensation device effect assessment target function additional corresponding is as follows:
Wherein, P
l, Q
lrepresent active power and the reactive power of bus l respectively, U
l, U
krepresent the voltage of bus l, k respectively,
with
represent upper voltage limit and the lower limit of bus k respectively, G
lk, B
lkrepresent the conductance between bus l, k and susceptance respectively, θ
lkrepresent the phase difference of voltage between bus l, k, I
krepresent the short circuit current of bus k,
represent the short circuit current upper limit of bus k, S
rrepresent the power of branch road r,
represent the power upper limit of branch road r.
8. the dynamic passive compensation reconnaissance method based on many feed-ins short circuit ratio according to claim 1, is characterized in that: in described step 4, adopts the preferred dynamic reactive compensation device mounting points of genetic algorithm, comprising:
(1) arranging evolutionary generation t is 0, and arranges maximum evolutionary generation T and individual in population number M;
(2) in n node, select s node as dynamic reactive compensation device mounting points, then have
the individual dynamic reactive compensation device reconnaissance strategy that may exist, and gene code is carried out to reactive power compensator reconnaissance strategy, gene code length L meets
the value of each bit of gene code length is 0 or 1;
(3) stochastic generation M individuality is as initial population P
0, each dynamic reactive compensation device reconnaissance strategy is as body one by one;
(4) calculate the fitness of each individuality, namely each website installs dynamic reactive compensation device effect assessment target function value;
(5) for each parent individuality distributes a random number, and according to corresponding random number, parent individuality is sorted according to order from big to small, adjacent two parent individualities are hybridized, certain bit in Stochastic choice gene order, 0 of bit or 1 exchange by two parent individualities, produce offspring individual, calculate the fitness of offspring individual, and offspring individual is joined in parent individuality composition parent colony;
(6) parent colony is screened, M individuality before retaining according to ideal adaptation degree size;
(7) certain individuality of Stochastic choice, and a bit in this genes of individuals sequence of Stochastic choice, overturn, obtain colony of future generation;
(8) if t=T, then the individuality in evolutionary process with maximum adaptation degree exports as optimal solution, namely completes the preferred of dynamic reactive compensation device mounting points; If t<T, then repeat (4) ~ (8).
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510794734.6A CN105356481B (en) | 2015-11-18 | 2015-11-18 | A kind of dynamic passive compensation reconnaissance method based on more feed-in short-circuit ratioes |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510794734.6A CN105356481B (en) | 2015-11-18 | 2015-11-18 | A kind of dynamic passive compensation reconnaissance method based on more feed-in short-circuit ratioes |
Publications (2)
Publication Number | Publication Date |
---|---|
CN105356481A true CN105356481A (en) | 2016-02-24 |
CN105356481B CN105356481B (en) | 2018-12-04 |
Family
ID=55332390
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510794734.6A Active CN105356481B (en) | 2015-11-18 | 2015-11-18 | A kind of dynamic passive compensation reconnaissance method based on more feed-in short-circuit ratioes |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105356481B (en) |
Cited By (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105846447A (en) * | 2016-05-19 | 2016-08-10 | 南方电网科学研究院有限责任公司 | Method for determining priority of dynamic reactive power compensation configuration of multi-loop direct current inverter station |
CN105939020A (en) * | 2016-05-14 | 2016-09-14 | 国家电网公司 | Station-configuration method for dynamic reactive compensation apparatus capable of improving multi-feed direct current restoration capability |
CN107104463A (en) * | 2017-06-28 | 2017-08-29 | 南京理工大学 | Current conversion station considers idle active input amount optimization method during black starting-up |
CN107134799A (en) * | 2017-04-19 | 2017-09-05 | 国家电网公司 | A kind of extra-high voltage grid transverter reactive configuration method based on short-circuit ratio |
CN107147107A (en) * | 2017-05-05 | 2017-09-08 | 国电南瑞科技股份有限公司 | A kind of phase modifier points distributing method for suppressing many chain commutation failures of direct current |
CN107482637A (en) * | 2017-08-03 | 2017-12-15 | 中国电力科学研究院 | It is a kind of count and active equipment difference direct-current short circuit than determination method and its system |
CN109167380A (en) * | 2018-10-31 | 2019-01-08 | 上海电力学院 | A kind of judgment method for the multi-infeed systems stability accessing voltage-source type converter station |
CN111446722A (en) * | 2020-03-20 | 2020-07-24 | 国网浙江省电力有限公司温州供电公司 | Measurement method for optimal placement of STATCOMS in multi-feed power electronic system |
CN112054515A (en) * | 2020-08-28 | 2020-12-08 | 武汉大学 | Receiving-end power grid DC receiving capacity detection method based on multi-objective optimization |
Citations (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2006014445A (en) * | 2004-06-24 | 2006-01-12 | Hitachi Ltd | Distribution line voltage fluctuation compensator |
CN103094905A (en) * | 2013-01-07 | 2013-05-08 | 广东电网公司电网规划研究中心 | Selection method of dynamic reactive power compensation configuration point |
CN103258299A (en) * | 2013-04-11 | 2013-08-21 | 国家电网公司 | Net rack optimization method for multi-direct-current centralized fed receiving end power grid |
CN103337864A (en) * | 2013-07-19 | 2013-10-02 | 国家电网公司 | Dynamic reactive compensation mounting point optimization method for multi-direct-current centralized infeed receiving end power grid |
CN103346576A (en) * | 2013-07-03 | 2013-10-09 | 杭州电子科技大学 | Power distribution network reactive compensation node sorting method based on second-order transmission loss sensitivity matrix |
CN103457279A (en) * | 2013-09-23 | 2013-12-18 | 广东电网公司电力调度控制中心 | Determination method for partition D-STATCOM integrated optimization configuration scheme of large-scale power grid |
CN103457275A (en) * | 2013-08-29 | 2013-12-18 | 国家电网公司 | Method for setting up reactive power compensation device model based on controlled alternating current source |
CN103972900A (en) * | 2014-05-21 | 2014-08-06 | 中国南方电网有限责任公司电网技术研究中心 | Method for determining distribution of reactive compensation devices of multi-feed direct current transmission system based on voltage control sensitive factors |
CN104333030A (en) * | 2014-11-03 | 2015-02-04 | 国网宁夏电力公司电力科学研究院 | Analysis method for multi-feed DC interaction factors based on reduced order Jacobian matrix |
CN104333019A (en) * | 2014-11-27 | 2015-02-04 | 宁波华浙电力技术有限公司 | Voltage optimizing point sorting method for 10KV overhead line reactive compensation device |
-
2015
- 2015-11-18 CN CN201510794734.6A patent/CN105356481B/en active Active
Patent Citations (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2006014445A (en) * | 2004-06-24 | 2006-01-12 | Hitachi Ltd | Distribution line voltage fluctuation compensator |
CN103094905A (en) * | 2013-01-07 | 2013-05-08 | 广东电网公司电网规划研究中心 | Selection method of dynamic reactive power compensation configuration point |
CN103258299A (en) * | 2013-04-11 | 2013-08-21 | 国家电网公司 | Net rack optimization method for multi-direct-current centralized fed receiving end power grid |
CN103346576A (en) * | 2013-07-03 | 2013-10-09 | 杭州电子科技大学 | Power distribution network reactive compensation node sorting method based on second-order transmission loss sensitivity matrix |
CN103337864A (en) * | 2013-07-19 | 2013-10-02 | 国家电网公司 | Dynamic reactive compensation mounting point optimization method for multi-direct-current centralized infeed receiving end power grid |
CN103457275A (en) * | 2013-08-29 | 2013-12-18 | 国家电网公司 | Method for setting up reactive power compensation device model based on controlled alternating current source |
CN103457279A (en) * | 2013-09-23 | 2013-12-18 | 广东电网公司电力调度控制中心 | Determination method for partition D-STATCOM integrated optimization configuration scheme of large-scale power grid |
CN103972900A (en) * | 2014-05-21 | 2014-08-06 | 中国南方电网有限责任公司电网技术研究中心 | Method for determining distribution of reactive compensation devices of multi-feed direct current transmission system based on voltage control sensitive factors |
CN104333030A (en) * | 2014-11-03 | 2015-02-04 | 国网宁夏电力公司电力科学研究院 | Analysis method for multi-feed DC interaction factors based on reduced order Jacobian matrix |
CN104333019A (en) * | 2014-11-27 | 2015-02-04 | 宁波华浙电力技术有限公司 | Voltage optimizing point sorting method for 10KV overhead line reactive compensation device |
Cited By (12)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105939020A (en) * | 2016-05-14 | 2016-09-14 | 国家电网公司 | Station-configuration method for dynamic reactive compensation apparatus capable of improving multi-feed direct current restoration capability |
CN105846447A (en) * | 2016-05-19 | 2016-08-10 | 南方电网科学研究院有限责任公司 | Method for determining priority of dynamic reactive power compensation configuration of multi-loop direct current inverter station |
CN107134799A (en) * | 2017-04-19 | 2017-09-05 | 国家电网公司 | A kind of extra-high voltage grid transverter reactive configuration method based on short-circuit ratio |
CN107147107A (en) * | 2017-05-05 | 2017-09-08 | 国电南瑞科技股份有限公司 | A kind of phase modifier points distributing method for suppressing many chain commutation failures of direct current |
CN107147107B (en) * | 2017-05-05 | 2020-05-15 | 国电南瑞科技股份有限公司 | Phase modulator point distribution method for inhibiting multi-direct-current cascading commutation failure |
CN107104463A (en) * | 2017-06-28 | 2017-08-29 | 南京理工大学 | Current conversion station considers idle active input amount optimization method during black starting-up |
CN107104463B (en) * | 2017-06-28 | 2019-10-18 | 南京理工大学 | Converter station considers idle active input amount optimization method during black starting-up |
CN107482637A (en) * | 2017-08-03 | 2017-12-15 | 中国电力科学研究院 | It is a kind of count and active equipment difference direct-current short circuit than determination method and its system |
CN107482637B (en) * | 2017-08-03 | 2022-06-03 | 中国电力科学研究院 | Method and system for determining difference direct current short circuit ratio of active equipment |
CN109167380A (en) * | 2018-10-31 | 2019-01-08 | 上海电力学院 | A kind of judgment method for the multi-infeed systems stability accessing voltage-source type converter station |
CN111446722A (en) * | 2020-03-20 | 2020-07-24 | 国网浙江省电力有限公司温州供电公司 | Measurement method for optimal placement of STATCOMS in multi-feed power electronic system |
CN112054515A (en) * | 2020-08-28 | 2020-12-08 | 武汉大学 | Receiving-end power grid DC receiving capacity detection method based on multi-objective optimization |
Also Published As
Publication number | Publication date |
---|---|
CN105356481B (en) | 2018-12-04 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN105356481A (en) | Multi-infeed-short-circuit-ratio-based dynamic reactive compensation point selection method | |
CN102570477B (en) | SVC and TCSC based on linearisation power flow algorithm model mixes planing method of layouting | |
CN107134799A (en) | A kind of extra-high voltage grid transverter reactive configuration method based on short-circuit ratio | |
CN104362642B (en) | Dynamic reactive reserved optimizing method for improving long-term voltage stabilization in AC/DC (Alternating Current/Direct Current) power grid | |
CN102856918A (en) | Power distribution network reactive power optimization method based on ecological niche particle swarm algorithm | |
CN102403716A (en) | Dynamic equalizing method for multi-infeed alternating/direct-current power grid | |
CN102522756B (en) | Inductive reactive compensation method for power grid for avoiding voltage off-normal risks | |
CN101609989B (en) | System for calculating power supply abundance of urban power network | |
CN104466984A (en) | Dynamic reactive standby optimization method for increasing safety level of direct current commutation | |
CN104598725A (en) | Evaluation method for power distribution network adaptability of distributed power supply after grid connection | |
CN107910887A (en) | A kind of black-start method for considering HVDC transmission system and participating in | |
CN107659194A (en) | A kind of optimal control collection model predictive control method of Modular multilevel converter | |
CN104899396A (en) | Fast decoupled flow calculation method of modified coefficient matrix | |
Barker et al. | Efficient ANFIS-driven power extraction and control strategies for PV-bess integrated electric vehicle charging station | |
CN109377020B (en) | Power transmission network planning method considering load transfer capacity of power distribution network | |
CN108599974B (en) | Micro-grid distributed communication topology design method based on graph theory connectivity | |
Bayati et al. | Accurate modeling of DC microgrid for fault and protection studies | |
CN103515964A (en) | Reactive compensation control method and reactive compensation control device | |
CN104009470B (en) | Electric power networks fault simulation method based on AC power flow | |
CN103208804B (en) | Based on the micro-capacitance sensor reactive power compensation point selection method of branch voltage stability analysis | |
CN110096767B (en) | Cascading failure simulation method for alternating current-direct current series-parallel power grid | |
CN107196313A (en) | A kind of extra-high-voltage alternating current power network high voltage reactor Optimal Configuration Method | |
CN104393777B (en) | Half-bridge module multi-level converter submodule voltage control method | |
CN105356457A (en) | Standby electrical network accident recovery capacity adequacy evaluation method | |
CN107482624B (en) | A kind of grid loss control strategy based on UPFC |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |