CN109245180A - A kind of scene fire storage coordination optimization operation method - Google Patents

A kind of scene fire storage coordination optimization operation method Download PDF

Info

Publication number
CN109245180A
CN109245180A CN201811376350.2A CN201811376350A CN109245180A CN 109245180 A CN109245180 A CN 109245180A CN 201811376350 A CN201811376350 A CN 201811376350A CN 109245180 A CN109245180 A CN 109245180A
Authority
CN
China
Prior art keywords
power
cogeneration units
simulation
type cogeneration
constraint
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
CN201811376350.2A
Other languages
Chinese (zh)
Other versions
CN109245180B (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.)
Xian Jiaotong University
Original Assignee
Xian Jiaotong University
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 Xian Jiaotong University filed Critical Xian Jiaotong University
Priority to CN201811376350.2A priority Critical patent/CN109245180B/en
Publication of CN109245180A publication Critical patent/CN109245180A/en
Application granted granted Critical
Publication of CN109245180B publication Critical patent/CN109245180B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • H—ELECTRICITY
    • H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
    • H02J3/00—Circuit arrangements for ac mains or ac distribution networks
    • H02J3/38—Arrangements for parallely feeding a single network by two or more generators, converters or transformers
    • H02J3/46—Controlling of the sharing of output between the generators, converters, or transformers
    • H—ELECTRICITY
    • H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
    • H02J3/00—Circuit arrangements for ac mains or ac distribution networks
    • H02J3/28—Arrangements for balancing of the load in a network by storage of energy
    • H—ELECTRICITY
    • H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
    • H02J2203/00—Indexing scheme relating to details of circuit arrangements for AC mains or AC distribution networks
    • H02J2203/20—Simulating, e g planning, reliability check, modelling or computer assisted design [CAD]
    • Y—GENERAL 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
    • Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02E—REDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
    • Y02E10/00—Energy generation through renewable energy sources
    • Y02E10/50—Photovoltaic [PV] energy
    • Y02E10/56—Power conversion systems, e.g. maximum power point trackers
    • Y—GENERAL 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
    • Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02E—REDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
    • Y02E10/00—Energy generation through renewable energy sources
    • Y02E10/70—Wind energy
    • Y02E10/76—Power conversion electric or electronic aspects
    • Y—GENERAL 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
    • Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02E—REDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
    • Y02E70/00—Other energy conversion or management systems reducing GHG emissions
    • Y02E70/30—Systems combining energy storage with energy generation of non-fossil origin

Landscapes

  • Engineering & Computer Science (AREA)
  • Power Engineering (AREA)
  • Supply And Distribution Of Alternating Current (AREA)

Abstract

A kind of scene fire storage coordination optimization operation method, fired power generating unit is obtained first, back pressure type cogeneration units and steam-extracting type cogeneration units open state sequence, secondly Wind turbines and photovoltaic unit available power sequences are obtained, then using Operation of Electric Systems cost minimization as target, consider electric system realtime power Constraints of Equilibrium, power of the assembling unit constraint, fired power generating unit Climing constant, cogeneration units are for thermal confinement, hydroenergy storage station constraint and electric system Reserve Constraint, establish scene fire storage coordinated operation model, finally, scene fire storage coordinated operation model is solved to POWER SYSTEM STATE in the different simulation periods, obtain power of the assembling unit arrangement；If feasible solution is not present in model, introduces abandonment power and abandon optical power, solved again after storing up coordinated operation model modification to honourable fire；The scene fire storage coordination optimization operation that the present invention is realized with Operation of Electric Systems cost minimization, can provide guidance for Operation of Electric Systems and scheduling.

Description

