IES scheduling method and device based on hydrogen production and low carbon willingness coordination

Abstract

The invention discloses an IES scheduling method and device based on coordinated hydrogen production and low carbon willingness, and relates to the field of comprehensive energy system optimization scheduling. According to a methane cracking hydrogen production method for kallaite hydrogen production, a kallaite hydrogen production mathematical model is constructed; the purposes of coordinating the proportions of green hydrogen and kallaite hydrogen prepared in different time periods and reducing the carbon emission of the system are achieved; analyzing the carbon footprint of the hydrogen-containing comprehensive energy system for supplying energy in a set time period according to the purchase willingness of a user for the zero-carbon energy, and dividing the energy supplied to the load side by the hydrogen-containing comprehensive energy system into zero-carbon electric energy heat energy and common electric energy heat energy; an extra price coefficient is set for the zero-carbon energy, and extra benefits of the zero-carbon energy are obtained; according to a set objective function, constructing a low-carbon economic dispatching model and constraint conditions; and inputting the set parameter data into the low-carbon economic dispatching model to output a hydrogen-containing comprehensive energy system dispatching strategy. The total operation cost and carbon emission of the system can be effectively reduced.

Description

Technical Field

[0001] The present invention relates to the technical field of integrated energy system optimal scheduling, and particularly relates to an IES scheduling method and device based on coordinated hydrogen production and low-carbon willingness. Background Art

[0002] With the increasingly severe global climate change problem, low-carbon economy and sustainable development have become important directions for the development of various countries. As an important means to achieve a low-carbon economy, the optimal scheduling problem of the integrated energy system (IES) has received extensive attention. The integrated energy system is an energy network that integrates multiple energy types (such as electricity, heat, gas, etc.). By efficiently integrating and optimizing various energy devices, such as wind power, energy storage, grid trading, gas turbines, gas boilers, electric chillers, etc., it realizes the coordinated supply of electric load, heat load, and cold load. The optimal scheduling of such a system is of great significance for improving energy utilization efficiency and reducing carbon emissions.

[0003] In the integrated energy system, hydrogen production, as an important energy conversion method, its efficiency and cost have a significant impact on the overall performance of the system. The selection and coordination of hydrogen production methods not only relate to the output and quality of hydrogen energy, but also directly affect the energy utilization efficiency and carbon emission level of the system. Hydrogen can be classified into green hydrogen, turquoise hydrogen, blue hydrogen, and grey hydrogen according to its source. Green hydrogen is hydrogen obtained by electrolyzing water using renewable energy for power generation, and there is no carbon dioxide (CO2) emission during the production process. Turquoise hydrogen is hydrogen produced by methane pyrolysis (MP), and there is also no carbon dioxide emission during the production process. Blue hydrogen is obtained by combining a carbon capture device with a methane steam reforming device. During the process of producing hydrogen by steam reforming methane, carbon capture, utilization, and storage (CCUS) technology is used to reduce CO2 emissions. Grey hydrogen is hydrogen obtained by steam reforming of fossil energy, and has a high carbon emission during the production process. Turquoise hydrogen is hydrogen produced by methane pyrolysis (MP) technology, and there is also no carbon dioxide emission during the production process. Currently, the integrated energy system usually meets the system hydrogen load demand by producing green hydrogen. However, the power consumption for green hydrogen production is relatively high. In order to meet the hydrogen load demand, the system needs to improve the operating efficiency of gas turbines at the peak of the electric load, resulting in a high carbon emission.

[0004] Therefore, in the process of optimal scheduling, it is necessary to fully consider the coordination problem of hydrogen production methods to achieve efficient energy utilization and low-carbon emissions. Summary of the Invention

[0005] To this end, the present invention provides an IES scheduling method and device based on coordinating hydrogen production and low-carbon willingness. A turquoise hydrogen production device is introduced into the integrated energy system, and a mathematical model for producing turquoise hydrogen using a methane cracking device is constructed. By coordinating the production ratios of green hydrogen and turquoise hydrogen in the system at different time periods, on the basis of meeting the system's hydrogen load, the total operating cost and carbon emissions of the system are effectively reduced. At the same time, the willingness of users to purchase zero-carbon energy is also considered, and by formulating a reasonable additional price coefficient for zero-carbon energy, the total operating cost of the system is effectively reduced.

[0006] To achieve the above object, the present invention provides the following technical solutions: An IES scheduling method based on coordinating hydrogen production and low-carbon willingness, including: Construct a mathematical model for producing turquoise hydrogen according to the methane cracking hydrogen production method for producing turquoise hydrogen; By introducing a coordinated hydrogen production method considering green hydrogen and turquoise hydrogen, construct a hydrogen-containing integrated energy system model; According to the willingness of users to purchase zero-carbon energy, by analyzing the carbon footprint of the energy supplied by the hydrogen-containing integrated energy system during the set time period, divide the energy supplied by the hydrogen-containing integrated energy system to the load side into zero-carbon electric energy, zero-carbon thermal energy, ordinary electric energy, and ordinary thermal energy; by setting an additional price coefficient for the zero-carbon energy, obtain the additional income of the zero-carbon energy; According to the set objective function, construct a low-carbon economic scheduling model and constraint conditions; Input the set parameter data into the low-carbon economic scheduling model; through the processing of the low-carbon economic scheduling model, output the scheduling strategy of the hydrogen-containing integrated energy system.

[0007] As a preferred solution of the IES scheduling method based on coordinating hydrogen production and low-carbon willingness, the expression of the mathematical model for producing turquoise hydrogen is: ; In the formula, is the volume of turquoise hydrogen produced by methane cracking in the time period; is the volume of natural gas consumed by the turquoise hydrogen production device in the time period; is the electric energy consumed by the turquoise hydrogen production device in the time period; is the electric energy consumed per unit volume of turquoise hydrogen produced in the time period; is the mass of carbon generated by the methane cracking reaction in the time period; is the efficiency of methane conversion to carbon; and are the upper and lower limits for hydrogen production by the methane cracking equipment during any period; is the low calorific value of carbon; is the low calorific value of methane.

[0008] As an optimal solution of the IES dispatching method based on coordinating hydrogen production and low-carbon willingness, in the hydrogen-integrated energy system model, the electrical load of the hydrogen-integrated energy system is supplied by purchasing electricity from the power grid, hydrogen fuel cells, wind turbines, and gas turbines; The thermal load of the hydrogen-integrated energy system is supplied by gas turbines, hydrogen fuel cells, electric boilers, and gas boilers; The hydrogen load of the hydrogen-integrated energy system is supplied by electrolytic hydrogen production equipment and methane cracking equipment; The energy storage equipment of the hydrogen-integrated energy system includes electricity storage equipment, heat storage equipment, and hydrogen storage equipment.

[0009] As an optimal solution of the IES dispatching method based on coordinating hydrogen production and low-carbon willingness, the mathematical expression for the user's willingness to purchase zero-carbon energy is: ; ; ; In the formula, are users who are very willing; are users who are relatively willing; are neutral users; are users who are relatively unwilling; are users who are completely unwilling; is the user willingness distribution matrix after the k-th iteration; is the willingness of a single user to purchase zero-carbon energy.

[0010] As an optimal solution of the IES dispatching method based on coordinating hydrogen production and low-carbon willingness, the mathematical model expression for dividing the energy supplied by the hydrogen-integrated energy system to the load side into the zero-carbon electric energy, the zero-carbon thermal energy, the ordinary electric energy, and the ordinary thermal energy is: ; In the formula, and are respectively the electrical load and the thermal load of the system during the and are respectively the zero-carbon electric energy and the zero-carbon thermal energy produced by the system during the is the volume of natural gas consumed by the gas turbine during the is the volume of natural gas consumed by the gas boiler during the is the main online power purchase quantity during a time period; is the electrical conversion efficiency of the gas turbine; is the thermal conversion efficiency of the gas turbine, is the thermal conversion efficiency of the gas boiler.

[0011] As an optimal solution of the IES scheduling method based on coordinated hydrogen production and low-carbon willingness, the objective function of the low-carbon economic scheduling model is: ; In the formula, is the total system operation cost; is the main online power purchase cost; is the equipment operation and maintenance cost; is the stepped carbon trading cost, is the carbon tax cost, is the green certificate income, is the wind curtailment penalty cost, is the additional income of zero-carbon energy; The constraint conditions include: wind power generation constraint, coupling constraint, power balance constraint and energy storage constraint.

[0012] The present invention also provides an IES scheduling device based on coordinated hydrogen production and low-carbon willingness. Based on the above IES scheduling method based on coordinated hydrogen production and low-carbon willingness, it includes: A mathematical model construction module for producing turquoise hydrogen, which is used to construct a mathematical model for producing turquoise hydrogen according to the methane cracking hydrogen production method for producing turquoise hydrogen; A hydrogen-containing integrated energy system model construction module, which is used to construct a hydrogen-containing integrated energy system model by introducing a coordinated hydrogen production method considering green hydrogen and turquoise hydrogen; An energy division module, which is used to divide the energy supplied by the hydrogen-containing integrated energy system to the load side into zero-carbon electric energy, zero-carbon thermal energy, ordinary electric energy and ordinary thermal energy by analyzing the carbon footprint of the energy supplied by the hydrogen-containing integrated energy system during a set time period according to the user's willingness to purchase zero-carbon energy; by setting an additional price coefficient for the zero-carbon energy, obtaining the additional income of the zero-carbon energy; A low-carbon economic scheduling model and constraint condition construction module, which is used to construct a low-carbon economic scheduling model and constraint conditions according to a set objective function; A scheduling strategy output module, which is used to input set parameter data into the low-carbon economic scheduling model; through processing by the low-carbon economic scheduling model, outputting a scheduling strategy for the hydrogen-containing integrated energy system.

[0013] As an optimal solution for the IES scheduling device based on coordinating hydrogen production and low-carbon willingness, in the turquoise hydrogen production mathematical model construction module, the expression of the turquoise hydrogen production mathematical model is: ; In the formula, is The volume of turquoise hydrogen produced by methane cracking in the is The volume of natural gas consumed by the turquoise hydrogen production equipment in the is The electric energy consumed by the turquoise hydrogen production equipment in the is The electric energy consumed per unit volume of turquoise hydrogen produced in the is The mass of carbon produced by the methane cracking reaction in the is the efficiency of methane conversion to carbon; is the efficiency of methane conversion to hydrogen; and are the upper and lower limits of hydrogen production by the methane cracking equipment at any time; is the lower calorific value of carbon; is the lower calorific value of methane.

[0014] As an optimal solution for the IES scheduling device based on coordinating hydrogen production and low-carbon willingness, in the hydrogen-containing integrated energy system model construction module, in the hydrogen-containing integrated energy system model, the electrical load of the hydrogen-containing integrated energy system is supplied by power grid power purchase, hydrogen fuel cells, wind turbines, and gas turbines; The thermal load of the hydrogen-containing integrated energy system is supplied by gas turbines, hydrogen fuel cells, electric boilers, and gas boilers; The hydrogen load of the hydrogen-containing integrated energy system is supplied by electrolytic hydrogen production equipment and methane cracking equipment; The energy storage equipment of the hydrogen-containing integrated energy system includes electricity storage equipment, heat storage equipment, and hydrogen storage equipment.

[0015] As an optimal solution for the IES scheduling device based on coordinating hydrogen production and low-carbon willingness, in the energy division module, the mathematical expression of the user's willingness to purchase zero-carbon energy is: ; ; ; In the formula, is the user who is very willing; is the user who is relatively willing; is the neutral user; is the user who is relatively unwilling; For users who are completely unwilling; is the user willingness distribution matrix after the k-th iteration; is the purchase willingness of a single user for zero-carbon energy.

[0016] As an optimal solution of the IES scheduling device based on coordinated hydrogen production and low-carbon willingness, in the energy division module, the mathematical model expression for dividing the energy supplied by the hydrogen-integrated energy system to the load side into the zero-carbon electric energy, the zero-carbon thermal energy, the ordinary electric energy, and the ordinary thermal energy is: ; In the formula, and are respectively the system's electrical load and thermal load during the period; and are respectively the zero-carbon electric energy and zero-carbon thermal energy produced by the system during the period; is the volume of natural gas consumed by the gas turbine during the period; is the volume of natural gas consumed by the gas boiler during the period; is the main network power purchase volume during the period; is the electrical conversion efficiency of the gas turbine; is the thermal conversion efficiency of the gas turbine, is the thermal conversion efficiency of the gas boiler.

[0017] As an optimal solution of the IES scheduling device based on coordinated hydrogen production and low-carbon willingness, in the low-carbon economic scheduling model and constraint condition construction module, the objective function of the low-carbon economic scheduling model is: ; In the formula, is the total system operation cost; is the main network power purchase cost; is the equipment operation and maintenance cost; is the stepped carbon trading cost, is the carbon tax cost, is the green certificate income, is the curtailment penalty cost, is the additional income of zero-carbon energy; The constraint conditions include: wind power generation constraint, coupling constraint, power balance constraint, and energy storage constraint.

[0018] The present invention has the following advantages: According to the hydrogen production method of methane cracking for producing turquoise hydrogen, the present invention constructs a mathematical model for producing turquoise hydrogen; by introducing a coordinated hydrogen production method considering green hydrogen and turquoise hydrogen, a hydrogen-containing integrated energy system model is constructed; according to the user's willingness to purchase zero-carbon energy, by analyzing the carbon footprint of the energy supplied by the hydrogen-containing integrated energy system during a set period, the energy supplied by the hydrogen-containing integrated energy system to the load side is divided into zero-carbon electric energy, zero-carbon heat energy, ordinary electric energy, and ordinary heat energy; by setting an additional price coefficient for the zero-carbon energy, the additional income of the zero-carbon energy is obtained; according to the set objective function, a low-carbon economic dispatch model and constraints are constructed; the set parameter data is input into the low-carbon economic dispatch model; through the processing of the low-carbon economic dispatch model, the dispatch strategy of the hydrogen-containing integrated energy system is output. The present invention takes into account that the power consumption for producing green hydrogen in the system is relatively high. When the system's electrical load reaches its peak, in order to meet the hydrogen load demand, the output of the gas turbine unit is increased, resulting in a significant increase in the system's carbon emissions. By introducing a turquoise hydrogen production device, the production volume of turquoise hydrogen during the electrical load peak is increased, reducing the power supply pressure of the system during the electrical load peak, and realizing the economic, environmental, safe, and reliable operation of the system. The physical meaning is clear, the method is scientific and reasonable, has strong applicability, and the dispatch strategy is reasonable and effective. The present invention introduces turquoise hydrogen into the integrated energy system, constructs a mathematical model for producing turquoise hydrogen by a methane cracking device, and by adjusting the output of the green hydrogen production device and the turquoise hydrogen production device at different times, on the basis of meeting the hydrogen load, the electrical energy consumed for producing hydrogen during the electrical load peak of the system is reduced, the carbon emissions of the system are reduced, and at the same time, the energy supplied by the system to the load side is divided into zero-carbon energy and ordinary energy, considering the influence of the user's low-carbon willingness, and by adjusting the price difference between zero-carbon energy and ordinary energy, the total operating cost of the system is reduced. BRIEF DESCRIPTION OF THE DRAWINGS

[0019] In order to more clearly illustrate the embodiments of the present invention or the technical solutions in the prior art, the following will briefly introduce the drawings required for use in the description of the embodiments or the prior art. Obviously, the drawings in the following description are only exemplary, and for those of ordinary skill in the art, without creative efforts, other implementation drawings can be obtained according to the provided drawings.

[0020] The structures, ratios, sizes, etc. shown in this specification are only used to cooperate with the content disclosed in the specification for those familiar with this technology to understand and read, and are not used to limit the limited conditions under which the present invention can be implemented. Therefore, they do not have a substantial technical meaning. Any modification of the structure, change in the proportional relationship, or adjustment of the size, without affecting the effects that the present invention can produce and the purposes that can be achieved, should still fall within the scope covered by the technical content disclosed in the present invention.

