Propuesta de ubicaci´onde la infraestructura de recarga r´apidapara veh´ıculosel´ectricosen

Universidad de Costa Rica Escuela de Ingenier´ıaEl´ectrica Electric Power and Energy Research Laboratory

Autores:

Ing. Jairo Quir´osTort´os,PhD. Ing. Luis Fernando Victor Gallardo Ing. Juan Mart´ınez Barboza Ing. Ricardo Bejarano Villachica Ing. Luis Ruiz Aguilar Ing. Jam Angulo Paniagua Sr. Daniel Fuentes Soto

Setiembre del 2018 ´Indice

Acr´onimosy siglas I

1. Introducci´on 1 1.1. Criterios para ubicar centros de recarga r´apida...... 3 1.2. Ubicaciones generales sugeridas durante el taller...... 4

2. Metodolog´ıade ubicaci´onde estaciones de recarga en Costa Rica8 2.1. Ubicaci´onde los puntos de recarga en la GAM...... 9 2.2. Ubicaci´onde los puntos de recarga fuera de la GAM...... 16

3. Verificaci´onde distancias y an´alisisde resultados 24 3.1. Rutas de la GAM a Paso Canoas...... 24 3.2. Rutas a ...... 27 3.3. Rutas a Guanacaste...... 30 3.4. Rutas hacia norte del pa´ıs...... 34

4. Recomendaciones finales para Ubicaci´on 38

5. Conclusiones 41

Referencias 42

i Propuesta de ubicaci´onde infraestructura de recarga r´apida EPERLab - UCR

´Indice de figuras

1. Cantidad de personas mayores de 15 a˜nosque trabajan fuera del cant´onen que habitan...... 10 2. Cantidad de veh´ıculosfamiliares y personas mayores de 15 a˜nosque trabajan fuera de la zona que habitan...... 10 3. Zonas de destino, aglomeraciones de trabajo y estudio...... 11 4. Rutas para la zona sur de la GAM...... 12 5. Rutas para la zona este de la GAM...... 12 6. Rutas para la zona noreste de la GAM...... 13 7. Rutas para la zona norte de la GAM...... 14 8. Rutas para la zona centro de la GAM...... 14 9. Rutas para la zona noroestede la GAM...... 15 10. Rutas para la zona suroeste de la GAM...... 15 11. Diagrama general del VE...... 16 12. Perfiles de manejo UDDS y HWFET de EPA...... 17 13. Consumo de la ruta: Limonal – Barranca – Orotina - San Jos´e...... 20 14. Consumo de la ruta: San Jos´e– Orotina – Barranca -Limonal...... 20 15. Consumo de la ruta: San Jos´e– Orotina – Barranca -Limonal...... 21 16. Consumo de la ruta: Limonal – Barranca – Orotina - San Jos´esin Orotina sin Orotina...... 21 17. Consumo de la ruta: San Jos´e– Orotina – Barranca -Limonal sin Barranca..... 22 18. Consumo de la ruta: Limonal – Barranca – Orotina - San Jos´esin Barranca.... 22 19. Rutas desde San Jos´ehasta ciudad Neily...... 25 20. C´alculosde la ruta 1 San Jos´e- Paso Canoas escenario 1...... 26 21. C´alculosde la ruta 1 San Jos´e- Paso Canoas escenario 2...... 26 22. C´alculosde la ruta 2 San Jos´e- Paso Canoas escenario 1...... 27 23. Rutas desde San Jos´ehasta Sixaola...... 28 24. C´alculosde la ruta 1 San Jos´e- Sixaola escenario 1...... 28 25. C´alculosde la ruta 1 San Jos´e- Sixaola escenario 2...... 29 26. C´alculosde la ruta 2 San Jos´e- Sixaola escenario 1...... 29 27. C´alculosde la ruta San Jos´e- Sixaola escenario 2...... 30 28. Rutas desde San Jos´ea los distintos puntos en Guanacaste...... 31 29. C´alculosde la ruta 1 San Jos´ea Limonal escenario 1...... 32 30. C´alculosde la ruta 1 San Jos´ea Limonal escenario 2...... 32

2018 ii Propuesta de ubicaci´onde infraestructura de recarga r´apida EPERLab - UCR

31. C´alculosde la ruta 2 San Jos´ea Limonal escenario 1...... 33 32. C´alculosde la ruta 2 San Jos´ea Limonal escenario 2...... 33 33. C´alculosde la ruta 1 Limonal a Pe˜nasBlancas escenario 1...... 33 34. C´alculosde la ruta 1 Limonal a Mal Pa´ısescenario 1...... 34 35. C´alculosde la ruta 1 Limonal a Tamarindo escenario 1...... 34 36. C´alculosde la ruta 2 Limonal a Tamarindo escenario 1...... 34 37. C´alculosde la ruta 1 San Jos´ea Los Chiles escenario 1...... 35 38. Rutas San Jos´ehacia el norte del pa´ıs...... 36 39. C´alculosde la ruta 2 San Jos´ea Los Chiles escenario 1...... 37 40. C´alculosde la ruta 2 San Jos´ea Los Chiles escenario 2...... 37 41. C´alculosde la ruta 3 San Jos´ea Los Chiles escenario 1...... 37

2018 iii Propuesta de ubicaci´onde infraestructura de recarga r´apida EPERLab - UCR

´Indice de tablas

1. Ubicaciones de inter´esmencionados en la mesa 1...... 4 2. Ubicaciones de inter´esmencionados en la mesa 2...... 5 3. Ubicaciones de inter´esmencionados en la mesa 3...... 6 4. Puntos de inter´esde todas las mesas de trabajo...... 7 5. Fuerzas involucradas en movimiento lineal del VE...... 16 6. Par´ametros que afectan al VE...... 16 7. Consideraciones sobre la eficiencia del sistema electro-mec´anicodel VE...... 18 8. Par´ametros utilizados para modelo de consumo del VE...... 19 9. Error entre rendimiento de fabricante y calculado con modelo...... 19 10. Capacidades nominales de bater´ıasde VEs y precios base de los modelos...... 23 11. Codificaci´onde colores para an´alisisde escenarios de ubicaci´on...... 24 12. Ubicaciones de la GAM - CNFL...... 38 13. Ubicaciones de la red de recarga m´ınima...... 39 14. Ubicaciones de la red de recarga media: a˜nadira la red m´ınima...... 40 15. Ubicaciones de la red de recarga robusta: a˜nadira la red media...... 40

2018 iv Propuesta de ubicaci´onde infraestructura de recarga r´apida EPERLab - UCR Acr´onimosy siglas

ARESEP Autoridad Reguladora de los Servicios P´ublicos.

ASOMOVE Asociaci´onCostarricense de Movilidad El´ectrica.

CMNUCC Convenci´onMarco de Naciones Unidas sobre el Cambio Clim´atico

EIE Escuela de Ingenier´ıaEl´ectrica.

EPERLab Laboratorio de Investigaci´onen Potencia y Energ´ıa.

GEI Gases de Efecto Invernadero.

GIZ Cooperaci´onT´ecnicaAlemana para el Desarrollo.

MIDEPLAN Ministerio de Planificaci´ony Pol´ıticaEcon´omica

MINAE Ministerio de Ambiente y Energ´ıa.

MOPT Ministerio de Obras P´ublicasy Transportes.

PND Plan Nacional de Desarrollo.

PNE Plan Nacional de Energ´ıa.

PNTE Plan Nacional de Transporte El´ectrico.

SEPSE Secretar´ıade Planificaci´ondel Subsector Energ´ıa.

UCR Universidad de Costa Rica.

VE Veh´ıculoEl´ectrico

2018 v 1. Introducci´on

La ubicaci´onde centros de recarga para VEs debe garantizar el desplazamiento de la tecno- log´ıapor el pa´ıs,evitando que este se quede sin energ´ıaen su bater´ıa.En particular, los centros de recarga r´apidabuscan brindar este servicio al menor tiempo de espera posible. El tema de ubicaci´ongeogr´aficafue un tema discutido en los talleres de infraestructura de recarga el´ectrica organizado por la Secretar´ıade Planificaci´ondel Subsector Energ´ıa(SEPSE) y la Cooperaci´on Alemana para el Desarrollo (GIZ). Algunas conclusiones importantes de los talleres son por ejem- plo, la definici´onen la Ley No 9518 de los 80 km en red vial nacional y 120 km en red cantonal para las distancias m´ınimas entre estaciones de recarga. A su vez, los participantes sugirieron diversos puntos de ubicaci´onde estaciones de recarga r´apida.

En el presente reporte se realiza un an´alisisde las ubicaciones de los centros de recarga r´apida considerando las recomendaciones de los participantes. Se considera adem´asque el ICE est´apor adquirir un total de 32 estaciones de recarga r´apidapara ubicar en el territorio nacional. Seg´unlas discusiones realizadas con los participantes, de esas 32, 4 se le entregar´ana la CNFL. Asimismo, el an´alisisconsidera que la CNFL tiene 5 estaciones adicionales.

A partir de estos antecedentes, el an´alisisse bas´oen las siguientes indagaciones:

⇒ ¿Puede un VE llegar a los puntos m´aslejanos del pa´ıs(respecto a la GAM) desde la GAM? ⇒ ¿Puede un VE devolverse de los mismos puntos (del apartado anterior) hasta la GAM? ⇒ ¿Cu´alesson las posibles ´areasde ubicaci´onseg´unlos criterios t´ecnicos?

Tal y como se mencion´o,las ubicaciones iniciales fueron sugeridas por los expertos de las empresas el´ectricasy dem´asparticipantes del taller. Superponiendo las sugerencias se obtiene un conjunto de 46 ubicaciones. La metodolog´ıapara encontrar la cantidad m´ınimade centros de recarga consiste en una eliminaci´onde las ubicaciones que menos benefician la extensi´onde la autonom´ıade un VE. La informaci´ongeogr´aficade la red vial nacional utilizada corresponde a la que el MOPT tiene disponible en l´ınea.Por otro lado, mediante los motores geogr´aficosGoogle Earth, Open Street Maps y Quantum GIS (QGIS), se realiza: i) una revisi´onde las distancias entre estaciones y ii) se extrae informaci´onde la topograf´ıa(elevaciones y asensos) para efectuar un an´alisisde autonom´ıaposteriormente. El an´alisispuede sintetizarse en los siguientes pasos:

1 Propuesta de ubicaci´onde infraestructura de recarga r´apida EPERLab - UCR

1 Cargar la red vial nacional en el software de manejo de informaci´ongeogr´aficaQGIS. 2 Utilizar los archivos de la red de Media Tensi´on(MT) dados por las empresas el´ectricas. 3 Ubicar los lugares iniciales (Tabla4) en el motor geogr´aficoOpen Street Maps. Descargar los archivos de este motor geogr´afico y cargarlos en QGIS. 4 Transponer las capas de los archivos de los pasos 1 y 2. Marcar en QGIS las coincidencias (´areas)entre la red vial y la red MT, cercanas a las zonas establecidas en el paso 3. 5 Ubicar las coincidencias en el motor geogr´aficoGoogle Earth. Marcar ´areasdel paso 4. 6 Calcular con Google Earth las rutas entre puntos lejanos. 7 Se listan las estaciones de cada ruta y se verifica el cumplimiento de las distancias. 8 Se extrae de Google Earth el perfil de elevaci´onentre una estaci´ony otra. 9 Utilizar modelo de consumo del VE (formalmente detallado en la seccion 2.2) para analizar el efecto de la topograf´ıaen su autonom´ıa.

La priorizaci´onse realiza en dos etapas: la primera consiste en utilizar las rutas de los puntos lejanos de la GAM para calcular la cantidad m´ınimade estaciones requeridos para llegar hasta o volver desde ellos. La segunda consiste en revisar las ubicaciones sugeridas para la GAM y verificar que estas atiendan la movilidad vehicular entre cantones de origen y destino de viajes: se utiliza informaci´ondel Censo 2011 y del Estado de la Naci´onpara identificar dichos or´ıgenes y destinos. La red de recarga inicial en la GAM ser´arobusta, debido a que existen 9 estaciones para CNFL. El ejercicio de simulaci´onde ubicaciones permite prever el efecto de las decisiones que las empresas distribuidoras pueden tomar en torno a la red de recarga r´apidade VEs en el corto plazo. La secci´on4 presenta una recomendaci´onfinal resultante del an´alisiscompleto.

2018 2 Propuesta de ubicaci´onde infraestructura de recarga r´apida EPERLab - UCR

1.1. Criterios para ubicar centros de recarga r´apida

Las mesas de trabajo determinaron cu´alesson los aspectos t´ecnicosque deben considerarse para escoger la ubicaci´onde los centros de recarga:

⇒ Considerar accidentes topogr´aficosque afectan autonom´ıade los EVs al elegir ubicaciones. ⇒ Tomar en cuenta el flujo vehicular para ubicar los centros de recarga. ⇒ Priorizar ubicaciones de los centros de recarga, cercanas a centros comerciales, tur´ısticos, lugares de descanso o de esparcimiento en general. Este criterio aplica para todas los centros de tipo comercial y de acceso p´ublico.Es posible que las alianzas entre empresas distribuidoras y otras instituciones o empresas facilite reunir amenidades deseables en los centros de recarga. Para destinos que impliquen una estad´ıamayor a 2 horas, los cargadores semi-r´apidosson una mejor alternativa debido a que se cuenta con m´astiempo para recargar el veh´ıculo. ⇒ Consideren rutas alternas para mejorar la accesibilidad al centro de recarga. ⇒ Proximidad a las rutas nacionales y cantonales y que exista facilidad de acceso entre las rutas y el punto de recarga. Esta proximidad a la ruta puede ser de 1 a 2 km. ⇒ Ubicaci´oncercana a los centros con mayor poblaci´on;por ejemplo, las cabeceras de cant´on. ⇒ Que los centros de recarga r´apidano se expongan a inundaciones.

Los criterios considerados para el desarrollo de la metodolog´ıade an´alisisde ubicaci´onse identifican en el cuadro anterior (mediante un ). Los dem´ascriterios pueden ser gu´ıaspara que las distribuidoras definan el punto de recarga dentro del ´areadefinida en este trabajo. Para todos los efectos, el criterio tecno-econ´omicoprincipal de ubicaci´ones la disponibilidad de red de distribuci´ontrif´asicaen la zona de inter´es,y su confluencia con la red vial nacional.

2018 3 Propuesta de ubicaci´onde infraestructura de recarga r´apida EPERLab - UCR

1.2. Ubicaciones generales sugeridas durante el taller

Se presentan las tres listas de ubicaci´ondiscutidas: una por mesa. Posteriormente, se presenta la propuesta acumulativa que corresponde a la intersecci´onde las tres propuestas.

Tabla 1: Ubicaciones de inter´esmencionados en la mesa 1

Distribuidora Provincia Lugar GAM/Zona Rural Escaz´u GAM Lindora GAM CNFL San Jos´e Desamparados GAM Curridabat GAM Guadalupe GAM CoopeAlfaroRuiz Alajuela Zarcero Rural Tamarindo Rural Guanacaste CoopeGuanacaste Nicoya Rural Puntarenas Paquera Rural La Fortuna Rural CoopeLesca Alajuela Ciudad Quesada Rural ESPH Heredia Heredia GAM Alajuela GAM Alajuela Bocas de San Carlos Rural Peaje de Naranjo Rural Cartago Turrialba Rural Liberia Rural Limonal Rural Guanacaste Tilar´an Rural Pe˜nasBlancas Rural ICE Lim´on Gu´apiles Rural Lim´on Rural Jac´o Rural Quepos Rural Puntarenas Palmar Norte Rural Ciudad Neily Rural Cruce Barranca Rural San Jos´e Perez Zeled´on Rural Heredia La Virgen de Sarapiqu´ı Rural JASEC Cartago Cartago GAM CoopeSantos Cartago El Empalme Rural

2018 4 Propuesta de ubicaci´onde infraestructura de recarga r´apida EPERLab - UCR

Tabla 2: Ubicaciones de inter´esmencionados en la mesa 2

Distribuidora Provincia Lugar GAM/Zona Rural Sabana GAM CNFL San Jos´e La Uruca GAM Curridabat GAM Huacas Rural CoopeGuanacaste Guanacaste Nicoya Rural San Carlos Rural CoopeLesca Alajuela Sarapiqu´ı Rural ESPH Heredia Heredia GAM Alajuela GAM Alajuela San Ram´on GAM Los Chiles Rural Cartago Turrialba Rural Liberia Rural Guanacaste Limonal Rural Tilar´an Rural Gu´apiles Rural ICE Lim´on Lim´on Rural Sixaola Rural Jac´o Rural Quepos Rural Puntarenas Palmar Norte Rural Barranca Rural San Jos´e San Isidro del General Rural Heredia La Virgen de Sarapiqu´ı Rural Cartago GAM JASEC Cartago Guarco Rural Cartago El Empalme Rural CoopeSantos San Jos´e Macho Gaff Rural

2018 5 Propuesta de ubicaci´onde infraestructura de recarga r´apida EPERLab - UCR

Tabla 3: Ubicaciones de inter´esmencionados en la mesa 3.

Distribuidora Provincia Lugar GAM/Zona Rural Escaz´u GAM La Uruca GAM CNFL San Jos´e Barrio Am´on GAM San Pedro GAM Huacas Rural CoopeGuanacaste Guanacaste Paquera Rural Nicoya Rural La Fortuna Rural Ciudad Quesada Rural CoopeLesca Alajuela San Rosa de Pocosol Rural Puerto Viejo de Sarapiqu´ı Rural Heredia GAM ESPH Heredia Vara Blanca Rural Alajuela - Aeropuerto GAM San Ram´on GAM Alajuela Atenas Rural Orotina Rural San Miguel de Naranjo GAM Cartago Turrialba Rural Liberia Rural Guanacaste Limonal Rural Tilar´an Rural ICE Gu´apiles Rural Lim´on Rural Lim´on Rural Jac´o Rural Quepos Rural Puntarenas Palmar Norte Rural Ciudad Neily Rural Barranca Rural P´erezZeled´on Rural San Jos´e La Georgina Rural Heredia La Virgen de Sarapiqu´ı Rural JASEC Cartago Cartago GAM CoopeSantos Cartago El Empalme Rural

2018 6 Propuesta de ubicaci´onde infraestructura de recarga r´apida EPERLab - UCR

Tabla 4: Puntos de inter´esde todas las mesas de trabajo.

Distribuidora Provincia Lugar GAM/Zona Rural La Uruca GAM Escaz´u GAM San Pedro GAM CNFL San Jos´e Barrio Am´on GAM Sabana GAM Curri GAM CoopeAlfaroRuiz Alajuela Zarcero Rural Huacas Rural Guanacaste Tamarindo Rural CoopeGuanacaste Nicoya Rural Puntarenas Paquera Rural La Fortuna Rural Alajuela Ciudad Quesada Rural CoopeLesca Santa Rosa de Pocosol Rural Heredia Sarapiqu´ı Rural Vara Blanca Rural ESPH Heredia Heredia GAM Alajuela-Aeropuerto GAM San Miguel de Naranjo GAM Atenas Rural Alajuela Orotina Rural San Ram´on GAM Los Chiles Rural Cartago Turrialba Rural Liberia Rural Tilar´an Rural Guanacaste Pe˜nasBlancas Rural S´amara Rural Gu´apiles Rural ICE Siquirres Rural Lim´on Lim´on Rural Puerto Viejo de Lim´on Rural Sixaola Rural Jac´o Rural Quepos Rural Palmar Norte Rural Puntarenas Ciudad Neily Rural Barranca Rural Limonal Rural La Georgina Rural San Jos´e Perez Zeled´on Rural San Isidro Rural Cartago GAM JASEC Cartago Guarco GAM Macho Gaff Rural CoopeSantos Cartago El Empalme Rural

2018 7 Propuesta de ubicaci´onde infraestructura de recarga r´apida EPERLab - UCR 2. Metodolog´ıade ubicaci´onde estaciones de recarga en Costa Rica

El fin de esta metodolog´ıaes responder a la necesidad de identificar una red de recarga r´apida en el corto plazo que pueda asegurar la autonom´ıade un VE modelo para alcanzar cualquier punto geogr´aficocubierto por la red vial nacional. Para lo anterior se utiliz´ola informaci´on georeferenciada del MOPT para identificar las rutas que m´asdemandan autonom´ıapara un VE, es decir, las que conectan a los puntos m´asalejados del pa´ısentre s´ı.Por la forma longitudinal en la que est´adesarrollada la red vial, las rutas suelen confluir hacia la GAM. Por ejemplo, para viajar de la costa Atl´antica a la costa Pac´ıfica,la GAM es actualmente la principal ´areade conexi´on(puede utilizarse tambi´enla ruta 4 entre Gu´apilesy Pto. Viejo de Sarapiqu´ı).Por esta raz´on,viajar de la GAM a Sixaola es una ruta utilizada en el an´alisisde la cobertura de la red de recarga m´ınima.Las rutas analizadas, as´ıcomo la lista de estaciones ubicadas durante los talleres por las cuales se viajar´ıa,se muestran a continuaci´on:

1 San Jos´e- Paso Canoas Empalme - Paso Macho - Georgina - PZ - Palmar Norte - CD Neily - Paso Canoas 2 San Jos´e- Paso Canoas Orotina - Jac´o- Quepos - Palmar Norte - CD Neily - Paso Canoas 3 San Jos´e- Sixaola Gu´apiles- Siquirres - Lim´on- Pto Viejo - Sixaola 4 San Jos´e- Sixaola Turrialba - Siquirres - Lim´on- Pto. Viejo - Sixaola 5 San Jos´e- Limonal Orotina - Barranca 6 San Jos´e- Limonal San Miguel - San Ram´on- Barranca 7 Limonal - Pe˜nasBlancas Limonal - Liberia - Pe˜nasBlancas 8 Limonal - Tamarindo Limonal - Nicoya - Tamarindo 9 Limonal - Tamarindo Limonal - Nicoya - Huacas - Tamarindo 10 Limonal - Mal Pa´ıs Limonal - Paquera - Mal Pa´ıs 11 San Jos´e- Los Chiles Pto Viejo - Santa Rosa 12 San Jos´e- Los Chiles San Miguel - San Ram´on- Los Chiles 13 Limonal - Los Chiles Vara Blanca - Los Chiles

La metodolog´ıa se encarga de identificar cu´alesestaciones eliminar y cu´alesagregar a lo largo de las rutas identificadas. En el an´alisis se identifica una red m´ınima,media y robusta: entre m´asrobusta, m´ascobertura tienen las rutas y menos tiempo de recarga o menos paradas deben realizar los usuarios de VEs. Adem´as,dado que la CNFL tiene 9 cargadores disponibles para suplir su ´areaservida, se realiza un an´alisisseparado para los cantones de la GAM en la secci´on 2.1, el cual tambi´enincluye ubicaciones en las cabeceras de provincia en Alajuela, Heredia y Cartago. La metodolog´ıade c´alculode consumo y de priorizaci´onde ubicaciones a lo largo de la red vial nacional se desarrolla en la secci´on 2.2 y responde al aseguramiento de la autonom´ıa de un veh´ıculoel´ectricomodelo.

2018 8 Propuesta de ubicaci´onde infraestructura de recarga r´apida EPERLab - UCR

2.1. Ubicaci´onde los puntos de recarga en la GAM

A diferencia del an´alisisde autonom´ıaque debe realizarse para ubicar las estaciones de recar- ga en los principales corredores del pa´ıs(presentados anteriormente), la ubicaci´onde las estaciones de recarga en la GAM debe de garantizar cobertura apropiada del ´areade la GAM y debe de considerar el flujo vehicular de las rutas que conectan los centros de poblaci´ony los centros de destino. Adem´as,una infraestructura robusta en la GAM, es tambi´enun buen punto de partida para incentivar los taxis el´ectricos.

La informaci´onde flujo vehicular disponible desde la Secretar´ıade Planificaci´ondel MOPT se revis´odetenidamente y se concluy´oque tiene un gran potencial de uso para la optimizaci´onde ubicaciones de estaciones de recarga cuando exista una alta penetraci´onde VEs. Esta informa- ci´ondebe analizarse con m´etodos estad´ısticosy geogr´aficospara utilizarse adecuadamente en el problema de ubicaci´onde estaciones de recarga del mediano y largo plazo. Sin embargo, para la red de corto plazo, el an´alisisse puede simplificar con informaci´onm´asgeneral. En este trabajo se utiliza un an´alisisde las zonas de origen y destino de viajes vehiculares.

Primeramente, se identificaron los puntos de recarga dentro de la GAM que son fundamenta- les para las personas que se dirigen hacia afuera de ella. Para la cobertura de la red vial nacional son prioritarios los siguientes puntos dentro de la GAM:

Cartago centro. Alajuela centro. Heredia centro.

Estos centros de recarga garantizan una cobertura apropiada al ubicarse en los extremos de la GAM, y al mismo tiempo, para el an´alisisde consumo planteado en la secci´on,eliminan la suposici´onde que todas las personas parten del mismo punto en la GAM. Para priorizar las otras ubicaciones de los puntos de recarga, se defini´ola correlaci´onen- tre los cantones que constituyen la GAM, tomando en cuenta que cada una de ellas presenta caracter´ısticasdistintas. Las caracter´ısticastomadas en cuenta fueron:

1 La cantidad de personas residentes de cada cant´onque trabajan fuera del mismo. 2 La cantidad de veh´ıculosfamiliares en cada cant´on. 3 Las regiones demandantes de movilidad urbana, caracterizadas por poseer los mayores cen- tros de educaci´ony empleadores de la GAM.

Los datos de los puntos 1 y 2 fueron obtenidos del Censo 2011. Estos datos ayudaron a determinar los puntos origen, de los que parte la mayor cantidad de personas, y a reconocer las caracter´ısticasde los mismos. Inicialmente, en la Figura1 se muestra la lista de los cantones con mayor n´umerode habitantes mayores de 15 a˜nosque trabajan fuera de ´el.

2018 9 Propuesta de ubicaci´onde infraestructura de recarga r´apida EPERLab - UCR

Cantones origen

Moravia 14,927 V´azquez 15,349 Tib´as 15,637 Curridabat 15,955 Alajuelita 18,782 Cartago 19,732 La Uni´on 24,204 Heredia 25,561 Goicochea 26,883 Alajuela 33,737 San Jos´e 38,637 Desamparados 45,336

1,5 2 2,5 3 3,5 4 4,5 N´umerode trabajadores(as) ·104

Figura 1: Cantidad de personas mayores de 15 a˜nosque trabajan fuera del cant´onen que habitan.

Seguidamente, todos los cantones se dividieron en siete grandes zonas: sur, norte, este, oeste, centro, noreste, noroeste y suroeste. A continuaci´onse detalla la informaci´onrecopilada para cada una de ellas:

Zonas origen

18,539 Suroeste Trabajadores(as) fuera del cant´on 20,549 Cantidad de veh´ıculosfamiliares 38,216 Noroeste 43,597 35,573 Centro 44,200 34,998 Noreste 62,441 49,178 Este 71,646 34,736 Sur 76,566 54,762 Norte 89,713

1 2 3 4 5 6 7 8 9 Datos seg´unzonas de la GAM ·104

Figura 2: Cantidad de veh´ıculosfamiliares y personas mayores de 15 a˜nosque trabajan fuera de la zona que habitan.

2018 10 Propuesta de ubicaci´onde infraestructura de recarga r´apida EPERLab - UCR

Esta separaci´onse realiz´ocon el objetivo de instalar los centros de recarga en ubicaciones que satisfagan a un grupo amplio de lugares a partir de sus necesidades. En este caso, las necesidades se relacionan en gran medida con la cantidad de veh´ıculosfamiliares en cada una de las zonas, y con las posibles rutas que pueden tomarse hacia los centros de educaci´ono de empleo. Estos puntos de destino se obtuvieron del informe final del estado de la naci´onen desarrollo sostenible: implicaciones sociales, econ´omicas y ambientales del modelo de ciudad vigente en la GAM. Con esto fue posible determinar los puntos destino, es decir, los principales lugares a los que se trasladan las personas. En la Figura3 se pueden apreciar estas zonas que representan las aglomeraciones m´asgrandes de centros de trabajo y estudio de la GAM.

Figura 3: Zonas de destino, aglomeraciones de trabajo y estudio.

Las zonas de destino tambi´ense agruparon en zonas tomando en cuenta su cercan´ıa,con el objetivo de facilitar el proceso de an´alisisde las posibles rutas. En QGIS se designaron pol´ıgonos correspondientes a las zonas de origen (en azul) y destino (en verde). Para la elecci´onde las ubicaciones se valoraron las posibles rutas designadas. Adem´as,se tom´oen cuenta el traslape entre zonas de origen y de destino. El objetivo de la metodolog´ıaes garantizar una cobertura geogr´aficaapropiada en la GAM, pero que esto se realice de forma estrat´egica para satisfacer la demanda de la mejor manera. En las im´agenesya se encuentran se˜naladaslas ubicaciones propuestas restantes para los centros de recarga, y a continuaci´onse justifica su elecci´on.Estos puntos restantes se ubican en la provincia de San Jos´e,y corresponden a la CNFL. Inicialmente, se muestra el an´alisispara la parte sur de la GAM (Desamparados, Alajuelita y Aserr´ı).En la Figura4 se aprecian las rutas que conectan a la zona sur de la GAM con los destinos frecuentes en la misma. Para esta zona se concluye que un punto de recarga es suficiente. Aunque la cantidad de personas que trabajan fuera de los cantones es alta, puede observarse que la cantidad de veh´ıculosfamiliares en cada zona es menor en relaci´oncon el resto de zonas. Al

2018 11 Propuesta de ubicaci´onde infraestructura de recarga r´apida EPERLab - UCR mismo tiempo, estos cantones no se traslapan con zonas de alta movilidad (destinos).

Figura 4: Rutas para la zona sur de la GAM.

De acuerdo con las rutas planteadas, la ubicaci´on´optimapara satisfacer a este grupo de cantones es cercano a Circunvalaci´on.De no ser posible, se recomienda el cant´onde Desamparados, al ser de hecho el cant´oncon la salida m´asgrande de personas a sus trabajos. Por su parte, la zona este abarca los cantones de La Uni´on,Curridabat y Montes de Oca. Esta zona es la segunda en la cantidad de veh´ıculosfamiliares. Para la misma, como se indic´o anteriormente, ya se cuenta con la ubicaci´onprioritaria de Cartago Centro, que abastece a este cant´ony a las personas que se dirigen afuera de la GAM por esta ruta. En la Figura5 se muestran las rutas salientes del este de la GAM.

Figura 5: Rutas para la zona este de la GAM.

2018 12 Propuesta de ubicaci´onde infraestructura de recarga r´apida EPERLab - UCR

Puede observarse que el punto de recarga ubicado en la zona sur es igualmente beneficiosos para este grupo de cantones. Por otro lado, los puntos cercanos a la Universidad de Costa Rica (UCR) y el Colegio Federado de Ingenieros y Arquitectos (CFIA) son valiosos pues se encuentran cercanos a las rutas y en zonas a los que se moviliza un grupo amplio de personas a sus trabajos y centros educativos, existe un traslape entre las rutas y las zonas de destino.

Figura 6: Rutas para la zona noreste de la GAM.

Continuando con el an´alisis,en la zona noreste se incluyeron los cantones de Goicoechea, V´azquezde Coronado, Moravia y San Isidro. En la Figura6 se muestran las rutas. Se nota que no existe coincidencias con las zonas con aglomeraciones de centros educativos y de trabajo. Adem´as, aunque la cantidad de personas que trabajan fuera de ella es alta, la presencia de autom´oviles es muy baja. Con el fin de responder a las necesidades de la zona se recomienda colocar una estaci´on de recarga sobre las rutas que se utilizan para llegar a esta zona de aglomeraci´on.De esta forma, se satisface la demanda de la zona noreste de la GAM, y al mismo tiempo los puntos de recarga pueden ser ´utilespara personas que no van a entrar a ella. Los puntos de la UCR y el CFIA, como puede verse en la imagen, pueden ser utilizados por personas que se dirigen al este de la GAM. As´ı,otro grupo de rutas estar´anacompa˜nadosdel cargador E4, ubicado en Barrio Am´on,y otro grupo de rutas por el cargador E5 en Tib´as.Estos cargadores tambi´enpodr´ıanser utilizados por personas en las zona central y norte respectivo. As´ımismo, el cargador E5 se coloc´oen esta zona porque funciona como punto de partida para las personas que utilizar´anla ruta 32 (San Jos´e- Lim´on). La zona norte inclu´ıaT´ıbasy 6 cantones heredianos. En la Figura7 se muestran las posibles rutas para las personas que salen de esta zona de la GAM. De hecho, esta agrupaci´onde cantones es la que posee el mayor n´umerode habitantes que trabajan fuera de ´el,y tambi´enla mayor cantidad de veh´ıculosfamiliares. Ac´aadem´asya se cuenta con el punto de Heredia Centro, como se mencion´oal inicio de esta secci´on.Este punto no es solamente estrat´egicoporque funciona para las personas que saldr´ande la GAM, sino tambi´enporque se ubica en un una zona que se traslapa con una ´areade empleadores y centros educativos importante.

2018 13 Propuesta de ubicaci´onde infraestructura de recarga r´apida EPERLab - UCR

Figura 7: Rutas para la zona norte de la GAM.

Como puede verse, el punto E5 es efectivamente ´utilpara quienes se dirigen al centro de San Jos´edesde la zona norte de la GAM. Al mismo tiempo, los puntos E6 y E7, ubicados en La Uruca y La Sabana se colocaron para abastecer la demanda en estas rutas que son altamente transitadas, como puede verse en el resto de im´agenes. Por su parte, los cantones de San Jos´ey Moras est´ancontenidos en la zona central de la GAM. Las rutas salientes de estos centros se pueden apreciar en la Figura8.

Figura 8: Rutas para la zona centro de la GAM.

La zona central y sus rutas se encuentran apropiadamente cubiertas. Con los puntos de

2018 14 Propuesta de ubicaci´onde infraestructura de recarga r´apida EPERLab - UCR recarga que se han establecido y discutido anteriormente, un veh´ıculopodr´aencontrar uno en la mayor´ıade sus rutas. Al noroeste se encuentran los cantones de Alajuela, Flores y Bel´en.En la Figura9 se encuentran las rutas.

Figura 9: Rutas para la zona noroestede la GAM.

En el caso del noroeste de la GAM, el punto E9 se coloca aprovechando el traslape entre estos cantones altamente poblados y una importante zona de empleadores. Este punto va a ser provechoso no solamente para las personas que se transporten del noroeste de la GAM hacia el resto de ella, sino tambi´enpara quienes trabajen en la zona o se dirijan hacia fuera de la GAM y necesiten recargar su veh´ıculo.

Figura 10: Rutas para la zona suroeste de la GAM.

2018 15 Propuesta de ubicaci´onde infraestructura de recarga r´apida EPERLab - UCR

Finalmente, los cantones de Escaz´uy Santa Ana componen la zona suroeste. Las rutas para esta zona pueden observarse en la Figura 10. Nuevamente, por el traslape de rutas, para garantizar una cobertura apropiada de la GAM y porque este punto sirve de paso para otra zona industrial, Lindora, se coloca el punto E8 en Escaz´u.La ubicaci´onde este punto es estrat´egicano solamente para las personas que habitan en la zona, sino tambi´enpara la gran cantidad de personas que viajan a esta zona para trabajar.

2.2. Ubicaci´onde los puntos de recarga fuera de la GAM

Para cada una de las rutas analizadas, se realiza un c´alculodel consumo energ´eticode un VE modelo. A partir de perfiles de elevaci´onextra´ıdosde Google Earth entre centros de recarga, se calcula el consumo del VE considerando su autonom´ıay se determina si este puede llegar a la siguiente ubicaci´oncon un cargador. De igual forma, se estima el consumo de una ruta y cu´alesser´ıanlas ubicaciones con m´asdemanda. La estimaci´onde consumo es muy dependiente de la forma de conducci´ondel VE, el cual es afectado por las fuerzas mostradas en la Tabla5 que determinan la potencia mec´anicarequerida para lograr su desplazamiento. El diagrama de la Figura 11 presenta el modelo simplificado del VE: una masa que se desplaza linealmente a una velocidad v. El significado de los par´ametrosque afectan las fuerzas se muestran en la Tabla6. El modelo de consumo explicado a continuaci´onest´abasado en [3] . Los resultados pueden mejorarse con mediciones de velocidad reales y afinaci´ondel modelo.

Tabla 5: Fuerzas involucradas en movimiento lineal del VE.

Frr = µrr m g Fuerza de Resistencia al Rodamiento 2 Fad = (1/2)ρACdv Arrastre Aerodin´amico

Fhc = m g sen(α) Componente del peso del VE.

Fla = 1,05 m a Fuerza requerida para la aceleraci´onlineal; el factor de 1.05 considera la fuerza requerida para dar la aceleraci´onrotacional al motor.

Tabla 6: Par´ametrosque afectan al VE.

g Aceleraci´ongravitacional ρ Densidad del aire

µrr Coeficiente de resistencia al rodamiento A Area´ frontal del VE

Cd Coeficiente de arrastre aerodin´amico α Angulo´ de elevaci´on m Masa del VE v Velocidad del VE Figura 11: Diagrama general del VE. a Aceleraci´onlineal del VE

2018 16 Propuesta de ubicaci´onde infraestructura de recarga r´apida EPERLab - UCR

El par´ametro α puede variar de acuerdo a la elevaci´onde la carretera por la que viaja el VE. Los dem´aspar´ametrosse ajustan a condiciones est´andarde modelado, excepto los marcados en color naranja de la Tabla6: la velocidad v y la aceleraci´on a se estiman a partir de los perfiles de pruebas que utiliza la Agencia de Protecci´onAmbiental de Estados Unidos (EPA del ingl´es Environmental Protection Agency). Se utiliza el perfil gen´erico Urban Dynamometer Driving Schedule (UDDS) [1] (ver Figura 12) para representar rutas de baja velocidad promedio y el perfil Highway Fuel Economy Driving Schedule (HWFET) [2] (ver Figura 12) para rutas de alta velocidad promedio. Con esta informaci´ones posible calcular la fuerza de tracci´ondel VE en cada segundo dada por la ecuaci´on(1).

Fte = Frr + Fad + Fhc + Fla (1)

Figura 12: Perfiles de manejo UDDS y HWFET de EPA.

La Figura 12 muestra las velocidades promedio (Avg Speed) que son 31.5 km/h para UDDS y

2018 17 Propuesta de ubicaci´onde infraestructura de recarga r´apida EPERLab - UCR

77.6 km/h para HWFET. Para el an´alisisde distancias mayores a las totales que se recorren en los perfiles, se extiende el alcance del perfil de conducci´on.Entonces, el perfil de conducci´ones c´ıclico.

En este trabajo se utiliza ´unicamente el perfil UDDS ya que representa el modo de conduc- ci´ont´ıpicoen Costa Rica, y que adem´asbrinda un an´alisism´asconservador. En general, los VEs consumen m´asenerg´ıaa mayor velocidad debido a la fuerza de arrastre del viento, la cual es directamente proporcional - y al cuadrado - de la velocidad del veh´ıculo.En pruebas realizadas para el ajuste del modelo, la autonom´ıadel veh´ıculodisminu´ıacuando la velocidad promedio del perfil HWFET se aumentaba a 93 km/h: aumento de 20 %. Manteniendo la velocidad promedio del perfil HWFET original, el modelo utilizado tiene un rendimiento superior al UDDS. Si se considera que la velocidad promedio en carreteras de alta velocidad ronda los 77.6 km/h para Costa Rica (igual a la del perfil del HWFET), utilizar el perfil UDDS permite asumir un peor caso al considerar consumos mayores y por lo tanto ser mejor para calcular el consumo entre rutas. En un futuro, pueden realizarse barridos de velocidades promedio o variaciones de los perfiles para determinar familias de tipos de conducci´on.

Al calcular la potencia que el VE necesita para desplazarse a una velocidad v con un perfil de conducci´onque permite obtener la aceleraci´on a en cada segundo, la energ´ıaque el VE requiere para su locomoci´onse calcula mediante la ecuaci´on(2), la cual depende del perfil de velocidad. Entonces, la velocidad a la que viaja el VE es un factor determinante en el consumo energ´etico:a mayores velocidades, mayor el consumo. Esta energ´ıaproviene de la bater´ıadel VE. Para llegar a utilizarse en la locomoci´ondel VE, la bater´ıaentrega energ´ıaa su motor y sistema mec´anico.

Entonces, existen las eficiencias mostradas en la Tabla7. La energ´ıa Ea corresponde a aquella utilizada para actividades auxiliares del VE. Finalmente, la ecuaci´on(3) muestra la forma de la estimaci´ondel consumo de la bater´ıasi el VE se impulsa o si tiene una regeneraci´ondebido al frenado o a una conducci´oncuesta abajo. Se considera una eficiencia adicional ηa para restringir la eficiencia neta de la regeneraci´ona un 20 %.

Tabla 7: Consideraciones sobre la eficiencia del sistema electro-mec´anicodel VE.

ηg Eficiencia de los sistemas mec´anicos Valor considerado: 0.9

ηm Eficiencia del motor el´ectrico Valor considerado: 0.9

Eac Energ´ıapara auxiliares Valor considerado: 15 % del consumo en superficie plana en perfil UDDS.

ηg · ηm · ηa Eficiencia neta de la regeneraci´on Valor considerado: 0.2

Z Ete = Fte · v(t)dt (2)

2018 18 Propuesta de ubicaci´onde infraestructura de recarga r´apida EPERLab - UCR

 Ete/(ηg · ηm), Impulsando Ebat = (3) Ete · (ηg · ηm · ηa) + Eac, Regenerando En el futuro puede indagarse con m´asdetalle los par´ametrosde cada fabricante de VE as´ı como las curvas de torque-eficiencia para tener un modelado m´asexacto que se pueda validar con mediciones experimentales de un VE. Los par´ametrosutilizados para este estudio asemejan las caracter´ısticasde un Nissan Leaf y se muestran en la Tabla8: la masa del VE de modelo es la misma que la de este veh´ıculo.Adem´as,se utiliza como referencia la autonom´ıade placa del Nissan Leaf para analizar la efectividad del modelo. Para este c´alculode comparaci´on,se toma valor te´oricoel dato de placa de rendimiento del VE. El valor calculado se obtiene con el modelo de consumo parametrizado con un ´angulode elevaci´on α = 0. El error se muestra en la Tabla9.

Tabla 8: Par´ametrosutilizados para modelo de consumo del VE.

Variable Descripci´on Valor Unidad Masa M 1557 kg

Coeficiente de resistencia al rodamiento Crr 0.02 Adimensional Densidad del aire ρ 1.225 kg/m3 Area´ frontal del VE A 2 m2

Coeficiente de arrastre aerodin´amico Cd 0.5 Adimensional

Tabla 9: Error entre rendimiento de fabricante y calculado con modelo.

Ejemplo: Nissan Leaf Conducci´on en Ciudad Conducci´on en Pista (Perfil (Perfil UDDS) HWFET con 20 % de aumen- to en velocidad promedio) Valor Te´orico 16.75 kWh/100km 20.94 kWh/100km Valor Calculado 17.46 kWh/100km 20.77 kWh/100km Valor Error 4.22 % 0.83 %

Mediante el modelo validado con los datos del fabricante, se realiz´ouna programaci´onde cada una de las rutas de inter´escon su respectiva elevaci´on.Conociendo el valor del ´angulo α en cada secci´onde la ruta, se logra calcular el consumo del VE en toda la ruta y entre cada ubicaci´on.Por ejemplo, la Figura 13 muestra el consumo del VE desde Limonal (como punto de origen) hasta San Jos´e,mientras que la Figura 14 muestra la misma ruta en el sentido contrario. Los consumos calculados en ambos sentidos difieren en 8.2 kWh, lo cual corresponde al 27.33 % de una carga completa de una bater´ıade 30 kWh. Por esta raz´on,el an´alisiscompleto de cada ruta debe realizarse en ambos sentidos y no basta con conocer la distancia entre centros de recarga r´apidade los VEs. Mediante estos valores de consumo, es posible analizar si se puede quitar o si hace falta agregar cargadores para cada ruta analizada.

2018 19 Propuesta de ubicaci´onde infraestructura de recarga r´apida EPERLab - UCR

Para efectos de este estudio, se utiliza un l´ımite de 22.5 kWh (75 % de la capacidad de un VE tipo Nissan Leaf modelo 2017) para identificar si se requiere un cargador entre dos ubicaciones: entonces, si el consumo del VE sobrepasa este valor, el VE presentar´aun riesgo de no alcanzar la pr´oximaubicaci´ondonde puede cargarse.

Figura 13: Consumo de la ruta: Limonal – Barranca – Orotina - San Jos´e.

Figura 14: Consumo de la ruta: San Jos´e– Orotina – Barranca -Limonal.

Con el fin de comprender la pertinencia de cada estaci´onde recarga, se procede a eliminar cada una de las estaciones ubicadas y realizar el an´alisisde consumo. Las Figuras 15y 16 muestran un escenario en el cual se eliminan la ubicaci´onde Orotina para las rutas San Jos´e- Limonal y viceversa, respectivamente. Luego, las Figuras 17y 18 muestran un escenario en el cual se

2018 20 Propuesta de ubicaci´onde infraestructura de recarga r´apida EPERLab - UCR eliminan la ubicaci´onde Barranca para las rutas San Jos´e- Limonal y viceversa, respectivamente. Al eliminar la estaci´onen Orotina, el VE puede desplazarse desde San Jos´ehasta Barranca, y viceversa, sin arriesgar quedarse sin carga. Cuando el VE viaja de Barranca a San Jos´econsume 19.51 kWh; este valor corresponde a 65 % de la energ´ıade una bater´ıatotalmente cargada y sin degeneraci´on.Al eliminar la estaci´onde Barranca, tampoco se presentan problemas.

Figura 15: Consumo de la ruta: San Jos´e– Orotina – Barranca -Limonal.

Figura 16: Consumo de la ruta: Limonal – Barranca – Orotina - San Jos´esin Orotina sin Orotina.

Estos escenarios muestran que Barranca u Orotina pueden omitirse para esta ruta. Sin embar- go, es necesario revisar otros caminos que contengan estas ubicaciones para no afectar a usuarios

2018 21 Propuesta de ubicaci´onde infraestructura de recarga r´apida EPERLab - UCR que viajen por las rutas alternativas. Adem´as,las estaciones que se pueden descartar son ne- cesarias si la otra se elimina. Entonces, el objetivo principal de la metodolog´ıaes cuantificar la cantidad de cargadores que se requieren para lograr una correcta cobertura de autonom´ıapar VEs y mostrar cu´alesson las opciones para una red m´ınima, una media y otra robusta. En algunos casos, como en el ejemplo (Figuras 15a 18), la decisi´onde mantener una opci´onde ubicaci´onse reduce a la cantidad de tiempo que el usuario requerir´ıapara cargar su VE, el cual es proporcional al consumo en kWh. El an´alisiscomprehensivo de los resultados se muestra en la secci´on3.

Figura 17: Consumo de la ruta: San Jos´e– Orotina – Barranca -Limonal sin Barranca.

Figura 18: Consumo de la ruta: Limonal – Barranca – Orotina - San Jos´esin Barranca.

Con el fin de analizar la ubicaci´onde las estaciones de recarga, tomando en cuenta la capa-

2018 22 Propuesta de ubicaci´onde infraestructura de recarga r´apida EPERLab - UCR cidad de las bater´ıasde los VEs que se comercializan actualmente (30 kWh) y de aquellos que ya circulan en el pa´ıs,se establece arbitrariamente un consumo m´aximoequivalente al 75 % de estas bater´ıas.De tal forma que el an´alisissiguiente busca asegurar que el consumo entre estaciones de recarga no supere los 22.5 kWh. Una red m´asrobusta puede servir a VEs que tienen una capaci- dad nominal menor o tienen una degradaci´on,lo cual disminuye la capacidad real de la bater´ıa. La Tabla 10 muestra una lista de modelos de VEs con sus respectivas capacidades nominales en kWh [4]. La mayor´ıade los VEs supera la marca de 22.5 kWh, especialmente si estos VEs son nuevos. Existe una tendencia creciente para la capacidad de almacenamiento de las bater´ıasde VEs, por lo que el 22.5 kWh es un valor que en futuro obligar´aa un sobre-dimensionamiento de la capacidad de la red de recarga, aunque cumple de forma responsable la necesidad actual de una red de corto plazo. Adem´as,el cargador r´apidopuede lograr alcanzar el estado de carga de la bater´ıaa 80 % de su capacidad real. Este valor corresponde a 24 kWh para bater´ıasde 30 kWh. Por esta raz´on,utilizar un l´ımitede 22.5 kWh es m´asexigente para el an´alisisde ubicaci´on.

Tabla 10: Capacidades nominales de bater´ıasde VEs y precios base de los modelos.

Modelo Bater´ıa Rango Precio (base y Promedios Promedios (kWh) (km) sin subsidios) Bater´ıa Precio Tesla Model S 100 D 100 539 $94,000 Tesla Model S P100D 100 507 $165,000 $99,857 Tesla Model X 100D 100 475 $96,000 89 kWh (base y sin Tesla Model X P100D 100 465 $140,000 subsidios) Tesla Model S 75 D 75 417 $74,500 Tesla Model X 75D 75 381 $79,500 Chevrolet Bolt EV 60 383 $36,620 Renault Zoe 41 299 $31,000 Nissan Leaf Gen 2 40 243 $29,990 $33,288 Volkswagen e-Golf 35.8 191 $30,495 38 kWh (base y sin Ford Focus Electric 33.5 185 $29,120 subsidios) BMW i3 33 183 $44,450 Kia Soul EV 30 179 $33,950 Nissan Leaf Gen 1 30 172 $30,680 Hyundai Ioniq 28 200 $29,500 $31,589 Honda Clarity 25.5 143 $37,150 23 kWh (base y sin Fiat 500 e 24 135 $32,995 subsidios) 22.5 Volkswagen e-Up! 18.7 159 $34,500 Smart ED 17.6 161 $23,800

2018 23 Propuesta de ubicaci´onde infraestructura de recarga r´apida EPERLab - UCR 3. Verificaci´onde distancias y an´alisisde resultados

Los puntos que se proponen en la Tabla4 son ubicaciones con un cargador r´apidopara suplir de energ´ıaa un VE suponiendo que pueda lograr cualquier viaje dentro del pa´ıs.Se consideran ´unicamente las rutas de la red vial nacional, lo cu´alobligar´ıaa que el viaje sea planeado con antelaci´onpara evitar que el VE se quede sin energ´ıaen medio de un trayecto alternativo. Para cubrir todas las rutas posibles entre caminos cantonales, ser´ıanecesario un extenso an´alisiscombi- nacional de rutas secundarias, lo cual puede responder a una red de recarga r´apidaen el mediano plazo. No obstante, para la red de corto plazo, desplazarse por la red vial nacional implica una cobertura total entre or´ıgenesy destinos de viaje principales del pa´ıs.En esta secci´onse muestran las principales rutas que se pueden definir para asegurar cobertura en todo el pa´ıscon los puntos de recarga. Adem´as,se analiza si estos cumplen con los requisitos m´ınimosde distancia de la ley actual y con las necesidades de autonom´ıaplanteadas en la secci´on 2.2.

Las Figuras 21a 41 muestran los c´alculosrealizados para cada una de las rutas analizadas. En las Secciones 3.1a 3.4 se describe el an´alisisde cada uno de los escenarios analizados: a partir de un an´alisiscompleto de cada escenario de colocaci´onde estaciones se logra deducir la necesidad de la ubicaci´oncorrespondiente. Los datos mostrados en las Figuras pueden tener una marca de color, la cual se codifica seg´unla Tabla 11.

Tabla 11: Codificaci´onde colores para an´alisisde escenarios de ubicaci´on.

C´odigode Color Significado Naranja Se elimina la estaci´onde la ubicaci´onmarcada. Esta ubicaci´onse puede catalogar como un destino o un punto final de la ruta analizada. Amarillo Se elimina la estaci´onde la ubicaci´onmarcada. Esta ubicaci´onse encuen- tra a lo largo de una ruta analizada, mas no corresponde a un destino o punto final. Rojo Indica que al eliminar una estaci´on,el VE se acerca al consumo l´ımite de los 22.5 kWh. Verde Indica que al eliminar una estaci´on, el VE se aleja del consumo l´ımite establecido en 22.5 kWh.

3.1. Rutas de la GAM a Paso Canoas

Para cubrir la ruta de San Jos´ea Paso Canoas se analizan dos rutas posibles:

Ruta 1 SJ - Paso Canoas Empalme - Paso Macho - Georgina - P´erezZeled´on- Palmar Norte - CD Neily - Paso Canoas Ruta 2 SJ - Paso Canoas Orotina - Jac´o- Quepos - Palmar Norte - CD Neily - Paso Canoas

2018 24 Propuesta de ubicaci´onde infraestructura de recarga r´apida EPERLab - UCR

3.4 km

54.4 km N 30.3 km

40.3 km 26.5 km 9.8 km 3.4 km

14.8 km 33.7 km

64.8 km

35.2 km

N 103 km

95.4 km

75.5 km

Figura 19: Rutas desde San Jos´ehasta ciudad Neily.

En la Figura 19 se puede ver que existen dos rutas para llegar desde la GAM hasta Ciu- dad Neily. Se puede ver que de P´erezZeled´ona Palmar Norte y que de Quepos a Palmar Norte la distancia entre estaciones supera la distancia m´aximapermitida por la ley, debido a que se encuentran a 95.4 km y 103 km respectivamente, por lo que seg´unla ley ser´ıanecesario que se coloque una cerca de Dominical. Sin embargo, no es posible determinar si es necesaria otra esta-

2018 25 Propuesta de ubicaci´onde infraestructura de recarga r´apida EPERLab - UCR ci´onsin un an´alisisde autonom´ıa(o sea, de consumo de energ´ıa)del VE.

En la Figura 20 se muestra el an´alisisde autonom´ıa del VE. Sobre la ruta por el cerro de la Muerte, se puede observar que al eliminar las estaciones de Paso Macho (amarillo) y de Ciudad Neily (naranja) la autonom´ıa del VE no se ve comprometida, ya que los valores de energ´ıaconsumidos son menores a 22.5 kWh. En el caso que se desee eliminar las estaciones de Paso Macho y de la Georgina el consumo del VE llega a 21 kWh (ver Figura 21).

RUTA 1 IDA SJ - EMPALME - PASO MACHO - GEORGINA - PZ - PALMAR NORTE - CIUDAD NEILY - PASO CANOAS ESTACIONES DISTANCIA (km) CONSUMO (kWh) (INICIO -ESTAC) CONSUMO (kWh) (ESTAC -ESTAC) CONSUMO (kWh) (ESTAC -ESTAC) Eliminación SAN JOSÉ 0 0 0 EMPALME 58.25 16.12 16.12 16.12 PASO MACHO 78.84 21.92 5.80 GEORGINA 96.54 25.05 3.13 8.93 PZ 130.81 25.58 0.53 0.53 PALMAR NORTE 222.04 40.06 14.48 14.48 CIUDAD NEILLY 294.45 52.50 12.44 PASO CANOAS 300.88 53.53 1.03 13.47

RUTA 1 VUELTA PASO CANOAS - CIUDAD NEILY - PALMAR NORTE - PZ- GEORGINA - PASO MACHO - EMPALME - SJ ESTACIONES DISTANCIA (km) CONSUMO (kWh) (INICIO -ESTAC) CONSUMO (kWh) (ESTAC -ESTAC) CONSUMO (kWh) (ESTAC -ESTAC) Eliminación PASO CANOAS 0.00 0.00 0 0 CIUDAD NEILLY 6.44 1.10 1.10 PALMAR NORTE 78.84 13.46 12.36 13.46 PZ 170.07 33.74 20.28 20.28 GEORGINA 204.34 48.86 15.12 PASO MACHO 222.04 52.89 4.03 15.12 EMPALME 242.64 54.73 1.84 5.87 SAN JOSÉ 300.88 59.79 5.05 5.05

Figura 20: C´alculosde la ruta 1 San Jos´e- Paso Canoas escenario 1.

RUTA 1 IDA EORGINA - PZ - PALMAR NORTE - CIUDAD NEILY - PASO CANOAS ESTACIONES DISTANCIA (km) CONSUMO (kWh) (INICIO -ESTAC) CONSUMO (kWh) (ESTAC -ESTAC) CONSUMO (kWh) (ESTAC -ESTAC) Eliminación SAN JOSÉ 0 0 0 0 EMPALME 58.25 16.12 16.12 16.12 PASO MACHO 78.84 21.92 5.80 GEORGINA 96.54 25.05 3.13 PZ 130.81 25.58 0.53 9.46 PALMAR NORTE 222.04 40.06 14.48 14.48 CIUDAD NEILLY 294.45 52.50 12.44 PASO CANOAS 300.88 53.53 1.03 13.47

RUTA 1 VUELTA ALMAR NORTE - PZ- GEORGINA - PASO MACHO - EMPALME - SJ ESTACIONES DISTANCIA (km) CONSUMO (kWh) (INICIO -ESTAC) CONSUMO (kWh) (ESTAC -ESTAC) CONSUMO (kWh) (ESTAC -ESTAC) Eliminación PASO CANOAS 0.00 0.00 0 0 CIUDAD NEILLY 6.44 1.10 1.10 PALMAR NORTE 78.84 13.46 12.36 13.46 PZ 170.07 33.74 20.28 20.28 GEORGINA 204.34 48.86 15.12 PASO MACHO 222.04 52.89 4.03 EMPALME 242.64 54.73 1.84 21.00 SAN JOSÉ 300.88 59.79 5.05 5.05

Figura 21: C´alculosde la ruta 1 San Jos´e- Paso Canoas escenario 2.

2018 26 Propuesta de ubicaci´onde infraestructura de recarga r´apida EPERLab - UCR

RUTA 2 IDA SJ - OROTINA - JACO - QUEPOS - PALMAR NORTE - CIUDAD NEILLY - PASO CANOAS ESTACIONES DISTANCIA (km) CONSUMO (kWh) (INICIO -ESTAC) CONSUMO (kWh) (ESTAC -ESTAC) CONSUMO (kWh) (ESTAC -ESTAC) Eliminación SAN JOSÉ 0 0 0.00 0 OROTINA 52.45 5.74 5.74 JACO 94.77 12.39 6.64 12.39 QUEPOS 158.16 23.38 11.00 11 PALMAR NORTE 260.98 40.80 17.42 17.42 CIUDAD NEILLY 334.67 53.37 12.57 PASO CANOAS 344.33 55.09 1.72 14.29

RUTA 2 VUELTA PASO CANOAS - CIUDAD NEILLY - PALMAR NORTE - QUEPOS - JACO - OROTINA -SJ ESTACIONES DISTANCIA (km) CONSUMO (kWh) (INICIO -ESTAC) CONSUMO (kWh) (ESTAC -ESTAC) CONSUMO (kWh) (ESTAC -ESTAC) Eliminación PASO CANOAS 0.00 0.00 0.00 0 CIUDAD NEILLY 9.65 1.65 1.65 PALMAR NORTE 83.35 14.18 12.53 14.18 QUEPOS 186.16 31.71 17.53 17.53 JACO 249.56 42.41 10.70 10.7 OROTINA 291.87 50.73 8.31 SAN JOSÉ 344.33 63.51 12.78 21.10

Figura 22: C´alculosde la ruta 2 San Jos´e- Paso Canoas escenario 1.

La segunda ruta solo tiene un escenario y es quitar la estaci´onde Orotina, debido a que si se eliminan el resto de estaciones los consumos son muy elevados. Al quitar la estaci´onde Orotina de regreso a la capital el consumo se vuelve de 21.1 kWh, estando por debajo del m´aximoconsumo permitido, pero muy cercano a este.

3.2. Rutas a Sixaola

Para cubrir la ruta de San Jos´ea Sixaola se analizan dos rutas posibles:

Ruta 1 SJ - Sixaola Gu´apiles- Siquirres - Lim´on- Pto Viejo - Sixaola Ruta 2 SJ - Sixaola Turrialba - Siquirres - Lim´on- Pto. Viejo - Sixaola

La Figura 23 muestra la cobertura en distancias de las rutas y los cargadores ubicados a lo largo de estas. El viaje desde San Jos´ehasta Sixaola se puede realizar por Gu´apileso por Turrial- ba. En ambos casos, las estaciones se encuentran a una distancia menor de la m´aximapermitida seg´unla ley. Desde Siquirres hasta Lim´ontambi´ense cumple el requisito de distancia m´axima.

Desde Lim´onhasta Sixaola, las estaciones se encuentran en lugares que permite al VE abas- tecerse cuando lo necesita, es decir, hay cobertura para recarga. Para llegar al centro de recarga de Puerto Viejo es necesario dar un desv´ıode 5 km desde la ruta 36, lo cual permite que el VE llegue sin problemas hasta su destino.

2018 27 Propuesta de ubicaci´onde infraestructura de recarga r´apida EPERLab - UCR

N 34.3 km

67.2 km

59.1 km

46.1 km

41.1 km

N

60 km 37 km 5 km

Figura 23: Rutas desde San Jos´ehasta Sixaola.

RUTA 1 IDA SJ - GUAPILES - SIQUIRRES - LIMON - PTO VIEJO - SIXAOLA CONSUMO (kWh) (ESTAC -ESTAC) ESTACIONES DISTANCIA (km) CONSUMO (kWh) (INICIO -ESTAC) CONSUMO (kWh) (ESTAC -ESTAC) ELIMINACIÓN GUAPILES 67.27 10.46 10.46 10.46 SIQUIRRES 100.91 15.52 5.05 LIMON 160.93 25.29 9.78 14.83 PTO VIEJO 217.26 34.76 9.47 SIXAOLA 246.23 39.75 4.99 14.46

RUTA 1 VUELTA SIXAOLA -PTO VIEJO - LIMON - SIQUIRRES - GUAPILES-SJ CONSUMO (kWh) (ESTAC -ESTAC) ESTACIONES DISTANCIA (km) CONSUMO (kWh) (INICIO -ESTAC) CONSUMO (kWh) (ESTAC -ESTAC) ELIMINACIÓN PTO VIEJO 28.97 5.03 5.03 LIMON 85.30 14.65 9.62 14.65 SIQUIRRES 145.32 25.38 10.72 GUAPILES 178.96 31.71 6.33 17.05 SAN JOSE 246.23 47.72 16.01 16.01

Figura 24: C´alculosde la ruta 1 San Jos´e- Sixaola escenario 1.

En la Figura 24 se puede ver que en la ruta de San Jos´ehasta Sixaola por Gu´apiles,se pueden eliminar las estaciones que se encuentran en Siquirres y en Puerto Viejo. El VE consume menos de los 22.5 kWh, por lo que en este escenario el VE no pierde autonom´ıaen su trayecto. Para esta misma ruta se considera un segundo escenario en el cual las estaciones de recarga que se eliminan son Puerto Viejo y Gu´apiles,pero para este caso si se elimina Gu´apiles el VE llega a

2018 28 Propuesta de ubicaci´onde infraestructura de recarga r´apida EPERLab - UCR un consumo de 22.23 kWh hasta San Jos´e,por lo que est´amuy cercano del consumo m´aximoy no ser´ıaconveniente quitar la estaci´onque es encuentra en Gu´apiles.

RUTA 1 IDA SJ - GUAPILES - SIQUIRRES - LIMON - PTO VIEJO - SIXAOLA CONSUMO (kWh) (ESTAC -ESTAC) ESTACIONES DISTANCIA (km) CONSUMO (kWh) (INICIO -ESTAC) CONSUMO (kWh) (ESTAC -ESTAC) ELIMINACIÓN GUAPILES 67.27 10.46 10.46 SIQUIRRES 100.91 15.52 5.05 15.52 LIMON 160.93 25.29 9.78 9.78 PTO VIEJO 217.26 34.76 9.47 SIXAOLA 246.23 39.75 4.99 14.46

RUTA 1 VUELTA SIXAOLA -PTO VIEJO - LIMON - SIQUIRRES - GUAPILES-SJ CONSUMO (kWh) (ESTAC -ESTAC) ESTACIONES DISTANCIA (km) CONSUMO (kWh) (INICIO -ESTAC) CONSUMO (kWh) (ESTAC -ESTAC) ELIMINACIÓN PTO VIEJO 28.97 5.03 5.03 LIMON 85.30 14.65 9.62 14.65 SIQUIRRES 145.32 25.38 10.72 10.72 GUAPILES 178.96 31.71 6.33 SAN JOSE 246.23 47.72 16.01 22.35

Figura 25: C´alculosde la ruta 1 San Jos´e- Sixaola escenario 2.

Para el an´alisisde la cantidad m´ınimade estaciones de recarga considerando la autonom´ıa del VE modelo en la ruta 2, se realizaron dos escenarios que se muestran en las Figuras 26y 27, respectivamente. En el escenario 1, se elimina la estaci´onSiquirres (marcada en amarillo) y las estaciones de destino Pto. Viejo y Sixaola (marcadas en naranja). El VE entre Turrialba y Lim´on consume una energ´ıaestimada de 18.92 kWh. Entonces, el VE consume menos de 22.5 kWh y no compromete su autonom´ıa.Las estaciones en los puntos de destino pueden eliminarse y el VE consumir´a14.4 kWh, por lo que no compromete su autonom´ıa.Al regreso de la ruta, dado que el VE parte con una carga llena desde Sixaola, llega a San Jos´esin inconvenientes. Entonces, es conveniente mantener tanto Gu´apilescomo Sixaola.

RUTA 2 IDA SJ - TURRIALBA - SIQUIRRES - LIMON - PTO VIEJO - SIXAOLA CONSUMO (kWh) (ESTAC -ESTAC) ESTACIONES DISTANCIA (km) CONSUMO (kWh) (INICIO -ESTAC) CONSUMO (kWh) (ESTAC -ESTAC) ELIMINACIÓN TURRIALBA 66.14 10.64 10.64 10.64 SIQUIRRES 111.21 19.37 8.72 LIMON 172.20 29.57 10.20 18.92 PTO VIEJO 223.70 38.33 8.76 SIXAOLA 257.49 44.01 5.68 14.44

RUTA 2 VUELTA SJ - TURRIALBA - SIQUIRRES - LIMON - PTO VIEJO - SIXAOLA CONSUMO (kWh) (ESTAC -ESTAC) ESTACIONES DISTANCIA (km) CONSUMO (kWh) (INICIO -ESTAC) CONSUMO (kWh) (ESTAC -ESTAC) ELIMINACIÓN PTO VIEJO 33.80 5.79 5.79 LIMON 85.30 14.68 8.89 14.68 SIQUIRRES 146.29 25.25 10.57 TURRIALBA 191.35 32.84 7.59 18.16 SAN JOSE 257.50 44.51 11.67 11.67

Figura 26: C´alculosde la ruta 2 San Jos´e- Sixaola escenario 1.

En el escenario 2, se elimina la estaci´onde Turrialba (marcada en amarillo), as´ıcomo las

2018 29 Propuesta de ubicaci´onde infraestructura de recarga r´apida EPERLab - UCR de destino marcadas en el escenario 1. En este caso, de San Jos´ea Siquirres se consume 19.37 kWh, inferior a la marca de 22.5 kWh para establecer satisfacci´onde autonom´ıa.A la vuelta, de Siquirres a San Jos´e,el VE consume 19.26 kWh, inferior a la marca de 22.5 kWh. No obstante, si se considera a estaci´onde Cartago, este consumo es menor. En general, se recomienda mantener la estaci´onde Siquirres y eliminar la de Turrialba considerando la ubicaci´onestrat´egicade Cartago. No obstante, para mejorar la robustez de la red de recarga, se puede recomendar incluir Turrialba como punto estrat´egicoy disminuir la ansiedad de recarga de los usuarios de VEs.

RUTA 2 IDA SJ - TURRIALBA - SIQUIRRES - LIMON - PTO VIEJO - SIXAOLA CONSUMO (kWh) (ESTAC -ESTAC) ESTACIONES DISTANCIA (km) CONSUMO (kWh) (INICIO -ESTAC) CONSUMO (kWh) (ESTAC -ESTAC) ELIMINACIÓN TURRIALBA 66.14 10.64 10.64 SIQUIRRES 111.21 19.37 8.72 19.37 LIMON 172.20 29.57 10.20 10.20 PTO VIEJO 223.70 38.33 8.76 SIXAOLA 257.49 44.01 5.68 14.44

RUTA 2 VUELTA SIXAOLA -PTO VIEJO - LIMON - SIQUIRRES - TURRIALBA-SJ CONSUMO (kWh) (ESTAC -ESTAC) ESTACIONES DISTANCIA (km) CONSUMO (kWh) (INICIO -ESTAC) CONSUMO (kWh) (ESTAC -ESTAC) ELIMINACIÓN PTO VIEJO 33.80 5.79 5.79 LIMON 85.30 14.68 8.89 14.68 SIQUIRRES 146.29 25.25 10.57 10.57 TURRIALBA 191.35 32.84 7.59 SAN JOSE 257.50 44.51 11.67 19.26

Figura 27: C´alculosde la ruta San Jos´e- Sixaola escenario 2.

3.3. Rutas a Guanacaste

Para cubrir la ruta de San Jos´ea Guanacaste se analizan varias rutas posibles:

Ruta 1 SJ - Limonal Orotina - Barranca Ruta 2 SJ - Limonal San Miguel - San Ram´on- Barranca Ruta 1 Limonal - Pe˜nasBlancas Limonal - Liberia - Pe˜nasBlancas Ruta 1 Limonal - Tamarindo Limonal - Nicoya - Tamarindo Ruta 2 Limonal - Tamarindo Limonal - Nicoya - Huacas - Tamarindo Ruta 1 Limonal - Mal Pa´ıs Limonal - Paquera - Mal Pa´ıs

En la Figura 28 se puede ver que de la capital del pa´ıshacia Limonal se puede llegar por dos rutas en las que un VE puede recargar sin problemas en varias estaciones durante todo su recorrido. Una de las rutas toma por Orotina y la otra por San Ram´on, pero en ambos casos existen suficiente estaciones para suplir de energ´ıay el adem´asse muestra la distancia del centro de recarga de Atenas por si fuera necesario su uso.

2018 30 Propuesta de ubicaci´onde infraestructura de recarga r´apida EPERLab - UCR

N

50 km 10.1 km

38.2 km 31.2 km

23 km 16.5 km 50.3 km 37 km

N

63.4 km 62.9 km

92.9 km

41 km 51.1 km 8.3 km

N

23 km

77.1 km 70.3 km

Figura 28: Rutas desde San Jos´ea los distintos puntos en Guanacaste.

Limonal es el punto de inflexi´onpara ir a las distintos destinos en Guanacaste. Para los casos de Tamarindo y de Pe˜nasBlancas las estaciones cumplen con la distancia requerida, pero si se

2018 31 Propuesta de ubicaci´onde infraestructura de recarga r´apida EPERLab - UCR deseara ir a Mal Pa´ıshabr´ıaun inconveniente en el trayecto de Limonal (o Nicoya) hasta Paquera, porque el recorrido que se debe hacer para llegar hasta Paquera supera los 80 km permitidos, siendo necesario que se coloque un centro de recarga cercano a Puerto Thiel (esto solo haciendo la consideraci´onde distancia m´axima).

Para llegar de San Jos´ea Limonal se analizan los dos rutas posibles, cuyos datos rendimiento seg´untopograf´ıay distancia se pueden ver en la Figura 29. En esta Figura se puede ver que a pesar de las distancias que se recorren la autonom´ıadel VE no se ven comprometida en el caso de que se elimine al estaci´onde recarga que se encuentra en barranca, debido a que tanto en la ida como en la vuelta del viaje el consumo m´aximoentre estaciones es de 14.35 kWh y 15 kWh respectivamente.

RUTA 1 IDA SJ - OROTINA - BARRANCA - LIMONAL ESTACIONES DISTANCIA (km) CONSUMO (kWh) (INICIO -ESTAC) CONSUMO (kWh) (ESTAC -ESTAC) CONSUMO (kWh) (ESTAC -ESTAC) ELIMINACION OROTINA 53.11 5.50 5.50 5.50 BARRANCA 87.87 11.08 5.58 LIMONAL 137.92 19.85 8.77 14.35

RUTA 1 VUELTA LIMONAL - BARRANCA - OROTINA - SJ ESTACIONES DISTANCIA (km) CONSUMO (kWh) (INICIO -ESTAC) CONSUMO (kWh) (ESTAC -ESTAC) CONSUMO (kWh) (ESTAC -ESTAC) ELIMINACION BARRANCA 50.05 8.54 8.54 OROTINA 84.81 15.00 6.46 15.00 SAN JOSÉ 137.92 28.05 13.05

Figura 29: C´alculosde la ruta 1 San Jos´ea Limonal escenario 1.

El segundo escenario de esta misma ruta es que se elimine la estaci´onde Orotina, pero en este caso en el viaje de vuelta el consumo del VE llega a los 19.51 kWh (ver Figura 30), lo cual es cercano a el m´aximode consumo que se desea, por lo que no es recomendable que se elimine ya que se estar´ıacomprometiendo la autonom´ıadel VE.

RUTA 1 IDA SJ - OROTINA - BARRANCA - LIMONAL ESTACIONES DISTANCIA (km) CONSUMO (kWh) (INICIO -ESTAC) CONSUMO (kWh) (ESTAC -ESTAC) CONSUMO (kWh) (ESTAC -ESTAC) ELIMINACION OROTINA 53.11 5.50 5.50 8.77 BARRANCA 87.87 11.08 5.58 11.08 LIMONAL 137.92 19.85 8.77

RUTA 1 VUELTA LIMONAL - BARRANCA - OROTINA - SJ ESTACIONES DISTANCIA (km) CONSUMO (kWh) (INICIO -ESTAC) CONSUMO (kWh) (ESTAC -ESTAC) CONSUMO (kWh) (ESTAC -ESTAC) ELIMINACION BARRANCA 50.05 8.54 8.54 8.54 OROTINA 84.81 15.00 6.46 SAN JOSÉ 137.92 28.05 13.05 19.51

Figura 30: C´alculosde la ruta 1 San Jos´ea Limonal escenario 2.

La segunda ruta para viajar de la capital hasta Limonal muestra en la Figura 31 que no es recomendable quitar la estaci´onde Barranca si se toma esta ruta y que la estaci´onde San Miguel podr´ıaser eliminada, debido a que el consumo m´aximoentre estaciones sin esta, llegar´ıaa un m´aximode 10.84 kWh.

2018 32 Propuesta de ubicaci´onde infraestructura de recarga r´apida EPERLab - UCR

RUTA 2 IDA SJ - S. MIGUEL - S. RAMÒN - BARRANCA - LIMONAL ESTACIONES DISTANCIA (km) CONSUMO (kWh) (INICIO -ESTAC) CONSUMO (kWh) (ESTAC -ESTAC) CONSUMO (kWh) (ESTAC -ESTAC) ELIMINACION S. MIGUEL 46.67 7.90 7.90 S. RAMÒN 56.81 10.76 2.86 10.76 BARRANCA 93.83 13.76 3.00 LIMONAL 144.52 22.30 8.54 11.54

RUTA 2 VUELTA LIMONAL - BARRANCA - S. RAMÒN - S. MIGUEL- SJ ESTACIONES DISTANCIA (km) CONSUMO (kWh) (INICIO -ESTAC) CONSUMO (kWh) (ESTAC -ESTAC) CONSUMO (kWh) (ESTAC -ESTAC) ELIMINACION BARRANCA 50.70 8.94 8.94 S. RAMÒN 87.71 19.48 10.54 19.48 S. MIGUEL 97.85 21.64 2.16 SJ 144.52 30.32 8.68 10.84

Figura 31: C´alculosde la ruta 2 San Jos´ea Limonal escenario 1.

RUTA 2 IDA SJ - S. MIGUEL - S. RAMÒN - BARRANCA - LIMONAL ESTACIONES DISTANCIA (km) CONSUMO (kWh) (INICIO -ESTAC) CONSUMO (kWh) (ESTAC -ESTAC) CONSUMO (kWh) (ESTAC -ESTAC) ELIMINACION S. MIGUEL 46.67 7.90 7.90 7.9 S. RAMÒN 56.81 10.76 2.86 BARRANCA 93.83 13.76 3.00 5.86 LIMONAL 144.52 22.30 8.54 8.54

RUTA 2 VUELTA LIMONAL - BARRANCA - S. RAMÒN - S. MIGUEL- SJ ESTACIONES DISTANCIA (km) CONSUMO (kWh) (INICIO -ESTAC) CONSUMO (kWh) (ESTAC -ESTAC) CONSUMO (kWh) (ESTAC -ESTAC) ELIMINACION BARRANCA 50.70 8.94 8.94 8.94 S. RAMÒN 87.71 19.48 10.54 S. MIGUEL 97.85 21.64 2.16 12.71 SJ 144.52 30.32 8.68 8.68

Figura 32: C´alculosde la ruta 2 San Jos´ea Limonal escenario 2.

El escenario 2 de la segunda ruta es eliminar la estaci´onque se encuentra en San Ram´on.En la Figura 32 se puede ver que al ser eliminada esta estaci´onel VE no pierde autonom´ıaen la ida ni en la vuelta, llegando a consumos de 5.86 kWh y 12.71 kWh respectivamente.

RUTA 1 IDA Limonal-Liberia-Peñas Blancas ESTACIONES DISTANCIA (km) CONSUMO (kWh) (INICIO -ESTAC) CONSUMO (kWh) (ESTAC -ESTAC) Liberia 69.20 11.97 11.97 Peñas Blancas 145.48 24.87 12.91

RUTA 1 VUELTA Peñas Blancas-Liberia-Limonal ESTACIONES DISTANCIA (km) CONSUMO (kWh) (INICIO -ESTAC) CONSUMO (kWh) (ESTAC -ESTAC) Liberia 76.28 13.44 13.44 Limonal 145.48 25.19 11.75

Figura 33: C´alculosde la ruta 1 Limonal a Pe˜nasBlancas escenario 1.

Para la ruta de Limonal a Pe˜nasBlancas se tiene solo un escenario y para este no se pueden eliminar estaciones de recarga, debido a que solo existe una intermedia y la distancia y el consumo no lo permiten (ver Figura 33). Para la ruta de Limonal hasta Mal Pa´ıstambi´enpresenta solo un escenario que se puede ver en la Figura 34, donde se muestra que el consumo entre las estaciones es menor a los 22.5 kWh permitidos y logrando una autonom´ıaen todo el trayecto. No se toma en consideraci´oneliminar alguna estaci´onde recarga por la gran distancia entre las mismas y que el consumo del VE se elevar´ıademasiado.

2018 33 Propuesta de ubicaci´onde infraestructura de recarga r´apida EPERLab - UCR

RUTA 1 IDA Limonal-Paquera-Mal país ESTACIONES DISTANCIA (km) CONSUMO (kWh) (INICIO -ESTAC) CONSUMO (kWh) (ESTAC -ESTAC) PAQUERA 127.14 21.69 21.69 MALPAÌS 177.03 30.79 9.10

RUTA 1 VUELTA Mal país-Paquera-Limonal ESTACIONES DISTANCIA (km) CONSUMO (kWh) (INICIO -ESTAC) CONSUMO (kWh) (ESTAC -ESTAC) PAQUERA 49.89 8.32 8.32 LIMONAL 177.03 30.57 22.26

Figura 34: C´alculosde la ruta 1 Limonal a Mal Pa´ısescenario 1.

Para ir de Limonal a Tamarindo se tienen dos rutas. En la primera se considera eliminar la estaci´onque se encuentra en Tamarindo, logrando que la autonom´ıadel VE no se vea perjudicada, tal y como se puede ver en la Figura 35. Adem´asse tiene que en la ruta 2 hacia Tamarindo se puede eliminar la estaci´onde recarga que se encuentra en Huacas, siendo esta muy cercana a la de Tamarindo y no necesaria para la autonom´ıadel VE podr´ıaser eliminada de la red sin mayores repercusiones (ver Figura 36).

RUTA 1 IDA LIMONAL-NICOYA-TAMARINDO ESTACIONES DISTANCIA (km) CONSUMO (kWh) (INICIO -ESTAC) CONSUMO (kWh) (ESTAC -ESTAC) CONSUMO (kWh) (ESTAC -ESTAC) ELIMINACION NICOYA 61.80 10.90 10.90 TAMARINDO 117.32 19.67 8.77 18.77

RUTA 1 VUELTA TAMARINDO-NICOYA-LIMONAL ESTACIONES DISTANCIA (km) CONSUMO (kWh) (INICIO -ESTAC) CONSUMO (kWh) (ESTAC -ESTAC) CONSUMO (kWh) (ESTAC -ESTAC) ELIMINACION NICOYA 55.52 10.19 10.19 LIMONAL 117.32 20.37 10.18

Figura 35: C´alculosde la ruta 1 Limonal a Tamarindo escenario 1.

RUTA 2 IDA LIMONAL-NICOYA-HUACAS-TAMARINDO ESTACIONES DISTANCIA (km) CONSUMO (kWh) (INICIO -ESTAC) CONSUMO (kWh) (ESTAC -ESTAC) CONSUMO (kWh) (ESTAC -ESTAC) ELIMINACION NICOYA 61.80 10.90 5.82 5.82 HUACAS 124.56 21.54 10.63 TAMARINDO 137.12 23.41 1.87 12.51

RUTA 2 VUELTA TAMARINDO-HUACAS-NICOYA-LIMONAL ESTACIONES DISTANCIA (km) CONSUMO (kWh) (INICIO -ESTAC) CONSUMO (kWh) (ESTAC -ESTAC) CONSUMO (kWh) (ESTAC -ESTAC) ELIMINACION HUACAS 12.55 2.56 2.56 NICOYA 75.32 13.66 11.10 13.66 LIMONAL 137.12 23.99 10.33 10.33

Figura 36: C´alculosde la ruta 2 Limonal a Tamarindo escenario 1.

3.4. Rutas hacia norte del pa´ıs

Uno de los puntos m´aslejanos en el norte del pa´ıses Los Chiles, al cual se puede llegar por distintas rutas. La primer ruta para llegar a Los Chiles es llegando a San Ram´ony llega hasta Santa Rosa de Pocosol, dicho trayecto es de 86 km, pero esta ruta no es por carretera nacional, por lo que la distancia entre estas estaciones respeta la ley, ya que esta indica que las estaciones en ruta cantonal deben de estar a un m´aximode 120 km. Tambi´ense puede llegar de San Jos´e

2018 34 Propuesta de ubicaci´onde infraestructura de recarga r´apida EPERLab - UCR hasta Los Chiles por la ruta de Zarcero y por esta ruta se cuenta con mayor cantidad de centros de recarga a disposici´ony por tanto con una menor distancia entre las mismas.

Ruta 1 SJ - Los Chiles Pto Viejo - Santa Rosa Ruta 2 SJ - Los Chiles San Miguel - San Ram´on- Los Chiles Ruta 3 Limonal - Los Chiles Vara Blanca - Los Chiles

Si se deseara ir desde la frontera de Pe˜nasBlancas hasta Los Chiles no se podr´ıarealizar el viaje de manera segura, debido a que la distancia entre los puntos supera por m´asdel doble a la distancia m´aximaque deben tener los centros de recarga, siendo necesario que se coloque uno cerca de Upala para disminuir las distancias.

RUTA 1 IDA SJ - PTO VIEJO SARAPIQUI- SANTA ROSA - LOS CHILES ESTACIONES DISTANCIA (km) CONSUMO (kWh) (INICIO -ESTAC) CONSUMO (kWh) (ESTAC -ESTAC) PTO VIEJO SARAPIQUI 85.30 12.76 12.76 SANTA ROSA 158.20 25.16 12.40 LOS CHILES 210.82 33.99 8.83

RUTA 1 VUELTA LOS CHILES -SANTA ROSA - PTO VIEJO SARAPIQUI- SJ ESTACIONES DISTANCIA (km) CONSUMO (kWh) (INICIO -ESTAC) CONSUMO (kWh) (ESTAC -ESTAC) SANTA ROSA 52.63 9.11 9.11 PTO VIEJO SARAPIQUI 125.53 21.67 12.56 SAN JOSE 210.82 41.69 20.02

Figura 37: C´alculosde la ruta 1 San Jos´ea Los Chiles escenario 1.

Los datos de la primer ruta (solo un escenario) se pueden ver en la Figura 37, en los cua- les se muestra que el viaje de vuelta de Pto Viejo de Sarapiqu´ıa la capital se hace un gasto de 20.02 kWh, a pesar de estar dentro del m´aximoconsumo permitido, ser´ıarecomendable que en esta ruta se agregue un centro de recarga para que el VE no corra peligro de perder autonom´ıa.

La segunda ruta cuenta con dos escenarios, el primero mostrado en la Figura 39 indica que la estaci´onde CD Quesada no es conveniente eliminarla, esto porque en el viaje hacia San Jos´e se puede llegar a 20.09 kWh de consumo en uno de los trayectos de esta ruta, por lo que es mejor mantener dicha estaci´onpara evitar comprometer la autonom´ıadel VE.

El segundo escenario se puede ver en la Figura 40, donde se muestra que la estaci´onde Zarcero es necesaria para que el VE no se quede sin energ´ıaen medio trayecto, esto porque su consumo m´aximoes de 19.75 kWh, pero dicho consumo est´acercano a los 22.5 kWh permitidos, por lo que es mejor mantener el centro de recarga, para que la red de recarga sea un poco m´as robusta. Con los datos de esta ruta y con los mostrados en otras rutas se puede ver que la existir la estaci´onen Naranjo y en San Ram´on,esta ´ultimapuede ser eliminada y en los distintos viajes de l VE este no perder´aautonom´ıa.

2018 35 Propuesta de ubicaci´onde infraestructura de recarga r´apida EPERLab - UCR

N

164 km

N

164 km

N

53.5 km

50 km 16.5 km 27 km

50.2 km

86.1 km 31.2 km 10.1 km

Figura 38: Rutas San Jos´ehacia el norte del pa´ıs.

2018 36 Propuesta de ubicaci´onde infraestructura de recarga r´apida EPERLab - UCR

RUTA 2 IDA SJ - SAN MIGUEL-SAN RAMON- SANTA ROSA - LOS CHILES ESTACIONES DISTANCIA (km) CONSUMO (kWh) (INICIO -ESTAC) CONSUMO (kWh) (ESTAC -ESTAC) CONSUMO (kWh) (ESTAC -ESTAC) ELIMINACIÓN ZARCERO 64.8 14.37 14.37 14.37 CD QUESADA 92.30 16.58 2.21 SANTA ROSA 135.80 22.38 5.80 8.01 LOS CHILES 187.80 31.14 8.75 8.75

RUTA 2 VUELTA LOS CHILES-SANTA ROSA - SAN RAMON - SAN MIGUEL -SJ ESTACIONES DISTANCIA (km) CONSUMO (kWh) (INICIO -ESTAC) CONSUMO (kWh) (ESTAC -ESTAC) CONSUMO (kWh) (ESTAC -ESTAC) ELIMINACIÓN SANTA ROSA 52.00 9.01 9.01 9.01 CD QUESADA 95.50 19.31 10.30 ZARCERO 123.00 29.10 9.79 20.09 SAN JOSE 187.80 39.06 9.96 9.96

Figura 39: C´alculosde la ruta 2 San Jos´ea Los Chiles escenario 1.

RUTA 2 IDA SJ - SAN MIGUEL-SAN RAMON- SANTA ROSA - LOS CHILES ESTACIONES DISTANCIA (km) CONSUMO (kWh) (INICIO -ESTAC) CONSUMO (kWh) (ESTAC -ESTAC) CONSUMO (kWh) (ESTAC -ESTAC) ELIMINACIÓN ZARCERO 64.8 14.37 14.37 CD QUESADA 92.30 16.58 2.21 16.58 SANTA ROSA 135.80 22.38 5.80 5.80 LOS CHILES 187.80 31.14 8.75 8.75

RUTA 2 VUELTA LOS CHILES-SANTA ROSA - SAN RAMON - SAN MIGUEL -SJ ESTACIONES DISTANCIA (km) CONSUMO (kWh) (INICIO -ESTAC) CONSUMO (kWh) (ESTAC -ESTAC) CONSUMO (kWh) (ESTAC -ESTAC) ELIMINACIÓN SANTA ROSA 52.00 9.01 9.01 9.01 CD QUESADA 95.50 19.31 10.30 10.30 ZARCERO 123.00 29.10 9.79 SAN JOSE 187.80 39.06 9.96 19.75

Figura 40: C´alculosde la ruta 2 San Jos´ea Los Chiles escenario 2.

RUTA 3 IDA SJ - VARA BLANCA- SANTA ROSA - LOS CHILES ESTACIONES DISTANCIA (km) CONSUMO (kWh) (INICIO -ESTAC) CONSUMO (kWh) (ESTAC -ESTAC) VARA BLANCA 40.23 10.87 10.87 SANTA ROSA 127.30 20.80 9.93 LOS CHILES 178.64 29.20 8.41

RUTA 3 VUELTA LOS CHILES -SANTA ROSA -VARA BLANCA- SJ ESTACIONES DISTANCIA (km) CONSUMO (kWh) (INICIO -ESTAC) CONSUMO (kWh) (ESTAC -ESTAC) SANTA ROSA 51.34 9.08 9.08 VARA BLANCA 138.40 32.12 23.04 SAN JOSE 178.64 37.24 5.12

Figura 41: C´alculosde la ruta 3 San Jos´ea Los Chiles escenario 1.

En la tercer ruta cuyos datos est´anen la Figura 41 se puede ver que el consumo que se da de Vara Blanca hacia la capital sobrepasa los 22.5 kWh, por lo que es necesario que exista una estaci´onentre estas dos estaciones.

2018 37 Propuesta de ubicaci´onde infraestructura de recarga r´apida EPERLab - UCR 4. Recomendaciones finales para Ubicaci´on

A partir del an´alisisde la secci´on3, se deducen las siguientes ubicaciones para incluirse como requeridas para lograr una red: m´ınima,media o robusta. Esta red se logra a partir de los insumos recogidos en el taller. Asimismo, se requerir´ıala red de distribuci´onen Media Tensi´ondel ICE, la cual no est´adisponible para este an´alisis.Se utiliza la siguiente simbolog´ıa:

B Requiere atenci´onpues aunque es una ubicaci´oneliminada para obtener la red m´ınima, el consumo entre esta ubicaci´ony la anterior es muy cercano al l´ımitede 22.5 kWh. BB Requiere atenci´onpues deber´ıaa˜nadirseuna ubicaci´onantes de esta o mover la ubicaci´on en una posici´onanterior, ya que el consumo entre esta ubicaci´ony la anterior es muy cercano al l´ımitede 22.5 kWh. F La ubicaci´ones un destino tur´ıstico.La premisa es que los destinos tienen cargadores semi-r´apidos.Sin embargo, algunos viajes pueden requerir un retorno inmediato desde el destino: estas ubicaciones pueden tener cargadores r´apidos,preferiblemente. v La ubicaci´ones un destino, aunque no necesariamente es tur´ıstico.La situaci´onde an´alisis es similar que la de los destinos tur´ısticos.

Todas las redes de carga (m´ınima, media o robusta) tiene las ubicaciones de la GAM, ya que existen 9 cargadores r´apidosdisponibles para el ´areaservida de la CNFL y las cabeceras de provincia aparte de San Jos´e.Desde el punto de vista de autonom´ıa,la cantidad de cargadores de la GAM no afectar´ıala robustez de la red de cargadores. Entonces, las ubicaciones de la GAM son (ver secci´on 2.1 para revisar c´omose determinan):

Tabla 12: Ubicaciones de la GAM - CNFL.

Provincia Empresa Distribuidora Lugar de Referencia Referencia en Mapa 1 San Jos´e CNFL Circunvalaci´on(sur) E1 2 San Jos´e CNFL Curridabat (CFIA) E2 3 San Jos´e CNFL San Pedro (UCR) E3 4 San Jos´e CNFL Barrio Am´on E4 5 San Jos´e CNFL Tib´as E5 6 San Jos´e CNFL Sabana E6 7 San Jos´e CNFL La Uruca E7 8 San Jos´e CNFL Escaz´u E8 9 San Jos´e CNFL General Ca˜nas E9

Adem´asde las 9 estaciones de la CNFL, a continuaci´onse muestran otras estaciones que conforman cada una de las redes identificadas.

2018 38 Propuesta de ubicaci´onde infraestructura de recarga r´apida EPERLab - UCR

Tabla 13: Ubicaciones de la red de recarga m´ınima.

Provincia Empresa Distribuidora Lugar de Referencia Indicador 10 Alajuela CoopeAlfaroRuiz Zarcero - 11 Guanacaste CoopeGuanacaste Nicoya - 12 Puntarenas CoopeGuanacaste Paquera BB 13 Alajuela CoopeLesca La Fortuna - 14 Alajuela CoopeLesca Santa Rosa de Pocosol - 15 Heredia CoopeLesca Sarapiqu´ı - 16 Heredia ESPH Heredia (centro) GAM 17 Cartago JASEC Cartago (centro) GAM 18 Cartago CoopeSantos El empalme - 19 Guanacaste ICE Tilar´an - 20 Alajuela ICE San Miguel de Naranjo - 21 Alajuela ICE Alajuela* GAM 22 Guanacaste ICE Liberia - 23 Lim´on ICE Siquirres - 24 Lim´on ICE Lim´on - 25 Puntarenas ICE Jac´o - 26 Puntarenas ICE Quepos - 27 Puntarenas ICE Palmar Norte - 28 Puntarenas ICE Barranca - 29 Puntarenas ICE Limonal - 30 San Jos´e ICE P´erezZeled´on - 31 Heredia ICE** Vara Blanca BB 32 Lim´on ICE Gu´apiles B 33 San Jos´e ICE La Georgina B 34 Alajuela ICE Orotina B *Se recomienda colocar el cargador de Alajuela en el Aeropuerto por su ubicaci´onestrat´egica. ** Se debe indagar con exactitud cu´ales la empresa distribuidora que corresponde a Varablanca.

2018 39 Propuesta de ubicaci´onde infraestructura de recarga r´apida EPERLab - UCR

Tabla 14: Ubicaciones de la red de recarga media: a˜nadira la red m´ınima.

Provincia Empresa Distribuidora Lugar de Referencia Indicador 35 Cartago ICE Turrialba B 36 Alajuela ICE Los Chiles v 37 Lim´on ICE Sixaola v 38 Guanacaste CoopeGuanacaste Tamarindo F 39 Guanacaste CoopeGuanacaste Huacas F 40 Guanacaste ICE Pe˜nasBlancas F 41 Guanacaste ICE S´amara F 42 Lim´on ICE Puerto Viejo F 43 Puntarenas ICE Ciudad Neily F

Tabla 15: Ubicaciones de la red de recarga robusta: a˜nadira la red media.

Provincia Empresa Distribuidora Lugar de Referencia Indicador 44 Lim´on ICE Sixaola v 45 Cartago JASEC Guarco - 46 Alajuela CoopeLesca Ciudad Quesada - 47 Cartago CoopeSantos Macho Gaff - 48 Alajuela ICE San Ram´on -

2018 40 Propuesta de ubicaci´onde infraestructura de recarga r´apida EPERLab - UCR 5. Conclusiones

Se desarroll´ouna metodolog´ıade an´alisisde las ubicaciones de los centros de recarga r´apida en la GAM que permite comprender la cobertura de los cargadores r´apidosen las rutas principales que conectan los cantones ubicados en esta parte del pa´ıs. Queda pendiente robustecer el an´alisiscon datos de flujo de transporte del MOPT, los cuales deben tratarse con m´etodos estad´ısticose integrarse a los sistemas de informaci´ongeogr´afica. Se desarroll´ouna metodolog´ıade an´alisisde los centros de recarga r´apidaque permiten ase- gurar la autonom´ıade los veh´ıculosel´ectricos.Se identifica que las ubicaciones se˜naladasen los talleres permiten una cobertura adecuada de las rutas de mayor distancia considerando el efecto de la elevaci´onen la autonom´ıa. Se requieren 34 centros de recarga r´apidaque asegure la suficiente autonom´ıade los veh´ıculos el´ectricospara llegar a los puntos m´aslejanos del pa´ıs. De estos 34, 9 se ubican en San Jos´ey algunos de ellos est´anpor la CNFL: estos cargadores responden a la demanda de transporte en la GAM m´asque al cumplimiento de autonom´ıa.Esta red de recarga corresponde a una infraestructura m´ınimade centros de recarga r´apida.La ubicaci´onde estas estaciones de determina con insumos de los talleres y de un an´alisisrealizado por el EPERLab. La red de infraestructura m´ınimade recarga r´apidapuede mejorarse con 9 cargadores para lograr desarrollar una infraestructura con robustez media que permitir´ıa a los usuarios cargar en los destinos alejados. Esta condici´onmejorar´ıala flexibilidad de los usuarios para regresar de sus destinos r´apidamente. La ubicaci´onde estas estaciones de determina con insumos de los talleres y de un an´alisisrealizado por el EPERLab. La red de infraestructura media de recarga r´apidapuede mejorarse con 5 cargadores, resul- tando en la conformaci´onde una red de recarga robusta. La ubicaci´onde estas estaciones de determina con insumos de los talleres y de un an´alisisrealizado por el EPERLab. El an´alisisde infraestructura de recarga r´apidadebe actualizarse cuando la penetraci´onde los veh´ıculosel´ectricosaumente. Adem´as,la informaci´ondeber´acomplementarse con infor- maci´onde tr´aficovehicular que deber´arecopilarse de forma sistematizada en coordinaci´on con los entes sectoriales. Una cobertura robusta de cargadores en la GAM puede facilitar la adopci´onde taxis el´ectri- cos en el sistema de transporte. No obstante, este cambio tecnol´ogicodepende de otros aspectos t´ecnicosy econ´omicospropios del gremio.

2018 41 Propuesta de ubicaci´onde infraestructura de recarga r´apida EPERLab - UCR Referencias

[1] Sitio web de la Agencia de Protecci´onAmbiental de EEUU (EPA). Emission Standards Reference Guide: EPA Urban Dynamometer Driving Schedule (UDDS). Recuperado el 10 de junio del 2018 de: https://www.epa.gov/emission-standards-reference-guide/ epa-urban-dynamometer-driving-schedule-udds

[2] Sitio web de la Agencia de Protecci´onAmbiental de EEUU (EPA). Vehicle and Fuel Emissions Testing: Dynamometer Drive Schedules. Recuperado el 10 de junio del 2018 de: https:// www.epa.gov/vehicle-and-fuel-emissions-testing/dynamometer-drive-schedules

[3] R. Maia, M. Silva, R. Ara´ujoand U. Nunes (2011). Electric vehicle simulator for energy consumption studies in electric mobility systems. IEEE Forum on Integrated and Sustainable Transportation Systems, Vienna, pp. 227-232. doi: 10.1109/FISTS.2011.5973655 http:// ieeexplore.ieee.org/stamp/stamp.jsp?tp=&arnumber=5973655&isnumber=5973576

[4] Sitio web evrater. The Electric Vehicle List. Recuperado el 9 de Julio del 2018 de https: //evrater.com/evs

2018 42