Scheduling control method and system for industrial integrated energy system, medium and processor

Abstract

The invention provides a scheduling control method for an industrial integrated energy system. The scheduling control method comprises the steps that S1, the electricity purchase price and the electricity purchase quantity of the integrated energy system are acquired; s2, acquiring a gas purchase price and a gas purchase amount of the integrated energy system; s3, acquiring a hydrogen purchase price and a hydrogen purchase amount of the comprehensive energy system; s4, the wind curtailment punishment price and the wind curtailment power of the integrated energy system are obtained; s5, calculating the carbon transaction cost of the integrated energy system; s6, calculating a translational power variable quantity and a replaceable demand response power variable quantity of the integrated energy system; and S7, establishing a scheduling control model of the integrated energy system, and calculating the total operation cost minimization of the system according to the data. By comprehensively and systematically analyzing the carbon emission process of the comprehensive energy industry park, deeply discussing the source and the flow direction of carbon emission and accurately estimating the capacity of IES by adopting a comprehensive carbon emission measurement method, the energy utilization efficiency is improved to the maximum extent, and meanwhile, the carbon emission is reduced to the maximum extent.

Description

Technical Field

[0001] The present invention relates to the technical field of carbon emission of energy systems, and in particular to a dispatching control method, system, medium and processor for an industrial integrated energy system. Background Art

[0002] In the context of global energy crisis and climate change, integrated energy system (IES) has attracted widespread attention as an important way to improve energy efficiency and reduce carbon emissions. Integrated energy industrial park integrates various forms of energy and industry, and its carbon emissions are directly related to the sustainable development of the region's energy structure. Comprehensive analysis of the formation mechanism of carbon flow within the park will help to accurately identify the main sources and distribution of carbon emissions, thereby providing a basis for the implementation of scientific and reasonable carbon reduction strategies. In addition, by optimizing the energy system within the park and improving energy efficiency, a solid foundation can be laid for achieving carbon reduction goals. Yu et al. proposed a robust emission inventory method based on enterprise-level data of Chinese industrial parks. Figure 1 A worldwide literature survey on carbon emissions from industrial parks is presented.

[0003] Europe and the United States have systematically revealed the characteristics and patterns of carbon emissions in comprehensive environmental strategies through methods such as life cycle assessment (LCA). Energy flow analysis and carbon accounting life cycle assessment can consider the entire life cycle and conduct a comprehensive and systematic assessment of carbon emissions. Srinath et al. proposed a method based on a calcium cycle carbon capture system to evaluate the life cycle energy, carbon and water footprint of a 660MWe supercritical coal-fired power plant in India. Wang et al. introduced a life cycle model to examine the impact of the circular economy on carbon emissions in the coal-fired power generation sector. Carbon energy flow has the advantage of quickly identifying the main sources and key stages of carbon emissions. Wang et al. proposed a low-carbon optimization scheduling model that aims to control carbon intensity by implementing demand response, thereby solving the problem of increased anthropogenic carbon emissions in the power system. Zhai et al. proposed a provincial energy-carbon interaction model to solve the problem of revealing provincial energy-carbon relationships. Liu et al. proposed a new method for real-time estimation of direct and indirect carbon emissions of a plant, which solves the limitations of traditional estimation methods in terms of comprehensiveness, timeliness and accuracy. The above studies have achieved gratifying practical results in managing energy carbon flows. However, there are few energy scheduling calculations and control methods involving carbon emissions in comprehensive energy industrial parks based on industry and integrating electricity, hydrogen, heat and gas.

[0004] In view of this, a dispatching control method, system, medium and processor for an industrial integrated energy system are needed. Summary of the invention

[0005] In view of the problem that there is no energy scheduling method involving carbon emissions in the integrated energy industrial park integrating electricity, hydrogen, heat and gas in the prior art, the present invention provides a scheduling control method, system, medium and processor for an industrial integrated energy system, which can use a comprehensive carbon emission measurement method to accurately estimate the capacity of IES, which helps to make wise decisions on the carbon emission reduction strategy of the industrial park. At the same time, it proposes an optimized energy scheduling strategy to minimize carbon emissions while maximizing energy utilization efficiency. The specific technical solution is as follows:

[0006] A dispatching control method for an industrial integrated energy system, comprising:

[0007] S1: Obtain the electricity purchase price and electricity purchase amount of the integrated energy system;

[0008] S2: Obtain the gas purchase price and gas purchase volume of the integrated energy system;

[0009] S3: Obtain the purchase price and purchase amount of hydrogen for the integrated energy system;

[0010] S4: Obtain the wind curtailment penalty price and wind curtailment power of the integrated energy system;

[0011] S5: Calculate the carbon trading costs of the integrated energy system;

[0012] S6: Calculate the shiftable power change of the integrated energy system and the replaceable demand response power change;

[0013] S7: Establish a dispatching and control model for the integrated energy system, and minimize the total system operation cost based on the above data.

[0014] Furthermore, the dispatching control model for establishing the integrated energy system is as follows:

[0015]

[0016] In the above formula, P e，buy (t), P g，buy (t) They represent the electricity, gas and hydrogen purchases at time t respectively; c e (t) is the electricity purchase price at time t; c g (t) is the gas purchase price at time t; is the price of hydrogen purchased at time t; cur is the penalty price for wind curtailment at time t; is the abandoned wind power; For carbon trading costs; are the shiftable and substitutable demand response prices at time t, respectively; is the amount of power variation that can be translated; is the amount of alternative demand response power change.

[0017] Furthermore, the calculation of the carbon trading cost of the integrated energy system includes the following steps:

[0018] S51: Calculate the carbon quota of the integrated energy system using the following formula:

[0019] E IES =E e，buy +E GT +E GB +E g，Load ;in,

[0020]

[0021] In the above formula, E IES 、E e，buy 、E GT 、E GB 、E g，Load They represent the carbon quota of the integrated energy system and the carbon quota caused by electricity purchased from the main grid, gas turbines, gas boilers and gas loads; μ e , μ h , μ g，Load represents the carbon quota coefficient calculated per unit of electricity, per unit of thermal power and per unit of gas load; P e，buy (t), P e，GT (t), P h，GT (t), P h，GB (t), P g，Load (t) represents the electricity purchase amount, the electric power output of the gas turbine, the thermal power output of the gas turbine, the thermal power output of the gas boiler and the gas load demand at time t respectively;

[0022] S52: Calculate the actual carbon emissions of the integrated energy system using the following formula:

[0023] E IES，a =E e，buy，a +E GT，a +E GB，a +E g，Load，a -E CCS，a ;

[0024]

[0025] In the above formula, E IES，a 、E e，buy，a 、E GT，a 、E GB，a 、E g，Load，a They represent the actual carbon emissions of the integrated energy system and the actual carbon emissions caused by purchasing electricity from the main grid, gas turbines, gas boilers and gas loads; E CCS，a is the amount of carbon actually captured; They are the first, second and third calculation parameters of the actual carbon emissions of coal-fired power units; are the actual carbon quota coefficients calculated based on unit electricity, unit heat energy, and unit gas load consumption; E CCS (t) is the total amount of CO2 captured during CCS operation at time t;

[0026] S53: The actual carbon difference of the comprehensive energy system is:

[0027] E IES，t =E IES，a -E IES ;

[0028] In the above formula, E IES，t Represents the carbon emission gap between actual carbon emissions and carbon quotas;

[0029] S54: From this, we can get the carbon trading cost of the comprehensive energy system for:

[0030]

[0031] Where ω is the benchmark price for carbon trading; τ represents the length of the carbon emission interval; and γ represents the price increase.

[0032] Furthermore, the total amount of carbon dioxide captured during the CCS operation at time t is E CCS (t) is calculated as follows:

[0033] P e，CCS (t) = P e，CCS，oper (t)+P CCS，base (t);

[0034]

[0035] E CC (t) = η c e G P GT (t);

[0036]

[0037] In the above formula, P e，CCS (t) is the total energy consumption of CCS; P CCS，base (t) is the CCS benchmark energy consumption; P e，CCS，oper (t) is the CCS operation energy consumption; is the carbon capture efficiency; P GT (t) represents the power generation of the gas turbine at time t; e G Indicates the carbon emission intensity per unit of electricity; ηc represents the carbon capture efficiency of CCS; Indicates the maximum energy consumption of CCS.

[0038] Furthermore, the calculation formula for calculating the shiftable power change amount and the replaceable demand response power change amount of the integrated energy system is as follows:

[0039]

[0040] In the above formula, is the amount of power variation that can be translated; is the amount of power change that can be substituted for demand response; They are respectively represented as the power of electricity, heat and gas loads participating in the shiftable demand response at time t; They are the upper and lower limits of power for electricity, heat and gas loads to participate in shiftable demand response; are the power of electricity, heat and gas loads participating in alternative demand response at time t respectively; They are the upper and lower power limits of electricity, heat and gas loads participating in alternative demand response.

[0041] Furthermore, when calculating the actual carbon emissions of the integrated energy system, each parameter must satisfy the following balance constraints:

[0042] In order to ensure the balance between power supply and demand, and thus maintain the stability and reliability of the power grid, the power balance equation is:

[0043]

[0044] Where: P w (t) is the wind power generation power at time t; P e，Load (t) is the power load power at time t; is the maximum power purchase at time t; is the maximum power generation of the wind turbine; are the charging and discharging power of the storage device at time t; P e，CHP (t) represents the electrical output power of the gas turbine and organic Rankine cycle power generation cogeneration system; P e，P2G (t) is the electric power input to the gas equipment during period t; P e，CCS (t) is the total energy consumption of CCS;

[0045] In order to achieve a balance between natural gas supply and demand and thus ensure the uninterrupted operation of the pipeline network, the natural gas balance equation is as follows:

[0046]

[0047] Where: P g，P2G (t) is the natural gas power output by P2G; Pg，Load (t) is the natural gas load power at time t; is the maximum gas purchasing power at time t; are the charging and deflating powers of the gas storage device at time t; P g，GT (t) is the natural gas power input to the gas turbine during period t; P g，GB (t) is the gas consumption of the gas boiler at time t;

[0048] In order to ensure the balance between heat supply and demand, and thus maintain the smooth operation of the district heating network, the heat balance equation is as follows:

[0049]

[0050] Where: P h，Load (t) is the heat load power at time t; are the charging and releasing powers of the heat storage device at time t; P h，CHP (t) represents the thermal output power of the gas turbine and the organic Rankine cycle heat and power cogeneration system; P h，GB (t) is the output thermal power of the gas boiler at time t;

[0051] In order to achieve a harmonious balance between hydrogen supply and demand, thereby ensuring seamless operation of hydrogen infrastructure, the hydrogen balance equation is as follows:

[0052]

[0053] in: is the hydrogen purchasing power at time t; is the maximum hydrogen purchasing power at time t; They are the charging and discharging power of the hydrogen storage device at time t respectively; It is the hydrogen energy power input from the electrolyzer EL to the hydrogen fuel cell HFC.

[0054] Furthermore, the natural gas power P output by the P2G g，P2G (t), the specific calculation formula is as follows:

[0055] P g，P2G (t) = η P2G P e，P2G (t);

[0056]

[0057] Where P e，P2G (t) is the electric power input to P2G during period t; P g，P2G (t) is the power of natural gas output by P2G during period t; η P2G is the conversion efficiency of P2G; V Gas(t) are the volumes of CO2 absorbed and natural gas produced by P2G during period t; H k is the calorific value of natural gas; is the mass of CO2 absorbed by P2G during period t; is the density of CO2.

[0058] A dispatching control system for an industrial integrated energy system, used to implement the dispatching control method for an industrial integrated energy system as described above, comprising:

[0059] A first acquisition module, which is used to obtain the electricity purchase price and the amount of electricity purchased by the comprehensive energy system;

[0060] The second acquisition module is used to obtain the gas purchase price and gas purchase volume of the integrated energy system;

[0061] The third acquisition module is used to obtain the purchase price and purchase amount of hydrogen for the integrated energy system;

[0062] A fourth acquisition module is used to obtain the wind abandonment penalty price and wind abandonment power of the integrated energy system;

[0063] A first calculation module, which is used to calculate the carbon trading cost of the comprehensive energy system;

[0064] A second calculation module is used to calculate the shiftable power variation and the replaceable demand response power variation of the integrated energy system;

[0065] The third calculation module is used to establish a dispatching and control model for the integrated energy system and minimize the total cost of calculating the system operation based on the above data.

[0066] A computer-readable storage medium includes a stored program, wherein when the program is executed, the device where the computer-readable storage medium is located is controlled to execute the scheduling control method for an industrial integrated energy system as described above.

[0067] A processor is used to run a program, wherein the program, when running, executes the scheduling control method for an industrial integrated energy system as described above.

[0068] Compared with the prior art, the present invention has the following beneficial effects:

[0069] 1. A comprehensive carbon emission measurement method can be used to accurately estimate the capacity of IES, which helps to make informed decisions on carbon emission reduction strategies for industrial parks. At the same time, an optimized energy scheduling strategy is proposed to minimize carbon emissions while maximizing energy efficiency. BRIEF DESCRIPTION OF THE DRAWINGS

[0070] In order to more clearly illustrate the specific embodiments of the present invention or the technical solutions in the prior art, the following is a brief introduction to the drawings required for the specific embodiments or the description of the prior art. In all the drawings, similar elements or parts are generally identified by similar reference numerals. In the drawings, the elements or parts are not necessarily drawn according to the actual scale.

[0071] Figure 1 A schematic flow chart of a dispatching control method for an industrial integrated energy system according to the present invention;

[0072] Figure 2 A schematic diagram of carbon emissions of different components in the IES integrated energy system of the present invention;

[0073] Figure 3 This is the Sankey diagram of carbon emission flow for the IES integrated energy system. DETAILED DESCRIPTION

[0074] The following will be combined with the drawings in the embodiments of the present invention to clearly and completely describe the technical solutions in the embodiments of the present invention. Obviously, the described embodiments are part of the embodiments of the present invention, not all of the embodiments. Based on the embodiments of the present invention, all other embodiments obtained by ordinary technicians in this field without creative work are within the scope of protection of the present invention.

[0075] It should be understood that when used in this specification and the appended claims, the terms "include" and "comprises" indicate the presence of described features, integers, steps, operations, elements and / or components, but do not exclude the presence or addition of one or more other features, integers, steps, operations, elements, components and / or combinations thereof.

[0076] It should also be understood that the terms used in the present specification are only for the purpose of describing specific embodiments and are not intended to limit the present invention. As used in the present specification and the appended claims, unless the context clearly indicates otherwise, the singular forms "a", "an" and "the" are intended to include plural forms.

[0077] It should be further understood that the term "and / or" used in the present description and the appended claims refers to any and all possible combinations of one or more of the associated listed items, and includes these combinations.

[0078] Embodiment 1

[0079] Industrial-based IES (Integrated Energy System) includes various components such as power-to-gas, hydrogen fuel cells, gas turbines, organic Rankine cycle, grid batteries, energy storage, etc. The synergistic advantages of all these components in industrial-based IES can improve energy utilization, reduce emissions, and improve energy efficiency.

[0080] like Figure 1 The figure shows a flow chart of a dispatching control method for an industrial integrated energy system, including:

[0081] S1: Obtain the electricity purchase price and electricity purchase amount of the integrated energy system;

[0082] S2: Obtain the gas purchase price and gas purchase volume of the integrated energy system;

[0083] S3: Obtain the purchase price and purchase amount of hydrogen for the integrated energy system;

[0084] S4: Obtain the wind curtailment penalty price and wind curtailment power of the integrated energy system;

[0085] S5: Calculate the carbon trading costs of the integrated energy system;

[0086] S6: Calculate the shiftable power change of the integrated energy system and the replaceable demand response power change;

[0087] S7: Establish a dispatching and control model for the integrated energy system, and minimize the total system operation cost based on the above data.

[0088] Furthermore, the dispatching control model for establishing the integrated energy system is as follows:

[0089]

[0090] In the above formula, P e，buy (t), P g，buy (t) They represent the electricity, gas and hydrogen purchases at time t respectively; c e (t) is the electricity purchase price at time t; c g (t) is the gas purchase price at time t; is the price of hydrogen purchased at time t; cur is the penalty price for wind curtailment at time t; is the abandoned wind power; For carbon trading costs; are the shiftable and substitutable demand response prices at time t, respectively; is the amount of power variation that can be translated; is the amount of alternative demand response power change.

[0091] Furthermore, the calculation of the carbon trading cost of the integrated energy system includes the following steps:

[0092] S51: Carbon quota is a mechanism to limit corporate carbon dioxide emissions by trading quotas. This policy is used to incentivize emission reduction and promote the development of a low-carbon economy. There are many methods for allocating carbon quotas, including historical allocation, auction allocation, benchmark allocation, and mixed allocation. Among the above carbon quota allocation methods, this study applies the benchmark allocation method to the IES industrial park. The benchmark allocation method sets emission standards for each industry and allocates carbon quotas based on their emissions. The carbon quota for the integrated energy system is calculated as follows:

[0093] E IES =E e，buy +E GT +E GB +E g，Load ;in,

[0094]

[0095] In the above formula, E IES 、E e，buy 、E GT 、E GB 、E g，Load They represent the carbon quota of the integrated energy system and the carbon quota caused by electricity purchased from the main grid, gas turbines, gas boilers and gas loads; μ e , μ h , μ g，Load represents the carbon quota coefficient calculated per unit of electricity, per unit of thermal power and per unit of gas load; P e，buy (t), P e，GT (t), P h，GT (t), P h，GB (t), P g，Load (t) represents the electricity purchase amount, the electric power output of the gas turbine, the thermal power output of the gas turbine, the thermal power output of the gas boiler and the gas load demand at time t respectively;

[0096] S52: Calculate the actual carbon emissions of the integrated energy system using the following formula:

[0097] E IES，a =E e，buy，a +E GT，a +E GB，a +E g，Load，a -E CCS，a ;

[0098]

[0099] In the above formula, E IES，a 、E e，buy，a 、E GT，a 、E GB，a 、E g，Load，aThey represent the actual carbon emissions of the integrated energy system and the actual carbon emissions caused by purchasing electricity from the main grid, gas turbines, gas boilers and gas loads; E CCS，a is the amount of carbon actually captured; They are the first, second and third calculation parameters of the actual carbon emissions of coal-fired power units; are the actual carbon quota coefficients calculated based on unit electricity, unit heat energy, and unit gas load consumption; E CCS (t) is the total amount of CO2 captured during CCS operation at time t;

[0100] S53: The actual carbon difference of the comprehensive energy system is:

[0101] E IES，t =E IES，a -E IES ;

[0102] In the above formula, E IES，t Represents the carbon emission gap between actual carbon emissions and carbon quotas;

[0103] S54: The tiered carbon trading mechanism adopts different trading methods according to the size of corporate emissions. Companies with lower emissions can obtain additional carbon emission quotas, while companies with higher emissions need to purchase more carbon emission rights to offset excess emissions. This tiered trading method encourages companies with lower emissions to continue to reduce their carbon footprint, thereby reducing the overall carbon emissions of the business sector. The carbon trading cost of the integrated energy system can be obtained from this for:

[0104]

[0105] Where ω is the benchmark price for carbon trading; τ represents the length of the carbon emission interval; and γ represents the price increase.

[0106] Furthermore, the total amount of carbon dioxide captured during the operation of COS at time t, ECCS(t), is calculated as follows:

[0107] P e，CCS (t) = P e，CCS，oper (t)+P CCS，base (t);

[0108]

[0109] E CC (t) = η c e G P GT (t);

[0110]

[0111] In the above formula, P e，CCS (t) is the total energy consumption of CCS; P CCS，base (t) is the CCS benchmark energy consumption; P e，COs，oper (t) is the CCS operation energy consumption; is the carbon capture efficiency; P GT (t) represents the power generation of the gas turbine at time t; e G Indicates the carbon emission intensity per unit of electricity; η c represents the carbon capture efficiency of CCS; Indicates the maximum energy consumption of CCS.

[0112] Carbon trading is a market-based carbon emission rights trading mechanism that aims to promote and reduce greenhouse gas emissions and address climate change through economic means. The basic principle of carbon trading is to treat carbon dioxide emissions as a tradable resource, and companies or countries can manage their emission levels by buying or selling carbon credits. Recently, carbon capture and storage (CCS) technology has brought great benefits to existing power plants, which can capture carbon dioxide from the atmosphere and use it for various purposes or store it to reduce greenhouse gas emissions and mitigate climate change.

[0113] Furthermore, demand response plays a vital role in reducing carbon emissions by flexibly adjusting energy consumption patterns and time. It helps to alleviate peak loads in industrial parks, thereby reducing dependence on traditional power generation facilities. On a horizontal time scale, demand response plays an important role in matching power system demand through load shifting and interruptible loads. On the alternative demand response scale, the supply characteristics of electricity, heating and gas loads can be flexibly selected to better adapt to demand fluctuations. The calculation formula for calculating the shiftable power change and the alternative demand response power change of the integrated energy system is as follows:

[0114]

[0115]

[0116] In the above formula, is the amount of power variation that can be translated; is the amount of power change that can be substituted for demand response; They are respectively represented as the power of electricity, heat and gas loads participating in the shiftable demand response at time t; They are the upper and lower limits of power for electricity, heat and gas loads to participate in shiftable demand response; are the power of electricity, heat and gas loads participating in alternative demand response at time t respectively; They are the upper and lower power limits of electricity, heat and gas loads participating in alternative demand response.

[0117] Furthermore, when calculating the actual carbon emissions of the integrated energy system, each parameter must satisfy the following balance constraints:

[0118] In order to ensure the balance between power supply and demand, and thus maintain the stability and reliability of the power grid, the power balance equation is:

[0119]

[0120] Where: P w (t) is the wind power generation power at time t; P e，Load (t) is the power load power at time t; is the maximum power purchase at time t; is the maximum power generation of the wind turbine; are the charging and discharging power of the storage device at time t; P e，CHP (t) represents the electrical output power of the gas turbine and organic Rankine cycle power generation cogeneration system; P e，P2G (t) is the electric power input to the gas equipment during period t; P e，CCS (t) is the total energy consumption of CCS;

[0121] In order to achieve a balance between natural gas supply and demand and thus ensure the uninterrupted operation of the pipeline network, the natural gas balance equation is as follows:

[0122]

[0123] Where: P g，P2G (t) is the natural gas power output by P2G; P g，Load (t) is the natural gas load power at time t; is the maximum gas purchasing power at time t; are the charging and deflating powers of the gas storage device at time t; P g，GT (t) is the natural gas power input to the gas turbine during period t; P g，GB (t) is the gas consumption of the gas boiler at time t;

[0124] In order to ensure the balance between heat supply and demand, and thus maintain the smooth operation of the district heating network, the heat balance equation is as follows:

[0125]

[0126] Where: P h，Load (t) is the heat load power at time t; are the charging and releasing powers of the heat storage device at time t; P h，CHP (t) represents the thermal output power of the gas turbine and the organic Rankine cycle heat and power cogeneration system; P h，GB (t) is the output thermal power of the gas boiler at time t;

[0127] In order to achieve a harmonious balance between hydrogen supply and demand, thereby ensuring seamless operation of hydrogen infrastructure, the hydrogen balance equation is as follows:

[0128]

[0129] in: is the hydrogen purchasing power at time t; is the maximum hydrogen purchasing power at time t; They are the charging and discharging power of the hydrogen storage device at time t respectively; It is the hydrogen energy power input from the electrolyzer EL to the hydrogen fuel cell HFC.

[0130] Furthermore, the P2G (power-to-gas) technology converts the surplus electricity generated by renewable energy into synthetic natural gas, and the CO2 emissions of the P2G device are calculated. The natural gas power P output by the P2G g，P2G (t), the specific calculation formula is as follows:

[0131] P g，P2G (t) = η P2G P e，P2G (t);

[0132]

[0133] Where P e，P2G (t) is the electric power input to P2G during period t; P g，P2G (t) is the power of natural gas output by P2G during period t; η P2G is the conversion efficiency of P2G; V Gas (t) are the volumes of CO2 absorbed and natural gas produced by P2G during period t; H k is the calorific value of natural gas; is the mass of CO2 absorbed by P2G during period t; is the density of CO2.

[0134] Furthermore, hydrogen fuel cells use hydrogen as fuel to generate electricity through an electrochemical process:

[0135]

[0136] Where: P is the hydrogen energy power input from the electrolyzer EL to the hydrogen fuel cell HFC; e，HFC (t), P h，HFC (t) are the electrical and thermal power output of HFC respectively; are the efficiency of HFC conversion into electricity and heat, respectively; They are the upper and lower limits of hydrogen power input into HFC respectively; They are the upper and lower limits of HFC hydrogen energy power climbing respectively.

[0137] Furthermore, the gas turbine converts natural gas into mechanical energy, which is then used to generate electricity. The specific calculation formula for the power generation and heat generation of the GT (gas turbine) device is as follows:

[0138]

[0139] Where: P g，GT (t) is the natural gas power input to GT during period t; P e，GT (t), P h，GT (t) are the electrical and thermal power output of GT in period t, respectively; They are the power generation and heat generation efficiencies of the GT equipment; are the upper and lower limits of GT input power respectively; They are the upper and lower limits of GT climbing power respectively.

[0140] Furthermore, the organic Rankine cycle uses organic liquids with a boiling point lower than that of water to generate electricity from waste heat sources. The specific calculation formula for the heat generation power of the hydrogen fuel cell HFC device is as follows:

[0141] P h，HFC (t) = P h，HFC，WHB (t)+P h，HFC，ORC (t);

[0142] P h，GT (t) = P h，GT，WHB (t)+P h，GT，ORC (t);

[0143] P h，in，WHB (t) = P h，HFC，WHB (t)+P h，GT，WHB (t);

[0144] P h，in，ORC (t) = P h，HFC，ORC (t)+P h，GT，ORC (t);

[0145] P h，ex，WHB (t) = η WHB P h，in，WHB (t);

[0146] P e，ex，ORC (t) = η ORC P e，in，ORC (t);

[0147]

[0148] Where: P h，HFC，WHB (t), P h，HFC，ORC(t) are the thermal power output of the hydrogen fuel cell HFC to the WHB (waste heat boiler) and ORC (organic Rankine cycle) during period t; P h，GT，WHB (t), P h，GT，ORC (t) are the thermal power output from GT to WHB and PRC during period t; P h，in，WHB (t), P h，in，ORC (t) are the input thermal power of WHB and ORC in period t; P h，ex，WHB (t) is the output thermal power of WHB during period t; P h，ex，ORC (t) is the output power of ORC during period t; η WHB , η ORC are the transformation efficiencies of WHB and ORC, respectively; They are the upper and lower limits of the input heat power of WHB respectively; are the upper and lower limits of the input thermal power of ORC respectively; They are the upper and lower limits of the climbing power of WHB respectively; They are the upper and lower limits of ORC climbing power respectively.

[0149] Among them, WHB stands for waste heat boiler; ORC stands for organic Rankine cycle.

[0150] Furthermore, the thermal power of the GT (gas turbine) is utilized through a combination of a waste heat boiler and an organic Rankine cycle. The specific calculation formula for the thermal power of the GT device is as follows:

[0151] P h，GT (t) = P h，WHB (t)+P h，ORC (t);

[0152] P e，ORC (t) = η ORC P h，ORC (t);

[0153] P h，WHB (t) = (1-η WHB ) h，WHB (t);

[0154]

[0155] Where: P h，WHB (t), P h，ORC (t) represents the thermal power input from GT device to WHB and ORC respectively; P h，WHB (t), P e，ORC (t) represents the output thermal power of WHB and ORC respectively; η ORC represents the conversion efficiency of ORC; η WHB It represents the heat loss of WHB; Respectively represent the upper and lower limits of the input power of ORC; They represent the upper and lower limits of ORC climbing respectively.

[0156] Furthermore, by combining gas turbine and organic Rankine cycle power generation, the overall power output of the heat and power unit is realized. The power generation power of the CHP (heat and power unit composed of gas turbine and organic Rankine cycle) device is calculated as follows:

[0157] P e，CHP (t) = P e，GT (t)+P e，ORC (t);

[0158] P h，CHP (t) = P h，WHB (t);

[0159] Where: P h，CHP (t), P e，CHP (t) represent the thermal and electrical output power of CHP, respectively.

[0160] Furthermore, the gas boiler generates heat by burning natural gas to provide thermal heating and hot water supply. The heat generation power of the GB device is calculated by formula (46). The specific calculation formula is as follows:

[0161] P h，GB (t) = η GB P g，GB (t);

[0162]

[0163] Where: η GB is the heat generation efficiency of the gas boiler; P g，GB (t) is the gas consumption of the gas boiler at time t; P is the lower and upper limits of the output thermal power of the gas boiler; h，GB (t) is the output thermal power of the gas boiler at time t; They are the lower and upper limits of gas boiler power ramp-up respectively.

[0164] Furthermore, energy storage adopts a multifunctional energy storage system, including electricity, heat, gas and hydrogen storage. The charging and discharging power and SOC of the energy storage system are calculated as follows:

[0165]

[0166] E j (1) = E j (T).

[0167]

[0168] Where: is the charging power of the j-th energy storage device at time t; is the energy release power of the j-th energy storage device at time t; is the maximum power of a single charge and discharge of the j-th energy storage device; are binary variables, which are the charging and discharging state parameters of the j-th energy storage device at time t; Indicates that the energy storage device is storing energy; Indicates that the energy storage device is discharging energy; are the charging and discharging efficiencies of the j-th energy storage device; E j (t) is the capacity of the j-th energy storage device at time t; are the upper and lower limits of the capacity of the j-th energy storage device; P j (t) is the power of the jth energy storage device at time t; E j (1) E j (T) are the capacities of the j-th energy storage device at the initial and end times, respectively.

[0169] In order to verify the effectiveness of the design, an industrial park with a multi-energy coupling system was selected for verification. Table 1 lists the technical parameters of each component of the system. The industrial park contains various types of equipment and energy consumption, such as electricity, heat, gas, and hydrogen. Through careful consideration of the technical parameters of each system component, the performance of the IES (integrated energy system) in the industrial park can be evaluated, including energy supply and demand balance, energy utilization efficiency, and cost-effectiveness.

[0170] Table 1 Technical parameters of different components

[0171]

[0172] Figure 2 (a) shows the change in carbon emission intensity of grid power generation (thermal power generation) within 24 hours. The vertical axis represents the carbon emission intensity in kilograms, and the horizontal axis represents different time points in a day, with a total of 24 points, each representing an hour, and carbon emissions reflect changes in usage. The highest carbon emission value is about 2,500 kilograms, and the lowest value is about 1,000 kilograms. The lowest carbon emissions occur at the beginning and end of each cycle, indicating that the demand for electricity is low during this period, resulting in a reduction in power generation from thermal power plants, thereby reducing carbon emissions.

[0173] The carbon emission intensity of gas turbines and gas boilers changes within 24 hours. Figure 2 As shown in (b), the carbon emission intensity of gas turbines fluctuates more compared with gas boilers. The peak carbon emission intensity of gas turbines is 244.9 kg, and the lowest point is about 122.5 kg. In contrast, the carbon emission intensity of gas boilers is relatively stable. Figure 2(c) shows the carbon consumption of the power-to-gas equipment and the carbon storage of the carbon capture and storage equipment within 24 hours. During 12:00-14:00 and 19:00-20:00, the carbon consumption of the power-to-gas equipment and the carbon storage of the carbon capture and storage equipment showed relatively opposite trends, which means that when the energy conversion efficiency of the power-to-gas equipment increases, the carbon storage capacity of the carbon capture and storage equipment will decrease, and vice versa. Figure 2 (d) shows the impact of two different types of demand response, namely vertical and horizontal demand response dispatch, on carbon emissions and absorption. Horizontal demand response generates carbon emissions during most hours, but also during low load hours in the early morning and at night, which may be related to load shifting or reduction measures, such as adjusting the operating time to avoid high carbon emission hours. Vertical demand response shows relatively stable carbon emissions throughout the day. This may be because vertical demand response focuses more on energy substitution rather than reducing overall energy consumption.

[0174] Figure 3 The Sankey diagram shown details the carbon flow of the IES in the industrial park, including carbon emissions and carbon sequestration at each stage. The largest carbon emission is generated by purchasing electricity from the main grid, with a total of 43,488 kg. This means that reducing carbon emissions from electricity purchased from the main grid is crucial to reducing the overall carbon emissions of the industrial park. Carbon capture technology plays an important role in this system, capturing 6,933 kg of carbon. The main focus should be on reducing carbon emissions from purchasing electricity from the main grid, which can be achieved by procuring renewable energy or investing in sustainable energy projects.

[0175] Therefore, this embodiment scheme adopts a comprehensive carbon emission measurement method to accurately estimate the capacity of IES, which helps to make wise decisions on the carbon emission reduction strategy of the industrial park. At the same time, an optimized energy scheduling strategy is proposed to minimize carbon emissions while maximizing energy utilization efficiency.

[0176] This study comprehensively and systematically analyzes the carbon emission process of the energy system of the integrated energy industrial park, and deeply explores the source and flow of carbon emissions. At the same time, combined with data simulation, real and reliable data support is obtained to ensure the scientificity and accuracy of the research results. The main contributions of this paper can be summarized in two aspects: (1) By considering the multi-energy utilization and conversion process, a comprehensive carbon emission measurement method is adopted to accurately estimate the capacity of IES. (2) An optimized energy scheduling strategy is proposed to minimize carbon emissions while maximizing energy utilization efficiency.

[0177] Embodiment 2

[0178] A dispatching control system for an industrial integrated energy system, used to implement the dispatching control method for an industrial integrated energy system as described above, comprising:

[0179] A first acquisition module, which is used to obtain the electricity purchase price and the amount of electricity purchased by the comprehensive energy system;

[0180] The second acquisition module is used to obtain the gas purchase price and gas purchase volume of the integrated energy system;

[0181] The third acquisition module is used to obtain the purchase price and purchase amount of hydrogen for the integrated energy system;

[0182] A fourth acquisition module is used to obtain the wind abandonment penalty price and wind abandonment power of the integrated energy system;

[0183] A first calculation module, which is used to calculate the carbon trading cost of the comprehensive energy system;

[0184] A second calculation module is used to calculate the shiftable power variation and the replaceable demand response power variation of the integrated energy system;

[0185] The third calculation module is used to establish a dispatching and control model for the integrated energy system and minimize the total cost of calculating the system operation based on the above data.

[0186] Embodiment 3

[0187] A computer-readable storage medium includes a stored program, wherein when the program is executed, the device where the computer-readable storage medium is located is controlled to execute the scheduling control method for an industrial integrated energy system as described above.

[0188] Embodiment 4

[0189] A processor is used to run a program, wherein the program, when running, executes the scheduling control method for an industrial integrated energy system as described above.

[0190] The present application provides a dispatching and control method for an industrial integrated energy system, including: S1: obtaining the electricity purchase price and the amount of electricity purchased by the integrated energy system; S2: obtaining the gas purchase price and the amount of gas purchased by the integrated energy system; S3: obtaining the hydrogen purchase price and the amount of hydrogen purchased by the integrated energy system; S4: obtaining the wind abandonment penalty price and wind abandonment power of the integrated energy system; S5: calculating the carbon trading cost of the integrated energy system; S6: calculating the shiftable power change and the replaceable demand response power change of the integrated energy system; S7: establishing a dispatching and control model for the integrated energy system, and minimizing the total cost of system operation based on the above data. Through a comprehensive and systematic analysis of the carbon emission process of the energy system of the integrated energy industrial park, in-depth exploration of the source and flow of carbon emissions, and the use of a comprehensive carbon emission measurement method, the capacity of IES is accurately estimated, while maximizing energy efficiency and minimizing carbon emissions.

[0191] Those of ordinary skill in the art will appreciate that the units of each example described in conjunction with the embodiments disclosed herein can be implemented in electronic hardware, computer software, or a combination of the two. In order to clearly illustrate the interchangeability of hardware and software, the composition of each example has been generally described in terms of function in the above description. Whether these functions are performed in hardware or software depends on the specific application and design constraints of the technical solution. Professional and technical personnel can use different methods to implement the described functions for each specific application, but such implementation should not be considered to be beyond the scope of the present invention.

[0192] In the embodiments provided by the present invention, it should be understood that the division of units is only a logical function division, and there may be other division methods in actual implementation, for example, multiple units can be combined into one unit, one unit can be split into multiple units, or some features can be ignored, etc.

[0193] In addition, each functional unit in each embodiment of the present invention may be integrated into one processing unit, or each unit may exist physically separately, or two or more units may be integrated into one unit. The above-mentioned integrated unit may be implemented in the form of hardware or in the form of software functional units.

[0194] If the integrated unit is implemented in the form of a software functional unit and sold or used as an independent product, it can be stored in a computer-readable storage medium. Based on this understanding, the technical solution of the present invention is essentially or the part that contributes to the prior art or all or part of the technical solution can be embodied in the form of a software product, and the computer software product is stored in a storage medium, including several instructions for a computer device (which can be a personal computer, a server or a network device, etc.) to perform all or part of the steps of the method described in each embodiment of the present invention. The aforementioned storage medium includes: U disk, read-only memory (ROM, Read-0nlyMemory), random access memory (RAM, RandomAccessMemory), mobile hard disk, disk or optical disk and other media that can store program codes.

[0195] Finally, it should be noted that the above embodiments are only used to illustrate the technical solutions of the present invention, rather than to limit it. Although the present invention has been described in detail with reference to the aforementioned embodiments, those skilled in the art should understand that they can still modify the technical solutions described in the aforementioned embodiments, or replace some or all of the technical features therein by equivalents. These modifications or replacements do not make the essence of the corresponding technical solutions deviate from the scope of the technical solutions of the embodiments of the present invention, and they should all be included in the scope of the claims and specification of the present invention.

Claims

1. A dispatching control method for an industrial integrated energy system, It is characterized in that include: S1: Obtain the electricity purchase price and electricity purchase amount of the integrated energy system; S2: Obtain the gas purchase price and gas purchase volume of the integrated energy system; S3: Obtain the purchase price and purchase amount of hydrogen for the integrated energy system; S4: Obtain the wind curtailment penalty price and wind curtailment power of the integrated energy system; S5: Calculate the carbon trading costs of the integrated energy system; S6: Calculate the shiftable power change of the integrated energy system and the replaceable demand response power change; S7: Establish a dispatching and control model for the integrated energy system, and minimize the total system operation cost based on the above data.

2. The dispatching control method for an industrial integrated energy system according to claim 1, It is characterized in that The dispatching control model for establishing the comprehensive energy system is as follows: In the above formula, P e,buy (t), P g,buy (t) They represent the electricity, gas and hydrogen purchases at time t respectively; c e (t) is the electricity purchase price at time t; c g (t) is the gas purchase price at time t; is the price of hydrogen purchased at time t; cur is the penalty price for wind curtailment at time t; is the abandoned wind power; For carbon trading costs; are the shiftable and substitutable demand response prices at time t, respectively; is the amount of translational power variation; is the amount of alternative demand response power change.

3. The dispatching control method for an industrial integrated energy system according to claim 2, It is characterized in that The calculation of the carbon trading cost of the comprehensive energy system includes the following steps: S51: Calculate the carbon quota of the integrated energy system using the following formula: E IES = E e，buy +E GT +E GB +E g，Load ; Among them, In the above formula, E IES 、E e，buy 、E GT 、E GB 、E g，Load They represent the carbon quota of the integrated energy system and the carbon quota caused by electricity purchased from the main grid, gas turbines, gas boilers and gas loads; μ e , μ h , μ g，Load represents the carbon quota coefficient calculated per unit of electricity, per unit of thermal power and per unit of gas load; P e，buy (t), P e，GT (t), P h，GT (t), P h，GB (t), P g，Load (t) represents the electricity purchase amount, the electric power output of the gas turbine, the thermal power output of the gas turbine, the thermal power output of the gas boiler and the gas load demand at time t respectively; S52: Calculate the actual carbon emissions of the integrated energy system using the following formula: AND IES，a =And e,buy，a +E GT，a +E GB，a +E g，Load，a -AND CCS，a ; In the above formula, E IES，a 、E e，buy，a 、E GT，a 、E GB，a 、E g，Load，a They represent the actual carbon emissions of the integrated energy system and the actual carbon emissions caused by purchasing electricity from the main grid, gas turbines, gas boilers and gas loads; E CCS，a is the amount of carbon actually captured; They are the first, second and third calculation parameters of the actual carbon emissions of coal-fired power units; are the actual carbon quota coefficients calculated based on unit electricity, unit heat energy, and unit gas load consumption; E CCS (t) is the total amount of CO2 captured during CCS operation at time t; S53: The actual carbon difference of the comprehensive energy system is: AND IES，t =And IES，a -AND IES ; In the above formula, E IES，t Represents the carbon emission gap between actual carbon emissions and carbon quotas; S54: From this, we can get the carbon trading cost of the comprehensive energy system for: Where ω is the benchmark price for carbon trading; τ represents the length of the carbon emission interval; and γ represents the price increase.

4. The dispatching control method for an industrial integrated energy system according to claim 3, It is characterized in that The total amount of carbon dioxide captured during CCS operation at time t is E CCS (t) is calculated as follows: P e,CCS (t)=P e,CCS，oper (t)+P CCS，base (t); AND CC (t)=η c and G P GT (t)? In the above formula, P e，CCS (t) is the total energy consumption of CCS; P CCS，base (t) is the CCS benchmark energy consumption; P e，CCS，oper (t) is the CCS operation energy consumption; is the carbon capture efficiency; P GT (t) represents the power generation of the gas turbine at time t; e G Indicates the carbon emission intensity per unit of electricity; η c represents the carbon capture efficiency of CCS; Indicates the maximum energy consumption of CCS.

5. The dispatching control method for an industrial integrated energy system according to claim 2, It is characterized in that The calculation formula for calculating the shiftable power change of the integrated energy system and the replaceable demand response power change is as follows: In the above formula, is the amount of power variation that can be translated; is the amount of power change that can be substituted for demand response; They are respectively represented as the power of electricity, heat and gas loads participating in the shiftable demand response at time t; They are the upper and lower limits of power for electricity, heat and gas loads to participate in shiftable demand response; are the power of electricity, heat and gas loads participating in alternative demand response at time t respectively; They are the upper and lower power limits of electricity, heat and gas loads participating in alternative demand response.

6. The dispatching control method for an industrial integrated energy system according to claim 3, It is characterized in that When calculating the actual carbon emissions of the integrated energy system, each parameter must satisfy the following balance constraints: In order to ensure the balance between power supply and demand, and thus maintain the stability and reliability of the power grid, the power balance equation is: Where: P w (t) is the wind power generation power at time t; P e，Load (t) is the power load power at time t; is the maximum power purchase at time t; is the maximum power generation of the wind turbine; are the charging and discharging power of the storage device at time t; P e，CHP (t) represents the electrical output power of the gas turbine and organic Rankine cycle power generation cogeneration system; P e，P2G (t) is the electric power input to the gas equipment during period t; P e，CCS (t) is the total energy consumption of CCS; In order to achieve a balance between natural gas supply and demand and thus ensure the uninterrupted operation of the pipeline network, the natural gas balance equation is as follows: Where: P g,P2G (t) is the natural gas power output by P2G; P g，Load (t) is the natural gas load power at time t; is the maximum gas purchasing power at time t; are the charging and deflating powers of the gas storage device at time t; P g，GT (t) is the natural gas power input to the gas turbine during period t; P g，GB (t) is the gas consumption of the gas boiler at time t; In order to ensure the balance between heat supply and demand, and thus maintain the smooth operation of the district heating network, the heat balance equation is as follows: Where: P h，Load (t) is the heat load power at time t; are the charging and releasing powers of the heat storage device at time t; P h，CHP (t) represents the thermal output power of the gas turbine and the organic Rankine cycle heat and power cogeneration system; P h，GB (t) is the output thermal power of the gas boiler at time t; In order to achieve a harmonious balance between hydrogen supply and demand, thereby ensuring seamless operation of hydrogen infrastructure, the hydrogen balance equation is as follows: in: is the hydrogen purchasing power at time t; is the maximum hydrogen purchasing power at time t; They are the charging and discharging power of the hydrogen storage device at time t respectively; It is the hydrogen energy power input from the electrolyzer EL to the hydrogen fuel cell HFC.

7. The dispatching control method for an industrial integrated energy system according to claim 6, It is characterized in that The natural gas power Pg output by the P2G is P2G (t), the specific calculation formula is as follows: P g，P2G (t)=η P2G P e，P2G (t); Where P e，P2G (t) is the electric power input to P2G during period t; P g，P2G (t) is the power of natural gas output by P2G during period t; η P2G is the conversion efficiency of P2G; V Gas (t) are the absorption of CO by P2G in period t 2 and the volume of natural gas produced; H k is the calorific value of natural gas; P2G absorbs CO during period t 2 quality; For CO 2 density.

8. A dispatching control system for an industrial integrated energy system, It is characterized in that A method for implementing the dispatching control method for an industrial integrated energy system according to any one of claims 1 to 7, comprising: A first acquisition module, which is used to obtain the electricity purchase price and the amount of electricity purchased by the comprehensive energy system; The second acquisition module is used to obtain the gas purchase price and gas purchase volume of the integrated energy system; The third acquisition module is used to obtain the purchase price and purchase amount of hydrogen for the integrated energy system; A fourth acquisition module is used to obtain the wind abandonment penalty price and wind abandonment power of the integrated energy system; A first calculation module, which is used to calculate the carbon trading cost of the comprehensive energy system; A second calculation module is used to calculate the shiftable power variation and the replaceable demand response power variation of the integrated energy system; The third calculation module is used to establish a dispatching and control model for the integrated energy system and minimize the total cost of calculating the system operation based on the above data.

9. A computer-readable storage medium, It is characterized in that The computer-readable storage medium includes a stored program, wherein when the program is executed, the device where the computer-readable storage medium is located is controlled to execute the dispatching control method for an industrial integrated energy system as described in any one of claims 1 to 7.

10. A processor, It is characterized in that The processor is used to run a program, wherein the program, when running, executes the dispatching control method for an industrial integrated energy system as described in any one of claims 1 to 7.

Images