[0021] Figure 1Schematic diagram of the IES scheduling method based on coordinating hydrogen production and low-carbon willingness provided in Embodiment 1 of the present invention; Figure 2 Schematic diagram of the integrated energy system framework in the IES scheduling method based on coordinating hydrogen production and low-carbon willingness provided in Embodiment 1 of the present invention; Figure 3 Schematic diagram of the turquoise hydrogen production equipment in the IES scheduling method based on coordinating hydrogen production and low-carbon willingness provided in Embodiment 1 of the present invention; Figure 4 Schematic diagram of the distribution of users' zero-carbon energy payment willingness in the IES scheduling method based on coordinating hydrogen production and low-carbon willingness provided in Embodiment 1 of the present invention; Figure 5 Schematic diagram of the change in users' payment willingness in the IES scheduling method based on coordinating hydrogen production and low-carbon willingness provided in Embodiment 1 of the present invention; Figure 6 Schematic diagram of the curve of users' zero-carbon energy payment willingness in the IES scheduling method based on coordinating hydrogen production and low-carbon willingness provided in Embodiment 1 of the present invention; Figure 7 Schematic diagram of the electrical power of the equipment in Scheme 1 in a possible embodiment provided in Embodiment 1 of the present invention, where E1MT is the output of the gas turbine, P1buy is the electricity purchased from the main grid by the system, WT is the output of wind power, E1SSD is the energy release of the energy storage device, E1HFC is the output of the hydrogen fuel cell, E1EL is the energy consumption of the electrolyzer, E1EB is the energy consumption of the electric boiler, E1SSC is the energy charging of the energy storage device, ELoad is the electrical load, and E1MP is the energy consumption of the methane cracking equipment; Figure 8 Schematic diagram of the thermal power of the equipment in Scheme 1 in a possible embodiment provided in Embodiment 1 of the present invention; where H1MT is the output of the gas turbine, H1SSD is the energy release of the heat storage device, H1HFC is the output of the hydrogen fuel cell, H1EB is the output of the electric boiler, H1GB is the output of the gas boiler, H1SSC is the energy charging of the heat storage device, and HLoad is the heat load; Figure 9 Schematic diagram of the gas power of the equipment in Scheme 1 in a possible embodiment provided in Embodiment 1 of the present invention; where G1buy is the gas purchased by the system, G1MT is the gas consumption of the gas turbine, and G1GB is the gas consumption of the gas boiler; Figure 10 Schematic diagram of the hydrogen power of the equipment in Scheme 1 in a possible embodiment provided in Embodiment 1 of the present invention; where H21EL is hydrogen production by electrolyzing water, H21SSD is the energy release of the hydrogen storage device, H21MP is hydrogen production by methane cracking, H21MT is the hydrogen consumption of the gas turbine, H21GB is the hydrogen consumption of the gas boiler, H21SSC is the energy charging of the hydrogen storage device, H2load is the hydrogen load, and H21HFC is the hydrogen consumption of the hydrogen fuel cell; Figure 11 Schematic diagram of the electrical power of the equipment in Scheme 2 of a possible embodiment provided in Embodiment 1 of the present invention; where E2MT is the output of the gas turbine, P2buy is the electricity purchased by the system from the main grid, WT is the output of wind power, E2SSD is the energy release of the energy storage device, E2HFC is the output of the hydrogen fuel cell, E2EL is the energy consumption of the electrolyzer, E2EB is the energy consumption of the electric boiler, E2SSC is the energy charging of the energy storage device, ELoad is the electrical load, and E2MP is the energy consumption of the methane cracking device; Figure 12 Schematic diagram of the thermal power of the equipment in Scheme 2 of a possible embodiment provided in Embodiment 1 of the present invention; where H2MT is the output of the gas turbine, H2SSD is the energy release of the heat storage device, H2HFC is the output of the hydrogen fuel cell, H2EB is the output of the electric boiler, H2GB is the output of the gas boiler, H2SSC is the energy charging of the heat storage device, and HLoad is the thermal load; Figure 13 Schematic diagram of the gas power of the equipment in Scheme 2 of a possible embodiment provided in Embodiment 1 of the present invention; where G2buy is the gas purchased by the system, G2MT is the gas consumption of the gas turbine, G2GB is the gas consumption of the gas boiler, and G2MP is the gas consumption of methane cracking; Figure 14 Schematic diagram of the hydrogen power of the equipment in Scheme 2 of a possible embodiment provided in Embodiment 1 of the present invention; where H22EL is hydrogen production by electrolyzing water, H22SSD is the energy release of the hydrogen storage device, H22MP is hydrogen production by methane cracking, H22MT is the hydrogen consumption of the gas turbine, H22GB is the hydrogen consumption of the gas boiler, H22SSC is the energy charging of the hydrogen storage device, H2load is the hydrogen load, and H22HFC is the hydrogen consumption of the hydrogen fuel cell; Figure 15 Schematic diagram of the electrical power of the equipment in Scheme 3 of a possible embodiment provided in Embodiment 1 of the present invention; where E3MT is the output of the gas turbine, P3buy is the electricity purchased by the system from the main grid, WT is the output of wind power, E3SSD is the energy release of the energy storage device, E3HFC is the output of the hydrogen fuel cell, E3EL is the energy consumption of the electrolyzer, E3EB is the energy consumption of the electric boiler, E3SSC is the energy charging of the energy storage device, ELoad is the electrical load, and E3MP is the energy consumption of the methane cracking device; Figure 16 Schematic diagram of the thermal power of the equipment in Scheme 3 of a possible embodiment provided in Embodiment 1 of the present invention; where H3MT is the output of the gas turbine, H3SSD is the energy release of the heat storage device, H3HFC is the output of the hydrogen fuel cell, H3EB is the output of the electric boiler, H3GB is the output of the gas boiler, H3SSC is the energy charging of the heat storage device, and HLoad is the thermal load; Figure 17Schematic diagram of the gas power of the equipment in the third solution of a possible embodiment provided in Embodiment 1 of the present invention; where G3buy is the system's purchased gas volume, G3MT is the gas consumption of the gas turbine, G3GB is the gas consumption of the gas boiler, and G3MP is the gas consumption for methane cracking; Figure 18 Schematic diagram of the hydrogen power of the equipment in the third solution of a possible embodiment provided in Embodiment 1 of the present invention; where H23EL is hydrogen production by electrolyzing water, H23SSD is the energy release of the hydrogen storage device, H23MP is hydrogen production by methane cracking, H23MT is the hydrogen consumption of the gas turbine, H23GB is the hydrogen consumption of the gas boiler, H23SSC is the energy charging of the hydrogen storage device, H2load is the hydrogen load, and H23HFC is the hydrogen consumption of the hydrogen fuel cell; Figure 19 Schematic diagram of the change in the total operating cost of the system in a possible embodiment provided in Embodiment 1 of the present invention; Figure 20 Schematic diagram of the change in carbon emissions in a possible embodiment provided in Embodiment 1 of the present invention; Figure 21 Schematic diagram of the architecture of the IES scheduling device based on coordinating hydrogen production and low-carbon willingness provided in Embodiment 2 of the present invention. Detailed implementation manners

[0022] The following specific embodiments illustrate the implementation manners of the present invention. Those skilled in the art can easily understand other advantages and effects of the present invention from the content disclosed in this specification. Obviously, the described embodiments are part of the embodiments of the present invention, rather than all of them. All other embodiments obtained by those of ordinary skill in the art based on the embodiments of the present invention without creative efforts shall fall within the protection scope of the present invention.

[0023] Embodiment 1

[0024] See Figure 1 , Embodiment 1 of the present invention provides an IES scheduling method based on coordinating hydrogen production and low-carbon willingness, including the following steps: S1. According to the methane cracking hydrogen production method for producing turquoise hydrogen, construct a mathematical model for producing turquoise hydrogen; S2. By introducing a coordinated hydrogen production method considering green hydrogen and turquoise hydrogen, construct a hydrogen-containing integrated energy system model; S3. According to the user's willingness to purchase zero-carbon energy, by analyzing the carbon footprint of the energy supplied by the hydrogen-containing integrated energy system during a set period, divide the energy supplied by the hydrogen-containing integrated energy system to the load side into zero-carbon electric energy, zero-carbon thermal energy, ordinary electric energy, and ordinary thermal energy; by setting an additional price coefficient for the zero-carbon energy, obtain the additional income of the zero-carbon energy; S4. According to the set objective function, construct a low-carbon economic scheduling model and constraint conditions; S5. Input the set parameter data into the low-carbon economic dispatch model; process it through the low-carbon economic dispatch model and output the dispatch strategy of the hydrogen-integrated energy system.

[0025] In this embodiment, in step S1, according to the methane pyrolysis method for producing turquoise hydrogen, a mathematical model for producing turquoise hydrogen is constructed. Specifically, as Figure 3 shown, the production method of turquoise hydrogen is to send methane into the fluidized bed, and under the action of a catalyst, it is thermally cracked to obtain solid carbon and hydrogen. The specific reaction equation is as follows: ; Compared with the production of green hydrogen, the advantage of producing turquoise hydrogen is that less electric energy is used to produce the same amount of hydrogen, and the production of turquoise hydrogen by methane pyrolysis is an efficient and clean production method, which can reduce carbon emissions. The current mainstream methane pyrolysis hydrogen production technologies can be divided into four types: high-temperature thermal cracking method, catalytic cracking method, plasma cracking method, and molten metal cracking method. The present invention selects the molten metal cracking method to produce turquoise hydrogen. Considering that the temperature of the molten metal cracking method needs to reach above 800 °C, the temperature condition required for methane pyrolysis is provided by electric heating at the fluidized bed equipment. The relevant mathematical model is as follows: ; In the formula, is the volume of turquoise hydrogen produced by methane pyrolysis in the time period; is the volume of natural gas consumed by the turquoise hydrogen production equipment in the time period; is the electric energy consumed by the turquoise hydrogen production equipment in the time period; is the electric energy consumed per unit volume of turquoise hydrogen produced in the time period; is the mass of carbon generated by the methane pyrolysis reaction in the time period; is the efficiency of methane conversion to carbon; is the efficiency of methane conversion to hydrogen; and are the upper and lower limits of hydrogen production by the methane pyrolysis equipment at any time period; is the low calorific value of carbon; is the low calorific value of methane.

