CN110070216B - Economic operation optimization method for industrial park comprehensive energy system - Google Patents

Economic operation optimization method for industrial park comprehensive energy system Download PDF

Info

Publication number
CN110070216B
CN110070216B CN201910287614.5A CN201910287614A CN110070216B CN 110070216 B CN110070216 B CN 110070216B CN 201910287614 A CN201910287614 A CN 201910287614A CN 110070216 B CN110070216 B CN 110070216B
Authority
CN
China
Prior art keywords
transmission
distribution
power
heat
park
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.)
Active
Application number
CN201910287614.5A
Other languages
Chinese (zh)
Other versions
CN110070216A (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

Images

Classifications

    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
    • G06Q10/00—Administration; Management
    • G06Q10/04—Forecasting or optimisation specially adapted for administrative or management purposes, e.g. linear programming or "cutting stock problem"
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
    • G06Q50/00—Systems or methods specially adapted for specific business sectors, e.g. utilities or tourism
    • G06Q50/06—Electricity, gas or water supply
    • Y—GENERAL 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
    • Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02E—REDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
    • Y02E40/00—Technologies for an efficient electrical power generation, transmission or distribution
    • Y02E40/70—Smart grids as climate change mitigation technology in the energy generation sector
    • Y—GENERAL 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
    • Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02P—CLIMATE CHANGE MITIGATION TECHNOLOGIES IN THE PRODUCTION OR PROCESSING OF GOODS
    • Y02P80/00—Climate change mitigation technologies for sector-wide applications
    • Y02P80/10—Efficient use of energy, e.g. using compressed air or pressurized fluid as energy carrier
    • Y—GENERAL 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
    • Y04—INFORMATION OR COMMUNICATION TECHNOLOGIES HAVING AN IMPACT ON OTHER TECHNOLOGY AREAS
    • Y04S—SYSTEMS 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/00—Systems supporting electrical power generation, transmission or distribution
    • Y04S10/50—Systems or methods supporting the power network operation or management, involving a certain degree of interaction with the load-side end user applications

Abstract

The invention discloses an economic operation optimization method of an industrial park comprehensive energy system, which comprises the steps of firstly establishing a thermal power plant, a gas boiler, an energy storage system, a photovoltaic power system, a heat transmission and distribution system and a power transmission and distribution system model of a novel industrial park comprehensive energy system on the basis of collecting park system data and information, and simultaneously setting operation constraints; then, a novel operation optimization method of the industrial park comprehensive energy system is provided, wherein the operation optimization method comprises a gas boiler operation method, an energy storage system operation method and an economic operation optimization method; and finally, solving the economic optimization operation model and outputting the novel industrial park comprehensive energy system information. The invention verifies the effectiveness and the rationality of the economic operation optimization method of the novel industrial park comprehensive energy system through example analysis, and has guiding significance for selecting the economic operation optimization method of the novel industrial park comprehensive energy system.

Description

Economic operation optimization method for industrial park comprehensive energy system
Technical Field
The invention relates to an industrial park comprehensive energy system, in particular to a novel industrial park comprehensive energy system economic operation optimization method.
Background
The comprehensive energy system is an important component of a new generation energy system, covers energy systems of power supply, heat supply, gas supply and the like, integrates energy supply, energy conversion and energy storage equipment in various forms, and realizes coupling of different types of energy in different links of source, network, load, storage and the like. The traditional industrial park energy system lacks the unified optimization of energy utilization, generally has the problems of energy waste, high operation cost and the like, and greatly influences the operation efficiency and the economic and environmental benefits of the system. In recent years, the novel industrial park comprehensive energy system taking the distributed energy and the combined heat and power system as main energy supply units has the beneficial interaction with a power grid by exploring the response potential of each main body of the industrial park, can better meet the energy demand of users, reduce the operation economic cost and improve the comprehensive energy utilization rate. The novel industrial park comprehensive energy system has the characteristics of clean energy utilization, energy utilization efficiency improvement, economy and environmental protection, and is rapidly developed in two years. Meanwhile, the novel industrial park comprehensive energy system can effectively relieve the contradiction between energy consumption increase and high-efficiency economy, plays a very important role in adjusting the electric power structure in China, has great significance for the development of the economic society, and is highly valued by people. The industrial park comprehensive energy system is used as a bottom layer coupling terminal of a multi-energy system, has the characteristics of high energy demand, concentrated energy consumption and the like, and has important practical significance in the aspects of realizing combined heat and power supply of a park, helping park enterprises save energy cost, reducing equipment operation investment, enhancing park recruitment attraction, consuming renewable energy on the spot, improving demand side scheduling flexibility, realizing multi-energy complementary cooperative utilization and the like.
Disclosure of Invention
The purpose of the invention is as follows: the invention aims to provide a method for optimizing the economic operation of an industrial park comprehensive energy system, which fully considers the economic requirements of a novel industrial park comprehensive energy system and provides method reference and theoretical guidance for the optimal economic optimization operation of the industrial park comprehensive energy system.
The technical scheme is as follows: the invention provides an economic operation optimization method of an industrial park comprehensive energy system, which comprises the following steps:
(1) collecting park system data and information, wherein the park system data and the information comprise heat load demand data, electric load demand data, gas boiler capacity, thermal power plant capacity, energy storage system capacity, photovoltaic capacity, heat transmission and distribution system topology information and equipment coupling relation information;
(2) establishing an operation model and setting operation constraints, wherein the operation model comprises a park equipment model and operation constraints thereof, and a park transmission and distribution system model and operation constraints thereof, the park equipment model comprises a thermal power plant model, a gas boiler model, an energy storage system model and a photovoltaic equipment model, and the park transmission and distribution system model comprises a park transmission and distribution thermal system model and a transmission and distribution system model;
(3) optimizing a park operation strategy, which comprises gas boiler operation optimization, energy storage system operation optimization and economic operation optimization;
(4) and outputting park system information which comprises the thermal output of a gas boiler, the electric thermal output of a thermal power plant, the charge and discharge power of an energy storage system, the photovoltaic power output, the hot water temperature in the pipeline of the transmission and distribution thermal system, the pipeline flow of the transmission and distribution thermal system, the node voltage of the transmission and distribution system, the electricity purchasing quantity, the gas consumption quantity and the economic cost.
Further, the step (2) of establishing the model and setting the operation constraint comprises:
A. park equipment model and operational constraints
A1 model and operation constraint of thermal power plant
The operating cost of the thermal power plant and the electric heating output have a certain functional relationship, and the model expression is as follows:
Figure BDA0002023822630000021
in the formula:
Figure BDA0002023822630000022
the coal consumption cost of the thermal power plant at the moment t; pCHP(t) is the electrical output of the thermal power plant at time t; hCHP(t) is the thermal output of the thermal power plant at time t; a isCHP、bCHP、cCHP、dCHP、eCHP、fCHPIs the cost characteristic parameter of the thermal power plant.
The operation constraint expression of the thermal power plant is as follows:
Figure BDA0002023822630000023
in the formula: alpha is alphai、βi、χiFor electric heating of thermal power plantsAn ith inequality constraint parameter of the output operation area, wherein i is 1,2,3 and 4;
Figure BDA0002023822630000024
Figure BDA0002023822630000025
respectively representing the lower limit and the upper limit of the thermal output of the thermal power plant;
Figure BDA0002023822630000026
the lower limit and the upper limit of the electric output of the thermal power plant are respectively.
A2, gas boiler model and operation constraint
The model expression of the gas boiler is as follows:
Figure BDA0002023822630000027
in the formula: hgas(t) is the thermal output of the gas boiler at time t; etagasPerformance efficiency for gas-fired boilers;
Figure BDA0002023822630000028
is the gas consumption of the gas boiler at time t.
The operation constraint expression of the gas boiler is as follows:
Figure BDA0002023822630000031
in the formula:
Figure BDA0002023822630000032
respectively the lower limit and the upper limit of the thermal output of the gas boiler.
A3 model and operation constraint of energy storage system
The model expression of the energy storage system is as follows:
Figure BDA0002023822630000033
in the formula: wbat(t) is the stored energy of the energy storage system at time t; deltabatIs an energy loss factor of the energy storage system; delta t is the economic optimization operation simulation step length; pbat(t) is the charge-discharge power quantity of the energy storage system at the moment t, wherein the positive value is charge, and the negative value is discharge;
Figure BDA0002023822630000034
the charge quantity of the energy storage system at the time t;
Figure BDA0002023822630000035
is the amount of discharge of the energy storage system at time t.
The operation constraint expression of the energy storage system is as follows:
Figure BDA0002023822630000036
in the formula:
Figure BDA0002023822630000037
respectively representing the lower limit and the upper limit of the stored energy of the energy storage system;
Figure BDA0002023822630000038
respectively is the lower limit and the upper limit of the charging power of the energy storage system;
Figure BDA0002023822630000039
respectively representing the lower limit and the upper limit of the discharge power of the energy storage system; and T is an economic optimization operation simulation period.
A4, photovoltaic model and operation constraint
The photovoltaic model expression is:
Figure BDA00020238226300000310
in the formula:
Figure BDA00020238226300000311
is light ofThe maximum predicted output electric power at the moment t; pstThe maximum test power generation power of the photovoltaic under the standard test condition; gpv(t)、GstThe actual solar radiation intensity of the photovoltaic at the moment t and the solar radiation intensity under the standard test condition are shown; epsilonTAdjusting a factor for the photovoltaic power temperature; t istem(t) is the actual operating temperature of the photovoltaic at time t;
Figure BDA0002023822630000041
is the photovoltaic surface temperature under standard test conditions.
The photovoltaic operation constraint expression is as follows:
Figure BDA0002023822630000042
in the formula:
Figure BDA0002023822630000043
the actual output electric power of the photovoltaic at the time t.
B. Park distribution system model and operational constraints
B1, heat transmission and distribution system model and operation constraint
The model expression of the heat transmission and distribution system is as follows:
Figure BDA0002023822630000044
in the formula:
Figure BDA0002023822630000045
respectively the thermal power required by hot users in the heat transmission and distribution system and the thermal power of a head end supply source;
Figure BDA0002023822630000046
Figure BDA0002023822630000047
respectively the hot water flow flowing through the hot user and the head end supply source in the transmission and distribution heat system;
Figure BDA0002023822630000048
the temperature of the water supply flowing into the heat user and the temperature of the hot water flowing out of the heat user in the transmission and distribution heat system are respectively set;
Figure BDA0002023822630000049
the water supply temperature and the water inlet temperature of the head end supply source are respectively; cpThe specific heat capacity parameter of the hot water of the transmission and distribution heat system;
Figure BDA00020238226300000410
respectively the hot water temperature at the head end and the terminal end of a pipeline in the transmission and distribution heat system;
Figure BDA00020238226300000411
the temperature of the environment where the heat transmission and distribution system is located; lambda [ alpha ]pIs the heat transfer coefficient of the pipe; l ispIs the length of the pipe; m ispIs the flow of hot water in the pipeline.
The operation constraint expression of the heat distribution and transmission system is as follows:
Figure BDA00020238226300000412
in the formula:
Figure BDA00020238226300000413
respectively the lower limit and the upper limit of the pipeline flow in the transmission and distribution heat system;
Figure BDA00020238226300000414
the lower limit and the upper limit of the temperature of hot water flowing out of a heat user of the transmission and distribution heat system are respectively set;
Figure BDA00020238226300000415
the lower limit and the upper limit of the water supply temperature of a head end supply source of the transmission and distribution heat system are respectively set; m isout、ToutRespectively the hot water flow and the hot water temperature of the nodes of the outflow heat distribution system; m isin、TinRespectively for the hot water flowing into the nodes of the heat distribution and transmission systemFlow, hot water temperature.
B2, power transmission and distribution system model and operation constraint
The model expression of the power transmission and distribution system is as follows:
Figure BDA0002023822630000051
in the formula: gmn、BmnRespectively the conductance and susceptance of the branch mn; delta Pm、ΔQmRespectively injecting active power and reactive power into the nodes m of the power transmission and distribution system; m and n are numbers of different nodes in the power transmission and distribution system; thetamnThe phase angle difference of the m and n voltages of the nodes of the power transmission and distribution system is obtained; u shapem、UnThe voltage amplitudes of the nodes of the m and n power transmission and distribution systems are respectively; n belongs to m and represents all branches connected with the power transmission and distribution system node m, and the branch end points are respectively the nodes m and n.
The operation constraint expression of the power transmission and distribution system is as follows:
Figure BDA0002023822630000052
in the formula: sb、
Figure BDA0002023822630000053
The branch capacity, the branch capacity lower limit and the branch capacity upper limit of the power transmission and distribution system are respectively;
Figure BDA0002023822630000054
Figure BDA0002023822630000055
respectively injecting a lower limit and an upper limit of active power into a node m of the power transmission and distribution system;
Figure BDA0002023822630000056
respectively injecting a lower limit and an upper limit of reactive power into a node m of the power transmission and distribution system;
Figure BDA0002023822630000057
the lower voltage limit and the upper voltage limit of the node m of the power transmission and distribution system are respectively.
Further, the park operation optimization method in the step (3) comprises the following steps:
A. method for operating gas boiler
The gas boiler is only used as a heat supplementing device in an industrial park comprehensive energy system, and is started to work only when the power output of a thermal power plant cannot meet the requirements of thermal users independently and completely, so that the peak shaving effect is achieved, and the expression of the operation method is as follows:
Figure BDA0002023822630000058
in the formula: sigma Hload(t) is the total thermal power required by the head end supply in the distributed heat system at time t.
B. Energy storage system operation method
Energy storage system among the industrial park comprehensive energy system mainly is for the running cost of reduction system, through charging at the valley moment, the peak moment discharges and reduces system economy running cost, and on the other hand, at the flat moment, energy storage system does not start, can effectual reduction charge-discharge number of times, extension energy storage equipment life, and its operation method expression is:
Figure BDA0002023822630000061
in the formula: t is tp、tt、tvThe time-of-use electricity price is the peak electricity price time interval, the flat electricity price time interval and the valley electricity price time interval of the time-of-use electricity price.
C. Economic operation optimization method
On the basis of a gas boiler operation method and an energy storage system operation method, the economic efficiency of the industrial park comprehensive energy system is mainly concerned, so that the economic cost of the industrial park comprehensive energy system is minimized. The economic cost of the industrial park comprehensive energy system mainly comprises the operation cost of a thermal power plant, the natural gas consumption cost of a gas boiler and the electricity purchasing cost from a large power grid, and the expression of the economic operation optimization method is as follows:
Figure BDA0002023822630000062
in the formula:
Figure BDA0002023822630000063
the total economic cost of the industrial park comprehensive energy system in the optimized operation period T is achieved; xigasA price to purchase natural gas for an industrial park energy system; xigrid(t) the time-of-use electricity price of the industrial park comprehensive energy system for purchasing electricity from the large power grid at the moment t;
Figure BDA0002023822630000064
in order to purchase power from the large grid at time t.
Has the advantages that: compared with the prior art, the novel industrial park type comprehensive energy system provided by the invention has the advantages of flexible operation mode, low carbon, high efficiency and the like; the gas boiler operation method, the energy storage system operation method and the economic optimization operation method provided by the invention are closer, accord with and meet the engineering requirements of an actual industrial park; the model of the equipment and the transmission and distribution system established by the invention is more refined, and the operation constraint of the equipment and the transmission and distribution system is more comprehensive.
Drawings
FIG. 1 is a flow chart of an embodiment of the present invention;
FIG. 2 is a schematic diagram of an exemplary embodiment of the new industrial park energy system;
FIG. 3 is a time-of-use electricity price diagram of the integrated energy system for purchasing electricity in the new industrial park;
fig. 4 is a graph of variation of electrical power injected into the power transmission and distribution system node 13 and photovoltaic electrical output power;
FIG. 5 is a diagram illustrating a variation of head end supply injected thermal power in a heat distribution and transmission system;
FIG. 6 is a heat power balance diagram of a comprehensive energy system of the novel industrial park;
FIG. 7 is a graph of an electrothermal output profile of a thermal power plant during economically optimized operation;
FIG. 8 is a graph illustrating cost changes during an economically optimized operation.
Detailed Description
The technical solution of the present invention is described in detail below with reference to the drawings and the specific embodiments, but the scope of the present invention is not limited to the embodiments.
An economic operation optimization method for an industrial park integrated energy system is shown in fig. 1, and comprises the following steps:
(1) collecting park system data and information
The method comprises the steps of collecting park system data and information, wherein the park system data and the information comprise heat load demand data, electric load demand data, gas boiler capacity, thermal power plant capacity, energy storage system capacity, photovoltaic capacity, heat transmission and distribution system topology information, equipment coupling relation information and the like.
(2) Modeling and setting operational constraints
A. Park equipment model and operational constraints
A1 model and operation constraint of thermal power plant
The thermal power plant can provide power and heat energy to the park as a core unit of the industrial park, the running cost of the thermal power plant and the electric heating output have a certain functional relationship, and the model expression is as follows:
Figure BDA0002023822630000071
in the formula:
Figure BDA0002023822630000072
the coal consumption cost of the thermal power plant at the moment t; pCHP(t) is the electrical output of the thermal power plant at time t; hCHP(t) is the thermal output of the thermal power plant at time t; a isCHP、bCHP、cCHP、dCHP、eCHP、fCHPIs the cost characteristic parameter of the thermal power plant.
The thermal power plant is physically restricted by a steam turbine and a generator when in normal operation, the electric heating output can only be in the safe operation range, and the operation restriction expression is as follows:
Figure BDA0002023822630000073
in the formula: alpha is alphai、βi、χiThe method comprises the following steps of (1) setting an i-th inequality constraint parameter of an electrothermal output operation area of the thermal power plant, wherein i is 1,2,3 and 4;
Figure BDA0002023822630000074
Figure BDA0002023822630000075
respectively representing the lower limit and the upper limit of the thermal output of the thermal power plant;
Figure BDA0002023822630000076
the lower limit and the upper limit of the electric output of the thermal power plant are respectively.
A2, gas boiler model and operation constraint
In the industrial park comprehensive energy system, the gas boiler is used as auxiliary heat supplementing equipment of a thermal power plant, hot water is prepared by consuming a certain amount of natural gas, the heat demand of a heat user is met, and the model expression is as follows:
Figure BDA0002023822630000081
in the formula: hgas(t) is the thermal output of the gas boiler at time t; etagasPerformance efficiency for gas-fired boilers;
Figure BDA0002023822630000082
is the gas consumption of the gas boiler at time t.
The thermal output of a gas boiler is limited by the capacity, and the operation constraint expression is as follows:
Figure BDA0002023822630000083
in the formula:
Figure BDA0002023822630000084
respectively the lower limit and the upper limit of the thermal output of the gas boiler.
A3 model and operation constraint of energy storage system
The energy storage system is a common equipment unit in an industrial park comprehensive energy system, and can play a role in peak clipping and valley filling by utilizing the time-shift characteristic of the energy storage system, and the model expression is as follows:
Figure BDA0002023822630000085
in the formula: wbat(t) is the stored energy of the energy storage system at time t; deltabatIs an energy loss factor of the energy storage system; delta t is the economic optimization operation simulation step length; pbat(t) is the charge-discharge power quantity of the energy storage system at the moment t, wherein the positive value is charge, and the negative value is discharge;
Figure BDA0002023822630000086
the charge quantity of the energy storage system at the time t;
Figure BDA0002023822630000087
is the amount of discharge of the energy storage system at time t.
In an integrated energy system of an actual industrial park, an energy storage system needs to meet strict operation requirements during normal operation, for example, when an economic operation cycle starts and ends, the energy storage is kept constant, the energy storage is kept in a certain range at any time, the charge and discharge power is also restrained to a certain extent, and the operation restraint expression is as follows:
Figure BDA0002023822630000088
in the formula:
Figure BDA0002023822630000089
respectively under the energy storage of the energy storage systemUpper and lower limits;
Figure BDA00020238226300000810
respectively is the lower limit and the upper limit of the charging power of the energy storage system;
Figure BDA00020238226300000811
respectively representing the lower limit and the upper limit of the discharge power of the energy storage system; and T is an economic optimization operation simulation period.
A4, photovoltaic model and operation constraint
In order to fully and effectively utilize solar energy, the maximum photovoltaic power output needs to be predicted, the photovoltaic power output is influenced by solar radiation intensity, photovoltaic working temperature and the like, and the model expression is as follows:
Figure BDA0002023822630000091
in the formula:
Figure BDA0002023822630000092
outputting electric power for maximum prediction of photovoltaic at the moment t; pstThe maximum test power generation power of the photovoltaic under the standard test condition; gpv(t)、GstThe actual solar radiation intensity of the photovoltaic at the moment t and the solar radiation intensity under the standard test condition are shown; epsilonTAdjusting a factor for the photovoltaic power temperature; t istem(t) is the actual operating temperature of the photovoltaic at time t;
Figure BDA0002023822630000093
is the photovoltaic surface temperature under standard test conditions.
In an industrial park comprehensive energy system, according to the requirements of different optimized operation methods, in a certain period of time, the actual output electric power of photovoltaic is often smaller than the maximum predicted electric output, and the operation constraint expression is as follows:
Figure BDA0002023822630000094
in the formula:
Figure BDA0002023822630000095
the actual output electric power of the photovoltaic at the time t.
B. Park distribution system model and operational constraints
B1, heat transmission and distribution system model and operation constraint
In an industrial park comprehensive energy system, a large amount of heat energy needs exist, a transmission and distribution heat system transmits and distributes heat energy generated by a heat source to different heat users, the transmission and distribution heat system can effectively reflect the change of hot water flow and hot water temperature, and the model expression is as follows:
Figure BDA0002023822630000096
in the formula:
Figure BDA0002023822630000097
respectively the thermal power required by hot users in the heat transmission and distribution system and the thermal power of a head end supply source;
Figure BDA0002023822630000098
Figure BDA0002023822630000099
respectively the hot water flow flowing through the hot user and the head end supply source in the transmission and distribution heat system;
Figure BDA00020238226300000910
the temperature of the water supply flowing into the heat user and the temperature of the hot water flowing out of the heat user in the transmission and distribution heat system are respectively set;
Figure BDA00020238226300000911
the water supply temperature and the water inlet temperature of the head end supply source are respectively; cpThe specific heat capacity parameter of the hot water of the transmission and distribution heat system;
Figure BDA0002023822630000101
are respectively asThe hot water temperature at the head end and the terminal end of a pipeline in the heat transmission and distribution system;
Figure BDA0002023822630000102
the temperature of the environment where the heat transmission and distribution system is located; lambda [ alpha ]pIs the heat transfer coefficient of the pipe; l ispIs the length of the pipe; m ispIs the flow of hot water in the pipeline.
The operation constraints of the heat transmission and distribution system mainly comprise pipeline flow constraints, water supply temperature constraints, heat load outflow temperature constraints and node energy conservation constraints, and the operation constraint expression is as follows:
Figure BDA0002023822630000103
in the formula:
Figure BDA0002023822630000104
respectively the lower limit and the upper limit of the pipeline flow in the transmission and distribution heat system;
Figure BDA0002023822630000105
the lower limit and the upper limit of the temperature of hot water flowing out of a heat user of the transmission and distribution heat system are respectively set;
Figure BDA0002023822630000106
the lower limit and the upper limit of the water supply temperature of a head end supply source of the transmission and distribution heat system are respectively set; m isout、ToutRespectively the hot water flow and the hot water temperature of the nodes of the outflow heat distribution system; m isin、TinRespectively the hot water flow and the hot water temperature flowing into the nodes of the heat transmission and distribution system.
B2, power transmission and distribution system model and operation constraint
The power transmission and distribution system adopts an alternating current model, and the model expression is as follows:
Figure BDA0002023822630000107
in the formula: gmn、BmnConductance of branches mn respectivelyAnd susceptance; delta Pm、ΔQmRespectively injecting active power and reactive power into the nodes m of the power transmission and distribution system; m and n are numbers of different nodes in the power transmission and distribution system; thetamnThe phase angle difference of the m and n voltages of the nodes of the power transmission and distribution system is obtained; u shapem、UnThe voltage amplitudes of the nodes of the m and n power transmission and distribution systems are respectively; n belongs to m and represents all branches connected with the power transmission and distribution system node m, and the branch end points are respectively the nodes m and n.
The branch capacity of the power transmission and distribution system is limited to a certain extent, the node voltage needs to be ensured within a certain quality, the active and reactive injection of the power supply is also within a certain range, and the operation constraint expression is as follows:
Figure BDA0002023822630000111
in the formula: sb、
Figure BDA0002023822630000112
The branch capacity, the branch capacity lower limit and the branch capacity upper limit of the power transmission and distribution system are respectively;
Figure BDA0002023822630000113
Figure BDA0002023822630000114
respectively injecting a lower limit and an upper limit of active power into a node m of the power transmission and distribution system;
Figure BDA0002023822630000115
respectively injecting a lower limit and an upper limit of reactive power into a node m of the power transmission and distribution system;
Figure BDA0002023822630000116
the lower voltage limit and the upper voltage limit of the node m of the power transmission and distribution system are respectively.
(3) Park operation optimization method
A. Method for operating gas boiler
The gas boiler is only used as a heat supplementing device in an industrial park comprehensive energy system, and is started to work only when the power output of a thermal power plant cannot meet the requirements of thermal users independently and completely, so that the peak shaving effect is achieved, and the expression of the operation method is as follows:
Figure BDA0002023822630000117
in the formula: sigma Hload(t) is the total thermal power required by the head end supply in the distributed heat system at time t.
B. Energy storage system operation method
Energy storage system among the industrial park comprehensive energy system mainly is for the running cost of reduction system, through charging at the valley moment, the peak moment discharges and reduces system economy running cost, and on the other hand, at the flat moment, energy storage system does not start, can effectual reduction charge-discharge number of times, extension energy storage equipment life, and its operation method expression is:
Figure BDA0002023822630000118
in the formula: t is tp、tt、tvThe time-of-use electricity price is the peak electricity price time interval, the flat electricity price time interval and the valley electricity price time interval of the time-of-use electricity price.
C. Economic operation optimization method
On the basis of a gas boiler operation method and an energy storage system operation method, the economic efficiency of the industrial park comprehensive energy system is mainly concerned, so that the economic cost of the industrial park comprehensive energy system is minimized. The economic cost of the industrial park comprehensive energy system mainly comprises the operation cost of a thermal power plant, the natural gas consumption cost of a gas boiler and the electricity purchasing cost from a large power grid, and the expression of the economic operation optimization method is as follows:
Figure BDA0002023822630000121
in the formula:
Figure BDA0002023822630000122
the total economic cost of the industrial park comprehensive energy system in the optimized operation period T is achieved; xigasA price to purchase natural gas for an industrial park energy system; xigrid(t) the time-of-use electricity price of the industrial park comprehensive energy system for purchasing electricity from the large power grid at the moment t;
Figure BDA0002023822630000123
in order to purchase power from the large grid at time t.
(4) Outputting park system information
And outputting park system information, which comprises information such as the heat output of a gas boiler, the electric heat output of a thermal power plant, the charge and discharge power of an energy storage system, the electric output of photovoltaic, the temperature of hot water in pipelines of a transmission and distribution heat system, the pipeline flow of the transmission and distribution heat system, the node voltage of the transmission and distribution system, the electricity purchasing quantity, the gas consumption quantity, the economic cost and the like.
(5) Example analysis
A. Introduction to the examples
The novel industrial park comprehensive energy system in the invention takes a typical winter day as a research object, the simulation step length is 1 hour, and the simulation period is 24 hours. The structure of the comprehensive energy system of the novel industrial park in the embodiment is shown in figure 2: the main equipment comprises a gas boiler, a thermal power plant, an energy storage system and a photovoltaic system; the main transmission and distribution system comprises a transmission and distribution heat system with 11 nodes and a transmission and distribution power system with 13 nodes. The gas boiler generates heat by consuming natural gas of the gas storage tank, and the output thermal power is injected into a node 1 of the transmission and distribution heat system; the thermal power plant generates electricity and supplies heat by consuming coal, the output heat energy is injected into a node 1 of a transmission and distribution heat system, and the output electric energy is injected into a node 13 of the transmission and distribution system; the energy storage system is connected with a large power grid, charging electric energy comes from the large power grid, and discharging electric energy is injected into a node 13 of the transmission and distribution system; the large power grid is also connected with a node 13 of the power transmission and distribution system; the photovoltaic is connected with a node 11 of the power transmission and distribution system; the heat transmission and distribution system has 13 nodes, wherein the node 1 is a head end supply source, the nodes 2 to 5 are pipeline connection cross points, and the nodes 6 to 11 are heat users; the power transmission and distribution system has 13 nodes in total, wherein the nodes 1 to 10 are electric load nodes, the node 12 is an electric load node, the node 11 is a photovoltaic power supply node, and the node 13 is a balance node; gas flow is formed between the gas storage tank and the gas boiler; coal flow between the coal and the thermal power plant; the hot water input into the transmission and distribution heat system and the transmission and distribution heat system is a heat energy flow; the electric energy input into the power transmission and distribution system and the power transmission and distribution system is electric energy flow; when the thermal output of the thermal power plant is not enough to meet the thermal load requirement, the gas-fired boiler is started to operate, otherwise, the gas-fired boiler is in a shutdown state; the energy storage system allows discharging at the peak value of the large power grid, allows charging at the valley value of the large power grid, and does not work at the average value of the large power grid.
The invention relates to a method for setting main parameters of a novel industrial park comprehensive energy system example, which comprises the following steps: the thermal performance efficiency of the gas boiler was 0.9; the energy loss rate of the energy storage system is 0.01; the natural gas price is 300 CNY/MWh; the maximum electric output of the thermal power plant is 21.1MW/h, and the maximum thermal output is 24.1 MW/h; the capacity of the energy storage system is 4 MWh; the voltage at node 13 is 1.1 pu; the voltage at node 11 is 1.05 pu; the time of use electricity price is shown in fig. 3.
B. Analysis of results
And compiling a model program based on a LINGO18.0 software platform and calling a global solver to solve the established economic operation optimization model of the novel industrial park comprehensive energy system. The electric power injected into the power transmission and distribution system node 13 and the electric output power of the photovoltaic are changed as shown in fig. 4, the electric power injected into the node 13 exceeds the maximum electric output of the thermal power plant in some time intervals, so the power transmission and distribution system inevitably discharges to the large power grid or the energy storage system to meet the electric power demand of the power transmission and distribution system, and as can be seen from fig. 4, the electric output of the photovoltaic is mainly distributed between the time interval 7 and the time interval 20 and bears a part of the electric load power of the power transmission and distribution system, therefore, the electric power injected into the node 13 is low in the period and is in negative correlation with the electric output of the photovoltaic; fig. 5 shows the variation of the heat power injected from the supply source at the head end of the heat distribution and transmission system, and it can also be found from fig. 5 that the heat power demand of the heat distribution and transmission system exceeds the maximum heat output of the thermal power plant in a part of the time period, which inevitably causes the gas-fired boiler to start to supplement the heat energy, and the gas-fired boiler plays a role of peak regulation; the heat power balance process of the comprehensive energy system of the novel industrial park is shown in fig. 6, and it can be found from fig. 6 that the thermal power plant bears the heat power requirement of the main heat transmission and distribution system, during the period from 9 to 17, the gas-fired boiler plays a role of heat compensation, although the heat output of the thermal power plant still has a little margin during the period, the heat output is indirectly limited due to the limitation of the climbing rate of the electric output, meanwhile, it can be found from fig. 4 that the electric output of the photovoltaic is mainly distributed between 7 and 20, the energy storage system is also in a state of allowing power generation, so that the electric heat output level of the thermal power plant is lower, so as to absorb the solar energy, and the electric heat output of the thermal power plant is relatively less; the distribution of the electric heating output of the thermal power plant during the economic optimization operation is shown in fig. 7, and according to the result of fig. 7, the distribution of the electric heating output is relatively concentrated on the upper right corner and the lower right corner, the adjustment of the thermal power plant is relatively stable, and the electric heating output is in a point safety range; the cost result of the economic optimization operation is shown in fig. 8, and it can be known from fig. 8 that the operating cost of the thermal power plant in the integrated energy system of the novel industrial park is the highest, which reflects that the thermal power plant is the core unit of the park, and the efficiency of the thermal power plant is improved, so that the operating cost of the system can be further reduced. The result of fig. 8 shows that the new industrial park comprehensive energy system mainly purchases electricity during the period from 1 to 7 and the period from 19 to 24, and during this period, the electricity price is mainly at the valley value, and the electricity purchasing cost is relatively low, so as to better meet the actual requirements of the new industrial park comprehensive energy system. From the above analysis, the model established by the invention and the proposed economic operation optimization method are reasonable and effective.

Claims (5)

1. An economic operation optimization method for an industrial park comprehensive energy system is characterized by comprising the following steps: the method comprises the following steps:
(1) collecting park system data and information, wherein the park system data and the information comprise heat load demand data, electric load demand data, gas boiler capacity, thermal power plant capacity, energy storage system capacity, photovoltaic capacity, topological information of a heat transmission and distribution system, topological information of a power transmission and distribution system and equipment coupling relation information;
(2) establishing an operation model and setting operation constraints, wherein the operation model comprises a park equipment model and operation constraints thereof, and a park transmission and distribution system model and operation constraints thereof, the park equipment model comprises a thermal power plant model, a gas boiler model, an energy storage system model and a photovoltaic equipment model, and the park transmission and distribution system model comprises a park transmission and distribution thermal system model and a transmission and distribution system model;
(3) optimizing a park operation strategy, which comprises gas boiler operation optimization, energy storage system operation optimization and economic operation optimization;
(4) outputting park system information, which comprises the thermal output of a gas boiler, the electric thermal output of a thermal power plant, the charge and discharge power of an energy storage system, the photovoltaic power output, the hot water temperature in pipelines of a transmission and distribution thermal system, the pipeline flow of the transmission and distribution thermal system, the node voltage of the transmission and distribution system, the electricity purchasing quantity, the gas consumption quantity and the economic cost;
the park equipment model and the operation constraint thereof are specifically as follows:
the expression of the thermal power plant model in the park equipment model is as follows:
Figure FDA0002891273890000011
in the formula:
Figure FDA0002891273890000012
the coal consumption cost of the thermal power plant at the moment t; pCHP(t) is the electrical output of the thermal power plant at time t; hCHP(t) is the thermal output of the thermal power plant at time t; a isCHP、bCHP、cCHP、dCHP、eCHP、fCHPThe cost characteristic parameter of the thermal power plant;
the thermal power plant operation constraint expression in the park equipment operation constraint is as follows:
Figure FDA0002891273890000013
in the formula: alpha is alphai、βi、χiThe method comprises the following steps of (1) setting an i-th inequality constraint parameter of an electrothermal output operation area of the thermal power plant, wherein i is 1,2,3 and 4;
Figure FDA0002891273890000014
Figure FDA0002891273890000015
respectively representing the lower limit and the upper limit of the thermal output of the thermal power plant;
Figure FDA0002891273890000016
respectively representing the lower limit and the upper limit of the electric output of the thermal power plant;
the expression of the gas boiler model in the park equipment model is as follows:
Figure FDA0002891273890000017
in the formula: hgas(t) is the thermal output of the gas boiler at time t; etagasPerformance efficiency for gas-fired boilers;
Figure FDA0002891273890000021
the gas consumption of the gas boiler at the moment t;
the operation constraint expression of the gas-fired boiler in the park equipment operation constraint is as follows:
Figure FDA0002891273890000022
in the formula:
Figure FDA0002891273890000023
respectively is the lower limit and the upper limit of the thermal output of the gas boiler;
the energy storage system model expression in the park equipment model is as follows:
Figure FDA0002891273890000024
in the formula: wbat(t) is the stored energy of the energy storage system at time t; deltabatIs an energy loss factor of the energy storage system; delta t is the economic optimization operation simulation step length; pbat(t) is the charge-discharge power quantity of the energy storage system at the moment t, wherein the positive value is charge, and the negative value is discharge;
Figure FDA0002891273890000025
the charge quantity of the energy storage system at the time t;
Figure FDA0002891273890000026
the discharge capacity of the energy storage system at the moment t;
the energy storage system operation constraint expression in the park equipment operation constraint is as follows:
Figure FDA0002891273890000027
in the formula:
Figure FDA0002891273890000028
respectively representing the lower limit and the upper limit of the stored energy of the energy storage system;
Figure FDA0002891273890000029
respectively is the lower limit and the upper limit of the charging power of the energy storage system;
Figure FDA00028912738900000210
respectively representing the lower limit and the upper limit of the discharge power of the energy storage system; t is an economic optimization operation simulation period;
the photovoltaic model expression in the park equipment model is as follows:
Figure FDA00028912738900000211
in the formula：
Figure FDA00028912738900000212
Outputting electric power for maximum prediction of photovoltaic at the moment t; pstThe maximum test power generation power of the photovoltaic under the standard test condition; gpv(t)、GstThe actual solar radiation intensity of the photovoltaic at the moment t and the solar radiation intensity under the standard test condition are shown; epsilonTAdjusting a factor for the photovoltaic power temperature; t istem(t) is the actual operating temperature of the photovoltaic at time t;
Figure FDA00028912738900000213
is the photovoltaic surface temperature under standard test conditions;
the photovoltaic operation constraint expression in the park equipment operation constraint is as follows:
Figure FDA0002891273890000031
in the formula:
Figure FDA0002891273890000032
the actual output electric power of the photovoltaic at the time t.
2. The method of claim 1 for optimizing the economic operation of the industrial park integrated energy system, wherein the method comprises the following steps: the park transmission and distribution system model and the operation constraint thereof are specifically as follows:
the expression of the heat transmission and distribution system model in the garden transmission and distribution system model is as follows:
Figure FDA0002891273890000033
in the formula:
Figure FDA0002891273890000034
respectively for users in the heat transmission and distribution systemThe required thermal power and the thermal power of the head end supply source;
Figure FDA0002891273890000035
Figure FDA0002891273890000036
respectively the hot water flow flowing through the hot user and the head end supply source in the transmission and distribution heat system;
Figure FDA0002891273890000037
the temperature of the water supply flowing into the heat user and the temperature of the hot water flowing out of the heat user in the transmission and distribution heat system are respectively set;
Figure FDA0002891273890000038
the water supply temperature and the water inlet temperature of the head end supply source are respectively; cpThe specific heat capacity parameter of the hot water of the transmission and distribution heat system;
Figure FDA0002891273890000039
respectively the hot water temperature at the head end and the terminal end of a pipeline in the transmission and distribution heat system;
Figure FDA00028912738900000310
the temperature of the environment where the heat transmission and distribution system is located; lambda [ alpha ]pIs the heat transfer coefficient of the pipe; l ispIs the length of the pipe; m ispThe flow rate of hot water in the pipeline;
the operation constraint expression of the heat transmission and distribution system in the operation constraint of the park transmission and distribution system is as follows:
Figure FDA00028912738900000311
in the formula:
Figure FDA00028912738900000312
respectively the lower limit and the upper limit of the pipeline flow in the transmission and distribution heat system;
Figure FDA00028912738900000313
the lower limit and the upper limit of the temperature of hot water flowing out of a heat user of the transmission and distribution heat system are respectively set;
Figure FDA00028912738900000314
the lower limit and the upper limit of the water supply temperature of a head end supply source of the transmission and distribution heat system are respectively set; m isout、ToutRespectively the hot water flow and the hot water temperature of the nodes of the outflow heat distribution system; m isin、TinRespectively the hot water flow and the hot water temperature flowing into the nodes of the heat transmission and distribution system;
the power transmission and distribution system model expression in the park power transmission and distribution system model is as follows:
Figure FDA0002891273890000041
in the formula: gmn、BmnRespectively the conductance and susceptance of the branch mn; delta Pm、ΔQmRespectively injecting active power and reactive power into the nodes m of the power transmission and distribution system; m and n are numbers of different nodes in the power transmission and distribution system; thetamnThe phase angle difference of the m and n voltages of the nodes of the power transmission and distribution system is obtained; u shapem、UnThe voltage amplitudes of the nodes of the m and n power transmission and distribution systems are respectively; n belongs to m and represents all branches connected with the nodes m of the power transmission and distribution system, and the end points of the branches are the nodes m and n respectively;
the operation constraint expression of the power transmission and distribution system in the operation constraint of the park power transmission and distribution system is as follows:
Figure FDA0002891273890000042
in the formula: sb、
Figure FDA0002891273890000043
The branch capacity, the branch capacity lower limit and the branch capacity upper limit of the power transmission and distribution system are respectively;
Figure FDA0002891273890000044
Figure FDA0002891273890000045
respectively injecting a lower limit and an upper limit of active power into a node m of the power transmission and distribution system;
Figure FDA0002891273890000046
respectively injecting a lower limit and an upper limit of reactive power into a node m of the power transmission and distribution system;
Figure FDA0002891273890000047
the lower voltage limit and the upper voltage limit of the node m of the power transmission and distribution system are respectively.
3. The method of claim 1 for optimizing the economic operation of the industrial park integrated energy system, wherein the method comprises the following steps: the corresponding mathematical expression of the operation optimization of the gas boiler is as follows:
Figure FDA0002891273890000048
in the formula: sigma Hload(t) is the total thermal power required by the head end supply in the distributed heat system at time t.
4. The method of claim 1 for optimizing the economic operation of the industrial park integrated energy system, wherein the method comprises the following steps: the corresponding mathematical expression of the operation optimization of the energy storage system is as follows:
Figure FDA0002891273890000051
in the formula: t is tp、tt、tvThe time-of-use electricity price is the peak electricity price time interval, the flat electricity price time interval and the valley electricity price time interval of the time-of-use electricity price.
5. The method of claim 1 for optimizing the economic operation of the industrial park integrated energy system, wherein the method comprises the following steps: the corresponding function expression of the economic operation optimization is as follows:
Figure FDA0002891273890000052
in the formula:
Figure FDA0002891273890000053
the total economic cost of the industrial park comprehensive energy system in the optimized operation period T is achieved; xigasA price to purchase natural gas for an industrial park energy system; xigrid(t) the time-of-use electricity price of the industrial park comprehensive energy system for purchasing electricity from the large power grid at the moment t;
Figure FDA0002891273890000054
in order to purchase power from the large grid at time 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 CN110070216A (en) 2019-07-30
CN110070216B true 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)

