Air traffic stability control method

Abstract

The invention discloses an air traffic stability control method, which belongs to the technical field of intelligent management of air traffic management, and comprises the following steps: constructing an aircraft flight space-time sequence based on high-altitude cruise aircraft information; calculating a relative flight position interval and a relative flight speed between the aircrafts based on the flight space-time sequence of the aircrafts; based on the relative flight position interval and the relative flight speed between the aircrafts, constructing an aircraft car-following model containing time-delay band bidirectional attention; based on the aircraft car-following model, constructing air traffic stability judgment conditions; and based on the air traffic stability judgment condition, the aircrafts in the unstable state are screened for feedback control, and the air traffic stability is guaranteed. According to the method, the air traffic stability condition can be judged, feedback control is performed by screening aircrafts in an unstable state, and the air traffic stability is guaranteed.

Description

Technical Field

[0001] This invention relates to a method for controlling air traffic stability, belonging to the field of intelligent air traffic management technology. Background Technology

[0002] With the rapid development of the global air transport industry, the number of aircraft in the high-altitude cruise phase continues to increase. In the field of intelligent air traffic management, effectively controlling the operational status of high-altitude cruise aircraft has become a research hotspot. Especially among multiple aircraft operating on the same route, which have high speeds and strong coupling, even a minor disturbance to one aircraft can easily trigger a chain reaction, causing an imbalance in overall air traffic stability. This further exacerbates the complexity and uncertainty of the air traffic system, affecting route operational efficiency and flight safety. Therefore, it is urgent to design highly coordinated control strategies to achieve dynamic collaboration and safe, efficient operation among multiple aircraft, ensuring the overall stability and safety of the air traffic system.

[0003] Existing research on air traffic stability control rarely considers modeling of car-following relationships between aircraft. Many studies rely on simplified assumptions to build car-following models, failing to capture the real-world time delay effects, attention allocation mechanisms, and two-way influence characteristics present in high-altitude cruise environments, resulting in insufficient model accuracy and adaptability. Consequently, air traffic stability assessment mechanisms are incomplete, lacking stability analysis methods based on the evolution of car-following behavior, making it difficult to identify potential systemic instability risks in a timely manner. Furthermore, the lack of closed-loop stability control strategies limits the means to control unstable aircraft, making it difficult to effectively ensure air traffic stability in complex operating environments. Summary of the Invention

[0004] The purpose of this invention is to provide an air traffic stability control method that can determine the air traffic stability status and perform feedback control on aircraft in unstable states to ensure air traffic stability.

[0005] To achieve the above objectives, the present invention provides the following technical solution: In a first aspect, the present invention provides an air traffic stability control method, comprising: Construct the spatiotemporal sequence of aircraft flight based on high-altitude cruise aircraft information; Based on the spatiotemporal sequence of aircraft flight, calculate the relative flight position interval and relative flight speed between aircraft; Based on the relative flight position interval and relative flight speed between aircraft, an aircraft car-following model with time delay and bidirectional attention is constructed. Based on the aircraft car-following model, conditions for determining air traffic stability are constructed. Based on the criteria for determining air traffic stability, aircraft in unstable states are selected for feedback control to ensure air traffic stability.

[0006] In conjunction with the first aspect, further, based on high-altitude cruise aircraft information, the construction of the aircraft's flight spatiotemporal sequence includes: Based on high-altitude cruise aircraft information, extract aircraft flight observation data streams; The flight observation data stream of aircraft is used to reconstruct the trajectory, project the route, and sort the queue to generate the spatiotemporal sequence of aircraft flight.

[0007] In conjunction with the first aspect, further, based on high-altitude cruise aircraft information, the extraction of aircraft flight observation data streams includes: Extract multi-source observation data based on information from high-altitude cruise aircraft; The timestamps of multi-source observation data of the same aircraft are identified and associated according to the aircraft call sign, forming an aircraft flight observation data stream that includes the route, three-dimensional position, speed and heading of each timestamp of the aircraft; The multi-source observation data includes flight plan messages, as well as route, segment, and aircraft information from flight information service systems, global navigation satellite systems, air-to-ground and air-to-air communication systems, secondary surveillance radar systems, and automatic dependent surveillance broadcast systems.

[0008] In conjunction with the first aspect, further steps include trajectory reconstruction, route projection, and queue sorting of the aircraft flight observation data stream to generate the aircraft's spatiotemporal flight sequence, including: The observation points for each aircraft are resampled, and the time step of the aircraft trajectory is standardized to obtain the continuous trajectory of the aircraft. The positions corresponding to each continuous trajectory at the time step are selected, and the arc length projection of the continuous trajectory is calculated along the centerline of the flight path to obtain the heading coordinates of the aircraft. The aircraft are then sorted in ascending order according to their heading coordinates to form a flight time-space sequence.

[0009] In conjunction with the first aspect, further, based on the aircraft's flight time-space sequence, the calculation of the relative flight position intervals and relative flight speeds between aircraft includes: The flight time-space sequence of aircraft is constructed into a formation sequence diagram. Based on the formation sequence diagram, adjacent aircraft are selected, and the relative flight position interval and relative flight speed between adjacent aircraft are calculated. The queue sequence diagram is a spatial topology constructed with aircraft as nodes and the adjacency relationships between aircraft as directed edges.

[0010] In conjunction with the first aspect, the relative flight position interval between adjacent aircraft is further defined as follows: ; in, express Time of the first The aircraft and the first The relative flight position interval between aircraft , They represent Time of the first , The flight positions of the aircraft; The relative flight speeds between adjacent aircraft are: ; in, express Time of the first The aircraft and the first The relative flight speed between aircraft , They represent Time of the first , The flight speed of an aircraft , They represent Time of the first , The flight positions of the aircraft.

[0011] Building upon the first aspect, further, based on the relative flight position intervals and relative flight speeds between aircraft, an aircraft car-following model with time delay and bidirectional attention is constructed, including: Calculate the desired positional spacing between aircraft and use the relative flight positional spacing between aircraft as the actual positional spacing between aircraft; Based on the expected and actual positional intervals between aircraft, calculate the forward and backward following strengths of the aircraft. Based on the forward and backward car-following strengths of the aircraft, calculate the forward and backward attention distributions of the pilots on the aircraft. A two-way attention distribution model is constructed from the forward and backward attention distributions of pilots on aircraft; Construct an air-to-air communication and response time delay model; Based on the relative flight speed between aircraft, a two-way attention distribution model and an air-to-air communication and response time delay model are combined to obtain an aircraft car-following model with time delay and two-way attention to describe the car-following characteristics of aircraft.

[0012] In conjunction with the first aspect, the desired positional spacing between aircraft is further as follows: ; in, , They represent Time of the first The aircraft and the one in front of it One aircraft, rear of the first The desired positional spacing between aircraft , , They represent Time of the first The aircraft, the The first aircraft in front The aircraft, the The aircraft behind The expected flight position of an aircraft; The actual positional spacing between aircraft is: ; in, , They represent Time of the first The aircraft and the first , The relative flight position interval between aircraft, i.e. Time of the first The aircraft and the one in front of it One aircraft, rear of the first The actual positional interval between aircraft , , They represent Time of the first , , The flight position of an aircraft, i.e. Time of the first The aircraft, the The first aircraft in front The aircraft, the The aircraft behind The actual flight position of the aircraft; The forward following strength of the aircraft is: ; in, express Time of the first The aircraft relative to the first one in front of it Forward following strength of an aircraft; The rearward following strength of the aircraft is: ; in, express Time of the first The aircraft relative to the one behind it The rearward following strength of an aircraft; The forward attention distribution of a pilot on an aircraft is as follows: ; in, express Time of the first The pilots on the aircraft were looking at the first... Forward attention value of an aircraft Indicates the first The total number of aircraft ahead of each aircraft; The distribution of a pilot's backward attention on an aircraft is as follows: ; in, express Time of the first The pilots on the aircraft were looking at the first... The backward attention value of an aircraft Indicates the first The total number of aircraft behind each aircraft; The air-to-air communication and response time delay model is as follows: ; in, This indicates air-to-air communication and response delay. This indicates the time required for air-to-air communication. Follows uniform distribution , , , They represent The minimum, maximum, and average values ​​of the value. This indicates the time it takes for the pilots on board the aircraft to respond to and manipulate the aircraft based on communication information. Follows uniform distribution , , , They represent Minimum, maximum, and average values; The aircraft following model is: ; in, express Time of the first The flight acceleration of an aircraft Indicates the first The sensitivity coefficient of an aircraft's response based on the status of the aircraft in front and behind it. , They represent Time of the first The aircraft and the first , The relative flight position interval between aircraft, i.e. Time of the first The aircraft and the one in front of it One aircraft, rear of the first The actual positional interval between aircraft express Time of the first The flight speed of an aircraft , They represent Time of the first The aircraft and the first , The relative flight speed between aircraft, i.e. Time of the first The aircraft and the one in front of it One aircraft, rear of the first The relative flight speed between aircraft , They represent Time of the first The aircraft and the one in front of it One aircraft, rear of the first Relative flight acceleration between aircraft Indicates the safe separation between aircraft. Represents the hyperbolic tangent function. Indicates according to and The calculated contents The The aircraft and the one in front of it The optimal flight speed for maintaining a safe distance between aircraft. Indicates according to and The calculated contents The The aircraft and the one behind it The optimal flight speed for maintaining a safe distance between aircraft. Indicates according to and The calculated first The aircraft and the one in front of it The optimal flight speed for maintaining a safe distance between aircraft. Indicates according to and The calculated first The aircraft and the one behind it The optimal flight speed for maintaining a safe distance between aircraft. , They represent , The first derivative, This represents the maximum value of the optimal flight speed.

[0013] In conjunction with the first aspect, the further criteria for determining air traffic stability are: ; in, Indicates the air traffic stability threshold. Indicates according to The first derivative of the determined optimal flight speed is the optimal flight acceleration.

[0014] Building upon the first aspect, further, based on the air traffic stability assessment criteria, aircraft in unstable states are selected for feedback control to ensure air traffic stability, including: Aircraft that do not meet the criteria for determining air traffic stability are considered to be in an unstable state. The system controls aircraft in unstable states to automatically change their flight acceleration according to the aircraft following model in order to adjust the flight position interval between aircraft and ensure air traffic stability. Among them, aircraft that meet the air traffic stability determination conditions are considered to be aircraft in a stable state; Control the aircraft to accelerate or maintain its current flight speed in a stable state.

[0015] In a second aspect, the present invention provides an air traffic stability control device, comprising: The aircraft flight spatiotemporal sequence construction module is used to construct the aircraft flight spatiotemporal sequence based on high-altitude cruise aircraft information; The aircraft relative relationship analysis module is used to calculate the relative flight position interval and relative flight speed between aircraft based on the aircraft flight time-space sequence; The car-following model building module is used to build an aircraft car-following model with time delay and bidirectional attention based on the relative flight position interval and relative flight speed between aircraft. The stability determination module is used to construct air traffic stability determination conditions based on the aircraft car-following model. The feedback control module is used to screen aircraft in unstable states based on air traffic stability determination conditions and perform feedback control to ensure air traffic stability.

[0016] Thirdly, the present invention provides a computer device, comprising: Storage medium used to store computer programs; A processor for executing the computer program to implement the air traffic stability control method of the first aspect.

[0017] Fourthly, the present invention provides a computer-readable storage medium having a computer program stored thereon, wherein the computer program, when executed by a processor, implements the air traffic stability control method described in the first aspect.

[0018] Fifthly, the present invention provides a computer program product, including a computer program that, when executed by a processor, implements the air traffic stability control method described in the first aspect.

[0019] Compared with the prior art, the beneficial effects of the present invention are: The air traffic stability control method provided by this invention acquires information on high-altitude cruising aircraft, constructs the aircraft flight time-space sequence, calculates the relative flight position interval and relative flight speed between aircraft, constructs an aircraft car-following model with time delay and bidirectional attention, and then constructs air traffic stability judgment conditions. It can determine the air traffic stability status and ensure air traffic stability by filtering out aircraft in unstable states for feedback control. Attached Figure Description

[0020] Figure 1 This is a flowchart of the air traffic stability control method provided in an embodiment of the present invention; Figure 2 This is a schematic diagram of the modeling process of the aircraft car-following model with time delay and bidirectional attention provided in an embodiment of the present invention; Figure 3 This is a schematic diagram of the stability region obtained based on the air traffic stability determination conditions provided in an embodiment of the present invention. Detailed Implementation

[0021] The technical solution of the present invention will be further described in detail below with reference to specific embodiments.

[0022] Embodiments of the present invention are described in detail below. Examples of these embodiments are illustrated in the accompanying drawings, wherein the same or similar reference numerals denote the same or similar elements or elements having the same or similar functions throughout. The embodiments described below with reference to the accompanying drawings are exemplary and are only used to explain the present invention, and should not be construed as limiting the present invention. Unless otherwise specified, embodiments of the present invention and the technical features thereof can be combined with each other.

[0023] This invention provides an air traffic stability control method, comprising: Construct the spatiotemporal sequence of aircraft flight based on high-altitude cruise aircraft information; Based on the spatiotemporal sequence of aircraft flight, calculate the relative flight position interval and relative flight speed between aircraft; Based on the relative flight position interval and relative flight speed between aircraft, an aircraft car-following model with time delay and bidirectional attention is constructed. Based on the aircraft car-following model, conditions for determining air traffic stability are constructed. Based on the criteria for determining air traffic stability, aircraft in unstable states are selected for feedback control to ensure air traffic stability.

[0024] The air traffic stability control method provided in this invention acquires information on high-altitude cruising aircraft, constructs the aircraft's flight time-space sequence, calculates the relative flight position interval and relative flight speed between aircraft, constructs an aircraft car-following model with time delay and bidirectional attention, and then constructs air traffic stability judgment conditions. This method can determine the air traffic stability status, and by filtering out aircraft in unstable states for feedback control, it ensures air traffic stability and provides a reference method for air traffic management.

[0025] Figure 1 This is a flowchart of an air traffic stability control method provided in an embodiment of the present invention. This flowchart only shows the logical sequence of the method in this embodiment. Provided there are no conflicts, different methods may be used. Figure 1 Complete the steps shown or described in the order indicated.

[0026] The air traffic stability control method provided in this embodiment of the invention can be applied to a terminal and can be executed by an air traffic stability control device. This device can be implemented by software and / or hardware and can be integrated into the terminal, such as any tablet computer or computer device with communication functions.

[0027] In one possible embodiment, constructing the aircraft's flight spatiotemporal sequence based on high-altitude cruise aircraft information specifically includes the following steps: Step 1: Extract the aircraft flight observation data stream based on high-altitude cruise aircraft information; In this embodiment, extracting the aircraft flight observation data stream based on high-altitude cruise aircraft information specifically includes the following steps: Step ①: Extract multi-source observation data based on high-altitude cruise aircraft information; Specifically, multi-source observation data includes flight plan messages, as well as route, segment, and aircraft information from flight information service systems, global navigation satellite systems, air-to-ground and air-to-air communication systems, secondary surveillance radar systems, and automatic dependent surveillance broadcast systems.

[0028] Step 2: Identify and associate the timestamps of multi-source observation data for the same aircraft according to the aircraft call sign, forming an aircraft flight observation data stream that includes the flight path, three-dimensional position, speed, and heading of each timestamp of the aircraft.

[0029] Step 2: Perform trajectory reconstruction, route projection, and queue sorting on the aircraft flight observation data stream to generate the aircraft flight spatiotemporal sequence.

[0030] In this embodiment, the process of reconstructing the trajectory, projecting the route, and sorting the queue of the aircraft flight observation data stream to generate the aircraft flight spatiotemporal sequence includes the following steps: Step ①: Resample the observation points for each aircraft to unify the time step of the aircraft's trajectory and obtain the continuous trajectory of the aircraft; Specifically, the aircraft flight observation data stream contains Same route at the time of observation All aircraft collection , for the Each aircraft has a discrete observation sequence This yields discrete observation sequences for all aircraft. This leads to the flight observation data streams of all aircraft. ,in, Indicates a set of routes. , … They represent the 1st, 2nd, ..., One aircraft, , … They represent the first The first, second, ..., of the aircraft A timestamp, , … They represent the first The aircraft , … Multi-source observation data.

[0031] The observation points for each aircraft are resampled, that is, the discrete observation sequence for each aircraft is sampled, with a uniform time step. , ,in, Indicates the start time of observation. Indicates the end time of observation, for any The continuous trajectory of the aircraft can be obtained by linear interpolation as follows: ; in, express Time of the first The continuous trajectory of an aircraft , Indicates the first The first aircraft , A timestamp, , Indicates the first The aircraft , Multi-source observation data.

[0032] Step 2: Select the position corresponding to each continuous trajectory at the time step, calculate the arc length projection of the continuous trajectory along the flight path centerline to obtain the heading coordinates of the aircraft, and sort the aircraft in ascending order of heading coordinates to form the aircraft flight time-space sequence.

[0033] Specifically, the spacetime sequence of aircraft flight is as follows: .

[0034] In one possible embodiment, routes W102, W45, W56, and W37 in the Beijing-Guangzhou corridor are selected. Flight schedule messages of high-altitude cruise aircraft along routes W102, W45, W56, and W37, as well as route, segment, and aircraft information from flight information service systems, global navigation satellite systems, air-to-ground and air-to-air communication systems, secondary surveillance radar systems, and automatic dependent surveillance broadcast systems are obtained. The multi-source observation timestamps of the same aircraft are identified and associated according to the aircraft call sign, forming a raw flight data stream containing the route where each timestamp of the aircraft is located, as well as its three-dimensional position, altitude, speed, and heading.

[0035] In this embodiment, the raw observation flight data streams of 1594 aircraft cruising at high altitudes along the W102, W45, W56, and W37 routes on May 7, 2025, obtained include... Same route at the time of observation All aircraft , for the Each aircraft has a discrete observation sequence .

[0036] The trajectory reconstruction, route projection, and queue sorting of 1,594 aircraft that cruised at high altitudes along routes W102, W45, W56, and W37 on May 7, 2025 were performed: the observation points of each aircraft were resampled to unify the time step and obtain continuous trajectories; the positions of each continuous trajectory at the time step were selected, the arc length projection along the route centerline was calculated to obtain the heading coordinates, and the trajectories were sorted in ascending order according to the route coordinates to form an ordered queue.

[0037] In one possible embodiment, calculating the relative flight position interval and relative flight speed between aircraft based on the aircraft flight time-space sequence specifically includes: constructing the aircraft flight time-space sequence into a formation sequence diagram, selecting adjacent aircraft based on the formation sequence diagram, and calculating the relative flight position interval and relative flight speed between adjacent aircraft.

[0038] In this embodiment, the queue sequence diagram is a spatial topology structure constructed with aircraft as nodes and the adjacency relationships between aircraft as directed edges.

[0039] In this embodiment, the spatial topology is defined as a directed graph. ,in, This represents the set of nodes in a directed graph, specifically the set of nodes consisting of aircraft. This represents the set of directed edges from one node to another, that is, the set of directed edges from one aircraft to another. express The Middle A directed edge.

[0040] Specifically, the relative flight position interval between adjacent aircraft is: ; in, express Time of the first The aircraft and the first The relative flight position interval between aircraft , They represent Time of the first , The flight positions of the aircraft.

[0041] The relative flight speeds between adjacent aircraft are: ; in, express Time of the first The aircraft and the first The relative flight speed between aircraft , They represent Time of the first , The flight speed of an aircraft , They represent Time of the first , The flight positions of the aircraft.

[0042] In one possible embodiment, such as Figure 2 As shown, based on the relative flight position interval and relative flight speed between aircraft, the specific steps for constructing an aircraft car-following model with time delay and bidirectional attention are as follows: Step 1: Calculate the desired positional separation between aircraft and use the relative flight positional separation between aircraft as the actual positional separation between aircraft; Specifically, the desired positional spacing between aircraft is: ; in, , They represent Time of the first The aircraft and the one in front of it One aircraft, rear of the first The desired positional spacing between aircraft , , They represent Time of the first The aircraft, the The first aircraft in front The aircraft, the The aircraft behind The desired flight location of an aircraft.

[0043] The actual positional spacing between aircraft is: ; in, , They represent Time of the first The aircraft and the first , The relative flight position interval between aircraft, i.e. Time of the first The aircraft and the one in front of it One aircraft, rear of the first The actual positional interval between aircraft , , They represent Time of the first , , The flight position of an aircraft, i.e. Time of the first The aircraft, the The first aircraft in front The aircraft, the The aircraft behind The actual flight position of the aircraft.

[0044] Step 2: Based on the expected and actual position intervals between aircraft, calculate the forward and backward following strengths of the aircraft; Specifically, the forward following strength of the aircraft is: ; in, express Time of the first The aircraft relative to the first one in front of it Forward following strength of an aircraft.

[0045] The rearward following strength of the aircraft is: ; in, express Time of the first The aircraft relative to the one behind it The rearward following strength of an aircraft.

[0046] Step 3: Based on the aircraft's forward and backward car-following strengths, calculate the forward and backward attention distributions of the pilots on the aircraft; Specifically, the distribution of a pilot's forward attention on an aircraft is as follows: ; in, express Time of the first The pilots on the aircraft were looking at the first... Forward attention value of an aircraft Indicates the first The total number of aircraft ahead of each aircraft.

[0047] The distribution of a pilot's backward attention on an aircraft is as follows: ; in, express Time of the first The pilots on the aircraft were looking at the first... The backward attention value of an aircraft Indicates the first The total number of aircraft behind each aircraft.

[0048] Step 4: Construct a two-way attention distribution model based on the forward and backward attention distributions of the pilots on the aircraft; The bidirectional attention distribution model is defined as the allocation of attention of the pilot on the current aircraft to the aircraft in front and behind based on the forward and backward car-following intensities of the current aircraft, thereby describing the degree of attention the pilot on the current aircraft pays to the aircraft in front and behind.

[0049] Step 5: Construct an air-to-air communication and response time delay model; Specifically, the air-to-air communication and response time delay model is as follows: ; in, This indicates air-to-air communication and response delay. This indicates the time required for air-to-air communication. Follows uniform distribution , , , They represent The minimum, maximum, and average values ​​of the value. This indicates the time it takes for the pilots on board the aircraft to respond to and manipulate the aircraft based on communication information. Follows uniform distribution , , , They represent The minimum, maximum, and average values.

[0050] Step 6: Based on the relative flight speed between aircraft, combine the two-way attention distribution model and the air-to-air communication and response time delay model to obtain an aircraft car-following model with time delay and two-way attention to describe the car-following characteristics of aircraft.

[0051] Specifically, the aircraft following model is as follows: ; in, express Time of the first The flight acceleration of an aircraft Indicates the first The sensitivity coefficient of an aircraft's response based on the status of the aircraft in front and behind it. , They represent Time of the first The aircraft and the first , The relative flight position interval between aircraft, i.e. Time of the first The aircraft and the one in front of it One aircraft, rear of the first The actual positional interval between aircraft express Time of the first The flight speed of an aircraft , They represent Time of the first The aircraft and the first , The relative flight speed between aircraft, i.e. Time of the first The aircraft and the one in front of it One aircraft, rear of the first The relative flight speed between aircraft , They represent Time of the first The aircraft and the one in front of it One aircraft, rear of the first Relative flight acceleration between aircraft Indicates the safe separation between aircraft. Represents the hyperbolic tangent function. Indicates according to and The calculated contents The The aircraft and the one in front of it The optimal flight speed for maintaining a safe distance between aircraft. Indicates according to and The calculated contents The The aircraft and the one behind it The optimal flight speed for maintaining a safe distance between aircraft. Indicates according to and The calculated first The aircraft and the one in front of it The optimal flight speed for maintaining a safe distance between aircraft. Indicates according to and The calculated first The aircraft and the one behind it The optimal flight speed for maintaining a safe distance between aircraft. , They represent , The first derivative, This represents the maximum value of the optimal flight speed.

[0052] Figure 2 middle, Indicates the region of stability. Indicates an unstable region. This indicates the total number of aircraft.

[0053] In one possible embodiment, the air traffic stability determination criterion is: ; in, Indicates the air traffic stability threshold. Indicates according to The first derivative of the determined optimal flight speed is the optimal flight acceleration.

[0054] In this embodiment, aircraft that meet the air traffic stability determination criteria are considered to be in a stable state, and are controlled to accelerate or maintain their current flight speed. Aircraft that do not meet the air traffic stability determination criteria are considered to be in an unstable state, and are controlled to automatically change their flight acceleration according to the aircraft following model to adjust the flight position interval between aircraft and ensure air traffic stability. The stability region obtained based on the air traffic stability determination criteria is as follows: Figure 3 As shown.

[0055] Figure 3 middle, This indicates the total number of aircraft.

[0056] This invention provides an air traffic stability control device, comprising: The aircraft flight spatiotemporal sequence construction module is used to construct the aircraft flight spatiotemporal sequence based on high-altitude cruise aircraft information; The aircraft relative relationship analysis module is used to calculate the relative flight position interval and relative flight speed between aircraft based on the aircraft flight time-space sequence; The car-following model building module is used to build an aircraft car-following model with time delay and bidirectional attention based on the relative flight position interval and relative flight speed between aircraft. The stability determination module is used to construct air traffic stability determination conditions based on the aircraft car-following model. The feedback control module is used to screen aircraft in unstable states based on air traffic stability determination conditions and perform feedback control to ensure air traffic stability.

[0057] The air traffic stability control device provided in this embodiment of the invention can execute the air traffic stability control method provided in this embodiment of the invention, and has the corresponding functional modules and beneficial effects of the method.

[0058] This invention provides a computer device, comprising: Storage medium used to store computer programs; A processor is used to execute computer programs to implement the air traffic stability control method provided in the embodiments of the present invention.

[0059] This invention provides a computer-readable storage medium storing a computer program thereon, which, when executed by a processor, implements the air traffic stability control method provided in this invention.

[0060] This embodiment provides a computer program product, including a computer program that, when executed by a processor, implements the air traffic stability control method provided in this embodiment of the invention.

[0061] Those skilled in the art will understand that embodiments of the present invention can be provided as methods, systems, or computer program products. Therefore, the present invention can take the form of a completely hardware embodiment, a completely software embodiment, or an embodiment combining software and hardware aspects. Furthermore, the present invention can take the form of a computer program product embodied on one or more computer-usable storage media (including, but not limited to, disk storage, CD-ROM, optical storage, etc.) containing computer-usable program code.

[0062] This invention is described with reference to flowchart illustrations and / or block diagrams of methods, apparatus (systems), and computer program products according to embodiments of the invention. It will be understood that each block of the flowchart illustrations and / or block diagrams, as well as combinations of blocks in the flowchart illustrations and / or block diagrams, can be implemented by computer program instructions. These computer program instructions can be provided to a processor of a general-purpose computer, special-purpose computer, embedded processor, or other programmable data processing apparatus to produce a machine, such that the instructions, which execute via the processor of the computer or other programmable data processing apparatus, generate instructions for implementing the flowchart. Figure 1 One or more processes and / or boxes Figure 1 A device that provides the functions specified in one or more boxes.

[0063] These computer program instructions may also be stored in a computer-readable storage medium that can direct a computer or other programmable data processing device to function in a particular manner, such that the instructions stored in the computer-readable storage medium produce an article of manufacture including instruction means, which are implemented in a process Figure 1 One or more processes and / or boxes Figure 1 The function specified in one or more boxes.

[0064] These computer program instructions may also be loaded onto a computer or other programmable data processing equipment to cause a series of operational steps to be performed on the computer or other programmable equipment to produce a computer-implemented process, thereby providing instructions that execute on the computer or other programmable equipment for implementing the process. Figure 1 One or more processes and / or boxes Figure 1 The steps of the function specified in one or more boxes.

[0065] The above are merely preferred embodiments of the present invention. It should be noted that those skilled in the art can make various improvements and modifications without departing from the technical principles of the present invention, and these improvements and modifications should also be considered within the scope of protection of the present invention.

Claims

1. A method for controlling air traffic stability, characterized in that, include: Construct the spatiotemporal sequence of aircraft flight based on high-altitude cruise aircraft information; Based on the spatiotemporal sequence of aircraft flight, calculate the relative flight position interval and relative flight speed between aircraft; Based on the relative flight position interval and relative flight speed between aircraft, an aircraft car-following model with time delay and bidirectional attention is constructed. Based on the aircraft car-following model, conditions for determining air traffic stability are constructed. Based on the criteria for determining air traffic stability, aircraft in unstable states are selected for feedback control to ensure air traffic stability.

2. The air traffic stability control method according to claim 1, characterized in that, Constructing the spatiotemporal sequence of an aircraft flight based on high-altitude cruise aircraft information includes: Based on high-altitude cruise aircraft information, extract aircraft flight observation data streams; The flight observation data stream of aircraft is used to reconstruct the trajectory, project the route, and sort the queue to generate the spatiotemporal sequence of aircraft flight.

3. The air traffic stability control method according to claim 2, characterized in that, Based on high-altitude cruise aircraft information, the extracted aircraft flight observation data stream includes: Extract multi-source observation data based on information from high-altitude cruise aircraft; The timestamps of multi-source observation data of the same aircraft are identified and associated according to the aircraft call sign, forming an aircraft flight observation data stream that includes the route, three-dimensional position, speed and heading of each timestamp of the aircraft; The multi-source observation data includes flight plan messages, as well as route, segment, and aircraft information from flight information service systems, global navigation satellite systems, air-to-ground and air-to-air communication systems, secondary surveillance radar systems, and automatic dependent surveillance broadcast systems.

4. The air traffic stability control method according to claim 2, characterized in that, The flight observation data stream is subjected to trajectory reconstruction, route projection, and queue sorting to generate the aircraft's spatiotemporal flight sequence, including: The observation points for each aircraft are resampled, and the time step of the aircraft trajectory is standardized to obtain the continuous trajectory of the aircraft. The positions corresponding to each continuous trajectory at the time step are selected, and the arc length projection of the continuous trajectory is calculated along the centerline of the flight path to obtain the heading coordinates of the aircraft. The aircraft are then sorted in ascending order according to their heading coordinates to form a spatiotemporal sequence of aircraft flight.

5. The air traffic stability control method according to claim 1, characterized in that, Based on the spatiotemporal sequence of aircraft flight, the calculation of the relative flight position interval and relative flight speed between aircraft includes: The flight time-space sequence of aircraft is constructed into a formation sequence diagram. Based on the formation sequence diagram, adjacent aircraft are selected, and the relative flight position interval and relative flight speed between adjacent aircraft are calculated. The queue sequence diagram is a spatial topology constructed with aircraft as nodes and the adjacency relationships between aircraft as directed edges.

6. The air traffic stability control method according to claim 5, characterized in that, The relative flight position interval between adjacent aircraft is: ； in, express Time of the first The aircraft and the first The relative flight position interval between aircraft , They represent Time of the first , The flight positions of the aircraft; The relative flight speeds between adjacent aircraft are: ； in, express Time of the first The aircraft and the first The relative flight speed between aircraft , They represent Time of the first , The flight speed of an aircraft , They represent Time of the first , The flight positions of the aircraft.

7. The air traffic stability control method according to claim 1, characterized in that, Based on the relative flight position intervals and relative flight speeds between aircraft, an aircraft car-following model with time delay and bidirectional attention is constructed, including: Calculate the desired positional spacing between aircraft and use the relative flight positional spacing between aircraft as the actual positional spacing between aircraft; Based on the expected and actual positional intervals between aircraft, calculate the forward and backward following strengths of the aircraft. Based on the forward and backward car-following strengths of the aircraft, calculate the forward and backward attention distributions of the pilots on the aircraft. A two-way attention distribution model is constructed from the forward and backward attention distributions of pilots on aircraft; Construct an air-to-air communication and response time delay model; Based on the relative flight speed between aircraft, a two-way attention distribution model and an air-to-air communication and response time delay model are combined to obtain an aircraft car-following model with time delay and two-way attention to describe the car-following characteristics of aircraft.

8. The air traffic stability control method according to claim 7, characterized in that, The desired spacing between aircraft is: ； in, , They represent Time of the first The aircraft and the one in front of it One aircraft, rear of the first The desired positional spacing between aircraft , , They represent Time of the first The aircraft, the The first aircraft in front The aircraft, the The aircraft behind The expected flight position of an aircraft; The actual positional spacing between aircraft is: ； in, , They represent Time of the first The aircraft and the first , The relative flight position interval between aircraft, i.e. Time of the first The aircraft and the one in front of it One aircraft, rear of the first The actual positional interval between aircraft , , They represent Time of the first , , The flight position of an aircraft, i.e. Time of the first The aircraft, the The first aircraft in front The aircraft, the The aircraft behind The actual flight position of the aircraft; The forward following strength of the aircraft is: ； in, express Time of the first The aircraft relative to the first one in front of it Forward following strength of an aircraft; The rearward following strength of the aircraft is: ； in, express Time of the first The aircraft relative to the one behind it The rearward following strength of an aircraft; The forward attention distribution of a pilot on an aircraft is as follows: ； in, express Time of the first The pilots on the aircraft were looking at the first... Forward attention value of an aircraft Indicates the first The total number of aircraft ahead of each aircraft; The distribution of a pilot's backward attention on an aircraft is as follows: ； in, express Time of the first The pilots on the aircraft were looking at the first... The backward attention value of an aircraft Indicates the first The total number of aircraft behind each aircraft; The air-to-air communication and response time delay model is as follows: ； in, This indicates air-to-air communication and response delay. This indicates the time required for air-to-air communication. Follows uniform distribution , , , They represent The minimum, maximum, and average values ​​of the value. This indicates the time it takes for the pilots on board the aircraft to respond to and manipulate the aircraft based on communication information. Follows uniform distribution , , , They represent Minimum, maximum, and average values; The aircraft following model is: ； in, express Time of the first The flight acceleration of an aircraft Indicates the first The sensitivity coefficient of an aircraft's response based on the status of the aircraft in front and behind it. , They represent Time of the first The aircraft and the first , The relative flight position interval between aircraft, i.e. Time of the first The aircraft and the one in front of it One aircraft, rear of the first The actual positional interval between aircraft express Time of the first The flight speed of an aircraft , They represent Time of the first The aircraft and the first , The relative flight speed between aircraft, i.e. Time of the first The aircraft and the one in front of it One aircraft, rear of the first The relative flight speed between aircraft , They represent Time of the first The aircraft and the one in front of it One aircraft, rear of the first Relative flight acceleration between aircraft Indicates the safe separation between aircraft. Represents the hyperbolic tangent function. Indicates according to and The calculated contents The The aircraft and the one in front of it The optimal flight speed for maintaining a safe distance between aircraft. Indicates according to and The calculated contents The The aircraft and the one behind it The optimal flight speed for maintaining a safe distance between aircraft. Indicates according to and The calculated first The aircraft and the one in front of it The optimal flight speed for maintaining a safe distance between aircraft. Indicates according to and The calculated first The aircraft and the one behind it The optimal flight speed for maintaining a safe distance between aircraft. , They represent , The first derivative, This represents the maximum value of the optimal flight speed.

9. The air traffic stability control method according to claim 8, characterized in that, The criteria for determining air traffic stability are: ； in, Indicates the air traffic stability threshold. Indicates according to The first derivative of the determined optimal flight speed is the optimal flight acceleration.

10. The air traffic stability control method according to claim 1, characterized in that, Based on air traffic stability assessment criteria, aircraft in unstable states are selected for feedback control to ensure air traffic stability, including: Aircraft that do not meet the criteria for determining air traffic stability are considered to be in an unstable state. The system controls aircraft in unstable states to automatically change their flight acceleration according to the aircraft following model in order to adjust the flight position interval between aircraft and ensure air traffic stability. Among them, aircraft that meet the air traffic stability determination conditions are considered to be aircraft in a stable state; Control the aircraft to accelerate or maintain its current flight speed in a stable state.

Images