scieee AI-readable full text Open interactive document viewer

Análisis y optimización de la oferta de transporte aéreo en España

Castro Díaz, Esteban

Abstract

En este trabajo se aborda el análisis de la desigualdad en el tráfico de pasajeros entre los diversos aeropuertos españoles, así como la distribución de estos pasajeros en las diferentes rutas ofrecidas por cada aeropuerto. Se considerará tanto el escenario nacional como internacional de pasajeros. Este análisis se amplía al comparar estos resultados con datos de aeropuertos de otros países como Inglaterra, Francia, Alemania y Estados Unidos. Posteriormente, se desarrolla y explica un modelo creado utilizando MATLAB® con el objetivo de optimizar el número de vuelos dentro de una red de aeropuertos españoles (tráfico nacional), con el propósito de minimizar el coste de operación en toda la red. En este modelo se consideran diferentes tipos de aeronaves y se aplican diversas restricciones con el fin de obtener resultados que tengan un sentido físico, coherente y práctico. Los datos de pasajeros y operaciones estarán basados en datos reales disponibles en AENA Estadísticas. Posteriormente se analizarán los resultados obtenidos para diferentes escenarios.

Full text

Equation Chapter 1 Section 1 Trabajo Fin de Grado Ingeniería Aeroespacial Análisis y optimización de la oferta de transporte aéreo en España Autor: Esteban Castro Díaz Tutor: José María del Castillo Granados Dpto. Ingeniería y Ciencia de los Materiales y del Transporte Escuela Técnica Superior de Ingeniería Universidad de Sevilla Sevilla, 2024 iii Trabajo Fin de Grado Ingeniería Aeroespacial Análisis y optimización de la oferta de transporte aéreo en España Autor: Esteban Castro Díaz Tutor: José María del Castillo Granados Catedrático de Universidad Dpto. Ingeniería y Ciencia de los Materiales y del Transporte Escuela Técnica Superior de Ingeniería Universidad de Sevilla Sevilla, 2024 v Trabajo Fin de Grado: Análisis y optimización de la oferta de transporte aéreo en España Autor: Esteban Castro Díaz Tutor: José María del Castillo Granados El tribunal nombrado para juzgar el Proyecto arriba indicado, compuesto por los siguientes miembros: Presidente: Vocales: Secretario: Acuerdan otorgarle la calificación de: Sevilla, 2024 El Secretario del Tribunal vii ix Agradecimientos A mis padres José Manuel y Esperanza por confiar siempre en mí, por su sacrificio, por darme la oportunidad de cumplir todos mis objetivos y pagarme la broma. A mi hermana Cecilia por aguantarme. A Blanca por estar a mi lado siempre, confiar en mí y quererme. A mi abuelo Manolo, que, aunque nunca entendiera que estaba haciendo, sé que estaría orgulloso de todos mis logros. A todos los que me han apoyado. Marchena, 2024 ÍNDICE DE TABLAS Tabla 1 Aeropuertos y población de diferentes países de Europa 10 Tabla 2 Parámetros del ajuste de las figuras 5 y 6 empleando el método de mínimos cuadrados ordinarios (MCO) y Theil-Sen 15 Tabla 3 Valores estadísticos más importantes de los ajustes realizados en la figura 10 19 Tabla 4 Porcentaje de pasajeros de los principales países de la Unión Europea 20 Tabla 5 Coeficiente de Gini del reparto de pasajeros entre los principales 27 aeropuertos de cada país 21 Tabla 6 Quince mayores aeropuertos por volumen de tráfico comercial de España, Estados Unidos y Reino Unido 23 Tabla 7 Número de pasajeros transportados por las cinco operadoras comerciales más importantes en España por volumen de tráfico y casos de Binter Canarias y Air Nostrum. [28] 32 Tabla 8 Flota conjunta de todas las aerolíneas contempladas en la tabla 5 [30] 32 Tabla 9 Principales características de las aeronaves seleccionadas [31] [32] 33 Tabla 10 Total de pasajeros y operaciones y media de pasajeros por operación de los 24 aeropuertos de la red. 35 Tabla 11 Tiempos de estacionamiento en el aeropuerto de las aeronaves empleadas 37 Tabla 12 Precios por aeronave de servicios adicionales contemplados 38 Tabla 13 Comparación de costes de operación para la aeronave Boeing 737-300 de diferentes operadoras estadounidenses (2003).[35] 40 Tabla 14 Estimación de consumo de combustible para vuelos de unos 930 Km [39] [40] 42 Tabla 15 CASM operativos asignados a las dos aeronaves empleadas 42 Tabla 16 Comparación del porcentaje de uso del Boeing 737-800 en para el modelo sin restricciones y con restricciones en algunas rutas en las cuales se empleaba mayoritariamente el Embraer 190 55 Tabla 17 Principales especificaciones necesarias dentro del modelo para las aeronaves escogidas A319 Neo y A321 Neo [43] [44] 70 Tabla 18 Eficiencia de combustible de todas las aeronaves empleadas [45] [46] 71 Tabla 19 Tiempos de ejecución del modelo en función de la variación relativa de CASM de las dos aeronaves empleadas sin restricciones 80 xvii ÍNDICE DE FIGURAS Figura 1 Recta de perfecta igualdad y curva Lorenz que determinan las áreas A y B. 8 Figura 2 Aeropuertos comerciales en España. En verde, aquellos de titularidad estatal operados por Aena. En naranja, los de titularidad autonómica. [13] 11 Figura 3: Evolución del número de pasajeros (Internacionales y nacionales) en la última década (datos numéricos en millones de pasajeros). 12 Figura 4: Fracción acumulada de pasajeros sobre el total en orden creciente del conjunto de todos los aeropuertos españoles. 13 Figura 5 : Coeficiente de Gini para los mayores 25 aeropuertos de España en volumen de tráfico nacional e internacional en función del número total de pasajeros de dicho aeropuerto 14 Figura 6 Coeficiente de Gini para los mayores 25 aeropuertos de España en volumen de tráfico exclusivamente nacional en función del número total de pasajeros de dicho aeropuerto 15 Figura 7 Coeficientes de Gini calculados para cada aeropuerto empleando tráfico nacional e internacional (azul) y exclusivamente nacional (rojo) 16 Figura 8: Comparación de los coeficientes de Gini reales y los calculados por el método de máxima entropía. 17 Figura 9 : Comparación del sumatorio de pasajeros normalizado tanto con datos reales como los calculados por el método de máxima entropía para los aeropuertos de Barcelona-El Prat, Palma de Mallorca, Ibiza y Sevilla. 18 Figura 10 Relación entre pasajeros reales y de máxima entropía (arriba) e influencia de la distancia en la diferencia de dichos valores (abajo) 19 Figura 11: Fracción acumulada de pasajeros sobre el total en orden creciente del conjunto de 27 mayores aeropuertos de Reino Unido, España, Alemania, Francia y Estados Unidos 22 Figura 12: Coeficientes de Gini ordenados de menor a mayor del reparto de los pasajeros entre los principales 10 destinos para cada uno de los 15 mayores aeropuertos de Estados Unidos (azul), España (rojo), Reino Unido (amarillo). Aeropuertos indicados en el eje X (España-Estados Unidos-Reino Unido) 24 Figura 13 Aeropuertos españoles empleados en el modelo. 31 Figura 14 Boeing 737-800 (RyanAir) 33 Figura 15 Embraer 190 (Binter) 33 Figura 16: Ejemplo de tarifas aplicadas por AENA Aeropuertos (enero 2024) por el aterrizaje y tránsito del aeródromo para los diferentes aeropuertos de España. 39 Figura 17: Determinación aproximada del coste por asiento disponible milla para diferentes modelos (1999) [38] 40 Figura 18: Diferentes costos por categoría para aeronaves de diferente tamaño. (Solo se consideran combustible, mantenimiento y tripulación). [25] 41 Figura 19: Número de rutas en función de su número de pasajeros en el mes de agosto de 2022 47 Figura 20: Porcentaje de vuelos realizado por la aeronave grande en cada ruta en función de la distancia 48 Figura 21: Puntos a partir del cual la aeronave grande pasa a ser más rentable de operar en función de la distancia para cada aeropuerto según sus tarifas aeroportuarias. 49 Figura 22: Variación del porcentaje de empleo de la aeronave grande en función de la distancia de la ruta para diferentes valores relativos del CASM de cada aeronave. 50 Figura 23: Resultados obtenidos de los vuelos necesarios a realizar con un modelo sin limitaciones de operación para el aeropuerto de Madrid-Barajas. 51 xix Figura 24 Comparación de los resultados obtenidos empleando (derecha) o no (izquierdas) las restricciones de operación y uso de aeronaves por distancia de la ruta 54 Figura 25 Porcentaje de vuelos realizado en cada ruta con la aeronave Boeing 737-800 en función de la distancia de la ruta. 56 Figura 26: Porcentaje de vuelos realizado en cada ruta con la aeronave Boeing 737-800 en función de la demanda de la ruta. 57 Figura 27 Resultados para los aeropuertos de Madrid-Barajas y Plama de Mallorca 58 Figura 28 Resultados obtenidos para los aeropuertos de Barcelona-El Prat y Gran Canaria 59 Figura 29 Resultados obtenidos para los aeropuertos de Tenerife Norte e Ibiza. 59 Figura 30 Resultados obtenidos para los aeropuertos Bilbao y Lanzarote-César Manrique 60 Figura 31 Resultados obtenidos para los aeropuertos de Vigo y A Coruña 61 Figura 32 Resultados obtenidos para los aeropuertos de Málaga-Costa del Sol y Santiago-Rosalía de Castro 61 Figura 33 Resultados obtenidos para los aeropuertos de Sevilla y Fuerteventura 62 Figura 34 Resultados obtenidos para los aeropuertos de Tenerife Sur y Menorca 63 Figura 35 Resultados obtenidos para los aeropuertos de Asturias y Valencia 63 Figura 36 Resultados obtenidos para los aeropuertos de Alicante-Elche Miguel Hernández y San Sebastián 64 Figura 37 Resultaos obtenidos para los aeropuertos de Santander-Severiano Ballesteros y Granada-Jaén 64 Figura 38 Resultados obtenidos para los aeropuertos de Zaragoza y Jerez de la Frontera-La Parra 65 Figura 39 Frecuencias diarias necesarias para cubrir la demanda dada en aquellas rutas que permitan tener al menos un vuelo diario. 67 Figura 40 Número de frecuencias semanales necesarias para cubrir la demanda dada para aquellas rutas que solo requieren varias conexiones semanales 67 Figura 41 Análisis de los precios por asiento de ambas aeronaves empleadas para todos vuelos a los aeropuertos de la red con origen en el aeropuerto de Madrid-Barajas 69 Figura 42 Análisis de los precios por asiento de ambas aeronaves empleadas para todos vuelos a los aeropuertos de la red con origen en los aeropuertos de Barcelona-El Prat (Izquierda) y Tenerife Norte (Derecha) 69 Figura 43 Airbus A321 Neo 70 Figura 44 Airbus A319 Neo 70 Figura 45 Porcentaje de uso de la aeronave grande (A321-Neo) en función de la distancia de la ruta en kilómetros (arriba) y en función de la demanda de la ruta (abajo) 72 Figura 46 Resultados para los aeropuertos de Madrid-Barajas, Palma de Mallorca, Barcelona-El Prat, Gran Canaria, Tenerife Norte e Ibiza 73 Figura 47 Resultados para los aeropuertos Bilbao, Lanzarote-César Manrique, Vigo, Coruña, Málaga-Costa del Sol y Santiago-Rosalía de Castro 74 Figura 48 Resultados para los aeropuertos Sevilla, Fuerteventura, Tenerife Sur, Menorca, Asturias, Valencia 75 Figura 49 Resultados para los aeropuertos Alicante-Elche, San Sebastián, Santander-Severiano Ballesteros, Granada-Jaén, Zaragoza y Jerez de la Frontera 76 Figura 50 Porcentaje de uso de cada aeronave y porcentaje de pasajeros transportados por cada aeronave en los dos modelos simulados. A319 Y A321 (arriba), Embraer 190 y B737 (abajo) 77 Figura 51 Rango en el que cada aeronave es más barata de emplear por pasajero según los datos contemplados en el trabajo 78 Figura 52 : Creación de la matriz (A) y vector de desigualdad (b) que genera la restricción de capacidad. 93 Figura 53: Creación de la matriz de desigualdad (A) que modela la restricción de operaciones débil. 94 Figura 54 Creación de la matriz y vector de igualdad que genera la restricción de operaciones de todas las aeronaves que no sean la más pequeña en vuelos de corta duración. 95 1 1 INTRODUCCIÓN l principal objetivo de este trabajo es desarrollar un modelo de asignación de aeronaves para atender la demanda de pasajeros dentro de una red de aeropuertos españoles, utilizando MATLAB® y datos estadísticos reales sobre el número de pasajeros y operaciones entre aeropuertos obtenidos de AENA Estadísticas como datos para entender el buen funcionamiento del modelo. Este enfoque busca encontrar la asignación óptima de aeronaves para satisfacer la demanda en diferentes rutas dentro de la red. El problema de asignación de aeronaves es de gran interés debido a su relevancia en la eficiencia operativa y económica de las aerolíneas y el transporte aéreo en general. Al encontrar la asignación óptima de aeronaves, se pueden minimizar los costes operativos, optimizar el uso de recursos y mejorar la calidad del servicio para los pasajeros. Lo que se busca es, seleccionando una serie de aeropuertos nacionales, y empleando datos reales de demanda y operaciones en un periodo de tiempo determinado entre dichos aeropuertos. Obtener cuál de las aeronaves que pueden emplearse es la óptima de seleccionar buscando que las operaciones sean lo más baratas posibles dentro de una serie de restricciones que limitan las operaciones que se pueden realizar a datos de operaciones reales de cada aeropuerto y ruta dentro de la red. Por otro lado, el análisis del número de pasajeros entre aeropuertos y la informidad del reparto de dichos pasajeros entre las diferentes rutas y destinos de cada aeropuerto para posteriormente realizar una comparación de resultados entre aeropuertos nacionales y de otros países, proporciona una visión detallada del funcionamiento de la aviación comercial en el mercado nacional. Esta uniformidad en la distribución de datos entre destinos ayuda a identificar patrones de demanda y comportamientos en las rutas aéreas, lo que permite a las aerolíneas y autoridades de aviación comprender mejor las necesidades de los pasajeros y optimizar la asignación de recursos, como aeronaves y rutas. 1.1. Breve historia de la aviación y la aviación comercial Tras la hazaña de los Hermanos Wright en 1903, los primeros años de evolución de esta nueva invención son debido principalmente a los acaudalados bolsillos de los denominados “sportsman”, nobles y millonarios que al igual que competían con automóviles, motos o practicaban la navegación ahora tenían una nueva disciplina donde batir marcas y récords. La aviación primitiva era más un deporte que un medio de transporte. Poco a poco entre diferentes grupos militares de las diversas potencias europeas comienzan a contemplar la idea del empleo de la aviación como un nuevo medio bélico. Más que como una máquina de guerra, como un nuevo instrumento para llevar a cabo labores de reconocimiento. Es difícil encontrar intentos de utilizar este nuevo método de transporte de forma comercial durante sus primeros años de existencia, como si fue la creación de DELAG en 1909 por el Conde Zeppellin, que permitió por primera vez el transporte de pasajeros usando un medio aéreo haciendo uso de dirigibles. En total, entre 1909 y 1913, facilitó el movimiento de unos 34000 pasajeros.[1] No fue hasta el final de la primera guerra mundial que la aviación demostró que había llegado para quedarse, pero no solo en el ámbito militar sino también en lo civil. Tras el final de la Gran Guerra tuvieron lugar las primeras convenciones en las que se buscaba regular el uso de la aviación, destacando el convenio de París en 1919 el cual da lugar a la fundación de la Comisión Internacional para la Navegación Aérea (ICAN). A su vez, las primeras rutas comerciales operadas de forma regular aparecen siendo empleadas principalmente E Introducción 2 para el transporte. La ruta Viena-Kiev (verano 1918), es considerada como la primera ruta regular de correo. Inicialmente, el principal motivo de su uso fue el deterioro que había producido la guerra en las redes ferroviarias de Europa central. Tras la derrota y caída del Imperio Alemán se da lugar en la ciudad alemana de Weimar una asamblea nacional donde se proclamaría la nueva república alemana. Huelgas y revueltas generales en todo el país impiden el transporte seguro de políticos a la ciudad lo que supuso el momento perfecto para el establecimiento de la primera ruta de pasajeros regulares Berlín-Weimar, creándose para tal efecto la compañía Deutsche Luft Reederei (DLR) el 5 de febrero de 1919. Este hecho también supuso la primera ruta subvencionada de la historia. Diferentes países europeos siguieron los mismos pasos, creando otras rutas aéreas directas y sendas compañías para operarlas como pueden ser KLM en Países Bajos y DDL en Dinamarca. Es en los primeros años de la década de 1920 donde Francia se consolidó como la nación con mayor actividad en el transporte de pasajeros. Su red de rutas aéreas tejía una conexión sin precedentes, uniendo París con ciudades como Londres, Ámsterdam, Bucarest, Varsovia, Ginebra y Marsella.[2] Se puede considerar a 1919 año clave en la historia de la aviación comercial puesto que, además de hitos antes mostrados, debemos sumar la fundación de IATA como organismo que buscaba la cooperación entre explotadores y operadores para el avance del sector. Antes del estallido de la Segunda Guerra Mundial, IATA concentró sus esfuerzos en simplificar los procesos, unificar los documentos de transporte y estandarizar las normas técnicas. Su labor en la estandarización fue especialmente notable, logrando acuerdos en los procedimientos necesarios para la gestión del tráfico. El objetivo primordial era hacer del transporte aéreo una experiencia más cómoda y atractiva para los usuarios. IATA sentó las bases para un futuro más eficiente y accesible en el transporte aéreo, allanando el camino para el crecimiento exponencial que experimentaría la aviación en la posguerra. Por otro lado, en Estados Unidos, sería Aeromarina Airways la primera compañía aérea del país, empleando hidroaviones de forma exclusiva, más adelantes nacerían otras aerolíneas más conocidas en la actualidad como TWA (1925), PAN-American Airways (1927) o Delta Airlines (1928).[3] La hazaña del Espíritu de San Luis a manos del piloto Charles Lindberrgh cruzando por primera vez el Atlántico en un vuelo sin escalas (1927) y la apertura de la primera línea regular de vuelos de Zeppelines entre Alemania y Brasil (1932), siembran la idea de futuras líneas directas de pasajeros entre países a ambos lados del Atlántico, siendo culminada cuando en 1939 Pan-American inaugura la primera línea regular directa de pasajeros entre Marsella y Nueva York usando un Boeing 314.[4] El transcurso de la segunda guerra mundial es, sin duda, un punto de inflexión para la evolución de la aviación. Procedimientos, tecnologías como el radar y avances significativos en las capacidades de operación de las aeronaves permitieron al final de esta entrar en la que se pueda considerar como la edad de oro de la aviación. A su vez, el amplio despliegue de bases aéreas militares que dejaron de ser utilizadas para tal efecto en las diferentes potencias beligerantes da lugar a una nueva red de aeropuertos con pistas cimentadas. Esto permitió el uso cada vez más extendido de aeronaves con tren de aterrizaje de ruedas a expensas del abandono del uso de hidroaviones que habían supuesto gran parte de la flota comercial mundial previa al gran conflicto. Se debe destacar en este período la fundación de la Organización de Aviación Civil Internacional (OACI), sucesora de la ICAN, que promueve el uso de regulaciones únicas de la aviación a nivel mundial. Además, el gran desarrollo de aeronaves de transporte por parte de Estados Unidos permitió una rápida conversión de estos diseños a aeronaves para el transporte comercial, lo que produjo una tónica general en las siguientes décadas de superioridad de las empresas aeronáuticas americanas. Destacan en estos 3 3 Análisis y optimización de la oferta de transporte aéreo en España años aeronaves como el Douglas DC-4 (1942) y el Lockheed Constellation (1943). Al respecto, otras potencias europeas como Inglaterra, que habían centrado sus esfuerzos de diseño y producción durante la guerra en bombarderos, buscaron remediar este nuevo desequilibrio tecnológico con el desarrollo del primer avión comercial a reacción. Esfuerzos que culminaron en el malogrado de Havilland Comet (1949). Se inicia así una nueva era ‘’jet age’’ que permite llegar más rápido, más lejos y con mayor comodidad a nuevos destinos. Rápidamente aparecen nuevos competidores como el Boeing 707 (1957) o su contraparte soviética Túpolev Tu-104 (1955). El desarrollo y prosperidad económica de las décadas entre 1950-1970, permiten la aparición de los primeros vuelos comerciales a menor precio. Volar estaba dejando de ser un lujo y comenzaba a llegar a las masas de los países desarrollados. Si bien la entrada de las aeronaves a reacción había permitido reducir a casi la mitad la duración de trayectos y el diseño de aeronaves con mayor capacidad respecto a las antiguas aeronaves propulsadas por hélices. El siguiente paso lógico debía ser la aviación supersónica que podría potencialmente reducir a la mitad de nuevo las duraciones de los trayectos. Intentos por parte de las diferentes naciones potencias en la industria aeronáuticas no tardaron en verse culminando con la aparición del Concorde (Inglaterra- Francia), el soviético Túpolev Tu-144 (1978) y el prototipo Boeing 2707, que vieron su futuro truncado tanto por un lado por regulaciones en diferentes países que prohibían vuelos supersónicos sobres sus espacios aéreos; como por la crisis del petróleo de 1973, que multiplicó en cuestión de meses los precios de todos los derivaos del petróleo como el queroseno. Estas aeronaves se vieron avocadas a un uso de nicho en futuras décadas. Fue sin dudas esta crisis del petróleo la que estableció a la fuerza el nuevo paradigma aeronáutico que, aún cincuenta años después, sigue dominando el mundo de la aviación y su economía. Se comienzan a reducir costes y optimizar vuelos apareciendo los primeros trimotores de fuselaje ancho como el Lockheed L-1011 (1972) o el McDonnell Douglas DC-10, más baratos de operar que los cuadrimotores aun cumpliendo regulaciones como la ETOPS de la FAA, que no permitía a aeronaves de dos motores operar a una distancia mayor de una hora de vuelo de un aeropuerto. También supone un punto de inflexión en la retirada de muchos modelos operados durante las dos últimas décadas que tenían consumos muy por encima de lo que, ya para entonces, permitían el precio de los combustibles operar de forma rentable. Casos significativos pueden ser modelos como el Convair 990-Coronado (1961) o el Boeing 747-100 que aún recién aparecido en el mercado las aerolíneas ven al principio poco rentable el uso de una aeronave tan grande. Esta norma mermaba el uso emergente de aeronaves de fuselaje ancho bimotores en vuelos transoceánicos. No fue hasta la introducción de aeronaves como el Airbus A-300, la primera aeronave bimotor certificada que cumplía la norma ETOPS, precedente que sembró la semilla que hoy en día permite ver que la mayor parte de la flota de aeronaves comerciales a nivel global sean aeronaves bimotores tanto de corto, medio y largo alcance. Las próximas décadas se destacan por la aparición de un mayor número de rutas aéreas tanto de corta como larga distancia. En el mismo orden, la fundación de nuevas aerolíneas comerciales de bajo coste como RyanAir (Irlanda), que opera a nivel europeo, o WOW air (Islandia) con conexiones entre Europa y Estados Unidos o Southwest Airlines (EE. UU.). El enfoque empresarial de estas compañías, caracterizado por la eliminación de servicios adicionales gratuitos, la reducción de los tiempos de escala y una estructura de costes renovada han transformado radicalmente la forma en que operan las compañías aéreas tradicionales. Este cambio ha facilitado un acceso más económico al transporte aéreo para un mayor número de personas. La aviación, en sus 120 años de historia, ha representado una de las mayores revoluciones debido a su impacto multidimensional en la sociedad, la economía y la tecnología. Desde la conectividad global y la Introducción 4 rapidez en los desplazamientos hasta la innovación tecnológica en materiales y sistemas de propulsión, la aviación ha transformado radicalmente la forma en que vivimos, trabajamos y nos comunicamos. Además, ha desempeñado un papel crucial en eventos históricos, desde conflictos bélicos hasta avances científicos y descubrimientos geográficos. En resumen, la aviación ha trascendido como una fuerza impulsora de progreso y cambio en todo el mundo. 1.2. Historia de la aviación comercial en España En los orígenes de la aviación española destacan personajes como Pedro Vives y Emilio Herrera, que, en 1910, a través de la Comisión de Experiencias del Ministerio de la Guerra, adquirieron los primeros aeroplanos en Francia. El primer vuelo en España se realiza en España el 11 de febrero de 1910 con una aeronave modelo Bleriot XI (1908) pilotado por el francés Julien Mamet. Solo meses después, en 1911, Benito Loigorri Pimentel realizó el primer vuelo desde Ciudad Lineal a Carabanchel, convirtiéndose en el primer español con licencia de piloto, obtenida en Francia. A partir de ese mismo año, se impartieron las primeras clases para obtener la licencia en nuestro territorio. [5] La primera gran carrera aérea a nivel mundial se trata de la carrera Paris-Madrid en 1911. La gesta de cubrir la distancia de más de 1200 kilómetros entre las dos ciudades con aeronaves que apenas solían permitir vuelos de más de 20 minutos no deja de ser una prueba más del peligro que asumían los pioneros de la aviación. Solo una aeronave de las más de 30 que comenzaron la carrera consiguió cubrir la distancia total de la carrera pilotado por Jules Védrines y aterrizando en Getafe en la Dehesa de Santa Quintería donde finalmente se asentaría el primer aeródromo en España.[6] Los principales avances en aviación primigenia en España se dan fundamentalmente por parte de la aviación militar. En 1911 se crea un grupo experimental para la aviación dentro de la sección del cuerpo de ingenieros del ejecito nacional. España puede considerarse como la primera potencia de la aviación militar a nivel mundial, solo dos años después de la primera compra de aeronaves por parte del ministerio de guerra, se realizan operaciones militares aéreas en el protectorado de Marruecos. Los pilotos españoles llevaron a cabo el 17 de diciembre de 1913 la primera acción de bombardeo exclusivamente aérea en la historia global, al emplear bombas diseñadas específicamente para su uso en aviones contra las tribus rifeñas sentando precedentes como el uso programado de bombas aéreas o el primer uso registrado de una mirilla de bombardeo. [7] Hasta el inicio de la Primera Guerra Mundial la mayor parte de material para aeronaves venían de parte de las potencias que ahora tenían centrados todos sus esfuerzos de producción para la propia guerra. La gran demanda de material por parte de las potencias beligerantes impulsó los primeros desarrollos de material aeronáutico en España como lo fue el motor Hispano-Suiza 8. Este motor fue probablemente el más empleado en aeronaves del eje, un motor pionero V8 con refrigeración líquida que suponía un gran avance en potencia y fiabilidad sobre los motores rotativos de estrella con cigüeñal fijo que solían emplearse al inicio del conflicto. El diseño se licencia en diferentes países como Estados Unidos, Francia, Italia y Japón. España a pesar de estar al margen de la guerra, se coloca a la cabeza en tecnología aeronáutica. Tras el final de la Gran Guerra miles de aeronaves de diferentes países en guerra quedan abandonados por no tener ya un uso real. El stock mundial de aeronaves impulsa que la mayoría de los países vendan todos sus modelos a estrenar a precio de coste. Este hecho supone uno de los mayores desastres para el sector aeronáutico español pues cualquier nueva tecnología debía competir con precios absurdamente bajos. El propio Ministerio de Guerra español, realiza la compra de más de doscientos modelos franceses e ingleses a precios muy bajos. Por poner en contexto, se compraron aeronaves a precios entre 20.000 y 30.000 pesetas cuando el presupuesto anual del servicio de aviación militar eran unos 10 millones de pesetas. El propio 5 5 Análisis y optimización de la oferta de transporte aéreo en España ejército abandonó la industria aeronáutica nacional, que se encontraba en proceso de diseño de diferentes aeronaves como bombarderos, de reconocimiento o cazas. Situación que llevo a la industria nacional a meras labores de mantenimiento. [8] Las primeras rutas comerciales se dan lugar a lo largo de 1919, con el inicio de la ruta aeropostal entre Toulouse y Casablanca que realizaban escalas técnicas tanto en Barcelona, Alicante y Málaga. Estas operaciones dieron lugar al primer ordenamiento del transporte aéreo español regulándose las escalas en nuestro país. Esto fue seguido rápidamente de otras regulaciones nacionales ese mismo año, como el primer Reglamento de la Navegación Aérea Civil y el primer organismo de aviación civil como seria el Negociado de Aerostación y Aviación Civiles. Así pues, es en la década de los años 20 cuando aparecen las primeras compañías aéreas españolas como fueron la Compañía Española de Tráfico Aéreo (C.E.T.A) (1921), dedicada al servicio postal, la Compañía Aeromarítima Mallorquina (CAMSA) (1923). Sin embargo, no fue hasta el 28 de junio de 1927 cuando se funda la principal aerolínea española, Iberia con un trayecto Barcelona-Madrid.[9] Probablemente una de las mayores hazañas de la aviación española podría ser el vuelo del Plus Ultra. Este significó la primera vez que se volaba entre España (Palos de la Frontera) y Latino América (Buenos Aires) en una fecha (1926) en la que pocas ocasiones con anterioridad se había cruzado el Atlántico y un año antes de que se cursase por primera vez en un vuelo sin escala. Durante este vuelo se dio la primera ocasión en la aviación del uso de un radiogoniómetro para la navegación a lo largo de la etapa más larga dentro del océano Atlántico sentós precedente y estándar para la navegación de largas distancias durante más de cincuenta años en muchas rutas comerciales. [10] La primera conexión internacional no se da hasta el final de la guerra civil española ese mismo año con la ruta Madrid-Lisboa. Es en 1946 cuando Iberia abre la primera ruta que conecta Europa con el continente americano haciendo uso de un Douglas DC-4. Las evoluciones en la aviación a nivel mundial en las próximas décadas dejan las flotas de los operadores españoles obsoletas y no es hasta 1960 cuando la compañía Aviaco realiza la compra de las primeras aeronaves a reacción, con el modelo Sud Aviation Caravelle. Seguida de la adquisición por parte de Iberia en 1961 de varios modelos de la aeronave Douglas DC-8 que permiten ampliar el número de destinos internacionales en nuestro país. Es en esta misma época es cuando se modernizan de forma considerable los grandes aeropuertos nacionales, además de la construcción de otros tantos, habiendo a fecha de 1962 un total de 33 aeropuertos civiles. Hasta finales de la década de los 60, la aviación nacional estaba predominada por el uso de vuelos chárter. No es hasta la década de los 70 cuando la tendencia cambia a la apertura de rutas regulares entre los principales aeropuertos. Las siguientes décadas, se ven marcadas principalmente por la entrada de España en el transporte aéreo comunitario y la liberación del espacio aéreo, además de la fundación en 1990 de AENA (Aeropuertos Nacionales y Navegación Aérea), en la que recaen las competencias de gestión de red de aeropuertos y las instalaciones de ayuda a la navegación aérea y control de la circulación aérea. Con la liberalización del espacio aéreo español aparecen en estas últimas dos décadas nuevas compañías de bajo coste, como podrían ser los casos de Vueling o Volotea, que han aumentado considerablemente el número de movimientos y operaciones, llegando a un total en la actualidad de 49 aeropuertos de interés general a fecha de 2024. Desigualdad entre aeropuertos españoles y comparación con otros países 12 Figura 3: Evolución del número de pasajeros (Internacionales y nacionales) en la última década (datos numéricos en millones de pasajeros). 2.3. Distribución del número de pasajeros totales entre todos los aeropuertos españoles A partir de los datos del número total de pasajeros (tanto nacionales como internacionales) de cada uno de los aeropuertos españoles en el año 2023, se tendrá que considerar un total de 48 aeropuertos con al menos algún pasajero. En la figura 4; donde se ordenan los aeropuertos de menor a mayor número y se traza la curva Lorenz de pasajeros anuales, se puede ver fácilmente como existe una gran desigualdad entre los grandes aeropuertos y los pequeños. Solo 10 aeropuertos gestionan el entorno del 80% de los pasajeros totales. Analíticamente esta desigualdad se puede medir con el coeficiente de Gini, que tiene un valor de 0.7978~0.8 muy cercano a la unidad. La primera conclusión que se puede extraer es que el conjunto de pasajeros que pasan por todos los aeropuertos españoles se reparte de forma muy poco uniforme. Existen numerosos aeropuertos que apenas son empleados como podrán ser el de Córdoba, Albacete, Salamanca, Burgos… Esto viene a reforzar la teoría enunciada anteriormente de que la infraestructura aeroportuaria nacional no refleja la realidad en cuestiones de demanda. Como es lógico, la existencia de aeropuertos más grandes o con mayor número de pasajeros es un fenómeno común. Difícilmente se pueda comprar un mercado como podrían ser los que acogen pequeños aeropuertos como el de Córdoba, Badajoz o La Gomera con aquellos que se encuentran en la mayores ciudades o destinos turísticos nacionales. 2013 2014 2015 2016 2017 2018 2019 2020 2021 2022 2023 INTERNACIONALES 129.481 137.128 144.815 162.461 175.654 182.981 189.194 42.040. 67.563. 160.882 190.024 NACIONALES 57.609. 58.734. 62.313. 67.769. 73.261. 80.773. 85.728. 33.825. 52.174. 82.506. 92.853. 0 50.000.000 100.000.000 150.000.000 200.000.000 250.000.000 300.000.000 Pajaeros Numero de pasajeros (Internacionales y nacionales) en el preriódo 2013-2023 13 13 Análisis y optimización de la oferta de transporte aéreo en España Figura 4: Fracción acumulada de pasajeros sobre el total en orden creciente del conjunto de todos los aeropuertos españoles. Es por ello, que para obtener conclusiones más profundas, se procede a analizar cómo se reparte en cada uno de los principales aeropuertos por volumen de tráfico del país el conjunto de pasajeros con origen en este. Para ello, a partir del conjunto de los principales destinos de cada uno de los mayores 25 aeropuertos españoles por volumen de pasajeros se calcula el coeficiente de Gini de todas las rutas del aeropuerto estudiado que tenga más de 100 pasajeros en 2023. La muestra para el cálculo del coeficiente de Gini de cada aeropuerto consiste en las 50 principales rutas por volumen de tráfico de dicho aeropuerto si tuvieron más de 100 pasajeros en el año 2023. En el caso de que el aeropuerto estudiado no disponga de tantas rutas (principalmente los de menor tamaño) se consideraran todas aquellas rutas en las que al menos si hubo 100 pasajeros. En la figura 5 se presenta la transformación logarítmica en el eje x del conjunto total de pasajeros de cada aeropuerto y el eje y, el cociente entre el coeficiente de Gini entre uno menos dicho coeficiente (Coeficiente de Gini: C.G, dentro de la figura). La figura muestra además sendas rectas de mínima dispersión de los coeficientes calculados. La línea verde realiza un ajuste mediante mínimos cuadrados ordinarios (MCO) de todos los puntos de la muestra, mientras que la roja realiza el mismo ajuste mediante el método de regresión Theil-Sen. El método de regresión de Theil-Sen, es un enfoque robusto para ajustar una línea de regresión a un conjunto de datos. A diferencia de los métodos de regresión lineal tradicionales, como el método de mínimos cuadrados ordinarios (MCO), que pueden ser sensibles a valores atípicos o datos ruidosos, el método de Theil-Sen es más resistente a estas influencias y proporciona estimaciones más estables. El método de Theil-Sen es robusto frente a valores atípicos y datos ruidosos en el conjunto de datos. Esto significa que su rendimiento no se ve tan afectado por puntos extremos como si lo haría el método de mínimos cuadrados ordinarios. En lugar de calcular una sola línea de regresión como en MCO, el método de Theil-Sen calcula todas las HSK LEC ABC RGS ODB QSA LES RJL SLM JCU LEN BJZ JCU GMZ PNA VLL VDE VIT EAS MLN ZAZ LEI RMU YRX GRX REU VGO SDR LCG GRO OVD SCQ MAH FUE TFN BIO SVQ ACE IBZ VLC TFS LPA ALC AGP PMI BCN MAD Aeropuertos 0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 Sumatorio de Pasajeros Normalizado Coeficiente de Gini: 0.79783 Curva Lorenz Línea de Igualdad Desigualdad entre aeropuertos españoles y comparación con otros países 14 posibles pendientes entre pares de puntos y luego toma la mediana de estas pendientes. Esta mediana de pendientes es robusta frente a valores extremos y proporciona una estimación más resistente al ruido. A diferencia de algunos métodos paramétricos que asumen normalidad en los datos, el método de Theil-Sen es no paramétrico y no requiere supuestos específicos sobre la distribución de los datos. [16] Este algoritmo proporciona tanto una intersección estimada (b), como una pendiente estimada (m), de forma similar a la que se obtiene realizando el ajuste por mínimos cuadrados ordinarios. Observando dichas rectas y los principales valores de los estadísticos de la regresión por mínimos cuadrados de la tabla 2, se puede extraer que el ajuste es significativo al 1% con un R2 de 0.636. Por lo que se puede concluir que, medida que el volumen de tráfico del aeropuerto disminuye peor uniformemente repartido está estarán sus pasajeros entre todas las conexiones que ofrece dicho aeropuerto. Dicho de otra manera, a medida que un aeropuerto tiene menor número de pasajeros, existirán menos rutas en dicho aeropuerto además de que existirán rutas con mucho tráfico y otras muchas más rutas con muy poco tráfico. Figura 5 : Coeficiente de Gini para los mayores 25 aeropuertos de España en volumen de tráfico nacional e internacional en función del número total de pasajeros de dicho aeropuerto Si se excluye de la muestra todas las conexiones internacionales y se pone el foco de atención únicamente en el mercado nacional de pasajeros se podrán obtener conclusiones más esclarecedoras sobre el mercado nacional de la aviación comercial. En la figura 6 se realiza el mismo estudio solo contemplando los pasajeros de origen nacional empleando todas las rutas nacionales que tuvieron al menos 100 pasajeros en dicho año. Al analizar los resultados obtenidos del ajuste por mínimos cuadrados de una muestra de 25 coeficientes de Gini del reparto de pasajeros entre los principales destinos de los 25 mayores aeropuertos a nivel nacional, se pueden extraer conclusiones similares al anterior caso estudiado, el ajuste es significativo al 1% con un R2 de 12.5 13 13.5 14 14.5 15 15.5 16 16.5 17 Logaritmo(Pasajeros nacional + internacional por aeropuerto) -1 -0.5 0 0.5 1 1.5 2 2.5 Logaritmo(G/(1-G)) Parámetros TheilSen Least Squares 15 15 Análisis y optimización de la oferta de transporte aéreo en España 0.187. Se puede concluir que para el tráfico nacional también ocurre que a medida que el tráfico del aeropuerto disminuye se reparte de forma menos uniforme dicho tráfico entre las diferentes rutas del aeropuerto. Figura 6 Coeficiente de Gini para los mayores 25 aeropuertos de España en volumen de tráfico exclusivamente nacional en función del número total de pasajeros de dicho aeropuerto Tabla 2 Parámetros del ajuste de las figuras 5 y 6 empleando el método de mínimos cuadrados ordinarios (MCO) y Theil-Sen MCO Theil-Sen Ajuste p R2 b m b m Nacional e Internacional 1.8e-06 0.636 7.9594 -0.52721 8.416954 -0.554425 Nacional 0.031 0.187 3.2075 -0.15897 3.098484 -0.158871 Dónde los parámetros de la tabla 2 m y b son respectivamente los valores de la pendiente de la recta de ajuste y su valor inicial. 9 10 11 12 13 14 15 16 Logaritmo(Pasajeros nacionales por aeropuerto) 0 0.5 1 1.5 2 2.5 Logaritmo(G/(1-G)) Parámetros TheilSen Least Squares Desigualdad entre aeropuertos españoles y comparación con otros países 16 Una vez analizado la relación entre el coeficiente de Gini y el número de pasajeros de los principales aeropuertos españoles. Es interesante añadir una comparativa de dichos coeficientes de Gini en función de si solo se contempla tráfico nacional o y se emplea el conjunto de pasajeros de dicho aeropuerto (figura 7), dónde se ordenan los aeropuertos de mayor a menor volumen de tráfico total. Por lo general el coeficiente de Gini de pasajeros nacionales es más alto, exceptuando casos como Bilbao, A Coruña o Vigo donde son muy similares o Santiago-Rosalía de Casto y Oviedo donde el nacional está muy por debajo del conjunto total. Destaca el aeropuerto de Santiago-Rosalía de Castro, que posee un coeficiente de Gini bastante bajo en comparación con el resto. Este posee conexiones importantes con la mayoría de los aeropuertos nacionales. En contraposición, los altos coeficientes obtenidos para aeropuertos como Girona o Vigo, que se deben principalmente a que ambos solo poseen conexiones regulares con Madrid-Barajas, ambos suelen recibir la mayor parte de vuelos desde destinos internacionales y que ambos se encuentran muy cerca de aeropuertos de mayor importancia, como lo son respectivamente Barcelona-El Prat o Santiago-Rosalía de Castro. También se puede destacar que para aeropuertos en zonas de alto interés turístico internacional sus coeficientes suelen ser muy bajos (Gran Canaria, Tenerife Sur o Alicante) en comparación con los obtenidos empleando el mercado nacional. En contraposición en aquellos en los que ambos coeficientes son muy similares suelen encontrarse en zonas donde el turismo internacional no tiene tanta relevancia. La mayoría de grandes aeropuertos españoles en tráfico total suelen tener coeficientes bajos (por debajo de 0.5) lo que puede mostrar que en España existe un número considerable de aeropuertos con una buena red de rutas a nivel internacional. Figura 7 Coeficientes de Gini calculados para cada aeropuerto empleando tráfico nacional e internacional (azul) y exclusivamente nacional (rojo) MAD BCN PMI AGP ALC LPA TFS VLC IBZ ACE SVQ BIO TNF FUE MAH SCQ OVD GRO LCG SDR VGO GRX YRX RMU ZAZ Aeropuertos 0.3 0.4 0.5 0.6 0.7 0.8 0.9 Coeficientes de Gini Coeficientes Gini tráfico internacional y nacional Coeficientes Gini tráfico nacional 17 17 Análisis y optimización de la oferta de transporte aéreo en España 2.4. Modelo de máxima entropía Se procede a emplear el método de máxima entropía explicado en el apartado 2.1.2 para ver, a partir de una matriz que engloba el conjunto de pasajeros que se desplazan de un aeropuerto i a otro aeropuerto j. Para ello se utilizará datos de pasajeros nacionales de los mismos aeropuertos que se emplearán posteriormente en el modelo del apartado 3. Una vez se tiene la matriz de pasajeros T y conociéndose los vectores Oi (número de viajes generados por la zona i) y Dj (que indica el número de viajes atraídos por cada zona de destino j). Se procede a calcular los coeficientes de Gini del número de pasajeros que se desplaza de cada aeropuerto al resto de aeropuertos de la red, como ya se hizo en la figura 7 y se comprara el resultado con los obtenidos mediante el método de máxima entropía. En la figura 8 se observa cómo los coeficientes de Gini para el conjunto de pasajeros generados por el método de máxima entropía suelen ser bastante menor para todos los aeropuertos estudiados y exceptuando el caso del aeropuerto Madrid-Barajas es más bajo que el coeficiente de Gini calculado para el conjunto de pasajeros reales de cada aeropuerto. Por lo general, a más alto sea el coeficiente de Gini real obtenido para un aeropuerto mayor será la diferencia con la distribución de máxima entropía. La razón principal por la que los coeficientes de Gini sean más similares entre sí en el escenario de máxima entropía es la uniformidad de la distribución de pasajeros. La máxima entropía busca eliminar cualquier preferencia o sesgo en la distribución, resultando en una dispersión equitativa de pasajeros; en contraste con los datos reales, donde factores externos crean distribuciones más dispares. De media, esta mejora del coeficiente de Gini si se comparan los resultados obtenidos mediante el método de máxima entropía en contraposición de los datos reales se encuentra en 0.15 puntos por debajo. Figura 8: Comparación de los coeficientes de Gini reales y los calculados por el método de máxima entropía. MAD PMI BCN LPA TFN IBZ BIO ACE VGO LCG AGP SCQ SVQ FUE TFS MAH OVD VLC ALC EAS SDR GRX ZAZ XRY Aeropuertos 0.45 0.5 0.55 0.6 0.65 0.7 0.75 0.8 0.85 0.9 0.95 Coeficiente de Gini Comparacion indices de gini operaciones reales vs operaciones de matriz de maxima entropia Coeficiente gini real Coeficiente gini maxima entropia Desigualdad entre aeropuertos españoles y comparación con otros países 18 Esta diferencia se puede observar de forma más detenida si se escogen algunos de los aeropuertos estudiados. En la figura 9 se muestra una comparación del sumatorio de pasajeros normalizado para los aeropuertos de Barcelona-El Prat, Mallorca, Ibiza y Sevilla. Como ya se pudo ver en la figura 8 el caso del aeropuerto de Palma de Mallorca es el que el caso de máxima entropía se asemeja más a la realidad. Para los otros aeropuertos mostrados, y en especial para el aeropuerto de Ibiza, la curva Lorenz en el caso de máxima entropía muestra un comportamiento mucho más uniforme que el que la curva Lorenz real. Figura 9 : Comparación del sumatorio de pasajeros normalizado tanto con datos reales como los calculados por el método de máxima entropía para los aeropuertos de Barcelona-El Prat, Palma de Mallorca, Ibiza y Sevilla. Si bien es interesante entender cómo se relacionan las distribuciones de pasajeros reales y de máxima, se puede profundizar más en el estudio si comparamos de forma cuantitativa la relación entre pasajeros reales y de máxima entropía. En la figura 10 se comparan dichos valores de forma logarítmica, donde para números de pasajeros por ruta grandes ( log*(𝑝𝑎𝑠𝑎𝑗𝑒𝑟𝑜𝑠)>9 ) existe una relación lineal entre los pasajeros reales por ruta y los calculados por el método de máxima entropía (figura 10 derecha). El coeficiente de determinación R2 indica que alrededor del 66.4% de la variabilidad del coeficiente de Gini de cada aeropuerto depende del volumen de tráfico de dicho aeropuerto. La estadística F es alta y su valor p es muy por debajo del umbral de confianza estándar de 0.05 (p=2.47e-72), lo que indica que el modelo en su conjunto es significativo. Es decir, para rutas con un número considerable de pasajeros la diferencia entre dicho número de pasajeros y los de máxima entropía es pequeña. Por otro lado, dicha relación entre pasajeros podría estar influenciada por la distancia entre los aeropuertos. Mediante la ecuación 2-5, es posible estudiar dicha relación. log ( 𝑃𝑎𝑠𝑎𝑗𝑒𝑟𝑜𝑠 ) =log ( 𝑃𝑎𝑠𝑎𝑗𝑒𝑟𝑜𝑠*𝑚á𝑥𝑖𝑚𝑎*𝑒𝑛𝑡𝑟𝑜𝑝í𝑎 ) + a ×𝐷𝑖𝑠𝑡𝑎𝑛𝑐𝑖𝑎*𝑒𝑛𝑡𝑟𝑒*𝑎𝑒𝑟𝑜𝑝𝑢𝑒𝑟𝑡𝑜𝑠** ( 2−5 ) 0 5 10 15 20 Aeropuertos 0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 Sumatorio de Pasajeros Normalizado Sumatorio de Pasajeros Normalizado aeropuerto de Palma de Mallorca Datos Palma de Mallorca Linea guía Maxima entropía Mallorca 0 5 10 15 20 Aeropuertos 0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 Sumatorio de Pasajeros Normalizado Sumatorio de Pasajeros Normalizado aeropuerto Barcelona-El Prat Datos reales Barcelona-El Prat Linea de Uniformidad Maxima entropia Barcelona-El Prat 0 5 10 15 20 Aeropuertos 0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 Sumatorio de Pasajeros Normalizado Sumatorio de Pasajeros Normalizado aeropuerto Ibiza Datos reales Ibiza Linea de Uniformidad Maxima entropia Ibiza 0 5 10 15 20 Aeropuertos 0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 Sumatorio de Pasajeros Normalizado Sumatorio de Pasajeros Normalizado aeropuerto Sevilla Datos reales Sevilla Linea de Uniformidad Maxima entropia Sevilla 19 19 Análisis y optimización de la oferta de transporte aéreo en España Los valores obtenidos de Alpha disminuyen de forma progresiva a medida que aumenta la distancia entre aeropuertos, pero en general son muy cercanos a 0. Un valor de Alpha cercano a cero sugiere que la distancia entre aeropuertos tiene poco impacto en la relación entre los pasajeros reales y los de máxima entropía. En otras palabras, la distancia no está contribuyendo significativamente a explicar la variabilidad en los pasajeros reales más allá de lo que ya explica la máxima entropía. El modelo de regresión lineal parece ser significativo en términos estadísticos, pero explica una proporción relativamente baja de la variabilidad en la variable Alpha y utilizando la distancia como variable predictora. Los coeficientes son significativos, el modelo tiene una capacidad predictiva limitada según el R2 Figura 10 Relación entre pasajeros reales y de máxima entropía (arriba) e influencia de la distancia en la diferencia de dichos valores (abajo) Tabla 3 Valores estadísticos más importantes de los ajustes realizados en la figura 10 MCO Theil-Sen Ajuste p R2 b m b m Pasajeros reales y máxima entropía 2.47e-72 0.664 2.984 0.72198 3.094796 0.719035 Alpha 2.44e-14 0.178 0.4123 -0.0593 0.475388 -0.067676 0 2 4 6 8 10 12 14 Logaritmo pasajeros reales 7 8 9 10 11 12 13 14 Logaritmo pasajeros de maxima entropía Compraración del numero de pasajeros reales con pasajeros de maxima entropía 9 9.5 10 10.5 11 11.5 12 12.5 13 13.5 14 Logaritmo pasajeros reales 8.5 9 9.5 10 10.5 11 11.5 12 12.5 13 13.5 14 Logaritmo pasajeros de maxima entropía Zona ampliada 0 500 1000 1500 2000 2500 Distancia entre aeropuertos -0.4 -0.3 -0.2 -0.1 0 0.1 0.2 0.3 0.4 Alpha Valores de los coeficientes alfa en función de la distancia de la ecuación: Log(Pax)=Log(Pax-Max-entropia)+alpha*Log(distancia) Desigualdad entre aeropuertos españoles y comparación con otros países 20 2.5. Comparación del ordenamiento de los pasajeros en los aeropuertos españoles con otros países Una vez conocida la distribución de pasajeros en España, es conveniente comparar los resultados y conclusiones obtenidos con anterioridad con otros países. Dentro del marco europeo, España se sitúa como el país que más pasajeros pasan por sus aeropuertos, seguido de Alemania, Francia e Italia. En la tabla 4 puede verse el porcentaje de pasajeros de los principales países de la Unión Europea del año 2023 que reportan datos a la agencia Eurostat, siendo estos un total de 32 países.[17] Tabla 4 Porcentaje de pasajeros de los principales países de la Unión Europea País % Del total de pasajeros de la Unión Europea España 17.276 Alemania 13.403 Francia 11.681 Italia 11.991 Países Bajos 5.1478 Portugal 4.9368 Suiza 3.7621 Polonia 3.6873 Irlanda 2.8719 Al igual que se hizo en la figura 3 podemos calcular cómo se distribuye el conjunto total de los pasajeros entre aeropuertos de diferentes países. Tendremos en consideración los principales 27 aeropuertos por volumen de pasajeros de 5 países diferentes, siendo estos España, Reino Unido, Francia, Alemania y Estados Unidos. Los datos empleados han sido recopilados de los diferentes organismos gubernamentales encargados de la distribución de estos. Se tiene pues, para España AENA [11], para Inglaterra “UK Civil Aviatión Authority” [18], De Francia ‘’Union des Aeroports Français” [19]. De Alemania “ADV Flighafenrband” [20] y finalmente para Estados Unidos “Bureau of transportation stadistics” [21]. Observando la figura 7 y la tabla 3, podemos distinguir tres comportamientos claramente diferenciables. El caso de Estados Unidos, que presenta el comportamiento mucho más uniforme que el resto de los países estudiados, España y Reino Unido en un término intermedio y los casos más desiguales, Francia y Alemania. 21 21 Análisis y optimización de la oferta de transporte aéreo en España Tabla 5 Coeficiente de Gini del reparto de pasajeros entre los principales 27 aeropuertos de cada país País Coeficiente de Gini (Top 27 aeropuertos por país) Estados Unidos 0.27357 España 0.66165 Reino Unido 0.70017 Francia 0.72744 Alemania 0.79892 Por un lado, Estados Unidos distribuye sus pasajeros entre aeropuertos de forma muy uniforme. El mercado de la aviación comercial americana es el mayor del mundo. Estados Unidos es un país que, por su propia geografía, la población se ve generalmente forzada al uso de la aviación para poder moverse tanto fuera y dentro del país. Es por ello por lo que gran parte de los grandes aeropuertos de Estados Unidos copan las listas de mayores pasajeros por volumen de pasajeros En el ranking de los 50 mayores aeropuertos mundiales por pasajeros, Estados Unidos posee 17. Mientras que España solo posee 2 (Madrid-Barajas y Barcelona-El Prat), Inglaterra 2 (Heathrow y Gatwick), Francia 1 (Paris-Charles de Gaulle) y Alemania 1 (Frankfurt). [21] Esto ya nos adelanta que, probablemente el conjunto de pasajeros totales, que, aunque sea el mayor del mundo, para los principales aeropuertos americanos se reparta de forma muy uniforme como se puede ver en la figura 11. Para Inglaterra y España parece que el reparto es relativamente similar. Aunque en Inglaterra los aeropuertos más pequeños parecen tener una menor importancia nacional que los españoles. Primeramente, tenemos el caso de Heathrow que absorbe en el entorno del 30% de tráfico de pasajeros del país, seguido de una seria de aeropuertos menores que copan cada uno en el entorno del 10-15% del tráfico total. Estos aeropuertos serian Gatwick, Manchester, Stanted y Luton. En el caso español, esto se puede ver extrapolado en el aeropuerto de Madrid-Barajas, seguido, pero mucho más de cerca que para los casos ingleses de Barcelona, Palma de Mallorca y en menor medida Málaga-Costa del Sol o Alicante. Este fenómeno también se ve reflejado en la figura 11, si observamos las líneas Azul (Reino Unido) y roja (España) vemos como el peso de Heathrow es mayor que el de Madrid-Barajas y para los aeropuertos de Barcelona-El Prat, Palma de Mallorca y Málaga- Costa del Sol mayores que los ingleses Gatwick, Manchester, Stanted y Luton. En el caso francés es muy similar al alemán, la uniformidad del reparto de pasajeros entre sus aeropuertos es muy baja. Esta desigualdad se amplifica aún más si entendemos la localización de sus aeropuertos. Francia posee principalmente dos grandes aeropuertos y ambos se encuentran en su capital (Paris-Charles de Gaulle y Paris–Orly). Si se suman los pasajeros de ambos aeropuertos obtendríamos que, en realidad solo París englobaría alrededor del 50% del conjunto total de pasajeros de todos los aeropuertos, cifra que en ninguna circunstancia se puede catalogar como equilibrada. Entre todos los aeropuertos de Alemania, el aeropuerto de Frankfurt se destaca por su alto tráfico de pasajeros, manejando alrededor del 34% del total de pasajeros que pasan por los aeropuertos alemanes. Esta cifra, representa una parte significativa del flujo aéreo nacional y refleja la importancia y el volumen de operaciones Modelo de estimación de vuelos necesarios para cubrir la demanda en una red de aeropuertos 28 demanda existente. Matemáticamente está descripción quedaría limitada por la siguiente expresión: 𝑓!" ≤* # 𝑏,*𝑦!" , ) ,&' * ( 3−4 ) • Máximo número de operaciones de cada aeropuerto: Es necesario limitar de alguna manera el número de operaciones de cada aeropuerto. La restricción matemática que modelaría esta restricción sería la siguiente: ##( 𝑦!" ,+𝑦"! , ) ≤*𝑣! ) ,&' 2 "&' "3! * ( 3−5 ) Donde 𝑣! sería el numero máximo de operaciones realizables entre cada par de aeropuertos. Aunque esta variable puede ser difícil de modelar, más aún, teniendo en cuenta que se van a realizar modelos en los que no se incluyen la totalidad de operaciones de cada aeropuerto. Pues se va a estudiar en función de la demanda de pasajeros nacionales y obviando la existencia de operaciones de tipo internacionales. Para la simplificación de los cálculos, en el conjunto de modelos realizados esta restricción se puede acotar de tres maneras diferentes. 1) Restricción débil Por un lado, el número de operaciones máximos que se realizan en un aeropuerto puede ser la suma de las componentes de cada columna de una matriz cuadrada de operaciones, esto generaría un vector 𝑣 que tendría longitud n. Por lo que tendríamos un total de n restricciones, las mismas que aeropuertos tenga el modelo. Esto queda expresado matemáticamente por la siguiente expresión: ##( 𝑦!" , ) ≤* # 𝑣!" 2 !&' ) ,&' 2 "&' * ( 3−6 ) 2) Restricción fuerte La restricción anterior se puede volver más estricta aún si en lugar de realizar la suma de cada columna, hacemos que para cada aeropuerto y cada ruta el número máximo de operaciones sea menor o igual a un nuevo de operaciones dado, es decir en lugar de realizar las sumas de las columnas comprar cada elemento de la matriz de restricciones con el valor de operaciones dado. Esto queda matemáticamente definido cómo: 𝑦!" '+𝑦!" 4+⋯+*𝑦!" )=𝑉!"* ( 3−7 ) Aplicada a la realidad supone que para una matriz demanda F y una matriz de operaciones V dadas a partir de 29 29 Análisis y optimización de la oferta de transporte aéreo en España datos reales. Nunca se podrá realizar para cada par origen-destino más operaciones de las que realmente se realizaron. En este caso se tendrían un total de n x n restricciones. 3) Restricción mixta Puede ocurrir que, para las dos restricciones de operación anteriormente desarrolladas, la primera sea demasiado abierta. Esta puede dar lugar a casos en los que se prefiera únicamente un tipo de aeronave pequeña para operar una ruta con mucha demanda lo que acarrearía un número de operaciones demasiado elevado. Por otro lado, la segunda restricción puede llegar a ser demasiado estricta. Para casos en los que el número de operaciones entre aeropuerto sea pequeño y deba emplear varias aeronaves diferentes para realizar un número total bajo de operaciones. Para equilibrar este desbalance se plantea la siguiente restricción. Por un lado, se usará la primera restricción de forma general, es decir, el número de operaciones totales en un aeropuerto no puede superar un valor establecido. A esta restricción se le añaden otras nuevas para aquellos casos en los que la demanda de pasajeros entre dos aeropuertos sea muy alta (Mayor que la demanda que usaremos como umbral), para estas rutas se especifica un número máximo de operaciones. Esta restricción se puede denominar como mixta, pues busca el balance de las otras dos. Esta nueva restricción queda definida matemáticamente de la siguiente forma: ##( 𝑦!" , ) ≤* # 𝑣!" 2 !&' ) ,&' 2 "&' 𝑦!" '+𝑦!" 4+⋯+*𝑦!" )=𝑣!"****𝑝𝑎𝑟𝑎*𝑓!" >𝐷𝑒𝑚𝑎𝑛𝑑𝑎*𝑢𝑚𝑏𝑟𝑎𝑙* ( 3−8 ) Esta será la restricción que se empleará en el modelo definitivo, pues dentro de su flexibilidad, en ningún caso se encontrarán resultados muy alejados de los que realmente ocurren. • Restricción del uso de la aeronave de mayor tamaño para la realización de vuelos de muy corta duración. En la ordenación del tráfico aéreo español nos encontramos con un fenómeno a destacar. Las diferentes islas de los archipiélagos nacionales suelen estar conectadas por rutas aéreas de muy alta frecuencia que permiten a los habitantes de estas islas una mayor facilidad para su movilidad en el archipiélago. El caso típico, si se piensa en este tipo de vuelos, serían los operados por la aerolínea Binter Canarias en cuya flotan dispone de un total de 27 aeronaves del tipo ATR-72, aeronaves con una capacidad entre los 68-78 pasajeros, algo menores a la de las aeronaves que emplearemos. Carecería de sentido emplear una aeronave como podría ser un Boeing- 737, Airbus A-320 o mayor. Es por ello por lo que, en el modelo se va a incluir una restricción adicional. Esta restricción elimina la posibilidad de uso de las aeronaves de mayor tamaño en este tipo de vuelos tan cortos. Por lo tanto, para una distancia entre aeropuertos menor a una distancia umbral se elimina la posibilidad del uso de una aeronave que no sea la más pequeña que contemple el modelo. Esta restricción queda expresada de forma matemática como: 𝑦!" 4+𝑦!" 5+⋯+*𝑦!" )=0********∀*𝑟𝑢𝑡𝑎*𝑐𝑜𝑛*𝑑!" <𝐷𝑖𝑠𝑡𝑎𝑛𝑐𝑖𝑎*𝑢𝑚𝑏𝑟𝑎𝑙* ( 3−9 ) Modelo de estimación de vuelos necesarios para cubrir la demanda en una red de aeropuertos 30 Es por ello por lo que si las aeronaves están ordenadas de menor a mayor tamaño solo se podrá emplear la primera (Que sería la más pequeña). 3.2 Implementación en MATLAB® La explicación más detallada de la creación del modelo y su funcionamiento se desarrolla en los ANEXOS A y B. En términos generales, el modelo se basa en la Programación Lineal de Enteros Mixtos (MILP), solo haciendo uso de valores enteros. Matlab, dispone de numerosas funciones que permiten la resolución de problemas de esta naturaleza como la función intlingprog. En el ANEXO B de detallan la implementación de diferentes restricciones desarrolladas en el apartado 3.1.2, así como el proceso de lectura de datos necesarios para el correcto funcionamiento y el tratamiento de los resultados finales.[27] 31 31 Análisis y optimización de la oferta de transporte aéreo en España 3.3 Datos empleados para la creación del modelo El conjunto de simulaciones que se realizarán más adelantes, están basadas en datos estadísticos reales de un conjunto de aeropuertos españoles escogidos siguiendo diferentes criterios. Estos datos estadísticos, serán tanto el número total de pasajeros que vuelan de cada aeropuerto a el resto de los aeropuertos de la red durante un periodo de tiempo determinado, así como, el conjunto de operaciones comerciales realizadas en ese mismo perdido. En este apartado se mostrarán con detenimiento todos los datos y variables necesarias para la correcta ejecución del modelo. 3.3.1 Aeropuertos empleados Se van a emplear, un conjunto total de 24 aeropuertos españoles, los cuales se detallan en la figura 13. Los aeropuertos han sido seleccionados buscándose un compromiso entre aeropuertos con un tráfico al menos moderadamente destacable y que ninguna región significativa se quede a grandes distancias de algún aeropuerto. Este compromiso se complica al intentar escoger algún aeropuerto del interior de la península. Como ya se explicó con anterioridad gran parte de los aeropuertos del interior de España suelen tener un número de pasajeros muy reducido. Esto ocurre con aeropuertos como Valladolid, Salamanca, Burgos, Teruel, Badajoz o Albacete. Su inclusión en el modelo, si queremos usar datos reales, no arrojaría resultados de gran interés. En cuanto a los aeropuertos insulares se ha buscado la mayor inclusión posible. En general estos tienen tres posibles usos. Por un lado, a lo largo de todo el año, ofrecen una conexión entre las islas de sendos archipiélagos además de permitir su conexión con la península. Por otro lado, en periodos estivales suelen ser aeropuertos que acumulan una gran cantidad de pasajeros al ser estas islas algunos de los puntos turísticos nacionales de mayor importancia. De todos los aeropuertos insulares, solo se han excluido los respectivos de las islas de El hierro, La Gomera y La Palma. La figura 13 puede compararse con la figura 2 para observar mejor cuáles han sido los aeropuertos españoles descartados para la realización del modelo. Figura 13 Aeropuertos españoles empleados en el modelo. Modelo de estimación de vuelos necesarios para cubrir la demanda en una red de aeropuertos 32 3.3.2 Aeronaves empleadas En cuanto a las aeronaves, empleadas, por limitaciones de computación solo será posible la elección entre dos tipos de aeronaves. Para realizar la selección de las dos aeronaves que serán empleadas en el modelo, en primer lugar, se procede a analizar la flota de las mayores aerolíneas que operan en España. Las principales cinco aerolíneas, en volumen de pasajeros en el año 2023 en España son las mostradas en la tabla 5: Tabla 7 Número de pasajeros transportados por las cinco operadoras comerciales más importantes en España por volumen de tráfico y casos de Binter Canarias y Air Nostrum. [28] Aerolínea Pasajeros transportados en 2023 (Millones) RyanAir 52,17 Vueling 40,35 Iberia 17,89 Air Europa 14,37 EasyJet 13,59 Binter (Especial) 8,65 Air Nostum (Especial) 7,49 Es interesante añadir a la lista anterior el caso de las aerolíneas Binter Canarias y Air Nostrum, las cuales poseen una flota a considerar de dos modelos de aeronaves muy interesantes para este trabajo. Estas aeronaves son el Embraer-190 y Bombardier CJR-1000, aeronaves cuya capacidad oscila un total de 100 pasajeros.[29] En la tabla 8 se muestra la flota combinada de aeronaves de estas 7 aerolíneas, con sus respectivas capacidades de pasajeros en configuración de una y dos clases: Tabla 8 Flota conjunta de todas las aerolíneas contempladas en la tabla 5 [30] Modelo Flota total Pasajeros 2 clases Pasajeros 1 clase Airbus A319-100 31 124 134 Airbus A320 Neo 76 150 164 Airbus A320-200 155 150 164 Airbus A321-200 104 185 199 Airbus A330-200 14 293 - Airbus A330-300 8 335 - Airbus A321 Neo 13 185 220 Boeing 737-800 229 160 189 ATR 72 26 - 68 Embraer 190 19 - 98 Bombardier CJR-1000 29 - 104 33 33 Análisis y optimización de la oferta de transporte aéreo en España Observando el total de aeronaves disponibles de cada modelo de las aerolíneas estudiadas, las aeronaves Airbus A320-200, Airbus A321-200 y Boeing 737-800, suponen la mayor parte de la flota de estas aerolíneas. No obstante, sus capacidades son muy similares entre si rondando cifras entre los 150-200 pasajeros, por lo que, a priori podríamos elegir cualquiera de ella como la aeronave de mayor tamaño. Para el caso de la aeronave pequeña, su capacidad debe oscilar unos 100 pasajeros, como lo es para las dos aeronaves de estas características (Embraer-190 y Bombardier CJR-1000). El uso del modelo ATR-72 se ve limitado desde un principio por su alcance, de unos 1500 km lo cual no permitiría su uso para las rutas más largas a cubrir. Además, no es una aeronave destinada al tipo de operación que se le exigirá en este modelo, pues suele emplearse en trayectos cortos, rápidos y de bajo coste. Dicho todo lo anterior las aeronaves que se van a emplear en este modelo serán los modelos Boeing 737-800 y Embraer 190, ambas en una configuración de pasajeros de una clase. En la siguiente tabla se detallan las características más importantes de ambas aeronaves necesarias para este estudio: Tabla 9 Principales características de las aeronaves seleccionadas [31] [32] El alcance máximo de ambas aeronaves nos permite tener la certeza de que ambas podrán realizar cualquiera ruta dentro del modelo pues la distancia máxima contemplada es de 2305 km que se da en la ruta entre los aeropuertos de Tenerife Sur y Menorca En el apartado de resultados se realizarán diferentes simulaciones donde se emplearán otras aeronaves. Los datos necesarios de dichas aeronaves se enumerarán y justificarán en el apartado donde se empleen. Modelo Capacidad de pasajeros MTOW (Toneladas) Alcance máximo (Kilómetros) Embraer 100 47.8 3334 Boeing 737-800 180 68 5665 Figura 14 Boeing 737-800 (RyanAir) Figura 15 Embraer 190 (Binter) Modelo de estimación de vuelos necesarios para cubrir la demanda en una red de aeropuertos 34 3.3.3 Datos de pasajeros (demanda), operaciones y distancias entre aeropuertos En cuanto a datos de demanda, para cada escenario simulado se han empleado datos reales de número de pasajeros y operaciones realizadas correspondientes a cada periodo estudiado, en el caso de este trabajo se emplearán datos de pasajeros y operaciones del mes de agosto de 2022. Estos estarán dispuestos en sendas matrices cuadradas de dimensión n (n =24), cuyas diagonales son nulas pues no tiene sentido realizar operaciones para transportar pasajeros de un aeropuerto a el mismo y donde cada columna corresponde con el número de pasajeros y operaciones con origen en el aeropuerto i al aeropuerto j. Esta matriz no conoce el número de pasajeros que quiere volar a un destino con escala previa en otros. Por ejemplo, supongamos que un pasajero haya volado desde Santander a Gran Canaria con escala en Madrid, en la matriz de demanda este pasajero aparecerá en dos ocasiones, tanto en un vuelo de Santander a Madrid, como en otro de Madrid a Gran Canaria. El número de operaciones reales servirá a su vez como la restricción del número de operaciones posibles a realizar para cada aeropuerto y servirán en cualquier caso como comparación de los resultados generados por el modelo. Es interesante realizar un pequeño análisis previo de los datos que se van a emplear. En la tabla 10, se muestran el conjunto total de pasajeros con origen en cada aeropuerto de la red, así como el total de operaciones que se realizaron en dicho aeropuerto, además del número de pasajeros medio de cada operación realizada en el aeropuerto. Como se puede observar, el número de pasajeros medio por operación de cada aeropuerto se encuentra en todos los casos entre las capacidades de las dos aeronaves empleadas. En términos globales en el mes de agosto de 2022 dentro de la red de aeropuertos empleadas hubo un total de 4.043.189 pasajeros y se realizaron un total de 32.200. Esto resulta en una media de pasajeros por operación de 125.56. A nivel local, si se observan estas cifras para cada aeropuerto. Destaca el hecho de que en los aeropuertos de las Islas Canarias el número de pasajeros por operación es bajo. Esto se debe principalmente al hecho de que muchas de las operaciones realizadas en dichos aeropuertos son de vuelos que posibilitan la conexión entre islas, haciendo uso de aeronaves pequeñas por lo que se reduce el valor medio. En contraposición, destacan aeropuertos como Santiago-Rosalía de Castro o Sevilla cuyos valores medios de pasajeros por operación son altos. Se puede adelantar que para aeropuertos como los mencionados, muy probablemente se emplee de forma mayoritaria el Boeing 737 o la aeronave grande en cuestión una vez que se incluyan restricciones en el número total de operaciones por aeropuerto. 35 35 Análisis y optimización de la oferta de transporte aéreo en España Tabla 10 Total de pasajeros y operaciones y media de pasajeros por operación de los 24 aeropuertos de la red. AEROPUERTO DE ORIGEN TOTAL PASAJEROS TOTAL OPERACIONES MEDIA DE PASAJEROS/OPERACION MADRID-BARAJAS 642.914 5.144 124,9 PALMA DE MALLORCA 431.915 3.385 127,5 BARCELON-EL PRAT 589.248 4.215 139,7 GRAN CANARIA 259.446 2.479 104,6 TENERIFE NORTE 216.296 2.051 105,4 IBIZA 230.157 2.012 114,3 BILBAO 175.631 1.323 132,7 LANZAROTE 139.692 1.127 123,9 VIGO 59.669 502 118,8 A CORUÑA 46.134 338 136,4 MÁLAGA 167.119 1.175 142,2 SANTIAGO 136.322 814 167,4 SEVILLA 171.347 1.039 164,9 FUERTEVENTURA 92.216 753 122,4 TENERIFE SUR 45.015 410 109,7 MENORCA 165.331 1.644 100,5 ASTURIAS 69.730 611 114,1 VALENCIA 141.654 1.161 122 ALICANTE-ELCHE 118.785 862 137,8 SAN SEBASTIÁN 20.409 169 120,7 SANTANDER 38.316 303 126,4 GRANADA-JAÉN 37.740 296 127,5 ZARAGOZA 15.565 122 127,6 JEREZ 32.538 265 122,8 Modelo de estimación de vuelos necesarios para cubrir la demanda en una red de aeropuertos 36 Para la distancia entre aeropuertos, se ha calculado previamente la distancia ortodrómica en millas náuticas entre cada par de aeropuertos y se genera una matriz cuadrada de dimensión n, simétrica y de diagonal nula (distancia de un aeropuerto al mismo aeropuerto es nula). La distancia ortodrómica es el camino más corto posible entre dos puntos. Por lo que, en cualquier caso, la estimación de distancia entre aeropuertos será siempre la más optimista posible y no se contemplan desviaciones en la ruta del vuelo o esperas durante la aproximación. Estas cifras son difíciles o casi imposibles de estimar con precisión y suelen variar para cada operación, su desarrollo con mayor exactitud serían motivo de un estudio diferente al que se procede a realizar. Por tanto, como ya se ha dado a entender, se desestiman. [33] 3.3.4 Costes operativos En cuanto a costes de operación, se procede a enumerar el conjunto de gastos que se han tenido en cuenta en la realización de un vuelo. Todos los gastos están expresados en este apartado están referidos en euros. Los costes, están separados en dos categorías principales. Los costes fijos, que son todos aquellos en los que se incurre al realizar la operación y los costes variables que son aquellos que depende principalmente de la distancia del vuelo. A continuación, se enumera cada uno de ellos: A) Costes fijos por vuelo Para el cálculo de costes fijos se emplean tarifas reales de costes aplicados por los diferentes aeropuertos españoles estipuladas por AENA dentro del marco de tarifas aeroportuarias establecidas por la misma compañía en enero de 2024. [34] Se tienen en cuenta el coste de aterrizaje, servicios de tránsito de aeródromo y coste de pasarelas. Estos dependerán tanto del peso máximo de despegue (MTOW) de la aeronave, como la duración de la actividad en cuestión. En cada operación realizada en el modelo estas tarifas están contempladas para el aeropuerto de destino, ya que no tendría sentido que se introdujesen costes que ya se contemplaron en la anterior operación. Es decir, las operaciones de aterrizaje, tránsito por el aeródromo y estacionamiento de cada vuelo se calculan en función de las tarifas aplicadas por el aeropuerto donde se van a realizar no del que proviene la aeronave. De donde provenga la aeronave, ya se aplicaron dichas tasas en la anterior operación. El cálculo del coste monetario de cada operación se detalla en las siguientes expresiones para un vuelo realizado con origen en el aeropuerto i con destino j por la aeronave de tipo l: • Coste aterrizaje ( 𝐶𝑎" ,) : Cada aeropuerto tiene un coste específico para el aterrizaje, a esta cuantía se le denominará 𝐴" y la expresión dependerá de peso máximo de despegue de la aeronave de la siguiente forma: 𝐶𝑜𝑠𝑡𝑒*𝑎𝑡𝑒𝑟𝑟𝑖𝑧𝑎𝑗𝑒*𝑎𝑒𝑟𝑜𝑝𝑢𝑒𝑟𝑡𝑜*𝑗*𝑑𝑒*𝑙𝑎*𝑎𝑒𝑟𝑜𝑛𝑎𝑣𝑒*𝑙=𝐶𝑎" ,=*𝐴"*𝑀𝑇𝑂𝑊,* ( 3−10 ) 37 37 Análisis y optimización de la oferta de transporte aéreo en España • Coste servicios de tránsito de aeródromo ( 𝑇𝑎" ): Se calcula de forma similar a la del coste del aterrizaje aplicando una tarifa dependiente del aeropuerto de destino denominada 𝑇" : 𝐶𝑜𝑠𝑡𝑒*𝑡𝑟𝑎𝑛𝑠𝑖𝑡𝑜*𝑎𝑒𝑟ó𝑑𝑟𝑜𝑚𝑜*𝑑𝑒𝑙*𝑎𝑟𝑜𝑝𝑢𝑒𝑟𝑡𝑜*𝑗*𝑑𝑒*𝑙𝑎*𝑎𝑒𝑟𝑜𝑛𝑎𝑣𝑒*𝑙=𝑇𝑎" ,=𝑇"*𝑀𝑇𝑂𝑊,** ( 3−11 ) • Coste servicio de estacionamiento de aeronave: Esta variable depende tanto del peso máximo de despegue de la aeronave como del tiempo expresado en fracciones de 15 minutos de la prestación del servicio. El tiempo de estacionamiento de cada aeronave se estipula en la siguiente tabla: Tabla 11 Tiempos de estacionamiento en el aeropuerto de las aeronaves empleadas Modelo Tiempo de estacionamiento (minutos) Fracciones de tiempo (Ft) Embraer 190 60 4 Boeing 737-800 105 7 Se contempla el hecho de que la aeronave pequeña pueda realizar todas las operaciones de salida y entrada de pasajeros, así como de equipaje y recarga de combustibles en un tiempo algo menor que la aeronave grande. En cuanto al cálculo del coste total de estacionamiento, la expresión dada por AENA depende del aeropuerto a realizar el estacionamiento siendo de forma general la siguiente: 𝐶𝑜𝑠𝑡𝑒*𝑒𝑠𝑡𝑎𝑐𝑖𝑜𝑛𝑎𝑚𝑖𝑒𝑛𝑡𝑜*𝑒𝑛*𝑒𝑙*𝑎𝑒𝑟𝑜𝑝𝑢𝑒𝑟𝑡𝑜*𝑗*𝑑𝑒*𝑙𝑎*𝑎𝑒𝑟𝑜𝑛𝑎𝑣𝑒*𝑙=*𝐸" ,=𝑒"#𝑀𝑇𝑂𝑊,*𝐹6 ,* ( 3−12 ) Donde 𝑒"# es el coeficiente que proporciona AENA para cada aeropuerto. Para los aeropuertos Vigo, A Coruña, San Sebastián, Santander-Severiano Ballesteros, F.G.L-Granada, Zaragoza y Jerez-La Parra, la expresión anterior se calcula de una manera diferente: 𝐶𝑜𝑠𝑡𝑒*𝑒𝑠𝑡𝑎𝑐𝑖𝑜𝑛𝑎𝑚𝑖𝑒𝑛𝑡𝑜*𝑒𝑛*𝑒𝑙*𝑎𝑒𝑟𝑜𝑝𝑢𝑒𝑟𝑡𝑜*𝑗*𝑑𝑒*𝑙𝑎*𝑎𝑒𝑟𝑜𝑛𝑎𝑣𝑒*𝑙=*𝐸" ,=𝑒" 96*𝑀𝑇𝑂𝑊,*𝐹6 ,* ( 3−13 ) • Coste de usos de pasarelas de embarque: El uso de pasarelas de embarque es en muchos casos opcional. En nuestro modelo su expresión es similar a la necesaria para calcular el coste de estacionamiento, pero no depende del peso máximo de despegue de la aeronave. Donde 𝑝" será el coeficiente proporcionado por AENA. Modelo de estimación de vuelos necesarios para cubrir la demanda en una red de aeropuertos 44 3.4 Principales limitaciones del modelo Este modelo se enfrenta a limitaciones importantes que afectan su precisión y aplicabilidad en entornos reales y que en ningún momento se pueden obviar, a continuación, se enumeran diferentes aspectos en los que se ha tenido que realizar simplificaciones que limitan al modelo. Una de las limitaciones fundamentales es la falta de información detallada sobre los viajes de los pasajeros. En nuestro caso al solo conocer como los pasajeros vuelan de forma directa de un aeropuerto a otro, es decir, si un pasajero vuela de un destino a otro mediante escala en un tercero, en nuestro modelo dicho pasajero aparecerá como dos pasajeros, uno que vuela del origen al lugar de escala y del lugar de escala al destino. Otra limitación de este modelo es la restricción del número de aeronaves utilizadas en la optimización. Aunque en la realidad las aerolíneas suelen tener flotas de aeronaves más diversas, con diferentes capacidades y costes operativos, el modelo desarrollado simplifica esta variable a solo dos posibilidades. Esta simplificación puede no ser excesivamente fiel a la realidad de la gestión de flotas que las aerolíneas disponen. En todo caso la resolución del problema de optimización incluyendo un mayor número de aeronaves requeriría de una potencia computacional de la cual no se dispone. Un modelo de mayor calidad debería contemplar al menos una aeronave aún menor a las empleadas como podría ser un ATR-72 de unas 70 plazas. Amén de diferentes modelos en el rango de 140-230 pasajeros como podrían ser las aeronaves Airbus A-319 o Airbus A-321 de menor y mayor capacidad que el Boeing 737-800 que se emplea. En cuanto a los costes asociados con los vuelos, el modelo de optimización utiliza un enfoque de coste por asiento-milla. Si bien esto puede ser útil para simplificar los cálculos, no refleja la realidad de los costes variables asociados con diferentes tipos de vuelos. Ignorar esta variabilidad puede llevar a decisiones subóptimas en la asignación de rutas y aeronaves. En el caso desarrollado en este trabajo, solo se contempla un valor de Coste por asiento-milla disponible único para cada aeronave. Además, estos valores empleados son en todo caso una aproximación a un posible valor real basado en algunos casos en datos antiguos y en otros a valores aproximado de variables como el consumo de las aeronaves o costes de tripulación. No se contemplan otros costos en los que recaen las aerolíneas a la hora de realizar un vuelo como catering, servicios en aeropuertos para pasajeros, etc. El modelo generado se centra en los costos operativos de un vuelo, excluyendo variables anteriormente mencionadas. Por último, es importante destacar la que probablemente sea la mayor limitación del modelo. A la hora de estimar que aeronave es la más optima a emplear el modelo solo calcula cual sería esta. A la hora de planificar las operaciones aéreas, es fundamental tener en cuenta la disponibilidad real de aviones en los aeropuertos de donde proceden. Este modelo no considera la disponibilidad de aviones ya que supone que los mejores aviones siempre están listos para volar y hacer su trabajo. Pero en la vida real, conseguir un avión en un determinado aeropuerto puede ser complicado debido a cosas como arreglar el avión, esperarlo, cambiar el horario o problemas inesperados que surgen. No se tiene en cuenta la disponibilidad real de aviones ni en tiempo ni en lugar. Por ejemplo, el modelo podría sugerir el uso de una aeronave en una ruta determinada debido a su eficiencia o capacidad, pero si esa aeronave no está disponible en el aeropuerto de origen a la hora requerida, la operación real podría enfrentar contratiempos temporales o no poder realizar la operación estipulada como óptima. Solo se selecciona cual es la aeronave óptima para realizar una operación dentro de las restricciones selecionadas. En resumen, la falta de información detallada sobre los viajes de los pasajeros, la simplificación de la flota de aeronaves, las simplificaciones realizadas en el los cálculos de los costes de operación y omisión de diferentes costes de las aerolíneas deben considerarse para desarrollar modelos más precisos y efectivos que reflejen con mayor precisión la complejidad del sector de la aviación comercial. El modelo creado solo busca dar 45 45 Análisis y optimización de la oferta de transporte aéreo en España una primera aproximación a como debería operarse una red para que esto se realice de la manera más optima posible a partir de los datos de los cuales se tiene disposición. 47 4 RESULTADOS DEL MODELO n este apartado, se van a discutir el conjunto de resultados obtenidos para simulaciones de diferente naturaleza. En primer lugar, se hará un conjunto de simulaciones que nos permitan conocer con mayor profundidad cómo funciona el modelo. Por un lado, es interesante conocer qué tipo de aeronave es óptima emplear en función de la distancia a partir de los datos anteriormente mostrados. Posteriormente, se realizarán una serie de simulaciones en la que se intentará justificar la elección del coste por asiento disponible por milla empleado. Esto permitirá extraer conclusiones, para más adelante crear un modelo con un mayor número de restricciones que se asemejen más a lo que realmente ocurre en la realidad. A continuación, se realizarán diferentes simulaciones para variados escenarios y usos de aeronaves haciendo uso de diversas restricciones, que se detallarán en cada apartado, estas serán usando datos de número de pasajeros reales y de operaciones realizadas, para el mes de agosto de 2022. Como análisis previo a la realización de cualquier simulación, es conveniente destacar como se organizan los datos de la muestra. En la figura 19 se detalla a modo de histograma el ordenamiento del número de pasajeros reales entre cada par origen-destino para el mes de alta demanda. Podemos observar cómo el conjunto de orígenes-destinos entre los que no hay pasajeros supone un 37.86% del total de los 552 pares (habiéndose eliminado con anterioridad los 24 pares que tendrían el mismo aeropuerto origen y destino). En general nos enfrentamos a una muestra muy desigual. En el mes de agosto parecen coexistir rutas altamente demandas con un número mayor de rutas que no se emplean. Esto tiene mayor sentido después de haber realizado el análisis del reparto de pasajeros entre los diferentes aeropuertos españoles. Si, por ejemplo, nos encontrásemos en un mercado más uniforme como Estados Unidos, la siguiente figura tendría un aspecto mucho más uniforme, desapareciendo gran parte de las rutas sin emplear. Figura 19: Número de rutas en función de su número de pasajeros en el mes de agosto de 2022 Número de pares origen-destino con los siguientes rangos de pasajeros 209 156 107 50 16 14 Nulos [1,5000] [5000,15000] [15000,30000] [30000,60000] >60000 Rangos de demanda 0 50 100 150 200 250 Valores E Resultados del modelo 48 48 4.1 Uso de la distancia como medida de selección de aeronaves El principal objetivo de este apartado es conocer que aeronave es, según el modelo, la más optima en función del coste de operación en los diferentes rangos de distancias que hay entre todos los aeropuertos. Se busca que el modelo tenga la mayor libertad de decisión a la hora de escoger la aeronave pues solo queremos conocer lo anteriormente nombrado y cualquier tipo de restricción operativa o de distancia para el uso de una aeronave condicionaría los resultados. Es por ello, por lo que, en este apartado solo se impondrá la restricción de capacidad de los vuelos (ecuación 3-4), que solo impone que la capacidad del total de vuelos necesarios sea mayor a la demanda real. En la figura 20 puede observarse el porcentaje de empleo de la aeronave grande en función de la distancia, es decir, de todas las rutas para las que el modelo ha calculado que al menos ser requiere una operación, que porcentaje de estas se realiza con la aeronave grande. Se puede interpretar fácilmente, que, para distancias cortas, el uso de la aeronave pequeña es predomínate. Como ya se adelantó anteriormente es el rango de operación óptimo que se estima para este tipo de aeronaves. A medida que la distancia de la ruta aumenta, su uso disminuye, prevaleciendo la aeronave grande. Figura 20: Porcentaje de vuelos realizado por la aeronave grande en cada ruta en función de la distancia Si se examina la figura 20 se observa que numerosos puntos se encuentran en un rango intermedio. Esto se debe principalmente a que muchas rutas operadas poseen baja demanda, es decir, el número de pajeros que desea realizar esa ruta es pequeño. La modelo está planteada de manera que debe dar cabida a todos los pasajeros existentes, como modela la restricción de capacidad de vuelo. Es preferible emplear algún resultado real para realizar una mejor explicación, la ruta entre Jerez y Santander en el mes de agosto de 2022 fue empleada por 696 pasajeros. Según el modelo la aeronave grande sería la más 0 500 1000 1500 2000 2500 Distancia entre aeropuertos 0 10 20 30 40 50 60 70 80 90 100 % de vuelos del avión grande Porcentaje de vuelos del avión grande vs. Distancia entre aeropuertos 49 49 Análisis y optimización de la oferta de transporte aéreo en España barata de emplear para cubrir esa ruta, sin embargo, los resultados que arroja es que deberían realizarse un total de 5 operaciones. Tres con el Boeing 737-800 y dos con el Embraer 190, esto es porque al ser el número de pasajeros totales no múltiplo de 180, puede que para acoger a los pasajeros sobrantes sea más rentable emplear la aeronave pequeña, o encontrar una proporción de uso de ambas aeronaves que minimice los pasajeros sobrantes y así minimizar costes. Si se atiende al número de vuelos para cubrir dicha ruta, se obtiene que solo el 60% de los vuelos serían con el Boeing 737 y el 40% por el Embraer 190. Es decir, no se usa exclusivamente una aeronave. Este fenómeno ocurre en la mayoría de las rutas, pero es más visible a medida que el número de operaciones es menor. Por poner otro ejemplo en la ruta que une Madrid-Barajas con Tenerife Norte se realizan 341 vuelos en total, 335 con la aeronave grande (98.24%) y 6 con la pequeña (1.76%). El rango en el cual, la aeronave pequeña es más rentable que la grande, es variable en función del aeropuerto al que se dirija el vuelo. Como se explicó en el apartado de datos, cada aeropuerto tiene sus tasas propias para el conjunto de operaciones que se realiza en él y como en el modelo creado, estas tasas se aplican según el aeropuerto de destino. Podemos calcular la distancia exacta a la cual la aeronave grande pasa a ser más rentable según las diferentes tarifas aeroportuarias. En la figura 21, se muestra este fenómeno para todos los aeropuertos con diferentes tasas. En la imagen se puede observar que, el rango de puntos de inflexión es grande. Desde 350 km para los aeropuertos en los cuales las tarifas son más baratas (Jerez de la Frontera, Zaragoza, Granada, etc.); hasta los 750 km para el caso del aeropuerto de Madrid-Barajas, que es el aeropuerto español con tarifas más caras. En general, para muchos de los aeropuertos peninsulares sería más económico emplear la aeronave pequeña para rutas con el resto de la península y Baleares y para los aeropuertos canarias el uso de la aeronave grande para realizar conexiones con la península y el pequeño para vuelos interinsulares. Figura 21: Puntos a partir del cual la aeronave grande pasa a ser más rentable de operar en función de la distancia para cada aeropuerto según sus tarifas aeroportuarias. Podemos concluir que a mayor sean los costos de las tasas aeroportuarias del aeropuerto de destino, mayor será la distancia en la que la aeronave pequeña será más rentable de operar. Generalmente estas tarifas son mayores en los grandes aeropuertos que en los que presentan menos tráfico. Esto presenta un conflicto para el modelo pues sin ningún tipo de restricción adicional, en los aeropuertos grandes las aeronaves pequeñas son Resultados del modelo 50 50 las que llegarían en mayor medida. Este suceso no es deseado pues puede provocar que aeropuertos como Madrid-Barajas o Barcelona reciban un número de vuelos muy superior a los que realmente ocurren en la realidad. No obstante, si se realizan sendas simulaciones en las que la variable CASM se modifica, podríamos observar cuál serían los valores óptimos para emplear de forma relativa. Sírvase este apartado como forma de justificar los valores Coste por Asiento-Milla Disponible que se han empleado para futuras simulaciones. En la figura 22 se detalla la variación del uso porcentual de la aeronave grande en función de la variación de del CASM de la aeronave pequeña con la distancia. En la imagen de arriba-izquierda donde el CASM de la aeronave pequeña es 0.0527 y el de la aeronave grande es 0.603 (es decir el de la aeronave pequeña mucho menor que el de la grande), el uso de la aeronave grande solo se da para pequeñas distancias. Siguiendo el orden arriba-derecha y abajo-izquierda, el Coste por Asiento- Milla Disponible de la aeronave pequeña aumenta respecto a la de la grande hasta hacerse igual. Mientras que el CASM de la pequeña sea menor que el de la grande su uso se dará principalmente en las grandes distancias y cuando se igualan comienza a prevalecer la grande. No es hasta que el CASM de la pequeña es mayor (por cierto margen) que el de la grande cuando se revierten sus funciones. Es decir, aparece el fenómeno buscado y cada aeronave se emplea en su rango óptimo de operación de forma aproximada. Figura 22: Variación del porcentaje de empleo de la aeronave grande en función de la distancia de la ruta para diferentes valores relativos del CASM de cada aeronave. 0 500 1000 1500 2000 2500 Distancia entre aeropuertos 0 10 20 30 40 50 60 70 80 90 100 % de vuelos del avión grande Porcentaje de vuelos del avión grande vs. Distancia entre aeropuertos: CASM EMBRAER-E190=0.0527, CASM BOEING 737=0.0603 0 500 1000 1500 2000 2500 Distancia entre aeropuertos 0 10 20 30 40 50 60 70 80 90 100 % de vuelos del avión grande Porcentaje de vuelos del avión grande vs. Distancia entre aeropuertos: CASM EMBRAER-E190=0.0573, CASM BOEING 737=0.0603 0 500 1000 1500 2000 2500 Distancia entre aeropuertos 0 10 20 30 40 50 60 70 80 90 100 % de vuelos del avión grande Porcentaje de vuelos del avión grande vs. Distancia entre aeropuertos: CASM EMBRAER-E190=0.0603, CASM BOEING 737=0.0603 0 500 1000 1500 2000 2500 Distancia entre aeropuertos 0 10 20 30 40 50 60 70 80 90 100 % de vuelos del avión grande Porcentaje de vuelos del avión grande vs. Distancia entre aeropuertos: CASM EMBRAER-E190=0.0672, CASM BOEING 737=0.063 51 51 Análisis y optimización de la oferta de transporte aéreo en España 4.2 Modelo de alta demanda aplicando restricciones operativas En el apartado anterior se han realizado sendas simulaciones muy poco restringidas. Puesto que se quería comprobar que aeronaves es la más adecuada en cada tipo de operación. No sin más, dichos modelos carecen en muchos sentidos de viabilidad operacional. Para poner en contexto este fenómeno, en la figura 23 se muestra el conjunto de operaciones que se debieron realizar en el mes de agosto para cubrir la demanda real que se produjo según el modelo anterior para el aeropuerto de Madrid-Barajas. Figura 23: Resultados obtenidos de los vuelos necesarios a realizar con un modelo sin limitaciones de operación para el aeropuerto de Madrid-Barajas. Si se observa con detenimiento la imagen se ve que para la mayoría de los destinos se deben realizar un número de operaciones muy elevado. Principalmente por el hecho de que dicho aeropuerto se encuentra en el interior de una península, y exceptuando los aeropuertos canarios y los más alejados de Madrid en la península, la mayoría de los aeropuertos se encuentran en el rango de distancias que es óptimo emplear la aeronave pequeña. Esto se ve con mayor claridad observando el caso de Palma de Mallorca, para el cuál según el modelo serían necesarias más de mil operaciones solo para este destino. Si se sumasen todas las operaciones que el modelo estima realizar serían un total de 6466, mientras que, si observamos los datos estadísticos reales de dicho mes, desde el aeropuerto de Madrid-Barajas al resto de aeropuertos de la red se realizaron un total exacto de 5144 operaciones. Esto supone un aumento aproximado del 25% de operaciones a realizar. Muchos aeropuertos españoles por dichas fechas tienen grandes volúmenes de operaciones y asumir un aumento tan brusco no sería una opción. Es por ello por lo que en este apartado se va a desarrollar un modelo más restringido, en especial en el tema operacional, para intentar que, este que es el modelo final presentado en este trabajo, arroje resultados viables operacionalmente. Aeropuerto MADRID-BARAJAS ADOLFO SUÁREZ 371 329 168 140 129 119 94 174 45 43 902 630 580 206 253 196 139 123 134 121 87 43 MAD PMI BCN LPA TFN IBZ BIO ACE VGO LCG AGP SCQ SVQ FUE TFS MAH OVD VLC ALC EAS SDR GRX ZAZ XRY Aeropuerto destino 0 100 200 300 400 500 600 700 800 900 1000 nº de vuelos Embraer 190 Boeing 737-800 Resultados del modelo 52 52 Las restricciones que se van a emplear están diseñadas buscando un consenso entre guiar al modelo para que arroje resultados de interés. Pero intentar dar la máxima flexibilidad a este para que los valores obtenidos no se vean creados de forma artificial. Las restricciones adicionales que se van a emplear son las siguientes: • Máximo número de operaciones de cada aeropuerto: Esta restricción estará formulada a su vez por dos restricciones, una fija para cada aeropuerto y otra dependiente del valor real de la demanda. En el apartado de desarrollo del modelo se definió esta restricción como de tipo mixto. Por un lado, se procede a generarse una restricción para cada aeropuerto (en total 24), la cual limita el número máximo de operaciones de salida de cada aeropuerto al total de operaciones que se realizaron en el mes de agosto de 2022 al conjunto del resto de aeropuertos de la red. Esta restricción garantiza que no se produzca un exceso de operaciones en determinadas rutas como ocurre en la figura 23. Pues si existe una solución viable esta debe ser menor o igual al dato real. Sin embargo, esta restricción todavía permite al modelo mucho margen en cuanto a las posibles opciones. Puede darse el caso de un aeropuerto que opte por el uso de la aeronave pequeña para una ruta y deban hacerse un gran número de operaciones, pero para el resto de las rutas opte por la grande. Este fenómeno no es deseable en ningún caso pues se busca que ninguna ruta tengo un número excesivamente elevado de operaciones. Para evitar este suceso se incluye una restricción adicional que solo afecta a las rutas de la red de mayor demanda. Se estipula que para el modelo de dicho mes una gran demanda, es toda aquella ruta, la cual tenga un total mayor a 45000 pasajeros para el mencionado periodo de tiempo. Si la ruta cumple la condición, el número de operaciones de esta ruta se limita a que pueda ser como muy grande igual al que ocurrió en la realidad. Esta doble condición puede generar resultados de gran interés, haciendo un balance entre libertad de decisión y restricción. Otra opción contemplable, seria restringir el número de operaciones de cada ruta, no solo de las de gran demanda si no para todas. Esta restricción no es completamente aconsejable pues al ser mucho más restrictiva los resultados pueden carecer de genuinidad. Dicho todo lo anterior, las restricciones descritas por las ecuaciones 3-8 tendrían la siguiente forma: # (𝑦!" 780/.1/ +*𝑦!" 9:;5;) 2 "&' ≤* # 𝑣!" 2 !&' 𝑦!" 780/.1/ +*𝑦!" 9:;5; ≤𝑣!"****𝑝𝑎𝑟𝑎*𝑓!" >45000*𝑝𝑎𝑥* ( 4−2 ) 53 53 Análisis y optimización de la oferta de transporte aéreo en España • Restricción del uso de la aeronave de mayor tamaño para la realización de vuelos de muy corta duración. Es extraño que una ruta de muy corta duración en vuelo se empleen aeronaves de mayor tamaño a los ATR-72 o en todo caso algo mayores, pero no demasiado. Carecería de sentido emplear una aeronave como podría ser un Boeing-737, Airbus A-320 o mayor. Es por ello por lo que, en el modelo se va a incluir una restricción adicional. Esta restricción elimina la posibilidad de uso de la aeronave de mayor tamaño en este tipo de vuelos tan cortos. Por lo tanto, para una distancia entre aeropuertos menor a una distancia umbral que delimitaremos a 270 kilómetros se elimina la posibilidad del uso de una aeronave que no sea la más pequeña que contemple el modelo. Esta distancia garantiza que dicha restricción actúe sobre todas las rutas que conectan todos los aeropuertos de las Islas Canarias que seleccionaremos se elimina la posibilidad del uso de una aeronave que no sea la más pequeña que contemple el modelo. Esta restricción tiene a su vez una excepción para el caso de las rutas que unen el aeropuerto de Barcelona-El Prat con las Islas Baleares (Palma de Mallorca, Menorca e Ibiza) pues son tres rutas que, encontrándose dentro del rango de 270 kilómetros, no son viables de operar haciendo uso exclusivo del Embraer-190 ya que su demanda de pasajeros en meses estivales es muy alta. Esta restricción basada en la ecuación 3-9 quedaría de la siguiente manera: 𝑦!" 9:;5; =0********∀*𝑟𝑢𝑡𝑎*𝑐𝑜𝑛*𝑑!" <270*𝑘𝑚* ( 4−3 ) Una vez se simula el modelo, es preferible comenzar comprobando las principales diferencias entre este y el anterior modelo. La manera más fácil de apreciar los cambios sería viendo cómo se reestructura las operaciones en los tres principales aeropuertos del país. La figura 24 muestra el número de vuelos a realizar con el nuevo modelo que incluye todas las restricciones planteadas (derecha) y los resultados que se obtendrían sin aplicar dichas restricciones (izquierda). Todas las figuras tienen la misma escala para facilitar ver las diferencias. Es posible percatarse de diferentes fenómenos. Por un lado, el número de operaciones en los aeropuertos disminuye considerablemente. En el conjunto completo de la red, estas restricciones consiguen una disminución total del 16.35% de operaciones necesarias. Si nos situamos en el contexto de cada aeropuerto, para el caso de Madrid-Barajas se eliminan 760 operaciones lo que supone una disminución del 12.71%, para el caso de Palma de Mallorca otras 838 operaciones o un 21.16% y para el caso de Barcelona-El Prat un total de 1144 operaciones o lo que es lo mismo un 20.96% de las operaciones totales del aeropuerto para dicho mes. Esta disminución se produce principalmente por el cambio en el uso de la aeronave pequeña por la aeronave grande para poder cumplir el máximo de operaciones posibles. Este cambio se da principalmente en los aeropuertos peninsulares y baleares pues la mayoría de sus destinos se encuentran dentro del rango óptimo de la aeronave pequeña. En la figura 24 se puede observar gráficamente este fenómeno, donde las grandes barras naranjas (Operaciones aeronave pequeña), se convierten en barras predominantemente azules (Operaciones aeronave grande). El modelo lo que realmente busca es emplear la aeronave pequeña pero la limitación de operaciones le obliga a emplear en mayor medida la grande. Resultados del modelo 60 60 Figura 30 Resultados obtenidos para los aeropuertos Bilbao y Lanzarote-César Manrique En las figuras 27 y 28 se muestran los resultados para algunos de los mayores aeropuertos de la red. Destaca por un lado el gran número de rutas que se combina el uso de ambas aeronaves, debido principalmente a que en estas son en las que se aplica la restricción de operaciones. Se emplea en Embraer 190 solo hasta que el número de operaciones lo permita. Puede observase este fenómeno en rutas como las que unen las Islas Baleares con Barcelona y Madrid. Por otro lado, Madrid-Barajas es el aeropuerto de la red en el que más conexiones y en un volumen de vuelos considerable tiene dentro de la red. También se puede destacar la gran polaridad en el uso de ambas aeronaves en los destinos del aeropuerto de Gran Canaria. Estas rutas o son muy cortas y se fuerza a el uso del Embraer 190 o son muy largas para entrar en el rango óptimo de operación de esta aeronave. Los aeropuertos de Tenerife Norte y Lanzarote tienen un comportamiento similar al que encontrábamos en el caso de Gran Canaria. A medida que aparecen aeropuertos menores, se observa como el conjunto de rutas que opera cada aeropuerto va disminuyendo. Destaca en ese sentido el aeropuerto de Ibiza, en el que encontramos por un lado rutas con una gran cantidad de operaciones (Madrid, Barcelona, Palma de Mallorca o Valencia), mientras que gran cantidad de posibles destinos no ofrece conexión. En este sentido el aeropuerto de Bilbao también parece tener un carácter muy polarizado la elección de aeronave. O se emplea una para cada ruta o la otra. Esto se debe principalmente a que, a pesar de ser un aeropuerto con un volumen de tráfico considerable, no posee ninguna ruta de las consideradas como de alta demanda (más de 45.000 pasajeros en el mes). Es por ello por lo que por la naturaleza del modelo existe a la práctica una casi libre elección de la aeronave a emplear. Aeropuerto BILBAO 44 52 65 51 78 85 87 216 199 185 51 117 MAD PMI BCN LPA TFN IBZ BIO ACE VGO LCG AGP SCQ SVQ FUE TFS MAH OVD VLC ALC EAS SDR GRX ZAZ XRY Aeropuerto destino 0 50 100 150 200 250 nº de vuelos Boeing 737 Embraer 190 Aeropuerto LANZAROTE CÉSAR MANRIQUE 174 67 50 36 357 229 MAD PMI BCN LPA TFN IBZ BIO ACE VGO LCG AGP SCQ SVQ FUE TFS MAH OVD VLC ALC EAS SDR GRX ZAZ XRY Aeropuerto destino 0 50 100 150 200 250 300 350 400 nº de vuelos Boeing 737 Embraer 190 61 61 Análisis y optimización de la oferta de transporte aéreo en España Figura 31 Resultados obtenidos para los aeropuertos de Vigo y A Coruña Figura 32 Resultados obtenidos para los aeropuertos de Málaga-Costa del Sol y Santiago-Rosalía de Castro Aeropuerto VIGO 76 308 MAD PMI BCN LPA TFN IBZ BIO ACE VGO LCG AGP SCQ SVQ FUE TFS MAH OVD VLC ALC EAS SDR GRX ZAZ XRY Aeropuerto destino 0 50 100 150 200 250 300 350 nº de vuelos Boeing 737 Embraer 190 Aeropuerto A CORUÑA 62 65 162 MAD PMI BCN LPA TFN IBZ BIO ACE VGO LCG AGP SCQ SVQ FUE TFS MAH OVD VLC ALC EAS SDR GRX ZAZ XRY Aeropuerto destino 0 50 100 150 200 250 nº de vuelos Boeing 737 Embraer 190 Aeropuerto MÁLAGA-COSTA DEL SOL 115 202 55 50 62 77 48 295 26 MAD PMI BCN LPA TFN IBZ BIO ACE VGO LCG AGP SCQ SVQ FUE TFS MAH OVD VLC ALC EAS SDR GRX ZAZ XRY Aeropuerto destino 0 50 100 150 200 250 300 nº de vuelos Boeing 737 Embraer 190 Aeropuerto SANTIAGO-ROSALÍA DE CASTRO 83 59 103 63 35 37 51 52 36 45 59 80 MAD PMI BCN LPA TFN IBZ BIO ACE VGO LCG AGP SCQ SVQ FUE TFS MAH OVD VLC ALC EAS SDR GRX ZAZ XRY Aeropuerto destino 0 20 40 60 80 100 120 140 160 180 nº de vuelos Boeing 737 Embraer 190 Resultados del modelo 62 62 En la figura 31, aparece un fenómeno que repetirá con mayor frecuencia a medida que se muestren los resultados de aeropuertos con menor volumen de tráfico. Los aeropuertos de Vigo y A Coruña, solo poseen conexiones de interés con especialmente el aeropuerto de Madrid-Barajas y en menor medida con Barcelona- El Prat. Como se vio en el apartado 2.3, en la figura 7, ambos son de los aeropuertos con mayor coeficiente de Gini dentro del mercado nacional. Por otro lado, en la misma comunidad, tenemos el caso de Santiago. Este aeropuerto es el principal dentro de Galicia y a diferencia de sus respectivos, este si posee conexiones importantes con la mayor parte de los aeropuertos nacionales. Al encontrarse en la zona más nor-este del país, muchos de los aeropuertos peninsulares como Sevilla, Málaga, Valencia, etc.; se encuentran ya fuera del rango de operación óptimo del Embraer 190. Solo se exceptúa el caso de Madrid donde si actúa la restricción por alta demanda. El aeropuerto de Málaga-Costa del Sol si es uno de los principales aeropuertos nacionales, aunque este más enfocado al tráfico internacional (Alemania e Inglaterra) posee numerosas conexiones dentro del país. El comportamiento de los resultados es similar al que aparece en el caso de Santiago-Rosalía de Castro. Figura 33 Resultados obtenidos para los aeropuertos de Sevilla y Fuerteventura Aeropuerto SEVILLA 125 196 68 74 54 82 55 39 149 60 MAD PMI BCN LPA TFN IBZ BIO ACE VGO LCG AGP SCQ SVQ FUE TFS MAH OVD VLC ALC EAS SDR GRX ZAZ XRY Aeropuerto destino 0 20 40 60 80 100 120 140 160 180 200 nº de vuelos Boeing 737 Embraer 190 Aeropuerto FUERTEVENTURA 123 40 274 186 MAD PMI BCN LPA TFN IBZ BIO ACE VGO LCG AGP SCQ SVQ FUE TFS MAH OVD VLC ALC EAS SDR GRX ZAZ XRY Aeropuerto destino 0 50 100 150 200 250 300 nº de vuelos Boeing 737 Embraer 190 63 63 Análisis y optimización de la oferta de transporte aéreo en España Figura 34 Resultados obtenidos para los aeropuertos de Tenerife Sur y Menorca Figura 35 Resultados obtenidos para los aeropuertos de Asturias y Valencia Aeropuerto TENERIFE SUR 95 35 79 30 MAD PMI BCN LPA TFN IBZ BIO ACE VGO LCG AGP SCQ SVQ FUE TFS MAH OVD VLC ALC EAS SDR GRX ZAZ XRY Aeropuerto destino 0 20 40 60 80 100 120 nº de vuelos Boeing 737 Embraer 190 Aeropuerto MENORCA 142 88 312 133 415 45 114 43 MAD PMI BCN LPA TFN IBZ BIO ACE VGO LCG AGP SCQ SVQ FUE TFS MAH OVD VLC ALC EAS SDR GRX ZAZ XRY Aeropuerto destino 0 100 200 300 400 500 600 nº de vuelos Boeing 737 Embraer 190 Aeropuerto ASTURIAS 38 33 28 39 32 135 118 MAD PMI BCN LPA TFN IBZ BIO ACE VGO LCG AGP SCQ SVQ FUE TFS MAH OVD VLC ALC EAS SDR GRX ZAZ XRY Aeropuerto destino 0 20 40 60 80 100 120 140 nº de vuelos Boeing 737 Embraer 190 Aeropuerto VALENCIA 32 39 29 45 33 62 155 411 30 177 MAD PMI BCN LPA TFN IBZ BIO ACE VGO LCG AGP SCQ SVQ FUE TFS MAH OVD VLC ALC EAS SDR GRX ZAZ XRY Aeropuerto destino 0 50 100 150 200 250 300 350 400 450 nº de vuelos Boeing 737 Embraer 190 Resultados del modelo 64 64 Figura 36 Resultados obtenidos para los aeropuertos de Alicante-Elche Miguel Hernández y San Sebastián Figura 37 Resultaos obtenidos para los aeropuertos de Santander-Severiano Ballesteros y Granada-Jaén Aeropuerto ALICANTE-ELCHE MIGUEL HDEZ. 89 38 66 58 28 35 135 244 105 MAD PMI BCN LPA TFN IBZ BIO ACE VGO LCG AGP SCQ SVQ FUE TFS MAH OVD VLC ALC EAS SDR GRX ZAZ XRY Aeropuerto destino 0 50 100 150 200 250 nº de vuelos Boeing 737 Embraer 190 Aeropuerto SAN SEBASTIÁN 87 MAD PMI BCN LPA TFN IBZ BIO ACE VGO LCG AGP SCQ SVQ FUE TFS MAH OVD VLC ALC EAS SDR GRX ZAZ XRY Aeropuerto destino 0 10 20 30 40 50 60 70 80 90 nº de vuelos Boeing 737 Embraer 190 Aeropuerto SANTANDER-SEVE BALLESTEROS 97 58 MAD PMI BCN LPA TFN IBZ BIO ACE VGO LCG AGP SCQ SVQ FUE TFS MAH OVD VLC ALC EAS SDR GRX ZAZ XRY Aeropuerto destino 0 10 20 30 40 50 60 70 80 90 100 nº de vuelos Boeing 737 Embraer 190 Aeropuerto F.G.L. GRANADA - JAÉN 65 50 76 50 MAD PMI BCN LPA TFN IBZ BIO ACE VGO LCG AGP SCQ SVQ FUE TFS MAH OVD VLC ALC EAS SDR GRX ZAZ XRY Aeropuerto destino 0 20 40 60 80 100 120 nº de vuelos Boeing 737 Embraer 190 65 65 Análisis y optimización de la oferta de transporte aéreo en España Figura 38 Resultados obtenidos para los aeropuertos de Zaragoza y Jerez de la Frontera-La Parra En las figuras 31,32 y 33 se muestran los resultados obtenidos para la mayoría de los aeropuertos menores considerados en la red. En ellos como ya se comentó para los casos de Vigo y A coruña apenas se ofrecen conexiones. Destaca en contraparte, el caso del aeropuerto de Santander-Severiano Ballesteros que, aunque requiere de pocos vuelos en conjunto si se operan un número considerable de rutas a más de la mitad de los aeropuertos de la red. En este sentido, otros aeropuertos poseen un carácter similar como Alicante-Elche o Valencia, estos aeropuertos si poseen un mayor volumen de tráfico internacional a diferencia de Santander, pero dentro del mercado nacional ofrecen numerosas rutas, pero con pocos vuelos. Aeropuerto ZARAGOZA 64 MAD PMI BCN LPA TFN IBZ BIO ACE VGO LCG AGP SCQ SVQ FUE TFS MAH OVD VLC ALC EAS SDR GRX ZAZ XRY Aeropuerto destino 0 10 20 30 40 50 60 70 nº de vuelos Boeing 737 Embraer 190 Aeropuerto JEREZ DE LA FRONTERA/ LA PARRA 28 64 87 MAD PMI BCN LPA TFN IBZ BIO ACE VGO LCG AGP SCQ SVQ FUE TFS MAH OVD VLC ALC EAS SDR GRX ZAZ XRY Aeropuerto destino 0 10 20 30 40 50 60 70 80 90 nº de vuelos Boeing 737 Embraer 190 Resultados del modelo 66 66 4.2.4 Frecuencias necesarias para cubrir la demanda Otra forma de analizar los resultados anteriormente mostrados podría ser el cálculo de qué frecuencias de vuelos serían necesarias para suplir la demanda dada según la aeronave empleada. Se van a estudiar dos casos principales, por un lado, para aquellas rutas que poseen un número necesario de vuelos, cuáles serían las frecuencias diarias necesarias (rutas más demandadas) y por otro lado para aquellos que no hay demanda suficiente cuantas frecuencias semanales, siempre menores a un vuelo diario son necesarias (rutas menos demandas pero que si permiten realizar vuelos de forma periódica). Dado el gran número de variables a analizar, y en pro de una mayor facilidad de comprensión de los resultados, se han generado sendas figuras. En estas, se muestran las frecuencias calculadas haciendo uso de mapas de calor. A grandes rasgos y para facilitar su comprensión, los valores en blanco son nulos y posteriormente los valores se distribuyen de forma entera entre 1 (Azul) a el máximo valor calculado (Rojo oscuro). Para el escenario de frecuencias diarias, basta que, a partir de los cálculos ya realizados, dividir el número total de vuelos necesarios de cada tipo de aeronave entre el total de los 31 días del mes de agosto, se han eliminado resultados menores de 0.5 y siempre se realiza un redondeo hacia arriba a partir de 0,5 que puede ayudar a obtener cifras de frecuencias que no están totalmente ajustadas al valor exacto. Este fenómeno en la realidad podría considerarse como la existencia de algunas plazas vacías en los distintos vuelos ofertados. En la figura 39 se muestran los cálculos realizados para el número de frecuencias diarias necesarias, tanto de la aeronave de mayor y menor tamaño. Las columnas verticales muestran el aeropuerto de origen (Representado por su código IATA de aeropuerto) y los destinos los elementos de cada fila. Las rutas de mayor demanda se dan entre los aeropuertos de Madrid-Barajas (MAD) y Barcelona-El Prat (BCN) con el aeropuerto de Palma de Mallorca. Siendo también las rutas de la mayoría de grandes aeropuertos con destinos vacacionales como pueden ser los de la mayoría de las islas o aeropuertos insulares como Valencia (VLC) y Málaga-Costa del Sol (AGP). Si nos centramos, por ejemplo, en la ruta entre Madrid-Barajas y Palma de Mallorca, esta requeriría unos 16 vuelos diarios de la aeronave grande o 112 vuelos semanales. Si compramos estos resultados con el registro de vuelos realizado en el mismo mes de agosto de 2022 (haciendo uso de la página Flightaware) [42], dependiendo del día se realizaron entre 17 y 22 vuelos diarios dependiendo del día de la semana. Dichas operaciones fueron realizadas en su mayoría por aeronaves como Boeing 737-800, Airbus A321 Neo o Airbus A319. Podemos concluir que para esta ruta la estimación es relativamente cercana a la realidad, a pesar de que se use en mayor medida una aeronave más pequeña de las que se usan realmente. Otra ruta como Barcelona-El Prat con Madrid-Barajas el modelo estipula que serían necesarias 3 vuelos de la aeronave grande y 14 con la pequeña diarios, es decir 17 diarios o 119 semanales en total. Realizando la misma comparación que en el caso anterior, esta ruta estuvo cubierta dependiendo del día entre 18 a 21 vuelos diarios. Por lo que nuestra estimación parece realizar de forma generalizada un menor número de vuelos a los que se realizaron en la realidad. 67 67 Análisis y optimización de la oferta de transporte aéreo en España Figura 39 Frecuencias diarias necesarias para cubrir la demanda dada en aquellas rutas que permitan tener al menos un vuelo diario. No obstante, existe un gran número de rutas que no requieren de si quiera un vuelo diario para satisfacer la demanda. En la figura 40 se muestran el número de frecuencias semanales necesarias para cada una de las rutas que poseen demanda suficiente para realizar de una a seis operaciones semanales. Figura 40 Número de frecuencias semanales necesarias para cubrir la demanda dada para aquellas rutas que solo requieren varias conexiones semanales Frecuencias diarias necesarias usando la aeronave pequeña MAD PMI BCN LPA TFN IBZ BIO ACE VGO LCG AGP SCQ SVQ FUE TFS MAH OVD VLC ALC EAS SDR GRX ZAZ XRY Aeropuertos Origen MAD PMI BCN LPA TFN IBZ BIO ACE VGO LCG AGP SCQ SVQ FUE TFS MAH OVD VLC ALC EAS SDR GRX ZAZ XRY 0 2 4 6 8 10 12 Frecuencias diarias necesarias usando la aeronave grande MAD PMI BCN LPA TFN IBZ BIO ACE VGO LCG AGP SCQ SVQ FUE TFS MAH OVD VLC ALC EAS SDR GRX ZAZ XRY Aeropuertos Origen MAD PMI BCN LPA TFN IBZ BIO ACE VGO LCG AGP SCQ SVQ FUE TFS MAH OVD VLC ALC EAS SDR GRX ZAZ XRY Aeropuertos Destino 0 2 4 6 8 10 12 14 16 Frecuencias semanales necesarias usando la aeronave pequeña MAD PMI BCN LPA TFN IBZ BIO ACE VGO LCG AGP SCQ SVQ FUE TFS MAH OVD VLC ALC EAS SDR GRX ZAZ XRY Aeropuertos Origen MAD PMI BCN LPA TFN IBZ BIO ACE VGO LCG AGP SCQ SVQ FUE TFS MAH OVD VLC ALC EAS SDR GRX ZAZ XRY Aeropuertos Destino 0 1 2 3 4 5 6Frecuencias semanales necesarias usando la aeronave grande MAD PMI BCN LPA TFN IBZ BIO ACE VGO LCG AGP SCQ SVQ FUE TFS MAH OVD VLC ALC EAS SDR GRX ZAZ XRY Aeropuertos Origen MAD PMI BCN LPA TFN IBZ BIO ACE VGO LCG AGP SCQ SVQ FUE TFS MAH OVD VLC ALC EAS SDR GRX ZAZ XRY Aeropuertos Destino 0 1 2 3 4 5 6 Resultados del modelo 68 68 4.2.5 Análisis de los costes de operación Una vez relizado el analisis de las diferentes distribuciones de operaciones para cada aeropuerto. Se procede a relizar una estimación de los costes de operación de diferentes vuelos por pasajero.Si se utiliza como referencia los vuelos que con origen en el aeropuerto de Madrid-Barajas. Este aeropuerto al encontrarse en el centro de la península, la mayoria de rutas nacionales se encuandran dentro del rango donde el uso de la aeronave pequeña (Embraer 190) es más económico que la aeronave grande (Boeing 737-800). En la figura 41 se muestran los precios por pasajero de cada ruta con origen en Madrid-Barajas, ordenados de menor a mayor precio. Se puede observar como los costes totales de la aeronave pequeña (Rojo) son menores para la mayoria de los destinos peninsulares y de las Islas Balerares. Sin embargo para las rutas de mayor distancia, es decir aquellas con destino a los aeropuertos de las Islas Canarias, la aeronave grande (Azul) es la más barata de emplear por pasajero. En la imagen salta a la vista la diferencia de precios que supone un vuelo de larga duracíon. Si bien en la mayoría de vuelos cortos, los costes fijos del vuelo y los cotes dependientes de la distancia son similares en escala, para los vuelos de larga duración los costes variables suponen alrededor del 90% del costo total del vuelo. También se puede observar como los precios por asiento para ambas aeronaves son muy similares, habiendo una diferencia general de pocos céntimos o porcentualmente entre un 3 a 5% en el precio total. Una idea importante que vienen a expresar estos gráficos es que los costes en la aviación comercial suelen estar muy optimizados. Pequeños cambios en el valor de una variable como puede ser el coste de operación por asiento-milla disponible pueden suponer una diferencia radical en los resultados finales. Además, una aeronave puede ser la más rentable para relizar una ruta pero diferentes aspectos como la disponiblidad de dicha aeronave o que estas no sea la más optima para el movimiento de grandes demandas pueden hacer que aunque sea la más barata de operar no resulte en la más óptima de emplear. Aunque la figura 41 muestre que para todos los vuelos de corta distancia (Penisulares e Islas Baleaers) sea más rentable emplear la aeronave pequeña (Excepto en el caso del aeropuerto de Menorca), el modelo creado haciendo uso de las diferentes restricciones operativas estima que es más optimo el uso mayoritario de un Boeing 737-800 para cubrir rutas que unen Madrid con aeropuertos como Palma de Mallorca, Ibiza ,etc. Si se relizase la misma comparativa de precios para otros aeropuertos en el que si exista mayores distancias entre rutas peninusalres, como Barcelona-El Prat. En la figura 42 (Derecha) se observa como si existe un mayor número de rutas dentro de la península en las que si es mas rentable el empleo de la aeronave grande. Para los aeropuertos peninsulares más lejanos a Barcelona como Vigo o Sevilla, es más rentable el uso de la aeronave grande. Por otro lado para el aeropuerto Tenerife Norte todos los vuelos hacia la península tienen un precio por pasajero entre 48.47€ para el aeropuerto peninsular más cercano (Jerez de la Frontera) y 79.14€ para el más lejano (Barcelona-El Prat) 69 69 Análisis y optimización de la oferta de transporte aéreo en España Figura 41 Análisis de los precios por asiento de ambas aeronaves empleadas para todos vuelos a los aeropuertos de la red con origen en el aeropuerto de Madrid-Barajas Figura 42 Análisis de los precios por asiento de ambas aeronaves empleadas para todos vuelos a los aeropuertos de la red con origen en los aeropuertos de Barcelona-El Prat (Izquierda) y Tenerife Norte (Derecha) ZAZ VLC EAS SDR BIO GRX OVD ALC SVQ VGO XRY IBZ AGP LCG SCQ BCN PMI MAH ACE FUE TFN LPA TFS Aeropuertos destino 0 10 20 30 40 50 60 70 Precio de vuelo por asiento aeropuerto Madrid-Barajas Coste total Boeing 737 Coste total Embraer 190 Variables Boeing 737 Variables Embraer 190 Fijos Boeing 737 Fijos Embraer 190 MAH PMI ZAZ IBZ VLC EAS ALC BIO SDR MAD GRX OVD AGP SVQ XRY LCG VGO SCQ ACE FUE TFN LPA TFS Aeropuertos destino 0 10 20 30 40 50 60 70 80 90 Precio de vuelo por asiento aeropuerto origen: Barcelona-El Prat Coste total Boeing 737 Coste total Embraer 190 Variables Boeing 737 Variables Embraer 190 Fijos Boeing 737 Fijos Embraer 190 TFS LPA BIO FUE ACE XRY SVQ AGP GRX VGO SCQ LCG MAD ALC OVD VLC SDR ZAZ IBZ EAS PMI BCN MAH Aeropuertos destino 0 10 20 30 40 50 60 70 80 90 Precio de vuelo por asiento aeropuerto origen: Tenerife Norte Coste total Boeing 737 Coste total Embraer 190 Variables Boeing 737 Variables Embraer 190 Fijos Boeing 737 Fijos Embraer 190 Resultados del modelo 76 76 Aeropuerto ALICANTE-ELCHE MIGUEL HDEZ. 35 63 55 31 108 189 135 81 40 MAD PMI BCN LPA TFN IBZ BIO ACE VGO LCG AGP SCQ SVQ FUE TFS MAH OVD VLC ALC EAS SDR GRX ZAZ XRY Aeropuerto destino 0 20 40 60 80 100 120 140 160 180 200 nº de vuelos A321 Neo A319-200 Aeropuerto SAN SEBASTIÁN 67 39 MAD PMI BCN LPA TFN IBZ BIO ACE VGO LCG AGP SCQ SVQ FUE TFS MAH OVD VLC ALC EAS SDR GRX ZAZ XRY Aeropuerto destino 0 10 20 30 40 50 60 70 nº de vuelos A321 Neo A319-200 Aeropuerto SANTANDER-SEVE BALLESTEROS 24 12 18 76 75 MAD PMI BCN LPA TFN IBZ BIO ACE VGO LCG AGP SCQ SVQ FUE TFS MAH OVD VLC ALC EAS SDR GRX ZAZ XRY Aeropuerto destino 0 10 20 30 40 50 60 70 80 nº de vuelos A321 Neo A319-200 Aeropuerto F.G.L. GRANADA - JAÉN 42 85 40 MAD PMI BCN LPA TFN IBZ BIO ACE VGO LCG AGP SCQ SVQ FUE TFS MAH OVD VLC ALC EAS SDR GRX ZAZ XRY Aeropuerto destino 0 10 20 30 40 50 60 70 80 90 nº de vuelos A321 Neo A319-200 Aeropuerto ZARAGOZA 13 14 52 MAD PMI BCN LPA TFN IBZ BIO ACE VGO LCG AGP SCQ SVQ FUE TFS MAH OVD VLC ALC EAS SDR GRX ZAZ XRY Aeropuerto destino 0 10 20 30 40 50 60 nº de vuelos A321 Neo A319-200 Aeropuerto JEREZ DE LA FRONTERA/ LA PARRA 28 64 15 67 MAD PMI BCN LPA TFN IBZ BIO ACE VGO LCG AGP SCQ SVQ FUE TFS MAH OVD VLC ALC EAS SDR GRX ZAZ XRY Aeropuerto destino 0 10 20 30 40 50 60 70 nº de vuelos A321 Neo A319-200 Figura 49 Resultados para los aeropuertos Alicante-Elche, San Sebastián, Santander-Severiano Ballesteros, Granada-Jaén, Zaragoza y Jerez de la Frontera 77 77 Análisis y optimización de la oferta de transporte aéreo en España 4.3.1 Comparación de los resultaos obtenidos empleando un modelo con aeronaves de muy diferente capacidad y aeronaves con capacidades más similares Una vez mostrados los resultados generales que se obtendrán con el uso respectivo del par de aeronave A319/A321 y Embraer 190/B737 es interesante realizar una comparación global del uso de cada aeronave para cada uno de los dos modelos simulados. Como es lógico el empleo de dos aeronaves de mayor tamaño que las empleadas con anterioridad dará como resultado una disminución del número de operaciones totales necesarias para suplir la demanda dada. Pasando de 28080 operaciones totales para el caso del empleo del par de aeronave Boeing 737 y Embraer a 24064 para el nuevo escenario. En la figura 50 se muestra el porcentaje de uso de cada tipo de aeronave en cada modelo, empleando las aeronaves A319 y A321 arriba y el Embraer 190 y Boeing 737 abajo. Es interesante observar como para el modelo anterior, el uso del Embraer se da en mayor medida dentro de su modelo que el del A319. Esto se debe principalmente a que al ser la aeronave menos pesada de todas y existe una mayor diferencia de peso con la grande (B737) que en el caso del uso de A319 y A321 es la aeronave más atractiva para operar el conjunto de rutas de muy corta distancia o menores a unos 500 kilómetros. Por otro lado, para el caso de las aeronaves grandes de cada modelo, el empleo de la aeronave A321 respecto al A319 se da en un rango de distancias más amplio que para el caso del Boeing 737 en comparación del Embraer. El porcentaje de pasajeros transportado por cada aeronave en ambos modelos en muy similar. Figura 50 Porcentaje de uso de cada aeronave y porcentaje de pasajeros transportados por cada aeronave en los dos modelos simulados. A319 Y A321 (arriba), Embraer 190 y B737 (abajo) Por ello se concluye que en general y a grandes rasgos no existe una gran diferencia entre ambos modelos. A nivel local si se pueden observar diferencias considerables, si se emplean las aeronaves Embraer 190 y B737 ocurre con mayor frecuencia el hecho de que ciertas rutas tengan que emplearse ambas aeronaves. Esto se debe a la inclusión de la restricción de operaciones en rutas con mucha demanda, si se estima que la aeronave pequeña es la más rentable se intenta emplear mientras que la restricción de operaciones lo permita. En el caso de empleo de las aeronaves A319 y A321 este fenómeno se da con menor frecuencia. Con las capacidades de Porcentaje de Uso de Aviones Operaciones totales A319 Neo: 9807 Operaciones totales A321 Neo: 14257 41% 59% % A319 Neo % A321 Neo Cantidad de Pasajeros por Tipo de Avión Pasajeros totales transportados por A319 Neo: 1274910 Pasajeros totales transportados por A321 Neo: 2780115 31% 69% % A319 Neo % A321 Neo Porcentaje de Uso de Aviones Operaciones totales Embraer 190: 1.257500e+04 Operaciones totales Boeing 737-800: 1.550500e+04 45% 55% % Embaer 190 % Boeing 737-800 Cantidad de Pasajeros por Tipo de Avión Pasajeros totales transportados Embraer 190: 1.257500e+06 Pasajeros totales transportados Boeing 737-800: 2.790900e+06 31% 69% % Embraer 190 % Boeing 737-800 Resultados del modelo 78 78 estos dos modelos el conjunto de restricciones empleadas no produce efectos tan agresivos como en el caso del primer modelo, pues en muchas rutas de alta demanda sea cual sea la aeronave más barata de emplear puede cubrir la demanda dada en dicha ruta con un menor número de operaciones que las que se dieron realmente. 4.4 Otros resultados de interés 4.4.1 Rango de operación óptimo de cada aeronave empleada Aunque por limitación del modelo sea imposible incluir más de dos aeronaves, a partir de los datos empleados si se puede calcular que aeronave de las cuatro contempladas es más rentable de operar en cada rango de distancia. Para ellos se calcula el precio del transporte de cada pasajero en función de los costes asociados a cada aeronave. Puesto que los precios fijos de cada destino varían dentro del modelo se va a emplear el aeropuerto de Barcelona-El Prat como referencia. Los rangos dependerán de dichas tarifas, pero el ejemplo mostrado serviría como una estimación. Por lo general el modelo Embraer 190 sería el más económico de emplear en rutas menores a unos 500 kilómetros, el Boeing 737 en el rango comprendido entre los 500 km a 1100 km y para distancia mayores a dicha cifra lo sería el Airbus A321. El Airbus A319 no es la más rentable para ninguna distancia de vuelo. No obstante, tampoco es la más cara de emplear para ninguna ruta. Figura 51 Rango en el que cada aeronave es más barata de emplear por pasajero según los datos contemplados en el trabajo En todo caso estos valores son orientativos, además, cualquier pequeño cambio en cualquier variable que afecte al coste de operación de una aeronave puede generar cambios importantes en las estimaciones realizadas. Por poner un ejemplo si se cambiase el CASM del Embraer 190 de 0.065 €/milla a 0.067 €/milla, o lo que es lo mismo un aumento únicamente del 3% del precio por milla. El rango en el que esta resulta más rentable de operar dicha aeronave se reduciría en el entorno de 200 km para solo ser la más rentable en 79 79 Análisis y optimización de la oferta de transporte aéreo en España operaciones de menos de unos 300 km. 4.4.2 Viabilidad computacional del problema de optimización de asignación de aeronaves El tiempo de ejecución de los algoritmos utilizados para obtenerlos diferentes resultados es una cuestión de vital importancia para el correcto empleo del modelo. Este tipo de problema se caracteriza por involucrar un gran número de variables, y el uso de algoritmos no deterministas de tipo polinómico. Puede parecer que a medida que aumenta el número de variables, el tamaño del problema crece de manera exponencial, lo que resultaría en un tiempo de ejecución del programa mayor, lo cual no sería cierto del todo. Además de la consideración del número de variables a optimizar, se debe tener en cuenta el conjunto de datos que modelan dichas variables. El problema de forma general busca minimizar el coste de operación, por lo que dicho coste puede tener un papel principal en la factibilidad de la obtención de una solución. En este contexto, al emplear el algoritmo creado sin restricciones de tipo operativo o de distancias, a medida que los costes de empleo de una aeronave u otra por pasajero se aproximan, mayor coste computacional requerirá la solución y por tanto más tiempo de ejecución requerirá. Por ejemplo, si el costo fijo de vuelo por pasajero para ambas aeronaves y el coste por asiento-milla disponible son iguales. A falta de restricciones el algoritmo determinará que existen numerosas soluciones viables. Se está generando un espacio de búsqueda demasiado grande o complejo para el algoritmo de optimización, lo que lleva a tiempos de ejecución prolongados o infinitos. En el problema planteado se ha visto que los costos de empleo de diferentes aeronaves por asiento son muy similares. Tanto referidos a fijos y variables por la distancia. Es por ello por lo que la elección de estas variables no puede ser aleatoria si se quiere obtener un resultado de forma garantizada. Si se consideran establecidos los costos fijos y solo se varía el CASM de cada aeronave dentro del modelo no se pueden tener en cuenta cualquier rango de valores. A medida que la diferencia relativa entre ambos es menor, mayor complejidad computacional tendrá el problema a optimizar. Si las pendientes de las rectas que estiman el coste de operación de cada aeronave son muy similares el algoritmo encargado de la optimización de la función (‘intlingprog’) tardará más en resolver el problema. Si se realizan sendas simulaciones variando de forma relativa el CASM de ambas aeronaves, obtenemos los siguientes tiempos de simulación para cada rango de valores. Para conocer dichas duraciones se emplea el comando tic y toc de MATLAB®. En la tabla 19 se muestran los tiempos de ejecución empleando diferentes rangos de diferencias relativas entre los valores de CASM de ambas aeronaves. A medida que dicha diferencia disminuye por lo general aumenta el tiempo de simulación. Esto puede generar grandes problemas si se desea realizar muchas simulaciones para comprobar diferentes resultados. Es por ello por lo que no se puede recomendar el uso de Costes por Asiento-Milla Disponible con una diferencia relativa por debajo del 5% con los costes fijos contemplados en el problema. En todo caso los valores que se van a emplear en el modelo tienen una diferencia porcentual en el entorno de 7-10% lo cual permite obtener simulaciones con una duración menor a 2 segundos en término medio. Valores aceptables en todo caso. En todo caso, estas simulaciones pueden mostrar la fragilidad del modelo a la hora de obtener resultados y cómo cualquier conjunto de datos a optimizar puede no tener una solución viable. La realización de un estudio previo de todos estos y como afectarán al modelo es siempre una recomendación, además de una correcta definición de restricciones dentro del problema que por lo general reducen los tiempos de ejecución. Resultados del modelo 80 80 Tabla 19 Tiempos de ejecución del modelo en función de la variación relativa de CASM de las dos aeronaves empleadas sin restricciones Valores del CASM Diferencia porcentual Tiempo de simulación (Segundos) [0.07-0.06] 14% 1.134 [0.068-0.06] 11% 1.482 [0.066-0.06] 9% 1.835 [0.064-0.06] 6% 27.003 [0.062-0.06] 3% 2013.07 81 81 Análisis y optimización de la oferta de transporte aéreo en España 5 CONCLUSIONES Y LÍNEAS DE TRABAJO FUTURAS 5.1 Conclusiones apartado 2 De los resultados obtenidos de los diferentes análisis estadísticos realizados en el apartado 2 se puede concluir que el tráfico en los aeropuertos españoles es irregular. Por un lado, un sector turístico líder a nivel mundial en España genera un conjunto de aeropuertos con un volumen de tráfico importante. Este fenómeno coexiste con el hecho de que el número de aeropuertos en España es elevado, y suele darse el caso de aeropuertos de poco o nulo tráfico cerca de aeropuertos mucho mayores que dificultan justificar su existencia. Si comparamos la situación de los aeropuertos españoles con la de otros países de naturaleza similar, podemos observar que España podría ser un caso fuera de la norma. En países de tamaño y población comparables, es común encontrar una estructura aeroportuaria donde solo unos pocos aeropuertos concentran gran parte del tráfico aéreo, mientras que el resto de los aeropuertos registran niveles de tráfico relativamente bajos. Este patrón es notable en países como Alemania o Francia, donde la disparidad en la distribución de pasajeros entre los aeropuertos es mucho más pronunciada que en España. La particularidad del caso español radica en que, además del significativo volumen de tráfico de los tres principales aeropuertos nacionales, los siguientes 8 a 10 aeropuertos en términos de tráfico mantienen una cercanía considerable con respecto a estos aeropuertos principales, una situación menos común en otros países donde la distancia entre el "top 3" de aeropuertos y los siguientes es más amplia. Los aeropuertos españoles exhiben una distribución del tráfico aéreo nacional más equilibrada en comparación con los del Reino Unido, y esto se debe en gran medida a varios factores distintivos del contexto español. En primer lugar, las distancias entre ciudades en España son, en promedio, mayores que las del Reino Unido, lo que fomenta el uso del transporte aéreo para conectar regiones distantes dentro del país. Además, la presencia de dos archipiélagos en España, los cuales dependen en gran medida del uso del avión para la comunicación con la península, contribuye significativamente a esta distribución más uniforme del tráfico aéreo nacional. En contraste, el panorama aeroportuario español se diferencia notablemente del estadounidense, donde las vastas distancias dentro del país y una población considerablemente mayor generan una gran demanda de vuelos nacionales. Esto se traduce en una distribución más uniforme de pasajeros entre los aeropuertos y las rutas nacionales, reflejando las necesidades de conectividad dentro de un territorio extenso y densamente poblado como el estadounidense. 5.2 Conclusiones apartados 3 y 4 Sin duda el fenómeno más interesante a destacar de todo el modelo es lo increíblemente variable que pueden ser los resultados en función de los datos empleados de precios de diferentes variables. Cambios sutiles en solo un par de toneladas en el peso máximo de despegue las aeronaves o décimas, o incluso centésimas, de céntimo en el coste operativo por milla pueden resultar en el abandono absoluto del modelo del uso de una aeronave u otra. Los datos empleados, por ejemplo, del CASM de cada aeronave son solo una estimación realizada por el autor. Emplear un valor u otro dentro del rango de la décima de céntimo resulta más una elección subjetiva que objetiva. A lo largo del desarrollo del modelo se ha dado un hecho a destacar, a medida que todos los datos empleados intentaban acercarse más a la realidad, más ajustado son los resultados obtenidos. Lo que se quiere decir con esto es que, por ejemplo, al contemplarse costes reales de tarifas aeroportuarias o más ajustado está el coste por asiento-milla disponible al coste real de variables como el precio del combustible y el consumo de cada Conclusiones y líneas de trabajo futuras 82 82 aeronave a valores de diferentes estudios; menor es el rango en el que se pueden variar estos valores para que se arrojen resultados de interés. No aporta valor a los resultados mostrados el hecho que se prefiera de forma casi exclusiva una aeronave u otra, pero un pequeño desajuste en los valores empleados acabaría en resultados como el fenómeno nombrado. Puede considerarse el modelo creado más como una herramienta de simulación de resultados o de ensayo y error, que como un modelo que arroja inequívocamente resultados completamente reales. El modelo admite la inclusión de cualquier tipo de restricción que el usuario desee contemplar, mientras que tengan sentido físico y la función de optimización de enteros encuentre un resultado. También se han ensayado simulaciones cambiando la demanda y operaciones por otros escenarios temporales como podría ser un mes de baja demanda o el conjunto de pasajeros y operaciones de un año completo. El cambio de escenario no afecta mucho a la naturaleza de los resultados. A iguales datos de costes para cada aeronave la más optima a emplear en cada ruta sería la misma y solo el cambio en el número de operaciones y por tanto en las restricciones ocasionaría cambios de interés en los resultados. En pro de una memoria más corta no se han incluido estos resultados. Por otro lado, la imposibilidad de incluir un mayor número de aeronaves dentro del modelo también ha supuesto un hándicap en la realización de este. La inclusión del Embraer 190 dentro de las primeras simulaciones se ha dado principalmente al hecho de que, si se empleaba una aeronave de algo menos de 200 pasajeros de capacidad, emplear otra aeronave de algo más de 100 pasajeros no reflejaría de forma correcta el hecho de que en España se realizan muchas operaciones empleado aeronaves mucho más pequeñas (véase el uso de ATR-72 en las Canarias). Intentar modela el conjunto de operaciones dentro del país haciendo uso de solo dos modelos es tarea imposible si lo que se quiere es obtener resultados que sean más parecidos a los que se dan en la realidad. Dicho todo esto se pasa a enunciar diferentes maneras en las que se podría proseguir con el trabajo. En el apartado 3.4 se hace referencia a la mayor parte de fenómenos que no se contemplan dentro del modelo. La inclusión de cualquiera de estos podría arrojar resultados de mayor interés que los obtenidos. Sin dudas, la posibilidad de incluir un mayor número de aeronaves sería un buen comienzo. La inclusión de una aeronave como el ATR-72 en detrimento del Embraer 190 y utilizar cuatro aeronaves con capacidades que oscilan los 80, 120 160 y 200 pasajeros sería incluso mejor. Por otro lado, el hecho de no disponer de datos reales de pasajeros con escala limita la posibilidad de obtener resultados más analizables. Saber esto permitiría realizar análisis de cuáles podrían ser los mejores aeropuertos dentro de la red a considerar como Hubs. Si consideramos fij como los pasajeros que vuelan del aeropuerto i al aeropuerto j, que es lo que se ha empleado en el modelo. Esta variable estaría comprendida por todos los pasajeros que vuelan directo de i a j, además de todos los pasajeros que vuelan de un aeropuerto h diferente a i y j con escala en i para llegar a j, los pasajeros que vuelan de i a h con escala en j y los que lo hacen de i a j con escala en h como se muestra en la ecuación 5-1: [23] 𝑥!" +* # 𝑥(!" 2 (&' (3!," +* # 𝑥!"( 2 (&' (3!," + # 𝑥!(" 2 (&' (3!," =𝑓!"* ( 5−1 ) * La determinación de dichas variables sin datos estadísticos certeros sería un problema estadístico complejo. Pero si pudiesen ser obtenidos, la restricción de capacidad de vuelos (ecuación 3-5), quedaría de la siguiente manera. Es decir, se deben eliminar del conjunto de pasajeros que vuelan de i a j con escala en h. 83 83 Análisis y optimización de la oferta de transporte aéreo en España 𝑥!" +* # 𝑥(!" 2 (&' (3!," +* # 𝑥!"( 2 (&' (3!," ≤* # 𝑏, ) ,&' *𝑦!" ,*↔*𝑓!" − # 𝑥!(" 2 (&' (3!," = # 𝑏, ) ,&' *𝑦!" ,* ( 5−2 ) Además, se podría realizar una estimación de cuáles serían los mayores hubs aeroportuarios dentro del mercado nacional. A partir de la ecuación 5-3, que determina el número de pasajeros de tránsito en un aeropuerto h en comparación con origen o destino en él. 𝑓(= ∑ ∑ 𝑥!(" 2 "&' "3(,! 2 !&' !3(," ∑ 𝑥"( 2 "&' "3( +* ∑ 𝑥!( 2 !&' !3( + ∑ ∑ ( 𝑥!"(+𝑥(!" ) 2 "&' "3(,! 2 !&' !3(," * ( 5−3 ) Otra manera interesante de continuar el estudio seria realizar una estimación de la flota de aeronaves mínima necesaria para cubrir el conjunto total de operaciones que calcula el modelo que se deben realizar para cubrir la demanda dada. Para ello se podrían considerar diferentes naturalezas de rutas, por ejemplo, para rutas de muy alta demanda estimar cuantas operaciones diarias podría realizar una aeronave y proporcionarlo al número de vuelos que serían necesarios. Para rutas de menor demanda, considerar en lugar de la ruta que une solo un par de aeropuertos, el conjunto de 3 o más rutas que unen diferentes aeropuertos y estimar cuantos vuelos del conjunto de rutas establecido son factibles a realizar por día y de la misma manera determinar el número de aeronaves que se requeriría para satisfacer el conjunto de rutas estipulados. Referencias 84 84 REFERENCIAS [1] Airships. DELAG: “The World’s First Airline”. [En línea] Disponible en: \https://www.airships.net/delagpassenger-zeppelins/ [Consultado el 22 de marzo de 2024] [2] Mulder, Rob. Europeanairlines. “Deutsche Luft-Reederei Gesellschaft mit beschränkter Haftung – DLR”. [En línea] Disponible en: \https://www.europeanairlines.no/wp-content/uploads/2010/06/A_DLR_300109_1- 2.pdf [Fecha de consulta: 22 marzo, 2024]. [3] Allan, John. “Britannica Money. Pan American World Airways, Inc.” [En línea] Disponible en: \https://www.britannica.com/money/Pan-American-World-Airways-Inc [Consultado el 22 de marzo de 2024] [4] Carryer, Simon. “El transporte aéreo en los años 30”. [En línea] Disponible en: \https://shadowlands.es/transporte-aereo-anos-30/ [Fecha de consulta: 23 marzo, 2024]. [5] Docsandroots. “Barcelona, febrero de 1910: el primer vuelo de aeroplano en España”. [En línea]. Disponible en: \https://www.docsandroots.com/2021/02/barcelona-1910-el-primer-vuelo-de-aeroplano- realizado-en-espana/ [Consultado el 22 de marzo de 2024] [6] Museo de Getafe. “Se solicita la Dehesa de Santa Quintería.” [En línea]. Disponible en: \https://museo.getafe.es/omeka/items/show/1561 [Fecha de consulta: 23 marzo, 2024]. [7] Archivo ABC. "Así revolucionó España la historia de los bombardeos aéreos durante la Guerra de Marruecos." [En línea]. Disponible en: \https://www.abc.es/archivo/abci-revoluciono-espana-historia- bombardeos-aereos-durante-guerra-marruecos-202207150352_noticia.html [Fecha de consulta: 25 marzo, 2024]. [8] Ivoox. "Motor y al Aire: Los inicios de la Aviación en España." [En línea]. Disponible en: \https://www.ivoox.com/podcast-motor-al-aire-historia-aviacion_sq_f1117313_1.html [Fecha de consulta: 20 de marzo, 2024]. [9] Ministerio de Transportes y Movilidad Sostenible. "100 años de transporte aéreo en España." [En línea] Disponible en: \https://www.transportes.gob.es/aviacion-civil/centenario-transporte-aereo- espania/historia/historia-de-los-100-anios-del-transporte-aereo [Fecha de consulta: 3 mayo, 2024]. [10] Wikipedia. "Plus Ultra (hidroavión)." [En línea]. Disponible en: \https://es.wikipedia.org/wiki/Plus_Ultra_(hidroavi%C3%B3n) [Consultado el 20 de marzo de 2024]. [11] Aena Estadísticas. [En línea]. Disponible en: \https://www.aena.es/es/estadisticas/inicio.html. [Consultado el 1 de marzo de 2024]. [12] Wikipedia. "Coeficiente de Gini." [En línea] Disponible en: \https://es.wikipedia.org/wiki/Coeficiente_de_Gini [Consultado el 30 de marzo de 2024]. [13] Del Castillo Granados, José María; Romero Pérez, Luis Miguel. "Ingeniería del transporte - 4º grado en Ingeniería Civil". Presentación, Universidad de Sevilla, Sevilla, España, 24 de octubre de 2023. [14] Wikipedia. "Anexo: Aeropuertos de España." [En línea] Disponible en: \https://es.wikipedia.org/wiki/Anexo:Aeropuertos_de_Espa%C3%B1a [Consultado el 30 de marzo de 2024]. [15] Comisión Nacional de los Mercados y la Competencia. "El sector aeroportuario en España: Situación actual y recomendaciones de liberalización." [En línea] 2014. Disponible en: \https://www.cnmc.es/sites/default/files/1185359_9.pdf [Consultado el 27 de marzo de 2024]. [16] Danziger, Zachary. "Theil-Sen Robust Linear Regression." [En línea]. MATLAB Central File Exchange. Disponible en: \https://www.mathworks.com/matlabcentral/fileexchange/48294-theil-sen-robust-linear- regression.[Fecha de consulta: 3 mayo, 2024]. [17] Eurostat. "Air passenger transport by reporting country." [En línea] Disponible en: \https://ec.europa.eu/eurostat/web/main/data/database. [Consultado: 10 de marzo de 2024]. 85 85 Análisis y optimización de la oferta de transporte aéreo en España [18] UK Civil Aviation Authority. "Annual airport data 2023." [En línea]. Disponible en: https://www.caa.co.uk/data-and-analysis/uk-aviation-market/airports/uk-airport-data/uk-airport-data- 2023/annual-2023/ [Consultado el 10 de marzo de 2024] [19] Union des Aeroports Français. "Résultat d'activité des aéroports français 2023." [En línea] Disponible en: https://www.aeroport.fr/public/page/uaf-134 [Consultado el 10 de marzo de 2024]. [20] ADV Flughafenverband. "Verkehrsstatistiken 2023." [En línea]. Disponible en: https://dashboard.adv.aero/statistiken. [Consultado el 10 de marzo de 2024]. [21] Bureau of Transportation Statistics. "TranStats." [En línea]. Disponible en: https://www.transtats.bts.gov/databases.asp?Z1qr_VQ=E&Z1qr_Qr5p=N8vn6v10&f7owrp6_VQF=D .[Consultado el 10 de marzo de 2024]. [22] Wikipedia. "List of busiest airports by passenger traffic." [En línea]. Disponible en: https://en.wikipedia.org/wiki/List_of_busiest_airports_by_passenger_traffic. [Consultado el 20 de abril de 2024]. [23] Del Castillo Granados, José María. "Aplicación de la programación lineal al transporte aéreo: problemas de explotación y gestión". Presentación sobre el tema 5 de "Explotación del Transporte Aéreo", curso 4º de Ingeniería Aeronáutica, Universidad de Sevilla [24] Fastercapital.com. "Coste por asiento-milla disponible: CASM: Analizando el impacto en RASM." [En línea]. Disponible en: https://fastercapital.com/es/contenido/Coste-por-asiento-milla-disponible--CASM-- Analizando-el-impacto-en-RASM.html#:~:text=que%20afectan%20CASM- ,El%20costo%20por%20asiento%2Dmilla%20disponible%20(CASM)%20es%20una,la%20rentabilidad%20 de%20una%20aerol%C3%A .[Consultado el 20 de abril de 2024]. [25] FAA. "Aircraft Operating Costs." [En línea]. Disponible en: https://www.faa.gov/sites/faa.gov/files/regulations_policies/policy_guidance/benefit_cost/econ-value-section- 4-op-costs.pdf [Consultado el 18 de abril de 2024] [26] Enaire. "Servicios de tránsito aéreo (ATS)." [En línea]. Disponible en: https://www.enaire.es/servicios/gestion_del_transito_aereo/servicios_de_transito_aereo . [Consultado el 10 de abril de 2024] [27] MathWorks. "intlinprog." [En línea]. Disponible en: https://es.mathworks.com/help/optim/ug/intlinprog.html [Consultado el 1 de marzo de 2024] [28] ESTATISTA. "Número de pasajeros de las principales compañías aéreas en España en 2022 y 2023." [En línea]. Disponible en: https://es.statista.com/estadisticas/568641/aerolineas-con-mayor-numero-de-pasajeros- en-espana/ [Consultado el 3 de marzo de 2024]. [29] Wikipedia. "Bombardier CRJ700/900/1000." [En línea]. Disponible en: https://es.wikipedia.org/wiki/Bombardier_CRJ700/900/1000 [Consultado el 17 de abril de 2024] [30] AIRFLEETS. "Búsqueda de flota por avión o compañía aérea." [En línea]. Disponible en: https://www.airfleets.es/home/ [Consultado el 3 de marzo de 2024] [31]Wikipedia. "Boeing 737." [En línea]. Disponible en: https://es.wikipedia.org/wiki/Boeing_737 [Consultado el 10 de marzo de 2014] [32]Embraer E-Jets. [Online]. Disponible en: https://es.wikipedia.org/wiki/Embraer_E-Jets#E-190/195. [Consultado el 10 de marzo de 2024] [33] Wikipedia. "Ortodrómica." [En línea]. Disponible en: https://es.wikipedia.org/wiki/Ortodr%C3%B3mica [Consultado el 11 de abril de 2024]. [34] AENA. "Guía de tarifas marzo 2024." [En línea]. Disponible en: https://www.aena.es/es/oportunidadesde-negocio/precios-y-facturacion.html [Consultado el 15 de marzo de 2024] [35] Doganis, Rigas. "The Airline Business." Google Books. [En línea].Disponible en: https://www.google.es/books/edition/The_Airline_Business/BlFJUrOIto0C?hl=es&gbpv=1&dq=operating%2 0cost%20of%20boeing%20737&pg=PA154&printsec=frontcover [Consultado el 17 de marzo de 2024]. Anexos 92 92 ANEXO B En este anexo se describe brevemente el funcionamiento y estructuras de los códigos empleados para la creación del modelo. A modo de resumen, el modelo creado consta de un código principal que se encargará de la minimización de la función objetivo, además de permitir seleccionar diferentes variables de interés para que puedan ser modificadas con facilidad. Este código hará uso de diferentes funciones que se encargarán de tareas más específicas, como la generación de restricciones, lectura de datos necesarios, creación de la función objetivo, etc. A) Lectura de datos El conjunto de datos necesarios para la correcta ejecución del modelo es amplio. Por un lado, son necesaria la matriz de demanda de pasajeros (F), la matriz de operaciones entre aeropuertos (V), y de distancia entre aeropuertos (D). Estas matrices estarán formadas por datos reales de pasajeros y operaciones entre aeropuertos de datos estadísticos reales obtenidos de AENA estadísticas. Se generan sendas matrices donde cada columna representa el aeropuerto de origen y cada fila el aeropuerto de destino. Para la correcta introducción de dichos datos se ha seguido el siguiente orden: 1. Obtención de datos de la página de estadísticas y creación de la matriz en EXCEL®. 2. Exportar las matrices generadas de EXCEL® a MATLAB®en formato ‘table’ y guardadas todas dentro de un archivo ‘.mat’ 3. Conversión de dichas matrices en formato ‘table’ a ‘array’ para facilitar la lectura de datos del programa haciendo uso de la función de MATLAB® table2array. Por otro lado, para la introducción de costes y tarifas aeroportuarias, en el apartado de datos del problema se hace alusión todas las tarifas contempladas. Las tarifas de cada aeropuerto son diferentes. Es por ello por lo que a partir de las tarifas contempladas se crea una matriz de costes fijos de dimensión 𝑛×𝑛 (Igual que el resto de las matrices de datos) en la que cada valor de la matriz engloba la suma de todos los costes contemplados para dicho aeropuerto. Puesto que los costes fijos dependen del aeropuerto de destino de la operación la matriz creada se traspone para que en cada coordenada de la matriz aparezcan los costos aeroportuarios del aeropuerto de destino. B) Generación de restricciones La función restricciones.m tiene la función de crear las matrices y vectores que modelan las restricciones de desigualdad e igualdad, es decir, A, b, Aeq y beq. Cada restricción debe plantearse de manera individual, es decir creando sus respectivas matrices y vectores que la modelan. Una vez se tengan todas ellas generadas, se unen en una sola. Cada restricción de forma individual se compone de un vector fila de 𝑛×𝑛×𝑘* que contempla todas las variables que aparecen en el problema. Si las variables no son afectadas por dicha restricción, la posición en la que colocaría esta aparecerá un cero. En este apartado se procede a explicar cómo se han implementado algunas restricciones formuladas en el apartado de restricciones (apartado 3.1.2). Siguiendo el mismo orden en el que se enumeraron. • Restricción de capacidad de los vuelos: Para este caso tendremos una matriz A de dimensiones ( 𝑛×𝑛 , * 𝑛×𝑛×𝑘 ) y un vector b de longitud * 𝑛×𝑛 . 93 93 Análisis y optimización de la oferta de transporte aéreo en España Eso se puede deducir atendiendo a como se formula la restricción, se deberá generar una para cada par origen destino de aeropuertos. Es decir, un total de 𝑛×𝑛 restricciones totales. Cada fila de la matriz A que hace referencia a una conexión entre aeropuertos y estará formada por ceros en todas aquellas posiciones que no hagan referencia a dicha conexión. Supongamos el ejemplo de la conexión entre el aeropuerto i de la matriz con el aeropuerto destino j. A partir de la matriz de demanda F, el valor correspondiente a Fij será la demanda existente de pasajeros entre el aeropuerto i y j del modelo. En la figura 52 se hace referencia a como se estructuraría cada fila de la matriz A. El vector columna b sería la matriz de demanda F reordenada como vector empleando el comando de MATLAB® reshape. La siguiente figura hace referencia a cómo se crearía cada fila de la matriz de restricciones A y el vector b en función de la matriz de demanda F. Figura 52 : Creación de la matriz (A) y vector de desigualdad (b) que genera la restricción de capacidad. Puesto que las restricciones de desigualdad deben introducirse de la forma 𝐴*𝑥≤𝑏 , esta restricción deberá implementarse negando ambas expresiones de la forma −𝐴*𝑥≤−𝑏 . • Restricción de máximo número de operaciones por aeropuerto Para la implementar la restricción que modela el hecho de que para cada aeropuerto se puedan realizar como máximo la suma de operaciones de dicho aeropuerto según la matriz de operaciones V. Esta restricción estará compuesta de tantas restricciones como aeropuertos tenga la red (n aeropuertos). Es decir, la matriz de restricciones A tendrá dimensiones ( 𝑛×𝑛 , * 𝑛×𝑛×𝑘 ) y el vector columna b longitud n. La figura 53 hace alusión a como se generaría cada fila de la matriz de restricciones A. Donde los unos de cada fila representan todas las conexiones del aeropuerto i al resto de la red. Se debe repetir el proceso k veces Anexos 94 94 en función del número de aeronaves que se emplee. Figura 53: Creación de la matriz de desigualdad (A) que modela la restricción de operaciones débil. En cuanto al vector b, a partir de la matriz de operaciones V tendrá la siguiente forma: 𝑏= ⎣ ⎢ ⎢ ⎢ ⎢ ⎢ ⎢ ⎡ # 𝑉'" 2 "&'.... # 𝑉2" 2 "&' ⎦ ⎥ ⎥ ⎥ ⎥ ⎥ ⎥ ⎤ * ( 𝐵−1 ) • Restricción del uso de la aeronave de mayor tamaño para la realización de vuelos de muy corta duración. Esta restricción que limita el uso de cualquier tipo de aeronave que no sea la más pequeña de las contempladas dentro del modelo. Supóngase que las aeronaves se ordenan de menor a mayor tamaño, esta restricción de igualdad de la forma 𝐴*𝑥=𝑏 , donde la matriz A incluye todas las aeronaves que no pueden emplearse en dicha ruta y el vector b es un vector columna cuyas componentes son todas nulas. La matriz de igualdad A tendrá dimensiones ( 𝑛×𝑛 , * 𝑛×𝑛×𝑘 ) y el vector columna b longitud 𝑛×𝑛 . En la figura 54 se hace referencia a como se modela cada fila de la matriz A en el caso de que si en la matriz de distancia entre aeropuertos D una ruta sea más corta que la distancia umbral a partir de la cual aparece la restricción. En el caso de que una distancia entre aeropuertos sea mayor que dicha distancia umbral la fila de dicha restricción será una fila completa de ceros. 95 95 Análisis y optimización de la oferta de transporte aéreo en España Figura 54 Creación de la matriz y vector de igualdad que genera la restricción de operaciones de todas las aeronaves que no sean la más pequeña en vuelos de corta duración. Todas las matrices de restricción son matrices de grandes dimensiones, las cuales la mayor parte de sus valores serán ceros, para obtener una mayor velocidad de ejecución del modelo todas las matrices descritas anteriormente se cómo matrices y vectores dispersos. A) Resolución del problema Una vez generados correctamente todos los datos y restricciones contempladas se genera la función objetivo como se contempla en la ecuación 3-1 y se minimiza haciendo uso de la función de MATLAB® intlingprog que devuelve tanto el coste total de operar la red completa como los valores del número de operaciones a realizar con cada tipo de aeronave para cada ruta del modelo. Dichos valores son devueltos por la función en forma de vector de longitud, * 𝑛×𝑛×𝑘 , haciendo uso de la función reshape, se convierte dicho vector en sendas matrices (k matrices) donde cada matriz muestra el número de operaciones con la aeronave de tipo l para cada par origen destino de la misma manera que se generaron las matrices de demanda u operaciones entre aeropuertos.