CN109167383A - Electric system peak regulation optimization method based on exact linearization method power network model - Google Patents
Electric system peak regulation optimization method based on exact linearization method power network model Download PDFInfo
- Publication number
- CN109167383A CN109167383A CN201810937487.4A CN201810937487A CN109167383A CN 109167383 A CN109167383 A CN 109167383A CN 201810937487 A CN201810937487 A CN 201810937487A CN 109167383 A CN109167383 A CN 109167383A
- Authority
- CN
- China
- Prior art keywords
- power
- unit
- peak regulation
- constraint
- peak
- 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
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
-
- 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
Landscapes
- Engineering & Computer Science (AREA)
- Power Engineering (AREA)
- Supply And Distribution Of Alternating Current (AREA)
- Feedback Control In General (AREA)
Abstract
The present invention relates to a kind of electric system peak regulation Optimization Scheduling based on exact linearization method power network model includes the following steps: S1: obtaining the conveying bound data of power system load demand data, each Generator Unit Operating Parameters and each route of cost coefficient, power grid;S2: according to acquired supplemental characteristic, the polymorphic type power supply combined adjusting peak Optimized model containing nuclear power is established;S3: carrying out exact linearization method processing to polymorphic type power supply combined adjusting peak Optimized model, obtains linearisation peak regulation scheduling optimization model;S4: initialization linearisation peak regulation scheduling optimization model parameters solve exact linearization method peak regulation Optimized model using branch and bound method, obtain unit output arrangement.The present invention comprehensively considers the power generation characteristics of different type power supply, establishes the peak-load regulating Optimized model for considering that nuclear power participates in peak regulation, provides Optimized Operation strategy for polymorphic type power supply in electric system.
Description
Technical field
The invention belongs to Economic Dispatchs to optimize field, and in particular to one kind is based on exact linearization method electric power networks
The electric system peak regulation optimization method of model.
Background technique
With the increase of China's power system load demand, power grid peak-valley difference is gradually expanded;On the other hand, all kinds of cleaning energy
Source installation scale accounting constantly increase, coastal area nuclear power installed capacity rapid growth, peak load regulation network face more challenges and
Pressure, the peak load regulation network optimizing research for carrying out the multiple power sources structure containing nuclear power are of great significance.
The reactor of the nuclear power unit in current China all has certain load-following capacity in design, in peak regulation technique
On be it is feasible, can satisfy peak load regulation network scheduling demand;On the other hand, nuclear power peak regulation bring fringe cost and its peak regulation
Participation mode and depth have direct relation.It runs and plans in layout process in dispatching of power netwoks, how to optimize nuclear power, wind-powered electricity generation, fire
The all types of power supply generation schedules such as electricity, water power realize that system is always sent out under the premise of meeting peak-load regulating demand and security constraint
The optimization of electric cost is Operation of Electric Systems and staff planners' urgent problem to be solved.
Summary of the invention
The object of the present invention is in the case where meeting the constraint of unit operational safety, with the minimum target letter of system total power production cost
Number, establishes Optimal Operation Model, reasonable arrangement unit output.
To achieve the above object, the technical solution adopted by the present invention is that:
A kind of combined adjusting peak Optimized Operation for considering nuclear power and participating in, it is characterized in that considering thermoelectricity fuel cost, thermoelectricity tune
Peak cost, thermoelectricity start-up and shut-down costs, nuclear power peak regulation cost and nuclear power peaking operation characteristic, Line Flow constraint, unit climbing speed
The various constraint conditions such as rate constraint, unit minimum start and stop constraint, system power Constraints of Equilibrium, establish combined adjusting peak model, to mould
Part Nonlinear Cost Function, electric power networks Nonlinear Constraints in type carry out linearization process, specifically include in following
Hold:
Step S1: with 24 hours for dispatching cycle, power system load demand data, each Generator Unit Operating Parameters are obtained
And the conveying constraint of each route of cost coefficient, power grid.
Step S2: according to acquired parameter, the polymorphic type power supply combined adjusting peak Optimized model containing nuclear power is established;
Step S3: carrying out exact linearization method processing to model, obtains linearisation peak regulation scheduling optimization model;
Step S4: initialization model parameters solve exact linearization method peak regulation Optimized model using branch and bound method, obtain
Take unit output arrangement.
Further, in step sl, acquired power system load demand data refer specifically to be in following 24 hours
It unites predicted load, is divided into 15 minutes between predicted time;Acquired Generator Unit Operating Parameters specifically include: nuclear power peak regulation fortune
Row characteristic, all kinds of unit ramping rate constraints such as water power, thermoelectricity, nuclear power, combustion gas, pumped storage, all kinds of unit minimum start and stop constraints;Hair
This coefficient of motor form specifically includes: thermoelectricity fuel cost, thermoelectricity peak regulation cost, thermoelectricity start-up and shut-down costs, nuclear power peak regulation cost;Institute
Each route conveying constraint of the power grid of acquisition includes the active power constraint of each route in power grid, transformer.
Further, in step s 2, the polymorphic type power supply combined adjusting peak Optimized model containing nuclear power established, is one
Nonlinear mixed-integer programming model, form are as follows:
Wherein x1、x2For the optimized variable of optimal load flow mathematical model, dimension is respectively n1、n2, wherein x1For continuous variable,
x2For discrete integer variable, f (x1,x2) be optimal load flow mathematical model objective function, H (x1,x2) it is equality constraint intersection, dimension
Degree is m, G (x1,x2) it is inequality constraints intersection, dimension s, Gmax、GminThe respectively upper lower limit value of inequality constraints.
Polymorphic type power supply combined adjusting peak Optimized model proposed by the present invention containing nuclear power, the table of objective function OF in (1) formula
Up to formula are as follows:
Wherein, T=24h is a dispatching cycle, m0、m1Respectively thermoelectricity, nuclear power unit number;For thermoelectricity
Fuel cost, whereinai、bi、ciFor cost coefficient, PG,it、PN,itFor thermoelectricity, nuclear power unit i
In the activity of force to be found out of t moment;For the opening of fired power generating unit, stop cost, wherein
STCG,it=stiyit, SDCG,it=sdizit, sti、sdiThe respectively starting of unit i, shutdown cost coefficient, yit、zitExpression machine
Beginning on/off state of the group i in t moment;For thermoelectricity peak regulation cost,
For nuclear power peak regulation cost, PN,iMaxFor the maximum output power of nuclear power unit i, CG、CNRespectively thermoelectricity depth peak regulation and nuclear power tune
Peak cost coefficient, for the nuclear power unit of CPR1000 heap-type, CN=71 yuan/MWh.
Polymorphic type power supply combined adjusting peak Optimized model proposed by the present invention containing nuclear power, optimized variable includes fire in (1) formula
Electricity, to be found out activity of force of the nuclear power unit i in t moment, startup-shutdown state.
Polymorphic type power supply combined adjusting peak Optimized model proposed by the present invention containing nuclear power, constraint condition includes: in (1) formula
S201: power-balance constraint
The sum of thermoelectricity, the nuclear power unit power output of t period whole, P are indicated in formula on the left of equationL,tFor the load of system t period
Demand.
S202: spinning reserve constraint
P in formulart, P 'rtRespectively positive and negative stand-by requirement of the system in t moment, P 'G,iMinIt is that fired power generating unit i is participated in substantially
Power output lower limit of the power when peak regulation, (4) formula mainly consider to be provided by fired power generating unit spare.
S203: the operation constraint of nuclear power
The power producing characteristics curve of nuclear power meets the peak regulation mode of " 12-3-6-3 ", i.e. 12 hours Operation at full power, and then 3
Hour downrating is kept for low power run 6 hours, then 3 hours power per liter are run to full power shape to low power state
State completes a peak regulation period.
Service capacity curve is expressed as:
PN,it=eitPN,iMax+fitPN,iMin+git(PN,iMin+ΔPN,i)+hit(PN,iMin+2ΔPN,i) (5)
In formula: eit, fit, git, hitIt indicates power operating states mark, is { 0,1 } variable;Ti e, Ti fIndicate most Grain Full function
Rate operation time and minimum low power run time value, PN,iMin、PN,iMaxIndicate minimum, the maximum output function of nuclear power unit i
Rate, Δ PN,iIndicate the power variation of nuclear power unit hour.
S204: the depth peak regulation constraint of nuclear power
ηit=(PN,iMax-PN,it)/PN,iMax≤ηmax (7)
The type of reactor of nuclear power unit, unit capacity are different, have different peak regulation depth threshold ηmax。
S205: fired power generating unit depth peak regulation constraint
Unit output bound constraint: PG,iMin≤PG,it≤PG,iMax
(8)
It needs to meet when fired power generating unit depth peak regulation:
P′G,iMin≤PG,it≤PG,iMax (9)
The peak regulation state of thermoelectricity is divided into basic peak regulation and depth peak regulation, only reach certain peak regulation depth just give peak regulation at
This compensation, P 'G,iMinIndicate power output lower limit when thermoelectricity participates in basic peak regulation, i.e., known paid peak regulation threshold value, PG,itIndicate to
Power output of the fired power generating unit i asked in t moment, PG,iMin、PG,iMaxFor known fired power generating unit i power output bound.
S206: Unit Ramp Rate constraint | PG,it-PG,it-1|≤Δi;Wherein, ΔiIndicate creep speed.
S207: Line Flow constraint
For a n node system, the node power flow equation of AC power flow is
Wherein, subscript i, j respectively indicates node i, j, Pi、QiThe respectively active power and reactive power of node i, Gij、
BijTransconductance and mutual susceptance respectively between node i and node j, Vi、VjIndicate node voltage amplitude, θijPhase angle between expression node
Difference.
Further, exact linearization method processing is carried out to model in step s3, obtains linearisation peak regulation optimizing scheduling mould
Type specifically includes following content:
S301: the piece-wise linearization of thermoelectricity fuel cost function, by fired power generating unit fuel cost argument of function unit
Changed power section [PG,iMin,PG,iMax] n equal part, nonlinear fuel cost function linearize in each equal segments close
Seemingly, principle is as shown in Figure 3
For the linear segmented slope of fuel cost,For power PG,itSegmentation variable, n is specified piecewise interval
Number, k be piecewise interval relative to call number, i.e., constant 1 arrive n.
S302: Unit Commitment and Climing constant linearization process
Recursive non-linear feature is presented in the constraint condition of minimum start-stop time, if directly there are difficulty for Optimization Solution, needs
This complicated nonlinear complementary problem is converted into simple linear restriction of equal value, then is solved.
Linearization process: uitIndicate that unit i indicates that unit is in operating status in the operating status of t moment, 1,0 indicates machine
Group is in shutdown status;yit、zitIndicate beginning on/off state of the unit i in t moment, yit=1 expression unit, which is in, to be opened
Beginning starting state, zit=1, which indicates that unit is in, starts shutdown status.
About above-mentioned constraint, it is described as follows:
1, the power of unit i remains less than power capacity, i.e.,
If 2, unit is shut down at next hour (t+1),Because of PG,it+1=0, PG,itNo more than
SDiValue;
If 3, unit is in (u in operating status of previous hourit-1=1), and continue to remain operational state, then
PG,itCompared to PG,it-1Growth can not be more than RUi.I.e.
If 4, unit is in close state (u in previous hourit-1=0), and in t moment start starting operation yit
=1, then PG,itIt can not be more than SUi, i.e.,
To sum up situation constitutes part (b) in constraint (13), and the part (c) similarly can be obtained.
The minimum of unit i stops machine time-constrain
The minimum runing time of unit i constrains
The minimum downtime of unit i constrains
In formula, RUi、RDiRise for the climbing under unit i normal operating condition and the lower rate of deceleration limits;SUi、SDiFor machine
Organize fall off rate limitation when climbing speed and the closing when i starting;UTi、DTiFor the minimum operation of unit i and downtime;It has initially been run for unit, downtime, uIt=0For unit initial operating state,P it、It indicates to consider creep speed
With the unit minimax operational limit after the constraint of start and stop operating status;ζiIt is that runing time and minimum are continuous about unit
Coupled relation between runing time, ζi=min { T, (UTi-Ui 0)uIt=0, ξiIt is about unit downtime and minimum
Coupled relation between continuous downtime, ξi=min { T, (DTi-Si 0)(1-uIt=0)}。
S303: trend constraint linearisation, in Unit Commitment and economic load dispatching, general use ignores network loss
DC power flow constraint, has the characteristics that linearisation and rapidity, but only account for the relationship of trend and voltage phase angle, accuracy is not
The description of height, especially voltage and reactive power.
Different from conventional DC flow model, reactive power and voltage constraint are considered here, is obtained using linearization approximate
Node voltage meets V ≈ 1.0p.u. to be had according to the expansion of reciprocal function: 1/V ≈ 2-V, substitutes into above formula equal sign
The linearisation of node active balance equation can be obtained in left side are as follows:
B′ijWith BijDifference be to be eliminated in element from susceptance bii;Similarly, node reactive power equilibrium equation linearizes
For
Linearisation Branch Power Flow equation can similarly be obtained
In step s 4, exact linearization method peak regulation Optimized model is solved using branch and bound method, obtains unit output arrangement,
Specifically include following content: initialization model parameter, including optimized variable initial value, objective function and constraint equation coefficient value,
Constraint condition upper limit value and lower limit value call the CPLEX solver based on branch and bound method to solve linearisation peak regulation optimization mould
Type obtains the power output arrangement of each unit.
Compared to the prior art, the beneficial effects of the present invention are:
1) present invention comprehensively considers the power generation characteristics of the different types power supply such as nuclear power, thermoelectricity, water power, pneumoelectric, pumped storage, establishes
Consider that nuclear power participates in the peak-load regulating Optimized model of peak regulation, provides Optimized Operation strategy for polymorphic type power supply in electric system;
2) present invention sufficiently takes into account model solvability and required precision, carries out to non-linear hybrid integration plan model accurate
Linearization process, significantly reduces the solution difficulty of model, while having ensured the optimality of model optimization result.
Detailed description of the invention
Fig. 1 is structure flow chart of the invention;
Fig. 2 is that nuclear power unit participates in peak regulation " 12-3-6-3 " operational mode;Abscissa represents whole day 24 hours in figure, indulges
Coordinate value 100% represents Operation at full power, and ordinate value 70% represents nuclear power unit and is in 70% Operation at full power state;
Fig. 3 is the piece-wise linearization schematic diagram of fuel cost function of the present invention;Abscissa represents unit in figure
Power is run, ordinate represents the fuel cost of unit, and abscissa unit is run power n equal part, and the combustion to each equal segments
Material curve is linearized.
Specific embodiment
With reference to the accompanying drawing, technical solution of the present invention is specifically described.
As shown in Figure 1, a kind of electric system peak regulation Optimization Scheduling based on exact linearization method power network model, packet
Include following steps:
Step S1: with 24 hours for dispatching cycle, power system load demand data, each Generator Unit Operating Parameters are obtained
And the conveying constraint of each route of cost coefficient, power grid.
Acquired power system load demand data refers specifically to system loading predicted value in 24 hours following, predicted time
Between be divided into 15 minutes;Acquired Generator Unit Operating Parameters specifically include: nuclear power peaking operation characteristic, water power, thermoelectricity, core
All kinds of unit ramping rate constraints such as electricity, combustion gas, pumped storage, all kinds of unit minimum start and stop constraints;Generating set cost coefficient is specific
It include: thermoelectricity fuel cost, thermoelectricity peak regulation cost, thermoelectricity start-up and shut-down costs, nuclear power peak regulation cost;Each route of acquired power grid is defeated
Sending constraint includes the active power constraint of each route in power grid, transformer.
Step S2: according to acquired parameter, the polymorphic type power supply combined adjusting peak Optimized model containing nuclear power is established;
The polymorphic type power supply combined adjusting peak Optimized model containing nuclear power established, is a nonlinear mixed-integer programming mould
Type, form are as follows:
Wherein x1、x2For the optimized variable of optimal load flow mathematical model, dimension is respectively n1、n2, wherein x1For continuous variable,
x2For discrete integer variable, f (x1,x2) be optimal load flow mathematical model objective function, H (x1,x2) it is equality constraint intersection, dimension
Degree is m, G (x1,x2) it is inequality constraints intersection, dimension s, Gmax、GminThe respectively upper lower limit value of inequality constraints.
Polymorphic type power supply combined adjusting peak Optimized model proposed by the present invention containing nuclear power, the table of objective function OF in (1) formula
Up to formula are as follows:
Wherein, T=24h is a dispatching cycle, m0、m1Respectively thermoelectricity, nuclear power unit number;For thermoelectricity combustion
Expect cost, whereinai、bi、ciFor cost coefficient, PG,it、PN,itIt is thermoelectricity, nuclear power unit i in t
The activity of force to be found out at moment;For the opening of fired power generating unit, stop cost, wherein STCG,it
=stiyit, SDCG,it=sdizit, sti、sdiThe respectively starting of unit i, shutdown cost coefficient, yit、zitIndicate unit i in t
The beginning on/off state at moment;For thermoelectricity peak regulation cost,For nuclear power
Peak regulation cost, PN,iMaxFor the maximum output power of nuclear power unit i, CG、CNRespectively thermoelectricity depth peak regulation and nuclear power peak regulation cost
Coefficient, for the nuclear power unit of CPR1000 heap-type, CN=71 yuan/MWh.
Polymorphic type power supply combined adjusting peak Optimized model proposed by the present invention containing nuclear power, optimized variable includes fire in (1) formula
Electricity, to be found out activity of force of the nuclear power unit i in t moment, startup-shutdown state.
Polymorphic type power supply combined adjusting peak Optimized model proposed by the present invention containing nuclear power, constraint condition includes: in (1) formula
S201: power-balance constraint
The sum of thermoelectricity, the nuclear power unit power output of t period whole, P are indicated in formula on the left of equationL,tFor the load of system t period
Demand.
S202: spinning reserve constraint
P in formulart, P 'rtRespectively positive and negative stand-by requirement of the system in t moment, P 'G,iMinIt is that fired power generating unit i is participated in substantially
Power output lower limit of the power when peak regulation, (4) formula mainly consider to be provided by fired power generating unit spare.
S203: the operation constraint of nuclear power
The power producing characteristics curve of nuclear power meets the peak regulation mode of " 12-3-6-3 ", i.e. 12 hours Operation at full power, and then 3
Hour downrating is kept for low power run 6 hours, then 3 hours power per liter are run to full power shape to low power state
State completes a peak regulation period.
Service capacity curve is expressed as:
PN,it=eitPN,iMax+fitPN,iMin+git(PN,iMin+ΔPN,i)+hit(PN,iMin+2ΔPN,i) (5)
In formula: eit, fit, git, hitIt indicates power operating states mark, is { 0,1 } variable;Ti e, Ti fIndicate most Grain Full function
Rate operation time and minimum low power run time value, PN,iMin、PN,iMaxIndicate minimum, the maximum output function of nuclear power unit i
Rate, Δ PN,iIndicate the power variation of nuclear power unit hour.
S204: the depth peak regulation constraint of nuclear power
ηit=(PN,iMax-PN,it)/PN,iMax≤ηmax (7)
The type of reactor of nuclear power unit, unit capacity are different, have different peak regulation depth threshold ηmax。
S205: fired power generating unit depth peak regulation constraint
Unit output bound constraint: PG,iMin≤PG,it≤PG,iMax
(8)
It needs to meet when fired power generating unit depth peak regulation:
P′G,iMin≤PG,it≤PG,iMax (9)
The peak regulation state of thermoelectricity is divided into basic peak regulation and depth peak regulation, only reach certain peak regulation depth just give peak regulation at
This compensation, P 'G,iMinIndicate power output lower limit when thermoelectricity participates in basic peak regulation, i.e., known paid peak regulation threshold value, PG,itIndicate to
Power output of the fired power generating unit i asked in t moment, PG,iMin、PG,iMaxFor known fired power generating unit i power output bound.
S206: Unit Ramp Rate constraint | PG,it-PG,it-1|≤Δi;Wherein, ΔiIndicate creep speed.
S207: Line Flow constraint
For a n node system, the node power flow equation of AC power flow is
Wherein, subscript i, j respectively indicates node i, j, Pi、QiThe respectively active power and reactive power of node i, Gij、
BijTransconductance and mutual susceptance respectively between node i and node j, Vi、VjIndicate node voltage amplitude, θijPhase angle between expression node
Difference.
Step S3: carrying out exact linearization method processing to model, obtains linearisation peak regulation scheduling optimization model;
Specifically include following content:
S301: the piece-wise linearization of thermoelectricity fuel cost function, by fired power generating unit fuel cost argument of function unit
Changed power section [PG,iMin,PG,iMax] n equal part, nonlinear fuel cost function linearize in each equal segments close
Seemingly, principle is as shown in Figure 3
For the linear segmented slope of fuel cost,For power PG,itSegmentation variable, n is specified piecewise interval
Number, k be piecewise interval relative to call number, i.e., constant 1 arrive n.
S302: Unit Commitment and Climing constant linearization process
Recursive non-linear feature is presented in the constraint condition of minimum start-stop time, if directly there are difficulty for Optimization Solution, needs
This complicated nonlinear complementary problem is converted into simple linear restriction of equal value, then is solved.
Linearization process: uitIndicate that unit i indicates that unit is in operating status in the operating status of t moment, 1,0 indicates machine
Group is in shutdown status;yit、zitIndicate beginning on/off state of the unit i in t moment, yit=1 expression unit, which is in, to be opened
Beginning starting state, zit=1, which indicates that unit is in, starts shutdown status.
About above-mentioned constraint, it is described as follows:
1, the power of unit i remains less than power capacity, i.e.,
If 2, unit is shut down at next hour (t+1),Because of PG,it+1=0, PG,itNo more than
SDiValue;
If 3, unit is in (u in operating status of previous hourit-1=1), and continue to remain operational state, then
PG,itCompared to PG,it-1Growth can not be more than RUi.I.e.
If 4, unit is in close state (u in previous hourit-1=0), and in t moment start starting operation yit
=1, then PG,itIt can not be more than SUi, i.e.,
To sum up situation constitutes part (b) in constraint (13), and the part (c) similarly can be obtained.
The minimum of unit i stops machine time-constrain
The minimum runing time of unit i constrains
The minimum downtime of unit i constrains
In formula, RUi、RDiRise for the climbing under unit i normal operating condition and the lower rate of deceleration limits;SUi、SDiFor machine
Organize fall off rate limitation when climbing speed and the closing when i starting;UTi、DTiFor the minimum operation of unit i and downtime;It has initially been run for unit, downtime, uIt=0For unit initial operating state,P it、It indicates to consider creep speed
With the unit minimax operational limit after the constraint of start and stop operating status;ζiIt is that runing time and minimum are continuous about unit
Coupled relation between runing time, ζi=min { T, (UTi-Ui 0)uIt=0, ξiIt is about unit downtime and minimum
Coupled relation between continuous downtime, ξi=min { T, (DTi-Si 0)(1-uIt=0)}。
S303: trend constraint linearisation, in Unit Commitment and economic load dispatching, general use ignores network loss
DC power flow constraint, has the characteristics that linearisation and rapidity, but only account for the relationship of trend and voltage phase angle, accuracy is not
The description of height, especially voltage and reactive power.
Different from conventional DC flow model, reactive power and voltage constraint are considered here, is obtained using linearization approximate
Node voltage meets V ≈ 1.0p.u. to be had according to the expansion of reciprocal function: 1/V ≈ 2-V, substitutes into above formula equal sign
The linearisation of node active balance equation can be obtained in left side are as follows:
B′ijWith BijDifference be to be eliminated in element from susceptance bii;Similarly, node reactive power equilibrium equation linearizes
For
Linearisation Branch Power Flow equation can similarly be obtained
Step S4: initialization model parameter, including optimized variable initial value, objective function and constraint equation coefficient value, about
Beam condition upper limit value and lower limit value call the CPLEX solver based on branch and bound method to solve the linearisation peak regulation Optimized model,
Obtain the power output arrangement of each unit.
The above are preferred embodiments of the present invention, all any changes made according to the technical solution of the present invention, and generated function is made
When with range without departing from technical solution of the present invention, all belong to the scope of protection of the present invention.
Claims (7)
1. a kind of electric system peak regulation Optimization Scheduling based on exact linearization method power network model, which is characterized in that packet
Include following steps:
Step S1: power system load demand data, each Generator Unit Operating Parameters and each route of cost coefficient, power grid are obtained
Convey bound data;
Step S2: according to acquired supplemental characteristic, the polymorphic type power supply combined adjusting peak Optimized model containing nuclear power is established;
Step S3: carrying out exact linearization method processing to polymorphic type power supply combined adjusting peak Optimized model, obtains linearisation peak regulation scheduling
Optimized model;
Step S4: initialization linearisation peak regulation scheduling optimization model parameters solve exact linearization method using branch and bound method
Peak regulation Optimized model obtains unit output arrangement.
2. a kind of linear tidal current computing method of three-phase polar coordinate system applied to active distribution network according to claim 1,
It is characterized by: acquired power system load demand data refers specifically to system loading predicted value in 24 hours following, prediction
Time interval is 15 minutes;The Generator Unit Operating Parameters specifically include: nuclear power peaking operation characteristic, water power, thermoelectricity, core
Electricity, combustion gas, pumped storage unit ramping rate constraints, unit minimum start and stop constraint;The generating set cost coefficient specifically includes: fire
Electric fuel cost, thermoelectricity peak regulation cost, thermoelectricity start-up and shut-down costs, nuclear power peak regulation cost;Each route conveying of power grid, which constrains, includes
The active power constraint of route, transformer in power grid.
3. a kind of linear tidal current computing method of three-phase polar coordinate system applied to active distribution network according to claim 1,
It is characterized by: in step S2 specifically: the polymorphic type power supply combined adjusting peak Optimized model containing nuclear power established, is one non-
Linear mixed integer programing model, form are as follows:
Wherein x1、x2For the optimized variable of optimal load flow mathematical model, dimension is respectively n1、n2, wherein x1For continuous variable, x2For
Discrete integer variable, f (x1,x2) be optimal load flow mathematical model objective function, H (x1,x2) it is equality constraint intersection, dimension
For m, G (x1,x2) it is inequality constraints intersection, dimension s, Gmax、GminThe respectively upper lower limit value of inequality constraints.
4. the polymorphic type power supply combined adjusting peak Optimized model containing nuclear power according to claim 3, it is characterised in that:
The expression formula of the objective function OF of the polymorphic type power supply combined adjusting peak Optimized model containing nuclear power are as follows:
Wherein, T=24h is a dispatching cycle, m0、m1Respectively thermoelectricity, nuclear power unit number;For thermoelectricity fuel at
This, whereinai、bi、ciFor cost coefficient, PG,it、PN,itIt is thermoelectricity, nuclear power unit i in t moment
Activity of force to be found out;For the opening of fired power generating unit, stop cost, wherein STCG,it=
stiyit, SDCG,it=sdizit, sti、sdiThe respectively starting of unit i, shutdown cost coefficient, yit、zitIndicate unit i in t
The beginning on/off state at quarter;For thermoelectricity peak regulation cost,For nuclear power tune
Peak cost, PN,iMaxFor the maximum output power of nuclear power unit i, CG、CNRespectively thermoelectricity depth peak regulation and nuclear power peak regulation cost system
Number, for the nuclear power unit of CPR1000 heap-type, CN=71 yuan/MWh;The optimized variable of model includes thermoelectricity, nuclear power unit i in t
The activity of force to be found out at moment, startup-shutdown state.
5. the polymorphic type power supply combined adjusting peak Optimized model containing nuclear power according to claim 3, it is characterised in that: described to contain core
The constraint condition of polymorphic type power supply combined adjusting peak Optimized model of electricity includes:
S501: power-balance constraint
The sum of thermoelectricity, the nuclear power unit power output of t period whole, P are indicated in formula on the left of equationL,tIt is needed for the load of system t period
It asks;
S502: spinning reserve constraint
P in formulart, P 'rtRespectively positive and negative stand-by requirement of the system in t moment, P 'G,iMinIt is that fired power generating unit i participates in basic peak regulation
When the power output lower limit of the power, (4) formula mainly considers to be provided by fired power generating unit spare;
S503: the operation constraint of nuclear power
The peak regulation mode of the power producing characteristics curve satisfaction " 12-3-6-3 " of nuclear power, i.e. 12 hours Operation at full power, then 3 hours
Downrating is kept for low power run 6 hours, then 3 hours power per liter are run to full power state, complete to low power state
At a peak regulation period;
Service capacity curve is expressed as:
PN,it=eitPN,iMax+fitPN,iMin+git(PN,iMin+ΔPN,i)+hit(PN,iMin+2ΔPN,i) (5)
In formula: eit, fit, git, hitIt indicates power operating states mark, is { 0,1 } variable;Ti e, Ti fIndicate minimum full power fortune
Row time value and minimum low power run time value, PN,iMin、PN,iMaxIndicate minimum, the maximum output power of nuclear power unit i, Δ
PN,iIndicate the power variation of nuclear power unit hour;
S504: the depth peak regulation constraint of nuclear power
ηit=(PN,iMax-PN,it)/PN,iMax≤ηmax (7)
The type of reactor of nuclear power unit, unit capacity are different, have different peak regulation depth threshold ηmax。
S505: fired power generating unit depth peak regulation constraint
Unit output bound constraint: PG,iMin≤PG,it≤PG,iMax (8)
It needs to meet when fired power generating unit depth peak regulation:
P′G,iMin≤PG,it≤PG,iMax (9)
The peak regulation state of thermoelectricity is divided into basic peak regulation and depth peak regulation, only reaches certain peak regulation depth and just gives peak regulation cost benefit
It repays, P 'G,iMinIndicate power output lower limit when thermoelectricity participates in basic peak regulation, i.e., known paid peak regulation threshold value, PG,itIt indicates to be asked
Power output of the fired power generating unit i in t moment, PG,iMin、PG,iMaxFor known fired power generating unit i power output bound;
S506: Unit Ramp Rate constraint | PG,it-PG,it-1|≤Δi;Wherein, ΔiIndicate creep speed;
S507: Line Flow constraint
For a n node system, the node power flow equation of AC power flow is
Wherein, subscript i, j respectively indicates node i, j, Pi、QiThe respectively active power and reactive power of node i, Gij、BijPoint
Transconductance and mutual susceptance not between node i and node j, Vi、VjIndicate node voltage amplitude, θijPhase angle difference between expression node.
6. a kind of linear tidal current computing method of three-phase polar coordinate system applied to active distribution network according to claim 1,
It is characterized by: the step S3 specifically:
S301: the piece-wise linearization of thermoelectricity fuel cost function, by the fired power generating unit fuel cost argument of function power of the assembling unit
Constant interval [PG,iMin,PG,iMax] n equal part, linearization approximate is carried out to nonlinear fuel cost function in each equal segments
For the linear segmented slope of fuel cost,For power PG,itSegmentation variable, n is specified piecewise interval number, k
For piecewise interval relative to call number, i.e., constant 1 arrive n;
S302: Unit Commitment and Climing constant linearization process
Recursive non-linear feature is presented in the constraint condition of minimum start-stop time, if directly Optimization Solution is needed there are difficulty by this
The nonlinear complementary problem of a complexity is converted into simple linear restriction of equal value, then is solved.
Linearization process: uitIndicate that unit i indicates that unit is in operating status in the operating status of t moment, 1,0 indicates at unit
In shutdown status;yit、zitIndicate beginning on/off state of the unit i in t moment, yit=1 expression unit, which is in, to be started to open
Dynamic state, zit=1, which indicates that unit is in, starts shutdown status.
The minimum of unit i stops machine time-constrain
The minimum runing time of unit i constrains
The minimum downtime of unit i constrains
In formula, RUi、RDiRise for the climbing under unit i normal operating condition and the lower rate of deceleration limits;SUi、SDiIt is opened for unit i
Fall off rate limitation when climbing speed when dynamic and closing;UTi、DTiFor the minimum operation of unit i and downtime;
It has initially been run for unit, downtime, uIt=0For unit initial operating state,P it、It indicates to consider creep speed and start and stop
Unit minimax operational limit after operating status constraint;ζiIt is about unit runing time and when minimum continuous operation
Between between coupled relation, ζi=min { T, (UTi-Ui 0)uIt=0, ξiBe about unit downtime with minimum continuously stop
Coupled relation between the machine time, ξi=min { T, (DTi-Si 0)(1-uIt=0)};
S303: trend constraint linearisation is general using the direct current for ignoring network loss in Unit Commitment and economic load dispatching
Trend constraint has the characteristics that linearisation and rapidity, but only accounts for the relationship of trend and voltage phase angle, and accuracy is not high,
The especially description of voltage and reactive power.
Different from conventional DC flow model, reactive power and voltage constraint are considered here, is obtained using linearization approximate:
Node voltage meets V ≈ 1.0p.u. to be had according to the expansion of reciprocal function: 1/V ≈ 2-V, is substituted on the left of above formula equal sign,
The linearisation of node active balance equation can be obtained are as follows:
B′ijWith BijDifference be to be eliminated in element from susceptance bii;Similarly, node reactive power equilibrium equation linearly turns to
Linearisation Branch Power Flow equation can similarly be obtained
7. a kind of linear tidal current computing method of three-phase polar coordinate system applied to active distribution network according to claim 1,
It is characterized by: the step S4 specifically:
Step S401: initialization model parameter, including optimized variable initial value, objective function and constraint equation coefficient value, constraint
Condition upper limit value and lower limit value;
Step S402: calling the CPLEX solver based on branch and bound method to solve the linearisation peak regulation Optimized model, obtains each
The power output arrangement of unit.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810937487.4A CN109167383B (en) | 2018-08-17 | 2018-08-17 | Power system peak regulation optimization method based on accurate linearized power network model |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810937487.4A CN109167383B (en) | 2018-08-17 | 2018-08-17 | Power system peak regulation optimization method based on accurate linearized power network model |
Publications (2)
Publication Number | Publication Date |
---|---|
CN109167383A true CN109167383A (en) | 2019-01-08 |
CN109167383B CN109167383B (en) | 2022-04-01 |
Family
ID=64896148
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201810937487.4A Active CN109167383B (en) | 2018-08-17 | 2018-08-17 | Power system peak regulation optimization method based on accurate linearized power network model |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109167383B (en) |
Cited By (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110649639A (en) * | 2019-10-30 | 2020-01-03 | 国网辽宁省电力有限公司阜新供电公司 | Regional power grid optimal scheduling method considering operation and loss cost of electric heating system |
CN111130118A (en) * | 2020-01-09 | 2020-05-08 | 清华大学 | Optimal power flow calculation method of power system based on piecewise linearization |
CN111934362A (en) * | 2020-07-29 | 2020-11-13 | 国网甘肃省电力公司电力科学研究院 | Multi-power-supply coordinated optimization peak regulation method for fluctuation characteristics of renewable energy |
CN112529648A (en) * | 2020-10-15 | 2021-03-19 | 中国南方电网有限责任公司 | Wind power considered combined clearing method and system for continuous quotation |
CN112803479A (en) * | 2021-01-06 | 2021-05-14 | 国网辽宁省电力有限公司 | Thermal power generating unit modeling method and model considering deep peak shaving |
CN112803500A (en) * | 2021-03-23 | 2021-05-14 | 国网山西省电力公司 | Method and system for constructing electric energy and deep peak regulation combined clearing model |
CN113344273A (en) * | 2021-06-08 | 2021-09-03 | 中国农业大学 | Building energy consumption based method and system for adjusting and optimizing peak-valley difference of regional distribution network |
CN113408648A (en) * | 2021-07-07 | 2021-09-17 | 华北电力大学 | Unit combination calculation method combined with deep learning |
CN113517691A (en) * | 2021-07-19 | 2021-10-19 | 海南电网有限责任公司 | Multi-type power supply cooperative scheduling method based on peak-valley time-of-use electricity price |
Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP0442756A1 (en) * | 1990-02-16 | 1991-08-21 | Westinghouse Electric Corporation | Electric power utility |
CN104166946A (en) * | 2014-08-15 | 2014-11-26 | 国家电网公司 | Standby and peak shaving auxiliary service cost allocation method facilitating new energy grid-connected consumption |
EP2818649A1 (en) * | 2013-06-27 | 2014-12-31 | Enrichment Technology Company Ltd. | Combination power plant |
CN104361403A (en) * | 2014-10-27 | 2015-02-18 | 国家电网公司 | Optimal grouping configuration method of distributed generations and microgrid |
CN104578183A (en) * | 2015-01-28 | 2015-04-29 | 国家电网公司 | Tie-line power transmission plan optimization method for improving low-carbon level of electric system |
CN104617591A (en) * | 2015-01-28 | 2015-05-13 | 西安交通大学 | Daily operation manner arranging and peak-load regulating method based on multi-scene new energy power generation simulation |
CN106786790A (en) * | 2016-11-19 | 2017-05-31 | 国网浙江省电力公司 | A kind of long-term many power supply coordinated scheduling methods of provincial power network of aqueous bottle coal nuclear power |
CN107104467A (en) * | 2017-06-05 | 2017-08-29 | 武汉大学 | A kind of Unit Combination optimization method for considering nuclear power peak regulation and its security constraint |
-
2018
- 2018-08-17 CN CN201810937487.4A patent/CN109167383B/en active Active
Patent Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP0442756A1 (en) * | 1990-02-16 | 1991-08-21 | Westinghouse Electric Corporation | Electric power utility |
EP2818649A1 (en) * | 2013-06-27 | 2014-12-31 | Enrichment Technology Company Ltd. | Combination power plant |
CN104166946A (en) * | 2014-08-15 | 2014-11-26 | 国家电网公司 | Standby and peak shaving auxiliary service cost allocation method facilitating new energy grid-connected consumption |
CN104361403A (en) * | 2014-10-27 | 2015-02-18 | 国家电网公司 | Optimal grouping configuration method of distributed generations and microgrid |
CN104578183A (en) * | 2015-01-28 | 2015-04-29 | 国家电网公司 | Tie-line power transmission plan optimization method for improving low-carbon level of electric system |
CN104617591A (en) * | 2015-01-28 | 2015-05-13 | 西安交通大学 | Daily operation manner arranging and peak-load regulating method based on multi-scene new energy power generation simulation |
CN106786790A (en) * | 2016-11-19 | 2017-05-31 | 国网浙江省电力公司 | A kind of long-term many power supply coordinated scheduling methods of provincial power network of aqueous bottle coal nuclear power |
CN107104467A (en) * | 2017-06-05 | 2017-08-29 | 武汉大学 | A kind of Unit Combination optimization method for considering nuclear power peak regulation and its security constraint |
Non-Patent Citations (2)
Title |
---|
侯文庭等: "考虑核电可调度性的风-光-核-水-火多源协调短期优化调度", 《电工技术学报》 * |
梁子鹏等: "考虑风电不确定度的风-火-水-气-核-抽水蓄能多源协同旋转备用优化", 《电网技术》 * |
Cited By (15)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110649639A (en) * | 2019-10-30 | 2020-01-03 | 国网辽宁省电力有限公司阜新供电公司 | Regional power grid optimal scheduling method considering operation and loss cost of electric heating system |
CN110649639B (en) * | 2019-10-30 | 2023-04-07 | 国网辽宁省电力有限公司阜新供电公司 | Regional power grid optimal scheduling method considering operation and loss cost of electric heating system |
CN111130118A (en) * | 2020-01-09 | 2020-05-08 | 清华大学 | Optimal power flow calculation method of power system based on piecewise linearization |
WO2022022099A1 (en) * | 2020-07-29 | 2022-02-03 | 国网甘肃省电力公司电力科学研究院 | Multi-power coordination optimization peak regulation method based on fluctuation characteristics of renewable energy |
CN111934362A (en) * | 2020-07-29 | 2020-11-13 | 国网甘肃省电力公司电力科学研究院 | Multi-power-supply coordinated optimization peak regulation method for fluctuation characteristics of renewable energy |
CN112529648A (en) * | 2020-10-15 | 2021-03-19 | 中国南方电网有限责任公司 | Wind power considered combined clearing method and system for continuous quotation |
CN112803479A (en) * | 2021-01-06 | 2021-05-14 | 国网辽宁省电力有限公司 | Thermal power generating unit modeling method and model considering deep peak shaving |
CN112803479B (en) * | 2021-01-06 | 2023-07-21 | 国网辽宁省电力有限公司 | Thermal power generating unit modeling method and model considering depth peak shaving |
CN112803500A (en) * | 2021-03-23 | 2021-05-14 | 国网山西省电力公司 | Method and system for constructing electric energy and deep peak regulation combined clearing model |
CN112803500B (en) * | 2021-03-23 | 2023-04-14 | 国网山西省电力公司 | Method and system for constructing electric energy and deep peak shaving combined clearing model |
CN113344273A (en) * | 2021-06-08 | 2021-09-03 | 中国农业大学 | Building energy consumption based method and system for adjusting and optimizing peak-valley difference of regional distribution network |
CN113344273B (en) * | 2021-06-08 | 2024-04-19 | 中国农业大学 | Regional power distribution network peak-valley difference adjustment optimization method and system based on building energy |
CN113408648A (en) * | 2021-07-07 | 2021-09-17 | 华北电力大学 | Unit combination calculation method combined with deep learning |
CN113517691A (en) * | 2021-07-19 | 2021-10-19 | 海南电网有限责任公司 | Multi-type power supply cooperative scheduling method based on peak-valley time-of-use electricity price |
CN113517691B (en) * | 2021-07-19 | 2023-01-10 | 海南电网有限责任公司 | Multi-type power supply cooperative scheduling method based on peak-valley time-of-use electricity price |
Also Published As
Publication number | Publication date |
---|---|
CN109167383B (en) | 2022-04-01 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109167383A (en) | Electric system peak regulation optimization method based on exact linearization method power network model | |
AU2019101104A4 (en) | An optimal dispatching method of multi-region power and gas coupled integrated energy system using tiered gas tariff | |
CN109149571B (en) | Energy storage optimal configuration method considering characteristics of system gas and thermal power generating unit | |
CN104242337B (en) | The real time coordination control method of photovoltaic microgrid system | |
CN102361330B (en) | Active power optimizing and distributing method for wind generator unit of wind power station | |
CN104504451B (en) | A kind of thermoelectricity system of multimode operation steam power plant Unit Combination optimization method a few days ago | |
CN111525627B (en) | Day-ahead scheduling method for flexible direct-current transmission system with pumped storage and new energy power generation | |
CN109284878A (en) | A kind of multi-source Optimization Scheduling for considering wind-powered electricity generation, nuclear power, water-storage and coordinating | |
CN108808740B (en) | Multi-time scale scheduling method and device for multi-region interconnected power grid and storage medium | |
CN109301818A (en) | A kind of wide-area distribution type energy-storage system coordinated scheduling method considering peak regulation, pressure regulation demand | |
CN108448646B (en) | Source network coordination peak regulation method considering direct current outgoing power regulation characteristic | |
CN108448619B (en) | AC/DC micro-grid robust scheduling method considering uncontrollable generator power tracking | |
CN105846456A (en) | Alternating current and direct current interconnected power grid wind and fire coordination dynamic economy scheduling optimization method | |
CN107451670B (en) | Day-ahead hierarchical coordination scheduling method for active power distribution network | |
CN107104467A (en) | A kind of Unit Combination optimization method for considering nuclear power peak regulation and its security constraint | |
CN111555366B (en) | Multi-time scale-based microgrid three-layer energy optimization management method | |
Tomin et al. | Deep reinforcement learning for energy microgrids management considering flexible energy sources | |
CN112886645B (en) | New energy power system operation simulation method based on ultrahigh hydrogen energy proportion | |
CN110336329A (en) | Receiving end peak load regulation network control method after extra-high voltage direct-current and new energy participation | |
Moshi et al. | Optimal generation scheduling of small diesel generators in a microgrid | |
CN108197766A (en) | A kind of active distribution network Optimal Operation Model for including micro-capacitance sensor group | |
CN109038532B (en) | Distributed transmission and distribution cooperative unit combination method based on direct current power flow equivalence | |
CN109992818A (en) | The Unit Combination model and method for solving of large-scale wind power participation primary frequency modulation | |
CN108336768B (en) | Wind power plant active power optimization control method | |
CN104659812B (en) | A kind of many microgrid control method for coordinating based on PREDICTIVE CONTROL |
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 |