CN109472413A - Optimal scheduling method of park comprehensive energy system considering heat pipe network transmission characteristics - Google Patents
Optimal scheduling method of park comprehensive energy system considering heat pipe network transmission characteristics Download PDFInfo
- Publication number
- CN109472413A CN109472413A CN201811353627.XA CN201811353627A CN109472413A CN 109472413 A CN109472413 A CN 109472413A CN 201811353627 A CN201811353627 A CN 201811353627A CN 109472413 A CN109472413 A CN 109472413A
- Authority
- CN
- China
- Prior art keywords
- pipe network
- pipeline
- energy
- source station
- moment
- 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.)
- Pending
Links
- 230000005540 biological transmission Effects 0.000 title claims abstract description 49
- 238000000034 method Methods 0.000 title abstract description 14
- 238000005457 optimization Methods 0.000 claims abstract description 18
- 230000008878 coupling Effects 0.000 claims abstract description 12
- 238000010168 coupling process Methods 0.000 claims abstract description 12
- 238000005859 coupling reaction Methods 0.000 claims abstract description 12
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims description 63
- 239000012141 concentrate Substances 0.000 claims description 28
- 238000004378 air conditioning Methods 0.000 claims description 9
- 230000005611 electricity Effects 0.000 claims description 9
- 238000003860 storage Methods 0.000 claims description 9
- 238000010438 heat treatment Methods 0.000 claims description 8
- 239000008400 supply water Substances 0.000 claims description 6
- 230000008859 change Effects 0.000 claims description 4
- 238000006243 chemical reaction Methods 0.000 claims description 4
- 238000004146 energy storage Methods 0.000 claims description 4
- 238000005286 illumination Methods 0.000 claims description 4
- 239000000203 mixture Substances 0.000 claims description 4
- 239000003795 chemical substances by application Substances 0.000 claims description 3
- 238000012546 transfer Methods 0.000 claims description 3
- 241000208340 Araliaceae Species 0.000 claims 1
- 235000005035 Panax pseudoginseng ssp. pseudoginseng Nutrition 0.000 claims 1
- 235000003140 Panax quinquefolius Nutrition 0.000 claims 1
- 235000008434 ginseng Nutrition 0.000 claims 1
- 238000001816 cooling Methods 0.000 abstract description 18
- 230000009466 transformation Effects 0.000 abstract 1
- 238000005057 refrigeration Methods 0.000 description 10
- 238000010586 diagram Methods 0.000 description 8
- 239000000498 cooling water Substances 0.000 description 3
- 239000013256 coordination polymer Substances 0.000 description 3
- LYCAIKOWRPUZTN-UHFFFAOYSA-N Ethylene glycol Chemical compound OCCO LYCAIKOWRPUZTN-UHFFFAOYSA-N 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 238000003912 environmental pollution Methods 0.000 description 2
- 230000006872 improvement Effects 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 238000011160 research Methods 0.000 description 2
- 238000009825 accumulation Methods 0.000 description 1
- 230000003111 delayed effect Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 239000006185 dispersion Substances 0.000 description 1
- 235000013399 edible fruits Nutrition 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 239000002803 fossil fuel Substances 0.000 description 1
- WGCNASOHLSPBMP-UHFFFAOYSA-N hydroxyacetaldehyde Natural products OCC=O WGCNASOHLSPBMP-UHFFFAOYSA-N 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 238000002360 preparation method Methods 0.000 description 1
- 230000008569 process Effects 0.000 description 1
- 230000000750 progressive effect Effects 0.000 description 1
- 230000009467 reduction Effects 0.000 description 1
- 230000009897 systematic effect Effects 0.000 description 1
- 238000010977 unit operation Methods 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/04—Forecasting or optimisation specially adapted for administrative or management purposes, e.g. linear programming or "cutting stock problem"
-
- 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
- G06Q10/06312—Adjustment or analysis of established resource schedule, e.g. resource or task levelling, or dynamic rescheduling
-
- 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
- G06Q10/06315—Needs-based resource requirements planning or 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
Landscapes
- Business, Economics & Management (AREA)
- Human Resources & Organizations (AREA)
- Engineering & Computer Science (AREA)
- Economics (AREA)
- Strategic Management (AREA)
- Entrepreneurship & Innovation (AREA)
- Tourism & Hospitality (AREA)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- General Business, Economics & Management (AREA)
- Marketing (AREA)
- Physics & Mathematics (AREA)
- Operations Research (AREA)
- Quality & Reliability (AREA)
- Game Theory and Decision Science (AREA)
- Development Economics (AREA)
- Educational Administration (AREA)
- Health & Medical Sciences (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)
Abstract
The invention discloses a campus comprehensive energy system optimization scheduling method considering heat pipe network transmission characteristics, which aims to solve the problem of energy supply optimization scheduling of a campus comprehensive energy system, fully considers the transmission loss and delay influence of a pipe network, the coordination of multiple devices of a centralized energy station and the coupling relation between the centralized energy station and an energy supply pipe network, establishes a campus comprehensive energy system optimization scheduling model comprising the centralized energy station and the energy supply pipe network, and calls a relevant mathematical solver to solve through nonlinear constraint linear transformation to obtain a day-ahead cooling supply scheduling plan.
Description
Technical field
The present invention relates to a kind of integrated energy system Optimization Schedulings, transmit more particularly to a kind of consideration hot pipe network special
The garden integrated energy system Optimization Scheduling of property.
Background technique
The energy is the basis that the mankind depend on for existence and development, with fossil fuel shortage and problem of environmental pollution increasingly
Seriously, how to improve the environmental pollution during social efficiency of energy utilization, reduction energy becomes if the whole world pays close attention to jointly
Topic.Garden integrated energy system (community integrated energy system, CIES) is met the tendency of in this background
And it gives birth to.Garden integrated energy system is integrated with the energy storage system and coupling dress of the various energy resources such as electric, hot and cold and multiple-energy-source form
It sets, by carrying out organic coordination and optimal control, it can be achieved that a variety of energy to multiple links such as production of energy, conversion, transmission, consumption
The high efficiency of source system entirety utilizes, economy with higher and energy supply reliability.
The importance that summer cooling, Winter heat supply are energized as regional complex energy resource system needs to formulate rationally effectively
Scheduling strategy.In the energy supply phase, is mostly generated by concentration energy source station and be used for air conditioner refrigerating/heat air conditioner cold water/hot water, pass through water supply
Pipeline to each building, then load side fan coil system work reaches cooling/heating effect；Air-conditioning after temperature change
Hot and cold water returns by water return pipeline and concentrates energy source station.Garden integrated energy system has the characteristics that source-lotus dispersion, optimization are adjusted
Degree should consider to energize the coordination optimization of side concentration energy source station, the influence that consideration energy supply network transmits again.At this time, it may be necessary to sufficiently
Study pipe network for thermal characteristics and transmission characteristic, being capable of accurate description pipe network shape based on thermodynamics and the foundation of pipe network basic principle
State, the mode for meeting actual motion feature.
Currently, the research of garden integrated energy system is mostly individually to concentrate energy source station as research object, less consideration is supplied
The influence of energy pipeline transmission, and have no the dynamic transmission model for considering pipe network energy loss and transmission delay.Therefore, it is badly in need of one kind
The Optimization Scheduling for concentrating energy source station and the multiple energy supply links of energy supply network is taken into account, considers pipe network energy loss and transmission delay
Dynamic characteristic, optimization concentrate energy source station and energy supply network coordinated operation, meet to economic and reliable garden integrated energy system
With energy demand.
Summary of the invention
The technical problem to be solved by the present invention is to the Optimized Operation for garden integrated energy system for the cold heat phase is asked
Topic can establish the garden integrated energy system Optimal Operation Model for considering hot pipe network transmission characteristic, comprehensively consider the concentration energy
It stands a variety of constraints such as operation, pipe network dynamic transmission, carries out including concentrating energy source station, the coordination for the energy resource system for energizing pipeline excellent
Change, it is final to formulate Optimized Operation scheme.
To achieve the above object, the embodiment of the invention provides following technical solutions:
The embodiment of the invention provides a kind of garden integrated energy system Optimized Operation sides for considering hot pipe network transmission characteristic
Method, comprising:
The selected garden integrated energy system of S1, basis, inputs electricity price information, reads electric load, cold heat load, illumination
The predicted value of intensity, energy source station equipment composition, equipment operating parameter, the currently stored heat of thermal storage equipment, unit open in input set
Shutdown expenses, pipe network structure and parameter, system call interval, heating medium parameter and concentrate the variables such as energy source station exit water temperature bound or
Parameter；
S2, according to the S1 provide the garden integrated energy system structure and parameter, establish consider hot pipe network pass
The garden integrated energy system Optimal Operation Model of defeated characteristic, comprising: set one scheduling week of the garden integrated energy system
The minimum objective function of the sum of operating cost and Unit Commitment expense in phase considers to concentrate energy source station operation constraint, heat pipe respectively
The constraint of net dynamic transmission concentrates energy source station-pipe network coupling constraint, pipe network-load coupling constraint, the cold heat load equilibrium of supply and demand about
Beam and the constraint of the electric load equilibrium of supply and demand；
Relevant nonlinear constraint carries out in S3, the garden integrated energy system Optimal Operation Model for obtaining the S2
Linearisation conversion, establishes mixed integer linear programming model and related solver is called to be solved；
The solving result of mixed integer linear programming model described in S4, the output S3, including operating cost, energy supply master
Machine start and stop instruction, operating condition, energy supply power, energy storage equipment energy supply instruction and energy supply power concentrate energy source station air-conditioning cold/hot water
Supply water temperature etc..
Preferably, the pipeline transmission delay model in the hot pipe network dynamic transmission constraint in the S2 specifically uses down
Formula indicates:
In above formula, ni,tFlow out the minimum scheduling interval of pipeline i passed through for the first time by the medium of t moment flow ipe i；
Si,t、Ri,tRespectively t moment is to t+ni,tMoment, t moment to t+ni,tThe cold water quality of -1 moment flow ipe i；ρ,Ai、LiPoint
It Wei not air conditioner cold water density, the length of the cross-sectional area of i-th pipeline and i-th pipeline；qi,tFor i-th pipeline delay of t moment
Weighting coefficient；mi,tFor the cold water flow of i-th pipeline of t moment；Δ t is system call interval；Respectively t
The head end temperature and terminal temperature of i-th pipeline of moment；X ∈ { s, r }, s, r respectively represent water supply line and water return pipeline；λiFor
The heat transfer coefficient of pipeline unit length；CwFor the specific heat capacity of heating agent；TaFor pipeline ambient temperature, mi,t+kThat represent is t+k
The cold water flow of i-th pipeline of moment.
Preferably, the concentration energy source station-pipe network coupling constraint in the S2 specifically uses following formula to indicate:
In above formula,Energy source station is concentrated to energize power for t moment；Respectively t moment concentrates energy source station
Air conditioner cold water return water temperature and concentration energy source station air-conditioning supply water temperature；Respectively it is connected with energy source station is concentrated
The temperature of water supply line headend node and the temperature of water return pipeline endpoint node；Tces,min、Tces,maxRespectively energy source station is concentrated to go out
Mouth air conditioner cold water lowest temperature and the upper limit；ΩcesSet for the pipeline being connected with concentration energy source station.
A kind of garden integrated energy system Optimized Operation side considering hot pipe network transmission characteristic provided in an embodiment of the present invention
Method fully considers pipe network transmission loss and delay shadow based on the energy supply Optimal Scheduling for solving garden integrated energy system
It rings, concentrate the cooperation of the more equipment of energy source station and concentrate energy source station-energy supply pipe network coupled relation, establish comprising concentrating the energy
It stands, energize the garden integrated energy system Optimal Operation Model of pipe network, pass through the linearisation conversion of nonlinear restriction, call related
Mathematical device is solved, and cooling supply operation plan a few days ago is obtained.
Detailed description of the invention
In order to more clearly explain the embodiment of the invention or the technical proposal in the existing technology, to embodiment or will show below
There is attached drawing needed in technical description to be briefly described, it should be apparent that, the accompanying drawings in the following description is only this
Some embodiments of invention for those of ordinary skill in the art without creative efforts, can be with
It obtains other drawings based on these drawings.
Fig. 1 is that a kind of garden integrated energy system for considering hot pipe network transmission characteristic disclosed by the embodiments of the present invention optimizes tune
Spend method flow schematic diagram；
Fig. 2 is pipeline transmission delay schematic diagram；
Fig. 3 is to concentrate energy source station energy supply structure chart；
Fig. 4 is system energy supply network structure chart；
Fig. 5 is building room temperature variation diagram when not considering pipeline transmission characteristic；
Fig. 6 is building room temperature variation diagram when only considering conduits of energy loss.
Specific embodiment
Following will be combined with the drawings in the embodiments of the present invention, and technical solution in the embodiment of the present invention carries out clear, complete
Site preparation description, it is clear that described embodiments are only a part of the embodiments of the present invention, instead of all the embodiments.It is based on
Embodiment in the present invention, it is obtained by those of ordinary skill in the art without making creative efforts every other
Embodiment shall fall within the protection scope of the present invention.
The embodiment of the invention discloses a kind of garden integrated energy system Optimized Operation sides for considering hot pipe network transmission characteristic
Method can establish the garden for considering hot pipe network transmission characteristic for garden integrated energy system for the Optimal Scheduling of cold heat phase
Area's integrated energy system Optimal Operation Model comprehensively considers and concentrates a variety of constraints such as energy source station operation, pipe network dynamic transmission, carries out
Coordination optimization including concentrating energy source station, the energy resource system for energizing pipeline, it is final to formulate Optimized Operation scheme.
Referring to Figure 1, Fig. 1 is that the embodiment of the invention discloses a kind of garden comprehensive energies for considering hot pipe network transmission characteristic
System optimization dispatching method flow diagram, this method comprises:
The selected garden integrated energy system of S1, basis, inputs electricity price information, reads electric load, cold heat load, illumination
The predicted value of intensity, energy source station equipment composition, equipment operating parameter, the currently stored heat of thermal storage equipment, unit open in input set
Shutdown expenses, pipe network structure and parameter, system call interval, heating medium parameter and concentrate the variables such as energy source station exit water temperature bound or
Parameter；
The structure and parameter for the garden integrated energy system that S2, foundation S1 are provided, establishes and considers hot pipe network transmission characteristic
Garden integrated energy system Optimal Operation Model, comprising: set operating cost in one dispatching cycle of garden integrated energy system
With the minimum objective function of the sum of Unit Commitment expense, consider to concentrate energy source station operation constraint, hot pipe network dynamic transmission about respectively
Beam concentrates energy source station-pipe network coupling constraint, pipe network-load coupling constraint, the constraint of the cold heat load equilibrium of supply and demand and electric load to supply
Need Constraints of Equilibrium.
Specifically, in the present embodiment, operating cost and Unit Commitment take in one dispatching cycle of garden integrated energy system
With the sum of minimum objective function specifically indicated using following formula:
In above formula, NTFor the total linear spacing number of a dispatching cycle,Indicate t moment purchase electricity price,For system and outside
The power of interconnecting ties, EiFor the switching cost of different hosts, S is cluster tool, comprising: earth source heat pump, conventional cold water master
Machine, Double-working-condition host,Respectively t moment jth platform earth source heat pump, conventional cold water main unit, Double-working-condition master
The start and stop state of machine, ΩHP、ΩWC、ΩDCThe respectively set of earth source heat pump, conventional cold water main unit, Double-working-condition host, ΩiInclude
ΩHP、ΩWC、ΩDC.The relevant binary variable of model represents powered mode/equipment for 1 and is in execution/starting state, is 0 generation
Table does not execute/shutdown status, similarly hereinafter.
Energy source station operation constraint is concentrated to specifically include following type:
First seed type: the operation constraint of earth source heat pump unit is specifically indicated using following formula:
In above formula,For i-th earth source heat pump cooling supply power of t moment；N HPFor ground-source heat pump host number；The respectively minimum of heat pump main frame, maximum refrigeration work consumption；For t moment heat pump unit power consumption；
For i-th performance coefficient of heat pump (COP)；PHP,CWPAnd PHP,CPThe respectively volume of heat pump main frame interlocking chilled water pump and cooling water pump
Determine electric power.
Second of type: chiller unit operation constraint is specifically indicated using following formula:
In above formula,For t moment i-th conventional cold water main unit cooling supply power；NWCFor conventional cold water main unit number；The respectively upper and lower limit of its refrigeration work consumption；For t moment chiller unit power consumption；It is i-th
The platform routine cold water main unit coefficient of performance；PWC,CWP、PWC,CPAnd PWC,CTRespectively conventional cold water main unit interlocks chilled water pump, cooling water
The specified electric power of pump and open cooling tower.
Third seed type: ice-storage system operation constraint is specifically indicated using following formula:
In above formula,For t moment ice-storage system, Ice Storage Tank refrigeration work consumption；When respectively t
Carve i-th Double-working-condition host refrigeration, ice making power；For the upper and lower limit of Double-working-condition host refrigeration work consumption；The upper and lower limit of power is made ice for it；NDCFor Double-working-condition host number； For i-th duplex of t moment
Condition host refrigeration, ice making operation mode；For t moment Double-working-condition unit refrigeration, ice making operation mode；For t
Moment Ice Storage Tank stores cooling capacity；W IT、The upper and lower limit of cooling capacity is stored for Ice Storage Tank；εITRate is let cool certainly for Ice Storage Tank；Δ t is
Scheduling interval；The upper limit of the power is let cool for Ice Storage Tank；For the separate unit chilled water pump refrigeration work consumption upper limit；When for t
Carve i-th ice-storage system chilled water pump start and stop state；For t moment ice-storage system power consumption；Respectively Double-working-condition host refrigeration, the ice making coefficient of performance, PEP、PDC,CP、PDC,CT、PIS,CWPRespectively second
Glycol solution pump, cooling water pump, open cooling tower, chilled water pump specified electric power.
The constraint of pipe network dynamic transmission specifically indicates (wherein to pass comprising the pipeline in the constraint of hot pipe network dynamic transmission using following formula
Defeated delay model):
In above formula, mi,tFor the cold water flow of i-th pipeline of t moment；Respectively i-th pipeline of t moment
Head end temperature and terminal temperature；X ∈ { s, r }, s, r respectively represent water supply line and water return pipeline；For the intersection of t moment i-th
The mixing temperature of node；Endpoint node, headend node are node i respectively in confession/return water network, that is, flow into,
Flow out the pipeline set of node, ni,tBetween the minimum scheduling passed through for flowing out pipeline i for the first time by the medium of t moment flow ipe i
Every；Si,t、Ri,tRespectively t moment is to t+ni,tMoment, t moment to t+ni,tThe cold water quality of -1 moment flow ipe i；ρ,Ai、Li
Respectively air conditioner cold water density, the length of the cross-sectional area of i-th pipeline and i-th pipeline；qi,tProlong for i-th pipeline of t moment
When weighting coefficient；Δ t is system call interval；λiFor the heat transfer coefficient of pipeline unit length；CwFor the specific heat capacity of heating agent；TaFor
Pipeline ambient temperature, mi,t+kWhat is represented is the cold water flow of i-th pipeline of t+k moment.In the constraint of hot pipe network dynamic transmission
Pipeline transmission delay model refer to Fig. 2, Fig. 2 is pipeline transmission delay schematic diagram.
Energy source station-pipe network coupling constraint is concentrated specifically following formula to be used to indicate:
In above formula,Energy source station is concentrated to energize power for t moment；Respectively t moment concentrates energy source station
Air conditioner cold water return water temperature and concentration energy source station air-conditioning supply water temperature；Respectively it is connected with energy source station is concentrated
The temperature of water supply line headend node and the temperature of water return pipeline endpoint node；Tces,min、Tces,maxRespectively energy source station is concentrated to go out
Mouth air conditioner cold water lowest temperature and the upper limit；ΩcesSet for the pipeline being connected with concentration energy source station.
Pipe network-load coupling constraint specifically uses following formula to indicate:
In above formula,It is t moment pipeline to the injecting power of building i,Respectively building i entrance, outlet
Locate heat medium temperature；Respectively the connected water supply line endpoint node of building, water return pipeline headend node temperature
Degree；Set for the pipeline being connected with building i.
The constraint of the cold heat load equilibrium of supply and demand is specifically indicated using following formula:
The constraint of the electric load equilibrium of supply and demand is specifically indicated using following formula:
In above formula,Respectively t moment photovoltaic system output power, dominant eigenvalues,For interconnection
Maximum allowable power value,For t moment system electric load.
S3, relevant nonlinear constraint in garden integrated energy system Optimal Operation Model that S2 is obtained is subjected to linearisation turns
It changes, establish mixed integer linear programming model and related solver is called to be solved.
For seizing the opportunity item, target side there are binary variable and continuous variable after formula (2), (5), (9), (10) expansion
Item is seized the opportunity there are binary variable and binary variable in journey (1), by introducing auxiliary variable and constraining these are non-linear
Item linearisation.After linearisation, this optimization problem is converted to a Mixed integer linear programming.
It is specifically divided into following steps:
(1) to nonlinear terms U, r, wherein U is binary variable, and r is continuous positive variable, and the upper limit of r isIt can introduce auxiliary
Variable R is helped to replace nonlinear terms, and R meets following constraint:
(2) to nonlinear terms U1、U2, wherein U1、U2For binary variable, it is non-linear to introduce assistant binary variable U replacement
, and U meets following constraint:
The solving result of mixed integer linear programming model in S4, output S3, including operating cost, energy supply host start and stop refer to
It enables, operating condition, energy supply power, energy storage equipment energy supply instructs and energy supply power, concentration energy source station air-conditioning cold/hot water supply water temperature
Deng.
It is adjusted as it can be seen that a kind of garden integrated energy system for considering hot pipe network transmission characteristic provided in an embodiment of the present invention optimizes
Degree method fully considers pipe network transmission loss and prolongs based on the energy supply Optimal Scheduling for solving garden integrated energy system
The cooperation of the more equipment of energy source station is influenced, concentrated late and concentrates energy source station-energy supply pipe network coupled relation, is established comprising concentrating
Energy source station, the garden integrated energy system Optimal Operation Model for energizing pipe network are converted by the linearisation of nonlinear restriction, are called
Related mathematics solver is solved, and cooling supply operation plan a few days ago is obtained.
Technical solution provided in an embodiment of the present invention is described in detail below with reference to practical application scene, is inputted first
Electricity price information, electric load predicted value dispatching cycle of system one, cooling load prediction value, intensity of illumination predicted value；Then input set
Middle energy source station equipment composition, equipment operating parameter, the currently stored cooling capacity of cold-storage apparatus, Unit Commitment expense, pipeline parameter, system
Scheduling interval, heating medium parameter and the initial value for concentrating the variables such as energy source station exit water temperature bound or parameter.In this system, by outside
Power grid and photovoltaic system meet electricity needs；Energy source station is concentrated to generate air conditioner cold water by energy supply pipeline to each building,
Meet cooling needs by fan coil.Concentrating energy source station includes: that 3 earth source heat pumps, 2 conventional cold water main units and one group of ice store
Refrigeration subsystem (two Double-working-condition hosts, an Ice Storage Tank).Concentrate energy source station for air-cooled structure as shown in figure 3, Fig. 3 is to concentrate the energy
It stands and energizes structure chart；Detail parameters are shown in Table 1, and energy supply network structure and parameter is shown in that Fig. 4 and table 2, Fig. 4 are system energy supply network structures
Figure, 2 energy supply network parameter information of table；Ice Storage Tank storage cooling capacity initial value is taken as 0；Ground source pump, conventional cold water main unit, Double-working-condition host
Starting expense is respectively 40.0 yuan/time, 120.0 yuan/time, 120.0 yuan/time；30min is divided between system call；Air-conditioning cold watertight
Degree and specific heat capacity are respectively 1000kg/m3,4.2kJ/ (kg DEG C)；Concentrate energy source station outlet air conditioner cold water at a temperature of, the upper limit difference
It is 5 DEG C, 8 DEG C；1.35 yuan/kWh of electricity price when peak (8:00-11:00,
18:00-23:00), 0.47 yuan/kWh of electricity price (00:00-7:00,23:00-00:00) when paddy, usually electricity price 0.89
Member/kWh (7:00-8:00,11:00-18:00).
The characteristics of for system by energy supply pipeline cooling supply, operating cost of the comparison using different scheduling strategy systems, knot
Fruit is shown in Table 3, wherein strategy 1 uses traditional Economic Scheduling Policy (dynamic transmission characteristics for not considering pipeline), strategy 2 considers
The energy loss of pipeline, but do not consider the influence of its delay, strategy 3 considers the energy loss and delayed impact of pipeline simultaneously.
Influence of the different scheduling strategies to running is mainly reflected in the variation of load side building room temperature, using strategy
The variation of load side (building) room temperature is shown in that Fig. 5, Fig. 5 are building room temperature variation diagrams when not considering pipeline transmission characteristic when 1；Using strategy
The variation of building room temperature is shown in that Fig. 6, Fig. 6 are building room temperature variation diagrams when only considering conduits of energy loss when 2.
Executing the computer hardware environment that optimization calculates is Intel (R) Xeon (R) CPU E5-2603, and dominant frequency is
1.60GHz inside saves as 8GB；Software environment is 10 operating system of Windows.
Compare systematic running cost under different scheduling strategies can be seen that consider pipeline after, due to the shadow of its transmission loss
It rings, operating cost can suitably increase, but due to considering pipe network transmission process in actual motion, Optimized Operation result will more
Meet practical application；After considering delay pipeline, then there is injustice between pipe network injecting power, building consumption power and loss power
It measures, there is pipe network certain accumulation of energy characteristic can reduce operating cost to a certain extent in conjunction with tou power price information at this time.
Comparing load side (building) indoor temperature change generated in case under different scheduling strategies can be seen that not while considering conduits of energy
Under the strategy of loss and transmission delay, building cooling supply power can deviate requirements (22 DEG C) when actual motion, these can be to using
Family is affected greatly with energy satisfaction；After the dynamic transmission characteristics for fully considering pipe network, practical building cooling supply power will be equal to
Theory demands value, room temperature are maintained at 22 DEG C of ideal temperature.As it can be seen that considering the garden comprehensive energy system of hot pipe network transmission characteristic
Optimization Scheduling unite close to running, can preferably meet user with energy demand, energy supply reliability with higher.
Table 1 concentrates energy station structure and parameter
2 energy supply network parameter information of table
The different scheduling strategy operating cost comparisons of table 3
Above to a kind of garden integrated energy system Optimized Operation of consideration hot pipe network transmission characteristic provided herein
Method is described in detail.Specific examples are used herein to illustrate the principle and implementation manner of the present application, with
The explanation of upper embodiment is merely used to help understand the present processes and its core concept.It should be pointed out that being led for this technology
For the those of ordinary skill in domain, under the premise of not departing from the application principle, can also to the application carry out it is several improvement and
Modification, these improvement and modification are also fallen into the protection scope of the claim of this application.
Each embodiment is described in a progressive manner in specification, the highlights of each of the examples are with other realities
The difference of example is applied, the same or similar parts in each embodiment may refer to each other.For device disclosed in embodiment
Speech, since it is corresponded to the methods disclosed in the examples, so being described relatively simple, related place is referring to method part illustration
?.
Claims (3)
1. a kind of garden integrated energy system Optimization Scheduling for considering hot pipe network transmission characteristic characterized by comprising
The selected garden integrated energy system of S1, basis, inputs electricity price information, reads electric load, cold heat load, intensity of illumination
Predicted value, energy source station equipment composition, equipment operating parameter, the currently stored heat of thermal storage equipment, Unit Commitment take in input set
With, pipe network structure and parameter, system call interval, heating medium parameter and concentrate variables or the ginsengs such as energy source station exit water temperature bound
Number；
S2, according to the structure and parameter of the S1 garden integrated energy system provided, establish and consider that hot pipe network transmission is special
The garden integrated energy system Optimal Operation Model of property, comprising: set in one dispatching cycle of the garden integrated energy system
The minimum objective function of the sum of operating cost and Unit Commitment expense considers to concentrate energy source station operation constraint, hot pipe network dynamic respectively
State transmission constraint, concentrate energy source station-pipe network coupling constraint, pipe network-load coupling constraint, the cold heat load equilibrium of supply and demand constraint and
The constraint of the electric load equilibrium of supply and demand；
Relevant nonlinear constraint carries out linear in S3, the garden integrated energy system Optimal Operation Model for obtaining the S2
Change conversion, establishes mixed integer linear programming model and related solver is called to be solved；
The solving result of mixed integer linear programming model described in S4, the output S3, including operating cost, energy supply host open
Stop instruction, operating condition, energy supply power, energy storage equipment energy supply instruction and energy supply power concentrate energy source station air-conditioning cold/hot water to supply water
Temperature etc..
2. the garden integrated energy system Optimization Scheduling according to claim 1 for considering hot pipe network transmission characteristic,
It is characterized in that, the pipeline transmission delay model in hot pipe network dynamic transmission constraint in the S2 specifically uses following formula table
Show:
In above formula, ni,tFlow out the minimum scheduling interval of pipeline i passed through for the first time by the medium of t moment flow ipe i；Si,t、
Ri,tRespectively t moment is to t+ni,tMoment, t moment to t+ni,tThe cold water quality of -1 moment flow ipe i；ρ,Ai、LiRespectively
The length of air conditioner cold water density, the cross-sectional area of i-th pipeline and i-th pipeline；qi,tFor i-th pipeline Delay weighting of t moment
Coefficient；mi,tFor the cold water flow of i-th pipeline of t moment；Δ t is system call interval；Respectively t moment
The head end temperature and terminal temperature of i-th pipeline；X ∈ { s, r }, respectively represents water supply line and water return pipeline；λiFor pipeline list
The heat transfer coefficient of bit length；CwFor the specific heat capacity of heating agent；TaFor pipeline ambient temperature, mi,t+kWhat is represented is the t+k moment i-th
The cold water flow of pipeline.
3. the garden integrated energy system Optimization Scheduling according to claim 1 for considering hot pipe network transmission characteristic,
It is characterized in that, the concentration energy source station-pipe network coupling constraint in the S2 specifically uses following formula to indicate:
Tces,min≤Tt ces,s≤Tces,max
In above formula,Energy source station is concentrated to energize power for t moment；Tt ces,r、Tt ces,sRespectively t moment concentrates energy source station air-conditioning
Cold-water return temperature and concentration energy source station air-conditioning supply water temperature；The water supply being respectively connected with concentration energy source station
The temperature of pipeline headend node and the temperature of water return pipeline endpoint node；
Tces,min、Tces,maxRespectively concentrate energy source station outlet air conditioner cold water lowest temperature and the upper limit；ΩcesFor with concentrate energy source station
The set for the pipeline being connected.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201811353627.XA CN109472413A (en) | 2018-11-14 | 2018-11-14 | Optimal scheduling method of park comprehensive energy system considering heat pipe network transmission characteristics |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201811353627.XA CN109472413A (en) | 2018-11-14 | 2018-11-14 | Optimal scheduling method of park comprehensive energy system considering heat pipe network transmission characteristics |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN109472413A true CN109472413A (en) | 2019-03-15 |
Family
ID=65672632
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201811353627.XA Pending CN109472413A (en) | 2018-11-14 | 2018-11-14 | Optimal scheduling method of park comprehensive energy system considering heat pipe network transmission characteristics |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN109472413A (en) |
Cited By (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109523137A (en) * | 2018-10-29 | 2019-03-26 | 天津大学 | Consider the garden comprehensive energy Optimization Scheduling of building thermal load demands response |
| CN111478326A (en) * | 2020-05-12 | 2020-07-31 | 南方电网科学研究院有限责任公司 | Comprehensive energy optimization method and device based on model-free reinforcement learning |
| CN111611690A (en) * | 2020-04-17 | 2020-09-01 | 清华大学 | Dynamic calculation method for operating parameters of heat pipe network in comprehensive energy network |
| CN112819191A (en) * | 2019-11-18 | 2021-05-18 | 国网能源研究院有限公司 | Garden comprehensive energy system operation optimization analysis method oriented to engineering application |
| CN113108353A (en) * | 2021-04-26 | 2021-07-13 | 杭州鸿晟电力设计咨询有限公司 | Variable flow variable temperature zone heating system scheduling method considering heat loss nonlinearity |
| CN114126365A (en) * | 2021-11-24 | 2022-03-01 | 华中科技大学 | High-energy-efficiency liquid cooling method and liquid cooling system for data center |
Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN105869075A (en) * | 2016-04-19 | 2016-08-17 | 东南大学 | Economic optimization scheduling method for cold, heat and electricity combined supply type miniature energy grid |
| CN106056251A (en) * | 2016-06-12 | 2016-10-26 | 清华大学 | Electric-thermal coupled multi-energy-flow system optimization scheduling method |
| CN106339772A (en) * | 2016-08-11 | 2017-01-18 | 清华大学 | Heat supply pipe network thermal storage efficiency-based thermal-electric joint optimal dispatching method |
| CN107807523A (en) * | 2017-10-18 | 2018-03-16 | 国网天津市电力公司电力科学研究院 | Consider the Regional Energy internet multi-source coordination optimization operation reserve of tou power price |
| CN108258679A (en) * | 2017-12-25 | 2018-07-06 | 国网浙江省电力有限公司经济技术研究院 | Consider the electric-thermal integrated energy system Optimization Scheduling of heating network heat accumulation characteristic |
| CN108446809A (en) * | 2018-04-09 | 2018-08-24 | 国网河南省电力公司经济技术研究院 | A kind of regional complex energy device and network dual-layer optimization configuration method |
-
2018
- 2018-11-14 CN CN201811353627.XA patent/CN109472413A/en active Pending
Patent Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN105869075A (en) * | 2016-04-19 | 2016-08-17 | 东南大学 | Economic optimization scheduling method for cold, heat and electricity combined supply type miniature energy grid |
| CN106056251A (en) * | 2016-06-12 | 2016-10-26 | 清华大学 | Electric-thermal coupled multi-energy-flow system optimization scheduling method |
| CN106339772A (en) * | 2016-08-11 | 2017-01-18 | 清华大学 | Heat supply pipe network thermal storage efficiency-based thermal-electric joint optimal dispatching method |
| CN107807523A (en) * | 2017-10-18 | 2018-03-16 | 国网天津市电力公司电力科学研究院 | Consider the Regional Energy internet multi-source coordination optimization operation reserve of tou power price |
| CN108258679A (en) * | 2017-12-25 | 2018-07-06 | 国网浙江省电力有限公司经济技术研究院 | Consider the electric-thermal integrated energy system Optimization Scheduling of heating network heat accumulation characteristic |
| CN108446809A (en) * | 2018-04-09 | 2018-08-24 | 国网河南省电力公司经济技术研究院 | A kind of regional complex energy device and network dual-layer optimization configuration method |
Non-Patent Citations (2)
| Title |
|---|
| 董帅等: "计及网络动态特性的电-气-热综合能源系统日前优化调度", 《电力系统自动化》 * |
| 顾泽鹏等: "考虑热网约束的电热能源集成系统运行优化及其风电消纳效益分析", 《中国电机工程学报》 * |
Cited By (9)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109523137A (en) * | 2018-10-29 | 2019-03-26 | 天津大学 | Consider the garden comprehensive energy Optimization Scheduling of building thermal load demands response |
| CN109523137B (en) * | 2018-10-29 | 2022-11-04 | 天津大学 | Garden comprehensive energy optimization scheduling method considering building thermal load demand response |
| CN112819191A (en) * | 2019-11-18 | 2021-05-18 | 国网能源研究院有限公司 | Garden comprehensive energy system operation optimization analysis method oriented to engineering application |
| CN111611690A (en) * | 2020-04-17 | 2020-09-01 | 清华大学 | Dynamic calculation method for operating parameters of heat pipe network in comprehensive energy network |
| CN111611690B (en) * | 2020-04-17 | 2021-06-22 | 清华大学 | Dynamic calculation method for operating parameters of heat pipe network in comprehensive energy network |
| CN111478326A (en) * | 2020-05-12 | 2020-07-31 | 南方电网科学研究院有限责任公司 | Comprehensive energy optimization method and device based on model-free reinforcement learning |
| CN111478326B (en) * | 2020-05-12 | 2021-09-03 | 南方电网科学研究院有限责任公司 | Comprehensive energy optimization method and device based on model-free reinforcement learning |
| CN113108353A (en) * | 2021-04-26 | 2021-07-13 | 杭州鸿晟电力设计咨询有限公司 | Variable flow variable temperature zone heating system scheduling method considering heat loss nonlinearity |
| CN114126365A (en) * | 2021-11-24 | 2022-03-01 | 华中科技大学 | High-energy-efficiency liquid cooling method and liquid cooling system for data center |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN109472413A (en) | Optimal scheduling method of park comprehensive energy system considering heat pipe network transmission characteristics | |
| CN111950809B (en) | Master-slave game-based hierarchical and partitioned optimized operation method for comprehensive energy system | |
| CN109523137B (en) | Garden comprehensive energy optimization scheduling method considering building thermal load demand response | |
| CN108932560A (en) | Garden integrated energy system Optimization Scheduling based on Model Predictive Control | |
| CN106960282A (en) | A kind of coordinated operation method of pneumoelectric comprehensive energy distribution network systems | |
| CN106447529A (en) | Distributed energy system modeling and running optimization method considering hot water pipe network | |
| CN106055773B (en) | In conjunction with the method for building up of the multizone integrated energy system allocation models of heat supply network model | |
| CN201297727Y (en) | Cold-storage heat-storage energy-saving synchronization system | |
| CN105241117A (en) | Data center CCHP energy supply system and method utilizing cold source | |
| CN112182887B (en) | Comprehensive energy system planning optimization simulation method | |
| CN108830743A (en) | Consider the garden integrated energy system Optimization Scheduling of a variety of cold-storage devices | |
| CN107358345A (en) | The distributed triple-generation system optimizing operation method of meter and dsm | |
| CN204827549U (en) | System is utilized to high efficiency that pipe network natural gas pressure energy and cold energy were retrieved | |
| CN112434402A (en) | Interval practical safety domain modeling method | |
| CN209068590U (en) | A kind of multi-energy complementation utilization system | |
| CN105226649A (en) | A kind of provincial power network power generation dispatching optimization method improved based on bus load prediction | |
| CN116911437A (en) | Data center energy system optimal scheduling method considering load response characteristics | |
| Mao et al. | Study on Multi-energy Flow Management for Multi-energy Network Systems. | |
| CN203432006U (en) | Cold storage and combined cold supplying type air-conditioning system | |
| CN106500397A (en) | A kind of cold, heat and power triple supply system | |
| Wang et al. | A district energy planning method with mutual interconnection and interchange of thermal grids | |
| CN111461429B (en) | Multi-energy complementary system optimization method for gas-electricity coordination | |
| CN110397979B (en) | Cold, heat and electricity triple supply new energy storage energy supply and peak shaving system | |
| Wirtz et al. | Bidirectional low temperature networks in urban districts: A novel design methodology based on mathematical optimization | |
| CN114186433A (en) | Optimized scheduling method, device, equipment and medium for comprehensive energy system |
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 | ||
| RJ01 | Rejection of invention patent application after publication | ||
| RJ01 | Rejection of invention patent application after publication |
Application publication date: 20190315 |