A kind of scene fire storage coordination optimization operation method
Technical field
The invention belongs to the electric system new energy of electric system new energy consumption field, the especially power supply containing polymorphic type to disappear Receive problem, and in particular to a kind of scene fire storage coordination optimization operation method.
Background technique
In recent years, generation of electricity by new energy is quickly grown, and proportion constantly increases in systems for new energy installed capacity.But by Have the characteristics that fluctuation, randomness and intermittent in generation of electricity by new energy, new energy large-scale grid connection gives system regulating power band Carry out huge challenge, causes fired power generating unit frequent start-stop and abandonment, abandons the problems such as light.Pump-storage generator is set as a kind of energy storage It is standby, it is played an important role in terms of system peak load shifting, frequency modulation and black starting-up, pump-storage generator combines fortune with wind-power electricity generation Row, it has also become solve the problems, such as that new energy dissolves important one of strategy in electric system.
In the northern area of China, the run-limiting of heat supply in winter phase cogeneration units " electricity determining by heat " peak regulation of unit Ability brings problem for new energy consumption.More power supplys and the coordinated operation model deposited do not consider cogeneration units to electricity at present The limitation of Force system bring, and existing electric system multi-source coordinated operation research only considers Wind turbines and the storage that draws water mostly Energy power station coordinated operation, with the continuous increase of photovoltaic plant installed capacity, photovoltaic power generation quantity ratio shared in total power generation Be continuously improved again, it is therefore necessary to study wind power plant, photovoltaic plant, thermal power plant and hydroenergy storage station coordinated operation power train The system method of operation.
Summary of the invention
For the deficiency of the existing more power supply coordinated operation methods of electric system, it is an object of the invention to propose a kind of scene Fire storage coordination optimization operation method, the present invention are directed to more power supplys electricity containing wind-powered electricity generation, photovoltaic, water-storage, conventional thermoelectricity and thermoelectricity Force system considers disappearing for the operation characteristic of hydroenergy storage station, the coupled thermomechanics relationship of cogeneration units and wind-powered electricity generation and photovoltaic Situation is received, to fired power generating unit power, back pressure type cogeneration units output power and heating power, steam-extracting type cogeneration units Output power and heating power, hydroenergy storage station draw water power or generated output is arranged, and realize with Operation of Electric Systems The scene fire storage coordination optimization operation of cost minimization, can provide guidance for Operation of Electric Systems and scheduling.
In order to achieve the above object, the present invention adopts the following technical scheme:
A kind of scene fire storage coordination optimization operation method, comprising the following steps:
1) input initial data, including electric system electric load data, thermic load data, fired power generating unit forced outage rate and Mean repair time, back pressure type cogeneration units forced outage rate and mean repair time, steam-extracting type cogeneration units are strong Compel outage rate and mean repair time, historical wind speed data and history solar irradiance data；Define maximum analog periodicity Dmax With number of segment T when maximum analog in a simulation cyclemax, and simulation cycle D=1 is set, simulate period T=1；
2) fired power generating unit open state sequence in simulation cycle D, back pressure type heat are obtained based on state duration sampling method Electricity Federation produces unit open state sequence and steam-extracting type cogeneration units open state sequence；
3) Wind turbines available power sequences and photovoltaic unit available power sequence in simulation cycle D are obtained based on inverse transformation method Column；
4) obtain simulation cycle D simulation period T in POWER SYSTEM STATE, including electric load data, thermic load data, Fired power generating unit open state, back pressure type cogeneration units open state and steam-extracting type cogeneration units open state, wind-powered electricity generation Unit available power and photovoltaic unit available power；
5) scene fire storage coordinated operation model is solved to the POWER SYSTEM STATE in step 4), if model there are feasible solution, Then illustrate that electric system can dissolve Wind turbines available power completely in the simulation period T of simulation cycle D and photovoltaic unit can With power, enter step 8)；If feasible solution is not present in model, illustrate that electric system cannot be complete in the simulation period T of simulation cycle D It totally disappeared and receive Wind turbines available power and photovoltaic unit available power, enter step 6)；Honourable fire storage coordinated operation model is by mesh Scalar functions and constraint condition composition, concrete model are as follows:
5.1) objective function
The objective function of honourable fire storage coordinated operation model is Operation of Electric Systems cost minimization in simulation cycle D, electric power System operation cost includes thermal power unit operation cost, back pressure type cogeneration units operating cost, steam-extracting type cogeneration of heat and power machine Group operating cost and pump-storage generator operating cost, specifically:
Operating cost of the fired power generating unit i in simulation period T are as follows:
Operating cost of the back pressure type cogeneration units m in simulation period T are as follows:
Operating cost of the steam-extracting type cogeneration units n in simulation period T are as follows:
Hydroenergy storage station operating cost are as follows:
In formula: NthFor fired power generating unit number,For back pressure type cogeneration units number,For steam-extracting type cogeneration of heat and power Unit number, NpFor hydroenergy storage station number,The power in period T is being simulated for fired power generating unit i,For back pressure type thermoelectricity Output power of the coproduction unit m in simulation period T,For confession of the steam-extracting type cogeneration units n in simulation period T Electrical power,For heating power of the steam-extracting type cogeneration units n in simulation period T, ai,biAnd ciFor fired power generating unit Consumption characterisitic parameter,For back pressure type cogeneration units consumption characterisitic parameter,For steam-extracting type cogeneration units consumption characterisitic parameter,It is opened for hydroenergy storage station u The primary cost of machine,Primary cost is shut down for hydroenergy storage station u,It is hydroenergy storage station u in the simulation period The operating status indicator variable that draws water in T,Refer to for generator operation state of the hydroenergy storage station u in simulation period T+1 Show variable, when hydroenergy storage station u is in simulation period T draws water operating status,Value is 1,Value is 0；Work as pumping When water accumulation of energy unit plant u is in generator operation state in simulation period T,Value is 1,Value is 0；；
5.2) constraint condition
The constraint condition of honourable fire storage coordinated operation model include electric system realtime power Constraints of Equilibrium, the power of the assembling unit about Beam, fired power generating unit Climing constant, cogeneration units for thermal confinement, hydroenergy storage station constraint and electric system Reserve Constraint, Specifically:
5.2.1) electric system realtime power Constraints of Equilibrium
Electric system realtime power Constraints of Equilibrium includes real-time electrical power Constraints of Equilibrium and the constraint of real-time heating power balance；；
Real-time electrical power Constraints of Equilibrium are as follows:
Real-time heating power balance constraint are as follows:
In formula:Draw water power or generated output in period T are being simulated for hydroenergy storage station u,For wind turbine Available power of the group j in simulation period T,For available power of the photovoltaic unit k in simulation period T, PL,TFor power train The electric load united in simulation period T, HL,TFor thermic load of the electric system in simulation period T, NwFor Wind turbines number, NsFor Photovoltaic unit number,For heating power of the back pressure type cogeneration units m in simulation period T；
5.2.2) power of the assembling unit constraint
Power of the assembling unit constraint includes fired power generating unit power constraint, the constraint of back pressure type cogeneration units output power, steam extraction The constraint of formula cogeneration units output power, the constraint of steam-extracting type cogeneration units heating power and hydroenergy storage station power are about Beam；
Fired power generating unit power constraint are as follows:
The constraint of back pressure type cogeneration units output power are as follows:
The constraint of steam-extracting type cogeneration units output power are as follows:
The constraint of steam-extracting type cogeneration units heating power are as follows:
Hydroenergy storage station power constraint are as follows:
In formula:For the minimum power of fired power generating unit i,For the maximum power of fired power generating unit i,For back pressure The minimum output power of formula cogeneration units m,For the maximum output power of back pressure type cogeneration units m, For the minimum output power of steam-extracting type cogeneration units n,For the maximum output power of steam-extracting type cogeneration units n, HnminFor the minimum heating power of steam-extracting type cogeneration units n, HnmaxFor the maximum for hot merit of steam-extracting type cogeneration units n Rate,It draws water power for the minimum of hydroenergy storage station u,It draws water power for the maximum of hydroenergy storage station u, For the minimum generated output of hydroenergy storage station u,For the maximum power generation of hydroenergy storage station u；
5.2.3) fired power generating unit Climing constant
In formula: Rup,iFor the ratio of slope of climbing of fired power generating unit i, Rdown,iFor the downward climbing rate of fired power generating unit i, Δ t is every A simulation period step-length；
5.2.4) cogeneration units are for thermal confinement
Cogeneration units include back pressure type cogeneration units for thermal confinement and steam-extracting type cogeneration of heat and power machine for thermal confinement Group is for thermal confinement；
The heat supply of back pressure type cogeneration units m is constrained to linear equality constraints, specifically:
The heat supply of steam-extracting type cogeneration units n is constrained to linear inequality constraint, specifically:
In formula:WithRespectively back pressure type cogeneration units electricity/linear Relation Parameters, by back pressure type heat The operating condition figure that Electricity Federation produces unit m is determining,AndRespectively steam-extracting type cogeneration units are electric/linear Relation Parameters are determined by the operating condition figure of steam-extracting type cogeneration units n；
5.2.5) hydroenergy storage station constraint
Hydroenergy storage station constraint includes the constraint of upper storage reservoir storage capacity, the constraint of the hydroenergy storage station conservation of energy, water-storage The constraint of power station last current state and economy constraint；
The constraint of upper storage reservoir storage capacity are as follows:
In formula: Eu minFor the upper storage reservoir minimum energy storage of hydroenergy storage station u, Eu maxFor the upper storage reservoir of hydroenergy storage station u Maximum energy storage,For upper storage reservoir energy storage of hydroenergy storage station u at the end of simulating period T, calculation formula are as follows:
In formula: Eu,0For the initial energy storage of hydroenergy storage station u upper storage reservoir,It is hydroenergy storage station u in operation of drawing water Energy conversion efficiency under state,For energy conversion efficiency of the hydroenergy storage station u under generator operation state；
The constraint of the hydroenergy storage station conservation of energy are as follows:
The constraint of hydroenergy storage station last current state are as follows:
Economy constraint are as follows:
W=1,2 ..., NpAnd v ≠ w
5.2.6) system reserve constraint
In formula: R is electric system spare capacity；
6) it introduces Wind turbines abandonment power and photovoltaic unit abandons optical power, to the scene fire storage coordinated operation in step 5) Model is modified；
6.1) amendment step 5.1) in objective function, it is modified scene fire storage coordinated operation model objective function Are as follows:
In formula:For abandonment penalty coefficient,To abandon light penalty coefficient,It is Wind turbines j in simulation period T Interior abandonment power,For abandoning optical power of the photovoltaic unit k in simulation period T；
6.2) increase Wind turbines abandonment power constraint on the basis of step 5.2) constraint condition and photovoltaic unit abandons light Power constraint；
Wind turbines j abandonment power constraint are as follows:
Photovoltaic unit k abandons optical power constraint are as follows:
6.3) increase Wind turbines available power, Wind turbines practical consumption function on the basis of step 5.2) constraint condition Relation constraint and the practical consumption power of photovoltaic unit available power, photovoltaic unit and light between rate and Wind turbines abandonment power Lie prostrate the relation constraint that unit is abandoned between optical power；
Relationship between the practical consumption power of Wind turbines j available power, Wind turbines j and Wind turbines j abandonment power Constraint are as follows:
The practical consumption power of photovoltaic unit k available power, photovoltaic unit k and photovoltaic unit k abandon the relationship between optical power Constraint are as follows:
Wherein:To simulate the practical consumption power of Wind turbines j in period T,For photovoltaic unit k in simulation period T Practical consumption power；
6.4) amendment step 5.2.1) in the real-time electrical power Constraints of Equilibrium of electric system, modified electric system is electric in real time Power-balance constraint are as follows:
7) to the POWER SYSTEM STATE solution procedure 6 in step 4)) in it is modified scene fire storage coordinated operation model, obtain To fired power generating unit power, back pressure type cogeneration units output power and heating power, steam-extracting type cogeneration units for electric work Rate and heating power, pump-storage generator draw water power or the practical consumption power of generated output, Wind turbines and photovoltaic unit it is real Border dissolves power；
8) fired power generating unit power in the simulation period T of statistical simulation period D, back pressure type cogeneration units output power and Heating power, steam-extracting type cogeneration units output power and heating power, pump-storage generator draw water power or generated output, Wind turbines available power and practical consumption power, photovoltaic unit available power and practical consumption power；
9) number of segment T when if the simulation period T of simulation cycle D is equal to maximum analogmax, then enter step 10), otherwise, return Step 4) is returned, into next simulation period；
10) if simulation cycle D is equal to maximum analog periodicity Dmax, then it enters step 11), otherwise, return step 2), Into next simulation cycle；
11) export fired power generating unit timing power in each simulation cycle, back pressure type cogeneration units timing output power and Timing heating power, steam-extracting type cogeneration units timing output power and timing heating power, pump-storage generator timing are taken out Water power and timing generated output, Wind turbines available power and reality dissolve power, photovoltaic unit available power and actually disappear Receive power.
Compared with prior art, the beneficial effects of the present invention are:
The present invention considers northern China heat supply in winter phase electric system practical operation situation and hydroenergy storage station in electric power The peak load shifting effect played in system, meter and electric system realtime power balance, conventional thermal power unit Climing constant, back pressure type And the constraint conditions such as steam-extracting type cogeneration units coupled thermomechanics, hydroenergy storage station constraint, coordinated by establishing scene fire storage Moving model calculates electric system fired power generating unit power, back pressure type cogeneration units output power and heating power, steam extraction Draw water power or generated output, Wind turbines of formula cogeneration units output power and heating power, pump-storage generator are available Power and practical consumption power, photovoltaic unit available power and practical consumption power, can formulate for Operation of Electric Systems personnel Power of the assembling unit arrangement provides guidance, more comprehensive compared to method before.In addition, present invention introduces Wind turbines abandonment power and Photovoltaic unit abandons optical power, is used not only for electric system all consumption Wind turbines available power and photovoltaic unit available work Scene fire storage coordination optimization operation in the case of rate, additionally it is possible to for considering that the scene that abandonment is abandoned in light situation occurs in electric system Fire storage coordination optimization operation.
Detailed description of the invention
Fig. 1 is present invention scene fire storage coordination optimization operation method flow chart.
Specific embodiment
A specific embodiment of the invention is described in further details with reference to the accompanying drawing.
Referring to Fig. 1, a kind of scene fire storage coordination optimization operation method, including state sampling, status assessment and result statistics Three parts are specifically divided into following steps:
1. inputting initial data, including electric system electric load data, thermic load data, fired power generating unit (only powering) are forced Outage rate and mean repair time, back pressure type cogeneration units forced outage rate and mean repair time, steam-extracting type thermoelectricity connection Produce unit forced outage rate and mean repair time, historical wind speed data and history solar irradiance data；Define maximum analog Periodicity DmaxWith number of segment T when maximum analog in a simulation cyclemax, and simulation cycle D=1 is set, simulate period T=1；
2. obtaining fired power generating unit open state sequence in simulation cycle D, back pressure type heat based on state duration sampling method Electricity Federation produces unit open state sequence and steam-extracting type cogeneration units open state sequence.Assuming that conventional thermal power unit, back pressure The operating status duration and malfunction duration of formula cogeneration units and steam-extracting type cogeneration units obey Exponential distribution, specific steps are as follows:
2.1) fired power generating unit open state sequence is obtained based on state duration sampling method
2.1.1 it) according to fired power generating unit forced outage rate and mean repair time in step 1, calculates fired power generating unit and is averaged nothing Fail operation time, calculation formula are as follows:
In formula: MTTFthFor fired power generating unit mean time between failures, FORthFor fired power generating unit forced outage rate, MTTRthFor fired power generating unit mean repair time；
2.1.2) defined nucleotide sequenceIndicate fired power generating unit open state in simulation cycle D Sequence,Value illustrates that conventional thermal power unit is in operating status in the simulation period T of simulation cycle D when being 1,When value is 0 Illustrate that conventional thermal power unit is in stoppage in transit state in the simulation period T of simulation cycle D；Define conventional thermal motor in simulation cycle D Group has determined that open state numberAnd it sets
2.1.3 it) generates and obeys equally distributed random number in (0,1) sectionWith
2.1.4 the thermal power unit operation duration) is calculatedAnd trouble duration
Thermal power unit operation duration calculation formula are as follows:
Fired power generating unit trouble duration calculation formula are as follows:
2.1.5 it) is determined according to thermal power unit operation duration and trouble durationValue, formula are as follows:
In formula:To be rounded symbol downwards；
2.1.6 it) updatesIfWhen illustrating whole simulations of simulation cycle D The open state of fired power generating unit determines in section；Otherwise, return step 2.1.3), continue other simulations for determining simulation cycle D Fired power generating unit open state in period；
2.2) back pressure type cogeneration units open state sequence is obtained based on state duration sampling method
2.2.1) according to back pressure type cogeneration units forced outage rate and mean repair time in step 1, back pressure is calculated Formula cogeneration units mean time between failures, calculation formula are as follows:
In formula:For the back pressure type cogeneration units mean time between failures,For back pressure type Cogeneration units forced outage rate,For back pressure type cogeneration units mean repair time；
2.2.2) defined nucleotide sequenceIndicate back pressure type thermoelectricity connection in simulation cycle D Unit open state sequence is produced,Value illustrates in the simulation period T of simulation cycle D at back pressure type cogeneration units when being 1 In operating status,Value illustrates that back pressure type cogeneration units are in stoppage in transit shape in the simulation period T of simulation cycle D when being 0 State；It defines back pressure type cogeneration units in simulation cycle D and has determined that open state numberAnd it sets
2.2.3 it) generates and obeys equally distributed random number in (0,1) sectionWith
2.2.4 the back pressure type cogeneration units duration of operation) is calculatedAnd trouble duration
Back pressure type cogeneration units duration of operation calculation formula are as follows:
Back pressure type cogeneration units trouble duration calculation formula are as follows:
2.2.5 it) is determined according to the back pressure type cogeneration units duration of operation and trouble durationValue, Formula are as follows:
In formula:To be rounded symbol downwards.
2.2.6 it) updatesIfIllustrate the whole of simulation cycle D The open state for simulating back pressure type cogeneration units in the period determines；Otherwise, return step 2.2.3), continue to determine simulation Back pressure type cogeneration units open state in other simulation periods of period D；
2.3) steam-extracting type cogeneration units open state sequence is obtained based on state duration sampling method
2.3.1) according to steam-extracting type cogeneration units forced outage rate and mean repair time in step 1, steam extraction is calculated Formula cogeneration units mean time between failures, calculation formula are as follows:
In formula:For the steam-extracting type cogeneration units mean time between failures,For steam extraction Formula cogeneration units forced outage rate,For steam-extracting type cogeneration units mean repair time；
2.3.2) defined nucleotide sequenceIndicate steam-extracting type thermoelectricity connection in simulation cycle D Unit open state sequence is produced,Value illustrates in the simulation period T of simulation cycle D at steam-extracting type cogeneration units when being 1 In operating status,Value illustrates that steam-extracting type cogeneration units are in stoppage in transit shape in the simulation period T of simulation cycle D when being 0 State；It defines steam-extracting type cogeneration units in simulation cycle D and has determined that open state numberAnd it sets
2.3.3 it) generates and obeys equally distributed random number in (0,1) sectionWith
2.3.4 the steam-extracting type cogeneration units duration of operation) is calculatedAnd trouble duration
Steam-extracting type cogeneration units duration of operation calculation formula are as follows:
Steam-extracting type cogeneration units trouble duration calculation formula are as follows:
2.3.5 it) is determined according to the steam-extracting type cogeneration units duration of operation and trouble durationValue, Formula are as follows:
In formula:To be rounded symbol downwards.
2.3.6 it) updatesIfIllustrate the whole of simulation cycle D The open state for simulating steam-extracting type cogeneration units in the period determines；Otherwise, return step 2.3.3), continue to determine simulation Steam-extracting type cogeneration units open state in other simulation periods of period D.
3. obtaining Wind turbines available power sequences and photovoltaic unit available power sequence in simulation cycle D based on inverse transformation method Column；Assuming that simulation wind speed obeys two parameter Weibull distribution, the distribution function of two parameter Weibull distribution are as follows:
In formula: v is simulation wind speed, awFor the scale parameter of two parameter Weibull distribution, bwFor two parameter Weibull distribution Form parameter；
Assuming that simulated solar irradiation level obeys β distribution, the distribution function of β distribution are as follows:
In formula: r is simulated solar irradiation level, rmaxFor history solar irradiance maximum value, asFor β distribution location parameter, bsFor the form parameter of β distribution；
Wind turbines available power sequences and photovoltaic unit available power sequences in simulation cycle D are obtained based on inverse transformation method Specific steps are as follows:
3.1) Wind turbines available power sequences in simulation cycle D are obtained
3.1.1) according to the historical wind speed data in step 1, the form parameter and scale parameter of Weibull distribution are calculated；
The form parameter calculation formula of Weibull distribution are as follows:
In formula: swFor the standard deviation of historical wind speed data,For the average value of historical wind speed data；
The scale parameter calculation formula of Weibull distribution are as follows:
In formula: Γ () is Gama function；
3.1.2 it) generates and obeys equally distributed random number sequence in (0,1) section
3.1.3) defined nucleotide sequenceIndicate the simulation wind series in simulation cycle D, Utilize the random number sequence calculating simulation wind series generated in step 3.1.2), calculation formula are as follows:
3.1.4) defined nucleotide sequenceIndicate that the Wind turbines in simulation cycle D are available Power sequence；According to the relationship of simulation wind speed and Wind turbines available power, wind-powered electricity generation in the simulation period T of calculating simulation period D Unit available powerSecondary segmenting function representation, tool can be used by simulating the relationship between wind speed and Wind turbines available power Body calculation formula are as follows:
In formula: PWnFor Wind turbines rated power, vinWind speed is cut for Wind turbines；voutWind is cut out for Wind turbines Speed；vnFor Wind turbines rated wind speed；WithFor Wind turbines power curve parameter, by vinAnd vnIt determines；
3.2) photovoltaic unit available power sequences in simulation cycle D are obtained
3.2.1) according to the history solar irradiance data in step 1, the location parameter and form parameter of β distribution are calculated.β The location parameter calculation formula of distribution are as follows:
In formula:For the average value of history solar irradiance data, ssFor the variance of history solar irradiance data；
The form parameter calculation formula of β distribution are as follows:
3.2.2 it) generates and obeys β (as,bs) distribution random number sequence
3.2.3) defined nucleotide sequence sD={ sD,T, T=1,2 ..., TmaxIndicate simulation cycle D in simulated solar irradiation level Sequence utilizes the random number sequence calculating simulation solar irradiation degree series generated in step 3.2.2), calculation formula are as follows:
3.2.4) defined nucleotide sequenceIndicate that the photovoltaic unit in simulation cycle D is available Power sequence.According to the relationship of simulated solar irradiation level and photovoltaic unit available power, period T is simulated in calculating simulation period D Interior photovoltaic unit available powerRelationship between simulated solar irradiation level and photovoltaic unit available power can use secondary segmenting Function representation, specific formula for calculation are as follows:
In formula: PsnFor photovoltaic unit rated power, sinSolar irradiance is cut for photovoltaic unit；soutIt is cut for photovoltaic unit Sunny irradiation level；snFor the specified solar irradiance of photovoltaic unit；WithFor photovoltaic power of the assembling unit parameter of curve, by sinWith snIt determines.
4. obtaining the POWER SYSTEM STATE in the simulation period T of simulation cycle D, including electric load number, thermic load number, thermoelectricity Unit open state, back pressure type cogeneration units open state and steam-extracting type cogeneration units open state, Wind turbines Available power and photovoltaic unit available power.
5. the POWER SYSTEM STATE in pair step 4 solves scene fire storage coordinated operation model, thus to the mould of simulation cycle D Conventional thermal power unit power, back pressure type cogeneration units output power and heating power, steam-extracting type cogeneration of heat and power in quasi- period T Unit output power and heating power, pump-storage generator draw water power or generated output optimizes, and fortune is coordinated in scene fire storage Row model belongs to the mixed-integer nonlinear programming model in Optimized model, can be solved using branch and bound method.
If model there are feasible solution, illustrates that electric system can dissolve wind turbine completely in the simulation period T of simulation cycle D Group available power and photovoltaic unit available power, the i.e. practical consumption power of Wind turbines are equal to Wind turbines available power, photovoltaic The practical consumption power of unit is equal to photovoltaic unit available power, enters step 8；If feasible solution is not present in model, illustrate simulation week Electric system cannot dissolve Wind turbines available power and photovoltaic unit available power completely in the simulation period T of phase D, into step Rapid 6；Honourable fire storage coordinated operation model is made of objective function and constraint condition, concrete model are as follows:
5.1) objective function
The objective function of honourable fire storage coordinated operation model is Operation of Electric Systems cost minimization in entire simulation cycle, electricity Force system operating cost includes conventional thermal power unit operating cost, back pressure type cogeneration units operating cost, steam-extracting type thermoelectricity Coproduction unit operating cost and pump-storage generator operating cost, specifically:
Operating cost of the fired power generating unit i in simulation period T are as follows:
Operating cost of the back pressure type cogeneration units m in simulation period T is only related with its output power, specific to calculate Formula are as follows:
Operating cost of the steam-extracting type cogeneration units n in simulation period T is not only related with its output power, and also and its Heating power is related, specific formula for calculation are as follows:
Hydroenergy storage station operating cost is
In formula: NthFor fired power generating unit number,For back pressure type cogeneration units number,For steam-extracting type thermoelectricity connection Produce unit number, NpFor hydroenergy storage station number,The power in period T is being simulated for fired power generating unit i,For back pressure type heat Electricity Federation produces output power of the unit m in simulation period T,It is steam-extracting type cogeneration units n in simulation period T Output power,For heating power of the steam-extracting type cogeneration units n in simulation period T, ai,biAnd ciFor fired power generating unit Consumption characterisitic parameter,For back pressure type cogeneration units consumption characterisitic parameter,For steam-extracting type cogeneration units consumption characterisitic parameter,It is opened for hydroenergy storage station u The primary cost of machine,Primary cost is shut down for hydroenergy storage station u,It is hydroenergy storage station u in simulation period T The interior operating status indicator variable that draws water,Refer to for generator operation state of the hydroenergy storage station u in simulation period T+1 Show variable；
5.2) constraint condition
The constraint condition of honourable fire storage coordinated operation model include electric system realtime power Constraints of Equilibrium, the power of the assembling unit about Beam, fired power generating unit Climing constant, cogeneration units for thermal confinement, hydroenergy storage station constraint and electric system Reserve Constraint, It is specific as follows:
5.2.1) electric system realtime power Constraints of Equilibrium
For the stable operation for guaranteeing electric system, electric system realtime power Constraints of Equilibrium need to be considered, electric system is real-time Power-balance constraint includes real-time electrical power Constraints of Equilibrium and the constraint of real-time heating power balance；
Real-time electrical power Constraints of Equilibrium are as follows:
Real-time heating power balance constraint are as follows:
In formula:Draw water power or generated output in period T are being simulated for hydroenergy storage station u,For wind turbine Available power of the group j in simulation period T,For available power of the photovoltaic unit k in simulation period T, PL,TFor power train The electric load united in simulation period T, HL,TFor thermic load of the electric system in simulation period T, NwFor Wind turbines number, NsFor Photovoltaic unit number,For heating power of the back pressure type cogeneration units m in simulation period T；
5.2.2) power of the assembling unit constraint
Power of the assembling unit constraint includes fired power generating unit power constraint, the constraint of back pressure type cogeneration units output power, steam extraction The constraint of formula cogeneration units output power, the constraint of steam-extracting type cogeneration units heating power and hydroenergy storage station power are about Beam；
Fired power generating unit power constraint are as follows:
The constraint of back pressure type cogeneration units output power are as follows:
The constraint of steam-extracting type cogeneration units output power are as follows:
The constraint of steam-extracting type cogeneration units heating power are as follows:
Hydroenergy storage station power constraint are as follows:
In formula:For the minimum power of fired power generating unit i,For the maximum power of fired power generating unit i,For back pressure The minimum output power of formula cogeneration units m,For the maximum output power of back pressure type cogeneration units m, For the minimum output power of steam-extracting type cogeneration units n,For the maximum output power of steam-extracting type cogeneration units n, HnminFor the minimum heating power of steam-extracting type cogeneration units n, HnmaxFor the maximum for hot merit of steam-extracting type cogeneration units n Rate,It draws water power for the minimum of hydroenergy storage station u,It draws water power for the maximum of hydroenergy storage station u, For the minimum generated output of hydroenergy storage station u,For the maximum power generation of hydroenergy storage station u；
5.2.3) fired power generating unit Climing constant
In formula: Rup,iFor the ratio of slope of climbing of fired power generating unit i, Rdown,iFor the downward climbing rate of fired power generating unit i, Δ t is every A simulation period step-length；
5.2.4) cogeneration units are for thermal confinement
Cogeneration units include back pressure type cogeneration units and steam-extracting type cogeneration units, back pressure type cogeneration of heat and power The heat supply of unit m is constrained to linear equality constraints, specifically:
The heat supply of steam-extracting type cogeneration units n is constrained to linear inequality constraint, specifically:
In formula:WithRespectively back pressure type cogeneration units electricity/linear Relation Parameters, by back pressure type heat The operating condition figure that Electricity Federation produces unit m is determining,AndRespectively steam-extracting type cogeneration units are electric/linear Relation Parameters are determined by the operating condition figure of steam-extracting type cogeneration units n；
5.2.5) hydroenergy storage station constraint
Hydroenergy storage station constraint includes the constraint of upper storage reservoir storage capacity, the constraint of the hydroenergy storage station conservation of energy, water-storage The constraint of power station last current state and economy constraint；
The constraint of upper storage reservoir storage capacity are as follows:
In formula: EuminFor the upper storage reservoir minimum energy storage of hydroenergy storage station u, EumaxFor hydroenergy storage station u upper storage reservoir most Big energy storage,For upper storage reservoir energy storage of hydroenergy storage station u at the end of simulating period T, calculation formula are as follows:
In formula: Eu,0For the initial energy storage of hydroenergy storage station u upper storage reservoir,It is hydroenergy storage station u in operation of drawing water Energy conversion efficiency under state,For energy conversion efficiency of the hydroenergy storage station u under generator operation state；
The constraint of the hydroenergy storage station conservation of energy are as follows:
The constraint of hydroenergy storage station last current state are as follows:
Economy constraint are as follows:
V=1,2 ..., Np, w=1,2 ..., NpAnd v ≠ w
5.2.6) system reserve constraint
In formula: R is electric system spare capacity.
6. introducing Wind turbines abandonment power and photovoltaic unit abandoning optical power, to the scene fire storage coordinated operation in step 5 Model is modified；
6.1) amendment step 5.1) in objective function.
To guarantee Wind turbines available power and the consumption as much as possible of photovoltaic unit available power, in Operation of Electric Systems At addition abandonment cost in this item and light cost is abandoned, modified scene fire stores up the objective function of coordinated operation model are as follows:
In formula:For abandonment penalty coefficient, size is determined by the urgent degree of consumption Wind turbines available power, To abandon light penalty coefficient, size is determined by the urgent degree of consumption photovoltaic unit available power,Exist for Wind turbines j The abandonment power in period T is simulated,For abandoning optical power of the photovoltaic unit k in simulation period T；
6.2) increase Wind turbines abandonment power constraint on the basis of step 5.2) constraint condition and photovoltaic unit abandons light Power constraint.Wind turbines j abandonment power constraint are as follows:
Photovoltaic unit k abandons optical power constraint are as follows:
6.3) increase Wind turbines available power, Wind turbines practical consumption function on the basis of step 5.2) constraint condition Relation constraint and the practical consumption power of photovoltaic unit available power, photovoltaic unit and light between rate and Wind turbines abandonment power Lie prostrate the relation constraint that unit is abandoned between optical power；
Relationship between the practical consumption power of Wind turbines j available power, Wind turbines j and Wind turbines j abandonment power Constraint are as follows:
The practical consumption power of photovoltaic unit k available power, photovoltaic unit k and photovoltaic unit k abandon the relationship between optical power Constraint are as follows:
Wherein:To simulate the practical consumption power of Wind turbines j in period T,For photovoltaic unit k in simulation period T Practical consumption power；
6.4) amendment step 5.2.1) in the real-time electrical power Constraints of Equilibrium of electric system
The modified real-time electrical power Constraints of Equilibrium of electric system are as follows:
7. modified scene fire storage coordinated operation model in the POWER SYSTEM STATE solution procedure 6 in pair step 4, thus To fired power generating unit power, back pressure type cogeneration units output power and heating power, pumping in the simulation period T of simulation cycle D Draw water power or generated output, Wind turbines of vapour formula cogeneration units output power and heating power, pump-storage generator are real Border consumption power and the practical consumption power of photovoltaic unit optimize, and modified scene fire storage coordinated operation model belongs to optimization Mixed-integer nonlinear programming model in model can be solved using branch and bound method.
8. fired power generating unit power in the simulation period T of statistical simulation period D, back pressure type cogeneration units output power and Heating power, steam-extracting type cogeneration units output power and heating power, pump-storage generator draw water power or generated output, The practical consumption power of Wind turbines and the practical consumption power of photovoltaic unit.
9. number of segment T when if the simulation period T of simulation cycle D is equal to maximum analogmax, then 10 are entered step, otherwise, is returned Step 4 is returned, into next simulation period.
10. if simulation cycle D is equal to maximum analog periodicity Dmax, then 11 are entered step, otherwise, return step 2, into Enter next simulation cycle.
11. export fired power generating unit power, back pressure type cogeneration units output power and heating power in each simulation cycle, Steam-extracting type cogeneration units output power and heating power, pump-storage generator draw water power and generated output, Wind turbines Practical consumption power, the practical consumption power of photovoltaic unit.

Claims (3)

1. a kind of scene fire storage coordination optimization operation method, which comprises the following steps:
1) initial data is inputted, including electric system electric load data, thermic load data, fired power generating unit forced outage rate and average Repair time, back pressure type cogeneration units forced outage rate and mean repair time, steam-extracting type cogeneration units, which are forced, stops Fortune rate and mean repair time, historical wind speed data and history solar irradiance data；Define maximum analog periodicity DmaxWith one Number of segment T when maximum analog in a simulation cyclemax, and simulation cycle D=1 is set, simulate period T=1；
2) fired power generating unit open state sequence in simulation cycle D, back pressure type thermoelectricity connection are obtained based on state duration sampling method Produce unit open state sequence and steam-extracting type cogeneration units open state sequence；
3) Wind turbines available power sequences and photovoltaic unit available power sequences in simulation cycle D are obtained based on inverse transformation method；
4) POWER SYSTEM STATE in the simulation period T of simulation cycle D, including electric load data, thermic load data, thermoelectricity are obtained Unit open state, back pressure type cogeneration units open state and steam-extracting type cogeneration units open state, Wind turbines Available power and photovoltaic unit available power；
5) scene fire storage coordinated operation model is solved to the POWER SYSTEM STATE in step 4) to say if there are feasible solutions for model Electric system can dissolve Wind turbines available power and photovoltaic unit available work completely in the simulation period T of bright simulation cycle D 8) rate enters step；If feasible solution is not present in model, illustrate that electric system cannot disappear completely in the simulation period T of simulation cycle D It receives Wind turbines available power and photovoltaic unit available power, enters step 6)；Honourable fire storage coordinated operation model is by target letter Several and constraint condition composition, concrete model are as follows:
5.1) objective function
The objective function of honourable fire storage coordinated operation model is Operation of Electric Systems cost minimization in simulation cycle D, electric system Operating cost includes thermal power unit operation cost, back pressure type cogeneration units operating cost, steam-extracting type cogeneration units fortune Row cost and pump-storage generator operating cost, specifically:
Operating cost of the fired power generating unit i in simulation period T are as follows:
Operating cost of the back pressure type cogeneration units m in simulation period T are as follows:
Operating cost of the steam-extracting type cogeneration units n in simulation period T are as follows:
Hydroenergy storage station operating cost are as follows:
In formula: NthFor fired power generating unit number,For back pressure type cogeneration units number,For steam-extracting type cogeneration units Number, NpFor hydroenergy storage station number,The power in period T is being simulated for fired power generating unit i,For back pressure type cogeneration of heat and power Output power of the unit m in simulation period T,Be steam-extracting type cogeneration units n simulation period T in for electric work Rate,For heating power of the steam-extracting type cogeneration units n in simulation period T, ai,biAnd ciFor fired power generating unit consumption Characterisitic parameter,For back pressure type cogeneration units consumption characterisitic parameter,For Steam-extracting type cogeneration units consumption characterisitic parameter,The primary cost of machine is opened for hydroenergy storage station u,It is stored to draw water Energy power station u shuts down primary cost,The operating status indicator variable that draws water in period T is being simulated for hydroenergy storage station u,The generator operation state instruction variable in period T+1 is being simulated for hydroenergy storage station u, when hydroenergy storage station u is in mould In quasi- period T in draw water operating status when,Value is 1,Value is 0；When pump-storage generator power station u is in simulation period T It is interior be in generator operation state when,Value is 1,Value is 0；
5.2) constraint condition
The constraint condition of honourable fire storage coordinated operation model include electric system realtime power Constraints of Equilibrium, power of the assembling unit constraint, Fired power generating unit Climing constant, cogeneration units are for thermal confinement, hydroenergy storage station constraint and electric system Reserve Constraint, specifically Are as follows:
5.2.1) electric system realtime power Constraints of Equilibrium
Electric system realtime power Constraints of Equilibrium includes real-time electrical power Constraints of Equilibrium and the constraint of real-time heating power balance；
Real-time electrical power Constraints of Equilibrium are as follows:
Real-time heating power balance constraint are as follows:
In formula:Draw water power or generated output in period T are being simulated for hydroenergy storage station u,For Wind turbines j Available power in simulation period T,For available power of the photovoltaic unit k in simulation period T, PL,TExist for electric system Simulate the electric load in period T, HL,TFor thermic load of the electric system in simulation period T, NwFor Wind turbines number, NsFor photovoltaic Unit number,For heating power of the back pressure type cogeneration units m in simulation period T；
5.2.2) power of the assembling unit constraint
Power of the assembling unit constraint includes fired power generating unit power constraint, the constraint of back pressure type cogeneration units output power, steam-extracting type heat Electricity Federation produces the constraint of unit output power, the constraint of steam-extracting type cogeneration units heating power and hydroenergy storage station power constraint；
Fired power generating unit power constraint are as follows:
The constraint of back pressure type cogeneration units output power are as follows:
The constraint of steam-extracting type cogeneration units output power are as follows:
The constraint of steam-extracting type cogeneration units heating power are as follows:
Hydroenergy storage station power constraint are as follows:
In formula:For the minimum power of fired power generating unit i,For the maximum power of fired power generating unit i,For back pressure type heat Electricity Federation produces the minimum output power of unit m,For the maximum output power of back pressure type cogeneration units m,To take out The minimum output power of vapour formula cogeneration units n,For the maximum output power of steam-extracting type cogeneration units n, Hnmin For the minimum heating power of steam-extracting type cogeneration units n, HnmaxFor the maximum heating power of steam-extracting type cogeneration units n,It draws water power for the minimum of hydroenergy storage station u,It draws water power for the maximum of hydroenergy storage station u,For The minimum generated output of hydroenergy storage station u,For the maximum power generation of hydroenergy storage station u；
5.2.3) fired power generating unit Climing constant
In formula: Rup,iFor the ratio of slope of climbing of fired power generating unit i, Rdown,iFor the downward climbing rate of fired power generating unit i, Δ t is each mould Quasi- period step-length；
5.2.4) cogeneration units are for thermal confinement
Cogeneration units include back pressure type cogeneration units for thermal confinement and the confession of steam-extracting type cogeneration units for thermal confinement Thermal confinement；
The heat supply of back pressure type cogeneration units m is constrained to linear equality constraints, specifically:
The heat supply of steam-extracting type cogeneration units n is constrained to linear inequality constraint, specifically:
In formula:WithRespectively back pressure type cogeneration units electricity/linear Relation Parameters, by back pressure type cogeneration of heat and power The operating condition figure of unit m is determining,AndRespectively steam-extracting type cogeneration units electricity/hot line sexual intercourse ginseng Number, is determined by the operating condition figure of steam-extracting type cogeneration units n；
5.2.5) hydroenergy storage station constraint
Hydroenergy storage station constraint includes the constraint of upper storage reservoir storage capacity, the constraint of the hydroenergy storage station conservation of energy, hydroenergy storage station Last current state constraint and economy constraint；
The constraint of upper storage reservoir storage capacity are as follows:
In formula: EuminFor the upper storage reservoir minimum energy storage of hydroenergy storage station u, EumaxIt is stored up for the upper storage reservoir maximum of hydroenergy storage station u Can,For upper storage reservoir energy storage of hydroenergy storage station u at the end of simulating period T, calculation formula are as follows:
In formula: Eu,0For the initial energy storage of hydroenergy storage station u upper storage reservoir,It is hydroenergy storage station u under operating status of drawing water Energy conversion efficiency,For energy conversion efficiency of the hydroenergy storage station u under generator operation state；
The constraint of the hydroenergy storage station conservation of energy are as follows:
The constraint of hydroenergy storage station last current state are as follows:
Economy constraint are as follows:
And v ≠ w
5.2.6) system reserve constraint
In formula: R is electric system spare capacity；
6) it introduces Wind turbines abandonment power and photovoltaic unit abandons optical power, to the scene fire storage coordinated operation model in step 5) It modifies；
6.1) amendment step 5.1) in objective function, it is modified scene fire storage coordinated operation model objective function are as follows:
In formula:For abandonment penalty coefficient,To abandon light penalty coefficient,For abandoning of the Wind turbines j in simulation period T Wind power,For abandoning optical power of the photovoltaic unit k in simulation period T；
6.2) increase Wind turbines abandonment power constraint on the basis of step 5.2) constraint condition and photovoltaic unit abandons optical power Constraint；
Wind turbines j abandonment power constraint are as follows:
Photovoltaic unit k abandons optical power constraint are as follows:
6.3) on the basis of step 5.2) constraint condition increase Wind turbines available power, the practical consumption power of Wind turbines and Relation constraint and the practical consumption power of photovoltaic unit available power, photovoltaic unit and photovoltaic machine between Wind turbines abandonment power Group abandons the relation constraint between optical power；
Relation constraint between the practical consumption power of Wind turbines j available power, Wind turbines j and Wind turbines j abandonment power Are as follows:
The practical consumption power of photovoltaic unit k available power, photovoltaic unit k and photovoltaic unit k abandon the relation constraint between optical power Are as follows:
Wherein:To simulate the practical consumption power of Wind turbines j in period T,It is practical for photovoltaic unit k in simulation period T Dissolve power；
6.4) amendment step 5.2.1) in the real-time electrical power Constraints of Equilibrium of electric system, the modified real-time electrical power of electric system Constraints of Equilibrium are as follows:
7) to the POWER SYSTEM STATE solution procedure 6 in step 4)) in it is modified scene fire storage coordinated operation model, obtain fire The electric power of the assembling unit, back pressure type cogeneration units output power and heating power, steam-extracting type cogeneration units output power and Heating power, pump-storage generator draw water power or the practical consumption power of generated output, Wind turbines and photovoltaic unit is practical disappears Receive power；
8) fired power generating unit power, back pressure type cogeneration units output power and heat supply in the simulation period T of statistical simulation period D Power, steam-extracting type cogeneration units output power and heating power, pump-storage generator draw water power or generated output, wind-powered electricity generation Unit available power and practical consumption power, photovoltaic unit available power and practical consumption power；
9) number of segment T when if the simulation period T of simulation cycle D is equal to maximum analogmax, then enter step 11), otherwise, return to step It is rapid 4), into next simulation period；
10) if simulation cycle D is equal to maximum analog periodicity Dmax, then it enters step 11), otherwise, return step 2), enter Next simulation cycle；
11) fired power generating unit timing power, back pressure type cogeneration units timing output power and timing in each simulation cycle are exported Heating power, steam-extracting type cogeneration units timing output power and timing heating power, pump-storage generator timing are drawn water function Rate and timing generated output, Wind turbines available power and practical consumption power, photovoltaic unit available power and practical consumption function Rate.
2. scene fire storage coordination optimization operation method according to claim 1, it is characterised in that: be based on shape in the step 2) State duration sampling method obtains fired power generating unit open state sequence in simulation cycle D, back pressure type cogeneration units booting shape State sequence and steam-extracting type cogeneration units open state sequence；Assuming that fired power generating unit, back pressure type cogeneration units and steam extraction The operating status duration and malfunction duration of formula cogeneration units obey exponential distribution, specific steps are as follows:
2.1) fired power generating unit open state sequence is obtained based on state duration sampling method
2.1.1) according to fired power generating unit forced outage rate and mean repair time in step 1), fired power generating unit Mean is calculated Working time, calculation formula are as follows:
In formula: MTTFthFor fired power generating unit mean time between failures, FORthFor fired power generating unit forced outage rate, MTTRthFor Fired power generating unit mean repair time；
2.1.2) defined nucleotide sequenceIndicate fired power generating unit open state sequence in simulation cycle D Column,Value illustrates that fired power generating unit is in operating status in the simulation period T of simulation cycle D when being 1,Value illustrates mould when being 0 Fired power generating unit is in stoppage in transit state in the simulation period T of quasi-periodicity D；It defines fired power generating unit in simulation cycle D and has determined that booting shape State numberAnd it sets
2.1.3 it) generates and obeys equally distributed random number in (0,1) sectionWith
2.1.4 the thermal power unit operation duration) is calculatedAnd trouble durationThe thermal power unit operation duration Calculation formula are as follows:
The calculation formula of fired power generating unit trouble duration are as follows:
2.1.5 it) according to the fired power generating unit duration of operation and trouble duration in step 2.1.4), determinesValue, formula Are as follows:
In formula:To be rounded symbol downwards；
2.1.6 it) updatesIfIllustrate in whole simulation periods of simulation cycle D The open state of fired power generating unit determines；Otherwise, return step 2.1.3), continue other simulation periods for determining simulation cycle D Interior fired power generating unit open state；
2.2) back pressure type cogeneration units open state sequence is obtained based on state duration sampling method
2.2.1) according to back pressure type cogeneration units forced outage rate and mean repair time in step 1), back pressure type heat is calculated Electricity Federation produces unit mean time between failures, calculation formula are as follows:
In formula:For the back pressure type cogeneration units mean time between failures,For back pressure type thermoelectricity Coproduction unit forced outage rate,For back pressure type cogeneration units mean repair time；
2.2.2) defined nucleotide sequenceIndicate back pressure type cogeneration of heat and power machine in simulation cycle D Group open state sequence,Value illustrates that back pressure type cogeneration units are in fortune in the simulation period T of simulation cycle D when being 1 Row state,Value illustrates that back pressure type cogeneration units are in stoppage in transit state in the simulation period T of simulation cycle D when being 0； It defines back pressure type cogeneration units in simulation cycle D and has determined that open state numberAnd it sets
2.2.3 it) generates and obeys equally distributed random number in (0,1) sectionWith
2.2.4 the back pressure type cogeneration units duration of operation) is calculatedAnd trouble duration
Back pressure type cogeneration units duration of operation calculation formula are as follows:
Back pressure type cogeneration units trouble duration calculation formula are as follows:
2.2.5 it) according to the back pressure type cogeneration units duration of operation and trouble duration in step 2.2.4), determinesValue, formula are as follows:
In formula:To be rounded symbol downwards；
2.2.6 it) updatesIfWhen illustrating whole simulations of simulation cycle D The open state of back pressure type cogeneration units determines in section；Otherwise, return step 2.2.3), continue to determine simulation cycle D Other simulation periods in back pressure type cogeneration units open state；
2.3) steam-extracting type cogeneration units open state sequence is obtained based on state duration sampling method
2.3.1) according to steam-extracting type cogeneration units forced outage rate and mean repair time in step 1), steam-extracting type heat is calculated Electricity Federation produces unit mean time between failures, calculation formula are as follows:
In formula:For the steam-extracting type cogeneration units mean time between failures,For steam-extracting type heat Electricity Federation produces unit forced outage rate,For steam-extracting type cogeneration units mean repair time；
2.3.2) defined nucleotide sequenceIndicate steam-extracting type cogeneration of heat and power machine in simulation cycle D Group open state sequence,Value illustrates that steam-extracting type cogeneration units are in fortune in the simulation period T of simulation cycle D when being 1 Row state,Value illustrates that steam-extracting type cogeneration units are in stoppage in transit state in the simulation period T of simulation cycle D when being 0； It defines steam-extracting type cogeneration units in simulation cycle D and has determined that open state numberAnd it sets
2.3.3 it) generates and obeys equally distributed random number in (0,1) sectionWith
2.3.4) calculate the steam-extracting type cogeneration units duration of operationAnd trouble duration
Steam-extracting type cogeneration units duration of operation calculation formula are as follows:
Steam-extracting type cogeneration units trouble duration calculation formula are as follows:
2.3.5 it) according to the steam-extracting type cogeneration units duration of operation and trouble duration in step 2.3.4), determinesValue, formula are as follows:
In formula:To be rounded symbol downwards；
2.3.6 it) updatesIfIllustrate whole simulations of simulation cycle D The open state of steam-extracting type cogeneration units determines in period；Otherwise, return step 2.3.3), continue to determine simulation cycle Steam-extracting type cogeneration units open state in other simulation periods of D.
3. according to claim 1 scene fire storage coordination optimization operation method, it is characterised in that: in the step 3) by pair The probability-distribution function of simulation wind speed and simulated solar irradiation level carries out sampling and obtains Wind turbines available power and photovoltaic unit Available power；
Assuming that simulation wind speed obeys two parameter Weibull distribution, the distribution function of two parameter Weibull distribution are as follows:
In formula: v is simulation wind speed, awFor the scale parameter of two parameter Weibull distribution, bwFor the shape of two parameter Weibull distribution Parameter；
Assuming that simulated solar irradiation level obeys β distribution, the distribution function of β distribution are as follows:
In formula: r is simulated solar irradiation level, rmaxFor history solar irradiance maximum value, asFor the location parameter of β distribution, bsFor β The form parameter of distribution；
3.1) Wind turbines available power in simulation cycle D is obtained
3.1.1) according to the historical wind speed data in step 1), the form parameter and scale parameter of Weibull distribution are calculated；
The form parameter calculation formula of Weibull distribution are as follows:
In formula: swFor the standard deviation of historical wind speed data,For the average value of historical wind speed data；
The scale parameter calculation formula of Weibull distribution are as follows:
In formula: Γ () is Gama function；
3.1.2 it) generates and obeys equally distributed random number sequence in (0,1) section
3.1.3) defined nucleotide sequenceIt indicates the simulation wind series in simulation cycle D, utilizes Step 3.1.2) in generate random number sequence calculating simulation wind series, calculation formula are as follows:
3.1.4) defined nucleotide sequenceIndicate the Wind turbines available power in simulation cycle D Sequence；According to the relationship of simulation wind speed and Wind turbines available power, Wind turbines in the simulation period T of calculating simulation period D Available powerSpecific formula for calculation are as follows:
In formula: PWnFor Wind turbines rated power, vinWind speed, v are cut for Wind turbinesoutFor Wind turbines cut-out wind speed, vnFor Wind turbines rated wind speed,WithFor Wind turbines power curve parameter, by vinAnd vnIt determines；
3.2.1) obtain simulation cycle D in photovoltaic unit available power sequences according to the history solar irradiation degree in step 1) According to the location parameter and form parameter of calculating β distribution；
The location parameter calculation formula of β distribution are as follows:
In formula:For the average value of history solar irradiance data, ssFor the variance of history solar irradiance data；
The form parameter calculation formula of β distribution are as follows:
3.2.2 it) generates and obeys β (as,bs) distribution random number sequence
3.2.3) defined nucleotide sequence sD={ sD,T, T=1,2 ..., TmaxIndicate that the simulated solar in simulation cycle D irradiates degree series, Utilize the random number sequence calculating simulation solar irradiation degree series generated in step 3.2.2), calculation formula are as follows:
3.2.4) defined nucleotide sequenceIndicate the photovoltaic unit available power sequence in simulation cycle D Column；According to the relationship of simulated solar irradiation level and photovoltaic unit available power, photovoltaic in the simulation period T of calculating simulation period D Unit available powerSpecific formula for calculation are as follows:
In formula: PsnFor photovoltaic unit rated power, sinSolar irradiance is cut for photovoltaic unit；soutIt is cut out too for photovoltaic unit Positive irradiation level；snFor the specified solar irradiance of photovoltaic unit；WithFor photovoltaic power of the assembling unit parameter of curve, by sinAnd snReally It is fixed.
CN201811376350.2A 2018-11-19 2018-11-19 Wind, light, fire and storage coordination optimization operation method Active CN109245180B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201811376350.2A CN109245180B (en) 2018-11-19 2018-11-19 Wind, light, fire and storage coordination optimization operation method

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201811376350.2A CN109245180B (en) 2018-11-19 2018-11-19 Wind, light, fire and storage coordination optimization operation method

