CN105226703B  Distributed windpowered electricity generation multiobjective planning method based on Intrusion Index and balance technology  Google Patents
Distributed windpowered electricity generation multiobjective planning method based on Intrusion Index and balance technology Download PDFInfo
 Publication number
 CN105226703B CN105226703B CN201510609613.XA CN201510609613A CN105226703B CN 105226703 B CN105226703 B CN 105226703B CN 201510609613 A CN201510609613 A CN 201510609613A CN 105226703 B CN105226703 B CN 105226703B
 Authority
 CN
 China
 Prior art keywords
 mrow
 msub
 dwg
 msubsup
 mfrac
 Prior art date
Links
 238000005516 engineering processes Methods 0.000 title claims abstract description 13
 230000003993 interaction Effects 0.000 claims abstract description 22
 238000004422 calculation algorithm Methods 0.000 claims abstract description 4
 230000000875 corresponding Effects 0.000 claims abstract description 4
 230000001010 compromised Effects 0.000 claims abstract description 3
 230000002068 genetic Effects 0.000 claims abstract description 3
 238000011068 load Methods 0.000 claims description 35
 230000005611 electricity Effects 0.000 claims description 15
 238000004364 calculation method Methods 0.000 claims description 10
 241001183012 Modified Vaccinia Ankara virus Species 0.000 claims description 9
 230000035969 Vmax Effects 0.000 claims description 7
 230000000694 effects Effects 0.000 claims description 4
 230000000670 limiting Effects 0.000 claims description 3
 230000001186 cumulative Effects 0.000 claims description 2
 230000002045 lasting Effects 0.000 claims description 2
 239000011159 matrix material Substances 0.000 claims description 2
 239000012466 permeate Substances 0.000 claims description 2
 238000011084 recovery Methods 0.000 claims description 2
 238000005457 optimization Methods 0.000 description 3
 230000002829 reduced Effects 0.000 description 3
 238000001764 infiltration Methods 0.000 description 2
 238000004088 simulation Methods 0.000 description 2
 235000006508 Nelumbo nucifera Nutrition 0.000 description 1
 240000002853 Nelumbo nucifera Species 0.000 description 1
 235000006510 Nelumbo pentapetala Nutrition 0.000 description 1
 241000220317 Rosa Species 0.000 description 1
 238000004458 analytical method Methods 0.000 description 1
 235000020127 ayran Nutrition 0.000 description 1
 230000015556 catabolic process Effects 0.000 description 1
 230000003247 decreasing Effects 0.000 description 1
 230000002950 deficient Effects 0.000 description 1
 230000001419 dependent Effects 0.000 description 1
 238000009795 derivation Methods 0.000 description 1
 238000005553 drilling Methods 0.000 description 1
 238000010348 incorporation Methods 0.000 description 1
 238000009434 installation Methods 0.000 description 1
 230000002452 interceptive Effects 0.000 description 1
 238000002955 isolation Methods 0.000 description 1
 238000000034 method Methods 0.000 description 1
 230000004048 modification Effects 0.000 description 1
 238000006011 modification reaction Methods 0.000 description 1
 230000000630 rising Effects 0.000 description 1
 238000006467 substitution reaction Methods 0.000 description 1
Classifications

 Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSSSECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSSREFERENCE 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
 Y02E10/00—Energy generation through renewable energy sources
 Y02E10/70—Wind energy
 Y02E10/76—Power conversion electric or electronic aspects
Abstract
The invention discloses the distributed windpowered electricity generation multiobjective planning method based on Intrusion Index and balance technology, step 1：Try to achieve five influence indexs for not accessing DWG；Step 2：The multiple objective function for determining DWG planing methods is and equation and inequality constraints condition；Setting interaction counts r=1；Step 3：One in five influence indexs of selection is used as major heading；Step 4：The corresponding deviation of each target setting is compromisedIt is determined that the limitation of maximum or minimum allowable, using remaining target as constraints；Step 5：Computing genetic algorithm calculates；Step 6：If result can compromise, into step 10；Otherwise, into next step；Step 7：Interaction counts increase by 1, r=r+1；Step 8：New deviation is set to compromiseAnd repeat step 5  step 8；Step 9：If major heading needs to change, return to step four；Step 10：Try to achieve most compromise solution.The present invention considers each target component, and designer can make decisions on one's own compromise nondomination solution.
Description
Technical field
The present invention relates to a kind of distributed windpowered electricity generation multiobjective planning method based on Intrusion Index and balance technology.
Background technology
In order to reduce power attenuation and peak value operating cost, voltage is preferably adjusted, improves integrality, the reliability of system
And efficiency, distributed windpower generator (distributed wind gerneration, DWG) obtained in recent years in power distribution network compared with
Used to be extensive.DWG is a kind of small winddriven generator that can access power distribution network.It is existing at present more to DG (distributions
Formula electric energy, distributed gerneration) and the DWG objects of planning research.
Kinematic nonlinearity planning problem is not known as multiple target, if hundreds and thousands of individual node accesses are distributed in planning region
Windpowered electricity generation, or in the practical application scenes such as distributed blower fan group polymorphic type, distributed power source variation, the optimal side of network arrangement
The determination of case can be highly difficult.In the Multiobjective programming models research of the power distribution network accessed to DWG, palpus consideration service reliability,
The systems organization target components such as running efficiency of system, cost of electricitygenerating, the quality of power supply, system safety.Included by the multiobjective Model
These target components contradiction and tradeoff relationship to each other be present.But traditional way does not account for some important key element mesh
Mark, such as the reliability in service and line security limit.Some are also not accounted for such as the practical situation of voltage rising.More
Document thinks that load is constant.But in fact, the load of system has voltagedependent, thus, it is supposed that load is constant
The result that analysis meeting must make mistake.So multiobjective problem is converted into singleobjective problem and draws a uncompromising solution by tradition
Certainly the optimization method of scheme is defective.
The content of the invention
The problem of existing present invention aim to address prior art, there is provided one kind is effectively based on Intrusion Index and balance
The distributed windpowered electricity generation multiobjective planning method of technology.
The technical scheme for realizing the object of the invention is that the distributed windpowered electricity generation multiple target based on Intrusion Index and balance technology is advised
The method of drawing,
Step 1：Try to achieve five influence indexs for not accessing DWG；Reliability index IR, active power loss coefficient ICL,
Electric power purchases coefficient correlation ICE, minimizes node voltage deviation effects index IVD, security implication index ISS；
Reliability index IR calculation formula is：
SAIDI is that system averagely has a power failure duration index, and SAIFI is system System average interruption frequency index, each for calculating
The lasting power failure cumulative time of user and frequency；
In formula, N_{i}It is the quantity for the user that ith powers off, N_{T}It is the total quantity of user, i and R are that ith is broken respectively
The recovery time of electricity and poweroff number；
w_{1}And w_{2}It is weight, respectively 0.4 and 0.6；Poweroff continuous weight shows reliably higher than poweroff frequency, IR value declines
Property improve.
The coefficient ICL calculation formula of active power loss is：
Active power loss cost containing DWG and without DWG is expressed as C_{LDWG}And C_{L}：
In formula, P_{LDWGj}And P_{Lj}It is the active power loss containing DWG and without DWG in jth of branch.λ_{t}Electricity when being t
Power price.Δ t is operation duration, and nbr is branch's number in network；ICL values be less than 1, show active power loss into
This reduction；ICL is equal to 1, shows that system is free of DWG.
Electric power buying coefficient correlation ICE calculation formula is：
Electric power purchase cost containing DWG and without DWG is respectively CE_{GDWG}And CE_{G}：
In formula, P_{Dit}The required load of node i when being t；λ_{t}Power price when being t；Δ t is operation duration；
P_{DWGi}It is the DWG of node i power；N is nodes；C_{lvDWG}It is the levelized power cost of DWG power supplies；
ICE numerical value is less than 1, shows that DWG operations are economical；When ICE be equal to 1, show that system is free of DWG；When ICE numerical value
More than 1, show that DWG operations are uneconomic.
The calculation formula for minimizing node voltage deviation IVD is：
System voltage containing DWG and without DWG is VDDWG, VD respectively compared with the average deviation of reference point：
In formula, V_{DWGi}And V_{i}It is the voltage of the node i containing DWG and without DWG respectively；V_{1}It is the voltage of reference mode；
IVD numerical value is less than 1, illustrates that node voltage is adjusted and is improved；IVD is equal to 1, illustrates that system is free of DWG.
The Intrusion Index ISS calculation formula of safety is：
MVA_{AV}And MVACAP_{AV}It is the average trend of circuit and the capacity of circuit respectively, is expressed as follows：
In formula, MVA_{j}And MVACAP_{j}It is circuit j actual trend and maximum capacity, nbr is branch's number；ISS numerical value reduces
Show the release of power system capacity.
Step 2：The multiple objective function of DWG planing methods is：
MinC (X)=min [IR (X), ICL (X), ICE (X), IVD (X), ISS (X)]^{T},
X is the vector of position and size in formula, and for determining the influence index of minimum, T is transposed matrix, related constraint bar
Part is as follows：
ε_{i}=ε_{i}*+Δε_{i}；
ε_{i}It is target i maximum allowable limiting value；ε_{i}* target i initial noninferior solution is represented, the compromise preference of designer is
Δε_{i}；
Equality constraint：Active and idle power protection must is fulfilled for the needs of network, including inputoutput power, owns
The DWG power of node, the algebraical sum of the line loss of whole power distribution network are necessarily equal to 0；
In formula, P_{i}And Q_{i}It is the active power and reactive power of the input of node i respectively；P_{Di}And Q_{Di}It is node i respectively
Active and reactive power load；P_{Lk}And Q_{Lk}It is branch road k active and reactive power loss；P_{DWGi}It is the DWG of node i wattful power
Rate, nbr are the branch road quantity of power distribution network；
Inequality constraints must be satisfied：
(a) the capacity MVA of circuit：The trend of any branch road does not allow more than the capacity of circuit in system：
Sij≤CSij
In formula, Sij and CSij are the actual trend capacity of connecting node i and j branches respectively；
(b) node voltage limits：It is electric containing meeting under DWG whether all loading condictions including offsetting voltage rise phenomenon
Press the constraints of limitation：
V0≤Vmax；
Vend ∣ max load, no DWG >=V min；
VDWG ∣ min load, max DWG≤V max；
VDWG ∣ min load, max DWG≤V0, i；
Vmin≤Vi≤Vmax；
In formula, V0 is the voltage of power supply point, and Vmin and Vmax are the minimum allowable value of node voltage respectively and maximum allowable
Value；Vi is the voltage of node i；Vend ∣ max load, no DWG are the feed voltages when load is maximum and is free of DWG；VDWG∣
Min load, max DWG are the voltage of DWG onpositions when minimum load and maximum DWG permeate；
Setting interaction counts r=1；
Step 3：One in five influence indexs of selection is used as major heading；
Step 4：The corresponding deviation of each target setting is compromisedIt is determined that the limitation of maximum or minimum allowable, by remaining
Target as constraints；
Step 5：Computing genetic algorithm calculates；Random generation DWG node location and size pair：
DWG_location=round (2+rand* (number_of_buses2))
DWG_size=(rand*maximum_DWG_size)
In formula, round is the function that rounds up, and non integer value tends to generate node location, rand during nearest integer value
The function for the Arbitrary Digit being randomly generated between 0~1；
Step 6：If result can compromise, into step 10；Otherwise, into next step；
Step 7：Interaction counts increase by 1, r=r+1；
Step 8：New deviation is set to compromiseAnd repeat step 5  step 8；
Step 9：If major heading needs to change, return to step four；
Step 10：Try to achieve most compromise solution.
After employing abovementioned technical proposal, the present invention has following positive effect：The present invention proposes distributed wind
The multiple objective function model of crucial pinpoint target in electricity planning, having considered 5 target components includes service reliability, system
Operational efficiency, cost of electricitygenerating, the quality of power supply and system safety, this method are solved more using interactive balance derivation algorithm
εconstrained procedure of objective optimisation problems.By should in this way, designer can make decisions on one's own compromise nondomination solution, nothing
By comprising load model whether, this method is effective.
Brief description of the drawings
In order that present disclosure is easier to be clearly understood, it is right below according to specific embodiment and with reference to accompanying drawing
The present invention is described in further detail, wherein
Fig. 1 is position containing RCS and the 33 meshed network figures for protecting fuse.
Embodiment
Test network is a 11KV, the power distribution network of 28 nodes, and related data is referred to shown in Tables 1 and 2.The base of the system
Quasi value is 1000KVA and 11KV.Consideration DWG range of capacity is 1.0884p.u. in 0~1p.u., system aggregate demand, and DWG is
With unity power factor operation.The levelized power cost of DWG power supplies includes the fund cost of installation, operation and maintenance expense
With run time and DWG life cycle etc., DWG rate for incorporation into the power network is about 0.51~0.58 yuan/kwh [19].Except relaxation
Outside node, each node can be used to configure the DWG in respective volume magnitude range in system.Outside 0.941.06p.u.
Node voltage be considered as collapse of voltage, if the node voltage that DWG is configured is higher than transformer substation voltage, what voltage rose shows
As will appear from.The capacity that the trend of any circuit exceedes the circuit will be considered as to violate circuit limitation.On any generation voltage
The DWG of phenomenon position and capacity is risen, violates voltage constraint, or violates circuit limitation, is all not distributed in power distribution network
Position.
As shown in figure 1, in addition to main circuit breaker performance DWG ability causes the reliability raising of system, two reclosings
The opening position in power distribution network acceptance of the bid arrow must be configured.To tackle failure, distribution network system is divided into 3 pieces of regions by these reclosings,
That is Z1, Z2 and Z3.The reclosing nearest from abort situation by Fault Isolation, must reduce the number of users of poweroff as far as possible.Originally grind
Study carefully and assume in emergency circumstances two reclosings while run work.If DWG is configured in nonfaulting region, and DWG disclosure satisfy that area
Domain load capacity requires that the user in region will not face poweroff.This method can reduce the System average interruption frequency index of system
(SAIFI) and averagely power failure duration index (SAIDI), so as to improve system reliability.The voltage constraint used in function can be protected
Demonstrate,proving any node does not have voltage outoflimit.Therefore, by using LTC transformers, transformer substation voltage is arranged to rated voltage
103%.
The system and load data of the node of table 1 28
SL=circuit maximum MVA (p.u)；PD=active power (p.u.)；Load power coefficient cos φ=0.7；It is idle negative
Lotus QL=tan φ PD；T=load types；I=commercial power；R=residential electricity consumptions；C=commercial powers.
The power failure frequency of the overhead line of table 2 and repair time
Two specific examples are carried out below to confirm that the method for the present invention is effective.
(example 1)：
Example 1 the results are shown in Table 4, since obtaining the initial scheme without DWG, the related Intrusion Index of reliability
(IR), the coefficient about power loss (ICL), electric power purchase cost index (ICE) are unit value 1.ICL and ICE unit values
Represent when DWG is not run, lose 1 times that cost and electric power purchase cost are corresponding costs after system optimization.IR unit value tables
Show that system reliability does not have any improvement, all users are likely to face power breakdown.Voltage relative influence index IVD is
0.0656, this represents that the average voltage deviations of node are 0.0656p.u., is now run without DWG.The system margin of safety influences to refer to
Number ISS is related to the marginal load of circuit, and when being run without DWG, ISS values are 0.9709, and this represents that load reaches capacity of trunk
97.09%, the only capacity of residue 2.91% loads available for other.In interaction 1, designer have selected loss cost correlation and refer to
ICL is marked as main target, and remaining target is as constraints.Balance preference is the maximum or most I compared with initial value
The deviation of receiving, this is designer's decision.This research is set as that the compromiseable preference of ICE target settings is 0.015p.u., this meaning
The reduction for ICE is no less than the 1.5% of initial value.And the electric power purchase cost containing DWG must not exceed reference scheme
98.5%.The compromiseable preference of other targets is arranged to 0, this explanation designer does not have any constraint to its these target.In mesh
In scalar functions, the limitation to voltage class and circuit sets constraint.The simulation result of interaction 1 is shown：DWG optimum position
It is node 7,0.48p.u. respectively with capacity.Distributing rationally for DWG reduces ICL indexs to 0.5530, shows to lose cost drop
Low 44.7%.Similarly, ICE is reduced to 0.9849, and the totle drilling cost of electric power buying reduces 1.51%.IVD from 0.0656 to
Reduce to 0.0361, it means that the average voltage deviations of node reduce 44.9%.ISS is reduced to 0.7902 by 0.9709,
It is 21%, rather than the 2.91% of reference scheme to show average line load margin.IR is reduced to 0.6271, illustrates 38.28%
Power supply trouble will not occur for user, therefore the reliability of system improves 38.28%.
The analog result of the example 1 of table 3
If designer's plan further reduces loss cost (ICL), and it is main target to select ICL again, then may be used
Capable scheme is to increase DWG infiltration.But the infiltration for improving DWG can increase electric power purchase cost, because electric power purchase cost
Include DWG power supply cost, this may be above power network power supply cost.Based on such a situation, designer selects in interaction 2
ICE compromiseable preference is 0.01p.u., it is meant that ICE reduces obtains 1% no less than initial value, and the limiting value of interaction 1 is
0.015p.u..Other constraints will be kept similar with interaction 1 in interaction 2.The solution of interaction 2, it is contemplated that new equilibrium valve,
Show that ICL is further decreased to 0.5419 from the 0.5530 of interaction 1.ICE rises to 0.9897 from the 0.9849 of interaction 1.IVD drops
As little as 0.0281, show the average voltage deviations of node has further reduction from the 0.0361 of interaction 1.ISS is from interaction 1
0.7902 increases to 0.7914.Compared with interaction 1, reliability index IR does not have significant change, because region has same time
Several disconnection faults occurs.In interaction 2, DWG position is still node 7, capacity from original 0.48p.u. improve to
0.62p.u.。
In the method, the selectable authority of designer, setting preference value is required according to optimization, other targets can be selected
Implemented as main target, and with this.It is assumed that designer is satisfied to the result of interaction 2, it is positioned at the capacity of node 7
0.62p.u. DWG allocation plan is optimal, because ICL is the minimum value of position and capacity logarithm.The shadow of some nodes
Snap number and DWG size are shown as 0, after this explanation DWG accesses the position of these nodes, it may appear that the phenomenon that voltage rises,
Either violate voltage or circuit constraint or be unsatisfactory for the status condition that designer determines in the interaction stage.Therefore, these are not
It is the position for being adapted to DWG accesses.In this case, node 2,3,4,11,12,13,14,15,16,21, and 28 be not to match somebody with somebody
Put DWG feasible location.
(example 2)
To make simulation test result more accurate, load model is included in interaction 2, its main target function and constraint bar
Part is consistent with example 1.House, industry, business, and four kinds of application scenarios of combined grade, respective analog result are considered respectively
As shown in table 4.
The analog result of the example 2 of table 4
As can be seen that load type does not influence DG position, for example node 7 can access the DWG of different load, and DWG
Capacity is changed because of the change of load.
Particular embodiments described above, the purpose of the present invention, technical scheme and beneficial effect are carried out further in detail
Describe in detail it is bright, should be understood that the foregoing is only the present invention specific embodiment, be not intended to limit the invention, it is all
Within the spirit and principles in the present invention, any modification, equivalent substitution and improvements done etc., it should be included in the guarantor of the present invention
Within the scope of shield.
Claims (6)
1. the distributed windpowered electricity generation multiobjective planning method based on Intrusion Index and balance technology, it is characterised in that including following step
Suddenly：
Step 1：Try to achieve five influence indexs for not accessing DWG；Reliability index IR, active power loss coefficient ICL, electric power
Coefficient correlation ICE is purchased, minimizes node voltage deviation effects index IVD, security implication index ISS；
Step 2：The multiple objective function of DWG planing methods is：
Min C (X)=min [IR (X), ICL (X), ICE (X), IVD (X), ISS (X)]^{T},
X is the vector of position and size in formula, and for determining the influence index of minimum, T is transposed matrix, and relevant constraint is such as
Under：
ε_{i}=ε_{i}*+Δε_{i}；
ε_{i}It is target i maximum allowable limiting value；ε_{i}* target i initial noninferior solution is represented, the compromise preference of designer is Δ ε_{i}；
Equality constraint：Active and idle power protection must is fulfilled for the needs of network, including inputoutput power, all nodes
DWG power, the algebraical sum of the line loss of whole power distribution network is necessarily equal to 0；
<mrow>
<msubsup>
<mo>&Sigma;</mo>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>n</mi>
</msubsup>
<msub>
<mi>P</mi>
<mi>i</mi>
</msub>
<mo></mo>
<msubsup>
<mo>&Sigma;</mo>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>2</mn>
</mrow>
<mi>n</mi>
</msubsup>
<mrow>
<mo>(</mo>
<msub>
<mi>P</mi>
<mrow>
<mi>D</mi>
<mi>i</mi>
</mrow>
</msub>
<mo></mo>
<msub>
<mi>P</mi>
<mrow>
<mi>D</mi>
<mi>W</mi>
<mi>G</mi>
<mi>i</mi>
</mrow>
</msub>
<mo>)</mo>
</mrow>
<mo></mo>
<msubsup>
<mo>&Sigma;</mo>
<mrow>
<mi>k</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mrow>
<mi>n</mi>
<mi>b</mi>
<mi>r</mi>
</mrow>
</msubsup>
<msub>
<mi>P</mi>
<mrow>
<mi>L</mi>
<mi>k</mi>
</mrow>
</msub>
<mo>=</mo>
<mn>0</mn>
<mo>;</mo>
</mrow>
<mrow>
<msubsup>
<mo>&Sigma;</mo>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>n</mi>
</msubsup>
<msub>
<mi>Q</mi>
<mi>i</mi>
</msub>
<mo></mo>
<msubsup>
<mo>&Sigma;</mo>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>2</mn>
</mrow>
<mi>n</mi>
</msubsup>
<msub>
<mi>Q</mi>
<mrow>
<mi>D</mi>
<mi>i</mi>
</mrow>
</msub>
<mo></mo>
<msubsup>
<mo>&Sigma;</mo>
<mrow>
<mi>k</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mrow>
<mi>n</mi>
<mi>b</mi>
<mi>r</mi>
</mrow>
</msubsup>
<msub>
<mi>Q</mi>
<mrow>
<mi>L</mi>
<mi>k</mi>
</mrow>
</msub>
<mo>=</mo>
<mn>0</mn>
<mo>;</mo>
</mrow>
In formula, P_{i}And Q_{i}It is the active power and reactive power of the input of node i respectively；P_{Di}And Q_{Di}It is the active of node i respectively
With reactive power load；P_{Lk}And Q_{Lk}It is branch road k active and reactive power loss；P_{DWGi}It is the DWG of node i active power,
Nbr is the branch road quantity of power distribution network；
Inequality constraints must be satisfied：
(a) the capacity MVA of circuit：The trend of any branch road does not allow more than the capacity of circuit in system：
Sij≤CSij
In formula, Sij and CSij are the actual trend capacity of connecting node i and j branches respectively；
(b) node voltage limits：Including offset voltage rise phenomenon, containing meet under DWG whether all loading condictions voltage limit
The constraints of system：
V0≤Vmax；
Vend ∣ max load, no DWG >=Vmin；
VDWG ∣ min load, max DWG≤Vmax；
VDWG ∣ min load, max DWG≤V0, i；
Vmin≤Vi≤Vmax；
In formula, V0 is the voltage of power supply point, and Vmin and Vmax are the minimum allowable value and maximum permissible value of node voltage respectively；Vi
It is the voltage of node i；Vend ∣ max load, no DWG are the feed voltages when load is maximum and is free of DWG；VDWG∣min
Load, max DWG are the voltage of DWG onpositions when minimum load and maximum DWG permeate；
Setting interaction counts r=1；
Step 3：One in five influence indexs of selection is used as major heading；
Step 4：The corresponding deviation of each target setting is compromisedIt is determined that the limitation of maximum or minimum allowable, by remaining target
As constraints；
Step 5：Computing genetic algorithm calculates；Random generation DWG node location and size pair：
DWG_location=round (2+rand* (number_of_buses2))
DWG_size=(rand*maximum_DWG_size)
In formula, round is the function that rounds up, and non integer value tends to generate node location during nearest integer value, rand be with
The function of Arbitrary Digit between machine generation 0~1；
Step 6：If result can compromise, into step 10；Otherwise, into next step；
Step 7：Interaction counts increase by 1, r=r+1；
Step 8：New deviation is set to compromiseAnd repeat step 5  step 8；
Step 9：If major heading needs to change, return to step four；
Step 10：Try to achieve most compromise solution.
2. the distributed windpowered electricity generation multiobjective planning method according to claim 1 based on Intrusion Index and balance technology, its
It is characterised by：
In the step 1, reliability index IR calculation formula is：
<mrow>
<mi>I</mi>
<mi>R</mi>
<mo>=</mo>
<msub>
<mi>W</mi>
<mn>1</mn>
</msub>
<mfrac>
<mrow>
<mi>S</mi>
<mi>A</mi>
<mi>I</mi>
<mi>F</mi>
<mi>I</mi>
</mrow>
<mrow>
<msub>
<mi>SAIFI</mi>
<mrow>
<mi>D</mi>
<mi>G</mi>
</mrow>
</msub>
</mrow>
</mfrac>
<mo>+</mo>
<msub>
<mi>W</mi>
<mn>2</mn>
</msub>
<mfrac>
<mrow>
<mi>S</mi>
<mi>A</mi>
<mi>I</mi>
<mi>F</mi>
<mi>I</mi>
</mrow>
<mrow>
<msub>
<mi>SAIFI</mi>
<mrow>
<mi>D</mi>
<mi>G</mi>
</mrow>
</msub>
</mrow>
</mfrac>
<mo>;</mo>
</mrow>
SAIDI is that system averagely has a power failure duration index, and SAIFI is system System average interruption frequency index, for calculating each user
The lasting power failure cumulative time and frequency；
<mrow>
<mi>S</mi>
<mi>A</mi>
<mi>I</mi>
<mi>F</mi>
<mi>I</mi>
<mo>=</mo>
<mfrac>
<mrow>
<msub>
<mo>&Sigma;</mo>
<mrow>
<mi>i</mi>
<mo>&Element;</mo>
<mi>R</mi>
</mrow>
</msub>
<msub>
<mi>N</mi>
<mi>i</mi>
</msub>
</mrow>
<msub>
<mi>N</mi>
<mi>T</mi>
</msub>
</mfrac>
</mrow>
<mrow>
<mi>S</mi>
<mi>A</mi>
<mi>I</mi>
<mi>D</mi>
<mi>I</mi>
<mo>=</mo>
<mfrac>
<mrow>
<msub>
<mo>&Sigma;</mo>
<mrow>
<mi>i</mi>
<mo>&Element;</mo>
<mi>R</mi>
</mrow>
</msub>
<msub>
<mi>&mu;</mi>
<mi>i</mi>
</msub>
<msub>
<mi>N</mi>
<mi>i</mi>
</msub>
</mrow>
<msub>
<mi>N</mi>
<mi>T</mi>
</msub>
</mfrac>
</mrow>
In formula, N_{i}It is the quantity for the user that ith powers off, N_{T}It is the total quantity of user, i and R are that ith powers off respectively
Recovery time and poweroff number；
W_{1}And W_{2}It is weight, respectively 0.4 and 0.6；Poweroff continuous weight shows that reliability carries higher than poweroff frequency, IR value declines
It is high.
3. the distributed windpowered electricity generation multiobjective planning method according to claim 1 based on Intrusion Index and balance technology, its
It is characterised by：
In the step 1, the coefficient ICL of active power loss calculation formula is：
<mrow>
<mi>I</mi>
<mi>C</mi>
<mi>L</mi>
<mo>=</mo>
<mfrac>
<msub>
<mi>C</mi>
<mrow>
<mi>L</mi>
<mi>D</mi>
<mi>W</mi>
<mi>G</mi>
</mrow>
</msub>
<msub>
<mi>C</mi>
<mi>L</mi>
</msub>
</mfrac>
<mo>;</mo>
</mrow>
Active power loss cost containing DWG and without DWG is expressed as C_{LDWG}And C_{L}：
<mrow>
<msub>
<mi>C</mi>
<mi>L</mi>
</msub>
<mo>=</mo>
<msubsup>
<mo>&Sigma;</mo>
<mrow>
<mi>j</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mrow>
<mi>n</mi>
<mi>b</mi>
<mi>r</mi>
</mrow>
</msubsup>
<msub>
<mi>P</mi>
<mrow>
<mi>L</mi>
<mi>j</mi>
</mrow>
</msub>
<msub>
<mi>&lambda;</mi>
<mi>t</mi>
</msub>
<mi>&Delta;</mi>
<mi>t</mi>
</mrow>
<mrow>
<msub>
<mi>C</mi>
<mrow>
<mi>L</mi>
<mi>D</mi>
<mi>E</mi>
<mi>G</mi>
</mrow>
</msub>
<mo>=</mo>
<msubsup>
<mo>&Sigma;</mo>
<mrow>
<mi>j</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mrow>
<mi>n</mi>
<mi>b</mi>
<mi>r</mi>
</mrow>
</msubsup>
<msub>
<mi>P</mi>
<mrow>
<mi>L</mi>
<mi>D</mi>
<mi>W</mi>
<mi>E</mi>
<mi>G</mi>
<mi>j</mi>
</mrow>
</msub>
<msub>
<mi>&lambda;</mi>
<mi>t</mi>
</msub>
<mi>&Delta;</mi>
<mi>t</mi>
</mrow>
In formula, P_{LDWGj}And P_{Lj}It is the active power loss containing DWG and without DWG in jth of branch；λ_{t}Electric power valency when being t
Lattice；Δ t is operation duration, and nbr is branch's number in network；ICL values are less than 1, show active power loss cost drop
It is low；ICL is equal to 1, shows that system is free of DWG.
4. the distributed windpowered electricity generation multiobjective planning method according to claim 1 based on Intrusion Index and balance technology, its
It is characterised by：
In the step 1, electric power buying coefficient correlation ICE calculation formula is：
<mrow>
<mi>I</mi>
<mi>C</mi>
<mi>E</mi>
<mo>=</mo>
<mfrac>
<mrow>
<msub>
<mi>CE</mi>
<mrow>
<mi>G</mi>
<mi>D</mi>
<mi>W</mi>
<mi>G</mi>
</mrow>
</msub>
</mrow>
<mrow>
<msub>
<mi>CE</mi>
<mi>G</mi>
</msub>
</mrow>
</mfrac>
</mrow>
Electric power purchase cost containing DWG and without DWG is respectively CE_{GDWG}And CE_{G}：
<mrow>
<msub>
<mi>CE</mi>
<mi>G</mi>
</msub>
<mo>=</mo>
<msubsup>
<mi>&Sigma;</mi>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>n</mi>
</msubsup>
<msub>
<mi>&lambda;</mi>
<mi>t</mi>
</msub>
<msub>
<mi>P</mi>
<mrow>
<mi>D</mi>
<mi>i</mi>
<mi>t</mi>
</mrow>
</msub>
<mi>&Delta;</mi>
<mi>t</mi>
</mrow>
<mrow>
<msub>
<mi>CE</mi>
<mrow>
<mi>G</mi>
<mi>D</mi>
<mi>W</mi>
<mi>G</mi>
</mrow>
</msub>
<mo>=</mo>
<msubsup>
<mi>&Sigma;</mi>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>n</mi>
</msubsup>
<msub>
<mi>&lambda;</mi>
<mi>t</mi>
</msub>
<mrow>
<mo>(</mo>
<msub>
<mi>P</mi>
<mrow>
<mi>D</mi>
<mi>i</mi>
<mi>t</mi>
</mrow>
</msub>
<mo></mo>
<msub>
<mi>P</mi>
<mrow>
<mi>D</mi>
<mi>W</mi>
<mi>G</mi>
<mi>i</mi>
</mrow>
</msub>
<mo>)</mo>
</mrow>
<mi>&Delta;</mi>
<mi>t</mi>
<mo>+</mo>
<msub>
<mi>C</mi>
<mrow>
<mi>l</mi>
<mi>v</mi>
<mi>D</mi>
<mi>W</mi>
<mi>G</mi>
</mrow>
</msub>
<msub>
<mi>P</mi>
<mrow>
<mi>D</mi>
<mi>W</mi>
<mi>G</mi>
<mi>i</mi>
</mrow>
</msub>
<mi>&Delta;</mi>
<mi>t</mi>
</mrow>
In formula, P_{Dit}The required load of node i when being t；λ_{t}Power price when being t；Δ t is operation duration；P_{DWGi}It is
The DWG of node i power；N is nodes；C_{lvDWG}It is the levelized power cost of DWG power supplies；
ICE numerical value is less than 1, shows that DWG operations are economical；When ICE be equal to 1, show that system is free of DWG；When ICE numerical value is more than
1, show that DWG operations are uneconomic.
5. the distributed windpowered electricity generation multiobjective planning method according to claim 1 based on Intrusion Index and balance technology, its
It is characterised by：
In the step 1, the calculation formula for minimizing node voltage deviation IVD is：
<mrow>
<mi>I</mi>
<mi>V</mi>
<mi>D</mi>
<mo>=</mo>
<mfrac>
<msub>
<mi>v</mi>
<mrow>
<mi>D</mi>
<mi>D</mi>
<mi>W</mi>
<mi>G</mi>
</mrow>
</msub>
<msub>
<mi>v</mi>
<mi>D</mi>
</msub>
</mfrac>
<mo>;</mo>
</mrow>
System voltage containing DWG and without DWG is VDDWG, VD respectively compared with the average deviation of reference point：
<mrow>
<mi>V</mi>
<mi>D</mi>
<mo>=</mo>
<mfrac>
<mn>1</mn>
<mrow>
<mi>n</mi>
<mo></mo>
<mn>1</mn>
</mrow>
</mfrac>
<msubsup>
<mo>&Sigma;</mo>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>2</mn>
</mrow>
<mi>n</mi>
</msubsup>
<mo></mo>
<mrow>
<mo>(</mo>
<msub>
<mi>V</mi>
<mi>i</mi>
</msub>
<mo></mo>
<msub>
<mi>V</mi>
<mn>1</mn>
</msub>
<mo>)</mo>
</mrow>
<mo></mo>
</mrow>
<mrow>
<mi>V</mi>
<mi>D</mi>
<mi>D</mi>
<mi>W</mi>
<mi>G</mi>
<mo>=</mo>
<mfrac>
<mn>1</mn>
<mrow>
<mi>n</mi>
<mo></mo>
<mn>1</mn>
</mrow>
</mfrac>
<msubsup>
<mo>&Sigma;</mo>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>2</mn>
</mrow>
<mi>n</mi>
</msubsup>
<mo></mo>
<mrow>
<mo>(</mo>
<msub>
<mi>V</mi>
<mrow>
<mi>D</mi>
<mi>W</mi>
<mi>G</mi>
<mi>i</mi>
</mrow>
</msub>
<mo></mo>
<msub>
<mi>V</mi>
<mn>1</mn>
</msub>
<mo>)</mo>
</mrow>
<mo></mo>
</mrow>
In formula, V_{DWGi}And V_{i}It is the voltage of the node i containing DWG and without DWG respectively；V_{1}It is the voltage of reference mode；
IVD numerical value is less than 1, illustrates that node voltage is adjusted and is improved；IVD is equal to 1, illustrates that system is free of DWG.
6. the distributed windpowered electricity generation multiobjective planning method according to claim 1 based on Intrusion Index and balance technology, its
It is characterised by：
In the step 1, safe Intrusion Index ISS calculation formula is：
<mrow>
<mi>I</mi>
<mi>S</mi>
<mi>S</mi>
<mo>=</mo>
<mfrac>
<mrow>
<msub>
<mi>MVA</mi>
<mrow>
<mi>A</mi>
<mi>V</mi>
</mrow>
</msub>
</mrow>
<mrow>
<msub>
<mi>MVACAP</mi>
<mrow>
<mi>A</mi>
<mi>V</mi>
</mrow>
</msub>
</mrow>
</mfrac>
<mo>;</mo>
</mrow>
MVA_{AV}And MVACAP_{AV}It is the average trend of circuit and the capacity of circuit respectively, is expressed as follows：
<mrow>
<msub>
<mi>MVA</mi>
<mrow>
<mi>A</mi>
<mi>V</mi>
</mrow>
</msub>
<mo>=</mo>
<mfrac>
<mn>1</mn>
<mrow>
<mi>n</mi>
<mi>b</mi>
<mi>r</mi>
</mrow>
</mfrac>
<msubsup>
<mo>&Sigma;</mo>
<mrow>
<mi>j</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mrow>
<mi>n</mi>
<mi>b</mi>
<mi>r</mi>
</mrow>
</msubsup>
<msub>
<mi>MVA</mi>
<mi>j</mi>
</msub>
</mrow>
<mrow>
<msub>
<mi>MVACAP</mi>
<mrow>
<mi>A</mi>
<mi>V</mi>
</mrow>
</msub>
<mo>=</mo>
<mfrac>
<mn>1</mn>
<mrow>
<mi>n</mi>
<mi>b</mi>
<mi>r</mi>
</mrow>
</mfrac>
<msubsup>
<mo>&Sigma;</mo>
<mrow>
<mi>j</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mrow>
<mi>n</mi>
<mi>b</mi>
<mi>r</mi>
</mrow>
</msubsup>
<msub>
<mi>MVACAP</mi>
<mi>j</mi>
</msub>
</mrow>
In formula, MVA_{j}And MVACAP_{j}It is circuit j actual trend and maximum capacity, nbr is branch's number；ISS numerical value, which reduces, to be shown
The release of power system capacity.
Priority Applications (1)
Application Number  Priority Date  Filing Date  Title 

CN201510609613.XA CN105226703B (en)  20150922  20150922  Distributed windpowered electricity generation multiobjective planning method based on Intrusion Index and balance technology 
Applications Claiming Priority (1)
Application Number  Priority Date  Filing Date  Title 

CN201510609613.XA CN105226703B (en)  20150922  20150922  Distributed windpowered electricity generation multiobjective planning method based on Intrusion Index and balance technology 
Publications (2)
Publication Number  Publication Date 

CN105226703A CN105226703A (en)  20160106 
CN105226703B true CN105226703B (en)  20171110 
Family
ID=54995493
Family Applications (1)
Application Number  Title  Priority Date  Filing Date 

CN201510609613.XA CN105226703B (en)  20150922  20150922  Distributed windpowered electricity generation multiobjective planning method based on Intrusion Index and balance technology 
Country Status (1)
Country  Link 

CN (1)  CN105226703B (en) 
Families Citing this family (3)
Publication number  Priority date  Publication date  Assignee  Title 

CN107392350A (en) *  20170608  20171124  国网宁夏电力公司电力科学研究院  Power distribution network Expansion Planning comprehensive optimization method containing distributed energy and charging station 
JP2019105973A (en) *  20171212  20190627  株式会社日立製作所  Maintenance plan generation system 
CN108988370A (en) *  20180822  20181211  广东电网有限责任公司  The capacity determining methods of energy storage device, equipment and storage medium in electric system 
Family Cites Families (5)
Publication number  Priority date  Publication date  Assignee  Title 

