Multi-node grid-connected complex tidal current electric energy settlement method and system for electric power system

Abstract

The invention discloses an electric power system multi-node grid-connected complex tidal current electric energy settlement method and system, and the method comprises the steps: obtaining the forward electric energy, reverse electric energy and power data of all electric energy meters, calculating the self-generating capacity, total load electric energy and net generating capacity of a new energy power generation main body, and achieving the virtual meter combination metering of the new energy power generation main body; constructing a distribution factor matrix for the nodes of the power grid side electric energy meter and the nodes of the new energy main body virtual combined table to represent the electric energy distribution between the bus nodes; calculating an electric energy contribution ratio of the power supply node to the power utilization node by using the distribution factor matrix, and calculating an electric energy distribution matrix according to the contribution ratio, the injection electric energy, the power supply electric energy and the power utilization electric energy; and according to the electric energy distribution matrix, for different transaction scenes of settlement between the new energy and the power grid and settlement between the new energy main bodies, performing differential settlement by adopting online and offline electricity prices and contract electricity prices. According to the method, the problems of unknown power flow path and unknown electric energy distribution during multi-main-body grid connection are solved, and the transaction fairness is ensured.

Description

Technical Field

[0001] This invention relates to the field of smart grids, specifically to a method and system for calculating complex power flow in a multi-node grid-connected power system. Background Technology

[0002] With the large-scale construction of new power systems and the increasing openness of the electricity market, accurate electricity settlement is a crucial link in supporting the transformation of the national energy structure and maintaining the order of the electricity market. However, due to natural constraints, new energy power generation exhibits significant intermittency, strong volatility, and unpredictable output, often manifesting as dynamic problems such as active power oscillations and short-term voltage amplitude and system frequency fluctuations in actual grid operation. In some cases, the power source of new energy power generation may be under maintenance or engaged in self-consumption, thus becoming a load. At different times, it may consume or generate electricity. This complexity and bidirectional nature of the power flow leads to secondary metering of self-generated and self-consumed electricity from new energy power plants at both the power source and load ends. Due to differences in grid and on-grid electricity prices, economic disputes arise between new energy power plants and the grid due to electricity settlement issues. Especially in scenarios where multiple new energy power generation entities adopt multi-point access and a shared bus topology for grid connection, the system power flow becomes even more complex, with the power transmission path and load consumption trajectory forming a complex network topology due to the energy interaction of multiple entities. If the traditional single-line energy metering method is still used for settling electricity consumption between grid connections, it will be impossible to accurately represent the power flow distribution characteristics during multi-path energy transmission, leading to a significant deviation between the settlement results and the actual amount of electricity transmitted. The limitations of this metering method become increasingly apparent in complex grid connection scenarios, causing serious economic losses to power generation companies. Summary of the Invention

[0003] The technical problem this invention aims to solve is to provide a method and system for settling complex power flow in multi-node grid-connected power systems, addressing the aforementioned issues in existing technologies. This method resolves the problems of unclear power flow paths and unclear power allocation when multiple entities are connected to the grid, ensuring fair transactions.

[0004] To solve the above-mentioned technical problems, the technical solution adopted by the present invention is as follows: A method for calculating power flow in complex multi-node grid-connected power systems includes the following steps: Acquire the forward energy, reverse energy, and power data of all energy meters within a specified time interval; The power flow direction of the corresponding node is determined based on the power data of the electricity meter corresponding to each new energy power generation entity. Based on the power flow direction of the corresponding node of each new energy power generation entity and the forward and reverse electricity of the corresponding electricity meter, the self-generated electricity, total load electricity and net power generation of each new energy power generation entity are calculated to realize virtual combined metering of each new energy entity. A distribution factor matrix is ​​constructed to characterize the power distribution among bus nodes for the nodes of the grid-side electricity meter and the nodes of the virtual combined meter of different new energy entities. One end of the grid-side electricity meter is a grid-side node and the other end is a bus node. One end of the virtual combined meter of new energy entities is a new energy entity node and the other end is a bus node. The grid-side node includes grid-side power supply nodes and grid-side power consumption nodes. The new energy entity node includes new energy entity power supply nodes and new energy entity power consumption nodes. The power contribution ratio of the power supply node to the power consumption node is calculated using the allocation factor matrix. Based on the power contribution ratio of the power supply node to the power consumption node, the injected power, the supplied power, and the consumed power, the power allocation matrix is ​​calculated. The power supply node includes grid-side power supply nodes and new energy main power supply nodes. The power consumption node includes grid-side power consumption nodes and new energy main power consumption nodes. The settlement of electricity between the power supply nodes of the new energy main body and the power consumption nodes of the grid side is carried out according to the power distribution matrix and the on-grid electricity price of the new energy power generation main body. The settlement of electricity between the power supply nodes of the grid side and the power consumption nodes of the new energy side is carried out according to the power distribution matrix and the off-grid electricity price of the grid. The settlement of electricity between the power supply nodes of the new energy main body and the power consumption nodes of the new energy main body is carried out according to the power distribution matrix and the contract electricity price.

[0005] Furthermore, when determining the power flow direction of the corresponding node based on the power data of the electricity meter corresponding to each new energy power generation entity, the specific steps include: If the electricity meter n q,c Power data If it is less than 0, then the electricity meter n q,c The power flow direction of the corresponding node is from the corresponding node to the bus node; If the electricity meter n q,c Power data If the value is greater than 0, then the electricity meter n q,c The power flow direction of the corresponding node is from the bus node to the corresponding node.

[0006] Furthermore, when calculating the self-generated power, total load, and net power generation of each new energy power generation entity based on the power flow direction of the corresponding node and the corresponding forward and reverse power of the electricity meter, the following steps are included: If the electricity meter nq,c The power flow direction of the corresponding node is from the corresponding node to the bus node. The status tag of the corresponding node is set to the first value, and the self-generated energy of the corresponding node is equal to that of the electricity meter. n q,c The reverse electrical energy increment; If the electricity meter n q,c The power flow direction of the corresponding node is from the bus node to the corresponding node. The status label of the corresponding node is set to the second value, and the load power of the corresponding node is equal to that of the power meter. n q,c The positive increase in electrical energy; The self-generated energy of all nodes with the first state label corresponding to the current new energy power generation entity is accumulated to obtain the self-generated energy of the current new energy power generation entity. The total load power of the current new energy power generation entity is obtained by summing the load power of all nodes whose status label is the second value. The net power generation of the current renewable energy power generation entity is obtained by subtracting the total load power from the self-generated power of the current renewable energy power generation entity.

[0007] Furthermore, before constructing the allocation factor matrix for the nodes of the grid-side electricity meters and the nodes of the virtual combined meters of different new energy entities, the process also includes the step of dividing the power supply nodes and the power consumption nodes, specifically including: If the power data of the power meter corresponding to the current power grid node is less than 0, then the current power grid node is a power supply node on the power grid side; if the power data of the power meter corresponding to the current power grid node is greater than 0, then the current power grid node is a power consumption node on the power grid side. If the net power generation of the new energy power generation entity corresponding to the current new energy main node is greater than 0, then the current new energy main node is a power supply node of the new energy main node; if the net power generation of the new energy power generation entity corresponding to the current new energy main node is less than 0, then the current new energy main node is a power consumption node of the new energy main node.

[0008] Furthermore, the mathematical expression for the allocation factor matrix is ​​as follows:

[0009] in, I For size M × M The identity matrix, A For size M × M The allocation factor matrix, whose elements A ( x , y Characterizing busbar nodes x For bus nodes yThe proportion of electricity contribution, x , y =1,2,…, M , and They represent vectors respectively F and E Z Construct a diagonal matrix by taking the reciprocal of each element. It is a matrix The transpose of the matrix, the matrix It is a tie-line correlation matrix constructed based on the main electrical wiring diagram and the positive and negative power ratings of the tie-line energy meters. B A matrix in which all non-zero elements are 1. It is a tie-line correlation matrix constructed based on the main electrical wiring diagram and the positive and negative power ratings of the tie-line energy meters. B A matrix in which all non-zero elements are -1, the connection line incidence matrix. B medium elements B d,m Indicates the first d The connecting line and the first m The connection relationship of each node. B d,m =1 represents the first d The starting node of the connecting line is m ; B d,m =-1 represents the first d The termination node of the connecting line is m ; B d,m =0 represents the first d Connecting lines and nodes m No connection; vector F The electrical energy vector of the tie line is expressed mathematically as follows: F =[ F 1, F 2,…, F d ,…, F D ] T in, d =1,2,…, D , No. d Transmission of electrical energy by the connecting line F d The formula is as follows:

[0010] in, Indicates the first d The positive energy increment of the electricity meter on the connecting line. and They are the first d The positive electrical energy of the electricity meter on the connecting line at the end and start times of a specified time interval. Indicates the first d The reverse energy increment of the electricity meter on the connecting line, and They are the first d The reverse energy of the electricity meter on the connecting line at the end and start times of a specified time interval. It is the first d The power data of the electricity meter on the connecting line at the end of the specified time interval; vector E Z The vector representing the injected electrical energy at the bus node is expressed mathematically as follows: E Z =[ ] T Among them, the bus node mid m The mathematical expression for the injected electrical energy is as follows:

[0011] In the formula, S m mid, the bus node m The collection of upstream connecting lines, For the first d m The power transmission of the connecting line To connect with the bus node mid m Connected power supply nodes v The value of the supplied electrical energy is the reverse electrical energy increment of the grid-side electricity meter corresponding to the grid-side electricity meter if the power supply node is a grid-side power supply node, and the value of the supplied electrical energy is the absolute value of the net power generation of the new energy entity corresponding to the new energy entity power supply node if the power supply node is a new energy entity power supply node.

[0012] Furthermore, the mathematical expression for the power distribution matrix is ​​as follows:

[0013] in, G V-diag Represented by vector G V Constructed diagonal matrix, vector G V This represents the vector of electrical energy supplied by the power supply node. G V The first vector Rg Each element represents the electrical energy supplied by the power grid-side power supply node. Q g Each element represents the electrical energy supplied by the main power supply node of the new energy source; A sub This represents the allocation factor submatrix between the downstream bus nodes of the power supply node and the upstream bus nodes of the power consumption node. A sub ( v , u ) represents the first element in the matrix. v Okay, number u The elements of the column are expressed mathematically as follows: A sub ( v , u )= A ( G ind ( v ), L ind ( u )) In the formula, A ( G ind ( v ), L ind ( u )) represents the allocation factor matrix G ind ( v )OK L ind ( u ) elements of column, G ind ( v ) represents the downstream bus node vector of the power supply node. G ind The Middle v Each bus node L ind ( u ) represents the upstream bus node vector of the power consumption node. L ind The Middle u Each bus node G ind The first vector R g Each element represents the bus node corresponding to the power supply node on the grid side. Q g Each element represents the bus node corresponding to the main power supply node of the new energy source. L ind The first vector Rl Each element represents the bus node corresponding to the power consumption node on the grid side. Q g Each element represents the bus node corresponding to the main electricity consumption node of new energy sources; , Represents a vector Copy in the vertical direction V This forms a new matrix and vector. Power consumption node power consumption vector L U transpose, L U The first vector R l Each element represents the electrical energy consumed by the power consumption nodes on the grid side. Q g Each element represents the electrical energy consumed by the main electricity consumption nodes of new energy sources; , Indicates to E Z The new vector obtained after vector indexing E Z ( G ind Copy horizontally V The resulting new matrix, where V represents the total number of power supply nodes. Matrix representation of a matrix E Z-sub The matrix obtained by taking the reciprocal of each element This represents the Hadamard product of a matrix. Furthermore, when settling electricity accounts between renewable energy power supply nodes and grid-side power consumption nodes based on the electricity distribution matrix and the on-grid tariff of renewable energy power generation, the mathematical expression is as follows:

[0014] in, Y up ( a ) indicates the main power supply node of new energy. a Revenue from electricity generation and grid connection for power consumption nodes on the grid side λ up The on-grid electricity price for electricity generated by the main source of renewable energy power generation. For the power distribution matrix, the first Line number Column elements represent the set of power supply nodes. S sup The Middle Power supply nodes For the set of power consumption nodes Suse The Middle individual electricity consumption nodes The amount of electrical energy distributed. S sup Before the set R g Each element represents a power supply node on the grid side. Q g Each element represents a new energy source power supply node. S use Before the set R l Each element represents a power consumption node on the grid side. Q g Each element represents a new energy main electricity consumption node.

[0015] Furthermore, when settling electricity accounts between grid-side power supply nodes and renewable energy-side power consumption nodes based on the electricity distribution matrix and the grid's off-grid electricity price, the data expression is as follows:

[0016] in, K down (b) Indicates the main electricity consumption node of new energy sources. b The cost of electricity drawn from power supply nodes on the grid side. λ down The off-grid electricity price of the power grid. For the power distribution matrix, the first Line number Column elements represent the set of power supply nodes. S sup The Middle Power supply nodes For the set of power consumption nodes S use The Middle individual electricity consumption nodes The amount of electrical energy allocated.

[0017] Furthermore, when settling electricity accounts between renewable energy power supply nodes and renewable energy consumption nodes based on the electricity allocation matrix and contracted electricity prices, the data expression is as follows:

[0018] in, Y PPA ( a ) indicates the main power supply node of new energy source. a Power supply to new energy main power consumption nodes b The revenue from electricity generation, K PPA ( b) represents the main electricity consumption node of new energy sources. b Power is drawn from the main power supply node of new energy sources. a Electricity costs, λ PPA This refers to the long-term contract electricity price signed between renewable energy power suppliers and renewable energy power consumers. For the power distribution matrix, the first Line number Column elements represent the set of power supply nodes. S sup The Middle Power supply nodes For the set of power consumption nodes S use The Middle individual electricity consumption nodes The amount of electrical energy allocated.

[0019] The present invention also proposes a power system for calculating complex power flow in a multi-node grid-connected power system, comprising a processor and a computer-readable storage medium, wherein the computer-readable storage medium stores a computer program, and the computer program is executed by the processor to implement the steps of the power system multi-node grid-connected complex power flow calculation method.

[0020] Compared with the prior art, the advantages of the present invention are as follows: This invention collects forward and reverse electrical energy and real-time power data, combines power to determine power flow direction, and accurately identifies the operating status of each renewable energy source at any given time. Simultaneously, by calculating the net power generation of each source, it achieves virtual combined metering for a single renewable energy source, integrating physically dispersed metering points into a logically unified node.

[0021] This invention establishes an allocation factor matrix to quantify the proportion of power contribution from each power supply node to the power consumption node. This, in turn, generates a power allocation matrix, clarifying the distribution of power from the power supply node to the power consumption node, thus solving the problems of unclear power flow paths and unclear power allocation in multi-entity grid interconnection.

[0022] Based on the power distribution matrix, this invention adopts differentiated settlement methods, namely grid connection price and contract price, for different transaction scenarios such as settlement between new energy sources and the power grid and settlement between new energy entities, to ensure fair transactions. Attached Figure Description

[0023] Figure 1 This is a flowchart of a method according to an embodiment of the present invention.

[0024] Figure 2 This is a schematic diagram of the flow of electrical energy. Detailed Implementation

[0025] The present invention will be further described below with reference to the accompanying drawings and specific preferred embodiments, but this does not limit the scope of protection of the present invention.

[0026] Due to the intermittent and fluctuating nature of renewable energy generation, individual renewable energy power plants may engage in self-consumption, meaning that the power source of a renewable energy power plant may supply power to its load end via a busbar. However, due to the geographically dispersed nature of the power generation and load units, combined metering is not possible. Existing settlement schemes use the meter readings of a single line to accumulate and settle electricity, failing to independently and effectively identify the flow path of self-generated and self-consumed electricity. This results in renewable energy power plants' self-generated and self-consumed electricity being measured as grid-connected electricity at the power source, while the electricity consumed at the load end is repeatedly measured as grid-connected electricity. Because of the price difference between grid-connected and grid-connected electricity, self-generated and self-consumed electricity undergoes two settlement processes with different pricing mechanisms, ultimately causing economic losses for the renewable energy power plant.

[0027] In multi-entity grid-connected renewable energy distribution areas, the electricity generated by renewable energy generators first flows into a common bus, and then is distributed by the bus. Part of this distribution may supply its own load, part may be directly consumed by other electricity consumers within the distribution area, and the remaining part will be transmitted to the upper-level grid, resulting in multiple distribution paths. Because electricity converges at the bus, the electricity from different entities mixes there, and traditional electricity meters can only record the electricity data of a single node. When settling electricity accounts, for the generator side, the electricity meter can only count the total amount of electricity injected into the bus by the generator unit; for tie lines, it can only count the electricity transmitted between two buses, unable to trace the actual source and destination of the electricity, and unable to distinguish the actual flow path of the electricity from the power supply and consumption entities in the network; for the consumption side, the electricity meter can only count the total amount of electricity taken from the bus by the consumption entity, unable to know whether this electricity comes from renewable energy generation within the distribution area or from the upper-level grid. This causes actual errors in electricity settlement.

[0028] In multi-new energy power plant grid-connected systems sharing a common busbar, a complex energy interaction network is formed among the various power generation entities, leading to direct power supply and consumption behaviors that are not entirely dependent on the power grid. The current settlement system only uses inter-station electricity meters to measure transactions between new energy power plants and the grid company. For the mutual power supply and consumption between different power plant entities via the busbar, there is a lack of accurate metering methods based on power flow paths, and no clear settlement rule system has been established. This dual deficiency in metering and rules leads to disputes over the ownership of supplied and consumed energy and cost allocation, seriously affecting the fairness and order of the electricity trading market.

[0029] To address the aforementioned issues, this embodiment proposes a method for settling complex power flow in multi-node grid-connected power systems. Based on the current status of power settlement, this research on a method for settling complex power flow in multi-node grid-connected new energy systems is of great significance for ensuring fairness and impartiality in power settlement, effectively safeguarding the interests of both power supply and demand sides, and promoting the construction and development of new power systems.

[0030] like Figure 1 As shown, the method includes the following steps: S1) Based on the main wiring diagram and the installation location of the electricity meters, classify the electricity meters into new energy side electricity meters, grid side electricity meters and tie line electricity meters, and obtain the forward energy, reverse energy and power data of all electricity meters in the specified time interval; S2) Determine the power flow direction of the corresponding node based on the power data of the electricity meter corresponding to each new energy power generation entity, and then accumulate the base code of the electricity meter of the collection line of the new energy entity. Specifically, calculate the self-generated power, total load power and net power generation of each new energy power generation entity based on the power flow direction of the corresponding node of each new energy power generation entity and the forward and reverse power of the corresponding electricity meter, so as to realize the virtual combined metering of each new energy entity. S3) The system nodes are divided into grid-side, new energy main body, and bus nodes according to their attributes. An allocation factor matrix is ​​constructed for the nodes of the grid-side electricity meter and the nodes of the virtual combined meter of different new energy main bodies to represent the power distribution among the bus nodes. One end of the grid-side electricity meter is a grid-side node and the other end is a bus node. One end of the virtual combined meter of the new energy main body is a new energy main body node and the other end is a bus node. The grid-side node includes grid-side power supply node and grid-side power consumption node. The new energy main body node includes new energy main body power supply node and new energy main body power consumption node. S4) Calculate the power contribution ratio of the power supply node to the power consumption node using the allocation factor matrix. Calculate the power allocation matrix based on the power contribution ratio of the power supply node to the power consumption node, the injected power, the supplied power, and the consumed power. The power supply node includes grid-side power supply nodes and new energy main power supply nodes. The power consumption node includes grid-side power consumption nodes and new energy main power consumption nodes. S5) For different trading scenarios involving settlement between renewable energy sources and the power grid, and settlement between renewable energy entities, differentiated settlement will be implemented using both grid connection and contract electricity prices. Specifically, this includes: The settlement of electricity between the power supply nodes of the new energy main body and the power consumption nodes of the grid side is carried out according to the power distribution matrix and the on-grid electricity price of the new energy power generation main body. The settlement of electricity between the power supply nodes of the grid side and the power consumption nodes of the new energy side is carried out according to the power distribution matrix and the off-grid electricity price of the grid. The settlement of electricity between the power supply nodes of the new energy main body and the power consumption nodes of the new energy main body is carried out according to the power distribution matrix and the contract electricity price.

[0031] Through the above steps, the method in this embodiment achieves virtual combined metering for a single renewable energy entity. Based on the contribution ratio of the power supply node to the power consumption node, the injected energy, the supplied energy, and the consumed energy, an energy distribution matrix is ​​calculated to clarify the distribution of energy from the power supply node to the power consumption node, thus resolving the problems of unclear power flow paths and unclear energy distribution when multiple entities are connected to the grid. During settlement, differentiated settlement is implemented using both grid connection prices and contract prices for different transaction scenarios, including settlement between renewable energy and the grid, and settlement between renewable energy entities, to ensure fair transactions.

[0032] The following is a detailed explanation of each step.

[0033] In step S1 of this embodiment, for the new energy power generation common bus grid connection system, one end of the new energy side energy meter is connected to the new energy side node, and the other end is connected to the bus node; one end of the grid side energy meter is connected to the grid side node, and the other end is connected to the bus node; the tie line energy meter is connected to two bus nodes.

[0034] Assuming the main wiring diagram for grid connection of new energy power generation contains Q One new energy power generation entity, N An electricity meter that only has bidirectional metering function, assuming that one of the meters is on the grid side. R One. Within the time interval 0~ T Within the transformer area, assuming the power flow direction remains unchanged, read... N Only the forward and reverse electrical energy at the start and end times of the electricity meter, and N The power at the end of the time interval of the electricity meter is obtained. N Only the corresponding electrical energy and power data of the electricity meter, including the forward electrical energy vector at the start time, the forward electrical energy vector at the end time, the reverse electrical energy vector at the start time, the forward electrical energy vector at the end time, and the power vector at the end time, are denoted as follows: W po-s , W po-e , W ne-s , W ne-e and P e The subscript po represents forward electrical energy, ne represents reverse electrical energy, s represents the start time, and e represents the end time. Each vector has a total of N The nth element, the nth element of each vector n element W po-s,n , W po-e,n , W ne-s,n , W ne-e,n and P e,n, respectively representing the n Only the forward energy at the start time, the forward energy at the end time, the reverse energy at the start time, the forward energy at the end time, and the power at the end time of the electricity meter. n =1,2,…, N .

[0035] In this embodiment, step S2 is used to perform virtual metering for a single new energy power generation entity, prioritizing that self-generated and self-consumed electricity is not included in the settlement.

[0036] When determining the power flow direction of a corresponding node based on the power data of the electricity meter corresponding to each new energy power generation entity, for the first... q One new energy power generation entity ( q =1,2,…, Q ), determine the common busbar collector line to which it belongs C q Only electricity meters, denoted as a set. S q ={ n q,1 , n q,2 ,…, n q,c ,…, n q,Cq},in n q,c express S q The Middle c One element, c =1,2,…, C q For sets S q Each electricity meter in the system, since the power flow of the distribution lines in the area did not change during the time period, calculates the power vector at the end of the specified time interval. P e Determine the set S q The power flow direction of the line where the electricity meter is located. Specifically, this includes: If the electricity meter n q,c Power data Less than 0 ( <0 indicates that the electricity meter n q,c The power of the line flows from the renewable energy node to the bus node, i.e., the electricity meter. n q,c If the power flow direction of the corresponding node is from the corresponding node to the bus node, then the new energy side node generates electricity. If the electricity meter n q,c Power data Greater than 0 ( >0), indicating that the electricity meter n q,c The power of the line flows from the bus node to the new energy side node, i.e., the electricity meter. n q,c If the power flow direction of the corresponding node is from the bus node to the corresponding node, then the new energy side node consumes electrical energy.

[0037] Correspondingly, when calculating the self-generated power, total load, and net power generation of each new energy power generation entity based on the power flow direction of the corresponding node and the corresponding forward and reverse power of the electricity meter, the following steps are included: If the electricity meter n q,c The power flow direction for the corresponding node is from the corresponding node to the bus node; set the status label for the corresponding node. δ q,c The first value ( δ q,c =1), and the self-generated energy of the corresponding node is equal to that of the electricity meter. n q,c The reverse electrical energy increment is = - ; If the electricity meter n q,c The power flow direction for the corresponding node is from the bus node to the corresponding node, and the status label for the corresponding node is set. δ q,c The second value ( δ q,c =0), and the load energy of the corresponding node is equal to the energy meter reading. n q,c The positive increase in electrical energy is = - ; The self-generated energy of the current renewable energy power generation entity is obtained by summing the self-generated energy of all nodes whose state label is the first value. The mathematical expression is as follows: ; The total load power of the current renewable energy generation entity is obtained by summing the load power of all nodes with the second state label corresponding to the current renewable energy generation entity. The mathematical expression is as follows: ; Subtracting the total load from the current self-generated electricity of the renewable energy power generation entity yields the net electricity generation of the current renewable energy power generation entity, expressed mathematically as follows: Enet-q = E gen-q - E load-q .

[0038] Step S2 is used to virtually combine the power collection lines of different new energy sources, and then the net power generation of each new energy source is used as the basis for the calculation. E net-q Settlement will be conducted with priority given to ensuring that new energy power generation entities can use their own generated electricity for their own consumption and that surplus electricity can be fed into the grid.

[0039] by Figure 2 For example, there are 3 new energy power generation sources in the diagram, namely main body 1 ( v 2) Main Body 2 ( u 2) and main body 3 ( v 3) Calculate the self-generated electricity, total load power and net power generation of each new energy power generation entity through step S2. The results are shown in Table 1.

[0040] Table 1. Total self-generated energy, total load energy, and net power generation of new energy power generation entities

[0041] In this embodiment, node classification, matrix construction, and allocation factor matrix calculation are achieved through step S3.

[0042] To clearly delineate power transmission paths and clarify the power interaction relationships between each meter and other nodes in a shared busbar scenario, avoiding issues such as unclear node affiliation and chaotic power transmission caused by the shared busbar, the first and last nodes of grid-side power meters and the first and last nodes of virtual combined meters for different renewable energy entities are recoded based on the actual main wiring diagram of renewable energy generation grid connection. One end of a grid-side power meter is a grid-side node, and the other end is a busbar node; one end of a virtual combined meter for renewable energy entities is a renewable energy entity node, and the other end is a busbar node; tie-line power meters connect two busbar nodes. Based on the aforementioned assumptions, the grid-side nodes in the network have a total of... R One, main nodes of new energy Q There are [number], assuming there are a total of [number] other necessary bus nodes. M These are denoted as mid1, mid2, ..., mid M .

[0043] Grid side R Let the set contain _ _ nodes. S grid ={grid1,grid2,…,grid R}, where grid r express S grid The Middle r One element, r =1,2,…,R New energy main nodes Q One, corresponding to Q There are several new energy power generation entities, and the nodes of these new energy entities are denoted as a set. S new ={new1,new2,…,new Q}, where new q express S new The Middle q One element, q =1,2,…, Q .

[0044] In step S3, before constructing the allocation factor matrix for the nodes of the grid-side electricity meters and the nodes of the virtual combined meters of different new energy entities, the step of dividing the power supply nodes and the power consumption nodes is also included, specifically: Based on the power direction of the grid-side energy meters at the grid-side nodes, the grid-side nodes are divided into grid-side power supply nodes and grid-side power consumption nodes. Specifically, if the current grid-side node is grid... r The corresponding power data of the grid-side electricity meter Less than 0 ( If <0, then the current grid-side node grid r For the power supply node on the grid side, if the current grid-side node is grid r The corresponding power data of the grid-side electricity meter Greater than ( If the value is greater than 0, then the current grid-side node is a grid-side power-consuming node. Assume there are grid-side power supply nodes. R g There are [number] power consumption nodes on the grid side. R l indivual, R g and R l satisfy R g + R l = R The set of power supply nodes on the grid side includes all nodes supplying power from the grid side, denoted as... G grid ={grid i | i ∈{ i 1, i 2,…, The set of power-consuming nodes on the grid side includes all power-consuming nodes on the grid side, denoted as... L grid ={grid j | j ∈{j 1, j 2,…, }}; Based on the net power generation of the corresponding renewable energy power generation entity, renewable energy power generation entities are divided into renewable energy power supply entities and renewable energy power consumption entities. Specifically, if the current renewable energy power generation entity is new... q The net power generation of the corresponding new energy power generation entities E net-q Greater than 0 ( If >0), then the current new energy main node new q For the main power supply node of new energy, if the current main power supply node is new q The net power generation of the corresponding new energy power generation entities E net-q Less than 0 ( If <0), then the current new energy main node new q This refers to the main electricity consumption nodes of new energy sources. Assume there are a total of [number missing] main power supply nodes of new energy sources. Q g There are [number] new energy main electricity consumption nodes. Q l indivual, Q g and Q l satisfy Q g + Q l =Q The set of power supply nodes for the main renewable energy source includes all nodes supplying power to the main renewable energy source, denoted as... G new ={new a | a ∈{ a 1, a 2,… The set of electricity-consuming nodes for new energy entities includes all electricity-consuming nodes of the new energy entities, denoted as... L new ={new b | b ∈{ b 1, b 2,… }}.

[0045] The system has a total of V Each power supply node U Each electricity consumption node, then V = R g + Q g , U =R l + Q l The power supply node is represented as S sup ={grid i new a | i ∈{ i 1, i 2,…, }, a ∈{ a 1, a 2,… } = {sup1,sup2,…,sup V}, the electricity consumption node is represented as S use ={ grid j new b | j ∈{ j 1, j 2,…, }, b ∈{ b 1, b 2,… } = {use1, use2, ..., use U}

[0046] Construct the power supply vector of the power supply node and the power consumption vector of the power consumption node: G V = [ ] T = [ ] T L U = [ ] T = [ ] T In the formula, G V The first vector R g Each element represents the electrical energy supplied by the power grid-side power supply node. Q g Each element represents the electrical energy supplied by the main power supply node of the new energy source. L U The first vector R l Each element represents the electrical energy consumed by the power consumption nodes on the grid side. Q gEach element represents the electrical energy consumed by the main electricity consumption nodes of new energy sources. G V and L U Each element is calculated using the following formula: ,

[0047] ,

[0048] That is, if the power supply node is a grid-side power supply node, the value of the supplied electrical energy is the reverse electrical energy increment of the grid-side electricity meter corresponding to the grid-side power supply node; if the power supply node is a renewable energy main body power supply node, the value of the supplied electrical energy is the absolute value of the net power generation of the renewable energy main body corresponding to the renewable energy main body power supply node; if the power consumption node is a grid-side power consumption node, the value of the consumed electrical energy is the forward electrical energy increment of the grid-side electricity meter corresponding to the grid-side power consumption node; if the power consumption node is a renewable energy main body power consumption node, the value of the supplied electrical energy is the absolute value of the net power generation of the renewable energy main body corresponding to the renewable energy main body power consumption node.

[0049] by Figure 2 For example, there are 3 power supply nodes in the diagram ( v 1. v 2. v 3) and 2 power consumption nodes ( u 1. u 2) The corresponding power supply vectors of the power supply nodes and the power consumption vectors of the power consumption nodes are as follows: G V =[82.4,54.5,9.7] T , L U =[103.7,42.9] T To clarify the power distribution path, it is necessary to construct the downstream bus node vector of the power supply node. G ind and the upstream bus node vector of the power consumption node L ind The mathematical expression is as follows: G ind =[ ] L ind =[ ] in, G ind The first vector R g Each element represents the bus node corresponding to the power supply node on the grid side.Q g Each element represents the bus node corresponding to the main power supply node of the new energy source. L ind The first vector R l Each element represents the bus node corresponding to the power consumption node on the grid side. Q g Each element represents a bus node corresponding to a main power consumption node of new energy sources.

[0050] by Figure 2 For example, power supply node v 1 and v The bus node corresponding to 2 is mid1, the power supply node. v The bus node corresponding to 3 is mid3, the power consumption node. u 1. u The bus nodes corresponding to 2 are mid2 and mid4, therefore, the downstream bus node vector of the power supply node is constructed. G ind and the upstream bus node vector of the power consumption node L ind for: G ind =[1,1,3], L ind =[2,4].

[0051] Assuming there are multiple [networks] within the network D Connecting lines, constructing the power vector of the connecting lines. F as follows: F =[ F 1, F 2,…, F d ,…, F D ] T in, d =1,2,…, D , No. d Transmission of electrical energy by the connecting line F d The cumulative calculation can be based on the meter readings on the connection line, using the following formula:

[0052] in, Indicates the first d The positive energy increment of the electricity meter on the connecting line. and They are the first d The positive electrical energy of the electricity meter on the connecting line at the end and start times of a specified time interval. Indicates the first d The reverse energy increment of the electricity meter on the connecting line, and They are the first d The reverse energy of the electricity meter on the connecting line at the end and start times of a specified time interval. It is the first d The power data of the electricity meter on the connecting line at the end of the specified time interval.

[0053] by Figure 2 For example, if there are tie-line energy meters between bus nodes mid1 and mid2, mid1 and mid3, mid1 and mid4, mid2 and mid3, and mid4 and mid3, then the tie-line energy vector is... F for: F =[105.3,20.2,11.4,1.6,22.7] T .

[0054] Based on the main electrical wiring diagram and the positive and negative power ratings of the tie line energy meters, the dimensions are constructed as follows: D × M Connection matrix of the connection lines B ,element B d,m Indicates the first d The connecting line and the first m The connection relationship of each node. B d,m =1 represents the first d The starting node of the connecting line is m ; B d,m =-1 represents the first d The termination node of the connecting line is m ; B d,m =0 represents the first d Connecting lines and nodes m No connections. (Identity matrix) B A matrix can be decomposed into a matrix whose non-zero elements are all -1. B u The matrix has all non-zero elements as 1. B d matrix, B d Corresponding to the starting node of the route. B u The corresponding line termination node satisfies B = B d + B u .

[0055] by Figure 2 For example, the connection matrix of the connecting lines B as follows:

[0056] Based on the conservation principle that the input and output electrical energy of a bus node are equal, the injected electrical energy at each node in the main electrical wiring diagram and the transmitted electrical energy along the lines can be calculated. This allows for the construction of a vector representing the injected electrical energy at each bus node. E Z The mathematical expression is as follows: E Z =[ ] T Bus node injected power vector E Z Length is M The elements represent the injected electrical energy at each bus node. Since the injected electrical energy at a bus node equals the input electrical energy from the upstream tie line plus the supplied electrical energy, the mid value of each bus node can be calculated. m The mathematical expression for the injected electrical energy is as follows:

[0057] In the formula, S m mid, the bus node m The collection of upstream connecting lines, For the first d m The power transmission of the connecting line To connect with the bus node mid m Connected power supply nodes v The power supply energy, if the power supply node is the first i For each grid-side power supply node, the value of the supplied electrical energy is the reverse energy increment of the grid-side energy meter corresponding to that node, i.e. ,and If the power supply node is the first a Each renewable energy power supply node supplies electricity at a given renewable energy source, and the value of the electricity supplied is the absolute value of the net power generation of the renewable energy source corresponding to that node. ,and .

[0058] by Figure 2 For example, the power vector injected into the bus node E Z for: E Z =[136.9,105.3,22.7,42.9] T .

[0059] Based on the aforementioned formula, an allocation factor matrix can be constructed to characterize the power distribution among bus nodes. The mathematical expression is as follows:

[0060] in, I For size M × M The identity matrix, A For size M × M The allocation factor matrix, whose elements A ( x , y Characterizing busbar nodes x For bus nodes y The proportion of electricity contribution, x , y =1,2,…, M , and They represent vectors respectively F and E Z Construct a diagonal matrix by taking the reciprocal of each element. It is a matrix The transpose of the matrix, the matrix This refers to the tie-line correlation matrix constructed based on the main electrical wiring diagram and the positive and negative power ratings of the tie-line energy meters. B A matrix in which all non-zero elements are 1. This refers to the tie-line correlation matrix constructed based on the main electrical wiring diagram and the positive and negative power ratings of the tie-line energy meters. B A matrix or vector whose non-zero elements are all -1. F This represents the aforementioned tie-line energy vector, the vector... E Z This represents the aforementioned power injection vector at the bus node.

[0061] by Figure 2 For example, the constructed allocation factor matrix is ​​as follows: .

[0062] In this embodiment, step S4 is used to process each power supply node sup v and electricity consumption nodes u , v =1,2,…, V , u =1,2,…, U Calculate the power distribution matrix.

[0063] In step S4, when calculating the energy contribution ratio of the power supply node to the power consumption node using the allocation factor matrix, the allocation factor sub-matrix between the downstream bus node of the power supply node and the upstream bus node of the power consumption node is constructed using the allocation factor matrix built in step S3. The mathematical expression is as follows: A sub ( v , u )= A ( G ind ( v ), L ind ( u )) In the formula, A sub ( v , u ) represents the th element in the allocation factor submatrix. v Okay, number u Column elements, A ( G ind ( v ), L ind ( u )) represents the allocation factor matrix G ind ( v )OK L ind ( u ) elements of column, G ind ( v ) represents the downstream bus node vector of the aforementioned power supply node. G ind In and power supply node sup v The corresponding number v Each bus node L ind ( u ) represents the upstream bus node vector of the aforementioned power consumption node. L ind In the middle and the electricity consumption node use u The corresponding number u Each bus node.

[0064] Based on this, the mathematical expression for calculating the energy distribution matrix according to the contribution ratio of the power supply node to the power consumption node, the injected energy, the supplied energy, and the consumed energy is as follows:

[0065] in, G V-diag Represented by vectorG V Constructed diagonal matrix, vector G V This represents the power vector supplied by the aforementioned power supply node; , Represents a vector Copy in the vertical direction V This forms a new matrix and vector. It is the aforementioned electrical energy vector of the electricity consumption node. L U transpose; , Indicates to E Z The new vector obtained after vector indexing E Z ( G ind Copy horizontally V The resulting new matrix, where V represents the total number of power supply nodes. Matrix representation of a matrix E Z-sub The matrix obtained by taking the reciprocal of each element This represents the Hadamard product of a matrix. Power distribution matrix DE No. v Line number u Column elements DE ( v , u ) indicates the first v The power supply node for the first u The power distribution at each power consumption node is used to clearly define the power supply and consumption path. Figure 2 For example, establish an energy distribution matrix. DE The power allocation from each power supply node to the power consumption node is shown in Table 2.

[0066] Table 2. Distribution of electrical energy between power supply nodes and power consumption nodes (unit: kWh)

[0067] In this embodiment, step S5 utilizes the power distribution matrix. DE Accurately settle the revenue and expenses of each power supply and consumption node.

[0068] For electricity settlement between renewable energy power supply nodes and grid-side power consumption nodes, due to the price difference between renewable energy sources and the grid, settlement is based on separate accounting for electricity from both sources: For new energy power supply nodes a That is, the set of power supply nodesS sup The Middle Power supply nodes When settling electricity accounts between renewable energy power supply nodes and grid-side power consumption nodes based on the electricity allocation matrix and the on-grid electricity price of renewable energy power generation entities, the mathematical expression is as follows:

[0069] in, Y up ( a ) indicates the main power supply node of new energy. a Revenue from electricity generation and grid connection for power consumption nodes on the grid side λ up The on-grid electricity price for electricity generated by the main source of renewable energy power generation. For the power distribution matrix, the first Line number Column elements represent the set of power supply nodes. S sup The Middle Power supply nodes For the set of power consumption nodes S use The Middle individual electricity consumption nodes The amount of electrical energy distributed, as can be seen from the preceding text, S sup Before the set R g Each element represents a power supply node on the grid side. Q g Each element represents a new energy source power supply node. S use Before the set R l Each element represents a power consumption node on the grid side. Q g Each element represents a new energy main electricity consumption node.

[0070] For new energy-side electricity consumption nodes b That is, the set of electricity consumption nodes S use The Middle individual electricity consumption nodes When settling electricity costs between grid-side power supply nodes and renewable energy-side power consumption nodes based on the electricity distribution matrix and the grid's off-grid electricity price, the data expression is as follows:

[0071] in, K down (b) Indicates the main electricity consumption node of new energy sources.b The cost of electricity drawn from power supply nodes on the grid side. λ down The off-grid electricity price of the power grid. For the power distribution matrix, the first Line number Column elements represent the set of power supply nodes. S sup The Middle Power supply nodes For the set of power consumption nodes S use The Middle individual electricity consumption nodes The amount of electrical energy allocated.

[0072] For new energy main power supply nodes a With new energy main electricity consumption nodes b When settling electricity payments between renewable energy power supply nodes and renewable energy consumption nodes based on the electricity allocation matrix and contract electricity price, the data expression is as follows:

[0073] in, Y PPA ( a ) indicates the main power supply node of new energy source. a Power supply to new energy main power consumption nodes b The revenue from electricity generation, K PPA ( b ) represents the main electricity consumption node of new energy sources. b Power is drawn from the main power supply node of new energy sources. a Electricity costs, λ PPA This refers to the long-term contract electricity price signed between renewable energy power suppliers and renewable energy power consumers. For the power distribution matrix, the first Line number Column elements represent the set of power supply nodes. S sup The Middle Power supply nodes For the set of power consumption nodes S use The Middle individual electricity consumption nodes The amount of electrical energy allocated.

[0074] Figure 2 In the example, the on-grid electricity price of the main renewable energy power generation unit λ up The off-grid electricity price of the power grid λdown Long-term contract electricity prices signed between renewable energy power supply entities and renewable energy consumption entities λ PPA Table 3 shows the power generation revenue and electricity costs of the three new energy entities under electricity price systems of 0.4 yuan / kWh, 0.6 yuan / kWh, and 0.35 yuan / kWh, respectively.

[0075] Table 3. Revenue from power generation and electricity costs of new energy sources (unit: yuan)

[0076] The data in the second row and first column of Table 3 is 16.4, representing that the revenue of renewable energy entity 1 from the power grid is 16.4 yuan. The other elements in the table are similar. The revenue from power supply to renewable energy entity 2 is 12.26 yuan, the revenue from power generation and grid connection of renewable energy entity 1 is 16.4 yuan, and the revenue from power supply to renewable energy entity 2 is 4.73 yuan. The revenue from power generation and grid connection of renewable energy entity 3 is 0.3 yuan, and the revenue from power supply to renewable energy entity 2 is 3.14 yuan. Similarly, the electricity cost for renewable energy entity 2 to draw power from the power grid is 12.26 yuan, the electricity cost to draw power from renewable energy entity 1 is 4.73 yuan, and the electricity cost to draw power from renewable energy entity 3 is 3.14 yuan. It can be seen that the formula in step S5 separates the revenue from the renewable energy supply node and the cost of the electricity consumption node, ensuring the fairness of the settlement between the renewable energy power plant and the power grid. At the same time, the above formula performs electricity settlement based on the electricity allocation results, ensuring the rationality of the settlement for different renewable energy entities.

[0077] Furthermore, this embodiment also proposes a power system for calculating complex power flow across multiple nodes in a power system, including a processor and a computer-readable storage medium. The computer-readable storage medium stores a computer program, which is executed by the processor to implement the steps of the power system power flow calculation method for complex power flow across multiple nodes in this embodiment.

[0078] In summary, this invention proposes a method and system for calculating complex power flow in a multi-node grid-connected power system. Based on the meter installation location and main wiring diagram, the meters are categorized into renewable energy-side meters and grid-side meters, and forward and reverse energy and power data are read. For a specific renewable energy source's collector line, the power flow direction is determined, the reverse energy increment is accumulated to calculate the total self-generated energy, and the forward energy increment is accumulated to calculate the total load, thus obtaining the net power generation and achieving virtual combined metering for a single renewable energy source. Then, the system nodes are divided into grid-side nodes, renewable energy source nodes, and bus nodes. Bus node injected energy vectors, tie-line energy vectors, and correlation matrices are constructed. Through correlation matrix decomposition and matrix operations, a distribution factor matrix is ​​obtained, and the energy contribution ratio of the power supply node to the power consumption node is calculated using the distribution factor matrix. Based on the contribution ratio of the power supply node to the power consumption node, the injected energy, the supplied energy, and the consumed energy, an energy distribution matrix is ​​calculated, clarifying the energy distribution from the power supply node to the power consumption node and solving the problems of unclear power flow paths and unclear energy distribution in multi-entity grid connection scenarios. During settlement, differentiated settlement methods are used, employing both on-grid and contract prices, for different transaction scenarios involving settlement between renewable energy sources and the grid, and between renewable energy entities, to ensure fair transactions. Compared to existing technologies, this method offers the following advantages: 1. Addressing the issue of self-generated and self-consumed electricity from renewable energy power plants being metered as grid-connected energy at the power source end, while the electricity consumed at the load end is repeatedly metered as off-grid energy, resulting in two different settlement processes with varying pricing mechanisms, this invention calculates the total self-generated energy and total load energy based on the power direction of the power meter on the individual renewable energy power plant's collection line. This yields the net power generation, enabling virtual combined metering for a single renewable energy power plant. This prioritizes self-consumption and surplus electricity fed into the grid by renewable energy power plants in energy calculations, facilitating accurate energy settlement between renewable energy power plants and the grid, as well as between renewable energy power plants and other renewable energy power plants.

[0079] 2. In scenarios involving multiple renewable energy generators connected to the grid, the electricity generated by these generators first flows into a common bus and is then distributed. A portion may supply its own load, another portion may be directly consumed by other users within the grid, and the remainder may be transmitted to the upper-level grid. This results in multiple distribution paths for the electricity, making it difficult to distinguish the actual power flow path. This invention divides the system nodes into grid-side nodes, renewable energy generator nodes, and bus nodes. It constructs an injected power vector for the bus nodes, a tie-line power vector, and an correlation matrix. Through correlation matrix decomposition and matrix operations, a distribution factor matrix is ​​obtained. This matrix is ​​then used to calculate the power contribution ratio of the power supply node to the power consumption node. Based on the contribution ratio of the power supply node to the power consumption node, the injected power, the supplied power, and the consumed power, a power distribution matrix is ​​calculated, clarifying the specific flow direction of electricity from the power supply node to the power consumption node and solving the problem of unclear power flow paths in multi-generator grid connection scenarios.

[0080] 3. The current settlement system only uses inter-station electricity meters to measure transactions between renewable energy power plants and grid companies, lacking a clear settlement rule system for the mutual power supply and consumption between different power plant entities. For the settlement of electricity costs between renewable energy power supply nodes and grid-side power consumption nodes, and between renewable energy power supply nodes and renewable energy power consumption nodes, this invention uses allocated electricity to accurately settle the revenue and costs of each power supply and consumption node, employing differentiated settlement using both grid-connected and contracted electricity prices to ensure fair transactions.

[0081] Those skilled in the art will understand that embodiments of this application can be provided as methods, systems, or computer program products. Therefore, this application can take the form of a completely hardware embodiment, a completely software embodiment, or an embodiment combining software and hardware aspects. Furthermore, this application can take the form of a computer program product embodied on one or more computer-readable storage media (including, but not limited to, disk storage, CD-ROM, optical storage, etc.) containing computer-usable program code. This application is described with reference to flowchart illustrations and / or block diagrams of methods, apparatus (systems), and computer program products according to embodiments of this application. It will be understood that each block of the flowchart illustrations and / or block diagrams, and combinations of blocks in the flowchart illustrations and / or block diagrams, can be implemented by computer program instructions. These computer program instructions can be provided to a processor of a general-purpose computer, special-purpose computer, embedded processor, or other programmable data processing apparatus to produce a machine, such that the instructions, which execute via the processor of the computer or other programmable data processing apparatus, create a machine for implementing the process. Figure 1 One or more processes and / or boxes Figure 1 The computer program instructions may also be stored in a computer-readable storage medium that can direct a computer or other programmable data processing device to operate in a particular manner, such that the instructions stored in the computer-readable storage medium produce an article of manufacture including instruction means, which are implemented in a process Figure 1 One or more processes and / or boxes Figure 1 The functions specified in one or more boxes. These computer program instructions may also be loaded onto a computer or other programmable data processing apparatus to cause a series of operational steps to be performed on the computer or other programmable apparatus to produce a computer-implemented process, thereby providing instructions that execute on the computer or other programmable apparatus for implementing the process. Figure 1 One or more processes and / or boxes Figure 1 The steps of the function specified in one or more boxes.

[0082] The above description is merely a preferred embodiment of the present invention. The scope of protection of the present invention is not limited to the above embodiments. All technical solutions falling within the scope of the present invention's concept are within the scope of protection of the present invention. It should be noted that for those skilled in the art, any improvements and modifications made without departing from the principles of the present invention should also be considered within the scope of protection of the present invention.

Claims

1. A method for power system multi-node grid-connected complex power flow energy settlement, characterized in that, The method comprises the following steps: obtaining forward power, reverse power and power data of all electric energy meters in a specified time interval; judging the power flow direction of the corresponding node according to the power data of the electric energy meter corresponding to each new energy power generation subject, calculating the self-generation capacity, total load capacity and net generation capacity of each new energy power generation subject according to the power flow direction of the corresponding node and the forward power and reverse power of the electric energy meter corresponding to each new energy power generation subject, and realizing the virtual metering of each new energy subject; constructing a distribution factor matrix between the nodes of the grid-side electric energy meter and the nodes of the virtual metering of different new energy subjects to represent the power distribution between bus nodes, one end of the grid-side electric energy meter being a grid-side node and the other end being a bus node, one end of the virtual metering of the new energy subject being a new energy subject node and the other end being a bus node, the grid-side node including a grid-side power supply node and a grid-side power consumption node, and the new energy subject node including a new energy subject power supply node and a new energy subject power consumption node; calculating the power contribution ratio of the power supply node to the power consumption node by using the distribution factor matrix, and calculating the power distribution matrix according to the contribution ratio of the power supply node to the power consumption node, the injected power, the power supply power and the power consumption power, the power supply node including the grid-side power supply node and the new energy subject power supply node, and the power consumption node including the grid-side power consumption node and the new energy subject power consumption node; carrying out power settlement between the new energy subject power supply node and the grid-side power consumption node according to the power distribution matrix and the on-grid price of the new energy power generation subject, carrying out power settlement between the grid-side power supply node and the new energy side power consumption node according to the power distribution matrix and the off-grid price of the grid, and carrying out power settlement between the new energy subject power supply node and the new energy subject power consumption node according to the power distribution matrix and the contract price.

2. The method of claim 1, wherein, When judging the power flow direction of the corresponding node according to the power data of the electric energy meter corresponding to each new energy power generation subject, the following steps are included: Electric energy meter n q,c Power data Less than 0, the electric energy meter n q,c The power flow direction of the corresponding node is from the corresponding node to the bus node; Electric energy meter n q,c Power data Greater than 0, the electric energy meter n q,c The power flow direction of the corresponding node is from the bus node to the corresponding node.

3. The method of claim 2, wherein, When calculating the self-generation capacity, total load and net generation capacity of each new energy power generation subject according to the power flow direction of the corresponding node and the forward power and reverse power of the electric energy meter corresponding to each new energy power generation subject, the following steps are included: Electric energy meter n q,c If the power flow direction of the corresponding node is from the corresponding node to the bus node, the state label of the corresponding node is set as the first value, and the self-generated electric energy of the corresponding node is equal to the reverse electric energy increment of the electric energy meter n q,c ​ Electric energy meter n q,c The power flow direction of the corresponding node is from the bus node to the corresponding node, the state label of the corresponding node is set to the second value, and the load electric energy of the corresponding node is equal to the positive forward electric energy increment of the electric energy meter n q,c ​ accumulating the self-generation capacity of all nodes corresponding to the current new energy power generation subject and having a state label of a first value to obtain the self-generation capacity of the current new energy power generation subject; accumulating the load capacity of all nodes corresponding to the current new energy power generation subject and having a state label of a second value to obtain the total load capacity of the current new energy power generation subject; subtracting the total load capacity from the self-generation capacity of the current new energy power generation subject to obtain the net generation capacity of the current new energy power generation subject.

4. The method of claim 3, wherein, Before constructing the distribution factor matrix between the nodes of the grid-side electric energy meter and the nodes of the virtual metering of different new energy subjects, the following step of dividing the power supply node and the power consumption node is included, which specifically includes: if the power data of the grid-side electric energy meter corresponding to the current grid-side node is less than 0, the current grid-side node is a grid-side power supply node, and if the power data of the grid-side electric energy meter corresponding to the current grid-side node is greater than 0, the current grid-side node is a grid-side power consumption node; If the net power generation of the new energy power generation subject corresponding to the current new energy subject node is greater than 0, the current new energy subject node is a new energy subject power supply node; if the net power generation of the new energy power generation subject corresponding to the current new energy subject node is less than 0, the current new energy subject node is a new energy subject power consumption node.

5. The method of claim 1, wherein, The mathematical expression of the distribution factor matrix is as follows: wherein, I is a unit matrix of size M × M , A is a distribution factor matrix of size M × M , whose elements A ( x , y ) represent the proportion of the bus node x to the bus node y , x , y =1,2,…, M , and represent the vectors F and E Z is a diagonal matrix constructed by taking the reciprocal of the elements, is the transpose matrix of the matrix , and the matrix is a tie-line incidence matrix B in which all non-zero elements are 1, and the matrix is a tie-line incidence matrix B in which all non-zero elements are -1, and the element B B d,m represents the connection relationship between the d th tie line and the m th node, B d,m =1 represents that the starting node of the d th tie line is m ; B d,m =-1 represents that the ending node of the d th tie line is m ; B d,m =0 represents that the d th tie line is not connected to the node m ;​ Vector F The vector of the liaison line electric energy is represented, the mathematical expression is as follows: F =[ F 1, F 2,…, F d ,…, F D ] T wherein d = 1, 2,... D , the d transmission of electrical energy F d The formula is as follows: wherein, denotes the positive energy increment of the energy meter of the d th tie line, and denote the positive energy of the energy meter of the d th tie line at the end and start time of the specified time interval, respectively, denotes the negative energy increment of the energy meter of the d th tie line, and denote the negative energy of the energy meter of the d th tie line at the end and start time of the specified time interval, respectively, is the power data of the energy meter of the d th tie line at the end time of the specified time interval; Vector E Z The vector of the power injection at a bus node is represented by the mathematical expression: E Z =[ ] T Wherein, the bus node mid m The mathematical expression of the injected power is as follows: In the formula, S m The upstream tie line set of the bus node mid m , The transmission electric energy of the i d m th tie line, The supply electric energy of the supply node sup m connected with the bus node mid v . If the supply node is a grid-side supply node, the value of the supply electric energy is the reverse electric energy increment of the grid-side electric energy meter corresponding to the grid-side supply node. If the supply node is a new energy main body supply node, the value of the supply electric energy is the absolute value of the net power generation of the new energy main body corresponding to the new energy main body supply node.

6. The method of claim 5, wherein, The mathematical expression of the power distribution matrix is as follows: wherein, G V-diag represents a vector G V constructed diagonal matrix, vector G V represents a power supply node power supply energy vector, G V the first R g elements of the vector are power supply energies of the grid side power supply nodes, and the last Q g elements are power supply energies of the new energy main body power supply nodes; A sub denotes the submatrix of distribution factors between the bus nodes downstream of the supply nodes and the bus nodes upstream of the load nodes, A sub ( v , u ) denotes the element of the matrix in row v , column u , which is given by the mathematical expression A sub ( v , u )= A ( G ind ( v ), L ind ( u )) In the formula, A G ind v L ind u represents the element in the distribution factor matrix G ind v row L ind u column, G ind v represents the mth bus node in the bus node vector G ind downstream of the power supply node, v L ind u represents the mth bus node in the bus node vector L ind upstream of the power consumption node, u G ind The first n elements in the vector are the bus nodes corresponding to the power supply nodes on the grid side, and the last m elements are the bus nodes corresponding to the power supply nodes of the new energy main body. R g Q g The first n elements in the vector are the bus nodes corresponding to the power consumption nodes on the grid side, and the last m elements are the bus nodes corresponding to the power consumption nodes of the new energy main body. L ind R l Q g ​​​​​​​​​​​​​​ , represents a vector Copy in the vertical direction V Form a new matrix, vector is the power consumption of the power node Power energy vector L U transpose, L U The first R l element in the vector is the power consumption of the power node Power energy, and the last Q g element is the power consumption of the new energy main body power node , represents the multiplication of E Z the new vector obtained by indexing the vector E Z ( G ind ) the matrix obtained by replicating horizontally V the new matrix obtained by indexing the matrix V, V is the total number of supply nodes, represents the multiplication of E Z-sub the matrix obtained by taking the reciprocal element-wise, represents the Hadamard product of 7. The method of claim 6, wherein, When the power settlement between the new energy subject power supply node and the grid side power consumption node is performed according to the power distribution matrix and the on-grid price of the new energy power generation subject, the mathematical expression is as follows: wherein, Y up a represents a new energy main power supply node new a power supply to the grid side power consumption node, λ up is the on-grid electricity price of the new energy power generation main body, is the element of the power distribution matrix in the row and the column, which represents the power supply node set S sup in the power supply node of the S use power consumption node in the power consumption node set S sup The first R g elements in the set are grid side power supply nodes, and the last Q g elements are new energy main power supply nodes, S use The first R l elements in the set are grid side power consumption nodes, and the last Q g elements are new energy main power consumption nodes.​ 8. The method of claim 7, wherein, When the power settlement between the grid side power supply node and the new energy side power consumption node is performed according to the power distribution matrix and the off-grid price of the grid, the data expression is as follows: wherein, K down (b) represents a new energy main body electricity node new b electricity consumption cost from the power grid side power supply node, λ down the lower grid electricity price of the power grid, is the element of the power distribution matrix row and column, representing the power supply node set S sup the th power supply node in the power supply node set S use the th electricity consumption node in the electricity consumption node set.

9. The method of claim 8, wherein, When the power settlement between the new energy subject power supply node and the new energy subject power consumption node is performed according to the power distribution matrix and the contract price, the data expression is as follows: in, Y PPA ( a ) indicates the main power supply node of new energy source. a Power supply to new energy main power consumption nodes b The revenue from electricity generation, K PPA ( b ) represents the main electricity consumption node of new energy sources. b Power is drawn from the main power supply node of new energy sources. a Electricity costs, λ PPA This refers to the long-term contract electricity price signed between renewable energy power suppliers and renewable energy power consumers. For the power distribution matrix, the first Line 1 Column elements represent the set of power supply nodes. S sup The Middle Power supply nodes For the set of power consumption nodes S use The Middle individual electricity consumption nodes The amount of electrical energy distributed.

10. A multi-node grid-connected complex power flow electricity settlement system for a power system, characterized in that, The computer readable storage medium stores a computer program, and the computer program is executed by the processor to realize the steps of the power system multi-node grid-connected complex power flow power settlement method in any one of claims 1-9.

Images