Line impedance online real-time accurate inversion calculation method and apparatus, device and medium

Abstract

A line impedance online real-time accurate inversion calculation method, comprising: selecting a plurality of measurement nodes on a main line and branch lines on the basis of sections to be measured; constructing a main line voltage matrix on the basis of the measurement nodes on the main line; constructing a branch line voltage matrix on the basis of the measurement nodes on each branch line; merging the main line voltage matrix and the branch line voltage matrix into a first voltage matrix (S5); constructing a second voltage matrix on the basis of a voltage difference between a line outgoing node and a measurement node closest thereto in a power distribution network (S6); merging the first voltage matrix and the second voltage matrix into a total voltage matrix (S7); merging the total voltage matrix, a total current column vector and a total admittance column vector into a total voltage matrix equation (S10); using the total voltage matrix equation to determine the admittance of all the sections to be measured (S11); and then determining the corresponding impedance (S12). According to the solution, more accurate line impedance can be calculated and inverted on line in real time, providing accurate line parameters and topological information for real-time setting and real-time load flow calculation of a power grid protection fixed value, and supporting applications such as power grid planning, scheduling, operation control and generation-grid-load-storage interaction.

Description

Line impedance online real-time accurate inversion calculation method, device, equipment and medium Technical Field

[0001] The present invention relates to the field of power electronics technology, and in particular to a method, device, equipment and medium for online real-time accurate inversion calculation of line impedance. Background Art

[0002] To serve the national energy transformation and "dual carbon" goals, a large number of distributed new energy and diversified loads are connected to the power grid, requiring the power grid to achieve multi-directional coordination and flexible interaction. Traditional power grids have shown problems of insufficient transparency and digitalization, especially the distribution network structure is complex and the equipment is numerous. There is an urgent need for transparency of power grid topology and electrical and physical parameters to truly realize the digital twin power grid.

[0003] Current technical means usually use power outages to measure line impedance, which affects power supply reliability and places a heavy workload on the distribution network. There is no effective method to perceive the entire power grid in real time and accurately obtain line impedance online.

[0004] Summary of the Invention

[0005] Based on this, it is necessary to propose an online real-time accurate inversion calculation method, device, equipment and medium for line impedance to address the above problems.

[0006] To achieve the above-mentioned objectives, the present application provides, in a first aspect, a method for online real-time accurate inversion calculation of line impedance, the method comprising:

[0007] Selecting several measurement nodes in the distribution network according to the section to be measured;

[0008] Determine a number of measurement nodes located on the main line and a number of measurement nodes located on each branch line;

[0009] Forming a main line voltage matrix according to the voltage difference between every two adjacent measurement nodes among the measurement nodes of the main line;

[0010] A branch line voltage matrix is ​​formed according to the voltage difference between the branch head node and the closest measurement node in the measurement nodes of each branch line, and the voltage difference between every two adjacent measurement nodes;

[0011] Combining the main line voltage matrix and the branch line voltage matrix to form a first voltage matrix;

[0012] A current column vector is formed according to the current flowing out of the measurement node of each main line in the distribution network and the current flowing out of the measurement node of each branch line along the direction of the distribution network line;

[0013] Construct an admittance column vector according to the admittance of each section to be measured;

[0014] Combining the first voltage matrix, the current column vector, and the admittance column vector to form a first voltage matrix equation for inversely calculating the first admittance;

[0015] Constructing a second voltage matrix equation based on the voltage difference between the node at the main line outlet and the closest measurement node in the distribution network, the outflow current of the node at the main line outlet, and the admittance between the node at the main line outlet and the closest measurement node;

[0016] Merging the second voltage matrix equation into the first voltage matrix equation to form a total voltage matrix equation for inversely calculating the admittance of the section to be measured;

[0017] Determining the admittance of all segments to be measured according to the total voltage matrix equation;

[0018] The corresponding impedance is determined according to the admittance of all the segments to be measured.

[0019] In some embodiments, selecting a plurality of measurement nodes in the distribution network according to the segment to be measured specifically includes:

[0020] Determine all monitoring nodes included in the segment to be tested;

[0021] Sorting all the monitoring nodes according to the direction of the distribution network line;

[0022] The head-end nodes or the end-end nodes among all the monitoring nodes after sorting are removed, and the remaining monitoring nodes are used as measurement nodes.

[0023] In some embodiments, the main line voltage matrix The elements of are:

[0024] Where i = 1 to n, n is the number of main line nodes, and j is 1 to n.

[0025] In some embodiments, the branch voltage matrix The elements of are:

[0026] Among them, if the i node has a branch, then k=i, k is the total number of branches, i=1 to p, p is the number of branch line nodes, and j is 1 to p.

[0027] In some implementations, merging the main line voltage matrix and the branch line voltage matrix to form a first voltage matrix specifically includes:

[0028] The main line voltage matrix is ​​combined with the diagonal line of the branch line voltage matrix of each branch line to form the first voltage matrix.

[0029] In some embodiments, the elements of the current column vector are:

[0030] In some embodiments, each element of the admittance column vector is: i1 =Y(i-1.i).

[0031] In some embodiments, the first voltage matrix equation is: in, is the main line voltage matrix, Y 主 is the principal line admittance vector, is the main line current vector, is the branch voltage matrix, Y 支k is the branch admittance vector, is the branch current vector.

[0032] In some embodiments, the second voltage matrix equation is merged into the first voltage matrix equation to form a total voltage matrix equation for inverse calculation of the admittance of the section to be measured, including:

[0033] Expand the voltage difference between the node at the main line outlet of the distribution network and the measurement node closest to it into a second voltage matrix in the second voltage matrix equation of 1*j, U11 is U0-U1, when j>1, U1j=0, and add it to the first row of the first voltage matrix in the first voltage matrix equation to form the total voltage matrix in the total voltage matrix equation;

[0034] The outflow current at the node of the main line outlet in the second voltage matrix equation is Added as the first row of the current column vector in the first voltage matrix equation, constituting the total current column vector in the total voltage matrix equation;

[0035] The total voltage matrix, the admittance column vector, and the total current column vector constitute the total voltage matrix equation used for inverse calculation of the admittance of the section to be measured.

[0036] In some embodiments, the total voltage matrix equation is: in, is the main line voltage matrix, Y 主 is the principal line admittance vector, is the main line current vector, is the branch voltage matrix, Y 支k is the branch admittance vector, is the branch current vector, The pressure difference between the node where the main line exits and the closest measurement node, It is the outflow current of the node where the main line exits.

[0037] To achieve the above-mentioned purpose, the second aspect of the present application provides an online real-time accurate inversion calculation device for line impedance, the device comprising:

[0038] A selection module is used to select several measurement nodes in the distribution network according to the section to be measured;

[0039] The selection module is further configured to determine a number of measurement nodes located on the main line and a number of measurement nodes located on each branch line;

[0040] A voltage module, configured to form a main line voltage matrix according to the voltage difference between every two adjacent measurement nodes of the main line;

[0041] The voltage module is further configured to form a branch voltage matrix based on the voltage difference between the branch head node and the closest measurement node in the measurement nodes of each branch line, and the voltage difference between every two adjacent measurement nodes;

[0042] The voltage module is further configured to combine the main line voltage matrix and the branch line voltage matrix to form a first voltage matrix;

[0043] A current module, configured to construct a current column vector based on the current flowing out of a measurement node of each main line in the distribution network and the current flowing out of a measurement node of each branch line along the direction of the distribution network line;

[0044] Admittance module, used to construct an admittance column vector according to the admittance of each section to be measured;

[0045] a calculation module, configured to combine the first voltage matrix, the current column vector, and the admittance column vector to form a first voltage matrix equation for inversely calculating the first admittance;

[0046] The calculation module is further configured to construct a second voltage matrix equation based on the voltage difference between the node at the main line outlet and the closest measurement node in the distribution network, the outflow current of the node at the main line outlet, and the admittance between the node at the main line outlet and the closest measurement node;

[0047] The calculation module is further configured to merge the second voltage matrix equation into the first voltage matrix equation to form a total voltage matrix equation for inverse calculation of the admittance of the section to be measured;

[0048] The calculation module is further configured to determine the admittance of all segments to be measured based on the total voltage matrix equation;

[0049] The calculation module is further configured to determine corresponding impedances according to the admittances of all the segments to be measured.

[0050] To achieve the above-mentioned objective, the third aspect of the present application provides a computer device, comprising a memory and a processor, wherein the memory stores a computer program, and when the computer program is executed by the processor, the processor performs the following steps:

[0051] Selecting several measurement nodes in the distribution network according to the section to be measured;

[0052] Determine a number of measurement nodes located on the main line and a number of measurement nodes located on each branch line;

[0053] Forming a main line voltage matrix according to the voltage difference between every two adjacent measurement nodes among the measurement nodes of the main line;

[0054] A branch line voltage matrix is ​​formed according to the voltage difference between the branch head node and the closest measurement node in the measurement nodes of each branch line, and the voltage difference between every two adjacent measurement nodes;

[0055] Combining the main line voltage matrix and the branch line voltage matrix to form a first voltage matrix;

[0056] A current column vector is formed according to the current flowing out of the measurement node of each main line in the distribution network and the current flowing out of the measurement node of each branch line along the direction of the distribution network line;

[0057] Construct an admittance column vector according to the admittance of each section to be measured;

[0058] Combining the first voltage matrix, the current column vector, and the admittance column vector to form a first voltage matrix equation for inversely calculating the first admittance;

[0059] Constructing a second voltage matrix equation based on the voltage difference between the node at the main line outlet and the closest measurement node in the distribution network, the outflow current of the node at the main line outlet, and the admittance between the node at the main line outlet and the closest measurement node;

[0060] Merging the second voltage matrix equation into the first voltage matrix equation to form a total voltage matrix equation for inversely calculating the admittance of the section to be measured;

[0061] Determining the admittance of all segments to be measured according to the total voltage matrix equation;

[0062] The corresponding impedance is determined according to the admittance of all the segments to be measured.

[0063] To achieve the above-mentioned object, the fourth aspect of the present application provides a computer-readable storage medium storing a computer program. When the computer program is executed by a processor, the processor performs the following steps:

[0064] Selecting several measurement nodes in the distribution network according to the section to be measured;

[0065] Determine a number of measurement nodes located on the main line and a number of measurement nodes located on each branch line;

[0066] Forming a main line voltage matrix according to the voltage difference between every two adjacent measurement nodes among the measurement nodes of the main line;

[0067] A branch line voltage matrix is ​​formed according to the voltage difference between the branch head node and the closest measurement node in the measurement nodes of each branch line, and the voltage difference between every two adjacent measurement nodes;

[0068] Combining the main line voltage matrix and the branch line voltage matrix to form a first voltage matrix;

[0069] A current column vector is formed according to the current flowing out of the measurement node of each main line in the distribution network and the current flowing out of the measurement node of each branch line along the direction of the distribution network line;

[0070] Construct an admittance column vector according to the admittance of each section to be measured;

[0071] Combining the first voltage matrix, the current column vector, and the admittance column vector to form a first voltage matrix equation for inversely calculating the first admittance;

[0072] Constructing a second voltage matrix equation based on the voltage difference between the node at the main line outlet and the closest measurement node in the distribution network, the outflow current of the node at the main line outlet, and the admittance between the node at the main line outlet and the closest measurement node;

[0073] Merging the second voltage matrix equation into the first voltage matrix equation to form a total voltage matrix equation for inversely calculating the admittance of the section to be measured;

[0074] Determining the admittance of all segments to be measured according to the total voltage matrix equation;

[0075] The corresponding impedance is determined according to the admittance of all the segments to be measured.

[0076] The embodiments of the present invention have the following beneficial effects:

[0077] In order to solve the problem of measuring the impedance of the entire power grid line in real time online without power outages, the present invention proposes a method for accurately measuring and reconstructing the line parameters of the distribution network online by using the voltage and current measured on the low-voltage side of the transformer. Several measurement nodes in the distribution network are selected according to the section to be measured; the total voltage matrix, total current column vector, and admittance column vector are determined according to the measurement nodes, and are combined to form a total voltage matrix equation for inverse calculation of admittance; the distribution network is reconstructed according to the total voltage matrix equation, and accurate line impedance is inverted online in real time, providing accurate line parameters and topology information for real-time setting of power grid protection values ​​and real-time power flow calculation, supporting applications such as power grid planning, scheduling, operation control, and source-grid-load-storage interaction, and facilitating the construction of new power systems. In addition, a second voltage matrix composed of the voltage difference between the outgoing node and the measurement node closest to it in the distribution network is added, which has stronger convergence and more accurate line impedance calculated by online real-time inversion. BRIEF DESCRIPTION OF THE DRAWINGS

[0078] In order to more clearly illustrate the embodiments of the present invention or the technical solutions in the prior art, the following briefly introduces the drawings required for use in the embodiments or the description of the prior art. Obviously, the drawings described below are only some embodiments of the present invention. For ordinary technicians in this field, other drawings can be obtained based on these drawings without paying any creative work.

[0079] in:

[0080] FIG1 is a flow chart of a method for online real-time accurate inversion calculation of line impedance according to one embodiment;

[0081] FIG2 is an example diagram of the basic principle of network reconstruction in one embodiment;

[0082] FIG3 is a structural diagram of a distribution network in one embodiment;

[0083] FIG4 is a schematic diagram of the main line of a distribution network structure diagram in one embodiment;

[0084] FIG5 is a schematic diagram of a vector analysis simulation model of a distribution transformer in one embodiment;

[0085] FIG6 is a positive sequence circuit diagram of a distribution transformer in one embodiment;

[0086] FIG7 is a negative sequence circuit diagram of a distribution transformer in one embodiment;

[0087] FIG8 is a structural diagram of a device for online real-time accurate inversion calculation of line impedance in one embodiment;

[0088] FIG9 is a schematic diagram of the structure of a computer device in one embodiment;

[0089] FIG10 is a schematic diagram of the structure of a computer-readable storage medium in one embodiment. DETAILED DESCRIPTION

[0090] The following will clearly and completely describe the technical solutions in the embodiments of the present invention in conjunction with the accompanying drawings. Obviously, the described embodiments are only part of the embodiments of the present invention, not all of the embodiments. All other embodiments obtained by ordinary technicians in this field based on the embodiments of the present invention without making any creative efforts shall fall within the scope of protection of the present invention.

[0091] In an embodiment of the present application, a method for online real-time accurate inversion calculation of line impedance is provided. Please refer to Figure 1, which is a flow chart of an online real-time accurate inversion calculation method for line impedance in one embodiment. The principle of the present application is to reconstruct the distribution network using matrix equations and inversely calculate the impedance of the line. Specifically: based on the vector voltage and current vector between two nodes, the line impedance between the nodes is calculated.

[0092] For example, according to the network reconstruction basic principle example diagram shown in Figure 2, write the equation Y1 and Y2 can be solved, and the corresponding line impedance between nodes can be calculated based on Y1 and Y2.

[0093] In an embodiment of the present application, the method for online real-time accurate inversion calculation of line impedance includes steps S1 to S12 shown in FIG1 .

[0094] Step S1 selects several measurement nodes in the distribution network according to the section to be measured;

[0095] In some embodiments, selecting a plurality of measurement nodes in the distribution network according to the segment to be measured specifically includes:

[0096] Determine all monitoring nodes included in the segment to be tested;

[0097] Sorting all the monitoring nodes according to the direction of the distribution network line;

[0098] The head-end nodes or the end-end nodes among all the monitoring nodes after sorting are removed, and the remaining monitoring nodes are used as measurement nodes.

[0099] Step S2, determining a number of measurement nodes located on the main line and a number of measurement nodes located on each branch line;

[0100] Step S3, forming a main line voltage matrix according to the voltage difference between every two adjacent measurement nodes among the measurement nodes of the main line;

[0101] In some embodiments, the main line voltage matrix The elements of are:

[0102] Where i = 1 to n, n is the number of main line nodes, and j is 1 to n.

[0103] Step S4, constructing a branch line voltage matrix according to the voltage difference between the branch head node and the closest measurement node among the measurement nodes of each branch line, and the voltage difference between every two adjacent measurement nodes;

[0104] In some embodiments, the branch voltage matrix The elements of are:

[0105] Among them, if the i node has a branch, then k=i, k is the total number of branches, i=1 to p, p is the number of branch line nodes, and j is 1 to p.

[0106] Step S5, merging the main line voltage matrix and the branch line voltage matrix to form a first voltage matrix;

[0107] In some implementations, merging the main line voltage matrix and the branch line voltage matrix to form a first voltage matrix specifically includes:

[0108] The main line voltage matrix is ​​combined with the diagonal line of the branch line voltage matrix of each branch line to form the first voltage matrix.

[0109] Specifically, the first element of the first row of the branch voltage matrix of the first branch is placed in the last column of the row below the last element in the main line voltage matrix;

[0110] Place the first element of the first row of the N-th branch voltage matrix in the remaining branches in the next row to the last element in the N-1-th branch voltage matrix and then in the next column;

[0111] Finally, the first voltage matrix is ​​constructed.

[0112] Step S6, forming a current column vector based on the current flowing out of the measurement node of each main line in the distribution network and the current flowing out of the measurement node of each branch line along the direction of the distribution network line;

[0113] In some embodiments, the elements of the current column vector are:

[0114] Specifically, the current column vector includes and Main line current column vector I 主 The elements of are: Where i = 1 to n, n is the number of main line nodes; the branch line current column vector I 支k The elements of are: Where k is the branch line number, i = 1 to p, and p is the number of branch line nodes.

[0115] Step S7, constructing an admittance column vector according to the admittance of each section to be measured;

[0116] In some embodiments, each element of the admittance column vector is: i1 =Y(i-1.i).

[0117] Specifically, the admittance column vector includes Y 主 and Y 支k , admittance column vector Y 主 The elements of Y are: i1 =Y(i-1.i), where i=1 to n, n is the number of main line nodes; the admittance column vector Y 支k The elements of are: Where k is the branch line number, i = 1 to p, and p is the number of branch line nodes.

[0118] Step S8, combining the first voltage matrix, current column vector, and admittance column vector to form a first voltage matrix equation for inverse calculation of admittance;

[0119] In some embodiments, the first voltage matrix equation is: in, is the main line voltage matrix, Y 主 is the principal line admittance vector, is the main line current vector, is the branch voltage matrix, Y 支k is the branch admittance vector, is the branch current vector.

[0120] Specifically, the main line voltage matrix Branch voltage matrix Combined to form the first voltage matrix is the admittance column vector, is the current column vector.

[0121] Step S9, constructing a second voltage matrix equation based on the voltage difference between the node at the main line outlet and the nearest measurement node in the distribution network, the outflow current of the node at the main line outlet, and the admittance between the node at the main line outlet and the nearest measurement node;

[0122] Step S10, merging the second voltage matrix equation into the first voltage matrix equation to form a total voltage matrix equation for inverse calculation of the admittance of the section to be measured;

[0123] In some embodiments, the second voltage matrix equation is merged into the first voltage matrix equation to form a total voltage matrix equation for inverse calculation of the admittance of the section to be measured, including:

[0124] Expand the voltage difference between the node at the main line outlet of the distribution network and the measurement node closest to it into a second voltage matrix in the second voltage matrix equation of 1*j, U11 is U0-U1, when j>1, U1j=0, and add it to the first row of the first voltage matrix in the first voltage matrix equation to form the total voltage matrix in the total voltage matrix equation;

[0125] The outflow current at the node of the main line outlet in the second voltage matrix equation is Added as the first row of the current column vector in the first voltage matrix equation, constituting the total current column vector in the total voltage matrix equation;

[0126] The total voltage matrix, the admittance column vector, and the total current column vector constitute the total voltage matrix equation used for inverse calculation of the admittance of the section to be measured.

[0127] It should be noted that the admittance column vector is composed of the admittance of each section to be measured, including the admittance between the node at the main line outlet and the measurement node closest to it in the second voltage matrix equation. Therefore, there is no need to add a row to the admittance column vector. The admittance column vector is the total admittance column vector.

[0128] In some embodiments, the total voltage matrix equation is: in, is the main line voltage matrix, Y 主 is the principal line admittance vector, is the main line current vector, is the branch voltage matrix, Y 支k is the branch admittance vector, is the branch current vector, The pressure difference between the node where the main line exits and the closest measurement node, It is the outflow current of the node where the main line exits.

[0129] It should be noted that: the first voltage matrix equation is a homogeneous linear equation group, and the total voltage matrix equation composed of merging the second voltage matrix equation into the first voltage matrix equation is a non-homogeneous linear equation group. With one more constraint equation, the calculated admittance and final impedance will be more accurate. If only the first voltage matrix equation is used, if the voltage measurement is wrong, it will cause inaccurate parameter calculation, which is difficult to detect. With one more constraint equation, it is possible to determine whether the measurement data is wrong by calculating whether the matrix converges. In this way, problems in the online real-time inversion calculation of the line impedance can be discovered in time, and the inversion results are more precise and the accuracy is more guaranteed.

[0130] Specifically, the main line voltage matrix Branch voltage matrix Combined to form the first voltage matrix First voltage matrix and the second voltage matrix Combined to form the total voltage matrix, is the admittance column vector, and is the total current column vector.

[0131] In some embodiments, Add one row as the first row to the first voltage matrix U and the current column vector I respectively, and the admittance column vector Y remains unchanged. The added row of the first voltage matrix U is: At this time, the first voltage matrix U becomes a (n+1)*n matrix;

[0132] The current column vector I adds one row to the first row, and the added row At this time, the current column vector I becomes a (n+1)*1 matrix.

[0133] in is the voltage at the outlet, is the current at the outlet, For the voltage of the main line node 1, which is the node closest to the outgoing line, the following equation is obtained:

[0134] In some embodiments, the following steps are constructed

[0135] Assume that the main line has m nodes, there are k branches, each branch has n nodes, then there are j=m+kn nodes in total.

[0136] The main line voltage matrix is ​​constructed as an m*m matrix, the main line current matrix is ​​constructed as an m*1 matrix, and the main line admittance matrix is ​​constructed as an m*1 matrix;

[0137] Construct each branch voltage matrix as an n*n matrix, each branch current matrix I as an n*1 matrix, and each branch admittance matrix as an n*1 matrix.

[0138] according to Combine all main line matrices and branch line matrices to form U 底 ×Y all =I 底 , where U 底 is the j*j matrix, which is the _th voltage matrix, Y all is a j*1 matrix, that is, the admittance column vector; I 底 is a j*1 matrix, which is the current column vector;

[0139] Construct the second voltage matrix U m0 is a 1*j matrix, the first element is The remaining elements are 0.

[0140] Will U m0 As the first line with U 底 Combined to form the complete total voltage matrix U all , then the total voltage matrix =(j+1)*j; the current at the outlet As the first line with I 底 Combined to form a complete current column vector I all , at this time the current column vector I all is a (j+1)*1 matrix;

[0141] Write the matrix equation as U all Y all =I all get

[0142] Specifically, for example, the voltage matrix equation is written to reconstruct the distribution network structure shown in Figure 3. At this time, the voltage matrix equation is:

[0143] Step S11, determining the admittance of all segments to be measured according to the voltage matrix equation;

[0144] Specifically, according to the voltage matrix equation After measuring the voltage and current of each node, we can calculate That is, the admittance of all segments to be measured.

[0145] Step S12: determining the corresponding impedance according to the admittance of all the segments to be measured.

[0146] Specifically, take The reciprocal of each element in the equation gives the corresponding impedance between each node.

[0147] By adopting the technical solution of this embodiment, a method is proposed to realize accurate online measurement and reconstruction of distribution network line parameters by using the voltage and current measured on the low-voltage side of the transformer, and the first voltage matrix and the second voltage matrix are merged into a total voltage matrix; the total voltage matrix, the total current column vector, and the admittance column vector are merged into a total voltage matrix equation; which is used to determine the admittance of all sections to be measured, and then determine the corresponding impedance. Such a voltage matrix equation has strong convergence due to the addition of the second matrix equation, and the line impedance inverted by online real-time calculation is more accurate, which can provide accurate line parameters and topology information for real-time setting of grid protection constants and real-time power flow calculation, support grid planning, scheduling, operation control, and source-grid-load-storage interaction, and help build new power systems.

[0148] In a specific embodiment, when reconstructing the distribution network using matrix equations, it is necessary to convert the voltage and current vectors on the low-voltage side of the distribution transformer to the high-voltage side line, and then obtain the voltage vectors of each measurement node in the section to be measured, for example And get the current vector of each measurement node of the section to be measured, for example It is used to subsequently substitute into the voltage matrix equation to calculate the admittance. Usually, the current vector of each measurement node of the measured section is obtained through current vector analysis, and the voltage vector of each measurement node of the measured section is obtained through voltage vector analysis.

[0149] Furthermore, the connection group of the power transformers in the distribution network is usually Dyn11. The following is based on the analysis of the Dyn11 distribution transformer and the conversion of the low-voltage side voltage and current vectors to the high-voltage side vectors. That is, in a specific embodiment, the current vector analysis and voltage vector analysis of each measurement node in the measured section can be referred to the following examples.

[0150] In a specific embodiment, the current vector is analyzed as follows:

[0151] First, the 0.4kV low-voltage side current is calculated using the formula The symmetrical component method decomposes the current into positive sequence, negative sequence and zero sequence current.

[0152] in, are the positive sequence, negative sequence and zero sequence components of the phase A current on the low voltage side of the distribution transformer, For the three-phase current of A, B, and C on the low-voltage side of the distribution transformer,

[0153] Since the high voltage side is a delta connection, the line current on the high voltage side does not include the zero sequence current. Ignoring the influence of the excitation impedance shunt of the distribution transformer (if the no-load current and no-load loss are known, they can be calculated), the current on the low voltage side of the distribution transformer is converted to the positive sequence current on the high voltage side. and negative sequence current They are

[0154] Further, the positive sequence current on the high voltage side of the distribution transformer and negative sequence Composite high-voltage side A, B, and C phase line currents for:

[0155] In summary, the high voltage side line current vector is completed. Conversion.

[0156] Furthermore, a simulation calculation is performed using the distribution transformer vector analysis simulation model schematic diagram shown in FIG5 to verify the high-voltage side line current vector Is the conversion correct?

[0157] The distribution transformer voltage ratio is 10kV / 0.4kV, the connection group is Dyn1, and an unbalanced load is set on the 0.4kV side of the distribution transformer. The distribution transformer parameters here mainly set the rated voltage, rated capacity, and leakage reactance on both sides, and ignore copper loss and iron loss. The distribution transformer parameters of the distribution transformer vector analysis simulation model shown in Figure 5 mainly set the rated voltage, rated capacity, and leakage reactance on both sides, and ignore copper loss and iron loss.

[0158] Assume that the rated voltage ratio of the distribution transformer is 25, the distribution transformer capacity is 50kVA and the leakage reactance is 0.1pu. The leakage reactance converted to the low voltage side is: X T2 =0.4×0.4÷50×1000×0.1=0.32jΩ. The current vector calculated by the theoretical formula and the PSCAD simulation is as follows:

[0159] (1) Calculation of positive sequence current on the high-voltage side

[0160] 1) Measure the positive sequence current on the low voltage side 369.8∠-130.2°A

[0161] 2) Calculate the positive sequence current on the high voltage side

[0162] 3) The positive sequence current on the high-voltage side calculated by PSCAD simulation is 14.8∠-160.25°A

[0163] (2) Calculation of negative sequence current on the high voltage side

[0164] 1) Measure the negative sequence current on the low voltage side 135.6∠106.9°A

[0165] 2) Calculate the negative sequence current on the high voltage side

[0166] 3) The negative sequence voltage on the high voltage side calculated by PSCAD simulation is 5.4∠136.9°

[0167] (3) Calculate the current of each phase on the high voltage side

[0168] (4) Phase currents on the high-voltage side obtained by PSCAD simulation

[0169] The maximum amplitude error between the calculated high-voltage side current and the high-voltage side current calculated by PSCAD simulation is 0.02A, and the maximum phase error is 0.06°, as shown in Table 1. This shows that the conclusion of the current vector analysis is correct. The current vector error obtained by this calculation is very small and has high accuracy. It can be used to substitute into the voltage matrix equation of this application to calculate the admittance and then calculate the corresponding impedance.

[0170] Table 1:

[0171] In a specific embodiment, the voltage vector is analyzed as follows:

[0172] The voltage on the 0.4kV low-voltage side of the distribution transformer is calculated using the formula The symmetrical component method decomposes the voltage into positive sequence, negative sequence and zero sequence voltage.

[0173] The positive-sequence power supply voltage and negative-sequence power supply voltage on the 0.4kV low-voltage side of the distribution transformer are further calculated based on the positive-sequence circuit and negative-sequence circuit on the 0.4kV low-voltage side of the distribution transformer. The positive-sequence circuit of the Dyn11 distribution transformer is shown in the positive-sequence circuit diagram of the distribution transformer in Figure 6.

[0174] is the positive sequence power supply voltage on the 0.4kV low voltage side of the distribution transformer, is the positive sequence voltage calculated by the symmetrical component method of the voltage measured at the 0.4kV low-voltage side terminal, is the positive sequence current calculated by the symmetrical component method of the current measured at the 0.4kV low-voltage side terminal, X T2+ It is the positive sequence impedance converted from the high voltage side of the distribution transformer to the low voltage side (only the leakage reactance is considered, and the excitation reactance and winding resistance are ignored).

[0175] According to Figure 6, the positive sequence power supply voltage on the 0.4kV low-voltage side of the distribution transformer is obtained

[0176] further The positive sequence voltage converted to the 10kV high voltage side of the distribution transformer is

[0177] in, is the positive sequence voltage on the high voltage side of the distribution transformer, and k is the distribution transformer ratio.

[0178] Furthermore, the negative sequence loop impedance of the Dyn11 distribution transformer is equal to the positive sequence loop impedance. The negative sequence loop of the distribution transformer is shown in the negative sequence loop diagram of the distribution transformer in Figure 7.

[0179] The negative sequence power supply voltage on the 0.4kV low voltage side of the distribution transformer, The negative sequence voltage is calculated by the symmetrical component method of the voltage measured at the 0.4kV low-voltage side terminal. is the negative sequence current calculated by the symmetrical component method of the current measured at the 0.4kV low-voltage side terminal, X T2- It is the negative sequence impedance from the high voltage side of the distribution transformer to the low voltage side.

[0180] According to Figure 7, the negative sequence power supply voltage on the 0.4kV low-voltage side of the distribution transformer is

[0181] further The positive sequence voltage converted to the high voltage side is

[0182] in, is the positive sequence voltage on the high voltage side of the distribution transformer, and k is the distribution transformer ratio.

[0183] Furthermore, the zero-sequence impedance on the high-voltage side of the distribution transformer Dyn11 is very small, and the zero-sequence voltage is approximately 0, so the positive-sequence voltage on the high-voltage side of the distribution transformer is and negative sequence voltage Synthesized high-voltage side A, B, and C phase voltages for:

[0184] In summary, the high-voltage side phase voltage is completed. Conversion.

[0185] Furthermore, a simulation calculation is performed using the distribution transformer vector analysis simulation model schematic diagram shown in FIG5 to verify the high-voltage side phase voltage vector Is the conversion correct?

[0186] The distribution transformer voltage ratio is 10kV / 0.4kV, the connection group is Dyn1, and an unbalanced load is set on the 0.4kV side of the distribution transformer. The distribution transformer parameters here mainly set the rated voltage, rated capacity, and leakage reactance on both sides, and ignore copper loss and iron loss. The distribution transformer parameters of the distribution transformer vector analysis simulation model shown in Figure 5 mainly set the rated voltage, rated capacity, and leakage reactance on both sides, and ignore copper loss and iron loss.

[0187] (1) Calculation of positive sequence voltage on the high-voltage side

[0188] 1) PSCAD simulation calculation of low voltage side positive sequence voltage 126.12∠-78.6°V

[0189] 2) PSCAD simulation calculation of low voltage side positive sequence current 369.8∠-130.2°A

[0190] 3) Calculate the positive sequence power supply voltage on the low voltage side

[0191] 4) Calculate the positive sequence voltage on the high voltage side

[0192] 5) The positive sequence voltage on the high-voltage side calculated by PSCAD simulation is 5773.12∠-90.04°

[0193] (2) Calculation of negative sequence voltage on the high voltage side

[0194] 1) PSCAD simulation calculation of low voltage side negative sequence voltage 43.4∠-16.9°V

[0195] 2) PSCAD simulation calculation of low voltage side negative sequence current 135.6∠106.9°A

[0196] 3) Calculate the open circuit negative sequence voltage on the low voltage side

[0197] 4) Calculate the negative sequence voltage on the high voltage side

[0198] 5) The negative sequence voltage on the high voltage side calculated by PSCAD simulation is 0.05∠-54.5°

[0199] (3) Calculate the voltage of each phase on the high voltage side

[0200] (4) Phase voltages on the high-voltage side obtained by PSCAD simulation

[0201] It can be seen that the maximum error between the calculated high-voltage side voltage and the high-voltage side voltage amplitude calculated by PSCAD simulation is 3V, the error is 0.05%, and the maximum phase error is 0.01°, as shown in Table 2. It can be seen that the conclusion of the voltage vector analysis is correct. The voltage vector error obtained by this calculation is very small and has high accuracy. It can be used to substitute into the voltage matrix equation of this application to calculate the admittance and then calculate the corresponding impedance.

[0202] Table 2

[0203] In an embodiment of the present application, an online real-time accurate inversion calculation device for line impedance is provided. Please refer to Figure 8, which is a structural diagram of the online real-time accurate inversion calculation device for line impedance in one embodiment. The online real-time accurate inversion calculation device for line impedance includes a selection module 201, a voltage module 202, a current module 203, an admittance module 204, and a calculation module 205.

[0204] The selection module 201 is configured to select a number of measurement nodes in the distribution network according to the section to be measured;

[0205] The selection module 201 is further configured to determine a number of measurement nodes located on the main line and measurement nodes located on each branch line;

[0206] The voltage module 202 is configured to form a main line voltage matrix according to the voltage difference between every two adjacent measurement nodes of the main line;

[0207] The voltage module 202 is further configured to form a branch voltage matrix based on the voltage difference between the branch head node and the closest measurement node in the measurement nodes of each branch, and the voltage difference between every two adjacent measurement nodes;

[0208] The voltage module 202 is further configured to combine the main line voltage matrix and the branch line voltage matrix to form a first voltage matrix;

[0209] The current module 203 is configured to form a current column vector according to the current flowing out of the measurement node of each main line in the distribution network and the current flowing out of the measurement node of each branch line along the direction of the distribution network line;

[0210] The admittance module 204 is configured to construct an admittance column vector according to the admittance of each segment to be measured;

[0211] A calculation module 205 is configured to combine the first voltage matrix, the current column vector, and the admittance column vector to form a first voltage matrix equation for inverse calculation of the first admittance;

[0212] The calculation module 205 is further configured to construct a second voltage matrix equation based on the voltage difference between the node at the main line outlet and the nearest measurement node in the distribution network, the outflow current of the node at the main line outlet, and the admittance between the node at the main line outlet and the nearest measurement node;

[0213] The calculation module 205 is further configured to merge the second voltage matrix equation into the first voltage matrix equation to form a total voltage matrix equation for inverse calculation of the admittance of the section to be measured;

[0214] The calculation module 205 is further configured to determine the admittance of all segments to be measured according to the total voltage matrix equation;

[0215] The calculation module 205 is further configured to determine corresponding impedances according to the admittances of all the segments to be measured.

[0216] In some embodiments, the selection module 201 is further configured to determine all monitoring nodes included in the segment to be measured;

[0217] Sorting all the monitoring nodes according to the direction of the distribution network line;

[0218] The head-end nodes or the end-end nodes among all the monitoring nodes after sorting are removed, and the remaining monitoring nodes are used as measurement nodes.

[0219] In some embodiments, the voltage module 202 is further configured to determine the main line voltage matrix The elements of are: Where i = 1 to n, n is the number of main line nodes, and j is 1 to n.

[0220] In some embodiments, the voltage module 202 is further configured to determine the branch voltage matrix The elements of Among them, if the i node has a branch, then k=i, k is the total number of branches, i=1 to p, p is the number of branch line nodes, and j is 1 to p.

[0221] In some embodiments, the voltage module 202 is further configured to combine the main line voltage matrix with the diagonal line of the branch line voltage matrix of each branch line to form the first voltage matrix.

[0222] In some embodiments, the current module 203 is further configured to determine each element of the current column vector as:

[0223] In some embodiments, the admittance module 204 is further configured to determine each element of the admittance column vector as: i1 =Y(i-1.i).

[0224] In some implementations, the calculation module 205 is further configured to determine the first voltage matrix equation as: in, is the main line voltage matrix, Y 主is the principal line admittance vector, is the main line current vector, is the branch voltage matrix, Y 支k is the branch admittance vector, is the branch current vector.

[0225] In some embodiments, the calculation module 205 is further configured to merge the second voltage matrix equation into the first voltage matrix equation to form a total voltage matrix equation for inverse calculation of the admittance of the section to be measured, including:

[0226] Expand the voltage difference between the node at the main line outlet of the distribution network and the measurement node closest to it into a second voltage matrix in the second voltage matrix equation of 1*j, U11 is U0-U1, when j>1, U1j=0, and add it to the first row of the first voltage matrix in the first voltage matrix equation to form the total voltage matrix in the total voltage matrix equation;

[0227] The outflow current at the node of the main line outlet in the second voltage matrix equation is Added as the first row of the current column vector in the first voltage matrix equation, constituting the total current column vector in the total voltage matrix equation;

[0228] The total voltage matrix, the admittance column vector, and the total current column vector constitute the total voltage matrix equation used for inverse calculation of the admittance of the section to be measured.

[0229] In some implementations, the calculation module 205 is further configured to determine the total voltage matrix equation as: in, is the main line voltage matrix, Y 主 is the principal line admittance vector, is the main line current vector, is the branch voltage matrix, Y 支k is the branch admittance vector, is the branch current vector, The pressure difference between the node where the main line exits and the closest measurement node, It is the outflow current of the node where the main line exits.

[0230] For other details about how the modules in the device for online real-time accurate inversion calculation of line impedance implement the above technical solution, please refer to the description of the device for online real-time accurate inversion calculation of line impedance provided above, which will not be repeated here.

[0231] In an embodiment of the present application, a computer device is provided. Please refer to FIG9 , which is a schematic diagram of the structure of a computer device in one embodiment. The device includes a memory 301 and a processor 302. The memory 301 stores a computer program. When the computer program is executed by the processor 302, the processor 302 performs the following steps:

[0232] Selecting several measurement nodes in the distribution network according to the section to be measured;

[0233] Determine a number of measurement nodes located on the main line and a number of measurement nodes located on each branch line;

[0234] Forming a main line voltage matrix according to the voltage difference between every two adjacent measurement nodes among the measurement nodes of the main line;

[0235] A branch line voltage matrix is ​​formed according to the voltage difference between the branch head node and the closest measurement node in the measurement nodes of each branch line, and the voltage difference between each two adjacent measurement nodes;

[0236] Combining the main line voltage matrix and the branch line voltage matrix to form a first voltage matrix;

[0237] A current column vector is formed according to the current flowing out of the measurement node of each main line in the distribution network and the current flowing out of the measurement node of each branch line along the direction of the distribution network line;

[0238] Construct an admittance column vector according to the admittance of each section to be measured;

[0239] Combining the first voltage matrix, the current column vector, and the admittance column vector to form a first voltage matrix equation for inversely calculating the first admittance;

[0240] Constructing a second voltage matrix equation based on the voltage difference between the node at the main line outlet and the closest measurement node in the distribution network, the outflow current of the node at the main line outlet, and the admittance between the node at the main line outlet and the closest measurement node;

[0241] Merging the second voltage matrix equation into the first voltage matrix equation to form a total voltage matrix equation for inversely calculating the admittance of the section to be measured;

[0242] Determining the admittance of all segments to be measured according to the total voltage matrix equation;

[0243] The corresponding impedance is determined according to the admittance of all the segments to be measured.

[0244] Among them, the processor 302 can also be called a CPU (Central Processing Unit), and the processor 302 may be an integrated circuit chip with signal processing capabilities; the processor 302 can also be a general-purpose processor, DSP (Digital Signal Processor), ASIC (Application Specific Integrated Circuit), FPGA (Field Programmable Gate Array) or other programmable logic devices, discrete gate or transistor logic devices, discrete hardware components, among which the general-purpose processor can be a microprocessor or the processor 302 can also be any conventional processor, etc.

[0245] In an embodiment of the present application, a computer-readable storage medium is provided. Please refer to FIG. 10 , which is a schematic diagram of the structure of a computer-readable storage medium in one embodiment. The storage medium stores a readable computer program 401. The computer program 401 may be stored in the storage medium in the form of a software product and includes several instructions for causing a computer device (which may be a personal computer, a server, or a network device, etc.) or a processor to perform the following steps:

[0246] Selecting several measurement nodes in the distribution network according to the section to be measured;

[0247] Determine a number of measurement nodes located on the main line and a number of measurement nodes located on each branch line;

[0248] Forming a main line voltage matrix according to the voltage difference between every two adjacent measurement nodes among the measurement nodes of the main line;

[0249] A branch line voltage matrix is ​​formed according to the voltage difference between the branch head node and the closest measurement node in the measurement nodes of each branch line, and the voltage difference between each two adjacent measurement nodes;

[0250] Combining the main line voltage matrix and the branch line voltage matrix to form a first voltage matrix;

[0251] A current column vector is formed according to the current flowing out of the measurement node of each main line in the distribution network and the current flowing out of the measurement node of each branch line along the direction of the distribution network line;

[0252] Construct an admittance column vector according to the admittance of each section to be measured;

[0253] Combining the first voltage matrix, the current column vector, and the admittance column vector to form a first voltage matrix equation for inversely calculating the first admittance;

[0254] Constructing a second voltage matrix equation based on the voltage difference between the node at the main line outlet and the closest measurement node in the distribution network, the outflow current of the node at the main line outlet, and the admittance between the node at the main line outlet and the closest measurement node;

[0255] Merging the second voltage matrix equation into the first voltage matrix equation to form a total voltage matrix equation for inversely calculating the admittance of the section to be measured;

[0256] Determining the admittance of all segments to be measured according to the total voltage matrix equation;

[0257] The corresponding impedance is determined according to the admittance of all the segments to be measured.

[0258] The aforementioned storage media include: USB flash drives, mobile hard drives, magnetic disks or optical disks, ROM (Read-Only Memory), RAM (Random Access Memory), and other media that can store program codes, or terminal devices such as computers, service machines, mobile phones, and tablets.

[0259] Those skilled in the art will appreciate that all or part of the processes in the above-mentioned embodiments can be implemented by instructing the relevant hardware through a computer program. The program can be stored in a non-volatile computer-readable storage medium. When the program is executed, it can include the processes of the embodiments of the above-mentioned methods. Among them, any reference to memory, storage, database or other media used in the embodiments provided in this application can include non-volatile and / or volatile memory. Non-volatile memory can include read-only memory (ROM), programmable ROM (PROM), electrically programmable ROM (EPROM), electrically erasable programmable ROM (EEPROM) or flash memory. Volatile memory can include random access memory (RAM) or external cache memory. By way of illustration and not limitation, RAM is available in various forms, such as static RAM (SRAM), dynamic RAM (DRAM), synchronous DRAM (SDRAM), double data rate SDRAM (DDRSDRAM), enhanced SDRAM (ESDRAM), synchronous link (Synchlink) DRAM (SLDRAM), memory bus (Rambus) direct RAM (RDRAM), direct memory bus dynamic RAM (DRDRAM), and memory bus dynamic RAM (RDRAM).

[0260] The technical features of the above embodiments can be combined arbitrarily. To make the description concise, not all possible combinations of the technical features in the above embodiments are described. However, as long as there is no contradiction in the combination of these technical features, they should be considered to be within the scope of this specification.

[0261] The above-described embodiments merely represent several implementation methods of the present application. While the descriptions are relatively specific and detailed, they should not be construed as limiting the scope of the present invention. It should be noted that a person of ordinary skill in the art may make various modifications and improvements without departing from the spirit of the present invention, all of which fall within the scope of protection of the present invention. Therefore, the scope of protection of the present invention shall be determined by the appended claims.

Claims

1. A method for online real-time accurate inversion calculation of line impedance, characterized in that: The method comprises: Selecting several measurement nodes in the distribution network according to the section to be measured; Determine a number of measurement nodes located on the main line and a number of measurement nodes located on each branch line; Forming a main line voltage matrix according to the voltage difference between every two adjacent measurement nodes among the measurement nodes of the main line; A branch line voltage matrix is formed according to the voltage difference between the branch head node and the closest measurement node in the measurement nodes of each branch line, and the voltage difference between every two adjacent measurement nodes; Combining the main line voltage matrix and the branch line voltage matrix to form a first voltage matrix; A current column vector is formed according to the current flowing out of the measurement node of each main line in the distribution network and the current flowing out of the measurement node of each branch line along the direction of the distribution network line; Construct an admittance column vector according to the admittance of each section to be measured; Combining the first voltage matrix, the current column vector, and the admittance column vector to form a first voltage matrix equation for inversely calculating the first admittance; A second voltage matrix equation is formed according to the voltage difference between the node at the main line outlet and the closest measurement node in the distribution network, the outflow current of the node at the main line outlet, and the admittance between the node at the main line outlet and the closest measurement node; Merging the second voltage matrix equation into the first voltage matrix equation to form a total voltage matrix equation for inversely calculating the admittance of the section to be measured; Determining the admittance of all segments to be measured according to the total voltage matrix equation; The corresponding impedance is determined according to the admittance of all the segments to be measured.

2. The method for online real-time accurate inversion calculation of line impedance according to claim 1, characterized in that: The step of selecting a plurality of measurement nodes in the distribution network according to the section to be measured specifically includes: Determine all monitoring nodes included in the segment to be tested; Sort all monitoring nodes according to the direction of the distribution network line; The head-end nodes or the end-end nodes among all the monitoring nodes after sorting are removed, and the remaining monitoring nodes are used as measurement nodes.

3. The method for online real-time accurate inversion calculation of line impedance according to claim 1 or 2, characterized in that: The main line voltage matrix The elements of are: Where i = 1 to n, n is the number of main line nodes, and j is 1 to n.

4. The method for online real-time accurate inversion calculation of line impedance according to claim 3, characterized in that: The branch voltage matrix The elements of are: Among them, if the i node has a branch, then k=i, k is the total number of branches, i=1 to p, p is the number of branch line nodes, and j is 1 to p.

5. The method for online real-time accurate inversion calculation of line impedance according to claim 4, characterized in that: The step of merging the main line voltage matrix and the branch line voltage matrix to form a first voltage matrix specifically includes: The main line voltage matrix is combined with the diagonal line of the branch line voltage matrix of each branch line to form the first voltage matrix.

6. The method for online real-time accurate inversion calculation of line impedance according to claim 5, characterized in that: The elements of the current column vector are:

7. The method for online real-time accurate inversion calculation of line impedance according to claim 6, characterized in that: The elements of the admittance column vector are: i1 =Y(i-1.i).

8. The method for online real-time accurate inversion calculation of line impedance according to claim 7, characterized in that: The first voltage matrix equation is: in, is the main line voltage matrix, Y 主 is the principal line admittance vector, is the main line current vector, is the branch voltage matrix, Y 支k is the branch admittance vector, is the branch current vector.

9. The method for online real-time accurate inversion calculation of line impedance according to claim 8, characterized in that: The second voltage matrix equation is merged into the first voltage matrix equation to form a total voltage matrix equation for inverse calculation of the admittance of the section to be measured, including: Expand the voltage difference between the node at the main line outlet of the distribution network and the measurement node closest to it into a second voltage matrix in the second voltage matrix equation of 1*j, U11 is U0-U1, when j>1, U1j=0, and add it to the first row of the first voltage matrix in the first voltage matrix equation to form the total voltage matrix in the total voltage matrix equation; The outflow current at the node of the main line outlet in the second voltage matrix equation is Added as the first row of the current column vector in the first voltage matrix equation, constituting the total current column vector in the total voltage matrix equation; The total voltage matrix, the admittance column vector, and the total current column vector constitute the total voltage matrix equation used for inverse calculation of the admittance of the section to be measured.

10. The method for online real-time accurate inversion calculation of line impedance according to claim 9, characterized in that: The total voltage matrix equation is: in, is the main line voltage matrix, Y 主 is the principal line admittance vector, is the main line current vector, is the branch voltage matrix, Y 支k is the branch admittance vector, is the branch current vector, The pressure difference between the node where the main line exits and the closest measurement node, It is the outflow current of the node where the main line exits.

11. A line impedance online real-time accurate inversion calculation device, characterized in that: The device comprises: A selection module is used to select several measurement nodes in the distribution network according to the section to be measured; The selection module is further configured to determine a number of measurement nodes located on the main line and a number of measurement nodes located on each branch line; A voltage module, configured to form a main line voltage matrix according to the voltage difference between every two adjacent measurement nodes of the main line; The voltage module is further configured to form a branch voltage matrix based on the voltage difference between the branch head node and the closest measurement node in the measurement nodes of each branch line, and the voltage difference between every two adjacent measurement nodes; The voltage module is further configured to combine the main line voltage matrix and the branch line voltage matrix to form a first voltage matrix; A current module, configured to construct a current column vector based on the current flowing out of a measurement node of each main line in the distribution network and the current flowing out of a measurement node of each branch line along the direction of the distribution network line; Admittance module, used to construct an admittance column vector according to the admittance of each section to be measured; a calculation module, configured to combine the first voltage matrix, the current column vector, and the admittance column vector to form a first voltage matrix equation for inversely calculating the first admittance; The calculation module is further configured to construct a second voltage matrix equation based on the voltage difference between the node at the main line outlet and the closest measurement node in the distribution network, the outflow current of the node at the main line outlet, and the admittance between the node at the main line outlet and the closest measurement node; The calculation module is further configured to merge the second voltage matrix equation into the first voltage matrix equation to form a total voltage matrix equation for inversely calculating the admittance of the section to be measured; The calculation module is further used to determine the admittance of all segments to be measured according to the total voltage matrix equation; The calculation module is further configured to determine corresponding impedances according to the admittances of all the segments to be measured.

12. A computer device comprising a memory and a processor, wherein the memory stores a computer program, and when the computer program is executed by the processor, the processor is caused to perform the steps of the method according to any one of claims 1 to 10.

13. A computer-readable storage medium storing a computer program, wherein when the computer program is executed by a processor, the processor is caused to perform the steps of the method according to any one of claims 1 to 10.

Images