Carbon emission metering method and device for district heating system, electronic equipment and medium

Abstract

The application discloses a regional heating system carbon emission metering method and device, electronic equipment and medium, wherein the method comprises: obtaining the steady-state carbon emission of the current regional heating system by using a pre-trained steady-state carbon emission flow model; obtaining the dynamic carbon emission of the current regional heating system by using a pre-trained dynamic carbon emission flow model; and calculating the carbon emission of the current regional heating system according to the steady-state carbon emission and the dynamic carbon emission. Thus, the application can effectively identify the carbon emission details of the source, network and load of the regional heating system, clarify the carbon emission responsibilities on both sides of the source and load, realize the accurate metering of the carbon emission of the regional heating system, and has high application value.

Description

TECHNICAL FIELD

[0001] The present application relates to the field of low-carbon technology of energy systems, and particularly relates to a carbon emission metering method and device for a regional heating system, an electronic device and a medium. BACKGROUND

[0002] At present, especially in China, the traditional heating method mainly burns coal, and the heating of residents is a relative "hard constraint", which leads to the fact that the regional heating is one of the industries with the largest coal consumption and carbon emission in China. The above characteristics determine that the "carbon reduction" of the regional heating system is not only the task of the source side, but also needs the cooperation of the "source-grid-load" whole chain. Real-time, accurate and comprehensive carbon emission metering is the basis and premise for mastering the carbon emission status and trend of the regional heating system, tapping the carbon emission reduction potential, guiding user interaction to reduce carbon, and promoting the low-carbon transformation of the regional heating system.

[0003] In order to better combine the characteristics of energy systems with the concept of low-carbon development and expand the research of low-carbon energy technology, it is necessary to understand and analyze the carbon emission problem in the energy system, especially the regional heating system, from the "carbon perspective". In recent years, in the international import and export trade network, scholars have found that the life cycle carbon emissions of most imported and exported goods are concentrated in the production link, and no carbon emissions are generated in the use process, which will be unfair to countries and regions that mainly export. Based on this, there should be a coupling carbon emission transfer from the exporting country to the importing country in the international import and export trade logistics, so as to clarify the carbon emission responsibility that the importing country should bear due to the consumption of goods to obtain utility, thereby deriving the new concept of carbon emission flow. In the regional heating system, different supply technologies have different carbon emission characteristics, but the generated carbon emission flow is not different. In this way, compared with the application in trade logistics, the carbon emission flow has a more simple and flexible application space in the regional heating system. Therefore, the regional heating system has a good foundation for building a theoretical system of carbon emission flow.

[0004] In the current research on low-carbon technology of energy systems, macro statistical method and life cycle method are usually used. However, this calculation method based on the source side has a large time scale and a too macro system space range, can only give the total carbon emission, cannot effectively identify the carbon emission details of each link of the source, grid and load, leads to unreasonable carbon emission responsibility allocation on the source and load sides, and is difficult to guide the low-carbon optimal operation of the regional heating system. SUMMARY

[0005] The carbon emission metering method and device for a regional heating system, the electronic device and the medium provided by the present application can effectively identify the carbon emission details of each link of the source, grid and load of the regional heating system, clarify the carbon emission responsibility on the source and load sides, realize the accurate metering of the carbon emission of the regional heating system, and have high application value.

[0006] The first aspect of this application provides a method for measuring carbon emissions from a district heating system, comprising the following steps: obtaining the steady-state carbon emissions of the current district heating system using a pre-trained steady-state carbon emission flow model, wherein the steady-state carbon emission flow model is constructed based on the pipe carbon flow rate, network loss carbon emissions, node carbon flow density, pipe carbon flow density, heat source carbon flow rate, and heat load carbon flow rate of the current district heating system; obtaining the dynamic carbon emissions of the current district heating system using a pre-trained dynamic carbon emission flow model, wherein the dynamic carbon emission flow model is constructed based on the multi-time period water micro-element carbon flow rate, the actual outlet carbon flow rate of the pipes, the network loss carbon flow rate, and the node carbon flow density of the current district heating system; and calculating the carbon emissions of the current district heating system based on the steady-state carbon emissions and the dynamic carbon emissions.

[0007] Optionally, in one embodiment of this application, it further includes:

[0008] A steady-state carbon emission flow model for the district heating system is constructed based on the carbon flow rate of the pipelines, network loss carbon emissions, node carbon flow density, pipeline carbon flow density, heat source carbon flow rate, and heat load carbon flow rate. Specifically:

[0009] Determine the carbon flow rate of the district heating system piping, including the carbon flow rate of the supply water network and the carbon flow rate of the return water network:

[0010] Wherein, the carbon flow rate of the water supply network is:

[0011]

[0012]

[0013] in, and These represent the inlet carbon flow rate and outlet carbon flow rate of pipe k in the water supply network, respectively, in units of tCO2 / h; Let be the carbon flux density of pipe k in the water supply network, in units of tCO2 / MWh; and c be the specific heat of water, in units of MWh / (kg·℃). The mass flow rate of pipe k in the water supply network is expressed in kg / h. and These represent the inlet and outlet temperatures of pipe k in the water supply network, respectively, in °C; Ω BH For the collection of pipes in a district heating system;

[0014] The carbon flow rate of the recirculation network is:

[0015]

[0016]

[0017] in, and These represent the inlet carbon flow rate and outlet carbon flow rate of pipe k in the return water network, respectively, in units of tCO2 / h; Here is the carbon flux density in pipe k of the return water network, expressed in tCO2 / h. The mass flow rate of pipe k in the return water network is expressed in kg / h. and These represent the inlet and outlet temperatures of pipe k in the return water network, respectively, in °C.

[0018] Determine the carbon emissions from the district heating system network losses, including the carbon flow rate of the supply water network losses and the carbon flow rate of the return water network losses, specifically as follows:

[0019] Determine the temperature difference in the pipes of the water supply and return networks:

[0020]

[0021]

[0022] in, and Let represent the temperature difference between the two ends of pipe k in the supply and return water networks, respectively, in °C; determine the carbon flow rate of the supply network and the carbon flow rate of the return water network:

[0023]

[0024]

[0025] in, and These represent the carbon flow rate of pipe k in the water supply network and the return network, respectively, in units of tCO2 / h;

[0026] Determine the carbon flux density at nodes in the district heating system:

[0027] For each node in a district heating system, the laws of conservation of matter and energy are satisfied at the node:

[0028]

[0029]

[0030] in, The total mass flow rate through node n in the water supply network is expressed in kg / h. Ω represents the set of injection pipes for node n in the water supply network. NH A set of nodes for a district heating system; The water flow temperature at node n in the water supply network is expressed in °C.

[0031] For each node in a district heating system, where carbon emission conservation is simultaneously satisfied, the carbon flow rate of node n is equal to the sum of the outlet carbon flow rates of all injected pipes and the network loss carbon flow rate allocated to the injected pipes:

[0032]

[0033] in, is the carbon flow rate of node n in the water supply network, in units of tCO2 / h; X is the allocation coefficient of the carbon flow rate of the injected pipeline network loss;

[0034] Determine the carbon flux density at node n in the water supply network:

[0035]

[0036] in, Here is the carbon flux density at node n in the water supply network, expressed in tCO2 / MWh.

[0037] Determine the carbon flux density at node n in the reclaimed water network:

[0038]

[0039] in, Here is the carbon flux density at node n in the reclaimed water network, expressed in tCO2 / MWh. Let n be the set of outflow pipes of node n in the water supply network.

[0040] Determine the carbon flow density in the district heating system pipes:

[0041]

[0042]

[0043] in, and Let these represent the inlet and outlet nodes of pipe k in the water supply network, respectively.

[0044] Determine the carbon flow rate of the heat source in the district heating system:

[0045] Determine the heat output of the heat source:

[0046]

[0047] Among them, Q i The thermal output of heat source i is expressed in MW. Ω represents the node where heat source i is located; GH A collection of heat sources;

[0048] Carbon emission conservation is satisfied at the heat source node:

[0049]

[0050] in, Let be the carbon flux density of heat source i, in units of tCO2 / MWh;

[0051] Determine the carbon flow rate of the heat load in the district heating system:

[0052] Determine heat load requirements:

[0053]

[0054] Where, q j The heat demand is expressed in MW for heat load j. Ω represents the node where heat load j is located; LH For heat load aggregation;

[0055] Determine the carbon flux density at the heat load nodes:

[0056]

[0057] Determine the carbon flow rate under heat load:

[0058]

[0059] in, denoted as carbon flow rate for heat load j, in units of tCO2 / h.

[0060] Optionally, in one embodiment of this application, the method further includes: calculating the matrix representation of the steady-state carbon emission flow model of the district heating system.

[0061] Optionally, in one embodiment of this application, the step of calculating the matrix representation of the steady-state carbon emission flow model of the district heating system includes: constructing the branch heat flow matrix, the branch network loss matrix, and the node heat flow matrix of the district heating system; and after calculating the carbon flow density vector of the heating network nodes, calculating the branch carbon flow rate matrix, the network loss carbon flow rate matrix, and the load carbon flow rate vector of the water supply network and the return water network.

[0062] Optionally, in one embodiment of this application, the calculation of the matrix representation of the steady-state carbon emission flow model of the district heating system includes:

[0063] Construct the branch heat flux matrix of the district heating system, including the heat flux matrix of the supply water network branches and the heat flux matrix of the return water network branches:

[0064] Determine the elements of the heat flux matrix for the water supply network branches:

[0065]

[0066] in, and Heat flux matrix for water supply network branches Elements in;

[0067] Determine the elements of the heat flux matrix for the return water network branches:

[0068]

[0069] in, and Heat flux matrix of return water network branches Elements in;

[0070] Construct the branch network loss matrix of the district heating system, including the branch network loss matrix of the water supply network and the branch network loss matrix of the return water network:

[0071] Determine the elements of the water supply network branch loss matrix:

[0072]

[0073] in, and Represents the network loss matrix of the water supply network branches Elements in;

[0074] Determine the elements of the network loss matrix for the return water network branches:

[0075]

[0076] in, and Represents the network loss matrix of the return water network branches Elements in;

[0077] Construct the nodal heat flux matrix of the district heating system, including the heat flux matrix of the water supply network nodes and the heat flux matrix of the return water network nodes:

[0078] For nodes not connected to a heat source, determine the heat flux matrix of the water supply network nodes:

[0079]

[0080] in, This is the heat flux matrix for the water supply network nodes. N is the coefficient matrix of the branch heat flux; NH The number of nodes in a district heating system;

[0081] For nodes not connected to the heat source, determine the heat flux matrix of the return water network nodes:

[0082]

[0083] in, The heat flux matrix of the return water network nodes;

[0084] For nodes connected to the heat source, determine the node's overall energy flux matrix:

[0085]

[0086] in, The node-wide energy flux matrix; The coefficient matrix for injecting heat flow into the heat source; N GH The number of heat sources in a district heating system;

[0087] For all nodes, the total injected carbon emissions of a node are equal to the sum of the injected carbon flow rates of all branches connected to that node:

[0088]

[0089]

[0090] Where, ρ NHS The carbon flux density at node n in the water supply network The matrix formed; ρ NHR The carbon flux density at node n in the reclaimed water network The matrix formed;

[0091] For the heat load node, the carbon flux density is equal in both the supply and return water networks:

[0092] Bρ NHS =Bρ NHR

[0093] Where B is the heat load-node correlation matrix, and when heat load j is connected to node n, B jn =1, otherwise 0; for heat source nodes, according to the carbon emission conservation law, the matrix relationship is:

[0094]

[0095] Where C is a 0-1 matrix associated with the heat source node, and C is the matrix associated with the heat source node when node n is connected to the heat source. nn =1, otherwise 0;

[0096] The carbon flux density vector at the nodes of the heating network is:

[0097]

[0098] Calculate the branch carbon flow rate matrix, network loss carbon flow rate matrix, and load carbon flow rate vector for the water supply and return networks.

[0099] Optionally, in one embodiment of this application, it further includes: constructing a dynamic carbon emission flow model of the district heating system based on the multi-time period water micro-element carbon flow rate, the actual outlet carbon flow rate of the pipeline, the network loss carbon flow rate, and the node carbon flow density, specifically:

[0100] Determine the carbon flow rate of water micro-elements over multiple time periods:

[0101]

[0102] in, Let M be the carbon flux of the water element during time period t; σ represents the density of water; M k This represents the volume of pipe k; δ represents the injection period of the earliest water element component contained in the water flow exiting pipe k at time period t; k,t A represents the injection period of the latest water element component contained in the water flow out of pipe k at time period t; k,t Indicates the period before time t Total injected water flow rate for each time period; B k,t Representing the time period t-δ k,t Total injected water flow rate up to time period t; A k,t and B k,t The expression is:

[0103]

[0104]

[0105] Determine the actual outlet carbon flow rate and network loss carbon flow rate of the pipeline over multiple time periods:

[0106] Calculate the actual temperature of the water flow at the pipe outlet based on the transmission loss of the pipe:

[0107]

[0108] in, This represents the actual temperature of the water flowing out of the pipe. λ represents the weighted average temperature of the injected water flow over the previous period; L represents the thermal conductivity of the pipe; k Indicates the length of pipe k;

[0109] Determine the actual outlet carbon flow rate of the pipeline during time period t:

[0110]

[0111] Determine the network loss carbon flow rate of the pipeline during time period t:

[0112]

[0113] Determine the carbon flux density at multiple time points:

[0114]

[0115] For the return water network, determine the carbon flow density at multiple time points:

[0116]

[0117] A second aspect of this application provides a carbon emission metering device for a district heating system, comprising: multiple heat source-side carbon meters connected to a heat source for measuring carbon emissions from the heat source; multiple heat pipeline carbon meters connected to a heat pipeline, the heat source-side carbon meters and the heat pipeline carbon meters being connected together for measuring carbon emissions from the heat pipeline; multiple user-side carbon meters, one end of which is connected to a heat user and the other end of which is connected to the heat pipeline carbon meters, for obtaining the carbon potential of the node where the user is located and obtaining the carbon emissions caused by the heat consumption of each user; and a central server connected to the multiple heat source-side carbon meters and the multiple heat pipeline carbon meters, for calculating the carbon emission intensity of the heat source based on the coal consumption data of the heat source and the emission coefficient of the coal used, calculating the distribution of carbon emission flow in the heat pipeline, and calculating the carbon emissions corresponding to the pipelines connected to the nodes based on the data from the heat pipeline carbon meters.

[0118] Optionally, in one embodiment of this application, the carbon meter is also used to display in real time the carbon emission results of the district heating system calculated by the central server.

[0119] A third aspect of this application provides an electronic device, including: a memory, a processor, and a computer program stored in the memory and executable on the processor, wherein the processor executes the program to perform the carbon emission metering method for a district heating system as described in the above embodiments.

[0120] A fourth aspect of this application provides a computer-readable storage medium having a computer program stored thereon, which is executed by a processor to perform the carbon emission metering method for a district heating system as described in the above embodiments.

[0121] The embodiments of this application have the following beneficial effects:

[0122] 1) The carbon steady-state emission flow model of the district heating system in this application can identify the details of carbon emissions in each link of the district heating system, including pipeline carbon flow rate, network loss carbon emissions, node carbon flow density, pipeline carbon flow density, heat source carbon flow rate and heat load carbon flow rate, thereby reasonably and effectively distributing the "carbon reduction" task on the source side to the entire "source-network-load" chain, and further clarifying the carbon emission responsibilities of each link of the district heating system.

[0123] 2) This application also considers the dynamic transmission characteristics of the district heating system and constructs a dynamic carbon emission flow model for the district heating system, which can grasp the current status and trend of carbon emissions of the system in real time, accurately and comprehensively.

[0124] 3) The carbon meter system of this application, with the help of a data acquisition, processing and display system, can intuitively display the calculation results of the carbon meter system to the heat users, which is conducive to further exploring the carbon emission reduction potential of the system and guiding users to interact in carbon reduction, and has high application value.

[0125] Additional aspects and advantages of this application will be set forth in part in the description which follows, and in part will be obvious from the description, or may be learned by practice of this application. Attached Figure Description

[0126] The above and / or additional aspects and advantages of this application will become apparent and readily understood from the following description of the embodiments taken in conjunction with the accompanying drawings, wherein:

[0127] Figure 1 This is a flowchart of a carbon emission metering method for a district heating system according to an embodiment of this application;

[0128] Figure 2 This is a schematic diagram of the framework structure of a carbon emission metering method for a district heating system according to an embodiment of this application;

[0129] Figure 3 This is a schematic diagram illustrating the solution results of steady-state carbon emission flow in a district heating system according to an embodiment of this application.

[0130] Figure 4 This is a schematic diagram of the solution results for the dynamic carbon emission flow of a district heating system according to an embodiment of this application;

[0131] Figure 5 This is an example diagram of a carbon emission metering device for a district heating system according to an embodiment of this application;

[0132] Figure 6 A schematic diagram of the structure of the electronic device provided in the application embodiment. Detailed Implementation

[0133] The embodiments of this application are described in detail below. Examples of these embodiments are shown in the accompanying drawings, wherein the same or similar reference numerals denote the same or similar elements or elements having the same or similar functions throughout. The embodiments described below with reference to the accompanying drawings are exemplary and intended to explain this application, and should not be construed as limiting this application.

[0134] Figure 1 This is a flowchart of a carbon emission metering method for a district heating system according to an embodiment of this application.

[0135] like Figure 1 As shown, the carbon emission metering method for this district heating system includes the following steps:

[0136] In step S101, the steady-state carbon emission flow model of the current district heating system is obtained using a pre-trained steady-state carbon emission flow model. The steady-state carbon emission flow model is constructed based on the carbon flow rate of the pipeline, the carbon emission of network loss, the carbon flow density of the nodes, the carbon flow density of the pipeline, the carbon flow rate of the heat source, and the carbon flow rate of the heat load of the current district heating system.

[0137] It is understood that the embodiments of this application establish a steady-state carbon emission flow model for a district heating system based on carbon emission flow theory, and use the constructed steady-state carbon emission flow model for a district heating system to obtain the current steady-state carbon emissions of the district heating system.

[0138] Specifically, the steady-state carbon emission flow model for a district heating system consists of the carbon flow rate of the system's pipelines, network loss carbon emissions, node carbon flow density, pipeline carbon flow density, heat source carbon flow rate, and heat load carbon flow rate. The specific steps are as follows:

[0139] 1-1) Under steady-state conditions, since transmission delay is not considered, there is no convergence or divergence of water flow on a single pipe. Therefore, the carbon flux density at the pipe inlet and outlet is equal and is uniformly defined as the pipe's carbon flux density. The inlet carbon flux rate and outlet carbon flux rate are equal to the products of the inlet energy flow and outlet energy flow, respectively, and the pipe's carbon flux density. The specific steps for determining the carbon flux rate of the district heating system pipes, including the carbon flux rate of the supply water network and the carbon flux rate of the return water network, are as follows:

[0140] 1-1-1) Determine the carbon flow rate of the water supply network, as shown in the following expression:

[0141]

[0142] in, and These represent the inlet carbon flow rate and outlet carbon flow rate of pipe k in the water supply network, respectively, in units of tCO2 / h; Let be the carbon flux density of pipe k in the water supply network, in units of tCO2 / MWh; and c be the specific heat of water, in units of MWh / (kg·℃). The mass flow rate of pipe k in the water supply network is expressed in kg / h. and These represent the inlet and outlet temperatures of pipe k in the water supply network, respectively, in °C; Ω BH For district heating system pipe collection.

[0143] 1-1-2) Determine the carbon flow rate of the return water network, as follows:

[0144]

[0145] in, and These represent the inlet carbon flow rate and outlet carbon flow rate of pipe k in the return water network, respectively, in units of tCO2 / h; Here is the carbon flux density in pipe k of the return water network, expressed in tCO2 / h. The mass flow rate of pipe k in the return water network is expressed in kg / h. and These represent the inlet and outlet temperatures of pipe k in the return water network, respectively, in °C.

[0146] 1-2) Significant transmission losses exist in the heating network, resulting in additional carbon emissions. The carbon emissions from network losses in the district heating system need to be determined, including the carbon flow rate of the supply water network and the carbon flow rate of the return water network. The specific steps are as follows:

[0147] 1-2-1) The losses in the heating network are mainly caused by the temperature loss of the water flowing in the pipes. Therefore, the carbon emissions from network losses mainly depend on the temperature difference between the two ends of the pipes. The temperature difference between the pipes in the supply and return water networks is determined as follows:

[0148]

[0149] in, and These represent the temperature difference between the two ends of pipe k in the water supply network and the return network, respectively, in °C.

[0150] 1-2-2) Determine the carbon flow rate of the water supply network and the carbon flow rate of the return water network, as follows:

[0151]

[0152] in, and These represent the carbon flow rate of pipe k in the water supply network and the return network, respectively, in units of tCO2 / h.

[0153] 1-3) For each node in the heating network, there may be a convergence of water flows from multiple pipes, where the conservation of mass and energy is satisfied. Determine the carbon flux density at each node of the district heating system; the specific steps are as follows:

[0154] 1-3-1) For each node in a district heating system, the laws of conservation of matter and energy are satisfied at the node, as follows:

[0155]

[0156] in, The total mass flow rate through node n in the water supply network is expressed in kg / h. Ω represents the set of injection pipes for node n in the water supply network. NH A set of nodes for a district heating system; The water temperature at node n in the water supply network is expressed in °C.

[0157] 1-3-2) For each node in the district heating system, carbon emission conservation is simultaneously satisfied at the node. Therefore, the carbon flow rate of node n is equal to the sum of the outlet carbon flow rates of all injection pipes and the network loss carbon flow rate allocated to the injection pipes, as follows:

[0158]

[0159] in, X represents the carbon flow rate of node n in the water supply network, expressed in tCO2 / h; X is the allocation coefficient for the carbon flow rate of the injected pipeline network loss.

[0160] 1-3-3) Determine the carbon flux density at node n in the water supply network, as follows:

[0161]

[0162] in, Let be the carbon flux density of node n in the water supply network, expressed in tCO2 / MWh.

[0163] 1-3-4) Determine the carbon flux density at node n in the reclaimed water network, as follows:

[0164]

[0165] in, Here is the carbon flux density at node n in the reclaimed water network, expressed in tCO2 / MWh. Let n be the set of outflow pipes of node n in the water supply network.

[0166] 1-4) Each pipe in the water supply and return networks satisfies the energy distribution criterion, that is, the carbon flow density of the pipe is equal to the carbon flow density of the pipe's injection node. The carbon flow density of the pipes in the district heating system is determined as follows:

[0167]

[0168] in, and These represent the inlet and outlet nodes of pipe k in the water supply network, respectively.

[0169] 1-5) The heat source and heat load are the connection points between the supply and return water networks. Through these two points, the carbon emission flows from the supply and return water networks can be linked, thus forming a complete heating network carbon emission flow. Determine the carbon flow rate of the heat source in the district heating system; the specific steps are as follows:

[0170]

[0171] Among them, Q i The thermal output of heat source i is expressed in MW. Ω represents the node where heat source i is located; GH It is a collection of heat sources.

[0172] 1-5-2) At the heat source node, carbon emission conservation is also satisfied, meaning the carbon flow rate of the supply water at the heat source node is equal to the sum of the carbon flow rate injected by the heat source and the carbon flow rate of the return water at the heat source node. This is equivalent to the heat source injecting corresponding carbon emissions while simultaneously injecting energy into the heating network; the ratio between the two is the carbon flux density of the heat source. Specifically:

[0173]

[0174] in, Let be the carbon flux density of heat source i, in units of tCO2 / MWh.

[0175] 1-6) The carbon flow rate of the heat load depends on the load magnitude and the carbon flow density at the load nodes. The heat load magnitude can be obtained from the mass flow rate at the load nodes and the supply and return water temperatures at the nodes; the specific steps are as follows:

[0176] 1-6-1) Determine the heat load requirement, as follows:

[0177]

[0178] Where, q j The heat demand is expressed in MW for heat load j. Ω represents the node where heat load j is located; LH For heat load aggregation;

[0179] 1-6-2) There is no external energy injection at the heat load node; therefore, the carbon flux density at the node connected to the heat load is equal in both the supply and return water networks. The carbon flux density at the heat load node is determined as follows:

[0180]

[0181] 1-6-3) This allows for the further determination of the carbon flow rate under heat load, as detailed below:

[0182]

[0183] in, denoted as carbon flow rate for heat load j, in units of tCO2 / h.

[0184] Optionally, in one embodiment of this application, the carbon emission metering method for a district heating system further includes: calculating the matrix representation of the steady-state carbon emission flow model of the district heating system.

[0185] In one embodiment of this application, the matrix representation of the steady-state carbon emission flow model of the district heating system is calculated, including: constructing the branch heat flow matrix, the branch network loss matrix, and the node heat flow matrix of the district heating system; after calculating the carbon flow density vector of the heating network nodes, calculating the branch carbon flow rate matrix, the network loss carbon flow rate matrix, and the load carbon flow rate vector of the water supply network and the return water network.

[0186] Specifically, the matrix expression for the steady-state carbon emission flow model of the district heating system is established, including constructing the branch heat flow matrix, branch network loss matrix, and nodal heat flux matrix of the district heating system; the specific steps are as follows:

[0187] 2-1) The elements of the branch energy flow matrix are related to the water flow rate and temperature in the pipes. Construct the branch heat flow matrix of the district heating system, including the heat flow matrix of the supply network branches and the heat flow matrix of the return network branches; the specific steps are as follows:

[0188] 2-1-1) Determine the elements of the heat flux matrix for the water supply network branches, as shown in the following expression:

[0189]

[0190] in, and Heat flux matrix for water supply network branches Elements in;

[0191] 2-1-2) Determine the elements of the heat flux matrix for the return water network branches, as shown in the following expression:

[0192]

[0193] in, and Heat flux matrix of return water network branches Elements in;

[0194] 2-2) The elements of the branch network loss matrix are related to the water flow rate and temperature in the pipeline. Construct the branch network loss matrix of the district heating system, including the water supply network branch network loss matrix and the return water network branch network loss matrix; the specific steps are as follows:

[0195] 2-2-1) Determine the elements of the water supply network branch loss matrix, as shown in the following expression:

[0196]

[0197] in, and Represents the network loss matrix of the water supply network branches Elements in;

[0198] 2-2-2) Determine the elements of the network loss matrix for the return water network branches, as shown in the following expression:

[0199]

[0200] in, and Represents the network loss matrix of the return water network branches. Elements in;

[0201] 2-3) The nodal energy flux matrix of the heating network is related to the branch energy flow matrix and the heat source injection matrix. Construct the nodal heat flux matrix of the district heating system, including the heat flux matrix of the supply network nodes and the heat flux matrix of the return water network nodes; the specific steps are as follows:

[0202] 2-3-1) For nodes not connected to a heat source, the energy flow of the node is equal to the sum of the energy flows of each injection branch. The heat flux matrix of the water supply network nodes is determined as follows:

[0203]

[0204] in, This is the heat flux matrix for the water supply network nodes; N is the coefficient matrix of the branch heat flux; NH The number of nodes in a district heating system;

[0205] 2-3-2) For nodes not connected to a heat source, the energy flow of the node is equal to the sum of the energy flows of each injection branch. The heat flux matrix of the return water network nodes is determined as follows:

[0206]

[0207] in, The heat flux matrix of the return water network nodes;

[0208] 2-3-3) For nodes connected to a heat source, the energy flow through the node is also related to the injection of the heat source. The comprehensive energy flux matrix of the node is determined as follows:

[0209]

[0210] in, The node-wide energy flux matrix; The coefficient matrix for injecting heat flow into the heat source; N GH The number of heat sources in a district heating system;

[0211] 2-4) For all nodes, the total injected carbon emissions of a node are equal to the sum of the injected carbon flow rates of all branches connected to the node, as follows:

[0212]

[0213] Where, ρ NHS The carbon flux density at node n in the water supply network (step 1-3-3) The matrix formed; ρ NHR The carbon flux density at node n in the reclaimed water network (steps 1-3-4) The matrix formed;

[0214] 2-4-1) For the heat load node, the carbon flux density is equal in both the supply and return water networks, as expressed below:

[0215] Bρ NHS =Bρ NHR (twenty three)

[0216] Where B is the heat load-node correlation matrix, and when heat load j is connected to node n, B jn =1, otherwise 0;

[0217] 2-4-2) For the heat source node, according to the carbon emission conservation relationship in step 1-5-2), the matrix relationship can be obtained as follows:

[0218]

[0219] Where C is a 0-1 matrix associated with the heat source node, and C is the matrix associated with the heat source node when node n is connected to the heat source. nn =1, otherwise 0;

[0220] 2-5) Further simplification of step 2-4) yields the following matrix expression:

[0221]

[0222] Therefore, the carbon flux density vector at the nodes of the heating network can be obtained by solving for it, and the expression is:

[0223]

[0224] Therefore, after obtaining the nodal carbon flow density vector of the heating network, the branch carbon flow rate matrix, network loss carbon flow rate matrix, and load carbon flow rate vector of the water supply network and return network can also be calculated.

[0225] In step S102, the dynamic carbon emissions of the current district heating system are obtained using a pre-trained dynamic carbon emission flow model. The dynamic carbon emission flow model is constructed based on the multi-period water micro-element carbon flow rate, the actual outlet carbon flow rate of the pipeline, the network loss carbon flow rate, and the node carbon flow density of the current district heating system.

[0226] Due to transmission delays in heating networks, the water flow exiting the network pipes at a given time interval differs from the water flow entering the pipes, resulting in differences between the inlet and outlet carbon flow rates and densities. A unit characterization method is used to describe the dynamic transmission characteristics of the district heating system. The water flow exiting the pipes is actually a mixture of the inlet water flow from several previous time intervals, which can be considered as a combination of several water micro-element components. A dynamic carbon emission flow model of the district heating system is established, including the carbon flow rates of water micro-element components at multiple time intervals, the actual outlet carbon flow rate of the pipes, the network loss carbon flow rate, and the nodal carbon flow density. The specific steps are as follows:

[0227] 3-1) A water element is specifically defined as the mass flow rate of water in a pipe per unit time. According to the conservation of carbon emissions, the carbon flow rate corresponding to the mixed hot water can also be expressed as a linear combination of the carbon flow rates corresponding to several water elements, as shown in the following expression:

[0228]

[0229] in, Let M be the carbon flux of the water element during time period t; σ represents the density of water; M k This represents the volume of pipe k; δ represents the injection period of the earliest water element component contained in the water flow exiting pipe k at time period t; k,t This represents the injection time of the latest water element component contained in the water flow out of pipe k during time period t, and is related to the transmission delay of the district heating system; A k,t Indicates the period before time t Total injected water flow rate for each time period; B k,t Representing the time period t-δ k,t Total injected water flow rate up to time period t; A k,t and B k,t The expression is as follows:

[0230]

[0231]

[0232] 3-2) Determine the actual outlet carbon flow rate and network loss carbon flow rate of the pipeline in multiple time periods; the specific steps are as follows:

[0233] 3-2-1) Considering the transmission loss of the pipeline, the actual temperature of the water flow at the pipeline outlet can be obtained, as expressed below:

[0234]

[0235] in, This represents the actual temperature of the water flowing out of the pipe. λ represents the weighted average temperature of the injected water flow over the previous period; L represents the thermal conductivity of the pipe; k Indicates the length of pipe k;

[0236] 3-2-2) Only the water temperature decreases along the pipe; there is no mixing or outflow process. Therefore, the carbon flow density remains constant during this process. The actual outlet carbon flow rate of the pipe in time period t is expressed as follows:

[0237]

[0238] 3-2-3) Determine the network loss carbon flow rate of the pipeline during time period t, as shown in the following expression:

[0239]

[0240] 3-3) Based on steps 1-3-2) and 1-3-3), determine the carbon flux density at multiple time periods, as follows:

[0241]

[0242] 3-4) For the return water network, determine the carbon flow density at multiple time points, as follows:

[0243]

[0244] In step S103, the carbon emissions of the current district heating system are calculated based on the steady-state carbon emissions and the dynamic carbon emissions.

[0245] like Figure 2 The diagram illustrates the framework structure of a carbon emission metering method for a district heating system according to an embodiment of this application. Figure 2 Based on the steady-state carbon emission flow model and dynamic carbon emission flow model constructed by the method shown, the steady-state and dynamic carbon emission indicators of the current district heating system are obtained.

[0246] In one specific embodiment of this application, based on a 6-node district heating system, the network structure and its steady-state carbon emission flow calculation results are as follows: Figure 3As shown. The heat loads of nodes 5 and 6 are 120MW and 100MW, respectively. Based on the steady-state carbon emission flow calculations, the carbon flux densities of the two heat sources HS1 and HS2 are 0.222tCO2 / MWh and 0.267tCO2 / MWh, respectively. The total carbon flux rates of the heat source injection, the total carbon flux rates of the heat load, and the total carbon flux rates caused by network losses are 68.63tCO2 / h, 53.14tCO2 / h, and 15.49tCO2 / h, respectively.

[0247] Based on the dynamic carbon emission flow model of the district heating system proposed in this application, the carbon flow density of each node in the district heating system's water supply network at different time periods is calculated, such as... Figure 4 As shown.

[0248] The carbon emission metering method for district heating systems proposed in this application establishes a steady-state carbon emission flow model based on carbon emission flow theory and provides corresponding matrix calculation expressions. Considering heat transfer delay, a dynamic carbon emission flow model for the district heating system is established, and the steady-state and dynamic carbon emission indicators of the district heating system are obtained by solving the model. This effectively identifies the carbon emission details at each stage of the district heating system—source, network, and load—clarifies the carbon emission responsibilities on both the source and load sides, and achieves accurate measurement of carbon emissions from the district heating system. Furthermore, with the help of a data acquisition, processing, and display system, the calculation results of the carbon meter system can be intuitively presented to heating users, guiding them to interact and reduce carbon emissions, thus possessing high application value.

[0249] Next, the carbon emission metering device for a district heating system proposed according to an embodiment of this application is described with reference to the accompanying drawings.

[0250] Figure 5 This is an example diagram of a carbon emission metering device for a district heating system according to an embodiment of this application.

[0251] like Figure 5 As shown, the carbon emission metering device for the regional heating system includes: multiple carbon meters on the heat source side, multiple carbon meters on the heat pipeline side, multiple carbon meters on the user side, and a central server.

[0252] Among them, the carbon meter on the heat source side is connected to the heat source and is used to measure carbon emissions on the heat source side.

[0253] The carbon meter for the heating pipeline is connected to the heating pipeline, and the carbon meter on the heat source side is connected to the carbon meter for the heating pipeline to measure the carbon emissions of the heating pipeline.

[0254] One end of the user-side carbon meter is connected to the heat user's end, and the other end is connected to the carbon meter in the heat pipeline. It is used to obtain the carbon potential of the node where the user is located and to obtain the carbon emissions caused by each user's heat consumption.

[0255] The central server, acting as the computing center, is responsible for calculating the carbon emission flow of the system and communicating with carbon meters distributed in various stages to obtain the carbon emission information of each stage. The central server is connected to multiple carbon meters on the heat source side and multiple carbon meters on the heat pipeline side. It is used to calculate the carbon emission intensity of the heat source based on the coal consumption data of the heat source and the emission coefficient of the coal used, and to calculate the distribution of carbon emission flow in the heat pipeline. Based on the data of the carbon meters on the heat pipeline, it calculates the carbon emission amount corresponding to the pipeline connected to the node.

[0256] Optionally, in one embodiment of this application, the carbon meter is also used to display in real time the carbon emission results of the district heating system calculated by the central server.

[0257] In the embodiments of this application, carbon meters are used to collect and calculate the information required for carbon emission measurement, and transmit the relevant information to a central server. The central server uses the collected information to calculate the carbon emissions of the district heating system, feeds back the calculation results to each carbon meter, and then each carbon meter displays the calculation results fed back by the central server in real time.

[0258] In the embodiments of this application, a carbon meter system is composed of multiple carbon meters. The carbon meter system is a hierarchical system, consisting of an upper layer and a lower layer. The upper layer system includes a central server, heat source-side carbon meters, and heat pipe carbon meters. It is used to calculate the carbon emission intensity of the heat source based on the coal consumption data of the heat source and the emission coefficient of the coal used, and then calculate the distribution of carbon emission flow in the heat pipe. Based on the data of the heat pipe carbon meters, it calculates the carbon emission amount corresponding to the pipe connected to the node. The lower layer system consists of user-side carbon meters. The user-side carbon meters communicate with the heat pipe carbon meters to obtain the carbon potential of the node where the user is located and to obtain the carbon emission amount caused by the heat consumption of each user.

[0259] It should be noted that the foregoing explanation of the embodiment of the carbon emission metering method for district heating systems also applies to the carbon emission metering device for district heating systems in this embodiment, and will not be repeated here.

[0260] The carbon emission metering device for a district heating system proposed in this application consists of carbon meters distributed throughout the district heating system for measuring carbon emissions, a central server, and communication lines. Specifically, it is divided into carbon meters on the heat source side, carbon meters on the pipeline side, and carbon meters on the heat user side. This effectively identifies the details of carbon emissions at each stage of the district heating system (source, network, and load), clarifies the carbon emission responsibilities on both the source and load sides, and achieves accurate measurement of carbon emissions in the district heating system, thus possessing high application value.

[0261] Figure 6 A schematic diagram of the structure of an electronic device provided in an embodiment of this application. The electronic device may include:

[0262] The memory 601, the processor 602, and the computer program stored on the memory 601 and capable of running on the processor 602.

[0263] When the processor 602 executes the program, it implements the carbon emission metering method for the district heating system provided in the above embodiments.

[0264] Furthermore, electronic devices also include:

[0265] Communication interface 603 is used for communication between memory 601 and processor 602.

[0266] The memory 601 is used to store computer programs that can run on the processor 602.

[0267] The memory 601 may include high-speed RAM memory, and may also include non-volatile memory, such as at least one disk storage device.

[0268] If the memory 601, processor 602, and communication interface 603 are implemented independently, then the communication interface 603, memory 601, and processor 602 can be interconnected via a bus to complete communication between them. The bus can be an Industry Standard Architecture (ISA) bus, a Peripheral Component Interconnect (PCI) bus, or an Extended Industry Standard Architecture (EISA) bus, etc. Buses can be categorized as address buses, data buses, control buses, etc. For ease of representation, Figure 6 The bus is represented by a single thick line, but this does not mean that there is only one bus or one type of bus.

[0269] Optionally, in a specific implementation, if the memory 601, processor 602, and communication interface 603 are integrated on a single chip, then the memory 601, processor 602, and communication interface 603 can communicate with each other through an internal interface.

[0270] The processor 602 may be a central processing unit (CPU), an application specific integrated circuit (ASIC), or one or more integrated circuits configured to implement the embodiments of this application.

[0271] This embodiment also provides a computer-readable storage medium storing a computer program, characterized in that the program, when executed by a processor, implements the above-described method for measuring carbon emissions in a district heating system.

[0272] In the description of this specification, the references to terms such as "one embodiment," "some embodiments," "example," "specific example," or "some examples," etc., indicate that a specific feature, structure, material, or characteristic described in connection with that embodiment or example is included in at least one embodiment or example of this application. In this specification, the illustrative expressions of the above terms do not necessarily refer to the same embodiment or example. Furthermore, the specific features, structures, materials, or characteristics described may be combined in any suitable manner in one or more embodiments or examples. Moreover, without contradiction, those skilled in the art can combine and integrate the different embodiments or examples described in this specification, as well as the features of different embodiments or examples.

[0273] Furthermore, the terms "first" and "second" are used for descriptive purposes only and should not be construed as indicating or implying relative importance or implicitly specifying the number of technical features indicated. Thus, a feature defined as "first" or "second" may explicitly or implicitly include at least one of that feature. In the description of this application, "N" means at least two, such as two, three, etc., unless otherwise explicitly specified.

[0274] Any process or method described in the flowchart or otherwise herein can be understood as representing a module, segment, or portion of code comprising one or more N executable instructions for implementing custom logic functions or processes, and the scope of the preferred embodiments of this application includes additional implementations in which functions may be performed not in the order shown or discussed, including substantially simultaneously or in reverse order depending on the functions involved, as should be understood by those skilled in the art to which embodiments of this application pertain.

[0275] It should be understood that the various parts of this application can be implemented using hardware, software, firmware, or a combination thereof. In the above embodiments, the N steps or methods can be implemented using software or firmware stored in memory and executed by a suitable instruction execution system. For example, if implemented in hardware as in another embodiment, it can be implemented using any one or a combination of the following techniques known in the art: discrete logic circuits having logic gates for implementing logical functions on data signals, application-specific integrated circuits (ASICs) having suitable combinational logic gates, programmable gate arrays (PGAs), field-programmable gate arrays (FPGAs), etc.

[0276] Those skilled in the art will understand that all or part of the steps of the methods described in the above embodiments can be implemented by a program instructing related hardware. The program can be stored in a computer-readable storage medium, and when executed, it includes one or a combination of the steps of the method embodiments.

Claims

1. A method for measuring carbon emissions in a district heating system, characterized in that, Includes the following steps: The steady-state carbon emissions of the current district heating system are obtained by using a pre-trained steady-state carbon emission flow model, wherein the steady-state carbon emission flow model is constructed based on the pipe carbon flow rate, network loss carbon emissions, node carbon flow density, pipe carbon flow density, heat source carbon flow rate and heat load carbon flow rate of the current district heating system. The dynamic carbon emissions of the current district heating system are obtained using a pre-trained dynamic carbon emission flow model, wherein the dynamic carbon emission flow model is constructed based on the multi-time period water element carbon flow rate, the actual outlet carbon flow rate of the pipeline, the network loss carbon flow rate, and the node carbon flow density of the current district heating system; and The carbon emissions of the current district heating system are calculated based on the steady-state carbon emissions and the dynamic carbon emissions. Based on the multi-time period water micro-element carbon flow rate, the actual outlet carbon flow rate of the pipeline, the network loss carbon flow rate, and the node carbon flow density, a dynamic carbon emission flow model for the district heating system is constructed, specifically as follows: Determine the carbon flow rate of water micro-elements over multiple time periods: in, For time period t Carbon flux of water microelement; Indicates the density of water; Indicates pipeline Volume; Indicates time period t Outflow pipe The earliest water element injection period contained in the water flow; Indicates time period t Outflow pipe The water flow includes the injection period of the latest water element component; Indicates time period t Previous Total water injection flow rate for each time period; Indicates time period When the time comes Total injected water flow rate; and The expression is: Determine the actual outlet carbon flow rate and network loss carbon flow rate of the pipeline over multiple time periods: Calculate the actual temperature of the water flow at the pipe outlet based on the transmission loss of the pipe: in, This represents the actual temperature of the water flowing out of the pipe. This represents the weighted average temperature of the injected water flow over the previous period; Indicates the thermal conductivity of the pipe; Indicates pipeline Length; Determine the time period t Actual outlet carbon flow rate of the pipeline: Determine the time period t Pipeline network loss carbon flow rate: Determine the carbon flux density at multiple time points: For the return water network, determine the carbon flow density at multiple time points: ； in, For water supply pipe k Water flow volume over a given time period; For water supply pipe k Water flow volume over a given time period; Let be the water flow volume in the return water pipe k during time period t.

2. The method according to claim 1, characterized in that, Also includes: A steady-state carbon emission flow model for the district heating system is constructed based on the carbon flow rate of the pipelines, network loss carbon emissions, node carbon flow density, pipeline carbon flow density, heat source carbon flow rate, and heat load carbon flow rate. Specifically: Determine the carbon flow rate of the district heating system piping, including the carbon flow rate of the supply water network and the carbon flow rate of the return water network: Wherein, the carbon flow rate of the water supply network is: in, and These represent the pipes in the water supply network. The inlet carbon flow rate and outlet carbon flow rate are expressed in tCO2 / h. Pipes in the water supply network The carbon flux density is expressed in tCO2 / MWh. Specific heat of water, expressed in MWh / (kg) ); Pipes in the water supply network The mass flow rate is expressed in kg / h. and These represent the pipes in the water supply network. The inlet and outlet temperatures, in units of ; For the collection of pipes in a district heating system; The carbon flow rate of the recirculation network is: in, and These represent the pipes in the return water network. The inlet carbon flow rate and outlet carbon flow rate are expressed in tCO2 / h. Pipes in the return water network The carbon flux density, expressed in tCO2 / h; Pipes in the return water network The mass flow rate is expressed in kg / h. and These represent the pipes in the return water network. The inlet and outlet temperatures, in units of ; Determine the carbon emissions from the district heating system network losses, including the carbon flow rate of the supply water network losses and the carbon flow rate of the return water network losses, specifically as follows: Determine the temperature difference in the pipes of the water supply and return networks: in, and These represent the pipes in the water supply network and the return network, respectively. The temperature difference between the two ends, in units of ; Determine the carbon flow rate of the supply network and the carbon flow rate of the return network: in, and These represent the pipes in the water supply network and the return network, respectively. The carbon flow rate of the network loss is expressed in tCO2 / h. Determine the carbon flux density at nodes in the district heating system: For each node in a district heating system, the laws of conservation of matter and energy are satisfied at the node: in, For nodes in the water supply network Total mass flow rate, expressed in kg / h; Nodes in the water supply network The set of injection pipelines; A set of nodes for a district heating system; Nodes in the water supply network Water flow temperature, in units of ; For each node in a district heating system, carbon emission conservation is simultaneously satisfied at the node. The carbon flow rate is equal to the sum of the outlet carbon flow rates of all injection pipes and the network loss carbon flow rate allocated to the injection pipes: in, Nodes in the water supply network The carbon flow rate, expressed in tCO2 / h; The allocation factor for the carbon flow rate of the injected pipeline network loss; Determine the nodes in the water supply network Carbon flux density: in, Nodes in the water supply network The carbon flux density is expressed in tCO2 / MWh. Determine the nodes in the backwater network Carbon flux density: in, For nodes in the water return network The carbon flux density is expressed in tCO2 / MWh. Nodes in the water supply network A collection of outflow pipes; Determine the carbon flow density in the district heating system pipes: in, and These represent the pipes in the water supply network. Injection nodes and outflow nodes; Determine the carbon flow rate of the heat source in the district heating system: Determine the heat output of the heat source: in, heat source Thermal output, measured in MW; Indicates heat source The node it is located at; A collection of heat sources; Carbon emission conservation is satisfied at the heat source node: in, heat source The carbon flux density is expressed in tCO2 / MWh. Determine the carbon flow rate of the heat load in the district heating system: Determine heat load requirements: in, For heat load The heat demand, in MW; Indicates heat load The node it is located at; For heat load aggregation; Determine the carbon flux density at the heat load nodes: Determine the carbon flow rate under heat load: in, For heat load The carbon flow rate is expressed in tCO2 / h.

3. The method according to claim 2, characterized in that, Also includes: Calculate the matrix representation of the steady-state carbon emission flow model of the district heating system.

4. The method according to claim 3, characterized in that, The matrix representation of the steady-state carbon emission flow model of the district heating system includes: Construct the branch heat flux matrix, branch network loss matrix, and node heat flux matrix of the district heating system; After calculating the carbon flow density vector of the heating network nodes, calculate the branch carbon flow rate matrix, network loss carbon flow rate matrix, and load carbon flow rate vector of the water supply network and return network.

5. The method according to claim 4, characterized in that, The matrix representation of the steady-state carbon emission flow model of the district heating system includes: Construct the branch heat flux matrix of the district heating system, including the heat flux matrix of the supply water network branches and the heat flux matrix of the return water network branches: Determine the elements of the heat flux matrix for the water supply network branches: in, and Heat flux matrix for water supply network branches Elements in; Determine the elements of the heat flux matrix for the return water network branches: in, and For the heat flux matrix of the return water network branches Elements in; Construct the branch network loss matrix of the district heating system, including the branch network loss matrix of the water supply network and the branch network loss matrix of the return water network: Determine the elements of the water supply network branch loss matrix: in, and Represents the network loss matrix of the water supply network branches Elements in; Determine the elements of the network loss matrix for the return water network branches: in, and Represents the network loss matrix of the return water network branches. Elements in; Construct the nodal heat flux matrix of the district heating system, including the heat flux matrix of the water supply network nodes and the heat flux matrix of the return water network nodes: For nodes not connected to a heat source, determine the heat flux matrix of the water supply network nodes: in, This is the heat flux matrix for the water supply network nodes; This is the coefficient matrix of the branch heat flux; The number of nodes in a district heating system; This is the energy flow distribution matrix of the water supply network branches; For nodes not connected to the heat source, determine the heat flux matrix of the return water network nodes: in, The heat flux matrix of the return water network nodes; The energy flow distribution matrix of the return water network branches; For nodes connected to the heat source, determine the node's overall energy flux matrix: in, The node-wide energy flux matrix; The coefficient matrix for injecting heat flow into the heat source; The number of heat sources in a district heating system; The distribution matrix of the heat flow injected into the heat source; For all nodes, the total injected carbon emissions of a node are equal to the sum of the injected carbon flow rates of all branches connected to that node: in, For the nodes in the water supply network carbon flux density The matrix formed; For the nodes in the water return network carbon flux density The matrix formed; For the heat load node, the carbon flux density is equal in both the supply and return water networks: in, The heat load-node correlation matrix is ​​given when the heat load... With nodes When connected, Otherwise, it is 0; For the heat source node, according to the carbon emission conservation law, the matrix relationship is as follows: in, The transpose of the distribution matrix of the heat flow injected into the heat source; The carbon flux density vector of the heat source; This is a 0-1 matrix associated with the heat source nodes. If the nodes... When connected to a heat source, Otherwise, it is 0; The carbon flux density vector at the nodes of the heating network is: ； Calculate the branch carbon flow rate matrix, network loss carbon flow rate matrix, and load carbon flow rate vector for the water supply and return networks.

6. An electronic device, characterized in that, include: A memory, a processor, and a computer program stored in the memory and executable on the processor, the processor executing the program to implement the carbon emission metering method for a district heating system as described in any one of claims 1-5.

7. A computer-readable storage medium having a computer program stored thereon, characterized in that, The program is executed by the processor to implement the carbon emission metering method for a district heating system as described in any one of claims 1-5.

Images