Chapter 3
Modeling of power system in PSCAD/EMTDC program

3.1 Introduction PSCAD/EMTDC

Electromagnetic Transient including DC (EMTDC) is an implementation of
the EMTP – type, initially designed for the solution of ac-dc power systems analysis.
It includes models of various physical components contained in functional modules.
It is a related software packages that provide the user with a very flexible power
systems electromagnetic transient simulation tool. EMTDC is the software, which
actually performs the electromagnetic transient analysis on the user defined power
system. While Power Systems Computer Aided Design (PSCAD) is the software that
uses the graphical interface. The interface between PSCAD and EMTDC is shown in
Figure 3.1 [21-22].

File Manager

Draft T-Line/Cable Run Time Multi Plot

Information
Graphic User
Interface PSCAD

EMTDC Progam
Data
EMTDC
DSDYN
Program
DSOUT

Output Files
 
Figure 3.1 PSCAD graphic user interface to EMTDC [21].
 

  37

The general steps that would be required to create a new power system model,
perform a simulation, analyze results and briefly description of each module is shown
as following:
- File Manager: This module allows the user organize many projects and cases
which ultimately will be created into the database.
- Draft: This module allows the user to graphically sketch power system circuit
to be simulated. And all parameters associated with the particular component are
entered into the special menu associated with the component’s icon.
- T-line/Cable: Two special modules have been written to handle the complex
computations which is required to generate data for the overhead transmission line
modules and cable modules.
- Run Time: This module allows the user to monitor and interact with the power
system simulation by two way, one is the EMTDC simulation that runs until
simulation finish time is reached and the other one is real time digital simulator
( RTDS). Simulation runs until the operator decides to stop it. Graphical icons of
meters, sliders and also the simulation by initiating a fault sequence through a push
buttons are used on interface. The run time cover the operation of EMTDC with on
line graphical output displays under run time. When the data input file and the
DSDYN and DSOUT subroutines have been assembled automatically by Draft or
manually by the user, Run time will compile and link DSDYN and DSOUT
subroutines to EMTDC, run, display output, allow for SNAPSHOT files to be written,
enable restarts from SNAPSOT files and provides for on line adjustment for
parameters during the run with graphical icons of potentiometers, push buttons and
switches.
- Multi Plot: This module is used to plot and directly analyze (by Fourier
analysis, rms values, harmonic distortion etc.) the data generated by the EMTDC
simulation and allows multiple graphs to be combined and displayed on a common
page [21].
3.1.1 Typical Studies Using PSCAD/EMTDC
EMTDC is software based electromagnetic transient simulation tool which can
be accessed through PSCAD. The user is able to assemble virtually and conceive
power system by using building blocks available from supplied libraries, or from user
 

  38

generated components. One of EMTDC’s strengths lines is its ability to model

complex power systems which may include HVDC transmission schemes and their
associated controls with relative ease. The simulation of studies routinely conducted
using EMTDC include [23-24]:
1. General power system electromagnetic transient studies
2. DC transmission configuration and controls
3. Effect of FACTS devises
4. Synchronous and induction machine torsion effects and self excitation
5. Static compensators
6. Non-linear control systems
7. Transmission line and cables
8. Transformer saturation studies such as ferroresonances and core saturation
instability
9. Development of new control system concepts
10. Steep front surge analysis.
3.1.2 What PSCAD/EMTDC can be used for
One of the ways understanding the behavior of a complicated system is to
study its response when subjected to disturbances or parametric variations. Simulation
is one way of producing these responses. In power systems, these responses can be
studied by observing either the time domain instantaneous values, time domain rms
values, or the frequency components of response.
PSCAD/EMTDC is the most suitable of simulating the time domain
instantaneous response, also popularly known as electromagnetic transients of
electrical systems. Following are some of the common components used in system
studied using PSCAD/EMTDC [22-24]:
• Resistors (R), inductors (L) and capacitors (C)
• Mutually coupled windings such as transformer
• Distributed frequency dependent transmission lines and cables

• Current and voltage sources

• Switches and breakers
• Diodes, thyristors and GTO’s
 

  39

• Analogue and digital control functions

• AC machines, exciters, governors, stabilizers and inertial models
• Meters and measuring functions
• Generic DC and AC controls
• HVDC, SVC, and other FACTS devices
The power of EMTDC is tremendously increased by its state of the art PSCAD
graphical user interface. PSCAD allows the user to assemble the circuit, run the
simulation, analyze the results, and manage the data in a completely integrated
graphical environment. Together, PSCAD/EMTDC is a very efficient tool for
electromagnetic transient simulation of power systems.
3.2 Equipments Modeling in PSCAD/EMTDC Program
3.2.1 Transmission Line Modeling
There are three basic transmission line modeling techniques in EMTDC: PI
sections, the Bergeron Model, and Frequency-Dependent Line Models. The
requirements of the study will determine which of the three models is suitable.
A simple Pi section model will give the correct fundamental impedance, but
cannot accurately represent other frequencies unless many sections are used (which is
inefficient). It also cannot accurately represent the frequency dependent parameters of
a line (such as skin effect). It is suitable for very short lines where the traveling wave
models cannot be used.
The Bergeron model represents the L and C elements of a Pi section in a
distributed manner (not using lumped parameters like Pi sections). It is roughly
equivalent to use an infinite number of Pi sections except that the resistance is lumped
(1/2 in the middle of the line, ¼ at each end). Like Pi sections, it also accurately
represents the fundamental frequency. It also represents impedances at other
frequencies, except that the losses do not change. This model is suitable for studies
where the fundamental frequency load-flow is most important (i.e. relay studies).
The Frequency-Dependent Line Model represents the frequency dependence of
all parameters. This model takes longer to run than the Bergeron model, but is
necessary for studies requiring a very detailed representation of the line over a wide
frequency range. Frequency dependent models can be solved using modal techniques
or using the more advanced phase domain techniques.
 

  40

The most accurate of these is the frequency dependent model, which represents
all frequency dependent effects of a transmission line, and should be used whenever
in doubt. When using Bergeron model, impedance/admittance data can also be
entered directly to define the transmission corridor, for frequency dependent models
detailed conductor information (i.e. line geometry, conductor radius) must be given.
This thesis have been chosen the frequency dependent model for the modeling
of 500 kV transmission system Nam Theun 2 – Roi Et 2 , the description of this
model shown in chapter 2 section 2.2.4.2 [22-24].
3.2.1.1 The Frequency Dependent Line Model
The Frequency Dependent Line Model is based on the theory developed in
[9]. In order to arrive at the time domain formation of the line equations, it is
convenient to first work in the frequency domain, where an exact solution for a given
frequency is possible.
Consider the following circuit of a transmission line as seen from the
terminations: as show in Figure 3.2.

i k   (ω  )   im  (ω )

V k  (ω )   V m  (ω )

Figure 3.2 The traveling wave frequency dependent line model [22].

For a given frequency, the voltage and currents at one end of the line may be
represented in terms of the voltage and current at other end by :
Vk (ω ) = cosh[γ (ω ).L ].Vm (ω ) − Z c (ω ).sinh[γ (ω ).L ].im (ω ) ,   (3.1)

sinh[γ (ω ). L]
ik (ω ) = .Vm (ω ) − cosh[γ (ω ).L ].im (ω ) ,   (3.2)
 Z c (ω )
 

  41

where:

γ (ω ) = Y (ω ).Z (ω )    : Propagation constant

 Z c (ω ) = Z (ω ) / Y (ω )   : Surge impedance

Y (ω ) = G + jω C     : Shunt admittance of the line

 Z (ω ) = R + jω L   : Series impedance of the line

G : Shunt conductance of the line

C : Capacitance of the line
By introducing the forward and back ward traveling wave functions F k   and  Bk   

in the node k  from Figure. 3.3

Fk (ω ) = Vk (ω ) + Z C (ω ).ik (ω )  ,   (3.3)

 Bk (ω ) = Vk (ω ) − ZC (ω ).ik (ω )  ,   (3.4)

and similarly at node m.

Fm (ω ) = Vm (ω ) + Z c (ω ).im (ω ) ,   (3.5)

 Bk (ω ) = Vm (ω ) − Z c (ω ).im (ω ) .   (3.6)

Substituting equation (3.4) into equation (3.3) gives:

Fk (ω) = 2.Vk (ω) − Bk  (ω ) , (3.7)

And similarly,
Fm (ω ) = 2.Vm (ω ) − Bm (ω ) , (3.8)

For equation ( 3.1) and (3.2) ( and their equivalents at node m ) can now be expressed
in terms of the forward and backward traveling functions.
 Bk (ω ) = A(ω ).F m (ω ) , (3.9)

 Bm (ω ) = A(ω ).F k  (ω ) , (3.10)

where :
1
 A(ω ) = = e −γ (ω ) L , (3.11)
cosh[γ (ω ). L].sinh[γ (ω ).L
 ]

 A( ω )  is known as the propagation constant and is a complex number. The real part

α  is the attenuation constant, and the imaginary part  β  is phase constant.
 

  42

im (ω )
 Z C   (ω )    Z C   (ω )

V k  (ω )  Bk  (ω )     Bm (ω ) V m (ω )

 
Figure 3.3 Frequency Dependent equivalent line model [22].

Equation ( 3.9 ) and ( 3.10 ) relate the source  Bk  (ω )  with the m side quantities

V m (ω )  and  Bm (ω ) .

Introducing in node k .

 Bk (ω ) = A(ω ).[2.Vm (ω ) − Bm (ω )] ,   ( 3.12 )

And similarly in node m 

 Bm (ω ) = A(ω ).[2.Vk (ω ) − Bk  (ω )] ,   ( 3.13 )

Equation ( 3.11 ) and ( 3. 13 ) are represented in Figure 3.3 in the time domain.
The problem however, is the multiplication in equation ( 3.13 ). A multiplication of
frequency domain becomes convolution in the time domain.
t 

 A(ω ).Bm (ω )  ↔ ∫

τ 
A(u ).Bm (t − u ).du . ( 3.14 )

The equation (3.14) is traveling time (t) because an impulse on one end of line
will not reach the other end until τ second. The traveling time is calculated using the
imaginary term  β    of the propagation constant. This equation is still in a very

inconvenient form for a time domain solution because with each time step, more and
more terms of the convolution integral must be evaluated. Fortunately, the
convolution can be computed in a recursive form suitable for time domain solution.
With the reference to the PSCAD/EMTDC line model, the selected transmission
line models are transposed frequency dependent phase model based on traveling time
and characteristic impedance of the line. The 500 kV double circuit lines between
Nam Theun 2 (NTN2) to Roi Et 2 (RE2) substations with 304.5 km length of the line
is divided into two parts. The 135 km line from NTN2 – SVJ with 795 MCM ACSR
of conductor is divided into 5 sections, the 169.5 km line from SVJ-RE2 with 1,272
 

  43

MCM ACSR of conductor is divided into 13 sections. It is illustrated in Appendix B,

Figure B.2.2 - B.2.3.
The input data required by each section for example; section length, phase
transformation data, tower configuration and ground conductivity are shown in  table:
3.1.
Table: 3.1 Data input for transmission line model.

♦Tower configuration ♦ Symbol

- Relative X position of tower center on right of way (m)
- Shunt conductance ( s/m)
- Number of conductors
- Show graphic of conductor sag ( Yes/No)
- Is this circuit ideal transpose ( Yes/No)
- How many ground wires ( Yes/No)
- Estimate ground wires ( Yes/No)
♦ Conductors data
- Conductor radius (m)
- Conductor DC resistance ( Ω /km)
- Sag of all conductor ( m)
- Number of sub-conductors in a bundle
- Bundle configuration ( Symmetrical or non symmetrical ) ( Yes/No)
- Show bundle graphics ( Yes/No)
- Horizontal coordinate and height distance ( X/Y)
- Phase/Node connection information
♦ Ground wires data
- Number of ground wires
- Ground wires radius (m)
- Ground wires DC resistance ( Ω /km)
- Sag of all ground wires (m)
- Horizontal coordinate and height distance ( X/Y)
♦ Other data
- Length of line (km)
- Tline interface component  and Tline info Component
- Steady state frequency (Hz)
- Ground resistivity (Ω)
 

  44

The data preparation and input data in the selected transmission model can be
explained. For the 500 kV double circuits from Nam Theun 2 substation to
Savannakhet junction (NTN2 – SVJ), the representation is shown in Figure 3.5 with
the following step.
-  Select the Tline Interface Component and TlineInfo Component from
Master Library
-  Name the circuit as NTN2-SVJ for both Tline Interface Component and
TlineInfo Component
-  Assign the line parameter for TlineInfo Component and tower
configuration according to the existing of 500 kV transmission line Nam
Theun 2 – Savannakhet junction (NTN 2 – RE 2).

Tline Interface Component at the sending end and the receiving end of the 500 kV
transmission line Nam Theun 2 – Savannakhet junction
 

  45

Edit parameter of Tline Info Component for the 500 kV

Figure 3.4 Tline Interface and Tline Info Conponent for Nam Theun 2 – Savannakhet
Junction model.
The 500 kV double circuits from Savannakhet junction to Roi Et 2 substation,
(SVJ – RE2) can also be represented as shown in Figure 3.5 with the similar step.

Tline Interface Component at sending end and receiving end of the 500 kV
double circuits transmission line Savannakhet junction – Roi Et 2
 

  46

Edit parameter for Tline Info SVJ-RE2

Figure 3.5 Tline Interface and Tline Info Component for Savannakhet Junction – Roi
Et 2 model.

3.2.2 Transformer modeling

The thesis there is two selected transformer models. The single phase
model is used for single phase step up power transformer at Nam Theun 2 substation
and step up power from generator 18 kV connecting with single phase step up power
transformer 18/ 525 kV. The three phases model is used for three phases step down
power autotransformer at Roi Et 2 substation and step down power from line 525 kV
- 242 kV-22 kV via the Roi Et 2 substation.
Transformer models as shown in Figure 3.6 are represented in detail with the
following information required for the PSCAD/ EMTDC transformer components:
- Transformer MVA rating.
- Winding configuration and winding voltage.
- Transformer tap change range and normal setting.
 

  47

- Leakage reactance between windings (this information can be obtained from

the load flow data used in system operating studies, if not readily available
from the name plate data ).
- Knee point of transformer core saturation characteristic in per unit of rate flux
or voltage.
- Estimated saturated air core reactance of transformer and winding it is based
on.
- Based operation frequency

Figure 3.6 Single and three phases transformer model in PSCAD/EMTDC.

3.2.3 Generator modeling

The generator can be represented by the infinite source series with
subtransient impedance matrix. The subtransient matrix can be built by the same
procedure as the three phases transformer as following:
- Positive sequence impedance (X1)
- Negative sequence impedance (X2)
- Zero sequence impedance (Xo)
- Base MVA (three phases)
- Base voltage (L-L rms)
- Source impedance type
 

  48

Build the matrix [X] from the short-circuit input reactance. The matrix contains
3×3 sub-matrixes of the form.

  X
  S
Xm X m 
  X XS X m
  , (3.15)
 m 
  X m Xm X S    

where  X S   is the self reactance of each phase, and  X m  is the mutual reactance among

the three phases. As in any other network three-phase network component, from
Figure 3.7, these self and mutual reactances are related to the positive and zero
sequence values  X 1 and  X 0 by:

 ( X 0 + 2 X 1 )
 X S  = , (3.16)
3
 ( X 0 − X 1 )
 X m = . (3.17)
3
 

Figure 3.7  Three phases X- matrix source.

3.2.4 Surge Arrester Modeling
Zinc oxide or metal oxide surge arresters are highly non-linear resistance,
with an almost infinite slope in the normal voltage region and almost horizontal slope
in the overvoltage protection region. The common form of the exponential function,
describing the protective characteristic of arrester, is as follows:
m
 V  
 I = k 
  .   , (3.18)
 V ref  
where V ref  is the reference voltage used to prevent the numerical overflow of the

exponential term for large voltage values, k   and m are constant values ( typical value
for m equal 20 to 30 ).
 

  49

Surge arrester is characterized with a nonlinear voltage versus resistance

characteristic. The model will be presented by entering points on the V-I
characteristic, or the default characteristic may be applied. It is suitable for designing
switching surge, temporary overvoltage protection as shown in Figure 3.8. 

Figure 3.8 Surge arrester model in PSCAD/EMTDC.

3.2.5 Circuit Breaker Modeling

The circuit breakers are modeled as simple time controlled switches in the
PSCAD/EMTDC program as shown in Figure 3.9. The program allows various
options to vary the closing time ranging from one-shot deterministic closings to
multiple shot statistical closings. For the statistical switching simulations, the
characteristics of the mechanical operating parts are modeled by a random pole
spreading with a standard deviation for the three phase breaker operating times.
Randomly generated switch closing times varying are also used to account for random
operation at any point on the power frequency waveform.
 

  50

The component allows simulation of single phase circuit breaker operation. The
“Open” and “Close” resistance of the breaker is specified along with its initial state. If
the input is 0, the breaker will close. If the input 1, the breaker will open. The breaker
current may be labeled and monitored via an output channel if desired. A pre-insertion
and post-insertion resistor may be represented if required. The data input for circuit
breaker model is shown in table: 3.2.

Table: 3.2 Parameter input for circuit breaker model.

♦ Configuration ♦ Symbol
- Breaker name ( the name should be the same as the output signal
label of the breaker logic component.
- Open possible if current flowing ( Yes/No)
- Use pr-insertion resistance ( Yes/No)

♦ Breaker main data

- Breaker close resistance (Ω)
- Breaker open resistance (Enter large value than close resistance) (Ω)
♦ Pre-Insertion data
- Pre-insertion resistance :
- Time delay for closing breaker (sec.)
- Time delay for bypassing pre-insertion (sec.)

Figure 3.9 Circuit breaker models in PSCAD/EMTDC.

 

  51

3.2.6 Shunt Reactor Modeling

Shunt Reactors in the studied 500 kV system are modeled by simple
lumped element. The locations of shunt reactors are identified. The neutral reactors
are installed in the neutral of the shunt reactors as presented in Figure 3.10.

Figure 3.10 Configuration of line connected shunt reactors with a neutral reactor.

The shunt reactance value of 3 ×55 Mvar line shunt reactor at Nam Theun 2 side
per circuit and 2×55 Mvar line shunt at Roi Et 2 side per circuit can be calculated
from the rate capacity with based voltage 525 kV as following equation 3.19:
For shunt reactor rate of Mvar.

(Voltage[kV ])2
Reactance  X  = . (3.19)
 Base( Mva r)

3.2.7 Load Modeling

The external grid connecting with autotransformer 525/242/22 kV at Roi
Et 2 substation are considered in load model. The external equivalence network is
employed to reduced the rest of the Roi Et 2 substation 500/230 kV. The complete
external network will be reduced to a voltage source behind a Thevenin equivalent
matrix of Zeq. The Zeq is obtained by initializing a three phases fault at the respective
bus and obtaining the sequence impedances at the faulted bus. The equivalent source
from Appendix A.9 is used for this study. The magnitude and phase angles of all
voltage sources are found from a load flow program.
3.3 Initial Condition of Program
There are generally two ways to start a simulation: start from time t = 0.0 with
no initial conditions or start with pre-calculated initial conditions imposed on some or
all elements. In this study, starting a simulation with initial conditions is achieved by
 

  52

using a Snapshot File. This can be accomplished in program by taking a snapshot at a

time just previous to the switch opening. Subsequent runs could be started from this
snapshot file, which would have stored, among other values, the current flowing in the
inductor at the time of the snapshot. All studied cases start at steady state condition by
the used of snapshot feature.
3.4 The Accuracy of the Models
To analysis of the studied system with the simulation, the models must be
verified. All these sophistications in the simulation techniques increase considerably
the complexity and potential value of the studies, but this of little practical use
without accurate “input” data related to the transient behavior of actual system
elements. Therefore, the caution of the model is very important and in this thesis, the
accuracy of the model is checked for validity with the following cases:
- Comparison with TNA case study for line energizing of the Mae Moh to Tha
Ta Ko circuit 1 [6], [8].
- Comparison with line energization overvoltage study report of EGAT on the
500 kV transmission line Nam Theun 2 to Roi Et 2 project [2].
3.4.1 Comparison with TNA case study for line energizing of the Mae Moh
to Tha Ta Ko circuit 1  
The switching overvoltage of 500 kV Mae Moh (MM 3) to Tha Ta Ko
(TTK) circuit 1 in case of line energization was taken from TNA case study to
compare with the result. The purpose is to check the accuracy of the selected models
and the reference case (MM3 – TTK) is quite similar to the studied system (NTN2 –
RE2) [6-8]. The network of MM3 – TTK circuit 1 is shown in Figure 3.11

Figure 3.11 Single line diagram Mae Moh – ThaTaKo circuit 1.

 

  53

All equipments are represented in 500 kV transmission line Mae Moh 3 to

ThaTaKo circuit 1 with the following detail:
- Generator model is the basic Thevenin voltage sources, three phases x matrix
and short circuit impedance.
- Transformers modeled: For each transformer is represented in detail with this
exiting information: MVA rating. Winding configuration and voltage, tap
change ranges and normal setting, leakage reactance between windings. Knee
point of transformer core saturation is characteristic in per unit of rate flux or
voltage, and estimated saturated air core reactance.
- Transmission line model: The selected transmission line models are distributed
parameter models based on the traveling time and characteristic impedance of
the line. The section of 500 kV line Mae Moh – ThaTaKo are presented by
untransposed frequency dependent line model.
- Circuit breaker model: the locations of the circuit breakers that will be
switched are identified on the studied system. Other parameters of the circuit
breaker are determined.
- Arresters models: the installed location and rating of surge arresters are
provided. The maximum ratings and particular the energy absorption
capability is determined with these studies.
- Shunt reactors models: the location of shunt reactors is identified. This is
included whether. They are line connected or bus connected.
3.4.2 Results of the Studies
The studied line energization case 500 kV transmission line Mae Moh to
Tha Ta Ko circuit 1 is shown in Figure 3.11.

Phase A
 

  54

1.750 p.u

Phase B

Phase C
( A) (B)
Figure 3.12 MM 3-TTK circuit 1 waveform from PSCAD/EMTDC (A) and TNA (B).

The TNA result and PSCAD/EMTDC result show relatively the similar
waveform with the same condition for breaker operation as shown in Figure 3.12.
3.4.3 Comparison with line energization overvoltage study report of EGAT
on the 500 kV transmission line Nam Theun 2 to Roi Et 2 project
The EMTP study is conducted to investigate the switching overvoltage
occurring during switching of the line circuit breaker on each side of the 500 kV
transmission line between Nam theun 2 project and EGAT’s network at Roi Et 2
substation. The initial study is for EGAT planning development before the project
starting [2]. The DigSILENT software program was used for the study. After
confirmation of the selected model from the previous section, the studied system is
represented with the existing parameter and information related to EGAT studied
case. The simulation is performed in three steps as follows:
Step 1 : Prior to the switching action commencing , EGAT’s system is in a
normal steady state operating condition. All circuit breakers at both ends of
these two circuits are opened and all generators at Nam Theun 2 Power Plant are
off line.
Step 2 : All generators of Nam Theun 2 Power Plant are running at no load
synchronous speed while all line circuit breakers at both ends of two circuits are
still open.
 

  55

Step 3 : One set of line circuit breakers on one circuit of the double circuit lines
at Nam Theun 2 end is closed.

Table: 3.3 Simulation results from DigSILENT and PSCAD/EMTDC program

DigSILENT simulation PSCAD/EMTDC
Detail
of EGAT simulation
Energized from Nam Theun 2 Power Plant
Maximum Overvoltage ( kV) -817.219 -817.030
Point Receiving end at RE 2 Receiving end at RE 2
Phase C C
Time closing of CB ( sec ) 0.020 0.020

The receiving end voltage from EMTDC simulation was recorded and compared
with the corresponding result of DigSILENT simulation. The comparison values are
shown in table: 3.3.
The waveform of DigSILENT and PSCAD/EMTDC simulation are shown in
Figure 3.13 and 3.14. The waveform of the voltage at the receiving end Roi Et 2
substation from both simulations are almost identical with similar highest voltage.

Figure 3.13 The voltage waveform of PSCAD/EMTDC simulation.

Figure 3.14 The voltage waveform of Dig SILENT simulation.

 

  56

3.5 Conclusion for the Accuracy of the Models

The accuracy of the PSCAD/EMTDC modeling program of the system
components as well as lines, transformers, generators, circuit breakers, shunt reactors,
surge arresters, loads are presented by taking into account the saturation of
transformer, shunt reactor and surge arresters. The raw data has to be converted to the
suitable value before input to the program.
The accuracy of the model with TNA case Mae Moh to Tha Ta Ko, the results
from simulation are similar. This comparison is taking the selected equipments model
and methodology of the studied for improvement of the selected model in case of
Nam Theun 2 to Roi Et 2.
The accuracy of the model with comparison of line energization overvoltage
EMTP study report of EGAT on the 500 kV transmission line Nam Theun 2 to Roi Et
2 project, the simulation from PSCAD/EMTDC and DigSILENT program are in good
agreement.
In case of the model verification, it can be concluded that the analysis of
transient overvoltage used PSCAD/EMTDC program in case Mae Moh to Tha Ta Ko
the maximum overvoltage are the same values and comparison with line energization
study report of EGAT on Nam Theun 2 to Roi Et 2. The maximum overvoltages from
both are similarly. The simulation had the identical system configuration as the
reference case. The voltage waveforms are almost identical, the comparisons give
satisfactory results, further confirming the accuracy of the simulation model and used
parameters. Therefore, the selected models in this thesis are accuracy and can be used
for analysis of the proposed system.