CN110110897A - A kind of integrated energy system optimization method considering different storage energy operation strategies - Google Patents
A kind of integrated energy system optimization method considering different storage energy operation strategies Download PDFInfo
- Publication number
- CN110110897A CN110110897A CN201910289594.5A CN201910289594A CN110110897A CN 110110897 A CN110110897 A CN 110110897A CN 201910289594 A CN201910289594 A CN 201910289594A CN 110110897 A CN110110897 A CN 110110897A
- Authority
- CN
- China
- Prior art keywords
- energy
- power
- integrated energy
- integrated
- electricity
- 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
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/0637—Strategic management or analysis, e.g. setting a goal or target of an organisation; Planning actions based on goals; Analysis or evaluation of effectiveness of goals
- G06Q10/06375—Prediction of business process outcome or impact based on a proposed change
-
- 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—Systems or methods specially adapted for specific business sectors, e.g. utilities or tourism
- G06Q50/06—Electricity, gas or water supply
Abstract
The invention discloses a kind of integrated energy system optimization methods for considering different storage energy operation strategies, comprising the following steps: S1, establishes a center of energy model；S2, the energy flow equation in center of energy is respectively obtained by center of energy model, energy flow equation includes power system equations, therrmodynamic system equation, cold supply system equation；S3, the balanced load operation reserve for formulating battery energy storage system respectively and combination peak valley operation reserve；S4, the cold supply system in integrated energy system, heating system, power supply system are planned respectively；S5, the Optimized model for establishing integrated energy system, Optimized model include annual total cost and constraint condition；S6, it is calculated and is solved using Optimized model of the genetic algorithm to integrated energy system.
Description
Technical field
The present invention relates to energy source optimization field more particularly to a kind of integrated energy systems for considering different storage energy operation strategies
Optimization method.
Background technique
Integrated energy system (IES) is a kind of multifunctional system, equipped with cold, heat, electricity and gas.In plan model, it is necessary to
Power grid, heat supply network, the model of combustion gas network and power constraint are considered simultaneously.It is asked individually however, most of currently existing schemes are all paid close attention to
Topic, there is no include solving the problems, such as a whole set of scheme.In addition, energy-storage battery plays an important role in IES, many rule of IES
The problem of drawing is to be planned under the premise of not considering energy storage based on certain boundaries；In addition, energy-storage system is that price is higher, storage
Energy battery all directly affects the preferred plan result and economic cost of IES for the selection of type and operation reserve, how to select to store up
Energy type and operation reserve are still urgent problem to be solved in IES optimization planning.
Summary of the invention
Object of the present invention is in view of the above-mentioned problems, providing a kind of comprehensive energy system based on different energy-storage system operation reserves
System optimization method.
To achieve the goals above, the technical scheme is that
A kind of integrated energy system optimization method considering different storage energy operation strategies, comprising the following steps:
S1, a center of energy model is established；
S2, the energy flow equation in center of energy is respectively obtained by center of energy model, energy flow equation includes electric power
System equation, therrmodynamic system equation, cold supply system equation；
S3, the balanced load operation reserve for formulating battery energy storage system respectively and combination peak valley operation reserve；
S4, the cold supply system in integrated energy system, heating system, power supply system are planned respectively；
S5, the Optimized model for establishing integrated energy system, Optimized model include annual total cost and constraint condition；
S6, it is calculated and is solved using Optimized model of the genetic algorithm to integrated energy system.
Further, in the step S1 center of energy model formula are as follows:
Wherein,Respectively represent electric load, thermic load, refrigeration duty；Pb、GbWith HbIt is energy respectively
Electric input quantity, the gas input quantity, heat input at center；C is the coefficient of coupling matrix, indicates electricity, gas, the cold coupling between heat
Relationship.
Further, power system equations in the step S2 are as follows:
Wherein, Pi pguIt is the output power of PGU, i ∈ b indicates that c-th of PGU unit is located at b-th of center of energy；
The therrmodynamic system equation are as follows:
Wherein,It is the heating output power of PGU；It is the power output situation of assist type electric boiler；frate(Pi pgu)
It is PGU exhaust mass flow；ftemp(Pi pgu) it is exhaust temperature of charge；ηauxIt is the transfer efficiency of donkey boiler；It is auxiliary
The fuel gas buring rate of boiler；
The cold supply system equation are as follows:
Wherein,WithIt is the output power of Absorption Refrigerator and electrical chillers；It is absorption refrigeration function
The consumption in source；copecIt is the Energy Efficiency Ratio of Absorption Refrigerator.
Further, in the step S3 battery energy storage system balanced load operation reserve are as follows: no matter electricity price be in use
Electric peak value, low ebb or flat section, when the output power of integrated energy system is greater than workload demand, energy-storage battery charges,
And extra electricity is sold to power grid after energy-storage battery charging；When the power output of integrated energy system is less than loading demand
When, energy-storage battery electric discharge；When energy-storage battery, which discharges into minimum and system, is still unable to satisfy demand, from power grid buy power with
Meet workload demand；
The combination peak valley operation reserve of the battery energy storage system are as follows: when the output power of integrated energy system is greater than load
Demand and when electricity price is in electricity consumption peak value, extra electricity sells to power grid；When electricity price is in level section, energy storage electricity
Half electric power is filled in pond；When electricity price is in low power consumption, energy-storage battery is fully charged；When the output power of integrated energy system
No more than workload demand and when electricity price is in electricity consumption peak value, energy-storage battery buys power after discharging into minimum state from power grid；When
When electricity price is in level section, buy power after energy-storage battery release half electric power from power grid；When to be in electricity consumption low for electricity price
Gu Shi, energy-storage battery do not discharge, and directly buy power from power grid.
Further, the cold supply system in integrated energy system is planned in the step S4, integrated energy system
Refrigeration duty by Absorption Refrigerator and electric refrigerating machine work complete, formula are as follows:
When the refrigeration duty of integrated energy system is less than the power output of Absorption Refrigerator, then absorption refrigeration is used completely
Machine carries out on-the-job satisfaction refrigeration duty, and electric refrigerating machine does not work, formula are as follows:
When the output power of Absorption Refrigerator is insufficient for the refrigeration duty of integrated energy system, Absorption Refrigerator
It is worked with rated power, electric refrigerating machine makes up remaining refrigeration duty, formula are as follows:
Further, the heating system in integrated energy system is planned in the step S4, integrated energy system
Thermic load by generating set and auxiliary electric boiler complete, formula are as follows:
When the thermic load of integrated energy system is less than the output power of generating set, then work is carried out by generating set completely
Make, auxiliary electric boiler does not work, formula are as follows:
When the output power of generating set is insufficient for the thermic load of integrated energy system, generating set will be with specified
Power work, assist type electric boiler make up remaining thermic load, formula are as follows:
Further, the power supply system in integrated energy system is planned in the step S4, integrated energy system
Power supply by photovoltaic panel, wind turbine, generating set, energy-storage battery, power grid complete, formula are as follows:
Wherein, Pi pv(t) be photovoltaic panel output electricity, Pi wt(t) be wind turbine output electricity；Pi pguIt is power generation
The output electricity of unit；Pi BESS(t) be energy-storage battery output electricity；Pgrid(t) power grid exports electricity.
Further, in the step S5 annual total cost calculation formula are as follows:
F=min ATC=Caic+Co&m+Cre+Ce+Cgrid；
Wherein, ATC is integrated energy system year minimum cost；CaicIt is the equipment investment cost of integrated energy system；Co&m
The operation expense of integrated energy system；CreIt is the replacement cost of integrated energy system；CeThe energy of integrated energy system at
This, CgridIt is the power grid interaction cost of integrated energy system；
Wherein, P and CpIt is the installed capacity and the specific investment cost cost of capacity of equipment；FcrFor capital recovery factor:
Wherein, Co&m,jIt is the unit operation and maintenance cost of equipment j；Wd,t,jIt is the production cost of equipment j；
Wherein, CrjIt is the replacement cost of equipment j, FrpIt is capital recovery factor, a is the service life of equipment, year；Y is comprehensive
Close the service life of energy resource system；
Wherein, cfuelIt is the energy cost of gas consumption；Fd,tThe gas consumption of integrated energy system；cgrid-PIt is
Electric energy purchasing price；Ed,t,pElectricity is bought from power grid；cgrid-sSale of electricity price；Ed,t,sIt is that integrated energy system to power grid sells electricity
Amount.
Further, constraint condition includes power constraint, natural gas grid constraint, heat supply network constraint, energy storage in the step S5
Battery constraint；Power constraint includes power balance constraint and the constraint of electric power stream；Natural gas grid constraint includes natural gas node
Flow equilibrium constraint and gas discharge-pressure confines；Heat supply network constraint includes thermodynamic branch characteristic equation, thermodynamics network
Equilibrium equation and equipment export-restriction；Energy-storage battery constraint includes the maximum constrained and energy-storage battery capacity of energy-storage battery capacity
Least commitment.
Further, when in the step S6 to the seismic responses calculated of integrated energy system, variable are as follows:
Wherein, 0,1 selects for the operation reserve of battery energy storage system, NGFor gas ductwork, NEFor power grid.
Compared with prior art, the advantages and positive effects of the present invention are:
The invention proposes two kinds of operation reserves of battery energy storage system, and have been respectively compared different operation strategy to comprehensive
Close the influence of energy system planning result；Its operation reserve based on each equipment, it is contemplated that power grid, heat supply network and combustion gas network and its
Related constraint establishes the optimum programming model of integrated energy system, with the annual total cost of integrated energy system for optimal mesh
Mark, solves the model using genetic algorithm, has obtained corresponding optimal solution, more realistically reflect in comprehensive energy physical planning
Consider the problems of, improves the economic effect of integrated energy system.
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 without any creative labor, may be used also for those of ordinary skill in the art
To obtain other drawings based on these drawings.
Fig. 1 is integrated energy system schematic diagram；
Fig. 2 is coupling process flow chart；
Fig. 3 is the first operation reserve figure of battery energy storage system；
Fig. 4 is second of operation reserve figure of battery energy storage system；
Fig. 5 is genetic algorithm flow chart；
Fig. 6 is integrated energy system network topology structure figure；
Fig. 7 is the power output figure of WT/PV；
Fig. 8 is multiple load chart；
Fig. 9 is the fitness curve graph of GA and PSO；
Figure 10 is the SOC state diagram one of battery energy storage system；
Figure 11 is distribution and official website result schematic diagram one；
Figure 12 is the SOC state diagram two of battery energy storage system；
Figure 13 is distribution and official website result schematic diagram two；
Figure 14 is battery price sensitivity analysis curve；
Figure 15 is load diagram hourly；
Figure 16 is the power output result figure of electric energy per hour；
The SOC state diagram of Figure 17 Different Strategies；
Figure 18 is the scheduling result figure of hot and cold per hour.
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, any modification, equivalent replacement, improvement and so on should all be included in the protection scope of the present invention.
As shown in Fig. 1 to Figure 18, integrated energy system coupled structure is illustrated first；
The cold-hot of integrated energy system-electric-gas structure is by central cooling system, heating system, electric system, natural gas
System and coupling element composition.And these systems and coupling element are distributed in the same area, as shown in Figure 1.It is coupled in multipotency
In mode, center of energy shown in FIG. 1 is integrated with energy conversion, distributes as its basic model unit, storage, the function such as transmission
Energy.It represents energy participant (supplier and consumer) and electric power, the energy interaction between combustion gas and therrmodynamic system.It is
The most basic and key link of high efficient and flexible energy supply is realized in integrated energy system.
In the present invention, generating set (PGU), wind turbine, photovoltaic and power grid are for meeting region electric load
Demand.Heating system is by heat source, heat supply network and thermic load composition.Thermal power is mentioned by PGU, gas fired-boiler (GB), electric boiler (EB) etc.
For.The structure of cold supply system is similar, and cooling power is converted by Absorption Refrigerator (AC) and electric refrigerating machine (EC).In IES
Natural gas system be mainly low-voltage network (7.5*103Pa and following), by gas source, the gentle load composition of steam line.Pressure
Draught control mechanism is typically mounted at gas source, and compressor and pressure regulator consider [30] not in piping network.In short, due to drawing
A series of coupling elements are entered, more energy streams have close coupling relationship, so needing to carry out Conjoint Analysis.
According to the coupled relation of multiple-energy-source, the present invention constructs center of energy (EH) model to describe its energy coupling
Mode.Energy conversion, the function of storage and distribution [31] may be implemented in center of energy, and the model can describe cold-hot-
Coupled relation between electric-gas.Concrete model is as follows.
Respectively represent electric load, thermic load, refrigeration duty.Pb,GbWith HbIt is the electricity of center of energy,
Gas and heat input.C is the coefficient of coupling matrix, indicates electricity, gas, the cold coupled relation between heat.
According to above-mentioned center of energy model, the energy flow equation in center of energy can pass through matrix-expand, such as equation
(2)-(9) abstraction matrix can be converted to particular equation and be added in model and as constraint.
(1) power subsystem
The electric output power of PGU is directly related with its fuel gas buring rate, and present invention introduces formula (3) expression.
Pi pguIt is the heat power output of PGU, i ∈ b indicates that c-th of PGU unit is located at b-th of center of energy.
(2) heating power subsystem
Thermic load is mainly met by the thermal output of PGU and the donkey boiler of b-th of center of energy.
It is the hot power output situation of PGU；It is the power output situation of assist type electric boiler；frate(Pi pgu) it is PGU exhaust
Mass flow .ftemp(Pi pgu) it is exhaust temperature of charge ηauxIt is the transfer efficiency of donkey boiler.It is the combustion gas of donkey boiler
Combustion rate.
(3) for refrigeration subsystem
Refrigeration duty is mainly met by AC and EC in the output of b-th of center of energy.
WithIt is the power output of Absorption Refrigerator and electrical chillers；It is the consumption of the Absorption Refrigerator energy
Amount；copecIt is the Energy Efficiency Ratio of Absorption Refrigerator.
In integrated energy system research of the invention, meets want to meeting cooling load of the air-conditioning system for refrigeration subsystem first
It asks.By the heat transfer needed for freezing, to generating set, and when being unable to satisfy refrigeration duty, starting electricity freezes Absorption Refrigerator
Machine.Secondly, thermic load subsystem meets the needs of thermic load.PGU works with FTL mode operation, and PGU needs to meet cooling supply and confession
Thermic load, and start electric refrigerating machine when being unable to satisfy refrigeration duty.Finally, electric energy subsystem is used to meet the electric load in system
With the needs of electrical power needed for cooling supply and heat supply.BESS plays the role of adjusting electric power output difference, and main power grid is comprehensive energy
The backup power source of source system compares the electric power and system loading of PGU, blower WT and photovoltaic PV generation based on this condition.Such as Fig. 2 institute
Show:
Obviously, the different operation reserves of IES lead to different system design schemes.In addition, suitable operation reserve is to being
The performance of system is critically important.Invention describes two kinds of different operation reserves of BESS, have obtained corresponding system configuration result.This
Invention selects different types of battery to compare, and analyzes different battery behaviors to the planning knot of BESS system difference operation reserve
The influence of fruit.
The first operation reserve of BESS is balanced load demand (as shown in Figure 3).No matter electricity price be in peak value, valley or
Person's usually electricity price level, when the output power of integrated energy system is greater than workload demand, BESS charging, and in energy-storage battery
Extra electricity is sold to power grid after charging.When the power output of integrated energy system is less than loading demand, BESS electric discharge.When
When BESS discharges into minimum and system and is still unable to satisfy demand, buy power from power grid to meet workload demand.
Second of configuration strategy of BESS is that workload demand is combined (as shown in Figure 4) with peak valley price.The of BESS
Two kinds of configuration modes are to carry out balanced load demand in conjunction with time-of-use tariffs.In this mode, it is necessary to judge whether the time is in use
Energy peak, flat sections or low ebb moment, and according to circumstances execute charge or discharge.
If Ppv(t)+Pwt(t)+Ppgu_e(t)-Peload> 0, when it is in electricity consumption peak value and BESS is not responded to,
Extra electricity sells to power grid；When it is in energy flush end, BESS charges part；When it is in energy low-valley interval
When, BESS executes charging.
If Ppv(t)+Pwt(t)+Ppgu_e(t)-Peload≤ 0, when it is in peak of power consumption, BESS discharges into minimum shape
State；When it is in the flat section of energy, BESS discharges some electrical power；When with the energy period at a low ebb, BESS will not discharge, directly
Insufficient partial power is bought from power grid.When above situation is unable to meet demand, bought power from power grid to meet load need
It asks.
Cold supply system planning
As described above, the refrigeration duties of IES are completed by Absorption Refrigerator and electric refrigerating machine:
Think Absorption Refrigerator better than electric refrigerating machine in refrigeration work.If refrigeration duty is less than Absorption Refrigerator
Calibration power output, then carried out meeting refrigeration duty completely with Absorption Refrigerator.Electric refrigerating machine does not work:
When the output power of Absorption Refrigerator is insufficient for refrigeration duty, it will will more with rated power operation, EC
Tonifying for the deficiency, EC operation carry out cooling supply:
After whole year is completed in simulation, the standard installed capacity of EC can be found:
Heating system planning
The thermal load demands of integrated energy system are as follows: thinking that AC keeps constant power output under different loads.
The heat supply work of integrated energy system is completed by PGU and auxiliary electric boiler.
Ph(t) be integrated energy system normal thermic load；copacIt is the transfer efficiency of Absorption Refrigerator.
If thermic load is less than the rated output of PGU, PGU is fully met, and auxiliary electric boiler does not work:
When the output power of PGU is insufficient for thermic load, it will will more with rated power operation, assist type electric boiler
Concurrent heating load difference:
After the completion of annual running simulation, the normal capacity of assist type electric boiler can be found:
Power supply system planning
The electricity needs of integrated energy system includes electric power needed for electric load and operation EC and EB.
The powered operation of IES is completed by PV, WT, PGU, BESS and power grid.
Pi pv(t) and Pi wt(t) be photovoltaic and blower output electricity；Pi pguIt is the output electricity of PGU；Pi BESSIt (t) is storage
The discharge capacity of energy battery；Pgrid(t) power grid exports electricity
In the present invention, Optimized model is established, plans integrated energy system based on year minimum total cost (ATC)
(IES).The life cycle of integrated energy system is 20 years.And the annual total cost of integrated energy system generally comprises the throwing of equipment
Provide cost, operation expense, replacement cost, power grid interaction cost and energy cost.Cost of investment mainly includes wind-power electricity generation
Machine, photovoltaic power generation equipment, BESS equipment, the acquisition cost of the equipment such as gas turbine.Operation expense is mainly equipment operation
Cost.Replacement cost is primarily referred to as the replacement cost of distinct device.It mainly includes electric power income from sales and electricity that power grid, which interacts cost,
Power buying expenses.Energy cost is the cost of gas consumption.Objective function is defined as:
F=min ATC=Caic+Co&m+Cre+Ce+Cgrid (25)
ATC is integrated energy system year minimum cost；CaicIt is the cost of integrated energy system；Co&mIntegrated energy system
Operation expense；CreIt is the replacement cost of comprehensive energy；CeIntegrated energy system buys gas cost, CgridIt is comprehensive energy
Source system interacts cost with power grid.
N and C is the initial outlay cost of some installation installation costs and equipment, respectively unit cost；P and CpIt is it
The installed capacity of his equipment and the specific investment cost cost of capacity；FcrIt can be defined as capital recovery factor:
The operation of different units and maintenance cost are always to produce relevant unit cost according to equipment to calculate.Co&m,jIt is to set
The unit operation of standby j and maintenance cost；Wd,t,jIt is the production cost of equipment j.
The replacement cost of distinct device depends on the service life of fund cost and equipment, CrjIt is the replacement cost of j, Frp
It is capital recovery factor, a is the service life of equipment, year；Y is the service life of integrated energy system, year.
cfuelIt is the energy cost of gas consumption；Fd,tThe gas consumption of integrated energy system；cgrid-PIt is electric energy
Purchasing price；Ed,t,pElectricity is bought from main power grid；cgrid-sSale of electricity price；Ed,t,sIt is that integrated energy system to main power grid sells electricity
Amount.
Constraint condition
(1) power grid limits
Power constraint includes: that power balance constraint (33)-(36) and electric power networks can flow constraint (37)-(38).
1) power-balance constraint
The constraint of power-balance mainly includes the gross electric capacity for being input to b-th of center of energy, start node and terminal section
Transimission power between point is poor and unbalanced electric energy.Their physical relationship is as follows:
U represents the polymerization of existing standard element；It is the electric power input of b-th of energy centre；It is line transmission function
Rate；R (l) and k (l) is respectively the start node of terminal node and route；It is the output power of the unit of power grid；For
Insufficient electric power；It is extra electric power；Power unbalance amount.If uneven electric flux is positive, show node electricity
It can be insufficient.If be negative, then it represents that node has remaining electric energy；Pi maxWith Pi minRespectively represent the upper limit of conventional elements output
And lower limit.
2) power flow constrains
Power flow constraint in power grid includes transimission power, circuit admittance etc..Their physical relationship is as follows:
The polymerization of Where L expression transmission line；BlAccess electricity；θ is the phase angle of two end node of line；M number is very big, specifically
Numerical value can refer to document [33]；Pl maxIt is the upper limit of line transmission power；It is the mode of operation of transmission line.
(2) limitation of natural gas network
In order to ensure the safe operation of gas distributing system, constraint includes that node flow balances, node pressure limitation and pipeline
Flow restriction.
1) the flow equilibrium constraint of natural gas node
The flow equilibrium constraint of natural gas node mainly includes the flow of i-th of injection network, upstream and downstream node
The natural gas load etc. of flow and i-th of node.Their physical relationship is as follows:
Fl,g,s,t,iIt is the process of the node i of injection network；J, k are the upstream and downstream nodes of node I；；Fg,s,t,j,
Fg,s,t,kIt is the flow of j, k node respectively；F2,g,s,t,iIt is the natural gas load of node I；It is the transmission flow pole of pipeline ij
Value.
2) gas discharge-pressure confines
Gas flow pressure constraint mainly includes gas flow direction, gas pressure amplitude, feed flow extreme value, gas function
Rate and gas heating value.Relationship is as follows:
Wherein, Fg,s,t,ijIt is the flow ij of pipeline；γs,t,ijIt is a symbolic variable, represents fuel gas flow direction；KijIt is
The transmission coefficient of pipeline；Ug,s,t,iIt is the pressure amplitude of node i；Ug,s,t,jIt is the pressure amplitude of node j；andIt is section
The upper and lower bound of the pressure amplitude of point i.
(3) heat supply network limits
1) thermodynamic branch characteristic equation
For pipeline k, inflow heat flux and outflow heat flux at time t are defined as follows:
Consider the transmission delay τ of the pipeline K of t momentk,t,WithRelationship it is as follows:
τk,t=Fk/Gk,t (47)
Gk,tIt is the flow of t moment pipeline K；C is refrigerant volume specific heat；WithBe respectively pipeline K entrance and
Outlet temperature；εkIt is the loss coefficient of pipeline K；FkIt is the characteristic quantity of pipeline k, it is determined by duct length, the parameter of area of section
It is fixed.
2) equilibrium equation of thermodynamics network
The equilibrium equation of thermodynamics network is as follows:
In formula, A1、A2It is starting and the termination incidence matrix of thermodynamics network respectively.It is each in time period t
The vector value of the inflow heat flux of pipeline.It is the vector value of the outflow heat flux of each pipeline in time period t.QtBe by when
Between in section t the hot-fluid inbound traffics composition of all nodes vector value, wherein heat source node inflow is positive, and the inflow of heat load node is
It is negative.
3) equipment export-restriction
Production of energy and energy transition equipment in integrated energy system specifically include that WT, PV, CHP, GB, ER and AR etc..
Its output constraint is as follows:
0≤Ppv≤Ppv_max (40)
0≤Pwt≤Pr (41)
Ppgu_min≤Ppgu≤Ppgu_max (42)
Peb_min≤Peb≤Peb_max (43)
Pec_min≤Pec≤Pec_max (44)
Pac_min≤Pac≤Pac_max (45)
Ppv_maxIt is the output upper limit of PV, PrIt is the output upper limit of WT, Ppgu_maxAnd Ppgu_minIt is the output upper limit of PGU under
Limit, Peb_maxAnd Peb_minIt is the output upper and lower bound of EB, Pec_maxAnd Pec_minIt is the output upper and lower bound of EC, Pac_maxWith
Pac_minIt is the bound output of AC.
(4) energy-storage battery limits
In the limitation of battery, the main maximum charge for considering battery and discharge power and SOC are constrained.
SOCmin≤SOC(t)≤SOCmax (46)
Wherein, SOCmin、SOCmaxThe minimum and maximum constraint of energy-storage battery capacity, Pch,Pdis,minIt is to be charged and discharged function
The minimum value of rate, Pch,Pdis,maxIt is the maximum value for being charged and discharged power.
The invention proposes for economic IES Optimized model.This is a typical nonlinear optimal problem, is related to being permitted
Multi-non-linear, discrete, random and uncertain factor.Formula (1)-(47) in Section 2 and Section 3 constitute integrated energy system
Coordinated planning Optimized model.The model includes many nonlinear components, such as Load flow calculation (37) and (38), and considers to use
In the calculating (41) of the combustion gas network pressure-flow relationships of traditional coordinated planning be non-linear variable.In addition, system equipment will be defeated
The scheduling factor that energy distributes to other equipment out makes plan model be nonlinear model.This is the non-linear more of CCHP system
The case where energy coupling, such as equation (4)-(9).The general type of MINLP mixed integer nonlinear programming problem is as follows.
Wherein, x, y are optimized variables；fn(x, y) is optimization object function；X is one group of variable；Y is all possible solutions
The set of scheme；J (x, y) is one group of inequality constraint；H (x, y) is one group of equality constraint.
There are two types of main Types for mixed integer nonlinear programming solution technology: mathematic programming methods and objective optimization side
Method (such as GA algorithm, PSO algorithm).The coupling being related between multiple energy and network constraint by the plan optimization problem proposed
Relationship, and many nonlinear equations and variable (including continuous variable and discrete variable) are also contemplated, therefore, it is difficult to pass through number
Programmed method is learned to solve the problems, such as this.Moreover, only locally optimal solution can be obtained by mathematics programmed method.It is mentioned in the present invention
In Optimized model out, primary variables is optimization time, place capacity and device type and pipe-line.A series of variables exist
In formula (49).
Wherein 0,1 variable represents the type selection of BESS, NG,NEIt is gas ductwork and power grid.
According to the Optimized model of integrated energy system, which is solved using genetic algorithm (GA).Genetic algorithm (GA) is
A kind of widely used heuristic value based on natural selection in wheel disc strategy.Since genetic algorithm has, precision is high,
Fast convergence rate is easy to the advantages that applying, and can efficiently solve mixing nonlinear optimal problem of the invention, prestissimo
Obtain globally optimal solution.There are four basic program, i.e. initial population, selections to intersect and make a variation for genetic algorithm.The stream of the GA of proposition
Journey figure is as shown in Figure 5.
Fig. 5 gives the optimization process of the IES system proposed in genetic algorithm (GA).Firstly, initial parameter is defined as follows:
Workload demand, key operating parameters and genetic algorithm parameter.Workload demand is primarily referred to as the refrigeration duty of IES, and thermic load and electricity are negative
Lotus.Crucial operating parameter is energy and economic aspect in IES.Genetic algorithm parameter mainly includes Population Size, is evolved generation
Number, select probability Ps, crossover probability PcWith mutation probability Pm.Next, initializaing variable is with 0/1 string encoding.In addition, hereditary
Random general sampling is respectively adopted in operator, and single point crossing and discrete variation determine Ps,PcWith Pm.Finally, when convergence satisfaction is wanted
When asking, stop search process, and determines optimum results by decoding.If being unsatisfactory for convergence, search process will continue,
Until obtaining optimum.
Analysis of cases
In order to verify proposed optimization method, using by 18 node power distribution networks and 11 node low pressure natural gas group of networks
At IES emulated and analyzed.Fig. 6 shows the topological structure of IES network.Distribution network is 20kV, shares 16 loads
Node, 2 power transformation tiny nodes and 24 planning routes.17 and 18 nodes are power transformation tiny nodes, other nodes are load nodes.It
There are 9 load nodes and 2 natural gas station nodes in right gas network.In addition, natural gas network includes 15 planning pipelines and one
A gas source (GS) center.10 and 11 nodes are natural gas station nodes, other nodes are load nodes.In energy coupling part
In, energy coupling is mainly executed by energy hub.Energy hub chief component be case study part CCHP and
Heat exchanger (HE).Since selected item has had complete ther mal network, ther mal network is not executed most in this case
Good planning.
In this case, time interval is one hour, therefore in planning process, there is 8760 points in 1 year.It is selected
Planned project include fixed photovoltaic and wind turbine, therefore the not capacity of re-optimization photovoltaic and wind turbine.Photovoltaic and
The initial outlay of wind turbine is not included in annual total cost.In the best configuration of integrated energy system, two kinds of situations are considered
It is analyzed.Three kinds of different types of batteries are selected, i.e. lithium battery, lead-lead carbon battery and sodium-sulfur battery is compared, and is analyzed
Influence of the different battery behaviors to the program results under two kinds of operation strategies of BESS.
In the present invention, renewable energy such as wind speed, radiation intensity are the important input data that China Meteorological Administration obtains, shadow
Ring the power output of WT and PV plate.The power output of WT/PV is annual 8760 points, as shown in Figure 7.In addition, the energy of the IES project
Source demand includes refrigeration duty, thermic load and electric load, and multiple load is as shown in Figure 8.
In the planning and operation of integrated energy system, electricity price and gas price directly affect the result of plan optimization.Cause
This, must clarify the electricity price and gas price of project before planning.In present case research, implement the day of comprehensive energy project
Right gas price lattice fixed value is 0.025 $/m3, and the electricity price of execution is the Peak-valley TOU power price of purchase and sale electric power.IES project
Peak-valley TOU power price is shown in Table 2.
The Peak-valley TOU power price of 2 IES project of table
In the planning process of present case, used equipment and relevant parameter mainly include initial outlay cost, operation
Maintenance cost, life cycle, efficiency, the economic technologies parameter such as rated power.In the best configuration of BESS, three types are selected
Battery be compared, including lithium battery, lead carbon battery and sodium-sulphur battery.The charge-discharge characteristic of these three batteries is different, technology
Parameter also has very big difference, and battery life is also influenced by depth of discharge.Design parameter is shown in Table 3 and table 4.
Each device parameter of table 3
The parameter of 4 BESS of table.
In above content of the present invention, derivation algorithm is described, and some typical algorithms can be used to solve to integrate
The configuration result of energy resource system.In order to study the feature of genetic algorithm, the model of foundation is respectively by genetic algorithm and differential evolution
Algorithm (DE) calculates.Obtained fitness Dependence Results are as shown in Figure 9.It is the optimization mould with lithium battery in the upper part of Fig. 9
The iterativecurve of type.It is the iterativecurve of the Optimized model with lead carbon battery in the middle section of Fig. 9.In the lower half of Fig. 9
Point, it is the iterativecurve using the Optimized model of sodium-sulphur battery.
From the upper part of Fig. 9 it can be clearly seen that when the computation model in GA, found after 15 iteration optimal
Solution, amounting to evaluation time is about 167.3734s, and the calculating time of optimal solution is 25.11s.However, PSO algorithm is in the 28th iteration
Decline into locally optimal solution, amounting to evaluation time is about 217.1427s.Similarly, in the middle section of Fig. 9 and lower part, GA is found
The number of iterations of optimal solution and time are less than number and the time of PSO algorithm.
Obviously, it determines in scheme in allocation optimum solution to model proposed by the present invention, genetic algorithm has higher convergence fast
Degree and better convergence, can search for globally optimal solution to the maximum extent, more suitable for solving model proposed by the present invention.
Energy-storage battery strategy one
Under first operation reserve of BESS, the battery complete response power budget difference of three types.In this feelings
Under condition, select the variation of three types battery in one week as analysis object.The SOC state of three kinds of batteries is as shown in Figure 10.
From fig. 10 it can be seen that the SOC curvilinear trend of three kinds of batteries is roughly the same, when power output overloads demand
When, BESS charges；Otherwise, BESS discharges.But since the power of three kinds of batteries is different, so the SOC of battery is when each
Between dotted state it is different.According to the operation reserve of BESS, the capacity configuration of IES is optimized.
According to solution shown in table 5 as a result, the total annual cost of the integrated energy system using lead carbon battery
(ATC) minimum, it is 223264 dollars.Total annual cost difference between lead carbon battery and lithium battery is very small, and only 3.16%.By
It is relatively small in the initial outlay cost of lead carbon battery, but its efficiency and service life also with depth of discharge increase and significantly drop
Low, corresponding replacement cost interacts increased costs with power grid.However, with the increase of depth of discharge, the efficiency of lithium battery and longevity
Life is relatively stable, and the initial outlay cost of lithium battery is relatively large, therefore the difference of total annual cost is relatively small.Due to initially at
This and replacement cost are higher, and the total annual cost of sodium-sulphur battery is maximum.
The energy-storage battery cost of 5 IES of table
The configuration result of 6 IES of table
According to the technical characterstic of different BESS, optimum results are as shown in table 6.By comparing three kinds of configuration results, because of lead
The transfer efficiency very little of carbon battery, it can be seen that the capacity of lead carbon battery is relatively large.However, sodium-sulphur battery is due to higher first
Begin to invest and there is relatively small capacity.In addition, the capacity of PGU is opposite with the capacity of electric boiler.Meanwhile Absorption Refrigerator
It is also opposite with the capacity of electronic cooling-water machine.Electric refrigerating machine installed capacity is significantly greater than Absorption Refrigerator, shows electric refrigerating machine more
It is economical.
In the distributing rationally of lead storage battery integrated energy system, the power of distributing rationally of BESS is 45kW, BESS capacity
For 190kwh；The capacity of PGU is 305 kilowatts；The capacity of AC is 81 kilowatts；The capacity of EB is 126 kilowatts；The capacity of EC is
315kW.Best configuration result and corresponding cost are shown in Table 7.
The corresponding cost of 7 best configuration of table
Under the operational mode of BESS, the capacity of PGU is 305kW.It is negative that the electric energy that PGU is generated can provide local electric energy
Lotus, thus the electric energy for reducing the transimission power of distribution network and being bought from power grid.In addition, power distribution network, which can choose, has smaller hold
The distribution line type of loading capability and more flexible electric network composition.PGU consumes natural gas for power supply and heat supply, and passes through absorption
Formula refrigeration machine provides cooling.Therefore, corresponding gas pipeline should be constructed according to the gas consumption of gas load and PGU.Specifically
Route planning result is as shown in figure 11.
Energy-storage battery strategy two
According to second of operation reserve of BESS, the battery response electric load difference and time-of-use tariffs of three types.Selection one
The variation of three kinds of batteries is as analysis object in week.The SOC state of three kinds of batteries is as shown in figure 12.
It can be recognized from fig. 12 that the SOC curvilinear trend of three kinds of batteries is roughly the same.When electric loading difference is greater than 0, on peak
Period electric power is directly sold to power grid, and BESS will not be responded.During paddy valence, BESS electricity is collected first and is sold to remainder
Power grid.During regular price, when SOC is greater than 0.5, electric power is directly sold to power grid, and BESS is not involved in response, when SOC is less than 0.5
When, SOC is charged to equal to 0.5 to BESS.Extra electric power is sold to power grid.When electric loading difference is less than 0, BESS is in peak period
Between discharge, directly buy power from power grid during paddy valence, BESS is not involved in response.During regular price, if SOC is greater than
0.5, BESS discharges into SOC equal to 0.5.It when SOC is less than 0.5, buys power from power grid, BESS is not involved in response.According to BESS
Operation reserve, optimize IES place capacity configuration.
8 IES energy-storage battery cost of table
The configuration result of 9 IES of table
According to solution shown in table 8 and table 9 as a result, under second of operation reserve of BESS, different batteries
The annual total cost of integrated energy system has different degrees of reduction, this is because the capacity of BESS it is smaller and reduce initial outlay and
Operation maintenance cost.According to this strategy, integrated energy system and interacting for power grid are more.During the price of peak, dump power
It is distributed directly to power grid, directly supplements short electric power from power grid during paddy valence, therefore, the BESS response time is relatively fewer, BESS
Capacity it is relatively small.In view of the influence of peak valley price, although integrated energy system and interacting for power grid are more, with power grid
Interaction cost but reduces.Do not change, the year cost using the integrated energy system of lead carbon battery is still the smallest.
In the distributing rationally of lead storage battery integrated energy system, the optimization power of BESS is 39kW, and BESS capacity is
106kwh；The capacity of PGU is 290 kilowatts；AC capacity is 78 kilowatts；The capacity of EB is 135 kilowatts；The capacity of EC is 318kW.Most
Good configuration result and corresponding cost are shown in Table 10.
The corresponding cost of 10 best configuration of table
Under the operational mode of BESS, the capacity of PGU is 285kW.Compared with the capacity of the lower PGU of strategy 1, the installation of PGU
The gas consumption that capacity reduces 9.84%, PGU is reduced, and the flow of required gas pipeline is correspondingly reduced.However, PGU
Power generation is also correspondingly reduced, and compared with strategy 1, the relatively large and required electric energy of the installed capacity of electric boiler is relatively large.
Therefore, it is necessary to the electric power for increasing the transimission power of power distribution network and being bought from power grid.Detailed programme path figure is as shown in figure 13.
In order to further analyze influence of the Parameters variation to model proposed by the present invention, this section has carried out battery price variation
The sensitivity analysis of rate.This section studies influence of the different battery prices variations to ATC and battery capacity.Due in BESS strategy 2
Lower ATC is lower, therefore selection scheme is calculated and analyzed in this mode.As the result is shown in Figure 14, wherein battery cost
Change rate is the ratio of each battery cost variation and its base cost, and base cost is battery cost of the present invention.
Influence of the battery price to ATC in the case of the left part of Figure 14 shows three kinds.With the rise of battery price,
ATC is rising in all three cases.In the case where sodium-sulphur battery, ATC increases at most, and the lead carbon battery the case where
Under, ATC increases minimum.The right part of Figure 14 it can be extrapolated that with price rise, the capacity of battery will reduce accordingly.But
The rising range of battery price has been more than the decline of capacity, leads to the ascendant trend of ATC in integrated energy system, and price compared with
The up-trend of high sodium-sulphur battery is the largest, and the up-trend of the lower lead carbon battery of price is the smallest.
Influence of the battery price to battery capacity in the case of the right part of Figure 14 shows three kinds.With battery price
Go up, battery capacity is all declining in the case of three kinds, but in order to keep the normal operation of integrated energy system and ensure system
Technical security and economy, it is necessary to configure the energy-storage battery of certain capacity as energy reserves.Therefore, this feelings are based on
Condition, the battery capacity in the case of three kinds will gradually tend to be roughly the same.For example, when battery price increases by 40%, three kinds of batteries
Decline is all tended to 83 kilowatt hours by capacity.With the variation of battery price, slight change can occur for the capacity of sodium-sulphur battery.Only
When price declines or increases significant proportion (30%, 40%), significant variation can just occur for the capacity of sodium-sulphur battery.When under price
When dropping 27%, the capacity of lithium battery starts the capacity greater than lead carbon battery.However, the year cost of lithium battery is still higher than lead carbon electricity
Pond, because difference in selling prices can not be made up by increasing capacity bring economic benefit Ji.When price decline 40%, although lithium battery
Configuration capacity is bigger than lead carbon battery, but the year cost of lithium battery is often identical as lead carbon battery because both batteries possess it is good
Difference in selling prices between good economical benefit and battery substantially reduces.
For the actual conditions that the device configuration examined obtained by above planning is held, selection has electricity, cold and thermic load
Typical case's day summer analyzes the IES situation of optimum capacity configuration.Whether the equipment that IES can be verified works according to operation strategy,
And examine the confidence level of Optimized model proposed by the present invention.In verification process, tested with lead carbon battery for BESS.One day
In load hourly it is as shown in figure 15.
The capacity of preferred plan optimization is analyzed, two stages of Bess are all tested on the day of, electricity, heat
With output power such as Figure 16 of equipment in refrigeration subsystem, shown in 17,18.What it is in first figure is solved based on BESS strategy 1
The optimal output of scheme is the power output of the optimal solution based on BESS strategy 2 in the second width figure.
Under two different BESS strategies, the power output of each equipment is as shown in figure 16.The different operation sides BESS
Formula has larger impact to the program results of integrated energy system, therefore running equipment power output is different.Due to wind-driven generator
Group and the quantity of photovoltaic generator are fixed, therefore the power output of two kinds of tactful leeward power generators and photovoltaic generator is
It is identical.The principle of " electricity determining by heat " is abided by PGU operation.Thermic load in the power output and integrated energy system of PGU is close
It is related.Higher in the capacity of the top half of Figure 16, PGU, power output is bigger.Under two kinds of different operational modes of BESS,
It directly charges to BESS during paddy valence of the remaining capacity between 0 to 8 points, and electric power is sold power grid until storage fills
It is full.When the power output of equipment is not able to satisfy loading demand, BESS makes an immediate response and discharges on the top of Figure 16.Electric power is from electricity
Online shopping is bought with meet demand.Until BESS discharges into minimum still, in second digit, system is directly bought from power grid with full
The loading demand of sufficient 4-8 point, and BESS is not responded to.During 8 points to 10 points, BESS electric discharge, and discharged into most in BESS
After lower bound degree, buy power from power grid with meet demand.This is also demonstrated under the operation reserve of the second energy storage strategy, integrates energy
The interaction of source system and power grid is more, and the energy storage response time is shorter.Importantly, in two kinds of operational modes of BESS
Under, main power grid maintains the power-balance of integrated energy system.
Since 24 hours operating statuses of integrated energy system can not reflect the difference between Different Strategies, Figure 17 completely
Show the SOC curve of BESS different operation reserves in one week.It can clearly be seen that BESS mode of operation from the figure
Difference.40 points during par, energy storage charging reaches maximum value under the first strategy.In second of strategy, BESS is carried out
Charging, but when the SOC state of BESS is 0.5, extra electric power is sold to power grid.Under 97 points, the first strategy, energy storage is put
Electricity is to minimum, and under the second strategy, when the SOC value of BESS is equal to 0.5, energy storage stops electric discharge, the shortage of integrated energy system
Electricity is supplemented by power grid purchase.
In addition, heating system of the invention includes generating set (PGU) and electric boiler, cold supply system includes absorption refrigeration
Machine and electric refrigerating machine.Figure 18 illustrates the operating status of heat supply heat and refrigeration system under BESS Different Strategies.It, can in from the graph
To be clear that the variation of different energy storage strategies.
Summer thermic load is smaller, and PGU is as main heating source, with thermic load (FTL) mode under BESS different running method
Work.In the upper left corner of Figure 18, the capacity of PGU is larger, and the output of PGU is larger.And in the lower left corner of Figure 18, due to peak-trough electricity
The capacity of valence, electric boiler increases, and the capacity of PGU is reduced, and the heating load of PGU is also accordingly reduced, but PGU is still used as main heat supply
Source.Electric boiler only works during high heat load and generates less electric power.
Cooling in summer load is big, and Absorption Refrigerator is system mainly for cold source.Absorption Refrigerator is on the day of with volume
Determine Power operation, the refrigeration work consumption of Absorption Refrigerator is fixed value hourly.However, the cooling from Absorption Refrigerator
When power is not still able to satisfy cooling load demand, cooling supply is supplemented using electric refrigerator.Equally, in the lower right corner of Figure 18, by
Increase in the capacity of time-of-use tariffs, electric cooler, the capacity of Absorption Refrigerator is reduced.However, can be with from two figures
Find out, Absorption Refrigerator is low capacity and electric refrigerating machine is large capacity, because electric refrigerating machine is more economical.
From what has been discussed above, the operating status result of BESS difference operation reserve demonstrates electricity, heat, the operation plan of cooling system
Slightly correctly execute.
The invention proposes using electric power as the integrated energy system capacity Optimized model of core, the system is by photovoltaic, wind-force
Generating set, CCHP, electric boiler, electric refrigerating machine and BESS composition.The system is by the annual total cost (ATC) of integrated energy system
As the most concerned optimum target of user.According to from the associated BESS of different operation reserves, genetic algorithm be also used for optimization close
Key design parameter, and obtain and distribute result rationally accordingly.The model is applied to IES project as case study, main to study
Achievement and contribution can be summarized as follows:
(1) in the present invention, the Optimized model of integrated energy system planning considers the different operational modes of BESS, and obtains
Corresponding allocation optimum result is arrived.The first operational mode of BESS is relatively common, and second of operational mode of BESS is examined
Time-of-use tariffs are considered, the annual total cost of integrated energy system reduces 12111 dollars with the configuration of identical lead carbon battery, accounts for
Than about 5.42%.
(2) three kinds of different types of batteries are compared in the optimization process of BESS, discovery lithium battery and lead carbon battery it
Between integrated energy system total annual cost difference it is relatively small.Under two kinds of different operation reserves of BESS, lead carbon battery is comprehensive
The total annual cost for closing energy resource system is slightly lower, and lower than lithium battery 3.16% and 2.77%.In addition, the cost to three kinds of batteries carries out spirit
Basis of sensitivity analysis, when battery cost reduces by 40% or more, the IES with lithium battery is more economical for discovery.
(3) the IES Optimal Planning Model that constructs of the present invention considers power grid, the physical constraint condition of fuel gas network and heat supply network,
And specific program results have been obtained by solving.In addition, network site is optimized.
(4) from the point of view of the optimum results of equipment, electronic water cooler is the cooling hair more more economical than absorption chiller
Electric equipment.After considering peak valley price, the capacity of BESS and PGU are reduced, and the capacity of electric boiler and electric cooler increases, comprehensive energy
Source system and interacting for power grid are more.
In order to facilitate understanding, special diacritic table is for reference:
。
Claims (10)
1. a kind of integrated energy system optimization method for considering different storage energy operation strategies, it is characterised in that: the following steps are included:
S1, a center of energy model is established；
S2, the energy flow equation in center of energy is respectively obtained by center of energy model, energy flow equation includes electric system
Equation, therrmodynamic system equation, cold supply system equation；
S3, the balanced load operation reserve for formulating battery energy storage system respectively and combination peak valley operation reserve；
S4, the cold supply system in integrated energy system, heating system, power supply system are planned respectively；
S5, the Optimized model for establishing integrated energy system, Optimized model include annual total cost and constraint condition；
S6, it is calculated and is solved using Optimized model of the genetic algorithm to integrated energy system.
2. considering the integrated energy system optimization method of different storage energy operation strategies as described in claim 1, it is characterised in that:
The formula of center of energy model in the step S1 are as follows:
Wherein,Respectively represent electric load, thermic load, refrigeration duty；Pb、GbWith HbIt is center of energy respectively
Electric input quantity, gas input quantity, heat input；C is the coefficient of coupling matrix, indicates that electricity, gas, the cold coupling between heat are closed
System.
3. considering the integrated energy system optimization method of different storage energy operation strategies as claimed in claim 2, it is characterised in that:
Power system equations in the step S2 are as follows:
Wherein, Pi pguIt is the output power of PGU, i ∈ b indicates that c-th of PGU unit is located at b-th of center of energy；
The therrmodynamic system equation are as follows:
Wherein,It is the heating output power of PGU；It is the power output situation of assist type electric boiler；It is PGU row
Gas mass flow；It is exhaust temperature of charge；ηauxIt is the transfer efficiency of donkey boiler；It is the combustion of donkey boiler
Gas combustion rate；
The cold supply system equation are as follows:
Wherein,WithIt is the output power of Absorption Refrigerator and electrical chillers；It is the Absorption Refrigerator energy
Consumption；copecIt is the Energy Efficiency Ratio of Absorption Refrigerator.
4. considering the integrated energy system optimization method of different storage energy operation strategies as claimed in claim 3, it is characterised in that:
The balanced load operation reserve of battery energy storage system in the step S3 are as follows: no matter electricity price is in electricity consumption peak value, low ebb or flat
Section, when the output power of integrated energy system is greater than workload demand, energy-storage battery charges, and charges in energy-storage battery
Electricity extra later is sold to power grid；When the power output of integrated energy system is less than loading demand, energy-storage battery electric discharge；When
When energy-storage battery discharges into minimum and system and is still unable to satisfy demand, buy power from power grid to meet workload demand；
The combination peak valley operation reserve of the battery energy storage system are as follows: when the output power of integrated energy system is greater than workload demand
And electricity price, when being in electricity consumption peak value, extra electricity sells to power grid；When electricity price is in level section, energy-storage battery fills
Half electric power；When electricity price is in low power consumption, energy-storage battery is fully charged；When the output power of integrated energy system is little
In workload demand and when electricity price is in electricity consumption peak value, energy-storage battery buys power after discharging into minimum state from power grid；Work as electricity price
When in level section, buy power after energy-storage battery release half electric power from power grid；When electricity price is in low power consumption,
Energy-storage battery does not discharge, and directly buys power from power grid.
5. considering the integrated energy system optimization method of different storage energy operation strategies as claimed in claim 4, it is characterised in that:
The cold supply system in integrated energy system is planned in the step S4, the refrigeration duty of integrated energy system is by absorption system
Cold and electric refrigerating machine work are completed, formula are as follows:
When the refrigeration duty of integrated energy system be less than Absorption Refrigerator power output when, then completely with Absorption Refrigerator into
Row on-the-job satisfaction refrigeration duty, electric refrigerating machine do not work, formula are as follows:
When the output power of Absorption Refrigerator is insufficient for the refrigeration duty of integrated energy system, Absorption Refrigerator is with volume
Determine power work, electric refrigerating machine makes up remaining refrigeration duty, formula are as follows:
6. considering the integrated energy system optimization method of different storage energy operation strategies as claimed in claim 5, it is characterised in that:
The heating system in integrated energy system is planned in the step S4, the thermic load of integrated energy system is by generating set
It is completed with auxiliary electric boiler, formula are as follows:
When the thermic load of integrated energy system is less than the output power of generating set, then worked completely by generating set,
Auxiliary electric boiler does not work, formula are as follows:
When the output power of generating set is insufficient for the thermic load of integrated energy system, generating set will be with rated power
Work, assist type electric boiler make up remaining thermic load, formula are as follows:
7. considering the integrated energy system optimization method of different storage energy operation strategies as claimed in claim 6, it is characterised in that:
The power supply system in integrated energy system is planned in the step S4, the power supply of integrated energy system is by photovoltaic panel, wind
Power turbine, generating set, energy-storage battery, power grid are completed, formula are as follows:
Wherein, Pi pv(t) be photovoltaic panel output electricity, Pi wt(t) be wind turbine output electricity；Pi pguIt is generating set
Output electricity；Pi BESS(t) be energy-storage battery output electricity；Pgrid(t) power grid exports electricity.
8. considering the integrated energy system optimization method of different storage energy operation strategies as claimed in claim 7, it is characterised in that:
The calculation formula of annual total cost in the step S5 are as follows:
F=minATC=Caic+Co&m+Cre+Ce+Cgrid；
Wherein, ATC is integrated energy system year minimum cost；CaicIt is the equipment investment cost of integrated energy system；Co&mIt is comprehensive
The operation expense of energy resource system；CreIt is the replacement cost of integrated energy system；CeThe energy cost of integrated energy system,
CgridIt is the power grid interaction cost of integrated energy system；
Wherein, P and CpIt is the installed capacity and the specific investment cost cost of capacity of equipment；FcrFor capital recovery factor:
Wherein, Co&m,jIt is the unit operation and maintenance cost of equipment j；Wd,t,jIt is the production cost of equipment j；
Wherein, CrjIt is the replacement cost of equipment j, FrpIt is capital recovery factor, a is the service life of equipment, year；Y is comprehensive energy
The service life of source system；
Wherein, cfuelIt is the energy cost of gas consumption；Fd,tThe gas consumption of integrated energy system；cgrid-PIt is electric energy
Purchasing price；Ed,t,pElectricity is bought from power grid；cgrid-sSale of electricity price；Ed,t,sIt is that integrated energy system to power grid sells electricity.
9. considering the integrated energy system optimization method of different storage energy operation strategies as claimed in claim 8, it is characterised in that:
Constraint condition includes power constraint, natural gas grid constraint, heat supply network constraint, energy-storage battery constraint in the step S5；Power constraint
Including power balance constraint and the constraint of electric power stream；Natural gas grid constraint includes flow equilibrium constraint and the day of natural gas node
Right throughput-pressure confines；Heat supply network constraint includes thermodynamic branch characteristic equation, thermodynamics network equilibrium equation and equipment output
Limitation；Energy-storage battery constraint includes the maximum constrained of energy-storage battery capacity and the least commitment of energy-storage battery capacity.
10. considering the integrated energy system optimization method of different storage energy operation strategies as claimed in claim 9, feature exists
In: when in the step S6 to the seismic responses calculated of integrated energy system, variable are as follows:
Wherein, 0,1 selects for the operation reserve of battery energy storage system, NGFor gas ductwork, NEFor power grid.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910289594.5A CN110110897A (en) | 2019-04-11 | 2019-04-11 | A kind of integrated energy system optimization method considering different storage energy operation strategies |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910289594.5A CN110110897A (en) | 2019-04-11 | 2019-04-11 | A kind of integrated energy system optimization method considering different storage energy operation strategies |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN110110897A true CN110110897A (en) | 2019-08-09 |
Family
ID=67485333
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201910289594.5A Pending CN110110897A (en) | 2019-04-11 | 2019-04-11 | A kind of integrated energy system optimization method considering different storage energy operation strategies |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN110110897A (en) |
Cited By (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110428346A (en) * | 2019-09-02 | 2019-11-08 | 云南电网有限责任公司电力科学研究院 | A kind of integrated energy system containing hydrogen energy storage |
| CN111799776A (en) * | 2020-06-09 | 2020-10-20 | 国网山东省电力公司经济技术研究院 | Electric-thermal comprehensive energy system optimization scheduling method considering heat storage characteristics of heat supply network |
| CN111987729A (en) * | 2020-08-10 | 2020-11-24 | 南瑞集团有限公司 | Power storage and gas storage capacity configuration and scheduling optimization method for P2G-containing microgrid |
| CN112465263A (en) * | 2020-12-11 | 2021-03-09 | 国网宁夏电力有限公司经济技术研究院 | Comprehensive energy operation optimization method suitable for multiple scenes |
| CN112636391A (en) * | 2020-12-09 | 2021-04-09 | 国网西藏电力有限公司 | Capacity configuration method for wind-solar-storage distributed energy system |
| CN112668202A (en) * | 2021-01-12 | 2021-04-16 | 山东大学 | Simulation model establishing method and system of comprehensive energy system based on comprehensive energy bus |
| CN112952915A (en) * | 2021-02-24 | 2021-06-11 | 国网山东省电力公司电力科学研究院 | Optimized scheduling method for stabilizing peak valley of power grid in comprehensive energy system |
| CN113793052A (en) * | 2021-09-23 | 2021-12-14 | 长沙理工大学 | Robust optimization scheduling method for regional comprehensive energy system |
Citations (9)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN106950840A (en) * | 2017-05-11 | 2017-07-14 | 山东理工大学 | Towards the integrated energy system layered distribution type control method for coordinating of power network peak clipping |
| CN107025519A (en) * | 2017-03-30 | 2017-08-08 | 厦门大学 | Area distribution formula cool and thermal power polygenerations systeme MIXED INTEGER nonlinear model optimization method |
| CN107665377A (en) * | 2017-09-20 | 2018-02-06 | 国网天津市电力公司 | A kind of multiple source-coupled integrated energy system planing method |
| CN107918919A (en) * | 2017-11-08 | 2018-04-17 | 华北电力大学 | A kind of industrial park integrated energy system Optimized Operation containing control strategy and evaluation system and method |
| CN108764737A (en) * | 2018-05-31 | 2018-11-06 | 南方电网科学研究院有限责任公司 | The assessment of economic benefit method and device of industrial park user side battery energy storage system |
| CN108832656A (en) * | 2018-06-20 | 2018-11-16 | 天津大学 | Turn gas and the micro- energy net multi-objective planning method of renewable energy utilization based on electricity |
| CN109146182A (en) * | 2018-08-24 | 2019-01-04 | 南京理工大学 | The economic load dispatching method of meter and the distributed triple-generation system of a variety of energy storage |
| CN109190824A (en) * | 2018-09-04 | 2019-01-11 | 东北大学 | Game with Coalitions optimizing operation method based on user side regional complex energy resource system |
| CN109409595A (en) * | 2018-10-19 | 2019-03-01 | 南京南瑞继保电气有限公司 | A kind of garden is provided multiple forms of energy to complement each other system dispatching method a few days ago |
-
2019
- 2019-04-11 CN CN201910289594.5A patent/CN110110897A/en active Pending
Patent Citations (9)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN107025519A (en) * | 2017-03-30 | 2017-08-08 | 厦门大学 | Area distribution formula cool and thermal power polygenerations systeme MIXED INTEGER nonlinear model optimization method |
| CN106950840A (en) * | 2017-05-11 | 2017-07-14 | 山东理工大学 | Towards the integrated energy system layered distribution type control method for coordinating of power network peak clipping |
| CN107665377A (en) * | 2017-09-20 | 2018-02-06 | 国网天津市电力公司 | A kind of multiple source-coupled integrated energy system planing method |
| CN107918919A (en) * | 2017-11-08 | 2018-04-17 | 华北电力大学 | A kind of industrial park integrated energy system Optimized Operation containing control strategy and evaluation system and method |
| CN108764737A (en) * | 2018-05-31 | 2018-11-06 | 南方电网科学研究院有限责任公司 | The assessment of economic benefit method and device of industrial park user side battery energy storage system |
| CN108832656A (en) * | 2018-06-20 | 2018-11-16 | 天津大学 | Turn gas and the micro- energy net multi-objective planning method of renewable energy utilization based on electricity |
| CN109146182A (en) * | 2018-08-24 | 2019-01-04 | 南京理工大学 | The economic load dispatching method of meter and the distributed triple-generation system of a variety of energy storage |
| CN109190824A (en) * | 2018-09-04 | 2019-01-11 | 东北大学 | Game with Coalitions optimizing operation method based on user side regional complex energy resource system |
| CN109409595A (en) * | 2018-10-19 | 2019-03-01 | 南京南瑞继保电气有限公司 | A kind of garden is provided multiple forms of energy to complement each other system dispatching method a few days ago |
Non-Patent Citations (2)
| Title |
|---|
| WANG YONGLI.ETC: "Short-term power load forecasting based on IVL-BP neural network technology", 《SYSTEMS ENGINEERING PROCEDIA 4》 * |
| 王永利: "区域客户端能源综合需求侧管理理论及应用研究", 《中国博士学位论文全文数据库》 * |
Cited By (10)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110428346A (en) * | 2019-09-02 | 2019-11-08 | 云南电网有限责任公司电力科学研究院 | A kind of integrated energy system containing hydrogen energy storage |
| CN111799776A (en) * | 2020-06-09 | 2020-10-20 | 国网山东省电力公司经济技术研究院 | Electric-thermal comprehensive energy system optimization scheduling method considering heat storage characteristics of heat supply network |
| CN111987729A (en) * | 2020-08-10 | 2020-11-24 | 南瑞集团有限公司 | Power storage and gas storage capacity configuration and scheduling optimization method for P2G-containing microgrid |
| CN112636391A (en) * | 2020-12-09 | 2021-04-09 | 国网西藏电力有限公司 | Capacity configuration method for wind-solar-storage distributed energy system |
| CN112636391B (en) * | 2020-12-09 | 2023-04-07 | 国网西藏电力有限公司 | Capacity configuration method for wind-solar-storage distributed energy system |
| CN112465263A (en) * | 2020-12-11 | 2021-03-09 | 国网宁夏电力有限公司经济技术研究院 | Comprehensive energy operation optimization method suitable for multiple scenes |
| CN112668202A (en) * | 2021-01-12 | 2021-04-16 | 山东大学 | Simulation model establishing method and system of comprehensive energy system based on comprehensive energy bus |
| CN112952915A (en) * | 2021-02-24 | 2021-06-11 | 国网山东省电力公司电力科学研究院 | Optimized scheduling method for stabilizing peak valley of power grid in comprehensive energy system |
| CN112952915B (en) * | 2021-02-24 | 2023-11-03 | 国网山东省电力公司电力科学研究院 | Optimized scheduling method for stabilizing peak and valley of power grid in comprehensive energy system |
| CN113793052A (en) * | 2021-09-23 | 2021-12-14 | 长沙理工大学 | Robust optimization scheduling method for regional comprehensive energy system |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN110110897A (en) | A kind of integrated energy system optimization method considering different storage energy operation strategies | |
| Guo et al. | Optimal operation of regional integrated energy system considering demand response | |
| CN109165788B (en) | Optimization method of combined cooling heating and power system | |
| Liu et al. | Two-phase collaborative optimization and operation strategy for a new distributed energy system that combines multi-energy storage for a nearly zero energy community | |
| CN104392286B (en) | Consider the micro-capacitance sensor running optimizatin method of supply of cooling, heating and electrical powers and storage energy operation strategy | |
| CN109659927B (en) | Energy storage capacity configuration method of comprehensive energy microgrid considering energy storage participation degree | |
| CN108734350A (en) | A kind of independent method for solving with combined dispatching of the power distribution network containing micro-capacitance sensor | |
| CN109740827A (en) | A kind of regional complex energy system planning optimization method based on dual-layer optimization | |
| CN103997062B (en) | A kind of microgrid energy management control method | |
| CN108009693A (en) | Grid-connected micro-capacitance sensor dual blank-holder based on two-stage demand response | |
| CN111737884B (en) | Multi-target random planning method for micro-energy network containing multiple clean energy sources | |
| Liu et al. | Energy management method of integrated energy system based on collaborative optimization of distributed flexible resources | |
| CN112036934A (en) | Quotation method for participation of load aggregators in demand response considering thermoelectric coordinated operation | |
| Wang et al. | Optimal design of integrated energy system considering different battery operation strategy | |
| CN105955931A (en) | High-density distributed photovoltaic absorption-oriented regional energy network optimizing and scheduling method | |
| CN112600253B (en) | Park comprehensive energy collaborative optimization method and equipment based on optimal energy utilization efficiency | |
| Zhai et al. | Optimization of integrated energy system considering photovoltaic uncertainty and multi-energy network | |
| CN112182887A (en) | Comprehensive energy system planning optimization simulation method | |
| CN109740824A (en) | A kind of multi-energy system Method for optimized planning considering heating network and thermic load | |
| CN114662752A (en) | Comprehensive energy system operation optimization method based on price type demand response model | |
| Lin et al. | Capacity planning of micro energy grid using double-level game model of environment-economic considering dynamic energy pricing strategy | |
| Yang et al. | Coordination and optimization of CCHP microgrid group game based on the interaction of electric and thermal energy considering conditional value at risk | |
| Jia et al. | Optimal operation of multi-energy systems in distributed energy network considering energy storage | |
| CN112883630A (en) | Day-ahead optimized economic dispatching method for multi-microgrid system for wind power consumption | |
| CN110594962B (en) | Distributed energy system optimal configuration method based on fuzzy demand response |
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: 20190809 |