INTRODUCCION
En Cuba existen 125 instalaciones hidroeléctricas a filo de agua. De ellas, 117 están diseñadas para que sólo funcionen de manera aislada.
Se han identificado proyectos a ser evaluados o ejecutados en el corto, mediano y largo plazo de construcción o rehabilitación de seis centrales a filo de agua de baja potencia para ser conectadas al sistema eléctrico nacional (SEN). (Ver tabla 1). El 50% son de baja carga y todas son de baja potencia.
Instalación | Carga de diseño (m3/s) | Flujo de diseño (m) | Potencia (kW) | Tipo | Acción | Plazo |
---|---|---|---|---|---|---|
Guaso | 180,0 | 1,300 | 1 870 | Agua fluente | Rehabilitación | Corto |
Guayabo | 211,5 | 0,116 | 196 | Agua fluente | Rehabilitación | Mediano/Largo |
Canal Yara | 4,9 | 6,500 | 256 | Canal | Construcción | Mediano/Largo |
Barranca | 7,1 | 3,750 | 214 | Canal | Rehabilitación | Mediano/Largo |
Guane | 10,0 | 0,650 | 52 | Agua fluente | Rehabilitación | Mediano/Largo |
Gota Blanca | 80,0 | 0,783 | 500 | Canal (trasvase) | Construcción | Largo Plazo |
Existe un potencial de nuevos proyectos de centrales a filo de agua por explorar, que se conectarían al SEN, fundamentalmente a partir de los trasvases construidos o en construcción. Se estima que estos nuevos desarrollos también serían de baja carga y potencia.
Los sitios de baja carga y potencia requieren de la mayor precisión en la determinación de la carga y el gasto de diseño para lograr que sean económicamente factibles, pues el costo de inversión por kW, de acuerdo a la evidencia nacional e internacional, aumenta en la medida que la potencia disminuye. La optimización de los proyectos de inversión, entonces, constituye una necesidad para que estos proyectos sean económicamente viables.
EL ENFOQUE TRADICIONAL PARA DETERMINAR EL FLUJO DE DISEÑO
El Dr. García Faure en su tesis de doctorado, recoge los métodos más utilizados para determinar el caudal de diseño de una central para los proyectos a filo de agua, y afirma que “…en Cuba se determina el caudal de proyecto (…) considerando la curva de caudales clasificados para el año seco con frecuencia de ocurrencia de alrededor del 85 % y de esa curva se obtiene el gasto que puede ser mantenido durante 310 días del año, o sea el 85% del tiempo”. La utilización de este método puede ser factible bajo determinadas condiciones específicas, sobre todo para algunas centrales aisladas, pero de ninguna manera puede generalizarse, so riesgo de subdimensionar los proyectos (García 2005).
El uso de esta regla se corroboró recientemente en la Solución Conceptual del proyecto de mini hidroeléctrica de Guane, con conexión al SEN (PROAGUA 2015). En el caso de la Solución Conceptual de Guayabo, igualmente con conexión al SEN, se utilizó la curva de caudales clasificados con frecuencia de ocurrencia del 85 % y como flujo de diseño el gasto promedio del periodo húmedo menos el flujo ecológico (RAUDAL 2018). Las curvas de caudales de años secos y húmedos son necesarias, pero a los efectos de revisar parámetros técnicos del sitio ante condiciones hidrológicas cambiantes o para realizar análisis de sensibilidad, entre otras aplicaciones.
Para el caso de Barranca, también a ser conectada al SEN, se utilizó la curva de caudales (medios) clasificados y como flujo de diseño el gasto mayor de los que más se repiten como caudal de diseño, obtenidos del histograma de caudales. (RAUDAL 2017). En estos casos, la producción de energía o el VAN son resultados pasivos de parámetros estadísticos fijados a priori.
En la práctica internacional actual se parte de los caudales promedio, y se recomienda que el diseño del flujo, sea ante todo un problema técnico-económico. Su fijación debe ser producto de un proceso iterativo de optimización (García 2005) (ESHA 2004). Las funciones objetivo más frecuentes son la energía (Mohsin et al. 2015) o el VAN (Rajšl et al. 2015).
EL ENFOQUE TRADICIONAL PARA DETERMINAR LA CARGA DE DISEÑO
El Dr. Pérez Franco define la carga de diseño como la carga… “por encima y por debajo de la cual, el promedio anual de generación por encima y por debajo es aproximadamente igual.” Recomienda utilizar como carga de diseño “la media aritmética de las cargas netas disponibles en cada periodo de tiempo” para estudios de pre-factibilidad y factibilidad, aunque precisa que esta carga obtenida puede variar en función de la turbina a utilizar (Pérez 1999). Bajo esta recomendación, compartida por numerosos autores, la carga de diseño (Hd) se calcularía de la siguiente forma:
donde:
Δt i : |
Periodo de tiempo que dura cada carga neta según dato i |
Hn i : |
Carga neta según el dato i |
n: |
número de datos de carga neta - periodo de tiempo |
Como que en los proyectos a filo de agua no hay acumulación de agua o esta no es significativa, no deben existir grandes variaciones de la carga bruta y, por lo tanto, la media aritmética ponderada puede ser una buena aproximación a la carga de diseño: cargas netas máximas y mínimas del sitio no deben estar fuera del rango permisible de operación de la turbina (ver tabla 2 y figura 1).
Turbina | Fuente | ||
---|---|---|---|
Kaplan | 1,20 | 0,70 | Promedio proveedores europeos y chinos en licitación en 2018 |
Francis | 1,20 | 0,70 | Promedio proveedores europeos y chinos en licitación en 2018 |
Pelton | 1,10 | 0,75 | Del software norteamericano ORNL-HEEA Tool (ORNL 2013) |
Propela | 1,10 | 0,80 | Del software norteamericano ORNL-HEEA Tool (ORNL 2013) |
Flujo cruzado | 1,25 | 0,55 | Área de operación de turbina del fabricante Ossberger |
Tornillo | 1,10 | 0,75 | Promedio de ofertas de varios proveedores europeos |
Leyenda:
lhmáx: Límite relativo superior de carga de la turbina
lhmín: Límite relativo inferior de carga de la turbina
En la práctica actual en Cuba, de acuerdo con las soluciones conceptuales de sitios a filo de agua revisadas, se han utilizado otras dos variantes:
En el caso de Barranca, a partir del rango de gastos con más posibilidades de turbinar, se evaluaron las pérdidas locales y de conducción de la instalación para esos caudales (condiciones de entrada y salida, rejillas, compuerta, conductora, contracciones, ampliaciones y cambio de dirección), y a la carga bruta correspondiente al flujo de diseño, se le restaron las pérdidas asociadas a ese flujo, y así se determinó la carga de diseño (RAUDAL 2017).
En el caso de Guayabo (RAUDAL 2018), se utilizó la ecuación:
En la práctica internacional también se utiliza la curva de duración de la carga para la determinación de la carga de diseño, en particular para los sitios de baja carga, en los cuales pequeñas variaciones absolutas del salto hidráulico pueden provocar significativos cambios en la energía a generar (AECOM 2011) (ORNL 2013).
Sin embargo, en trabajos recientes se también se recomienda que el cálculo de la carga de diseño debe ser producto de un proceso iterativo de optimización (Rajšl et al. 2015).
MÉTODOS ALTERNATIVOS PARA LA DETERMINACIÓN DEL FLUJO Y LA CARGA DE DISEÑO PRELIMINAR
El modelo matemático para optimización tiene tres componentes: la función objetivo, que es la que expresa lo que debe ser optimizado; el conjunto de variables de decisión del problema cuyos valores deben determinarse para que la función objetivo alcanza su valor óptimo; y las restricciones que acotan los valores que pueden tomar las variables.
Los avances en los medios de cómputo y de las comunicaciones desde fines del siglo pasado, han permitido la “democratización” de la optimización, de manera que puede ser realizada en prácticamente cualquier sector en tiempos relativamente breves y con bajos costos.
La explicación de los métodos de optimización para los proyectos hidroenergéticos a filo de agua se realizará a través de la aplicación informática “Optim-Filo agua 8.1 xlms” desarrollada por profesionales de la Empresa de Hidroenergía.
El Dr. García Faure, en su pionera incursión en el tema de la optimización de los proyectos de centrales hidroeléctricas en Cuba, define como la función objetivo la maximización de la relación Ingreso/Costo (García 2015).
En el presente trabajo se aborda la optimización a partir de dos modelos:
PLANTEAMIENTO DE LA OPTIMIZACION DE LA GENERACION DE ENERGIA
En determinadas circunstancias, lo que se desea es maximizar la generación eléctrica de un proyecto. Los indicadores económico-financieros serían un resultado de la energía obtenida como máxima.
El universo de información con los cuales se resolverá el problema son los datos i de la curva de caudales clasificados turbinables y las cargas netas correspondientes a cada gasto. El planteamiento muy simplificado del problema es el siguiente (ver tabla 3).
Función objetivo | |
| |
Variables de decisión: | |
|
No se aprovecha la carga |
|
No se turbina si el flujo |
|
La eficiencia de la turbina depende de |
|
Según el tipo de turbina existe una función para la eficiencia para |
|
La eficiencia del generador depende del cociente de la potencia |
|
Se asumen pérdidas de transformación de un 3%. |
Si hay inversor y reductor se asumen entre un 6-8%. | |
Se asume una disponibilidad técnica de alrededor del 96%. | |
Los límites operativos de carga dependen del tipo de turbina (Ver tabla 4) | |
Los límites operativos de flujo dependen del tipo de turbina. (Ver tabla 5). | |
- | |
- | |
|
- |
|
- |
Donde:
E: Generación de energía anual; E i : Generación de energía en el momento i en MWh
P i : Potencia total en el momento i en kW; t i : Tiempo en horas en el momento i
Qi: Flujo en el momento i en m3/s; Hd; Carga de diseño en m
Qd: Flujo de diseño en m3/s; Hn i : Carga neta en el momento i en m
ηt i : Efic. de la turbina en el momento i en %; P i : Potencia en el momento i en kW
T: Tipo de turbina; ηg i : Eficiencia del generador en el momento i en %
Pt i : Pot. de la turbina en el momento i en kW; Pt: Pot. máxima de la turbina en kW
ηf: Pérdidas de transformación en %; p: Pérdidas varias en % (inversor/reductor)
d: Disponibilidad técnica en % del tiempo; Qmáx: Flujo máx. de la curva de caudales clasificados en m3/s
Turbina | Fuente | ||
---|---|---|---|
Kaplan | El flujo máximo que pasa por la turbina es |
0,30 | Promedio proveedores europeos y chinos en licitación en 2018 |
Francis | 0,45 | Promedio proveedores europeos y chinos en licitación en 2018 | |
Pelton | 0,15 | Promedio proveedores europeos y chinos en licitación en 2018 | |
Propela | 0,35 | Del |
|
Flujo cruzado | 0,20 | Información de fabricantes |
|
Tornillo | 0,15 | Promedio de ofertas de varios proveedores europeos | |
Leyenda |
|
Año 2018 | Kaplan/Francis /Pelton | Propela (Proporción muestra de |
Flujo Cruzado (Proporción muestra de |
Tornillo | |
---|---|---|---|---|---|
| |||||
|
|||||
|
|
|
|||
|
|
|
|||
|
|
|
|||
|
|||||
|
|||||
|
|||||
PLANTEAMIENTO DE LA OPTIMIZACION DEL VALOR ACTUAL NETO
En planteamiento en forma también simplificada del problema, agrega varias restricciones al anterior (ver tabla 6).
Función objetivo | |
| |
Restricciones adicionales | Comentarios |
|
El costo del equipamiento y de la ingeniería básica asociada en moneda extranjera es una función de la potencia de las turbinas (o de la capacidad instalada) y del tipo de turbina. Es necesario tener las curvas de costos (ver tabla 6). |
|
Costo del equipamiento y de la ingeniería básica asociada en pesos. |
Los costos de inversión en obras y civiles y los gastos previos son específicos del sitio. Se obtienen de los presupuestos de los proyectos. | |
Se calcula como un porcentaje de |
|
La tarifa del kWh generado es $0.27 para estudios de factibilidad. | |
Los gastos anuales de salarios y seguridad social | |
Los gastos de mantenimiento se calculan enter $0.02 y $0,04/kWh . | |
|
Impuestos |
|
La depreciación en el año |
Se asume un 35% de la utilidad antes de impuestos | |
Usualmente se asume una tasa entre el 3% y el 15%. |
Donde:
K d : Costo del equipamiento y de la ingeniería básica (importado) en moneda extranjera
K: Costo del equipamiento y de la ingeniería básica (importado) en pesos
e: Tasa de cambio en moneda extranjera/pesos
g: Margen comercial de la empresa importadora
I: El resto de los costos de inversión en pesos
k: Capital de trabajo en pesos
n: Vida útil de la instalación en años
E i : Generación de energía en el momento i en MWh
p: Tarifa por kWh en pesos
s: Gastos anuales de salario y seguridad social en pesos
m: Gastos de mantenimiento en pesos
d i : Depreciación en el año j en pesos
t: Tasa impositiva en %
r: Tasa de descuento en %
RESULTADOS DE LA APLICACION DEL MÉTODO TRADICIONAL
Una condición importante para aplicar los métodos de optimización es programar el modelo (la función objetivo, las variables y las restricciones), en Excel, Matlab u otro software. La realización de cálculos manualmente puede ser tedioso, ocupar una enorme cantidad de tiempo y existe, además, mayor posibilidad de errores de cálculo. Con una aplicación como “Optim-Filo agua 8.1 xlms”, la tarea se simplifica.
A continuación se expondrá el método de solución mediante iteraciones y gráficos. Para una descripción más ilustrativa de los métodos de optimización mediante iteraciones y gráficos, y mediante el uso del Gradiente Reducido Generalizado, será utilizada para la información del proyecto de Barranca (RAUDAL 2018): la curva de caudales clasificados y las cargas netas asociadas a cada valor de flujo, así como las restricciones específicas relativas al proyecto. Asimismo, sólo se expondrá la aplicación de ambos métodos para la optimización del VAN, pues en esencia son similares. El modelo para la maximización de la energía sólo se diferencia en la función objetivo y en que utiliza menos restricciones.
En la figura 2, tomada de una de las hojas de la aplicación “Optim-Filo agua 8.1 xlms”, son mostrados los resultados de la generación y el VAN, a partir de la carga y el flujo de diseño propuestos en la Solución Conceptual (Hd=7,13 m, Qd=3,75 m3/s). También se muestran en la tabla 7, valores estadísticos mínimo, medio y máximo del flujo y la carga neta de dicho sitio.
SOLUCION MEDIANTE ITERACIONES Y GRÁFICOS
Este método es sumamente sencillo si se dispone de medios de cómputo:
a) De la tabla 7 de toman los flujos y cargas netas mínimos y máximos del sitio.
b) Se determinan a priori para cuántos posibles valores de flujo y carga neta de diseño se van a evaluar el VAN. En el ejemplo, son 21 de flujo de diseño y 5 de carga de diseño.
c) Se determinan ΔQ y ΔHn. Ellos serán los incrementos que se adicionarán a los flujos y cargas mínimos del sitio.
d) Se organizan en la tabla los flujos y las cargas sumando ΔQ y ΔHn a los flujos y cargas precedentes.
e) Se determina el VAN para cada una de las combinaciones de flujo y carga de diseño, con ello queda rellena la tabla 8. El máximo VAN es producto de la combinación Hd= 7,01 y Qd= 5,25 m3/s (valores en negritas en la tabla 8)
f) Con los datos de dicha tabla se crea un gráfico 3D. (ver figura 3)
g) Se observan en la parte superior del gráfico los rangos de las combinaciones Hd y Qd que maximizan el VAN.
h) Se repite el proceso desde el acápite a), pero se toman nuevos flujos y cargas de diseño máximos y mínimos para cuyo rango se determinará el VAN (se va cerrando la zona de optimización), se rellena de nuevo la tabla y se construye el nuevo gráfico. (ver figura 4). Obsérvese que con el segundo gráfico se ha obtenido una combinación Hd/Qd que aumenta el VAN. Este proceso se repite hasta que no existan cambios en el VAN.
SOLUCION MEDIANTE EL METODO DE GRG
Otra forma se solución es la siguiente: una vez expresado el problema de optimización en hojas de cálculo, se puede aplicar el método GRG no lineal. Debe ser aplicado varias veces para tener la mayor certeza de que ha obtenido un óptimo global. (Ver figuras 5 y 6)
COMPARACION DE LOS DIFERENTES MÉTODOS
A continuación se muestran los resultados de aplicar el método tradicional y los métodos de optimización a los datos del proyecto Barranca (Ver tabla 9). También se realiza una evaluación cualitativa de los métodos. (Ver tabla 10).
Tasa de descuento 10% / Proveedor A/B Proy-Eur. | Optimización de energía | ||||
---|---|---|---|---|---|
Indicadores | Método tradicional | Iteraciones y gráficos | GRG | % | |
A | B | C | D | C/B - 1 | D/B - 1 |
7,13 | 7,25 | 7,25 | 1,7 | 1,7 | |
|
3,75 | 5,70 | 5,70 | 52.0 | 52.0 |
Energía (MWh) | 1 402 | 1 604 | 1 604 | 14,4 | 14,4 |
Capacidad instalada (kW) | 213 | 331 | 331 | 55,4 | 55,4 |
VAN |
852 631 | 1 014 078 | 1 014 078 | 18,9 | 18,9 |
Tasa de descuento 10% / Proveedor A/B Proy-Eur. | Optimización del VAN | ||||
Indicadores | Método tradicional | Iteraciones y gráficos | GRG | % | |
A | B | C | D | C/B - 1 | D/B - 1 |
7,13 | 7,05 | 7,02 | -1,1 | -1,5 | |
|
3,75 | 5,40 | 5,40 | 44,0 | 44,0 |
Energía (MWh) | 1 402 | 1 595 | 1 594 | 13,8 | 13,9 |
Capacidad instalada (kW) | 213 | 304 | 303 | 42,7 | 42,3 |
VAN |
852 631 | 1 032 125 | 1 033 309 | 21,1 | 21,2 |
Criterio | Método tradicional | Métodos que optimizan la producción de energía o el VAN | |
---|---|---|---|
Iteraciones y gráficos | Gradiente Reducido Generalizado | ||
Forma de cálculo | Valores estadísticos | Mediante iteraciones sucesivas, reduciendo el rango de las pruebas de |
Aplicando el método de GRG de la herramienta Solver (Excel) |
Precisión | Exacto, pues son valores estadísticos | Con este método no se garantiza el óptimo global, pues las variables sólo toman valores discretos predeterminados. No obstante, cambiando rangos se puede tener la certeza de que se ha logrado o se está cerca del óptimo global. Los gráficos nos dan orientaciones para determinar los rangos. | Con el método GRG no se garantiza un óptimo global. Hay que ejecutar varias veces para tener certeza de que se ha obtenido el máximo global o los valores obtenidos están muy cercanos. |
Relación entre |
Tanto la energía como el VAN son resultados pasivos | Se obtienen los valores máximos o muy cercanos de generación de energía o VAN, según se desee. | |
Potencia resultante | La potencia es un resultado pasivo. | La potencia puede ser mayor o menor con relación a la obtenida a partir del método tradicional. |
CONCLUSIONES
Como se ha probado con este trabajo:
El método tradicional de determinación de la carga y el flujo de diseño para las instalaciones a filo de agua no garantizan que se optimicen los proyectos.
La carga y el flujo de diseño determinados de manera preliminar por los dos métodos de optimización descritos aumentan la potencial generación de energía o el VAN, con relación a los que se obtendrían a partir del método tradicional basado en la estadística,
los métodos aquí propuestos, de ser utilizados por los proyectistas, coadyuvarán a un cálculo más preciso de este importante parámetro.
Los métodos de optimización mediante iteraciones y gráficos o aplicando el método del Gradiente Reducido Generalizado son complementarios, es decir, es recomendable que ambos sean utilizados: el segundo es más preciso, pero el primero ayuda mediante las representaciones gráficas a identificar la zona de posible máximo y permite corrobar los parámetros de diseño obtenidos por el primero.
Se recomienda, además, que las instituciones nacionales de proyectos, universidades, el Instituto Nacional de Recursos Hidráulicos y la Unión Eléctrica, revisen y posiblemente mejoren y mantengan actualizada la herramienta “Optim-Filo agua 8.1 xlms”, desarrollado en la Empresa de Hidroenergía, como un modelo de optimización para la proyección futura de centrales hidroeléctricas más eficientes.
Finalmente, es necesario que se actualice el levantamiento hidroenergético de Cuba, y que se incluyan en dicho trabajo aquellos sitios potenciales que se deriven de los trasvases en proyecto o construcción.