Method and device for calculating electric power carbon emission of areas and enterprises under full coverage of green evidence

Abstract

The invention provides a method and a device for calculating electric power carbon emission of a region and an enterprise under full coverage of green certificates. The method comprises the following steps: acquiring a power flow matrix and a transaction contract of a target power grid; contract unit nodes and contract load nodes involved in the transaction contract in the target power grid are determined, and a transaction flow branch load flow distribution matrix is constructed; subtracting the transaction flow branch load flow distribution matrix from the load flow matrix; determining a carbon emission factor of a contract unit node and a contract electric quantity of a contract load node in the transaction flow branch power flow distribution matrix, and multiplying the carbon emission factor by the contract electric quantity; determining the carbon potential of a unit node and the natural power consumption of a load node in the natural distribution branch load flow distribution matrix, and multiplying the carbon potential of the load node by the natural power consumption; and adding the carbon emissions belonging to the same load node to obtain the carbon emissions of the load nodes in the target power grid. Therefore, the carbon emission of each load node is calculated according to the actual power consumption condition, and the accuracy of carbon emission calculation is improved.

Description

TECHNICAL FIELD

[0001] The present application relates to the technical field of carbon emissions, and particularly relates to a method and device for calculating regional and enterprise electric power carbon emissions under green certificate full coverage, computer equipment, a computer readable storage medium, and a computer program product. BACKGROUND

[0002] Calculating enterprise carbon emissions is a core step for quantifying the influence of human activities on climate change. Through accurate measurement, scientific emission reduction plans can be developed based on data, and technology upgrading and energy efficiency improvement can be promoted.

[0003] At present, the carbon emissions of enterprises in the power grid are mainly calculated by the formula: enterprise electricity carbon emissions = electricity consumption x power grid carbon emission factor. The power grid carbon emission factor is the key to the entire calculation. The power grid carbon emission factor represents the average amount of carbon dioxide emissions per unit of electricity consumed. Generally, the power grid carbon emission factor can be obtained from the official data released by official authoritative agencies every year.

[0004] However, the power grid carbon emission factor obtained in this way is an average value. Using the same power grid carbon emission factor to calculate the carbon emissions of enterprises of different electricity types cannot achieve accurate calculation of the carbon emissions of individual enterprises, and thus reduces the accuracy of enterprise carbon emission calculation. SUMMARY

[0005] The purpose of the embodiments of the present application is to provide a method and device for calculating regional and enterprise electric power carbon emissions under green certificate full coverage, computer equipment, a computer readable storage medium, and a computer program product, so as to improve the accuracy of enterprise carbon emission calculation.

[0006] To solve the above technical problems, the embodiments of the present application provide the following technical solutions: The first aspect of the application provides a method for calculating regional and enterprise power carbon emissions under green certificate full coverage, which comprises: obtaining a power flow matrix and a transaction contract of a target power grid, wherein the target power grid at least includes unit nodes and load nodes, the power flow matrix is used to represent the power flow between nodes in the target power grid, and the transaction contract is used to represent the power transaction between the unit nodes and the load nodes in the target power grid; determining the contract unit nodes and the contract load nodes involved in the transaction contract in the target power grid, and constructing a transaction flow branch power distribution matrix based on the connection relationship between the contract unit nodes and the contract load nodes in the target power grid and the power flow matrix; subtracting the transaction flow branch power distribution matrix from the power flow matrix to obtain a natural distribution branch power distribution matrix; determining the carbon emission factor of the contract unit node and the contract power of the contract load node in the transaction flow branch power distribution matrix, and multiplying the carbon emission factor by the contract power to obtain the contract carbon emission of the contract load node in the transaction flow branch power distribution matrix, wherein the carbon emission factor is used to represent the carbon emission generated by the unit node producing unit power; determining the carbon potential of the unit node and the natural power consumption of the load node in the natural distribution branch power distribution matrix, and determining the carbon potential of the load node based on the carbon potential of the unit node according to the node internal carbon concentration proportion sharing principle, and multiplying the carbon potential of the load node by the natural power consumption to obtain the natural carbon emission of the load node in the natural distribution branch power distribution matrix, wherein the carbon potential is used to represent the carbon emission generated by the corresponding node producing unit power; adding the carbon emission of the contract load node in the transaction flow branch power distribution matrix and the carbon emission of the load node in the natural distribution branch power distribution matrix belonging to the same load node to obtain the carbon emission of each load node in the target power grid.

[0007] Compared with the prior art, the method for calculating regional and enterprise power carbon emissions under green certificate full coverage provided by the first aspect of the application distinguishes the power involved in the transaction contract and the natural flow in the power grid, generates the corresponding branch power distribution matrix, respectively calculates the carbon emission of each load node under the transaction contract and the natural flow by multiplying the power consumption by the carbon emission factor on the power generation side and the carbon potential recursion and multiplying the carbon potential by the power consumption at the load node, so that the carbon emission of each load node is calculated according to the actual power consumption, thereby improving the accuracy of the calculation of the carbon emission of each enterprise in the power grid.

[0008] In other embodiments provided in the application, before constructing the transaction flow branch power distribution matrix based on the connection relationship between the contract generating unit nodes and the contract load nodes in the target power grid and the power flow matrix, the method further comprises: generating a branch network loss matrix based on the power flow matrix, the branch network loss matrix being used to represent the power loss of the branches between nodes in the target power grid; subtracting the power loss of the corresponding branch in the branch network loss matrix from the power of the node in the power flow matrix to obtain a lossless power flow matrix of the branches without loss in the ideal state of the target power grid; and constructing the transaction flow branch power distribution matrix based on the connection relationship between the contract generating unit nodes and the contract load nodes in the target power grid and the power flow matrix, comprising: constructing the transaction flow branch power distribution matrix based on the connection relationship between the contract generating unit nodes and the contract load nodes in the target power grid and the lossless power flow matrix; the method further comprises: taking the branch power loss in the branch network loss matrix as the generated power of the branch starting node; and distributing the generated power of each node to each load node according to the proportion of the power consumption of each load node in the target power grid.

[0009] Considering the branch network loss in the power grid and constructing a virtual load at the branch starting node, the power of the virtual load can be distributed to the corresponding load node, which can fairly distribute the branch network loss to the corresponding load node for carbon emission calculation, further improving the accuracy of carbon emission calculation of each enterprise in the power grid.

[0010] In other embodiments provided in the application, before distributing the generated power of each node to each load node according to the proportion of the power consumption of each load node in the target power grid, the method further comprises: determining a first branch network loss matrix corresponding to the transaction flow branch power distribution matrix based on the power flow matrix and the transaction flow branch power distribution matrix; subtracting the first branch network loss matrix from the branch network loss matrix to obtain a second branch network loss matrix corresponding to the natural distribution branch power distribution matrix; and distributing the generated power of each node to each load node according to the proportion of the power consumption of each load node in the target power grid, comprising: distributing the generated power of each node corresponding to the first branch network loss matrix to the load node corresponding to the transaction contract; and distributing the generated power of each node corresponding to the second branch network loss matrix to the corresponding load node according to the proportion of the natural power consumption of each load node.

[0011] Distributing the network loss based on the transaction contract first and then distributing the remaining network loss according to the power consumption proportion of the load node can realize correct distribution of the network loss, further improving the accuracy of carbon emission calculation of each enterprise in the power grid.

[0012] In other embodiments provided in this application, the transaction contract also includes a green certificate transaction contract for the purchase of green certificates between various generating nodes; the method further includes: determining the purchasing node and the selling node from each generating node based on the green certificate transaction contract; determining the carbon potential of the generating nodes and the natural electricity consumption of the load nodes in the natural distribution branch power flow distribution matrix, including: determining the original carbon emission factors of the purchasing node and the selling node in the generating nodes in the natural distribution branch power flow distribution matrix; diluting the original carbon emission factor corresponding to the purchasing node using the green electricity purchased by the purchasing node in the green certificate transaction contract to obtain the new carbon emission factor of the purchasing node; and concentrating the original carbon emission factor corresponding to the selling node using the green electricity sold by the selling node in the green certificate transaction contract to obtain the new carbon emission factor of the selling node; and determining the new carbon emission factors corresponding to the purchasing node and the selling node respectively as the carbon emission factors of the generating nodes in the natural distribution branch power flow distribution matrix.

[0013] In cases where generating units purchase green certificates from each other, the carbon potential of the generating unit purchasing green certificates is diluted according to the amount purchased, while the carbon potential of the generating unit selling green certificates is concentrated according to the amount sold. This approach fully considers the actual trading situation between generating units, thereby improving the accuracy of determining the carbon potential of each generating unit and thus improving the accuracy of carbon emission calculations for each enterprise in the power grid.

[0014] In other embodiments provided in this application, the transaction contract also includes a green certificate transaction contract for load nodes to purchase green certificates; determining the carbon potential of the generating nodes and the natural electricity consumption of the load nodes in the natural distribution branch power flow distribution matrix includes: determining the original natural electricity consumption of the load nodes in the natural distribution branch power flow distribution matrix; subtracting the purchased electricity of the corresponding load node in the green certificate transaction contract from the original natural electricity consumption to obtain the natural electricity consumption of the load nodes in the natural distribution branch power flow distribution matrix.

[0015] When load nodes purchase green certificates, the electricity consumption of the load nodes that purchase green certificates is offset, which can fully take into account the green certificate trading situation of load nodes, more accurately measure the carbon emissions of load nodes, and thus improve the accuracy of carbon emission calculations for various enterprises in the power grid.

[0016] The second aspect of this application provides a regional and enterprise electricity carbon emission calculation device under full green certificate coverage. The device includes: an acquisition module for acquiring the power flow matrix and trading contracts of a target power grid, wherein the target power grid includes at least generating unit nodes and load nodes, the power flow matrix is ​​used to characterize the power flow between nodes in the target power grid, and the trading contracts are used to characterize the power trading between generating unit nodes and load nodes in the target power grid; a determination module for determining the contract generating unit nodes and contract load nodes involved in the trading contracts in the target power grid, and constructing a power flow distribution matrix of trading flow branches based on the connection relationship of the contract generating unit nodes and contract load nodes in the target power grid and the power flow matrix; a subtraction module for subtracting the power flow distribution matrix of trading flow branches from the power flow matrix to obtain a power flow distribution matrix of naturally distributed branches; and a contract module for determining the carbon emission factor of the contract generating unit nodes and the contract electricity of the contract load nodes in the power flow distribution matrix of trading flow branches. The system calculates the carbon emissions of each load node in the power flow distribution matrix. The carbon emission factor is multiplied by the contracted electricity volume to obtain the contracted carbon emissions of the load nodes in the power flow distribution matrix. The carbon emission factor characterizes the carbon emissions generated per unit of electricity produced by the generating nodes. The system also includes a natural module to determine the carbon potential of generating nodes and the natural electricity consumption of load nodes in the power flow distribution matrix. Based on the carbon potential of generating nodes and the principle of sharing carbon concentration ratios within nodes, the system determines the carbon potential of load nodes and multiplies it by the natural electricity consumption to obtain the natural carbon emissions of load nodes in the power flow distribution matrix. The carbon potential characterizes the carbon emissions generated per unit of electricity produced by the corresponding node. Finally, the system calculates the carbon emissions of each load node in the target power grid by adding the contracted carbon emissions of the load nodes in the power flow distribution matrix to the carbon emissions of the load nodes belonging to the same load node within the natural carbon emissions of the load nodes in the power flow distribution matrix.

[0017] A third aspect of this application provides a computer device including a memory, a processor, and a computer program stored in the memory, wherein the processor executes the computer program to implement the method of the first aspect.

[0018] A fourth aspect of this application provides a computer-readable storage medium having a computer program stored thereon, which, when executed by a processor, implements the method of the first aspect.

[0019] The fifth aspect of this application provides a computer program product, including a computer program that, when executed by a processor, implements the method of the first aspect.

[0020] The regional and enterprise electricity carbon emission calculation device under full green certificate coverage provided in the second aspect of this application, the computer equipment provided in the third aspect, the computer-readable storage medium provided in the fourth aspect, and the computer program product provided in the fifth aspect have the same or similar beneficial effects as the regional and enterprise electricity carbon emission calculation method under full green certificate coverage provided in the first aspect. Attached Figure Description

[0021] The above and other objects, features, and advantages of exemplary embodiments of this application will become readily understood by reading the following detailed description with reference to the accompanying drawings. In the drawings, several embodiments of this application are illustrated by way of example and not limitation, with the same or corresponding reference numerals denoteing the same or corresponding parts, wherein: Figure 1 This is a flowchart illustrating the method for calculating regional and enterprise electricity carbon emissions under full green certificate coverage in this application embodiment. Figure 1 ; Figure 2 This is a flowchart illustrating the method for calculating regional and enterprise electricity carbon emissions under full green certificate coverage in this application embodiment. Figure 2 ; Figure 3 This is a schematic diagram of the structure of the regional and enterprise electricity carbon emission calculation device under the full coverage of green certificates in this application embodiment. Figure 1 ; Figure 4 This is a schematic diagram of the structure of the regional and enterprise electricity carbon emission calculation device under the full coverage of green certificates in this application embodiment. Figure 2 ; Figure 5 This is a schematic diagram of the structure of the computer device in the embodiments of this application. Detailed Implementation

[0022] Exemplary embodiments of this application will now be described in more detail with reference to the accompanying drawings. While exemplary embodiments of this application are shown in the drawings, it should be understood that this application may be implemented in various forms and should not be limited to the embodiments set forth herein. Rather, these embodiments are provided to enable a more thorough understanding of this application and to fully convey the scope of this application to those skilled in the art.

[0023] It should be noted that, unless otherwise stated, the technical or scientific terms used in this application shall have the ordinary meaning as understood by one of ordinary skill in the art to which this application pertains.

[0024] Currently, a uniform carbon emission factor is used to calculate the carbon emissions of various enterprises in the power grid. This is unfair to enterprises with different electricity consumption types and different positions in the power grid, thereby reducing the accuracy of carbon emission calculations for each enterprise in the power grid.

[0025] In view of this, embodiments of this application provide a method for calculating regional and enterprise electricity carbon emissions under full green certificate coverage, a device for calculating regional and enterprise electricity carbon emissions under full green certificate coverage, computer equipment, computer-readable storage medium, and computer program products. Through transaction contracts involving green electricity and green certificates, the power flow matrix of the power grid is separated according to the presence of transaction contracts and natural electricity consumption. Then, under different matrices, different methods are used to calculate the carbon emissions of corresponding load nodes, and the carbon emissions are then combined based on the same load node. In this way, the carbon emission factor calculated based on carbon emissions is sufficiently subdivided within the power grid. Using different carbon emission factors to calculate the carbon emissions of corresponding load nodes can improve the accuracy of carbon emission calculations for each load node, thereby improving the accuracy of carbon emission calculations for each enterprise in the power grid.

[0026] It should be noted that all components, data, and related processing methods involved in this application are authorized by the user or fully authorized by all parties, and the collection, use, and processing of related data comply with the relevant laws, regulations, and standards of the relevant countries and regions.

[0027] First, the calculation method for regional and enterprise electricity carbon emissions under the full coverage of green certificates provided in this application embodiment will be described in detail.

[0028] Figure 1 This is a flowchart illustrating the method for calculating regional and enterprise electricity carbon emissions under full green certificate coverage in this application embodiment. Figure 1 See Figure 1 As shown, the method may include: S11: Obtain the power flow matrix and transaction contracts of the target power grid, wherein the target power grid includes at least generating nodes and load nodes. The power flow matrix is ​​used to characterize the power flow between nodes in the target power grid, and the transaction contracts are used to characterize the power transaction between generating nodes and load nodes in the target power grid.

[0029] The target power grid refers to a specific power system that is abstracted into a network model for detailed calculation of electricity carbon emissions.

[0030] The target power grid can include generating nodes, intermediate nodes, load nodes, and their interconnections. Generating nodes can be various types of power plants, such as thermal power plants, wind farms, and photovoltaic power plants. Intermediate nodes are grid nodes in the target power grid that are neither sources of carbon emissions (generating nodes) nor the ultimate bearers of carbon emission responsibility (load nodes), such as 500kV / 220kV hub substations, 220kV / 110kV regional substations, and high-voltage direct current (HVDC) converter stations. Load nodes are electricity users, which can be a factory, a commercial building, a residential area, or an aggregated electricity sales company; they are the ultimate bearers of carbon emission responsibility.

[0031] A power flow matrix is ​​a mathematical square matrix used to precisely describe the net flow of active power through branches in a target power grid at a specific time (or time period). The power flow matrix is ​​an N×N square matrix, where N is the total number of nodes in the target power grid. Each element P_ij in the matrix represents the active power flowing from node i to node j.

[0032] The power flow matrix can be obtained directly from the state estimator of the power grid dispatch center, or it can be obtained by performing power flow calculations on the target power grid through a power flow calculation program.

[0033] A power purchase agreement (PPA) is a legally binding agreement reached between power generation companies (generator nodes) and power users (load nodes) or electricity retailers in the electricity market environment, regarding the quantity, price, time period, and energy attributes (such as whether it is green electricity) of the electricity traded. For example, load node L2 (data center) and generator node G2 (wind farm) signed a "May 2024 Green Electricity Purchase Agreement." The traded electricity volume is 5000 MWh, the trading period is from 00:00 on May 1, 2024 to 24:00 on May 31, 2024, the electricity type is green electricity, and the settlement point is node L2.

[0034] Trading contracts can be obtained directly from power trading centers, green certificate issuance and trading institutions, and power grid company marketing and settlement systems.

[0035] S12: Determine the contract generator nodes and contract load nodes involved in the transaction contracts in the target power grid, and construct the transaction flow branch power flow distribution matrix based on the connection relationship of the contract generator nodes and contract load nodes in the target power grid and the power flow matrix.

[0036] Determining the contracted generating unit nodes and contracted load nodes involves locating the generating unit nodes and load nodes involved in the transaction contract within the target power grid. Specifically, this can be done using node identifiers or a node identifier lookup table to map the nodes in the transaction contract to the nodes in the target power grid.

[0037] The power flow distribution matrix for trading branches is an NxN square matrix with the same dimensions as the global power flow matrix PB. It quantitatively describes, at the branch level, the virtual power flow caused by a single or all power trading contracts. It depicts a scenario where contracted power does not participate in the unified power flow distribution across the entire network, but rather flows "directly" from the contracted generator node to the contracted load node, like a "dedicated line." The matrix shows the power distribution along this "dedicated line" path. For each element PB_trade(i,j) in the matrix, it represents the active power flowing from node i to node j in the trading flow scenario.

[0038] When constructing the transaction flow branch power distribution matrix, the first step is to determine the contribution of each transaction to the branch power flow. Using the PTDF matrix, the power flow F_l, trade generated by a single transaction on each branch l is calculated: Fl, trade =(Φ l , s -Φ l , b )× Ts → b Where (Φ_l,s - Φ_l,b) is the sensitivity factor of the transaction pair branch l from node s to node b, that is, what proportion of the power supplied from s to b will flow through branch l. T_{s->b} is the power of the transaction.

[0039] Next, the power flow distribution matrix PB_trad for the transaction flow branches is constructed. This is a process of transforming the branch power flow F_l, trade back into the node-view matrix PB_trade. An N x N zero matrix PB_trade is initialized. For each branch l in the power grid, it connects nodes i and j. Based on the previously calculated F_l, trade: if F_l, trade > 0, it means the transaction flow is from node i to j, then: PB_trade(i, j) = PB_trade(i, j) + F_l, trade, PB_trade(j, i) = PB_trade(j, i) - F_l, trade (maintaining antisymmetry). If F_l, trade < 0, it means the transaction flow is from node j to i, then the opposite operation is performed. After traversing all branches, the PB_trade matrix corresponding to this transaction is obtained.

[0040] Finally, process multiple transactions. Calculate the transaction flow matrix PB_trade_k for each transaction contract k, and then superimpose them to obtain the total transaction flow branch distribution matrix PB_trade_total.

[0041] S13: Subtract the transaction flow branch flow distribution matrix from the flow matrix to obtain the natural distribution branch flow distribution matrix.

[0042] This involves directly performing matrix subtraction. S14: Determine the carbon emission factor of the contract unit node and the contract electricity of the contract load node in the transaction flow branch power distribution matrix, and multiply the carbon emission factor by the contract electricity to obtain the contract carbon emission of the contract load node in the transaction flow branch power distribution matrix. The carbon emission factor is used to characterize the carbon emission generated by the unit node producing a unit of electricity.

[0043] The core logic here is "whoever signs the contract bears the responsibility," meaning that load nodes are responsible for the carbon emissions of the electricity they purchase under contract, and the responsibility is directly traced back to the power source with which they signed the contract.

[0044] First, determine the carbon emission factor of the contracted generator node. Locate the carbon emission factor of the generator unit (contracted generator node) that has signed a contract with the load node. For fossil fuel units such as thermal power and gas-fired power: the carbon emission factor is based on parameters such as fuel type and power generation efficiency, using national standards or measured values. For example, the carbon emission factor of a typical coal-fired power unit might be between 0.8 tCO2 / MWh and 1.0 tCO2 / MWh. For renewable energy units such as wind power and solar power: since their power generation process does not directly generate carbon emissions, their carbon emission factor is 0 tCO2 / MWh. This is the fundamental reason why green electricity trading achieves carbon emission reduction.

[0045] Next, determine the contracted electricity volume of the load node. Find the actual electricity volume settled by the load node under this contract within the accounting period. This can be obtained directly from the settlement data of the power trading center. For example, a monthly green electricity contract with a settled electricity volume of 10,000 MWh.

[0046] Finally, calculate the contract carbon emissions. Multiply the two values ​​above. Contract carbon emissions = Carbon emission factor of contract generator nodes × Contract electricity consumption of contract load nodes.

[0047] For example, consider the following scenarios: Contract A (Green Electricity Trading): Contract generator node: G2, carbon emission factor: 0 tCO2 / MWh (because it is wind power), contract load node: L3, contract power volume: 5000 MWh (settlement power this month). Contract B (Conventional Thermal Power Trading): Contract generator node: G1, carbon emission factor: 0.9 tCO2 / MWh (typical coal power value), contract load node: L3, contract power volume: 2000 MWh (settlement power this month).

[0048] For Contract A (Green Electricity): Contract carbon emissions = G2 carbon emission factor × L3 contract electricity = 0 tCO2 / MWh × 5000 MWh = 0 tCO2.

[0049] For Contract B (thermal power): Contract carbon emissions = G1 carbon emission factor × L3 contract electricity = 0.9 tCO2 / MWh × 2000 MWh = 1800 tCO2 The total contract carbon emissions of load node L3 = Contract A carbon emissions + Contract B carbon emissions = 0 tCO2 + 1800 tCO2 = 1800 tCO2.

[0050] S15: Determine the carbon potential of the generating unit nodes and the natural electricity consumption of the load nodes in the natural distribution branch power flow distribution matrix. Based on the carbon potential of the generating unit nodes, determine the carbon potential of the load nodes according to the principle of sharing the carbon concentration ratio within the node. Multiply the carbon potential of the load nodes by the natural electricity consumption to obtain the natural carbon emissions of the load nodes in the natural distribution branch power flow distribution matrix. The carbon potential is used to characterize the carbon emissions generated by the corresponding node for producing a unit of electricity.

[0051] Specifically, first, determine the unit's carbon potential and the load's natural power consumption. The carbon potential of a unit node under naturally distributed power flow is an inherent property of the generating unit. For thermal power units, it is the carbon emission per unit of electricity generated (tCO2 / MWh); for wind power, photovoltaic, etc., it is 0. The natural power consumption of a load node can be obtained from the naturally distributed branch power flow distribution matrix. The natural power consumption of a load node is equal to the sum of all branch power flows flowing into that node.

[0052] Next, the nodal carbon potential is calculated. The "proportional sharing principle" needs to be applied. That is, assuming all electricity flowing into a node is completely mixed at that node, then any kilowatt-hour of electricity flowing out of that node has the same carbon footprint (i.e., carbon potential), equal to the weighted average of the carbon potentials of all incoming electricity. This is a recursive process from source to load. It must start from the generating node (where the carbon potential is known) and proceed step-by-step along the direction of the power flow to calculate the carbon potential of each downstream node.

[0053] Finally, the natural carbon emissions of the load nodes are calculated. For each load node l l Its natural carbon emissions are: Enatural , l = CEl × Pnatural , l .in, CEl It is the carbon potential at the loading node. Pnatural , l It is the natural electricity consumption of the load node.

[0054] For example, suppose the scenario is: Node 1: Thermal power unit G1, whose power generation carbon emission factor is... CEG 1 = 0.8 tCO2 / MWh. Node 2: Load L2. Branch: Connecting nodes 1 and 2. Power flowing from node 1 to node 2 is... PBnatural (1,2)=100 MW Next, we calculate the natural carbon emissions of load node L2.

[0055] Unit node carbon potential: CE 1 = 0.8 tCO2 / MWh.

[0056] Natural power consumption at load nodes: Pnatural ,2= PBnatural (1,2)=100 MWh (Assume the time period is 1 hour).

[0057] For node 2 (the load node), there is only one injection source, namely node 1. The carbon potential of node 2 is: CE 2= P 2, inPBnatural (1,2)× CE 1 = 100 × 0.8 = 0.8 tCO2 / MWh. (If node 2 has a next-level node, continue the calculation in this way.) Natural carbon emissions at load nodes: Enatural ,2= CE 2× Pnatural ,2=0.8×100=80 tCO2.

[0058] It should be noted that, in order to improve the efficiency of carbon emission calculation, steps S14 and S15 can be performed simultaneously.

[0059] S16: Add the contract carbon emissions of the contract load nodes in the transaction flow branch power flow distribution matrix to the carbon emissions of the load nodes in the natural distribution branch power flow distribution matrix that belong to the same load node, to obtain the carbon emissions of each load node in the target power grid.

[0060] In other words, for each unique load node, its carbon emission responsibility in two parallel worlds (the transaction flow world and the natural flow world) is merged.

[0061] First, input two independent carbon accounting results: List A (Transaction Flow Results): A list recording all load nodes participating in the transaction contracts and their corresponding contract carbon emissions. For example: [Load Node L1: 50 tCO2, Load Node L3: 0 tCO2]. List B (Natural Flow Results): A list recording all existing load nodes and their corresponding natural carbon emissions. For example: [Load Node L1: 30 tCO2, Load Node L2: 80 tCO2, Load Node L3: 20 tCO2]. Then, iterate through each load node in List B. For each node (e.g., L1), search for a record with the same node number in List A. If found, add the contract carbon emissions in List A to the natural carbon emissions in List B to obtain the node's total carbon emissions. If not found (e.g., L2), it means that the node has not participated in any trading contracts, and its total carbon emissions are equal to the natural carbon emissions in List B.

[0062] Finally, a complete inventory of total carbon emissions from all load nodes is output.

[0063] As can be seen from the above, the method for calculating regional and enterprise electricity carbon emissions under full green certificate coverage provided in this application distinguishes between electricity in the power grid involving trading contracts and natural flow, and generates corresponding branch power flow distribution matrices. The carbon emissions of each load node under trading contracts and natural flow are calculated by multiplying the carbon emission factor on the generation side by the electricity consumption measured by the power consumption, and by recursively extrapolating the carbon potential of the nodes, and multiplying the carbon potential by the electricity consumption at the load nodes. This ensures that the carbon emissions of each load node are calculated based on the actual electricity consumption, thereby improving the accuracy of carbon emission calculations for each enterprise in the power grid.

[0064] Furthermore, as Figure 1 In a refinement and extension of the method shown, this application also provides a method for calculating regional and enterprise electricity carbon emissions under full coverage of green certificates.

[0065] Figure 2 This is a flowchart illustrating the method for calculating regional and enterprise electricity carbon emissions under full green certificate coverage in this application embodiment. Figure 2 See Figure 2 As shown, the method may include: S21: Obtain the power flow matrix and transaction contracts of the target power grid, wherein the target power grid includes at least generating unit nodes and load nodes. The power flow matrix is ​​used to characterize the power flow between nodes in the target power grid, and the transaction contracts are used to characterize the power transaction between generating unit nodes and load nodes in the target power grid.

[0066] Step S21 here is implemented in the same way as step S11 in the previous embodiment, and can be found in the description in the previous embodiment. It will not be repeated here.

[0067] S22: Generate a branch network loss matrix based on the power flow matrix. The branch network loss matrix is ​​used to characterize the power loss of branches between nodes in the target power grid. Subtract the power loss of the corresponding branch in the branch network loss matrix from the power of the nodes in the power flow matrix to obtain the lossless power flow matrix with no branch loss under the ideal state of the target power grid.

[0068] In physically lossy networks, power flow is nonlinear. If the transaction flow is directly subtracted from the total power flow, the resulting "natural flow" may not satisfy Kirchhoff's laws and is an illegal, non-conservative power flow distribution that cannot be used for calculation.

[0069] By equating network losses to node loads, an "equivalent lossless network" is created. In this network, the power flow model can be approximated as a linear system, satisfying the superposition theorem. This makes it possible to precisely decompose the overall power flow into "transactional flow" and "natural flow".

[0070] The branch network loss matrix PBL is an N x N square matrix, where the element PBL(i,j) represents the active power loss on branch ij.

[0071] Network loss cannot be directly "generated" from the power flow matrix; it needs to be obtained through the results of power flow calculations.

[0072] In specific generation, for the branch connecting nodes i and j, its active power loss is... Ploss , ij This can be obtained from power flow calculations. In AC power flow, the formula for calculating network loss is: Ploss , ij = Iij 2× Rij .in, Iij It is the branch current. Rij This represents the branch resistance. Then, fill in the matrix. P_loss(i, j) = P_loss, ij (representing the loss associated with the power flow from i to j). P_loss(j, i) = P_loss, ij (network loss is a scalar and has no direction, so the same value is filled in symmetrical positions). Diagonal elements are 0.

[0073] Next, the network loss is represented as the additional load on each node. This is not matrix subtraction, but a redistribution process. In principle, the network loss of each branch is attached to the starting node of that branch as an additional "virtual load." For branch ij, its network loss is... Ploss , ijAdd an equivalent load to node i: Δ Loadi += Ploss , ij By traversing all branches, the total "equivalent network loss load" for each node k is calculated. Loadloss , k .

[0074] The lossless power flow matrix is ​​an N x N square matrix that describes the power flow distribution in a virtual power grid without branch losses.

[0075] In the specific calculation, the original lossy power flow matrix PB is used as the initial value for the lossless matrix PB_lossless, and then the power balance of all nodes is ensured. In most research and applications, when using a DC power flow model for this type of analysis, resistance is ignored by default, thus naturally resulting in a lossless power flow matrix.

[0076] S23: Determine the contract generator nodes and contract load nodes involved in the transaction contracts in the target power grid, and construct the transaction flow branch power flow distribution matrix based on the connection relationship of the contract generator nodes and contract load nodes in the target power grid and the lossless power flow matrix.

[0077] S24: Subtract the transaction flow branch flow distribution matrix from the lossless flow matrix to obtain the natural distribution branch flow distribution matrix.

[0078] The difference between this step and steps S12 and S13 in the previous embodiments is that the power flow matrix is ​​replaced with an ideal lossless power flow matrix. The specific methods for determining the contract unit nodes, contract load nodes, and the construction of the transaction flow branch power flow distribution matrix are the same as in steps S12 and S13 in the previous embodiments, and can be found in the relevant descriptions in the previous embodiments; they will not be repeated here.

[0079] S25: Use the branch power loss in the branch network loss matrix as the generated power of the branch starting node; according to the proportion of power consumption of each load node in the target power grid, allocate the generated power of each node to each load node.

[0080] In other words, the energy loss on each line is equivalent to the extra electricity generated at the "starting node" of the line. Then, based on the principle of "whoever uses the electricity, bears the cost," this loss is fairly distributed among the load nodes according to their actual electricity consumption ratio.

[0081] In active power flow, for a branch ij, if the power flow is from i to j, then i is the "starting node" of that branch. The power loss P_loss(i,j) on branch ij is entirely included in its starting node i, which is regarded as an additional virtual power generated by that node, specifically used to compensate for the loss on this line.

[0082] When allocating this portion of electricity, a power allocation principle can be adopted. Electricity generated by one node is transmitted along the power grid to various loads, and network losses are distributed according to the proportion of electricity each load node receives from the entire network. First, calculate the sum of the total active power consumption P_total_load of all load nodes in the target power grid. Then, for each load node l, calculate the proportion αl of its power consumption P_load,l to the total power consumption P_total_load. For the virtual generation G_loss,k of each node k, it is distributed to each load node l according to the proportion α_l.

[0083] Since both market transactions and natural use involve carbon emissions, in order to fairly distribute network losses and their carbon emissions among market transaction users and natural use users, the total loss can first be divided into "market transaction loss" and "natural use loss". Then, the "market transaction loss" is directly borne by the corresponding load node. The "natural use loss" is then shared by all using load nodes according to their usage ratio.

[0084] Specifically, step S25 above may include: S251: Based on the power flow matrix and the power flow distribution matrix of the transaction flow branches, determine the network loss matrix of the first branch corresponding to the power flow distribution matrix of the transaction flow branches.

[0085] For each branch l in the target power grid, calculate the proportion factor K_l of the transaction flow in the total power flow of that branch. Based on the calculated proportion factor K_l, calculate the transaction flow network loss for each branch l. Here, the data in each matrix can be regarded as the power in the power grid. After calculating the network loss of the transaction flow, it is converted into a matrix to obtain the network loss matrix of the first branch.

[0086] S252: Subtract the first branch network loss matrix from the branch network loss matrix to obtain the second branch network loss matrix corresponding to the natural distribution branch power flow distribution matrix.

[0087] First, confirm that the positions of elements in both matrices represent the same physical branch (e.g., (i,j) both represent the branch from node i to node j). Then, subtract element by element. Perform the subtraction operation once for each element (i.e., each branch) in the matrix. Finally, obtain the second branch network loss matrix.

[0088] S253: Distribute the generated electricity of each node corresponding to the first branch network loss matrix to the corresponding load node with a transaction contract.

[0089] In principle, whoever makes the transaction bears the responsibility. The network loss matrix of the first branch is equivalent to the virtual power generation of the nodes. The virtual power generation is then allocated to the corresponding contract load nodes.

[0090] S254: Distribute the generated electricity of each node corresponding to the second branch network loss matrix to the corresponding load nodes according to the proportion of the natural electricity consumption of each load node.

[0091] In principle, the losses are allocated fairly according to the proportion of electricity consumption. The second branch network loss matrix is ​​equivalent to the virtual generation of nodes. The natural electricity consumption ratio of each load node is calculated. The virtual generation is allocated proportionally. Ultimately, each load node (regardless of whether it participates in the transaction) receives a portion of the "natural flow network loss electricity".

[0092] Once the ideal lossless power flow matrix has been obtained and each node has borne its corresponding losses, we can continue to trade the flow and the naturally distributed carbon emissions separately.

[0093] S26: Determine the carbon emission factor of the contract unit node and the contract electricity of the contract load node in the transaction flow branch power distribution matrix, and multiply the carbon emission factor by the contract electricity to obtain the contract carbon emission of the contract load node in the transaction flow branch power distribution matrix. The carbon emission factor is used to characterize the carbon emission generated by the unit node producing a unit of electricity.

[0094] S27: Determine the carbon potential of the generating unit nodes and the natural electricity consumption of the load nodes in the natural distribution branch power flow distribution matrix. Based on the carbon potential of the generating unit nodes, determine the carbon potential of the load nodes according to the principle of sharing the carbon concentration ratio within the node. Multiply the carbon potential of the load nodes by the natural electricity consumption to obtain the natural carbon emissions of the load nodes in the natural distribution branch power flow distribution matrix. The carbon potential is used to characterize the carbon emissions generated by the corresponding node for producing a unit of electricity.

[0095] Under the "certificate-electricity separation" green certificate trading system, without proper handling, both the wind farm selling green certificates and the thermal power plant buying them might simultaneously claim ownership of the zero-carbon benefits of 100 MWh of green electricity. The wind farm would believe it produced green electricity, while the thermal power plant would believe it purchased environmental rights. This would result in the same low-carbon benefit being calculated twice.

[0096] Therefore, through a "dilution" and "concentration" mechanism, it is clearly stipulated that once a wind farm (the selling node) sells its green certificate, the corresponding portion of its physical electricity loses its zero-carbon attribute (is "concentrated"). The corresponding zero-carbon attribute is transferred to the thermal power plant (the purchasing node), causing a portion of its high-carbon electricity to be "greened" ("diluted"). This ensures that each environmental right is calculated only once within an accounting cycle.

[0097] On the one hand, the trading contracts also include green certificate trading contracts between various generating unit nodes for the purchase of green certificates. This is a special type of transaction that occurs between power plants, namely green electricity certificate (green certificate) trading contracts.

[0098] Specifically, determining the carbon potential of the unit node in step S27 above (the natural power consumption of the load node has now been factored in with branch power losses) may include: Step A1: Determine the buying and selling nodes from each generator group node based on the green certificate trading contract.

[0099] The green certificate trading contract stores the identifiers of each power plant. In the target power grid, each generating unit node also has a corresponding identifier. The power plants in the green certificate trading contract and some generating unit nodes in the target power grid belong to the same object. Therefore, a mapping relationship between power plant identifiers and generating unit node identifiers is established in advance. Then, the power plant identifiers and the buy / sell relationship are extracted from the green certificate trading contract. Finally, the buying and selling nodes in the target power grid are determined from the mapping relationship based on the power plant identifiers and the buy / sell relationship.

[0100] Step A2: Determine the original carbon emission factors of the purchasing and selling nodes in the unit nodes of the natural distribution branch power flow distribution matrix.

[0101] Specifically, it can be obtained through calculation: Raw Emission Factor (EF) = (Fuel consumption × Carbon content per unit calorific value of fuel × Carbon oxidation rate × 44 / 12) / Electricity generation. Alternatively, it can be obtained directly from officially published data.

[0102] Step A3: Dilute the original carbon emission factor corresponding to the purchasing node by using the green electricity purchased by the purchasing node in the green certificate trading contract to obtain the new carbon emission factor of the purchasing node; and concentrate the original carbon emission factor corresponding to the selling node by using the green electricity sold by the selling node in the green certificate trading contract to obtain the new carbon emission factor of the selling node.

[0103] In the dilution calculation, the purchased "zero-carbon" electricity (green certificate electricity) is virtually mixed into the high-carbon electricity, thereby lowering the overall average carbon concentration. New carbon potential of the purchasing node = (original carbon potential of the purchasing node × total power generation of the purchasing node + 0 × traded green electricity) / (total power generation of the purchasing node + traded green electricity).

[0104] In the enrichment calculation, the "zero-carbon" electricity sold by the node is virtually replaced with the system's average high-carbon electricity, thereby increasing the average carbon concentration of its remaining electricity. New carbon potential of the selling node = (Original carbon potential of the selling node × Total electricity generation of the selling node + System default carbon potential × Traded green electricity) / (Total electricity generation of the selling node + Traded green electricity).

[0105] Step A4: Determine the new carbon emission factors corresponding to the purchasing node and the selling node as the carbon emission factors of the unit nodes in the natural distribution branch power flow distribution matrix.

[0106] It should be noted that steps A1-A4 above also apply to step S26 above.

[0107] On the other hand, the trading contracts also include green certificate trading contracts where load nodes purchase green certificates. That is, electricity users (load nodes) directly purchase green electricity certificates (green certificates) from renewable energy generators (generator nodes) without the direct trading of physical electricity.

[0108] Specifically, determining the natural electricity consumption of the load node in step S27 above may include: Step B1: Determine the original natural power consumption of the load nodes in the natural distribution branch power flow distribution matrix.

[0109] The original natural electricity consumption is a value calculated from the natural distribution branch power flow distribution matrix. It represents the actual electricity consumed by a power user (load node) purely through the natural distribution of the physical power grid after removing the influence of all market transaction contracts.

[0110] Step B2: Subtract the purchased electricity of the corresponding load node in the green certificate trading contract from the original natural electricity consumption to obtain the natural electricity consumption of the load node in the natural distribution branch power flow distribution matrix.

[0111] The electricity purchased by the load node in the green certificate trading contract refers to the electricity corresponding to the green electricity certificate (green certificate) purchased by a specific power user (load node) through the official green certificate trading platform.

[0112] The physical electrical energy consumed by users is still hybrid electricity, which cannot be changed. However, "despite physically using so much electricity, because 'zero carbon certificates' (green certificates) have been purchased for a portion of it, when calculating carbon emissions, only the remaining 'uncertified' electricity is accounted for." In the target power grid, the carbon emission factors of each generating unit node have been diluted and enriched accordingly, and the electricity consumption of each load node has also been adjusted by adding branch losses and deducting based on green certificates. The carbon emissions of each load node can then be calculated separately for transaction flows and natural distribution.

[0113] S28: Add the contract carbon emissions of the contract load nodes in the transaction flow branch power flow distribution matrix to the carbon emissions of the load nodes in the natural distribution branch power flow distribution matrix that belong to the same load node, to obtain the carbon emissions of each load node in the target power grid.

[0114] Step S28 here is implemented in the same way as step S16 in the previous embodiment, and can be found in the description in the previous embodiment, which will not be repeated here.

[0115] This concludes the explanation of the regional and enterprise electricity carbon emission calculation methods under the full coverage of green certificates provided in this application.

[0116] Based on the same inventive concept, this application also provides a regional and enterprise electricity carbon emission calculation device under full coverage of green certificates.

[0117] Figure 3 This is a schematic diagram of the structure of the regional and enterprise electricity carbon emission calculation device under the full coverage of green certificates in this application embodiment. Figure 1 See Figure 3 As shown, the device may include: The acquisition module 31 is used to acquire the power flow matrix and transaction contracts of the target power grid. The target power grid includes at least generating nodes and load nodes. The power flow matrix is ​​used to characterize the power flow between nodes in the target power grid, and the transaction contracts are used to characterize the power transaction between generating nodes and load nodes in the target power grid.

[0118] The determination module 32 is used to determine the contract generator nodes and contract load nodes involved in the transaction contract in the target power grid, and construct the transaction flow branch power flow distribution matrix based on the connection relationship of the contract generator nodes and contract load nodes in the target power grid and the power flow matrix.

[0119] Subtraction module 33 is used to subtract the transaction flow branch flow distribution matrix from the flow matrix to obtain the natural distribution branch flow distribution matrix.

[0120] Contract module 34 is used to determine the carbon emission factor of the contract unit node and the contract electricity of the contract load node in the transaction flow branch power distribution matrix, and multiply the carbon emission factor by the contract electricity to obtain the contract carbon emission of the contract load node in the transaction flow branch power distribution matrix. The carbon emission factor is used to characterize the carbon emission generated by the unit node producing a unit of electricity.

[0121] Natural module 35 is used to determine the carbon potential of the unit nodes and the natural electricity consumption of the load nodes in the natural distribution branch power flow distribution matrix. Based on the carbon potential of the unit nodes, the carbon potential of the load nodes is determined according to the principle of sharing the carbon concentration ratio within the node. The carbon potential of the load nodes is multiplied by the natural electricity consumption to obtain the natural carbon emissions of the load nodes in the natural distribution branch power flow distribution matrix. The carbon potential is used to characterize the carbon emissions generated by the corresponding node for producing a unit of electricity.

[0122] The calculation module 36 is used to add the contract carbon emissions of the contract load nodes in the transaction flow branch power flow distribution matrix to the carbon emissions of the load nodes in the natural distribution branch power flow distribution matrix that belong to the same load node, so as to obtain the carbon emissions of each load node in the target power grid.

[0123] Furthermore, as Figure 3 In addition to the refinement and expansion of the device shown, this application also provides a regional and enterprise electricity carbon emission calculation device under full coverage of green certificates.

[0124] Figure 4 This is a schematic diagram of the structure of the regional and enterprise electricity carbon emission calculation device under the full coverage of green certificates in this application embodiment. Figure 2 See Figure 4 As shown, the device may include: The acquisition module 41 is used to acquire the power flow matrix and transaction contracts of the target power grid. The target power grid includes at least generating nodes and load nodes. The power flow matrix is ​​used to characterize the power flow between nodes in the target power grid, and the transaction contracts are used to characterize the power transaction between generating nodes and load nodes in the target power grid.

[0125] The network loss module 42 is used to generate a branch network loss matrix based on the power flow matrix. The branch network loss matrix is ​​used to characterize the power loss of branches between nodes in the target power grid. The power loss of the corresponding branch in the branch network loss matrix is ​​subtracted from the power of the nodes in the power flow matrix to obtain the lossless power flow matrix with no branch loss under the ideal state of the target power grid.

[0126] The network loss module 42 is also used to take the branch power loss in the branch network loss matrix as the generated power of the branch starting node; and to allocate the generated power of each node to each load node according to the proportion of power consumption of each load node in the target power grid.

[0127] The network loss module 42 is specifically used to determine the first branch network loss matrix corresponding to the transaction flow branch flow distribution matrix based on the power flow matrix and the transaction flow branch power flow distribution matrix; subtract the first branch network loss matrix from the branch network loss matrix to obtain the second branch network loss matrix corresponding to the natural distribution branch power flow distribution matrix; allocate the generated electricity of each node corresponding to the first branch network loss matrix to the corresponding load node with a transaction contract; and allocate the generated electricity of each node corresponding to the second branch network loss matrix to the corresponding load node according to the proportion of the natural power consumption of each load node.

[0128] The determination module 43 is used to determine the contract generator nodes and contract load nodes involved in the transaction contract in the target power grid, and to construct the transaction flow branch power flow distribution matrix based on the connection relationship of the contract generator nodes and contract load nodes in the target power grid and the lossless power flow matrix.

[0129] Subtraction module 44 is used to subtract the transaction flow branch flow distribution matrix from the lossless flow matrix to obtain the natural distribution branch flow distribution matrix.

[0130] Contract module 45 is used to determine the carbon emission factor of the contract unit node and the contract electricity of the contract load node in the transaction flow branch power distribution matrix, and multiply the carbon emission factor by the contract electricity to obtain the contract carbon emission of the contract load node in the transaction flow branch power distribution matrix. The carbon emission factor is used to characterize the carbon emission generated by the unit node producing a unit of electricity.

[0131] Natural module 46 is used to determine the carbon potential of the unit nodes and the natural electricity consumption of the load nodes in the natural distribution branch power flow distribution matrix. Based on the carbon potential of the unit nodes, the carbon potential of the load nodes is determined according to the principle of sharing the carbon concentration ratio within the node. The carbon potential of the load nodes is multiplied by the natural electricity consumption to obtain the natural carbon emissions of the load nodes in the natural distribution branch power flow distribution matrix. The carbon potential is used to characterize the carbon emissions generated by the corresponding node for producing a unit of electricity.

[0132] The calculation module 47 is used to add the contract carbon emissions of the contract load nodes in the transaction flow branch power flow distribution matrix to the carbon emissions of the load nodes in the natural distribution branch power flow distribution matrix that belong to the same load node, so as to obtain the carbon emissions of each load node in the target power grid.

[0133] In cases where the transaction contract also includes a green certificate transaction contract for the purchase of green certificates between various unit nodes, the calculation module 47 is specifically used to determine the purchasing node and the selling node from each unit node based on the green certificate transaction contract; determine the original carbon emission factors of the purchasing node and the selling node in the unit node in the natural distribution branch power flow distribution matrix; dilute the original carbon emission factor corresponding to the purchasing node using the green electricity purchased by the purchasing node in the green certificate transaction contract to obtain the new carbon emission factor of the purchasing node; and concentrate the original carbon emission factor corresponding to the selling node using the green electricity sold by the selling node in the green certificate transaction contract to obtain the new carbon emission factor of the selling node; and determine the new carbon emission factors corresponding to the purchasing node and the selling node respectively as the carbon emission factors of the unit nodes in the natural distribution branch power flow distribution matrix.

[0134] In the case where the transaction contract also includes a green certificate transaction contract for load nodes to purchase green certificates, the calculation module 47 is specifically used to determine the original natural electricity consumption of load nodes in the natural distribution branch power flow distribution matrix; and to subtract the purchased electricity of the corresponding load node in the green certificate transaction contract from the original natural electricity consumption to obtain the natural electricity consumption of load nodes in the natural distribution branch power flow distribution matrix.

[0135] It should be noted that the description of the above device embodiments is similar to the description of the above method embodiments, and has similar beneficial effects. For technical details not disclosed in the device embodiments of this application, please refer to the description of the method embodiments of this application for understanding.

[0136] Based on the same inventive concept, this application also provides a computer device.

[0137] Figure 5 This is a schematic diagram of the structure of the computer device in an embodiment of this application. See also... Figure 5 As shown, the computer device may include: a memory 51, a processor 52, and a computer program stored on the memory 51, wherein the processor 52 executes the computer program to implement the methods described in the foregoing embodiments.

[0138] It should be noted that the description of the above computer device embodiments is similar to the description of the above method embodiments, and has similar beneficial effects. For technical details not disclosed in the computer device embodiments of this application, please refer to the description of the method embodiments of this application for understanding.

[0139] Based on the same inventive concept, embodiments of this application also provide a computer-readable storage medium. The computer-readable storage medium stores a computer program, which, when executed by a processor, implements the methods described in the foregoing embodiments.

[0140] It should be noted that the description of the above computer-readable storage medium embodiments is similar to the description of the above method embodiments, and has similar beneficial effects. For technical details not disclosed in the computer-readable storage medium embodiments of this application, please refer to the description of the method embodiments of this application for understanding.

[0141] Based on the same inventive concept, this application also provides a computer program product. The computer program product includes a computer program that, when executed by a processor, implements the methods described in the foregoing embodiments.

[0142] It should be noted that the descriptions of the above computer program product embodiments are similar to those of the above method embodiments, and have similar beneficial effects. For technical details not disclosed in the computer program product embodiments of this application, please refer to the descriptions of the method embodiments of this application for understanding.

[0143] The above description is merely a specific embodiment of this application, but the scope of protection of this application is not limited thereto. Any variations or substitutions that can be easily conceived by those skilled in the art within the scope of the technology disclosed in this application should be included within the scope of protection of this application. Therefore, the scope of protection of this application should be determined by the scope of the claims.

Claims

1. A method for calculating regional and enterprise electric power carbon emissions under green certificate full coverage, characterized in that, The method comprises: obtaining a power flow matrix and a transaction contract of a target power grid, wherein the target power grid comprises at least a unit node and a load node, the power flow matrix is used to represent the power flow between nodes in the target power grid, and the transaction contract is used to represent the power transaction between the unit node and the load node in the target power grid; determining the contract unit node and the contract load node involved in the transaction contract in the target power grid, and constructing a transaction flow branch power distribution matrix based on the connection relationship between the contract unit node and the contract load node in the target power grid and the power flow matrix; subtracting the transaction flow branch power distribution matrix from the power flow matrix to obtain a natural distribution branch power distribution matrix; determining the carbon emission factor of the contract unit node and the contract power of the contract load node in the transaction flow branch power distribution matrix, and multiplying the carbon emission factor by the contract power to obtain the contract carbon emission of the contract load node in the transaction flow branch power distribution matrix, wherein the carbon emission factor is used to represent the carbon emission generated by the unit node in the production of unit power; determining the carbon potential of the unit node and the natural power consumption of the load node in the natural distribution branch power distribution matrix, and determining the carbon potential of the load node based on the carbon potential of the unit node according to the node internal carbon concentration proportion sharing principle, and multiplying the carbon potential of the load node by the natural power consumption to obtain the natural carbon emission of the load node in the natural distribution branch power distribution matrix, wherein the carbon potential is used to represent the carbon emission generated by the corresponding node in the production of unit power; adding the contract carbon emission of the contract load node in the transaction flow branch power distribution matrix to the carbon emission of the same load node in the natural distribution branch power distribution matrix to obtain the carbon emission of each load node in the target power grid.

2. The method of claim 1, wherein, Before constructing the transaction flow branch power distribution matrix based on the connection relationship between the contract unit node and the contract load node in the target power grid and the power flow matrix, the method further comprises: generating a branch network loss matrix based on the power flow matrix, wherein the branch network loss matrix is used to represent the power loss of the branch between nodes in the target power grid; subtracting the power loss of the corresponding branch in the branch network loss matrix from the power of the node in the power flow matrix to obtain a lossless power flow matrix of the branch without loss in the ideal state of the target power grid; constructing the transaction flow branch power distribution matrix based on the connection relationship between the contract unit node and the contract load node in the target power grid and the lossless power flow matrix; the method further comprises: regarding the branch power loss in the branch network loss matrix as the generated power of the branch starting node; allocating the generated power of each node to each load node according to the proportion of the power consumption of each load node in the target power grid. ​ 3. The method of claim 2, wherein, Before distributing the generated power of each node to each load node according to the proportion of the power consumption of each load node in the target power grid, the method further comprises: determining a first branch loss matrix corresponding to the transaction flow branch power distribution matrix based on the power flow matrix and the transaction flow branch power distribution matrix; subtracting the first branch loss matrix from the branch loss matrix to obtain a second branch loss matrix corresponding to the natural distribution branch power distribution matrix; the method further comprises: distributing the generated power of each node corresponding to the first branch loss matrix to the load nodes corresponding to the transaction contract; distributing the generated power of each node corresponding to the second branch loss matrix to the corresponding load nodes according to the proportion of the natural power consumption of each load node.

4. The method according to any one of claims 1 to 3, characterized in that, The transaction contract further comprises a green certificate transaction contract for purchasing green certificates among each unit node; the method further comprises: determining the purchase nodes and the sale nodes from each unit node based on the green certificate transaction contract; the method further comprises: determining the original carbon emission factors of the purchase nodes and the sale nodes in the natural distribution branch power distribution matrix; diluting the original carbon emission factor corresponding to the purchase node by using the green power quantity purchased by the purchase node in the green certificate transaction contract to obtain a new carbon emission factor of the purchase node; and concentrating the original carbon emission factor corresponding to the sale node by using the green power quantity sold by the sale node in the green certificate transaction contract to obtain a new carbon emission factor of the sale node; determining the new carbon emission factors corresponding to the purchase node and the sale node respectively as the carbon emission factors of the unit nodes in the natural distribution branch power distribution matrix.

5. The method according to any one of claims 1 to 3, characterized in that, The transaction contract further comprises a green certificate transaction contract for purchasing green certificates by load nodes; the method further comprises: determining the original natural power consumption of the load nodes in the natural distribution branch power distribution matrix; subtracting the purchase power of the corresponding load nodes in the green certificate transaction contract from the original natural power consumption to obtain the natural power consumption of the load nodes in the natural distribution branch power distribution matrix.

6. A device for calculating regional and enterprise electric power carbon emissions under green certificate full coverage, characterized in that, The device comprises: an acquisition module configured to acquire a power flow matrix and a transaction contract of a target power grid, wherein the target power grid comprises at least unit nodes and load nodes, the power flow matrix is used to represent the power flow among the nodes in the target power grid, and the transaction contract is used to represent the power transaction between the unit nodes and the load nodes in the target power grid; a determination module configured to determine contract unit nodes and contract load nodes involved in the transaction contract in the target power grid, and construct a transaction flow branch power distribution matrix based on the connection relationship between the contract unit nodes and the contract load nodes in the target power grid and the power flow matrix; a subtraction module, configured to subtract the transaction flow branch power flow distribution matrix from the power flow matrix to obtain a natural distribution branch power flow distribution matrix; a contract module, configured to determine a carbon emission factor of a contract unit node and a contract electricity of a contract load node in the transaction flow branch power flow distribution matrix, and multiply the carbon emission factor by the contract electricity to obtain a contract carbon emission of the contract load node in the transaction flow branch power flow distribution matrix, the carbon emission factor being used to represent carbon emission generated by a unit of electricity generated by the unit node; a natural module, configured to determine a carbon potential of a unit node and a natural electricity of a load node in the natural distribution branch power flow distribution matrix, and determine a carbon potential of the load node based on the carbon potential of the unit node according to a principle of sharing in proportion to carbon concentration in a node, and multiply the carbon potential of the load node by the natural electricity to obtain a natural carbon emission of the load node in the natural distribution branch power flow distribution matrix, the carbon potential being used to represent carbon emission generated by a unit of electricity generated by the corresponding node; a calculation module, configured to add the contract carbon emission of the contract load node in the transaction flow branch power flow distribution matrix and the natural carbon emission of the load node in the natural distribution branch power flow distribution matrix that belongs to the same load node to obtain carbon emissions of each load node in the target power grid.

7. The apparatus of claim 6, wherein, The device further comprises: a network loss module, configured to generate a branch network loss matrix based on the power flow matrix, the branch network loss matrix being used to represent electricity loss of a branch between nodes in the target power grid; and subtract the electricity loss of the corresponding branch in the branch network loss matrix from the electricity of the node in the power flow matrix to obtain a lossless power flow matrix without loss of the branch in an ideal state of the target power grid; the determination module is specifically configured to construct a transaction flow branch power flow distribution matrix based on a connection relationship between the contract unit node and the contract load node in the target power grid, and the lossless power flow matrix; the network loss module is further configured to take the electricity loss of the branch in the branch network loss matrix as generated electricity of a branch starting node; and distribute the generated electricity of each node to each load node according to a proportion of the electricity consumption of each load node in the target power grid.

8. A computer device comprising a memory, a processor, and a computer program stored on the memory, wherein the computer program comprises instructions that, when executed by the processor, cause the processor to perform the method of any one of claims 1-7. The processor executes the computer program to implement the steps of the method in any one of claims 1 to 5.

9. A computer readable storage medium having stored thereon a computer program, characterized in that, The computer program is executed by the processor to implement the steps of the method in any one of claims 1 to 5.

10. A computer program product comprising a computer program, characterized in that, The computer program is executed by the processor to implement the steps of the method in any one of claims 1 to 5. The computer program is executed by the processor to implement the steps of the method in any one of claims 1 to 5.

Images