Families Citing this family (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111864805B (en) * 2020-08-26 2021-11-09 合肥工业大学 Energy management method of industrial park micro-grid light storage combined power generation device
CN112966873B (en) * 2021-03-18 2023-01-24 南方电网科学研究院有限责任公司 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
Economic optimization for configuration and sizing of micro integrated energy systems;Haibo YU et al.;《J. Mod. Power Syst. Clean Energy》;20170512;全文 *
考虑运行约束的区域电力–天然气–热力综合能源系统能量流优化分析;王伟亮 等;《中国电机工程学报》;20171220;全文 *

Also Published As

Publication number Publication date
CN110070216A (en) 2019-07-30

Similar Documents

Publication Publication Date Title
CN109193815B (en) Combined heat and power dispatching method for improving wind power consumption
CN107204632A (en) A kind of flexible load Optimization Scheduling for lifting wind electricity digestion
CN110991000B (en) Modeling method for energy hub considering solid oxide fuel cell and electric conversion gas
CN110676847A (en) Optimal scheduling method considering wind power-heat storage unit-electric boiler combined operation
CN110163415B (en) Multi-objective fuzzy cooperative optimization method for multi-energy flow system under variable working condition characteristic
CN112836882B (en) Regional comprehensive energy system operation optimization method considering equipment load rate change
CN110543966A (en) day-ahead scheduling optimization method for micro-energy grid with combined supply of electricity, heat and gas
CN110070216B (en) Economic operation optimization method for industrial park comprehensive energy system
CN115238987A (en) Energy efficiency improvement scheduling method considering multi-type low-carbon factors and demand response
CN112036652A (en) Photovoltaic-energy storage integrated energy system planning method based on opportunity constraint planning
CN114784797A (en) Thermoelectric optimization day-ahead scheduling method for residential comprehensive energy system considering multiple uncertainties
CN110957722A (en) Day-ahead optimized scheduling method for micro energy grid with electricity-to-gas conversion equipment
Li et al. Planning model of integrated energy system considering P2G and energy storage
CN114386256A (en) Regional electric heating system optimal scheduling method considering flexibility constraint of electric heating equipment and heat supply network characteristics
CN114626624A (en) Comprehensive energy system electric-heat combined optimization scheduling method considering heat load conversion
CN110752610B (en) Electric power regulation method for cogeneration unit containing heat storage tank
CN114255137A (en) Low-carbon comprehensive energy system optimization planning method and system considering clean energy
CN113313305A (en) Non-cooperative game-based comprehensive energy system optimization scheduling method
CN113078689A (en) Load tracking management operation strategy for comprehensive energy optimization configuration
CN113313329A (en) Optimal scheduling method for power distribution network containing comprehensive energy system
CN112202203A (en) Cascade utilization optimization method for energy coupling of electric heating microgrid
Cui et al. Review of Optimal Dispatching of Combined Heat-Power System for Wind Power Accommodation
CN110793093B (en) Electric power control method for cogeneration unit comprising electric boiler and heat storage tank
Zhong et al. Economic optimization operation strategy for industrial park integrated energy system
CN113131513B (en) Method for optimizing operation of electric, thermal and gas conversion system with consideration of carbon emission and storage medium

Legal Events

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