Instituto Tecnológico de

Morelia
“José María Morelos y Pavón”

Departamento de
Ingeniería Eléctrica
Modelado de líneas
De transmisión

Juan Francisco Yrena Heredia

N° de control: 09120788

Morelia Michoacán a 8 de julio de 2013

 1.0 Introducción
Podemos definir una línea de transmisión como el elemento que toma potencia desde la
generación y la entrega a la subestación de gran potencia. Esta línea de transmisión está formada
por básicamente por conductores sostenidos por una estructura metálica solida cuidadosamente
diseñada, a diferencia de un conductor como el que podría ser el que está usando para conectar
su ordenador al toma corriente en una línea de transmisión los puntos están tan distantes que ya
no pueden ser consideradas un nodo o un bus como se le llama en sistemas eléctricos de potencia,
en una línea de transición se manifiestan tres parámetros los cuales son: resistencia, inductancia y
capacitancia, lo cual nos puede conducir a modelos de circuitos de tercer orden que se resuelven
como ecuaciones diferenciales, para evitar estos cálculos una línea de transmisión se representa
como una red de dos puertos cuya solución es relativamente fácil. Al momento de modelar una
línea hay que tomar en consideración varias cosas, las más importante s son: Distancia y
frecuencia. Los modelos principalmente usados son el modelo de línea corta que son distancias
menores a 80 km a 60 Hz, modelo PI nominal que comprende distancias entre 80 y 240 km y
modelo de exacto o de línea larga para distancias mayores a 240 km. Ante esto puede surgir la
pregunta ¿Qué es lo que hace que se tomen diferentes modelos y no se generaliza en uno solo?, la
respuesta a esto es principalmente la capacitancia, en una línea corta su efecto es tan poco
significante que puede despreciarse, y como veremos las ecuaciones para una línea corta son las
más sencillas de todas, en línea media la capacitancia puede tomarse como una concentración
solida (modelo T), o distribuida en dos (modelo pi) y tener una respuesta que no pierda
credibilidad por error, pero en una línea larga se debe tomar una distribución uniforme a lo largo
de toda la línea

2.0 Desarrollo:

2.1 Línea corta.

Como se mencionó anteriormente la línea corta tiene una longitud menor a 80 millas, para poder
llegar a sus ecuaciones correspondientes tomaremos en cuenta las siguientes consideraciones: La
relación X/R es muy grande por lo que se puede ignorar la resistencia serie de la línea, además la
admitancia capacitiva en derivación no se toma en cuenta, en el presente documento se tomara
en cuenta la resistencia serie, a continuación de muestra el desarrollo de las ecuaciones para el
modelo de línea corta:

La línea se buscara llevar al equivalente de una red de dos puertos como la que se puede apreciar
en la figura 2.1.1, dicha red de dos puertos tiene una solución de la forma:

…………………………. (1)
……………………..……..(2)

Las ecuaciones (1) y (2) pueden reescribirse en forma matricial como:

* +* +* +……………………………(3)
Los parámetros A, B, C, D dependen de los parámetros de línea (R, L y C) y además cumplen con las
siguientes condiciones:
…………………………………(4)
……………………………………………..(5)

Figura 2.1.1. Red de dos puertos

Reiterando lo que hemos venido comentando, en la modelación de la línea corta no tomaremos

en cuenta el efecto capacitivo en derivación por lo que el circuito de línea quedara como esta
representado en la figura 2.1.2

Figura 2.1.2 Modelo de línea corta.

Definiendo las siguientes variables.

Z= R+jῳL
Y=G+jῳC
Z=zl
Y=gl

Analizando la malla por LTK y uno de los nodos por LCK obtenemos que:

A=D=1 p.u.
B=Z Ω.
C=0 S.

2.2 Modelo de linea media:

2.2.1 Modelo PI nominal.

El modelo PI nominal toma en cuenta el efecto capasitivo y lo divide en dos como se muetra en la
figura 2.2.1.1

Figura 2.2.1.1. Modelo pi nominal

Al igual que el modelo de linea corta buscaremos llevar el modelo a un equivalente de red de dos
puertos. Analizando el equivalente del modelo pi nominal obtenemos las ecuciones (6) y (7):

( ) ………………………………(6)
( ) ( ) ………..………(7)

O en su forma matricial como se denota en (8):

( )
* +[ ]* +……………(8)
( ) ( )

De lo cual resulta evidente que:

A=D=( )
B=Z Ω
C= ( )

2.2.2 Modelo T
También existe el llamado modelo T. Donde la admitancia capacitiva en derivación se concentra en
un en el centro de la línea de transmisión y la impedancia serie se divide en dos:
Figura 2.2.2. Modelo T
La figura 2.2.2 muestra un equivalente monofásico al modelo T del cual tras realizar un análisis
puede demostrarse que:

A=D=ZY/4+1
B=Z(ZY/4+1)
C=Y

2.3 Modelo de línea larga.

Para la solución de la línea de transmisión de cualquier distancia al calcular líneas de 60Hz puesto
a que los parámetros de la línea no están agrupados sino que están distribuidos a lo largo de la
línea de forma uniforme. En la figura 2.3.1 se muestra el equivalente monofásico de la línea con
parámetros distribuidos y se considera un elemento diferencial dx a una distancia x del punto de
recepción de la línea, entonces x dx y x dx son la impedancia en serie y la admitancia en paralelo
de la línea y V e I son fasores que varían con respecto a x.

Figura 2.3.1 Diagrama monofásico de la línea para el modelo exacto o de línea larga.

Al realizar el análisis obtenemos la ecuación (9) que es la solución al voltaje:

…………………………….(9)

La solución de (9) tiene la forma:

Dónde:
A1 y A2 son constantes de integración
=√
Similarmente la solución a la corriente se muerta en (10)

………….……….(10)
De la ecuación (10) vemos que:
√ ……………………..…………(11)
√
Entonces definimos:
√ …………………………………..…….(12)
Y a esta nueva constante se le conoce como impedancia característica:
Solucionando las ecuaciones y aplicando las identidades de Euler vemos el resultado expresado en
su forma matricial en la ecuación (13)

[ ][ ]* +…….(14)

3.0 Método computacional

Los análisis mediante el uso de un ordenador mejoran la precisión y dan resultados rápidos, estos
programas incluso más complejos pueden obtener desde el cálculo de parámetros, modelados de
línea de alta precisión, cálculos de flujos, Ybus Zbus y fallas de cualquier tipo en cualquier punto
del sistema con solo ingresar datos como lo son configuraciones de línea, distancias y la forma en
que están interconectadas ya sea acoplamientos eléctricos o magnéticos. En este caso se dará uso
al tan conocido software MATLAB. El código del programa así como su uso se ilustran a
continuación.

%Solución de modelado de lineas

clc, clear all;
%%Datos de la línea
r=0.036; %Resistencia por unidad de longitud (Ohm/Km)
L=0.8E-3; %Inductancia por unidad de longitud (Henry/Km)
C=0.0112E-6; %Capacitancia por unidad de longitud (Farad/Km)
lon=input('Cual es la longitud de la linea en km: ');
%Longitud de la línea (Km)
f=60; %Frecuencia de operación en Hertz
%%Datos de la carga
Sr3=270E6; %Potencia trifásica consumida por la carga
FP=0.8; %Factor de potencia en atraso en la carga
Vr=325E3; %Voltaje de línea a línea en la carga
%%Cálculo de los parámetros ABCD

Z=(r+j*2*pi*f*L)*lon;
Y=1/(1/(j*2*pi*f*C)/lon);
if lon<=80
disp('El modelo usado es: Modelo de linea corta')

A=1;
 B=Z;
C=0;
D=A;
ABCD=[A B
C D];1
elseif lon<=240
asd=input('Que modelo desea usar Pi nominal (1); Modelo T (2): ');
if asd==1
disp('El modelo usado es: Modelo Pi nominal')

%modelo pi
A=1+Z*Y/2;
B=Z;
C=Y*(1+Z*Y/4);
D=A;
ABCD=[A B
C D];
else %modelo T
disp('El modelo usado es: Modelo T')

A=Z*Y/2+1;
B=Y*(Z*Y/4+1);
C=Y;
D=A;
ABCD=[A B
C D];
end
else
disp('El modelo usado es: Modelo de linea larga')

gamma=sqrt(Z*Y/lon^2);
Zc=sqrt(Z/Y);
A=cosh(gamma*lon);
B=Zc*sinh(gamma*lon);
C=1/Zc*sinh(gamma*lon);
D=A;
ABCD=[A B
C D];
end
%%%% de aki pa abajo es en espera_________________________________
%%Cálculo del voltaje y la corriente de envío
Vrp=Vr/sqrt(3);
Srp=Sr3/3*(FP+j*sin(acos(FP)));
Irp=conj(Srp)/Vrp;
VsIs=ABCD*[Vrp; Irp];
disp('El voltaje de envío de línea a neutro (en Volts)')
Vsp=VsIs(1)
disp('Magnitud del voltaje de envío de línea a línea (en Volts)')
mag_Vsll=sqrt(3)*abs(Vsp)
disp('La magnitud de la corriente de envío (en Amperes)')
Isp=VsIs(2);
mag_Isp=abs(Isp)
disp('La potencia aparente de envío (en VA)')
Ss3=3*Vsp*conj(Isp);
mag_Ss3=abs(Ss3)
disp('Factor de potencia de envío')
FPs=cos(angle(Ss3))
disp('Regulacion de voltaje (%)')
RV=(abs(Vsp/A)-abs(Vrp))/abs(Vrp)*100
disp('Eficiencia (%)')
n=real(3*Srp)/real(Ss3)

Uso del programa:

El presente programa está diseñado para el voltaje de envió (de línea a neutro), Corriente de
envió, potencia aparente de envió, factor de potencia de envió, regulación de voltaje y eficiencia
de la línea a partir de los siguientes datos: r=0.036 Ω/km, L=0.8x10^-3 H/km, C=0.0112x10^-6
F/km, f=60 Hz, Sr3=270x10^6 VA, FP=0.8 y Vr=325x10^3 v.
Una vez se corre el programa este solicita se ingrese la distancia de la línea como lo muestra la
figura 3.1

Figura 3.1. Ejecución del programa para el cálculo de parámetros.

El segundo dato que nos pide es cuando la distancia se encuentra entre los 80 y 240 km (longitud
media de línea) y esta vez solicita al operador que indique el método que se decía utilizar como se
indica en la figura 3.2.
Figura 3.2. Selección del modelo de línea media.

Se realizaran varias corridas con diferentes distancias y métodos para ilustrar los resultados que el
programa nos ofrece, así como realizarán comparaciones.

Distancia Método Voltaje de envió Corriente de Potencia aparente de Factor de Regulacion de Eficiencia error
(km) usado (magnitud) envió envió mag. (VA) potencia de voltaje (%) (%)
(v) magnitud envio
(A)
10 Línea 3.2675e+005 479.6448 2.7145e+008 0.7966 0.5377 99.8851 6.3418e-5
corta
20 Línea 3.2851e+005 479.6448 2.7291e+008 0.7933 1.0786 99.7705 2.5269e-
corta 004
50 Línea 3.3384e+005 479.6448 2.7734e+008 0.7833 2.7193 99.4281 0.0016
corta
70 Línea 3.3744e+005 479.6448 2.8033e+008 0.7767 3.8277 99.2012 0.0030
corta
80 Línea 3.3925e+005 479.6448 2.8184e+008 0.7734 4.3862 99.0882 0.0040
corta
85 Pi 3.3868e+005 440.4341 2.5836e+008 0.8436 4.6910 99.1078 0.0045
nominal
85 T 3.3864e+005 440.3731 2.5830e+008 0.8438 4.6787 99.1071 4.2459e-
005
100 Pi 3.4086e+005 433.8581 2.5614e+008 0.8521 5.5517 98.9666 0.0061
nominal
100 T 3.4079e+005 433.7645 2.5604e+008 0.8524 5.5315 98.9653 7.0021e-
005
150 Pi 3.4761e+005 413.0088 2.4866e+008 0.8816 8.5103 98.5266 0.0136
nominal
150 T 3.4738e+005 412.7627 2.4835e+008 0.8828 8.4404 98.5206 2.4621e-
004
200 Pi 3.5357e+005 394.2100 2.4141e+008 0.9118 11.6326 98.1283 0.0239
nominal
200 T 3.5303e+005 393.8151 2.4080e+008 0.9143 11.4621 98.1106 6.0705e-
004
235 Pi 3.5727e+005 382.5542 2.3673e+008 0.9323 13.9357 97.8709 0.0328
nominal
235 T 3.5639e+005 382.1501 2.3589e+008 0.9359 13.6531 97.8382 0.0010
240 Pi 3.5777e+005 381.0026 2.3610e+008 0.9351 14.2736 97.8354 0.0342
nominal
240 T 3.5683e+005 380.6097 2.3523e+008 0.9389 13.9717 97.8000 0.0011
241 Línea 3.5731e+005 380.6502 2.3557e+008 0.9371 14.1342 97.8431 0
larga
300 Línea 3.6207e+005 365.0798 2.2895e+008 0.9681 18.1080 97.4528 0
 larga
400 Línea 3.6718e+005 351.2403 2.2338e+008 0.9984 25.5455 96.8522 0
larga
500 Línea 3.6852e+005 355.3638 2.2683e+008 0.9891 34.1685 96.2752 0
larga
700 Línea 3.6009e+005 412.4462 2.5724e+008 0.8843 57.2755 94.9577 0
larga
1000 Línea 3.2260e+005 559.3044 3.1252e+008 0.7523 130.0290 91.8773 0
larga
1000 Línea 5.3157e+005 479.6448 4.4161e+008 0.5454 63.5595 89.6837 0.6477
corta

Nota para el cálculo del error será el error relativo en el voltaje magnitud descrito por la formula
(15).

……………………….(15)

Tomaremos como valor calculado como el valor que nos provea matlab para cada corrida realizada
y como valor verdadero el valor que nos provea matlab para el cálculo de esa misma distancia
usando el modelo de línea larga (lo cual dará un error de cero para todos los casos donde se
calcule con el método exacto).

Conclusiones:
El uso del programa ilustra cómo es que la línea se comporta para las distancias y observando los
datos se puede uno dar cuenta porque las distancias son estandarizadas en los modelos que
usamos para una línea a 60 Hz, la variación al principio no es mucha entre un modelo aproximado
y el exacto, pero a medida que aumenta la distancia el error generado se hace mayor en especial
en los voltajes como se muestra en la última corrida de la tabla 3.1 tanto el voltaje como el error
son realmente considerables. Despreciar el efecto capacitivo o considerarlo concentrado en ya sea
un punto (modelo T) o en dos (modelo PI) genera errores y aunque en el modelo de línea corta al
llegar a los 80 km no se nota un error tan grande como en el que tenemos al llegar a 240 km en
modelo de línea media recordemos que este modelo también recomienda ignorar la resistencia
serie lo cual generaría un mayor error. Este mismo programa con procedimientos parecidos se
puede usar para determinar también a que distancias deja de ser viable un modelo y se tiene que
cambiar a otro mejor a diferentes frecuencias con solo cambiar la frecuencia de operación, realizar
algunas corridas en el programa y observar donde los errores comienzan a ser consistentes, y por
ultimo concluyo que las distancias son solo consideraciones y estándares, se puede usar un
modelo PI nominal en una distancia de 75 km para aumentar la precisión sin que la dificultad
aumente considerablemente o se puede usar un modelo igualmente PI para una distancia de 245
km que pueden facilitar un análisis rápido cuando así se requiera.

Fuentes de información:
[1] Power system analysis John J. Graiger, William D. Stevenson Jr. Ed Mc Graw hill
[2] Cálculo de parámetros operación de estado estacionario de líneas de transmisión. Lino Coria Cisneros
[3] Sistemas de distribución y análisis de cargas. J Lemus Duarte.