Energy island planning and designing method based on mathematical optimization model

Abstract

The invention relates to an energy island planning and designing method based on a mathematical optimization model, which comprises the following steps: firstly, establishing a topological structure diagram of an energy island system, and then constructing an energy island system mathematical model according to the topological structure diagram of the energy island system by taking the maximum capital internal return rate and the minimum abandoned electric quantity as targets; and finally, inputting the mathematical model of the energy island system into an optimization solver, and calculating and solving to obtain an optimal configuration and operation scheme of the energy island system. According to the method, an energy flow topology model is constructed, an optimization target and a constraint condition are established, and a commercial solver is combined for solving, so that the optimal planning design of the energy island under the multi-energy cooperative condition is realized.

Description

Technical Field

[0001] This invention relates to the field of integrated energy system planning technology, and in particular to an energy island planning and design method based on a mathematical optimization model. Background Technology

[0002] With the acceleration of the global energy transition, the large-scale development and utilization of clean energy has gradually become an important direction for energy structure adjustment. In recent years, energy islands, as a new type of integrated and multi-energy complementary energy system, have gradually attracted widespread attention. Energy islands are usually located in the sea, and by aggregating renewable energy sources such as offshore wind power, photovoltaics, and ocean energy, and combining them with diversified loads such as energy storage, hydrogen production, seawater desalination, and computing centers, they achieve local consumption and efficient utilization of energy.

[0003] However, existing energy island planning and design methods mostly adopt decentralized or experience-driven approaches, which have the following shortcomings: lack of holistic modeling: existing studies often focus on planning for a single energy form such as electricity or hydrogen, failing to fully consider the coupling relationship between multiple energy flows such as electricity, hydrogen, and energy storage; difficulty in balancing multiple objectives: in practical applications, energy islands must meet the economic requirements of energy supply while ensuring the reliability and environmental friendliness of system operation, and existing methods are insufficient in balancing multiple objectives; lack of scalable optimization methods: traditional heuristic or static planning methods are difficult to cope with the complex equipment configuration, energy flow, and dynamic load characteristics of energy islands, and cannot guarantee the global optimality and adaptability of the solution. Summary of the Invention

[0004] Based on this, an energy island planning and design method based on a mathematical optimization model is provided. By constructing an energy flow topology model, establishing optimization objectives and constraints, and combining it with a commercial solver, the optimal planning and design of the energy island under multi-energy synergy conditions can be achieved.

[0005] This invention provides a method for planning and designing energy islands based on a mathematical optimization model, comprising the following steps:

[0006] Establish the topology diagram of the energy island system;

[0007] Based on the topology diagram of the energy island system, a mathematical model of the energy island system is constructed with the objectives of maximizing the internal rate of return on capital and minimizing the amount of abandoned electricity.

[0008] The optimal configuration and operation scheme of the energy island system are obtained by inputting the mathematical model of the energy island system into the optimization solver.

[0009] The mathematical model of the energy island system is as follows:

[0010] maxIRR = f IRR (cash)

[0011]

[0012] st{electric balance constraints, hydrogen balance constraints, ammonia balance constraints, alcohol balance constraints, electricity purchase constraints, electrical energy storage constraints, hydrogen equipment constraints, synthetic ammonia equipment constraints, alcohol production equipment constraints, energy island constraints, hydrogen storage constraints, ammonia storage constraints, and alcohol storage constraints};

[0013] In the formula, IRR is the internal rate of return on capital, cash is the total cash flow of the energy island system, and f IRR The formula for calculating the internal rate of return (IRR) of capital is e. abandon For the abandoned electricity of the energy island system, yPG ki Let xRPV be the electricity purchased in the i-th hour of year k, xRPV be the installed photovoltaic capacity on the island, and effPV be the electricity purchased in the i-th hour of year k. i Let capPV be the photovoltaic power generation coefficient in the i-th hour. k Let yBES be the remaining capacity factor of the wind turbine in year k. ki Let OE be the charge / discharge amount of electrical energy storage in the i-th hour of year k. ki Let lE be the amount of electricity generated by ocean energy in the i-th hour of year k. i Let ySA be the electrical load for the i-th hour. ki Let y_cSA be the ammonia production capacity of the ammonia synthesis unit in the i-th hour of year k, and y_ECU be the unit power consumption of the ammonia synthesis unit. ki Let be the ammonia production output of the hydrogen production unit in the i-th hour of year k, and ecECU be the unit power consumption of the hydrogen production unit, yMMW. ki Let x represent the alcohol production capacity of the alcohol production equipment in the i-th hour of year k, ecMMW represent the unit power consumption of the ammonia production equipment, xCPC represent the installed capacity of the computing center, ecCPC represent the unit power consumption of the computing center, xDAS represent the installed capacity of the seawater desalination equipment, ecDAS represent the unit power consumption of the seawater desalination equipment, xFE represent the installed capacity of the refueling equipment, ecFE represent the unit power consumption of the refueling equipment, xOMH represent the construction scale of the maintenance home port, ecOMH represent the unit power consumption of the maintenance home port, xMRC represent the construction scale of the rescue center, ecMRC represent the unit power consumption of the rescue center, xCAT represent the construction scale of the cultural tourism area, ecCAT represent the unit power consumption of the cultural tourism area, xMCT represent the construction scale of the aquaculture area, ecMCT represent the unit power consumption of the aquaculture area, and eCELD represent the unit power consumption of the energy island.

[0014] In one embodiment of the formula, the total cash flow of the energy island system is calculated as follows:

[0015]

[0016] cash = [cash1, cash2, ..., cash]k ,...,cash K ]

[0017] In the formula, cash k For the cash flow in year k, reward k For the income in year k, Let the investment cost be in year k. The maintenance cost in year k is... Let K be the electricity cost in year k, and K be the investment period of the project.

[0018] In one embodiment of the formula, the investment cost calculation formula for year k is:

[0019]

[0020] In the formula, Let x_k be the investment cost coefficient for the photovoltaic equipment on the island in year k, and x_ECU be the installed capacity of the hydrogen production equipment. Let xSA be the investment cost coefficient for hydrogen production equipment in year k, and xSA be the installed capacity of ammonia synthesis equipment. Let xMMW be the investment cost coefficient for the ammonia synthesis equipment in year k, and xMMW be the installed capacity of the methanol production equipment. Let be the investment cost coefficient for the alcohol production equipment in year k. Let be the investment cost coefficient of the computing center in year k. Let xELD be the investment cost coefficient for the hydrogen production equipment in year k, and xELD be the total area of ​​the energy island. Let xOE be the investment cost coefficient for the energy island in year k, and xOE be the installed capacity of the marine hub. Let xPH be the investment cost coefficient for ocean energy in year k, and xPH be the installed capacity of the power hub. Let be the investment cost coefficient of the power hub in year k. This represents the investment cost coefficient for the equipment in year k. This represents the investment cost coefficient for the maintenance of the home port in year k. Let be the investment cost coefficient of the rescue center in year k. Let be the investment cost coefficient for cultural tourism in year k. Let xSBAT be the investment cost coefficient for aquaculture in year k, and let xSBAT be the installed capacity of electrical energy storage. Let xSH be the investment cost coefficient for electric energy storage in year k, and xSH be the installed capacity of hydrogen storage equipment. Let x be the investment cost coefficient for hydrogen storage in year k, and xHSA be the installed capacity of ammonia storage equipment. Let xHMW be the investment cost coefficient for the ammonia storage equipment in year k, and let xHMW be the installed capacity of the methanol storage equipment. Let be the investment cost coefficient for alcohol storage in year k.

[0021] In one embodiment of the formula, the maintenance cost calculation formula for year k is:

[0022]

[0023] In the formula, Let be the maintenance cost coefficient for the photovoltaic equipment on the island in year k. Let be the maintenance cost coefficient for the hydrogen production equipment in year k. Let be the maintenance cost coefficient for the ammonia synthesis equipment in year k. Let be the maintenance cost coefficient for the alcohol production equipment in year k. Let be the maintenance cost coefficient of the computing center in year k. Let be the maintenance cost coefficient for the hydrogen production equipment in year k. Let be the maintenance cost coefficient of the energy island in year k. Let be the maintenance cost coefficient for ocean energy in year k. Let be the maintenance cost coefficient of the power hub in year k. This represents the maintenance cost coefficient for the refueling equipment in year k. This represents the maintenance cost coefficient for the home port in year k. Let this be the maintenance cost coefficient of the rescue center in year k. Let be the maintenance cost coefficient for cultural tourism in year k. Let this be the maintenance cost coefficient for aquaculture in year k. Let be the maintenance cost coefficient for energy storage in year k. Let be the maintenance cost coefficient for hydrogen storage in year k. Let be the maintenance cost coefficient for the ammonia storage equipment in year k. Let be the maintenance cost coefficient for alcohol storage in year k.

[0024] In one embodiment of the formula, the electricity cost calculation formula for year k is:

[0025]

[0026] In the formula, pe i The electricity cost for the i-th hour;

[0027] The formula for calculating the return in year k is:

[0028]

[0029] In the formula, se i sh i ssa i smmwi Let LE be the selling price of electricity, hydrogen, ammonia, and alcohol in the i-th hour. i For the electrical load in the i-th hour, lH i For the hydrogen load in the i-th hour, lSA i For the ammonia load in the i-th hour, lMME i For the alcohol load in the i-th hour, Let the unit revenue of the computing center be at the k-th position. Let $\frac{ ... Let the unit revenue of the refueling equipment be in year k. For the unit revenue of the maintenance home port in year k, For the unit revenue of the rescue center in year k, Let the unit revenue of cultural tourism in year k be... Let $\frac{ ...

[0030] In one embodiment of the formula, the electrical balance constraint is:

[0031]

[0032] The hydrogen balance constraint is:

[0033] yECU ki +ySH ki ≥

[0034] lH i +ySA ki ×hcSA+yMMW ki ×hcMMW

[0035] In the formula, ySH ki Let hcSA be the amount of hydrogen stored and released by the hydrogen storage equipment in the i-th hour of year k, hcMA be the unit hydrogen consumption of the ammonia synthesis equipment, and hcMMW be the unit hydrogen consumption of the alcohol production equipment.

[0036] The ammonia balance constraint is:

[0037] ySA ki +yHSA ki ≥lSA i

[0038] In the formula, yHSA ki Let be the amount of ammonia stored and released by the ammonia storage equipment in the i-th hour of year k.

[0039] The equilibrium constraint for alcohols is:

[0040] yMMW ki +yHMW ki ≥lMMW i

[0041] In the formula, yHMW ki Let be the amount of ammonia stored and released by the methanol storage device in the i-th hour of year k.

[0042] The power purchase limit is:

[0043] yPG ki ≤xSPV×effSPV i ×capSPV k +xSWT×effSWT i

[0044] In the formula, xSPV represents the installed capacity of offshore photovoltaic power, and effSPV represents the installed capacity of offshore photovoltaic power. i Let capSPV be the power generation coefficient of offshore photovoltaic power generation in the i-th hour. k Let xSWT be the remaining capacity factor of offshore photovoltaic power in year k, xSWT be the installed capacity of offshore wind turbines, and effSWT be the remaining capacity factor. i Let be the power generation coefficient of offshore photovoltaic power in the i-th hour;

[0045] The constraints of electrical energy storage are:

[0046] yBES ki ≤zBES ki

[0047] -yBES ki ≤xBES-zBES ki

[0048] zBES k,i+1 =zBES ki -yBES ki

[0049]

[0050] In the formula, zBES ki Let xBES be the remaining energy storage capacity at time i in year k, and zBES be the installed capacity of energy storage. k,i+1 Let be the remaining electrical energy stored at time i+1 in year k. The charging rate for electrical energy storage. The discharge rate of the stored energy;

[0051] Hydrogen equipment constraints are: yECU ki ≤xECU;

[0052] The constraints for ammonia synthesis equipment are: ySA ki ≤xSA;

[0053] The constraints for the alcohol production equipment are: yMMW ki ≤xMMW;

[0054] In one embodiment of the formula, the energy island constraint is:

[0055] xELD=xRPV×aRPV+xOE×aOE+xECU×aECU+

[0056] xSA×aSA+xMMW×aMMW+xSH×aSH+

[0057] xHSA×aHSA+xHMW×aHMW+xCPC×aCPC+

[0058] xDAS×aDAS+xPH×aPH+xFE×aFE+

[0059] xOMH×aOMH+xMRC×aMRC+xCAT×aCAT+

[0060] xMCT×aMCT

[0061] In the formula, aRPV represents the area occupied by the unit installed capacity of onshore photovoltaics, aOE represents the area occupied by the unit installed capacity of ocean energy, aECU represents the area occupied by the unit installed capacity of hydrogen production equipment, aSA represents the area occupied by the unit installed capacity of ammonia synthesis equipment, aMMW represents the area occupied by the unit installed capacity of methanol production equipment, aSH represents the area occupied by the unit installed capacity of hydrogen storage equipment, aHSA represents the area occupied by the unit installed capacity of ammonia storage equipment, aHMW represents the area occupied by the unit installed capacity of methanol storage, aCPC represents the area occupied by the unit installed capacity of computing center, aDAS represents the area occupied by the unit installed capacity of seawater desalination, aPH represents the area occupied by the unit installed capacity of power hub, aFE represents the area occupied by the unit installed capacity of refueling equipment, aOMH represents the area occupied by the unit construction capacity of operation and maintenance home port, aMRC represents the area occupied by the unit construction capacity of rescue center, aCAT represents the area occupied by the unit construction capacity of cultural tourism, and aMCT represents the area occupied by the unit construction capacity of aquaculture.

[0062] In one embodiment of the formula, the hydrogen storage constraint is:

[0063] ySH ki ≤zSH ki

[0064] -ySH ki ≤xSH-zSH ki

[0065] zSH k,i+1 =zSH ki -ySH ki

[0066]

[0067] In the formula, zSH kiLet zSH be the remaining hydrogen quantity of the hydrogen storage device at time i in year k. k,i+1 Let be the remaining amount of hydrogen in the hydrogen storage device at time i+1 in year k. For hydrogen storage ratio, This refers to the hydrogen release rate;

[0068] Ammonia storage constraints are:

[0069] yHSA ki ≤zHSA ki

[0070] -yHSA ki ≤xHSA-zHSA ki

[0071] zHSA k,i+1 =zHSA ki -yHSA ki

[0072]

[0073] In the formula, zHSA ki Let zHSA be the amount of remaining ammonia in the ammonia storage device at time i in year k. k,i+1 Let be the amount of ammonia remaining in the ammonia storage device at time i+1 in year k. This refers to the ammonia storage ratio. This refers to the ammonia release ratio;

[0074] The alcohol storage constraints are:

[0075] yHMW ki ≤zHMW ki

[0076] -yHMW ki ≤xHMW-zHMW ki

[0077] zHMW k,i+1 =zHMW ki -yHMW ki

[0078]

[0079] In the formula, zHMW ki Let zHMW be the amount of alcohol remaining in the alcohol storage device at time i in year k. k,i+1 Let be the amount of alcohol remaining in the alcohol storage device at time i+1 in year k. This refers to the ammonia storage ratio. This represents the alcohol release ratio.

[0080] In one embodiment of the formula, the optimal configuration and operation scheme of the energy island system includes:

[0081] The capacity configuration of various production, storage, and consumption devices, as well as the optimal construction scale of various infrastructures in the topology;

[0082] Energy island energy flow and operation scheduling strategies, including the production, storage and consumption of electricity and hydrogen energy at different time scales;

[0083] The system's economic and environmental indicators include total life cycle cost, return on equity, carbon emission levels, and renewable energy utilization rate.

[0084] In one embodiment of the formula, the optimization solver is a linear programming solver, a mixed integer programming solver, or a nonlinear programming solver.

[0085] The beneficial effects of this invention are:

[0086] This invention is the first to apply mathematical optimization methods to the planning and design of energy islands. By establishing a unified coupled model for electricity, hydrogen energy, and energy storage, it achieves multi-energy collaborative modeling and optimization, overcoming the shortcomings of traditional methods that rely on experience-based configuration and single-energy analysis. This method can achieve a trade-off between system economy, energy supply reliability, and environmental friendliness, avoiding resource waste and improving energy utilization efficiency. Simultaneously, the use of commercial solvers for model solving ensures global optimality and computational efficiency. The proposed method has good versatility and scalability, applicable to the construction and operation of energy islands of different scales and in various application scenarios, providing scientific decision support for the future development of integrated marine energy systems. Attached Figure Description

[0087] Figure 1 This is one of the flowcharts illustrating the energy island planning and design method based on a mathematical optimization model provided in an embodiment of the present invention.

[0088] Figure 2 This is a topology diagram of an energy island system provided in an embodiment of the present invention. Detailed Implementation

[0089] To make the objectives, technical solutions, and advantages of this invention clearer, the invention will be further described in detail below with reference to the accompanying drawings and embodiments. It should be understood that the specific embodiments described herein are merely illustrative and not intended to limit the invention.

[0090] In one embodiment, such as Figure 1 As shown, Figure 1 This is one of the flowcharts illustrating the energy island planning and design method based on a mathematical optimization model provided in this invention. When applied to computer equipment, this method includes the following steps:

[0091] S101. Establish the topology diagram of the energy island system. Figure 2 This is a topology diagram of an energy island system provided in an embodiment of the present invention. The topology diagram of the energy island system includes a computing center, seawater desalination, power hub, refueling equipment, energy storage, ocean energy, onshore photovoltaic, offshore photovoltaic, offshore wind power, operation and maintenance home port, rescue center, cultural tourism, aquaculture, basic hydrogen production, hydrogen storage, hydrogen production, ammonia storage, ammonia production, alcohol storage, and alcohol production.

[0092] S102. Based on the topology diagram of the energy island system, construct a mathematical model of the energy island system with the objectives of maximizing the internal rate of return on capital and minimizing the amount of abandoned electricity.

[0093] S103. Input the mathematical model of the energy island system into the optimization solver to obtain the optimal configuration and operation scheme of the energy island system. The optimal configuration and operation scheme of the energy island system includes: the capacity configuration of various production, storage, and consumption devices, as well as the optimal construction scale of various infrastructures in the topology; the energy flow and operation scheduling strategy of the energy island, including the production, storage, and consumption of electricity and hydrogen energy at different time scales; and the economic and environmental indicators of the system, including the total life cycle cost, return on equity, carbon emission level, and renewable energy utilization rate. This scheme can be directly used for the design and construction of energy island projects, and at the same time provides a reference for scheduling optimization in the operation phase, realizing a comprehensive improvement in the economic efficiency and energy supply reliability of the energy island.

[0094] Specifically, the optimization solver can be a linear programming solver, a mixed integer programming solver, or a nonlinear programming solver.

[0095] The energy island planning and design method based on mathematical optimization models in this embodiment constructs a unified energy flow topology covering power generation, energy storage, hydrogen production, and various types of loads. It establishes an optimization model that includes supply and demand balance, capacity constraints, energy conversion efficiency, and environmental constraints. By combining multi-objective functions with commercial solvers, it obtains the optimal scheme for energy island capacity configuration and operation scheduling. This improves energy utilization efficiency and environmental friendliness while ensuring economy and reliability, providing scientific decision support for the construction and operation of energy islands in different scenarios.

[0096] The mathematical model of the energy island system is as follows:

[0097] maxIRR = f IRR (cash)

[0098]

[0099] st{electric balance constraints, hydrogen balance constraints, ammonia balance constraints, alcohol balance constraints, electricity purchase constraints, electrical energy storage constraints, hydrogen equipment constraints, synthetic ammonia equipment constraints, alcohol production equipment constraints, energy island constraints, hydrogen storage constraints, ammonia storage constraints, and alcohol storage constraints};

[0100] In the formula, IRR is the internal rate of return on capital, cash is the total cash flow of the energy island system, and f IRR The formula for calculating the internal rate of return (IRR) of capital is e. abandon For the abandoned electricity of the energy island system, yPG ki Let xRPV be the electricity purchased in the i-th hour of year k, xRPV be the installed photovoltaic capacity on the island, and effPV be the electricity purchased in the i-th hour of year k. i Let capPV be the photovoltaic power generation coefficient in the i-th hour. k Let yBES be the remaining capacity factor of the wind turbine in year k. ki Let OE be the charge / discharge amount of electrical energy storage in the i-th hour of year k. ki Let lE be the amount of electricity generated by ocean energy in the i-th hour of year k. i Let ySA be the electrical load for the i-th hour. ki Let y_cSA be the ammonia production capacity of the ammonia synthesis unit in the i-th hour of year k, and y_ECU be the unit power consumption of the ammonia synthesis unit. ki Let be the ammonia production output of the hydrogen production unit in the i-th hour of year k, and ecECU be the unit power consumption of the hydrogen production unit, yMMW. ki Let x represent the alcohol production capacity of the alcohol production equipment in the i-th hour of year k, ecMMW represent the unit power consumption of the ammonia production equipment, xCPC represent the installed capacity of the computing center, ecCPC represent the unit power consumption of the computing center, xDAS represent the installed capacity of the seawater desalination equipment, ecDAS represent the unit power consumption of the seawater desalination equipment, xFE represent the installed capacity of the refueling equipment, ecFE represent the unit power consumption of the refueling equipment, xOMH represent the construction scale of the maintenance home port, ecOMH represent the unit power consumption of the maintenance home port, xMRC represent the construction scale of the rescue center, ecMRC represent the unit power consumption of the rescue center, xCAT represent the construction scale of the cultural tourism area, ecCAT represent the unit power consumption of the cultural tourism area, xMCT represent the construction scale of the aquaculture area, ecMCT represent the unit power consumption of the aquaculture area, and eCELD represent the unit power consumption of the energy island.

[0101] In one embodiment of the formula, the total cash flow of the energy island system is calculated as follows:

[0102]

[0103] cash = [cash1, cash2, ..., cash] k ,...,cash K ]

[0104] In the formula, cash k For the cash flow in year k, reward k For the income in year k, Let the investment cost be in year k. The maintenance cost in year k is... Let K be the electricity cost in year k, and K be the investment period of the project.

[0105] In one embodiment of the formula, the investment cost calculation formula for year k is:

[0106]

[0107] In the formula, Let x_k be the investment cost coefficient for the photovoltaic equipment on the island in year k, and x_ECU be the installed capacity of the hydrogen production equipment. Let xSA be the investment cost coefficient for hydrogen production equipment in year k, and xSA be the installed capacity of ammonia synthesis equipment. Let xMMW be the investment cost coefficient for the ammonia synthesis equipment in year k, and xMMW be the installed capacity of the methanol production equipment. Let be the investment cost coefficient for the alcohol production equipment in year k. Let be the investment cost coefficient of the computing center in year k. Let xELD be the investment cost coefficient for the hydrogen production equipment in year k, and xELD be the total area of ​​the energy island. Let xOE be the investment cost coefficient for the energy island in year k, and xOE be the installed capacity of the marine hub. Let xPH be the investment cost coefficient for ocean energy in year k, and xPH be the installed capacity of the power hub. Let be the investment cost coefficient of the power hub in year k. This represents the investment cost coefficient for the equipment in year k. This represents the investment cost coefficient for the maintenance of the home port in year k. Let be the investment cost coefficient of the rescue center in year k. Let be the investment cost coefficient for cultural tourism in year k. Let xSBAT be the investment cost coefficient for aquaculture in year k, and let xSBAT be the installed capacity of electrical energy storage. Let xSH be the investment cost coefficient for electric energy storage in year k, and xSH be the installed capacity of hydrogen storage equipment. Let x be the investment cost coefficient for hydrogen storage in year k, and xHSA be the installed capacity of ammonia storage equipment. Let xHMW be the investment cost coefficient for the ammonia storage equipment in year k, and let xHMW be the installed capacity of the methanol storage equipment. Let be the investment cost coefficient for alcohol storage in year k.

[0108] In one embodiment of the formula, the maintenance cost calculation formula for year k is:

[0109]

[0110] In the formula, Let be the maintenance cost coefficient for the photovoltaic equipment on the island in year k. Let be the maintenance cost coefficient for the hydrogen production equipment in year k. Let be the maintenance cost coefficient for the ammonia synthesis equipment in year k. Let be the maintenance cost coefficient for the alcohol production equipment in year k. Let be the maintenance cost coefficient of the computing center in year k. Let be the maintenance cost coefficient for the hydrogen production equipment in year k. Let be the maintenance cost coefficient of the energy island in year k. Let be the maintenance cost coefficient for ocean energy in year k. Let be the maintenance cost coefficient of the power hub in year k. This represents the maintenance cost coefficient for the refueling equipment in year k. This represents the maintenance cost coefficient for the home port in year k. Let this be the maintenance cost coefficient of the rescue center in year k. Let be the maintenance cost coefficient for cultural tourism in year k. Let this be the maintenance cost coefficient for aquaculture in year k. Let be the maintenance cost coefficient for energy storage in year k. Let be the maintenance cost coefficient for hydrogen storage in year k. Let be the maintenance cost coefficient for the ammonia storage equipment in year k. Let be the maintenance cost coefficient for alcohol storage in year k.

[0111] In one embodiment of the formula, the electricity cost calculation formula for year k is:

[0112]

[0113] In the formula, pe i The electricity cost for the i-th hour;

[0114] The formula for calculating the return in year k is:

[0115]

[0116] In the formula, se i sh i ssa i smmw i Let LE be the selling price of electricity, hydrogen, ammonia, and alcohol in the i-th hour. i For the electrical load in the i-th hour, lH i For the hydrogen load in the i-th hour, lSA i For the ammonia load in the i-th hour, lMME i For the alcohol load in the i-th hour, Let the unit revenue of the computing center be at the k-th position. Let $\frac{ ... Let the unit revenue of the refueling equipment be in year k. For the unit revenue of the maintenance home port in year k, For the unit revenue of the rescue center in year k, Let the unit revenue of cultural tourism in year k be... Let $\frac{ ...

[0117] In one embodiment of the formula, the electrical balance constraint is:

[0118]

[0119] The hydrogen balance constraint is:

[0120] yECU ki +ySH ki ≥

[0121] lH i +ySA ki ×hcSA+yMMW ki ×hcMMW

[0122] In the formula, ySH ki Let hcSA be the amount of hydrogen stored and released by the hydrogen storage equipment in the i-th hour of year k, hcMA be the unit hydrogen consumption of the ammonia synthesis equipment, and hcMMW be the unit hydrogen consumption of the alcohol production equipment.

[0123] The ammonia balance constraint is:

[0124] ySA ki +yHSA ki ≥lSA i

[0125] In the formula, yHSA ki Let be the amount of ammonia stored and released by the ammonia storage equipment in the i-th hour of year k.

[0126] The equilibrium constraint for alcohols is:

[0127] yMMW ki +yHMW ki ≥lMMW i

[0128] In the formula, yHMW ki Let be the amount of ammonia stored and released by the methanol storage device in the i-th hour of year k.

[0129] The power purchase limit is:

[0130] yPG ki ≤xSPV×effSPV i ×capSPVk +xSWT×effSWT i

[0131] In the formula, xSPV represents the installed capacity of offshore photovoltaic power, and effSPV represents the installed capacity of offshore photovoltaic power. i Let capSPV be the power generation coefficient of offshore photovoltaic power generation in the i-th hour. k Let xSWT be the remaining capacity factor of offshore photovoltaic power in year k, xSWT be the installed capacity of offshore wind turbines, and effSWT be the remaining capacity factor. i Let be the power generation coefficient of offshore photovoltaic power in the i-th hour;

[0132] The constraints of electrical energy storage are:

[0133] yBES ki ≤zBES ki

[0134] -yBES ki ≤xBES-zBES ki

[0135] zBES k,i+1 =zBES ki -yBES ki

[0136]

[0137] In the formula, zBES ki Let xBES be the remaining energy storage capacity at time i in year k, and zBES be the installed capacity of energy storage. k,i+1 Let be the remaining electrical energy stored at time i+1 in year k. The charging rate for electrical energy storage. The discharge rate of the stored energy;

[0138] Hydrogen equipment constraints are: yECU ki ≤xECU;

[0139] The constraints for ammonia synthesis equipment are: ySA ki ≤xSA;

[0140] The constraints for the alcohol production equipment are: yMMW ki ≤xMMW;

[0141] In one embodiment of the formula, the energy island constraint is:

[0142] xELD=xRPV×aRPV+xOE×aOE+xECU×aECU+

[0143] xSA×aSA+xMMW×aMMW+xSH×aSH+

[0144] xHSA×aHSA+xHMW×aHMW+xCPC×aCPC+

[0145] xDAS×aDAS+xPH×aPH+xFE×aFE+

[0146] xOMH×aOMH+xMRC×aMRC+xCAT×aCAT+

[0147] xMCT×aMCT

[0148] In the formula, aRPV represents the area occupied by the unit installed capacity of onshore photovoltaics, aOE represents the area occupied by the unit installed capacity of ocean energy, aECU represents the area occupied by the unit installed capacity of hydrogen production equipment, aSA represents the area occupied by the unit installed capacity of ammonia synthesis equipment, aMMW represents the area occupied by the unit installed capacity of methanol production equipment, aSH represents the area occupied by the unit installed capacity of hydrogen storage equipment, aHSA represents the area occupied by the unit installed capacity of ammonia storage equipment, aHMW represents the area occupied by the unit installed capacity of methanol storage, aCPC represents the area occupied by the unit installed capacity of computing center, aDAS represents the area occupied by the unit installed capacity of seawater desalination, aPH represents the area occupied by the unit installed capacity of power hub, aFE represents the area occupied by the unit installed capacity of refueling equipment, aOMH represents the area occupied by the unit construction capacity of operation and maintenance home port, aMRC represents the area occupied by the unit construction capacity of rescue center, aCAT represents the area occupied by the unit construction capacity of cultural tourism, and aMCT represents the area occupied by the unit construction capacity of aquaculture.

[0149] In one embodiment of the formula, the hydrogen storage constraint is:

[0150] ySH ki ≤zSH ki

[0151] -ySH ki ≤xSH-zSH ki

[0152] zSH k,i+1 =zSH ki -ySH ki

[0153]

[0154] In the formula, zSH ki Let zSH be the remaining hydrogen quantity of the hydrogen storage device at time i in year k. k,i+1 Let be the remaining amount of hydrogen in the hydrogen storage device at time i+1 in year k. For hydrogen storage ratio, This refers to the hydrogen release rate;

[0155] Ammonia storage constraints are:

[0156] yHSA ki≤zHSA ki

[0157] -yHSA ki ≤xHSA-zHSA ki

[0158] zHSA k,i+1 =zHSA ki -yHSA ki

[0159]

[0160] In the formula, zHSA ki Let zHSA be the amount of remaining ammonia in the ammonia storage device at time i in year k. k,i+1 Let be the amount of ammonia remaining in the ammonia storage device at time i+1 in year k. This refers to the ammonia storage ratio. This refers to the ammonia release ratio;

[0161] The alcohol storage constraints are:

[0162] yHMW ki ≤zHMW ki

[0163] -yHMW ki ≤xHMW-zHMW ki

[0164] zHMW k,i+1 =zHMW ki -yHMW ki

[0165]

[0166] In the formula, zHMW ki Let zHMW be the amount of alcohol remaining in the alcohol storage device at time i in year k. k,i+1 Let be the amount of alcohol remaining in the alcohol storage device at time i+1 in year k. This refers to the ammonia storage ratio. This represents the alcohol release ratio.

[0167] In one embodiment of the formula,

[0168] It should be understood that although the steps in the flowcharts of the embodiments described above are shown sequentially according to the arrows, these steps are not necessarily executed in the order indicated by the arrows. Unless explicitly stated herein, there is no strict order restriction on the execution of these steps, and they can be executed in other orders. Moreover, at least some steps in the flowcharts of the embodiments described above may include multiple steps or multiple stages. These steps or stages are not necessarily completed at the same time, but can be executed at different times. The execution order of these steps or stages is not necessarily sequential, but can be performed alternately or in turn with other steps or at least some of the steps or stages of other steps.

[0169] The technical features of the above embodiments can be combined in any way. For the sake of brevity, not all possible combinations of the technical features in the above embodiments are described. However, as long as there is no contradiction in the combination of these technical features, they should be considered to be within the scope of this specification.

[0170] The embodiments described above are merely illustrative of several implementations of the present invention, and while the descriptions are specific and detailed, they should not be construed as limiting the scope of the present invention. It should be noted that those skilled in the art can make various modifications and improvements without departing from the concept of the present invention, and these modifications and improvements all fall within the scope of protection of the present invention. Therefore, the scope of protection of the present invention should be determined by the appended claims.

Claims

1. A method for planning and designing energy islands based on a mathematical optimization model, characterized in that, Includes the following steps: Establish the topology diagram of the energy island system; Based on the topology diagram of the energy island system, a mathematical model of the energy island system is constructed with the objectives of maximizing the internal rate of return on capital and minimizing the amount of abandoned electricity. The mathematical model of the energy island system is input into the optimization solver to obtain the optimal configuration and operation scheme of the energy island system; The mathematical model of the energy island system is as follows: maxIRR=f IRR (cash) st{electric balance constraints, hydrogen balance constraints, ammonia balance constraints, alcohol balance constraints, electricity purchase constraints, electrical energy storage constraints, hydrogen equipment constraints, synthetic ammonia equipment constraints, alcohol production equipment constraints, energy island constraints, hydrogen storage constraints, ammonia storage constraints, and alcohol storage constraints}; In the formula, IRR is the internal rate of return on capital, cash is the total cash flow of the energy island system, and f IRR The formula for calculating the internal rate of return (IRR) of capital is e. abandon For the abandoned electricity of the energy island system, yPG ki Let xRPV be the electricity purchased in the i-th hour of year k, xRPV be the installed photovoltaic capacity on the island, and effPV be the electricity purchased in the i-th hour of year k. i Let capPV be the photovoltaic power generation coefficient in the i-th hour. k Let yBES be the remaining capacity factor of the wind turbine in year k. ki Let OE be the charge / discharge amount of electrical energy storage in the i-th hour of year k. ki Let lE be the amount of electricity generated by ocean energy in the i-th hour of year k. i Let ySA be the electrical load for the i-th hour. ki Let y_cSA be the ammonia production capacity of the ammonia synthesis unit in the i-th hour of year k, and y_ECU be the unit power consumption of the ammonia synthesis unit. ki Let be the ammonia production output of the hydrogen production unit in the i-th hour of year k, and ecECU be the unit power consumption of the hydrogen production unit, yMMW. ki Let x represent the alcohol production capacity of the alcohol production equipment in the i-th hour of year k, ecMMW represent the unit power consumption of the ammonia production equipment, xCPC represent the installed capacity of the computing center, ecCPC represent the unit power consumption of the computing center, xDAS represent the installed capacity of the seawater desalination equipment, ecDAS represent the unit power consumption of the seawater desalination equipment, xFE represent the installed capacity of the refueling equipment, ecFE represent the unit power consumption of the refueling equipment, xOMH represent the construction scale of the maintenance home port, ecOMH represent the unit power consumption of the maintenance home port, xMRC represent the construction scale of the rescue center, ecMRC represent the unit power consumption of the rescue center, xCAT represent the construction scale of the cultural tourism area, ecCAT represent the unit power consumption of the cultural tourism area, xMCT represent the construction scale of the aquaculture area, ecMCT represent the unit power consumption of the aquaculture area, and eCELD represent the unit power consumption of the energy island.

2. The energy island planning and design method based on a mathematical optimization model according to claim 1, characterized in that, The formula for calculating the total cash flow of the energy island system is: cash=[cash1,cash2,...,cash k ,...,cash K ] In the formula, cash k For the cash flow in year k, reward k For the income in year k, Let the investment cost be in year k. The maintenance cost in year k is... Let K be the electricity cost in year k, and K be the investment period of the project.

3. The energy island planning and design method based on a mathematical optimization model according to claim 2, characterized in that, The formula for calculating the investment cost in year k is: In the formula, Let x_k be the investment cost coefficient for the photovoltaic equipment on the island in year k, and x_ECU be the installed capacity of the hydrogen production equipment. Let xSA be the investment cost coefficient for hydrogen production equipment in year k, and xSA be the installed capacity of ammonia synthesis equipment. Let xMMW be the investment cost coefficient for the ammonia synthesis equipment in year k, and xMMW be the installed capacity of the methanol production equipment. Let be the investment cost coefficient for the alcohol production equipment in year k. Let be the investment cost coefficient of the computing center in year k. Let xELD be the investment cost coefficient for the hydrogen production equipment in year k, and xELD be the total area of ​​the energy island. Let xOE be the investment cost coefficient for the energy island in year k, and xOE be the installed capacity of the marine hub. Let xPH be the investment cost coefficient for ocean energy in year k, and xPH be the installed capacity of the power hub. Let be the investment cost coefficient of the power hub in year k. This represents the investment cost coefficient for the equipment in year k. This represents the investment cost coefficient for the maintenance of the home port in year k. Let be the investment cost coefficient of the rescue center in year k. Let be the investment cost coefficient for cultural tourism in year k. Let xSBES be the investment cost coefficient for aquaculture in year k, and let xSBES be the installed capacity of electrical energy storage. Let xSH be the investment cost coefficient for electric energy storage in year k, and xSH be the installed capacity of hydrogen storage equipment. Let x be the investment cost coefficient for hydrogen storage in year k, and xHSA be the installed capacity of ammonia storage equipment. Let xHMW be the investment cost coefficient for the ammonia storage equipment in year k, and let xHMW be the installed capacity of the methanol storage equipment. Let be the investment cost coefficient for alcohol storage in year k.

4. The energy island planning and design method based on a mathematical optimization model according to claim 2, characterized in that, The formula for calculating the maintenance cost in year k is: In the formula, Let be the maintenance cost coefficient for the photovoltaic equipment on the island in year k. Let be the maintenance cost coefficient for the hydrogen production equipment in year k. Let be the maintenance cost coefficient for the ammonia synthesis equipment in year k. Let be the maintenance cost coefficient for the alcohol production equipment in year k. Let be the maintenance cost coefficient of the computing center in year k. Let be the maintenance cost coefficient for the hydrogen production equipment in year k. Let be the maintenance cost coefficient of the energy island in year k. Let be the maintenance cost coefficient for ocean energy in year k. Let be the maintenance cost coefficient of the power hub in year k. This represents the maintenance cost coefficient for the refueling equipment in year k. This represents the maintenance cost coefficient for the home port in year k. Let this be the maintenance cost coefficient of the rescue center in year k. Let be the maintenance cost coefficient for cultural tourism in year k. Let this be the maintenance cost coefficient for aquaculture in year k. Let be the maintenance cost coefficient for energy storage in year k. Let be the maintenance cost coefficient for hydrogen storage in year k. Let be the maintenance cost coefficient for the ammonia storage equipment in year k. Let be the maintenance cost coefficient for alcohol storage in year k.

5. The energy island planning and design method based on a mathematical optimization model according to claim 2, characterized in that, The formula for calculating the electricity cost in year k is: In the formula, pe i The electricity cost for the i-th hour; The formula for calculating the return in year k is: In the formula, se i sh i ssa i smmw i Let LE be the selling price of electricity, hydrogen, ammonia, and alcohol in the i-th hour. i For the electrical load in the i-th hour, lH i For the hydrogen load in the i-th hour, lSA i For the ammonia load in the i-th hour, lMME i For the alcohol load in the i-th hour, Let the unit revenue of the computing center be at the k-th position. Let $\frac{ ... Let the unit revenue of the refueling equipment be in year k. For the unit revenue of the maintenance home port in year k, For the unit revenue of the rescue center in year k, Let the unit revenue of cultural tourism in year k be... Let $\frac{ ...

6. The energy island planning and design method based on a mathematical optimization model according to claim 1, characterized in that, The electrical balance constraint is: The hydrogen balance constraint is: yECU ki +ySH ki ≥ lH i +ySA ki ×hcSA+yMMW ki ×hcMMW In the formula, ySH ki Let hcSA be the amount of hydrogen stored and released by the hydrogen storage equipment in the i-th hour of year k, hcMA be the unit hydrogen consumption of the ammonia synthesis equipment, and hcMMW be the unit hydrogen consumption of the alcohol production equipment. The ammonia balance constraint is: ySA ki +yHSA ki ≥lSA i In the formula, yHSA ki Let be the amount of ammonia stored and released by the ammonia storage equipment in the i-th hour of year k. The alcohol equilibrium constraint is: yMMW ki +yHMW ki ≥lMMW i In the formula, yHMW ki Let be the amount of ammonia stored and released by the methanol storage device in the i-th hour of year k. The constraint on the amount of electricity purchased is: yPG ki ≤xSPV×effSPV i ×capSPV k +xSWT×effSWT i In the formula, xSPV represents the installed capacity of offshore photovoltaic power, and effSPV represents the installed capacity of offshore photovoltaic power. i Let capSPV be the power generation coefficient of offshore photovoltaic power generation in the i-th hour. k Let xSWT be the remaining capacity factor of offshore photovoltaic power in year k, xSWT be the installed capacity of offshore wind turbines, and effSWT be the remaining capacity factor. i Let be the power generation coefficient of offshore photovoltaic power in the i-th hour; The energy storage constraint is: yBES ki ≤zBES ki -yBES ki ≤xBES-zBES ki zBES k,i+1 =zBES ki -yBES ki In the formula, zBES ki Let xBES be the remaining energy storage capacity at time i in year k, and zBES be the installed capacity of energy storage. k,i+1 Let be the remaining electrical energy stored at time i+1 in year k. The charging rate for electrical energy storage. The discharge rate of the stored energy; The hydrogen equipment constraint is: yECU ki ≤xECU; The constraints of the ammonia synthesis equipment are: ySA ki ≤xSA; The constraints of the alcohol production equipment are: yMMW ki ≤xMMW.

7. The energy island planning and design method based on a mathematical optimization model according to claim 1, characterized in that, The constraints of the energy island are: xELD=xRPV×aRPV+xOE×aOE+xECU×aECU+ xSA×aSA+xMMW×aMMW+xSH×aSH+ xHSA×aHSA+xHMW×aHMW+xCPC×aCPC+ xDAS×aDAS+xPH×aPH+xFE×aFE+ xOMH×aOMH+xMRC×aMRC+xCAT×aCAT+ xMCT×aMCT In the formula, aRPV represents the area occupied by the unit installed capacity of onshore photovoltaics, aOE represents the area occupied by the unit installed capacity of ocean energy, aECU represents the area occupied by the unit installed capacity of hydrogen production equipment, aSA represents the area occupied by the unit installed capacity of ammonia synthesis equipment, aMMW represents the area occupied by the unit installed capacity of methanol production equipment, aSH represents the area occupied by the unit installed capacity of hydrogen storage equipment, aHSA represents the area occupied by the unit installed capacity of ammonia storage equipment, aHMW represents the area occupied by the unit installed capacity of methanol storage, aCPC represents the area occupied by the unit installed capacity of computing center, aDAS represents the area occupied by the unit installed capacity of seawater desalination, aPH represents the area occupied by the unit installed capacity of power hub, aFE represents the area occupied by the unit installed capacity of refueling equipment, aOMH represents the area occupied by the unit construction capacity of operation and maintenance home port, aMRC represents the area occupied by the unit construction capacity of rescue center, aCAT represents the area occupied by the unit construction capacity of cultural tourism, and aMCT represents the area occupied by the unit construction capacity of aquaculture.

8. The energy island planning and design method based on a mathematical optimization model according to claim 1, characterized in that, The hydrogen storage constraint is: ySH ki ≤zSH ki -ySH ki ≤xSH-zSH ki zSH k,i+1 =zSH ki -ySH ki In the formula, zSH ki Let zSH be the remaining hydrogen quantity of the hydrogen storage device at time i in year k. k,i+1 Let be the remaining amount of hydrogen in the hydrogen storage device at time i+1 in year k. For hydrogen storage ratio, This refers to the hydrogen release rate; The ammonia storage constraint is: yHSA ki ≤zHSA ki -yHSA ki ≤xHSA-zHSA ki zHSA k,i+1 =zHSA ki -yHSA ki In the formula, zHSA ki Let zHSA be the amount of remaining ammonia in the ammonia storage device at time i in year k. k,i+1 Let be the amount of ammonia remaining in the ammonia storage device at time i+1 in year k. This refers to the ammonia storage ratio. This refers to the ammonia release ratio; The alcohol storage constraint is: yHMW ki ≤zHMW ki -yHMW ki ≤xHMW-zHMW ki zHMW k,i+1 =zHMW ki -yHMW ki In the formula, zHMW ki Let zHMW be the amount of alcohol remaining in the alcohol storage device at time i in year k. k,i+1 Let be the amount of alcohol remaining in the alcohol storage device at time i+1 in year k. This refers to the ammonia storage ratio. This represents the alcohol release ratio.

9. The energy island planning and design method based on a mathematical optimization model according to claim 1, characterized in that, The optimal configuration and operation scheme of the energy island system includes: The capacity configuration of various production, storage, and consumption devices, as well as the optimal construction scale of various infrastructures in the topology; Energy island energy flow and operation scheduling strategies, including the production, storage and consumption of electricity and hydrogen energy at different time scales; The system's economic and environmental indicators include total life cycle cost, return on equity, carbon emission levels, and renewable energy utilization rate.

10. The energy island planning and design method based on a mathematical optimization model according to claim 1, characterized in that, The optimization solver is a linear programming solver, a mixed integer programming solver, or a nonlinear programming solver.

Images