Three-phase distribution transformer connections modeling based on matrix operation method by phase-coordinates

This paper proposes a matrix operation method for modeling the three-phase transformer by phase-coordinates. Based on decoupling theory, the 12 × 12 dimension primitive admittance matrix is obtained at first employing the coupling configuration of the windings. Under the condition of asymmetric magnetic circuits, according to the boundary conditions for transformer connections, the transformers in different connections enable to be modeling by the matrix operation method from the primitive admittance matrix. Another purpose of this paper is to explain the differences of the phase-coordinates and the positive sequence parameters in the impedances of the transformers. The numerical testing results in IEEE-4 system show that the proposed method is valid and efficient.

proposed in [4][5][6][7][8][9][10]. Reference [4] developed an approach from the KCL and KVL, which was able to generate the 6 × 6 matrix in different transformer connections according to single-phase transformer symmetrical lattice equivalent circuits as the units grouping up. In the paper, the authors described the transformer in model relationship between the phase-coordinates and component coordinates as well. However, the interphase coupling did not be considered in this approach. Later, the improved models proposed in [5][6][7] derived from this approach of assembling singlephase transformer equivalent circuits by different connections. Another representative approach generated a primitive matrix by the six coins equivalent circuit of the transformer in (YN, yn0) connection [8,9].The main characteristic of this approach was that the models were obtained from product relations between the primitive matrix and the incidence matrix in different connections, while it was not accurate description that the left incidence matrix and the right one were same. Reference [10] accounted for a method of handling matrix singularity in the use of the transformer power-flow models. And the modified augmented nodal analysis (MANA) [11] was proposed, which enriched the model application of transformers.
There were several flaws in the previous approaches for modeling the three-phase transformers by phase-coordinates. For one thing, the phase self-impedances and mutual-impedances come from the transformer positive impedances directly, though the phase impedances are much closed to the positive impedances. For another, the most of models described the parameters of 6 coins by 6 × 6 dimension to 7 × 7 dimension matrix in the models for transformer, only consider the injection currents at a-b-c phase between the primary and the secondary sides, but not all the currents. In this case, the models did not enable to cover the use of both in power-flow and short circuit calculations under the asymmetric magnetic circuits.
The three-phase AC transmission theory used to be generated from the single-phase models of the system equipment by the AC circuit theory based on the symmetrical characteristics of the three-phase voltages and currents, which the three-phase symmetry is the characteristic of traditional power system analysis. And the unavoidable asymmetry conditions of the transformers are even more serious than the lines. The purpose of this paper is to model the three-phase transformers in different connections by phasecoordinates by matrix operation method. The main work for studying is show as follows: • Generated the primitive admittance matrix from the coupling windings by decoupling method. • Modeling the transformers in different connections by matrix operation method based on the asymmetric magnetic circuits. • Method for obtaining phase-coordinate modified parameters for the models of the transformers.

Methods/experimental
The aim of this paper is to solve the problem for winding connections modeling of three-phase transformer by phase-coordinates based on matrix operation method.
Firstly, we introduced the steps to obtain transformer admittance matrix by matrix operation method. Then we analyzed the modeling for the connections of three-phase distribution transformer, and also we analyzed the differences of the impedance parameters between phase-coordinates and sequence-coordinates. Finally, we verified the effectiveness of the modeling methods by simulation.

Modeling methodology
In this section, the modeling approach for a transformer is described by the matrix operation method from the 12 × 12 dimension primitive admittance matrix. The complex variables and the values of parameters are given in per-unit system. Firstly, the coupling configuration used prefers to describe and analyze the coupling phenomenon by comparing with the two circuit topologies for a single-phase double-winding transformer. And then the phase-coordinates construction methodology for the three-phase model is described by the following steps.
• Definition of the transformer primitive admittance matrix Y P . • Definition of the transformer admittance matrix Y T by the matrix operation method in different connections. Figure 1 shows the main steps of the conceptual scheme of the construction methodology [12]. Firstly, the primitive matrix is generated from the coupling windings shown as the inner blue frame. And then, according to the boundary conditions of the connections at the windings, the transformers are modeled by the matrix operation method.

Primitive admittance matrix Y P
The two kinds of single-phase equivalent circuits are shown in Fig. 2. Figure 2a is the T-configuration, and Fig. 2b is the coupled configuration. In Fig. 2, "A", "X" are the two buses at the primary side, "a", "x" are the two buses at the secondary side respectively. I 1 , I 2 (I A , I a ) are the Injection currents, and V 1 , V 2 (V A , V a ) stand for the node voltages at the windings. In Fig. 2a, Z 1 , Z 2 are the windings self-impedances, and jwM is the mutual-impedance between windings. w is the angular acceleration. In Fig. 2b, R A , R a , L A , L a are the resistances and self-inductances at the windings. M is the mutual-inductance. And g 0 + jb 0 is the noload admittance. The circuit relationship in Fig. 2a (in p.u system) can be given by (1) Fig. 2 Single-phase equivalent circuits of the double-winding transformer; a is the model of T-configuration, Z 1 and Z 2 are connected together, and The mutual impedance is connected between them. They form a T-shaped connection. b Model of Coupling configuration, which mutual impedance is generated by coupling There is potential difference between Bus "X" and "x" in transformer testing experiment. The T-configuration fails to express the electrical characteristics of the transformer, although the calculation enables to equipotential potential. When running in symmetrical operation, the two buses (X, x) are the same at "zero" potential point, so they can be connected in the form of equipotential. Generally, the two buses are not the same at the zero potential point when they run in three-phase asymmetric operation. In order to reflect the potential offset phenomenon of neutral point and consider the various connections of transformers, three-phase transformers for modeling should adopt the form of Fig. 2b.
Consequently, there is no electrical connection between Bus X and Bus x in Fig. 2b, which can indicate the potential difference between the two buses. The floating phenomenon without buses grounding and the different connections of the transformer enable to explain as well. While the T-configuration equivalent circuit of transformer in Fig. 2a utilizes the equipotential characteristics at the symmetrical operation, which is not suitable for asymmetrical operation analysis.
The application of Fig. 2b at Fig. 1, the three-phase transformer contains 6 coins and 12 buses (A-B-C buses and X-Y-Z buses at the primary side; a-b-c buses and x-y-z buses at the secondary side) shown as Fig. 3, and the generalized primitive model can be given by a 12 × 12 matrix using the decoupling methodology. The branch current equation can be given as Z AA Z AB Z AC Z Aa Z Ab Z Ac Z BA Z BB Z BC Z Ba Z Bb Z Bb Z CA Z CB Z CC Z Ca Z Cb Z Cc Z aA Z aB Z aC Z aa Z ab Z ac Z bA Z bB Z bC Z ba Z bb Z bc Z cA Z cB Z cC Z ca Z cb Z cc , which presents the impedance relationship of the transformer, and its admittance is y prim = Z −1 prim . The node voltage equation of the transformer can be shown as

Transformer admittance matrix Y T
According to different connections, we obtain the admittance models based on the matrix operation method from the derivation of the initial admittance matrix Y P . Figure 4 is used to describe the (YN, d11) connection of transformer. According to potential relations, we obtains The Eq. (4) presents the boundary conditions of the transformer in YN,d11 connection. The admittance matrix can be calculated by the matrix operation method according to the connected relationship of the transformer in YN,d11 connection.
Steps to obtain transformer admittance matrix by matrix operation method can be shown as follows: (1) YN connection of primary side: Bus X, Y and Z are grounding, and the rows and columns of them are retained in the primitive admittance matrix Y P . (2) d11 connection of secondary side: (3) Block 1: Unchanged (the changes of secondary winding connection do not affect the array of the primary side at Y P ); (4) Block 2: For the ja column in Y T , sum with the elements of the j.jy column, (i.e. j.ja = j.ja + j.jy ) in Y P ; (3) Non-diagonal elements in Y P : a-y: column element: Y P * (jb,jy) → Y T Δ(jb,ja),row element: Ring network without grounding branch); Preserve three-phase voltage variables, and delete Buses (X, Y, Z, x, y, z) by needed, which can be obtained a 6 × 6 standard matrix (shown as Fig. 5) to a 7 × 7 matrix by retention of neutral buses.
Based on matrix operation method instead of scanning the branch, the models of the transformer are derived from the relationship of the connections, which directly forms the nodal admittance matrix. The analysis method can be used to three-winding transformers as well.
In accordance with the above rules, we can obtain two incidence matrixes (C Y and C d11 ). The Y-bus model relationship (containing neutral voltage variable) of the transformer in YN,d11 connection between (C Y ,C d11 ) and Y P can be given by The model in the Eq. (5) is the complete full transformer admittance model. Generally, we need to retain the related parameters of the three-phase voltage variables. Repeat Step 3, we can enable to obtain a 6 × 6 matrix.

Modeling of three-phase transformers
In this section, the generalized modeling methodology is applied to represent the threephase constructions of transformers. The aim is to demonstrate the matrices of the models in derivation process. All the three-phase two winding transformers are defined in magnetic circuits of asymmetry configurations in this paper, and symmetric configuration is a special kind of asymmetric configurations.
There are the magnetic circuits connecting closely of the transformer in three-phase three-limb core, besides the magnetic coupling of the primary and secondary windings of each phase, as well as the magnetic circuits coupling of the different phase windings as shown Fig. 6. The effects of the coupling of the inter-phase windings are obvious when the transformer runs asymmetrically.
Considering mutual inductance between the windings, the branch current equation for conveniently analyzing mutual inductance is Bus X and x, Bus Y and y, Bus Z and z are checked in unequal potentials, and those buses do not connect together. The self-impedance of each winding is nearly equal in no-load test, and the relationship of the self-impedance can be given as According to the nodal voltage Eqs. (6) and (7), the nodal admittance matrix can be shown by 12 × 12 dimensions as z AA = z BB = z CC = z aa = z bb = z cc = z s where y ′ m is the mutual admittance between the primary and secondary sides of windings on the same iron-core. ẏ m is the mutual admittance among different phases primary (secondary) sides of windings on the different iron-cores.ÿ m is the mutual admittance among different phases from the primary sides to secondary sides of windings on different iron-cores. y s is the self-admittance on the primary sides or the secondary sides of windings.
The Bus X, Y, Z and Bus x, y, z of transformer can be connected according to the connection of transformer to express the corresponding voltages in neutral points, based on the method of Sect. 2.
(1) YN,yn0 connection The network topology of the transformer in YN,yn0 connection can be shown as Fig. 7. And its 6 × 6 Y-bus matrix can be given by the matrix operation method as Three-phase two winding transformer in magnetic circuit asymmetry configuration; Three-phase two winding transformer in magnetic circuit asymmetry configuration shows that there are coupling in windings, as well as the coupling between phases Ignoring mutual-inductances among the three phases ( ẏ = 0,ÿ = 0, y ′ m = y m ), the 6 × 6 admittance matrix can be given by   Fig. 7 Topology of the three-phase three-limb core transformer; the three-phase three-limb core transformer has asymmetric magnetic circuits. The topology shows the relationship between equivalent reactance and connection of the transformer The network topology of the transformer in the Eq. (10) can be shown as Fig. 8.
(2) YN,d11 connection In Fig. 9, the connected topology of the transformer in YN,d11 connection is presented. For the transformer in the asymmetric magnetic circuits (such as the transformer in threephase three-limb cores), the node admittance matrix is full because of the coupling among the windings.
The boundary condition is given by the Eq. (4). And the Buses X, Y, Z are grounded. We enable to obtain the 6 × 6 Y-bus matrix by the matrix operation method retaining the threephase variables, shown as  (1) Ignoring mutual-inductances among the three phases, the model is able to be seen as the connection by three single-phase transformers assembling. The 6 × 6 admittance matrix can be given as The network topology of the transformer in the Eq. (12) can be shown as Fig. 10.

Impedances obtained in the symmetry
Considering the symmetry of structural parameters for three-phase Transformer, the relationship of the Eq. (6) in per-unit system enables to be expressed by where Z P * is the self-admittance on the primary winding; Z s * is the self-admittance on the secondary winding; Z m * is the mutual-admittance between the primary and the secondary windings at the same iron-core; Z ′ m * is the mutual-admittance between the primary windings. Z ′′ m * is the mutual-admittance between the primary and the secondary windings at the different iron-cores.
Z ′′′ m * is the mutual-admittance between the secondary windings.
(1) No-load test: The test enables to obtain the two groups of data, and they are no-load current I 0 and no-load loss P 0 .The relationships of parameters in no-load state can be given by where y0* is the no-load admittance; z0* is no-load impedance; g0* is the no-load conductance; b0* is the no-load susceptance; STN is the transformer's rated capacity. The no-load impedance can be expressed by (2) Short circuit test: The test enables to obtain the two groups of data, they are the percentage of short-circuit voltage V s % and short circuit power loss P k . The relationships of parameters at short-circuit state can be given by where R T* is the transformer winding resistance; |z T * | is the leakage reactance; x T* is the winding reactance. And the impedance in short-circuit state can be given by The Eqs. (10)-(20) deduced under the condition of the symmetrical currents or voltages, the above parameters stand for the positive sequence parameters. The steps for extrapolating the phase-coordinate parameters can be shown as follows: Considering the same values of leakage reactance on the primary side and secondary side, the relationship (in value system) can be shown by where Z σ 1 ,Z σ 2 are the leakage reactances on the primary and secondary sides as Fig. 2b shown. The basic parameters in the symmetry in no load test can be given by The reluctances of transformer's magnetic circuit are mainly from the air gap between iron cores. According to the principle of magnetic circuit of phase separation test method, there is the relationship shown as follows: The leakage fluxes run in the air different from those run in the iron-cores. The three-phase symmetrical currents have the effect of magnetization. Similar to threephase transmission lines, they have the basic relationships of positive sequence impedance and zero sequence impedance. The relationship between positive sequence leakage reactance and phase separation leakage reactance of transformer obtained from short circuit test is as follows: where x σ s is self-inductance leakage flux; x σ m is mutual-inductance leakage flux.
The self-impedance and mutual impedance of the leakage reactance from the short current test can be given by By substituting the Eqs. (14)-(25), the relationship of the impedance matrix by phasecoordinates in YN,yn0 connection can be given by the modified equation in per-unit system shown as Figure 11 presents the relationships of the vectors in the parametric systems. The relationships for the parameters in References are shown in Fig. 11a. While the relationships for our modified parameters are shown in Fig. 11b, which enable to explain the work in this section intuitively.

Impedances and admittances obtained in the asymmetry
The magnetic circuits of the three-phase three-limb core transformer can be shown in Fig. 12. The branch Eq. (6) enables to express the relationships of the transformer conveniently.
Considering the relations of the linear circuits without magnetic saturation, the impedance parameters in the Eq. (6) can be obtained by the open circuit test. Figure 13 shows the coupling relations between the primary side and the secondary side. Fig. 11 Vector relations in parametric systems; Unmodified parametric system is used to describe vector relations by sequence-coordinates, which only the coupling relationship among windings is considered in (a).While Modified parametric system is used to describe vector relations by phase-coordinates, which coupling relationship among windings is considered containing the tree-phase coupling relations in (b) (1) At wye side for the open circuit test, when the voltage V AX = V N (where V N stand for the rated voltage) is applied at the winding "AX", and the rest windings in opening. We enable to obtain test data for the Eq. (6) shown as cz are the testing voltages; I 0 A0 is the testing current. We can obtain the impedances by per-unit system from the Eq. (20) expressed as z AA z AB z AC z Aa z Ab z Ac z BA z BB z BC z Ba z Bb z Bc z CA z CB z CC z Ca z Cb z Cc z aA z aB z aC z aa z ab z ac z bA z bB z bC z ba z bb z bc z cA z cB z cC z ca z cb z cc 12 Magnetic circuits of the three-phase three-limb core transformer; Magnetic circuits of the three-phase three-limb core transformer is showed when The current I A0 is bring to bear on A-phase winding in primary side Fig. 13 Coupling relations between the primary side and the secondary side by phase-coordinates; Vector graph of the relationship between voltage and coupling is used to describe the Coupling relations between the primary side and the secondary side by phase-coordinates There are the transformers given as (Type I: Capacity: 60 k VA, Ratio: 400 V/110 V; Type II: Capacity: 1000 k VA, Ratio: 6300 V/400 V). In the Table 1, the mutual inductance coefficient is inversely proportional to the length of magnetic circuit. Ignore measurement errors, the mutual inductance (reactance) coefficients can be approximated as: a-phase [1 0.66 0.33]; b-phase [0.5 1 0.5]; c-phase [0.33 0.66 1] (shown as Fig. 12).
(2) Simulation calculation The IEEE 4 node feeder test system network is shown in the Fig. 14. The parameters of the test system are as follows: transformer ratio is 12.47(kV):24.9 (kV), and the parameters of transformer and line present in [13]. The tolerance for calculation is 10 −5 for testing. The unbalanced loads in bus 4 of the test system are 1250 kW, 1800 kW, 2375 kW, and the power factors are − 0.85, − 0.9 and − 0.95 respectively (Complex power of load marked as S 1 ).
There are 6 types and parameters of the transformers for testing by power flow calculations. Table 2 presents the types and the parameters of the transformers. In Table 2, The admittance matrix in No. 1 and No. 2 groups of the transformers employ the method in [8] and [12]. And the rest of the transformers use the method in this paper. Table 3 lists the results of power flow calculation for the transformers in D,yn1 connection. And  Table 4 presents the results of power flow calculation for the transformers in YN,yn0 connection. Table 5 lists the power flow results for the different power supply in the transformer in YN,yn0 connection.

Discussion
The calculated results of No. 1 groups at Tables 3 and 4 are the same as [13].Comparing with the results in No. 1 groups of Tables 3 and 4, The results of No. 2 groups are merely smaller, but less errors considering no-load parameters. While both the admittance   matrix parameters in the two groups of transformers are calculated by the method in [8] and [12], which the parameters essentially are the positive sequence ones. The method employed in Sect. 5 of this paper is used to calculate the three-phase admittance matrix parameters, and the experiments also show the differences in the test. Comparing with   Tables 3 and 4 are a little smaller and closer in the same condition. Analysis from the physical definition, the true phase admittances are different from the positive sequence ones due to the demagnetization by three-phase symmetrical currents. And the results of different magnetic circuits are given in Tables 3 and 4 as well. Those results also are different from each other.
In Table 5, comparing with the No. 1 group, c-phase voltages are increasing gradually by the changes of loads in No. 6 group, due to the increased unbalances of loads' currents. Comparing with symmetric component parameters, the effect of node voltages is more different by the phase-coordinates parameters.
Consequently, the three-phase symmetrical currents have the demagnetization, which make the phase-impedances slightly smaller than the synthetic impedances (positive sequence impedance). In addition, the mutual-impedances effect significantly by asymmetry.

Conclusions
In this paper, the authors propose the phase-coordinate modeling method by a 12 × 12 primitive full admittance matrix to structure the 6 × 6 admittance matrix (7 × 7 matrix by retention of 1 neutral point) for three-phase transformers in different connections from asymmetric magnetic circuits by the matrix operation method. The phase-coordinate impedances have been corrected by considering the magnetization of three-phase symmetrical currents. There are the advantages shown as follows: • Modeling the transformers in different connections and asymmetric magnetic circuits; • Real phase-coordinates parameters, which enable calculations accurately.
Finally, the matrix operation method also enables to use for modeling three-phase three-winding transformers usefully, besides the three-phase double-winding transformers. And the method applies to model the non-standard transformers as well.