CN103490449A - Optimization method for multi-energy combined power generation system operation simulation - Google Patents

Optimization method for multi-energy combined power generation system operation simulation Download PDF

Info

Publication number
CN103490449A
CN103490449A CN201310469621.XA CN201310469621A CN103490449A CN 103490449 A CN103490449 A CN 103490449A CN 201310469621 A CN201310469621 A CN 201310469621A CN 103490449 A CN103490449 A CN 103490449A
Authority
CN
China
Prior art keywords
adjustable type
power generating
fired power
generating unit
period
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
CN201310469621.XA
Other languages
Chinese (zh)
Other versions
CN103490449B (en
Inventor
赵冬梅
李星宇
牟澎涛
俞勤政
李龙龙
王建锋
郭威
金小明
张东辉
魏国清
胡剑琛
祁永福
吴锋
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
China South Power Grid International Co ltd
North China Electric Power University
Hainan Power Grid Co Ltd
Original Assignee
China South Power Grid International Co ltd
North China Electric Power University
Hainan Power Grid Co Ltd
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by China South Power Grid International Co ltd, North China Electric Power University, Hainan Power Grid Co Ltd filed Critical China South Power Grid International Co ltd
Priority to CN201310469621.XA priority Critical patent/CN103490449B/en
Publication of CN103490449A publication Critical patent/CN103490449A/en
Application granted granted Critical
Publication of CN103490449B publication Critical patent/CN103490449B/en
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • YGENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02EREDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
    • Y02E10/00Energy generation through renewable energy sources
    • Y02E10/50Photovoltaic [PV] energy
    • Y02E10/56Power conversion systems, e.g. maximum power point trackers
    • YGENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02EREDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
    • Y02E10/00Energy generation through renewable energy sources
    • Y02E10/70Wind energy
    • Y02E10/76Power conversion electric or electronic aspects

Landscapes

  • Supply And Distribution Of Alternating Current (AREA)

Abstract

The invention discloses the technical field of operation and scheduling of power systems, and particularly relates to an optimization method for operation simulation of a multi-energy combined power generation system. Selecting days or a week as a scheduling period, and establishing an expected scene of the output of wind power and photovoltaic energy sources according to the predicted deviation value of the output of intermittent power sources such as wind power, photovoltaic and the like; the method comprises the steps of firstly determining the working position of a hydroelectric generating set with flexibly adjustable unit output in a load curve in a multi-energy combined power generation system, then starting from the initial period of a dispatching cycle, arranging the combination mode of a thermal power generating set, reasonably distributing load to each type of set by taking optimized operation cost as a local target, further forming a load regulation adequacy interval aiming at the unit operated in each period, and finally realizing unit combination arrangement and load distribution among the units in the whole dispatching cycle by continuously selecting and updating the load regulation adequacy interval of the unit in each period. The invention ensures the power balance and stable operation of the system.

Description

A kind of optimization method of multi-energy resource combined electric generating system operation simulation
Technical field
The invention belongs to power system operation and dispatching technique field, relate in particular to a kind of optimization method of multi-energy resource combined electric generating system operation simulation.
Background technology
In recent years, along with the tremendous development of domestic power industry and the carrying out of market-oriented reform, China's electric power system fast development, system scale is day by day huge, and running environment is complicated all the more.Meanwhile, the ratio of the new forms of energy such as wind energy turbine set, photovoltaic plant unit in electric power system progressively increases, this makes power supply architecture occur obviously to change, and the electric power system accelerating transition is the combined generating system that contains the polymorphic type power supplys such as thermoelectricity, water power, nuclear power, wind-powered electricity generation, photovoltaic, energy storage.The permeability of the wind-powered electricity generation of some areas can not be ignored, and randomness and intermittence because wind-powered electricity generation, photovoltaic are exerted oneself, after connecting system, can make the more uncertain factor of systems face.Therefore, the operation of traditional power system dispatching is theoretical needs in addition Improvement and perfection, with the change of adaptive system power supply architecture.Under these circumstances, a kind of optimization method of multi-energy resource combined electric generating system operation simulation is proposed, Study system Unit Combination and scheduling problem within the dispatching cycle of a few days or one week, make exerting oneself of all types of power supplys in system can meet the demand of system power load, guarantee system power balance and stability operation.
Summary of the invention
The object of the invention is to, a kind of optimization method of multi-energy resource combined electric generating system operation simulation is provided, for after solving wind-powered electricity generation unit and photovoltaic generation unit connecting system, due to randomness and the intermittence of the two, the unsettled problem of the electricity generation system caused, thus system power balance and stability operation guaranteed.
To achieve these goals, the technical scheme that the present invention proposes is that a kind of optimization method of multi-energy resource combined electric generating system operation simulation is characterized in that described method comprises:
Step 1: the configuration scheduling cycle also will be divided into T period described dispatching cycle;
Step 2: set each period sampling number of scenes N, and determine that each period wind-powered electricity generation unit and photoelectricity unit of each sampling scene blends the power value
Figure BDA0000393281780000021
wherein, m=1,2 ..., N, j=1,2 ..., T;
Step 3: according to formula P l1, j=P l, j-P f, jcalculate the load P that in each period system, the adjustable type generating set be should bear l1, j, and by P l1, jbe designated as the first duty value such as grade;
In described system, the adjustable type generating set comprises adjustable type Hydropower Unit and adjustable type fired power generating unit;
P l, jbe j the total load that the period system be should bear;
P f, jfor the load that in j period system, the generating set except the adjustable type generating set provides;
Step 4: according to formula
Figure BDA0000393281780000022
calculate the load P that in each period system, the adjustable type fired power generating unit be should bear l2, j, and by P l2, jbe designated as the second duty value such as grade; Simultaneously, calculate each period adjustable type Hydropower Unit accommodation limit value that makes progress with downward accommodation limit value
Figure BDA0000393281780000024
μ hi, jbe the operating state of j period Hi adjustable type Hydropower Unit, when j period, Hi adjustable type Hydropower Unit put into operation, μ hi, j=1; When j period, Hi adjustable type Hydropower Unit do not put into operation, μ hi, j=0;
P hi, jbe power output and the P of j period Hi adjustable type Hydropower Unit hi, jhi* g * (H up, Hi, j-H down, Hi, j) * Q hi, j;
η hiit is the generating efficiency coefficient of Hi adjustable type Hydropower Unit;
G is gravity constant;
H up, Hi, jit is water level on the dam of j period Hi adjustable type Hydropower Unit;
H down, Hi, jit is the level of tail water of j period Hi adjustable type Hydropower Unit;
Q hi, jit is the generating average discharge of j period Hi adjustable type Hydropower Unit;
N hquantity for the adjustable type Hydropower Unit;
μ hp, jbe the state that draws water of j period pumped storage machine, when j the period, pumped storage machine put into operation, μ hp, j=1, when j the period, pumped storage machine did not put into operation, μ hp, j=0;
P hp, jbe the power that draws water of j period pumped storage machine;
Step 5: the state that initially puts into operation of determining the adjustable type fired power generating unit; Specifically, calculate respectively the duty value minimum value such as second of the second duty value maximums such as grade of front 5 periods of dispatching cycle and front 5 periods, sequentially drop into the adjustable type fired power generating unit until the total installation of generating capacity of the adjustable type fired power generating unit put into operation is greater than the second maximum such as duty value such as grade of front 5 periods of dispatching cycle according to the first setting, and the minimum technology of the adjustable type fired power generating unit put into operation the summation of exerting oneself is less than the duty value minimum values such as second of front 5 periods of dispatching cycle;
Described adjustable type fired power generating unit comprises the adjustable type fired power generating unit that can not shut down, adjustable type fired power generating unit and the gas turbine unit that can shut down;
Step 6: set initial period t=5;
Step 7: according to formula
Figure BDA0000393281780000034
calculate the load that current period t adjustable type fired power generating unit be should bear under each sampling scene
Figure BDA0000393281780000031
and will
Figure BDA0000393281780000032
be designated as the C grade duty value;
Step 8: the adjustable type fired power generating unit distribution load put into operation under each the sampling scene for current period t
Figure BDA0000393281780000033
determine the Load Regulation abundant intensity interval of the adjustable type fired power generating unit put into operation under each sampling scene of current period t simultaneously;
Step 9: calculate the operating cost of the adjustable type fired power generating unit that each sampling of current period t put into operation under scene, select the adjustable type fired power generating unit that put into operation under sampling scene τ corresponding to the minimum value of operating cost and the load of distribution thereof
Figure BDA0000393281780000041
as the best of current period t adjustable type fired power generating unit scheme that puts into operation; Simultaneously, select the Load Regulation abundant intensity interval of the adjustable type fired power generating unit that put into operation under sampling scene τ corresponding to the minimum value of operating cost to regulate the abundant intensity interval as the optimum load of current period t adjustable type fired power generating unit;
Step 10: make t=t+1, calculate the duty value maximum such as second of continuous 4 periods before current period and current period, and the duty value minimum value such as second of continuous 4 periods before current period and current period;
If the total installation of generating capacity that the adjustable type fired power generating unit of the previous period of current period puts into operation is less than the duty value maximum such as second of continuous 4 periods before current period and current period, according to the first setting, sequentially drop into the adjustable type fired power generating unit also do not put into operation, until the total installation of generating capacity of the adjustable type fired power generating unit put into operation is greater than the duty value maximum such as second of continuous 4 periods before current period and current period;
The summation if the minimum technology of adjustable type fired power generating unit that the duty value maximum such as second of continuous 4 periods is less than the previous period of current period before current period and current period is exerted oneself, the adjustable type fired power generating unit that stops having put into operation according to the inverted order of the first setting order, until the minimum technology of the adjustable type fired power generating unit put into operation is exerted oneself, summation is less than the duty value minimum value such as second of continuous 4 periods before current period and current period;
Step 11: judge that whether the adjustable type fired power generating unit that adjustable type fired power generating unit that current time puts into operation puts into operation with the previous moment of current time is identical, if the adjustable type fired power generating unit that the adjustable type fired power generating unit that current time puts into operation puts into operation with the previous moment of current time is identical, perform step 12; Otherwise, perform step 15;
Step 12: calculate adjustable type fired power generating unit load variations amount under each sampling scene of current period t
Figure BDA0000393281780000051
if Δ P>0, perform step 13; If Δ P<0, perform step 14; If Δ P=0, perform step 15;
Step 13: sequentially adjust according to the second setting the adjustable type generating set that the current period put into operation and adjustable capacity is still arranged, until the adjustable type generating set that participates in adjusting meets &Delta;P = &Sigma; Hi = 1 N H &lambda; Hi P Hi , t , u + &Sigma; Gk = 1 N G &lambda; Gk P Gk , t , u m , Then perform step 15;
Wherein, λ hifor current period Hi platform adjustable type Hydropower Unit participates in the state of adjusting, when current period Hi platform adjustable type Hydropower Unit participates in adjusting, λ hi=1; Current period Hi platform adjustable type Hydropower Unit has neither part nor lot in while adjusting, λ hi=0;
P hi, t, uthe increase variable quantity of exerting oneself for current period Hi platform adjustable type Hydropower Unit;
λ gkin the fired power generating unit put into operation for the current period, Gk platform fired power generating unit participates in the state of adjusting, when in the fired power generating unit that puts into operation of current period, Gk platform fired power generating unit participates in adjusting, and λ gK=1; In the fired power generating unit that puts into operation of current period, Gk platform fired power generating unit has neither part nor lot in while adjusting, λ gk=0;
Figure BDA0000393281780000053
for the increase variable quantity of exerting oneself of Gk platform fired power generating unit in the fired power generating unit put into operation under m of current period sampling background;
N hquantity for the adjustable type Hydropower Unit;
N gquantity for the adjustable type fired power generating unit that put into operation;
Step 14: sequentially adjust according to the second setting the adjustable type generating set that the current period put into operation and adjustable capacity is still arranged, until the adjustable type generating set that participates in adjusting meets
| &Delta;P | = &Sigma; Hi = 1 N H &lambda; Hi P Hi , t , d + &Sigma; Gk = 1 N G &lambda; Gk P Gk , t , d m ;
Wherein, λ hifor current period Hi platform adjustable type Hydropower Unit participates in the state of adjusting, when current period Hi platform adjustable type Hydropower Unit participates in adjusting, λ hi=1; Current period Hi platform adjustable type Hydropower Unit has neither part nor lot in while adjusting, λ hi=0;
P hi, t, dthe minimizing variable quantity of exerting oneself for current period Hi platform adjustable type Hydropower Unit;
λ gkin the fired power generating unit put into operation for the current period, Gk platform fired power generating unit participates in the state of adjusting, when in the fired power generating unit that puts into operation of current period, Gk platform fired power generating unit participates in adjusting, and λ gK=1; In the fired power generating unit that puts into operation of current period, Gk platform fired power generating unit has neither part nor lot in while adjusting, λ gk=0;
Figure BDA0000393281780000061
for the minimizing variable quantity of exerting oneself of Gk platform fired power generating unit in the fired power generating unit put into operation under m of current period sampling background;
N hquantity for the adjustable type Hydropower Unit;
N gquantity for the adjustable type fired power generating unit that put into operation;
Step 15: the adjustable type fired power generating unit distribution load put into operation under each the sampling scene for the current period
Figure BDA0000393281780000062
calculate the Load Regulation abundant intensity interval of the adjustable type fired power generating unit put into operation under each sampling scene of current period simultaneously;
Step 16: calculate the operating cost of the adjustable type fired power generating unit that each sampling of current period t put into operation under scene, select the adjustable type fired power generating unit that put into operation under sampling scene τ corresponding to the minimum value of operating cost and the load of distribution thereof
Figure BDA0000393281780000063
as the best of current period t adjustable type fired power generating unit scheme that puts into operation; Simultaneously, select the Load Regulation abundant intensity interval of the adjustable type fired power generating unit that put into operation under sampling scene τ corresponding to the minimum value of operating cost to regulate the abundant intensity interval as the optimum load of current period t adjustable type fired power generating unit;
Step 17: judge t > whether N sets up, if t > N, perform step 18; Otherwise, return to step 10;
Step 18: finish.
Wind-powered electricity generation unit and the photoelectricity unit of described each sampling scene of definite each period blend the power value comprise following sub-step:
Sub-step 101: respectively according to the historical data of wind-powered electricity generation unit and photoelectricity unit output, make the histogram of wind-powered electricity generation unit output prediction deviation and the histogram of photoelectricity unit output prediction deviation;
Sub-step 102: by approximating method, obtain respectively each period wind-powered electricity generation unit output prediction deviation probability density function and photoelectricity unit output prediction deviation probability density function;
Sub-step 103: each period wind-powered electricity generation unit output prediction deviation probability density function and photoelectricity unit output prediction deviation probability density function are carried out respectively to integral operation, obtain wind-powered electricity generation unit output prediction deviation probability density function and the photoelectricity unit output prediction deviation probability density function of each period;
Sub-step 104: utilize formula calculate the basic sampled value under each sampling scene of each period; Wherein, m=1,2 ..., N, N is each period sampling number of scenes, U is the random number on interval (0,1);
Sub-step 105: utilize formula
Figure BDA0000393281780000073
calculate each sampling scene leeward group of motors of each period prediction deviation of exerting oneself, and utilize formula
Figure BDA0000393281780000074
calculate photoelectricity unit output prediction deviation under each sampling scene of each period; Wherein, for the inverse function of wind-powered electricity generation unit output prediction deviation probability density function,
Figure BDA0000393281780000076
inverse function for photoelectricity unit output prediction deviation probability density function; J=1,2 ..., T, the time hop count that T is dispatching cycle;
Sub-step 106: according to formula wind-powered electricity generation unit and photoelectricity unit blend the power value;
Wherein, P w, jfor each period wind-powered electricity generation unit output predicted value; P p, jfor each period photoelectricity unit output predicted value.
Each period adjustable type Hydropower Unit of described calculating accommodation limit value that makes progress
Figure BDA0000393281780000081
adopt formula
P H , j op = &Sigma; Hi = 1 N H &mu; Hi , j P Hi , j op ;
Wherein,
Figure BDA0000393281780000083
be the upwards pondage of j period Hi adjustable type Hydropower Unit, and
P Hi , j op = P Hi , j max - P Hi , j ;
Figure BDA0000393281780000085
it is the peak power output of j period Hi adjustable type Hydropower Unit;
P hi, jit is the power output of j period Hi adjustable type Hydropower Unit;
μ hi, jbe the operating state of j period Hi adjustable type Hydropower Unit, when j period, Hi adjustable type Hydropower Unit put into operation, μ hi, j=1; When j period, Hi adjustable type Hydropower Unit do not put into operation, μ hi, j=0.
The downward accommodation limit value of each period adjustable type Hydropower Unit of described calculating
Figure BDA0000393281780000086
adopt formula
P H , j ne = &Sigma; Hi = 1 N H &mu; Hi , j P Hi , j ne ;
Wherein,
Figure BDA0000393281780000088
be the downward pondage of j period Hi adjustable type Hydropower Unit, and
P Hi , j ne = P Hi , j - P Hi , j min ;
Figure BDA00003932817800000810
the minimum that is j period Hi adjustable type Hydropower Unit is forced power output;
P hi, jit is the power output of j period Hi adjustable type Hydropower Unit;
μ hi, jbe the operating state of j period Hi adjustable type Hydropower Unit, when j period, Hi adjustable type Hydropower Unit put into operation, μ hi, j=1; When j period, Hi adjustable type Hydropower Unit do not put into operation, μ hi, j=0.
Describedly according to the first setting, sequentially drop into the adjustable type fired power generating unit specifically; first by minimum operating ratio consumption, order from low to high drops into the adjustable type fired power generating unit that can not shut down; by minimum specific consumption, order from low to high drops into the adjustable type fired power generating unit that can shut down again, and finally by the power output regulations speed, order from high to low drops into gas turbine unit.
The adjustable type fired power generating unit distribution load put into operation under described each sampling scene for current period t
Figure BDA0000393281780000094
specifically comprise following sub-step:
Sub-step 201: set the adjustable type fired power generating unit put into operation the out-of-limit sign nFlag=1 that exerts oneself;
Sub-step 202: according to formula
Figure BDA0000393281780000091
calculate the micro-gaining rate of consumption of fired power generating unit; Wherein,
Figure BDA0000393281780000092
for C grade duty value and m=1,2 ..., N;
N is each period sampling number of scenes;
N ' gfor current period t put into operation and the exerted oneself quantity of the adjustable type fired power generating unit that do not transfinite;
A gk 'for current period t the Gk ' platform put into operation and the exerted oneself quadratic term coefficient of consumption characteristic quadratic function of the adjustable type fired power generating unit that do not transfinite;
B gk 'for current period t the Gk ' platform put into operation and the exerted oneself Monomial coefficient of consumption characteristic quadratic function of the adjustable type fired power generating unit that do not transfinite;
Sub-step 203: make Gk=1;
Sub-step 204: according to formula
Figure BDA0000393281780000093
calculate adjustable type fired power generating unit value of exerting oneself of having put into operation of Gk platform under each sampling scene of current period t;
A gkthe quadratic term coefficient of the consumption characteristic quadratic function of the adjustable type fired power generating unit put into operation for current period t Gk platform;
B gkthe Monomial coefficient of the consumption characteristic quadratic function of the adjustable type fired power generating unit put into operation for current period t Gk platform;
Sub-step 205: if
Figure BDA0000393281780000101
transfinite, order
Figure BDA0000393281780000102
maximum output value during for this adjustable type fired power generating unit put into operation not out-of-limit, and make nFlag=0;
Sub-step 206: judgement Gk>N gwhether set up, if Gk N g, carry out sub-step 207; Otherwise, make Gk=Gk+1, return to sub-step 204;
N gthe quantity of the adjustable type fired power generating unit put into operation for current period t;
Sub-step 207: judge whether nFlag=0 sets up, if nFlag=0 returns to sub-step 201; Otherwise, carry out sub-step 208;
Sub-step 208: make the load that under each sampling scene of current period t, the adjustable type fired power generating unit that put into operation of Gk platform provides be
Figure BDA0000393281780000103
The operating cost of the adjustable type fired power generating unit put into operation under each sampling scene of the current period t of described calculating adopts formula F cos t ( P G , t m ) = &Sigma; Gk = 1 N T [ &mu; Gk , t m &times; f k ( P Gk , t m ) + &mu; Gk , t m ( 1 - &mu; Gk , t - 1 m ) S Gk ] ;
Wherein,
Figure BDA0000393281780000105
the operating cost of the adjustable type fired power generating unit put into operation under m the sampling scene for current period t, m=1,2 ..., N, N is the sampling number of scenes in each period;
N tfor adjustable type fired power generating unit quantity;
Figure BDA0000393281780000106
for the current period t Gk platform adjustable type fired power generating unit state that puts into operation, current period t Gk platform adjustable type fired power generating unit puts into operation,
Figure BDA0000393281780000107
current period t Gk platform adjustable type fired power generating unit does not put into operation,
Figure BDA0000393281780000108
Figure BDA0000393281780000109
for the consumption characteristic quadratic function of current period t Gk platform adjustable type fired power generating unit, and
f k ( P Gk , t m ) = A Gk &times; ( P Gk , t m ) 2 + B Gk &times; P Gk , t m + C Gk ;
A gkquadratic term coefficient for the consumption characteristic quadratic function of current period t Gk platform adjustable type fired power generating unit;
B gkmonomial coefficient for the consumption characteristic quadratic function of current period t Gk platform adjustable type fired power generating unit;
C gkconstant term coefficient for the consumption characteristic quadratic function of current period t Gk platform adjustable type fired power generating unit;
S gkfor the payment for initiation of current period t Gk platform adjustable type fired power generating unit is used.
The Load Regulation abundant intensity interval of the adjustable type fired power generating unit put into operation under each sampling scene of described current period t is [ P L 3 , t m - P H , t ne - P G , t ne , m , P L 3 , t m + P H , t op + P G , t op , m ] ;
Wherein,
Figure BDA0000393281780000113
for the downward accommodation limit value of current period t adjustable type Hydropower Unit;
Figure BDA0000393281780000114
for adjustable type thermal power unit operation lower bound under each sampling scene of current period t, and
P G , t ne , m = &Sigma; Gk = 1 N T &mu; Gk , t m P Gk , t ne , m ;
Figure BDA0000393281780000116
for the operation lower bound of current period t Gk platform adjustable type fired power generating unit, and
P Gk , t ne , m = min ( P Gk , t m - P Gk , min , - R Gk , d &times; &Delta;t ) ;
P gk, minit is the Gk platform adjustable type fired power generating unit lower limit of exerting oneself;
R gk, dit is the downward creep speed of Gk platform adjustable type fired power generating unit;
for the current period t Gk platform adjustable type fired power generating unit state that puts into operation, current period t Gk platform adjustable type fired power generating unit puts into operation,
Figure BDA0000393281780000119
current period t Gk platform adjustable type fired power generating unit does not put into operation,
Figure BDA0000393281780000121
The time span that Δ t is current period t;
N tquantity for the adjustable type fired power generating unit;
for the current period t adjustable type Hydropower Unit accommodation limit value that makes progress;
Figure BDA0000393281780000123
for the adjustable type thermal power unit operation upper bound under each sampling scene of current period t, and
P G , t op , m = &Sigma; Gk = 1 N T &mu; Gk , t m P Gk , t op , m ;
for the current period t Gk platform adjustable type thermal power unit operation upper bound, and
P Gk , t op , m = min ( P Gk , max - P Gk , t m , - R Gk , u &times; &Delta;t ) ;
P gk, maxit is the upper limit of exerting oneself of Gk platform adjustable type fired power generating unit;
R gk, uit is the Gk platform adjustable type fired power generating unit creep speed that makes progress;
M=1,2 ..., N, N is the sampling number of scenes in each period.
Describedly sequentially adjust the current period according to the second setting and put into operation and still had the adjustable type generating set of adjustable capacity to be specially, first adjust the adjustable type Hydropower Unit that the current period put into operation and adjustable capacity is still arranged, when all adjustable type Hydropower Unit of having put into operation of current period do not have adjustable capacity, the gas turbine unit that order adjustment has from high to low put into operation by the power output regulations speed, when all gas turbine units of having put into operation of current period do not have adjustable capacity, the adjustable type fired power generating unit of shutting down that order adjustment from low to high puts into operation by minimum specific consumption, when all adjustable type fired power generating unit of shutting down of having put into operation of current period do not have adjustable capacity, the adjustable type fired power generating unit that order adjustment from low to high can not be shut down by minimum operating ratio consumption.
Method provided by the invention, by research system Unit Combination and scheduling problem within the dispatching cycle of a few days or one week, make exerting oneself of all types of power supplys in system can meet the demand of system power load, guarantees system power balance and stability operation.
The accompanying drawing explanation
Fig. 1 is the optimization method flow chart of a kind of multi-energy resource combined electric generating system operation simulation provided by the invention;
Fig. 2 is each sampling scene apparatus for lower wind generating set of each period schematic diagram of exerting oneself;
Fig. 3 is adjustable type Hydropower Unit service position schematic diagram;
Fig. 4 is the adjustable type fired power generating unit distribution load flow chart put into operation under each sampling scene of each period.
Embodiment
Below in conjunction with accompanying drawing, preferred embodiment is elaborated.Should be emphasized that, following explanation is only exemplary, rather than in order to limit the scope of the invention and to apply.
Fig. 1 is the optimization method flow chart of a kind of multi-energy resource combined electric generating system operation simulation provided by the invention.As shown in Figure 1, the optimization method of multi-energy resource combined electric generating system operation simulation provided by the invention comprises
Following steps:
Step 1: the configuration scheduling cycle also will be divided into T period described dispatching cycle.
Step 2: set each period sampling number of scenes N, and determine that each period wind-powered electricity generation unit and photoelectricity unit of each sampling scene blends the power value
Figure BDA0000393281780000131
wherein, m=1,2 ..., N, j=1,2 ..., T.
Wind power generation and photovoltaic generation are all generations of electricity by new energy, usually according to wind turbine generator and photovoltaic generation unit dope power value and prediction deviation value thereof (be the actual value of exerting oneself with prediction exert oneself poor) statistics, generate some scenes of exerting oneself within each period, with the intermittence of exerting oneself and the fluctuation of the generation of electricity by new energy units such as reflection wind-force, photovoltaic.Therefore, determine that each period wind-powered electricity generation unit and photoelectricity unit of each sampling scene blends the power value
Figure BDA0000393281780000132
comprise following sub-step:
Sub-step 101: respectively according to the historical data of wind-powered electricity generation unit and photoelectricity unit output, make the histogram of wind-powered electricity generation unit output prediction deviation and the histogram of photoelectricity unit output prediction deviation.
Sub-step 102: by approximating method, obtain respectively each period wind-powered electricity generation unit output prediction deviation probability density function and photoelectricity unit output prediction deviation probability density function.
Sub-step 103: each period wind-powered electricity generation unit output prediction deviation probability density function and photoelectricity unit output prediction deviation probability density function are carried out respectively to integral operation, obtain wind-powered electricity generation unit output prediction deviation probability density function and the photoelectricity unit output prediction deviation probability density function of each period.
Above-mentioned sub-step 101-103 is those skilled in the art's generating set commonly used prediction deviation computational methods of exerting oneself.Because in each period, each sampling can only have a predicted value under scene, therefore above-mentioned two probability density functions are all invertible functions one to one.
Sub-step 104: calculate the basic sampled value under each sampling scene of each period, its computing formula is as follows:
U m = U N + m - 1 N - - - ( 1 )
In formula (1), m=1,2 ..., N, N is each period sampling number of scenes, U is the random number on interval (0,1).
Sub-step 105: calculate each sampling scene leeward group of motors of each period exert oneself prediction deviation and photoelectricity unit output prediction deviation.Wherein, the exert oneself computing formula of prediction deviation of each sampling scene leeward group of motors of each period is:
X W , j m = F W , j - 1 ( U m ) - - - ( 2 )
Under each sampling scene of each period, the computing formula of photoelectricity unit output prediction deviation is:
X P , j m = F P , j - 1 ( U m ) - - - ( 3 )
In formula (2), inverse function for wind-powered electricity generation unit output prediction deviation probability density function.In formula (3),
Figure BDA0000393281780000145
inverse function for photoelectricity unit output prediction deviation probability density function.J=1,2 ..., T, the time hop count that T is dispatching cycle.
Sub-step 106: the wind-powered electricity generation unit and the photoelectricity unit that calculate each sampling scene of each period blend the power value, and the computing formula of this mixing value of exerting oneself is:
P New , j m = P W , j + X W , j m + P P , j + X P , j m - - - ( 4 )
In formula (4), P w, jfor each period wind-powered electricity generation unit output predicted value, P p, jfor each period photoelectricity unit output predicted value.Fig. 2 is each sampling scene apparatus for lower wind generating set of each period schematic diagram of exerting oneself.Each period, each sampling scene apparatus for lower wind generating set value of exerting oneself should be under each sampling scene of each period, photoelectricity unit output value should be
Figure BDA0000393281780000153
Step 3: calculate the load P that in each period system, the adjustable type generating set be should bear l1, j, and by P l1, jbe designated as the first duty value such as grade.
In the multi-energy resource combined electric generating system, generating set is divided into unadjustable generator group and adjustable type generating set usually.Exerting oneself in whole power generation process of unadjustable generator group is all that can't change stable exerted oneself.The exerting oneself of unadjustable generator group comprise without forcing of the exerting oneself of adjustment type Hydropower Unit, adjustable type Hydropower Unit exert oneself, forcing of thermoelectricity unit exert oneself and the exerting oneself of nuclear power generating sets, these are exerted oneself owing to being steadily exerting oneself of system operation, service position is at the load curve base portion, and it is exerted oneself and does not participate in the adjustment process of power-balance.Adjustable type generating set in system comprises adjustable type Hydropower Unit and adjustable type fired power generating unit.
Therefore, the load P that in each period system, the adjustable type generating set be should bear l1, jcomputing formula be:
P L1,j=P L,j-P f,j (5)
In formula (5), P l, jbe j the total load that the period system be should bear, P f, jfor the load that in j period system, the generating set except the adjustable type generating set provides.
Step 4: calculate the load P that in each period system, the adjustable type fired power generating unit be should bear l2, j, and by P l2, jbe designated as the second duty value such as grade.
Due to P l1, jit is the load that in each period system, the adjustable type generating set be should bear, therefore the load that in each period system, the adjustable type fired power generating unit be should bear, the load that namely in each period system, the adjustable type generating set be should bear deducts the load that in each period system, the adjustable type Hydropower Unit be should bear.
Adjustable type Hydropower Unit and pumped storage machine are in load peak time generating.Pumped storage machine is arranged in the low ebb load period of every day in dispatching cycle and is drawn water, in the next adjustable type Hydropower Unit that is considered as of generating operation mode.According to the hydrologic condition of adjustable type Hydropower Unit and the plan energy output of every day, determine the service position of adjustable type Hydropower Unit on load curve, as shown in Figure 3.In Fig. 3, P h, upfor the service position upper limit of corresponding period adjustable type Hydropower Unit on load curve, P h, downfor the service position lower limit of corresponding period adjustable type Hydropower Unit on load curve, in real system, expection value of exerting oneself of adjustable type Hydropower Unit is: P H , j = P H , up - P H , down , P L 1 , j &GreaterEqual; P H , up P L 1 , j - P H , down , P H , down < P L 1 , j < P H , up 0 , P L 1 , j &le; P H , down , It is also expection value of exerting oneself of the adjustable type Hydropower Unit of j period.
In dispatching cycle, the expection energy output of adjustable type Hydropower Unit, be according to power station hydrologic condition, the definite definite value of storage capacity control plan.Therefore, can utilize the mode moved down from the daily load peak value fixed step size of load curve, determine the position of the work bound of adjustable type Hydropower Unit, the difference of assurance bound power approaches the anticipation of adjustable type Hydropower Unit as far as possible and exerts oneself.
So, for load P that in each period system, the adjustable type fired power generating unit be should bear l2, j, its computing formula should be:
P L 2 , j = P L 1 , j - &Sigma; Hi = 1 N H &mu; Hi , j P Hi , j + &mu; Hp , j P Hp , j - - - ( 6 )
In formula (6), μ hi, jbe the operating state of j period Hi adjustable type Hydropower Unit, when j period, Hi adjustable type Hydropower Unit put into operation, μ hi, j=1; When j period, Hi adjustable type Hydropower Unit do not put into operation, μ hi, j=0.N hquantity for the adjustable type Hydropower Unit.μ hp, jbe the state that draws water of j period pumped storage machine, when j the period, pumped storage machine put into operation, μ hp, j=1, when j the period, pumped storage machine did not put into operation, μ hp, j=0.P hp, jbe the power that draws water of j period pumped storage machine.P hi, jbe the power output of j period Hi adjustable type Hydropower Unit, its computing formula is:
P Hi,jHi×g×(H up,Hi,j-H down,Hi,j)×Q Hi,j (7)
In formula (7), η hibe the generating efficiency coefficient of Hi adjustable type Hydropower Unit, g is gravity constant and gets g=9.81, unit: newton/kilogram, H up, Hi, jbe water level on the dam of j period Hi adjustable type Hydropower Unit, H down, Hi, jbe the level of tail water of j period Hi adjustable type Hydropower Unit, Q hi, jit is the generating average discharge of j period Hi adjustable type Hydropower Unit.
Meanwhile, calculate respectively each period adjustable type Hydropower Unit accommodation limit value that makes progress
Figure BDA0000393281780000171
with downward accommodation limit value
Figure BDA0000393281780000172
above-mentioned two values can be used in the calculating of subsequent step.
Calculate each period adjustable type Hydropower Unit accommodation limit value that makes progress
Figure BDA0000393281780000173
the formula adopted is:
P H , j op = &Sigma; Hi = 1 N H &mu; Hi , j P Hi , j op - - - ( 8 )
In formula (8),
Figure BDA0000393281780000175
be j period Hi adjustable type Hydropower Unit upwards pondage and
Figure BDA0000393281780000176
Figure BDA0000393281780000177
be the peak power output of j period Hi adjustable type Hydropower Unit, P hi, jit is the power output of j period Hi adjustable type Hydropower Unit.μ hi, jbe the operating state of j period Hi adjustable type Hydropower Unit, when j period, Hi adjustable type Hydropower Unit put into operation, μ hi, j=1; When j period, Hi adjustable type Hydropower Unit do not put into operation, μ hi, j=0.
Calculate the downward accommodation limit value of each period adjustable type Hydropower Unit
Figure BDA0000393281780000181
the formula adopted is:
P H , j ne = &Sigma; Hi = 1 N H &mu; Hi , j P Hi , j ne - - - ( 9 )
In formula (9),
Figure BDA0000393281780000183
be j period Hi adjustable type Hydropower Unit downward pondage and the minimum that is j period Hi adjustable type Hydropower Unit is forced power output, P hi, jit is the power output of j period Hi adjustable type Hydropower Unit.μ hi, jbe the operating state of j period Hi adjustable type Hydropower Unit, when j period, Hi adjustable type Hydropower Unit put into operation, μ hi, j=1; When j period, Hi adjustable type Hydropower Unit do not put into operation, μ hi, j=0.
Step 5: the state that initially puts into operation of determining the adjustable type fired power generating unit.
In hybrid power system, the adjustable type fired power generating unit generally comprises the adjustable type fired power generating unit that can not shut down, adjustable type fired power generating unit and the gas turbine unit that can shut down.The adjustable type fired power generating unit that can not shut down is not carried out start and stop peak regulation (not considering the trouble hunting special case) within dispatching cycle.The adjustable type fired power generating unit that can shut down can be subject to dispatching command to arrange start and stop.The gas turbine unit capacity is little, fast response time, does and regulates the power supply use.
Consider the small-sized thermoelectricity and the combustion gas unit that participate in the start and stop peak regulation, the restriction of minimum start-stop time is arranged, while therefore considering the Unit Combination scheme of current period, need associating prediction load condition afterwards, in order to allow unit have certain buffer time started or close down, this method is selected continuous 5 periods backward that comprise the current period.So, when determining the state that initially puts into operation of adjustable type fired power generating unit, first calculate respectively the duty value minimum value such as second of the second duty value maximums such as grade of front 5 periods of dispatching cycle and front 5 periods, then according to the first setting, sequentially drop into the adjustable type fired power generating unit.
Herein; sequentially dropping into the adjustable type fired power generating unit according to the first setting is specially; first by minimum operating ratio consumption, order from low to high drops into the adjustable type fired power generating unit that can not shut down; by minimum specific consumption, order from low to high drops into the adjustable type fired power generating unit that can shut down again, and finally by the power output regulations speed, order from high to low drops into gas turbine unit.Drop into successively the adjustable type fired power generating unit, when the total installation of generating capacity of the adjustable type fired power generating unit put into operation is greater than the second maximums such as duty value such as grade of front 5 periods of dispatching cycle, and the minimum technology of the adjustable type fired power generating unit put into operation is exerted oneself summation while being less than the duty value minimum value such as second of front 5 periods of dispatching cycle, stop dropping into the adjustable type fired power generating unit, the initially state of putting into operation of adjustable type fired power generating unit is determined complete.
Above-mentioned, the minimum technology of adjustable type fired power generating unit is exerted oneself and is referred to, guarantees that boiler maintains in the smooth combustion situation of (be unlikely to cease fire or depart from temperature and the pressure of stable state), the minimal power values that unit can send.The minimum technology of every fired power generating unit provides after exerting oneself and all having been manufactured by unit producer.
Step 6: set initial period t=5.
Step 7: calculate the load that current period t adjustable type fired power generating unit be should bear under each sampling scene
Figure BDA0000393281780000191
and will
Figure BDA0000393281780000192
be designated as the C grade duty value.The C grade duty value
Figure BDA0000393281780000193
computing formula be:
P L 3 , t m = P L 2 , t - P New , t m - - - ( 10 )
Step 8: the adjustable type fired power generating unit distribution load put into operation under each the sampling scene for current period t
Figure BDA0000393281780000195
determine the Load Regulation abundant intensity interval of the adjustable type fired power generating unit put into operation under each sampling scene of current period t simultaneously.
The C grade duty value
Figure BDA0000393281780000196
be the total load that current period t adjustable type fired power generating unit be should bear under each sampling scene, specific to every adjustable type fired power generating unit, how many loads of should bear under each sampling scene of current period t, should calculate by following method.
At first, every fired power generating unit has consumption characteristic quadratic function separately, and this consumption characteristic quadratic function manufacturer of fired power generating unit sometimes can directly provide.By also obtaining fired power generating unit consumption characteristic curve fitting.The consumption characteristic quadratic function that we set Gk platform fired power generating unit is
Figure BDA0000393281780000197
wherein, A kbe the quadratic term coefficient of the consumption characteristic quadratic function of k platform fired power generating unit, B kbe the Monomial coefficient of the consumption characteristic quadratic function of k platform fired power generating unit, C kbe the constant term coefficient of the consumption characteristic quadratic function of k platform fired power generating unit, P kit is the load of k platform fired power generating unit.
When fired power generating unit, when exerting oneself variation, if the micro-gaining rate equalization of consumption, the regulating effect optimum that load-responsive changes, describe by the math equation group, is
f op = 2 A 1 P 1 + B 1 = 2 A 2 P 2 + B 2 = &CenterDot; &CenterDot; &CenterDot; 2 A M P M + B M &Sigma; k = 1 M P k = P L - - - ( 11 )
In formula (11), the quantity that M is fired power generating unit, P lfor the total load of all fired power generating unit, f opfor the micro-gaining rate of consumption.According to formula (11), can solve:
f op = P L + &Sigma; k = 1 M ( 0.5 B k / A k ) &Sigma; k = 1 M ( 2 A k ) - 1 - - - ( 12 )
Every load that the adjustable type fired power generating unit be should bear under each sampling scene of current period, should be not out-of-limit value.Therefore, the present invention adopts the mode of loop iteration, by the micro-gaining rate of adjustable type fired power generating unit consumption put into operation, regulates the value of exerting oneself of the adjustable type fired power generating unit put into operation, thereby guarantees that exerting oneself of every adjustable type fired power generating unit is all not out-of-limit.Therefore, according to formula (12), the concrete computational process of the load that every adjustable type fired power generating unit be should bear under each sampling scene of current period t as shown in Figure 4, comprising:
Sub-step 201: set the adjustable type fired power generating unit put into operation the out-of-limit sign nFlag=1 that exerts oneself.
Sub-step 202: calculate the micro-gaining rate of adjustable type fired power generating unit consumption put into operation.During beginning, may have the value of exerting oneself of the adjustable type fired power generating unit put into operation in out-of-limit state.While calculating the micro-gaining rate of consumption, should only consider the adjustable type fired power generating unit put into operation that those values of exerting oneself are not out-of-limit.
Adopt formula (13) to calculate the micro-gaining rate of adjustable type fired power generating unit consumption put into operation:
f op = P L 3 , t m + &Sigma; Gk &prime; = 1 N G &prime; ( 0.5 B Gk &prime; / A G k &prime; ) &Sigma; Gk &prime; = 1 N G &prime; ( 2 A Gk &prime; ) - 1 - - - ( 13 )
In formula (13),
Figure BDA0000393281780000212
for C grade duty value and m=1,2 ..., N, N is each period sampling number of scenes, N ' gfor current period t put into operation and the exerted oneself quantity of the adjustable type fired power generating unit that do not transfinite, A gk 'for current period t the Gk ' platform put into operation and the exerted oneself quadratic term coefficient of consumption characteristic quadratic function of the adjustable type fired power generating unit that do not transfinite, B gk 'for current period t the Gk ' platform put into operation and the exerted oneself Monomial coefficient of consumption characteristic quadratic function of the adjustable type fired power generating unit that do not transfinite.
Sub-step 203: make Gk=1.
Sub-step 204: utilize the micro-gaining rate f of consumption op, upgrading adjustable type fired power generating unit value of exerting oneself that each has put into operation, its formula is:
P Gk , t m = f op - B Gk 2 A Gk - - - ( 14 )
In formula (14), A gkthe quadratic term coefficient of the consumption characteristic quadratic function of the adjustable type fired power generating unit put into operation for current period t Gk platform, B gkthe Monomial coefficient of the consumption characteristic quadratic function of the adjustable type fired power generating unit put into operation for current period t Gk platform.
Sub-step 205: if
Figure BDA0000393281780000214
transfinite, order
Figure BDA0000393281780000215
maximum output value during for this adjustable type fired power generating unit put into operation not out-of-limit, and make nFlag=0.
Sub-step 206: judgement Gk>N gwhether set up, if Gk N g, carry out sub-step 207; Otherwise, make Gk=Gk+1, return to sub-step 204.N gthe quantity of the adjustable type fired power generating unit put into operation for current period t.
After the value of exerting oneself of all adjustable type fired power generating unit that put into operation has been upgraded, may also have situation about transfiniting and exist, now out-of-limit sign nFlag is set to 0, carry out further iteration by sub-step 207.
Sub-step 207: judge whether nFlag=0 sets up, if nFlag=0 returns to sub-step 201; Otherwise, carry out sub-step 208.
Sub-step 208: make the load that under each sampling scene of current period t, the adjustable type fired power generating unit that put into operation of Gk platform provides be
Figure BDA0000393281780000221
If in certain renewal process, out-of-limit sign nFlag is 1, the value of exerting oneself that all adjustable type fired power generating unit that put into operation are described is not out-of-limit, and the load that the adjustable type fired power generating unit that now under each sampling scene of t of current period, the Gk platform has put into operation provides is for being its value of exerting oneself
Figure BDA0000393281780000222
Step 9: calculate the operating cost of the adjustable type fired power generating unit put into operation under each sampling scene of current period t, its computing formula is:
F cos t ( P G , t m ) = &Sigma; Gk = 1 N T [ &mu; Gk , t m &times; f k ( P Gk , t m ) + &mu; Gk , t m ( 1 - &mu; Gk , t - 1 m ) S Gk ] - - - ( 15 )
In formula (15), the operating cost of the adjustable type fired power generating unit put into operation under m the sampling scene for current period t, m=1,2 ..., N, N is the sampling number of scenes in each period, N tfor adjustable type fired power generating unit quantity.
Figure BDA0000393281780000225
for the current period t Gk platform adjustable type fired power generating unit state that puts into operation, current period t Gk platform adjustable type fired power generating unit puts into operation,
Figure BDA0000393281780000226
current period t Gk platform adjustable type fired power generating unit does not put into operation,
Figure BDA0000393281780000227
for the consumption characteristic quadratic function of current period t Gk platform adjustable type fired power generating unit, and
Figure BDA0000393281780000228
Figure BDA0000393281780000229
for the quadratic term coefficient of the consumption characteristic quadratic function of current period t Gk platform adjustable type fired power generating unit, B gkfor the Monomial coefficient of the consumption characteristic quadratic function of current period t Gk platform adjustable type fired power generating unit, C gkfor the constant term coefficient of the consumption characteristic quadratic function of current period t Gk platform adjustable type fired power generating unit, S gkfor the payment for initiation of current period t Gk platform adjustable type fired power generating unit is used.
Simultaneously, determine the Load Regulation abundant intensity interval of the adjustable type fired power generating unit put into operation under each sampling scene of current period t.This interval is actual is
Figure BDA0000393281780000231
wherein,
Figure BDA0000393281780000232
for the downward accommodation limit value of current period t adjustable type Hydropower Unit,
Figure BDA0000393281780000233
for current period t each the sampling scene under adjustable type thermal power unit operation lower bound and
Figure BDA0000393281780000234
Figure BDA0000393281780000235
for the operation lower bound of current period t Gk platform adjustable type fired power generating unit and
Figure BDA0000393281780000236
p gk, minbe the Gk platform adjustable type fired power generating unit lower limit of exerting oneself, RGk, d is the downward creep speed of Gk platform adjustable type fired power generating unit.
Figure BDA0000393281780000237
for the current period t adjustable type Hydropower Unit accommodation limit value that makes progress,
Figure BDA0000393281780000238
for current period t each the sampling scene under the adjustable type thermal power unit operation upper bound and
Figure BDA0000393281780000239
Figure BDA00003932817800002310
for the current period t Gk platform adjustable type thermal power unit operation upper bound and p gk, maxbe the upper limit of exerting oneself of Gk platform adjustable type fired power generating unit, R gk, uit is the Gk platform adjustable type fired power generating unit creep speed that makes progress. for the current period t Gk platform adjustable type fired power generating unit state that puts into operation, current period t Gk platform adjustable type fired power generating unit puts into operation,
Figure BDA00003932817800002313
current period t Gk platform adjustable type fired power generating unit does not put into operation,
Figure BDA00003932817800002314
the time span that Δ t is current period t, N tfor the quantity of adjustable type fired power generating unit, m=1,2 ..., N, N is the sampling number of scenes in each period.
Under each sampling scene of current period t, all can calculate one value, select all
Figure BDA00003932817800002316
value (m=1,2 ..., N) middle minimum value, when
Figure BDA00003932817800002317
sampling scene value while getting minimum value is made as τ, the adjustable type fired power generating unit put into operation under the scene of sampling τ and the load of distribution thereof
Figure BDA00003932817800002318
the best of current exactly period t adjustable type fired power generating unit scheme that puts into operation.The Load Regulation abundant intensity interval of the adjustable type fired power generating unit simultaneously, put into operation under sampling scene τ
Figure BDA0000393281780000241
the optimum load that is current period t adjustable type fired power generating unit is regulated the abundant intensity interval.
Step 10: make t=t+1, calculate the duty value maximum such as second of continuous 4 periods before current period and current period, and the duty value minimum value such as second of continuous 4 periods before current period and current period.
If the total installation of generating capacity that the adjustable type fired power generating unit of the previous period of current period puts into operation is less than the duty value maximum such as second of continuous 4 periods before current period and current period, according to the first setting, sequentially drop into the adjustable type fired power generating unit also do not put into operation, until the total installation of generating capacity of the adjustable type fired power generating unit put into operation is greater than the duty value maximum such as second of continuous 4 periods before current period and current period.
The summation if the minimum technology of adjustable type fired power generating unit that the duty value maximum such as second of continuous 4 periods is less than the previous period of current period before current period and current period is exerted oneself, the adjustable type fired power generating unit that stops having put into operation according to the inverted order of the first setting order, until the minimum technology of the adjustable type fired power generating unit put into operation is exerted oneself, summation is less than the duty value minimum value such as second of continuous 4 periods before current period and current period.
The adjustable type fired power generating unit that stops having put into operation according to the inverted order of the first setting order is actual is exactly; first by power output regulations speed order stopping gas turbine unit from low to high; by minimum specific consumption, order from high to low stops the adjustable type fired power generating unit that can shutting down again, finally by minimum operating ratio consumption order from high to low, stops the adjustable type fired power generating unit that can not shutting down.
Step 11: judge that whether the adjustable type fired power generating unit that adjustable type fired power generating unit that current time puts into operation puts into operation with the previous moment of current time is identical.
If the adjustable type fired power generating unit that the adjustable type fired power generating unit that current time puts into operation puts into operation with the previous moment of current time is identical, the generating set that explanation does not newly increase, the generating set of also newly not closing, now need the adjustable type generating set to put into operation is finely tuned, jump directly to step 12.
If the adjustable type fired power generating unit that the adjustable type fired power generating unit that current time puts into operation puts into operation with the previous moment of current time is not identical, explanation has the unit newly increased or the unit of newly closing, direct step 15, select put into operation scheme and optimum load of current time the best to regulate the abundant intensity interval.
Step 12: calculate adjustable type fired power generating unit load variations amount under each sampling scene of current period t
&Delta;P = P L 3 , t m - P L 3 , t - 1 &tau; .
If Δ P > 0, the explanation load has increased, and need to be finely tuned the adjustable type generating set put into operation, now skips to step 13, makes the adjustable type generating set put into operation bear the load variations amount of increase.
If Δ P<0, the explanation load has reduced, and need to be finely tuned the adjustable type generating set put into operation, now skips to step 14, makes the adjustable type generating set put into operation bear the load variations amount of minimizing.
If Δ P=0, finely tuned without the adjustable type generating set to putting into operation, jump directly to step 15, select put into operation scheme and optimum load of current time the best to regulate the abundant intensity interval.
Step 13: sequentially adjust according to the second setting the adjustable type generating set that the current period put into operation and adjustable capacity is still arranged, until the adjustable type generating set that participates in adjusting meets &Delta;P = &Sigma; Hi = 1 N H &lambda; Hi P Hi , t , u + &Sigma; Gk = 1 N G &lambda; Gk P Gk , t , u m , Then perform step 15.
Wherein, λ hifor current period Hi platform adjustable type Hydropower Unit participates in the state of adjusting, when current period Hi platform adjustable type Hydropower Unit participates in adjusting, λ hi=1; Current period Hi platform adjustable type Hydropower Unit has neither part nor lot in while adjusting, λ hi=0.P hi, t, uthe increase variable quantity of exerting oneself for current period Hi platform adjustable type Hydropower Unit.λ gkin the fired power generating unit put into operation for the current period, Gk platform fired power generating unit participates in the state of adjusting, when in the fired power generating unit that puts into operation of current period, Gk platform fired power generating unit participates in adjusting, and λ gK=1; In the fired power generating unit that puts into operation of current period, Gk platform fired power generating unit has neither part nor lot in while adjusting, λ gk=0.
Figure BDA0000393281780000262
for the increase variable quantity of exerting oneself of Gk platform fired power generating unit in the fired power generating unit put into operation under m of current period sampling background, N hfor the quantity of adjustable type Hydropower Unit, N gquantity for the adjustable type fired power generating unit that put into operation.
Sequentially adjusting the current period according to the second setting has put into operation and has still had the adjustable type generating set of adjustable capacity to be specially, first adjust the adjustable type Hydropower Unit that the current period put into operation and adjustable capacity is still arranged, when all adjustable type Hydropower Unit of having put into operation of current period do not have adjustable capacity, the gas turbine unit that order adjustment has from high to low put into operation by the power output regulations speed, when all gas turbine units of having put into operation of current period do not have adjustable capacity, the adjustable type fired power generating unit of shutting down that order adjustment from low to high puts into operation by minimum specific consumption, when all adjustable type fired power generating unit of shutting down of having put into operation of current period do not have adjustable capacity, the adjustable type fired power generating unit that order adjustment from low to high can not be shut down by minimum operating ratio consumption.
Step 14: sequentially adjust according to the second setting the adjustable type generating set that the current period put into operation and adjustable capacity is still arranged, until the adjustable type generating set that participates in adjusting meets
| &Delta;P | = &Sigma; Hi = 1 N H &lambda; Hi P Hi , t , d + &Sigma; Gk = 1 N G &lambda; Gk P Gk , t , d m .
Wherein, λ hifor current period Hi platform adjustable type Hydropower Unit participates in the state of adjusting, when current period Hi platform adjustable type Hydropower Unit participates in adjusting, λ hi=1; Current period Hi platform adjustable type Hydropower Unit has neither part nor lot in while adjusting, λ hi=0.P hi, t, dthe minimizing variable quantity of exerting oneself for current period Hi platform adjustable type Hydropower Unit.λ gkin the fired power generating unit put into operation for the current period, Gk platform fired power generating unit participates in the state of adjusting, when in the fired power generating unit that puts into operation of current period, Gk platform fired power generating unit participates in adjusting, and λ gK=1; In the fired power generating unit that puts into operation of current period, Gk platform fired power generating unit has neither part nor lot in while adjusting, λ gk=0.
Figure BDA0000393281780000263
for the minimizing variable quantity of exerting oneself of Gk platform fired power generating unit in the fired power generating unit put into operation under m of current period sampling background, N hfor the quantity of adjustable type Hydropower Unit, N gquantity for the adjustable type fired power generating unit that put into operation.
Step 15: the adjustable type fired power generating unit distribution load put into operation under each the sampling scene for the current period
Figure BDA0000393281780000271
calculate the Load Regulation abundant intensity interval of the adjustable type fired power generating unit put into operation under each sampling scene of current period simultaneously.
The detailed process of this step is identical with step 8, repeats no more herein.
Step 16: calculate the operating cost of the adjustable type fired power generating unit that each sampling of current period t put into operation under scene, select the adjustable type fired power generating unit that put into operation under sampling scene τ corresponding to the minimum value of operating cost and the load of distribution thereof
Figure BDA0000393281780000272
as the best of current period t adjustable type fired power generating unit scheme that puts into operation.Simultaneously, select the Load Regulation abundant intensity interval of the adjustable type fired power generating unit that put into operation under sampling scene τ corresponding to the minimum value of operating cost to regulate the abundant intensity interval as the optimum load of current period t adjustable type fired power generating unit.
The detailed process of this step is identical with step 9, repeats no more herein.
Step 17: judge t > whether N sets up, if t > N, and illustrate and finish whole dispatching cycle, now perform step 18; Otherwise, return to step 10, select put into operation scheme and optimum load of next the best constantly to regulate the abundant intensity interval.
Step 18: finish.
The above; be only the present invention's embodiment preferably, but protection scope of the present invention is not limited to this, anyly is familiar with in technical scope that those skilled in the art disclose in the present invention; the variation that can expect easily or replacement, within all should being encompassed in protection scope of the present invention.Therefore, protection scope of the present invention should be as the criterion with the protection range of claim.

Claims (9)

1. a multi-energy resource combined electric generating system is moved the optimization method of simulating, and it is characterized in that described method comprises:
Step 1: the configuration scheduling cycle also will be divided into T period described dispatching cycle;
Step 2: set each period sampling number of scenes N, and determine that each period wind-powered electricity generation unit and photoelectricity unit of each sampling scene blends the power value
Figure FDA0000393281770000011
wherein, m=1,2 ..., N, j=1,2 ..., T;
Step 3: according to formula P l1, j=P l, j-P f, jcalculate the load P that in each period system, the adjustable type generating set be should bear l1, j, and by P l1, jbe designated as the first duty value such as grade;
In described system, the adjustable type generating set comprises adjustable type Hydropower Unit and adjustable type fired power generating unit;
P l, jbe j the total load that the period system be should bear;
P f, jfor the load that in j period system, the generating set except the adjustable type generating set provides;
Step 4: according to formula
Figure FDA0000393281770000012
calculate the load P that in each period system, the adjustable type fired power generating unit be should bear l2, j, and by P l2, jbe designated as the second duty value such as grade; Simultaneously, calculate each period adjustable type Hydropower Unit accommodation limit value that makes progress
Figure FDA0000393281770000013
with downward accommodation limit value
Figure FDA0000393281770000014
μ hi, jbe the operating state of j period Hi adjustable type Hydropower Unit, when j period, Hi adjustable type Hydropower Unit put into operation, μ hi, j=1; When j period, Hi adjustable type Hydropower Unit do not put into operation, μ hi, j=0;
P hi, jbe power output and the P of j period Hi adjustable type Hydropower Unit hi, jhi* g * (H up, Hi, j-H down, Hi, j) * Q hi, j;
η hiit is the generating efficiency coefficient of Hi adjustable type Hydropower Unit;
G is gravity constant;
H up, Hi, jit is water level on the dam of j period Hi adjustable type Hydropower Unit;
H down, Hi, jit is the level of tail water of j period Hi adjustable type Hydropower Unit;
Q hi, jit is the generating average discharge of j period Hi adjustable type Hydropower Unit;
N hquantity for the adjustable type Hydropower Unit;
μ hp, jbe the state that draws water of j period pumped storage machine, when j the period, pumped storage machine put into operation, μ hp, j=1, when j the period, pumped storage machine did not put into operation, μ hp, j=0;
P hp, jbe the power that draws water of j period pumped storage machine;
Step 5: the state that initially puts into operation of determining the adjustable type fired power generating unit; Specifically, calculate respectively the duty value minimum value such as second of the second duty value maximums such as grade of front 5 periods of dispatching cycle and front 5 periods, sequentially drop into the adjustable type fired power generating unit until the total installation of generating capacity of the adjustable type fired power generating unit put into operation is greater than the second maximum such as duty value such as grade of front 5 periods of dispatching cycle according to the first setting, and the minimum technology of the adjustable type fired power generating unit put into operation the summation of exerting oneself is less than the duty value minimum values such as second of front 5 periods of dispatching cycle;
Described adjustable type fired power generating unit comprises the adjustable type fired power generating unit that can not shut down, adjustable type fired power generating unit and the gas turbine unit that can shut down;
Step 6: set initial period t=5;
Step 7: according to formula
Figure FDA0000393281770000021
calculate the load that current period t adjustable type fired power generating unit be should bear under each sampling scene
Figure FDA0000393281770000022
and will
Figure FDA0000393281770000023
be designated as the C grade duty value;
Step 8: the adjustable type fired power generating unit distribution load put into operation under each the sampling scene for current period t
Figure FDA0000393281770000031
determine the Load Regulation abundant intensity interval of the adjustable type fired power generating unit put into operation under each sampling scene of current period t simultaneously;
Step 9: calculate the operating cost of the adjustable type fired power generating unit that each sampling of current period t put into operation under scene, select the adjustable type fired power generating unit that put into operation under sampling scene τ corresponding to the minimum value of operating cost and the load of distribution thereof
Figure FDA0000393281770000032
as the best of current period t adjustable type fired power generating unit scheme that puts into operation; Simultaneously, select the Load Regulation abundant intensity interval of the adjustable type fired power generating unit that put into operation under sampling scene τ corresponding to the minimum value of operating cost to regulate the abundant intensity interval as the optimum load of current period t adjustable type fired power generating unit;
Step 10: make t=t+1, calculate the duty value maximum such as second of continuous 4 periods before current period and current period, and the duty value minimum value such as second of continuous 4 periods before current period and current period;
If the total installation of generating capacity that the adjustable type fired power generating unit of the previous period of current period puts into operation is less than the duty value maximum such as second of continuous 4 periods before current period and current period, according to the first setting, sequentially drop into the adjustable type fired power generating unit also do not put into operation, until the total installation of generating capacity of the adjustable type fired power generating unit put into operation is greater than the duty value maximum such as second of continuous 4 periods before current period and current period;
The summation if the minimum technology of adjustable type fired power generating unit that the duty value maximum such as second of continuous 4 periods is less than the previous period of current period before current period and current period is exerted oneself, the adjustable type fired power generating unit that stops having put into operation according to the inverted order of the first setting order, until the minimum technology of the adjustable type fired power generating unit put into operation is exerted oneself, summation is less than the duty value minimum value such as second of continuous 4 periods before current period and current period;
Step 11: judge that whether the adjustable type fired power generating unit that adjustable type fired power generating unit that current time puts into operation puts into operation with the previous moment of current time is identical, if the adjustable type fired power generating unit that the adjustable type fired power generating unit that current time puts into operation puts into operation with the previous moment of current time is identical, perform step 12; Otherwise, perform step 15;
Step 12: calculate adjustable type fired power generating unit load variations amount under each sampling scene of current period t
Figure FDA0000393281770000041
if Δ P>0, perform step 13; If Δ P<0, perform step 14; If Δ P=0, perform step 15;
Step 13: sequentially adjust according to the second setting the adjustable type generating set that the current period put into operation and adjustable capacity is still arranged, until the adjustable type generating set that participates in adjusting meets &Delta;P = &Sigma; Hi = 1 N H &lambda; Hi P Hi , t , u + &Sigma; Gk = 1 N G &lambda; Gk P Gk , t , u m ; Then perform step 15;
Wherein, λ hifor current period Hi platform adjustable type Hydropower Unit participates in the state of adjusting, when current period Hi platform adjustable type Hydropower Unit participates in adjusting, λ hi=1; Current period Hi platform adjustable type Hydropower Unit has neither part nor lot in while adjusting, λ hi=0;
P hi, t, uthe increase variable quantity of exerting oneself for current period Hi platform adjustable type Hydropower Unit;
λ gkin the fired power generating unit put into operation for the current period, Gk platform fired power generating unit participates in the state of adjusting, when in the fired power generating unit that puts into operation of current period, Gk platform fired power generating unit participates in adjusting, and λ gK=1; In the fired power generating unit that puts into operation of current period, Gk platform fired power generating unit has neither part nor lot in while adjusting, λ gk=0;
for the increase variable quantity of exerting oneself of Gk platform fired power generating unit in the fired power generating unit put into operation under m of current period sampling background;
N hquantity for the adjustable type Hydropower Unit;
N gquantity for the adjustable type fired power generating unit that put into operation;
Step 14: sequentially adjust according to the second setting the adjustable type generating set that the current period put into operation and adjustable capacity is still arranged, until the adjustable type generating set that participates in adjusting meets
| &Delta;P | = &Sigma; Hi = 1 N H &lambda; Hi P Hi , t , d + &Sigma; Gk = 1 N G &lambda; Gk P Gk , t , d m ;
Wherein, λ hifor current period Hi platform adjustable type Hydropower Unit participates in the state of adjusting, when current period Hi platform adjustable type Hydropower Unit participates in adjusting, λ hi=1; Current period Hi platform adjustable type Hydropower Unit has neither part nor lot in while adjusting, λ hi=0;
P hi, t, dthe minimizing variable quantity of exerting oneself for current period Hi platform adjustable type Hydropower Unit;
λ gkin the fired power generating unit put into operation for the current period, Gk platform fired power generating unit participates in the state of adjusting, when in the fired power generating unit that puts into operation of current period, Gk platform fired power generating unit participates in adjusting, and λ gK=1; In the fired power generating unit that puts into operation of current period, Gk platform fired power generating unit has neither part nor lot in while adjusting, λ gk=0;
Figure FDA0000393281770000052
for the minimizing variable quantity of exerting oneself of Gk platform fired power generating unit in the fired power generating unit put into operation under m of current period sampling background;
N hquantity for the adjustable type Hydropower Unit;
N gquantity for the adjustable type fired power generating unit that put into operation;
Step 15: the adjustable type fired power generating unit distribution load put into operation under each the sampling scene for the current period
Figure FDA0000393281770000054
calculate the Load Regulation abundant intensity interval of the adjustable type fired power generating unit put into operation under each sampling scene of current period simultaneously;
Step 16: calculate the operating cost of the adjustable type fired power generating unit that each sampling of current period t put into operation under scene, select the adjustable type fired power generating unit that put into operation under sampling scene τ corresponding to the minimum value of operating cost and the load of distribution thereof
Figure FDA0000393281770000053
as the best of current period t adjustable type fired power generating unit scheme that puts into operation; Simultaneously, select the Load Regulation abundant intensity interval of the adjustable type fired power generating unit that put into operation under sampling scene τ corresponding to the minimum value of operating cost to regulate the abundant intensity interval as the optimum load of current period t adjustable type fired power generating unit;
Step 17: judge t > whether N sets up, if t > N, perform step 18; Otherwise, return to step 10;
Step 18: finish.
2. method according to claim 1, is characterized in that wind-powered electricity generation unit and the photoelectricity unit of described each sampling scene of definite each period blends the power value comprise following sub-step:
Sub-step 101: respectively according to the historical data of wind-powered electricity generation unit and photoelectricity unit output, make the histogram of wind-powered electricity generation unit output prediction deviation and the histogram of photoelectricity unit output prediction deviation;
Sub-step 102: by approximating method, obtain respectively each period wind-powered electricity generation unit output prediction deviation probability density function and photoelectricity unit output prediction deviation probability density function;
Sub-step 103: each period wind-powered electricity generation unit output prediction deviation probability density function and photoelectricity unit output prediction deviation probability density function are carried out respectively to integral operation, obtain wind-powered electricity generation unit output prediction deviation probability density function and the photoelectricity unit output prediction deviation probability density function of each period;
Sub-step 104: utilize formula
Figure FDA0000393281770000062
calculate the basic sampled value under each sampling scene of each period; Wherein, m=1,2 ..., N, N is each period sampling number of scenes, U is the random number on interval (0,1);
Sub-step 105: utilize formula
Figure FDA0000393281770000063
calculate each sampling scene leeward group of motors of each period prediction deviation of exerting oneself, and utilize formula
Figure FDA0000393281770000064
calculate photoelectricity unit output prediction deviation under each sampling scene of each period; Wherein,
Figure FDA0000393281770000065
for the inverse function of wind-powered electricity generation unit output prediction deviation probability density function,
Figure FDA0000393281770000066
inverse function for photoelectricity unit output prediction deviation probability density function; J=1,2 ..., T, the time hop count that T is dispatching cycle;
Sub-step 106: according to formula
Figure FDA0000393281770000067
the wind-powered electricity generation unit and the photoelectricity unit that calculate each sampling scene of each period blend the power value;
Wherein, P w, jfor each period wind-powered electricity generation unit output predicted value; P p, jfor each period photoelectricity unit output predicted value.
3. method according to claim 1, is characterized in that each period adjustable type Hydropower Unit of described calculating accommodation limit value that makes progress
Figure FDA0000393281770000071
adopt formula
Figure FDA0000393281770000072
Wherein, be the upwards pondage of j period Hi adjustable type Hydropower Unit, and
P Hi , j op = P Hi , j max - P Hi , j ;
Figure FDA0000393281770000075
it is the peak power output of j period Hi adjustable type Hydropower Unit;
P hi, jit is the power output of j period Hi adjustable type Hydropower Unit;
μ hi, jbe the operating state of j period Hi adjustable type Hydropower Unit, when j period, Hi adjustable type Hydropower Unit put into operation, μ hi, j=1; When j period, Hi adjustable type Hydropower Unit do not put into operation, μ hi, j=0.
4. method according to claim 1, is characterized in that the downward accommodation limit value of each period adjustable type Hydropower Unit of described calculating
Figure FDA0000393281770000076
adopt formula
Figure FDA0000393281770000077
Wherein, be the downward pondage of j period Hi adjustable type Hydropower Unit, and
P Hi , j ne = P Hi , j - P Hi , j min ;
Figure FDA00003932817700000710
the minimum that is j period Hi adjustable type Hydropower Unit is forced power output;
P hi, jit is the power output of j period Hi adjustable type Hydropower Unit;
μ hi, jbe the operating state of j period Hi adjustable type Hydropower Unit, when j period, Hi adjustable type Hydropower Unit put into operation, μ hi, j=1; When j period, Hi adjustable type Hydropower Unit do not put into operation, μ hi, j=0.
5. method according to claim 1; it is characterized in that describedly according to the first setting, sequentially dropping into the adjustable type fired power generating unit specifically; first by minimum operating ratio consumption, order from low to high drops into the adjustable type fired power generating unit that can not shut down; by minimum specific consumption, order from low to high drops into the adjustable type fired power generating unit that can shut down again, and finally by the power output regulations speed, order from high to low drops into gas turbine unit.
6. method according to claim 1, is characterized in that the adjustable type fired power generating unit distribution load put into operation under described each sampling scene for current period t
Figure FDA0000393281770000081
specifically comprise following sub-step:
Sub-step 201: set the adjustable type fired power generating unit put into operation the out-of-limit sign nFlag=1 that exerts oneself;
Sub-step 202: according to formula calculate the micro-gaining rate of consumption of fired power generating unit; Wherein,
Figure FDA0000393281770000084
for C grade duty value and m=1,2 ..., N;
N is each period sampling number of scenes;
N ' gfor current period t put into operation and the exerted oneself quantity of the adjustable type fired power generating unit that do not transfinite;
A gk 'for current period t the Gk ' platform put into operation and the exerted oneself quadratic term coefficient of consumption characteristic quadratic function of the adjustable type fired power generating unit that do not transfinite;
B gk 'for current period t the Gk ' platform put into operation and the exerted oneself Monomial coefficient of consumption characteristic quadratic function of the adjustable type fired power generating unit that do not transfinite;
Sub-step 203: make Gk=1;
Sub-step 204: according to formula
Figure FDA0000393281770000083
calculate adjustable type fired power generating unit value of exerting oneself of having put into operation of Gk platform under each sampling scene of current period t;
A gkthe quadratic term coefficient of the consumption characteristic quadratic function of the adjustable type fired power generating unit put into operation for current period t Gk platform;
B gkthe Monomial coefficient of the consumption characteristic quadratic function of the adjustable type fired power generating unit put into operation for current period t Gk platform;
Sub-step 205: if
Figure FDA0000393281770000091
transfinite, order
Figure FDA0000393281770000092
maximum output value during for this adjustable type fired power generating unit put into operation not out-of-limit, and make nFlag=0;
Sub-step 206: judgement Gk>N gwhether set up, if Gk N g, carry out sub-step 207; Otherwise, make Gk=Gk+1, return to sub-step 204;
N gthe quantity of the adjustable type fired power generating unit put into operation for current period t;
Sub-step 207: judge whether nFlag=0 sets up, if nFlag=0 returns to sub-step 201; Otherwise, carry out sub-step 208;
Sub-step 208: make the load that under each sampling scene of current period t, the adjustable type fired power generating unit that put into operation of Gk platform provides be
Figure FDA0000393281770000093
7. method according to claim 1, is characterized in that the operating cost of the adjustable type fired power generating unit that put into operation under each sampling scene of the current period t of described calculating adopts formula
F cos t ( P G , t m ) = &Sigma; Gk = 1 N T [ &mu; Gk , t m &times; f k ( P Gk , t m ) + &mu; Gk , t m ( 1 - &mu; Gk , t - 1 m ) S Gk ] ;
Wherein, the operating cost of the adjustable type fired power generating unit put into operation under m the sampling scene for current period t, m=1,2 ..., N, N is the sampling number of scenes in each period;
N tfor adjustable type fired power generating unit quantity;
Figure FDA0000393281770000096
for the current period t Gk platform adjustable type fired power generating unit state that puts into operation, current period t Gk platform adjustable type fired power generating unit puts into operation,
Figure FDA0000393281770000101
current period t Gk platform adjustable type fired power generating unit does not put into operation,
Figure FDA0000393281770000103
for the consumption characteristic quadratic function of current period t Gk platform adjustable type fired power generating unit, and
f k ( P Gk , t m ) = A Gk &times; ( P Gk , t m ) 2 + B Gk &times; P Gk , t m + C Gk ;
A gkquadratic term coefficient for the consumption characteristic quadratic function of current period t Gk platform adjustable type fired power generating unit;
B gkmonomial coefficient for the consumption characteristic quadratic function of current period t Gk platform adjustable type fired power generating unit;
C gkconstant term coefficient for the consumption characteristic quadratic function of current period t Gk platform adjustable type fired power generating unit;
S gkfor the payment for initiation of current period t Gk platform adjustable type fired power generating unit is used.
8. method according to claim 1, is characterized in that the Load Regulation abundant intensity interval of the adjustable type fired power generating unit that put into operation under each sampling scene of described current period t is
[ P L 3 , t m - P H , t ne - P G , t ne , m , P L 3 , t m + P H , t op + P G , t op , m ] ;
Wherein,
Figure FDA0000393281770000106
for the downward accommodation limit value of current period t adjustable type Hydropower Unit;
Figure FDA0000393281770000107
for adjustable type thermal power unit operation lower bound under each sampling scene of current period t, and
P G , t ne , m = &Sigma; Gk = 1 N T &mu; Gk , t m P Gk , t ne , m ;
Figure FDA0000393281770000109
for the operation lower bound of current period t Gk platform adjustable type fired power generating unit, and
P Gk , t ne , m = min ( P Gk , t m - P Gk , min , - R Gk , d &times; &Delta;t ) ;
P gk, minit is the Gk platform adjustable type fired power generating unit lower limit of exerting oneself;
R gk, dit is the downward creep speed of Gk platform adjustable type fired power generating unit;
Figure FDA0000393281770000111
for the current period t Gk platform adjustable type fired power generating unit state that puts into operation, current period t Gk platform adjustable type fired power generating unit puts into operation, current period t Gk platform adjustable type fired power generating unit does not put into operation,
Figure FDA0000393281770000113
The time span that Δ t is current period t;
N tquantity for the adjustable type fired power generating unit;
for the current period t adjustable type Hydropower Unit accommodation limit value that makes progress;
Figure FDA0000393281770000115
for the adjustable type thermal power unit operation upper bound under each sampling scene of current period t, and
P G , t op , m = &Sigma; Gk = 1 N T &mu; Gk , t m P Gk , t op , m ;
Figure FDA0000393281770000117
for the current period t Gk platform adjustable type thermal power unit operation upper bound, and
P Gk , t op , m = min ( P Gk , max - P Gk , t m , - R Gk , u &times; &Delta;t ) ;
P gk, maxit is the upper limit of exerting oneself of Gk platform adjustable type fired power generating unit;
R gk, uit is the Gk platform adjustable type fired power generating unit creep speed that makes progress;
M=1,2 ..., N, N is the sampling number of scenes in each period.
9. method according to claim 1, it is characterized in that describedly sequentially adjusting the current period according to the second setting and having put into operation and still had the adjustable type generating set of adjustable capacity to be specially, first adjust the adjustable type Hydropower Unit that the current period put into operation and adjustable capacity is still arranged, when all adjustable type Hydropower Unit of having put into operation of current period do not have adjustable capacity, the gas turbine unit that order adjustment has from high to low put into operation by the power output regulations speed, when all gas turbine units of having put into operation of current period do not have adjustable capacity, the adjustable type fired power generating unit of shutting down that order adjustment from low to high puts into operation by minimum specific consumption, when all adjustable type fired power generating unit of shutting down of having put into operation of current period do not have adjustable capacity, the adjustable type fired power generating unit that order adjustment from low to high can not be shut down by minimum operating ratio consumption.
CN201310469621.XA 2013-10-10 2013-10-10 Optimization method for multi-energy combined power generation system operation simulation Expired - Fee Related CN103490449B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201310469621.XA CN103490449B (en) 2013-10-10 2013-10-10 Optimization method for multi-energy combined power generation system operation simulation

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201310469621.XA CN103490449B (en) 2013-10-10 2013-10-10 Optimization method for multi-energy combined power generation system operation simulation

Publications (2)

Publication Number Publication Date
CN103490449A true CN103490449A (en) 2014-01-01
CN103490449B CN103490449B (en) 2015-04-15

Family

ID=49830481

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201310469621.XA Expired - Fee Related CN103490449B (en) 2013-10-10 2013-10-10 Optimization method for multi-energy combined power generation system operation simulation

Country Status (1)

Country Link
CN (1) CN103490449B (en)

Cited By (13)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103887826A (en) * 2014-03-05 2014-06-25 中国东方电气集团有限公司 Active power distribution method for achieving least fuel consumption in system of a plurality of conventional generating sets
CN103956742A (en) * 2014-04-14 2014-07-30 国家电网公司 Determining method of operating mode of electric power system under complex energy source environment
CN104124707A (en) * 2014-07-23 2014-10-29 北京四方继保自动化股份有限公司 Hydropower station power quick adjustment system and implementation method based on intensive small hydropower station groups
CN104156886A (en) * 2014-08-11 2014-11-19 上海交通大学 Method of evaluating flexibility of power supply with renewable energy power system
CN104467029A (en) * 2014-11-25 2015-03-25 国家电网公司 Regional small hydropower and wind power day combined optimization and dispatching method
CN104600747A (en) * 2015-01-21 2015-05-06 西安交通大学 Operation optimizing method capable of coordinating operation risk and wind energy consumption of power system
CN105048516A (en) * 2015-08-18 2015-11-11 四川大学 Wind-light-water-fire multi-source complementary optimization scheduling method
CN106130079A (en) * 2016-07-08 2016-11-16 华北电力大学 A kind of edema due to wind pathogen fire short-term joint optimal operation method
CN107482656A (en) * 2017-08-11 2017-12-15 电力规划总院有限公司 Power planning method, apparatus, equipment and the computer-readable recording medium of micro-capacitance sensor
CN107634547A (en) * 2017-11-22 2018-01-26 国家电网公司 Contributed based on new energy and predict that the electric association system of error goes out electric control method
CN108258732A (en) * 2018-01-26 2018-07-06 国网辽宁省电力有限公司 Control method of the New-energy power system under the method for operation of normal regulation domain
CN109936170A (en) * 2019-04-08 2019-06-25 东北电力大学 Consider the honourable extreme misery complementation coordination optimization dispatching method of power supply flexibility nargin
CN110492478A (en) * 2019-08-20 2019-11-22 云南电网有限责任公司 A kind of water power and wind-powered electricity generation combine power producing characteristics analysis method

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102738835A (en) * 2012-07-17 2012-10-17 太原电力高等专科学校 Wind-fire-water co-scheduling method on basis of multi-agent system
WO2013013656A2 (en) * 2011-07-22 2013-01-31 ct.e Controltechnology Engineering GmbH Method for adaptively controlling the supply voltage in local distribution systems
CN103066626A (en) * 2012-10-31 2013-04-24 中国电力科学研究院 Multi-source coordinating day-ahead generation scheduling method
CN103151803A (en) * 2013-03-14 2013-06-12 吉林省电力有限公司电力科学研究院 Method for optimizing wind power system-contained unit and backup configuration

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2013013656A2 (en) * 2011-07-22 2013-01-31 ct.e Controltechnology Engineering GmbH Method for adaptively controlling the supply voltage in local distribution systems
CN102738835A (en) * 2012-07-17 2012-10-17 太原电力高等专科学校 Wind-fire-water co-scheduling method on basis of multi-agent system
CN103066626A (en) * 2012-10-31 2013-04-24 中国电力科学研究院 Multi-source coordinating day-ahead generation scheduling method
CN103151803A (en) * 2013-03-14 2013-06-12 吉林省电力有限公司电力科学研究院 Method for optimizing wind power system-contained unit and backup configuration

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
刘可真等: "一种实用的机组最优启停计划方法", 《昆明理工大学学报》 *

Cited By (21)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103887826A (en) * 2014-03-05 2014-06-25 中国东方电气集团有限公司 Active power distribution method for achieving least fuel consumption in system of a plurality of conventional generating sets
CN103956742B (en) * 2014-04-14 2015-12-02 国家电网公司 The defining method of power system operation mode under a kind of complicated energy environment
CN103956742A (en) * 2014-04-14 2014-07-30 国家电网公司 Determining method of operating mode of electric power system under complex energy source environment
CN104124707A (en) * 2014-07-23 2014-10-29 北京四方继保自动化股份有限公司 Hydropower station power quick adjustment system and implementation method based on intensive small hydropower station groups
CN104156886A (en) * 2014-08-11 2014-11-19 上海交通大学 Method of evaluating flexibility of power supply with renewable energy power system
CN104156886B (en) * 2014-08-11 2017-05-24 上海交通大学 Method of evaluating flexibility of power supply with renewable energy power system
CN104467029A (en) * 2014-11-25 2015-03-25 国家电网公司 Regional small hydropower and wind power day combined optimization and dispatching method
CN104600747A (en) * 2015-01-21 2015-05-06 西安交通大学 Operation optimizing method capable of coordinating operation risk and wind energy consumption of power system
CN104600747B (en) * 2015-01-21 2016-09-28 西安交通大学 The Operation of Electric Systems optimization method that coordinated operation risk and wind energy are dissolved
CN105048516A (en) * 2015-08-18 2015-11-11 四川大学 Wind-light-water-fire multi-source complementary optimization scheduling method
CN106130079B (en) * 2016-07-08 2018-11-09 华北电力大学 A kind of short-term joint optimal operation method of geomantic omen fire
CN106130079A (en) * 2016-07-08 2016-11-16 华北电力大学 A kind of edema due to wind pathogen fire short-term joint optimal operation method
CN107482656A (en) * 2017-08-11 2017-12-15 电力规划总院有限公司 Power planning method, apparatus, equipment and the computer-readable recording medium of micro-capacitance sensor
CN107482656B (en) * 2017-08-11 2019-12-27 电力规划总院有限公司 Power planning method, device and equipment for microgrid and computer readable storage medium
CN107634547A (en) * 2017-11-22 2018-01-26 国家电网公司 Contributed based on new energy and predict that the electric association system of error goes out electric control method
CN108258732A (en) * 2018-01-26 2018-07-06 国网辽宁省电力有限公司 Control method of the New-energy power system under the method for operation of normal regulation domain
CN108258732B (en) * 2018-01-26 2022-06-21 国网辽宁省电力有限公司 Control method of new energy power system in normal regulation and control domain operation mode
CN109936170A (en) * 2019-04-08 2019-06-25 东北电力大学 Consider the honourable extreme misery complementation coordination optimization dispatching method of power supply flexibility nargin
CN109936170B (en) * 2019-04-08 2022-02-18 东北电力大学 Wind, light, water and fire complementary coordination optimization scheduling method considering power supply flexibility margin
CN110492478A (en) * 2019-08-20 2019-11-22 云南电网有限责任公司 A kind of water power and wind-powered electricity generation combine power producing characteristics analysis method
CN110492478B (en) * 2019-08-20 2023-04-11 云南电网有限责任公司 Water, electricity and wind electricity combined output characteristic analysis method

Also Published As

Publication number Publication date
CN103490449B (en) 2015-04-15

Similar Documents

Publication Publication Date Title
CN103490449A (en) Optimization method for multi-energy combined power generation system operation simulation
CN112104007B (en) Centralized control method for generalized source storage system scheduling
CN109767105B (en) Wind, light, water, fire and storage combined system-based multi-energy complementary coordinated power generation scheduling method
CN103580063B (en) A kind of method of large-scale grid connection wind-powered electricity generation of dissolving based on demanding party&#39;s response
CN107248751A (en) A kind of energy storage station dispatch control method for realizing distribution network load power peak load shifting
CN105140936B (en) A kind of micro-capacitance sensor frequency modulation control method for minimizing cost of electricity-generating
CN115018260A (en) Peak regulation optimal scheduling method of cascade water-light complementary power generation system
CN108039737B (en) Source-grid-load coordinated operation simulation system
CN109038689B (en) Ultra-short-term optimal scheduling method for power system
CN105225022A (en) A kind of economy optimizing operation method of cogeneration of heat and power type micro-capacitance sensor
CN108521132B (en) Multi-time scale optimization control method for frequency adjustment of multi-energy complementary support power grid
CN109543889A (en) A kind of regional complex energy resource system cooperates with optimizing operation method a few days ago
CN104167765A (en) Admitting ability distribution-based maximum wind power installed capacity calculation method
CN103050998B (en) Thermal power system dynamic scheduling method of wind power integration
CN107730032B (en) Universal energy network electricity-cooling-heating triple co-generation optimized scheduling system and method
CN112134271A (en) Installed capacity optimization method and system for multi-energy complementary system
CN116979578A (en) Electric and thermal triple generation optimal scheduling method and system for wind, light, water and fire storage
CN112134307A (en) Electric energy storage/heat energy storage capacity joint optimization method and system of multi-energy complementary system
CN114465226A (en) Method for establishing multi-level standby acquisition joint optimization model of power system
CN117060396A (en) Day-ahead optimal operation method of wind-solar-fire-storage multi-energy power system
CN115936336B (en) Virtual power plant capacity configuration and regulation operation optimization method
CN115577549B (en) Virtual power plant closed-loop direct control method containing multiple types of resources
CN116562554A (en) Long-short-period energy storage joint planning method considering flexible supply-demand balance
CN114357725A (en) Source-to-load double-end uncertainty modeling method considering carbon capture emission
CN109617135B (en) Power scheduling multi-objective decision method for hybrid power generation system

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
C14 Grant of patent or utility model
GR01 Patent grant
CF01 Termination of patent right due to non-payment of annual fee
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20150415

Termination date: 20181010