METHOD AND APPARATUS FOR MEASURING CARBON EMISSION OF DISTRICT HEATING SYSTEM, ELECTRONIC DEVICE, AND MEDIUM

Abstract

The present disclosure discloses a method and apparatus for measuring carbon emission of a district heating system, an electronic device, and a medium. The method includes: obtaining a steady carbon emission amount of a current district heating system using a pre-trained steady carbon emission flow model; obtaining a dynamic carbon emission amount of the current district heating system using a pre-trained dynamic carbon emission flow model; and counting a carbon emission amount of the current district heating system based on the steady carbon emission amount and the dynamic carbon emission amount. Therefore, the present disclosure can effectively identify carbon emission details of each link of a source, a grid, and a load of the district heating system, clarify carbon emission responsibilities on both a source side and a load side, and realize accurate measurement of carbon emission of the district heating system, which has high application value.

Claims

1. A method for measuring carbon emission of a district heating system, the method comprising: obtaining a steady carbon emission amount of a current district heating system using a pre-trained steady carbon emission flow model, wherein the pre-trained steady carbon emission flow model is constructed based on a pipeline carbon flow rate, network loss carbon emission, a nodal carbon flow density, a pipeline carbon flow density, a heat source carbon flow rate, and a heat load carbon flow rate of the current district heating system; obtaining a dynamic carbon emission amount of the current district heating system using a pre-trained dynamic carbon emission flow model, wherein the pre-trained dynamic carbon emission flow model is constructed based on water element carbon flow rates at a plurality of time periods, actual outlet carbon flow rates of a pipeline at the plurality of time periods, network loss carbon flow rates at the plurality of time periods, and nodal carbon flow densities at the plurality of time periods of the current district heating system; and counting a carbon emission amount of the current district heating system based on the steady carbon emission amount and the dynamic carbon emission amount.

2. The method according to claim 1, further comprising: constructing the steady carbon emission flow model of the district heating system based on the pipeline carbon flow rate, the network loss carbon emission, the nodal carbon flow density, the pipeline carbon flow density, the heat source carbon flow rate, and the heat load carbon flow rate of the district heating system, wherein said constructing the steady carbon emission flow model of the district heating system based on the pipeline carbon flow rate, the network loss carbon emission, the nodal carbon flow density, the pipeline carbon flow density, the heat source carbon flow rate, and the heat load carbon flow rate of the district heating system is: determining the pipeline carbon flow rate of the district heating system comprising a carbon flow rate of a water supply network and a carbon flow rate of a water return network, the carbon flow rate of the water supply network being:
R.sub.k.sup.BHS,in=ρ.sub.k.sup.BHScm.sub.k.sup.ST.sub.k.sup.S,in,∀k∈Ω.sup.BH
R.sub.k.sup.BHS,out=ρ.sub.k.sup.BHScm.sub.k.sup.ST.sub.k.sup.S,out,∀k∈Ω.sup.BH, where R.sub.k.sup.BHS,in and R.sub.k.sup.BHS,out represent an inlet carbon flow rate and an outlet carbon flow rate of a pipeline k in the water supply network, respectively, unit: tCO.sub.2/h; ρ.sub.k.sup.BHS represents a carbon flow density of the pipeline k in the water supply network, unit: tCO.sub.2/MWh; c represents a specific heat capacity of water, unit: MWh/(kg.Math.° C.); m.sub.k.sup.S represents a mass flow of the pipeline k in the water supply network, unit: kg/h; T.sub.k.sup.S,in and T.sub.k.sup.S,out represent an inlet temperature and an outlet temperature of the pipeline k in the water supply network, unit: ° C.; and Ω.sup.BH represents a set of pipelines in the district heating system; and the carbon flow rate of the water return network being:
R.sub.k.sup.BHR,in=ρ.sub.k.sup.BHRcm.sub.k.sup.RT.sub.k.sup.R,in,∀k∈Ω.sup.BH
R.sub.k.sup.BHR,out=ρ.sub.k.sup.BHRcm.sub.k.sup.RT.sub.k.sup.R,out,∀k∈Ω.sup.BH, where R.sub.k.sup.BHR,in and R.sub.k.sup.BHR,out represent an inlet carbon flow rate and an outlet carbon flow rate of a pipeline k in the water return network, respectively, unit: tCO.sub.2/h; ρ.sub.k.sup.BHR represents a carbon flow density of the pipeline k in the water return network, unit: tCO.sub.2/h; m.sub.k.sup.R represents a mass flow of the pipeline k in the water return network, unit: kg/h; and T.sub.k.sup.R,in and T.sub.k.sup.R,out represent an inlet temperature and an outlet temperature of the pipeline k in the water return network, unit: ° C.; determining the network loss carbon emission of the district heating system comprising a network loss carbon flow rate of the water supply network and a network loss carbon flow rate of the water return network, said determining the network loss carbon emission of the district heating system being: determining a temperature difference of a pipeline in the water supply network and a temperature difference of a pipeline in the water return network:
T.sub.k.sup.S,Loss=T.sub.k.sup.S,in−T.sub.k.sup.S,out
T.sub.k.sup.R,Loss=T.sub.k.sup.R,in−T.sub.k.sup.R,out, where T.sub.k.sup.S,Loss and T.sub.k.sup.R,Loss represent a temperature difference between both ends of the pipeline k in the water supply network and a temperature difference between both ends of the pipeline k in the water return network, respectively, unit: ° C.; and determining the network loss carbon flow rate of the water supply network and the network loss carbon flow rate of the water return network:
R.sub.k.sup.BHS,Loss=ρ.sub.k.sup.BHScm.sub.k.sup.ST.sub.k.sup.S,Loss,∀k∈Ω.sup.BH
R.sub.k.sup.BHR,Loss=ρ.sub.k.sup.BHRcm.sub.k.sup.RT.sub.k.sup.R,Loss,∀k∈Ω.sup.BH, where R.sub.k.sup.BHS,Loss and R.sub.k.sup.BHR,Loss represent a network loss carbon flow rate of the pipeline k in the water supply network and a network loss carbon flow rate of the pipeline k in the water return network, respectively, unit: tCO.sub.2/h; determining the nodal carbon flow density of the district heating system, said determining the nodal carbon flow density of the district heating system comprising: determining, for each node in the district heating system, that conservation of mass and conservation of energy are satisfied at the node: $m_n^S=\underset{k\in {\Omega }_n^{BH+}}{.Math.}{m}_k^S,\forall n\in {\Omega }^{NH}$ $T_n^S{m}_n^S=\underset{k\in {\Omega }_n^{BH+}}{.Math.}{T}_k^{S,out}{m}_k^S,\forall n\in {\Omega }^{NH},$ where m.sub.n.sup.S represents a total mass flow flowing through a node n in the water supply network, unit: kg/h; Ω.sub.n.sup.BH+ represents a set of injection pipelines at the node n in the water supply network; Ω.sup.NH represents a set of nodes in the district heating system; and T.sub.n.sup.S represents a water flow temperature of the node n in the water supply network, unit: ° C.; determining, for each node in the district heating system, that conservation of carbon emission is satisfied at the node, a carbon flow rate of the node n being equal to a sum of outlet carbon flow rates of all the injection pipelines and network loss carbon flow rates allocated to the injection pipelines; $R_n^{NHS}=\underset{k\in {\Omega }_n^{BH+}}{.Math.}\left(R_k^{BHS,out}+{\mathrm{XR}}_k^{\mathrm{BHS},\mathrm{Loss}}\right)=\underset{k\in {\Omega }_n^{BH+}}{.Math.}{\rho }_k^{\mathrm{BHS}}c{m}_k^S\left(T_k^{S,out}+X{T}_k^{S,\mathrm{Loss}}\right),\forall n\in {\Omega }^{NH},$ where R.sub.n.sup.NHS represents the carbon flow rate of the node n in the water supply network, unit: tCO.sub.2/h; and X represents an allocating coefficient of a network loss carbon flow rate of an injection pipeline; determining a carbon flow density of the node n in the water supply network: ${\rho }_n^{NHS}=\frac{R_n^{NHS}}{c{m}_n^S{T}_n^S}=\frac{\underset{k\in {\Omega }_n^{BH+}}{.Math.}{\rho }_k^{BHS}{m}_k^S\left(T_k^{S,\mathrm{out}}+X{T}_k^{S,\mathrm{Loss}}\right)}{\underset{k\in {\Omega }_n^{BH+}}{.Math.}{m}_k^S{T}_k^{S,\mathrm{out}}},\forall n\in {\Omega }^{NH},$ where ρ.sub.n.sup.NHS represents the carbon flow density of the node n in the water supply network, unit: tCO.sub.2/MWh; and determining a carbon flow density of the node n in the water return network: ${\rho }_n^{NHR}=\frac{\underset{k\in {\Omega }_n^{\mathrm{BH}-}}{.Math.}{\rho }_k^{BHR}{m}_k^R\left(T_k^{R,\mathrm{out}}+{\mathrm{XT}}_k^{R,\mathrm{Loss}}\right)}{\underset{k\in {\Omega }_n^{\mathrm{BH}-}}{.Math.}{m}_k^R{T}_k^{R,\mathrm{out}}},\forall n\in {\Omega }^{NH},$ where ρ.sub.n.sup.NHR represents the carbon flow density of the node n in the water return network, unit: tCO.sub.2/MWh; and Ω.sub.n.sup.BH− represents a set of outflow pipelines at the node n in the water supply network; determining the pipeline carbon flow density of the district heating system:
ρ.sub.k.sup.BHS=ρ.sub.n.sup.NHS,n=Γ.sub.k.sup.NH+,∀k∈Ω.sup.BH
ρ.sub.k.sup.BHR=ρ.sub.n.sup.NHR,n=Γ.sub.k.sup.NH−,∀k∈Ω.sup.BH, where Γ.sub.k.sup.NH+ and Γ.sub.k.sup.NH+ represent an injection node and an outflow node of the pipeline k in the water supply network, respectively; determining the heat source carbon flow rate of the district heating system, said determining the heat source carbon flow rate of the district heating system comprising: determining a heat output of a heat source:
Q.sub.i=cm.sub.n.sup.S(T.sub.n.sup.S−T.sub.n.sup.R),n=Γ.sub.i.sup.NH,∀i∈Ω.sup.GH, where Q.sub.i represents a heat output of a heat source i, unit: MW; Γ.sub.i.sup.NH represents a node where the heat source i is located; and Ω.sup.GH represents a set of heat sources; and determining that conservation of carbon emission at a heat source node is satisfied:
ρ.sub.n.sup.NHScm.sub.n.sup.ST.sub.n.sup.S=ρ.sub.i.sup.GHQ.sub.i+ρ.sub.n.sup.NHRcm.sub.k.sup.ST.sub.k.sup.R,n=Γ.sub.i.sup.NH,∀i∈Ω.sup.GH, where ρ.sub.i.sup.GH represents a carbon flow density of the heat source i, unit: tCO.sub.2/MWh; and determining the heat load carbon flow rate of the district heating system, said determining the heat load carbon flow rate of the district heating system comprising: determining a heat load demand:
q.sub.j=cm.sub.n.sup.S(T.sub.n.sup.S−T.sub.n.sup.R),n=Γ.sub.j.sup.NH,∀j∈Ω.sup.LH, where q.sub.j represents a heat demand of a heat load j, unit: MW; Γ.sub.j.sup.NH represents a node where the heat load j is located; and Ω.sup.LH represents a set of heat loads; determining a carbon flow density of a heat load node:
ρ.sub.n.sup.NHR=ρ.sub.n.sup.NHS,n=Γ.sub.j.sup.NH,∀j∈Ω.sup.LH; and determining the heat load carbon flow rate:
R.sub.j.sup.LH=ρ.sub.n.sup.NHSq.sub.j=ρ.sub.n.sup.NHScm.sub.n.sup.S(T.sub.n.sup.S−T.sub.n.sup.R),n=Γ.sub.j.sup.NH,∀j∈Ω.sup.LH, where R.sub.j.sup.LH represents a carbon flow rate of the heat load j, unit: tCO.sub.2/h.

3. The method according to claim 2, further comprising: calculating a matrix representation of the steady carbon emission flow model of the district heating system.

4. The method according to claim 3, wherein said calculating the matrix representation of the steady carbon emission flow model of the district heating system comprises: constructing a branch heat flow matrix, a branch network loss matrix, and a nodal heat flow flux matrix of the district heating system; and calculating a nodal carbon flow density vector of a heat network, and then calculating a branch carbon flow rate matrix, a network loss carbon flow rate matrix, and a load carbon flow rate vector of each of the water supply network and the water return network.

5. The method according to claim 4, wherein said calculating the matrix representation of the steady carbon emission flow model of the district heating system comprises: constructing the branch heat flow matrix of the district heating system comprising a branch heat flow matrix of the water supply network and a branch heat flow matrix of the water return network, wherein said constructing the branch heat flow matrix of the district heating system comprises: determining elements of the branch heat flow matrix of the water supply network:
Q.sub.B,ij.sup.S=cm.sub.k.sup.ST.sub.k.sup.S,out,Q.sub.B,ji.sup.S=0, where Q.sub.B,ij.sup.S and Q.sub.B,ji.sup.S represent elements in a branch heat flow matrix Q.sub.B.sup.S of the water supply network; and determining elements of the branch heat flow matrix of the water return network:
Q.sub.B,ij.sup.R=cm.sub.k.sup.RT.sub.k.sup.R,out,Q.sub.B,ji.sup.R=0, where Q.sub.B,ij.sup.R and Q.sub.B,ji.sup.R represent elements in a branch heat flow matrix Q.sub.B.sup.R of the water return network; constructing the branch network loss matrix of the district heating system comprising a branch network loss matrix of the water supply network and a branch network loss matrix of the water return network, wherein said constructing the branch network loss matrix of the district heating system comprises: determining elements of the branch network loss matrix of the water supply network:
Q.sub.BL,ij.sup.S=cm.sub.k.sup.ST.sub.k.sup.S,Loss,Q.sub.BL,ji.sup.S=0 where Q.sub.BL,ij.sup.S and Q.sub.BL,ji.sup.S represent elements in a branch network loss matrix Q.sub.BL.sup.S of the water supply network; and determining elements of the branch network loss matrix of the water return network:
Q.sub.BL,ij.sup.R=cm.sub.k.sup.RT.sub.k.sup.R,Loss,Q.sub.BL,ji.sup.R=0, where Q.sub.BL,ij.sup.R and Q.sub.BL,ji.sup.R represent elements in a branch network loss matrix Q.sub.BL.sup.R of the water return network; constructing the nodal heat flow flux matrix of the district heating system comprising a nodal heat flow flux matrix of the water supply network and a nodal heat flow flux matrix of the water return network, wherein said constructing the nodal heat flow flux matrix of the district heating system comprises: determining, for a node in no connection to the heat source, the nodal heat flow flux matrix of the water supply network:
Q.sub.N.sup.S=diag{ζ.sub.N.sub.NHQ.sub.B.sup.S}, where Q.sub.N.sup.S represents the nodal heat flow flux matrix of the water supply network; ζ.sub.N.sub.NH represents a coefficient matrix of a branch heat flow; and N.sub.NH represents a number of nodes of the district heating system; determining, for a node in no connection to the heat source, the nodal heat flow flux matrix of the water return network:
Q.sub.N.sup.R=diag{ζ.sub.N.sub.NHQ.sub.B.sup.R}, where Q.sub.N.sup.R represents the nodal heat flow flux matrix of the water return network; determining, for a node connected to the heat source, a nodal integrated energy flow flux matrix:
{circumflex over (Q)}.sub.N.sup.S=diag{ζ.sub.N.sub.NHQ.sub.B.sup.S+ζ.sub.N.sub.GHQ.sub.G}, where {circumflex over (Q)}.sub.N.sup.S represents the nodal integrated energy flow flux matrix; ζ.sub.N.sub.GH represents a coefficient matrix of a heat flow injected by the heat source; and N.sub.GH represents a number of heat sources of the district heating system; determining, for all nodes, that a total injected carbon emission of the nodes is equal to a sum of injected carbon flow rates of all branches connected to the nodes:
Q.sub.N.sup.Sρ.sup.NHS=(Q.sub.B.sup.S+XQ.sub.BL.sup.S).sup.Tρ.sup.NHS
Q.sub.N.sup.Rρ.sup.NHR=(Q.sub.B.sup.R+XQ.sub.BL.sup.R).sup.Tρ.sup.NHR, where ρ.sup.NHS represents a matrix formed by the carbon flow density ρ.sub.n.sup.NHS of the node n in the water supply network; and ρ.sup.NHR represents a matrix formed by the carbon flow density ρ.sub.n.sup.NHR of the node n in the water return network; determining, for the heat load node, that the heat load node has an equal carbon flow density in the water supply network and the water return network:
Bρ.sup.NHS=Bρ.sup.NHR, where B represents a heat load-node association matrix, when the heat load j is connected to the node n, B.sub.jn=1, otherwise B.sub.jn=0; determining, for the heat source node, a matrix relation of the heat source node based on a conservation of carbon emission as:
C{circumflex over (Q)}.sub.N.sup.Sρ.sup.NHS=CQ.sub.N.sup.Rρ.sup.NHR+Q.sub.G.sup.Tρ.sup.GH, where C represents a 0-1 matrix associated with the heat source node, when the node n is connected to the heat source, C.sub.nn=1, otherwise C.sub.nn=0; determining the nodal carbon flow density vector of the heat network as: $\left[\begin{array}{c}{\rho }^{\mathrm{NHS}}\\ {\rho }^{NHR}\end{array}\right]={\left[\begin{array}{cc}{Q}_N^S-{\left(Q_B^S+X{Q}_{BL}^S\right)}^T& 0\\ 0& Q_N^R-{\left(Q_B^R+{\mathrm{XQ}}_{BL}^R\right)}^T\\ B& -B\\ C{\hat{Q}}_N^S& -C{Q}_N^R\end{array}\right]}^{-1}\left[\begin{array}{c}0\\ 0\\ 0\\ Q_G^T{\rho }^{GH}\end{array}\right];$ and calculating the branch carbon flow rate matrix, the network loss carbon flow rate matrix, and the load carbon flow rate vector of each of the water supply network and the water return network.

6. The method according to claim 5, further comprising: constructing the dynamic carbon emission flow model of the district heating system based on the water element carbon flow rates at the plurality of time periods, the actual outlet carbon flow rates of the pipeline at the plurality of time periods, the network loss carbon flow rates at the plurality of time periods, and the nodal carbon flow densities at the plurality of time periods, wherein said constructing the dynamic carbon emission flow model of the district heating system based on the water element carbon flow rates at the plurality of time periods, the actual outlet carbon flow rates of the pipeline at the plurality of time periods, the network loss carbon flow rates at the plurality of time periods, and the nodal carbon flow densities at the plurality of time periods is: determining the water element carbon flow rates at the plurality of time periods: ${\tilde{R}}_{k,t}^{\mathrm{BHS},\mathrm{out}}=\frac{R_{k,t-{\delta }_{k,t}}^{\mathrm{BHS},in}\left(B_{k,t}-\sigma M_k\right)}{\left(m_{k,t-{\delta }_{k,t}}^S\Delta t\right)}+\underset{\tau =t-{\phi }_{k,t}+1}{\overset{t-{\delta }_{k,t}-1}{.Math.}}{R}_{k,\tau }^{\mathrm{BHS},in}+\frac{R_{k,t-{\delta }_{k,t}}^{\mathrm{BHS},in}\left(m_{k,t}^S\Delta t+\sigma M_k-A_{k,t}\right)}{\left(m_{k,t-{\phi }_{k,t}}^S\Delta t\right)},$ where {tilde over (R)}.sub.k,t.sup.BHS,out represents a water element carbon flow rate at a time period t; σ represents a density of water; M.sub.k represents a volume of the pipeline k; φ.sub.k,t represents an injection time period of an earliest water element component contained in a water flow flowing out of the pipeline k at the time period t; δ.sub.k,t represents an injection time period of a latest water element component contained in the water flow flowing out of the pipeline k at the time period t; A.sub.k,t represents a total injection water flow amount for φ.sub.k,t time periods before the time period t; B.sub.k,t represents a total injection water flow amount from a time period t−δ.sub.k,t to the time period t; and expressions for A.sub.k,t and B.sub.k,t are: $A_{k,t}=\{\begin{array}{cc}\underset{\tau =t-{\phi }_{k,t}+1}{\overset{t}{.Math.}}{m}_{k,\tau }^S\Delta t,& {\phi }_{k,t}\ge {\delta }_{k,t}+1\\ B_{k,t},& {\phi }_{k,t}<{\delta }_{k,t}+1\end{array},$ $B_{k,t}=\underset{\tau =t-{\delta }_{k,t}}{\overset{t}{.Math.}}{m}_{k,\tau }^S\Delta t;$ and determining the actual outlet carbon flow rates of the pipeline at the plurality of time periods and the network loss carbon flow rates at the plurality of time periods, said determining the actual outlet carbon flow rates of the pipeline at the plurality of time periods and the network loss carbon flow rates at the plurality of time periods comprising: calculating an actual temperature of an outlet water flow of the pipeline based on a conveying loss of the pipeline: $T_{k,t}^{S,\mathrm{out}}={\tilde{T}}_{k,t}^{S,\mathrm{out}}\exp \left(-\frac{\lambda L_k}{{\mathrm{cm}}_{k,t}^S}\right),$ where T.sub.k,t.sup.S,out represents the actual temperature of the outlet water flow of the pipeline; {tilde over (T)}.sub.k,t.sup.S,out represents a weighted average of temperatures of injection water flows at previous time periods; λ represents a thermal conductivity coefficient of the pipeline; and L.sub.k represents a length of the pipeline k; determining an actual outlet carbon flow rate of the pipeline at the time period t: $R_{k,t}^{\mathrm{BHS},\mathrm{out}}={\tilde{R}}_{k,t}^{{\mathrm{BHS},\mathrm{out}}_{}}\frac{T_{k,t}^{S,\mathrm{out}}}{{\tilde{T}}_{k,t}^{{S,out}_{}}};$ and determining a network loss carbon flow rate of the pipeline at the time period t: $R_{k,t}^{\mathrm{BHS},\mathrm{Loss}}={\tilde{R}}_{k,t}^{\mathrm{BHS},\mathrm{Loss}}\left(1-\frac{T_{k,t}^{S,\mathrm{out}}}{{\tilde{T}}_{k,t}^{S,\mathrm{out}}}\right);$ determining the nodal carbon flow densities at the plurality of time periods: ${\rho }_{n,t}^{NHS}=\frac{\underset{k\in {\Omega }_n^{BH+}}{.Math.}\left(R_{k,t}^{\mathrm{BHS},\mathrm{out}}+X{R}_{k,t}^{\mathrm{BHS},\mathrm{Loss}}\right)}{\underset{k\in {\Omega }_n^{BH+}}{.Math.}c{m}_{k,t}^S{T}_{k,t}^{S,\mathrm{out}}};$ and determining, for the water return network, the nodal carbon flow densities at the plurality of time periods: ${\rho }_{n,t}^{\mathrm{NHR}}=\frac{\underset{k\in {\Omega }_n^{\mathrm{BH}-}}{.Math.}\left(R_{k,t}^{\mathrm{BHR},\mathrm{out}}+X{R}_{k,t}^{\mathrm{BHR},\mathrm{Loss}}\right)}{\underset{k\in {\Omega }_n^{\mathrm{BH}-}}{.Math.}c{m}_{k,t}^R{T}_{k,t}^{R,\mathrm{out}}}.$

7. An apparatus for measuring carbon emission of a district heating system, the apparatus comprising: a plurality of heat source-side carbon meters, wherein each of the plurality of heat source-side carbon meters is connected to a heat source and configured to measure heat source-side carbon emission; a plurality of heat pipeline carbon meters, wherein each of the plurality of heat pipeline carbon meters is connected to a heat pipeline and the heat source-side carbon meter, and the heat pipeline carbon emission meter is configured to measure heat pipeline carbon emission; a plurality of user-side carbon meters, wherein each of the plurality of user-side carbon meters has one end connected to a heat user end and another end connected to the heat pipeline carbon meter, and is configured to obtain a carbon potential of a node where a user is located and obtain a carbon emission amount resulted from heat consumption of the user; and a central server, wherein the central server is connected to the plurality of heat source-side carbon meters and the plurality of heat pipeline carbon meters, and configured to calculate a carbon emission intensity of the heat source based on coal consumption data of the heat source and an emission coefficient of the coal, calculate a distribution of a carbon emission flow in the heat pipeline, and calculate a carbon emission amount corresponding to a pipeline connected to the node based on data of the heat pipeline carbon meter.

8. The apparatus according to claim 7, wherein each of the heat source-side carbon meter, the heat pipeline carbon meter, and the user-side carbon meter is further configured to display in real time a result of the carbon emission of the district heating system calculated by the central server.

9. An electronic device, comprising: a memory; a processor; and a computer program stored in the memory and executable on the processor, wherein the processor, when executing the computer program, implements: obtaining a steady carbon emission amount of a current district heating system using a pre-trained steady carbon emission flow model, wherein the pre-trained steady carbon emission flow model is constructed based on a pipeline carbon flow rate, network loss carbon emission, a nodal carbon flow density, a pipeline carbon flow density, a heat source carbon flow rate, and a heat load carbon flow rate of the current district heating system; obtaining a dynamic carbon emission amount of the current district heating system using a pre-trained dynamic carbon emission flow model, wherein the pre-trained dynamic carbon emission flow model is constructed based on water element carbon flow rates at a plurality of time periods, actual outlet carbon flow rates of a pipeline at the plurality of time periods, network loss carbon flow rates at the plurality of time periods, and nodal carbon flow densities at the plurality of time periods of the current district heating system; and counting a carbon emission amount of the current district heating system based on the steady carbon emission amount and the dynamic carbon emission amount.

10. The electronic device according to claim 9, further comprising: constructing the steady carbon emission flow model of the district heating system based on the pipeline carbon flow rate, the network loss carbon emission, the nodal carbon flow density, the pipeline carbon flow density, the heat source carbon flow rate, and the heat load carbon flow rate of the district heating system, wherein said constructing the steady carbon emission flow model of the district heating system based on the pipeline carbon flow rate, the network loss carbon emission, the nodal carbon flow density, the pipeline carbon flow density, the heat source carbon flow rate, and the heat load carbon flow rate of the district heating system is: determining the pipeline carbon flow rate of the district heating system comprising a carbon flow rate of a water supply network and a carbon flow rate of a water return network, the carbon flow rate of the water supply network being:
R.sub.k.sup.BHS,in=ρ.sub.k.sup.BHScm.sub.k.sup.ST.sub.k.sup.S,in,∀k∈Ω.sup.BH
R.sub.k.sup.BHS,out=ρ.sub.k.sup.BHScm.sub.k.sup.ST.sub.k.sup.S,out,∀k∈Ω.sup.BH, where R.sub.k.sup.BHS,in and R.sub.k.sup.BHS,out represent an inlet carbon flow rate and an outlet carbon flow rate of a pipeline k in the water supply network, respectively, unit: tCO.sub.2/h; ρ.sub.k.sup.BHS represents a carbon flow density of the pipeline k in the water supply network, unit: tCO.sub.2/MWh; c represents a specific heat capacity of water, unit: MWh/(kg.Math.° C.); m.sub.k.sup.S represents a mass flow of the pipeline k in the water supply network, unit: kg/h; T.sub.k.sup.S,in and T.sub.k.sup.S,out represent an inlet temperature and an outlet temperature of the pipeline k in the water supply network, unit: ° C.; and Ω.sup.BH represents a set of pipelines in the district heating system; and the carbon flow rate of the water return network being:
R.sub.k.sup.BHR,in=ρ.sub.k.sup.BHRcm.sub.k.sup.RT.sub.k.sup.R,in,∀k∈Ω.sup.BH
R.sub.k.sup.BHR,out=ρ.sub.k.sup.BHRcm.sub.k.sup.RT.sub.k.sup.R,out,∀k∈Ω.sup.BH, where R.sub.k.sup.BHR,in and R.sub.k.sup.BHR,out represent an inlet carbon flow rate and an outlet carbon flow rate of a pipeline k in the water return network, respectively, unit: tCO.sub.2/h; ρ.sub.k.sup.BHR represents a carbon flow density of the pipeline k in the water return network, unit: tCO.sub.2/h; m.sub.k.sup.R represents a mass flow of the pipeline k in the water return network, unit: kg/h; and T.sub.k.sup.R,in and T.sub.k.sup.R,out represent an inlet temperature and an outlet temperature of the pipeline k in the water return network, unit: ° C.; determining the network loss carbon emission of the district heating system comprising a network loss carbon flow rate of the water supply network and a network loss carbon flow rate of the water return network, said determining the network loss carbon emission of the district heating system being: determining a temperature difference of a pipeline in the water supply network and a temperature difference of a pipeline in the water return network:
T.sub.k.sup.S,Loss=T.sub.k.sup.S,in−T.sub.k.sup.S,out
T.sub.k.sup.R,Loss=T.sub.k.sup.R,in−T.sub.k.sup.R,out, where T.sub.k.sup.S,Loss and T.sub.k.sup.R,Loss represent a temperature difference between both ends of the pipeline k in the water supply network and a temperature difference between both ends of the pipeline k in the water return network, respectively, unit: ° C.; and determining the network loss carbon flow rate of the water supply network and the network loss carbon flow rate of the water return network:
R.sub.k.sup.BHS,Loss=ρ.sub.k.sup.BHScm.sub.k.sup.ST.sub.k.sup.S,Loss,∀k∈Ω.sup.BH
R.sub.k.sup.BHR,Loss=ρ.sub.k.sup.BHRcm.sub.k.sup.RT.sub.k.sup.R,Loss,∀k∈Ω.sup.BH, where R.sub.k.sup.BHS,Loss and R.sub.k.sup.BHR,Loss represent a network loss carbon flow rate of the pipeline k in the water supply network and a network loss carbon flow rate of the pipeline k in the water return network, respectively, unit: tCO.sub.2/h; determining the nodal carbon flow density of the district heating system, said determining the nodal carbon flow density of the district heating system comprising: determining, for each node in the district heating system, that conservation of mass and conservation of energy are satisfied at the node: $m_n^S=\underset{k\in {\Omega }_n^{BH+}}{.Math.}{m}_k^S,\forall n\in {\Omega }^{NH}$ $T_n^S{m}_n^S=\underset{k\in {\Omega }_n^{BH+}}{.Math.}{T}_k^{S,out}{m}_k^S,\forall n\in {\Omega }^{NH},$ where m.sub.n.sup.S represents a total mass flow flowing through a node n in the water supply network, unit: kg/h; Ω.sub.n.sup.BH+ represents a set of injection pipelines at the node n in the water supply network; Ω.sup.NH represents a set of nodes in the district heating system; and T.sub.n.sup.S represents a water flow temperature of the node n in the water supply network, unit: ° C.; determining, for each node in the district heating system, that conservation of carbon emission is satisfied at the node, a carbon flow rate of the node n being equal to a sum of outlet carbon flow rates of all the injection pipelines and network loss carbon flow rates allocated to the injection pipelines; $R_n^{NHS}=\underset{k\in {\Omega }_n^{BH+}}{.Math.}\left(R_k^{BHS,out}+{\mathrm{XR}}_k^{\mathrm{BHS},\mathrm{Loss}}\right)=\underset{k\in {\Omega }_n^{BH+}}{.Math.}{\rho }_k^{\mathrm{BHS}}{\mathrm{cm}}_k^S\left(T_k^{S,out}+X{T}_k^{S,\mathrm{Loss}}\right),$ $\forall n\in {\Omega }^{NH},$ where R.sub.n.sup.NHS represents the carbon flow rate of the node n in the water supply network, unit: tCO.sub.2/h; and X represents an allocating coefficient of a network loss carbon flow rate of an injection pipeline; determining a carbon flow density of the node n in the water supply ${\rho }_n^{NHS}=\frac{R_n^{NHS}}{c{m}_n^S{T}_n^S}=\frac{\underset{k\in {\Omega }_n^{BH+}}{.Math.}{\rho }_k^{BHS}{m}_k^S\left(T_k^{S,\mathrm{out}}+X{T}_k^{S,\mathrm{Loss}}\right)}{\underset{k\in {\Omega }_n^{BH+}}{.Math.}{m}_k^S{T}_k^{S,\mathrm{out}}},\forall n\in {\Omega }^{NH},$ where R.sub.n.sup.NHS represents the carbon flow density of the node n in the water supply network, unit: tCO.sub.2/MWh; and determining a carbon flow density of the node n in the water return network: ${\rho }_n^{NHR}=\frac{\underset{k\in {\Omega }_n^{\mathrm{BH}-}}{.Math.}{\rho }_k^{BHR}{m}_k^R\left(T_k^{R,\mathrm{out}}+{\mathrm{XT}}_k^{R,\mathrm{Loss}}\right)}{\underset{k\in {\Omega }_n^{\mathrm{BH}-}}{.Math.}{m}_k^R{T}_k^{R,\mathrm{out}}},\forall n\in {\Omega }^{NH},$ where ρ.sub.n.sup.NHR represents the carbon flow density of the node n in the water return network, unit: tCO.sub.2/MWh; and Ω.sub.n.sup.BH− represents a set of outflow pipelines at the node n in the water supply network; determining the pipeline carbon flow density of the district heating system:
ρ.sub.k.sup.BHS=ρ.sub.n.sup.NHS,n=Γ.sub.k.sup.NH+,∀k∈Ω.sup.BH
ρ.sub.k.sup.BHR=ρ.sub.n.sup.NHR,n=Γ.sub.k.sup.NH−,∀k∈Ω.sup.BH, where Γ.sub.k.sup.NH+ and Γ.sub.k.sup.NH+ represent an injection node and an outflow node of the pipeline k in the water supply network, respectively; determining the heat source carbon flow rate of the district heating system, said determining the heat source carbon flow rate of the district heating system comprising: determining a heat output of a heat source:
Q.sub.i=cm.sub.n.sup.S(T.sub.n.sup.S−T.sub.n.sup.R),n=Γ.sub.i.sup.NH,∀i∈Ω.sup.GH, where Q.sub.i represents a heat output of a heat source i, unit: MW; Γ.sub.i.sup.NH represents a node where the heat source i is located; and Ω.sup.GH represents a set of heat sources; and determining that conservation of carbon emission at a heat source node is satisfied:
ρ.sub.n.sup.NHScm.sub.n.sup.ST.sub.n.sup.S=ρ.sub.i.sup.GHQ.sub.i+ρ.sub.n.sup.NHRcm.sub.n.sup.ST.sub.n.sup.R,n=Γ.sub.i.sup.NH,∀i∈ΩGH, where ρ.sub.i.sup.GH represents a carbon flow density of the heat source i, unit: tCO.sub.2/MWh; and determining the heat load carbon flow rate of the district heating system, said determining the heat load carbon flow rate of the district heating system comprising: determining a heat load demand:
q.sub.j=cm.sub.n.sup.S(T.sub.n.sup.S−T.sub.n.sup.R),n=Γ.sub.j.sup.NH,∀j∈Ω.sup.LH, where q.sub.j represents a heat demand of a heat load j, unit: MW; Γ.sub.j.sup.NH represents a node where the heat load j is located; and Ω.sup.LH represents a set of heat loads; determining a carbon flow density of a heat load node:
ρ.sub.n.sup.NHR=ρ.sub.n.sup.NHS,n=Γ.sub.j.sup.NH,∀j∈Ω.sup.LH; and determining the heat load carbon flow rate:
R.sub.j.sup.LH=ρ.sub.n.sup.NHSq.sub.j=ρ.sub.n.sup.NHScm.sub.n.sup.S(T.sub.n.sup.S−T.sub.n.sup.R),n=Γ.sub.j.sup.NH,∀j∈Ω.sup.LH, where R.sub.j.sup.LH represents a carbon flow rate of the heat load j, unit: tCO.sub.2/h.

11. The electronic device according to claim 10, further comprising: calculating a matrix representation of the steady carbon emission flow model of the district heating system.

12. The electronic device according to claim 11, wherein said calculating the matrix representation of the steady carbon emission flow model of the district heating system comprises: constructing a branch heat flow matrix, a branch network loss matrix, and a nodal heat flow flux matrix of the district heating system; and calculating a nodal carbon flow density vector of a heat network, and then calculating a branch carbon flow rate matrix, a network loss carbon flow rate matrix, and a load carbon flow rate vector of each of the water supply network and the water return network.

13. The electronic device according to claim 12, wherein said calculating the matrix representation of the steady carbon emission flow model of the district heating system comprises: constructing the branch heat flow matrix of the district heating system comprising a branch heat flow matrix of the water supply network and a branch heat flow matrix of the water return network, wherein said constructing the branch heat flow matrix of the district heating system comprises: determining elements of the branch heat flow matrix of the water supply network:
Q.sub.B,ij.sup.S=cm.sub.k.sup.ST.sub.k.sup.S,out,Q.sub.B,ji.sup.S=0, where Q.sub.B,ij.sup.S and Q.sub.B,ji.sup.S represent elements in a branch heat flow matrix Q.sub.B.sup.S of the water supply network; and determining elements of the branch heat flow matrix of the water return network:
Q.sub.B,ij.sup.R=cm.sub.k.sup.RT.sub.k.sup.R,out,Q.sub.B,ji.sup.R=0, where Q.sub.B,ij.sup.R and Q.sub.B,ji.sup.R represent elements in a branch heat flow matrix Q.sub.B.sup.R of the water return network; constructing the branch network loss matrix of the district heating system comprising a branch network loss matrix of the water supply network and a branch network loss matrix of the water return network, wherein said constructing the branch network loss matrix of the district heating system comprises: determining elements of the branch network loss matrix of the water supply network:
Q.sub.BL,ij.sup.S=cm.sub.k.sup.ST.sub.k.sup.S,Loss,Q.sub.BL,ji.sup.S=0, where Q.sub.BL,ij.sup.S and Q.sub.BL,ji.sup.S represent elements in a branch network loss matrix Q.sub.BL.sup.S of the water supply network; and determining elements of the branch network loss matrix of the water return network:
Q.sub.BL,ij.sup.R=cm.sub.k.sup.RT.sub.k.sup.R,Loss,Q.sub.BL,ji.sup.R=0, where Q.sub.BL,ij.sup.R and Q.sub.BL,ji.sup.R represent elements in a branch network loss matrix Q.sub.BL.sup.R of the water return network; constructing the nodal heat flow flux matrix of the district heating system comprising a nodal heat flow flux matrix of the water supply network and a nodal heat flow flux matrix of the water return network, wherein said constructing the nodal heat flow flux matrix of the district heating system comprises: determining, for a node in no connection to the heat source, the nodal heat flow flux matrix of the water supply network:
Q.sub.N.sup.S=diag{ζ.sub.N.sub.NHQ.sub.B.sup.S}, where Q.sub.N.sup.S represents the nodal heat flow flux matrix of the water supply network; ζ.sub.N.sub.NH represents a coefficient matrix of a branch heat flow; and N.sub.NH represents a number of nodes of the district heating system; determining, for a node in no connection to the heat source, the nodal heat flow flux matrix of the water return network:
Q.sub.N.sup.R=diag{ζ.sub.N.sub.NHQ.sub.B.sup.R}, where Q.sub.N.sup.R represents the nodal heat flow flux matrix of the water return network; determining, for a node connected to the heat source, a nodal integrated energy flow flux matrix:
{circumflex over (Q)}.sub.N.sup.S=diag{ζ.sub.N.sub.NHQ.sub.B.sup.S+ζ.sub.N.sub.GHQ.sub.G}, where {circumflex over (Q)}.sub.N.sup.S represents the nodal integrated energy flow flux matrix; ζ.sub.N.sub.GH represents a coefficient matrix of a heat flow injected by the heat source; and N.sub.GH represents a number of heat sources of the district heating system; determining, for all nodes, that a total injected carbon emission of the nodes is equal to a sum of injected carbon flow rates of all branches connected to the nodes:
Q.sub.N.sup.Sρ.sup.NHS=(Q.sub.B.sup.S+XQ.sub.BL.sup.S).sup.Tρ.sup.NHS
Q.sub.N.sup.Rρ.sup.NHR=(Q.sub.B.sup.R+XQ.sub.BL.sup.R).sup.Tρ.sup.NHR, where ρ.sup.NHS represents a matrix formed by the carbon flow density ρ.sub.n.sup.NHS of the node n in the water supply network; and ρ.sup.NHR represents a matrix formed by the carbon flow density ρ.sub.n.sup.NHR of the node n in the water return network; determining, for the heat load node, that the heat load node has an equal carbon flow density in the water supply network and the water return network:
Bρ.sup.NHS=Bρ.sup.NHR, where B represents a heat load-node association matrix, when the heat load j is connected to the node n, B.sub.jn=1, otherwise B.sub.jn=0; determining, for the heat source node, a matrix relation of the heat source node based on a conservation of carbon emission as:
C{circumflex over (Q)}.sub.N.sup.Sρ.sup.NHS=CQ.sub.N.sup.Rρ.sup.NHR+Q.sub.G.sup.Tρ.sup.GH, where C represents a 0-1 matrix associated with the heat source node, when the node n is connected to the heat source, C.sub.nn=1, otherwise C.sub.nn=0; determining the nodal carbon flow density vector of the heat network as: $\left[\begin{array}{c}{\rho }^{NHS}\\ {\rho }^{NHR}\end{array}\right]={\left[\begin{array}{cc}{Q}_N^S-{\left(Q_B^S+X{Q}_{BL}^S\right)}^T& 0\\ 0& Q_N^R-{\left(Q_B^R+{\mathrm{XQ}}_{BL}^R\right)}^T\\ B& -B\\ C{\hat{Q}}_N^S& -C{Q}_N^R\end{array}\right]}^{-1}\left[\begin{array}{c}0\\ 0\\ 0\\ Q_G^T{\rho }^{GH}\end{array}\right];$ calculating the branch carbon flow rate matrix, the network loss carbon flow rate matrix, and the load carbon flow rate vector of each of the water supply network and the water return network.

14. The electronic device according to claim 13, further comprising: constructing the dynamic carbon emission flow model of the district heating system based on the water element carbon flow rates at the plurality of time periods, the actual outlet carbon flow rates of the pipeline at the plurality of time periods, the network loss carbon flow rates at the plurality of time periods, and the nodal carbon flow densities at the plurality of time periods, wherein said constructing the dynamic carbon emission flow model of the district heating system based on the water element carbon flow rates at the plurality of time periods, the actual outlet carbon flow rates of the pipeline at the plurality of time periods, the network loss carbon flow rates at the plurality of time periods, and the nodal carbon flow densities at the plurality of time periods is: determining the water element carbon flow rates at the plurality of time periods: ${\tilde{R}}_{k,t}^{\mathrm{BHS},\mathrm{out}}=\frac{R_{k,t-{\delta }_{k,t}}^{\mathrm{BHS},\mathrm{in}}\left(B_{k,t}-\sigma M_k\right)}{\left(m_{k,t-{\delta }_{k,t}}^S\Delta t\right)}+\underset{\tau =t-{\phi }_{k,t}+1}{\overset{t-{\delta }_{k,t}-1}{.Math.}}{R}_{k,\tau }^{\mathrm{BHS},\mathrm{in}}+\frac{R_{k,t-{\delta }_{k,t}}^{\mathrm{BHS},\mathrm{in}}\left(m_{k,t}^S\Delta t+\sigma M_k-A_{k,t}\right)}{\left(m_{k,t-{\phi }_{k,t}}^S\Delta t\right)},$ where {tilde over (R)}.sub.k,t.sup.BHS,out represents a water element carbon flow rate at a time period t; σ represents a density of water; M.sub.k represents a volume of the pipeline k; φ.sub.k,t represents an injection time period of an earliest water element component contained in a water flow flowing out of the pipeline k at the time period t; δ.sub.k,t represents an injection time period of a latest water element component contained in the water flow flowing out of the pipeline k at the time period t; A.sub.k,t represents a total injection water flow amount for φ.sub.k,t time periods before the time period t; B.sub.k,t represents a total injection water flow amount from a time period t−δ.sub.k,t to the time period t; and expressions for A.sub.k,t and B.sub.k,t are: $A_{k,t}=\{\begin{array}{cc}\underset{\tau =t-{\phi }_{k,t}+1}{\overset{t}{.Math.}}{m}_{k,\tau }^S\Delta t,& {\phi }_{k,t}\ge {\delta }_{k,t}+1\\ B_{k,t},& {\phi }_{k,t}<{\delta }_{k,t}+1\end{array},$ $B_{k,t}=\underset{\tau =t-{\delta }_{k,t}}{\overset{t}{.Math.}}{m}_{k,\tau }^S\Delta t;$ determining the actual outlet carbon flow rates of the pipeline at the plurality of time periods and the network loss carbon flow rates at the plurality of time periods, said determining the actual outlet carbon flow rates of the pipeline at the plurality of time periods and the network loss carbon flow rates at the plurality of time periods comprising: calculating an actual temperature of an outlet water flow of the pipeline based on a conveying loss of the pipeline: $T_{k,t}^{S,\mathrm{out}}={\tilde{T}}_{k,t}^{S,\mathrm{out}}\exp \left(-\frac{\lambda L_k}{{\mathrm{cm}}_{k,t}^S}\right),$ where T.sub.k,t.sup.S,out represents the actual temperature of the outlet water flow of the pipeline; T.sub.k,t.sup.S,out represents a weighted average of temperatures of injection water flows at previous time periods; λ represents a thermal conductivity coefficient of the pipeline; and L.sub.k represents a length of the pipeline k; determining an actual outlet carbon flow rate of the pipeline at the time period t: $R_{k,t}^{\mathrm{BHS},\mathrm{out}}={\tilde{R}}_{k,t}^{{\mathrm{BHS},\mathrm{out}}_{}}\frac{T_{k,t}^{S,\mathrm{out}}}{{\tilde{T}}_{k,t}^{S,\mathrm{out}}};$ and determining a network loss carbon flow rate of the pipeline at the time period t: $R_{k,t}^{\mathrm{BHS},\mathrm{Loss}}={\tilde{R}}_{k,t}^{\mathrm{BHS},\mathrm{Loss}}\left(1-\frac{T_{k,t}^{S,\mathrm{out}}}{{\tilde{T}}_{k,t}^{S,\mathrm{out}}}\right);$ determining the nodal carbon flow densities at the plurality of time periods: ${\rho }_{n,t}^{\mathrm{NHS}}=\frac{\underset{k\in {\Omega }_n^{BH+}}{.Math.}\left(R_{k,t}^{\mathrm{BHS},\mathrm{out}}+X{R}_{k,t}^{\mathrm{BHS},\mathrm{Loss}}\right)}{\underset{k\in {\Omega }_n^{BH+}}{.Math.}c{m}_{k,t}^S{T}_{k,t}^{S,\mathrm{out}}};$ determining, for the water return network, the nodal carbon flow densities at the plurality of time periods: ${\rho }_{n,t}^{\mathrm{NHS}}=\frac{\underset{k\in {\Omega }_n^{\mathrm{BH}-}}{.Math.}\left(R_{k,t}^{\mathrm{BHS},\mathrm{out}}+{\mathrm{XR}}_{k,t}^{\mathrm{BHS},\mathrm{Loss}}\right)}{\underset{k\in {\Omega }_n^{\mathrm{BH}-}}{.Math.}c{m}_{k,t}^R{T}_{k,t}^{R,\mathrm{out}}}.$

15. A computer-readable storage medium, having a computer program stored thereon, wherein the computer program, when executed by a processor, implements the method for measuring the carbon emission of the district heating system according to claim 1.

16. The computer-readable storage medium according to claim 15, wherein the method further comprises: constructing the steady carbon emission flow model of the district heating system based on the pipeline carbon flow rate, the network loss carbon emission, the nodal carbon flow density, the pipeline carbon flow density, the heat source carbon flow rate, and the heat load carbon flow rate of the district heating system, wherein said constructing the steady carbon emission flow model of the district heating system based on the pipeline carbon flow rate, the network loss carbon emission, the nodal carbon flow density, the pipeline carbon flow density, the heat source carbon flow rate, and the heat load carbon flow rate of the district heating system is: determining the pipeline carbon flow rate of the district heating system comprising a carbon flow rate of a water supply network and a carbon flow rate of a water return network, the carbon flow rate of the water supply network being:
R.sub.k.sup.BHS,in=ρ.sub.k.sup.BHScm.sub.k.sup.ST.sub.k.sup.S,in,∀k∈Ω.sup.BH
R.sub.k.sup.BHS,out=ρ.sub.k.sup.BHScm.sub.k.sup.ST.sub.k.sup.S,out,∀k∈Ω.sup.BH, where R.sub.k.sup.BHS,in and R.sub.k.sup.BHS,out represent an inlet carbon flow rate and an outlet carbon flow rate of a pipeline k in the water supply network, respectively, unit: tCO.sub.2/h; ρ.sub.k.sup.BHS represents a carbon flow density of the pipeline k in the water supply network, unit: tCO.sub.2/MWh; c represents a specific heat capacity of water, unit: MWh/(kg.Math.° C.); m.sub.k.sup.S represents a mass flow of the pipeline k in the water supply network, unit: kg/h; T.sub.k.sup.S,in and T.sub.k.sup.S,out represent an inlet temperature and an outlet temperature of the pipeline k in the water supply network, unit: ° C.; and Ω.sup.BH represents a set of pipelines in the district heating system; and the carbon flow rate of the water return network being:
R.sub.k.sup.BHR,in=ρ.sub.k.sup.BHRcm.sub.k.sup.RT.sub.k.sup.R,in,∀k∈Ω.sup.BH
R.sub.k.sup.BHR,out=ρ.sub.k.sup.BHRcm.sub.k.sup.RT.sub.k.sup.R,out,∀k∈Ω.sup.BH, where R.sub.k.sup.BHR,in and R.sub.k.sup.BHR,out represent an inlet carbon flow rate and an outlet carbon flow rate of a pipeline k in the water return network, respectively, unit: tCO.sub.2/h; ρ.sub.k.sup.BHR represents a carbon flow density of the pipeline k in the water return network, unit: tCO.sub.2/h; m.sub.k.sup.R represents a mass flow of the pipeline k in the water return network, unit: kg/h; and T.sub.k.sup.R,in and T.sub.k.sup.R,out represent an inlet temperature and an outlet temperature of the pipeline k in the water return network, unit: ° C.; determining the network loss carbon emission of the district heating system comprising a network loss carbon flow rate of the water supply network and a network loss carbon flow rate of the water return network, said determining the network loss carbon emission of the district heating system being: determining a temperature difference of a pipeline in the water supply network and a temperature difference of a pipeline in the water return network:
T.sub.k.sup.S,Loss=T.sub.k.sup.S,in−T.sub.k.sup.S,out
T.sub.k.sup.R,Loss=T.sub.k.sup.R,in−T.sub.k.sup.R,out where T.sub.k.sup.S,Loss and T.sub.k.sup.R,Loss represent a temperature difference between both ends of the pipeline k in the water supply network and a temperature difference between both ends of the pipeline k in the water return network, respectively, unit: ° C.; and determining the network loss carbon flow rate of the water supply network and the network loss carbon flow rate of the water return network:
R.sub.k.sup.BHS,Loss=ρ.sub.k.sup.BHScm.sub.k.sup.ST.sub.k.sup.S,Loss,∀k∈Ω.sup.BH
R.sub.k.sup.BHR,Loss=ρ.sub.k.sup.BHRcm.sub.k.sup.RT.sub.k.sup.R,Loss,∀k∈Ω.sup.BH, where R.sub.k.sup.BHS,Loss and R.sub.k.sup.BHR,Loss represent a network loss carbon flow rate of the pipeline k in the water supply network and a network loss carbon flow rate of the pipeline k in the water return network, respectively, unit: tCO.sub.2/h; determining the nodal carbon flow density of the district heating system, said determining the nodal carbon flow density of the district heating system comprising: determining, for each node in the district heating system, that conservation of mass and conservation of energy are satisfied at the node: $m_n^S=\underset{k\in {\Omega }_n^{BH+}}{.Math.}{m}_k^S,\forall n\in {\Omega }^{\mathrm{NH}}$ $T_n^S{m}_n^S=\underset{k\in {\Omega }_n^{BH+}}{.Math.}{T}_k^{S,out}{m}_k^S,\forall n\in {\Omega }^{\mathrm{NH}},$ where m.sub.n.sup.S represents a total mass flow flowing through a node n in the water supply network, unit: kg/h; Ω.sub.n.sup.BH+ represents a set of injection pipelines at the node n in the water supply network; Ω.sup.NH represents a set of nodes in the district heating system; and T.sub.n.sup.S represents a water flow temperature of the node n in the water supply network, unit: ° C.; determining, for each node in the district heating system, that conservation of carbon emission is satisfied at the node, a carbon flow rate of the node n being equal to a sum of outlet carbon flow rates of all the injection pipelines and network loss carbon flow rates allocated to the injection pipelines; $R_n^{\mathrm{NHS}}=\underset{k\in {\Omega }_n^{BH+}}{.Math.}\left(R_k^{\mathrm{BHS},out}+\lambda R_k^{\mathrm{BHS},\mathrm{Loss}}\right)=\underset{k\in {\Omega }_n^{BH+}}{.Math.}{\rho }_k^{\mathrm{BHS}}{\mathrm{cm}}_k^S\left(T_k^{S,out}+X{T}_k^{S,\mathrm{Loss}}\right),\forall n\in {\Omega }^{\mathrm{NH}},$ where R.sub.n.sup.NHS represents the carbon flow rate of the node n in the water supply network, unit: tCO.sub.2/h; and X represents an allocating coefficient of a network loss carbon flow rate of an injection pipeline; determining a carbon flow density of the node n in the water supply network: ${\rho }_n^{\mathrm{NHS}}=\frac{R_n^{\mathrm{NHS}}}{c{m}_n^S{T}_n^S}=\frac{\underset{k\in {\Omega }_n^{BH+}}{.Math.}{\rho }_k^{\mathrm{BHS}}{m}_k^S\left(T_k^{S,\mathrm{out}}+X{T}_k^{S,\mathrm{Loss}}\right)}{\underset{k\in {\Omega }_n^{BH+}}{.Math.}{m}_k^S{T}_k^{S,\mathrm{out}}},\forall n\in {\Omega }^{\mathrm{NH}},$ where ρ.sub.n.sup.NHS represents the carbon flow density of the node n in the water supply network, unit: tCO.sub.2/MWh; and determining a carbon flow density of the node n in the water return network: ${\rho }_n^{\mathrm{NHR}}=\frac{\underset{k\in {\Omega }_n^{\mathrm{BH}-}}{.Math.}{\rho }_k^{\mathrm{BHR}}{m}_k^R\left(T_k^{R,\mathrm{out}}+{\mathrm{XT}}_k^{R,\mathrm{Loss}}\right)}{\underset{k\in {\Omega }_n^{\mathrm{BH}-}}{.Math.}{m}_k^R{T}_k^{R,\mathrm{out}}},\forall n\in {\Omega }^{\mathrm{NH}},$ where ρ.sub.n.sup.NHR represents the carbon flow density of the node n in the water return network, unit: tCO.sub.2/MWh; and Ω.sub.n.sup.BH− represents a set of outflow pipelines at the node n in the water supply network; determining the pipeline carbon flow density of the district heating system:
ρ.sub.k.sup.BHS=ρ.sub.n.sup.NHS,n=Γ.sub.k.sup.NH+,∀k∈Ω.sup.BH
ρ.sub.k.sup.BHR=ρ.sub.n.sup.NHR,n=Γ.sub.k.sup.NH−,∀k∈Ω.sup.BH, where Γ.sub.k.sup.NH+ and Γ.sub.k.sup.NH+ represent an injection node and an outflow node of the pipeline k in the water supply network, respectively; determining the heat source carbon flow rate of the district heating system, said determining the heat source carbon flow rate of the district heating system comprising: determining a heat output of a heat source:
Q.sub.i=cm.sub.n.sup.S(T.sub.n.sup.S−T.sub.n.sup.R),n=Γ.sub.i.sup.NH,∀i∈Ω.sup.GH, where Q.sub.i represents a heat output of a heat source i, unit: MW; Γ.sub.i.sup.NH represents a node where the heat source i is located; and Ω.sup.GH represents a set of heat sources; and determining that conservation of carbon emission at a heat source node is satisfied:
ρ.sub.n.sup.NHScm.sub.n.sup.ST.sub.n.sup.S=ρ.sub.i.sup.GHQ.sub.i+ρ.sub.n.sup.NHRcm.sub.n.sup.ST.sub.n.sup.R,n=Γ.sub.i.sup.NH,∀i∈Ω.sup.GH, where ρ.sub.i.sup.GH represents a carbon flow density of the heat source i, unit: tCO.sub.2/MWh; and determining the heat load carbon flow rate of the district heating system, said determining the heat load carbon flow rate of the district heating system comprising: determining a heat load demand:
q.sub.j=cm.sub.n.sup.S(T.sub.n.sup.S−T.sub.n.sup.R),n=Γ.sub.j.sup.NH,∀j∈Ω.sup.LH, where q.sub.j represents a heat demand of a heat load j, unit: MW; Γ.sub.j.sup.NH represents a node where the heat load j is located; and Ω.sup.LH represents a set of heat loads; determining a carbon flow density of a heat load node:
ρ.sub.n.sup.NHR=ρ.sub.n.sup.NHS,n=Γ.sub.j.sup.NH,∀j∈Ω.sup.LH; and determining the heat load carbon flow rate:
R.sub.j.sup.LH=ρ.sub.n.sup.NHSq.sub.j=ρ.sub.n.sup.NHScm.sub.n.sup.S(T.sub.n.sup.S−T.sub.n.sup.R),n=Γ.sub.j.sup.NH,∀j∈Ω.sup.LH, where R.sub.j.sup.LH represents a carbon flow rate of the heat load j, unit: tCO.sub.2/h.

17. The computer-readable storage medium according to claim 16, wherein the method further comprises: calculating a matrix representation of the steady carbon emission flow model of the district heating system.

18. The computer-readable storage medium according to claim 17, wherein said calculating the matrix representation of the steady carbon emission flow model of the district heating system comprises: constructing a branch heat flow matrix, a branch network loss matrix, and a nodal heat flow flux matrix of the district heating system; and calculating a nodal carbon flow density vector of a heat network, and then calculating a branch carbon flow rate matrix, a network loss carbon flow rate matrix, and a load carbon flow rate vector of each of the water supply network and the water return network.

19. The computer-readable storage medium according to claim 18, wherein said calculating the matrix representation of the steady carbon emission flow model of the district heating system comprises: constructing the branch heat flow matrix of the district heating system comprising a branch heat flow matrix of the water supply network and a branch heat flow matrix of the water return network, wherein said constructing the branch heat flow matrix of the district heating system comprises: determining elements of the branch heat flow matrix of the water supply network:
Q.sub.B,ij.sup.S=cm.sub.k.sup.ST.sub.k.sup.S,out,Q.sub.B,ji.sup.S=0, where Q.sub.B,ij.sup.S and Q.sub.B,ji.sup.S represent elements in a branch heat flow matrix Q.sub.B.sup.S of the water supply network; and determining elements of the branch heat flow matrix of the water return network:
Q.sub.B,ij.sup.R=cm.sub.k.sup.RT.sub.k.sup.R,out,Q.sub.B,ji.sup.R=0, where Q.sub.B,ij.sup.R and Q.sub.B,ji.sup.R represent elements in a branch heat flow matrix Q.sub.B.sup.R of the water return network; constructing the branch network loss matrix of the district heating system comprising a branch network loss matrix of the water supply network and a branch network loss matrix of the water return network, wherein said constructing the branch network loss matrix of the district heating system comprises: determining elements of the branch network loss matrix of the water supply network:
Q.sub.BL,ij.sup.S=cm.sub.k.sup.ST.sub.k.sup.S,Loss,Q.sub.BL,ji.sup.S=0, where Q.sub.BL,ij.sup.S and Q.sub.BL,ji.sup.S represent elements in a branch network loss matrix Q.sub.BL.sup.S of the water supply network; and determining elements of the branch network loss matrix of the water return network:
Q.sub.BL,ij.sup.R=cm.sub.k.sup.RT.sub.k.sup.R,Loss,Q.sub.BL,ji.sup.R=0, where Q.sub.BL,ij.sup.R and Q.sub.BL,ji.sup.R represent elements in a branch network loss matrix Q.sub.BL.sup.R of the water return network; constructing the nodal heat flow flux matrix of the district heating system comprising a nodal heat flow flux matrix of the water supply network and a nodal heat flow flux matrix of the water return network, wherein said constructing the nodal heat flow flux matrix of the district heating system comprises: determining, for a node in no connection to the heat source, the nodal heat flow flux matrix of the water supply network:
Q.sub.N.sup.S=diag{ζ.sub.N.sub.NHQ.sub.B.sup.S}, where Q.sub.N.sup.S represents the nodal heat flow flux matrix of the water supply network; ζ.sub.N.sub.NH represents a coefficient matrix of a branch heat flow; and N.sub.NH represents a number of nodes of the district heating system; determining, for a node in no connection to the heat source, the nodal heat flow flux matrix of the water return network:
Q.sub.N.sup.R=diag{ζ.sub.N.sub.NHQ.sub.B.sup.R}, where Q.sub.N.sup.R represents the nodal heat flow flux matrix of the water return network; determining, for a node connected to the heat source, a nodal integrated energy flow flux matrix:
{circumflex over (Q)}.sub.N.sup.S=diag{ζ.sub.N.sub.NHQ.sub.B.sup.S+ζ.sub.N.sub.GHQ.sub.G}, where {circumflex over (Q)}.sub.N.sup.S represents the nodal integrated energy flow flux matrix; ζ.sub.N.sub.GH represents a coefficient matrix of a heat flow injected by the heat source; and N.sub.GH represents a number of heat sources of the district heating system; determining, for all nodes, that a total injected carbon emission of the nodes is equal to a sum of injected carbon flow rates of all branches connected to the nodes:
Q.sub.N.sup.Sρ.sup.NHS=(Q.sub.B.sup.S+XQ.sub.BL.sup.S).sup.Tρ.sup.NHS
Q.sub.N.sup.Rρ.sup.NHR=(Q.sub.B.sup.R+XQ.sub.BL.sup.R).sup.Tρ.sup.NHR, where ρ.sup.NHS represents a matrix formed by the carbon flow density ρ.sub.n.sup.NHS of the node n in the water supply network; and ρ.sup.NHR represents a matrix formed by the carbon flow density ρ.sub.n.sup.NHR of the node n in the water return network; determining, for the heat load node, that the heat load node has an equal carbon flow density in the water supply network and the water return network:
Bρ.sup.NHS=Bρ.sup.NHR, where B represents a heat load-node association matrix, when the heat load j is connected to the node n, B.sub.jn=1, otherwise B.sub.jn=0; determining, for the heat source node, a matrix relation of the heat source node based on a conservation of carbon emission as:
C{circumflex over (Q)}.sub.N.sup.Sρ.sup.NHS=CQ.sub.N.sup.Rρ.sup.NHR+Q.sub.G.sup.Tσ.sup.GH, where C represents a 0-1 matrix associated with the heat source node, when the node n is connected to the heat source, C.sub.nn=1, otherwise C.sub.nn=0; determining the nodal carbon flow density vector of the heat network as: $\left[\begin{array}{c}{\rho }^{\mathrm{NHS}}\\ {\rho }^{\mathrm{NHR}}\end{array}\right]={\left[\begin{array}{ccc}{Q}_N^S-(Q_B^S& {+X{Q}_{BL}^S)}^T& 0\\ & 0& Q_N^R-{\left(Q_B^R+X{Q}_{BL}^R\right)}^T\\ & B& -B\\ & C{\hat{Q}}_N^S& -C{Q}_N^R\end{array}\right]}^{-1}\left[\begin{array}{c}0\\ 0\\ 0\\ Q_G^T{\rho }^{GH}\end{array}\right];$ and calculating the branch carbon flow rate matrix, the network loss carbon flow rate matrix, and the load carbon flow rate vector of each of the water supply network and the water return network.

20. The computer-readable storage medium according to claim 19, wherein the method further comprises: constructing the dynamic carbon emission flow model of the district heating system based on the water element carbon flow rates at the plurality of time periods, the actual outlet carbon flow rates of the pipeline at the plurality of time periods, the network loss carbon flow rates at the plurality of time periods, and the nodal carbon flow densities at the plurality of time periods, wherein said constructing the dynamic carbon emission flow model of the district heating system based on the water element carbon flow rates at the plurality of time periods, the actual outlet carbon flow rates of the pipeline at the plurality of time periods, the network loss carbon flow rates at the plurality of time periods, and the nodal carbon flow densities at the plurality of time periods is: determining the water element carbon flow rates at the plurality of time periods: ${\tilde{R}}_{k,t}^{\mathrm{BHS},\mathrm{out}}=\frac{R_{k,t-{\delta }_{k,t}}^{\mathrm{BHS},\mathrm{in}}\left(B_{k,t}-\sigma M_k\right)}{\left(m_{k,t-{\delta }_{k,t}}^S\Delta t\right)}+\underset{\tau =l-{\phi }_{k,t}+1}{\overset{t-{\delta }_{k,t}-1}{.Math.}}{R}_{k,\tau }^{\mathrm{BHS},\mathrm{in}}+\frac{R_{k,t-{\delta }_{k,t}}^{\mathrm{BHS},\mathrm{in}}\left(m_{k,t}^S\Delta t+\sigma M_k-A_{k,t}\right)}{\left(m_{k,t-{\phi }_{k,t}}^S\Delta t\right)},$ where {tilde over (R)}.sub.k,t.sup.BHS,out represents a water element carbon flow rate at a time period t; σ represents a density of water; M.sub.k represents a volume of the pipeline k; φ.sub.k,t represents an injection time period of an earliest water element component contained in a water flow flowing out of the pipeline k at the time period t; δ.sub.k,t represents an injection time period of a latest water element component contained in the water flow flowing out of the pipeline k at the time period t; A.sub.k,t represents a total injection water flow amount for φ.sub.k,t time periods before the time period t; B.sub.k,t represents a total injection water flow amount from a time period t−δ.sub.k,t to the time period t; and expressions for A.sub.k,t and B.sub.k,t are: $A_{k,t}=\{\begin{array}{cc}\underset{\tau =t-{\phi }_{k,1}+1}{\overset{t}{.Math.}}{m}_{k,\tau }^S\Delta t,& {\phi }_{k,t}\ge {\delta }_{k,t}+1\\ B_{k,t},& {\phi }_{k,t}<{\delta }_{k,t}+1\end{array},$ $B_{k,t}=\underset{\tau =t-{\delta }_{k,t}}{\overset{t}{.Math.}}{m}_{k,\tau }^S\Delta t;$ determining the actual outlet carbon flow rates of the pipeline at the plurality of time periods and the network loss carbon flow rates at the plurality of time periods, said determining the actual outlet carbon flow rates of the pipeline at the plurality of time periods and the network loss carbon flow rates at the plurality of time periods comprising: calculating an actual temperature of an outlet water flow of the pipeline based on a conveying loss of the pipeline: $T_{k,t}^{S,\mathrm{out}}={\tilde{T}}_{k,t}^{S,\mathrm{out}}\exp \left(-\frac{\lambda L_k}{c{m}_{k,t}^S}\right),$ where T.sub.k,t.sup.S,out represents the actual temperature of the outlet water flow of the pipeline; {acute over (T)}.sub.k,t.sup.S,out represents a weighted average of temperatures of injection water flows at previous time periods; λ represents a thermal conductivity coefficient of the pipeline; and L.sub.k represents a length of the pipeline k; determining an actual outlet carbon flow rate of the pipeline at the time period t: $R_{k,t}^{\mathrm{BHS},\mathrm{out}}={\tilde{R}}_{k,t}^{\mathrm{BHS},\mathrm{out}}\frac{T_{k,t}^{S,\mathrm{out}}}{{\tilde{T}}_{k,t}^{S,\mathrm{out}}};$ and determining a network loss carbon flow rate of the pipeline at the time period t: $R_{k,t}^{\mathrm{BHS},\mathrm{Loss}}={\tilde{R}}_{k,t}^{\mathrm{BHS},\mathrm{Loss}}\left(1-\frac{T_{k,t}^{S,\mathrm{out}}}{{\tilde{T}}_{k,t}^{S,\mathrm{out}}}\right);$ determining the nodal carbon flow densities at the plurality of time periods: ${\rho }_{n,t}^{\mathrm{NHS}}=\frac{\underset{k\in {\Omega }_n^{BH+}}{.Math.}\left(R_{k,t}^{\mathrm{BHS},\mathrm{out}}+X{R}_{k,t}^{\mathrm{BHS},\mathrm{Loss}}\right)}{\underset{k\in {\Omega }_n^{BH+}}{.Math.}c{m}_{k,t}^S{T}_{k,t}^{S,\mathrm{out}}};$ and determining, for the water return network, the nodal carbon flow densities at the plurality of time periods: ${\rho }_{n,t}^{\mathrm{NHS}}=\frac{\underset{k\in {\Omega }_n^{\mathrm{BH}-}}{.Math.}\left(R_{k,t}^{\mathrm{BHS},\mathrm{out}}+{\mathrm{XR}}_{k,t}^{\mathrm{BHS},\mathrm{Loss}}\right)}{\underset{k\in {\Omega }_n^{\mathrm{BH}-}}{.Math.}c{m}_{k,t}^R{T}_{k,t}^{R,\mathrm{out}}}.$

Description

BRIEF DESCRIPTION OF THE DRAWINGS

[0090] The above and/or additional aspects and advantages of the present disclosure will become more apparent and more understandable from the following description of embodiments taken in conjunction with the accompanying drawings.

[0091] FIG. 1 is a flowchart illustrating a method for measuring carbon emission of a district heating system according to an embodiment of the present disclosure.

[0092] FIG. 2 is a schematic diagram showing a framework structure of a method for measuring carbon emission of a district heating system according to an embodiment of the present disclosure.

[0093] FIG. 3 is a schematic diagram showing a result of solving a steady carbon emission flow of a district heating system according to an embodiment of the present disclosure.

[0094] FIG. 4 is a schematic diagram showing a result of solving a dynamic carbon emission flow of a district heating system according to an embodiment of the present disclosure.

[0095] FIG. 5 is an example diagram of an apparatus for measuring carbon emission of a district heating system according to an embodiment of the present disclosure.

[0096] FIG. 6 is a schematic diagram showing a structure of an electronic device according to an embodiment of the present disclosure.

DETAILED DESCRIPTION

[0097] The embodiments of the present disclosure will be described in detail below with reference to examples thereof as illustrated in the accompanying drawings, throughout which same or similar elements, or elements having same or similar functions, are denoted by same or similar reference numerals. The embodiments described below with reference to the drawings are illustrative only, and are intended to explain, rather than limiting, the present disclosure.

[0098] FIG. 1 is a flowchart illustrating a method for measuring carbon emission of a district heating system according to an embodiment of the present disclosure.

[0099] As illustrated in FIG. 1, the method for measuring the carbon emission of the district heating system includes actions at blocks S101 to S103.

[0100] At block S101, a steady carbon emission amount of a current district heating system is obtained using a pre-trained steady carbon emission flow model. The pre-trained steady carbon emission flow model is constructed based on a pipeline carbon flow rate, network loss carbon emission, a nodal carbon flow density, a pipeline carbon flow density, a heat source carbon flow rate, and a heat load carbon flow rate of the current district heating system.

[0101] It should be understood that, according to the embodiments of the present disclosure, the steady carbon emission flow model of the district heating system is constructed based on the carbon emission flow theory. The constructed steady carbon emission flow model of the district heating system is used to obtain the steady carbon emission amount of the current district heating system.

[0102] In some embodiments, the steady carbon emission flow model of the district heating system is formed by the pipeline carbon flow rate, the network loss carbon emission, the nodal carbon flow density, the pipeline carbon flow density, the heat source carbon flow rate, and the heat load carbon flow rate of the district heating system. The operation of forming the steady carbon emission flow model of the district heating system includes the following actions.

[0103] 1-1) In a case of a steady carbon emission flow, no convergence or divergence of a water flow is found in a single pipeline since a conveying delay is not considered. Therefore, carbon flow densities at a pipeline inlet and a pipeline outlet are equal and uniformly defined as a carbon flow density of the pipeline. An inlet carbon flow rate of the pipeline is equal to a product of an inlet energy flow and the pipeline carbon flow density. An outlet carbon flow rate of the pipeline is equal to a product of an outlet energy flow and the pipeline carbon flow density. The pipeline carbon flow rate of the district heating system is determined. The pipeline carbon flow rate of the district heating system includes a carbon flow rate of a water supply network and a carbon flow rate of a water return network. The operation of determining the pipeline carbon flow rate of the district heating system includes the following actions at 1-1-1) and 1-1-2).

[0104] 1-1-1) The carbon flow rate of the water supply network is determined. An expression of the carbon flow rate of the water supply network is:

R.sub.k.sup.BHS,in=ρ.sub.k.sup.BHScm.sub.k.sup.ST.sub.k.sup.S,in,∀k∈Ω.sup.BH

R.sub.k.sup.BHS,out=ρ.sub.k.sup.BHScm.sub.k.sup.ST.sub.k.sup.S,out,∀k∈Ω.sup.BH  (1)

where R.sub.k.sup.BHS,in and R.sub.k.sup.BHS,out represent an inlet carbon flow rate and an outlet carbon flow rate of a pipeline k in the water supply network, respectively, unit: tCO.sub.2/h; ρ.sub.k.sup.BHS represents a carbon flow density of the pipeline k in the water supply network, unit: tCO.sub.2/MWh; c represents a specific heat capacity of water, unit: MWh/(kg.Math.° C.); m.sub.k.sup.S represents a mass flow of the pipeline k in the water supply network, unit: kg/h; T.sub.k.sup.S,in and T.sub.k.sup.S,out represent an inlet temperature and an outlet temperature of the pipeline k in the water supply network, unit: ° C.; and Ω.sup.BH represents a set of pipelines in the district heating system.

[0105] 1-1-2) The carbon flow rate of the water return network is determined as follows:

R.sub.k.sup.BHR,in=ρ.sub.k.sup.BHRcm.sub.k.sup.RT.sub.k.sup.R,in,∀k∈Ω.sup.BH

R.sub.k.sup.BHR,out=ρ.sub.k.sup.BHRcm.sub.k.sup.RT.sub.k.sup.R,out,∀k∈Ω.sup.BH  (2),

where R.sub.k.sup.BHR,in and R.sub.k.sup.BHR,out represent an inlet carbon flow rate and an outlet carbon flow rate of a pipeline k in the water return network, respectively, unit: tCO.sub.2/h; ρ.sub.k.sup.BHR represents a carbon flow density of the pipeline k in the water return network, unit: tCO.sub.2/h; m.sub.k.sup.R represents a mass flow of the pipeline k in the water return network, unit: kg/h; and T.sub.k.sup.R,in and T.sub.k.sup.R,out represent an inlet temperature and an outlet temperature of the pipeline k in the water return network, unit: ° C.

[0106] 1-2) A heat network bears a large conveying loss, which causes additional carbon emission. The network loss carbon emission of the district heating system is determined. The network loss carbon emission of the district heating system includes a network loss carbon flow rate of the water supply network and a network loss carbon flow rate of the water return network. The operation of determining the network loss carbon emission of the district heating system includes the following actions at 1-2-1) and 1-2-2).

[0107] 1-2-1) A loss in the heat network is mainly resulted from a temperature loss of a water flow in a pipeline. Therefore, the network loss carbon emission mainly depends on a temperature difference between two ends of the pipeline. A temperature difference of a pipeline in the water supply network and a temperature difference of a pipeline in the water return network are determined as follows:

T.sub.k.sup.S,Loss=T.sub.k.sup.S,in−T.sub.k.sup.S,out

T.sub.k.sup.R,Loss=T.sub.k.sup.R,in−T.sub.k.sup.R,out  (3),

where T.sub.k.sup.S,Loss and T.sub.k.sup.R,Loss represent a temperature difference between both ends of the pipeline k in the water supply network and a temperature difference between both ends of the pipeline k in the water return network, respectively, unit: ° C.

[0108] 1-2-2) The network loss carbon flow rate of the water supply network and the network loss carbon flow rate of the water return network are determined as follows:

R.sub.k.sup.BHS,Loss=ρ.sub.k.sup.BHScm.sub.k.sup.ST.sub.k.sup.S,Loss,∀k∈Ω.sup.BH

R.sub.k.sup.BHR,Loss=ρ.sub.k.sup.BHRcm.sub.k.sup.RT.sub.k.sup.R,Loss,∀k∈Ω.sup.BH  (4),

where R.sub.k.sup.BHS,Loss and R.sub.k.sup.BHR,Loss represent a network loss carbon flow rate of the pipeline k in the water supply network and a network loss carbon flow rate of the pipeline k in the water return network, respectively, unit: tCO.sub.2/h.

[0109] 1-3) For each node in the heat network, convergence of water flows from several pipelines may exist. Conservation of mass and conservation of energy are satisfied at the node. The nodal carbon flow density of the district heating system is determined. The operation of determining the nodal carbon flow density of the district heating system includes the following actions at 1-3-1), 1-3-2), 1-3-3) and 1-3-4).

[0110] 1-3-1) For each node in the district heating system, it is determined that conservation of mass and conservation of energy are satisfied at the node:

[00013]$\begin{array}{cc}\begin{array}{c}{m}_n^S=\underset{k\in {\Omega }_n^{BH+}}{.Math.}{m}_k^S,\forall n\in {\Omega }^{NH}\\ T_n^S{m}_n^S=\underset{k\in {\Omega }_n^{BH+}}{.Math.}{T}_k^{S,out}{m}_k^S,\forall n\in {\Omega }^{NH}\end{array},& \left(5\right)\end{array}$

where m.sub.n.sup.S represents a total mass flow flowing through a node n in the water supply network, unit: kg/h; Ω.sub.n.sup.BH+ represents a set of injection pipelines at the node n in the water supply network; Ω.sup.NH represents a set of nodes in the district heating system; T.sub.n.sup.S represents a water flow temperature of the node n in the water supply network, unit: ° C.

[0111] 1-3-2) For each node in the district heating system, conservation of carbon emission is satisfied at the node. In this way, a carbon flow rate of the node n is equal to a sum of outlet carbon flow rates of all the injection pipelines and network loss carbon flow rates allocated to the injection pipelines, which is expressed as:

[00014]$\begin{array}{cc}{R}_n^{NHS}=\underset{k\in {\Omega }_n^{BH+}}{.Math.}\left(R_k^{BHS,out}+X{R}_k^{\mathrm{BHS},\mathrm{Loss}}\right)=\underset{k\in {\Omega }_n^{BH+}}{.Math.}{\rho }_k^{\mathrm{BHS}}{\mathrm{cm}}_k^S\left(T_k^{S,out}+X{T}_k^{S,\mathrm{Loss}}\right),\forall n\in {\Omega }^{NH},& \left(6\right)\end{array}$

where R.sub.n.sup.NHS represents the carbon flow rate of the node n in the water supply network, unit: tCO.sub.2/h; and X represents an allocating coefficient of a network loss carbon flow rate of an injection pipeline.

[0112] 1-3-3) A carbon flow density of the node n in the water supply network is determined as follows:

[00015]$\begin{array}{cc}{\rho }_n^{NHS}=\frac{R_n^{\mathrm{NHS}}}{c{m}_n^S{T}_n^S}=\frac{\underset{k\in {\Omega }_n^{BH+}}{.Math.}{\rho }_k^{\mathrm{BHS}}{m}_k^S\left(T_k^{S,\mathrm{out}}+X{T}_k^{S,\mathrm{Loss}}\right)}{\underset{k\in {\Omega }_n^{BH+}}{.Math.}{m}_k^S{T}_k^{S,\mathrm{out}}},\forall n\in {\Omega }^{NH}& \left(7\right)\end{array}$

where ρ.sub.n.sup.NHS represents the carbon flow density of the node n in the water supply network, unit: tCO.sub.2/MWh.

[0113] 1-3-4) A carbon flow density of the node n in the water return network is determined as follows:

[00016]$\begin{array}{cc}{\rho }_n^{NHR}=\frac{\underset{k\in {\Omega }_n^{\mathrm{BH}-}}{.Math.}{\rho }_k^{\mathrm{BHR}}{m}_k^R\left(T_k^{R,\mathrm{out}}+{\mathrm{XT}}_k^{R,\mathrm{Loss}}\right)}{\underset{k\in {\Omega }_n^{\mathrm{BH}-}}{.Math.}{m}_k^R{T}_k^{R,\mathrm{out}}},\forall n\in {\Omega }^{NH},& \left(8\right)\end{array}$

where ρ.sub.n.sup.NHR represents the carbon flow density of the node n in the water return network, unit: tCO.sub.2/MWh; and Ω.sup.BH− represents a set of outflow pipelines at the node n in the water supply network.

[0114] 1-4) Each pipeline in the water supply network and the water return network satisfies an energy allocation criterion. That is, the pipeline carbon flow density is equal to a carbon flow density of an injection node of a pipeline. The pipeline carbon flow density of the district heating system is determined as follows:

ρ.sub.k.sup.BHS=ρ.sub.n.sup.NHS,n=Γ.sub.k.sup.NH+,∀k∈ΩBH

ρ.sub.k.sup.BHR=ρ.sub.n.sup.NHR,n=θ.sub.k.sup.NH−,∀k∈Ω.sup.BH  (9),

where Γ.sub.k.sup.NH+ and Γ.sub.k.sup.NH+ represent an injection node and an outflow node of the pipeline k in the water supply network, respectively.

[0115] 1-5) The heat source and the heat load are connection points of the water supply network and the water return network, through which a carbon emission flow of the water supply network and a carbon emission flow of the water return network can be connected to form a complete carbon emission flow of the heat network. The heat source carbon flow rate of the district heating system is determined. The operation of determining the heat source carbon flow rate of the district heating system includes the following actions at 1-5-1) and 1-5-2).

[0116] 1-5-1) A heat output of a heat source is determined as follows:

Q.sub.i=cm.sub.n.sup.S(T.sub.n.sup.S=T.sub.n.sup.R),n=Γ.sub.i.sup.NH,∀i∈Ω.sup.GH  (10),

where Q.sub.i represents a heat output of a heat source i, unit: MW; Γ.sub.i.sup.NH represents a node where the heat source i is located; and Ω.sup.GH represents a set of heat sources.

[0117] 1-5-2) At a heat source node, conservation of carbon emission is also satisfied. That is, a carbon flow rate of water supply at the heat source node is equal to a sum of an injected carbon flow rate of the heat source and a carbon flow rate of water return at the heat source node. It is equivalent to that the heat source injects energy into the heat network while injecting corresponding carbon emission. A ratio between the injected energy and the injected carbon emission is the carbon flow density of the heat source. In some embodiments,

ρ.sub.n.sup.NHScm.sub.n.sup.ST.sub.n.sup.S=ρ.sub.i.sup.GHQ.sub.i+ρ.sub.n.sup.NHRcm.sub.n.sup.ST.sub.n.sup.R,n=Γ.sub.i.sup.NH,∀i∈Ω.sup.GH  (11),

where ρ.sub.i.sup.GH represents a carbon flow density of the heat source i, unit: tCO.sub.2/MWh.

[0118] 1-6) The carbon flow rate of the heat load depends on a heat load size and a carbon flow density of a heat load node. The heat load size may be obtained from a mass flow of the heat load node and supplied and returned water temperatures of the heat load node. The operation of determining the heat load carbon flow rate includes the following actions at 1-6-1), 1-6-2) and 1-6-3).

[0119] 1-6-1) A heat load demand is determined as follows:

q.sub.j=cm.sub.n.sup.S(T.sub.n.sup.S−T.sub.n.sup.R),n=Γ.sub.j.sup.NH,∀j∈Ω.sup.LH  (12),

where q.sub.j represents a heat demand of a heat load j, unit: MW; Γ.sub.j.sup.NH represents a node where the heat load j is located; and Ω.sup.LH represents a set of heat loads.

[0120] 1-6-2) Since there is no external energy injection at the heat load node, a node connected to the heat load have equal carbon flow densities in the water supply network and the water return network. A carbon flow density of a heat load node is determined as follows:

ρ.sub.n.sup.NHR=ρ.sub.n.sup.NHS,n=Γ.sub.j.sup.NH,∀j∈Ω.sup.LH  (13).

[0121] 1-6-3) The heat load carbon flow rate may be determined as follows:

R.sub.j.sup.LH=ρ.sub.n.sup.NHSq.sub.j=ρ.sub.n.sup.NHScm.sub.n.sup.S(T.sub.n.sup.S−T.sub.n.sup.R),n=Γ.sub.j.sup.NH,∀j∈Ω.sup.LH  (14),

where R.sub.j.sup.LH represents a carbon flow rate of the heat load j, unit: tCO.sub.2/h.

[0122] Optionally, in an embodiment of the present disclosure, the method for measuring the carbon emission of the district heating system further includes: calculating a matrix representation of the steady carbon emission flow model of the district heating system.

[0123] In an embodiment of the present disclosure, the operation of calculating the matrix representation of the steady carbon emission flow model of the district heating system includes: constructing a branch heat flow matrix, a branch network loss matrix, and a nodal heat flow flux matrix of the district heating system; and calculating a nodal carbon flow density vector of a heat network, and then calculating a branch carbon flow rate matrix, a network loss carbon flow rate matrix, and a load carbon flow rate vector of each of the water supply network and the water return network.

[0124] In some embodiments, an operation of constructing a matrix representation expression of the steady carbon emission flow model of the district heating system includes: constructing the branch heat flow matrix, the branch network loss matrix, and the nodal heat flow flux matrix of the district heating system. The operation of constructing the branch heat flow matrix, the branch network loss matrix, and the nodal heat flow flux matrix of the district heating system includes the following actions at 2-1), 2-2) and 2-3).

[0125] 2-1) Elements of a branch energy flow matrix are related to a water flow and a temperature in a pipeline. The branch heat flow matrix of the district heating system is constructed. The branch heat flow matrix of the district heating system includes a branch heat flow matrix of the water supply network and a branch heat flow matrix of the water return network. The operation of constructing the branch heat flow matrix of the district heating system includes the following actions at 2-1-1) and 2-1-2).

[0126] 2-1-1) Elements of the branch heat flow matrix of the water supply network are determined. Expressions of the elements are:

Q.sub.B,ij.sup.S=cm.sub.k.sup.ST.sub.k.sup.S,out,Q.sub.B,ji.sup.S=0  (15),

where Q.sub.B,ij.sup.S and Q.sub.B,ji.sup.S represent elements in a branch heat flow matrix Q.sub.B.sup.S of the water supply network.

[0127] 2-1-2) Elements of the branch heat flow matrix of the water return network are determined. Expressions of the elements are:

Q.sub.B,ij.sup.R=cm.sub.k.sup.RT.sub.k.sup.R,out,Q.sub.B,ji.sup.R=0  (16),

where Q.sub.B,ij.sup.R and Q.sub.B,ji.sup.R represent elements in a branch heat flow matrix Q.sub.B.sup.R of the water return network.

[0128] 2-2) Elements of the branch network loss matrix are related to a water flow and a temperature in a pipeline. The branch network loss matrix of the district heating system is constructed. The branch network loss matrix of the district heating system includes a branch network loss matrix of the water supply network and a branch network loss matrix of the water return network. The operation of constructing the branch network loss matrix of the district heating system includes the following actions at 2-2-1) and 2-2-2).

[0129] 2-2-1) Elements of the branch network loss matrix of the water supply network are determined. Expressions of the elements are:

Q.sub.BL,ij.sup.S=cm.sub.k.sup.ST.sub.k.sup.S,Loss,Q.sub.BL,ji.sup.S=0  (17),

where Q.sub.BL,ij.sup.S and Q.sub.BL,ji.sup.S represent elements in a branch network loss matrix Q.sub.BL.sup.S of the water supply network.

[0130] 2-2-2) Elements of the branch network loss matrix of the water return network are determined. Expressions of the elements are:

Q.sub.BL,ij.sup.R=cm.sub.k.sup.RT.sub.k.sup.R,Loss,Q.sub.BL,ji.sup.R=0  (18),

where Q.sub.BL,ij.sup.R and Q.sub.BL,ji.sup.R represent elements in a branch network loss matrix Q.sub.BL.sup.R of the water return network.

[0131] 2-3) A nodal energy flow flux matrix of the heat network is related to a branch energy flow matrix and a heat source injection matrix. The nodal heat flow flux matrix of the district heating system is constructed. The nodal heat flow flux matrix of the district heating system includes a nodal heat flow flux matrix of the water supply network and a nodal heat flow flux matrix of the water return network. The operation of constructing the nodal heat flow flux matrix of the district heating system includes the following actions at 2-3-1), 2-3-2) and 2-3-3).

[0132] 2-3-1) For a node in no connection to the heat source, an energy flow at the node is equal to a sum of energy flows of respective injection branches. The nodal heat flow flux matrix of the water supply network is determined. An expression of the nodal heat flow flux matrix of the water supply network is as follows:

Q.sub.N.sup.S=diag{ζ.sub.N.sub.NHQ.sub.B.sup.S}  (19),

where Q.sub.N.sup.S represents the nodal heat flow flux matrix of the water supply network; ζ.sub.N.sub.NH represents a coefficient matrix of a branch heat flow; and N.sub.NH represents a number of nodes of the district heating system.

[0133] 2-3-2) For a node in no connection to the heat source, an energy flow at the node is equal to a sum of energy flows of respective injection branches. The nodal heat flow flux matrix of the water return network is determined. An expression of the nodal heat flow flux matrix of the water return network is as follows:

Q.sub.N.sup.R=diag{ζ.sub.N.sub.NHQ.sub.B.sup.R}  (20),

where Q.sub.N.sup.R represents the nodal heat flow flux matrix of the water return network.

[0134] 2-3-3) For a node connected to the heat source, an energy flow flowing through the node is also related to an injection of the heat source. A nodal integrated energy flow flux matrix is determined. An expression of the nodal integrated energy flow flux matrix is as follows:

{circumflex over (Q)}.sub.N.sup.S=diag{ζ.sub.N.sub.NHQ.sub.B.sup.S+ζ.sub.N.sub.GHQ.sub.G}  (21),

where {circumflex over (Q)}.sub.N.sup.S represents the nodal integrated energy flow flux matrix; ζ.sub.N.sub.GH represents a coefficient matrix of a heat flow injected by the heat source; and N.sub.GH represents a number of heat sources of the district heating system.

[0135] 2-4) For all nodes, a total injected carbon emission of the nodes is determined to be equal to a sum of injected carbon flow rates of all branches connected to the nodes. In some embodiments,

Q.sub.N.sup.Sρ.sup.NHS=(Q.sub.B.sup.S+XQ.sub.BL.sup.S).sup.Tρ.sup.NHS

Q.sub.N.sup.Rρ.sup.NHR=(Q.sub.B.sup.R+XQ.sub.BL.sup.R).sup.Tρ.sup.NHR,  (22),

where ρ.sup.NHS represents a matrix formed by the carbon flow density ρ.sub.n.sup.NHS of the node n in the water supply network determined at the action of 1-3-3); and ρ.sup.NHR represents a matrix formed by the carbon flow density ρ.sup.NHR of the node n in the water return network determined at the action of 1-3-4).

[0136] 2-4-1) For the heat load node, the heat load node is determined to have an equal carbon flow density in the water supply network and the water return network, an expression of which is as follows:

Bρ.sup.NHS=Bρ.sup.NHR  (23),

where B represents a heat load-node association matrix, when the heat load j is connected to the node n, B.sub.jn=1, otherwise B.sub.jn=0.

[0137] 2-4-2) For the heat source node, a matrix relation of the heat source node is determined based on the conservation of carbon emission obtained at the action of 1-5-2) as:

C{circumflex over (Q)}.sub.N.sup.Sρ.sup.NHS=CQ.sub.N.sup.Rρ.sup.NHR+Q.sub.G.sup.Tρ.sup.GH  (24),

where C represents a 0-1 matrix associated with the heat source node, when the node n is connected to the heat source, C.sub.nn=1, otherwise C.sub.nn=0.

[0138] 2-5) The action of 2-4) is rearranged to obtain a matrix expression as:

[00017]$\begin{array}{cc}\left[\begin{array}{cc}{Q}_N^S-{\left(Q_B^S+X{Q}_{BL}^S\right)}^T& 0\\ 0& Q_N^R-{\left(Q_B^R+{\mathrm{XQ}}_{BL}^R\right)}^T\\ B& -B\\ C{\hat{Q}}_N^S& -C{Q}_N^R\end{array}\right]\left[\begin{array}{c}{\rho }^{\mathrm{NHS}}\\ {\rho }^{NHR}\end{array}\right]=\left[\begin{array}{c}0\\ 0\\ 0\\ Q_G^T{\rho }^{GH}\end{array}\right].& \left(25\right)\end{array}$

[0139] Therefore, the nodal carbon flow density vector of the heat network may be obtained. An expression of the nodal carbon flow density vector of the heat network is:

[00018]$\begin{array}{cc}\left[\begin{array}{c}{\rho }^{\mathrm{NHS}}\\ {\rho }^{NHR}\end{array}\right]={\left[\begin{array}{cc}{Q}_N^S-{\left(Q_B^S+X{Q}_{BL}^S\right)}^T& 0\\ 0& Q_N^R-{\left(Q_B^R+{\mathrm{XQ}}_{BL}^R\right)}^T\\ B& -B\\ C{\hat{Q}}_N^S& -C{Q}_N^R\end{array}\right]}^{-1}\left[\begin{array}{c}0\\ 0\\ 0\\ Q_G^T{\rho }^{GH}\end{array}\right].& \left(26\right)\end{array}$

[0140] Therefore, after the nodal carbon flow density vector of the heat network is obtained, the branch carbon flow rate matrix, the network loss carbon flow rate matrix, and the load carbon flow rate vector of each of the water supply network and the water return network may be calculated.

[0141] At block S102, a dynamic carbon emission amount of the current district heating system is obtained using a pre-trained dynamic carbon emission flow model. The pre-trained dynamic carbon emission flow model is constructed based on water element carbon flow rates at a plurality of time periods, actual outlet carbon flow rates of a pipeline at the plurality of time periods, network loss carbon flow rates at the plurality of time periods, and nodal carbon flow densities at the plurality of time periods of the current district heating system.

[0142] Due to a conveying delay in the heat network, a water flow flowing out of a pipeline in the heat network at a time period is different from a water flow injecting into the pipeline at the time period, which leads to a difference between the inlet carbon flow rate and an inlet carbon flow density, and the outlet carbon flow rate and an outlet carbon flow density of the pipeline in the heat network. Cell characterization is used to describe characteristics of dynamic conveying of the district heating system. The water flow flowing out of the pipeline is actually a mixture of injected water flows at several previous time periods, which may be seen as a combination of several water element components. The dynamic carbon emission flow model of the district heating system is established, which includes the water element carbon flow rates at the plurality of time periods, the actual outlet carbon flow rates of the pipeline at the plurality of time periods, the network loss carbon flow rates at the plurality of time periods, and the nodal carbon flow densities at the plurality of time periods. The operation of establishing the dynamic carbon emission flow model of the district heating system includes the following actions at 3-1), 3-2), 3-3) and 3-4).

[0143] 3-1) The water element is defined as a mass flow of water in a pipeline per unit time. According to the conservation of carbon emission, a carbon flow rate corresponding to mixed hot water may be similarly expressed as a linear combination of carbon flow rates corresponding to several water elements. An expression of the carbon flow rate corresponding to mixed hot water is:

[00019]$\begin{array}{cc}{\tilde{R}}_{k,t}^{\mathrm{BHS},\mathrm{out}}=\frac{R_{k,t-{\delta }_{k,t}}^{\mathrm{BHS},in}\left(B_{k,t}-\sigma M_k\right)}{\left(m_{k,t-{\delta }_{k,t}}^S\Delta t\right)}+\underset{\tau =t-{\phi }_{k,t}+1}{\overset{t-{\delta }_{k,t}-1}{.Math.}}{R}_{k,\tau }^{\mathrm{BHS},in}+\frac{R_{k,t-{\delta }_{k,t}}^{\mathrm{BHS},in}\left(m_{k,t}^S\Delta t+\sigma M_k-A_{k,t}\right)}{\left(m_{k,t-{\phi }_{k,t}}^S\Delta t\right)},& \left(27\right)\end{array}$

where {tilde over (R)}.sub.k,t.sup.BHS,out represents a water element carbon flow rate at a time period t; σ represents a density of water; M.sub.k represents a volume of the pipeline k; φ.sub.k,t represents an injection time period of an earliest water element component contained in a water flow flowing out of the pipeline k at the time period t; δ.sub.k,t represents an injection time period of a latest water element component contained in the water flow flowing out of the pipeline k at the time period t, and is related to the conveying delay of the district heating system; A.sub.k,t represents a total injection water flow amount for φ.sub.k,t time periods before the time period t; B.sub.k,t represents a total injection water flow amount from a time period t−δ.sub.k,t to the time period t; and expressions for A.sub.k,t and B.sub.k,t are:

[00020]$\begin{array}{cc}{A}_{k,t}=\{\begin{array}{cc}\underset{\tau =t-{\phi }_{k,t}+1}{\overset{t}{.Math.}}{m}_{k,\tau }^S\Delta t,& {\phi }_{k,t}\ge {\delta }_{k,t}+1\\ B_{k,t},& {\phi }_{k,t}<{\delta }_{k,t}+1\end{array},& \left(28\right)\end{array}$$\begin{array}{cc}{B}_{k,t}=\underset{\tau =t-{\delta }_{k,t}}{\overset{t}{.Math.}}{m}_{k,\tau }^S\Delta t.& \left(29\right)\end{array}$

and

[0144] 3-2) The actual outlet carbon flow rates of the pipeline at the plurality of time periods and the network loss carbon flow rates at the plurality of time periods are determined. The operation of determining the actual outlet carbon flow rates of the pipeline at the plurality of time periods and the network loss carbon flow rates at the plurality of time periods includes the following actions at 3-2-1), 3-2-2) and 3-2-3).

[0145] 3-2-1) An actual temperature of an outlet water flow of the pipeline may be obtained by taking a conveying loss of the pipeline into consideration. An expression of the actual temperature of the outlet water flow of the pipeline is as follows:

[00021]$\begin{array}{cc}{T}_{k,t}^{S,\mathrm{out}}={\tilde{T}}_{k,t}^{S,\mathrm{out}}\exp \left(-\frac{\lambda L_k}{c{m}_{k,t}^S}\right),& \left(30\right)\end{array}$

where T.sub.k,t.sup.S,out represents the actual temperature of the outlet water flow of the pipeline; {tilde over (T)}.sub.k,t.sup.S,out represents a weighted average of temperatures of injection water flows at previous time periods; λ represents a thermal conductivity coefficient of the pipeline; and L.sub.k represents a length of the pipeline k.

[0146] 3-2-2) The carbon flow density remains constant in a process in which only a temperature of a water flow decreases along a pipeline, without any mixing or an outflow of the water flow. An expression of an actual outlet carbon flow rate of the pipeline at the time period t is as follows:

[00022]$\begin{array}{cc}{R}_{k,t}^{\mathrm{BHS},\mathrm{out}}={\tilde{R}}_{k,t}^{\mathrm{BHS},\mathrm{out}}\frac{T_{k,t}^{S,\mathrm{out}}}{{\tilde{T}}_{k,t}^{S,\mathrm{out}}}.& \left(31\right)\end{array}$

[0147] 3-2-3) A network loss carbon flow rate of the pipeline at the time period t is determined. An expression of the network loss carbon flow rate of the pipeline at the time period t is as follows:

[00023]$\begin{array}{cc}{R}_{k,t}^{\mathrm{BHS},\mathrm{loss}}={\tilde{R}}_{k,t}^{\mathrm{BHS},\mathrm{loss}}\left(1-\frac{T_{k,t}^{S,\mathrm{out}}}{{\tilde{T}}_{k,t}^{S,\mathrm{out}}}\right).& \left(32\right)\end{array}$

[0148] 3-3) The nodal carbon flow densities at the plurality of time periods are determined based on the action of 1-3-2) and the action of 1-3-3). An expression of the nodal carbon flow densities at the plurality of time periods is:

[00024]$\begin{array}{cc}{\rho }_{n,t}^{NHS}=\frac{\underset{k\in {\Omega }_n^{\mathrm{BH}+}}{.Math.}\left(R_{k,t}^{\mathrm{BHS},\mathrm{out}}+X{R}_{k,t}^{\mathrm{BHS},\mathrm{Loss}}\right)}{\underset{k\in {\Omega }_n^{BH+}}{.Math.}c{m}_{k,t}^S{T}_{k,t}^{S,\mathrm{out}}}.& \left(33\right)\end{array}$

[0149] 3-4) The nodal carbon flow densities at the plurality of time periods are determined for the water return network. An expression of the nodal carbon flow densities at the plurality of time periods is:

[00025]$\begin{array}{cc}{\rho }_{n,t}^{\mathrm{NHR}}=\frac{\underset{k\in {\Omega }_n^{\mathrm{BH}-}}{.Math.}\left(R_{k,t}^{\mathrm{BHR},\mathrm{out}}+X{R}_{k,t}^{\mathrm{BHR},\mathrm{Loss}}\right)}{\underset{k\in {\Omega }_n^{\mathrm{BH}-}}{.Math.}c{m}_{k,t}^R{T}_{k,t}^{R,\mathrm{out}}}.& \left(34\right)\end{array}$

[0150] At block S103, a carbon emission amount of the current district heating system is counted based on the steady carbon emission amount and the dynamic carbon emission amount.

[0151] As illustrated in FIG. 2, a framework structure of a method for measuring carbon emission of a district heating system according to an embodiment of the present disclosure is displayed. Based on the steady carbon emission flow model and the dynamic carbon emission flow model that are constructed by the method illustrated in FIG. 2, steady and dynamic carbon emission indicators of the current district heating system are obtained.

[0152] In an embodiment of the present disclosure, a network structure and a calculation result of a steady carbon emission flow of a 6-node district heating system are as illustrated in FIG. 3. A heat load size of node 5 and node 6 are 120 MW and 100 MW, respectively. According to the calculation result of the steady carbon emission flow, carbon flow densities of two heat sources HS1 and HS2 are 0.222 tCO.sub.2/MWh and 0.267 tCO.sub.2/MWh, respectively. A total carbon flow rate of an injection of the heat source, a total carbon flow rate of the heat load, and a total carbon flow rate due to a network loss are 68.63 tCO.sub.2/h, 53.14 tCO.sub.2/h, and 15.49 tCO.sub.2/h, respectively.

[0153] According to the dynamic carbon emission flow model of the district heating system provided in the present disclosure, carbon flow densities of each node in the water supply network of the district heating system at different time periods are calculated and illustrated in FIG. 4.

[0154] With the method for measuring the carbon emission of the district heating system according to the embodiments of the present disclosure, the steady carbon emission flow model of the district heating system is established based on the carbon emission flow theory, and a corresponding matrix calculation expression is provided. The dynamic carbon emission flow model of the district heating system is established with a conveying delay in a heat flow taken into consideration. The steady and dynamic carbon emission indicators of the district heating system are obtained. In this way, carbon emission details of each link of a source, a grid, and a load of the district heating system can be effectively identified, carbon emission responsibilities on both a source side and a load side can be clarified, and accurate measurement of the carbon emission of the district heating system can be realized. In addition, with the help of a data collection, processing, and display system, a calculation result of a carbon meter system can be displayed intuitively to heat users, which can guide the heat users to reduce the carbon emission interactively, and have high application value.

[0155] An apparatus for measuring carbon emission of a district heating system according to an embodiment of the present disclosure is described below with reference to the accompanying drawings.

[0156] FIG. 5 is an example diagram of an apparatus for measuring carbon emission of a district heating system according to an embodiment of the present disclosure.

[0157] As illustrated in FIG. 5, the apparatus for measuring the carbon emission of the district heating system includes a plurality of heat source-side carbon meters, a plurality of heat pipeline carbon meters, a plurality of user-side carbon meters, and a central server.

[0158] Each of the plurality of heat source-side carbon meters is connected to a heat source and configured to measure heat source-side carbon emission.

[0159] Each of the plurality of heat pipeline carbon meters is connected to a heat pipeline and the heat source-side carbon meter. The heat pipeline carbon emission meter is configured to measure heat pipeline carbon emission.

[0160] Each of the plurality of user-side carbon meters has one end connected to a heat user end and another end connected to the heat pipeline carbon meter, and is configured to obtain a carbon potential of a node where a user is located and obtain a carbon emission amount resulted from heat consumption of the user.

[0161] As a calculation center, the central server is responsible for calculation of the carbon emission flow of the district heating system, and for communicating with the carbon meter disposed in each link to obtain carbon emission of each link. The central server is connected to the plurality of heat source-side carbon meters and the plurality of heat pipeline carbon meters, and configured to calculate a carbon emission intensity of the heat source based on coal consumption data of the heat source and an emission coefficient of the coal, calculate a distribution of a carbon emission flow in the heat pipeline, and calculate a carbon emission amount corresponding to a pipeline connected to the node based on data of the heat pipeline carbon meter.

[0162] Optionally, in an embodiment of the present disclosure, each of the heat source-side carbon meter, the heat pipeline carbon meter, and the user-side carbon meter is further configured to display in real time a result of the carbon emission of the district heating system calculated by the central server.

[0163] In the embodiments of the present disclosure, each of the heat source-side carbon meter, the heat pipeline carbon meter, and the user-side carbon meter is configured to collect and calculate information needed for measuring the carbon emission, and transmit relevant information to the central server. The central server is configured to calculate the carbon emission of the district heating system using the collected information, and feed a calculation result back to each carbon meter. The carbon meter is configured to display the calculation result fed back by the central server in real time.

[0164] In the embodiments of the present disclosure, the plurality of carbon meters forms a carbon meter system. The carbon meter system is a hierarchical system divided into an upper level and a lower level. The upper level of the carbon meter system includes the central server, the heat source-side carbon meter, and the heat pipeline carbon meter, and is configured 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, and then calculate the distribution of the carbon emission flow in the heat pipeline; and calculate the carbon emission amount corresponding to the pipeline connected to the node based on the data of the heat pipeline carbon meter. The lower level of the carbon meter system consists of user-side carbon meters. The user-side carbon meter is configured to communicate with the heat pipeline carbon meter to obtain the carbon potential of the node where the user is located and obtain the carbon emission amount resulted from the heat consumption of the user.

[0165] It should be noted that the above explanatory description of the embodiments of the method for measuring the carbon emission of the district heating system is also applicable to the apparatus for measuring the carbon emission of the district heating system, and details thereof will be omitted here.

[0166] The apparatus for measuring the carbon emission of the district heating system according to the embodiments of the present disclosure is formed by carbon meters scattered in various links of the district heating system for measuring the carbon emission, the central server, and a communication line. In some embodiments, the carbon meters scattered in various links of the district heating system include the heat source-side carbon meter, a pipeline-side carbon meter, and a heat-user side carbon meter. In this way, the carbon emission details of each link of the source, the grid, and the load of the district heating system can be effectively identified, the carbon emission responsibilities on both the source side and the load side can be clarified, and the accurate measurement of the carbon emission of the district heating system can be realized, which has high application value.

[0167] FIG. 6 is a schematic diagram showing a structure of an electronic device according to an embodiment of the present disclosure. The electronic device may include a memory 601, a processor 602, and a computer program stored in the memory 601 and executable on the processor 602. The processor 602, when executing the computer program, implements the method for measuring the carbon emission of the district heating system according to any of the above embodiments.

[0168] Further, the electronic device includes: a communication interface 603 configured for communication between the memory 601 and the processor 602; and the memory 601 stores a computer program executable on the processor 602. The memory 601 may include a high-speed Random Access Memory (RAM), or may also include a non-volatile memory, such as at least one disk memory.

[0169] When the memory 601, the processor 602, and the communication interface 603 are implemented independently, the communication interface 603, the memory 601, and the processor 602 may be interconnected and communicate with each other via a bus. The bus may be an Industry Standard Architecture (ISA) bus, a Peripheral Component (PCI) bus, or an Extended Industry Standard Architecture (EISA) bus. The buses may be divided into an address bus, a data bus, a control bus, etc. For the convenience of description, only one thick line is used in FIG. 6, but it does not mean that there is only one bus or one type of bus.

[0170] Optionally, in a specific implementation, when the memory 601, the processor 602, and the communication interface 603 are integrated on a single chip for an implementation, the memory 601, the processor 602, and the communication interface 603 may communicate with each other through an internal interface.

[0171] The processor 602 may be a Central Processing Unit (CPU), or an Application Specific Integrated Circuit (ASIC), or one or more integrated circuits configured to implement an embodiment of the present disclosure.

[0172] A computer-readable storage medium is further provided according to an embodiment. The computer-readable storage medium stores a computer program. The computer program, when executed by a processor, implements the method for measuring the carbon emission of the district heating system as described above.

[0173] In the description of this specification, descriptions with reference to the terms “an embodiment”, “some embodiments”, “examples”, “specific examples”, or “some examples” etc., mean that specific features, structure, materials or characteristics described in conjunction with the embodiment or example are included in at least one embodiment or example of the present disclosure. In this specification, the schematic representations of the above terms do not necessarily refer to the same embodiment or example. Moreover, the described specific features, structures, materials or characteristics may be combined in any one or more embodiments or examples in a suitable manner. In addition, those skilled in the art can combine the different embodiments or examples and the features of the different embodiments or examples described in this specification without contradicting each other.

[0174] In addition, the terms “first” and “second” are only used for descriptive purposes, and cannot be understood as indicating or implying relative importance or implicitly indicating the number of indicated technical features. Therefore, the features defined with “first” and “second” may explicitly or implicitly include at least one of the features. In the description of the present disclosure, “N” means at least two, such as two, three, etc., unless otherwise specifically defined.

[0175] Any process or method described in a flowchart or described herein in other ways may be understood to include one or N modules, segments, or portions of codes of executable instructions for achieving specific logical functions or steps in the process. The scope of a preferred embodiment of the present disclosure includes other implementations, in which a function may be performed not in a sequence shown or discussed, including a substantially simultaneous manner or a reverse sequence based on the function involved, which should be understood by those skilled in the art to which the embodiments of the present disclosure belong.

[0176] It should be understood that each part of the present disclosure may be realized by hardware, software, firmware, or a combination thereof. In the above embodiments, N steps or methods may be realized by software or firmware stored in the memory and executed by an appropriate instruction execution system. For example, when it is realized by the hardware, likewise in another embodiment, the steps or methods may be realized by one or a combination of the following techniques known in the art: a discrete logic circuit having a logic gate circuit for realizing a logic function of a data signal, an application-specific integrated circuit having an appropriate combination logic gate circuit, a Programmable Gate Array (PGA), a Field Programmable Gate Array (FPGA), etc.

[0177] It should be understood by those skilled in the art that all or a part of the steps carried by the method in the above-described embodiments may be completed by relevant hardware instructed by a program. The program may be stored in a computer-readable storage medium. When the program is executed, one or a combination of the steps of the method in the above-described embodiments may be included.