Publications (2)

Publication Number Publication Date
CN109245180A true CN109245180A (en) 2019-01-18
CN109245180B CN109245180B (en) 2020-08-04

Family

ID=65075933

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201811376350.2A Active CN109245180B (en) 2018-11-19 2018-11-19 Wind, light, fire and storage coordination optimization operation method

Country Status (1)

Country Link
CN (1) CN109245180B (en)

Cited By (11)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109724269A (en) * 2019-02-20 2019-05-07 清华四川能源互联网研究院 The full spectrum cogeneration system of solar energy and energy storage configuration method
CN110474352A (en) * 2019-07-05 2019-11-19 国家电网有限公司 Operating condition evaluation method, device and its electronic equipment of generator motor
CN112072709A (en) * 2020-07-29 2020-12-11 国网甘肃省电力公司电力科学研究院 Coordinated optimization method for joint peak regulation system under given outgoing power curve
CN113098071A (en) * 2021-03-17 2021-07-09 华中科技大学 Method and device for establishing operation model of wind power photovoltaic system
CN114142536A (en) * 2021-12-03 2022-03-04 国家电网有限公司西北分部 Multi-type unit coordination method considering capacity reserve
CN114285088A (en) * 2021-11-29 2022-04-05 中国华能集团清洁能源技术研究院有限公司 Energy system optimal configuration method and system considering thermal power generating unit constraints
CN114285089A (en) * 2021-11-29 2022-04-05 中国华能集团清洁能源技术研究院有限公司 Method and system for optimizing starting and stopping of thermal power generating unit in wind-solar-energy fire storage system
CN114336778A (en) * 2021-11-29 2022-04-12 中国华能集团清洁能源技术研究院有限公司 Method and device for determining starting sequence of thermoelectric generator set in wind-light-fire storage system
CN114336776A (en) * 2021-11-29 2022-04-12 中国华能集团清洁能源技术研究院有限公司 Thermal power generating unit starting sequence determination method and system based on optimized binary tree
CN114336777A (en) * 2021-11-29 2022-04-12 中国华能集团清洁能源技术研究院有限公司 Thermal power generating unit starting sequence determination method and system considering energy utilization sequence
CN114938011A (en) * 2022-07-07 2022-08-23 中国长江三峡集团有限公司 Wind-solar-fire-storage system combined operation method considering energy storage optimization configuration

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2015036950A1 (en) * 2013-09-12 2015-03-19 Frezite - Equipamentos Energéticos & Ambiente, Lda. Power distribution hub
US8996187B2 (en) * 2011-07-18 2015-03-31 Nec Laboratories America, Inc. Optimal storage sizing for integrating wind and load forecast uncertainties
CN104795846A (en) * 2015-04-24 2015-07-22 清华大学 Optimized operation method of pumped-storage power station and wind power combined system
CN107171367A (en) * 2017-05-19 2017-09-15 山东大学 Fired power generating unit decision making of combinatorial optimization method of the wind-powered electricity generation with energy-storage system under complementary
CN107332286A (en) * 2017-07-27 2017-11-07 清华大学 A kind of cogeneration of heat and power containing heat accumulation and wind-powered electricity generation coordinated scheduling method

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US8996187B2 (en) * 2011-07-18 2015-03-31 Nec Laboratories America, Inc. Optimal storage sizing for integrating wind and load forecast uncertainties
WO2015036950A1 (en) * 2013-09-12 2015-03-19 Frezite - Equipamentos Energéticos & Ambiente, Lda. Power distribution hub
CN104795846A (en) * 2015-04-24 2015-07-22 清华大学 Optimized operation method of pumped-storage power station and wind power combined system
CN107171367A (en) * 2017-05-19 2017-09-15 山东大学 Fired power generating unit decision making of combinatorial optimization method of the wind-powered electricity generation with energy-storage system under complementary
CN107332286A (en) * 2017-07-27 2017-11-07 清华大学 A kind of cogeneration of heat and power containing heat accumulation and wind-powered electricity generation coordinated scheduling method

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
CAN WANG 等: "Evaluation ofPower Grids" Renewable Energy", 《2016 IEEE PES ASIA-PACIFIC POWER AND ENERGY CONFERENCE - XI" AN - CHINA》 *
黄亭 等: "风火联合系统的多目标随机优化调度方法的研究", 《西安理工大学学报》 *

