CN117217008B - Digital twin model construction method of power transmission network transformer - Google Patents
Digital twin model construction method of power transmission network transformer Download PDFInfo
- Publication number
- CN117217008B CN117217008B CN202311197441.0A CN202311197441A CN117217008B CN 117217008 B CN117217008 B CN 117217008B CN 202311197441 A CN202311197441 A CN 202311197441A CN 117217008 B CN117217008 B CN 117217008B
- Authority
- CN
- China
- Prior art keywords
- transformer
- phase
- voltage
- loop
- shifting
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
- 230000005540 biological transmission Effects 0.000 title claims abstract description 49
- 238000010276 construction Methods 0.000 title claims abstract description 11
- 239000011159 matrix material Substances 0.000 claims abstract description 63
- 239000013598 vector Substances 0.000 claims abstract description 37
- 238000000034 method Methods 0.000 claims abstract description 34
- 238000012795 verification Methods 0.000 claims abstract description 13
- 238000012546 transfer Methods 0.000 claims abstract description 7
- 238000004804 winding Methods 0.000 claims description 56
- 230000010363 phase shift Effects 0.000 claims description 53
- 230000008030 elimination Effects 0.000 claims description 7
- 238000003379 elimination reaction Methods 0.000 claims description 7
- 230000002441 reversible effect Effects 0.000 claims description 6
- 230000008569 process Effects 0.000 claims description 4
- 230000009466 transformation Effects 0.000 claims description 4
- 230000007704 transition Effects 0.000 claims description 4
- 238000006467 substitution reaction Methods 0.000 claims description 3
- 230000001052 transient effect Effects 0.000 description 11
- 238000004364 calculation method Methods 0.000 description 4
- 238000010586 diagram Methods 0.000 description 4
- 238000012937 correction Methods 0.000 description 3
- 238000004088 simulation Methods 0.000 description 3
- 101100063069 Caenorhabditis elegans deg-1 gene Proteins 0.000 description 2
- 230000001360 synchronised effect Effects 0.000 description 2
- 238000004422 calculation algorithm Methods 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 230000003993 interaction Effects 0.000 description 1
- 230000009191 jumping Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000007935 neutral effect Effects 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
Landscapes
- Supply And Distribution Of Alternating Current (AREA)
Abstract
The invention discloses a digital twin model construction method of a power transmission network transformer, which comprises the following steps: obtaining short-circuit impedance parameters of each transformer in a power transmission network, and constructing a digital twin model of the power transmission network transformer based on the short-circuit impedance parameters of each transformer; based on a digital twin model of the power transmission network transformer, constructing an adjacent matrix and branch transfer item vectors, and determining a phase-shifting loop verification result in the digital twin model by using the adjacent matrix and the branch transfer item vectors, wherein the phase-shifting loop verification result comprises no phase-shifting loop or a phase-shifting loop; if the phase-shifting ring verification result is that the phase-shifting ring exists, positioning the phase-shifting ring to obtain a phase-shifting ring positioning result; and correcting the transformer connection method on the branch in the phase-shifting ring according to the positioning result of the phase-shifting ring so as to eliminate the phase-shifting ring and obtain the final digital twin module. The invention searches the phase shifting ring in the transmission network digital twin model and eliminates the phase shifting ring, thereby ensuring the correctness of the transmission network digital twin model.
Description
Technical Field
The invention belongs to the field of digital twin modeling of a power transmission network, and particularly relates to a digital twin model construction method of a power transmission network transformer.
Background
Along with the evolution of a novel power system, a digital twin model of a power transmission network is built, and the online synchronous evolution of the digital twin model and the power transmission network state is ensured. The construction of the digital twin model of the power transmission network transformer often adopts the existing electromechanical transient data as a basis, however, the direct conversion of the current manually maintained online electromechanical transient data into the electromagnetic transient model can be problematic, so that the digital twin model based on the electromagnetic transient refined simulation is difficult to normally operate: 1) The three-winding transformer can generate negative resistance after the resistance is converted, and the negative resistance can cause the electromagnetic transient simulation to generate a phenomenon of divergence of results; 2) A large number of electromagnetic looped networks exist in an alternating current transmission network, under the condition of error of a main transformer connection method, a phase shift loop possibly exists in the network, and a large reactive loop current is generated in digital twin simulation and deviates from a steady-state operating point, so that the analysis result of a digital twin model of the network is incorrect. The problems of transformer parameters and wiring often need to be searched in the whole system, if the manual searching is performed, time and labor are wasted, and automatic synchronous updating evolution of the digital twin model of the power transmission network cannot be achieved.
Therefore, when constructing a power transmission network transformer model, it is necessary to develop an automatic searching and correcting method for the two problems, so as to ensure the correctness of the digital twin model of the power transmission network.
Disclosure of Invention
The invention aims to provide a digital twin model construction method of a power transmission network transformer, which solves the problems existing in the prior art.
The invention is realized by the following technical scheme:
a digital twin model construction method of a power transmission network transformer comprises the following steps:
obtaining short-circuit impedance parameters of each transformer in a power transmission network, and constructing a digital twin model of the power transmission network transformer based on the short-circuit impedance parameters of each transformer;
Constructing an adjacency matrix and branch transition term vectors based on a digital twin model of a power transmission network transformer, and determining a phase-shifting loop verification result in the digital twin model according to the adjacency matrix and the branch transition term vectors, wherein the phase-shifting loop verification result comprises no phase-shifting loop or a phase-shifting loop;
If the phase-shifting ring verification result is that the phase-shifting ring exists, positioning the phase-shifting ring to obtain a phase-shifting ring positioning result;
And correcting a transformer connection method on a branch in the phase-shifting ring according to the positioning result of the phase-shifting ring so as to eliminate the phase-shifting ring and obtain a final digital twin module.
In one possible embodiment, obtaining a short circuit impedance parameter for each transformer in a power transmission network includes:
A1, acquiring the type of a transformer and nameplate parameters, judging whether the type of the transformer is a two-winding transformer, if so, entering a step A2, otherwise, judging that the transformer is a three-winding transformer, and entering a step A3;
A2, determining short-circuit impedance parameters of the two-winding transformer according to nameplate parameters of the transformer, wherein the short-circuit impedance parameters are as follows:
Wherein U 0% represents the transformer short-circuit voltage percentage, P 0 represents the transformer short-circuit loss, U B represents the transformer high-voltage side rated voltage, S B represents the transformer high-voltage side rated capacity, X represents the transformer winding leakage reactance, and R represents the transformer winding resistance;
a3, determining short-circuit impedance parameters of the three-winding transformer according to nameplate parameters of the transformer, wherein the short-circuit impedance parameters are as follows:
Wherein, X 1 represents a transformer high-voltage winding leakage reactance, X 2 represents a transformer medium-voltage winding leakage reactance, X 3 represents a transformer low-voltage winding leakage reactance, U 120% represents a transformer high-medium voltage short-circuit voltage percentage, U 230% represents a transformer medium-voltage low-voltage short-circuit voltage percentage, U 130% represents a transformer high-voltage low-voltage short-circuit voltage percentage, U 1B represents a transformer high-voltage side rated voltage, U 2B represents a transformer medium-voltage side rated voltage, U 3B represents a transformer low-voltage side rated voltage, S 1B represents a transformer high-voltage side rated capacity, S 2B represents a transformer medium-voltage side rated capacity, S 3B represents a transformer low-voltage side rated capacity, P 120 represents a transformer high-medium voltage short-circuit loss, P 230 represents a transformer medium-voltage low voltage short-circuit loss, P 130 represents a transformer high-voltage winding resistance, R 1 represents a transformer medium-voltage winding resistance, R 2 represents a transformer medium-voltage winding resistance, and R 3 represents a transformer low-voltage winding resistance;
a4, obtaining short-circuit impedance parameters of each transformer according to the methods in the step A2 and the step A3.
In one possible implementation, a digital twinning model of a grid transformer is constructed based on short circuit impedance parameters of the respective transformers, comprising: and constructing a digital twin model of the power transmission network transformer in an equivalent substitution mode according to the short-circuit impedance parameters of each transformer.
In one possible implementation, the adjacency matrix a is:
Wherein A (i) represents the ith row of the adjacent matrix, i is more than or equal to 1 and less than or equal to b, and b is the number of power network branches; corresponding branch i, k is the starting node of branch i, m represents the ending node of branch i, i.e. the ith row A (i: of matrix A) of which the elements except k columns and m columns are all 0; n is the number of electrical nodes of the power network.
In one possible implementation, the branch term vector Δθ is:
Δθ＝[Δθ1 Δθ2 … Δθm]T
Wherein, delta theta i represents the phase difference of the ith branch caused by the connection of the transformer, T represents the transpose.
In one possible implementation manner, positioning the phase shift ring to obtain a positioning result of the phase shift ring includes:
adopting Gaussian elimination method to process the adjacent matrix A, converting the adjacent matrix A into an upper triangular matrix U through series elementary line transformation, and obtaining the following steps:
KA＝U
Wherein K represents a reversible matrix;
According to the branch term transfer vector, the loop phase angle difference vector delta theta * is obtained as follows:
Δθ*＝KΔθ
Wherein Δθ represents a branch term vector;
determining rows which are all 0 in the upper triangular matrix U to obtain an end p rows; when no row which is 0 exists in the upper triangular matrix U, determining that no loop exists, namely the positioning result of the phase-shifting loop is no phase-shifting loop, otherwise, determining the phase-shifting loop;
judging whether the last p rows in the loop phase angle difference vector delta theta * are all 0, if yes, determining that no loop exists, namely the positioning result of the phase shifting loop is that no phase shifting loop exists, otherwise, determining that non-zero elements exist;
and determining that the number of non-zero elements in the last p rows of the loop phase angle difference vector delta theta * is q, and determining that q phase shift loops exist as a phase shift loop check result.
In one possible implementation manner, positioning the phase shift ring to obtain a positioning result of the phase shift ring includes:
Determining the row where q non-zero elements are located in the last p rows of the loop phase angle difference vector delta theta * to obtain a target row number;
determining the columns of non-zero elements in the row corresponding to the target row number in the reversible matrix K to obtain the target column number;
and taking the branch number corresponding to the target column number as a branch for forming the phase shifting ring to obtain a phase shifting ring positioning result.
In one possible implementation manner, according to the positioning result of the phase shift ring, correcting the transformer connection on the branch in the phase shift ring to eliminate the phase shift ring, including:
determining whether a transformer connection corresponding to a branch circuit forming the phase-shifting ring in the phase-shifting ring positioning result meets standard requirements, if so, determining that the transformer connection is not required to be corrected, otherwise, determining that the transformer connection is required to be corrected;
And obtaining all target transformers which need to be corrected in the branches forming the phase-shifting ring, correcting the transformer connection method of the target transformers, and eliminating the phase-shifting ring.
In one possible implementation, the standard requirements include:
Boosting power plant: low-voltage side delta connection and high-voltage side Y connection;
500kV and 220kV transformer substation: low-pressure side delta connection and high-and medium-pressure side Y connection;
User load specialization: low-side delta connection and high-side Y connection.
In one possible implementation, after the phase shift loop is eliminated, the last p rows of the loop phase angle difference vector Δθ * are all 0.
According to the digital twin model construction method of the power transmission network transformer, provided by the embodiment of the invention, the digital twin model is constructed, the phase shifting ring in the power transmission network digital twin model is searched, the phase shifting ring is eliminated, the correctness of the power transmission network digital twin model is ensured, and the problem that manual searching is time-consuming and labor-consuming in the prior art is solved.
Drawings
In order to more clearly illustrate the technical solutions of the exemplary embodiments of the present invention, the drawings that are needed in the examples will be briefly described below, it being understood that the following drawings only illustrate some examples of the present invention and therefore should not be considered as limiting the scope, and that other related drawings may be obtained from these drawings without inventive effort for a person skilled in the art. In the drawings:
Fig. 1 is a flowchart of a method for constructing a digital twin model of a power transmission network transformer according to an embodiment of the present invention.
Fig. 2 is a topology diagram of an example system according to an embodiment of the present invention.
Fig. 3 is a schematic diagram of a digital twin model of a two-winding transformer according to an embodiment of the present invention.
Fig. 4 is a schematic diagram of a digital twin model of a three-winding transformer according to an embodiment of the present invention.
Detailed Description
For the purpose of making apparent the objects, technical solutions and advantages of the present invention, the present invention will be further described in detail with reference to the following examples and the accompanying drawings, wherein the exemplary embodiments of the present invention and the descriptions thereof are for illustrating the present invention only and are not to be construed as limiting the present invention.
Example 1
As shown in fig. 1, a digital twin model construction method of a power transmission network transformer includes:
S1, obtaining short-circuit impedance parameters of each transformer in a power transmission network, and constructing a digital twin model of the power transmission network transformer based on the short-circuit impedance parameters of each transformer.
The transformers in the power transmission network are divided into a two-winding transformer and a three-winding transformer, so that the number of windings of the transformers needs to be determined before the short-circuit impedance parameters of the transformers are determined, and the digital twin model of the transformers is constructed more accurately.
S2, constructing an adjacency matrix and branch line transfer item vectors based on a digital twin model of the transmission network transformer, and determining a phase shift loop verification result in the digital twin model according to the adjacency matrix and the branch line transfer item vectors, wherein the phase shift loop verification result comprises no phase shift loop or a phase shift loop.
Alternatively, each node and branch needs to be assigned a continuous and unique number, based on which an adjacency matrix and branch term vectors can be constructed, so that the branch with the phase shift loop can be determined.
And S3, if the phase-shifting ring verification result is that the phase-shifting ring exists, positioning the phase-shifting ring to obtain a phase-shifting ring positioning result.
After the branch with the phase-shifting ring is determined, the fault of the connection method of the transformers can be rapidly determined, so that the connection method of the transformers is corrected, the phase-shifting ring is eliminated, and the digital twin model of the power transmission network is constructed more accurately.
S4, correcting a transformer connection method on a branch in the phase-shifting ring according to the positioning result of the phase-shifting ring so as to eliminate the phase-shifting ring and obtain a final digital twin module.
In one possible embodiment, obtaining a short circuit impedance parameter for each transformer in a power transmission network includes:
A1, acquiring the type of the transformer and nameplate parameters, judging whether the type of the transformer is a two-winding transformer, if so, entering a step A2, otherwise, judging that the transformer is a three-winding transformer, and entering a step A3.
Optionally, the type of the transformer and the nameplate parameters may be data pre-stored in a database, or may be data generated by man-machine interaction.
A2, determining short-circuit impedance parameters of the two-winding transformer according to nameplate parameters of the transformer, wherein the short-circuit impedance parameters are as follows:
Wherein U 0% represents a transformer short-circuit voltage percentage, P 0 represents a transformer short-circuit loss, U B represents a transformer high-voltage side rated voltage, S B represents a transformer high-voltage side rated capacity, X represents a transformer winding leakage reactance, and R represents a transformer winding resistance.
A3, determining short-circuit impedance parameters of the three-winding transformer according to nameplate parameters of the transformer, wherein the short-circuit impedance parameters are as follows:
Wherein X 1 represents a transformer high voltage winding leakage reactance, X 2 represents a transformer medium voltage winding leakage reactance, X 3 represents a transformer low voltage winding leakage reactance, U 120% represents a transformer high-medium voltage short circuit voltage percentage, U 230% represents a transformer medium voltage-low voltage short circuit voltage percentage, U 130% represents a transformer high voltage-low voltage short circuit voltage percentage, U 1B represents a transformer high voltage side rated voltage, U 2B represents a transformer medium voltage side rated voltage, U 3B represents a transformer low voltage side rated voltage, S 1B represents a transformer high voltage side rated capacity, S 2B represents a transformer medium voltage side rated capacity, S 3B represents a transformer low voltage side rated capacity, P 120 represents a transformer high voltage-medium voltage short circuit loss, P 230 represents a transformer medium voltage-low voltage short circuit loss, P 130 represents a transformer high voltage-low voltage short circuit loss, R 1 represents a transformer high voltage winding resistance, R 2 represents a transformer medium voltage winding resistance, and R 3 represents a transformer low voltage winding resistance.
A4, obtaining short-circuit impedance parameters of each transformer according to the methods in the step A2 and the step A3.
Optionally, to ensure the stability of electromagnetic transient calculation, R 1,R2,R3 is more than or equal to 0, and if a negative value exists in the electromagnetic transient calculation, 0 is used for replacing the electromagnetic transient calculation.
In one possible implementation, a digital twinning model of a grid transformer is constructed based on short circuit impedance parameters of the respective transformers, comprising: and constructing a digital twin model of the power transmission network transformer in an equivalent substitution mode according to the short-circuit impedance parameters of each transformer.
As shown in fig. 3, for a two-winding transformer, it is replaced with R, X series branches; as shown in fig. 4, for the three-winding transformer, a virtual neutral node is added, and is replaced by a three-port network formed by R 1、R2、R3、X1、X2、X3.
In one possible implementation, the adjacency matrix a is:
Wherein A (i) represents the ith row of the adjacent matrix, i is more than or equal to 1 and less than or equal to b, and b is the number of power network branches; corresponding branch i, k is the starting node of branch i, m represents the ending node of branch i, i.e. the ith row A (i: of matrix A) of which the elements except k columns and m columns are all 0; n is the number of electrical nodes of the power network.
Where i is the branch number, and assuming that it is a branch from k to m, the phase difference caused by the transformer connection is Δθ i: for the Y-delta connection, the phase difference is 30.
In one possible implementation, the branch term vector Δθ is:
Δθ＝[Δθ1 Δθ2 … Δθm]T
Wherein, delta theta i represents the phase difference of the ith branch caused by the connection of the transformer, T represents the transpose.
In one possible implementation manner, positioning the phase shift ring to obtain a positioning result of the phase shift ring includes:
adopting Gaussian elimination method to process the adjacent matrix A, converting the adjacent matrix A into an upper triangular matrix U through series elementary line transformation, and obtaining the following steps:
KA＝U
Where K represents the invertible matrix.
The method for obtaining the reversible matrix K can be as follows:
1. initializing B: let B 1 be the B-order identity matrix:
a ij is the ith row and jth column element of A (for any 1.ltoreq.i.ltoreq.b, 1.ltoreq.j.ltoreq.n)
2. And (3) circulation: for i continuously varying from 1 to n
2.1 Judging whether a ii、ai+1i…、abi is all 0:
2.1.1 if yes, jumping to the next cycle of i;
2.1.2 if no:
2.1.2.1 swaps the non-zero elements to the a ii position by way of the i-B row swap of B.
2.1.2.2 Cycles: for j continuously varying from i+1 to n
Wherein B (i,:), B (j,:), and B '(i,:) are respectively the ith row of matrix B, the jth row of matrix B, and the ith row of matrix B'.
3. B= [ U K ], where U is the matrix of the first n columns of transformed matrix B (this is the upper triangular matrix according to the gaussian elimination step), and K is the matrix of the last B columns of transformed matrix B.
According to the branch term transfer vector, the loop phase angle difference vector delta theta * is obtained as follows:
Δθ*＝KΔθ
Wherein Δθ represents a branch term vector;
consider all rows of 0 in matrix U (assuming that there are no rows of 0, then the network is shown to have no loops, and no phase shift loop is required), assuming p rows at the end.
Determining rows which are all 0 in the upper triangular matrix U to obtain an end p rows; i.e. the row where U is all 0 is determined to determine the value of q, and then check how many non-0 elements there are in the last q row in Δθ *. And when no row which is 0 exists in the upper triangular matrix U, determining that no loop exists, namely the positioning result of the phase-shifting loop is no phase-shifting loop, and otherwise, determining the phase-shifting loop.
Judging whether the last p rows in the loop phase angle difference vector delta theta * are all 0, if yes, determining that no loop exists, namely the positioning result of the phase-shifting loop is that no phase-shifting loop exists, otherwise, determining that non-zero elements exist.
And determining that the number of non-zero elements in the last p rows of the loop phase angle difference vector delta theta * is q, and determining that q phase shift loops exist as a phase shift loop check result.
In one possible implementation manner, positioning the phase shift ring to obtain a positioning result of the phase shift ring includes:
And determining the row where q non-zero elements are located in the last p rows of the loop phase angle difference vector delta theta *, and obtaining the target row number.
And determining the columns of non-zero elements in the row corresponding to the target row number in the reversible matrix K to obtain the target column number.
And taking the branch number corresponding to the target column number as a branch for forming the phase shifting ring to obtain a phase shifting ring positioning result.
In one possible implementation manner, according to the positioning result of the phase shift ring, correcting the transformer connection on the branch in the phase shift ring to eliminate the phase shift ring, including:
Determining whether a transformer connection corresponding to a branch circuit forming the phase-shifting ring in the phase-shifting ring positioning result meets standard requirements, if so, determining that the transformer connection is not required to be corrected, otherwise, determining that the transformer connection is required to be corrected.
And obtaining all target transformers which need to be corrected in the branches forming the phase-shifting ring, correcting the transformer connection method of the target transformers, and eliminating the phase-shifting ring.
Optionally, when the transformer connection corresponding to the branch circuit forming the phase-shifting ring does not meet the standard requirement, the corresponding transformer connection is modified to the connection corresponding to the standard requirement, so as to eliminate the phase-shifting ring.
In one possible implementation, the standard requirements include:
boosting power plant: low-side delta connection and high-side Y connection.
500KV and 220kV transformer substation: low-side delta connection and high-and medium-side Y connection.
User load specialization: low-side delta connection and high-side Y connection.
In one possible implementation, after the phase shift loop is eliminated, the last p rows of the loop phase angle difference vector Δθ * are all 0.
According to the digital twin model construction method of the power transmission network transformer, provided by the embodiment of the invention, the digital twin model is constructed, the phase shifting ring in the power transmission network digital twin model is searched, the phase shifting ring is eliminated, the correctness of the power transmission network digital twin model is ensured, and the problem that manual searching is time-consuming and labor-consuming in the prior art is solved.
Example 2
The embodiment discloses an automatic verification and correction method for a digital twin model of a transformer, which comprises the following steps: (1) Calculating the short circuit impedance according to the nameplate parameters of the transformer; (2) constructing an adjacency matrix and branch phase shift vectors; (3) processing the adjacency matrix by adopting a Gaussian elimination method; (4) positioning a phase shift ring; (5) For the branches constituting the phase-shifting loop, a transformer connection which may have errors is corrected.
1) Calculating the short circuit impedance according to the nameplate parameters of the transformer
For a two-winding transformer (reduced to the high-voltage side):
U 0% -the short-circuit voltage percentage of the transformer;
P 0 -transformer short-circuit loss;
u B -the high side rated voltage of the transformer;
s B -high-side rated capacity of the transformer;
X-transformer winding leakage reactance (Ω);
r-transformer winding resistance (Ω);
For a three-winding transformer (reduced to three sides):
U 120% -percent transformer high-medium voltage short circuit voltage; u 230% -the percentage of medium-voltage-low-voltage short-circuit voltage of the transformer; u 130% -the percentage of high-voltage-low-voltage short-circuit voltage of the transformer; p 120 -transformer high-medium voltage short-circuit loss;
P 230 -transformer medium-low voltage short circuit loss;
P 130 -transformer high-low voltage short-circuit loss;
U 1B -the high side rated voltage of the transformer;
u 2B -the voltage rating on the medium voltage side of the transformer;
u 3B -the rated voltage of the low-voltage side of the transformer;
s 1B -high-side rated capacity of the transformer;
s 2B -rated capacity of the medium voltage side of the transformer;
S 3B -the rated capacity of the low-voltage side of the transformer;
X 1 —transformer high voltage winding leakage reactance (Ω);
x 2 —the transformer medium voltage winding leakage reactance (Ω);
X 3 —transformer low voltage winding leakage reactance (Ω);
r 1 -high voltage winding resistance of transformer (Ω);
R 2 -the voltage winding resistance (Ω) in the transformer;
R 3 -the low-voltage winding resistance (Ω) of the transformer.
In order to ensure the calculation stability of the electromagnetic transient state, R 1,R2,R3 is more than or equal to 0, and if a negative value exists in the electromagnetic transient state, 0 is used for replacing the electromagnetic transient state.
2) Constructing an adjacency matrix A and a branch phase shift vector delta theta;
Each node, branch, needs to be assigned a continuous, unique number.
Let the number of electrical nodes be n and the number of branches be b.
Δθ＝[Δθ1 Δθ2 … Δθm]T
Where i is the branch number, and assuming that it is a branch from k to m, the phase difference caused by the transformer connection is Δθ i. For the Y-delta connection, the phase difference is 30.
3) Adopting Gaussian elimination method to process the adjacent matrix A, and converting the adjacent matrix A into an upper triangular matrix U through series elementary line transformation, namely:
KA＝U
k is the invertible matrix.
Defining a loop phase angle difference vector delta theta *
Δθ*＝KΔθ
Consider all rows of 0 in matrix U (assuming that there are no rows of 0, then the network is shown to have no loops, and no phase shift loop is required), assuming p rows at the end.
If the last p rows of elements of delta theta * are all 0s, then there is no phase shift loop in the network.
If there are q non-zero elements in the last p rows of Δθ *, there are q phase-shifting loops in the network.
As shown in fig. 2, the embodiment of the present invention provides a topological connection diagram of a power grid system, and it can be seen that three phase shift loops exist.
Table 1 example system branch and node, branch phase shift relationship table
| Branch numbering | Head-end node | End node | Phase shift/deg |
| 1 | 1 | 2 | 0 |
| 2 | 2 | 3 | 0 |
| 3 | 4 | 5 | 0 |
| 4 | 4 | 6 | 0 |
| 5 | 2 | 11 | 0 |
| 6 | 4 | 11 | -30 |
| 7 | 8 | 11 | -30 |
| 8 | 3 | 12 | 0 |
| 9 | 9 | 12 | -30 |
| 10 | 5 | 12 | 0 |
| 11 | 6 | 7 | 30 |
| 12 | 6 | 7 | 0 |
| 13 | 1 | 10 | 30 |
| 14 | 1 | 10 | 0 |
According to table 1, an adjacency matrix a can be constructed, the branch phase shift vector delta theta is obtained by Gaussian elimination, and the matrix conditions of U, K and delta theta * are obtained:
the U matrix is shown in Table 2.
TABLE 2
The last 3 rows are all 0.
K. The Δθ * matrix is shown in table 3.
TABLE 3 Table 3
The algorithm finds 3 phase shift loops:
(1) Branch 11, branch 12.
(2) Branch 2, branch 3, branch 5, branch 6, branch 8, branch 10.
(3) Branch 13, branch 14.
According to the rule check, the transformer branch circuit to be corrected is: the specific correction results of branch 12, branch 6, and branch 14 are shown in table 4.
TABLE 4 Table 4
| Branch numbering | Head-end node | End node | Pre-correction phase shift/deg | Corrected phase shift/deg |
| 1 | 1 | 2 | 0 | 0 |
| 2 | 2 | 3 | 0 | 0 |
| 3 | 4 | 5 | 0 | 0 |
| 4 | 4 | 6 | 0 | 0 |
| 5 | 2 | 11 | 0 | 0 |
| 6 | 4 | 11 | -30 | 0 |
| 7 | 8 | 11 | -30 | -30 |
| 8 | 3 | 12 | 0 | 0 |
| 9 | 9 | 12 | -30 | -30 |
| 10 | 5 | 12 | 0 | 0 |
| 11 | 6 | 7 | 30 | 30 |
| 12 | 6 | 7 | 0 | 30 |
| 13 | 1 | 10 | 30 | 30 |
| 14 | 1 | 10 | 0 | 30 |
The corrected matrices of K, Δθ * are shown in table 5.
TABLE 5
No phase shift loop is present.
4) Positioning a phase shift ring:
Looking at the q non-zero elements of the last p rows of Δθ * are in the row: let k 1、…、kq be the line.
The branch numbers corresponding to the columns of the non-zero elements in the K j row in the matrix K are the branches forming the phase shift ring.
5) For branches forming a phase-shifting loop, a possibly erroneous transformer connection is corrected
Checking whether the branches in the phase shift loop meet the following requirements:
6-1) power plant boost variation: low-side delta connection and high-side Y connection.
6-2) 500KV, 220kV transformer substation: low-side delta connection and high-and medium-side Y connection.
6-3) User load specialization: low-side delta connection and high-side Y connection.
And adjusting the transformer branch which does not meet the requirements in the phase-shifting loop, eliminating the phase-shifting loop, and ensuring that the last p rows of delta theta * are all 0.
The foregoing description of the embodiments has been provided for the purpose of illustrating the general principles of the invention, and is not meant to limit the scope of the invention, but to limit the invention to the particular embodiments, and any modifications, equivalents, improvements, etc. that fall within the spirit and principles of the invention are intended to be included within the scope of the invention.
Claims (9)
1. The digital twin model construction method of the power transmission network transformer is characterized by comprising the following steps of:
obtaining short-circuit impedance parameters of each transformer in a power transmission network, and constructing a digital twin model of the power transmission network transformer based on the short-circuit impedance parameters of each transformer;
Constructing an adjacency matrix and branch transition term vectors based on a digital twin model of a power transmission network transformer, and determining a phase-shifting loop verification result in the digital twin model according to the adjacency matrix and the branch transition term vectors, wherein the phase-shifting loop verification result comprises no phase-shifting loop or a phase-shifting loop;
If the phase-shifting ring verification result is that the phase-shifting ring exists, positioning the phase-shifting ring to obtain a phase-shifting ring positioning result;
Correcting a transformer connection method on a branch in the phase-shifting ring according to the positioning result of the phase-shifting ring so as to eliminate the phase-shifting ring and obtain a final digital twin module;
Obtaining short-circuit impedance parameters of each transformer in a power transmission network comprises the following steps:
A1, acquiring the type of a transformer and nameplate parameters, judging whether the type of the transformer is a two-winding transformer, if so, entering a step A2, otherwise, judging that the transformer is a three-winding transformer, and entering a step A3;
A2, determining short-circuit impedance parameters of the two-winding transformer according to nameplate parameters of the transformer, wherein the short-circuit impedance parameters are as follows:
Wherein U 0% represents the transformer short-circuit voltage percentage, P 0 represents the transformer short-circuit loss, U B represents the transformer high-voltage side rated voltage, S B represents the transformer high-voltage side rated capacity, X represents the transformer winding leakage reactance, and R represents the transformer winding resistance;
a3, determining short-circuit impedance parameters of the three-winding transformer according to nameplate parameters of the transformer, wherein the short-circuit impedance parameters are as follows:
Wherein, X 1 represents a transformer high-voltage winding leakage reactance, X 2 represents a transformer medium-voltage winding leakage reactance, X 3 represents a transformer low-voltage winding leakage reactance, U 120% represents a transformer high-medium voltage short-circuit voltage percentage, U 230% represents a transformer medium-voltage low-voltage short-circuit voltage percentage, U 130% represents a transformer high-voltage low-voltage short-circuit voltage percentage, U 1B represents a transformer high-voltage side rated voltage, U 2B represents a transformer medium-voltage side rated voltage, U 3B represents a transformer low-voltage side rated voltage, S 1B represents a transformer high-voltage side rated capacity, S 2B represents a transformer medium-voltage side rated capacity, S 3B represents a transformer low-voltage side rated capacity, P 120 represents a transformer high-medium voltage short-circuit loss, P 230 represents a transformer medium-voltage low voltage short-circuit loss, P 130 represents a transformer high-voltage winding resistance, R 1 represents a transformer medium-voltage winding resistance, R 2 represents a transformer medium-voltage winding resistance, and R 3 represents a transformer low-voltage winding resistance;
a4, obtaining short-circuit impedance parameters of each transformer according to the methods in the step A2 and the step A3.
2. The method of constructing a digital twin model of a power transmission network transformer according to claim 1, wherein constructing the digital twin model of the power transmission network transformer based on the short-circuit impedance parameters of the respective transformers comprises:
And constructing a digital twin model of the power transmission network transformer in an equivalent substitution mode according to the short-circuit impedance parameters of each transformer.
3. The method for constructing a digital twin model of a grid transformer according to claim 2, wherein the adjacency matrix a is:
Wherein A (i) represents the ith row of the adjacent matrix, i is more than or equal to 1 and less than or equal to b, and b is the number of power network branches; corresponding branch i, k is the starting node of branch i, m represents the ending node of branch i, i.e. the ith row A (i: of matrix A) of which the elements except k columns and m columns are all 0; n is the number of electrical nodes of the power network.
4. A method of constructing a digital twin model of a grid transformer according to claim 3, wherein the branch term vector Δθ is:
Δθ＝[Δθ1 Δθ2 … Δθm]T
Wherein, delta theta i represents the phase difference of the ith branch caused by the connection of the transformer, T represents the transpose.
5. The method for constructing a digital twin model of a grid transformer according to claim 4, wherein positioning the phase shift ring to obtain a phase shift ring positioning result comprises:
adopting Gaussian elimination method to process the adjacent matrix A, converting the adjacent matrix A into an upper triangular matrix U through series elementary line transformation, and obtaining the following steps:
KA＝U
Wherein K represents a reversible matrix;
According to the branch term transfer vector, the loop phase angle difference vector delta theta * is obtained as follows:
Δθ*＝KΔθ
Wherein Δθ represents a branch term vector;
determining rows which are all 0 in the upper triangular matrix U to obtain an end p rows; when no row which is 0 exists in the upper triangular matrix U, determining that no loop exists, namely the positioning result of the phase-shifting loop is no phase-shifting loop, otherwise, determining the phase-shifting loop;
judging whether the last p rows in the loop phase angle difference vector delta theta * are all 0, if yes, determining that no loop exists, namely the positioning result of the phase shifting loop is that no phase shifting loop exists, otherwise, determining that non-zero elements exist;
and determining that the number of non-zero elements in the last p rows of the loop phase angle difference vector delta theta * is q, and determining that q phase shift loops exist as a phase shift loop check result.
6. The method for constructing a digital twin model of a grid transformer according to claim 5, wherein positioning the phase shift ring to obtain a phase shift ring positioning result comprises:
Determining the row where q non-zero elements are located in the last p rows of the loop phase angle difference vector delta theta * to obtain a target row number;
determining the columns of non-zero elements in the row corresponding to the target row number in the reversible matrix K to obtain the target column number;
and taking the branch number corresponding to the target column number as a branch for forming the phase shifting ring to obtain a phase shifting ring positioning result.
7. The method of constructing a digital twin model of a grid transformer according to claim 6, wherein correcting the transformer connection on the branch in the phase shift loop to eliminate the phase shift loop according to the positioning result of the phase shift loop comprises:
determining whether a transformer connection corresponding to a branch circuit forming the phase-shifting ring in the phase-shifting ring positioning result meets standard requirements, if so, determining that the transformer connection is not required to be corrected, otherwise, determining that the transformer connection is required to be corrected;
And obtaining all target transformers which need to be corrected in the branches forming the phase-shifting ring, correcting the transformer connection method of the target transformers, and eliminating the phase-shifting ring.
8. The method of constructing a digital twinning model of a power transmission network transformer of claim 7, wherein the standard requirements include:
Boosting power plant: low-voltage side delta connection and high-voltage side Y connection;
500kV and 220kV transformer substation: low-pressure side delta connection and high-and medium-pressure side Y connection;
User load specialization: low-side delta connection and high-side Y connection.
9. The method of constructing a digital twinning model for a power transmission network transformer according to claim 7, wherein the last p rows of the loop phase angle difference vector Δθ * are all 0 after the phase shift loop is eliminated.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202311197441.0A CN117217008B (en) | 2023-09-15 | 2023-09-15 | Digital twin model construction method of power transmission network transformer |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202311197441.0A CN117217008B (en) | 2023-09-15 | 2023-09-15 | Digital twin model construction method of power transmission network transformer |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN117217008A CN117217008A (en) | 2023-12-12 |
| CN117217008B true CN117217008B (en) | 2024-05-07 |
Family
ID=89036698
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202311197441.0A Active CN117217008B (en) | 2023-09-15 | 2023-09-15 | Digital twin model construction method of power transmission network transformer |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN117217008B (en) |
Citations (11)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN112306658A (en) * | 2020-10-31 | 2021-02-02 | 贵州电网有限责任公司 | Digital twin application management scheduling method for multi-energy system |
| CN112685949A (en) * | 2020-11-25 | 2021-04-20 | 江苏科技大学 | Transformer health prediction method based on digital twinning |
| CN113591173A (en) * | 2021-07-30 | 2021-11-02 | 贵州电网有限责任公司 | Data visualization interaction method for multi-energy system digital twin |
| CN114217571A (en) * | 2021-12-15 | 2022-03-22 | 苏州华光智控电子科技有限公司 | Laser cutting follow-up system health monitoring method based on digital twinning |
| CN114595528A (en) * | 2022-02-28 | 2022-06-07 | 国网陕西省电力有限公司电力科学研究院 | Visualization system and method of digital twin-based transformer outgoing line device |
| CN114744764A (en) * | 2022-04-25 | 2022-07-12 | 南方电网科学研究院有限责任公司 | Digital twin terminal, system and terminal control method for power distribution network |
| CN115224698A (en) * | 2022-07-20 | 2022-10-21 | 南京理工大学 | Reactive power-voltage optimization control method for new energy power system based on digital twinning |
| CN115310227A (en) * | 2022-08-09 | 2022-11-08 | 重庆大学 | Method for constructing gear power propagation digital twin model |
| CN115563897A (en) * | 2022-10-13 | 2023-01-03 | 南京电力设计研究院有限公司 | Transformer state evaluation method based on digital twinning |
| WO2023078566A1 (en) * | 2021-11-05 | 2023-05-11 | Fraunhofer-Gesellschaft zur Förderung der angewandten Forschung e.V. | Operating a wireless communication network using a digital twin |
| CN116224931A (en) * | 2023-02-10 | 2023-06-06 | 上海益实智能科技有限公司 | Real-time data driving-based digital twin simulation operation method for cigarette factory |
Family Cites Families (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109884459B (en) * | 2019-01-30 | 2020-12-01 | 国网浙江省电力有限公司电力科学研究院 | Intelligent online diagnosis and positioning method for winding deformation of power transformer |
| US11755689B2 (en) * | 2019-06-24 | 2023-09-12 | Intel Corporation | Methods, systems, articles of manufacture and apparatus to manage process excursions |
-
2023
- 2023-09-15 CN CN202311197441.0A patent/CN117217008B/en active Active
Patent Citations (11)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN112306658A (en) * | 2020-10-31 | 2021-02-02 | 贵州电网有限责任公司 | Digital twin application management scheduling method for multi-energy system |
| CN112685949A (en) * | 2020-11-25 | 2021-04-20 | 江苏科技大学 | Transformer health prediction method based on digital twinning |
| CN113591173A (en) * | 2021-07-30 | 2021-11-02 | 贵州电网有限责任公司 | Data visualization interaction method for multi-energy system digital twin |
| WO2023078566A1 (en) * | 2021-11-05 | 2023-05-11 | Fraunhofer-Gesellschaft zur Förderung der angewandten Forschung e.V. | Operating a wireless communication network using a digital twin |
| CN114217571A (en) * | 2021-12-15 | 2022-03-22 | 苏州华光智控电子科技有限公司 | Laser cutting follow-up system health monitoring method based on digital twinning |
| CN114595528A (en) * | 2022-02-28 | 2022-06-07 | 国网陕西省电力有限公司电力科学研究院 | Visualization system and method of digital twin-based transformer outgoing line device |
| CN114744764A (en) * | 2022-04-25 | 2022-07-12 | 南方电网科学研究院有限责任公司 | Digital twin terminal, system and terminal control method for power distribution network |
| CN115224698A (en) * | 2022-07-20 | 2022-10-21 | 南京理工大学 | Reactive power-voltage optimization control method for new energy power system based on digital twinning |
| CN115310227A (en) * | 2022-08-09 | 2022-11-08 | 重庆大学 | Method for constructing gear power propagation digital twin model |
| CN115563897A (en) * | 2022-10-13 | 2023-01-03 | 南京电力设计研究院有限公司 | Transformer state evaluation method based on digital twinning |
| CN116224931A (en) * | 2023-02-10 | 2023-06-06 | 上海益实智能科技有限公司 | Real-time data driving-based digital twin simulation operation method for cigarette factory |
Non-Patent Citations (2)
| Title |
|---|
| 基于Simulink配电变压器模型的建立及仿真;王希平;李文才;李燕;路文梅;;自动化技术与应用;20090925(09);第124-126, 132页 * |
| 基于电力系统全数字实时仿真装置的大电网机电暂态电磁暂态混合仿真;朱旭凯;周孝信;田芳;徐得超;黄柳强;李刚;;电网技术;20110305(03);第26-31页 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN117217008A (en) | 2023-12-12 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Haque | Efficient load flow method for distribution systems with radial or mesh configuration | |
| CN108429257B (en) | Distribution network frame planning method based on load point clustering partitioning and considering geographic factors | |
| Lin et al. | Heuristic rule-based phase balancing of distribution systems by considering customer load patterns | |
| CN112255567B (en) | Short-circuit current rapid determination method for power distribution network containing photovoltaic power supply | |
| CN108270248B (en) | Active safety correction control method based on quasi-steady-state sensitivity analysis | |
| CN110880764B (en) | Fault processing method for unbalanced distribution network containing inversion type distributed power supply | |
| CN105320810A (en) | Short-circuit current scanning optimization method | |
| Pujiantara et al. | Optimization technique based adaptive overcurrent protection in radial system with DG using genetic algorithm | |
| CN113363987B (en) | Master-slave fault self-healing control method and system for power distribution network | |
| CN117217008B (en) | Digital twin model construction method of power transmission network transformer | |
| CN111324997B (en) | Transformer three-phase model establishing method based on incidence matrix considering neutral point | |
| CN110427635B (en) | LCC-HVDC optimized electromagnetic transient simulation method | |
| CN110176755B (en) | Fault current limiter configuration method and device | |
| CN108595819B (en) | Transformer three-phase model establishing method based on conversion matrix | |
| CN116722583A (en) | Method and system for determining negative sequence current of doubly-fed wind turbine under grid short circuit fault | |
| CN115455737A (en) | Power system short-circuit current calculation method, device, terminal and medium | |
| CN108565866B (en) | Transformer three-phase model establishing method based on unified transformation matrix | |
| CN110968973B (en) | Simulation method, control method, electronic equipment and storage medium of transformer model | |
| CN111965478B (en) | Method and system for determining short-circuit current by considering phase shifting effect of transformer | |
| CN113991652A (en) | Data-driven multi-output calculation method for short-circuit current of IIDG-containing power distribution network | |
| Li et al. | Explicit linear model of bus impedance with unknown branch status | |
| CN111965477A (en) | Method and system for determining short-circuit current by considering zero sequence mutual inductance influence | |
| CN114928035B (en) | Generation method and system based on strategy for restraining short-circuit current by multiple limiters | |
| CN112803420B (en) | SCOPF model active power flow constraint reduction method based on full rank discrimination | |
| Chen et al. | Modelling asymmetrical HVDC transfer network for multi-vendor-multi-terminal interoperability |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |