A carbon flow calculation method for integrated energy system pool trading network

Abstract

A carbon flow calculation method for a joint trading network of an integrated energy system includes the following steps: converting the joint trading in the mixed trading of the integrated energy system into multiple virtual bilateral tradings to obtain an integrated energy system containing only bilateral trading; calculating the network loss shared by each bilateral trading, and calculating the total network loss that the joint trading needs to bear; evenly distributing the total network loss to both sides of the virtual bilateral trading, and obtaining the DC parameters of the generator and load in the adjusted joint trading network; calculating the node active flux of the joint trading network, and calculating the node carbon potential of each node according to the node active flux; calculating the carbon flow distribution of the joint trading network and the load carbon flow rate according to the node carbon potential. The present invention proposes to realize the decoupling of multiple trading networks in the mixed trading of the integrated energy system, accurately calculate the carbon flow distribution of the entire integrated energy system, and clarify the influence of different types of power trading on the carbon flow distribution of the entire integrated energy system.

Description

Technical Field

[0001] The present invention relates to the field of integrated energy systems, and in particular to a carbon flow calculation method for an integrated energy system joint venture trading network. Background Art

[0002] In the process of the power industry's transformation towards low-carbonization, the calculation method of carbon emission flows between and within various integrated energy systems plays a vital role. In the previous hybrid power trading process, it is generally believed that carbon emissions only exist in the power generation link, and the transmission, distribution and use of electricity will not produce carbon emissions. It is very unfair to ask the power generation link to bear all carbon emission responsibilities.

[0003] At present, the electricity market is mainly divided into joint trading mode and bilateral trading mode. If there is a mixed electricity market of joint trading and bilateral trading in the network, the entire integrated energy system is not linear due to the existence of network losses, that is, the sum of the network losses generated when each transaction acts on the system alone is not equal to the total network losses generated by all transactions acting on the system together. Therefore, when calculating the carbon flow distribution of the power system involving the mixed trading mode of the integrated energy system, the existence of network loss coupling terms between the transactions makes the calculation difficult to achieve.

[0004] Therefore, in order to clarify the impact of joint trading on the carbon flow distribution of the system in the integrated energy system, it is urgent to propose a carbon flow calculation method for the joint trading network of the integrated energy system to solve the problems of unclear allocation of coupling network losses between different transactions and unquantifiable carbon flow distribution. Summary of the invention

[0005] The purpose of the present invention is to overcome the problems in the prior art of unclear allocation of coupling network losses between different transactions in a hybrid trading network of an integrated energy system and inability to quantify carbon flow distribution, and to provide a carbon flow calculation method for a joint trading network of an integrated energy system.

[0006] To achieve the above objectives, the technical solution of the present invention is: a carbon flow calculation method for an integrated energy system joint trading network, comprising the following steps:

[0007] Convert the pooled transactions in the mixed transactions of the integrated energy system into multiple virtual bilateral transactions, and obtain an integrated energy system that only includes bilateral transactions;

[0008] Calculate the network loss shared by each bilateral transaction and the total network loss that needs to be borne by the joint transaction;

[0009] The total network loss is evenly distributed to both sides of the virtual bilateral transaction, and the DC parameters of the generators and loads in the adjusted pooled transaction network are obtained;

[0010] Calculate the node active flux of the pooled trading network, and calculate the node carbon potential of each node based on the node active flux;

[0011] The carbon flow distribution and load carbon flow rate of the joint trading network are calculated based on the node carbon potential.

[0012] Assume that there are n bilateral transactions and pooled transactions in the hybrid transaction of the integrated energy system, where the DC parameters of the pooled transaction network are:

[0013]

[0014] Where: G mL , LD mL are the power generation and load of node m in the comprehensive energy joint trading network; G m , LD m are the power generation and load of node m in the hybrid transaction of the integrated energy system; m = 1, ..., a, a is the number of nodes in the hybrid transaction of the integrated energy system; G mB is the power generation of node m in the comprehensive energy bilateral trading network; LD mB is the load of node m in the comprehensive energy bilateral trading network;

[0015] If the joint venture transaction is converted into s virtual bilateral transactions, then there are n+s bilateral transactions in the integrated energy system;

[0016] The transaction volume of each virtual bilateral transaction is:

[0017]

[0018] Where: B n+i is the transaction volume of the i-th virtual bilateral transaction; is a column vector whose kth component is 1 and the rest are 0; e i is a column vector with the i-th component being 1 and the other components being 0, k and i represent the node positions of the load and the generator; P k is the power flowing through load node k; LD k is the load of load node k in the mixed transaction of the integrated energy system; LD kB is the load of load node k in the bilateral trading network of the comprehensive energy system; G i G is the power generation of generator node i in the mixed transaction of the integrated energy system; iB is the power generation of generator node i in the bilateral trading network of the integrated energy system; is the upstream distribution matrix.

[0019] The integrated energy system is split into r sub-integrated energy systems, and the network equivalents of the r sub-integrated energy systems are converted into a connected graph containing only tree branches;

[0020] Remove the redundant interconnection line flows between different sub-integrated energy systems, so that there is only one interconnection line between different sub-integrated energy systems, and then an equivalent circulation is formed between different sub-integrated energy systems;

[0021] The relationship between each branch flow and each bilateral transaction in the integrated energy system 1 for:

[0022] F 1 =C cb B;

[0023] Where: C cb =C a C b is the branch transaction association matrix, which is an l×(n+s)-order matrix, C a is the branch incidence matrix, C b is the transaction association matrix; l is the number of branches of the mixed transaction of the integrated energy system.

[0024] The network loss shared by each bilateral transaction is:

[0025]

[0026] Where: C cb，j C cb The (n×s)-order row vector of the j-th row of the matrix, j=1,...,l, where l is the number of branches of the hybrid transaction of the integrated energy system; P loss，j is the branch loss on the jth branch; B i is a (n×s)-order column vector whose ith element is the transaction amount of the ith bilateral transaction and whose other elements are 0. i =[0, ..., B i , ...0];

[0027] The total network loss that the joint venture transaction needs to bear is:

[0028]

[0029] The DC parameters of generators and loads in the pooled trading network in hybrid trading are:

[0030]

[0031] Where: are respectively the active output of the balanced units in the joint trading network, the active output of the unbalanced units and the active input of the load; P G1，j , P G2 , j , LD iThey are respectively the active output of the balancing units, the active output of the unbalanced units and the active input of the loads in the mixed transaction of the integrated energy system; The carbon emission intensity of the units in the pool trading network; They are respectively the active output of the unbalanced units and the active input of the loads in the bilateral trading network;

[0032] The total network loss is evenly distributed to both sides of the virtual bilateral transaction, and the DC parameters of the generator and load in the adjusted pooled transaction network are obtained as follows:

[0033]

[0034] Where: are the DC parameters of the generator and load in the adjusted pooling trading network; B s represents the set of virtual bilateral transactions; Gi is the i-th generator; They are used to adjust the DC parameters of generators and loads in the pre-pooling trading network respectively.

[0035] Calculate the DC branch power flow distribution matrix, unit injection distribution matrix and load distribution matrix of the pool trading network;

[0036] The node active flux is obtained based on the splicing matrix of the DC branch power flow distribution matrix and the unit injection distribution matrix. for:

[0037]

[0038] Where: ξ is a row vector of order l+h, and all elements are 1, l is the number of branches of the integrated energy system hybrid trading, and h is the number of generator sets of the integrated energy system hybrid trading; It is the concatenation matrix of the DC branch power flow distribution matrix and the unit injection distribution matrix. is the DC branch power flow distribution matrix, Inject distribution matrix into the unit;

[0039] The node carbon potential of each node is obtained through the node active flux of the joint trading network, the DC branch power flow distribution matrix and the unit injection distribution matrix. for:

[0040]

[0041] The carbon flow distribution and load carbon flow rate of the joint trading network are:

[0042]

[0043] Where: They are the load carbon flow rate distribution matrix and branch carbon flow rate distribution matrix of the joint trading network respectively; LD T，after is the load distribution matrix.

[0044] A carbon flow calculation device for a comprehensive energy system joint trading network, the device is applied to the above-mentioned method, the device comprises:

[0045] A transaction conversion module is used to convert the joint transaction in the mixed transaction of the integrated energy system into multiple virtual bilateral transactions to obtain an integrated energy system containing only bilateral transactions;

[0046] The network loss calculation module is used to calculate the network loss shared by each bilateral transaction and the total network loss that needs to be borne by the joint transaction;

[0047] The DC parameter calculation module is used to evenly distribute the total network loss to both sides of the virtual bilateral transaction and obtain the DC parameters of the generators and loads in the adjusted pooled transaction network;

[0048] The carbon flow distribution calculation module is used to calculate the node active flux of the joint trading network, and calculate the node carbon potential of each node based on the node active flux, and calculate the carbon flow distribution and load carbon flow rate of the joint trading network based on the node carbon potential.

[0049] A carbon flow computing device for an integrated energy system pool trading network, comprising a memory and a processor;

[0050] The memory is used to store computer program code and transmit the computer program code to the processor;

[0051] The processor is used to execute the method as described above according to the instructions in the computer program code.

[0052] A computer-readable storage medium stores a computer program, and when the computer program is executed by a processor, the method described above is implemented.

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

[0054] The present invention provides a carbon flow calculation method for an integrated energy system joint trading network. Aiming at the problems of unclear allocation of coupled network losses between different transactions and unquantifiable carbon flow distribution in a hybrid trading network of an integrated energy system, it is proposed that in an integrated energy system with an expanding scale of transactions between multiple regions, the carbon flow distribution of the entire integrated energy system can be accurately calculated, the influence of different types of power transactions on the carbon flow distribution of the entire integrated energy system can be clarified, the carbon flow direction of transactions in the system can be clearly reflected, and the influence of transaction order on the system is ignored, so that the network loss is allocated to the hybrid transaction in a more fair and reasonable manner. BRIEF DESCRIPTION OF THE DRAWINGS

[0055] Figure 1 It is a flow chart of a carbon flow calculation method for a comprehensive energy system joint trading network of the present invention.

[0056] Figure 2 It is a structural block diagram of a carbon flow calculation device for a comprehensive energy system joint trading network according to the present invention.

[0057] Figure 3 It is a structural block diagram of a carbon flow calculation device for a comprehensive energy system joint trading network according to the present invention. DETAILED DESCRIPTION

[0058] The present invention is further described in detail below in conjunction with the accompanying drawings and specific implementation methods.

[0059] Embodiment 1:

[0060] See also Figure 1 , a carbon flow calculation method for an integrated energy system pool trading network, comprising the following steps:

[0061] S1. Convert the joint transaction in the mixed transaction of the integrated energy system into multiple virtual bilateral transactions to obtain an integrated energy system containing only bilateral transactions;

[0062] Assume that there are joint transactions and n bilateral transactions in the mixed transaction of the integrated energy system:

[0063] B1: (G 1B LD 1B ),...,Bn:(G nB LD nB );

[0064] That is, there are n+1 transactions in the integrated energy system. In the DC power flow of the mixed transaction of the integrated energy system, after removing the bilateral transactions, the DC parameters of the joint transaction network are:

[0065]

[0066] Where: G mL , LD mL are the power generation and load of node m in the comprehensive energy joint trading network; G m , LD m are the power generation and load of node m in the hybrid transaction of the integrated energy system; m = 1, ..., a, a is the number of nodes in the hybrid transaction of the integrated energy system; G mB is the power generation of node m in the comprehensive energy bilateral trading network; LD mB is the load of node m in the comprehensive energy bilateral trading network; B is the transaction volume of bilateral transactions;

[0067] If the joint venture transaction is converted into S virtual bilateral transactions, there will be n+s bilateral transactions in the integrated energy system;

[0068] The transaction volume of each virtual bilateral transaction is:

[0069]

[0070] Where: B n+i is the transaction volume of the i-th virtual bilateral transaction; is a column vector whose kth component is 1 and the rest are 0; e i is a column vector with the i-th component being 1 and the other components being 0, k and i represent the node positions of the load and the generator; P k is the power flowing through load node k; LD k is the load of load node k in the mixed transaction of the integrated energy system; LD kB is the load of load node k in the bilateral trading network of the comprehensive energy system; G i G is the power generation of generator node i in the mixed transaction of the integrated energy system; iB is the power generation of generator node i in the bilateral trading network of the integrated energy system; is the upstream distribution matrix;

[0071] The supply and demand electricity of the virtual bilateral transaction is determined by the transaction electricity of each virtual bilateral transaction, and the original system is changed into a system with only n+s bilateral transactions (among which, the n+1th to n+sth transactions are virtual bilateral transactions). Since the processed integrated energy system only has bilateral transactions, it is easier to distribute the network loss to each bilateral transaction, thereby realizing the decoupling of the hybrid transaction mode.

[0072] S2. Calculate the network loss shared by each bilateral transaction and the total network loss that needs to be borne by the joint transaction;

[0073] The integrated energy system is split into r sub-integrated energy systems, and the network equivalent of the r sub-integrated energy systems is converted into a connected graph containing only tree branches. After the internal processing of each sub-region, the redundant connecting line flow between different sub-integrated energy systems is removed, so that there is only a unique connecting line between different sub-integrated energy systems. The place where the redundant connecting line is removed forms a corresponding equivalent circulation in the corresponding loop. Removing the connecting line between any two integrated energy systems will form an equivalent circulation between different integrated energy systems.

[0074] The relationship between the original branch flow and the flow distribution caused by transactions is F 1 :

[0075] F 1 =C a F2 ;

[0076] Where: C a is the branch correlation matrix; F 2 It is the power flow through the branches of the system after the network split, which only includes bilateral transactions within the sub-integrated energy system and bilateral transactions between different integrated energy systems;

[0077] Assume that there are n bilateral transactions in the integrated energy system B = [B 1 , B 2 , ..., B n ], since only tree branches exist in the sub-integrated energy system after the network splits, and only a unique connection line exists between different integrated energy systems, the flow of transaction volume is determined, and the relationship between the two is proposed;

[0078] F 2 =C b B;

[0079] Where: C b represents the transaction association matrix;

[0080] Since the processed integrated energy system only contains (n+s) bilateral transactions, the relationship between each branch flow and each bilateral transaction in the integrated energy system is 1 for:

[0081] F 1 =C cb B;

[0082] Where: C cb =C a C b is the branch transaction association matrix, which is an l×(n+s)-order matrix; l is the number of branches in the mixed transaction of the integrated energy system;

[0083] The network loss shared by each bilateral transaction is:

[0084]

[0085] Where: C cb，j C cb The (n×s)-order row vector of the j-th row of the matrix, j=1,...,l, where l is the number of branches of the hybrid transaction of the integrated energy system; P loss，j is the branch loss on the jth branch; B i is a (n×s)-order column vector whose jth element is the transaction amount of the i-th bilateral transaction and whose other elements are 0. i =[0, ..., B i , ...0];

[0086] The total network loss that the joint venture transaction needs to bear is:

[0087]

[0088] S3, the total network loss is evenly distributed to both sides of the virtual bilateral transaction, and the DC parameters of the generators and loads in the adjusted pooled transaction network are obtained;

[0089] The parameters of each line and network structure of the joint trading network are consistent with the mixed trading of the integrated energy system. The DC parameters of the generator and load in the joint trading network are:

[0090]

[0091] Where: are respectively the active output of the balanced units in the joint trading network, the active output of the unbalanced units and the active input of the load; P G1，j , P G2，j , LD i They are respectively the active output of the balancing units, the active output of the unbalanced units and the active input of the loads in the mixed transaction of the integrated energy system; The carbon emission intensity of the units in the pool trading network; They are respectively the active output of the unbalanced units and the active input of the loads in the bilateral trading network;

[0092] Since the joint trading is conducted through the power pool, the two parties of the transaction are not clear. According to the method of virtual bilateral trading, the network loss is evenly distributed to both sides of the virtual bilateral transaction, and the DC parameters of the generator and load in the joint trading network after adjustment are obtained as follows:

[0093]

[0094] Where: are the DC parameters of the generator and load in the adjusted pooling trading network; B s represents the set of virtual bilateral transactions; Gi is the i-th generator; They are used to adjust the DC parameters of generators and loads in the pre-pooling trading network respectively.

[0095] S4, calculating the node active flux of the pooled trading network, and calculating the node carbon potential of each node according to the node active flux;

[0096] The net power generation and net load input of the pooling trading network are known, and the line parameters in the network remain unchanged. The DC branch power flow distribution matrix, unit injection distribution matrix and load distribution matrix of the pooling trading network are calculated by the DC power flow calculation method.

[0097] The node active flux is obtained based on the splicing matrix of the DC branch power flow distribution matrix and the unit injection distribution matrix. for:

[0098]

[0099] Where: ξ is a row vector of order l+h, and all elements are 1, l is the number of branches of the integrated energy system hybrid trading, and h is the number of generator sets of the integrated energy system hybrid trading; It is the concatenation matrix of the DC branch power flow distribution matrix and the unit injection distribution matrix. is the DC branch power flow distribution matrix, Inject distribution matrix into the unit;

[0100] The node carbon potential of each node is obtained through the node active flux of the joint trading network, the DC branch power flow distribution matrix and the unit injection distribution matrix. for:

[0101]

[0102] S5. Calculate the carbon flow distribution and load carbon flow rate of the joint trading network based on the node carbon potential.

[0103] The carbon flow distribution and load carbon flow rate of the joint trading network are:

[0104]

[0105] Where: They are the load carbon flow rate distribution matrix and branch carbon flow rate distribution matrix of the joint trading network respectively; LD T，ofter is the load distribution matrix.

[0106] Embodiment 2:

[0107] See also Figure 2 , a carbon flow calculation device for a comprehensive energy system joint trading network, the device is applied to the method described in Example 1, the device includes;

[0108] A transaction conversion module is used to convert the joint transaction in the mixed transaction of the integrated energy system into multiple virtual bilateral transactions to obtain an integrated energy system containing only bilateral transactions;

[0109] The network loss calculation module is used to calculate the network loss shared by each bilateral transaction and the total network loss that needs to be borne by the joint transaction;

[0110] The DC parameter calculation module is used to evenly distribute the total network loss to both sides of the virtual bilateral transaction and obtain the DC parameters of the generators and loads in the adjusted pooled transaction network;

[0111] The carbon flow distribution calculation module is used to calculate the node active flux of the joint trading network, and calculate the node carbon potential of each node based on the node active flux, and calculate the carbon flow distribution and load carbon flow rate of the joint trading network based on the node carbon potential.

[0112] Embodiment 3:

[0113] See also Figure 3 , a carbon flow calculation device for an integrated energy system joint trading network, comprising a memory and a processor; the memory is used to store computer program code and transmit the computer program code to the processor; the processor is used to execute the method described in Example 1 according to the instructions in the computer program code.

[0114] A computer-readable storage medium having a computer program stored thereon, wherein the computer program implements the method described in Example 1 when executed by a processor.

[0115] Generally speaking, the computer instructions for implementing the method of the present invention may be carried in any combination of one or more computer-readable storage media. Non-transitory computer-readable storage media may include any computer-readable media, except for the signal itself that is temporarily propagating.

[0116] The computer-readable storage medium may be, for example, but not limited to, an electrical, magnetic, optical, electromagnetic, infrared, or semiconductor system, device or device, or any combination thereof. More specific examples of computer-readable storage media (a non-exhaustive list) include: an electrical connection with one or more wires, a portable computer disk, a hard disk, a random access memory (RAM), a read-only memory (ROM), an erasable programmable read-only memory (EKROM or flash memory), an optical fiber, a portable compact disk read-only memory (CD-ROM), an optical storage device, a magnetic storage device, or any suitable combination thereof. In the present invention, a computer-readable storage medium may be any tangible medium containing or storing a program that may be used by or in conjunction with an instruction execution system, device, or device.

[0117] Computer program code for performing the operation of the present invention can be written in one or more programming languages ​​or a combination thereof, including object-oriented programming languages ​​such as Java, Smalltalk, C++, and conventional procedural programming languages ​​such as "C" or similar programming languages, in particular, Python suitable for neural network computing and platform frameworks based on TensorFlow, PyTorch, etc. can be used. The program code can be executed entirely on the user's computer, partially on the user's computer, as an independent software package, partially on the user's computer and partially on a remote computer, or entirely on a remote computer or server. In the case of a remote computer, the remote computer can be connected to the user's computer or to an external computer (for example, using an Internet service provider to connect via the Internet) through any type of network, including a local area network (LAN) or a wide area network (WAN).

[0118] The above-mentioned device and non-temporary computer-readable storage medium can be found in the specific description of a carbon flow calculation method for an integrated energy system joint trading network and its beneficial effects, which will not be repeated here.

[0119] Although the embodiments of the present invention have been shown and described above, it should be understood that the above embodiments are exemplary and should not be construed as limitations on the present invention. A person skilled in the art may change, modify, replace and vary the above embodiments within the scope of the present invention.

Claims

1. A carbon flow calculation method for a comprehensive energy system joint trading network, characterized in that: The method comprises the following steps: Convert the joint trading in the mixed trading of the integrated energy system into multiple virtual bilateral trading to obtain an integrated energy system that only includes bilateral trading; Assume that there are n bilateral transactions and pooled transactions in the hybrid transaction of the integrated energy system, where the DC parameters of the pooled transaction network are: Where: G mL , LD mL are the power generation and load of node m in the comprehensive energy joint trading network; G m , LD m are the power generation and load of node m in the hybrid transaction of the integrated energy system; m = 1, ..., a, a is the number of nodes in the hybrid transaction of the integrated energy system; G mB is the power generation of node m in the comprehensive energy bilateral trading network; LD mB is the load of node m in the comprehensive energy bilateral trading network; If the joint venture transaction is converted into S virtual bilateral transactions, there will be n+s bilateral transactions in the integrated energy system; The transaction volume of each virtual bilateral transaction is: Where: B n+i is the transaction volume of the i-th virtual bilateral transaction; is a column vector whose kth component is 1 and the rest are 0; e i is a column vector with the i-th component being 1 and the other components being 0, k and i represent the node positions of the load and the generator; P k is the power flowing through load node k; LD k is the load of load node k in the mixed transaction of the integrated energy system; LD kB is the load of load node k in the bilateral trading network of the comprehensive energy system; G i G is the power generation of generator node i in the mixed transaction of the integrated energy system; iB is the power generation of generator node i in the bilateral trading network of the integrated energy system; is the upstream distribution matrix; The integrated energy system is split into r sub-integrated energy systems, and the network equivalents of the r sub-integrated energy systems are converted into a connected graph containing only tree branches; Remove the redundant interconnection line flows between different sub-integrated energy systems, so that there is only one interconnection line between different sub-integrated energy systems, and then an equivalent circulation is formed between different sub-integrated energy systems; The relationship F1 between each branch flow and each bilateral transaction in the integrated energy system is: F1=C cb B; Where: C cb =C a C b is the branch transaction association matrix, which is an l×(n+s)-order matrix, C a is the branch incidence matrix, C b is the transaction correlation matrix; l is the number of branches of the mixed transaction of the integrated energy system; B is the transaction volume of bilateral transactions; Calculate the network loss shared by each bilateral transaction and the total network loss that needs to be borne by the joint transaction; The network loss shared by each bilateral transaction is: Where: C cb，j C cb The (n×s)-order row vector of the j-th row of the matrix, j=1,...,l, where l is the number of branches of the hybrid transaction of the integrated energy system; P loss，j is the branch loss on the jth branch; B i is a (n×s)-order column vector whose ith element is the transaction amount of the ith bilateral transaction and whose other elements are 0. i =[0, ..., B i , ...0]; The total network loss that the joint venture transaction needs to bear is: The total network loss is evenly distributed to both sides of the virtual bilateral transaction, and the DC parameters of the generators and loads in the adjusted pooled transaction network are obtained; The DC parameters of generators and loads in the pooled trading network in hybrid trading are: Where: are respectively the active output of the balanced units, the active output of the unbalanced units and the active input of the load in the joint trading network; P G1，j , P G2，j , LD i They are respectively the active output of the balancing units, the active output of the unbalanced units and the active input of the loads in the mixed transaction of the integrated energy system; The carbon emission intensity of the units in the pool trading network; They are respectively the active output of the unbalanced units and the active input of the loads in the bilateral trading network; The total network loss is evenly distributed to both sides of the virtual bilateral transaction, and the DC parameters of the generator and load in the adjusted pooled transaction network are obtained as follows: Where: are the DC parameters of the generator and load in the adjusted pooling trading network; B s is the set of virtual bilateral transactions; Gi is the i-th generator set; To adjust the DC parameters of generators and loads in the pre-pool trading network respectively; Calculate the node active flux of the pooled trading network, and calculate the node carbon potential of each node based on the node active flux; The carbon flow distribution and load carbon flow rate of the joint trading network are calculated based on the node carbon potential.

2. A carbon flow calculation method for a comprehensive energy system joint trading network according to claim 1, characterized in that: Calculate the DC branch power flow distribution matrix, unit injection distribution matrix and load distribution matrix of the pool trading network; The node active flux is obtained based on the splicing matrix of the DC branch power flow distribution matrix and the unit injection distribution matrix. for: Where: ξ is a row vector of order l+h, and all elements are 1, l is the number of branches of the integrated energy system hybrid trading, and h is the number of generator sets of the integrated energy system hybrid trading; It is the concatenation matrix of the DC branch power flow distribution matrix and the unit injection distribution matrix. is the DC branch power flow distribution matrix, Inject distribution matrix into the unit; The node carbon potential of each node is obtained through the node active flux of the joint trading network, the DC branch power flow distribution matrix and the unit injection distribution matrix. for: Where: is the carbon emission intensity of the units in the pool trading network.

3. A carbon flow calculation method for a comprehensive energy system joint trading network according to claim 2, characterized in that: The carbon flow distribution and load carbon flow rate of the joint trading network are: Where: are the load carbon flow rate distribution matrix and the branch carbon flow rate distribution matrix of the joint trading network respectively; LD T，after is the load distribution matrix.

4. A carbon flow calculation device for a comprehensive energy system joint trading network, characterized in that: The device is applied to the method described in any one of claims 1 to 3, and the device comprises: A transaction conversion module is used to convert the joint transaction in the mixed transaction of the integrated energy system into multiple virtual bilateral transactions to obtain an integrated energy system containing only bilateral transactions; The network loss calculation module is used to calculate the network loss shared by each bilateral transaction and the total network loss that needs to be borne by the joint transaction; The DC parameter calculation module is used to evenly distribute the total network loss to both sides of the virtual bilateral transaction and obtain the DC parameters of the generators and loads in the adjusted pooled transaction network; The carbon flow distribution calculation module is used to calculate the node active flux of the joint trading network, and calculate the node carbon potential of each node based on the node active flux, and calculate the carbon flow distribution and load carbon flow rate of the joint trading network based on the node carbon potential.

5. A carbon flow calculation device for a comprehensive energy system joint trading network, characterized in that: including memory and processor; The memory is used to store computer program code and transmit the computer program code to the processor; The processor is configured to execute the method according to any one of claims 1 to 3 according to instructions in the computer program code.

6. A computer-readable storage medium, characterized in that: The computer-readable storage medium stores a computer program, and when the computer program is executed by a processor, the method according to any one of claims 1 to 3 is implemented.

Images