ESTUDIO E IMPLEMENTACIÓN DE

MODELOS DE CARGA EN SISTEMAS

ELÉCTRICOS DE DISTRIBUCIÓN

Sabela González Castro

JULIO 2018

TRABAJO FIN DE GRADO PARA Sabela González Castro

LA OBTENCIÓN DEL TÍTULO DE
GRADUADO EN INGENIERÍA EN DIRECTOR DEL TRABAJO FIN DE GRADO:

TECNOLOGÍAS INDUSTRIALES Rosa María de Castro

Fernández
Agradecimientos

A mi madre, por ser mi pilar fundamental e infalible siempre y haber luchado con cariño
a mi lado en cada paso; a mi padre, por su impecable fe en mi, ánimo y admiración
constantes; a mi hermana, por tornar cualquier sentimiento siempre en alegrı́a (incluso
en las largas horas de biblioteca) y a mis Bio y Bia por su cariño, preocupación y
ánimos. Os quiero y os debo todo.

También a todas esas personas y amigos que me han ayudado sin dudarlo y han hecho
que esta etapa haya sido soportable.
Gracias por acompañarme a ser quien hoy soy.

Por último agradecer a mi tutora Rosa Marı́a su disposición y ayuda al realizar este
Proyecto de forma tan eficiente, y por todo lo que he aprendido.

I
Resumen ejecutivo.

La energı́a eléctrica que se usa a diario es generada en centrales, transportada a conti-

nuación por lı́neas de transmisión y finalmente distribuida a los centros receptores. Para
controlar el sistema eléctrico en su totalidad hay que analizar cada una de estas partes
integradoras.

Las redes eléctricas son desequilibradas y acopladas. Por ello es más complejo su reso-
lución, puesto que el método de componentes simétricas no aporta ventajas a la hora del
estudio y se realiza por componentes de fase.

Con el objetivo de analizar y controlar el comportamiento de los sistemas de distribu-

ción, en este Trabajo de Fin de Grado se van a implementar y estudiar modelos que
permitan trabajar a las redes desequilibradas con un comportamiento determinado.
Conseguir controlar la red para que aporte energı́a de la manera más exacta posible en
relación a lo demandado supone unos beneficios increı́bles. De esta manera se trabajarı́a
con mayor eficiencia, eficacia y minimizando las pérdidas.

Los cargas se pueden clasificar en estáticas y dinámicas. En este Proyecto se trabajará

con cargas estáticas para los tres tipos de análisis que se realizarán.
Ası́ se desarrollarán unos modelos de cargas estáticas que permitan que el sistema trabaje
consumiendo potencia, impedancia o intensidad constante para diferentes situaciones y
análisis.

Para poder establecer relaciones y comparaciones entre los resultados obtenidos de los
modelos, estos se van a implementar sobre una lı́nea de cuatro nudos y un transformador.
Los parámetros de esta lı́nea se calculan con anterioridad y se mantendrán constantes
durante todo el Trabajo.
La lı́nea consta de cuatro nudos, pero los resultados recogidos se medirán en el nudo 4 al
ser el que está conectado a la carga.

II
 Estudio de modelos de carga en sistemas eléctricos de distribución

En primer lugar se hará un análisis de flujo de cargas correspondiente a un determinado

instante. Para este tipo de simulación, Simulink ofrece una carga que permite escoger si
trabaja consumiendo potencia, intensidad o impedancia constante. Se seleccionará cada
opción y se recogerán los datos de tensiones de fase medidas en el nudo 4 con el propósito
de su posterior comparación.

A continuación se procederá al análisis fasorial, que trabaja de manera similar al ejecutar

infinitos flujos de carga para cada instante a lo largo del tiempo. Es decir, como antes se
tenı́a el estudio de un punto de tiempo determinado, en este caso se tiene el de todos los
puntos en régimen permanente.
Para este análisis, como ya no es válida la herramienta de Simulink de la carga, se mode-
larán unas cargas de manera que trabajen según el objetivo de cada modelo.
Partiendo de la base teórica de cada uno de ellos, se obtendrá una intensidad correspon-
diente, la cual se inyectará a la lı́nea. Esto es ası́ ya que en Simulink una carga se simula
como tres fuentes de intensidad más una resistencia en paralelo para evitar indetermina-
ciones. Por eso, se calculará la intensidad para cada modelo de potencia, impedancia o
intensidad constante.
Los resultados de estos modelos se compararán con los del flujo de cargas, donde se com-
probará que funcionan correctamente si los valores obtenidos son idénticos.

En último lugar se procederá al análisis temporal. En este caso los modelos se modifican
de manera que permitan estudiar el régimen transitorio del sistema.
En principio las simulaciones se harán con valores de potencia constantes. Más adelante se
obtendrán curvas de carga, que se utilizarán para estudiar y comparar el funcionamiento
del sistema. Para demostrar el correcto funcionamiento de los modelos, se variarán las
condiciones de trabajo y se comprobará si, en cada caso, se mantiene constante la variable
(potencia, impedancia o intensidad) necesaria.
Se trabajará con tres curvas de carga. Las curvas de carga representan el consumo de
potencia frente al tiempo. Las curvas empleadas de consumo industrial corresponden a
tres dı́as en un año separadas entre sı́ por varios meses. La localización será la misma para
las tres.
Estos perfiles se aplicarán a los tres modelos temporales, obteniendo un total de nueve
simulaciones además de las tres para valores de potencia constante. Se comprobará si la
constante de tiempo del sistema es lo suficientemente pequeña como para que entre cada
intervalo correspondiente a un valor de potencia, se llegue a alcanzar el régimen perma-
nente. Si se observan incongruencias y las gráficas representan oscilaciones, implicará que

Sabela González Castro III

 Estudio de modelos de carga en sistemas eléctricos de distribución

todavı́a se mantiene en el régimen transitorio.

Por interés se empleará un programa que permita generar y programar curvas de car-
ga de uso residencial, pero que no se aplicarán a los modelos.

Palabras clave: redes de distribución, modelos de carga, flujo de cargas, análisis fasorial,
análisis temporal, curvas de carga, potencia, intensidad, impedancia.

Código UNESCO: 330609 - Transmisión y distribución.

Sabela González Castro IV

Índice general

Resumen IV

1. Introducción 1
1.1. Sistemas de distribución. . . . . . . . . . . . . . . . . . . . . . . . . . . . . 1
1.1.1. Componentes de una red de distribución. . . . . . . . . . . . . . . . 2
1.1.2. Redes con cargas desequilibradas. . . . . . . . . . . . . . . . . . . . 3
1.2. Modelos de cargas. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 5
1.2.1. Modelos de cargas estáticas. . . . . . . . . . . . . . . . . . . . . . . 5
1.2.2. Modelos de cargas dinámicas. . . . . . . . . . . . . . . . . . . . . . 6
1.3. Modelos de curvas de carga. . . . . . . . . . . . . . . . . . . . . . . . . . . 7
1.3.1. Modelos Top-Down. . . . . . . . . . . . . . . . . . . . . . . . . . . 8
1.3.2. Modelos Bottom-Up. . . . . . . . . . . . . . . . . . . . . . . . . . . 8

2. Objetivos y Alcance del trabajo 10

3. Metodologı́a 12
3.1. Análisis en régimen estacionario . . . . . . . . . . . . . . . . . . . . . . . . 13
3.1.1. Flujo de cargas. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 13
3.1.2. Régimen fasorial. . . . . . . . . . . . . . . . . . . . . . . . . . . . . 15
3.2. Análisis temporal. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 29
3.2.1. Modelo potencia constante. . . . . . . . . . . . . . . . . . . . . . . 29
3.2.2. Modelo impedancia constante. . . . . . . . . . . . . . . . . . . . . . 37
3.2.3. Modelo intensidad constante. . . . . . . . . . . . . . . . . . . . . . 42
3.3. Generación y aplicación de curvas de carga . . . . . . . . . . . . . . . . . . 47

4. Resultados 50
4.1. Valores obtenidos y representaciones. . . . . . . . . . . . . . . . . . . . . . 50
4.1.1. Resultados del flujo de cargas. . . . . . . . . . . . . . . . . . . . . . 50
4.1.2. Resultados del modelo fasorial . . . . . . . . . . . . . . . . . . . . . 54
4.1.3. Resultados del modelo temporal. . . . . . . . . . . . . . . . . . . . 59

V
 Estudio de modelos de carga en sistemas eléctricos de distribución

4.1.4. Impedancia constante. . . . . . . . . . . . . . . . . . . . . . . . . . 67

4.1.5. Intensidad constante. . . . . . . . . . . . . . . . . . . . . . . . . . . 71
4.2. Esquema y comparación de resultados. . . . . . . . . . . . . . . . . . . . . 77
4.2.1. Flujo de cargas y análisis fasorial. . . . . . . . . . . . . . . . . . . . 77
4.2.2. Análisis temporal. . . . . . . . . . . . . . . . . . . . . . . . . . . . . 80

5. Conclusiones 84

6. Impacto Social y Ambiental. 85

7. Planificación Temporal y Presupuesto 87

7.1. Planificación temporal . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 87
7.2. Presupuesto . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 89

8. Trabajos futuros. 91

Trabajos futuros 91

Lista de figuras 94

Lista de tablas 97

Sabela González Castro VI

Capı́tulo 1

Introducción

1.1. Sistemas de distribución.

Las redes de distribución son una parte del sistema de suministro de energı́a eléctrica que
se encarga de transportar la energı́a desde la subestación hasta los receptores.
En la Figura 1.1 se presenta un esquema de los componentes de los sistemas de distribu-
ción.

Figura 1.1: Componentes de sistemas de energı́a eléctrica.

A la hora de diseñar las redes eléctricas es fundamental conocer la carga con la que se
va a funcionar y trabajarlas de acorde a ello. Ası́ se puede operar con mayor eficiencia y
menores pérdidas.
El sistema de distribución parte de la subestación de distribución. Esta está alimentada
por las lı́neas de subtransmisión, aunque a veces directamente de una lı́nea de tranmisión
de alta tensión.
Por lo general, los sistemas de distribución son radiales, lo que implica que solo hay un
sentido de recorrido de la potencia: desde la subestación de distribución hasta el consumo.
La carga de los sistemas de distribución es inevitablemente desequilibrada debido al gran
numero de cargas monofásicas que deben ser alimentadas.
Esto complica muchı́simo el análisis y el estudio frente a cargas equilibradas, donde se

Sabela González Castro 1

 Estudio de modelos de carga en sistemas eléctricos de distribución

analiza por fase. Por ello es necesario diseñar las lı́neas de la manera más precisa posible
empleando modelos.
No hay un diseño fijo de subestaciones, pero la gran mayorı́a están contectados a una o
más fuentes de alimentación radiales.

1.1.1. Componentes de una red de distribución.

Los componentes que forman la red de distribución son los siguientes:

Alimentación primaria trifásica

Laterales trifásicos, bifásicos y monofásicos

Reguladores de tensión (Estos solo se emplean en E.E.U.U, en España se utilizan

transformadores con tomas)

Transformadores de lı́nea

Bancos de condensadores

Transformadores de distribución

Cargas trifásicas, bifásicas y monofásicas

Estos elementos se representan en la Figura 1.2. [2]

Sabela González Castro 2

 Estudio de modelos de carga en sistemas eléctricos de distribución

Figura 1.2: Componentes de una red de distribución

1.1.2. Redes con cargas desequilibradas.

A la hora del estudio de circuitos, lo más básico es la resolución de los equilibrados. Esto
se presenta cuando todos los elementos del circuito en cada fase son simétricos y equili-
brados, presentando las mismas caracterı́sticas para todas las fases. Ası́, se puede estudiar
el circuito con el equivalente monofásico, mediante componentes de fase, simplificando
mucho su análisis.

El desequilibrio es debido a factores como la diferencia de impedancias, transformado-

res conectados en estrella o triángulo abiertos y sobrecarga en una de las fases, entre
muchos otros.

Los desequilibrios son perjudiciales para el sistema eléctrico, pues con que la tensión
varı́e un poco puede originar diferencias enormes en la intensidad. Esto es indeseable tan-
to porque reduce la seguridad y estibilidad como que perjudica los aislamientos.

Sabela González Castro 3

 Estudio de modelos de carga en sistemas eléctricos de distribución

Para poder resolver también los circuitos con cargas desequilibradas, Fortescue desarrolló
el método de las componentes simétricas. Este permite descomponer un sistema trifásico
asimétrico y desacoplado en la suma de tres sistemas simétricos y resolubles.
Esto es aplicable a cualquier sistema polifásico y no únicamente trifásico, siempre y cu-
nado esté desacoplado.
Los sistemas en los que se descompone el circuito desequilibrado son:

Sistema de secuencia 1, positiva o directa

Sistema, antagónico al primero, de secuencia 2, negativa o inversa

Sistema, de tres fasores monofásicos de igual módulo, de secuencia 0, homopolar o

nula.

Los sistemas citados aparecen representados en la Figura 1.3.

Figura 1.3: Descomposición en componentes simétricas.

Las redes eléctricas son en sı́ desequilibradas y las fases están acopladas. Lo que implica
que en este caso no supondrı́a ninguna ventaja el descomponer el sistema en tres subsiste-
mas resolubles de manera intependiente. Es decir, no ayudarı́a el resolverlo por el método
de las componentes simétricas mencionado.
Por este motivo se resuelve cada una de las fases de la red por separado, mediante com-
ponentes de fase, sin descomponerlas según los sistemas mostrados.

Sabela González Castro 4

 Estudio de modelos de carga en sistemas eléctricos de distribución

1.2. Modelos de cargas.

Hasta ahora se ha expuesto la importancia del estudio de las redes y su diseño. A conti-
nuación se proceden a explicar tipos de cargas que pueden ser alimentadas por la red.

La estabilidad de un sistema de potencia viene determinada por la relación entre la energı́a

ofertada y la demanda de la carga. Por ello, las cargas juegan un papel crucial a la hora de
analizar y estudiar un sistema. Esto es tan complicado como importante, ya que depende
de la hora del dı́a, estación, equipo receptor y tiempo de uso, entre muchas otras.

Los modelos de carga se pueden dividir en estáticos y dinámicos. En este proyecto se

trabajará con modelos estáticos, pero antes se explicará cada uno contrastando las dife-
rencias.

1.2.1. Modelos de cargas estáticas.

Estos modelos representan las caracterı́sticas de la carga como funciones algebraicas de la
magnitud y frecuencia del nudo de tensión en cada instante. Para ello, se opera indepen-
dientemente con la potencia activa y reactiva.

El primer modelo es el que se ha usado tradicionalmente: el Modelo exponencial.

Viene dado por las siguients expresiones [3]:

P = P0 (V̄ )a (1.1)

Q = Q0 (V̄ )b (1.2)

V
V̄ = (1.3)
V0

Como parámetros del Modelo exponencial están a y b. Estos pueden tomar los valores 0,
1 o 2 de manera que representan potencia constante, corriente constante e impedancia
constante, respectivamente.

Para cargas compuestas los valores dependen de las caracterı́sticas de cada parte. Es
decir, los parámetros tomarán valores entre 0, 1 y 2 obtenidos según las caracterı́sticas
de la carga, al no operar simplemente con potencia, impedancia o intensidad constante.
Por ejemplo: si la carga trabaja en mayor medida consumiendo una potencia constante

Sabela González Castro 5

 Estudio de modelos de carga en sistemas eléctricos de distribución

frente al resto de caracterı́sticas, el valor que tomarán a o b será más próximo a 0 (pero
no igual) que al resto, ya que corresponde a potencia constante.

Este Proyecto se centrará en los valores simples y se trabajará con los tres modelos citados.

Otro modelo alternativo es el que representa la dependencia de la tensión y las cargas: el

Modelo polinómico, también llamado ZIP pues está compuesto de impedancia constante
(Z), intensidad constante (I) y potencia constante (P).

P = P0 [p1 V̄ 2 + p2 V̄ + p3 ] (1.4)

Los parámetros son los coeficientes pi , los cuales definen la proporción de cada componente.

La dependencia de la frecuencia de la carga es normalmente representada por un fac-

tor que multiplica a las expresiones exponencial o polinómica.

Otra manera es combinar ambos, teniendo:

P = P0 [PZIP + PEX1 + PEX2 ]. (1.5)

Donde:
PZIP = p1 V̄ 2 + p2 V̄ + p3 (1.6)

PEX1 = p4 (V̄ )a1 (1 + Kpf 1 ∆f ) (1.7)

PEX2 = p5 (V̄ )a2 (1 + Kpf 2 ∆f ) (1.8)

Con Kpf 1 un factor relacionado con la frecuencia.

Siendo las expresiones similares para la potencia reactiva.

1.2.2. Modelos de cargas dinámicas.

Para cambios moderados de voltaje y frecuencia, los modelos estáticos descritos con an-
terioridad son adecuados. Sin embargo, cuando la respuesta del sistema no se realiza
lo suficientemente rápido se requieren modelos dinámicos de cargas. Este es el caso del
estudio de oscilaciones especı́ficas, modelado de motores de inducción y estabilidad de
tensiones, entre otros [4].

Sabela González Castro 6

 Estudio de modelos de carga en sistemas eléctricos de distribución

Además se pueden estudiar modelos en los que las cargas sean compuestas. Sin embargo,
esto no está dentro de los lı́mites de este trabajo y, por ello, no se profundizará más.

1.3. Modelos de curvas de carga.

Las redes de distribución transportan la energı́a eléctrica para que esta pueda ser consu-
mida. Para una mayor eficiencia es indispensable conocer el comportamiento de la carga,
que no es más que el consumo de esta energı́a.
Si se conoce con anterioridad la demanda de potencia en cada instante o intervalo de
tiempo, la red puede configurarse de manera trabaje con menores pérdidas y proporcione
lo necesario. Sin embargo, a la hora de estudiar el consumo, es inmediato pensar que
intervienen millones de parámetros variables en cada vivienda o centro receptor.

Las curvas de carga representan el consumo de potencia respecto del tiempo, de ma-
nera que da información de la demanda. Por ello, se han ido realizando estudios a lo largo
de los últimos años con el fin de obtener un algoritmo o metodologı́a que permita obtener
unos modelos de trabajo eficiente. Estos modelos han de ser paramétricos, ya que deben
permitir simular diferentes situaciones.

A continuación se presentan los Modelos de curvas de carga, divididos en dos grandes

bloques: los modelos “Top-down” y los modelos “Bottom-up”.

Figura 1.4: Esquema de representación de los modelos de curvas de carga.

Sabela González Castro 7

 Estudio de modelos de carga en sistemas eléctricos de distribución

1.3.1. Modelos Top-Down.

Se atribuye un consumo de energı́a a cada vivienda de acuerdo a las caracterı́sticas que
esta presenta. Son modelos estadı́sticos mediante desagregación.
Se toma como dato la curva de carga total de una residencia (incluyendo todo fin último)
y se propone separar de esta cada contribución. Es decir, se miden curvas de carga y se
dividen para identificar aparatos y equipos.

Dentro de estos modelos destaca el modelo de Aigner et al [5]. Los autores partieron
de un desarrollo matemático llamado Análisis de demanda condicional (CDA) y lo
mejoraron.
Como datos de entrada toman curvas de carga medidas en intervalos de 15 minutos.
Las variables que presentan la demanda de potencia son: la temperatura deseada del
hogar, el tamaño de la vivienda y otras que resumen la presencia o no de equipos
receptores.
El modelo propuesto es adecuado pero algo limitado en cuanto a la precisión de los
resultados. Puede ser mejorado dependiendo de la información de la que parte como
entrada, es decir, los datos sobre el uso real de los electrodomésticos.

El segundo y último que se comenta en este proyecto es el de Bartels et al [6]. Este

requiere más información como entrada y consiste en construir una curva de carga
DELMOD desarrollada a partir del modelo CDA.
El objetivo de los autores era obtener una herramienta que fuese capaz de simular
el impacto de diferentes escenarios sobre la demanda regional de potencia.
DELMOD está compuesta por dos módulos. El primero es DELMOD B (base): cal-
cula en un dı́a especı́fico de trabajo la curva del sector residencial.
El segundo es DELMOD W (clima): estudia cómo el tiempo atmosférico influye en
los resultados de DELMOD B. Como resultados de este modelo se tienen curvas de
carga horarias.

1.3.2. Modelos Bottom-Up.

Se calcula el consumo de energı́a de una vivienda o conjunto de viviendas y se extrapola
para todas las residencias, a partir de datos microscópicos.
Los datos de entrada son el consumo individual de electrodomésticos, sus propiedades
o caracterı́sticas técnicas, las propiedades geométricas y térmicas de las residencias a
estudiar, el clima, las facturas de electricidad de los residentes y el comportamiento y
acciones de estos.
Ası́ como el modelo anterior era estadı́stico mediante desagregación, este puede dividir en

Sabela González Castro 8

 Estudio de modelos de carga en sistemas eléctricos de distribución

modelo estadı́stico aleatorio, modelo empı́rico probabilı́stico y modelo basado en el uso

del tiempo.

El Modelo aleatorio estadı́stico fue presentado por Yao y Steemers [7], llamado Méto-
do Simple de fomulación de Curvas de Carga (SMLP). Con él se elaboran Curvas de
Carga con la consideración de varios escenarios de personas. Para realizarlas toman
la peor de las situaciones en cada residencia, con intervlalos de 1, 5, 15 o 30 minutos.

En Modelo empı́rico probabilı́stico destacan Stokes y Paatero y Lund. Sin embargo

se comentará únicamente y brevemente el modelo de Stokes, sin entrar en cálculos y
teorı́as matemáticas. Stokes [8] desarrolló un modelo aplicable a un sector doméstico
determinado capaz de generar tres tipos de curvas de carga globales: el primer nivel
o 30 minutos de demanda para un residente de la media, el segundo nivel o 30
minutos de demanda para una vivienda especı́fica y el tercer nivel o 1 minuto de
demanda para una vivienda seleccionada de un cliente. Calcula el factor de potencia
de cada elemento en cada instante.

Los Modelos basados en el uso del tiempo se corresponden a unas curvas de uso de
frecuencia diaria para cada electrodoméstico/equipo.

Sabela González Castro 9

Capı́tulo 2

Objetivos y Alcance del trabajo

Si se reflexiona sobre el dı́a a dı́a, es inmediato pensar que las cargas alimentadas por los
sistemas de distribución varı́an. No todos los interruptores, por ejemplo, de cada vivienda
se encienden y se apagan al mismo tiempo. Lo mismo ocurre con dispositivos electrónicos
y electrodomésticos.
Por ello, es imposible que las redes de distribución trabajen con cargas equilibradas.

Como se puede ver es realmente importante estudiar y analizar los sistemas desequili-
brados de manera que permita aumentar el rendimiento de los sistemas de distribución
y mejorar considerablemente su trabajo. Actualmente se siguen buscando soluciones que
permitan aumentar la capacidad del sistema.

Las redes inteligentes o Smart Grids son aquellas que pueden integrar de forma eficiente el
comportamiento y las acciones de todos los usuarios que están conectados a ellas. De esta
manera, implican un sistema con menores pérdidas y mayores seguridad y capacidad. [1]
Para poder optimizar la producción y transporte de electricidad requiere de tecnologı́a
eléctrica y medidores inteligentes. Es decir, para el eficiente funcionamiento de la red se
tiene en cuenta a la demanda.

Sabela González Castro 10

 Estudio de modelos de carga en sistemas eléctricos de distribución

Figura 2.1: Relación de Smart Grid.

Para poder estudiar una red desequilibrada y ”modificarla”de manera que opere con ge-
neración distribuida como una Smart Grid, hay que realizar una serie de tareas que son
el objeto de estudio de este Proyecto de Fin de Carrera.

Como es lógico, son muchos los parámetros y elementos que influyen en los sistemas de
distribución que hay que tener en cuenta a la hora de diseñarlos y que se pueden abordar
como objeto de estudio e investigación.
Este trabajo se centra en los modelos de carga con el fin de poder analizar el comporta-
miento de una lı́nea para cada uno de ellos. Ası́, se realiza el estudio con los tres modelos
estáticos principales: potencia constante, impedancia constante e intensidad constante;
todos ellos elaborados y trabajados en Simulink de Matlab.

Los tres modelos citados se desarrollan en las dos grandes partes en las que se divide
este Proyecto de Fin de Grado. En la primera mitad, se aborda el estudio en régimen
estacionario.

En la segunda mitad, se aborda el estudio en régimen temporal. Ası́ como en la primera

parte se opera con fasores que representan el punto de trabajo en régimen permanente,
en este caso se tiene en cuenta el transitorio y los valores que toma el sistema en cada
instante. Se realizarán las comparaciones pertinentes con el régimen estacionario.

Sabela González Castro 11

Capı́tulo 3

Metodologı́a

Para realizar los modelos descritos con anterioridad se ha empleado la herramienta Simu-
link de Matlab.
Lo primero que se ha desarrollado es un análisis de flujo de cargas para un modelo de
cuatro nudos, calculando previemente los parámetros de la lı́nea (las impedancias de la
lı́nea a ambos lados del transformador).
La herramienta del flujo de cargas permite conocer los valores medidos en la lı́nea para
un determinado instante. Ası́, presenta una opción que permite seleccionar una carga que
trabaje a potencia, impedancia o intensidad constante. Esto se selecciona mediante el
bloque powergui : continuous.

A continuación se han realizado los modelos en régimen estacionario operando con fa-
sores y comparando con el caso anterior.
Gracias al análisis fasorial se estudian varios puntos a lo largo del tiempo. Es decir, serı́a
como hacer infinitos flujos de carga.
Para esta opción, Simulink no presenta la herramienta anterior de una carga que te permi-
ta seleccionar el tipo de funcionamiento. Es por esto que en esta parte se realizan modelos
que simulen cada una de las tres cargas.

El último paso ha sido modificar los modelos estacionarios para poder estudiar el transito-
rio y el comportamiento temporal. En este caso ya no se trabaja en un punto determinado,
sino que las variaciones han comprendido valores instantáneos. En este caso se selecciona
en el bloque powergui: phasor 60 Hz.

Para finalizar, se ha usado el programa Load Profile Generator. Gracias a él se han ob-
tenido curvas de carga muy particulares, ya que permite elegir y combinar viviendas con
gran cantidad de variables para proporcionar curvas de carga que represente un consumo

Sabela González Castro 12

 Estudio de modelos de carga en sistemas eléctricos de distribución

de potencia adecuado.

3.1. Análisis en régimen estacionario

Para comenzar se trabaja en dominio fasorial, con una resolución permamente a lo largo
del tiempo.

El primer análisis se hace para un flujo de cargas. Es decir, sobre un punto determinado en
el tiempo. Esto es sencillo puesto que Simulink permite obtener los cálculos rápidamente.
El segundo análisis es el fasorial. En este caso se estudia la solución en regimen perma-
nente a lo largo del tiempo. Es decir, serı́a equiparable a hacer infinitos flujos de carga.
Este caso es más complicado, ya que habrá que diseñar los modelos.

Se empieza realizando una serie de modelos brevemente mencionados al principio y, con

menos complejidad que en la siguiente parte, se aborda el estudio del comportamiento de
la lı́nea en régimen estacionario.
De esta manera se encuentra limitado a trabajar con valores correspondientes a un instante
de tiempo y poder operar mediante fasores que representen a cada una de las variables
eléctricas: tensión, intensidad y potencia. El modelo general, que es sobre el cual se trabaja
en este Proyecto, se presenta más adelante. Representa una lı́nea de cuatro nudos con un
transformador y una fuente de tensión de alterna.

3.1.1. Flujo de cargas.

Se implementa en Simulink una lı́nea de cuatro nudos con un transformador. Los datos
se calculan, a partir de un enunciado obtenido del IEEE [9], en Matlab y se exportan a
Simulink.
La lı́nea es desequilibrada, lo que se comprueba si se observan las matrices de los paráme-
tros de la lı́nea: las componentes fuera de la diagonal no son simétricas.

Cada bloque naranja representa un nudo (el circuito presenta cuatro en total) y, al ejecutar
como flujo de cargas, cada uno debe programarse como nudo balance, carga o generador.

Sabela González Castro 13

 Estudio de modelos de carga en sistemas eléctricos de distribución

Figura 3.1: Modelo de 4 nudos para el flujo de cargas.

El transformador presente es un Dyg11, es decir con la estrella conectada a tierra. La

frecuencia de trabajo son 60 Hz y las tensiones nominales son 12470 V para el primario y
4160 V para el secundario (ambos valores de lı́nea).
Los bloques LY y LD representan las impedancias de la lı́nea calculadas.

Para poder establecer una comparación con los modelos siguientes, se emplea la opción
de Simulink de carga modificable para flujos de cargas. Se puede escoger que la carga
trifásica consuma potencia constante, impedancia constante o intensidad constante.
Se selecciona en el bloque de la carga de la siguiente manera:
Three-Phase Series RLC Load > Load Flow > Load type: constant Z, constant PQ y
constant I.

Figura 3.2: Desplegable de la carga.

Es decir, para cada tipo de carga del modelo de flujo de cargas se obtienen unos valores
que se compararán más adelante con los resultados obtenidos en la siguiente sección.
Para ejecutar el flujo de cargas se acciona de la manera siguiente:

Sabela González Castro 14

 Estudio de modelos de carga en sistemas eléctricos de distribución

Figura 3.3: Obtención de valores del flujo de cargas.

3.1.2. Régimen fasorial.

Para el análisis fasorial no se encuentra la opción de la carga que se escoge para que
funcione a potencia, impedancia o intensidad constantes. Es por esto que se desarrollan
los modelos. Para todos los modelos realizados en el proyecto se cuenta con una carga
trifásica con fases a,b y c conectada en estrella.

Figura 3.4: Carga conectada en estrella.

Sabela González Castro 15

 Estudio de modelos de carga en sistemas eléctricos de distribución

El sistema trabaja realizando iteraciones hasta que converge a un resultado que cumple
las condiciones de cada modelo: potencia, impedancia o corriente constantes. A partir de
esto se desarrolla cada caso.

Antes de presentar la construcción de cada modelo en Simulink, se expone la base teórica

de cada uno.

Para la carga conectada en estrella, se cumple en los tres modelos que las potencias
de fase son:
Sa = |Sa | θa (3.1)

Sb = |Sb | θb (3.2)

Sc = |Sc | θc (3.3)
En el modelo de potencia constante se tiene en cuenta lo siguiente:
Sa ∗
 
|Sa |
ILa = = δa − θa = |ILa | αa (3.4)
Van |Van |
 ∗
Sb |Sb |
ILb = = δb − θb = |ILb | αb (3.5)
Vbn |Vbn |
 ∗
Sc |Sc |
ILc = = δc − θc = |ILc | αc (3.6)
Vcn |Vcn |
Donde la potencia consumida será constante si se varı́an las condiciones nominales
o no, mientras que la tensión y la intensidad varı́an.

En el modelo de impedancia constante se tiene en cuenta que:

|Van |2 |Van |2
Za = = θa = |Za | θa (3.7)
Sa ∗ |Sa |

|Vbn |2 |Vbn |2
Zb = = θb = |Zb | θb (3.8)
Sb ∗ |Sb |
|Vcn |2 |Vcn |2
Zc = = θc = |Zc | θc (3.9)
Sc ∗ |Sc |

Sabela González Castro 16

 Estudio de modelos de carga en sistemas eléctricos de distribución

A partir de ello se calculan las intensidades de lı́nea según las expresiones:

 
Van
ILa = (3.10)
Za
 
Vbn
ILb = (3.11)
Zb
 
Vcn
ILc = (3.12)
Zc
Todo esto aparece representado en los siguientes apartados, donde se explican los
modelos en Simulink.
En este caso mientras que al modificar las condiciones nominales la tensión consu-
mida varı́a, la impedancia se mantiene constante.

En el modelo de intensidad constante se cumple:

ILa = |ILa | δa − θa (3.13)

ILb = |ILb | δb − θb (3.14)

ILc = |ILc | δc − θc (3.15)

Donde δi representa los ángulos de las tensiones de fase y θi los ángulos de la po-
tencia, para i= a, b y c.
En este caso se mantendrán constantes el módulo de la intensidad y la componente
de la fase correspondiente al factor de potencia (θ).

[Link]. Modelo potencia constante.

El primer modelo de análisis fasorial desarrollado en Simulink es el de potencia constante.

Al implementarlo se tuvo en cuenta lo mostrado en la sección anterior para las ecua-

ciones de potencia (3.1) a (3.3) y (3.4) a (3.6).

A continuación, se muestran imágenes del modelo para la mejor comprensión de la expli-

cación.

Sabela González Castro 17

Figura 3.5: Modelo fasorial de potencia constante.
 Estudio de modelos de carga en sistemas eléctricos de distribución

A primera vista se observa el mismo circuito eléctrico que en el modelo del flujo de cargas
expuesto con anterioridad. Antes, tal y como se explicó al principio, se trabajaba con
powergui : continuous y ahora con P hasor 60.

Como se explicó al comienzo del capı́tulo, para la opción de análisis fasorial no existe
la opción de una carga que trabaje a potencia, impedancia o intensidad constantes.

La carga a configurar en Simulink se representa con tres fuentes de intensidad y una

resistencia en paralelo. La función de esta es evitar posibles indeterminaciones del progra-
ma.
De esta manera se puede trabajar con esta carga con potencia PQ constante, ya que la
intensidad que “define” la carga es la correspondiente a un cálculo en el que la potencia
no varı́a.

Se observan dos bloques a la derecha del nudo 4. La función que tienen es de medidores
de tensión e intensidad. Para ello, se les asigna dos “tags” V abc, y V abc4 respectivamente,
de manera que se pueden estudiar los valores más arriba, tal y como se ve.

Se comprueba que en ambos puntos se obtienen los mismos valores. Los destacados en
azul son los valores obtenidos de Vabc y los rosas de V abc4.

Sabela González Castro 19

Figura 3.6: Comparación de valores obtenidos en dos puntos.
 Estudio de modelos de carga en sistemas eléctricos de distribución

Los módulos de la tensión son valores eficaces y los ángulos están en grados.

El objetivo era calcular para cada modelo una intensidad que, efectivamente, implica-
se un funcionamiento a potencia PQ constante. Como se puede observar, la corriente que
se obtiene como salida del bloque Subsystem P Q cte se inyecta en las tres fuentes de
intensidad.

Figura 3.7: Ampliación de la parte modificada.

Para posteriores simulaciones, se da la opción de introducir los valores de trabajo deseados

de potencia activa y el factor de potencia. En este caso, la potencia de trabajo es de 1800
kW y el factor de potencia 0,9. Observando los medidores de la etiqueta Vabc4 se comprueba
que, efectivamente, la potencia medida es la que se ha impuesto.

Sabela González Castro 21

 Estudio de modelos de carga en sistemas eléctricos de distribución

Figura 3.8: Máscara del Subsystem P Q cte.

Adentrándonos en el subsistema:

Figura 3.9: Interior del Subsystem P Q cte.

El bloque M emory1 aplica un pequeño retraso para que no haya indeterminaciones en el

sistema. √
La ganancia K tiene valor - 2. El signo negativo es por el sentido de la intensidad, y el
valor es porque se transforma de valor eficaz a amplitud de la intensidad. (La ganancia K

Sabela González Castro 22

 Estudio de modelos de carga en sistemas eléctricos de distribución

a la derecha de V abc4 transforma a la tensión en valor eficaz.

Como se puede observar, a partir de las P y Q fijadas y la tensión medida en el nudo 4 (la
cual se pasa a valor eficaz), se obtiene el valor de la intensidad mediante los cálculos que
representan el modelo y que se han expuesto al inicio de esta sección (ecuaciones (3.4),
(3.5) y (3.6)).

[Link]. Modelo impedancia constante.

El siguiente modelo a desarrollar es el de impedancia constante. Como su propio nombre
indica, en este modelo se trabaja de manera que la lı́nea opera a impedancia constante.
Para ello se aplica una metodologı́a similar a la del apartado anterior. Partiendo de las
potencias de fase indicadas en las ecuaciones (3.1) a (3.3) se obtienen las impedancias me-
diante las ecuaciones (3.7) a (3.9) indicadas en el inicio de este capı́tulo. Las intensidades
de lı́nea se calculan a partir de la impedancia constante gracias a las ecuaciones (3.10) a
(3.12).

En este modelo las tensiones fase-neutro varı́an en cada iteración, pero las impedancias
expuestas anteriormente se mantendrán constantes.

A continuación, se muestran imágenes del modelo que se explican para su mejor com-
prensión.

Sabela González Castro 23

Figura 3.10: Modelo fasorial de impedancia constante.
 Estudio de modelos de carga en sistemas eléctricos de distribución

La distinción de este modelo del anterior se basa en el interior del bloque Subsystem Z cte.

Exterior a este bloque el modelo es exactamente igual y funciona de la misma forma:

se obtiene como salida del bloque una intensidad caracterı́stica del funcionamiento a im-
pedancia constante la cual se inyecta en la carga simulada por tres fuentes de intensidad,
como ya se explicó en el modelo anterior.

Para posteriores simulaciones, se da la opción de introducir los valores de trabajo deseados

de tensión nominal de fase, potencia activa y factor de potencia.

Figura 3.11: Máscara del Subsystem Z cte.

Como el exterior al subsistema no varı́a en nada más respecto del caso anterior, se procede
a comentar directamente el Subsystem Z cte.

Entrando en el modelo se tiene lo siguiente, que viene representado en la Figura 3.12:

a partir de las P y Q nominales se obtiene la conjugada de la potencia aparente.
Con la relación expuesta al principio de esta sección, se utiliza la tensión nominal para
calcular la impedancia.
El tag F rom2 mide la tensión en el nudo 4 y calcula la intensidad a inyectar a partir de
la impedancia anterior (constante).
El bloque M emory1 aplica un pequeño retraso para que no haya indeterminaciones en el
sistema.

Sabela González Castro 25

 Estudio de modelos de carga en sistemas eléctricos de distribución

Figura 3.12: Subsystem Z cte.

Como se puede observar al hacer la simulación, se obtienen valores de impedancia positi-

vos, lo que es lógico.

[Link]. Modelo intensidad constante.

El objetivo ahora es simular el comportamiento de la lı́nea trabajando a intensidad cons-
tante.
Para ello se ha de tener en cuenta lo mencionado a partir de la carga en estrella y las
potencias de fase indicadas al inicio del capı́tulo (ecuaciones (3.1) a (3.3).
Las intensidades a calcular vienen representadas en las ecuaciones (3.13) a (3.15). Como
es lógico, el módulo de la intensidad se mantendrá constante y la componente de la fase
debido al factor de potencia también.

De nuevo la base es la misma, variando solo el cálculo de la intensidad resumido en

el bloque Subsystem I cte.
El modelo se representa en la Figura 3.13.

Sabela González Castro 26

Figura 3.13: Modelo fasorial de intensidad constante.
 Estudio de modelos de carga en sistemas eléctricos de distribución

Con el mismo objetivo que en los dos modelos anteriores, se configura una máscara de
manera que da la opción de simular distintas situaciones variando la potencia activa, el
factor de potencia y la tensión nominal.

Figura 3.14: Máscara del Subsystem I cte.

El subsistema presenta de nuevo el cálculo de la intensidad. En este caso juegan un papel

muy importante las fases y, por ello, se procederá a su explicación más detenidamente.

Figura 3.15: Subsystem I cte.

Como se puede ver, se parte de tres constantes: potencia activa P, potencia reactiva Q y
tensión nominal fase-neutro; asignadas previamente en la máscara externa.

Por un lado, con el ángulo de la potencia aparente y la fase de la tensión medida en

Sabela González Castro 28

 Estudio de modelos de carga en sistemas eléctricos de distribución

el nudo 4 se obtiene, mediante la resta, el ángulo de la intensidad buscada.

Por otro lado, con el módulo de la potencia aparente y el de la tensión nominal se obtiene,
mediante la división, el módulo de la intensidad a inyectar. La ganancia negativa se añade
ya que la intensidad circula en sentido contrario.

En este modelo varı́a el ángulo de la corriente que, desglosado, lo forman el ángulo de la

tensión medida (variable) y el ángulo del factor de potencia (constante).

3.2. Análisis temporal.

En la primera parte se ha estudiado el comportamiento de una lı́nea para análisis en régi-
men estacionario, donde la resolución es permanente a lo largo del tiempo.
En esta segunda parte, se considerará una solución en cada instante del tiempo. Para ello,
se parte de los modelos con los que se ha trabajado previamente y se modifican de manera
que permitan el análisis temporal.

Los tres modelos son trifásicos donde cada fase está conectada a la entrada de una fuente.
Para su mejor comprensión se explica estudiando una fase, al ser las otras dos restantes
idénticas en la construcción, que no en los resultados.

3.2.1. Modelo potencia constante.

La base teórica de este modelo es la misma que la expuesta en la sección 3.1.2. Sin
embargo, al tener en cuenta el tiempo para cada instante hay que realizar una serie de
modificaciones.
Para su mejor comprensión se explicará con ayuda del modelo realizado. En la Figura
3.16 se muestra el modelo temporal de potencia constante.

Sabela González Castro 29

Figura 3.16: Modelo temporal de potencia constante.
 Estudio de modelos de carga en sistemas eléctricos de distribución

Se observa a la izquierda de todo una fuente trifásica de corriente alterna conectada al

subsistema Subsystem P Q.
La salida de este subsistema la componen las tensiones e intensidades de cada fase. A
continuación se transforman en potencia activa y reactiva de cada fase, que aparecerán
representadas y medidas junto con las tensiones e intensidades a la derecha.

Para futuras simulaciones se crean dos máscaras en la fuente y en el Subsystem P Q.

Figura 3.17: Máscaras de fuente y Subsystem P Q.

La máscara de la fuente permite ajustar la tensión de fase en valor eficaz de la fuente y

la frecuencia de trabajo.
La máscara del subsistema Subsystem P Q permite ajustar la tensión de fase en valor
eficaz de la fuente, la frecuencia de trabajo, el ángulo de fase inicial de la tensión y el
valor inicial de la intensidad en por unidad.
Entrando en el Subsystem P Q, se tiene lo mostrado en la siguiente figura. El propósito
de esta es mostrar las tres fases.

Sabela González Castro 31

Figura 3.18: Máscaras de fuente y Subsystem P Q.
 Estudio de modelos de carga en sistemas eléctricos de distribución

Como ya se indicó previamente, se usará una fase para su explicación, al ser la construcción
de las dos restantes idéntica.

Figura 3.19: Fase A del Subsystem P Q.

Debe tenerse en cuenta lo siguiente. En todos los modelos se parte de unos valores de
potencia. Estos pueden ser constantes o valores correspondientes a curvas de carga.
En la Figura 3.20 se aprecian los bloques rosas P 4 y Q4 para operar con constantes. Los
bloques blancos que sı́ están conectados son bloques F rom W orkspace gracias a los cuales
se importan a Simulink los valores de curvas de carga desde Matlab.

Figura 3.20: Bloques P,Q constantes (rosas) y From Workspace.

La selección de las curvas de carga con las que se opera se explica más detallademente al
final del capı́tulo.

Sabela González Castro 33

 Estudio de modelos de carga en sistemas eléctricos de distribución

Se recibe como entrada la tensión fase-neutro de la fuente. Como ahora hay que tener
en cuenta el valor instantáneo y no se opera con fasores, es necesario calcular en cada
instante el valor eficaz y la componente wt.
Teniendo estos valores, podemos llegar a conocer el valor instantáneo de intensidad:

i(t) = Iamplitud cos(wt + φ) (3.16)

La Figura 3.21 es una ampliación del modelo anterior.

Figura 3.21: Ampliación de la Figura 3.19.

Para poder calcular el valor eficaz, se parte de un valor incial de tensión de fase (bloque
V f ). Este valor pasa por el RMS, que proporciona el verdadero valor eficaz de la tensión
al no poder, simplemente, dividir entre raı́z de dos. Esto se debe a que al trabajar ins-
tantáneamente, se calcula el valor eficaz mediante un muestreo de resultados anteriores.
De esta manera obtenemos para cada instante, el valor eficaz de la tensión.

A continuación: los ángulos. Para calcular el ángulo de fase de la tensión en cada ins-
tante, se divide al valor inicial de tensión V f entre el valor nominal, obteniendo el valor
de la tensiń en por unidad. Con este valor se entra en el bloque P LL (se explicará con
detalle en el siguiente párrafo).
Como primera salida del PLL se obtiene el ángulo α en radianes. De esta manera se ha

Sabela González Castro 34

 Estudio de modelos de carga en sistemas eléctricos de distribución

calculado el fasor de tensión para cada instante con el módulo en valor eficaz y el ángulo
de fase en radianes.

El bloque PLL (Phase-locked loop, o Lazo de seguimiento de fase) no es más que un

circuito que permite controlar la frecuencia y la fase a partir de una señal de referencia
externa. [10].
En la siguiente figura se muestra el interior del subsistema PLL:

Figura 3.22: Circuito PLL.

Con el valor de tensión de entrada en por unidad, el bloque superior obtiene el ángulo de
π
fase en grados y con la ganancia de valor se transforma en radianes.
180
A partir del valor de frecuencia de entrada en por unidad, multiplicando por la ganancia
de valor f rec (frec es el valor de frecuencia nominal definida en la máscara del Subsystem
P Q externo) y por 2π se tiene un valor que, integrado, proporciona el valor de wt.

Explicado esto, se pasa a la segunda parte del diagrama, representado en la Figura 3.23.

Sabela González Castro 35

Figura 3.23: Ampliación de la Figura 3.19.
 Estudio de modelos de carga en sistemas eléctricos de distribución

Se parte del fasor de tensión y el valor de wt. Tal y como viene expresado en las ecuaciones
(3.4), (3.5) y (3.6):  ∗
S |S|
IL = = δ − θ = |IL | α (3.17)
Vn |Vn |
Tras dividir la potencia compleja (conector sobre el que está S = P + jQ ) entre la tensión
y hacer su conjugada (bloque V̄ ) se tiene la intensidad. Los bloque real e imag son un
retardo para evitar indeterminaciones en el caso de que la tensión tome valor nulo en un
instante determinado.

A continuación, el valor complejo de intensidad obtenido se divide en módulo y ángu-

lo, con el proposito de añadir a la fase el valor de wt.
El módulo de I se multiplica por la ganancia raı́z de dos, para tener la amplitud de la
intensidad correspondiente a la ecuación (3.16). En este caso y al haber separado la com-
ponente wt se pueden relacionar ası́ valor eficaz y amplitud. Al ángulo de la intensidad se
le suma el wt y se crea la función seno que, multiplicada por la amplitud de la tensión,
porporcionaria la fución de la ecuación 3.16.
La intensidad final se multiplica por una ganancia negativa debido al sentido real de la
corriente, análogo a lo que ocurrı́a en los modelos estáticos (sección 3.1)

Como ya se mencionó antes, en las primeras simulaciones se usaron valores de poten-

cia constantes. Sin embargo lo interesante es analizar el comportamiento de la lı́nea en
situaciones lo más próximas a la realidad.
De esta manera, se modificó el modelo de manera que emplease como datos de entrada
perfiles de curvas de carga. Parte de ellas fueron creadas gracias al Load Profile Generator,
del que se hablará al final de este capı́tulo. También se utilizaron curvas de carga de uso
industrial.
Para poder trabajar con estos datos desde Simulink, se modificó el bloque constante por
un F rom W orkspace, gracias al cual se importan desde Matlab.

3.2.2. Modelo impedancia constante.

De nuevo, la base teórica de este modelo es la del de impedancia constante en régimen
fasorial (sección 3.1.3), con la diferencia de la influencia temporal.
En la Figura 3.24 viene representado el modelo desarrollado en Simulink.

Sabela González Castro 37

Figura 3.24: Modelo temporal de impedancia constante.
 Estudio de modelos de carga en sistemas eléctricos de distribución

En primer lugar se observa la misma fuente de tensión alterna trifásica que en el modelo
anterior. Los valores de la tensión son la entrada al subsistema Subsystem Z, que tiene
por salida los valores en módulo y complejos de la impedancia en cada fase.
Los valores de la impedancia vienen representados en gráficos (naranjas para el módulo,
azules para valor complejo) y por medidores que van indicando su valor en cada instante.

Igual que para el caso anterior, se hacen máscaras tanto en la fuente como en el Subsystem
Z para poder modificar los valores de entrada y de trabajo.

Figura 3.25: Máscaras de fuente y Subsystem Z.

Los valores a introducir en la fuente son la tensión de la fuente de fase en valor eficaz y
la frecuencia de trabajo.
Las variables a añadir para el subsistema son la tensión de la fuente de fase en valor eficaz,
la frecuencia de trabajo, el desfase inicial de la tensión y el valor inicial de la intensidad
en p.u.

Entrando en el Subsystem Z se tiene lo siguiente:

Sabela González Castro 39

Figura 3.26: Subsystem Z.
 Estudio de modelos de carga en sistemas eléctricos de distribución

Se observan las tres fases que de construcción son idénticas aunque operen con valores
ligeramente distintos (desfases).
Para su mejor comprensión se estudia a continuación la primera de las fases, representada
en la Figura 3.27.

Figura 3.27: Ampliación del Subsystem Z.

Como se puede ver, el modelo es similar al expuesto de potencia constante en la sección

3.2.1. A continuación se van a destacar y razonar las diferencias estructurales, dejando la
comparación de los valores obtenidos para los capı́tulos 4 y 5.

La obtención del valor eficaz, el ángulo de fase y la componente wt de la tensión de

entrada es análoga a la sección anterior. La diferencia está en la forma de calcular la
intensidad.
Ası́ como antes, sin considerar wt, se obtenı́a dividiendo la potencia aparente entre la ten-
sión (ecuación 3.17), ahora se calcula obteniendo previamente los valores de impedancia
según las ecuaciones (3.7), (3.8), y (3.9). Es decir:
∗
|Vin |2 |Vin |2
Zi = = θi = |Zi | θi (3.18)
Si |Si |
Este cálculo se ve representado en esta sección del subsistema:

Sabela González Castro 41

 Estudio de modelos de carga en sistemas eléctricos de distribución

Figura 3.28: Ampliación 2 del Subsystem Z.

En este caso se utiliza el módulo de la tensión no medido de la fuente, sino el nominal y

asignado en la máscara externa.

Con el valor de la tensión se hace el cuadrado y el conjugado, tal y como aparece en

la ecuación, para poder dividirlo entre el conjugado de la potencia aparente: S=P-jQ.

Ası́ se obtiene la impedancia constante. A partir de ella y con la tensión medida (sin
su componente wt) se obtiene la intensidad, siendo el resto análogo al modelo anterior.

3.2.3. Modelo intensidad constante.

Como se explicó en el modelo de intensidad constante en régimen fasorial,en este caso
juegan un papel importante los ángulos. La base teórica es la de las ecuaciones (3.13),
(3.14) y (3.15); es decir:
ILi = |ILi | δi − θi (3.19)

 ∗
Si
ILi = (3.20)
Ui
El modelo viene representado en Simulink como aparece en la Figura 3.29:

Sabela González Castro 42

 Estudio de modelos de carga en sistemas eléctricos de distribución

Figura 3.29: Modelo temporal de intensidad constante.

Sabela González Castro 43

 Estudio de modelos de carga en sistemas eléctricos de distribución

Se realizan las máscaras correspondientes de la fuente y subsistema Subsystem I para

posteriores simulaciones. En el primero de la fuente se escogen como variables de entrada
la tensión nominal eficaz de fase de la fuente y la frecuencia de trabajo.
En el segundo se añaden la tensión de la fuente de fase en valor eficaz, la frecuencia de
trabajo, el desfase inicial de la tensión y el valor inicial de la intensidad en p.u.

Figura 3.30: Máscaras de fuente y Subsystem I.

La entrada al Subsystem I son las tres tensiones de fase de la fuente. Como salida se ob-
tienen el módulo de la intensidad y el ángulo de la potencia aparente, φ, que se mantienen
constantes a lo largo del tiempo.

En este caso se obtendrá una onda de la forma:

i(t) = Iamplitud cos(wt + θ − φ) (3.21)
Donde wt en radianes es función de cada instante (variable) , θ es el desfase de la tensión
(variable) y φ es el ángulo de la potencia aparente (constante).

Sabela González Castro 44

 Estudio de modelos de carga en sistemas eléctricos de distribución

Entrando en el Subsystem I se tiene:

Figura 3.31: Máscaras de fuente y Subsystem I.

Sabela González Castro 45

 Estudio de modelos de carga en sistemas eléctricos de distribución

Como en lo dos casos anteriores, la construcción de las 3 fases es idéntica aunque los
resultados sean distintos. Por ello se estudiará, de nuevo, una sola fase.

Ampliando a la primera de las fases se tiene lo representado en la figura 3.32:

Figura 3.32: Máscaras de fuente y Subsystem I.

La parte de la obtención del valor eficaz de la tensión, el desfase y la componente wt es

análoga a la explicada en la sección 3.2.1.

Centrándo el análisis en el cálculo de la intensidad, se tiene:

Figura 3.33: Máscaras de fuente y Subsystem I.

El módulo de la intensidad se obtiene dividiendo el conjugado de la potencia aparente por

el valor nominal de la tensión eficaz asignado en la máscara exterior. Multiplicando por
la ganancia raı́z de 2 se obtiene la amplitud de la intensidad. De esta se mantiene constante.

El ángulo phi (φ) se toma a partir del conjugado de la potencia aparente, por lo que
es en negativo. Este -phi(−φ) se suma al wt y el theta (θ) para hacer la función seno y
mutliplicarla por el módulo de la intensidad.

Sabela González Castro 46

 Estudio de modelos de carga en sistemas eléctricos de distribución

La Figura 3.34 muestra la salida del subsistema. Para obtener el ángulo phi φ se ha
operado como se muestra.

Figura 3.34: Salida de la primera fase del subsistema.

3.3. Generación y aplicación de curvas de carga

Como se mencionó con anterioridad, las primeras simulaciones de los modelos en régimen
temporal se hacı́an para valores de P y Q constantes como en los modelos de régimen
fasorial (ver Figura 3.20).

Para poder acercar estos modelos a la realidad y poder analizar el comportamiento, se ha

estudiado un programa que genera curvas de carga de consumos residenciales.
Como se mencionó en la sección 1.3, Modelos de curvas de carga, las Curvas de carga son
una representación del consumo de potencia respecto del tiempo.

El programa Load Profile Generator [11] se encarga de proporcionar, entre otros, unas
tablas en las que vienen valores de potencia activa y reactiva para cada instante de tiem-
po (ajustable). Este programa trabaja en Windows y es gratuito.
Lo útil y particular del programa es que las curvas se generan atendiendo a muchı́simos
parámetros y variables a programar. Entre ellos se encuentran el tipo de vivienda, la edad
de los residentes, el paı́s (con sus correspondientes dı́as festivos, clima y ciudad), el periodo
vacacional, la situación de los residentes (parado/estudiante/empleado) y el tipo de gas
que consumen.

Sabela González Castro 47

 Estudio de modelos de carga en sistemas eléctricos de distribución

En este Proyecto se ha optado por emplear Curvas de carga correspondientes a consumos

industriales, obtenidas de la web de Red Eléctrica de España [12] .
Se han estudiado tres situaciones para P y Q de uso industrial en la isla del Hierro, Ca-
narias.
En el apartado de resultados se presentan los datos y conclusiones obtenidas para cada
una de ellas.

A continuación aparecen representadas cada una de las situaciones.

En el eje vertical se representa la potencia (W o Var) y en el eje horizontal los intervalos
de tiempo a lo largo de un dı́a.
Se han seleccionado los dás 12 de marzo, 13 e junio y 13 de septiembre. No son fehcas
muy alejadas pero suficientes para observar diferencias en cuanto a temperatura y horas
de sol.

Figura 3.35: Curva de carga, 13 de junio.

El máximo de la primera situación son 20,6 MW de potencia activa y 93,6 MVar de

reactiva.

Sabela González Castro 48

 Estudio de modelos de carga en sistemas eléctricos de distribución

Figura 3.36: Curva de carga, 13 de septiembre.

El máximo de potencia activa es de 21,6 MW, y el de reactiva es de 93,6 MVar.

Figura 3.37: Curva de carga, 12 de marzo.

El máximo se da para 19 MW de potencia activa y 92,02 MVar de reactiva, hacia el final

del dı́a.

Sabela González Castro 49

Capı́tulo 4

Resultados

En este capı́tulo se muestran los resultados y comprobaciones numéricas de los modelos

expuestos en el capı́tulo anterior. Se dividen por secciones según el análisis, pero al final
del capı́tulo en la sección 4.2 se recoge una tabla resumen con todos los valores obtenidos.

4.1. Valores obtenidos y representaciones.

4.1.1. Resultados del flujo de cargas.
Para los mismos valores de entrada y nominales se procede a simular el flujo de cargas tal
y como se indicó en la sección 3.1.1.

Sabela González Castro 50

 Estudio de modelos de carga en sistemas eléctricos de distribución

[Link]. Potencia constante.

Figura 4.1: Resultados para el flujo de cargas con PQ constante.

Las tensiones de fase medidas en el nudo 4 son:

Fase a: U4a = 0,799 −39,08o pu
Fase b: U4b = 0,855 −158,32o pu
Fase c: U4c = 0,8267 80,84o pu
41600
Como la tensión base en el lado de baja tensión del transformador es √ , se tie-
3
nen las siguientes tensiones de fase del nudo 4 medidas en voltios:
Fase a: U4a = 1919 −39,08o V
Fase b: U4b = 2054 −158,32o V
Fase c: U4c = 1986 80,84o V

Sabela González Castro 51

 Estudio de modelos de carga en sistemas eléctricos de distribución

[Link]. Impedancia constante.

Figura 4.2: Resultados para el flujo de cargas con Z constante.

Las tensiones de fase medidas en el nudo 4 son:

Fase a: U4a = 0, 8656 −35, 99o pu
Fase b: U4b = 0, 8859 −156, 29o pu
Fase c: U4c = 0, 8783 83, 13o pu
41600
Como la tensión base en el lado de baja tensión del transformador es √ , se tienen las
3
siguientes tensiones de fase del nudo 4 medidas en voltios:
Fase a: U4a = 2079 −35, 99o V
Fase b: U4b = 2128 −156, 29o V
Fase c: U4c = 2110 83, 13o V

Sabela González Castro 52

 Estudio de modelos de carga en sistemas eléctricos de distribución

[Link]. Intensidad constante.

Figura 4.3: Resultados para el flujo de cargas con I constante.

Se obtienen unas tensiones de fase para el nudo 4:

Fase a: U4a = 0, 8432 −37, 00o pu

Fase b: U4b = 0, 8723 −157, 10o pu
Fase c: U4c = 0, 8601 82, 26o pu
41600
Como la tensión base en el lado de baja tensión del transformador es √ , se tie-
3
nen las siguientes tensiones de fase del nudo 4 medidas en voltios:
Fase a: U4a = 2025 −37, 00o V
Fase b: U4b = 2095 −157, 10o V
Fase c: U4c = 2066 82, 26o V

Sabela González Castro 53

 Estudio de modelos de carga en sistemas eléctricos de distribución

4.1.2. Resultados del modelo fasorial

[Link]. Potencia constante.

Figura 4.4: Resultados para el modelo fasorial con PQ constante.

Fase a: U4a = 1919 −39, 07o V

Fase b: U4b = 2054 −158, 3o V
Fase c: U4c = 1986 80, 85o V

Para ver que funciona consumiendo una potencia impuesta en la máscara, se varı́a el
valor nominal de potencia (que eran 1800 kW) a, por ejemplo, 2000kW y se observa cómo
es ese valor el que se mide en el nudo 4:

Sabela González Castro 54

 Estudio de modelos de carga en sistemas eléctricos de distribución

Figura 4.5: Potencia consumida en la simulación.

Si se modifica el valor de tensión de la fuente, se observa que la potencia consumida se

mantiene constante y tal cual se ha impuesto, aunque varı́e la tensión en el nudo 4.
En la siguiente figura se muestran los resultados de potencia consumidos y tensión medida
en el nudo 4 (esta última rodeada en rojo).

Figura 4.6: Resultados para el modelo fasorial en condiciones nominales para el modelo
de potencia constante.

A continuación se modifica el valor de tensión de la fuente de alimentación, observando

cambios en el valor de tensión medido en el nudo 4, y comprobando que los valores de
potencia consumidos son constantes.

Sabela González Castro 55

 Estudio de modelos de carga en sistemas eléctricos de distribución

Figura 4.7: Tensión en condiciones no nominales para el modelo de potencia constante.

[Link]. Impedancia constante.

Al simular en condiciones nominales el modelo fasorial de impedancia constante, se ob-
tienen los siguientes resultados:

Figura 4.8: Resultados para el modelo fasorial con Z constante.

Sabela González Castro 56

 Estudio de modelos de carga en sistemas eléctricos de distribución

Fase a: U4a = 2079 −35, 98o V

Fase b: U4b = 2128 −156, 3o V
Fase c: U4c = 2110 83, 14o V

El tag Vabc mide la tensión en el nudo 4 para cada fase. Los medidores M odulo y AnguloV
obtienen lo valores de módulo y ángulo por cada fase.

Los valores obtenidos de la impedancia al operar con condiciones nominales son los que
aparecen en la siguiente figura:

Figura 4.9: Resultados para el modelo fasorial en condiciones nominales para el modelo
de impedancia constante.

Z=2.882 25, 84o Ω para cada una de las tres fases.

Si se simula para condiciones no nominales, por ejemplo variando el valor de tensión de la
fuente de tensión, se observa que la tensión medida en el nudo 4 varı́a pero la impedancia
se mantiene tal y como la muestra el medidor Z de la Figura 4.9: constante. Esto no se
añade pues es muy sencillo comprobarlo y la imagen serı́a exactamente la misma.

[Link]. Intensidad constante.

Al simular en condiciones nominales el modelo fasorial de intensidad constante se obtienen
los siguientes valores de tensión medidos en el nudo 4.

Sabela González Castro 57

 Estudio de modelos de carga en sistemas eléctricos de distribución

Figura 4.10: Resultados para el modelo fasorial para el modelo de itensidad constante.

Fase a: U4a = 2025 −37, 0o V

Fase b: U4b = 2095 −157, 1o V
Fase c: U4c = 2066 82, 26o V

La intensidad de trabajo en condiciones nominales es la mostrada en al Figura 4.11.

Figura 4.11: Resultados para el modelo fasorial en condiciones nominales para el modelo
de intensidad constante.

Sabela González Castro 58

 Estudio de modelos de carga en sistemas eléctricos de distribución

Fase a: I4a = 833 −62, 84o A

Fase b: I4b = 833 177, 1o A
Fase c: I4c = 833 56, 42o A

Si se modifican los valores de la fuente de tensión para operar en condiciones no no-

minales, se comprueba que el módulo de la corriente se mantiene constante y el ángulo
varı́a. Esto se muestra en la Figura 4.12.

Figura 4.12: Intensidad en condiciones no nominales para el modelo de intensidad cons-

tante.

Como se puede observar, el ángulo de la fase A se modifica de -62,84o a -59,71o , el de la

fase B varı́a de 177,1o a -179,8o y el de la fase C de 56,42o a 59,86o . En ambos casos el
módulo de la intensidad para las tres fases es de 833 A. El ángulo varı́a ya que solo la
componente debida al factor de potencia se mantiene constante, pero la otra de la tensión
varı́a.

4.1.3. Resultados del modelo temporal.

A continuación se muestran los valores obtenidos de las simulaciones de los modelos. Se
simularán para potencia constante y para curvas de carga. Estas últimas se harán en el
mismo intervalo de tiempo para su posterior análisis.

[Link]. Potencia constante.

El valor de la potencia consumida para valores de P y Q constantes (50 MW y 25 Mvar)
es, efectivamente 50 MW y 25 MVar.

Sabela González Castro 59

 Estudio de modelos de carga en sistemas eléctricos de distribución

Se puede observar que se mantienen constante en el tiempo en la siguiente representación,

donde también se ve el pequeño transitorio del inicio.

Figura 4.13: P(arriba) y Q (abajo) respecto del tiempo.

Resultados para valores de P y Q de los 3 perfiles de carga:

Los resultados obtenidos para la primera situación (dı́a 13 de junio) de las curvas
de carga P11 y Q11 son los que se muestran a continuación.
Las curvas de carga aparecen en intervalos de tiempo modificados para la simplifi-
cación del análisis. Como Simulink estudia cada instante, podrı́a interpolarse para
obtener una mayor exactitud. Sin embargo, como interesa una visión más global y
la comprobación del modelo, no se profundizará y se trabajará con intervalos.

Simulando para 5s, se tienen los siguientes resultados para una de las fases (se
obtiene el mismo en las otras dos restantes ya que la curva de carga es la misma
para las tres fases).

Sabela González Castro 60

 Estudio de modelos de carga en sistemas eléctricos de distribución

Figura 4.14: Simulación de 5 segundos para el perfil 1.

P= 1,518 MW ; Q=0,696 MVar.

Figura 4.15: Valores de P (izquierda) y Q (derecha) de la curva de carga a los 5s.

Como se puede observar los valores medidos en el modelo (Figura 4.13), no son
exactamente iguales que los indicados por el perfil de carga para ese intervalo (Figura
4.14). Se puede explicar si se observan las gráficas donde vienen representadas las
potencias respecto del tiempo. Se tiene:

Sabela González Castro 61

 Estudio de modelos de carga en sistemas eléctricos de distribución

Figura 4.16: P y Q en función del tiempo.

Si se amplia en torno a a los 5 segundos, se puede comprobar que al alcanzar ese

valor, al que le corresponde un nuevo valor de potencia, todavı́a no ha alcanzado el
régimen permanente. Esto explica por qué los valores obtenidos en el medidor no
coinciden exactamente con los marcados en verde en la Figura 4.14, porque el sis-
tema no ha alcanzado el régimen permanente correspondiente al valor de esa Figura.

Sabela González Castro 62

 Estudio de modelos de carga en sistemas eléctricos de distribución

Figura 4.17: P y Q en función del tiempo ampliado.

Es decir, que la constante de tiempo es mayor que el tiempo en el que actúa cada
valor de potencia, por lo que no termina de salir del transitorio.

Los resultados obtenidos para la segunda situación (dı́a 13 de septiembre) son, para
una simulación de 5s:

Figura 4.18: Simulación de 5 segundos para el perfil 2.

P= 1,82 MW ; Q=0,841 MVar.

Sabela González Castro 63

 Estudio de modelos de carga en sistemas eléctricos de distribución

Figura 4.19: Valores de P (izquierda) y Q (derecha) de la curva de carga a los 5s.

Es decir, de nuevo el sistema presenta una constante de tiempo mayor que el tiempo
de trabajo. Lo que se puede observar en las siguientes figuras.

Figura 4.20: P y Q respecto del tiempo.

Sabela González Castro 64

 Estudio de modelos de carga en sistemas eléctricos de distribución

Figura 4.21: P y Q respecto del tiempo ampliado.

Los resultados obtenidos para la tercera situación (dı́a 12 de marzo) son, para una
simulación de 5s:

Figura 4.22: Simulación de 5 segundos para el perfil 3.

P= 1,35 MW ; Q=0,614 MVar.

Sabela González Castro 65

 Estudio de modelos de carga en sistemas eléctricos de distribución

Figura 4.23: Valores de P (izquierda) y Q (derecha) de la curva de carga a los 5s.

Al igual que en los dos casos anteriores, no se consigue alcanzar el régimen perma-
nente en cada intervalo.
Las siguientes figuras lo representan:

Figura 4.24: P y Q respecto del tiempo.

Sabela González Castro 66

 Estudio de modelos de carga en sistemas eléctricos de distribución

Figura 4.25: P y Q respecto del tiempo ampliado.

A lo largo de la simulación se observan valores de potencia constantes y de tensión e

intensidad variables entre cada intervalo. Es decir, entre los segundos 0 y 1 para valores
de (0,1; 0,4; etc) se miden tensiones y corrientes variables pero las potencias activa y
reactiva se mantienen en el valor de la curva de carga correspondiente al intervalo 0.

4.1.4. Impedancia constante.

El valor del módulo de la impedancia obtenido para valores de P y Q constantes (50 MW
y 25 Mvar) es: 4472 Ω.
Se puede observar que se mantiene constante en el tiempo en la siguiente figura:

Sabela González Castro 67

 Estudio de modelos de carga en sistemas eléctricos de distribución

Figura 4.26: Módulo de impedancia respecto del tiempo.

Resultados para valores de P y Q de los 3 perfiles de carga: Simulando para 5s, se tienen los
siguientes resultados para una de las fases (se obtiene el mismo en las otras dos restantes,
ya que la curva de carga es la misma en cada fase).

Los resultados obtenidos para la primera situación (dı́a 13 de junio) son:

Figura 4.27: Valores de impedancia para una simulación de 5s.

La impedancia es de 157 kΩ. Si se observa las gráfica donde viene representadas

la impedancia respecto del tiempo se tiene que no alcanza un valor en régimen
permanente en todos los intervalos (solo en el (1,2) y (3,4)). Para la simulación de
5 segundos se mantiene en régimen transitorio.

Sabela González Castro 68

 Estudio de modelos de carga en sistemas eléctricos de distribución

Figura 4.28: Impedancia respecto del tiempo.

Para la simulación de 5s para el perfil de carga correspondiente a la situación 2 (dı́a

13 de septiembre), se obtienen los siguientes valores:

Figura 4.29: Valores de impedancia para una simulación de 5s.

La impedancia es de 129 kΩ. Si se observa las gráfica donde viene representadas

la impedancia respecto del tiempo se tiene que no alcanza un valor en régimen
permanente, es decir que la constante de tiempo es mayor que el intervalo de trabajo
y no alcanza un valor estable.

Sabela González Castro 69

 Estudio de modelos de carga en sistemas eléctricos de distribución

Figura 4.30: Impedancia respecto del tiempo.

Los resultados obtenidos para la tercera situación (dı́a 12 de marzo) :

Figura 4.31: Valores de impedancia para una simulación de 5s.

La impedancia es de 177 kΩ. Si se observa de nuevo la gráfica donde viene repre-

sentada la impedancia respecto del tiempo se tiene que tampoco alcanza un valor
en régimen permanente en los intervalos (1,2), (2,3) (3,4) y (4,5) es decir que la
constante de tiempo es mayor que el intervalo y no alcanza un valor estable. En los
intervalos (0,1) y (5,6) se alcanza un valor que se mantiene constante, equivalente
al régimen permanente.

Sabela González Castro 70

 Estudio de modelos de carga en sistemas eléctricos de distribución

Figura 4.32: Impedancia respecto del tiempo.

4.1.5. Intensidad constante.

En este caso el módulo de la intensidad y la componente de la fase debida al factor de
potencia se mantienen constantes.

El valor del módulo de la intensidad obtenido para valores de P y Q constantes (50

MW y 25 Mvar) es: 159 A.
Como se observa en la siguiente figura, se mantiene constante.

Sabela González Castro 71

 Estudio de modelos de carga en sistemas eléctricos de distribución

Figura 4.33: Módulo de la intensidad respecto del tiempo.

A continuación se muestran los valores obtenidos con las simulaciones, de nuevo de 5

segundos y para las tres situaciones presentadas de las curvas de carga. El factor de
potencia aparece con el signo negativo ya que aparece restando al resto de componentes
del ángulo de fase de la intensidad, por lo que su valor es al contrario.

Simulación de la situación 1 correspondiente al dı́a 13 de junio:

Figura 4.34: Valores de intensidad y factor de potencia para una simulación de 5s.

Intensidad de 4,95 A y factor de potencia 0,451. Si se observa de nuevo la gráfica

donde viene representada la intensidad respecto del tiempo se tiene que no alcanza
un valor en régimen permanente (salvo los intervalos (1,2) y (3,4)), es decir que la
constante de tiempo es mayor que el intervalo de trabajo.

Sabela González Castro 72

 Estudio de modelos de carga en sistemas eléctricos de distribución

Figura 4.35: Intensidad respecto del tiempo.

Y más concretamente ampliando en torno a los 5 segundos:

Figura 4.36: Intensidad respecto del tiempo.

Simulación 2 del dı́a 13 de septiembre:

Sabela González Castro 73

 Estudio de modelos de carga en sistemas eléctricos de distribución

Figura 4.37: Valores de intensidad y factor de potencia para una simulación de 5s.

Intensidad de 5,45 A y factor de potencia 0,[Link] se observa de nuevo la gráfica

donde viene representada la intensidad respecto del tiempo se tiene que no alcanza
un valor en régimen permanente, es decir que la constante de tiempo es mayor que
el intervalo de trabajo y no alcanza un valor estable.

Figura 4.38: Intensidad respecto del tiempo.

Y más concretamente ampliando en torno a los 5 segundos:

Sabela González Castro 74

 Estudio de modelos de carga en sistemas eléctricos de distribución

Figura 4.39: Intensidad respecto del tiempo.

Simulación 3 del dı́a 12 de marzo:

Figura 4.40: Valores de intensidad y factor de potencia para una simulación de 5s.

Intensidad de 4,19 A y factor de potencia 0,451. Como en los casos anteriores, se

muestra la intensidad respecto del tiempo.

Sabela González Castro 75

 Estudio de modelos de carga en sistemas eléctricos de distribución

Figura 4.41: Intensidad respecto del tiempo.

Y más concretamente ampliando en torno a los 5 segundos:

Figura 4.42: Intensidad respecto del tiempo.

Se observa que se obtiene un valor permanente a partir de los 5 segundos.

Sabela González Castro 76

 Estudio de modelos de carga en sistemas eléctricos de distribución

4.2. Esquema y comparación de resultados.

4.2.1. Flujo de cargas y análisis fasorial.
En esta sección se incluyen unas tablas resumen de los resultados recopilados en las simu-
laciones.

En primer lugar se presentan las tensiones medidas en el nudo 4 de la carga corres-

pondientes a cada modelo ejecutando el flujo de cargas y el análisis fasorial.

Tabla 4.1: Valores obtenidos del flujo de cargas y análisis fasorial.

Como se puede observar, los valores coinciden.

A continuación se incluye una tabla recopilatoria de los valores del análisis fasorial en
distintas condiciones de trabajo.
Cada modelo se ha ejecutado con dos valores de tensión de lı́nea de alimentación: 12470
V (nominal) y 22470 V.
El objetivo de esta tabla es mostrar que, aún variando las condiciones de trabajo, pa-
ra cada modelo se mantiene constante la variable deseada. Es decir, para el modelo de
potencia constante se obtienen distintas intensidades y tensiones medidas si se varı́a la
tensión de alimentación. Sin embargo la potencia (coloreada en naranja) se mantiene igual.

Esto mismo ocurre con los otros dos modelos de impedancia e intensidad constante. Al
trabajar variando la tensión de la fuente de tensión, la impedancia se mantiene constante
de valor 25,8 Ω y la tensión en el nudo 4 varı́a.
En el caso de intensidad constante se mantienen el valor del módulo (833 A, como se
observa) y la componente de la fase del ángulo correspondiente al factor de potencia. En

Sabela González Castro 77

 Estudio de modelos de carga en sistemas eléctricos de distribución

la tabla se ve que el ángulo de fase varı́a, y es porque está incluida la componente debida
a la tensión.

Sabela González Castro 78

Tabla 4.2: Valores obtenidos del análisis fasorial.
 Estudio de modelos de carga en sistemas eléctricos de distribución

4.2.2. Análisis temporal.

La siguiente y última tabla que se muestra resume las simulaciones del análisis temporal.
El funcionamiento es el mismo, para cada modelo se realizan simualaciones donde se
comprueba que cada variable se mantiene constante.
Sin embargo hay una diferencia con la tabla anterior. En los modelos fasoriales se ha
trabajado con valores de potencia constantes. Para el análisis temporal se han empleado
los perfiles de carga, en concreto el correspondiente a la situación 3 para la realización de
la tabla.
Se han realizado por cada modelo cuatro simulaciones. Dos de ellas de 1 segundo (para
tensiones de alimentación de 500 y 600 kV) y el resto de 2 segundos, para las mismas
tensiones.
Ası́, se ha comprobado cómo se mantienen constantes las variables en cada intervalo de
simulación. Es en cada intervalo ya que al tratarse de perfiles de carga, los valores de
potencia van cambiando.

Sabela González Castro 80

Tabla 4.3: Valores obtenidos del análisis temporal.
 Estudio de modelos de carga en sistemas eléctricos de distribución

A continuación se incluye una posible reflexión respecto a los valores obtenidos.

Se observa que el máximo de potencia se da para el dı́a 13 de septiembre (1,82 MW

y 0,69 Mvar), situación 2. Esto puede ser debido a que en la isla del Hierro hay altas
temperaturas y el consumo puede destinarse al aire acondicionado y también a la luz, ya
que en el mes de septiembre las horas de sol son más reducidas que en junio.
El siguiente dı́a en el que se consume mayor potencia es el correspondiente a la situación
1, el dı́a 13 de junio (1,518 MW y 0,69 MVar). En este caso las temperaturas también
son muy elevadas pero las horas de sol son mayores, por lo que prácticamente el consumo
podrı́a ir destinado al aire acondicionado.
El dı́a que se consume menos es el 12 de marzo, situación 3 (1,35 MW y 0,614 MVar). La
isla del Hierro presenta altas temperaturas en verano y moderadas en invierno. Por este
motivo se puede pensar que no hay consumo prácticamente dedicado a la calefacción y
que casi todo se emplea en la luz por el menor número de horas solares.

Este análisis son meras reflexiones de acorde a los resultados obtenidos. Esto es debi-
do a que el clima es muy aleatorio y depende del dı́a en que se tomen las mediciones y
del momento de la simulación. Sin embargo, para el propósito de este Proyecto de Fin de
Grado no interesa profundizar más sobre ello.

Habiendo realizado las tres simulaciones para el mismo intervalo de tiempo, se tiene que
para el dı́a 13 de junio la impedancia es de 157 kΩ , para el 13 de septiembre de 129 kΩ
y para el 12 de marzo de 177 kΩ.
Antes se tenı́a que la mayor potencia se consumı́a para la situación 2, luego la 1 y final-
mente la 3. En este caso, la impedancia mayor se obtiene para la situación 3, frente a la
1 y, finalmente la 2.
Esto es lo coherente, ya que la potencia y la impedancia son inversamente proporcionales,
de acuerdo con:
|Vin |2
Zi = (4.1)
Si ∗
La justificación del orden serı́a la misma que la de la potencia, pero a la inversa.

Se comprueba que modificando la tensión de la alimentación, la impedancia sigue siendo

la misma.

Se obtiene que la intensidad mayor se da para la situación 2 (13 de septiembre) con

un valor de 5,54 A. A continuación se encuenytra la situación 1 (13 de junio) con 4,96 A
y finalmente la situación 3 (12 de marzo) con 4,19 A.

Sabela González Castro 82

 Estudio de modelos de carga en sistemas eléctricos de distribución

Como la intensidad es directamente proporcional a la potencia según la ecuación:

Sc ∗
 
ILc = (4.2)
Vcn
Se aplica la justificación de la sección 5.2.1 en el modelo de potencia constante.

Se comprueba que variando la tensión de alimentación se obtienen valores de intensi-

dad iguales para cada intervalo.

Sabela González Castro 83

Capı́tulo 5

Conclusiones

En este capı́tulo se comentan las conclusiones del Proyecto.

Se ha simulado un sistema de distribución con sus desequilibrios correspondientes en la
red. Para ello se ha empleado la herramienta Simulink de Matlab, la cual se ha compro-
bado que es versátil para estos estudios.

El primer análisis que se ha realizado ha sido el del flujo de cargas para el sistema de
cuatro nudos. En este análisis Simulink aporta una carga que permite seleccionar si se
quiere que funcione consumiendo potencia, impedancia o intensidad constante.
Una vez obtenidos los datos, se han implementado los modelos fasoriales para el mismo
sistema de cuatro nudos. La simulación de estos se corresponde con la ejecución de in-
finitos flujos de carga en el tiempo. En este caso, Simulink no presenta la opción de la
carga anterior. Por eso se han desarrollado modelos de carga que permiten simular ese
comportamiento.
Con ellos se ha podido comprobar que los resultados obtenidos para un punto por ambos
métodos (flujo de cargas y análisis fasorial) coinciden.

Una vez se ha terminado de analizar los modelos anteriores se han desarrollado tam-
bién en Simulink 3 modelos de carga. Estos son similares a los del modelo fasorial pero
adaptándolos de manera que permiten estudiar el régimen transitorio. Se han analizado
los resultados para valores de potencia constantes y para perfiles de carga en cada uno de
los tres modelos.

Sabela González Castro 84

Capı́tulo 6

Impacto Social y Ambiental.

El estudio del diseño y trabajo de las redes desequilibradas y Smart Grids es algo continuo.
Está claro que es debido a que su mejora supone un acercamiento y mayores beneficios.

Esto presenta consecuencias a nivel social, económico y ambiental.

El hecho de dominar la red correctamente no tiene un impacto ambiental ”directo”. Se

podrı́a decir que un buen funcionamiento podrı́a aumentar el uso de energı́as renovables
frente a comubstibles fósiles, reduciendo ası́ la contaminación y acelerando la transición
energética hacia un futuro de mayor limpieza y sostenibilidad.

El mayor uso de energı́a renovable supone un impacto económico. Ahora mismo, los com-
bustibles fósiles son lı́deres en el consumo de su energı́a. Además, las energı́as renovables
no se ven sometidas a las fluctuaciones de precio como los combustibles fósiles, lo que
supone mayor estabilidad.
Una vez construidas las infraestructuras, el uso de la energı́a eléctrica no implica la emi-
sión de contaminantes o gases que contribuyen al efecto invernadero. Si es cierto que tiene
sus efectos adversos, sobre todo ante la fauna. Además la aportación es más complicada
de predecir debido a que dependen del clima y la demanda. Sin embargo, la inteligencia
artificial está poniendo fin a estos inconvenientes. [13]

Por otro lado, es mucho más viable proporcionar energı́a renovable frente a comprar com-
bustibles fósiles en los paı́ses en vı́as de desarrollo.

Si se tuviese la opción de proporcionar a los usuarios la energı́a eléctrica adecuada en

el tiempo necesario, se podrı́an minimizar las pérdidas y, por tanto, reducir las tarifas

Sabela González Castro 85

 Estudio de modelos de carga en sistemas eléctricos de distribución

de la energı́a consumida. Además se presentarı́a una mayor libertad de elección para los
usuarios. El hecho de no emitir gases tóxicos mejorarı́a la salud de los habitantes expues-
tos a ellos.

Sabela González Castro 86

Capı́tulo 7

Planificación Temporal y
Presupuesto

7.1. Planificación temporal

En esta sección se muestra la estructuración del proyecto y la distribución temporal del
mismo.

En primer lugar se expone la Estructura de Descomposición del Proyecto (EDP). Es-

ta divide en los principales bloques fundamentales del Proyecto, que aparecen en verde.
Debajo de los bloques verdes estan los sub-bloques azules que, dentro de cada bloque fun-
damental representan el primer paso. Los sub-bloques naranjas son las tareas realizadas.

Sabela González Castro 87

 Estudio de modelos de carga en sistemas eléctricos de distribución

Figura 7.1: Estructura de Descomposición del Proyecto.

A continuación se procede a presentar el diagrama de Gantt.

Gracias a el se puede entender de manera visual la duración del Proyecto. Esta comienza
el dı́a de la asignación del Proyecto, el 16 de octubre de 2017, y finaliza el dı́a 18 de julio
de 2018.
Entre estas dos fechas se ha dedicado al Proyecto un total de 320 horas, aproximadamente.

La siguiente figura muestra el diagrama de Gantt, que se ha podido realizar considerando

lo anterior.

Sabela González Castro 88

 Estudio de modelos de carga en sistemas eléctricos de distribución

Figura 7.2: Diagram de Gantt.

7.2. Presupuesto
En esta sección se indicarán brevemente los costes del proyecto y el cálculo del presupues-
to. Al final consta una tabla donde vienen representados los costes con IVA (21 %) sin
IVA .
Para calcular el presupuesto se ha tenido en cuenta lo siguiente. En primer lugar, el uso del
ordenador: se ha estimado que, de 320 horas de proyecto, solo 20 han sido sin ordenador
correspondientes a una fase de documentación en papel. La vida útil del ordenador son
10 años y el coste del ordenador 1000e.
El tiempo que se ha estado usando han sido 8 meses en los que se ha desarrollado el
trabajo, con un factor de uso aproximadamente del 45 %. Teniendo en cuenta que el IVA
es del 21 % se tendrı́a un coste con IVA de 30e.

En segundo lugar, los sueldos: se ha considerado que el tutor tiene un contrato de 3500e/mes
y el alumno 900e/mes. Estimando un trabajo diario de 8 h que implica aproximadamente
22e/h y 5.6e/h para tutor y alumno respectivamente. El alumno ha trabajado un total
de 320 horas, por lo que le corresponden 1792e. El tutor ha trabajado 45h, por lo que le
corresponden 990e.

En tercer lugar, la energı́a eléctrica consumida: solo se ha considerado lo que consume

el ordenador en las 300 horas de trabajo para un consumo de 0.125kW/h y un coste de

Sabela González Castro 89

 Estudio de modelos de carga en sistemas eléctricos de distribución

0.1e/kW, lo que viene a resultar un coste de 3.75e.

En cuarto y último lugar, los costes de desplazamiento: siendo de 20e/mes para 8 meses
de trabajo, resultando 160e. No se han considerado los costes de programas por tener la
licencia MATLAB de estudiante y ser el TeXstudio y el Load Profile Generator gratuitos.

Tabla 7.1: Tabla de presupuesto.

Sabela González Castro 90

Capı́tulo 8

Trabajos futuros.

En primer lugar podrı́a realizarse un estudio más exhaustivo del modelo temporal, ya
que en este Proyecto ha estado limitado por la constante de tiempo. Como se veı́a en el
Capı́tulo 4 para las simulaciones con perfiles de carga, la constante de tiempo era mayor
que el intervalo de trabajo y, por ello, no se llegaba a alcanzar el régimen permanente
para cada valor.

Como se comentó en la introducción, además de las cargas estáticas se encuentran las

cargas dinámicas. En este proyecto se ha trabajado solo con modelos de las primeras.
Por ello una posible continuación serı́a desarrollar unos modelos que permitiesen estudiar
el comportamiento ante este tipo de cargas más complejas. Además ası́ podrı́a no solo
compararse modelos dentro de las cargas dinámicas, como se hace en este Proyecto para
los estáticos: potencia, impedancia e intensidad; sino comparar el comportamiento de las
cargas dinámicas frente a las estáticas.

Sabela González Castro 91

Bibliografı́a

[1] Red Eléctrica de España. ”Smart Grids”. Internet: [Link]

inteligentes/que-son-las-smartgrid
[2] William H. Kersting. Ïntroduction to Distribution System.”Distribution System Mo-
deling and Analysis. Las Cruces, New Mexico: CRC Press.
[3] Prabha Kundur. ”Power System Loads”. Power System Stability And Control. CA:
McGraw-Hill, pp 271-274.
[4] Prabha Kundur. ”Power System Loads”. Power System Stability And Control. CA:
McGraw-Hill, pp 274-278.
[5] Aigner DJ, Sorooshian C, Kerwin P. Conditional demand analysis for estimating re-
sidential end-use load profiles. The Energy Journal 1984, pp 81-97.
[6] Bartels R, Fiebig DG, Garben M, Lumsdaine R. An end-use electicity load simulation
model. Delomd. Utilities Policy 1992, pp 71-82.
[7] Yao R, Steemers K. A method of formulating energy load profile for domestic buildings
in the UK. Energy and Buildings 2005, pp 663-71.
[8] Stokes M. Removing barriers to embedded generation: a fine-grained load model to
support low voltage networl performance analyisis. PhD [Link] Montfort University;
2005.
[9] Roger Dugan. ÏEE 4 Node Test Feeder”. Internet: [Link], Sept. 2006.
[10] Daniel Rabinovich y Oscar Santa Cruz. ”PLL: Phase Loop Locked Ïnternet:
[Link]
2010.
[11] ”Load Profile Generator”. Internet: [Link]
[12] Red Eléctrica de España. Internet: [Link]

Sabela González Castro 92

 Estudio de modelos de carga en sistemas eléctricos de distribución

[13] Javier Martı́n Arroyo, ”U na solución para dominar los picos de la energı́a verde”,
[Link] Jun.
22, 2018.

Sabela González Castro 93

Índice de figuras

1.1. Componentes de sistemas de energı́a eléctrica. . . . . . . . . . . . . . . . . 1

1.2. Componentes de una red de distribución . . . . . . . . . . . . . . . . . . . 3
1.3. Descomposición en componentes simétricas. . . . . . . . . . . . . . . . . . 4
1.4. Esquema de representación de los modelos de curvas de carga. . . . . . . . 7

2.1. Relación de Smart Grid. . . . . . . . . . . . . . . . . . . . . . . . . . . . . 11

3.1. Modelo de 4 nudos para el flujo de cargas. . . . . . . . . . . . . . . . . . . 14

3.2. Desplegable de la carga. . . . . . . . . . . . . . . . . . . . . . . . . . . . . 14
3.3. Obtención de valores del flujo de cargas. . . . . . . . . . . . . . . . . . . . 15
3.4. Carga conectada en estrella. . . . . . . . . . . . . . . . . . . . . . . . . . . 15
3.5. Modelo fasorial de potencia constante. . . . . . . . . . . . . . . . . . . . . 18
3.6. Comparación de valores obtenidos en dos puntos. . . . . . . . . . . . . . . 20
3.7. Ampliación de la parte modificada. . . . . . . . . . . . . . . . . . . . . . . 21
3.8. Máscara del Subsystem P Q cte. . . . . . . . . . . . . . . . . . . . . . . . . 22
3.9. Interior del Subsystem P Q cte. . . . . . . . . . . . . . . . . . . . . . . . . 22
3.10. Modelo fasorial de impedancia constante. . . . . . . . . . . . . . . . . . . . 24
3.11. Máscara del Subsystem Z cte. . . . . . . . . . . . . . . . . . . . . . . . . . 25
3.12. Subsystem Z cte. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 26
3.13. Modelo fasorial de intensidad constante. . . . . . . . . . . . . . . . . . . . 27
3.14. Máscara del Subsystem I cte. . . . . . . . . . . . . . . . . . . . . . . . . . 28
3.15. Subsystem I cte. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 28
3.16. Modelo temporal de potencia constante. . . . . . . . . . . . . . . . . . . . 30
3.17. Máscaras de fuente y Subsystem P Q. . . . . . . . . . . . . . . . . . . . . . 31
3.18. Máscaras de fuente y Subsystem P Q. . . . . . . . . . . . . . . . . . . . . . 32
3.19. Fase A del Subsystem P Q. . . . . . . . . . . . . . . . . . . . . . . . . . . . 33
3.20. Bloques P,Q constantes (rosas) y From Workspace. . . . . . . . . . . . . . 33
3.21. Ampliación de la Figura 3.19. . . . . . . . . . . . . . . . . . . . . . . . . . 34
3.22. Circuito PLL. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 35
3.23. Ampliación de la Figura 3.19. . . . . . . . . . . . . . . . . . . . . . . . . . 36

Sabela González Castro 94

 Estudio de modelos de carga en sistemas eléctricos de distribución

3.24. Modelo temporal de impedancia constante. . . . . . . . . . . . . . . . . . . 38

3.25. Máscaras de fuente y Subsystem Z. . . . . . . . . . . . . . . . . . . . . . . 39
3.26. Subsystem Z. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 40
3.27. Ampliación del Subsystem Z. . . . . . . . . . . . . . . . . . . . . . . . . . 41
3.28. Ampliación 2 del Subsystem Z. . . . . . . . . . . . . . . . . . . . . . . . . 42
3.29. Modelo temporal de intensidad constante. . . . . . . . . . . . . . . . . . . 43
3.30. Máscaras de fuente y Subsystem I. . . . . . . . . . . . . . . . . . . . . . . 44
3.31. Máscaras de fuente y Subsystem I. . . . . . . . . . . . . . . . . . . . . . . 45
3.32. Máscaras de fuente y Subsystem I. . . . . . . . . . . . . . . . . . . . . . . 46
3.33. Máscaras de fuente y Subsystem I. . . . . . . . . . . . . . . . . . . . . . . 46
3.34. Salida de la primera fase del subsistema. . . . . . . . . . . . . . . . . . . . 47
3.35. Curva de carga, 13 de junio. . . . . . . . . . . . . . . . . . . . . . . . . . . 48
3.36. Curva de carga, 13 de septiembre. . . . . . . . . . . . . . . . . . . . . . . . 49
3.37. Curva de carga, 12 de marzo. . . . . . . . . . . . . . . . . . . . . . . . . . 49

4.1. Resultados para el flujo de cargas con PQ constante. . . . . . . . . . . . . 51

4.2. Resultados para el flujo de cargas con Z constante. . . . . . . . . . . . . . 52
4.3. Resultados para el flujo de cargas con I constante. . . . . . . . . . . . . . . 53
4.4. Resultados para el modelo fasorial con PQ constante. . . . . . . . . . . . . 54
4.5. Potencia consumida en la simulación. . . . . . . . . . . . . . . . . . . . . . 55
4.6. Resultados para el modelo fasorial en condiciones nominales para el modelo
de potencia constante. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 55
4.7. Tensión en condiciones no nominales para el modelo de potencia constante. 56
4.8. Resultados para el modelo fasorial con Z constante. . . . . . . . . . . . . . 56
4.9. Resultados para el modelo fasorial en condiciones nominales para el modelo
de impedancia constante. . . . . . . . . . . . . . . . . . . . . . . . . . . . . 57
4.10. Resultados para el modelo fasorial para el modelo de itensidad constante. . 58
4.11. Resultados para el modelo fasorial en condiciones nominales para el modelo
de intensidad constante. . . . . . . . . . . . . . . . . . . . . . . . . . . . . 58
4.12. Intensidad en condiciones no nominales para el modelo de intensidad cons-
tante. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 59
4.13. P(arriba) y Q (abajo) respecto del tiempo. . . . . . . . . . . . . . . . . . . 60
4.14. Simulación de 5 segundos para el perfil 1. . . . . . . . . . . . . . . . . . . . 61
4.15. Valores de P (izquierda) y Q (derecha) de la curva de carga a los 5s. . . . . 61
4.16. P y Q en función del tiempo. . . . . . . . . . . . . . . . . . . . . . . . . . 62
4.17. P y Q en función del tiempo ampliado. . . . . . . . . . . . . . . . . . . . . 63
4.18. Simulación de 5 segundos para el perfil 2. . . . . . . . . . . . . . . . . . . . 63
4.19. Valores de P (izquierda) y Q (derecha) de la curva de carga a los 5s. . . . . 64
4.20. P y Q respecto del tiempo. . . . . . . . . . . . . . . . . . . . . . . . . . . . 64

Sabela González Castro 95

 Estudio de modelos de carga en sistemas eléctricos de distribución

4.21. P y Q respecto del tiempo ampliado. . . . . . . . . . . . . . . . . . . . . . 65

4.22. Simulación de 5 segundos para el perfil 3. . . . . . . . . . . . . . . . . . . . 65
4.23. Valores de P (izquierda) y Q (derecha) de la curva de carga a los 5s. . . . . 66
4.24. P y Q respecto del tiempo. . . . . . . . . . . . . . . . . . . . . . . . . . . . 66
4.25. P y Q respecto del tiempo ampliado. . . . . . . . . . . . . . . . . . . . . . 67
4.26. Módulo de impedancia respecto del tiempo. . . . . . . . . . . . . . . . . . 68
4.27. Valores de impedancia para una simulación de 5s. . . . . . . . . . . . . . . 68
4.28. Impedancia respecto del tiempo. . . . . . . . . . . . . . . . . . . . . . . . . 69
4.29. Valores de impedancia para una simulación de 5s. . . . . . . . . . . . . . . 69
4.30. Impedancia respecto del tiempo. . . . . . . . . . . . . . . . . . . . . . . . . 70
4.31. Valores de impedancia para una simulación de 5s. . . . . . . . . . . . . . . 70
4.32. Impedancia respecto del tiempo. . . . . . . . . . . . . . . . . . . . . . . . . 71
4.33. Módulo de la intensidad respecto del tiempo. . . . . . . . . . . . . . . . . . 72
4.34. Valores de intensidad y factor de potencia para una simulación de 5s. . . . 72
4.35. Intensidad respecto del tiempo. . . . . . . . . . . . . . . . . . . . . . . . . 73
4.36. Intensidad respecto del tiempo. . . . . . . . . . . . . . . . . . . . . . . . . 73
4.37. Valores de intensidad y factor de potencia para una simulación de 5s. . . . 74
4.38. Intensidad respecto del tiempo. . . . . . . . . . . . . . . . . . . . . . . . . 74
4.39. Intensidad respecto del tiempo. . . . . . . . . . . . . . . . . . . . . . . . . 75
4.40. Valores de intensidad y factor de potencia para una simulación de 5s. . . . 75
4.41. Intensidad respecto del tiempo. . . . . . . . . . . . . . . . . . . . . . . . . 76
4.42. Intensidad respecto del tiempo. . . . . . . . . . . . . . . . . . . . . . . . . 76

7.1. Estructura de Descomposición del Proyecto. . . . . . . . . . . . . . . . . . 88

7.2. Diagram de Gantt. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 89

Sabela González Castro 96

Índice de tablas

4.1. Valores obtenidos del flujo de cargas y análisis fasorial. . . . . . . . . . . . 77

4.2. Valores obtenidos del análisis fasorial. . . . . . . . . . . . . . . . . . . . . . 79
4.3. Valores obtenidos del análisis temporal. . . . . . . . . . . . . . . . . . . . . 81

7.1. Tabla de presupuesto. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 90

Sabela González Castro 97