CN115329518A - Energy flow-based electrothermal coupling network carbon emission flow metering method, electronic equipment and storage medium - Google Patents
Energy flow-based electrothermal coupling network carbon emission flow metering method, electronic equipment and storage medium Download PDFInfo
- Publication number
- CN115329518A CN115329518A CN202211026388.3A CN202211026388A CN115329518A CN 115329518 A CN115329518 A CN 115329518A CN 202211026388 A CN202211026388 A CN 202211026388A CN 115329518 A CN115329518 A CN 115329518A
- Authority
- CN
- China
- Prior art keywords
- network
- node
- energy
- flow
- carbon
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Pending
Links
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/10—Geometric CAD
- G06F30/18—Network design, e.g. design based on topological or interconnect aspects of utility systems, piping, heating ventilation air conditioning [HVAC] or cabling
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/16—Matrix or vector computation, e.g. matrix-matrix or matrix-vector multiplication, matrix factorization
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q10/00—Administration; Management
- G06Q10/06—Resources, workflows, human or project management; Enterprise or organisation planning; Enterprise or organisation modelling
- G06Q10/063—Operations research, analysis or management
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q50/00—Systems or methods specially adapted for specific business sectors, e.g. utilities or tourism
- G06Q50/06—Electricity, gas or water supply
Abstract
An electrothermal coupling network carbon emission flow metering method based on energy flow, electronic equipment and a storage medium belong to the field of carbon emission flow metering methods. The accuracy of the carbon emission measurement of the whole process is improved. The method comprises the steps of constructing a power network carbon emission flow model; constructing a thermal network carbon emission flow model; constructing a carbon emission flow model of the energy conversion equipment; calculating the energy flow of the comprehensive energy network according to the power network carbon emission flow model and the thermal network carbon emission flow model; collecting parameters of an electrothermal coupling network, calculating the carbon emission flow of the electrothermal coupling network according to a power network carbon emission flow model, a thermal network carbon emission flow model, an energy conversion equipment carbon emission flow model and the calculation result of energy flow, generating a node load matrix, an energy matrix of an energy supply node injection network, an energy flow distribution matrix and a node energy flux matrix in the electrothermal coupling network, calculating the node carbon potential of the electrothermal coupling network, and calculating the carbon flow rate of load nodes. The invention has accurate measurement.
Description
Technical Field
The invention belongs to the field of carbon emission flow metering methods, and particularly relates to an energy flow-based electrothermal coupling network carbon emission flow metering method, electronic equipment and a storage medium.
Background
With the increasing greenhouse effect and the decreasing fossil energy, people pay more attention to the emission of greenhouse gases and the use of fossil energy. Meanwhile, under the requirement of sustainable development of environment and energy, a comprehensive energy network is also required to be developed rapidly so as to improve the use efficiency of energy and promote the application of an efficient energy-saving energy utilization mode.
Reasonable and accurate carbon measurement is an important technical support for measuring carbon emission and is also a data support. Electric power and heat supply are important industries of carbon dioxide emission, play an important role in realizing energy conservation and emission reduction, and therefore, the research on carbon emission measurement of an electric heating network is very important.
In the current stage, the carbon emission measurement of the electric heating network mainly adopts a macroscopic statistical method, namely, the total carbon emission of the electric power system and the heat energy system is calculated according to the carbon emission factors of different energy sources by counting the sum of the use of various energy sources. Problems with such methods include: on the one hand, both the power network and the heat network are network systems which consume primary energy and provide secondary energy, and are characterized in that the generation of carbon emissions is entirely from the energy supply side. The carbon emission responsibility is attributed to the energy supply side during carbon metering, the carbon emission caused by the supply and demand relationship is neglected, the carbon emission responsibility of the user side cannot be defined, and meanwhile, the energy and the energy consumption of regions are not distributed uniformly, so that the difficulty in defining the carbon emission responsibility of the user side is increased. On the other hand, the access of new energy and a macroscopic carbon emission metering method cannot meet the real-time carbon emission metering problem of an electric heating network comprising a plurality of distributed energy supply systems.
Disclosure of Invention
The invention aims to solve the problem that an electrothermal coupling network ignores carbon emission measurement on a user side, and provides an electrothermal coupling network carbon emission flow measurement method based on energy flow, electronic equipment and a storage medium in order to improve the accuracy of the carbon emission measurement of a whole flow.
In order to achieve the purpose, the invention is realized by the following technical scheme:
an electro-thermally coupled network carbon emission flow metering method based on energy flow, comprising the steps of:
an electro-thermally coupled network carbon emission flow metering method based on energy flow, comprising the steps of:
s1, constructing a power network carbon emission flow model;
s2, constructing a thermal network carbon emission flow model;
s3, constructing a carbon emission flow model of the energy conversion equipment;
s4, calculating the energy flow of the comprehensive energy network according to the power network carbon emission flow model in the step S1 and the heat power network carbon emission flow model in the step S2;
s5, collecting parameters of an electrothermal coupling network, calculating the carbon emission flow of the electrothermal coupling network according to the power network carbon emission flow model in the step S1, the thermal network carbon emission flow model in the step S2, the energy conversion equipment carbon emission flow model in the step S3 and the calculation result of the energy flow in the step S4, generating a node load matrix, an energy matrix of an energy supply node injection network, an energy flow distribution matrix and a node energy flux matrix in the electrothermal coupling network, and calculating a node carbon potential of the electrothermal coupling network;
and S6, calculating the carbon flow rate of the load node according to the node carbon potential of the electrothermal coupling network obtained in the step S5.
Further, the specific implementation method of step S1 includes the following steps:
s1.1, setting active power P and reactive power P' in the power network, as shown in a formula (1):
P i ＝V i ∑V j (G ij cosθ ij +B ij sinθ ij )
wherein, P i Active power of node i, P' i Is the reactive power of a node i, i and j respectively represent the node serial number, ij is a branch circuit connected with the node i and the node j, V is the node voltage, G ij As the real part of the nodal admittance matrix, B ij Is the imaginary part of the node admittance matrix, and theta is the included angle of the voltage and the current at the branch;
s1.2, setting carbon flow F of power network node i ei As shown in equation (2):
wherein λ is eik Carbon emission factor, P, of kth generator connected to node i eik The active power of the kth generator accessed to the node i;
s1.3, setting branch carbon flow rate R of power network branch ij eij As shown in equation (3):
wherein the carbon flow of the power network node j is F ej T is time;
s1.4, setting branch carbon flow density rho of power network branch ij eij As shown in equation (4):
wherein, P eij Branch active power for power network branch ij;
s1.5, setting the carbon potential of a power network node i, as shown in a formula (5):
wherein, P j Is the active power, rho, of a node j in the power network ej Is the carbon flow density of power network node j;
s1.6, setting active loss of a power network branch to be P l And then the carbon current rate R of active loss of the branch es As shown in equation (6):
R es ＝diag(E ei )·P l (6)
wherein diag is the diagonal matrix.
Further, the specific implementation method of step S2 includes the following steps:
s2.1, setting hydraulic balance of nodes of the thermodynamic network, as shown in a formula (7):
wherein the content of the first and second substances,
in order to be able to take the mass flow into account,
in order to be able to flow out of the mass flow,
is the mass flow consumed;
s2.2, setting the active power of the thermodynamic network, as shown in a formula (8):
where Φ is the active power of the node, T s For supply of heat, T o Is the outlet temperature, C p Is the specific heat capacity of water;
s2.3, setting the temperature of the tail end of the pipeline of the thermal network, as shown in a formula (9):
wherein, T end Is the pipe end temperature, T start Is the initial temperature of the pipeline, T a Lambda is the heat transfer coefficient of the unit pipe, and L is the length of the pipe;
s2.4, setting the carbon flow F of the node i of the thermodynamic network hi As shown in equation (10):
wherein λ is hik Carbon emission factor, Φ, of kth heat-producing plant connected to node i hik Active power of the kth heat production equipment accessed for the node i;
s2.5 setting the Branch carbon flow Rate R of the Branch ij of the thermodynamic network hij As shown in equation (11):
s2.6, setting branch carbon flow density rho of branch ij of thermal power network hij As shown in equation (12):
wherein phi hij Branch active power for the thermodynamic network branch ij;
s2.7, setting the carbon potential of the thermodynamic network node i as shown in a formula (13):
wherein phi j Active power of node j, ρ hj Is the carbon flow density of the thermodynamic network node j;
s2.8, setting the energy loss of the branch of the thermodynamic network to phi l Then branch active loss carbon current rate R hs As shown in equation (14):
R hs ＝diag(E hi )·Φ l (14)。
further, the step S3 is that the energy conversion device is a CHP electric heating network coupling system, and the specific implementation method for establishing the CHP carbon emission flow model includes the following steps:
s3.1, setting the power supply quantity A of the CHP electric heating network coupling system E The heat supply amount B is shown in the formula (15) H As shown in the formula (16)
A E ＝Q·η e (15)
B H ＝Q·η h ·λ a (16)
Wherein Q is consumption of natural gas of the CHP electric heating network coupling system, eta e Efficiency of electric power generation of prime mover [ ] h For waste heat recovery efficiency, lambda a Is the performance coefficient of the absorption heat pump;
s3.2, setting a carbon potential E of a power network node of a CHP electric heating network coupling system e The thermodynamic network node carbon potential E is shown in equation (17) h As shown in equation (18):
wherein the content of the first and second substances,
is the carbon emission factor of natural gas.
Further, step S4 is a method for calculating an energy flow of the integrated energy network, where each node has a known quantity and an unknown quantity according to the parameters of the integrated energy network, the known quantity is input into the power network carbon emission flow model and the thermal network carbon emission flow model constructed in steps S1-S2, the unknown quantity of each node is calculated by using a newton-raphson algorithm, and an obtained calculation result is the integrated energy network energy flow, and the specific implementation method includes the following steps:
s4.1, the known quantity of a balance node of the power network of the comprehensive energy network is V and theta, the unknown quantity of the balance node is P and P ', the calculation formula is formula (1), V represents the node voltage, theta represents the included angle of the voltage and the current at the branch, P represents the active power of the node, and P' represents the reactive power of the node;
s4.2, the known quantities of PV nodes of the power network of the comprehensive energy network are V and P, the unknown quantities are theta and P', and the calculation formula is a formula (1);
s4.3, the known quantity of the PQ node of the power network of the comprehensive energy network is P, P', the unknown quantity is V, theta, and the calculation formula is formula (1);
s4.4, the known quantity of balancing nodes of the thermodynamic network of the integrated energy network is T s Unknown quantity is phi, T r ,
The calculation formula is formula (7, 8, 9), T r The temperature of the return water is the temperature of the return water,
is the mass flow;
s4.5, phi T of thermodynamic network of integrated energy network s The known quantity of a node is phi, T s Unknown quantity is T r ,
The calculation formula is formula (7, 8, 9);
s4.6 phi T of thermodynamic network of comprehensive energy network r The known quantity of a node is phi, T r Unknown quantity is T s ,
The calculation formula is formula (7, 8, 9).
Further, the specific implementation method of step S5 includes the following steps:
s5.1, setting nodes of the electric-thermal coupling network into two types, namely an energy supply node and an energy utilization node, wherein the position and the carbon emission intensity of the energy supply node are known, and calculating the carbon emission intensity in the power network by using a formula (2, 3, 4, 5 and 6) according to the calculation result of the energy flow; calculating the carbon emission intensity of the thermodynamic network by using formulas (10, 11, 12, 13 and 14); the carbon emission intensity of the CHP node is obtained through calculation of the formulas (15, 16, 17 and 18), and the carbon emission intensity is the node carbon potential;
s5.2, setting a power network in the electrothermal coupling network to be provided with K generators, and setting K heating equipment in the heating network, wherein the carbon potential matrix E of the nodes of the power network of the electrothermal coupling network G A carbon potential matrix E of the nodes of the thermodynamic network, as shown in equation (19) B As shown in equation (20):
E G ＝[e e1 e e2 … e eK ] T (19)
E B ＝[e h1 e h2 … e hK ] T (20)
s5.3, setting the load node matrixes of the electrothermal coupling networks to be P respectively L ,Φ L The energy matrixes of the energy supply node injection networks are respectively P G ,Φ B The energy flow distribution matrix is respectively P e ,Φ h The node energy flux matrix is P N ,Φ N Then, obtaining the node carbon potential of the electrothermal coupling network, which is expressed as the formulas (21) and (22):
further, the formula for calculating the carbon flow rate of the load node in step S6 is:
R eL ＝P L E e (23)
R hL ＝Φ L E h (24)
wherein R is eL Carbon flow rate, R, of load node of electric power network for electro-thermal coupling network hL The carbon flow rate of a thermal network load node of the electro-thermal coupling network.
Electronic apparatus comprising a memory storing a computer program and a processor implementing the steps of the energy flow-based method of electrothermal coupled network carbon emission flow metering, when the computer program is executed by the processor.
A computer-readable storage medium having stored thereon a computer program which, when executed by a processor, implements a method of energy flow-based electrothermal coupled network carbon exhaust flow metering.
The invention has the beneficial effects that:
the invention relates to an energy flow-based carbon emission flow measurement method of an electrothermal coupling network, wherein the electrothermal coupling network is a network which takes CHP as a representative and can simultaneously generate electric power and heat power, or takes electric boilers and other devices which can convert electric energy and heat energy into each other and connects a power supply network and a heat supply network with each other. Due to the coupling of different energy networks, energy conversion is generated, and the time-space scales and the time scales of electric power and heating power are different, so that the decoupling analysis of different networks and coupling equipment is needed for the analysis of the electric-thermal coupling network. The invention provides a full-process carbon emission metering model for an electrothermal coupling network. The model is a carbon emission allocation model based on carbon emission flow and network energy conversion and loss, takes the energy loss of each link into consideration, and effectively allocates the responsibility of carbon emission.
The energy flow-based carbon emission flow metering method of the electric heating coupling network can reasonably define carbon emission responsibility and effectively promote the development of clean energy.
Drawings
FIG. 1 is a flow chart of a method for measuring carbon emission flow in an energy flow-based electro-thermally coupled network according to the present invention;
FIG. 2 is a schematic view of the overall flow of the electro-thermal coupling network according to the present invention;
FIG. 3 is a schematic diagram of the CHP energization relationship of the electrothermal coupling network according to the present invention;
FIG. 4 is a schematic diagram of an exemplary power network according to a first embodiment of the present invention;
fig. 5 is a schematic diagram of an exemplary thermal network structure according to a first embodiment of the present invention.
Detailed Description
In order to make the objects, technical solutions and advantages of the present invention more apparent, the present invention will be described in detail with reference to the accompanying drawings and the detailed description. It is to be understood that the embodiments described herein are illustrative only and are not limiting, i.e., that the embodiments described are only a few embodiments, rather than all, of the present invention. While the components of the embodiments of the present invention generally described and illustrated in the figures herein may be arranged and designed in a wide variety of different configurations, the present invention is capable of other embodiments.
Thus, the following detailed description of specific embodiments of the present invention presented in the accompanying drawings is not intended to limit the scope of the invention as claimed, but is merely representative of selected embodiments of the invention. All other embodiments, which can be derived by a person skilled in the art from the detailed description of the invention without inventive step, are within the scope of protection of the invention.
For a further understanding of the contents, features and effects of the present invention, the following embodiments will be illustrated in detail with reference to the accompanying drawings 1-5:
the first specific implementation way is as follows:
the whole energy transfer process of the electric-thermal coupling network can be divided into four parts, namely an energy production link, an energy transmission link, an energy conversion link and an energy use link. In the energy production link, primary energy is consumed to generate heat energy and electric energy along with carbon emission; the energy is transmitted through the power network and the heat network, and the loss of the energy in the process is equal to the consumption of carbon emission; CHP is mainly considered in the energy conversion link; although the carbon emission is not directly generated in the energy use link, the carbon emission corresponding to the energy also needs to be calculated in consideration of the responsibility of the carbon emission corresponding to the energy demand. The schematic flow diagram of the electrothermal coupling network is shown in fig. 2. Relationship between the relative index of carbon emission flow and energy flow in the electrothermal coupling network: relevant indexes of carbon emission flow of the electrothermal coupling network are called carbon flow indexes and comprise branch carbon flow, branch carbon flow rate, branch carbon flow density and node carbon potential; the branch carbon flow corresponds to the electric quantity transmitted by the branch in the power flow, and corresponds to the heat transmitted by the branch in the thermal power flow. The branch carbon flow rate corresponds to the active power of the branch in the power flow and corresponds to the thermal power transmitted by the branch in the thermal flow.
An energy flow based electro-thermal coupling network carbon emission flow metering method, comprising the steps of:
s1, constructing a power network carbon emission flow model;
further, the specific implementation method of step S1 includes the following steps:
s1.1, setting active power P and reactive power P' in the power network, as shown in formula (1):
P i ＝V i ∑V j (G ij cosθ ij +B ij sinθ ij )
wherein, P i Active power of node i, P' i I and j represent node serial numbers respectively, ij is the connection of the node i and the node jConnected branch, V being the node voltage, G ij Being the real part of the nodal admittance matrix, B ij Is the imaginary part of the node admittance matrix, and theta is the included angle of the voltage and the current at the branch;
s1.2, setting carbon flow F of power network node i ei As shown in equation (2):
wherein λ is eik Carbon emission factor, P, of kth generator connected to node i eik The active power of a kth generator accessed to the node i;
s1.3, setting branch carbon flow rate R of power network branch ij eij As shown in equation (3):
wherein the carbon flow of the power network node j is F ej T is time;
s1.4, setting branch carbon flow density rho of power network branch ij eij As shown in equation (4):
wherein, P eij Branch active power for power network branch ij;
s1.5, setting the carbon potential of a power network node i, as shown in a formula (5):
wherein, P j Active power, rho, for node j in the power network ej Is the carbon flow density of power network node j;
s1.6, setting the active loss of a power network branch asP l Then branch active loss carbon current rate R es As shown in equation (6):
R es ＝diag(E ei )·P l (6)
wherein diag is to solve a diagonal matrix;
s2, constructing a thermal network carbon emission flow model;
furthermore, the heat supply network transfers heat through the hot water in the pipeline, and the structure of the heat supply network consists of a water supply pipeline and a water return pipeline. A thermodynamic network differs from an electrical network in that the flow of water is of a type such that, in modeling the thermodynamic network, in addition to the thermodynamic model, the hydraulic balance of the fluid is also taken into account. The flow of liquid in the pipe should obey kirchhoff's law. The outlet temperature of the heat source and the inlet temperature of the load are known in the thermodynamic model, and the return water temperature depends on the heating temperature, the room temperature, and the heat load. Heat loss in the pipe is affected by the length of the pipe and the ambient temperature.
Further, the specific implementation method of step S2 includes the following steps:
s2.1, setting hydraulic balance of nodes of the thermodynamic network, as shown in a formula (7):
wherein, the first and the second end of the pipe are connected with each other,
in order to be able to take the mass flow into account,
in order to be able to flow out of the mass flow,
is the mass flow consumed;
s2.2, setting the active power of the thermodynamic network, as shown in a formula (8):
where Φ is the active power of the node, T s For heating temperature, T o Is the outlet temperature, C p Is the specific heat capacity of water;
s2.3, setting the temperature of the tail end of the pipeline of the heat distribution network, as shown in a formula (9):
wherein, T end Is the temperature at the end of the pipe, T start Is the initial temperature of the pipeline, T a Lambda is the heat transfer coefficient of the unit pipe, L is the length of the pipe;
s2.4, setting the carbon flow F of the node i of the thermodynamic network hi As shown in equation (10):
wherein λ is hik Carbon emission factor, Φ, of kth heat-producing plant connected to node i hik Active power of the kth heat production equipment accessed for the node i;
s2.5 setting the Branch carbon flow Rate R of the Branch ij of the thermodynamic network hij As shown in equation (11):
s2.6, setting branch carbon flow density rho of branch ij of thermal network hij As shown in equation (12):
wherein phi is hij Branch active power for the thermodynamic network branch ij;
s2.7, setting the carbon potential of the thermodynamic network node i as shown in a formula (13):
wherein phi j Active power of node j, ρ hj Is the carbon flow density of the thermodynamic network node j;
s2.8, setting the energy loss of the branch of the thermodynamic network to phi l And then the carbon current rate R of active loss of the branch hs As shown in equation (14):
R hs ＝diag(E hi )·Φ l (14)；
s3, constructing a carbon emission flow model of the energy conversion equipment;
further, the CHP is a coupling device in the electric heating network, which can generate electric energy and heat simultaneously, and is an important device for energy conversion. Cogeneration units can increase the efficiency of heat generation, but their overall efficiency is constant. The energy supply relation of the CHP is shown in figure 3, the CHP consumes natural gas, and the electricity generation link directly generates electric energy through a gas generator; the heat generation is realized by utilizing the waste heat through an absorption heat pump and a refrigeration technology.
Further, the step S3 is that the energy conversion device is a CHP electric heating network coupling system, and the specific implementation method for establishing the CHP carbon emission flow model includes the following steps:
s3.1, setting the power supply quantity A of the CHP electric heating network coupling system E As shown in the formula (15), the heat supply amount B H As shown in equation (16):
A E ＝Q·η e (15)
B H ＝Q·η h ·λ a (16)
wherein Q is consumption of natural gas of CHP electric heating network coupling system, eta e Efficiency of electric power generation of prime mover [ ] h For waste heat recovery device efficiency, λ a Is the absorption heat pump coefficient of performance;
s3.2, setting a carbon potential E of a power network node of a CHP electric heating network coupling system e As shown in equation (17), heatForce network node carbon potential E h As shown in equation (18):
wherein, the first and the second end of the pipe are connected with each other,
a carbon emission factor for natural gas;
s4, calculating the energy flow of the comprehensive energy network according to the power network carbon emission flow model in the step S1 and the thermal network carbon emission flow model in the step S2;
further, step S4 is a method for calculating the energy flow of the integrated energy network, where each node has a known quantity and an unknown quantity according to the parameters of the integrated energy network, the known quantity is input into the power network carbon emission flow model and the thermal network carbon emission flow model constructed in steps S1-S2, the unknown quantity of each node is calculated by using a newton-raphson algorithm, and the obtained calculation result is the integrated energy network energy flow, and the specific implementation method includes the following steps:
s4.1, the known quantity of a balance node of the power network of the comprehensive energy network is V and theta, the unknown quantity of the balance node is P and P ', the calculation formula is formula (1), V represents the node voltage, theta represents the included angle of the voltage and the current at the branch, P represents the active power of the node, and P' represents the reactive power of the node;
s4.2, the known quantities of PV nodes of the power network of the comprehensive energy network are V and P, the unknown quantities are theta and P', and the calculation formula is a formula (1);
s4.3, the known quantity of the PQ node of the power network of the comprehensive energy network is P, P', the unknown quantity is V, theta, and the calculation formula is a formula (1);
s4.4, the known quantity of balancing nodes of the thermodynamic network of the integrated energy network is T s Unknown quantity is phi, T r ,
The calculation formula is formula (7, 8, 9), T r The temperature of the return water is;
s4.5, phi T of thermodynamic network of comprehensive energy network s The known quantity of a node is phi, T s Unknown quantity is T r ,
The calculation formula is formula (7, 8, 9);
s4.6 phi T of thermodynamic network of integrated energy network r The known quantity of a node is phi, T r Unknown quantity is T s ,
The calculation formula is formula (7, 8, 9);
further, the energy flow calculation requirements are shown in table 1:
TABLE 1 requirements for energy flow calculation
S5, acquiring parameters of an electrothermal coupling network, calculating the carbon emission flow of the electrothermal coupling network according to the electric power network carbon emission flow model in the step S1, the heating power network carbon emission flow model in the step S2, the energy conversion equipment carbon emission flow model in the step S3 and the calculation result of the energy flow in the step S4, generating a node load matrix, an energy matrix of an energy supply node injection network, an energy flow distribution matrix and a node energy flux matrix in the electrothermal coupling network, and calculating a node carbon potential of the electrothermal coupling network;
further, the specific implementation method of step S5 includes the following steps:
s5.1, setting nodes of the electric-thermal coupling network into two types, namely an energy supply node and an energy utilization node, wherein the position and the carbon emission intensity of the energy supply node are known, and calculating the carbon emission intensity in the power network by using a formula (2, 3, 4, 5 and 6) according to the calculation result of the energy flow; calculating the carbon emission intensity of the thermodynamic network by using formulas (10, 11, 12, 13 and 14); the carbon emission intensity of the CHP node is obtained through calculation of the formulas (15, 16, 17 and 18), and the carbon emission intensity is the node carbon potential;
s5.2, setting a power network in the electrothermal coupling network to be provided with K generators, and setting K heating equipment in the heating network, wherein the carbon potential matrix E of the nodes of the power network of the electrothermal coupling network G The carbon potential matrix E of the nodes of the thermodynamic network, as shown in equation (19) B As shown in equation (20):
E G ＝[e e1 e e2 ... e eK ] T (19)
E B ＝[e h1 e h2 ... e hK ] T (20)
s5.3, setting the load node matrixes of the electrothermal coupling networks to be P respectively L ,Φ L The energy matrixes of the energy supply nodes injected into the network are respectively P G ,Φ B The energy flow distribution matrix is respectively P e ,Φ h The node energy flux matrix is P N ,Φ N Then, the node carbon potential of the electrothermal coupling network can be expressed as shown in formulas (21), (22):
the load node matrixes of the electrothermal coupling network are respectively P L ,Φ L The energy matrixes of the energy supply nodes injected into the network are respectively P G ,Φ B The energy flow distribution matrix is respectively P e ,Φ h The node energy flux matrix is P N ,Φ N ；
The node load matrix and the energy matrix of the energy supply node injection network are matrixes formed by arranging known quantities in a table II, and the energy flow distribution matrix and the node flux matrix are obtained by calculation according to formulas 1 to 9.
Further, the carbon emission flow calculation requirements are shown in table 2:
TABLE 2 carbon emissions stream calculated requirements
S6, calculating the carbon flow rate of the load node according to the node carbon potential of the electrothermal coupling network obtained in the step S5;
further, the formula for calculating the carbon flow rate of the load node in step S6 is as follows:
R eL ＝P L E e (23)
R hL ＝Φ L E h (24)
wherein R is eL Carbon flow rate, R, of load node of electric power network for electro-thermal coupling network hL The carbon flow rate of a thermal network load node of the electro-thermal coupling network.
In order to verify the accuracy of the energy flow-based electrothermal coupling network carbon emission flow metering method proposed by the embodiment, a power network as shown in fig. 4 is adopted. The network is a standard IEEE-14 node power grid, wherein nodes 1 and 2 are connected with coal-fired power generating units, and the carbon emission intensity is 875 and 525]gCO 2 /kWh. The nodes 3 and 7 are connected with the wind turbine generator set, and the carbon emission intensity of the wind turbine generator set is zero. The power grid node 6 is connected with the heat supply network node 5 and has CHP electric heating carbon emission intensity of 520,480]gCO 2 /kWh. The heat network model is shown in FIG. 5, wherein the node 1 is a coal-fired central heating plant and has carbon emission intensity of 965gCO 2 /kWh. The load at each node and the length of the pipeline are indicated. A total of 14 nodes can be seen from the figure.
Tables 3 and 4 show the carbon emissions of the entire process of the electrothermal coupling network calculated by the method of the present embodiment.
TABLE 3 electrothermal coupling network node carbon potential and carbon flow rate
TABLE 4 carbon flow rate loss for branches of the electrothermal coupling network
Table 3 shows the carbon potential and the load carbon flow rate of each node of the power grid and the heat supply grid, and it can be seen from table 3 that the carbon potential of node 1 is equal to the carbon emission intensity of the node unit, because no energy is injected into node 1 except for the electric energy of the generator. The node 7 is connected with a wind generating set, and carbon emission of the wind generating set is zero, so that the node carbon potential is zero. The flow rate of carbon injected by the unit corresponds to the energy value generated in unit time, and the flow rate of loaded carbon corresponds to the active power consumed by the node. In addition, the difference exists between the sum of the carbon flow rates injected by the units and the sum of the carbon flow rates of the loads, and the difference is the carbon flow rate of the system loss. The carbon flow rate consumed by the system is given in table 6, the carbon flow rate of the branch loss is related to the energy loss of the branch, the power network is mainly represented by the impedance of the line, and the heat network is directly related to the length of the pipeline and the heat conduction coefficient of the pipeline. The data in tables 3 and 4 are consistent with the law of conservation of energy, i.e., conservation of carbon emissions.
Comparative example 1 is a carbon emission flow calculation using a conventional algorithm on the constructed electro-thermal coupling networks of fig. 4 and 5, and the results are shown in table 5: as can be seen from table 5, in contrast to the method proposed in the present embodiment, the method of table 5 does not consider the situation of the full flow apportionment, and the carbon emissions of the network loss are all calculated to the user side, which is obviously unreasonable.
TABLE 5 calculation of system node carbon potential and carbon flow rate for comparative example 1
Comparative example 2 is that the carbon emission measurement at the user side can be relatively simple under the condition of adopting centralized power supply and centralized heat supply, but at the same time, the energy conservation and emission reduction cannot be well promoted, so in the above example, the distributed photovoltaic power generation system is connected to the grid 11 node, and the carbon emission calculation result is shown in table 6:
TABLE 6 calculation of system node carbon potential and carbon flow rate in comparative example 2
As can be seen from comparison between table 6 and the method of the present embodiment, table 6 adds a distributed photovoltaic power generation system to node 11, the carbon emission intensity of the photovoltaic power generation system is zero, the generated energy can be supplied to itself, and the remaining power is on line. The carbon potential at node 11 is 0 and the carbon flow rate also becomes zero. At the same time, the carbon flow rate of the three nodes is also reduced due to the transmission of the clean power generated by the clean power to the nodes 6, 9 and 10.
Therefore, the electrothermal coupling network carbon emission flow metering method based on the energy flow can reasonably define carbon emission responsibility and effectively promote the development of clean energy.
The second embodiment is as follows:
an electronic device comprising a memory storing a computer program and a processor that when executed implements the steps of a method for energy flow based electrothermal coupled network carbon effluent flow metering described in the detailed description.
The computer device of the present invention may be a device including a processor, a memory, and the like, for example, a single chip microcomputer including a central processing unit and the like. And the processor is used for implementing the steps of the recommendation method capable of modifying the relationship-driven recommendation data based on the CREO software when executing the computer program stored in the memory.
The Processor may be a Central Processing Unit (CPU), other general purpose Processor, a Digital Signal Processor (DSP), an Application Specific Integrated Circuit (ASIC), an off-the-shelf Programmable Gate Array (FPGA) or other Programmable logic device, discrete Gate or transistor logic, discrete hardware components, etc. A general purpose processor may be a microprocessor or the processor may be any conventional processor or the like.
The memory may mainly include a storage program area and a storage data area, wherein the storage program area may store an operating system, an application program required by at least one function (such as a sound playing function, an image playing function, etc.), and the like; the storage data area may store data (such as audio data, a phonebook, etc.) created according to the use of the cellular phone, etc. In addition, the memory may include high speed random access memory, and may also include non-volatile memory, such as a hard disk, a memory, a plug-in hard disk, a Smart Media Card (SMC), a Secure Digital (SD) Card, a Flash memory Card (Flash Card), at least one magnetic disk storage device, a Flash memory device, or other volatile solid state storage device.
The third concrete implementation mode:
a computer-readable storage medium, on which a computer program is stored, which, when being executed by a processor, carries out a method of energy flow based electrothermal coupled network carbon emission flow metering according to embodiments.
The computer readable storage medium of the present invention may be any form of storage medium read by a processor of a computer device, including but not limited to non-volatile memory, ferroelectric memory, etc., on which a computer program is stored, which when read and executed by the processor of the computer device, may implement the steps of the above-described CREO software-based modeling method that can modify relationship-driven modeling data.
The computer program comprises computer program code which may be in the form of source code, object code, an executable file or some intermediate form, etc. The computer-readable medium may include: any entity or device capable of carrying the computer program code, recording medium, usb disk, removable hard disk, magnetic disk, optical disk, computer Memory, read-Only Memory (ROM), random Access Memory (RAM), electrical carrier wave signals, telecommunications signals, software distribution medium, and the like. It should be noted that the computer-readable medium may contain suitable additions or subtractions depending on the requirements of legislation and patent practice in jurisdictions, for example, in some jurisdictions, computer-readable media may not include electrical carrier signals or telecommunication signals in accordance with legislation and patent practice.
It is noted that relational terms such as "first" and "second," and the like, may be used solely to distinguish one entity or action from another entity or action without necessarily requiring or implying any actual such relationship or order between such entities or actions. Also, the terms "comprises," "comprising," or any other variation thereof, are intended to cover a non-exclusive inclusion, such that a process, method, article, or apparatus that comprises a list of elements does not include only those elements but may include other elements not expressly listed or inherent to such process, method, article, or apparatus. Without further limitation, an element defined by the phrases "comprising one of 8230; \8230;" 8230; "does not exclude the presence of additional like elements in a process, method, article, or apparatus that comprises the element.
While the application has been described above with reference to specific embodiments, various modifications may be made and equivalents may be substituted for elements thereof without departing from the scope of the application. In particular, the various features of the embodiments disclosed herein may be used in any combination that is not inconsistent with the structure, and the failure to exhaustively describe such combinations in this specification is merely for brevity and resource conservation. Therefore, it is intended that the application not be limited to the particular embodiments disclosed, but that the application will include all embodiments falling within the scope of the appended claims.
Claims (9)
1. An electrothermal coupling network carbon emission flow metering method based on energy flow is characterized in that: the method comprises the following steps:
s1, constructing a power network carbon emission flow model;
s2, constructing a thermal network carbon emission flow model;
s3, constructing a carbon emission flow model of the energy conversion equipment;
s4, calculating the energy flow of the comprehensive energy network according to the power network carbon emission flow model in the step S1 and the heat power network carbon emission flow model in the step S2;
s5, collecting parameters of an electrothermal coupling network, calculating the carbon emission flow of the electrothermal coupling network according to the power network carbon emission flow model in the step S1, the thermal network carbon emission flow model in the step S2, the energy conversion equipment carbon emission flow model in the step S3 and the calculation result of the energy flow in the step S4, generating a node load matrix, an energy matrix of an energy supply node injection network, an energy flow distribution matrix and a node energy flux matrix in the electrothermal coupling network, and calculating a node carbon potential of the electrothermal coupling network;
and S6, calculating the carbon flow rate of the load node according to the node carbon potential of the electric-thermal coupling network obtained in the step S5.
2. An energy flow based electrothermal coupling network carbon emission flow metering method of claim 1, wherein: the specific implementation method of the step S1 comprises the following steps:
s1.1, setting active power P and reactive power P' in the power network, as shown in formula (1):
P i ＝V i ∑V j (G ij cosθ ij +B ij sinθ ij )
wherein, P i Active power of node i, P' i I and j represent node serial numbers respectively, ij is a branch circuit connected with the node i and the node j, V is node voltage, and G is node voltage ij Being the real part of the nodal admittance matrix, B ij Is the imaginary part of the node admittance matrix, and theta is the included angle of the voltage and the current at the branch;
s1.2, setting carbon flow F of power network node i ei As shown in equation (2):
wherein λ is eik Carbon emission factor, P, of kth generator connected to node i eik The active power of a kth generator accessed to the node i;
s1.3, setting branch carbon flow rate R of power network branch ij eij As shown in equation (3):
wherein the carbon flow of the power network node j is F ej T is time;
s1.4, setting branch carbon current density rho of power network branch ij eij As shown in equation (4):
wherein, P eij Branch active power for power network branch ij;
s1.5, setting the carbon potential of the power network node i as shown in a formula (5):
wherein, P j Active power, rho, for node j in the power network ej Is the carbon flow density of power network node j;
s1.6, setting active loss of a power network branch to be P l Then branch active loss carbon current rate R es As shown in equation (6):
R es ＝diag(E ei )·P l (6)
wherein diag is the diagonal matrix.
3. A method of energy flow based electrothermal coupled network carbon emission flow metering according to claim 1 or 2, wherein: the specific implementation method of the step S2 comprises the following steps:
s2.1, setting hydraulic balance of nodes of the thermodynamic network, as shown in a formula (7):
wherein, the first and the second end of the pipe are connected with each other,
in order to be able to take the mass flow into account,
in order to be able to take out the mass flow,
is the mass flow consumed;
s2.2, setting the active power of the heat distribution network, as shown in a formula (8):
where Φ is the active power of the node, T s For supply of heat, T o Is the outlet temperature, C p Is the specific heat capacity of water;
s2.3, setting the temperature of the tail end of the pipeline of the thermal network, as shown in a formula (9):
wherein, T end Is the temperature at the end of the pipe, T start Is the initial temperature of the pipeline, T a Lambda is the heat transfer coefficient of the unit pipe, L is the length of the pipe;
s2.4, setting the carbon flow F of the node i of the thermodynamic network hi As shown in equation (10):
wherein λ is hik Carbon emission factor, Φ, of kth heat producing equipment connected to node i hik Active power of the kth heat generating equipment accessed for the node i;
s2.5 setting the Branch carbon flow Rate R of the Branch ij of the thermodynamic network hij As shown in equation (11):
s2.6, setting branch carbon flow density rho of branch ij of thermal network hij As shown in equation (12):
wherein phi is hij Branch active power for the thermodynamic network branch ij;
s2.7, setting the carbon potential of the thermodynamic network node i as shown in a formula (13):
wherein phi j Active power of node j, ρ hj Is the carbon flow density of the thermodynamic network node j;
s2.8, setting the energy loss of the branch of the thermodynamic network to phi l Then branch active loss carbon current rate R hs As shown in equation (14):
R hs ＝diag(E hi )·Φ l (14)。
4. the energy flow-based method for metering carbon emission flow in an electrothermal coupling network according to claim 3, wherein: s3, the energy conversion equipment is a CHP electric heating network coupling system, and the specific implementation method for establishing the CHP carbon emission flow model comprises the following steps:
s3.1, setting the power supply quantity A of the CHP electric heating network coupling system E The heat supply amount B is shown in the formula (15) H As shown in equation (16):
A E ＝Q·η e (15)
B H ＝Q·η h ·λ a (16)
wherein Q is consumption of natural gas of CHP electric heating network coupling system, eta e For the efficiency of the prime mover h For waste heat recovery device efficiency, λ a Is the performance coefficient of the absorption heat pump;
s3.2, setting a carbon potential E of a power network node of a CHP electric heating network coupling system e As shown in equation (17), the thermal network node carbon potential E h As shown in equation (18):
wherein the content of the first and second substances,
is the carbon emission factor of natural gas.
5. The energy flow-based method for metering carbon emission flow in an electrothermal coupling network according to claim 4, wherein: step S4, the method for calculating the energy flow of the comprehensive energy network comprises the steps of inputting the known quantity into the power network carbon emission flow model and the thermal network carbon emission flow model which are constructed in the steps S1-S2 according to the parameters of the comprehensive energy network, wherein each node has the known quantity and the unknown quantity, calculating the unknown quantity of each node by using a Newton-Raphson algorithm, and obtaining a calculation result as the energy flow of the comprehensive energy network, wherein the specific implementation method comprises the following steps:
s4.1, the known quantity of a balance node of the power network of the comprehensive energy network is V and theta, the unknown quantity of the balance node is P and P ', the calculation formula is formula (1), V represents the node voltage, theta represents the included angle of the voltage and the current at the branch, P represents the active power of the node, and P' represents the reactive power of the node;
s4.2, the known quantities of PV nodes of the power network of the comprehensive energy network are V and P, the unknown quantities are theta and P', and the calculation formula is a formula (1);
s4.3, the known quantity of the PQ node of the power network of the comprehensive energy network is P, P', the unknown quantity is V, theta, and the calculation formula is formula (1);
s4.4, the known quantity of the balance nodes of the thermodynamic network of the integrated energy network is T s Unknown quantity is phi, T r ,
The calculation formula is formula (7, 8, 9), T r The temperature of the return water is set as the temperature of the return water,
mass flow (please give meaning);
s4.5, phi T of thermodynamic network of comprehensive energy network s The known quantity of a node is phi, T s Unknown quantity is T r ,
The calculation formula is formula (7, 8, 9);
s4.6 phi T of thermodynamic network of comprehensive energy network r The known quantity of a node is phi, T r Unknown quantity is T s ,
The calculation formula is formula (7, 8, 9).
6. The energy flow-based method for metering carbon emission flow in an electrothermal coupling network according to claim 5, wherein: the specific implementation method of the step S5 comprises the following steps:
s5.1, setting nodes of the electric-thermal coupling network into two types, namely energy supply nodes and energy utilization nodes, wherein the positions of the energy supply nodes and the carbon emission intensity are known, and calculating the carbon emission intensity in the power network by using formulas (2, 3, 4, 5 and 6) according to the calculation result of the energy flow; calculating the carbon emission intensity of the heat power network by formulas (10, 11, 12, 13 and 14); the carbon emission intensity of the CHP node is obtained through calculation of the formulas (15, 16, 17 and 18), and the carbon emission intensity is the node carbon potential;
s5.2, setting a power network in the electrothermal coupling network to be provided with K generators, and setting K heating equipment in the heating network, wherein the carbon potential matrix E of the nodes of the power network of the electrothermal coupling network G A carbon potential matrix E of the nodes of the thermodynamic network, as shown in equation (19) B As shown in equation (20):
E G ＝[e e1 e e2 … e eK ] T (19)
E B ＝[e h1 e h2 … e hK ] T (20)
s5.3, setting the load node matrixes of the electrothermal coupling networks to be P respectively L ,Φ L The energy matrixes of the energy supply node injection networks are respectively P G ,Φ B The energy flow distribution matrix is P e ,Φ h The node energy flux matrix is P N ,Φ N Then, obtaining the node carbon potential of the electrothermal coupling network, which is expressed as the formulas (21) and (22):
7. the energy flow-based method for metering carbon emission flow in an electrothermal coupling network according to claim 6, wherein: the formula for calculating the carbon flow rate of the load node in the step S6 is as follows:
R eL ＝P L E e (23)
R hL ＝Φ L E h (24)
wherein R is eL Carbon flow rate, R, of load node of electric power network for electro-thermal coupling network hL The carbon flow rate is the thermal network load node of the electro-thermal coupling network.
8. Electronic device, characterized in that it comprises a memory storing a computer program and a processor implementing the steps of an energy flow based electro-thermal coupled network carbon emission flow metering method according to any one of claims 1-7 when executing said computer program.
9. A computer-readable storage medium, on which a computer program is stored, which, when being executed by a processor, carries out a method of energy flow based electrothermal coupled network carbon exhaust flow metering according to any one of claims 1 to 7.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202211026388.3A CN115329518A (en) | 2022-08-25 | 2022-08-25 | Energy flow-based electrothermal coupling network carbon emission flow metering method, electronic equipment and storage medium |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202211026388.3A CN115329518A (en) | 2022-08-25 | 2022-08-25 | Energy flow-based electrothermal coupling network carbon emission flow metering method, electronic equipment and storage medium |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN115329518A true CN115329518A (en) | 2022-11-11 |
Family
ID=83927398
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202211026388.3A Pending CN115329518A (en) | 2022-08-25 | 2022-08-25 | Energy flow-based electrothermal coupling network carbon emission flow metering method, electronic equipment and storage medium |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN115329518A (en) |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN117112978A (en) * | 2023-10-19 | 2023-11-24 | 华北电力大学 | Method, equipment and medium for tracking carbon flow of electric power system |
-
2022
- 2022-08-25 CN CN202211026388.3A patent/CN115329518A/en active Pending
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN117112978A (en) * | 2023-10-19 | 2023-11-24 | 华北电力大学 | Method, equipment and medium for tracking carbon flow of electric power system |
| CN117112978B (en) * | 2023-10-19 | 2024-01-09 | 华北电力大学 | Method, equipment and medium for tracking carbon flow of electric power system |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Zhou et al. | Integrated power and heat dispatch considering available reserve of combined heat and power units | |
| Chen et al. | Increasing the flexibility of combined heat and power for wind power integration in China: Modeling and implications | |
| Lu et al. | Optimization model for the short-term joint operation of a grid-connected wind-photovoltaic-hydro hybrid energy system with cascade hydropower plants | |
| CN106777708B (en) | Steady state analysis method of electric power-natural gas regional comprehensive energy system | |
| CN109711601A (en) | The hot integrated energy system distributed optimization dispatching method of electric-gas-and device | |
| Bai et al. | Interdependence of electricity and heat distribution systems coupled by an AA‐CAES‐based energy hub | |
| CN109272353B (en) | System dynamic probability energy flow analysis method considering comprehensive demand response uncertainty | |
| CN115329518A (en) | Energy flow-based electrothermal coupling network carbon emission flow metering method, electronic equipment and storage medium | |
| Zhao et al. | Assessment of distributed photovoltaic hosting capacity in integrated electricity and heat systems considering uncertainty | |
| CN116822683A (en) | Comprehensive energy system optimization operation method based on carbon transaction meter and uncertainty | |
| Qin et al. | An adaptive piecewise linearized weighted directed graph for the modeling and operational optimization of integrated energy systems | |
| CN113268698B (en) | New energy consumption capacity optimization method and device of comprehensive energy system and terminal equipment | |
| Wang et al. | Integrated modeling of regional and park-level multi-heterogeneous energy systems | |
| Li et al. | Planning model of integrated energy system considering P2G and energy storage | |
| CN110457745B (en) | Method for establishing water-energy contact network model with close coupling relation | |
| Sun et al. | Multi-energy flow calculation method for we-energy based energy internet | |
| CN111724026A (en) | Optimization method for coupling operation of multi-energy network and water distribution network | |
| Liu et al. | Influence Evaluation of Integrated Energy System on the Unit Commitment in Power System | |
| CN115906456A (en) | Hydrogen-containing energy IES scheduling optimization model considering response uncertainty of demand side | |
| CN111695235B (en) | Regional decomposition method, device and system for renewable energy source absorption responsibility weight | |
| Tang et al. | Benefits allocation method for demand response in smart distribution network with high PV penetration | |
| CN112418732A (en) | Multi-energy hub-containing comprehensive energy system planning method and system based on maximum flow and minimum cut theorem | |
| Chen et al. | Multi-energy flow calculation considering the convexification network constraints for the integrated energy system | |
| CN113131513B (en) | Method for optimizing operation of electric, thermal and gas conversion system with consideration of carbon emission and storage medium | |
| Chen et al. | A joint solution framework for the optimal operation and evaluation of integrated energy systems based on exergy |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination |