CN109325679A - Consider that the integrated energy system of integration requirement response linearizes accidental scheduling method - Google Patents
Consider that the integrated energy system of integration requirement response linearizes accidental scheduling method Download PDFInfo
- Publication number
- CN109325679A CN109325679A CN201811057112.5A CN201811057112A CN109325679A CN 109325679 A CN109325679 A CN 109325679A CN 201811057112 A CN201811057112 A CN 201811057112A CN 109325679 A CN109325679 A CN 109325679A
- Authority
- CN
- China
- Prior art keywords
- energy
- scene
- gas
- constraint
- hinge
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
- 230000004044 response Effects 0.000 title claims abstract description 47
- 238000000034 method Methods 0.000 title claims abstract description 22
- 230000010354 integration Effects 0.000 title claims abstract description 17
- 230000005611 electricity Effects 0.000 claims abstract description 30
- 238000004146 energy storage Methods 0.000 claims abstract description 28
- 230000011218 segmentation Effects 0.000 claims abstract description 23
- 230000008878 coupling Effects 0.000 claims abstract description 7
- 238000010168 coupling process Methods 0.000 claims abstract description 7
- 238000005859 coupling reaction Methods 0.000 claims abstract description 7
- VNWKTOKETHGBQD-UHFFFAOYSA-N methane Chemical group C VNWKTOKETHGBQD-UHFFFAOYSA-N 0.000 claims description 63
- 239000007789 gas Substances 0.000 claims description 52
- 238000009825 accumulation Methods 0.000 claims description 45
- 239000003345 natural gas Substances 0.000 claims description 29
- 238000012797 qualification Methods 0.000 claims description 15
- 239000011159 matrix material Substances 0.000 claims description 9
- 230000008859 change Effects 0.000 claims description 7
- 238000006073 displacement reaction Methods 0.000 claims description 6
- 238000010943 off-gassing Methods 0.000 claims description 6
- 238000006243 chemical reaction Methods 0.000 claims description 5
- 238000005273 aeration Methods 0.000 claims description 3
- 238000005338 heat storage Methods 0.000 claims description 3
- 238000010438 heat treatment Methods 0.000 claims description 3
- 230000007704 transition Effects 0.000 claims description 3
- 239000013256 coordination polymer Substances 0.000 claims 1
- 230000000694 effects Effects 0.000 claims 1
- 230000009467 reduction Effects 0.000 description 3
- 230000005540 biological transmission Effects 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000008901 benefit Effects 0.000 description 1
- 238000004140 cleaning Methods 0.000 description 1
- 239000000567 combustion gas Substances 0.000 description 1
- 230000005684 electric field Effects 0.000 description 1
- 230000008569 process Effects 0.000 description 1
- 230000000087 stabilizing effect Effects 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q10/00—Administration; Management
- G06Q10/06—Resources, workflows, human or project management; Enterprise or organisation planning; Enterprise or organisation modelling
- G06Q10/063—Operations research, analysis or management
- G06Q10/0639—Performance analysis of employees; Performance analysis of enterprise or organisation operations
- G06Q10/06393—Score-carding, benchmarking or key performance indicator [KPI] analysis
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q50/00—Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
- G06Q50/06—Energy or water supply
-
- 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
- Y04—INFORMATION OR COMMUNICATION TECHNOLOGIES HAVING AN IMPACT ON OTHER TECHNOLOGY AREAS
- Y04S—SYSTEMS INTEGRATING TECHNOLOGIES RELATED TO POWER NETWORK OPERATION, COMMUNICATION OR INFORMATION TECHNOLOGIES FOR IMPROVING THE ELECTRICAL POWER GENERATION, TRANSMISSION, DISTRIBUTION, MANAGEMENT OR USAGE, i.e. SMART GRIDS
- Y04S10/00—Systems supporting electrical power generation, transmission or distribution
- Y04S10/50—Systems or methods supporting the power network operation or management, involving a certain degree of interaction with the load-side end user applications
-
- 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
- Y04—INFORMATION OR COMMUNICATION TECHNOLOGIES HAVING AN IMPACT ON OTHER TECHNOLOGY AREAS
- Y04S—SYSTEMS INTEGRATING TECHNOLOGIES RELATED TO POWER NETWORK OPERATION, COMMUNICATION OR INFORMATION TECHNOLOGIES FOR IMPROVING THE ELECTRICAL POWER GENERATION, TRANSMISSION, DISTRIBUTION, MANAGEMENT OR USAGE, i.e. SMART GRIDS
- Y04S50/00—Market activities related to the operation of systems integrating technologies related to power network operation or related to communication or information technologies
- Y04S50/16—Energy services, e.g. dispersed generation or demand or load or energy savings aggregation
Landscapes
- Business, Economics & Management (AREA)
- Human Resources & Organizations (AREA)
- Engineering & Computer Science (AREA)
- Economics (AREA)
- Strategic Management (AREA)
- Theoretical Computer Science (AREA)
- Entrepreneurship & Innovation (AREA)
- Educational Administration (AREA)
- Marketing (AREA)
- Development Economics (AREA)
- Health & Medical Sciences (AREA)
- Tourism & Hospitality (AREA)
- Physics & Mathematics (AREA)
- General Business, Economics & Management (AREA)
- General Physics & Mathematics (AREA)
- Public Health (AREA)
- Primary Health Care (AREA)
- Water Supply & Treatment (AREA)
- General Health & Medical Sciences (AREA)
- Game Theory and Decision Science (AREA)
- Operations Research (AREA)
- Quality & Reliability (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
The invention discloses a kind of integrated energy systems of consideration integration requirement response to linearize accidental scheduling method.The present invention determines objective function first, comprehensive energy network constraint is set again, including Unit Combination constraint, electric network constraint, gas network constraint, distributed energy storage constraint and coupling constraint, then integration requirement response model is set, including setting energy hinge internal demands response model and external demand response model, wherein external demand response model include can time shift load and interruptible load；Method is cut down with reversed wind-powered electricity generation again and generates N number of random scene, is simplified the quadratic term in Weymouth equation using contiguous segmentation linearized algorithm, is finally solved entire model with commercialization solver YALMIP.The present invention establishes the model of the hot integrated energy system of electric-gas-, formulates a unit output plan；There is provided it is a kind of more comprehensively, economic, reliable comprehensive energy network operation dispatching method.
Description
Technical field
The present invention relates to the accidental scheduling methods of integrated energy system, belong to integrated energy system field, and in particular to one
Kind considers that the integrated energy system of integration requirement response linearizes accidental scheduling method.
Background technique
Natural gas have the characteristics that cleaning, at a low price, it is environmental-friendly, and natural gas unit have rapid starting/stopping, be easy to store
Characteristic, this make the hot integrated energy system of the electric-gas-for having merged natural gas system have boundless application prospect.?
Couple in one region multiple energy hinges integrated energy system can both make full use of cogeneration unit (CHP) in hinge,
Electric boiler (EB), heat accumulation element (HS) etc. carry out the conversion of hinge internal power source, and can by energy source route carry out hinge it
Between energy transmitting, greatly improve utilization rate of each form energy in integrated energy system, improve running efficiency of system.
The bootable user's rationality energy of demand response rationally shifts Peak power use to translate load peak valley, to promotion electric power
Running efficiency of system is of great significance.
Wind power output has stronger uncertainty, and electric power system dispatching need to fully consider this fluctuation to electric system
Reliability bring influences, and accurately portrays the process that wind-powered electricity generation is contributed at random.
In the integrated energy system linearisation accidental scheduling method for considering integration requirement response, different wind are not only featured
Under electric field scape the case where wind power output, it is also contemplated that energy hinge inside and outside demand response is to the work for stabilizing load curve
With, operation plan a few days ago determine on the basis of, the advantage of the flexible start and stop of natural gas unit is made full use of, with a kind of piecewise linearity
Change method solves the model,
Summary of the invention
The problem to be solved in the present invention is, uncertain in face of wind power output, in power transmission network level, formulates a unit
Power output plan, so that including Unit Commitment expense, coal-fired expense, combustion gas expense, demand response expense, wind-powered electricity generation is cut down and cutting load
Rejection penalty including system overall running cost it is minimum.The present invention determines objective function first, then comprehensive energy network is arranged
Constraint, including Unit Combination constraint, electric network constraint, gas network constraint, distributed energy storage constraint and coupling constraint, are then arranged
Integration requirement response model, including setting energy hinge internal demands response model and external demand response model, wherein external
Demand response model includes can time shift load and interruptible load；Method is cut down with reversed wind-powered electricity generation again and generates N number of random scene, is used
Contiguous segmentation linearized algorithm simplifies the quadratic term in Weymouth equation, finally solves entire mould with commercialization solver YALMIP
Type.
Determine objective function, specific as follows:
Wherein, PrSIt is the probability of wind-powered electricity generation scene s,It is the expense of demand response in scene s,It is machine in scene s
Group i t moment unlatching expense,It is that unit i in t moment shuts down expense in scene s,It is that unit i exists in scene s
The generated energy of t moment,It isCorresponding cost of electricity-generating, ρGASIt is natural gas unit price,It is gas well sp in scene s
Displacement under t moment, ω are the rejection penalties that unit wind-powered electricity generation is cut down,It is wind of the blower w under t moment in scene s
Electric reduction, α are the rejection penalties of unit cutting load, and { } form of energy can represent electricity, natural gas or heat,It is scene
In s in energy hinge r all kinds of form of energy cutting load amount.
Comprehensive energy network constraint is set, is implemented as follows:
2-1. is arranged Unit Combination and constrains；
Climbed downwards about according to decision unit output bound, the upward Climing constant of unit, unit the characteristics of Unit Combination
Beam, unit starting expense, unit shut down expense, the constraint of unit minimum continuous operating time and unit minimum continuous idle time
Constraint；
2-2. is arranged electric network and constrains；
Gone out according to decision network trend constraint, node power Constraints of Equilibrium, phase angle bound and wind-powered electricity generation the characteristics of network
Power bound；
Gas network constraint is arranged in 2-3.；
1.
It is the displacement of gas well sp in scene s under t moment；
2.
It is the pressure of moral point m in scene s under t moment；
3.
It is the average value of the both ends gas pipeline mn throughput in scene s under t moment,It is t moment in scene s
Under；
4.
It is the natural gas load of moral point under t moment in scene s,WithIt is t moment in scene s respectively
The outlet end of lower pipeline mn and the throughput of arrival end,WithIt is that natural gas storage fills under t moment in scene s respectively
Set the release quantity and aeration quantity of q；The gas trend balance of constraint qualification moral point；
5.
It is the gas quantity of Gas Generator Set i under t moment in scene s,WithIt is unit i in scene s respectively
T moment unlatching and shut down air consumption, γ is that electricity is converted to the energy conversion factor of gas, and unit is m3/MW,It is scene
In s under t moment energy hinge r air inflow, the constraint qualification composition of moral point load；
6.
It is the gas-storing capacity of pipeline mn under t moment in scene s,It is the gas storage coefficient of discharge of pipeline mn,It is field
In scape s under t moment, the average value of pipeline mn two-port pressure, the constraint qualification pass of pipeline mn gas-storing capacity and both ends pressure
System；
7.
Constraint qualification pipeline mn gas-storing capacity and both ends enter and leave the time coupled relation of throughput；
2-4. is arranged distributed storage and constrains；
①.
It is the storage state of energy-storage travelling wave tube q under t moment in scene s, 0 expression is taken to release in energy, takes in 1 expression energy storage；It is the maximum energy storage value of energy-storage travelling wave tube q；The restrict energy storage upper limit of energy-storage travelling wave tube；
②.
It is that release can value for the maximum of energy-storage travelling wave tube q；Releasing for restrict energy-storage travelling wave tube can the upper limit；
③.
It is the gas-storing capacity of reservation component q under t moment in scene s,WithIt is reservation component q inflation and outgassing respectively
Efficiency；Time coupled relation of the constraint qualification between reservation component gas-storing capacity and gas storage outgassing；
④.
WithIt is the minimum value and maximum value of reservation component q gas-storing capacity respectively；Restrict reservation component
The bound of gas-storing capacity；
⑤.
It is the gas-storing capacity of reservation component q under the last one period NT in scene s, SGq,0It is that reservation component q is in office
Gas-storing capacity under one scene under first period 0；Constraint qualification reservation component is before and after the cycle of operation
Gas-storing capacity should be constant；
Coupling constraint is arranged in 2-5.；
⑴.
The operating condition for indicating the unit i under t moment in scene s takes in the operation of 1 expression unit, 0 expression unit is taken to close
In stopping；Unit Combination arrangement under the different wind-powered electricity generation scenes of the restrict should be the same；
⑵.
It is Gas Generator Set energy conversion efficiency, unit is m3/MW；Constraint qualification Gas Generator Set can turn fossil
It is changed to the energy relationship of electric energy.
Integration requirement response model is set, is implemented as follows:
Energy hinge internal demands response model is arranged in 3-1.
Inside energy hinge, there are cogeneration unit, electric boiler and heat to store the energy transition equipment of these types of form；
IfWithIt is the input quantity of energy hinge electric power and natural gas respectively,WithBe respectively energy hinge electric power,
The output quantity η of natural gas and heating powerCEAnd ηCHIt is the efficiency that natural gas is converted to electric energy and thermal energy by cogeneration unit respectively,
ηHCAnd ηHDIt is the efficiency that heat accumulation element fills heat and heat release, η respectivelyEBIt is the transfer efficiency that electric boiler converts the electricity into thermal energy；The energy
Hinge byWithEnergy input forms v inside hinge1~v11Energy stream, it is respective between there is also the constraint relationships；
Energy hinge energy input and output relation are constrained:
①
WithIt is the input and output value of energy hinge r energy stream in s scene, C respectivelyrIt is the energy of energy hinge r
Amount stream input/output relation matrix；
②
Indicate the change value of the heat accumulation element quantity of heat storage under t moment in scene s；The matrix indicates that energy hinge r's is outer
Portion's input power flow and internal energy stream together constitute the energy stream input of energy hinge；
③
The matrix indicates that energy hinge output quantity is made of various forms of energy streams；
④
Indicate in scene s under t moment in the energy hinge r heat accumulation element storage state, take 1 expression to fill and hanker,
It takes in 0 expression heat release；Indicate that the maximum of heat accumulation element in energy hinge r fills heat；The restrict is by v6And v10Group
At energy stream as heat accumulation element HS input the upper limit；
⑤
Indicate the maximum heat release amount of heat accumulation element in energy hinge r；The restrict is by v11The energy stream of composition
The upper limit of output as heat accumulation element；
⑥
Indicate in scene s under t moment in the energy hinge r heat accumulation element heat accumulation value；Constraint qualification energy pivot
The time coupled relation of heat accumulation element heat accumulation value and thermal energy change value in knob r；
⑦
Indicate the maximum heat accumulation value of heat accumulation element in energy hinge r；Heat accumulation member in restrict energy hinge r
The bound of part heat accumulation value；
⑧
It is the maximum input level of cogeneration unit in energy hinge r；The restrict is by v5The energy of composition
Flow the upper limit of the input as cogeneration unit；
⑨
It is the maximum input level of electric boiler in energy hinge r；The restrict is by v2And v5The energy stream of composition
The upper limit as water-tube boiler input；
Energy hinge external demand response model is arranged in 3-2.
External demand response is carried out in terminal user side outside energy hinge, external demand response model is set:
For in the t period, energy hinge r responded without external demand before all kinds of energy demands,For the t at scene s
Period, all kinds of energy form demands of the energy hinge r after external demand responds.
Wherein, can time shift load have following constraint:
Wherein, interruptible load has following constraint:
Contiguous segmentation linearized algorithm is set, is implemented as follows:
The initial precision A of algorithm is arranged in 5-1.f=1, Ap=1, accuracy threshold values is set
The segments k of 5-2. setting natural gas trendfWith node air pressure kp；
5-3. is in existing solution segmentation, by natural gas trend quadratic termLinearisation；
Existing natural gas trend section is further subdivided into Kf sections；Wherein cmn,lAnd bmn,lIt is in first of solution segmentation respectively
Slope and intercept, δmn,lIt is the existing selection variable for solving segmentation, takes 1 expression solution to fall in first of section, take 0 expression solution not
Fall in first of section；
5-4. calculates the model problem after linearisation with solver YALMIP, finds rough in the Linear Segmentation
Solution, and be that new existing solution is segmented by the subsection setup where the rough solution；
5-5. updates arithmetic accuracy Af=kf*Af, Ap=kp*Ap, and check whether the precision is more than precision threshold values, i.e.,Such as the inequality is set up, then this contiguous segmentation linearized algorithm terminates；It is such as invalid, then it repeats
Step 5-2 to 5-5.
The present invention has the beneficial effect that:
The present invention establishes the model of the hot integrated energy system of electric-gas-, comprehensive energy hinge is added in a model, in energy
Setting cogeneration unit, electric boiler and heat-storing device inside source hinge, be arranged outside the energy hinge electricity, gas transmission line and
Natural gas storage device；The internal demands response that the model is realized by the internal energy conversion of energy hinge, passes through energy pivot
Knob terminal user setting can time shift load and interruptible load realize the external demand response of the model；There is provided a kind of continuous point
Section linearized algorithm solves the model.To sum up, this method can provide it is a kind of more comprehensively, economic, reliable comprehensive energy network fortune
Row dispatching method.
Detailed description of the invention
Fig. 1 is flow chart of the present invention.
Specific embodiment
Present invention will be further explained below with reference to the attached drawings and examples.
As shown in Figure 1, a kind of integrated energy system for considering integration requirement response linearizes accidental scheduling method.First really
Set the goal function, and Unit Combination constraint, electric network constraint, gas network constraint, distributed energy storage constraint and coupling is then arranged about
Beam, be then arranged internal demands respond and including can time shift load and interruptible load external demand response, then with reversed wind
Electric reduction method generates N number of random scene, simplifies the quadratic term in Weymouth equation using contiguous segmentation linearized algorithm, finally
Entire model is solved with commercial solver YALMIP.
In order to solve the technical problem, solution of the invention the following steps are included:
Objective function is arranged in step (1)；
Wherein, PrSIt is the probability of wind-powered electricity generation scene s,It is the expense of demand response in scene s,It is machine in scene s
Group i t moment unlatching expense,It is that unit i in t moment shuts down expense in scene s,It is that unit i exists in scene s
The generated energy of t moment,It isCorresponding cost of electricity-generating, ρGASIt is natural gas unit price,It is gas well sp in scene s
Displacement under t moment, ω are the rejection penalties that unit wind-powered electricity generation is cut down,It is wind of the blower w under t moment in scene s
Electric reduction, α are the rejection penalties of unit cutting load, and { } form of energy can represent electricity, natural gas or heat,It is scene
In s in energy hinge r all kinds of form of energy cutting load amount.
Comprehensive energy network constraint is arranged in step (2)
2-1. is arranged Unit Combination and constrains；
Climbed downwards about according to decision unit output bound, the upward Climing constant of unit, unit the characteristics of Unit Combination
Beam, unit starting expense, unit shut down expense, the constraint of unit minimum continuous operating time and unit minimum continuous idle time
Constraint.
2-2. is arranged electric network and constrains；
Gone out according to decision network trend constraint, node power Constraints of Equilibrium, phase angle bound and wind-powered electricity generation the characteristics of network
Power bound.
Gas network constraint is arranged in 2-3.；
1.
It is the displacement of gas well sp in scene s under t moment.Restrict gas well gas output.
2.
It is the pressure of moral point m in scene s under t moment.The restrict value of moral point pressure.
3.
It is the average value of the both ends gas pipeline mn throughput in scene s under t moment,It is t moment in scene s
Under, the characterisitic parameter of gas pipeline mn is related to duct length.This be constrained to Weymouth equation, it is specified that chimneying with
Relationship between node air pressure.
4.
It is the natural gas load of moral point under t moment in scene s,WithIt is t moment in scene s respectively
The outlet end of lower pipeline mn and the throughput of arrival end,WithIt is that natural gas storage fills under t moment in scene s respectively
Set the release quantity and aeration quantity of q.The gas trend balance of constraint representation moral point.
5.
It is the gas quantity of Gas Generator Set i under t moment in scene s,WithIt is unit i in scene s respectively
T moment unlatching and shut down air consumption, γ is that electricity is converted to the energy conversion factor of gas, and unit is m3/MW,It is scene
In s under t moment energy hinge r air inflow.The constraint representation composition of moral point load.
6.
It is the gas-storing capacity of pipeline mn under t moment in scene s,It is the gas storage coefficient of discharge of pipeline mn,It is field
In scape s under t moment, the average value of pipeline mn two-port pressure.The constraint representation pass of pipeline mn gas-storing capacity and both ends pressure
System.
7.
Constraint representation pipeline mn gas-storing capacity and both ends enter and leave the time coupled relation of throughput.
2-4. is arranged distributed storage and constrains；
①.
It is the storage state of energy-storage travelling wave tube q under t moment in scene s, 0 expression is taken to release in energy, takes in 1 expression energy storage.It is the maximum energy storage value of energy-storage travelling wave tube q.The restrict energy storage upper limit of the power of energy-storage travelling wave tube.
②.
It is that release can value for the maximum of energy-storage travelling wave tube q.Releasing for restrict energy-storage travelling wave tube can the upper limit.
③.
It is the gas-storing capacity of reservation component q under t moment in scene s,WithIt is reservation component q inflation and outgassing respectively
Efficiency.Time coupled relation of the constraint qualification between reservation component gas-storing capacity and gas storage outgassing.
④.
WithIt is the minimum value and maximum value of reservation component q gas-storing capacity respectively.Restrict reservation component
The bound of gas-storing capacity.
⑤.
It is the gas-storing capacity of reservation component q under the last one period NT in scene s, SGq,0It is that reservation component q is in office
Gas-storing capacity under one scene under first period 0.Constraint qualification reservation component is before and after the cycle of operation
Gas-storing capacity should be constant.
Coupling constraint is arranged in 2-5.；
⑴.
The operating condition for indicating the unit i under t moment in scene s takes in the operation of 1 expression unit, 0 expression unit is taken to close
In stopping.Unit Combination arrangement under the different wind-powered electricity generation scenes of the restrict should be the same.
⑵.
It is Gas Generator Set energy conversion efficiency, unit is m3/MW.Constraint representation Gas Generator Set can turn fossil
It is changed to the energy relationship of electric energy.
Integration requirement response model is arranged in step (3)；
Energy hinge internal demands response model is arranged in 3-1.
As shown in Fig. 1, inside energy hinge, cogeneration unit, electric boiler and heat store these types of form
Energy transition equipment.In figure,WithIt is the input quantity of energy hinge electric power and natural gas respectively,WithPoint
It is not the output quantity of energy hinge electric power, natural gas and heating power, ηCEAnd ηCHIt is that cogeneration unit is converted to natural gas respectively
The efficiency of electric energy and thermal energy, ηHCAnd ηHDIt is the efficiency that heat accumulation element fills heat and heat release, η respectivelyEBIt is that electric boiler converts the electricity into heat
The transfer efficiency of energy.As shown in figure, energy hinge byWithEnergy input forms v inside hinge1~v11Energy stream,
Its respectively between there is also the constraint relationships.
Energy hinge energy input and output relation are constrained:
①
WithIt is the input and output value of energy hinge r energy stream in s scene, C respectivelyrIt is the energy of energy hinge r
Amount stream input/output relation matrix.
②
Indicate the change value of the heat accumulation element quantity of heat storage under t moment in scene s.The matrix indicates that energy hinge r's is outer
Portion's input power flow and internal energy stream together constitute the energy stream input of energy hinge.
③
The matrix indicates that energy hinge output quantity is made of various forms of energy streams.
④
Indicate in scene s under t moment in the energy hinge r heat accumulation element storage state, take 1 expression to fill and hanker,
It takes in 0 expression heat release.Indicate that the maximum of heat accumulation element in energy hinge r fills heat.The restrict is by v6And v10Group
At energy stream as heat accumulation element HS input the upper limit.
⑤
Indicate the maximum heat release amount of heat accumulation element in energy hinge r.The restrict is by v11The energy stream of composition
The upper limit of output as heat accumulation element.
⑥
Indicate in scene s under t moment in the energy hinge r heat accumulation element heat accumulation value.The constraint qualification energy
The time coupled relation of heat accumulation element heat accumulation value and thermal energy change value in hinge r.
⑦
Indicate the maximum heat accumulation value of heat accumulation element in energy hinge r.Heat accumulation member in restrict energy hinge r
The bound of part heat accumulation value.
⑧
It is the maximum input level of cogeneration unit in energy hinge r.The restrict is by v5The energy of composition
Flow the upper limit of the input as cogeneration unit.
⑨
It is the maximum input level of electric boiler in energy hinge r.The restrict is by v2And v5The energy stream of composition is made
For the upper limit of water-tube boiler input.
Energy hinge external demand response model is arranged in 3-2.
External demand response is carried out in terminal user side outside energy hinge, external demand response model is set:
For in the t period, energy hinge r responded without external demand before all kinds of energy demands,For the t at scene s
Period, all kinds of energy form demands of the energy hinge r after external demand responds.
Wherein, can time shift load have following constraint:
Wherein, interruptible load has following constraint:
Wind-powered electricity generation random scene is arranged in step (4)；
According to the uncertain characteristic of wind-powered electricity generation, N number of wind-powered electricity generation random scene is generated with reversed scene "flop-out" method.
Contiguous segmentation linearized algorithm is arranged in step (5)；
The initial precision A of algorithm is arranged in 5-1.f=1, Ap=1, accuracy threshold values is set
The segments k of 5-2. setting natural gas trendfWith node air pressure kp。
5-3. is in existing solution segmentation, by natural gas trend quadratic termLinearisation.
Existing natural gas trend section is further subdivided into Kf sections；Wherein cmn,lAnd bmn,lIt is in first of solution segmentation respectively
Slope and intercept, δmn,lIt is the existing selection variable for solving segmentation, takes 1 expression solution to fall in first of section, take 0 expression Xie Weiluo
In first of section.
5-4. calculates the model problem after linearisation with solver YALMIP, finds rough in the Linear Segmentation
Solution, and be that new existing solution is segmented by the subsection setup where the rough solution.
5-5. updates arithmetic accuracy Af=kf*Af, Ap=kp*Ap, and check whether the precision is more than precision threshold values, i.e.,Such as the inequality is set up, then this contiguous segmentation linearized algorithm terminates；It is such as invalid, then it repeats
Step 5-2 to 5-5.
Step (6) solves entire model.
Due to the linearisation by step (5), entire model has been converted into a MIXED INTEGER linear problem, can use commercialization
Solver YALMIP, which is solved, to be calculated.
Claims (5)
1. considering that the integrated energy system of integration requirement response linearizes accidental scheduling method, it is characterised in that determine target first
Function, then comprehensive energy network constraint is set, including Unit Combination constraint, electric network constraint, gas network constraint, distributed energy storage
Constraint and coupling constraint, are then arranged integration requirement response model, including setting energy hinge internal demands response model and outer
Portion's demand response model, wherein external demand response model include can time shift load and interruptible load；It is cut again with reversed wind-powered electricity generation
Subtraction generates N number of random scene, simplifies the quadratic term in Weymouth equation using contiguous segmentation linearized algorithm, finally uses quotient
Entire model is solved with solver YALMIP.
2. the integrated energy system according to claim 1 for considering integration requirement response linearizes accidental scheduling method,
It is characterized in that determining objective function, specific as follows:
Wherein, PrSIt is the probability of wind-powered electricity generation scene s,It is the expense of demand response in scene s,It is unit i in scene s
In the unlatching expense of t moment,It is that unit i in t moment shuts down expense in scene s,Be in scene s unit i in t
The generated energy at quarter,It isCorresponding cost of electricity-generating, ρGASIt is natural gas unit price,Be in scene s gas well sp in t
The displacement inscribed, ω are the rejection penalties that unit wind-powered electricity generation is cut down,It is that wind-powered electricity generation of the blower w under t moment is cut down in scene s
Amount, α are the rejection penalties of unit cutting load, and { } form of energy can represent electricity, natural gas or heat,It is energy in scene s
The cutting load amount of all kinds of form of energy in the hinge r of source.
3. the integrated energy system according to claim 2 for considering integration requirement response linearizes accidental scheduling method,
It is characterized in that the setting comprehensive energy network constraint, is implemented as follows:
2-1. is arranged Unit Combination and constrains；
According to decision unit output bound, the upward Climing constant of unit, the downward Climing constant of unit, machine the characteristics of Unit Combination
Group starting expense, unit shut down expense, the constraint of unit minimum continuous operating time and the minimum continuous idle time constraint of unit；
2-2. is arranged electric network and constrains；
It is determined in network trend constraint, node power Constraints of Equilibrium, phase angle bound and wind power output according to the characteristics of network
Lower limit；
Gas network constraint is arranged in 2-3.；
1.
It is the displacement of gas well sp in scene s under t moment；Restrict gas well gas output；
2.
It is the pressure of moral point m in scene s under t moment；The restrict value of moral point pressure；
3.
It is the average value of the both ends gas pipeline mn throughput in scene s under t moment,It is in scene s under t moment, it is defeated
The characterisitic parameter of feed channel mn is related to duct length；This is constrained to Weymouth equation, it is specified that chimneying and node gas
Relationship between pressure；
4.
It is the natural gas load of moral point under t moment in scene s,WithIt is pipeline under t moment in scene s respectively
The outlet end of mn and the throughput of arrival end,WithIt is that natural gas storage device q is released under t moment in scene s respectively
Tolerance and aeration quantity；The gas trend balance of constraint qualification moral point；
5.
It is the gas quantity of Gas Generator Set i under t moment in scene s,WithBe respectively in scene s unit i in t
The unlatching at quarter and air consumption is shut down, γ is the energy conversion factor that electricity is converted to gas, and unit is m3/MW,It is t in scene s
When inscribe the air inflow of energy hinge r, the constraint representation composition of moral point load；
6.
It is the gas-storing capacity of pipeline mn under t moment in scene s,It is the gas storage coefficient of discharge of pipeline mn,It is t in scene s
When inscribe, the average value of pipeline mn two-port pressure, the constraint representation relationship of pipeline mn gas-storing capacity and both ends pressure；
7.
Constraint representation pipeline mn gas-storing capacity and both ends enter and leave the time coupled relation of throughput；
2-4. is arranged distributed storage and constrains；
①.
It is the storage state of energy-storage travelling wave tube q under t moment in scene s, 0 expression is taken to release in energy, takes in 1 expression energy storage；
It is the maximum energy storage value of energy-storage travelling wave tube q；The restrict energy storage upper limit of the power of energy-storage travelling wave tube；
②.
It is that release can value for the maximum of energy-storage travelling wave tube q；Releasing for restrict energy-storage travelling wave tube can the upper limit；
③.
It is the gas-storing capacity of reservation component q under t moment in scene s,WithIt is the effect of reservation component q inflation and outgassing respectively
Rate；Time coupled relation of the constraint qualification between reservation component gas-storing capacity and gas storage outgassing；
④.
WithIt is the minimum value and maximum value of reservation component q gas-storing capacity respectively；Restrict reservation component gas storage
The bound of amount；
⑤.
It is the gas-storing capacity of reservation component q under the last one period NT in scene s, SGq,0It is reservation component q in any field
Gas-storing capacity under scape under first period 0；Constraint qualification gas storage of the reservation component before and after a cycle of operation
Amount should be constant；
Coupling constraint is arranged in 2-5.；
⑴.
The operating condition for indicating the unit i under t moment in scene s takes in the operation of 1 expression unit, 0 expression unit is taken to shut down
In；Unit Combination arrangement under the different wind-powered electricity generation scenes of the restrict should be the same；
⑵.
It is Gas Generator Set energy conversion efficiency, unit is m3/MW；Constraint representation Gas Generator Set can be converted to fossil
The energy relationship of electric energy.
4. the integrated energy system according to claim 3 for considering integration requirement response linearizes accidental scheduling method,
It is characterized in that the setting integration requirement response model, is implemented as follows:
Energy hinge internal demands response model is arranged in 3-1.
Inside energy hinge, there are cogeneration unit, electric boiler and heat to store the energy transition equipment of these types of form；If
WithIt is the input quantity of energy hinge electric power and natural gas respectively,WithIt is energy hinge electric power respectively, natural
The output quantity of gas and heating power, ηCEAnd ηCHIt is the efficiency that natural gas is converted to electric energy and thermal energy by cogeneration unit, η respectivelyHCWith
ηHDIt is the efficiency that heat accumulation element fills heat and heat release, η respectivelyEBIt is the transfer efficiency that electric boiler converts the electricity into thermal energy；Energy hinge
ByWithEnergy input forms v inside hinge1~v11Energy stream, it is respective between there is also the constraint relationships；
Energy hinge energy input and output relation are constrained:
①
WithIt is the input and output value of energy hinge r energy stream in s scene, C respectivelyrBe energy hinge r energy stream it is defeated
Enter output relation matrix；
②
Indicate the change value of the heat accumulation element quantity of heat storage under t moment in scene s；The matrix indicates that the outside of energy hinge r is defeated
Enter energy stream and internal energy stream together constitutes the energy stream input of energy hinge；
③
The matrix indicates that energy hinge output quantity is made of various forms of energy streams；
④
Indicate in scene s under t moment in the energy hinge r heat accumulation element storage state, take 1 expression to fill and hanker, take 0 table
Show in heat release；Indicate that the maximum of heat accumulation element in energy hinge r fills heat；The restrict is by v6And v10The energy of composition
The upper limit of the amount stream as the input of heat accumulation element HS；
⑤
Indicate the maximum heat release amount of heat accumulation element in energy hinge r；The restrict is by v11The energy stream conduct of composition
The upper limit of the output of heat accumulation element；
⑥
Indicate in scene s under t moment in the energy hinge r heat accumulation element heat accumulation value；Constraint qualification energy hinge r
The time coupled relation of middle heat accumulation element heat accumulation value and thermal energy change value；
⑦
Indicate the maximum heat accumulation value of heat accumulation element in energy hinge r；Heat accumulation element storage in restrict energy hinge r
The bound of calorific value；
⑧
CPr maxIt is the maximum input level of cogeneration unit in energy hinge r；The restrict is by v5The energy stream of composition is made
For the upper limit of the input of cogeneration unit；
⑨
It is the maximum input level of electric boiler in energy hinge r；The restrict is by v2And v5The energy stream of composition is as heat
The upper limit of electric boiler input；
Energy hinge external demand response model is arranged in 3-2.
External demand response is carried out in terminal user side outside energy hinge, external demand response model is set:
For in the t period, energy hinge r responded without external demand before all kinds of energy demands,For at scene s
The t period, all kinds of energy form demands of the energy hinge r after external demand responds.
Wherein, can time shift load have following constraint:
Wherein, interruptible load has following constraint:
5. the integrated energy system according to claim 4 for considering integration requirement response linearizes accidental scheduling method,
It is characterized in that setting contiguous segmentation linearized algorithm, is implemented as follows:
The initial precision A of algorithm is arranged in 5-1.f=1, Ap=1, accuracy threshold values is set
The segments k of 5-2. setting natural gas trendfWith node air pressure kp；
5-3. is in existing solution segmentation, by natural gas trend quadratic termLinearisation；
Existing natural gas trend section is further subdivided into Kf sections；Wherein cmn,lAnd bmn,lIt is oblique in first of solution segmentation respectively
Rate and intercept, δmn,lIt is the existing selection variable for solving segmentation, takes 1 expression solution to fall in first of section, 0 expression solution is taken not fall within
First of section；
5-4. calculates the model problem after linearisation with solver YALMIP, finds the rough solution in the Linear Segmentation, and
It is new existing solution segmentation by the subsection setup where the rough solution；
5-5. updates arithmetic accuracy Af=kf*Af, Ap=kp*Ap, and check whether the precision is more than precision threshold values, i.e.,Such as the inequality is set up, then this contiguous segmentation linearized algorithm terminates；It is such as invalid, then it repeats
Step 5-2 to 5-5.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201811057112.5A CN109325679B (en) | 2018-09-11 | 2018-09-11 | Linear random scheduling method of comprehensive energy system considering comprehensive demand response |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201811057112.5A CN109325679B (en) | 2018-09-11 | 2018-09-11 | Linear random scheduling method of comprehensive energy system considering comprehensive demand response |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN109325679A true CN109325679A (en) | 2019-02-12 |
| CN109325679B CN109325679B (en) | 2022-04-05 |
Family
ID=65264952
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201811057112.5A Active CN109325679B (en) | 2018-09-11 | 2018-09-11 | Linear random scheduling method of comprehensive energy system considering comprehensive demand response |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN109325679B (en) |
Cited By (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110163443A (en) * | 2019-05-27 | 2019-08-23 | 西南石油大学 | Consider the micro- energy net Optimization Scheduling in the natural gas pressure regulating station of electric-gas integration requirement response |
| CN110377973A (en) * | 2019-06-28 | 2019-10-25 | 西安交通大学 | A kind of construction method of normal linearity integrated energy system model |
| CN110941799A (en) * | 2019-11-29 | 2020-03-31 | 国网辽宁省电力有限公司经济技术研究院 | Energy hub random planning method considering system comprehensive uncertainty factors |
| CN111831963A (en) * | 2020-07-14 | 2020-10-27 | 中国南方电网有限责任公司 | Method for evaluating controllability of comprehensive energy service provider under power market background |
| CN112381375A (en) * | 2020-11-09 | 2021-02-19 | 浙江大学 | Power grid economic operation domain rapid generation method based on power flow distribution matrix |
| CN113343497A (en) * | 2021-06-30 | 2021-09-03 | 广东电网有限责任公司 | Decoupling method and device for cooling, heating and power multifunctional coupling system |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN106356902A (en) * | 2016-11-01 | 2017-01-25 | 东南大学 | Energy management method applicable to multi-energy cooperation comprehensive energy network in Internet of energy |
| CN107257132A (en) * | 2017-05-17 | 2017-10-17 | 重庆大学 | A kind of electrical interconnection system synthetic load for considering that wind-powered electricity generation is abandoned cuts down model building method |
| CN107368927A (en) * | 2017-08-01 | 2017-11-21 | 重庆大学 | Electrical energy flow point cloth collaboration optimized calculation method based on target cascade analysis |
| CN108494015A (en) * | 2018-02-09 | 2018-09-04 | 中国科学院电工研究所 | The integrated energy system design method of one introduces a collection-lotus-storage coordination and interaction |
-
2018
- 2018-09-11 CN CN201811057112.5A patent/CN109325679B/en active Active
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN106356902A (en) * | 2016-11-01 | 2017-01-25 | 东南大学 | Energy management method applicable to multi-energy cooperation comprehensive energy network in Internet of energy |
| CN107257132A (en) * | 2017-05-17 | 2017-10-17 | 重庆大学 | A kind of electrical interconnection system synthetic load for considering that wind-powered electricity generation is abandoned cuts down model building method |
| CN107368927A (en) * | 2017-08-01 | 2017-11-21 | 重庆大学 | Electrical energy flow point cloth collaboration optimized calculation method based on target cascade analysis |
| CN108494015A (en) * | 2018-02-09 | 2018-09-04 | 中国科学院电工研究所 | The integrated energy system design method of one introduces a collection-lotus-storage coordination and interaction |
Non-Patent Citations (2)
| Title |
|---|
| MICHAEL URBINA ET.AL.: ""A Combined Model for Analyzing the Interdependency of Electrical and Gas Systems"", 《2007 39TH NORTH AMERICAN POWER SYMPOSIUM》 * |
| 仉梦林等: ""基于动态场景集和需求响应的二阶段随机规划调度模型"", 《电力系统自动化》 * |
Cited By (10)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110163443A (en) * | 2019-05-27 | 2019-08-23 | 西南石油大学 | Consider the micro- energy net Optimization Scheduling in the natural gas pressure regulating station of electric-gas integration requirement response |
| CN110163443B (en) * | 2019-05-27 | 2022-09-09 | 西南石油大学 | Natural gas pressure regulating station micro-energy network optimization scheduling method considering electricity-gas comprehensive demand response |
| CN110377973A (en) * | 2019-06-28 | 2019-10-25 | 西安交通大学 | A kind of construction method of normal linearity integrated energy system model |
| CN110377973B (en) * | 2019-06-28 | 2021-08-13 | 西安交通大学 | Construction method of standard linear comprehensive energy system model |
| CN110941799A (en) * | 2019-11-29 | 2020-03-31 | 国网辽宁省电力有限公司经济技术研究院 | Energy hub random planning method considering system comprehensive uncertainty factors |
| CN110941799B (en) * | 2019-11-29 | 2023-08-08 | 国网辽宁省电力有限公司经济技术研究院 | Energy hub stochastic programming method considering comprehensive uncertainty factors of system |
| CN111831963A (en) * | 2020-07-14 | 2020-10-27 | 中国南方电网有限责任公司 | Method for evaluating controllability of comprehensive energy service provider under power market background |
| CN112381375A (en) * | 2020-11-09 | 2021-02-19 | 浙江大学 | Power grid economic operation domain rapid generation method based on power flow distribution matrix |
| CN112381375B (en) * | 2020-11-09 | 2024-03-29 | 浙江大学 | Rapid generation method for power grid economic operation domain based on tide distribution matrix |
| CN113343497A (en) * | 2021-06-30 | 2021-09-03 | 广东电网有限责任公司 | Decoupling method and device for cooling, heating and power multifunctional coupling system |
Also Published As
| Publication number | Publication date |
|---|---|
| CN109325679B (en) | 2022-04-05 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN109325679A (en) | Consider that the integrated energy system of integration requirement response linearizes accidental scheduling method | |
| WO2019233134A1 (en) | Data-driven three-stage scheduling method for power-heat-gas grid based on wind power uncertainty | |
| CN105576657B (en) | Random fuzzy continuous power flow modeling method for active power distribution network and power transmission network integrated system | |
| CN109193815A (en) | A kind of combined heat and power dispatching method improving wind electricity digestion | |
| CN107944757A (en) | Electric power interacted system regenerative resource digestion capability analysis and assessment method | |
| CN108628176B (en) | Thermoelectric combined optimization scheduling method considering heat supply network flexibility | |
| CN105305419A (en) | Compressed air energy storage-containing independent microgrid capacity optimal configuration method | |
| CN111030094B (en) | Method for evaluating reliability of energy supply of multi-energy complementary park | |
| CN107257132B (en) | Method for constructing comprehensive load reduction model of electricity-gas interconnection system considering wind power abandonment | |
| CN104504451A (en) | Optimization method of day-ahead unit commitment of thermal power system of thermal power plant capable of operating under multiple states | |
| CN109523076A (en) | A kind of energy storage participate under integrated energy system steady-state optimization analysis method a few days ago | |
| CN104935017B (en) | Based on the wind-powered electricity generation and fired power generating unit combined method for improving light Robust Optimization Model | |
| CN113313351B (en) | Method for evaluating flexibility of electric-gas-heat system considering multi-energy coupling influence | |
| CN111985165B (en) | Electric heating comprehensive energy system scheduling method considering heat storage characteristics of intelligent building | |
| CN111564873B (en) | Integrated gas-electric energy system optimization method considering unit start and stop | |
| CN111047069B (en) | Thermoelectric combined scheduling flexibility assessment method considering energy storage characteristics of heat supply network | |
| CN107871052A (en) | A kind of meter and regenerative resource and the energy hub Optimal Operation Model of energy storage | |
| CN115936336B (en) | Virtual power plant capacity configuration and regulation operation optimization method | |
| CN112862158A (en) | Day-ahead optimization scheduling method and system for comprehensive energy system | |
| CN114387129B (en) | Comprehensive energy scheduling method, device, equipment and medium considering natural gas pipe network | |
| CN115169128A (en) | Comprehensive energy system unit combination method and system considering network dynamic characteristics | |
| Yuan et al. | A multi-energy flow calculation method considering multiple energy coupling operation modes | |
| CN114492917A (en) | Comprehensive energy optimization scheduling method and system based on dynamic characteristics of electric and thermal network | |
| CN113837589A (en) | N-1 scheduling method and device for comprehensive energy system | |
| CN113011714A (en) | Combined configuration method and system for cogeneration unit |
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 |