Island and ship integrated energy system two-stage robust optimization method

Abstract

The application discloses a kind of island and ship integrated energy system two-stage robust optimization method, which comprises the following steps: 1) establishing the objective function of island and ship integrated energy system two-stage robust optimization model;2) establish the uncertain set of island and ship integrated energy system two-stage robust optimization model;3) establish the constraint condition of ship integrated energy system in island and ship integrated energy system two-stage robust optimization model;4) establish the constraint condition of energy equipment in island and ship integrated energy system two-stage robust optimization model, including fuel cell, electrolytic cell, gas turbine, energy storage device constraint condition;5) establish the power balance constraint condition of each stage in island and ship integrated energy system two-stage robust optimization model;6) solve island and ship integrated energy system two-stage robust optimization model using nested column constraint generation algorithm, obtain island and ship integrated energy system joint optimization scheme.

Description

Technical Field

[0001] This invention belongs to the field of integrated energy technology, and specifically relates to a two-stage robust optimization method for integrated energy systems of islands and ships. Background Technology

[0002] Driven by global pressures for energy transition and carbon emission reduction in the shipping industry, the large-scale application of electricity and hydrogen energy equipment in island scenarios has become a trend. However, the development of renewable energy on islands still faces multiple challenges: on the one hand, limited local load capacity leads to high curtailment rates, with diesel generator sets still dominating; on the other hand, uneven spatial distribution of energy sources and loads among multiple islands exacerbates the complexity of energy transmission. Hydrogen energy, as a high-energy-density energy storage carrier, plays a crucial role in promoting multi-energy flow coupling and renewable energy consumption in integrated island energy systems, but its large-scale application introduces new challenges such as increased probability of auxiliary equipment failure and increased difficulty in dynamic system matching.

[0003] As a novel mobile energy carrier, the integrated ship energy system offers a technological solution to the aforementioned problems by participating in the conversion and allocation of electricity and hydrogen energy between islands. However, its operation is constrained by uncertainties such as complex sea conditions and equipment failures. For example, fluctuations in speed directly affect the timeliness of energy transmission, while uneven distribution of energy sources and loads further amplifies the impact of failures. Traditional independent planning models for single islands or ships are no longer sufficient to meet the system's resilience requirements. Although existing research has proposed a multi-island cluster collaborative energy supply architecture and a multi-energy flow interaction mechanism for ships, there are still theoretical gaps in the collaborative planning and operational strategy integration of ship and island systems. In particular, a complete framework has not yet been established for scenarios with a high proportion of hydrogen energy access, resulting in limited system adaptability under extreme events such as equipment failures and uncertainties in navigation time. In addition, the considerable speed loss caused by the uncertainty of ship navigation hinders the timely response of emergency resources. These challenges will lead to severe energy supply disruptions, requiring multi-energy islands and ship microgrids to optimize energy equipment capacity to ensure safe and stable operation.

[0004] Against this backdrop, a two-stage robust optimization method offers an innovative approach to enhancing system robustness through the coordinated optimization of equipment capacity configuration and ship navigation scheduling. This method combines the strong adaptability of robust optimization to uncertainty with the phased flexibility of two-stage decision-making. In the first stage, it achieves optimal planning of equipment capacity, while in the second stage, it coordinates ship energy interaction and navigation scheduling strategies, effectively balancing economy and reliability. Through the coordinated optimization of equipment capacity and ship scheduling, it effectively improves the system's ability to cope with equipment failures and navigation time uncertainties, thereby enhancing the overall robustness of island and ship integrated energy systems. Summary of the Invention

[0005] Technical Solution: To address the technical problems mentioned in the background section, this invention proposes a two-stage robust optimization method for integrated energy systems of islands and ships, which includes the following steps:

[0006] (1) Establish the objective function of a two-stage robust optimization model for integrated energy systems of islands and ships;

[0007] (2) Establish the uncertainty set of a two-stage robust optimization model for integrated energy systems of islands and ships;

[0008] (3) Establish the constraints of the integrated energy system of the ship in the two-stage robust optimization model of the island and ship integrated energy system, including the navigation scheduling, propulsion system and energy interaction constraints of the integrated energy system of the ship;

[0009] (4) Establish constraints on energy equipment in the two-stage robust optimization model of the integrated energy system of islands and ships, including constraints on fuel cells, electrolyzers, gas turbines and energy storage equipment;

[0010] (5) Establish the power balance constraints for each stage in the two-stage robust optimization model of the island and ship integrated energy system, including the power balance constraints for the first stage and the second stage.

[0011] (6) The nested column constraint generation algorithm is used to solve the two-stage robust optimization model of the island and ship integrated energy system, and the joint optimization scheme of the island and ship integrated energy system is obtained.

[0012] Furthermore, in step (1), the objective function of the two-stage robust optimization model for the integrated energy system of the island and the ship is established, as follows:

[0013]

[0014] In the formula, x represents the planning stage variable; κ represents the depreciation rate; n represents the system's operating life; y represents the scenario; θ represents the energy equipment on the island, which is one of the following: wind turbine, electrolyzer, fuel cell, hydrogen storage equipment, or gas turbine; i represents the island. Let θ be the unit investment cost of the energy equipment; γ represents the planned capacity of equipment θ on island i; γ represents the ship's energy equipment, which is one of fuel cells, energy storage equipment, and hydrogen storage equipment. The unit investment cost of ship energy equipment γ; The planned capacity for ship energy equipment γ; Ω NOPE N is a set of scenarios during normal operation. y c is the total number of days in scenario y; OPE The cost of the operation phase; t is time point; The electrical power of energy equipment θ on island i at time t during normal operation; These represent the charging power and discharging power of the energy device θ on island i at time t during normal operation, respectively. Let γ be the electrical power of the ship's energy equipment at time t during normal operation. These represent the charging and discharging power of the ship's energy equipment γ at time t during normal operation; Δt is the interval time; T is a time period; c WT,CUR The unit cost of wind curtailment for wind turbine units; The wind curtailment power of the wind turbines on island i at time t during normal operation. Ω represents the uncertainty set of a two-stage robust optimization model for integrated energy systems of islands and ships; z represents the variables during the emergency operation phase; EOPE This is a set of scenarios for the emergency operation phase, which takes into account the uncertainty of equipment failure and ship sailing time. The electrical power of the energy equipment θ on island i at time t during the emergency operation phase; These represent the charging power and discharging power of the energy equipment θ on island i at time t during the emergency operation phase, respectively. The electrical power of the ship's energy equipment γ at time t during the emergency operation phase; These represent the charging power and discharging power of the ship's energy equipment γ at time t during the emergency operation phase, respectively. c represents the wind curtailment power of the wind turbines on island i at time t during the emergency operation phase. ELS,EOPE Cost per unit of power load shedding; Let c be the power load shedding amount on island i at time t; HLS,EOPE Cost per unit of hydrogen energy load shedding; Let be the hydrogen energy load shedding amount on island i at time t.

[0015] Furthermore, in step (2), the uncertainty set of the two-stage robust optimization model for the island and ship integrated energy system is established, as follows:

[0016]

[0017]

[0018] In the formula, For a set of uncertain sailing times; For equipment failure uncertainty set; σ represents the ship's actual speed at time t; k1 and k2 are the power calculation coefficients for the ship's propulsion system; PS To improve system efficiency; P t EV,PS Let t be the output power of the all-electric ship's motor. For resistance during still water navigation; Adds resistance to the waves; W V R is the beam of the ship. V L is the block coefficient of the ship. V ρ is the length of the ship; g is the density of seawater; and g is the acceleration due to gravity. These are the actual and predicted values ​​of the wave height, respectively. This is the positive offset of the wave height; This is the negative offset of the wave height; This is a binary variable, indicating whether the wave height has shifted positively at time t; a positive shift is 1, otherwise it is 0. This is a binary variable representing whether the wave height has a negative shift at time t; a value of 1 indicates a negative shift, and a value of 0 indicates otherwise. SV For uncertain budgets related to ship navigation; This is a binary variable indicating whether the electrolytic cell is faulty; 1 indicates a fault, and 0 indicates otherwise. This is a binary variable indicating whether the fuel cell is faulty; 1 indicates a fault, and 0 indicates otherwise. This is a binary variable representing whether the gas turbine is faulty; 1 indicates a fault, and 0 indicates otherwise. EC Uncertain budget for emergency equipment use.

[0019] Furthermore, the specific process of step (3) is as follows:

[0020] (301) Establish navigation scheduling constraints for the ship's integrated energy system

[0021]

[0022] X i,t =1,t=T i LA (6)

[0023]

[0024] In the formula, i+i' represents all islands except the starting and ending points; X i,t Let X be the position variable of the ship relative to island i at time t. i,t =1 indicates that the ship is docked on island i at time t, X i,t =0 indicates that the ship departs from island i at time t and is in a sailing state; X i,t+1 Let X be the position variable of the ship relative to island i at time t+1. i,t+1 =1 means the ship is docked on island i at time t+1, X i,t+1 =0 indicates that the ship departs from island i at time t+1 and is in a sailing state; T MBD This refers to the minimum berthing time for a vessel on an island. These represent the latest times when ships arrive at islands i-1, i, and i+1, respectively; R i,t R i,t+1 These represent the real-time sailing distances of the ship on its route starting from island i at times t and t+1, respectively. Let be the ship's speed at time t; This represents the total distance traveled by the vessel along the route originating from island i.

[0025] (302) Establishing propulsion system constraints for integrated ship energy systems

[0026]

[0027] In the formula, This indicates the start / stop status of the motor; 1 indicates start, and 0 indicates stop.

[0028] (303) Establishing energy interaction constraints for integrated ship energy systems

[0029]

[0030] In the formula, P represents the power exchange between the island's i and the ship's electrical energy. EEXC,max This represents the upper limit of the electrical power exchange between island i and the ship; For the hydrogen energy interaction power between island i and ship; P HEXC,max This represents the upper limit of the hydrogen energy interaction power between the island i and the ship.

[0031] Furthermore, the specific process of step (4) is as follows:

[0032] (401) Establish fuel cell constraints

[0033]

[0034] P FC,max =l1·N FC ·(α FC ·C FC,max +β FC (18)

[0035] In the formula, P t FC L1 represents the output power of the fuel cell at time t; L1 and L2 are the unit conversion factors; α FC With β FC The fitting coefficients for the operating power of the fuel cell are denoted as . N is the load current of the fuel cell at time t; FC This represents the number of fuel cell stacks; ω represents the start / stop state variable of the fuel cell; 1 indicates start-up and 0 indicates otherwise. FCThis is the operating state variable of the fuel cell; it is 1 for normal operation and 0 otherwise. Let be the hydrogen power consumed by the fuel cell at time t; F be the Faraday constant; C be the hydrogen power consumed by the fuel cell at time t. FC,min C FC,max These are the minimum and maximum load currents of the fuel cell, respectively; P FC,max This refers to the rated capacity of the fuel cell;

[0036] (402) Establishing constraints for the electrolytic cell

[0037]

[0038] P ED,max =l1·N ED ·(α ED ·C ED,max +β ED ) (twenty two)

[0039] In the formula, P t ED α is the input power of the electrolyzer at time t; ED With β ED The fitting coefficients for the working power of the electrolytic cell; N represents the load current of the electrolytic cell at time t. ED This refers to the number of electrolysis units; This represents the start / stop status of the electrolytic cell; 1 indicates start-up, and 0 indicates otherwise. ED This is the operating status variable of the electrolysis unit; it is 1 for normal operation and 0 otherwise. σ is the output power of the electrolytic cell; ED For the working efficiency of the electrolytic cell; C ED,min C ED,max These are the minimum and maximum load currents of the electrolytic cell, respectively; P ED,max This refers to the rated capacity of the electrolytic cell;

[0040] (403) Establish gas turbine constraints

[0041]

[0042] In the formula, σ GT,max P represents the maximum gradeability of the gas turbine. GT,max P represents the rated capacity of the gas turbine. t GT , These represent the output power of the gas turbine at times t and t+1, respectively; ζ GT,min ζ GT,max These are the minimum and maximum operating ratios of the gas turbine, respectively. ω represents the start-stop state variable of the gas turbine; 1 indicates start-up and 0 indicates shutdown. GTThis is the operating state variable of the gas turbine; it is 1 for normal operation and 0 otherwise. The hydrogen consumption power of the gas turbine; σ GT For the working efficiency of the gas turbine;

[0043] (404) Establish constraints for energy storage devices:

[0044]

[0045] In the formula, π represents the energy storage device, which can be either an electrical storage device or a hydrogen storage device; V π,t V π,t+1 Let π be the self-discharge of the energy storage device at times t and t+1, respectively; σ π π represents the self-discharge rate of the energy storage device. These represent the charging power and discharging power of the energy storage device at time t, respectively. These represent the charging efficiency and discharging efficiency of the energy storage device π, respectively. π represents the maximum ramp rate of the energy storage device. The rated capacity of the energy storage device π; These represent the minimum and maximum operating coefficients (π) of the energy storage device, respectively; V π,0 V represents the initial capacity of the energy storage device π; π,T The capacity of the energy storage device π during the termination period.

[0046] Furthermore, the specific process of step (5) is as follows:

[0047] (501) Establish the power balance constraints in the first stage of the two-stage robust optimization model of the integrated energy system of islands and ships.

[0048]

[0049]

[0050] In the formula, The maximum installed capacity of energy equipment θ on island i; The maximum installed capacity of the ship's energy equipment γ; These represent the hydrogen charging power and hydrogen discharging power of the hydrogen storage system on island i at time t during normal operation; H Load,BV This serves as a benchmark value for ship hydrogen energy load. Hydrogen load for ships; The hydrogen consumption power of fuel cell i on island at time t during normal operation; The hydrogen consumption power of the gas turbine on island i at time t during normal operation; P represents the hydrogen energy interaction power between island i and the ship at time t during normal operation. t NOPE,VFCEP represents the discharge power of the fuel cell on the ship at time t during normal operation. t NOPE,VESS,CH P t NOPE,VESS,DIS These represent the charging power and discharging power of the ship's energy storage system at time t during normal operation, respectively. The power transmitted from island i to the ship during energy interaction at time t; P AVL,BV P is the baseline value for ship auxiliary loads. t AVL P represents the ship's auxiliary load at time t; t NOPE,PS P represents the operating power of the ship's propulsion system at time t during normal operation. t NOPE,VFCH The hydrogen consumption power of the fuel cell on island i at time t during normal operation; The effective output of the wind turbine on island i at time t during normal operation; P represents the power transferred from the wind turbine on island i to the electrolyzer at time t. t WT P represents the output of a single wind turbine unit on island i at time t. i WT,max Let t be the planned capacity of the wind turbines on island i at time t; The hydrogen release power of the electrolyzer on island i at time t during normal operation; The discharge power of the fuel cell on island i at time t during normal operation; E represents the discharge power of the gas turbine on island i at time t during normal operation. Load,BV This is the baseline value for ship electrical load; Let t be the electrical load of island i at time t;

[0051] (502) Establish the power constraint balance condition in the second stage of the two-stage robust optimization model of the integrated energy system of islands and ships.

[0052]

[0053]

[0054] In the formula, The hydrogen charging and discharging power of the hydrogen storage system on island i at time t during the emergency operation phase; P represents the hydrogen consumption power of the fuel cell on island i at time t during the emergency operation phase. t EOPE,HEC The hydrogen energy interaction power between island i and the ship at time t during the emergency operation phase; The hydrogen consumption power of the gas turbine on island i at time t during the emergency operation phase; This represents the output power of the wind turbine on island i at time t during the emergency operation phase. The electrical interaction power between island i and the ship at time t during the emergency operation phase; The electrical power transmitted from the wind turbine on island i to the electrolyzer at time t during the emergency operation phase; The hydrogen power transmitted from the wind turbine on island i to the electrolyzer at time t during the emergency operation phase; To ensure the effective output of the wind turbines on island i at time t during the emergency operation phase; The discharge power of the fuel cell on island i at time t during the emergency operation phase; ζ represents the discharge power of the gas turbine on island i at time t during the emergency operation phase; EOPE,HSS,min ζ EOPE,HSS,max These are the minimum and maximum operating coefficients of the hydrogen storage system during the emergency operation phase, respectively; P i HSS,max This represents the upper limit of the power of the hydrogen storage system on island i; σ represents the capacity of the hydrogen storage system on island i at time T during the emergency operation phase; HSS The self-leakage rate of the hydrogen storage system; These represent the hydrogen charging and discharging power of the hydrogen storage system on island i at time T during the emergency operation phase; σ HSS,CH σ HSS,DIS These are the hydrogen charging efficiency and hydrogen discharging efficiency of the hydrogen storage system, respectively.

[0055] Furthermore, in step (6), a nested column constraint generation algorithm is used to solve the two-stage robust optimization model of the island and ship integrated energy system to obtain the joint optimization scheme of the island and ship integrated energy system.

[0056] Beneficial effects: Compared with the prior art, the technical solution of the present invention has the following beneficial technical effects:

[0057] (1) This invention considers the participation of the ship's integrated energy system as a mobile energy storage device in the scheduling of the island integrated energy system. When the ship is anchored between the energy island and the load island, the energy interaction system enables energy interconnection between the ship and the island: it can both supply energy to the island and store excess electrical energy in the island system. Other equipment in the ship's integrated energy system, such as fuel cells and energy storage devices, can also participate in the island's power supply during this period, realizing the spatiotemporal optimization and sharing of energy between islands, thereby improving the collaborative scheduling capability and overall operating efficiency of multiple islands and the ship's integrated energy system. Through a two-stage robust optimization model of the island and the ship's integrated energy system, the ship's integrated energy system can adjust its navigation plan, effectively avoiding energy shortages caused by navigation delays. Through the flexible transfer of hydrogen energy, the reliability of power supply under extreme scenarios can be improved, which not only enhances the renewable energy absorption capacity of the integrated energy system, but also significantly reduces investment and operating costs.

[0058] (2) This invention considers the synergistic optimization of the island and ship integrated energy system. By adjusting the ship's navigation plan and equipment capacity, the impact of equipment failure and ship navigation delay on the integrated energy system is reduced. At the same time, the effect of ocean wave fluctuation on ship navigation is considered to further improve the reliability of ship navigation, thereby enhancing the overall robustness of the island and ship integrated energy system. Attached Figure Description

[0059] Figure 1 This is a flowchart of the method of the present invention. Detailed Implementation

[0060] The technical solution of the present invention will be described in detail below with reference to the accompanying drawings.

[0061] like Figure 1 As shown, this invention proposes a two-stage robust optimization method for integrated energy systems of islands and ships, comprising the following steps:

[0062] (1) Establish the objective function of a two-stage robust optimization model for integrated energy systems of islands and ships;

[0063] (2) Establish the uncertainty set of a two-stage robust optimization model for integrated energy systems of islands and ships;

[0064] (3) Establish the constraints of the integrated energy system of the ship in the two-stage robust optimization model of the island and ship integrated energy system, including the navigation scheduling, propulsion system and energy interaction constraints of the integrated energy system of the ship;

[0065] (4) Establish constraints on energy equipment in the two-stage robust optimization model of the integrated energy system of islands and ships, including constraints on fuel cells, electrolyzers, gas turbines and energy storage equipment;

[0066] (5) Establish the power balance constraints for each stage in the two-stage robust optimization model of the island and ship integrated energy system, including the power balance constraints for the first stage and the second stage.

[0067] (6) The nested column constraint generation algorithm is used to solve the two-stage robust optimization model of the island and ship integrated energy system, and the joint optimization scheme of the island and ship integrated energy system is obtained.

[0068] Using a combined energy system of an island and a ship as an example, a numerical test was conducted on the GAMS optimization software, and the constructed mixed integer linear programming problem was solved using the CPLEX solver.

[0069] To illustrate the advantages of the two-stage robust optimization method for integrated island and ship energy systems proposed in this invention, the impact of different optimization methods for integrated island and ship energy systems on cost, curtailment rate, and hydrogen curtailment rate is explored through comparison.

[0070] Option 1: An optimization model for integrated energy systems of islands and ships that does not consider the worst-case operating scenarios;

[0071] Option 2: Resilience optimization model for integrated energy systems of islands and ships considering only equipment failures;

[0072] Option 3: The proposed two-stage robust resilience optimization model for integrated energy systems of islands and ships.

[0073] The comparison results of the three schemes are shown in Table 1. Compared with Method 1, the total cost of Method 2 and Method 3 increased by 8.1% and 8.2% respectively under ideal conditions. However, due to the effective handling of failure uncertainty by the robust optimization method, the total cost and load loss under the worst-case scenario are both lower than those of Model 1. Further comparison of Method 2 and Method 3 shows that the total cost of Method 2 under ideal conditions is similar to that of Method 3, but there is still a risk of load loss under the worst-case scenario. This is because the speed reduction caused by random fluctuations in sea waves is not considered, which prevents energy from being delivered to the load island as planned. The above results indicate that considering the uncertainty of ship sailing time does not significantly increase the normal operating cost of the system and can effectively reduce the risk of load loss under extreme scenarios, thereby improving the overall economy and stability of the system.

[0074] Table 1 Comparison of Scheme Results

[0075] plan <![CDATA[Cost(10 6 $)]]> Curtailed electricity (MW) <![CDATA[hydrogen discard (10 3 m 3 )]]> Option 1 9.341 67.316 0 Option 2 5.087 0 0.743 Option 3 4.911 0.064 0

[0076] The embodiments are merely illustrative of the technical concept of the present invention and should not be construed as limiting the scope of protection of the present invention. Any modifications made to the technical solution based on the technical concept proposed in this invention shall fall within the scope of protection of this invention.

Claims

1. A two-stage robust optimization method for integrated energy systems of islands and ships, characterized in that, The method includes the following steps: (1) Establish the objective function of a two-stage robust optimization model for integrated energy systems of islands and ships; (2) Establish the uncertainty set of a two-stage robust optimization model for integrated energy systems of islands and ships; (3) Establish the constraints of the integrated energy system of the ship in the two-stage robust optimization model of the island and ship integrated energy system, including the navigation scheduling, propulsion system and energy interaction constraints of the integrated energy system of the ship; (4) Establish the constraints of energy equipment in the two-stage robust optimization model of the integrated energy system of islands and ships, including the constraints of fuel cells, electrolyzers, gas turbines and energy storage equipment; (5) Establish the power balance constraints for each stage in the two-stage robust optimization model of the island and ship integrated energy system, including the power balance constraints for the first stage and the second stage. (6) The nested column constraint generation algorithm is used to solve the two-stage robust optimization model of the island and ship integrated energy system, and the joint optimization scheme of the island and ship integrated energy system is obtained; In step (1), the objective function of the two-stage robust optimization model for the integrated energy system of the island and the ship is established, as follows: （1） In the formula, For variables in the planning phase; This refers to the depreciation rate; The system's operating lifespan; For the scene; Energy equipment for islands includes wind turbines, electrolyzers, fuel cells, hydrogen storage devices, and gas turbines. It is an island; For energy equipment The unit investment cost; For islands Up equipment The planned capacity; It is a marine energy equipment, which is one of the fuel cells, energy storage devices, and hydrogen storage devices; For ship energy equipment The unit investment cost; For ship energy equipment The planned capacity; A collection of scenarios during normal operation; For the scene Total number of days; Costs incurred during the operational phase; For a specific moment; Normal operation phase Time Island Upper energy equipment The electrical power; , These are the normal operation phases. Time Island Upper energy equipment The charging power and discharging power; Normal operation phase Ship energy equipment The electrical power; , These are the normal operation phases. Shipboard energy equipment The charging power and discharging power; The interval time; For one time period; The unit cost of wind curtailment for wind turbine units; Normal operation phase Time Island Wind curtailment power of upstream wind turbines; The uncertainty set for a two-stage robust optimization model of integrated energy systems for islands and ships; Variables for the emergency operation phase; This is a set of scenarios for the emergency operation phase, which takes into account the uncertainty of equipment failure and ship sailing time. Emergency Operation Phase Time Island Upper energy equipment The electrical power; , They are divided into emergency operation phases. Time Island Upper energy equipment The charging power and discharging power; Emergency Operation Phase Ship energy equipment The electrical power; , They are divided into emergency operation phases. Ship energy equipment The charging power and discharging power; Emergency Operation Phase Time Island Wind curtailment power of upstream wind turbines; Cost per unit of power load shedding; for Time Island Power load shedding; Cost per unit of hydrogen energy load shedding; for Time Island Hydrogen energy load shedding.

2. The two-stage robust optimization method for an integrated energy system for islands and ships according to claim 1, characterized in that, In step (2), the uncertainty set of the two-stage robust optimization model for the integrated energy system of the island and the ship is established, as follows: (2) (3) (4) In the formula, For a set of uncertain sailing times; For equipment failure uncertainty set; For ships in The actual sailing speed at any given moment; and For calculating the power of the ship's propulsion system; To improve system efficiency; for The output power of the all-electric ship's electric motor at all times; For resistance during still water navigation; Adding resistance to the waves; The beam of the ship; This is the block coefficient of the ship; The length of the ship; The density of seawater; It is the acceleration due to gravity; , These are the actual and predicted values ​​of the wave height, respectively. This is the positive offset of the wave height; This is the negative offset of the wave height; A binary variable, representing in The value is 1 if the wave height deviates positively at any given moment, and 0 otherwise. A binary variable, representing in The value is 1 if the wave height has a negative offset at any given time, and 0 otherwise. For uncertain budgets related to ship navigation; This is a binary variable indicating whether the electrolytic cell is faulty; 1 indicates a fault, and 0 indicates otherwise. This is a binary variable indicating whether the fuel cell is faulty; 1 indicates a fault, and 0 indicates otherwise. This is a binary variable indicating whether the gas turbine is malfunctioning; 1 indicates a malfunction, and 0 indicates otherwise. Uncertain budget for emergency equipment use.

3. The two-stage robust optimization method for an integrated energy system for islands and ships according to claim 2, characterized in that, The specific process of step (3) is as follows: (301) Establishing navigation scheduling constraints for integrated ship energy systems (5) (6) (7) (8) (9) (10) In the formula, This includes all islands except for the starting and ending points; In order to be in Ships relative to islands at any time Position variables, For ships in Always docked on the island superior, Characterizing the ship in Always leaving the island Currently in navigation mode; In order to be in Ships relative to islands at any time Position variables, For ships in Always docked on the island superior, Characterizing the ship in Always leaving the island Currently in navigation mode; This refers to the minimum berthing time for a vessel on an island. , , The arrival of ships at the island , and The latest time; , The ships are respectively in and Always on the island Real-time distance traveled along the route originating from the current point; For ships in The sailing speed at any given moment; For ships in islands The total distance along the route originating from the given point; (302) Establishing propulsion system constraints for integrated ship energy systems (11) (12) In the formula, This indicates the start / stop status of the motor; 1 indicates start, and 0 indicates stop. (303) Establishing energy interaction constraints for integrated ship energy systems (13) (14) In the formula, For islands Power exchanged with ship electrical power; For islands Upper limit of power exchange with ship electrical energy; For islands Hydrogen energy interaction power with ships; For islands Upper limit of the power interaction between hydrogen energy and ships.

4. The two-stage robust optimization method for an integrated energy system for islands and ships according to claim 3, characterized in that, The specific process of step (4) is as follows: (401) Establishing constraints for fuel cells (15) (16) (17) (18) In the formula, For fuel cells in Output power at any given moment; and This is the unit conversion factor; and The fitting coefficients for the operating power of the fuel cell are denoted as . For fuel cells in Load current at any given moment; This represents the number of fuel cell stacks; This is the start / stop state variable for the fuel cell; it is 1 when the fuel cell is running and 0 otherwise. This is the operating state variable of the fuel cell; it is 1 for normal operation and 0 otherwise. for The hydrogen consumption power of the fuel cell at all times; It is Faraday's constant; , These are the minimum and maximum load currents of the fuel cell, respectively. This refers to the rated capacity of the fuel cell; (402) Establishing the constraints of the electrolytic cell (19) (20) (21) (22) In the formula, For the electrolytic cell in Input power at any given time; and The fitting coefficients for the working power of the electrolytic cell; For the electrolytic cell in Load current at any given moment; This refers to the number of electrolysis units; This indicates the start / stop status of the electrolytic cell; 1 indicates start-up, and 0 indicates otherwise. This is the operating status variable of the electrolysis unit; it is 1 for normal operation and 0 otherwise. This refers to the output power of the electrolytic cell; To improve the working efficiency of the electrolytic cell; , These are the minimum and maximum load currents of the electrolytic cell, respectively. This refers to the rated capacity of the electrolytic cell; (403) Establishing constraints for the gas turbine (23) (24) (25) In the formula, This represents the maximum gradeability of the gas turbine. This refers to the rated capacity of the gas turbine. , Gas turbines and Output power at time +1; , These are the minimum and maximum operating ratios of the gas turbine, respectively. This is the start-stop state variable for the gas turbine; it is 1 when the turbine is running and 0 otherwise. This is the operating state variable of the gas turbine; it is 1 for normal operation and 0 otherwise. This refers to the hydrogen consumption power of the gas turbine. For the working efficiency of the gas turbine; (404) Establish constraints for energy storage devices: (26) (27) (28) (29) (30) In the formula, It is an energy storage device, which is one of the following: electrical storage devices and hydrogen storage devices. , Energy storage devices at times t and t+1 respectively The amount of self-discharge; For energy storage devices The self-discharge rate; , These represent the charging power and discharging power of the energy storage device at time t, respectively. , Energy storage devices The charging efficiency and discharging efficiency; For energy storage devices Maximum gradeability; Energy storage devices Rated capacity; , Energy storage devices The minimum and maximum operating coefficients; Energy storage device for the initial period The capacity; Energy storage devices for the end of the period The capacity.

5. The two-stage robust optimization method for an integrated energy system for islands and ships according to claim 4, characterized in that, The specific process of step (5) is as follows: (501) Establish the power balance constraints in the first stage of the two-stage robust optimization model for integrated energy systems of islands and ships. (31) (32) (33) (34) (35) (36) (37) (38) (39) In the formula, For islands Upper energy equipment Maximum installed capacity; Energy equipment for ships Maximum installed capacity; , These are the normal operation phases. Time Island The hydrogen charging and discharging power of the hydrogen storage system; This serves as a benchmark value for ship hydrogen energy load. Hydrogen load for ships; Normal operation phase Time Island Fuel cell hydrogen consumption power; Normal operation phase Time Island Hydrogen power consumption of gas turbine; Normal operation phase Time Island Hydrogen energy interaction power with ships; Normal operation phase The discharge power of fuel cells on the ship at any given time; , These are the normal operation phases. The charging and discharging power of the ship's energy storage system at all times; Time Island Power transmitted to the ship during energy exchange; This is the baseline value for ship auxiliary loads; for Ship auxiliary loads at all times; Normal operation phase The operating power of the ship's propulsion system at all times; Normal operation phase Time Island The hydrogen consumption power of the fuel cell; Normal operation phase Time Island Effective output of the wind turbine unit; for Time Island The power transmitted from the wind turbine to the electrolytic cell; for Time Island The output of the wind turbine unit in the upper unit; for Time Island The planned capacity of the wind turbine units; Normal operation phase Time Island Hydrogen release power of the upper electrolyzer; Normal operation phase Time Island Upgrade fuel cell discharge power; Normal operation phase Time Island Gas turbine discharge power; This is the baseline value for ship electrical load; for Time Island Electrical load; (502) Establish the power constraint balance condition in the second stage of the two-stage robust optimization model of the integrated energy system of islands and ships. (40) (41) (42) (43) (44) (45) (46) (47) In the formula, , Emergency Operation Phase Time Island The hydrogen charging and discharging power of the hydrogen storage system; Emergency Operation Phase Time Island The hydrogen consumption power of the fuel cell; Emergency Operation Phase Time Island Hydrogen energy interaction power with ships; Emergency Operation Phase Time Island Hydrogen power consumption of gas turbine; Emergency Operation Phase Time Island The output power of the wind turbine unit; Emergency Operation Phase Time Island Electrical power exchange with the ship; Emergency Operation Phase Time Island The electrical power transmitted from the wind turbine to the electrolytic cell; Emergency Operation Phase Time Island Hydrogen power transferred from the wind turbine to the electrolyzer; Emergency Operation Phase Time Island Effective output of the wind turbine unit; Emergency Operation Phase Time Island The discharge power of the fuel cell; Emergency Operation Phase Time Island The discharge power of the gas turbine; , These are the minimum and maximum operating coefficients of the hydrogen storage system during the emergency operation phase, respectively. For islands The upper power limit of the hydrogen storage system; Emergency Operation Phase Time Island The capacity of the hydrogen storage system; The self-leakage rate of the hydrogen storage system; , They are divided into emergency operation phases. Time Island The hydrogen charging and discharging power of the hydrogen storage system; , These are the hydrogen charging efficiency and hydrogen discharging efficiency of the hydrogen storage system, respectively.

Images