CN117764798B - Method and system for checking carbon meter measurement data of user - Google Patents
Method and system for checking carbon meter measurement data of user Download PDFInfo
- Publication number
- CN117764798B CN117764798B CN202410194801.XA CN202410194801A CN117764798B CN 117764798 B CN117764798 B CN 117764798B CN 202410194801 A CN202410194801 A CN 202410194801A CN 117764798 B CN117764798 B CN 117764798B
- Authority
- CN
- China
- Prior art keywords
- carbon emission
- user
- carbon
- tolerance
- power grid
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
- OKTJSMMVPCPJKN-UHFFFAOYSA-N Carbon Chemical compound [C] OKTJSMMVPCPJKN-UHFFFAOYSA-N 0.000 title claims abstract description 185
- 229910052799 carbon Inorganic materials 0.000 title claims abstract description 185
- 238000005259 measurement Methods 0.000 title claims abstract description 84
- 238000000034 method Methods 0.000 title claims abstract description 39
- 238000004364 calculation method Methods 0.000 claims description 23
- 239000011159 matrix material Substances 0.000 claims description 18
- 230000002159 abnormal effect Effects 0.000 claims description 13
- 238000012549 training Methods 0.000 claims description 8
- 238000010276 construction Methods 0.000 claims description 6
- 230000008030 elimination Effects 0.000 claims description 6
- 238000003379 elimination reaction Methods 0.000 claims description 6
- 238000012937 correction Methods 0.000 claims description 5
- 238000012706 support-vector machine Methods 0.000 claims description 5
- 238000004891 communication Methods 0.000 abstract description 6
- 239000003245 coal Substances 0.000 description 4
- 230000005611 electricity Effects 0.000 description 4
- 238000010586 diagram Methods 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 238000005303 weighing Methods 0.000 description 2
- 238000002485 combustion reaction Methods 0.000 description 1
- 238000009795 derivation Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000005265 energy consumption Methods 0.000 description 1
- 239000000843 powder Substances 0.000 description 1
- 238000010248 power generation Methods 0.000 description 1
- 230000035945 sensitivity Effects 0.000 description 1
Landscapes
- Investigating, Analyzing Materials By Fluorescence Or Luminescence (AREA)
Abstract
The invention provides a method and a system for checking measurement data of a user carbon table, wherein the method comprises the following steps: acquiring carbon emission measurement values of a current updating period of all user carbon tables of a current regional power grid; constructing a set of carbon emission measurement values of each group; correspondingly forming an equation set aiming at each group of carbon emission measurement value sets, and solving the intermediate quantity of each group through each equation set; and judging whether the intermediate quantity of each group is out of tolerance, if so, correcting a carbon emission metering value set corresponding to the out of tolerance intermediate quantity, thereby realizing the purpose of checking the metering result of the carbon table of the user, identifying the out of tolerance result in the carbon emission quantity of the carbon table of the user caused by network, communication, time delay and the like, correcting the out of tolerance result, and making up for the blank of the lack of a checking method for the carbon emission quantity of the carbon table of the current user.
Description
Technical Field
The invention relates to the technical field of carbon emission measurement, in particular to a method and a system for checking carbon meter measurement data of a user.
Background
Accurate measurement of carbon emission is an important basis for grasping emission contribution and clear emission reduction responsibility of each participant in an energy system. Along with the development of a carbon flow theory, the carbon emission factors of all nodes of the whole network can be calculated by obtaining the carbon emission intensity of a power supply, so that the carbon emission factors can be combined with an electric energy meter of a user to form a carbon meter, and the carbon emission amount generated by the energy consumption of the user is measured.
The carbon emission measurement of the user carbon table needs to obtain the carbon emission factor of the current access node, and the factor can be an annual average published by related organizations, or a dynamic area average issued and updated in real time by a platform area or a dispatching center. However, due to the influence of various aspects such as network, communication and time delay, the carbon emission factor value adopted when the carbon table updates the carbon emission result may have the conditions of coding error, bad data, untimely updating and the like, thereby further causing the output carbon emission result error.
At present, the calculation and the check of the carbon emission are mainly aimed at determining the emission by a coal and electricity enterprise in a flowmeter or weighing mode, and the check is mainly used for judging whether the emission of the coal and electricity enterprise is in an allowable range or not, and the lack of a method for checking a carbon meter measurement result based on a carbon flow theory may lead to inaccurate carbon emission calculation result of a user.
The calculation and the acquisition of the carbon emission are mainly divided into two types at present, one type is to calculate the emission of coal-electricity enterprises, and the two types comprise weighing coal powder and detecting and calculating the emission by a flowmeter for combustion emission; one is to distribute the carbon emission of the power generation part to the user by the carbon flow theory. The carbon meter calculates the carbon emission of the user by using the node carbon emission factor and the user electricity consumption, so that the emission contribution of each user can be mastered, and the emission reduction responsibility can be definitely realized.
In the prior art, the calculation and the check of the carbon emission are mainly aimed at the calculation of the emission of coal and electricity enterprises and the check of whether the emission exceeds the standard, and a check method for the measurement result of a carbon meter of a user is lacked.
Disclosure of Invention
The invention aims to solve the technical problem of providing a method and a system for checking the metering data of a user carbon meter, wherein the method and the system calculate an intermediate value through the carbon emission of each user and check whether the carbon emission exceeds the standard, and calculate a new carbon emission to correct the original result if the carbon emission exceeds the standard so as to check the metering result of the user carbon meter.
In a first aspect, the present invention provides a method for checking measurement data of a user carbon table, including the following steps:
S101, acquiring carbon emission measurement values C i, i=1, 2, … and n of all user carbon tables of a current regional power grid in n current update periods;
S102, constructing a set Z of n-3m+1 groups of carbon emission measurement values p, , P=1, 2, …, n-3m+1, m is the number of power supplies generating carbon emission in the current regional power grid, m < < n, i.e. the 1 st to 3m carbon emission measurement values are selected as a1 st group carbon emission measurement value set Z 1, and the 2 nd to 3m+1, 3 rd to 3m+2, …, n-3m+1-n carbon emission measurement values are sequentially selected as a 2 nd, 3 rd, … … th, n-3m+1 th group carbon emission measurement value set Z 2、Z3、……、Zn-3m+1;
S103, correspondingly forming an equation set theta p, p=1, 2, … and n-3m+1 for each group of carbon emission measurement value sets Z p;
；
Wherein:
C p+j-1 is the carbon emission measurement value of the p+j-1 th user;
t is the carbon emission factor update period;
K itself is not meant, but simply the sign of the reduced expression, and the superscript p, j represents a variable in the j-th equation in the set of equations Θp; the 3m variables are generalized to 4K as follows:
，，/>，；
G p+j-1 and B p+j-1 are respectively equivalent conductance and susceptance of the (p+j) -1 th user, and are obtained through calculation of a regional power grid impedance matrix, wherein the regional power grid impedance matrix is obtained from a power grid dispatching center;
And/> The equivalent resistance and the equivalent reactance between the (p+j-1) th user and the kth power supply are obtained through calculation of a region power grid impedance matrix Kron reduced order elimination passive node;
the carbon emission intensity of the kth power supply is calculated for the equation set theta p;
And/> D and q axis components of the kth power supply output current calculated by the equation set Θ p respectively;
S104, solving a p-th group intermediate quantity X p through each equation set Θ p: ；
each group of intermediate quantity X p has 、/>And/>A total of 3m variables (k=1, 2, … m), and 3m equations, the unique solution can be found by numerical or least squares methods;
S105, judging whether the intermediate quantities X p are out of tolerance, if so, correcting a carbon emission metering value set Z s corresponding to the out of tolerance intermediate quantity X s to be Z s'; wherein X s is the s-th out-of-tolerance intermediate quantity.
In a second aspect, the present invention provides a system for checking measurement data of a carbon table of a user, including:
The acquisition module is used for acquiring carbon emission measurement values C i, i=1, 2, … and n of all user carbon tables of the current regional power grid in n current update periods;
A set construction module for constructing a set Z p of n-3m+1 groups of carbon emission measurement values, P=1, 2, …, n-3m+1, m is the number of power supplies generating carbon emission in the current regional power grid, m < < n, namely, the 1 st to 3m carbon emission measurement values are selected as a 1 st group carbon emission measurement value set Z 1, and the 2 nd to 3m+1, 3 rd to 3m+2, …, n-3m+1-n carbon emission measurement values are sequentially selected as a 2 nd, 3 rd, … … th and n-3m+1 th group carbon emission measurement value set Z 2、Z3、……、Zn-3m+1;
The equation set construction module is used for correspondingly forming an equation set theta p, p=1, 2, … and n-3m+1 for each set of carbon emission measurement value sets Z p;
；
Wherein:
C p+j-1 is the carbon emission measurement value of the p+j-1 th user;
t is the carbon emission factor update period;
K itself is not meant, but simply the sign of the reduced expression, and the superscript p, j represents a variable in the j-th equation in the set of equations Θp; the 3m variables are generalized to 4K as follows:
，，/>，；
G p+j-1 and B p+j-1 are respectively equivalent conductance and susceptance of the (p+j) -1 th user, and are obtained through calculation of a regional power grid impedance matrix, wherein the regional power grid impedance matrix is obtained from a power grid dispatching center;
And/> The equivalent resistance and the equivalent reactance between the (p+j-1) th user and the kth power supply are obtained through calculation of a region power grid impedance matrix Kron reduced order elimination passive node;
the carbon emission intensity of the kth power supply is calculated for the equation set theta p;
And/> D and q axis components of the kth power supply output current calculated by the equation set Θ p respectively;
a solving module for solving the p-th group intermediate quantity X p by each equation set Θ p: ；
each group of intermediate quantity X p has 、/>And/>A total of 3m variables (k=1, 2, … m), and 3m equations, the unique solution can be found by numerical or least squares methods;
The correction module is used for judging whether the intermediate quantity X p is out of tolerance, if so, correcting a carbon emission metering value set Z s corresponding to the out of tolerance intermediate quantity X s to be Z s'; wherein X s is the s-th out-of-tolerance intermediate quantity.
One or more technical solutions provided in the embodiments of the present invention at least have the following technical effects or advantages: the method has the advantages of being convenient to implement, high in sensitivity and the like, can identify and correct out-of-tolerance results in the carbon emission of the user carbon table caused by network, communication, delay factors and the like, and makes up for the blank of the lack of a checking method for the carbon emission of the current user carbon table.
The foregoing description is only an overview of the present invention, and is intended to be implemented in accordance with the teachings of the present invention in order that the same may be more clearly understood and to make the same and other objects, features and advantages of the present invention more readily apparent.
Drawings
The invention will be further described with reference to examples of embodiments with reference to the accompanying drawings.
FIG. 1 is a schematic diagram of a frame of a system of the present invention;
FIG. 2 is a flow chart of a method according to a first embodiment of the invention;
fig. 3 is a schematic structural diagram of a system in a second embodiment of the present invention.
Detailed Description
The embodiment of the application provides a method and a system for checking the metering data of a user carbon meter, wherein the method and the system calculate the intermediate value through the carbon emission of each user and check whether the out-of-tolerance occurs, and calculate the new carbon emission to correct the original result if the out-of-tolerance occurs, so as to check the metering result of the user carbon meter.
The technical scheme in the embodiment of the application has the following overall thought: and collecting the middle value of the carbon emission amount calculation of each user and checking whether the deviation occurs, if so, calculating the new carbon emission amount to correct the original result, and checking the carbon meter calculation result of the user, so that the deviation result in the carbon emission amount of the user caused by network, communication, time delay and the like can be identified and corrected, and the blank of the current carbon emission amount lack checking method of the user carbon meter is made up. Wherein, collecting the intermediate value of the carbon emission amount calculation of each user and checking whether the out-of-tolerance occurs means that: reversely calculating a plurality of groups of intermediate quantities (comprising carbon emission intensity of each emission power supply, current values of d axis and q axis) according to a carbon flow theory by using network parameters and carbon emission quantity of each user in one period; if the abnormal emission quantity value exists, abnormal intermediate quantity calculation results are caused, so that a part of intermediate quantity normal values (or a pile of intermediate quantity far away from the normal values) and a part of intermediate quantity abnormal values exist, the intermediate quantity calculation results are used as multidimensional characteristic parameters, the class of the intermediate quantity normal values and the abnormal values can be separated by adopting a support vector machine, the abnormal values are corrected by adopting the average value of the intermediate quantity normal values, and the abnormal carbon emission quantity of users is corrected.
Before describing a specific embodiment, a system framework corresponding to the method of the embodiment of the present application is described, as shown in fig. 1, the system is roughly divided into two parts:
The user carbon table is used for metering the carbon emission of the user, and due to the influences of various aspects such as network, communication and time delay, the carbon emission factor value adopted when the carbon table updates the carbon emission result can have the conditions of coding error, bad data, untimely updating and the like, so that the output carbon emission result is wrong;
And the checking system of the user carbon meter measurement data calculates the intermediate value through the carbon emission of each user and checks whether the deviation occurs, if so, calculates the new carbon emission to correct the original result, and realizes the checking of the user carbon meter measurement result.
Example 1
As shown in fig. 2, the present embodiment provides a method for checking the measurement data of the carbon table of the user, which includes the following steps:
S101, acquiring carbon emission measurement values C i, i=1, 2, … and n of all user carbon tables of a current regional power grid in n current update periods;
the current regional power grid is a region of interest to a technical application, and is as large as a certain province or a certain city, or as small as a certain community or a certain unit, and the like. Because the whole power grids are interconnected, the carbon emission factor of one user is related to the thermal power plants in the whole country theoretically, and the calculated amount is overlarge, the regional power grids need to be limited, and the equivalent is carried out on the points or connecting lines of the current regional power grids and other regional power grids, so that the carbon emission intensity of the equivalent power source is obtained.
S102, constructing a set Z of n-3m+1 groups of carbon emission measurement values p, , P=1, 2, …, n-3m+1, m is the number of power supplies generating carbon emission in the current regional power grid, m < < n, i.e. the 1 st to 3m carbon emission measurement values are selected as a1 st group carbon emission measurement value set Z 1, and the 2 nd to 3m+1, 3 rd to 3m+2, …, n-3m+1-n carbon emission measurement values are sequentially selected as a 2 nd, 3 rd, … … th, n-3m+1 th group carbon emission measurement value set Z 2、Z3、……、Zn-3m+1;
S103, correspondingly forming an equation set theta p, p=1, 2, … and n-3m+1 for each group of carbon emission measurement value sets Z p;
；
Wherein:
C p+j-1 is the carbon emission measurement value of the p+j-1 th user;
t is the carbon emission factor update period;
K itself is not meant, but simply the sign of the reduced expression, and the superscript p, j represents a variable in the j-th equation in the set of equations Θp; the 3m variables can be generalized to 4K in the derivation of the formula:
，，，/>，；
G p+j-1 and B p+j-1 are respectively equivalent conductance and susceptance of the (p+j) -1 th user, and are obtained through calculation of a regional power grid impedance matrix, wherein the regional power grid impedance matrix is obtained from a power grid dispatching center;
And/> The equivalent resistance and the equivalent reactance between the (p+j-1) th user and the kth power supply are obtained through calculation of a region power grid impedance matrix Kron reduced order elimination passive node;
the carbon emission intensity of the kth power supply is calculated for the equation set theta p;
And/> D and q axis components of the kth power supply output current calculated by the equation set Θ p respectively;
S104, solving a p-th group intermediate quantity X p through each equation set Θ p: ；
each group of intermediate quantity X p has 、/>And/>A total of 3m variables (k=1, 2, … m), and 3m equations, the unique solution can be found by numerical or least squares methods;
S105, judging whether the intermediate quantities X p are out of tolerance, if so, correcting a carbon emission metering value set Z s corresponding to the out of tolerance intermediate quantity X s to be Z s'; wherein X s is the s-th out-of-tolerance intermediate quantity.
By out-of-tolerance is meant an abnormal emission measurement, and whether or not the emission measurement is abnormal depends on the requirement of the error, such as for a carbon table with an error of 1%, calculating a result 1% (or setting a certain margin, such as 1.1%) from the average value of the normal values of the emission measurement even if the emission measurement is abnormal.
Taking m=3 as an example, if m=3, each set of carbon emission measurement values includes 9 carbon emission measurement values, then a total of n-3m+1=n-8 sets of carbon emission measurement values may be constructed, and the last set of sets of carbon emission measurement values may be denoted as Z n-8; and each system of equations consists of 9 equations.
In the step S102, the 1 st to 9 th carbon emission measurement values are selected as the 1 st group carbon emission measurement value set Z 1, the 2 nd to 10 th carbon emission measurement values are selected as the 2 nd group carbon emission measurement value set Z 2, the 3 rd to 11 th carbon emission measurement values are selected as the 3 rd group carbon emission measurement value set Z 3, and the n-8~n th carbon emission measurement values are selected as the n-8 th group carbon emission measurement value set Z n-8, then:；/>；/>；；
in the step S103, the 1 st equation set Θ 1 is configured as: ；
Equation set Θ 2 of FIG. 2 is constructed as: ；
equation set Θ 3 is constructed as: ；
The n-3m+1 th equation set Θ n-3m+1 is constructed as: ；
wherein, ，/>，，/>，
，，/>，；
In the step S104, the intermediate amounts X 1、X2、X3、……、Xn-3m+1 are respectively:
；
；
；
。
further, the step S105 specifically includes:
training the sample set { X p } by using a support vector machine, wherein p=1, 2, …, n-3m+1;
After training, the classification function f (X s) <0, s epsilon omega, wherein omega is an out-of-tolerance variable set, and the number of elements is r;
Judging whether an out-of-tolerance variable X s exists or not, if so, judging that the corresponding Z s exceeds the out-of-tolerance Z s', and determining that Corrected to/>Wherein:
；
，/>，，/>；
；/>；；/>；
In the method, in the process of the invention, The carbon emission intensity of the kth power supply is calculated for the equation set theta j; /(I)、/>D and q axis components of the kth power supply output current calculated by the equation set Θ j respectively; /(I)、/>、/>Mean value of all normal values after excluding abnormal values,/>, respectivelyRepresents the intermediate quantity of the mean value of each element of all normal values after the exception value is removed.
Based on the same inventive concept, the application also provides a system corresponding to the method in the first embodiment, and the details of the second embodiment are described in detail.
Example two
As shown in fig. 3, in this embodiment, a system for checking the measurement data of the carbon table of the user is provided, which includes:
The acquisition module is used for acquiring carbon emission measurement values C i, i=1, 2, … and n of all user carbon tables of the current regional power grid in n current update periods;
A set construction module for constructing a set Z p of n-3m+1 groups of carbon emission measurement values, P=1, 2, …, n-3m+1, m is the number of power supplies generating carbon emission in the current regional power grid, m < < n, namely, the 1 st to 3m carbon emission measurement values are selected as a 1 st group carbon emission measurement value set Z 1, and the 2 nd to 3m+1, 3 rd to 3m+2, …, n-3m+1-n carbon emission measurement values are sequentially selected as a 2 nd, 3 rd, … … th and n-3m+1 th group carbon emission measurement value set Z 2、Z3、……、Zn-3m+1;
The equation set construction module is used for correspondingly forming an equation set theta p, p=1, 2, … and n-3m+1 for each set of carbon emission measurement value sets Z p;
；
Wherein:
C p+j-1 is the carbon emission measurement value of the p+j-1 th user;
t is the carbon emission factor update period;
K itself is not meant, but simply the sign of the reduced expression, and the superscript p, j represents a variable in the j-th equation in the set of equations Θp; the 3m variables are generalized to 4K as follows:
，，/>，；
G p+j-1 and B p+j-1 are respectively equivalent conductance and susceptance of the (p+j) -1 th user, and are obtained through calculation of a regional power grid impedance matrix, wherein the regional power grid impedance matrix is obtained from a power grid dispatching center;
And/> The equivalent resistance and the equivalent reactance between the (p+j-1) th user and the kth power supply are obtained through calculation of a region power grid impedance matrix Kron reduced order elimination passive node;
the carbon emission intensity of the kth power supply is calculated for the equation set theta p;
And/> D and q axis components of the kth power supply output current calculated by the equation set Θ p respectively;
a solving module for solving the p-th group intermediate quantity X p by each equation set Θ p: ；
each group of intermediate quantity X p has 、/>And/>A total of 3m variables (k=1, 2, … m), and 3m equations, the unique solution can be found by numerical or least squares methods;
The correction module is used for judging whether the intermediate quantity X p is out of tolerance, if so, correcting a carbon emission metering value set Z s corresponding to the out of tolerance intermediate quantity X s to be Z s'; wherein X s is the s-th out-of-tolerance intermediate quantity.
Further, the correction module specifically performs the following process:
training the sample set { X p } by using a support vector machine, wherein p=1, 2, …, n-3m+1;
After training, the classification function f (X s) <0, s epsilon omega, wherein omega is an out-of-tolerance variable set, and the number of elements is r;
Judging whether an out-of-tolerance variable X s exists or not, if so, judging that the corresponding Z s exceeds the out-of-tolerance Z s', and determining that Corrected to/>Wherein:；
，/>，，/>；
；/>；；/>；
In the method, in the process of the invention, The carbon emission intensity of the kth power supply is calculated for the equation set theta j; /(I)、/>D and q axis components of the kth power supply output current calculated by the equation set Θ j respectively; /(I)、/>、/>Mean value of all normal values after excluding abnormal values,/>, respectivelyRepresents the intermediate quantity of the mean value of each element of all normal values after the exception value is removed.
Since the system described in the second embodiment of the present invention is a device for implementing the method described in the first embodiment of the present invention, based on the method described in the first embodiment of the present invention, a person skilled in the art can understand the specific structure and the modification of the device, and therefore, the description thereof is omitted herein. All devices used in the method according to the first embodiment of the present invention are within the scope of the present invention.
The technical scheme provided by the embodiment of the application has at least the following technical effects or advantages: and acquiring the middle value of the carbon emission amount calculation of each user, checking whether the middle value is out of tolerance, and calculating new carbon emission amount to correct the original result if the middle value is out of tolerance, so that the out-of-tolerance result in the carbon emission amount of the user carbon table caused by network, communication, delay factors and the like can be identified and corrected, and the blank of the lack of a checking method for the carbon emission amount of the current user carbon table is made up.
While specific embodiments of the invention have been described above, it will be appreciated by those skilled in the art that the specific embodiments described are illustrative only and not intended to limit the scope of the invention, and that equivalent modifications and variations of the invention in light of the spirit of the invention will be covered by the claims of the present invention.
Claims (2)
1. A checking method of user carbon meter measurement data is characterized in that: comprising the following steps:
S101, acquiring carbon emission measurement values C i, i=1, 2, … and n of all user carbon tables of a current regional power grid in n current update periods;
S102, constructing a set Z of n-3m+1 groups of carbon emission measurement values p, , P=1, 2, …, n-3m+1, m is the number of power supplies generating carbon emission in the current regional power grid, m < < n, i.e. the 1 st to 3m carbon emission measurement values are selected as a1 st group carbon emission measurement value set Z 1, and the 2 nd to 3m+1, 3 rd to 3m+2, …, n-3m+1-n carbon emission measurement values are sequentially selected as a 2 nd, 3 rd, … … th, n-3m+1 th group carbon emission measurement value set Z 2、Z3、……、Zn-3m+1;
S103, correspondingly forming an equation set theta p, p=1, 2, … and n-3m+1 for each group of carbon emission measurement value sets Z p;
；
Wherein:
C p+j-1 is the carbon emission measurement value of the p+j-1 th user;
t is the carbon emission factor update period;
K itself is not meant, but simply the sign of the simplified expression, the superscript p, j represents the variable in the j-th equation in the set of equations Θp, and 3m variables are generalized to 4K as follows:
，/>，，/>；
G p+j-1 and B p+j-1 are respectively equivalent conductance and susceptance of the (p+j) -1 th user, and are obtained through calculation of a regional power grid impedance matrix, wherein the regional power grid impedance matrix is obtained from a power grid dispatching center;
And/> The equivalent resistance and the equivalent reactance between the (p+j-1) th user and the kth power supply are obtained through calculation of a region power grid impedance matrix Kron reduced order elimination passive node;
the carbon emission intensity of the kth power supply is calculated for the equation set theta p;
And/> D and q axis components of the kth power supply output current calculated by the equation set Θ p respectively;
S104, solving a p-th group intermediate quantity X p through each equation set Θ p: ；
each group of intermediate quantity X p has 、/>And/>A total of 3m variables (k=1, 2, … m), and 3m equations, and a unique solution for each equation set Θ p is found;
S105, judging whether the intermediate quantities X p are out of tolerance, if yes, collecting carbon emission metering values corresponding to the out of tolerance intermediate quantity X s Modified as/>；
Wherein X s is the s-th out-of-tolerance intermediate quantity;
the step S105 specifically includes:
training the sample set { X p } by using a support vector machine, wherein p=1, 2, …, n-3m+1;
After training, the classification function f (X s) <0, s epsilon omega, wherein omega is an out-of-tolerance variable set, and the number of elements is r;
Judging whether the out-of-tolerance variable intermediate quantity X s exists or not, if so, judging that the corresponding Z s is out-of-tolerance Z s', and determining that Corrected to/>Wherein:
；
，/>，，/>；
；/>；；/>；
In the method, in the process of the invention, The carbon emission intensity of the kth power supply is calculated for the equation set theta j; /(I)、/>D and q axis components of the kth power supply output current calculated by the equation set Θ j respectively; /(I)、/>、/>Mean value of all normal values after excluding abnormal values,/>, respectivelyRepresents the intermediate quantity of the mean value of each element of all normal values after the exception value is removed.
2.A checking system of user carbon meter measurement data is characterized in that: comprising the following steps:
The acquisition module is used for acquiring carbon emission measurement values C i, i=1, 2, … and n of all user carbon tables of the current regional power grid in n current update periods;
A set construction module for constructing a set Z p of n-3m+1 groups of carbon emission measurement values, P=1, 2, …, n-3m+1, m is the number of power supplies generating carbon emission in the current regional power grid, m < < n, namely, the 1 st to 3m carbon emission measurement values are selected as a 1 st group carbon emission measurement value set Z 1, and the 2 nd to 3m+1, 3 rd to 3m+2, …, n-3m+1-n carbon emission measurement values are sequentially selected as a 2 nd, 3 rd, … … th and n-3m+1 th group carbon emission measurement value set Z 2、Z3、……、Zn-3m+1;
The equation set construction module is used for correspondingly forming an equation set theta p, p=1, 2, … and n-3m+1 for each set of carbon emission measurement value sets Z p;
；
Wherein:
C p+j-1 is the carbon emission measurement value of the p+j-1 th user;
t is the carbon emission factor update period;
K itself is not meant, but simply the sign of the reduced expression, and the superscript p, j represents a variable in the j-th equation in the set of equations Θp; the 3m variables are generalized to 4K as follows:
，/>，，/>；
G p+j-1 and B p+j-1 are respectively equivalent conductance and susceptance of the (p+j) -1 th user, and are obtained through calculation of a regional power grid impedance matrix, wherein the regional power grid impedance matrix is obtained from a power grid dispatching center;
And/> The equivalent resistance and the equivalent reactance between the (p+j-1) th user and the kth power supply are obtained through calculation of a region power grid impedance matrix Kron reduced order elimination passive node;
the carbon emission intensity of the kth power supply is calculated for the equation set theta p;
And/> D and q axis components of the kth power supply output current calculated by the equation set Θ p respectively;
a solving module for solving the p-th group intermediate quantity X p by each equation set Θ p: ；
each group of intermediate quantity X p has 、/>And/>A total of 3m variables (k=1, 2, … m), and 3m equations, and a unique solution for each equation set Θ p is found;
The correction module is used for judging whether the intermediate quantities X p are out of tolerance, if so, correcting a carbon emission metering value set Z s corresponding to the out of tolerance intermediate quantity X s to be Z s', wherein X s is the s-th out of tolerance intermediate quantity;
The correction module specifically executes the following processes:
training the sample set { X p } by using a support vector machine, wherein p=1, 2, …, n-3m+1;
After training, the classification function f (X s) <0, s epsilon omega, wherein omega is an out-of-tolerance variable set, and the number of elements is r;
Judging whether the out-of-tolerance intermediate quantity X s exists or not, if so, judging that the corresponding Z s is out-of-tolerance Z s', and Corrected to/>Wherein:；
，/>，，/>；
；
；/>；；
In the method, in the process of the invention, The carbon emission intensity of the kth power supply is calculated for the equation set theta j; /(I)、/>D and q axis components of the kth power supply output current calculated by the equation set Θ j respectively; /(I)、/>、/>Mean value of all normal values after excluding abnormal values,/>, respectivelyRepresents the intermediate quantity of the mean value of each element of all normal values after the exception value is removed.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202410194801.XA CN117764798B (en) | 2024-02-22 | 2024-02-22 | Method and system for checking carbon meter measurement data of user |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202410194801.XA CN117764798B (en) | 2024-02-22 | 2024-02-22 | Method and system for checking carbon meter measurement data of user |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN117764798A CN117764798A (en) | 2024-03-26 |
| CN117764798B true CN117764798B (en) | 2024-05-24 |
Family
ID=90322395
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202410194801.XA Active CN117764798B (en) | 2024-02-22 | 2024-02-22 | Method and system for checking carbon meter measurement data of user |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN117764798B (en) |
Citations (13)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN106251095A (en) * | 2016-09-06 | 2016-12-21 | 清华大学 | A kind of method of power system carbon emission real time measure and carbon table system |
| CN112816652A (en) * | 2021-02-23 | 2021-05-18 | 国电环境保护研究院有限公司 | Device for testing fuel oxidation rate and carbon emission of combustion engine |
| CN114358971A (en) * | 2022-01-07 | 2022-04-15 | 清华大学 | Carbon emission metering method and device for natural gas system, electronic equipment and storage medium |
| CN114662781A (en) * | 2022-04-07 | 2022-06-24 | 南京师范大学 | Cloud-based intelligent carbon footprint evaluation management system and evaluation method thereof |
| CN115392528A (en) * | 2022-06-06 | 2022-11-25 | 安徽大学 | Carbon emission fine metering method for electricity utilization side based on carbon emission flow theory |
| WO2023160283A1 (en) * | 2022-02-24 | 2023-08-31 | 清华大学 | Carbon emission response calculation method and apparatus based on carbon emission flow |
| CN117252490A (en) * | 2023-07-31 | 2023-12-19 | 安徽南瑞中天电力电子有限公司 | Carbon emission accounting method, system and device based on special transformer acquisition terminal |
| CN117273276A (en) * | 2023-10-10 | 2023-12-22 | 南方电网能源发展研究院有限责任公司 | Carbon emission monitoring method and device based on electric power data |
| CN117272362A (en) * | 2023-06-12 | 2023-12-22 | 中亿丰数字科技集团有限公司 | Block chain-based carbon emission perfect management method and system |
| CN117371650A (en) * | 2023-10-09 | 2024-01-09 | 国网江苏省电力有限公司连云港供电分公司 | Accurate carbon metering method and system for power distribution network considering load side electric energy substitution |
| CN117412201A (en) * | 2023-10-16 | 2024-01-16 | 国网江苏省电力有限公司营销服务中心 | Carbon metering network system and carbon metering data optimization acquisition method |
| WO2024015014A1 (en) * | 2022-07-12 | 2024-01-18 | Envision Digital International Pte. Ltd. | Method and apparatus for displaying carbon intensities, and device, storage medium, and program product thereof |
| CN117557111A (en) * | 2023-09-28 | 2024-02-13 | 国网四川省电力公司电力科学研究院 | Method and system for calculating carbon emission factor of coal-fired power plant for real-time carbon transaction |
-
2024
- 2024-02-22 CN CN202410194801.XA patent/CN117764798B/en active Active
Patent Citations (13)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN106251095A (en) * | 2016-09-06 | 2016-12-21 | 清华大学 | A kind of method of power system carbon emission real time measure and carbon table system |
| CN112816652A (en) * | 2021-02-23 | 2021-05-18 | 国电环境保护研究院有限公司 | Device for testing fuel oxidation rate and carbon emission of combustion engine |
| CN114358971A (en) * | 2022-01-07 | 2022-04-15 | 清华大学 | Carbon emission metering method and device for natural gas system, electronic equipment and storage medium |
| WO2023160283A1 (en) * | 2022-02-24 | 2023-08-31 | 清华大学 | Carbon emission response calculation method and apparatus based on carbon emission flow |
| CN114662781A (en) * | 2022-04-07 | 2022-06-24 | 南京师范大学 | Cloud-based intelligent carbon footprint evaluation management system and evaluation method thereof |
| CN115392528A (en) * | 2022-06-06 | 2022-11-25 | 安徽大学 | Carbon emission fine metering method for electricity utilization side based on carbon emission flow theory |
| WO2024015014A1 (en) * | 2022-07-12 | 2024-01-18 | Envision Digital International Pte. Ltd. | Method and apparatus for displaying carbon intensities, and device, storage medium, and program product thereof |
| CN117272362A (en) * | 2023-06-12 | 2023-12-22 | 中亿丰数字科技集团有限公司 | Block chain-based carbon emission perfect management method and system |
| CN117252490A (en) * | 2023-07-31 | 2023-12-19 | 安徽南瑞中天电力电子有限公司 | Carbon emission accounting method, system and device based on special transformer acquisition terminal |
| CN117557111A (en) * | 2023-09-28 | 2024-02-13 | 国网四川省电力公司电力科学研究院 | Method and system for calculating carbon emission factor of coal-fired power plant for real-time carbon transaction |
| CN117371650A (en) * | 2023-10-09 | 2024-01-09 | 国网江苏省电力有限公司连云港供电分公司 | Accurate carbon metering method and system for power distribution network considering load side electric energy substitution |
| CN117273276A (en) * | 2023-10-10 | 2023-12-22 | 南方电网能源发展研究院有限责任公司 | Carbon emission monitoring method and device based on electric power data |
| CN117412201A (en) * | 2023-10-16 | 2024-01-16 | 国网江苏省电力有限公司营销服务中心 | Carbon metering network system and carbon metering data optimization acquisition method |
Non-Patent Citations (3)
| Title |
|---|
| Modeling carbon dioxide emissions with a system of differential equations;Tsokos, Chris P.;《NONLINEAR ANALYSIS-THEORY METHODS & APPLICATIONS》;20091215;第71卷(第12期);第1182-1197页 * |
| 基于碳排放流迭代算法的分布式碳表系统（一）：理论方法与分析;李业辉 等;《电网技术》;20230630;第47卷(第06期);第2165-2173页 * |
| 电力系统全环节碳计量方法与碳表系统;张宁 等;《电力系统自动化》;20230510;第47卷(第09期);第2-12页 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN117764798A (en) | 2024-03-26 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN109754013B (en) | Electric power system hybrid measurement fusion method based on unscented Kalman filtering | |
| CN102521080B (en) | Computer data recovery method for electricity-consumption information collecting system for power consumers | |
| CN113985339B (en) | Error diagnosis method and system for intelligent ammeter, equipment and storage medium | |
| CN105207204B (en) | One kind meter and the probabilistic Probabilistic Load Flow analysis method of primary frequency modulation | |
| CN102692615A (en) | System capable of automatically acquiring electric quantity data | |
| CN108336739A (en) | A kind of Probabilistic Load Flow on-line calculation method based on RBF neural | |
| CN108448568A (en) | Power distribution network admixture method of estimation based on a variety of time cycle measurement data | |
| CN112130109A (en) | Method for detecting metering performance abnormity of intelligent electric energy meter | |
| CN103900666A (en) | Long-time running measuring precision control method of membrane type gas meter | |
| CN110348720A (en) | The electricity quality evaluation method of rural area photovoltaic access system | |
| CN115616473B (en) | Identification method, device, equipment and storage medium of out-of-tolerance electric energy meter | |
| CN105262146B (en) | Method and system for calculating reserve capacity of power system containing wind power | |
| CN112508338A (en) | New energy station grid-connected performance comprehensive evaluation system and evaluation method thereof | |
| CN115618654A (en) | Identification method and device for out-of-tolerance electric energy meter | |
| CN117764798B (en) | Method and system for checking carbon meter measurement data of user | |
| CN102780220A (en) | Power flow calculation method for electric distribution network comprising PV constant distributed generation | |
| CN105956755A (en) | Method and system for establishing quantitative relationship of general power line loss rate influencing factors | |
| CN110518627A (en) | Meter and the probabilistic power system steady state voltage stability probability evaluation method of failure of wind-powered electricity generation | |
| CN117913974A (en) | Daily electric quantity monitoring and analyzing system | |
| CN107146166B (en) | Distribution network planning auxiliary system | |
| CN117033356A (en) | Real-time data quality optimization method for power enterprise production | |
| CN114707860B (en) | ECI-based user electricity consumption information management system | |
| CN111461297A (en) | Solar irradiation quantity optimization prediction algorithm based on MPC and E L M neural network | |
| CN115776135A (en) | Distributed power supply sensing method based on state estimation | |
| CN116805029A (en) | Earth potential distribution calculation method, device and system based on actual measurement of direct current magnetic bias current of power grid and storage medium |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |