Research article

An en route capacity optimization model based on air traffic control process

  • Received: 01 January 2022 Revised: 03 February 2022 Accepted: 14 February 2022 Published: 25 February 2022
  • In order to resolve the imbalance of demand-capacity and airspace congestion, improve the performance of the en route air traffic management, promote the development of air traffic control automation system in the future, this paper proposes an En route air traffic control process model from the perspective of operation requirements. Taking the minimization of operation time, instantaneous density, maximum lateral displacement and air traffic controllers' workload as the optimization objectives, the commonly used air traffic control instructions such as climb and descent and speed restriction are set as constraints, the algorithm is designed based on the air traffic control scheme, and a complete air traffic control process are modeled which outputs instructions for each aircraft. Finally, the model is applied to a case study in the northwest region of China. The simulation results show that compared with the actual operation process, the total operation time is reduced by 18.6%, the variance of the lateral displacement and the vertical separation are efficiently reduced, and the en route air traffic capacity is improved. The proposed model envisages the following two innovations: (ⅰ) the whole process of air traffic controllers' command is considered in the model, especially the control scheme and different types of instructions, and (ⅱ) the en route historical trajectory data of aircraft is used to as the key parameters of the input data to efficiently yield the acceptable results of the model. By quantifying the operation requirements of air traffic control, this model can also balance the distribution of traffic flow density, reduce the utilization rate of horizontal airspace, alleviate flight conflicts on air routes, and lessen the workload of controllers.

    Citation: Jie Ren, Shiru Qu, Lili Wang, Yu Wang. An en route capacity optimization model based on air traffic control process[J]. Mathematical Biosciences and Engineering, 2022, 19(4): 4277-4299. doi: 10.3934/mbe.2022198

    Related Papers:

    [1] Bo Sun, Ming Wei, Binbin Jing . Optimal model for the aircraft arrival and departure scheduling problem with fuzzy runway incursion time. Mathematical Biosciences and Engineering, 2021, 18(5): 6724-6738. doi: 10.3934/mbe.2021334
    [2] Lijing Ma, Shiru Qu, Jie Ren, Xiangzhou Zhang . Mixed traffic flow of human-driven vehicles and connected autonomous vehicles: String stability and fundamental diagram. Mathematical Biosciences and Engineering, 2023, 20(2): 2280-2295. doi: 10.3934/mbe.2023107
    [3] Suresh Y, Kalaivani T, Senthilkumar J, Mohanraj V . BF2 VHDR based dynamic routing with hydrodynamics for QoS development in WSN. Mathematical Biosciences and Engineering, 2020, 17(1): 930-947. doi: 10.3934/mbe.2020050
    [4] Lei Jin, Sixiang Lin, Binglei Xie, Lin Liu . A vulnerability-based vehicle routing approach for solving capacitated arc routing problem in urban snow plowing operations. Mathematical Biosciences and Engineering, 2021, 18(1): 166-181. doi: 10.3934/mbe.2021009
    [5] Jin Zhang, Nan Ma, Zhixuan Wu, Cheng Wang, Yongqiang Yao . Intelligent control of self-driving vehicles based on adaptive sampling supervised actor-critic and human driving experience. Mathematical Biosciences and Engineering, 2024, 21(5): 6077-6096. doi: 10.3934/mbe.2024267
    [6] Qixin Zhu, Mengyuan Liu, Hongli Liu, Yonghong Zhu . Application of machine learning and its improvement technology in modeling of total energy consumption of air conditioning water system. Mathematical Biosciences and Engineering, 2022, 19(5): 4841-4855. doi: 10.3934/mbe.2022226
    [7] Matthew Hayden, Bryce Morrow, Wesley Yang, Jin Wang . Quantifying the role of airborne transmission in the spread of COVID-19. Mathematical Biosciences and Engineering, 2023, 20(1): 587-612. doi: 10.3934/mbe.2023027
    [8] Simone Göttlich, Stephan Knapp . Modeling random traffic accidents by conservation laws. Mathematical Biosciences and Engineering, 2020, 17(2): 1677-1701. doi: 10.3934/mbe.2020088
    [9] Haiquan Wang, Hans-Dietrich Haasis, Menghao Su, Jianhua Wei, Xiaobin Xu, Shengjun Wen, Juntao Li, Wenxuan Yue . Improved artificial bee colony algorithm for air freight station scheduling. Mathematical Biosciences and Engineering, 2022, 19(12): 13007-13027. doi: 10.3934/mbe.2022607
    [10] Hongwen Hu, Miao Ye, Chenwei Zhao, Qiuxiang Jiang, Xingsi Xue . Intelligent multicast routing method based on multi-agent deep reinforcement learning in SDWN. Mathematical Biosciences and Engineering, 2023, 20(9): 17158-17196. doi: 10.3934/mbe.2023765
  • In order to resolve the imbalance of demand-capacity and airspace congestion, improve the performance of the en route air traffic management, promote the development of air traffic control automation system in the future, this paper proposes an En route air traffic control process model from the perspective of operation requirements. Taking the minimization of operation time, instantaneous density, maximum lateral displacement and air traffic controllers' workload as the optimization objectives, the commonly used air traffic control instructions such as climb and descent and speed restriction are set as constraints, the algorithm is designed based on the air traffic control scheme, and a complete air traffic control process are modeled which outputs instructions for each aircraft. Finally, the model is applied to a case study in the northwest region of China. The simulation results show that compared with the actual operation process, the total operation time is reduced by 18.6%, the variance of the lateral displacement and the vertical separation are efficiently reduced, and the en route air traffic capacity is improved. The proposed model envisages the following two innovations: (ⅰ) the whole process of air traffic controllers' command is considered in the model, especially the control scheme and different types of instructions, and (ⅱ) the en route historical trajectory data of aircraft is used to as the key parameters of the input data to efficiently yield the acceptable results of the model. By quantifying the operation requirements of air traffic control, this model can also balance the distribution of traffic flow density, reduce the utilization rate of horizontal airspace, alleviate flight conflicts on air routes, and lessen the workload of controllers.



    According to an air traffic forecast scenario that Aviation Industry Corporation of China identified as the most-likely scenario for China, the revenue passenger kilometer (RPK) will increase to 3.3 trillion passenger kilometers by the year 2020, which is 175% more than in 2019 with an average annual growth of 5.3%. It is also predicted that the growth rate of air traffic demand for transport services will outpace the growth rate of the airport infrastructure development and of the available airspace resources [1]. This will result in the imbalance of demand-capacity and flight delays which may cause safety and congestion issues in the air traffic management (ATM) system. Therefore, it is critical to improve the performance of ATM system from the aspects of increase capacity and efficiency of the airspace, such as terminal maneuvering area (TMA) and area control center (ACC). Compared with TMA, the en route traffic congestion in ACC becomes more and more prominent in crowded airspace around the world [2]. Hence, the issue of en route capacity has attracted considerable attention in the academic literature recently. However, most scholars address the evaluation of en route capacity and the route network structure optimization, the issues of en route capacity improvement from the perspective of air traffic control operation are rarely discussed in existing studies. Another motivation for this paper is to study the command process of Air traffic controller (ATCO) to establish an en route capacity optimization model. Air traffic control is currently performed by human controllers using radar displays and voice communication with pilots [3]. The U.S. NextGen (Next Generation Air Transportation System) envisions an ATM system that allocates functions to the human and automation and between different stakeholders such as pilots and controllers in a significantly different manner than today's system. These concepts promise increases in capacity while simultaneously improving safety. For example, in an automation controlled or managed airspace, a computer can perform certain tasks that are needed to separate aircraft and hence the airspace may be able to accommodate a larger amount of traffic [4]. To improve the performance of ATC, mainly in anticipation of increasing levels of air traffic, research effort has been devoted over the last decade on creating tools to assist ATC with trajectory prediction, conflict detection and resolution tasks [5]. However, the best of the author's knowledge, no air traffic automation model with aircraft historical trajectory data has been comprehensively analyzed yet.

    The primary goal of our research is to present an en route capacity optimization model to improve the route capacity and decrease the workload of ATCO, and promote the development of ATC automation system in the long run. The main developments are summarized as follows: (ⅰ) development of an en route capacity optimization model considering the whole process of ATCO command; (ⅱ) the en route historical trajectory data of aircraft is used to as the key parameters of the input data. Finally, a real-world case study is conducted to demonstrate the validity of the proposed approach.

    The remainder of this paper is organized as follows: Section 2 discusses previous related work and establishes the contributions of this paper, Section 3 presents the proposed model, while Section 4 describes the results of the application of the proposed model. Finally, Section 5 summarizes the research conclusions and provides directions for future research.

    Despite the goal of improving route capacity, the focus of this paper is on ATCO workflow-based optimization model, so the literature on air route capacity evaluation and route network optimization will not be discussed in this section.

    Today, the ATC automation system only plays an auxiliary role, and the performance of the air traffic management system mainly depends on the level of the service provided by air traffic controller [6]. The operation state of air traffic flow is directly related to the capacity of aviation system. Before takeoff, aircraft file flight plans which cover the entire flight. During the flight, ATCO sends additional instructions to them, depending on the actual traffic, to improve traffic flow and avoid dangerous encounters. The primary concern of ATC is to maintain safe separation between the aircraft [7]. The control command of the controller is to provide control services, coordinate the flight conflicts caused by the operational requirements of different aircraft, and at the same time ensure the effective use of airspace resources. According to the process of air traffic control, the research of ATC automation system mainly focuses on trajectory prediction, conflict detection and conflict resolution.

    Inseok Hwang et al. designed interactive multiple model (IMM) and residual mean interactive multiple model (RMIMM), and analyzed the conflict detection problem in multi aircraft scenarios based on the model [8]. Liu et al. predicted the transition between flight states by stochastic linear hybrid system [9]. Porretta et al. proposed a kinematic model, which considered the influence of wind and aircraft intention information, and verified the feasibility of the model through real data [10]. Kaneshige et al. also proposed a track prediction method based on kinematic model [11]. The flight process is divided into two stages: horizontal plane and altitude profile. Based on the flight time, fuel consumption and other indicators, the difference between no trajectory prediction and with trajectory prediction was compared. It is proved that trajectory prediction can improve the efficiency of trajectory operation. Tang proposed a flight contour clustering method based on tweed and K-means algorithm [12]. Tastamberkov et al. achieved 4D track fitting prediction by regression method [13]. In view of the fact that neural network can infinitely approximate any complex nonlinear relationship, Shi proposed a track prediction model based on long-term memory (LSTM) network [14]. The four interaction layers of the repeat module in LSTM enable it to connect the long-term dependence to the current prediction task. Zhang et al. proposed to use ant lion optimization (ALO) algorithm to optimize LSTM neural network, which provides a good reference for trajectory prediction [15]. Cui et al. proposed a trajectory prediction model with higher accuracy on the basis of previous track prediction research, aiming at the defects of traditional technology, recurrent and multilayer neural network structure [16]. The algorithm can provide an ideal conflict solution for each aircraft in a short time.

    There has been a lot of research to solve the problem of aircraft conflict detection in the presence of uncertainty. The most common way is to treat the expected positions of the aircraft as random variables and assume that they follow a probability distribution. Paielli and Erzberger modeled the lateral error and longitudinal error of the track as a normal distribution with a constant increase in the variance over the entire time range, and estimated the conflict probability on this basis [17]. Krozel and Peters proposed a non-deterministic conflict detection model for two aircraft, in which the relative position, speed and heading between them are considered to obey a Gaussian distribution [18]. Prandini et al. proposed a probabilistic model for mid-term conflict detection, in which the tracking error obeys a normal distribution, and its variance increases over time [19]. Prandini et al. also proposed a short-term conflict detection model, in which the motion of the aircraft is treated as a deterministic motion plus Brownian motion disturbance. Jacquemart and Morio considered Brownian motion and the variance that grows over time, not only can estimate the probability of aircraft collision, but also use an interactive particle system algorithm suitable for the estimation of the probability of a small probability event to estimate the possibility of aircraft collision [20]. Shi et al. proposed a fast algorithm based on the idea of position space discretization to calculate the conflict probability. And on this basis, a method based on three-dimensional Brownian motion using coordinate transformation and the Bachelier-Levy theorem to estimate the probability of conflict is proposed [21]. Most probabilistic conflict detection methods are based on calculating the conflict probability, that is, the probability that the distance between a pair of aircraft is less than the specified minimum [22,23]. Blom and Bakker proposed and analyzed feature matrices, such as crossover probability and overlap probability [24,25]. However, it is difficult to obtain an accurate calculation result of the conflict rate through analysis, so some studies have adopted geometric approximation or hypothetical methods to solve this problem. For example, the approximation proposed by Paielli and Erzberger is based on a constant speed cross trajectory, but the stability of this method is not good. If one of the aircraft is turning, it may cause the calculation to fail [17,26]. In addition, Liu and Hwang used the classical approximation results for the Gaussian quadratic form through series expansion [27].

    The optimization method is the most studied method, usually combining the kinematics model of the aircraft with a set of cost matrices, and determining the best conflict resolution strategy by solving the trajectory of the maximum cost. The core of this optimization method is a reasonable definition of the cost function. Typical costs include fuel consumption, time consumption, and work load. Mixed integer linear and nonlinear programming are commonly used methods to resolve conflict resolution problems. In order to obtain the optimal trajectory for conflict resolution, Cafieri et al. established two mixed-integer linear programming models with the goal of minimizing the total flight time, each of which allows only one factor of speed or heading to change [28]. Alonso Ayuso et al. proposed a hybrid 0-1 nonlinear model based on geometric structure in order to obtain the best strategy for conflict avoidance of different aircraft configurations. With the goal of minimizing aircraft acceleration changes, the model is solved by linear approximation using iterative Taylor polynomials [29]. Omer proposed a mixed integer linear model based on spatial discretization. With the optimization goal of minimizing fuel consumption and delay time, the continuous speed vector is considered comprehensively, as well as the limitation of speed, acceleration and yaw rate, to ensure the feasibility of the actual planned trajectory [30]. Intelligent optimization algorithms have also been widely used in trajectory planning. Zhou et al. improved the genetic algorithm considering the arrival time constraints, and realized the UAV trajectory planning in a dynamic environment. Because of its strong parallel computing capability and global search capability, it can quickly obtain the global optimal solution [31]. The rolling time domain method completes the real-time update of planning sub-objectives through the advancement of rolling windows, and is often combined with fuzzy logic and swarm intelligence algorithms to improve the efficiency of local trajectory planning. Hu improved the ant colony algorithm for UAV static trajectory pre-planning, and proposed an artificial bee colony algorithm to realize dynamic trajectory planning considering sudden threats [32]. The study found that the introduction of random ant colony can expand the search space and obtain more diverse solutions, while the bee colony algorithm can well complete the local trajectory correction. In addition, Chiang et al. used the spatiotemporal flow method to solve the conflict problem, which uses the graphic search method in discrete time and space to find a conflict free trajectory for the aircraft through iteration [33]. Tomlin et al. proposed a conflict resolution method based on Lie algebra and Hamilton Jacobian equation [34]. Bicchi used optimal control and game theory to solve the problem of aircraft conflict. The model assumes that the linear speed of aircraft is constant, and takes airspeed, course, longitude and latitude as optimization parameters [35]. Goodchild et al. proposed a collaborative optimal conflict resolution algorithm based on distributed artificial intelligence by using dynamic optimization algorithm [36]. Bayen et al. proposed a Lagrange model based on simple conflict resolution instruction, which allows aircraft to fly at different altitudes [37]. Durand et al. established an unsupervised learning neural network to calculate the approximate optimal trajectory of two aircraft in conflict. When collision is detected, the three-layer neural network only changes the course. By using genetic algorithm to help train neural network, this method can solve the conflict of no more than three aircraft [38].

    In summary, the current research primarily paid attention to the efficiency and accuracy of conflict detection, trajectory prediction and conflict resolution in the field of ATC automation system study. However, very little of them considered the key factors of the whole ATC process in actual operation, such as the issue of different restriction instructions and control scheme. In addition, the historical trajectory data of aircraft are not commonly used to deal with the problem of ATC automation system. To this end, this paper aims to propose an en route capacity optimization model from the aspect of ATCO's instruction on account of the characteristics of the historical trajectory data of aircraft.

    In this section, we present an en route capacity optimization model with the consideration of ATC command process to fulfill the following operational requirements: there is no flight conflict on the route; the traffic flow is less than the capacity; the ATCO workload is reduced; the density distribution of aircraft is balanced; and the lateral occupied space is minimized. For the reason that better meet the actual en route operation, the historical trajectory data of aircraft are processed and used as the input of this model.

    In most cases, the airspace of en route flight belongs to controlled airspace in which all the maneuvers of aircraft, i.e., climb, descent, acceleration and deceleration require the instructions of the ATCO. Usually, there is only one ATCO in the designated airspace unit. Before the aircraft arrives at the airspace unit under his/her jurisdiction, ATCO will prepare a preliminary control scheme based on the scheduled arrival time and its operational needs. When the aircraft actually arrives at the airspace unit, ATCO will issue more detailed control instructions according to the current airspace conditions.

    In the whole process of air traffic control, the ATCO needs to consider many factors, coordinate different operational objectives, and avoid potential risks between aircraft. For the convenience of research, it is necessary to simplify the complexity of the real- life problem and facilitate the calculation. Refer to [3,4,7], the following assumptions are made:

    1). Unless otherwise specified, "route" in this paper refers to a route within the airspace under the jurisdiction of the ATCO. And there is no alternative route for the aircraft to fly.

    2). Aircraft should fly along the route direction when entering or exiting the route. The entry and exit process to the route is not considered in this paper.

    3). If ATCO does not issue maneuver instructions, the aircraft will only move in a straight line at a constant speed.

    4). Aircraft will not change its heading and flight level simultaneously.

    5). There is no convective weather or military activities.

    6). The clearance interval is 2 minutes at the exit waypoint of the route in radar control circumstances.

    7). Repeated climb and descend, or repeated acceleration and deceleration instructions are not allowed.

    8). The rate of climb, rate of descent, acceleration and deceleration of aircraft are limited respectively. The average value, maximum value and minimum value are extracted from the historical radar data of the aircraft on the route. The average value is taken as the default value, and the maximum value and minimum value are used to set the limit interval.

    Although the above assumptions were made during the modeling, it did not have much impact on the practicability and validity of the model, mainly because most of the assumptions occurred infrequently during the en-route operation phase of the aircraft. Taking the military activities as an example, if the military activities such as artillery fire around the route are considered, the instructions issued by ATCO will become much more complicated, i.e., lateral displacement is forbidden, or the en route holding procedure will be conducted. However, these situations are not common.

    To facilitate the model elaboration, the notations and definitions used hereafter are summed up in Table 1.

    Table 1.  Parameters and variables in the proposed model.
    Parameters
    Δti Operation time of the ith aircraft on the route
    Ymax Maximum value of lateral displacement of all aircraft
    M Number of ATC instructions issued by ATCO which are different from the control scheme
    var(Sepai) Separation between the ith aircraft and the preceding aircraft
    amax Maximum value of acceleration of aircraft
    amean Average value of acceleration of aircraft
    bmean Average value of deceleration of aircraft
    bmin Minimum value of deceleration of aircraft
    v1 Speed of aircraft before speed adjustment, i.e., initial speed
    v2 Target speed
    gradCmax Maximum value of climb gradient of aircraft
    gradCmean Average value of climb gradient of aircraft
    gradDmean Average value of descent gradient of aircraft
    gradDmin Minimum value of descent gradient of aircraft
    h1 Flight altitude of aircraft before the flight level adjustment, i.e., initial altitude
    h2 Target altitude
    v Velocity during the flight level adjustment
    Ldel Distance the aircraft is delayed
    tdel Time the aircraft is delayed
    ΔL Actual delayed distance of aircraft
    α1 Angle of first turn
    α1+α2 Angle of second turn
    LMK Distance of MK segment
    Dslope Slope distance between the aircraft and other aircraft on the route
    Dlen Longitudinal distance between the aircraft and other aircraft on the route
    Dver Vertical distance between the aircraft and other aircraft on the route
    Xhci Position of the ith aircraft receiving flight level adjustment instructions
    Xvci Position of the ith aircraft receiving speed adjustment instructions
    ai Acceleration of the ith aircraft receiving speed adjustment instructions
    gradi Gradient of the ith aircraft receiving flight level adjustment instructions
    Si Conflict resolution strategy used by the ith aircraft

     | Show Table
    DownLoad: CSV

    In order to improve the efficiency of the route as much as possible on the premise of ensuring safety, the process of the ATCO commanding the aircraft in the en route phase is analyzed firstly when proposing this optimization model.

    Assuming that there are N aircraft entering the route in turn, ATCO takes over the aircraft from the upstream route one by one, command the aircraft on the principle of first come, first serve (FCFS), and then transfers them to the downstream route at the exit waypoint. The ATCO should not only ensure the safety separation between aircraft, but also guarantee the operation efficiency of the route, that is, reduce the total time of all aircraft occupying the route. In some circumstances, if the longitudinal separation cannot be met, ATCO will issue instructions to aircraft to execute lateral displacement to avoid conflicts, which leads to the increase of the airspace resources usage. Thus, the shorter distance of lateral displacement which aircraft conduct on the route, the better. Furthermore, because aircraft is unevenly distributed on the route, the instantaneous density distribution of the route is also inhomogeneous which raises the safety risk. Therefore, the smaller the variance of longitudinal separation is, the better. In addition, the less instructions issued by ATCO that are different from the initial control scheme, the lower workload.

    In the light of the discussions above, we propose the following mathematical model which considers the total operation time of all aircraft on the route, instantaneous density of the route, maximum lateral displacement, and ATCO's workload.

    minZ=NΔti+Ymax+M+var(Sepai) (1)

    where Δti is the operation time of the ith aircraft on the route, in seconds; Ymax is the maximum lateral displacement of all aircraft, in meters; M is the number of instructions issued by ATCO different from the control scheme which represents the workload of ATCO; var(Sepai) is the variance function of longitudinal separation between the ith aircraft and the preceding aircraft that quantifies the instantaneous density of the route.

    Since the dimensions and magnitudes of the above four indicators are not the same, for the convenience of modeling, normalization is adopted to make the value of them between 0 and 1. The normalization formula used in this section is as follows:

    x=11+x (2)

    Hence, the newly proposed en route capacity optimization model is given by Eq (3).

    minZ=11+NΔti+11+Ymax+11+M+11+var(Sepai) (3)

    In order to be more in line with the actual operation, the proposed route capacity optimization model is implemented by taking the consideration and judgment of the ATCO when issuing the control instructions as the process. Taking into account the main responsibilities of the ATCO in commanding the aircraft, it mainly includes the formulation of control scheme, flight conflict detection, conflict resolution, and the final issue of instructions. Since the process of conflict detection and conflict resolution directly affects the speed, heading and flight altitude of the commanded aircraft, the ATCO instruction modelization can be divided into six parts: speed adjustment instruction, flight level adjustment instruction, lateral displacement instruction, the establishment of control scheme, flight conflict detection and flight conflict resolution.

    In the actual operation process, the speed of the aircraft is difficult to maintain at a constant value which often fluctuates in small amplitudes. And speed changes are more complicated when aircraft in the stage of acceleration or deceleration. To simplify the problem, it is assumed that the aircraft fly at a constant speed or conduct a constant acceleration and deceleration on the route.

    Thus, the acceleration value of the aircraft on the route can be extracted from the aircraft trajectory data which are derived from the ATC automation system. The maximum value of acceleration amax, the average value of acceleration amean, the average value of deceleration bmean and the minimum value of deceleration bmin can be calculated by mathematical statistics. Obviously, the minimum value of acceleration and the maximum value of deceleration are zero.

    Suppose the speed of the aircraft before ATCO issue the speed adjustment instruction is the initial speed v1 and the target speed after the adjustment is v2, then:

    t1={v2v1ameanv1<v2v2v1bmeanv1>v2 (4)
    l1={v22v212ameanv1<v2v22v212bmeanv1>v2 (5)

    where t1 is the time required for speed adjustment; l1 is the distance required for speed adjustment.

    If the calculated time and distance cannot meet the requirement, the acceleration can be calculated inversely. Then judge whether it is within the range of values, otherwise use other instructions.

    In the phase of en route flight, aircraft generally maintains at a fixed altitude, i.e., flight level. When aircraft climbs or descends, the change of speed and altitude is complicated. For the sake of convenience, it is assumed that the aircraft climbs or descends at a constant speed and gradient.

    In the same way, the climb and descent gradient values of the aircraft on the route can also be extracted from the aircraft trajectory data which are derived from the ATC automation system. The maximum value of the climb gradient gradCmax, the average value of the climb gradient gradCmean, the average value of the descent gradient gradDmean and the minimum value of the descent gradient gradDmincan be calculated by mathematical statistics. Obviously, both the minimum value of the climb gradient and the maximum value of the descent gradient are zero.

    Suppose the flight level before adjustment is the initial flight level h1, the target flight level after the adjustment is h2, and the speed of the aircraft during the flight level adjustment is v, then:

    l2={h2h1gradCmeanh1<h2h2h1gradDmeanh1>h2 (6)
    t2=l2v (7)

    where t2 is the time required for flight level adjustment; l2 is the distance required for flight level adjustment.

    In order to alleviate the congestion of downstream route and resolve potential flight conflicts, it is necessary for the ATCO to adjust the separation between aircraft before they converge. In addition to issue speed reduction instructions, the ATCO can also command the aircraft to make a lateral displacement to guarantee the required separation. As shown in Figure 1, there is a route with its entry waypoint and exit waypoint E and F. The green aircraft follows the yellow one. In order to ensure the required separation, the ATCO issue lateral displacement instructions to guide the green aircraft turns left from point M and flies to N, and then turns right at N to go back to the route at K. Meanwhile, the yellow aircraft maintains on the route and flies from point P to Q.

    Figure 1.  Schematic diagram of lateral displacement.

    The constraints of lateral displacement instruction are as follows:

    1) The purpose of the lateral displacement instruction is to delay the following aircraft by a certain distance Ldel or a certain time tdel.

    2) When selecting point N, the length of MN should be greater than two times of actual delayed distance, i.e., 2ΔL and less than 10 km of the pilot's visual range. The initial length of MN can be set to (5+Ldel) km.

    3) The angle of the first turn α1 is generally set to 15°, 30°, 45°or 60°; the angle of the second turn α1+α2 should not be greater than 120°; α1 is generally greater than or equal to α2.

    Taking the delay of a certain distance Ldel as an example, suppose the overall speed of the aircraft is v, and the equation set formed by the constraint of the lateral displacement instruction is:

    {t1=(5+Ldel)/vLMNK=v(t1+t2)LMK=vt1cosα1+vt2cosα2ΔL=LMNKLMKvt2/sinα1=vt1/sinα2α1α2δ=ΔLLdelα1=15304560α1+α2120δLdel×10% (8)

    where t1 and t2 are the flight time of the MN and MK segments respectively, and ΔL is the actual delayed distance. Try different values of α1 and α2 to make sure δ is less than 10% of Ldel.

    For the convenience of calculation, the arc of the turn is not considered in the above formula. However, when the turning angle is less than 60 degrees, the distance error will increase by about 4%, which cannot be ignored. To this end, modify the equations as follows:

    ~ΔL=ΔL+rα12rtan(α1/2)+r(α1+α2)2rtan[(α1+α2)/2]+rα22rtan(α2/2)r=90v2281πtan25 (9)

    Control scheme is always prepared in advance by ATCO before commanding the aircraft and the potential conflicts are not considered. In order to shorten the operation time of each aircraft and increase the efficiency of the route, aircraft that needs to accelerate will be commanded to accelerate as soon as it enters the route, and the one that needs to decelerate will be commanded to decelerate until it leaves the route. In most cases, flight level adjustment instruction is issued in the middle of the route segment, and it may be modified according to the subsequent situation. The process of control scheme is shown in Figure 2.

    Figure 2.  Flowchart of the control scheme.

    Firstly, by establishing the kinematic equation of the aircraft, the slope distance between any two aircrafts or the projection distance in three directions at any time can be calculated by using the distance formula between two points.

    When an aircraft joins the route, its kinematic equation is as follows:

    Xi=xi(t);Yi=yi(t);Zi=zi(t)t(i)0tt(i)L (10)

    The slope distance, longitudinal distance, and vertical distance between this aircraft and other aircraft on the route are:

    {Dslope=(XiXm)2+(YiYm)2+(ZiZm)2Dlen=(XiXm)2Dver=(ZiZm)2 (11)

    According to the criteria of ICAO Doc4444 "Air traffic management" and the national regulation "Civil aviation air traffic management regulations" in China, in radar control circumstances, the en route longitudinal separation of aircraft should not be less than 20 km, and the vertical separation should not be less than 300 m. In order to simplify the calculation, for each moment in the time range, the first step of flight conflict detection is to judge whether the slope distance of two aircraft Dslope is greater than the required minimum value which is 200002+3002 meters. If yes, it means that the required safety separation is met; otherwise, the second step is to judge whether the longitudinal distance between the aircraft and other aircraft Dlen is greater than 20 km or the vertical distance between the aircraft and other aircraft Dver is greater than 300 m, and if yes, the required safety separation is met.

    For the purpose of lower cost, the priority order when the ATCO formulate conflict resolution strategies is as follows:

    1) Issue the speed adjustment instruction which can be satisfied by adjusting the position in the control scheme when the instruction is issued, i.e., Xv.

    2) Issue the flight level adjustment instruction which can be satisfied by adjusting the position in the control scheme when the instruction is issued, i.e., Xh.

    3) Issue lateral displacement instruction

    The lateral displacement in the conflict resolution plan is only to allow the front and rear aircraft to meet the required safe separation, that is, the longitudinal distance between the two aircraft at least 20 km.

    Suppose the speeds of the front and rear aircraft are v2 and v1 respectively, and the new equations are as follows:

    {LPQ=v2tLMNK=v1t1+v1t2t=t1+t2LMK=v1t1cosα1+v1t2cosα2vt2/sinα1=vt1/sinα2α1α2α1=15304560α1+α2120LKQ=LMP+LPQLMKδ2=LKQ2×1040δ22×104×10% (12)

    In addition, when it is necessary to optimize the trajectory and adjust the separation, these strategies of the conflict resolution solutions mentioned above can also be used.

    The uncertain part of the whole model is whether to modify the control scheme and how to modify it, the key parameters of the conflict resolution strategy and the unknown variables in the objective function. Therefore, the decision variables are added as shown in Table 2.

    Table 2.  Values of variables.
    Expression Connotation Value (EMPL = 1) Value (EMPL = 0)
    EMPL 0–1 variable, equal to 1 means the ith aircraft executes the control scheme, equal to 0, the opposite 1 0
    Xhci The position of the ith aircraft receiving the flight level adjustment instruction 50 km [0,L]
    Xvci The position of the ith aircraft receiving the speed adjustment instruction See 2.4 [0,L]
    ai The acceleration of the ith aircraft receiving the speed adjustment instruction See 2.4 [bmin,amax]
    gradi The gradient of the ith aircraft receiving the flight level adjustment instruction See 2.4 [gradDmin,gradCmax]
    Si The conflict resolution strategy used by the ith aircraft, equal to 0 means not to use, equal to 1 means to modify the control scheme, and equal to 2 means to use the lateral displacement strategy. 0 1 or 2

     | Show Table
    DownLoad: CSV

    The solving steps of the model in this paper are as follows in Figure 3.

    Figure 3.  Flowchart of the solving steps of the proposed model.

    This paper takes one route of Xi'an control area in northwest of China as the case study. The route selected is part of air route G212 from waypoints WJC to OKVUM whose length is 128 kilometers and magnetic course is 57 degrees, as shown in Figure 4.

    Figure 4.  Route selected in the case study.

    The aircraft trajectory data used in this paper is the flight information of the aircraft recorded by the control radar in the airspace, with a total of 6 days. The Radar data is derived from the air traffic control automation system and is updated every 4 seconds. After parsing and transcribing, it is saved in a format that is easy to read and manipulate, such as txt format. Figure 5 shows the data header and part of the trajectory data of one aircraft. The data header contains the track number, the aircraft call sign and the number of recorded trajectory points. Trajectory data shows the flight time, longitude, latitude, altitude, airspeed and magnetic heading of the aircraft from left to right.

    Figure 5.  Historical trajectory data of the aircraft.

    In the modeling process, it is necessary to set parameters for the speed adjustment and flight level adjustment instructions issued by ATCO. The actual operational data of the route are extracted by mathematical statistics method, and the average value of the obtained parameters is shown in Table 3.

    Table 3.  Values of the model input parameters.
    Parameters gradCmax gradCmean gradDmean gradDmin amax amean bmean bmin
    Value 0.147 0.019 −0.022 −0.119 89.2 13.9 −12.5 −89.4

     | Show Table
    DownLoad: CSV

    In order to facilitate model simulation, the trajectory data of aircraft is simplified. For example, the speed of the aircraft is rounded up to the next 10 knots increment, and the altitude of the aircraft is rounded according to the flight level. Part of the data is shown in Table 4, where ActTF is the actual departure time of each flight.

    Table 4.  Input data of the model (first 40 flights).
    S.N. TE ActTF hE hF vE vF
    1 0:14:02 0:25:02 11300 10100 200 220
    2 0:14:34 0:32:06 10100 10100 210 200
    3 0:18:46 0:29:34 10100 10100 220 250
    4 0:24:14 0:35:06 10100 10100 220 280
    5 1:51:54 2:2:38 10700 10700 200 290
    6 1:53:14 2:4:02 11600 11300 180 270
    7 2:20:38 2:32:10 8900 8900 190 220
    8 2:27:42 2:38:50 8900 10100 200 350
    9 2:46:54 2:57:54 8900 10700 160 250
    10 3:13:22 3:24:38 10700 10700 200 320
    11 7:11:38 7:23:06 8900 8900 160 350
    12 7:33:54 7:44:50 11300 10100 250 270
    13 7:43:18 7:54:38 8900 8900 200 270
    14 7:46:02 8:02:10 10100 8600 190 140
    15 8:15:06 8:26:10 10700 10100 200 220
    16 8:24:30 8:36:19 7200 8900 230 280
    17 8:44:10 8:56:58 8400 7800 −210 −190
    18 8:44:46 8:56:03 6900 8900 210 200
    19 8:48:58 9:00:10 10100 10100 160 320
    20 8:58:34 9:10:07 7500 10700 250 250

     | Show Table
    DownLoad: CSV

    According to the model proposed in Section 3, the results are obtained in Table 5.

    Table 5.  Output data of the model (first 40 flights).
    S.N. EMPL Xchi Xcvi ai gradi Si Δti TF α1/2
    1 1 50000 0 13.9627 0.0153 0 581.8 00:23
    2 0 127836.6 −12.5463 2 134 00:25 60°/65°
    3 0 0 13.9627 2 188 00:35 60°/65°
    4 0 0 13.9627 2 200 00:37 60°/65°
    5 1 0 13.9627 0 442 01:59
    6 1 50000 0 13.9627 0.0038 0 475 02:01
    7 1 0 13.9627 0 582 02:30
    8 1 50000 0 13.9627 −0.0154 0 368 02: 34
    9 1 50000 0 13.9627 0.0199 0 513 02:55
    10 1 0 13.9627 0 402 03:20
    11 1 0 13.9627 0 369 07:18
    12 1 50000 0 13.9627 0.0154 0 474 07:42
    13 1 0 13.9627 0 475 07:51
    14 1 50000 127342.4 −12.5463 −0.0218 0 674 07:57
    15 1 50000 0 13.9627 0.0077 0 582 08:25
    16 1 50000 0 13.9627 0.0199 0 457 08:32
    17 1 50000 0 13.9627 0.0077 0 674 08:33
    18 1 50000 127836.6 −12.5463 −0.0256 0 609 08:55
    19 0 0 13.9627 1 423 08:56
    20 1 50000 0.0410 0 512 09:07

     | Show Table
    DownLoad: CSV

    Statistical analysis of the input data according to the aircraft trajectory and the simulation results of the optimization model is shown in Table 6. According to the comparison in this table, it can be seen that the total operation time is reduced by 18.6%; the maximum lateral displacement distance is magnificently reduced by 57.4%; the variance of the separation of aircraft are also reduced by 20.3%; and the en route capacity is improved by 23%.

    Table 6.  Comparison of simulation results and actual trajectory data.
    Average operation time
    (Seconds)
    Maximum lateral displacement distance (Meters) var(Sepa) (Meters) Capacity
    (Flights/hour)
    Actual trajectory 691.65 36790.8 1030 20.8
    Simulation results 562.78 15669.6 820 25.6

     | Show Table
    DownLoad: CSV

    Figures 69 demonstrate the plane view and profile view of the aircraft trajectory on this route before and after optimization. It can be seen intuitively from the above figures that the en route capacity optimization model established according to the ATCO's command process can effectively reduce the en route lateral displacement distance and the amount of flight level changes of aircraft, so as to improve the route utilization efficiency and reduce the workload of ATCO and pilots.

    Figure 6.  Plane view of the historical trajectory.
    Figure 7.  Plane view of the simulated trajectory.
    Figure 8.  Profile view of the historical trajectory.
    Figure 9.  Profile view of the simulated trajectory.

    The improvement of the en route capacity problem has increasing drawn the attention of the research community, policy makers and professionals in the ATC community. Existing research has placed the main focus on the capacity evaluation and route structural optimization, but there is still ample room for further investigation of the en route capacity optimization model. This paper presents an en route capacity optimization model to improve the route capacity and decrease the workload of ATCO. Compared with the available approaches, the proposed model envisages the following two innovations: (ⅰ) the whole process of ATCO command is considered in the model, especially the control scheme and different types of instructions, and (ⅱ) the en route historical trajectory data of aircraft is used to as the key parameters of the input data to efficiently yield the acceptable results of the model. A case study is conducted to validate the applicability and feasibility of the model. Results show that the total en route operation time, maximum lateral displacement distance and variance of the separation of aircraft will be considerably reduced, while the capacity and the utilization efficiency will be increased as compared to the conventional model. Although, the experimental results are encouraging in terms of the ability of the proposed capacity optimization model to tackle en route problem instances, further experimentation is required to: (ⅰ) ensure that the proposed model performs well on the route network instances that resemble more closely to real-world cases, and (ⅱ) generalize the results regarding the implications of aircraft classifications.

    The work reported in this paper has been jointly funded by the National Natural Science Foundation of China and Civil Aviation Administration of China through Programme Grant U1633124 "Research on Traffic Capacity of Airway Intersections Based on TAAM Simulation Platform and Traffic Flow Feature Recognition".

    The authors declare there is no conflict of interest.



    [1] ADR in Beijing: Civil aircraft market forecast annual report 2020-2039, 2020. Available from: https://www.avic.com/c/2020-12-22/513568.shtml.
    [2] Y. Lin, L. Li, P. Ren, Y. Wang, W. Szeto, From aircraft tracking data to network delay model: A data-driven approach considering en-route congestion, Transp. Res. Part C, 131 (2021), 103329. https://doi.org/10.1016/j.trc.2021.103329 doi: 10.1016/j.trc.2021.103329
    [3] R. A. Paielli, H. Erzberger, Trajectory specification for terminal air traffic: pairwise conflict detection and resolution, in AIAA Aviat. Technol. Integr. Oper. Conf, 2017,312-318. https://doi.org/10.2514/6.2017-4256
    [4] N. Shen, H. Idris, Tradeoff between airspace capacity and risk mitigation, in AIAA Guid. Navig. Control Conf., 2013,512-526. https://doi.org/10.2514/6.2013-5031
    [5] J. Kuchar, L. Yang, A review of conflict detection and resolution modeling methods, IEEE Trans. Intell. Transp. Syst., 1 (2000), 179-189. https://doi.org/10.1109/6979.898217 doi: 10.1109/6979.898217
    [6] Y. Wang, The analysis and design of emergency ATC system, Master thesis, Sichuan University, 2005.
    [7] A. Lecchini, W. Glover, J. Lygeros, J. Maciejowski, Monte Carlo optimization strategies for Air Traffic Control, in AIAA Guid. Navig. Control Conf., 2005. https://doi.org/10.2514/6.2005-5821
    [8] I. Hwang, J. Hwang, C. Tomlin, Flight-mode-based aircraft conflict detection using a residual-mean interacting multiple model algorithm, in AIAA Guid. Navig. Control Conf., 2003,312-318. https://doi.org/10.2514/6.2003-5340
    [9] W. Liu, I. Hwang, Probabilistic trajectory prediction and conflict detection for air traffic control, J. Guid. Control Dyn., 34 (2011), 1183-1188. https://doi.org/10.2514/1.53645 doi: 10.2514/1.53645
    [10] M. Porretta, M. Dupuy, W. Schuster, A. Majumdar, W. Ochieng, Performance Evaluation of a Novel 4D Trajectory Prediction Model for Civil Aircraft, J. Navig., 61, (2008), 393-420. https://doi.org/10.1017/S0373463308004761 doi: 10.1017/S0373463308004761
    [11] J. Benavides, J. Kaneshige, S. Sharma, R. Panda, M. Steglinski, Implementation of a trajectory prediction function for trajectory based operations, in AIAA. Atmos. Flight Mech. Conf., 2014,189-203. https://doi.org/10.2514/6.2014-2198
    [12] X. Tang, J. Gu, Z. Shen, P. Chen, A flight profile clustering method combining TWED with K-means algorithm for 4D trajectory prediction, in 2015 Integr. Commun. Navig. Surv. Conf., 2015,489-501. https://doi.org/10.1109/ICNSURV.2015.7121260
    [13] K. Tastambekov, S. Puechmorel, D. Delahaye, C. Rabut, Aircraft trajectory forecasting using local functional regression in sobolev space, Transp. Res. Part C, 39 (2014), 1-22. https://doi.org/10.1016/j.trc.2013.11.013 doi: 10.1016/j.trc.2013.11.013
    [14] Z. Shi, M. Xu, Q. Pan, B. Yan, H. Zhang, LSTM-based flight trajectory prediction, in 2018 Int. Joint Conf. Neural Network (IJCNN), 2018, 1156-1162. https://doi.org/10.1109/IJCNN.2018.8489734
    [15] Z. Zhang, R. Yang, Y. Fang, LSTM network based on antlion optimization and its application in flight trajectory prediction, IEEE Adv. Inf. Manage. Comun. Electron. Auto. Control Conf, 2018,356-465. https://doi.org/10.1109/IMCEC.2018.8469476 doi: 10.1109/IMCEC.2018.8469476
    [16] Y. Cui, X. Wei, H. You, Adaptive forecast model for uncertain track, ACTA. Aerosp. Astro. Sinica., 40 (2019), 241-250.
    [17] R. A. Paielli, H. Erzberger, Conflict probability estimation for free flight, J. Guid. Control Dyn., 20 (1997), 588-596. https://doi.org/10.2514/2.4081 doi: 10.2514/2.4081
    [18] J. Krozel, M. Peters, Strategic conflict detection and resolution for free flight, in IEEE Conf. Decis. Control, 1997, 45-59. https://doi.org/10.1109/CDC.1997.657844
    [19] M. Prandini, J. Hu, J. Lygeros, S. Sastry, A probabilistic approach to aircraft conflict detection, IEEE Trans. Intell. Transp. Sys., 1 (2000), 199-220. https://doi.org/10.1109/6979.898224 doi: 10.1109/6979.898224
    [20] D. Jacquemart, J. Morio, Adaptive interacting particle system algorithm for aircraft conflict probability estimation, Aerosp. Sci. Technol., 55 (2016), 431-438. https://doi.org/10.1016/j.ast.2016.05.027 doi: 10.1016/j.ast.2016.05.027
    [21] L. Shi, Research on probabilistic flight conflict detection algorithms in air traffic management, Master thesis, Tianjin University, 2014.
    [22] I. Hwang, C. E. Seah, Intent-based probabilistic conflict detection for the next generation air transportation system, Proc. IEEE, 96 (2008), 2040-2059. https://doi.org/10.1109/JPROC.2008.2006138 doi: 10.1109/JPROC.2008.2006138
    [23] R. A. Paielli, H. Erzberger, D. Chiu, K. R. Heere, Tactical conflict alerting aid for air traffic controllers, J. Guid. Control Dyn., 32 (2009), 184-193. https://doi.org/10.2514/1.36449 doi: 10.2514/1.36449
    [24] H. A. P. Blom, G. J. Bakker, Conflict probability and in-crossing probability in air traffic management, in IEEE Conf. Decis. Control, 2002, 2421-2426.
    [25] G. J. Bakker, H. J. Kremer, H. A. P. Blom, Geometric and probabilistic approaches towards conflict prediction, USA/Eur. ATM R & D Semi., 2000, 1-10.
    [26] H. Erzberger, R. A. Paielli, D. R. Isaacson, M. M. Eshow, Conflict detection and resolution in the presence of prediction error, USA/Eur. ATM R & D Semi., 1997, 17-20.
    [27] W. Liu, C. E. Seah, I. Hwang, Aircraft 4D trajectory prediction and conflict detection for air traffic control, IEEE Conf. Decis. Control, 3 (2010), 1-7.
    [28] S. Cafieri, R. Omheni, Mixed-integer nonlinear programming for aircraft conflict avoidance by sequentially applying velocity and heading angle changes, Eur. J. Oper. Res., 260 (2017), 283-290. https://doi.org/10.1016/j.ejor.2016.12.010 doi: 10.1016/j.ejor.2016.12.010
    [29] A. A. Alonso-Ayuso, L. F. Escudero, F. J. Martín-Campo, A mixed 0-1 nonlinear optimization model and algorithmic approach for the collision avoidance in ATM: velocity changes through a time horizon, Comput. Oper. Res., 39 (2012), 3136-3146. https://doi.org/10.1016/j.cor.2012.03.015 doi: 10.1016/j.cor.2012.03.015
    [30] J. Omer, A space-discretized mixed-integer linear model for air conflict resolution with speed and heading maneuvers, Comput. Oper. Res., 58 (2015), 75-86. https://doi.org/10.1016/j.cor.2014.12.012 doi: 10.1016/j.cor.2014.12.012
    [31] Q. Zhou, R. Zhang, X. Suo, K. Qiang, UAV genetic algorithm trajectory planning with time constraints, Aerosp. Comput. Technol., 46 (2011), 93-97.
    [32] Z. Hu, Research on some key techniques of UAV path planning based on intelligent optimization algorithm, Ph.D thesis, Nanjing University of Aeronautics and Astronautics, 2011.
    [33] Y. Chiang, J. T. Klosowski, C. Lee, J. S. B. Mitchell, Geometric algorithms for conflict detection/resolution in air traffic management, in IEEE Conf. Decis. Control, 1997, 1-12. https://doi.org/10.1109/CDC.1997.657848
    [34] C. Tomlin, G. J. Pappas, S. Sastry, Conflict resolution for air traffic management: A study in multi agent hybrid systems, IEEE Trans. Autom. Control, 43 (1998), 509-521. https://doi.org/10.1109/9.664154 doi: 10.1109/9.664154
    [35] A. Bicchi, L. Pallottino, On optimal cooperative conflict resolution for air traffic management systems, IEEE Trans. Intell. Transp Syst., 1 (2000), 221-231. https://doi.org/10.1109/6979.898228 doi: 10.1109/6979.898228
    [36] C. Goodchild, M. A. Vilaplana, S. Elefante, Co-operative optimal airborne separation assurance in free flight airspace, USA/Eur. ATM R & D Semin., 2000.
    [37] A. M. Bayen, P. Grieder, G. Meyer, C. J. Tomlin, Lagrangian delay predictive model for sector-based air traffic flow, J. Guid. Control Dyn., 28 (2005), 1015-1026. https://doi.org/10.2514/1.15242 doi: 10.2514/1.15242
    [38] N. Durand, J. M. Alliot, F. Edioni, Neural nets trained by genetic algorithms for collision avoidance, Appl. Intell., 13 (2000), 205-213. https://doi.org/10.1023/A:1026507809196 doi: 10.1023/A:1026507809196
  • This article has been cited by:

    1. A V Mikhailov, S S Kumkov, Air-traffic controller and autopilot joint control of aircraft as mixed integer quadratic problem, 2023, 2514, 1742-6588, 012013, 10.1088/1742-6596/2514/1/012013
    2. Honghao Chen, Xinping Zhu, Yajun Zou, Zhongkun Li, Tianxiong Zhang, Indrajit Ghosh, Optimized ST-Moran’s I Model for Characterizing the Dynamic Evolution of Terminal Airspace Congestion, 2024, 2024, 2042-3195, 1, 10.1155/2024/7151746
    3. Verma Aditya, Dande Sureshkumar Aswin, Somasundaram Vanitha Dhaneesh, Sakthivelan Chakravarthy, Bhukya Shanmuk Kumar, Marimuthu Venkadavarahan, A review on air traffic flow management optimization: trends, challenges, and future directions, 2024, 5, 2662-9984, 10.1007/s43621-024-00781-7
    4. Rafik Amara, Maamar Bougherara, Brahim Gouraya, Mostefa Belhadj Aissa, 2024, Air Traffic Optimization in Training Environment: A Comparison Between Genetic Algorithms and Simulated Annealing, 979-8-3503-6756-0, 1, 10.1109/ICCA62237.2024.10927963
  • Reader Comments
  • © 2022 the Author(s), licensee AIMS Press. This is an open access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0)
通讯作者: 陈斌, bchen63@163.com
  • 1. 

    沈阳化工大学材料科学与工程学院 沈阳 110142

  1. 本站搜索
  2. 百度学术搜索
  3. 万方数据库搜索
  4. CNKI搜索

Metrics

Article views(2264) PDF downloads(119) Cited by(4)

Figures and Tables

Figures(9)  /  Tables(6)

Other Articles By Authors

/

DownLoad:  Full-Size Img  PowerPoint
Return
Return

Catalog