An inner-outer cooperative port comprehensive energy system low-carbon economic optimization scheduling method

Abstract

The application discloses a kind of inside and outside synergies port comprehensive energy system low-carbon economic optimization scheduling method, comprising the following steps: S1, obtain the output prediction data of photovoltaic, wind power and cold and heat load, time-of-use electricity price, the capacity size of various equipment and output constraint;S2, based on carbon flow characteristics build green certificate-carbon trading mechanism model and comprehensive demand response model;S3, with the minimum operation cost after inside and outside synergies as objective function, port objective function and port area objective function are constructed;S4, port objective function solves the coordinated optimization problem outside port, and the interactive power obtained by solving is passed to basic energy supply unit, and port area objective function solves internal optimization problem of port area.The present application fully considers the carbon trading mechanism of port area and the energy interaction between multiple port areas, which can effectively improve the economic benefit and low-carbon characteristics of port comprehensive energy system.

Description

TECHNICAL FIELD

[0001] The present application relates to the field of integrated energy system optimal operation, in particular to a port integrated energy system low-carbon economic optimal scheduling method with internal and external coordination. BACKGROUND

[0002] With the increasingly serious problem of carbon emissions in the shipping industry, as the hub of shipping, it is urgent to realize the green and low-carbon construction of the port. And the coordinated development between various energy systems is the current mainstream trend, so the port integrated energy system which can realize the coupling and interconnection of multiple energies has gradually attracted the attention of the energy industry.

[0003] In the existing research on the optimization of port integrated energy systems, most of them adopt the mode of separate planning, separate design and independent operation, without considering the source-load economic regulation and mutual aid between various types of energy, the problems such as different load types and power characteristics, large configuration redundancy of supply-side energy supply, strong uncertainty of new energy output, etc. have not been effectively solved. At the same time, carbon emission constraints further reduce the effective way of mutual aid between various types of energy. In order to realize the low-carbon economic operation of the port integrated energy system, on the one hand, the interconnection and coupling complementation of various energies are needed; on the other hand, it is also necessary to adopt carbon trading and demand response mechanism to improve energy use and improve the flexibility of the system. Therefore, it is of great research significance to consider the economy and environmental protection of the system in the internal and external coordinated optimal scheduling of IPES. SUMMARY

[0004] The purpose of the present application is to provide a port integrated energy system low-carbon economic optimal scheduling method with internal and external coordination, which can effectively reduce the operation cost of the system and improve the low-carbon characteristics by fully considering the green-carbon trading mechanism of the port area and the energy interaction between multiple port areas.

[0005] Technical scheme: A port integrated energy system low-carbon economic optimal scheduling method with internal and external coordination, wherein the port integrated energy system realizes transaction with the upper energy market through the port power distribution network and the natural gas network, the port is composed of multiple port areas, each port area is a basic energy supply unit, and the port areas interact with each other through the tie line to exchange electric energy and thermal energy; the equipment in the port area mainly includes photovoltaic, wind turbine, gas boiler, electric refrigerator, combined heat and power and electric storage and heat storage device; comprising the following steps:

[0006] S1, obtaining the output prediction data of photovoltaic, wind power and combined heat and power load, time-of-use electricity price, capacity size and output constraint of various devices;

[0007] S2, building a green certificate-carbon trading mechanism model and a comprehensive demand response model based on carbon flow characteristics;

[0008] S3, taking the minimum operation cost after internal and external coordination as the objective function, a port objective function and a port area objective function are constructed;

[0009] S4, the port objective function solves the coordination optimization problem outside the port and passes the obtained interaction power to the basic energy supply unit, and the port area objective function solves the internal optimization problem.

[0010] Further, in step S2, based on the carbon flow characteristics of the port comprehensive energy system, a green certificate-carbon trading mechanism is built, and the carbon trading interval and green certificate benefit are calculated;

[0011] The carbon flow model of the input side of the port comprehensive energy system is:

[0012] E e,t = P e,t ρ e,t

[0013] E g,t = P g,t ρ g,t

[0014] In the formula, E e,t , E g,t are carbon emissions brought by electricity and gas purchase at time t; P e,t , P g,t are power supply of the power grid and the gas grid at time t; ρ e,t , ρ g,t are carbon flow densities corresponding to the power supply of the power grid and the gas grid.

[0015] According to the carbon flow conservation of device input and output, the carbon flow model of the energy hub is:

[0016]

[0017] In the formula, are electric power, heat power and cold power generated by the CCHP at time t; are carbon flow densities of the electric power, heat power and cold power generated by the CCHP; are heat power and corresponding carbon flow density generated by the GB at time t; are electric power and cold power generated by the EC at time t; are carbon flow densities of the electric power and cold power generated by the EC;

[0018] The carbon emissions of the output side of the port comprehensive energy system are:

[0019]

[0020] In the formula, E Le,t , E Lh,t, E Lc,t , E PV,t , E WT,t , E se,t , E sh,t , E PV,t , E WT,t , E se,t , E sh,t , E , E

[0021]

[0022] , E , E , E , E , E , E

[0023] The green certificate trading mechanism of the port comprehensive energy system is:

[0024]

[0025] , E , E GCT , E , E , E , E

[0026] The port comprehensive demand response model is:

[0027]

[0028] , E e,dr,t , E , E , E

[0029] The port electric energy Pe,port,t The demand is:

[0030]

[0031] In the formula, is the port power demand without considering comprehensive demand response, in units of MW;

[0032] The controllable load constraint is:

[0033]

[0034] In the formula, is the proportion of the port power load that can be reduced in the total power load; are the proportions of the power load transferred from the t period to other periods and from other periods to the t period, respectively, in the total power load; are whether to perform power load transfer, 0 indicating no transfer and 1 indicating transfer.

[0035] Further, in step S4, the objective function established with the minimum total daily operation cost is a cost function composed of energy interaction cost energy storage cost , then the expression of the port objective function is:

[0036]

[0037] In the formula, N is the total number of port areas; C1 is the total daily operation cost, are the energy interaction cost and energy storage cost between the nth port area, respectively; g e , g h is the electricity and heat transaction cost; is the transmission electric and heat power of the ith EH device in the port area, respectively represent the electric power delivered to the EH i from the EH j , the EH j to the EH i at t time, i and j represent the number of EH devices; respectively represent whether the EH i transmits power to the EH j , the EH j to the EH i at t time, 0 indicates no power transmission and 1 indicates power transmission; d i,j is the distance between the EH i and the EH j ; σ is the tie line loss rate, and energy cannot be bidirectionally transmitted in the same line at the same time; R ε is the operation and maintenance coefficient of the ε device, ε ∈ {ES, HS}. The charging and discharging power corresponding to the ES and HS devices;

[0038] The energy storage and device constraints are as follows:

[0039]

[0040] The tie-line interaction power constraint:

[0041]

[0042] In the formula, is the upper and lower limits of the energy storage of the ε device; is the upper and lower limits of the energy discharge of the ε device; P i,j,max is the maximum transmissible capacity, in kW;

[0043] The port area establishes an objective function with the minimum all-day operation cost, and the cost function is composed of the energy purchase cost C BUY,t , the device operation cost C OP,t , the demand response cost C n,dr,t , the carbon trading cost and the green certificate trading cost C GCT,t , so the expression of the port area objective function is as follows:

[0044]

[0045] C BUY,t = P e,t g e,t + P g,t g g,t

[0046]

[0047] In the formula, C BUY,t , C OP,t , C dr,t , C GCT,t , are the energy purchase cost, the device operation cost, the demand response cost, the carbon trading cost and the green certificate trading cost of the port area respectively; g e,t , g g,t are the time-of-use electricity price and the gas price; ξ ∈ {CCHP, PV, WT, GB, ES, HS}; P ξ t is the operation capacity power of each unit device, in MW; R ξ is the operation cost coefficient of each device; are the load reduction and transfer costs, yuan / MWh respectively; T is the length of a dispatching period, and Δt is the length of the considered time period;

[0048] The energy hub coupling constraints are as follows:

[0049]

[0050] L e,t , L ship,t , L h,t , L c,t respectively represent total demand of electric load, ship shore power load, heat load and cold load; P PV,t , P WT,t respectively represent photovoltaic and wind power input; H g represents natural gas heat value; and respectively represent energy conversion efficiency of CCHP electricity, heat and cold; alpha, beta and gamma represent distribution coefficients corresponding to CCHP electricity, heat and cold; P ship,t represents total ship shore power load at t moment; m represents number of ships at t moment; represents ship shore power load of the i-th ship at t moment; represents state of the i-th ship connecting to port shore power at t moment, when the ship connects to port shore power when the ship does not connect to port shore power

[0051] The electricity and gas purchase and equipment power constraints are as follows:

[0052]

[0053] P e,max , P e,min respectively represent upper and lower limits of electricity purchase; P g,max , P g,min respectively represent upper and lower limits of gas purchase; P i represents power of each equipment; P i,min represents lower limit of power of each equipment; P i,max represents upper limit of power of each equipment.

[0054] Further, the obtained interactive power is the interactive power of electric energy and heat energy between ports.

[0055] Compared with the prior art, the present application has the following remarkable effects:

[0056] 1. The present application adopts a hierarchical optimization method based on target cascade analysis for solving, the upper target function solves the coordination optimization problem outside the port, and the obtained interactive power is transmitted to the lower layer, the lower target function solves the internal optimization problem of the port area, fully considers the carbon trading mechanism of the port area and the energy interaction between multiple port areas, and can effectively improve the economic benefit of the port comprehensive energy system.

[0057] 2、The green-carbon trading mechanism and the incentive comprehensive demand response are considered, so that the optimal scheduling of the port integrated energy system can be realized through the energy interaction between the ports on the basis of the coordinated operation of a single port, thereby effectively reducing the operation cost of the system and improving the low-carbon characteristics. BRIEF DESCRIPTION OF DRAWINGS

[0058] Figure 1 The port integrated energy system operation model in the embodiment of the present application is shown in the figure.

[0059] Figure 2 The flow chart of the low-carbon economic optimization method of the port integrated energy system in the embodiment of the present application is shown in the figure.

[0060] Figure 3 (a) is the carbon emission interval of the bulk cargo ship port, (b) is the carbon emission interval of the cruise ship port, and (c) is the carbon emission interval of the container ship port.

[0061] Figure 4 (a) is the electric energy scheduling strategy of the bulk cargo port, (b) is the electric energy scheduling strategy of the cruise ship port, and (c) is the electric energy scheduling strategy of the container ship port.

[0062] Figure 5 (a) is the thermal energy scheduling strategy of the bulk cargo port, (b) is the thermal energy scheduling strategy of the cruise ship port, and (c) is the thermal energy scheduling strategy of the container ship port. DETAILED DESCRIPTION

[0063] The present application will be further described in detail in combination with the drawings and specific embodiments of the present application.

[0064] In the embodiment, first, the framework of the port integrated energy system is built, as shown in the figure. Figure 1 The port integrated energy system realizes the transaction with the superior energy market through the port power distribution network and the natural gas network, the port is composed of multiple port districts, each port district is a basic energy supply unit, and the port districts interact with each other through the tie lines to exchange electric energy and thermal energy, thereby meeting the load demand of the user side. The equipment in the energy hub mainly includes photovoltaic, fan, gas boiler (GB), electric chiller (EC), combined cooling, heating and power (CCHP), and electric storage and heat storage devices. Based on the carbon energy flow characteristics of the port integrated energy system (PIES), the green-carbon trading mechanism is introduced to calculate the carbon trading interval and the green certificate benefit; and according to the multi-energy characteristics of the port integrated energy system, an incentive comprehensive demand response model is established to guide the change of the load energy consumption mode.

[0065] Figure 2 This is a flowchart of a low-carbon economic optimization method for a port integrated energy system in an embodiment of the present invention.

[0066] according to Figure 2 The hierarchical optimization method based on target cascade analysis shown for system optimization scheduling includes the following steps:

[0067] Step 1: Obtain output forecast data for photovoltaic, wind power, and cooling, heating and power loads, time-of-use electricity prices, capacity size of various equipment, and output constraints.

[0068] Step 2: Build a green certificate-carbon trading mechanism model and a comprehensive demand response model;

[0069] (21) PIES Green-Carbon Trading Model

[0070] In this embodiment, based on the carbon flow characteristics of the Port Integrated Energy System (PISE), a green certificate-carbon trading mechanism is introduced to calculate the carbon trading range and the benefits of green certificates.

[0071] The carbon flow model on the PIES input side is:

[0072] E e,t =P e,t ρ e,t (1)

[0073] E g,t =P g,t ρ g,t (2)

[0074] In the formula, E e,t E g,t These represent the carbon emissions from purchasing electricity and gas at time t, respectively; P e,t P g,t The power supplied by the power grid and gas grid at time t are respectively; ρ e,t ρ g,t These are the carbon flux densities corresponding to the power supplied by the power grid and gas grid, respectively.

[0075] Based on the conservation of carbon flow at the equipment input and output, the carbon flow model for the energy hub is derived as follows:

[0076]

[0077] In the formula, These represent the electrical, thermal, and cooling power produced by CCHP at time t, respectively. These are the carbon flux densities for the electrical, thermal, and cold power output of CCHP, respectively. These represent the power of the thermal energy produced by GB at time t and the corresponding carbon flux density, respectively. These represent the electrical and cooling power produced by EC at time t, respectively. Carbon flow density of EC output electricity and cold power, respectively.

[0078] Carbon emission of PIES output side, respectively:

[0079]

[0080] In the formula, E Le,t , E Lh,t , E Lc,t Carbon emission of electricity, heat, and cold load at t moment, respectively; e PV,t , e WT,t Carbon potential of photovoltaic and wind turbine generator, respectively; e se,t , e sh,t Carbon potential of electricity and heat energy storage device when discharging; P PV,t , P WT,t Electricity power of photovoltaic and wind power at t moment, respectively; P se,t , P sh,t Discharging power of electricity and heat energy storage device at t moment, respectively; Thermal power of GB at t moment.

[0081]

[0082] In the formula: Carbon trading cost of load i at t moment; Carbon trading interval range of marginal action obtained by shapley value method, respectively; Carbon emission of load i at t moment; Carbon flow density of load i at t moment; Load amount of load i at t moment; λ1, λ2 and λ3 are three-grade carbon trading unit prices, unit: yuan / t.

[0083] Green certificate trading mechanism of PIES is:

[0084]

[0085] In the formula, Green certificate trading cost at t moment; α GCT Green certificate trading price; Green certificate amount purchased / sold at t moment; Total electricity consumption at t moment; Renewable energy consumption amount at t moment; k is a quantitative coefficient, 1 green certificate corresponds to 1 MW·h renewable energy electricity; ω is the weight of power consumption responsibility.

[0086] (22) Port comprehensive demand response model is:

[0087]

[0088] In the formula, P e,dr,t The power load participating in integrated demand response (IDR) at time t, in MW; The load reduction at time t is expressed in MW. These represent the load transferred from time t to other time periods and the load transferred from other time periods to time t, respectively.

[0089] Port power P after integrated demand response e,port,t The requirement is:

[0090]

[0091] In the formula: The port's power demand is measured in MW, without considering the overall demand response.

[0092] The controllable load constraint is:

[0093]

[0094] In the formula: This can reduce the proportion of the port's electricity load in the total electricity load; These represent the proportions of electricity load that can be transferred from time period t to other time periods and the proportions of electricity load that can be transferred from other time periods to time period t, respectively, to the total electricity load. These indicate whether to perform power load transfer, with 0 indicating no transfer and 1 indicating transfer.

[0095] Step 3: Using the minimum operating cost after internal and external coordination as the objective function, construct a two-layer optimization scheduling model for the port integrated energy system that takes into account internal and external coordination;

[0096] In the two-layer optimization scheduling model of the port's integrated energy system, which takes into account both internal and external coordination, the objective function established by the port is to minimize the daily operating cost.

[0097] The port optimization objective is to minimize the total daily operating cost, and the cost function is determined by the energy interaction cost. Energy storage costs The objective function of the port is composed of the following:

[0098]

[0099] In the formula: N is the total number of port areas, n represents the nth port area; C1 is the total daily operating cost. For the energy exchange cost and energy storage cost between the nth port area; g e g h Costs associated with electricity and heat energy transactions; The transmission electric and thermal power of the i-th EH device in the port area, EH i to the EH j , EH j to the EH i The electric power delivered to the EH device, i and j are the numbers of the EH devices respectively; EH i to the EH j , EH j to the EH i Whether to transmit power, 0 represents no power transmission, and 1 represents power transmission; d i,j The distance between the EH i and the EH j ; σ is the tie line loss rate; In the same line at the same time, energy cannot be transmitted bidirectionally. ε ε is the operation and maintenance coefficient of the ε device, and ε ∈ {ES, HS}; ES (electricity storage) and HS (heat storage) devices correspond to the charging and discharging power.

[0100] The constraints of the port objective function are: energy storage and device constraints and tie line interaction power constraints;

[0101] The energy storage and device constraints are as follows:

[0102]

[0103] In the formula: are the upper and lower limits of the energy storage of the ε device respectively; are the upper and lower limits of the energy storage of the ε device respectively.

[0104] The tie line interaction power constraints are as follows:

[0105]

[0106] In the formula: P i,j,max is the maximum transmissible capacity, with the unit of kW.

[0107] In the port comprehensive energy system double-layer optimization scheduling model considering internal and external collaboration, the objective function established by the port area to minimize the whole day operation cost is:

[0108] The optimization objective of the port area is to minimize the whole day operation cost, and the cost function is composed of the energy purchase cost C BUY,t , the device operation cost C OP,t , the demand response cost C n,dr,t , the carbon trading cost and the green certificate trading cost C GCT,tThe expression of the port objective function is as follows:

[0109]

[0110] C BUY,t = P e,t g e,t + P g,t g g,t (22)

[0111]

[0112] In the formula, C BUY,t , C OP,t , C dr,t , C GCT,t are the energy purchasing cost, equipment operation cost, demand response cost, carbon trading cost and green certificate trading cost of the port respectively; g e,t , g g,t are the time-of-use electricity price and gas price; ξ ∈ {CCHP, PV, WT, GB, ES, HS}; P ξ t is the operation power of each unit equipment, in MW; R ξ is the operation cost coefficient of each equipment; are the load reduction and transfer costs respectively, in yuan / MWh; T is the length of a dispatching period, and Δt is the length of the considered time period.

[0113] The constraint conditions of the port objective function are: energy hub coupling constraint and electricity and gas purchasing and equipment power constraint;

[0114] The energy hub coupling constraint is as follows:

[0115]

[0116] In the formula, L e,t , L ship,t , L h,t , L c,t are the total demand of the electric load, ship shore power load, heat load and cold load respectively; P PV,t , P WT,t are the photovoltaic and wind power inputs respectively; H g is the natural gas heat value; and are the energy conversion efficiencies of CCHP electricity, heat and cold respectively; α, β, γ are the distribution coefficients of CCHP electricity, heat and cold.

[0117]

[0118] In the formula, P ship,tis the total ship shore power at time t; m is the number of ships in port at time t; is the shore power of the i-th ship at time t; is the shore power of the i-th ship at time t; is the shore power of the i-th ship at time t;

[0119] The electricity and gas purchase constraints are as follows:

[0120]

[0121] wherein: P e,max , P e,min are the upper and lower limits of electricity purchase; P g,max , P g,min are the upper and lower limits of gas purchase; P i is the power of each device; P i,min is the lower limit of the power of each device; P i,max is the upper limit of the power of each device.

[0122] Step 4: solving by a hierarchical optimization method based on target cascade analysis.

[0123] Equation (15) is used as the upper layer (i.e. port) target function to solve the coordination optimization problem outside the port, and the obtained interactive power (i.e. the interactive power of electric energy and thermal energy between the port and the port area) is passed to the lower layer (i.e. the port area), and equation (21) is used as the lower layer target function to solve the optimization problem inside the port area.

[0124] In summary, the port integrated energy system is optimized by the present application, fully considering the carbon trading mechanism of the port area and the energy interaction between multiple port areas, which can effectively improve the economic benefit and low-carbon characteristics of the port integrated energy system.

[0125] The effectiveness of the proposed method is verified by simulation analysis of the inland river port in Suzhou. The IPES composed of bulk cargo port area, cruise ship port area and container port area is optimized and dispatched, and six scenarios are set as shown in Table 1.

[0126] Table 1 Scenario setting

[0127] Scenario Carbon trading Clean Development Mechanism IDR EHs in synergy Ⅰ Ⅱ √ Ⅲ √ √ Ⅳ √ Ⅴ √ √ √ Ⅵ √ √ √ √

[0128] (1) Carbon trading result analysis

[0129] The step carbon trading results of each port area are analyzed. The carbon emission quota of the bulk cargo ship port area at t=1h is taken as an example to set the interval, and the carbon trading unit price of three grades is taken as 5, 10 and 20 yuan / t respectively, and the carbon quota of 24h is obtained EH(s=1, 2, 3; respectively denoted as EH1, EH2 and EH3) as shown in Table 2.Figure 3 As shown in (a), (b), and (c) in the figure.

[0130] It can be seen that the carbon emission allowance for EH3 is larger than that for EH1 and EH2, because the container port area has the highest load and the largest carbon emissions. Since the port's demand for electricity and gas loads is high during the day, and electricity and gas loads account for a large proportion, the carbon emission allowance ranges for each EH area from 9:00 to 20:00 are all high, belonging to the high-carbon period, while those from 0:00 to 8:00 are low, belonging to the low-carbon period. Compared to using a single range to calculate the carbon emission liability cost for each time period, the IPES carbon trading model based on the Shapley value method is more reasonable.

[0131] (2) Scheduling Result Analysis

[0132] The optimized scheduling results for different scenarios are shown in Table 2.

[0133] Table 2 Optimization scheduling results under different scenarios

[0134]

[0135] Comparing scenarios I, IV, II, and V, when considering the hierarchical coordination and optimization of IPES, the operating costs are all reduced compared to the independent operation of each port area. Therefore, the internal and external coordination and optimization of IPES, by realizing energy sharing among the port areas, makes energy use more rational. However, without considering low-carbon economic optimization, compared to scenarios I and IV, the port's electricity purchases increase, leading to a corresponding increase in carbon emissions. Therefore, although the internal and external coordination strategy of IPES has certain economic benefits, it will have a certain negative impact on carbon emissions.

[0136] When only considering the carbon trading mechanism and not considering the green certificate mechanism, compared with scenarios I and II, the operation cost, carbon emission and electricity purchase amount are all reduced, and the gas purchase amount is increased by 10.9 MW·h. This is because under the carbon trading mechanism, the port will be more inclined to use cleaner CCHP to supply power. When considering both the carbon trading mechanism and the green certificate mechanism, compared with scenarios IV and V, the operation cost is reduced by 2086 yuan, the carbon emission is reduced by 8.2 tons, the electricity purchase amount is reduced by 8.3 MW·h, and the gas purchase amount is increased by 16.3 MW·h. Therefore, when considering carbon trading and green certificate trading, although the carbon income is reduced, the system can reduce the carbon emission and bring green certificate income. Compared with scenarios II and III, after introducing the green certificate mechanism, the carbon emission is reduced by 1.1 tons, and the operation cost is reduced by 979 yuan. It shows that the green certificate mechanism further promotes the reduction of carbon emission and helps to balance the economic benefit of the system. Compared with scenarios V and VI, under the condition of considering carbon trading, after considering demand response in scenario VI, the total cost of the port energy system is reduced by 3255 yuan, the electricity purchase amount is reduced, the natural gas purchase cost is increased, and the carbon emission is reduced. It shows that the introduction of IDR makes IPES increase the utilization of gas network and reduce the power supply of power grid, thereby reducing the carbon emission. Taking scenario VI as an example, the dispatching strategies of the three port areas of IPES are shown in (a), (b) and (c) of Figure 4 , (a), (b), (c) of Figure 5 , (a), (b), (c) of

[0137] The application provides a port comprehensive energy system low-carbon economic optimization dispatching method based on internal-external coordination. By considering the carbon trading mechanism of green power certificates and incentive comprehensive demand response, a double-layer optimization dispatching model of the port comprehensive energy system considering internal-external coordination is established; a hierarchical optimization method based on target cascade analysis is used for solving, so as to realize the low-carbon economic dispatching of the port. The optimization model fully considers the carbon trading mechanism of the port area and the energy interaction between multiple port areas, and effectively improves the economic benefit and low-carbon characteristics of the port comprehensive energy system.

[0138] The above only describes the preferred embodiments of the application, and it should be noted that for ordinary skilled persons in the art, without departing from the principles of the application, some improvements and refinements can be made, and these improvements and refinements should also be considered as the protection scope of the application.

Claims

1. A method for low-carbon economic optimal dispatch of an internal-external coordinated port integrated energy system, wherein, The port comprehensive energy system realizes transaction with the superior energy market through a port power distribution network and a natural gas network, the port is composed of multiple port districts, each port district is a basic energy supply unit, and the port districts interact with each other through tie lines; the equipment in the port district mainly includes photovoltaic, wind turbine, gas boiler, electric refrigerator, combined heat and power and storage and heat storage device; characterized in that, comprising the following steps: S1, obtain the output prediction data of photovoltaic, wind power and cold heat and power load, time-of-use electricity price, capacity size and output constraint of various equipment; S2, build a green certificate-carbon trading mechanism model and a comprehensive demand response model based on carbon flow characteristics; S3, taking the minimum operation cost after internal and external cooperation as the objective function, a port objective function and a port district objective function are constructed; S4, the port objective function solves the coordination and optimization problem of the outside of the port, and the obtained interactive power is transmitted to the basic energy supply unit, and the port district objective function solves the internal optimization problem of the port district; In step S4, the target function is established with the minimum total daily operation cost, and the cost function is composed of the energy interaction cost and the energy storage cost , and the expression of the port target function is as follows: In the formula, N is the total number of ports; is the total operating cost of the day, , is the energy interaction cost and energy storage cost between the nth port, respectively; g e , g h is the electricity and heat transaction cost; , is the transmission electric and heat power of the ith EH device in the port, , respectively represent the electric power delivered by EH i to EH j , EH j to EH i at time t , i, j represent the number of EH devices; , ∈{0,1}, respectively represent whether EH i transmits power to EH j , EH j transmits power to EH i at time t , 0 represents no power transmission, and 1 represents power transmission; d i,j is the distance between EH i and EH j ; is the tie line loss rate, and energy cannot be transmitted bidirectionally in the same line at the same time; is the operating and maintenance coefficient of the device, ; , is the charging and discharging power corresponding to the ES and HS devices; The energy storage and equipment constraints are as follows: The tie line interactive power constraint is as follows: In the formula, , is Upper and lower limits of the energy storage of the device; , is Upper and lower limits of the energy release of the device; P i,j,max is the maximum transmissible capacity, in kW; The port area establishes a target function with the minimum total operation cost in a day, and the cost function is composed of the cost of purchasing energy , the operation cost of equipment , the cost of demand response , the cost of carbon trading and the cost of green certificate trading , and the expression of the target function of the port area is as follows: In the formula, C BUY,t , C OP,t , C dr,t , C co2,t , C GCT,t , are the port area's purchase cost, equipment operation cost, demand response cost, carbon trading cost, and green certificate trading cost, respectively; g e,t , g g,t is the time-of-use electricity price and gas price; ; is the operating capacity power of each unit equipment, with the unit being MW; is the operating cost coefficient of each equipment; , are the load reduction and transfer costs, respectively, in yuan / MWh; T is the length of a dispatching period, is the length of the considered time period; The energy hub coupling constraint is as follows: In the formula: L e,t , L ship,t , L h,t , L c,t These are the total demand for electrical load, ship-to-shore power load, heat load, and cooling load, respectively. P PV,t , P WT,t The inputs are solar power and wind power, respectively. H g This refers to the calorific value of natural gas. , and These are the energy conversion efficiencies of CCHP for electricity, heat, and cold, respectively. , , The distribution coefficients for electrical, thermal, and cold energy in CCHP; P ship,t for t Total ship shore power load at any given time; m for t Number of ships berthing at port at any given time; for t Time of the first i Only the shore power load capacity of ships; for t Time of the first i Only when the ship is connected to the port shore power status, when the ship is connected to the port shore power status. When the ship is not connected to the port shore power The power and gas purchase and equipment power constraint is as follows: In the formula, P e,max , P e,min are upper and lower limits of electricity purchase amount, respectively; P g,max , P g,min are upper and lower limits of gas purchase amount, respectively; P i is power of each device; P i,min is lower limit of power of each device; P i,max is upper limit of power of each device; The obtained interactive power is the interactive power of the port district.

2. The method of claim 1, wherein, In step S2, based on the carbon flow characteristics of the port comprehensive energy system, a green certificate-carbon trading mechanism is built, and the carbon trading interval and the green certificate benefit are calculated; The carbon flow model of the input side of the port comprehensive energy system is: In the formula, E e,t , E g,t respectively are t carbon emissions brought by electricity and gas purchase at the moment; P e,t , P g,t respectively are t power supply of the power grid and the gas grid at the moment; , respectively are carbon flow density corresponding to power supply of the power grid and the gas grid. According to the carbon flow conservation of equipment input and output, the carbon flow model of the energy hub is obtained as follows: In the formula, , , are respectively the electric, heat and cold power output by the CCHP at the moment t; t , , , are respectively the carbon flow density of the electric, heat and cold power output by the CCHP; , are respectively the power of the heat energy output by the GB at the moment t and the corresponding carbon flow density; t , , are respectively the electric and cold power output by the EC at the moment t; t , , are respectively the carbon flow density of the electric and cold power output by the EC; The carbon emissions of the output side of the port comprehensive energy system are as follows: In the formula, , , They are respectively t Carbon emissions from electricity, heat, and cooling loads at all times; , These are the carbon potentials of photovoltaic and wind turbine units, respectively. , These are the carbon potentials of the electric and thermal energy storage devices when they release energy, respectively. , They are respectively t The electrical power output of photovoltaic and wind power at any given time; , They are respectively t The energy output of the electric and thermal energy storage devices at all times; for t The thermal power of GB at any given time; In the formula, Load i At t Carbon trading cost at the moment; , , The marginal effect of carbon trading interval range obtained by shapley value method respectively; Load i At t Carbon emissions at the moment; Load i At t Carbon flow density at the moment; Load i At t Load at the moment; , And Three carbon trading unit prices, unit: yuan / t; The green certificate trading mechanism of the port comprehensive energy system is as follows: In the formula, is t the green certificate transaction cost at the moment; is the green certificate transaction price; is t the green certificate quantity purchased / sold at the moment; is t the total power consumption at the moment; is t the renewable energy consumption at the moment; k is a quantitative coefficient, 1 green certificate corresponds to 1 MW·h of renewable energy power; is the power consumption responsibility weight; The port comprehensive demand response model is as follows: In the formula, P e,dr,t for t The power load that participates in integrated demand response at all times, measured in MW; for t Load reduction amount at any time; , They are respectively from t Time shift to other time slots and from other time slots to t The workload of the moment; Port power after integrated demand response Demand is: In the formula, Pport is the port electric energy demand without considering integrated demand response, in MW. The controllable load constraint is as follows: In the formula, The proportion of the power load that can be reduced by the port in the total power load; , The power load transferred from the other time period to the t The proportion of the power load transferred from the other time period to the t The proportion of the power load transferred from the other time period to the , ∈{0, 1}, respectively, whether to transfer the power load, 0 means not to transfer, 1 means to transfer.

Images