A method for enhancing the resilience of integrated energy systems in seaports considering multi-energy network failures

Abstract

This paper discloses a method for enhancing the resilience of a port integrated energy system considering multi-energy network failures. The method includes the following steps: 1) establishing the objective function of the port integrated energy system resilience enhancement model; 2) establishing multi-energy network failure constraints for the port integrated energy system resilience enhancement model; 3) establishing combined cooling, heating, and power (CCHP) constraints for the port integrated energy system resilience enhancement model; 4) establishing container transportation constraints for the port integrated energy system resilience enhancement model; 5) solving the port integrated energy system resilience enhancement model to obtain a resilient operation strategy for the port integrated energy system. This invention considers the coupling characteristics of the multi-energy network (electricity, gas, heat, and cooling) of the port integrated energy system during extreme events. Through collaborative optimization of the multi-energy network and the port logistics system, it improves the flexibility, reliability, and resilience of the port integrated energy system.

Description

Technical Field

[0001] This invention belongs to the field of integrated energy technology for seaports, and specifically relates to a method for improving the resilience of integrated energy systems for seaports that takes into account multi-energy network failures. Background Technology

[0002] Against the backdrop of continuous expansion in global trade, the energy supply and demand system of seaports, as international logistics hubs, is exhibiting multi-dimensional and complex characteristics, encompassing a complex supply network involving the coupling of multiple energy flows such as electricity, gas, heat, and cooling. As core nodes of the global supply chain and regional energy distribution centers, modern seaport energy systems must simultaneously meet high concurrent load demands, the coupling characteristics of multiple energy sources, and the requirements for all-day operational stability. They also face multi-dimensional risks and challenges, including extreme weather conditions, aging infrastructure, and cybersecurity threats. The structural vulnerabilities of the system easily trigger cross-media propagation effects of multi-energy network failures, forming cascading failures through energy flow-information flow coupling channels. This can not only lead to the collapse of critical energy nodes and the paralysis of supply chain operations, but also threaten the regional energy security architecture, resulting in significant economic losses and social impacts. This risk is particularly prominent in the context of carbon neutrality transition, as traditional single-energy dispatch models are no longer sufficient to balance the multi-dimensional contradictions between improving system resilience, achieving low-carbon goals, and ensuring seaport operational efficiency.

[0003] Currently, research on improving the resilience of integrated energy systems mainly focuses on traditional power systems or urban area energy networks, while insufficient attention is paid to the unique application scenario of seaports. Multi-energy networks are tightly coupled through energy conversion equipment, and a single network failure can easily trigger a cascading failure across networks through energy flow transmission. Existing research lacks comprehensive modeling and suppression strategies for dynamic fault propagation paths. Therefore, improving the resilience of integrated energy systems in seaports that consider multi-energy network failures has become a pressing technical challenge in this field. Summary of the Invention

[0004] The purpose of this invention is to address the aforementioned problems by proposing a method for enhancing the resilience of integrated energy systems in seaports that consider multi-energy network failures, thereby improving the resilience of integrated energy systems in seaports that consider multi-energy network failures and solving a technical problem that urgently needs to be overcome in the current field.

[0005] Technical Solution: To address the aforementioned technical issues, this invention proposes a method for enhancing the resilience of a port integrated energy system considering multi-energy network failures. This method includes the following steps:

[0006] Step 1: Establish the objective function of the resilience enhancement model for the integrated energy system of the seaport;

[0007] Step 2: Establish multi-energy network fault constraints for the resilience enhancement model of the integrated energy system of the seaport;

[0008] Step 3: Establish combined cooling, heating and power (CCHP) constraints for the resilience enhancement model of the integrated energy system of the seaport;

[0009] Step 4: Establish container transport constraints for the resilience enhancement model of the integrated energy system of the seaport;

[0010] Step 5: Based on the objective function of Step 1 and the constraints of Steps 2-4, solve the resilience enhancement model of the integrated energy system of the seaport to obtain the resilience operation strategy of the integrated energy system of the seaport.

[0011] Furthermore, in step (1), the objective function for establishing the resilience enhancement model of the integrated energy system of the seaport is expressed as follows:

[0012]

[0013] In the formula, a represents a combined cooling, heating and power (CCHP) unit; t represents a time period; s represents a ship; e represents a grid node; n represents a new energy unit; d represents an electrolytic cell; b represents a boiler; p represents a heat pump; r represents an electric chiller; g represents a gas grid node; h represents a heating grid node; and c represents a cooling grid node. This represents the fixed cost of combined cooling, heating and power (CCHP) unit a; This represents the startup cost of combined cooling, heating and power (CCHP) unit a; This represents the shutdown cost of combined cooling, heating and power (CCHP) unit a; This indicates whether the combined cooling, heating and power (CCHP) unit a is fixed during time period t; if yes, set to 1, otherwise set to 0. This indicates whether the combined cooling, heating and power (CCHP) unit a is started during time period t; if yes, set to 1, otherwise set to 0. This indicates whether the combined cooling, heating and power (CCHP) unit a is off during time period t; set to 1 if yes, otherwise set to 0. This represents the power generation cost of combined cooling, heating and power (CCHP) unit a; This represents the active power output of the combined cooling, heating and power (CCHP) unit a during time period t. This represents the operation and maintenance cost of combined cooling, heating and power (CCHP) unit a; This indicates the penalty for late departure of vessel s; Indicates the actual departure time of vessel s; Indicates the expected departure time of vessel s; This represents the electricity price during time period t; This represents the active power injected by the upstream power grid into grid node e during time period t; This represents the operation and maintenance cost of the new energy unit n; This represents the active power output of the new energy unit n during time period t; This represents the operation and maintenance cost of electrolytic cell d; This represents the active power consumption of electrolytic cell d during time period t; This represents the operation and maintenance cost of boiler b; This represents the gas consumption of boiler b during time period t; This represents the operation and maintenance cost of heat pump p; This represents the active power consumption of the heat pump p during time period t; This represents the maintenance cost of the electric chiller r; This represents the active power consumption of the electric chiller r during time period t; Indicates the cost of electrical load loss; Indicates the cost of gas loss of load; Indicates the cost of thermal load loss; Indicates the cost of cooling load loss; This represents the load shedding rate of grid node e during time period t; This represents the total active load of grid node e during time period t; This represents the gas load loss rate of gas network node g during time period t; This represents the total gas load of gas network node g during time period t; This represents the heat load loss rate of heating network node h during time period t; This represents the total heat load of heat network node h during time period t; This represents the cooling load loss rate of cold network node c during time period t; This represents the total cooling load of cold network node c during time period t.

[0014] Furthermore, the specific process of step (2) is as follows:

[0015] (201) Establish grid fault constraints:

[0016]

[0017] In the formula, f represents a power grid node; o represents a berth; q represents a quay crane; y represents a yard crane; v represents an automated guided container truck; and ef represents the power grid branch connecting power grid nodes e and f. This represents the set of combined cooling, heating and power (CCHP) units connected to grid node e; This represents the set of new energy generating units connected to grid node e; This represents the set of electrolytic cells connected to grid node e; This represents the set of heat pumps connected to grid node e; This represents the set of electric chillers connected to the power grid node e; Represents the set of power grid branches with terminal node e; Represents the set of power grid branches with the first node being e; This represents the set of berths connected to the power grid node e; This represents the set of quay bridges connected to the power grid node e; This represents the set of field bridges connected to grid node e; This represents the set of autonomously navigated container trucks connected to grid node e; Let f represent the set of power grid branches with terminal node f; This represents the active power transmitted by the power grid branch ef during time period t; This represents the active power consumption of berth o during time period t; This represents the active power consumption of the quay crane q during time period t; This represents the active power consumption of the field bridge y during time period t; This represents the active power that the automated guided container truck v receives from the power grid during time period t; This represents the reactive power injected by the upstream power grid into grid node e during time period t; This indicates the reactive power output of the combined cooling, heating and power (CCHP) unit a during time period t; This represents the reactive power output of the new energy unit n during time period t; This represents the reactive power consumption of electrolytic cell d during time period t; This represents the reactive power consumption of the heat pump p during time period t; This represents the reactive power consumption of the electric chiller r during time period t; This represents the reactive power transmitted by the power grid branch ef during time period t; This represents the reactive power consumption of berth o during time period t; This represents the reactive power consumption of the quay crane q during time period t; This represents the reactive power consumption of the field bridge y during time period t; This represents the reactive power that the automated guided container truck v receives from the power grid during time period t; This represents the total reactive load of grid node e during time period t; This represents a very large constant; This indicates whether the power grid branch ef is faulty during time period t; if yes, set to 0, otherwise set to 1. This represents the voltage amplitude at grid node e during time period t; This represents the voltage amplitude at grid node f during time period t; This represents the resistance of branch ef; V represents the reactance of branch ef; EB,0 Indicates the reference voltage;

[0018] (202) Establish gas network fault constraints:

[0019]

[0020] In the formula, k represents a gas network node; l represents a gas compressor; gk represents the gas network pipeline connecting gas network nodes g and k; This represents the set of gas pipelines with the first node g. This represents the set of gas compressors connected to gas network node g. This represents the set of combined cooling, heating and power (CCHP) units connected to the gas network node g. This represents the set of electrolytic cells connected to the gas network node g; This represents the set of boilers connected to the gas network node g; σ represents the hydrogen gas transported by the gas pipeline gk during time period t; l This represents the compression coefficient of the gas compressor l; This represents the hydrogen output by gas compressor l during time period t; This represents the hydrogen consumption of combined cooling, heating and power unit a during time period t; This represents the conversion efficiency of electrolytic cell d during time period t; This indicates the gas transmitted by the gas pipeline gk during time period t; This represents the gas output of gas compressor l during time period t; This represents the amount of gas injected into gas network node g during time period t; This represents the gas consumption of a combined cooling, heating, and power (CCHP) unit a during time period t; κ GP,H Indicates the hydrogen blending ratio in the gas pipeline network; This indicates whether the gas pipeline gk is faulty during time period t; if yes, set to 0, otherwise set to 1. This indicates the maximum gas flow rate in the gas pipeline gk during time period t; This indicates the maximum hydrogen flow rate in the gas pipeline gk during time period t;

[0021] (203) Establish heating network fault constraints:

[0022]

[0023] In the formula, i represents a heat network node; hi represents the heat network pipe connecting heat network nodes h and i; This represents the set of combined cooling, heating and power (CCHP) units connected to the heating network node h; This represents the set of boilers connected to the heating network node h; This represents the set of heat pumps connected to the heat network node h. This represents the set of heating network pipes with terminal node h; This represents the set of heat network pipelines with the first node being h; This indicates the heat output of the combined cooling, heating and power (CCHP) unit a during time period t; This indicates the heat conversion efficiency of boiler b; ρ represents the heat conversion efficiency of the heat pump p; W This indicates the specific heat capacity of water; This represents the heat source mass flow rate of heat network node h during time period t; This represents the water supply temperature at node h in the heating network during time period t; This represents the return water temperature of heating network node h during time period t; This represents the heat load mass flow rate of heat network node h during time period t; This indicates the outlet temperature of the heating network pipe hi during time period t; This indicates whether the heating network pipeline hi is faulty during time period t; set to 0 if yes, otherwise set to 1. This represents the ambient temperature during time period t; This indicates the inlet temperature of the heating network pipe hi during time period t; This represents the heat transfer coefficient of the heating network pipe hi; Indicates the length of the heating network pipe hi; This represents the hot water mass flow rate of the heating network pipe hi during time period t; This represents the mixing temperature of heat network node h during time period t;

[0024] (204) Establish cold network fault constraints:

[0025]

[0026]

[0027] In the formula, x represents a cold network node; cx represents the cold network pipe connecting cold network nodes c and x; This represents the set of combined cooling, heating, and power (CCHP) units connected to the cold network node c; This represents the set of electric chillers connected to node c of the cold network; Represents the set of cold network pipes with terminal node c; This represents the set of cold network pipes with the first end node being c; This represents the cooling output of the combined cooling, heating and power (CCHP) unit a during time period t. This indicates the refrigeration conversion efficiency of the electric chiller r; This represents the cold source mass flow rate of cold network node c during time period t; This represents the return water temperature of cold network node c during time period t; This represents the water supply temperature of cold network node c during time period t; This represents the mass flow rate of the cooling load at node c in the cold network during time period t. This indicates whether the cold network pipe CX is faulty during time period t; if yes, set to 0, otherwise set to 1. This represents the outlet temperature of the cold network pipe cx during time period t; This represents the inlet temperature of the cold network pipe cx during time period t; This represents the heat transfer coefficient of the cold network pipe cx; This indicates the length of the cold network pipe cx; This represents the mass flow rate of chilled water in the cold network pipe cx during time period t; This represents the mixing temperature of cold network node c during time period t.

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

[0029] (301) Establishing the combined cooling, heating and power (CCHP) constraints for the resilience enhancement model of a seaport integrated energy system:

[0030]

[0031] In the formula, This indicates the power generation efficiency of combined cooling, heating and power unit a; This indicates the heat production efficiency of combined cooling, heating and power unit a; This indicates the maximum downward ramp rate of the combined cooling, heating and power (CCHP) unit a; This indicates the maximum upward ramp rate of the combined cooling, heating and power (CCHP) unit a; This indicates the maximum downhill ramp rate when the combined cooling, heating and power (CCHP) unit a is shut down. This indicates the maximum uphill ramp rate at startup for combined cooling, heating and power (CCHP) unit a; This represents the active power output of the combined cooling, heating and power (CCHP) unit a during the time period t-1; This indicates whether the combined cooling, heating and power (CCHP) unit a is fixed during the t-1 time period; if yes, set it to 1, otherwise set it to 0. This represents the minimum active power output of the combined cooling, heating and power (CCHP) unit a; This indicates the maximum active power output of the combined cooling, heating and power (CCHP) unit a.

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

[0033] (401) Establishing container transport constraints for a resilience enhancement model of a seaport integrated energy system:

[0034]

[0035] In the formula: This represents the number of containers loaded and unloaded by quay crane q from ship s during time period t; This indicates whether the quay crane q serves vessel s during time period t; if yes, it is set to 1, otherwise it is set to 0. This indicates the number of containers loaded and unloaded from ship s by the yard crane y during time period t; This indicates whether the field bridge y serves vessel s during time period t; if yes, set to 1, otherwise set to 0. This represents the number of containers loaded and unloaded by the automated guided vehicle (AGV) from the ship (S) during time period t. This indicates whether the automated navigation container truck v serves the vessel s during time period t; set to 1 if yes, otherwise set to 0. This indicates the number of containers loaded on ship s; This represents the electrical energy storage of the automated guided container truck v during time period t; This represents the electrical energy storage of the automated guided container truck v during the time period t-1; This indicates the charging efficiency of the automated guided container truck v. This indicates the discharge efficiency of the automatically guided container truck v. This represents the active power consumption of the automated guided container truck v during time period t.

[0036] Furthermore, in step (5), a resilience enhancement model for the integrated energy system of a port with multi-energy network failure is written in GAMS, and the resilience operation strategy of the integrated energy system of the port is obtained by solving the model.

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

[0038] This invention establishes a dynamic fault constraint model for a multi-energy network encompassing electricity, gas, heat, and cooling, accurately characterizing fault propagation paths and energy coupling characteristics to effectively suppress the risk of cross-energy cascading failures. Simultaneously, it dynamically correlates the energy consumption of container transport equipment with grid load, coordinating the charging and discharging strategies of quay cranes, yard cranes, and automated guided vehicles (AGVs) to ensure port operational efficiency and energy supply-demand balance. By solving for and generating a multi-objective resilient operation strategy that considers the loss-of-load costs, maintenance costs, and late vessel return penalties of various energy sources, it addresses complex scenarios such as extreme weather and equipment aging, reducing economic losses. Ultimately, it achieves a comprehensive improvement in the safety, economy, and sustainability of the port's integrated energy system, providing systematic support for the green and efficient operation of the port. Attached Figure Description

