A collaborative optimization method for industrial park integrated energy system and factory production plan

Abstract

The invention discloses a method for collaborative optimization of an integrated energy system and a factory production plan in an industrial park. Aiming at the structure of the integrated energy system and industrial users in the industrial park, a hydrogen production model, a carbon dioxide recovery and utilization and storage model, and a production raw material storage area model are established on the basis of the original energy supply, conversion, and storage equipment models. A multi-energy coupling model of energy consumption in different workshops of industrial users and a factory production constraint model are constructed. Taking the operation efficiency of all equipment in the integrated energy system in the park, the multi-energy coupling model of different workshops of industrial users, and the factory production process as constraint conditions, and taking the total energy consumption cost of the industrial park in one day as the objective function, a collaborative optimization model of the operation of the integrated energy system in the park and the industrial users is constructed. The collaborative optimization model is solved to obtain the output plan of each equipment in the integrated energy system in the park and the factory production scheduling arrangement, thereby reducing the total daily operation energy consumption cost of the integrated energy system and industrial users in the park.

Description

Technical Field

[0001] The present invention relates to the technical field of energy configuration optimization, and in particular to a method for collaborative optimization of an industrial park integrated energy system and a factory production plan. Background Art

[0002] With the increasing demand for various energy sources such as electricity, natural gas, heating and cooling, modern power systems are developing towards integrated energy systems. Integrated energy systems are considered to be an effective model for reducing the use of fossil fuels, reducing energy costs, and improving the efficiency and flexibility of traditional energy systems. Its energy input includes natural gas and electricity, and its output types can include multiple energy sources such as electricity, heating, cooling, and natural gas at the same time. Since industrial loads require a large amount of electricity, heat, and cooling energy, they play an important role among end users; however, the traditional factory operation model lacks the management and optimization of multiple energy sources, which seriously affects the economic benefits and energy efficiency of enterprises; therefore, it is necessary to study the energy consumption optimization strategy of the integrated energy system of industrial parks and industrial users.

[0003] With the increasing attention paid to environmental issues, many research works have taken carbon emissions into account in integrated energy systems, and proposed various frameworks to optimize the operation of integrated energy systems in parks to reduce energy costs. In addition, in demand response, the management and optimization of industrial loads are widely used to improve economic benefits and energy efficiency. However, there are few studies on the coordinated optimization of energy conversion in integrated energy systems in parks and production plans of factory users. In addition, the recycling and utilization of carbon dioxide in integrated energy systems in parks and the coupling relationship between different energy sources consumed by factory production are not considered, and the daily production tasks of industrial users are not considered for demand response, which seriously affects the enthusiasm of users to participate. Summary of the invention

[0004] The technical problem to be solved by the present invention is to provide a method for collaborative optimization of the integrated energy system of an industrial park and the factory production plan, thereby reducing the total daily operating energy consumption cost of the integrated energy system and industrial users in the park, improving the comprehensive utilization rate of energy, and reducing carbon emissions.

[0005] In order to solve the above technical problems, the technical solution adopted by the present invention is: a method for collaborative optimization of an industrial park integrated energy system and a factory production plan, comprising the following steps:

[0006] Step S1: Based on the existing energy supply, conversion and storage equipment models, a hydrogen production model, a carbon dioxide recycling and storage model, and a production raw material storage area model are established for the integrated energy system and the structure of industrial users in the industrial park;

[0007] Step S2: The park integrated energy management center receives the production tasks, workshop energy consumption types, and production constraints of industrial users, and builds a multi-energy coupling model for energy consumption of different workshops of industrial users and a factory production constraint model;

[0008] Step S3: The park integrated energy management center integrates the received industrial user information, equipment information of the integrated energy system, price and power information of the external gas network and power grid, and builds a collaborative optimization model of the integrated energy system and industrial user operation on the premise of completing the daily production task of the factory, with the operating efficiency of all equipment in the integrated energy system in the park, the multi-energy coupling model of different workshops of industrial users and the factory production process as constraints, and the total energy consumption cost of the industrial park in one day as the objective function;

[0009] Step S4: The park integrated energy management center solves the collaborative optimization model of the integrated energy system and the operation of industrial users to obtain the output plan of each device in the park integrated energy system and the production scheduling arrangement of the factory;

[0010] Step S5: The park integrated energy management center controls the operation of the integrated energy system according to the solution results, and sends the obtained results to the factory users so that they can carry out production activities according to the corresponding arrangements.

