A dynamic aggregation method for integrated energy service providers considering node load profitability

Abstract

The present invention provides a dynamic aggregation method for integrated energy service providers that considers node load yields. The method constructs an IESP operating structure based on the IEDN. The operating structure includes a network infrastructure and EH within the service provider's jurisdiction. The network infrastructure is centered on the distribution network, and the IEDN couples the gas and heat distribution networks through the EH. The present invention also constructs an EH input-output model, performs aggregation based on the IESP aggregation principle of node yields, obtains the aggregated yield of node resources based on the initial aggregation incentive, constructs a node yield calculation model, establishes a dynamic correction model for the aggregation incentive, and uses the correction model to establish a node load aggregation potential calculation model. Finally, an IESP aggregation yield maximization model is established to obtain a method for dynamic node load aggregation. The present invention fully aggregates the electricity, gas, and heat resources within the IEDN, quantifies the real-time contribution during the aggregation process, identifies key nodes that contribute to improved efficiency, and increases aggregation yield.

Description

Technical Field

[0001] The present invention belongs to the field of integrated energy technology, and in particular relates to a dynamic aggregation method for integrated energy service providers considering node load yield rates. Background Art

[0002] With the increasing penetration of renewable energy in power grids and the maturation of technologies like distributed energy and combined heat and power, traditional power grids are evolving into local integrated energy systems that integrate distributed energy and multiple energy flows. Within the distribution network, this local integrated energy system is called an Integrated Energy Distribution Network (IEDN). Due to the geographical dispersion of distributed energy and the independence of different energy flow systems, IEDNs exhibit abundant diverse energy resources but are difficult to centrally regulate. With the development and construction of the Energy Internet, the Integrated Energy System (IES) has gradually developed as a specialized form of the Energy Internet. As managers of regional and local integrated energy systems, Integrated Energy System Providers (IESPs) can manage the equipment and loads within the IES and make appropriate adjustments to ensure optimal system operation. For IESP, since users have some flexible loads, they can flexibly adjust their energy consumption strategies based on the original simple energy purchase to obtain higher profits or reduce their own energy purchase costs. IESP can also participate in integrated energy market transactions by aggregating users' flexible loads, transforming their identity as price acceptors into market price leaders. Therefore, improving their ability to aggregate adjustable resources in the region is the key to improving operating profits.

[0003] In recent years, many advances have been made in the research on aggregation optimization of IESP. However, traditional research generally focuses on the power system, mainly on its aggregation optimization in the power system. However, in the integrated energy system environment with multi-energy integration, especially in the IEDN environment with rich distributed resources, there is a lack of research on the aggregation optimization of key node user resources by IESP, and there is also little research on the dynamic aggregation of real-time contribution of resources at different nodes.

[0004] Therefore, how to screen user-side adjustable resources with aggregation value in the complex IEDN environment, select high-quality load users based on their contribution to the trading profits of integrated energy service providers after participating in aggregation, provide them with attractive aggregation incentives, maximize user potential, and improve aggregation capabilities, is an urgent problem to be solved in the process of IESP participating in market transactions. Summary of the Invention

[0005] In response to the shortcomings of the existing technology, the present invention provides a dynamic aggregation method for integrated energy service providers that takes into account the node load yield rate, quantifies the real-time contribution of the aggregation process of electricity, gas and heat resources, identifies the node loads with the most aggregation value in different energy systems, and effectively improves the aggregation benefit.

[0006] The present invention achieves the above technical objectives through the following technical means.

[0007] A dynamic aggregation method for integrated energy service providers considering node load profitability includes the following steps:

[0008] Step 1: Build an IESP operating structure based on IEDN, including the network foundation and EH within the jurisdiction of the service provider. The network foundation is centered on the distribution network, and the IEDN is coupled to the gas and heat distribution networks through EH.

[0009] Step 2: Construct the basic EH structure. The EH converts natural gas into electricity and heat through cogeneration, and interacts with the distribution network, gas distribution network, and heat distribution network through the electricity bus, natural gas bus, and heat bus. Build the EH input and output model.

[0010] Step 3: Aggregate according to the IESP aggregation principle based on node yield;

[0011] Step 4: Obtain the aggregated income of node resources based on the initial aggregated incentives in the aggregation process of step 3, build a node yield calculation model, establish an aggregated incentive dynamic correction model, and correct the aggregated incentives in the aggregation process of step 3;

[0012] Step 5: Combined with the dynamic correction model of aggregation incentives, a node load aggregation potential calculation model is established to obtain key nodes that can improve efficiency;

[0013] Step 6: Taking the maximum aggregated benefit of IESP within the aggregation period as the goal, establish an IESP aggregated benefit maximization model;

[0014] Step 7: Perform aggregation simulation and analyze the aggregation effect.

[0015] Furthermore, in step 2, the input and output model of EH is as follows:

[0016] P EH,t =f EH,t H vg λ e

[0017] Q EH,t =f EH,t H vg λ h

[0018] Among them, PEH,t Indicates the output electrical power; f EH, t represents input energy, including electrical energy, thermal energy, and natural gas flow; H vg Indicates the lower calorific value of natural gas; e Indicates the power conversion coefficient; Q EH,t Indicates output heat energy; λ h Represents the thermal conversion coefficient.

[0019] Furthermore, in step 3, when aggregation is performed according to the IESP aggregation principle based on node yield, the electricity, gas and heat resources aggregated based on IESP are P re 、P rg 、P rh , combined with IESP's initial aggregation incentives for electricity, gas and heat resources within the IEDN, are I re,0 , I rg,0 , I rh,0 The aggregate benefits of electricity, gas and heat resources are W re 、W rg 、W rh , and then calculate the node yields of electricity, gas and heat resources as r ie 、r ig 、r ih ;

[0020] IESP modifies the aggregate incentives of electricity, gas and heat resources based on the calculated node yield and publishes them to users in IEDN. The modified aggregate incentives of electricity, gas and heat resources are I re , I rg , I rh After the user receives the modified aggregation incentive signal, the aggregation potential of the electricity, gas and heat resources generated in response is P re,max 、P rg,max 、P rh,max ; Combined with the original load L corresponding to electricity, gas and heat resources e , L g , L h , forming a power of aggregation participation in IESP.

[0021] Furthermore, in step 4, the node yield calculation model is constructed as follows:

[0022] The ratio of the single aggregation income of the aggregated unit to the overall regional aggregation income is defined as the node yield of the aggregated unit. The calculation formula for the node yield of the distribution network, gas distribution network, and heat distribution network (i.e., the node yield calculation model) is as follows:

[0023]

[0024]

[0025]

[0026] Among them, r ie,i,t 、W re,i,t They represent the node yield and aggregation benefit of the power aggregated unit located at the distribution network node i in period t; r ig,k,t 、W rg,k,t They represent the node yield and aggregated yield of the natural gas aggregated unit at node k in the gas distribution network during period t; r ih,m,t 、W rh,m,t They represent the node yield and aggregated yield of the heat aggregated unit at the node m of the heat distribution network in period t; W r,t Represents the overall aggregated return within the region during period t.

[0027] Furthermore, the W r,t The calculation formula is as follows:

[0028]

[0029] Among them, W re,i,t 、W rg,k,t 、W rh,m,t The specific calculation formula is as follows:

[0030] W re,i,t =P re,i,t (p sle,t -I re,i,t )

[0031] W rg,k,t =P rg,k,t (p slg,t -I rg,k,t )

[0032] W rh,m,t =P rh,m,t (p slh,t -I rh,m,t )

[0033] Among them, P re,i,t , I re,i,t They represent the aggregated electric power and aggregated incentive of the aggregated unit on the distribution network node i in time period t; P rg,k,t , I rg,k,t They represent the aggregated natural gas flow and aggregated incentive of the aggregated unit at the gas distribution network node k in time period t; P rh,m,t , I rh,m,t They represent the thermal power and aggregated excitation of the aggregated unit at the node m of the heat distribution network during the period t; p sle,t 、p slg,t 、pslh,t They represent the electricity price, gas price and heat price in which IESP participates in market transactions.

[0034] Furthermore, in step 4, the dynamic correction model of the aggregate incentive is as follows:

[0035] I re,i,t =I re,i,t-1 (1+r ie,i,t-1 )

[0036] I rg,k,t =I rg,k,t-1 (1+r ig,k,t-1 )

[0037] I rh,m,t =I rh,m,t-1 (1+r ih,m,t-1 )

[0038] in, I represents the aggregated incentive of the aggregated unit on the distribution network node i in the t-1 period; rg,k,t-1 I represents the aggregated excitation of the aggregated unit on the gas distribution network node k in the t-1 period; rh,m,t-1 represents the aggregated excitation of the aggregated unit at the node m of the heat distribution network in the period t-1; r ie,i,t-1 represents the node yield of the power aggregation unit located at the distribution network node i in the t-1 period; r ig,k,t-1 represents the node yield of the natural gas aggregation unit at node k in the gas distribution network during period t-1; r ih,m,t-1 It represents the node yield of the heat aggregated unit located at the node m of the heat distribution network in the period t-1.

[0039] Furthermore, the node load aggregation potential calculation model in step 5 is as follows:

[0040] Aggregation incentives will directly affect the aggregation potential of the aggregated units. Taking into account the influence of the aggregation willingness of the aggregated units, the aggregation potential is expressed as follows:

[0041] P re,i,t,max =α e I re,i,t L e,i,t

[0042] P rg,k,t,max =α g I rg,k,t L g,k,t

[0043] P rh,m,t,max =α h I rh,m,t L h,m,t

[0044] Among them, Pre,i,t,max P represents the maximum value of the aggregated electric power of the aggregated unit at the distribution network node i in the t period; rg,k,t,max P represents the maximum value of the natural gas flow that can be aggregated by the aggregated unit at the gas distribution network node k during the period t; rh,m,t,max represents the maximum value of the thermal power that can be aggregated by the aggregated unit at the node m of the heat distribution network during the period t; α e , α g , α h Respectively represent the calculation coefficients of the aggregation potential of electricity, gas and heat resources; L e,i,t represents the electrical load at node i during period t; L g,k,t represents the gas load at node k during period t; L h,m,t represents the heat load at node m during period t; I re,i,t I represents the aggregated incentive of the aggregated unit on the distribution network node i during the period t; rg,k,t I represents the aggregation incentive of the aggregated units on the gas distribution network node k in the time period t; rh,m,t It represents the aggregated excitation of the aggregated units on the heat distribution network node m during the period t.

[0045] Furthermore, in step 6, the objective is to maximize the aggregated profit F of the IESP within the aggregation period, and the objective function is established as follows:

[0046]

[0047] Among them, W r,t Represents the overall aggregated return within the region during period t.

[0048] Furthermore, the distribution network considers node power balance, unit output constraints, ramp constraints, and branch flow constraints; the gas distribution network considers pipeline flow constraints, gas source point constraints, flow balance constraints, compressor constraints, and node pressure constraints; the heat distribution network considers node flow balance, node power integration, load utilization characteristics, supply and return water temperature constraints, and pipe section heat transfer characteristics; the IESP based on the network aggregates the adjustable electricity, gas, and heat resources of different nodes within the IEDN range under the condition of meeting the IEDN operation constraints, and the aggregated energy participates in energy market transactions.

[0049] The present invention has the following beneficial effects:

[0050] The present invention can fully aggregate the electricity, gas and heat resources within the IEDN, quantify the real-time contribution of the aggregation process, identify the key nodes of the electricity, gas and heat systems within the IEDN that are more conducive to improving efficiency, effectively increase the aggregation benefits, and provide a basis for resource configuration and allocation. BRIEF DESCRIPTION OF THE DRAWINGS

[0051] Figure 1 Schematic diagram of the IESP operation structure based on IEDN according to the present invention;

[0052] Figure 2 This is a schematic diagram of the basic structure of the EH of the present invention;

[0053] Figure 3 Schematic diagram of the IESP aggregation principle based on node yield according to the present invention. DETAILED DESCRIPTION

[0054] The present invention will be further described below with reference to the accompanying drawings and specific embodiments, but the protection scope of the present invention is not limited thereto.

[0055] The method for dynamic aggregation of virtual power plants considering real-time contribution of the present invention comprises the following steps:

[0056] Step 1: Build Figure 1 The IEDN-based IESP operating structure shown in the figure includes the network infrastructure and energy hub (EH) within the service provider's jurisdiction. The network infrastructure is centered on the Electric Distribution Network (EDN), and couples the IEDN with the Natural Gas Distribution Network (NGDN) and the Heat Distribution Network (HDN) through the EH. Under this network infrastructure, the IESP mainly aggregates the adjustable electricity, gas, and heat resources at different nodes within the IEDN, while meeting the IEDN's operating constraints. The aggregated energy participates in energy market transactions.

[0057] Among them, the distribution network considers node power balance, unit output constraints, climbing constraints and branch flow constraints; the gas distribution network considers pipeline flow constraints, gas source point constraints, flow balance constraints, compressor constraints and node pressure constraints; the heating distribution network considers node flow balance, node power fusion, load utilization characteristics, supply and return water temperature constraints and pipe section heat transfer characteristics.

[0058] Step 2: Build Figure 2The basic structure of the EH is shown in the figure. The EH is an important support for the IEDN. The EH mainly converts natural gas into electricity and heat energy through combined heat and power (CHP). It also exchanges energy with the distribution network, gas distribution network, and heat distribution network through the electrical bus (EB), gas bus (GB), and heat bus (HB). On the one hand, the EH supplies the electric load (EL), gas load (GL), and heat load (HL) of the node where it is located. On the other hand, it supports energy exchange between the distribution network, gas distribution network, and heat distribution network.

[0059] Construct the input and output model of EH, which is specifically expressed as follows:

[0060] P EH,t =f EH,t H vg λ e

[0061] Q EH,t =f EH,t H vg λ h

[0062] Among them, P EH,t Indicates the output electrical power; f EH,t Represents input energy, including electrical energy, thermal energy, and natural gas flow; H vg Indicates the lower calorific value of natural gas; e Indicates the power conversion coefficient; Q EH,t Indicates output heat energy; h Represents the thermal conversion coefficient.

[0063] Step 3: Follow the steps below Figure 3 The IESP aggregation principle based on node yield is shown in the figure. The electricity, gas and heat resources based on IESP aggregation are P re 、P rg 、P rh , combined with IESP's initial aggregation incentives for electricity, gas and heat resources within the IEDN, are I re,0 , I rg,0 , I rh,0 , respectively obtain the aggregate benefits W of electricity, gas and heat resources re 、W rg 、W rh , and then calculate the node load yield of electricity, gas and heat resources as r ie 、r ig 、r ih ;

[0064] IESP modifies the aggregate incentives of electricity, gas and heat resources based on the calculated node yield and publishes them to users in IEDN. The modified aggregate incentives of electricity, gas and heat resources are I re , I rg , I rh After the user receives the modified aggregation incentive signal, the aggregation potential of the electricity, gas and heat resources generated in response is P re,max 、P rg,max 、P rh,max ; Combined with the original load L corresponding to electricity, gas and heat resources e , L g , L h , forming the power of aggregation participating in IESP; the above process must meet the operation constraints of IEDN.

[0065] Step 4: Obtain the aggregated income of node resources based on the initial aggregate incentive, build a node yield calculation model, and establish a dynamic correction model for aggregate incentives. The specific process is as follows:

[0066] The ratio of the single aggregation income of the aggregated unit to the overall regional aggregation income is defined as the node load yield rate of the aggregated unit. The calculation formula for the node load yield rate of the distribution network, gas distribution network, and heat distribution network (i.e., the node yield rate calculation model) is as follows:

[0067]

[0068]

[0069]

[0070] Among them, r ie,i,t 、W re,i,t They represent the node yield and aggregation benefit of the power aggregated unit located at the distribution network node i in period t; r ig,k,t 、W rg,k,t They represent the node yield and aggregated yield of the natural gas aggregated unit at node k in the gas distribution network during period t; r ih,m,t 、W rh,m,t They represent the node yield and aggregation yield of the heat aggregated unit at the node m of the heat distribution network in period t; W r,t Represents the overall aggregated return within the region during period t.

[0071] The W r,t The calculation formula is as follows:

[0072]

[0073] Among them, W re,i,t 、Wrg,k,t 、W rh,m,t The specific calculation formula is as follows:

[0074] W re,i,t =P re,i,t (p sle,t -I re,i,t )

[0075] W rg,k,t =P rg,k,t (p slg,t -I rg,k,t )

[0076] W rh,m,t =P rh,m,t (p slh,t -I rh,m,t )

[0077] Among them, P re,i,t , I re,i,t They represent the aggregated electric power and aggregated incentive of the aggregated unit on the distribution network node i in time period t; P rg,k,t , I rg,k,t They represent the aggregated natural gas flow and aggregated incentive of the aggregated unit at the gas distribution network node k in time period t; P rh,m,t , I rh,m,t They represent the thermal power and aggregated excitation of the aggregated unit at the node m of the heat distribution network during the period t; p sle,t 、p slg,t 、p slh,t They represent the electricity price, gas price and heat price in which IESP participates in market transactions;

[0078] Based on the distribution network, gas distribution network, and heat distribution network, the calculation formula for the aggregated incentive based on the node yield rate (i.e., the dynamic correction model for the aggregated incentive) is as follows:

[0079] I re,i,t =I re,i,t-1 (1+r ie,i,t-1 )

[0080] I rg,k,t =I rg,k,t-1 (1+r ig,k,t-1 )

[0081] I rh,m,t =I rh,m,t-1 (1+r ih,m,t-1 ).

[0082] Step 5: Combined with the dynamic correction model of aggregation incentives, a node load aggregation potential calculation model is established. The specific process is as follows:

[0083] Aggregation incentives will directly affect the aggregation potential of the aggregated units. Taking into account the influence of the aggregation willingness of the aggregated units, this application expresses the aggregation potential as follows:

[0084] P re,i,t,max =α e I re,i,t L e,i,t

[0085] P rg,k,t,max =α g I rg,k,t L g,k,t

[0086] P rh,m,t,max =α h I rh,m,t L h,m,t

[0087] Among them, P re,i,t,max P represents the maximum value of the aggregated electric power of the aggregated unit at the distribution network node i in the t period; rg,k,t,max P represents the maximum value of the natural gas flow that can be aggregated by the aggregated unit at the gas distribution network node k during the period t; rh,m,t,max represents the maximum value of the thermal power that can be aggregated by the aggregated unit at the node m of the heat distribution network during the period t; α e , α g , α h Respectively represent the calculation coefficients of the aggregation potential of electricity, gas and heat resources; L e,i,t represents the electrical load at node i during period t; L g,k,t represents the gas load at node k during period t; L h,m,t represents the heat load at node m during period t.

[0088] Step 6: Establish an IESP aggregation benefit maximization model to obtain a node load dynamic aggregation method. The goal is to maximize the IESP aggregation benefit F within the aggregation period. The established objective function is as follows:

[0089]

[0090] Step 7: Conduct aggregation simulation on the dynamic aggregation method of integrated energy service providers considering node load yield and analyze the aggregation effect.

[0091] The embodiments described are preferred implementations of the present invention, but the present invention is not limited to the above implementations. Any obvious improvements, substitutions or modifications that can be made by those skilled in the art without departing from the essence of the present invention are within the scope of protection of the present invention.

Claims

1. A dynamic aggregation method for integrated energy service providers considering node load profitability, characterized in that: The steps include: Step 1: Build an IESP operating structure based on IEDN, including the network foundation and EH within the jurisdiction of the service provider. The network foundation is centered on the distribution network, and the IEDN is coupled to the gas and heat distribution networks through EH. Step 2: Construct the basic EH structure. The EH converts natural gas into electricity and heat through cogeneration, and interacts with the distribution network, gas distribution network, and heat distribution network through the electricity bus, natural gas bus, and heat bus. Build the EH input and output model. Step 3: Aggregate according to the IESP aggregation principle based on node yield; Step 4: Obtain the aggregated income of node resources based on the initial aggregated incentives in the aggregation process of step 3, build a node yield calculation model, establish an aggregated incentive dynamic correction model, and correct the aggregated incentives in the aggregation process of step 3; Step 5: Combined with the dynamic correction model of aggregation incentives, a node load aggregation potential calculation model is established to obtain key nodes that can improve efficiency; Step 6: Taking the maximum aggregated benefit of IESP within the aggregation period as the goal, establish an IESP aggregated benefit maximization model; Step 7: Perform aggregation simulation and analyze the aggregation effect.

2. The method for dynamic aggregation of integrated energy service providers considering node load yield according to claim 1 is characterized in that: In step 2, the input and output model of EH is as follows: P EH,t =f EH,t H vg λ e Q EH,t =f EH,t H vg λ h Among them, P EH,t Indicates the output electrical power; f EH,t Represents input energy, including electrical energy, thermal energy, and natural gas flow; H vg Indicates the lower calorific value of natural gas; e Indicates the power conversion coefficient; Q EH,t Indicates output heat energy; h Represents the thermal conversion coefficient.

3. The method for dynamic aggregation of integrated energy service providers considering node load yield according to claim 1 is characterized in that: In step 3, when aggregation is performed according to the IESP aggregation principle based on node yield, the electricity, gas and heat resources aggregated based on IESP are P re 、P rg 、P rh , combined with IESP's initial aggregation incentives for electricity, gas and heat resources within the IEDN, are I re,0 , I rg,0 , I rh,0 The aggregate benefits of electricity, gas and heat resources are W re 、W rg 、W rh , and then calculate the node yields of electricity, gas and heat resources as r ie 、r ig 、r ih ; IESP modifies the aggregate incentives of electricity, gas and heat resources based on the calculated node yield and publishes them to users in IEDN. The modified aggregate incentives of electricity, gas and heat resources are I re , I rg , I rh After the user receives the modified aggregation incentive signal, the aggregation potential of the electricity, gas and heat resources generated in response is P re,max 、P rg,max 、P rh,max ; Combined with the original load L corresponding to electricity, gas and heat resources e , L g , L h , forming a power of aggregation participation in IESP.

4. The method for dynamic aggregation of integrated energy service providers considering node load yield according to claim 1 is characterized in that: In step 4, the method for constructing the node yield calculation model is as follows: The ratio of the single aggregation income of the aggregated unit to the overall regional aggregation income is defined as the node yield of the aggregated unit. The calculation formula of the node yield based on the distribution network, gas distribution network, and heat distribution network is as follows: Among them, r ie,i,t 、W re,i,t They represent the node yield and aggregation benefit of the power aggregated unit located at the distribution network node i in period t; r ig,k,t 、W rg,k,t They represent the node yield and aggregated yield of the natural gas aggregated unit at node k in the gas distribution network during period t; r ih,m,t 、W rh,m,t They represent the node yield and aggregated yield of the heat aggregated unit at the node m of the heat distribution network in period t; W r,t Represents the overall aggregated return within the region during period t.

5. The method for dynamic aggregation of integrated energy service providers considering node load yield according to claim 4 is characterized in that: The W r,t The calculation formula is as follows: Among them, W re,i,t 、W rg,k,t 、W rh,m,t The specific calculation formula is as follows: W re,i,t =P re,i,t (p sle,t -I re,i,t ) W rg,k,t =P rg,k,t (p slg,t -I rg,k,t ) W rh,m,t =P rh,m,t (p slh,t -I rh,m,t ) Among them, P re,i,t , I re,i,t They represent the aggregated electric power and aggregated incentive of the aggregated unit on the distribution network node i in time period t; P rg,k,t , I rg,k,t They represent the aggregated natural gas flow and aggregated incentive of the aggregated unit at the gas distribution network node k in time period t; P rh,m,t , I rh,m,t They represent the thermal power and aggregated excitation of the aggregated unit at the node m of the heat distribution network during the period t; p sle,t 、p slg,t 、p slh,t They represent the electricity price, gas price and heat price in which IESP participates in market transactions.

6. The method for dynamic aggregation of integrated energy service providers considering node load yield according to claim 5 is characterized in that: In step 4, the dynamic correction model of the aggregate incentive is as follows: I re,i,t =I re,i,t-1 (1+r ie,i,t-1 ) I rg,k,t =I rg,k,t-1 (1+r ig,k,t-1 ) I rh,m,t =I rh,m,t-1 (1+r ih,m,t-1 ) Among them, I re,i,t-1 I represents the aggregated incentive of the aggregated unit on the distribution network node i in the t-1 period; rg,k,t-1 I represents the aggregated excitation of the aggregated unit on the gas distribution network node k in the t-1 period; rh,m,t-1 represents the aggregated excitation of the aggregated unit at the node m of the heat distribution network in the period t-1; r ie,i,t-1 represents the node yield of the power aggregation unit located at the distribution network node i in the t-1 period; r ig,k,t-1 represents the node yield of the natural gas aggregation unit at the gas distribution network node k in the t-1 period; r ih,m,t-1 It represents the node yield of the heat aggregated unit located at the node m of the heat distribution network in the t-1 period.

7. The method for dynamic aggregation of integrated energy service providers considering node load profitability according to claim 1, characterized in that: The calculation model of node load aggregation potential in step 5 is as follows: Aggregation incentives will directly affect the aggregation potential of the aggregated units. Taking into account the influence of the aggregation willingness of the aggregated units, the aggregation potential is expressed as follows: P re,i,t,max =a e I re,i,t L e,i,t P rg,k,t,max =a g I rg,k,t L g,k,t P rh,m,t,max =a h I rh,m,t L h,m,t Among them, P re,i,t,max P represents the maximum value of the aggregated electric power of the aggregated unit at the distribution network node i in the t period; rg,k,t,max P represents the maximum value of the natural gas flow that can be aggregated by the aggregated unit at the gas distribution network node k during the period t; rh,m,t,max represents the maximum value of the thermal power that can be aggregated by the aggregated unit at the node m of the heat distribution network during the period t; α e , α g , α h Respectively represent the calculation coefficients of the aggregation potential of electricity, gas and heat resources; L e,i,t represents the electrical load at node i during period t; L g,k,t represents the gas load at node k during period t; L h,m,t represents the heat load at node m during period t; I re,i,t I represents the aggregated incentive of the aggregated unit on the distribution network node i during the period t; rg,k,t I represents the aggregation incentive of the aggregated units on the gas distribution network node k in the time period t; rh,m,t It represents the aggregated excitation of the aggregated units on the heat distribution network node m during the period t.

8. The method for dynamic aggregation of integrated energy service providers considering node load profitability according to claim 1, characterized in that: In step 6, the objective is to maximize the aggregated benefit F of the IESP within the aggregation period, and the objective function established is as follows: Among them, W r,t Represents the overall aggregated return within the region during period t.

9. The method for dynamic aggregation of integrated energy service providers considering node load profitability according to claim 1, characterized in that: The distribution network considers node power balance, unit output constraints, ramp constraints, and branch flow constraints; the gas distribution network considers pipeline flow constraints, gas source point constraints, flow balance constraints, compressor constraints, and node pressure constraints; the heat distribution network considers node flow balance, node power integration, load utilization characteristics, supply and return water temperature constraints, and pipe section heat transfer characteristics; the IESP based on the network aggregates the adjustable electricity, gas, and heat resources at different nodes within the IEDN range under the condition of meeting the IEDN operation constraints, and the aggregated energy participates in energy market transactions.

Images