Dispatching and trading strategy of integrated power generation, grid, load and storage system considering electricity-carbon market synergy

Abstract

Considering the dispatching and trading strategy of the integrated source-grid-load-storage system with the coordination of the electricity-carbon market, a mathematical model of the integrated source-grid-load-storage system is established; a voluntary carbon emission reduction mechanism including day-ahead-day carbon offset, intraday carbon transfer and utilization, and intraday carbon emission reduction subsidies is established; a set of wind-solar uncertainty scenarios is generated based on non-parametric kernel density estimation and Wasserstein distance, and a multi-stage electricity-carbon market bidding strategy and benefit distribution model is established under the worst distribution of wind and solar in the wind-solar uncertainty scenario set; asymmetric Nash bargaining theory is introduced to solve the multi-stage electricity-carbon market bidding strategy and benefit distribution model; the Myerson value of each integrated source-grid-load-storage system is calculated according to the solution results, and the benefit distribution is carried out based on the Myerson value. The present invention takes the multi-subject integrated source-grid-load-storage system as the object, carries out the coordination of the electricity-carbon market and multi-stage optimal dispatching, so that the dispatching strategy can take into account both low-carbon and economic efficiency, and efficiently handle the regulation and dispatch between different energy flows.

Description

Technical Field

[0001] The present invention belongs to the technical field of electric power dispatching, and specifically relates to a dispatching and trading strategy for a source-grid-load-storage integrated system taking into account the synergy of the electricity-carbon market. Background Art

[0002] The emergence of an integrated source-grid-load-storage system with multi-energy complementary technology as its core has played a key supporting role in building a new type of power system that is clean, low-carbon, safe, controllable, flexible, efficient, intelligent, friendly, open and interactive. It can also improve the energy supply and demand system from various links such as production, conversion, transmission, storage and consumption, enhance the self-consistency of source-grid-load-storage, and help reduce carbon emissions.

[0003] However, with the increase in the penetration rate of intermittent renewable energy such as wind and light, the operational reliability and flexibility of the integrated source-grid-load-storage system will inevitably be affected. First, the randomness of renewable energy output leads to fluctuations in the energy supply of the integrated source-grid-load-storage system, making it difficult to meet load demand. Secondly, the increase in the penetration rate of renewable energy reduces the proportion of reliable energy supply units in the structure of the integrated source-grid-load-storage system, thereby weakening the energy scheduling flexibility of the integrated source-grid-load-storage system. It can be seen that eliminating the influence of uncertain factors is the key to ensuring the stability of the integrated source-grid-load-storage system and improving the utilization rate of renewable energy.

[0004] In addition, as my country's carbon market begins to operate, it is expected to play an important role in promoting the development of renewable energy power generation, the construction of new power systems, and the realization of carbon emission reduction targets. However, due to the lack of transparent and widely available market information, the current trading varieties and pricing methods of my country's carbon market are relatively simple, and there is a lack of a fair and reasonable cooperative profit distribution mechanism, which makes it difficult for market participants to make sophisticated trading decisions.

[0005] Therefore, in order to improve the self-consistency of the integrated power generation, grid, load and storage system and promote the realization of the goal of reducing carbon emissions, it is necessary to systematically study the economic characteristics and carbon emission characteristics of the integrated power generation, grid, load and storage system under uncertain scenarios, the interaction between the electricity and carbon markets, and the coordination strategies of multiple subjects participating in the electricity and carbon markets. Summary of the invention

[0006] The purpose of the present invention is to provide a scheduling and trading strategy for a source-grid-load-storage integrated system that takes into account both low-carbon and economic efficiency and the synergy of the electricity-carbon market in order to address the above-mentioned problems existing in the prior art.

[0007] To achieve the above objectives, the technical solution of the present invention is as follows:

[0008] Considering the coordination of electricity and carbon markets, a dispatching and trading strategy for a source-grid-load-storage integrated system is proposed, wherein the dispatching and trading strategy comprises the following steps:

[0009] S1. Establish a mathematical model of the source-grid-load-storage integrated system;

[0010] S2. Establish a voluntary carbon emission reduction mechanism including day-ahead and intraday carbon offset, intraday carbon transfer and utilization, and intraday carbon emission reduction subsidies;

[0011] S3, generate a set of wind and solar uncertainty scenarios based on non-parametric kernel density estimation and Wasserstein distance, and determine the worst wind and solar distribution;

[0012] S4. Multi-stage electricity-carbon market bidding strategy and benefit allocation model under the worst distribution of wind and solar power;

[0013] S5. Introduce asymmetric Nash bargaining theory to solve the bidding strategy and benefit distribution model of multi-stage electricity-carbon market;

[0014] S6. Calculate the Myerson value of each source-grid-load-storage integrated system according to the solution results, and distribute benefits based on the Myerson value.

[0015] The mathematical model of the source-grid-load-storage integrated system includes:

[0016] CHP unit model:

[0017] ; ;

[0018] ;

[0019] In the above formula, is the output power of the gas turbine GT at time t; is the natural gas consumption of gas turbine GT at time t; The lower heating value of the gas turbine GT; is the thermal power output of the gas turbine GT at time t; , They are the power generation efficiency and heat loss rate of the gas turbine GT respectively; is the total heat energy recovered and converted by the waste heat boiler WHB at time t; is the operating efficiency of the waste heat boiler WHB;

[0020] P2G device model:

[0021] ; ; ;

[0022] In the above formula, is the hydrogen output power of electrolyzer EL at time t; is the input power of electrolyzer EL at time t; is the operating efficiency of the electrolyzer EL; is the hydrogen output power of the P2G equipment used to meet the hydrogen energy supply at time t; is the hydrogen input power of the methane reactor MR used to produce natural gas at time t; is the amount of natural gas produced by the methane reactor MR at time t; is the operating efficiency of the methane reactor MR;

[0023] HFC hydrogen fuel cell model:

[0024] ; ;

[0025] In the above formula, is the electrical output power of the hydrogen fuel cell HFC at time t; is the thermal output power of the hydrogen fuel cell HFC at time t; is the hydrogen input power of the hydrogen fuel cell HFC at time t; is the electricity conversion efficiency of the hydrogen fuel cell HFC; is the thermal conversion efficiency of hydrogen fuel cell HFC;

[0026] Carbon capture model:

[0027] ; ;

[0028] ; ;

[0029] ;

[0030] ; ;

[0031] In the above formula, is the carbon emission of the CHP system at time t; is the carbon emission intensity coefficient of the CHP system; is the carbon emissions not captured by the storage equipment CCUS at time t; The carbon capture efficiency of CCUS for storage equipment; The total carbon emissions of the integrated source-grid-load-storage system; is the carbon emission consumed by the P2G equipment to produce methane at time t; is the carbon emissions stored by the storage equipment CCUS at time t; is the carbon emissions captured by the storage device CCUS at time t; Carbon emissions transferred from the integrated source-grid-load-storage system; is the carbon consumption coefficient of methane produced by the methane reactor MR; , , They are the total energy consumption of the storage equipment CCUS at time t, the energy consumption generated by carbon capture, and the energy consumption generated by carbon storage; is the total energy consumption of the storage equipment CCUS at time t+1; , They are the unit carbon capture power consumption coefficient and carbon storage power consumption coefficient of the storage equipment CCUS respectively; , They are the upper and lower limits of CCUS operating power for storage equipment respectively; , They are the upper and lower limits of CCUS ramp rate for storage equipment, respectively;

[0032] The SOC state of charge of the energy storage device during charging:

[0033] ;

[0034] In the above formula, is the charge state of the energy storage device o at time t, and the energy storage device refers to the energy storage system ESS that is composed of the energy storage system; is the charge state of the energy storage device o at time t-1; Charging power of each energy storage device o; is the charging efficiency of the energy storage device o; is the capacity of the energy storage device o;

[0035] The SOC state of charge of the energy storage device when discharging energy:

[0036] ;

[0037] In the above formula, is the energy release power of energy storage device o; is the energy release efficiency of the energy storage device o.

[0038] The voluntary carbon emission reduction mechanism includes:

[0039] ;

[0040] ;

[0041] ; ; ;

[0042] ; ; ;

[0043] ;

[0044] In the above formula, is the day-ahead pre-bid electricity of the integrated power-grid-load-storage system i at time t; , They are the predicted power consumption of the electrolyzer EL and the storage equipment CCUS in the source-grid-load-storage integrated system i at time t; , They are respectively the day-ahead predicted charge and discharge capacity of the energy storage system i in the source-grid-load-storage integrated system at time t; is the power consumption of the electric load in the source-grid-load-storage integrated system i at time t; , , , They are the day-ahead predicted output values ​​of photovoltaic, wind power, hydrogen fuel cell and gas turbine in the source-grid-load-storage integrated system i at time t; is the electrical output power of the CHP unit; is the carbon emission of the CHP system in the integrated source-grid-load-storage system i; , They are the carbon quota benchmark values ​​for power supply and heating supply of CHP system respectively; is the energy supply correction factor of the CHP system; is the heat energy generated by the CHP system in the source-grid-load-storage integrated system i at time t; The carbon quota obtained by the integrated power generation, grid-load and storage system i in the pre-bid stage; Carbon quota obtained by using green certificates for the integrated source-grid-load-storage system; Indicates the amount of carbon quota that can be bid for the integrated source-grid-load-storage system i days ago; The surplus carbon quota of the integrated power grid-load-storage system i; represents the carbon transfer amount of the integrated source-grid-load-storage system i through intraday scheduling; , , , They are the total carbon emissions after the source-grid-load-storage integrated system i participates in carbon transfer, the carbon emissions of the CHP unit, the carbon quotas exchanged for green certificates, and the carbon quotas sold; , , They are the green certificates obtained by the source-grid-load-storage integrated system i through wind power, photovoltaic and green certificate trading markets; Represents the green certificate trading income of the source-grid-load-storage integrated system i; Represents the unit green certificate transaction price of the source-grid-load-storage integrated system i; is the carbon transfer cost of the integrated source-grid-load-storage system i; is the carbon transfer value of the integrated source-grid-load-storage system i and the integrated source-grid-load-storage system j at time t; is the carbon transfer cost coefficient; , are subsidy income and subsidy coefficient respectively.

[0045] The step S3 is specifically as follows:

[0046] S31. Establish the expression of total prediction error of wind power and photovoltaic power:

[0047] ; ; ;

[0048] In the above formula, represents the total prediction error of wind power and photovoltaic power at time t in the source-grid-load-storage integrated system i; , are the prediction errors of wind power and photovoltaic power at time t in the source-grid-load-storage integrated system i; , They are the actual outputs of wind power and photovoltaic power at time t in the source-grid-load-storage integrated system i; , They are the predicted power of wind power and photovoltaic power at time t in the source-grid-load-storage integrated system i;

[0049] S32. Construct the probability density function of the wind and solar random variable based on non-parametric kernel density estimation:

[0050] ;

[0051] In the above formula, for The probability density function of Indicates that the scenery changes randomly quantity ; Indicates the number of typical daily data samples of wind power and photovoltaic power; is bandwidth; is the kernel function;

[0052] S33. Construct a cumulative distribution function based on the probability density function:

[0053] ;

[0054] In the above formula, is the cumulative distribution function, also representing the empirical distribution;

[0055] S34. Determine the uncertainty set of the cumulative distribution function through the quantile of the wind and solar random variable:

[0056] ;

[0057] ;

[0058] In the above formula, , Respectively The upper and lower boundaries are determined by the quantile of the wind and solar random variable prediction error distribution, and the confidence level is set to ; is the uncertainty set;

[0059] S35, As an empirical distribution, the Wasserstein distance is used to describe the true distribution F and the empirical distribution Distance between:

[0060] ; ;

[0061] In the above formula, Represents the true distribution F and the empirical distribution Wasserstein distance between them; is the supremum function; is the intermediate calculation amount; Represents the empirical distribution value of wind power and photovoltaic forecast data;

[0062] S36, constructing a distributed fuzzy uncertainty set based on Wasserstein distance, wherein the distributed fuzzy uncertainty set can be regarded as an empirical distribution The Wasserstein sphere with the center as , and the true distribution F as the radius, is expressed as:

[0063] ;

[0064] In the above formula, Represents a distributed fuzzy uncertainty set; represents the set of true distributions F; represents the radius of the Wasserstein sphere, satisfying .

[0065] The multi-stage electricity-carbon market bidding strategy and benefit allocation model includes a first-stage objective function and its constraints, a second-stage objective function and its constraints; the first-stage objective function includes:

[0066] ;

[0067] ;

[0068] ; ; ;

[0069] ; ;

[0070] In the above formula, is the objective function of one stage; represents the day-ahead period return; is the objective function of the second stage; is the day-ahead electricity bidding price of the integrated power grid-load-storage system i; The estimated winning bid power of the source-grid-load-storage integrated system i; is the day-ahead carbon bidding price of the integrated power grid-load-storage system i; The expected winning carbon quota for the integrated source-grid-load-storage system i; is the pre-dispatch cost of the source-grid-load-storage integrated system i; , , , They are the operation and maintenance cost, gas purchase cost, carbon capture cost, and demand response cost of the source-grid-load-storage integrated system i; The amount of natural gas purchased for the integrated power grid, load and storage system i; is the price of natural gas; The green certificate trading volume of the integrated source-grid-load-storage system i; is the green certificate price; is the total cost of carbon capture and storage for the integrated source-grid-load-storage system i at time t; is the cost of carbon capture for the integrated source-grid-load-storage system i at time t; is the cost of carbon storage for the integrated source-grid-load-storage system i at time t; is the interruption load factor; is the elastic load coefficient; is the interruptible power of the integrated source-grid-load-storage system i at time t Load capacity; is the elastically transferable energy of the integrated source-grid-load-storage system i at time t Load capacity; It means electricity / heat / hydrogen.

[0071] The constraints of the first stage include:

[0072] Exchange power constraints for energy trading:

[0073] ;

[0074] In the above formula, is the day-ahead pre-bid electricity of the integrated power-grid-load-storage system i at time t; , They are the upper and lower limits of the power of the source-grid-load-storage integrated system i respectively; is the gas consumption of the integrated source-grid-load-storage system i at time t; , are the upper and lower limits of gas consumption of the source-grid-load-storage integrated system i respectively;

[0075] Energy supply and demand balance constraints:

[0076] ;

[0077] In the above formula, , , , They are the day-ahead predicted outputs of photovoltaic power, wind power, storage equipment CCUS, and electrolyzer EL in the source-grid-load-storage integrated system i at time t; is the electrical output power of the hydrogen fuel cell HFC at time t; , are the heat outputs of the waste heat boiler WHB and the hydrogen fuel cell HFC at time t respectively; , are the heat storage and release of the thermal energy storage system at time t respectively; , are the storage and discharge amounts of the electric energy storage system at time t respectively; , are the storage and release amounts of hydrogen in the hydrogen energy storage system at time t, respectively; , , They are heat load, electricity load and hydrogen load at time t respectively; , , are the hydrogen input powers of the electrolyzer EL, hydrogen fuel cell HFC, and methane reactor MR at time t, respectively; is the amount of natural gas produced by the methane reactor MR at time t; is the natural gas consumption of gas turbine GT at time t;

[0078] Carbon emission flow balance constraint:

[0079] ;

[0080] In the above formula, is the carbon emission of the CHP unit at time t; is the carbon emission consumed by the P2G equipment to produce methane after carbon capture at time t; is the carbon emissions not captured by the storage equipment CCUS at time t; is the carbon emissions captured by the storage equipment CCUS at time t; is the carbon emissions consumed by the storage equipment CCUS for carbon storage at time t;

[0081] Demand response constraints:

[0082] ;

[0083] ; ;

[0084] In the above formula, is the time t load, They represent electricity / heat / hydrogen respectively; , They represent the elastic transferability at time t. Load, interruptible load; , are elastically transferable at time t Load, interruptible The upper limit of the load; is the fixed value at time t load.

[0085] The objective function of the two stages includes:

[0086] ;

[0087] ;

[0088] ; ; ;

[0089] ;

[0090] ;

[0091] ;

[0092] ; ;

[0093] ;

[0094] In the above formula, is the objective function of the second stage, expressed in Wasserstein fuzzy sets The worst distribution of scenery The intraday scheduling cost under ; is the intraday re-bidding cost of the integrated source-grid-load-storage system i at time t; , , , They are respectively the operation and maintenance cost, gas purchase cost, carbon capture cost, and demand response cost of the source-grid-load-storage integrated system i after adjusting the operating status of the unit; is the adjusted power of device k in the source-grid-load-storage integrated system i at time t; Adjust the power of gas purchase for the integrated source-grid-load-storage system i; , They are the interruptible electric / heat / hydrogen load adjustment power and the elastically transferable electric / heat / hydrogen load adjustment power of the source-grid-load-storage integrated system i at time t, respectively; The cost of the integrated source-grid-load-storage system i consuming green certificates in exchange for carbon quotas at time t; is the penalty cost of the source-grid-load-storage integrated system i at time t; is the voluntary carbon emission reduction subsidy for the integrated power grid-load-storage system i at time t; The carbon quota obtained by the integrated source-grid-load-storage system i in the pre-bidding stage before time t; is the actual carbon emission of the integrated source-grid-load-storage system i during the intraday period at time t; is the subsidy coefficient; The cost incurred by the integrated source-grid-load-storage system i when transferring carbon emissions to other systems at time t; is the energy sharing income of the integrated source-grid-load-storage system i at time t; is the carbon offset adjustment of the integrated source-grid-load-storage system i at time t; is the voluntary carbon emission reduction benefit of the integrated source-grid-load-storage system i at time t; is the operation and maintenance cost coefficient of equipment k; is the price of natural gas; is the carbon capture adjustment amount of the integrated source-grid-load-storage system i at time t; is the carbon storage adjustment of the integrated source-grid-load-storage system i at time t; is the interruption load factor; is the elastic load coefficient; The green certificate quota adjustment amount for the integrated power grid-load-storage system i; is the carbon offset of the integrated source-grid-load-storage system i at time t; The actual adjustment amount of the green certificate of the integrated power grid-load-storage system i; is the green certificate price; is the wind abandonment factor; is the abandoned wind volume; is the light abandonment coefficient; is the amount of abandoned light; is the load shedding factor; is the load adjustment amount; is the carbon transfer cost coefficient; is the carbon transfer from source-grid-load-storage integrated system i to source-grid-load-storage integrated system j; The total number of integrated source-grid-load-storage systems that receive energy.

[0095] The constraints of the second stage include:

[0096] Energy supply and demand balance compliance re-dispatch constraints:

[0097] ;

[0098] In the above formula, , , , , They are the adjusted outputs of photovoltaic power, wind power, storage equipment CCUS, electrolyzer EL, and gas turbine GT in the integrated source-grid-load-storage system i at time t during the intraday phase; is the adjusted power purchase amount of system i during the intraday period; is the adjusted electric power of the hydrogen fuel cell HFC at time t during the intra-day phase; , They are the adjusted thermal power of the waste heat boiler WHB and the hydrogen fuel cell HFC at time t in the intra-day stage; , They are respectively the heat storage and heat release adjustment of the thermal energy storage system at time t in the intra-day stage; , They are the storage and discharge adjustment of the electric energy storage system at time t during the intra-day phase; , They are the storage and release adjustment amounts of hydrogen in the hydrogen energy storage system at time t during the intra-day stage; , , They are the heat load, electricity load and hydrogen load adjustment at time t in the day phase respectively; , They are the adjusted hydrogen power of the hydrogen fuel cell HFC and the methane reactor MR at time t in the intra-day stage; is the adjustment amount of natural gas produced by the methane reactor MR at time t during the intra-day period; is the adjusted natural gas consumption of the gas turbine GT at time t during the intra-day phase; , , They represent the amount of electric energy, thermal energy and hydrogen used to prepare methane transferred from the integrated source-grid-load-storage system i to the integrated source-grid-load-storage system j during the day; It represents the adjustment amount of natural gas purchased by the integrated source-grid-load-storage system i during the intra-day period; represents the amount of carbon dioxide for methane production transferred from the integrated source-grid-load-storage system i to the integrated source-grid-load-storage system j; It represents the hydrogen adjustment amount of the integrated source-grid-load-storage system i in the intra-day stage;

[0099] Energy sharing balance constraints for compliance redispatch:

[0100] ;

[0101] In the above formula, , , , It represents the amount of electricity, heat, hydrogen energy and natural gas output from the integrated source-grid-load-storage system i to the integrated source-grid-load-storage system j; , , , It represents the amount of electricity, heat, hydrogen energy and natural gas received by the integrated source-grid-load-storage system j; The total number of integrated source-grid-load-storage systems that output energy; Indicates the total number of integrated source-grid-load-storage systems receiving energy;

[0102] Energy sharing price constraints for compliance and redispatch:

[0103] ;

[0104] In the above formula, Indicates the lower limit of the selling price of electricity / heat / hydrogen energy; Indicates the shared price of electricity / heat / hydrogen energy; Indicates the lower limit of the purchase price of electricity / heat / hydrogen energy;

[0105] Carbon emission flow balance constraints for compliance redispatching:

[0106] ;

[0107] In the above formula, It represents the carbon emission of CHP unit in the intra-day period; It indicates that CCUS equipment adjusts the amount of carbon captured during the intra-day period; It represents the adjusted carbon emissions consumed by the P2G equipment to produce methane after carbon capture during the day; Indicates the amount of carbon stored by CCUS equipment during the intra-day period; represents the amount of carbon dioxide transferred from the integrated source-grid-load-storage system i to the integrated source-grid-load-storage system j; Indicates the carbon emission adjustment amount during the intraday period.

[0108] The S5 is specifically:

[0109] S51. Introduce asymmetric Nash bargaining theory to describe the cooperative game within the source-grid-load-storage integrated system alliance, and convert the multi-stage electricity-carbon market bidding strategy and benefit distribution model into the following asymmetric Nash bargaining model:

[0110] ;

[0111] In the above formula, represents the benefits obtained by the integrated source-grid-load-storage system i participating in Nash bargaining; It represents the benefits of the integrated power grid-load-storage system when it operates independently, i.e. the negotiation breakdown point; It represents the energy sharing contribution index of the source-grid-load-storage integrated system i; Indicates the total number of source-grid-load-storage integrated systems;

[0112] S52. Convert the asymmetric Nash bargaining model into the alliance benefit maximization problem, the energy sharing contribution quantification problem, and the cooperation benefit allocation problem, and solve these three sub-problems in turn; the alliance benefit maximization problem includes:

[0113] ;

[0114] In the above formula, is the objective function of the second stage; The cost of transferring carbon emissions from the integrated power grid, load and storage system i to other integrated power grid, load and storage systems; Energy sharing benefits of the integrated source-grid-load-storage system i; It indicated that it would not participate in energy sharing and carbon transfer;

[0115] The energy sharing contribution quantification issues include:

[0116] ;

[0117] ;

[0118] In the above formula, , They are the energy sharing adjustment coefficient and the carbon transfer adjustment coefficient; The power sold / maximum power sold of the integrated power-grid-load-storage system i; The purchased power / maximum purchased power of the integrated power-grid-load-storage system i; is the carbon output / maximum carbon output of the integrated source-grid-load-storage system i; is the carbon input / maximum carbon input of the source-grid-load-storage integrated system i; The energy traded from source-grid-load-storage integrated system i to source-grid-load-storage integrated system j; The energy traded from the integrated source-grid-load-storage system j to the integrated source-grid-load-storage system i; is the amount of carbon transferred from the integrated source-grid-load-storage system i to the integrated source-grid-load-storage system j; is the amount of carbon transferred from the integrated source-grid-load-storage system j to the integrated source-grid-load-storage system i;

[0119] The cooperation benefit distribution issues include:

[0120] ;

[0121] In the above formula, , They represent the energy sharing benefits and carbon transfer costs when the optimal solution is obtained for the alliance benefit maximization problem; is the redispatching cost.

[0122] The calculation formula of the Myerson value is:

[0123] ;

[0124] ;

[0125] In the above formula, is the Myerson value of the source-grid-load-storage integrated system i; Indicates transaction status; alliance It is a collection of source-grid-load-storage integrated system i. yes A subset of is the total number of samples; For the Alliance The number of Zhongyuan grid-load-storage integrated systems; Indicates the marginal contribution in cooperative games; It represents the marginal contribution in non-cooperative game, which is directly calculated without considering carbon transfer and energy sharing strategy transactions; Indicates that the source-grid-load-storage integrated system forms a cooperative game alliance Total revenue; Represents the Alliance for Independent Operation of Source-Grid-Load-Storage Integrated System total revenue.

[0126] Compared with the prior art, the present invention has the following beneficial effects:

[0127] The scheduling and trading strategy of the source-grid-load-storage integrated system considering the coordination of the electricity-carbon market described in the present invention establishes a mathematical model of the source-grid-load-storage integrated system; establishes a voluntary carbon emission reduction mechanism including day-ahead-intraday carbon offset, intraday carbon transfer and utilization, and intraday carbon emission reduction subsidies; generates a set of wind-solar uncertainty scenarios based on non-parametric kernel density estimation and Wasserstein distance, and determines the worst distribution of wind-solar; establishes a multi-stage electricity-carbon market bidding strategy and benefit distribution model under the worst distribution of wind-solar; introduces asymmetric Nash bargaining theory to solve the multi-stage electricity-carbon market bidding strategy and benefit distribution model; calculates the Myerson value of each source-grid-load-storage integrated system according to the solution results, and distributes benefits based on the Myerson value. The present invention takes the multi-subject source-grid-load-storage integrated system as the object, carries out electricity-carbon market coordination and multi-stage optimization scheduling, so that the scheduling strategy can take into account both low-carbon and economic efficiency, and efficiently handle the regulation and scheduling between different energy flows. BRIEF DESCRIPTION OF THE DRAWINGS

[0128] Figure 1 It is a flow chart of the present invention. DETAILED DESCRIPTION

[0129] The present invention is further described in detail below in conjunction with specific implementations and drawings.

[0130] See also Figure 1 , a dispatching and trading strategy for the integrated source-grid-load-storage system considering the synergy of the electricity-carbon market is carried out in the following steps:

[0131] S1. Establish mathematical model of source-grid-load-storage integrated system

[0132] The source-grid-load-storage integrated system in this embodiment includes a combined heat and power unit CHP, P2G equipment, a hydrogen fuel cell HFC, a carbon capture and storage equipment CCUS, an energy storage facility ESS, etc. The CHP unit produces electricity through a gas turbine GT, and the waste heat generated when the gas turbine generates electricity is recovered through a waste heat boiler WHB, thereby realizing the combined supply of heat and power; the P2G equipment can realize the conversion of electricity and hydrogen energy, including a methane reactor MR and an electrolyzer EL. The electrolyzer is used to convert electricity into hydrogen. Part of the generated hydrogen is used to supply hydrogen energy, and the other part is used to react with CO in the methane reactor. 2 The reaction produces methane and supplies natural gas; hydrogen energy can be converted into electricity and heat energy through hydrogen fuel cells HFC; the carbon capture and storage equipment CCUS adopts a post-combustion capture system to capture CO emitted by the CHP cogeneration unit 2 ; The energy storage facility ESS considers three types of energy storage devices: electricity, heat, and hydrogen, and introduces the SOC state of charge variable to describe the state of each energy storage device;

[0133] The mathematical model of the source-grid-load-storage integrated system constructed includes:

[0134] CHP unit model:

[0135] ; ;

[0136] ;

[0137] In the above formula, is the output power of the gas turbine GT at time t; is the natural gas consumption of gas turbine GT at time t; The lower heating value of the gas turbine GT; is the thermal power output of the gas turbine GT at time t; , They are the power generation efficiency and heat loss rate of the gas turbine GT respectively; is the total heat energy recovered and converted by the waste heat boiler WHB at time t; is the operating efficiency of the waste heat boiler WHB;

[0138] P2G device model:

[0139] ; ; ;

[0140] In the above formula, is the hydrogen output power of electrolyzer EL at time t; is the input power of electrolyzer EL at time t; is the operating efficiency of the electrolyzer EL; is the hydrogen output power of the P2G equipment used to meet the hydrogen energy supply at time t; is the hydrogen input power of the methane reactor MR used to produce natural gas at time t; is the amount of natural gas produced by the methane reactor MR at time t; is the operating efficiency of the methane reactor MR;

[0141] HFC hydrogen fuel cell model:

[0142] ; ;

[0143] In the above formula, is the electrical output power of the hydrogen fuel cell HFC at time t; is the thermal output power of the hydrogen fuel cell HFC at time t; is the hydrogen input power of the hydrogen fuel cell HFC at time t; is the electricity conversion efficiency of the hydrogen fuel cell HFC; is the thermal conversion efficiency of hydrogen fuel cell HFC;

[0144] Carbon capture model:

[0145] ; ;

[0146] ; ;

[0147] ;

[0148] ; ;

[0149] In the above formula, is the carbon emission of the CHP system at time t; is the carbon emission intensity coefficient of the CHP system; is the carbon emissions not captured by the storage equipment CCUS at time t; The carbon capture efficiency of CCUS for storage equipment; The total carbon emissions of the integrated source-grid-load-storage system; is the carbon emission consumed by the P2G equipment to produce methane at time t; is the carbon emissions stored by the storage equipment CCUS at time t; is the carbon emissions captured by the storage device CCUS at time t; Carbon emissions transferred from the integrated source-grid-load-storage system; is the carbon consumption coefficient of methane produced by the methane reactor MR; , , They are the total energy consumption of the storage equipment CCUS at time t, the energy consumption generated by carbon capture, and the energy consumption generated by carbon storage; is the total energy consumption of the storage equipment CCUS at time t+1; , They are the unit carbon capture power consumption coefficient and carbon storage power consumption coefficient of the storage equipment CCUS respectively; , They are the upper and lower limits of CCUS operating power for storage equipment respectively; , They are the upper and lower limits of CCUS ramp rate for storage equipment, respectively;

[0150] The SOC state of charge of the energy storage device during charging:

[0151] ;

[0152] In the above formula, is the charge state of the energy storage device o at time t. The energy storage device refers to energy storage devices such as batteries, hydrogen storage tanks, and heat storage tanks, which together constitute the energy storage system ESS; is the state of charge of the energy storage device o at time t-1; is the charging power of the energy storage device o; is the charging efficiency of the energy storage device o; is the capacity of the energy storage device o;

[0153] The SOC state of charge of the energy storage device when discharging energy:

[0154] ;

[0155] In the above formula, is the energy release power of the energy storage device o; is the energy release efficiency of the energy storage device o;

[0156] S2. Establish a voluntary carbon emission reduction mechanism including day-ahead and intraday carbon offset, intraday carbon transfer and utilization, and intraday carbon emission reduction subsidies;

[0157] In the day-ahead pre-bid stage, the integrated source-grid-load-storage system gives priority to ensuring the stable supply of loads, and optimizes its own system energy supply status based on the predicted values ​​of wind and solar random variables and the day-ahead power trading price. The integrated source-grid-load-storage system calculates its own surplus electricity and independently submits the pre-bid electricity for each time period to the dispatch center to generate pre-bid revenue; the calculation formula for the day-ahead pre-bid electricity is:

[0158] ;

[0159] In the above formula, is the day-ahead pre-bid electricity of the integrated power-grid-load-storage system i at time t; , They are the predicted power consumption of the electrolyzer EL and the storage equipment CCUS in the source-grid-load-storage integrated system i at time t; , They are respectively the day-ahead predicted charge and discharge capacity of the energy storage system i in the source-grid-load-storage integrated system at time t; is the power consumption of the electric load in the source-grid-load-storage integrated system i at time t; , , , They are the day-ahead predicted output values ​​of photovoltaic, wind power, hydrogen fuel cell and gas turbine in the source-grid-load-storage integrated system i at time t;

[0160] Before the integrated source-grid-load-storage system participates in the carbon quota bidding, the regulatory center allocates the initial carbon emission quota; the allocation method can be described as follows:

[0161] ;

[0162] In the above formula, in the above formula, is the electrical output power of the CHP unit; is the carbon emission of the CHP system in the integrated source-grid-load-storage system i; , They are the carbon quota benchmark values ​​for power supply and heating supply of CHP system respectively; is the energy supply correction factor of the CHP system; is the heat energy generated by the CHP system in the source-grid-load-storage integrated system i at time t;

[0163] In the day-ahead stage, the integrated power grid, load and storage system adjusts the unit output on the premise of giving priority to its own day-ahead power supply, and uses CCUS for carbon capture to reduce system carbon emissions and increase system carbon quotas; surplus electricity and carbon quotas can participate in bidding in the electricity market and carbon market; each integrated power grid, load and storage system bids independently and calculates its own expected bid electricity and expected bid carbon quota; at the same time, the integrated power grid, load and storage system can exchange green certificates for more carbon quotas; the carbon quotas obtained in the pre-bidding stage can be expressed as follows:

[0164] ;

[0165] In the above formula, The carbon quota obtained by the integrated power generation, grid-load and storage system i in the pre-bid stage; Carbon quota obtained by using green certificates for the integrated source-grid-load-storage system;

[0166] The total carbon quota obtained in the pre-bidding stage shall be used for its own carbon emission consumption first, and the surplus carbon quota can be used for bidding in the day-ahead carbon trading center; the day-ahead carbon quota that can be bid for the integrated source-grid-load-storage system can be expressed as:

[0167] ;

[0168] In the above formula, Indicates the amount of carbon quota that can be bid for the integrated source-grid-load-storage system i days ago; The surplus carbon quota of the integrated power grid-load-storage system i;

[0169] In the intraday stage, multiple integrated source-grid-load-storage systems form an alliance. The alliance optimizes the system output based on energy sharing, carbon transfer, etc. to reduce the actual carbon emissions of the system and obtain more carbon quotas through green certificate-carbon offset. The biddable carbon quota after the intraday stage of cooperative game can be expressed as:

[0170] ;

[0171] In the above formula, represents the carbon transfer amount of the integrated source-grid-load-storage system i through intraday scheduling; , , , They are the total carbon emissions of the integrated source-grid-load-storage system i after participating in carbon transfer in the daily stage, the carbon emissions of the CHP unit, the carbon quotas exchanged for green certificates, and the carbon quotas sold;

[0172] In order to increase the biddable carbon quota of the integrated source-grid-load-storage system, a green certificate-carbon offset mechanism is proposed, that is, the integrated source-grid-load-storage system can consume part of the green certificates in exchange for carbon quotas when calculating the expected carbon emissions in the day-ahead pre-bidding stage and the actual carbon emissions in the intraday real-time bidding stage; green certificates can be obtained through clean energy power generation certification or purchased from the green certificate trading market; the number of green certificates obtained by the integrated source-grid-load-storage system can be expressed as:

[0173] ;

[0174] In the above formula, Indicates the number of green certificates obtained by the source-grid-load-storage integrated system; , , They are the green certificates obtained by the source-grid-load-storage integrated system i through wind power, photovoltaic and green certificate trading markets;

[0175] In the day-ahead stage, the source-grid-load-storage integrated system is certified for green certificates through the predicted value of wind-solar random variables, and one green certificate is exchanged for each MW of clean energy power generation; in the intraday stage, green certificate certification is carried out considering the wind-solar output error penalty, referring to the average historical transaction price in the "China Green Power Certificate Subscription and Trading Platform", assuming that the green certificate transaction prices of wind power and photovoltaic power are equal, both are 130 yuan / book; the green certificate transaction income of the source-grid-load-storage integrated system can be expressed as:

[0176] ; ;

[0177] In the above formula, Indicates the amount of green certificates obtained by the integrated source-grid-load-storage system i through wind power / photovoltaic; represents the predicted value of wind power / photovoltaic output in the source-grid-load-storage integrated system i; Represents the green certificate trading income of the source-grid-load-storage integrated system i; Represents the unit green certificate transaction price of the source-grid-load-storage integrated system i;

[0178] In the intraday stage, the alliance composed of multiple integrated power grid, load and storage systems also optimizes scheduling through cooperative games. That is, the high-carbon emission integrated power grid, load and storage system will use the carbon dioxide captured by CCUS to transfer and transport it to the low-carbon emission integrated power grid, load and storage system for use. The carbon transferor bears the transfer cost and the carbon transferor obtains the transfer income. Through cooperative games and carbon transfer, the carbon emissions of the integrated power grid, load and storage system will be significantly reduced, and the regulatory agency will provide voluntary carbon emission reduction subsidies to the integrated power grid, load and storage system. This process can be expressed as:

[0179] ; ;

[0180] In the above formula, is the carbon transfer cost of the integrated source-grid-load-storage system i; is the carbon transfer value of the integrated source-grid-load-storage system i and the integrated source-grid-load-storage system j at time t; is the carbon transfer cost coefficient; , are subsidy income and subsidy coefficient respectively;

[0181] S3. Generate scenery uncertainty scene set based on nonparametric kernel density estimation and Wasserstein distance

[0182] S31. Since the electricity-carbon market bidding strategy of the integrated source-grid-load-storage system will be affected by the uncertainty of wind power and photovoltaic power, the prediction error of wind power and photovoltaic power can be expressed as follows:

[0183] ; ; ;

[0184] In the above formula, represents the total prediction error of wind power and photovoltaic power at time t in the source-grid-load-storage integrated system i; , are the prediction errors of wind power and photovoltaic power at time t in the source-grid-load-storage integrated system i; , They are the actual outputs of wind power and photovoltaic power at time t in the source-grid-load-storage integrated system i; , They are the predicted power of wind power and photovoltaic power at time t in the source-grid-load-storage integrated system i;

[0185] S32. Construct the probability density function of the wind and solar random variable based on non-parametric kernel density estimation:

[0186] ;

[0187] In the above formula, for The probability density function of Indicates that the scenery changes randomly quantity ; Indicates the number of typical daily data samples of wind power and photovoltaic power; is bandwidth; is the kernel function;

[0188] S33. Construct a cumulative distribution function based on the probability density function:

[0189] ;

[0190] In the above formula, is the cumulative distribution function, also representing the empirical distribution;

[0191] S34. The uncertainty set of the cumulative distribution function is determined by the quantile of the wind and solar random variable, which reduces the robustness of the uncertainty set without losing real information:

[0192] ; ;

[0193] In the above formula, , Respectively The upper and lower boundaries are determined by the quantile of the wind and solar random variable prediction error distribution, and the confidence level is set to ; is the uncertainty set;

[0194] S35, As an empirical distribution, the Wasserstein distance is used to describe the empirical distribution The distance from the true distribution F is:

[0195] ; ;

[0196] In the above formula, Represents the true distribution F and the empirical distribution Wasserstein distance between them; is the supremum function; is the intermediate calculation amount; Represents the empirical distribution value of wind power and photovoltaic forecast data;

[0197] S36, constructing a distributed fuzzy uncertainty set based on Wasserstein distance, wherein the distributed fuzzy uncertainty set can be regarded as an empirical distribution The Wasserstein sphere with the center as , and the true distribution F as the radius, is expressed as:

[0198] ;

[0199] In the above formula, Represents a distributed fuzzy uncertainty set; represents the set of true distributions F; represents the radius of the Wasserstein sphere, satisfying ;

[0200] S4. Multi-stage electricity-carbon market bidding strategy and benefit allocation model under the worst distribution of wind and solar power

[0201] The multi-stage electricity-carbon market bidding strategy and benefit distribution model includes the objective function and constraints of the first stage and the objective function and constraints of the second stage. The objective function of the first stage is to maximize the sum of the day-ahead phase revenue and the intraday phase scheduling cost under the worst distribution of wind and solar forecast errors, which is specifically expressed as follows:

[0202] ;

[0203] ; ; ;

[0204] ; ; ;

[0205] In the above formula, is the objective function of one stage; represents the day-ahead period return; is the objective function of the second stage, i.e., minimizing the scheduling cost during the intraday period under the worst distribution of wind and solar forecast errors; is the day-ahead electricity bidding price of the integrated power grid-load-storage system i; The estimated winning bid power of the source-grid-load-storage integrated system i; is the day-ahead carbon bidding price of the integrated power grid-load-storage system i; The expected winning carbon quota for the integrated source-grid-load-storage system i; is the pre-dispatch cost of the source-grid-load-storage integrated system i; , , , They are the operation and maintenance cost, gas purchase cost, carbon capture cost, and demand response cost of the source-grid-load-storage integrated system i; The amount of natural gas purchased for the integrated power grid, load and storage system i; is the price of natural gas; The green certificate trading volume of the integrated source-grid-load-storage system i; is the green certificate price; is the total cost of carbon capture and storage for the integrated source-grid-load-storage system i at time t; is the cost of carbon capture for the integrated source-grid-load-storage system i at time t; is the cost of carbon storage for the integrated source-grid-load-storage system i at time t; is the interruption load factor; is the elastic load coefficient; is the interruptible power of the integrated source-grid-load-storage system i at time t Load capacity; is the elastically transferable energy of the integrated source-grid-load-storage system i at time t Load capacity; It means electricity / heat / hydrogen.

[0206] The constraints of the first stage include:

[0207] Exchange power constraints for energy trading:

[0208] ;

[0209] In the above formula, is the day-ahead pre-bid electricity of the integrated power-grid-load-storage system i at time t; , They are the upper and lower limits of the power of the source-grid-load-storage integrated system i respectively; is the gas consumption of the integrated source-grid-load-storage system i at time t; , are the upper and lower limits of gas consumption of the source-grid-load-storage integrated system i respectively;

[0210] Energy supply and demand balance constraints:

[0211] ;

[0212] In the above formula, , , , They are the day-ahead predicted outputs of photovoltaic power, wind power, storage equipment CCUS, and electrolyzer EL in the source-grid-load-storage integrated system i at time t; is the electrical output power of the hydrogen fuel cell HFC at time t; , are the heat outputs of the waste heat boiler WHB and the hydrogen fuel cell HFC at time t respectively; , are the heat storage and release of the thermal energy storage system at time t respectively; , are the storage and discharge amounts of the electric energy storage system at time t respectively; , are the storage and release amounts of hydrogen in the hydrogen energy storage system at time t, respectively; , , They are heat load, electricity load and hydrogen load at time t respectively; , , are the hydrogen input powers of the electrolyzer EL, hydrogen fuel cell HFC, and methane reactor MR at time t, respectively; is the amount of natural gas produced by the methane reactor MR at time t; is the natural gas consumption of gas turbine GT at time t;

[0213] Carbon emission flow balance constraint:

[0214] ;

[0215] In the above formula, is the carbon emission of the CHP unit at time t; is the carbon emission consumed by the P2G equipment to produce methane after carbon capture at time t; is the carbon emissions not captured by the storage equipment CCUS at time t; is the carbon emissions captured by the storage equipment CCUS at time t; is the carbon emissions consumed by the storage equipment CCUS for carbon storage at time t;

[0216] Demand response constraints:

[0217] ;

[0218] ; ;

[0219] In the above formula, is the time t load, They represent electricity / heat / hydrogen respectively; , They represent the elastic transferability at time t. Load, interruptible load; , are elastically transferable at time t Load, interruptible The upper limit of the load; is the fixed value at time t load.

[0220] In the intraday stage, the total forecast error of wind power and photovoltaic power is taken into account, and the pre-dispatch plan of the day-ahead stage is adjusted by means of load reduction, energy sharing, voluntary carbon emission reduction, etc. to effectively deal with the forecast error; under the worst distribution of wind and photovoltaic forecast errors, a real-time dispatch plan is formulated with the goal of minimizing the dispatch cost of the intraday stage; the objective functions of the two stages include:

[0221] ;

[0222] ; ;

[0223] ; ;

[0224] ; ;

[0225] ;

[0226] ; ;

[0227] ;

[0228] In the above formula, is the objective function of the second stage, expressed in Wasserstein fuzzy sets The worst distribution of scenery Minimize the intra-day scheduling cost under is the intraday re-bidding cost of the integrated source-grid-load-storage system i at time t; , , , They are respectively the operation and maintenance cost, gas purchase cost, carbon capture cost, and demand response cost of the source-grid-load-storage integrated system i after adjusting the operating status of the unit; is the adjusted power of device k in the source-grid-load-storage integrated system i at time t; Adjust the power of gas purchase for the integrated source-grid-load-storage system i; , They are the interruptible electric / heat / hydrogen load adjustment power and the elastically transferable electric / heat / hydrogen load adjustment power of the source-grid-load-storage integrated system i at time t, respectively; The cost of the integrated source-grid-load-storage system i consuming green certificates in exchange for carbon quotas at time t; is the penalty cost of the source-grid-load-storage integrated system i at time t; The carbon quota obtained by the integrated source-grid-load-storage system i in the pre-bidding stage before time t; is the actual carbon emission of the integrated source-grid-load-storage system i during the intraday period at time t; is the subsidy coefficient; The cost incurred by the integrated source-grid-load-storage system i when transferring carbon emissions to other systems at time t; is the energy sharing income of the integrated source-grid-load-storage system i at time t; is the carbon offset adjustment of the integrated source-grid-load-storage system i at time t; is the voluntary carbon emission reduction benefit of the integrated source-grid-load-storage system i at time t; is the operation and maintenance cost coefficient of equipment k; is the price of natural gas; is the carbon capture adjustment amount of the integrated source-grid-load-storage system i at time t; is the carbon storage adjustment of the integrated source-grid-load-storage system i at time t; is the interruption load factor; is the elastic load coefficient; The green certificate quota adjustment amount for the integrated power grid-load-storage system i; is the carbon offset of the integrated source-grid-load-storage system i at time t; The actual adjustment amount of the green certificate of the integrated power grid-load-storage system i; is the green certificate price; is the wind abandonment factor; is the abandoned wind volume; is the light abandonment coefficient; is the amount of abandoned light; is the load shedding factor; is the load adjustment amount; is the carbon transfer cost coefficient; is the carbon transfer from source-grid-load-storage integrated system i to source-grid-load-storage integrated system j; The number of integrated source-grid-load-storage systems receiving energy;

[0229] The constraints of the second stage include:

[0230] Energy supply and demand balance compliance re-dispatch constraints:

[0231] ;

[0232] In the above formula, , , , , They are the adjusted outputs of photovoltaic power, wind power, storage equipment CCUS, electrolyzer EL, and gas turbine GT in the integrated source-grid-load-storage system i at time t during the intraday phase; is the adjusted power purchase amount of system i during the intraday period; is the adjusted electric power of the hydrogen fuel cell HFC at time t during the intra-day phase; , They are the adjusted thermal power of the waste heat boiler WHB and the hydrogen fuel cell HFC at time t in the intra-day stage; , They are respectively the heat storage and heat release adjustment of the thermal energy storage system at time t in the intra-day stage; , They are the storage and discharge adjustment of the electric energy storage system at time t during the intra-day phase; , They are the storage and release adjustment amounts of hydrogen in the hydrogen energy storage system at time t during the intra-day stage; , , They are the heat load, electricity load and hydrogen load adjustment at time t in the day phase respectively; , They are the adjusted hydrogen power of the hydrogen fuel cell HFC and the methane reactor MR at time t in the intra-day stage; is the adjustment amount of natural gas produced by the methane reactor MR at time t during the intra-day period; is the adjusted natural gas consumption of the gas turbine GT at time t during the intra-day phase; , , They represent the amount of electric energy, thermal energy and hydrogen used to prepare methane transferred from the integrated source-grid-load-storage system i to the integrated source-grid-load-storage system j during the day; It represents the adjustment amount of natural gas purchased by the integrated source-grid-load-storage system i during the intra-day period; represents the amount of carbon dioxide for methane production transferred from the integrated source-grid-load-storage system i to the integrated source-grid-load-storage system j; It represents the hydrogen adjustment amount of the integrated source-grid-load-storage system i in the intra-day stage;

[0233] Energy sharing balance constraints for compliance redispatch:

[0234] ;

[0235] In the above formula, , , , It represents the amount of electricity, heat, hydrogen energy and natural gas output from the integrated source-grid-load-storage system i to the integrated source-grid-load-storage system j; , , , It represents the amount of electricity, heat, hydrogen energy and natural gas received by the integrated source-grid-load-storage system j; The total number of integrated source-grid-load-storage systems that output energy; Indicates the total number of integrated source-grid-load-storage systems receiving energy;

[0236] Energy sharing price constraints for compliance and redispatch:

[0237] ;

[0238] In the above formula, Indicates the lower limit of the selling price of electricity / heat / hydrogen energy; Indicates the shared price of electricity / heat / hydrogen energy; Indicates the lower limit of the purchase price of electricity / heat / hydrogen energy;

[0239] Carbon emission flow balance constraints for compliance redispatching:

[0240] ;

[0241] In the above formula, It represents the carbon emission of CHP unit in the intra-day period; It indicates that CCUS equipment adjusts the amount of carbon captured during the intra-day period; It represents the adjusted carbon emissions consumed by the P2G equipment to produce methane after carbon capture during the day; Indicates the amount of carbon stored by CCUS equipment during the intra-day period; represents the amount of carbon dioxide transferred from the integrated source-grid-load-storage system i to the integrated source-grid-load-storage system j; Indicates the carbon emission adjustment amount during the day;

[0242] S5. In order to reduce the intraday dispatching cost as much as possible, the asymmetric Nash bargaining theory is introduced to describe the cooperative game of the alliance composed of multiple integrated power grid, load and storage systems; the total re-bid electricity and carbon quota of the integrated power grid, load and storage system alliance is equal to the sum of the expected bid electricity and expected bid carbon quota of each integrated power grid, load and storage system in the previous bid; in the alliance, the integrated power grid, load and storage systems cooperate with each other and are independent of each other, so as to minimize the intraday dispatching cost of the alliance while ensuring that their own benefits will not be damaged; the specific steps are as follows:

[0243] S51. Introduce asymmetric Nash bargaining theory to describe the cooperative game within the source-grid-load-storage integrated system alliance, and convert the multi-stage electricity-carbon market bidding strategy and benefit distribution model into the following asymmetric Nash bargaining model:

[0244] ;

[0245] In the above formula, represents the benefits obtained by the integrated source-grid-load-storage system i participating in Nash bargaining; It represents the benefits of the integrated power grid-load-storage system when it operates independently, i.e. the negotiation breakdown point; It represents the energy sharing contribution index of the source-grid-load-storage integrated system i; Indicates the total number of source-grid-load-storage integrated systems;

[0246] S52. Convert the above asymmetric Nash bargaining model into the alliance benefit maximization problem, energy sharing contribution quantification problem, and cooperation benefit distribution problem, and solve these three sub-problems in turn; the alliance benefit maximization problem includes:

[0247] ;

[0248] In the above formula, is the objective function of the second stage; The cost of transferring carbon emissions from the integrated power grid, load and storage system i to other integrated power grid, load and storage systems; Energy sharing benefits of the integrated source-grid-load-storage system i; It indicated that it would not participate in energy sharing and carbon transfer;

[0249] The energy sharing contribution quantification issues include:

[0250] ;

[0251] ;

[0252] In the above formula, represents the energy sharing contribution of the source-grid-load-storage integrated system i; , They are the energy sharing adjustment coefficient and the carbon transfer adjustment coefficient; The energy sold / maximum energy sold of the integrated source-grid-load-storage system i, including electrical energy, thermal energy, and hydrogen energy; is the purchased energy / maximum purchased energy of the integrated power grid-load-storage system i; is the carbon output / maximum carbon output of the integrated source-grid-load-storage system i; is the carbon input / maximum carbon input of the source-grid-load-storage integrated system i; The energy traded from source-grid-load-storage integrated system i to source-grid-load-storage integrated system j; The energy traded from the source-grid-load-storage integrated system j to the source-grid-load-storage integrated system i at time t; is the amount of carbon transferred from the integrated source-grid-load-storage system i to the integrated source-grid-load-storage system j at time t; is the amount of carbon transferred from the integrated source-grid-load-storage system j to the integrated source-grid-load-storage system i at time t;

[0253] The cooperation benefit distribution issues include:

[0254] ;

[0255] In the above formula, , They represent the energy sharing benefits and carbon transfer costs when the optimal solution is obtained for the alliance benefit maximization problem; is the redispatching cost.

[0256] S6. Introduce the concept of Myerson value to quantify the contribution of the integrated power generation, grid, load and storage system to the profits obtained through electricity-carbon joint bidding; calculate the Myerson value of each integrated power generation, grid, load and storage system according to the solution results, and distribute the benefits of each integrated power generation, grid, load and storage system based on the Myerson value.

[0257] The calculation formula of the Myerson value is:

[0258] ;

[0259] ;

[0260] In the above formula, is the Myerson value of the source-grid-load-storage integrated system i; Indicates transaction status; alliance It is a collection of source-grid-load-storage integrated system i. yes A subset of is the total number of samples; For the Alliance The number of Zhongyuan grid-load-storage integrated systems; Indicates the marginal contribution in cooperative games; It represents the marginal contribution in non-cooperative game, which is directly calculated without considering carbon transfer and energy sharing strategy transactions; Indicates that the source-grid-load-storage integrated system forms a cooperative game alliance Total revenue Represents the Alliance for Independent Operation of Source-Grid-Load-Storage Integrated System total revenue.

Claims

1. A dispatching and trading strategy for the integrated source-grid-load-storage system considering the synergy of the electricity-carbon market, characterized by: The scheduling trading strategy includes the following steps: S1. Establish a mathematical model of the source-grid-load-storage integrated system; The mathematical model of the source-grid-load-storage integrated system includes a CHP unit model, a P2G equipment model, an HFC hydrogen fuel cell model, a carbon capture model, the SOC state of charge of the energy storage device during charging, and the SOC state of charge of the energy storage device during discharging; S2. Establish a voluntary carbon emission reduction mechanism including day-ahead and intraday carbon offset, intraday carbon transfer and utilization, and intraday carbon emission reduction subsidies; The voluntary carbon emission reduction mechanism includes: ； ； ； ； ； ； ； ； ； In the above formula, is the day-ahead pre-bid electricity of the integrated power-grid-load-storage system i at time t; , They are the predicted power consumption of the electrolyzer EL and the storage equipment CCUS in the source-grid-load-storage integrated system i at time t; , They are respectively the day-ahead predicted charge and discharge capacity of the energy storage system i in the source-grid-load-storage integrated system at time t; is the power consumption of the electric load in the source-grid-load-storage integrated system i at time t; , , , They are the day-ahead predicted output values ​​of photovoltaic, wind power, hydrogen fuel cell and gas turbine in the source-grid-load-storage integrated system i at time t; is the electrical output power of the CHP unit; is the carbon emission of the CHP system in the integrated source-grid-load-storage system i; , They are the carbon quota benchmark values ​​for power supply and heating supply of CHP system respectively; is the energy supply correction factor of the CHP system; is the heat energy generated by the CHP system in the source-grid-load-storage integrated system i at time t; The carbon quota obtained by the integrated power generation, grid-load and storage system i in the pre-bid stage; Carbon quota obtained by using green certificates for the integrated source-grid-load-storage system; Indicates the amount of carbon quota that can be bid for the integrated source-grid-load-storage system i days ago; The surplus carbon quota of the integrated power grid-load-storage system i; represents the carbon transfer amount of the integrated source-grid-load-storage system i through intraday scheduling; , , , They are the total carbon emissions after the source-grid-load-storage integrated system i participates in carbon transfer, the carbon emissions of the CHP unit, the carbon quotas exchanged for green certificates, and the carbon quotas sold; , , They are the green certificates obtained by the source-grid-load-storage integrated system i through wind power, photovoltaic and green certificate trading markets; Represents the green certificate trading income of the source-grid-load-storage integrated system i; Represents the unit green certificate transaction price of the source-grid-load-storage integrated system i; is the carbon transfer cost of the integrated source-grid-load-storage system i; is the carbon transfer value of the integrated source-grid-load-storage system i and the integrated source-grid-load-storage system j at time t; is the carbon transfer cost coefficient; , are subsidy income and subsidy coefficient respectively; S3. Generate a set of wind and solar uncertainty scenarios based on nonparametric kernel density estimation and Wasserstein distance to determine the worst wind and solar distribution. The specific steps include: S31. Establish the expression of total prediction error of wind power and photovoltaic power; S32. Construct the probability density function of wind and solar random variables based on non-parametric kernel density estimation; S33, constructing a cumulative distribution function based on the probability density function; S34, determining the uncertainty set of the cumulative distribution function through the quantile of the wind and solar random variable; S35. Treat the cumulative distribution function as an empirical distribution and use the Wasserstein distance to describe the distance between the true distribution and the empirical distribution: S36, constructing a distributed fuzzy uncertainty set based on the Wasserstein distance, wherein the distributed fuzzy uncertainty set can be regarded as a Wasserstein ball with the empirical distribution as the center and the true distribution as the radius; S4. Multi-stage electricity-carbon market bidding strategy and benefit allocation model under the worst distribution of wind and solar power; The multi-stage electricity-carbon market bidding strategy and benefit distribution model includes the objective function and constraints of the first stage and the objective function and constraints of the second stage; the objective function of the first stage is to maximize the sum of the day-ahead stage revenue and the intraday stage scheduling cost under the worst distribution of wind and solar forecast errors, and the constraints of the first stage include the exchange power constraints of energy transactions, energy supply and demand balance constraints, carbon emission flow balance constraints, and demand response constraints; the objective function of the second stage is to minimize the intraday stage scheduling cost under the worst distribution of wind and solar forecast errors, and the constraints of the second stage include energy supply and demand balance compliance re-dispatch constraints, compliance re-dispatch energy sharing balance constraints, compliance re-dispatch energy sharing price constraints, and compliance re-dispatch carbon emission flow balance constraints; S5. Introduce asymmetric Nash bargaining theory to solve the bidding strategy and benefit distribution model of multi-stage electricity-carbon market; S6. Calculate the Myerson value of each source-grid-load-storage integrated system according to the solution results, and distribute benefits based on the Myerson value.

2. The dispatching and trading strategy of the integrated source-grid-load-storage system considering the synergy of the electricity-carbon market according to claim 1 is characterized by: The CHP unit model includes: ； ； ； In the above formula, is the output power of the gas turbine GT at time t; is the natural gas consumption of gas turbine GT at time t; The lower heating value of the gas turbine GT; is the thermal power output of the gas turbine GT at time t; , They are the power generation efficiency and heat loss rate of the gas turbine GT respectively; is the total heat energy recovered and converted by the waste heat boiler WHB at time t; is the operating efficiency of the waste heat boiler WHB; The P2G device model includes: ； ； ； In the above formula, is the hydrogen output power of electrolyzer EL at time t; is the input power of electrolyzer EL at time t; is the operating efficiency of the electrolyzer EL; is the hydrogen output power of the P2G equipment used to meet the hydrogen energy supply at time t; is the hydrogen input power of the methane reactor MR used to produce natural gas at time t; is the amount of natural gas produced by the methane reactor MR at time t; is the operating efficiency of the methane reactor MR; The HFC hydrogen fuel cell model includes: ； ； In the above formula, is the electrical output power of the hydrogen fuel cell HFC at time t; is the thermal output power of the hydrogen fuel cell HFC at time t; is the hydrogen input power of the hydrogen fuel cell HFC at time t; is the electricity conversion efficiency of the hydrogen fuel cell HFC; is the thermal conversion efficiency of hydrogen fuel cell HFC; The carbon capture model includes: ； ； ； ； ； ； ； In the above formula, is the carbon emission of the CHP system at time t; is the carbon emission intensity coefficient of the CHP system; is the carbon emissions not captured by the storage equipment CCUS at time t; The carbon capture efficiency of CCUS for storage equipment; The total carbon emissions of the integrated source-grid-load-storage system; is the carbon emission consumed by the P2G equipment to produce methane at time t; is the carbon emissions stored by the storage equipment CCUS at time t; is the carbon emissions captured by the storage device CCUS at time t; Carbon emissions transferred from the integrated source-grid-load-storage system; is the carbon consumption coefficient of methane produced by the methane reactor MR; , , They are the total energy consumption of the storage equipment CCUS at time t, the energy consumption generated by carbon capture, and the energy consumption generated by carbon storage; is the total energy consumption of the storage equipment CCUS at time t+1; , They are the unit carbon capture power consumption coefficient and carbon storage power consumption coefficient of the storage equipment CCUS respectively; , They are the upper and lower limits of CCUS operating power for storage equipment respectively; , They are the upper and lower limits of CCUS ramp rate for storage equipment, respectively; The SOC state of charge of the energy storage device during charging includes: ； In the above formula, is the charge state of the energy storage device o at time t; is the charge state of the energy storage device o at time t-1; Charging power of each energy storage device o; is the charging efficiency of the energy storage device o; is the capacity of the energy storage device o; The SOC state of charge of the energy storage device during energy discharge includes: ； In the above formula, is the energy release power of energy storage device o; is the energy release efficiency of the energy storage device o.

3. The dispatching and trading strategy of the integrated source-grid-load-storage system considering the synergy of the electricity-carbon market according to claim 2 is characterized by: The step S3 is specifically as follows: S31. Establish the expression of total prediction error of wind power and photovoltaic power: ； ； ； In the above formula, represents the total prediction error of wind power and photovoltaic power at time t in the source-grid-load-storage integrated system i; , are the prediction errors of wind power and photovoltaic power at time t in the source-grid-load-storage integrated system i; , They are the actual outputs of wind power and photovoltaic power at time t in the source-grid-load-storage integrated system i; , They are the predicted power of wind power and photovoltaic power at time t in the source-grid-load-storage integrated system i; S32. Construct the probability density function of the wind and solar random variable based on non-parametric kernel density estimation: ； In the above formula, for The probability density function of Indicates that the scenery changes randomly quantity ; Indicates the number of typical daily data samples of wind power and photovoltaic power; is bandwidth; is the kernel function; Represents the empirical distribution value of wind power and photovoltaic forecast data; S33. Construct a cumulative distribution function based on the probability density function: ； In the above formula, is the cumulative distribution function, also representing the empirical distribution; S34. Determine the uncertainty set of the cumulative distribution function through the quantile of the wind and solar random variable: ； ； In the above formula, , Respectively The upper and lower boundaries are determined by the quantile of the wind and solar random variable prediction error distribution, and the confidence level is set to ; is the uncertainty set; S35, As an empirical distribution, the Wasserstein distance is used to describe the true distribution F and the empirical distribution Distance between: ； ； In the above formula, Represents the true distribution F and the empirical distribution Wasserstein distance between them; is the supremum function; is the intermediate calculation amount; S36, constructing a distributed fuzzy uncertainty set based on Wasserstein distance, wherein the distributed fuzzy uncertainty set can be regarded as an empirical distribution The Wasserstein sphere with center and true distribution F as radius is expressed as: ； In the above formula, Represents a distributed fuzzy uncertainty set; represents the set of true distributions F; represents the radius of the Wasserstein sphere, satisfying .

4. The dispatching and trading strategy of the integrated source-grid-load-storage system considering the synergy of the electricity-carbon market according to claim 1 is characterized by: The objective function of the first stage includes: ； ； ； ； ； ； ； In the above formula, is the objective function of one stage; represents the day-ahead period return; is the objective function of the second stage; is the day-ahead electricity bidding price of the integrated power grid-load-storage system i; The estimated winning bid power of the source-grid-load-storage integrated system i; is the day-ahead carbon bidding price of the integrated power grid-load-storage system i; The expected winning carbon quota for the integrated source-grid-load-storage system i; is the pre-dispatch cost of the source-grid-load-storage integrated system i; , , , They are the operation and maintenance cost, gas purchase cost, carbon capture cost, and demand response cost of the source-grid-load-storage integrated system i; The amount of natural gas purchased for the integrated power grid, load and storage system i; is the price of natural gas; The green certificate trading volume of the integrated source-grid-load-storage system i; is the green certificate price; is the total cost of carbon capture and storage for the integrated source-grid-load-storage system i at time t; is the cost of carbon capture for the integrated source-grid-load-storage system i at time t; is the cost of carbon storage for the integrated source-grid-load-storage system i at time t; is the interruption load factor; is the elastic load coefficient; is the interruptible power of the integrated source-grid-load-storage system i at time t Load capacity; is the elastically transferable energy of the integrated source-grid-load-storage system i at time t Load capacity; It refers to electricity / heat / hydrogen; is the operation and maintenance cost coefficient of device k.

5. The dispatching and trading strategy of the integrated source-grid-load-storage system considering the synergy of the electricity-carbon market according to claim 1 is characterized by: The exchange power constraints of the energy transaction include: ； In the above formula, is the day-ahead pre-bid electricity of the integrated power-grid-load-storage system i at time t; , They are the upper and lower limits of the power of the source-grid-load-storage integrated system i respectively; is the gas consumption of the source-grid-load-storage integrated system i at time t; , are the upper and lower limits of gas consumption of the source-grid-load-storage integrated system i, respectively; The energy supply and demand balance constraints include: ； In the above formula, , , , They are the day-ahead predicted outputs of photovoltaic, wind power, storage equipment CCUS, and electrolyzer EL in the source-grid-load-storage integrated system i at time t; is the electrical output power of the hydrogen fuel cell HFC at time t; , are the heat outputs of the waste heat boiler WHB and the hydrogen fuel cell HFC at time t; , are the heat storage and release of the thermal energy storage system at time t respectively; , are the storage and discharge amounts of the electric energy storage system at time t respectively; , are the storage and release amounts of hydrogen in the hydrogen energy storage system at time t, respectively; , , They are heat load, electricity load and hydrogen load at time t respectively; , , are the hydrogen input powers of the electrolyzer EL, hydrogen fuel cell HFC, and methane reactor MR at time t, respectively; is the amount of natural gas produced by the methane reactor MR at time t; is the natural gas consumption of gas turbine GT at time t; is the output power of the gas turbine GT at time t; The carbon emission flow balance constraints include: ； In the above formula, is the carbon emission of the CHP unit at time t; is the carbon emission consumed by the P2G equipment to produce methane after carbon capture at time t; is the carbon emissions not captured by the storage equipment CCUS at time t; is the carbon emissions captured by the storage equipment CCUS at time t; is the carbon emissions consumed by the storage equipment CCUS for carbon storage at time t; The demand response constraints include: ； ； ； In the above formula, is the time t load, They represent electricity / heat / hydrogen respectively; , They represent the elastic transferability at time t. Load, interruptible load; , are elastically transferable at time t Load, interruptible The upper limit of the load; is the fixed value at time t load.

6. The dispatching and trading strategy of the integrated source-grid-load-storage system considering the synergy of the electricity-carbon market according to claim 5 is characterized by: The objective function of the two stages includes: ； ； ； ； ； ； ； ； ； ； ； In the above formula, is the objective function of the second stage, expressed in Wasserstein fuzzy sets The worst distribution of scenery The intraday scheduling cost under ; is the intraday re-bidding cost of the integrated source-grid-load-storage system i at time t; , , , They are respectively the operation and maintenance cost, gas purchase cost, carbon capture cost, and demand response cost of the source-grid-load-storage integrated system i after adjusting the operating status of the unit; is the adjusted power of device k in the source-grid-load-storage integrated system i at time t; Adjust the power of gas purchase for the integrated source-grid-load-storage system i; , They are the interruptible electric / heat / hydrogen load adjustment power and the elastically transferable electric / heat / hydrogen load adjustment power of the source-grid-load-storage integrated system i at time t, respectively; The cost of the integrated source-grid-load-storage system i consuming green certificates in exchange for carbon quotas at time t; is the penalty cost of the source-grid-load-storage integrated system i at time t; is the voluntary carbon emission reduction subsidy for the integrated power grid-load-storage system i at time t; The carbon quota obtained by the integrated source-grid-load-storage system i in the pre-bidding stage before time t; is the actual carbon emission of the integrated source-grid-load-storage system i during the intraday period at time t; is the subsidy coefficient; The cost incurred by the integrated source-grid-load-storage system i when transferring carbon emissions to other systems at time t; is the energy sharing income of the integrated source-grid-load-storage system i at time t; is the carbon offset adjustment of the integrated source-grid-load-storage system i at time t; is the voluntary carbon emission reduction benefit of the integrated source-grid-load-storage system i at time t; is the operation and maintenance cost coefficient of equipment k; is the price of natural gas; is the carbon capture adjustment amount of the integrated source-grid-load-storage system i at time t; is the carbon storage adjustment of the integrated source-grid-load-storage system i at time t; is the interruption load factor; is the elastic load coefficient; The green certificate quota adjustment amount for the integrated power grid-load-storage system i; is the carbon offset of the integrated source-grid-load-storage system i at time t; The actual adjustment amount of the green certificate of the integrated power grid-load-storage system i; is the green certificate price; is the wind abandonment factor; is the abandoned wind volume; is the light abandonment coefficient; is the amount of abandoned light; is the load shedding factor; is the load adjustment amount; is the carbon transfer cost coefficient; is the carbon transfer from source-grid-load-storage integrated system i to source-grid-load-storage integrated system j; The total number of integrated source-grid-load-storage systems that receive energy.

7. The dispatching and trading strategy of the integrated source-grid-load-storage system considering the synergy of the electricity-carbon market according to claim 6 is characterized by: The energy supply and demand balance compliance re-dispatch constraints include: ； In the above formula, , , , , They are the adjusted outputs of photovoltaic power, wind power, storage equipment CCUS, electrolyzer EL, and gas turbine GT in the integrated source-grid-load-storage system i at time t during the intraday phase; is the adjusted power purchase amount of system i during the intraday period; is the adjusted electric power of the hydrogen fuel cell HFC at time t during the intra-day phase; , They are the adjusted thermal power of the waste heat boiler WHB and the hydrogen fuel cell HFC at time t in the intra-day stage; , are the heat storage and release adjustment of the thermal energy storage system at time t in the intra-day stage respectively; , They are the storage and discharge adjustment of the electric energy storage system at time t during the intra-day phase; , They are the storage and release adjustment amounts of hydrogen in the hydrogen energy storage system at time t during the intra-day stage; , , They are the heat load, electricity load and hydrogen load adjustment at time t in the day phase respectively; , They are the adjusted hydrogen power of the hydrogen fuel cell HFC and the methane reactor MR at time t in the intra-day stage; is the adjustment amount of natural gas produced by the methane reactor MR at time t during the intra-day period; is the adjusted natural gas consumption of the gas turbine GT at time t during the intra-day phase; , , They represent the amount of electric energy, thermal energy and hydrogen used to prepare methane transferred from the integrated source-grid-load-storage system i to the integrated source-grid-load-storage system j during the day; It represents the adjustment amount of natural gas purchased by the integrated source-grid-load-storage system i during the intra-day period; represents the amount of carbon dioxide for methane production transferred from the integrated source-grid-load-storage system i to the integrated source-grid-load-storage system j; It represents the hydrogen adjustment amount of the integrated source-grid-load-storage system i in the intra-day stage; The energy sharing balance constraints of compliance redispatch include: ； In the above formula, , , , It represents the amount of electricity, heat, hydrogen energy and natural gas output from the integrated source-grid-load-storage system i to the integrated source-grid-load-storage system j; , , , It represents the amount of electricity, heat, hydrogen energy and natural gas received by the integrated source-grid-load-storage system j; The total number of integrated source-grid-load-storage systems that output energy; Indicates the total number of integrated source-grid-load-storage systems receiving energy; The compliance redispatch energy sharing price constraints include: ； In the above formula, Indicates the lower limit of the selling price of electricity / heat / hydrogen energy; Indicates the shared price of electricity / heat / hydrogen energy; Indicates the lower limit of the purchase price of electricity / heat / hydrogen energy; The compliance rescheduling carbon emission flow balance constraints include: ； In the above formula, It represents the carbon emission of CHP unit in the intra-day period; It indicates that CCUS equipment adjusts the amount of carbon captured during the intra-day period; It represents the adjusted carbon emissions consumed by the P2G equipment to produce methane after carbon capture during the day; Indicates the amount of carbon stored by CCUS equipment during the intra-day period; represents the amount of carbon dioxide transferred from the integrated source-grid-load-storage system i to the integrated source-grid-load-storage system j; Indicates the carbon emission adjustment amount during the intraday period.

8. The dispatching and trading strategy of the integrated source-grid-load-storage system considering the synergy of the electricity-carbon market according to claim 4 is characterized by: The S5 is specifically: S51. Introduce asymmetric Nash bargaining theory to describe the cooperative game within the source-grid-load-storage integrated system alliance, and convert the multi-stage electricity-carbon market bidding strategy and benefit distribution model into the following asymmetric Nash bargaining model: ； In the above formula, represents the benefits obtained by the integrated source-grid-load-storage system i participating in Nash bargaining; It represents the benefits of the integrated power grid-load-storage system when it operates independently, i.e. the negotiation breakdown point; It represents the energy sharing contribution index of the source-grid-load-storage integrated system i; Indicates the total number of source-grid-load-storage integrated systems; S52. Convert the asymmetric Nash bargaining model into the problem of maximizing alliance benefits, the problem of quantifying energy sharing contributions, and the problem of allocating cooperative benefits, and solve these three sub-problems in turn; The alliance revenue maximization problem includes: ； In the above formula, is the objective function of the second stage; The cost of transferring carbon emissions from the integrated power grid, load and storage system i to other integrated power grid, load and storage systems; Energy sharing benefits of the integrated source-grid-load-storage system i; It indicated that it would not participate in energy sharing and carbon transfer; The energy sharing contribution quantification issues include: ； ； In the above formula, , They are the energy sharing adjustment coefficient and the carbon transfer adjustment coefficient; The power sold / maximum power sold of the integrated power-grid-load-storage system i; The purchased power / maximum purchased power of the integrated power-grid-load-storage system i; is the carbon output / maximum carbon output of the integrated source-grid-load-storage system i; is the carbon input / maximum carbon input of the source-grid-load-storage integrated system i; The energy traded from source-grid-load-storage integrated system i to source-grid-load-storage integrated system j; The energy traded from the integrated source-grid-load-storage system j to the integrated source-grid-load-storage system i; is the amount of carbon transferred from the integrated source-grid-load-storage system i to the integrated source-grid-load-storage system j; is the amount of carbon transferred from the integrated source-grid-load-storage system j to the integrated source-grid-load-storage system i; The cooperation benefit distribution issues include: ； In the above formula, , They represent the energy sharing benefits and carbon transfer costs when the optimal solution is obtained for the alliance benefit maximization problem; is the redispatching cost.

9. The dispatching and trading strategy of the integrated source-grid-load-storage system considering the synergy of the electricity-carbon market according to claim 1 is characterized by: The calculation formula of the Myerson value is: ； ； In the above formula, is the Myerson value of the source-grid-load-storage integrated system i; Indicates transaction status; alliance It is a collection of source-grid-load-storage integrated systems. yes A subset of is the total number of samples of all trading strategies; For the Alliance The number of Zhongyuan grid-load-storage integrated systems; Indicates the marginal contribution in cooperative games; It represents the marginal contribution in non-cooperative game, which is directly calculated without considering carbon transfer and energy sharing strategy transactions; Indicates that the source-grid-load-storage integrated system forms a cooperative game alliance Total revenue; Represents the Alliance for Independent Operation of Source-Grid-Load-Storage Integrated System of total revenue.

Images