CN108695905A - Safety constraint optimization scheduling method for 110 kV-level-contained power grid - Google Patents
Safety constraint optimization scheduling method for 110 kV-level-contained power grid Download PDFInfo
- Publication number
- CN108695905A CN108695905A CN201810345203.2A CN201810345203A CN108695905A CN 108695905 A CN108695905 A CN 108695905A CN 201810345203 A CN201810345203 A CN 201810345203A CN 108695905 A CN108695905 A CN 108695905A
- Authority
- CN
- China
- Prior art keywords
- power
- unit
- period
- constraint
- 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
Links
- 238000005457 optimization Methods 0.000 title claims abstract description 23
- 238000000034 method Methods 0.000 title claims abstract description 16
- 238000010248 power generation Methods 0.000 claims abstract description 12
- 238000003860 storage Methods 0.000 claims description 29
- 238000009987 spinning Methods 0.000 claims description 23
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims description 12
- 230000005611 electricity Effects 0.000 claims description 9
- 230000009194 climbing Effects 0.000 claims description 7
- 239000000446 fuel Substances 0.000 claims description 6
- 230000006870 function Effects 0.000 claims description 5
- 230000006872 improvement Effects 0.000 claims description 5
- 230000005540 biological transmission Effects 0.000 claims description 4
- 238000002347 injection Methods 0.000 claims description 4
- 239000007924 injection Substances 0.000 claims description 4
- 230000001186 cumulative effect Effects 0.000 claims description 3
- 230000008569 process Effects 0.000 claims description 3
- 230000004044 response Effects 0.000 claims description 2
- 239000004576 sand Substances 0.000 claims description 2
- 238000009825 accumulation Methods 0.000 claims 1
- 238000005086 pumping Methods 0.000 claims 1
- 239000000243 solution Substances 0.000 description 15
- 238000004364 calculation method Methods 0.000 description 7
- 238000004458 analytical method Methods 0.000 description 3
- 230000008859 change Effects 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 238000009826 distribution Methods 0.000 description 2
- 235000013399 edible fruits Nutrition 0.000 description 2
- 238000009472 formulation Methods 0.000 description 2
- 239000000203 mixture Substances 0.000 description 2
- 241001269238 Data Species 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000000052 comparative effect Effects 0.000 description 1
- 230000007812 deficiency Effects 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000013178 mathematical model Methods 0.000 description 1
- 230000015654 memory Effects 0.000 description 1
- 238000003012 network analysis Methods 0.000 description 1
- 238000003825 pressing Methods 0.000 description 1
- 230000003068 static effect Effects 0.000 description 1
- 238000012795 verification Methods 0.000 description 1
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
- 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/04—Circuit arrangements for ac mains or ac distribution networks for connecting networks of the same frequency but supplied from different sources
- H02J3/06—Controlling transfer of power between connected networks; Controlling sharing of load between connected networks
-
- 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]
Landscapes
- Engineering & Computer Science (AREA)
- Power Engineering (AREA)
- Supply And Distribution Of Alternating Current (AREA)
Abstract
The invention discloses a safety constraint optimization scheduling method for a 110 kV-level-contained power grid, which comprises the steps of firstly, improving a direct current power flow model, aiming at the simplified condition that the ratio of the reactance and the resistance of a line in the 110 kV-level power grid is relatively small and does not meet X > R, so that the resistance of the line cannot be ignored, the voltage amplitude of each node is approximately considered to be near a rated value, the phase angle difference of the voltage of the nodes at two ends of the line is very small and is approximately 0, and a grounding branch in the power grid is ignored, so that the alternating current power flow model is simplified, and the accuracy of the calculated result of the line power is higher; the improved direct current power flow model is applied to the optimization scheduling problem of power generation before the day of safety constraint of the power grid with the 110kV voltage level for solving, the solving speed is high, and the solving result is more reasonable and accurate.
Description
Technical field
The present invention relates to technical field of electric power, and in particular to a kind of security constraint of level Grid containing 110kV Optimized Operation side
Method.
Background technology
The Plan Curve of power generation dispatching a few days ago for formulating each generating set in power grid is the key that in Operation of Electric Systems scheduling
Link, and safety is the most basic requirement of operation of power networks, it is therefore desirable to net is carried out to the power generation dispatching Plan Curve of formulation
Network Security Checking ensures that the active power of the corresponding all transmission lines of electricity of day part is no more than its maximum transmission power.Electric power
System load flow calculating is the grid structure in given system, the service condition of branch parameters and load and generator node power
Under, calculate transimission power, the power attenuation of network of the voltage magnitude of each node and phase angle in electric system, each circuit
Etc. states, can be used for checking network element whether overload, each node voltage whether meet the requirements, the distribution of power and point
With whether reasonable etc., therefore, can be applied to carry out network security check in power system optimal dispatch problem.Load flow calculation can
To be divided into AC power flow and DC power flow according to the difference for solving mathematical model, AC power flow equation is nonlinear, calculates knot
The precision of fruit is higher, can comprehensively reflect the truth of system, but its solution efficiency is low, therefore is generally not particularly suited for
Bulk power grid optimizes the line power Security Checking in traffic control problem.DC power flow is simplified on the basis of AC power flow
It obtains, following simplified process has been done on the basis of AC power flow:1) ignore the ground branch in power grid;2) due to normally transporting
In the case of row, each node voltage amplitude is near rated value in electric system, therefore is approximately considered Vi≈1;3) two end segment of circuit
The voltage phase angle of point differs very little, i.e. θij≈ 0, therefore it is approximately considered sin θij≈θij,cosθij≈1;4) in supergrid
Meet branch reactance and is far longer than branch resistance, i.e. X>>R, therefore branch resistance can be ignored.DC flow model does not consider net
Reactive power distribution in network and loss are one group of linear equations, can effectively reduce calculation amount, and solution efficiency has obtained very big
It improves, is widely used in the electric system Security Checking of generation schedule, electricity market congestion management, security constraint a few days ago at present
In the problems such as economic load dispatching, static security analysis and security constraint Unit Combination.
Currently, China there are multiple key city power grids independent from provincial power network, it is directly subordinate to regional power grid management, example
Such as, Guangzhou, Shenzhen Power Grid are independent from Guangdong Power Grid, are directly subordinate to south electric network management.Key city grid dispatching center
Since the range of administration is smaller, in addition to the operation shape of local 220kV and ratings above power grid belonging to concern in day-to-day operation
Outside state, it is also necessary to the operating status of local 110kV level Grids belonging to concern, to the power generation dispatching Plan Curve of formulation into
When row network security is checked, the safe operation for considering 110kV level Grids is needed.At present in the major network for considering Network Security Constraints
In Optimal Scheduling, it is generally used DC flow model, this is because it is generally acknowledged that major network is typically high voltage network, including
The grid structure of 220kV or more meets above-mentioned X>>4 simplified conditions such as R, computational accuracy are also relatively high.But due to
The ratio of line reactance and resistance is relatively small in 110kV power grids, is unsatisfactory for X>>This simplified condition of R, therefore, when needs pair
Including when 110kV level Grids carry out security constraint Optimized Operation, it, may if still being analyzed using DC flow model
Cause the line power being calculated inaccurate, it is larger inclined to occur with actual conditions when carrying out line power safety check
Difference, the accuracy of the optimal dispatch result influenced, the Optimized Operation result that can cause when serious are discontented with full
The Network Security Constraints on border.Therefore, in comprising 110kV level Grid security constraint Optimal Schedulings, to DC power flow mould
It is necessary that type, which is improved, to keep the result of calculation of line power more accurate, the optimal dispatch knot of acquisition
Fruit is also more in line with actual conditions.
For the power grid containing 110kV grades, since the ratio of line reactance and resistance is relatively small in 110kV racks,
It is unsatisfactory for X>>This simplified condition of R, if directly using DC flow model to calculate line power carries out network security check,
The line power being calculated can be caused inaccurate, it is larger to occur with actual conditions when carrying out line power safety check
Deviation influences the Optimized Operation operation result of power grid.
At present both at home and abroad in power system security constraints Optimal Scheduling, typically directly DC flow model being used to calculate
Line power carries out the Security Checking of network;Also, (it is unsatisfactory for X> to containing 110kV level Grids;>R conditions) security constraint
Optimal Scheduling should use any network analysis model that could accurately calculate line power, and there are no carry out correlative study.
Invention content
In view of the deficiencies of the prior art, it is excellent to be intended to provide a kind of security constraint of level Grid containing 110kV for the purpose of the present invention
Change dispatching method, to obtain the optimal dispatch result for being more in line with actual conditions.
To achieve the above object, the present invention adopts the following technical scheme that:
A kind of security constraint Optimization Scheduling of level Grid containing 110kV, the method includes:
Ignore the ground branch in power grid, the active power for obtaining circuit under AC Ioad flow model is:
Pij=ViVj(Gijcosθij+Bijsinθij)-Vi 2Gij (1)
In formula:PijFor the active power of circuit ij;ViAnd VjThe respectively voltage magnitude of node i and node j;θijFor node
Phase difference of voltage between i and node j;Gij,BijThe respectively transconductance between node i and node j and mutual susceptance;
The active power of circuit under AC Ioad flow model is simplified, obtains improving direct current tide model, the improvement direct current
Damp model includes the active power equation of arbitrary circuit ij:
And the active power balance equation of injection of arbitrary node i:
Security constraint a few days ago generation optimization scheduling of the direct current tide model applied to the power grid of voltage class containing 110kV will be improved to ask
Topic, to obtain optimal dispatch result.
The beneficial effects of the present invention are:
This method first improves DC flow model, for line reactance in 110kV voltage class power grids with
The ratio of resistance is relatively small, is unsatisfactory for X>>This simplified condition of R, therefore line resistance cannot be ignored, and each node voltage
Amplitude is approximately considered near rated value, and circuit both ends node voltage phase angle differs very little, is approximately 0, and ignore in power grid
Ground branch, to simplify to AC Ioad flow model, the accuracy of obtained circuit active power calculating result is higher;It will
Improved DC flow model is applied to the power system security constraints of voltage class containing 110kV to carry out in generation optimization scheduling problem a few days ago
It solves, solving speed is quickly and solving result is more rationally accurate.
Description of the drawings
Fig. 1 is IEEE-9 node system grid structure wiring diagrams;
Fig. 2 is total load prediction graph;
Fig. 3 is the trend chart that out-of-limit line power counts;
Fig. 4 a-4c are the output Plan Curve comparison diagram of unit.
Specific implementation mode
In the following, in conjunction with attached drawing and specific implementation mode, the present invention is described further:
A kind of security constraint Optimization Scheduling of level Grid containing 110kV is present embodiments provided, this method includes
(1) it establishes and improves DC flow model
Ignore the ground branch in power grid, then the active power of circuit is under AC Ioad flow model:
Pij=ViVj(Gijcosθij+Bijsinθij)-Vi 2Gij (1)
In formula:PijFor the active power of circuit ij;ViAnd VjThe respectively voltage magnitude of node i and node j;θijFor node
Phase difference of voltage between i and node j;Gij,BijThe respectively transconductance between node i and node j and mutual susceptance.
The process of simplification:
1) due under normally operating condition, each node voltage amplitude is therefore close near rated value in electric system
Seemingly think Vi≈1,Vj≈1。
2) voltage phase angle of two end node of circuit differs very little, i.e. θij≈ 0, therefore it is approximately considered sin (θij/2)≈θij/
2,cos(θij/2)≈1。
Therefore, circuit active power formula (1) can be reduced to:
And then obtain the injecting power equilibrium equation of arbitrary node i:
(2) security constraint of level Grid containing 110kV Optimal Operation Model
A, object function
For the sum of operating cost with all units in minimum system for target, the operating cost of unit includes startup-shutdown expense
With with fuel used to generate electricity consumption costs, such as formula (4):
In formula, T is total when hop count dispatching cycle, is divided into 96 periods by one day herein, then each period is 15min,
N1For conventional thermal power unit sum, CiU,tAnd CiD,tBootings and idleness expense of the respectively conventional thermal power unit i in period t, Fi,t
For conventional thermal power unit i period t power generation expense;N2For pumped storage unit sum, CsU,tAnd CsD,tRespectively pumped storage unit s exists
The booting of period t and idleness expense, since the operation of pumped storage unit will not consume fuel, so power generation expense is zero.
For conventional thermal power unit, generate electricity expense such as formula (5):
Fi,t=Ai×Pi,t (5)
In formula, AiThe fuel consumption characteristic coefficient for the expense that generates electricity for unit i, Pi,tFor unit i period t output.
B, constraints
1) the startup-shutdown expense restriction of unit
Be switched on expense restriction:
Idleness expense constrains:
In formula, KiAnd JiThe single of respectively conventional power unit i is switched on and idleness expense, KsAnd JsRespectively pumped storage unit s's
Single is switched on and idleness expense.Ii,t/Zs,tIt is conventional power unit i/ pumped storage unit s in the start and stop state of period t, starting state value is
1, shutdown status value is 0.
2) system power Constraints of Equilibrium
In formula, PLoad,tFor the system total load of period t, PLoss,tFor the system losses of period t, by being taken as total load work(
The certain percentage of rate is come approximate;Ppg,s,tAnd Ppp,s,tGenerated outputs and draw water power of the respectively pumped storage unit s in period t,
Generated output is just that the power that draws water is negative.
3) the conventional power unit output upper limit constrains
Ii,tPi,min≤Pi,t≤Ii,tPi,max (9)
In formula, Pi,minFor the minimum load of conventional power unit i, Pi,maxFor the maximum output of conventional power unit i.
4) unit climbing/landslide constraint
Consider that unit is no more than the minimum load of unit in first period of booting or the last one period of shutdown, then
Unit climbing/landslide constraint can indicate as follows:
In formula, ruiAnd rdiRespectively the climbing rate of unit i and landslide rate, T15For a run the period 15min.
5) unit minimum startup-shutdown time-constrain
Minimum available machine time constraint:
Minimum downtime constraint:
In formula, Ui/DiIt indicates the period that unit i starts to be switched on/shut down in dispatching cycle, is by last dispatching cycle
At the end of unit state determine, Ton_i/Toff_iFor minimum available machine time/minimum downtime of unit i, Xon_i,0/
Xoff_i,0For unit i when dispatching cycle starts continuous time of booting/shutdown.
6) the operation constraint of pump-storage generator
Output bound constrains:
In formula, Ppg,s,maxAnd Ppp,s,maxThe respectively maximum power generation and maximum of pumped storage unit s is drawn water power, Zpg,s,t
And Zpp,s,tGenerating states and draw water state of the respectively pumped storage unit s in period t, value are that 1 expression is in corresponding states, are taken
Value is that 0 expression is not at corresponding states.
Operating condition Constraints, i.e., same period pumped storage unit can not possibly be in and draw water and generating operation mode simultaneously, as follows:
Zs,t=Zpg,s,t+Zpp,s,t≤1 (14)
In actual operation, it is also necessary to meet balance of daily electric energy constraint, it is as follows:
In formula, ξ is the transfer efficiency of pumped storage unit, usually takes 75%.
In order to extend the service life of pumped storage unit, in actual motion the switching of two kinds of operating conditions of pumped storage unit need satisfaction one
Fixed time restriction, the present embodiment define switching time as half an hour, that is, need the switching time of 2 periods, it is therefore desirable to full
It is enough lower constraint:
7) system spinning reserve constrains
By reserved enough system spinning reserve capacities with cope with load prediction error band come influence.Utilize positive rotation
Spare capacity compensates the influence brought because underestimating system loading, is compensated using negative spinning reserve capacity and is brought because over-evaluating system loading
Influence.
Conventional power unit spinning reserve constrains:
In formula, sui,tAnd sdi,tThe positive and negative spinning reserve capacity that respectively unit i is capable of providing in period t, T10For unit
The spinning reserve response time, the present embodiment takes 10min.
Pumped storage unit spinning reserve constrains:
The spinning reserve constraint of demand that the spinning reserve of cumulative all units obtains system is as follows:
In formula, Su,tAnd Sd,tRespectively system is in the positive and negative spinning reserve capacity of period t, Lu% and Ld% is respectively load
Service demand factor of the prediction deviation to the positive and negative spinning reserve capacity of system.
8) node injects active power balance constraint
In formula,It is injection active power of the node i in period t;Pij,tFor the active power of period t circuit ij;θij,t
For the phase difference of voltage between period t node i and node j.
9) Network Security Constraints
According to DC flow model is improved, line transmission power constraint is as follows:
In formula, l is the circuit for needing to carry out safety check, Pl,km,tFor the transimission power of period t circuit l (k, m);;θkm,t
For the phase difference of voltage between period t node k and node m;For the transimission power limit value of circuit l (k, m).
Number of lines is very more in practical bulk power grid, if by the transimission power security constraint of all circuits in all periods
It is directly placed into model and is solved, the scale of model can be caused excessive, solve and take very much.In view of the operation of actual electric network
Situation, many circuits are not out-of-limit, therefore many line power security constraints are inoperative.In order to improve the solution of model
Efficiency is solved using the method for " verification-addition-is verified-added again again ":First to being added without line power security constraint
Model solved, obtain Unit Commitment and the optimum results of output;Then network security verification is carried out to this group solution, not
It is added in model by the line power constraint of safety check, solves again and once obtain optimum results;If being not met by
All safety checks will be then added in model again by the constraint of the line power of safety check and be solved again, verified, until
Until obtained optimum results meet all safety checks.
It can be seen that the security constraint Optimal Operation Model of level Grid containing 110kV of above-mentioned formula (4)~(21) description is mixed
Integral nonlinear program-ming model is closed, DICOPT the or SBB solvers that may be used in GAMS softwares carry out solution
The validity of this method is verified with reference to a sample calculation analysis:
Sample calculation analysis is carried out by taking IEEE9 node systems as an example, the grid structure of system is as shown in Figure 1.The load of the system
Number of nodes is 3, and transformer branch number is 3, and circuit number is 6, and 6 line parameter circuit values are shown in Table 1, it can be seen that its reactance and electricity
The ratio of resistance and practical 110kV circuits are similar;There are 3 generating sets, wherein 1 pump-storage generator is connected on BUS-1, and 2
Platform conventional thermal power unit is connected on BUS-2 and BUS-3 respectively, and the characterisitic parameter of 3 generating sets is shown in Table 2.Total load prediction curve
As shown in Fig. 2, maximum power is 385MW, minimum power 226.1875MW, due to the specific load power prediction of each node
Value can not obtain, so pressing percentage shared by each node load in BPA operation datas using by each moment total load predicted value
Than come the mode distributed, obtaining the predicted load of each node of each moment, system reference power is 100MW.It is respectively adopted
Improved DC flow model carries out Load flow calculation with DC flow model, and is compared relatively accurately with AC power flow result
Property;DC flow model will accurately be improved it is applied to power system security constraints and will ask in generation optimization scheduling problem a few days ago
The start and stop of unit a few days ago and the plan of contributing are formulated in solution, consideration 96 periods of whole day, optimization.Used computer is Intel (R)
Xeon (R) CPU E3-1270v5@3.60GHz, 32GB memories.
1 line parameter circuit value of table
The characterisitic parameter of 2 unit of table
It is calculated using the active power for improving DC flow model and DC flow model calculating circuit, and with AC power flow
As a result it is compared, the results are shown in Table 3.The excessive situation of the relative error brought in order to avoid the too small branch of transmitted power,
The branch rejecting that AC power flow line power is less than to circuit maximum power value 10% in network counts.
The different tide model comparison of computational results of table 3
From table 3 it is observed that the line power relative error that DC flow model is calculated has more than 20%, accidentally
Difference is larger, and it is more accurate relative to DC flow model to improve the power flow solutions that DC flow model is calculated, each circuit
Power relative error is both less than 5%, it is therefore proposed that improvement DC flow model be more suitable for calculate voltage class containing 110kV
The circuit active power of power grid, accuracy are higher.
In the Optimal Operation Model of foundation, the positive and negative spinning reserve service demand factor L of load prediction erroru% takes 3%,
Ld% takes 1%, and the percentage that network loss accounts for total load power is taken as 2.2%.Using in GAMS softwares CPLEX and DICOPT ask
Solution device solves the security constraint Optimal Operation Model based on DC flow model respectively, and (model is mixed integer linear programming mould
Type MILP) and the present embodiment establish based on improve DC flow model security constraint Optimal Operation Model (model is mixed
Close Integral nonlinear program-ming model M INLP).The target function value and calculate time comparison such as 4 institute of table that two kinds of model solutions obtain
Show.It can be seen that the solving speed of two kinds of models is quickly, based on the security constraint Optimized Operation mould for improving DC flow model
The target function value that type solves is slightly below the security constraint Optimal Operation Model based on DC flow model, straight due to improving
It is more accurate to flow the line power that tide model calculates, therefore the Optimized Operation result obtained by MINLP models is also more in line with
Actual conditions.
4 two kinds of model solution Comparative results of table
In the solution procedure for being based respectively on two kinds of models, the out-of-limit counting of line power becomes with the variation of iterations
Gesture is as shown in figure 3, out-of-limit in different periods for same circuit, then counting is cumulative.It can be seen that based on DC power flow mould
In the solution procedure of the security constraint Optimal Operation Model of type, line power security constraint has " verification-addition " 4 times altogether;And
In the solution procedure based on the security constraint Optimal Operation Model for improving DC flow model, line power security constraint has altogether
" verification-addition " 3 times.This is because the raising of line power computational accuracy so that the safety check number of needs is reduced.
The unit output Plan Curve comparison solved based on two kinds of tide models is as shown in Figure 4.It can be seen that
The fluctuation for the power generation dispatching Plan Curve that MINLP model solutions obtain is smaller, more meets practical operation situation.
It follows that the security constraint Optimization Scheduling of level Grid containing 110kV provided in this embodiment and the prior art
Compared to having the following advantages:
1) DC flow model is improved so that the result of calculation of line power is more in level Grid containing 110kV
Accurately.
2) the higher improvement DC flow model of accuracy is applied to the power system security constraints of voltage class containing 110kV a few days ago
It is solved in optimal hydropower operation model, solving speed is quickly and solving result is more rationally accurate.
It will be apparent to those skilled in the art that technical solution that can be as described above and design, make various other
Corresponding change and deformation, and all these changes and deformation should all belong to the protection domain of the claims in the present invention
Within.
Claims (10)
1. a kind of security constraint Optimization Scheduling of level Grid containing 110kV, which is characterized in that the method includes:
Ignore the ground branch in power grid, the active power for obtaining circuit under AC Ioad flow model is:
Pij=ViVj(Gijcosθij+Bijsinθij)-Vi 2Gij (1)
In formula:PijFor the active power of circuit ij;ViAnd VjThe respectively voltage magnitude of node i and node j;θijFor node i and
Phase difference of voltage between node j;Gij,BijThe respectively transconductance between node i and node j and mutual susceptance;
The active power of circuit under AC Ioad flow model is simplified, obtains improving direct current tide model, the improvement direct current tide mould
Type includes the active power equation of arbitrary circuit ij:
And the active power balance equation of injection of arbitrary node i:
Direct current tide model will be improved it is applied to the security constraint generation optimization scheduling problem a few days ago of the power grid of voltage class containing 110kV,
To obtain optimal dispatch result.
2. the security constraint Optimization Scheduling of level Grid containing 110kV as described in claim 1, which is characterized in that improving
Direct current tide model is applied to the security constraint of the power grid of voltage class containing 110kV a few days ago in generation optimization scheduling problem, Optimized Operation
Using the total generator operation cost minimization of power grid as object function, constraints includes the start-stop expense restriction of unit, system power
Constraints of Equilibrium, the constraint of the conventional power unit output upper limit, unit climbing/landslide constraint, unit minimum startup-shutdown time-constrain, draw water storage
The operation constraint of energy unit, the constraint of system spinning reserve, node inject active power balance constraint, Network Security Constraints.
3. the security constraint Optimization Scheduling of level Grid containing 110kV as claimed in claim 1 or 2, which is characterized in that handing over
The active power of circuit is simplified under stream tide model, is obtained improvement direct current tide model process and is:
1) due under normally operating condition, each node voltage amplitude is therefore close near rated value in grid power system
Seemingly think Vi≈1,Vj≈1;
2) voltage phase angle of two end node of circuit differs very little, i.e. θij≈ 0, therefore it is approximately considered sin (θij/2)≈θij/2,cos
(θij/2)≈1;
Therefore, formula (1) can be reduced to:
And then obtain improving the injecting power equilibrium equation of direct current tide model arbitrary node i:
4. the security constraint Optimization Scheduling of level Grid containing 110kV as claimed in claim 2, which is characterized in that the mesh
Scalar functions:
With the sum of operating cost for minimizing i.e. all units of the total generator operation cost of power grid for target, the operating cost packet of unit
Startup-shutdown expense and fuel used to generate electricity consumption costs are included, such as formula (4):
In formula, T is total when hop count dispatching cycle, is divided into 96 periods by one day herein, then each period is 15min, N1It is normal
Advise fired power generating unit sum, CiU,tAnd CiD,tBootings and idleness expense of the respectively conventional thermal power unit i in period t, Fi,tFor routine
Power generation expenses of the fired power generating unit i in period t;N2For pumped storage unit sum, CsU,tAnd CsD,tRespectively pumped storage unit s is period t's
Booting and idleness expense, since the operation of pumped storage unit will not consume fuel, so power generation expense is zero;
For conventional thermal power unit, generate electricity expense such as formula (5):
Fi,t=Ai×Pi,t
In formula, AiThe fuel consumption characteristic coefficient for the expense that generates electricity for unit i, Pi,tFor unit i period t output.
5. the security constraint Optimization Scheduling of level Grid containing 110kV as described in claim 2 or 4, which is characterized in that described
The startup-shutdown expense restriction of unit is:
Be switched on expense restriction:
Idleness expense constrains:
In formula, KiAnd JiThe single of respectively conventional power unit i is switched on and idleness expense, KsAnd JsThe respectively single of pumped storage unit s
Booting and idleness expense.Ii,t/Zs,tIt is start and stop states of the conventional power unit i/ pumped storage unit s in period t, starting state value is 1, is stopped
Machine state value is 0.
6. the security constraint Optimization Scheduling of level Grid containing 110kV as claimed in claim 2, which is characterized in that the system
The power-balance constraint of system is:
In formula, PLoad,tFor the system total load of period t, PLoss,tFor the system losses of period t, by being taken as total load power
Certain percentage is come approximate;Ppg,s,tAnd Ppp,s,tGenerated outputs and draw water power of the respectively pumped storage unit s in period t, power generation
Power is just that the power that draws water is negative;
The conventional power unit output upper limit is constrained to:
Ii,tPi,min≤Pi,t≤Ii,tPi,max (9)
In formula, Pi,minFor the minimum load of conventional power unit i, Pi,maxFor the maximum output of conventional power unit i.
7. the security constraint Optimization Scheduling of level Grid containing 110kV as claimed in claim 2, which is characterized in that the machine
Group climbing/landslide is constrained to:
Consider that unit is no more than the minimum load of unit in first period of booting or the last one period of shutdown, then unit
Climbing/landslide constraint representation is as follows:
In formula, ruiAnd rdiRespectively the climbing rate of unit i and landslide rate, T15For a run the period 15min.
The unit minimum startup-shutdown time-constrain is:
Minimum available machine time constraint:
Minimum downtime constraint:
In formula, Ui/DiIndicate the period that unit i starts to be switched on/shut down in dispatching cycle, was terminated by last dispatching cycle
When unit state determine, Ton_i/Toff_iFor minimum available machine time/minimum downtime of unit i, Xon_i,0/Xoff_i,0For
Unit i continuous times of booting/shutdown when dispatching cycle starts.
8. the security constraint Optimization Scheduling of level Grid containing 110kV as claimed in claim 2, which is characterized in that the pumping
The operation of water accumulation of energy unit is constrained to:
Output bound constrains:
In formula, Ppg,s,maxAnd Ppp,s,maxThe respectively maximum power generation and maximum of pumped storage unit s is drawn water power, Zpg,s,tWith
Zpp,s,tGenerating states and draw water state of the respectively pumped storage unit s in period t, value are that 1 expression is in corresponding states, value
It indicates to be not at corresponding states for 0;
Operating condition Constraints, i.e., same period pumped storage unit can not possibly be in and draw water and generating operation mode simultaneously, as follows:
Zs,t=Zpg,s,t+Zpp,s,t≤1 (14)
In actual operation, it is also necessary to meet balance of daily electric energy constraint, it is as follows:
In formula, ξ is the transfer efficiency of pumped storage unit, usually takes 75%.
In order to extend the service life of pumped storage unit, in actual motion the switching of two kinds of operating conditions of pumped storage unit need to meet time limit
System, definition switching time are half an hour, that is, need the switching time of 2 periods, it is therefore desirable to meet following constraint:
9. the security constraint Optimization Scheduling of level Grid containing 110kV as claimed in claim 2, which is characterized in that the system
System spinning reserve is constrained to:
By reserved enough system spinning reserve capacities with cope with load prediction error band come influence, it is spare using positive rotation
The influence that capacity compensation is brought because underestimating system loading compensates the shadow brought because over-evaluating system loading using negative spinning reserve capacity
It rings;
Conventional power unit spinning reserve constrains:
In formula, sui,tAnd sdi,tThe positive and negative spinning reserve capacity that respectively unit i is capable of providing in period t, T10For the rotation of unit
Turn the spare response time;
Pumped storage unit spinning reserve constrains:
The spinning reserve constraint of demand that the spinning reserve of cumulative all units obtains system is as follows:
In formula, Su,tAnd Sd,tRespectively system is in the positive and negative spinning reserve capacity of period t, Lu% and Ld% is respectively load prediction
Service demand factor of the deviation to the positive and negative spinning reserve capacity of system.
10. the security constraint Optimization Scheduling of level Grid containing 110kV as claimed in claim 2, which is characterized in that the section
Point injects active power balance constraint:
In formula,It is injection active power of the node i in period t;Pij,tFor the active power of period t circuit ij;θij,tFor when
Phase difference of voltage between section t node is and node j;
The Network Security Constraints
According to DC flow model is improved, line transmission power constraint is as follows:
In formula, l is the circuit for needing to carry out safety check, Pl,km,tFor the transimission power of period t circuit l (k, m);θkm,tFor the period
Phase difference of voltage between t nodes k and node m;For the transimission power limit value of circuit l (k, m).
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810345203.2A CN108695905B (en) | 2018-04-17 | 2018-04-17 | Safety constraint optimization scheduling method for 110 kV-level-contained power grid |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810345203.2A CN108695905B (en) | 2018-04-17 | 2018-04-17 | Safety constraint optimization scheduling method for 110 kV-level-contained power grid |
Publications (2)
Publication Number | Publication Date |
---|---|
CN108695905A true CN108695905A (en) | 2018-10-23 |
CN108695905B CN108695905B (en) | 2021-04-06 |
Family
ID=63845000
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201810345203.2A Active CN108695905B (en) | 2018-04-17 | 2018-04-17 | Safety constraint optimization scheduling method for 110 kV-level-contained power grid |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN108695905B (en) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111711185A (en) * | 2020-05-25 | 2020-09-25 | 国网青海省电力公司 | Day safety checking and blocking management method based on linearized alternating current power flow |
Family Cites Families (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101976836B (en) * | 2010-09-30 | 2013-03-06 | 河海大学 | Method for statically analyzing voltage stabilization of VSC-HVDC (Voltage-Sourced Converter-High Voltage Director Current) containing AC and DC system |
CN102623988B (en) * | 2012-03-27 | 2014-01-29 | 清华大学 | Automatic generating method for day-ahead plan power flow in power grid |
CN102983573B (en) * | 2012-11-09 | 2014-10-15 | 天津大学 | Security constraint economic dispatch method based on security domains |
EP3035476B1 (en) * | 2014-12-18 | 2019-07-24 | Siemens Aktiengesellschaft | Method, controller, and computer program product for controlling a voltage in a DC power grid |
CN105552915B (en) * | 2016-02-25 | 2018-02-06 | 国网福建省电力有限公司 | A kind of alternating current-direct current mixing network optimization tidal current computing methods of HVDC containing MMC for considering Static Security Constraints |
CN105976065A (en) * | 2016-05-20 | 2016-09-28 | 燕山大学 | Environment economic dispatching two-level dispatching solution method including centralized charging stations |
CN106385035B (en) * | 2016-12-05 | 2019-02-19 | 国网江西省电力公司经济技术研究院 | It is a kind of meter and charging capacitor DC power flow calculation method |
CN107069741A (en) * | 2017-04-20 | 2017-08-18 | 河海大学 | A kind of novel linear tidal current computing method |
CN106899022A (en) * | 2017-04-25 | 2017-06-27 | 大连理工大学 | A kind of out-of-limit control analysis method of interconnection |
-
2018
- 2018-04-17 CN CN201810345203.2A patent/CN108695905B/en active Active
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111711185A (en) * | 2020-05-25 | 2020-09-25 | 国网青海省电力公司 | Day safety checking and blocking management method based on linearized alternating current power flow |
CN111711185B (en) * | 2020-05-25 | 2021-11-19 | 国网青海省电力公司 | Day safety checking and blocking management method based on linearized alternating current power flow |
Also Published As
Publication number | Publication date |
---|---|
CN108695905B (en) | 2021-04-06 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN108110804B (en) | Power grid optimal scheduling method with pumped storage unit | |
KR20190073340A (en) | Intelligent grid operation system for managing distributed energy resources in the grid | |
Fu et al. | Research on joint optimal dispatching method for hybrid power system considering system security | |
JP2015080401A (en) | Methods and systems for controlling electric network | |
CN102074954A (en) | Comprehensive energy-conservation evaluation and decision-making method for urban and rural distribution network | |
CN105470985B (en) | A kind of wind stores up the flexible self-start method of Isolated Network System | |
CN107749635A (en) | A kind of closed loop network running optimizatin method | |
CN103187735A (en) | Bidirectional intelligent gateway device for distributed new energy grid connection | |
CN109494727A (en) | Consider the active and idle coordination optimization operation method of power distribution network of demand response | |
CN108551176A (en) | A kind of energy storage battery system capacity collocation method of combination energy storage balancing technique | |
CN106712103B (en) | A kind of micro-capacitance sensor voltage stabilizing control system and control method | |
CN106684931A (en) | Active control method and system of new-energy power plant | |
CN108695905A (en) | Safety constraint optimization scheduling method for 110 kV-level-contained power grid | |
CN103928944B (en) | Photovoltaic DC-to-AC converter power stage combined adjuster | |
CN105429172B (en) | The power generation dispatching method of strengthening system operational flexibility | |
CN107612018B (en) | The control method and device of inverter, inverter, computer readable storage medium | |
CN107086586B (en) | The reactive-load compensation method of wind farm group based on the idle generating ability of double-fed blower | |
CN113572191B (en) | Distributed camera setting configuration method and system for large-scale new energy direct current delivery system | |
CN115689187A (en) | Distributed photovoltaic open capacity calculation method and system for power distribution network | |
CN115912496A (en) | Power distribution network source-storage-load multi-mode power scheduling method | |
CN113131528A (en) | Method, device, equipment and storage medium for determining optimal capacity of wind fire bundling | |
CN110311947A (en) | A kind of new energy electricity consumption percent information method for pushing and system | |
Li et al. | Integrated Frequency-constrained Scheduling Considering Coordination of Frequency Regulation Capabilities from Multi-source Converters | |
CN105207249A (en) | Wind power plant constant power factor control system and control method thereof | |
CN117996861B (en) | Scheduling method for light-water coupling residual electricity hydrogen production of power distribution network and energy management device |
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 |