Next Article in Journal
Electromagnetically Induced Transparency in Media with Rydberg Excitons 2: Cross-Kerr Modulation
Previous Article in Journal
A Novel Color Image Encryption Algorithm Based on Hyperchaotic Maps and Mitochondrial DNA Sequences
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Probabilistic Strategic Conflict-Management for 4D Trajectories in Free-Route Airspace

by
Javier Alberto Pérez-Castán
*,
Álvaro Rodríguez-Sanz
,
Luis Pérez Sanz
,
Rosa M. Arnaldo Valdés
,
V. Fernando Gómez Comendador
,
Clemence Greatti
and
Lidia Serrano-Mira
Higher Technical School of Aeronautical and Space Engineering (ETSIAE), Polytechnic University of Madrid, Plaza Cardenal Cisneros, 28040 Madrid, Spain
*
Author to whom correspondence should be addressed.
Entropy 2020, 22(2), 159; https://doi.org/10.3390/e22020159
Submission received: 26 November 2019 / Revised: 18 January 2020 / Accepted: 27 January 2020 / Published: 30 January 2020
(This article belongs to the Section Multidisciplinary Applications)

Abstract

:
The expected growth of air traffic in the following decades demands the implementation of new operational concepts to avoid current limitations of the air traffic management system. This paper focuses on the strategic conflict management for four-dimensional trajectories (4DT) in free-route airspace. 4DT has been proposed as the future operational concept to manage air traffic. Thus, aircraft must fulfil temporary restrictions at specific waypoints in the airspace based on time windows. Based on the temporary restrictions, a strategic conflict management method is proposed to calculate the conflict probability of an aircraft pair (that intersects in the air) and to calculate temporary-blocking windows that quantify the time span at which an aircraft cannot depart because one conflict could occur. This methodology was applied in a case-study for an aircraft pair, including the uncertainty associated with 4DT. Moreover, a sensitivity analysis was performed to characterise the impact of wind conditions and speed control on the temporary-blocking windows. The results concluded that it is feasible to propose 4DT strategic de-confliction based on temporary-blocking windows. Although, uncertainty variables such as wind and speed control impact on the conflict probability and the size of the temporary-blocking windows.

1. Introduction

The continuous growth of air traffic is expected to double their present values around 2040 [1]. The current system of air traffic management (ATM) presents several inefficiencies and will not be able to manage the air traffic increase. The European macro programme Single European Sky ATM Research (SESAR) addresses the development and implementation of the future ATM system [2]. Among other operational concepts, this paper focuses on the trajectory-based operations (TBO) or four-dimensional trajectories (4DT) and their implementation on a free-route airspace (FRA).
The 4D trajectory of an aircraft consists of the three spatial dimensions plus time as a fourth dimension, defining each waypoint by position (latitude, longitude and flight level) and time [3]. This operational framework foresees trajectory restrictions regarding time: users will follow trajectories without practically any geographical limitations (according to free-route trajectories), as long as they comply with the planned times. Thus, the main objective of this approach is to ensure an optimal trajectory (free-route trajectory whenever would be possible) by “forcing” the aircraft to accurately meet the time of arrival at designated waypoints in exchange [4]. Most of the proposed methods to model 4DT and solve the aircraft trajectory prediction problem can be categorised into two approaches: deterministic and probabilistic [5]. The traditional approach is based on a deterministic formulation, a mathematical problem that describes the aircraft motion. If some external forces or parameters (e.g., aircraft performance, weather conditions, navigation systems accuracy, traffic regulations) are unknown or cannot be precisely appraised, the probabilistic approach transforms the problem into a stochastic one [6,7].
The TBO concept involves separating aircraft via strategic (long-term) trajectory definition, rather than the currently-practised tactical (short-term) conflict resolution [8]. This aims to increase the air traffic capacity by reducing the controllers’ workload. Nevertheless, real-time measures (over the trajectory) will be required to improve reliability, react to unplanned conditions and thus maintain the expected capacity [9].
Numerous studies focus on conflict detection and resolution in a short-term horizon based on collision avoidance manoeuvres. Several techniques were developed for conflict detection and the projection of the further status of both aircraft; the authors recommend [10,11] for a full understanding of the techniques. Moreover, more complex solutions were developed by different authors based on mathematical techniques [12,13,14]. Particularly, conflict detection in free-flight is not a novel concept [15,16]. Paielli and Erzberger [17] developed the concept of trajectory specification based on pairwise conflict detection and resolution. The authors paved the way with this work for further assessment of 4DT. Barnier and Allignol [18] presented an optimisation problem to estimate the cost of directly solving conflicts in the upper airspace with delays. However, there is few research that focuses on mathematical criteria and probabilistic assessment to appraise the viability of the introduction of 4DT in a FRA.
Typically, conflict detection was tackled based on a static approach. Typically, there exists a cylinder-shaped “forbidden volume or protected airspace zone” around the aircraft defined longitudinally by 5 Nautical Miles (NM) and vertically by 1000 feet (ft.). The conflict area or critical section refers to the whole airspace near an intersection where an aircraft can infringe the separation minima depending on the relative position of other aircraft. Geisinger [19] considered an airspace intersection geometry without any geometrical boundaries. International Civil Aviation Organisation ICAO [20] denoted a conflict area around the intersecting point shaped like a rectangle. Netjasov [21] transformed this conflict area to a critical section as a circle around the intersecting points. Although every geometrical concept is valid, the problem arises when uncertainties are associated with the aircraft movement (stochastic approach). Among different solutions, one of them at the present time is the blocking-area concept [22,23]. The blocking-area allows knowing if an aircraft influences the trajectory of other aircraft and may imply a modification of the other trajectory. When the blocking-area transforms its application from location to time, it evolves to a temporary-blocking window [24,25,26]. In this case, the temporary-blocking window provides the time span an aircraft cannot enter in the critical section because one conflict will arise. Then, it provides temporary restrictions to the operation of an aircraft. The main question arises when it is demanded to know if it is possible to specify temporary restrictions between aircraft departing from different airports based on 4DT for the future scheduling process.
This research provides a probabilistic assessment of the incompatibilities between 4DT depending on the city-pairs. The goal is to provide the temporary-blocking windows on the strategic horizon that limits the departures between different airports. This concept tries to solve the typical approach of solving the conflicts in-flight in order to avoid departure delays [17,27]. However, air traffic is expected to double their values in the following 20 years and, in addition, there is a lack of capacity because of air traffic control (ATC) limitations [28]. The integration of 4DT together with temporary-blocking windows can ensure conflict-free trajectories and, in turn, increase capacity because of a reduction of ATC workload. Therefore, this paper tackles the safe integration of 4DT by strategic conflict management. This approach analyses the impact of uncertainty of 4DT in terms of conflict probability according to temporary restrictions that must be fulfilled by aircraft. A further goal is to calculate the temporary-blocking windows that block the departure of an aircraft depending on another aircraft operating 4D trajectories.
The remainder of this work is structured as follows. Section 3 introduces the methodology to tackle the problem. Section 3 details the probabilistic assessment for the strategic conflict management and Section 4 describes the process to analyse the conflict probability. Section 5 presents the results in a real case evaluating ideal results with the inclusion of uncertainty parameters. Lastly, Section 6 provides the primary conclusions.

2. Methods

2.1. Path Modelling for Free-Route Airspace

Currently, airspace users (AUs) define the flight profile based on the predefined routes and waypoints. The operational concept of FRA seeks the path optimisation depending on the AU needs. In this way, the AUs can select the trajectory that best fulfils its operational requirement according to fuel consumptions, aviation charges or flight duration. Herein, it is assumed that full airspace is FRA and, therefore, aircraft will fly according to the shortest route (orthodromic routes) between city-pairs. Equation (1) provides the orthodromic route:
s i n ( λ λ 0 ) = t a n ( θ 0 ) t a n ( φ )
where λ and φ are the longitude and latitude of each point, θ 0 is the path course at the equator, and λ 0 the latitude at the equator. In this way, the lateral path is calculated based on the geographical coordinates of the city-pair. Although it is considered the shortest route, the straight route is not the flight with the shortest flying-time in aviation because the shortest flying-time flight is the maximum circle arc. Optimised or loxodromic routes could be as well introduced depending on the AU priorities because the information required for the conflict-detection model is just the 4DT.
According to the vertical profile, the trajectory is defined based on a simple model that entails: a continuous climbing up to Top of Climb (ToC), a cruise level phase (constant speed is assumed) and a continuous descending phase from the Top of Descent (ToD).
Moreover, the 4D trajectory is not characterised by pre-defined waypoints, but it demands new waypoints or “control points” based on distance-separation criterion [3]. Herein, N waypoints are fixed for time-windows requirements of 4DT. They are separated by a constant distance l that is being currently considered from 150 to 200 NM because the path degradation can be assumed [4]. The N waypoints are uniformly distributed from the ToC to the ToD.
The geographical coordinates of the new waypoints are calculated based on Equation (2).
s i n ( l ) = s i n ( φ j φ i ) c o s ( θ i )
where the indicators i and j denote two consecutive waypoints.
Figure 1 shows that l is the distance between the two consecutive waypoints, t is the time, and t T o C , t W P i and t T o D are the time of the aircraft at the ToC, waypoint i and ToD. The generation of trajectories is based on BADA performance model [29].

2.2. Conflict-Detection Model

The goal of this module is to develop the conflict-detection model for 4DT by an aircraft pair flying intersecting city-pairs. This model determines the temporary-blocking windows that should apply between departures, according to both trajectories. The temporary-blocking windows calculate the required on-ground time-based separation to ensure separation assurance at the crossing point. This model only considers the case of aircraft departing from different airports. This limitation implies that trailing conflicts are out of the scope of this research. Nonetheless, the authors recommend the work of [21] and [30] that propose different methods to solve this type of conflicts.
Conflict is the safety metric considered herein because it is one of the most used worldwide. However, to detect the number of conflicts and then the conflict probability means several issues because different features of the conflicts are not considered. Other safety metrics considered by different authors [20,26,31,32,33] are out of the scope of this problem. In future work, to bear in mind other safety metrics could improve the risk assessment of 4DT.

2.2.1. Critical Section

A conflict is defined as the situation where a separation minima infringement will occur. This infringement is modelled herein based on the critical section [34,35]. The critical section ( d c r i ) is located around the intersection of two trajectories depending on the separation minima S m i n and the crossing angle γ .
d c r i = S m i n s i n ( γ )
The best situation is the geometry with a crossing angle γ = 90 ° . As the crossing angle deviates from 90º, both the critical section and the conflict probability increases. Netjasov [21] already investigated the way it can affect and concluded that the crossing angle is a crucial variable to take into account.

2.2.2. Temporary-Blocking Window

The operational concept of the temporary-blocking window is it calculates on-ground time-based separation on-ground for an aircraft pair ensuring the separation assurance at the crossing point. Based on the position of the aircraft m , there are two critical situations regarding the initial or final points of the temporary-blocking windows. If the aircraft pierces into these two critical situations, a separation minima infringement will occur:
  • The initial critical situation occurs when the aircraft m is at the beginning of the critical section, and the n aircraft is leaving.
  • The final critical situation occurs when the aircraft m is leaving the critical section, and the aircraft n is entering.
Figure 2 represents critical situations.
The initial critical situation is characterised because the aircraft m is located at the instant t e m at the entry point of the critical section. It should be noted that the origin time is the departure time of the aircraft n ( t d e p n ). According to the conflict geometry and the operational parameters of the aircraft n :
t i n t m = t e n + d c r i v n
where t i n t n is the expected time for the aircraft n to be at the intersection, and v n is the cruise ground speed. Going backwards in the trajectory of the aircraft m , the initial departing time ( t d e p i n i m ) can be calculated at which the aircraft m should not depart to not generate a conflict.
t d e p i n i m = t e m t T o C m t c r m t c r m = d ( t e m ) d T o C m v m
where t T o C m is the expected time of the aircraft m to be located at the ToC, t c r m is the time required from the aircraft m to cover the separation between the ToC and the entry point of the critical section ( d ( t e m ) ).
The final critical situation is characterised because the aircraft m leaves the critical section at the time t f m while the aircraft n is piercing into the critical section. According to the geometry of both trajectories, the following condition must be fulfilled:
t f m = t i n t n d c r i v n
The final departing time ( t d e p f i n m ) at which the aircraft m should not depart in order not to generate conflict is:
t d e p f i n m = t f m t T o C m t c r m 2 d c r i v m t c r m = d ( t e m ) d T o C m v m
In this way, the time span [ t d e p m ; t d e p f i n m ] provides the temporary-blocking windows for the aircraft m between the aircraft m and n . In other words, in the case the aircraft m departs inside the temporary-blocking windows, a conflict can occur between both aircrafts at the critical section. The crossing angle affects the size of temporary-blocking windows by d c r i . The greater the value of d c r i , the temporary-blocking windows increases. One issue that should be investigated in further research is the operational application of the temporary-blocking windows. The operational application is crucial because if it considers the application prior to the push-back, the aircraft reduces fuel consumption but new uncertainties are considered. If it considers the location before reaching the departure, it will have a great impact on the queue because of an unfair burden. This research only considers the strategic horizon in which the temporary-blocking windows applies into the allocation of flight plans

3. Probabilistic Models

Once the trajectories for both city-pairs are calculated based on FRA, the uncertainties associated with different operational factors are introduced. First, the trajectories are going to be modelled based on the uncertainty associated with 4DT time windows. This is the most complex step because the aircraft must ensure they are going to ensure time windows at specific waypoints based on predefined trajectories. Second, the aircraft evolution suffers uncertainty throughout the stretch between both waypoints, mainly by the wind and speed. The goal of this work is strategic conflict management; therefore, it only has been considered the uncertainty from the last waypoint before the intersection to the next waypoint. In other words, the previous deviations to the last waypoint before the intersection are introduced in the uncertainty at these waypoints.

3.1. Uncertainty Associated With 4D Trajectories

The introduction of 4DT implies the fulfilment of an agreement between the AUs and the air navigation service providers (ANSPs). This agreement entails the definition of a set of waypoints throughout the trajectories, in which the pilot must pass through them fulfilling a predefined time window of +/- 2 min the 95% of the time [36]. The main implications of this concept are there is a fair distribution of the responsibility between AUs, ANSPs and the Network Manager (NM) in Europe [37]. Initial validation activities and exercises about 4DT concept confirmed that the trajectory must be adjusted between NM and AUs at the strategic level. Once the 4DT is defined, pilots must ensure that they fulfilled the time-windows requirement supported by ATC. Then, the workload can be increased both pilot and ATC at a tactical level depending on the proper development in the strategic and pre-tactical operational phases.
Herein, the previous waypoints to the intersection are denoted as X for the aircraft m and Y for the aircraft n . 4DT uncertainty at the waypoints can be modelled as μ = 0 and σ = 60 to ensure 95% ( 2 σ ) the time required at a specific waypoint. For instance, the time span the aircraft m must fulfil at waypoint X is [ t X m 2 σ ;   t X m + 2 σ ] . Where t X m is the time of the aircraft m at the waypoint X .

3.2. Speed Uncertainty

When both aircraft reach the waypoints X and Y respectively, there is speed uncertainty until the intersection. The planned speed in advance can be different from the operating speed. Once the aircraft crosses a waypoint, it should adjust the speed to fulfil the time windows at the next waypoints. With the aim of seeking different possibilities, two different speed uncertainty modes have been analysed: adjusted and updated speed. These modes are applied independently to both aircraft.

3.2.1. Mode of Adjusted Speed

This mode adjusts the speed to the average value of the 4DT time windows for the next waypoint. In the case the aircraft reaches the waypoint X beforehand or delayed to the estimated average value ( t X ¯ ), the speed up to the next waypoint is modified to achieve the next estimated average value. Therefore, in the case the aircraft reaches before the expected time (( Δ t < 0 ), it will speed down; otherwise, the aircraft will speed up. The adjusted speed v is calculated as follows:
v = l Δ t + l v Δ t = t X ¯ t X
where t X is the real-time of the aircraft at the waypoint X . Moreover, the speed modification is constrained up to 5% nominal speed to avoid disproportionate and non-operating speed modifications during the en-route phase.

3.2.2. Mode of Updated Speed

This mode modifies the planned speed with an updating factor α . This updating factor relates the expected average time at the waypoint with the real crossing time, i.e., it updates the speed ( v ) to the real speed the aircraft has operated, and it is expected to continue flying.
v = α v α = ( t X ¯ t T o C ) ( t X t T o C )
One operational problem of this concept is that speed adjustments increase operational costs for airlines. However, 4DT concept may involve this type of requirements to ensure the achievement of time windows.

3.3. Wind Uncertainty

The wind is the most important source of uncertainty in the planning of trajectory because of the inaccuracy of the forecasts. The wind is modelled following two assumptions: (1) The wind remains constant throughout the trajectory, and (2) the wind varies at different geographical areas. However, one of the assumptions of this study is that we only consider the effect of the wind in the ground speed ( G S ) and not in the path adhesion. The aircraft deviation with regard to the lateral path due to the wind is not considered. Then, the real speed is calculated as:
G S = v W c o s ( β )
where W is the wind and β is the angle of the wind concerning the aircraft heading. Convective weather can generate hardest delays or non-compliance of time windows because aircraft would have to avoid these areas. However, this type of non-standard events are not considered and requires further research.

4. Probabilistic Assessment

The final result of this work is two-fold: to determine the conflict probability in 4DT and to estimate the temporary-blocking windows at which an aircraft should not depart. Figure 3 shows the process to calculate the temporary-blocking windows that allow strategic conflict management with 4D trajectories.
The process to detect conflicts and to determine the temporary-blocking windows for the aircraft m is as follows:
  • It is known the departing time of the aircraft n ( t d e p n ) and the planned 4DT of both aircraft.
  • It is calculated the 4DT time windows of the aircraft n at the closest waypoint ( Y ) to the intersection [ t Y n 2 σ ;   t Y n + 2 σ ] .
  • According to the equations described in Section 2.2, the temporary-blocking window associated with aircraft m at the waypoint X is calculated [ t X m 2 σ ;   t X m + 2 σ ] . This temporary-blocking window is calculated for the probabilistic distribution of the aircraft n based on its 4DT temporary windows.
  • Once the whole temporary-blocking window at the waypoint X is calculated for the aircraft m , this value is extrapolated to the departing time: [ t d e p i n i m ; t d e p f i n m ] .
The temporary-blocking windows for departures do not mean that it is not possible for aircraft m to depart, but it means that there exists conflict probability associated with aircraft n . This process allows determining the conflict probability between both aircraft for specific departing times considering 4DT time windows. The conflict probability is calculated by the number of infringements detected based on the total number of simulations:
P = N o   c o n f l i c t s N o   s i m u l a t i o n s

5. Results and Discussion

5.1. Description of Case Study

Strategic conflict management is applied to a case study between two city-pairs in Spain. The city-pairs are Sevilla (LESV)–Barcelona (LEBL) and Vigo (LEVX)–Murcia (LELC), which are separated by similar distance and aircraft operationally and statistically fly at similar flight levels. Table 1 provides the main parameters of the case study. The value of the different times is provided as minutes: seconds.
Figure 4 shows the 2D representation of the city-pairs Sevilla-Barcelona and Vigo-Murcia.
It can be observed at both trajectories that they are constituted by four waypoints at which the 4DT temporary restrictions must be fulfilled. The distance between the waypoints is l = 100 NM. Waypoints X and Y (the last waypoints before the intersection) are highlighted in Figure 4, as well as the intersecting point. The geographic coordinates of these three points are:
X { λ X = 3.9461 φ X = 38.4566 Y { λ Y = 3.9471 φ Y = 39.6966 I n t { λ i n t = 2.822228 φ i n t = 39.028519

5.2. Conflict-Detection Management without Uncertainty

For the case study of Section 5.1, the temporary-blocking windows were calculated to the departures of the city-pair Sevilla–Barcelona. In this case, no 4DT uncertainty was considered, either wind or speed. The static temporary-blocking window for departures was [21:43; 24:34], a time span of 171 s. This was the baseline to be compared with the introduction of the probabilistic assessment.

5.3. Probabilistic Assessment

The next step was to introduce the uncertainty models described in Section 3 and to analyse how they modified the static temporary-blocking window. 10 4 combinations assessed the integration of both trajectories because of the introduction of 4DT time windows.

5.3.1. 4DT Uncertainty

In this section, the 4DT time windows were applied to calculate the new temporary-blocking window. By applying the process described in Section 4, the new temporary-blocking window based on 4DT uncertainty was [20:43; 25:31] for the 95% (2 σ ), an increase of 112 s with the static temporary-blocking window. Figure 5 represents the probabilistic distribution obtained for the MC simulations.
The conflict probability between both aircraft depended on the combination of both departing times. The average value of the conflict probability distribution is very high, exceeding 80% conflict probability. Therefore, it can be concluded that the introduction of 4DT trajectories (fulfilling temporary restrictions of 2 min at specific waypoints) can be introduced with a temporary-blocking window of 283 s to avoid conflict between this city-pair. This technique allows calculating the conflict probability based on 4DT and temporary requirements and the calculation of temporary-blocking windows for de-conflicting 4DT based on their departing times. The main concern of this study will be to ensure the application on the real operation and the validation of these results.

5.3.2. Speed Uncertainty

As was explained in Section 3.2, two different modes of speed uncertainty were assumed to analyse their impact on 4DT: adjusted and updated speed. Table 2 provides the temporary-blocking windows associated with each mode.
The following conclusions could be extracted:
  • The time span of the temporary-blocking window for the updated speed was 302 s, very similar to the 4DT uncertainty without speed uncertainty. Therefore, this technique increased the size of the temporary-blocking window and deteriorated the previous value of conflict-probability.
  • The time span of the adjusted speed decreased to 166 s. This value was similar to the static case where no uncertainty was considered. The range of the speeds needed to correct the time windows were minor than 5%. This mode provided the most accurate results for conflict probability considering uncertainty. In this case, the pilot would adjust the speed once the aircraft reached a waypoint to ensure the time requirement of the following waypoint. The major limitation of this mode is that the aircraft had to modify its speed throughout the trajectory.
Therefore, the introduction of the adjusted-speed mode with 4DT improves the path predictability to the static temporary-blocking window. In other words, the requirement of adjusting the speed once a waypoint is reached provided the assurance of the same temporal requirements of the case without uncertainty. It could be possible the strategic conflict management depending on the temporary-blocking windows for departures by ensuring a requirement of adjusting the speed at the waypoints. However, the validation of these theoretical results for the path adherence could not be confirmed until 4DT will fly with these requirements.

5.3.3. Wind Uncertainty

The last analysis focused on the way the wind uncertainty affects 4DT in terms of strategic conflict management. As was mentioned before, two different evaluations were performed depending on the wind that was considered as uniform throughout the scenario, or it varied in different geographical areas.

Constant Wind

The first analysis assumed a constant wind throughout the trajectories. This simplification can be used in typical non-convective days where the wind variation is reduced. It was simulated a uniform wind for every direction, and the intensity varied from 0 to 40 knots (kts). These values were typical wind values that could appear for similar altitude conditions. Besides, wind direction varied every 20°. Table 3 shows the temporary-blocking windows for different cases.
From previous results, the temporary-blocking windows varied their time span and location for each combination of intensity and direction:
  • The variation of the time span increased with the wind intensity, although the increase was limited until 3.1% (9 s). Therefore, the variation of the wind intensity below 40 kts can be neglected with respect to the size of the temporary-blocking windows with 4DT.
  • However, the temporal displacement of the temporary-blocking windows cannot be neglected. Those variations could generate a delay or advance for the same wind direction up to 168 s of delay (wind direction 140°) and 138 s of advancement (direction 320°). The size of the 4DT temporary-blocking windows without speed uncertainty is 283 s (Section 5.3.1), this implied that there were temporal displacements over 40%. In the case the static temporary-blocking windows were considered, the displacement could reach up to 100%.
Lastly, Figure 6 shows the variations of the temporary-blocking windows for the different combinations of intensity and direction compared with the reference case (0 kts). The temporary-blocking windows were delayed when the points are located outside of the envelope of 0 kts and advanced when they were inside of it. It was observed that the cross from advancement to delay was roughly obtained for 50° and 240°. Therefore, winds with directions from 50° to 240° advanced the temporary-blocking window and winds with directions from 240° to 50° delayed the temporary-blocking window. Nonetheless, these wind directions mainly depended on the city-pair selected and will affect differently for each city-pair.
According to the above results, the wind uncertainty affected the location of the temporary-blocking window but not their time span. To minimise the wind impact to the temporary-blocking windows, the wind direction should be known in advance with an accuracy of +/-20°. Therefore, the wind direction is a crucial factor for the calculation of temporary-blocking windows in strategic conflict management.

Variable Wind

The analysis for the influence of variable wind throughout the trajectories depended on the geographical areas they have flown. The airspace was divided into four geographical areas (see Figure 7) with a similar length around 100 NM where the wind typically does not exceed 10% variations. Hence, geographical areas 1 and 4 did not take part in the calculation of the uncertainty because the path stretch to study were located at geographical areas 2 and 3. Throughout geographical areas 1 and 4, the aircraft should take the necessary actions to ensure the 4DT time windows.
The variation of the wind in geographical areas 2 and 3 was modelled with a variation of the direction of ± 5° and intensity of ± 10% from the geographical area 1. Although the wind variations were the same in module, the uncertainty was calculated independently for geographical area 2 and 3.
Table 4 shows the results of the temporary-blocking windows for variable wind considering different wind directions. For the sake of clarity, Table 4 only shows the results for a nominal wind of 30 kts (which represents the worst scenario). The column difference indicates the variation in module of the case of variable wind compared with constant wind, and the column displacement indicates the variation of the location between the cases of variable wind compared with constant wind.
The variation of the difference in module between the time span of the temporary-blocking windows did not exceed three seconds. Then, the influence of 10% wind variation between consecutive areas could be neglected. This result agreed with the results of the constant wind effect. About the displacement of the temporary-blocking windows, the effect of the variable wind implied a movement to both sides depending on the uniform wind. Moreover, this displacement did not exceed 5 s, the displacement could be neglected because it modified the size of the temporary-blocking windows less than 2%. However, it should be studied in further work the inclusion of a constant factor in order to take the wind variation into consideration.
Finally, these results concluded that the introduction of temporary-blocking windows could be a reliable and trustworthy system for strategic conflict management based on 4DT. The temporary-blocking windows were a powerful tool for air traffic management to ascertain departing restrictions between city-pairs. The implications could imply a reduction in the ATC workload by the reduction of conflict probability, but, it will demand to pilot modifications to ensure 4DT time windows that will mean pilot workload increase. Then, the model developed herein and the results achieved require further study. There are levels of error in the calculation of the temporary-blocking windows that should be reduced. On the one hand, the temporary-blocking windows were conservative because there were situations assumed as conflicts, but a more detailed assessment could indicate that there was no conflict. The strategic conflict management required knowledge of wind that cannot be ascertained at this strategic horizon; then, it is necessary to introduce new techniques to include more real wind conditions. Lastly further work will include the integration of temporary-blocking windows in the whole European path network; in other words, to tackle the multi-aircraft conflict-detection problem based on 4DT and the optimisation of the results.

6. Conclusions

In this paper, the authors presented the results of a new methodology to deal with the strategic conflict management based on 4D trajectories, while taking a free-route scenario into account. The proposed methodology consisted of the calculation of conflict probability and temporary-blocking windows between two city-pairs. The temporary-blocking windows provided temporary time spans that blocked the departure of one aircraft because of the potential separation infringement in the air. The introduction of 4DT implied that aircraft must fulfil temporary restrictions at specific waypoints throughout the trajectories (4DT time windows). Moreover, the probabilistic assessment took into account the main uncertainties because of speed and wind effect. Speed was modelled based on two models: (1) Updated speed that recalculated the nominal speed based on the time the aircraft passes through the waypoint, and (2) adjusted speed that recalculated the nominal speed to fulfil the estimated average time at the next waypoint. The wind was modelled based on two assumptions as well: (1) The wind was constant throughout the trajectories, and (2) the wind varied throughout the trajectories depending on different geographical areas. The results of the probabilistic assessment were founded in 10 9 Monte Carlo simulations for each experiment. Particularly, this research only considered conflict probability because it is the most widespread safety metric. The main concern of this safety metric is that it does not take into account the specific features of the separation infringements as the severity, duration or minimum distance. Although there is no unique and accepted safety metric, different safety metrics should be considered in further work to deepen in the characteristics of the conflicts, and to bring to the light new insights and methods to tune up conflict management.
First, the length (time span) of the temporary-blocking windows increased with the introduction of uncertainty. The introduction of 4DT time windows implied a temporary-blocking window of almost 300 s while the static approach was 180 s. The introduction of the speed uncertainty provided different results depending on the mode analysed. Updated speed provided an increase of the temporary-blocking windows because of the uncertainty and modifications of the speed. However, adjusted mode reduced the temporary-blocking windows because the aircraft adjusted the speed to fulfil the estimated average time at the next waypoint. This case improved the results and reduced the uncertainty until similar values to the static case. Second, the introduction of wind did not influence the length of the temporary-blocking windows, although the variation of the intensity encompassed from 0 to 40 kts. However, the wind displaced the temporary-blocking windows by the advancement or delay from its original time location. This displacement was bigger depending on the wind direction and the wind intensity.
Therefore, the solution proposed in this work was the implementation of temporary-blocking windows between aircraft of different city-pairs. This solution would imply scheduling restrictions during the strategic phase based on 4DT, and it would demand a collaboration decision-making procedure between the airspace users, the ANSPs and the network manager. About the length of the temporary-blocking windows, the best solution demanded to minimise the conflict probability and the strategic conflict-management by the usage of the speed adjusted mode. This concept ensures that the pilot fulfils the temporary requirements for the 4D time windows, which, in turn, will increase the pilot workload because of the path modifications demanded. The main concern will be to ensure the application of this concept on the future real operation and the validation of these results. Further research should focus on the analysis of this study in the whole European network and the specifications to comply between AUs, ATC, ANSPs and NM.

Author Contributions

J.A.P.-C. was the leader of this research by developing the framework and ideas that will underlie the study. In addition, he was the main contributor to this paper both its redaction and reviewing process. Á.R.-S. helped on the concept development about the integration of 4D trajectories in real aerospace and supported also the writing process. L.P.S. worked during the concept development about the potential issues that could arise by the integration of 4DT in free-rout airspace and during the reviewing process. V.F.G.C. and R.M.A.V. supported the development of the project with their conceptual design. C.G. developed the software and the different technical aspects of the study, helped by L.S.-M. that supported the reviewing process. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Eurocontrol. EUROCONTROL Seven-Year Forecast February 2018. Available online: https://www.eurocontrol.int/publication/seven-year-forecast-flight-movements-and-service-units-february-2019 (accessed on 30 March 2019).
  2. SESAR JU. SESAR Solutions Catalogue 2019, 3rd ed. 2019. Available online: https://www.sesarju.eu/index.php/newsroom/brochures-publications/sesar-solutions-catalogue (accessed on 11 June 2019).
  3. SESAR Joint Undertaking. SESAR 2020 Concept of Operations. 2017. Available online: https://www.sesardeploymentmanager.eu/wp-content/uploads/2016/12/421968-SESAR-Deployment-Programme-edition-2017-Hfdst-1.pdf (accessed on 22 February 2019).
  4. Rodríguez-Sanz, Á.; Comendador, F.G.; Valdés, R.M.A.; Pérez-Castán, J.A.; García, P.G.; Godoy, M.N.G.N. 4D-trajectory time windows: Definition and uncertainty management. Aircr. Eng. Aerosp. Technol. 2019, 91, 761–782. [Google Scholar] [CrossRef]
  5. Matsuno, Y.; Tsuchiya, T. Stochastic 4D Trajectory Optimization for Aircraft Conflict Resolution. In Proceedings of the 2014 IEEE Aerospace Conference, Big Sky, MT, USA, 1–8 March 2014. [Google Scholar]
  6. Margellos, K.; Lygeros, J. Toward 4-D Trajectory Management in Air Traffic Control: A Study Based on Monte Carlo Simulation and Reachability Analysis. IEEE Trans. Control Syst. Technol. 2013, 21, 1820–1833. [Google Scholar] [CrossRef]
  7. Casado, E.; Goodchild, C.; Vilaplana, M. Identification and Initial Characterization of Sources of Uncertainty Affecting the Performance of Future Trajectory Management Automation Systems. Presented at 2nd International Conference on Application and Theory of Automation in Command and Control Systems, London, UK, 29–31 May 2012. [Google Scholar]
  8. Chaimatanan, S.; Delahaye, D.; Mongeau, M. Aircraft 4D Trajectories Planning under Uncertainties. In Proceedings of the 2015 IEEE Symposium Series on Computational Intelligence (SSCI 2015), Cape Town, South Africa, 7–10 December 2015; pp. 51–58. [Google Scholar] [CrossRef]
  9. SESAR. SESAR Concept of Operations Step 2 Edition 2014 (Ed. 01.01.00); SESAR Joint Undertaking: Brussels, Belgium, 2014. [Google Scholar]
  10. Krozel, J.; Peters, M.E.; Hunter, G. Conflict Detection and Resolution for Future Air Transportation; TR97138-01; NASA: Washington, DC, USA, 1997.
  11. Kuchar, J.K.; Yang, L.C. Conflict Detection and Resolution, Air Traffic Control, Alerting Systems, Warning Systems. IEEE Trans. Intell. Transp. Syst. 2000, 1, 179–189. [Google Scholar] [CrossRef] [Green Version]
  12. Erzberger, H.; Heere, K. Algorithm and Operational Concept for Resolving Short-Range Conflicts. Proc. Inst. Mech. Eng. Part G J. Aerosp. Eng. 2010, 224, 225–243. [Google Scholar] [CrossRef]
  13. Lehouillier, T.; Nasri, M.I.; Soumis, F.; Desaulniers, G.; Omer, J. Solving the Air Conflict Resolution Problem Under Uncertainty Using an Iterative Biobjective Mixed Integer Programming Approach. Transp. Sci. 2017, 51, 1242–1258. [Google Scholar] [CrossRef] [Green Version]
  14. Sunil, E.; Ellerbroek, J.; Hoekstra, J.M.; Maas, J. Three-Dimensional Conflict Count Models for Unstructured and Layered Airspace Designs. Transp. Res. Part C Emerg. Technol. 2018, 95, 295–319. [Google Scholar] [CrossRef]
  15. Alam, S.; Shafi, K.; Abbass, H.A.; Barlow, M. An Ensemble Approach for Conflict Detection in Free Flight by Data Mining. Transp. Res. Part C Emerg. Technol. 2009, 17, 298–317. [Google Scholar] [CrossRef]
  16. Paielli, R.A.; Erzberger, H. Conflict Probability Estimation for Free Flight. J. Guid. Control. Dyn. 1997, 20, 588–596. [Google Scholar] [CrossRef]
  17. Paielli, R.A.; Erzberger, H. Trajectory Specification for Terminal Air Traffic: Pairwise Conflict Detection and Resolution. In Proceedings of the 17th AIAA Aviation Technology, Integration, and Operations Conference, Denver, CO, USA, 5–9 June 2017. [Google Scholar]
  18. Barnier, N.; Allignol, C. 4D-Trajectory Deconfliction Through Departure Time Adjustment. Available online: https://hal-enac.archives-ouvertes.fr/hal-00938231/ (accessed on 28 January 2009).
  19. Geisinger, K.E. Airspace Conflict Equations. Transp. Sci. 1985, 19, 139–153. [Google Scholar] [CrossRef]
  20. ICAO. Doc 9689-AN/953 - Manual on Airspace Planning Methodology for the Determination of Separation Minima. Available online: https://searchworks.stanford.edu/view/4026380 (accessed on 18 January 2019).
  21. Netjasov, F. Framework for Airspace Planning and Design Based on Conflict Risk Assessment Part 1: Conflict Risk Assessment Model for Airspace Strategic Planning. Transp. Res. Part C Emerg. Technol. 2012, 24, 190–212. [Google Scholar] [CrossRef]
  22. Speijker, L.J.P.; Blom, H.A.P.; Bakker, G.J.; Karwal, A.K.; Van baren, G.B.; Klompstra, M.B.; Kruijsen, E.A.C. RASMAR Final Report: Risk Analysis of Simultaneous Missed Approaches on Schiphol Converging Runways I9R and 22; Nationaal Nationaal Lucht- en Ruimtevaartlaboratorium: Amsterdam, The Netherlands, 2001. [Google Scholar]
  23. Henry, M.; Schmitz, S.; Revenko, N.; Kelbaugh, K. A Monte Carlo Simulation for Evaluating Airborne Collision Risk in Intersecting Runways. In Proceedings of the AIAA Modeling and Simulation Technologies (MST) Conference, Boston, MA, USA, 19–22 August 2013. [Google Scholar]
  24. Pérez-Castán, J.A.; Gómez Comendador, F.; Rodríguez-Sanz, A.; Barragán Montes, R.; Arnaldo Valdés, R.; Pérez Sanz, L. Impact of Continuous Climb Operations on Airport Capacity. Transp. Res. Part C Emerg. Technol. 2018, 96, 231–250. [Google Scholar] [CrossRef] [Green Version]
  25. Pérez-Castán, J.A.; Comendador, F.G.; Rodríguez-Sanz, Á.; Barragán, R.; Arnaldo-Valdés, R.M. Design of a Conflict-Detection Air Traffic Control Tool for the Implementation of Continuous Climb Operations: A Case Study at Palma TMA. Proc. Inst. Mech. Eng. Part G J. Aerosp. Eng. 2019, 233, 095441001983119. [Google Scholar] [CrossRef]
  26. Pérez-Castán, J.A.; Comendador, F.G.; Rodriguez-Sanz, Á.; Arnaldo Valdés, R.M.; Agueda, G. RPAS Integration in Non-Segregated Airspace: Safety Metrics for Tactical Planning. Proc. Inst. Mech. Eng. Part G J. Aerosp. Eng. 2019, 233. [Google Scholar] [CrossRef]
  27. Erzberger, H.; Lauderdale, T.A.; Chu, Y.-C. Automated conflict resolution, arrival management, and weather avoidance for air traffic management. Proc. Inst. Mech. Eng. Part G J. Aerosp. Eng. 2012, 226, 930–949. [Google Scholar] [CrossRef]
  28. EUROCONTROL. European Aviation in 2040—Challenges of Growth. 2018. Available online: https://www.eurocontrol.int/publication/challenges-growth-2018 (accessed on 26 January 2020).
  29. EUROCONTROL. User Manual for the Base of Aircraft DAta (BADA), Family 4. 2014. Available online: https://www.eurocontrol.int/publication/user-manual-base-aircraft-data-bada (accessed on 26 January 2020).
  30. Rey, D.; Rapine, C.; Fondacci, R.; Faouzi, N.-E. Subliminal Speed Control in Air Traffic Management: Optimization and Simulation. Transp. Sci. 2016, 50, 240–262. [Google Scholar] [CrossRef] [Green Version]
  31. Piric, S.; De Boer, R.J.; Roelen, A.; Karanikas, N.; Kaspers, S. How Does Aviation Industry Measure Safety Performance Current Practice and Limitations. Int. J. Aviat. Manag. 2019, 4, 224. [Google Scholar] [CrossRef] [Green Version]
  32. Campos, L.M.B.C.; Marques, J.M.G. On a Dimensionless Alternative to the ICAO Target Level of Safety. Proc. Inst. Mech. Eng. Part G J. Aerosp. Eng. 2016, 230, 1548–1557. [Google Scholar] [CrossRef]
  33. Netjasov, F.; Janic, M. A Review of Research on Risk and Safety Modelling in Civil Aviation. J. Air Transp. Manag. 2008, 14, 213–220. [Google Scholar] [CrossRef]
  34. Netjasov, F. Framework for Airspace Planning and Design Based on Conflict Risk Assessment. Part 2: Conflict Risk Assessment Model for Airspace Tactical Planning. Transp. Res. Part C Emerg. Technol. 2012, 24, 213–226. [Google Scholar] [CrossRef]
  35. Pérez-Castán, J.A.; Gómez Comendador, F.; Rodríguez-Sanz, Á.; Arnaldo Valdés, R.M. Conflict-Risk Assessment Model for Continuous Climb Operations. Aerosp. Sci. Technol. 2019, 84, 812–820. [Google Scholar] [CrossRef]
  36. EUROCONTROL. “RECAT-EU” European Wake Turbulence Categorisation and Separation Minima on Approach and Departure. 2018. Available online: https://www.eurocontrol.int/publication/european-wake-turbulence-categorisation-and-separation-minima-approach-and-departure (accessed on 26 January 2020).
  37. SESAR Joint Undertaking. Optimised Airspace User Operations. 2016. Available online: https://www.eurocontrol.int/project/optimised-airspace-users-operations (accessed on 26 January 2020).
Figure 1. Vertical trajectory profile.
Figure 1. Vertical trajectory profile.
Entropy 22 00159 g001
Figure 2. At the left initial critical situation; at the right final critical situation.
Figure 2. At the left initial critical situation; at the right final critical situation.
Entropy 22 00159 g002
Figure 3. Process of the strategic conflict management based on 4D trajectories.
Figure 3. Process of the strategic conflict management based on 4D trajectories.
Entropy 22 00159 g003
Figure 4. 2D representation of case-study.
Figure 4. 2D representation of case-study.
Entropy 22 00159 g004
Figure 5. Conflict probability as a function of t d e p S e v i l l a B a r c e l o n a .
Figure 5. Conflict probability as a function of t d e p S e v i l l a B a r c e l o n a .
Entropy 22 00159 g005
Figure 6. Envelope of temporary-blocking windows with 4DT and constant wind uncertainty.
Figure 6. Envelope of temporary-blocking windows with 4DT and constant wind uncertainty.
Entropy 22 00159 g006
Figure 7. Limits of the areas where the wind varies.
Figure 7. Limits of the areas where the wind varies.
Entropy 22 00159 g007
Table 1. Parameters of the case-study.
Table 1. Parameters of the case-study.
City-PairSevilla-BarcelonaVigo-Murcia
Departure5°53′35″ W, 37°25′05″ N8°37′36″ W, 42°13′54″ N
Arrival2°04′42″ E, 41°17′49″ N0°48′45″ W, 37°46′30″ N
Mach v 1 = 0.78 v 2 = 0.78
Flight Level350350
Climbing time t T o C 1 = 25:00 t T o C 2 = 26:00
Departure time- t d e p 2 = 00:00
Time at intersection- t i n t 2 = 47:00
Table 2. Temporary-blocking windows with speed uncertainty.
Table 2. Temporary-blocking windows with speed uncertainty.
ModeTemporary-Blocking WindowsTime Span (sec)
Updated speed[20:36; 25:38]302
Adjusted speed[21:44; 24:30]166
Table 3. Temporary-blocking windows for constant wind.
Table 3. Temporary-blocking windows for constant wind.
Wind DirectionIntensity (kts)
051015203040
[20:43; 25:31][20:29; 25:15][20:16; 25:03][19:54; 24:43][19:50; 24:38][19:18; 24:06][18:42; 23:32]
20°[20:31; 25:19][20:26; 25:15][20:15; 25:05][20:05; 24:55][19:36; 24:26][19:13; 24:05]
40°[20:40; 25:28][20:38; 25:28][20:31; 25:20][20:27; 25:18][20:14; 25:06][20:03; 24:55]
60°[20:48; 25:37][20:49; 25:39][20:49; 25:41][20:50; 25:42][20:48; 25:42][20:55; 25:51]
80°[20:54; 25:44][20:53; 25:45][21:12; 25:36][21:18; 26:10][21:35; 26:30][21:52; 26:49]
100°[21:02; 25:50][21:11; 26:01][21:23; 26:14][21:45; 26:39][22:12; 27:07][22:42; 27:38]
120°[20:57; 25:47][21:15; 26:06][21:43; 26:33][21:57; 26:49][22:32; 27:26][23:15; 28:12]
140°[21:05; 25:53][21:32; 26:24][21:44; 26:35][21:59; 26:51][22:47; 27:40][23:31; 28:24]
160°[20:56; 25:46][21:23; 26:12][21:49; 26:41][21:59; 26:49][22:41; 27:31][23:21; 28:13]
180°[20:59; 25:47][21:18; 26:06][21:32; 26:20][21:51; 26:39][22:18; 27:07][22:53; 27:41]
200°[21:00; 25:50][21:07; 25:54][21:17; 26:05][21:23; 26:11][21:49; 26:35][22:09; 26:55]
220°[20:52; 25:40][21:03; 25:15][20:57; 25:43][20:58; 25:44][21:15; 26:00][21:20; 26:03]
240°[20:47; 25:35][20:46; 25:32][20:41; 25:27][20:37; 25:21][20:38; 25:21][20:29; 25:10]
260°[20:36; 25:23][20:34; 25:20][20:19; 25:04][20:09; 24:52][19:52; 24:36][19:42; 24:23]
280°[20:29; 25:17][20:23; 25:09][20:04; 24:50][19:58; 24:42][19:24; 24:08][19:05; 23:46]
300°[20:20; 25:10][20:06; 24:52][19:54; 24:39][19:43; 24:29][19:09; 24:51][18:38; 23:20]
320°[20:27; 25:15][20:06; 24:52][19:51; 24:37][19:30; 24:16][18:56; 24:43][18:25; 23:09]
340°[20:30; 25:18][20:09; 24:56][19:55; 24:42][19:34; 24:20][19:08; 24:56][18:26; 23:12]
Table 4. Temporary-blocking windows for constant and variable wind.
Table 4. Temporary-blocking windows for constant and variable wind.
Variable WindDifferenceDisplacement
[19:23; 24:02]2 sec(0.7%)−5 sec
20°[19:54; 24:28]0 sec(0%)2 sec
40°[20:22; 24:04]1 sec(0.3%)−2.5 sec
60°[20:25; 24:46]0 sec(0%)4 sec
80°[21:26; 24:29]0 sec(0%)−1 sec
100°[22:32; 24:04]−1 sec(−0.3%)−2.5 sec
120°[22:28; 24:28]−2 sec(−0.7%)3 sec
140°[22:27; 24:38]−1 sec(−0.3%)−1.5 sec
160°[22:27; 24:31]0 sec(0%)0 sec
180°[22:27; 24:07]−1 sec(−0.3%)0.5 sec
200°[21:26; 24:35]0 sec(0%)0 sec
220°[21:25; 24:55]−1 sec(−0.4%)−4.5 sec
240°[20:25; 24:16]1 sec(0.4%)−5.5 sec
260°[19:54; 24:39]−3 sec(–1.1%)4.5 sec
280°[19:44; 24:11]−2 sec(−0.7%)4 sec
300°[19:23; 24:51]1 sec(0.4%)−0.5 sec
320°[18:53; 24:42]0 sec(0%)−1 sec
340°[18:23; 24:45]−2 sec(−0.7%)−5 sec

Share and Cite

MDPI and ACS Style

Pérez-Castán, J.A.; Rodríguez-Sanz, Á.; Pérez Sanz, L.; Arnaldo Valdés, R.M.; Gómez Comendador, V.F.; Greatti, C.; Serrano-Mira, L. Probabilistic Strategic Conflict-Management for 4D Trajectories in Free-Route Airspace. Entropy 2020, 22, 159. https://doi.org/10.3390/e22020159

AMA Style

Pérez-Castán JA, Rodríguez-Sanz Á, Pérez Sanz L, Arnaldo Valdés RM, Gómez Comendador VF, Greatti C, Serrano-Mira L. Probabilistic Strategic Conflict-Management for 4D Trajectories in Free-Route Airspace. Entropy. 2020; 22(2):159. https://doi.org/10.3390/e22020159

Chicago/Turabian Style

Pérez-Castán, Javier Alberto, Álvaro Rodríguez-Sanz, Luis Pérez Sanz, Rosa M. Arnaldo Valdés, V. Fernando Gómez Comendador, Clemence Greatti, and Lidia Serrano-Mira. 2020. "Probabilistic Strategic Conflict-Management for 4D Trajectories in Free-Route Airspace" Entropy 22, no. 2: 159. https://doi.org/10.3390/e22020159

APA Style

Pérez-Castán, J. A., Rodríguez-Sanz, Á., Pérez Sanz, L., Arnaldo Valdés, R. M., Gómez Comendador, V. F., Greatti, C., & Serrano-Mira, L. (2020). Probabilistic Strategic Conflict-Management for 4D Trajectories in Free-Route Airspace. Entropy, 22(2), 159. https://doi.org/10.3390/e22020159

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop