CN114818269A - Industrial park comprehensive energy system and factory production plan collaborative optimization method - Google Patents
Industrial park comprehensive energy system and factory production plan collaborative optimization method Download PDFInfo
- Publication number
- CN114818269A CN114818269A CN202210294463.8A CN202210294463A CN114818269A CN 114818269 A CN114818269 A CN 114818269A CN 202210294463 A CN202210294463 A CN 202210294463A CN 114818269 A CN114818269 A CN 114818269A
- Authority
- CN
- China
- Prior art keywords
- representing
- energy
- production
- workshop
- power
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Pending
Links
- 238000004519 manufacturing process Methods 0.000 title claims abstract description 126
- 238000000034 method Methods 0.000 title claims abstract description 34
- 238000005457 optimization Methods 0.000 title claims abstract description 34
- CURLTUGMZLYLDI-UHFFFAOYSA-N Carbon dioxide Chemical compound O=C=O CURLTUGMZLYLDI-UHFFFAOYSA-N 0.000 claims abstract description 164
- 229910002092 carbon dioxide Inorganic materials 0.000 claims abstract description 82
- 239000001569 carbon dioxide Substances 0.000 claims abstract description 82
- 238000003860 storage Methods 0.000 claims abstract description 77
- 230000008878 coupling Effects 0.000 claims abstract description 58
- 238000010168 coupling process Methods 0.000 claims abstract description 58
- 238000005859 coupling reaction Methods 0.000 claims abstract description 58
- 238000005265 energy consumption Methods 0.000 claims abstract description 51
- 238000006243 chemical reaction Methods 0.000 claims abstract description 35
- UFHFLCQGNIYNRP-UHFFFAOYSA-N Hydrogen Chemical compound [H][H] UFHFLCQGNIYNRP-UHFFFAOYSA-N 0.000 claims abstract description 30
- 239000001257 hydrogen Substances 0.000 claims abstract description 30
- 229910052739 hydrogen Inorganic materials 0.000 claims abstract description 30
- 238000004064 recycling Methods 0.000 claims abstract description 13
- 239000002994 raw material Substances 0.000 claims abstract description 12
- 238000002360 preparation method Methods 0.000 claims abstract description 6
- VNWKTOKETHGBQD-UHFFFAOYSA-N methane Chemical compound C VNWKTOKETHGBQD-UHFFFAOYSA-N 0.000 claims description 66
- 239000007789 gas Substances 0.000 claims description 50
- 239000003345 natural gas Substances 0.000 claims description 33
- 239000013067 intermediate product Substances 0.000 claims description 30
- 230000005611 electricity Effects 0.000 claims description 27
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims description 20
- QVGXLLKOCUKJST-UHFFFAOYSA-N atomic oxygen Chemical compound [O] QVGXLLKOCUKJST-UHFFFAOYSA-N 0.000 claims description 15
- 239000001301 oxygen Substances 0.000 claims description 15
- 229910052760 oxygen Inorganic materials 0.000 claims description 15
- 238000010521 absorption reaction Methods 0.000 claims description 6
- 238000001816 cooling Methods 0.000 claims description 6
- 238000004146 energy storage Methods 0.000 claims description 6
- 230000003993 interaction Effects 0.000 claims description 6
- 238000005868 electrolysis reaction Methods 0.000 claims description 5
- 238000010248 power generation Methods 0.000 claims description 5
- 229910002056 binary alloy Inorganic materials 0.000 claims description 3
- 239000006227 byproduct Substances 0.000 claims description 3
- 230000000694 effects Effects 0.000 claims description 3
- 238000005485 electric heating Methods 0.000 claims description 3
- 238000005516 engineering process Methods 0.000 claims description 3
- 230000020169 heat generation Effects 0.000 claims description 3
- 230000008569 process Effects 0.000 claims description 3
- OKTJSMMVPCPJKN-UHFFFAOYSA-N Carbon Chemical compound [C] OKTJSMMVPCPJKN-UHFFFAOYSA-N 0.000 description 4
- 229910052799 carbon Inorganic materials 0.000 description 4
- 238000010586 diagram Methods 0.000 description 3
- 238000005057 refrigeration Methods 0.000 description 2
- 230000004044 response Effects 0.000 description 2
- 230000008901 benefit Effects 0.000 description 1
- 230000007613 environmental effect Effects 0.000 description 1
- 239000002803 fossil fuel Substances 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q10/00—Administration; Management
- G06Q10/06—Resources, workflows, human or project management; Enterprise or organisation planning; Enterprise or organisation modelling
- G06Q10/063—Operations research, analysis or management
- G06Q10/0631—Resource planning, allocation, distributing or scheduling for enterprises or organisations
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q50/00—Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
- G06Q50/06—Energy or water supply
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2111/00—Details relating to CAD techniques
- G06F2111/04—Constraint-based CAD
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2111/00—Details relating to CAD techniques
- G06F2111/06—Multi-objective optimisation, e.g. Pareto optimisation using simulated annealing [SA], ant colony algorithms or genetic algorithms [GA]
Landscapes
- Engineering & Computer Science (AREA)
- Business, Economics & Management (AREA)
- Human Resources & Organizations (AREA)
- Theoretical Computer Science (AREA)
- Economics (AREA)
- Physics & Mathematics (AREA)
- Strategic Management (AREA)
- General Physics & Mathematics (AREA)
- Tourism & Hospitality (AREA)
- Health & Medical Sciences (AREA)
- Marketing (AREA)
- General Business, Economics & Management (AREA)
- Entrepreneurship & Innovation (AREA)
- Water Supply & Treatment (AREA)
- Game Theory and Decision Science (AREA)
- Educational Administration (AREA)
- Development Economics (AREA)
- Operations Research (AREA)
- Public Health (AREA)
- Quality & Reliability (AREA)
- General Health & Medical Sciences (AREA)
- Evolutionary Computation (AREA)
- Computer Hardware Design (AREA)
- Primary Health Care (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
The invention discloses a method for collaborative optimization of an industrial park comprehensive energy system and a factory production plan, which aims at the structures of the industrial park comprehensive energy system and industrial users, and establishes a comprehensive energy system hydrogen preparation model, a carbon dioxide recycling and storage model and a production raw material storage area model on the basis of the original energy supply, conversion and storage equipment model; constructing a multi-energy coupling model and a factory production constraint model of different workshop energy consumption of an industrial user; the method comprises the steps of constructing a collaborative optimization model of the operation of the comprehensive energy system and industrial users in the park by taking the operation efficiency of all devices of the comprehensive energy system in the park, a multi-energy coupling model of different workshops of the industrial users and a factory production flow as constraint conditions and taking the total energy consumption cost of one day of the industrial park as an objective function; and solving the collaborative optimization model to obtain the output plan of each device of the park integrated energy system and the production scheduling arrangement of the factory, thereby reducing the total daily operation energy consumption cost of the park integrated energy system and industrial users.
Description
Technical Field
The invention relates to the technical field of energy configuration optimization, in particular to a collaborative optimization method for an industrial park comprehensive energy system and a factory production plan.
Background
With the increasing demand for various energy sources such as electricity, natural gas, heat supply, refrigeration and the like, modern power systems are developing towards comprehensive energy systems. Integrated energy systems are considered an effective model to reduce fossil fuel usage, reduce energy costs, and improve the efficiency and flexibility of traditional energy systems. The energy input comprises natural gas and electric power, and the output type can simultaneously comprise various energy sources such as electric power, heat supply, refrigeration, natural gas and the like. Industrial loads require a large amount of electrical, thermal and cold energy, and thus play an important role in end users; however, the traditional factory operation mode lacks of management and optimization of multiple energy sources, and the economic benefit and the energy efficiency of an enterprise are seriously influenced; therefore, there is a need to develop energy consumption optimization strategies for industrial park integrated energy systems and industrial users.
With increasing concern for environmental issues, much research effort has taken carbon emissions into account in integrated energy systems and various frameworks have been proposed to optimize the operation of integrated energy systems in parks to reduce energy costs; in addition, in demand response, management and optimization of industrial loads are widely used to improve economic efficiency and energy efficiency. However, the research on the energy conversion of the park integrated energy system and the production plan of the factory user for collaborative optimization scheduling is few; and the recycling of the carbon dioxide of the park comprehensive energy system and the coupling relation among different energy sources consumed by factory production are not considered, and the daily production task of an industrial user is not considered to carry out demand response, so that the user participation is seriously influenced.
Disclosure of Invention
The technical problem to be solved by the invention is to provide a method for the collaborative optimization of the industrial park comprehensive energy system and the factory production plan, so that the total daily operation energy consumption cost of the industrial user and the comprehensive energy system in the park is reduced, the comprehensive energy utilization rate is improved, and the carbon emission is reduced.
In order to solve the technical problems, the technical scheme adopted by the invention is as follows: a method for collaborative optimization of an industrial park comprehensive energy system and a factory production plan comprises the following steps:
step S1, aiming at the structures of the comprehensive energy system and the industrial users in the industrial park, on the basis of the original energy supply, conversion and storage equipment model, a hydrogen preparation model, a carbon dioxide recycling and storage model and a production raw material storage area model of the comprehensive energy system are established;
step S2, the park integrated energy management center receives the production tasks, the workshop energy consumption types and the production constraints of the industrial users, and a multi-energy coupling model and a factory production constraint model of different workshop energy consumptions of the industrial users are constructed;
step S3, the park integrated energy management center integrates the received industrial user information, the equipment information of the integrated energy system, the price and the power information of the external air network and the power grid, and under the premise of finishing the daily production task of the factory, the operation efficiency of all the equipment of the integrated energy system in the park, the multi-energy coupling model of different workshops of the industrial user and the factory production flow are taken as constraint conditions, the total daily energy consumption cost of the industrial park is taken as a target function, and a collaborative optimization model of the operation of the integrated energy system and the industrial user is constructed;
step S4, the park integrated energy management center solves the collaborative optimization model of the integrated energy system and the industrial user operation to obtain the output plan of each device and the factory production scheduling arrangement of the park integrated energy system;
and step S5, the park integrated energy management center controls the operation of the integrated energy system according to the solving result, and sends the obtained result to the factory user to enable the factory user to engage in the production activities according to the corresponding arrangement.
The technical scheme of the invention is further improved as follows: the energy supply equipment model in the step S1 comprises a natural gas network, a power grid, a photovoltaic generator set and a solar water heater; the energy conversion equipment model comprises a gas boiler, a gas turbine, a P2G unit, an absorption refrigerator and an electric refrigerator; the storage device model of energy includes an electrical energy storage device; the comprehensive energy system hydrogen production model comprises an electrolytic water unit; the carbon dioxide recycling and storing model comprises carbon dioxide capturing equipment; the production raw material storage area model comprises an initial raw material warehouse, a workshop 1 intermediate product storage warehouse, a workshop 2 intermediate product storage warehouse and a workshop z intermediate product storage warehouse.
The technical scheme of the invention is further improved as follows: the model of the P2G unit and the electrolyzed water unit is as follows:
in the above formula:representing the rate at which the electric gas conversion unit produces natural gas; δ represents a proportionality coefficient between the rate of hydrogen consumption and the rate of natural gas production;representing the rate at which hydrogen is produced by the electrolysis cell of water;
in the above formula:representing the rate at which carbon dioxide is consumed by the electrical conversion unit; λ represents the proportionality coefficient between the rate of hydrogen consumption and the rate of carbon dioxide consumption for natural gas production;
in the above formula:power representing electrical energy consumed by the electrical conversion unit;representing a conversion factor between a rate of hydrogen consumption and power consumption;
in the above formula:represents the maximum power of the electric energy consumed by the electric gas conversion unit;
in the above formula:representing the rate of oxygen production from the electrolytic cell; τ represents a proportionality coefficient between the rate of hydrogen production and the rate of oxygen production;
in the above formula:representing the rate at which the electric gas conversion unit produces natural gas; θ represents a coefficient for converting the production rate of natural gas into natural gas power.
The technical scheme of the invention is further improved as follows: the carbon dioxide recycling and storing model of the integrated energy system in the step S1 is as follows:
in the above formula:representing the rate of capture of carbon dioxide produced by the cogeneration unit and the gas boiler;indicating the rate of carbon dioxide being discharged into the air due to capture technology limitations; ρ represents a conversion coefficient;representing the heat production power of the gas boiler;representing the heat production efficiency of the gas boiler;respectively representing the power generation and heat production of the cogeneration unit;respectively representing the electricity generation efficiency and the heat generation efficiency of the cogeneration unit;
in the above formula:represents the rate at which the recovered carbon dioxide is used for consumption by the P2G unit;indicating the rate of storage to the carbon dioxide tank.
The technical scheme of the invention is further improved as follows: the carbon dioxide storage model of the integrated energy system in the step S1 is as follows:
in the above formula: CO 2 2s,t Represents the carbon dioxide storage amount;representing an initial storage amount of carbon dioxide;represents the carbon dioxide sales rate;
the above formula shows that the carbon dioxide storage amount is kept unchanged at the beginning and end of carbon dioxide;
0≤CO 2s,t ≤CO 2s,max ,
in the above formula: CO 2 2s,max Represents a maximum storage capacity of carbon dioxide;
The technical scheme of the invention is further improved as follows: the types of the energy consumption of the workshop in the step S2 comprise a pure electric energy workshop, a pure cold energy workshop, a pure heat energy workshop, an electric-heat coupling workshop, an electric-cold coupling workshop and a hot-cold coupling workshop.
The technical scheme of the invention is further improved as follows: in step S2, the multi-energy coupling model of the energy consumption of different workshops of the industrial user is:
in the above formula: p E,t Represents the total power consumed by the plant; z represents a shop mark; e' represents a set of pure electric energy consumption workshops; x is a binary system and represents the running state of the workshop, and the workshop works when 1 is taken out, otherwise the workshop does not work; p E,z Represents the electric power consumed by the z plant; EH' represents a hot spot coupled plant set;represents the electric power consumed by the electro-thermal coupling plant z; EC' represents an electric cold coupling workshop set;represents the electric power consumed by the electric cold coupling workshop z;
in the above formula: p H,t Representing the total power of the plant consuming thermal energy; h' represents a set of pure thermal energy consumption plants; p H,z Represents the thermal power consumed by the z plant;represents the total thermal power consumed by the electro-thermal coupling workshop; HC' represents a set of hot and cold coupled plants;represents the thermal power consumed by the hot-cold coupling workshop z;
in the above formula: p C,t Represents the total power of the plant consuming cooling energy; c' represents a pure cold energy consumption workshop set;representing the total cold power consumed by the electric cold coupling workshop;representing the total cold power consumed by the hot-cold coupling workshop;
wherein:
in the above formula:representing the z workshop electrical thermal coupling coefficient;representing the z workshop electric cold coupling coefficient;represents the z workshop heat-cold coupling coefficient.
The technical scheme of the invention is further improved as follows: the production constraints in the step S2 include the production sequence, the production speed, the operation state of different workshops, and the initial number and the remaining capacity of the storage warehouses;
P (EHC),z representing various coupling energy consumption powers of a workshop z; υ ═ υ 1 …υ Z Indicates the production rates of the different plants;an initial storage amount of intermediate products representing the plant z;representing a plant z intermediate product storage capacity;
specifically, the intermediate product storage warehouse model:
in the above formula:representing the number of intermediate product stores in shop z;representing the initial storage amount of the intermediate products in the workshop z; upsilon is z,t Representing the processing rate of the workshop z at different moments;
the above expression shows that the storage capacity of the intermediate product is kept unchanged at the beginning and the end;
in the above formula:representing the maximum storage capacity of the intermediate product in the workshop z;
υ z,t =x z,t ·υ z ,
in the above formula: upsilon is z Representing the plant z processing rate.
The technical scheme of the invention is further improved as follows: the specific process of step S3 is as follows:
on the premise of finishing daily production tasks of a factory, taking total energy consumption cost of one day in an industrial park as a target function, specifically comprising electricity purchasing cost to a power grid, gas purchasing cost to a natural gas station, carbon dioxide emission cost, income of carbon dioxide sale capture and income of oxygen sale as a byproduct of hydrogen production by electrolyzing water;
specifically, the objective function expression of the total energy consumption cost of the industrial park in one day is as follows:
in the above formula: c total Represents the total energy consumption cost of one day of the industrial park; c t Represents the total cost of energy consumption at a certain time in one day of the industrial park;
wherein:
in the above formula:representing the power purchased by the power grid; p is a radical of E,t Representing the electricity purchase price of the power grid;representing the gas purchasing power of the natural gas station; p is a radical of g,t Representing the gas purchase price of the natural gas station;represents the carbon dioxide emission rate;represents a carbon dioxide emission price;represents a carbon dioxide purchase rate;represents a carbon dioxide purchase price;representing the rate of oxygen production from the electrolysis of water to produce hydrogen;represents the price for oxygen sale;represents the rate at which carbon dioxide is sold;represents the price for selling carbon dioxide;
the constraint conditions include: the method comprises the following steps of electric heating and cooling power balance constraint, power grid and gas grid interaction power constraint, operation constraint of each energy device, energy storage constraint and factory production constraint;
the relation expression among electricity, heat and cold produced by the comprehensive energy system is as follows:
in the above formula:representing the power of electricity, heat and cold provided by the park comprehensive energy system; represents the charge and discharge power of the storage battery;representing the photovoltaic power generation;representing the electricity converted into the electricity consumption power;indicating electricityGas conversion and gas production power;the heat production power of the solar water heater is represented;
wherein:
in the above formula: ω represents a percentage;representing the heat production efficiency and the electricity production efficiency of the cogeneration unit;representing gas boiler efficiency; COP EC Representing the conversion coefficient of the electric refrigerator; COP AC Representing the conversion coefficient of the absorption refrigerator;
specifically, the balance constraints of the comprehensive energy system and the electricity, heat and cold powers of the industrial users in the park are as follows:
specifically, the power grid and gas grid interaction power constraint is as follows:
Due to the adoption of the technical scheme, the invention has the technical progress that:
the invention provides a method for collaborative optimization of an industrial park comprehensive energy system and a factory production plan, which aims at the structures of the industrial park comprehensive energy system and industrial users, and establishes a comprehensive energy system hydrogen preparation model, a carbon dioxide recycling and storage model and a production raw material storage area model on the basis of the original energy supply, conversion and storage equipment model; constructing a multi-energy coupling model and a factory production constraint model of different workshop energy consumption of an industrial user; on the premise of finishing daily production tasks of a factory, constructing a collaborative optimization model of the operation of the comprehensive energy system and industrial users in the park by taking the operation efficiency of all devices of the comprehensive energy system in the park, a multi-energy coupling model of different workshops of the industrial users and the production flow of the factory as constraint conditions and taking the total energy consumption cost of one day of the industrial park as an objective function; and solving the established collaborative optimization model to obtain the output plan of each device of the park integrated energy system and the production scheduling arrangement of the factory, so that the total daily operation energy consumption cost of the park integrated energy system and industrial users can be reduced, the energy integrated utilization rate is improved, and the carbon emission is reduced.
Drawings
FIG. 1 is a block diagram of the collaborative optimization steps for operation of an integrated energy system and industrial users on a campus in accordance with the present invention;
FIG. 2 is a system block diagram of a collaborative optimization model of the operation of the integrated energy system and the industrial users on the campus in accordance with the present invention;
figure 3 is a schematic diagram of the production of the various plants of the plant of the present invention on a campus.
Detailed Description
The present invention will be described in further detail with reference to the following examples:
as shown in fig. 1, a method for collaborative optimization of an industrial park integrated energy system and a factory production plan includes the following steps:
and step S1, aiming at the structures of the comprehensive energy system and the industrial users in the industrial park, establishing a hydrogen preparation model, a carbon dioxide recycling and storing model and a production raw material storage area model of the comprehensive energy system on the basis of the original energy supply, conversion and storage equipment models. As shown in fig. 2, the relation between the supply and demand of the electric energy, the heat energy and the cold energy produced by the park integrated energy system and the consumption of the electric energy, the heat energy and the cold energy in different workshops of the factory is described, and the optimization of the park integrated energy system and the production scheduling arrangement of the factory are realized through the cooperative optimization between the supply and demand, so that the total daily production energy consumption cost of the park integrated energy system and the factory is reduced; the supply equipment model of the energy of the comprehensive energy system comprises a natural gas network, a power grid, a photovoltaic generator set and a solar water heater; the energy conversion equipment model comprises a gas boiler, a gas turbine, a P2G unit, an absorption refrigerator and an electric refrigerator; the storage device model of energy includes an electrical energy storage device; the comprehensive energy system hydrogen production model comprises an electrolytic water unit; the carbon dioxide recycling and storing model comprises carbon dioxide capturing equipment; the factory is provided with N production workshops, and the production raw material storage area model comprises an initial raw material warehouse, a workshop 1 intermediate product storage warehouse, a workshop 2 intermediate product storage warehouse and a workshop z intermediate product storage warehouse.
The model of the P2G cell and the electrolyzed water cell is as follows:
in the above formula:representing the rate at which the electric gas conversion unit produces natural gas; δ represents a proportionality coefficient between the rate of hydrogen consumption and the rate of natural gas production;representing the rate at which hydrogen is produced by the electrolysis cell of water;
in the above formula:representing the rate at which carbon dioxide is consumed by the electrical conversion unit; λ represents the proportionality coefficient between the rate of hydrogen consumption and the rate of carbon dioxide consumption for natural gas production;
in the above formula:power representing electrical energy consumed by the electrical conversion unit;representing a conversion factor between a rate of hydrogen consumption and power consumption;
in the above formula:represents the maximum power of the electric energy consumed by the electric gas conversion unit;
in the above formula:representing the rate of oxygen production from the electrolytic cell; τ represents the rate of hydrogen production and oxygen productionA proportionality coefficient between rates of;
in the above formula:representing the rate at which the electric gas conversion unit produces natural gas; θ represents a coefficient for converting the production rate of natural gas into natural gas power.
The comprehensive energy system carbon dioxide recycling model comprises the following steps:
in the above formula:representing the rate of capture of carbon dioxide produced by the cogeneration unit and the gas boiler;indicating the rate of carbon dioxide being discharged into the air due to capture technology limitations; ρ represents a conversion coefficient;representing the heat production power of the gas boiler;representing the heat production efficiency of the gas boiler;respectively representing the power generation and heat production of the cogeneration unit;respectively representing the electricity generation efficiency and the heat generation efficiency of the cogeneration unit;
in the above formula:represents the rate at which the recovered carbon dioxide is used for consumption by the P2G unit;indicating the rate of storage to the carbon dioxide tank.
in the above formula: CO 2 2s,t Represents the carbon dioxide storage amount;representing an initial storage amount of carbon dioxide;represents the carbon dioxide sales rate;
the above formula shows that the carbon dioxide storage amount is kept unchanged at the beginning and end of carbon dioxide;
0≤CO 2s,t ≤CO 2s,max ,
in the above formula: CO 2 2s,max Represents a maximum storage capacity of carbon dioxide;
And S2, the park comprehensive energy management center receives the production tasks, the workshop energy consumption types and the production constraints of the industrial users, and a multi-energy coupling model and a factory production constraint model of different workshop energy consumptions of the industrial users are constructed, wherein the workshop energy consumption types comprise a pure electric energy workshop, a pure cold energy workshop, a pure heat energy workshop, an electric-heat coupling workshop, an electric-cold coupling workshop and a hot-cold coupling workshop. The multi-energy coupling model of the energy consumption of different workshops of the industrial user is as follows:
in the above formula: p E,t Represents the total power consumed by the plant; z represents a shop mark; e' represents a set of pure electric energy consumption plants; x is a binary system and represents the running state of a workshop, and the workshop works when 1 is taken out, otherwise the workshop does not work; p E,z Represents the electric power consumed by the z plant; EH' represents a hot spot coupled plant set;represents the electric power consumed by the electro-thermal coupling plant z; EC' represents an electric cold coupling workshop set;represents the electric power consumed by the electric cold coupling workshop z;
in the above formula: p H,t Representing the total power of the plant consuming thermal energy; h' represents a set of pure thermal energy consumption plants; p H,z Represents the thermal power consumed by the z plant;represents the total thermal power consumed by the electro-thermal coupling workshop; HC' represents a set of hot and cold coupled plants;represents the thermal power consumed by the hot-cold coupling workshop z;
in the above formula: p C,t Represents the total power of the plant consuming cooling energy; c' represents a pure cold energy consumption workshop set;representing the total cold power consumed by the electric cold coupling workshop;representing the total cold power consumed by the hot-cold coupling workshop;
wherein:
in the above formula:representing the z workshop electrical thermal coupling coefficient;representing the z workshop electric cold coupling coefficient;represents the z workshop heat-cold coupling coefficient.
The production constraints comprise the production sequence, the production speed and the running state of different workshops, and the initial number and the residual capacity of the storage warehouses;
as shown in FIG. 3, P (EHC),z Representing various coupling energy consumption powers of a workshop z; υ ═ υ 1 …υ Z Indicates the production rates of the different plants;an initial storage amount of intermediate products representing the plant z;representing a plant z intermediate product storage capacity;
specifically, the intermediate product storage warehouse model:
in the above formula:representing the number of intermediate product stores in shop z;representing the initial storage amount of the intermediate products in the workshop z; upsilon is z,t Representing the processing rate of the workshop z at different moments;
the above expression shows that the storage capacity of the intermediate product is kept unchanged at the beginning and the end;
in the above formula:representing the maximum storage capacity of the intermediate product in the workshop z;
υ z,t =x z,t ·υ z ,
in the above formula: upsilon is z Representing the plant z processing rate.
Step S3, the park integrated energy management center integrates the received industrial user information, the equipment information of the integrated energy system, the price and the power information of the external air network and the power grid, and under the premise of finishing the daily production task of the factory, the operation efficiency of all the equipment of the integrated energy system in the park, the multi-energy coupling model of different workshops of the industrial user and the factory production flow are taken as constraint conditions, the total daily energy consumption cost of the industrial park is taken as a target function, and a collaborative optimization model of the operation of the integrated energy system and the industrial user is constructed; the specific process is as follows:
on the premise of finishing daily production tasks of a factory, taking total energy consumption cost of one day in an industrial park as a target function, specifically comprising electricity purchasing cost to a power grid, gas purchasing cost to a natural gas station, carbon dioxide emission cost, income of carbon dioxide sale capture and income of oxygen sale as a byproduct of hydrogen production by electrolyzing water;
specifically, the objective function expression of the total energy consumption cost of the industrial park in one day is as follows:
in the above formula: c total Represents the total energy consumption cost of one day of the industrial park; c t The total energy consumption cost at a certain time in one day of the industrial park is represented;
wherein:
in the above formula:representing the power purchased by the power grid; p is a radical of E,t Representing the electricity purchase price of the power grid;representing the gas purchasing power of the natural gas station; p is a radical of g,t Representing the gas purchase price of the natural gas station;represents the carbon dioxide emission rate;represents a carbon dioxide emission price;representing a carbon dioxide purchase rate;represents a carbon dioxide purchase price;representing the rate of oxygen production from the electrolysis of water to produce hydrogen;represents the price for oxygen sale;represents the rate at which carbon dioxide is sold;represents the price for selling carbon dioxide;
the constraint conditions include: the method comprises the following steps of electric heating and cooling power balance constraint, power grid and gas grid interaction power constraint, operation constraint of each energy device, energy storage constraint and factory production constraint;
the relation expression among the electricity, heat and cold produced by the comprehensive energy system is as follows:
in the above formula:representing the power of electricity, heat and cold provided by the park comprehensive energy system; represents the charge and discharge power of the storage battery;representing the photovoltaic power generation;representing electricity to gas consumption power;the power of the electricity-to-gas production is shown;the heat production power of the solar water heater is represented;
wherein:
in the above formula: ω represents a percentage;representing the heat production efficiency and the electricity production efficiency of the cogeneration unit;representing gas boiler efficiency; COP EC Representing the conversion coefficient of the electric refrigerator; COP AC Representing the conversion coefficient of the absorption refrigerator;
specifically, the balance constraints of the comprehensive energy system and the electricity, heat and cold powers of the industrial users in the park are as follows:
specifically, the power grid and gas grid interaction power constraint is as follows:
Step S4, the park integrated energy management center solves the collaborative optimization model of the integrated energy system and the industrial user operation to obtain the output plan of each device and the factory production scheduling arrangement of the park integrated energy system;
and step S5, the park integrated energy management center controls the operation of the integrated energy system according to the solving result, and sends the obtained result to the factory user to enable the factory user to engage in the production activities according to the corresponding arrangement.
The invention provides a method for collaborative optimization of an industrial park comprehensive energy system and a factory production plan, which aims at the structures of the industrial park comprehensive energy system and industrial users, and establishes a comprehensive energy system hydrogen preparation model, a carbon dioxide recycling and storage model and a production raw material storage area model on the basis of the original energy supply, conversion and storage equipment model; constructing a multi-energy coupling model and a factory production constraint model of different workshop energy consumption of an industrial user; on the premise of finishing daily production tasks of a factory, constructing a collaborative optimization model of the operation of the comprehensive energy system and industrial users in the park by taking the operation efficiency of all devices of the comprehensive energy system in the park, a multi-energy coupling model of different workshops of the industrial users and the production flow of the factory as constraint conditions and taking the total energy consumption cost of one day of the industrial park as an objective function; and solving the established collaborative optimization model to obtain the output plan of each device of the park integrated energy system and the production scheduling arrangement of the factory, so that the total daily operation energy consumption cost of the park integrated energy system and industrial users can be reduced, the energy integrated utilization rate is improved, and the carbon emission is reduced.
Claims (9)
1. A method for collaborative optimization of an industrial park comprehensive energy system and a factory production plan is characterized by comprising the following steps: the method comprises the following steps:
step S1, aiming at the structures of the comprehensive energy system and the industrial users in the industrial park, on the basis of the original energy supply, conversion and storage equipment model, a hydrogen preparation model, a carbon dioxide recycling and storage model and a production raw material storage area model of the comprehensive energy system are established;
step S2, the park integrated energy management center receives the production tasks, the workshop energy consumption types and the production constraints of the industrial users, and a multi-energy coupling model and a factory production constraint model of different workshop energy consumptions of the industrial users are constructed;
step S3, the park integrated energy management center integrates the received industrial user information, the equipment information of the integrated energy system, the price and the power information of the external air network and the power grid, and under the premise of finishing the daily production task of the factory, the operation efficiency of all the equipment of the integrated energy system in the park, the multi-energy coupling model of different workshops of the industrial user and the factory production flow are taken as constraint conditions, the total daily energy consumption cost of the industrial park is taken as a target function, and a collaborative optimization model of the operation of the integrated energy system and the industrial user is constructed;
step S4, the park integrated energy management center solves the collaborative optimization model of the integrated energy system and the industrial user operation to obtain the output plan of each device and the factory production scheduling arrangement of the park integrated energy system;
and step S5, the park integrated energy management center controls the operation of the integrated energy system according to the solving result, and sends the obtained result to the factory user to enable the factory user to engage in the production activities according to the corresponding arrangement.
2. The method of claim 1 for collaborative optimization of an industrial park integrated energy system and a factory production plan, wherein the method comprises the following steps: the energy supply equipment model in the step S1 comprises a natural gas network, a power grid, a photovoltaic generator set and a solar water heater; the energy conversion equipment model comprises a gas boiler, a gas turbine, a P2G unit, an absorption refrigerator and an electric refrigerator; the storage device model of energy includes an electrical energy storage device; the comprehensive energy system hydrogen production model comprises an electrolytic water unit; the carbon dioxide recycling and storing model comprises carbon dioxide capturing equipment; the production raw material storage area model comprises an initial raw material warehouse, a workshop 1 intermediate product storage warehouse, a workshop 2 intermediate product storage warehouse and a workshop z intermediate product storage warehouse.
3. The method of claim 2, wherein the method comprises the steps of: the model of the P2G unit and the electrolyzed water unit is as follows:
in the above formula:representing the rate at which the electric gas conversion unit produces natural gas; δ represents a proportionality coefficient between the rate of hydrogen consumption and the rate of natural gas production;indicating waterThe rate at which hydrogen is produced by the electrolyzer of (a);
in the above formula:representing the rate at which carbon dioxide is consumed by the electrical conversion unit; λ represents the proportionality coefficient between the rate of hydrogen consumption and the rate of carbon dioxide consumption for natural gas production;
in the above formula:power representing the electric energy consumed by the electric gas conversion unit;representing a conversion factor between a rate of hydrogen consumption and power consumption;
in the above formula:represents the maximum power of the electric energy consumed by the electric gas conversion unit;
in the above formula:representing the rate of oxygen production from the electrolytic cell; τ represents a proportionality coefficient between the rate of hydrogen production and the rate of oxygen production;
4. The method of claim 3, wherein the method comprises: the carbon dioxide recycling model of the integrated energy system in the step S1 is as follows:
in the above formula:representing the rate of capture of carbon dioxide produced by the cogeneration unit and the gas boiler;indicating the rate of carbon dioxide being discharged into the air due to capture technology limitations; ρ represents a conversion coefficient;representing the heat production power of the gas boiler;representing the heat production efficiency of the gas boiler;respectively representing the power generation and heat production of the cogeneration unit;respectively representing the electricity generation efficiency and the heat generation efficiency of the cogeneration unit;
5. The method of claim 4 for collaborative optimization of an industrial park energy system and a factory production plan, wherein the method comprises the following steps: the carbon dioxide storage model of the integrated energy system in the step S1 is as follows:
in the above formula: CO 2 2s,t Represents the carbon dioxide storage amount;representing an initial storage amount of carbon dioxide;represents the carbon dioxide sales rate;
the above formula shows that the carbon dioxide storage amount is kept unchanged at the beginning and end of carbon dioxide;
0≤CO 2s,t ≤CO 2s,max ,
in the above formula: CO 2 2s,max Represents a maximum storage capacity of carbon dioxide;
6. The method of claim 5 for collaborative optimization of an industrial park energy system and a factory production plan, wherein the method comprises the following steps: the types of the energy consumption of the workshop in the step S2 comprise a pure electric energy workshop, a pure cold energy workshop, a pure heat energy workshop, an electric-heat coupling workshop, an electric-cold coupling workshop and a hot-cold coupling workshop.
7. The method of claim 6, wherein the method comprises: in step S2, the multi-energy coupling model of the energy consumption of different workshops of the industrial user is:
in the above formula: p E,t Represents the total power consumed by the plant; z represents a shop mark; e' represents a set of pure electric energy consumption workshops; x is a binary system and represents the running state of the workshop, and the workshop works when 1 is taken out, otherwise the workshop does not work; p E,z Represents the electric power consumed by the z plant; EH' represents a hot spot coupled plant set;represents the electric power consumed by the electro-thermal coupling plant z; EC' represents an electric cold coupling workshop set;represents the electric power consumed by the electric cold coupling workshop z;
in the above formula: p is H,t Representing the total power of the plant consuming thermal energy; h' represents a set of pure thermal energy consumption plants; p H,z Represents the thermal power consumed by the z plant;represents the total thermal power consumed by the electro-thermal coupling plant; HC' represents a set of hot and cold coupled plants;represents the thermal power consumed by the hot-cold coupling workshop z;
in the above formula: p C,t Represents the total power of the plant consuming cooling energy; c' represents a pure cold energy consumption workshop set;representing the total cold power consumed by the electric cold coupling workshop;representing the total cold power consumed by the hot-cold coupling workshop;
wherein:
8. The method of claim 7, wherein the method comprises: the production constraints in the step S2 include the production sequence, the production speed, the operation state of different workshops, and the initial number and the remaining capacity of the storage warehouses;
P (EHC),z representing various coupling energy consumption powers of a workshop z; υ ═ υ 1 …υ Z Indicates the production rates of the different plants;an initial storage amount of intermediate products representing the plant z;representing a plant z intermediate product storage capacity;
specifically, the intermediate product storage warehouse model:
in the above formula:representing the number of intermediate product stores in shop z;representing the initial storage amount of the intermediate products in the workshop z; v is a cell z,t Representing the processing rate of the workshop z at different moments;
the above expression shows that the storage capacity of the intermediate product is kept unchanged at the beginning and the end;
in the above formula:representing the maximum storage capacity of the intermediate product in the workshop z;
υ z,t =x z,t ·υ z ,
in the above formula: v is a cell z Representing the plant z processing rate.
9. The method of claim 8, wherein the method comprises: the specific process of step S3 is as follows:
on the premise of finishing daily production tasks of a factory, taking total energy consumption cost of one day in an industrial park as a target function, specifically comprising electricity purchasing cost to a power grid, gas purchasing cost to a natural gas station, carbon dioxide emission cost, income of carbon dioxide sale capture and income of oxygen sale as a byproduct of hydrogen production by electrolyzing water;
specifically, the objective function expression of the total energy consumption cost of the industrial park in one day is as follows:
in the above formula: c total Represents the total energy consumption cost of one day of the industrial park; c t Represents the total cost of energy consumption at a certain time in one day of the industrial park;
wherein:
in the above formula:representing the power purchased by the power grid; p is a radical of E,t Representing the electricity purchase price of the power grid;representing the gas purchasing power of the natural gas station; p is a radical of g,t Representing the gas purchase price of the natural gas station;represents the carbon dioxide emission rate;represents a carbon dioxide emission price;represents a carbon dioxide purchase rate;represents a carbon dioxide purchase price;representing the rate of oxygen production from the electrolysis of water to produce hydrogen;represents the price at which oxygen is sold;represents the rate at which carbon dioxide is sold;represents the price for selling carbon dioxide;
the constraint conditions include: the method comprises the following steps of electric heating and cooling power balance constraint, power grid and gas grid interaction power constraint, operation constraint of each energy device, energy storage constraint and factory production constraint;
the relation expression among the electricity, heat and cold produced by the comprehensive energy system is as follows:
in the above formula:representing the power of electricity, heat and cold provided by the park comprehensive energy system; represents the charge and discharge power of the storage battery;representing photovoltaic productsElectrical power;representing the electricity converted into the electricity consumption power;representing the power of electricity-to-gas production;the heat production power of the solar water heater is represented;
wherein:
in the above formula: ω represents a percentage;representing the heat production efficiency and the electricity production efficiency of the cogeneration unit;representing gas boiler efficiency; COP EC Representing the conversion coefficient of the electric refrigerator; COP (coefficient of Performance) AC Representing the conversion coefficient of the absorption refrigerator;
specifically, the balance constraints of the comprehensive energy system and the electricity, heat and cold powers of the industrial users in the park are as follows:
specifically, the power grid and gas grid interaction power constraint is as follows:
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210294463.8A CN114818269A (en) | 2022-03-23 | 2022-03-23 | Industrial park comprehensive energy system and factory production plan collaborative optimization method |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210294463.8A CN114818269A (en) | 2022-03-23 | 2022-03-23 | Industrial park comprehensive energy system and factory production plan collaborative optimization method |
Publications (1)
Publication Number | Publication Date |
---|---|
CN114818269A true CN114818269A (en) | 2022-07-29 |
Family
ID=82531402
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202210294463.8A Pending CN114818269A (en) | 2022-03-23 | 2022-03-23 | Industrial park comprehensive energy system and factory production plan collaborative optimization method |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN114818269A (en) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN117911055A (en) * | 2024-03-19 | 2024-04-19 | 国网辽宁省电力有限公司技能培训中心 | Carbon emission optimizing system based on regional comprehensive energy coupling characteristics |
-
2022
- 2022-03-23 CN CN202210294463.8A patent/CN114818269A/en active Pending
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN117911055A (en) * | 2024-03-19 | 2024-04-19 | 国网辽宁省电力有限公司技能培训中心 | Carbon emission optimizing system based on regional comprehensive energy coupling characteristics |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN111738502B (en) | Multi-energy complementary system demand response operation optimization method for promoting surplus wind power consumption | |
Chen et al. | An optimization on an integrated energy system of combined heat and power, carbon capture system and power to gas by considering flexible load | |
CN111463836B (en) | Comprehensive energy system optimal scheduling method | |
CN113315242B (en) | Virtual wind abandoning-hydrogen production combination for promoting wind abandoning consumption based on hydrogen energy economy | |
CN109146182A (en) | The economic load dispatching method of meter and the distributed triple-generation system of a variety of energy storage | |
CN111639824A (en) | Thermoelectric optimization scheduling method for regional comprehensive energy system with electric-to-gas conversion function | |
CN114091913B (en) | Low-carbon economic dispatching method considering heat supply network and P2G multi-park comprehensive energy system | |
CN114580863A (en) | Regional comprehensive energy system economic dispatching method of carbon-containing capture equipment and photo-thermal power station considering comprehensive demand response | |
CN116542370A (en) | Park low-carbon economic operation method considering carbon capture and carbon transaction | |
CN112270433A (en) | Micro-grid optimization method considering renewable energy uncertainty and user satisfaction | |
CN112085263A (en) | User side distributed energy system hybrid energy storage optimal configuration method and system | |
CN117081143A (en) | Method for promoting coordination and optimization operation of park comprehensive energy system for distributed photovoltaic on-site digestion | |
CN114818269A (en) | Industrial park comprehensive energy system and factory production plan collaborative optimization method | |
Wang et al. | Low carbon optimal operation of integrated energy system based on concentrating solar power plant and power to hydrogen | |
CN113806952A (en) | Source-load-storage-considered cooling, heating and power comprehensive energy system and optimized operation method thereof | |
CN113128868A (en) | Regional comprehensive energy system scheduling optimization method and device | |
CN110516863A (en) | A kind of more microgrid active distribution system dual blank-holders of supply of cooling, heating and electrical powers type | |
CN114676897A (en) | Optimal scheduling method for comprehensive energy system of park containing CHP-P2G-hydrogen energy | |
Shi et al. | Research on energy management strategy of integrated energy system | |
Miao et al. | Low carbon and economic operation of integrated energy system considering electricity-thermal flexible load | |
Han et al. | Optimized Dispatch of Integrated Energy System with Hydrogen Energy and Carbon Capture Under Demand Side Response Mechanisms | |
Chen et al. | Flexible transformation cost analysis of Thermal power in Jiangxi Province for depth peak-shaving | |
CN114362152B (en) | Multi-time scale scheduling method for comprehensive energy system | |
Yan et al. | Optimal Scheduling of Integrated Energy System considering Hydrogen and Integrated Demand Response | |
Zhang et al. | Distributed Bargaining Method of a Multi-integrated Energy System Based on Nash Theory |
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 |