CN102623989B (en) *  20120328  20131127  湖南大学  Method for optimization and configuration of intermittent distributed generation (DG) 
CN103353979B (en) *  20130531  20161228  国家电网公司  The Optimizing Site Selection of a kind of distributed power source and constant volume method 
CN103903073B (en) *  20140423  20170308  河海大学  A kind of microcapacitance sensor Method for optimized planning containing distributed power source and energy storage and system 
CN104852373B (en) *  20150429  20180316  南方电网科学研究院有限责任公司  A kind of dynamic economic dispatch optimization method of the microcapacitance sensor of independence containing intermittent energy source 
CN104881716A (en) *  20150528  20150902  贵州电网公司电网规划研究中心  Optimization programming and evaluation method of microgrid power supply 

2015
 20150922 CN CN201510609613.XA patent/CN105226703B/en not_active IP Right Cessation
Also Published As
Publication number  Publication date 

CN105226703A (en)  20160106 
Similar Documents
Publication  Publication Date  Title 

Wang et al.  Coordinated control of distributed energystorage systems for voltage regulation in distribution networks  
Biggar et al.  The economics of electricity markets  
Capitanescu et al.  Assessing the potential of network reconfiguration to improve distributed generation hosting capacity in active distribution systems  
Collins et al.  Real and reactive power control of distributed PV inverters for overvoltage prevention and increased renewable generation hosting capacity  
Ding et al.  On distributed PV hosting capacity estimation, sensitivity study, and improvement  
Long et al.  Voltage control of PVrich LV networks: OLTCfitted transformer and capacitor banks  
Zhao et al.  A review of active management for distribution networks: current status and future development trends  
Mohammadi et al.  Challenges of PV integration in lowvoltage secondary networks  
Chen et al.  Analysis of voltage profile problems due to the penetration of distributed generation in lowvoltage secondary distribution networks  
PadilhaFeltrin et al.  VoltVAR multiobjective optimization to peakload relief and energy efficiency in distribution networks  
Jayasekara et al.  Optimal operation of distributed energy storage systems to improve distribution network load and generation hosting capability  
Abdi et al.  Application of IPSOMonte Carlo for optimal distributed generation allocation and sizing  
AbdelGhany et al.  Optimizing DG penetration in distribution networks concerning protection schemes and technical impact  
Shahnia et al.  Voltage imbalance analysis in residential low voltage distribution networks with rooftop PVs  
KavousiFard et al.  Optimal distribution feeder reconfiguration for reliability improvement considering uncertainty  
Abapour et al.  Dynamic planning of distributed generation units in active distribution network  
Kumar et al.  DG integrated approach for service restoration under cold load pickup  
Wang et al.  Integrating electrical energy storage into coordinated voltage control schemes for distribution networks  
Khan et al.  Implementation of Distributed Generation (IDG) algorithm for performance enhancement of distribution feeder under extreme load growth  
Sakar et al.  Increasing PV hosting capacity in distorted distribution systems using passive harmonic filtering  
Liu et al.  Probabilistic estimation of voltage unbalance in MV distribution networks with unbalanced load  
Fan et al.  Probabilistic power flow analysis with generation dispatch including photovoltaic resources  
Deilami et al.  Voltage profile and THD distortion of residential network with high penetration of Plugin Electrical Vehicles  
Dondi et al.  Network integration of distributed power generation  
Biswas et al.  Optimum distributed generation placement with voltage sag effect minimization 
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  
CF01  Termination of patent right due to nonpayment of annual fee  
CF01  Termination of patent right due to nonpayment of annual fee 
Granted publication date: 20171110 Termination date: 20180922 