[0011] A further improvement of the technical solution of the present invention is that: the energy supply equipment model in step S1 includes a natural gas grid, a power grid, a photovoltaic power generation unit, and a solar water heater; the energy conversion equipment model includes a gas boiler, a gas turbine, a P2G unit, an absorption refrigerator, and an electric refrigerator; the energy storage equipment model includes an electric energy storage device; the integrated energy system hydrogen production model includes a water electrolysis unit; the carbon dioxide recycling and storage model includes a carbon dioxide capture device; the production raw material storage area model includes 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.

[0012] A further improvement of the technical solution of the present invention is that the P2G unit and the water electrolysis unit model are:

[0013]

[0014] In the above formula: Indicates the rate at which the power-to-gas unit produces natural gas; δ indicates the proportionality coefficient between the rate at which hydrogen is consumed and the rate at which natural gas is produced; represents the rate at which hydrogen is produced by the water electrolyzer;

[0015]

[0016] In the above formula: represents the rate at which the power-to-gas unit consumes carbon dioxide; λ represents the proportional coefficient between the rate at which hydrogen is consumed and the rate at which carbon dioxide is consumed in the production of natural gas;

[0017]

[0018] In the above formula: Indicates the power consumed by the power-to-gas unit; Indicates the conversion coefficient between the rate of hydrogen consumption and power consumption;

[0019]

[0020] In the above formula: Indicates the maximum power consumed by the power-to-gas unit;

[0021]

[0022] In the above formula: Indicates the rate at which the electrolytic water tank produces oxygen; τ indicates the proportional coefficient between the rate at which hydrogen is produced and the rate at which oxygen is produced;

[0023]

[0024] In the above formula: represents the rate at which the power-to-gas unit produces natural gas; θ represents the coefficient of converting the natural gas production rate into natural gas power.

[0025] A further improvement of the technical solution of the present invention is that the carbon dioxide recovery, utilization and storage model of the integrated energy system in step S1 is:

[0026]

[0027] In the above formula: represents the rate of capturing carbon dioxide produced by cogeneration units and gas boilers; represents the rate of carbon dioxide being emitted into the air due to the limitation of capture technology; ρ represents the conversion coefficient; Indicates the heating power of the gas boiler; Indicates the heat generation efficiency of gas boiler; They represent the electricity and heat production of the cogeneration unit respectively; They represent the electricity and heat production efficiencies of the cogeneration unit respectively;

[0028]

[0029] In the above formula: represents the rate at which recovered CO2 is consumed by the P2G unit; Indicates the rate at which CO2 is stored in the CO2 tank.

[0030] A further improvement of the technical solution of the present invention is that the carbon dioxide storage model of the integrated energy system in step S1 is:

[0031]

[0032] In the above formula: CO 2s,t represents the amount of carbon dioxide stored; represents the initial storage amount of carbon dioxide; represents the CO2 sales rate;

[0033]

[0034] The above formula indicates that the amount of carbon dioxide stored at the beginning and end of the carbon dioxide period remains unchanged;

[0035] 0≤CO 2s,t ≤CO 2s,max ,

[0036] In the above formula: CO 2s,max Indicates the maximum storage capacity of carbon dioxide;

[0037]

[0038] In the above formula: Indicates the maximum storage rate of carbon dioxide.

[0039] A further improvement of the technical solution of the present invention is that the workshop energy consumption types in step S2 include pure electric energy workshop, pure cold energy workshop, pure heat energy workshop, electric-heat coupling workshop, electric-cold coupling workshop, and heat-cold coupling workshop.

[0040] A further improvement of the technical solution of the present invention is that the multi-energy coupling model of energy consumption of different workshops of industrial users in step S2 is:

[0041]

[0042] In the above formula: P E,t represents the total power consumed by the factory; z represents the workshop number; E' represents the set of workshops that consume pure power; x is binary, representing the operation status of the workshop, 1 represents the workshop working, otherwise it does not work; P E,z represents the electric power consumed by workshop z; EH' represents the set of hotspot coupled workshops; represents the electric power consumed by the electric-heat coupling workshop z; EC' represents the set of electric-cooling coupling workshops; It represents the electric power consumed by the electric-cooling coupling workshop z;

[0043]

[0044] In the above formula: P H,trepresents the total power of thermal energy consumed by the factory; H' represents the set of workshops consuming pure thermal energy; P H,z represents the thermal power consumed by workshop z; represents the total thermal power consumed by the electric-thermal coupling workshop; HC' represents the set of thermal-cold coupling workshops; represents the thermal power consumed by the heat-cold coupling workshop z;

[0045]

[0046] In the above formula: P C,t represents the total power of cooling energy consumed by the factory; C' represents the set of workshops consuming pure cooling energy; It indicates the total cooling power consumed by the electric-cooling coupled workshop; It indicates the total cooling power consumed by the hot-cold coupled workshop;

[0047] in:

[0048]

[0049]

[0050]

[0051] In the above formula: represents the electrothermal coupling coefficient of workshop z; It represents the electric-cooling coupling coefficient of workshop z; Represents the thermal-cold coupling coefficient of workshop z.

[0052] A further improvement of the technical solution of the present invention is that: the production constraints in step S2 include the production sequence, production speed, operation status of different workshops, and the initial quantity and remaining capacity of the storage warehouse;

[0053] P (EHC),z Indicates the various coupled energy consumption power of workshop z; υ={υ 1 … Z} represents the production rate of different workshops; represents the initial storage quantity of the intermediate products in workshop z; represents the storage capacity of intermediate products in workshop z;

[0054] Specifically, the intermediate product storage warehouse model:

[0055]

[0056] In the above formula: represents the storage quantity of intermediate products in workshop z; represents the initial storage volume of intermediate products in workshop z; z,t It represents the processing rate of workshop z at different times;

[0057]

[0058] The above formula indicates that the storage volume of intermediate products remains unchanged at the initial and final moments;

[0059]

[0060] In the above formula: represents the maximum storage capacity of the intermediate products of workshop z;

[0061] υ z,t =x z,t ·υ z ,

[0062] In the above formula: z Represents the processing rate of workshop z.

[0063] A further improvement of the technical solution of the present invention is that the specific process of step S3 is as follows:

[0064] On the premise of completing the daily production task of the factory, the total energy consumption cost of the industrial park in one day is taken as the objective function, which specifically includes the cost of purchasing electricity from the power grid, the cost of purchasing gas from the natural gas station, the cost of carbon dioxide emissions, the income from the sale of captured carbon dioxide, and the income from the sale of oxygen, a byproduct of hydrogen production from electrolysis of water;

[0065] Specifically, the objective function expression of the total energy consumption cost of the industrial park in one day is:

[0066]

[0067] In the above formula: C total represents the total energy consumption cost of the industrial park in one day; C t It represents the total energy consumption cost of the industrial park at a certain time of the day;

[0068] in:

[0069]

[0070] In the above formula: Indicates the power purchased by the power grid; p E,t Indicates the price of electricity purchased from the power grid; Indicates the gas purchasing power of the natural gas station; p g,t Indicates the gas purchase price at the natural gas station; represents the carbon dioxide emission rate; represents the price of carbon dioxide emissions; represents the CO2 purchase rate; represents the purchase price of CO2; It indicates the rate at which oxygen is produced by electrolyzing water to produce hydrogen; represents the price at which oxygen is sold; represents the rate at which CO2 is sold; represents the price at which CO2 is sold;

[0071] Constraints include: power balance constraints for electricity, heating and cooling, interactive power constraints for the power grid and gas grid, operation constraints for various energy equipment, energy storage constraints, and factory production constraints;

[0072] The relationship between electricity, heat and cold produced by the integrated energy system is expressed as:

[0073]

[0074] In the above formula: Indicates the electricity, heat and cooling power provided by the park's integrated energy system; Indicates the battery charging and discharging power; Indicates photovoltaic power generation; Indicates the power consumption of power-to-gas conversion; Indicates the power of electricity-to-gas conversion; Indicates the heating power of the solar water heater;

[0075] in:

[0076]

[0077] In the above formula: ω represents a percentage; It represents the heat production efficiency and power production efficiency of the cogeneration unit; Indicates gas boiler efficiency; COP EC Indicates the conversion coefficient of electric refrigerator; COP AC It represents the conversion factor of absorption chiller;

[0078] Specifically, the power balance constraints of the integrated energy system in the park and industrial users in terms of electricity, heat and cooling are:

[0079]

[0080]

[0081]

[0082] Specifically, the interactive power constraints of the power grid and gas grid are:

[0083]

[0084]

[0085] In the above formula: Indicates the maximum power of electricity purchased from the power grid and gas purchased from the gas grid.

[0086] Due to the adoption of the above technical solution, the technical progress achieved by the present invention is:

[0087] The present invention proposes a method for collaborative optimization of an integrated energy system of an industrial park and a factory production plan. Aiming at the structure of the integrated energy system and industrial users in the industrial park, a hydrogen production model, a carbon dioxide recovery and utilization and storage model, and a production raw material storage area model are established on the basis of the original energy supply, conversion, and storage equipment models; a multi-energy coupling model of energy consumption in different workshops of industrial users and a factory production constraint model are constructed; on the premise of completing the daily production task of the factory, the operating efficiency of all equipment in the integrated energy system in the park, the multi-energy coupling model of different workshops of industrial users, and the factory production process are used as constraints, and the total energy consumption cost of the industrial park in one day is used as the objective function to construct a collaborative optimization model of the operation of the integrated energy system in the park and industrial users; the constructed collaborative optimization model is solved to obtain the output plan of each equipment in the integrated energy system of the park and the production scheduling arrangement of the factory, which can reduce the total daily operating energy consumption cost of the integrated energy system and industrial users in the park, improve the comprehensive energy utilization rate, and reduce carbon emissions. BRIEF DESCRIPTION OF THE DRAWINGS

[0088] Figure 1 It is a block diagram of the steps of collaborative optimization of the integrated energy system and industrial user operation in the park of the present invention;

[0089] Figure 2 It is a system block diagram of the collaborative optimization model of the integrated energy system and industrial users in the park of the present invention;

[0090] Figure 3 It is a production schematic diagram of different workshops in the factory park of the present invention. DETAILED DESCRIPTION

[0091] The present invention is further described in detail below in conjunction with embodiments:

[0092] like Figure 1 As shown, a method for collaborative optimization of an industrial park integrated energy system and a factory production plan includes the following steps:

[0093] Step S1: Based on the existing energy supply, conversion and storage equipment models, a comprehensive energy system hydrogen production model, a carbon dioxide recycling and storage model, and a production raw material storage area model are established for the structure of the comprehensive energy system and industrial users in the industrial park. Figure 2As shown in the figure, it describes the supply and demand relationship between the production of electricity, heat and cold energy by the park's comprehensive energy system and the consumption of electricity, heat and cold energy by different workshops in the factory, and realizes the optimization of the comprehensive energy system and the production scheduling of the factory through the coordinated optimization between the two, thereby reducing the total daily production energy consumption cost of the park's comprehensive energy system and the factory; among them, the energy supply equipment model of the comprehensive energy system includes natural gas grid, power grid, photovoltaic power generation unit, solar water heater; the energy conversion equipment model includes gas boiler, gas turbine, P2G unit, absorption refrigerator, electric refrigerator; the energy storage equipment model includes electric energy storage equipment; the hydrogen production model of the comprehensive energy system includes water electrolysis unit; the carbon dioxide recycling and storage model includes carbon dioxide capture equipment; there are N production workshops in the factory, and the production raw material storage area model includes initial raw material warehouse, workshop 1 intermediate product storage warehouse, workshop 2 intermediate product storage warehouse, workshop z intermediate product storage warehouse.

[0094] The P2G unit and water electrolysis unit models are:

[0095]

[0096] In the above formula: Indicates the rate at which the power-to-gas unit produces natural gas; δ indicates the proportionality coefficient between the rate at which hydrogen is consumed and the rate at which natural gas is produced; represents the rate at which hydrogen is produced by the water electrolyzer;

[0097]

[0098] In the above formula: represents the rate at which the power-to-gas unit consumes carbon dioxide; λ represents the proportional coefficient between the rate at which hydrogen is consumed and the rate at which carbon dioxide is consumed in the production of natural gas;

[0099]

[0100] In the above formula: Indicates the power consumed by the power-to-gas unit; Indicates the conversion coefficient between the rate of hydrogen consumption and power consumption;

[0101]

[0102] In the above formula: Indicates the maximum power consumed by the power-to-gas unit;

[0103]

[0104] In the above formula: Indicates the rate at which the electrolytic water tank produces oxygen; τ indicates the proportional coefficient between the rate at which hydrogen is produced and the rate at which oxygen is produced;

[0105]

[0106] In the above formula: represents the rate at which the power-to-gas unit produces natural gas; θ represents the coefficient of converting the natural gas production rate into natural gas power.

[0107] The carbon dioxide recovery and utilization model of the integrated energy system is:

[0108]

[0109] In the above formula: represents the rate of capturing carbon dioxide produced by cogeneration units and gas boilers; represents the rate of carbon dioxide being emitted into the air due to the limitation of capture technology; ρ represents the conversion coefficient; Indicates the heating power of the gas boiler; Indicates the heat generation efficiency of gas boiler; They represent the electricity and heat production of the cogeneration unit respectively; They represent the electricity and heat production efficiencies of the cogeneration unit respectively;

[0110]

[0111] In the above formula: represents the rate at which recovered CO2 is consumed by the P2G unit; Indicates the rate at which CO2 is stored in the CO2 tank.

[0112] The carbon dioxide storage model of the integrated energy system is:

[0113] In the above formula: CO 2s,t represents the amount of carbon dioxide stored; represents the initial storage amount of carbon dioxide; represents the CO2 sales rate;

[0114]

[0115] The above formula indicates that the amount of carbon dioxide stored at the beginning and end of the carbon dioxide period remains unchanged;

[0116] 0≤CO 2s,t ≤CO 2s,max ,

[0117] In the above formula: CO 2s,max Indicates the maximum storage capacity of carbon dioxide;

[0118]

[0119] In the above formula: Indicates the maximum storage rate of carbon dioxide.

[0120] Step S2: The park integrated energy management center receives the production tasks, workshop energy consumption types, and production constraints of industrial users, and constructs a multi-energy coupling model of the energy consumption of different workshops of industrial users and a factory production constraint model, wherein the workshop energy consumption types include pure electric energy workshops, pure cold energy workshops, pure heat energy workshops, electric-heat coupling workshops, electric-cold coupling workshops, and heat-cold coupling workshops. The multi-energy coupling model of the energy consumption of different workshops of industrial users is:

[0121]

[0122] In the above formula: P E,t represents the total power consumed by the factory; z represents the workshop number; E' represents the set of workshops that consume pure power; x is binary, representing the operation status of the workshop, 1 represents the workshop working, otherwise it does not work; P E,z represents the electric power consumed by workshop z; EH' represents the set of hotspot coupled workshops; represents the electric power consumed by the electric-heat coupling workshop z; EC' represents the set of electric-cooling coupling workshops; It represents the electric power consumed by the electric-cooling coupling workshop z;

[0123]

[0124] In the above formula: P H,t represents the total power of thermal energy consumed by the factory; H' represents the set of workshops consuming pure thermal energy; P H,z represents the thermal power consumed by workshop z; represents the total thermal power consumed by the electric-thermal coupling workshop; HC' represents the set of thermal-cold coupling workshops; represents the thermal power consumed by the heat-cold coupling workshop z;

[0125]

[0126] In the above formula: P C,t represents the total power of cooling energy consumed by the factory; C' represents the set of workshops consuming pure cooling energy; It indicates the total cooling power consumed by the electric-cooling coupled workshop; It indicates the total cooling power consumed by the hot-cold coupled workshop;

[0127] in:

[0128]

[0129]

[0130]

[0131] In the above formula: represents the electrothermal coupling coefficient of workshop z; It represents the electric-cooling coupling coefficient of workshop z; Represents the thermal-cold coupling coefficient of workshop z.

[0132] Production constraints include the production sequence, production speed, and operation status of different workshops, as well as the initial quantity and remaining capacity of the storage warehouse;

[0133] like Figure 3 As shown, P (EHC),z Indicates the various coupled energy consumption power of workshop z; υ={υ 1 … Z} represents the production rate of different workshops; represents the initial storage quantity of the intermediate products in workshop z; represents the storage capacity of intermediate products in workshop z;

[0134] Specifically, the intermediate product storage warehouse model:

[0135]

[0136] In the above formula: represents the storage quantity of intermediate products in workshop z; represents the initial storage volume of intermediate products in workshop z; z,t It represents the processing rate of workshop z at different times;

[0137]

[0138] The above formula indicates that the storage volume of intermediate products remains unchanged at the initial and final moments;

[0139]

[0140] In the above formula: represents the maximum storage capacity of the intermediate products of workshop z;

[0141] υ z,t =x z,t ·υ z ,

[0142] In the above formula: z Represents the processing rate of workshop z.

[0143] Step S3: The park integrated energy management center will receive industrial user information, equipment information of the integrated energy system, price and power information of the external gas network and power grid, and integrate them. On the premise of completing the daily production task of the factory, the operation efficiency of all equipment in the integrated energy system in the park, the multi-energy coupling model of different workshops of industrial users and the factory production process are used as constraints, and the total energy consumption cost of the industrial park in one day is used as the objective function, and a collaborative optimization model of the operation of the integrated energy system and industrial users is constructed; the specific process is as follows:

[0144] On the premise of completing the daily production task of the factory, the total energy consumption cost of the industrial park in one day is taken as the objective function, which specifically includes the cost of purchasing electricity from the power grid, the cost of purchasing gas from the natural gas station, the cost of carbon dioxide emissions, the income from the sale of captured carbon dioxide, and the income from the sale of oxygen, a byproduct of hydrogen production from electrolysis of water;

[0145] Specifically, the objective function expression of the total energy consumption cost of the industrial park in one day is:

[0146]

[0147] In the above formula: C total represents the total energy consumption cost of the industrial park in one day; C t It represents the total energy consumption cost of the industrial park at a certain time of the day;

[0148] in:

[0149]

[0150] In the above formula: Indicates the power purchased by the power grid; p E,t Indicates the price of electricity purchased from the power grid; Indicates the gas purchasing power of the natural gas station; p g,t Indicates the gas purchase price at the natural gas station; represents the carbon dioxide emission rate; represents the price of carbon dioxide emissions; represents the CO2 purchase rate; represents the purchase price of CO2; It indicates the rate at which oxygen is produced by electrolyzing water to produce hydrogen; represents the price at which oxygen is sold; represents the rate at which CO2 is sold; represents the price at which CO2 is sold;

[0151] Constraints include: power balance constraints for electricity, heating and cooling, interactive power constraints for the power grid and gas grid, operation constraints for various energy equipment, energy storage constraints, and factory production constraints;

[0152] The relationship between electricity, heat and cold produced by the integrated energy system is expressed as:

[0153]

[0154] In the above formula: Indicates the electricity, heat and cooling power provided by the park's integrated energy system; Indicates the battery charging and discharging power; Indicates photovoltaic power generation; Indicates the power consumption of power-to-gas conversion; Indicates the power of electricity-to-gas conversion; Indicates the heating power of the solar water heater;

[0155] in:

[0156]

[0157] In the above formula: ω represents a percentage; It represents the heat production efficiency and power production efficiency of the cogeneration unit; Indicates gas boiler efficiency; COP EC Indicates the conversion coefficient of electric refrigerator; COP AC It represents the conversion factor of absorption chiller;

[0158] Specifically, the power balance constraints of the integrated energy system in the park and industrial users in terms of electricity, heat and cooling are:

[0159]

[0160]

[0161]

[0162] Specifically, the interactive power constraints of the power grid and gas grid are:

[0163]

[0164]

[0165] In the above formula: Indicates the maximum power of electricity purchased from the power grid and gas purchased from the gas grid.

[0166] Step S4: The park integrated energy management center solves the collaborative optimization model of the integrated energy system and the operation of industrial users to obtain the output plan of each device in the park integrated energy system and the production scheduling arrangement of the factory;

[0167] Step S5: The park integrated energy management center controls the operation of the integrated energy system according to the solution results, and sends the obtained results to the factory users so that they can carry out production activities according to the corresponding arrangements.

[0168] The present invention proposes a method for collaborative optimization of an integrated energy system of an industrial park and a factory production plan. Aiming at the structure of the integrated energy system and industrial users in the industrial park, a hydrogen production model, a carbon dioxide recovery and utilization and storage model, and a production raw material storage area model are established on the basis of the original energy supply, conversion, and storage equipment models; a multi-energy coupling model of energy consumption in different workshops of industrial users and a factory production constraint model are constructed; on the premise of completing the daily production task of the factory, the operating efficiency of all equipment in the integrated energy system in the park, the multi-energy coupling model of different workshops of industrial users, and the factory production process are used as constraints, and the total energy consumption cost of the industrial park in one day is used as the objective function to construct a collaborative optimization model of the operation of the integrated energy system in the park and industrial users; the constructed collaborative optimization model is solved to obtain the output plan of each equipment in the integrated energy system of the park and the production scheduling arrangement of the factory, which can reduce the total daily operating energy consumption cost of the integrated energy system and industrial users in the park, improve the comprehensive energy utilization rate, and reduce carbon emissions.

Claims

1. A method for collaborative optimization of an industrial park integrated energy system and a factory production plan, characterized in that: The steps include: Step S1: Based on the existing energy supply, conversion and storage equipment models, a hydrogen production model, a carbon dioxide recycling and storage model and a production raw material storage area model are established for the integrated energy system and the structure of industrial users in the industrial park; Energy supply equipment models include natural gas grid, power grid, photovoltaic generators, and solar water heaters; The energy conversion equipment model includes gas boilers, gas turbines, P2G units, absorption refrigerators, and electric refrigerators; the energy storage equipment model includes electric energy storage equipment; the integrated energy system hydrogen production model includes a water electrolysis unit; the carbon dioxide recycling and storage model includes a carbon dioxide capture device; the production raw material storage area model includes the initial raw material warehouse, the workshop 1 intermediate product storage warehouse, the workshop 2 intermediate product storage warehouse, and the workshop z intermediate product storage warehouse; The P2G unit and water electrolysis unit models are: In the above formula: Indicates the rate at which the power-to-gas unit produces natural gas; δ indicates the proportionality coefficient between the rate at which hydrogen is consumed and the rate at which natural gas is produced; represents the rate at which hydrogen is produced by the water electrolyzer; In the above formula: represents the rate at which the power-to-gas unit consumes carbon dioxide; λ represents the proportional coefficient between the rate at which hydrogen is consumed and the rate at which carbon dioxide is consumed in the production of natural gas; In the above formula: Indicates the power consumed by the power-to-gas unit; Indicates the conversion coefficient between the rate of hydrogen consumption and power consumption; In the above formula: Indicates the maximum power consumed by the power-to-gas unit; In the above formula: Indicates the rate at which the electrolytic water tank produces oxygen; τ indicates the proportional coefficient between the rate at which hydrogen is produced and the rate at which oxygen is produced; In the above formula: represents the rate at which the power-to-gas unit produces natural gas; θ represents the coefficient of converting the natural gas production rate into natural gas power; The carbon dioxide recovery and utilization model of the integrated energy system is: In the above formula: represents the rate of capturing carbon dioxide produced by cogeneration units and gas boilers; represents the rate of carbon dioxide being emitted into the air due to the limitation of capture technology; ρ represents the conversion coefficient; Indicates the heating power of the gas boiler; Indicates the heat generation efficiency of gas boiler; They represent the electricity and heat production of the cogeneration unit respectively; They represent the electricity and heat production efficiencies of the cogeneration unit respectively; In the above formula: represents the rate at which recovered CO2 is consumed by the P2G unit; Indicates the rate of storage into the CO2 tank; The carbon dioxide storage model of the integrated energy system is: In the above formula: CO 2s,t represents the amount of carbon dioxide stored; represents the initial storage amount of carbon dioxide; represents the CO2 sales rate; The above formula indicates that the amount of carbon dioxide stored at the beginning and end of the carbon dioxide period remains unchanged; 0≤CO 2s,t ≤CO 2s,max , In the above formula: CO 2s,max Indicates the maximum storage capacity of carbon dioxide; In the above formula: Indicates the maximum storage rate of carbon dioxide; Step S2: The park comprehensive energy management center receives the production tasks, workshop energy consumption types, and production constraints of industrial users, and constructs a multi-energy coupling model of energy consumption of different workshops of industrial users and a factory production constraint model; the workshop energy consumption types include pure electric energy workshop, pure cold energy workshop, pure heat energy workshop, electric-heat coupling workshop, electric-cold coupling workshop, and heat-cold coupling workshop; The multi-energy coupling model of energy consumption in different workshops of industrial users is: In the above formula: P E,t represents the total power consumed by the factory; z represents the workshop number; E' represents the set of workshops that consume pure power; x is binary, representing the operation status of the workshop, 1 represents the workshop working, otherwise it does not work; P E,z represents the electric power consumed by workshop z; EH' represents the set of hotspot coupled workshops; represents the electric power consumed by the electrothermal coupling workshop z; EC' represents the electrical-cooling coupled workshop set; It represents the electric power consumed by the electric-cooling coupling workshop z; In the above formula: P H,t represents the total power of thermal energy consumed by the factory; H' represents the set of workshops consuming pure thermal energy; P H,z represents the thermal power consumed by workshop z; represents the total thermal power consumed by the electric-thermal coupling workshop; HC' represents the set of thermal-cold coupling workshops; represents the thermal power consumed by the heat-cold coupling workshop z; In the above formula: P C,t represents the total power of cooling energy consumed by the factory; C' represents the set of workshops consuming pure cooling energy; It indicates the total cooling power consumed by the electric-cooling coupled workshop; It represents the total cooling power consumed by the hot-cold coupled workshop; in: In the above formula: represents the electrothermal coupling coefficient of workshop z; It represents the electric-cooling coupling coefficient of workshop z; It represents the thermal-cold coupling coefficient of workshop z; Production constraints include the production sequence, production speed, and operation status of different workshops, as well as the initial quantity and remaining capacity of the storage warehouse; P (EHC),z Indicates the various coupled energy consumption power of workshop z; υ={υ1…υ Z } represents the production rate of different workshops; represents the initial storage quantity of the intermediate products in workshop z; represents the storage capacity of intermediate products in workshop z; Specifically, the intermediate product storage warehouse model: In the above formula: represents the storage quantity of intermediate products in workshop z; represents the initial storage volume of intermediate products in workshop z; z,t It represents the processing rate of workshop z at different times; The above formula indicates that the storage volume of intermediate products remains unchanged at the initial and final moments; In the above formula: represents the maximum storage capacity of the intermediate products of workshop z; u z,t =x z,t ·u z , In the above formula: z Indicates the processing rate of workshop z; Step S3, the park integrated energy management center will receive industrial user information, equipment information of the integrated energy system, price and power information of the external gas network and power grid, and integrate them. On the premise of completing the daily production task of the factory, the operating efficiency of all equipment in the integrated energy system in the park, the multi-energy coupling model of different workshops of industrial users and the factory production process are used as constraints, and the total energy consumption cost of the industrial park in one day is used as the objective function, and a collaborative optimization model of the integrated energy system and industrial user operation is constructed; Step S4: The park integrated energy management center solves the collaborative optimization model of the integrated energy system and the operation of industrial users to obtain the output plan of each device in the park integrated energy system and the production scheduling arrangement of the factory; Step S5: The park integrated energy management center controls the operation of the integrated energy system according to the solution results, and sends the obtained results to the factory users so that they can carry out production activities according to the corresponding arrangements.

2. The method for collaborative optimization of an industrial park integrated energy system and factory production plan according to claim 1, characterized in that: The specific process of step S3 is as follows: On the premise of completing the daily production task of the factory, the total energy consumption cost of the industrial park in one day is taken as the objective function, which specifically includes the cost of purchasing electricity from the power grid, the cost of purchasing gas from the natural gas station, the cost of carbon dioxide emissions, the income from the sale of captured carbon dioxide, and the income from the sale of oxygen, a byproduct of hydrogen production from electrolysis of water; Specifically, the objective function expression of the total energy consumption cost of the industrial park in one day is: In the above formula: C total represents the total energy consumption cost of the industrial park in one day; C t It represents the total energy consumption cost of the industrial park at a certain time of the day; in: In the above formula: Indicates the power purchased by the power grid; p E,t Indicates the price of electricity purchased from the power grid; Indicates the gas purchasing power of the natural gas station; p g,t Indicates the gas purchase price at the natural gas station; represents the carbon dioxide emission rate; represents the price of carbon dioxide emissions; represents the CO2 purchase rate; represents the purchase price of CO2; Indicates the rate at which oxygen is produced by electrolyzing water to produce hydrogen; represents the price at which oxygen is sold; represents the rate at which CO2 is sold; represents the price at which CO2 is sold; Constraints include: power balance constraints for electricity, heating and cooling, interactive power constraints for the power grid and gas grid, operation constraints for various energy equipment, energy storage constraints, and factory production constraints; The relationship between electricity, heat and cold produced by the integrated energy system is expressed as: In the above formula: Indicates the electricity, heat and cooling power provided by the park's integrated energy system; Indicates the battery charging and discharging power; Indicates photovoltaic power generation; Indicates the power consumption of electricity-to-gas conversion; Indicates the power of electricity-to-gas conversion; Indicates the heating power of the solar water heater; in: In the above formula: ω represents a percentage; It represents the heat production efficiency and power production efficiency of the cogeneration unit; Indicates gas boiler efficiency; COP EC Indicates the conversion coefficient of electric refrigerator; COP AC It represents the conversion factor of absorption chiller; Specifically, the power balance constraints of the integrated energy system in the park and industrial users in terms of electricity, heat and cooling are: Specifically, the interactive power constraints of the power grid and gas grid are: In the above formula: Indicates the maximum power of electricity purchased from the power grid and gas purchased from the gas grid.

Images