CN109066808A - A kind of uncertain active distribution network running optimizatin method of adaptation power supply power output - Google Patents

A kind of uncertain active distribution network running optimizatin method of adaptation power supply power output Download PDF

Info

Publication number
CN109066808A
CN109066808A CN201810979860.2A CN201810979860A CN109066808A CN 109066808 A CN109066808 A CN 109066808A CN 201810979860 A CN201810979860 A CN 201810979860A CN 109066808 A CN109066808 A CN 109066808A
Authority
CN
China
Prior art keywords
scene
power
under
active
node
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
CN201810979860.2A
Other languages
Chinese (zh)
Other versions
CN109066808B (en
Inventor
刘志坚
陈潇雅
边居政
王洪亮
罗凌灵
宋琪
谢静
刘亚锦
王雁红
周于尧
徐慧
王妃
王一妃
余进
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Kunming University of Science and Technology
Original Assignee
Kunming University of Science and Technology
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 Kunming University of Science and Technology filed Critical Kunming University of Science and Technology
Priority to CN201810979860.2A priority Critical patent/CN109066808B/en
Publication of CN109066808A publication Critical patent/CN109066808A/en
Application granted granted Critical
Publication of CN109066808B publication Critical patent/CN109066808B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • H—ELECTRICITY
    • H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
    • H02J3/00—Circuit arrangements for ac mains or ac distribution networks
    • H02J3/38—Arrangements for parallely feeding a single network by two or more generators, converters or transformers
    • H02J3/46—Controlling of the sharing of output between the generators, converters, or transformers
    • H—ELECTRICITY
    • H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
    • H02J3/00—Circuit arrangements for ac mains or ac distribution networks
    • H02J3/38—Arrangements for parallely feeding a single network by two or more generators, converters or transformers
    • H02J3/46—Controlling of the sharing of output between the generators, converters, or transformers
    • H02J3/48—Controlling the sharing of the in-phase component
    • H—ELECTRICITY
    • H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
    • H02J3/00—Circuit arrangements for ac mains or ac distribution networks
    • H02J3/38—Arrangements for parallely feeding a single network by two or more generators, converters or transformers
    • H02J3/46—Controlling of the sharing of output between the generators, converters, or transformers
    • H02J3/50—Controlling the sharing of the out-of-phase component
    • H—ELECTRICITY
    • H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
    • H02J2203/00—Indexing scheme relating to details of circuit arrangements for AC mains or AC distribution networks
    • H02J2203/20—Simulating, e g planning, reliability check, modelling or computer assisted design [CAD]

Abstract

The present invention relates to a kind of uncertain active distribution network running optimizatin methods of adaptation power supply power output, belong to power system active power distribution network technical field.The stochastic uncertainty of power supply power output is portrayed in the form of scene, consider that standby configuration, spare response participate in consumption stochastic uncertainty mechanism, with the minimum target of active distribution network operating cost, in the case where meeting the power distribution network safe operation technical requirements under stochastic uncertainty scene, building adapts to the probabilistic active distribution network optimal operation model of power supply power output, and selects CONOPT solver to solve mentioned Optimized model based on GAMS Optimization Platform.The invention proposes adapt to the probabilistic active distribution network optimal operation model of power supply power output, illustrate that voltage character of load exchanges reduction distribution system operating cost, reduction power with higher level's power transmission network and uncertain consumption is promoted to have positive influence, and load voltage response characteristic ratio is bigger, positive influence is more significant.

Description

A kind of uncertain active distribution network running optimizatin method of adaptation power supply power output
Technical field
The present invention relates to a kind of uncertain active distribution network running optimizatin methods of adaptation power supply power output, belong to electric system Active distribution network technical field.
Background technique
The output power of power supply of the cleaning form such as wind-powered electricity generation, photovoltaic power generation has the characteristics that stochastic uncertainty, extensive Feed-in power distribution network brings stern challenge to dispatching of power netwoks operation.Therefore, how flexible resource in power distribution network is sufficiently called, protected Under the premise of hindering electric power netting safe running, the development of renewable energies such as wind-powered electricity generation, the photovoltaic power generation of stochastic uncertainty are as often as possible dissolved Electricity is current active distribution network running optimizatin hot issue urgently to be resolved.
Prediction is to hold the probabilistic important means of renewable energy power generation, and there are many prediction techniques, such as area at present Between prediction, probabilistic forecasting etc., but its precision of prediction is still very low.As a result, only only in accordance with can in power distribution network optimal operation model The Optimal Decision-making of the being determined property of desired value of renewable source of energy generation will make the biggish operation risk of systems face, contain high permeability The active distribution network running optimizatin of renewable energy power generation must move towards unascertained decision.Based on prediction to probabilistic Modeling is that active distribution network rationally dissolves probabilistic prerequisite.It is usually selected about probabilistic mathematical expression fuzzy Number, interval number and random number form.Fuzzy number expression needs to use complicated subordinating degree function, and interval number expression is needed by robust Optimization realizes, to meet in the fluctuation range of section whole feasible, and the way is overly conservative, and random number then needs to use probability point Cloth functional form, convenient for being mutually connected with prediction techniques such as existing probabilistic forecasting, interval predictions.How to call main in power distribution network It is that active distribution network dissolves probabilistic key that dynamic flexible resource, which is reduced to the pressure regulation of higher level's power transmission network, the dependence of frequency modulation resource,.
Summary of the invention
The technical problem to be solved in the present invention is to provide a kind of uncertain active distribution network operation of adaptation power supply power output is excellent Change method, by calling initiative and flexible resource in power distribution network to reduce to the pressure regulation of higher level's power transmission network, the dependence of frequency modulation resource.
The technical solution adopted by the present invention is that: a kind of uncertain active distribution network running optimizatin side of adaptation power supply power output Method includes the following steps:
Step 1: the stochastic uncertainty of power supply power output is portrayed in the form of scene；
Step 2: considering that standby configuration, spare response participate in the mechanism of consumption stochastic uncertainty, transported with active distribution network Row cost minimization is target；
Step 3: in the case where meeting the power distribution network safe operation technical requirements under stochastic uncertainty scene, building adapts to power supply It contributes probabilistic active distribution network optimal operation model；
Step 4: selecting CONOPT solver Optimized model mentioned to step 3 to solve based on GAMS Optimization Platform.
Specifically, specific step is as follows for the step 1:
Step1: scene production
In entire distribution system the uncertainty of renewable energy power generation be exactly it is consistent with renewable energy quantity with The combination of machine time series, if renewable power supply number is r, look-ahead time range includes t period, then renewable power supply is contributed It is represented byWhereinIndicate that r-th of renewable power supply exists The output power of period t, a specific implementation of Random time sequence are renewable power supply power within the scope of look-ahead time One scene, each scene within the scope of obtained prediction desired value and prediction deviation will be endowed certain probability value, such as S-th of sceneIts probability value note For ρs, each renewable power supply is above and below the fluctuation range that accordingly there are a desired value and prediction deviation characterization in the given period Limit, random sampling corresponds to a specific implementation in the range as a result, and all renewable power supplies are taken out at random in all periods The combination of sample constitutes a scene, generates Ns scene at random using which, then the probability of each scene is
Step2: scene reduction
If scene i is labeled as P(i), scene j is labeled as P(j), probability of happening is denoted as ρ respectivelyiAnd ρj, define two scenes it Between distance be two scene vector differences second order norm:
Dij=| | P(i)-P(j)||2 (1)
In formula, DijAs the distance between scene i and scene j, the purpose of scene reduction, which is to select, can most represent the whole audience It is minimum to pursue formula (1) under conditions of the scene number J of deletion is given for the scene subclass of scape set.
Specifically, active distribution network operating cost minimum target function in the step 2 are as follows:
In formula, JGThe set constituted for adjustable synchronous units all in distribution system；JWFor double-fed wind all in distribution system The set that motor group is constituted；JVThe set constituted for photovoltaic generating systems all in distribution system；JSIt is all defeated in distribution system The set that electric device is constituted；Pex,0Active power for active distribution network from higher level's power transmission network, Cex,0Indicate the limit of root node Electricity price；PgFor the active power of gas turbine g output；ρsIndicate scene probability, Indicate gas turbine Non-firm power；Indicate the non-firm power of s scene higher level's power transmission network；Indicate s scene abandonment electrical power； S scene photovoltaic generating system is marked to abandon optical quantum；Cg() is the cost of electricity-generating characterisitic function of gas turbine g；For combustion gas wheel The stand-by cost characterisitic function of machine；For the abandoning light cost function of photovoltaic generating system；For the abandoning of double-fed fan motor unit Eolian characterisitic function；For the stand-by cost characterisitic function of higher level's power transmission network；Δ τ indicates period lasts duration.
Specifically, the electrical physical constraint and safety specifications that the step 3 need to be met with distribution system operation are about Beam condition, constraint condition is divided into equality constraint and inequality constraints, specific as follows:
1. equality constraint
(1) basic point power-balance constraint:
Wherein, PwFor the active power output of Wind turbines scheduling, PvFor the active power output of photovoltaic generating system scheduling, PdFor electricity The active power demand of power load, JDThe set constituted for electric loads all in distribution system；
(2) node active power, reactive power equilibrium constraint:
In formula (4), Pl(s) and QlIt (s) is respectively the active power and reactive power transmitted on transmission of electricity element l under s scene, It is represented by formula (5), formula (6) respectively；JS,iFor the set constituted using node i as all transmission of electricity elements of first node；JE,iFor The set constituted using node i as all transmission of electricity elements of end-node；JNThe set constituted for nodes all in distribution system；JD,i The set constituted for electric loads all in node i；Pi(s) and Qi(s) be respectively s scene bet ingress i active power and Reactive power is represented by formula (7), formula (8) respectively；Pd(s) and QdIt (s) is respectively electric load d active at scene s Power, reactive power demand；
Pl(s)=Vi 2(s)gl-Vi(s)Vj(s)·(glcosθij(s)+blsinθij(s)) (5)
Ql(s)=- Vi 2(s)bl+Vi(s)Vj(s)·(blcosθij(s)-glsinθij(s)) (6)
In formula, θij(s) phase angle difference of scene s lower node i and node j voltage phasor are indicated；Vi(s) it indicates to prop up under scene s The first node i-node voltage magnitude of road l；Vj(s) the end-node j node voltage amplitude of branch l under scene s is indicated；glAnd blRespectively For the electric conductivity value and susceptance value of the branch l that transmits electricity；
Wherein, JG,iThe set constituted for synchronous gas turbines all in node i；JR,iFor renewable electricity all in node i The set that source is constituted；JC,iThe set constituted for reactive-load compensation equipments all in node i；Pg(s) and QgIt (s) is respectively under s scene The active power and reactive power of synchronous gas turbine g output；Pr(s) and QrIt (s) is respectively renewable power supply r output under s scene Active power and reactive power；Qc(s) lagging reactive power exported at scene s for reactive-load compensation equipment c；
(3) schedulable synchronous unit relevant equations constraint:
Wherein,For unit g stator reactive current synchronous under s scene；Vi(s) it indicates under s scene where synchronous unit g The voltage magnitude of node i；It is expressed as the floating voltage of synchronous unit g excitation con-trol setting；KgIndicate the voltage of synchronous unit g Difference coefficient；Eg(s) built-in potential of synchronous unit g under s scene is indicated；δg(s) the generator rotor angle value of synchronous unit g under s scene is indicated；For synchronous unit g d-axis reactance；
(4) voltage character of load:
In formula, QdFor the reactive power demand of electric load；WithRespectively electric load d is in voltage rating level Under active power, reactive power demand；WithRespectively the constant impedance active power of electric load d, reactive power portion Point；WithRespectively constant current active power, the reactance capacity of electric load d；WithRespectively electric load Firm power active power, the reactance capacity of d；V (s) is system voltage under s scene；V0For system nominal voltage；
2. inequality constraints
(1) the active basic point power bound constraint of schedulable synchronous unit (small power station's unit, gas turbine):
Wherein, PgFor the basic point power of schedulable synchronous machine unit scheduling；WithWhat respectively synchronous unit g allowed has Function power bound；
(2) Wind turbines active power output constrains under basic point mode:
Wherein,For the desired active power output of Wind turbines；
(3) photovoltaic generating system active power output constrains under basic point mode:
Wherein,For the desired active power output of photovoltaic generating system；
(4) the marginal capacity range constraint that higher level's power transmission network provides:
Wherein,The maximum marginal capacity provided for higher level's power transmission network；
(5) schedulable synchronous unit spare capacity range constraint:
Wherein,The spare capacity upper limit value provided for synchronous unit；
(6) schedulable synchronous unit output power range constraint:
Wherein,WithRespectively schedulable synchronous unit g excitation potential bound；δg(s) schedulable synchronization is indicated Generator rotor angle value of the unit g under s scene；For synchronous unit g stator current maximum value；
(7) double-fed fan motor unit range of operation constrains
In formula, PwIt (s) is the active power output of double-fed fan motor unit w Wind turbines under s scene；QwIt (s) is double-fed fan motor machine The idle power output of group w Wind turbines under s scene；Indicate that double-fed fan motor unit w is limited under s scene by natural conditions And the maximum active power that can be exported；VwIt (s) is double-fed fan motor unit w stator side set end voltage under s scene；For double-fed The stator reactance of Wind turbines w；For the excitation reactance of double-fed fan motor unit w；It is maximum for double-fed fan motor unit w rotor-side Electric current；
(8) photovoltaic generating system range of operation constrains
Wherein: Pv(s) active power output of the photovoltaic generating system v under s scene is indicated；Qv(s) photovoltaic generating system v is indicated Idle power output under s scene；Indicate that photovoltaic generating system v is limited by natural conditions under s scene and can be exported Maximum active power；Vv(s) indicate photovoltaic generating system v under s scene and net side node voltage；Indicate photovoltaic hair The maximum carrying electric current of electric system v inverter；
(9) node voltage amplitude bound constrains:
In formula (26), Vi maxAnd Vi minRespectively indicate the bound of node i voltage magnitude；Vi(s) for s scene lower node i's Voltage magnitude；
(10) the thermocurrent range constraint that transmission of electricity element allows:
In formula (27),Indicate the maximum current of transmission of electricity element l, Ilij(s) electric current of transmission of electricity element l under s scene is indicated Amplitude may be expressed as:
In formula (28), YlIndicate the admittance modulus value of transmission of electricity element l；
(11) reactive-load compensation equipment capacity bound constrains:
In formula (29), Qc(s) reactive power that reactive-load compensation equipment c is compensated under s scene is indicated,WithPoint It Wei not reactive-load compensation equipment c maximum compensating power and minimum compensating power under s scene；JCFor distribution system institute whether there is or not Function compensates the set that equipment is constituted.
The beneficial effects of the present invention are:
(1) stochastic uncertainty of power supply power output, including scene production and scene reduction are portrayed in the form of scene；
(2) mechanism of standby configuration, spare response participation consumption stochastic uncertainty is considered, running with distribution system needs The electrical physical constraint and safety specifications met is constraint condition, meets the power distribution network safety under stochastic uncertainty scene Running technology requirement fully considers the influence that voltage characteristic is distributed power balance system and trend, on the basis of random scene To pursue the minimum target of distribution system operating cost, building adapts to the uncertain active distribution network running optimizatin mould of power supply power output Type；
(3) present invention considers voltage character of load, dissolves the influence of uncertain mechanism etc., expands power distribution network operation The feasible zone of Optimized model or optimizing space, promote grid operation quality to a certain extent, avoid the operation of conventional electrical distribution net excellent The deficiencies of changing the conservative and limitation of decision, realization carry out in advance distribution system active balance mode and voltage support mode Optimal Decision-making.
Detailed description of the invention
Fig. 1 is overall flow figure of the present invention；
Fig. 2 scene reduces flow chart；
The IEEE33 Node power distribution system electric hookup of Fig. 3 modification.
Specific embodiment
In order to illustrate the objectives, technical solutions, and advantages of the present invention, with reference to the accompanying drawings and examples, to the present invention make into One step is described in detail.
Embodiment 1: a kind of uncertain active distribution network running optimizatin method as shown in Figure 1-3, adaptation power supply is contributed, packet Include following steps:
Step 1: the stochastic uncertainty of power supply power output is portrayed in the form of scene；
Step 2: considering that standby configuration, spare response participate in the mechanism of consumption stochastic uncertainty, transported with active distribution network Row cost minimization is target；
Step 3: in the case where meeting the power distribution network safe operation technical requirements under stochastic uncertainty scene, building adapts to power supply It contributes probabilistic active distribution network optimal operation model；
Step 4: selecting CONOPT solver Optimized model mentioned to step 3 to solve based on GAMS Optimization Platform.
Specifically, specific step is as follows for the step 1:
Step1: scene production
In entire distribution system the uncertainty of renewable energy power generation be exactly it is consistent with renewable energy quantity with The combination of machine time series, if renewable power supply number is r, look-ahead time range includes t period, then renewable power supply is contributed It is represented byWhereinIndicate that r-th of renewable power supply exists The output power of period t, a specific implementation of Random time sequence are renewable power supply power within the scope of look-ahead time One scene, each scene within the scope of obtained prediction desired value and prediction deviation will be endowed certain probability value, such as S-th of sceneIts probability value note For ρs, each renewable power supply is above and below the fluctuation range that accordingly there are a desired value and prediction deviation characterization in the given period Limit, random sampling corresponds to a specific implementation in the range as a result, and all renewable power supplies are taken out at random in all periods The combination of sample constitutes a scene, generates Ns scene at random using which, then the probability of each scene is
Step2: scene reduction
If scene i is labeled as P(i), scene j is labeled as P(j), probability of happening is denoted as ρ respectivelyiAnd ρj, define two scenes it Between distance be two scene vector differences second order norm:
Dij=| | P(i)-P(j)||2 (1)
In formula, DijAs the distance between scene i and scene j, the purpose of scene reduction, which is to select, can most represent the whole audience It is minimum to pursue formula (1) under conditions of the scene number J of deletion is given for the scene subclass of scape set.
Scene reduction realization comprising the following steps:
1) the number of iterations k=0, the set C that reduction scene is constituted are initializedkIt is initialized as empty set；
2) the scene P for calculating kth time iteration reduction is traversed according to formula (4-2)k；
3) set that reduction scene is constituted, C are updatedk=Ck-1∪{Pk}；
4) the number of iterations updates, k=k+1.
5) judge whether to reach maximum number of iterations, k < α?
6) update reduction after scene set probability distribution, by the scene probability of reduction be superimposed to reservation with a distance from it most On close scene probability；
7) the scene set after obtaining final reduction.
Adapting to the uncertain active distribution network optimal operation model of power supply power output is on the basis of random scene with power distribution system The electrical physical constraint and safety specifications that system operation need to meet are constraint condition, fully consider that distribution system response is uncertain Property mechanism, the influence that voltage characteristic is distributed power balance system and trend is considered, to pursue distribution system operating cost most Small is target, and realization carries out advanced Optimal Decision-making to distribution system active balance mode and voltage support mode.
Further, adapting to the uncertain active distribution network optimal operation model of power supply power output is on the basis of random scene To pursue the minimum target of distribution system operating cost, active distribution network operating cost minimum target function in the step 2 are as follows:
In formula, JGThe set constituted for adjustable synchronous units (small power station's unit, gas turbine) all in distribution system；JW The set constituted for double-fed fan motor units all in distribution system；JVThe collection constituted for photovoltaic generating systems all in distribution system It closes；JSThe set constituted for transmission of electricity elements all in distribution system；Pex,0Wattful power for active distribution network from higher level's power transmission network Rate, Cex,0Indicate the Marginal Pricing of root node；PgFor the active power of gas turbine g output；ρsIndicate scene probability, Indicate gas turbine non-firm power；Indicate the non-firm power of s scene higher level's power transmission network；Indicate s scene abandonment electrical power；S scene photovoltaic generating system is marked to abandon optical quantum；Cg() is gas turbine g Cost of electricity-generating characterisitic function；For the stand-by cost characterisitic function of gas turbine；For the abandoning light of photovoltaic generating system Cost function；For the abandonment cost behavior function of double-fed fan motor unit；For the stand-by cost characteristic of higher level's power transmission network Function；Δ τ indicates period lasts duration.
Further, which fully considers that distribution system responds probabilistic mechanism, considers voltage characteristic pair The influence of power balance system and trend distribution runs the electrical physical constraint and safety specifications that need to meet with distribution system For constraint condition, realization carries out advanced Optimal Decision-making to distribution system active balance mode and voltage support mode.
Equality constraint and inequality constraints in the step 3, specific as follows:
1. equality constraint
(1) basic point power-balance constraint:
Wherein, PwFor the active power output of Wind turbines scheduling, PvFor the active power output of photovoltaic generating system scheduling, PdFor electricity The active power demand of power load, JDThe set constituted for electric loads all in distribution system；
(2) node active power, reactive power equilibrium constraint (power flow equation):
In formula (4), Pl(s) and QlIt (s) is respectively the active power and reactive power transmitted on transmission of electricity element l under s scene, It is represented by formula (5), formula (6) respectively；JS,iFor the set constituted using node i as all transmission of electricity elements of first node；JE,iFor The set constituted using node i as all transmission of electricity elements of end-node；JNThe set constituted for nodes all in distribution system；JD,i The set constituted for electric loads all in node i；Pi(s) and Qi(s) be respectively s scene bet ingress i active power and Reactive power is represented by formula (7), formula (8) respectively；Pd(s) and QdIt (s) is respectively electric load d active at scene s Power, reactive power demand；
Pl(s)=Vi 2(s)gl-Vi(s)Vj(s)·(glcosθij(s)+blsinθij(s)) (5)
Ql(s)=- Vi 2(s)bl+Vi(s)Vj(s)·(blcosθij(s)-glsinθij(s)) (6)
In formula, θij(s) phase angle difference of scene s lower node i and node j voltage phasor are indicated；Vi(s) it indicates to prop up under scene s The first node i-node voltage magnitude of road l；Vj(s) the end-node j node voltage amplitude of branch l under scene s is indicated；glAnd blRespectively For the electric conductivity value and susceptance value of the branch l that transmits electricity；
Wherein, JG,iThe set constituted for synchronous gas turbines all in node i；JR,iFor renewable electricity all in node i The set that source is constituted；JC,iThe set constituted for reactive-load compensation equipments all in node i；Pg(s) and QgIt (s) is respectively under s scene The active power and reactive power of synchronous gas turbine g output；Pr(s) and QrIt (s) is respectively renewable power supply r output under s scene Active power and reactive power；Qc(s) lagging reactive power exported at scene s for reactive-load compensation equipment c；
(3) schedulable synchronous unit (small power station's unit, gas turbine) relevant equations constraint:
Wherein,For unit g stator reactive current synchronous under s scene；Vi(s) it indicates under s scene where synchronous unit g The voltage magnitude of node i；It is expressed as the floating voltage of synchronous unit g excitation con-trol setting；KgIndicate the voltage of synchronous unit g Difference coefficient；Eg(s) built-in potential of synchronous unit g under s scene is indicated；δg(s) the generator rotor angle value of synchronous unit g under s scene is indicated；For synchronous unit g d-axis reactance；
(4) voltage character of load:
In formula, QdFor the reactive power demand of electric load；WithRespectively electric load d is in voltage rating level Under active power, reactive power demand；WithRespectively the constant impedance active power of electric load d, reactive power portion Point；WithRespectively constant current active power, the reactance capacity of electric load d；WithRespectively electric load Firm power active power, the reactance capacity of d；V (s) is system voltage under s scene；V0For system nominal voltage；
2. inequality constraints
(1) the active basic point power bound constraint of schedulable synchronous unit (small power station's unit, gas turbine):
Wherein, PgFor the basic point power of schedulable synchronous machine unit scheduling；WithWhat respectively synchronous unit g allowed has Function power bound；
(2) Wind turbines active power output constrains under basic point mode:
Wherein,For the desired active power output of Wind turbines；
(3) photovoltaic generating system active power output constrains under basic point mode:
Wherein,For the desired active power output of photovoltaic generating system；
(4) the marginal capacity range constraint that higher level's power transmission network provides:
Wherein,The maximum marginal capacity provided for higher level's power transmission network；
(5) schedulable synchronous unit spare capacity range constraint:
Wherein,The spare capacity upper limit value provided for synchronous unit；
(6) schedulable synchronous unit output power range constraint:
Wherein,WithRespectively schedulable synchronous unit g excitation potential bound；δg(s) schedulable synchronization is indicated Generator rotor angle value of the unit g under s scene；For synchronous unit g stator current maximum value；
(7) double-fed fan motor unit range of operation constrains
In formula, PwIt (s) is the active power output of double-fed fan motor unit w Wind turbines under s scene；QwIt (s) is double-fed fan motor machine The idle power output of group w Wind turbines under s scene；Indicate that double-fed fan motor unit w is limited under s scene by natural conditions And the maximum active power that can be exported；VwIt (s) is double-fed fan motor unit w stator side set end voltage under s scene；For double-fed The stator reactance of Wind turbines w；For the excitation reactance of double-fed fan motor unit w；Most for double-fed fan motor unit w rotor-side High current；
(8) photovoltaic generating system range of operation constrains
Wherein: Pv(s) active power output of the photovoltaic generating system v under s scene is indicated；Qv(s) photovoltaic generating system v is indicated Idle power output under s scene；Indicate that photovoltaic generating system v is limited by natural conditions under s scene and can be exported Maximum active power；Vv(s) indicate photovoltaic generating system v under s scene and net side node voltage；Indicate photovoltaic hair The maximum carrying electric current of electric system v inverter；
(9) node voltage amplitude bound constrains:
In formula (26), Vi maxAnd Vi minRespectively indicate the bound of node i voltage magnitude；Vi(s) for s scene lower node i's Voltage magnitude；
(10) the thermocurrent range constraint that transmission of electricity element allows:
In formula (27),Indicate the maximum current of transmission of electricity element l, Il,ij(s) electric current of transmission of electricity element l under s scene is indicated Amplitude may be expressed as:
In formula (28), YlIndicate the admittance modulus value of transmission of electricity element l；
(11) reactive-load compensation equipment capacity bound constrains:
In formula (29), Qc(s) reactive power that reactive-load compensation equipment c is compensated under s scene is indicated,WithPoint It Wei not reactive-load compensation equipment c maximum compensating power and minimum compensating power under s scene；JCFor distribution system institute whether there is or not Function compensates the set that equipment is constituted.
The present invention is further illustrated With reference to embodiment.
Based on the present invention utilizes IEEE33 Node power distribution system, wind-powered electricity generation, photovoltaic distributed development of renewable energy are added Electricity, the synchronous form of power source such as addition small power station, miniature gas turbine is to verify effectiveness of the invention.The IEEE33 node of modification For distribution system electric hookup as shown in figure 3, in Simulation Example analysis, which is selected as 10MVA, benchmark electricity Pressure is selected as 12.66kV, it is assumed that the deviation range of node voltage is ± the 5% of voltage rating, and optimizing cycle is selected as 15min.Through imitative It is true to calculate, unit, Wind turbines and photovoltaic generating system running optimizatin result are synchronized respectively as shown in table 1 to table 3, the two target Functional value exchanges power and the comparison of spare capacity value from higher level's power transmission network as shown in table 4.
Table 1 synchronizes unit running optimization result
2 double-fed fan motor unit running optimizatin result of table
3 photovoltaic generating system running optimizatin result of table
Table 4 considers and does not consider the running optimizatin comparing result of voltage characteristic
Through comparison it is found that compared with not using power distribution network randomized optimization process of the invention, using distribution of the invention Net running optimizatin method decision optimization operation totle drilling cost, with higher level's power transmission network exchange power be intended to low and total abandonment, Abandoning light desired value reduces, and thus illustrates to consider that voltage regulation properties participate in responding probabilistic economic benefit.
The invention proposes the probabilistic active distribution network optimal operation model of power supply power output is adapted to, load electricity is illustrated Pressure characteristic exchanges reduction distribution system operating cost, reduction power with higher level's power transmission network and uncertain consumption is promoted to have Positive influence, and load voltage response characteristic ratio is bigger, positive influence is more significant, it expands for mathematical optimization angle Feasible zone range is conducive to promote distribution system on-road efficiency.
In conjunction with attached drawing, the embodiment of the present invention is explained in detail above, but the present invention is not limited to above-mentioned Embodiment within the knowledge of a person skilled in the art can also be before not departing from present inventive concept Put that various changes can be made.

Claims (4)

  1. A kind of uncertain active distribution network running optimizatin method 1. adaptation power supply is contributed, characterized by the following steps:
    Step 1: the stochastic uncertainty of power supply power output is portrayed in the form of scene；
    Step 2: consider standby configuration, it is spare response participate in consumption stochastic uncertainty mechanism, with active distribution network operation at This minimum target；
    Step 3: in the case where meeting the power distribution network safe operation technical requirements under stochastic uncertainty scene, building adapts to power supply power output Uncertain active distribution network optimal operation model；
    Step 4: selecting CONOPT solver Optimized model mentioned to step 3 to solve based on GAMS Optimization Platform.
  2. 2. a kind of adaptation power supply according to claim 1 is contributed, uncertain active distribution network runs randomized optimization process, It is characterized by: specific step is as follows for step 1:
    Step1: scene production
    When the uncertainty of renewable energy power generation is exactly consistent with renewable energy quantity random in entire distribution system Between sequence combination, if renewable power supply number be r, look-ahead time range include t period, then renewable power supply contribute can table It is shown asWhereinIndicate r-th of renewable power supply in the period The output power of t, a specific implementation of Random time sequence are one of renewable power supply power within the scope of look-ahead time Scene, each scene within the scope of obtained prediction desired value and prediction deviation will be endowed certain probability value, and such as s-th SceneIts probability value is denoted as ρs, Each renewable power supply accordingly has the fluctuation range bound an of desired value and prediction deviation characterization in the given period, by This, random sampling corresponds to a specific implementation in the range, and all renewable power supplies are in all period random samplings Combination constitutes a scene, generates Ns scene at random using which, then the probability of each scene is
    Step2: scene reduction
    If scene i is labeled as P(i), scene j is labeled as P(j), probability of happening is denoted as ρ respectivelyiAnd ρj, define between two scenes Distance is the second order norm of two scene vector differences:
    Dij=| | P(i)-P(j)||2 (1)
    In formula, DijAs the distance between scene i and scene j, the purpose of scene reduction, which is to select, can most represent whole scene collection It is minimum to pursue formula (1) under conditions of the scene number J of deletion is given for the scene subclass of conjunction.
  3. The uncertain active distribution network running optimizatin method 3. a kind of adaptation power supply according to claim 1 is contributed, it is special Sign is: active distribution network operating cost minimum target function in step 2 are as follows:
    In formula, JGThe set constituted for adjustable synchronous units all in distribution system；JWFor double-fed fan motor machines all in distribution system The set that group is constituted；JVThe set constituted for photovoltaic generating systems all in distribution system；JSFor transmission of electricity members all in distribution system The set that part is constituted；Pex,0Active power for active distribution network from higher level's power transmission network, Cex,0Indicate the limit electricity of root node Valence；PgFor the active power of gas turbine g output；ρsIndicate scene probability, Indicate that gas turbine is spare Power；Indicate the non-firm power of s scene higher level's power transmission network；Indicate s scene abandonment electrical power；Mark s Scene photovoltaic generating system abandons optical quantum；Cg() is the cost of electricity-generating characterisitic function of gas turbine g；For the standby of gas turbine With cost behavior function；For the abandoning light cost function of photovoltaic generating system；For the abandonment cost of double-fed fan motor unit Characterisitic function；For the stand-by cost characterisitic function of higher level's power transmission network；Δ τ indicates period lasts duration.
  4. The uncertain active distribution network running optimizatin method 4. a kind of adaptation power supply according to claim 3 is contributed, it is special Sign is: the electrical physical constraint and safety specifications that step 3 need to be met using distribution system operation is constraint conditions, constraint Condition is divided into equality constraint and inequality constraints, specific as follows:
    1. equality constraint
    (1) basic point power-balance constraint:
    Wherein, PwFor the active power output of Wind turbines scheduling, PvFor the active power output of photovoltaic generating system scheduling, PdFor power load The active power demand of lotus, JDThe set constituted for electric loads all in distribution system；
    (2) node active power, reactive power equilibrium constraint:
    In formula (4), Pl(s) and QlIt (s) is respectively the active power and reactive power transmitted on transmission of electricity element l under s scene, point It is not represented by formula (5), formula (6)；JS,iFor the set constituted using node i as all transmission of electricity elements of first node；JE,iFor with section Point i is the set that all transmission of electricity elements of end-node are constituted；JNThe set constituted for nodes all in distribution system；JD,iFor section The set that all electric loads are constituted on point i；Pi(s) and QiIt (s) is respectively the s scene bet active power of ingress i and idle Power is represented by formula (7), formula (8) respectively；Pd(s) and Qd(s) be respectively active power of the electric load d at scene s, Reactive power demand；
    Pl(s)=Vi 2(s)gl-Vi(s)Vj(s)·(glcosθij(s)+blsinθij(s)) (5)
    Ql(s)=- Vi 2(s)bl+Vi(s)Vj(s)·(blcosθij(s)-glsinθij(s)) (6)
    In formula, θij(s) phase angle difference of scene s lower node i and node j voltage phasor are indicated；Vi(s) branch l under scene s is indicated First node i-node voltage magnitude；Vj(s) the end-node j node voltage amplitude of branch l under scene s is indicated；glAnd blIt is respectively defeated The electric conductivity value and susceptance value of electric branch l；
    Wherein, JG,iThe set constituted for synchronous gas turbines all in node i；JR,iFor renewable power supply structures all in node i At set；JC,iThe set constituted for reactive-load compensation equipments all in node i；Pg(s) and QgIt (s) is respectively synchronous under s scene The active power and reactive power of gas turbine g output；Pr(s) and Qr(s) be respectively under s scene renewable power supply r export have Function power and reactive power；Qc(s) lagging reactive power exported at scene s for reactive-load compensation equipment c；
    (3) schedulable synchronous unit relevant equations constraint:
    Wherein,For unit g stator reactive current synchronous under s scene；Vi(s) node where synchronous unit g under s scene is indicated The voltage magnitude of i；It is expressed as the floating voltage of synchronous unit g excitation con-trol setting；KgIndicate that the voltage tune of synchronous unit g is poor Coefficient；Eg(s) built-in potential of synchronous unit g under s scene is indicated；δg(s) the generator rotor angle value of synchronous unit g under s scene is indicated；For Synchronous unit g d-axis reactance；
    (4) voltage character of load:
    In formula, QdFor the reactive power demand of electric load；WithRespectively electric load d having in the case where voltage rating is horizontal Function power, reactive power demand；WithRespectively constant impedance active power, the reactance capacity of electric load d； WithRespectively constant current active power, the reactance capacity of electric load d；WithThe respectively perseverance of electric load d Determine power active power, reactance capacity；V (s) is system voltage under s scene；V0For system nominal voltage；
    2. inequality constraints
    (1) the active basic point power bound constraint of schedulable synchronous unit:
    Wherein, PgFor the basic point power of schedulable synchronous machine unit scheduling；WithThe wattful power that respectively synchronous unit g allows Rate bound；
    (2) Wind turbines active power output constrains under basic point mode:
    Wherein,For the desired active power output of Wind turbines；
    (3) photovoltaic generating system active power output constrains under basic point mode:
    Wherein,For the desired active power output of photovoltaic generating system；
    (4) the marginal capacity range constraint that higher level's power transmission network provides:
    Wherein,The maximum marginal capacity provided for higher level's power transmission network；
    (5) schedulable synchronous unit spare capacity range constraint:
    Wherein,The spare capacity upper limit value provided for synchronous unit；
    (6) schedulable synchronous unit output power range constraint:
    Wherein,WithRespectively schedulable synchronous unit g excitation potential bound；δg(s) schedulable synchronous unit g is indicated Generator rotor angle value under s scene；For synchronous unit g stator current maximum value；
    (7) double-fed fan motor unit range of operation constrains
    In formula, PwIt (s) is the active power output of double-fed fan motor unit w Wind turbines under s scene；QwIt (s) is double-fed fan motor unit w The idle power output of Wind turbines under s scene；Indicate that double-fed fan motor unit w is limited under s scene by natural conditions and The maximum active power that can be exported；VwIt (s) is double-fed fan motor unit w stator side set end voltage under s scene；For double-fed wind The stator reactance of motor group w；For the excitation reactance of double-fed fan motor unit w；For the maximum electricity of double-fed fan motor unit w rotor-side Stream；
    (8) photovoltaic generating system range of operation constrains
    Wherein: Pv(s) active power output of the photovoltaic generating system v under s scene is indicated；Qv(s) indicate photovoltaic generating system v at s Idle power output under scape；Indicate the maximum that photovoltaic generating system v is limited under s scene by natural conditions and can be exported Active power；Vv(s) indicate photovoltaic generating system v under s scene and net side node voltage；Indicate photovoltaic generating system v The maximum carrying electric current of inverter；
    (9) node voltage amplitude bound constrains:
    In formula (26),WithRespectively indicate the bound of node i voltage magnitude；Vi(s) voltage for being s scene lower node i Amplitude；
    (10) the thermocurrent range constraint that transmission of electricity element allows:
    In formula (27),Indicate the maximum current of transmission of electricity element l, Il,ij(s) current amplitude of transmission of electricity element l under s scene is indicated, It may be expressed as:
    In formula (28), YlIndicate the admittance modulus value of transmission of electricity element l；
    (11) reactive-load compensation equipment capacity bound constrains:
    In formula (29), Qc(s) reactive power that reactive-load compensation equipment c is compensated under s scene is indicated,WithRespectively s Reactive-load compensation equipment c maximum compensating power and minimum compensating power under scene；JCFor all idle benefits of distribution system Repay the set of equipment composition.
CN201810979860.2A 2018-08-27 2018-08-27 Active power distribution network operation optimization method adaptive to uncertainty of power output Active CN109066808B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201810979860.2A CN109066808B (en) 2018-08-27 2018-08-27 Active power distribution network operation optimization method adaptive to uncertainty of power output

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201810979860.2A CN109066808B (en) 2018-08-27 2018-08-27 Active power distribution network operation optimization method adaptive to uncertainty of power output

Publications (2)

Publication Number Publication Date
CN109066808A true CN109066808A (en) 2018-12-21
CN109066808B CN109066808B (en) 2021-07-16

Family

ID=64757023

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201810979860.2A Active CN109066808B (en) 2018-08-27 2018-08-27 Active power distribution network operation optimization method adaptive to uncertainty of power output

Country Status (1)

Country Link
CN (1) CN109066808B (en)

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2015077196A1 (en) * 2013-11-19 2015-05-28 Siemens Aktiengesellschaft Holistic optimization of distribution automation using survivability modeling to support storm hardening
CN106684889A (en) * 2017-03-24 2017-05-17 河海大学 Random reactive optimization method of active distribution network based on scenario method
CN106953354A (en) * 2017-03-10 2017-07-14 国网山东省电力公司经济技术研究院 Consider the method for Unit Commitment containing wind-powered electricity generation of voltage support

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2015077196A1 (en) * 2013-11-19 2015-05-28 Siemens Aktiengesellschaft Holistic optimization of distribution automation using survivability modeling to support storm hardening
CN106953354A (en) * 2017-03-10 2017-07-14 国网山东省电力公司经济技术研究院 Consider the method for Unit Commitment containing wind-powered electricity generation of voltage support
CN106684889A (en) * 2017-03-24 2017-05-17 河海大学 Random reactive optimization method of active distribution network based on scenario method

Also Published As

Publication number Publication date
CN109066808B (en) 2021-07-16

Similar Documents

Publication Publication Date Title
CN106485352B (en) A kind of multiclass power supply generation schedule combination decision-making method
CN110970906B (en) Method and system for coordinated control of reactive voltage of wind power plant
García et al. Energy management system based on techno-economic optimization for microgrids
Li et al. Stability-constrained two-stage robust optimization for integrated hydrogen hybrid energy system
CN111400937B (en) Comprehensive energy system tide calculation method
CN109103929A (en) Based on the power distribution network economic optimization dispatching method for improving dynamic gram Li Sijin model
CN110007600A (en) A kind of constrained multipotency stream coordinated scheduling aid decision-making system of tool
CN106529737A (en) Planning and distribution method for peak load regulation power source on supply side of power distribution network
CN104167765A (en) Admitting ability distribution-based maximum wind power installed capacity calculation method
CN105790291B (en) A kind of electric system energy storage Optimal Configuration Method based on flexibility assessment
Saadatmand et al. PMU-based FOPID controller of large-scale wind-PV farms for LFO damping in smart grid
CN109524955A (en) A kind of active distribution network optimization method of consideration source lotus voltage characteristic
CN104868468A (en) UPFC optimization configuration method based on life cycle cost
CN109615125A (en) A kind of multizone random production analog method and application considering extra-high voltage peak regulation
CN117154746A (en) Micro-grid day-ahead optimal scheduling method considering uncertainty of photovoltaic output
CN109066808A (en) A kind of uncertain active distribution network running optimizatin method of adaptation power supply power output
CN116436040A (en) Power transmission and distribution collaborative optimization scheduling method and system considering frequency dynamic security constraint
CN111064231A (en) New energy graded interactive consumption method and system
CN109995084A (en) A kind of step power station-thermal power plant&#39;s joint optimal operation method and system
CN109149630A (en) A kind of new energy electricity digestion capability analysis planing method
CN114583762A (en) Daily power and electric quantity balance analysis method considering multi-partition thermal power balance scheduling
CN104466956B (en) The tidal current computing method and device of a kind of distributed generation resource
Liu et al. Sensitivity region based optimization for maximiz-ing renewable generation hosting capacity of an islanded microgrid
CN110198053A (en) It is a kind of to concentrate with the micro-capacitance sensor real-time voltage control method and system combined on the spot
Foote et al. Ubiquitous energy storage

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant