2025-11-14T17:49:11.423932

Bilevel Optimization for Improved Flexibility Aggregation Models of Electric Vehicle Fleets

Härtel, von Bonin
Electric vehicle (EV) fleets are expected to become an increasingly important source of flexibility for power system operations. However, accurately capturing the flexibility potential of numerous and heterogeneous EVs remains a significant challenge. We propose a bilevel optimization formulation to enhance flexibility aggregations of electric vehicle fleets. The outer level minimizes scheduling deviations between the aggregated and reference EV units, while the inner level maximizes the aggregated unit's profits. Our approach introduces hourly to daily scaling factor mappings to parameterize the aggregated EV units. Compared to simple aggregation methods, the proposed framework reduces the root-mean-square error of charging power by 78~per cent, providing more accurate flexibility representations. The proposed framework also provides a foundation for several potential extensions in future work.
academic

Optimización Bilineal para Modelos Mejorados de Agregación de Flexibilidad de Flotas de Vehículos Eléctricos

Información Básica

  • ID del Artículo: 2506.04843
  • Título: Optimización Bilineal para Modelos Mejorados de Agregación de Flexibilidad de Flotas de Vehículos Eléctricos
  • Autores: Philipp Härtel (Fraunhofer IEE), Michael von Bonin (Fraunhofer IEE)
  • Clasificación: math.OC cs.SY eess.SY
  • Fecha de Publicación: 5 de junio de 2025 (preimpresión en arXiv)
  • Enlace del Artículo: https://arxiv.org/abs/2506.04843

Resumen

Se espera que las flotas de vehículos eléctricos (VE) se conviertan en un recurso de flexibilidad cada vez más importante en la operación de sistemas eléctricos. Sin embargo, capturar con precisión el potencial de flexibilidad de numerosos vehículos eléctricos heterogéneos sigue siendo un desafío significativo. Este artículo propone una formulación de optimización bilineal para mejorar la agregación de flexibilidad de flotas de VE. El nivel superior minimiza la desviación de programación entre la unidad VE agregada y la unidad VE de referencia, mientras que el nivel inferior maximiza la ganancia de la unidad agregada. El método introduce un mapeo de factores de escala de hora a día para parametrizar la unidad VE agregada. En comparación con métodos de agregación simple, el marco propuesto reduce el error cuadrático medio de la potencia de carga en un 78%, proporcionando una representación más precisa de la flexibilidad.

Antecedentes y Motivación de la Investigación

Definición del Problema

Con la adopción generalizada de vehículos eléctricos, las flotas de VE se convertirán en un recurso de flexibilidad importante para los sistemas eléctricos. Sin embargo, modelar con precisión la flexibilidad de carga y descarga de flotas de VE heterogéneos a gran escala presenta desafíos significativos:

  1. Complejidad Computacional: Modelar individualmente cada VE en la planificación de sistemas eléctricos a gran escala es computacionalmente inviable
  2. Desafíos de Heterogeneidad: Diferentes VE poseen características distintas como patrones de conducción, capacidad de batería, comportamiento de carga, etc.
  3. Sobrestimación de Flexibilidad: Los métodos de agregación simple existentes tienden a sobrestimar el potencial real de flexibilidad de la flota

Importancia de la Investigación

Las predicciones del mercado eléctrico alemán para 2045 muestran que la demanda de energía de VE mantiene una mediana inferior a 10 GW durante todo el año, pero el cuartil superior frecuentemente supera 40 GW, exhibiendo una variabilidad enorme. La introducción de capacidades de vehículo a red (V2G) amplifica aún más la necesidad de una representación precisa de la flexibilidad.

Limitaciones de Métodos Existentes

  1. Métodos de Agregación Ingenua: Estiman la flexibilidad de la flota mediante suma simple de disponibilidad de carga y demanda de conducción de VE individuales, ignorando interacciones y dependencias entre vehículos
  2. Agregación de Batería Virtual: Modelan la flota de VE como un sistema de almacenamiento de energía virtual, pero aún dependen de factores de escala genéricos ajustados manualmente
  3. Técnicas de Agrupamiento: Agrupan VE con patrones de carga similares, pero carecen de estrategias de agregación optimizadas

Contribuciones Principales

  1. Propone Marco de Optimización Bilineal: Aplica por primera vez la optimización bilineal a la agregación de flexibilidad de flotas de VE, optimizando factores de escala en el nivel superior para minimizar desviación de programación y optimizando ganancias de participación en el mercado de la unidad agregada en el nivel inferior
  2. Diseña Mecanismo de Mapeo de Factores de Escala: Introduce mapeo de factores de escala de escala temporal de hora a día, capaz de capturar características de flexibilidad en diferentes granularidades temporales
  3. Mejora Significativa de Precisión de Agregación: Logra una reducción del 78% en RMSE de potencia de carga en comparación con métodos de agregación simple, proporcionando una representación más precisa de la flexibilidad
  4. Proporciona Marco Escalable: Establece infraestructura para futuras extensiones (como V2G, diferentes tipos de usuarios, carga rápida, etc.)

Explicación Detallada del Método

Definición de la Tarea

Entrada: Disponibilidad de carga/descarga de unidades VE individuales, trayectorias de SOC, precios de electricidad, demandas de conducción Salida: Factores de escala optimizados y estrategia de programación de la unidad VE agregada (AEV) Objetivo: Minimizar la desviación de programación entre la unidad AEV y la unidad VE de referencia individual

Arquitectura de Optimización Bilineal

Problema de Nivel Superior (Nivel Externo)

Minimiza la desviación de programación entre la unidad agregada y la unidad VE de referencia:

minΞOLtT[γC(xu,tCX^Vu,tC)2+γD(xu,tDX^Vu,tD)2+γS(xu,tSX^Vu,tS)2]\min_{\Xi^{OL}} \sum_{t \in T} \left[ \gamma^C (x^C_{u,t} - \hat{X}^C_{V_u,t})^2 + \gamma^D (x^D_{u,t} - \hat{X}^D_{V_u,t})^2 + \gamma^S (x^S_{u,t} - \hat{X}^S_{V_u,t})^2 \right]

Las restricciones incluyen:

  • Mapeo de factores de escala: xu,tC=κu,τCvVuXv,tCx^C_{u,t} = \kappa^C_{u,\tau} \sum_{v \in V_u} X^C_{v,t}, donde τ=fnC(t)\tau = f^C_n(t)
  • No negatividad de factores de escala: κu,τC,κu,τCR+\kappa^C_{u,\tau}, \overline{\kappa}^C_{u,\tau} \in \mathbb{R}_+

Problema de Nivel Inferior (Nivel Interno)

Maximiza la ganancia de participación en el mercado de la unidad agregada:

minΞILtTΠt(xu,tCxu,tD)\min_{\Xi^{IL}} \sum_{t \in T} \Pi_t (x^C_{u,t} - x^D_{u,t})

Restricciones:

  • Límites de potencia de carga/descarga: xu,tCxu,tCxu,tC\underline{x}^C_{u,t} \leq x^C_{u,t} \leq \overline{x}^C_{u,t}
  • Continuidad de SOC: xu,t+1S=ρuxu,tSΦu,tDRΦu,tTH+ηuCxu,tC1ηuDxu,tDx^S_{u,t+1} = \rho_u x^S_{u,t} - \Phi^{DR}_{u,t} - \Phi^{TH}_{u,t} + \eta^C_u x^C_{u,t} - \frac{1}{\eta^D_u} x^D_{u,t}

Mecanismo de Mapeo de Factores de Escala

Se diseña una función de mapeo periódica: fnC(t)=fnD(t)=fnS(t)=tmod(7×24)nf^C_n(t) = f^D_n(t) = f^S_n(t) = \left\lfloor \frac{t \bmod (7 \times 24)}{n} \right\rfloor

Este mapeo asegura un ciclo semanal (168 horas), mientras agrupa cada n horas dentro del ciclo, procesando cada día laboral.

Reconstrucción de Condiciones KKT

Debido a la convexidad del problema de nivel inferior, se reconstruye el problema bilineal como un programa lineal entero mixto de un solo nivel:

  • Introduce variables duales λu,t,μu,t()\lambda_{u,t}, \mu^{(\cdot)}_{u,t}
  • Añade condiciones de optimalidad KKT
  • Reconstruye condiciones de holgura complementaria usando el método de gran M

Configuración Experimental

Conjunto de Datos

  • Rango Temporal: Período de 3 semanas, resolución horaria (datos alemanes de enero de 2012)
  • Configuración de VE: 900 perfiles de conductores commuters, basados en datos de encuestas de viajes históricos
  • Precios de Electricidad: Precios eléctricos alemanes predichos para 2035
  • Especificaciones de Vehículos: Incluye capacidad de batería, eficiencia de conducción, potencia de carga/descarga doméstica, límites de SOC

Métricas de Evaluación

  • Métrica Principal: Error Cuadrático Medio (RMSE) de potencia de carga
  • Métricas de Optimización: Valor de función objetivo, brecha de optimalidad
  • Dimensiones de Comparación: Trayectorias de SOC, desviación de programación de potencia de carga

Métodos de Comparación

  • Agregación Simple (SA): Método heurístico de agregación basado en almacenamiento de energía virtual
  • Variantes de AEV: Diferentes mapeos de factores de escala (1h, 2h, 4h, 6h, 24h)

Detalles de Implementación

  • Solucionador: Gurobi Optimizer
  • Configuración Simplificada: No considera capacidades V2G actualmente, enfocándose en verificación del método principal
  • Tipos de Usuarios: Solo se consideran commuters, reduciendo complejidad de heterogeneidad

Resultados Experimentales

Resultados Principales

MétodoValor Objetivo ÓptimoCota Inferior ÓptimaBrecha RelativaRMSE (MW)
AEV 1h2891.691957.9532.29%2.926
AEV 2h3734.222921.8921.75%3.125
AEV 4h4130.893709.8910.19%3.164
AEV 6h4300.004162.433.20%3.256
AEV 24h5666.595609.961.00%3.827
SA---18.084

Hallazgos Clave:

  • El método AEV logra una reducción del 78% en RMSE en comparación con SA
  • Los mapeos de granularidad más fina (1h, 2h) producen desviaciones de programación más pequeñas
  • El mapeo de 24h alcanza una brecha de optimalidad del 1%

Análisis de Factores de Escala

  • La disponibilidad de potencia de carga máxima agregada se reduce entre 0%-18%
  • La trayectoria máxima de SOC se reduce hasta el 60%
  • El rango de trayectoria mínima de SOC es 0%-15% del SOC máximo agregado
  • Existen patrones de variación entre diferentes días laborales

Análisis de Comportamiento de Programación

De la Figura 4 se puede observar:

  • La programación de la unidad AEV es altamente consistente con la programación de referencia de VE individuales
  • El método SA produce picos de carga más pronunciados durante períodos de bajo precio
  • El método AEV captura mejor las restricciones reales de flexibilidad

Trabajo Relacionado

Métodos de Agregación de VE

  1. Agregación de Batería Virtual: Modelan la flota de VE como un sistema de almacenamiento de energía virtual, considerando límites de carga/descarga y capacidad de batería
  2. Técnicas de Agrupamiento: Agrupan VE con patrones de carga similares, evitando mezcla de vehículos heterogéneos
  3. Modelos Equivalentes: Adoptan conceptos de modelos compuestos o equivalentes de sistemas hidroeléctricos

Aplicaciones de Optimización Bilineal

  • Identificación de modelos simplificados precisos en sistemas hidroeléctricos
  • Despeje de mercados eléctricos y planificación de capacidad
  • Respuesta a demanda y optimización de almacenamiento

Conclusiones y Discusión

Conclusiones Principales

  1. Efectividad del Método: El marco de optimización bilineal mejora significativamente la precisión de agregación de flexibilidad de flotas de VE
  2. Viabilidad Computacional: Se logra solución eficiente mediante reconstrucción KKT
  3. Valor Práctico: Proporciona una representación precisa de flexibilidad de VE para planificación de sistemas eléctricos a gran escala

Limitaciones

  1. Supuestos Simplificados: La versión actual no considera V2G, carga en lugares de trabajo, carga rápida y otros escenarios
  2. Heterogeneidad de Usuarios: Solo considera commuters, sin cubrir diferentes tipos de usuarios
  3. Complejidad Computacional: Algunas instancias no alcanzan optimalidad global (brecha relativa >1%)

Direcciones Futuras

  1. Extensiones Funcionales: Integrar capacidades V2G, múltiples tipos de usuarios, infraestructura de carga rápida
  2. Análisis de Sensibilidad: Diferentes escenarios de precios, factores de escala específicos por día laboral
  3. Rol del Agregador: Utilizar estructura bilineal para analizar mecanismos de incentivos de agregadores de flexibilidad
  4. Aplicaciones en Tiempo Real: Extender a programación en tiempo real y participación en mercados

Evaluación Profunda

Fortalezas

  1. Innovación Metodológica: Primera aplicación sistemática de optimización bilineal a agregación de flotas de VE, con base teórica sólida
  2. Suficiencia Experimental: Comparación de múltiples escalas temporales de mapeo, resultados convincentes
  3. Alto Valor Práctico: La mejora de precisión del 78% es significativa para planificación de sistemas eléctricos
  4. Marco Extensible: Proporciona buena infraestructura para investigación futura

Deficiencias

  1. Limitaciones de Escenarios: Los escenarios experimentales actuales son relativamente simplificados, careciendo de verificación de funciones importantes como V2G
  2. Eficiencia Computacional: Algunas instancias presentan brechas de solución mayores, que podrían afectar aplicaciones prácticas
  3. Sensibilidad de Parámetros: Falta análisis de sensibilidad para parámetros clave (como pesos de desviación γ)
  4. Verificación de Escalabilidad: No verifica suficientemente el desempeño del método en flotas de mayor escala

Impacto

  1. Contribución Académica: Proporciona nuevo paradigma de optimización para el campo de agregación de VE
  2. Valor de Ingeniería: Puede aplicarse directamente en herramientas de planificación de sistemas eléctricos
  3. Apoyo Político: Proporciona base técnica para formular mecanismos de mercado de flexibilidad de VE

Escenarios Aplicables

  1. Planificación de Sistemas Eléctricos: Evaluación de flexibilidad de VE en decisiones de inversión a largo plazo
  2. Diseño de Mercados: Diseño de mecanismos para participación de agregadores de VE en mercados eléctricos
  3. Optimización Operativa: Control coordinado de flotas de VE en redes de distribución

Referencias

El artículo cita 16 referencias relacionadas, cubriendo principalmente:

  • Métodos de agregación de VE: agregación de almacenamiento virtual, técnicas de agrupamiento
  • Teoría de optimización bilineal: condiciones KKT, reconstrucción de enteros mixtos
  • Modelado equivalente de sistemas hidroeléctricos: modelos compuestos, aplicaciones bilineales
  • Datos de viajes alemanes: encuesta MiD, estadísticas de transporte

Evaluación General: Este es un artículo de alta calidad con contribuciones importantes en el campo de agregación de flexibilidad de flotas de VE. El método es novedoso, los experimentos son exhaustivos, los resultados son significativos y proporciona una solución efectiva para abordar desafíos tecnológicos clave en la integración de VE a gran escala. Aunque existen algunas limitaciones, su valor innovador y potencial práctico lo convierten en un progreso importante en este campo.