CN117421542A - Method, system and device for calculating sliding speed of aircraft in and out of service area - Google Patents

Method, system and device for calculating sliding speed of aircraft in and out of service area Download PDF

Info

Publication number
CN117421542A
CN117421542A CN202311744250.1A CN202311744250A CN117421542A CN 117421542 A CN117421542 A CN 117421542A CN 202311744250 A CN202311744250 A CN 202311744250A CN 117421542 A CN117421542 A CN 117421542A
Authority
CN
China
Prior art keywords
aircraft
stand
speed
distance
time
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
CN202311744250.1A
Other languages
Chinese (zh)
Other versions
CN117421542B (en
Inventor
胡婕
王志锋
刘松
赵辰旭
李明
王燕
杨实
王治宇
黄智豪
郝德月
吴熙泰
钟雨泉
吴坤永
郭浩彬
苏剑飞
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Zhuhai Xiangyi Aviation Technology Co Ltd
Original Assignee
Zhuhai Xiangyi Aviation Technology Co Ltd
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Zhuhai Xiangyi Aviation Technology Co Ltd filed Critical Zhuhai Xiangyi Aviation Technology Co Ltd
Priority to CN202311744250.1A priority Critical patent/CN117421542B/en
Publication of CN117421542A publication Critical patent/CN117421542A/en
Application granted granted Critical
Publication of CN117421542B publication Critical patent/CN117421542B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F18/00—Pattern recognition
    • G06F18/10—Pre-processing; Data cleansing
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F18/00—Pattern recognition
    • G06F18/20—Analysing
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F18/00—Pattern recognition
    • G06F18/20—Analysing
    • G06F18/24—Classification techniques
    • G06F18/243—Classification techniques relating to the number of classes
    • G06F18/2433—Single-class perspective, e.g. one-against-all classification; Novelty detection; Outlier detection
    • G—PHYSICS
    • G07—CHECKING-DEVICES
    • G07C—TIME OR ATTENDANCE REGISTERS; REGISTERING OR INDICATING THE WORKING OF MACHINES; GENERATING RANDOM NUMBERS; VOTING OR LOTTERY APPARATUS; ARRANGEMENTS, SYSTEMS OR APPARATUS FOR CHECKING NOT PROVIDED FOR ELSEWHERE
    • G07C5/00—Registering or indicating the working of vehicles
    • G07C5/08—Registering or indicating performance data other than driving, working, idle, or waiting time, with or without registering driving, working, idle or waiting time
    • G07C5/0841—Registering performance data

Landscapes

  • Engineering & Computer Science (AREA)
  • Theoretical Computer Science (AREA)
  • Data Mining & Analysis (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Bioinformatics & Computational Biology (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Bioinformatics & Cheminformatics (AREA)
  • Evolutionary Biology (AREA)
  • Evolutionary Computation (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • General Engineering & Computer Science (AREA)
  • Artificial Intelligence (AREA)
  • Traffic Control Systems (AREA)

Abstract

The invention belongs to the field of data analysis, in particular relates to a method, a system and a device for calculating the sliding speed of an airplane in a stop-entering area, and aims to solve the problem that the evaluation dimension of the airplane in the stop-entering process is single in the prior art, and the defect that the actual safe running condition of the airplane in the stop-entering stage cannot be accurately evaluated exists. The invention comprises the following steps: acquiring an original operation parameter data set of a target aircraft in a landing taxiing stage; removing abnormal values in the original operation parameter data set, and performing standardization processing to obtain standardized operation parameters, so as to calculate the moment when the aircraft arrives at the stand, and further calculate the speed of each moment in the process of entering the stand area of the aircraft. The invention can calculate and monitor the sliding speed of the airplane at the final airplane entering and stopping stage to obtain quantitative evaluation so as to adjust the driving habit of the pilot and eliminate potential safety hazards.

Description

Method, system and device for calculating sliding speed of aircraft in and out of service area
Technical Field
The invention belongs to the field of data analysis, and particularly relates to a method, a system and a device for calculating the sliding speed of an aircraft in and out of a parking area.
Background
In recent years, the development of aeronautical technology is rapid, and an airplane becomes one of ideal transportation means in daily travel of people. The taxiing speed of an aircraft during the entering of a stand is a critical parameter in the operation of the aircraft, and has important significance for the safety assessment of whether the aircraft can accurately stop at a specified stand. At present, the safety research on the taxi speed of the aircraft in the airplane station process is less, even if the track of the aircraft in the airplane station process is drawn simply through longitude and latitude parameters, whether the situation of deviating from a taxi guiding line exists before the aircraft is in place is judged simply through the track, however, the taxi speed of the aircraft in the airplane station process is ignored to be an important factor affecting the safety of the flying operation, and if the speed of the aircraft in the airplane station process is too high or the speed change rate is too high, unexpected stopping of the aircraft or the aircraft rushing through a stopping line to collide with a corridor or ground facilities can be caused.
There is a need for an objective and reliable method for accurately counting the speed of the aircraft during the last aircraft stand stage.
Disclosure of Invention
In order to solve the above problems in the prior art, that is, the problem that the dimension of the airplane estimated in the process of entering and stopping a station is single in the prior art, and the defect that the actual safe running condition of the airplane in the stage of entering and stopping the station cannot be accurately estimated, so that the problem that potential safety hazards exist is solved, the invention provides a method for calculating the sliding speed of the airplane in the area of entering and stopping the station, which comprises the following steps:
step S100, acquiring an original operation parameter data set of a target aircraft in a landing taxiing stage;
step S200, eliminating abnormal values in the original operation parameter data set and performing standardization processing to obtain a standardized operation parameter data set;
step S300, calculating the moment of the aircraft to the stand according to the operation parameters in the standardized operation parameter data set；
The step S300 specifically includes:
step S310, according to the parking brake parameter, the space gate starts to traverse forward from the point where AIR changes to GROUND, and finds that the parking brake parameter changes from mark 0The time of the mark 1 is recorded as the set stay timeThe method comprises the steps of carrying out a first treatment on the surface of the When the stay braking parameter is 0, the brake is released, and when the stay braking parameter is 1, the stay braking is set;
step S320, setting the stay time according to the longitude and latitude of the airport standCalculating the arc length between the coordinate point of the aircraft and the coordinate point of the stand>：
；
；
；
；
Wherein,is the longitude of the assigned airport stand, +.>Is the latitude of the designated airport stand; />Is to set the stay time +.>Longitude, latitude of current aircraft>Is to set the stay time +.>Current aircraft latitude,/->Indicating deflection angle +.>Is the difference of dimensions, +.>Is the difference in longitude;
judging the arc lengthWith a preset accurate berth threshold value, if the arc length is + ->If the preset accurate berth threshold value is smaller than or equal to the preset accurate berth threshold value, judging the current momenttThe aircraft is correctly stopped at the stand, and the current time is output +.>As aircraft to stand time->The method comprises the steps of carrying out a first treatment on the surface of the If the arc length is +>If the data is larger than a preset accurate berth threshold value, judging that the data acquired by the current aircraft cannot accurately judge that the aircraft stops on a stand;
step S330, determining the time of arrival at the standIs>；
Step S340, if the stand time is setIs>0, go to step S400; if the stand is at time->Is the taxiing speed V of the aircraft t If the current value is not 0, the point moment of the stand is recorded>Is>Reverse traversal would be within 30 seconds of the corresponding taxi speed +.>All minus the first aircraft speed magnitudeEnsure stand time +.>The taxiing speed is 0, and an airplane taxiing speed sequence is obtained;
step S400, according to the moment when the aircraft arrives at the standThe operation parameter data are standardized, and the speed of each moment in the process of entering the stand area of the aircraft is calculated;
the step S400 specifically includes:
step S410, obtaining standardized operation parameter data, including:
Time[n]: distance standnThe moment of the rice is initiallyNone；
LIST_GS[n]: distance standnThe taxiing speed of the rice aircraft is initiallyNone；
GS[i]: second at the present timeiThe taxiing speed of the aircraft is initiallyNone；
: cumulative distance, representing the distance moved by the aircraft per second, initially 0;
: the target observation distance from the preset airplane to the stand is set to be 1 meter;
step S420, from the aircraft to the standStart reverse cycle traverse of coasting speedGSThe iteration mode is as follows:
；
wherein the method comprises the steps ofFor the current time secondiAnd the last time secondi-1, i is the initial value fromCurrent time of start second, 1 sea/hour = 1852/3600 m/s;
step S430, the steps are performedAccumulated to accumulated distance +.>In (I)>；
Step S440, if the accumulated distanceLess than->Advancing the time second i by 1 second, and repeating the steps S420 and S430; when accumulating distance->First time greater than or equal to the target viewing distance +.>Then it is stated that +.>Distance from rice, i.e.)>Jump to S450;
step S450, ifGS[i]=GS[i-1]Distance to stand target observation distanceThe sliding speed +.>Is thatGS[i-1]；
If it isGS[i]≠GS[i-1]According to the accumulated distance、/>And i-1 and the current second i find the viewing distance from the stand target>Speed at rice:
a=GS[i-1]-GS[i]；
b=GS[i]*2；
；
；
；
by->Rounding to take a decimal place as the observation distance from the shutdown target>A sliding speed at the rice;
step S460, object observation distanceIncreasing the viewing distance of the new target by 1 meterIf the new target is observed at a distance +>Marking the new target observation distance as the target observation distance, repeating the steps S420-S450 until +.>；
Step S470, output，LIST_The value in GS sequentially corresponds to the ground speed of 1 to 10 meters from the stand as the speed of each moment in the process of entering the stand area of the airplane.
In some preferred embodiments, the step S200 specifically includes:
step S210, eliminating or replacing abnormal operation parameter data based on the original operation parameter data set to obtain an operation parameter data set with abnormal values eliminated;
the abnormal operation parameter data comprise frame dropping parameters and operation parameter data with parameter records exceeding a preset abnormal threshold;
step S220, data cleaning and data standardization are carried out on the operation parameter data set with abnormal values removed, and a standardized operation parameter data set is obtained.
In some preferred embodiments, the standardized set of operating parameter data specifically includes:
the aircraft stand longitude and latitude, the aircraft longitude and latitude, the sliding speed, the stay braking parameter and the magnetic heading parameter;
the longitude and latitude of the parking apron and the longitude and latitude of the airplane use a unified coordinate system;
the stay brake parameter includes a brake-on flag of 1 and a brake-off flag of 0.
The method further comprises step S500:
and comparing the speed of the aircraft at each moment in the aircraft entering and exiting station area with the interpolation of the standard value to obtain the safe operation evaluation result of the aircraft entering and exiting station stage.
In another aspect of the invention, a system for calculating the taxiing speed of an aircraft in and out of a parking area is provided, the system comprising:
the original parameter acquisition module is configured to acquire an original operation parameter data set of the target aircraft in a landing taxiing stage;
the preprocessing module is configured to remove abnormal values in the original operation parameter data set and perform standardized processing to obtain a standardized operation parameter data set;
the stand time calculation module is configured to calculate the time when the aircraft arrives at the stand according to the operation parameters in the standardized operation parameter data set；
The method specifically comprises the following steps:
according to the stay braking parameters, the point of the AIR-GROUND gate changing from AIR to GROUND starts to go forward, and the stay braking parameters are found outThe time when the number changes from mark 0 to mark 1 is recorded as the set stay timeThe method comprises the steps of carrying out a first treatment on the surface of the When the stay braking parameter is 0, the brake is released, and when the stay braking parameter is 1, the stay braking is set;
according to the longitude and latitude of the airport stand and the setting of the stay timeCalculating the arc length between the coordinate point of the aircraft and the coordinate point of the stand>：
；
；
；
；
Wherein,is the longitude of the assigned airport stand, +.>Is the latitude of the designated airport stand; />Is to set the stay time +.>Longitude, latitude of current aircraft>Is to set the stay time +.>Current aircraft latitude,/->Indicating deflection angle +.>Is the difference of dimensions, +.>Is the difference in longitude;
judging the arc lengthWith a preset accurate berth threshold value, if the arc length is + ->If the preset accurate berth threshold value is smaller than or equal to the preset accurate berth threshold value, judging the current momenttThe aircraft is correctly stopped at the stand, and the current time is output +.>As aircraft to stand time->The method comprises the steps of carrying out a first treatment on the surface of the If the arc length is +>If the data is larger than a preset accurate berth threshold value, judging that the data acquired by the current aircraft cannot accurately judge that the aircraft stops on a stand;
determining the arrival time of the machine at the standIs>；
If the stand time is setIs>0, carrying out a speed change condition acquisition module; if the stand is at time->Is the taxiing speed V of the aircraft t If the current value is not 0, the point moment of the stand is recorded>Is>Reverse traversal would be within 30 seconds of the corresponding taxi speed +.>The first aircraft speed magnitude is subtracted>Ensure stand time +.>The taxiing speed is 0, and an airplane taxiing speed sequence is obtained;
a speed change condition acquisition module configured to, according to the aircraft arrival timeThe operation parameter data are standardized, and the speed of each moment in the process of entering the stand area of the aircraft is calculated;
an initializing unit for obtaining standardized operation parameter data, including:
Time[n]: distance standnThe moment of the rice is initiallyNone；
LIST_GS[n]: distance standnThe taxiing speed of the rice aircraft is initiallyNone；
GS[i]: second at the present timeiThe taxiing speed of the aircraft at the time,initial asNone；
: cumulative distance, representing the distance moved by the aircraft per second, initially 0;
: the target observation distance from the preset airplane to the stand is set to be 1 meter;
taxiing speed iteration unit for time from airplane to standStart reverse cycle traverse of coasting speedGSThe iteration mode is as follows:
；
wherein the method comprises the steps ofFor the current time secondiAnd the last time secondi-1, i is the initial value fromCurrent time of start second, 1 sea/hour = 1852/3600 m/s;
an accumulation unit for accumulating the aboveAccumulated to accumulated distance +.>In,；
an accumulated distance iteration unit for calculating the accumulated distanceLess than->Advancing the time second i for 1 second, and repeating the sliding speed iteration unit and the accumulation unit; when accumulating distance->First greater than or equal to the target viewing distanceThen it is stated that +.>Distance from rice, i.e.)>Jumping to a sliding speed calculating unit;
a sliding speed calculating unit, ifGS[i]=GS[i-1]Distance to stand target observation distanceThe sliding speed +.>Is thatGS[i-1]；
If it isGS[i]≠GS[i-1]According to the accumulated distance、/>And i-1 and the current second i find the viewing distance from the stand target>Speed at rice:
a=GS[i-1]-GS[i]；
b=GS[i]*2；
；
；
；
by->Rounding to take a decimal place as the observation distance from the shutdown target>A sliding speed at the rice;
target observation distance updating unit for updating target observation distanceIncreasing the viewing distance of the new target by 1 meterIf the new target is observed at a distance +>Marking the new target observation distance as the target observation distance, and repeating the sliding speed iteration unit to the sliding speed calculation unit until；
Data output unit for outputting，LIST_The value in GS sequentially corresponds to the ground speed of 1 to 10 meters from the stand as the speed of each moment in the process of entering the stand area of the airplane.
In a third aspect of the present invention, an electronic device is provided, including:
at least one processor; and
a memory communicatively coupled to at least one of the processors; wherein,
the memory stores instructions executable by the processor for execution by the processor to implement the method of calculating aircraft taxiing speed into a stopping area described above.
In a fourth aspect of the present invention, a computer readable storage medium is provided, the computer readable storage medium storing computer instructions for execution by the computer to implement the above-described method of calculating aircraft taxiing speed in a stopping zone.
The invention has the beneficial effects that:
(1) The method for calculating the taxi speed of the airplane in the stop zone can calculate and monitor the taxi speed of the airplane in the final stop station stage, can obtain quantitative evaluation and provide analysis basis so as to adjust the driving habit of a pilot and eliminate potential safety hazards.
(2) And calculating the speed of 10 meters before the aircraft enters the stand and 1 meter at intervals based on the data set. On the one hand, the method can calculate and monitor the sliding speed of the aircraft in the final aircraft entering and stopping stage, on the other hand, the pilot skill shortboard is quantized by carrying out the processing on each execution flight of the pilot, the speed and throttle control habit of the aircraft in the aircraft entering and stopping stage are mastered, and the personnel with poor speed and throttle control are found out to carry out technical improvement. On the other hand, the method can realize accurate assessment of the safe running condition of the aircraft in the station entering stage and provide basis for aviation security event investigation.
Drawings
Other features, objects and advantages of the present application will become more apparent upon reading of the detailed description of non-limiting embodiments, made with reference to the following drawings, in which:
FIG. 1 is a flow chart of a method of calculating taxi speeds in an aircraft on-off region according to an embodiment of the invention;
FIG. 2 is a schematic illustration of a calculation flow of aircraft to stand time in an embodiment of the present invention;
FIG. 3 is a schematic view of a calculation flow of speed at each moment in the process of entering a stand area of an aircraft according to an embodiment of the invention;
fig. 4 is a schematic diagram of the calculation of arc length in an embodiment of the invention.
Detailed Description
The present application is described in further detail below with reference to the drawings and examples. It is to be understood that the specific embodiments described herein are merely illustrative of the invention and are not limiting of the invention. It should be noted that, for convenience of description, only the portions related to the present invention are shown in the drawings.
It should be noted that, in the case of no conflict, the embodiments and features in the embodiments may be combined with each other. The present application will be described in detail below with reference to the accompanying drawings in conjunction with embodiments.
In order to more clearly describe the method for calculating the taxiing speed of an aircraft in a stopping area according to the present invention, each step of the embodiment of the present invention will be described in detail with reference to fig. 1.
The method for calculating the taxiing speed of the aircraft in the stopping area according to the first embodiment of the invention comprises the following steps of step S100 to step S400:
step S100, acquiring an original operation parameter data set of a target aircraft in a landing taxiing stage; the present invention is implemented based on QAR data and airport data.
Step S200, eliminating abnormal values in the original operation parameter data set and performing standardization processing to obtain a standardized operation parameter data set;
in this embodiment, the step S200 specifically includes:
step S210, eliminating or replacing abnormal operation parameter data based on the original operation parameter data set to obtain an operation parameter data set with abnormal values eliminated;
the abnormal operation parameter data comprise frame dropping parameters and operation parameter data with parameter records exceeding a preset abnormal threshold;
step S220, data cleaning and data standardization are carried out on the operation parameter data set with abnormal values removed, and a standardized operation parameter data set is obtained. In this embodiment, the requirement on the latitude and longitude data is high, and the accuracy and smoothness of the latitude and longitude data must be ensured, so that the steps including interpolation, fitting and smoothing processing are required in standardization.
In this embodiment, the standardized operation parameter data set specifically includes:
the aircraft stand longitude and latitude, the aircraft longitude and latitude, the sliding speed, the stay braking parameter and the magnetic heading parameter;
the longitude and latitude of the parking apron and the longitude and latitude of the airplane use a unified coordinate system; the present embodiment uses the WGS84 ellipsoid coordinate system of the earth;
wherein, the standardized parameter process comprises the following steps: taking the model a320 as an example, in the data set of the invention, a320 has a total of 10410, 10411, 10412, 103251, 10413, 320020, 320232 and other 7 code libraries, and the same logical parameters of each code library are mapped into the same standard parameter name, namely the standardized parameter process, and the standardized parameters such as longitude and latitude, speed, stop brake and the like of the aircraft are shown in table 1.
Table 1 standardizes aircraft longitude and latitude, speed, and stopping brake
Step S300, calculating the moment of the aircraft to the stand according to the operation parameters in the standardized operation parameter data set；
In this embodiment, the calculation process of step S300 is shown in fig. 2, and specifically includes:
step S310, according to the stop brake parameter, starting forward traversal from the flight data landing point (i.e. the AIR gate is changed from AIR to GROUND), and recording the time when the stop brake parameter is changed from tag 0 to tag 1 as the set stop timeThe method comprises the steps of carrying out a first treatment on the surface of the When the stay braking parameter is 0, the brake is released, and when the stay braking parameter is 1, the stay braking is set;
step S320, setting the stay time according to the longitude and latitude of the airport standCalculating the arc length between the coordinate point of the aircraft and the coordinate point of the stand>：
；
；
；
；
Wherein,is the longitude of the assigned airport stand, +.>Is the latitude of the designated airport stand; />Is to set the stay time +.>Longitude, latitude of current aircraft>Is to set the stay time +.>Current aircraft latitude,/->Indicating deflection angle +.>Is the difference of dimensions, +.>Is the difference in longitude;
judging the arc lengthWith a preset accurate berth threshold value, if the arc length is + ->If the current time is smaller than or equal to the preset accurate berth threshold value, judging the current time +.>The aircraft is correctly stopped at the stand, and the current time is output +.>As aircraft to stand time->The method comprises the steps of carrying out a first treatment on the surface of the If the arc length is +>If the data is larger than a preset accurate berth threshold value, judging that the data acquired by the current aircraft cannot accurately judge that the aircraft stops on a stand;
step S330, determining the time of arrival at the standIs>；
Step S340, if setStand timeIs>0, go to step S340; if the stand is at time->Is>If the current value is not 0, the point moment of the stand is recorded>Is>Reverse traversal would be within 30 seconds of the corresponding taxi speed +.>All minus the first aircraft speed magnitudeEnsure stand time +.>And the taxiing speed is 0, and an airplane taxiing speed sequence is obtained.
If the arc lengthIf the arc length is smaller than or equal to the preset arc length threshold value, determining that the current aircraft is stopped at the aircraft stand, wherein the arc length is represented as shown in fig. 4, O in fig. 4 represents the sphere center of the earth, AB is the arc length, and outputting the current momenttAs aircraft to stand time->The method comprises the steps of carrying out a first treatment on the surface of the If the arc length is +>And if the data is larger than the preset arc length threshold value, judging that the current aircraft data cannot accurately judge that the aircraft stops at the aircraft stand. In this embodiment, the arc length threshold may be selected to be 0.00011, that is, the arc length is smaller than 0.00011, and the current aircraft is considered to stop at the stand, if the arc length is larger than 0.00012, the data collected by the current aircraft is considered to be unable to accurately support the judgment of the stop of the aircraft at the stand, and if the data is unable to accurately support the judgment of the stop of the aircraft at the stand, the corresponding data is marked as problem data and is returned to the processing center.
Step S400, according to the moment when the aircraft arrives at the standThe operation parameter data are standardized, and the speed of each moment in the process of entering the stand area of the aircraft is calculated;
in this embodiment, the calculation process of step S400 is shown in fig. 3, and specifically includes:
step S410, obtaining the standardized operation parameter data, including:
Time[n]: distance standnThe moment of the rice is initiallyNone；
LIST_GS[n]: distance standnThe taxiing speed of the rice aircraft is initiallyNone；
GS[i]: second at the present timeiThe taxiing speed of the aircraft is initiallyNone；
: cumulative distance, representing the distance moved by the aircraft per second, initially 0;
: the target observation distance from the preset airplane to the stand is 1 meter initially;
step S420, from the aircraft to the standStart reverse cycle traverse of coasting speedGSThe iteration mode is as follows:
；
；
wherein the method comprises the steps ofFor the current time secondiAnd the last time secondi-1, i is the initial value fromCurrent time of start second, 1 sea/hour = 1852/3600 m/s;
step S430, the steps are performedAccumulated to accumulated distance +.>In,；
step S440, if the accumulated distanceLess than->Advancing the time second i by 1 second, and repeating the steps S420 and S430; when accumulating distance->First time greater than or equal to the target viewing distance +.>Then it is stated that +.>Distance from rice, i.e.)>Jump to S450;
step S450, accumulating distanceSet to the target observation distance +.>Rice, if GS [ i ]]=GS[i-1]Then distance to stand target observation distance +.>The sliding speed +.>；
If it isGS[i]≠GS[i-1]According to the accumulated distanceThe taxi speed at time i-1 and the taxi speed at current time second i find the viewing distance from the stand target +.>Speed at rice:
a=GS[i-1]-GS[i]；
b=GS[i]*2；
；
；
；
by->Rounding to take a decimal place as the observation distance from the shutdown target>A sliding speed at the rice;
step S460, object observation distanceIncreasing the viewing distance of the new target by 1 meterIf the new target is observed at a distance +>Marking the new target observation distance as target observation distance, repeating steps S410-S450 until +.>；
Step S470, output and， LIST_the value in GS sequentially corresponds to the ground speed of 1 to 10 meters from the stand as the speed of each moment in the process of entering the stand area of the airplane. In this embodiment, the method further includes: and S500, carrying out early warning according to a set early warning rule based on the speed of each moment in the process of entering the airplane station area.
In this embodiment, the step S500,
and comparing the speed of the aircraft at each moment in the aircraft entering and exiting station area with the interpolation of the standard value to obtain the safe operation evaluation result of the aircraft entering and exiting station stage.
The invention can obtain the speed of 10 meters before the aircraft enters the stand and the speed of 1 meter at intervals. On the one hand, the method can calculate and monitor the sliding speed of the aircraft in the final aircraft entering and stopping stage, on the other hand, the pilot skill shortboard is quantized by carrying out the processing on each execution flight of the pilot, the speed and throttle control habit of the aircraft in the aircraft entering and stopping stage are mastered, and the personnel with poor speed and throttle control are found out to carry out technical improvement. On the other hand, the method can realize accurate assessment of the safe running condition of the aircraft in the station entering stage and provide basis for aviation security event investigation.
Although the steps are described in the above-described sequential order in the above-described embodiments, it will be appreciated by those skilled in the art that in order to achieve the effects of the present embodiments, the steps need not be performed in such order, and may be performed simultaneously (in parallel) or in reverse order, and such simple variations are within the scope of the present invention.
A system for calculating a taxiing speed of an aircraft in a stopping zone according to a second embodiment of the invention, the system comprising:
the original parameter acquisition module is configured to acquire an original operation parameter data set of the target aircraft in a landing taxiing stage;
the preprocessing module is configured to remove abnormal values in the original operation parameter data set and perform standardized processing to obtain a standardized operation parameter data set;
the stand time calculation module is configured to calculate the time when the aircraft arrives at the stand according to the operation parameters in the standardized operation parameter data set；
A speed change condition acquisition module configured to, according to the aircraft arrival timeThe operation parameter data are standardized, and the speed of each moment in the process of entering the stand area of the aircraft is calculated;
an initializing unit for obtaining standardized operation parameter data, including:
Time[n]: distance standnAt the moment when the rice is positioned,initial asNone；
LIST_GS[n]: distance standnThe taxiing speed of the rice aircraft is initiallyNone；
GS[i]: second at the present timeiThe taxiing speed of the aircraft is initiallyNone；
: cumulative distance, representing the distance moved by the aircraft per second, initially 0;
: the target observation distance from the preset airplane to the stand is set to be 1 meter;
taxiing speed iteration unit for time from airplane to standStart reverse cycle traverse of coasting speedGSThe iteration mode is as follows:
；
wherein the method comprises the steps ofFor the current time secondiAnd the last time secondi-1, i is the initial value from +.>Current time of start second, 1 sea/hour = 1852/3600 m/s;
an accumulation unit for accumulating the aboveAccumulated to accumulated distance +.>In,；
an accumulated distance iteration unit for calculating the accumulated distanceLess than->Advancing the time second i for 1 second, and repeating the sliding speed iteration unit and the accumulation unit; when accumulating distance->First greater than or equal to the target viewing distanceThen it is stated that +.>Distance from rice, i.e.)>Jumping to a sliding speed calculating unit;
a sliding speed calculating unit, ifGS[i]=GS[i-1]Distance to stand target observation distanceThe sliding speed +.>Is thatGS[i-1]；
If it isGS[i]≠GS[i-1]According to the accumulated distance、/>And i-1 and the current second i find the viewing distance from the stand target>Speed at rice：
a=GS[i-1]-GS[i]；
b=GS[i]*2；
；
；
；
By->Rounding to take a decimal place as the observation distance from the shutdown target>A sliding speed at the rice;
target observation distance updating unit for updating target observation distanceIncreasing the viewing distance of the new target by 1 meterIf the new target is observed at a distance +>Marking the new target observation distance as the target observation distance, and repeating the sliding speed iteration unit to the sliding speed calculation unit until；
Data output unit for outputting，LIST_The value in GS sequentially corresponds to the ground speed of 1 to 10 meters from the stand as the speed of each moment in the process of entering the stand area of the airplane.
It will be clear to those skilled in the art that, for convenience and brevity of description, the specific working process of the system described above and the related description may refer to the corresponding process in the foregoing method embodiment, which is not repeated here.
It should be noted that, in the system for calculating the taxiing speed of the aircraft in the stop zone provided in the above embodiment, only the division of the above functional modules is used as an example, in practical application, the above functional allocation may be performed by different functional modules according to needs, that is, the modules or steps in the embodiment of the present invention are further decomposed or combined, for example, the modules in the embodiment may be combined into one module, or may be further decomposed into a plurality of sub-modules, so as to complete all or part of the functions described above. The names of the modules and steps related to the embodiments of the present invention are merely for distinguishing the respective modules or steps, and are not to be construed as unduly limiting the present invention.
An electronic device of a third embodiment of the present invention includes:
at least one processor; and
a memory communicatively coupled to at least one of the processors; wherein,
the memory stores instructions executable by the processor for execution by the processor to implement a method of calculating taxi speed of an aircraft into and out of a region of aircraft as described above.
A fourth embodiment of the invention is a computer readable storage medium storing computer instructions for execution by the computer to perform a method of calculating an aircraft taxiing speed in an approach and stop area as described above.
It will be clearly understood by those skilled in the art that, for convenience and brevity of description, the specific working process of the storage device and the processing device described above and the related description may refer to the corresponding process in the foregoing method embodiment, which is not repeated herein.
Those of skill in the art will appreciate that the various illustrative modules, method steps described in connection with the embodiments disclosed herein may be implemented as electronic hardware, computer software, or combinations of both, and that the program(s) corresponding to the software modules, method steps, may be embodied in Random Access Memory (RAM), memory, read Only Memory (ROM), electrically programmable ROM, electrically erasable programmable ROM, registers, hard disk, removable disk, CD-ROM, or any other form of storage medium known in the art. To clearly illustrate this interchangeability of electronic hardware and software, various illustrative components and steps have been described above generally in terms of their functionality. Whether such functionality is implemented as electronic hardware or software depends upon the particular application and design constraints imposed on the solution. Those skilled in the art may implement the described functionality using different approaches for each particular application, but such implementation is not intended to be limiting.
Computer program code for carrying out operations of the present application may be written in one or more programming languages, including an object oriented programming language such as Java, smalltalk, C ++ and conventional procedural programming languages, such as the "C" programming language or similar programming languages. The program code may execute entirely on the user's computer, partly on the user's computer, as a stand-alone software package, partly on the user's computer and partly on a remote computer or entirely on the remote computer or server. In the case of a remote computer, the remote computer may be connected to the user's computer through any kind of network, including a Local Area Network (LAN) or a Wide Area Network (WAN), or may be connected to an external computer (for example, through the Internet using an Internet service provider).
The flowcharts and block diagrams in the figures illustrate the architecture, functionality, and operation of possible implementations of systems, methods and computer program products according to various embodiments of the present application. In this regard, each block in the flowchart or block diagrams may represent a module, segment, or portion of code, which comprises one or more executable instructions for implementing the specified logical function(s). It should also be noted that, in some alternative implementations, the functions noted in the block may occur out of the order noted in the figures. For example, two blocks shown in succession may, in fact, be executed substantially concurrently, or the blocks may sometimes be executed in the reverse order, depending upon the functionality involved. It will also be noted that each block of the block diagrams and/or flowchart illustration, and combinations of blocks in the block diagrams and/or flowchart illustration, can be implemented by special purpose hardware-based systems which perform the specified functions or acts, or combinations of special purpose hardware and computer instructions.
The terms "first," "second," and the like, are used for distinguishing between similar objects and not for describing a particular sequential or chronological order.
The terms "comprises," "comprising," or any other variation thereof, are intended to cover a non-exclusive inclusion, such that a process, method, article, or apparatus/apparatus that comprises a list of elements does not include only those elements but may include other elements not expressly listed or inherent to such process, method, article, or apparatus/apparatus.
Thus far, the technical solution of the present invention has been described in connection with the preferred embodiments shown in the drawings, but it is easily understood by those skilled in the art that the scope of protection of the present invention is not limited to these specific embodiments. Equivalent modifications and substitutions for related technical features may be made by those skilled in the art without departing from the principles of the present invention, and such modifications and substitutions will be within the scope of the present invention.

Claims (7)

1. A method of calculating a taxiing speed of an aircraft in a stopping zone, the method comprising:
step S100, acquiring an original operation parameter data set of a target aircraft in a landing taxiing stage;
step S200, eliminating abnormal values in the original operation parameter data set and performing standardization processing to obtain a standardized operation parameter data set;
step S300, calculating the moment of the aircraft to the stand according to the operation parameters in the standardized operation parameter data set；
The step S300 specifically includes:
step S310, according to the stop brake parameter, the space gate starts to traverse forward from the point when the AIR changes to GROUND, and the time when the stop brake parameter changes from the mark 0 to the mark 1 is found to be the set stop timeThe method comprises the steps of carrying out a first treatment on the surface of the When the stay braking parameter is 0, the brake is released, and when the stay braking parameter is 1, the stay braking is set;
step S320, setting the stay time according to the longitude and latitude of the airport standCalculating the arc length between the coordinate point of the aircraft and the coordinate point of the stand>：
；
；
；
；
Wherein,is the longitude of the assigned airport stand, +.>Is the latitude of the designated airport stand; />Is to set the stay timeLongitude, latitude of current aircraft>Is to set the stay time +.>Current aircraft latitude,/->Indicating deflection angle +.>Is the difference of dimensions, +.>Is the difference in longitude;
judging the arc lengthWith a preset accurate berth threshold value, if the arc length is + ->Less than or equal to a preset valueDetermining the accurate berth threshold value, and judging the current momenttThe aircraft is correctly stopped at the stand, and the current time is output +.>As aircraft to stand time->The method comprises the steps of carrying out a first treatment on the surface of the If the arc length is +>If the data is larger than a preset accurate berth threshold value, judging that the data acquired by the current aircraft cannot accurately judge that the aircraft stops on a stand;
step S330, determining the time of arrival at the standIs>；
Step S340, if the stand time is setIs>0, go to step S400; if the stand is at time->Is the taxiing speed V of the aircraft t If the current value is not 0, the point moment of the stand is recorded>Is>Reverse traversal would be within 30 seconds of the corresponding taxi speed +.>The first aircraft speed magnitude is subtracted>Ensure stand time +.>The taxiing speed is 0, and an airplane taxiing speed sequence is obtained;
step S400, according to the moment when the aircraft arrives at the standThe standard operation parameter data are used for calculating the speed of the equal difference distance from the stand of the aircraft in the process of entering the stand area of the aircraft;
the step S400 specifically includes:
step S410, obtaining standardized operation parameter data, including:
Time[n]: distance standnThe moment of the rice is initiallyNone；
LIST_GS[n]: distance standnThe taxiing speed of the rice aircraft is initiallyNone；
GS[i]: second at the present timeiThe taxiing speed of the aircraft is initiallyNone；
: cumulative distance, representing the distance moved by the aircraft per second, initially 0;
: the target observation distance from the preset airplane to the stand is set to be 1 meter;
step S420, from the aircraft to the standStart reverse cycle traverse of coasting speedGSThe iteration mode is as follows:
；
wherein the method comprises the steps ofFor the current time secondiAnd the last time secondi-1, i is the initial value from +.>Current time of start second, 1 sea/hour = 1852/3600 m/s;
step S430, the steps are performedAccumulated to accumulated distance +.>In (I)>；
Step S440, if the accumulated distanceLess than->Advancing the time second i by 1 second, and repeating the steps S420 and S430; when accumulating distance->First time greater than or equal to the target viewing distance +.>Then it is stated that +.>Distance from rice, i.e.)>Jump to S450;
step S450, ifGS[i]=GS[i-1]Distance to stand target observation distanceCoasting speed atIs thatGS[i-1]；
If it isGS[i]≠GS[i-1]According to the accumulated distance、/>And i-1 and the current second i find the viewing distance from the stand target>Speed at rice:
a=GS[i-1]-GS[i] ；
b=GS[i]*2；
；
；
；
by->Rounding to take a decimal place as the observation distance from the shutdown target>A sliding speed at the rice;
step S460, object observation distanceIncreasing the viewing distance of the new target by 1 meterIf the new target is observed at a distance +>Marking the new target observation distance as the target observation distance, repeating the steps S420-S450 until +.>；
Step S470, output，LIST_The value in GS sequentially corresponds to the ground speed of 1 to 10 meters from the stand as the speed of each moment in the process of entering the stand area of the airplane.
2. The method of calculating the taxiing speed of an aircraft in a stopping area according to claim 1, wherein the step S200 specifically includes:
step S210, eliminating or replacing abnormal operation parameter data based on the original operation parameter data set to obtain an operation parameter data set with abnormal values eliminated;
the abnormal operation parameter data comprise frame dropping parameters and operation parameter data with parameter records exceeding a preset abnormal threshold;
step S220, data cleaning and data standardization are carried out on the operation parameter data set with abnormal values removed, and a standardized operation parameter data set is obtained.
3. A method of calculating the taxiing speed of an aircraft in a stopping area according to claim 1, characterized in that said standardized operating parameter data set comprises in particular:
the aircraft stand longitude and latitude, the aircraft longitude and latitude, the sliding speed, the stay brake and the air-ground electric door flight recording parameters;
the longitude and latitude of the stand and the longitude and latitude of the airplane use a unified coordinate system;
the stay brake parameter includes a brake-on flag of 1 and a brake-off flag of 0.
4. The method of calculating the taxiing speed of an aircraft in a stopping area according to claim 1, further comprising step S500:
and comparing the speed of the aircraft at each moment in the aircraft entering and exiting station area with the interpolation of the standard value to obtain the safe operation evaluation result of the aircraft entering and exiting station stage.
5. A system for calculating taxi speeds of an aircraft in and out of a region, said system comprising:
the original parameter acquisition module is configured to acquire an original operation parameter data set of the target aircraft in a landing taxiing stage;
the preprocessing module is configured to remove abnormal values in the original operation parameter data set and perform standardized processing to obtain a standardized operation parameter data set;
the stand time calculation module is configured to calculate the time when the aircraft arrives at the stand according to the operation parameters in the standardized operation parameter data set；
The method specifically comprises the following steps:
according to the stop brake parameter, the space electric gate starts to go forward from the point that the AIR changes to the GROUND state, and the moment when the stop brake parameter changes from the mark 0 to the mark 1 is found and recorded as the set stop momentThe method comprises the steps of carrying out a first treatment on the surface of the When the stay braking parameter is 0, the brake is released, and when the stay braking parameter is 1, the stay braking is set;
according to the longitude and latitude of the airport stand and the setting of the stay timeCalculating the arc length between the coordinate point of the aircraft and the coordinate point of the stand>：
；
；
；
；
Wherein,is the longitude of the assigned airport stand, +.>Is the latitude of the designated airport stand; />Is to set the stay timeLongitude, latitude of current aircraft>Is to set the stay time +.>Current aircraft latitude,/->Indicating deflection angle +.>Is the difference of dimensions, +.>Is the difference in longitude;
judging the arc lengthWith a preset accurate berth threshold value, if the arc length is + ->If the preset accurate berth threshold value is smaller than or equal to the preset accurate berth threshold value, judging the current momenttThe aircraft is correctly stopped at the stand, and the current time is output +.>As aircraft to stand time->The method comprises the steps of carrying out a first treatment on the surface of the If the arc length is +>If the data is larger than a preset accurate berth threshold value, judging that the data acquired by the current aircraft cannot accurately judge that the aircraft stops on a stand;
determining the arrival time of the machine at the standIs>；
If the stand time is setIs>0, carrying out a speed change condition acquisition module; if the stand is at time->Is the taxiing speed V of the aircraft t If the current value is not 0, the point moment of the stand is recorded>Is>Reverse traversal would be within 30 seconds of the corresponding taxi speed +.>All minus the first aircraft speed magnitudeEnsure stand time +.>The taxiing speed is 0, and an airplane taxiing speed sequence is obtained;
a speed change condition acquisition module configured to, according to the aircraft arrival timeThe operation parameter data are standardized, and the speed of each moment in the process of entering the stand area of the aircraft is calculated;
an initializing unit for obtaining standardized operation parameter data, including:
Time[n]: distance standnThe moment of the rice is initiallyNone；
LIST_GS[n]: distance standnThe taxiing speed of the rice aircraft is initiallyNone；
GS[i]: second at the present timeiThe taxiing speed of the aircraft is initiallyNone；
: cumulative distance, representing the distance moved by the aircraft per second, initially 0;
: the target observation distance from the preset airplane to the stand is set to be 1 meter;
taxiing speed iteration unit for time from airplane to standStart reverse cycle traverse of coasting speedGSThe iteration mode is as follows:
；
wherein the method comprises the steps ofFor the current time secondiAnd the last time secondi-1, distance betweenFrom increment, i is the initial value from +.>Current time of start second, 1 sea/hour = 1852/3600 m/s;
an accumulation unit for accumulating the aboveAccumulated to accumulated distance +.>In (I)>；
An accumulated distance iteration unit for calculating the accumulated distanceLess than->Advancing the time second i for 1 second, and repeating the sliding speed iteration unit and the accumulation unit; when accumulating distance->First greater than or equal to the target viewing distanceThen it is stated that +.>Distance from rice, i.e.)>Jumping to a sliding speed calculating unit;
a sliding speed calculating unit, ifGS[i]=GS[i-1]Distance to stand target observation distanceThe sliding speed +.>Is thatGS[i-1]；
If it isGS[i]≠GS[i-1]According to the accumulated distance、/>And i-1 and the current second i find the viewing distance from the stand target>Speed at rice:
a=GS[i-1]-GS[i] ；
b=GS[i]*2；
；
；
；
by->Rounding to take a decimal place as the observation distance from the shutdown target>A sliding speed at the rice;
target observation distance updating unit for updating target observation distanceIncreasing the viewing distance of the new target by 1 meterIf the new target is observed at a distance +>Marking the new target observation distance as the target observation distance, and repeating the sliding speed iteration unit to the sliding speed calculation unit until；
Data output unit for outputting，LIST_The value in GS sequentially corresponds to the ground speed of 1 to 10 meters from the stand as the speed of each moment in the process of entering the stand area of the airplane.
6. An electronic device, comprising:
at least one processor; and
a memory communicatively coupled to at least one of the processors; wherein,
the memory stores instructions executable by the processor for performing a method of calculating aircraft taxi speed into a taxi stop zone as claimed in any one of claims 1 to 4.
7. A computer readable storage medium storing computer instructions for execution by the computer to implement a method of calculating aircraft taxiing speed in a stopping area as claimed in any one of claims 1 to 4.
CN202311744250.1A 2023-12-19 2023-12-19 Method, system and device for calculating sliding speed of aircraft in and out of service area Active CN117421542B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202311744250.1A CN117421542B (en) 2023-12-19 2023-12-19 Method, system and device for calculating sliding speed of aircraft in and out of service area

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202311744250.1A CN117421542B (en) 2023-12-19 2023-12-19 Method, system and device for calculating sliding speed of aircraft in and out of service area

Publications (2)

Publication Number Publication Date
CN117421542A true CN117421542A (en) 2024-01-19
CN117421542B CN117421542B (en) 2024-03-29

Family

ID=89530673

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202311744250.1A Active CN117421542B (en) 2023-12-19 2023-12-19 Method, system and device for calculating sliding speed of aircraft in and out of service area

Country Status (1)

Country Link
CN (1) CN117421542B (en)

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107944701A (en) * 2017-11-23 2018-04-20 北京航空航天大学 A kind of detection method and device for the risk that guns off the runway during aircraft landing
CN110288857A (en) * 2019-06-26 2019-09-27 中国民航大学 It is a kind of based on when Spatial Dimension decompose airport aircraft gate fast dispatch method
CN114895705A (en) * 2022-05-07 2022-08-12 中国民航科学技术研究院 Method and device for calculating turning radius of airplane in ground taxiing stage
EP4109433A1 (en) * 2021-06-22 2022-12-28 ADB Safegate Sweden AB Method for monitoring backward movement of an aircraft at an airport stand

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107944701A (en) * 2017-11-23 2018-04-20 北京航空航天大学 A kind of detection method and device for the risk that guns off the runway during aircraft landing
CN110288857A (en) * 2019-06-26 2019-09-27 中国民航大学 It is a kind of based on when Spatial Dimension decompose airport aircraft gate fast dispatch method
EP4109433A1 (en) * 2021-06-22 2022-12-28 ADB Safegate Sweden AB Method for monitoring backward movement of an aircraft at an airport stand
CN114895705A (en) * 2022-05-07 2022-08-12 中国民航科学技术研究院 Method and device for calculating turning radius of airplane in ground taxiing stage

Also Published As

Publication number Publication date
CN117421542B (en) 2024-03-29

Similar Documents

Publication Publication Date Title
Medeossi et al. A method for using stochastic blocking times to improve timetable planning
RU2573735C2 (en) Method and system for analysis of flight data recorded during aircraft flight
US20150032418A1 (en) Aggregating and filtering data from mobile devices
CN109493449B (en) Freight train cargo state estimation method based on freight train GPS track data and high-speed transaction data
CN112732687A (en) Aviation flight data visualization processing system and analysis method based on data cleaning
CN109478350B (en) System and method for determining maintenance-focused aircraft data record frame configuration
CN103473436A (en) Method and device for assisting the mission tracking of an aircraft
EP4020425A2 (en) Method and apparatus for determining green wave speed, electronic device and storage medium
Cunillera et al. A literature review on train motion model calibration
CN104977922A (en) Equipment monitoring device and method thereof
Wang et al. Airspace risk management using surveillance track data: Stabilized approaches
Cunillera et al. Real-time train motion parameter estimation using an unscented Kalman filter
CN117421542B (en) Method, system and device for calculating sliding speed of aircraft in and out of service area
CN102495305A (en) Grey theory-based fault trend analysis method of track traffic operation electromechanical equipment
Figuet et al. Data-driven mid-air collision risk modelling using extreme-value theory
Bilimoria Methodology for the performance evaluation of a conflict probe
Hrastovec et al. Machine learning model for aircraft performances
CN117953731A (en) Incoming flight flow flight plan prediction method for terminal area traffic simulation
Hu et al. Runway occupancy time constraint and runway throughput estimation under reduced arrival wake separation rules
CN111967676A (en) Method and system for predicting risk of aircraft tail rubbing during takeoff based on stepwise regression
CN116415818A (en) Method and system for confirming risk points in aircraft approach stage based on clustering algorithm
CN106096152A (en) A kind of aircraft safety landing processing method based on function resonance analyzing
Kistler et al. Relationship between airport efficiency and surface traffic
Yan et al. Precognitive maintenance and probabilistic assessment of tool wear using particle filters
CN109190968B (en) Empty pipe operation efficiency grade evaluation method and device

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant