Implementation of an Alternative Frequency-Dependent Three-Phase Transmission Line Model Based on the Folded Line Equivalent Model in MatLab-Simulink

Abstract

This paper proposes an alternative multiconductor transmission line model that combines the folded line equivalent with the modal transformation. The folded line equivalent decomposes the nodal admittance matrix of a transmission line into its open-circuit and short-circuit contributions. These contributions are fitted to rational functions, which are associated with Norton equivalent circuits based on their state space models. The proposed model uses an orthogonal matrix to transform voltages and currents from the phase domain to the folded line equivalent domain and vice versa. Because the transformation matrix is orthogonal, we represent it using ideal transformers in simulation software. First, we use a circuit representation of Clarke’s matrix to decompose a transmission line into its modes. Then, each mode is decomposed into its open-circuit and short-circuit contributions using a circuit implementation of the proposed matrix. The proposed approach can accurately represent short lines in simulations with time steps equal to or greater than the propagation time of the transmission line. We compare the results obtained with the proposed approach to those obtained with power systems computer-aided design/electromagnetic transients including the DC universal line model.

1. Introduction and Literature Review

Frequency-dependent cable and transmission line (TL) models have been widely used in the study of electromagnetic transients in power systems. These models can be developed in the modal domain or directly in the phase domain. For instance, the JMarti model, introduced in [1], assumes a constant matrix to transform electrical parameters to the mode domain and vice versa. This assumption is valid for ideally transposed TLs and gives good results for non-transposed TLs with vertical symmetry. For this reason, the JMarti model is available as a built-in model in simulation software such as the Alternative Transient Program-Electromagnetic Transients Program (ATP-EMTP) [2]. However, assuming frequency-independent transformation matrices is not valid for non-transposed asymmetric lines or parallel circuits. Other TL models that do not make this assumption are preferred for these cases e.g., the universal line model (ULM) [3]. The ULM computes voltages and currents directly in the phase domain and is considered one of the most accurate and general models for the simulation of overhead transmission lines and cables. The ULM is available as a built-in model in simulation software such as the Power Systems Computer Aided Design (PSCAD) [4] and the Electromagnetic Transients Program—Restructured Version (EMTP-RV) [5].

Numerous transmission line models (including the ULM and the JMarti model) derive from the method of characteristics. The method of characteristics requires (simulation) time steps to be smaller than the propagation time of the transmission line. Therefore, the time step used in the simulation of complex power systems is limited to the propagation time of the shortest line [6,7]. The folded line equivalent (FLE) [8] circumvents this problem. The FLE decouples a TL into two blocks that represent the short-circuit and open-circuit contributions of the TL. These blocks are fitted to rational functions using the vector fitting (VF) algorithm [9,10]. The state-space representations of the open-circuit and short-circuit contributions are then associated with equivalent Norton circuits, which are compatible with simulation software such as the EMTP. The main advantage of the FLE is that it provides accurate results for short lines when the time step is greater than the propagation time of the transmission line, reducing execution times in the simulation of complex power systems. In contrast, other transmission line models require smaller time steps and thereby take longer to compute simulations.

This paper proposes an EMTP-compatible multi-conductor TL model based on the FLE. The proposed model first transforms voltages and currents to the mode domain using a circuit representation of Clarke’s matrix. Then, voltages and currents are transformed to the FLE domain using a new orthogonal transformation matrix, different from the one originally proposed in [8]. The proposed transformation matrix has a simpler interpretation and allows a more intuitive implementation using ideal transformers, similar to the arrangement used to represent Clarke’s matrix. We implement the proposed model in the MatLab/Simulink simulation software [11].

The proposed approach inherits the advantages of the FLE and thus can be simulated with time steps greater than the propagation time of the TL. Because our approach is combined with modal transformation, it provides a better understanding of mode behavior during high-frequency phenomena such as lightning strikes or switching operations, which could cause outages and damage equipment. Moreover, analyzing three independent modes is simpler and computationally more efficient than analyzing a coupled three-phase TL. Because the proposed model uses basic circuit elements, it can be implemented in simulation software e.g., EMTP-RV, PSCAD/EMTDC, and ATP. Results show that the proposed approach outputs accurate results when the step size of the simulations is greater than the propagation time of the TL.

2. The Folded Line Equivalent

The nodal admittance matrix of the TL of Figure 1 relates voltages ${\mathit{V}}_1=$ ${\left[\begin{array}{ccc}{V}_{1,1}& \dots & V_{\mathrm{n},1}\end{array}\right]}^T$ and currents ${\mathit{I}}_1={\left[\begin{array}{ccc}{I}_{1,1}& \dots & I_{2,1}\end{array}\right]}^T$ at one of the line terminals to the voltages ${\mathit{V}}_2={\left[\begin{array}{ccc}{V}_{1,2}& \dots & V_{\mathrm{n},2}\end{array}\right]}^T$ and currents ${\mathit{I}}_2=$ ${\left[\begin{array}{ccc}{I}_{1,2}& \dots & I_{2,2}\end{array}\right]}^T$ at the other terminal, as follows:

(1)$\left[\begin{array}{c}{\mathit{I}}_1\\ {\mathit{I}}_2\end{array}\right]=\left[\begin{array}{cc}{\mathit{Y}}_{\mathrm{s}}& {\mathit{Y}}_{\mathrm{m}}\\ {\mathit{Y}}_{\mathrm{m}}& {\mathit{Y}}_{\mathrm{s}}\end{array}\right]\left[\begin{array}{c}{\mathit{V}}_1\\ {\mathit{V}}_2\end{array}\right].$

Refs. [12,13,14] compute the frequency-dependent nodal admittance matrix in (1) from the per unit length (p.u.l.) TL parameters. The FLE takes advantage of the symmetry of the nodal admittance matrix in (1) and decouples it into its open-circuit and short-circuit contributions as follows:

(2)$\left[\begin{array}{c}{\mathit{I}}_{\mathrm{oc}}\\ {\mathit{I}}_{\mathrm{sc}}\end{array}\right]=\left[\begin{array}{cc}{\mathit{Y}}_{\mathrm{oc}}& \mathbf{0}\\ \mathbf{0}& {\mathit{Y}}_{\mathrm{sc}}\end{array}\right]\left[\begin{array}{c}{\mathit{V}}_{\mathrm{oc}}\\ {\mathit{V}}_{\mathrm{sc}}\end{array}\right]$

where the open-circuit admittance ${\mathit{Y}}_{\mathrm{oc}}$ and the short-circuit admittance ${\mathit{Y}}_{\mathrm{sc}}$ are, respectively, given by

(3)${\mathit{Y}}_{\mathrm{oc}}={\mathit{Y}}_{\mathrm{s}}+{\mathit{Y}}_{\mathrm{m}}$

(4)${\mathit{Y}}_{\mathrm{sc}}={\mathit{Y}}_{\mathrm{s}}-{\mathit{Y}}_{\mathrm{m}}.$

Voltages and currents are transformed to and from the FLE domain as follows:

(5)$\left[\begin{array}{c}{\mathit{V}}_{\mathrm{oc}}\\ {\mathit{V}}_{\mathrm{sc}}\end{array}\right]={\mathit{K}}^{-1}\left[\begin{array}{c}{\mathit{V}}_1\\ {\mathit{V}}_2\end{array}\right]$

(6)$\left[\begin{array}{c}{\mathit{I}}_{\mathrm{oc}}\\ {\mathit{I}}_{\mathrm{sc}}\end{array}\right]={\mathit{K}}^{-1}\left[\begin{array}{c}{\mathit{I}}_1\\ {\mathit{I}}_2\end{array}\right]$

where the transformation matrix ${\mathit{K}}^{-1}$ is given by

(7)$\mathit{K}=\left[\begin{array}{cc}\mathbf{U}& \mathbf{U}\\ \mathbf{U}& -\mathbf{U}\end{array}\right],$

and $\mathbf{U}$ is the identity matrix.

3. Our Approach: Proposed Folded Line Equivalent

The transformation of voltages $\mathit{V}$ and currents $\mathit{I}$ at line terminals, longitudinal impedances $\mathit{Z}$, and shunt admittances $\mathit{Y}$ of the TL of Figure 1 from the phase domain to the mode domain and vice versa is given by

(8)$\mathit{V}={\mathit{T}}_{\mathrm{V}}{\mathit{V}}_{\alpha \beta 0}$

(9)$\mathit{I}={\mathit{T}}_{\mathrm{I}}{\mathit{I}}_{\alpha \beta 0}$

(10)${\mathit{Z}}_{\alpha \beta 0}={\mathit{T}}_{\mathrm{I}}^{\mathrm{T}}\mathit{Z}{\mathit{T}}_{\mathrm{I}}$

(11)${\mathit{Y}}_{\alpha \beta 0}={\mathit{T}}_{\mathrm{I}}^{-1}\mathit{Y}{\mathit{T}}_{\mathrm{I}}^{-\mathrm{T}}.$

In (8) and (9), the transformation matrices ${\mathit{T}}_{\mathrm{I}}$ and ${\mathit{T}}_{\mathrm{V}}$ are frequency-dependent.

Clarke’s matrix (12) is a frequency-constant orthogonal transformation matrix that gives accurate results for perfectly transposed TLs and relatively accurate results for TLs with vertical symmetry [15].

(12)${\mathit{T}}_{\mathrm{clk}}=\left[\begin{array}{ccc}2/\sqrt{6}& -1/\sqrt{6}& -1/\sqrt{6}\\ 0& 1/\sqrt{2}& -1/\sqrt{2}\\ 1/\sqrt{3}& 1/\sqrt{3}& 1/\sqrt{3}\end{array}\right]$

We propose to use Clarke’s matrix to decouple three-phase TLs into its modes, as shown in Figure 2. Then, we transform each independent mode to the FLE domain, as shown in Figure 3. We represent each of these transformations using circuit elements, as shown in Figure 4.

3.1. Modal Transformation: Circuit Representation of Clarke’s Matrix

The block labeled as transformation matrix in Figure 2 transforms voltages and currents from the phase domain to the mode domain and vice versa. This block contains the arrangement of ideal transformers of Figure 4a [16], which relates voltages at the left and right terminals Figure 4a as follows:

(13)$v_{\mathrm{a}}=v_{\mathrm{a}1}+v_{\mathrm{a}3}=\frac{2}{\sqrt{6}}{v}_{\alpha }+\frac{1}{\sqrt{3}}{v}_0$

(14)$v_{\mathrm{b}}=v_{\mathrm{b}1}+v_{\mathrm{b}2}+v_{\mathrm{b}3}=-\frac{1}{\sqrt{6}}{v}_{\alpha }+\frac{1}{\sqrt{2}}{v}_{\beta }+\frac{1}{\sqrt{3}}{v}_0$

(15)$v_{\mathrm{c}}=v_{\mathrm{c}1}+v_{\mathrm{c}2}+v_{\mathrm{c}3}=-\frac{1}{\sqrt{6}}{v}_{\alpha }+\frac{1}{\sqrt{2}}{v}_{\beta }+\frac{1}{\sqrt{3}}{v}_0$

The relationship between currents at right and left terminals of Figure 4a is given by

(16)$i_{\alpha }=i_{\alpha 1}+i_{\alpha 2}+i_{\alpha 3}=\frac{2}{\sqrt{6}}{i}_{\mathrm{a}}-\frac{1}{\sqrt{6}}{i}_{\mathrm{b}}-\frac{1}{\sqrt{6}}{i}_{\mathrm{c}}$

(17)$i_{\beta }=i_{\beta 2}+i_{\beta 3}=\frac{1}{\sqrt{2}}{i}_{\mathrm{b}}-\frac{1}{\sqrt{2}}{i}_{\mathrm{c}}$

(18)$i_0=i_{01}+i_{02}+i_{03}=\frac{1}{\sqrt{3}}{i}_{\mathrm{a}}+\frac{1}{\sqrt{3}}{i}_{\mathrm{b}}+\frac{1}{\sqrt{3}}{i}_{\mathrm{c}}$

In matrix form, the relationships presented in (13)–(18) can be written as follows:

(19)${\mathit{V}}_{\mathrm{abc}}={\mathit{T}}_{\mathrm{clk}}^{-1}{\mathit{V}}_{\alpha \beta 0}$

(20)${\mathit{I}}_{\alpha \beta 0}={\mathit{T}}_{\mathrm{clk}}{\mathit{V}}_{abc}$

which is equivalent to (9) when ${\mathit{T}}_{\mathrm{I}}={\mathit{T}}_{\mathrm{clk}}$ (Clarke’s transformation).

3.2. FLE Transformation: Circuit Representation of the Proposed Transformation Matrix

We combine the circuit representation of the $\alpha \beta 0$ decomposition of Section 3.1 with the circuit representation of the FLE transformation presented in this section. After the three-phase TL is decomposed into its modes, as shown in Figure 2, each mode is individually transformed to the FLE domain, as shown in Figure 3. Each mode has its open-circuit (OC) and short-circuit (SC) contributions. The block labeled as Matrix K in Figure 3 transforms voltages and currents from the mode domain to the FLE domain and vice versa. This block contains the arrangement of ideal transformers of Figure 4b, which relates voltages at terminals 1 and 2 to voltages at terminals OC and SC as follows:

(21)$V_{1,\mathrm{mk}}=V_{1,\mathrm{mk}1}+V_{1,\mathrm{mk}2}=\frac{1}{\sqrt{2}}{V}_{\mathrm{oc},\mathrm{mk}}+\frac{1}{\sqrt{2}}{V}_{\mathrm{sc},\mathrm{mk}}$

(22)$V_{2,\mathrm{mk}}=V_{2,\mathrm{mk}1}+V_{2,\mathrm{mk}2}=\frac{1}{\sqrt{2}}{V}_{\mathrm{oc},\mathrm{mk}}-\frac{1}{\sqrt{2}}{V}_{\mathrm{sc},\mathrm{mk}}$

where subscript mk = $\left\{\alpha ,\beta ,0\right\}$. The relationship between currents at terminals 1 and 2 to voltages at terminals OC and SC is given by

(23)$I_{\mathrm{oc},\mathrm{mk}}=I_{\mathrm{oc},\mathrm{mk}1}+I_{\mathrm{oc},\mathrm{mk}2}=\frac{1}{\sqrt{2}}{I}_{1,\mathrm{mk}}+\frac{1}{\sqrt{2}}{I}_{2,\mathrm{mk}}$

(24)$I_{\mathrm{sc},\mathrm{mk}}=I_{\mathrm{sc},\mathrm{mk}1}+I_{\mathrm{sc},\mathrm{mk}2}=\frac{1}{\sqrt{2}}{I}_{1,\mathrm{mk}}-\frac{1}{\sqrt{2}}{I}_{2,\mathrm{mk}}$

In matrix form, the relationships presented in (21)–(24) can be written as follows:

(25)$\left[\begin{array}{c}{V}_{1,\mathrm{mk}}\\ V_{2,\mathrm{mk}}\end{array}\right]=\mathit{K}\left[\begin{array}{c}{V}_{\mathrm{oc},\mathrm{mk}}\\ V_{\mathrm{sc},\mathrm{mk}}\end{array}\right]$

(26)$\left[\begin{array}{c}{I}_{\mathrm{oc},\mathrm{mk}}\\ I_{\mathrm{sc},\mathrm{mk}}\end{array}\right]={\mathit{K}}^{-1}\left[\begin{array}{c}{I}_{1,\mathrm{mk}}\\ I_{2,\mathrm{mk}}\end{array}\right].$

In (25) and (26), instead of using the transformation Matrix K introduced in [8], we propose using the following transformation matrix to decouple single-phase TL into its open-circuit and short-circuit contributions.

(27)$\mathit{K}=\left[\begin{array}{cc}\frac{1}{\sqrt{2}}& \frac{1}{\sqrt{2}}\\ \frac{1}{\sqrt{2}}& -\frac{1}{\sqrt{2}}\end{array}\right]$

The proposed transformation matrix preserves the properties and relationships of the transformation matrix introduced in [8], yet it allows a simpler and more intuitive implementation. Figure 4b shows the circuit implementation of (27) using an arrangement of ideal transformers.

3.3. Computation of the Nodal Admittance Matrix of Each Mode

Each of the modes resulting from the transformation presented in Section 3.1 is considered an independent transmission line and thus has a nodal admittance matrix given by

(28)${\mathit{Y}}_{\mathrm{mk}}=\left[\begin{array}{cc}{Y}_{\mathrm{mk},\mathrm{s}}& Y_{\mathrm{mk},\mathrm{m}}\\ Y_{\mathrm{mk},\mathrm{m}}& Y_{\mathrm{mk},\mathrm{s}}\end{array}\right].$

The nodal admittance matrix of the mk-th mode can be calculated from the propagation function $H_{\mathrm{mk}}$ and characteristic admittance $Y_{\mathrm{c},\mathrm{mk}}$ of the mode as follows

(29)$Y_{\mathrm{mk},\mathrm{s}}={\left(1-H_{\mathrm{mk}}^2\right)}^{-1}\left(H_{\mathrm{mk}}^2+1\right)Y_{\mathrm{c},\mathrm{mk}}$

(30)$Y_{\mathrm{mk},\mathrm{m}}=-2{\left(1-H_{\mathrm{mk}}^2\right)}^{-1}{H}_{\mathrm{mk}}{Y}_{\mathrm{c},\mathrm{mk}}$

where $H_{\mathrm{mk}}$ and $Y_{\mathrm{c},\mathrm{mk}}$ are computed from the mode’s p.u.l. impedance $Z_{\mathrm{mk}}$, p.u.l. admittance $Y_{\mathrm{mk}}$, and the length of the line ℓ, as shown below

(31)$Y_{\mathrm{c},\mathrm{k}}=Z_{\mathrm{mk}}^{-1}\sqrt{Z_{\mathrm{mk}}{Y}_{\mathrm{mk}}}$

(32)$H_{\mathrm{k}}={\mathrm{e}}^{-\ell \sqrt{Y_{\mathrm{mk}}{Z}_{\mathrm{mk}}}}.$

There are alternative procedures of computing $Y_{\mathrm{mk},\mathrm{s}}$ and $Y_{\mathrm{mk},\mathrm{m}}$. In [12], Morched et al. compute the elements of the modal admittance matrix of a single-phase TL as follows:

(33)$Y_{\mathrm{mk},\mathrm{s}}=Y_{\mathrm{mk}}{\Lambda }_{\mathrm{mk},\mathrm{zy}}^{-1}coth\left({\Lambda }_{\mathrm{mk},\mathrm{zy}}\ell \right)$

(34)$Y_{\mathrm{k},\mathrm{m}}=-Y_{\mathrm{mk}}{\Lambda }_{\mathrm{mk},\mathrm{zy}}^{-1}\csc \mathrm{h}\left({\Lambda }_{\mathrm{mk},\mathrm{zy}}\ell \right)$

where propagation constant ${\Lambda }_{\mathrm{mk},\mathrm{zy}}=\sqrt{Z_{\mathrm{mk}}{Y}_{\mathrm{mk}}}$. Alternatively, in [14], Gustavsen avoids hyperbolic functions and computes $Y_{\mathrm{mk},\mathrm{s}}$ and $Y_{\mathrm{mk},\mathrm{m}}$ as follows:

(35)$Y_{\mathrm{mk},\mathrm{s}}=Z_{\mathrm{mk}}^{-1}{a}_{\mathrm{mk}}$

(36)$Y_{\mathrm{mk},\mathrm{m}}=Z_{\mathrm{mk}}^{-1}{b}_{\mathrm{mk}}$

where

(37)$a_{\mathrm{mk}}=\frac{\sqrt{d_{\mathrm{mk}}}\left(1+h_{\mathrm{mk}}^2\right)}{1-h_{\mathrm{mk}}^2}$

(38)$b_{\mathrm{mk}}=\frac{-2\sqrt{d_{\mathrm{mk}}}.h_{\mathrm{mk}}}{1-h_{\mathrm{mk}}^2}$

(39)$h_{\mathrm{mk}}=\exp \left(-\sqrt{d_{\mathrm{mk}}}\ell \right)$

(40)$d_{\mathrm{mk}}={\left(\sqrt{Z_{\mathrm{mk}}{Y}_{\mathrm{mk}}}\right)}^{-1}$

3.4. Circuit Representation of the FLE Parameters of Each Mode

The FLE parameters (open-circuit and short-circuit contributions) of each mode can be represented by Norton equivalent circuits [17] or an arrangement of resistors, capacitors, and inductors [18]. In this paper, we employ Norton equivalent circuits to represent the open-circuit and short-circuit contributions for each mode in simulation software.

The VF algorithm [9,10] fits the open-circuit $Y_{\mathrm{oc},\mathrm{mk}}$ and short-circuit $Y_{\mathrm{sc},\mathrm{mk}}$ contributions of each mode to rational functions. We use both the VF algorithm and the fast modal perturbation (FMP) [19] to ensure passivity in the fitted rational functions. The state-space model of each fitted matrix is then generated following the procedure described in [20]. The state-space model that describes the open-circuit contribution of the mk-th mode is given by

(41)${\dot{\mathit{x}}}_{\mathrm{mk},\mathrm{oc}}\left(t\right)={\mathit{A}}_{\mathrm{mk},\mathrm{oc}}{\mathit{x}}_{\mathrm{mk},\mathrm{oc}}\left(t\right)+{\mathit{B}}_{\mathrm{mk},\mathrm{oc}}{\mathit{v}}_{\mathrm{mk},\mathrm{oc}}\left(t\right)$

(42)$i_{\mathrm{mk},\mathrm{oc}}\left(t\right)={\mathit{C}}_{\mathrm{mk},\mathrm{oc}}{\mathit{x}}_{\mathrm{mk},\mathrm{oc}}\left(t\right)+D_{\mathrm{mk},\mathrm{oc}}{v}_{\mathrm{mk},\mathrm{oc}}\left(t\right)$

Similarly, the state-space model that describes the short-circuit contribution of the mk-th mode is given by

(43)${\dot{\mathit{x}}}_{\mathrm{mk},\mathrm{sc}}\left(t\right)={\mathit{A}}_{\mathrm{mk},\mathrm{sc}}{\mathit{x}}_{\mathrm{mk},\mathrm{sc}}\left(t\right)+{\mathit{B}}_{\mathrm{mk},\mathrm{sc}}{\mathit{v}}_{\mathrm{mk},\mathrm{sc}}\left(t\right)$

(44)$i_{\mathrm{mk},\mathrm{sc}}\left(t\right)={\mathit{C}}_{\mathrm{mk},\mathrm{sc}}{\mathit{x}}_{\mathrm{mk},\mathrm{sc}}\left(t\right)+D_{\mathrm{mk},\mathrm{sc}}{v}_{\mathrm{mk},\mathrm{sc}}\left(t\right)$

In (41)–(44), state matrices $\mathit{A}$, $\mathit{B}$, $\mathit{C}$, and D are constant.

The Norton equivalent circuit of a state-space model consists of a constant conductance G and a historic current $i_{\mathrm{hist}}$ that updates its value in each time step, as portrayed in Figure 3. The historic current sources and the constant parameters of the Norton equivalent circuit allow the proposed model to consider the wave phenomena. Application of the trapezoidal rule of integration to (41)–(44) results in the following Norton parameters [17]:

(45)$G_{\mathrm{ck},\mathrm{mk}}=D_{\mathrm{ck},\mathrm{mk}}+{\mathit{C}}_{\mathrm{ck},\mathrm{mk}}{\lambda }_{\mathrm{ck},\mathrm{mk}}{\mathit{B}}_{\mathrm{ck},\mathrm{mk}}$

(46)$i_{\mathrm{hist},\mathrm{ck},\mathrm{mk}}\left[n\right]=-{\mathit{C}}_{\mathrm{ck},\mathrm{mk}}\left({\mathit{\alpha }}_{\mathrm{ck},\mathrm{mk}}+1\right){\mathit{\lambda }}_{\mathrm{ck},\mathrm{mk}}{\mathit{x}}_{\mathrm{ck},\mathrm{mk}}\left[n\right]$

where

(47)${\mathit{\lambda }}_{\mathrm{ck},\mathrm{mk}}={\left(1-{\mathit{A}}_{\mathrm{ck},\mathrm{mk}}\frac{\Delta t}{2}\right)}^{-1}\frac{\Delta t}{2}$

(48)${\mathit{\alpha }}_{\mathrm{ck},\mathrm{mk}}={\left(1-{\mathit{A}}_{\mathrm{ck},\mathrm{mk}}\frac{\Delta t}{2}\right)}^{-1}\left(1+{\mathit{A}}_{\mathrm{ck},\mathrm{mk}}\frac{\Delta t}{2}\right)$

(49)${\mathit{x}}_{\mathrm{ck},\mathrm{mk}}\left[n\right]={\mathit{\alpha }}_{\mathrm{ck},\mathrm{mk}}{\mathit{x}}_{\mathrm{ck},\mathrm{mk}}[n-1]+{\mathit{B}}_{\mathrm{ck},\mathrm{mk}}{v}_{\mathrm{ck},\mathrm{mk}}[n-1]$

and ck = $\left\{\mathrm{oc},\mathrm{sc}\right\}$, mk = $\left\{\alpha ,\beta ,0\right\}$. The Norton parameters in (45) and (46) are EMTP-compatible because all passive components are represented by Norton equivalent circuits [6].

4. Results

4.1. Line Configuration and Simulation Approach

To show results, we

• Computed the p.u.l. line parameters of the TL of Figure 5 using the PSCAD/EMTDC’s Line Constants Program;

• Simulated the TL using the proposed approach, implemented in Matlab/Simulink, as shown in Figure 4;

• Simulated the TL using PSCAD/EMTDC’s ULM;

• Compared the outputs of both models.

Figure 5

Three-phase transmission line geometry.

[Figure omitted. See PDF]

The flowchart of Figure 6 shows all the steps mentioned above, needed to implement the proposed model in MatLab’s Simulink (ver. 2021b).

One of the benefits of the proposed approach is that it provides accurate results when the time step is greater than the propagation delay of the TL. This is particularly useful in networks containing short and long lines because the time step of the entire network depends on the propagation delay of the shortest line. To show this advantage, we considered the TL of Figure 5 300 m long, resulting in a propagation delay of approximately 1 $\mathsf{\mu }$s. We show that the PFLE outputs accurate results in simulations with time steps greater than 1 $\mathsf{\mu }$s.

4.2. Fitting Procedure and Passivity Enforcement

The TL of Figure 5 has three modes of propagation: $\alpha$, $\beta$, and 0. Each mode is composed of an open-circuit and short-circuit nodal admittance.

We employed the VF routine [21] to fit the open-circuit and short-circuit admittances of each mode using twenty poles and five iterations. Passivity is independently enforced for each mode using the “RPdriver” routine [20]. Note that the “RPdriver” routine is applied to a frequency range that is twice the fitting frequency. Therefore, it is necessary to compute parameters up to 2 MHz.

Figure 7, Figure 8 and Figure 9 show the magnitude of the open-circuit and short-circuit nodal admittance of each mode (solid blue line), the fitted rational function after enforcing passivity (dotted red line), and the deviation between them (solid green line). Note that the deviation curve has a behavior similar to the data because the weight used in the fitting is inversely proportional to the magnitude of the data. This weighting scheme allows the accurate fitting of small and large values.

The number of poles required to fit a given function usually depends on the number of peaks present in the curve. The number of peaks in the frequency range of power systems transients increases with the length of the line. As a result, short lines require fewer poles to achieve a certain level of accuracy. This paper shows that the proposed model outputs accurate results in the simulation of short lines even when the time step of the simulation is larger than the propagation delay of the TL.

A Norton equivalent circuit is considered passive when the eigenvalues of the real part of the open-circuit nodal admittance $Y_{\mathrm{oc},\mathrm{mk}}$ and short-circuit nodal admittance $Y_{\mathrm{sc},\mathrm{mk}}$ are positive for all frequencies [8,22].

Figure 10, Figure 11 and Figure 12 show that the eigenvalues of the open-circuit and short-circuit conductances of each mode are positive after passivity is enforced.

4.3. Open-Circuit, Capacitive Load, and Single-Phase Fault

4.3.1. Simulation Layout

Figure 13 shows the TL of Figure 5 when one of its terminals is energized at 0 s, and the other terminal is

• 1.. Left open-circuited when the line is energized;

• 2.. Connected to a capacitive load at 4 ms;

• 3.. A single-phase fault occurs at 8 ms.

We compared the results obtained by the proposed approach (labeled as PFLE) to those obtained by PSCAD’s ULM. To show the advantages of the proposed approach, the simulations were performed considering different time steps.

4.3.2. Results

Results show that the proposed approach outputs accurate results even when the time step is four times the propagation delay of the TL. In contrast, the maximum integration time step that the ULM supports is equal to the propagation delay of the TL, which is approximately 1 $\mathsf{\mu }$s. This requirement must be met for all line models based on the method of characteristics. We use a time step equal to 1/10 of the propagation delay of the TL to compute results using the ULM.

In all simulations, the red, blue, and green curves represent phases A, B, and C, respectively. The results are organized as follows:

• 1.. Figure 14 and Figure 15 show, respectively, the voltages at the receiving terminal and the currents at the sending terminal during line energization;

• 2.. Figure 16 and Figure 17 show, respectively, the voltages at the receiving terminal and the currents at the sending terminal when the capacitive load is connected;

• 3.. Figure 18 and Figure 19 show, respectively, the voltages at the receiving terminal and the currents at the sending terminal when a single-phase fault occurs at 8 ms.

Figure 14, Figure 15, Figure 16, Figure 17, Figure 18 and Figure 19 show that the PFLE outputs accurate results with time steps greater than the propagation delay of the TL. This can greatly reduce execution times in the simulation of large complex networks where the time step is limited to the propagation delay of the shortest line when TL models based on the methods of characteristics are used.

4.4. Line Energization

4.4.1. Simulation Layout

Figure 20 shows the TL of Figure 5 when one of its terminals is energized at 0 s and the other terminal is connected to a resistive-inductive load.

4.4.2. Results

Figure 21 and Figure 22 compare the voltages at the receiving terminal simulated with the proposed approach (labeled as PFLE) and PSCAD’s ULM for several time steps. Figure 21 shows the transient voltages right after the line is energized, and Figure 22 shows the voltages after the first high-frequency voltages disappear. Note that using a time step greater than the propagation delay of the TL causes a small delay during the first high-frequency oscillations, which disappears as the TL reaches steady-state conditions, as shown in Figure 22.

Figure 23 shows the currents at the sending terminal of the TL. Note that the proposed approach agrees with the ULM for all time steps.

The execution time of a simulation depends on the time required to solve a single time step and the total number of time steps in the simulation. Reducing the average time to solve one time step is possible through optimization techniques and lightweight models. However, the most effective way of reducing the execution time of simulation is by reducing the total number of time steps that need to be computed. TL models based on the method of characteristics are limited to time steps equal to the propagation delay of the shortest TL in the network. The proposed approach supports time steps greater than the propagation delay of the shortest TL without compromising accuracy. Our results show that the proposed approach outputs accurate results even with time steps that are four times the propagation delay of the shortest TL, resulting in execution times that are up to 400% faster.

4.5. Error

Section 4.3.2 and Section 4.4.2 show the results obtained by implementing the proposed approach, the PFLE, in MatLab’s Simulink. We measure the error of the PFLE using the normalized root-mean-square deviation (NRMSD), which is given by

(50)$\mathrm{NRMSD}=\frac{\sqrt{\frac{1}{N}{\sum }_{i=1}^N{\left(y_{PM,i}-y_{\mathrm{ULM}}\right)}^2}}{max\left\{y_{\mathrm{ULM}}\right\}-min\left\{y_{\mathrm{ULM}}\right\}}$

where N is the number of points of y.

In all simulations, we compare the results obtained by the PFLE to the results obtained by the ULM when the time step is 1/10 of the propagation delay of the TL (1 $\mathsf{\mu }$s). We use PSCAD/EMTDC’s ULM as a reference model because this model is widely accepted in technical literature and professional software due to its accuracy and numerical stability. The ULM computes results directly in the phase domain and time domain while considering the frequency dependence and distributed nature of line parameters.

Our numerical comparisons are organized as follows:

• Table 1 shows the NRMSD of the proposed approach (PFLE) in Figure 14 and Figure 15.

• Table 2 shows the NRMSD of the proposed approach (PFLE) in Figure 16 and Figure 17.

• Table 3 shows the NRMSDof the proposed approach (PFLE) in Figure 18 and Figure 19.

• Table 4 shows the NRMSD of the proposed approach (PFLE) in Figure 22 and Figure 23.

5. Conclusions

In this paper, we propose an alternative frequency-dependent transmission line (TL) model based on the folded line equivalent (FLE) line model to represent three-phase TLs. The proposed approach uses a circuit representation of Clarke’s matrix to transform parameters from the phase domain to the mode domain. Then, instead of using the transformation matrix proposed in [8], we use a new transformation matrix to transform parameters from the mode domain to the FLE domain. The proposed matrix allows the proposed model to be implemented in simulation software using ideal transformers. The open-circuit and short-circuit contributions of the proposed approach are represented in simulation software using Norton equivalent circuits. We implement the proposed approach in MatLab’s Simulink multidomain simulation software.

Traditional TL models based on the method of characteristics can only run in simulations where the time step is equal to or less than the propagation delay of the shortest TL in the network. The main advantage of the proposed approach is that it can run on simulations where the step size is greater than the propagation time of the TL without compromising accuracy. Figure 14 through Figure 23 show good agreement between the results obtained with the ULM and the PFLE. Furthermore, they show that the PFLE provides accurate results even when the time step is four times greater than the propagation delay of the TL. Table 1, Table 2, Table 3 and Table 4 show that the NRMSD of the results obtained by the PFLE, when compared to the results obtained by the ULM with a time step of 1/10 of the propagation delay of the TL, is less than 0.005. These results prove that the proposed model outputs accurate results in simulations where the time step is 10%, 100%, 200%, and 400% of the propagation time of the TL.

Author Contributions

Conceptualization, J.S.L.C., S.K., L.T., P.T.C. and J.P.F.; methodology, J.S.L.C. and L.T.; software, J.S.L.C. and L.T.; validation, J.S.L.C. and L.T.; formal analysis, J.S.L.C., L.T. and P.T.C.; investigation, J.S.L.C. and L.T.; writing—original draft preparation, J.S.L.C., L.T. and P.T.C.; visualization, J.S.L.C. and L.T.; supervision, P.T.C., S.K. and J.P.F.; project administration, S.K. and J.P.F.; funding, S.K. All authors have read and agreed to the published version of the manuscript.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

Abbreviations

The following abbreviations are used in this manuscript:

TL	Transmission line
ATP-EMTP	Alternative Transient Program-Electromagnetic Transients Program
ULM	Universal Line Model
PSCAD	Power Systems Computer Aided Design
EMRP-RV	Electromagnetic Transients Program-Restructured Version
FLE	Folded Line Equivalent
VF	Vector Fitting
p.u.l.	per unit length
FMP	Fast Modal Perturbation
PFLE	Proposed Folded Line Equivalent

Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Figures and Tables

Figure 1. Multiconductor transmission line.

Figure 2. Representation of a three-phase transmission line in the mode domain: abc [Forumla omitted. See PDF.] abc.

Figure 3. Mode representation for the proposed FLE (PFLE) model.

Figure 4. Circuit representation of (a) Clarke’s matrix ([Forumla omitted. See PDF.] transformation) and (b) the proposed FLE transformation matrix (FLE transformation).

Figure 6. Steps required to implement the PFLE in MatLab’s Simulink.

Figure 7. Fitting of the open-circuit admittance [Forumla omitted. See PDF.] and short-circuit admittance [Forumla omitted. See PDF.] of mode [Forumla omitted. See PDF.].

Figure 8. Fitting of the open-circuit admittance [Forumla omitted. See PDF.] and short-circuit admittance [Forumla omitted. See PDF.] of mode [Forumla omitted. See PDF.].

Figure 9. Fitting of the open-circuit admittance [Forumla omitted. See PDF.] and short-circuit admittance [Forumla omitted. See PDF.] of mode 0.

Figure 10. Eigenvalues of the nodal admittance matrix of mode [Forumla omitted. See PDF.].

Figure 11. Eigenvalues of the nodal admittance matrix of mode [Forumla omitted. See PDF.].

Figure 12. Eigenvalues of the nodal admittance matrix of mode 0.

Figure 13. Simulation configuration.

Figure 14. Voltage at the receiving end of the line when the line is energized.

Figure 15. Current at the sending end of the line when the line is energized.

Figure 16. Voltage at the receiving end of the line when the cable is connected.

Figure 17. Current at the sending end of the line when the cable is connected.

Figure 18. Voltage at the receiving end of the line after a single-phase fault occurs.

Figure 19. Current at the sending end of the line after a single-phase fault occurs.

Figure 20. Transmission line connected to a resistive-inductive load.

Figure 21. Voltages at the receiving terminal of the line right after it is energized.

Figure 22. Voltages at the receiving end of the line in the first 20 ms of the simulation.

Figure 23. Current at the sending terminal of the line.

Table 1

NRMSD of the voltage and current curves obtained with the PFLE with respect to the voltage and current curves obtained with the ULM and a time step of 1 $\mathsf{\mu }$s: Figure 14 and Figure 15.

Voltage			Current
Phases	Model: $\Delta$ t	NRMSD	Phases	Model: $\Delta$ t	NRMSD
1(Figure 14)	PFLE: 0.1 $\mathsf{\mu }$s	0.004	1 (Figure 15)	PFLE: 0.1 $\mathsf{\mu }$s	1 × 10${}^{-4}$
	PFLE: 1 $\mathsf{\mu }$s	0.003		PFLE: 1 $\mathsf{\mu }$s	3 × 10${}^{-4}$
	PFLE: 2 $\mathsf{\mu }$s	0.003		PFLE: 2 $\mathsf{\mu }$s	1 × 10${}^{-4}$
	PFLE: 4 $\mathsf{\mu }$s	0.005		PFLE: 4 $\mathsf{\mu }$s	2 × 10${}^{-4}$
2 and 3 (Figure 14)	PFLE: 0.1 $\mathsf{\mu }$s	0.005	2 and 3 (Figure 15)	PFLE: 0.1 $\mathsf{\mu }$s	1 × 10${}^{-4}$
	PFLE: 1 $\mathsf{\mu }$s	0.003		PFLE: 1 $\mathsf{\mu }$s	2 × 10${}^{-4}$
	PFLE: 2 $\mathsf{\mu }$s	0.003		PFLE: 2 $\mathsf{\mu }$s	4 × 10${}^{-4}$
	PFLE: 4 $\mathsf{\mu }$s	0.005		PFLE: 4 $\mathsf{\mu }$s	1 × 10${}^{-4}$

Table 2

NRMSD of the voltage and current curves obtained with the PFLE with respect to the voltage and current curves obtained with the ULM and a time step of 1 $\mathsf{\mu }$s: Figure 16 and Figure 17.

Voltage			Current
Phases	Model: $\Delta$ t	NRMSD	Phases	Model: $\Delta$ t	NRMSD
1 (Figure 16)	PFLE 0.1 $\mathsf{\mu }$s	0.004	1 (Figure 17)	PFLE 0.1 $\mathsf{\mu }$s	2 × 10${}^{-4}$
	PFLE 1 $\mathsf{\mu }$s	0.003		PFLE 1 $\mathsf{\mu }$s	3 × 10${}^{-4}$
	PFLE 2 $\mathsf{\mu }$s	0.003		PFLE 2 $\mathsf{\mu }$s	4 × 10${}^{-4}$
	PFLE 4 $\mathsf{\mu }$s	0.005		PFLE 4 $\mathsf{\mu }$s	1 × 10${}^{-4}$
2 and 3 (Figure 16)	2	0.005	2 and 3 (Figure 17)	PFLE 0.1 $\mathsf{\mu }$s	2 × 10${}^{-4}$
	PFLE 1 $\mathsf{\mu }$s	0.003		PFLE 1 $\mathsf{\mu }$s	1 × 10${}^{-4}$
	PFLE 2 $\mathsf{\mu }$s	0.003		PFLE 2 $\mathsf{\mu }$s	2 × 10${}^{-4}$
	PFLE 4 $\mathsf{\mu }$s	0.005		PFLE 4 $\mathsf{\mu }$s	3 × 10${}^{-4}$

Table 3

NRMSD of the voltage and current curves obtained with the PFLE with respect to the voltage and current curves obtained with the ULM and a time step of 1 $\mathsf{\mu }$s: Figure 18 and Figure 19.

Voltage			Current
Phases	Model: $\Delta$ t	NRMSD	Phases	Model: $\Delta$ t	NRMSD
1 (Figure 18)	PFLE 0.1 $\mathsf{\mu }$s	0.004	1 (Figure 19)	PFLE 0.1 $\mathsf{\mu }$s	2 × 10${}^{-4}$
	PFLE 1 $\mathsf{\mu }$s	0.003		PFLE 1 $\mathsf{\mu }$s	3 × 10${}^{-4}$
	PFLE 2 $\mathsf{\mu }$s	0.003		PFLE 2 $\mathsf{\mu }$s	2 × 10${}^{-4}$
	PFLE 4 $\mathsf{\mu }$s	0.005		PFLE 4 $\mathsf{\mu }$s	1 × 10${}^{-4}$
2 and 3 (Figure 18)	2	0.005	2 and 3 (Figure 19)	PFLE 0.1 $\mathsf{\mu }$s	42 × 10${}^{-4}$
	PFLE 1 $\mathsf{\mu }$s	0.003		PFLE 1 $\mathsf{\mu }$s	1 × 10${}^{-4}$
	PFLE 2 $\mathsf{\mu }$s	0.003		PFLE 2 $\mathsf{\mu }$s	2 × 10${}^{-4}$
	PFLE 4 $\mathsf{\mu }$s	0.005		PFLE 4 $\mathsf{\mu }$s	4 × 10${}^{-4}$

Table 4

NRMSD of the voltage and current curves obtained with the PFLE with respect to the voltage and current curves obtained with the ULM and a time step of 1 $\mathsf{\mu }$s: Figure 22 and Figure 23.

Voltage			Current
Phases	Model	NRMSD	Phases	Model	NRMSD
1 (Figure 22)	PFLE 0.1 $\mathsf{\mu }$s	0.004	1 (Figure 23)	PFLE 0.1 $\mathsf{\mu }$s	3 × 10${}^{-4}$
	PFLE 1 $\mathsf{\mu }$s	0.003		PFLE 1 $\mathsf{\mu }$s	3 × 10${}^{-4}$
	PFLE 2 $\mathsf{\mu }$s	0.003		PFLE 2 $\mathsf{\mu }$s	1 × 10${}^{-4}$
	PFLE 4 $\mathsf{\mu }$s	0.005		PFLE 4 $\mathsf{\mu }$s	1 × 10${}^{-4}$
2 and 3 (Figure 22)	PFLE 0.1 $\mathsf{\mu }$s	0.005	2 and 3 (Figure 23)	PFLE 0.1 $\mathsf{\mu }$s	2 × 10${}^{-4}$
	PFLE 1 $\mathsf{\mu }$s	0.003		PFLE 1 $\mathsf{\mu }$s	2 × 10${}^{-4}$
	PFLE 2 $\mathsf{\mu }$s	0.003		PFLE 2 $\mathsf{\mu }$s	1 × 10${}^{-4}$
	PFLE 4 $\mathsf{\mu }$s	0.005		PFLE 4 $\mathsf{\mu }$s	4 × 10${}^{-4}$

References

1. Marti, J.R. Accurate Modeling of Frequency-Dependent Transmission Lines in Electromagnetic Transient Simulations. IEEE Power Eng. Rev.; 1982; 2, pp. 29-30. [DOI: https://dx.doi.org/10.1109/MPER.1982.5519686]

2. Laurent Dubé European EMTP-ATP Users Group e.V. Users Guide to Models in ATP. April 1999; Available online: https://pdfcoffee.com/users-guide-to-models-in-atp-acknowledgments-pdf-free.html (accessed on 30 January 2022).

3. Morched, A.S.; Gustavsen, B.; Tartibi, M. A universal model for accurate calculation of electromagnetic transients on overhead lines and underground cables. IEEE Trans. Power Deliv.; 1999; 14, pp. 1032-1037. [DOI: https://dx.doi.org/10.1109/61.772350]

4. Manitoba HVDC Research Centre Inc. PSCAD-EMTDC Version 4.5 The Professional’s Tool for Electromagnetic Transients Simulation; Manitoba HVDC Research Centre, Inc.: Winnipeg, MB, Canada, 2008.

5. EMTP-RV. 2021; Available online: http://emtp.com/ (accessed on 1 February 2022).

6. Dommel, H.W. Digital Computer Solution of Electromagnetic Transients in Single-and Multiphase Networks. IEEE Trans. Power Appar. Syst.; 1969; PAS-88, pp. 388-399. [DOI: https://dx.doi.org/10.1109/TPAS.1969.292459]

7. De Silva, H.M.J.; Kariyawasam, K.K.M.A.; Gole, A.M.; Nordstrom, J.E. Accurate time domain simulation of frequency dependent transmission line models for large time steps. Proceedings of the International Conference on Power Systems Transients (IPST); Vancouver, BC, Canada, 18–20 July 2013.

8. Gustavsen, B.; Semlyen, A. Admittance-Based Modeling of Transmission Lines by a Folded Line Equivalent. IEEE Trans. Power Deliv.; 2009; 24, pp. 231-239. [DOI: https://dx.doi.org/10.1109/TPWRD.2008.2002960]

9. Gustavsen, B.; Semlyen, A. Rational approximation of frequency domain responses by vector fitting. IEEE Trans. Power Deliv.; 1999; 14, pp. 1052-1061. [DOI: https://dx.doi.org/10.1109/61.772353]

10. Gustavsen, B. Improving the pole relocating properties of vector fitting. IEEE Trans. Power Deliv.; 2006; 21, pp. 1587-1592. [DOI: https://dx.doi.org/10.1109/TPWRD.2005.860281]

11. Simscape Electrical User’s Guide (Specialized Power Systems). Hydro-Quebec and Math Works. 2022; Available online: https://www.mathworks.com/help/physmod/sps/ (accessed on 1 February 2022).

12. Morched, A.S.; Ottevangers, J.H.; Marti, L. Multi-port frequency dependent network equivalents for the EMTP. IEEE Trans. Power Deliv.; 1993; 8, pp. 1402-1412. [DOI: https://dx.doi.org/10.1109/61.252667]

13. Semlyen, A.; Gustavsen, B. Phase-Domain Transmission-Line Modeling With Enforcement of Symmetry Via The Propagated Characteristic Admittance Matrix. IEEE Trans. Power Deliv.; 2012; 27, pp. 626-631. [DOI: https://dx.doi.org/10.1109/TPWRD.2011.2180405]

14. Gustavsen, B. Validation of frequency-dependent transmission line models. IEEE Trans. Power Deliv.; 2005; 20, pp. 925-933. [DOI: https://dx.doi.org/10.1109/TPWRD.2004.837676]

15. Tavares, M.C.; Pissolato, J.; Portela, C.M. New mode domain multiphase transmission line model-transformation matrix modeling. Proceedings of the POWERCON, International Conference on Power System Technology, Proceedings (Cat. No.98EX151); Beijing, China, 18–21 August 1998; Volume 2, pp. 855-859. [DOI: https://dx.doi.org/10.1109/ICPST.1998.729206]

16. Colqui, J.S.L.; Timaná, L.; Caballero, P.T.; Kurokawa, S.; Filho, J.P. Implementation of modal domain transmission line models in the ATP software. IEEE Access; 2022; 10, pp. 15924-15934. [DOI: https://dx.doi.org/10.1109/ACCESS.2022.3146880]

17. Gustavsen, B.; De Silva, H.M.J. Inclusion of Rational Models in an Electromagnetic Transients Program: Y-Parameters, Z-Parameters, S-Parameters, Transfer Functions. IEEE Trans. Power Deliv.; 2013; 28, pp. 1164-1174. [DOI: https://dx.doi.org/10.1109/TPWRD.2013.2247067]

18. Gustavsen, B. Computer code for rational approximation of frequency dependent admittance matrices. IEEE Trans. Power Deliv.; 2002; 17, pp. 1093-1098. [DOI: https://dx.doi.org/10.1109/TPWRD.2002.803829]

19. Gustavsen, B. Fast Passivity Enforcement for Pole-Residue Models by Perturbation of Residue Matrix Eigenvalues. IEEE Trans. Power Deliv.; 2008; 23, pp. 2278-2285. [DOI: https://dx.doi.org/10.1109/TPWRD.2008.919027]

20. Gustavsen, B. Matrix Fitting Toolbox (User’s Guide and Reference); SINTEF Energy Research: Trondheim, Norway, 2009.

21. Gustavsen, B. User’s Guide for Vectfit3.m (Fast Relaxed Vector Fitting); SINTEF Energy Research: Trondheim, Norway, 2008.

22. Boyd, S.; Chua, L.O. On the passivity criterion for LTI n-ports. Circuit Theory Appl.; 1982; 10, pp. 323-333. [DOI: https://dx.doi.org/10.1002/cta.4490100404]