CN116596330A - Accurate carbon emission flow acquisition method considering loss correction of transmission line - Google Patents

Accurate carbon emission flow acquisition method considering loss correction of transmission line Download PDF

Info

Publication number
CN116596330A
CN116596330A CN202310455923.5A CN202310455923A CN116596330A CN 116596330 A CN116596330 A CN 116596330A CN 202310455923 A CN202310455923 A CN 202310455923A CN 116596330 A CN116596330 A CN 116596330A
Authority
CN
China
Prior art keywords
loss
line
node
network
power
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
Application number
CN202310455923.5A
Other languages
Chinese (zh)
Inventor
申冉
唐爱红
李黄强
刘子伟
张宇
谢琼瑶
金潇
王庆铭
游欣宇
余文晗
毛鲁洁
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Wuhan University of Technology WUT
Yichang Power Supply Co of State Grid Hubei Electric Power Co Ltd
Original Assignee
Wuhan University of Technology WUT
Yichang Power Supply Co of State Grid Hubei Electric Power Co Ltd
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Wuhan University of Technology WUT, Yichang Power Supply Co of State Grid Hubei Electric Power Co Ltd filed Critical Wuhan University of Technology WUT
Priority to CN202310455923.5A priority Critical patent/CN116596330A/en
Publication of CN116596330A publication Critical patent/CN116596330A/en
Pending legal-status Critical Current

Links

Classifications

    • HELECTRICITY
    • H02GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02JCIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
    • H02J3/00Circuit arrangements for ac mains or ac distribution networks
    • H02J3/04Circuit arrangements for ac mains or ac distribution networks for connecting networks of the same frequency but supplied from different sources
    • H02J3/06Controlling transfer of power between connected networks; Controlling sharing of load between connected networks
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06QINFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
    • G06Q10/00Administration; Management
    • G06Q10/06Resources, workflows, human or project management; Enterprise or organisation planning; Enterprise or organisation modelling
    • G06Q10/063Operations research, analysis or management
    • G06Q10/0637Strategic management or analysis, e.g. setting a goal or target of an organisation; Planning actions based on goals; Analysis or evaluation of effectiveness of goals
    • G06Q10/06375Prediction of business process outcome or impact based on a proposed change
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06QINFORMATION 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/00Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
    • G06Q50/06Energy or water supply
    • HELECTRICITY
    • H02GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02JCIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
    • H02J3/00Circuit arrangements for ac mains or ac distribution networks
    • HELECTRICITY
    • H02GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02JCIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
    • H02J2203/00Indexing scheme relating to details of circuit arrangements for AC mains or AC distribution networks
    • H02J2203/10Power transmission or distribution systems management focussing at grid-level, e.g. load flow analysis, node profile computation, meshed network optimisation, active network management or spinning reserve management
    • HELECTRICITY
    • H02GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02JCIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
    • H02J2203/00Indexing scheme relating to details of circuit arrangements for AC mains or AC distribution networks
    • H02J2203/20Simulating, e g planning, reliability check, modelling or computer assisted design [CAD]
    • YGENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02PCLIMATE CHANGE MITIGATION TECHNOLOGIES IN THE PRODUCTION OR PROCESSING OF GOODS
    • Y02P90/00Enabling technologies with a potential contribution to greenhouse gas [GHG] emissions mitigation
    • Y02P90/80Management or planning
    • Y02P90/84Greenhouse gas [GHG] management systems

Landscapes

  • Business, Economics & Management (AREA)
  • Engineering & Computer Science (AREA)
  • Human Resources & Organizations (AREA)
  • Economics (AREA)
  • Strategic Management (AREA)
  • Health & Medical Sciences (AREA)
  • Physics & Mathematics (AREA)
  • Entrepreneurship & Innovation (AREA)
  • Power Engineering (AREA)
  • Educational Administration (AREA)
  • Marketing (AREA)
  • Theoretical Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Tourism & Hospitality (AREA)
  • General Business, Economics & Management (AREA)
  • Development Economics (AREA)
  • Quality & Reliability (AREA)
  • Operations Research (AREA)
  • Game Theory and Decision Science (AREA)
  • Public Health (AREA)
  • Water Supply & Treatment (AREA)
  • General Health & Medical Sciences (AREA)
  • Primary Health Care (AREA)
  • Supply And Distribution Of Alternating Current (AREA)

Abstract

A method for accurately acquiring carbon emission flow considering loss correction of a power transmission line comprises the following steps: step 1: carrying out correction of transmission line loss; step 2: performing lossless network conversion for correcting the loss of the transmission line; step 3: the obtained line loss correction amounts of all the lines are equivalent to the load of the line non-node; step 4: and acquiring a carbon emission stream for transmission line loss correction. The invention aims to solve the technical problem that the carbon flow rate of the power generation side and the carbon flow rate of the load side obtained by the existing carbon flow calculation method have a certain difference, so that the monitoring effect of the carbon emission of a power grid can be affected.

Description

Accurate carbon emission flow acquisition method considering loss correction of transmission line
Technical Field
The invention belongs to the technical field of power systems, in particular to the fields of accurate calculation, carbon reduction and carbon reduction of carbon emission, and particularly relates to an accurate calculation method of carbon emission flow considering power transmission line loss correction.
Background
According to the existing carbon flow obtaining method, the calculated power generation side carbon flow rate should be equal to the load side carbon flow rate, but in practice, there is a certain difference between the power generation side carbon flow rate and the load side carbon flow rate calculated by the existing carbon flow calculating method. The actual situation is different from the theoretical calculation result, and the research difficulty of the problems about line loss carbon rate allocation and the like is improved.
Patent document with the authority of publication number CN115907307A discloses an online analysis method of carbon emission flow of a power system facing real-time data interaction of a power grid, and based on a carbon emission flow calculation method, real-time data in the power grid are analyzed and calculated to obtain the real-time carbon emission flow of the power system. Although the method can obtain the real-time carbon emission flow of the electric power system, the calculated line loss carbon rate still has a certain error with the actual line loss carbon flow.
When the line loss carbon rate calculation is carried out, the theoretical line loss carbon rate is not only included, but also the actual line carbon rate generated by factors such as temperature compensation loss, corona loss and the like are included.
Therefore, the applicant proposes a method for accurately calculating and acquiring a carbon emission flow taking account of transmission line loss correction.
Disclosure of Invention
The invention aims to solve the technical problem that the carbon flow rate of the power generation side and the carbon flow rate of the load side obtained by the existing carbon flow calculation method have a certain difference, so that the monitoring effect of the carbon emission of a power grid can be affected.
In order to solve the technical problems, the invention adopts the following technical scheme:
a method for accurately acquiring carbon emission flow considering loss correction of a power transmission line comprises the following steps:
step 1: carrying out correction of transmission line loss;
step 2: performing lossless network conversion for correcting the loss of the transmission line;
step 3: the obtained line loss correction amounts of all the lines are equivalent to the load of the line non-node;
step 4: and acquiring a carbon emission stream for transmission line loss correction.
In step 1, correction includes overhead line temperature compensation, corona loss compensation and overhead ground wire power loss compensation;
step 1.1: calculating the temperature compensation of the overhead line:
ΔS′ L =k w ×ΔS L (1)
wherein, deltaS L Line loss before temperature compensation; ΔS' L Line loss after temperature compensation is carried out on the basis of load flow calculation; k (k) w Is the temperature rise coefficient;
step 1.2: calculating corona loss of overhead line
The calculation method of the corona loss is shown as a formula (2).
In the method, in the process of the invention,to account for line losses after corona losses.
Step 1.3: the electric energy loss of the overhead ground wire is calculated as follows:
(1) Double overhead ground wire
Loss Δs generated by induced current between two ground wires H Calculated from equation (3).
Wherein L is the length of the overhead ground wire; d, d 1A 、d 1C And d 12 The space between the overhead ground wire, the lead and the two-phase line is respectively; t is the line loss calculation period; i jf Calculating root mean square current in a period for the line loss; r is R i The resistance of the overhead ground wire is 20 ℃; r is the radius of the overhead ground wire; x is x i The unit inductive reactance is the double-overhead ground wire;
loss Δs due to induced current between surrounding ground and earth D
Wherein f is the current frequency; d (D) 0 The equivalent depth of the induced current is;
so that the double-frame ground wire and the two sides are all grounded for loss,
ΔS 2 =ΔS H +ΔS D (5)
(2) Single overhead ground wire
Single overhead ground wire energy consumption (or double overhead ground wire single side grounding energy consumption) delta S 1
Wherein S is 1 Is complex induced electromotive force.
In step 2, the network loss obtained by calculating the power flow is required to be distributed to the load side.
For any node in the net rack, the inflow power S of the node i is obtained i The following is shown:
wherein i is - A downstream node set for node i; s is S ij Power flowing out for node i on line ij; s is S Li If the node i is not loaded, the variable is 0.
The invention adopts the countercurrent tracking method to solve the network loss of each generator caused by the power injected into the net rack, thus S in the above formula can be obtained ij Power |s to become injection node j ji Net loss S on i and line ij ij,loss And, if the sum is the same, the formula (7) can be:
wherein S is i-,loss Is the sum of the network losses of all downstream node branches connected with the node i. Decreasing the homogeneous phase, the above formula can be rewritten as:
the above formula is converted into a matrix form, and can be expressed as:
A d S=S L +S Σ-,loss (10)
wherein S is Σ-,loss A network loss set on a branch consisting of downstream nodes of all nodes in the network frame; s is a vector consisting of injection power of all nodes in the net rack; s is S L The node load power vector; a is that d The elements of the downstream distribution matrix are shown as follows:
let the injection power on the generator bus i in the network be S Gi The parameter may be rewritten as shown in the following equation:
in the formula e i T ∈R n Column vectors with the ith variable being 1 and the remaining variables being 0, where n represents the number of nodes in the grid. From equation (12), the net loss caused by the generator node i is S Gi,loss
Thus, the output S 'of each generator after conversion to a lossless network' Gi Can be calculated from equation (14):
S′ Gi =S Gi -S Gi,loss (14)
in summary, the network can be converted into a substantially lossless network.
In order to not increase the calculation dimension and reduce the number of nodes in the grid, the method equates the calculated correction quantity of the power transmission line loss to the load at the end node of the line, as shown in fig. 1, so that a lossless network for calculating and correcting the line loss can be obtained.
In the step (3) of the process,
step 3.1: the direction of the current of each line is determined as follows:
t bus =[t 1 ,t 2 ,t 3 ,…,t i ,…,t nl ] (15)
wherein t is bus Injecting a number vector of the node for the line tide; n is n l The number of the lines is the number of the lines; t is t i The number of the power flow injection node on the ith line is given.
Step 3.2: the line loss correction quantity of each line calculated in the step 1 is equivalent to t bus Node load at.
In step 4, when obtaining a carbon emission stream for transmission line loss correction;
step 4.1: firstly, a direct current tide distribution matrix, a unit injection distribution matrix and a generator unit carbon emission intensity vector are calculated according to tide results. Wherein P is B Is a direct current power flow distribution matrix of the network, and(n l the number of branches in the net rack); p (P) G Injecting a distribution matrix for units of the network, and +.>(n g The number of the generators in the net rack is n b Is the number of nodes in the net rack); e (E) G Is the carbon emission intensity of the generator set of the network, and +.>
Step 4.2: according to the variables, let P Z =[P B P G ] T The active power flux P can be obtained N The method comprises the following steps:
in the method, in the process of the invention,is n b +n g The column vector of the rank, and all elements are 1.
Step 4.3: calculating the node carbon potential E of the network N Branch carbon flow rate distribution matrix R B
And (3) obtaining the carbon emission flow considering the loss correction of the transmission line.
Compared with the prior art, the invention has the following technical effects:
according to the accurate acquisition method for the carbon emission flow considering the power transmission line loss correction, which is provided by the invention, the actual line loss of each line in the grid can be accurately calculated, the network is converted into a lossless network through a tide tracking method, then the line loss correction value of each power transmission line is equivalent to a load, and finally the carbon emission rate of each node in the actual operation of the power grid is calculated through the carbon emission flow, so that a foundation is laid for the accurate calculation of the regional power grid carbon emission flow, and the actual carbon emission rate of each node in the obtained grid can be ensured.
Drawings
The invention is further illustrated by the following examples in conjunction with the accompanying drawings:
FIG. 1 is a flow chart of the present invention;
FIG. 2 is a schematic diagram of a lossless network for calculating and correcting line loss according to the present invention.
Detailed Description
A method for accurately acquiring carbon emission flow considering loss correction of a power transmission line comprises the following steps:
step 1: carrying out correction of transmission line loss;
step 2: performing lossless network conversion for correcting the loss of the transmission line;
step 3: the obtained line loss correction amounts of all the lines are equivalent to the load of the line non-node; step 4: and acquiring a carbon emission stream for transmission line loss correction.
In step 1, correction includes overhead line temperature compensation, corona loss compensation and overhead ground wire power loss compensation;
step 1.1: calculating the temperature compensation of the overhead line:
ΔS′ L =k w ×ΔS L (1)
wherein, deltaS L Line loss before temperature compensation; deltaS′ L Line loss after temperature compensation is carried out on the basis of load flow calculation; k (k) w Is the temperature rise coefficient;
step 1.2: calculating corona loss of overhead line
The calculation method of the corona loss is shown as a formula (2).
In the method, in the process of the invention,to account for line losses after corona losses.
Step 1.3: the electric energy loss of the overhead ground wire is calculated as follows:
(1) Double overhead ground wire
Loss Δs generated by induced current between two ground wires H Calculated from equation (3).
Wherein L is the length of the overhead ground wire; d, d 1A 、d 1C And d 12 The space between the overhead ground wire, the lead and the two-phase line is respectively; t is the line loss calculation period; i jf Calculating root mean square current in a period for the line loss; r is R i The resistance of the overhead ground wire is 20 ℃; r is the radius of the overhead ground wire; x is x i The unit inductive reactance is the double-overhead ground wire;
loss Δs due to induced current between surrounding ground and earth D
Wherein f is the current frequency; d (D) 0 The equivalent depth of the induced current is;
so that the double-frame ground wire and the two sides are all grounded for loss,
ΔS 2 =ΔS H +ΔS D (5)
(2) Single overhead ground wire
Single overhead ground wire energy consumption (or double overhead ground wire single side grounding energy consumption) delta S 1
Wherein S is 1 Is complex induced electromotive force.
In step 2, the network loss obtained by calculating the power flow is required to be distributed to the load side.
For any node in the net rack, the inflow power S of the node i is obtained i The following is shown:
wherein i is - A downstream node set for node i; s is S ij Power flowing out for node i on line ij; s is S Li If the node i is not loaded, the variable is 0.
The invention adopts the countercurrent tracking method to solve the network loss of each generator caused by the power injected into the net rack, thus S in the above formula can be obtained ij Power |s to become injection node j ji Net loss S on i and line ij ij,loss And, if the sum is the same, the formula (7) can be:
wherein S is i-,loss Is the sum of the network losses of all downstream node branches connected with the node i. Decreasing the homogeneous phase, the above formula can be rewritten as:
the above formula is converted into a matrix form, and can be expressed as:
A d S=S L +S Σ-,loss (10)
wherein S is Σ-,loss A network loss set on a branch consisting of downstream nodes of all nodes in the network frame; s is a vector consisting of injection power of all nodes in the net rack; s is S L The node load power vector; a is that d The elements of the downstream distribution matrix are shown as follows:
let the injection power on the generator bus i in the network be S Gi The parameter may be rewritten as shown in the following equation:
in the method, in the process of the invention,column vectors for the ith variable of 1 and the remaining variables of 0, where n represents the number of nodes in the grid. From equation (12), the net loss caused by the generator node i is S Gi,loss
Thus, the output S 'of each generator after conversion to a lossless network' Gi Can be calculated from equation (14):
S′ Gi =S Gi -S Gi,loss (14)
in summary, the network can be converted into a substantially lossless network.
In order to not increase the calculation dimension and reduce the number of nodes in the grid, the method equates the calculated correction quantity of the power transmission line loss to the load at the end node of the line, as shown in fig. 1, so that a lossless network for calculating and correcting the line loss can be obtained.
In step 3:
step 3.1: the direction of the current of each line is determined as follows:
t bus =[t 1 ,t 2 ,t 3 ,…,t i ,…,t nl ] (15)
wherein t is bus Injecting a number vector of the node for the line tide; n is n l The number of the lines is the number of the lines; t is t i The number of the power flow injection node on the ith line is given.
Step 3.2: the line loss correction quantity of each line calculated in the step 1 is equivalent to t bus Node load at.
In step 4, when obtaining a carbon emission stream for transmission line loss correction;
step 4.1: firstly, a direct current tide distribution matrix, a unit injection distribution matrix and a generator unit carbon emission intensity vector are calculated according to tide results. Wherein P is B Is a direct current power flow distribution matrix of the network, and(n l the number of branches in the net rack); p (P) G Injecting a distribution matrix for units of the network, and +.>(n g The number of the generators in the net rack is n b Is the number of nodes in the net rack); e (E) G Is the carbon emission intensity of the generator set of the network, and +.>
Step 4.2: according to the variables, let P Z =[P B P G ] T The active power flux P can be obtained N The method comprises the following steps:
in the method, in the process of the invention,is n b +n g The column vector of the rank, and all elements are 1.
Step 4.3: calculating the node carbon potential E of the network N Branch carbon flow rate distribution matrix R B
And (3) obtaining the carbon emission flow considering the loss correction of the transmission line.

Claims (5)

1. The accurate acquisition method of the carbon emission flow considering the loss correction of the power transmission line is characterized by comprising the following steps of:
step 1: carrying out correction of transmission line loss;
step 2: performing lossless network conversion for correcting the loss of the transmission line;
step 3: the obtained line loss correction amounts of all the lines are equivalent to the load of the line non-node;
step 4: and acquiring a carbon emission stream for transmission line loss correction.
2. The method of claim 1, wherein in step 1, the correction includes overhead line temperature compensation, corona loss compensation, and overhead ground wire power loss compensation;
step 1.1: calculating the temperature compensation of the overhead line:
ΔS′ L =k w ×ΔS L (1)
wherein, deltaS L Line loss before temperature compensation; ΔS' L Line loss after temperature compensation is carried out on the basis of load flow calculation; k (k) w Is the temperature rise coefficient;
step 1.2: calculating corona loss of overhead line
The calculation method of the corona loss is shown in the formula (2);
in the method, in the process of the invention,to account for line losses after corona losses;
step 1.3: the electric energy loss of the overhead ground wire is calculated as follows:
(1) Double overhead ground wire
Loss Δs generated by induced current between two ground wires H Calculated from formula (3);
wherein L is the length of the overhead ground wire; d, d 1A 、d 1C And d 12 The space between the overhead ground wire, the lead and the two-phase line is respectively; t is the line loss calculation period; i jf Calculating root mean square current in a period for the line loss; r is R i The resistance of the overhead ground wire is 20 ℃; r is the radius of the overhead ground wire; x is x i The unit inductive reactance is the double-overhead ground wire;
loss Δs due to induced current between surrounding ground and earth D
Wherein f is the current frequency; d (D) 0 The equivalent depth of the induced current is;
so that the double-frame ground wire and the two sides are all grounded for loss,
ΔS 2 =ΔS H +ΔS D (5)
(2) Single overhead ground wire
Single overhead ground wire energy consumption (or double overhead ground wire single side grounding energy consumption) delta S 1
Wherein S is 1 Is complex induced electromotive force.
3. The method according to claim 1, characterized in that in step 2, the power flow calculation is performed to obtain a net loss and the net loss is distributed to the load side;
for any node in the net rack, the inflow power S of the node i is obtained i The following is shown:
wherein i is - A downstream node set for node i; s is S ij Power flowing out for node i on line ij; s is S Li For the load on node i, if node i is not loaded, the variable is 0;
the invention adopts the countercurrent tracking method to solve the network loss of each generator caused by the power injected into the net rack, thus S in the above formula can be obtained ij Power |s to become injection node j ji Net loss S on i and line ij ij,loss And, if the sum is the same, the formula (7) can be:
in the method, in the process of the invention,the sum of the network losses of all downstream node branches connected with the node i is obtained; decreasing the homogeneous phase, the above formula can be rewritten as:
the above formula is converted into a matrix form, and can be expressed as:
in the method, in the process of the invention,a network loss set on a branch consisting of downstream nodes of all nodes in the network frame; s is a vector consisting of injection power of all nodes in the net rack; s is S L The node load power vector; a is that d The elements of the downstream distribution matrix are shown as follows:
let the injection power on the generator bus i in the network be S Gi The parameter may be rewritten as shown in the following equation:
in the method, in the process of the invention,column vectors with the ith variable being 1 and the rest variables being 0, wherein n represents the number of nodes in the net rack;
from equation (12), the net loss caused by the generator node i is S Gi,loss
Thus, the output S 'of each generator after conversion to a lossless network' Gi Can be calculated from equation (14):
S′ Gi =S Gi -S Gi,loss (14)
to sum up, the network can be converted into a basically lossless network;
in order to not increase the calculation dimension and reduce the number of nodes in the grid, the method equates the calculated correction quantity of the power transmission line loss to the load at the end node of the line, as shown in fig. 1, so that a lossless network for calculating and correcting the line loss can be obtained.
4. The method according to claim 1, wherein, in step 3,
step 3.1: the direction of the current of each line is determined as follows:
t bus =[t 1 ,t 2 ,t 3 ,…,t i ,…,t nl ] (15)
wherein t is bus Injecting a number vector of the node for the line tide; n is n l The number of the lines is the number of the lines; t is t i Numbering the power flow injection nodes on the ith line;
step 3.2: the line loss correction quantity of each line calculated in the step 1 is equivalent to t bus Node load at.
5. The method according to claim 1, characterized in that in step 4, when obtaining a transmission line loss corrected carbon emission stream;
step 4.1: firstly, calculating a direct current tide distribution matrix, a unit injection distribution matrix and a carbon emission intensity vector of a generator unit according to tide results; wherein P is B Is a direct current power flow distribution matrix of the network, and(n l the number of branches in the net rack); p (P) G Injecting a distribution matrix for units of the network, and +.>(n g The number of the generators in the net rack is n b Is the number of nodes in the net rack); e (E) G Is the carbon emission intensity of the generator set of the network, and +.>
Step 4.2: according to the variables, let P Z =[P B P G ] T The active power flux P can be obtained N The method comprises the following steps:
in the method, in the process of the invention,is n b +n g A row vector of the order, and all elements are 1;
step 4.3: calculating the node carbon potential E of the network N Branch carbon flow rate distribution matrix R B
And (3) obtaining the carbon emission flow considering the loss correction of the transmission line.
CN202310455923.5A 2023-04-25 2023-04-25 Accurate carbon emission flow acquisition method considering loss correction of transmission line Pending CN116596330A (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202310455923.5A CN116596330A (en) 2023-04-25 2023-04-25 Accurate carbon emission flow acquisition method considering loss correction of transmission line

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202310455923.5A CN116596330A (en) 2023-04-25 2023-04-25 Accurate carbon emission flow acquisition method considering loss correction of transmission line

Publications (1)

Publication Number Publication Date
CN116596330A true CN116596330A (en) 2023-08-15

Family

ID=87607052

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202310455923.5A Pending CN116596330A (en) 2023-04-25 2023-04-25 Accurate carbon emission flow acquisition method considering loss correction of transmission line

Country Status (1)

Country Link
CN (1) CN116596330A (en)

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20180067089A1 (en) * 2016-09-06 2018-03-08 Tsinghua University System for measuring carbon emission in power system
CN113886767A (en) * 2021-08-30 2022-01-04 国网江苏省电力有限公司 Real-time calculation method for carbon emission flow of power system
CN115795889A (en) * 2022-12-06 2023-03-14 国家电网有限公司华中分部 Power transmission line loss optimization calculation method based on multi-parameter correction
CN115906347A (en) * 2022-10-25 2023-04-04 国网天津市电力公司电力科学研究院 Power system carbon flow tracking method based on carbon emission flow theory

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20180067089A1 (en) * 2016-09-06 2018-03-08 Tsinghua University System for measuring carbon emission in power system
CN113886767A (en) * 2021-08-30 2022-01-04 国网江苏省电力有限公司 Real-time calculation method for carbon emission flow of power system
CN115906347A (en) * 2022-10-25 2023-04-04 国网天津市电力公司电力科学研究院 Power system carbon flow tracking method based on carbon emission flow theory
CN115795889A (en) * 2022-12-06 2023-03-14 国家电网有限公司华中分部 Power transmission line loss optimization calculation method based on multi-parameter correction

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
冯欣等: "考虑网络损耗的碳排放流理论改进与完善", 电力自动化设备, vol. 36, no. 05, 31 May 2016 (2016-05-31), pages 81 - 86 *

Similar Documents

Publication Publication Date Title
CN108270248B (en) Active safety correction control method based on quasi-steady-state sensitivity analysis
CN107204617B (en) Linear programming based interval power flow calculation method in rectangular coordinate form
CN104636829B (en) A kind of decoupling algorithm improving electric system temperature optimal load flow computational efficiency
CN108023364A (en) Power distribution network distributed generation resource maximum access capability computational methods based on convex difference planning
CN110489806B (en) Electromagnetic transient modeling and calculating method for power grid comprising multi-voltage-source converter
CN107689772A (en) The evaluation method of the submatrix parallel connection mismatch of photovoltaic plant
CN103969552A (en) Harmonic source positioning and analyzing method for distributed power generation system
EP2628224B1 (en) Method and apparatus for status signal generation
CN111900718B (en) Active power distribution network dynamic simulation method based on multi-stage optimization catch-up variational iteration method
CN113777440A (en) Method for positioning single-phase earth fault of current collection circuit of photovoltaic power station
CN104242308A (en) Multisource system partitioning method considering voltage stabilization
CN107482639B (en) dynamic reactive power equipment reactive power optimization method suitable for wind power system in power grid
CN112434408A (en) Photovoltaic power station cable length optimization method, cable model selection method and device
CN107732904B (en) Active static security domain construction method considering single TCSC equivalent reactance parameter
CN116596330A (en) Accurate carbon emission flow acquisition method considering loss correction of transmission line
CN106991229B (en) Wind power plant equivalent modeling method for complex topology
CN109149584A (en) Improve IPFC power injection model constringent method in Load flow calculation
CN112836364A (en) Unified calculation method for urban rail transit stray current
CN107271941A (en) A kind of zero-sequence mutual inductance impedance computation method on multiple-loop line transmission line of electricity
CN101841154A (en) Voltage stability margin real-time evaluation and optimum control method after grid major failure
CN110011313A (en) A kind of flexible direct current power grid load flow calculation method and system
CN108183487A (en) A kind of distribution power flow rapid analysis method based on linear mathematical model
Ćetković et al. Analysis of justification for using capacitor banks in distribution networks with low power demand
CN111239511B (en) Online measurement method for generalized short-circuit ratio of multi-feed-in new energy grid-connected system
CN110543739A (en) Circuit simulation model of overhead power transmission line

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