CN116797245A - Power utilization side carbon emission coefficient allocation method and system - Google Patents
Power utilization side carbon emission coefficient allocation method and system Download PDFInfo
- Publication number
- CN116797245A CN116797245A CN202310793108.XA CN202310793108A CN116797245A CN 116797245 A CN116797245 A CN 116797245A CN 202310793108 A CN202310793108 A CN 202310793108A CN 116797245 A CN116797245 A CN 116797245A
- Authority
- CN
- China
- Prior art keywords
- carbon
- carbon emission
- load
- power
- node
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Pending
Links
- OKTJSMMVPCPJKN-UHFFFAOYSA-N Carbon Chemical compound [C] OKTJSMMVPCPJKN-UHFFFAOYSA-N 0.000 title claims abstract description 291
- 229910052799 carbon Inorganic materials 0.000 title claims abstract description 291
- 238000000034 method Methods 0.000 title claims abstract description 39
- 239000011159 matrix material Substances 0.000 claims abstract description 78
- 238000002347 injection Methods 0.000 claims abstract description 47
- 239000007924 injection Substances 0.000 claims abstract description 47
- 230000005611 electricity Effects 0.000 claims abstract description 37
- 238000004401 flow injection analysis Methods 0.000 claims abstract description 26
- 238000010248 power generation Methods 0.000 claims abstract description 18
- 239000013598 vector Substances 0.000 claims description 17
- 238000004458 analytical method Methods 0.000 claims description 16
- 230000008878 coupling Effects 0.000 abstract description 2
- 238000010168 coupling process Methods 0.000 abstract description 2
- 238000005859 coupling reaction Methods 0.000 abstract description 2
- 230000006399 behavior Effects 0.000 description 7
- 238000005457 optimization Methods 0.000 description 3
- 239000000243 solution Substances 0.000 description 3
- CURLTUGMZLYLDI-UHFFFAOYSA-N Carbon dioxide Chemical compound O=C=O CURLTUGMZLYLDI-UHFFFAOYSA-N 0.000 description 2
- 230000005540 biological transmission Effects 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- VNWKTOKETHGBQD-UHFFFAOYSA-N methane Chemical compound C VNWKTOKETHGBQD-UHFFFAOYSA-N 0.000 description 2
- 230000000903 blocking effect Effects 0.000 description 1
- 229910002092 carbon dioxide Inorganic materials 0.000 description 1
- 239000001569 carbon dioxide Substances 0.000 description 1
- 239000003245 coal Substances 0.000 description 1
- 238000002485 combustion reaction Methods 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 238000005265 energy consumption Methods 0.000 description 1
- 230000007613 environmental effect Effects 0.000 description 1
- 230000000977 initiatory effect Effects 0.000 description 1
- 238000002372 labelling Methods 0.000 description 1
- 239000003345 natural gas Substances 0.000 description 1
- 230000000750 progressive effect Effects 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q30/00—Commerce
- G06Q30/018—Certifying business or products
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q50/00—Systems or methods specially adapted for specific business sectors, e.g. utilities or tourism
- G06Q50/06—Electricity, gas or water supply
Landscapes
- Business, Economics & Management (AREA)
- Engineering & Computer Science (AREA)
- Economics (AREA)
- General Physics & Mathematics (AREA)
- Health & Medical Sciences (AREA)
- Theoretical Computer Science (AREA)
- Marketing (AREA)
- Strategic Management (AREA)
- Physics & Mathematics (AREA)
- General Business, Economics & Management (AREA)
- Finance (AREA)
- Development Economics (AREA)
- Accounting & Taxation (AREA)
- Entrepreneurship & Innovation (AREA)
- Public Health (AREA)
- Water Supply & Treatment (AREA)
- General Health & Medical Sciences (AREA)
- Human Resources & Organizations (AREA)
- Primary Health Care (AREA)
- Tourism & Hospitality (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
The invention discloses a method and a system for allocating carbon emission coefficients of a power utilization side, and relates to the field of electric carbon coupling markets, wherein the method comprises the steps of obtaining injection power of each node in each period according to a power generation side quotation of a target area and historical power of each node in each period; the nodes comprise generator nodes and load nodes; obtaining a load side carbon flow injection matrix according to the injection power of all the nodes in each period; obtaining a total carbon quota coefficient of the target area according to the load side carbon flow injection matrix, and obtaining carbon emission influence factors of each load node; and obtaining the carbon emission coefficient of the load node according to the carbon quota total coefficient and the carbon emission influence factor of each load node. According to the invention, the influence degree of the user on the total carbon emission amount of the target area is considered, and the accurate tracking and tracing of the indirect carbon emission of the user side are realized, so that the fair allocation of the carbon emission coefficient is realized, and the user is stimulated to improve the self electricity consumption behavior.
Description
Technical Field
The invention relates to the field of electric carbon coupling markets, in particular to a power consumption side carbon emission coefficient allocation method and system based on electric power network carbon traceability analysis.
Background
When non-renewable energy sources such as coal, natural gas and the like are adopted for power generation, a large amount of carbon dioxide is generated after the combustion, and carbon emission is formed. At the same time of consuming electric energy, the user side has different carbon emission responsibilities allocated to different load nodes due to factors such as space distribution difference, electricity consumption behavior difference, power transmission blocking problem and the like of each load node, and the carbon emission coefficients corresponding to the electricity consumption behavior are different. It is therefore necessary to consider how to fairly and reasonably meter the carbon emission responsibilities of the user side and achieve a reasonable apportionment of the carbon emission coefficients. In recent years, with the formal initiation of carbon market trading, power generation enterprises have been required to purchase carbon quotas from the carbon market to meet power generation demands. The environmental coefficient of carbon emissions has hysteresis to users by electricity price conduction, and although there is a certain effect in reducing carbon emissions, it is difficult to guide the user-side electricity behavior. Considering that the power system 'source follow-up' feature makes the load side the main responsible person of the carbon emission of the power system, it is necessary to consider the coefficient allocation method of the carbon emission from the viewpoint of the electricity utilization side.
Disclosure of Invention
The invention aims to provide a method and a system for distributing carbon emission coefficients at an electricity utilization side, which can realize fair distribution of the carbon emission coefficients.
In order to achieve the above object, the present invention provides the following solutions:
a method for apportioning carbon emission coefficients on an electricity utilization side, the method comprising:
calculating the injection power of each node in each period according to the power generation side quotation of the target area and the historical power of each node in each period; the nodes comprise generator nodes and load nodes;
calculating a load side carbon flow injection matrix according to the injection power of all the nodes in each period; the load side carbon stream injection matrix includes carbon emission streams for each of the load nodes at each time period;
calculating a total carbon quota coefficient of the target area and a carbon emission influence factor of each load node according to the load side carbon flow injection matrix; the carbon emission influence factor characterizes the influence degree of the electricity consumption of a user on the carbon emission of the electric power system;
for each load node, calculating a carbon emission coefficient of the load node according to the carbon quota total coefficient and a carbon emission influence factor of the load node.
Optionally, the calculating the load side carbon flow injection matrix according to the injection power of all the nodes in each period specifically includes:
and taking the injection power of all the nodes in each period as input, and calculating a load side carbon flow injection matrix by using a carbon flow tracking formula.
Optionally, the method for establishing the carbon flow tracking formula is as follows:
establishing a direct current power flow model and a power phase angle incidence matrix of a load node; the expression of the direct current power flow model is as follows:
wherein P is G Is connected with the generatorGenerator power train vector, P, consisting of injection power at each time period L A load power column vector consisting of the injection power of the load node at each time period; b is a node susceptance matrix; b (B) GG ，B LG ，B GL And B LL The node susceptance matrix B is a partitioned matrix which is distinguished according to generator nodes and load nodes; θ G A phase angle column vector for an engine node; θ L Phase angle column vectors for load nodes;
the expression of the power phase angle incidence matrix of the load node is thatWherein (1)>For the elements of the ith row and jth column of the power phase angle correlation matrix of the load node, +.>And->Respectively P L And theta L Is the i-th element of (a);
and deducing the carbon flow tracking formula according to the direct current power flow model and the power phase angle incidence matrix of the load node.
Optionally, the carbon flow tracking formula is:
wherein C is L Representing a load side carbon stream injection matrix; q (Q) LG For theta L And theta G Is a correlation matrix of (a); e (E) G For a generator carbon emission factor matrix composed of average carbon emission factors of the generator nodes in each period, elements on a diagonal line represent the average carbon emission factors of the generator nodes; mu (mu) G Is P G And theta G Is a correlation matrix of (a); Γ -shaped structure L Is the power phase angle correlation matrix of the load node.
Optionally, the expression of the carbon quota total coefficient is:
wherein pi s Is the total coefficient of carbon quota;representing the carbon emission flow of load node l during the t-th period; t is all time periods; l is the total number of load nodes; r is R s Representing the free carbon quota allocated to the target area power industry.
Optionally, the expression of the carbon emission influencing factor is:
wherein τ l Carbon emission influencing factors for load nodes l;representing the carbon emission flow of load node l during the t-th period; t is all time periods; l is the total number of load nodes.
Optionally, the calculating the carbon emission coefficient of the load node according to the carbon quota total coefficient and the carbon emission influence factor of the load node specifically includes:
and taking the total carbon quota coefficient and the carbon emission influence factor of the load node as inputs, and calculating the carbon emission coefficient of the load node by using a carbon emission coefficient analysis formula.
Optionally, the carbon emission coefficient analysis formula is:
u l ＝τ l π s ；
wherein u is l Carbon emission coefficient for load node l; τ l Carbon emission influencing factor for load node l；π s Is the carbon quota total coefficient.
Optionally, the method for establishing the carbon emission coefficient analysis formula comprises the following steps:
to be used forFor the objective function, in->Establishing a carbon emission coefficient allocation model for constraint conditions; wherein u is l Carbon emission coefficient for load node l; τ 1 Carbon emission influencing factors for load nodes l; pi s Is the total coefficient of carbon quota; l is the total number of load nodes;
and deducing the carbon emission coefficient allocation model to obtain a carbon emission coefficient analysis formula.
The invention also provides a system for apportioning the carbon emission coefficient of the electricity utilization side, which comprises:
the power prediction module is used for calculating the injection power of each node in each period according to the power generation side quotation of the target area and the historical power of each node in each period; the nodes comprise generator nodes and load nodes;
the load side carbon flow injection matrix determining module is used for calculating a load side carbon flow injection matrix according to the injection power of all the nodes in each period; the load side carbon stream injection matrix includes carbon emission streams for each of the load nodes at each time period;
the carbon emission parameter determining module is used for calculating the total carbon quota coefficient of the target area and the carbon emission influence factor of each load node according to the load side carbon flow injection matrix; the carbon emission influence factor characterizes the influence degree of the electricity consumption of a user on the carbon emission of the electric power system;
and the carbon transaction coefficient allocation module is used for calculating the carbon emission coefficient of each load node according to the total carbon quota coefficient and the carbon emission influence factor of the load node.
According to the specific embodiment provided by the invention, the invention discloses the following technical effects: the invention provides a method and a system for allocating carbon emission coefficients of a power utilization side, wherein the method comprises the following steps: obtaining the injection power of each node in each period according to the power generation side quotation of the target area and the historical power of each node in each period; the nodes comprise generator nodes and load nodes; obtaining a load side carbon flow injection matrix according to the injection power of all the nodes in each period; the load side carbon stream injection matrix comprises carbon emission streams of each load node in each period; obtaining a total carbon quota coefficient of the target area according to the load side carbon flow injection matrix, and obtaining carbon emission influence factors of each load node; the carbon emission influence factor characterizes the influence degree of the electricity consumption of the user on the carbon emission of the electric power system; and for each load node, obtaining the carbon emission coefficient of the load node according to the total carbon quota coefficient and the carbon emission influence factor of the load node. According to the invention, the influence degree (carbon emission influence factor) of the user on the total carbon emission amount of the target area is considered, so that the indirect carbon emission of the user side is tracked and traced in a refined way, the fair allocation of the carbon emission coefficient is realized, and the user is stimulated to improve the self electricity consumption behavior.
Drawings
In order to more clearly illustrate the embodiments of the present invention or the technical solutions of the prior art, the drawings that are needed in the embodiments will be briefly described below, it being obvious that the drawings in the following description are only some embodiments of the present invention, and that other drawings may be obtained according to these drawings without inventive effort for a person skilled in the art.
FIG. 1 is a schematic flow chart of a method for allocating carbon emission coefficients on an electricity consumption side according to an embodiment of the present invention;
FIG. 2 is a schematic diagram of an implementation flow of a method for allocating carbon emission coefficients of an electricity consumption side according to an embodiment of the present invention;
fig. 3 is a block diagram of a power consumption side carbon emission coefficient allocation system according to an embodiment of the present invention.
Detailed Description
The following description of the embodiments of the present invention will be made clearly and completely with reference to the accompanying drawings, in which it is apparent that the embodiments described are only some embodiments of the present invention, but not all embodiments. All other embodiments, which can be made by those skilled in the art based on the embodiments of the invention without making any inventive effort, are intended to be within the scope of the invention.
The carbon emission flow theory of the electric power system realizes the responsibility allocation of the carbon emission of the electric power by labeling the electric power flow with carbon labels. The theory can distribute the direct carbon emission generated by the source side power generation to the network side and the load side in real time according to the source side power generation information, the local new energy consumption information, the line tide information, the line network loss information and the like, and distinguish the electricity consumption carbon emission responsibility difference of different electricity consumption behaviors of different users.
The invention aims to provide a power consumption side carbon emission coefficient allocation method and system, a tracking and tracing method (carbon flow tracking formula) of a user side carbon emission flow is established based on deduction of a direct current trend model in a circuit theory, and the coefficient allocation method considering user side carbon tracing analysis is provided. The carbon flow tracking method provided by the invention does not depend on a proportional distribution principle, and accords with the actual physical operation rule of the power grid. The influence degree (carbon emission influence factor) of the user on the total carbon emission amount of the target area is considered, the indirect carbon emission caused by the electricity consumption behavior of the user side can be accurately metered through tracking the carbon flow in the electric power system, and the accurate tracking and tracing of the indirect carbon emission of the user side are realized, so that the fair allocation of the carbon emission coefficient is realized, the user is further stimulated to improve the electricity consumption behavior of the user, the carbon emission amount of the user is reduced, and theoretical guidance is provided for realizing the carbon reduction of the electric power system.
In order that the above-recited objects, features and advantages of the present invention will become more readily apparent, a more particular description of the invention will be rendered by reference to the appended drawings and appended detailed description.
As shown in fig. 1, the invention provides a method for allocating carbon emission coefficients of an electricity utilization side, which comprises the following steps:
s1: calculating the injection power of each node in each period according to the power generation side quotation of the target area and the historical power of each node in each period; the nodes include generator nodes and load nodes. The power generation side quotation includes a declared electricity price and declared electricity quantity of the power generation side.
S2: calculating a load side carbon flow injection matrix according to the injection power of all the nodes in each period; the load side carbon stream injection matrix includes a carbon emission stream for each of the load nodes at each time period.
S3: calculating a total carbon quota coefficient of the target area and a carbon emission influence factor of each load node according to the load side carbon flow injection matrix; the carbon emission influence factor characterizes influence degree of electricity consumption of a user on carbon emission of the electric power system.
S4: for each load node, calculating a carbon emission coefficient of the load node according to the carbon quota total coefficient and a carbon emission influence factor of the load node.
The method for calculating the load side carbon flow injection matrix according to the injection power of all the nodes in each period specifically comprises the following steps:
and taking the injection power of all the nodes in each period as input, and calculating a load side carbon flow injection matrix by using a carbon flow tracking formula. Specifically:
in this embodiment, first, system optimization scheduling needs to be performed according to the declared electricity price and declared electricity quantity of the power generation side and the load power scene prediction of each period, that is, the injection power of each node in each period is calculated according to the power generation side quotation of the target area and the historical power of each node (including the generator node and the load node) in each period, specifically, the injection power of each load node in each period is predicted through the historical power of the load node in each period, and then the injection power of each generator node in each period is obtained through the system optimization scheduling calculation. The present embodiment divides each of 24 hours a day into 24 periods in total, as one period. The injection power of each node in each period is predicted by taking the power of each node in a certain day in the history as the historical power and the quotation of the power generation side. Treating the net inflow node as a generator node, and treating the net outflowThe node is regarded as a load node, and a positive power column vector P of the generator node injection system of each period is obtained G And a negative power column vector P for each period load node injection system L . The power load prediction is one of important works of the power department, accurate load prediction can be achieved, starting and stopping of a generator set in the power grid can be economically and reasonably arranged, and safety and stability of operation of the power grid are maintained.
In the process of establishing the carbon flow tracking model of the electric power system, the carbon flow tracking formula is established according to the obtained power system optimization scheduling result, the power and phase angle correlation expression under the direct current power flow model is firstly established, then the correlation expression of the injection power and the phase angle of the load node is constructed, the dynamic carbon emission factor of the generator node is introduced, and the carbon flow tracking expression (carbon flow tracking formula) is deduced.
The establishment method of the carbon flow tracking formula comprises the following steps:
establishing a direct current power flow model and a power phase angle incidence matrix of a load node; the direct current power flow model is established through a power grid topological structure, and the expression of the direct current power flow model is as follows:
wherein P is G For a generator power train vector consisting of the injection power of the generator node at each time period, P L A load power column vector consisting of the injection power of the load node at each time period; b is a node susceptance matrix; b (B) GG ，B LG ，B GL And B LL For a partitioned matrix in node susceptance matrix B, which is distinguished by generator nodes and load nodes, elements of M rows and M columns before matrix B form B GG Elements of the first M rows and the last N columns form B GL The elements of the first M columns of the last N rows form B LG Element constitution B of N columns and N rows LL M is the number of generator nodes, and N is the number of load nodes. θ G A phase angle column vector for an engine node; θ L Is the loadPhase angle column vector of the node.
The expression of the power phase angle incidence matrix of the load node is thatWherein the power phase angle associated matrix of the load node is a pair of angle matrix, < >>For the elements of the ith row and jth column of the power phase angle correlation matrix of the load node, +.>And->Respectively P L And theta L Is the i-th element of the load node, the power phase angle correlation matrix Γ L Is P L And theta L Is used for the correlation matrix of the (a).
And deducing the carbon flow tracking formula according to the direct current power flow model and the power phase angle incidence matrix of the load node. Substituting the power phase angle incidence matrix of the load node into a direct current power flow model to obtain a carbon flow tracking formula, wherein the carbon flow tracking formula is as follows:
wherein C is L Representing a load side carbon stream injection matrix; q (Q) LG For theta L And theta G Is a correlation matrix of (a); e (E) G For a generator carbon emission factor matrix composed of average carbon emission factors of generator nodes in each period, the generator carbon emission factor matrix is a pair of angle matrixes, and elements on the diagonal represent the average carbon emission factors (dynamic carbon emission factors) of the generator nodes, and the carbon emission factors E of the generator nodes in each period G Positive power column vector P with generator node injection system G The carbon emission of the generator in each period can be obtained through multiplication; mu (mu) G Is P G And theta G Is a correlation matrix of (a); Γ -shaped structure L Is the power phase angle correlation matrix of the load node.
According to the carbon flow tracking formula, the power column vector P of the generator G Load power column vector P L And susceptance parameters (node susceptance matrix) of the power transmission line, and a load side carbon flow injection matrix C can be calculated and obtained L Which includes a carbon emission stream per time period for each of the load nodes.
According to the load side carbon flow injection matrix C L The coefficient allocation model (namely the carbon emission coefficient allocation model) based on Nash bargaining theory is constructed, and the concrete steps are as follows:
in the construction of the carbon emission coefficient allocation model, the calculation result of the carbon emission flow at the user side is considered, the influence factor (carbon emission influence factor of the load node) of the electricity consumption of the user on the system carbon emission is embedded into the carbon emission coefficient allocation model, and the analysis solution (carbon emission coefficient analysis formula) of the carbon emission coefficient allocation model is obtained, so that the influence degree of the user on the total regional carbon emission amount can be distinguished when the carbon emission coefficient is allocated, and the fair carbon emission coefficient allocation result is obtained.
The expression of the total carbon quota coefficient corresponding to the power consumption of the target area is as follows:
wherein pi s Is the total coefficient of carbon quota;representing the carbon emission flow of load node l during the t-th period; t is all time periods; l is the total number of load nodes; r is R s Representing the free carbon quota allocated to the target area power industry.
The expression of the carbon emission influencing factor is:
wherein τ l Carbon emission influencing factors for load nodes l;representing the carbon emission flow of load node l during the t-th period; t is all time periods; l is the total number of load nodes.
In step S4, calculating a carbon emission coefficient of the load node according to the carbon quota total coefficient and the carbon emission influence factor of the load node specifically includes:
and taking the total carbon quota coefficient and the carbon emission influence factor of the load node as inputs, and calculating the carbon emission coefficient of the load node by using a carbon emission coefficient analysis formula.
The establishment method of the carbon emission coefficient analysis formula is as follows:
to be used forFor the objective function, in->Establishing a carbon emission coefficient allocation model for constraint conditions; wherein u is l Carbon emission coefficient for load node l; τ l Carbon emission influencing factors for load nodes l; pi s Is the total coefficient of carbon quota; l is the total number of load nodes;
and deducing the carbon emission coefficient allocation model to obtain a carbon emission coefficient analysis formula.
The analytical formula of the carbon emission coefficient is as follows:
u l ＝τ l π s 。
wherein u is l Carbon emission coefficient for load node l; τ l Carbon emission influencing factors for load nodes l; pi s Is the carbon quota total coefficient.
And obtaining the carbon emission coefficient allocated to the user where each load node is located through the carbon emission coefficient analysis formula, the obtained carbon quota total coefficient and the carbon emission influence factor of each load node.
As shown in fig. 3, the present invention also provides a power-consumption-side carbon emission coefficient allocation system, which includes:
the power prediction module T1 is used for calculating the injection power of each node in each period according to the power generation side quotation of the target area and the historical power of each node in each period; the nodes include generator nodes and load nodes.
A load side carbon flow injection matrix determining module T2, configured to calculate a load side carbon flow injection matrix according to injection powers of all the nodes in each period; the load side carbon stream injection matrix includes a carbon emission stream for each of the load nodes at each time period.
A carbon emission parameter determining module T3, configured to calculate a total carbon quota coefficient of the target area and a carbon emission influence factor of each load node according to the load side carbon flow injection matrix; the carbon emission influence factor characterizes influence degree of electricity consumption of a user on carbon emission of the electric power system.
And the carbon transaction coefficient allocation module T4 is used for calculating the carbon emission coefficient of each load node according to the total carbon quota coefficient and the carbon emission influence factor of the load node.
In the present specification, each embodiment is described in a progressive manner, and each embodiment is mainly described in a different point from other embodiments, and identical and similar parts between the embodiments are all enough to refer to each other. For the power system disclosed in the embodiment, since the power system corresponds to the method disclosed in the embodiment, the description is relatively simple, and the relevant points are referred to in the description of the method section.
The principles and embodiments of the present invention have been described herein with reference to specific examples, the description of which is intended only to assist in understanding the methods of the present invention and the core ideas thereof; also, it is within the scope of the present invention to be modified by those of ordinary skill in the art in light of the present teachings. In view of the foregoing, this description should not be construed as limiting the invention.
Claims (10)
1. A method for apportioning carbon emission coefficients of an electricity utilization side, the method comprising:
calculating the injection power of each node in each period according to the power generation side quotation of the target area and the historical power of each node in each period; the nodes comprise generator nodes and load nodes;
calculating a load side carbon flow injection matrix according to the injection power of all the nodes in each period; the load side carbon stream injection matrix includes carbon emission streams for each of the load nodes at each time period;
calculating a total carbon quota coefficient of the target area and a carbon emission influence factor of each load node according to the load side carbon flow injection matrix; the carbon emission influence factor characterizes the influence degree of the electricity consumption of a user on the carbon emission of the electric power system;
for each load node, calculating a carbon emission coefficient of the load node according to the carbon quota total coefficient and a carbon emission influence factor of the load node.
2. The electricity side carbon emission factor allocation method according to claim 1, wherein the calculating the load side carbon stream injection matrix according to the injection power of all the nodes in each period specifically comprises:
and taking the injection power of all the nodes in each period as input, and calculating a load side carbon flow injection matrix by using a carbon flow tracking formula.
3. The electricity consumption side carbon emission coefficient allocation method according to claim 2, wherein the establishment method of the carbon flow tracking formula is as follows:
establishing a direct current power flow model and a power phase angle incidence matrix of a load node; the expression of the direct current power flow model is as follows:
wherein P is G For a generator power train vector consisting of the injection power of the generator node at each time period, P L A load power column vector consisting of the injection power of the load node at each time period; b is a node susceptance matrix; b (B) GG ，B LG ，B GL And B LL The node susceptance matrix B is a partitioned matrix which is distinguished according to generator nodes and load nodes; θ G A phase angle column vector for an engine node; θ L Phase angle column vectors for load nodes;
the expression of the power phase angle incidence matrix of the load node is thatWherein (1)>For the elements of the ith row and jth column of the power phase angle correlation matrix of the load node, +.>And->Respectively P L And theta L Is the i-th element of (a);
and deducing the carbon flow tracking formula according to the direct current power flow model and the power phase angle incidence matrix of the load node.
4. The electricity side carbon emission factor allocation method according to claim 3, wherein the carbon flow tracking formula is:
wherein C is L Representing a load side carbon stream injection matrix; q (Q) LG For theta L And theta G Is a correlation matrix of (a); e (E) G For a generator carbon emission factor matrix composed of average carbon emission factors of the generator nodes in each period, elements on a diagonal line represent the average carbon emission factors of the generator nodes; mu (mu) G Is P G And theta G Is a correlation matrix of (a); Γ -shaped structure L Is the power phase angle correlation matrix of the load node.
5. The electricity consumption side carbon emission coefficient allocation method according to claim 1, wherein the expression of the carbon quota total coefficient is:
wherein pi s Is the total coefficient of carbon quota;representing the carbon emission flow of load node l during the t-th period; t is all time periods; l is the total number of load nodes; r is R s Representing the free carbon quota allocated to the target area power industry.
6. The electricity consumption side carbon emission factor allocation method according to claim 1, wherein the expression of the carbon emission influence factor is:
wherein τ l Carbon emission influencing factors for load nodes l;representing the carbon emission flow of load node l during the t-th period; t is all time periods; l is the total number of load nodes.
7. The electricity-side carbon emission factor allocation method according to claim 1, wherein the calculating the carbon emission factor of the load node according to the carbon quota total factor and the carbon emission influence factor of the load node specifically includes:
and taking the total carbon quota coefficient and the carbon emission influence factor of the load node as inputs, and calculating the carbon emission coefficient of the load node by using a carbon emission coefficient analysis formula.
8. The electricity consumption side carbon emission coefficient allocation method according to claim 7, wherein the carbon emission coefficient analysis formula is:
u l ＝τ l π s ；
wherein u is l Carbon emission coefficient for load node l; τ l Carbon emission influencing factors for load nodes l; pi s Is the carbon quota total coefficient.
9. The electricity consumption side carbon emission coefficient allocation method according to claim 7, wherein the establishing method of the carbon emission coefficient analysis formula is as follows:
to be used forFor the objective function, in->Establishing a carbon emission coefficient allocation model for constraint conditions; wherein u is l Carbon emission coefficient for load node l; τ l Carbon emission influencing factors for load nodes l; pi s Is the total coefficient of carbon quota; l is the total number of load nodes;
and deducing the carbon emission coefficient allocation model to obtain a carbon emission coefficient analysis formula.
10. A power side carbon emission factor apportionment system, the system comprising:
the power prediction module is used for calculating the injection power of each node in each period according to the power generation side quotation of the target area and the historical power of each node in each period; the nodes comprise generator nodes and load nodes;
the load side carbon flow injection matrix determining module is used for calculating a load side carbon flow injection matrix according to the injection power of all the nodes in each period; the load side carbon stream injection matrix includes carbon emission streams for each of the load nodes at each time period;
the carbon emission parameter determining module is used for calculating the total carbon quota coefficient of the target area and the carbon emission influence factor of each load node according to the load side carbon flow injection matrix; the carbon emission influence factor characterizes the influence degree of the electricity consumption of a user on the carbon emission of the electric power system;
and the carbon transaction coefficient allocation module is used for calculating the carbon emission coefficient of each load node according to the total carbon quota coefficient and the carbon emission influence factor of the load node.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202310793108.XA CN116797245A (en) | 2023-06-30 | 2023-06-30 | Power utilization side carbon emission coefficient allocation method and system |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202310793108.XA CN116797245A (en) | 2023-06-30 | 2023-06-30 | Power utilization side carbon emission coefficient allocation method and system |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN116797245A true CN116797245A (en) | 2023-09-22 |
Family
ID=88034592
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202310793108.XA Pending CN116797245A (en) | 2023-06-30 | 2023-06-30 | Power utilization side carbon emission coefficient allocation method and system |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN116797245A (en) |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN117239844A (en) * | 2023-11-15 | 2023-12-15 | 广东电网有限责任公司广州供电局 | Power system scheduling method, device and storage medium based on carbon emission responsibility |
-
2023
- 2023-06-30 CN CN202310793108.XA patent/CN116797245A/en active Pending
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN117239844A (en) * | 2023-11-15 | 2023-12-15 | 广东电网有限责任公司广州供电局 | Power system scheduling method, device and storage medium based on carbon emission responsibility |
| CN117239844B (en) * | 2023-11-15 | 2024-04-05 | 广东电网有限责任公司广州供电局 | Power system scheduling method, device and storage medium based on carbon emission responsibility |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN105844367A (en) | Sensitivity factor based calculating method for power grid node marginal pricing | |
| CN106230028A (en) | A kind of Multipurpose Optimal Method of wind-powered electricity generation water-storage association system | |
| CN104537428B (en) | One kind meter and the probabilistic economical operation appraisal procedure of wind power integration | |
| Dias et al. | Parallel computing applied to the stochastic dynamic programming for long term operation planning of hydrothermal power systems | |
| CN110490363A (en) | More days Unit Combination optimization methods of one kind and system | |
| CN116797245A (en) | Power utilization side carbon emission coefficient allocation method and system | |
| Muthirayan et al. | Mechanism design for self-reporting baselines in demand response | |
| Rong et al. | An effective heuristic for combined heat-and-power production planning with power ramp constraints | |
| CN105956247A (en) | Open, fair and impartial scheduling model building method based on Gini coefficient | |
| CN112001576A (en) | Accounting method for electric power consumption of renewable energy source | |
| Sharma et al. | Metrics-based assessment of sustainability in demand response | |
| Teng et al. | Availability estimation of wind power forecasting and optimization of day-ahead unit commitment | |
| CN105825348A (en) | Method for quantitatively analyzing influence exerted by bilateral trade on power grid power transmission cost | |
| CN105226707A (en) | A kind of methodology based on Shapley value wind-electricity integration system fixed cost of power transmission | |
| Chen et al. | Multi-objective robust optimal bidding strategy for a data center operator based on bi-level optimization | |
| CN112101607A (en) | Active power distribution network rolling optimization scheduling method considering demand response time effect | |
| CN108767855A (en) | A kind of electric system random production analog method that sequential persistently mixes | |
| Reinschmidt et al. | Construction cash flow management system | |
| CN111242438A (en) | Evaluation method and system for flexibility adjustment capability of power generation and utilization resources of self-contained power plant | |
| CN102609824B (en) | The acquisition methods of critical point power trade capacity index under polynary electrification structure | |
| Brodsky et al. | A Decision-guided advisor to maximize ROI in local generation & utility contracts | |
| Celik et al. | A comparison of three parallel processing methods for a resource allocation problem in the smart grid | |
| Chen et al. | Overview of transmission expansion planning in the market environment | |
| CN109740178B (en) | Multi-tenant data center energy efficiency optimization method and system and combined modeling method | |
| de Oca et al. | Demand response program for supercomputing and datacenters providing ancillary services in the electricity market |
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 |