Method and device for correcting node carbon emission factor of power system

Abstract

The invention discloses an electric power system node carbon emission factor correction method and device, and belongs to the technical field of electric power systems, and the method comprises the steps: collecting SCADA data of an electric power system in real time; collecting carbon emission data of generator nodes of the power system in real time; acquiring power grid node topological data of the power system; constructing a nonlinear measurement model based on SCADA data and power grid node topological data; solving the nonlinear measurement model to obtain voltage correction values of all nodes of the power system; calculating the complex power of each branch in the power system based on the voltage correction value; decomposing the complex power to obtain the active power correction of each branch; and based on the active power correction, the carbon emission data and the power grid node topological data, correcting a carbon emission factor of each node in the power system. Therefore, by implementing the method, the problem of low calculation accuracy of the carbon emission factor in the prior art can be solved.

Description

TECHNICAL FIELD

[0001] The present application relates to the technical field of power systems, in particular to a power system node carbon emission factor correction method and device. BACKGROUND

[0002] The current carbon emission factor calculation method includes the following several kinds: 1. Through the installation of sensors and the collection of carbon emission data of power users, and through the data law after preprocessing to predict the carbon emission factor. 2. Through the load sequence and the carbon emission intensity of the generator to calculate the corresponding node electric carbon factor, and use the graph attention network method to predict the electric carbon factor. 3. Through the coal quality parameters and the efficiency of the power plant to calculate the carbon emission factor. 4. Through the construction of multiple linear regression model and time series prediction model, using federal learning algorithm to protect data privacy, realizing the comprehensive prediction of carbon emission factor of remote data.

[0003] However, in practical application, the prior art has the problems of large amount of data information, too many required calculation conditions, etc. And due to cost factors, a large number of nodes and branches in the system are not installed with sensors, and the power flow distribution of the entire power system cannot be measured, and sometimes it may be difficult to measure, cannot be directly obtained or data packet loss, etc. Therefore, the carbon emission factor calculated by the existing method has low accuracy. SUMMARY

[0004] The present application provides a power system node carbon emission factor correction method and device, which can solve the problem of low accuracy of carbon emission factor calculation in the prior art.

[0005] In order to solve the above technical problems, the present application provides a power system node carbon emission factor correction method, comprising:

[0006] Real-time collection of SCADA data of the power system; wherein the SCADA data includes node voltage amplitude, node voltage phase, node active power, node reactive power, branch active power and branch reactive power;

[0007] Real-time collection of carbon emission data of the generator node of the power system;

[0008] Obtaining the power grid node topology data of the power system; wherein the power grid node topology data includes the number of node voltages, the number of node injected powers and the number of branch powers;

[0009] Based on the SCADA data and the power grid node topology data, a nonlinear measurement model is constructed;

[0010] solving the nonlinear measurement model to obtain voltage correction values of all nodes of the power system; wherein the voltage correction values include node voltage amplitude correction values and node voltage phase correction values;

[0011] calculating complex power of each branch in the power system based on the voltage correction values;

[0012] decomposing the complex power to obtain active power correction values of each branch;

[0013] correcting carbon emission factors of each node in the power system based on the active power correction values, the carbon emission data and the grid node topology data.

[0014] As a preferred solution, the constructing a nonlinear measurement model based on the SCADA data and the grid node topology data comprises:

[0015] determining node topology relationships and branch types of the power system according to the grid node topology data;

[0016] calculating branch admittance data based on the branch types and element parameters of the power system;

[0017] respectively constructing node power balance equations and branch power equations according to the SCADA data, the node topology relationships and the branch admittance data;

[0018] constructing state equations of the nonlinear measurement model based on the node power balance equations and the branch power equations;

[0019] adopting a weighted least square method to construct an objective function of the nonlinear measurement model with the objective of minimizing residual sum of squares.

[0020] As a preferred solution, the node power balance equation is specifically:

[0021]

[0022] In the formula, P i is node active power of node i; Q i is node reactive power of node i; V i is node voltage amplitude of node i; V j is node voltage amplitude of node j; θ ij = θ i - θ j , θ i is node voltage phase of node i, θ j is node voltage phase of node j; G ij is electric conductance of branch ij; B ij is electric susceptance of branch ij; n = 2Nv +2N p +2B S , N v is the number of node voltages, N p is the number of node injection powers, B S is the number of branch powers.

[0023] As a preferred solution, the branch power equation is specifically:

[0024]

[0025] In the formula, P ij is the branch active power of branch ij; Q ij is the branch reactive power of branch ij; V i is the node voltage amplitude of node i; V j is the node voltage amplitude of node j; θ ij = θ i - θ j , θ i is the node voltage phase of node i, θ j is the node voltage phase of node j; G ij is the conductance of branch ij; B ij is the susceptance of branch ij; g i0 is the ground conductance of the branch; and b i0 is the ground susceptance of the branch.

[0026] As a preferred solution, the state equation is specifically:

[0027] z = h(x) + w

[0028] In the formula, z is the state equation; h(x) is a nonlinear function, including the node power balance equation and the branch power equation; and w is a preset variable.

[0029] As a preferred solution, the objective function is specifically:

[0030] J(x) = [z - h(x)] T W[z - h(x)]

[0031] In the formula, J(x) is the objective function; h(x) is a nonlinear function; z is the state equation; and W is a weight matrix.

[0032] As a preferred solution, the calculation of the complex power of each branch in the power system based on the voltage correction value includes:

[0033] The complex power of each branch in the power system is calculated by the following formula:

[0034]

[0035] Sij=V ij is the complex power of branch ij; Y ij is the branch admittance, Y ij =G ij +jB ij ; V i 'is the node voltage magnitude correction value of node i; V j 'is the node voltage magnitude correction value of node j; θ ij =θ i -θ j , θ i is the node voltage phase correction value of node i, θ j is the node voltage phase correction value of node j.

[0036] As a preferred solution, the decomposing the complex power to obtain the active power correction quantity of each branch comprises:

[0037] The active power correction quantity of each branch is calculated by using the following formula:

[0038]

[0039] P ij 'is the active power correction quantity of branch ij; Q ij 'is the reactive power correction quantity of branch ij; V i 'is the node voltage magnitude correction value of node i; V j 'is the node voltage magnitude correction value of node j; θ ij =θ i -θ j , θ i is the node voltage phase correction value of node i, θ j is the node voltage phase correction value of node j; G ij is the conductance of branch ij; B ij is the susceptance of branch ij.

[0040] As a preferred solution, the correcting the carbon emission factor of each node in the power system based on the active power correction quantity, the carbon emission data and the power grid node topology data comprises:

[0041] For each node, a plurality of power input nodes of the node are determined according to the node topology relationship;

[0042] The carbon emission factor before correction of each power input node is determined according to the carbon emission data;

[0043] The carbon emission factor of the node is calculated based on the carbon emission factor before correction of each power input node by using the following formula:

[0044]

[0045] wherein e i is the carbon emission factor of node i; P ij is the active power correction of branch ij; e j is the carbon emission factor of node j before correction.

[0046] Correspondingly, the application provides a device for correcting carbon emission factors of nodes in a power system, comprising a SCADA data acquisition module, a carbon emission data acquisition module, a topology data acquisition module, a model construction module, a model solution module, a complex power calculation module, a power decomposition module and a correction module.

[0047] The SCADA data acquisition module is configured to acquire SCADA data of the power system in real time; wherein the SCADA data comprises node voltage amplitude, node voltage phase, node active power, node reactive power, branch active power and branch reactive power.

[0048] The carbon emission data acquisition module is configured to acquire carbon emission data of generator nodes of the power system in real time.

[0049] The topology data acquisition module is configured to acquire grid node topology data of the power system; wherein the grid node topology data comprises node voltage quantity, node injected power quantity and branch power quantity.

[0050] The model construction module is configured to construct a nonlinear measurement model based on the SCADA data and the grid node topology data.

[0051] The model solution module is configured to solve the nonlinear measurement model to obtain voltage correction values of all nodes of the power system; wherein the voltage correction values comprise node voltage amplitude correction values and node voltage phase correction values.

[0052] The complex power calculation module is configured to calculate complex power of each branch in the power system based on the voltage correction values.

[0053] The power decomposition module is configured to decompose the complex power to obtain active power correction of each branch.

[0054] The correction module is configured to correct carbon emission factors of each node in the power system based on the active power correction, the carbon emission data and the grid node topology data.

[0055] Compared with the prior art, the embodiments of the application have the following beneficial effects:

[0056] The application provides a power system node carbon emission factor correction method, which comprises the following steps: collecting SCADA data of a power system in real time; collecting carbon emission data of a generator node of the power system in real time; obtaining power grid node topology data of the power system; constructing a nonlinear measurement model based on the SCADA data and the power grid node topology data; solving the nonlinear measurement model to obtain voltage correction values of all nodes of the power system; calculating complex power of each branch in the power system based on the voltage correction values; decomposing the complex power to obtain active power correction values of each branch; and correcting carbon emission factors of each node in the power system based on the active power correction values, the carbon emission data and the power grid node topology data. BRIEF DESCRIPTION OF DRAWINGS

[0057] In order to more clearly illustrate the technical solutions of the present application, the following will briefly introduce the drawings needed to be used in the embodiments. Obviously, the drawings described in the following are only some of the embodiments of the present application, and other drawings can also be obtained by those skilled in the art without any creative effort based on these drawings.

[0058] Figure 1 The flowchart of an embodiment of the power system node carbon emission factor correction method provided by the present application is shown in the figure.

[0059] Figure 2 The structural schematic diagram of an embodiment of the power system provided by the present application is shown in the figure.

[0060] Figure 3 The structural schematic diagram of an embodiment of the power system node carbon emission factor correction device provided by the present application is shown in the figure. DETAILED DESCRIPTION

[0061] In order to make the purpose, technical solutions and advantages of the present application more clear, the technical solutions in the present application will be described clearly and completely in the following with reference to the drawings in the embodiments of the present application. Obviously, the described embodiments are only some of the embodiments of the present application, but not all the embodiments. Based on the embodiments in the present application, all other embodiments obtained by those skilled in the art without any creative effort fall within the scope of protection of the present application.

[0062] Unless otherwise defined, all technical and scientific terms used herein have the same meaning as commonly understood by one of ordinary skill in the art to which this application belongs; the terminology used herein is for the purpose of describing particular embodiments only and is not intended to be limiting of this application; the use of the terms "including," "comprising," or "having" and variations thereof herein is intended to be broad and encompass the terms "consisting of" and "consisting essentially of" and variations thereof. Unless otherwise required by context, singular terms shall include pluralities and vice versa. Unless otherwise required by context, the use herein of the singular is also to be construed as a use of the plural and vice versa.

[0063] In the description of the embodiments of the present application, the technical terms "first", "second", etc. are only used to distinguish different objects, and cannot be understood as indicating or implying relative importance or implicitly indicating the number, specific order or primary and secondary relationship of the indicated technical features. In the description of the embodiments of the present application, the meaning of "multiple" is more than two, unless otherwise explicitly specified.

[0064] Reference herein to "an embodiment" means that a particular feature, structure, or characteristic described in connection with the embodiment can be included in at least one embodiment of the application. The appearance of the phrase in various places in the specification does not necessarily all refer to the same embodiment, nor is it necessarily independent or alternative embodiments to each other. It is explicitly and implicitly understood by those skilled in the art that the embodiments described herein can be combined with other embodiments.

[0065] In the description of the embodiments of the present application, the term "and / or" is only a description of the association relationship of the associated objects, which means that there can be three relationships, for example, A and / or B, which can represent the three cases of A alone, A and B together, and B alone. In addition, the character " / " in this paper generally represents that the front and rear associated objects are a "or" relationship.

[0066] In the description of the embodiments of the present application, the term "multiple" refers to more than two (including two), and similarly, "multiple groups" refers to more than two groups (including two groups), and "multiple pieces" refers to more than two pieces (including two pieces).

[0067] In the description of the embodiments of the present application, unless otherwise explicitly specified and limited, the technical terms "mounting", "connecting", "connecting", "fixing" and the like should be broadly understood, for example, it can be fixedly connected, or it can be detachably connected, or it can be integrated; it can be mechanical connection, or it can be electrical connection; it can be directly connected, or it can be indirectly connected through an intermediate medium; it can be the internal communication of two elements or the interaction relationship between two elements. For those skilled in the art, the specific meaning of the above terms in the embodiments of the present application can be understood according to the specific circumstances.

[0068] Reference Figure 1To solve the problem of low accuracy of carbon emission factor calculation in the prior art, an embodiment of the present application provides a power system node carbon emission factor correction method, which comprises steps 101 to 108, and each step is as follows:

[0069] Step 101: Real-time collection of SCADA data of the power system; wherein the SCADA data comprises node voltage amplitude, node voltage phase, node active power, node reactive power, branch active power and branch reactive power.

[0070] In the embodiment of the present application, SCADA (Supervisory Control And Data Acquisition) is a data acquisition and monitoring control system for industrial field process control and data acquisition. The SCADA data refers to various types of data related to industrial production processes collected, processed and stored by the system during operation, which is of great significance for monitoring production processes, optimizing operation and ensuring safety. The SCADA data collected in the present application includes node voltage amplitude, node voltage phase, node active power, node reactive power, branch active power and branch reactive power. After obtaining the SCADA data, abnormality detection is performed on the data, and non-abnormal data is extracted from the collected SCADA data as the basis for subsequent carbon emission factor correction.

[0071] Step 102: Real-time collection of carbon emission data of the generator node of the power system.

[0072] In the embodiment of the present application, the carbon emission data of the generator node is the carbon emission factor calculated based on the power flow data of the power system at the current time, which is the carbon emission data to be corrected.

[0073] Step 103: Obtain grid node topology data of the power system; wherein the grid node topology data comprises node voltage quantity, node injection power quantity and branch power quantity.

[0074] In the embodiment of the present application, the grid node topology data of the power system contains the node topology relationship and branch type of the power system, as well as the node quantity and branch quantity of the grid topology. According to the grid node topology data, the node voltage quantity, node injection power quantity and branch power quantity can be determined. Among them, the node voltage quantity and the node injection power quantity are equal to the node quantity in the grid topology, and the branch power quantity is equal to the branch quantity in the grid topology.

[0075] Step 104: Construction of a nonlinear measurement model based on the SCADA data and the grid node topology data.

[0076] In the embodiment of the present application, the SCADA data of the power system does not necessarily include data of all nodes on the power system, and by constructing the nonlinear measurement model, all data of the nodes on the power system can be corrected based on the SCADA data containing partial node data.

[0077] As a preferred scheme of the embodiment, the nonlinear measurement model is constructed based on the SCADA data and the grid node topology data, and the nonlinear measurement model includes:

[0078] The node topology relationship and the branch type of the power system are determined according to the grid node topology data;

[0079] The branch admittance data are calculated based on the branch type and the element parameters of the power system;

[0080] The node power balance equation and the branch power equation are respectively constructed according to the SCADA data, the node topology relationship and the branch admittance data;

[0081] The state equation of the nonlinear measurement model is constructed based on the node power balance equation and the branch power equation;

[0082] The objective function of the nonlinear measurement model is constructed by using the weighted least square method with the objective of minimizing the residual sum of squares.

[0083] In the embodiment of the present application, the known input of the nonlinear measurement model is the SCADA data and the grid node topology data, and the unknown output is the voltage correction value of all nodes of the power system. The state equation of the nonlinear measurement model includes two types of equations, i.e., the node power balance equation and the branch power equation. The objective of the nonlinear measurement model is to minimize the residual sum of squares, and the objective function can be constructed by using the weighted least square method.

[0084] As a preferred scheme of the embodiment, the state equation is specifically:

[0085] z = h(x) + w

[0086] In the formula, z is the state equation, h(x) is a nonlinear function including the node power balance equation and the branch power equation, w is a preset variable, and x represents the voltage correction value of the node.

[0087] In the embodiment of the present application, the nonlinear function h(x) includes the node power balance equation and the branch power equation.

[0088] As a preferred scheme of the embodiment, the node power balance equation is specifically:

[0089]

[0090] In the formula, P iPi is the node active power of node i; Q i Qi is the node reactive power of node i; V i Vi is the node voltage amplitude of node i; V j Vj is the node voltage amplitude of node j; θ ij = θ i - θ j , θ i is the node voltage phase of node i, θ j is the node voltage phase of node j; G ij is the conductance of branch ij; B ij is the susceptance of branch ij; n is the dimension, n = 2N v + 2N p + 2B S , N v is the number of node voltages, N p is the number of node injection powers, B S is the number of branch powers.

[0091] In the embodiment of the present application, in the node power balance equation constructed above, the node voltage amplitude and the node voltage phase are obtained from the SCADA data, the number of node voltages, the number of node injection powers and the number of branch powers are obtained according to the node topological relationship in the power grid node topological data, and the branch conductance and the branch susceptance are calculated according to the branch type and the element parameters of the power system.

[0092] As a preferred scheme of the embodiment, the branch power equation is specifically:

[0093]

[0094] In the formula, P ij is the branch active power of branch ij; Q ij is the branch reactive power of branch ij; V i is the node voltage amplitude of node i; V j is the node voltage amplitude of node j; θ ij = θ i - θ j , θ i is the node voltage phase of node i, θ j is the node voltage phase of node j; G ij is the conductance of branch ij; B ij is the susceptance of branch ij; g i0 is the ground conductance of branch; and b i0 is the ground susceptance of branch.

[0095] In the embodiment of the present application, the branch-to-ground conductance and the branch-to-ground susceptance can be calculated according to the element parameters of the power system, for example, the device nodes such as transformers can be calculated by equivalent circuits and no-load test parameters (such as no-load loss and excitation reactive power). The branch-to-ground conductance and the branch-to-ground susceptance can also be calculated by using a power frequency parameter tester (such as an off-frequency method) to apply a power frequency voltage to the line, measure the current and phase difference, and calculate the susceptance.

[0096] As a preferred scheme of the embodiment, the objective function, specifically,

[0097] J(x) = [z - h(x)] T W[z - h(x)]

[0098] In the formula, J(x) is the objective function; h(x) is a nonlinear function; z is a state equation; and W is a weight matrix.

[0099] In the embodiment of the present application, the weight matrix W is a diagonal matrix, and the elements are:

[0100]

[0101] After the state equation and the objective function of the nonlinear measurement model are determined, the Newton-Raphson iteration method can be used to solve the optimal correction value, and the voltage correction value of all nodes of the power system is obtained.

[0102] Step 105: solving the nonlinear measurement model to obtain the voltage correction value of all nodes of the power system; wherein the voltage correction value includes a node voltage amplitude correction value and a node voltage phase correction value.

[0103] In the embodiment of the present application, the nonlinear measurement model can be solved by the Newton-Raphson iteration method, and the initial voltage correction value is first set as x (0) , and the iteration formula is calculated, which is expressed as follows:

[0104] Δx (k) = (H T WH) -1 H T W(Z - h(x (k) ))

[0105] x (k+1) = x (k) + Δx (k)

[0106] In the formula, H is a Jacobian matrix: x (k) is the voltage correction value at the kth iteration; Δx (k) is the kth voltage correction difference; and x (k+1) is the voltage correction value at the (k+1)th iteration.

[0107] In the solving process, whether the result converges is judged by residual norm and maximum iteration number. The residual norm is expressed as:

[0108] ||E-h(x (k) )||<ε

[0109] In the formula, E represents an expected value, ε represents a preset convergence threshold, and ||·|| represents a norm.

[0110] Step 106: calculating complex power of each branch in the power system based on the voltage correction value.

[0111] As a preferred scheme of the embodiment, the complex power of each branch in the power system is calculated based on the voltage correction value, and the calculation includes:

[0112] The complex power of each branch in the power system is calculated by using the following formula:

[0113]

[0114] In the formula, S ij represents the complex power of the branch ij; Y ij represents branch admittance, Y ij =G ij +jB ij ; V i ' represents the node voltage amplitude correction value of node i; V j ' represents the node voltage amplitude correction value of node j; θ ij =θ i -θ j , θ i represents the node voltage phase correction value of node i, and θ j represents the node voltage phase correction value of node j.

[0115] In the embodiment, when the voltage correction values of all nodes in the power system are obtained by solving the nonlinear measurement model, the complex power of the branch is calculated by using the following formula:

[0116]

[0117] In the formula, I ij represents branch current, which is calculated by using branch admittance and node voltage, and the calculation formula is as follows:

[0118] Step 107: decomposing the complex power to obtain the active power correction amount of each branch.

[0119] As a preferred scheme of the embodiment, the complex power is decomposed to obtain the active power correction amount of each branch, and the decomposition includes:

[0120] The active power correction amount of each branch is calculated by the following formula:

[0121]

[0122] In the formula, P ij is the active power correction amount of branch ij; Q ij is the reactive power correction amount of branch ij; V i is the node voltage amplitude correction value of node i; V j is the node voltage amplitude correction value of node j; θ ij = θ i - θ j , θ i is the node voltage phase correction value of node i; θ j is the node voltage phase correction value of node j; G ij is the conductance of branch ij; and B ij is the susceptance of branch ij.

[0123] In the embodiment of the present application, after the complex power of each branch in the power system is calculated according to the voltage correction values of all nodes in the power system, the active power correction amount of each branch can be obtained by analyzing the complex power.

[0124] Step 108: Correct the carbon emission factor of each node in the power system based on the active power correction amount, the carbon emission data and the grid node topology data.

[0125] As a preferred scheme of the embodiment, the carbon emission factor of each node in the power system is corrected based on the active power correction amount, the carbon emission data and the grid node topology data, and the method comprises the following steps.

[0126] For each node, a plurality of power input nodes of the node are determined according to the node topology relationship;

[0127] The carbon emission factor of each power input node before correction is determined according to the carbon emission data;

[0128] The carbon emission factor of the node is calculated based on the carbon emission factor of each power input node before correction by the following formula:

[0129]

[0130] In the formula, e i is the carbon emission factor of node i; P ij is the active power correction amount of branch ij; e j is the carbon emission factor of node j before correction.

[0131] In an embodiment of the present invention, in the power system, the active power flow carries carbon emissions starting from the power source and flows to the load node via the intermediate node, thereby forming a "carbon flow". The carbon flows from different power sources are mixed at each node, resulting in different "carbon content" of the injected power at different nodes. The power flow path of each node can be determined through the node topology relationship of the power system, so the power input node of each node can be obtained. The power input node is the upstream node of its corresponding node and is used to input power to the node. After determining the active power correction amount of each branch of the power system and the power input node of each node, the carbon emission factor of each node can be corrected according to the above formula.

[0132] In an embodiment of the present invention, based on the carbon emission factor correction method for each node described above, a carbon emission matrix of the power system can be obtained:

[0133]

[0134] Where, E N is the carbon emission matrix; P N is the active flux matrix; P B is the branch power flow distribution matrix; P G Inject the distribution matrix into the generator set; E G is the carbon emission intensity distribution matrix of the generator set.

[0135] In an embodiment of the present invention, after obtaining the corrected carbon emission factor for each node in the power system, the carbon emission factor of the generator set node is compared with a preset carbon emission threshold. When the sum of the carbon emission factors of the generator set node is greater than the preset carbon emission threshold, adjustments are made to each generator set. For example, for coal-fired power units, the carbon emission factor can be reduced by adjusting the boiler combustion temperature and the steam turbine vacuum to reduce coal consumption for power generation, and by controlling the oxygen content of the flue gas at the boiler outlet. For gas-fired units, the natural gas consumption per unit of electricity can be reduced by adjusting the fuel nozzle angle and the compressor inlet guide vane opening, thereby improving the thermal efficiency of the gas turbine; and the number of unit starts and stops can be reduced, thereby reducing the carbon emission factor by extending the continuous operation time.

[0136] In an embodiment of the present invention, after the carbon emission factor of each node in the power system is corrected, the uncertainty of the carbon emission factor of each node can be calculated using the Monte Carlo method. Specifically:

[0137] A mathematical model for calculating the uncertainty of carbon emission factors (output) of a power system is established; the input quantities that have an impact on the calculation results of carbon emission intensity are determined according to the established output mathematical model, and the distribution and probability density function of the input quantities are determined through a large number of tests; the number of Monte Carlo method (MCM method) experiments M is determined; M times of random sampling are performed in the input quantities to obtain corresponding sample values; the sample values are brought into the mathematical model to obtain the model values of the output quantities; the M model values of the output quantities are incrementally sorted to obtain a discrete representation of the distribution function of the output quantities; the average value and the standard uncertainty of the output quantities are determined, and the given containing probability and the containing interval of the measured quantities are determined from the discrete representation of the distribution function of the output quantities.

[0138] As an example of an embodiment of the present application, refer to Figure 2 , which is a structural schematic diagram of an embodiment of the power system provided by the present application, which is an EEE 3 machine 9 node system structure, containing three types of machine groups G1, G2 and G3, G1 connects node 2, G2 connects node 4, and G3 connects node 3; nodes are connected through lines, node 2 has power 8 flowing out, node 3 has power 3 flowing in, node 4 has power 1 flowing in, node 5 has power flowing out to load 1, node 7 has power flowing out to load 2, node 9 has power flowing out to load 3, each node and line together constitute a power transmission network, realizing power transmission from machine groups to loads. The parameter settings of each generator group are shown in Table 1 below:

[0139] Generator node number Carbon emission factor 1 1.06 kg / kWh 2 0.5 kg / kWh 3 0 kg / kWh

[0140] Table 1: Parameter settings of generator groups

[0141] SCADA is used to collect node data and branch data of the power system, wherein nodes 1, 2, 6 and 9 are data collection nodes, and the collected data are used as observations, and the measured node data and branch data at a certain time are shown in Tables 2 and 3 below:

[0142]

[0143]

[0144] Table 2: Node data of some nodes

[0145] Node number Active injection (MV) Reactive injection (MV) 4 5 0.3070 0.103 3 6 0.85 -0.1086 6 7 0.2418 0.0312 8 9 0.8662 -0.0838

[0146] Table 3: Branch active power flow

[0147] Based on the above collected part of the node data and branch data, a nonlinear measurement model is constructed, and the Newton-Raphson iteration method is used to solve the nonlinear measurement model to obtain the voltage correction value of all nodes in the power system, and then the carbon emission factors of each node in the power system are corrected.

[0148] The above embodiment is implemented, and has the following effects:

[0149] The application provides a power system node carbon emission factor correction method, real-time SCADA data of a power system is collected, carbon emission data of a generator node of the power system is collected in real time, power grid node topology data of the power system is acquired, a nonlinear measurement model is constructed based on the SCADA data and the power grid node topology data, the nonlinear measurement model is solved, voltage correction values of all nodes of the power system are obtained, complex power of each branch in the power system is calculated based on the voltage correction values, active power correction values of each branch are obtained by decomposing the complex power, and carbon emission factors of each node in the power system are corrected based on the active power correction values, the carbon emission data and the power grid node topology data. The SCADA data and the power grid node topology data of the generator set are collected, the voltage correction values of all nodes of the power system are solved by using the weighted least square method, the power flow distribution of the entire power system can be measured, the carbon emission data of the generator node are corrected according to the voltage correction values, the corrected carbon emission factors of each node in the power system are obtained, and the accuracy of the carbon emission factors can be effectively improved.

[0150] As Figure 3 shown, on the basis of the above method embodiment, a corresponding device embodiment is provided;

[0151] An embodiment of the application provides a power system node carbon emission factor correction device, which comprises a SCADA data acquisition module, a carbon emission data acquisition module, a topology data acquisition module, a model construction module, a model solving module, a complex power calculation module, a power decomposition module and a correction module.

[0152] The SCADA data acquisition module is used for collecting SCADA data of a power system in real time, wherein the SCADA data comprises node voltage amplitude, node voltage phase, node active power, node reactive power, branch active power and branch reactive power.

[0153] The carbon emission data acquisition module is used for collecting carbon emission data of a generator node of the power system in real time.

[0154] The topology data acquisition module is used for acquiring power grid node topology data of the power system, wherein the power grid node topology data comprises node voltage quantity, node injected power quantity and branch power quantity.

[0155] The model construction module is used for constructing a nonlinear measurement model based on the SCADA data and the power grid node topology data.

[0156] The model solving module is configured to solve the nonlinear measurement model to obtain voltage correction values of all nodes of the power system; wherein the voltage correction values include node voltage amplitude correction values and node voltage phase correction values;

[0157] The complex power calculation module is configured to calculate complex power of each branch in the power system based on the voltage correction values;

[0158] The power decomposition module is configured to decompose the complex power to obtain active power correction values of each branch;

[0159] The correction module is configured to correct carbon emission factors of each node in the power system based on the active power correction values, the carbon emission data and the grid node topology data.

[0160] It can be understood that the above-mentioned device item embodiments are corresponding to the method item embodiments of the present application, and can realize the power system node carbon emission factor correction method provided by any one of the above-mentioned method item embodiments of the present application.

[0161] It should be noted that the device embodiments described above are only schematic, and part or all of the modules thereof can be selected to achieve the purpose of the present embodiment scheme. In addition, in the device embodiment provided by the present application, the connection relationship between the modules indicates that there is a communication connection between them, which can be realized as one or more communication buses or signal lines. Those skilled in the art can understand and implement it without creative labor.

[0162] The above-described specific embodiments further illustrate the purpose, technical solutions and beneficial effects of the present application. It should be understood that the above-described specific embodiments are only for the purpose of the present application and do not limit the protection scope of the present application. It is particularly pointed out that any modification, equivalent replacement, improvement, etc. made by those skilled in the art within the spirit and principles of the present application shall be included in the protection scope of the present application.

Claims

1. A method for correcting carbon emission factors for nodes of an electric power system, characterized by, The method comprises: collecting SCADA data of a power system in real time; wherein the SCADA data comprises node voltage amplitude, node voltage phase, node active power, node reactive power, branch active power and branch reactive power; collecting carbon emission data of a generator node of the power system in real time; obtaining power grid node topology data of the power system; wherein the power grid node topology data comprises node voltage quantity, node injected power quantity and branch power quantity; constructing a nonlinear measurement model based on the SCADA data and the power grid node topology data; solving the nonlinear measurement model to obtain voltage correction values of all nodes of the power system; wherein the voltage correction values comprise node voltage amplitude correction values and node voltage phase correction values; calculating complex power of each branch in the power system based on the voltage correction values; decomposing the complex power to obtain active power correction values of each branch; correcting carbon emission factors of each node in the power system based on the active power correction values, the carbon emission data and the power grid node topology data.

2. The method of claim 1, wherein, The method of constructing the nonlinear measurement model based on the SCADA data and the power grid node topology data comprises: determining node topology relationships and branch types of the power system according to the power grid node topology data; calculating branch admittance data based on the branch types and element parameters of the power system; respectively constructing node power balance equations and branch power equations according to the SCADA data, the node topology relationships and the branch admittance data; constructing a state equation of the nonlinear measurement model based on the node power balance equations and the branch power equations; constructing an objective function of the nonlinear measurement model by using a weighted least squares method with the objective of minimizing residual sum of squares.

3. The power system node carbon emission factor revision method of claim 2, wherein, The node power balance equation is specifically: where P i is the node active power of node i; Q i is the node reactive power of node i; V i is the node voltage magnitude of node i; V j is the node voltage magnitude of node j; θ ij = θ i - θ j , θ i is the phase of the node voltage at node i, θ j is the phase of the node voltage at node j; G ij is the conductance of branch ij; B ij is the susceptance of branch ij; n = 2N v + 2N p + 2B S , N v is the number of node voltages, N p is the number of node injections, B S is the number of branch powers.

4. The power system node carbon emission factor revision method of claim 3, wherein, The branch power equation is specifically: Where, P ij is the active power of branch ij; Q ij is the branch reactive power of branch ij; V i is the node voltage amplitude of node i; V j is the node voltage amplitude of node j; θ ij = θ i - θ j , θ i is the phase of the node voltage at node i, θ j is the phase of the node voltage at node j; G ij is the conductance of branch ij; B ij is the susceptance of branch ij; g i0 is the conductance to ground of branch; and b i0 is the susceptance to ground of branch.

5. The power system node carbon emission factor revision method of claim 4, wherein, The state equation is specifically: z=h(x)+w In the formula, z is a state equation; h(x) is a nonlinear function, which comprises node power balance equations and branch power equations; and w is a preset variable.

6. The power system node carbon emission factor revision method of claim 5, wherein, The objective function is specifically: J(x) = [z - h(x)] T W[z - h(x)] In the formula, J(x) is an objective function; h(x) is a nonlinear function; z is a state equation; and W is a weight matrix.

7. The power system node carbon emission factor revision method of claim 6, wherein, The method of calculating complex power of each branch in the power system based on the voltage correction values comprises: The method of decomposing the complex power to obtain active power correction values of each branch comprises: where S ij is the complex power of branch ij; Y ij is the branch admittance, Y ij = G ij + jB ij ; V' i is the node voltage magnitude correction value of node i; V j ' is the node voltage magnitude correction value of node j; θ' ij = θ' i - θ' j , θ' i is the node voltage phase correction value of node i; and θ' j is the node voltage phase correction value of node j.

8. The power system node carbon emission factor revision method of claim 7, wherein, The method of calculating active power correction values of each branch comprises: The method of correcting carbon emission factors of each node in the power system based on the active power correction values, the carbon emission data and the power grid node topology data comprises: Where P′ ij is the active power correction value of branch ij; Q′ ij is the reactive power correction value of branch ij; V′ i is the node voltage amplitude correction value of node i; V′ j is the node voltage amplitude correction value of node j; θ' ij =θ' i -θ' j ,θ' i is the node voltage phase correction value of node i, θ' j is the node voltage phase correction value of node j; G ij is the conductance of branch ij; B ij is the susceptance of branch ij.

9. The power system node carbon emission factor revision method of claim 8, wherein, For each node, a plurality of power input nodes of the node are determined according to the node topology relationships; correction-before carbon emission factors of each power input node are determined according to the carbon emission data; carbon emission factors of the node are calculated based on the correction-before carbon emission factors of each power input node by using the following formula: The method comprises: wherein e i is the carbon emission factor of node i; P′ ij is the active power correction of branch ij; e j is the carbon emission factor of node j before correction.

10. An apparatus for correcting carbon emission factors for power system nodes, characterized by, ​ The SCADA data acquisition module, the carbon emission data acquisition module, the topology data acquisition module, the model construction module, the model solving module, the complex power calculation module, the power decomposition module and the correction module; The SCADA data acquisition module is configured to acquire SCADA data of the power system in real time; wherein the SCADA data comprises node voltage amplitude, node voltage phase, node active power, node reactive power, branch active power and branch reactive power; The carbon emission data acquisition module is configured to acquire carbon emission data of a generator node of the power system in real time; The topology data acquisition module is configured to acquire power grid node topology data of the power system; wherein the power grid node topology data comprises node voltage quantity, node injection power quantity and branch power quantity; The model construction module is configured to construct a nonlinear measurement model based on the SCADA data and the power grid node topology data; The model solving module is configured to solve the nonlinear measurement model to obtain voltage correction values of all nodes of the power system; wherein the voltage correction values comprise node voltage amplitude correction values and node voltage phase correction values; The complex power calculation module is configured to calculate complex power of each branch in the power system based on the voltage correction values; The power decomposition module is configured to decompose the complex power to obtain active power correction values of each branch; The correction module is configured to correct carbon emission factors of each node in the power system based on the active power correction values, the carbon emission data and the power grid node topology data.

Images