CN117217008A - 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
- CN117217008A CN117217008A CN202311197441.0A CN202311197441A CN117217008A CN 117217008 A CN117217008 A CN 117217008A CN 202311197441 A CN202311197441 A CN 202311197441A CN 117217008 A CN117217008 A CN 117217008A
- Authority
- CN
- China
- Prior art keywords
- transformer
- phase
- voltage
- loop
- digital twin
- 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.)
- Granted
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 64
- 238000000034 method Methods 0.000 claims abstract description 38
- 239000013598 vector Substances 0.000 claims abstract description 37
- 238000012795 verification Methods 0.000 claims abstract description 16
- 230000010363 phase shift Effects 0.000 claims description 55
- 238000004804 winding Methods 0.000 claims description 54
- 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
- 238000012546 transfer Methods 0.000 abstract description 4
- 230000001052 transient effect Effects 0.000 description 7
- 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
- 238000004364 calculation method Methods 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 is 0 % represents the percent of short-circuit voltage of the transformer, P 0 Indicating transformer short-circuit loss, U B Represents the rated voltage of the high-voltage side of the transformer, S B The high-voltage side rated capacity of the transformer is represented, X represents the leakage reactance of the transformer winding, and R represents the resistance of the transformer winding;
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 is 1 Represents leakage reactance of high-voltage winding of transformer, X 2 Represents the leakage reactance of the medium-voltage winding of the transformer, X 3 Indicating leakage reactance of low-voltage winding of transformer, U 120 % represents the percentage of the short circuit voltage between the high voltage and the medium voltage of the transformer, U 230 % represents the percentage of the medium-voltage-low-voltage short-circuit voltage of the transformer, U 130 % represents the percentage of high-voltage-low-voltage short-circuit voltage of the transformer, U 1B Indicating the rated voltage of the high-voltage side of the transformer, U 2B Indicating the rated voltage of the medium voltage side of the transformer, U 3B Representing the rated voltage of the low-voltage side of the transformer, S 1B Representing the rated capacity of the high-voltage side of the transformer, S 2B Representing the rated capacity of the medium voltage side of the transformer, S 3B Representing the rated capacity of the low-voltage side of the transformer, P 120 Representing the high-voltage and medium-voltage short-circuit loss of the transformer, P 230 Representing the medium-voltage-low-voltage short-circuit loss of the transformer, P 130 Representing high-voltage-low-voltage short-circuit loss of transformer, R 1 R represents the resistance of the high-voltage winding of the transformer 2 Representing the resistance of the medium-voltage winding of the transformer, R 3 Representing the low-voltage winding resistance of the transformer;
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 Indicating the phase difference of the ith branch caused by transformer connection,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 shift vector, a loop phase angle difference vector delta theta is obtained * The method comprises the following steps:
Δθ * ＝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;
determining a loop phase angle difference vector delta theta * Whether the last p rows are all 0 or not, if yes, determining no loop, namely, determining that the positioning result of the phase shift loop is no phase shift loop, otherwise, determining non-zero elements;
determining a loop phase angle difference vector delta theta * If the number of non-zero elements in the last p rows of (1) is q, then determining that q phase shift loops exist as a result of phase shift loop verification.
In one possible implementation manner, positioning the phase shift ring to obtain a positioning result of the phase shift ring includes:
determining a loop phase angle difference vector delta theta * Q non-zero element rows in the last p rows 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-shifting loop is eliminated, the loop phase angle difference vector Δθ * All 0 s in the last p rows of (a).
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 is 0 % represents the percent of short-circuit voltage of the transformer, P 0 Indicating transformer short-circuit loss, U B Represents the rated voltage of the high-voltage side of the transformer, S B The rated capacity of the high-voltage side of the transformer is represented, X represents the leakage reactance of the winding of the transformer, and R represents the variationVoltage 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 is 1 Represents leakage reactance of high-voltage winding of transformer, X 2 Represents the leakage reactance of the medium-voltage winding of the transformer, X 3 Indicating leakage reactance of low-voltage winding of transformer, U 120 % represents the percentage of the short circuit voltage between the high voltage and the medium voltage of the transformer, U 230 % represents the percentage of the medium-voltage-low-voltage short-circuit voltage of the transformer, U 130 % represents the percentage of high-voltage-low-voltage short-circuit voltage of the transformer, U 1B Indicating the rated voltage of the high-voltage side of the transformer, U 2B Indicating the rated voltage of the medium voltage side of the transformer, U 3B Representing the rated voltage of the low-voltage side of the transformer, S 1B Representing the rated capacity of the high-voltage side of the transformer, S 2B Representing the rated capacity of the medium voltage side of the transformer, S 3B Representing the rated capacity of the low-voltage side of the transformer, P 120 Representing the high-voltage and medium-voltage short-circuit loss of the transformer, P 230 Representing the medium-voltage-low-voltage short-circuit loss of the transformer, P 130 Representing high-voltage-low-voltage short-circuit loss of transformer, R 1 R represents the resistance of the high-voltage winding of the transformer 2 Representing the resistance of the medium-voltage winding of the transformer, R 3 Representing the 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 electromagnetic transient calculation stability, R is required to be satisfied 1 ,R 2 ,R 3 0 or more, if negative values are present therein, 0 is substituted.
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. 3For a two-winding transformer, it is replaced with a R, X series branch; as shown in fig. 4, for a three-winding transformer, a virtual neutral point node is added and used as R 1 、R 2 、R 3 、X 1 、X 2 、X 3 The three-port network is constructed instead.
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 due to 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 Indicating the phase difference of the ith branch caused by transformer connection,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 For the b-order identity matrix:
a ij is the element of the ith row and the jth column 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 a ii 、a i+1i …、a bi Whether all are 0:
2.1.1 if yes, jumping to the next cycle of i;
2.1.2 if no:
2.1.2.1 non-zero elements are swapped to a by way of the i-B row swap of B ii Position.
2.1.2.2 cycle: 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 a matrix of the first n columns of the transformed matrix B (this is the upper triangular matrix according to the gaussian elimination step), and K is a matrix of the last B columns of the transformed matrix B.
According to the branch term shift vector, a loop phase angle difference vector delta theta is obtained * The method comprises the following steps:
Δθ * ＝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 Δθ is checked * How many non-0 elements there are in the last q rows in (a). When on the shelfAnd if no row which is 0 exists in the triangular matrix U, determining that no loop exists, namely, determining that the phase-shifting loop positioning result is no phase-shifting loop, and otherwise, determining that the phase-shifting loop exists.
Determining a loop phase angle difference vector delta theta * And (3) determining whether the last p rows are all 0, if so, determining that no loop exists, namely, determining that the positioning result of the phase-shifting loop is that no phase-shifting loop exists, and otherwise, determining that non-zero elements exist.
Determining a loop phase angle difference vector delta theta * If the number of non-zero elements in the last p rows of (1) is q, then determining that q phase shift loops exist as a result of phase shift loop verification.
In one possible implementation manner, positioning the phase shift ring to obtain a positioning result of the phase shift ring includes:
determining a loop phase angle difference vector delta theta * The row of q non-zero elements in the last p rows of (a) gets 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-shifting loop is eliminated, the loop phase angle difference vector Δθ * All 0 s in the last p rows of (a).
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 percent-transformer short circuit voltage;
P 0 -transformer short-circuit loss;
U B -the high-side rated voltage of the transformer;
S B -the 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 of transformer high-medium voltage short circuit voltage; u (U) 230 Percent of transformer medium-voltage-low-voltage short-circuit voltage; u (U) 130 Percent of high-low short circuit voltage of the transformer; p (P) 120 -transformer high-medium voltage short-circuit loss;
P 230 -transformer medium-low voltage short-circuit losses;
P 130 -transformer high-low voltage short-circuit loss;
U 1B -the high-side rated voltage of the transformer;
U 2B -voltage rating of the medium voltage side of the transformer;
U 3B -the low-side rated voltage of the transformer;
S 1B -the high-side rated capacity of the transformer;
S 2B -medium voltage side rated capacity of the transformer;
S 3B -the low-side rated capacity 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 -a transformer high voltage winding resistance (Ω);
R 2 -the voltage winding resistance (Ω) in the transformer;
R 3 -transformer low voltage winding resistance (Ω).
To ensure the electromagnetic transient calculation stability, R needs to be satisfied 1 ,R 2 ,R 3 0 or more, if negative values are present therein, 0 is substituted.
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 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 delta theta * The last p rows of elements are all 0, and no phase shift loop exists in the network.
If delta theta * The last p rows have q non-zero elements, and q phase shift loops are present 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, a branch can be constructedThe path phase shift vector delta theta is subjected to Gaussian elimination to obtain U, K and delta theta * Matrix case:
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 |
Corrected K, delta theta * The matrix is shown in table 5.
TABLE 5
No phase shift loop is present.
4) Positioning a phase shift ring:
look at Δθ * The q non-zero elements of the last p rows are in the row: let k be 1 、…、k q And (3) row.
The kth of the matrix K j The branch numbers corresponding to the columns where the non-zero elements of the rows are located are the branches forming the phase shifting 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.
The transformer branch circuit which does not meet the requirements in the phase-shifting ring is adjusted, the phase-shifting ring is eliminated, and delta theta is ensured * The last p rows 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 (10)
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;
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.
2. The method of constructing a digital twin model of a power transmission network transformer according to claim 1, wherein obtaining short-circuit impedance parameters of each transformer in the power transmission network comprises:
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 is 0 % represents the percent of short-circuit voltage of the transformer, P 0 Indicating a short circuit of the transformerLoss, U B Represents the rated voltage of the high-voltage side of the transformer, S B The high-voltage side rated capacity of the transformer is represented, X represents the leakage reactance of the transformer winding, and R represents the resistance of the transformer winding;
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 is 1 Represents leakage reactance of high-voltage winding of transformer, X 2 Represents the leakage reactance of the medium-voltage winding of the transformer, X 3 Indicating leakage reactance of low-voltage winding of transformer, U 120 % represents the percentage of the short circuit voltage between the high voltage and the medium voltage of the transformer, U 230 % represents the percentage of the medium-voltage-low-voltage short-circuit voltage of the transformer, U 130 % represents the percentage of high-voltage-low-voltage short-circuit voltage of the transformer, U 1B Indicating the rated voltage of the high-voltage side of the transformer, U 2B Indicating the rated voltage of the medium voltage side of the transformer, U 3B Representing the rated voltage of the low-voltage side of the transformer, S 1B Representing the rated capacity of the high-voltage side of the transformer, S 2B Representing the rated capacity of the medium voltage side of the transformer, S 3B Representing the rated capacity of the low-voltage side of the transformer, P 120 Representing the high-voltage and medium-voltage short-circuit loss of the transformer, P 230 Representing the medium-voltage-low-voltage short-circuit loss of the transformer, P 130 Representing high-voltage-low-voltage short-circuit loss of transformer, R 1 R represents the resistance of the high-voltage winding of the transformer 2 Representing the resistance of the medium-voltage winding of the transformer, R 3 Representing the low-voltage winding resistance of the transformer;
a4, obtaining short-circuit impedance parameters of each transformer according to the methods in the step A2 and the step A3.
3. The method of constructing a digital twin model of a power transmission network transformer according to claim 2, 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.
4. A method of constructing a digital twin model of a grid transformer according to claim 3, in which 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.
5. The method of constructing a digital twinning model for a power transmission network transformer of claim 4, wherein the branch term vector Δθ is:
Δθ＝[Δθ 1 Δθ 2 … Δθ m ] T
wherein, delta theta i Indicating the phase difference of the ith branch caused by transformer connection,t represents the transpose.
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:
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 shift vector, a loop phase angle difference vector delta theta is obtained * The method comprises the following steps:
Δθ * ＝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;
determining a loop phase angle difference vector delta theta * Whether the last p rows are all 0 or not, if yes, determining no loop, namely, determining that the positioning result of the phase shift loop is no phase shift loop, otherwise, determining non-zero elements;
determining a loop phase angle difference vector delta theta * If the number of non-zero elements in the last p rows of (1) is q, then determining that q phase shift loops exist as a result of phase shift loop verification.
7. The method for constructing a digital twin model of a grid transformer according to claim 6, wherein positioning the phase shift ring to obtain a phase shift ring positioning result comprises:
determining a loop phase angle difference vector delta theta * Q non-zero element rows in the last p rows 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.
8. The method of constructing a digital twin model of a grid transformer according to claim 7, 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.
9. The method of constructing a digital twin model of a grid transformer of claim 8, 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.
10. The method of constructing a digital twin model of a grid transformer according to claim 8, wherein after eliminating the phase shift loop, a loop phase angle difference vector Δθ * All 0 s in the last p rows of (a).
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 true CN117217008A (en) | 2023-12-12 |
| CN117217008B 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 (13)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20190340843A1 (en) * | 2019-06-24 | 2019-11-07 | Intel Corporation | Methods, systems, articles of manufacture and apparatus to manage process excursions |
| US20200241066A1 (en) * | 2019-01-30 | 2020-07-30 | State Grid Zhejiang Electric Power Company Limited Electric Power Research Institute | Intelligent on-line diagnosis and positioning method for winding deformation of power transformers |
| 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 |
-
2023
- 2023-09-15 CN CN202311197441.0A patent/CN117217008B/en active Active
Patent Citations (13)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20200241066A1 (en) * | 2019-01-30 | 2020-07-30 | State Grid Zhejiang Electric Power Company Limited Electric Power Research Institute | Intelligent on-line diagnosis and positioning method for winding deformation of power transformers |
| US20190340843A1 (en) * | 2019-06-24 | 2019-11-07 | Intel Corporation | Methods, systems, articles of manufacture and apparatus to manage process excursions |
| 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 |
|---|
| 朱旭凯;周孝信;田芳;徐得超;黄柳强;李刚;: "基于电力系统全数字实时仿真装置的大电网机电暂态电磁暂态混合仿真", 电网技术, no. 03, 5 March 2011 (2011-03-05), pages 26 - 31 * |
| 王希平;李文才;李燕;路文梅;: "基于Simulink配电变压器模型的建立及仿真", 自动化技术与应用, no. 09, 25 September 2009 (2009-09-25), pages 124 - 126 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN117217008B (en) | 2024-05-07 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Haque | Efficient load flow method for distribution systems with radial or mesh configuration | |
| Lin et al. | Heuristic rule-based phase balancing of distribution systems by considering customer load patterns | |
| CN100367620C (en) | Power network topology error identification method based on mixed state estimation | |
| CN112255567B (en) | Short-circuit current rapid determination method for power distribution network containing photovoltaic power supply | |
| CN110880764B (en) | Fault processing method for unbalanced distribution network containing inversion type distributed power supply | |
| Quintero-Duran et al. | A modified backward/forward sweep-based method for reconfiguration of unbalanced distribution networks. | |
| Pujiantara et al. | Optimization technique based adaptive overcurrent protection in radial system with DG using genetic algorithm | |
| CN117217008B (en) | Digital twin model construction method of power transmission network transformer | |
| CN105320810A (en) | Short-circuit current scanning optimization method | |
| CN111324997B (en) | Transformer three-phase model establishing method based on incidence matrix considering neutral point | |
| CN110176755B (en) | Fault current limiter configuration method and device | |
| CN110427635B (en) | LCC-HVDC optimized electromagnetic transient simulation method | |
| CN110208634B (en) | Method for acquiring asymmetric short-circuit current direct-current component of complex power system | |
| CN111651889A (en) | High-risk event screening method, device, equipment and storage medium for power transmission system | |
| CN110968973B (en) | Simulation method, control method, electronic equipment and storage medium of transformer model | |
| CN108565866B (en) | Transformer three-phase model establishing method based on unified transformation matrix | |
| Medina et al. | A three phase load flow algorithm for shipboard power systems (SPS) | |
| CN113991652A (en) | Data-driven multi-output calculation method for short-circuit current of IIDG-containing power distribution network | |
| CN114928035B (en) | Generation method and system based on strategy for restraining short-circuit current by multiple limiters | |
| Chen et al. | Modelling asymmetrical HVDC transfer network for multi-vendor-multi-terminal interoperability | |
| CN111965478A (en) | Method and system for determining short-circuit current by considering transformer phase shifting effect | |
| CN112803420B (en) | SCOPF model active power flow constraint reduction method based on full rank discrimination | |
| CN111965477A (en) | Method and system for determining short-circuit current by considering zero sequence mutual inductance influence | |
| Liu et al. | DC Fault Current Calculation and Fault Level Analysis in MMC-MVDC System | |
| CN112287500B (en) | Grid planning method for power distribution network in grid based on optimal cutting |
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 |