CN107769215A - Garden energy mix system optimization dispatching method based on energy hub - Google Patents
Garden energy mix system optimization dispatching method based on energy hub Download PDFInfo
- Publication number
- CN107769215A CN107769215A CN201711152889.5A CN201711152889A CN107769215A CN 107769215 A CN107769215 A CN 107769215A CN 201711152889 A CN201711152889 A CN 201711152889A CN 107769215 A CN107769215 A CN 107769215A
- Authority
- CN
- China
- Prior art keywords
- mrow
- msub
- energy
- natural gas
- power
- 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
- 238000000034 method Methods 0.000 title claims abstract description 29
- 238000005457 optimization Methods 0.000 title claims abstract description 16
- 230000005611 electricity Effects 0.000 claims abstract description 22
- 238000004146 energy storage Methods 0.000 claims abstract description 12
- VNWKTOKETHGBQD-UHFFFAOYSA-N methane Chemical compound C VNWKTOKETHGBQD-UHFFFAOYSA-N 0.000 claims description 220
- 239000003345 natural gas Substances 0.000 claims description 110
- 239000000203 mixture Substances 0.000 claims description 35
- 239000007789 gas Substances 0.000 claims description 31
- 238000006243 chemical reaction Methods 0.000 claims description 19
- 238000010438 heat treatment Methods 0.000 claims description 15
- 238000004378 air conditioning Methods 0.000 claims description 12
- 230000008878 coupling Effects 0.000 claims description 11
- 238000010168 coupling process Methods 0.000 claims description 11
- 238000005859 coupling reaction Methods 0.000 claims description 11
- 238000009826 distribution Methods 0.000 claims description 11
- 238000003860 storage Methods 0.000 claims description 11
- 239000000567 combustion gas Substances 0.000 claims description 8
- 239000011159 matrix material Substances 0.000 claims description 8
- 230000005540 biological transmission Effects 0.000 claims description 4
- 238000004364 calculation method Methods 0.000 claims description 4
- 238000005057 refrigeration Methods 0.000 claims description 4
- 238000011160 research Methods 0.000 claims description 4
- 238000002485 combustion reaction Methods 0.000 claims description 3
- 238000005265 energy consumption Methods 0.000 claims description 3
- 239000012530 fluid Substances 0.000 claims description 3
- 230000003993 interaction Effects 0.000 claims description 3
- 230000009897 systematic effect Effects 0.000 claims description 3
- 230000009466 transformation Effects 0.000 claims description 3
- 230000006835 compression Effects 0.000 claims 2
- 238000007906 compression Methods 0.000 claims 2
- 238000011144 upstream manufacturing Methods 0.000 claims 1
- 230000009467 reduction Effects 0.000 abstract description 2
- 230000006870 function Effects 0.000 description 11
- 238000010586 diagram Methods 0.000 description 10
- 238000001816 cooling Methods 0.000 description 5
- 238000010248 power generation Methods 0.000 description 3
- 238000004422 calculation algorithm Methods 0.000 description 2
- 229910052799 carbon Inorganic materials 0.000 description 2
- 230000000295 complement effect Effects 0.000 description 2
- 238000011161 development Methods 0.000 description 2
- 230000018109 developmental process Effects 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 238000011156 evaluation Methods 0.000 description 2
- 238000007726 management method Methods 0.000 description 2
- OKTJSMMVPCPJKN-UHFFFAOYSA-N Carbon Chemical compound [C] OKTJSMMVPCPJKN-UHFFFAOYSA-N 0.000 description 1
- 235000006508 Nelumbo nucifera Nutrition 0.000 description 1
- 240000002853 Nelumbo nucifera Species 0.000 description 1
- 235000006510 Nelumbo pentapetala Nutrition 0.000 description 1
- 230000008901 benefit Effects 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 238000011217 control strategy Methods 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 238000005868 electrolysis reaction Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 230000007613 environmental effect Effects 0.000 description 1
- 230000006855 networking Effects 0.000 description 1
- 238000013486 operation strategy Methods 0.000 description 1
- 239000002245 particle Substances 0.000 description 1
- 230000008092 positive effect Effects 0.000 description 1
- 230000001172 regenerating effect Effects 0.000 description 1
- 230000004044 response Effects 0.000 description 1
- 238000004088 simulation Methods 0.000 description 1
- 230000033772 system development Effects 0.000 description 1
- 230000007474 system interaction Effects 0.000 description 1
- 238000012546 transfer Methods 0.000 description 1
- 230000007704 transition Effects 0.000 description 1
Classifications
-
- H—ELECTRICITY
- H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
- H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
- H02J3/00—Circuit arrangements for ac mains or ac distribution networks
- H02J3/04—Circuit arrangements for ac mains or ac distribution networks for connecting networks of the same frequency but supplied from different sources
- H02J3/06—Controlling transfer of power between connected networks; Controlling sharing of load between connected networks
-
- H—ELECTRICITY
- H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
- H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
- H02J3/00—Circuit arrangements for ac mains or ac distribution networks
- H02J3/008—Circuit arrangements for ac mains or ac distribution networks involving trading of energy or energy transmission rights
-
- H—ELECTRICITY
- H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
- H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
- H02J3/00—Circuit arrangements for ac mains or ac distribution networks
- H02J3/28—Arrangements for balancing of the load in a network by storage of energy
-
- H—ELECTRICITY
- H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
- H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
- H02J2203/00—Indexing scheme relating to details of circuit arrangements for AC mains or AC distribution networks
- H02J2203/20—Simulating, e g planning, reliability check, modelling or computer assisted design [CAD]
-
- H—ELECTRICITY
- H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
- H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
- H02J3/00—Circuit arrangements for ac mains or ac distribution networks
- H02J3/003—Load forecast, e.g. methods or systems for forecasting future load demand
Landscapes
- Engineering & Computer Science (AREA)
- Power Engineering (AREA)
- Engine Equipment That Uses Special Cycles (AREA)
Abstract
The present invention relates to a kind of garden energy mix system optimization dispatching method based on energy hub, method has including step:(1) the Power Exchange modeling of energy hub internal;(2) energy-storage system models in energy hub;(3) the Power Exchange modeling between energy hub;(4) garden energy mix system modelling;(5) Optimal Operation Model a few days ago;(6) acquisition of scheduling scheme.The Optimization Scheduling of the present invention can fully excavate the mutual enabling capabilities between electricity/cold heat different energy sources form, it is favorably improved the economy of HEP operations, dispatching method realizes the obvious reduction of garden energy mix system energy supply cost, method application is simple, quick, has good practical value.
Description
Technical field
The invention belongs to integrated energy system technical field, particularly a kind of garden energy mix based on energy hub
System optimization dispatching method.
Background technology
With increase of the modern society to energy demand, energy and environment shortage problem is increasingly serious, to human society
Sustainable development brings great harm.In a distributed manner power supply, cooling heating and power generation system CCHP (combined cooling,
Heating and power) be main supply unit multiple-energy-source garden energy mix system HEP (hybrid energy
Park), due to the characteristic of efficiency of energy utilization, low-carbon environment-friendly is utilized, improved with clean energy resource, hair at full speed is obtained within nearly 2 years
Exhibition.The system can effectively alleviate the contradiction between energy-consuming growth and environmental protection, in China's electric power structural adjustment
Very important effect is played, it is significant to socio-economic development.
HEP is effective from that can carry out the devices such as distributed generation unit, CCHP systems, load, energy storage and control system
It is integrated, meet that user combines the needs of energy supply for electric energy, heat energy and refrigeration, and be incorporated into the power networks by power distribution network, ultimately form
One can be incorporated into the power networks can also isolated operation flexible system.Formulated most according to the operation characteristic of each unit inside HEP
Excellent scheduling scheme, scheduling is optimized to multiple-energy-source in HEP, the complementary abundant consumption with regenerative resource of various energy resources can be achieved
Utilize, reduce system operation cost.Turn however, HEP is one containing various energy resources input, multi output of procucts and various energy resources
The non-same sex complexity body of unit is changed, the mutual conversion for being related to electricity/gas/cold heat link optimizes with complementary, is badly in need of new dispatching party
Method.
Research about garden energy mix system optimization dispatching method has made some progress.Associated specialist is never
The Optimization Scheduling of garden energy mix system is proposed with angle.Compared with conventional electric power type garden energy mix system, HEP
In the presence of by multiple-energy-source link coordinated scheduling come improve reduce energy supply cost active demand, currently in HEP electrically it is cold and hot not
Still lack effective model at present with the transition form between energy form and carry out Unify legislation, cause scheduling scheme to be guarded.
By to following pertinent literature:1st, intelligent grid innovative demonstration area energy internet evaluation index and evaluation method
[J] Power System and its Automation journals, 2016,28 (1):39-45;2、A heuristic operation strategy
for commercial building Parks containing EVs and PV system[J].IEEE
Transactions on Industrial Electronics, 2015,62 (4):2560-2570;3rd, to China's comprehensive energy
Thinking [J] power constructions of System Development, 2015,36 (1):16-25;4th, under the method for operation of electrolysis coupling containing Cogeneration Heat
Microgrid energy complex optimum [J] Power System and its Automation journals, 2016,28 (1):51-57;5th, combined cooling and power system
Multiple target running optimizatin [J] Power System and its Automation journals, 2016,28 (5):62-68;6th, supply of cooling, heating and electrical powers garden mixes
Close energy resource system Optimized Operation universal modeling method [J] Proceedings of the CSEEs, 2013,33 (31):26-33;7、
Hierarchical energy management system for multi-source multi-product Parks
[J] .Renewable Energy, 2015,78:621-630;8th, the microgrid energy management system control strategy of hybrid energy-storing
[J] Power System and its Automation journals, 2016,28 (10):85-91;9th, arranged using the carbon of fuzzy self-tuning particle cluster algorithm
Delegate power transaction cool and thermal power Multiobjective Scheduling [J] Proceedings of the CSEEs, 2014,34 (34):6119-6126;10th, response peak
Cooling heating and power generation system Optimized Operation [J] Power System and its Automation journals of paddy electricity valency, 2016,28 (4):25-30;11、
Optimal scheduling of buildings with energy generation and thermal energy
storage under dynamic electricity pricing using mixed-integer nonlinear
Programming [J] .Applied Energy, 2015,147:49-58;12nd, the cold and hot of Demand-side virtual energy storage system is merged
The energy mix system optimization of CCHP building garden dispatching method [J] Proceedings of the CSEEs, 2017,37 (2):581-
590, this patent content has the difference of essence with pertinent literature, has the creativeness described in Patent Law.
The content of the invention
The purpose of the present invention be in view of the shortcomings of the prior art, and propose it is a kind of based on hotspot stress regulation urban energy it is mutual
Networking tidal current computing method.
The present invention solves its technical problem and takes following technical scheme to realize:
A kind of garden energy mix system optimization dispatching method based on energy hub, the energy collection described in this method
Line device is the control unit of conversion, transmission and a storage comprising diversified forms energy;Garden described in this method mixes
Energy resource system is made up of three energy similarly configured hubs, each energy hub by power network and natural gas network connection,
Power input end mouth P of the power network by power network gas turbine and distributed photovoltaic generation for energy hube1, Pe2, Pe3For
Electricity;Natural gas network is supplied by a natural gas network N and is sent to the natural gas input mouth P of energy hubg1, Pg2, Pg3;
Compressor C is equiped with natural gas contact pipeline between node 1-2 and 1-312And compressor C13, carried for the flowing of natural gas
Voltage supply power, the compressor of natural gas network pipeline is by gas turbine drives, during load equipment that cold and hot source of the gas is provided in power network is
Air-conditioning is entreated, while using heat-storing device come energy storage;HEP is formed with a power network and a natural gas network in this method
HEP is modeler model, it is characterised in that:It is as follows that the method comprising the steps of:
(1) the Power Exchange modeling of energy hub internal
For the energy converter of a single-input single-output, input is with the relation exported:
Lβ=CαβPα (1)
In formula:PαAnd LβThe respectively input and output of systematic steady state;CαβFor the coefficient of coup between inputting and exporting, bag
HEP containing multiple energy converters and various energy resources form, then the coupling of input and output is described by a coupling matrix C
Conjunction relation, i.e.,:
Vectorial P and L are respectively HEP input and output, and the coefficient of coup in coupling matrix C not only turns with conversion equipment
Change that efficiency is relevant, also the distribution coefficient with the energy in different switching device is relevant, introduces an energy distribution coefficient ν here, and 0
≤ ν≤1, e.g., ν P represent directly to supply the electric energy of electric load, (1- ν) PeThen represent the electric energy of supply central air-conditioning;
The conversion efficiency of the electric energy of gas turbine and central air-conditioning is set as constant,WithRespectively natural gas passes through
Cross the conversion efficiency that gas turbine is converted into electric power and heat energy;ηACThe Energy Efficiency Ratio of refrigeration and heating for central air-conditioning, and then
Arrive
Wherein:PeAnd PgIt is energy hub and power network and the energy interaction value of natural gas grid respectively;LeAnd LhRespectively can
The electric load and thermic load that source hub is supplied,
The form for being write as matrix is:
(2) energy-storage system models in energy hub
HEP is using heat-storing device come energy storage, and at the k moment, electricity exchanges power M with hoth,i(k) it is actual and in heat-storing device
The ENERGY E of storageh,i(k) relation is:
In formula,WithStorage energy respectively between heat-storing device and system and the efficiency to release energy;Eh,i(k)
For the energy of heat-storing device physical holding of the stock in moment k;The energy loss of heat-storing device in per a period of time is represented,
Storage device is considered to enter in formula (5), formula (8) can be obtained:
(3) the Power Exchange modeling between energy hub
Trend between energy cluster system, the i.e. trend on the interconnection of connection different energy sources cluster system
Described by steady-state equation, for power network and natural gas network, tide model is established based on node power balance,
1. power network
Direction of energy model is balanced by node complex power to establish, and in node m, node complex power balance can represent such as
Under:
In formula:SmTo inject node m complex power;SmnTo flow to the trend with the node m all nodes being associated, circuit
On trend represented by node voltage amplitude U, vector and line parameter circuit value,
In formula:ymnFor circuit mn transadmittance;ym0For node m self-admittance;
2. natural gas line network
The tide model of the piping network of natural gas grid is also to be established according to node flow balance, following trend side
Journey is applied to the Load flow calculation of all types of isothermal pipelines, and node m volume flow equation of equilibrium is as follows:
In formula:QmTo inject node m natural gas volume flow;QmnIt is the volume flow of pipeline;PmAnd PnRepresent on pipeline
Trip and the pressure in downstream;kmnTo characterize the parameter of pipeline and natural gas fluid;TbFor the temperature of standard state, K;pbFor standard shape
The pressure of state, kPa;DmnFor internal diameter of the pipeline, mm;TfFor the temperature of combustion gas in pipeline, K;G is the relative ratio with air of natural gas
Weight;Z is the compressibility factor of combustion gas;LmmFor gas pipeline length, km;fmnFor the coefficient of friction of gas pipeline, dimensionless,
smnThe direction of air in pipeline flowing is characterized, it is specifically calculated as follows:
For the compressor of natural gas line by gas turbine drives, corresponding power consumption is considered as the extra energy of flow ipe,
The energy consumption of compressor is
Qcom=kcomQmn(Pm-Pk) (15)
In formula:kcomFor compressor pressure ratios;PkFor the pressure of suction port of compressor side;PmFor the pressure of compressor outlet side,
The volume flow Q of natural gas linemnCorresponding direction of energy Pmn, relation between the two is:
In formula:K is heating value of natural gas and the conversion coefficient of electrical power;GHV be natural gas high heating value, MJ/Nm3, Qmn's
Unit is m3/ h, PmnUnit be kW;Because 1MJ=0.278kWh, a transformation ratio 0.278 in formula (16) be present;
(4) garden energy mix system modelling
Garden energy mix system is made up of three energy similarly configured hubs, each energy hub by power network and
Natural gas network connection, power input end mouth P of the power network by power network and distributed photovoltaic generation for energy hube1,
Pe2, Pe3Power supply;Natural gas grid is supplied by a natural gas network N and is sent to the natural gas input mouth P of energy hubg1,
Pg2, Pg3;Compressor C has been installed on natural gas contact pipeline between node 1-2 and 1-312, C13, provided for the flowing of natural gas
Pressure;
(5) Optimal Operation Model a few days ago
1. object function
Because the garden multiple-energy-source hybrid system model that this patent is established is provided by power network, natural gas network and photovoltaic respectively
The energy, so object function uses total purchases strategies and gas consumption cost and minimum object function, such as formula (17)
Shown, electricity price uses tou power price data,
In formula:Cph,iRepresent the Research on electricity price prediction value of i-th hour;PiRepresent i-th hour electrical power bought;Cgas,iRepresent the
The natural gas price predicted value of i hours;PMT,iI-th hour miniature combustion engine electromotive power output, by formula (16) by the body of natural gas network
Product equivalent flow is the trend of electric power networks, therefore the consumption characteristics of natural gas network are equal into electric power networks to handle,
2. constraints
The equality constraint of Optimal Operation Model is the power flow equation (9) of electric power networks, the flow side of natural gas network a few days ago
The equality constraint composition of the equilibrium equation (5) of journey (16) and energy cluster system composition;
Inequality constraints by energy hub input Pi, the flow F of electric power networks and natural gas networka, distribution factor
vi, generator voltage amplitude UmAnd phase angle thetam, generated power output PeiWith idle output Qei, natural gas line pressure pmAnd
The ratio k of compressor delivery pressure and inlet pressurecpLimitation composition;
0≤νi≤1 (20)
(6) acquisition of scheduling scheme
Above-mentioned equation group is solved, obtains scheduling scheme.
The advantages and positive effects of the present invention are:
1st, Optimization Scheduling of the invention can fully excavate the mutual support energy between electricity/cold heat different energy sources form
Power, it is favorably improved the economy of HEP operations;
2nd, dispatching method of the invention realizes the obvious reduction of garden energy mix system energy supply cost, method application letter
It is single, quick, there is good practical value.
Brief description of the drawings
Fig. 1 is energy hub model schematic in the inventive method;
Fig. 2 is the natural gas line model being made up of in the inventive method natural gas compressor (C) and pipeline (P);
Fig. 3 is the garden energy mix system structure diagram that three interconnection energy hubs are formed in the present invention;
Fig. 4 is load model in January schematic diagram in hotel's in the present invention;
Fig. 5 is hotel's August part load model schematic diagram in the present invention;
Fig. 6 is school's load model in January schematic diagram in the present invention;
Fig. 7 is school's August part load model schematic diagram in the present invention;
Fig. 8 is commercial center's load model in January schematic diagram in the present invention;
Fig. 9 is commercial center's August load model schematic diagram in the present invention;
Figure 10 is scheduling scheme schematic diagram when January only has supply of electric power in the present invention;
Scheduling scheme schematic diagram when Figure 11 is energy mix supply in January in the present invention;
Figure 12 is scheduling scheme schematic diagram when August part only has supply of electric power in the present invention;
Figure 13 is the energy mix system schematic of electric power in the present invention, natural gas and photovoltaic.
Embodiment
The embodiment of the present invention is further described below:It is emphasized that embodiment of the present invention is explanation
Property, rather than it is limited, therefore the present invention is not limited to the embodiment described in embodiment, it is every by this area
The other embodiment that technical staff's technique according to the invention scheme is drawn, also belongs to the scope of protection of the invention.
A kind of garden energy mix system optimization dispatching method based on energy hub, its method and step are as follows:
Nominal definition and explanation:
Energy hub:Energy hub is conversion, transmission and a storage that can include diversified forms energy
Control unit, it is the interface platform between different energy facilities and different demands load, and energy hub is from macroscopically
It is connection micro-capacitance sensor and a control centre of bulk power grid or control platform, ultra-short term and real-time online can be passed through
Distributed energy, every state of power distribution network are monitored, optimizes control to each Generation Side and controlled-load, in microcosmic point,
Many basic energy resource facilities similar to factory, big groups of building, rural area, urban area and train etc. can be regarded as energy
Source hub.
The HEP definition and explanation of present patent application:HEP is with a power network and a natural gas network in present patent application
The HEP of composition is modeler model, and garden multiple-energy-source hybrid system model provides energy by power network, natural gas network and photovoltaic respectively
Source, while using heat-storing device come energy storage, the phase in HEP between the coupling of different-energy trend and resulting system
Interaction can be described by the concept of energy hub, as shown in figure 1, a typical HEP system can be abstracted as one
Or multiple energy hub models, comprising input and output, the unit for changing and storing various energy carrier functions forms.
(1) the Power Exchange modeling of energy hub internal
For the energy converter of a single-input single-output, input is with the relation exported:
Lβ=CαβPα (1)
In formula:PαAnd LβThe respectively input and output of systematic steady state;CαβFor the coefficient of coup between inputting and exporting, bag
HEP containing multiple energy converters and various energy resources form, then the coupling of input and output is described by a coupling matrix C
Conjunction relation, i.e.,:
Vectorial P and L are respectively HEP input and output, because a form of energy may enter different energy
Measure in conversion equipment, such as electric power is entered in miniature gas turbine and central air-conditioning simultaneously in HEP in Fig. 1, coupling moment
The coefficient of coup in battle array C is not only relevant with the conversion efficiency of conversion equipment, the distribution system also with the energy in different switching device
Number is relevant, introduces an energy distribution coefficient ν here, and 0≤ν≤1, e.g., ν P represent directly to supply the electric energy of electric load, (1-
ν)PeThen represent the electric energy of supply central air-conditioning;
This patent sets the conversion efficiency of gas turbine and central air-conditioning as constant,WithRespectively natural gas passes through
Cross the conversion efficiency that miniature gas turbine is converted into electric power and heat energy;ηACThe Energy Efficiency Ratio of refrigeration and heating for central air-conditioning, enters
And obtain
Wherein:PeAnd PgIt is energy hub and power network and the energy interaction value of natural gas grid respectively;LeAnd LhRespectively can
The electric load and thermic load that source hub is supplied,
The form for being write as matrix is:
(2) energy-storage system models in energy hub
This patent HEP is using heat-storing device come energy storage, and at the k moment, electricity exchanges power M with hoth,iAnd heat-storing device (k)
The ENERGY E of middle physical holding of the stockh,i(k) relation is:
In formula,WithStorage energy respectively between heat-storing device and system and the efficiency to release energy;Eh,i(k)
For the energy of heat-storing device physical holding of the stock in moment k;The energy loss of heat-storing device in per a period of time is represented,
Storage device is considered to enter in formula (5), formula (8) can be obtained:
(3) the Power Exchange modeling between energy hub
Trend between energy cluster system, the i.e. trend on the interconnection of connection different energy sources cluster system
It can be described by steady-state equation, for power network and natural gas grid, tide model is established based on node power balance,
1. power network
Direction of energy model is balanced by node complex power to establish, and in node m, node complex power balance can represent such as
Under:
In formula:SmTo inject node m complex power;SmnTo flow to the trend with the node m all nodes being associated, circuit
On trend represented by node voltage amplitude U, vector and line parameter circuit value,
In formula:ymnFor circuit mn transadmittance;ym0For node m self-admittance;
2. natural gas line network
The tide model of the piping network of natural gas grid is also to be established according to node flow balance, following trend side
Journey is applied to the Load flow calculation of all types of isothermal pipelines, and node m volume flow equation of equilibrium is as follows:
In formula:QmTo inject node m natural gas volume flow;QmnIt is the volume flow of pipeline;PmAnd PnRepresent on pipeline
Trip and the pressure in downstream;kmnTo characterize the parameter of pipeline and natural gas fluid;TbFor the temperature of standard state, K;pbFor standard shape
The pressure of state, kPa;DmnFor internal diameter of the pipeline, mm;TfFor the temperature of combustion gas in pipeline, K;G is the relative ratio with air of natural gas
Weight;Z is the compressibility factor of combustion gas;LmmFor gas pipeline length, km;fmnFor the coefficient of friction of gas pipeline, dimensionless,
smnThe direction of air in pipeline flowing is characterized, it is specifically calculated as follows:
The compressor of natural gas line needs energy to drive, if driving compressor by miniature gas turbine, accordingly
Power consumption can be considered as the extra energy of flow ipe, as shown in Fig. 2 the energy consumption of compressor is
Qcom=kcomQmn(Pm-Pk) (15)
In formula:kcomFor compressor pressure ratios;PkFor the pressure of suction port of compressor side;PmFor the pressure of compressor outlet side,
The volume flow Q of natural gas linemnCorresponding direction of energy Pmn, relation between the two is:
In formula:K is heating value of natural gas and the conversion coefficient of electrical power;GHV be natural gas high heating value, MJ/Nm3, Qmn's
Unit is m3/ h, PmnUnit be kW;Because 1MJ=0.278kWh, a transformation ratio 0.278 in formula (16) be present;
(4) garden energy mix system modelling
Garden hybrid system containing power network and natural gas grid in this patent is based on several energy line concentrations connected each other
Device system models, accordingly, it is capable to source hub represents energy producers, the interface between consumer and transmission facilities, such as Fig. 3
It is shown, it is the composition of garden energy mix system, system is made up of three energy similarly configured hubs, each energy line concentration
Device is by AC network and natural gas network connection, electricity of the power network by power network and distributed photovoltaic generation for energy hub
Power input port Pe1, Pe2, Pe3Power supply;Natural gas grid is supplied by a natural gas network N and is sent to the natural gas of energy hub
Input port Pg1, Pg2, Pg3;Compressor C has been installed on natural gas contact pipeline between node 1-2 and 1-312, C13, it is natural
The flowing of gas provides pressure;
(5) Optimal Operation Model a few days ago
1. object function
Because the garden multiple-energy-source hybrid system model that this patent is established is provided by power network, natural gas network and photovoltaic respectively
The energy, so object function uses total purchases strategies and gas consumption cost and minimum object function, such as formula (17)
Shown, electricity price uses tou power price data,
In formula:Cph,iRepresent the Research on electricity price prediction value of i-th hour;PiRepresent i-th hour electrical power bought;Cgas,iRepresent the
The natural gas price predicted value of i hours;PMT,iI-th hour miniature combustion engine electromotive power output, by formula (16) by the body of natural gas network
Product equivalent flow is the trend of electric power networks, therefore the consumption characteristics of natural gas network are equal into electric power networks to handle,
2. constraints
The equality constraint of Optimal Operation Model is the power flow equation (9) of electric power networks, the flow side of natural gas network a few days ago
The equality constraint composition of the equilibrium equation (5) of journey (16) and energy cluster system composition;
Inequality constraints by energy hub input Pi, the flow F of electric power networks and natural gas networka, distribution factor
vi, generator voltage amplitude UmAnd phase angle thetam, generated power output PeiWith idle output Qei, natural gas line pressure pmAnd
The ratio k of compressor delivery pressure and inlet pressurecpLimitation composition;
0≤νi≤1 (20)
(6) acquisition of scheduling scheme
Above-mentioned equation group is solved, obtains scheduling scheme.
Instantiation and Contrast on effect
(1) garden multiple-energy-source hybrid system parameter
1. electric power networks parameter
This patent includes three nodes altogether using the electric power networks of model, and node 1 is power network and multiple-energy-source hybrid system
Interface, it is set to balance nodes;Node 2 has accessed distributed photovoltaic power generation, is PV node when photovoltaic is contributed, in photovoltaic not
It is PQ nodes when output;Node 3 is PQ nodes, and the voltage of power network is 10.5kV.Design parameter is shown in Table 1 and table 2.
The power node parameter of table 1
The Electrical Power Line Parameter of table 2
2. natural gas network parameter
Natural gas line parameter is shown in Table 3:
The natural gas line parameter of table 3
This patent selects 8500kcal/Nm3=35.56MJ/Nm3As heating value of natural gas transfer standard, therefore according to formula
(16) there is Pmn=127.91Qmn, this patent by the flow of natural gas network by heating value of natural gas be converted into the direction of energy go into
Row optimization, the reference power of whole system is SB=1MVA, by formula PmnIt is as follows to be converted to perunit value:
It is as shown in table 4 that specific natural gas network calculating parameter can be obtained according to formula (13) and table 3;Garden multiple-energy-source hybrid system
Other specification it is as shown in table 5.
The natural gas network calculating parameter of table 4
The other specification of the garden energy mix system of table 5
3. electric price parameter
This patent electricity price data reference, uses the model of tou power price in optimization calculates, and the price in natural gas January is
0.0425/kWh, the price of August part are 0.0405/kWh, and the specific data of the electricity price of January and August part are shown in Table 6.Natural gas
Data changed significantly with seasonality, relevant with the collection capacity of source of the gas, this patent is integrated energy system, can pass through combustion
Gas generates electricity, and by UTILIZATION OF VESIDUAL HEAT IN, therefore the price of natural gas is contained, optimized together with remaining electricity,
The January of table 6 and the electricity price data of August part
4. load parameter
Multiple-energy-source hybrid network in this patent shares 3 energy cluster systems, the load of energy hub 1 in Fig. 3
For hotel's load, the load of energy hub 2 is school's load, and the load of energy hub 3 is commercial center's load, and three kinds negative
Lotus model is the load model of some typical day.Hotel's load January and August part model are shown in Fig. 4 and Fig. 5 respectively;School
Load January and August part model are shown in Fig. 6 and Fig. 7 respectively;Commercial center's load January and August part model are shown in Fig. 8 and figure respectively
9。
(2) simulation result
Optimization scheduling algorithm a few days ago is entered by taking the garden energy mix system of three interconnection energy hub compositions in Fig. 3 as an example
Row checking.
1. scheduling result in January
1) there was only supply of electric power
Now the energy supply of energy cluster system only has electric power, and without natural gas, scheduling result is as shown in Figure 10, this
When system goal function value be:2411.2, natural gas network is not contributed in system as we can see from the figure.
2) electric power, natural gas and photovoltaic joint supply
System includes electric energy, natural gas and photovoltaic, as shown in figure 3,2 nodes have accessed photo-voltaic power supply, goes out in photo-voltaic power supply
The node is considered as PV node when power, the node is considered as PQ nodes when photovoltaic is not contributed, by optimization obtain as
Figure 11 scheduling scheme, now system goal function value be:2308.4$.
2. August part scheduling result
August part only has that scheduling result during supply of electric power is as shown in figure 12, and system goal function value is:2854.5 $, multipotency
Source joint supply scheduling scheme is as shown in figure 13, and system goal function value is:2467.6$.
(3) scheduling scheme contrasts
Scheduling Cost comparisons under different two kinds of powered modes of month are as shown in table 7,
Scheduling Cost comparisons under the different months of table 7 two kinds of powered modes
It can be seen that only supply of electric power when with multiple-energy-source joint supply when compare, it is higher that the former significantly energizes cost,
It can be seen that comprehensive energy supply cost is substantially reduced for the energy resource system of the more single form of garden energy mix hybrid system.
Claims (1)
1. a kind of garden energy mix system optimization dispatching method based on energy hub, the energy line concentration described in this method
Device is the control unit of conversion, transmission and a storage comprising diversified forms energy;Garden hybrid energy described in this method
Source system is made up of three energy similarly configured hubs, and each energy hub is by power network and natural gas network connection, electricity
Power input end mouth P of the net by power network gas turbine and distributed photovoltaic generation for energy hube1, Pe2, Pe3Power supply;
Natural gas grid is supplied by a natural gas network N and is sent to the natural gas input mouth P of energy hubg1, Pg2, Pg3;Node 1-
Compressor C is equiped with natural gas contact pipeline between 2 and 1-312And compressor C13, pressure is provided for the flowing of natural gas,
For the compressor of natural gas network pipeline by gas turbine drives, the load equipment that cold and hot source of the gas is provided in power network is central air-conditioning,
Simultaneously using heat-storing device come energy storage;HEP is using the HEP that a power network and a natural gas network are formed as modeling in this method
Model, it is characterised in that:It is as follows that the method comprising the steps of:
(1) the Power Exchange modeling of energy hub internal
For the energy converter of a single-input single-output, input is with the relation exported:
Lβ=CαβPα (1)
In formula:PαAnd LβThe respectively input and output of systematic steady state;CαβFor the coefficient of coup between inputting and exporting, comprising more
The HEP of individual energy converter and various energy resources form, then closed by a coupling matrix C to describe the coupling of input and output
System, i.e.,:
Vectorial P and L are respectively HEP input and output, and the coefficient of coup in coupling matrix C is not only imitated with the conversion of conversion equipment
Rate is relevant, and also the distribution coefficient with the energy in different switching device is relevant, introduces the energy distribution coefficient ν, 0≤ν here
≤ 1, e.g., ν P represent directly to supply the electric energy of electric load, (1- ν) PeThen represent the electric energy of supply central air-conditioning;
The conversion efficiency of the electric energy of gas turbine and central air-conditioning is set as constant,WithRespectively natural gas passes through combustion gas
Turbine is converted into the conversion efficiency of electric power and heat energy;ηACThe Energy Efficiency Ratio of refrigeration and heating for central air-conditioning, and then obtain
<mrow>
<msub>
<mi>L</mi>
<mi>e</mi>
</msub>
<mo>=</mo>
<msub>
<mi>vP</mi>
<mi>e</mi>
</msub>
<mo>+</mo>
<msubsup>
<mi>&eta;</mi>
<mrow>
<mi>g</mi>
<mi>e</mi>
</mrow>
<mrow>
<mi>C</mi>
<mi>H</mi>
<mi>P</mi>
</mrow>
</msubsup>
<msub>
<mi>P</mi>
<mi>g</mi>
</msub>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>3</mn>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<msub>
<mi>L</mi>
<mi>h</mi>
</msub>
<mo>=</mo>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>-</mo>
<mi>v</mi>
<mo>)</mo>
</mrow>
<msup>
<mi>&eta;</mi>
<mrow>
<mi>A</mi>
<mi>C</mi>
</mrow>
</msup>
<msub>
<mi>P</mi>
<mi>e</mi>
</msub>
<mo>+</mo>
<msubsup>
<mi>&eta;</mi>
<mrow>
<mi>g</mi>
<mi>h</mi>
</mrow>
<mrow>
<mi>C</mi>
<mi>H</mi>
<mi>P</mi>
</mrow>
</msubsup>
<msub>
<mi>P</mi>
<mi>g</mi>
</msub>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>4</mn>
<mo>)</mo>
</mrow>
</mrow>
Wherein:PeAnd PgIt is energy hub and power network and the energy interaction value of natural gas grid respectively;LeAnd LhRespectively energy collection
The electric load and thermic load that line device is supplied,
The form for being write as matrix is:
(2) energy-storage system models in energy hub
HEP is using heat-storing device come energy storage, and at the k moment, electricity exchanges power M with hoth,i(k) and heat-storing device in physical holding of the stock
ENERGY Ehi(k) relation is:
<mrow>
<mtable>
<mtr>
<mtd>
<mrow>
<msub>
<mi>M</mi>
<mrow>
<mi>h</mi>
<mo>,</mo>
<mi>i</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mfrac>
<msub>
<mover>
<mi>E</mi>
<mo>&CenterDot;</mo>
</mover>
<mrow>
<mi>h</mi>
<mo>,</mo>
<mi>i</mi>
</mrow>
</msub>
<msub>
<mi>e</mi>
<mrow>
<mi>h</mi>
<mo>,</mo>
<mi>i</mi>
</mrow>
</msub>
</mfrac>
<mo>=</mo>
<mfrac>
<mn>1</mn>
<msub>
<mi>e</mi>
<mrow>
<mi>h</mi>
<mo>,</mo>
<mi>i</mi>
</mrow>
</msub>
</mfrac>
<mfrac>
<mrow>
<msub>
<mi>dE</mi>
<mrow>
<mi>h</mi>
<mo>,</mo>
<mi>i</mi>
</mrow>
</msub>
</mrow>
<mrow>
<mi>d</mi>
<mi>t</mi>
</mrow>
</mfrac>
<mo>&ap;</mo>
<mfrac>
<mn>1</mn>
<msub>
<mi>e</mi>
<mrow>
<mi>h</mi>
<mo>,</mo>
<mi>i</mi>
</mrow>
</msub>
</mfrac>
<mfrac>
<mrow>
<msub>
<mi>&Delta;E</mi>
<mrow>
<mi>h</mi>
<mo>,</mo>
<mi>i</mi>
</mrow>
</msub>
</mrow>
<mrow>
<mi>&Delta;</mi>
<mi>t</mi>
</mrow>
</mfrac>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mo>=</mo>
<mfrac>
<mn>1</mn>
<msub>
<mi>e</mi>
<mrow>
<mi>h</mi>
<mo>,</mo>
<mi>i</mi>
</mrow>
</msub>
</mfrac>
<mrow>
<mo>(</mo>
<msub>
<mi>E</mi>
<mrow>
<mi>h</mi>
<mo>,</mo>
<mi>i</mi>
</mrow>
</msub>
<mo>(</mo>
<mi>k</mi>
<mo>)</mo>
<mo>-</mo>
<msub>
<mi>E</mi>
<mrow>
<mi>h</mi>
<mo>,</mo>
<mi>i</mi>
</mrow>
</msub>
<mo>(</mo>
<mrow>
<mi>k</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
<mo>)</mo>
<mo>+</mo>
<msubsup>
<mi>E</mi>
<mrow>
<mi>h</mi>
<mo>,</mo>
<mi>i</mi>
</mrow>
<mrow>
<mi>s</mi>
<mi>t</mi>
<mi>b</mi>
</mrow>
</msubsup>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
</mtable>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>6</mn>
<mo>)</mo>
</mrow>
</mrow>
In formula,WithStorage energy respectively between heat-storing device and system and the efficiency to release energy;Eh,i(k) it is storage
The energy of thermal physical holding of the stock in moment k;The energy loss of heat-storing device in per a period of time is represented,
Storage device is considered to enter in formula (5), formula (8) can be obtained:
(3) the Power Exchange modeling between energy hub
Trend between energy cluster system, the i.e. trend on the interconnection of connection different energy sources cluster system pass through
Steady-state equation is described, and for power network and natural gas network, tide model is established based on node power balance,
1. power network
Direction of energy model is balanced by node complex power to establish, and in node m, node complex power balance can represent as follows:
<mrow>
<msub>
<mi>S</mi>
<mi>m</mi>
</msub>
<mo>-</mo>
<munder>
<mo>&Sigma;</mo>
<mrow>
<mi>n</mi>
<mo>&Element;</mo>
<msub>
<mi>N</mi>
<mi>m</mi>
</msub>
</mrow>
</munder>
<msub>
<mi>S</mi>
<mrow>
<mi>m</mi>
<mi>n</mi>
</mrow>
</msub>
<mo>=</mo>
<mn>0</mn>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>9</mn>
<mo>)</mo>
</mrow>
</mrow>
In formula:SmTo inject node m complex power;SmnTo flow to and the trend of all nodes that node m is associated, on circuit
Trend represents by node voltage amplitude U, vector and line parameter circuit value,
<mrow>
<msub>
<mi>S</mi>
<mrow>
<mi>m</mi>
<mi>n</mi>
</mrow>
</msub>
<mo>=</mo>
<msubsup>
<mi>U</mi>
<mi>m</mi>
<mn>2</mn>
</msubsup>
<mo>&CenterDot;</mo>
<mrow>
<mo>(</mo>
<msub>
<mover>
<mi>y</mi>
<mo>*</mo>
</mover>
<mrow>
<mi>m</mi>
<mn>0</mn>
</mrow>
</msub>
<mo>+</mo>
<msub>
<mover>
<mi>y</mi>
<mo>*</mo>
</mover>
<mrow>
<mi>m</mi>
<mi>n</mi>
</mrow>
</msub>
<mo>)</mo>
</mrow>
<mo>-</mo>
<msub>
<mover>
<mi>U</mi>
<mo>&CenterDot;</mo>
</mover>
<mi>m</mi>
</msub>
<mover>
<msub>
<mi>U</mi>
<mi>n</mi>
</msub>
<mo>*</mo>
</mover>
<msub>
<mover>
<mi>y</mi>
<mo>*</mo>
</mover>
<mrow>
<mi>m</mi>
<mi>n</mi>
</mrow>
</msub>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>10</mn>
<mo>)</mo>
</mrow>
</mrow>
In formula:ymnFor circuit mn transadmittance;ym0For node m self-admittance;
2. natural gas line network
The tide model of the piping network of natural gas grid is also to be established according to node flow balance, and following power flow equation is fitted
For the Load flow calculation of all types of isothermal pipelines, node m volume flow equation of equilibrium is as follows:
<mrow>
<msub>
<mi>Q</mi>
<mi>m</mi>
</msub>
<mo>-</mo>
<munder>
<mo>&Sigma;</mo>
<mrow>
<mi>n</mi>
<mo>&Element;</mo>
<msub>
<mi>N</mi>
<mi>m</mi>
</msub>
</mrow>
</munder>
<msub>
<mi>Q</mi>
<mrow>
<mi>m</mi>
<mi>n</mi>
</mrow>
</msub>
<mo>=</mo>
<mn>0</mn>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>11</mn>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<msub>
<mi>Q</mi>
<mrow>
<mi>m</mi>
<mi>n</mi>
</mrow>
</msub>
<mo>=</mo>
<msub>
<mi>k</mi>
<mrow>
<mi>m</mi>
<mi>n</mi>
</mrow>
</msub>
<msub>
<mi>s</mi>
<mrow>
<mi>m</mi>
<mi>n</mi>
</mrow>
</msub>
<msqrt>
<mrow>
<msub>
<mi>s</mi>
<mrow>
<mi>m</mi>
<mi>n</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<msubsup>
<mi>P</mi>
<mi>m</mi>
<mn>2</mn>
</msubsup>
<mo>-</mo>
<msubsup>
<mi>P</mi>
<mi>n</mi>
<mn>2</mn>
</msubsup>
<mo>)</mo>
</mrow>
</mrow>
</msqrt>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>12</mn>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<msub>
<mi>k</mi>
<mrow>
<mi>m</mi>
<mi>n</mi>
</mrow>
</msub>
<mo>=</mo>
<mn>4.7892</mn>
<mo>&times;</mo>
<msup>
<mn>10</mn>
<mrow>
<mo>-</mo>
<mn>5</mn>
</mrow>
</msup>
<mo>&CenterDot;</mo>
<mrow>
<mo>(</mo>
<mfrac>
<msub>
<mi>T</mi>
<mi>b</mi>
</msub>
<msub>
<mi>p</mi>
<mi>b</mi>
</msub>
</mfrac>
<mo>)</mo>
</mrow>
<mo>&CenterDot;</mo>
<msqrt>
<mfrac>
<msubsup>
<mi>D</mi>
<mrow>
<mi>m</mi>
<mi>n</mi>
</mrow>
<mn>5</mn>
</msubsup>
<mrow>
<msub>
<mi>T</mi>
<mi>f</mi>
</msub>
<msub>
<mi>GZL</mi>
<mrow>
<mi>m</mi>
<mi>m</mi>
</mrow>
</msub>
<msub>
<mi>f</mi>
<mrow>
<mi>m</mi>
<mi>n</mi>
</mrow>
</msub>
</mrow>
</mfrac>
</msqrt>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>13</mn>
<mo>)</mo>
</mrow>
</mrow>
In formula:QmTo inject node m natural gas volume flow;QmnIt is the volume flow of pipeline;PmAnd PnRepresent ducts upstream and
The pressure in downstream;kmnTo characterize the parameter of pipeline and natural gas fluid;TbFor the temperature of standard state, K;pbFor standard state
Pressure, kPa;DmnFor internal diameter of the pipeline, mm;TfFor the temperature of combustion gas in pipeline, K;G is the relative proportion with air of natural gas;Z is
The compressibility factor of combustion gas;LmmFor gas pipeline length, km;fmnFor the coefficient of friction of gas pipeline, dimensionless,
smnThe direction of air in pipeline flowing is characterized, it is specifically calculated as follows:
<mrow>
<msub>
<mi>s</mi>
<mrow>
<mi>m</mi>
<mi>n</mi>
</mrow>
</msub>
<mo>=</mo>
<mfenced open = "{" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<mo>+</mo>
<mn>1</mn>
</mrow>
</mtd>
<mtd>
<mrow>
<mi>i</mi>
<mi>f</mi>
<mi> </mi>
<msub>
<mi>P</mi>
<mi>m</mi>
</msub>
<mo>&GreaterEqual;</mo>
<msub>
<mi>P</mi>
<mi>n</mi>
</msub>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mo>-</mo>
<mn>1</mn>
</mrow>
</mtd>
<mtd>
<mrow>
<mi>e</mi>
<mi>l</mi>
<mi>s</mi>
<mi>e</mi>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>14</mn>
<mo>)</mo>
</mrow>
</mrow>
For the compressor of natural gas line by gas turbine drives, corresponding power consumption is considered as the extra energy of flow ipe, compression
The energy consumption of machine is
Qcom=kcomQmn(Pm-Pk) (15)
In formula:kcomFor compressor pressure ratios;PkFor the pressure of suction port of compressor side;PmFor the pressure of compressor outlet side,
The volume flow Q of natural gas linemnCorresponding direction of energy Pmn, relation between the two is:
<mrow>
<msub>
<mi>P</mi>
<mrow>
<mi>m</mi>
<mi>n</mi>
</mrow>
</msub>
<mo>=</mo>
<mi>k</mi>
<mo>&CenterDot;</mo>
<msub>
<mi>Q</mi>
<mrow>
<mi>m</mi>
<mi>n</mi>
</mrow>
</msub>
<mo>=</mo>
<mfrac>
<mrow>
<mi>G</mi>
<mi>H</mi>
<mi>V</mi>
</mrow>
<mn>0.278</mn>
</mfrac>
<mo>&CenterDot;</mo>
<msub>
<mi>Q</mi>
<mrow>
<mi>m</mi>
<mi>n</mi>
</mrow>
</msub>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>16</mn>
<mo>)</mo>
</mrow>
</mrow>
In formula:K is heating value of natural gas and the conversion coefficient of electrical power;GHV be natural gas high heating value, MJ/Nm3, QmnUnit
For m3/ h, PmnUnit be kW;Because 1MJ=0.278kWh, a transformation ratio 0.278 in formula (16) be present;
(4) garden energy mix system modelling
Garden energy mix system is made up of three energy similarly configured hubs, and each energy hub is by power network and naturally
Gas network connection, power input end mouth P of the power network by power network and distributed photovoltaic generation for energy hube1, Pe2, Pe3
Power supply;Natural gas grid is supplied by a natural gas network N and is sent to the natural gas input mouth P of energy hubg1, Pg2, Pg3;
Compressor C has been installed on natural gas contact pipeline between node 1-2 and 1-312, C13, pressure is provided for the flowing of natural gas;
(5) Optimal Operation Model a few days ago
1. object function
Because the garden multiple-energy-source hybrid system model that this patent is established provides energy by power network, natural gas network and photovoltaic respectively
Source, so object function uses total purchases strategies and gas consumption cost and minimum object function, such as formula (17) institute
To show, electricity price uses tou power price data,
<mrow>
<msub>
<mi>C</mi>
<mn>1</mn>
</msub>
<mo>=</mo>
<mi>m</mi>
<mi>i</mi>
<mi>n</mi>
<munderover>
<mo>&Sigma;</mo>
<mi>i</mi>
<mi>h</mi>
</munderover>
<msub>
<mi>C</mi>
<mrow>
<mi>p</mi>
<mi>h</mi>
<mo>,</mo>
<mi>i</mi>
</mrow>
</msub>
<msub>
<mi>P</mi>
<mi>i</mi>
</msub>
<mo>+</mo>
<msub>
<mi>C</mi>
<mrow>
<mi>g</mi>
<mi>a</mi>
<mi>s</mi>
<mo>,</mo>
<mi>i</mi>
</mrow>
</msub>
<mfrac>
<msub>
<mi>P</mi>
<mrow>
<mi>M</mi>
<mi>T</mi>
<mo>,</mo>
<mi>i</mi>
</mrow>
</msub>
<msubsup>
<mi>&eta;</mi>
<mrow>
<mi>g</mi>
<mi>e</mi>
</mrow>
<mrow>
<mi>C</mi>
<mi>H</mi>
<mi>P</mi>
</mrow>
</msubsup>
</mfrac>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>17</mn>
<mo>)</mo>
</mrow>
</mrow>
In formula:Cph,iRepresent the Research on electricity price prediction value of i-th hour;PiRepresent i-th hour electrical power bought;Cgas,iRepresent i-th hour
Natural gas price predicted value;PMT,iI-th hour miniature combustion engine electromotive power output, by formula (16) by the volume flow of natural gas network
The trend of electric power networks is equivalent to, therefore the consumption characteristics of natural gas network are equal to electric power networks to handle,
2. constraints
The equality constraint of Optimal Operation Model is power flow equation (9), the flow equation of natural gas network of electric power networks a few days ago
(16) and energy cluster system equilibrium equation (5) composition equality constraint composition;
Inequality constraints by energy hub input Pi, the flow F of electric power networks and natural gas networka, distribution factor vi, hair
Electric moter voltage amplitude UmAnd phase angle thetam, generated power output PeiWith idle output Qei, natural gas line pressure pmAnd compression
The ratio k of machine outlet pressure and inlet pressurecpLimitation composition;
<mrow>
<msub>
<munder>
<mi>P</mi>
<mo>&OverBar;</mo>
</munder>
<mi>i</mi>
</msub>
<mo>&le;</mo>
<msub>
<mi>P</mi>
<mi>i</mi>
</msub>
<mo>&le;</mo>
<msub>
<mover>
<mi>P</mi>
<mo>&OverBar;</mo>
</mover>
<mi>i</mi>
</msub>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>18</mn>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<msub>
<munder>
<mi>F</mi>
<mo>&OverBar;</mo>
</munder>
<mi>a</mi>
</msub>
<mo>&le;</mo>
<msub>
<mi>F</mi>
<mi>a</mi>
</msub>
<mo>&le;</mo>
<msub>
<mover>
<mi>F</mi>
<mo>&OverBar;</mo>
</mover>
<mi>a</mi>
</msub>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>19</mn>
<mo>)</mo>
</mrow>
</mrow>
0≤νi≤1 (20)
<mrow>
<msub>
<munder>
<mi>U</mi>
<mo>&OverBar;</mo>
</munder>
<mi>m</mi>
</msub>
<mo>&le;</mo>
<msub>
<mi>U</mi>
<mi>m</mi>
</msub>
<mo>&le;</mo>
<msub>
<mover>
<mi>U</mi>
<mo>&OverBar;</mo>
</mover>
<mi>m</mi>
</msub>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>21</mn>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<msub>
<munder>
<mi>&theta;</mi>
<mo>&OverBar;</mo>
</munder>
<mi>m</mi>
</msub>
<mo>&le;</mo>
<msub>
<mi>&theta;</mi>
<mi>m</mi>
</msub>
<mo>&le;</mo>
<msub>
<mover>
<mi>&theta;</mi>
<mo>&OverBar;</mo>
</mover>
<mi>m</mi>
</msub>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>22</mn>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<msub>
<munder>
<mi>P</mi>
<mo>&OverBar;</mo>
</munder>
<mrow>
<mi>e</mi>
<mi>i</mi>
</mrow>
</msub>
<mo>&le;</mo>
<msub>
<mi>P</mi>
<mrow>
<mi>e</mi>
<mi>i</mi>
</mrow>
</msub>
<mo>&le;</mo>
<msub>
<mover>
<mi>P</mi>
<mo>&OverBar;</mo>
</mover>
<mrow>
<mi>e</mi>
<mi>i</mi>
</mrow>
</msub>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>23</mn>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<msub>
<munder>
<mi>Q</mi>
<mo>&OverBar;</mo>
</munder>
<mrow>
<mi>e</mi>
<mi>i</mi>
</mrow>
</msub>
<mo>&le;</mo>
<msub>
<mi>Q</mi>
<mrow>
<mi>e</mi>
<mi>i</mi>
</mrow>
</msub>
<mo>&le;</mo>
<msub>
<mover>
<mi>Q</mi>
<mo>&OverBar;</mo>
</mover>
<mrow>
<mi>e</mi>
<mi>i</mi>
</mrow>
</msub>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>24</mn>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<msub>
<munder>
<mi>p</mi>
<mo>&OverBar;</mo>
</munder>
<mi>m</mi>
</msub>
<mo>&le;</mo>
<msub>
<mi>p</mi>
<mi>m</mi>
</msub>
<mo>&le;</mo>
<msub>
<mover>
<mi>p</mi>
<mo>&OverBar;</mo>
</mover>
<mi>m</mi>
</msub>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>25</mn>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<msub>
<munder>
<mi>k</mi>
<mo>&OverBar;</mo>
</munder>
<mrow>
<mi>c</mi>
<mi>p</mi>
</mrow>
</msub>
<mo>&le;</mo>
<msub>
<mi>k</mi>
<mrow>
<mi>c</mi>
<mi>p</mi>
</mrow>
</msub>
<mo>=</mo>
<mfrac>
<msub>
<mi>p</mi>
<mi>m</mi>
</msub>
<msub>
<mi>p</mi>
<mi>k</mi>
</msub>
</mfrac>
<mo>&le;</mo>
<msub>
<mover>
<mi>k</mi>
<mo>&OverBar;</mo>
</mover>
<mrow>
<mi>c</mi>
<mi>p</mi>
</mrow>
</msub>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>26</mn>
<mo>)</mo>
</mrow>
<mo>;</mo>
</mrow>
(6) acquisition of scheduling scheme
Above-mentioned equation group is solved, obtains scheduling scheme.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201711152889.5A CN107769215A (en) | 2018-01-19 | 2018-01-19 | Garden energy mix system optimization dispatching method based on energy hub |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201711152889.5A CN107769215A (en) | 2018-01-19 | 2018-01-19 | Garden energy mix system optimization dispatching method based on energy hub |
Publications (1)
Publication Number | Publication Date |
---|---|
CN107769215A true CN107769215A (en) | 2018-03-06 |
Family
ID=61278463
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201711152889.5A Pending CN107769215A (en) | 2018-01-19 | 2018-01-19 | Garden energy mix system optimization dispatching method based on energy hub |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107769215A (en) |
Cited By (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108846506A (en) * | 2018-05-29 | 2018-11-20 | 清华-伯克利深圳学院筹备办公室 | A kind of planing method of energy resource system, device, equipment and storage medium |
CN109783916A (en) * | 2019-01-02 | 2019-05-21 | 大连理工大学 | A kind of air compressor group Optimized Operation decision-making technique based on emulation technology |
CN109919450A (en) * | 2019-02-14 | 2019-06-21 | 国核电力规划设计研究院有限公司 | Solve the game optimization method of comprehensive intelligent energy system scheduling |
CN110163411A (en) * | 2019-04-11 | 2019-08-23 | 华北电力大学 | A kind of regional complex energy system operation optimization method |
CN110263387A (en) * | 2019-05-29 | 2019-09-20 | 国网福建省电力有限公司 | The integrated energy system optimal operation model and its linear processing methods that power distribution network, natural gas grid based on energy hub are established |
CN110263966A (en) * | 2019-05-06 | 2019-09-20 | 天津大学 | Consider the electric-thermal integrated energy system Optimization Scheduling of dynamic heat transfer process |
CN110377973A (en) * | 2019-06-28 | 2019-10-25 | 西安交通大学 | A kind of construction method of normal linearity integrated energy system model |
CN111681130A (en) * | 2020-06-15 | 2020-09-18 | 西安交通大学 | Comprehensive energy system optimization scheduling method considering condition risk value |
CN111768036A (en) * | 2020-06-29 | 2020-10-13 | 国网上海市电力公司 | Power optimization method for interactive operation of comprehensive energy power distribution system and upper-level power grid |
CN117335429A (en) * | 2023-09-28 | 2024-01-02 | 华北电力大学 | Optimal scheduling method and system for power transmission channel and electronic equipment |
-
2018
- 2018-01-19 CN CN201711152889.5A patent/CN107769215A/en active Pending
Cited By (19)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108846506A (en) * | 2018-05-29 | 2018-11-20 | 清华-伯克利深圳学院筹备办公室 | A kind of planing method of energy resource system, device, equipment and storage medium |
CN108846506B (en) * | 2018-05-29 | 2021-02-09 | 清华-伯克利深圳学院筹备办公室 | Energy system planning method, device, equipment and storage medium |
CN109783916B (en) * | 2019-01-02 | 2021-06-18 | 大连理工大学 | Air compressor group optimal scheduling decision method based on simulation technology |
CN109783916A (en) * | 2019-01-02 | 2019-05-21 | 大连理工大学 | A kind of air compressor group Optimized Operation decision-making technique based on emulation technology |
CN109919450A (en) * | 2019-02-14 | 2019-06-21 | 国核电力规划设计研究院有限公司 | Solve the game optimization method of comprehensive intelligent energy system scheduling |
CN109919450B (en) * | 2019-02-14 | 2020-11-24 | 国核电力规划设计研究院有限公司 | Game optimization method for solving scheduling of comprehensive intelligent energy system |
CN110163411A (en) * | 2019-04-11 | 2019-08-23 | 华北电力大学 | A kind of regional complex energy system operation optimization method |
CN110263966A (en) * | 2019-05-06 | 2019-09-20 | 天津大学 | Consider the electric-thermal integrated energy system Optimization Scheduling of dynamic heat transfer process |
CN110263966B (en) * | 2019-05-06 | 2023-04-18 | 天津大学 | Electric-thermal comprehensive energy system optimal scheduling method considering dynamic heat transfer process |
CN110263387A (en) * | 2019-05-29 | 2019-09-20 | 国网福建省电力有限公司 | The integrated energy system optimal operation model and its linear processing methods that power distribution network, natural gas grid based on energy hub are established |
CN110263387B (en) * | 2019-05-29 | 2023-04-07 | 国网福建省电力有限公司 | Energy system operation device based on power distribution network and natural gas network |
CN110377973B (en) * | 2019-06-28 | 2021-08-13 | 西安交通大学 | Construction method of standard linear comprehensive energy system model |
CN110377973A (en) * | 2019-06-28 | 2019-10-25 | 西安交通大学 | A kind of construction method of normal linearity integrated energy system model |
CN111681130A (en) * | 2020-06-15 | 2020-09-18 | 西安交通大学 | Comprehensive energy system optimization scheduling method considering condition risk value |
CN111681130B (en) * | 2020-06-15 | 2024-04-16 | 西安交通大学 | Comprehensive energy system optimal scheduling method considering conditional risk value |
CN111768036A (en) * | 2020-06-29 | 2020-10-13 | 国网上海市电力公司 | Power optimization method for interactive operation of comprehensive energy power distribution system and upper-level power grid |
CN111768036B (en) * | 2020-06-29 | 2023-11-03 | 国网上海市电力公司 | Power optimization method for interactive operation of comprehensive energy distribution system and superior power grid |
CN117335429A (en) * | 2023-09-28 | 2024-01-02 | 华北电力大学 | Optimal scheduling method and system for power transmission channel and electronic equipment |
CN117335429B (en) * | 2023-09-28 | 2024-05-24 | 华北电力大学 | Optimal scheduling method and system for power transmission channel and electronic equipment |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107769215A (en) | Garden energy mix system optimization dispatching method based on energy hub | |
CN108229025B (en) | Economic optimization scheduling method for cooling, heating and power combined supply type multi-microgrid active power distribution system | |
CN107808218A (en) | Urban energy internet tidal current computing method based on hotspot stress regulation | |
CN111738502B (en) | Multi-energy complementary system demand response operation optimization method for promoting surplus wind power consumption | |
CN110866627B (en) | Multi-zone electricity-gas coupling comprehensive energy system optimal scheduling method considering step gas price | |
CN104734168B (en) | Microgrid running optimization system and method based on power and heat combined dispatching | |
CN106447113B (en) | A kind of multizone integrated energy system operation method based on optimal operation model | |
CN109063992A (en) | Consider the power distribution network Expansion Planning method of regional complex energy resource system optimization operation | |
CN110163411A (en) | A kind of regional complex energy system operation optimization method | |
CN112186755B (en) | Flexible load energy storage modeling method for regional comprehensive energy system | |
CN109523076A (en) | A kind of energy storage participate under integrated energy system steady-state optimization analysis method a few days ago | |
Lu et al. | Optimal dispatch of community integrated energy system based on Stackelberg game and integrated demand response under carbon trading mechanism | |
CN111969657B (en) | Multi-energy station system economic optimization operation method based on negotiation game | |
CN105955931A (en) | High-density distributed photovoltaic absorption-oriented regional energy network optimizing and scheduling method | |
CN115170343A (en) | Distributed resource and energy storage collaborative planning method for regional comprehensive energy system | |
CN113141005B (en) | New energy consumption-oriented comprehensive energy system multi-time scale scheduling method | |
CN111461417A (en) | Comprehensive energy system collaborative optimization integrated modeling method based on large system theory | |
CN115906456A (en) | Hydrogen-containing energy IES scheduling optimization model considering response uncertainty of demand side | |
CN113131513B (en) | Method for optimizing operation of electric, thermal and gas conversion system with consideration of carbon emission and storage medium | |
Zhang et al. | Modeling of CHP-EHP coupled energy station considering load side flexibility | |
Li et al. | Integrated energy system model with multi-time scale optimal dispatch method based on a demand response mechanism | |
Sun et al. | Efficiency and benefit evaluation of multi-energy microgrid considering multi-type heat pump | |
Chen et al. | A multi-energy microgrid modelling and optimization method based on exergy theory | |
Zheng et al. | Optimal operation of regional integrated energy system considering flexible heating load | |
Liu et al. | Expansion Planning of Community-Scale Regional Integrated Energy System Considering Grid-Source Coordination: A Cooperative Game Approach |
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 |
Application publication date: 20180306 |
|
RJ01 | Rejection of invention patent application after publication |