CN110070216A - A kind of industrial park integrated energy system economic operation optimization method - Google Patents

A kind of industrial park integrated energy system economic operation optimization method Download PDF

Info

Publication number
CN110070216A
CN110070216A CN201910287614.5A CN201910287614A CN110070216A CN 110070216 A CN110070216 A CN 110070216A CN 201910287614 A CN201910287614 A CN 201910287614A CN 110070216 A CN110070216 A CN 110070216A
Authority
CN
China
Prior art keywords
distribution
energy
garden
transmission
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.)
Granted
Application number
CN201910287614.5A
Other languages
Chinese (zh)
Other versions
CN110070216B (en
Inventor
钟永洁
孙永辉
王建喜
许周
翟苏巍
熊俊杰
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Hohai University HHU
Original Assignee
Hohai University HHU
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Hohai University HHU filed Critical Hohai University HHU
Priority to CN201910287614.5A priority Critical patent/CN110070216B/en
Publication of CN110070216A publication Critical patent/CN110070216A/en
Application granted granted Critical
Publication of CN110070216B publication Critical patent/CN110070216B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06QINFORMATION 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/00Administration; Management
    • G06Q10/04Forecasting or optimisation specially adapted for administrative or management purposes, e.g. linear programming or "cutting stock problem"
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06QINFORMATION 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/00Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
    • G06Q50/06Energy or water supply
    • YGENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02EREDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
    • Y02E40/00Technologies for an efficient electrical power generation, transmission or distribution
    • Y02E40/70Smart grids as climate change mitigation technology in the energy generation sector
    • YGENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02PCLIMATE CHANGE MITIGATION TECHNOLOGIES IN THE PRODUCTION OR PROCESSING OF GOODS
    • Y02P80/00Climate change mitigation technologies for sector-wide applications
    • Y02P80/10Efficient use of energy, e.g. using compressed air or pressurized fluid as energy carrier
    • YGENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y04INFORMATION OR COMMUNICATION TECHNOLOGIES HAVING AN IMPACT ON OTHER TECHNOLOGY AREAS
    • Y04SSYSTEMS INTEGRATING TECHNOLOGIES RELATED TO POWER NETWORK OPERATION, COMMUNICATION OR INFORMATION TECHNOLOGIES FOR IMPROVING THE ELECTRICAL POWER GENERATION, TRANSMISSION, DISTRIBUTION, MANAGEMENT OR USAGE, i.e. SMART GRIDS
    • Y04S10/00Systems supporting electrical power generation, transmission or distribution
    • Y04S10/50Systems or methods supporting the power network operation or management, involving a certain degree of interaction with the load-side end user applications

Landscapes

  • Business, Economics & Management (AREA)
  • Engineering & Computer Science (AREA)
  • Economics (AREA)
  • Human Resources & Organizations (AREA)
  • Strategic Management (AREA)
  • Theoretical Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Health & Medical Sciences (AREA)
  • Marketing (AREA)
  • General Physics & Mathematics (AREA)
  • General Business, Economics & Management (AREA)
  • Tourism & Hospitality (AREA)
  • Quality & Reliability (AREA)
  • Game Theory and Decision Science (AREA)
  • Operations Research (AREA)
  • Development Economics (AREA)
  • Entrepreneurship & Innovation (AREA)
  • Public Health (AREA)
  • Water Supply & Treatment (AREA)
  • General Health & Medical Sciences (AREA)
  • Primary Health Care (AREA)
  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)
  • Supply And Distribution Of Alternating Current (AREA)

Abstract

The invention discloses a kind of industrial park integrated energy system economic operation optimization methods, first on the basis of acquiring garden system data and information, steam power plant, the gas fired-boiler, energy-storage system, photovoltaic, transmission & distribution hot systems, electrical power trans mission/distribution system model of infant industry garden integrated energy system are established, while operation constraint is set;Then infant industry garden integrated energy system running optimizatin method, including gas fired-boiler operation method, energy-storage system operation method, economic operation optimization method are proposed;Finally economic optimization moving model is solved and exports infant industry garden integrated energy system information.The present invention is verified the validity and reasonability of infant industry garden integrated energy system economic operation optimization method by sample calculation analysis, has directive significance for the selection of infant industry garden integrated energy system economic operation optimization method.

Description

A kind of industrial park integrated energy system economic operation optimization method
Technical field
The present invention relates to industrial park integrated energy systems, and in particular to a kind of infant industry garden integrated energy system warp Ji running optimizatin method.
Background technique
Important composition of the integrated energy system as energy resource system of new generation, covers the energy such as power supply, heat supply, gas supply system System is integrated with the energy supply, energy conversion and energy storage device of diversified forms, realizes difference in the different links such as source, net, lotus, storage The coupling of the type energy.Traditional industrial park energy resource system lacks the unified optimization with energy, generally existing energy waste, operation The problems such as with high costs, greatly affected the operational efficiency and economic and environmental benefit of system.In recent years, in a distributed manner the energy, Cogeneration system is the infant industry garden integrated energy system of main supply unit, by excavating each main body in industrial park Respond potentiality, be allowed to power grid carry out good interaction, can preferably meet user with can demand, reduce economical operation cost and Improve comprehensive energy utilization rate.Since there is infant industry garden integrated energy system clean energy resource to utilize, raising using energy source Efficiency, economic and environment-friendly characteristic are developed rapidly for nearly 2 years.Meanwhile infant industry garden integrated energy system can have The contradiction between energy-consuming growth and high-efficiency and economic is alleviated on effect ground, plays in China's electric power structural adjustment very important Effect, great attention by people significant to socio-economic development.Industrial park integrated energy system is as multiple-energy-source System bottom coupling terminal has the characteristics that energy demand is big, consumes energy and concentrate, in the cogeneration for realizing garden, helps garden Enterprise saves with energy cost, reduction equipment operation investment, enhancing garden trade and investment promotion attraction, renewable energy on-site elimination, improves The collaboration of providing multiple forms of energy to complement each other of Demand-side dispatching flexibility, realization is utilized etc. and to be had important practical significance.
Summary of the invention
Goal of the invention: the present invention is directed to a kind of industrial park integrated energy system economic operation optimization methods, fully consider Infant industry garden integrated energy system economy demand mentions for the optimal economic optimization operation of industrial park integrated energy system For method reference and theoretical direction.
Technical solution: a kind of industrial park integrated energy system economic operation optimization method provided by the present invention, including Following steps:
(1) garden system data and information, including thermal load demands data, electrical load requirement data, gas fired-boiler are acquired Capacity, steam power plant's capacity, energy storage system capacity, photovoltaic capacity, transmission & distribution hot systems topology information, transmission & distribution hot systems topology information, Equipment coupled relation information;
(2) moving model and setting operation constraint, including garden device model and its operation constraint and garden transmission & distribution are established System model and its operation constraint, the garden device model includes steam power plant's model, gas fired-boiler model, energy-storage system model With photovoltaic apparatus model, the garden distributing system model includes garden transmission & distribution heat system model and electrical power trans mission/distribution system model;
(3) optimization garden operation reserve, including gas fired-boiler running optimizatin, energy-storage system running optimizatin and economical operation are excellent Change;
(4) garden system information is exported, heat power output, the electric heating power output of steam power plant, energy-storage system including gas fired-boiler Charge-discharge electric power, photovoltaic electricity power output, transmission & distribution hot systems pipeline in hot water temperature, the pipeline flow of transmission & distribution hot systems, transmission & distribution Node voltage, purchase of electricity, air consumption and the economic cost of electric system.
Further, step (2), which is established model and constrained with setting operation, includes:
A, garden device model and operation constrain
A1, steam power plant's model and operation constrain
There are certain functional relation, model expressions for operating cost and the electric heating power output of steam power plant are as follows:
In formula:For steam power plant moment t consumption coal cost;PCHP(t) it contributes for steam power plant in the electricity of moment t; HCHP(t) it contributes for steam power plant in the heat of moment t;aCHP、bCHP、cCHP、dCHP、eCHP、fCHPFor the cost feature parameter of steam power plant.
The operation constraint expression formula of steam power plant are as follows:
In formula: αi、βi、χiIt contributes the i-th inequality constraints parameter of operation area for steam power plant's electric heating, wherein i=1,2,3, 4; Respectively heat the power output lower limit, the upper limit of steam power plant;Respectively steam power plant electricity power output lower limit, The upper limit.
A2, gas fired-boiler model and operation constrain
Gas fired-boiler makees model expression are as follows:
In formula: Hgas(t) heat for gas fired-boiler in moment t is contributed;ηgasFor the effectiveness of performance of gas fired-boiler;For Air consumption of the gas fired-boiler in moment t.
The operation constraint expression formula of gas fired-boiler are as follows:
In formula:Respectively heat the power output lower limit, the upper limit of gas fired-boiler.
A3, energy-storage system model and operation constrain
The model expression of energy-storage system are as follows:
In formula: WbatIt (t) is energy storage capacity of the energy-storage system in moment t;δbatFor the energy-storage system energy loss factor;Δ t is Economic optimization runs simulation step length;PbatIt (t) is energy-storage system in the charge-discharge electric power amount of moment t, positive value is charging, negative value is Electric discharge;For energy-storage system moment t charge volume;For energy-storage system moment t discharge capacity.
The operation constraint expression formula of energy-storage system are as follows:
In formula:Respectively the energy storage capacity lower limit of energy-storage system, the upper limit;Respectively store up The charge power lower limit of energy system, the upper limit;Respectively the discharge power lower limit of energy-storage system, the upper limit;T is Economic optimization runs emulation cycle.
A4, photovoltage model and operation constrain
The model expression of photovoltaic are as follows:
In formula:For photovoltaic moment t maximum predicted electromotive power output;PstFor the photovoltaic under standard test condition Full test generated output;Gpv(t)、GstIt is photovoltaic under the practical intensity of solar radiation of moment t, standard test condition Intensity of solar radiation;εTFor photovoltaic power temperature Dynamic gene;TtemIt (t) is actual work temperature of the photovoltaic in moment t;For Photovoltaic surface temperature under standard test condition.
The operation constraint expression formula of photovoltaic are as follows:
In formula:For photovoltaic moment t reality output electrical power.
B, garden distributing system model and operation constrain
B1, transmission & distribution heat system model and operation constrain
The model expression of transmission & distribution hot systems are as follows:
In formula:The thermal power of heat user demand, the hot merit of head end source of supply respectively in transmission & distribution hot systems Rate; The hot water flow of heat user, head end source of supply is respectively flowed through in transmission & distribution hot systems;It is respectively defeated Hot water temperature with the supply water temperature of inflow heat user, outflow heat user in hot systems;Respectively head end source of supply Supply water temperature, inflow temperature;CpFor the hot water specific heat capacity parameter of transmission & distribution hot systems;Respectively transmission & distribution heat system The head end of pipeline, terminal hot water temperature in system;For environment temperature locating for transmission & distribution hot systems;λpFor the heat transfer coefficient of pipeline; LpFor the length of pipeline;mpFor the hot water flow in pipeline.
The operation constraint expression formula of transmission & distribution hot systems are as follows:
In formula:Pipeline flow lower limit, the upper limit respectively in transmission & distribution hot systems;Respectively Hot water temperature's lower limit, the upper limit are flowed out for the heat user of transmission & distribution hot systems;The respectively head end of transmission & distribution hot systems Source of supply supply water temperature lower limit, the upper limit;mout、ToutRespectively flow out hot water flow, the hot water temperature of transmission & distribution hot systems node; min、TinRespectively flow into hot water flow, the hot water temperature of transmission & distribution hot systems node.
B2, electrical power trans mission/distribution system model and operation constrain
The model expression that electrical power trans mission/distribution system is adopted are as follows:
In formula: Gmn、BmnThe respectively conductance of branch mn, susceptance;ΔPm、ΔQmThe respectively note of electrical power trans mission/distribution system node m Enter active and reactive power;M, n is different node serial numbers in electrical power trans mission/distribution system;θmnFor the phase angle of electrical power trans mission/distribution system node m, n voltage Difference;Um、UnThe respectively voltage magnitude of m, n electrical power trans mission/distribution system node;N ∈ m indicates all and is connected with electrical power trans mission/distribution system node m Branch, and leg endpoint is respectively node m, n.
The operation constraint expression formula of electrical power trans mission/distribution system are as follows:
In formula: SbRespectively the tributary capacity of electrical power trans mission/distribution system, branch hold limit, the tributary capacity upper limit; Respectively the node m of electrical power trans mission/distribution system injects the active lower limit of the power, the upper limit;Respectively power transmission and distribution system The node m of system injects reactive power lower limit, the upper limit;Respectively node m lower voltage limit, the upper limit of electrical power trans mission/distribution system.
Further, step (3) garden running optimizatin method includes:
A, gas fired-boiler operation method
Gas fired-boiler in the integrated energy system of industrial park as just concurrent heating equipment, only when steam power plant electricity power output Just start work when can not individually fully meet the demand of heat user, play Peak Load Adjustment, operation method expression formula are as follows:
In formula: ∑ HloadIt (t) is total thermal power needed for the head end source of supply in moment t transmission & distribution hot systems.
B, energy-storage system operation method
Energy-storage system in the integrated energy system of industrial park is mainly the operating cost for being reduction system, by valley Moment charging, peak value moment electric discharge are to reduce systematic economy operating cost, and on the other hand, at the level values moment, energy-storage system is not opened It is dynamic, charge and discharge number can be effectively reduced, energy storage device service life, operation method expression formula are extended are as follows:
In formula: tp、tt、tvRespectively the peak value rate period of tou power price, level values rate period, valley rate period.
C, economic operation optimization method
On the basis of gas fired-boiler operation method and energy-storage system operation method, it is primarily upon industrial park comprehensive energy system System economy, so that industrial park integrated energy system economic cost minimizes.Industrial park integrated energy system economic cost It mainly include steam power plant's operating cost, gas fired-boiler consumption gas cost, from bulk power grid power purchase expense, economical operation optimizes Method expression formula are as follows:
In formula:For industrial park integrated energy system economic cost total in optimization cycle of operation T;ξgasFor work The price of industry garden integrated energy system purchase natural gas;ξgrid(t) for industrial park integrated energy system in moment t from big electricity The tou power price of online shopping electricity;To buy electricity from bulk power grid in moment t.
The utility model has the advantages that compared with prior art, the mentioned infant industry garden pattern synthesis energy resource system of the present invention has operation Mode flexibly, low-carbon, efficiently etc. advantages;The mentioned gas fired-boiler operation method of the present invention, energy-storage system operation method, economic optimization Operation method more close to, meet, meet actual industrial garden engine request;Equipment established by the present invention and distributing system Model more refines, and the operation constraint of set equipment and distributing system is more comprehensive.
Detailed description of the invention
Fig. 1 is implementation flow chart of the present invention;
Fig. 2 is the example structure chart of infant industry garden integrated energy system;
Fig. 3 is the tou power price figure of infant industry garden integrated energy system power purchase;
Fig. 4 is the electricity power output changed power figure of the electrical power for injecting electrical power trans mission/distribution system node 13 and photovoltaic;
Fig. 5 is that head end source of supply injects thermal power variation diagram in transmission & distribution hot systems;
Fig. 6 is infant industry garden integrated energy system heating power balance figure;
The electric heating power output distribution map of steam power plant during Fig. 7 is economic optimization operation;
General expenses variation diagram during Fig. 8 is economic optimization operation.
Specific embodiment
Technical solution of the present invention is described in detail with specific embodiment with reference to the accompanying drawings of the specification, but it is of the invention Protection scope be not limited to the embodiment.
A kind of industrial park integrated energy system economic operation optimization method, as shown in Figure 1, comprising the following steps:
(1) garden system data and information are acquired
Acquire garden system data and information, including thermal load demands data, electrical load requirement data, gas-fired boiler heat size Amount, energy storage system capacity, photovoltaic capacity, transmission & distribution hot systems topology information, transmission & distribution hot systems topology information, is set steam power plant's capacity Standby coupled relation information etc..
(2) model and setting operation constraint are established
A, garden device model and operation constrain
A1, steam power plant's model and operation constrain
Core cell of the steam power plant as industrial park can provide electricity, thermal energy to garden simultaneously, the operation of steam power plant at There are certain functional relation, model expressions with electric heating power output for this are as follows:
In formula:For steam power plant moment t consumption coal cost;PCHP(t) it contributes for steam power plant in the electricity of moment t; HCHP(t) it contributes for steam power plant in the heat of moment t;aCHP、bCHP、cCHP、dCHP、eCHP、fCHPFor the cost feature parameter of steam power plant.
Physical constraint of the steam power plant when operating normally by steam turbine, generator, electric heating power output may only be at it In operational envelope, constraint expression formula is run are as follows:
In formula: αi、βi、χiIt contributes the i-th inequality constraints parameter of operation area for steam power plant's electric heating, wherein i=1,2,3, 4; Respectively heat the power output lower limit, the upper limit of steam power plant;Respectively steam power plant electricity power output lower limit, The upper limit.
A2, gas fired-boiler model and operation constrain
In the integrated energy system of industrial park, concurrent heating equipment of the gas fired-boiler as steam power plant passes through consumption one Quantitative natural gas carrys out hot water preparing, meets the heat demand of heat user, model expression are as follows:
In formula: Hgas(t) heat for gas fired-boiler in moment t is contributed;ηgasFor the effectiveness of performance of gas fired-boiler;For Air consumption of the gas fired-boiler in moment t.
The heat power output of gas fired-boiler is limited by capacity, runs constraint expression formula are as follows:
In formula:Respectively heat the power output lower limit, the upper limit of gas fired-boiler.
A3, energy-storage system model and operation constrain
Energy-storage system is unit common in the integrated energy system of industrial park, utilizes its time shift feature, Ke Yiqi To the effect of peak load shifting, model expression are as follows:
In formula: WbatIt (t) is energy storage capacity of the energy-storage system in moment t;δbatFor the energy-storage system energy loss factor;Δ t is Economic optimization runs simulation step length;PbatIt (t) is energy-storage system in the charge-discharge electric power amount of moment t, positive value is charging, negative value is Electric discharge;For energy-storage system moment t charge volume;For energy-storage system moment t discharge capacity.
In actual industrial garden integrated energy system, energy-storage system, which needs to meet stringent operation when operating normally, to be wanted Ask, as the economical operation period start be with end energy storage capacity maintain it is constant, any time energy storage capacity should be maintained within a certain range, Charge-discharge electric power is also centainly constrained, and constraint expression formula is run are as follows:
In formula:Respectively the energy storage capacity lower limit of energy-storage system, the upper limit;Respectively store up The charge power lower limit of energy system, the upper limit;Respectively the discharge power lower limit of energy-storage system, the upper limit;T is Economic optimization runs emulation cycle.
A4, photovoltage model and operation constrain
Solar energy is utilized to be sufficiently effective, needs to predict the maximum electricity power output of photovoltaic, the electricity power output of photovoltaic is by the sun Radiation intensity, temperature of photovoltaic work etc. influence, model expression are as follows:
In formula:For photovoltaic moment t maximum predicted electromotive power output;PstFor the photovoltaic under standard test condition Full test generated output;Gpv(t)、GstIt is photovoltaic under the practical intensity of solar radiation of moment t, standard test condition Intensity of solar radiation;εTFor photovoltaic power temperature Dynamic gene;TtemIt (t) is actual work temperature of the photovoltaic in moment t;For Photovoltaic surface temperature under standard test condition.
In the integrated energy system of industrial park, according to the requirement of Different Optimization operation method, in some period, photovoltaic Reality output electrical power is often below maximum predicted electricity power output, runs constraint expression formula are as follows:
In formula:For photovoltaic moment t reality output electrical power.
B, garden distributing system model and operation constrain
B1, transmission & distribution heat system model and operation constrain
In an industrial park integrated energy system, there are a large amount of thermal demand, the heat that transmission & distribution hot systems generate heat source In place of energy transmission & distribution to different heat users, transmission & distribution hot systems can effectively react hot water flow and the variation of hot water temperature, Model expression are as follows:
In formula:The thermal power of heat user demand, the hot merit of head end source of supply respectively in transmission & distribution hot systems Rate; The hot water flow of heat user, head end source of supply is respectively flowed through in transmission & distribution hot systems;Respectively transmission & distribution The supply water temperature of heat user is flowed into hot systems, flows out the hot water temperature of heat user;Respectively head end source of supply Supply water temperature, inflow temperature;CpFor the hot water specific heat capacity parameter of transmission & distribution hot systems;Respectively transmission & distribution hot systems The head end of middle pipeline, terminal hot water temperature;For environment temperature locating for transmission & distribution hot systems;λpFor the heat transfer coefficient of pipeline;Lp For the length of pipeline;mpFor the hot water flow in pipeline.
The operation constraint of transmission & distribution hot systems mainly includes pipeline flow constraint, supply water temperature constrains, thermic load flows out temperature Constraint, node energy conservation constraints run constraint expression formula are as follows:
In formula:Pipeline flow lower limit, the upper limit respectively in transmission & distribution hot systems;Respectively Hot water temperature's lower limit, the upper limit are flowed out for the heat user of transmission & distribution hot systems;The respectively head end of transmission & distribution hot systems Source of supply supply water temperature lower limit, the upper limit;mout、ToutRespectively flow out hot water flow, the hot water temperature of transmission & distribution hot systems node; min、TinRespectively flow into hot water flow, the hot water temperature of transmission & distribution hot systems node.
B2, electrical power trans mission/distribution system model and operation constrain
Electrical power trans mission/distribution system uses AC model, model expression are as follows:
In formula: Gmn、BmnThe respectively conductance of branch mn, susceptance;ΔPm、ΔQmThe respectively note of electrical power trans mission/distribution system node m Enter active and reactive power;M, n is different node serial numbers in electrical power trans mission/distribution system;θmnFor the phase angle of electrical power trans mission/distribution system node m, n voltage Difference;Um、UnThe respectively voltage magnitude of m, n electrical power trans mission/distribution system node;N ∈ m indicates all and is connected with electrical power trans mission/distribution system node m Branch, and leg endpoint is respectively node m, n.
The tributary capacity of electrical power trans mission/distribution system has certain limitations, and node voltage needs guarantee in certain mass quality, power supply Active and reactive injection also in a certain range, run constraint expression formula are as follows:
In formula: SbRespectively the tributary capacity of electrical power trans mission/distribution system, branch hold limit, the tributary capacity upper limit; Respectively the node m of electrical power trans mission/distribution system injects the active lower limit of the power, the upper limit;Respectively power transmission and distribution system The node m of system injects reactive power lower limit, the upper limit;Respectively node m lower voltage limit, the upper limit of electrical power trans mission/distribution system.
(3) garden running optimizatin method
A, gas fired-boiler operation method
Gas fired-boiler in the integrated energy system of industrial park as just concurrent heating equipment, only when steam power plant electricity power output Just start work when can not individually fully meet the demand of heat user, play Peak Load Adjustment, operation method expression formula are as follows:
In formula: ∑ HloadIt (t) is total thermal power needed for the head end source of supply in moment t transmission & distribution hot systems.
B, energy-storage system operation method
Energy-storage system in the integrated energy system of industrial park is mainly the operating cost for being reduction system, by valley Moment charging, peak value moment electric discharge are to reduce systematic economy operating cost, and on the other hand, at the level values moment, energy-storage system is not opened It is dynamic, charge and discharge number can be effectively reduced, energy storage device service life, operation method expression formula are extended are as follows:
In formula: tp、tt、tvRespectively the peak value rate period of tou power price, level values rate period, valley rate period.
C, economic operation optimization method
On the basis of gas fired-boiler operation method and energy-storage system operation method, it is primarily upon industrial park comprehensive energy system System economy, so that industrial park integrated energy system economic cost minimizes.Industrial park integrated energy system economic cost It mainly include steam power plant's operating cost, gas fired-boiler consumption gas cost, from bulk power grid power purchase expense, economical operation optimizes Method expression formula are as follows:
In formula:For industrial park integrated energy system economic cost total in optimization cycle of operation T;ξgasFor work The price of industry garden integrated energy system purchase natural gas;ξgrid(t) for industrial park integrated energy system in moment t from big electricity The tou power price of online shopping electricity;To buy electricity from bulk power grid in moment t.
(4) garden system information is exported
Garden system information is exported, the charge and discharge of the electric heating power output, energy-storage system of heat power output, steam power plant including gas fired-boiler Electrical power, photovoltaic electricity power output, transmission & distribution hot systems pipeline in hot water temperature, the pipeline flow of transmission & distribution hot systems, power transmission and distribution system The information such as node voltage, purchase of electricity, air consumption, the economic cost of system.
(5) sample calculation analysis
A, example introduction
Using typical day in winter as research object, simulation step length is infant industry garden integrated energy system example in the present invention 1 hour, emulation cycle was 24 hours.Infant industry garden integrated energy system structure is as shown in Figure 2 in example: capital equipment packet Include gas fired-boiler, steam power plant, energy-storage system, photovoltaic;Main distributing system include the transmission & distribution hot systems of 11 nodes, 13 nodes it is defeated Distribution system.Gas fired-boiler is heated by consuming the natural gas of air accumulator, and the thermal power of output is injected into transmission & distribution hot systems Node 1 in;Steam power plant is generated electricity by consumption coal and heat supply, and the thermal energy of output is injected into the node 1 of transmission & distribution hot systems In, output electric energy be injected into the node 13 of electrical power trans mission/distribution system;Energy-storage system is connect with bulk power grid, and rechargeable electrical energy carrys out arrogant electricity Net, property are injected into the node 13 of distributing system;Bulk power grid is also connected with the node of electrical power trans mission/distribution system 13;Photovoltaic with The node 11 of electrical power trans mission/distribution system is connected;Transmission & distribution hot systems share 13 nodes, and node 1 is head end source of supply, node 2 to section Point 5 is pipeline connection crosspoint, node 6 to node 11 is heat user;Electrical power trans mission/distribution system shares 13 nodes, its node 1 to section Point 10 is electric load node, node 12 is electric load node, node 11 is photo-voltaic power supply node, node 13 is balance nodes;Storage Gas tank between gas fired-boiler be combustion gas stream;Coal between steam power plant be coal stream;It is input to transmission & distribution hot systems and transmission & distribution heat Hot water in system is thermal energy stream;The electric energy being input in electrical power trans mission/distribution system and electrical power trans mission/distribution system is electric energy streams;When steam power plant When hot undercapacity is to meet thermal load demands, gas fired-boiler just starts operation, is otherwise in shutdown status;Energy-storage system is big Allow to discharge when power grid peak value electricity price, allows to charge in bulk power grid valley electricity price, not work in bulk power grid level values electricity price.
The major parameter setting of infant industry garden integrated energy system example in the present invention: the hot efficiency of gas fired-boiler Rate is 0.9;Energy-storage system rate of energy loss is 0.01;Gas Prices are 300CNY/MWh;The maximum electricity of steam power plant, which is contributed, is 21.1MW/h, maximum heat power output are 24.1MW/h;The capacity of energy-storage system is 4MWh;The voltage of node 13 is 1.1pu;Node 11 Voltage be 1.05pu;Tou power price is as shown in Figure 3.
B, interpretation of result
Model program is write based on LINGO18.0 software platform and calls global solver to above-mentioned established infant industry Garden integrated energy system economical operation Optimized model is solved.The electrical power and photovoltaic injected to electrical power trans mission/distribution system node 13 Electricity power output changed power as shown in figure 4, being injected into the maximum electricity that the electrical power of node 13 has been more than steam power plant in some periods Power output, therefore electrical power trans mission/distribution system will necessarily meet the electric work of electrical power trans mission/distribution system to bulk power grid power purchase or energy-storage system to its electric discharge Rate demand, from Fig. 4 it can be found that photovoltaic electricity power output be mainly distributed on period 7 to the period 20, assume responsibility for a part of defeated Distribution system electric load power, therefore, during this period, the electrical power for being injected into node 13 is lower, presents with the electricity power output of photovoltaic It is negatively correlated;Head end source of supply injection thermal power variation in transmission & distribution hot systems from Fig. 5 as shown in figure 5, can also find, transmission & distribution heat For the thermal power demand of system in the maximum heat power output that a part of period has been more than steam power plant, this just necessarily causes gas fired-boiler to be switched on To supplement thermal energy, gas fired-boiler plays the role of peak regulation;Infant industry garden integrated energy system heating power balance process is such as Shown in Fig. 6, it can be found that steam power plant assumes responsibility for the thermal power demand of main transmission & distribution hot systems, in the period 9 to period from Fig. 6 During 17, gas fired-boiler plays the role of concurrent heating, although the heat power output of steam power plant still has some nargin during this period, due to electricity The limitation for climbing rate of contributing indirectly causes heat power output to receive limitation, meanwhile, it can be found that the electricity of photovoltaic goes out to advocate from Fig. 4 It is distributed between period 7 to the period 20, energy-storage system is also at permission generating state, causes the electric heating power output of steam power plant horizontal Lower, to dissolve solar energy, the electric heating power output of steam power plant is relatively less;The electric heating of steam power plant during economic optimization is run Power output distribution is as shown in fig. 7, according to the result of Fig. 7 it can be found that electric heating power output distribution is relatively concentrated in the upper right corner and the lower right corner The adjusting of part, steam power plant is relatively steady, and electric heating power output is on one point in safe range;The cost knot of economic optimization operation Fruit is as shown in figure 8, as it can be observed in the picture that steam power plant's operating cost accounting highest in infant industry garden integrated energy system, this is anti- The core cell that steam power plant is garden should be gone out, the efficiency for improving steam power plant can further decrease system operation cost.Fig. 8 result Show infant industry garden integrated energy system mainly in period 1 to period 7,24 period of period 19 to period power purchase, during this Predominantly valley electricity price, and power purchase cost is relatively fewer, preferably meets being actually subjected to for infant industry garden integrated energy system It asks.From the above analysis, the model that the present invention establishes and the economic operation optimization method of proposition are rationally effective.

Claims (6)

1. a kind of industrial park integrated energy system economic operation optimization method, it is characterised in that: the following steps are included:
(1) acquire garden system data and information, including thermal load demands data, electrical load requirement data, gas fired-boiler capacity, Steam power plant's capacity, energy storage system capacity, photovoltaic capacity, transmission & distribution hot systems topology information, transmission & distribution hot systems topology information, equipment coupling Close relation information;
(2) moving model and setting operation constraint, including garden device model and its operation constraint and garden distributing system are established Model and its operation constraint, the garden device model includes steam power plant's model, gas fired-boiler model, energy-storage system model and light Device model is lied prostrate, the garden distributing system model includes garden transmission & distribution heat system model and electrical power trans mission/distribution system model;
(3) optimize garden operation reserve, including gas fired-boiler running optimizatin, energy-storage system running optimizatin and economical operation optimization;
(4) garden system information is exported, the charge and discharge of the electric heating power output, energy-storage system of heat power output, steam power plant including gas fired-boiler Electrical power, photovoltaic electricity power output, transmission & distribution hot systems pipeline in hot water temperature, the pipeline flow of transmission & distribution hot systems, power transmission and distribution system Node voltage, purchase of electricity, air consumption and the economic cost of system.
2. a kind of industrial park integrated energy system economic operation optimization method according to claim 1, it is characterised in that: The garden device model and its operation constraint is specific as follows:
Steam power plant's model expression in the garden device model are as follows:
In formula:For steam power plant moment t consumption coal cost;PCHP(t) it contributes for steam power plant in the electricity of moment t;HCHP(t) It contributes for steam power plant in the heat of moment t;aCHP、bCHP、cCHP、dCHP、eCHP、fCHPFor the cost feature parameter of steam power plant;
Steam power plant in the garden equipment operation constraint runs constraint expression formula are as follows:
In formula: αi、βi、χiIt contributes the i-th inequality constraints parameter of operation area for steam power plant's electric heating, wherein i=1,2,3,4; Respectively heat the power output lower limit, the upper limit of steam power plant;Respectively steam power plant electricity power output lower limit, on Limit;
Gas fired-boiler model expression in the garden device model are as follows:
In formula: Hgas(t) heat for gas fired-boiler in moment t is contributed;ηgasFor the effectiveness of performance of gas fired-boiler;For combustion gas Air consumption of the boiler in moment t;
Gas fired-boiler in the garden equipment operation constraint runs constraint expression formula are as follows:
In formula:Respectively heat the power output lower limit, the upper limit of gas fired-boiler;
Energy-storage system model expression in the garden device model are as follows:
In formula: WbatIt (t) is energy storage capacity of the energy-storage system in moment t;δbatFor the energy-storage system energy loss factor;Δ t is economy Optimization operation simulation step length;PbatIt (t) is energy-storage system in the charge-discharge electric power amount of moment t, positive value is charging, negative value is electric discharge;For energy-storage system moment t charge volume;For energy-storage system moment t discharge capacity;
Energy-storage system in the garden equipment operation constraint runs constraint expression formula are as follows:
In formula:Respectively the energy storage capacity lower limit of energy-storage system, the upper limit;Respectively energy storage system The charge power lower limit of system, the upper limit;Respectively the discharge power lower limit of energy-storage system, the upper limit;T is economy Optimization operation emulation cycle;
Photovoltage model expression formula in the garden device model are as follows:
In formula:For photovoltaic moment t maximum predicted electromotive power output;PstMost for the photovoltaic under standard test condition Big test generated output;Gpv(t)、GstFor the sun of the photovoltaic under the practical intensity of solar radiation of moment t, standard test condition Radiation intensity;εTFor photovoltaic power temperature Dynamic gene;TtemIt (t) is actual work temperature of the photovoltaic in moment t;For standard Photovoltaic surface temperature under test condition;
Photovoltaic in the garden equipment operation constraint runs constraint expression formula are as follows:
In formula:For photovoltaic moment t reality output electrical power.
3. a kind of industrial park integrated energy system economic operation optimization method according to claim 1, it is characterised in that: The garden distributing system model and its operation constraint is specific as follows:
Transmission & distribution heat system model expression formula in the garden distributing system model are as follows:
In formula:The thermal power of heat user demand, the thermal power of head end source of supply respectively in transmission & distribution hot systems; The hot water flow of heat user, head end source of supply is respectively flowed through in transmission & distribution hot systems;Respectively transmission & distribution heat The supply water temperature of heat user is flowed into system, flows out the hot water temperature of heat user;The respectively confession of head end source of supply Coolant-temperature gage, inflow temperature;CpFor the hot water specific heat capacity parameter of transmission & distribution hot systems;Respectively in transmission & distribution hot systems The head end of pipeline, terminal hot water temperature;For environment temperature locating for transmission & distribution hot systems;λpFor the heat transfer coefficient of pipeline;LpFor The length of pipeline;mpFor the hot water flow in pipeline;
Transmission & distribution hot systems in the garden distributing system operation constraint run constraint expression formula are as follows:
In formula:Pipeline flow lower limit, the upper limit respectively in transmission & distribution hot systems;It is respectively defeated Heat user with hot systems flows out hot water temperature's lower limit, the upper limit;The respectively head end supply of transmission & distribution hot systems Source supply water temperature lower limit, the upper limit;mout、ToutRespectively flow out hot water flow, the hot water temperature of transmission & distribution hot systems node;min、Tin Respectively flow into hot water flow, the hot water temperature of transmission & distribution hot systems node;
Electrical power trans mission/distribution system model expression in the garden distributing system model are as follows:
In formula: Gmn、BmnThe respectively conductance of branch mn, susceptance;ΔPm、ΔQmRespectively electrical power trans mission/distribution system node m's is injected with Function, reactive power;M, n is different node serial numbers in electrical power trans mission/distribution system;θmnFor the phase angle difference of electrical power trans mission/distribution system node m, n voltage; Um、UnThe respectively voltage magnitude of m, n electrical power trans mission/distribution system node;N ∈ m indicates all branch being connected with electrical power trans mission/distribution system node m Road, and leg endpoint is respectively node m, n;
Electrical power trans mission/distribution system in the garden distributing system operation constraint runs constraint expression formula are as follows:
In formula: SbRespectively the tributary capacity of electrical power trans mission/distribution system, branch hold limit, the tributary capacity upper limit; Respectively the node m of electrical power trans mission/distribution system injects the active lower limit of the power, the upper limit;Respectively electrical power trans mission/distribution system Node m injects reactive power lower limit, the upper limit;Respectively node m lower voltage limit, the upper limit of electrical power trans mission/distribution system.
4. a kind of industrial park integrated energy system economic operation optimization method according to claim 1, it is characterised in that: The corresponding mathematic(al) representation of gas fired-boiler running optimizatin is as follows:
In formula: ∑ HloadIt (t) is total thermal power needed for the head end source of supply in moment t transmission & distribution hot systems.
5. a kind of industrial park integrated energy system economic operation optimization method according to claim 1, it is characterised in that: The corresponding mathematic(al) representation of energy-storage system running optimizatin is as follows:
In formula: tp、tt、tvRespectively the peak value rate period of tou power price, level values rate period, valley rate period.
6. a kind of industrial park integrated energy system economic operation optimization method according to claim 1, it is characterised in that: It is as follows that the economical operation optimizes corresponding function expression:
In formula:For industrial park integrated energy system economic cost total in optimization cycle of operation T;ξgasFor industry park The price of area's integrated energy system purchase natural gas;ξgrid(t) it is purchased in moment t from bulk power grid for industrial park integrated energy system The tou power price of electricity;To buy electricity from bulk power grid in moment t.
CN201910287614.5A 2019-04-11 2019-04-11 Economic operation optimization method for industrial park comprehensive energy system Active CN110070216B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201910287614.5A CN110070216B (en) 2019-04-11 2019-04-11 Economic operation optimization method for industrial park comprehensive energy system

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201910287614.5A CN110070216B (en) 2019-04-11 2019-04-11 Economic operation optimization method for industrial park comprehensive energy system

