CN105226703B - Distributed wind-powered electricity generation multi-objective planning method based on Intrusion Index and balance technology - Google Patents

Distributed wind-powered electricity generation multi-objective planning method based on Intrusion Index and balance technology Download PDF

Info

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
Application number
CN201510609613.XA
Other languages
Chinese (zh)
Other versions
CN105226703A (en
Inventor
蔡浩
朱熀秋
唐静
Original Assignee
江苏大学
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by 江苏大学 filed Critical 江苏大学
Priority to CN201510609613.XA priority Critical patent/CN105226703B/en
Publication of CN105226703A publication Critical patent/CN105226703A/en
Application granted granted Critical
Publication of CN105226703B publication Critical patent/CN105226703B/en

Links

Classifications

    • YGENERAL 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
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02EREDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
    • Y02E10/00Energy generation through renewable energy sources
    • Y02E10/70Wind energy
    • Y02E10/76Power conversion electric or electronic aspects

Abstract

The invention discloses the distributed wind-powered electricity generation multi-objective 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 non-domination solution.

Description

Distributed wind-powered electricity generation multi-objective planning method based on Intrusion Index and balance technology
Technical field
The present invention relates to a kind of distributed wind-powered electricity generation multi-objective 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 wind-power 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 wind-driven 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 Wind-powered 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 electricity-generating, the quality of power supply, system safety.Included by the multi-objective Model These target components contradiction and trade-off 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 voltage-dependent, thus, it is supposed that load is constant The result that analysis meeting must make mistake.So multi-objective problem is converted into single-objective 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 wind-powered electricity generation multi-objective planning method of technology.
The technical scheme for realizing the object of the invention is that the distributed wind-powered 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, NiIt is the quantity for the user that ith powers off, NTIt is the total quantity of user, i and R are that ith is broken respectively The recovery time of electricity and power-off number;
w1And w2It is weight, respectively 0.4 and 0.6;Power-off continuous weight shows reliably higher than power-off 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 CLDWGAnd CL
In formula, PLDWGjAnd PLjIt is the active power loss containing DWG and without DWG in j-th of branch.λtElectricity 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 CEGDWGAnd CEG
In formula, PDitThe required load of node i when being t;λtPower price when being t;Δ t is operation duration; PDWGiIt is the DWG of node i power;N is nodes;ClvDWGIt 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, VDWGiAnd ViIt is the voltage of the node i containing DWG and without DWG respectively;V1It 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:
MVAAVAnd MVACAPAVIt is the average trend of circuit and the capacity of circuit respectively, is expressed as follows:
In formula, MVAjAnd MVACAPjIt 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:
εii*+Δεi
εiIt 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 input-output 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, PiAnd QiIt is the active power and reactive power of the input of node i respectively;PDiAnd QDiIt is node i respectively Active and reactive power load;PLkAnd QLkIt is branch road k active and reactive power loss;PDWGiIt 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 on-positions 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_buses-2))
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 above-mentioned 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 electricity-generating, 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 non-domination 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.94-1.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 power-off as far as possible.Originally grind Study carefully and assume in emergency circumstances two reclosings while run work.If DWG is configured in non-faulting region, and DWG disclosure satisfy that area Domain load capacity requires that the user in region will not face power-off.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 out-of-limit.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 wind-powered electricity generation multi-objective 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:
εii*+Δεi
εiIt 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 input-output 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>&amp;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>&amp;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>&amp;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>&amp;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>&amp;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>&amp;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, PiAnd QiIt is the active power and reactive power of the input of node i respectively;PDiAnd QDiIt is the active of node i respectively With reactive power load;PLkAnd QLkIt is branch road k active and reactive power loss;PDWGiIt 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 on-positions 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_buses-2))
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 wind-powered electricity generation multi-objective 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>&amp;Sigma;</mo> <mrow> <mi>i</mi> <mo>&amp;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>&amp;Sigma;</mo> <mrow> <mi>i</mi> <mo>&amp;Element;</mo> <mi>R</mi> </mrow> </msub> <msub> <mi>&amp;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, NiIt is the quantity for the user that ith powers off, NTIt is the total quantity of user, i and R are that ith powers off respectively Recovery time and power-off number;
W1And W2It is weight, respectively 0.4 and 0.6;Power-off continuous weight shows that reliability carries higher than power-off frequency, IR value declines It is high.
3. the distributed wind-powered electricity generation multi-objective 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 CLDWGAnd CL
<mrow> <msub> <mi>C</mi> <mi>L</mi> </msub> <mo>=</mo> <msubsup> <mo>&amp;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>&amp;lambda;</mi> <mi>t</mi> </msub> <mi>&amp;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>&amp;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>&amp;lambda;</mi> <mi>t</mi> </msub> <mi>&amp;Delta;</mi> <mi>t</mi> </mrow>
In formula, PLDWGjAnd PLjIt is the active power loss containing DWG and without DWG in j-th of branch;λtElectric 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 wind-powered electricity generation multi-objective 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 CEGDWGAnd CEG
<mrow> <msub> <mi>CE</mi> <mi>G</mi> </msub> <mo>=</mo> <msubsup> <mi>&amp;Sigma;</mi> <mrow> <mi>i</mi> <mo>=</mo> <mn>1</mn> </mrow> <mi>n</mi> </msubsup> <msub> <mi>&amp;lambda;</mi> <mi>t</mi> </msub> <msub> <mi>P</mi> <mrow> <mi>D</mi> <mi>i</mi> <mi>t</mi> </mrow> </msub> <mi>&amp;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>&amp;Sigma;</mi> <mrow> <mi>i</mi> <mo>=</mo> <mn>1</mn> </mrow> <mi>n</mi> </msubsup> <msub> <mi>&amp;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>&amp;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>&amp;Delta;</mi> <mi>t</mi> </mrow>
In formula, PDitThe required load of node i when being t;λtPower price when being t;Δ t is operation duration;PDWGiIt is The DWG of node i power;N is nodes;ClvDWGIt 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 wind-powered electricity generation multi-objective 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>&amp;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>&amp;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, VDWGiAnd ViIt is the voltage of the node i containing DWG and without DWG respectively;V1It 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 wind-powered electricity generation multi-objective 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>
MVAAVAnd MVACAPAVIt 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>&amp;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>&amp;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, MVAjAnd MVACAPjIt 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.
CN201510609613.XA 2015-09-22 2015-09-22 Distributed wind-powered electricity generation multi-objective planning method based on Intrusion Index and balance technology CN105226703B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201510609613.XA CN105226703B (en) 2015-09-22 2015-09-22 Distributed wind-powered electricity generation multi-objective planning method based on Intrusion Index and balance technology

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201510609613.XA CN105226703B (en) 2015-09-22 2015-09-22 Distributed wind-powered electricity generation multi-objective planning method based on Intrusion Index and balance technology

Publications (2)

Publication Number Publication Date
CN105226703A CN105226703A (en) 2016-01-06
CN105226703B true CN105226703B (en) 2017-11-10

Family

ID=54995493

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201510609613.XA CN105226703B (en) 2015-09-22 2015-09-22 Distributed wind-powered electricity generation multi-objective planning method based on Intrusion Index and balance technology

Country Status (1)

Country Link
CN (1) CN105226703B (en)

Families Citing this family (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107392350A (en) * 2017-06-08 2017-11-24 国网宁夏电力公司电力科学研究院 Power distribution network Expansion Planning comprehensive optimization method containing distributed energy and charging station
JP2019105973A (en) * 2017-12-12 2019-06-27 株式会社日立製作所 Maintenance plan generation system
CN108988370A (en) * 2018-08-22 2018-12-11 广东电网有限责任公司 The capacity determining methods of energy storage device, equipment and storage medium in electric system

Family Cites Families (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102623989B (en) * 2012-03-28 2013-11-27 湖南大学 Method for optimization and configuration of intermittent distributed generation (DG)
CN103353979B (en) * 2013-05-31 2016-12-28 国家电网公司 The Optimizing Site Selection of a kind of distributed power source and constant volume method
CN103903073B (en) * 2014-04-23 2017-03-08 河海大学 A kind of micro-capacitance sensor Method for optimized planning containing distributed power source and energy storage and system
CN104852373B (en) * 2015-04-29 2018-03-16 南方电网科学研究院有限责任公司 A kind of dynamic economic dispatch optimization method of the micro-capacitance sensor of independence containing intermittent energy source
CN104881716A (en) * 2015-05-28 2015-09-02 贵州电网公司电网规划研究中心 Optimization programming and evaluation method of micro-grid power supply

Also Published As

Publication number Publication date
CN105226703A (en) 2016-01-06

Similar Documents

Publication Publication Date Title
Wang et al. Coordinated control of distributed energy-storage 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 PV-rich LV networks: OLTC-fitted 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 low-voltage secondary networks
Chen et al. Analysis of voltage profile problems due to the penetration of distributed generation in low-voltage secondary distribution networks
Padilha-Feltrin et al. Volt-VAR multiobjective optimization to peak-load 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 IPSO-Monte Carlo for optimal distributed generation allocation and sizing
Abdel-Ghany 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
Kavousi-Fard 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 Plug-in 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 non-payment of annual fee
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20171110

Termination date: 20180922