HM 000772
HM 000772
HM 000772
Diciembre, 2020
ISSN 2602-8042 impreso / 2631-2522 digital
Colectores solares de
placa plana en Uruguay.
Eficiencia en la producción
de energía eléctrica,
Bolivia.
Instalación solar de un
proceso productivo,
Argentina.
Permeabilidad al aire
de edificios en Uruguay.
Catalizadores para
producción de gas natural.
Residuos eléctricos
y electrónicos.
Cadena productiva de la
energía y emisión de GEI,
Argentina.
COMITÉ AD-HONOREM
Andrés Romero C.
Pontificia Universidad Católica de Chile.
Leonardo Beltrán.
Institute of the Americas. México.
Manlio Coviello.
Pontificia Universidad Católica de Chile.
Mauricio Medinaceli.
Investigador independiente. Bolivia.
COORDINADORES DE LA EDICIÓN
DIRECTOR GENERAL
Alfonso Blanco
DIRECTORES EJECUTIVOS
Pablo Garcés
Marcelo Vega
COORDINADORA DE PRODUCCIÓN
Blanca Guanocunga. Organización Latinoamericana
de Energía (OLADE).
COLABORADORES
Henry Milton Espada Romero. Rodrigo Alonso Suárez. Universidad de la República (UdelaR).
Universidad Pública de El Alto. Bolivia. Facultad de Ingeniería. Uruguay.
© Copyright Organización Latinoamericana de Energía Diseño de la portada y contraportada Ana María Arroyo.
(OLADE) 2020. Todos los derechos reservados.
NOTA DE RESPONSABILIDAD DE CONTENIDO
ISSN: 2602-8042 (Impresa) Las ideas expresadas en este documento son responsabilidad
ISSN: 2631-2522 (Electrónica) de los autores y no comprometen a las organizaciones
mencionadas.
Dirección: Av. Mariscal Antonio José de Sucre N58-63 y
Fernández Salvador.
Quito - Ecuador
SARS-COV2
2020
“ Según estimaciones
del FMI, América Latina
2020 ha sido un año complejo en el que debemos
lamentar las aproximadamente 1.4 millones de
víctimas del SARS-COV2 a nivel mundial, la crisis
sanitaria derivada de la transmisión acelerada
y El Caribe ha sido una de del virus ha impactado de manera relevante el
comportamiento social y la situación económica
las regiones más afectadas a nivel mundial.
”
7.9% en el PIB nominal 2020 respecto a 2019.
2020 respecto a 2019. Recuperar el nivel precrisis tomará a la región
varios años.
El sector energético se ha visto sometido a El sector energético apalancará la recuperación
fuertes tensiones, por un lado, debiendo apa- económica de la región. Es imprescindible acelerar
lancar el suministro eléctrico desplazado hacia las transiciones energéticas y la innovación en
el sector residencial y servicios de sanidad, el sector, la digitalización en la cadena eléctrica,
desde la demanda industrial, transporte y co- la aceleración del recambio de combustibles en
mercial drásticamente contraída y, por otro, el trasporte, el acceso universal más necesario
las medidas de restricciones de movilidad en esta crisis sanitaria, son ejes necesarios de
provocaron una dramática disminución de la impulsar en esta coyuntura. Estos temas fueron
demanda de combustibles, aumentando los tratados extensamente en la reciente edición de
stocks de crudo que presionaron los precios la Semana de la Energía, evento desarrollado
a la baja pronunciadamente, obligando a la de manera virtual del 16 al 25 de noviembre
OPEP, conjuntamente con productores no OPEP, que se ha convertido en el principal espacio de
a tomar medidas de recorte de la oferta para encuentro de los actores energéticos en la región.
intentar revertir el desplome de los precios. El material grabado del evento se encuentra en el
Éstas entre otras consecuencias de las medidas sitio web: https://semanadelaenergia.olade.org/
adoptadas para intentar la contención de la
transmisión de la enfermedad. El futuro inmediato está lleno de desafíos para
el sector, por lo que este espacio de divulgación
OLADE, en su panorama energético publicado mantiene abiertas sus puertas para apoyar a la
recientemente1, estima una contracción de la comunidad energética en la tarea de análisis y
demanda total de energía cercana al 8% res- búsqueda de soluciones para afrontar el período
pecto a los niveles de 2019. La generación postcrisis.
eléctrica al cierre del año presentará una
disminución de 5% en relación al año anterior, Para finalizar, una buena noticia: en agosto de
disminuyendo principalmente la generación tér- este año, la revista ENERLAC ha ingresado
mica con GN, CM y derivados de petróleo. en el índice LATINDEX gracias a sus procesos
que aseguran la calidad de sus publicaciones.
Las transiciones energéticas en marcha en la Queremos extender las gracias a todos los co-
región y los compromisos climáticos se verán laboradores: autores, revisores, comité editorial,
ante la incertidumbre de un escenario de con- edición y toda la gente que hace posible mantener
tracción económica y precios relativamente bajos este espacio creado con el objetivo de contribuir
de los hidrocarburos. Las emisiones de gases a la difusión de resultados de investigaciones
de efecto invernadero del sector energético en y trabajos científicos de carácter técnico del
este período han reducido drásticamente (11% sector energético de América Latina y el Caribe,
respecto a 20192), explicado principalmente por elaborados por académicos e investigadores
las restricciones de movilidad y desaceleración de pertenecientes a instituciones vinculadas con el
los sectores industrial y comercial, un resultado sector de los 27 países miembros de OLADE.
positivo para el ambiente.
p. 10
ENSAYO CUASI-DINÁMICO DE COLECTORES SOLARES DE
PLACA PLANA EN URUGUAY DE ACUERDO A LA NORMA
ISO 9806:2017
Juan Manuel Rodríguez Muñoz, Andrés Monetta, Italo Carlos Bove Vanzulli,
Rodrigo Alonso-Suárez
p. 28
ANÁLISIS DE LA EFICIENCIA EN LA PRODUCCIÓN DE ENERGÍA
ELÉCTRICA PARA LA CENTRAL SOLAR FOTOVOLTAICA DE 50
kW CONECTADA A LA RED, EN EL CAMPUS UNIVERSITARIO DE
COTA COTA, UBICADA EN LA CIUDAD DE LA PAZ, BOLIVIA
p. 40
ANÁLISIS EXPERIMENTAL DE LA TRANSMISIÓN SOLAR DE
DOS TIPOS DE PROTECCIONES SOLARES PARA DISTINTAS
POSICIONES EN MONTEVIDEO, URUGUAY
p. 76
ESTUDIO EXPERIMENTAL DE LA PERMEABILIDAD
AL AIRE DE EDIFICIOS EN URUGUAY
p. 92
DESARROLLO DE NUEVOS CATALIZADORES BASADOS
EN NÍQUEL PARA LA PRODUCCIÓN DE GAS NATURAL
SINTÉTICO A PARTIR DE LA HIDROGENACIÓN DE CO2
p. 132
CONTRIBUCIÓN DE LA CADENA PRODUCTIVA DE LA ENERGÍA
EN LAS EMISIONES DE GASES DE EFECTO INVERNADERO EN
ARGENTINA
Maxime Le Bail
p. 146
HACIA UNA JUSTICIA ENERGÉTICA EN MÉXICO. UN ANÁLISIS
DEL RÉGIMEN TARIFARIO DEL SECTOR ELÉCTRICO
COLOMBIANO Y SU POSIBLE IMPLEMENTACIÓN EN LA
REPÚBLICA MEXICANA
Juan Manuel Rodríguez Muñoz 1, Andrés Monetta 2, Italo Carlos Bove Vanzulli 3, Rodrigo Alonso-Suárez 4
10
ENERLAC • Volumen IV. Número 2. Diciembre, 2020. ISSN: 2602-8042 (impreso) / 2631-2522 (digital).
OLADE – AUGM
RESUMEN ABSTRACT
Los ensayos de colectores solares en estado Transient-state efficiency tests for solar thermal
transitorio permiten aumentar la capacidad de collectors allow increasing the laboratories
ensayo de los laboratorios respecto a los ensayos annual test capacity in comparison to steady-
en estado estacionario. Esto ha sido validado en state efficiency tests. This characteristic has
laboratorios internacionales, especialmente en been evaluated by international laboratories,
Europa. En Latinoamérica este método no ha especially in Europe. In Latin America this
tenido mucha difusión y a la fecha ningún labo- method has not been widely disseminated and
ratorio de esta región había validado (ni imple- to date no regional laboratory had validated
mentado) su uso en climas locales, por lo que los (or implemented) its utilization under the local
laboratorios de la región que están desarrollando climates, so laboratories in the region that are
capacidades de ensayo de colectores solares no developing solar collector testing capabilities
consideran aún este método en forma comercial. do not yet consider this method commercially.
En este trabajo se presentan los resultados de This article presents the results of a transient-
un ensayo en estado transitorio (cuasi dinámico) state (quasi-dynamic) test of a flat plate solar
de un colector solar de placa plana de acuerdo collector according to ISO 9806:2017 standard
a la norma ISO 9806:2017 en Uruguay, región in Uruguay, a part of the broader Pampa
de la Pampa Húmeda en Sudamérica. Se compara Húmeda region in South America. Test results are
sus resultados con la metodología de ensayo compared to those of the traditional steady-state
tradicional en estado estacionario y se demues- methodology, demonstrating that the method
tra que es posible aplicarlo en la región incluso can be applied in the region and has also the
con la ventaja de ampliar la cantidad anual de advantage of increasing the annual test capacity,
ensayos, como se ha observado en Europa. as observed in Europe.
Palabras clave: Colectores Solares Térmicos, Keywords: Solar Thermal Collectors, Performance
Ensayos de Desempeño, Modelos Dinámicos, Test, Dynamic Model, Steady State, ISO 9806,
Estado Estacionario, ISO 9806, Uruguay Uruguay.
11
ENSAYO CUASI-DINÁMICO DE COLECTORES SOLARES DE PLACA PLANA EN URUGUAY DE ACUERDO A LA NORMA ISO 9806:2017
Rodríguez Muñoz, Juan Manuel; Monetta, Andrés; Bove Vanzulli, Italo Carlos; Alonso-Suárez, Rodrigo
INTRODUCCIÓN
El aumento sostenido del consumo final de aumente, dado que esto dependerá en definitiva
energía mundial sumado al hecho de que el de las políticas energéticas que adopte cada país.
85% proviene de combustibles fósiles (IEA, No obstante, dicha capacidad es necesaria para
2016a) genera gran preocupación respecto a el desarrollo del sector, ya que difícilmente una
las emisiones de CO2 y su efecto en el cambio política energética tendrá éxito si la calidad de
climático. Las dos acciones más importantes los colectores comercializados no es la adecuada.
para reducir estas emisiones son la incorpora- Además, los ensayos de desempeño térmico
ción de fuentes de energías renovables al uso final constituyen un insumo básico para la evaluación
de energía y la reducción de la demanda a través económico financiera de los proyectos.
de medidas de eficiencia energética. En 2016,
el sector edilicio fue responsable del 30 % del Tradicionalmente los ensayos de eficiencia de
consumo final de energía mundial (IEA, 2016a) colectores solares (desempeño térmico) se
mientras que en Uruguay este sector represento han realizado en estado estacionario (Steady
el 25 % (IEA, 2016b). Una parte importante de State Testing, SST). Para asegurar el estado
este consumo se debe al calentamiento de agua estacionario se requiere mantener con gran
sanitaria (ACS) y a los sistemas de calefacción y estabilidad diversas variables durante el ensayo.
refrigeración. Algunas de ellas son controlables, como la
temperatura del fluido a la entrada del colector
La energía solar térmica de baja temperatura y el caudal, y otras son no controlables, como
se presenta como una candidata para sustituir la irradiancia solar o las condiciones de viento
el uso tradicional de energía para ACS y para sobre el colector. Es posible mitigar la acción
los sistemas de acondicionamiento térmico en del viento en los alrededores forzando viento
el sector edilicio. Entre 2007 y 2017, la capaci- de forma artificial a una velocidad determinada
dad instalada de colectores solares térmicos sobre el colector, elemento que es considerado
aumentó de forma sostenida pasando de 145 GW por todos los laboratorios al exterior. El control
en 2007 a 472 GW en 2017 (REN21, 2018). del ángulo de incidencia es también una limi-
No obstante, si bien la energía solar térmica tante, dado que el movimiento aparente del sol
ha crecido de forma notoria, la participación varía para cada momento del año. Por ende,
de esta tecnología en la matriz energética dada todas las restricciones geométricas, no
mundial sigue siendo pequeña. Para lograr una es sencillo conseguir los ángulos de incidencia
mayor participación se deben abordar distintas necesarios para un ensayo bajo norma, incluso
problemáticas, una de ellas es la calidad de utilizando sistemas de seguimiento solar. Estas
los productos que se comercializan. En este condicionantes, en especial, las relacionadas con
sentido, varias organizaciones internacionales irradiancia solar y ángulo de incidencia, limitan
de reconocido prestigio, como la International la cantidad de ensayos que un laboratorio puede
Organization for Standarization (ISO) y la realizar en un año.
American Society of Heating, Refrigeration and
Air-Conditioning Enginners (ASHRAE), han elabo- Buscando solucionar esta situación, en la década
rado normativas que describen procedimientos de los 90’ se realizaron los primeros ensayos
para evaluar el desempeño térmico de colectores de eficiencia en estado transitorio (Quasi
solares y la calidad de los materiales con los Dynamic Testing, QDT) para colectores solares
que se fabrican. El hecho que un país cuente (Muschaweck y Spirkl, 1993; Perers, 1993) y
con estas capacidades de ensayo no garantiza se comenzaron a evidenciar las ventajas de
que la capacidad instalada de colectores solares este método en sitios con nubosidad variable,
12
ENERLAC • Volumen IV. Número 2. Diciembre, 2020. ISSN: 2602-8042 (impreso) / 2631-2522 (digital).
OLADE – AUGM
donde las condiciones meteorológicas para el con la metodología tradicional. El trabajo busca
ensayo SST son difíciles de conseguir. Algunos además darle mayor difusión a la metodología de
laboratorios en Europa ubicados en sitios con ensayo QDT en Latinoamérica, describiendo una
estas características climáticas adoptaron implementación del método realizada y validada
esta metodología (Fischer et al., 2004; Osório por un laboratorio local.
y Carvalho, 2014; García de Jalón, 2011;
Zambolin, 2010) y también en Estados Unidos Este artículo se organiza de la siguiente forma.
(Rojas, 2008), pero en Latinoamérica, por En primer lugar se describe el modelado
distintos motivos, esta metodología no ha sido matemático de la norma ISO 9806:2017 para
adoptada por la mayoría de los laboratorios, describir el comportamiento térmico de un
existiendo muy pocos antecedentes a la fecha colector solar con cubierta, incluyendo una
(Kratzenberg, 2006). Dado que el método SST descripción de las variables de ensayo rele-
requiere condiciones exigentes de cielo des- vantes. En la siguiente sección (p.17) se pre-
pejado para su correcta implementación, el senta brevemente la plataforma de ensayo, a
método QDT se presenta como una alternativa saber, el Banco de Ensayos de Calentadores
más factible de realizarse en climas con mayor Solares (BECS) del Laboratorio de Energía
presencia de nubosidad, donde los días de Solar (LES) de la Udelar, y se describen las
cielo despejados son escasos. Dadas sus carac- características comunes a ambos ensayos.
terísticas, la metodología QDT no sólo Posteriormente, (p.20) se describe el procedi-
permite, sino que es un requisito de la misma, miento de ensayo utilizado para ambas
que durante el ensayo hayan ocurrido dife- metodologías (SST y QDT) y la siguiente sección
rentes condiciones de cielo (cielo despejado, (p.23) presenta y compara sus resultados. Se
parcialmente nublado y nublado), siendo difícil presenta al final de la misma el análisis de la
de implementar en sitios con predominancias cantidad de ensayos anuales que se podrían
de cielos despejados (en climas áridos o desér- obtener con ambas metodologías, evidenciando
ticos, por ejemplo). Los laboratorios pueden las ventajas del método QDT en la región.
optar por utilizar una u otra metodología según Finalmente, se resumen las conclusiones (p.25).
sean las condiciones meteorológicas locales
(utilizando la misma infraestructura de ensayo).
13
ENSAYO CUASI-DINÁMICO DE COLECTORES SOLARES DE PLACA PLANA EN URUGUAY DE ACUERDO A LA NORMA ISO 9806:2017
Rodríguez Muñoz, Juan Manuel; Monetta, Andrés; Bove Vanzulli, Italo Carlos; Alonso-Suárez, Rodrigo
La norma ISO 9806:2017 (ISO, 2017) propone cubierta. El primer modelo es el más simple y se
dos modelos matemáticos para describir el describe según:
desempeño térmico de un colector solar con
dTm
C = Ag [ ηo,hem . Khem (θ) . G a₁ (Tm Ta ) a₂ (Tm Ta ) ] ḿ́cp ( To Ti ) (1)
dt
Lo que da lugar a:
dTm
C = Ag [ ηo,b . Kb (θ) . Gb + Kd . Gd ] a₁ ( Tm Ta ) a₂ ( Tm Ta ) ] ḿ́cp ( To Ti ) (3)
dt
En el caso del método SST se determinan los destacar que los primeros dos ensayos son en
coeficientes η0,hem, a1, a2, Khem(θ) y C de la estado estacionario pero el tercero no. En el caso
Ecuación (1) en forma separada a través de tres del método QDT se determinan los coeficientes
ensayos independientes: (i) Ensayo de rendi- η0,b, a1, a2, Kb(θ), Kd y C de la Ecuación (3) en
miento en donde se determinan los coeficientes un único ensayo que integra los tres ensayos
η0,hem, a1 y a2; (ii) Ensayo de modificador por independientes del método SST.
ángulo de incidencia en donde se determina
Khem(θ); y (iii) Ensayo de capacidad térmica Para el modificador por ángulo de incidencia
en donde se determina la constante C. Cabe Khem(θ) se utiliza la siguiente relación:
Khem ( θ ) = 1 b0 ( 1
cosθ
1 ) (4)
90
∫0 kb ( θ )sinθcosθdθ
Kd = (5)
90
∫0 kb sinθcosθdθ
14
ENERLAC • Volumen IV. Número 2. Diciembre, 2020. ISSN: 2602-8042 (impreso) / 2631-2522 (digital).
OLADE – AUGM
ηo,hem
ηo,b = (6)
(1 d ) + fd Kd
Las pruebas se realizaron en el Banco de Ensayos Figura 1. Ubicación del LES dentro de la región de
de Calentadores Solares (BECS) del Laboratorio la pampa húmeda (verde)
de Energía Solar (LES, http://les.edu.uy/) de la
Universidad de la República (Udelar) ubicado
en la localidad de Salto (31° sur), Uruguay. En la
Figura 1 se muestra la ubicación del LES dentro
de la región climática afín de la pampa húmeda
(verde claro). Dicha región se ubica al sureste
de Sudamérica y está clasificada como Cfa
(clima templado-cálido, húmedo, con veranos
calurosos) según la clasificación actualizada
de Köppen-Geiger (Peel et al., 2007). Las zonas
identificadas en verde oscuro corresponden
a la clasificación Cfb, muy similar a Cfa, con la
excepción de que los veranos son más frescos
por la influencia de la costa.
15
ENSAYO CUASI-DINÁMICO DE COLECTORES SOLARES DE PLACA PLANA EN URUGUAY DE ACUERDO A LA NORMA ISO 9806:2017
Rodríguez Muñoz, Juan Manuel; Monetta, Andrés; Bove Vanzulli, Italo Carlos; Alonso-Suárez, Rodrigo
Para este trabajo se ensayó un colector solar utilizando los ventiladores axiales que se mues-
térmico de placa plana de 2.01 m2 de área bruta tra en Figura 2. El caudal de ensayo se fijó en
entre el 17 y el 25 de noviembre de 2019. Los 2.4 lts/min y durante todas las pruebas el ángulo
datos se procesaron de diferente forma según de inclinación del seguidor se fijó en 45°, el
el método empleado, SST o QDT. El colector se azimut se ajustó de forma manual o automática
montó sobre un seguidor móvil como se muestra según el ensayo. En todos los casos, antes del
en la Figura 2. La inclinación horizontal de este periodo de medición, al colector se lo hizo pasar
seguidor se puede ajustar manualmente entre por un periodo de acondicionamiento previo
0° y 75° y el azimut se puede ajustar de forma de 15 minutos a la temperatura de ensayo
manual o automática con un paso de 2 minutos correspondiente; dicho período no fue incluido
entre -90° y 90°. Se impuso una velocidad de en la identificación de parámetros de ensayo en
3 m/s (promedio espacial) sobre el colector ninguna de las metodologías.
16
ENERLAC • Volumen IV. Número 2. Diciembre, 2020. ISSN: 2602-8042 (impreso) / 2631-2522 (digital).
OLADE – AUGM
Tamb Visor
T2 Anemómetro
V53
GTI Piranómetro
V Piranómetro de difusa
T4
Sensor de temperatura ambiente
IC1
GV
T1 Sensor de temperatura
Caudalímetro
Circuito primario
Circuito de calentamiento
V51 T3 C1 Circuito de enfriamiento
Cañería flexible
TTE1 Válvula de seguridad
IC2
Intercambiador de calor
T5
17
ENSAYO CUASI-DINÁMICO DE COLECTORES SOLARES DE PLACA PLANA EN URUGUAY DE ACUERDO A LA NORMA ISO 9806:2017
Rodríguez Muñoz, Juan Manuel; Monetta, Andrés; Bove Vanzulli, Italo Carlos; Alonso-Suárez, Rodrigo
Para medir la temperatura a la entrada y a del colector se obtuvo Gb. La medida de irra-
la salida del colector (T1 y T2) se utilizaron diancia difusa en plano inclinado con la banda
PT100 de 3 hilos con transmisores de 4-20 mA de sombra es corregida con un factor provisto
de la marca Herten. Dichos sensores fueron por el fabricante. Dichos piranómetros son
calibrados bajo normativa ISO vigente en el calibrados anualmente en el LES de acuerdo a
LES reportándose una incertidumbre estándar la norma ISO 9747:1992 contra un estándar
de 0.02 °C. La temperatura ambiente (Tamb) secundario Kipp & Zonen CMP22 que es enviado
se registró con un sensor PT1000 de 2 hilos cada cinco años a calibrar contra la referencia
Honeywell también calibrado en el LES con una radiométrica mundial en World Radiation
incertidumbre estándar de 0.02 °C. La medida Center en Davos, Suiza. Todas las medidas
de caudal (C1) se realizó con un caudalímetro se registraron cada 10 segundos usando un
electromagnético de Endress & Hauser con una adquisidor de datos Fischer Scientific DT85.
incertidumbre estándar de 0.5 % de la medida.
La velocidad del viento paralelo al plano del PROCEDIMIENTO
colector (V) se midió con un anemómetro
de copa NGR con incertidumbre estándar de Condiciones requeridas
0.25 m/s. La irradiancia global en el plano del
colector (GTI) se midió con un piranómetro La Tabla 1 muestra las condiciones requeridas
Kipp & Zonen CMP 10. La irradiancia directa por cada metodología para cada variable,
Gb en el plano del colector se estimó de la incluyendo la variabilidad permitida. Estas son
siguiente manera: por un lado, se midió la las condiciones que deben cumplir las medidas
irradiancia global en plano horizontal (GHI) registradas durante el ensayo para poder ser
con un piranómetro Kipp & Zonen CMP11 y la utilizadas para determinar los parámetros con
irradiancia difusa en plano horizontal (DHI) cada metodología. Como se puede observar, la
con un piranómetro Kipp & Zonen CMP6 mon- metodología SST presenta requerimientos más
tado con una banda de sombra del mismo estrictos respecto a los valores requeridos y
fabricante, luego, la irradiancia en incidencia su variabilidad. Por otro lado, la metodología
normal (DNI) se estimó por diferencia de estas QDT requiere que se representen diferentes
dos medidas, y, multiplicando la DNI por el condiciones meteorológicas durante el ensayo,
coseno del ángulo de incidencia en el plano como se explica en la página 24.
18
ENERLAC • Volumen IV. Número 2. Diciembre, 2020. ISSN: 2602-8042 (impreso) / 2631-2522 (digital).
OLADE – AUGM
Ensayo de Rendimiento: Para este ensayo se estacionario, para cada punto, se tomaron
midió la energía útil por unidad de área para promedios en períodos sucesivos de 30 segundos
cuatro temperaturas de entrada diferentes, y y se verificaron que las variaciones de dichos
la temperatura de ensayo más alta fue 85 ⁰C. promedios respecto al promedio 10 minutal se
A su vez, por cada temperatura de entrada, encontraran dentro de los límites dados por la
se obtuvieron cuatro puntos independientes Tabla 1. Finalmente, los parámetros η0,hem, a1
dando así un total de 16 puntos. Cada punto y a2 se determinaron por regresión lineal
consistió en un promedio sobre un período de usando el siguiente modelo:
medición 10 minutos. Para verificar el estado
Qu = η0,hem . G - a₁ (Tm Ta ) a₂ (Tm Ta ) (7)
ḿcp (To Ti )
Qu = (8)
Ag
Qu + a₁ (Tm Ta ) + a₂ (Tm Ta )²
Khem (Ѳ) = (9)
ηo,hem . G
19
ENSAYO CUASI-DINÁMICO DE COLECTORES SOLARES DE PLACA PLANA EN URUGUAY DE ACUERDO A LA NORMA ISO 9806:2017
Rodríguez Muñoz, Juan Manuel; Monetta, Andrés; Bove Vanzulli, Italo Carlos; Alonso-Suárez, Rodrigo
Ensayo de Capacidad Térmica: Este ensayo se se removió y se dejó que el colector alcance un
hizo de acuerdo a la sección 25.2 de la norma nuevo punto en estado estacionario, diferente
ISO 9806:2017 pero teniendo en consideración al primero por el efecto de la irradiancia solar.
el ajuste de segundo orden en las pérdidas La capacidad térmica efectiva se determinó
térmicas, es decir el coeficiente a2. Al inicio del integrado la Ecuación (1) sobre el periodo de
ensayo se fijó una temperatura de entrada igual tiempo entre los dos puntos de funcionamiento
a la temperatura ambiente y luego se cubrió el estacionarios, como se muestra en la Ecuación
colector con una manta reflectiva dejando que (10). Para esto se utilizaron los datos
alcance el estado estacionario. Luego, la cubierta promediados cada 30 segundos.
C
=
ʃ t1
t2
(η 0, hem . G a₁ (Tm Ta ) a₂ (Tm Ta )² Q dt
u ) (10)
Ag Tm2 Tm1
En este caso, también se realizó un seguimiento tomó una serie de datos realizando seguimiento
acimutal y por este motivo en la Ecuación (10) se acimutal, de modo de obtener medidas con
asume incidencia normal (Khem = 1). bajo ángulo de incidencia, y otra serie de datos
fijando el azimut del seguidor en 0° (orientación
Ensayo en estado Cuasi-Dinámico QDT norte) para obtener ángulos de incidencia altos.
Ambas series se tomaron bajo condiciones
Para este ensayo se midió la energía útil por de cielo claro. Las temperaturas restantes se
unidad de área para cuatro temperaturas de en- realizaron bajo distintas condiciones de cielo
trada diferentes. Los datos fueron promediados realizando seguimiento acimutal para trabajar
cada 5 minutos y para cada temperatura el con ángulos de incidencia pequeños (Khem = 1).
periodo de medición fue de al menos 3 horas. Si bien esto no es requisito de la norma,
Se verificó además que durante cada período permite un mayor desacople de las variables
de medición se cumplieran las condiciones de dependientes y mejora la identificación de
la Tabla 1, comparando los promedios cada 30 parámetros. Finalmente, las constantes η0,b, a1,
segundos con los promedios de todo el período a2, b0, Kd y C se determinaron todas a la vez por
de medición. Para la temperatura más baja se regresión lineal utilizando el modelo siguiente:
(Ec. Q
11) u
= η0,b . Gb η0,b b0
1
cosθ ( )
- 1 Gb + η0,b . Kd . Gd a₁ (Tm Ta ) a₂ (Tm Ta )² C dTm
Ag dt (11)
20
ENERLAC • Volumen IV. Número 2. Diciembre, 2020. ISSN: 2602-8042 (impreso) / 2631-2522 (digital).
OLADE – AUGM
RESULTADOS
En la Tabla 2 se presentan los resultados de los cada parámetro (cociente entre el valor de un
ensayos para ambos métodos. En dicha tabla parámetro y su incertidumbre). No se reportó la
se presentan los coeficientes de los modelos incertidumbre de los parámetros η0,b y Kd para el
térmicos de las Ecuaciones (1) y (3), sus res- caso SST por tratarse de magnitudes derivadas
pectivas incertidumbres y el estadístico t para (Ecuaciones (5) y (6)).
Tabla 2. Resultados
En todos los casos se obtuvo un estadístico y Carvalho, 2014). Una forma más adecuada
t mayor 3, que indica que los parámetros son de comparar los coeficientes de pérdidas es
estadísticamente significativos. Los parámetros utilizando un coeficiente de pérdida linealizado
η0,b y b0 difieren en menos de 1% (cociente a una diferencia de temperatura dada ∆T, eso es:
entre diferencia de los valores y el valor
promedio) lo cual da cuenta de una buena a (ΔT) = a₁ + a₂ ΔT (13)
concordancia entre estos parámetros. En el caso
de los coeficientes a1 y a2 la diferencia es mayor, Para una diferencia de 50 K se tienen
3% y 18% respectivamente. Ahora, si se presta a(50 K) = 4.585 W/m2.K para el método SST y
atención a la Tabla 2 se puede observar que el a(50 K) = 4.553 W/m2.K para el caso QDT, lo
término a1 en el caso SST es mayor que en el QDT que da una diferencia menor al 1%. Las mayores
y el término a2 en el caso SST es menor que en diferencias se presentan en el modificador de
el caso QDT, es decir, que la diferencia tiende a ángulo para la irradiación difusa Kd y para la
compensarse. El signo negativo de la covarianza capacidad térmica por unidad de área C/Ag, 5%
entre estos dos parámetros, obtenido en los y 12% respectivamente. Entendemos que la
dos casos, confirma esta relación inversa. Este diferencia en Kd puede deberse en parte al factor
comportamiento no es nuevo; lo mismo se obser- de corrección empleado para corregir la media
va en las referencias (Fischer et al., 2004; Osório de difusa en plano horizontal por el efecto de
21
ENSAYO CUASI-DINÁMICO DE COLECTORES SOLARES DE PLACA PLANA EN URUGUAY DE ACUERDO A LA NORMA ISO 9806:2017
Rodríguez Muñoz, Juan Manuel; Monetta, Andrés; Bove Vanzulli, Italo Carlos; Alonso-Suárez, Rodrigo
la banda de sombra. Este factor supone que la demasiado grande y por este motivo tal vez no
difusa es isotrópica, por lo tanto, este factor se se capturó correctamente el comportamiento
puede mejorar utilizando factores de corrección dinámico del colector. Todos estos aspectos son
que contemplen la anisotropía de la atmósfera motivo de nuestro trabajo actual, en particular,
(Batlle et al., 1995). Otro causante de esta en lo que refiere a los factores de corrección de
diferencia puede ser que las horas bajo la radiación difusa.
condiciones de cielo parcialmente y totalmente
cubierto durante el ensayo realizado estén Finalmente, para estudiar el efecto combinado de
sub-representadas, dado que una adecuada las diferencias entre los parámetros, se calculó
representación de estas condiciones es im- para cada caso la energía útil producida por el
portante para la determinación de este colector utilizando la Ecuación (11), asumiendo
parámetro bajo la metodología QDT. La diferencia incidencia normal y estado estacionario, para
en la capacidad térmica efectiva puede deberse distintas condiciones de temperatura y de
a dos motivos: (i) puede ser que la cantidad de cielo. Para las distintas condiciones de cielo se
horas en condiciones de cielo variable (variable utilizaron las condiciones de reporte estándar
dTm/dt alta en valor absoluto) estén sub- (Standard Reporting Conditions, SRC) dadas
representadas, y (ii) que el intervalo de tiempo por la norma ISO 9806:2017. Los resultados se
5 minutos para promediar los datos sea muestran en la Tabla 3.
Tabla 3. Potencia útil (W) producida por el colector bajo distintas condiciones de cielo y temperatura
Para condiciones de cielo claro la diferencia en el modificador por ángulo de incidencia para
en la energía útil es poco significativa (menor la irradiancia difusa Kd, lo que es motivo actual
al 1%). En condiciones de cielo parcialmente de estudio.
nublado la diferencia comienza crecer a medida
que crece la diferencia de temperatura, si bien Estimación de cantidad de ensayos anuales
sigue siendo poco significativa. En condiciones
de cielo cubierto la diferencia comienza a ser Para mostrar la aplicabilidad del método QDT en
significativa alcanzando el 15 % a una diferencia el clima de la región, en la Tabla 4 se muestra
de temperatura de 40 K. Entendemos que esa un estimativo de la cantidad de ensayos posibles
diferencia se debe principalmente a la diferencia por mes en la localidad de Salto (Uruguay)
22
ENERLAC • Volumen IV. Número 2. Diciembre, 2020. ISSN: 2602-8042 (impreso) / 2631-2522 (digital).
OLADE – AUGM
tanto para seguir el método SST realizando QDT, tiempo estimado en el cual se logra medir
un seguimiento acimutal como para seguir el las variables necesarias para cada ensayo si se
método QDT con una orientación norte, en cumplen las condiciones de la Tabla 1. Además,
ambos casos, con una inclinación horizontal fija para el caso QDT se impuso que al menos un
de 45°. Una hora se considera válida para un tipo 30% de las horas fueran con cielo claro,
de ensayo si cumple las condiciones presentadas debido a que es importante para este ensayo
en la Tabla 1. Para el caso SST se verificó capturar cierta dispersión en las condiciones
además que las horas fueran en condiciones meteorológicas, y las condiciones de cielo
de cielo claro. El análisis se realizó utilizando despejado no deben ser sub-representadas.
el año meteorológico típico de la localidad de Para verificar esta condición se impuso que el
Salto (latitud = 31.28° S, longitud = 57.92° O, índice de cielo claro fuera mayor a 0.8 y que la
altitud = 56 m sobre el nivel del mar), Uruguay, irradiación horaria en incidencia normal (DNI)
a escala horaria (Alonso-Suárez et al., 2019). fuera mayor a 200 Wh/m2. Para el cálculo del
Se supuso una duración de los ensayos de 10 índice de cielo claro se utilizó la expresión dada
horas para el caso SST y de 15 horas para el caso por Pérez et al. (1990):
kt
k´t = (14)
1 + 1.031e ( )
-1.4m
9.4 + 0.9m
MES ENE FEB MAR ABR MAY JUN JUL AGO SEP OCT NOV DIC Anual
SST 1.0 1.2 3.5 5.0 2.3 2.6 3.4 3.4 4.6 2.7 1.9 1.5 33
QDT 2.0 2.4 5.8 10.2 6.8 8.8 10.2 9.2 9.2 5.4 3.8 3.0 76
En la Tabla 4 se puede ver que la cantidad de llegar a realizarse si se siguiera este método.
ensayos posibles en un año en el caso QDT es más Esta cantidad es mayor al doble de la cantidad
del doble que la cantidad de ensayos en el caso de ensayos que podrían obtenerse con la
SST. Esta proporción coincide con lo estimado metodología tradicional en estado estacionario,
en otras localidades (Rojas, 2008; Osório y lo cual muestra la ventaja de la metodología
Carvalho, 2014). QDT en Uruguay y en la región climática
homogénea de la Pampa Húmeda del sureste
CONCLUSIONES de Sudamérica.
Se estudió la aplicabilidad del método cuasi Por otra parte, se concluye que el método QDT
dinámico en Uruguay y se determinó la se implementó correctamente para el colector
cantidad de ensayos a nivel anual que podrían estudiado y se puede utilizar en el clima de
23
ENSAYO CUASI-DINÁMICO DE COLECTORES SOLARES DE PLACA PLANA EN URUGUAY DE ACUERDO A LA NORMA ISO 9806:2017
Rodríguez Muñoz, Juan Manuel; Monetta, Andrés; Bove Vanzulli, Italo Carlos; Alonso-Suárez, Rodrigo
•••
AGRADECIMIENTOS
Los autores desean agradecer al Ministerio de para el desarrollo del BECS y haber impulsado este
Industria, Energía y Minería (MIEM, Uruguay), en proyecto con capacidades nacionales. Los autores
especial a su Dirección Nacional de Energía (DNE), agradecen también al PTB de Alemania por im-
al Fideicomiso Uruguayo de Ahorro y Eficiencia pulsar y financiar el inter-laboratorio de ensayos
Energética (Fudaee, Uruguay) y a la Corporación de eficiencia de colectores solares, lo que nos ha
Nacional para el Desarrollo (CND, Uruguay), por brindado certezas técnicas sobre la capacidad de
haber proporcionado apoyo financiero y logístico ensayo local.
24
ENERLAC • Volumen IV. Número 2. Diciembre, 2020. ISSN: 2602-8042 (impreso) / 2631-2522 (digital).
OLADE – AUGM
NOMENCLATURA
25
ENSAYO CUASI-DINÁMICO DE COLECTORES SOLARES DE PLACA PLANA EN URUGUAY DE ACUERDO A LA NORMA ISO 9806:2017
Rodríguez Muñoz, Juan Manuel; Monetta, Andrés; Bove Vanzulli, Italo Carlos; Alonso-Suárez, Rodrigo
REFERENCIAS
Alonso-Suárez R., Bidegain M., Abal G., Modernell P. (2019). Kovács, P., Petterson, U., Persson, M., Perers, B., Fischer, S.
Año Meteorológico Típico para aplicaciones de Energía Solar (2011). Improving the accurancy in performance prediction
(AMTUes): series horarias típicas para 5 sitios del Uruguay. for new collector deings. Proceedings of ISES Solar World
Ediciones Universitarias (UCUR) de la Universidad de la Congress, 28 august - 2 september, Kassel, Germany.
República (Udelar), ISBN 978-9974-0-1674-7.
Kratzenberg, M.G., Beyer, H.G., Colle, S. (2006). Uncertainty
Batlles, F.J., Olmo, F.J. y Alados-Armboledas, L. (1995). On calculation applied to different regression methods in the
shadowband correction methods for diffuse irradiance quasi-dynamic collector test. Solar Energy, 80, 1453-1462.
measurements. Solar Energy, 54(2): 105-114.
Muschaweck, J. y Spirkl, W. (1993). Dynamic solar collector
Duffie, J.A. y Beckman, W.A. (2013). Solar Engineering of performance testing. Solar Energy Materials and Solar Cells,
Thermal Processes. Wiley and Sons, Hoboken, New Jersey, 30(2):95-105.
third edition, 2006.
Osório, T. y Carvalho, M. (2014). Testing of solar thermal
Fischer, S., Heidemann, W., Müller-Steinhagen, H., Perers, collectors under transient conditions. Solar Energy, 104,
B., Bergquist, P., Hellström, B. (2004). Collector test method 71-81.
under quasi-dynamic conditions according to the European
Standard EN 12975-2. Solar Energy, 76(1-3):117-123. Peel, M.C., Finlayson, B.L. y Mcmahon, T.A. (2007). Update
world mapo of the köppen-geiger climate classification.
Fisher, S. (2020). Expert Report - Quality Infrastructure for Hydrology and Earth System Sciences Discussions, 11, 1633-
Energy Efficiency and Renewable Energy in Latin America 1644.
and the Caribbean. + Report #95309. Recuperado de
http://les.edu.uy/papers/Expert_Report-RR-Test_R1.pdf Perers, B. (1993). Dynamic method for solar collector
array testing and evaluation with standard database and
García de Jalón, A., Sallaberry, F., Olano, X., Mateu, E., Astiz simulation programas. Solar Energy, 50(6), 517-526.
R., Ezcurra, M., Ramírez, L. (2011). Comparison of thermal
efficiency curves of solar collectors tested in outdoor Pérez, R., Ineichen, P., Seals, R. y Zelenka, A. (1990). Making
conditions. Proceedings of ISES Solar World Congress, full use of the clearness index for parameterizing hourly
28 august - 2 september, Kassel, Germany. insolation conditions. Solar Energy 45(2), 111-114.
International Energy Agency (IEA). (2016a). Final Energy Renewable Energy Policy Newtwork for the 21st Century
Consumption in the World. Recuperado de https://www.iea. (REN21). (2018). Renewables 2018 Global Status Report.
org/Sankey/#?c=World&s=Final%20consumption Recuperado de http://www.ren21.net/gsr-2018/pages/
foreword/foreword/
International Energy Agency (IEA). (2016b). Final Energy
Consumption in Uruguay. Recuperado de https://www.iea. Rojas, D., Beermann, S.A., Klein, S.A. y Reindl, D.T. (2008).
org/Sankey/#?c=Uruguay&s=Final20consumption Thermal performance testing of flat-plate collectors. Solar
Energy, 82, 746-757.
ISO 9806. (2017). Solar Energy – Solar thermal collectors
– Test methods. Switzerland: International Organization of Zambolin, E. y Del Col D. (2010). Experimental analysis of
Standarization. thermal performance of flat plate and evacuated tube solar
collectors in stationary standard and daily conditions. Solar
Energy, 84, 1382-1396.
26
ANÁLISIS DE LA EFICIENCIA EN LA PRODUCCIÓN
DE ENERGÍA ELÉCTRICA PARA LA CENTRAL SOLAR
FOTOVOLTAICA DE 50 kW CONECTADA A LA RED,
EN EL CAMPUS UNIVERSITARIO DE COTA COTA,
UBICADA EN LA CIUDAD DE LA PAZ, BOLIVIA
EFFICIENCY ANALYSIS OF THE PRODUCTION OF ELECTRICITY IN THE PHOTOVOLTAIC
SOLAR CENTRAL OF 50 kW, IN THE UNIVERSITY CAMPUS OF COTA COTA,
LOCATED IN LA PAZ, BOLIVIA
Álvaro Montano Saavedra1, María Belén Aruquipa Saravia2, Cristian Limachi Ochoa3
28
ENERLAC • Volumen IV. Número 2. Diciembre, 2020. ISSN: 2602-8042 (impreso) / 2631-2522 (digital).
OLADE – AUGM
29
ANÁLISIS DE LA EFICIENCIA EN LA PRODUCCIÓN DE ENERGÍA ELÉCTRICA PARA LA CENTRAL SOLAR FOTOVOLTAICA DE 50 kW CONECTADA A LA RED, EN EL CAMPUS
UNIVERSITARIO DE COTA COTA, UBICADA EN LA CIUDAD DE LA PAZ, BOLIVIA. Montano Saavedra, Álvaro; Aruquipa Saravia, María Belén; Limachi Ochoa, Cristian
ABSTRACT INTRODUCCIÓN
The article has as objective to perform the analysis Bolivia es un país que cuenta con grandes
of electrical energy on the 50 kW Photovoltaic posibilidades de implementación de sistemas
Solar Central located in Cota Cota’s University de generación con energías alternativas y
Campus, which was built thanks to a donation renovables, y que ha tenido gran relevancia en
from Japan International Cooperation Agency los últimos años con la construcción de pequeñas
(JICA), that is now under the charge of Universidad centrales fotovoltaicas, parques eólicos y pe-
Mayor de San Andrés’ Electric Engineering queñas centrales hidroeléctricas conectadas al
Investigation Institute. This Central began to Sistema Interconectado Nacional. Sin embargo,
operate in early 2018, and till the date has stored el concepto de generación distribuida todavía
periodic information regarding radiation and no fue colocado en práctica, debido a varios
temperature on photovoltaic panels, generated factores, como ser: normativos, regulatorios,
energy of the Photovoltaic Solar Central and técnicos, culturales y otros.
injected in the electric distributor company’s
grid circuit of La Paz (DELAPAZ). The Central En ese sentido, la energía solar fotovoltaica se
was designed under the concept of distributed presenta como una opción para iniciar el de-
generation, which means that all the generated sarrollo del concepto de generación distribuida;
energy is injected to the grid. It is considered that entonces, el realizar investigaciones específicas
due to its geographical location, the photovoltaic acerca del desempeño de Sistemas Solares
panels’ efficiency should be greater than stipulated Fotovoltaicos conectados a la Red de Energía
on standard conditions; however, the results of Eléctrica de la empresa distribuidora en la
the analysis made to collected data, shows that ciudad de La Paz (que por sus características
photovoltaic panels’ efficiency is around 17%, topográficas y climáticas cuenta con aspectos
and it must be considered other factors affecting especiales que deben ser tomados en cuenta)
efficiency: cleanliness and maintenance of permitirá conocer información relevante acerca
photovoltaic panels, shadow effect, etc. de: la radiación solar en la zona, la eficiencia de los
paneles fotovoltaicos en la altura, el rendimiento
del conjunto de todo el sistema fotovoltaico,
Keywords: Photovoltaic Solar Central, Energy además de otros aspectos importantes; así
Generation, Energy Efficiency, Photovoltaic Panels, mismo, posibilitará estimar la cantidad de
Irradiation, Bolivia. energía eléctrica que se podrá generar en
ciertos períodos de tiempo y realizar la difusión
y transferencia de conocimientos de éste tipo
de sistemas, entre los alumnos y docentes de la
Bolivia es un país con Facultad de Ingeniería de la Universidad Mayor
grandes posibilidades de de San Andrés.
implementación de sistemas
de generación con energías Al respecto, se debe mencionar que actualmente
alternativas y renovables, se cuenta con escasa bibliografía relacionada
sin embargo, el concepto con estudios realizados al desempeño de sis-
temas fotovoltaicos en localidades situadas a
de generación distribuida
alturas mayores a 3,000 m.s.n.m.; por otra parte,
todavía no fue colocado
en Bolivia se cuenta con un mapa solar, el cual
en práctica, debido a promedia extensas regiones de 1° de latitud por
varios factores. 1° de longitud en base a imágenes satelitales
de la NASA y un método de interpolación, para
30
ENERLAC • Volumen IV. Número 2. Diciembre, 2020. ISSN: 2602-8042 (impreso) / 2631-2522 (digital).
OLADE – AUGM
31
ANÁLISIS DE LA EFICIENCIA EN LA PRODUCCIÓN DE ENERGÍA ELÉCTRICA PARA LA CENTRAL SOLAR FOTOVOLTAICA DE 50 kW CONECTADA A LA RED, EN EL CAMPUS
UNIVERSITARIO DE COTA COTA, UBICADA EN LA CIUDAD DE LA PAZ, BOLIVIA. Montano Saavedra, Álvaro; Aruquipa Saravia, María Belén; Limachi Ochoa, Cristian
RESULTADOS
Para el análisis del desempeño de la CSFV y la acerca de las variables: irradiación, temperatura
determinación de la eficiencia de los PFVs, se y energía eléctrica generada. En la Figura 1 se
recopiló información desde el momento de puede observar la relación existente entre la
entrada en operación de la Central (marzo energía solar y la energía eléctrica producida
de 2018) hasta el mes de mayo del 2019; en la CFV, durante una semana del mes de
debido a la gran cantidad de información octubre de la gestión 2018. La Figura 2 presenta
disponible, en primera instancia se identifi- la información recopilada que consta de más
caron los períodos de tiempo en los cuales de 400 datos, correspondientes al día 12 de
se disponía de la totalidad de información octubre.
Figura 1. Energía Solar recibida y Energía Eléctrica producida en la CFV evaluadas en la semana del
8 al 12 de octubre de la gestión 2018
700
Energía Solar Energía Eléctrica
600
Energía [kW-h/día]
500
400
300
200
100
0
8/10/2018 9/10/2018 10/10/2018 11/10/2018 12/10/2018 13/10/2018 14/10/2018
Días
Fuente: Elaboración propia con información del Sistema de adquisición de datos de la Central Fotovoltaica (Campus
Universitario de Cota Cota, 2018).
32
ENERLAC • Volumen IV. Número 2. Diciembre, 2020. ISSN: 2602-8042 (impreso) / 2631-2522 (digital).
OLADE – AUGM
1600,0 20,0
18,0
1400,0
16,0
1200,0
TEMPERATURA [˚C]
14,0
IRRADIACIÓN [W/m2]
1000,0
12,0
800,0 10,0
8,0
600,0
6,0
400,0
4,0
200,0
2,0
0,0 0,0
06:02:50
06:16:20
06:26:20
06:36:30
06:46:30
06:56:40
07:06:40
07:16:50
07:26:50
07:36:00
07:47:00
07:57:10
08:07:10
08:18:20
08:27:40
08:39:40
08:47:30
08:57:40
09:10:10
09:17:40
09:27:50
09:37:50
09:50:10
10:00:20
10:08:10
10:20:10
10:28:20
10:40:40
10:48:30
10:58:30
11:10:00
11:18:50
11:28:50
11:38:50
11:51:10
12:00:00
12:12:20
12:22:00
12:32:30
12:42:30
12:52:30
13:01:20
13:12:50
13:22:50
13:33:00
13:43:00
13:53:30
14:02:40
14:13:10
14:23:20
14:33:20
14:42:00
14:53:30
15:02:50
15:13:40
15:23:50
15:33:50
15:43:20
15:53:50
16:02:00
16:14:10
16:24:20
16:34:10
16:44:30
16:53:50
17:03:50
17:14:30
17:24:40
17:34:50
17:45:00
17:55:00
18:03:20
18:12:20
18:24:10
HORAS
IRRADIACIÓN W/m2 TEMPERATURA °C
Fuente: Elaboración propia con datos del estudio de irradiación y la temperatura ambiente (Campus Universitario de
Cota Cota, 2018).
33
ANÁLISIS DE LA EFICIENCIA EN LA PRODUCCIÓN DE ENERGÍA ELÉCTRICA PARA LA CENTRAL SOLAR FOTOVOLTAICA DE 50 kW CONECTADA A LA RED, EN EL CAMPUS
UNIVERSITARIO DE COTA COTA, UBICADA EN LA CIUDAD DE LA PAZ, BOLIVIA. Montano Saavedra, Álvaro; Aruquipa Saravia, María Belén; Limachi Ochoa, Cristian
170 Energía
160 Solar
150
Energía
140 Eléctrica
130
120
110
100
90
kW
80
70
60
50
40
30
20
10
0
00:00:00
01:00:00
02:00:00
03:00:00
04:00:00
05:00:00
06:00:00
07:00:00
08:00:00
09:00:00
10:00:00
11:00:00
12:00:00
13:00:00
14:00:00
15:00:00
16:00:00
17:00:00
18:00:00
19:00:00
20:00:00
21:00:00
22:00:00
23:00:00
00:00:00
HORAS
Fuente: Elaboración propia con datos del estudio de energía solar recibida y energía eléctrica generada (Campus Universitario
de Cota Cota, 2018).
34
ENERLAC • Volumen IV. Número 2. Diciembre, 2020. ISSN: 2602-8042 (impreso) / 2631-2522 (digital).
OLADE – AUGM
Figura 4. Irradiación solar en el Campus de Cota Cota. 16 ° 32´ 17´´S, 68 ° 03´ 44´´ O 3,492 m.s.n.m.
1600
1400
IRRADIACIÓN - W/m2
1200
1000
800
600
400
200
0
ENE FEB MAR ABR MAY JUN JUL AGO SEP OCT NOV DIC
Fuente: Elaboración propia con datos del estudio de irradiación (Campus Universitario de Cota Cota, 2018).
35
IRRADIACIÓN [W/m2] IRRADIACIÓN [W/m2]
1000
1200
1400
1600
200
400
600
800
0
200
400
600
800
1000
1200
1400
00:02:00 00:12:10
00:18:40 00:29:00
00:35:30 00:45:40
00:52:20 01:02:30
01:09:10 01:19:20
01:26:00 01:36:10
IRRADIACIÓN W/m2
IRRADIACIÓN W/m2
07:35:30 08:00:10
07:52:20 08:02:30
08:09:10 08:24:40
08:37:10 08:38:50
08:50:10 09:01:00
09:10:30 09:11:30
09:16:30 09:38:10
09:41:30 09:43:20
09:52:20 10:00:10
10:17:00 10:24:40
10:26:00 10:33:40
10:49:00 10:02:40
10:57:10 11:10:40
11:14:30 11:24:50
11:31:40 11:50:00
11:54:40
36
12:04:30 12:07:00
12:18:40
Figura 5. Irradiación y Temperatura en la CFV, Verano
12:27:50 12:31:20
HORAS
HORAS
12:41:50
TEMPERATURA °C
TEMPERATURA °C
22:26:00 22:19:20
22:42:40 22:36:10
22:59:30 22:53:00
23:16:20 23:09:40
23:33:10 23:26:30
23:50:00 23:43:20
0,0
2,0
4,0
6,0
8,0
0,0
5,0
UNIVERSITARIO DE COTA COTA, UBICADA EN LA CIUDAD DE LA PAZ, BOLIVIA. Montano Saavedra, Álvaro; Aruquipa Saravia, María Belén; Limachi Ochoa, Cristian
10,0
15,0
20,0
25,0
10,0
12,0
14,0
16,0
18,0
20,0
Fuente: Elaboración propia con datos del estudio de irradiación y temperatura en invierno (Campus Universitario de
Fuente: Elaboración propia con datos del estudio de irradiación y temperatura en verano (Campus Universitario de
ANÁLISIS DE LA EFICIENCIA EN LA PRODUCCIÓN DE ENERGÍA ELÉCTRICA PARA LA CENTRAL SOLAR FOTOVOLTAICA DE 50 kW CONECTADA A LA RED, EN EL CAMPUS
Figura 7. Irradiación solar en la CFV del campus de Cota Cota, solsticio de verano y solsticio de invierno
1600
1400
IRRADIACIÓN [W/m2]
1200
1000
800
600
400
200
0
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24
HORAS
VERANO INVIERNO
Fuente: Elaboración propia con datos del estudio de irradiación en Verano e invierno (Campus Universitario de Cota
Cota, 2018)
DISCUSIÓN
Los resultados encontrados luego de la También se deben tomar en cuenta que existen
recopilación de datos y análisis realizados, otros factores externos que influencian en
muestran que el desempeño de la Central Solar el desempeño de la CSFV, como el estado de
Fotovoltaica es aceptable, llegando a obtener limpieza de los paneles y el efecto sombra que
una eficiencia de los PFVs cercana al 17% para el puede ser producido por los árboles presentes
periodo analizado. Se debe tomar en cuenta que en la zona.
la CSFV opera en condiciones muy diferentes a
las estipuladas en las denominadas condiciones CONCLUSIONES
estándar y que, para poder determinar criterios
de evaluación más certeros, de deben analizar El valor de la eficiencia encontrado en relación
la mayor cantidad de datos posible, en mayores a la energía solar recibida por los paneles
períodos de tiempo (semestrales, anuales). fotovoltaicos y la energía eléctrica generada por
37
ANÁLISIS DE LA EFICIENCIA EN LA PRODUCCIÓN DE ENERGÍA ELÉCTRICA PARA LA CENTRAL SOLAR FOTOVOLTAICA DE 50 kW CONECTADA A LA RED, EN EL CAMPUS
UNIVERSITARIO DE COTA COTA, UBICADA EN LA CIUDAD DE LA PAZ, BOLIVIA. Montano Saavedra, Álvaro; Aruquipa Saravia, María Belén; Limachi Ochoa, Cristian
•••
AGRADECIMIENTOS
38
ConvocatorIa
para POSTULACIÓN de artículos
Más InformacIÓN
CALL FOR
SUBMISSION OF ARTICLES
MORE INFORMATION
[email protected]
ANÁLISIS EXPERIMENTAL DE LA TRANSMISIÓN
SOLAR DE DOS TIPOS DE PROTECCIONES SOLARES
PARA DISTINTAS POSICIONES EN MONTEVIDEO,
URUGUAY
EXPERIMENTAL ANALYSIS OF THE SOLAR TRANSMISSION OF TWO TYPES OF SOLAR
PROTECTIONS FOR DIFFERENT POSITIONS IN MONTEVIDEO, URUGUAY
40
ENERLAC • Volumen IV. Número 2. Diciembre, 2020. ISSN: 2602-8042 (impreso) / 2631-2522 (digital).
OLADE – AUGM
41
ANÁLISIS EXPERIMENTAL DE LA TRANSMISIÓN SOLAR DE DOS TIPOS DE PROTECCIONES SOLARES PARA DISTINTAS POSICIONES EN MONTEVIDEO, URUGUAY.
Pereira Ruchansky, Lucía; Gutiérrez Bazterrica, Lucía; Sosa Ibarra, Daniel; López Salgado, María Noel; Camacho Roberts, Magdalena
position for three orientations (east, north and la principal fuente de pérdida de calor en
northwest) and in a position perpendicular to the invierno, por su alta transmitancia térmica, y la
incident direct solar radiation . principal fuente de ganancia de calor en verano
por su comportamiento frente a la radiación solar
The results of solar transmission obtained (transparente a la radiación solar de onda corta
experimentally, allow us to identify that the y opaca a la radiación solar de onda larga
behavior of the solar protections is closely linked (Serra y Coch, 1995; Maraveas, 2019).
to their position relative to the sun - angle of
incident radiation - and to the geometry of the En particular, para optimizar el desempeño
solar protection; putting into question the solar térmico y prevenir el sobrecalentamiento en
factor as a representative parameter of the el interior de los locales durante el periodo
behavior of solar protections under conditions of caluroso es fundamental la correcta utilización de
use and its usefulness as a design tool. estrategias de diseño que incorporen el control de
ganancias de energía por radiación solar, así como
considerar desde el proyecto las propiedades
Keywords: Solar Transmission, Solar Factor, térmicas y lumínicas de los componentes de estos
Solar Protection, Solar Radiation, Energy cerramientos.
Efficiency, Fenestration System, Uruguay.
Existen dos indicadores que se relacionan con
el manejo de la energía térmica y se utilizan pa-
ra evaluar el desempeño de los cerramientos
INTRODUCCIÓN transparentes: transmitancia térmica U (W/m2.K)
y factor solar (FS). El FS, o coeficiente de
ganancia solar (por sus siglas en inglés SHGC),
La tendencia de la arquitectura en las últimas es la fracción de irradiancia incidente que
décadas promueve la utilización del vidrio y la ingresa por una ventana, e incluye la energía
transparencia como lenguaje expresivo frente transmitida directamente y la reemitida, ver
a una fascinación por la luz y la ligereza en el Figura 1. Es un valor entre 0 y 1, donde 0 indica
diseño. En muchos casos, estas expresiones que no se transmite radiación y 1 significa que se
estéticas se anteponen a la posibilidad de transmite la totalidad de la radiación incidente;
brindar condiciones de confort térmico y visual puede especificarse para cada componente del
a los ocupantes de forma natural, condicionando cerramiento por separado (vidrio y dispositivo
al uso de ventilación mecánica y luz artificial. de protección solar) o para todo el conjunto del
cerramiento transparente.
Para que un edificio pueda ser eficiente y
proporcionar confort a sus usuarios debe existir Investigaciones recientes demuestran que estos
una correcta relación entre su forma, orientación, indicadores, tradicionalmente utilizados, no son
materialidad y los factores ambientales externos, suficientes para evaluar el desempeño energético de
logrando que las estrategias de diseño sean los cerramientos transparentes. Stazi et al. (2014)
apropiadas al clima donde se apliquen. identifican la incidencia de algunos factores
claves en el comportamiento físico-térmico
En este sentido, la envolvente debe estar de las protecciones solares como la geometría
cuidadosamente diseñada de modo de controlar (inclinación, ancho y separación entre lamas),
los intercambios de energía, mereciendo materialidad (calor específico y conductividad),
especial atención los cerramientos transparentes y su colocación (interior y exterior). Dentro de
(conjunto conformado por: marco, vidrio y esta misma línea, Mainini et al. (2014) estudian
protección solar), los cuales pueden representar la efectividad de estos dispositivos en función
42
ENERLAC • Volumen IV. Número 2. Diciembre, 2020. ISSN: 2602-8042 (impreso) / 2631-2522 (digital).
OLADE – AUGM
Radiación incidente
FS
- Radiación absorvida
EXTERIOR INTERIOR
43
ANÁLISIS EXPERIMENTAL DE LA TRANSMISIÓN SOLAR DE DOS TIPOS DE PROTECCIONES SOLARES PARA DISTINTAS POSICIONES EN MONTEVIDEO, URUGUAY.
Pereira Ruchansky, Lucía; Gutiérrez Bazterrica, Lucía; Sosa Ibarra, Daniel; López Salgado, María Noel; Camacho Roberts, Magdalena
Figura 2. Características geométricas y textura de protecciones solares ensayadas: (a) PS1; (b) PS2
a) b)
44
ENERLAC • Volumen IV. Número 2. Diciembre, 2020. ISSN: 2602-8042 (impreso) / 2631-2522 (digital).
OLADE – AUGM
a) b)
Las protecciones se ensayaron para los planos ciones características de las fachadas de
verticales de orientación este, norte y noroeste los edificios. Adicionalmente, se realizaron
abarcando las horas en que el plano de estudio mediciones durante el mediodía solar con el
recibe los mayores niveles de radiación solar dispositivo de ensayo inclinado, orientando la
directa, de acuerdo a la distribución teórica protección solar normal a la radiación solar
de la radiación por plano para el solsticio de directa incidente y considerando una desvia-
verano en Montevideo, ver Figura 4. Con- ción de hasta 5 grados, según ANSI/NFRC
siderando que las fachadas este y oeste son 200-2014 [5]. En la Figura 5 se pueden
simétricas a la radiación, las mediciones observar las posiciones del dispositivo de ensayo
realizadas abarcan un amplio rango de posi- durante las mediciones.
45
ANÁLISIS EXPERIMENTAL DE LA TRANSMISIÓN SOLAR DE DOS TIPOS DE PROTECCIONES SOLARES PARA DISTINTAS POSICIONES EN MONTEVIDEO, URUGUAY.
Pereira Ruchansky, Lucía; Gutiérrez Bazterrica, Lucía; Sosa Ibarra, Daniel; López Salgado, María Noel; Camacho Roberts, Magdalena
Figura 4. Niveles de radiación solar incidente por plano (W/m2) para solsticio de verano en Montevideo
1000
PH
900
800
RADIACIÓN SOLAR (W/m)
700 O
E
600
SE SO
500
NE NO
400
300
200 N
S S
100
0
4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20
HORAS
Fuente: Cátedra Acondicionamiento térmico FADU-UdelaR, 2000
Figura 5. Dispositivo de ensayo de protecciones solares: (a) medición con PS1 en posición vertical;
(b) medición con PS2 en posición inclinada, normal a la radiación incidente
a) b)
46
ENERLAC • Volumen IV. Número 2. Diciembre, 2020. ISSN: 2602-8042 (impreso) / 2631-2522 (digital).
OLADE – AUGM
Las protecciones solares se ensayaron sin la un lado, las interacciones existentes entre el
utilización de vidrio, considerando que los vidrio y la protección solar que no permitirían
valores de factor solar, reflexión lumínica identificar el comportamiento aislado de la
y transmitancia térmica son conocidos y misma, además del efecto invernadero al
ensayados en laboratorio según norma interior de la caja que generaría la presencia de
ISO 15099:2003 para los distintos tipos de vidrio, afectando las mediciones y los resultados.
vidrios. La supresión del vidrio simplificó
las mediciones. De modo que se evitaron por RESULTADOS Y ANÁLISIS
47
ANÁLISIS EXPERIMENTAL DE LA TRANSMISIÓN SOLAR DE DOS TIPOS DE PROTECCIONES SOLARES PARA DISTINTAS POSICIONES EN MONTEVIDEO, URUGUAY.
Pereira Ruchansky, Lucía; Gutiérrez Bazterrica, Lucía; Sosa Ibarra, Daniel; López Salgado, María Noel; Camacho Roberts, Magdalena
Figura 6. Radiación solar incidente (W/m2), transmitida (W/m2) y transmisión solar (%) en PS1
para diferentes posiciones: (a) normal a la radiación solar directa incidente; (b) plano vertical
-orientación norte-
a) b)
1200 50% 1200 50%
45%
1000 1000
40% 40%
35% 35%
s (%)
800 800
s (%)
30% 30%
600 25% 600 25%
20% 20%
400 400
15% 15%
10% 10%
200 200
5% 5%
0 0% 0 0%
12:31
12:51
13:11
11:00
11:20
11:40
12:00
12:20
12:40
13:00
13:20
13:40
13:00
14:20
Radiación Incidente (W/m2) Radiación Trasmitida (W/m2) Radiación solar - TS (%)
Debido a la trayectoria del sol, que presenta orientación noroeste, Figura 7(b), con medi-
un comportamiento simétrico con eje en la ciones realizadas entre las 15:00h (H=58.6°,
hora del mediodía solar, se esperaba registrar Z=296.3°) y las 18:30h (H=16.2°, Z=257.6°),
un comportamiento simétrico en la radiación se puede observar un comportamiento más
incidente y en la radiación medida al interior, en estable, con variaciones durante toda la medición
la posición vertical orientada al norte. Los datos entre el 16.0% y el 18.3%.
obtenidos del piranómetro exterior muestran
esta simetría en su curva, sin embargo, se Si bien la tela screen no presenta a priori una
encuentra desplazada del eje correspondiente fuerte condicionante geométrica, se constató
al mediodía solar, lo que podría estar indicando que a medida que varía su posición relativa
que la ubicación del dispositivo de medición con respecto al sol, varía el porcentaje de
tiene una desviación con respecto al norte solar radiación transmitida. El valor medido en con-
o que la posición de los instrumentos presenta dición normal a la radiación directa incidente
una pequeña inclinación respecto a la protección (21.5%) se aproxima únicamente al registrado a
solar; provocando que los datos obtenidos no las 7:00h en la orientación este, cuando el azimut
reflejen esa simetría. relativo es menor a 20°, y representa el mayor
valor de transmisión solar registrado para esta
Para la PS1 en la orientación este, Figura 7(a), protección. Para el resto de las mediciones se
la transmisión solar es variable. Esta relación registraron valores que presentan una diferencia
tiende a disminuir a medida que avanza la respecto a la medición en posición normal de
medición, variando entre 20.5% al inicio hasta un 7.2% menos (diferencia entre el valor
(7:00h, H=12.0°, Z=105.9°), y un 14.3% al representativo en posición normal de 21.5% y
finalizar (12:00h, H=70.3°, Z=42.8°). En la el valor mínimo en orientación este de 14.3%).
48
ENERLAC • Volumen IV. Número 2. Diciembre, 2020. ISSN: 2602-8042 (impreso) / 2631-2522 (digital).
OLADE – AUGM
Figura 7. Radiación solar incidente (W/m2), transmitida (W/m²) y transmisión solar (%) en PS1 para
plano vertical en diferentes orientaciones: (a) este; (b) noroeste
a) b)
1200 50% 1200 50%
45%
4 45%
1000 1000
40%
4 40%
35% 35%
s (%)
800 800
s (%)
30% 30%
600 25% 600 25%
20% 20%
400 400
15% 15%
10% 10%
200 200
5% 5%
0 0% 0 0%
7:00
7:20
7:40
8:00
8:20
8:40
9:00
9:20
9:40
10:00
10:20
10:40
11:00
11:20
11:40
12:00
15:00
15:20
15:40
16:00
16:20
16:40
17:00
17:20
17:40
18:00
18:20
Radiación Incidente (W/m2) Radiación Trasmitida (W/m2) Radiación solar - TS (%)
En la Figura 8(a), correspondiente a la PS2 en res a los obtenidos en posición normal. Al igual
condición normal a la radiación incidente se que en las figuras anteriores, la línea vertical
identifica un comportamiento fluctuante de la corresponde al mediodía solar. La ausencia de
radiación transmitida, con picos y valles, que se un comportamiento simétrico en los resultados,
refleja en valores de transmisión solar que oscilan respecto a esta hora, es posible atribuirla a la
entre el 15.1% y el 24.2%. La transmisión solar ubicación relativa del piranómetro interior en
puede variar en un 8.9% al considerar un pico relación a la geometría de la protección; el cual
(24.2%, 12:31 h) y un valle (15.3%, 12:38 h), quedó desplazado del eje de simetría del trapecio
aún con un movimiento relativo del sol muy que conforma la chapa plegada. No obstante,
acotado. Para la orientación norte, Figura 8(b), se considera que los resultados obtenidos pre-
el porcentaje de radiación que pasa a través sentan una coherencia global en relación a lo
de la protección mantiene un comportamiento esperado desde el punto de vista teórico.
fluctuante, con valores máximos del 14.7%
(11:29h, H=60.4°, Z=48.8°) y mínimos del
3.4%, 16 minutos después (11:45h, H=62.7°,
Z=42.3°). La mediciones de las Figura 8(a) y
8(b) presentan el mismo azimut, variando la
altura y la extensión del periodo de medición,
sin embargo, se puede observar que los valores
obtenidos para la posición vertical son inferio-
49
ANÁLISIS EXPERIMENTAL DE LA TRANSMISIÓN SOLAR DE DOS TIPOS DE PROTECCIONES SOLARES PARA DISTINTAS POSICIONES EN MONTEVIDEO, URUGUAY.
Pereira Ruchansky, Lucía; Gutiérrez Bazterrica, Lucía; Sosa Ibarra, Daniel; López Salgado, María Noel; Camacho Roberts, Magdalena
Figura 8. Radiación solar incidente (W/m2), transmitida (W/m2) y transmisión solar (%) de PS2
para diferentes posiciones: (a) normal a la radiación solar directa incidente; (b) plano vertical
-orientación norte-
a) b)
1200 50%
45%
45%
1000
40% 1000
40%
35% 35%
s (%)
800 800
s (%)
30% 30%
600 25% 600 25%
20% 20%
400 400
15% 15%
10% 10%
200 200
5% 5%
0 0% 0 11:00 0%
11:30
12:00
12:30
13:00
13:30
14:00
14:30
12:21
12:31
12:41
12:51
13:01
b)
Radiación Incidente (W/m2) Radiación Trasmitida (W/m2) Radiación solar - TS (%)
50
ENERLAC • Volumen IV. Número 2. Diciembre, 2020. ISSN: 2602-8042 (impreso) / 2631-2522 (digital).
OLADE – AUGM
Figura 9. Radiación solar incidente (W/m2), transmitida (W/m2) y transmisión solar (%) de PS2 para
plano vertical en diferentes orientaciones: (a) este; (b) noroeste
a) b)
1200 50% 1200 50%
45% 45%
1000 40% 1000
40%
35% 35%
s (%)
800 800
s (%)
30% 30%
0 0% 0 0%
13:00
13:20
13:40
14:00
14:20
14:40
15:00
15:20
15:40
16:00
16:20
16:40
17:00
17:20
17:40
18:00
18:20
7:00
7:20
7:40
8:00
8:20
8:40
9:00
9:20
9:40
10:00
10:20
10:40
11:00
11:20
11:40
12:00
Los resultados para la PS2, en todos los casos El plegado trapezoidal de la chapa por su parte,
presentan un comportamiento fluctuante en genera un auto-sombreado, mientras que las
magnitud y amplitud, lo cual permite constatar perforaciones presentes en la misma permiten
que las variaciones en la transmisión solar están mayor o menor captación de radiación directa
condicionadas por su geometría. dependiendo del ángulo de incidencia, ver Figura
10(a) y (b).
Figura 10. Detalle de incidencia de radiación solar directa sobre protección “Chapa perforada
trapezoidal” y presentación esquemática
MÁXIMA
MÁXIMA
CAPTACIÓN CAPTACIÓN
51
ANÁLISIS EXPERIMENTAL DE LA TRANSMISIÓN SOLAR DE DOS TIPOS DE PROTECCIONES SOLARES PARA DISTINTAS POSICIONES EN MONTEVIDEO, URUGUAY.
Pereira Ruchansky, Lucía; Gutiérrez Bazterrica, Lucía; Sosa Ibarra, Daniel; López Salgado, María Noel; Camacho Roberts, Magdalena
Análisis comparativo
Figura 11. Transmisión solar (%) de PS1 y PS2 para todas las posiciones y orientaciones
25%
20%
15%
S (%)
10%
5%
0%
7:00
7:20
7:40
8:00
8:20
8:40
9:00
9:20
9:40
10:00
10:20
10:40
11:00
11:20
11:40
12:00
12:20
12:40
13:00
13:20
13:40
14:00
14:20
14:40
15:00
15:20
15:40
16:00
16:20
16:40
17:00
17:20
17:40
18:00
18:20
Estos datos permiten afirmar que, aun protec- ciones para todas las orientaciones.
ciones de geometría homogénea -PS1- presentan La PS1 presenta un rango de transmisión solar
desempeños variables según el ángulo de inci- más acotado (14.0% a 21.0%), mientras que
dencia de la radiación solar. Mientras que, para la PS2 presenta un rango más extenso (3.4% a
el caso de la PS2, se constata la incidencia de 25.0%), ver Figura 12. De esto se desprende
la geometría compleja en su desempeño como que, para todas las orientaciones y posiciones,
protección solar, registrándose resultados aún el valor medio de transmitancia solar es más
más variables. En la Figura 12, se puede ob- representativo para la PS1 que para la PS2,
servar la caracterización de los valores de debido a que en la primera los resultados
transmisión solar para cada protección solar presentan una menor desviación respecto a
si consideramos los resultados de las medi- la media.
52
ENERLAC • Volumen IV. Número 2. Diciembre, 2020. ISSN: 2602-8042 (impreso) / 2631-2522 (digital).
OLADE – AUGM
Figura 12. Distribución normal de coeficientes de transmisión solar para las protecciones PS1 y PS2
Caracterización del desempeño de las Figura 13. Lectura de ángulos de incidencia para
protecciones el plano de medición
La Figura 14 presenta los resultados obtenidos niveles de transmitancia (20% o más); así
para la PS1 donde se observa que en ángulos mismo los menores valores (entre 14% y 17%)
cercanos a la perpendicular, tanto en azimut se presentan únicamente a partir de los 55°
como en altura, se presentan los mayores de altura.
53
ANÁLISIS EXPERIMENTAL DE LA TRANSMISIÓN SOLAR DE DOS TIPOS DE PROTECCIONES SOLARES PARA DISTINTAS POSICIONES EN MONTEVIDEO, URUGUAY.
Pereira Ruchansky, Lucía; Gutiérrez Bazterrica, Lucía; Sosa Ibarra, Daniel; López Salgado, María Noel; Camacho Roberts, Magdalena
Figura 14. Transmisión solar (%) de PS1 en función de posición relativa del sol
90̄
80̄
70̄
60̄
ALTURA
50̄
40̄
30̄
20̄
10̄
0̄
0̄
90̄
80̄
70̄
60̄
50̄
40̄
30̄
20̄
10̄
10̄
20̄
30̄
40̄
50̄
60̄
70̄
80̄
90̄
AZIMUT RELATIVO
14 a 17% 17 a 20% Mayor a 20%
Figura 15. Transmisión solar (%) de PS2 en función de posición relativa del sol
90̄
80̄
70̄
60̄
ALTURA
50̄
40̄
30̄
20̄
10̄
0̄
90̄
80̄
70̄
60̄
50̄
40̄
30̄
20̄
10̄
0̄
10̄
20̄
30̄
40̄
50̄
60̄
70̄
80̄
90̄
AZIMUT RELATIVO
Menor a 8% 8% a 11% 11% a 14% 14% a 17% 17% a 20% Mayor a 20%
Fuente: Elaboración propia
54
ENERLAC • Volumen IV. Número 2. Diciembre, 2020. ISSN: 2602-8042 (impreso) / 2631-2522 (digital).
OLADE – AUGM
Para la PS2, Figura 15, es posible identificar diferencia podría, por ejemplo, llevar a un sobre
sobre los 50.0° de altura un cambio en los rangos dimensionamiento de los cálculos de demanda
de transmitancia independiente del azimut. de refrigeración o condicionar los niveles de
Mientras que, para alturas menores a 50.0° el iluminación interior, poniendo en discusión la
rango oscila entre 14.0% y 20.0%, para alturas utilidad del factor solar como herramienta de
mayores la transmisión varía entre 3.4% y 17.0%. diseño.
A su vez, dentro de cada zona se identifica que El estudio realizado permitió caracterizar el
los niveles más altos se presentan en diferentes desempeño de las protecciones estudiadas para
ángulos de incidencia – altura y azimut relativo al posibles orientaciones de fachadas de edificios
plano de la protección -, lo cual tiene relación con en Montevideo, Uruguay. Adicionalmente, se
el comportamiento presentado en la Figura 10. presentó un análisis preliminar que permite
caracterizar el desempeño de las protecciones
CONCLUSIONES solares para cualquier ubicación y orientación,
a través del análisis realizado que expresa la
Este trabajo presenta los estudios realizados transmisión solar en función del ángulo de
sobre dos tipos de protecciones solares, cortina radiación incidente. Para profundizar el estudio
de tela screen y chapa perforada trapezoidal, de las protecciones solares, es necesario un
para evaluar la transmisión de radiación solar mayor número de mediciones frente a otros
de onda corta a través de ellas, en diferentes ángulos de incidencia de radiación para obte-
posiciones y orientaciones. En las mediciones ner una caracterización global por protección
realizadas se pudo constatar que el desempeño solar, así como para independizar esta variable
de las protecciones solares es variable según de la incidencia de la geometría que la carac-
el ángulo de la radiación incidente, o sea, a su teriza. El desarrollo de esta metodología podría
posición relativa al sol, y la geometría propia de aportar un indicador más adecuado para la
la protección. toma de decisiones en el diseño y evaluación de
cerramientos transparentes.
Los resultados obtenidos de forma experimen-
tal para la ciudad de Montevideo, demuestran Los datos de este artículo forman parte del
que las protecciones ensayadas presentan una proyecto Comisión Sectorial de Investigación
diferencia en su desempeño según la posición Científica (CSIC), I+D “Efectos de las
en que se orienten, lo cual revela la dificultad características de las ventanas, los vidrios y
de utilizar un indicador único por dispositivo los dispositivos de protección solar sobre los
de protección solar. Considerando siempre su requerimientos energéticos del edificio y las
evaluación para el período de verano, ambas condiciones de confort térmico y visual”.
protecciones presentan un mejor desempeño,
es decir menores niveles de transmisión solar,
cuando se orientan hacia el norte. •••
55
ANÁLISIS EXPERIMENTAL DE LA TRANSMISIÓN SOLAR DE DOS TIPOS DE PROTECCIONES SOLARES PARA DISTINTAS POSICIONES EN MONTEVIDEO, URUGUAY.
Pereira Ruchansky, Lucía; Gutiérrez Bazterrica, Lucía; Sosa Ibarra, Daniel; López Salgado, María Noel; Camacho Roberts, Magdalena
REFERENCIAS
56
ISSN 2602-8042 [Impresa]
ISSN 2631-2522 [Electrónica]
EVALUACIÓN TÉRMICA Y EFECTOS DE LA
DEMANDA EN UNA INSTALACIÓN SOLAR DE UN
PROCESO PRODUCTIVO
THERMAL ASSESSMENT AND EFFECTS OF DEMAND ON A
SOLAR INSTALLATION OF A PRODUCTIVE PROCESS
58
ENERLAC • Volumen IV. Número 2. Diciembre, 2020. ISSN: 2602-8042 (impreso) / 2631-2522 (digital).
OLADE – AUGM
59
EVALUACIÓN TÉRMICA Y EFECTOS DE LA DEMANDA EN UNA INSTALACIÓN SOLAR DE UN PROCESO PRODUCTIVO.
Albizzati, Enrique Donato
ABSTRACT INTRODUCCIÓN
60
ENERLAC • Volumen IV. Número 2. Diciembre, 2020. ISSN: 2602-8042 (impreso) / 2631-2522 (digital).
OLADE – AUGM
al rango entre 150-400°C, y el 48% supera conveniente como primer paso analizar su
los 400°C de temperatura (Kumar et al., 2019; performance para satisfacer diferentes deman-
Solar Payback, 2017). En principio se han das térmicas en aplicaciones relacionadas con el
priorizado las aplicaciones de la energía solar sector productivo de la región latinoamericana.
en la producción de agua caliente para lavar y Dado que el fin es mejorar la productividad
limpiar, calentamiento de agua para genera- y competitividad industrial, en posteriores
ción de vapor, calor para procesos en baños y etapas será necesario evaluar ambiental y
recipientes, y secado por convección con aire económicamente la instalación solar propuesta.
caliente (IEA-ETSAP e IRENA, 2015; Lauterbach
et al., 2014). En este trabajo se simula una instalación solar
de mediana escala para calentamiento de agua
Particularmente en la provincia de Santa Fe, que abastece un proceso productivo cuya de-
Argentina se encuentran radicadas numerosas manda térmica es variable durante el día y a
empresas pequeñas y medianas relacionadas con lo largo del año. El agua caliente suministrada
el sector agropecuario y agroindustrial. En las se emplea como fluido de intercambio térmico
mismas se fabrican equipamientos y producen en el equipo donde se lleva a cabo el proceso.
alimentos, a la vez que se ejecutan programas Con el software System Advisor Model (SAM
para fortalecer la gestión de la energía y promover versión 2020.2.29) se determinan la radiación
la incorporación de fuentes renovables, siendo solar incidente y aprovechada, la temperatura
por ello un campo donde puede aprovecharse la del agua suministrada, la energía termosolar
energía solar térmica. aportada y los consumos energéticos comple-
mentarios. Luego se analiza el efecto de la
Por otra parte se dispone de técnicas de mode- demanda térmica sobre la performance de la
lado y simulación en programas como TRNSYS, instalación considerando casos donde solo
T-SOL, WATSUN y Polysun, para predecir el aumenta el caudal de agua caliente, y otros
comportamiento de sistemas solares térmicos donde además se necesita energía adicional
en diferentes condiciones climáticas, mejorar para mantenimiento térmico cuando no se lleva
el rendimiento de los elementos componentes y a cabo el proceso. Este suministro de calor para
analizar su sensibilidad frente a cambios en los mantenimiento evita que la carga contenida
parámetros de operación. En este marco muchas en el equipo se enfríe, ya que permanece a una
publicaciones han presentado estudios sobre el temperatura prefijada.
uso de la energía termosolar en áreas productivas
empleando variadas herramientas de cálculo y
tienen el respaldo de software con confiabilidad MATERIALES Y MÉTODO
comprobada (Abdel-Dayem et al., 2001; Frein et
al., 2014; Kalogirou, 2003; Quijera et al., 2011; Integración de la energía solar térmica
Suresh y Rao, 2017). En otros artículos se ha
incluido la determinación de los efectos de la En la Tabla 1 se muestran procesos y servicios
variación en el tiempo de la demanda térmica industriales donde se aplica la energía solar
sobre la performance de instalaciones solares térmica, observándose que los mismos abarcan
industriales (Mousa y Taylor, 2016; Fenercom, varios sectores. Las industrias de alimentos y
2014; Lauterbach et al., 2014). bebidas, fabricación de materiales o piezas metá-
licas, procesos químicos y plantas para lavado
Debe tenerse en cuenta que es factible y limpieza son las que poseen mayor potencial
construir plantas solares modulares estándar (Farjana, 2018; Kumar et al., 2019; Sharma
para ofrecerlas en el mercado, por lo cual es et al., 2017).
61
EVALUACIÓN TÉRMICA Y EFECTOS DE LA DEMANDA EN UNA INSTALACIÓN SOLAR DE UN PROCESO PRODUCTIVO.
Albizzati, Enrique Donato
La aplicación del calor solar está en constante porado 567 MW térmicos. En ese mismo año
crecimiento y se han implementado plantas 108 plantas fueron proyectadas, y a los países
que van desde la pequeña escala hasta grandes que figuran en la Tabla 2 se deben agregar
sistemas que producen 100 MW térmicos. Argentina, Chipre, Grecia, Israel, Italia, Corea
Según datos publicados, al menos 741 plantas del Sur, Sudáfrica, Turquía, la Unión Europea
con 662,648 m2 de colectores fueron operadas y Estados Unidos (Weiss y Spörk- Dür, 2019).
en todo el mundo en 2018, habiéndose incor-
62
ENERLAC • Volumen IV. Número 2. Diciembre, 2020. ISSN: 2602-8042 (impreso) / 2631-2522 (digital).
OLADE – AUGM
Una vez identificado el proceso o servicio in- tienen mayor durabilidad y los últimos generan
dustrial, se requiere integrar la energía solar energía térmica más eficientemente. Los procesos
aplicando una metodología para dimensionar y servicios industriales en los que ya se usa la
y configurar las instalaciones en función de las radiación solar específicamente para calentar
necesidades propias de la producción y el uso de agua son el precalentamiento de la alimentación
las otras fuentes de energía complementarias. a calderas, higiene de equipos e instalaciones
Por ello se están adecuando continuamente pro- y calor, para operaciones en tinas, baños y
cedimientos, según la localización geográfica y recipientes (Farjana, 2018; Sharma et al., 2017).
el tipo de aplicación que se planifica. Se pueden
contemplar las industrias que presentan el mayor Dado que el calentamiento solar de agua puede
consumo de energía, así como casos en que se abastecer procesos térmicos donde se opera du-
dispone de fluidos que retornan de un proceso rante un tiempo predeterminado, se pueden
o servicio a una temperatura relativamente aprovechar corrientes de retorno y el nivel de
alta. También se debe tener en cuenta que las temperatura no es muy elevado, en el caso elegido
tecnologías más confiables son las desarrolladas como base para este estudio el proceso se lleva
para bajas temperaturas y utilizan colectores a cabo durante seis horas al día y cinco días de
solares de placa plana y de tubos evacuados, la semana, excepto en la primer quincena de
habiéndose ya reconocido que los primeros enero (Figura 1).
100
50
0
1 2 3 4 5 6 7 8 9 10 11 12 1 2 3 4 5 6 7 1 5 9 13 18 23
La instalación solar estará localizada cerca de entre éste y el tanque acumulador. El circuito
la ciudad de Santa Fe, Argentina. En la misma hidráulico que posee bombas, y un controlador
el agua circula por convección forzada entre que gobierna la captación de la radiación solar
los colectores solares y el intercambiador, y (Figura 2).
63
EVALUACIÓN TÉRMICA Y EFECTOS DE LA DEMANDA EN UNA INSTALACIÓN SOLAR DE UN PROCESO PRODUCTIVO.
Albizzati, Enrique Donato
Tanque
Colectores Intercambiador acumulado
solares de calor
El controlador posibilita que el agua circule Los colectores solares adoptados son de placa
por el circuito de captación de la radiación plana. En los mismos una parte de la radiación
solar, solo si la temperatura de salida del agua solar se transmite por una cubierta de vidrio
de los colectores es varios grados mayor que semitransparente, se absorbe en una placa de
la temperatura del agua en el fondo del tanque aluminio recubierta con un material selectivo, y
acumulador. El sistema posee dos fuentes al convertirse en energía térmica se transfiere al
de energía complementarias: suministro de agua que circula por tubos de cobre. Entre la placa
energía auxiliar con gas natural que completa y la cubierta hay una cámara de aire, se utiliza un
los requerimientos térmicos, y electricidad para aislante térmico para evitar las pérdidas de calor,
bombeo en el circuito hidráulico. y una caja de chapa metálica contiene todos los
componentes (Albizzati, 2012; Kalogirou, 2014).
La temperatura del proceso es 45°C, la
temperatura fijada para suministro del agua El sistema modular está configurado con cuatro
caliente es 75°C, la temperatura de la corriente baterías de cinco colectores cada una, todas en
de retorno 50°C, y el caudal circulante por el paralelo (Figura 3).
equipo de proceso 260 kg/h. El agua caliente
circula por un intercambiador de calor que Figura 3. Configuración de los colectores solares
alimenta térmicamente al equipo donde trans-
curre el proceso, siendo la temperatura de la AGUA
corriente de salida 5°C más alta que la necesa- DEL
TANQUE
ria en el mismo. Si el sistema solar suministra
agua a menos de 75°C, esta temperatura se
alcanza mediante la fuente térmica auxiliar.
Cuando la temperatura de la corriente que
sale del tanque es mayor a 75°C la corriente se
mezcla con una parte de la corriente de retorno
para lograr el nivel térmico fijado. Si even-
tualmente la temperatura del tanque supera los
100 °C, se descarga calor para evitar los pro- AGUA
DEL
blemas originados por el sobrecalentamiento. TANQUE
64
ENERLAC • Volumen IV. Número 2. Diciembre, 2020. ISSN: 2602-8042 (impreso) / 2631-2522 (digital).
OLADE – AUGM
La energía solar captada por unidad de tiempo (qc) Para calcular la radiación solar incidente sobre
se computa con los parámetros característicos, los colectores y luego evaluar la performance
siendo la diferencia entre la energía absorbida y del sistema solar se emplea el software SAM.
la perdida en los colectores. Intervienen también Este simulador fue desarrollado por el National
en el cálculo la radiación total incidente (IT), el Renewable Energy Laboratory (USA) para la
área neta total (Ac), la temperatura de la corriente planificación y evaluación de proyectos de ener-
que ingresa (Te) y la temperatura ambiente (Ta) gías renovables, así como para la investigación
(Duffie y Beckman, 2013) académica y ha sido verificado mediante la
aplicación de otros modelos más complejos
qc = (FR EO ) IT Ac Nc − (FR Uc ) Ac Nc (Te − Ta ) (1) como TRNSYS (DiOrio et al., 2014).
65
EVALUACIÓN TÉRMICA Y EFECTOS DE LA DEMANDA EN UNA INSTALACIÓN SOLAR DE UN PROCESO PRODUCTIVO.
Albizzati, Enrique Donato
El programa para calentamiento de agua permite el año, seleccionando en forma apropiada las
modificar la localización geográfica, los perfiles entradas y salidas.
de la demanda, las características del colector,
y muchas otras variables. Para su utilización se SAM tiene varias opciones para hallar la radiación
deben cargar valores horarios de la radiación recibida por los colectores solares, partiendo de la
solar en el plano horizontal, y datos climáticos radiación solar global (I), la radiación directa (ID)
correspondientes al Typical Meteorological Year y la radiación difusa (Id) en el plano horizontal.
del lugar donde se encuentra la instalación. El En este caso se calcula con el modelo anisotrópico
software es dinámico y genera 8,760 valores HDRK que contempla el índice de anisotropía (Ai),
horarios, relacionados fundamentalmente con el cociente entre la radiación directa en planos
el aporte de energía térmica de origen solar y la inclinados y el horizontal (RD), el ángulo de
temperatura del agua. Se puede obtener también inclinación de los colectores (β) y la reflectividad
otra información horaria, mensual y para todo de los alrededores (r) (Duffie y Beckman, 2013):
IT = (ID + Id Ai) RD + Id (1 − Ai) [(1+ cos β)/2] [1 +(ID/I)0,5 . sen3 (β/2)] + I r [(1− cos β)/2] (2)
El modelo físico-matemático aplicado supone que intercambiador y egresa otra corriente más fría
hay dos situaciones claves en el funcionamiento hacia ese equipo (Figura 4). La energía ganada
del sistema solar. En una de ellas se produce es menor que la energía colectada debido a las
la captación de la radiación y ganancia de pérdidas en las cañerías y a que la eficiencia del
energía en el tanque perfectamente mezclado, intercambiado de calor es menor al 100 %.
ya que ingresa la corriente que proviene del
Figura 4. Modelos del tanque acumulador: con mezclado perfecto (a) y estratificado en dos zonas (b)
Energía
Energía perdida
Energía Zona
caliente acumulada
Energía Energía de reposición Energía
ganada acumulada de reposición
Energía
perdida Zona Energía
fría Energía
acumulada
perdida
(a) (b)
66
ENERLAC • Volumen IV. Número 2. Diciembre, 2020. ISSN: 2602-8042 (impreso) / 2631-2522 (digital).
OLADE – AUGM
Las pérdidas de energía térmica se calculan con el Se debe tener en cuenta que también se consume
coeficiente total de transferencia de calor (Ut) y el energía eléctrica en las bombas (Pb), y entonces la
área (At) del tanque: fracción solar efectiva es:
Qd = M Cp (Tc – Tr ) (6)
Fs = Qe / Qd (8)
67
EVALUACIÓN TÉRMICA Y EFECTOS DE LA DEMANDA EN UNA INSTALACIÓN SOLAR DE UN PROCESO PRODUCTIVO.
Albizzati, Enrique Donato
Los procesos industriales pueden requerir sola- da una energía adicional para mantenimiento,
mente del suministro tipo A o del suministro tipo que es el 15% de la originalmente fijada para el
B, y en algunos casos es posible que se pueda suministro tipo A. Se cubren de este modo las
elegir entre las dos alternativas. Debe observarse pérdidas de calor del equipo en momentos en
que en los casos II-IV no se modifica el perfil de que no se lleva a cabo el proceso (Figura 5). La
la demanda propuesto en el caso I, se requieren operación implica que aunque no se modifiquen
diferentes caudales y se opera durante seis horas. el caudal y el tiempo del proceso, la energía
En los suministros tipo B (casos V-VIII) se deman- anual que se debe proveer es mayor.
100
DEMANDA TÉRMICA (%)
50
0
1 2 3 4 5 6 7 8 9 10 11 12 1 2 3 4 5 6 7 1 5 9 13 18 23
68
ENERLAC • Volumen IV. Número 2. Diciembre, 2020. ISSN: 2602-8042 (impreso) / 2631-2522 (digital).
OLADE – AUGM
800
RAD SOLAR (W/m)
600
TOTAL
400
DIRECTA
200
DIFUSA
0
0 4 8 12 16 20 24
HORA DEL DÍA
Fuente: Elaboración propia
Siendo la radiación total incidente anual invierno. En los meses fríos la radiación solar
1,899.9 kWh/m2 y la demanda térmica anual igual recibida por los colectores disminuye sensible-
a 11,326.3 kWh, al aplicar el software se obtiene mente, y la energía térmica ganada que llega al
que la fracción solar es 0.835 y la eficiencia térmica tanque es mucho menor que en el resto del año.
de la instalación es 0.124. Se encuentra también
que como son necesarios 1,870.5 kWh de energía
para completar el suministro de calor y 517.4 kWh
son consumidos para bombear el fluido, la frac-
ción solar efectiva es 0.789.
69
EVALUACIÓN TÉRMICA Y EFECTOS DE LA DEMANDA EN UNA INSTALACIÓN SOLAR DE UN PROCESO PRODUCTIVO.
Albizzati, Enrique Donato
90
80
T TAN (C)
70
60
50
40
ENE FEB MAR ABR MAY JUN JUL AGO SEP OCT NOV DIC
900
Q (kWh)
600
300
0
ENE FEB MAR ABR MAY JUN JUL AGO SEP OCT NOV DIC
Fuente: Elaboración propia
70
ENERLAC • Volumen IV. Número 2. Diciembre, 2020. ISSN: 2602-8042 (impreso) / 2631-2522 (digital).
OLADE – AUGM
Al determinar anualmente el efecto de la deman- La Tabla 6 muestra que los mayores valores de
da sobre la performance de la instalación la fracción solar son 0.835 para suministros A
modular, los cálculos indican que con suministro (caso I) y 0.670 para suministros B (caso V). Las
tipo A (casos I-V) al crecer la misma se reduce eficiencias térmicas máximas son 0.322 y 0.267
la fracción de energía solar, pero aumentan la para los casos IV y VIII respectivamente.
energía termosolar entregada y la eficiencia. Lo
mismo sucede al analizar los suministros tipo B
(casos V-VIII).
El desfasaje entre el inicio del proceso y la casos. En V-VIII, se observa que al multiplicarse
ganancia real de energía provoca que exista un la demanda baja la fracción solar, creciendo la
consumo de energía auxiliar a partir de las 7 de eficiencia y la energía termosolar que se entrega.
la mañana, que es relativamente más importante
cuando crece el caudal que se suministra. Al Otros efectos se ven en la referida tabla
comparar los casos I y II se comprueba que cuando no cambia el caudal de agua caliente
aunque se duplica la demanda, el aumento de la y se consume energía para mantenimiento
energía termosolar suministrada es del 76 %, por térmico, aunque siempre hay disminución de
lo que la energía térmica auxiliar prácticamente la fracción solar. Si por ejemplo se comparan
se triplica. El suministro crece porque la energía los casos I y V el suministro de energía y la
ganada en el sistema es mayor, la temperatura en eficiencia anual aumentan, mientras que esos
el tanque es menor, se capta más radiación solar indicadores disminuyen al comparar los casos
y las pérdidas de calor disminuyen. III y VII. Entonces existe un valor de la demanda
a partir del cual se tiene una menor eficiencia
Efectos similares se producen cuando se al incorporar mantenimiento térmico a la
aumenta el caudal del agua en los restantes instalación (Figuras 9 y 10).
71
EVALUACIÓN TÉRMICA Y EFECTOS DE LA DEMANDA EN UNA INSTALACIÓN SOLAR DE UN PROCESO PRODUCTIVO.
Albizzati, Enrique Donato
1000
800
energía termosolar (kWh/m2)
A
Suuministro anual de
600
400
B
200
Demanda anual de
energía térmica(kWh/m2)
Fuente: Elaboración propia
0,8
A
Fracción anual de energía
0,6
termosolar
0,4
B
0,2
0,0
72
ENERLAC • Volumen IV. Número 2. Diciembre, 2020. ISSN: 2602-8042 (impreso) / 2631-2522 (digital).
OLADE – AUGM
Las figuras anteriores muestran que si no hay Los indicadores calculados en los casos tipo
mantenimiento térmico la energía termosolar A son inferiores a los obtenidos en Madrid,
anual máxima entregada es 611.5 kWh/m2 de España para instalaciones industriales cuya
colector y la fracción solar mínima es 0.540, pero demanda energética diaria promedio por uni-
cuando se requiere mantenimiento dichos valo- dad de área de colección se encuentra en el
res son 506.8 kWh/m2 de colector y 0.259. La mismo rango. En ellas ha verificado que al calen-
energía suministrada aumenta un 38% entre I y V, tar agua desde 15°C hasta a 60°C en forma
mientras que disminuye un 17% entre IV y VIII. continua durante todo el año, se obtiene una
alta eficiencia y suministros anuales entre 700
Cuando se requiere mantenimiento térmico si y 1.000 kWh/m2 de colector, notándose que
crece el caudal de fluido suministrado, aunque la menor temperatura del agua en el sistema
la ganancia de energía termosolar es mayor, provoca un crecimiento de esos valores
el consumo de energía auxiliar al inicio del (Fenercom, 2010).
proceso es más importante. Por esto en los casos
VII y VIII donde la demanda es más elevada, Para los suministros tipo B los indicadores
la eficiencia térmica es más baja que en III y IV se pueden comparar con los de instalaciones
respectivamente. de igual demanda energética diaria promedio
por unidad de área de colección localizadas en
Los resultados obtenidos muestran también en Madrid, España. En ellas el baño está a 45°C
los casos con suministro tipo A, que la mayor y opera doce horas por día, excepto sábado y
demanda disminuye el cociente entre la potencia domingo (Fenercom, 2010). Los colectores pla-
necesaria para bombeo y la energía solar térmica nos están inclinados a 35° mirando al sur, la
suministrada. Una tendencia similar se verifica acumulación es de 0.070 m3/m2, la temperatura
en los casos tipo B, aunque en cada instalación de la corriente de suministro es 70°C, la del
el mantenimiento térmico provoca un consumo retorno es 50°C y el mantenimiento es el 20%
mayor de energía eléctrica. del consumo en el proceso (Tabla 7).
En los casos F1-F4 que están en la citada tabla 1,769.0 kWh/m2, casi un 7% menor que en Santa
se alcanzan valores más elevados de todos los Fe. Debe tenerse en cuenta que hay diferencias
indicadores que en los casos estudiados V-VIII, en el modo de operación, y que en la performance
aunque la radiación solar anual sobre Madrid es influyen la distribución de la radiación durante
73
EVALUACIÓN TÉRMICA Y EFECTOS DE LA DEMANDA EN UNA INSTALACIÓN SOLAR DE UN PROCESO PRODUCTIVO.
Albizzati, Enrique Donato
el año, los parámetros de los colectores, el perfil por año los resultados indican que la fracción
diario de consumo del agua caliente, la tem- de energía solar es 0.835 y la eficiencia térmica
peratura y el salto térmico del fluido requerido en es 0.124, siendo la fracción efectiva es 0.789 al
el proceso. considerar el consumo de energía eléctrica para
bombeo.
La energía térmica solar anual suministrada
en las instalaciones españolas alcanza el valor Luego se comprueba que al cambiar la demanda
máximo de 615.0 kWh/m2 de colector que debido solo a aumentos del caudal de agua,
supera en un 20% al máximo calculado para tanto en casos con mantenimiento térmico como
Santa Fe, ambos para la misma demanda por sin mantenimiento, se reduce la fracción de
unidad de área de colector. Puede decirse que energía solar, pero crecen la energía termosolar
un efecto muy significativo produce la reducción suministrada y la eficiencia.
en 5°C del salto térmico del fluido calefactor,
pues porcentualmente su caída es similar al Asimismo, al adicionar el mantenimiento a la
crecimiento del suministro. instalación en todos los casos disminuye la
fracción solar, pero para baja demanda térmica
Luego, para el sistema que se ha estudiado se el suministro y eficiencia anuales aumentan,
han calculado rendimientos térmicos similares mientras que para una más alta demanda esos
a aquellos encontrados en lugares donde se indicadores disminuyen. Finalmente, se estima
están realizando intensos programas para la que los indicadores de la performance posibilitan
incorporación de sistemas solares térmicos a las instalaciones abastecer satisfactoriamente
en el sector productivo. Sin embargo, para un diferentes demandas térmicas. Entonces, el
mejor aprovechamiento local deberá avanzarse aprovechamiento de la energía solar térmica es
principalmente en lo que respecta a la disponi- un recurso a tener en cuenta para contribuir a
bilidad de tecnologías para la captación de la mejorar la productividad y competitividad de
radiación solar y recursos humanos capacitados pequeñas y medianas empresas industriales que
para el diseño y montaje, además de los aspectos están radicadas en la región central de Argentina.
ambientales y económicos involucrados en las
instalaciones. •••
CONCLUSIONES AGRADECIMIENTOS
El software System Advisor Model (SAM) posibilita Se agradece a la Facultad de Ingeniería Química
analizar detalladamente el comportamiento de de la Universidad Nacional del Litoral por el apoyo
instalaciones solares para calentar agua. Los recibido para realizar este trabajo.
datos horarios de la radiación solar, variables
climáticas locales, cantidad y temperatura del
agua requerida, características de los colectores
adoptados y del tanque acumulador, se encuen- REFERENCIAS
tran entre los entre los más importantes para
posibilitar el uso del programa. Abdel-Dayem, A. M., Mohamad, M. (2001). Potential of
solar energy utilization in the textile industry-a case study.
Al simular el sistema solar para calentar agua Renewable Energy, 23, 3-4, 685-94.
que abastece un proceso productivo que tiene
20 colectores planos selectivos con 40 m2 de Albizzati, E. (2012). Energía solar térmica: principios básicos
superficie de captación y un tanque acumulador y aplicaciones tecnológicas. Centro de Publicaciones de la
de 3.0 m3, cuando la demanda es de 11,326.3 kWh Universidad Nacional del Litoral, Santa Fe, Argentina.
74
ENERLAC • Volumen IV. Número 2. Diciembre, 2020. ISSN: 2602-8042 (impreso) / 2631-2522 (digital).
OLADE – AUGM
ANSI/ASHRAE 93. (2010). Methods of testing to determine Kalogirou, S. (2014). Solar energy engineering: processes
the thermal performance of solar collectors. American and systems. 2nd ed., Elsevier/Academic Press, Burlington,
National Standards Institute /American Society of Heating, MA, USA.
Refrigerating and Air Conditioning Engineers.
Kumar, L., Hasanuzzaman, M., Rahim, N. A. (2019). Global
Christensen, C., Maguire, J., Burch J., DiOrio, N. (2014). advancement of solar thermal energy technologies for
Simplified solar water heater simulation using a multi-mode industrial process heat and its future prospects: A review.
tank model. https://sam.nrel.gov/images/web_page_files/ Energy Conversion and Management, 195, 885-908
christensen-2014-simplified-swh-simulation-multi-mode-
tank.pdf. Lauterbach, C., Schmitt, B., Vajen, K. (2014). System analysis
of a low-temperature solar process heat system. Solar
DiOrio, N., Christensen, C., Burch, J., Dobos, A. (2014). Energy, 101, 117–130.
Technical Manual for the SAM Solar Water Heating Model.
https://sam.nrel.gov/resources. Mousa, O. B., Taylor, R. (2016). Solar thermal sterilization:
a TRNSYS performance analysis. Asia Pacific Solar Research
Duffie, J. A., Beckman, W. A. (2013). Solar engineering of Conference. http://apvi.org.au/solar-research-conference/
thermal processes. 4th ed., J. Wiley & Sons, Hoboken, NJ, USA. wp-content/uploads/2017/02/O-Bany-Mousa-and-R-
Taylor-Solar-Thermal-Sterilization-A-TRNSYS-Performance-
EnergyPlus. (2017). EnergyPlus Weather Data. https:// Analysis.pdf.
energyplus.net/weather.
Quijera, J. A., González Alriols, M., Labidi, J. (2011). Usage of
Farjana, S. H., Huda, N., Mahmud, M. A. P., Saidur, R. (2018). solar energy in an industrial process. Chemical Engineering
Solar process heat in industrial systems-A global review. Transactions, 25, 875–80.
Renewable and Sustainable Energy Reviews, 82, Part 3,
February, 2270-2286. REN21. (2019). Renewables 2019 Global status report. REN21.
http://www.ren21.net/status-of-renewables/global-
Fenercom. (2010). Manual Técnico de Energía Solar para status-report.
Procesos Industriales. https://www.fenercom.com/wp-
content/uploads/2010/11. Sharma, A.K., Sharma, C., Mullic, S. C., Kandpal, T. C. (2017).
Solar industrial process heating: a review. Renewable and
Frein, A., Calderoni, M., Mott, M. (2014). Solar thermal plant Sustainable Energy Reviews, 78, October, 124-137.
integration into an industrial process. Energy Procedia, 8,
152-1163. Solar Payback. (2017). Solar heat for industry.
https://www.international-climate-initiative.com/file
IEA-ETSAP e IRENA. (2015). Solar heat for industrial admin/Dokumente/2017/170530_Brochure_EN_Solar
processes technology Brief. https://www.irena.org/pub _Payback_digital.pd.pdf.
lications/2015/Jan/Solar-Heat-for-Industrial-Processes.
Suresh, N. S., Rao, B. S. (2017). Solar Energy for Process
IRENA. (2019). Transforming the energy system-and holding Heating: a Case Study of Select Indian Industries. Journal of
the line on the rise of global temperatures. https://www. Cleaner Production, 151, 439–451.
irena.org/publications/2019/Sep/Transforming-the-
energy-system. Weiss, W. and Spörk-Dür, M. (2019). Solar heat worldwide.
The Solar and Heating Programme, International Energy
Kalogirou, S. (2003). The potential of solar industrial Agency. https://www.iea-shc.org/Data/Sites/1/publica
process heat applications. Applied Energy, 76, 4, 337-361. tions/Solar-Heat-Worldwide-2019.pdf.
75
ESTUDIO EXPERIMENTAL DE LA PERMEABILIDAD
AL AIRE DE EDIFICIOS EN URUGUAY
EXPERIMENTAL STUDY OF THE AIR PERMEABILITY OF BUILDINGS IN URUGUAY
Juan Manuel Rodríguez Muñoz 1, Andrea Sevrini 2, Gunnar Gil 3, Italo Carlos Bove Vanzulli 4
76
ENERLAC • Volumen IV. Número 2. Diciembre, 2020. ISSN: 2602-8042 (impreso) / 2631-2522 (digital).
OLADE – AUGM
77
ESTUDIO EXPERIMENTAL DE LA PERMEABILIDAD AL AIRE DE EDIFICIOS EN URUGUAY.
Rodríguez Muñoz, Juan Manuel; Sevrini, Andrea; Gil, Gunnar; Bove Vanzulli, Italo Carlos
RESUMEN
ABSTRACT
78
ENERLAC • Volumen IV. Número 2. Diciembre, 2020. ISSN: 2602-8042 (impreso) / 2631-2522 (digital).
OLADE – AUGM
79
ESTUDIO EXPERIMENTAL DE LA PERMEABILIDAD AL AIRE DE EDIFICIOS EN URUGUAY.
Rodríguez Muñoz, Juan Manuel; Sevrini, Andrea; Gil, Gunnar; Bove Vanzulli, Italo Carlos
Es por todo ello que a través de este trabajo se Dada la no linealidad de la Ecuación (1) resulta
busca obtener, como objetivo primario, datos conveniente realizar el siguiente cambio de
empíricos de infiltraciones de aire en edificios en variable:
Uruguay: dichos datos son medidos de acuerdo
al método que se describe en la norma ISO Y = ln qenv (2)
9972:2015. En segundo lugar, se busca a través
de la comparación con los datos presentados X = ln Δp (3)
en publicaciones de otros países, responder a la
interrogante antes planteada de si es adecuado Entonces la Ecuación (1) se puede rescribir:
usar datos empíricos de otros países, o por el
contrario, es necesario invertir en continuar Y = ln Cenv + nX (4)
realizando medidas para lograr una mejor
caracterización y comprensión del parque Para determinar los coeficientes Cenv y n se deben
edilicio nacional. A la vez, a través de la tomar diferentes medidas de qenv y Δp para
comparación con estándares internacionales, luego realizar un ajuste lineal de la Ecuación (4).
se trata de inferir si es aceptable el nivel de per- El cálculo de Cenv y n y sus respectivas incer-
meabilidad al aire de los casos considerados. Por tidumbres se realiza de acuerdo a anexo C de la
último, dada la falta de datos en esta materia norma ISO 9972:2015.
a nivel regional, se busca difundir el método de
ensayo. Por esta razón, a continuación, se pone El coeficiente de caudal depende de la tem-
un énfasis especial en la descripción del método y peratura del aire durante el ensayo. Para
se explicita el alcance de éste, luego se presentan independizarnos de esto se define el coeficiente
los casos de estudio, se analizan sus resultados y, de fuga de aire CL, el cual se calcula corrigiendo
por último, se comparan los mismos con los el coeficiente de caudal Cenv como se indica en
obtenidos en otros trabajos a nivel internacional. las fórmulas (5) y (6), para despresurización y
presurización respectivamente.
MÉTODO DE ENSAYO
CL ≅ Cenv ( T0 ⁄ Te ) 1-n (5)
Fundamentos
CL ≅ Cenv ( T0 ⁄ Tint ) 1-n (6)
En la norma ISO 9972:2015 la permeabilidad al
aire del edificio se describe mediante la Ecuación: en donde Te y Tint son las temperatura exterior e
interior durante el ensayo (expresadas en Kelvin)
qenv = Cenv Δpn (1) y T0 la temperatura del aire en condiciones
estándar (20 °C).
en donde qenv (m3/h) es el caudal de aire a
través de la envolvente del edificio y Δp (Pa)
es la diferencia de presiones entre el aire
exterior e interior, Cenv (m3/h.Pan) se denomina En Uruguay en los
coeficiente de caudal y n coeficiente de presión. últimos años se ha
El coeficiente de presión está asociado al tipo
comenzado a medir la
de flujo y normalmente se encuentra entre los
valores n = 1/2 (turbulento) y n = 1 (laminar). La permeabilidad al aire de
determinación de la permeabilidad al aire de un ventanas en banco de
edificio se reduce a determinar los coeficientes de ensayo fijo.
caudal y presión.
80
ENERLAC • Volumen IV. Número 2. Diciembre, 2020. ISSN: 2602-8042 (impreso) / 2631-2522 (digital).
OLADE – AUGM
Definición de parámetros derivados a usar para su cálculo, pero existen varios modelos en
la literatura, uno muy sencillo consiste en dividir
Si bien la permeabilidad al aire de un edificio n50 por veinte (d’Ambrosio Alfano et al., 2012;
queda dada a través de los coeficientes de fuga Feijó-Muñoz et al., 2019); es decir:
y presión (CL y n), estos no resultan prácticos a
la hora establecer comparaciones entre distintos n0 = n50 ⁄ 20 (9)
edificios. Por este motivo, se introducen las
renovaciones hora a una diferencia de presión Desafortunadamente, en los artículos fuente de
exterior-interior de 50 Pascales valor deno- la información no se señala el autor.
minado n50. Esta cantidad se utiliza como
característica de la permeabilidad al aire de un La ASRHAE (2019) provee de otro método, el
edificio y es apta para establecer comparaciones. cual está basado en el modelo físico de Sherman
Las exigencias dispuestas en varios países son y Grimsrud (1980). Este modelo fue validado
establecidas en función de este parámetro. experimentalmente en la referencia (Sherman y
Además, esta magnitud es independiente de Modera, 1986) y da una expresión para el caudal
la implantación del edificio, lo que permite de infiltraciones q (m3/h) que ingresa a un edificio
comparar edificios situados en distintos lugares. por causa de la diferencia de temperatura ΔT (°C)
entre el aire exterior y el interior, y la velocidad
A partir del coeficiente de fuga CL y del coeficiente del viento v (m/s) a diez metros de altura:
de presión n, se calcula un caudal normalizado
q50 (m3/h):
AL
q = Cs | ∆T | + Cwʋ2 (10)
n 1000
q50 = CL 50 (7)
81
ESTUDIO EXPERIMENTAL DE LA PERMEABILIDAD AL AIRE DE EDIFICIOS EN URUGUAY.
Rodríguez Muñoz, Juan Manuel; Sevrini, Andrea; Gil, Gunnar; Bove Vanzulli, Italo Carlos
82
ENERLAC • Volumen IV. Número 2. Diciembre, 2020. ISSN: 2602-8042 (impreso) / 2631-2522 (digital).
OLADE – AUGM
En primer lugar, se debe preparar el edificio estar espaciadas 10 Pa máximo, el valor mínimo
según el método de ensayos escogido. Poste- inicial debe ser de 10 Pa y el valor máximo
riormente, se hace un relevamiento de las debe ser de al menos 50 Pa, pero se recomien-
condiciones climáticas: temperatura y velocidad da extender el rango hasta 100\Pa. La norma
del viento. Luego, se procede a determinar la recomienda realizar dos series de mediciones,
diferencia de presión a caudal cero. Para este una para presurización y otra para despresu-
paso se deben apagar y tapar los ventiladores. rización, sin embargo, se admite que se realice
Luego se debe registrar la diferencia de presión solo una serie de medidas.
exterior-interior durante un periodo de al menos
30 segundos. El promedio de los valores posi- Al final de estas series de medidas se debe
tivos Δp01+, el de los valores negativos Δp01- determinar nuevamente la diferencia de presión
y el de todos los valores Δp01 sobre este a caudal cero. Y nuevamente, el promedio de
periodo no deben superar los 5 Pa, en caso con- los valores positivos Δp02+, el de los valores
trario el ensayo se declara no válido. Para negativos Δp02- y el de todos los valores Δp02
alcanzar una diferencia de presión exterior- sobre este periodo no deben superar los 5 Pa,
interior a caudal cero aceptable (menor a 5 Pa) en caso contrario el ensayo se declara no válido.
se recomienda que la velocidad del viento Las medidas de presión en la Ecuación (1) son
durante el ensayo no alcance 3 en la escala corregidas con la diferencia de presión a flujo
Beaufort y que la diferencia de temperatura cero, es decir, a todas las medidas de Δp se les
exterior-interior multiplicada por la altura del resta el promedio entre Δp01 y Δp02.
edificio sea menor a 250 m.K.
Finalmente, los datos son procesados de
Una vez determinada la diferencia de presión acuerdo a la sección 2 y son determinados: el
exterior-interior a flujo cero, se destapan y se coeficiente de fuga CL, el coeficiente de presión
encienden los ventiladores y se procede a medir n, las renovaciones horas a una diferencia de
diferentes caudales de aire a distintas diferencias 50 Pa n50 y el promedio anual de la tasa de in-
de presión exterior-interior. Las medidas de filtraciones del edificio n0.
diferencia de presión exterior-interior deben
83
ESTUDIO EXPERIMENTAL DE LA PERMEABILIDAD AL AIRE DE EDIFICIOS EN URUGUAY.
Rodríguez Muñoz, Juan Manuel; Sevrini, Andrea; Gil, Gunnar; Bove Vanzulli, Italo Carlos
CASOS DE ESTUDIO
Se ensayó un total de siete edificios con diversas por el año de construcción, y en ellas se resumen
características. En las Tablas 2 y 3 se presentan las principales características de los edificios, las
los casos de estudio ordenados en forma creciente cuales se describen en los párrafos siguientes.
84
ENERLAC • Volumen IV. Número 2. Diciembre, 2020. ISSN: 2602-8042 (impreso) / 2631-2522 (digital).
OLADE – AUGM
De acuerdo a sus funciones, los casos de estudio Todos los edificios ensayados tienen en común
se dividen en dos grupos: servicios y residencial el sistema constructivo tradicional pesado en
(colectivo). Los edificios 1 al 5 pertenecen al mampostería de cerámica cocida, donde los casos
sector público, estando relacionados a funciones 6 y 7 presentan además estructura indepen-
educativas los casos 1, 2 y 3. El edificio del caso diente en hormigón armado. Las cubiertas son
4 tiene como función principal los servicios, pero de hormigón armado en los casos 2, 4, 5, 6 y 7.
el módulo ensayado pertenece al grupo resi- Los casos 1 y 3 cuentan con cubierta de cerámica
dencial, dado que es un edificio de uso mixto. armada abovedada con directriz catenaria. Las
El caso 5 es un edificio público destinado a la aberturas del grupo de edificios más antiguos
investigación científica. Dentro del grupo resi- son de diversa materialidad, encontrándose de
dencial, los casos 6 y 7 son unidades de edificios hierro, madera, aluminio y PVC. La forma
de apartamentos en altura. de apertura más frecuente es corrediza para
85
ESTUDIO EXPERIMENTAL DE LA PERMEABILIDAD AL AIRE DE EDIFICIOS EN URUGUAY.
Rodríguez Muñoz, Juan Manuel; Sevrini, Andrea; Gil, Gunnar; Bove Vanzulli, Italo Carlos
ventanas y batiente para puertas, aunque en Para el caso 4 se realizó una serie de medidas
el caso 3 las ventanas principales cuentan con de despresurización y otra de presurización.
un sector oscilante y el resto de la abertura Para el resto de los casos, las series de medidas
con vidrio fijo. En el caso 4 hay banderolas fueron siempre de despresurización.
proyectantes en los servicios sanitarios. En
la Tabla 3 se indica el material y la forma En todos los casos la velocidad del viento
de apertura de las aberturas de cada caso. durante el ensayo fue menor a 3 en la escala
El caso 2 cuenta con una chimenea como parte Beaufort. Además, en todos ellos la diferencia
del equipamiento de calefacción la cual no fue de temperatura exterior-interior multiplicada
sellada durante el ensayo. El caso 5 tiene cañe- por la altura del edificio fue menor a 250 K.m,
rías de salida de conductores eléctricos por a excepción del caso 6 donde no se midió
su cubierta, necesarias para el equipamiento temperatura. Sin embargo, en todos los casos
de medición ubicado sobre la misma. Estos se logró una diferencia de presión exterior-
elementos influyen en la hermeticidad de interior a caudal cero aceptable (< 5 Pa).
los edificios considerados. Los casos 1 y 3
tienen una particularidad importante: origi- RESULTADOS Y ANÁLISIS
nalmente eran construcciones similares, pero
posteriormente el caso 1 fue modificado Para comparar los resultados entre sí y a su vez
reemplazando las ventanas por otras de también con los resultados de otras publica-
diferente material y mejor calidad. Esto per- ciones, se determinó el valor de n50 para cada
mitió comparar dos edificios esencialmente caso. Los resultados se exponen en la Tabla 4
iguales, casos 1 y 3, con aberturas de distinta y se presentan de forma gráfica en la Figura 3 solo
calidad. los de despresurización.
Tabla 4. Resultados
1 - D – Despresurización, P – Presurización
En la Tabla 4 se muestra además el coeficiente n varía entre 0.541 y 0.830 y el de fuga de aire
de fuga de aire CL, el coeficiente de presión n y entre 20.3 y 1,562 m3/(h.Pan). El error relativo
el coeficiente de correlación r2. El coeficiente promedio del coeficiente de presión es del 2% y
de correlación fue mayor a 0.99 en todos los del de fuga de aire 8%. El promedio general de
casos excepto en el 7. El coeficiente de presión n50 (tomando en cuenta para el caso 4 solo el
86
ENERLAC • Volumen IV. Número 2. Diciembre, 2020. ISSN: 2602-8042 (impreso) / 2631-2522 (digital).
OLADE – AUGM
valor de despresurización) fue de 11.0 h-1 con lo sugerido por la norma de realizar únicamente
una gran variabilidad, que se traduce en una una serie de medidas (presurización o des-
alta desviación estándar, la cual es de 10.6 h-1. presurización).
En el caso 4, como fue mencionado ante- Resulta interesante notar que, como se dijo
riormente, se realizó una serie de medidas de antes, los casos 1 y 3 consisten en edificios
despresurización y otra de presurización. Como esencialmente iguales, pero con aberturas de
se puede ver en la Tabla 4 los resultados de diferente calidad. Se puede ver en la Tabla 4
ambos ensayos son muy parecidos, difiriendo que a través del cambio de aberturas en este
menos de un 6%, lo cual coincide con lo visto caso se logra una mejora de casi el 50% en el
en otras publicaciones (Targo, 2007) y valida valor de n50.
Figura 3. Valores de n50 para los casos de estudio y niveles de referencia internacional
35
Mediciones
n = . h
30
25
20
n (h)
15
10
0
1 2 3 4 5 6 7
Casos de estudio
87
ESTUDIO EXPERIMENTAL DE LA PERMEABILIDAD AL AIRE DE EDIFICIOS EN URUGUAY.
Rodríguez Muñoz, Juan Manuel; Sevrini, Andrea; Gil, Gunnar; Bove Vanzulli, Italo Carlos
88
ENERLAC • Volumen IV. Número 2. Diciembre, 2020. ISSN: 2602-8042 (impreso) / 2631-2522 (digital).
OLADE – AUGM
89
ESTUDIO EXPERIMENTAL DE LA PERMEABILIDAD AL AIRE DE EDIFICIOS EN URUGUAY.
Rodríguez Muñoz, Juan Manuel; Sevrini, Andrea; Gil, Gunnar; Bove Vanzulli, Italo Carlos
Al comparar los datos que aquí se presentan, por la Ecuación (1) ajusta muy bien los datos
se puede observar que los resultados de los experimentales y que la metodología ha sido
trabajos de los países mediterráneos son bas- correctamente implementada. A esto último,
tante similares entre sí. Si bien los tres trabajos se suma el hecho de que se utilizaron instru-
evalúan edificios en cada país de diferentes an- mentos calibrados y que se siguió un método
tigüedades, las desviaciones estándar de los estandarizado.
promedios presentados son relativamente bajas,
mostrando una cierta homogeneidad, lo cual El valor promedio de n50 de los casos estudiados
habla del buen estado de conservación y mante- fue de 11.0 h-1 con una desviación estándar
nimiento de los edificios evaluados. Por otro de 10.6 h-1 Al comparar los resultados con
lado, el trabajo chileno, si bien es sobre una la exigencia del Código Técnico de la Edifica-
amplia muestra de edificaciones relativamente ción de España (n50 < 6 h-1) se encontró que los
nuevas, muestra un promedio del orden de un edificios nuevos y/o reacondicionados satisfa-
50% mayor que el de las edificaciones europeas cen dicha exigencia en tanto que, los antiguos
acá presentadas. Respecto a las medidas realiza- sin reacondicionar, no lo hacen. Además, el valor
das en el presente trabajo, una primera mirada de n50 de los edificios del primer grupo casi
muestra una similitud entre el promedio de ellas cuadriplica el valor del primer grupo.
y el promedio del trabajo chileno. Sin embargo,
la alta desviación estándar de nuestros resulta- Partiendo del modelo de Sherman (1980) se
dos muestra que hay una alta diversidad de pudo calcular la tasa de infiltraciones promedio
resultados entre sí. Por un lado, están los anual n0, el valor el valor promedio general de n0
edificios antiguos sin reacondicionados, con un fue de 0.55 h-1, el de los edificios de construcción
promedio de n50 de 19.7 h-1 y desviación estándar nueva más el caso reacondicionado 0.22 h-1, y
de 11.7 h-1 y, por otro, los edificios nuevos y/o el de construcción antigua sin reacondicionar
reacondicionados, con un promedio de 4.4 h-1 0.98 h-1, el cual es más del cuádruple que el
y desviación estándar de 0.7 h-1. Por ende, los primero. La comparación de estos resultados
resultados acá medidos muestran una gran con el criterio utilizado en FADU deja en eviden-
diferencia tanto con los resultados de Chile para cia que este criterio debe ser revisado y
nuevas construcciones como con los resulta- actualizado, fundamentalmente para los edificios
dos europeos. Se observa entonces que utilizar de reciente construcción.
valores empíricos de infiltraciones de otros
países, a la hora de modelar construcciones lo- Finalmente, al comparar los valores de n50
cales, conlleva a un error. Por lo tanto, será obtenidos con los resultados de otras publica-
necesario continuar midiendo las infiltraciones ciones, se encuentran diferencias importantes
de edificaciones en Uruguay, ya sean nuevas, en los promedios y en la variabilidad de los re-
antiguas o recicladas, en cantidad suficiente como sultados (desviación estándar). A partir de esto,
para tener valores medios característicos que y de lo expuesto en el párrafo anterior, se
puedan ser usados de referencia. concluye que es importante que se sigan rea-
lizando ensayos en Uruguay que permitan tener
CONCLUSIONES valores más acordes al momento de modelar
térmicamente los edificios locales. A futuro, sería
Se determinó con éxito la permeabilidad al aire deseable ensayar una cantidad significativa de
de siete edificios en Uruguay de acuerdo al edificios, con diversas características, como para
estándar ISO 9972:2015. Se hallaron los inferir conclusiones más generales de la situación
coeficientes CL y n con un muy buen nivel de actual del parque edilicio nacional. Ello podría
incertidumbre (error relativo de 8%). El coe- servir como insumo para diseñar normativas
ficiente de correlación en 6 casos fue mayor adecuadas a nuestro país con el objetivo de limi-
a 0.99. Se deduce que el modelo teórico dado tar el consumo de energía de los edificios.
90
ENERLAC • Volumen IV. Número 2. Diciembre, 2020. ISSN: 2602-8042 (impreso) / 2631-2522 (digital).
OLADE – AUGM
REFERENCIAS
Alonso-Suárez, R., Bidegain, M., Abal, G. y Modernell, P. IEA. (2016). Final Energy Consumption in the World.
(2019). Año Meteorológico Típico para aplicaciones de Obtenido de https://www.iea.org/Sankey/#?c=World&s=
Energía Solar (AMTUes): series horarias típicas para 5 sitios Final%20consumption
del Uruguay. Montevideo, Uruguay: Ediciones Universitarias
(UCUR) de la Universidad de la República (Udelar). ISO. (2015). ISO 9972. Thermal performance of buildings
– determination of the air permeability of buildings – fan
ASHRAE. (2009). 2009 ASHRAE Handbook Fundamentals. pressurization method. Geneva, Switzerland.
Atlanta.
MIEM. (2018). Balance Energético Nacional de Uruguay.
Casañas, V., Pereira, L., Liston, M. y Cataldo, J. (2016). Obtenido de https://ben.miem.gub.uy/index.html
Caracterización y valoración energética de infiltraciones
por aberturas en el parque habitacional uruguayo. XXXIX Ministerio de Fomento. (2019). Código Técnico de
Reunión de Trabajo de la Asociación Argentina de Energías la Edificación de España. Obtenido de https://www.
Renovables y Medio Ambiente, 4, págs. 07.11-07.18. codigotecnico.org/images/stories/pdf/ahorroEnergia/
Argentina. DBHE.pdf
Centro de Investigación en Tecnologías de la Construcción Passive House. (s.f.). Passive House Institute. Obtenido de
de la Universidad del Bío-Bío. (2014). Manual de https://passivehouse.com/
hermeticidad al aire de edificaciones. Obtenido de
http://construccionsustentable.uc.cl/images/Documentos Ramos, N. M., Almeida, R. M., Curado, A., Pereira, P. F.,
/Manual_de_hermeticidad_al_aire_de_edificaciones.pdf Manuel, S., and Maia, J. (2015). Airtghtness and ventilation
in mild climate counrty rehabilitated social housing
d’Ambrosio Alfano, F., Dell’Isola, M., Ficco, G. and Tassini, buildings - What users want and what they get. Building and
F. (2012). Experimental analysis of air tightness in Environment (92), 97-110.
Mediterranean buildings using the fan pressuritation
method. Building and Environment, 53, 16-25. Sherman, M.H. and Modera, M.P. (1986). Comparation
of Measured and Predicted infiltration Using the LBL
Facultad de Arquitectura, Diseño y Urbanismo. (2002). Infiltration Model. Measured Air Leakage of Buildings,
AT02 Repartido. Cátedra de Acondicionamiento Térmico. ASTM STP 904, H.R. Trechsel and P.L Lagus, Eds., American
Montevideo, Uruguay. Obtenido de http://www.fadu. Society for Testing and Materials, Philadelpphia, 325-347.
edu.uy/acondicionamiento-termico/files/2012/02/
repartido-02.pdf Sherman, M.H., Grimsrud, D.T. (1980). Infiltration-
presurization correlation: Simplify physical modeling.
Feijó-Muñoz, J., González-Lezcano, R. A., Poza-Casado, I., ASHRAE Transactions, 86(2), 778
Padilla-Marcos, M. Á. and Meiss, A. (2019). Airtightness
of residential buildings in the Continental area of Spain. Targo, K. (2007). Air tightness and air leakages of new
Building and Environment, 148, 299-308. lightweigth single-family datached houses in Estonia.
Building and Environment, 42, 2369-2377.
Feijó-Muñoz, J., Pardal, C., Echarri, V., Fernandéz-Agüera, J.,
Assiego de Larriva, R., Montesdeoca Calderín, M., . . . Meiss,
A. (2019). Energy impact of the air infiltration in residential
buildings in the Mediterranean area of Spain and the Canary
inslans. Energy & Buildings, 226-238.
91
DESARROLLO DE NUEVOS CATALIZADORES BASADOS
EN NÍQUEL PARA LA PRODUCCIÓN DE GAS NATURAL
SINTÉTICO A PARTIR DE LA HIDROGENACIÓN DE CO2
DEVELOPMENT OF NEW NICKEL-BASED CATALYSTS FOR SYNTHETIC NATURAL GAS
PRODUCTION FROM CO2 HYDROGENATION
92
ENERLAC • Volumen IV. Número 2. Diciembre, 2020. ISSN: 2602-8042 (impreso) / 2631-2522 (digital).
OLADE – AUGM
RESUMEN ABSTRACT
La hidrogenación de CO2 puro fue estudiada, The hydrogenation of pure CO2 was studied
para la producción de CH4 (gas natural sintético), to produce CH4 (synthetic natural gas), on
sobre catalizadores NiLaMe (Me = Ce, Ti, Zr). NiLaMe catalysts (Me = Ce, Ti, Zr). The catalysts
Los catalizadores fueron sintetizados por el were synthesized by the modified polymer
método del precursor polimérico modificado precursor method and characterized by different
y caracterizados mediante diferentes técnicas physicochemical techniques. In all cases, the
fisicoquímicas. En todos los casos, los catali- reduced catalysts are composed of Ni (12% w/w)
zadores reducidos están compuestos por Ni and a mixed oxide of defined stoichiometry.
(12% p/p) y un óxido mixto de estequiometría The NiLaZr catalyst that exhibits higher dis-
definida. El catalizador NiLaZr que presenta persion of the active phase, and a higher
mayor dispersión de la fase activa, y una mayor concentration of the basic surface sites leads to
concentración de los sitios básicos superfi- a higher conversion of CO2 and greater selec-
ciales conduce a una mayor conversión de CO2 tivity to CH4. This catalyst exhibits good stability
y mayor selectividad a CH4. Este catalizador and selectivity for 100 hours at 400 °C and
presenta buena estabilidad y selectividad a then also evaluated in the direct hydrogenation
CH4 durante 100 horas de reacción a 400 °C. of a model biogas demonstrating to have good
Este mismo catalizador fue evaluado en la catalytic properties to efficiently transform CO2
hidrogenación directa de un biogás modelo into CH4 and thereby to increase its heating value
demostrando poseer buenas propiedades cata- in several applications.
líticas para transformar de manera eficiente
el CO2 en CH4 y, de esa forma, aumentar su
contenido energético para varias aplicaciones.
93
DESARROLLO DE NUEVOS CATALIZADORES BASADOS EN NÍQUEL PARA LA PRODUCCIÓN DE GAS NATURAL SINTÉTICO A PARTIR DE LA HIDROGENACIÓN DE CO2.
Veiga Bacci, Santiago German; Bussi, Juan
94
ENERLAC • Volumen IV. Número 2. Diciembre, 2020. ISSN: 2602-8042 (impreso) / 2631-2522 (digital).
OLADE – AUGM
Para la reacción de Sabatier puede utilizarse se está investigando esta reacción como una
CO2 recuperado de diferentes procesos indus- alternativa para la valorización de biogás
triales (Ashok et al., 2017), aunque también (Witte et al., 2019).
CAPTURA
ENERGÍA
CO
HIDROELÉCTRICA
RENOVABLE
ELECTRICIDAD
REACTOR
METANACIÓN
GAS NATURAL
H CATALIZADORES SINTÉTICO
ELECTRÓLISIS Ni - La - Me
DE AGUA
(Me - Ce, Zr, Ti)
TRANSPORTE CALOR
Este biocombustible disponible en Uruguay y en electricidad por combustión (Cipiti et al., 2016).
la región, se obtiene a través de la degradación Sin embargo, el uso directo del biogás tiene dos
anaerobia de materiales orgánicos y está com- inconvenientes principales debido al contenido
puesto principalmente por CH4 y CO2 con algunos relativamente alto de CO2 (40-60%) que
componentes minoritarios (H2S, NH3, siloxanos, disminuye su poder calorífico y la generación de
etc.). Su uso directo incluye la generación de mayores emisiones de CO2 a la atmósfera.
95
DESARROLLO DE NUEVOS CATALIZADORES BASADOS EN NÍQUEL PARA LA PRODUCCIÓN DE GAS NATURAL SINTÉTICO A PARTIR DE LA HIDROGENACIÓN DE CO2.
Veiga Bacci, Santiago German; Bussi, Juan
Por lo anterior, la hidrogenación de esta fracción En un trabajo anterior se ha sintetizado una serie
de CO2 presente en el biogás es una opción de catalizadores conteniendo 12% de níquel que
atractiva para su valorización. Para ello se podría se encuentra disperso sobre un óxido de lantano-
utilizar el CO2 que se obtuvo en etapas previas cerio y lantano-zirconio que fueron evaluados en
de separación o realizar la hidrogenación directa el reformado con vapor de glicerina cruda para la
del biogás. En este último caso se eliminaría la producción de hidrógeno (Veiga et al., 2020).
necesidad de retirar el CO2 que normalmente
se realiza mediante la absorción con aminas El objetivo del presente trabajo es sumar a los
consumiendo una gran cantidad de energía. catalizadores antes mencionados uno en el cual
el soporte esté constituido por lantano-titanio y
Sin embargo, existen importantes limitaciones evaluar estos catalizadores en la hidrogenación
cinéticas que dificultan la hidrogenación directa de CO2 puro y en la hidrogenación de un biogás
de CO2 a CH4, por lo que es necesario utilizar modelo para la producción de CH4 (GNS). En forma
catalizadores que presenten elevada actividad y adicional, con fines comparativos, se prepararon
selectividad a CH4 (Pastor-Pérez et al., 2018). y evaluaron los catalizadores Ni/Al2O3 y Rh/Al2O3
que están ampliamente documentados en la
En los últimos años, se han publicado una gran literatura existente. (Swalus et al., 2012).
cantidad de trabajos (Gödde et al., 2021; Zhou
et al., 2016; Zhang et al., 2020) con el objetivo MATERIALES Y MÉTODOS
de desarrollar un catalizador que presente
una elevada estabilidad térmica con una alta Preparación de catalizadores
conversión de CO2 y elevada selectividad a CH4.
Los catalizadores de níquel son los más La síntesis de los catalizadores se realizó
investigados para aplicaciones en procesos mediante el método del precursor polimérico
de hidrogenación a escala industrial debido a modificado que es una variación del método de
su elevada actividad y relativo bajo costo. Sin Pechini (Gaura et al., 2011). En este se polimeriza
embargo, el catalizador convencional Ni/Al2O3 una solución de etilenglicol (E), ácido cítrico
sufre desactivación incluso a temperaturas (AC) y los precursores metálicos. Las sales
de reacción bajas debido principalmente a la precursoras utilizadas fueron Ni(NO3)2.5H2O,
sinterización de partículas de Ni y la deposición La(NO3)3.6H2O, Ce(NO3)3.6H2O, Zr(OiPro)4 70%
de carbón (Falbo et al., 2017). Además de la en isopropanol y Ti(OCH2CH2CH3)4. La relación
fase activa, la elección de un soporte adecuado molar E:AC utilizada fue 4:1 y la relación molar
es fundamental para obtener una importante de los metales para los cinco catalizadores
actividad catalítica y estabilidad. Los soportes preparados fue: La:Ce=1:1, La:Ti=1:1, La:Zr=1:1.
más reportados son Al2O3, TiO2, SiO2, La2O3, En todos los casos el contenido de níquel fue
CeO2 y ZrO2. Para resolver los problemas 12% en peso. La mezcla resultante se agita por 1
asociados a los fenómenos de desactivación se hora a 70 °C. Luego se calienta hasta 130 °C para
vienen estudiando nuevos catalizadores prepa- favorecer la polimerización y remover el exceso
rados a partir de óxidos mixtos tipo perovskita de solvente. La resina resultante se pre-calcina a
o pirocloro que han demostrado tener buena 300 °C por 2 horas para posteriormente calcinar
actividad y elevada estabilidad térmica (Sekine el sólido resultante a 850 °C (temperatura
et al., 2013). Algunos de estos óxidos son necesaria para obtener el sistema bifásico NiO-
utilizados por esta última propiedad como La2Me2O7 por 4 horas. Para su identificación
recubrimiento de barrera térmica en turbinas de en el resto del trabajo los catalizadores serán
gas (Ramachandran et al., 2013). citados como NiLaCe, NiLaTi, y NiLaZr.
96
ENERLAC • Volumen IV. Número 2. Diciembre, 2020. ISSN: 2602-8042 (impreso) / 2631-2522 (digital).
OLADE – AUGM
Ensayos catalíticos
97
DESARROLLO DE NUEVOS CATALIZADORES BASADOS EN NÍQUEL PARA LA PRODUCCIÓN DE GAS NATURAL SINTÉTICO A PARTIR DE LA HIDROGENACIÓN DE CO2.
Veiga Bacci, Santiago German; Bussi, Juan
Los datos termodinámicos fueron calculados conversión de CO2, la selectividad hacia CH4 y CO
utilizando el software Fact-Web. Se trabajó en se determinaron por las fórmulas siguientes y el
el rango 250-500 °C. Para cada temperatura, T50 se define como la temperatura necesaria para
el análisis de la fase gaseosa se realizó luego convertir el 50% de CO2.
de alcanzar el estado estacionario (1 hora). La
moles CH₄ o CO
Selectividad CH₄ o CO (%) = * 100 (4)
moles CH₄ + CO
RESULTADOS Y DISCUSIÓN
Tabla 1. Parámetros texturales y estructurales de catalizadores Ni-La-Me (Me = Ce, Zr, Ti)
Varios trabajos han demostrado que los sitios uno de los reactivos es el CO2, que es ligeramente
básicos son muy importantes para la metanación ácido, la presencia de sitios básicos es útil en
de CO2 (Vita et al., 2018), por lo tanto, se realiza- términos de adsorción y activación de este
ron los perfiles de desorción térmica programada gas. Los sitios básicos más débiles causan la
de CO2 para los catalizadores reducidos para desorción de CO2 adsorbido a temperaturas más
cuantificar la basicidad superficial. Dado que bajas, mientras que los sitios básicos más fuertes
98
ENERLAC • Volumen IV. Número 2. Diciembre, 2020. ISSN: 2602-8042 (impreso) / 2631-2522 (digital).
OLADE – AUGM
NiLaZr
NiLaCe
TEMPERATURA (C)
La conversión de CO2 entre 200 – 500 °C sobre utilizó una relación H2 /CO2 = 4 y una velocidad
los catalizadores NiLaCe, NiLaTi y NiLaZr se espacial (GHSV) de 21000 mL h-1 gcat-1.
presenta en la Figura 3. Para estos ensayos se
99
DESARROLLO DE NUEVOS CATALIZADORES BASADOS EN NÍQUEL PARA LA PRODUCCIÓN DE GAS NATURAL SINTÉTICO A PARTIR DE LA HIDROGENACIÓN DE CO2.
Veiga Bacci, Santiago German; Bussi, Juan
Figura 3. Conversión de CO2 para los catalizadores Ni-La-Me (Me = Ce, Ti, Zr)
100 100
80 Equ 80
ilib rio
CONVERSIÓN CO₂ (%)
60 60
55 C 99 C
40 40
NiLaTi
20 20
NiLaCe
NiLaZr
0 0
250 300 350 400 450 500
TEMPERATURA ( C)
La línea continua representa la conversión de 250 a 400 °C, la conversión de CO2 se incre-
CO2 predicha por el equilibrio termodinámico. mentó gradualmente, mientras que, a partir
En las condiciones de trabajo la conversión de de los 400 °C, la conversión disminuye con
CO2 está lejos de los valores de equilibrio por la temperatura (debido al hecho que la meta-
debajo de 400 °C. Aunque la hidrogenación de nación de CO2 es exotérmica) con valores
CO2 es termodinámicamente favorable a bajas cercanos al equilibrio termodinámico. Para
temperaturas, lograr una elevada conversión a los otros dos catalizadores se observa un
tales temperaturas es un desafío. Esto se debe aumento en la conversión de CO2 con la
principalmente a la dificultad asociada con la temperatura, pero ambos presentan valores
activación de CO2 y cinética de la reacción a bajas de T50 inferiores que el catalizador NiLaZr:
temperaturas. Para el catalizador NiLaZr, desde (55 °C vs. NiLaCe) y (154 °C vs. NiLaTi).
100
ENERLAC • Volumen IV. Número 2. Diciembre, 2020. ISSN: 2602-8042 (impreso) / 2631-2522 (digital).
OLADE – AUGM
Figura 4. Selectividad de catalizadores Ni-La-Me (Me = Ce, Ti, Zr) en la hidrogenación de CO2 puro
100 100
SELECTIVIDAD CH₄ (%)
90 90
80 80
Equilibrio
70 70
NiLaTi
60 60
NiLaCe
50 50
NiLaZr
40 40
250 300 350 400 450 500
60 60
SELECTIVIDAD CO (%)
50 50
40 40
30 30
20 20
10 10
0 0
250 300 350 400 450 500
TEMPERATURA (C)
La selectividad hacia CH4 y CO para los tres cata- presenta diferencias significativas con los dos
lizadores ensayados se muestra en la Figura 4. anteriores. Este catalizador convierte CO2 vía
Como se comentó anteriormente la metanación RDAR pero no es eficaz en la metanación. Otros
de CO2 es exotérmica, sin embargo, la RDAR es posibles productos como etano, etileno y propano
endotérmica. Por lo tanto, a bajas temperaturas no se detectaron durante las experiencias.
la selectividad a CH4 es cercana al 100%
mientras que al ir aumentando la temperatura La mejor performance catalítica del catalizador
la selectividad de CH4 comienza a disminuir NiLaZr está atribuida a una mejor dispersión
mientras que la de CO aumenta. Los valores para de Ni observada en los ensayos de quimisorción
NiLaZr y NiLaCe están cercanos a lo que predice de H2 y a una mayor densidad de sitios básicos
el equilibrio termodinámico mientras que NiLaTi observada en los ensayos de TPD-CO2.
101
DESARROLLO DE NUEVOS CATALIZADORES BASADOS EN NÍQUEL PARA LA PRODUCCIÓN DE GAS NATURAL SINTÉTICO A PARTIR DE LA HIDROGENACIÓN DE CO2.
Veiga Bacci, Santiago German; Bussi, Juan
100 100
80 80
CONVERSIÓN Y SELECTIVIDAD (%)
60 60
Conversión CO₂
40 40
Selectividad CH₄
Selectividad CO
20 20
0 0
0 10 20 30 40 50 60 70 80 90 100
TIEMPO (h)
La actividad del catalizador NiLaZr fue luego la Figura 6. En ambos difractogramas pueden
evaluada en una experiencia durante 100 horas identificarse los picos asociados a la fase piro-
a 400 °C y los datos de conversión de CO2 y cloro La2Zr2O7 en 2θ=28.7, 33.4, 47.8 y 56.6° y los
selectividad a CH4 y CO se muestran en la Figura picos en 2θ=44.5° y 51.8° son asociados al níquel
4. La conversión se mantuvo estable alrededor de metálico. Utilizando la ecuación de Scherrer se
72% y se observa un leve descenso en la calculó el tamaño de cristalita de níquel y del
selectividad a CH4 que pasa de 97.1% a 95.4% pirocloro antes y después del ensayo de 100
luego de las 100 horas de ensayo demostrando horas. No existen diferencias significativas en
una desactivación prácticamente nula en el tiem- ambos parámetros indicando que no existe sin-
po de reacción estudiado. terización apreciable de las partículas de níquel
y confirmando la buena estabilidad térmica y
Luego de este ensayo se procedió a estudiar química que presenta el compuesto La2Zr2O7
posibles modificaciones estructurales del cata- utilizado como soporte en este catalizador su-
lizador mediante difracción de rayos X. Con fines mado a la ausencia de otros picos como La2O3
comparativos también se presenta el catalizador y ZrO2 que podrían estar asociados a modifica-
sin usar. Ambos difractogramas se presentan en ciones estructurales del óxido mixto.
102
ENERLAC • Volumen IV. Número 2. Diciembre, 2020. ISSN: 2602-8042 (impreso) / 2631-2522 (digital).
OLADE – AUGM
Figura 6. Difractogramas del catalizador NiLaZr sin usar y usado durante 100 h a 400 °C
NiLaZr usado
DLa2Zr2O7 = 15,3 nm
DNi = 18,6 nm
INTENSIDAD (u.a)
DLa2Zr2O7 = 15,1 nm
DNi = 18,2 nm
20 30 40 50 60
2 THETA ()
En la Figura 7 se muestran los resultados de dad para este proceso posiblemente asociado a
la hidrogenación directa de un biogás modelo menor basicidad superficial.
(CH4 = 70% , CO2 = 30%) con una relación
CH4/CO2/H2 = 2.33/1/4 y una velocidad espacial La selectividad para los catalizadores ensayados
de 21000 mL h-1 gcat-1. En este ensayo el catali- se presenta en la Figura 8. No se observan dife-
zador NiLaZr fue comparado con el catalizador rencias significativas en los valores de selectivi-
Ni(12%)/Al2O3 que está ampliamente reportado dad alcanzando valores cercanos al equilibrio
en la literatura (Stangeland et al., 2018) y otro termodinámico para temperaturas menores a
catalizador basado en un metal noble como el 400 °C. Sin embargo, el gas obtenido no alcanza la
rodio soportado también sobre Al2O3 (Rh(1%) concentración necesaria para poder ser inyecta-
Al2O3). A temperaturas bajas este último cata- da directamente en la red de gas natural
lizador presenta la mayor conversión de CO2 (CH4 > 96%, H2 < 2% y CO2 < 4%) por lo que sería ne-
pero a partir de 400 °C los valores son similares cesario adicionar otra etapa de acondicionamiento
a los obtenidos por el catalizador NiLaZr. Los como etapas de separación de H2 y CO2 o utilizar
valores de conversión de CO2 para el catalizador un segundo reactor de metanación para cumplir
Ni(12%)/Al2O3 son inferiores a los que presenta con las especificaciones. Esta opción ya fue con-
el catalizador NiLaZr indicando una menor activi- templada por otros autores (Witte et al., 2018).
103
DESARROLLO DE NUEVOS CATALIZADORES BASADOS EN NÍQUEL PARA LA PRODUCCIÓN DE GAS NATURAL SINTÉTICO A PARTIR DE LA HIDROGENACIÓN DE CO2.
Veiga Bacci, Santiago German; Bussi, Juan
Equilibr
io
80
80
CONVERSIÓN CO₂ (%)
60
60
Equilibrio 40
40 NiLaZr
Ni/Al₂O₃
20
Rh/Al₂O₃
20
0
250 300 350 400 450 500
TEMPERATURA (C)
104
ENERLAC • Volumen IV. Número 2. Diciembre, 2020. ISSN: 2602-8042 (impreso) / 2631-2522 (digital).
OLADE – AUGM
100 100
98 98
Equilibrio
96 96
NiLaZr
94 Ni/Al₂O₃ 94
92 Rh/Al₂O₃ 92
90 90
250 300 350 400 450 500
10 10
SELECTIVIDAD CO (%)
8 8
6 6
4 4
2 2
0 0
250 300 350 400 450 500
TEMPERATURA (C)
105
DESARROLLO DE NUEVOS CATALIZADORES BASADOS EN NÍQUEL PARA LA PRODUCCIÓN DE GAS NATURAL SINTÉTICO A PARTIR DE LA HIDROGENACIÓN DE CO2.
Veiga Bacci, Santiago German; Bussi, Juan
AGRADECIMIENTOS
REFERENCIAS
Ashok, J., Ang, M. L. and Kawi, S. (2015). Enhanced activity of Hussain, I., Jalil, A.A., Hassan, N.S., Hambali, H.U. and Jusoh,
CO2 methanation over Ni/CeO2-ZrO2 catalysts: Influence of N.W.C. (2020). Fabrication and characterization of highly
preparation methods. Catalysis Today, 281, 2–9. active fibrous silica-mordenite (FS@SiO2-MOR) cockscomb
shaped catalyst for enhanced CO2 methanation. Chemical
Cipiti, F., Barbera, O., Briguglio, N., Giacoppo, G., Italiano, C. and Engineering Science, 228, 115978.
Vita, A. (2016). Design of a biogas steam reforming reactor:
A modelling and experimental approach. International Iaquaniello, G., Setini, S., Salladini, A. and Falco, M. De. (2018).
Journal of Hydrogen Energy, 41(27), 11577–11583. CO2 valorization through direct methanation of flue gas
and renewable hydrogen : A technical and economic
Falbo, L., Martinelli, M., Visconti, C. G., Lietti, L., Bassano, assessment. International Journal of Hydrogen Energy, 43,
C. and Deiana, P. (2017). Kinetics of CO2 methanation on a 17069-17081.
Ru-based catalyst at process conditions relevant for Power-
to-Gas applications. Applied Catalysis B: Environmental, Janke, C., Duyar, M.S., Hoskins, M. and Farrauto, R.J. (2014).
225(17), 354–363. Catalytic and adsorption studies for the hydrogenation of
CO2 to methane. Applied Catalysis. B Environmental. 152–
Gaura, S., Haynes, D. and Spivey, J. (2011). Rh, Ni, and 153.
Ca substituted pyrochlore catalysts for dry reforming of
methane. Applied Catalysis A: General, 403, 142-151. Milani, D., Khalilpour, R., Zahedi, G. and Abbas, A. (2015).
A model-based analysis of CO2 utilization in methanol
Gödde, J., Merko, M., Xia, W. and Muhler, M. (2021). Nickel synthesis plant. Journal of CO2 Utilization, 10, 12–22.
nanoparticles supported on nitrogen–doped carbon
nanotubes are a highly active, selective and stable CO2 Odumugbo, C.A. (2010). Natural gas utilisation in Nigeria:
methanation catalyst. Journal of Energy Chemistry, 54, Challenges and opportunities. Journal of Natural Gas Science
323–331. and Engineering, 2, 310-316.
Haynes, D. J., Shekhawat, D., Berry, D. A., Zondlo, J., Roy, A. Pastor-Pérez, L., Saché, E. L., Jones, C., Gu, S., Arellano-García,
and Spivey, J. J. (2017). Characterization of calcination H. and Reina, T. R. (2018). Synthetic natural gas production
temperature on a Ni-substituted lanthanum-strontium- from CO2 over Ni-x/CeO2 -ZrO2 (x = Fe, Co) catalysts:
zirconate pyrochlore. Ceramics International, 43(18), Influence of promoters and space velocity. Catalysis Today,
16744–16752. 317, 108-113.
106
ENERLAC • Volumen IV. Número 2. Diciembre, 2020. ISSN: 2602-8042 (impreso) / 2631-2522 (digital).
OLADE – AUGM
Ramachandran, C.S., Balasubramanian, V. and Witte, J., Settino, J., Biollaz, S., Tilman, J. and Schildhauer,
Ananthapadmanabhan, P.V. (2013). Thermal cycling T. (2018). Direct catalytic methanation of biogas – Part I:
behaviour of plasma sprayed lanthanum zirconate based New insights into biomethane production using rate-based
coatings under concurrent infiltration by a molten glass modelling and detailed process análisis. Energy Conversion
concoction. Ceramics International, 39(2), 1413-1431. and Management, 171, 750–768.
Sekine, Y., Mukai, D., Murai, Y., Tochiya, S., Izutsu, Y., Sekiguchi, Xu, L., Wang F., Chen, M., Nie, D., Lian, X., Lu, Z., Chen, H.,
K. & Sugiura, Y. (2013). Steam reforming of toluene over Zhang, K., Ge, P. (2017). CO2 methanation over rare earth
perovskite-supported Ni catalysts. Applied Catalysis A: doped Ni based mesoporous catalysts with intensified
General, 451, 160–167. low-temperature activity. International Journal of
Hydrogen Energy, 42, 15523-15539.
Spivey, J.J. and Hutchings, G. (2014). Catalytic aromatization
of methane. Chemical Society Reviews, 43, 792–803. Zhang, F., Lu, B.and Sun, P. (2020), Highly stable Ni-based
catalysts derived from LDHs supported on zeolite for CO2
Stangeland, K., Kalai, D. Y., Li, H. and Yu, Z. (2018). Active and methanation. International Journal of Hydrogen Energy,
stable Ni based catalysts and processes for biogas upgrading: 45, 16183-16192.
The effect of temperature and initial methane concentration
on CO2 methanation. Applied Energy, 227, 206-212. Zhou, R., Rui, N., Fan, Z. and Liu, C. (2016). Effect of
the structure of Ni/TiO2 catalyst on CO2 methanation.
Swalus, C., Jacquemina, M., Poleunis, C., Bertrand, P. and Ruiz, International Journal of Hydrogen Energy, 41(47), 22017–
P. (2012). CO2 methanation on Rh/g-Al2O3 catalyst at low 22025.
temperature: “In situ” supply of hydrogen by Ni/activated
carbon catalyst. Applied Catalysis B: Environmental, 125,
41–50.
Tada, S., Ikeda, S., Shimoda, N., Honma, T., Takahashi, M.,
Nariyuki, A. and Satokawa, S. (2017). Sponge Ni catalyst
with high activity in CO2 methanation. International Journal
of Hydrogen Energy, 42(51), 30126–30134.
Veiga, S., Romero, M., Faccio, R., Segobia, D., Duarte, H.,
Apesteguía, C., and Bussi, J. (2020). Hydrogen-rich gas
production by steam and oxidative steam reforming of
crude glycerol over Ni-La-Me mixed oxide catalysts (Me= Ce
and/or Zr). Catalysis Today, 344, 190–198.
107
RESIDUOS DE APARATOS ELÉCTRICOS Y
ELECTRÓNICOS (RAEE): IMPACTO SOCIAL,
AMBIENTAL, GESTIÓN Y METODOLOGÍAS
SOBRE SU MANEJO
WASTE OF ELECTRICAL AND ELECTRONIC DEVICES (WEEE):
SOCIAL, ENVIRONMENTAL IMPACT, MANAGEMENT AND
METHODOLOGIES ON ITS HANDLING
Dulce Kristal Becerra Paniagua 1, Araceli Hernández Granados 2, Evelyn Betsabe Díaz Cruz 3,
Karla Graciela Cedano Villavicencio 4, Horacio Martínez Valencia 5
108
ENERLAC • Volumen IV. Número 2. Diciembre, 2020. ISSN: 2602-8042 (impreso) / 2631-2522 (digital).
OLADE – AUGM
RESUMEN
109
RESIDUOS DE APARATOS ELÉCTRICOS Y ELECTRÓNICOS (RAEE): IMPACTO SOCIAL, AMBIENTAL, GESTIÓN Y METODOLOGÍAS SOBRE SU MANEJO.
Becerra Paniagua, Dulce Kristal; Hernández Granados, Araceli; Díaz Cruz, Evelyn Betsabe; Cedano Villavicencio, Karla Graciela; Martínez Valencia, Horacio
ABSTRACT INTRODUCCIÓN
As the technology market continues to evolve, Se definen como RAEE a los residuos o desechos
the flow waste from the electrical and electronic de aparatos eléctricos y electrónicos, o WEEE
devices (WEEE) generated may constitute an por sus siglas en inglés (Waste Electrical and
environmental and social threat, due in large Electronic Equipment) a cualquier aparato
part to the inadequate management of these eléctrico o electrónico cuyo periodo de vida
wastes, as well the insufficient collection services, útil ha concluido. Los dispositivos eléctricos o
and the existing disinformation on this topic. electrónicos (DE) son por definición aquellos
In this article its present a systematic review of que contienen circuitos y/o componentes eléc-
the literature published about WEEE, including tricos o que para su funcionamiento necesitan
information on their generation, destination and corriente eléctrica, campos electromagnéticos
environmentally friendly methodologies that o baterías (Blaisdell-Vidal, 2019). Los RAEE son
have been proposed for its management- generados cuando el usuario decide dejar de
recycling, as well as on environmental and utilizar sus DE como resultado de la creciente
health impacts derived from WEEE and its demanda de nuevos productos, falla física, las
handling. This review aims to sensitize the constantes actualizaciones o la llamada obso-
population about the danger of inadequate lescencia programada (problemas de compatibi-
treatment of WEEE and in turn to consider lidad) (Ongondo et al., 2011). Por ejemplo, en
the implementation of management strategies 1992 la vida útil (EoL, End of Life, en inglés) de
and actions to be taken for these, through las computadoras personales (PC, personal
environmentally friendly recycling processes. The computer, en inglés) era en promedio de 4.5 años,
first mentions of WEEE in scientific production esta decreció a 3 años en 1999 y pasó a ser de
date from 1996, evidencing that the topic has sólo 2 años en el 2005 (Culver, 2005; Widmer
continued to grow steadily until 2019, when 735 et al., 2005). Entre los RAEE más generados se
documents were published. The categories with encuentran teléfonos celulares, computadoras,
the greatest number of documents where the impresoras, televisores, tabletas, paneles fotovol-
greatest knowledge and high-impact scientific taicos, refrigeradores, entre otros más (Cucchiella
contributions have been generated include et al., 2015; Wong et al., 2007). Esto se debe en
those focused on the environment or energy sus- gran parte al deseo de los consumidores por
tainability, as well as the areas of metallurgy, adquirir los modelos más recientes en materia
chemistry, electricity and fossil fuels. In addition, tecnológica.
the mapping of electronic waste produced in
millions of tons that is produced worldwide is Datos recientes informan que anualmente en el
presented and the regions that send WEEE to mundo se generan de 20 a 50 millones de tonela-
other countries, mostly illegally, are shown. das aproximadamente de estos residuos (Tuncuk
On the other hand, the recycling of WEEE is a et al., 2012). Tan solo en el 2016 se generaron
complicated process, but with a high potential 44.7 millones de toneladas métricas (Mt) de
to generate a green circular economy. desechos, que significa 6.1 kilogramos por habi-
tante en el planeta (kg/hab), lo cual equivale
al peso de 9 pirámides de Giza o 4,500 torres
Keywords: WEEE, ODS, Electronic Waste, Eiffel; y se estima que esta cifra pueda incre-
Recycling, Electrochemistry, Metal Recovery. mentar a 52.2 Mt para el 2021 (Baldé et al., 2017);
esta estimación se puede observar en la Figura 1.
110
ENERLAC • Volumen IV. Número 2. Diciembre, 2020. ISSN: 2602-8042 (impreso) / 2631-2522 (digital).
OLADE – AUGM
Este incremento de RAEE no sólo origina la aluminio, bario, berilio, cadmio, cobalto, plomo,
generación de desechos sólidos, si no que otro níquel, entre otros más (Cesaro et al., 2018;
problema que se produce es el inadecuado ma- Yang et al., 2018). Debido a la magnitud de las
nejo de estos, ocasionando problemas como la problemáticas anteriormente mencionadas,
contaminación del suelo y agua; asimismo si los en este artículo se presenta una revisión
residuos se incineran, también se contamina el sistemática de la literatura publicada sobre los
aire 1, impactando de manera negativa la salud RAEE. Además, se proporciona información
humana y al ambiente (Dimitrakakis et al., 2009; sobre su generación, destino, impacto a la salud
Stenvall et al., 2013). Cuando los contaminantes y a la calidad del ambiente, con el objetivo de
penetran en el suelo, pueden filtrarse al sub- sensibilizar e informar sobre el peligro que exis-
suelo, de manera que alteran la calidad de las te en el tratamiento inadecuado de los RAEE.
aguas subterráneas. Las personas que habitan Adicionalmente, se muestra la importancia de
cerca de los sitios de reciclaje al tomar agua, formular e implementar estrategias de gestión
usarla en la comida o respirar aire contaminado de residuos mediante procesos de reciclaje am-
quedan expuestos a compuestos tóxicos como bientalmente amigables.
Figura 1. Generación mundial de residuos electrónicos por habitante y total mundial, periodo
2014 a 2016
50
40
44.7
(Mt)
30
20
10
0
2014 2015 2016 2017 2018 2019 2020 2021
0
1
por habitante
2
(kg/hab)
3
4
5
6.1
6
7
Nota: Los valores de 2017 a 2021 son estimaciones del estudio de Baldé et al. (2017)
Fuente: Elaboración propia con datos basados del estudio de Baldé et al. (2017)
111
RESIDUOS DE APARATOS ELÉCTRICOS Y ELECTRÓNICOS (RAEE): IMPACTO SOCIAL, AMBIENTAL, GESTIÓN Y METODOLOGÍAS SOBRE SU MANEJO.
Becerra Paniagua, Dulce Kristal; Hernández Granados, Araceli; Díaz Cruz, Evelyn Betsabe; Cedano Villavicencio, Karla Graciela; Martínez Valencia, Horacio
112
ENERLAC • Volumen IV. Número 2. Diciembre, 2020. ISSN: 2602-8042 (impreso) / 2631-2522 (digital).
OLADE – AUGM
Figura 2. Producción científica por año durante el periodo de 1996-2020, para los criterios de
búsqueda establecidos
800
735 716 Número de publicaciones
700
628
600 586
500 493
400
346 347 330
300 254 236
200 181 179
122 132
100 88 87 92
46 31 22 28
10 10 13 6
0
20 019 018 017 016 015 014 013 012 011 010 009 008 007 006 005 004 003 002 001 000 999 998 997 996
20 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 1 1 1 1
6 7 2
1,380,004,385
1,439,323,776
CHINA INDIA
CHINA EUA INDIA
Fuente: Elaboración propia con información de Web of Science (WOS, 2020) y cifras de la población mundial por país
(Worldometer, 2020).
113
RESIDUOS DE APARATOS ELÉCTRICOS Y ELECTRÓNICOS (RAEE): IMPACTO SOCIAL, AMBIENTAL, GESTIÓN Y METODOLOGÍAS SOBRE SU MANEJO.
Becerra Paniagua, Dulce Kristal; Hernández Granados, Araceli; Díaz Cruz, Evelyn Betsabe; Cedano Villavicencio, Karla Graciela; Martínez Valencia, Horacio
En la Figura 4, se muestran las categorías con destacan las enfocadas al medio ambiente o a
mayor número de documentos donde se ha la sustentabilidad energética, así también las
generado mayor conocimiento y aportaciones áreas de metalurgia, química, eléctrica y de
científicas de alto impacto. Dentro de estas, combustibles fósiles.
ORIGEN DE LOS RAEE Y MATERIALES London lo acuñó en uno de sus ensayos en 1932
CONTENIDOS (London, 1932) y se popularizó en 1954 por
Brooks Stevens, el cual lo llevó a la definición
Entre los diferentes tipos de dispositivos de obsolescencia planeada como: “Inculcar en el
eléctricos y electrónicos que han llamado más comprador el deseo de poseer algo un poco más
la atención por su amplia popularidad de uso nuevo, un poco mejor, un poco antes de lo nece-
y desechos que generan, se encuentran los sario”. De acuerdo con los datos de GSMA2, hasta
celulares, computadoras, pantallas y monitores el día de hoy existen 5.8 mil millones de perso-
de CRT, LCD y LED, tabletas, placas de circuitos nas que poseen un dispositivo móvil en el mundo,
impresos (PCB, Printed Circuit Board, en lo cual significa que al 67.95% de la pobla-
inglés) discos duros de estado sólido y paneles ción mundial posee un dispositivo de este tipo
fotovoltaicos (Cucchiella et al., 2015). Por (Turner, 2020). En el 2004, se estimó que se
ejemplo, los celulares son dispositivos que vendieron más de 180 millones de computadoras
poseen diferentes tipos de baterías como personales (PC), aproximadamente 100 millones
níquel-cadmio (NiCd), níquel-metal hidruro y de PC’s fueron desechadas y en el mismo año
de iones de litio (Li-ion), estas baterías poseen alrededor de 315 millones de PC’s quedaron
elementos altamente contaminantes. En 2006 obsoletas. Dejando a su paso una generación de
se estimó que se desecharon mil millones de aproximadamente 544,310 toneladas de plomo,
celulares (Rahmani et al., 2014). Sin embargo, es 907 toneladas de cadmio y 1,814,369 toneladas
difícil calcular con certeza los RAEE derivados de de plástico (Babu et al., 2007; Daven & Klein,
esta industria ya que estos pueden permanecer 2008). Otros materiales contaminantes presen-
en casa de los usuarios por años, incluso tes en su gran mayoría en las computadoras
cuando dejan de funcionar, o dejan de usarse son resinas como ABS (acrylonitrile butadiene
(aunque sean funcionales) o por obsolescencia styrene), HIPS (high-impact polystyrene), PPO
programada (Cucchiella et al., 2015). Este tér-
mino se utilizó por primera vez cuando Bernard
2 Global System for Mobile Communications
114
ENERLAC • Volumen IV. Número 2. Diciembre, 2020. ISSN: 2602-8042 (impreso) / 2631-2522 (digital).
OLADE – AUGM
(polyphenylene oxide), PPE (polyphenylene ether) (Adachi et al., 2013). La industria solar FV se
y PVC (polyvinyl chloride) (Babu et al., 2007) que encuentra exclusivamente dominada por la
pueden resultar altamente tóxicos. Las panta- tecnología convencional de semiconductores
llas y monitores de CRT, LCD y LED son inorgánicos como el silicio policristalino,
consideradas una sección altamente generadora monocristalino y las celdas de película delgada.
de RAEE. La composición de las pantallas o Esta tecnología actualmente representa un
televisores de CRT son de mayor volumen en porcentaje bajo de generación de RAEE a nivel
comparación con las de LCD. Por ejemplo, una mundial. De acuerdo con la Agencia Internacional
TV y un monitor de CRT pesan aproximadamente de Energía Renovable (Stephanie Weckend et
25 y 16 kilos respectivamente, en cambio las de al., 2016) menos del 1% de los paneles FV que
LCD tienen un peso menor de aproximadamente produjeron una potencia energética de 222 GW
60% (10 y 5 kilos, para TV y CRT respectiva- a finales del 2015 cumplieron su vida útil.
mente). Esto se puede observar al hacer Sin embargo, para el 2050 se proyecta una
la comparación en gramos/unidad de tres instalación de paneles FV que producirá
de los principales materiales (vidrio, plástico, 4,500 GW y por ende los desechos de paneles
acero/hierro) contenidos en los monitores solares se incrementará. La vida útil promedio
y/o pantallas de CRT y LCD. Por ejemplo, vidrio de un panel fotovoltaico es de hasta 30 años.
15,760 vs 590, plástico 8,755 vs 1,780, acero/ Tomando en cuenta este último dato y la
hierro 2,088 vs 2,530 en gramos/unidad para capacidad de desechos de paneles solares
CRT vs LCD (Cucchiella et al., 2015). Los acumulados para el 2050, se tendrán alrededor
paneles fotovoltaicos (FV) se presentan como de 60 a 70 millones de toneladas generadas de
una tecnología alterna viable para la genera- residuos de paneles FV (Stephanie Weckend et al.,
ción de energía limpia, debido a la intensa 2016); está estimación se puede observar en la
radiación que incide en la tierra, casi Figura 5, donde el Escenario de pérdida-regular:
127 mil TW de energía solar llegan a la tierra supone una vida útil de 30 años para los paneles
anualmente, de los cuales sólo 1.2 % se afirma solares, sin desgaste anticipado, mientras que el
que son aprovechables, es decir, 1,600 TW (Kabir Escenario de pérdida-temprana: tiene en cuenta
et al., 2018) o aproximadamente 104 veces más una vida útil “temprana”, “mediana edad” y fallas
de lo que la humanidad consume actualmente de “desgaste” antes de los 30 años.
100
90 4,500 GW 4500
78 millón 4000
80
70 3500
60 millón 3000
60
1,630 GW 2500
50
40 2000
30 1500
20 1000
270 GW
10 8 millón 500
43,500 250,000 1,7 millón
0 0
2016 2030 2050
Escenario de pérdida-regular Escenario de pérdida-temprana Lineal (Capacidad acumulada de paneles FV)
Fuente: Elaboración propia con datos basados del estudio de Stephanie Weckend et al. (2016) Definir pérdida-regular y
pérdida temprana
115
RESIDUOS DE APARATOS ELÉCTRICOS Y ELECTRÓNICOS (RAEE): IMPACTO SOCIAL, AMBIENTAL, GESTIÓN Y METODOLOGÍAS SOBRE SU MANEJO.
Becerra Paniagua, Dulce Kristal; Hernández Granados, Araceli; Díaz Cruz, Evelyn Betsabe; Cedano Villavicencio, Karla Graciela; Martínez Valencia, Horacio
Los RAEE contienen una combinación de mate- bromados (BFR, brominated flame retardants,
riales para hacer posible el funcionamiento y en inglés), que contienen grupos orgánicos
el diseño que el consumidor espera. Los tipos como: los bifenilos policlorados (BPCs), éteres
de RAEE se han clasificado en cinco categorías bifenílicos polibromados (PBDEs) e hidrocarburos
principalmente (Figura 6), las cuales son: i) me- aromáticos policíclicos (HAPs), entre otros. Así
tales-férreos, ii) plásticos, iii) metales no férreos, también algunos elementos tóxicos como: bario
iv) vidrio y v) otros materiales como madera, (Ba), berilio (Be), cadmio (Cd), cobalto (Co), co-
caucho, cartón. Los metales-férreos como el bre (Cu), cromo (Cu), hierro (Fe), plomo (Pb), litio
hierro y el acero son los materiales más comunes (Li), mercurio (Hg), níquel (Ni) y plata (Ag) (Kiddee
de encontrarse en los RAEE y ambos pueden et al., 2013). Lo cual genera una alta preocupa-
contribuir al 50% del peso total de los desechos; ción entre la comunidad científica y por lo que
el segundo grupo que más peso aporta es el sugiere documentarse con la literatura generada
plástico, constituye aproximadamente el 21%; sobre los RAEE, para conocer las propuestas que
el tercer grupo incluye a los metales preciosos se han hecho y poder sugerir procesos alternos
y representa un 13%; finalmente el porcentaje a los existentes para la gestión de estos residuos
total lo conforman otros materiales como el con metodologías y tecnologías que consideren
vidrio. Además, se han encontrado más de 1,000 principios de producción más limpia y economía
substancias tóxicas en los RAEE (Puckett et al., circular de los residuos que incluyan procesos
2002), entre ellos los retardadores de llamas ecológicos o que minimicen el daño generado.
Plásticos 15%
Metales 60%
Mezcla plástico-metal 5%
Contaminantes 3%
Cables 2%
Circuitos impresos 2%
Otros 1%
Fuente: Elaboración propia con datos adaptados del estudio de Widmer et al. (2005)
CIFRAS INTERESANTES SOBRE LOS RAEE cionan algunos datos que son preocupantes
y que tienen un potencial educativo para
Parte fundamental para conocer el daño concienciar acerca de los desechos gene-
ocasionado por los RAEE, es conocer en cifras el rados por el uso de aparatos eléctricos y
impacto que generan. A continuación, se men- electrónicos.
116
ENERLAC • Volumen IV. Número 2. Diciembre, 2020. ISSN: 2602-8042 (impreso) / 2631-2522 (digital).
OLADE – AUGM
Los RAEE representan el 2% de los desechos El reciclaje de los RAEE ofrece un ahorro
sólidos en el mundo, pero pueden llegar a de energía, ya que al no producir estos la
ser el 70% de los residuos peligrosos que se disminución de gases de efecto de invernadero
encuentran en los rellenos (ONU, 2019) y se es- será de uno de los efectos positivos más impor-
tima que para el 2040 los RAEE representen el tantes. Por ejemplo, al reciclar 10 kg de aluminio
14% de las emisiones totales de carbono (Ellen se estima que el ahorro de energía es del 90%,
MacArthur Foundation, 2018). Además, sólo se evitan aproximadamente 13kg de bauxita3,
el 20% de los RAEE son reciclados y el 80% 20 kg de dióxido de carbono (CO2) y 0.11 kg de
restante se incinera o termina en los rellenos dióxido de azufre (SO2) (Kumar et al., 2017).
sanitarios (World Economic Forum, 2019). La Otra ventaja de los RAEE es que el grado de
importancia del reciclaje de los RAEE tiene concentración de metales en los electrónicos
grandes beneficios ambientales y de salud, pero es mayor que los extraídos por la actividad
también su valor económico es considerable. minera. Por ejemplo, en promedio varios RAEE
Los RAEE al ser compuestos por más de 60 pueden contener 1,009 ppm de plata, 127 ppm
metales diferentes, se estima que su valor es de oro y 51.6 ppm de paladio, para el caso de
de 48 mil millones de euros. Bullion Street, la minería los valores son 215.5 ppm de plata,
analizó que la industria electrónica consume 1.01 ppm de oro y 2.7 ppm de paladio, lo cual es
320 toneladas de oro y 7,500 toneladas de plata una diferencia significativa (Kumar et al., 2017).
cada año y que la extracción urbana de desechos
electrónicos podría generar $ 21 mil millones DESTINO DE LOS RAEE
cada año (United Nations University, 2012). El
valor calculado de los RAEE con un tratamiento El código HS o código arancelario es un grupo
de reciclaje adecuado rebasa los 62,500,000,000 de números que se asignan de acuerdo con
de dólares anuales, lo cual supera al producto el sistema de clasificación estandarizado
interno bruto (PIB) de algunos países (UN, 2019). internacionalmente a las importaciones y ex-
Otro beneficio económico es que al reciclar portaciones entre países. Oficialmente se le
10,000 toneladas de desechos pueden crear los conoce como Sistema Armonizado o en inglés
siguientes empleos: 1 empleo para incineración, como Harmonized Commodity Description and
6 empleos en rellenos sanitarios, 36 empleos en Coding System o Harmonized System o sim-
centros de reciclaje y 296 empleos en el reciclaje plemente HS. Los códigos HS clasifican
de computadoras (ILSR, 2002). Aun cuando aproximadamente 8,000 productos, pero nin-
este es un problema global, actualmente sólo 67 guno de estos incluye una categoría para los
países cuentan con una legislación ambiental RAEE. Cuando un cargamento de este tipo
para el manejo adecuado de los RAEE (World ocurre, este solo tiene la opción de usar la
Economic Forum, 2019). clasificación de “electrónicos nuevos”. Alrededor
de 67 países cuentan con una legislación
De acuerdo con la Agencia de Protección Ambien- referente a los RAEE, lo cual comprende aproxi-
tal (EPA, Environmental Protection Agency, en madamente 2/3 del total de países del mundo
inglés) en EUA se desechan diariamente 416,000 (World Economic Forum, 2019). No obstante,
celulares y 142,000 computadoras, mientras existen regiones como África, América Latina y
que anualmente se producen 1,000,000,000 parte de Asia que no cuentan ninguna restricción
de celulares y 300 millones de computadoras sobre los RAEE. Por lo tanto, se vuelve más
(Chandler, 2017). Se ha estimado que al reciclar complicado hacer una estimación del volumen
un millón de celulares (148.4 toneladas) se
pueden obtener 34 kilos de oro, 350 kilos de 3 Roca compuesta en su mayoría por alúmina y, en menores
plata, 15 kilos de paladio y más de 16,000 kilos proporciones óxido de hierro y sílice. Es la fuente principal
de cobre (Kumar et al., 2017). de donde se extrae el aluminio.
117
RESIDUOS DE APARATOS ELÉCTRICOS Y ELECTRÓNICOS (RAEE): IMPACTO SOCIAL, AMBIENTAL, GESTIÓN Y METODOLOGÍAS SOBRE SU MANEJO.
Becerra Paniagua, Dulce Kristal; Hernández Granados, Araceli; Díaz Cruz, Evelyn Betsabe; Cedano Villavicencio, Karla Graciela; Martínez Valencia, Horacio
de intercambio o comercio de los RAEE de (P. Prakashan, 2019)] ocupe el primer lugar.
forma correcta. Por este motivo se han hecho Lo que es preocupante es observar que Europa
diversas investigaciones al respecto, por ejemplo, [10,180,000 km2 (Cañellas Blanco y Castellano
en el trabajo de Hopson y Puckett mencionan Alcaide, 2019)] siendo 4.3 veces menor del
que se rastrearon 205 DE con un sistema de tamaño de Asia, posea el segundo lugar en
posicionamiento global (GPS, Global Positioning generación de RAEE. También se muestra en
System, en inglés). Del total de estos, el 74% la tabla indexada el top 10 de países (gris con
(152 DE) se envió a un centro de reciclaje, el puntos negros) que más RAEE producen en
25% (51 DE) fueron donaciones y el 1% restante millones de toneladas métricas y la aportación
(2 DE) a reventa, se observó que de los 152 DE, calculada de kg/habitante (Baldé et al., 2017;
el 34% (69 DE) salieron del país en envíos ile- Prasad, 2020). Asimismo, se muestran (color
gales a países como Hong Kong, China o Taiwán. morado) las regiones que envían RAEE a otros
países en su mayoría en forma ilegal (flechas
En la Figura 7 se presenta el mapeo de los negras). Por ejemplo, América del Norte envía
desechos electrónicos, se describe (círculos a algunas regiones de América Latina y África;
amarillos) la cantidad en millones de toneladas Europa a regiones como China, Corea del Sur,
que se produce mundialmente, así como las la India, África. Finalmente, se ilustra también
aportaciones que generan cada continente. Si se la contraparte los países que más reciben RAEE
correlaciona el área de los continentes en km2 y en (color naranja) los cuales son países como
función de la cantidad de desechos que generan, México, Brasil y algunas regiones en Europa,
es de esperarse que Asia [44,579,000 km2 China, India y Tailandia.
Fuente: Elaboración propia con datos de World Economic Forum (2019); Baldé et al. (2017) y el uso de mapas de
MapChart (2020)
118
ENERLAC • Volumen IV. Número 2. Diciembre, 2020. ISSN: 2602-8042 (impreso) / 2631-2522 (digital).
OLADE – AUGM
IMPACTOS A LA SALUD Y AL AMBIENTE refleja en una pequeña escala todos los recur-
sos energéticos no renovables y renovables que
Existen investigaciones (Babu et al., 2007; se requieren para la fabricación de estos
Grant et al., 2013; Hopson y Puckett, 2016) que dispositivos (Rahmani et al., 2014).
han comprobado la enorme problemática que
conlleva el manejo inadecuado de desechos de Los contaminantes tóxicos encontrados en
DE en rellenos sanitarios, debido a que estos los RAEE se pueden clasificar en dos grupos:
se componen de elementos tóxicos como el Cd, 1) éteres de bifenilos poli bromados (PBDEs) y
Ni, Li, Cr, Cu, Pb, Sn, Hg, que generan problemas 2) metales pesados. Los PBDEs son popular-
ambientales y a la salud del ser humano al no mente usados como retardantes de llama.
saberlos manipular. Los RAEE impactan de ma- Diversas investigaciones han señalado a lo largo
nera directa a las personas que se encuentran de los años la presencia de PBDEs en sangre,
trabajando o viviendo cerca de los sitios de cabello y orina en personas por exposición
reciclaje o separación de éstos. La ingesta de estos directa o indirecta (Kiddee et al., 2013; Lechuga
contaminantes puede ocurrir por inhalación, Vázquez y Paredes Rizo, 2014). Este hallazgo
absorción dérmica o digestiva de vapores o es considerado un riesgo potencial de salud
partículas de compuestos tóxicos producidas debido a que los PBDEs están asociados a
cuando se incineran o cuando los contaminan- problemas de desestabilización del sistema
tes penetran en el suelo, ya que pueden filtrarse hormonal-tiroideo, tumores y toxicidad neuro-
al subsuelo, de manera que pueden llegar a las lógica, por mencionar algunos (Qu et al., 2007;
fuentes de consumo de agua (Kiddee et al., 2013). Rocha-Gutiérrez et al., 2015). Por otra parte, se
La población infantil, las mujeres embarazadas encuentran los metales pesados, se estima que
y los adultos mayores son los más vulnerables el 70% de los metales pesados que llegan a los
a este tipo de sustancias, ocasionando diver- rellenos en EUA provienen de los RAEE (Lines
sas enfermedades mencionadas en la Tabla 1. et al., 2016). Estos pueden provocar serias
afectaciones, como daños al bazo, riñones,
Por ejemplo, una computadora en promedio pesa dermatitis, asma y afectaciones neurodege-
29.6 kg, de lo cual el porcentaje en peso es 43.7% nerativas como el Alzheimer (aluminio) y el
metal, 23.3% es plástico, 17.3% componentes Parkinson (plomo) (Grant et al., 2013). En la
electrónicos y 15.7% es vidrio. Para producir una Tabla 1 se presenta un resumen de los tipos de
computadora, así como un monitor se necesitan contaminantes comúnmente encontrados en
1.5 toneladas de agua, 22 kilos de componentes RAEE, su aplicación en la industria y los impactos
y 240 kilos de combustibles fósiles, lo cual que pueden ocasionar en la salud.
119
RESIDUOS DE APARATOS ELÉCTRICOS Y ELECTRÓNICOS (RAEE): IMPACTO SOCIAL, AMBIENTAL, GESTIÓN Y METODOLOGÍAS SOBRE SU MANEJO.
Becerra Paniagua, Dulce Kristal; Hernández Granados, Araceli; Díaz Cruz, Evelyn Betsabe; Cedano Villavicencio, Karla Graciela; Martínez Valencia, Horacio
Fuente: Elaboración propia con datos del estudio de Grant (2013); Ramesh (2007); Puckett (2002); Cucchiella (2015)
120
ENERLAC • Volumen IV. Número 2. Diciembre, 2020. ISSN: 2602-8042 (impreso) / 2631-2522 (digital).
OLADE – AUGM
Tabla 2. Objetivos de desarrollo sostenible (ODS) y metas con relación a los RAEE
121
RESIDUOS DE APARATOS ELÉCTRICOS Y ELECTRÓNICOS (RAEE): IMPACTO SOCIAL, AMBIENTAL, GESTIÓN Y METODOLOGÍAS SOBRE SU MANEJO.
Becerra Paniagua, Dulce Kristal; Hernández Granados, Araceli; Díaz Cruz, Evelyn Betsabe; Cedano Villavicencio, Karla Graciela; Martínez Valencia, Horacio
ESTRATEGIAS PARA COMBATIR LA GENE- de 25 a 250 veces más oro comparado con la
RACIÓN DE RAEE cantidad encontrada en los minerales de los que
se extraen estos metales en la actividad minera
Como se mencionó anteriormente, la impor- (Tuncuk et al., 2012). Además, el costo de extraer
tancia del manejo correcto de los RAEE es una estos metales a partir de minerales es 13 veces
pieza fundamental para evitar los daños a la más alto en yacimientos naturales que en RAEE
salud pública y al ambiente (Dalrymple et al., (Zeng et al., 2018). Aun cuando esta actividad
2007). Este gran problema del manejo de RAEE tiene un alto impacto económico, los metales
puede tener diversas razones, una de ellas es en RAEE son reciclados y extraídos mediante
que los países desarrollados poseen un costo de herramientas de uso común como desarmado-
jornada laboral alta y varias regulaciones para res, martillos y cinceles. Algunos otros ejemplos
disposición final de estos residuos, por lo cual de extracción son las placas de circuitos im-
se exportan a países subdesarrollados. Estos a presos que son calentadas y los componentes
su vez, poseen herramientas y tecnologías más son removidos; el oro y otros metales se
rudimentarias para la extracción/reciclaje de extraen con baños de ácido; los plásticos son
los RAEE (Perkins et al., 2014) y es de esperarse triturados y se funden, y en el caso del cobre,
que estos sean recibidos en casi cualquier este se obtiene quemando el plástico que tiene
vertedero por su alto valor comercial (Tuncuk como cubierta. Todo lo anterior mencionado
et al., 2012). Por otra parte, se estima que el con escasa o nula protección personal, lo cual
40% de la composición de los RAEE son metales conlleva a una actividad de alto riesgo para la
de interés comercial, tales como cobre, plomo, integridad y salud del reciclador (Perkins et al.,
estaño, mercurio, plata, oro, platino y paladio. 2014). Es por ello, que en los últimos tiempos
Estos últimos, los metales preciosos, atraen la comunidad científica ha enfocado su atención
gran interés por su alto valor comercial. en el desarrollo de estrategias que incluyan
(Hagelüken y Corti, 2010). Aunque su extrac- métodos, procesos, protocolos y técnicas que
ción está limitada debido a la compleja no sólo permitan recuperar estos metales de
composición de metales en los desechos los desechos electrónicos, principalmente los
electrónicos, más del 80% del valor total de metales preciosos, sino, que consideren un
recuperación comprende metales preciosos buen uso de los desechos de DE, un manejo ade-
(Ag, Au y Pt) y más del 95% del total de metales cuado de estos y un reciclaje óptimo de RAEE
son básicos (Cu, Sn, Hg, Pb) (Diaz et al., 2017). con un enfoque de economía circular y amigables
Aunado a esto, los RAEE pueden contener con el ambiente.
122
ENERLAC • Volumen IV. Número 2. Diciembre, 2020. ISSN: 2602-8042 (impreso) / 2631-2522 (digital).
OLADE – AUGM
empresas
Diversas
reemplazarse
si ¿Se no si • Donar
puede reparar? ¿Sirve? • Vender a un
menor costo
no
Trasladar a un Clasiicar de
Desarmar y separar
Reciclar centro de reciclaje acuerdo al tipo
de componentes
de RAEE de RAEE
Fuente: Elaboración propia con base a los estudios de Cordova-Pizarro et al. (2019); Estrada-Ayub & Kahhat (2014)
123
RESIDUOS DE APARATOS ELÉCTRICOS Y ELECTRÓNICOS (RAEE): IMPACTO SOCIAL, AMBIENTAL, GESTIÓN Y METODOLOGÍAS SOBRE SU MANEJO.
Becerra Paniagua, Dulce Kristal; Hernández Granados, Araceli; Díaz Cruz, Evelyn Betsabe; Cedano Villavicencio, Karla Graciela; Martínez Valencia, Horacio
124
ENERLAC • Volumen IV. Número 2. Diciembre, 2020. ISSN: 2602-8042 (impreso) / 2631-2522 (digital).
OLADE – AUGM
125
RESIDUOS DE APARATOS ELÉCTRICOS Y ELECTRÓNICOS (RAEE): IMPACTO SOCIAL, AMBIENTAL, GESTIÓN Y METODOLOGÍAS SOBRE SU MANEJO.
Becerra Paniagua, Dulce Kristal; Hernández Granados, Araceli; Díaz Cruz, Evelyn Betsabe; Cedano Villavicencio, Karla Graciela; Martínez Valencia, Horacio
AGRADECIMIENTOS
CONFLICTO DE INTERÉS
126
ENERLAC • Volumen IV. Número 2. Diciembre, 2020. ISSN: 2602-8042 (impreso) / 2631-2522 (digital).
OLADE – AUGM
REFERENCIAS
Adachi, M., Tanino, R., Adachi, J., Mori, Y., Tsuchiya, K.,
Isoda, S. and Uchida, F. (2013). Verification of necessity
of one-dimensional titania nanoscale materials for dye-
sensitized solar cells. Journal of Power Sources, 226, 94–100.
https://doi.org/10.1016/j.jpowsour.2012.10.024
Baldé, C. P., Forti, V., Gray, V., Kuehr, R. and Stegmann, P. (2017).
The Global E-waste Monitor 2017: Quantities, Flows, and
Resources. United Nations University (UNU), International
Telecommunication Union (ITU) & International Solid
Waste Association (ISWA), Bonn/Geneva/Vienna. https://
www.itu.int/en/ITU-D/Climate-Change/Documents/GEM
2017/Global-E-waste Monitor 2017 .pdf
127
RESIDUOS DE APARATOS ELÉCTRICOS Y ELECTRÓNICOS (RAEE): IMPACTO SOCIAL, AMBIENTAL, GESTIÓN Y METODOLOGÍAS SOBRE SU MANEJO.
Becerra Paniagua, Dulce Kristal; Hernández Granados, Araceli; Díaz Cruz, Evelyn Betsabe; Cedano Villavicencio, Karla Graciela; Martínez Valencia, Horacio
Cañellas Blanco, A. M. y Castellano Alcaide, J. M. (2019). Daven, J. I. and Klein, R. N. (2008). Progress in waste
Ámbito Lingüístico y Social. Geografía Humana. Editex, S.A. management research. Nova Science, Publishers, Inc, pp.
ISBN: 978-84-9161-998-7 93–102.
Cesaro, A., Belgiorno, V., Vaccari, M., Jandric, A., Chung, T. D., Dimitrakakis, E., Janz, A., Bilitewski, B. and Gidarakos,
Dias, M. I., Hursthouse, A. and Salhofer, S. (2018). A device- E. (2009). Small WEEE: Determining recyclables and
specific prioritization strategy based on the potential hazardous substances in plastics. Journal of Hazardous
for harm to human health in informal WEEE recycling. Materials, 161(2–3), 913–919. https://doi.org/10.1016/j.
Environmental Science and Pollution Research, 25(1), jhazmat.2008.04.054
683–692. https://doi.org/10.1007/s11356-017-0390-7
Ellen MacArthur Foundation. (2018). The circular economy
Chandler, D. (2017). Strategic Corporate Social Respon- opportunity for urban & industrial innovation in China.
sibility: Sustainable Value Creation. https://books.google. Circular Economy Perspectives Series, 1–166. https://
com.mx/books?id=q4CzDAAAQBAJ&printsec=frontcover www.ellenmacarthurfoundation.org/assets/downloads/
&hl=es&source=gbs_ge_summary_r&cad=0#v=onepage The-circular-economy-opportunity-for-urban-industrial-
&q&f=false innovation-in-China_19-9-18_1.pdf
Chen, M., Wang, J., Chen, H., Ogunseitan, O. A., Zhang, M., Estrada-Ayub, J. A. and Kahhat, R. (2014). Decision
Zang, H. and Hu, J. (2013). Electronic waste disassembly with factors for e-waste in Northern Mexico: To waste or trade.
industrial waste heat. Environmental Science and Technology, Resources, Conservation and Recycling, 86, 93–106. https://
47(21), 12409–12416. https://doi.org/10.1021/es402102t doi.org/10.1016/j.resconrec.2014.02.012
Chu, Y., Chen, M., Chen, S., Wang, B., Fu, K. and Chen, H. (2015). Estrada-Ruiz, R. H., Flores-Campos, R., Gámez-Altamirano,
Micro-copper powders recovered from waste printed circuit H. A. and Velarde-Sánchez, E. J. (2016). Separation of the
boards by electrolysis. Hydrometallurgy, 156, 152–157. metallic and non-metallic fraction from printed circuit
https://doi.org/10.1016/j.hydromet.2015.06.006 boards employing green technology. Journal of Hazardous
Materials, 311, 91–99. https://doi.org/10.1016/j.jhazmat.
Cordova-Pizarro, D., Aguilar-Barajas, I., Romero, D. and 2016.02.061
Rodriguez, C. (2019). Circular economy in the electronic
products sector: material flow analysis and economic Fogarasi, S., Imre-Lucaci, F., Ilea, P. and Imre-Lucaci, Á. (2013).
impact of cellphone e-waste in Mexico. Sustainability, 11(5), The environmental assessment of two new copper recovery
1361. https://doi.org/10.3390/su11051361 processes from Waste Printed Circuit Boards. Journal of
Cleaner Production, 54, 264–269. https://doi.org/10.1016
Cucchiella, F., D’Adamo, I., Lenny Koh, S. C. and Rosa, P. (2015). /j.jclepro.2013.04.044
Recycling of WEEEs: An economic assessment of present
and future e-waste streams. Renewable and Sustainable Fondo Monetario Internacional. (2019). Perspectivas de la
Energy Reviews, 51, 263–272. https://doi.org/10.1016/j. economía mundial. Desaceleración mundial de la actividad
rser.2015.06.010 manufacturera, crecientes barreras comerciales. ISSN 1564-
5215 (online)
Culver, J. (2015). The life cycle of the CPU. http://www.
cpushack.com/life-cycle-of-cpu.html Grant, K., Goldizen, F. C., Sly, P. D., Brune, M. N., Neira, M., van
den Berg, M. and Norman, R. E. (2013). Health consequences
Dalrymple, I., Wright, N., Kellner, R., Bains, N., Geraghty, K., of exposure to e-waste: A systematic review. The Lancet
Goosey, M. and Lightfoot, L. (2007). An integrated approach Global Health, 1(6), e350–e361. https://doi.org/10.1016/
to electronic waste (WEEE) recycling. Circuit World, 33(2), S2214-109X(13)70101-3
52–58. https://doi.org/10.1108/03056120710750256
128
ENERLAC • Volumen IV. Número 2. Diciembre, 2020. ISSN: 2602-8042 (impreso) / 2631-2522 (digital).
OLADE – AUGM
Hagelüken, C. and Corti, C. W. (2010). Recycling of gold from Kaya, M. (2016). Recovery of metals and nonmetals from
electronics: Cost-effective use through “design for recycling.” electronic waste by physical and chemical recycling
Gold Bulletin, 43(3), 209–220. https://doi.org/10.1007/ processes. In Waste Management, 57, 64–90. Elsevier Ltd.
BF03214988 https://doi.org/10.1016/j.wasman.2016.08.004
Havlik, T., Orac, D., Petranikova, M., Miskufova, A., Kukurugya, Kiddee, P., Naidu, R. and Wong, M. H. (2013). Electronic
F. and Takacova, Z. (2010). Leaching of copper and tin from waste management approaches: An overview. Waste
used printed circuit boards after thermal treatment. Journal Management, 33(5), 1237–1250. https://doi.org/10.1016/j.
of Hazardous Materials, 183(1–3), 866–873. https://doi. wasman.2013.01.006
org/10.1016/j.jhazmat.2010.07.107
Kim, M., Jang, Y. C. and Lee, S. (2013). Application of Delphi-
Hopson, E. and Puckett, J. (2016). Scam Recycling: e-Dumping AHP methods to select the priorities of WEEE for recycling
on Asia by US Recyclers, Basel Action Network, USA. in a waste management decision-making tool. Journal of
https://wiki.ban.org/images/1/12/ScamRecycling Environmental Management, 128, 941–948. https://doi.
Report-web.pdf org/10.1016/j.jenvman.2013.06.049
ILSR. (2002, February). Recycling Means Business - Institute Kumar, A., Holuszko, M. and Espinosa, D. C. R. (2017).
for Local Self-Reliance. https://ilsr.org/recycling-means- E-waste: An overview on generation, collection, legislation
business/ and recycling practices. Resources, Conservation and
Recycling, 122, 32–42. https://doi.org/10.1016/j.resconrec.
Imre-Lucaci, F., Dorneanu, S. A. and Ilea, P. (2009). 2017.01.018
Electrochemical metals recovery from electronic wastes.
Part. I. copper recovery from synthetic solutions. Studia Lechuga Vázquez, P. y Paredes Rizo, M. . L. (2014).
Universitatis Babes-Bolyai Chemia, SPEC. ISSUE 1, 105–111. Efectos biológicos derivados de la exposición a PBDES en
https://www.researchgate.net/profile/Sorin_Dorneanu/ trabajadores del reciclaje de e-waste: revisión sistemática.
publication/294494630_Electrochemical_metals_recovery_ Medicina y Seguridad Del Trabajo, 60(237), 685–713.
from_electronic_wastes_Part_I_copper_recovery_from_ https://doi.org/10.4321/s0465-546x2014000400008
synthetic_solutions/links/56d5cd5208aee1aa5f 730b22/
Electrochemical-metals-recovery-from-electronic-wastes- Leirnes, J. S. and Lundstrom, M. S. (1983). Method for
Part-I-copper-recovery-from-synthetic-solutions.pdf working-up metal-containing waste products. https://www.
google.com/patents/US4415360
Instituto Nacional de Ecología INECC. (2006). Diagnóstico
sobre la Generación de Residuos Electrónicos en México. Lines, K., Garside, B., Sinha, S. and Fedorenko, I. (2016). Clean
SEMARNAT. http://www2.inecc.gob.mx/descargas/sqre/ and inclusive? Recycling e-waste in China and India Issue
res_electronicos_borrador_final.pdf paper Sustainable markets (Issue March). https://pubs.iied.
org/pdfs/16611IIED.pdf
Jadhao, P., Chauhan, G., Pant, K. K. and Nigam, K. D. P. (2016).
Greener approach for the extraction of copper metal from London, B. (1932). Ending the depression through
electronic waste. Waste Management, 57, 102–112. https:// planned obsolescence. Self-Published, 1–8. https://faculty.
doi.org/10.1016/j.wasman.2015.11.023 bemidjistate.edu/mlawrence/London.pdf
Kabir, E., Kumar, P., Kumar, S., Adelodun, A. A. and Kim, K. Long, E., Kokke, S., Lundie, D., Shaw, N., Ijomah, W. and Kao,
H. (2018). Solar energy: Potential and future prospects. In C. chuan. (2016). Technical solutions to improve global sus-
Renewable and Sustainable Energy Reviews, 82, 894–900. tainable management of waste electrical and electronic
Elsevier Ltd. https://doi.org/10.1016/j.rser.2017.09.094 equipment (WEEE) in the EU and China. Journal of Rema-
nufacturing, 6(1), 1–27. https://doi.org/10.1186/s13243-
015-0023-6
129
RESIDUOS DE APARATOS ELÉCTRICOS Y ELECTRÓNICOS (RAEE): IMPACTO SOCIAL, AMBIENTAL, GESTIÓN Y METODOLOGÍAS SOBRE SU MANEJO.
Becerra Paniagua, Dulce Kristal; Hernández Granados, Araceli; Díaz Cruz, Evelyn Betsabe; Cedano Villavicencio, Karla Graciela; Martínez Valencia, Horacio
MapChart. (2020). World Map - Detailed - MapChart. https:// Rahmani, M., Nabizadeh, R., Yaghmaeian, K., Mahvi, A. H. and
mapchart.net/detworld.html Yunesian, M. (2014). Estimation of waste from computers
and mobile phones in Iran. Resources, Conservation and
Ongondo, F. O., Williams, I. D. and Cherrett, T. J. (2011). How Recycling, 87, 21–29. https://doi.org/10.1016/j.resconrec.
are WEEE doing? A global review of the management of 2014.03.009
electrical and electronic wastes. Waste Management, 31(4),
714–730. https://doi.org/10.1016/j.wasman.2010.10.023 Rocha-Gutiérrez, B. A., Peralta-Pérez, M. D. R. y Zavala-Díaz
de la Serna, F. J. (2015). Revisión global de los contaminantes
ONU. (2019). Los desechos electrónicos, una oportunidad de emergentes PBDE y el caso particular de México. In
oro para el trabajo decente. https://news.un.org/es/story/ Revista Internacional de Contaminación Ambiental, 31(3),
2019/04/1455621 311–320. http://www.scielo.org.mx/scielo.php?script=
sci_arttext&pid=S0188-49992015000300010
Perkins, D. N., Brune Drisse, M. N., Nxele, T. and Sly, P. D. (2014).
E-waste: A global hazard. Annals of Global Health, 80(4), Rubin, R. S., Castro, M. A. S. De, Brandão, D., Schalch, V. and
286–295. https://doi.org/10.1016/j.aogh.2014.10.001 Ometto, A. R. (2014). Utilization of Life Cycle Assessment
methodology to compare two strategies for recovery of
Petter, P. M. H., Veit, H. M. and Bernardes, A. M. (2014). copper from printed circuit board scrap. Journal of Cleaner
Evaluation of gold and silver leaching from printed circuit Production, 64, 297–305. https://doi.org/10.1016/j.
board of cellphones. Waste Management, 34(2), 475–482. jclepro.2013.07.051
https://doi.org/10.1016/j.wasman.2013.10.032
Smedley, T. (2020). How to mine precious metals in
Prakashan, P. (2019). General Knowledge 2019 (Prabhat your home. BBC Future. https://www.bbc.com/future/
Prakashan (ed.)). ISBN : 9789352665976 article/20200407-urban-mining-how-your-home-may-be-
a-gold-mine
Prasad, M. N. V. (2020). Handbook of Electronic Waste
Management International Best Practices and Case Studies. Stenvall, E., Tostar, S., Boldizar, A., Foreman, M. R. S. J. and
Butterworth-Heinemann. eBook ISBN: 9780128170311 Möller, K. (2013). An analysis of the composition and
metal contamination of plastics from waste electrical and
Puckett, J., Byster, L., Westervelt, S., Gutierrez, R., Davis, S., electronic equipment (WEEE). Waste Management, 33(4),
Puckett, J. and Smith, T. (2002). Exporting Harm The High- 915–922. https://doi.org/10.1016/j.wasman.2012.12.022
Tech Trashing of Asia. Svt C. http://svtc.org/wp-content/
uploads/technotrash.pdf Stephanie Weckend, Andreas Wade and Garvin Heath.
(2016). End of Life Management Solar PV Panels.
Qu, W., Bi, X., Sheng, G., Lu, S., Fu, J., Yuan, J. and Li, L. (2007). https://www.irena.org/documentdownloads/publications
Exposure to polybrominated diphenyl ethers among /irena_ieapvps_end-of-life_solar_pv_panels_2016.pdf
workers at an electronic waste dismantling region in Guang-
dong, China. Environment International, 33(8), 1029–1034. Swain, B., Mishra, C., Kang, L., Park, K. S., Lee, C. G. and Hong,
https://doi.org/10.1016/j.envint.2007.05.009 H. S. (2015). Recycling process for recovery of gallium
from GaN an e-waste of LED industry through ball milling,
Queiruga, D., Walther, G., González-Benito, J. and Spengler, annealing and leaching. Environmental Research, 138
T. (2008). Evaluation of sites for the location of WEEE (2015), 401–408. https://doi.org/10.1016/j.envres.2015.
recycling plants in Spain. Waste Management, 28(1), 181– 02.027
190. https://doi.org/10.1016/j.wasman.2006.11.001
130
ENERLAC • Volumen IV. Número 2. Diciembre, 2020. ISSN: 2602-8042 (impreso) / 2631-2522 (digital).
OLADE – AUGM
Tuncuk, A., Stazi, V., Akcil, A., Yazici, E. Y. and Deveci, H. Worldometer. (2020). Population by Country (2020) -
(2012). Aqueous metal recovery techniques from e-scrap: Worldometer. https://www.worldometers.info/world-
Hydrometallurgy in recycling. Minerals Engineering, 25(1), population/population-by-country/
28–37. https://doi.org/10.1016/j.mineng.2011.09.019
WOS. (2020). Web of Science [v.5.33] - Colección principal de
Turner, A. (2020). How many people have smartphones in Web of Science Búsqueda básica. www.webofknowledge.com
2020? Bank My Cell. https://www.bankmycell.com/blog/
how-many-phones-are-in-the-world Yang, H., Ma, M., Thompson, J. R. and Flower, R. J. (2018).
Waste management, informal recycling, environmental
UN. (2019). UN report: Time to seize opportunity, tackle pollution and public health. Journal of Epidemiology
challenge of e-waste | UNEP - UN Environment Programme. and Community Health, 72(3), 237–243. https://doi.
https://www.unenvironment.org/news-and-stories/ org/10.1136/jech-2016-208597
press-release/un-report-time-seize-opportunity-tackle-
challenge-e-waste Yue, C., Sun, H., Liu, W.-J., Guan, B., Deng, X., Zhang, X. and
Yang, P. (2017). Environmentally benign, rapid, and selective
United Nations University. (2012). E-waste: Annual Gold, extraction of gold from ores and waste electronic materials.
Silver “Deposits” in New High-Tech Goods Worth $21 Billion+; Angewandte Chemie, 129(32), 9459–9463. https://doi.
Less Than 15% Recovered - United Nations University. org/10.1002/ange.201703412
https://unu.edu/media-relations/releases/step-news-
release-6-july-2012-e-waste-precious-metals-recovery.html
Wong, M. H., Wu, S. C., Deng, W. J., Yu, X. Z., Luo, Q., Leung,
A. O. W., Wong, C. S. C., Luksemburg, W. J. and Wong, A. S.
(2007). Export of toxic chemicals - A review of the case of
uncontrolled electronic-waste recycling. Environmental
Pollution, 149(2), 131–140. https://doi.org/10.1016/j.
envpol.2007.01.044
131
CONTRIBUCIÓN DE LA CADENA PRODUCTIVA DE LA
ENERGÍA EN LAS EMISIONES DE GASES DE EFECTO
INVERNADERO EN ARGENTINA
CONTRIBUTION OF THE ENERGY VALUE CHAIN IN
ARGENTINA’S GREENHOUSE GAS EMISSIONS
Maxime Le Bail 1
132
ENERLAC • Volumen IV. Número 2. Diciembre, 2020. ISSN: 2602-8042 (impreso) / 2631-2522 (digital).
OLADE – AUGM
RESUMEN
133
CONTRIBUCIÓN DE LA CADENA PRODUCTIVA DE LA ENERGÍA EN LAS EMISIONES DE GASES DE EFECTO INVERNADERO EN ARGENTINA.
Le Bail, Maxime
ABSTRACT INTRODUCCIÓN
This paper gives an overview of the emissions of Argentina ha asumido compromisos ante la
greenhouse gases (GHG) in Argentina, detailing comunidad internacional en materia de combate
their origin by sector; then, it analyzes the al cambio climático al presentar su Contribución
Intended Nationally Determined Contributions Prevista y Determinada a Nivel Nacional (INDC,
(INDC) presented to the United Nations Frame- por sus siglas en inglés) en octubre de 2015, antes
work Convention on Climate Change (UNFCCC) de la vigésimo primera reunión de la Convención
and the feasibility of their compliance, and de las Partes celebrada en París donde se firmó el
studies the implications of the coronavirus Acuerdo homónimo, posteriormente ratificado y
pandemic Covid-19 on the country’s emissions. promulgado por el país mediante la Ley 27270.
In a transversal way, it emphasizes on the energy
sector and demonstrates its wide contribution Siendo un país con una matriz energética do-
to the country’s total GHG. This is done by anali- minada casi exclusivamente por los recursos
zing the official communications of the country, fósiles, y cuyo consumo energético final es
of the UNFCCC, of the International Energy compuesto en un 75% de derivados de petróleo
Agency and of the United Nations Environment y gas natural (Agencia Internacional de Energía,
Programme, as well as several research papers 2020), las emisiones del sector energético en
and articles in the medias. Argentina representan poco más de la mitad de las
emisiones totales (Moreira Muzio, M., 2019), por
lo que el cumplimiento de dichos compromisos
Keywords: Climate Change, Greenhouse Gases, implica inevitablemente la reducción de las
Intended Nationally Determined Contributions, emisiones de gases de efecto invernadero (GEI)
Energy, INDC, Covid-19, Argentina. en ese sector.
134
ENERLAC • Volumen IV. Número 2. Diciembre, 2020. ISSN: 2602-8042 (impreso) / 2631-2522 (digital).
OLADE – AUGM
Entonces, ¿qué medidas considera Argentina más GEI emite. Entre 1990 y 2016, sus emi-
para disminuir sus emisiones en el sector ener- siones de GEI aumentaron un 35%.
gético? ¿Serán suficientes para cumplir a cabalidad
con los compromisos adquiridos? ¿La pandemia De esos 364 millones de toneladas, 232 MtCO2e
de Covid-19 habrá contribuido lo suficiente en la son de dióxido de carbono (CO2), representando
reducción de emisiones? cerca del 64% de las emisiones totales del país,
mientras que 79 MtCO2e son de metano (CH4),
Este ensayo presenta las emisiones de GEI de 48 MtCO2e de óxido nitroso (N2O), 5 MtCO2e de
Argentina, analiza la INDC comprometida por Hidrofluorocarbono, y menos de un millón de
el país ante la Convención Marco de Naciones toneladas de CO2e de perfluorocarbonos y
Unidas sobre Cambio Climático (CMNUCC) y hexafluoruro de azufre (Secretaría de Ambiente
la factibilidad de que sea cumplida, y estudia y Desarrollo Sustentable, 2019).
las implicaciones de la actual pandemia de
coronavirus Covid-19 sobre las emisiones del Con una población total de 43.51 millones de
país. De manera transversal, hace énfasis sobre habitantes (Ourworldindata.org, 2020) y un
la contribución del sector energético. Producto Interno Bruto (PIB) de alrededor de
447 mil millones de dólares (Ourworldindata.
EMISIONES DE GASES DE EFECTO INVER- org, 2020), ello representa la emisión de 8.37
NADERO EN ARGENTINA toneladas de CO2 equivalente per cápita y 814
tCO2e por millón de dólares de PIB.
La Secretaría de Ambiente y Desarrollo Susten-
table de Argentina reporta que 364 millones de Poco más de la mitad de dichas emisiones (el
toneladas de dióxido de carbono equivalente 53%) fueron emitidas por el sector energía con
(MtCO2e) de GEI fueron emitidas en 2016 en 193 MtCO2e emitidas, mientras que el 37% fue
el país (Moreira Muzio, M., 2019). Lo anterior emitido por el sector agropecuario (136 MtCO2e),
representa cerca del 0.7% de las emisiones el 6% por los procesos industriales y el uso de
totales a nivel mundial (UNEP, 2016) y posiciona productos (20 MtCO2e), y el 4% por los residuos
a Argentina como el vigésimo noveno país que (15 MtCO2e).
Figura 1. Emisiones de gases de efecto invernadero (GEI) por sector en Argentina en 2016
15 MtCO2e
20 MtCO2e 4%
6%
Energía
Agropecuario
Residuos
135
CONTRIBUCIÓN DE LA CADENA PRODUCTIVA DE LA ENERGÍA EN LAS EMISIONES DE GASES DE EFECTO INVERNADERO EN ARGENTINA.
Le Bail, Maxime
El sector energía emitió 193 MtCO2e en 2016, El sector agropecuario emitió 136 MtCO2e en
representando poco más de la mitad de las emi- 2016, representando el 37% de las emisiones tota-
siones totales del país durante dicho año (Moreira les del país en ese año (Moreira Muzio, M., 2019).
Muzio, M., 2019).
En este sector, se incluyen la gestión de ganado
Este sector incluye todas las emisiones de GEI vivo y de estiércol, que representan el 41% de las
que emanan de la combustión de combustibles emisiones de este sector, así como las emisiones
fósiles, que representan el 94% de las emisiones de los suelos gestionados y las emisiones de la
de este sector, y las emisiones fugitivas, que repre- aplicación de fertilizantes, que representan el
sentan el 6% restante. 36% de las emisiones, y las emisiones y absor-
ciones de tierras forestales, tierras de cultivo,
Dentro de las emisiones que emanan de la com- pastizales, y otros tipos de uso de la tierra, que re-
bustión de combustibles fósiles, cabe resaltar presentan el 23% de las emisiones en este sector.
la fuerte participación de la industria de la
energía: la refinación del petróleo, la fabricación Sus emisiones han oscilado constantemente,
de combustibles sólidos y otras industrias ener- conociendo un aumento de +33% entre 1990
géticas, y la producción de electricidad y calor y 2000 al pasar de emitir 158.8 MtCO2e a
representan el 33% de las emisiones en este 211 MtCO2e, debido particularmente a una
sector; por su parte, el transporte representa el fuerte intensificación de la agricultura en el país.
26% de las emisiones de este sector, mientras A raíz de la crisis de 2001, la producción bajó
que la quema de combustibles de uso residencial, hasta llegar a emitir 177.7 MtCO2e en 2010,
comercial y agropecuario representa el 18%, antes de aumentar nuevamente en 2012 a
delante de la quema de combustibles de la in- 194.8 MtCO2e. Desde 2012, las emisiones de este
dustria manufacturera y de la construcción que sector están disminuyendo de manera continua
representa el 17%. A lo definido por Moreira (Moreira Muzio, M., 2019). Lo anterior revela, de
Muzio, M. (2019), faltaría agregar la extracción de manera global, una disminución de 22.3 MtCO2e,
hidrocarburos que genera amplias emisiones. cerca del 15%, en las emisiones de este sector
entre 1990 y 2016.
Este sector duplicó sus emisiones en las tres
últimas décadas, conociendo un aumento de Los GEI en este sector son, en orden de im-
97.9 MtCO2e en ese periodo, al pasar de portancia, el CH4 (43%) debido a la crianza bovina
95.5 MtCO2e en 1990 a 193 MtCO2e en 2016. principalmente, el N2O (34%) debido al uso de
El aumento fue particularmente fuerte a partir del fertilizantes, y el CO2 (23%) debido al cambio
año 2000, año que corresponde tanto al inicio de uso de suelo para los cultivos (Secretaría de
del aumento de la producción de gas natural en Ambiente y Desarrollo Sustentable, 2019).
el país, como del consumo de derivados de petró-
leo, gas natural y electricidad (Moreira Muzio, Procesos industriales y uso de productos
M., 2019).
Los procesos industriales y el uso de productos
Según la Secretaría de Ambiente y Desarrollo emitieron 20 MtCO2e en 2016, representando
Sustentable (2019), el 96% de las emisiones en aproximadamente el 6% de las emisiones totales
este sector fueron de CO2, debido a la quema de de GEI del país en ese año (Moreira Muzio,
los combustibles fósiles. M., 2019).
136
ENERLAC • Volumen IV. Número 2. Diciembre, 2020. ISSN: 2602-8042 (impreso) / 2631-2522 (digital).
OLADE – AUGM
Este sector incluye todas las emisiones de Este sector duplicó sus emisiones en las tres
GEI generadas como resultado de la reacción últimas décadas, pasando de 7.6 MtCO2e en 1990,
entre materias primas empleadas en diferentes a 15.4 MtCO2e en 2016. Su aumento ha sido
procesos químicos. constante sobre este periodo, a una velocidad de
0.3 MtCO2e por año (Moreira Muzio, M., 2019).
Las emisiones de este sector provienen en un
35% de la industria de los minerales, así como Los GEI generados en este sector son princi-
de la industria de los metales que representa el palmente el CH4 (94%) y el N2O (6%) (Secretaría
26% de las emisiones en este sector, mientras de Ambiente y Desarrollo Sustentable, 2019).
que el uso de productos y la industria química
representan respectivamente el 26% y el 13%.
El sector multiplicó sus emisiones por 2.4 desde Así, si bien el sector procesos industriales y el
1990, conociendo un aumento de 11.7 MtCO2e uso de productos es el que conoció el mayor
en ese periodo al pasar de emitir 8.4 MtCO2e a aumento, en cifras brutas es el sector energía que
20.1 MtCO2e. Dicho aumento fue relativamente mayor aumento conoció en el periodo 1990-2016.
constante, a una tasa de 0.45 MtCO2e por año, Ante este panorama, Argentina ha desarrollado
conociendo su mayor aceleración en el periodo algunos compromisos de reducción de emisiones
post-crisis 2001 debido a la reactivación de la que fueron presentados algunas semanas antes
economía y del sector secundario (Moreira Muzio, de la Vigésimo Primera Conferencia de las Partes
M., 2019). (COP 21).
137
CONTRIBUCIÓN DE LA CADENA PRODUCTIVA DE LA ENERGÍA EN LAS EMISIONES DE GASES DE EFECTO INVERNADERO EN ARGENTINA.
Le Bail, Maxime
Entrega de la INDC
En la primera entrega de la INDC (República 2030 con respecto de su BAU, es decir que
Argentina, 2015), Argentina presentó un esce- Argentina llegaría a emitir 570 MtCO2e en ese
nario de línea de base Business as usual (BAU) año. Aunado a la meta incondicional, Argentina
con una proyección que inicia en el año 2005 propuso reducir de manera condicionada el 30%
con 450 MtCO2e y llega al año 2030 con de las emisiones de GEI para 2030 de recibir
emisiones de GEI equivalentes a 670 MtCO2e. financiamiento internacional, transferencia de
La meta incondicional de reducción pretendía tecnologías y capacitación adecuada, entre otros;
la reducción de 15% de emisiones para el año el país estaría entonces emitiendo 470 MtCO2e.
Figura 2. Proyección de emisiones de Argentina al año 2030 según un escenario de línea de base BAU,
y la reducción de emisiones gracias a la meta incondicional y la meta condicionada
700
650
600
550
500
450
400
350
2005
2006
2007
2008
2009
2010
2011
2012
2013
2014
2015
2016
2017
2018
2019
2020
2021
2022
2023
2024
2025
2026
2027
2028
2029
Línea Base BAU (MtCO2e) 2030
En la primera revisión de su INDC (República significaría que el país no deberá emitir más que
Argentina, 2016), Argentina presenta una nueva 483 MtCO2e en 2030. Aunado a lo anterior, de
línea de base BAU que empieza en alrededor de manera condicionada, propone una reducción
300 MtCO2e en 1990 y alcanza 592 MtCO2e en de hasta 37%, equivalente a una reducción de
2030. De manera incondicional, se compromete 223 MtCO2e; el país estaría emitiendo entonces
a realizar una reducción de 18% equivalente 369 MtCO2e en 2030, una cifra similar a lo que se
a 109 MtCO2e sobre esta proyección, lo que emitió en 2016.
138
ENERLAC • Volumen IV. Número 2. Diciembre, 2020. ISSN: 2602-8042 (impreso) / 2631-2522 (digital).
OLADE – AUGM
Figura 3. Proyección de emisiones de Argentina al año 2030 según un escenario de línea de base BAU,
y la reducción de emisiones gracias a la meta incondicional y la meta condicionada, derivado de la
primera revisión de su INDC
700
600 592
Reducción por medidas
incondicionales
109 MtCO2e
400 369
114 MtCO2e
1
2005
2030
1990
1992
1994
1996
1998
2000
2002
2004
2006
2008
2010
2012
2014
2016
2018
2020
2022
2024
2026
2028
Ahorros proyectados de modo Ahorros adicionales proyectados por medidas Emisiones proyectadas como resultado de la
incondicional sobre línea de base condicionales sobre línea de base implementación efectiva de todos los ahorros
proyectados (condicional e incondicional)
Cabe resaltar las grandes diferencias entre medidas más ambiciosas” habrán resultado
el primer escenario presentado en 2015 y el únicamente en la reducción adicional de 9
segundo escenario presentado en 2016 donde MtCO2e. De igual manera, la revisión de las metas
existe una discrepancia de 50 MtCO2e para el condicionadas aportó únicamente la reducción
año 2005 (datos históricos) y de hasta 78 MtCO2e adicional de 19 MtCO2e.
para la proyección al año 2030. Por ende, la
declaración relativa a que su contribución fue Aunada a la tímida revisión de sus compromisos,
mejorada “planificando medidas de mitigación Argentina no explicita las medidas que tomará
incondicionales que logran bajar su meta al para lograr dicha reducción; sólo menciona
2030 de 570 a 483 millones de tCO2eq” debe lo siguiente: “la meta se logrará a través de la
ser moderada; Argentina reconoce que ello es implementación de una serie de medidas a lo
principalmente debido al “cambio a metodología” largo de la economía, focalizando en los sectores
utilizada para estimar las emisiones, entre la de energía, agricultura, bosques, transporte,
INDC y su primera revisión, lo que permitió industria y residuos”.
“mejorar la calidad del inventario evitando la
sobreestimación de las emisiones agrícolas”. En cambio, existe mucho mayor grado de detalle
Finalmente, la “revisión de más de 50 medidas respecto a las acciones de adaptación que el
incondicionales y la incorporación de nuevas país se compromete a emprender: “ampliación
139
CONTRIBUCIÓN DE LA CADENA PRODUCTIVA DE LA ENERGÍA EN LAS EMISIONES DE GASES DE EFECTO INVERNADERO EN ARGENTINA.
Le Bail, Maxime
de las redes de monitoreo, el fortalecimiento de 2017) permitió no sólo compensarlo sino también
los sistemas de alerta temprana y los servicios dar un rumbo para lograr el compromiso de no
climáticos para la salud, seguridad alimentaria, rebasar 483 MtCO2e en 2030.
el agua, la energía y la reducción de riesgos de
desastres”, “impulso a proyectos de I+D y de Dicho Plan presenta una serie de políticas,
tecnologías como así también de buenas prácticas medidas y acciones para a) abastecer de energía
productivas”, “mapeo de vulnerabilidades y limpia y sostenible, b) acompañar el crecimiento
riesgos climáticos como herramienta de diag- productivo y poblacional, c) fomentar el uso
nóstico del estado de situación y apoyo a la responsable de la energía, y d) promover la
gestión de la adaptación al cambio climático”, eficiencia energética. Tiene la pretensión de
“cuantificación económica de los impactos del evitar la emisión de 77 MtCO2e en 2030, lo que
cambio climático”, “creación de capacidades en representa el 70% de su meta no condicionada
recursos humanos y la mejora en la coordinación total, así como de adaptar la infraestructura de
interinstitucional para la planificación y gestión, generación, transporte y distribución de energía
en temas de adaptación al cambio climático”, etc. para disminuir su vulnerabilidad ante los efectos
adversos del cambio climático.
El Gobierno de Argentina justifica su postura
explicando que el país sólo emite el 0.7% de las El Plan pretende actuar tanto sobre la oferta
emisiones mundiales de GEI, mientras que su de energía como sobre la demanda, e implica
contribución de mitigación revisada representa medidas de mitigación relativas a la eficiencia
una participación de 2.7% sobre el total de energética, al desarrollo de energías renovables,
las reducciones no condicionadas para 2030 la promoción de los biocombustibles y la gene-
comunicadas por las Partes de la CMNUCC, las ración a gran escala de energía a partir de centrales
cuales se elevan a 4,000 MtCO2e (UNEP, 2016). nucleares e hidroeléctricas, la sustitución de
petróleo por gas natural y el incremento de la
No obstante, es menester resaltar el esfuerzo del eficiencia de las centrales térmicas. Pretende
Gobierno de Argentina para desarrollar planes la disminución de los combustibles fósiles en
sectoriales para reducir emisiones de GEI y la matriz energética del país, los cuales serían
adaptarse al cambio climático. En particular, el sustituidos por energía procedente de las
Plan de Acción Nacional de Energía y Cambio renovables y de centrales nucleares. En materia
Climático (Gabinete Nacional de Cambio de generación eléctrica, pretende aumentar con-
Climático, Ministerio de Energía y Ministerio de siderablemente la proporción de las energías
Ambiente y Desarrollo Sustentable, 2017) merece renovables (de 2% a 25%), así como la del
ser destacado. nuclear (de 6% a 15%), a la par con una fuerte
disminución de la proporción de la energía tér-
mica (de 66% actualmente a 31%) en 2030.
Plan de Acción Nacional de Energía y Cambio
Climático Es de reconocer que dicho Plan está muy bien
elaborado y contiene medidas pertinentes y
La INDC presentada por Argentina ante la ambiciosas para reducir las emisiones de aquí
CMNUCC no fue muy explícita respecto a las a 2030. Sin embargo, la intensificación de la
medidas que se tomarán para disminuir las producción de gas natural, así como la explota-
emisiones en el sector energético del país. La ción de yacimientos no convencionales en Vaca
publicación del Plan de Acción Nacional de Muerta, ambos promovidos en la política ener-
Energía y Cambio Climático (Gabinete Nacional gética actual del país, se contraponen con el
de Cambio Climático, Ministerio de Energía y objetivo establecido en el presente Plan de
Ministerio de Ambiente y Desarrollo Sustentable, disminuir la importancia de los combustibles
140
ENERLAC • Volumen IV. Número 2. Diciembre, 2020. ISSN: 2602-8042 (impreso) / 2631-2522 (digital).
OLADE – AUGM
fósiles en la matriz energética del país y gene- se pueden controlar las emisiones fugitivas, a
rarán invariablemente un aumento de las mayor intensidad de la producción de gas natural,
emisiones de CH4 en el sector energético. existen mayores emisiones de CH4.
En efecto, varios estudios recientes revelan la Las emisiones de CH4 serían todavía más
importancia de las emisiones de CH4 a lo largo importantes en la explotación de yacimientos
de la cadena de valor del sector hidrocarburos. no convencionales. En un estudio de Howarth
Le Bail et al. (2019) explican que la contribución et al. (2011), en el cual estiman las emisiones de
del sector hidrocarburos se debe a que “el CH4 CH4 del ciclo de vida del gas natural, los autores
es uno de los principales componentes del gas llegan a la conclusión que estas emisiones
natural, mismo que es liberado a la atmósfera pueden ser especialmente elevadas en la ex-
-de manera deliberada o accidental- durante tracción de recursos no convencionales: estiman
la producción, procesamiento y transporte del que “entre 3.6% y 7.9% del metano obtenido
petróleo y gas natural” por lo que determinan que gracias a la producción de gas no convencional
“el papel que pueda tener el gas natural dentro escapa a la atmósfera debido al venteo y fugas
de las estrategias globales de descarbonización a lo largo del ciclo de la vida del pozo”, lo que
de la matriz energética depende en gran medida representaría “un aumento de 30% y quizás más
de la minimización de las emisiones de CH4 del doble de la producción de gas convencional”.
de la cadena de valor del sector hidrocarburos”. Explican que las mayores emisiones de CH4 en la
producción de gas shale se debe principalmente
al escape de las emisiones durante la perforación,
que es mucho más intensa que en la producción
convencional, y por los fluidos de retorno,
Alrededor de 1/4 de las también en volúmenes más importantes que en
emisiones totales serían la producción convencional. El estudio de Honty,
de CH4, las cuales G. (2014) coincide con esta postura, y agrega que
provendrían en un 40% las emisiones fugitivas pueden ocurrir durante
de emisiones fugitivas la fase de perforación, durante el fracturamiento
durante la producción hidráulico, durante la fase de producción y,
también, durante la fase de abandono si los
de gas natural.
pozos abandonados no están sellados de manera
adecuada.
141
CONTRIBUCIÓN DE LA CADENA PRODUCTIVA DE LA ENERGÍA EN LAS EMISIONES DE GASES DE EFECTO INVERNADERO EN ARGENTINA.
Le Bail, Maxime
142
ENERLAC • Volumen IV. Número 2. Diciembre, 2020. ISSN: 2602-8042 (impreso) / 2631-2522 (digital).
OLADE – AUGM
CONCLUSIÓN
En 2014, las emisiones de Argentina fueron de y extracción de los combustibles fósiles y, parti-
368 MtCO2e. Lo anterior representó cerca del cularmente, el abandono de la explotación de
0.7% de las emisiones totales a nivel mundial y recursos no convencionales con el fracturamiento
posicionó a Argentina como el vigésimo noveno hidráulico en Vaca Muerta, y su sustitución pro-
país que más GEI emite. gresiva por energías renovables y más limpias.
143
CONTRIBUCIÓN DE LA CADENA PRODUCTIVA DE LA ENERGÍA EN LAS EMISIONES DE GASES DE EFECTO INVERNADERO EN ARGENTINA.
Le Bail, Maxime
REFERENCIAS
Agencia Internacional de Energía. (2020). Key energy Joselow, M. (2020). IEA chief Fatih Birol on trauma
statistics, 2018. Total primary energy supply (TPES) by across the energy sector. E&E News. 22 de mayo de 2020.
source, Argentina 1990-2017. Recuperado el 25/07/2020 en: Recuperado el 13/06/2020 en: https://www.eenews.net/
https://www.iea.org/countries/argentina stories/1063205029?stream=top
Álvarez, J. M. (2017). Emisiones de gases de efecto invernadero Le Bail, M. Zavala-Araiza, D. y Gómez Solares, A. M. (2019).
en el sector de gas y petróleo en Argentina. Situación y La mitigación de las emisiones de metano en el sector
Perspectivas (Tesis de Maestría). Centro de Estudios de hidrocarburos: la medida ausente del Reporte del IPCC sobre
la Actividad Regulatoria Energética (CEARE). Argentina. las consecuencias de un aumento de la temperatura media
Recuperado el 25/07/2020 en: https://www.ceare.org/ global de 1.5°C. En: José Clemente Rueda Abad (Ed.) 2019.
tesis/2017/tes18.pdf ¿Aún estamos a tiempo para el 1.5°C? Voces y Visiones sobre
el Reporte Especial del IPCC. Universidad Nacional Autónoma
Bidegaray, M. (2020). Coronavirus en Argentina: cae 23% el De México - Programa de Investigación en Cambio Climático.
consumo de electricidad en la primera semana de cuarentena. ISBN: 978-607-30-2099-2
Clarín. 01/04/2020. Recuperado el 13/06/2020 en:
https://www.clarin.com/economia/economia/coronavirus- Le Quéré, C., Jackson, R.B., Jones, M.W. et al. (2020).
argentina-cae-23-consumo-electricidad-primera-semana- Temporary reduction in daily global CO2 emissions during
cuarentena_0_968Atlwxi.html the Covid-19 forced confinement. Nat. Clim. Chang. https://doi.
org/10.1038/s41558-020-0797-x
Convención Marco de Naciones Unidas sobre el Cambio
Climático (CMNUCC). (Sin fecha). INDC Registry - Argentina. Moreira Muzio, M., Gaioli, F., y Galbusera, S. (2019). Inventario
Recuperado el 08/06/2020 en: https://www4.unfccc.int/ nacional de gases de efecto invernadero: Argentina, 2019.
sites/INDCStaging/Pages/All.aspx Ciudad Autónoma de Buenos Aires: Secretaría de Ambiente y
Desarrollo Sustentable de la Nación.
Gabinete Nacional de Cambio Climático, Ministerio de Energía,
Ministerio de Ambiente y Desarrollo Sustentable (2017). Plan Ourworldindata.org. (2020). Argentina - Population, 1800
de Acción Nacional de Energía y Cambio Climático. Versión I - to 2100. Historical estimates of population combined with
2017. the projected population to 2100 based on the UN’s medium
variant scenario. Recuperado el 07/06/2020 en: https://
Garrison, C. (2020). Refinería de argentina YPF reduce ourworldindata.org/grapher/projected-population-by-
personal por caída del consumo. Infobae. 29 de abril de 2020. country?country=~ARG
Recuperado el 13/06/2020 en: https://www.infobae.com/
america/agencias/2020/04/29/refineria-de-argentina-ypf- Ourworldindata.org. (2020 bis). Argentina -Gross Domestic
reduce-personal-por-caida-del-consumo/ Product, 1960 to 2017. Gross domestic product adjusted for
price changes over time (inflation) and expressed in US-Dollars.
Honty, G. (2014). Potenciales emisiones de la explotación de Recuperado el 07/06/2020 en: https://ourworldindata.org/
hidrocarburos no convencionales. Centro Latino Americano de grapher/gross-domestic-product?tab=chart&country=~ARG
Ecología Social. Informe preparado especialmente para 350.
org. Octubre, 2014. República Argentina. (2015). Contribución Prevista y
Determinada a Nivel Nacional. 01 de octubre de 2015. Recu-
Howarth, R. W., Santoro, R. and Ingraffea, A. (2011). Methane perado el 08/06/2020 en: https://www4.unfccc.int/sites/
and the greenhouse-gas footprint of natural gas from shale submissions/INDC/Published%20Documents/Argentina/1/
formations. Climatic change, 106(4), 679. INDC%20Argentina.pdf
144
ENERLAC • Volumen IV. Número 2. Diciembre, 2020. ISSN: 2602-8042 (impreso) / 2631-2522 (digital).
OLADE – AUGM
145
HACIA UNA JUSTICIA ENERGÉTICA EN MÉXICO.
UN ANÁLISIS DEL RÉGIMEN TARIFARIO DEL
SECTOR ELÉCTRICO COLOMBIANO Y SU POSIBLE
IMPLEMENTACIÓN EN LA REPÚBLICA MEXICANA
TOWARDS AN ENERGY JUSTICE IN MEXICO. AN ANALYSIS OF THE TARIFF REGIME
OF THE COLOMBIAN ELECTRICITY SECTOR AND ITS POSSIBLE IMPLEMENTATION
IN THE MEXICAN REPUBLIC
146
ENERLAC • Volumen IV. Número 2. Diciembre, 2020. ISSN: 2602-8042 (impreso) / 2631-2522 (digital).
OLADE – AUGM
RESUMEN ABSTRACT
El estudio tiene como objetivo analizar y The study aim is to compare and analyse the
comparar el régimen tarifario de la energía Colombian and Mexican electricity tariff regime.
eléctrica en Colombia y México. El hecho de que The current scenario in which ninety-nine
en México exista una cobertura del (99%) de percent (99%) of the population has access to
la población da pauta a buscar esquemas tarifa- electricity is the perfect scenario to review the
rios en el que los recursos públicos destinados electricity tariff and the subsidy schemes to
al subsidio universal sean destinados a otros redefine the public resources provided for that
rubros con el objetivo de dar un suministro de purpose to improve a more reliable, continuous,
energía más confiable, continuo y sostenible en and sustainable energy supply to achieve an
aras de lograr una mayor justicia energética. energy justice.
El esquema de subsidios en los servicios públicos The subsidy scheme for domestic public services
domésticos en Colombia es un modelo qué de in Colombia is a model that, if it is similarly im-
implementarlo de forma similar, podría brindar plemented in Mexico, it could provide a solution to
una solución a los problemas actuales en México. the current lack of public resources for electricity
infrastructure in Mexico.
Palabras clave: Energía Eléctrica, Tarifas Keywords: Electric Power, Electric Tariff, Subsidy,
Eléctricas, Subsidio, Justicia Energética, México, Energy Justice, Mexico, Colombia.
Colombia.
147
HACIA UNA JUSTICIA ENERGÉTICA EN MÉXICO. UN ANÁLISIS DEL RÉGIMEN TARIFARIO DEL SECTOR ELÉCTRICO COLOMBIANO Y SU POSIBLE IMPLEMENTACIÓN
EN LA REPÚBLICA MEXICANA. Beltrán Nishizaki, Ernesto
Justicia energética
Ante el reto de mayor inversión en el sector
eléctrico mexicano como detonante para No existe una definición única aceptada por
el desarrollo y la falta de recursos públicos la academia sobre qué es justicia energética,
para lograrlo, existe una imperiosa necesidad sin embargo, este concepto se desarrolló con
de buscar nuevos esquemas que permitan tres objetivos fundamentales que son el acceso
refocalizar los recursos destinados al subsidio universal a una energía segura, sostenible y un
de energía eléctrica a través de un uso eficiente precio asequible (Sourav, 2018).
de los recursos sin dejar a un lado a las clases
menos favorecidas. Para ello, la doctrina ha señalado 8 elementos
que deben ser considerados para conseguirlo.
Este artículo analizará, a través de una meto- Estos son: disponibilidad, asequibilidad, debido
dología hipotética-deductiva, el modelo tarifario proceso, transparencia y rendición de cuentas,
eléctrico colombiano y su posible implementación sostenibilidad, equidad intergeneracional,
en México con el objetivo de reducir el gasto equidad social y, finalmente, responsabilidad
público destinado al subsidio universal de la tarifa (Sovacool y Dworkin, 2015). Entendidos estos
eléctrica en casa habitación, para que de esta cómo:
forma se focalicen los recursos públicos en otros
proyectos prioritarios que el Estado considere 1. Disponibilidad. El garantizar un suministro
estratégicos para el desarrollo de la nación y un seguro, confiable y suficiente.
sistema energético más justo.
2. Asequibilidad. El costo que permita el acceso
En la primera parte del artículo analizaremos universal a la energía eléctrica. Esto no quiere
el esquema de subsidio (directo y cruzado) en decir que para ello se requiera únicamente
el servicio público de energía eléctrica, imple- precios bajos sino también que la facturación no
mentado en Colombia. Al mismo tiempo, se sea excesiva (Sovacool y Dworkin, 2015).
contrastarán las ventajas y desventajas del
mecanismo utilizado para determinar la apli- 3. Debido proceso. Actualmente la regulación
cabilidad del subsidio, las reglas bajo el cual busca cumplir con esta condición ya que la
opera dicho esquema y el posible cumplimento de regulación obliga a presentar una licencia de
los elementos que conforman el principio de impacto social, consulta indígena y manifestación
justicia energética. de impacto ambiental para los proyectos
eléctricos1. Con ello se busca que los interesados
Posteriormente, se realizará un análisis jurídico puedan conocer las características del proyecto y,
del sistema eléctrico mexicano y el subsidio en su caso, oponerse a través de los mecanismos
universal existente en la tarifa eléctrica resi- previstos en los ordenamientos jurídicos.
dencial. El costo de mantener un sistema de
subsidio que beneficia -en mayor medida- a los 4. Transparencia y rendición de cuentas.
que más consumen, no incentivan el uso eficiente La transparencia de la información contable y
(Sánchez et al., 2018) y sostenible de los recursos. financiera que permita conocer el costo real que
Por último, se analiza si el cambio de mecanismo 1 Con excepción de aquellos que se encuentren exento de
para el otorgamiento del subsidio en México permiso por parte de la CRE. Como plantas de emergencia y
a través de la implementación de un régimen generación por debajo del .5MW.
148
ENERLAC • Volumen IV. Número 2. Diciembre, 2020. ISSN: 2602-8042 (impreso) / 2631-2522 (digital).
OLADE – AUGM
COLOMBIA
149
HACIA UNA JUSTICIA ENERGÉTICA EN MÉXICO. UN ANÁLISIS DEL RÉGIMEN TARIFARIO DEL SECTOR ELÉCTRICO COLOMBIANO Y SU POSIBLE IMPLEMENTACIÓN
EN LA REPÚBLICA MEXICANA. Beltrán Nishizaki, Ernesto
departamentos, los distritos, los municipios y las garantizar la calidad, eficiencia, libre competen-
entidades descentralizadas conceder subsidios, cia, acceso universal y solidaridad (siendo este
en sus respectivos presupuestos, para lograr una último la ratio legis del ordenamiento jurídico).
mayor equidad social y las personas de menores
ingresos puedan pagar las tarifas de suministro Para su consecución, contempla el otorgamiento
eléctrico que cubran sus necesidades básicas. de subsidios para las personas de menores
ingresos y en la obligación de realizar una
Incluso, a través de la sentencia número C-580/92 contraprestación extraordinaria a los estratos
de su Corte Constitucional, la CPRC reconoce el más altos.
principio fundamental de “solidaridad” aplicable
al régimen tarifario, el cual considera que los La energía eléctrica como servicio público
servicios púbicos se prestarán en beneficio de
la comunidad y elimina la posibilidad que se Mediante la Ley 143 de 1994, se establece el
establezca la prestación de estos en condiciones régimen para la generación, interconexión,
favorables a un segmento de la sociedad. Al transmisión, distribución y comercialización
mismo tiempo, este principio impone la obli- de electricidad en el territorio nacional. La
gación de tomar en cuenta los distintos estratos ley prevé principios que regirán el servicio de
sociales que participan como usuarios de los energía eléctrica. Estos son eficiencia, calidad,
servicios públicos, su capacidad económica y continuidad, adaptabilidad, neutralidad, solida-
como consecuencia el establecimiento de tarifas ridad y equidad, estipulados en el artículo 6
diferenciales para impedir que la mayoría de la de la Ley 143 de 1994. Para efecto del artículo,
población asuma el mayor costo de los servicios. únicamente haremos un análisis de los últi-
mos dos.
Aunado a lo anterior, el artículo 367 de la CPRC
determina que en todo momento el esquema El principio de solidaridad, bajo el contexto
tarifario debe buscar cumplir con otro de los de la ley 143 de 1994, se debe entender como
objetivos planteados: solucionar las necesidades el mandato de diseñar un régimen tarifario
básicas insatisfechas y el acceso a los servicios que tome en consideración factores para que
públicos domiciliados de las clases menos los sectores de consumo de mayores ingresos
favorecidas. En otras palabras, que exista ayuden a que las personas de menos ingresos
disponibilidad, asequibilidad y equidad social. puedan pagar las tarifas de consumo eléctrico
que cubra las necesidades básicas y con ello
Con la finalidad de cumplir con el mandato lograr cumplir con la asequibilidad del servicio.
constitucional, las autoridades colombianas En cambio, el principio de equidad se refiere al
crearon la normatividad para garantizar el objetivo de alcanzar una cobertura equilibrada
equilibrio económico y financiero de las empresas, y adecuada en todo el país que garantice la
así como la cobertura futura de los servicios y el satisfacción de las necesidades básicas (Casas,
acceso universal al servicio de electricidad. Medina y Meléndez, 2005).
150
ENERLAC • Volumen IV. Número 2. Diciembre, 2020. ISSN: 2602-8042 (impreso) / 2631-2522 (digital).
OLADE – AUGM
151
HACIA UNA JUSTICIA ENERGÉTICA EN MÉXICO. UN ANÁLISIS DEL RÉGIMEN TARIFARIO DEL SECTOR ELÉCTRICO COLOMBIANO Y SU POSIBLE IMPLEMENTACIÓN
EN LA REPÚBLICA MEXICANA. Beltrán Nishizaki, Ernesto
152
ENERLAC • Volumen IV. Número 2. Diciembre, 2020. ISSN: 2602-8042 (impreso) / 2631-2522 (digital).
OLADE – AUGM
2019). El contar con una metodología para el a) La aprobación de la tarifa por la Comisión
cobro de las contraprestaciones cumple en todo Reguladora de Energía (CRE), órgano regulador
momento con el principio de transparencia y coordinado con autonomía técnica operativa y de
rendición de cuentas consagrado como otro gestión.
elemento fundamental para lograr una justicia
energética. b) Existen 12 tarifas4 regionalizadas5.
153
HACIA UNA JUSTICIA ENERGÉTICA EN MÉXICO. UN ANÁLISIS DEL RÉGIMEN TARIFARIO DEL SECTOR ELÉCTRICO COLOMBIANO Y SU POSIBLE IMPLEMENTACIÓN
EN LA REPÚBLICA MEXICANA. Beltrán Nishizaki, Ernesto
154
ENERLAC • Volumen IV. Número 2. Diciembre, 2020. ISSN: 2602-8042 (impreso) / 2631-2522 (digital).
OLADE – AUGM
(Sánchez et al., 2018). De igual forma, prueba A diferencia de Colombia en donde el pago de
que beneficia a las clases con mayores recursos y la cuota extraordinaria depende del estrato
no a las menos favorecidas, ya que las primeras social al que se forme parte, en México se busca
de ellas consumen mayor cantidad de energía que sea a través del monto de consumo men-
eléctrica. sual total, ya que no existe una regulación que
permita una estratificación uniforme, no refleja
Un cambio de régimen de subsidio por un -en general- la capacidad de pago de los
principio de justicia energética usuarios finales y no incentiva un consumo
racional de los recursos.
Como se ha mencionado, el régimen de subsidio
a la tarifa eléctrica en México se rige bajo el En ese sentido, se propone que la estructura para
principio de universalidad el cual tiene como el cobro de la cuota extraordinaria sea similar
objetivo incentivar el acceso y uso de la energía a la actual, es decir, que la tarifa residencial
eléctrica a la población en general; sin embargo, autorizada para el suministro de energía eléctrica
en el escenario actual en el que la cobertura en casa habitación sea segmentada conforme a
eléctrica es de 99% (INEGI, 2018), los efectos de los límites de consumo dependiendo la región
mantener una política de subsidio similar generan en la que se encuentre la localidad y tomando en
un sistema menos justo. Por un lado, un uso consideración un consumo mínimo para cubrir
indiscriminado e irresponsable de los recursos y, las necesidades básicas. Todo consumo mayor
por otro, una infraestructura cada vez más vieja deberá ser tomado como base para el pago de la
que pone entredicho la prestación de un servicio contraprestación extraordinaria.
continuo, eficiente y sostenible.
155
HACIA UNA JUSTICIA ENERGÉTICA EN MÉXICO. UN ANÁLISIS DEL RÉGIMEN TARIFARIO DEL SECTOR ELÉCTRICO COLOMBIANO Y SU POSIBLE IMPLEMENTACIÓN
EN LA REPÚBLICA MEXICANA. Beltrán Nishizaki, Ernesto
II. Asequibilidad. Con un subsidio focalizado VI. Equidad intergeneracional. Para efecto de
se lograría el apoyo real y directo al sector este documento se entiende bajo lo estipulado por
con menores ingresos y evitar el pago de una Fernández Santos (2005) como aquella que hace
contraprestación desproporcional o destinar referencia a la igualdad en el ámbito ambiental
gran parte de su presupuesto para ello. entre distintas generaciones. En otras palabras, es
el derecho que tienen las generaciones futuras de
III. Debido proceso. Actualmente la regulación desarrollarse o vivir en un lugar al menos con las
cumple con esta condición ya que la regulación mismas características ambientales a las actuales.
obliga a presentar una licencia de impacto social,
consulta indígena y manifestación de impacto VII. Equidad social. Se debe entender como un
ambiental para los proyectos eléctricos9. Con ello acceso imparcial a la energía eléctrica. El acceso
se busca que los interesados puedan conocer universal a la energía eléctrica es un derecho
las características del proyecto y, en su caso, humano, ya que sin ese insumo básico no se pue-
oponerse a través de los mecanismos previstos de cubrir las necesidades básicas. El focalizar
en los ordenamientos jurídicos. los subsidios a grupos vulnerables permitiría
redireccionar recursos públicos para mejor ac-
IV. Transparencia y rendición de cuentas. ceso a una energía segura y confiable, así como
El cambio en el mecanismo de subsidio debe programas sociales para mejorar la calidad de
aparejar mayor transparencia en los recursos vida (como podría ser el uso de estufas eléctricas
públicos y la aportación extraordinaria realizada que permitan reducir el consumo de leña o
por aquellos consumidores con mayor consumo carbón para calentar o cocinar alimentos). De
y capacidad de pago que sean destinados para acuerdo con cifras oficiales, el 11% de la población
apoyar al sector población con menos ingresos. sigue utilizando como principal combustible la
Actualmente se conoce el monto total destinado leña (INEGI – ENIGH, 2018).
a subsidiar el suministro eléctrico, pero no hay
claridad ni información de fácil acceso para VIII. Responsabilidad. Es la corresponsabilidad
conocer los montos efectivamente erogados. que tiene el gobierno y ciudadanía de proteger el
medio ambiente y reducir los impactos negativos
V. Sostenibilidad. La explotación de los recursos derivado del suministro de energía eléctrica.
debe tener como limitante el no comprometer Al igual que en los otros casos, el focalizar los
su disponibilidad para las futuras generaciones. subsidios permitirá destinar recursos públicos
La aportación extraordinaria propuesta sería para el desarrollo de infraestructura eléctrica
un aliciente para reducir el consumo excesivo que mitigue el impacto ambiental como líneas
de los recursos. El hecho de no recibir apoyo de distribución y transmisión que han registrado
gubernamental para el consumo aunado a una hasta un 21.2% de pérdidas de energía en el país
aportación extraordinaria lograría desincentivar (CEPAL, 2018), sustitución de equipos obsoletos
prácticas que generen mayor consumo. para reducir la pérdida de energía eléctrica.
156
ENERLAC • Volumen IV. Número 2. Diciembre, 2020. ISSN: 2602-8042 (impreso) / 2631-2522 (digital).
OLADE – AUGM
REFERENCIAS
Beato, P. (1999). Cross subsidies in public services: Some CONUEE-SENER (2020). Clasificación de climas y su
issues. Estados Unidos: Inter- American Development Bank. aplicación a la norma para envolvente de viviendas: notas
para discusión. México: Cuadernos de la CONUEE. No. 6/
Casas L., C., Medina F., P., Meléndez A., M. (2005). Nuevo Ciclo. Disponible en https://www.conuee.gob.
Subsidios al consumo de los servicios públicos en Colombia: mx/transparencia/boletines/Cuadernos/cuaderno6Nvo
¿hacia dónde nos movemos? Disponible en: https:// Ciclo.pdf
www.repository.fedesarrollo.org.co/bitstream/handle
/11445/1085/Co_So_Diciembre_2005_Melendez.pdf; Corte Constitucional de la República de Colombia.
jsessionid=E0DCA0908A0F963954E0548CC0AF7E5E Sentencia número C. -580/92. Colombia. 5 de noviembre
?sequence=2 de 1992. Disponible en: https://www.corteconstitucional.
gov.co/relatoria/1992/C-580-92.htm
CEPAL (2018). Informe nacional de monitoreo de eficiencia
energética de México 2018. México: Naciones Unidas. Fernández Santos, Y. (2005). Evaluación de la equidad
Disponible en: https://repositorio.cepal.org/bitstream/ intrageneracional e intergeneracional en una entidad
handle/11362/43612/1/S1800496_es.pdf pública local a través de indicadores de gestión. IX
Congresso Internacional de Custos – Florianópolis, SC.
Constitución Política de la República de Colombia. Artículos Brasil, 28 a 30 de noviembre de 2005. Disponible en:
365, 366 y 367. Colombia. Disponible en: http://www. https://anaiscbc.emnuvens.com.br/anais/article/view
secretariasenado.gov.co/senado/basedoc/constitucion_ File/1944/1944
politica_1991.html
157
HACIA UNA JUSTICIA ENERGÉTICA EN MÉXICO. UN ANÁLISIS DEL RÉGIMEN TARIFARIO DEL SECTOR ELÉCTRICO COLOMBIANO Y SU POSIBLE IMPLEMENTACIÓN
EN LA REPÚBLICA MEXICANA. Beltrán Nishizaki, Ernesto
Gordillo, A. (2013). Tratado de derecho administrativo y Ministerio de la Información y Tecnologías (2013 al 2016).
obras selectas. Teoría general del derecho administrativo, Cobertura del servicio de energía eléctrica por municipio.
tomo 8. Buenos Aires: FDA. Disponible en: http://gordillo. Colombia. Disponible en: https://www.datos.gov.co/Minas-
com/tomo8.php y-Energ-a/Cobertura-del-servicio-de-energia-electrica-por-
mu/fuk9-m94c/data
INEGI (2018). Primera encuesta nacional sobre consumo
de energéticos en viviendas particulares (encevi). 7 de PRODESEN - Centro Nacional de Control de Energía (2019).
noviembre de 2018. Disponible en: https://www.inegi. Programa de Ampliación y Modernización de la Red Nacional
org.mx/contenidos/saladeprensa/boletines/2018/ de Transmisión y Redes Generales de Distribución del
EstSociodemo/ENCEVI2018.pdf Mercado Eléctrico Mayorista PRODESEN 2019-2033, Centro
Nacional de Control de Energía. Disponible en: https://
INEGI – ENIGH (2018). Encuesta Nacional de Ingresos www.cenace.gob.mx/Docs/Planeacion/ProgramaRNT/
y Gastos de los Hogares. México: ENIGH. Disponible en: Programa%20de%20Ampliaci%C3%B3n%20y%20
https://www.inegi.org.mx/contenidos/programas/enigh/ Modernizaci%C3%B3n%20de%20la%20RNT%20y%20
nc/2018/doc/enigh2018_ns_presentacion_resultados.pdf RGD%202019%20-%202033.pdf
Ley 142 de 1994. Artículo 2. Por la cual se establece el Rojas López, J.G. (2017). El régimen de contratación de las
régimen de los servicios públicos domiciliarios y se dictan empresas de servicios públicos, domiciliarios de carácter
otras disposiciones. 11 de julio de 1994. Colombia: Diario oficial. Ratio Juris, Vol. 1 Núm. 2 (2005): Enero – Junio, 23-
Oficial No. 41.433, de 12 de julio de 1994. Disponible en: 24. Disponible en: https://publicaciones.unaula.edu.co/
http://www.bogotajuridica.gov.co/sisjur/normas/Norma1. index.php/ratiojuris/article/view/294
jsp?i=2752
Sánchez, L., Echeverría, D., Wooders, P., Kuehne, K., Lean,
Ley 142 de 1994. Artículo 3, Inciso 3.9. Colombia: Diario T., Beaton, C., Shama, S., Oharenko, Y. (2018). Mejorando
Oficial No. 41.433, de 12 de julio de 1994. Por la cual se y refocalizando los subsidios a la electricidad. Opciones
establece el régimen de los servicios públicos domiciliarios para su optimización. México - Alemania: International
y se dictan otras disposiciones. 11 de julio de 1994. Institute for Sustainable Development (ISSD) Disponible
Disponible en: http://www.bogotajuridica.gov.co/sisjur/ en: https://www.energypartnership.mx/fileadmin/user_
normas/Norma1.jsp?i=2752 upload/mexico/media_elements/reports/Subsidios
Electricos-MEX.pdf
Ley 143 de 1994. Por la cual se establece el régimen para
la generación, interconexión, transmisión, distribución y Santa María, M., Von Der Fehr, N.H., Millán, J., Benavides,
comercialización de electricidad en el territorio nacional, se J., Gracia, O., Schutt, E. (2009). El mercado de la energía
conceden unas autorizaciones y se dictan otras disposiciones eléctrica en Colombia: características, evolución e impacto
en materia energética. 12 de julio de 1994. Colombia: sobre otros sectores. Cuaderno de Fedesarrollo No. 30.
Diario Oficial No. 41.434. Disponible en: http://www. Colombia: Fedesarrollo.
secretariasenado.gov.co/senado/basedoc/ley_0143_1994.
html Senado de la República (2013). Proyecto de Dictamen en
las Comisiones Unidas de Energía y de Estudios Legislativos,
Ley 143 de 1994. Por la cual se establece el régimen para primera con proyecto de Decreto por el que se expiden la
la generación, interconexión, transmisión, distribución y Ley de la Industria Eléctrica, Ley de Energía Geotérmica y se
comercialización de electricidad en el territorio nacional, adiciona y reforman diversas disposiciones de la Ley de Aguas
se conceden unas autorizaciones y se dictan otras Nacionales. México. Disponible en: https://www.senado.
disposiciones en materia energética. 12 de julio de 1994. gob.mx/comisiones/energia/docs/reforma_energetica/
Colombia: Diario Oficial No. 41.434. Disponible en: dictamen2.pdf
http://www.secretariasenado.gov.co/senado/basedoc/
ley_0143_1994.html
158
ENERLAC • Volumen IV. Número 2. Diciembre, 2020. ISSN: 2602-8042 (impreso) / 2631-2522 (digital).
OLADE – AUGM
159
Asociación de Universidades
GRUPO MONTEVIDEO