Publications (2)

Publication Number Publication Date
CN110070216A true CN110070216A (en) 2019-07-30
CN110070216B CN110070216B (en) 2021-02-26

Family

ID=67367352

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201910287614.5A Active CN110070216B (en) 2019-04-11 2019-04-11 Economic operation optimization method for industrial park comprehensive energy system

Country Status (1)

Country Link
CN (1) CN110070216B (en)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111864805A (en) * 2020-08-26 2020-10-30 合肥工业大学 Energy management method of industrial park micro-grid light storage combined power generation device
CN112966873A (en) * 2021-03-18 2021-06-15 南方电网科学研究院有限责任公司 Industrial park comprehensive energy system optimization operation method, device and storage medium

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107609684A (en) * 2017-08-24 2018-01-19 浙江万克新能源科技有限公司 A kind of integrated energy system economic optimization dispatching method based on micro-capacitance sensor
CN108932560A (en) * 2018-06-13 2018-12-04 天津大学 Garden integrated energy system Optimization Scheduling based on Model Predictive Control
CN108960503A (en) * 2018-07-02 2018-12-07 河海大学 The more scene optimization analysis methods of integrated energy system based on interior point method
CN109345045A (en) * 2018-11-29 2019-02-15 东北大学 Electric heating integrated energy system economic load dispatching method based on double multiplier iterative algorithms
CN109523076A (en) * 2018-11-15 2019-03-26 河海大学 A kind of energy storage participate under integrated energy system steady-state optimization analysis method a few days ago

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107609684A (en) * 2017-08-24 2018-01-19 浙江万克新能源科技有限公司 A kind of integrated energy system economic optimization dispatching method based on micro-capacitance sensor
CN108932560A (en) * 2018-06-13 2018-12-04 天津大学 Garden integrated energy system Optimization Scheduling based on Model Predictive Control
CN108960503A (en) * 2018-07-02 2018-12-07 河海大学 The more scene optimization analysis methods of integrated energy system based on interior point method
CN109523076A (en) * 2018-11-15 2019-03-26 河海大学 A kind of energy storage participate under integrated energy system steady-state optimization analysis method a few days ago
CN109345045A (en) * 2018-11-29 2019-02-15 东北大学 Electric heating integrated energy system economic load dispatching method based on double multiplier iterative algorithms

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
HAIBO YU ET AL.: "Economic optimization for configuration and sizing of micro integrated energy systems", 《J. MOD. POWER SYST. CLEAN ENERGY》 *
王伟亮 等: "考虑运行约束的区域电力–天然气–热力综合能源系统能量流优化分析", 《中国电机工程学报》 *

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111864805A (en) * 2020-08-26 2020-10-30 合肥工业大学 Energy management method of industrial park micro-grid light storage combined power generation device
CN112966873A (en) * 2021-03-18 2021-06-15 南方电网科学研究院有限责任公司 Industrial park comprehensive energy system optimization operation method, device and storage medium

Also Published As

Publication number Publication date
CN110070216B (en) 2021-02-26

Similar Documents

Publication Publication Date Title
Jiang et al. Optimal integrated demand response scheduling in regional integrated energy system with concentrating solar power
CN110210747B (en) Flexible scheduling method for electric-heating-gas interconnection energy system
CN109242350A (en) The cooling heating and power generation system capacity configuration optimizing method of meter and translatable load
CN103593711B (en) A kind of distributed power source Optimal Configuration Method
CN109993345B (en) Garden-oriented dynamic economic dispatching method for multi-energy complementary system for island operation
CN110543966A (en) day-ahead scheduling optimization method for micro-energy grid with combined supply of electricity, heat and gas
CN110163415A (en) A kind of multipotency streaming system multi objective fuzzy cooperative optimization method under Study on Variable Condition Features
Yuan et al. A control strategy for distributed energy system considering the state of thermal energy storage
CN110070216A (en) A kind of industrial park integrated energy system economic operation optimization method
CN115330021A (en) Comprehensive energy operation optimization system and method considering methane electric heat utilization ratio
CN115146868A (en) Electricity-heat gas coupled rural micro energy system collaborative planning method and terminal equipment
Li et al. Planning model of integrated energy system considering P2G and energy storage
CN109255487A (en) A kind of integrated energy system optimization method based on normalized matrix model
CN108960564A (en) Supply of cooling, heating and electrical powers type microgrid distribution system dispatching method based on centralization modeling
CN115906456A (en) Hydrogen-containing energy IES scheduling optimization model considering response uncertainty of demand side
Zhu et al. Demand response strategy for promoting distributed power consumption based on resident load participation
Zheng et al. Optimal rural integrated energy system configuration against the background of the rural energy transformation strategy
CN112270481A (en) Multi-target planning method and system for power and natural gas coupling system and storage medium
Wang et al. Source-load-storage coordinated optimal scheduling of virtual power plant participating demand response
Gong et al. Operation mode of integrated energy system considering behavior characteristics of multiple subsystems
Ma et al. Optimization of Integrated Energy System Scheduling Considering Stepped Carbon Emission Trading Mechanism and Electricity Heat Demand Response
Zhong et al. Economic optimization operation strategy for industrial park integrated energy system
Yang et al. Multi-objective Day-ahead Scheduling of Regional Integrated Energy Systems Considering Energy Sharing
Rommyonegge et al. Generation Optimization and Planning To Boost Affordable Electrification: A Case Study of Nias Island System as an Outcome of Power System Operation and Control Course at Institut Teknologi Bandung
Li et al. Integrated demand response model of microgrid considering electro-thermal coupling in market environment

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant