CN105356481B - A kind of dynamic passive compensation reconnaissance method based on more feed-in short-circuit ratioes - Google Patents
A kind of dynamic passive compensation reconnaissance method based on more feed-in short-circuit ratioes Download PDFInfo
- Publication number
- CN105356481B CN105356481B CN201510794734.6A CN201510794734A CN105356481B CN 105356481 B CN105356481 B CN 105356481B CN 201510794734 A CN201510794734 A CN 201510794734A CN 105356481 B CN105356481 B CN 105356481B
- Authority
- CN
- China
- Prior art keywords
- indicate
- bus
- node
- compensation device
- voltage
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
- 238000000034 method Methods 0.000 title claims abstract description 28
- 238000009434 installation Methods 0.000 claims abstract description 32
- 230000000694 effects Effects 0.000 claims abstract description 22
- 239000011159 matrix material Substances 0.000 claims abstract description 13
- 230000002068 genetic effect Effects 0.000 claims abstract description 10
- 108090000623 proteins and genes Proteins 0.000 claims description 15
- 230000008859 change Effects 0.000 claims description 14
- 238000002347 injection Methods 0.000 claims description 12
- 239000007924 injection Substances 0.000 claims description 12
- 230000005611 electricity Effects 0.000 claims description 9
- 230000006978 adaptation Effects 0.000 claims description 6
- 239000000243 solution Substances 0.000 claims description 5
- 238000006243 chemical reaction Methods 0.000 claims description 4
- 239000003990 capacitor Substances 0.000 claims description 3
- 230000010429 evolutionary process Effects 0.000 claims description 3
- 230000001052 transient effect Effects 0.000 claims description 3
- 235000006508 Nelumbo nucifera Nutrition 0.000 claims 1
- 240000002853 Nelumbo nucifera Species 0.000 claims 1
- 235000006510 Nelumbo pentapetala Nutrition 0.000 claims 1
- 238000005457 optimization Methods 0.000 abstract description 5
- 230000008569 process Effects 0.000 abstract description 2
- 238000011217 control strategy Methods 0.000 description 3
- 230000008901 benefit Effects 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 238000005259 measurement Methods 0.000 description 2
- 239000000203 mixture Substances 0.000 description 2
- 230000003068 static effect Effects 0.000 description 2
- 238000012885 constant function Methods 0.000 description 1
- 230000007423 decrease Effects 0.000 description 1
- 230000007812 deficiency Effects 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 238000011084 recovery Methods 0.000 description 1
- 230000009467 reduction Effects 0.000 description 1
- 230000004044 response Effects 0.000 description 1
- 230000001568 sexual effect Effects 0.000 description 1
- 238000004088 simulation Methods 0.000 description 1
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 present invention provides a kind of dynamic passive compensation reconnaissance method based on more feed-in short-circuit ratioes, comprising the following steps: establishes Jacobian matrix and solves the voltage influence factor;More feed-in short-circuit ratioes are determined according to the voltage influence factor;Establish each website installation dynamic reactive compensation device effect assessment objective function and corresponding constraint condition;Using the preferred dynamic reactive compensation device installation point of genetic algorithm.The present invention considers the reciprocal effect between multiple-circuit line route, and the selection of the dynamic reactive compensation device installation point in feed-in areas more for China's direct current provides effective ways and technical support;And influence of each element responds characteristic to system stability in dynamic process is considered, reconnaissance method is more realistic, has engineering application value;The present invention carries out dynamic reconnaissance optimization using genetic algorithm, and computational efficiency is high.
Description
Technical field
The invention belongs to technical field of power systems, and in particular to a kind of dynamic passive compensation based on more feed-in short-circuit ratioes
Reconnaissance method.
Background technique
The load centers areas such as China East China Yangtze River Delta are limited by the development of local normal power supplies, and more direct currents will be presented in future
Feed-in structure is concentrated, the system failure may cause multiple-circuit line and chain reaction occurs, and cause Voltage-stabilizing Problems.Domestic external power grid
Operating experience shows all being related to REACTIVE POWER/VOLTAGE CONTROL problem, and dynamic reactive compensation device is more satisfactory solution party
Case, dynamic reactive compensation device control strategy is flexible, fast response time, and the Voltage-stabilizing Problems of multi-infeed DC receiving end power grid can
It is considered as dynamic reactive compensation device, improves the voltage support ability of extra-high voltage grid, is especially disturbed in system failure
After dynamic, the fast quick-recovery for realizing system weak spot voltage is helped.
In view of receiving end Net Frame of Electric Network is intensive, land resource scarcity is installed dynamic reactive compensation device in substation and is held
Vulnerable to the limitation of site, so, dynamic passive compensation allocation plan can should be selected in the case where reaching effect same
Capacity minimum or the least scheme of configuration site, need to determine by optimizing.
More feed-in short-circuit ratioes are relatively common for measuring exchange in multi-infeed DC power grid for DC support ability
Index, be widely used, but it is defined as static index, be only capable of reflection network topology structure, can not count and dynamic element model.
And stability of power system focuses on the system performance in dynamic process under each element responds, various dynamic elements such as generator mould
Type, load model, DC control model, dynamic passive compensation model etc. are affected to voltage stability.Currently, both at home and abroad
The weakness zone that existing dynamic reactive compensation device configuration method is generally basede on static voltage stability carries out reconnaissance, not yet comprehensive
Consider multi-infeed DC system performance and dynamic element model characteristics.
Summary of the invention
In order to overcome the above-mentioned deficiencies of the prior art, the present invention provides a kind of dynamic reactive benefit based on more feed-in short-circuit ratioes
Reconnaissance method is repaid, reconnaissance method is more realistic, has engineering application value.
In order to achieve the above-mentioned object of the invention, the present invention adopts the following technical scheme that:
The present invention provides a kind of dynamic passive compensation reconnaissance method based on more feed-in short-circuit ratioes, and the method includes following
Step:
Step 1: establishing Jacobian matrix and solve the voltage influence factor;
Step 2: more feed-in short-circuit ratioes are determined according to the voltage influence factor;
Step 3: establishing each website installation dynamic reactive compensation device effect assessment objective function and corresponding constraint condition;
Step 4: using the preferred dynamic reactive compensation device installation point of genetic algorithm.
The step 1 the following steps are included:
Step 1-1: consider dynamic element model, establish power balance equation;
Step 1-2: establishing Jacobian matrix, and solves the voltage influence factor.
In the step 1-1, dynamic element model includes generator model, load model, DC control model and dynamic
Reactive power compensation device model.
In the step 1-1, following power balance equation is established:
Wherein, Δ Pi、ΔQiRespectively indicate the active power variable quantity and reactive power variable quantity of node i injection, PGi、QGi
Respectively indicate the active power output and idle power output of generator injection node i, PLi、QLiRespectively indicate the burden with power and nothing of node i
Workload, PDiIndicate the dc power of node i, QDiIndicate the reactive power of DC filter capacitor injection node i, Ui、UjRespectively
Indicate the voltage of node i, j, QSiIndicate the idle power output of dynamic reactive compensation device injection node i, Gij、BijRespectively indicate section
Conductance and susceptance between point i, j, θijPhase difference of voltage between expression node i, j, i=1,2 ... ..., n, j=1,
2 ... ..., n, n are node total number.
In the step 1-2, following Jacobian matrix equation is established:
Wherein, MIIFlkIndicate the voltage influence factor of bus l opposing busbars k, and MIIFlk=Δ Ul/ΔUk, Δ UlTable
Show the voltage variety of bus l, Δ UkIndicate the voltage variety of bus k;
Jacobian matrix element Hii、Nii、Mii、LiiIt calculates according to the following formula:
Wherein, GiiIndicate the conductance of node i, BiiIndicate the susceptance of node i;
It respectively indicates are as follows:
Wherein, Ei" indicate generator potential, θδiIndicate Ei" and UiPhase angle difference, X "diIndicate the super transient state electricity of generator d axis
It is anti-;If node i is constant current load bus,It respectively indicates are as follows:
Wherein, IPi、IQiRespectively indicate the watt current and reactive current of constant current load bus;
If node i is constant impedance load bus,It respectively indicates are as follows:
Wherein, Gi、BiRespectively indicate the conductance and susceptance of constant impedance load bus;
If node i is DC line access node,It respectively indicates are as follows:
Wherein, IdIndicate DC current, ntIndicate six pulse conversion devices series connection number, kTIndicate converter power transformer no-load voltage ratio, kγ
Indicate the equivalent no-load voltage ratio of converter power transformer, θdIndicate the direct current angle of overlap of rectification side or the blow-out angle of inverter side, XcIndicate equivalent commutation
Reactance;Indicate Equivalent Power Factor angle, and
If node i is that dynamic reactive compensation device installs node,It indicates are as follows:
Wherein, BiIndicate the susceptance of dynamic passive compensation installation node, and Bi=-K Δ Ui=-K (Ui-Ui0), Δ UiIt indicates
Voltage deviation before and after dynamic reactive compensation device, U are installedi0Indicate the initial voltage of dynamic reactive compensation device installation node, K
Indicate proportionality coefficient.
In the step 2, the drop point that pth returns DC line is bus l, and the drop point that q returns DC line is bus k, root
More feed-in short-circuit ratioes are determined according to the voltage influence factor, are had:
Wherein, MISCR 'pIndicate that pth returns more feed-in short-circuit ratioes of DC line, SlIndicate the system short circuit capacity of bus l,
PpIndicate that pth returns the power of DC line, PqIndicate that q returns the power of DC line, MIIFlkIndicate bus l opposing busbars k's
The voltage influence factor, and MIIFlk=Δ Ul/ΔUk, Δ UlIndicate the voltage variety of bus l, Δ UkIndicate the voltage of bus k
Variable quantity, m indicate that returning for DC line counts in the more feed-in power grids of direct current.
In the step 3, each website installation dynamic reactive compensation device effect assessment objective function is established, is had:
Wherein, f indicates that each website installs dynamic reactive compensation device effect assessment objective function,Indicate that pth returns direct current
The equivalent weight coefficient of route, has:
Wherein, ωpIndicate that pth returns the weight coefficient of DC line, ωqIndicate that q returns the weight coefficient of DC line,
ωp、ωqReflect that pth, q return influence of the DC line to other direct currents in the more feed-in power grids of direct current respectively, have:
Wherein, ZpqIndicate the equivalent resistance that pth is returned between DC line change of current bus and q times DC line change of current buses
It is anti-, ZppIndicate that pth returns the equivalent impedance of DC line change of current bus, ZqqIndicate that q returns the equivalent resistance of DC line change of current bus
It is anti-;
It is as follows that each website installs the corresponding constraint condition of dynamic reactive compensation device effect assessment objective function additional:
Wherein, Pl、QlRespectively indicate the active power and reactive power of bus l, Ul、UkRespectively indicate the electricity of bus l, k
Pressure,WithRespectively indicate the upper voltage limit and lower limit of bus k, Glk、BlkRespectively indicate conductance between bus l, k and
Susceptance, θlkIndicate the phase difference of voltage between bus l, k, IkIndicate the short circuit current of bus k,Indicate the short circuit of bus k
Upper current limit, SrIndicate the power of branch r,Indicate the upper limit of the power of branch r.
In the step 4, using the preferred dynamic reactive compensation device installation point of genetic algorithm, comprising:
(1) setting evolutionary generation t is 0, and maximum evolutionary generation T and individual in population number M is arranged;
(2) it selects s node as dynamic reactive compensation device installation point in n node, then sharesIt is a to deposit
Dynamic reactive compensation device reconnaissance strategy, and to reactive power compensator reconnaissance strategy carry out gene coding, gene coding length
L is spent to meetThe value of each bit of gene code length is 0 or 1;
(3) random to generate M individual as initial population P0, each dynamic reactive compensation device reconnaissance strategy is as one
Individual;
(4) fitness of each individual is calculated, i.e., each website installs dynamic reactive compensation device effect assessment objective function
Value;
(5) distribute a random number for each parent individuality, and according to corresponding random number to parent individuality according to from big
It being ranked up to small sequence, two neighboring parent individuality is hybridized, some bit in random selection gene order, and two
The 0 of bit or 1 are interchangeable by parent individuality, are generated offspring individual, are calculated the fitness of offspring individual, and offspring individual is added
Enter into parent individuality composition parent group;
(6) parent group is screened, M individual before retaining according to individual adaptation degree size;
(7) some individual is randomly choosed, and randomly chooses a bit in the genes of individuals sequence, is overturn, is obtained
Next-generation group;
(8) if t=T, the individual using in evolutionary process with maximum adaptation degree is exported as optimal solution, that is, completes dynamic
Reactive power compensator installation point it is preferred;If t < T, (4)~(8) are repeated.
Compared with the immediate prior art, technical solution provided by the invention is had the advantages that
1. the present invention is using more feed-in short-circuit ratioes as the measurement index of idle investment, dynamic in feed-in power grids more for direct current
State reactive power compensator installation point carries out preferred, it is contemplated that and the reciprocal effect between multiple-circuit line route is more for China's direct current
The selection of the dynamic reactive compensation device installation point in feed-in area provides effective ways and technical support;
2. the present invention comprehensively considers multi-infeed DC system performance and dynamic element model characteristics, by generator, load, straight
The dynamic elements model such as stream, dynamic reactive compensation device is brought into the measurement index of more feed-in short-circuit ratioes, it is contemplated that dynamic mistake
Influence of each element responds characteristic to system stability in journey, reconnaissance method is more realistic, has engineering application value;
3. the present invention carries out dynamic reconnaissance optimization using genetic algorithm, it is not based on time-domain-simulation, computational efficiency is high.
Detailed description of the invention
Fig. 1 is the dynamic passive compensation reconnaissance method flow diagram in the embodiment of the present invention based on more feed-in short-circuit ratioes;
Fig. 2 is direct current more feed-in power grid geography wiring schematic diagrams in the year two thousand twenty East China in the embodiment of the present invention 2.
Specific embodiment
The present invention is described in further detail below in conjunction with the accompanying drawings.
The present invention provides a kind of dynamic passive compensation reconnaissance method based on more feed-in short-circuit ratioes, such as Fig. 1, the method packet
Include following steps:
Step 1: establishing Jacobian matrix and solve the voltage influence factor;
Step 2: more feed-in short-circuit ratioes are determined according to the voltage influence factor;
Step 3: establishing each website installation dynamic reactive compensation device effect assessment objective function and corresponding constraint condition;
Step 4: using the preferred dynamic reactive compensation device installation point of genetic algorithm.
In the step 1, dynamic element model includes generator model, load model, DC control model and dynamic nothing
Reactive power compensation installations model.
The step 1 the following steps are included:
Step 1-1: consider dynamic element model, establish power balance equation;
Step 1-2: establishing Jacobian matrix, and solves the voltage influence factor.
In the step 1-1, following power balance equation is established:
Wherein, Δ Pi、ΔQiRespectively indicate the active power variable quantity and reactive power variable quantity of node i injection, PGi、QGi
Respectively indicate the active power output and idle power output of generator injection node i, PLi、QLiRespectively indicate the burden with power and nothing of node i
Workload, PDiIndicate the dc power of node i, QDiIndicate the reactive power of DC filter capacitor injection node i, Ui、UjRespectively
Indicate the voltage of node i, j, QSiIndicate the idle power output of dynamic reactive compensation device injection node i, Gij、BijRespectively indicate section
Conductance and susceptance between point i, j, θijPhase difference of voltage between expression node i, j, i=1,2 ... ..., n, j=1,
2 ... ..., n, n are node total number.
In the step 1-2, following Jacobian matrix equation is established:
Wherein, MIIFlkIndicate the voltage influence factor of bus l opposing busbars k, and MIIFlk=Δ Ul/ΔUk, Δ UlTable
Show the voltage variety of bus l, Δ UkIndicate the voltage variety of bus k;
Jacobian matrix element Hii、Nii、Mii、LiiIt calculates according to the following formula:
Wherein, GiiIndicate the conductance of node i, BiiIndicate the susceptance of node i;
It respectively indicates are as follows:
Wherein, Ei" indicate generator potential, θδiIndicate Ei" and UiPhase angle difference, X "diIndicate the super transient state electricity of generator d axis
It is anti-;
If node i is constant current load bus,It respectively indicates are as follows:
Wherein, IPi、IQiRespectively indicate the watt current and reactive current of constant current load bus;
If node i is constant impedance load bus,It respectively indicates are as follows:
Wherein, Gi、BiRespectively indicate the conductance and susceptance of constant impedance load bus;
If node i is DC line power node,It respectively indicates are as follows:
Wherein, IdIndicate DC current, ntIndicate six pulse conversion devices series connection number, kTIndicate converter power transformer no-load voltage ratio, kγ
Indicate the equivalent no-load voltage ratio of converter power transformer, θdIndicate the direct current angle of overlap of rectification side or the blow-out angle of inverter side, XcIndicate equivalent commutation
Reactance;Indicate Equivalent Power Factor angle, and
If node i is that dynamic reactive compensation device installs node,It indicates are as follows:
Wherein, BiIndicate the susceptance of dynamic passive compensation installation node, and Bi=-K Δ Ui=-K (Ui-Ui0), Δ UiIt indicates
Voltage deviation before and after dynamic reactive compensation device, U are installedi0Indicate the initial voltage of dynamic reactive compensation device installation node, K
Indicate proportionality coefficient.
In the step 2, the drop point that pth returns DC line is bus l, and the drop point that q returns DC line is bus k, root
More feed-in short-circuit ratioes are determined according to the voltage influence factor, are had:
Wherein, MISCR 'pIndicate that pth returns more feed-in short-circuit ratioes of DC line, SlIndicate the system short circuit capacity of bus l,
PpIndicate that pth returns the power of DC line, PqIndicate that q returns the power of DC line, MIIFlkIndicate bus l opposing busbars k's
The voltage influence factor, and MIIFlk=Δ Ul/ΔUk, Δ UlIndicate the voltage variety of bus l, Δ UkIndicate the voltage of bus k
Variable quantity, m indicate that returning for DC line counts in the more feed-in power grids of direct current.
In the step 3, each website installation dynamic reactive compensation device effect assessment objective function is established, is had:
Wherein, f indicates that each website installs dynamic reactive compensation device effect assessment objective function,Indicate that pth returns direct current
The equivalent weight coefficient of route, has:
Wherein, ωpIndicate that pth returns the weight coefficient of DC line, ωqIndicate that q returns the weight coefficient of DC line,
ωp、ωqReflect that pth, q return influence of the DC line to other direct currents in the more feed-in power grids of direct current respectively, have:
Wherein, ZpqIndicate the equivalent resistance that pth is returned between DC line change of current bus and q times DC line change of current buses
It is anti-, ZppIndicate that pth returns the equivalent impedance of DC line change of current bus, ZqqIndicate that q returns the equivalent resistance of DC line change of current bus
It is anti-;
It is as follows that each website installs the corresponding constraint condition of dynamic reactive compensation device effect assessment objective function additional:
Wherein, Pl、QlRespectively indicate the active power and reactive power of bus l, Ul、UkRespectively indicate the electricity of bus l, k
Pressure,WithRespectively indicate the upper voltage limit and lower limit of bus k, Glk、BlkRespectively indicate conductance between bus l, k and
Susceptance, θlkIndicate the phase difference of voltage between bus l, k, IkIndicate the short circuit current of bus k,Indicate the short circuit of bus k
Upper current limit, SrIndicate the power of branch r,Indicate the upper limit of the power of branch r.
In the step 4, using the preferred dynamic reactive compensation device installation point of genetic algorithm, comprising:
(1) setting evolutionary generation t is 0, and maximum evolutionary generation T and individual in population number M is arranged;
(2) it selects s node as dynamic reactive compensation device installation point in n node, then sharesIt is a to deposit
Dynamic reactive compensation device reconnaissance strategy, and to reactive power compensator reconnaissance strategy carry out gene coding, gene coding length
L is spent to meetThe value of each bit of gene code length is 0 or 1;
(3) random to generate M individual as initial population P0, each dynamic reactive compensation device reconnaissance strategy is as one
Individual;
(4) fitness of each individual is calculated, i.e., each website installs dynamic reactive compensation device effect assessment objective function
Value;
(5) distribute a random number for each parent individuality, and according to corresponding random number to parent individuality according to from big
It being ranked up to small sequence, two neighboring parent individuality is hybridized, some bit in random selection gene order, and two
The 0 of bit or 1 are interchangeable by parent individuality, are generated offspring individual, are calculated the fitness of offspring individual, and offspring individual is added
Enter into parent individuality composition parent group;
(6) parent group is screened, M individual before retaining according to individual adaptation degree size;
(7) some individual is randomly choosed, and randomly chooses a bit in the genes of individuals sequence, is overturn, is obtained
Next-generation group;
(8) if t=T, the individual using in evolutionary process with maximum adaptation degree is exported as optimal solution, that is, completes dynamic
Reactive power compensator installation point it is preferred;If t < T, (4)~(8) are repeated.
Embodiment 1
Reconnaissance method provided by the invention is applied to present East China Power Grid, carries out more feed-in short-circuit ratioes and calculates.If with
Traditional reconnaissance method acquired results are MISCRp, reconnaissance method acquired results provided by the invention are MISCRp'.East China Power Grid feedback
Enter direct current 9 to return, shown in the more feed-in short-circuit ratio calculated results 1 of East China Power Grid direct current in 2016:
Table 1
In table 1, MISCR' does not account for hvdc control mode, and load model is " the constant function of 40% constant impedance+60%
Rate ".
As can be seen from the results in the table that two class method calculated results have different, but according to more feed-in short-circuit ratio sizes
Sequence substantially there is no variation.
To same 9 times direct currents, the more feed-in short-circuit ratioes calculated under different DC control strategies are calculated, East China Power Grid in 2016
More feed-in short-circuit ratio calculated result such as tables 2 under direct current difference control strategy:
Table 2
2 calculated result of table shows: in most cases, more feed-in short-circuit ratio ratios under constant current control mode ignore direct current control
It is bigger when processed;The more feed-in short-circuit ratioes of constant dc power control mode are small.This is because, constant dc power control mode supports reactive voltage
The reason of situation is deteriorated is: when voltage drop is low on inverter side change of current bus, inverter makes blow-out angle by adjusting Trigger Angle
It remains unchanged, inverter side DC voltage reduces, so as to cause DC current raising;Due to the influence of commutating reactance, inverter side is straight
Galvanic electricity pressure further decreases, and leads to the reduction of inverter side converter power factor, and reactive requirement increases.For rectification side, due to straight
Galvanic electricity stream increases, in order to keep power invariability, it is necessary to DC voltage is reduced, to also reactive requirement be made to increase.
Illustrated by above-mentioned calculated result, the more feed-in short-circuit ratio calculated results and former short-circuit ratio derived using the method for the present invention
Calculated result is consistent sexual intercourse, reflects simultaneously, and the model of all kinds of dynamic elements is different, to the meter of index and its stability
Certain influence will be generated by calculating result.
Embodiment 2
Reconnaissance method provided by the invention is applied to the year two thousand twenty Jiangsu planning power grid, as shown in Fig. 2, shared political affairs are flat, same
In, Liyang, Taizhou, Nanjing, Changshu, 7, Changzhou direct current drop point, form typical multi-infeed DC power grid, Jiangsu Power Grid is by electricity
Ratio 40%.Since feed-in direct current scale is excessive, when not installing dynamic reactive compensation device, there are 14 backcrossing Flow Lines in Jiangsu
Three-phase permanent short failure will lead to Voltage Instability.
The 500kV substation of southern area of Jiangsu Province 25 alternately point is chosen, as shown in table 3, passes through procedure selection 5 stations of optimization
The SVC of point installation 2 × 240Mvar capacity;
Table 3
Different evolutionary generation T, group individual amount M are set, and optimizes calculating, southern Jiangsu dynamic passive compensation reconnaissance
Scenario outcomes are as shown in table 4, from the point of view of calculated result, influenced not using different evolutionary generation and individual amount on calculated result
Greatly, calculated result can stablize the scheme relatively determined at one substantially.Work as T=30, when M=30, result and remaining scheme
Slightly difference, but from the point of view of target function value, the two difference very little.
Table 4
By result in table it is found that solving the above problems using genetic algorithm, preferable convergence can be reached.To above-mentioned two
A different reactive compensation allocation plan carries out analysis safety and stability and checks, the results show that N-1 leads to system under two class schemes
The route of unstability is 7 times, obvious to the promotion effect of system stability.
It can be seen that the dynamic reactive compensation device optimization reconnaissance method in the more feed-in areas of direct current mentioned herein calculates effect
Rate is high, and effect of optimization is good, has very strong engineering application value.
Finally it should be noted that: the above embodiments are merely illustrative of the technical scheme of the present invention and are not intended to be limiting thereof, institute
The those of ordinary skill in category field can still modify to a specific embodiment of the invention referring to above-described embodiment or
Equivalent replacement, these are applying for this pending hair without departing from any modification of spirit and scope of the invention or equivalent replacement
Within bright claims.
Claims (2)
1. a kind of dynamic passive compensation reconnaissance method based on more feed-in short-circuit ratioes, it is characterised in that: the method includes following
Step:
Step 1: establishing Jacobian matrix and solve the voltage influence factor;
Step 2: more feed-in short-circuit ratioes are determined according to the voltage influence factor;
Step 3: establishing each website installation dynamic reactive compensation device effect assessment objective function and corresponding constraint condition;
Step 4: using the preferred dynamic reactive compensation device installation point of genetic algorithm;
The step 1 the following steps are included:
Step 1-1: consider dynamic element model, establish power balance equation;
Step 1-2: establishing Jacobian matrix, and solves the voltage influence factor;
In the step 1-1, dynamic element model includes generator model, load model, DC control model and dynamic reactive
Compensation device model;
In the step 1-1, following power balance equation is established:
Wherein, Δ Pi、ΔQiRespectively indicate the active power variable quantity and reactive power variable quantity of node i injection, PGi、QGiRespectively
Indicate the active power output and idle power output of generator injection node i, PLi、QLiRespectively indicate the burden with power of node i and idle negative
Lotus, PDiIndicate the dc power of node i, QDiIndicate the reactive power of DC filter capacitor injection node i, Ui、UjIt respectively indicates
The voltage of node i, j, QSiIndicate the idle power output of dynamic reactive compensation device injection node i, Gij、BijRespectively indicate node i, j
Between conductance and susceptance, θijIndicate the phase difference of voltage between node i, j, i=1,2 ... ..., n, j=1,2 ... ..., n, n
For node total number;
In the step 1-2, following Jacobian matrix equation is established:
Wherein, MIIFlkIndicate the voltage influence factor of bus l opposing busbars k, and MIIFlk=Δ Ul/ΔUk, Δ UlIndicate bus
The voltage variety of l, Δ UkIndicate the voltage variety of bus k;
Jacobian matrix element Hii、Nii、Mii、LiiIt calculates according to the following formula:
Wherein, GiiIndicate the conductance of node i, BiiIndicate the susceptance of node i;
It respectively indicates are as follows:
Wherein, E "iIndicate generator potential, θδiIndicate E "iWith UiPhase angle difference, X "diIndicate the super transient state reactance of generator d axis;If
Node i is constant current load bus,It respectively indicates are as follows:
Wherein, IPi、IQiRespectively indicate the watt current and reactive current of constant current load bus;
If node i is constant impedance load bus,It respectively indicates are as follows:
Wherein, Gi、BiRespectively indicate the conductance and susceptance of constant impedance load bus;
If node i is DC line access node,It respectively indicates are as follows:
Wherein, IdIndicate DC current, ntIndicate six pulse conversion devices series connection number, kTIndicate converter power transformer no-load voltage ratio, kγIt indicates
The equivalent no-load voltage ratio of converter power transformer, θdIndicate the direct current angle of overlap of rectification side or the blow-out angle of inverter side, XcIndicate equivalent commutation electricity
It is anti-;Indicate Equivalent Power Factor angle, and
If node i is that dynamic reactive compensation device installs node,It indicates are as follows:
Wherein, BiIndicate the susceptance of dynamic passive compensation installation node, and Bi=-K Δ Ui=-K (Ui-Ui0), Δ Ui indicates installation
Voltage deviation before and after dynamic reactive compensation device, Ui0Indicate that the initial voltage of dynamic reactive compensation device installation node, K indicate
Proportionality coefficient;
In the step 2, the drop point that pth returns DC line is bus l, and the drop point that q returns DC line is bus k, according to electricity
Pressure impact factor determines more feed-in short-circuit ratioes, has:
Wherein, MISCRp' indicate that pth returns more feed-in short-circuit ratioes of DC line, SlIndicate the system short circuit capacity of bus l, PpTable
Show that pth returns the power of DC line, PqIndicate that q returns the power of DC line, MIIFlkIndicate the voltage of bus l opposing busbars k
Impact factor, and MIIFlk=Δ Ul/ΔUk, Δ UlIndicate the voltage variety of bus l, Δ UkIndicate the voltage change of bus k
Amount, m indicate that returning for DC line counts in the more feed-in power grids of direct current;
In the step 3, each website installation dynamic reactive compensation device effect assessment objective function is established, is had:
Wherein, f indicates that each website installs dynamic reactive compensation device effect assessment objective function,Indicate that pth returns DC line
Equivalent weight coefficient, have:
Wherein, ωpIndicate that pth returns the weight coefficient of DC line, ωqIndicate that q returns the weight coefficient of DC line, ωp、ωq
Reflect that pth, q return influence of the DC line to other direct currents in the more feed-in power grids of direct current respectively, have:
Wherein, ZpqIndicate the equivalent impedance that pth is returned between DC line change of current bus and q times DC line change of current buses, Zpp
Indicate that pth returns the equivalent impedance of DC line change of current bus, ZqqIndicate that q returns the equivalent impedance of DC line change of current bus;
It is as follows that each website installs the corresponding constraint condition of dynamic reactive compensation device effect assessment objective function additional:
Wherein, Pl、QlRespectively indicate the active power and reactive power of bus l, Ul、UkThe voltage of bus l, k are respectively indicated,
WithRespectively indicate the upper voltage limit and lower limit of bus k, Glk、BlkRespectively indicate the conductance and susceptance between bus l, k, θlk
Indicate the phase difference of voltage between bus l, k, IkIndicate the short circuit current of bus k,Indicate the short circuit current upper limit of bus k,
SrIndicate the power of branch r,Indicate the upper limit of the power of branch r.
2. the dynamic passive compensation reconnaissance method according to claim 1 based on more feed-in short-circuit ratioes, it is characterised in that: institute
It states in step 4, using the preferred dynamic reactive compensation device installation point of genetic algorithm, comprising:
(1) setting evolutionary generation t is 0, and maximum evolutionary generation T and individual in population number M is arranged;
(2) it selects s node as dynamic reactive compensation device installation point in n node, then sharesIt is a that may be present
Dynamic reactive compensation device reconnaissance strategy, and gene coding is carried out to reactive power compensator reconnaissance strategy, gene code length L is full
FootThe value of each bit of gene code length is 0 or 1;
(3) random to generate M individual as initial population P0, each dynamic reactive compensation device reconnaissance strategy is as an individual;
(4) fitness of each individual is calculated, i.e., each website installs dynamic reactive compensation device effect assessment target function value;
(5) for each parent individuality distribute a random number, and according to corresponding random number to parent individuality according to from big to small
Sequence be ranked up, two neighboring parent individuality is hybridized, randomly choose gene order on some bit, two parents
The 0 of bit or 1 are interchangeable by individual, are generated offspring individual, are calculated the fitness of offspring individual, and offspring individual is added to
Parent individuality forms in parent group;
(6) parent group is screened, M individual before retaining according to individual adaptation degree size;
(7) some individual is randomly choosed, and randomly chooses a bit in the genes of individuals sequence, is overturn, is obtained next
For group;
(8) if t=T, the individual using in evolutionary process with maximum adaptation degree is exported as optimal solution, i.e. completion dynamic reactive
Compensation device installation point it is preferred;If t < T, (4)~(8) are repeated.
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 CN105356481A (en) | 2016-02-24 |
CN105356481B true 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) |
Families Citing this family (9)
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 |
CN105846447B (en) * | 2016-05-19 | 2018-09-18 | 南方电网科学研究院有限责任公司 | 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 |
CN107147107B (en) * | 2017-05-05 | 2020-05-15 | 国电南瑞科技股份有限公司 | Phase modulator point distribution method for inhibiting multi-direct-current cascading commutation failure |
CN107104463B (en) * | 2017-06-28 | 2019-10-18 | 南京理工大学 | Converter station considers idle active input amount optimization method during black starting-up |
CN107482637B (en) * | 2017-08-03 | 2022-06-03 | 中国电力科学研究院 | Method and system for determining difference direct current short circuit ratio of active equipment |
CN109167380B (en) * | 2018-10-31 | 2020-04-07 | 上海电力学院 | Method for judging stability of multi-feed-in system accessed to 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 |
CN112054515B (en) * | 2020-08-28 | 2022-12-06 | 武汉大学 | Receiving-end power grid DC receiving capacity evaluation method based on multi-objective optimization |
Family Cites Families (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2006014445A (en) * | 2004-06-24 | 2006-01-12 | Hitachi Ltd | Distribution line voltage fluctuation compensator |
CN103094905B (en) * | 2013-01-07 | 2015-03-25 | 广东电网公司电网规划研究中心 | Selection method of dynamic reactive power compensation configuration point |
CN103258299B (en) * | 2013-04-11 | 2016-09-21 | 国家电网公司 | A kind of many direct currents concentrate the receiving end Net Frame of Electric Network optimization method of feed-in |
CN103346576B (en) * | 2013-07-03 | 2015-04-22 | 杭州电子科技大学 | Power distribution network reactive compensation node sorting method based on second-order transmission loss sensitivity matrix |
CN103337864B (en) * | 2013-07-19 | 2015-11-11 | 国家电网公司 | Many direct currents concentrate the dynamic passive compensation mounting points optimization method of feed-in receiving end electrical network |
CN103457275B (en) * | 2013-08-29 | 2015-08-05 | 国家电网公司 | Based on the modeling method of the reactive power compensation device model of controlled AC 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 |
CN103972900B (en) * | 2014-05-21 | 2016-03-02 | 中国南方电网有限责任公司电网技术研究中心 | Utilize the method that voltage control sensitive factor determination Multi-infeed HVDC transmission system reactive power compensator is layouted |
CN104333030B (en) * | 2014-11-03 | 2016-08-17 | 国网宁夏电力公司电力科学研究院 | A kind of analysis method of multi-infeed HVDC interaction factor |
CN104333019B (en) * | 2014-11-27 | 2016-06-29 | 宁波华浙电力技术有限公司 | 10KV overhead transmission line reactive power compensator voltage optimization reconnaissance method |
-
2015
- 2015-11-18 CN CN201510794734.6A patent/CN105356481B/en active Active
Also Published As
Publication number | Publication date |
---|---|
CN105356481A (en) | 2016-02-24 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN105356481B (en) | A kind of dynamic passive compensation reconnaissance method based on more feed-in short-circuit ratioes | |
CN107039993A (en) | Symmetrical bipolar flexible direct-current transmission converter power turns band control method | |
Papanikolaou | Low‐voltage ride‐through concept in flyback inverter‐based alternating current–photovoltaic modules | |
CN110601198B (en) | Hybrid micro-grid optimized operation method considering harmonic and voltage unbalance constraints | |
CN108808715B (en) | Multi-terminal flexible direct system static safety analysis method considering direct current network fault power | |
CN104362642B (en) | Dynamic reactive reserved optimizing method for improving long-term voltage stabilization in AC/DC (Alternating Current/Direct Current) power grid | |
CN108429252B (en) | Method for calculating contribution short-circuit current of alternating current system during direct current fault of multi-terminal alternating current-direct current hybrid power distribution network | |
CN107276110B (en) | Multi-infeed direct current transmission receiving end alternating current system equivalent evaluation method | |
CN109638839B (en) | Load flow calculation method of bipolar flexible direct-current transmission system | |
CN108599227B (en) | MMC direct-current voltage balance control method for forming direct-current converter station based on MMC cascade connection | |
CN108306324A (en) | Modularization centralized formula energy-storage system | |
CN110212558A (en) | A kind of addressing constant volume configuration method of distributed energy storage system | |
CN104899396A (en) | Fast decoupled flow calculation method of modified coefficient matrix | |
CN107332290B (en) | Regional load transfer method based on direct current circuit | |
CN111242389A (en) | Intelligent energy storage soft switch planning method, system, equipment and medium | |
CN109377020B (en) | Power transmission network planning method considering load transfer capacity of power distribution network | |
Fotoohabadi et al. | Evaluating the technical benefits of AC–DC hybrid distribution systems consisting of solid-state transformers using a multiobjective index | |
CN108988401B (en) | Method for evaluating power transmission capacity and constraint conditions of multi-channel alternating current-direct current series-parallel connection section | |
CN105262077B (en) | Control method for optimizing power flow of direct-current power distribution network | |
CN104462821B (en) | Method for quickly calculating maximum length of collecting lines of wind farm | |
CN109617079A (en) | A kind of existence and method for analyzing stability of direct-flow distribution system flow solution | |
CN108988336B (en) | Optimization planning method for charging pile system with nested micro-grid | |
CN110096767A (en) | A kind of AC-DC hybrid power grid cascading failure emulation mode | |
CN107069703B (en) | AC/DC distribution network planning method considering new energy access | |
CN109921406A (en) | True bipolar flexible DC grid is layered tidal current computing method |
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 |