Research on Flight Delay Based on Fuzzy Evaluation Algorithm

Show more

1. Introduction

Since the reform and opening-up, China’s civil aviation develops rapidly and has made remarkable achievements. Fixed-point flight traffic has been ranked second place in the International Civil Aviation Organization. China Civil Aviation’s average annual growth is 19.5% and it develops much quicker than any other transportation ways [1] . With the rapid developments, flight delays which influence civil aviation are also constantly exacerbated and upgraded. In recent years, about 25% of the flights are not normal according to China’s civil aviation, for which complaints caused by flight delay accounted for more than 45% of the total civil aviation complaints. Currently, the most prominent problem in China’s air transport industry is the growing flight delays caused by uneven distribution of aviation resources and airspace congestion. For all the reasons, study on flight delays is very important.

Research on reducing flight delays began much earlier in abroad. For example, in the optimization of flight plans, as early as the middle of the last century, how to optimize the flight time has already attracted the attention of the optimization experts of the airline flight plan and they carried out some researches. Mathaisel and Etshmaier [2] summed up internationally renowned scholars for the flight planning research from the 60 s of last century to the mid-80 s. They believed that the airline should use the structural optimization of the flight plan assessment and establishment of flight plan these two means to produce the best flight scheduling which can satisfy the actual demand. And Barnhart C and Corn A [3] , the two scholars have made a summary of the research results on existing airline flight scheduling. Researches on time optimization and flight frequency of flight scheduling are also rich in the foreign countries. For example, Wei and Hanse [4] proposed to establish the target planning model to maximize the profitability of the entire route and minimize passenger delay rate. Then it can determine the actual flight frequency and simultaneously carry out market share analysis to optimize the flight running time. In addition, Lohatepanont M and Barnhart C [5] also analyzed the relationship between the airline supply and the passengers demand and put forward the establishment of a comprehensive model optimization flight plan from the aspects of model assignment and route selection, etc. What’s more, one of the effective solutions is to implement air traffic flow management.

Research on flow control began in 1970s overseas. The traffic management system it used can effectively alleviate the partial traffic jam. In America, the United States Federal Aviation Administration has adopted an enhanced traffic management system. The system played a huge role in the system integration, solutions, monitoring and controlling, expert assessment and easy expansion, etc. In Europe, the central traffic management system was established in 1989 to provide ATFM services for all EU member states. In terms of flow control, foreign scholars also made a lot of researches. For instance, Sherali HD, Bish EK and Zhu X. [6] discussed airline fleet assignment concepts, models, and algorithms. Peter B. Vranas, Dimitris J. Bertsimas, Amedeo R. Odoni [7] and some other scholars [8] [9] [10] [11] [12] studied on the real-time problem of terrestrial waiting strategy and the multi-airport limited ground waiting strategy. In addition, the management of airspace is also one of the hotspots in air traffic management, and especially the collaborative decision-making mechanism has attracted much attention from many experts. The collaborative decision-making mechanism changed the previous centralized control model, and combined with the traditional optimization algorithm to improve the mathematical model. In the meantime, it enabled decentralized management collaborative decision-making model to the process of air traffic management. Initially, it implemented the CDM program in the United States, and used information sharing of air traffic control, airlines and airports three subsystems to achieve the effective resources integration. It has reached the purpose of controlling traffic, and the program has reached a certain effect.

Compared with foreign countries, domestic research on the problem of flight delays is inadequate. Most researchers paid attention to the flight delay service recovery strategy. They have discussed abundantly on service recovery, governance model and legislation of the flight delay. Although service remediation after flight delays is a solution to flight delays, it’s better to develop a more optimized flight plan by harmonizing the relationship among the airport sector, regulators, air traffic control departments, airlines and other aspects, so that it can better control the flow and can effectively solve the problem of flight delays. Scholars did this research. They studied flight Planning and optimization and provided mathematical models. Pita JP, Adler N, and Antunes AP [13] , based on the concept of collaborative cooperation, used collaborative mechanism to improve the efficiency of air transport system by integrating air traffic control, airlines and airport resources.

But there is no systematic research inland, from the perspective of the civil aviation departments and government departments, focusing on the main causes of flight delays and its improvement measures. It can be said that the study of domestic flight delays is in the ascendant. In this article, we will analyze the reasons for the delay of several airports which have higher delay rates and give suggestions. This is the innovation point of the paper.

2. Models and Questions

2.1. Problems’ Proposing and Analyzing

Nowadays, flight delay caused some troubles to travelers. In the article we collect related data and build the corresponding mathematical model to solve the following problems:

Question 1: Establish a model to determine the extent of an airport delay, and list 10 airports which has the most severe flight delays in China (or some other country).

Question 2: What is the main reason for the flight delay in our country (or some other country)?

Question 3: What are the corresponding improvements?

For Question 1, this paper collects the airport flight delay data first. It uses the fuzzy comprehensive evaluation method to construct the flight delay level assessment model, and determines the weight of the six indexes of the airport delay degree. Then the membership is calculated and the degree of airport delay will be determined.

For Question 2, this paper analysis the causes of flight delay from the aspects of its impact and the frequency. In the aspect of the influence degree, it builds grey correlation model to evaluate. And its influence can be analyzed by comparing the flight delay trends and the number of delays. Using the grey correlation model to assess the degree of impacts, you can determine the main reason for flight delays from these two angles.

For Question 3, this paper starts from the main impact-flow. On one hand, it qualitatively analysis some factors which influence airport traffic, including airport infrastructure construction, modes of transportation and other factors. On the other hand, quantitatively, it establishes a regression analysis model and tests the model. And finally it adds dummy variables to study the impact of policy factors on airport traffic.

2.2. Model Assumption

1) Each airport is consistent with the criteria of flight delays.

2) The source of the data is real and effective.

3) Reasons for the flight delay have no influences with each other.

4) Assuming that under normal circumstances, all flight indexes meet the flight requirements.

2.3. Symbols Description

$\rho $ : Distinguishing coefficient; ${x}_{i}$ : Correlation data array; ${x}_{0}$ :Standard data array;

$k$ : Moment; $\xi $ : Correlation coefficient; $r$ : Correlation;

$y$ : Sequence initialization results; ${\beta}_{i}$ : Parameter point estimation;

${R}^{2}$ : The square of the correlation coefficient R of the regression equation;

$F$ : Regression equation statistics; ${s}^{2}$ and ${\sigma}^{2}$ : estimated value;

${Y}_{i}$ : Air transport passenger turnover; ${X}_{2i}$ : Railway passenger turnover;

${X}_{3i}$ : Investment of Civil Aviation Capital Construction and Technical Renovation;

$D$ : Dummy variable.

3. Model’s Established and Solved

3.1. Flight Delay Assessment Model’s Established and Solved

Flight delay assessment index system is an organic whole which is composed of a series of interrelated indexes that intuitively reflect real-time flight delays in airspace. From the flight delay time and sphere of influence these two angles, this paper respectively selects time-related indexes and scope-related indexes to establish a real-time assessment index system for flight delays. According to spatial and temporal distribution feature, the article chooses the time-related indexes and scope-related indexes of the flight delays, and builds the real-time evaluation system for flight delays. It is shown in Table 1.

We determine the weight of the index by AHP (analytic hierarchy process)

Table 1. Flight delay real-time evaluation system.

1) Establish index set. According to the established index system, all indexes are shown below:

First class index:

$X=\left\{{X}_{1},{X}_{2}\right\}=\left\{\text{time-relatedindex},\text{scope-relatedindex}\right\}$

Second class index:

${X}_{1}=\left\{{X}_{11},{X}_{12}\right\}=\left\{\text{flightdelayaveragetime},\text{flightdelayratechange}\right\}$

2) Solution of index weight. Compare its impact of the same level index with each other in the index system by using of analytic hierarchy process. And it uses 1 - 9 scale method. Meanings of each index are shown in Table 2. Forming decision matrix according to the Table 1, and marking it as $W=\left({w}_{ij}\right)$ and

${w}_{ij}>0$ . Feature of the comparative matrix is ${a}_{ij}>0$ , ${a}_{ii}=1$ , ${a}_{ji}=\frac{1}{{a}_{ij}}$ .

Solving the eigenvalue problem ${W}_{n}={n}_{n}$ and doing a consistency check for calculation, after normalization processing, weight of index $X=\left\{{X}_{1},{X}_{2}\right\}$ can be calculated and it is marked as $U=\left\{{u}_{1},{u}_{2}\right\}$ . Weight of index ${X}_{1}=\left\{{X}_{11},{X}_{12}\right\}$ is marked as ${U}_{1}=\left\{{u}_{11},{u}_{12}\right\}$ . Weight of index ${X}_{2}=\left\{{X}_{21},{X}_{22}{X}_{23},{X}_{24}\right\}$ is marked as

${U}_{2}=\left\{{u}_{21},{u}_{22},{u}_{23},{u}_{24}\right\}$

3) Determine each index weight. First-level and second-level judgment matrices are as follow:

$W=\left[\begin{array}{cc}1& 1\\ 1& 1\end{array}\right],\text{\hspace{0.17em}}\text{\hspace{0.17em}}W=\left[\begin{array}{cc}1& 2\\ 1/2& 1\end{array}\right],\text{\hspace{0.17em}}\text{\hspace{0.17em}}W=\left[\begin{array}{cccc}1& 2& 2& 3\\ 1/2& 1& 1& 2\\ 1/2& 1& 1& 2\\ 1/2& 1/2& 1/2& 1\end{array}\right]$

After consistency test, we can get

$U=\left(0.500,0.500\right),\text{\hspace{0.17em}}\text{\hspace{0.17em}}{U}_{1}=\left(0.667,0.333\right),\text{\hspace{0.17em}}\text{\hspace{0.17em}}{U}_{2}=\left(0.423,0.227,0.227,0.123\right)$

In the paper, the fuzzy comprehensive evaluation method is used to construct the flight delay level assessment model.

According to the features of flight delay, we divide flight delay degree into 5 levels. Value of each level is shown in Table 3.

Table 2. Scale and its definition.

Table 3. Flight delay real-time assessment indicators and evaluative standards.

The fuzzy level of each evaluation index is determined by the triangular membership function. The function is shown in Figure 1.

The membership function of flight average delay time is

$C\left({x}_{1}\right)=\{\begin{array}{l}0\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}{x}_{1}\le 10\\ \frac{{x}_{1}-10}{40-10}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}10<{x}_{1}<40\\ 1\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}{x}_{1}\ge 40\end{array}$ (3.1)

Membership function of flight delay rate change is

$C\left({x}_{2}\right)=\{\begin{array}{l}0\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}{x}_{2}\le -9\\ \frac{{x}_{2}-\left(-9\right)}{18-\left(-9\right)}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}-9<{x}_{2}<9\\ 1\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}{x}_{2}\ge 18\end{array}$ (3.2)

Membership function of flight delay rate, airport delay rate, air control and weather is:

$C\left({x}_{i}\right)=\{\begin{array}{l}0\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}{x}_{i}\le 20\\ \frac{{x}_{i}-20}{80-20}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}20<{x}_{i}<80\\ 1\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}{x}_{i}\ge 80\end{array}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\left(i=3,4,5,6\right)$ (3.3)

Substituting data of delays to Formulas (3.1) - (3.3), it will get C. Then we can get the membership of each delay level. Thus we can get fuzzy membership matrix ${R}_{1}$ and ${R}_{2}$ of time-related indexes and scope-related indexes, where

${R}_{1}=\left[\begin{array}{ccccc}{a}_{11}& {a}_{12}& {a}_{13}& {a}_{14}& {a}_{15}\\ {a}_{21}& {a}_{22}& {a}_{23}& {a}_{24}& {a}_{25}\end{array}\right],\text{\hspace{0.17em}}{R}_{2}=\left[\begin{array}{ccccc}{b}_{11}& {b}_{12}& {b}_{13}& {b}_{14}& {b}_{15}\\ {b}_{21}& {b}_{22}& {b}_{23}& {b}_{24}& {b}_{25}\\ {b}_{31}& {b}_{32}& {b}_{33}& {b}_{34}& {b}_{35}\\ {b}_{41}& {b}_{42}& {b}_{43}& {b}_{44}& {b}_{45}\end{array}\right]$ (3.4)

We use the operator $M\left(\u2022,\oplus \right)$ to calculate the results by combining the weight of the secondary index and the fuzzy membership matrix. Then we can get the evaluation result of time-related index

$A={U}_{1}\u2022{R}_{1}=\left({a}_{1},{a}_{2},{a}_{3},{a}_{4},{a}_{5}\right),\text{\hspace{0.17em}}\text{\hspace{0.17em}}B={U}_{2}\u2022{R}_{2}=\left({b}_{1},{b}_{2},{b}_{3},{b}_{4},{b}_{5}\right)$ (3.5)

Then we can get the evaluation result of flight delay levels

Figure 1. Triangular membership function and membership of delay level.

$D=U\u2022R=\left({u}_{1},{u}_{2}\right)\u2022\left[\begin{array}{ccccc}{a}_{1}& {a}_{2}& {a}_{3}& {a}_{4}& {a}_{5}\\ {b}_{1}& {b}_{2}& {b}_{3}& {b}_{4}& {b}_{5}\end{array}\right]=\left({d}_{1},{d}_{2},{d}_{3},{d}_{4},{d}_{5}\right)$ (3.6)

According to the principle of maximum membership

$V=\mathrm{max}\left({d}_{1},{d}_{2},{d}_{3},{d}_{4},{d}_{5}\right)$ (3.7)

We get the flight delay levels in real time, and that is delay level.

Next, we check the average delay time and the punctuality rate and other data in 2014 to analyze the actually data according to the theory above. (Data are from Flight stats and Civil aviation resource network) the summary data of the Shanghai Pudong airport is shown in Table 4.

Substituting data of delays to Formulas (3.1) - (3.3), we get

$C\left({x}_{1}\right)=0.8077,\text{\hspace{0.17em}}\text{\hspace{0.17em}}C\left({x}_{2}\right)=0.8685,\text{\hspace{0.17em}}\text{\hspace{0.17em}}C\left({x}_{3}\right)=0.2833$

$C\left({x}_{4}\right)=0.8667,\text{\hspace{0.17em}}\text{\hspace{0.17em}}C\left({x}_{5}\right)=0.29,\text{\hspace{0.17em}}\text{\hspace{0.17em}}C\left({x}_{6}\right)=0.1067$

Then we get the membership of each delay level. Thus we can get fuzzy membership matrix ${R}_{1}$ and ${R}_{2}$ of time-related indexes and scope-related indexes.

${R}_{1}=\left[\begin{array}{ccccc}0& 0& 0& 0.4615& 0.5385\\ 0& 0& 0& 0.1575& 0.8425\end{array}\right],$

${R}_{2}=\left[\begin{array}{ccccc}0.0835& 0.9615& 0& 0& 0\\ 0& 0& 0& 0.1615& 0.8335\\ 0.05& 0.95& 0& 0& 0\\ 0.9665& 0.0335& 0& 0& 0\end{array}\right]$

Substituting ${R}_{1},{R}_{2}$ to Formulas (3.5) and (3.6), we can get:

$A=0\text{\hspace{0.17em}}\text{\hspace{0.17em}}0\text{\hspace{0.17em}}\text{\hspace{0.17em}}0\text{\hspace{0.17em}}\text{\hspace{0.17em}}0.3603\text{\hspace{0.17em}}\text{\hspace{0.17em}}0.6397$

$B=0.1656\text{\hspace{0.17em}}\text{\hspace{0.17em}}0.6444\text{\hspace{0.17em}}\text{\hspace{0.17em}}0\text{}0.0367\text{\hspace{0.17em}}\text{\hspace{0.17em}}0.1892$

$D=0.0828\text{\hspace{0.17em}}\text{\hspace{0.17em}}0.32220\text{\hspace{0.17em}}\text{\hspace{0.17em}}0.1985\text{\hspace{0.17em}}\text{\hspace{0.17em}}0.4145$

According to the principle of maximum membership (3.7), we get the result that flight delay level of Shanghai Pudong Airport is level 5. It’s in serious delay state, and this state is in accord with the data from flight stats and also consistent with the actual situation of the month.

Table 4. Flight delay data of Shanghai Pudong airport in June.

Statistics of the domestic large-scale airport punctuality rate from civil aviation administration website in 2014 is displayed that 10 airports have the most serious flight delays: Nanjing Lukou International Airport, Changzhou Airport, Wuxi Airport, Hangzhou Xiaoshan International Airport, Shanghai Hongqiao International Airport, Shanghai Pudong International Airport, Ningbo Lishe Airport, Tianjin Binhai International Airport, Shenzhen Baoan International Airport, Guangzhou Baiyun Airport.

3.2. Flight Delay Reasons

There are a variety of reasons for the flight delay. But in general reasons can be broadly divided into five categories: airline; flow control; bad weather; other reasons, including passenger reasons. At present, the major media published statistics can only reflect the frequency of each occurrence, but cannot reflect the influence extent of various reasons. For this situation, we use the grey relational analysis method to analyze the relationship between the various reasons and flight delays, and distinguish the impacts of these reasons. Thus, grey relational analysis should be introduced firstly.

For a sequence $x=\left(x\left(1\right),x\left(2\right),\cdots ,x\left(n\right)\right)$ , the mapping $f\left(x\left(k\right)\right)=y\left(k\right)$ , $k=1,2,\cdots ,n$ is called a data transformation from sequence $x$ to sequence $y$ . $f$ is called initialization transformation if

$f\left(x\left(k\right)\right)=\frac{x\left(k\right)}{x\left(1\right)}=y\left(k\right),\text{\hspace{0.17em}}x\left(1\right)\ne 0$

The concrete steps of solving relevancy are: select reference numeral.

${y}_{0}=\left\{{y}_{0}\left(k\right)|k=1,2,\cdots ,n\right\}=\left({x}_{0}\left(1\right),{x}_{0}\left(2\right),\cdots ,{x}_{0}\left(n\right)\right)$

Here $k$ indicates the moment. Suppose it has $m$ correlation data arrays.

${y}_{i}=\left\{{y}_{i}\left(k\right)|k=1,2,\cdots ,n\right\}=\left({y}_{i}\left(1\right),{y}_{i}\left(2\right),\cdots ,{y}_{i}\left(n\right)\right),\text{\hspace{0.17em}}i=1,2,\cdots ,m$

Then

${\xi}_{i}\left(k\right)=\frac{\underset{s}{\mathrm{min}}\underset{t}{\mathrm{min}}\left|{y}_{0}\left(t\right)-{y}_{s}\left(t\right)\right|+\rho \underset{s}{\mathrm{max}}\underset{t}{\mathrm{max}}\left|{y}_{0}\left(t\right)-{y}_{s}\left(t\right)\right|}{\left|{y}_{0}\left(k\right)-{y}_{i}\left(k\right)\right|+\rho \underset{s}{\mathrm{max}}\underset{t}{\mathrm{max}}\left|{y}_{0}\left(t\right)-{y}_{s}\left(t\right)\right|}$ (3.8)

is called correlation coefficient between correlation data array ${x}_{i}$ and standard data array ${x}_{0}$ at the $k$ moment. Also $\rho \in \left[0,1\right]$ is distinguishing coefficient, and in Formula (3.8), $\underset{s}{\mathrm{min}}\underset{t}{\mathrm{min}}\left|{y}_{0}\left(t\right)-{y}_{s}\left(t\right)\right|$ and $\underset{s}{\mathrm{max}}\underset{t}{\mathrm{max}}\left|{y}_{0}\left(t\right)-{y}_{s}\left(t\right)\right|$ is two-level minimum difference and two-level maximum difference respectively.

In general, the resolution coefficient and the resolution are of the same trend. In this case, we take resolution coefficient as 0.5. The correlation coefficient is an index that describes the degree of correlation between correlation data array and standard data array at a certain time. Since there is an association number at each moment, the information becomes too scattered and is not easy to compare.

${r}_{i}=\frac{1}{n}{\displaystyle \underset{k=1}{\overset{n}{\sum}}{\xi}_{i}\left(k\right)}$ (3.9)

is called the degree of correlation between ${x}_{i}$ and ${x}_{0}$ . It is easy to see the degree of correlation is the concentration of each moment correlation coefficient to focus on an average, that is, to achieve a scattered centralized information processing. Using the concept of correlation, we can analyze the various factors on the flight problems. At present, China Civil Aviation Authority has made statistics on the frequency of abnormal flights caused by various reasons. It is shown in Table 5.

The Figure 2 below shows the trends in the number of occurrences of the four causes from 2009 to 2014.

Initialization of each sequence is showed in Table 6.

Using gray system software to solve it, then we get the degree of correlation ( $\rho =0.5$ ) of each sequence.

Figure 2. Frequency changing trend of flight delays caused by various reasons.

Table 5. Nationwide frequency of abnormal flights caused by various reasons in 2009-2014.

Table 6. Initialization results of each sequence.

$r=0.686125\text{\hspace{0.17em}}\text{\hspace{0.17em}}0.721406\text{\hspace{0.17em}}\text{\hspace{0.17em}}0.681394\text{\hspace{0.17em}}\text{\hspace{0.17em}}0.64157$

Comprehensive analysis of all the reasons is shown in Table 7.

From Table 7, we can conclude that improper control flow is the main impact, and followed by airline reasons and weather reasons, and finally for other reasons.

3.3. Flow Factors Analysis

In order to ensure the most efficient using of airspace, the method to adjusting traffic flow of a given airspace and a given route to an airport is flow controlling. From the view of airports conditions: airport infrastructure construction also has a great impact on airport traffic. In addition, flight flow is different times. With the rapid development of China’s aviation and the rapid increase in the number of aircraft, the load on these roads is also growing. Moreover, in order to ensure the needs of national defense, the country restricts the airspace strictly, and there are many restricted areas in the air. And the military is responsible for organizing the implementation of the national flight control, while the adjustable margin of civil aviation is very small. From this point of view, the airport is passive control of traffic.

There are three modes of transportation available at the moment: air, rail and highway. The road transport is suitable for short trips. Its price fluctuation has little effect on air transport. And the proportion of waterways is small, so there is not much analysis about it. The biggest competitor of air transport is rail transport, especially the high-speed train. The market situation is their price and the actual time consumption. Insiders generally believe that high-speed train may cause a certain traffic divergence to airline traffic in a short-term. But over time, this diversion phenomenon will gradually weaken, thus this short supply situation of airport traffic won’t change too much.

Due to the strong uncontrollability of weather conditions and uneasy quantization of the airport and the airline these indexes, this article builds a model to analyze the impact on airport traffic only from the airport infrastructure and other transportation modes these two aspects. The paper takes air transport traveler turnover ${Y}_{i}$ as the standard to measure the airport traffic, takes rail transport traveler turnover ${X}_{2i}$ as the standard to measure the railway transport capacity and takes civil aviation capital construction and technical renovation investment ${X}_{3i}$ as yard stick. Then we can get regression equation:

Table 7. Frequency and extent of the factors that cause flight delays.

${Y}_{i}={\beta}_{0}+{\beta}_{1}{X}_{2i}+{\beta}_{2}{X}_{3i}+{\epsilon}_{i}$ (3.10)

Based on the relevant published data in China Statistical Yearbook and 2010 Civil Aviation Industry Development Statistics Bulletin, using Matlab to calculate the parameters and draw the residual map 3-3, then we can get an initial regression model:

${\stackrel{^}{Y}}_{i}=-50.1667+0.0825{X}_{2i}-0.218{X}_{3i}$

It can be seen that the second point in Figure 3 is an outlier. Then deleting the point and using regression again, one obtains coefficients, the coefficient confidence interval, the statistic, and the residual graph of the improved regression model. Details can be seen in Table 8 and Figure 4.

At this moment, confidence interval does not include null point. When F-statistics increases, coefficient of determination will increase from 0.9875 to 1.0000, and then we can get a regression model:

${\stackrel{^}{Y}}_{i}=-44.5673+0.0790{X}_{2i}-0.1682{X}_{3i}$

Having been tested, the model is well fitted to the sample. Both the model and the parameters of the test are passed. So it can be concluded that the passengers increase in railway transport won’t lead to the passenger reduction in air transport. Thus it proves the insiders’ view point that high-speed train has limited influence. It should be noted that no one person chooses the train as a travel tool on the corresponding 0.0790 million people choose the aircraft as a travel tool. It only indicates that an overall increasing in passenger traffic. When the investment in civil aviation infrastructure is reduced by 100 million RMB, the average passenger traffic volume is reduced by 0.168 billion tons/km. Therefore, strengthen investment on airport infrastructure is effective to improve the airport traffic.

From the above model, we can see that the impact on airport traffic comes from many aspects, such as weather, air control, policy factors and other uncertain factors, but also some strong controllable factors, such as airport facilities, other transportation modes. For the uncertain factors, generally cannot be controlled. For controlling factors, first of all, should increase the airport infrastructure investment, and secondly, should be appropriate to adjust the fare to attract more passengers to take the plane.

4. Model’s Improvement and Prospect

After 2010, China’s Civil Aviation Administration has implemented a number of policies to encourage the development of the aviation industry, which makes

Figure 3. Initial residual plot.

Figure 4. Final residual plot.

Table 8. Parameters value of the final regression model.

R^{2} = 1.000; F = 1.8934e+005; p < 0.0016; s^{2} = 0.0630.

significant changes in the development of aviation industry. The relevant data shows the air transport passenger turnover rate dramatically improved compared with previous years since 2010. Therefore, the above model may not be able to completely describe the air transport passenger turnover changes in the trend and the model should be improved. For the introduction of “policy fac-

tors” into the model and dummy variable $D=\{\begin{array}{l}1,\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{after}\text{\hspace{0.17em}}\text{2010}\text{\hspace{0.17em}}\text{year}\\ 0,\text{\hspace{1em}}\text{before}\text{\hspace{0.17em}}2010\text{\hspace{0.17em}}\text{year}\end{array}$ will be introduced, the model becomes

${Y}_{i}={\beta}_{1}+{\beta}_{2}{X}_{2i}+{\beta}_{3}{X}_{3i}+{\beta}_{4}D{X}_{2i}+{\beta}_{5}D{X}_{3i}+{\epsilon}_{i}$

After computed by Matlab, we obtain an interpretation report below in the case of considering policy factors:

${\stackrel{^}{Y}}_{i}=-654.528+0.258{X}_{2i}+2.345{X}_{3i}-0.303D{X}_{2i}-1.958D{X}_{3i}$

Having been tested, all tests of the model can pass, and the model is well fitted. Therefore, the model after adding the dummy variable is a better model, and the improvement of the model is successful.

It can be seen that the model has a great deal of flexibility, and it can coordinate the arrival and departure processes in each time interval on the airport and make full use of the airport capacity. At the same time, the model also consider the weather factors, which can facilitate decision-makers to combine the meteorological department information with the expected flow in the airport within a next period of time to initially forecast future feasible traffic allocation. This provides a reliable basis for its future scheduling strategy.

5. Conclusions

This paper firstly establishes a method to assess the degree of airport delays and gets all factors which caused the flight delays quantification, and ultimately offers a proposal to deal with the flow factor, which is the principal contradiction for flight delays.

For Question 1, this paper uses the fuzzy comprehensive evaluation method to construct a flight degree evaluation model due to different cognition for different flight level and different people. We choose six indicators to evaluate the degree of airport delay consisting of delay time, the value of flight delay rate change, rate of flight delays, airport delays, weather, and air control. Next we use the analytic hierarchy process to determine the index weight and use triangular membership functions to determine the fuzzy level evaluation indexes and calculate the membership degree of each index. Finally, we get the degree of Shanghai Pudong airport flight delays is level 5, which belongs to the serious delay condition.

For Question 2, based on the grey correlation model, we analyze the various reasons which are correlated with flight delays and divide them into primary and secondary influence on these reasons. Firstly, we initialize the data from 2009 to 2014 and calculate the correlation of sequence based on grey system software. The closer correlation of it, the greater contact it is. Finally, we get the first three main reasons which influence the flight delay are flow, airlines and weather conditions.

For Question 3, starting from the flow factor, which is the most important factor from problem 2, we analyze the airport infrastructure, other transportation ways and policy factors which influence the flow of airport from the aspect of qualitative .We also establish a regression analysis model from the perspective of quantitative and get it tested. Finally, we offer some suggestions to improve the airport traffic such as increasing investment on airport infrastructure, selling ticket at a discount, streamlining government and delegating authorities.

Acknowledgements

This work was supported by the Key project of Education Department for Hubei Province, China (No: D20172905).

References

[1] Civil Aviation Administration of China (2012) Civil Aviation Industry Statistical Development Bulletin. China Civil Aviation Press, Beijing.

[2] Etschmaier, M.M. and Mathaisel, D.F.X. (1985) Airline Scheduling: An Overview. Transportation Science, 5, 36-38.

https://doi.org/10.1287/trsc.19.2.127

[3] Barnhart, C. and Com, A. (2004) Commissioned Paper Airline Schedule Planning. Manufacturing Service Operations Management, 1, 56-66.

[4] Wei, W. and Hanse, M. (2005) Impact of Aircraft Size and Seat Availability on Airlines’ Demand and Market Share in Duopoly Markets. Transportation Research Part E, 41, 315-327.

[5] Lohatepanont, M. and Barnhart, C. (2004) Airline Schedule Planning: Integrated Models and Algorithms for Schedule Design and Fleet Assignment. Transportation Science, 38, 19-32.

https://doi.org/10.1287/trsc.1030.0026

[6] Sherali, H.D., Bish, E.K. and Zhu, X. (2006) Airline Fleet Assignment Concepts, Models, and Algorithms. European Journal of Operational Research, 172, 1-30.

[7] Vranas, P.B., Bertsimas, D.J., Odoni, A.R. (1994) Dynamic Ground-Holding Policies for a Network of Airports. Transportation Science, 28, 275-291.

https://doi.org/10.1287/trsc.28.4.275

[8] Brunetta, L., Guastalla, G. and Navazio, L. (1998) Solving the Multi-Airport Ground Holding Problem. Annals of Operations Research, 81, 271-288.

https://doi.org/10.1023/A:1018909224543

[9] Zou, B. and Hansen, M. (2005) Flight Delay Impact on Airfare and Flight Frequency: A Comprehensive Assessment. Transportation Research Part E, 69, 54-74.

[10] Gopalan, R. and Talluri, K.T. (1998) Mathematical models in Airline Schedule Planning: A survey. Annals of Operations Research, 76, 155-185.

https://doi.org/10.1023/A:1018988203220

[11] Cheung, R.K.-M. and Powell, W.B. (1996) Models and Algorithms for Distribution Problems with Uncertain Demands. Transportation Science, 30, 43-59.

https://doi.org/10.1287/trsc.30.1.43

[12] Barnhart, C. and Cohn, A. (2004) Airline Schedule Planning: Accomplishments and Opportunities. Manufacturing & Service Operations Management, 6, 3-22.

https://doi.org/10.1287/msom.1030.0018

[13] Pita, J.P., Adler, N. and Antunes, A.P. (2014) Socially-Oriented Flight Scheduling and Fleet Assignment Model with an Application to Norway. Transportation Research Part B: Methodological, 61, 17-32.