[0026] In this embodiment, in step S2, by introducing a coordinated hydrogen production method considering green hydrogen and turquoise hydrogen, a hydrogen-integrated energy system model is constructed. Specifically, as Figure 2As shown, in the hydrogen-integrated energy system model, the electrical load of the hydrogen-integrated energy system is supplied by purchasing electricity from the grid, hydrogen fuel cells, wind turbines, and gas turbines; The thermal load of the hydrogen-integrated energy system is supplied by gas turbines, hydrogen fuel cells, electric boilers, and gas boilers; The hydrogen load of the hydrogen-integrated energy system is supplied by electrolytic hydrogen production equipment and methane cracking equipment; The energy storage equipment of the hydrogen-integrated energy system includes electrical energy storage equipment, thermal energy storage equipment, and hydrogen energy storage equipment.

[0027] Among them, the mathematical model of the electrolyzer is as follows: ; In the formula, is the electrical energy consumed by electrolyzing water in the is the volume of hydrogen gas produced by the electrolyzer in the is the efficiency of converting electrical energy into hydrogen energy in the electrolyzer; is the mass of water consumed by the electrolyzer in the is the mass of water consumed per unit volume of hydrogen produced by the electrolyzer in the and are respectively the upper and lower limits of the power injected into the electrolyzer in the is the lower calorific value of hydrogen.

[0028] Among them, the electric boiler equipment can convert electrical energy into thermal energy, and its model is as follows: ; In the formula, is the electrical energy consumed by the electric boiler in the is the thermal energy produced by the electric boiler in the and are the upper and lower limits of the electrical energy consumed by the electric boiler at any time; is the efficiency of the electric boiler converting electrical energy into thermal energy.

[0029] Among them, the mathematical model of the hydrogen-blended gas turbine is as follows: ; In the formula, is the electrical energy output by the gas turbine in the is the thermal energy output by the gas turbine in the is the volume of natural gas consumed by the gas turbine in the is the volume of hydrogen consumed by the gas turbine during the is the mass of CO2 emitted by the gas turbine during the the mass of CO2 emitted per unit volume of natural gas burned by the gas turbine; is the total power output by the gas turbine during the the electrical energy conversion efficiency of the gas turbine; the thermal energy conversion efficiency of the gas turbine; the maximum ramp power of the gas turbine; and is the upper and lower limits of the electrical power output by the gas turbine during the and is the upper and lower limits of the thermal power output by the gas turbine during the

[0030] Among them, the mathematical model of the hydrogen - blended gas boiler is as follows: ; In the formula, is the thermal energy output by the gas boiler during the is the volume of natural gas consumed by the gas boiler during the is the volume of hydrogen consumed by the gas boiler during the is the mass of CO2 emitted by the gas boiler during the the thermal efficiency of the gas boiler; and is and is the upper and lower limits of the ramp power of the gas boiler;

[0031] Among them, the mathematical model of the hydrogen fuel cell is as follows: ; In the formula, is the volume of hydrogen consumed by the hydrogen fuel cell during the is the electrical power generated by the hydrogen fuel cell during the is the thermal power generated by the hydrogen fuel cell during the the electrical efficiency of the hydrogen fuel cell; the thermal efficiency of the hydrogen fuel cell; and are the upper and lower limits of the hydrogen fuel cell's ramp power during any period; and are the upper and lower limits of the hydrogen volume consumed by the hydrogen fuel cell during any period; and are the upper and lower limits of the thermoelectric ratio of the hydrogen fuel cell during any period.

[0032] Among them, the mathematical model related to the energy storage device is as follows: ; In the formula, is the electric energy in the electricity storage device at time and are respectively the charging and discharging powers of the heat storage device during the and are respectively the charging and heat release powers of the heat storage device during the and are respectively the hydrogen charging and discharging powers of the hydrogen storage device during the and are respectively the self-loss coefficients of the electricity storage, heat storage, and hydrogen storage devices; and are respectively the charging and discharging efficiencies of the electric energy of the electricity storage device; and are respectively the charging and discharging efficiencies of the heat energy of the heat storage device; and are respectively the charging and discharging efficiencies of the hydrogen energy of the hydrogen storage device.

[0033] In this embodiment, in step S3, according to the user's willingness to purchase zero-carbon energy, by analyzing the carbon footprint of the hydrogen-integrated energy system supplying energy during the set period, the energy supplied by the hydrogen-integrated energy system to the load side is divided into zero-carbon electric energy, zero-carbon thermal energy, ordinary electric energy, and ordinary thermal energy; by setting an additional price coefficient for the zero-carbon energy, the additional income of the zero-carbon energy is obtained; Specifically, as Figure 4 shows, different users have different attitudes towards zero-carbon energy, and users with different low-carbon willingness have different willingness to purchase zero-carbon energy. Extreme low-carbon enthusiasts tend to be willing to pay a higher price to purchase zero-carbon energy when faced with zero-carbon energy and ordinary energy. In addition, the willingness of users to pay for zero-carbon electric energy is also affected by the price of zero-carbon electric energy. As the price of zero-carbon electric energy increases, the number of users willing to actively purchase zero-carbon electric energy will gradually decrease.

[0034] Users' willingness to purchase zero-carbon energy can also change due to the decisions of people around them. Users with a high willingness to purchase zero-carbon energy tend to promote their low-carbon concepts to people around them, thus driving those neutral people with ambiguous willingness to purchase to become biased again in their willingness to purchase.

[0035] As Figure 5 shown, considering that at the current price of zero-carbon energy, users who are very willing will actively promote to their neighbors, so the influence factor of them on the people around is set to 2, the influence factor of relatively willing users is set to 1, the influence factor of neutral users is set to 0, the influence factor of relatively unwilling users is set to -1, and the influence factor of completely unwilling users is set to -2. Add up the influence factors of the four users around a neutral user. When the result is greater than 0, the neutral user changes to a relatively willing user; when it is equal to 0, the user remains neutral; when it is less than 0, the neutral user changes to a relatively unwilling user. Repeatedly iterate until the users no longer change. There may still be a small number of neutral users remaining neutral due to the sum of the influence factors of the people around them being 0. Considering that relevant statistics show that when more than half of the people around are willing to pay an extra fee to purchase zero-carbon energy, neutral users will be more inclined to pay an extra fee to purchase zero-carbon energy. Therefore, finally, the people who remain neutral are changed to users who are relatively willing to pay an extra fee for zero-carbon electricity.

[0036] The mathematical expression for users' willingness to purchase zero-carbon energy is: ; ; ; In the formula, is the user who is very willing; is the relatively willing user; is the neutral user; is the relatively unwilling user; is the completely unwilling user; is the user willingness distribution matrix after the k-th iteration; is the willingness of a single user to purchase zero-carbon energy.

[0037] ; At each iteration, if , then sum the four values around . After summing, is converted into the corresponding user willingness.

[0038] In this embodiment, regarding the user's willingness to purchase zero-carbon energy, by analyzing the carbon footprint of the energy supplied by the system at different times, the electric energy and heat energy supplied by the system to the load side are classified into zero-carbon electric energy, zero-carbon heat energy, ordinary electric energy, and ordinary heat energy. The relevant mathematical model is as follows: ; In the formula, and are respectively the system's electrical load and heat load during the period; and are respectively the zero-carbon electric energy and zero-carbon heat energy produced by the system during the period; is the volume of natural gas consumed by the gas turbine during the period; is the volume of natural gas consumed by the gas boiler during the period; is the main online power purchase volume during the period; is the electrical conversion efficiency of the gas turbine; is the heat conversion efficiency of the gas turbine, is the heat conversion efficiency of the gas boiler.

