Definición 3 (PPISEP) Un g-ésimo Portafolio de Proyectos de Inversión en un Sistema Eléctrico de Potencia (PPISEP) es un conjunto finito P_g de q-ésimos proyectos de inversión, de acuerdo a la Definición 2. Según las definiciones y notaciones empleadas previamente, un P_g puede ser expresado en términos matemáticos como:

donde:

• n: cantidad total de proyectos de inversión,

• t_p = (t_p1,t_p2,...t_pn)^t: tipos o clases de cada uno de los proyectos de inversión (proyectos de generación o de transmisión),

• N_i = (N_i1,N_i2,...N_in)^t Î {1,2,...,N_g+N_t+1}ⁿ: Nodos de Cola de los arcos nominales de la Red R que son alterados por cada uno de los proyectos de inversión,

• N_f = (N_f1,N_f2,...N_fn)^t Î {2,3,...,N_g+N_t+N_c+1}ⁿ: Nodos de Cabeza de los arcos nominales de la Red R que son alterados por cada uno de los proyectos de inversión,

• e_p = (e_p1,e_p2,...e_pn)^t Î Rⁿ: energías activas promedio que introducen cada uno de los proyectos de inversión a la Red R,

• t_i = (t_i1,t_i2,...t_in)^t Î Zⁿ: intervalos de tiempo en los cuales se inicia la construcción de cada uno de los proyectos de inversión.

Definición 4 (Demanda insatisfecha del SEP) Para el presente artículo la demanda insatisfecha viene representada por la diferencia entre la energía activa promedio efectivamente entregada a los consumidores menos el consumo de energía activa promedio proyectado o pronosticado de éstos. Considerando la representación del SEP como una Red R, la diferencia de los valores de los arcos efectivos de consumo respecto a los valores de los arcos nominales de consumo representa la demanda insatisfecha. Para poder representar el efecto que la demanda insatisfecha tiene sobre la empresa eléctrica, su valor debe ser convertido en dinero mediante el uso de factores de penalización de demanda insatisfecha. Cada tipo de consumidor tiene asociado su respectivo factor de penalización de demanda insatisfecha. Para cada intervalo de tiempo se debe calcular el valor de demanda insatisfecha de cada consumidor (en términos de dinero) para luego sumarlos y, de esta manera, obtener el valor de demanda insatisfecha total (en términos de dinero) de cada intervalo de tiempo. Los valores de demanda insatisfecha total (en términos de dinero) de cada intervalo de tiempo deben ser llevados a valor presente (usando tasas fijas de interés y de inflación promedio por intervalo de tiempo) donde pueden ser sumados y, de esta forma obtener el valor de demanda insatisfecha total (en términos de dinero) para el Horizonte de Planificación.

Definición 5 (Balance económico de la empresa eléctrica) La evaluación del balance económico de la empresa eléctrica consiste en determinar los ingresos menos los egresos de la misma durante el Horizonte de Planificación establecido. Los ingresos para cada intervalo de tiempo se obtienen mediante tarifas aplicadas a la energía activa promedio efectivamente entregada en cada intervalo de tiempo a los consumidores (valores de los arcos efectivos de consumo), donde cada tipo de consumidor tiene su tarifa respectiva. Por su parte, los egresos para cada intervalo de tiempo se obtienen mediante la adición de los costos de los proyectos de inversión, los costos variables y los costos fijos. Los costos de los proyectos de inversión se consideran dependientes de manera lineal de la cantidad de energía activa promedio que éstos introducen al SEP y son considerados únicamente en el intervalo de tiempo de inicio de su construcción. Por otra parte, los costos variables dependen de la energía activa efectivamente generada por el SEP y están presentes en cada intervalo de tiempo. Considerando la representación del SEP como una Red R, el cálculo de los costos variables para cada intervalo de tiempo se obtiene mediante la suma de la aplicación de un costo por unidad de energía activa generada a cada uno de los valores de los arcos efectivos de generación. Por último, los costos fijos, los cuales representan los costos de operación y mantenimiento, los compromisos adquiridos y los cargos por depreciación son constantes y están presentes en cada intervalo de tiempo del Horizonte de Planificación. Luego de determinar los ingresos y los egresos en cada uno de los intervalos de tiempo que componen el Horizonte de Planificación, se debe realizar una suma aritmética de los mismos en valor presente, utilizando tasas fijas de interés mínimas atractivas e inflación promedio por intervalo de tiempo, obteniéndose de esta manera el valor del balance económico de la empresa eléctrica para el Horizonte de Planificación.

Sean e_p Î Rⁿ el conjunto de n variables reales independientes que representan las energías activas promedio que introducen a la Red R cada uno de los proyectos de inversión de un PPISEP, y t_i Î Zⁿ el conjunto de n variables enteras independientes que representan los intervalos de tiempo en los cuales se inicia la construcción de cada uno de los proyectos de inversión del PPISEP. Sean además B_p el valor promedio de las muestras de balance económico de la empresa eléctrica para el Horizonte de Planificación H_p, y D_p el valor promedio de las muestras de demanda insatisfecha total (en términos de dinero) para el Horizonte de Planificación H_p, siendo este último menor o igual a cero. Entonces, el presente problema de optimización mediante la simulación puede ser expresado de la siguiente manera:

sujeto a:

donde x₁ Î (0,1), x₂ Î (0,1) son los factores de ponderación, tales que x₁ + x₂ = 1, R_d es el presupuesto que se dispone en valor presente para los proyectos de inversión,(e_p[q])|_p es el costo de inversión del q-ésimo proyecto de inversión referido a su valor presente, siendo el costo por unidad de energía activa promedio que introduce el proyecto de inversión y, = si el q-ésimo proyecto de inversión es de generación o = si el q-ésimo proyecto de inversión es de transmisión;

donde t_i[q] + [e_p[q]] es el intervalo de tiempo de puesta en marcha o ejecución del q-ésimo proyecto de inversión, siendo un factor que multiplica la energía activa promedio que introduce el proyecto de inversión y, = si el q-ésimo proyecto de inversión es de generación o = si el q-ésimo proyecto de inversión es de transmisión. El operador [·] es utilizado para aproximar el resultado a su próximo superior entero;

e_p[q] ≥ 0,     " q Î {1,2,...,n},                (7)

Como se puede observar, las variables de decisión del presente problema de optimización son las energías activas promedio e_p que introducen cada uno de los proyectos de inversión a la Red R (conjunto de variables reales independientes) y los intervalos de tiempos ti en los cuales dichos proyectos se empiezan a construir (conjunto de variables enteras independientes).

Algoritmo Simplex Entero Mixto - Una de las técnicas de optimización utilizadas para la resolución del problema antes planteado es el Algoritmo Simplex Entero Mixto (ASEM) puesto que éste permite considerar el campo numérico entero del intervalo de tiempo de inicio de construcción y el campo numérico real de la cantidad de energía activa promedio de cada uno de los proyectos de inversión que constituyen el PPISEP. El ASEM, denominado de esta manera por el autor Brea (2013), es una extensión del método algorítmico simplex de Nelder Mead (Nelder et al.,1965) para la identificación de al menos un óptimo local, cuando es usado en problemas no lineales enteros mixtos irrestrictos. El ASEM se basa en una doble estructura de símpleces, la cual está compuesta de una estructura simplex real de dimensión n (simplex real), y otra estructura simplex entera de igual dimensión n (simplex entero). Las operaciones propuestas originalmente por Nelder et al., (1965) son ejecutadas sobre el simplex real mientras que otro grupo de nuevas operaciones desarrolladas por Brea (2013) son empleadas en el simplex entero. Este conjunto de nuevas operaciones junto con las operaciones originales del método de Nelder Mead garantizan que en cada iteración del ASEM se arroje un nuevo punto de prueba definido en el campo numérico de dimensión 2n enteros mixtos Rⁿ x Zⁿ, con el fin de garantizar la identificación del óptimo local entero mixto, sin convertir las variables enteras en variables reales. Es muy importante destacar que debido a que el ASEM es un algoritmo de minimización, el problema de la mejora del PPISEP debe ser replanteado como un problema de minimización.

Optimización por camino de muestras - La otra técnica o estrategia de optimización empleada para la resolución del problema antes planteado es la optimización por camino de muestras. En la optimización por camino de muestras el valor esperado de la función objetivo se calcula mediante el valor promedio de una muestra relativamente grande de observaciones en cada punto (camino de muestras) (Carson et al., 1997), (Azadivar, 1999), para luego hacer uso de métodos deterministas de programación no lineal (Carson et al., 1997). De esta manera, esta técnica de optimización permite considerar las incertidumbres aleatorias asociadas a la función objetivo del Problema.

El presente artículo utiliza el ASEM como método determinista de programación no lineal en la optimización por camino de muestras, de esta forma se considera la naturaleza entera-mixta del problema y las incertidumbres asociadas al mismo.

REPLANTEAMIENTO DEL PROBLEMA

Con el propósito de redefinir el problema como un problema de minimización irrestricto y, de esta manera poder utilizar el algoritmo de optimización ASEM, la ecuación (3) es replanteada de la siguiente manera:

y cada una de las restricciones expresadas en las ecuaciones (4), (5), (6) y (7) es replanteada como una función de penalización:

donde j_p es el coeficiente de las funciones de penalización.

MODELO COMPUTACIONAL PARA RESOLVER EL PROBLEMA PLANTEADO

A continuación se presentan los algoritmos, procedimientos y funciones que conforman el modelo computacional diseñado para resolver el Problema anteriormente planteado.

Algoritmo 1. Metodología que permite mejorar los Portafolios de Proyectos de Inversión en un Sistema Eléctrico de Potencia (PPISEP). Este es el algoritmo principal, en el cual se presentan las instrucciones necesarias para emplear la metodología.

Como se puede observar, el Algoritmo 1 agrupa en un solo vector v₀ ÎRⁿ x Zⁿ las cantidades iniciales de energía activa promedio y de intervalo de tiempo de inicio de construcción de cada uno de los proyectos inversión que componen el PPISEP inicial, y establece una matriz Z donde se representan todas las características fijas de cada uno de dichos proyectos de inversión. El vector v₀ representa el vector inicial del algoritmo de optimización ASEM, el cual comienza a generar vectores v hasta determinar un vector final v_f ÎRⁿ x Zⁿ que genera una disminución notoria en el valor de su función objetivo. Este último vector representa la energía activa promedio y el intervalo de tiempo de inicio de construcción de cada uno de los proyectos de inversión que conforman el PPISEP mejorado.

Algoritmo 2. Metodología que permite evaluar la expansión de un Sistema Eléctrico de Potencia (SEP). Este algoritmo constituye la función objetivo del algoritmo de optimización ASEM puesto que permite evaluar un PPISEP.

En este último algoritmo existen dos ciclos iterativos principales, el primero calcula la cantidad de muestras que se deseen y el segundo simula el paso de los intervalos de tiempo durante el Horizonte de Planificación considerado. A continuación se explican los procedimientos que componen este algoritmo.

Mediante este procedimiento se construye una Matriz W en la que se guarda toda la información asociada a cada uno de los proyectos de inversión que componen al PPISEP que desea ser evaluado. El vector nulo de dimensión n es representado mediante la notación 0_n.

Procedimiento 2. Asignación de los valores de los arcos nominales de consumo de la Red R.

De acuerdo a los tipos de clientes t_ih presentes en los Nodos Terminales de la Red R con sus respectivas distribuciones de probabilidad tipo triangular, se determinan los valores de los arcos nominales de consumo para el intervalo de tiempo considerado.

Procedimiento 3. Adición de los costos de proyectos de inversión.

Procedimiento 4. Modificación de los valores de los arcos nominales de generación y de los arcos nominales máximos de transmisión de la Red R debido a la entrada en funcionamiento de los proyectos de inversión.

Procedimiento 5. Determinación de los valores de los arcos nominales de transmisión de la Red R.

Mediante este procedimiento se calculan los valores de los arcos nominales de transmisión de la Red R para el intervalo de tiempo considerado como una disminución aleatoria de los valores de los arcos nominales máximos de transmisión , mediante el uso de una distribución de probabilidad tipo triangular. En los procedimientos siguientes se emplea la notación S(i) para representar los Nodos sucesores al Nodo i y, la notación P(i) para representar los Nodos predecesores al Nodo i (Strang, 2005).

Procedimiento 6. Determinación de los valores de los arcos efectivos x_ij de la Red R. Mediante este procedimiento se resuelve el problema de flujo máximo (también denominado despacho óptimo de energía activa promedio) de la Red R, obteniéndose los valores de todos los arcos efectivos x_ij para el intervalo de tiempo considerado. Dicho problema es resuelto mediante programación lineal puesto que tanto la función objetivo como las restricciones del problema son expresiones lineales.

Procedimiento 7. Contabilización de los ingresos. Mediante el uso de las tarifas k_tipo y de la energía efectivamente entregada a cada uno de los consumidores, este procedimiento permite determinar los ingresos de la empresa eléctrica para el intervalo de tiempo considerado.

Procedimiento 8. Contabilización de los egresos.

Procedimiento 9. Cálculo del balance económico en valor presente. Mediante el método de valor presente (Sullivan et al., 2004), este procedimiento determina el balance económico de la empresa eléctrica para el intervalo de tiempo considerado.

Procedimiento 10. Cálculo de la demanda insatisfecha total en valor presente. Mediante el método de valor presente (Sullivan et al., 2004), este procedimiento determina la demanda insatisfecha total (en términos de dinero) para el intervalo de tiempo considerado.

Procedimiento 11. Cálculo de los valores promedio de balance económico y de demanda insatisfecha.

Procedimiento 12. Cálculo de los valores de las funciones de penalización. Este procedimiento permite evaluar cada una de las funciones de penalización asociadas al PPISEP que desee ser evaluado.

Función 1. Generador de distribución triangular a partir de una uniforme.

IMPLEMENTACIÓN DEL EJEMPLO NUMÉRICO

El ejemplo numérico considerado se basa en un Sistema Eléctrico de Potencia (SEP) cuyo diagrama unifilar se presenta en la figura 1. Como se puede observar en la figura, se tienen dos generadores, once líneas de transmisión y seis cargas en el SEP considerado para el ejemplo numérico. Para el presente ejemplo numérico, los intervalos de tiempo fueron establecidos de una duración igual a dos años.

En la tabla 1 se muestran las capacidades nominales de generación de energía de los generadores G₁ y G₂ y las capacidades nominales máximas de transmisión de energía de las líneas de transmisión L₁ a L₁₁ del SEP considerado. En cuanto a las líneas de transmisión, se estableció que éstas podían tener una capacidad nominal típica o más probable igual al 90% de su capacidad nominal máxima y una capacidad nominal mínima igual al 80% de su capacidad nominal máxima, es decir, mín_falla = 0,8; mod_falla = 0,9; máx_falla = 1. Finalmente, a cada una de las cargas C₁, C₂ y C₃ se le asignó un cliente residencial, mientras que a cada una de las cargas C₄ y C₅ se le asignó un cliente comercial y a la carga C₆ se le asignó un cliente industrial.

Descripción del PPISEP - Para el presente ejemplo numérico se estableció estudiar el SEP antes descrito durante un Horizonte de Planificación Hp de veinte años, compuesto de 10 intervalos de tiempo de 2 años de duración (= 10, t= 2 años).

Dentro de este Horizonte de Planificación, se consideró el siguiente PPISEP inicial de tres proyectos inversión (n3=):

• Proyecto de transmisión N°1: incremento de la capacidad nominal máxima de transmisión de energía activa promedio por intervalo de la línea de transmisión de energía L₇ en 10 unidades de energía. Intervalo de tiempo de inicio de su construcción: N°2 (cuarto año).

• Proyecto de transmisión N°2: construcción de una nueva línea de transmisión de energía que conecte el transformador T₁ con la Carga C₄ y cuya capacidad nominal máxima de transmisión de energía activa promedio por intervalo de tiempo es de 10 unidades de energía. Intervalo de tiempo de inicio de su construcción: N°2 (cuarto año).

• Proyecto de generación N°1: incremento de la capacidad nominal de generación de energía activa promedio por intervalo de tiempo del Generador G₂ en 10 unidades de energía. Intervalo de tiempo de inicio de su construcción: N°2 (cuarto año).

Por otra parte, en la Red R considerada se tienen dos arcos de generación (arcos que unen al Nodo Fuente con los Nodos Iniciales), doce arcos de transmisión (arcos que unen a los Nodos Iniciales con los Nodos Intermedios y arcos que unen a los Nodos Intermedios con los Nodos Terminales) y seis arcos de consumo (arcos que unen a los Nodos Terminales con el Nodo Hundido).

Todos estos arcos presentan una dirección establecida desde el Nodo Fuente hasta el Nodo Hundido. Ahora bien, el arco de la Red R que une al Nodo Intermedio 4 con el Nodo Terminal 10 no representa alguna línea de transmisión en el diagrama unifilar del SEP mostrado en la figura 1. Esto se debe a que dicho arco hace referencia a la nueva línea de transmisión que puede introducir el proyecto de transmisión N°2.

PPISEP inicial en términos matemáticos - Es importante destacar que para definir un proyecto de generación, se seleccionó la variable t_p = 1, mientras que para definir un proyecto de transmisión, se seleccionó la variable t_p = 2. De acuerdo con la descripción expuesta en la definición 3, a continuación se presenta el PPISEP inicial:

Escenarios estudiados - La metodología fue implementada en dos escenarios de consumo pronosticado en alza de energía activa promedio en el tiempo por parte de los clientes presentes en el SEP, los cuales difieren únicamente en los valores de las tasas de interés atractiva e inflación promedio por intervalo de tiempo. El primer escenario toma en consideración una tasa de interés atractiva por intervalo de tiempo () igual a 17% y una tasa de inflación promedio por intervalo de tiempo () igual a 10 %, mientras que en el segundo escenario se considera una tasa de interés atractiva por intervalo de tiempo () igual a 10% y una tasa de inflación promedio por intervalo de tiempo () igual a 17 %. En la tabla 2 se presentan los valores de los parámetros a_tipo, b_tipo, g_tipo, mín_tipo, mod_tipo y máx_tipo para cada tipo de consumidor (residencial, comercial e industrial) y los valores fijos de los parámetros restantes.

El ejemplo numérico fue programado en MATLAB 7 ®, estableciendo un total de 30 muestras (=30) o ejecuciones del modelo de simulación para obtener los valores de la función objetivo del ASEM en cada punto. El número de muestras o ejecuciones del modelo de simulación fue establecido en 30 para obtener una mínima muestra significativa de valores en cada punto de evaluación del ASEM que permitiera obtener la mejora del PPISEP inicial en un tiempo razonable. Por último, para el algoritmo ASEM se consideró como criterio de parada un número máximo de etapas igual a 12, estableciendo de esta manera que los valores que brinda el ASEM de tiempo de inicio de construcción y de cantidad de energía activa promedio que introduce cada proyecto de inversión del PPISEP considerado constituyen unos valores mejores a los iniciales mas no los valores óptimos. A medida que se considere un número mayor de etapas del ASEM, los valores obtenidos estarán más cercanos a los valores óptimos.

RESULTADOS DEL EJEMPLO NUMÉRICO

• Escenario N°1. Consumo en alza. Interés mayor que inflación.

El PPISEP mejorado arrojado por el ASEM en la etapa 12 luego de 1024 evaluaciones de la función objetivo es el siguiente:

Al evaluar el rendimiento del PPISEP mejorado se obtuvo un intervalo de confianza al 95% para la media de la demanda insatisfecha total (en términos de dinero) μ_d igual a: -18017,25506 ≤ μ_d ≤ -3782,208917 y un intervalo de confianza al 95% para la media del balance económico μ_b igual a: 385369,923 ≤ μ_b ≤ 386862,6541 Las figuras 5 y 6 presentan los histogramas de este escenario con el PPISEP mejorado.

• Escenario N°2. Consumo en alza. Inflación mayor que interés.

Al evaluar el rendimiento del PPISEP inicial se obtuvo un intervalo de confianza al 95% para la media de la demanda insatisfecha total (en términos de dinero) μ_d igual a: -2843822,939 ≤ μ_d ≤ -2629141,657 y un intervalo de confianza al 95% para la media del balance económico μ_b igual a: 1029475,9 ≤ μ_b ≤ 1032349,327. Las figuras 7 y 8 presentan los histogramas de este escenario con el PPISEP inicial.

El PPISEP mejorado arrojado por el ASEM en la etapa 12, luego de 1424 evaluaciones de la función objetivo, es el siguiente:

Al evaluar el rendimiento del PPISEP mejorado se obtuvo un intervalo de confianza al 95% para la media del balance económico μ_b igual a: 837336,4542 ≤ μ_b ≤ 840096,5465 y una media de demanda insatisfecha total (en términos de dinero) μ_d igual a cero. Las figuras 9 y 10 presentan los histogramas de este escenario con el PPISEP mejorado.

ANÁLISIS DE RESULTADOS

CONCLUSIONES

La propuesta metodológica desarrollada constituye un aporte inédito dentro del campo de la planificación de la expansión de los SEP. Partiendo de un PPISEP inicial y haciendo uso de la técnica de optimización ASEM, esta metodología identifica unos valores mejores de tiempo de inicio de construcción (variable entera) y de energía activa promedio introducida al SEP (variable real) de cada uno de los proyectos de inversión que componen el PPISEP. Éstos constituyen una disminución del consumo pronosticado no cubierto de energía activa promedio y un incremento del balance económico de la empresa eléctrica. Por otra parte, mediante el uso de la simulación por Monte Carlo, la metodología diseñada considera las incertidumbres aleatorias asociadas al consumo energético y las fallas en las líneas de transmisión del SEP. Finalmente, la virtud más importante de la propuesta metodológica desarrollada es la de permitir considerar diferentes topologías de SEP y distintos tipos de PPISEP, con diferente cantidad y tipos de proyectos de inversión. Aun cuando en el presente artículo se muestra un SEP simplificado para ilustrar la metodología desarrollada, el modelo computacional expuesto permite considerar diversas topologías radiales de SEP. Es relevante destacar que aun cuando no se toma en cuenta la ampliación de la capacidad energética de los transformadores en la metodología descrita, los mismos pueden ser considerados como arcos de transmisión cuyos valores representan las capacidades de transmisión de energía activa promedio de los transformadores.

ANEXO 1. NOMENCLATURA

m, t, r: contador de muestras, de intervalos de tiempo y de proyectos de inversión respectivamente;

i: contador de los Nodos de Cola de los arcos de la Red R;

j: contador de Nodos de Cabeza de los arcos de la Red R;

k: contador de Nodos de Cola y Nodos de Cabeza de los arcos de la Red R;

b_m, d_m: balance económico de la empresa eléctrica y demanda insatisfecha total (en términos de dinero) para el Horizonte de Planificación Hp en valor presente de la muestra m;

: sumador de balance económico de la empresa eléctrica y sumador de demanda insatisfecha total (en términos de dinero) en valor presente de cada muestra m;

: arco nominal de generación de la Red R cuyo Nodo de Cola es el Nodo 1 (Nodo Fuente) y cuyo Nodo de Cabeza es el j-ésimo Nodo Inicial, es decir, tiene una dirección asignada hacia el Nodo Inicial, y su valor representa la capacidad nominal de generación de energía activa promedio en el intervalo de tiempo t del generador asociado al j-ésimo Nodo Inicial,

: arco nominal máximo de transmisión de la Red R cuyo Nodo de Cola es el i-ésimo Nodo Inicial o el i-ésimo Nodo Intermedio y cuyo Nodo de Cabeza es el j-ésimo Nodo Intermedio o el j-ésimo Nodo Terminal, y su valor representa la capacidad nominal máxima de transmisión de energía activa promedio en el intervalo de tiempo t de la línea de transmisión ubicada entre dichos Nodos (sin considerar las fallas aleatorias asociadas a dicha línea en el intervalo de tiempo t), "i Î {2, 3,...,N_g + N_t + 1}, "j Î {N_g +  2, N_g + 3,...,N_g + N_t + N_c + 1}. Si alguno de estos arcos es añadido a la Red R por parte de algún proyecto de transmisión, su valor debe ser inicialmente igual a cero hasta que dicho proyecto de transmisión comience a ejecutarse, momento en el cual pasa a ser igual a la energía activa promedio que introduce el proyecto de transmisión.

: arco nominal de transmisión de la Red R cuyo Nodo de Cola es el i-ésimo Nodo Inicial o el i-ésimo Nodo Intermedio y cuyo Nodo de Cabeza es el j-ésimo Nodo Intermedio o el j-ésimo Nodo Terminal, y su valor representa la capacidad nominal de transmisión de energía activa promedio en el intervalo de tiempo t de la línea de transmisión ubicada entre dichos Nodos (considerando las fallas aleatorias asociadas a dicha línea en el intervalo de tiempo t), "i Î {2, 3,...,N_g + N_t + 1}, "j Î {N_g +  2, N_g + 3,...,N_g + N_t + N_c + 1};

: arco nominal de consumo de la Red R cuyo Nodo de Cola es el i-ésimo Nodo Terminal y cuyo Nodo de Cabeza es el Nodo h = N_g + N_t + N_c + 2 (Nodo Hundido), y su valor representa el consumo pronosticado de energía activa promedio en el intervalo de tiempo t de la carga asociada al i-ésimo Nodo Terminal, "i Î {N_g + N_t + 2, N_g + N_t + 3,...,N_g + N_t + N_c + 1}.

x_ij: arco efectivo de generación, transmisión o consumo de la Red R cuyo valor representa el mayor flujo de energía activa promedio (energía efectiva) entre el i-ésimo Nodo y el j-ésimo Nodo en cada intervalo de tiempo t, "i Î {1, 2,...,N_g + N_t + N_c + 1}, "j Î {2, 3,...,N_g + N_t + N_c + 2};

i_tipo: sumador de los ingresos obtenidos por la empresa eléctrica, de acuerdo al tipo de consumidor, en cada intervalo de tiempo t;

i_g: ingresos totales obtenidos por la empresa eléctrica en cada intervalo de tiempo t;

: sumador de los costos variables (costos que varían de acuerdo a la energía activa promedio generada por intervalo de tiempo) que tiene la empresa eléctrica en cada intervalo de tiempo t;

e_c: sumador de los egresos de la empresa eléctrica en cada intervalo de tiempo t. Dichos egresos abarcan los costos fijos y los costos variables en el intervalo de tiempo t, además de los costos de los proyectos de inversión, en el caso de que éstos inicien su construcción en el intervalo de tiempo t;

b_p, d_p: balance económico de la empresa eléctrica y demanda insatisfecha total (en términos de dinero) en cada intervalo de tiempo t;

: sumador de balance económico de la empresa eléctrica y sumador de demanda insatisfecha total (en términos de dinero) para el Horizonte de Planificación Hp en valor presente de cada muestra m;

B_p: promedio de los valores de balance económico de la empresa eléctrica para el Horizonte de Planificación Hp en valor presente obtenidos en cada una de las muestras m;

D_p: promedio de los valores de demanda insatisfecha total (en términos de dinero) para el Horizonte de Planificación Hp en valor presente obtenidos en cada una de las muestras m;

: sumadores asociados a las funciones de penalización;

: funciones de penalización.

N_g, N_t, N_c: número de Nodos Iniciales, Nodos Intermedios y Nodos Terminales de la Red R respectivamente. Para el caso de N_g y N_t, dicho número considera la posible adición de generadores o líneas de transmisión al SEP por parte de algún proyecto de inversión;

G_1j: arco inicial de generación de la Red R cuyo Nodo de Cola es el Nodo 1 (Nodo Fuente) y cuyo Nodo de Cabeza es el j-ésimo Nodo Inicial, y su valor representa la capacidad nominal inicial de generación de energía activa promedio en el intervalo de tiempo t del generador asociado al j-ésimo Nodo Inicial, "j Î {2, 3,...,N_g + 1}. Si alguno de estos arcos es añadido a la Red R por parte de algún proyecto de generación, su valor debe ser igual a cero;

T_ij: arco inicial de transmisión de la Red R cuyo Nodo de Cola es el i-ésimo Nodo Inicial o el i-ésimo Nodo Intermedio y cuyo Nodo de Cabeza es el j-ésimo Nodo Intermedio o el j-ésimo Nodo Terminal, y su valor representa la capacidad nominal inicial máxima de transmisión de energía activa promedio en el intervalo de tiempo t de la línea de transmisión ubicada entre dichos Nodos, "i Î {2, 3,...,N_g + N_t + 1}, "j Î {N_g +  2, N_g + 3,...,N_g + N_t + N_c + 1}. Si alguno de estos arcos es añadido a la Red R por parte de algún proyecto de transmisión, su valor debe ser igual a cero;

: número total de muestras del modelo de simulación;

: número de intervalos de tiempo que conforman el Horizonte de Planificación Hp;

n: número total de proyectos que componen el PPISEP;

t_ih: tipo de cliente asociado a cada Nodo Terminal i de la Red R, "i Î {N_g + N_t + 2, N_g + N_t + 3,...,N_g + N_t + N_c + 1};

a_tipo, β_tipo, g_tipo: parámetros que definen la tendencia del consumo pronosticado de un tipo de consumidor;

mín_tipo, mod_tipo, máx_tipo: parámetros que definen la distribución de probabilidad tipo triangular (valor mínimo, moda y valor máximo) asociada al consumo pronosticado de un tipo de consumidor;

mín_falla, mod_falla, máx_falla: parámetros que definen la distribución de probabilidad tipo triangular (valor mínimo, moda y valor máximo) asociada al cálculo de las capacidades nominales de las líneas de transmisión de energía (valores de los arcos );

k_tipo: tarifa a la cual se factura la energía efectivamente entregada a los consumidores (valores de los arcos efectivos de consumo) según el tipo de consumidor;

c_u: costo por unidad de energía activa promedio generada (valores de los arcos efectivos de generación) que tiene la empresa eléctrica;

c_f: costo fijo de cada intervalo de tiempo que tiene la empresa eléctrica;

l_tipo: factor de penalización de demanda insatisfecha según el tipo de consumidor;

,: tasa de interés mínima atractiva y tasa de inflación por intervalo de tiempo;

REFERENCIAS

1. Azadivar F. (1999). Simulation Optimization Methodologies. Proceedings of the 1999 Winter Simulation Conference, pp. 93 - 100.        [ Links ]

2. Brea E. (2013). Una extensión del método de Nelder Mead a problemas de optimización no lineales enteros mixtos. Revista Internacional de Métodos Numéricos para Cálculo y Diseño en Ingeniería, vol. 29, no. 3, July- September 2013, Pages 163-174, ISSN_0213-1315, http://dx.doi.org/10.1016/j.rimni.2013.06.005        [ Links ]

3. Carson Y. & Maria A. (1997). Simulation Optimization: Methods and Applications. Proceedings of the 1997 Winter Simulation Conference, pp. 118 - 126.        [ Links ]

4. Hreinsson E. & Jónasson K. (2003). Monte Carlo based risk analysis in Hydroelectric Power System Expansion Planning in the presence of uncertainty in project cost and capacity. 38th International Universities Power Engineering Conference (UPEC), University of Thessaloniki, septiembre 2003.        [ Links ]

5. Nelder J. & Mead R. (1965). A simplex method for function minimization. The Computer Journal, vol. 7, pp. 308 - 313.        [ Links ]

6. Strang G. (2005). Linear Algebra and its Applications. India: 4ta ed. Brooks Cole. 128 - 140.        [ Links ]

7. Su C. & Lii G. (1997). Power System Capacity Expansion Planning Using Monte Carlo Simulation and Genetic Algorithm. MODSIM 1997 International Congress on Modelling and Simulation. Modelling and Simulation Society of Australia and New Zealand, pp. 1450 - 1455, diciembre 1997.        [ Links ]

8. Sullivan W., Wicks E., Luxhoj J. (2004). Ingeniería Económica de DeGarmo (12). México: Editorial Pearson Educación. 70 - 148.        [ Links ] ]]> 1 1999 93 - 100 2 2013 20 13 29 3 3 163-174 3 1997 118 - 126 4 2003 septiembre 2003 5 1965 7 308 - 313 6 2005 4ta 128 - 140 7 1997 1450 - 1455 8 2004 70 - 148