A cross-border integrated energy system collaborative operation optimization method considering carbon-green certificate combined transaction

Abstract

The present application relates to a kind of cross-border integrated energy system collaborative operation optimization method considering carbon-green certificate combined transaction, first, build cross-border integrated energy system CBIES, cross-border integrated energy system CBIES including gas turbine, gas boiler, renewable energy unit, energy storage device, electric boiler, wherein renewable energy unit includes wind power generation and photovoltaic power generation, energy storage device includes battery and heat storage tank;Second, propose the cross-border integrated energy system model considering heat network, power grid loss, on the basis of fully considering cross-border energy flow and conversion, introduce cross-border green certificate, carbon transaction, establish the cross-border integrated energy system collaborative optimization model considering carbon-green certificate combined transaction;Finally, the collaborative operation optimization result is obtained by the model simulation established.This application provides optimization method, which can reduce the operation cost of integrated energy system in different countries, and also improve energy optimization configuration and renewable energy consumption amount.

Description

TECHNICAL FIELD

[0001] The application discloses a cross-border comprehensive energy system collaborative operation optimization method considering carbon-green certificate joint transaction, belongs to the field of electric power engineering research, and particularly relates to the cross-border comprehensive energy system collaborative operation optimization method based on consideration of carbon-green certificate joint transaction. BACKGROUND

[0002] The border areas of China's southwest region and Southeast Asian countries have high population density, rapid growth of trade and energy demand, rapid development of cooperation in cross-border trade, tourism, processing industry and agriculture, and great future development potential. Energy interconnection is an important part of China's foreign cooperation, and cross-border energy cooperation has become an important link to promote rapid economic development, plays a crucial role in realizing energy cooperation benefits and destiny, and ensuring national energy security. However, as an important part of China's foreign cooperation, energy interconnection faces problems such as cross-border market transaction risks in the aspects of laws, policies and finance, great differences in bilateral dispatching mechanisms, and limited information exchange. Incomplete system construction and technical means cannot guarantee the safe and stable operation of the cross-border comprehensive energy system, and in addition, may reduce the enthusiasm of cross-border parties, thereby affecting the optimal allocation of cross-border resources and the improvement of overall welfare.

[0003] The existing cross-border energy cooperation mainly takes the form of single energy, and is mainly point-to-point transmission of electricity and natural gas pipelines, which has a long construction period, low channel utilization rate, and lacks bidirectional supply capacity, local new energy consumption capacity and resource complementary utilization capacity. SUMMARY

[0004] The application proposes a cross-border comprehensive energy system collaborative operation optimization method considering carbon-green certificate joint transaction to solve the problems of different energy market mechanisms and great difficulty in multilateral collaborative dispatching between China and neighboring countries, and is used for determining the optimal scale of existing cross-border collaborative operation and large-scale utilization of new energy.

[0005] A cross-border comprehensive energy system collaborative operation optimization method considering carbon-green certificate joint transaction, characterized by comprising the following steps:

[0006] S1, a cross-border comprehensive energy system CBIES model is constructed, the cross-border comprehensive energy system CBIES comprises a gas turbine, a gas boiler, a renewable energy unit, an energy storage device and an electric boiler, wherein the renewable energy unit comprises wind power generation and photovoltaic power generation, and the energy storage device comprises a battery and a heat storage tank;

[0007] S2, a cross-border green certificate transaction model and a carbon transaction model are introduced on the basis of the cross-border comprehensive energy system CBIES model in step S1, and a cross-border comprehensive energy system collaborative operation optimization model considering carbon transaction and cross-border green certificate transaction is established;

[0008] S3, solving the cross-border integrated energy system collaborative operation optimization model considering carbon trading and cross-border green certificate trading established in step S2 using a solver to obtain an optimization result.

[0009] S11, the CBIES model constructed in step S1 specifically includes a gas turbine model, a gas boiler model, a renewable energy unit model, an energy storage device model, and an electric boiler model, and is specifically as follows:

[0010] The gas turbine model is specifically as follows:

[0011]

[0012]

[0013] wherein represents the output electric power of the gas turbine of the i-th country integrated energy system at time t, represents the output thermal power of the gas turbine of the i-th country integrated energy system at time t, L NG represents the low heat value of natural gas, represents the natural gas consumption of the gas turbine at time t, η GT represents the power generation efficiency of the gas turbine, N GT represents the waste heat recovery coefficient;

[0014] The gas boiler model is specifically as follows:

[0015]

[0016] wherein represents the output thermal power of the gas boiler of the i-th country integrated energy system at time t, L NG represents the low heat value of natural gas, represents the natural gas consumption of the gas boiler at time t, η GB represents the heat generation efficiency of the gas boiler;

[0017] The renewable energy unit model includes a wind power generation model and a photovoltaic power generation model, and is specifically as follows:

[0018] The wind power generation model is:

[0019]

[0020] wherein is the wind power of the i-th country integrated energy system at time t, ρ is the air density, A is the swept area of the wind turbine blade, v is the wind speed, c wt is the wind energy utilization coefficient, which is the ratio of the wind energy absorbed by the wind turbine to the total wind energy passing through the wind turbine rotating surface per unit time, λ wtFor tip speed ratio;

[0021] Photovoltaic power generation model:

[0022]

[0023] T s = T a + 0.0138 · (1 + 0.031T a ) · (1 - 0.042v) · G (6)

[0024] Wherein is the i th country comprehensive energy system t time photovoltaic power generation, G is the light intensity, T s is the photovoltaic cell surface temperature, P stc , G stc , T stc is the maximum output power, light intensity and photovoltaic cell surface temperature under standard test conditions, respectively, ε is the temperature coefficient of photovoltaic cell; T a is the ambient temperature; v is the wind speed;

[0025] Energy storage device model includes battery model and heat storage tank model, as follows:

[0026] Battery model:

[0027] Discharge

[0028]

[0029] Charging

[0030]

[0031] Wherein is the i th country comprehensive energy system t+1 time storage in the battery power, W t e,i is the i th country comprehensive energy system t time storage in the battery power, and is the i th country comprehensive energy system t time battery discharge power and charging power, respectively; and is the discharge efficiency and charging efficiency of the battery itself, respectively; and is the discharge loss and charging loss of the battery itself, respectively;

[0032] Heat storage tank model:

[0033] Heat release

[0034]

[0035] Charging

[0036]

[0037] wherein is the thermal energy stored in the thermal storage tank of the i-th national integrated energy system at time t+1, W t h,i is the thermal energy stored in the thermal storage tank of the i-th national integrated energy system at time t, W and are the discharging power and charging power of the thermal storage tank of the i-th national integrated energy system at time t, respectively; and are the discharging efficiency and charging efficiency of the thermal storage tank itself, respectively; and are the discharging loss and charging loss of the thermal storage tank itself, respectively;

[0038] The electric boiler model is specifically as follows:

[0039]

[0040] wherein, is the heating power of the electric boiler of the i-th national integrated energy system at time t, W is the electric power required by the electric boiler of the i-th national integrated energy system at time t, η EB is the conversion efficiency of the electric boiler.

[0041] S21, introducing a carbon trading model on the basis of the cross-border integrated energy system CBIES model described in step S1, specifically including a carbon trading cost, and specifically as follows:

[0042]

[0043] wherein is the carbon trading cost of the cross-border integrated energy system, N is the number of integrated energy systems, is the carbon emission of the i-th national integrated energy system, is the initial carbon quota of the i-th national integrated energy system, is the price of carbon emission right on the market; d is the length of the carbon emission interval; σ is the growth rate of the carbon emission right price of each step, and the carbon emission right price grows When the carbon trading cost is negative, it means that the carbon emission right is sold to obtain a profit.

[0044] The carbon emission of the i-th national integrated energy system is calculated as follows:

[0045]

[0046] wherein Direct carbon emissions of cross-border integrated energy system, Carbon emissions of purchased electricity, Carbon emissions of purchased heat;

[0047] Direct carbon emissions of cross-border integrated energy system The calculation is as follows:

[0048]

[0049] Where N is the number of integrated energy systems; T is the number of hours, taken as 24 hours; is the electric power output of the gas turbine at time t in the i th country integrated energy system, is the heat power output of the gas boiler at time t in the i th country integrated energy system; and are the power generation efficiency of the gas turbine and the heat generation efficiency of the gas boiler in the i th country integrated energy system, respectively;

[0050] Carbon emissions of purchased electricity The calculation is as follows:

[0051]

[0052] Where represents the purchased electricity consumption of the i th country integrated energy system, with units of ten thousand kilowatt-hours; is the grid baseline emission factor, which represents the indirect emissions per unit of electricity consumption, taken as 0.8367 in this paper, with units of tons of CO2 / kilowatt-hour;

[0053] Carbon emissions of purchased heat The calculation is as follows:

[0054]

[0055] Where represents the purchased heat consumption of the i th country integrated energy system, with units of million kilojoules; is the heat network baseline emission factor, which represents the indirect emissions per unit of heat consumption, taken as 0.096 in this paper, with units of tons of CO2 / million kilojoules;

[0056] Carbon initial quota of the i th country integrated energy system The calculation is as follows:

[0057]

[0058] Where is the power supply of the gas unit at time t in the i th country integrated energy system, with units of MWh, P e,bThe power supply reference value for the gas turbine unit is 0.3791, which is expressed in tCO2 / MWh. h This is the correction factor for the heat supply of the gas turbine unit, and the correction factor for the heat supply of the gas turbine unit is 1. The heat supplied by the gas turbine units of the integrated energy system of the i-th country is expressed in GJ and P. h,b The reference value for heating by the gas turbine unit is 0.0600, and the unit is tCO2 / GJ.

[0059] S22. Based on the CBIES (Cross-border Integrated Energy System) model described in step S1, a cross-border green certificate trading model is introduced, specifically including cross-border green certificate trading costs, as follows:

[0060]

[0061] in Let N be the cost of cross-border green certificate transactions for the integrated energy system of the i-th country, and N be the number of integrated energy systems. Let i be the actual amount of renewable energy consumed by the integrated energy system of the i-th country. Let i be the renewable energy quota of the integrated energy system of the i-th country. and The prices for buying and selling green certificates are λ, respectively. f This is the penalty coefficient;

[0062] The renewable energy quota for the i-th country's integrated energy system is calculated as follows:

[0063]

[0064] in For the i-th country's integrated energy system, the region is projected to absorb renewable energy in a given year. For the i-th country's integrated energy system, the predicted annual renewable energy access volume for the region is given. Let N be the estimated annual electricity consumption of the entire society in the region of the integrated energy system of the i-th country, and N be the number of integrated energy systems.

[0065] S23. In step S2, a collaborative operation optimization model for the cross-border integrated energy system considering carbon trading and cross-border green certificate trading is established, as follows:

[0066]

[0067] Where β is the conversion factor between cross-border green certificate trading and carbon trading rights, and α z The transaction price adjustment factor is determined by the type of renewable energy and the local gas price; The price of carbon emission rights in the market. The price for selling green certificates;

[0068] The cross-border integrated energy system collaborative operation optimization model considering carbon trading and cross-border green certificate trading includes system operation total cost, specifically as follows:

[0069]

[0070] Wherein C i , and are the total cost of the i-th country integrated energy system, fuel cost, energy purchase cost, energy sale revenue, system maintenance cost, interaction cost between integrated energy systems of different countries, carbon trading cost and green certificate trading cost, unit yuan;

[0071] The fuel cost of the i-th country integrated energy system is:

[0072]

[0073] Wherein c fuel is the fuel cost, unit yuan, represents the natural gas consumption of the gas turbine at t time, represents the natural gas consumption of the gas boiler at t time;

[0074] The energy purchase cost of the i-th country integrated energy system is:

[0075]

[0076] Wherein and are the electricity purchase cost and heat purchase cost of the integrated energy system of each country at t time, unit yuan; and are the purchased electricity power and purchased heat power of the i-th country integrated energy system at t time;

[0077] The energy sale revenue of the i-th country integrated energy system is:

[0078]

[0079] Wherein and are the electricity sale cost and heat sale cost, unit yuan; and are the sold electricity power and sold heat power of the i-th country integrated energy system at t time;

[0080] The system maintenance cost of the i-th country integrated energy system is:

[0081]

[0082] wherein c GT , c GB , c EB , c HS , c ES , c PV , c WT are the unit power maintenance costs of the gas turbine, the gas boiler, the electric boiler, the heat storage tank, the battery, the photovoltaic, and the wind turbine, respectively, unit yuan, is the output electric power of the gas turbine of the i-th national integrated energy system at time t, is the output heat power of the gas boiler of the i-th national integrated energy system at time t, is the heating power of the electric boiler of the i-th national integrated energy system at time t, is the discharging or charging power of the heat storage tank of the i-th national integrated energy system at time t, is the charging or discharging power of the battery of the i-th national integrated energy system at time t, is the photovoltaic power of the i-th national integrated energy system at time t, is the wind power of the i-th national integrated energy system at time t;

[0083] Interactions between national integrated energy systems:

[0084]

[0085] wherein is the electricity purchase cost between the integrated energy systems, is the heat purchase cost between the integrated energy systems, is the electric power transmitted between the i-th and j-th national integrated energy systems, is the heat power transmitted between the i-th and j-th national integrated energy systems;

[0086] The constraint conditions of the cross-border integrated energy system collaborative operation optimization model considering carbon trading and cross-border green certificate trading specifically include the internal power constraints of the national integrated energy systems, the interaction power constraints between the national integrated energy systems and the power grid, the interaction power constraints between the national integrated energy systems and the heat grid, the interaction power constraints between the power grids of the national integrated energy systems, and the interaction power constraints between the heat grids of the national integrated energy systems, and are specifically as follows:

[0087] The internal power constraints of the national integrated energy systems specifically include the electric power balance constraint, the heat power balance constraint, the upper and lower limits of the output of the devices of the national integrated energy systems, the power constraints of the batteries, and the power constraints of the heat storage tanks, and are specifically as follows:

[0088] Electric power balance constraint:

[0089]

[0090] wherein, is the purchased electricity power of the i-th country comprehensive energy system at time t interacting with the power grid, is the sold electricity power of the i-th country comprehensive energy system at time t interacting with the power grid, is the output electricity power of the gas turbine of the i-th country comprehensive energy system at time t, is the electricity power loss of the power grid between the i-th and j-th country comprehensive energy systems, is the electricity power transmitted between the i-th and j-th country comprehensive energy systems, is the electricity power required by the electric boiler of the i-th country comprehensive energy system at time t, are respectively the charging power and discharging power of the battery of the i-th country comprehensive energy system at time t, is the photovoltaic power of the i-th country comprehensive energy system at time t, is the wind power of the i-th country comprehensive energy system at time t, is the electricity load of the i-th country comprehensive energy system;

[0091] thermal power balance constraint:

[0092]

[0093] wherein is the purchased heat power of the i-th country comprehensive energy system at time t interacting with the heat grid, is the sold heat power of the i-th country comprehensive energy system at time t interacting with the heat grid, is the output heat power of the gas boiler of the i-th country comprehensive energy system at time t, is the heating power of the electric boiler of the i-th country comprehensive energy system at time t, is the heat power loss of the heat grid between the i-th and j-th country comprehensive energy systems, is the heat power transmitted between the i-th and j-th country comprehensive energy systems, are respectively the discharging power and charging power of the heat storage tank of the i-th country comprehensive energy system at time t, is the heat load of the i-th country;

[0094] upper and lower limit constraints of the output power of the devices of the country comprehensive energy systems:

[0095]

[0096] wherein is the electricity power of the device m in the i-th country comprehensive energy system, are the upper and lower limits of the electricity power of the device m in the i-th country comprehensive energy system; is the electricity power of the device m, The upper and lower limits of the thermal power of the equipment m in the i-th national comprehensive energy system;

[0097] Battery power constraint:

[0098]

[0099]

[0100]

[0101]

[0102] The battery capacity of the i-th national comprehensive energy system, The maximum charging rate, The maximum discharging rate, The battery charging power and discharging power of the i-th national comprehensive energy system at time t, respectively; The state bit of charging at time t, The state bit of discharging at time t, is a 0-1 variable, indicating that the charging and discharging state of the same equipment at the same time is unique, The maximum and minimum energy storage of the battery of the i-th national comprehensive energy system;

[0103] Heat storage tank power constraint:

[0104]

[0105]

[0106]

[0107]

[0108] Wherein The heat storage tank capacity of the i-th national comprehensive energy system, The maximum charging rate, The maximum discharging rate, The heat storage tank discharging power and charging power of the i-th national comprehensive energy system at time t, respectively; The maximum and minimum heat storage of the heat storage tank of the i-th national comprehensive energy system;

[0109] The interaction power constraint between the national comprehensive energy system and the electric and heat network is as follows:

[0110] The interaction power constraint between the national comprehensive energy system and the power grid is as follows:

[0111]

[0112] Wherein is the maximum allowed purchased power of the i-th country integrated energy system from the power grid, is the maximum allowed sold power of the i-th country integrated energy system to the power grid, is the purchased power of the i-th country integrated energy system from the power grid at time t, is the sold power of the i-th country integrated energy system to the power grid at time t;

[0113] Power exchange constraints between the integrated energy systems and the heat grid of each country:

[0114]

[0115] wherein is the maximum allowed purchased power of the i-th country integrated energy system from the heat grid, is the maximum allowed sold power of the i-th country integrated energy system to the heat grid, is the purchased heat power of the i-th country integrated energy system from the heat grid at time t, is the sold heat power of the i-th country integrated energy system to the heat grid at time t;

[0116] Power exchange constraints between the integrated energy systems and the heat grid of each country:

[0117] Power exchange constraints between the integrated energy systems and the power grid of each country:

[0118]

[0119] wherein is the maximum value of the transmitted electric power between the i-th and j-th country integrated energy systems, is the maximum value of the transmitted heat power between the i-th and j-th country integrated energy systems, is the transmitted electric power between the i-th and j-th country integrated energy systems, is the transmitted heat power between the i-th and j-th country integrated energy systems;

[0120] Power exchange constraints between the integrated energy systems and the heat grid of each country:

[0121]

[0122] wherein is the maximum value of the electric power loss of the power grid between the i-th and j-th country integrated energy systems, is the maximum value of the heat power loss of the heat grid between the i-th and j-th country integrated energy systems, is the electric power loss of the power grid between the i-th and j-th country integrated energy systems, is the heat power loss of the heat grid between the i-th and j-th country integrated energy systems;

[0123] The solver is a CPLEX solver or a GUROBI solver.

[0124] The application discloses a cross-border integrated energy system collaborative operation optimization method considering carbon-green certificate combined transaction, which comprises the following steps. DETAILED DESCRIPTION

[0125] In order to make the purpose, technical scheme and advantages of the present application clearer, the present application is described in detail below in combination with the drawings and examples. It should be understood that the specific examples described herein are only used to explain the present application and do not limit the present application.

[0126] The application provides a cross-border integrated energy system collaborative operation optimization method considering carbon-green certificate combined transaction, which comprises the following steps.

[0127] S1, a cross-border integrated energy system CBIES model is constructed, the cross-border integrated energy system CBIES comprises a gas turbine, a gas boiler, a renewable energy unit, energy storage equipment and an electric boiler, wherein the renewable energy unit comprises wind power generation and photovoltaic power generation, and the energy storage equipment comprises a battery and a heat storage tank;

[0128] S2, a cross-border green certificate transaction model and a carbon transaction model are introduced on the basis of the cross-border integrated energy system CBIES model in step S1, and a cross-border integrated energy system collaborative operation optimization model considering carbon transaction and cross-border green certificate transaction is established;

[0129] S3, the cross-border integrated energy system collaborative operation optimization model considering carbon transaction and cross-border green certificate transaction established in step S2 is solved by using a solver to obtain an optimization result.

[0130] S11, the cross-border integrated energy system CBIES model constructed in step S1 specifically comprises a gas turbine model, a gas boiler model, a renewable energy unit model, an energy storage equipment model and an electric boiler model, and specifically comprises the following:

[0131] The gas turbine model specifically comprises the following:

[0132]

[0133]

[0134] wherein represents the output electric power of the gas turbine of the comprehensive energy system of the i-th country at time t, represents the output thermal power of the gas turbine of the comprehensive energy system of the i-th country at time t, L NG represents the low calorific value of natural gas, represents the natural gas consumption of the gas turbine at time t, η GT represents the power generation efficiency of the gas turbine, N GT represents the waste heat recovery coefficient;

[0135] The gas boiler model is specifically as follows:

[0136]

[0137] wherein represents the output thermal power of the gas boiler of the comprehensive energy system of the i-th country at time t, L NG represents the low calorific value of natural gas, represents the natural gas consumption of the gas boiler at time t, η GB represents the heat generation efficiency of the gas boiler;

[0138] The renewable energy unit model includes a wind power generation model and a photovoltaic power generation model, and is specifically as follows:

[0139] Wind power generation model:

[0140]

[0141] wherein is the wind power of the comprehensive energy system of the i-th country at time t, ρ is the air density, A is the swept area of the wind turbine blade, v is the wind speed, c wt is the wind energy utilization coefficient, which is the ratio of the wind energy absorbed by the wind turbine to the total wind energy passing through the wind turbine rotating surface per unit time, λ wt is the tip speed ratio;

[0142] Photovoltaic power generation model:

[0143]

[0144] T s = T a + 0.0138 · (1 + 0.031 T a ) · (1 - 0.042 v) · G (47)

[0145] wherein is the photovoltaic power of the i th country comprehensive energy system at time t, G is the light intensity, T s is the photovoltaic cell surface temperature, P stc , G stc , T stc are the maximum output power, light intensity and photovoltaic cell surface temperature under standard test conditions respectively, and ε is the temperature coefficient of the photovoltaic cell; Ta is the ambient temperature; and v is the wind speed;

[0146] The energy storage device model includes a battery model and a heat storage tank model, and is specifically as follows:

[0147] Battery model:

[0148] Discharge

[0149]

[0150] Charge

[0151]

[0152] wherein is the amount of electricity stored in the battery of the i th country comprehensive energy system at time t+1, is the amount of electricity stored in the battery of the i th country comprehensive energy system at time t, and are the discharge power and charge power of the battery of the i th country comprehensive energy system at time t, respectively; and are the discharge efficiency and charge efficiency of the battery itself, respectively; and are the discharge loss and charge loss of the battery itself, respectively;

[0153] Heat storage tank model:

[0154] Heat release

[0155]

[0156] Heat charging

[0157]

[0158] wherein is the amount of heat stored in the heat storage tank of the i th country comprehensive energy system at time t+1, W t h,i is the amount of heat stored in the heat storage tank of the i th country comprehensive energy system at time t, and are the heat release power and heat charging power of the heat storage tank of the i th country comprehensive energy system at time t, respectively; and These are the heat release efficiency and heat charging efficiency of the heat storage tank itself, respectively. and These are the heat release loss and heat charging loss of the heat storage tank itself, respectively.

[0159] The electric boiler model is as follows:

[0160]

[0161] in, Let be the heating capacity of the electric boiler in the integrated energy system of the i-th country at time t; Let η be the electrical power required by the electric boiler in the integrated energy system of the i-th country at time t. EB This refers to the conversion efficiency of the electric boiler.

[0162] S21. Based on the CBIES (Cross-border Integrated Energy System) model described in step S1, a carbon trading model is introduced, specifically including carbon trading costs, as follows:

[0163]

[0164] in The carbon trading cost of a cross-border integrated energy system is given by N, where N is the number of integrated energy systems. Let i be the carbon emissions of the integrated energy system of country i. For the i-th country's integrated energy system initial carbon allowance, d represents the market price of carbon emission rights; d is the length of the carbon emission range; σ is the price increase for each tier of carbon emission rights, with the price increasing for each tier. When carbon trading costs are negative, it means that the sale of carbon emission rights generates revenue.

[0165] Carbon emissions from the integrated energy system of country i The calculation is as follows:

[0166]

[0167] in For direct carbon emissions from cross-border integrated energy systems, Carbon emissions from purchased electricity Carbon emissions from purchased heat;

[0168] Direct carbon emissions from cross-border integrated energy systems The calculation is as follows:

[0169]

[0170] Where N is the number of integrated energy systems; T is the number of hours, taken as 24 hours; is the electric power output of the gas turbine in the i th country comprehensive energy system at time t, is the heat power output of the gas boiler in the i th country comprehensive energy system at time t; and are the power generation efficiency of the gas turbine and the heat generation efficiency of the gas boiler in the i th country comprehensive energy system, respectively;

[0171] Carbon emissions of purchased electricity is calculated as follows:

[0172]

[0173] wherein represents the consumption of purchased electricity in the i th country comprehensive energy system, in units of ten thousand kilowatt-hours; is the grid baseline emission factor, representing the indirect emissions per unit of electricity consumption, and is taken as 0.8367 in this paper, in units of tons of CO2 / kilowatt-hour;

[0174] Carbon emissions of purchased heat is calculated as follows:

[0175]

[0176] wherein represents the consumption of purchased heat in the i th country comprehensive energy system, in units of million kilojoules; is the heat network baseline emission factor, representing the indirect emissions per unit of heat consumption, and is taken as 0.096 in this paper, in units of tons of CO2 / million kilojoules;

[0177] Carbon initial quota of the i th country comprehensive energy system is calculated as follows:

[0178]

[0179] wherein is the power supply of the gas unit in the i th country comprehensive energy system at time t, in units of MWh, P e,b is the power supply baseline value of the gas unit, taken as 0.3791 in this paper, in units of tCO2 / MWh, γ h is the heat supply correction coefficient of the gas unit, and the heat supply correction coefficient of the gas unit is 1; is the heat supply of the gas unit in the i th country comprehensive energy system, in units of GJ, P h,b is the heat supply baseline value of the gas unit, taken as 0.0600 in this paper, in units of tCO2 / GJ;

[0180] S22, introducing a cross-border green certificate transaction model based on the cross-border comprehensive energy system CBIES model described in step S1, specifically including cross-border green certificate transaction cost, specifically as follows:

[0181]

[0182] wherein is the cross-border green certificate transaction cost of the i th country comprehensive energy system, and N is the number of comprehensive energy systems, is the actual renewable energy consumption of the i th country comprehensive energy system, is the renewable energy quota of the i th country comprehensive energy system, and are the prices of buying and selling green certificates, respectively, and λ f is the penalty coefficient;

[0183] The renewable energy quota of the i th country comprehensive energy system is calculated as follows:

[0184]

[0185] wherein is the regional predicted annual renewable energy consumption of the i th country comprehensive energy system, is the regional predicted annual renewable energy access of the i th country comprehensive energy system, is the regional predicted annual electricity consumption of the i th country comprehensive energy system, and N is the number of comprehensive energy systems;

[0186] S23, the step S2 of establishing a cross-border comprehensive energy system collaborative operation optimization model considering carbon trading and cross-border green certificate trading, specifically as follows:

[0187]

[0188] wherein β is the cross-border green certificate transaction and carbon trading right conversion coefficient, and α z is the transaction price adjustment coefficient, determined by the type of renewable energy and the local gas price; is the market carbon emission right price, is the price of selling green certificates;

[0189] The cross-border comprehensive energy system collaborative operation optimization model considering carbon trading and cross-border green certificate trading includes the total system operation cost, specifically as follows:

[0190]

[0191] wherein C i , and The total cost of the i-th country's comprehensive energy system, fuel cost, energy purchase cost, energy sale revenue, system maintenance cost, interaction cost between comprehensive energy systems, carbon trading cost and green certificate trading cost, unit yuan;

[0192] The fuel cost of the i-th country's comprehensive energy system:

[0193]

[0194] Wherein c fuel is the fuel cost, unit yuan, represents the natural gas consumption of the gas turbine at time t, represents the natural gas consumption of the gas boiler at time t;

[0195] The energy purchase cost of the i-th country's comprehensive energy system:

[0196]

[0197] Wherein and are the electricity purchase cost and heat purchase cost of the i-th country's comprehensive energy system at time t, unit yuan; and are the electricity purchase power and heat purchase power of the i-th country's comprehensive energy system at time t, unit MW;

[0198] The energy sale revenue of the i-th country's comprehensive energy system:

[0199]

[0200] Wherein and are the electricity sale cost and heat sale cost, unit yuan; and are the electricity sale power and heat sale power of the i-th country's comprehensive energy system at time t;

[0201] The system maintenance cost of the i-th country's comprehensive energy system:

[0202]

[0203] Wherein c GT , c GB , c EB , c HS , c ES , c PV , c WT are the unit power maintenance costs of the gas turbine, gas boiler, electric boiler, heat storage tank, battery, photovoltaic and wind turbine, unit yuan, is the output electric power of the gas turbine of the i th national integrated energy system at time t, is the output heat power of the gas boiler of the i th national integrated energy system at time t, is the heating power of the electric boiler of the i th national integrated energy system at time t, is the heat release power or heat charging power of the heat storage tank of the i th national integrated energy system at time t, is the charging power or discharging power of the battery of the i th national integrated energy system at time t, is the photovoltaic power generation power of the i th national integrated energy system at time t, is the wind power of the i th national integrated energy system at time t;

[0204] Interactions between national integrated energy systems:

[0205]

[0206] wherein He is the electricity purchase cost and heat purchase cost between integrated energy systems, is the electric power transmitted between the i th and j th national integrated energy systems, is the heat power transmitted between the i th and j th national integrated energy systems;

[0207] The constraint conditions of the cross-border integrated energy system coordinated operation optimization model considering carbon trading and cross-border green certificate trading specifically include the internal power constraints of the national integrated energy system, the interaction power constraints between the national integrated energy system and the power grid, the interaction power constraints between the national integrated energy system and the heat grid, the interaction power constraints between the power grids of the national integrated energy systems, and the interaction power constraints between the heat grids of the national integrated energy systems, and are specifically as follows:

[0208] The internal power constraints of the national integrated energy system specifically include electric power balance constraints, heat power balance constraints, upper and lower limits of the output of the devices of the national integrated energy system, battery power constraints, and heat storage tank power constraints, and are specifically as follows:

[0209] Electric power balance constraints:

[0210]

[0211] wherein, is the electricity purchase power of the i th national integrated energy system interacting with the power grid at time t, is the electricity sale power of the i th national integrated energy system interacting with the power grid at time t, is the output electric power of the gas turbine of the i th national integrated energy system at time t, is the electric power loss of the power grid transmitted between the i th and j th national integrated energy systems, Pij(t) is the electric power transmitted between the two countries i, j, Pbi(t) is the electric power required by the electric boiler of the i-th country at time t, Pbi(t) is the electric power required by the electric boiler of the i-th country at time t, Ppv(t) is the photovoltaic power of the i-th country at time t, Pw(t) is the wind power of the i-th country at time t, L i E Pdi(t) is the electric load of the i-th country comprehensive energy system;

[0212] Thermal power balance constraints:

[0213]

[0214] Pbi(t) is the electric power required by the electric boiler of the i-th country at time t, Pbi(t) is the electric power required by the electric boiler of the i-th country at time t, Pbi(t) is the electric power required by the electric boiler of the i-th country at time t, Pbi(t) is the electric power required by the electric boiler of the i-th country at time t, Pbi(t) is the electric power required by the electric boiler of the i-th country at time t, Pij(t) is the electric power transmitted between the two countries i, j, Pij(t) is the electric power transmitted between the two countries i, j, Pbi(t) is the electric power required by the electric boiler of the i-th country at time t, Pdi(t) is the electric load of the i-th country comprehensive energy system;

[0215] Upper and lower limits of the output of the equipment of the country comprehensive energy system:

[0216]

[0217] Pbi(t) is the electric power required by the electric boiler of the i-th country at time t, Pbi(t) is the electric power required by the electric boiler of the i-th country at time t, Pbi(t) is the electric power required by the electric boiler of the i-th country at time t, Pbi(t) is the electric power required by the electric boiler of the i-th country at time t, Pbi(t) is the electric power required by the electric boiler of the i-th country at time t,

[0218] Battery power constraints:

[0219]

[0220]

[0221]

[0222]

[0223] is the maximum charging rate of the i-th country comprehensive energy system, is the maximum discharging rate of the i-th country comprehensive energy system, respectively represent the charging and discharging power of the i-th country comprehensive energy system at time t; is the state bit of charging at time t, is the state bit of discharging at time t, is a 0-1 variable, indicating that the charging and discharging state of the same device at the same time is unique, is the maximum and minimum energy storage of the i-th country comprehensive energy system battery,

[0224] Regulation of regenerative tank power:

[0225]

[0226]

[0227]

[0228]

[0229] wherein is the maximum charging rate of the i-th country comprehensive energy system, is the maximum discharging rate of the i-th country comprehensive energy system, respectively represent the charging and discharging power of the i-th country comprehensive energy system at time t; is the maximum and minimum energy storage of the i-th country comprehensive energy system battery,

[0230] The interaction power constraints between the comprehensive energy system of each country and the electricity and heat network are as follows:

[0231] The interaction power constraints between the comprehensive energy system of each country and the electricity network are as follows:

[0232]

[0233] wherein is the maximum allowed power purchase of the i-th country comprehensive energy system and the electricity network, is the maximum allowed power sale of the i-th country comprehensive energy system and the electricity network, is the power purchase of the i-th country comprehensive energy system and the electricity network at time t, is the power sale of the i-th country comprehensive energy system and the electricity network at time t;​​

[0234] The maximum allowed power exchange between the i-th country's integrated energy system and the heat network,

[0235]

[0236] wherein is the maximum allowed power purchase of the i-th country's integrated energy system and the heat network, is the maximum allowed power sale of the i-th country's integrated energy system and the heat network, is the heat purchase of the i-th country's integrated energy system and the heat network at time t, is the heat sale of the i-th country's integrated energy system and the heat network at time t.

[0237] The power exchange constraints between the integrated energy systems of different countries via the power grid and the heat network are as follows:

[0238] The power exchange constraints between the integrated energy systems of different countries via the power grid are as follows:

[0239]

[0240] wherein is the maximum value of the power transmission between the integrated energy systems of the i-th and j-th countries, is the maximum value of the heat transmission between the integrated energy systems of the i-th and j-th countries, is the power transmission between the integrated energy systems of the i-th and j-th countries, is the heat transmission between the integrated energy systems of the i-th and j-th countries.

[0241] The power exchange constraints between the integrated energy systems of different countries via the heat network are as follows:

[0242]

[0243] wherein is the maximum value of the power loss of the power grid between the integrated energy systems of the i-th and j-th countries, is the maximum value of the heat loss of the heat network between the integrated energy systems of the i-th and j-th countries, is the power loss of the power grid between the integrated energy systems of the i-th and j-th countries, is the heat loss of the heat network between the integrated energy systems of the i-th and j-th countries.

[0244] The solver is a CPLEX solver or a GUROBI solver.

[0245] Embodiment

[0246] The cross-border integrated energy system CBIES interacts through the power grid and the heating pipe network, and supplies electric energy and heat energy to the loads of the respective national integrated energy systems; the cross-border integrated energy system is equipped with a micro gas turbine, a gas boiler, a gas turbine, an electric boiler, electricity storage, heat storage, photovoltaic, a fan and a waste heat recovery device, and different countries set different parameters, and the model built is simulated on a computer equipped with an i7 CPU and 16 GB RAM, and the simulation platform is Matlab2020b.

[0247] (1) Parameters and simulation settings of the cross-border integrated energy system

[0248] Table 1 gives the device capacity and parameters used in the simulation, and Table 2 sets six groups of cases for the cross-border integrated energy system operation scenarios. The first three groups of cases do not consider cross-border interaction of the energy system, and the last three groups of cases consider cross-border interaction of the energy system, and whether to introduce carbon and green certificate trading is considered to compare the differences in the planning results of the cross-border integrated energy system.

[0249] Table 1 Device capacity and parameters

[0250]

[0251] Table 2 Cross-border integrated energy system operation scenarios

[0252]

[0253] The three national integrated energy systems in case1, case2 and case3 run independently without energy transmission between systems through the electric and heat networks; the three national integrated energy systems in case4, case5 and case6 run cooperatively with energy transmission between systems through the electric and heat networks, and when energy is transmitted between systems, the price of purchased energy is less than the price of purchased electricity from the power grid, and both ways aim to minimize the total operating cost of the three integrated energy systems.

[0254] (2) Optimization analysis

[0255] The optimization results of the six groups of cases are shown in Tables 3 and 4.

[0256] Table 3 Operating cost of each energy system without considering cross-border interaction

[0257]

[0258] As can be seen from Table 3, compared with Case1, the operating costs of the integrated energy systems of each country in Case2 considering green certificate trading are reduced by 57680, 50260 and 41660 yuan respectively; compared with Case2, the operating costs of the integrated energy systems of each country in Case3 considering green certificate and carbon trading are reduced by 46440, 43970 and 43160 yuan respectively.

[0259] Table 4 considers the operation cost of each energy system when cross-border interaction

[0260]

[0261] As shown in Table 4, compared with Case 4, the operation cost of the comprehensive energy system of each country in Case 5 considering green certificate transaction is respectively reduced by 57640, 49560 and 52890 yuan; compared with Case 5, the operation cost of the comprehensive energy system of each country in Case 6 considering green certificate and carbon transaction is respectively reduced by 43230, 42360 and 45623 yuan.

[0262] From the comparison of Case 1, 4, Case 2, 5, Case 3, 6, it can be seen that after considering the cross-border interaction of the system, the operation cost of the comprehensive energy system of A country is respectively reduced by 50580, 65830 and 66410 yuan, the operation cost of the comprehensive energy system of B country is respectively reduced by 50540, 66120 and 77640 yuan, and the operation cost of the comprehensive energy system of C country is respectively reduced by 47330, 64510 and 80103 yuan.

[0263] Considering the influence of carbon transaction and green certificate on the annual operation cost and wind-solar consumption of each country, the following important conclusions can be obtained:

[0264] a) The cross-border comprehensive energy system with carbon and green certificate transaction mechanism in the cross-border region can effectively reduce the operation cost of the system, guide each energy system to reduce its carbon emission, and make the cross-border comprehensive energy system give priority to renewable energy consumption for maximizing the income;

[0265] b) Considering the cross-border interaction between energy systems can regulate the imbalance of resource allocation between countries, effectively reduce the power purchase quantity of the grid, increase the proportion of natural gas and renewable energy consumption, and has important significance for reducing the carbon emission of system units, which is an important way to realize the collaborative operation of the cross-border comprehensive energy system.

[0266] In summary, the simulation results verify the effectiveness of the proposed carbon-green certificate optimization planning model in the cross-border comprehensive energy system CBIES, and the results show that the setting of the carbon-green certificate optimization planning model in the cross-border comprehensive energy system CBIES only brings a little investment cost, but brings considerable income to the system operation, therefore, the cross-border comprehensive energy system collaborative operation optimization method considering carbon-green certificate joint transaction provided by the application has potential engineering application value.

[0267] The above description of the embodiments is only used to help understand the method of the present application and its core idea; meanwhile, for those skilled in the art, according to the idea of the present application, there will be changes in the specific implementation and application range, and the above description should not be understood as a limitation on the present application.

Claims

1. A method for optimizing the coordinated operation of a cross-border integrated energy system considering carbon-green certificate joint trading, characterized in that, Includes the following steps: S1. Construct a cross-border integrated energy system (CBIES) model. The cross-border integrated energy system (CBIES) includes gas turbines, gas boilers, renewable energy units, energy storage equipment, and electric boilers. The renewable energy units include wind power generation and photovoltaic power generation, and the energy storage equipment includes batteries and thermal storage tanks. S2. Based on the CBIES model of the cross-border integrated energy system described in step S1, introduce the cross-border green certificate trading model and the carbon trading model to establish a collaborative operation optimization model of the cross-border integrated energy system that considers carbon trading and cross-border green certificate trading. S21. Based on the CBIES (Cross-border Integrated Energy System) model described in step S1, a carbon trading model is introduced, specifically including carbon trading costs, as follows: in The carbon trading cost of a cross-border integrated energy system is given by N, where N is the number of integrated energy systems. Let i be the carbon emissions of the integrated energy system of country i. For the i-th country's integrated energy system initial carbon allowance, d represents the market price of carbon emission rights; d is the length of the carbon emission range; σ is the price increase for each tier of carbon emission rights, with the price increasing for each tier. When carbon trading costs are negative, it means that the sale of carbon emission rights generates revenue. Carbon emissions from the integrated energy system of country i The calculation is as follows: in For direct carbon emissions from cross-border integrated energy systems, Carbon emissions from purchased electricity Carbon emissions from purchased heat; Direct carbon emissions from cross-border integrated energy systems The calculation is as follows: Where N is the number of integrated energy systems; T is the number of hours, taken as 24 hours; It is the electrical power output of the gas turbine at time t in the integrated energy system of the i-th country. It is the thermal power output of the gas-fired boiler in the i-th national integrated energy system at time t; and These are the power generation efficiency of the gas turbine and the heating efficiency of the gas boiler in the integrated energy system of the i-th country, respectively. Carbon emissions from purchased electricity The calculation is as follows: in This represents the amount of electricity purchased from outside the integrated energy system of the i-th country, in ten thousand kilowatt-hours. The baseline emission factor for the power grid represents the indirect emissions per unit of electricity consumed. In this paper, the value is taken as 0.8367, and the unit is tons of CO2 per kilowatt-hour. Carbon emissions from purchased heat The calculation is as follows: in This represents the amount of purchased heat consumed by the integrated energy system of the i-th country, in millions of kilojoules. The baseline emission factor for the heating network represents the indirect emissions per unit of heat consumed. In this paper, it is taken as 0.096, and the unit is tons of CO2 per million kilojoules. Initial carbon quota for the integrated energy system of the i-th country The calculation is as follows: in The power supply of the gas turbine unit in the integrated energy system of the i-th country at time t, in MWh, P e,b The power supply reference value for the gas turbine unit is 0.3791, which is expressed in tCO2 / MWh. h This is the correction factor for the heat supply of the gas turbine unit, and the correction factor for the heat supply of the gas turbine unit is 1. The heat supplied by the gas turbine units of the integrated energy system of the i-th country is expressed in GJ and P. h,b The reference value for heating by the gas turbine unit is 0.0600, and the unit is tCO2 / GJ. S22. Based on the CBIES (Cross-border Integrated Energy System) model described in step S1, a cross-border green certificate trading model is introduced, specifically including cross-border green certificate trading costs, as follows: in Let N be the cost of cross-border green certificate transactions for the integrated energy system of the i-th country, and N be the number of integrated energy systems. Let i be the actual amount of renewable energy consumed by the integrated energy system of the i-th country. Let i be the renewable energy quota of the integrated energy system of the i-th country. and The prices for buying and selling green certificates are λ, respectively. f This is the penalty coefficient; The renewable energy quota for the i-th country's integrated energy system is calculated as follows: in For the i-th country's integrated energy system, the region is projected to absorb renewable energy in a given year. For the i-th country's integrated energy system, the predicted annual renewable energy access volume for the region is given. Let N be the estimated annual electricity consumption of the entire society in the region of the integrated energy system of the i-th country, and N be the number of integrated energy systems. S23. In step S2, a collaborative operation optimization model for the cross-border integrated energy system considering carbon trading and cross-border green certificate trading is established, as follows: Where β is the conversion factor between cross-border green certificate trading and carbon trading rights, and α z The transaction price adjustment factor is determined by the type of renewable energy and the local gas price; The price of carbon emission rights in the market. The price for selling green certificates; The cross-border integrated energy system collaborative operation optimization model considering carbon trading and cross-border green certificate trading includes the total system operating cost, as detailed below: Where C i , and These represent the total operating cost of the integrated energy system of country i, fuel cost, energy purchase cost, energy sales revenue, system maintenance cost, interaction cost between integrated energy systems of various countries, carbon trading cost, and green certificate trading cost, respectively, in yuan; Fuel cost of operating the integrated energy system of country i: Where c fuel Fuel cost, in yuan. This represents the amount of natural gas consumed by the gas turbine at time t. This represents the amount of natural gas consumed by the gas-fired boiler at time t. Energy purchase cost for the integrated energy system operation of country i: in and These represent the costs of purchasing electricity and heat from the grid by each country's integrated energy system at time t, in yuan. and These represent the electricity purchased by the integrated energy system of the i-th country at time t and the heat purchased by the integrated energy system with the grid. Revenue from energy sales during the operation of the integrated energy system in country i: in and These are the electricity sales fees and heat sales fees, respectively, in yuan; and These represent the electricity sales power and heat sales power of the integrated energy system of the i-th country at time t, which interact with the power grid and the heat network, respectively. The operation and maintenance cost of the integrated energy system of the i-th country: Where c GT c GB c EB c HS c ES c PV c WT The unit power maintenance cost for gas turbines, gas boilers, electric boilers, thermal storage tanks, batteries, photovoltaic systems, and wind turbines are listed below, in yuan. Let be the output electrical power of the gas turbine in the integrated energy system of the i-th country at time t. Let be the output thermal power of the gas-fired boiler in the integrated energy system of the i-th country at time t. Let be the heating capacity of the electric boiler in the integrated energy system of the i-th country at time t. Let t be the heat release or charging power of the thermal storage tank in the integrated energy system of the i-th country at time t. The charging or discharging power of the battery in the integrated energy system of the i-th country at time t. Let be the photovoltaic power generation capacity of the integrated energy system of the i-th country at time t. Let be the wind power output of the integrated energy system of the i-th country at time t; Interaction costs between national integrated energy systems: in what These are the costs for purchasing electricity and heat between integrated energy systems. The electrical power transmitted between the integrated energy systems of countries i and j The heat power transferred between the integrated energy systems of countries i and j; The constraints of the cross-border integrated energy system collaborative operation optimization model considering carbon trading and cross-border green certificate trading specifically include the power constraints within each country's integrated energy system, the power constraints between each country's integrated energy system and the power grid, the power constraints between each country's integrated energy system and the heating network, the power constraints between each country's integrated energy system and the power grid, and the power constraints between each country's integrated energy system and the heating network, as detailed below: The power constraints within the integrated energy systems of various countries specifically include electrical power balance constraints, thermal power balance constraints, upper and lower limits of equipment output constraints in the integrated energy systems of various countries, battery power constraints, and thermal storage tank power constraints, as detailed below: Electric power balance constraints: in, Let be the power purchased by the integrated energy system of country i at time t when interacting with the power grid. Let be the electricity sales power of the integrated energy system of country i interacting with the power grid at time t. Let be the output electrical power of the gas turbine in the integrated energy system of the i-th country at time t. The power loss during power transmission between the integrated energy systems of countries i and j. The electrical power transmitted between the integrated energy systems of countries i and j Let be the electrical power required by the electric boiler in the integrated energy system of the i-th country at time t. These represent the charging power and discharging power of the battery in the integrated energy system of the i-th country at time t, respectively. Let be the photovoltaic power generation capacity of the integrated energy system of the i-th country at time t. Let be the wind power output of the integrated energy system of the i-th country at time t. Let i be the electrical load of the integrated energy system of the i-th country; Thermal power balance constraint: in Let t be the heat purchase power of the i-th country's integrated energy system interacting with the heating network at time t. Let be the heat sales power of the i-th country's integrated energy system interacting with the heating network at time t. Let be the output thermal power of the gas-fired boiler in the integrated energy system of the i-th country at time t. Let be the heating capacity of the electric boiler in the integrated energy system of the i-th country at time t. The heat power loss in the heat network transmission between the integrated energy systems of countries i and j For the heat power transferred between the integrated energy systems of countries i and j, These represent the heat release power and heat charge power of the thermal storage tank at time t in the integrated energy system of the i-th country. Let the heat load be that of the i-th country. Constraints on the upper and lower limits of integrated energy system equipment output in various countries: in Let m be the electrical power of equipment m in the integrated energy system of the i-th country. Let m be the upper and lower limits of the electrical power of equipment m in the integrated energy system of the i-th country; Let m be the electrical power of the device. Let m be the upper and lower limits of the thermal power of equipment m in the integrated energy system of the i-th country. Battery power constraints: For the battery capacity of the integrated energy system of the i-th country, Maximum charging rate, This is the maximum discharge rate. These represent the charging power and discharging power of the battery in the integrated energy system of the i-th country at time t, respectively. The state bit for charging at time t. This is the state bit for energy release at time t. It is a 0-1 variable, indicating that the charging / discharging state of the same device at the same time is unique. Let the maximum and minimum energy storage capacity of the battery in the integrated energy system of the i-th country be denoted by . Thermal storage tank power constraints: in Let the thermal storage tank capacity of the integrated energy system of the i-th country be , For the maximum heat charging ratio, For maximum heat release rate, These represent the heat release power and heat charge power of the thermal storage tank at time t in the integrated energy system of the i-th country. Let be the maximum and minimum heat storage capacity of the thermal storage tank in the integrated energy system of the i-th country; The power constraints of the interaction between the integrated energy systems and electricity and heat networks in various countries are as follows: Power constraints between integrated energy systems and power grids in various countries: in Let i be the maximum allowable power purchase capacity for the integrated energy system and power grid interaction of the i-th country. Let i be the maximum permissible power sales capacity for interaction between the integrated energy system and the power grid in the i-th country. Let be the power purchased by the integrated energy system of country i at time t when interacting with the power grid. Let be the electricity sales power of the integrated energy system of the i-th country at time t, which interacts with the power grid. Interactive power constraints between integrated energy systems and heating networks in various countries: in The maximum permissible power purchase capacity for the interaction between the integrated energy system and the heating network of the i-th country. Let i be the maximum allowable power output for interaction between the integrated energy system and the heating network of the i-th country. Let t be the heat purchase power of the i-th country's integrated energy system interacting with the heating network at time t. Let be the heat sales power of the i-th national integrated energy system interacting with the heating network at time t; The power constraints of electricity and heat network interactions between national integrated energy systems are as follows: Power constraints of grid interaction between national integrated energy systems: in This represents the maximum electrical power transmitted between the integrated energy systems of countries i and j. This represents the maximum heat transfer capacity between the integrated energy systems of countries i and j. The electrical power transmitted between the integrated energy systems of countries i and j The heat power transferred between the integrated energy systems of countries i and j; Inter-regional power constraints of heat networks in national integrated energy systems: in This represents the maximum power loss during power transmission between the integrated energy systems of countries i and j. This represents the maximum heat power loss transmitted between the integrated energy systems of countries i and j via the heat network. The power loss during power transmission between the integrated energy systems of countries i and j. The heat power loss transmitted between the integrated energy systems of countries i and j in the heat network; S3. Based on the cross-border integrated energy system collaborative operation optimization model considering carbon trading and cross-border green certificate trading established in step S2, use the solver to obtain the optimization results.

2. The method for optimizing the coordinated operation of a cross-border integrated energy system considering carbon-green certificate joint trading as described in claim 1, characterized in that: S11. The Cross-Border Integrated Energy System (CBIES) model constructed in step S1 specifically includes a gas turbine model, a gas boiler model, a renewable energy unit model, an energy storage device model, and an electric boiler model, as detailed below: The gas turbine model is as follows: in This represents the output electrical power of the gas turbine in the integrated energy system of the i-th country at time t. L represents the output thermal power of the gas turbine in the integrated energy system of the i-th country at time t. NG This indicates the lower heating value of natural gas. η represents the amount of natural gas consumed by the gas turbine at time t. GT N represents the power generation efficiency of a gas turbine. GT Indicates the waste heat recovery coefficient; The specific model of the gas-fired boiler is as follows: in L represents the output thermal power of the gas-fired boiler in the integrated energy system of the i-th country at time t. NG This indicates the lower heating value of natural gas. η represents the natural gas consumption of the gas-fired boiler at time t. GB This indicates the heating efficiency of the gas-fired boiler; Renewable energy unit models include wind power generation models and photovoltaic power generation models, as detailed below: Wind power generation model: in Let be the wind power output of the integrated energy system of the i-th country at time t, ρ be the air density, A be the swept area of ​​the wind turbine blades, v be the wind speed, and c be the wind power output of the i-th country. wt The wind energy utilization coefficient is the ratio of the wind energy absorbed by the wind turbine per unit time to the total wind energy passing through the rotor's rotating surface, λ. wt The ratio of leaf tip speed; Photovoltaic power generation model: T s =T a +0.0138·(1+0.031T a )·(1-0.042v)·G(6) where Let G be the photovoltaic power generation of the integrated energy system of the i-th country at time t, and let G be the solar irradiance. s P is the surface temperature of the photovoltaic cell. stc G stc T stc These represent the maximum output power, light intensity, and photovoltaic cell surface temperature under standard test conditions, respectively, where ε is the photovoltaic cell temperature coefficient; T a v represents ambient temperature; v represents wind speed. Energy storage device models include battery models and thermal storage tank models, as detailed below: Battery model: Discharge Charge in Let W be the amount of electricity stored in the battery of the i-th country's integrated energy system at time t+1. t e,i Let be the amount of electricity stored in the battery of the i-th country's integrated energy system at time t. and These represent the battery discharge power and charging power of the i-th country's integrated energy system at time t, respectively. and These are the battery's own discharge efficiency and charging efficiency, respectively. and These are the battery's own discharge loss and charging loss, respectively. Thermal storage tank model: Exothermic heating in W represents the thermal energy stored in the thermal storage tank of the i-th country's integrated energy system at time t+1. t h,i The thermal energy stored in the thermal storage tank of the i-th country's integrated energy system at time t. and These are the heat release power and heat charging power of the thermal storage tank at time t in the integrated energy system of the i-th country, respectively. and These are the heat release efficiency and heat charging efficiency of the heat storage tank itself, respectively. and These are the heat release loss and heat charging loss of the heat storage tank itself, respectively. The electric boiler model is as follows: in, Let be the heating capacity of the electric boiler in the integrated energy system of the i-th country at time t; Let η be the electrical power required by the electric boiler in the integrated energy system of the i-th country at time t. EB This refers to the conversion efficiency of the electric boiler.

3. The method for optimizing the coordinated operation of a cross-border integrated energy system considering carbon-green certificate joint trading as described in claim 2, characterized in that: The solver is either the CPLEX solver or the GUROBI solver.

Images