[0039] The system divides the electric energy and heat energy supplied to users in each period into zero-carbon electric energy, ordinary electric energy, zero-carbon heat energy, and non-zero-carbon heat energy. However, considering the impact of price on the user's willingness to pay, not all zero-carbon electric energy and zero-carbon heat energy can be sold at an increased price. It is necessary to consider the proportion of users on the user side who are willing to pay for zero-carbon energy. Therefore, for the pricing of zero-carbon energy, referring to the user's zero-carbon energy willingness-to-pay curve, as Figure 6 shown, the load side is divided into a part that can obtain zero-carbon energy benefits and a part that cannot obtain zero-carbon energy according to the corresponding proportion. The relevant constraints are as follows: ; In the formula, and are respectively the zero-carbon electric energy and zero-carbon heat energy that obtain additional benefits during the period; and are respectively the maximum proportions of the energy supplied to the user side at the additional price coefficients of the current zero-carbon electric energy and zero-carbon heat energy that can be sold at the zero-carbon energy price.

[0040] In this embodiment, the stepped carbon trading mechanism, carbon tax, and green certificate trading mechanism are combined and introduced into the electric-heat-hydrogen integrated energy system as a carbon constraint mechanism.

[0041] The actual carbon emission calculation model is as follows: ; In the formula, is the actual carbon emissions of the system; is the carbon emissions of the main online power purchase; is the carbon emissions of the gas turbine; is the carbon emissions of the gas boiler; and are the carbon emission coefficients of the main online power purchase, gas turbine and gas boiler respectively.

[0042] The main sources of carbon emissions in the system are as follows: the main online power purchase, gas turbine and gas boiler. For the calculation of carbon emission quota for stepped carbon trading, based on the baseline method, considering the part of zero-carbon electricity and zero-carbon heat in the electric energy and heat energy supplied to the load side, it is converted into carbon emission quota according to a certain proportion and added to the calculation of the total carbon emission quota. The relevant mathematical model is as follows: ; In the formula, is the total carbon emission quota of the system; is the carbon emission quota of the main online power purchase; is the carbon emission quota of the gas turbine; is the carbon emission quota of the gas boiler; is the carbon emission quota obtained by the system for producing zero-carbon energy; is the carbon emission quota coefficient corresponding to the main online power purchase; is the carbon emission quota coefficient corresponding to the gas turbine; is the carbon emission quota coefficient corresponding to the gas boiler.

[0043] The calculation model of the carbon emission trading volume actually participating in the carbon trading market is as follows: ; In the formula, is the carbon emission trading volume.

[0044] The stepped carbon trading mechanism encourages enterprises to reduce carbon emissions by setting different carbon emission quotas and charging stepped fees for the part exceeding the quota. Different from the fixed-price carbon trading mechanism, the stepped carbon trading mechanism is flexible and can be adjusted according to market demand and policy goals to improve the economic benefits of carbon emission reduction. The stepped carbon trading mechanism model is as follows: ; In the formula, is the base price of carbon trading; is the interval length; is the price growth rate; is the total cost of stepped carbon trading.

[0045] The carbon tax model is as follows: ; Wherein, is the total carbon tax cost; is the carbon tax coefficient.

[0046] The green certificate trading mechanism model is as follows: The full name of the green certificate is the renewable energy green power certificate, which is used for the accounting of renewable energy power consumption and the certification of renewable energy power consumption. One unit of green certificate corresponds to 1000 kWh of renewable energy electricity. The specific model of the green certificate trading mechanism is as follows: ; Wherein, is the green certificate quota index, is the actual number of green certificates obtained by the system, is the green certificate quota coefficient, is the wind power consumed by the system during the period, is the green certificate income coefficient, is the income of the green certificate trading mechanism.

[0047] In this embodiment, in step S4, according to the set objective function, a low-carbon economic dispatch model and constraint conditions are constructed; Specifically, the objective function of the low-carbon economic dispatch model is: ; Wherein, is the total system operation cost; is the main network power purchase cost; is the equipment operation and maintenance cost; is the stepped carbon trading cost, is the carbon tax cost, is the green certificate income, is the curtailment penalty cost, is the additional income of zero-carbon energy; Among them, the energy purchase cost is as follows: ; Wherein, is the power purchase price, is the gas purchase price, is the volume of natural gas purchased during the period.

[0048] The equipment operation and maintenance cost is as follows: ; Wherein, and are the operation and maintenance costs per unit power of a gas turbine, a gas boiler, a hydrogen fuel cell, an electric boiler, and a methane cracking device, respectively, is the price per kilogram of water, is the price per kilogram of carbon.

[0049] Additional revenue from zero-carbon energy is as follows: ; In the formula, and are the additional price coefficients of zero-carbon electric energy and zero-carbon thermal energy, respectively.

[0050] Penalty cost for wind curtailment is as follows: ; In the formula, is the predicted wind power, is the wind power actually used by the system, is the unit wind curtailment penalty coefficient.

[0051] In this embodiment, the constraint conditions include: wind power generation constraint, coupling constraint, power balance constraint, and energy storage constraint.

[0052] Among them, the electric power balance is as follows: ; In the formula, is the wind power used in the is time period, and

[0053] the electric load of the system in the ; In the formula, is the thermal load of the system in the

[0054] The gas power balance is as follows: ; In the formula, is the volume of natural gas purchased by the system in the is the volume of natural gas consumed for producing turquoise hydrogen in the is the volume of natural gas consumed by the gas turbine in the is the volume of natural gas consumed by the gas boiler in the

[0055] The hydrogen power balance is as follows: ; In the formula, for The volume of hydrogen produced by the electrolyzer during the period; for The volume of hydrogen produced by the methane cracking equipment during the period; for The volume of hydrogen released by the hydrogen storage device during the period; for The volume of hydrogen consumed by the hydrogen fuel cell during the period; for The volume of hydrogen consumed by the gas turbine during the period; for The volume of hydrogen consumed by the gas boiler during the period; for The volume of hydrogen injected into the hydrogen storage device during a period of time; for The hydrogen load of the system during the period.

[0056] In this embodiment, the low-carbon economic dispatch model is solved using MATLAB software and CPLEX solver.

[0057] In this embodiment, in step S5, the set parameter data is input into the low-carbon economic dispatch model; the low-carbon economic dispatch model is used to process the data and output a dispatch strategy for the hydrogen-containing integrated energy system.

[0058] Specifically, the forecast data of wind power and load; the correlation coefficients and operation and maintenance costs of energy storage equipment, gas turbines and other equipment are input into the low-carbon economic dispatch model; the low-carbon economic dispatch model is used to process the data to obtain the optimal output curve and operation strategy of each equipment and unit in the system, and the changes in the total operating cost and carbon emissions of the system caused by the additional price of zero-carbon energy are analyzed.

[0059] In a possible embodiment, a verification example is provided as follows: In order to verify the impact of the coordination of hydrogen production methods and users' low-carbon willingness on the system's economy and low-carbon performance, three schemes were set up to verify the effectiveness of the proposed model, and the impact of different zero-carbon energy additional price coefficients on the system's total operating cost and carbon emissions was analyzed.

[0060] Option 1: Optimal scheduling of the electricity-heat-hydrogen integrated energy system without considering the equipment for turquoise hydrogen production; Option 2: Optimal scheduling of the electricity-heat-hydrogen integrated energy system considering turquoise hydrogen production; Option 3: Optimal scheduling of the electricity-heat-hydrogen integrated energy system considering turquoise hydrogen production and users' low-carbon willingness.

[0061] The scheduling results of the three scenarios are compared as shown in Table 1: Table 1 Comparison of Scheduling Results of Scenarios 1 - 3

[0062] It can be seen from Table 1 that after introducing the turquoise hydrogen production equipment in Scenario 2, compared with Scenario 1, the total operating cost of the system is reduced by 1496.89 yuan, the energy purchase cost is reduced by 400.07 yuan, and the carbon constraint cost is reduced by 865.12 yuan. From Figure 7 and Figure 11 it can be seen that this is because in Scenario 1, more electric energy is consumed in green hydrogen production. During the period from 7 to 21, wind power is difficult to meet the electrical load demand of the system, and natural gas can only be purchased and electricity is generated through a gas turbine to maintain the stable operation of the system. After introducing the turquoise hydrogen production equipment in Scenario 2, during the period from 0 to 4, considering that the wind power is relatively abundant, at this time the system mainly produces green hydrogen to meet the hydrogen load demand of the system. From 5 to 21, considering that the wind power is difficult to meet the system demand at this time, the turquoise hydrogen production equipment starts to operate. Compared with green hydrogen, the electric energy consumed per unit mass of turquoise hydrogen produced by the system is significantly less than that of producing green hydrogen. Therefore, during the period from 10 to 21, the power generation of the gas turbine in Scenario 2 is significantly less than that in Scenario 1. However, since the production of turquoise hydrogen requires the consumption of natural gas, from Figure 9 and Figure 13 it can be seen that during this period, the total amount of natural gas consumed in Scenario 2 is not significantly less than that in Scenario 1, but the natural gas consumed by the gas boiler and gas turbine in Scenario 2 from 10 to 12 is significantly higher than that in Scenario 1. In addition, during the period from 10 to 12 in Scenario 2, the power of the electric boiler is significantly higher than that in Scenario 1, which is also because less electric energy is consumed in producing turquoise hydrogen, so the system can convert more electric energy into heat energy through the electric boiler and supply it to users. Comparing Figure 8 and Figure 12 it can be seen that after introducing the turquoise hydrogen production equipment, the thermal energy stored in the system during the periods from 1 to 4 and from 17 to 20 is significantly reduced. This is because the turquoise hydrogen equipment consumes less electric energy in producing hydrogen, effectively reducing the power supply pressure of the system at the peak of the electrical load. The system reduces the output of the gas turbine and increases the output of the gas boiler. Comparing Figure 10 and Figure 14 it can be found that less hydrogen is stored in Scenario 2 than in Scenario 1 at 1 o'clock. This is because after considering the production of turquoise hydrogen in Scenario 2, the problem of more electric energy consumption in hydrogen production at the peak of the electrical load is effectively reduced. Therefore, at 1 o'clock, more wind energy can be converted into heat energy to meet the thermal load demand of the system while reducing the carbon emissions generated by the gas boiler for heating. Through the comparison of Scenario 1 and Scenario 2, it can be seen that introducing the turquoise hydrogen production equipment into the integrated energy system can effectively reduce the system carbon emissions and total operating cost while meeting the hydrogen load demand.

[0063] As can be seen from Table 1, compared with Scheme 2, after considering the low-carbon willingness of users in Scheme 3, the total operating cost of the system decreased by 3,366.07 yuan, and the total carbon emissions decreased by 6.91 kg. Comparing Figure 11 and Figure 15 it can be seen that when the system is at the peak of the electrical load, the use of gas equipment is slightly reduced, the output of the energy storage equipment at the peak of the load is increased, and the proportion of zero-carbon electric energy supplied to the user side by the system during the peak load period is reduced. Comparing Figure 12 and Figure 16 it can be seen that the system mainly adjusts the output time period of the heat storage equipment. During the 1-3 time periods, considering that during this time period, the zero-carbon heat energy supplied by the system to the load side exceeds the user's demand for zero-carbon heat energy, and part of the zero-carbon heat energy cannot obtain benefits. Therefore, it is considered to charge the heat storage equipment during this time period, and at the same time increase the output of the gas boiler to meet the heat load demand of the system. During the 13-14, 16, and 18-19 time periods, at this time, the zero-carbon heat energy supplied by the system to the load side has not exceeded the total amount of zero-carbon heat energy that users are willing to purchase. Therefore, the system reduces the heat storage power of the heat storage equipment during the 13-14 and 18 time periods, reduces the heat release power of the heat storage equipment during the 17 time period, and increases the heat release power of the heat storage equipment during the 18-19 time period, making full use of the user's willingness to purchase zero-carbon energy and further reducing the operating cost of the system. Comparing Figure 13 and Figure 17 it can be seen that in order to use the user's willingness to purchase zero-carbon energy to reduce the operating cost of the system, the system slightly adjusts the natural gas purchased from the main grid in each time period, so as to adjust the proportion of zero-carbon energy in the energy supplied to the user side in each stage, so as to maximize the zero-carbon energy income. Comparing Figure 14 and Figure 18 it can be seen that the system reduces the input power of the energy storage equipment in the 1 time period, increases the hydrogen blending ratio of the gas boiler during the 1-3 time periods, and increases the proportion of zero-carbon heat energy supplied by the system to the user side during the 1-3 time periods, so as to obtain more zero-carbon energy income.

[0064] The relationship between the additional price coefficient of zero-carbon energy and the total operating cost and carbon emissions of the system is as Figure 19 and 20 shown. From Figure 19It can be seen that the minimum total operating cost of the system is 12,092.96 yuan. When the additional price coefficient of zero-carbon heat energy is fixed, as the additional price coefficient of zero-carbon electrical energy increases, the total operating cost of the system shows a trend of first decreasing and then increasing. Without considering the impact of the additional price of zero-carbon heat energy, when the additional price coefficient of zero-carbon electrical energy is in the range of 0.25 - 0.375, the total operating cost of the system is the lowest. This is because the additional price of zero-carbon energy is related to the price of ordinary energy. The system adopts time-of-use electricity prices. During the peak electricity load, the price of ordinary electrical energy is significantly higher than that during the valley electricity load. The profit obtained from selling zero-carbon electrical energy with the same capacity during the peak electricity load is significantly higher than that during the valley electricity load. Therefore, when the additional price coefficient of zero-carbon electrical energy starts to increase, although the supply of zero-carbon electrical energy at night exceeds the zero-carbon electrical energy that users are willing to pay an additional price for, the system loses a part of the profit that can be received from zero-carbon electrical energy at night. However, during the peak electricity load, the profit obtained from the zero-carbon electrical energy supplied by the system increases. After comprehensive calculation, the total profit obtained by the system from zero-carbon electrical energy is higher. When the additional price coefficient of zero-carbon electrical energy exceeds 0.375, the total operating cost of the system shows an upward trend. This is because as the additional price coefficient increases, the number of users willing to pay an additional fee to purchase zero-carbon electrical energy during the peak electricity load gradually decreases. When the additional price coefficient of zero-carbon electrical energy exceeds 0.375, the zero-carbon electrical energy provided by the system during the peak electricity load exceeds the capacity of zero-carbon electrical energy that users are willing to purchase. At this time, some zero-carbon electrical energy cannot obtain benefits, resulting in an increase in the operating cost of the system.

[0065] When the additional price coefficient of zero-carbon electrical energy in the system is fixed, the total operating cost of the system increases as the additional price coefficient of zero-carbon heat energy rises. This is because the wind power resources are relatively abundant at night, and zero-carbon electrical energy can be converted into zero-carbon heat energy through an electric boiler and supplied to users. Moreover, the heat load on the user side is relatively high at night. The increase in the additional price coefficient of zero-carbon heat energy will lead to a decrease in the capacity of zero-carbon heat energy that users are willing to purchase, resulting in a decrease in the income from zero-carbon heat energy provided at night and thus a decrease in the total operating cost of the system. In addition, since the heat load is relatively low during the day and the heat load is mainly supplied through gas equipment, the amount of zero-carbon heat energy that can be provided is relatively small. Therefore, the change in the additional price coefficient of zero-carbon heat energy has little impact on the profit obtained from supplying zero-carbon heat energy by the system during the day.

[0066] For Figure 20It can be seen that the minimum carbon emission of the system is 1586.03 kg. The carbon emission of the system is less affected by the change of the additional price coefficient of zero-carbon energy. Only when the additional price coefficient of zero-carbon electricity is in the range of 0.1 - 0.375, the carbon emission of the system changes significantly. When the additional price coefficient of zero-carbon electricity of the system is 0.15 and the additional price coefficient of zero-carbon heat is 0.2, the carbon emission of the system reaches the minimum. When the additional price coefficient of zero-carbon electricity is 0.15 and the additional price coefficient of zero-carbon heat is 0.45, the carbon emission of the system reaches the maximum. This is because when the additional price coefficient of zero-carbon electricity increases from 0 to 0.25, the system converts part of the zero-carbon electricity into zero-carbon heat, reduces the thermal power of the gas boiler at night, and reduces the carbon emission of the system. With the further increase of the additional price coefficient of zero-carbon electricity, in order to obtain the profit of zero-carbon electricity at the peak of the electric load, the system increases the output of gas equipment during the valley and normal periods of the electric load, increases the zero-carbon energy stored in the energy storage equipment in this way, and then the energy storage equipment releases it at the peak of the electric load to obtain the profit of zero-carbon energy. However, there is energy loss in the energy storage equipment, and the consumption of natural gas by gas equipment increases, resulting in an increase in the carbon emission of the system. The same is true for the increase in the carbon emission of the system with the increase of the additional price coefficient of zero-carbon heat.

[0067] In summary, according to the hydrogen production method by methane cracking for producing turquoise hydrogen, the present invention constructs a mathematical model for producing turquoise hydrogen; by introducing a coordinated hydrogen production method considering green hydrogen and turquoise hydrogen, a hydrogen-containing integrated energy system model is constructed; according to the user's willingness to purchase zero-carbon energy, by analyzing the carbon footprint of the energy supplied by the hydrogen-containing integrated energy system during a set period, the energy supplied by the hydrogen-containing integrated energy system to the load side is divided into zero-carbon electric energy, zero-carbon thermal energy, ordinary electric energy, and ordinary thermal energy; by setting an additional price coefficient for the zero-carbon energy, the additional income of the zero-carbon energy is obtained; according to the set objective function, a low-carbon economic dispatch model and constraint conditions are constructed; the set parameter data is input into the low-carbon economic dispatch model; through the processing of the low-carbon economic dispatch model, the dispatch strategy of the hydrogen-containing integrated energy system is output. The present invention takes into account that the power consumption for producing green hydrogen in the system is relatively high. When the system's electrical load reaches its peak, in order to meet the hydrogen load demand, the output of the gas turbine unit is increased, resulting in a significant increase in the system's carbon emissions. By introducing a turquoise hydrogen production device, the production volume of turquoise hydrogen during the electrical load peak is increased, reducing the power supply pressure of the system during the electrical load peak, and realizing the economic, environmental, safe, and reliable operation of the system. The physical meaning is clear, the method is scientific and reasonable, has strong applicability, and the dispatch strategy is reasonable and effective. The present invention introduces turquoise hydrogen into the integrated energy system, constructs a mathematical model for producing turquoise hydrogen by a methane cracking device, and by adjusting the output of the green hydrogen production device and the turquoise hydrogen production device at different times, on the basis of meeting the hydrogen load, reduces the electrical energy consumed for producing hydrogen during the electrical load peak of the system, reduces the carbon emissions of the system, and at the same time differentiates the energy supplied by the system to the load side into zero-carbon energy and ordinary energy, considering the influence of the user's low-carbon willingness, and reduces the total operating cost of the system by adjusting the price difference between zero-carbon energy and ordinary energy.

[0068] It should be noted that the method of the embodiments of the present disclosure can be executed by a single device, such as a computer or a server. The method of this embodiment can also be applied to a distributed scenario, and completed by multiple devices cooperating with each other. In this case of a distributed scenario, one of the multiple devices can only execute one or more steps of the method of the embodiments of the present disclosure, and these multiple devices will interact with each other to complete the described method.

[0069] It should be noted that some embodiments of the present disclosure have been described. In some cases, the recorded actions or steps can be executed in a different order from that in the above embodiments and still achieve the desired results. Additionally, the processes depicted in the drawings do not necessarily require the specific order or continuous order shown to achieve the desired results. In certain embodiments, multitasking and parallel processing are also possible or may be advantageous.

[0070] Embodiment 2 See Figure 21, Embodiment 2 of the present invention further provides an IES scheduling device based on coordinating hydrogen production and low-carbon willingness, including: A mathematical model construction module 001 for producing turquoise hydrogen, which is used to construct a mathematical model for producing turquoise hydrogen according to the methane cracking hydrogen production method for producing turquoise hydrogen; A hydrogen-containing integrated energy system model construction module 002, which is used to construct a hydrogen-containing integrated energy system model by introducing a coordinated hydrogen production method considering green hydrogen and turquoise hydrogen; An energy division module 003, which is used to divide the energy supplied by the hydrogen-containing integrated energy system to the load side into zero-carbon electric energy, zero-carbon heat energy, ordinary electric energy, and ordinary heat energy according to the user's willingness to purchase zero-carbon energy by analyzing the carbon footprint of the energy supplied by the hydrogen-containing integrated energy system during a set period; by setting an additional price coefficient for the zero-carbon energy, obtaining the additional income of the zero-carbon energy; A low-carbon economic scheduling model and constraint condition construction module 004, which is used to construct a low-carbon economic scheduling model and constraint conditions according to a set objective function; A scheduling strategy output module 005, which is used to input set parameter data into the low-carbon economic scheduling model; through processing by the low-carbon economic scheduling model, outputting a scheduling strategy for the hydrogen-containing integrated energy system.

[0071] In this embodiment, in the mathematical model construction module 001 for producing turquoise hydrogen, the expression of the mathematical model for producing turquoise hydrogen is: ; In the formula, is the volume of turquoise hydrogen produced by methane cracking in the time period; is the volume of natural gas consumed by the turquoise hydrogen production equipment in the time period; is the electric energy consumed by the turquoise hydrogen production equipment per unit volume of turquoise hydrogen produced in the time period; is the mass of carbon produced by the methane cracking reaction in the time period; is the efficiency of methane conversion to carbon; is the efficiency of methane conversion to hydrogen; is the lower heating value of carbon; is the lower heating value of methane.

[0072] In this embodiment, in the hydrogen-integrated energy system model construction module 002, in the hydrogen-integrated energy system model, the electrical load of the hydrogen-integrated energy system is supplied by purchasing electricity from the power grid, hydrogen fuel cells, wind turbines, and gas turbines; The thermal load of the hydrogen-integrated energy system is supplied by gas turbines, hydrogen fuel cells, electric boilers, and gas boilers; The hydrogen load of the hydrogen-integrated energy system is supplied by electrolytic hydrogen production equipment and methane cracking equipment; The energy storage equipment of the hydrogen-integrated energy system includes electricity storage equipment, heat storage equipment, and hydrogen storage equipment.

[0073] In this embodiment, in the energy division module 003, the mathematical expression for the user's willingness to purchase zero-carbon energy is: ; ; ; In the formula, is the user who is very willing; is the user who is relatively willing; is the neutral user; is the user who is relatively unwilling; is the user who is completely unwilling; is the user willingness distribution matrix after the k-th iteration; is the willingness of a single user to purchase zero-carbon energy.

[0074] In this embodiment, in the energy division module 003, the mathematical model expression for dividing the energy supplied by the hydrogen-integrated energy system to the load side into the zero-carbon electric energy, the zero-carbon thermal energy, the ordinary electric energy, and the ordinary thermal energy is: ; In the formula, and are respectively the system electrical load and thermal load during the and are respectively the zero-carbon electric energy and zero-carbon thermal energy produced by the system during the is the volume of natural gas consumed by the gas turbine during the is the volume of natural gas consumed by the gas boiler during the is the main grid power purchase volume during the is the electrical conversion efficiency of the gas turbine; is the thermal conversion efficiency of the gas turbine, is the thermal conversion efficiency of the gas boiler.

[0075] In this embodiment, in the low-carbon economic dispatch model and constraint construction module 004, the objective function of the low-carbon economic dispatch model is as follows: ; In the formula, is the total system operation cost; is the main network power purchase cost; is the equipment operation and maintenance cost; is the stepped carbon trading cost, is the carbon tax cost, is the green certificate income, is the wind curtailment penalty cost, is the additional income of zero-carbon energy; The constraints include: wind power generation constraint, coupling constraint, power balance constraint and energy storage constraint.

[0076] It should be noted that the information interaction, execution process, etc. between the above system modules, due to being based on the same concept as the method embodiment in Embodiment 1 of this application, have the same technical effects as the method embodiment of this application. The specific content can be referred to the description in the method embodiment shown above in this application and will not be elaborated here.

[0077] Embodiment 3 Embodiment 3 of the present invention provides a non-transitory computer-readable storage medium, in which program codes of an IES dispatch method based on coordinated hydrogen production and low-carbon willingness are stored, and the program codes include instructions for executing the IES dispatch method based on coordinated hydrogen production and low-carbon willingness in Embodiment 1 or any possible implementation manner thereof.

[0078] The computer-readable storage medium can be any available medium that can be accessed by a computer or a data storage device such as a server or a data center integrating one or more available media. The available media can be magnetic media (for example, floppy disks, hard disks, magnetic tapes), optical media (for example, DVDs), or semiconductor media (for example, solid state disks (Solid State Disk, SSD)), etc.

[0079] Embodiment 4 Embodiment 4 of the present invention provides an electronic device, including: a memory and a processor; The processor and the memory communicate with each other through a bus; the memory stores program instructions executable by the processor, and the processor can execute the IES dispatch method based on coordinated hydrogen production and low-carbon willingness in Embodiment 1 or any possible implementation manner thereof by calling the program instructions.

[0080] Specifically, the processor can be implemented by hardware or software. When implemented by hardware, the processor can be a logic circuit, an integrated circuit, etc.; when implemented by software, the processor can be a general-purpose processor that realizes its functions by reading software codes stored in a memory. The memory can be integrated in the processor or exist independently outside the processor.

[0081] In the above embodiments, they can be implemented in whole or in part by software, hardware, firmware, or any combination thereof. When implemented using software, they can be implemented in whole or in part in the form of a computer program product. The computer program product includes one or more computer instructions. When the computer program instructions are loaded and executed on a computer, the processes or functions described in the embodiments of the present invention are generated in whole or in part. The computer can be a general-purpose computer, a special-purpose computer, a computer network, or other programmable systems. The computer instructions can be stored in a computer-readable storage medium or transmitted from one computer-readable storage medium to another. For example, the computer instructions can be transmitted from a website, a computer, a server, or a data center to another website, a computer, a server, or a data center by wire (such as coaxial cable, optical fiber, digital subscriber line (DSL)) or wirelessly (such as infrared, wireless, microwave, etc.).

[0082] Obviously, those skilled in the art should understand that the above-mentioned modules or steps of the present invention can be implemented by a general computing system. They can be concentrated on a single computing system or distributed on a network composed of multiple computing systems. Optionally, they can be implemented by program codes executable by the computing system. Thus, they can be stored in a storage system and executed by the computing system. And in some cases, the steps shown or described can be executed in a sequence different from that here, or they can be separately fabricated into individual integrated circuit modules, or multiple modules or steps among them can be fabricated into a single integrated circuit module for implementation. In this way, the present invention is not limited to any specific combination of hardware and software.

[0083] Although the present invention has been described in detail above with general descriptions and specific embodiments, modifications or improvements can be made to it on the basis of the present invention, which are obvious to those skilled in the art. Therefore, these modifications or improvements made without departing from the spirit of the present invention all fall within the scope of protection required by the present invention.

Claims

1. An IES scheduling method based on coordinating hydrogen production and low-carbon willingness, characterized in that Including: According to the hydrogen production method by methane pyrolysis for producing turquoise hydrogen, a mathematical model for producing turquoise hydrogen is constructed; By introducing a coordinated hydrogen production method considering green hydrogen and turquoise hydrogen, a hydrogen-integrated energy system model is constructed; According to the user's willingness to purchase zero-carbon energy, by analyzing the carbon footprint of the energy supplied by the hydrogen-integrated energy system during a set period, the energy supplied by the hydrogen-integrated energy system to the load side is divided into zero-carbon electric energy, zero-carbon thermal energy, ordinary electric energy, and ordinary thermal energy; by setting an additional price coefficient for the zero-carbon energy, the additional income of the zero-carbon energy is obtained; According to the set objective function, a low-carbon economic dispatch model and constraint conditions are constructed; The set parameter data is input into the low-carbon economic dispatch model; through the processing of the low-carbon economic dispatch model, the dispatch strategy of the hydrogen-integrated energy system is output.

2. The IES scheduling method based on coordinating hydrogen production and low-carbon willingness according to claim 1, wherein The expression of the mathematical model for producing turquoise hydrogen is: ； In the formula, is the volume of turquoise hydrogen produced by methane pyrolysis during the is the volume of natural gas consumed by the turquoise hydrogen production equipment during the is the electric energy consumed by the turquoise hydrogen production equipment during the is the electric energy consumed to produce a unit volume of turquoise hydrogen during the is the mass of carbon produced by the methane pyrolysis reaction during the is the efficiency of methane conversion to carbon; is the efficiency of methane conversion to hydrogen; and are the upper and lower limits of hydrogen production by the methane pyrolysis equipment at any time; is the low calorific value of carbon; is the low calorific value of methane.

3. The IES scheduling method based on coordinating hydrogen production and low-carbon willingness according to claim 2, characterized in that, In the hydrogen-integrated energy system model, the electrical load of the hydrogen-integrated energy system is supplied by purchasing electricity from the grid, hydrogen fuel cells, wind turbines, and gas turbines; The thermal load of the hydrogen-integrated energy system is supplied by gas turbines, hydrogen fuel cells, electric boilers, and gas boilers; The hydrogen load of the hydrogen-integrated energy system is supplied by electrolytic hydrogen production equipment and methane pyrolysis equipment; The energy storage equipment of the hydrogen-integrated energy system includes electricity storage equipment, heat storage equipment, and hydrogen storage equipment.

4. The IES scheduling method based on coordinating hydrogen production and low-carbon willingness according to claim 3, wherein The mathematical expression of the user's willingness to purchase zero-carbon energy is: ； ； ； wherein, represents users who are very willing; represents users who are relatively willing; represents neutral users; represents users who are relatively unwilling; represents users who are completely unwilling; represents the user willingness distribution matrix after the k-th iteration; represents the purchase willingness of a single user for zero-carbon energy.

5. The IES scheduling method based on coordinating hydrogen production and low-carbon willingness according to claim 4, characterized in that, The mathematical model expression for dividing the energy supplied by the hydrogen-integrated energy system to the load side into the zero-carbon electric energy, the zero-carbon thermal energy, the ordinary electric energy, and the ordinary thermal energy is: ； In the formula, and are respectively the system's electrical load and heat load during a period; and are respectively the zero-carbon electric energy and zero-carbon heat energy produced by the system during a period; is the volume of natural gas consumed by the gas turbine during a period; is the volume of natural gas consumed by the gas boiler during a period; is the main network power purchase during a period; is the electrical conversion efficiency of the gas turbine; is the heat conversion efficiency of the gas turbine, is the heat conversion efficiency of the gas boiler.

6. The IES scheduling method based on coordinating hydrogen production and low-carbon willingness according to claim 5, characterized in that The objective function of the low-carbon economic dispatch model is: ； In the formula, is the total operating cost of the system; is the main online power purchase cost; is the equipment operation and maintenance cost; is the stepped carbon trading cost, is the carbon tax cost, is the green certificate income, is the penalty cost for wind curtailment, is the additional income of zero-carbon energy; The constraint conditions include: wind power generation constraint, coupling constraint, power balance constraint, and energy storage constraint.

7. The IES dispatching device based on coordinating hydrogen production and low-carbon willingness adopts the IES dispatching method based on coordinating hydrogen production and low-carbon willingness according to any one of claims 1-6, characterized in that Including: A module for constructing a mathematical model for producing turquoise hydrogen, which is used to construct a mathematical model for producing turquoise hydrogen according to the hydrogen production method by methane pyrolysis for producing turquoise hydrogen; A module for constructing a hydrogen-integrated energy system model, which is used to construct a hydrogen-integrated energy system model by introducing a coordinated hydrogen production method considering green hydrogen and turquoise hydrogen; An energy division module, which is used to divide the energy supplied by the hydrogen-integrated energy system to the load side into zero-carbon electric energy, zero-carbon thermal energy, ordinary electric energy, and ordinary thermal energy according to the user's willingness to purchase zero-carbon energy by analyzing the carbon footprint of the energy supplied by the hydrogen-integrated energy system during a set period; by setting an additional price coefficient for the zero-carbon energy, the additional income of the zero-carbon energy is obtained; A module for constructing a low-carbon economic dispatch model and constraint conditions, which is used to construct a low-carbon economic dispatch model and constraint conditions according to the set objective function; A dispatch strategy output module, which is used to input the set parameter data into the low-carbon economic dispatch model; through the processing of the low-carbon economic dispatch model, the dispatch strategy of the hydrogen-integrated energy system is output.

8. The IES scheduling device based on coordinating hydrogen production and low-carbon willingness according to claim 7, characterized in that, In the module for constructing a mathematical model for producing turquoise hydrogen, the expression of the mathematical model for producing turquoise hydrogen is: ； Wherein, is the volume of turquoise hydrogen produced by methane cracking during the is the volume of natural gas consumed by the turquoise hydrogen production equipment during the is the electric energy consumed by the turquoise hydrogen production equipment during the is the electric energy consumed to produce a unit volume of turquoise hydrogen during the is the mass of carbon produced by the methane cracking reaction during the is the efficiency of methane conversion to carbon; is the efficiency of methane conversion to hydrogen; and are the upper and lower limits of hydrogen production by the methane cracking equipment at any time; is the low calorific value of carbon; is the low calorific value of methane.

9. The IES scheduling device based on coordinating hydrogen production and low-carbon willingness according to claim 8, characterized in that, In the module for constructing a hydrogen-integrated energy system model, in the hydrogen-integrated energy system model, the electrical load of the hydrogen-integrated energy system is supplied by purchasing electricity from the grid, hydrogen fuel cells, wind turbines, and gas turbines; The thermal load of the hydrogen-integrated energy system is supplied by gas turbines, hydrogen fuel cells, electric boilers, and gas boilers; The hydrogen load of the hydrogen-integrated energy system is supplied by electrolytic hydrogen production equipment and methane cracking equipment; The energy storage equipment of the hydrogen-integrated energy system includes electricity storage equipment, heat storage equipment, and hydrogen storage equipment.

10. The IES scheduling device based on coordinating hydrogen production and low-carbon willingness according to claim 9, characterized in that, In the said energy division module, the mathematical expression for the user's willingness to purchase zero-carbon energy is: ； ； ； In the formula, represents users who are very willing; represents users who are relatively willing; represents neutral users; represents users who are relatively unwilling; represents users who are completely unwilling; represents the user willingness distribution matrix after the k-th iteration; represents the purchase willingness of a single user for zero-carbon energy.

Images