CN116502961A - Power grid carbon flow calculation method and device considering system loss and wind power uncertainty - Google Patents
Power grid carbon flow calculation method and device considering system loss and wind power uncertainty Download PDFInfo
- Publication number
- CN116502961A CN116502961A CN202310770811.9A CN202310770811A CN116502961A CN 116502961 A CN116502961 A CN 116502961A CN 202310770811 A CN202310770811 A CN 202310770811A CN 116502961 A CN116502961 A CN 116502961A
- Authority
- CN
- China
- Prior art keywords
- power
- wind
- carbon flow
- wind speed
- 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 282
- 229910052799 carbon Inorganic materials 0.000 title claims abstract description 282
- 238000004364 calculation method Methods 0.000 title claims abstract description 113
- 238000000034 method Methods 0.000 claims abstract description 154
- 238000002347 injection Methods 0.000 claims abstract description 99
- 239000007924 injection Substances 0.000 claims abstract description 99
- 238000004458 analytical method Methods 0.000 claims abstract description 28
- 238000005206 flow analysis Methods 0.000 claims abstract description 21
- 239000011159 matrix material Substances 0.000 claims description 114
- 238000009826 distribution Methods 0.000 claims description 92
- 230000008569 process Effects 0.000 claims description 60
- 230000006870 function Effects 0.000 claims description 23
- 230000004907 flux Effects 0.000 claims description 21
- 230000000694 effects Effects 0.000 claims description 20
- 238000005315 distribution function Methods 0.000 claims description 15
- 230000001186 cumulative effect Effects 0.000 claims description 13
- 238000011144 upstream manufacturing Methods 0.000 claims description 11
- 238000009827 uniform distribution Methods 0.000 claims description 10
- 238000003860 storage Methods 0.000 claims description 8
- 238000004590 computer program Methods 0.000 claims description 6
- 238000000605 extraction Methods 0.000 claims description 4
- 238000004519 manufacturing process Methods 0.000 claims description 4
- 238000010276 construction Methods 0.000 claims description 3
- 239000000243 solution Substances 0.000 description 7
- 238000010205 computational analysis Methods 0.000 description 3
- 230000008878 coupling Effects 0.000 description 3
- 238000010168 coupling process Methods 0.000 description 3
- 238000005859 coupling reaction Methods 0.000 description 3
- 238000004401 flow injection analysis Methods 0.000 description 3
- 238000010248 power generation Methods 0.000 description 3
- 230000009286 beneficial effect Effects 0.000 description 2
- 238000006243 chemical reaction Methods 0.000 description 2
- 238000004891 communication Methods 0.000 description 2
- 238000005314 correlation function Methods 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 239000004973 liquid crystal related substance Substances 0.000 description 2
- 238000011160 research Methods 0.000 description 2
- 238000011166 aliquoting Methods 0.000 description 1
- 230000005540 biological transmission Effects 0.000 description 1
- 238000004422 calculation algorithm Methods 0.000 description 1
- 238000005265 energy consumption Methods 0.000 description 1
- 230000007246 mechanism Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 238000005070 sampling Methods 0.000 description 1
- 238000007619 statistical method Methods 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- 238000013076 uncertainty analysis Methods 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
- G06Q10/00—Administration; Management
- G06Q10/06—Resources, workflows, human or project management; Enterprise or organisation planning; Enterprise or organisation modelling
- G06Q10/063—Operations research, analysis or management
- G06Q10/0639—Performance analysis of employees; Performance analysis of enterprise or organisation operations
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/11—Complex mathematical operations for solving equations, e.g. nonlinear equations, general mathematical optimization problems
- G06F17/12—Simultaneous equations, e.g. systems of linear equations
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/16—Matrix or vector computation, e.g. matrix-matrix or matrix-vector multiplication, matrix factorization
-
- 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—Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
- G06Q50/06—Energy or water supply
-
- H—ELECTRICITY
- H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
- H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
- H02J3/00—Circuit arrangements for ac mains or ac distribution networks
- H02J3/38—Arrangements for parallely feeding a single network by two or more generators, converters or transformers
- H02J3/381—Dispersed generators
-
- H—ELECTRICITY
- H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
- H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
- H02J2203/00—Indexing scheme relating to details of circuit arrangements for AC mains or AC distribution networks
- H02J2203/10—Power 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
-
- H—ELECTRICITY
- H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
- H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
- H02J2300/00—Systems for supplying or distributing electric power characterised by decentralized, dispersed, or local generation
- H02J2300/20—The dispersed energy generation being of renewable origin
- H02J2300/28—The renewable source being wind energy
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Business, Economics & Management (AREA)
- Mathematical Physics (AREA)
- Theoretical Computer Science (AREA)
- Mathematical Optimization (AREA)
- Pure & Applied Mathematics (AREA)
- Mathematical Analysis (AREA)
- Computational Mathematics (AREA)
- Human Resources & Organizations (AREA)
- Economics (AREA)
- Data Mining & Analysis (AREA)
- Strategic Management (AREA)
- General Business, Economics & Management (AREA)
- Software Systems (AREA)
- Entrepreneurship & Innovation (AREA)
- Health & Medical Sciences (AREA)
- General Engineering & Computer Science (AREA)
- Databases & Information Systems (AREA)
- Marketing (AREA)
- Algebra (AREA)
- Tourism & Hospitality (AREA)
- Operations Research (AREA)
- Educational Administration (AREA)
- Development Economics (AREA)
- Quality & Reliability (AREA)
- Primary Health Care (AREA)
- General Health & Medical Sciences (AREA)
- Water Supply & Treatment (AREA)
- Game Theory and Decision Science (AREA)
- Public Health (AREA)
- Computing Systems (AREA)
- Power Engineering (AREA)
- Wind Motors (AREA)
Abstract
The invention relates to the technical field of analysis of carbon emission flows of power systems, and discloses a power grid carbon flow calculation method and device considering system loss and wind power uncertainty. The method comprises the steps of constructing a random wind speed sample set aiming at a target power system; based on each group of random wind speed samples in the random wind speed sample set, carrying out deterministic power flow calculation accounting for active loss, and determining a functional expression between wind power injection power considering wind speed uncertainty and the total carbon flow rate of the target power system; calculating related carbon flow indexes according to the obtained deterministic trend calculation result, and carrying out carbon flow analysis based on the related carbon flow indexes; the carbon flow analysis includes calculating contributions of balancing units to active power flows of portions of the target power system and calculating carbon flow rate influence factors of balancing units to nodes and branches. The method and the system realize the calculation of the carbon flow of the power grid by considering the system loss and wind power uncertainty.
Description
Technical Field
The invention relates to the technical field of analysis of carbon emission flows of power systems, in particular to a method and a device for calculating a power grid carbon flow by considering system loss and wind power uncertainty.
Background
Researches show that the wind power can be connected to replace the traditional energy source to reduce emission, and the tide distribution in the system can be effectively changed. Because the carbon emission flow of the system is closely related to tide characteristics, the influence and analysis of the random and intermittent wind power injection power on the carbon emission flow of the system are considered, and the method has important significance for quantitatively evaluating the contribution of wind power generation to the low carbon of the system.
Currently, the carbon emission calculation of the electric power system mainly comprises a macroscopic statistical method and a carbon flow analysis method. The macro statistics method can start from macro data and count according to total energy consumption in a period of time, but the calculation result of the method is too extensive, and cannot track the concrete flow direction of carbon emission, so that the division and the identification of carbon emission responsibility are not facilitated. The carbon flow analysis method is a carbon flow tracking method based on power distribution, can clearly reveal the distribution characteristics and transmission consumption mechanisms of carbon flows in an electric power network, realizes accurate tracking and tracing of specific flow directions of carbon emission, and has important significance in various carbon emission analysis and statistics of an electric power system.
Based on the carbon flow analysis method, ma Rui et al in journal entitled power system automation, present a method for calculating and analyzing the influence of wind power injection power considering wind speed uncertainty on the carbon emission flow of a power system. Firstly, acquiring a system tide based on a random direct current tide of uncertainty of wind power injection power, and acquiring a correlation function of wind speed for determining wind power output and total carbon injection flow rate of the system through a calculation relation between a carbon emission stream of a power system and the system tide. And secondly, acquiring a system node-node path output distribution factor by adopting a directional path algorithm, and solving the average influence factor of wind power injection power on other nodes and branches by combining a conventional unit and a system node and branch correlation matrix calculation model on the premise of assuming that wind power fluctuation is borne by a balance unit, thereby acquiring the uncertain characteristics of the carbon emission flow of the power system under random and intermittent wind power injection.
However, the calculation and analysis methods are performed on the basis of direct current power flow, and larger calculation errors are generated when the actual lossy network is faced, while in the prior art, the research on wind power uncertainty is lacking in a carbon flow analysis method for considering system loss, and most of carbon flow calculation methods for considering system loss are complex.
Therefore, it is necessary to provide a power grid carbon flow calculation scheme that takes into account system losses and wind power uncertainty.
Disclosure of Invention
The invention provides a power grid carbon flow calculation method and device considering system loss and wind power uncertainty, and solves the technical problem of how to realize power grid carbon flow calculation considering system loss and wind power uncertainty.
The first aspect of the invention provides a power grid carbon flow calculation method considering system loss and wind power uncertainty, comprising the following steps:
constructing a random wind speed sample set aiming at a target power system;
based on each group of random wind speed samples in the random wind speed sample set, carrying out deterministic power flow calculation accounting for active loss, and determining a functional expression between wind power injection power considering wind speed uncertainty and the total carbon flow rate of the target power system;
calculating related carbon flow indexes according to the obtained deterministic power flow calculation result; the related carbon flow indexes comprise a balance unit-node carbon flow incidence matrix, a balance unit-branch carbon flow incidence matrix, a unit-branch loss carbon flow incidence matrix, a node output distribution matrix and a path output distribution matrix;
Performing carbon flow analysis based on the relevant carbon flow index; the carbon flow analysis includes calculating contributions of balancing units to active power flows of portions of the target power system and calculating carbon flow rate influence factors of balancing units to nodes and branches.
According to one implementation manner of the first aspect of the present invention, the constructing a random wind speed sample set for a target power system includes:
estimating a probability density function of the wind speed by considering the influence of the wind speed uncertainty on the wind power injection power of the target power system;
obtaining a cumulative distribution function from the probability density function of the wind speed; the cumulative distribution function obeys uniform distribution between 0 and 1;
dividing the value range of the cumulative distribution function into a plurality of disjoint intervals, extracting a random number from each interval according to uniform distribution, and obtaining the random wind speed sample set without repeated extraction in each interval.
According to one implementation manner of the first aspect of the present invention, the determining a functional expression between wind power injection power considering wind speed uncertainty and total carbon flow rate of the target power system based on the deterministic power flow calculation accounting for active loss of each group of random wind speed samples in the random wind speed sample set includes:
Calculating the mechanical power of a fan under each group of random wind speed samples, and determining the relation between the mechanical power of the wind turbine and the wind speed according to the obtained mechanical power calculation result;
obtaining a function expression of wind power injection power and wind speed based on the relation between the mechanical power and wind speed of the wind turbine;
based on each load flow calculation equation, carrying out deterministic load flow calculation accounting for active loss on the target power system; the load flow calculation equations comprise that the injected active power of each node except the balance node is equal to the sum of the power difference value and the active loss of the corresponding node, the sum of the injected active power of each node except the balance node is equal to the active output of the corresponding node, and the total carbon flow rate of the target power system is equal to the sum of the carbon emission generated in unit time of all units; the power difference value is the difference value between the injection power of the access unit and the load injection power;
the simultaneous tidal current calculation equations result in a functional expression between wind power injection power taking into account wind speed uncertainty and the total carbon flow rate of the target power system.
According to one possible implementation manner of the first aspect of the present invention, the determining a relationship between mechanical power and wind speed of a wind turbine according to a result of obtaining the mechanical power calculation includes:
Describing a relationship between mechanical power and wind speed of the wind turbine based on a power curve of the wind turbine; the power curve of the wind turbine is expressed as an equation of the form:
；
in the method, in the process of the invention,P v indicating wind speedvThe mechanical power at the time of the manufacture,ρfor air density, a is the wind turbine area,C p is the power coefficient of the wind turbine,v cut-in for the cut-in wind speed,v _r for the rated wind speed,P _r for the rated power of the electric motor,v cut-out to cut out wind speed.
According to one possible implementation manner of the first aspect of the present invention, the simultaneous tidal current calculation equations obtain a functional expression between wind power injection power considering wind speed uncertainty and a total carbon flow rate of the target power system, including:
and combining the tide calculation equations to obtain a functional expression between wind power injection power considering wind speed uncertainty and the total carbon flow rate of the target power system, wherein the functional expression is as follows:
；
in the method, in the process of the invention,representing a target power systemTotal carbon flow rate,/->In order to balance the carbon emission intensity of the unit,P i is the target power systemiThe injected active power of the individual nodes,P i_loss is the firstiThe active power loss of the individual nodes,P Lm wind power access nodemIs used for the active load of the (a),srepresenting a balancing node in the target power system, +.>Is the number of the conventional units and sets,Nfor the number of nodes of the target power system, +. >Is the carbon emission intensity of the conventional unit,P Gk is the firstkThe injection power of the connected unit on each node,vfor the wind speed of the wind,v cut-in for the cut-in wind speed,v cut-out in order to cut out the wind speed,athe number of the wind power machines is the number of the wind power machines,ρfor air density, a is the wind turbine area,C p is the power coefficient of the wind turbine,P _r for the rated power of the electric motor,v _r is rated wind speed.
According to one implementation manner of the first aspect of the present invention, the calculating the contribution of the balancing unit to the active power flow of each part of the target power system includes:
defining an upstream distribution matrix based on a countercurrent flow method;
and calculating the contribution quantity of each unit corresponding to the active power flow consumed by the consumer in the target power system by combining the upstream distribution matrix, wherein the active power output provided by each unit in the system network loss is calculated according to the following formula:
；
in the method, in the process of the invention,P Gk-l for units in system network losskThe active force to be provided is that,P Gk is the firstkThe injection power of the connected unit on each node,to represent the first in the path output distribution matrixkAll the elements of the column are listed,P Nii is a unitkWind power access node of (a)iAt the position ofNActive flux matrix of order diagonal array nodeP N Is a function of the corresponding element of the (c),ξ N is thatNThe rank vector and all elements in the vector are 1,P l for branches in a target power systemlIs marked by the active loss of (2) TRepresenting a transpose;
wind power access node in calculation target power systemiIs the first of (2)kThe station set provides active output force for loads connected with system nodes:
；
in the method, in the process of the invention,P Gk-L wind power access node for target power systemiIs the first of (2)kThe station set provides active output force for loads connected with system nodes,ξ M is thatMThe rank vector and all elements in the vector are 1,P L is thatM×NThe order-load distribution matrix is used to determine,P L the element in (2) is the connection relation between all the electric loads and the target electric power system and the active load quantity.
According to one implementation manner of the first aspect of the present invention, the calculating the carbon flow rate influence factor of the balancing unit on each node and branch includes:
in the random tide analysis, only wind power uncertainty is considered temporarily and wind power fluctuation power is all borne by a balance unit, namely
；
In the method, in the process of the invention,P m target power system node for wind turbine generator system accessmIs used for the wind power injection power of the (1),P Gs balancing the output power of the unit;
on the basis of the association relation between the balance unit and nodes and branches of the target power system, analyzing and calculating the target power system to determine carbon flow rate influence factors of the balance unit on each node and branch:
is provided withR us-N To balance the unit-node carbon flow correlation matrix, namely:
；
In the method, in the process of the invention,to balance the carbon emission intensity of the unit +.>Output a distribution matrix for the path, ">Is 1 XN-dimensional row vector, corresponding to balancing machine setsThe element value of (2) is 1, and the other element values are 0;
is provided withZ us-N To balance the influence factor vector of the unit-node, namely:
；
the effect of wind power injection power on the carbon flow rate flux of each node of the target power system is expressed as:
；
in the method, in the process of the invention,representing the carbon flow rate flux of each node of the target power system;
the effect of wind power injection power on the carbon flow rate distribution of each branch of the target power system is expressed as:
；
in the method, in the process of the invention,representing the carbon flow rate distribution for each leg of the target electrical power system,R us-B in order to balance the set-branch carbon flow correlation matrix,Houtputting a distribution matrix for the nodes;
the effect of wind power injection power on the carbon flow rate distribution of each load of the target power system is expressed as:
；
in the method, in the process of the invention,representing the carbon flow rate distribution for each load of the target electrical power system,R us-L for the unit-branch loss carbon flow correlation matrix,ξ N is thatNThe rank vector and all elements in the vector are 1,P L is thatM×NThe order-load distribution matrix is used to determine,P L the elements in the system are the connection relation between all the electric loads and the target electric power system and the active load quantity,P N is thatNAn active flux matrix of the order diagonal array node;
the effect of wind power injection power on the carbon flow rate distribution of the system active loss is expressed as:
；
In the method, in the process of the invention,a carbon flow rate distribution representing the active loss of the system,P l for branches in a target power systemlIs not limited to the active loss of the (a).
A second aspect of the present invention provides a power grid carbon flow calculation device taking into account system loss and wind power uncertainty, comprising:
the construction module is used for constructing a random wind speed sample set aiming at a target power system;
the first calculation analysis module is used for carrying out deterministic power flow calculation accounting for active loss based on each group of random wind speed samples in the random wind speed sample set, and determining a functional expression between wind power injection power considering wind speed uncertainty and the total carbon flow rate of the target power system;
the second calculation analysis module is used for calculating the related carbon flow index according to the obtained deterministic power flow calculation result; the related carbon flow indexes comprise a balance unit-node carbon flow incidence matrix, a balance unit-branch carbon flow incidence matrix, a unit-branch loss carbon flow incidence matrix, a node output distribution matrix and a path output distribution matrix;
the third calculation analysis module is used for carrying out carbon flow analysis based on the related carbon flow index; the carbon flow analysis includes calculating contributions of balancing units to active power flows of portions of the target power system and calculating carbon flow rate influence factors of balancing units to nodes and branches.
According to one manner in which the second aspect of the present invention can be implemented, the building block includes:
the estimating unit is used for estimating a probability density function of the wind speed in consideration of the influence of the wind speed uncertainty on the wind power injection power of the target power system;
the function conversion unit is used for obtaining a cumulative distribution function from the probability density function of the wind speed; the cumulative distribution function obeys uniform distribution between 0 and 1;
the sample set determining unit is used for dividing the value range of the cumulative distribution function into a plurality of disjoint intervals, extracting a random number from each interval according to uniform distribution, and obtaining the random wind speed sample set without repeated extraction in each interval.
According to one implementation manner of the second aspect of the present invention, the first calculation and analysis module includes:
the first determining unit is used for calculating the mechanical power of the fan under each group of random wind speed samples and determining the relation between the mechanical power of the wind turbine and the wind speed according to the obtained mechanical power calculation result;
the second determining unit is used for obtaining a functional expression of wind power injection power and wind speed based on the relation between the mechanical power and wind speed of the wind turbine;
A calculation unit for performing deterministic power flow calculation accounting for active loss on the target power system based on each power flow calculation equation; the load flow calculation equations comprise that the injected active power of each node except the balance node is equal to the sum of the power difference value and the active loss of the corresponding node, the sum of the injected active power of each node except the balance node is equal to the active output of the corresponding node, and the total carbon flow rate of the target power system is equal to the sum of the carbon emission generated in unit time of all units; the power difference value is the difference value between the injection power of the access unit and the load injection power;
and the third determining unit is used for combining the tide calculation equations to obtain a functional expression between the wind power injection power considering the wind speed uncertainty and the total carbon flow rate of the target power system.
According to one possible implementation manner of the second aspect of the present invention, the first determining unit is specifically configured to:
describing a relationship between mechanical power and wind speed of the wind turbine based on a power curve of the wind turbine; the power curve of the wind turbine is expressed as an equation of the form:
；
in the method, in the process of the invention,P v indicating wind speedvThe mechanical power at the time of the manufacture, ρFor air density, a is the wind turbine area,C p is the power coefficient of the wind turbine,v cut-in for the cut-in wind speed,v _r for the rated wind speed,P _r for the rated power of the electric motor,v cut-out to cut out wind speed.
According to one possible implementation manner of the second aspect of the present invention, the third determining unit is specifically configured to:
and combining the tide calculation equations to obtain a functional expression between wind power injection power considering wind speed uncertainty and the total carbon flow rate of the target power system, wherein the functional expression is as follows:
；
in the method, in the process of the invention,representing the total carbon flow rate of the target power system, +.>In order to balance the carbon emission intensity of the unit,P i is the target power systemiThe injected active power of the individual nodes,P i_loss is the firstiThe active power loss of the individual nodes,P Lm wind power access nodemIs used for the active load of the (a),srepresenting a balancing node in the target power system, +.>Is the number of the conventional units and sets,Nfor the number of nodes of the target power system, +.>Is the carbon emission intensity of the conventional unit,P Gk is the firstkThe injection power of the connected unit on each node,vfor the wind speed of the wind,v cut-in for the cut-in wind speed,v cut-out in order to cut out the wind speed,athe number of the wind power machines is the number of the wind power machines,ρfor air density, a is the wind turbine area,C p is the power coefficient of the wind turbine,P _r for the rated power of the electric motor,v _r is rated wind speed.
According to one possible manner of the second aspect of the present invention, when the third calculation and analysis module calculates the contribution of the balancing unit to the active power flow of each part of the target power system, the third calculation and analysis module is specifically configured to:
Defining an upstream distribution matrix based on a countercurrent flow method;
and calculating the contribution quantity of each unit corresponding to the active power flow consumed by the consumer in the target power system by combining the upstream distribution matrix, wherein the active power output provided by each unit in the system network loss is calculated according to the following formula:
；
in the method, in the process of the invention,P Gk-l for units in system network losskThe active force to be provided is that,P Gk is the firstkThe injection power of the connected unit on each node,to represent the first in the path output distribution matrixkAll the elements of the column are listed,P Nii is a unitkWind power access node of (a)iAt the position ofNActive flux matrix of order diagonal array nodeP N Is a function of the corresponding element of the (c),ξ N is thatNThe rank vector and all elements in the vector are 1,P l for branches in a target power systemlIs marked by the active loss of (2)TRepresenting a transpose;
wind power access node in calculation target power systemiIs the first of (2)kThe station set provides active output force for loads connected with system nodes:
；
in the method, in the process of the invention,P Gk-L wind power access node for target power systemiIs the first of (2)kThe station set provides active output force for loads connected with system nodes,ξ M is thatMThe rank vector and all elements in the vector are 1,P L is thatM×NThe order-load distribution matrix is used to determine,P L the element in (2) is the connection relation between all the electric loads and the target electric power system and the active load quantity.
According to one implementation manner of the second aspect of the present invention, when the third calculation and analysis module calculates the carbon flow rate influence factor of the balancing unit on each node and branch, the third calculation and analysis module is specifically configured to:
in the random tide analysis, only wind power uncertainty is considered temporarily and wind power fluctuation power is all borne by a balance unit, namely
；
In the method, in the process of the invention,P m target power system node for wind turbine generator system accessmIs used for the wind power injection power of the (1),P Gs balancing the output power of the unit;
on the basis of the association relation between the balance unit and nodes and branches of the target power system, analyzing and calculating the target power system to determine carbon flow rate influence factors of the balance unit on each node and branch:
is provided withR us-N To balance the unit-node carbon flow correlation matrix, namely:
；
in the method, in the process of the invention,to balance the carbon emission intensity of the unit +.>Output a distribution matrix for the path, ">Is 1 XN-dimensional row vector, corresponding to balancing machine setsThe element value of (2) is 1, and the other element values are 0;
is provided withZ us-N To balance the influence factor vector of the unit-node, namely:
；
the effect of wind power injection power on the carbon flow rate flux of each node of the target power system is expressed as:
；
in the method, in the process of the invention,representing the carbon flow rate flux of each node of the target power system;
The effect of wind power injection power on the carbon flow rate distribution of each branch of the target power system is expressed as:
；
in the method, in the process of the invention,representing the carbon flow rate distribution for each leg of the target electrical power system,R us-B in order to balance the set-branch carbon flow correlation matrix,Houtputting a distribution matrix for the nodes;
the effect of wind power injection power on the carbon flow rate distribution of each load of the target power system is expressed as:
；
in the method, in the process of the invention,representing the carbon flow rate distribution for each load of the target electrical power system,R us-L for the unit-branch loss carbon flow correlation matrix,ξ N is thatNThe rank vector and all elements in the vector are 1,P L is thatM×NThe order-load distribution matrix is used to determine,P L the elements in the system are the connection relation between all the electric loads and the target electric power system and the active load quantity,P N is thatNAn active flux matrix of the order diagonal array node;
the effect of wind power injection power on the carbon flow rate distribution of the system active loss is expressed as:
；
in the method, in the process of the invention,a carbon flow rate distribution representing the active loss of the system,P l for branches in a target power systemlIs not limited to the active loss of the (a).
A third aspect of the present invention provides a power grid carbon flow calculation device taking into account system loss and wind power uncertainty, comprising:
a memory for storing instructions; the instructions are used for realizing the power grid carbon flow calculation method considering the system loss and wind power uncertainty in the mode that any one of the above can be realized;
And the processor is used for executing the instructions in the memory.
A fourth aspect of the present invention is a computer readable storage medium, where a computer program is stored, where the computer program when executed by a processor implements a method for calculating a grid carbon flow, where the system loss and the wind power uncertainty are considered as described in any one of the above modes.
From the above technical scheme, the invention has the following advantages:
the method comprises the steps of constructing a random wind speed sample set aiming at a target power system; based on each group of random wind speed samples in the random wind speed sample set, carrying out deterministic power flow calculation accounting for active loss, and determining a functional expression between wind power injection power considering wind speed uncertainty and the total carbon flow rate of the target power system; calculating related carbon flow indexes according to the obtained deterministic power flow calculation result; the related carbon flow indexes comprise a balance unit-node carbon flow incidence matrix, a balance unit-branch carbon flow incidence matrix, a unit-branch loss carbon flow incidence matrix, a node output distribution matrix and a path output distribution matrix; performing carbon flow analysis based on the relevant carbon flow index; the carbon flow analysis comprises the steps of calculating the contribution of a balance unit to active power flow of each part of the target power system and calculating the carbon flow rate influence factors of the balance unit to each node and each branch; the method realizes the calculation of the carbon flow of the power grid taking the system loss and wind power uncertainty into consideration, and can provide a theoretical basis for quantitatively evaluating the influence of wind power access on the carbon emission of the system.
Drawings
In order to more clearly illustrate the embodiments of the invention or the technical solutions in the prior art, the drawings that are required in the embodiments or the description of the prior art will be briefly described, it being obvious that the drawings in the following description are only some embodiments of the 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 flowchart of a method for calculating a power grid carbon flow taking into account system loss and wind power uncertainty according to an alternative embodiment of the present invention;
FIG. 2 is a block diagram of a power grid carbon flow calculation device with consideration of system loss and wind power uncertainty according to an alternative embodiment of the present invention.
Reference numerals:
1-building a module; 2-a first computational analysis module; 3-a second computational analysis module; 4-a third calculation analysis module.
Detailed Description
The embodiment of the invention provides a power grid carbon flow calculation method and device considering system loss and wind power uncertainty, which are used for solving the technical problem of how to realize power grid carbon flow calculation considering system loss and wind power uncertainty.
In order to make the objects, features and advantages of the present invention more comprehensible, the technical solutions in the embodiments of the present invention are described in detail below with reference to the accompanying drawings, and it is apparent that the embodiments described below are only some embodiments of the present invention, but not all embodiments of the present invention. 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 invention provides a power grid carbon flow calculation method considering system loss and wind power uncertainty.
Referring to fig. 1, fig. 1 shows a flowchart of a power grid carbon flow calculation method considering system loss and wind power uncertainty according to an embodiment of the present invention.
The power grid carbon flow calculation method considering system loss and wind power uncertainty provided by the embodiment of the invention comprises the steps S1-S4.
Step S1, constructing a random wind speed sample set aiming at a target power system.
Wherein a Latin hypercube sampling method can be used to generate a wind power sample set related to wind speed from uniform distributionU sample As a random wind speed sample set for a target power system, wherein the number of samplesnThe setting can be performed according to actual conditions.
As an embodiment, the constructing a random wind speed sample set for a target power system includes:
taking into account wind speedvInfluence of uncertainty on wind power injection power, and probability density function of estimated wind speed isf(v)；
From probability density function of wind speedf(v) Obtaining a cumulative distribution functionU=F(v) Then the distribution function is accumulatedU=F(v) Obeys [0,1 ]]The space is uniformly distributed;
integrating the distribution functionUThe value range [0,1 ]]Aliquoting intonSubintervals, i.e. divided into nThe disjoint intervals [0,1 ]n），[1/n，2/n)，…，[(n-1)/n，1]Extracting a random number in each interval according to uniform distribution, and not repeatedly extracting each subinterval to obtain a sample matrixU sample WhereinnIs the number of samples.
And S2, carrying out deterministic power flow calculation accounting for active loss based on each group of random wind speed samples in the random wind speed sample set, and determining a functional expression between wind power injection power considering wind speed uncertainty and the total carbon flow rate of the target power system.
In one implementation, the determining a functional expression between wind power injection power accounting for wind speed uncertainty and total carbon flow rate of the target power system based on deterministic power flow calculation accounting for active losses for each set of random wind speed samples in the set of random wind speed samples includes:
calculating the mechanical power of a fan under each group of random wind speed samples, and determining the relation between the mechanical power of the wind turbine and the wind speed according to the obtained mechanical power calculation result;
obtaining a function expression of wind power injection power and wind speed based on the relation between the mechanical power and wind speed of the wind turbine;
based on each load flow calculation equation, carrying out deterministic load flow calculation accounting for active loss on the target power system; the load flow calculation equations comprise that the injected active power of each node except the balance node is equal to the sum of the power difference value and the active loss of the corresponding node, the sum of the injected active power of each node except the balance node is equal to the active output of the corresponding node, and the total carbon flow rate of the target power system is equal to the sum of the carbon emission generated in unit time of all units; the power difference value is the difference value between the injection power of the access unit and the load injection power;
The simultaneous tidal current calculation equations result in a functional expression between wind power injection power taking into account wind speed uncertainty and the total carbon flow rate of the target power system.
In one possible implementation, the determining the relationship between the mechanical power and the wind speed of the wind turbine according to the obtained mechanical power calculation result includes:
describing a relationship between mechanical power and wind speed of the wind turbine based on a power curve of the wind turbine; the power curve of the wind turbine is expressed as an equation of the form:
（1）
in the method, in the process of the invention,P v indicating wind speedvTime machineThe mechanical power of the device is calculated,ρfor air density, a is the wind turbine area,C p is the power coefficient of the wind turbine,v cut-in for the cut-in wind speed,v _r for the rated wind speed,P _r for the rated power of the electric motor,v cut-out to cut out wind speed.
Assume that the node number of the system isNWind turbine generator access system nodem(m=1,2,…,N) Wind power injection powerP m The expression of the function with wind speed is
（2）
Wherein, the liquid crystal display device comprises a liquid crystal display device,athe number of wind turbines is the number.
The conventional units and loads of the system are all known quantities, and the injection power of each node except the balance node is the sum of the injection power of the access unit and the injection power of the load of the node, namely:
（3）
in the method, in the process of the invention,P i is the target power systemi(i=1,2,…,N) The injected active power of the individual nodes, P Gi Is the firstiInjection power of the connected units on the individual nodes (if no unit is connectedP Gi =0），P Li For the active load of the node (if no load is accessedP Li =0），P i_loss Is the firstiActive loss of individual nodes;
suppose a nodesIs a balance node in the target power system, then the sum of node injection powers of all other nodes is the active power output of that node, i.e., the nodesThe injection power of (2) is:
（4）
in the method, in the process of the invention,P m wind power access nodemIs used for the injection of active power,P m_loss wind power access nodemActive loss of (2);
for the entire target power system, the sum of the carbon emissions per unit time, i.e., the node carbon flow rates, is equal to the carbon emissions produced per unit time for each unit, i.e., the total system carbon flow rate is:
（5）
in the method, in the process of the invention,Kis the number of the conventional units,P Gk is the firstk(k=1,2,…,N) The injection power of the conventional unit on the individual nodes,e Gs ande Gk the carbon emission intensity of the balancing unit and the conventional unit is respectively;
the simultaneous equation may yield a relationship of injection power to total node carbon flow rate of the system taking into account wind speed uncertainty,
（6）
in the method, in the process of the invention,representing the total carbon flow rate of the target power system, +.>In order to balance the carbon emission intensity of the unit,P i is the target power systemiThe injected active power of the individual nodes,P i_loss is the firstiThe active power loss of the individual nodes, P Lm Wind power access nodemIs used for the active load of the (a),srepresenting a balancing node in the target power system, +.>Is the number of the conventional units and sets,Nfor the number of nodes of the target power system, +.>Is the carbon emission intensity of the conventional unit,P Gk is the firstkThe injection power of the connected unit on each node,vfor the wind speed of the wind,v cut-in for the cut-in wind speed,v cut-out in order to cut out the wind speed,athe number of the wind power machines is the number of the wind power machines,ρfor air density, a is the wind turbine area,C p is the power coefficient of the wind turbine,P _r for the rated power of the electric motor,v _r is rated wind speed.
The above shows that the wind power injection power can generate fluctuation on the total node carbon flow rate of the system under different conditions.
Step S3, calculating related carbon flow indexes according to the obtained deterministic power flow calculation result; the related carbon flow indexes comprise a balance unit-node carbon flow correlation matrix, a balance unit-branch carbon flow correlation matrix, a unit-branch loss carbon flow correlation matrix, a node output distribution matrix and a path output distribution matrix.
To give boundary conditions for carbon emission flow distribution from the power network level, a branch tidal current distribution matrix is definedP B If nodeiAnd nodej(i,j=1,2,…,N) Is connected by a branch and is connected from a node by the branchiTo the nodejThe forward active power flow flowing through is pThenP Bij =p，P Bji =0; in other casesP Bij =P Bji =0；
Defining a set injection distribution matrixP G To describe the connection relation between all generator sets and the power system and define the active power injected into the system by the generator setsP L To describe the connection relation between all the electric loads and the electric system and the active load quantity, define the active loss carbon flow rate matrixR l To represent the active loss of each branch in the systemP l The resulting carbon flow rate:
（7）
set the firsti(i=1,2,…,N) The carbon potential of each node ise Ni Then the node carbon potential vectorE N =[e N1 e N2 …e NN ] T ；
Definition of the definitionNActive flux matrix of order diagonal array nodeP N When the system isP B AndP G when the matrix is known to be a function of the matrix,P N can pass throughP B AndP G matrix is directly generated, orderP Z =[P B P G ] T Then:
（8）
in the method, in the process of the invention,ξ N K+ is thatN+KAn order vector, wherein all elements in the vector are 1 (the same applies below);
further describing distribution information of tide and carbon emission stream of the system under given running steady state, and defining node output distribution matrixH：
（9）
In the method, in the process of the invention,Iis a unit matrix;
in order to characterize the flow of the power generation unit injection system and the path information of the carbon emission flow flowing from the node where the power generation unit is positioned to the target node, a path output distribution matrix of the system injection is definedDBasic calculation method and matrix for system node carbon potentialDCan be defined as follows:
（10）
determining contribution of carbon flow injection to carbon flow rate into another node for all gensets in a system Artificial machine set-node carbon flow association matrixR U-N ：
（11）
In the method, in the process of the invention,E G is a carbon emission intensity vector;
s47, the first system is available in the same waykUnit-branch carbon flow correlation matrix of contribution condition of carbon flow rate of all branches in carbon flow injection system of generator unitR U-B k, And the unit-node carbon flow correlation matrix of the contribution condition of the carbon flow injection of all the generator units to the load of other nodes corresponding to the carbon flow rate in the systemR U-L ：
（12）
In the method, in the process of the invention,to represent a matrixDMiddle (f)kAll elements are listed.
S4, carrying out carbon flow analysis based on the related carbon flow index; the carbon flow analysis includes calculating contributions of balancing units to active power flows of portions of the target power system and calculating carbon flow rate influence factors of balancing units to nodes and branches.
In one manner that can be implemented, the calculating the contribution of the balancing unit to the active power flow of each part of the target power system includes:
the first of the system is defined as the unit-branch carbon flow associated matrixkUnit-branch loss carbon flow correlation matrix for distribution of carbon flow rates caused by active loss of each branch in a station generator unitR Uk-l ：
（13）
The countercurrent flow method is to add virtual nodes in the branch and equivalent the branch loss as virtual negative The load, first define the upstream distribution matrixA u The elements are as follows:
（14）
in the method, in the process of the invention,i,j=1,2,…,N，U i is a nodeiUpstream node sets of (a);
since the branch loss is equivalent to a virtual load, the total loss of the network is obtained by adding the power. By combining the upstream distribution matrix, the contribution of each unit corresponding to the active power flow consumed by the consumer in the system can be obtained, and the active power output provided by each unit in the system network loss can be represented by the following formula:
（15）/>
in the method, in the process of the invention,P Gk-l for units in system network losskThe active force to be provided is that,P Gk is the firstkThe injection power of the connected unit on each node,to represent the first in the path output distribution matrixkAll the elements of the column are listed,P Nii is a unitkWind power access node of (a)iAt the position ofNActive flux matrix of order diagonal array nodeP N Is a function of the corresponding element of the (c),ξ N is thatNThe rank vector and all elements in the vector are 1,P l for branches in a target power systemlIs marked by the active loss of (2)TRepresenting a transpose;
similarly, the data of the active power flow of the system unit to each load of the system can be obtained, and the access node in the systemiIs the first of (2)kThe active output force provided by the station set for the load connected with the system node can be expressed by the following formula:
（16）
in the method, in the process of the invention,P Gk-L wind power access node for target power systemiIs the first of (2)kThe station set provides active output force for loads connected with system nodes, ξ M Is thatMThe rank vector and all elements in the vector are 1,P L is thatM×NThe order-load distribution matrix is used to determine,P L the element in (2) is the connection relation between all the electric loads and the target electric power system and the active load quantity.
In one manner that can be implemented, the calculating the carbon flow rate impact factor of the balancing unit on each node and branch includes:
in the random tide analysis, only wind power uncertainty is considered temporarily and wind power fluctuation power is all borne by a balance unit, namely
；
In the method, in the process of the invention,P m target power system node for wind turbine generator system accessmIs used for the wind power injection power of the (1),P Gs balancing the output power of the unit;
on the basis of the association relation between the balance unit and nodes and branches of the target power system, analyzing and calculating the target power system to determine carbon flow rate influence factors of the balance unit on each node and branch:
is provided withR us-N To balance the unit-node carbon flow correlation matrix, namely:
；
in the method, in the process of the invention,to balance the carbon emission intensity of the unit +.>Output a distribution matrix for the path, ">Is 1 XN-dimensional row vector, corresponding to balancing machine setsThe element value of (2) is 1, and the other element values are 0;
is provided withZ us-N To balance the influence factor vector of the unit-node, namely:
；
the effect of wind power injection power on the carbon flow rate flux of each node of the target power system is expressed as:
；
In the method, in the process of the invention,representing the carbon flow rate flux of each node of the target power system;
the effect of wind power injection power on the carbon flow rate distribution of each branch of the target power system is expressed as:
；
in the method, in the process of the invention,representing the carbon flow rate distribution for each leg of the target electrical power system,R us-B in order to balance the set-branch carbon flow correlation matrix,Houtputting a distribution matrix for the nodes;
the effect of wind power injection power on the carbon flow rate distribution of each load of the target power system is expressed as:
；
in the method, in the process of the invention,representing the carbon flow rate distribution for each load of the target electrical power system,R us-L for the unit-branch loss carbon flow correlation matrix,ξ N is thatNThe rank vector and all elements in the vector are 1,P L is thatM×NThe order-load distribution matrix is used to determine,P L the elements in the system are the connection relation between all the electric loads and the target electric power system and the active load quantity,P N is thatNAn active flux matrix of the order diagonal array node;
the effect of wind power injection power on the carbon flow rate distribution of the system active loss is expressed as:
；
in the method, in the process of the invention,a carbon flow rate distribution representing the active loss of the system,P l for branches in a target power systemlIs not limited to the active loss of the (a).
According to the embodiment of the invention, a carbon emission flow distribution analysis method of random alternating current flow considering wind speed uncertainty is provided, a correlation function of wind speed, wind direction and total carbon flow rate of a system is obtained, and an uncertainty analysis environment is developed for deterministic carbon emission flow analysis; the method for analyzing the carbon emission flow of the power system based on the countercurrent flow method and considering the active loss of the line is provided, and on the basis, the average influence factors of balance nodes on the carbon flow rates of other nodes and branches under the wind power injection power are provided; and taking the wind speed uncertainty characteristic into consideration, the average and variance statistical characteristics of the node carbon potential and the average influence factor of wind power injection power on the carbon flow rate can be obtained, so that the carbon emission flow distribution characteristic considering the wind power uncertainty influence is obtained.
The invention also provides a power grid carbon flow calculation device considering the system loss and the wind power uncertainty, which can be used for executing the power grid carbon flow calculation method considering the system loss and the wind power uncertainty according to any one of the embodiments of the invention.
Referring to fig. 2, fig. 2 shows a block diagram of a structural connection of a power grid carbon flow calculation device taking system loss and wind power uncertainty into consideration according to an embodiment of the invention.
The embodiment of the invention provides a power grid carbon flow calculation device considering system loss and wind power uncertainty, which comprises the following components:
a construction module 1 for constructing a random wind speed sample set for a target power system;
a first calculation and analysis module 2, configured to perform deterministic power flow calculation accounting for active loss based on each group of random wind speed samples in the random wind speed sample set, and determine a functional expression between wind power injection power considering wind speed uncertainty and a total carbon flow rate of the target power system;
the second calculation and analysis module 3 is used for calculating the related carbon flow index according to the obtained deterministic power flow calculation result; the related carbon flow indexes comprise a balance unit-node carbon flow incidence matrix, a balance unit-branch carbon flow incidence matrix, a unit-branch loss carbon flow incidence matrix, a node output distribution matrix and a path output distribution matrix;
A third calculation and analysis module 4, configured to perform carbon flow analysis based on the related carbon flow index; the carbon flow analysis includes calculating contributions of balancing units to active power flows of portions of the target power system and calculating carbon flow rate influence factors of balancing units to nodes and branches.
In one possible implementation, the building block 1 comprises:
the estimating unit is used for estimating a probability density function of the wind speed in consideration of the influence of the wind speed uncertainty on the wind power injection power of the target power system;
the function conversion unit is used for obtaining a cumulative distribution function from the probability density function of the wind speed; the cumulative distribution function obeys uniform distribution between 0 and 1;
the sample set determining unit is used for dividing the value range of the cumulative distribution function into a plurality of disjoint intervals, extracting a random number from each interval according to uniform distribution, and obtaining the random wind speed sample set without repeated extraction in each interval.
In one possible implementation, the first computational analysis module 2 comprises:
the first determining unit is used for calculating the mechanical power of the fan under each group of random wind speed samples and determining the relation between the mechanical power of the wind turbine and the wind speed according to the obtained mechanical power calculation result;
The second determining unit is used for obtaining a functional expression of wind power injection power and wind speed based on the relation between the mechanical power and wind speed of the wind turbine;
a calculation unit for performing deterministic power flow calculation accounting for active loss on the target power system based on each power flow calculation equation; the load flow calculation equations comprise that the injected active power of each node except the balance node is equal to the sum of the power difference value and the active loss of the corresponding node, the sum of the injected active power of each node except the balance node is equal to the active output of the corresponding node, and the total carbon flow rate of the target power system is equal to the sum of the carbon emission generated in unit time of all units; the power difference value is the difference value between the injection power of the access unit and the load injection power;
and the third determining unit is used for combining the tide calculation equations to obtain a functional expression between the wind power injection power considering the wind speed uncertainty and the total carbon flow rate of the target power system.
In one implementation manner, the first determining unit is specifically configured to:
describing a relationship between mechanical power and wind speed of the wind turbine based on a power curve of the wind turbine; the power curve of the wind turbine is expressed as an equation of the form:
；
In the method, in the process of the invention,P v indicating wind speedvThe mechanical power at the time of the manufacture,ρfor air density, a is the wind turbine area,C p is the power coefficient of the wind turbine,v cut-in for the cut-in wind speed,v _r for the rated wind speed,P _r for the rated power of the electric motor,v cut-out to cut out wind speed.
In one implementation manner, the third determining unit is specifically configured to:
and combining the tide calculation equations to obtain a functional expression between wind power injection power considering wind speed uncertainty and the total carbon flow rate of the target power system, wherein the functional expression is as follows:
；/>
in the method, in the process of the invention,representing the total carbon flow rate of the target power system, +.>In order to balance the carbon emission intensity of the unit,P i is the target power systemiThe injected active power of the individual nodes,P i_loss is the firstiThe active power loss of the individual nodes,P Lm wind power access nodemIs used for the active load of the (a),srepresenting a balancing node in the target power system, +.>Is the number of the conventional units and sets,Nfor the number of nodes of the target power system, +.>Is the carbon emission intensity of the conventional unit,P Gk is the firstkThe injection power of the connected unit on each node,vfor the wind speed of the wind,v cut-in for the cut-in wind speed,v cut-out in order to cut out the wind speed,athe number of the wind power machines is the number of the wind power machines,ρfor air density, a is the wind turbine area,C p is the power coefficient of the wind turbine,P _r for the rated power of the electric motor,v _r is rated wind speed.
In one possible implementation manner, when the third calculation and analysis module 4 calculates the contribution of the balancing unit to the active power flow of each part of the target power system, the method is specifically used for:
defining an upstream distribution matrix based on a countercurrent flow method;
and calculating the contribution quantity of each unit corresponding to the active power flow consumed by the consumer in the target power system by combining the upstream distribution matrix, wherein the active power output provided by each unit in the system network loss is calculated according to the following formula:
；
in the method, in the process of the invention,P Gk-l for units in system network losskThe active force to be provided is that,P Gk is the firstkThe injection power of the connected unit on each node,to represent the first in the path output distribution matrixkAll the elements of the column are listed,P Nii is a unitkWind power access node of (a)iAt the position ofNActive flux matrix of order diagonal array nodeP N Is a function of the corresponding element of the (c),ξ N is thatNThe rank vector and all elements in the vector are 1,P l for branches in a target power systemlIs marked by the active loss of (2)TRepresenting a transpose;
wind power access node in calculation target power systemiIs the first of (2)kThe station set provides active output force for loads connected with system nodes:
；
in the method, in the process of the invention,P Gk-L wind power access node for target power systemiIs the first of (2)kThe station set provides active output force for loads connected with system nodes, ξ M Is thatMThe rank vector and all elements in the vector are 1,P L is thatM×NThe order-load distribution matrix is used to determine,P L the elements in the system are the connection relation between all the electric loads and the target power system and the active powerLoad amount.
In one possible implementation, the third calculation and analysis module 4 is specifically configured to, when calculating the carbon flow rate influence factor of the balancing unit on each node and branch:
in the random tide analysis, only wind power uncertainty is considered temporarily and wind power fluctuation power is all borne by a balance unit, namely
；
In the method, in the process of the invention,P m target power system node for wind turbine generator system accessmIs used for the wind power injection power of the (1),P Gs balancing the output power of the unit;
on the basis of the association relation between the balance unit and nodes and branches of the target power system, analyzing and calculating the target power system to determine carbon flow rate influence factors of the balance unit on each node and branch:
is provided withR us-N To balance the unit-node carbon flow correlation matrix, namely:
；/>
in the method, in the process of the invention,to balance the carbon emission intensity of the unit +.>Output a distribution matrix for the path, ">Is 1 XN-dimensional row vector, corresponding to balancing machine setsThe element value of (2) is 1, and the other element values are 0;
is provided withZ us-N To balance the influence factor vector of the unit-node, namely:
；
the effect of wind power injection power on the carbon flow rate flux of each node of the target power system is expressed as:
；
In the method, in the process of the invention,representing the carbon flow rate flux of each node of the target power system;
the effect of wind power injection power on the carbon flow rate distribution of each branch of the target power system is expressed as:
；
in the method, in the process of the invention,representing the carbon flow rate distribution for each leg of the target electrical power system,R us-B in order to balance the set-branch carbon flow correlation matrix,Houtputting a distribution matrix for the nodes;
the effect of wind power injection power on the carbon flow rate distribution of each load of the target power system is expressed as:
；
in the method, in the process of the invention,representing the carbon flow rate distribution for each load of the target electrical power system,R us-L for the unit-branch loss carbon flow correlation matrix,ξ N is thatNThe rank vector and all elements in the vector are 1,P L is thatM×NThe order-load distribution matrix is used to determine,P L the elements in the system are the connection relation between all the electric loads and the target electric power system and the active load quantity,P N is thatNAn active flux matrix of the order diagonal array node;
the effect of wind power injection power on the carbon flow rate distribution of the system active loss is expressed as:
；
in the method, in the process of the invention,a carbon flow rate distribution representing the active loss of the system,P l for branches in a target power systemlIs not limited to the active loss of the (a).
The invention also provides a power grid carbon flow calculation device considering system loss and wind power uncertainty, which comprises:
a memory for storing instructions; the instructions are used for realizing the power grid carbon flow calculation method considering system loss and wind power uncertainty according to any one of the embodiments;
And the processor is used for executing the instructions in the memory.
The invention also provides a computer readable storage medium, wherein the computer readable storage medium is stored with a computer program, and the computer program is executed by a processor to realize the grid carbon flow calculation method considering the system loss and wind power uncertainty according to any one of the embodiments.
It will be clearly understood by those skilled in the art that, for convenience and brevity of description, specific working processes of the above-described apparatus, modules and units may refer to corresponding processes in the foregoing method embodiments, and specific beneficial effects of the above-described apparatus, modules and units may refer to corresponding beneficial effects in the foregoing method embodiments, which are not repeated herein.
In the several embodiments provided in this application, it should be understood that the disclosed apparatus and method may be implemented in other ways. For example, the apparatus embodiments described above are merely illustrative, and for example, the division of the modules is merely a logical function division, and there may be additional divisions when actually implemented, for example, multiple modules or components may be combined or integrated into another apparatus, or some features may be omitted or not performed. Alternatively, the coupling or direct coupling or communication connection shown or discussed with each other may be an indirect coupling or communication connection via some interfaces, devices or modules, which may be in electrical, mechanical, or other forms.
The modules described as separate components may or may not be physically separate, and components shown as modules may or may not be physical modules, i.e., may be located in one place, or may be distributed over a plurality of network modules. Some or all of the modules may be selected according to actual needs to achieve the purpose of the solution of this embodiment.
In addition, each functional module in each embodiment of the present invention may be integrated into one processing module, or each module may exist alone physically, or two or more modules may be integrated into one module. The integrated modules may be implemented in hardware or in software functional modules.
The integrated modules, if implemented in the form of software functional modules and sold or used as a stand-alone product, may be stored in a computer readable storage medium. Based on such understanding, the technical solution of the present invention may be embodied essentially or in part or all of the technical solution or in part in the form of a software product stored in a storage medium, including instructions for causing a computer device (which may be a personal computer, a server, or a network device, etc.) to perform all or part of the steps of the method according to the embodiments of the present invention. And the aforementioned storage medium includes: a U-disk, a removable hard disk, a Read-Only Memory (ROM), a random access Memory (RAM, random Access Memory), a magnetic disk, or an optical disk, or other various media capable of storing program codes.
The above embodiments are only for illustrating the technical solution of the present invention, and not for limiting the same; although the invention has been described in detail with reference to the foregoing embodiments, it will be understood by those of ordinary skill in the art that: the technical scheme described in the foregoing embodiments can be modified or some technical features thereof can be replaced by equivalents; such modifications and substitutions do not depart from the spirit and scope of the technical solutions of the embodiments of the present invention.
Claims (10)
1. A power grid carbon flow calculation method considering system loss and wind power uncertainty is characterized by comprising the following steps:
constructing a random wind speed sample set aiming at a target power system;
based on each group of random wind speed samples in the random wind speed sample set, carrying out deterministic power flow calculation accounting for active loss, and determining a functional expression between wind power injection power considering wind speed uncertainty and the total carbon flow rate of the target power system;
calculating related carbon flow indexes according to the obtained deterministic power flow calculation result; the related carbon flow indexes comprise a balance unit-node carbon flow incidence matrix, a balance unit-branch carbon flow incidence matrix, a unit-branch loss carbon flow incidence matrix, a node output distribution matrix and a path output distribution matrix;
Performing carbon flow analysis based on the relevant carbon flow index; the carbon flow analysis includes calculating contributions of balancing units to active power flows of portions of the target power system and calculating carbon flow rate influence factors of balancing units to nodes and branches.
2. The method for calculating the carbon flow of the power grid taking into account the system loss and the wind power uncertainty as claimed in claim 1, wherein the constructing a random wind speed sample set for the target power system comprises:
estimating a probability density function of the wind speed by considering the influence of the wind speed uncertainty on the wind power injection power of the target power system;
obtaining a cumulative distribution function from the probability density function of the wind speed; the cumulative distribution function obeys uniform distribution between 0 and 1;
dividing the value range of the cumulative distribution function into a plurality of disjoint intervals, extracting a random number from each interval according to uniform distribution, and obtaining the random wind speed sample set without repeated extraction in each interval.
3. The method for calculating the carbon flow of the power grid taking into account the system loss and the wind power uncertainty as claimed in claim 1, wherein the step of calculating the deterministic power flow accounting for the active loss based on each group of random wind speed samples in the random wind speed sample set, determining a functional expression between the wind power injection power taking into account the wind speed uncertainty and the total carbon flow rate of the target power system, comprises the following steps:
Calculating the mechanical power of a fan under each group of random wind speed samples, and determining the relation between the mechanical power of the wind turbine and the wind speed according to the obtained mechanical power calculation result;
obtaining a function expression of wind power injection power and wind speed based on the relation between the mechanical power and wind speed of the wind turbine;
based on each load flow calculation equation, carrying out deterministic load flow calculation accounting for active loss on the target power system; the load flow calculation equations comprise that the injected active power of each node except the balance node is equal to the sum of the power difference value and the active loss of the corresponding node, the sum of the injected active power of each node except the balance node is equal to the active output of the corresponding node, and the total carbon flow rate of the target power system is equal to the sum of the carbon emission generated in unit time of all units; the power difference value is the difference value between the injection power of the access unit and the load injection power;
the simultaneous tidal current calculation equations result in a functional expression between wind power injection power taking into account wind speed uncertainty and the total carbon flow rate of the target power system.
4. A method of grid carbon flow calculation taking into account system losses and wind power uncertainty as defined in claim 3, wherein said determining a relationship between mechanical power and wind speed of a wind turbine from the obtained mechanical power calculations comprises:
Describing a relationship between mechanical power and wind speed of the wind turbine based on a power curve of the wind turbine; the power curve of the wind turbine is expressed as an equation of the form:
；
in the method, in the process of the invention,P v indicating wind speedvThe mechanical power at the time of the manufacture,ρfor air density, a is the wind turbine area,C p is the power coefficient of the wind turbine,v cut-in for the cut-in wind speed,v _r for the rated wind speed,P _r for the rated power of the electric motor,v cut-out to cut out wind speed.
5. A method of calculating a grid carbon flow taking into account system losses and wind power uncertainty as defined in claim 3, wherein said simultaneous tidal current calculation equations yield a functional expression between wind power injection power taking into account wind speed uncertainty and a total carbon flow rate of said target power system, comprising:
and combining the tide calculation equations to obtain a functional expression between wind power injection power considering wind speed uncertainty and the total carbon flow rate of the target power system, wherein the functional expression is as follows:
；
in the method, in the process of the invention,representing the total carbon flow rate of the target power system, +.>In order to balance the carbon emission intensity of the unit, P i is the target power systemiThe injected active power of the individual nodes,P i_loss is the firstiThe active power loss of the individual nodes,P Lm wind power access nodemIs used for the active load of the (a),srepresenting a balancing node in the target power system, +. >Is the number of the conventional units and sets,Nfor the number of nodes of the target power system, +.>Is the carbon emission intensity of the conventional unit,P Gk is the firstkThe injection power of the connected unit on each node,vfor the wind speed of the wind,v cut-in for the cut-in wind speed,v cut-out in order to cut out the wind speed,athe number of the wind power machines is the number of the wind power machines,ρfor air density, a is the wind turbine area,C p is the power coefficient of the wind turbine,P _r for the rated power of the electric motor,v _r is rated wind speed.
6. The method for calculating the carbon flow of the power grid taking system loss and wind power uncertainty into consideration as set forth in claim 1, wherein the calculating the contribution of the balancing unit to the active power flow of each part of the target power system includes:
defining an upstream distribution matrix based on a countercurrent flow method;
and calculating the contribution quantity of each unit corresponding to the active power flow consumed by the consumer in the target power system by combining the upstream distribution matrix, wherein the active power output provided by each unit in the system network loss is calculated according to the following formula:
；
in the method, in the process of the invention,P Gk-l for units in system network losskThe active force to be provided is that,P Gk is the firstkThe injection power of the connected unit on each node,to represent the first in the path output distribution matrixkAll the elements of the column are listed,P Nii is a unitkWind power access node of (a)iAt the position ofNActive flux matrix of order diagonal array node P N Is a function of the corresponding element of the (c),ξ N is thatNThe rank vector and all elements in the vector are 1,P l for branches in a target power systemlIs marked by the active loss of (2)TRepresenting a transpose;
wind power access node in calculation target power systemiIs the first of (2)kThe station set provides active output force for loads connected with system nodes:
；
in the method, in the process of the invention,P Gk-L wind power access node for target power systemiIs the first of (2)kThe station set provides active output force for loads connected with system nodes,ξ M is thatMThe rank vector and all elements in the vector are 1,P L is thatM×NThe order-load distribution matrix is used to determine,P L the element in (2) is the connection relation between all the electric loads and the target electric power system and the active load quantity.
7. The method for calculating the carbon flow of the power grid taking system loss and wind power uncertainty into consideration according to claim 1, wherein the calculating the carbon flow rate influence factor of the balancing unit on each node and branch comprises:
in the random tide analysis, only wind power uncertainty is considered temporarily and wind power fluctuation power is all borne by a balance unit, namely
；
In the method, in the process of the invention,P m target power system node for wind turbine generator system accessmIs used for the wind power injection power of the (1),P Gs balancing the output power of the unit;
on the basis of the association relation between the balance unit and nodes and branches of the target power system, analyzing and calculating the target power system to determine carbon flow rate influence factors of the balance unit on each node and branch:
Is provided withR us-N To balance the unit-node carbon flow correlation matrix, namely:
；
in the method, in the process of the invention,to balance the carbon emission intensity of the unit +.>Output a distribution matrix for the path, ">Is 1 XN-dimensional row vector, corresponding to balancing machine setsThe element value of (2) is 1, and the other element values are 0;
is provided withZ us-N To balance the influence factor vector of the unit-node, namely:
；
the effect of wind power injection power on the carbon flow rate flux of each node of the target power system is expressed as:
；
in the method, in the process of the invention,representing the carbon flow rate flux of each node of the target power system;
the effect of wind power injection power on the carbon flow rate distribution of each branch of the target power system is expressed as:
；
in the method, in the process of the invention,representing the carbon flow rate distribution for each leg of the target electrical power system,R us-B in order to balance the set-branch carbon flow correlation matrix,Houtputting a distribution matrix for the nodes;
the effect of wind power injection power on the carbon flow rate distribution of each load of the target power system is expressed as:
；
in the method, in the process of the invention,representing the carbon flow rate distribution for each load of the target electrical power system,R us-L for the unit-branch loss carbon flow correlation matrix,ξ N is thatNThe rank vector and all elements in the vector are 1,P L is thatM×NThe order-load distribution matrix is used to determine,P L the elements in the system are the connection relation between all the electric loads and the target electric power system and the active load quantity, P N Is thatNAn active flux matrix of the order diagonal array node;
the effect of wind power injection power on the carbon flow rate distribution of the system active loss is expressed as:
；
in the method, in the process of the invention,a carbon flow rate distribution representing the active loss of the system,P l for branches in a target power systemlIs not limited to the active loss of the (a).
8. A grid carbon flow calculation device taking system loss and wind power uncertainty into account, comprising:
the construction module is used for constructing a random wind speed sample set aiming at a target power system;
the first calculation analysis module is used for carrying out deterministic power flow calculation accounting for active loss based on each group of random wind speed samples in the random wind speed sample set, and determining a functional expression between wind power injection power considering wind speed uncertainty and the total carbon flow rate of the target power system;
the second calculation analysis module is used for calculating the related carbon flow index according to the obtained deterministic power flow calculation result; the related carbon flow indexes comprise a balance unit-node carbon flow incidence matrix, a balance unit-branch carbon flow incidence matrix, a unit-branch loss carbon flow incidence matrix, a node output distribution matrix and a path output distribution matrix;
the third calculation analysis module is used for carrying out carbon flow analysis based on the related carbon flow index; the carbon flow analysis includes calculating contributions of balancing units to active power flows of portions of the target power system and calculating carbon flow rate influence factors of balancing units to nodes and branches.
9. A grid carbon flow calculation device taking system loss and wind power uncertainty into account, comprising:
a memory for storing instructions; the instructions are used for realizing the power grid carbon flow calculation method taking system loss and wind power uncertainty into consideration according to any one of claims 1-7;
and the processor is used for executing the instructions in the memory.
10. A computer readable storage medium, wherein a computer program is stored on the computer readable storage medium, and when executed by a processor, the computer program implements the method for calculating the carbon flow of the power grid taking into account the system loss and wind power uncertainty as claimed in any one of claims 1 to 7.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202310770811.9A CN116502961A (en) | 2023-06-28 | 2023-06-28 | Power grid carbon flow calculation method and device considering system loss and wind power uncertainty |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202310770811.9A CN116502961A (en) | 2023-06-28 | 2023-06-28 | Power grid carbon flow calculation method and device considering system loss and wind power uncertainty |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN116502961A true CN116502961A (en) | 2023-07-28 |
Family
ID=87325234
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202310770811.9A Pending CN116502961A (en) | 2023-06-28 | 2023-06-28 | Power grid carbon flow calculation method and device considering system loss and wind power uncertainty |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN116502961A (en) |
Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN113886767A (en) * | 2021-08-30 | 2022-01-04 | 国网江苏省电力有限公司 | Real-time calculation method for carbon emission flow of power system |
| CN114707121A (en) * | 2022-04-29 | 2022-07-05 | 深圳供电局有限公司 | Carbon emission calculation method for distributed power supply random output power distribution network |
| CN115882461A (en) * | 2022-12-23 | 2023-03-31 | 国网辽宁省电力有限公司电力科学研究院 | Grid loss considered power system node marginal carbon potential calculation method and system |
| CN115907281A (en) * | 2022-09-23 | 2023-04-04 | 国网甘肃省电力公司 | Power system carbon emission flow analysis method considering photovoltaic uncertainty |
| CN116231657A (en) * | 2023-05-09 | 2023-06-06 | 国网浙江省电力有限公司 | Global carbon flow distributed determination method and device for transmission and distribution network |
| US20230198258A1 (en) * | 2021-12-16 | 2023-06-22 | Singularity Energy, Inc. | Apparatus and method for optimizing carbon emissions in a power grid |
-
2023
- 2023-06-28 CN CN202310770811.9A patent/CN116502961A/en active Pending
Patent Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN113886767A (en) * | 2021-08-30 | 2022-01-04 | 国网江苏省电力有限公司 | Real-time calculation method for carbon emission flow of power system |
| US20230198258A1 (en) * | 2021-12-16 | 2023-06-22 | Singularity Energy, Inc. | Apparatus and method for optimizing carbon emissions in a power grid |
| CN114707121A (en) * | 2022-04-29 | 2022-07-05 | 深圳供电局有限公司 | Carbon emission calculation method for distributed power supply random output power distribution network |
| CN115907281A (en) * | 2022-09-23 | 2023-04-04 | 国网甘肃省电力公司 | Power system carbon emission flow analysis method considering photovoltaic uncertainty |
| CN115882461A (en) * | 2022-12-23 | 2023-03-31 | 国网辽宁省电力有限公司电力科学研究院 | Grid loss considered power system node marginal carbon potential calculation method and system |
| CN116231657A (en) * | 2023-05-09 | 2023-06-06 | 国网浙江省电力有限公司 | Global carbon flow distributed determination method and device for transmission and distribution network |
Non-Patent Citations (2)
| Title |
|---|
| 袁书林 等: "基于电力系统碳排放流理论的碳排放分摊模型研究", 现代电力, vol. 31, no. 6, pages 70 - 75 * |
| 马瑞 等: "考虑风电不确定性的电力系统碳排放流分析", 电力系统自动化, vol. 38, no. 17, pages 124 - 129 * |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN107133695A (en) | A kind of wind power forecasting method and system | |
| Wang et al. | Two-stage multi-objective unit commitment optimization under hybrid uncertainties | |
| Fang et al. | Deep reinforcement learning for scenario-based robust economic dispatch strategy in internet of energy | |
| CN111709569A (en) | Method and device for predicting and correcting output power of wind power plant | |
| CN105634018A (en) | Random response surface method and interior point method based wind-power-plant-contained random optimal power flow solving method | |
| CN114782217A (en) | Indirect carbon emission refined accounting method and system for power system | |
| CN106372440B (en) | A kind of adaptive robust state estimation method of the power distribution network of parallel computation and device | |
| CN113935588A (en) | Digital currency/carbon emission multithread computing network | |
| CN108957352A (en) | A kind of capacity life consumption calculation method based on state-of-charge | |
| CN109149559A (en) | A kind of Demand-side interconnection reliability estimation method and system | |
| CN116502961A (en) | Power grid carbon flow calculation method and device considering system loss and wind power uncertainty | |
| Wiesner et al. | A testbed for carbon-aware applications and systems | |
| CN111786391A (en) | Probabilistic power flow analysis method and system considering energy uncertainty | |
| CN110059840B (en) | Method and system for selecting address of battery energy storage system in receiving-end power grid | |
| CN113627655B (en) | Method and device for simulating and predicting pre-disaster fault scene of power distribution network | |
| CN116418000A (en) | Bayesian state estimation method for unobservable distribution network based on embedded tide neural network | |
| Deshmukh et al. | Application of genetic algorithm to neural network model for estimation of wind power potential | |
| CN116561502A (en) | Interpolation method, device and equipment for wind speed data of wind power plant and storage medium | |
| CN115296345A (en) | Start-stop-output-standby combined optimization method and device for gas generator set | |
| CN103326396A (en) | Method for testing wind power service capacity upper limit value of power system | |
| CN115660187B (en) | Low-carbon town ground source heat pump capacity optimization configuration method | |
| CN116050306B (en) | Power frequency power grid reliability assessment method and system considering offshore wind power frequency division access | |
| CN113138899B (en) | Server fan vibration management and control method, system, terminal and storage medium | |
| Sun et al. | Credible Capacity Evaluation Method for Wind Farms Considering Output Correlation | |
| CN117833351A (en) | New energy power system operation risk online assessment method and system |
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 |