# An Improved Runway Operation Capacity Model for V-Open Multirunway Airports in China

Authors: Zhiyuan Shen; Xinyu Xu; Yong He; YongGang Yan; Lin Zhou; Yingying Hu
Journal: Journal of Advanced Transportation
 (2022)

---

## Abstract

Diverging runways are receiving considerable attention as replacements for parallel runways when constructing three or more runways in China. This has been prompted by the limitations associated with a single-direction airspace and associated regulations when using parallel runways. However, the operating efficiency of diverging runways is still greatly constrained by the separation standard for departing aircraft in CCAR-93TM-R5 and Decree of CAAC NO.123. To address these problems, an improved runway operation capacity model for a V-open multirunway airport is proposed to reduce the separation standard based on the equivalent lateral separation operation (ELSO) standard. A collision risk module is first proposed to assess the feasibility and safety of the proposed policy. Also, a novel runway capacity model for V-open multiple runways is then built to improve the operating efficiency for diverging runway. A series of comparative experiments based on actual flight data from Chengdu Tianfu International Airport demonstrated that reducing the divergence angle from 15° to 10° improved the runway operation capacity by about 8% based on the proposed improved policy. These research results could play a major role in facilitating appropriate amendments to the associated CAAC operation regulations.


---

## Body

## 1. Introduction

The International Air Transport Association (IATA) has predicted that China will become the largest civil aviation market in the world by around 2025 and involve the flow of 1.6 billion passengers by around 2037 [1]. Then, there is no doubt that COVID-19 significantly impacted civil aviation market and considered as the main factor in much recent studies [2]. Research studies on runway capacity estimation usually focus on the maximum of capacity in peak hour or one day. So, airline and airport pay more attention to theoretical maximum flight number under the priority of operation safety. By reducing the impact of the diverging angle, the departure separation can be decreased. Thereby, the numbers of airport flights permitted to taking-off and landing at peak hour increase. The ticket price at peak hour is usually higher than those at other time. Consequently, it leads to a higher airline revenue.In the past, large airports at China tended to build additional parallel runways to improve their operation capacity. However, it has been difficult to reach estimated maximum operation capacities due to limitations associated with a single-direction airspace and its associated regulations. The diverging runway has advantages of low building cost and flexible operation procedures; therefore, it is the preferred choice for configuration of large domestic hub airports. It has led to diverging runways being considered for replacing parallel runways when constructing three or more runways in China, such as at Beijing Daxing International Airport (ZBAD) and Chengdu Tianfu International Airport (ZUTF). While configuring runways in different directions facilitates the flexible utilization of airport airspace and operation procedures, such diverging runways have mainly been used for aircraft takeoffs [3].However, the departure efficiency for diverging runways is still greatly affected by the separation standard for departing aircraft in CCAR-93TM-R5 [4] and Decree of CAAC NO.123 [5]. Specifically, aircraft departing simultaneously from parallel or diverging runways need to comply with the separation standard that their takeoff courses diverge by 15° or more immediately after departure. The initial purpose of implementing diverging departures is to improve the operating efficiency and reduce flight delays [6]. However, numerous studies have shown that the strict separation standard greatly restricts the use of the available airspace, especially for multiple runways and diverging runways [7, 8]. For example, a larger diverging airspace is needed when applying the separation standard to three parallel runways, and separated departure flight paths could overlap the diverging runways to result in increased waiting time for aircraft departures or reducing the takeoff runway distance, resulting in a decreased departure capacity. Therefore, many studies have attempted to reduce the divergence angle while still conforming to safety standards [9, 10].The recent rapid development and application of performance based on navigation (PBN) has significantly improved navigation precision in the terminal area. A novel concept of the equivalent lateral separation operation (ELSO) standard [11] was proposed, which is able to reduce divergence angles while maintaining established minimum lateral spacing between departure paths. A series of related theorical and practical evaluation experiments were conducted at Atlanta Airport over a 1-year period in 2011 [12]. The Federal Aviation Administration (FAA) finally approved for the divergence angle to be reduced to 10° at Atlanta Airport based on the safety of the experimental results, which facilitated the development of a new departure procedure. Nowadays, any airport in the USA can operate with the standard 10° divergence angle when both aircraft are flying RNAV (area navigation) standard instrument departures (SID) from parallel runways according to the policy in FAAO.7110.65W [13].All of the international airports in China have been designed and operated under RNAV SID procedures, and more runways are planned to be built in order to meet the increasing number of flights at the busiest Chinese airports (e.g., Shanghai Pudong, Hangzhou Xiaoshan, and Yunnan Changshui), including to stimulate the economy affected by COVID-19 with greater infrastructure investment [14–16]. Reducing delays and improving operation efficiency are two important steps to reduce the impact of COVID-19 on the aviation industry. However, the superiority of the higher navigation precision of RNAV has not been performed to improve the operating efficiency for the applicable standard and additional runways [17, 18]. Most importantly, the existing separation standard, including the minimum requirement for a 15° divergence angle between independent parallel departure operations, significantly constrains the potential of multirunway operations [19].In accordance with the development strategies of the CAAC and the Central Military Commission (CMC) of ATC that pertain to higher flight punctuality rates and better civil airspace planning [20], based on the actual aircraft departure path data from ZUTF parallel runways, we proposed modifying the regulations associated with the separation standard by referring to the ELSO standard. We further proposed a capacity estimation model for V-open multiple runways and tested the operating efficiency of the proposed new separation standard relative to the conventional one. Specifically, taking ZUTF as an example, two main modules are developed to test the assumption: (1) the collision risk module based on ELSO for assessing the feasibility and safety of the proposed policy, and (2) the V-open multirunway capacity module for measuring the increases in runway capacity. This study provides theoretical support and reference data for the development of future regulations to reduce the divergence angle standard and the operation procedure of V-open runways. All terminology used in this paper is listed in Table 1.Table 1
List of terminology used in this paper.
TerminologyReferred toELSOEquivalent lateral separation operationCAACCivil Aviation Administration of ChinaIATAInternational Air Transport AssociationZBADBeijing Daxing International AirportZUTFChengdu Tianfu International AirportPBNPerformance based on navigationFAAFederal Aviation AdministrationRNAVArea navigationSIDStandard instrument departuresCMCCentral Military CommissionThe remainder of this paper is organized as follows: Section2 describes the basic theory of V-open runway configuration and ELSO, Section 3 details the proposed collision risk module and V-open multirunway capacity module, Section 4 reports on the theoretical analyses and experimental simulations performed to test our proposed framework, and Section 5 summarizes the obtained findings.


## 2. Basic Theory

The current divergence requirements for departing aircraft are defined in CCAR-93TM-R5 and Decree of CAAC NO.123. There are three key rules involved in separating aircraft departing from the same runway or parallel runways based on using either conventional SID or RNAV SID. In all cases, aircraft must establish radar identification within 1 mile of the end of the takeoff runway and travel along courses that diverge by 15° or more immediately after departure.


### 2.1. Layout of V-Open Runways

The layout of the runways at two airports (ZBAD and ZUTF) with diverging runways served as the research model of V-open runways in this study. As illustrated in Figure1, including three runways totally in our study scene, two runways are parallel, and the other runway is V-open with the parallel runways. For parallel runways, we try to use ELSO to reduce the angle of course divergence for improving the efficiency of airspace usage. For the V-open runway, the main effort is to reduce the coverage impact of the parallel runway departure courses on the third runway. The diverging runway in the figure is used only for aircraft departures toward the east direction. One departure course coincides with the runway centerline extension; another one diverged by alpha angle. In this case, pilots are rather easier to master departing course whatever runway is used. On the other hand, the location of the runway direction is significantly affected by the surrounding terrain, so departure courses along the extension line of the runway centerline could minimize obstacle restrictions. Therefore, there is no need to simultaneously consider mutual diversion and considering one runway’s trajectory divergence could better estimate the runway.Figure 1
Layout of a V-open runway.In accordance with the regulations in CCAR-93TM-R5, aircraft can depart from parallel runways simultaneously when the centerlines of the parallel runways are separated by at least 760 m and the courses of the aircraft diverge by 15° or more immediately after departure. In terms of approach operations, correlated parallel approaches and independent parallel approaches are authorized for minimum parallel runway spacing of 915 m and 1035 m, respectively.The departure course relative to the runway centerline might cover the takeoff zone of the diverging runway in Figure1. In other words, the departure interval for diverging runways is highly dependent on the time separation between aircraft flying over the intersection zone. There is no doubt that this situation will adversely affect the airport operating efficiency and runway capacity [21]. The concept of ELSO is introduced to solve this problem in this study.


### 2.2. Principle of ELSO

ELSO was proposed to meet the minimum spacing of the general divergence standard defined in CCAR-93TM-R5. In other words, the divergence angle can be reduced as long as the established minimum lateral spacing between departure paths is maintained. After a 1-year safety evaluation experiment was conducted at Atlanta Airport in 2011, the FAA authorized that the course divergence angle could be reduced to 10° when both aircraft are departing from parallel runways and are flying RNAV SID.For the independent parallel departure operation, the general standard defines that the minimum spacing between parallel runwaysrmin should be 760 m and that the departure-course divergence αmin should be a minimum of 15°. Figure 2 illustrates some key parameters of the general departure operation. In the figure, dashed lines represent the nominal departure path and the distance (d) from the runway end will determine nominal spacing (n(d)) between the nominal departure path, which can be obtained as follows: [22](1)nd=rmin+d×tanαmin.Figure 2
Conventional departure path and associated parameters.The actual path deviation due to various factors such as the wind direction is also illustrated in Figure2. We can assume that the resulting width of the departure path can be estimated by single or multiple standard deviations σ of the related dispersion distribution. Hence, σSTO,C represents the track width of the conventional straight departure procedure along the centerline of the runway, and σDIV,C indicates the track width of the diverging conventional departure procedure. The lateral spacing of the diverging departure path is given as follows:(2)sd=rmin+d×tanαmin−σDIV,C−d×tanσSTO,C.Therefore, distancesd can be regarded as the minimum lateral spacing between actual departure flight paths under the conventional divergence standard and serves as the reference distance when using the ELSO standard.When parallel runways are separated by more than 760 m, lateral spacing requirementsd can also be achieved for a smaller divergence angle, as shown in Figure 3. In the figure, σSTO,R indicates the width of the RNAV departure path following the extended centerline, and the width of the RNAV departure path whose centerline is angled relative to the runway centerline is σDIV,R.Figure 3
Improved departure path based on ELSO.Given reference distancesd and staggered distance t at the end of the runway, the runway distance can be expressed as follows:(3)r=sd+d·tanσSTO,Rd−d+ttanβ−σDIV,Rd+t,where σSTO,Rd and σDIV,Rd+t represent the actual track divergence characteristics of RNAV straight-line departure and offset departure, respectively. Therefore, based on the specific runway configuration, equivalent divergence angle β that can implement independent parallel departures under the ELSO standard is as follows:(4)β=arctan1d+t·sd−r+d·tanσSTO,Rd+σDIV,Rd+t.Equivalent divergence angleβ can ensure that the departure path distance at the end of the runway d is the same as or even greater than the reference distance. If the nominal departure path fails to establish a lateral separation of d nautical miles from the end of the runway, the ELSO divergence standard will assume that the departure will be implemented with the traditional 15° divergence angle after this distance.


## 2.1. Layout of V-Open Runways

The layout of the runways at two airports (ZBAD and ZUTF) with diverging runways served as the research model of V-open runways in this study. As illustrated in Figure1, including three runways totally in our study scene, two runways are parallel, and the other runway is V-open with the parallel runways. For parallel runways, we try to use ELSO to reduce the angle of course divergence for improving the efficiency of airspace usage. For the V-open runway, the main effort is to reduce the coverage impact of the parallel runway departure courses on the third runway. The diverging runway in the figure is used only for aircraft departures toward the east direction. One departure course coincides with the runway centerline extension; another one diverged by alpha angle. In this case, pilots are rather easier to master departing course whatever runway is used. On the other hand, the location of the runway direction is significantly affected by the surrounding terrain, so departure courses along the extension line of the runway centerline could minimize obstacle restrictions. Therefore, there is no need to simultaneously consider mutual diversion and considering one runway’s trajectory divergence could better estimate the runway.Figure 1
Layout of a V-open runway.In accordance with the regulations in CCAR-93TM-R5, aircraft can depart from parallel runways simultaneously when the centerlines of the parallel runways are separated by at least 760 m and the courses of the aircraft diverge by 15° or more immediately after departure. In terms of approach operations, correlated parallel approaches and independent parallel approaches are authorized for minimum parallel runway spacing of 915 m and 1035 m, respectively.The departure course relative to the runway centerline might cover the takeoff zone of the diverging runway in Figure1. In other words, the departure interval for diverging runways is highly dependent on the time separation between aircraft flying over the intersection zone. There is no doubt that this situation will adversely affect the airport operating efficiency and runway capacity [21]. The concept of ELSO is introduced to solve this problem in this study.


## 2.2. Principle of ELSO

ELSO was proposed to meet the minimum spacing of the general divergence standard defined in CCAR-93TM-R5. In other words, the divergence angle can be reduced as long as the established minimum lateral spacing between departure paths is maintained. After a 1-year safety evaluation experiment was conducted at Atlanta Airport in 2011, the FAA authorized that the course divergence angle could be reduced to 10° when both aircraft are departing from parallel runways and are flying RNAV SID.For the independent parallel departure operation, the general standard defines that the minimum spacing between parallel runwaysrmin should be 760 m and that the departure-course divergence αmin should be a minimum of 15°. Figure 2 illustrates some key parameters of the general departure operation. In the figure, dashed lines represent the nominal departure path and the distance (d) from the runway end will determine nominal spacing (n(d)) between the nominal departure path, which can be obtained as follows: [22](1)nd=rmin+d×tanαmin.Figure 2
Conventional departure path and associated parameters.The actual path deviation due to various factors such as the wind direction is also illustrated in Figure2. We can assume that the resulting width of the departure path can be estimated by single or multiple standard deviations σ of the related dispersion distribution. Hence, σSTO,C represents the track width of the conventional straight departure procedure along the centerline of the runway, and σDIV,C indicates the track width of the diverging conventional departure procedure. The lateral spacing of the diverging departure path is given as follows:(2)sd=rmin+d×tanαmin−σDIV,C−d×tanσSTO,C.Therefore, distancesd can be regarded as the minimum lateral spacing between actual departure flight paths under the conventional divergence standard and serves as the reference distance when using the ELSO standard.When parallel runways are separated by more than 760 m, lateral spacing requirementsd can also be achieved for a smaller divergence angle, as shown in Figure 3. In the figure, σSTO,R indicates the width of the RNAV departure path following the extended centerline, and the width of the RNAV departure path whose centerline is angled relative to the runway centerline is σDIV,R.Figure 3
Improved departure path based on ELSO.Given reference distancesd and staggered distance t at the end of the runway, the runway distance can be expressed as follows:(3)r=sd+d·tanσSTO,Rd−d+ttanβ−σDIV,Rd+t,where σSTO,Rd and σDIV,Rd+t represent the actual track divergence characteristics of RNAV straight-line departure and offset departure, respectively. Therefore, based on the specific runway configuration, equivalent divergence angle β that can implement independent parallel departures under the ELSO standard is as follows:(4)β=arctan1d+t·sd−r+d·tanσSTO,Rd+σDIV,Rd+t.Equivalent divergence angleβ can ensure that the departure path distance at the end of the runway d is the same as or even greater than the reference distance. If the nominal departure path fails to establish a lateral separation of d nautical miles from the end of the runway, the ELSO divergence standard will assume that the departure will be implemented with the traditional 15° divergence angle after this distance.


## 3. Proposed Modules

A collision risk model was proposed to quantitatively evaluate safety when the divergence angle is reduced based on ELSO. The minimum divergence angle can also be obtained under the safety requirement. In order to further assess the benefits of the new standard, we determined the operating efficiency of a V-open multirunway capacity model.


### 3.1. Collision Risk Model

The collision risk model mainly considered the path deviation of the actual flight caused by the precision of navigation equipment, which is greatly affected by the aircraft speed. We can assume that both the navigation error and speed error conform to a Gaussian distribution and are given as follows:(5)εny∼Nμn,σn2,(6)εay∼N12μat2,14σa2t4,where μn and μa indicate the mean values and σn and σa are the standard deviations. Subscripts ny and ay denote the aircraft lateral deviation caused by two factors. The lateral collision risk can be computed using equations (7) to (10):(7)εy∼Nμy,σy2μy=μn+12μat2,σy2=σn2+14σa2t4,where εy denotes the distribution of the lateral deviation. The actual lateral distance between simultaneously departing aircraft is given as follows:(8)Dy′t=L2y′t−L1y′t=Dyt+ε2y−ε1y∼NDyt+μ2y−μ1y,σ2y2+σ1y2L1y′t=L1yt+ε1y,L2y′t=L2yt+ε2y,where L1y′t and L2y′t represent the actual path positions of two different aircraft. The probability density function for the lateral direction can be expressed as follows:(9)fy=12πσ2y2+σ1y2exp−y−Dyt+μ2y−μ1y22σ2y2+σ1y2.The probability density function of the actual distance between the two aircraft is integrated to obtain the lateral collision risk:(10)Py=py1<Dy′t<y2=∫y1y2fydy,where y1 and y2 represent the minimum and maximum separations between the aircraft, respectively, and can be calculated from the aircraft wing spans. Similarly, the probability density functions in the vertical and horizontal directions are given as follows:(11)fx=12πσ2x2+σ1x2exp−x−Dxt+μ2x−μ1x22σ2x2+σ1x2,(12)fz=12πσ2z2+σ1z2exp−z−Dzt+μ2z−μ1z22σ2z2+σ1z2,Finally, the collision risk model is obtained as follows:(13)P=∫x1x2fxdx∫y1y2fydy∫z1z2fzdz.


### 3.2. V-Open Multirunway Capacity Model

Wake turbulence was one of the main factors influencing the design of the V-open multirunway capacity model [23, 24]. At domestic hub airports, the number of medium aircraft accounts for 99% and heavy aircraft accounts for 1% at Chinese domestic hub airports. Therefore, in most cases, it is necessary to consider the wake impact for medium aircraft following heavy or medium aircraft, and detailed wake interval and accounts of different categories are shown in Table 2. Mainly, the wake generated from taking-off aircraft on the diverging runway affects the operation of departing aircraft on parallel runways and vice versa. A typical layout of V-open multiple runways is illustrated in Figure 1, in which there are two parallel runways and one diverging runway. There are generally two approach operations for parallel runways depending on whether or not they have different widths: independent parallel approach and correlated parallel approach [25]. Small airports in China currently use independent parallel approaches due to the high requirements for equipment and working skills. However, another diverging path of 30° for aircraft to safely go around also needs to be considered [26]. Therefore, the capacity evaluation model was proposed here based only on the correlated parallel approach.Table 2
Wake influence parameters.
CategorySH-SHSH-SMSM-SHSM-SMWake interval8 km6 km10 km6 kmAccounts4%16%16%64%Remarks: SH/SM denotes heavy and medium aircraft separately.An operation schematic diagram of parallel runways using the correlated parallel-approach-independent parallel departure mode is shown in Figure4. Both runways are used for aircraft arrivals and departures. Aircraft arrive and depart on each runway independently of each other. However, in order to ensure the safe operation of aircraft, the aircraft must land in the same direction on Runways 1 and 2. A safe distance needs to be maintained between aircraft on the parallel runway approaches δa.Figure 4
Parallel approach operations.Assuming that the distance between the two runways isD, as shown in Figure 4, the distance between two successively approaching aircraft on Runway 1 is the following:(14)lmin=2δa2−D2.Then, the time interval between two successively approaching aircraft a and aircraftc on Runway 1 is as follows:(15)tmin=maxlminv.The controller can use the insertion method to arrange for aircraftb to land on Runway 2 between aircraft a and aircraft c, but the time interval between aircraft a and aircraft c must be consistent with the following formula [27]:(16)Tac+Bac≥tmin.The operation sequence diagram of three runways for Runways 1 and 2 using related parallel-approach-independent parallel departure procedures is shown in Figure5. Assuming that approaching aircraft a on Runway 1 arrives at the runway threshold at 00:00, then within the interval between aircraft a and aircraft c, the departing aircraft can be released on Runway 1 or the approaching aircraft can be arranged to land on Runway 2. During time gap G2 of the continuous-approach aircraft on Runway 2, the departing aircraft can continue to be released, but the controller can also release the departing aircraft on Runway 3 after the departing aircraft on Runway 2 has been released for time tz.Figure 5
Time sequence model of V-open three-runway operation.Time intervalG1 between two consecutively approaching aircraft on Runway 1 is as follows:(17)G1=T1−2−T1−1=Tac+Bac−ROTa−maxtbc,tmin.WhenG1>0 (i.e., when the following conditions are satisfied), n departing aircraft can be inserted in gap G1:(18)ETac+Bac≥EROTa+n−1⋅Eεac+Emaxtbc,tmin+Eτ.Then, the entry and exit capacities of Runway 1 are as follows:(19)C1−arrival=3600ETac+Bac,(20)C1−departure=∑i=1n∑j=1nTG1⋅Pij⋅nij,where TG1 is the total number of clearances on Runway 1 within 1 h and nij is the number of departing aircraft that can be released between two adjacent approaching aircraft N/M on Runway 1.Time intervalG2 between two successively approaching aircraft on Runway 2 is as follows:(21)G2=T1−2−T1−1=Tij+Bij−maxROTb,th−maxt2,δvj.When the following conditions are met,k departing aircraft can be inserted in gap G2:(22)ETij+Bij≥EmaxROTb,th+n−1⋅Emaxtz+th,εij+Emaxt2,δvj+Eτ.Then, the entry and exit capacities of Runway 2 are as follows:(23)C2−arrival=∑i=1n∑j=1n∑xTG1GxETij+Bij,(24)C2−departure=∑i=1n∑j=1nTG2⋅Pij⋅kij,where Gx is the time interval between the first successive approach aircraft on Runway 1, TG1 and TG2 are the total numbers of clearances on Runways 1 and 2, respectively, within 1 h, and kij is the number of takeoff and departing aircraft that can be inserted in the time interval between two consecutively arriving aircraft i and j on Runway 2.In gapG2 of Runway 2, assuming that Runway 3 can release m aircraft, according to the operation sequence diagram, the following equation can be derived:(25)Tj−t2−T2−1≥m−1⋅maxTDd,tz+th,εij+τ,Therefore, the departure capacity of Runway 3 is as follows:(26)C3−departure=∑i=1n∑j=1nTG2⋅Pij⋅mij,where mij is the number of takeoff and departing aircraft that can be inserted in the interval between two consecutively arriving aircraft i and j on Runway 2.


## 3.1. Collision Risk Model

The collision risk model mainly considered the path deviation of the actual flight caused by the precision of navigation equipment, which is greatly affected by the aircraft speed. We can assume that both the navigation error and speed error conform to a Gaussian distribution and are given as follows:(5)εny∼Nμn,σn2,(6)εay∼N12μat2,14σa2t4,where μn and μa indicate the mean values and σn and σa are the standard deviations. Subscripts ny and ay denote the aircraft lateral deviation caused by two factors. The lateral collision risk can be computed using equations (7) to (10):(7)εy∼Nμy,σy2μy=μn+12μat2,σy2=σn2+14σa2t4,where εy denotes the distribution of the lateral deviation. The actual lateral distance between simultaneously departing aircraft is given as follows:(8)Dy′t=L2y′t−L1y′t=Dyt+ε2y−ε1y∼NDyt+μ2y−μ1y,σ2y2+σ1y2L1y′t=L1yt+ε1y,L2y′t=L2yt+ε2y,where L1y′t and L2y′t represent the actual path positions of two different aircraft. The probability density function for the lateral direction can be expressed as follows:(9)fy=12πσ2y2+σ1y2exp−y−Dyt+μ2y−μ1y22σ2y2+σ1y2.The probability density function of the actual distance between the two aircraft is integrated to obtain the lateral collision risk:(10)Py=py1<Dy′t<y2=∫y1y2fydy,where y1 and y2 represent the minimum and maximum separations between the aircraft, respectively, and can be calculated from the aircraft wing spans. Similarly, the probability density functions in the vertical and horizontal directions are given as follows:(11)fx=12πσ2x2+σ1x2exp−x−Dxt+μ2x−μ1x22σ2x2+σ1x2,(12)fz=12πσ2z2+σ1z2exp−z−Dzt+μ2z−μ1z22σ2z2+σ1z2,Finally, the collision risk model is obtained as follows:(13)P=∫x1x2fxdx∫y1y2fydy∫z1z2fzdz.


## 3.2. V-Open Multirunway Capacity Model

Wake turbulence was one of the main factors influencing the design of the V-open multirunway capacity model [23, 24]. At domestic hub airports, the number of medium aircraft accounts for 99% and heavy aircraft accounts for 1% at Chinese domestic hub airports. Therefore, in most cases, it is necessary to consider the wake impact for medium aircraft following heavy or medium aircraft, and detailed wake interval and accounts of different categories are shown in Table 2. Mainly, the wake generated from taking-off aircraft on the diverging runway affects the operation of departing aircraft on parallel runways and vice versa. A typical layout of V-open multiple runways is illustrated in Figure 1, in which there are two parallel runways and one diverging runway. There are generally two approach operations for parallel runways depending on whether or not they have different widths: independent parallel approach and correlated parallel approach [25]. Small airports in China currently use independent parallel approaches due to the high requirements for equipment and working skills. However, another diverging path of 30° for aircraft to safely go around also needs to be considered [26]. Therefore, the capacity evaluation model was proposed here based only on the correlated parallel approach.Table 2
Wake influence parameters.
CategorySH-SHSH-SMSM-SHSM-SMWake interval8 km6 km10 km6 kmAccounts4%16%16%64%Remarks: SH/SM denotes heavy and medium aircraft separately.An operation schematic diagram of parallel runways using the correlated parallel-approach-independent parallel departure mode is shown in Figure4. Both runways are used for aircraft arrivals and departures. Aircraft arrive and depart on each runway independently of each other. However, in order to ensure the safe operation of aircraft, the aircraft must land in the same direction on Runways 1 and 2. A safe distance needs to be maintained between aircraft on the parallel runway approaches δa.Figure 4
Parallel approach operations.Assuming that the distance between the two runways isD, as shown in Figure 4, the distance between two successively approaching aircraft on Runway 1 is the following:(14)lmin=2δa2−D2.Then, the time interval between two successively approaching aircraft a and aircraftc on Runway 1 is as follows:(15)tmin=maxlminv.The controller can use the insertion method to arrange for aircraftb to land on Runway 2 between aircraft a and aircraft c, but the time interval between aircraft a and aircraft c must be consistent with the following formula [27]:(16)Tac+Bac≥tmin.The operation sequence diagram of three runways for Runways 1 and 2 using related parallel-approach-independent parallel departure procedures is shown in Figure5. Assuming that approaching aircraft a on Runway 1 arrives at the runway threshold at 00:00, then within the interval between aircraft a and aircraft c, the departing aircraft can be released on Runway 1 or the approaching aircraft can be arranged to land on Runway 2. During time gap G2 of the continuous-approach aircraft on Runway 2, the departing aircraft can continue to be released, but the controller can also release the departing aircraft on Runway 3 after the departing aircraft on Runway 2 has been released for time tz.Figure 5
Time sequence model of V-open three-runway operation.Time intervalG1 between two consecutively approaching aircraft on Runway 1 is as follows:(17)G1=T1−2−T1−1=Tac+Bac−ROTa−maxtbc,tmin.WhenG1>0 (i.e., when the following conditions are satisfied), n departing aircraft can be inserted in gap G1:(18)ETac+Bac≥EROTa+n−1⋅Eεac+Emaxtbc,tmin+Eτ.Then, the entry and exit capacities of Runway 1 are as follows:(19)C1−arrival=3600ETac+Bac,(20)C1−departure=∑i=1n∑j=1nTG1⋅Pij⋅nij,where TG1 is the total number of clearances on Runway 1 within 1 h and nij is the number of departing aircraft that can be released between two adjacent approaching aircraft N/M on Runway 1.Time intervalG2 between two successively approaching aircraft on Runway 2 is as follows:(21)G2=T1−2−T1−1=Tij+Bij−maxROTb,th−maxt2,δvj.When the following conditions are met,k departing aircraft can be inserted in gap G2:(22)ETij+Bij≥EmaxROTb,th+n−1⋅Emaxtz+th,εij+Emaxt2,δvj+Eτ.Then, the entry and exit capacities of Runway 2 are as follows:(23)C2−arrival=∑i=1n∑j=1n∑xTG1GxETij+Bij,(24)C2−departure=∑i=1n∑j=1nTG2⋅Pij⋅kij,where Gx is the time interval between the first successive approach aircraft on Runway 1, TG1 and TG2 are the total numbers of clearances on Runways 1 and 2, respectively, within 1 h, and kij is the number of takeoff and departing aircraft that can be inserted in the time interval between two consecutively arriving aircraft i and j on Runway 2.In gapG2 of Runway 2, assuming that Runway 3 can release m aircraft, according to the operation sequence diagram, the following equation can be derived:(25)Tj−t2−T2−1≥m−1⋅maxTDd,tz+th,εij+τ,Therefore, the departure capacity of Runway 3 is as follows:(26)C3−departure=∑i=1n∑j=1nTG2⋅Pij⋅mij,where mij is the number of takeoff and departing aircraft that can be inserted in the interval between two consecutively arriving aircraft i and j on Runway 2.


## 4. Experiments and Analysis

This section describes associated experiments and data analysis for the runway layout at ZUTF. As presented in Figure6, this airport has two parallel runways separated by 2400 m and one diverging runway whose centerline is angled at 90°. The West Runway and East Runway are used for northbound departures and landings, and the North Runway is only used for eastward departures. The collision risk for a reduced divergence angle is evaluated to determine whether it satisfies the minimum requirement. The increase in the runway capacity is determined by a theoretical analysis and further tested using simulation experiments.Figure 6
V-open runway layout at ZUTF.


### 4.1. Analysis of Collision Risk

The reference distance from the runway end is generally set as 10 nm (18.6 km) to calculate the reference spacing when applying the ELSO standard. Then, by analyzing the actual departure path at ZUTF (see Figure7), which contains 998 flight data all departing from parallel runways on December 14 and 15, 2021, the path variance angle is obtained as σSTO,C = 0.4° and σDIV,C = 0.6°. In this case, minimum lateral spacing sd can be easily obtained using equation (2) as the conventional divergence standard.Figure 7
Analysis of actual departure paths at ZUTF.The width of the parallel runways at ZUTF is 2400 m and the runways are staggered by 1350 m. In this case, equation (4) can be used to evaluate the effect of reducing the divergence angle to 10° according to the ELSO standard.Experiments on the collision risk were performed for the proposed risk evaluation model using MATLAB software (version 2020). The experimental parameters for aircraft departing from the parallel runways with a separation path angle of 10° and three types of aircraft (heavy, medium, and light) are listed in Table3.Table 3
Performance parameters of different aircraft.
Aircraft typeLength offuselage (m)Wing span (m)Height offuselage (m)Final approachspeed (km/h)Runway thresholdspeed (km/h)Takeoff speed(km/h)Normal cruisingspeed (km/h)Heavy68.564.312.55296287293916Medium39.534.319.7258249269839Light16.5516.974.86213204189333When the crosswind speed is 6 kt, the aircraft wake moves at 15 ft/s. The time needed from aircraft takeoff to the first navigation point during the departure procedure is about 120 s. Research has shown that the lateral movement distance of the wake at a wind speed of 6 kt is 640 m.The parameters of the collision risk model are listed in Table4. The upper limit of the collision probability integration interval is given by the following equation:(27)y1=λ1y+λ2y2.Table 4
Collision risk calculation parameters.
ParameterValueμn0σn3μa0σa1.27R2400σg50Dzt300The lower limits of the collision probability integration interval in the lateral, longitudinal, and vertical directions are given by the following equation:(28)y2=miny1+640,2400−y1x2=x1+5+v¯former.flight×120z2=minz1+2+2×120,300−z1.According to the collision risk model based on the calculation parameters, the collision risk of departing aircraft on the parallel West Runway and East Runway can be calculated as presented in Table5.Table 5
Probabilities of collision risk with different divergence angles.
Divergence angle (degree)Aircraft typeSide collision riskLongitudinal collision riskVertical collision riskTotal collision risk14Heavy-heavy7.4268 × 10−62.7528 × 10−101.6612 × 10−53.3872 × 10−20Heavy-medium7.3578 × 10−62.7415 × 10−101.6844 × 10−53.3833 × 10−20Medium-medium7.1581 × 10−62.6924 × 10−101.5512 × 10−52.9811 × 10−2013Heavy-heavy9.6952 × 10−62.6852 × 10−101.6512 × 10−54.2849 × 10−20Heavy-medium9.5586 × 10−62.6964 × 10−101.6300 × 10−54.1873 × 10−20Medium-medium9.3254 × 10−62.5368 × 10−101.5946 × 10−53.7491 × 10−2012Heavy-heavy2.7485 × 10−52.7256 × 10−101.6102 × 10−51.1736 × 10−20Heavy-medium2.5263 × 10−52.7415 × 10−101.5968 × 10−51.0732 × 10−20Medium-medium2.1452 × 10−52.7238 × 10−101.5317 × 10−58.6751 × 10−1911Heavy-heavy4.5685 × 10−52.6852 × 10−101.6012 × 10−51.9553 × 10−19Heavy-medium4.3947 × 10−52.6732 × 10−101.5858 × 10−51.7664 × 10−19Medium-medium4.0359 × 10−52.5845 × 10−101.5577 × 10−51.6075 × 10−1910Heavy-heavy7.3195 × 10−52.8601 × 10−101.6302 × 10−53.4127 × 10−19Heavy-medium7.2873 × 10−52.8586 × 10−101.6400 × 10−53.4163 × 10−19Medium-medium7.2647 × 10−52.7985 × 10−101.32852 × 10−52.7009 × 10−19The results in Table5 indicate that when the West Runway and East Runway at ZUTF are operated in an independent parallel departure mode with a reduced divergence angle from 14 degree to 10°, the probability of aircraft collision risk is all less than the overall safety target level (1.0 × 10−8). Therefore, it is considered that the operation spacing for the departing aircraft is safe and that the divergence angle for the simultaneous use of the two long-distance West Runway and East Runway at ZUTF can be reduced to 10°.


### 4.2. V-Open Multirunway Capacity Evaluation

#### 4.2.1. Theoretical Analysis

When the departing aircraft on the East Runway has a divergence angle of 15°, the aircraft departure protection area is 724 m from the North Runway. Therefore, the aircraft on the East Runway will affect the departure of aircraft on the North Runway.After the aircraft on the North Runway has taxied forth=25s, the air traffic controller can allow the release of the departing aircraft on the East Runway. The departing aircraft on the East Runway is released at tz=70s, and the departing aircraft on the North Runway can then be aligned to the runway and allowed to taxi.When applying the ELSO standard, the departure divergence angle is reduced to 10°. In this case, the departure protection zone of the departing aircraft on the East Runway still overlaps with the head of the North Runway, but the overlap distance is reduced to 494 m. Because departing aircraft on the North Runway have a taxiing distance of 600 m behind the runway head, the departing aircraft on East Runway are not limited by the divergence angle of 10°.The improvements in the runway capacity when using the proposed model and associated parameters are given in Table6. The capacity results for the conventional standard (divergence angle of 15°) and the ELSO standard (divergence angle of 10°) are presented in Table 7. It is clear that the number of flights per hour can increase by six with the reduced divergence angle.Table 6
Runway capacity calculation parameters.
ParameterValueDD interval100 sDA interval180 sAA interval100 sMinimum distance between landing aircraft and entrance95 sAverage runway occupation time of approaching aircraftROTi.50 sAverage runway occupation time of departing aircraftTDi.55 sFinal approach speed of heavy aircraft300 km/hFinal approach speed of medium aircraft270 km/hσVheavy0.14 (km/h)σVmedium0.14 (km/h)Table 7
Calculation results for runway capacity.
Runway capacity (aircraft/h)Related parallel-approach-independent parallel departureExistingOptimizedWest Runway3232East Runway2629North Runway3033Total capacity8894


#### 4.2.2. Simulation Results

In order to test the theoretical results, several simulation experiments were performed for ZUTF based on the general airport simulation software jointly developed using the Visual Basic programming language. This airport surface simulation software has been approved by the ATC office of the CAAC. The airport runway layout was first set to summarize the measurement data of the airport ground flight area, and AutoCAD software was used to build a ground simulation model. After standardizing the positioning point information, two-dimensional coordinate information, taxi-path constraint rules, and parking space information, these data were imported into the ground capacity evaluation system to establish the ground simulation capacity evaluation model at ZUTF. The simulation interface is shown in Figure8.Figure 8
Ground simulation interface.In accordance with the regulations for ZUTF operation, the North Runway, East Runway, and West Runway account for 70%, 5%, and 25%, respectively, of all departing flights, while the East Runway and West Runway account for 83% and 17%, respectively, of all arriving flights.The flight schedule is selected as the data source of the simulation evaluation model, and the ground parameters used in the simulation program are listed in Table8. Considering the characteristics of runway operation and operating rules, and under the restriction conditions of the percentages of the numbers of arriving and departing flights, the flight flow distribution proportion for each hour, and other characteristics, the original flight schedule of a typical day was randomly selected to generate pressurized flight and aircraft flows for use in the simulations.Table 8
Simulation operating parameters.
Simulation parameterValueTakeoff runway occupation time50 sLanding runway occupation time (fast departure lane)50 sLanding runway occupation time (end of runway)75 sMinimum time interval between consecutive takeoffs100 sMinimum time interval for continuous landings100 sMinimum distance between the landing aircraft and the entrance required to issue the takeoff permission95 sTime to cross the taxiway starting from standstill40 sRunway threshold speed150 ktRapid departure speed30 ktMain slide and vertical connecting track speed20 ktGeneral taxiway speed10 ktTaxi speed in apron area5 ktTime from the sliding position withdrawing the wheel to the completion of driving150 sTable9 presents the results from the simulations of the air-ground joint operation capacity of the runway system at ZUTF under the conventional standard (divergence angle of 15°) and the ELSO standard (divergence angle of 10°). The runway capacity under the conventional standard was 74 per hour, which is less than that under the proposed ELSO standard (81 per hour). There were 8.1% more departing flights for the proposed policy during the peak hour and 11.1% more during the peak quarter than the conventional one. Moreover, the average flight delay time reduced significantly by 30% when using the proposed policy. Furthermore, given that 95% of flights involved medium aircraft and 5% involved heavy aircraft with an average passenger load factor of 85%, it can be inferred that passenger turnover at ZUTF would increase by about 6.3%.Table 9
Simulation results under the conventional standard and the ELSO standard.
Operation modeFlight number in peak hourTotal flight numberAverage delay (min)Delay in the air (min)Delay on the ground (min)Peak hour delay (min)Flight numberin peak quarterConventional standard7411897.94.24.713.918ELSO standard8112646.83.03.813.620Improvement+8.1%+6.3%–13.9%–30%–19%+2.1%+11.1%The present theoretical analysis and simulations have demonstrated that the efficiency of ZUTF can be markedly improved when the divergence angle of aircraft on the East Runway is reduced from 15° to 10°. This results in no overlap in the ground starting area, and the departure operations on the two runways are independent of each other. During operations, it is only necessary to consider that the go-around aircraft on the East Runway will release aircraft on the North Runway before 4.2 km from the runway entrance, which can effectively improve runway capacity by about 10%.In summary, the proposed policy can improve (1) airport capacity for either parallel runways or diverging runways, (2) the travel experience in civil aviation with a higher flight punctuality rate, and (3) airspace resource optimization with more flexible operation procedures. The present research results can provide theoretical guidance for the future development strategies of the CAAC.


## 4.1. Analysis of Collision Risk

The reference distance from the runway end is generally set as 10 nm (18.6 km) to calculate the reference spacing when applying the ELSO standard. Then, by analyzing the actual departure path at ZUTF (see Figure7), which contains 998 flight data all departing from parallel runways on December 14 and 15, 2021, the path variance angle is obtained as σSTO,C = 0.4° and σDIV,C = 0.6°. In this case, minimum lateral spacing sd can be easily obtained using equation (2) as the conventional divergence standard.Figure 7
Analysis of actual departure paths at ZUTF.The width of the parallel runways at ZUTF is 2400 m and the runways are staggered by 1350 m. In this case, equation (4) can be used to evaluate the effect of reducing the divergence angle to 10° according to the ELSO standard.Experiments on the collision risk were performed for the proposed risk evaluation model using MATLAB software (version 2020). The experimental parameters for aircraft departing from the parallel runways with a separation path angle of 10° and three types of aircraft (heavy, medium, and light) are listed in Table3.Table 3
Performance parameters of different aircraft.
Aircraft typeLength offuselage (m)Wing span (m)Height offuselage (m)Final approachspeed (km/h)Runway thresholdspeed (km/h)Takeoff speed(km/h)Normal cruisingspeed (km/h)Heavy68.564.312.55296287293916Medium39.534.319.7258249269839Light16.5516.974.86213204189333When the crosswind speed is 6 kt, the aircraft wake moves at 15 ft/s. The time needed from aircraft takeoff to the first navigation point during the departure procedure is about 120 s. Research has shown that the lateral movement distance of the wake at a wind speed of 6 kt is 640 m.The parameters of the collision risk model are listed in Table4. The upper limit of the collision probability integration interval is given by the following equation:(27)y1=λ1y+λ2y2.Table 4
Collision risk calculation parameters.
ParameterValueμn0σn3μa0σa1.27R2400σg50Dzt300The lower limits of the collision probability integration interval in the lateral, longitudinal, and vertical directions are given by the following equation:(28)y2=miny1+640,2400−y1x2=x1+5+v¯former.flight×120z2=minz1+2+2×120,300−z1.According to the collision risk model based on the calculation parameters, the collision risk of departing aircraft on the parallel West Runway and East Runway can be calculated as presented in Table5.Table 5
Probabilities of collision risk with different divergence angles.
Divergence angle (degree)Aircraft typeSide collision riskLongitudinal collision riskVertical collision riskTotal collision risk14Heavy-heavy7.4268 × 10−62.7528 × 10−101.6612 × 10−53.3872 × 10−20Heavy-medium7.3578 × 10−62.7415 × 10−101.6844 × 10−53.3833 × 10−20Medium-medium7.1581 × 10−62.6924 × 10−101.5512 × 10−52.9811 × 10−2013Heavy-heavy9.6952 × 10−62.6852 × 10−101.6512 × 10−54.2849 × 10−20Heavy-medium9.5586 × 10−62.6964 × 10−101.6300 × 10−54.1873 × 10−20Medium-medium9.3254 × 10−62.5368 × 10−101.5946 × 10−53.7491 × 10−2012Heavy-heavy2.7485 × 10−52.7256 × 10−101.6102 × 10−51.1736 × 10−20Heavy-medium2.5263 × 10−52.7415 × 10−101.5968 × 10−51.0732 × 10−20Medium-medium2.1452 × 10−52.7238 × 10−101.5317 × 10−58.6751 × 10−1911Heavy-heavy4.5685 × 10−52.6852 × 10−101.6012 × 10−51.9553 × 10−19Heavy-medium4.3947 × 10−52.6732 × 10−101.5858 × 10−51.7664 × 10−19Medium-medium4.0359 × 10−52.5845 × 10−101.5577 × 10−51.6075 × 10−1910Heavy-heavy7.3195 × 10−52.8601 × 10−101.6302 × 10−53.4127 × 10−19Heavy-medium7.2873 × 10−52.8586 × 10−101.6400 × 10−53.4163 × 10−19Medium-medium7.2647 × 10−52.7985 × 10−101.32852 × 10−52.7009 × 10−19The results in Table5 indicate that when the West Runway and East Runway at ZUTF are operated in an independent parallel departure mode with a reduced divergence angle from 14 degree to 10°, the probability of aircraft collision risk is all less than the overall safety target level (1.0 × 10−8). Therefore, it is considered that the operation spacing for the departing aircraft is safe and that the divergence angle for the simultaneous use of the two long-distance West Runway and East Runway at ZUTF can be reduced to 10°.


## 4.2. V-Open Multirunway Capacity Evaluation

### 4.2.1. Theoretical Analysis

When the departing aircraft on the East Runway has a divergence angle of 15°, the aircraft departure protection area is 724 m from the North Runway. Therefore, the aircraft on the East Runway will affect the departure of aircraft on the North Runway.After the aircraft on the North Runway has taxied forth=25s, the air traffic controller can allow the release of the departing aircraft on the East Runway. The departing aircraft on the East Runway is released at tz=70s, and the departing aircraft on the North Runway can then be aligned to the runway and allowed to taxi.When applying the ELSO standard, the departure divergence angle is reduced to 10°. In this case, the departure protection zone of the departing aircraft on the East Runway still overlaps with the head of the North Runway, but the overlap distance is reduced to 494 m. Because departing aircraft on the North Runway have a taxiing distance of 600 m behind the runway head, the departing aircraft on East Runway are not limited by the divergence angle of 10°.The improvements in the runway capacity when using the proposed model and associated parameters are given in Table6. The capacity results for the conventional standard (divergence angle of 15°) and the ELSO standard (divergence angle of 10°) are presented in Table 7. It is clear that the number of flights per hour can increase by six with the reduced divergence angle.Table 6
Runway capacity calculation parameters.
ParameterValueDD interval100 sDA interval180 sAA interval100 sMinimum distance between landing aircraft and entrance95 sAverage runway occupation time of approaching aircraftROTi.50 sAverage runway occupation time of departing aircraftTDi.55 sFinal approach speed of heavy aircraft300 km/hFinal approach speed of medium aircraft270 km/hσVheavy0.14 (km/h)σVmedium0.14 (km/h)Table 7
Calculation results for runway capacity.
Runway capacity (aircraft/h)Related parallel-approach-independent parallel departureExistingOptimizedWest Runway3232East Runway2629North Runway3033Total capacity8894


### 4.2.2. Simulation Results

In order to test the theoretical results, several simulation experiments were performed for ZUTF based on the general airport simulation software jointly developed using the Visual Basic programming language. This airport surface simulation software has been approved by the ATC office of the CAAC. The airport runway layout was first set to summarize the measurement data of the airport ground flight area, and AutoCAD software was used to build a ground simulation model. After standardizing the positioning point information, two-dimensional coordinate information, taxi-path constraint rules, and parking space information, these data were imported into the ground capacity evaluation system to establish the ground simulation capacity evaluation model at ZUTF. The simulation interface is shown in Figure8.Figure 8
Ground simulation interface.In accordance with the regulations for ZUTF operation, the North Runway, East Runway, and West Runway account for 70%, 5%, and 25%, respectively, of all departing flights, while the East Runway and West Runway account for 83% and 17%, respectively, of all arriving flights.The flight schedule is selected as the data source of the simulation evaluation model, and the ground parameters used in the simulation program are listed in Table8. Considering the characteristics of runway operation and operating rules, and under the restriction conditions of the percentages of the numbers of arriving and departing flights, the flight flow distribution proportion for each hour, and other characteristics, the original flight schedule of a typical day was randomly selected to generate pressurized flight and aircraft flows for use in the simulations.Table 8
Simulation operating parameters.
Simulation parameterValueTakeoff runway occupation time50 sLanding runway occupation time (fast departure lane)50 sLanding runway occupation time (end of runway)75 sMinimum time interval between consecutive takeoffs100 sMinimum time interval for continuous landings100 sMinimum distance between the landing aircraft and the entrance required to issue the takeoff permission95 sTime to cross the taxiway starting from standstill40 sRunway threshold speed150 ktRapid departure speed30 ktMain slide and vertical connecting track speed20 ktGeneral taxiway speed10 ktTaxi speed in apron area5 ktTime from the sliding position withdrawing the wheel to the completion of driving150 sTable9 presents the results from the simulations of the air-ground joint operation capacity of the runway system at ZUTF under the conventional standard (divergence angle of 15°) and the ELSO standard (divergence angle of 10°). The runway capacity under the conventional standard was 74 per hour, which is less than that under the proposed ELSO standard (81 per hour). There were 8.1% more departing flights for the proposed policy during the peak hour and 11.1% more during the peak quarter than the conventional one. Moreover, the average flight delay time reduced significantly by 30% when using the proposed policy. Furthermore, given that 95% of flights involved medium aircraft and 5% involved heavy aircraft with an average passenger load factor of 85%, it can be inferred that passenger turnover at ZUTF would increase by about 6.3%.Table 9
Simulation results under the conventional standard and the ELSO standard.
Operation modeFlight number in peak hourTotal flight numberAverage delay (min)Delay in the air (min)Delay on the ground (min)Peak hour delay (min)Flight numberin peak quarterConventional standard7411897.94.24.713.918ELSO standard8112646.83.03.813.620Improvement+8.1%+6.3%–13.9%–30%–19%+2.1%+11.1%The present theoretical analysis and simulations have demonstrated that the efficiency of ZUTF can be markedly improved when the divergence angle of aircraft on the East Runway is reduced from 15° to 10°. This results in no overlap in the ground starting area, and the departure operations on the two runways are independent of each other. During operations, it is only necessary to consider that the go-around aircraft on the East Runway will release aircraft on the North Runway before 4.2 km from the runway entrance, which can effectively improve runway capacity by about 10%.In summary, the proposed policy can improve (1) airport capacity for either parallel runways or diverging runways, (2) the travel experience in civil aviation with a higher flight punctuality rate, and (3) airspace resource optimization with more flexible operation procedures. The present research results can provide theoretical guidance for the future development strategies of the CAAC.


## 4.2.1. Theoretical Analysis

When the departing aircraft on the East Runway has a divergence angle of 15°, the aircraft departure protection area is 724 m from the North Runway. Therefore, the aircraft on the East Runway will affect the departure of aircraft on the North Runway.After the aircraft on the North Runway has taxied forth=25s, the air traffic controller can allow the release of the departing aircraft on the East Runway. The departing aircraft on the East Runway is released at tz=70s, and the departing aircraft on the North Runway can then be aligned to the runway and allowed to taxi.When applying the ELSO standard, the departure divergence angle is reduced to 10°. In this case, the departure protection zone of the departing aircraft on the East Runway still overlaps with the head of the North Runway, but the overlap distance is reduced to 494 m. Because departing aircraft on the North Runway have a taxiing distance of 600 m behind the runway head, the departing aircraft on East Runway are not limited by the divergence angle of 10°.The improvements in the runway capacity when using the proposed model and associated parameters are given in Table6. The capacity results for the conventional standard (divergence angle of 15°) and the ELSO standard (divergence angle of 10°) are presented in Table 7. It is clear that the number of flights per hour can increase by six with the reduced divergence angle.Table 6
Runway capacity calculation parameters.
ParameterValueDD interval100 sDA interval180 sAA interval100 sMinimum distance between landing aircraft and entrance95 sAverage runway occupation time of approaching aircraftROTi.50 sAverage runway occupation time of departing aircraftTDi.55 sFinal approach speed of heavy aircraft300 km/hFinal approach speed of medium aircraft270 km/hσVheavy0.14 (km/h)σVmedium0.14 (km/h)Table 7
Calculation results for runway capacity.
Runway capacity (aircraft/h)Related parallel-approach-independent parallel departureExistingOptimizedWest Runway3232East Runway2629North Runway3033Total capacity8894


## 4.2.2. Simulation Results

In order to test the theoretical results, several simulation experiments were performed for ZUTF based on the general airport simulation software jointly developed using the Visual Basic programming language. This airport surface simulation software has been approved by the ATC office of the CAAC. The airport runway layout was first set to summarize the measurement data of the airport ground flight area, and AutoCAD software was used to build a ground simulation model. After standardizing the positioning point information, two-dimensional coordinate information, taxi-path constraint rules, and parking space information, these data were imported into the ground capacity evaluation system to establish the ground simulation capacity evaluation model at ZUTF. The simulation interface is shown in Figure8.Figure 8
Ground simulation interface.In accordance with the regulations for ZUTF operation, the North Runway, East Runway, and West Runway account for 70%, 5%, and 25%, respectively, of all departing flights, while the East Runway and West Runway account for 83% and 17%, respectively, of all arriving flights.The flight schedule is selected as the data source of the simulation evaluation model, and the ground parameters used in the simulation program are listed in Table8. Considering the characteristics of runway operation and operating rules, and under the restriction conditions of the percentages of the numbers of arriving and departing flights, the flight flow distribution proportion for each hour, and other characteristics, the original flight schedule of a typical day was randomly selected to generate pressurized flight and aircraft flows for use in the simulations.Table 8
Simulation operating parameters.
Simulation parameterValueTakeoff runway occupation time50 sLanding runway occupation time (fast departure lane)50 sLanding runway occupation time (end of runway)75 sMinimum time interval between consecutive takeoffs100 sMinimum time interval for continuous landings100 sMinimum distance between the landing aircraft and the entrance required to issue the takeoff permission95 sTime to cross the taxiway starting from standstill40 sRunway threshold speed150 ktRapid departure speed30 ktMain slide and vertical connecting track speed20 ktGeneral taxiway speed10 ktTaxi speed in apron area5 ktTime from the sliding position withdrawing the wheel to the completion of driving150 sTable9 presents the results from the simulations of the air-ground joint operation capacity of the runway system at ZUTF under the conventional standard (divergence angle of 15°) and the ELSO standard (divergence angle of 10°). The runway capacity under the conventional standard was 74 per hour, which is less than that under the proposed ELSO standard (81 per hour). There were 8.1% more departing flights for the proposed policy during the peak hour and 11.1% more during the peak quarter than the conventional one. Moreover, the average flight delay time reduced significantly by 30% when using the proposed policy. Furthermore, given that 95% of flights involved medium aircraft and 5% involved heavy aircraft with an average passenger load factor of 85%, it can be inferred that passenger turnover at ZUTF would increase by about 6.3%.Table 9
Simulation results under the conventional standard and the ELSO standard.
Operation modeFlight number in peak hourTotal flight numberAverage delay (min)Delay in the air (min)Delay on the ground (min)Peak hour delay (min)Flight numberin peak quarterConventional standard7411897.94.24.713.918ELSO standard8112646.83.03.813.620Improvement+8.1%+6.3%–13.9%–30%–19%+2.1%+11.1%The present theoretical analysis and simulations have demonstrated that the efficiency of ZUTF can be markedly improved when the divergence angle of aircraft on the East Runway is reduced from 15° to 10°. This results in no overlap in the ground starting area, and the departure operations on the two runways are independent of each other. During operations, it is only necessary to consider that the go-around aircraft on the East Runway will release aircraft on the North Runway before 4.2 km from the runway entrance, which can effectively improve runway capacity by about 10%.In summary, the proposed policy can improve (1) airport capacity for either parallel runways or diverging runways, (2) the travel experience in civil aviation with a higher flight punctuality rate, and (3) airspace resource optimization with more flexible operation procedures. The present research results can provide theoretical guidance for the future development strategies of the CAAC.


## 5. Conclusion

Limited by the conventional separation standard of 15 degrees, V-open multiple runways cannot perform its maximum operation function. A collision risk model is proposed to assess the feasibility and safety of the proposed policy at ZUTF airport. A novel capacity estimation model for V-open multiple runways is built to test the operating efficiency of the proposed new separation standard relative to the conventional one. Comparative experiments based on actual flight data from Chengdu Tianfu International Airport showed that the probability of aircraft collision risk is all less than the overall safety target level when divergence angle reduced from 15 degree to 10 degree and new capacity model effectively improves runway capacity by about 10%.The present research findings can provide theoretical guidance for the CAAC to apply an operation standard to future V-open multiple runways. They might also be useful as a reference for capacity evaluations at airports with diverging runways used for civil aviation. In the future, when divergence angle is reducing from 15 degree to 10 degree, the proposed ELSO standard can also be applied in other airports with three parallel runways to improve the efficiency of airspace in China. Furthermore, the trajectories of go-around aircraft could also cover the diverging runway area for independent parallel approaching operation. Therefore, we will attempt to modify our capacity model with such considerations in our future work.



---

*Source: 5869101-2022-10-12.xml*