[0039] Figure 1 A flowchart illustrating a method for enhancing the resilience of a port integrated energy system in the event of multi-energy network failures;

[0040] Figure 2 This is a cost result diagram for a multi-energy network integrated energy system in a seaport. Detailed Implementation

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

[0042] like Figure 1 As shown, this invention proposes a method for enhancing the resilience of a port integrated energy system considering multi-energy network failures. This method includes the following steps:

[0043] Step 1: Establish the objective function of the resilience enhancement model for the integrated energy system of the seaport;

[0044] Step 2: Establish multi-energy network fault constraints for the resilience enhancement model of the integrated energy system of the seaport;

[0045] Step 3: Establish combined cooling, heating and power (CCHP) constraints for the resilience enhancement model of the integrated energy system of the seaport;

[0046] Step 4: Establish container transport constraints for the resilience enhancement model of the integrated energy system of the seaport;

[0047] Step 5: Based on the objective function of Step 1 and the constraints of Steps 2-4, solve the resilience enhancement model of the integrated energy system of the seaport to obtain the resilience operation strategy of the integrated energy system of the seaport.

[0048] Furthermore, in step (1), the objective function for establishing the resilience enhancement model of the integrated energy system of the seaport is expressed as follows:

[0049]

[0050] In the formula, a represents a combined cooling, heating and power (CCHP) unit; t represents a time period; s represents a ship; e represents a grid node; n represents a new energy unit; d represents an electrolytic cell; b represents a boiler; p represents a heat pump; r represents an electric chiller; g represents a gas grid node; h represents a heating grid node; and c represents a cooling grid node. This represents the fixed cost of combined cooling, heating and power (CCHP) unit a; This represents the startup cost of combined cooling, heating and power (CCHP) unit a; This represents the shutdown cost of combined cooling, heating and power (CCHP) unit a; This indicates whether the combined cooling, heating and power (CCHP) unit a is fixed during time period t; if yes, set to 1, otherwise set to 0. This indicates whether the combined cooling, heating and power (CCHP) unit a is started during time period t; if yes, set to 1, otherwise set to 0. This indicates whether the combined cooling, heating and power (CCHP) unit a is off during time period t; set to 1 if yes, otherwise set to 0. This represents the power generation cost of combined cooling, heating and power (CCHP) unit a; This represents the active power output of the combined cooling, heating and power (CCHP) unit a during time period t. This represents the operation and maintenance cost of combined cooling, heating and power (CCHP) unit a; This indicates the penalty for late departure of vessel s; Indicates the actual departure time of vessel s; Indicates the expected departure time of vessel s; This represents the electricity price during time period t; This represents the active power injected by the upstream power grid into grid node e during time period t; This represents the operation and maintenance cost of the new energy unit n; This represents the active power output of the new energy unit n during time period t; This represents the operation and maintenance cost of electrolytic cell d; This represents the active power consumption of electrolytic cell d during time period t; This represents the operation and maintenance cost of boiler b; This represents the gas consumption of boiler b during time period t; This represents the operation and maintenance cost of heat pump p; This represents the active power consumption of the heat pump p during time period t; This represents the maintenance cost of the electric chiller r; This represents the active power consumption of the electric chiller r during time period t; Indicates the cost of electrical load loss; Indicates the cost of gas loss of load; Indicates the cost of thermal load loss; Indicates the cost of cooling load loss; This represents the load shedding rate of grid node e during time period t; This represents the total active load of grid node e during time period t; This represents the gas load loss rate of gas network node g during time period t; This represents the total gas load of gas network node g during time period t; This represents the heat load loss rate of heating network node h during time period t; This represents the total heat load of heat network node h during time period t; This represents the cooling load loss rate of cold network node c during time period t; This represents the total cooling load of cold network node c during time period t.

[0051] Furthermore, the specific process of step (2) is as follows:

[0052] (201) Establish grid fault constraints:

[0053]

[0054] In the formula, f represents a power grid node; o represents a berth; q represents a quay crane; y represents a yard crane; v represents an automated guided container truck; and ef represents the power grid branch connecting power grid nodes e and f. This represents the set of combined cooling, heating and power (CCHP) units connected to grid node e; This represents the set of new energy generating units connected to grid node e; This represents the set of electrolytic cells connected to grid node e; This represents the set of heat pumps connected to grid node e; This represents the set of electric chillers connected to the power grid node e; Represents the set of power grid branches with terminal node e; Represents the set of power grid branches with the first node being e; This represents the set of berths connected to the power grid node e; This represents the set of quay bridges connected to the power grid node e; This represents the set of field bridges connected to grid node e; This represents the set of autonomously navigated container trucks connected to grid node e; Let f represent the set of power grid branches with terminal node f; This represents the active power transmitted by the power grid branch ef during time period t; This represents the active power consumption of berth o during time period t; This represents the active power consumption of the quay crane q during time period t; This represents the active power consumption of the field bridge y during time period t; This represents the active power that the automated guided container truck v receives from the power grid during time period t; This represents the reactive power injected by the upstream power grid into grid node e during time period t; This indicates the reactive power output of the combined cooling, heating and power (CCHP) unit a during time period t; This represents the reactive power output of the new energy unit n during time period t; This represents the reactive power consumption of electrolytic cell d during time period t; This represents the reactive power consumption of the heat pump p during time period t; This represents the reactive power consumption of the electric chiller r during time period t; This represents the reactive power transmitted by the power grid branch ef during time period t; This represents the reactive power consumption of berth o during time period t; This represents the reactive power consumption of the quay crane q during time period t; This represents the reactive power consumption of the field bridge y during time period t; This represents the reactive power that the automated guided container truck v receives from the power grid during time period t; This represents the total reactive load of grid node e during time period t; This represents a very large constant; This indicates whether the power grid branch ef is faulty during time period t; if yes, set to 0, otherwise set to 1. This represents the voltage amplitude at grid node e during time period t; This represents the voltage amplitude at grid node f during time period t; This represents the resistance of branch ef; V represents the reactance of branch ef; EB,0 Indicates the reference voltage;

[0055] (202) Establish gas network fault constraints:

[0056]

[0057] In the formula, k represents a gas network node; l represents a gas compressor; gk represents the gas network pipeline connecting gas network nodes g and k; This represents the set of gas pipelines with the first node g. This represents the set of gas compressors connected to gas network node g. This represents the set of combined cooling, heating and power (CCHP) units connected to the gas network node g. This represents the set of electrolytic cells connected to the gas network node g; This represents the set of boilers connected to the gas network node g; σ represents the hydrogen gas transported by the gas pipeline gk during time period t; l This represents the compression coefficient of the gas compressor l; This represents the hydrogen output by gas compressor l during time period t; This represents the hydrogen consumption of combined cooling, heating and power unit a during time period t; This represents the conversion efficiency of electrolytic cell d during time period t; This indicates the gas transmitted by the gas pipeline gk during time period t; This represents the gas output of gas compressor l during time period t; This represents the amount of gas injected into gas network node g during time period t; This represents the gas consumption of a combined cooling, heating, and power (CCHP) unit a during time period t; κ GP,H Indicates the hydrogen blending ratio in the gas pipeline network; This indicates whether the gas pipeline gk is faulty during time period t; if yes, set to 0, otherwise set to 1. This indicates the maximum gas flow rate in the gas pipeline gk during time period t; This indicates the maximum hydrogen flow rate in the gas pipeline gk during time period t;

[0058] (203) Establish heating network fault constraints:

[0059]

[0060] In the formula, i represents a heat network node; hi represents the heat network pipe connecting heat network nodes h and i; This represents the set of combined cooling, heating and power (CCHP) units connected to the heating network node h; This represents the set of boilers connected to the heating network node h; This represents the set of heat pumps connected to the heat network node h. This represents the set of heating network pipes with terminal node h; This represents the set of heat network pipelines with the first node being h; This indicates the heat output of the combined cooling, heating and power (CCHP) unit a during time period t; This indicates the heat conversion efficiency of boiler b; ρ represents the heat conversion efficiency of the heat pump p; W This indicates the specific heat capacity of water; This represents the heat source mass flow rate of heat network node h during time period t; This represents the water supply temperature at node h in the heating network during time period t; This represents the return water temperature of heating network node h during time period t; This represents the heat load mass flow rate of heat network node h during time period t; This indicates the outlet temperature of the heating network pipe hi during time period t; This indicates whether the heating network pipeline hi is faulty during time period t; set to 0 if yes, otherwise set to 1. This represents the ambient temperature during time period t; This indicates the inlet temperature of the heating network pipe hi during time period t; This represents the heat transfer coefficient of the heating network pipe hi; Indicates the length of the heating network pipe hi; This represents the hot water mass flow rate of the heating network pipe hi during time period t; This represents the mixing temperature of heat network node h during time period t;

[0061] (204) Establish cold network fault constraints:

[0062]

[0063] In the formula, x represents a cold network node; cx represents the cold network pipe connecting cold network nodes c and x; This represents the set of combined cooling, heating, and power (CCHP) units connected to the cold network node c; This represents the set of electric chillers connected to node c of the cold network; Represents the set of cold network pipes with terminal node c; This represents the set of cold network pipes with the first end node being c; This represents the cooling output of the combined cooling, heating and power (CCHP) unit a during time period t. This indicates the refrigeration conversion efficiency of the electric chiller r; This represents the cold source mass flow rate of cold network node c during time period t; This represents the return water temperature of cold network node c during time period t; This represents the water supply temperature of cold network node c during time period t; This represents the mass flow rate of the cooling load at node c in the cold network during time period t. This indicates whether the cold network pipe CX is faulty during time period t; if yes, set to 0, otherwise set to 1. This represents the outlet temperature of the cold network pipe cx during time period t; This represents the inlet temperature of the cold network pipe cx during time period t; This represents the heat transfer coefficient of the cold network pipe cx; This indicates the length of the cold network pipe cx; This represents the mass flow rate of chilled water in the cold network pipe cx during time period t; This represents the mixing temperature of cold network node c during time period t.

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

[0065] (301) Establishing the combined cooling, heating and power (CCHP) constraints for the resilience enhancement model of a seaport integrated energy system:

[0066]

[0067]

[0068] In the formula, This indicates the power generation efficiency of combined cooling, heating and power unit a; This indicates the heat production efficiency of combined cooling, heating and power unit a; This indicates the maximum downward ramp rate of the combined cooling, heating and power (CCHP) unit a; This indicates the maximum upward ramp rate of the combined cooling, heating and power (CCHP) unit a; This indicates the maximum downhill ramp rate when the combined cooling, heating and power (CCHP) unit a is shut down. This indicates the maximum uphill ramp rate at startup for combined cooling, heating and power (CCHP) unit a; This represents the active power output of the combined cooling, heating and power (CCHP) unit a during the time period t-1; This indicates whether the combined cooling, heating and power (CCHP) unit a is fixed during the t-1 time period; if yes, set it to 1, otherwise set it to 0. This represents the minimum active power output of the combined cooling, heating and power (CCHP) unit a; This indicates the maximum active power output of the combined cooling, heating and power (CCHP) unit a.

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

[0070] (401) Establishing container transport constraints for a resilience enhancement model of a seaport integrated energy system:

[0071]

[0072] In the formula: This represents the number of containers loaded and unloaded by quay crane q from ship s during time period t; This indicates whether the quay crane q serves vessel s during time period t; if yes, it is set to 1, otherwise it is set to 0. This indicates the number of containers loaded and unloaded from ship s by the yard crane y during time period t; This indicates whether the field bridge y serves vessel s during time period t; if yes, set to 1, otherwise set to 0. This represents the number of containers loaded and unloaded by the automated guided vehicle (AGV) from the ship (S) during time period t. This indicates whether the automated navigation container truck v serves the vessel s during time period t; set to 1 if yes, otherwise set to 0. This indicates the number of containers loaded on ship s; This represents the electrical energy storage of the automated guided container truck v during time period t; This represents the electrical energy storage of the automated guided container truck v during the time period t-1; This indicates the charging efficiency of the automated guided container truck v. This indicates the discharge efficiency of the automatically guided container truck v. This represents the active power consumption of the automated guided container truck v during time period t.

[0073] Furthermore, in step (5), a resilience enhancement model for the integrated energy system of a port with multi-energy network failure is written in GAMS, and the resilience operation strategy of the integrated energy system of the port is obtained by solving the model.

[0074] A multi-energy network system of a typical seaport is selected as a case study. Using the ship berthing schedule and container loading / unloading tasks for a specific day in actual operation as input data, a 24-hour continuous energy-transportation joint optimization model is constructed. The simulation is implemented on the GAMS platform, and the mixed-integer linear programming problem in the model is solved by calling the CPLEX solver.

[0075] To illustrate the advantages of the proposed method for enhancing the resilience of a port integrated energy system in the event of multi-energy network failures, the impact on total cost is explored by employing the following different schemes.

[0076] Option 1: Ignore multi-energy network failures in the integrated energy system of the seaport;

[0077] Option 2: Only consider power network failures in the integrated energy system of the seaport;

[0078] Option 3: Consider the multi-energy network failure of the integrated energy system in the seaport;

[0079] The scheduling results of the three schemes are shown in Table 1. Scheme 1 does not consider multi-energy network failures. Although the total cost is relatively low, the system's resilience and reliability are poor in the event of a system failure. Scheme 2 only considers power network failures. Although it has some countermeasures for power system failures, it ignores the impact of failures in other energy networks such as gas, heating, and cooling networks, resulting in higher costs for ship delays, electricity purchases, and total load shedding, leading to poor overall system performance and economy. Scheme 3 considers multi-energy network failures. Through dynamic fault constraint modeling and collaborative optimization of the power, gas, heating, and cooling multi-energy networks, it can effectively suppress the risk of cross-energy cascading failures and ensure port operation efficiency and energy supply and demand balance. Although its total cost is slightly higher than Scheme 1, it has significant advantages in system resilience, reliability, and economy, making it an effective solution for dealing with complex scenarios and achieving sustainable development of the port's integrated energy system.

[0080] Table 1 Comparison of System Costs ($) for Different Solutions

[0081] Option 1 90.00 7090.00 2578.80 2784.74 1543.06 / 14086.6 Option 2 240.00 4650.00 1978.45 4285.51 1256.00 9856.32 22266.28 Option 3 150.00 5763.00 2175.98 1972.96 1345.58 4978.45 16385.97

[0082] 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 method for enhancing the resilience of a port integrated energy system considering multi-energy network failures, characterized in that, The method includes the following steps: Step 1: Establish the objective function of the resilience enhancement model for the integrated energy system of the seaport; Step 2: Establish multi-energy network fault constraints for the resilience enhancement model of the integrated energy system of the seaport; Step 3: Establish combined cooling, heating and power (CCHP) constraints for the resilience enhancement model of the integrated energy system of the seaport; Step 4: Establish container transport constraints for the resilience enhancement model of the integrated energy system of the seaport; Step 5: Based on the objective function of Step 1 and the constraints of Steps 2-4, solve the resilience enhancement model of the integrated energy system of the seaport to obtain the resilience operation strategy of the integrated energy system of the seaport. In step (1), the objective function for establishing the resilience enhancement model of the integrated energy system of the seaport is expressed as follows: （1） In the formula, a represents a combined cooling, heating and power (CCHP) unit; t represents a time period; s represents a ship; e represents a grid node; n represents a new energy unit; d represents an electrolytic cell; b represents a boiler; p represents a heat pump; r represents an electric chiller; g represents a gas grid node; h represents a heating grid node; and c represents a cooling grid node. This represents the fixed cost of combined cooling, heating and power (CCHP) unit a; This represents the startup cost of combined cooling, heating and power (CCHP) unit a; This represents the shutdown cost of combined cooling, heating and power (CCHP) unit a; This indicates whether the combined cooling, heating and power (CCHP) unit a is fixed during time period t; if yes, set to 1, otherwise set to 0. This indicates whether the combined cooling, heating and power (CCHP) unit a is started during time period t; if yes, set to 1, otherwise set to 0. This indicates whether the combined cooling, heating and power (CCHP) unit a is off during time period t; set to 1 if yes, otherwise set to 0. This represents the power generation cost of combined cooling, heating and power (CCHP) unit a; This represents the active power output of the combined cooling, heating and power (CCHP) unit a during time period t. This represents the operation and maintenance cost of combined cooling, heating and power (CCHP) unit a; This indicates the penalty for late departure of vessel s; Indicates the actual departure time of vessel s; Indicates the expected departure time of vessel s; This represents the electricity price during time period t; This represents the active power injected by the upstream power grid into grid node e during time period t; This represents the operation and maintenance cost of the new energy unit n; This represents the active power output of the new energy unit n during time period t; This represents the operation and maintenance cost of electrolytic cell d; This represents the active power consumption of electrolytic cell d during time period t; This represents the operation and maintenance cost of boiler b; This represents the gas consumption of boiler b during time period t; This represents the operation and maintenance cost of heat pump p; This represents the active power consumption of the heat pump p during time period t; This represents the maintenance cost of the electric chiller r; This represents the active power consumption of the electric chiller r during time period t; Indicates the cost of electrical load loss; Indicates the cost of gas loss of load; Indicates the cost of thermal load loss; Indicates the cost of cooling load loss; This represents the load shedding rate of grid node e during time period t; This represents the total active load of grid node e during time period t; This represents the gas load loss rate of gas network node g during time period t; This represents the total gas load of gas network node g during time period t; This represents the heat load loss rate of heating network node h during time period t; This represents the total heat load of heat network node h during time period t; This represents the cooling load loss rate of cold network node c during time period t; This represents the total cooling load of cold network node c during time period t.

2. The method for enhancing the resilience of a port integrated energy system considering multi-energy network failures as described in claim 1, characterized in that, The specific process of step (2) is as follows: (201) Establish grid fault constraints: （2） （3） （4） In the formula, f represents a power grid node; o represents a berth; q represents a quay crane; y represents a yard crane; v represents an automated guided container truck; and ef represents the power grid branch connecting power grid nodes e and f. This represents the set of combined cooling, heating and power (CCHP) units connected to grid node e; This represents the set of new energy generating units connected to grid node e; This represents the set of electrolytic cells connected to grid node e; This represents the set of heat pumps connected to grid node e; This represents the set of electric chillers connected to the power grid node e; Represents the set of power grid branches with terminal node e; Represents the set of power grid branches with the first node being e; This represents the set of berths connected to the power grid node e; This represents the set of quay bridges connected to the power grid node e; This represents the set of field bridges connected to grid node e; This represents the set of autonomously navigated container trucks connected to grid node e; Let f represent the set of power grid branches with terminal node f; This represents the active power transmitted by the power grid branch ef during time period t; This represents the active power consumption of berth o during time period t; This represents the active power consumption of the quay crane q during time period t; This represents the active power consumption of the field bridge y during time period t; This represents the active power that the automated guided container truck v receives from the power grid during time period t; This represents the reactive power injected by the upstream power grid into grid node e during time period t; This indicates the reactive power output of the combined cooling, heating and power (CCHP) unit a during time period t; This represents the reactive power output of the new energy unit n during time period t; This represents the reactive power consumption of electrolytic cell d during time period t; This represents the reactive power consumption of the heat pump p during time period t; This represents the reactive power consumption of the electric chiller r during time period t; This represents the reactive power transmitted by the power grid branch ef during time period t; This represents the reactive power consumption of berth o during time period t; This represents the reactive power consumption of the quay crane q during time period t; This represents the reactive power consumption of the field bridge y during time period t; This represents the reactive power that the automated guided container truck v receives from the power grid during time period t; This represents the total reactive load of grid node e during time period t; It is a constant; This indicates whether the power grid branch ef is faulty during time period t; if yes, set to 0, otherwise set to 1. This represents the voltage amplitude at grid node e during time period t; This represents the voltage amplitude at grid node f during time period t; This represents the resistance of branch ef; Indicates the reactance of branch ef; Indicates the reference voltage; (202) Establish gas network fault constraints: （5） （6） （7） （8） In the formula, k represents a gas network node; l represents a gas compressor; gk represents the gas network pipeline connecting gas network nodes g and k; This represents the set of gas pipelines with the first node g. This represents the set of gas compressors connected to gas network node g. This represents the set of combined cooling, heating and power (CCHP) units connected to the gas network node g. This represents the set of electrolytic cells connected to the gas network node g; This represents the set of boilers connected to the gas network node g; This represents the hydrogen gas transported by the gas pipeline gk during time period t; This represents the compression coefficient of the gas compressor l; This represents the hydrogen output by gas compressor l during time period t; This represents the hydrogen consumption of combined cooling, heating and power unit a during time period t; This represents the conversion efficiency of electrolytic cell d during time period t; This indicates the gas transmitted by the gas pipeline gk during time period t; This represents the gas output of gas compressor l during time period t; This represents the amount of gas injected into gas network node g during time period t; This represents the gas consumption of combined cooling, heating and power unit a during time period t; Indicates the hydrogen blending ratio in the gas pipeline network; This indicates whether the gas pipeline gk is faulty during time period t; if yes, set to 0, otherwise set to 1. This indicates the maximum gas flow rate in the gas pipeline gk during time period t; This indicates the maximum hydrogen flow rate in the gas pipeline gk during time period t; (203) Establish heating network fault constraints: （9） （10） （11） （12） （13） In the formula, i represents a heat network node; hi represents the heat network pipe connecting heat network nodes h and i; This represents the set of combined cooling, heating and power (CCHP) units connected to the heating network node h; This represents the set of boilers connected to the heating network node h; This represents the set of heat pumps connected to the heat network node h. This represents the set of heating network pipes with terminal node h; This represents the set of heat network pipelines with the first node being h; This indicates the heat output of the combined cooling, heating and power (CCHP) unit a during time period t; This indicates the heat conversion efficiency of boiler b; This indicates the heat conversion efficiency of the heat pump p; This indicates the specific heat capacity of water; This represents the heat source mass flow rate of heat network node h during time period t; This represents the water supply temperature at node h in the heating network during time period t; This represents the return water temperature of heating network node h during time period t; This represents the heat load mass flow rate of heat network node h during time period t; This indicates the outlet temperature of the heating network pipe hi during time period t; This indicates whether the heating network pipeline hi is faulty during time period t; set to 0 if yes, otherwise set to 1. This represents the ambient temperature during time period t; This indicates the inlet temperature of the heating network pipe hi during time period t; This represents the heat transfer coefficient of the heating network pipe hi; Indicates the length of the heating network pipe hi; This represents the hot water mass flow rate of the heating network pipe hi during time period t; This represents the mixing temperature of heat network node h during time period t; (204) Establish cold network fault constraints: （14） （15） （16） （17） （18） In the formula, x represents a cold network node; cx represents the cold network pipe connecting cold network nodes c and x; This represents the set of combined cooling, heating, and power (CCHP) units connected to the cold network node c; This represents the set of electric chillers connected to node c of the cold network; Represents the set of cold network pipes with terminal node c; This represents the set of cold network pipes with the first end node being c; This represents the cooling output of the combined cooling, heating and power (CCHP) unit a during time period t. This indicates the refrigeration conversion efficiency of the electric chiller r; This represents the cold source mass flow rate of cold network node c during time period t; This represents the return water temperature of cold network node c during time period t; This represents the water supply temperature of cold network node c during time period t; This represents the mass flow rate of the cooling load at node c in the cold network during time period t. This indicates whether the cold network pipe CX is faulty during time period t; if yes, set to 0, otherwise set to 1. This represents the outlet temperature of the cold network pipe cx during time period t; This represents the inlet temperature of the cold network pipe cx during time period t; This represents the heat transfer coefficient of the cold network pipe cx; This indicates the length of the cold network pipe cx; This represents the mass flow rate of chilled water in the cold network pipe cx during time period t; This represents the mixing temperature of cold network node c during time period t.

3. The method for enhancing the resilience of a port integrated energy system considering multi-energy network failures according to claim 1, characterized in that, The specific process of step (3) is as follows: (301) Establishing the combined cooling, heating and power (CCHP) constraints for the resilience enhancement model of a seaport integrated energy system: （19） （20） （21） （22） （23） （24） In the formula, This indicates the power generation efficiency of combined cooling, heating and power unit a; This indicates the heat production efficiency of combined cooling, heating and power unit a; This indicates the maximum downward ramp rate of the combined cooling, heating and power (CCHP) unit a; This indicates the maximum upward ramp rate of the combined cooling, heating and power (CCHP) unit a; This indicates the maximum downhill ramp rate when the combined cooling, heating and power (CCHP) unit a is shut down. This indicates the maximum uphill ramp rate at startup for combined cooling, heating and power (CCHP) unit a; This represents the active power output of the combined cooling, heating and power (CCHP) unit a during the time period t-1; This indicates whether the combined cooling, heating and power (CCHP) unit a is fixed during the t-1 time period; if yes, set it to 1, otherwise set it to 0. This represents the minimum active power output of the combined cooling, heating and power (CCHP) unit a; This indicates the maximum active power output of the combined cooling, heating and power (CCHP) unit a.

4. The method for enhancing the resilience of a port integrated energy system considering multi-energy network failures according to claim 1, characterized in that, The specific process of step (4) is as follows: (401) Establishing container transport constraints for a resilience enhancement model of a seaport integrated energy system: （25） （26） In the formula: This represents the number of containers loaded and unloaded by quay crane q from ship s during time period t; This indicates whether the quay crane q serves vessel s during time period t; if yes, it is set to 1, otherwise it is set to 0. This indicates the number of containers loaded and unloaded from ship s by the yard crane y during time period t; This indicates whether the field bridge y serves vessel s during time period t; if yes, set to 1, otherwise set to 0. This represents the number of containers loaded and unloaded by the automated guided vehicle (AGV) from the ship (S) during time period t. This indicates whether the automated navigation container truck v serves the vessel s during time period t; set to 1 if yes, otherwise set to 0. This indicates the number of containers loaded on ship s; This represents the electrical energy storage of the automated guided container truck v during time period t; This represents the electrical energy storage of the automated guided container truck v during the time period t-1; This indicates the charging efficiency of the automated guided container truck v. This indicates the discharge efficiency of the automatically guided container truck v. This represents the active power consumption of the automated guided container truck v during time period t.

5. The method for enhancing the resilience of a port integrated energy system considering multi-energy network failures according to claim 1, characterized in that, In step (5), a resilience enhancement model for the integrated energy system of a port with multi-energy network failure is written in GAMS, and the resilience operation strategy of the integrated energy system of the port is obtained by solving the model.

Images