Cited By (19)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109724269A (en) * 2019-02-20 2019-05-07 清华四川能源互联网研究院 The full spectrum cogeneration system of solar energy and energy storage configuration method
CN109724269B (en) * 2019-02-20 2024-03-22 清华四川能源互联网研究院 Solar full-spectrum cogeneration system and energy storage configuration method
CN110474352A (en) * 2019-07-05 2019-11-19 国家电网有限公司 Operating condition evaluation method, device and its electronic equipment of generator motor
CN112072709A (en) * 2020-07-29 2020-12-11 国网甘肃省电力公司电力科学研究院 Coordinated optimization method for joint peak regulation system under given outgoing power curve
CN113098071A (en) * 2021-03-17 2021-07-09 华中科技大学 Method and device for establishing operation model of wind power photovoltaic system
CN114336777B (en) * 2021-11-29 2023-08-25 中国华能集团清洁能源技术研究院有限公司 Thermal power generating unit startup sequence determination method and system considering energy utilization sequence
CN114285088A (en) * 2021-11-29 2022-04-05 中国华能集团清洁能源技术研究院有限公司 Energy system optimal configuration method and system considering thermal power generating unit constraints
CN114336778A (en) * 2021-11-29 2022-04-12 中国华能集团清洁能源技术研究院有限公司 Method and device for determining starting sequence of thermoelectric generator set in wind-light-fire storage system
CN114336776A (en) * 2021-11-29 2022-04-12 中国华能集团清洁能源技术研究院有限公司 Thermal power generating unit starting sequence determination method and system based on optimized binary tree
CN114336777A (en) * 2021-11-29 2022-04-12 中国华能集团清洁能源技术研究院有限公司 Thermal power generating unit starting sequence determination method and system considering energy utilization sequence
CN114285088B (en) * 2021-11-29 2023-10-10 中国华能集团清洁能源技术研究院有限公司 Energy system optimal configuration method and system considering thermal power unit constraint
CN114336778B (en) * 2021-11-29 2023-09-22 中国华能集团清洁能源技术研究院有限公司 Method and device for determining starting sequence of thermal power generating unit in wind, light and fire storage system
CN114285089A (en) * 2021-11-29 2022-04-05 中国华能集团清洁能源技术研究院有限公司 Method and system for optimizing starting and stopping of thermal power generating unit in wind-solar-energy fire storage system
CN114336776B (en) * 2021-11-29 2023-08-25 中国华能集团清洁能源技术研究院有限公司 Thermal power generating unit startup sequence determining method and system based on optimized binary tree
CN114285089B (en) * 2021-11-29 2023-08-29 中国华能集团清洁能源技术研究院有限公司 Method and system for optimizing start-stop of thermal power generating unit in wind, light and fire storage system
CN114142536B (en) * 2021-12-03 2023-09-26 国家电网有限公司西北分部 Multi-type unit coordination method considering capacity reserve
CN114142536A (en) * 2021-12-03 2022-03-04 国家电网有限公司西北分部 Multi-type unit coordination method considering capacity reserve
CN114938011B (en) * 2022-07-07 2022-10-14 中国长江三峡集团有限公司 Wind-solar-fire-storage system combined operation method considering energy storage optimization configuration
CN114938011A (en) * 2022-07-07 2022-08-23 中国长江三峡集团有限公司 Wind-solar-fire-storage system combined operation method considering energy storage optimization configuration

Also Published As

Publication number Publication date
CN109245180B (en) 2020-08-04

Similar Documents

Publication Publication Date Title
CN109245180A (en) A kind of scene fire storage coordination optimization operation method
Yang et al. Optimal capacity and operation strategy of a solar-wind hybrid renewable energy system
CN103580063B (en) A kind of method of large-scale grid connection wind-powered electricity generation of dissolving based on demanding party's response
CN106099993B (en) A kind of power source planning method for adapting to new energy and accessing on a large scale
CN109103929B (en) Power distribution network economic optimization scheduling method based on improved dynamic kriging model
CN103762589A (en) Method for optimizing new energy capacity ratio in layers in power grid
CN107332286A (en) A kind of cogeneration of heat and power containing heat accumulation and wind-powered electricity generation coordinated scheduling method
CN111064229A (en) Wind-light-gas-storage combined dynamic economic dispatching optimization method based on Q learning
Zhu et al. Complementary operational research for a hydro-wind-solar hybrid power system on the upper Jinsha River
CN103944175A (en) Wind-solar-storage combined power generation system output characteristic optimization method
CN111585310B (en) Method and device for predicting output of distributed power supply
CN106712033A (en) Wind curtailment absorption method in thermal power plant
CN112217232B (en) Photovoltaic and photo-thermal coupling power generation system and power generation control method
CN114066204A (en) Integrated optimization planning and operation method and device of comprehensive energy system
CN106600022B (en) Wind-light-gas-seawater pumped storage isolated power system capacity optimal configuration method based on multi-objective optimization
CN108667077A (en) A kind of wind storage association system Optimization Scheduling
CN110198052B (en) Photo-thermal-wind power combined grid-connected power generation coordination control method
CN107102545B (en) Load optimization distribution and residual gas amount correction calculation method and system for shutdown process of combined cycle unit
CN112214883A (en) Micro-energy-grid comprehensive energy system operation optimization method and system
CN109376406B (en) Energy supply system superstructure model, modeling method, computer device and storage medium
CN112686425B (en) Energy internet optimal scheduling method and system based on cooperative game
CN114386256A (en) Regional electric heating system optimal scheduling method considering flexibility constraint of electric heating equipment and heat supply network characteristics
Qian et al. Research on source-load coordinated dispatch based on multi-scenario wind power consumption
CN113410850A (en) Light-heat wind-power combined frequency modulation model and frequency modulation strategy based on MPC
Xiao et al. Complementary coordinated dispatch of multi-energy systems considering power flexibility margin

Legal Events

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