CN110516951A - A kind of integrated energy system dispatching method of dynamic interval - Google Patents
A kind of integrated energy system dispatching method of dynamic interval Download PDFInfo
- Publication number
- CN110516951A CN110516951A CN201910776326.6A CN201910776326A CN110516951A CN 110516951 A CN110516951 A CN 110516951A CN 201910776326 A CN201910776326 A CN 201910776326A CN 110516951 A CN110516951 A CN 110516951A
- Authority
- CN
- China
- Prior art keywords
- few days
- scheduling
- control
- integrated energy
- reference locus
- 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
- 238000000034 method Methods 0.000 title claims abstract description 52
- 238000006243 chemical reaction Methods 0.000 claims abstract description 13
- 239000011159 matrix material Substances 0.000 claims description 37
- 238000005457 optimization Methods 0.000 claims description 29
- 238000012937 correction Methods 0.000 claims description 24
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims description 17
- 230000008859 change Effects 0.000 claims description 15
- 238000005096 rolling process Methods 0.000 claims description 10
- 230000001276 controlling effect Effects 0.000 claims description 9
- 238000012546 transfer Methods 0.000 claims description 8
- 238000004519 manufacturing process Methods 0.000 claims description 7
- 230000033228 biological regulation Effects 0.000 claims description 6
- 230000001105 regulatory effect Effects 0.000 claims description 6
- 230000000694 effects Effects 0.000 claims description 4
- 230000004927 fusion Effects 0.000 claims description 3
- 230000001960 triggered effect Effects 0.000 claims description 2
- 238000010248 power generation Methods 0.000 claims 1
- 230000032258 transport Effects 0.000 claims 1
- 230000009286 beneficial effect Effects 0.000 abstract description 2
- VNWKTOKETHGBQD-UHFFFAOYSA-N methane Chemical compound C VNWKTOKETHGBQD-UHFFFAOYSA-N 0.000 description 38
- 239000007789 gas Substances 0.000 description 36
- 239000003345 natural gas Substances 0.000 description 19
- 238000010438 heat treatment Methods 0.000 description 9
- 230000005611 electricity Effects 0.000 description 8
- 238000011160 research Methods 0.000 description 6
- 238000012360 testing method Methods 0.000 description 6
- 230000008569 process Effects 0.000 description 5
- 230000001052 transient effect Effects 0.000 description 5
- 238000010586 diagram Methods 0.000 description 4
- 238000009826 distribution Methods 0.000 description 4
- 230000004044 response Effects 0.000 description 4
- 230000005540 biological transmission Effects 0.000 description 3
- 230000003321 amplification Effects 0.000 description 2
- 238000004364 calculation method Methods 0.000 description 2
- 238000001816 cooling Methods 0.000 description 2
- 230000008878 coupling Effects 0.000 description 2
- 238000010168 coupling process Methods 0.000 description 2
- 238000005859 coupling reaction Methods 0.000 description 2
- 235000013399 edible fruits Nutrition 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 238000003199 nucleic acid amplification method Methods 0.000 description 2
- 230000036961 partial effect Effects 0.000 description 2
- 238000005057 refrigeration Methods 0.000 description 2
- 241001672694 Citrus reticulata Species 0.000 description 1
- 241001269238 Data Species 0.000 description 1
- 238000010521 absorption reaction Methods 0.000 description 1
- 239000003795 chemical substances by application Substances 0.000 description 1
- 230000009194 climbing Effects 0.000 description 1
- 239000000567 combustion gas Substances 0.000 description 1
- 238000010835 comparative analysis Methods 0.000 description 1
- 238000005094 computer simulation Methods 0.000 description 1
- 230000010485 coping Effects 0.000 description 1
- 230000007423 decrease Effects 0.000 description 1
- 230000007812 deficiency Effects 0.000 description 1
- 230000001934 delay Effects 0.000 description 1
- 210000004262 dental pulp cavity Anatomy 0.000 description 1
- 230000001419 dependent effect Effects 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 238000004146 energy storage Methods 0.000 description 1
- 230000007613 environmental effect Effects 0.000 description 1
- 239000012530 fluid Substances 0.000 description 1
- 238000010348 incorporation Methods 0.000 description 1
- 238000009434 installation Methods 0.000 description 1
- 238000007726 management method Methods 0.000 description 1
- 238000002156 mixing Methods 0.000 description 1
- 230000000750 progressive effect Effects 0.000 description 1
- 230000009467 reduction Effects 0.000 description 1
- 239000007787 solid Substances 0.000 description 1
- 230000003068 static effect Effects 0.000 description 1
- 239000008400 supply water Substances 0.000 description 1
- 230000033772 system development Effects 0.000 description 1
- 230000009897 systematic effect Effects 0.000 description 1
- 230000007704 transition 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/0631—Resource planning, allocation, distributing or scheduling for enterprises or organisations
-
- 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
Landscapes
- Business, Economics & Management (AREA)
- Human Resources & Organizations (AREA)
- Engineering & Computer Science (AREA)
- Economics (AREA)
- Strategic Management (AREA)
- Tourism & Hospitality (AREA)
- Physics & Mathematics (AREA)
- Health & Medical Sciences (AREA)
- Theoretical Computer Science (AREA)
- Marketing (AREA)
- Entrepreneurship & Innovation (AREA)
- General Physics & Mathematics (AREA)
- General Business, Economics & Management (AREA)
- Operations Research (AREA)
- Development Economics (AREA)
- Quality & Reliability (AREA)
- Educational Administration (AREA)
- Game Theory and Decision Science (AREA)
- Public Health (AREA)
- Water Supply & Treatment (AREA)
- General Health & Medical Sciences (AREA)
- Primary Health Care (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
- Feedback Control In General (AREA)
Abstract
The invention discloses a kind of integrated energy system dispatching methods of dynamic interval, include the following steps, establish the integrated energy system Scheduling Framework based on dynamic interval Model Predictive Control；Establish integrated energy system scheduling model a few days ago；Establish in a few days scheduling " trajector deviation-control energy-control cost " performance indicator based on Model Predictive Control；Establish dynamic interval decision index system；Establish energy conversion unit, heat supply network model, gas pessimistic concurrency control, electric network model.Beneficial effects of the present invention: consider that being conducive to performance indicator has enough adjustable nargin, copes with different scene demands with the integrated energy system dispatching method for the dynamic interval that " trajector deviation-control energy-control cost " is performance indicator；Dynamic interval decision index system is considered, the time of day for even having occurred and that large error will be occurred by being conducive to reflection system, provide decision index system for whether system needs to assign dispatch command.
Description
Technical field
The present invention relates to the technical fields of integrated energy system scheduling, more particularly to one kind to be based on dynamic interval model
The integrated energy system dispatching method of dynamic interval.
Background technique
With a large amount of accesses of renewable energy, integrated energy system is the important channel realizing high efficiency of energy and utilizing.It is comprehensive
The Optimized Operation for closing energy resource system is the key technology realizing all kinds of energy harmonizings and utilizing.In Optimized Operation, dispatching party
The accuracy of case is heavily dependent on the levels of precision of internal loading dispatching cycle prediction, and net load prediction inevitably exists certain
Error, and increase with the increase of dispatching cycle, predict that influence of the error to Optimized Operation can be by dispatching cycle
It divides to reduce.However the response speed of electric power networks, natural gas network and the instruction of heating power network degree of exchanging is different, electric power
System inertia is minimum, and adjustment speed is fast, and natural gas system takes second place, and hot (cold) system is most slow, since comprehensive energy subsystem is dynamic
There are significant differences for step response, workload demand/response characteristic etc., when adding same on discontinuity surface, different sub-systems
Schedulable resource is also different, these factors bring challenge to the selection of integrated energy system dispatching cycle, if dispatching cycle
It is too long, lead to that there are large errors with the actual operational process of system, if dispatching cycle is too small, meaningless increase system development and
Operating cost.Therefore, the rational management time how is selected, comprehensive energy is realized by scheduling in necessary and correct time section
The resources effective utilization of system is current the problem of being badly in need of research.
It is at this stage based on Fixed Time Interval, the day of optimization system mostly about the scheduling research of integrated energy system
Before/in a few days scheduling strategy be with Fixed Time Interval and to consider rail mostly in terms of the dispatching method based on Model Predictive Control
The performance indicator of mark deviation, but these researchs seldom consider the dispatching method of the Model Predictive Control of dynamic interval.
Summary of the invention
The purpose of this section is to summarize some aspects of the embodiment of the present invention and briefly introduce some preferable implementations
Example.It may do a little simplified or be omitted to avoid our department is made in this section and the description of the application and the title of the invention
Point, the purpose of abstract of description and denomination of invention it is fuzzy, and this simplification or omit and cannot be used for limiting the scope of the invention.
In view of above-mentioned existing problem, the present invention is proposed.
Therefore, it is an object of the present invention to provide a kind of integrated energy system dispatching methods of dynamic interval.
In order to solve the above technical problems, the invention provides the following technical scheme: a kind of comprehensive energy of dynamic interval
System scheduling method includes the following steps, establishes the integrated energy system scheduling based on dynamic interval Model Predictive Control
Frame obtains the scheduling principle of the integrated energy system based on dynamic interval Model Predictive Control；Establish comprehensive energy system
System scheduling model a few days ago, obtain a few days ago operation plan as the reference locus in a few days dispatched；It establishes based on Model Predictive Control
In a few days dispatch " trajector deviation-control energy-control cost " performance indicator, the objective function as in a few days optimizing scheduling；It establishes
Dynamic interval decision index system, judges whether scheduling interval needs to change；Reference locus correction index is established, is judged with reference to rail
Whether mark, which needs, corrects；Energy conversion unit, heat supply network model, gas pessimistic concurrency control, electric network model are established, integrated energy system tune is obtained
Degree constraint condition obtains the in a few days scheduling scheme of integrated energy system in the case where meeting constraint condition.
A kind of preferred embodiment of integrated energy system dispatching method as dynamic interval of the present invention,
In: the integrated energy system Scheduling Framework of the foundation based on dynamic interval Model Predictive Control, when obtaining based on dynamic
Between gap model PREDICTIVE CONTROL integrated energy system scheduling principle it is further comprising the steps of, dispatched a few days ago to run minimized
Cost is target making hour service capacity plan in lower day and issues.And in a few days rolling optimal dispatching is then in terms of to dispatch a few days ago
It divides reference locus w (k+i) into, carries out rolling optimization by initial value of running power output.Shape is run according to the system at k moment
State sequence u (k), the system mode sequence based on the prediction model prediction k+1 momentSystem mode sequenceWith
The system mode error sequence e (k) at k moment collectively forms output prediction linkOn the one hand output prediction link passes through
Cross dynamic interval decision index system ξt, judge whether the instruction for needing to assign amendment operation plan, on the other hand by a few days
Reference locus corrects indexJudge whether to need to be modified reference locus w (k+i).
A kind of preferred embodiment of integrated energy system dispatching method as dynamic interval of the present invention,
In: it is described to establish integrated energy system scheduling model a few days ago, obtain a few days ago operation plan as the reference locus in a few days dispatched also
Include the following steps, dispatching a few days ago and optimizing scheduling by target of economic optimum, objective function is as follows:
In formula: T is the Optimized Operation period a few days ago, is taken for 24 hours, time interval 1h；U is generating set set；V is natural gas
Source set；Eu,tFor the active power output of t moment generating set u, Gv,tFor the gas production of t moment gas source v；For generating set u's
Cost of electricity-generating,For the Gas Prices of gas source v.
A kind of preferred embodiment of integrated energy system dispatching method as dynamic interval of the present invention,
In: described in a few days scheduling " trajector deviation-control energy-control cost " performance indicator of foundation based on Model Predictive Control is made
Objective function in a few days optimizing scheduling is further comprising the steps of, the reference locus w (k+i) dispatched a few days ago and output predictionAfter scene Recognition, optimized according to performance indicator, performance indicator is the objective function of rolling optimization.Performance
Target function considers trajector deviation, control energy, control cost, as shown in formula:
Being expressed as matrix form is formula:
Wherein, T is prediction time domain, and M is control time domain, Q=diag [q1, q2..., qT], R=diag [r1, r2..., rM]、
S=diag [s1, s2..., sM] be different scenes under weight matrix,To predict that output matrix, W (k+1) are reference
Track matrix, Δ p (k) are controlling increment matrix, and C (k) is control Cost matrix.Predict output matrix and reference locus matrix packet
The active power output and gas production of equipment are included, controlling increment matrix is regulated quantity of the equipment in plan a few days ago, controls Cost matrix
The cost of load is supplied for equipment.
A kind of preferred embodiment of integrated energy system dispatching method as dynamic interval of the present invention,
In: it is described to establish dynamic interval decision index system, judge scheduling interval whether need to change it is further comprising the steps of, based on to
Decision-making period predict time domain in integrated energy system overall situation operating cost in a few days under reference locus the overall situation operating cost deviation
Rate establishes dynamic interval decision index system ξt.Global cost is shown below:
S in formular,tIndicate that decision period and the global cost for predicting time domain, A indicate the unit of 1 row T column to t moment
Battle array, P indicate that decision period and prediction time domain regulate and control moment matrix, p in t momentlT,tIndicate t moment Demand-side and system side resource
Regulation amount of the l in prediction time domain T, C expression Cost matrix, clIndicate the cost coefficient regulated and controled to l.
Dynamic interval decision index system ξtShown in being expressed as follows:
Wherein, Sw,tIndicate the global cost of reference locus.
When the deviation between the actual path and reference locus after the optimizing scheduling based on Model Predictive Control is to complete
Office's operating cost will trigger dynamic interval decision when producing bigger effect, change control room every, in advance to control time domain into
Row, which optimizes, simultaneously executes dispatch command, when influencing smaller, does not then trigger decision, and control room reduces unnecessary holds every will be late by
Row instruction.
A kind of preferred embodiment of integrated energy system dispatching method as dynamic interval of the present invention,
In: it is described to establish reference locus correction index, judge that whether reference locus needs to correct further comprising the steps of, works as actual motion
When the deviation that scene and prediction scene are likely to occur, the reference significance for being intended to be reference locus a few days ago is just reduced.Therefore will
Actual path is compared with reference locus, establishes reference locus correction index.Operation plan essentially dictates unit a few days ago
Start and stop and combination plan, if therefore current time occur relatively large deviation, Unit Commitment or combination plan should be adjusted in time, then
In a few days reference locus corrects indexIt can indicate are as follows:
Wherein,For actual path,For in a few days reference locus.
A kind of preferred embodiment of integrated energy system dispatching method as dynamic interval of the present invention,
In: it is described to establish energy conversion unit, heat supply network model, gas pessimistic concurrency control, electric network model, obtain integrated energy system schedule constraints item
Part, in the case where meeting constraint condition, the in a few days scheduling scheme for obtaining integrated energy system is further comprising the steps of, the energy
It measures converting unit and considers energy balance constraint；The heat supply network model considers that node flow balance, node power fusion, load are taken
Characteristic, supply and return water temperature constraint and pipeline section heat-transfer character；The gas pessimistic concurrency control consideration pipeline flow constrains, gas source point constrains,
Flow equilibrium constraint, compressor constraint and node pressure constraint；The electric network model considers node power balance, unit output
Constraint, Climing constant and Branch Power Flow constraint.
Beneficial effects of the present invention: when considering the dynamic with " trajector deviation-control energy-control cost " for performance indicator
Between the integrated energy system dispatching method that is spaced, being conducive to performance indicator has enough adjustable nargin, and coping with different scenes needs
It asks；Dynamic interval decision index system is considered, being conducive to reflection system will occur even to have occurred and that the true of large error
Real state provides decision index system for whether system needs to assign dispatch command；In a few days reference locus correction index is considered, favorably
In the deviation of reply actual scene and prediction scene, real now forecast time domain global optimization scheduling facilitates the popularization of this programme
With implementation.
Detailed description of the invention
In order to illustrate the technical solution of the embodiments of the present invention more clearly, required use in being described below to embodiment
Attached drawing be briefly described, it should be apparent that, drawings in the following description are only some embodiments of the invention, for this
For the those of ordinary skill of field, without any creative labor, it can also be obtained according to these attached drawings other
Attached drawing.Wherein:
Fig. 1 is the bulk flow of the integrated energy system dispatching method of dynamic interval described in first embodiment of the invention
Cheng Tu；
Fig. 2 is the integrated energy system based on dynamic interval Model Predictive Control described in first practical example of the invention
Scheduling Framework figure；
Fig. 3 is the test example figure of integrated energy system described in first embodiment of the invention；
Fig. 4 is that the integrated energy system dispatching method of dynamic interval described in first embodiment of the invention is calculated in test
Scheduling scheme a few days ago in example；
Fig. 5 is that the integrated energy system dispatching method of dynamic interval described in first embodiment of the invention is calculated in test
In a few days scheduling scheme in example；
Fig. 6 is the integrated energy system dispatching method of dynamic interval described in second embodiment of the invention and is based on solid
It fixes time gap model comparing result schematic diagram；
Fig. 7 is that performance indicator described in second embodiment of the invention is shown with the result of variations of dynamic interval decision index system
It is intended to；
Fig. 8 be second embodiment of the invention described in dynamic interval integrated energy system dispatching method with do not consider
The dynamic interval model comparing result schematic diagram of reference locus correction index；
Fig. 9 is that performance indicator described in second embodiment of the invention is illustrated with the result of variations that reference locus corrects index
Figure.
Specific embodiment
In order to make the foregoing objectives, features and advantages of the present invention clearer and more comprehensible, right with reference to the accompanying drawings of the specification
A specific embodiment of the invention is described in detail, it is clear that and described embodiment is a part of the embodiments of the present invention, and
It is not all of embodiment.Based on the embodiments of the present invention, ordinary people in the field is without making creative work
Every other embodiment obtained, all should belong to the range of protection of the invention.
In the following description, numerous specific details are set forth in order to facilitate a full understanding of the present invention, but the present invention can be with
Implemented using other than the one described here other way, those skilled in the art can be without prejudice to intension of the present invention
In the case of do similar popularization, therefore the present invention is not limited by the specific embodiments disclosed below.
Secondly, " one embodiment " or " embodiment " referred to herein, which refers to, may be included at least one realization side of the invention
A particular feature, structure, or characteristic in formula." in one embodiment " that different places occur in the present specification not refers both to
The same embodiment, nor the individual or selective embodiment mutually exclusive with other embodiments.
Combination schematic diagram of the present invention is described in detail, when describing the embodiments of the present invention, for purposes of illustration only, indicating device
The sectional view of structure can disobey general proportion and make partial enlargement, and the schematic diagram is example, should not limit this herein
Invent the range of protection.In addition, the three-dimensional space of length, width and depth should be included in actual fabrication.
Simultaneously in the description of the present invention, it should be noted that the orientation of the instructions such as " upper and lower, inner and outer " in term
Or positional relationship is to be based on the orientation or positional relationship shown in the drawings, and is merely for convenience of description of the present invention and simplification of the description, and
It is not that the device of indication or suggestion meaning or element must have a particular orientation, be constructed and operated in a specific orientation, therefore
It is not considered as limiting the invention.In addition, term " first, second or third " is used for description purposes only, and cannot understand
For indication or suggestion relative importance.
In the present invention unless otherwise clearly defined and limited, term " installation is connected, connection " shall be understood in a broad sense, example
Such as: may be a fixed connection, be detachably connected or integral type connection；It equally can be mechanical connection, be electrically connected or be directly connected to,
Can also indirectly connected through an intermediary, the connection being also possible to inside two elements.For the ordinary skill people of this field
For member, the concrete meaning of above-mentioned term in the present invention can be understood with concrete condition.
Embodiment 1
In terms of integrated energy system Optimized Operation, current research can establish setting for separate unit and coupling unit
On the basis of standby physical model, economic model and electric power, natural gas, heating power energy subsystem steady-state load flow model, needle is realized
To the comprehensive energy Optimized Operation of the particular systems such as electric-gas, electric-thermal, also it is capable of forming excellent including energy transition equipment
Change and adjust strategy, robust cooperative scheduling optimization can also be carried out for the uncertain problem of the renewable energy such as wind-powered electricity generation, also can
The responding ability of user is excavated by pricing information and formulates a few days ago/in a few days operation plan, but these Optimization Schedulings are most
It is to optimize scheduling controlling for sometime section or more periods, still falls within static optimization.For electric system, naturally
Gas system and therrmodynamic system also have correlative study can be to electricity-in the fast dynamics and slow motion step response of different time scales
The particular systems such as gas, electric-thermal carry out in a few days/ultra-short term rolling optimal dispatching, but still are open loop Optimization Scheduling, by right
Load prediction precision can be improved in the division of time scale, but current research is not counted and real system operation is to optimal control
The influence of process.It is with one day for 24 hours for scheduling time domain, 15min is most small step but the present invention is based on Model Predictive Control principle
Long, the optimal value of current time variable has an impact the variate-value of following instant, belongs to dynamic optimization, setting output prediction loop
Section is fed back, and the optimization of moment variable considered the influence to following instant variable at that time, belonged to closed-loop optimization.Last is defeated
Out the result is that one day for 24 hours in multiple moment dispatch value.
Influence of the real system operation to optimal control process is that " it is pre- that load can be improved by the division to time scale
Survey precision ", time scale can be divided into the several seconds, but so will not frequently be controlled when real system operation, can increase
Add very big scheduling cost.By temperature when actual motion, many factors such as pressure influence, and there are certain loss, effective supplies
There are certain deviations for value and operation plan value.
The deficiency on anti-interference ability and robustness, Model Predictive Control (Model are dispatched for open loop
Predictive Control, MPC) method is gradually applied on Optimized Operation.MPC is as a kind of system optimization controlling party
Method, it is different from refinement time scale Optimization Scheduling, quantity of state feedback compensation link is introduced, prediction error etc. can be corrected
Optimized Operation deviation caused by factor obtains the optimum control variable of current and future period, so that following export and refer to rail
Mark deviation is minimum, has extremely strong anti-interference ability and robustness.MPC method has applied to family lan, micro-capacitance sensor, has matched
All kinds of grid load frequency controls such as power grid and supply of cooling, heating and electrical powers type microgrid, cool and thermal power association system and electric power, comprehensive energy
The Optimized Operation in field, and good effect is obtained in terms of stability contorting and system robust, but mostly between the set time
It is divided into and rolls the period and using trajector deviation as performance indicator rolling optimization is successively carried out to current time, it is difficult to accurate description load
The influence to the flowing of integrated energy system self-energy such as prediction, renewable energy power output, in the scheduling track of multistage time scale
Amendment, the scene correction of different time section, control room study deep not enough, less consideration every adjustment etc..This implementation
Example can only describe the system mode at given at the time of point by Fixed Time Interval, that is, can only carry out in specific point
Scheduling, and these points not necessarily (assigning dispatch command or not assigning dispatch command and scheduling is tied there is an urgent need to scheduling
The influence of fruit is little).Dispatch command, shadow are exactly assigned in the scheduling of dynamic interval when system is there is an urgent need to dispatch
Ringing little point, there is no need to reschedule.
Signal referring to Fig.1 illustrates that is proposed in the present embodiment considers the integrated energy system multiclass of system reserve value
The overall flow figure of type energy storage configuration method, the also as comprehensive energy of dynamic interval described in first embodiment of the invention
The overall flow figure of system scheduling method is made of in specific implementation procedure following steps:
S100: the integrated energy system Scheduling Framework based on dynamic interval Model Predictive Control is established, is referring to Fig. 2
Integrated energy system Scheduling Framework figure described in the present embodiment based on dynamic interval Model Predictive Control is obtained based on dynamic
The scheduling principle of the integrated energy system of time interval Model Predictive Control.Specifically, it is further comprising the steps of in this step,
It is dispatched a few days ago using the cost that runs minimized as target making hour service capacity plan in lower day and is issued.And in a few days
Rolling optimal dispatching is then to be rolled using operation plan a few days ago as reference locus w (k+i) using running power output as initial value
Dynamic optimization.
The scheduling principle for being namely based on Model Predictive Control of scheduling principle description.This method of Model Predictive Control itself
It is not intended to do and dispatch, be a kind of control theory.Had one in electric power system dispatching, the scheduling of supply of cooling, heating and electrical powers type microgrid later
A little applications, then we also apply it in the scheduling of integrated energy system.
Dispatching method of the invention is based on " Model Predictive Control " first, and feature or innovative point are " dynamic times
Interval ".Scheduling principle is exactly the principle that Model Predictive Control is applied to integrated energy system scheduling, is with operation plan a few days ago
Reference locus.System in future state is predicted by prediction model, and collectively forms output in advance with the state error of system
Link is surveyed, output prediction link enters dynamic interval decision index system, judges whether to need to correct scheduling interval, then carry out
The optimization of performance indicator；Another aspect output prediction link also enters in a few days reference locus and corrects index, judges whether needs pair
Reference locus is modified.
It is specific as follows: according to the system running state sequence u (k) at k moment, to be based on what prediction model predicted the k+1 moment
System status switch, system mode sequence and the system mode error sequence e (k) at k moment collectively form output prediction link, output
On the one hand prediction link passes through dynamic interval decision index system ξt, judge whether the instruction for needing to assign amendment operation plan,
This judgement herein is to need to provide a judgment criteria according to system requirements, and the result that standard difference then optimizes is different, quite
Then an input parameter.I can list 5 groups of dynamic interval decision index system ξ t of setting in example below, and respectively 0.04,
0.08,0.12,0.16,0.20, scheduling result is compared, it is the parameter of an input that judgment criteria is provided according to the demand of system,
Different reference locus correction indexs is considered in example, and 5 groups of reference locus are set and correct indexRespectively 0.04,0.08,
0.12,0.16,0.20, it is scheduled the comparative analysis of result.On the other hand index is corrected by a few days reference locusJudgement
Whether need to be modified reference locus w (k+i).
The instruction for correcting operation plan is according to the calculation formula of index ξ t, this index at computing system current time
Value is then adjusted on the basis of needing in the operation plan that reference locus provides at current time when reaching given standard
Degree adjustment, assigns new dispatch command, if not up to given standard, current time according to script reference locus scheduling
Plan, does not need additionally to assign dispatch command.
Track correct is according to indexCalculation formula, this index value at computing system current time, when reach to
It when fixed standard, then needs to be modified the operation plan that reference locus provides at current time, changes reference locus, specifically
Operation is the scheduling result at the moment before, re-starts optimizing scheduling a few days ago, updates current time and tune a few days ago later
Degree plan, that is, reference locus.If not up to given standard, current time still use script reference locus.
Difference with dynamic interval decision index system is:
Reference locus corrects index, and correction is reference locus itself, and dynamic interval decision index system is modified to be
No needs are in a few days dispatched at current time according to reference locus.The Model Predictive Control mainly includes three parts prediction mould
Type, feedback compensation and rolling optimization.Prediction model is exactly the result according to scheduled good period prediction following sessions.
It should also be noted that, prediction result is a sequence, it is not a value.For example, there are within one day 96 15min,
8 points of (k+1 moment) this moment are scheduled, just 8 points are once predicted at 7: 45 (k moment).Based on 8 points it
Preceding 32 data (this 32 data are k moment, that is, 7: 45 already existing operation datas) predict 64 numbers below
According to this sequence (the system mode sequence at k+1 moment) of output includes 96 that 32 data and 64 data form together
Data.This is to predict the sequence come.And system mode error sequence only has 7: 45 and its data before, that is, 32
Error information, this 32 error informations constitute the state error sequence at 7: 45 (k moment).That is, 7: 45 to not
Come at 8 points to be predicted, 8 points of prediction result is original existing 32 data+subsequent 96 data, wherein containing at 8 points
Predicted value.But there was only 32 data in 7: 45 error information, there is no 8 points of error informations.So state error sequence
Column are the k moment, and the system mode sequence of prediction is the k+1 moment.It is by 96 data when then carrying out performance index optimization
Optimize together (7: 45 and scheduling numerical value before will not change, mainly to the scheduling numerical value after 8 points and 8 points
Optimization), guarantee global optimum, but only carry out the optimizing scheduling at 8 o'clock as a result, although later point is optimized, but
It is not execute, the waiting time is rolled to following time interval.It to be scheduled to 8: 15,8: 15 become k+1, and 8 points become
K is once predicted again 8: 15, that is, based on 33 data before, predict subsequent 63 data, behind with 8 points
Scheduling it is identical.
S200: establishing integrated energy system scheduling model a few days ago, obtain a few days ago operation plan as the reference in a few days dispatched
Track.It is specific further comprising the steps of,
It dispatches a few days ago and optimizes scheduling by target of economic optimum, objective function is as follows:
In formula: T is the Optimized Operation period a few days ago, is taken for 24 hours, time interval 1h；U is generating set set；V is natural gas
Source set；Eu,tFor the active power output of t moment generating set u, Gv,tFor the gas production of t moment gas source v；For generating set u's
Cost of electricity-generating,For the Gas Prices of gas source v.Scheduling model a few days ago is solved, tune hourly under obtaining 1 day for 24 hours
Degree plan, that is, the power output of each equipment per hour.Operation plan is equivalent to be the rank as unit of hour a few days ago for this
Jump signal.Then in a few days scheduling is exactly using this step signal dispatched a few days ago as reference locus, in the base of this step signal
It is adjusted on plinth.
S300: it establishes in a few days scheduling " trajector deviation-control energy-control cost " performance based on Model Predictive Control and refers to
Mark, the objective function as in a few days optimizing scheduling.It is specific further comprising the steps of,
The reference locus w (k+i) dispatched a few days ago and output predictionReference locus is adjusted according to optimization a few days ago
The result of degree obtains, and output prediction is obtained by system mode sequence and system mode error sequence, after scene Recognition, foundation
Performance indicator optimizes, which refers to the minimum of performance indicator, and min F can embody in formula, and performance indicator is just
It is F, comprising trajector deviation, control energy and control cost three parts, performance indicator is the objective function of rolling optimization.Performance refers to
Scalar functions consider trajector deviation, control energy and control cost, as shown in formula:
Being expressed as matrix form is formula:
Wherein, T is prediction time domain, and M is control time domain, Q=diag [q1, q2..., qT], R=diag [r1, r2..., rM]、
S=diag [s1, s2..., sM] be different scenes under weight matrix,To predict that output matrix, W (k+1) are reference
Track matrix, Δ p (k) are controlling increment matrix, and C (k) is control Cost matrix.Predict output matrix and reference locus matrix packet
The active power output and gas production of equipment are included, controlling increment matrix is regulated quantity of the equipment in plan a few days ago, controls Cost matrix
The cost of load is supplied for equipment.
S400: establishing dynamic interval decision index system, judges whether scheduling interval needs to change；It further include following step
Suddenly,
Based on to integrated energy system in decision period forecasting time domain, (" be based on dynamic interval Model Predictive Control " be
The method of scheduling, " to decision period forecasting time domain " are the global operating cost under this method, when carrying out in a few days rolling optimization
An expression) global operating cost and operating cost global under in a few days reference locus deviation ratio, establish between dynamic time
Every decision index system ξt.Global cost (the global cost of t moment decision period and prediction time domain, i.e. Sr, t) is shown below:
S in formular,tIndicate that decision period and the global cost for predicting time domain, A indicate the unit of 1 row T column to t moment
Battle array, P indicate that decision period and prediction time domain regulate and control moment matrix, p in t momentlT,tIndicate t moment Demand-side and system side resource
Regulation amount of the l in prediction time domain T, C expression Cost matrix, clIndicate the cost coefficient regulated and controled to l.
Dynamic interval decision index system ξtShown in being expressed as follows:
Wherein, Sw,tIndicate the global cost of reference locus.
When the deviation between the actual path and reference locus after the optimizing scheduling based on Model Predictive Control is to complete
Office's operating cost will trigger dynamic interval decision when producing bigger effect, change control room every, in advance to control time domain into
Row, which optimizes, simultaneously executes dispatch command, when influencing smaller, does not then trigger decision, and control room reduces unnecessary holds every will be late by
Row instruction.
S500: reference locus correction index is established, judges whether reference locus needs to correct.It is further comprising the steps of,
When the deviation that actual motion scene and prediction scene are likely to occur, it is intended to be the reference meaning of reference locus a few days ago
Justice just reduces.Therefore actual path and reference locus are compared, establishes reference locus correction index.
Operation plan essentially dictates the start and stop and combination plan of unit a few days ago, if therefore occurring at current time larger inclined
Difference should then adjust Unit Commitment or combination plan in time, then in a few days reference locus corrects indexIt can indicate are as follows:
Wherein,For actual path,For in a few days reference locus.
S600: establishing energy conversion unit, heat supply network model, gas pessimistic concurrency control and electric network model, obtains integrated energy system tune
Degree constraint condition obtains the in a few days scheduling scheme of integrated energy system in the case where meeting constraint condition.Above-mentioned establishes three
A model, these models are mainly physical model, that is, need to meet energy conversion unit, heat supply network, gas net, power grid object
Reason constraint, when being then scheduled to integrated energy system, needs to meet these constraint conditions.The constraint condition refers to energy
Converting unit, heat supply network, gas net and power grid need the physical constraint met, and the boundary that different systems needs to meet can be different,
Which boundary I can indicate in example below.
In a few days scheduling scheme refer mainly to all fired power generating units power output and natural air-air source air demand (heat be also pass through
Electricity, natural gas are by heating equipment supply, and the energy input of most original is exactly that electricity is gentle, so final scheduling scheme sums up
For the power output of fired power generating unit and the air demand of natural air-air source).
It is further more specifically, wherein further comprising the steps of in S600,
S601: energy conversion model of element is established.
The energy conversion unit of integrated energy system passes through Coupling device for electric system, natural gas system and heating power system
System connects, and the energy conversion unit of the present embodiment building mainly includes cogeneration of heat and power (CHP), power transformer (T), heat exchange
Device (HE), heat regenerator (HR), electricity turn gas (P2G), electricity refrigeration (EC) and absorption refrigeration (AC) equipment, energy-balance equation
(energy conversion model of element) are as follows:
Wherein Le、Lg、Lh、LeRespectively electricity, air and heat, cold load consumption, ηCHP、ηHE、ηT、ηHR、ηAC、ηP2G、ηECFor
The efficiency of CHP, HE, T, HR, AC, P2G, EC, φCHPFor the hotspot stress of CHP, λe,1、λe,2、λe,3For the electric power distribution ratio of input
Example, λg,1、λg,2For the natural gas distribution ratio of input, PgFor gas consumption, PeFor electric power consumption.
S602: establishing the heat supply network model of integrated energy system, herein heat supply network model be not a simple expression formula or
Figure, but need to meet certain physical constraint in pipeline, the hot water flow of node, temperature, it is following to specifically describe:
Heating system often uses steam and hot water as heat-carrying agent.Regard heat supply network as fluid network, mainly considers node-flow
Amount balance, node power fusion, load take characteristic, supply and return water temperature constraint and pipeline section heat-transfer character.Specifically, as follows:
1. node flow balances, for any node in heat supply network, the sum of hot water flow of inflow be equal to outflow flow it
With that is,In formula,WithRespectively it is connected with node i and from node i starting and ending pipeline
Set；For hot water quality's flow in period t pipeline j.
2. node temperature merges, the hot water of different temperatures is mixed after different pipeline flow-direction same node points, after mixing
The hot water temperature for flowing into different pipelines from same node is identical, i.e.,In formula,For period t pipe
Hot water outlet temperature in road j；For hot water inlet's temperature in period t pipeline k.
3. load takes characteristic, for the heat supply network branch comprising heat user, load bus i is in period t consumption of calorie
Flowing through load bus mass flow isWater by supply water temperatureIt is down to return water temperatureI.e.
4. supply and return water temperature constrains, in order to guarantee the heating quality of heat source and heat user, need to heat source and heat user
It is limited for, return water temperature, i.e.,
5. pipeline section heat-transfer character, the flowing that heat supply network relies on hot water realizes that energy transmission, the sluggishness of delivery may be led
Several minutes of power transmission delays not waited to a few hours are caused, this dispatches the ultra-short term of integrated energy system, and there may be significant
It influences.Steady state heat transfer characteristic is expressed as:In formula, x be on pipeline section certain point and pipeline section head end away from
From；R is the thermal resistance of pipeline section unit length, and Ts, Te, Ta is respectively the head end temperature of a root canal section, temperature and ambient temperature at x, f
For hot water flow.Transient state heat-transfer character is as follows: for apart from the closer place of heat source, transient process is shorter, Primary regulation it
Afterwards, before adjusting next time, pipeline section temperature has reached stable state, and the transient state heat-transfer character of pipeline section can indicate at these points are as follows:
In formula, Ti(x, t) is temperature of the heat-net-pipeline in the i-th period, at heat source x in t moment；WithPoint
It Wei not ti-1Moment and ti-2The temperature of moment heat source；I=1,2,3 ....
For apart from the farther away place of heat source, transient process is longer, and after Primary regulation, stable state is had not yet been reached in temperature, under
Primary regulation has begun, and the transient state heat-transfer character of pipeline section can indicate at these points are as follows:
In formula, Ti(x, t) is temperature of the heat-net-pipeline in the (i-1)-th period, at heat source x at the ti-1 moment；I=1,
2,3 ....
S603: establishing the gas pessimistic concurrency control of integrated energy system, be on natural gas network, it is natural in pipeline, node
The physical constraint that throughput, pressure need to meet, described below:
Gas pessimistic concurrency control mainly includes pipeline flow constraint, gas source point constraint, flow equilibrium constraint, compressor constraint and section
Press force constraint.
1. pipeline flow constrains, natural gas line flow equation is related with pipe ends pressure and many physical characteristics of pipeline,
General form is had no, the gas flow under particular condition is usually described with nonlinear equation.Gas pipeline is insulated for ideal, is examined
Consider natural gas two-way flow, flow equation may be expressed as:In formula,
Indicate that t moment flows through the average flow rate of pipeline ij, whereinThe respectively first section natural gas filling of t moment pipeline ij becomes a mandarin
Amount and end natural gas output flow；Cij is the related constants such as pipeline ij efficiency, temperature, length, internal diameter, compressibility factor；
pi,t、pj,tRespectively t moment first and last node i, the pressure value of j.
2. gas source point constrains,Wherein,The natural gas of respectively gas source point n supplies
Answer flow bound.
3. flow equilibrium constrains,Wherein,For t moment node
Gas source feed flow on i；Turn gas gas supply flow for the electricity in t moment node i；For combustion gas wheel in t moment node i
The gas discharge of machine consumption；For the natural gas load in t moment node i；It is consumed for CCHP in t moment node i
Gas discharge.
4. compressor constrains, using simplified compressor model are as follows: pl,t≤βcompi,t, in formula, pl,tFor the pressure of compressor
Contracting coefficient.
5. node pressure constrains,Wherein,Respectively node i pressure value is upper and lower
Limit.
S604: establish integrated energy system electric network model (at the route, node of power grid, power, voltage, phase angle need
Meet certain physical constraint).
Power system modeling mainly include node power balance, unit output constraint, Climing constant and Branch Power Flow about
Beam.
1. node power balances,Wherein,For
The active power output of fired power generating unit in t moment node i；For the active power output of Wind turbines in t moment node i；When for t
Carve the active power output of gas turbine in node i；For the burden with power in t moment node i；Pij,tFor on t moment route ij
Active power；For the idle power output of fired power generating unit in t moment node i；For the load or burden without work in t moment node i；Qij,t
For the reactive power on t moment route ij.
2. unit output constrains,Wherein,Fired power generating unit is contributed respectively in node i
Bound.
3. Climing constant,Wherein, RUi、RDiFor fired power generating unit in node i
Climbing and descending grade amplitude.
4. Branch Power Flow constrains,
Wherein, Vi,tFor the voltage magnitude in t moment node i；Voltage magnitude is upper and lower respectively in node i
Limit；θijFor t moment node i, the phase difference of voltage of j；The respectively bound of t moment node i j phase difference of voltage；
GijFor the conductance between node i j；BijFor the susceptance between node i j；Route is active between respectively node i j
The bound of power transmission.
Embodiment 2
Referring to the signal of Fig. 3, it is illustrated as based on improved IEEE14 node power distribution net system, 11 node natural gas distribution systems
It unites, the example test macro for the integrated energy system that the therrmodynamic system of 11 nodes is built, i.e. test example figure.By gas discharge
It is converted into power unit according to heating value of natural gas, operation plan is illustrated as the comprehensive of dynamic interval referring to the signal of Fig. 4 a few days ago
Close a few days ago scheduling scheme of the energy resource system dispatching method in test example.This few days ago scheduling scheme this be also in a few days operation plan
Reference locus, in a few days scheduling scheme is obtained based on dispatching method, referring to the signal of Fig. 5, takes the day of 18:00-20:45 period
Interior scheduling result is as displaying.
System is when state error is lesser, it is possible to reduce dispatch command issues, reduce scheduling cost, error in
When plan difference is larger a few days ago, then need to be modified the reference locus planned a few days ago.
Specifically, the in a few days scheduling of script is the fixed intervals of 15min, in 19:30-20 referring again to the signal of Fig. 5:
00 this 30min is a scheduling interval, has ignored the scheduling interval of script 19:45.This is because load fluctuation and ring at this time
Border variation is smaller, and the error generated to system is smaller.In the electricity power output reference locus and 20:00-20:45 of 18:00-19:00
Natural gas power output track corrected, this is because the load of the two periods and environmental factor and the feelings predicted a few days ago
Relatively large deviation has occurred in condition, is affected to system cost.Final scheduling cost is
Dispatching method of the invention is compared with based on Fixed Time Interval model, that is, does not consider dynamic interval
The scene of decision index system.As a result excellent although the system mode in this 30min of 19:30-20:00 is highly stable referring to the signal of Fig. 6
The deviation very little of power and reference locus is dissolved, system is still to carry out successively dispatch command every 15 minutes, and dispatch each time
Instruction has corresponding energy variation, causes to control energy and controls increased costs, final scheduling cost is
Compared with dispatching method of the invention, scheduling cost is increased.
Influence of the different dynamic interval decision index systems to scheduling result is compared, considers different dynamic intervals
5 groups of dynamic interval decision index system ξ are arranged in decision index systemt, the value of triggering is respectively 0.04,0.08,0.12,0.16,
0.20.If the dynamic interval decision index system value of current state is less than the value, dynamic interval decision is not triggered, such as
Fruit is greater than, then triggers.Performance indicator with dynamic interval decision index system result of variations referring to Fig. 7 signal.Dynamic time
The different numerical value for being spaced decision index system can generate different influences to the optimization of system performance index.When dynamic interval decision
It is unobvious to the variation detection of system mode error when index is excessive, cause trajector deviation to increase, and then performance indicator increases.
It is higher to the precision of system mode error when dynamic interval decision index system is too small, cause scheduling excessively frequent, although having
Conducive to reduction trajector deviation, but control energy and control cost are increased, and then performance indicator increases.According to Fig. 7, with dynamic
The increase of time interval decision index system, performance indicator first reduce, then increase.First reduce is because of dynamic interval decision index system
Amplification, reduce control energy and control cost.Increase again is because dynamic interval decision index system is further amplified
So that system trajectory deviation greatly increases.And the part increased is more, it is seen that dynamic interval decision index system size is more
Ground affects this partial properties index of trajector deviation.
Dispatching method of the invention is compared with the dynamic interval model for not considering reference locus correction index.As a result
Referring to the signal of Fig. 8, do not consider that reference locus corrects index, system is then excellent according to performance indicator progress based on reference locus always
Change scheduling, it is difficult to guarantee global optimum of the system in prediction time domain.According to Fig. 8, the electric generating optimization of 18:00-19:00 and
On the natural gas generating optimization of 20:00-20:45, due to being corrected without reference to trajectory corrector index so that actual path with
The deviation of reference locus is larger, causes system cost larger, and final scheduling cost isCompared with S1, increase
Scheduling cost.
Influence of the different reference locus correction indexs to scheduling result is compared, considers different reference locus correction indexs,
5 groups of reference locus correction indexs are setThe value of triggering is respectively 0.04,0.08,0.12,0.16,0.20.If current state
Reference locus correction index value be less than the value, then do not trigger dynamic interval decision, if it is greater, then triggering.Performance refers to
Mark corrects signal of the result of variations referring to Fig. 9 of index with reference locus.The different numerical value of reference locus correction index will affect
In a few days to the correction situation planned a few days ago when Optimized Operation.When reference locus correction index is larger, to actual path and reference
The deviation nargin of track is larger, may lose certain necessary corrections, increases the trajector deviation in performance indicator.When reference rail
It is smaller to the deviation nargin of actual path and reference locus when mark correction index is smaller, some unnecessary schools may be executed
Just, increase the control energy and control cost in performance indicator.According to Fig. 9, with the increase of reference locus correction index, performance
Index first increases and then decreases.First reducing is the amplification because of deviation nargin, and control energy and control cost are minimized.After increase
It is because deviation nargin is further amplified so that the trajector deviation of system considerably increases.
The present embodiment is based on Model Predictive Control Theory, the integrated energy system dispatching method of research trends time interval.
Using trajector deviation-control energy-control cost as performance indicator, dynamic interval decision index system is established, it is ensured that based on being
The promptly and accurately scheduling of system time of day.Time of day and the accurate scheduling following embody:
1, the output prediction link based on Model Predictive Control feeds back system call, is based on system mode progressive
The optimization of energy index；
2, dynamic interval decision index system is capable of the size of reaction system state change, small and to system shadow when changing
When ringing little, it is possible to reduce scheduling times, when variation is big and big, the then publication dispatch command in time to systematic influence.
Establish reference locus correction index, it is ensured that system realizes the global optimization tune of incorporation time scale in prediction time domain
Degree.Elaborate the integrated energy system scheduling flow based on dynamic interval Model Predictive Control.
Example gives the scheduling scheme based on dispatching method of the present invention.
The technical solution of the present embodiment is used, it can be achieved that following the utility model has the advantages that considering with " trajector deviation-control energy-control
Be made this " be performance indicator dynamic interval integrated energy system dispatching method, be conducive to performance indicator have it is enough
Adjustable nargin copes with different scene demands；Dynamic interval decision index system is considered, being conducive to reflection system will occur
The time of day for even having occurred and that large error, provides decision index system for whether system needs to assign dispatch command；It considers
In a few days reference locus corrects index, is conducive to cope with actual scene and predicts that the deviation of scene, the real now forecast time domain overall situation are excellent
Change scheduling, facilitates the popularization and implementation of this programme.Adjustable nargin is that performance indicator includes trajector deviation, control energy and control
Cost three parts, requirement having differences property of the different scenes to three, some scenes are more demanding to trajector deviation, some fields
Scape is higher to control cost requirement, and the optimization range of tuning performance index can be carried out by adjusting the weight coefficient between three,
To adapt to more scene demands.
It should be noted that the above examples are only used to illustrate the technical scheme of the present invention and are not limiting, although referring to preferable
Embodiment describes the invention in detail, those skilled in the art should understand that, it can be to technology of the invention
Scheme is modified or replaced equivalently, and without departing from the spirit and scope of the technical solution of the present invention, should all be covered in this hair
In bright scope of the claims.
Claims (7)
1. a kind of integrated energy system dispatching method of dynamic interval, it is characterised in that: include the following steps,
The integrated energy system Scheduling Framework based on dynamic interval Model Predictive Control is established, the comprehensive energy system is obtained
The scheduling principle of system；
Establish the scheduling model a few days ago of the integrated energy system, obtain a few days ago operation plan as the reference rail in a few days dispatched
Mark；
In a few days scheduling " trajector deviation-control energy-control cost " performance indicator based on Model Predictive Control is established, as day
The objective function of interior optimizing scheduling；
Dynamic interval decision index system is established, judges whether scheduling interval needs to change；
Reference locus correction index is established, judges whether reference locus needs to correct；
Energy conversion unit, heat supply network model, gas pessimistic concurrency control and electric network model are established, integrated energy system schedule constraints item is obtained
Part obtains the in a few days scheduling scheme of integrated energy system in the case where meeting constraint condition.
2. the integrated energy system dispatching method of dynamic interval as described in claim 1, it is characterised in that: the foundation
Integrated energy system Scheduling Framework based on dynamic interval Model Predictive Control obtains pre- based on dynamic interval model
The scheduling principle of the integrated energy system of observing and controlling is further comprising the steps of,
It is dispatched a few days ago using the cost that runs minimized as target making hour service capacity plan in lower day and is issued.And it in a few days rolls
Optimized Operation is then roll using running power output as initial value excellent using operation plan a few days ago as reference locus w (k+i)
Change.According to the system running state sequence u (k) at k moment, the system mode sequence based on the prediction model prediction k+1 momentSystem mode sequenceOutput prediction link is collectively formed with the system mode error sequence e (k) at k momentOn the one hand output prediction link passes through dynamic interval decision index system ξt, judge whether that needing to assign amendment adjusts
The instruction of plan is spent, on the other hand corrects index by a few days reference locusJudge whether to need to reference locus w (k+i)
It is modified.
3. the integrated energy system dispatching method of dynamic interval as described in claim 1, it is characterised in that: the foundation
Integrated energy system scheduling model a few days ago, obtaining operation plan a few days ago as the reference locus in a few days dispatched further includes following step
Suddenly,
It dispatches a few days ago and optimizes scheduling by target of economic optimum, objective function is as follows:
In formula: T is the Optimized Operation period a few days ago, is taken for 24 hours, time interval 1h；U is generating set set；V is gas source collection
It closes；Eu,tFor the active power output of t moment generating set u, Gv,tFor the gas production of t moment gas source v；For the power generation of generating set u
Cost,For the Gas Prices of gas source v.
4. the integrated energy system dispatching method of dynamic interval as described in claim 1, it is characterised in that: the foundation
In a few days scheduling " trajector deviation-control energy-control cost " performance indicator based on Model Predictive Control, it is excellent as in a few days dispatching
The objective function of change is further comprising the steps of,
The reference locus w (k+i) dispatched a few days ago and output predictionAfter scene Recognition, according to performance indicator
It optimizes, performance indicator is the objective function of rolling optimization.Performance index function considers trajector deviation, control energy, control
Cost, as shown in formula:
Being expressed as matrix form is formula:
Wherein, T is prediction time domain, and M is control time domain, Q=diag [q1, q2..., qT], R=diag [r1, r2..., rM], S=
diag[s1, s2..., sM] be different scenes under weight matrix,To predict that output matrix, W (k+1) are with reference to rail
Mark matrix, Δ p (k) are controlling increment matrix, and C (k) is control Cost matrix.Prediction output matrix and reference locus matrix include
The active power output and gas production of equipment, controlling increment matrix are regulated quantity of the equipment in plan a few days ago, and control Cost matrix is
The cost of equipment supply load.
5. the integrated energy system dispatching method of dynamic interval as claimed in claim 2, it is characterised in that: the foundation
Dynamic interval decision index system, judge scheduling interval whether need to change it is further comprising the steps of,
Based on to integrated energy system overall situation operating cost in decision period forecasting time domain and operation global under in a few days reference locus
The deviation ratio of cost establishes dynamic interval decision index system ξt.Global cost is shown below:
S in formular,tIndicate that decision period and the global cost for predicting time domain, A indicate the unit matrix of 1 row T column, P to t moment
Indicate that decision period and prediction time domain regulate and control moment matrix, p in t momentlT,tIndicate that t moment Demand-side and system side resource l exist
Predict the regulation amount of time domain T, C indicates Cost matrix, clIndicate the cost coefficient regulated and controled to l.
Dynamic interval decision index system ξtShown in being expressed as follows:
Wherein, Sw,tIndicate the global cost of reference locus.
When the deviation between the actual path and reference locus after the optimizing scheduling based on Model Predictive Control transports the overall situation
When row cost produces bigger effect, dynamic interval decision will be triggered, changes control room every excellent to control time domain progress in advance
Change and execute dispatch command, when influencing smaller, then do not trigger decision, control room reduces unnecessary execution and refers to every will be late by
It enables.
6. the integrated energy system dispatching method of dynamic interval as described in claim 1, it is characterised in that: the foundation
Reference locus correct index, judge reference locus whether need to correct it is further comprising the steps of,
When the deviation that actual motion scene and prediction scene are likely to occur, it is intended to be the reference significance of reference locus a few days ago just
It reduces.Therefore actual path and reference locus are compared, establishes reference locus correction index.Operation plan is main a few days ago
Determine the start and stop and combination plan of unit, if therefore current time occur relatively large deviation, Unit Commitment should be adjusted in time
Or combination plan, then in a few days reference locus corrects indexIt can indicate are as follows:
Wherein,For actual path,For in a few days reference locus.
7. the integrated energy system dispatching method of dynamic interval as described in claim 1, it is characterised in that: the foundation
Energy conversion unit, heat supply network model, gas pessimistic concurrency control, electric network model obtain integrated energy system scheduling constraint, are meeting about
In the case where beam condition, the in a few days scheduling scheme for obtaining integrated energy system is further comprising the steps of,
The energy conversion unit considers energy balance constraint；
The heat supply network model consider node flow balance, node power fusion, load take characteristic, supply and return water temperature constraint and
Pipeline section heat-transfer character；
The gas pessimistic concurrency control considers pipeline flow constraint, gas source point constraint, flow equilibrium constraint, compressor constraint and node pressure
Force constraint；
The electric network model considers node power balance, unit output constraint, Climing constant and Branch Power Flow constraint.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910776326.6A CN110516951B (en) | 2019-08-22 | 2019-08-22 | Dynamic time interval comprehensive energy system scheduling method |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910776326.6A CN110516951B (en) | 2019-08-22 | 2019-08-22 | Dynamic time interval comprehensive energy system scheduling method |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN110516951A true CN110516951A (en) | 2019-11-29 |
| CN110516951B CN110516951B (en) | 2021-01-05 |
Family
ID=68627055
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201910776326.6A Active CN110516951B (en) | 2019-08-22 | 2019-08-22 | Dynamic time interval comprehensive energy system scheduling method |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN110516951B (en) |
Cited By (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN111898876A (en) * | 2020-07-13 | 2020-11-06 | 江苏方天电力技术有限公司 | Comprehensive energy regulation and control method considering air pipe network storage |
| CN112186752A (en) * | 2020-09-24 | 2021-01-05 | 国网辽宁省电力有限公司葫芦岛供电公司 | Single-target multi-time-period accurate adjustment method |
| CN112434420A (en) * | 2020-11-20 | 2021-03-02 | 国网山东省电力公司电力科学研究院 | Time synchronization and data interaction method for hybrid simulation of comprehensive energy system |
| CN113486532A (en) * | 2021-07-22 | 2021-10-08 | 东南大学 | Dynamic safety control method for electric heating comprehensive energy system |
| CN113840384A (en) * | 2021-11-29 | 2021-12-24 | 成都成电光信科技股份有限公司 | Variable step length scheduling method for time trigger message in TT-FC network |
| CN115562029A (en) * | 2022-10-17 | 2023-01-03 | 杭州天然气有限公司 | Intelligent control method and system for natural gas turbine expansion generator set |
Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20170060109A1 (en) * | 2013-06-06 | 2017-03-02 | International Business Machines Corporation | Managing time-substitutable electricity usage using dynamic controls |
| CN107732982A (en) * | 2017-10-20 | 2018-02-23 | 河海大学 | Consider the integrated energy system Multiple Time Scales dispatching method of Model Predictive Control |
| CN108932560A (en) * | 2018-06-13 | 2018-12-04 | 天津大学 | Garden integrated energy system Optimization Scheduling based on Model Predictive Control |
| CN109787227A (en) * | 2019-02-02 | 2019-05-21 | 国网江苏省电力有限公司南京供电分公司 | One kind is provided multiple forms of energy to complement each other system Multiple Time Scales Optimization Scheduling |
| CN110137942A (en) * | 2019-04-23 | 2019-08-16 | 河海大学 | Multiple Time Scales flexible load rolling scheduling method and system based on Model Predictive Control |
-
2019
- 2019-08-22 CN CN201910776326.6A patent/CN110516951B/en active Active
Patent Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20170060109A1 (en) * | 2013-06-06 | 2017-03-02 | International Business Machines Corporation | Managing time-substitutable electricity usage using dynamic controls |
| CN107732982A (en) * | 2017-10-20 | 2018-02-23 | 河海大学 | Consider the integrated energy system Multiple Time Scales dispatching method of Model Predictive Control |
| CN108932560A (en) * | 2018-06-13 | 2018-12-04 | 天津大学 | Garden integrated energy system Optimization Scheduling based on Model Predictive Control |
| CN109787227A (en) * | 2019-02-02 | 2019-05-21 | 国网江苏省电力有限公司南京供电分公司 | One kind is provided multiple forms of energy to complement each other system Multiple Time Scales Optimization Scheduling |
| CN110137942A (en) * | 2019-04-23 | 2019-08-16 | 河海大学 | Multiple Time Scales flexible load rolling scheduling method and system based on Model Predictive Control |
Cited By (9)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN111898876A (en) * | 2020-07-13 | 2020-11-06 | 江苏方天电力技术有限公司 | Comprehensive energy regulation and control method considering air pipe network storage |
| CN112186752A (en) * | 2020-09-24 | 2021-01-05 | 国网辽宁省电力有限公司葫芦岛供电公司 | Single-target multi-time-period accurate adjustment method |
| CN112434420A (en) * | 2020-11-20 | 2021-03-02 | 国网山东省电力公司电力科学研究院 | Time synchronization and data interaction method for hybrid simulation of comprehensive energy system |
| CN112434420B (en) * | 2020-11-20 | 2023-09-05 | 国网山东省电力公司电力科学研究院 | Time synchronization and data interaction method for hybrid simulation of comprehensive energy system |
| CN113486532A (en) * | 2021-07-22 | 2021-10-08 | 东南大学 | Dynamic safety control method for electric heating comprehensive energy system |
| CN113486532B (en) * | 2021-07-22 | 2023-12-19 | 东南大学 | Dynamic safety control method for electric heating comprehensive energy system |
| CN113840384A (en) * | 2021-11-29 | 2021-12-24 | 成都成电光信科技股份有限公司 | Variable step length scheduling method for time trigger message in TT-FC network |
| CN113840384B (en) * | 2021-11-29 | 2022-03-08 | 成都成电光信科技股份有限公司 | Variable step length scheduling method for time trigger message in TT-FC network |
| CN115562029A (en) * | 2022-10-17 | 2023-01-03 | 杭州天然气有限公司 | Intelligent control method and system for natural gas turbine expansion generator set |
Also Published As
| Publication number | Publication date |
|---|---|
| CN110516951B (en) | 2021-01-05 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN110516951A (en) | A kind of integrated energy system dispatching method of dynamic interval | |
| Zhang et al. | Optimal operation of integrated electricity and heat system: A review of modeling and solution methods | |
| Yang et al. | Economic-emission dispatch problem in integrated electricity and heat system considering multi-energy demand response and carbon capture Technologies | |
| CN108009693A (en) | Grid-connected micro-capacitance sensor dual blank-holder based on two-stage demand response | |
| CN106384168A (en) | Multi-objective coordinated optimization scheduling model for different power sources | |
| CN106505635A (en) | Abandon the minimum active power dispatch model of wind and scheduling system | |
| CN106682808A (en) | Online rolling optimization scheduling model | |
| CN109617142A (en) | A kind of CCHP type micro-capacitance sensor Multiple Time Scales Optimization Scheduling and system | |
| CN108876066A (en) | Heat supply network off-the-line operating scheme real-time optimization method and its system | |
| Ma et al. | A compact model to coordinate flexibility and efficiency for decomposed scheduling of integrated energy system | |
| CN109447323A (en) | It is a kind of meter and node caloric value integrated energy system two stages capacity collocation method | |
| CN104573875A (en) | Low-carbon power source and power grid optimization planning method | |
| Wu et al. | A two-stage rolling optimization strategy for park-level integrated energy system considering multi-energy flexibility | |
| Zhao et al. | Model predictive based frequency control of power system incorporating air-conditioning loads with communication delay | |
| CN109286208A (en) | A kind of integrated energy system dispatching method and system | |
| CN110400056A (en) | Cascade hydropower based on honourable space-time characterisation Optimization Scheduling and device a few days ago | |
| CN115498668A (en) | Optimization method of comprehensive energy system | |
| CN116308881A (en) | Multi-time scale scheduling method for comprehensive energy system utilizing heat supply pipe network for heat storage | |
| Liu et al. | Multi-objective optimization strategy of integrated electric-heat system based on energy storage situation division | |
| Wang et al. | Event-triggered online energy flow control strategy for regional integrated energy system using Lyapunov optimization | |
| Wang et al. | A time-scale adaptive dispatching strategy considering the matching of time characteristics and dispatching periods of the integrated energy system | |
| Hu et al. | A hierarchical transactive energy management framework for optimizing the reserve profile in district energy systems | |
| Li et al. | Stochastic-robust planning optimization method based on tracking-economy extreme scenario tradeoff for CCHP multi-energy system | |
| CN116502921A (en) | Park comprehensive energy system optimization management system and coordination scheduling method thereof | |
| CN115713151A (en) | Data center source-load collaborative optimization scheduling method, system, equipment and medium |
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 |