Международный информационно-аналитический журнал «Crede Experto: транспорт, общество, образование, язык». № 1 (36). Март 2023 (http://ce.if-mstuca.ru)
УДК 629.7.058.6 : 629.7.016.2 DOI 10.51955/2312-1327_2023_1_118
ПОДТВЕРЖДЕНИЕ ДОСТОВЕРНОСТИ СООБЩЕНИЙ АЗН-В ПУТЕМ ОЦЕНКИ ВЫСОТЫ ПОЛЕТА ВОЗДУШНОГО СУДНА
Александр Петрович Плясовских, orcid.org/0000-0003-2250-8852, доктор технических наук, доцент АО «ВНИИРА», пр. Обуховской обороны, д. 120 Санкт-Петербург, 192012, Россия info@vniira. ru
Евгений Андреевич Рубцов, orcid.org/0000-0003-2126-0015, кандидат технических наук Российский университет транспорта (МИИТ),
ул. Образцова, д. 9, стр. 9 Москва, 127994, Россия rubtsov.rut.miit@gmail. com
Андрей Сергеевич Калинцев, orcid.org/0000-0002-2578-2892, ФГУП «Госкорпорация по ОрВД», Архангельский центр ОВД филиала «Аэронавигация Северо-Запада»,
ул. Северная, д. 20 Мезень, 164751, Россия kas4job@gmail.com
Владислав Юрьевич Давиденко, orcid.org/0009-0001-9707-166X, ФГУП «Госкорпорация по ОрВД» филиал «Аэронавигация Северо-Запада», 10 линия В.О., д. 53А Санкт-Петербург, 199048, Россия david071521@mail. ru
Аннотация. Предложена методика подтверждения данных АЗН-В, основанная на анализе информации о высоте воздушного судна и параметрах атмосферы. Для подтверждения достоверности сообщений АЗН-В предлагается сравнивать геометрическую высоту, полученную по данным спутниковой системы навигации, с высотой, рассчитанной для реальной атмосферы (по информации от метеослужбы о давлении и температуре). Для реализации предлагаемой методики необходимы следующие входные данные: давление на уровне метеостанции, температура на уровне метеостанции и на высотах выполнения полетов и передаваемые в сообщении АЗН данные о геометрической и барометрической высотах. Установлено, что наибольшую погрешность при выполнении расчетов вносит температура на высоте полета ВС. С целью снижения влияния ошибок определения температуры предложено использование данных от метеозондов. С помощью разработанной методики выполнен расчет с использованием реальных данных о высотах, передаваемых в сообщении АЗН-В и требуемых данных от метеослужбы. Были использованы сообщения, полученные от наземной станции АЗН-В НС-1А производства АО «ВНИИРА», установленной на аэродроме Мезень.
© А.П.Плясовских, Е.А.Рубцов, А.С.Калинцев, В.Ю. Давиденко, 2023
Разработанную методику предлагается применять для подтверждения данных АЗН-В на региональных аэродромах с малой интенсивностью полетов.
Ключевые слова: АЗН-В, 1090 ES, метеозонд, барометрическая высота, геометрическая высота, давление, температура, формула Лапласа, ISA.
CONFIRMATION OF ADS-B MESSAGES BY AIRCRAFT FLIGHT
ALTITUDE
Alexander P. Plyasovskih, orcid.org/0000-0003-2250-8852, Doctor of Technical Sciences JSC VNIIRA, 120, Obukhovskoy Oborony av.
Saint-Petersburg, 192012, Russia info@vniira. ru
Evgeny A. Rubtsov, orcid.org/0000-0003-2126-0015, Candidate of Technical Sciences Russian university of transport (MIIT), 9, buid. 9, Obrazcova str. Moscow, 127994, Russia rubtsov.rut.miit@gmail. com
Andrey S. Kalintsev, orcid.org/0000-0002-2578-2892, State ATM Corporation of the Russian Federation, Arkhangelsk ATM center,
North-West Air Navigation branch 20, Severnaya str.
Mezen, 164751, Russia kas4job@gmail.com
Vladislav Y. Davidenko, orcid.org/0009-0001-9707-166X State ATM Corporation of the Russian Federation, North-West Air Navigation branch,
53A, 10 Line of V. I.
Saint-Petersburg, 199048, Russia david071521@mail. ru
Abstract. Proposed the method for confirming ADS-B data based on the analysis of information about the aircraft altitude and atmospheric parameters. To confirm the reliability of ADS-B messages, proposed to compare the geometric height obtained from satellite navigation system data with the height calculated for the real atmosphere (based on information from the meteorological service on pressure and temperature). To implement the proposed methodology, required the following input data: pressure at the level of the weather station, temperature at the level of the weather station and at flight altitudes, and data on geometric and barometric altitudes transmitted in the ADS message. It has been established that the greatest error in the calculations is introduced by the temperature at the aircraft flight altitude. In order to reduce the influence of errors in determining the temperature, proposed to use data from weather balloons. With developed methodology, a calculation was performed using real altitude data transmitted in the ADS-B message and the required data from the meteorological service. Messages received from the ADS-B ground station NS-1A manufactured by JSC VNIIRA, which is installed at the Mezen airfield, were used. The developed methodology is proposed to be used to confirm ADS-B data at regional aerodromes with low flight intensity.
Key words: ADS-B, weather balloon, barometric altitude, geometric altitude, pressure, temperature, Laplace formula, ISA.
Введение
Автоматическое зависимое наблюдение вещательного типа (АЗН-В) отнесено Международной организацией гражданской авиации (ИКАО) к перспективным технологиям наблюдения [Aeronautical Surveillance..., 2020; Global air., 2016]. К достоинствам АЗН-В относят небольшую стоимость внедрения наземных станций по сравнению с затратами на развертывание вторичных радиолокаторов (ВРЛ) и многопозиционных систем наблюдения (МПСН). Эксплуатационная стоимость также значительно ниже. При наличии на борту высокоточных навигационных датчиков (как правило, приемников глобальных навигационных спутниковых систем - ГНСС) АЗН-В обеспечивает точность наблюдения не хуже, чем у радиолокационных средств.
Однако, в Руководстве по авиационному наблюдению ИКАО [Aeronautical Surveillance., 2020] отмечается необходимость подтверждать достоверность получаемых от АЗН-В данных с помощью ВРЛ или МПСН. Это связано с возможностью возрастания погрешности навигации (при использовании автономных и азимутально-дальномерных систем навигации). Таким образом, необходимость подтверждения данных нивелирует экономическую привлекательность технологии АЗН-В.
Задача подтверждения данных АЗН-В без применения дополнительных средств наблюдения может быть решена с применением алгоритмических методов. В работах [Калинцев и др., 2021; Плясовских и др., 2019; Плясовских и др., 2020] были предложены методы подтверждения данных АЗН-В для обеспечения наземного движения на аэродроме и полетов в аэродромной зоне. Эти методы основаны на накоплении измерений координат воздушного судна (ВС) на станции АЗН-В с последующим статистическим анализом полученных данных.
В гражданской авиации данные о высоте ВС имеют большое значение: на их основе выполняется эшелонирование ВС и управление воздушным движением (УВД). Сообщение АЗН-В содержит информацию о горизонтальном и вертикальном местоположении ВС.
Вертикальная точность определения местоположения ГНСС значительно уступает горизонтальной. Согласно [Taib et al., 2016] погрешности при определении геометрической высоты могут превышать погрешности определения горизонтального местоположения в 3 раза. В нормальных условиях погрешность определения геометрической высоты находится в пределах 30-50 м [Taib et al., 2016]. Очевидно, что неоптимальное расположение спутников, отказ или нарушение работоспособности оборудования определения местоположения ВС приведут к большим погрешностям определения вертикального местоположения. Учитывая значение данных о высоте ВС и вероятность значительных погрешностей при определении геометрической высоты АЗН-В, можно заключить, что разработка методики подтверждения данных АЗН-В путем оценки высоты полета ВС является важной и актуальной задачей. Решение
этой задачи сделает возможным применение АЗН-В в качестве единственного источника наблюдения, что, принимая во внимание экономическую привлекательность данной технологии, позволит снизить стоимость развертывания и эксплуатации инфраструктуры наблюдения, что согласуется с документами ИКАО по планированию развития гражданской авиации [Global air..., 2016].
В сообщении АЗН-В передаются данные о геометрической и барометрической высоте. Барометрическая высота передается в виде значения занимаемого ВС эшелона полета относительно среднего уровня моря (mean sea level, MSL). Данные о барометрической высоте передаются по показаниям бортового высотомера, являющегося независимым от системы ГНСС источником данных о высоте полета ВС. Барометрический высотомер откалиброван в соответствии с Международной стандартной атмосферой (International Standard Atmosphere, ISA), предполагающей стандартное давление на уровне моря 1013,25 гПа при 15° [Taib et al., 2016; Manual of the ICAO..., 1993]. Сообщения расширенного сквиттера о положении ВС в воздухе помимо кодированных значений широты и долготы содержат данные о геометрической высоте. Данные ГНСС о геометрической высоте представляют собой высоту над эллипсоидом в системе координат WGS-84 [Minimum operational., 2009].
Барометрический высотомер и приемник сигналов ГНСС измеряют принципиально разные величины, однако независимость этих источников позволяет использовать измерения барометрического высотомера для проверки геометрической высоты ВС и реализовать методику подтверждения данных АЗН-В путем оценки высоты полета ВС.
Материалы и методы
Для реализации методики подтверждения данных АЗН-В путем оценки высоты полета ВС необходимы следующие исходные данные:
1) Прогнозные значения температуры на высотах, данные метеозонда, данные давления от метеостанции.
2) Высота превышения ВПП;
3) Данные, передаваемые в сообщении АЗН-В:
- Барометрическая высота;
- Геометрическая высота, координаты ВС.
Рассмотрим передаваемые значения о высоте их различия.
Геометрическая высота, передаваемая в сообщении АЗН-В [Global air.,
2016], вычисляется бортовым приемником ГНСС в системе координат WGS84. В сообщении стандарта 1090ES значение геометрической высоты передается в поле «GNSS Height (НАЕ)», то есть как превышение над эллипсоидом в системе координат WGS84 [Minimum operational..., 2009]. Стоит отметить, что возможны два варианта передачи геометрической высоты [Minimum operational., 2009]:
1) как непосредственное значение высоты «GNSS Height (НАЕ)»;
2) как разницы между барометрической и геометрической высотами «Difference between "Baro Altitude" and "GNSS Height (НАЕ)"».
В обоих случаях значения геометрической высоты определяются как превышения над эллипсоидом WGS84 [Global air., 2016].
Эллипсоид WGS84 в вертикальной плоскости в качестве приближения к MSL, согласно [Aeronautical Surveillance., 2020], может иметь ошибки от -100 до +70 м по отношению к геоиду (рис. 1), конкретные значения ошибок зависят от местоположения на земном шаре. Чтобы сравнить данные рассчитанной высоты для реальной атмосферы и данных о высоте полета ГНСС, их необходимо привести к одному уровню. Предлагается приводить все высоты к среднему уровню моря (MSL).
Рисунок 1 - Схематичное изображение уровней: геоида, эллипсоида WGS84 и высоты, определяемой приемником ГНСС
На рис. 1 обозначено: H - высота над геоидом, которая определяется гравитацией Земли и аппроксимируется средним уровнем моря (MSL) [Use International Standard., s.a.]. По координатам от ВС находим волну геоида N (разность между поверхностью эллипсоида WGS84 и геоида). Для показанного на рис. 1 случая, для нахождения H нужно сложить высоту ВС над эллипсоидом h и волну геоида N в точке с координатами ВС:
H = h + N (1)
В случае, если поверхность эллипсоида проходит над геоидом, необходимо для нахождения H вычесть из высоты над эллипсоидом h волну геоида N в точке с координатами ВС, которые получаем из сообщения АЗН-В:
H = h - N (2)
Барометрический высотомер в кабине ВС откалиброван по геопотенциальной высоте. Таким образом, конкретной высоте (эшелону полета) соответствует определенное значение давления по стандартной атмосфере. Геопотенциальная высота выражается в геопотенциальных метрах (гп. м). Согласно документу ГОСТ 4401-81 «Атмосфера стандартная, параметры» геопотенциальная высота может быть представлена как:
H = —= — "hg (h)dh (3)
gc gc о
где H - геопотенциальная высота; — - геопотенциал;
gc - стандартное ускорение свободного падения;
g(h) - ускорение свободного падения на высоте h, определяемое как:
Я = gc
r
(4)
V г + к у
где г - условный радиус Земли (равен 6356767 м); к - геометрическая высота.
Таким образом, геопотенциальная высота (Н) и геометрическая высота (к) преобразуются в соответствии с документом ГОСТ 4401-81 «Атмосфера стандартная, параметры»:
гк
(5)
H =
к =
r + к
rH r - H
(6)
Для реализации методики подтверждения данных АЗН-В путем оценки высоты полета ВС необходимы следующие исходные данные:
- превышение метеостанции над уровнем моря - hi;
- барометрическая высота полета ВС - Явебаро;
- геометрическая высота полета ВС - ЯвегМ;
- температура на уровне метеостанции - to;
- температура на высоте полета ВС - fa;
- атмосферное давление на высоте полета ВС - Рвс;
- атмосферное давление на уровне метеостанции - Po.
Данные АЗН-В включают барометрическую высоту НВСбаро и геометрическую высоту Нвсгм полета ВС [Minimum operational..., 2009].
Предоставляемые метеослужбой данные включают: прогнозные и фактические значения давления Po и температуры to и tH, которые получают от метеосервера или от другого источника метеоданных. Барометрическую и геометрическую высоту полета получают из сообщения АЗН-В от воздушного судна. Данные о превышении ВПП (давление, приведенное к уровню порога ВПП) получают из сборника аэронавигационной информации (AIP) в разделе географические и административные данные по аэродрому.
Согласно разработанной методике, геометрическая высота воздушного судна находится как сумма высоты превышения метеостанции над уровнем моря к\ и высоты от уровня метеостанции до ВС кг (рис. 2).
Высота ВС
Уровень метеостанции
h.
Уровень моря (MSL)
К
Рисунок 2 - К задаче определения геометрической высоты АЗН-В
Подтверждение достоверности данных АЗН-В осуществляют следующим образом:
1. Определяют высоту превышения ВПП над уровнем моря - hi с использованием сборника аэронавигационной информации (AIP).
2. Определяют фактические значения давления Po и температуры to на уровне метеостанции. Данные получают из сводок METAR, которые обновляются каждые полчаса и содержат фактические значения давления (QFE) и температуры на уровне метеостанции.
3. Определяют барометрическую высоту полета ВС Нвсбаро, значение которой содержится в сообщении АЗН-В.
4. Находят атмосферное давление Рвс на высоте полета ВС Нвсбаро, согласно данным документа ГОСТ 4401-81 «Атмосфера стандартная, параметры». Для высот до 11 000 м давление определяется как:
Р = Ро ^ьк;, (7)
где Ро - давление на уровне моря (101325 Па);
Ь - градиент температуры в соответствии с параметрами стандартной атмосферы (0,0065 К/м);
Я - газовая постоянная (287,0531 );
кг - К
е = T
To
T = To - L • h ;
To - температура поверхностного слоя стандартной атмосферы (288,15 К);
h - барометрическая (геопотенциальная) высота ВС.
5. Определяют температуру на высоте полета воздушного судна tH с использованием прогнозных значений, получаемых от метеослужбы. Прогнозные значения корректируются с использованием данных метеозонда, которые содержат фактическое значение температуры по высотам. Метеозонд запускают каждые 12 часов. Данные метеозонда находятся в открытом доступе (табл. 1) [Department of Atmospheric..., 2023]. Для их получения необходимо знать код метеостанции, от которой запускался радиозонд, в районе, где находится ВС (например, для Архангельска и Архангельской области необходимо применять данные от метеостанции с кодом 22543). Данные о значениях температуры записываются с разной дискретностью, поэтому для получения значения температуры на искомой высоте воспользуемся выражением:
t =
i HВС -
v H+1 - H
(tj+1- tj) + tj' (8)
где £ - температура в градусах на высоте Нвс, по данным метеозонда; Нвс - геометрическая высота АЗН-В, полученная с борта ВС; Н - ближайшее значение высоты к Нвс, полученное из таблицы показаний радиозонда в меньшую сторону;
Н+1 - ближайшее значение высоты к Нвс, полученное из таблицы показаний метеозонда в большую сторону;
^ - температура, соответствующая высоте Н/, полученная по данным метеозонда;
- температура, соответствующая высоте Н+1, полученная по данным метеозонда.
Таблица 1 - Пример данных метеозонда для метеостанции №22543 (Архангельск)
Атмосферное Геопотенциальная Температура, Относительная Скорость
давление, высота, м °С влажность, % ветра,
ГПа узлы
1000 6 2.8 91 8
974 217 1.8 92 25
925 629 -0,3 94 33
850 1301 -3,9 94 39
746 2322 -8,7 96 43
725 2545 -9,7 96 42
700 2815 -12,3 93 41
669 3156 -15,3 94 39
633 3572 -18,9 94 43
574 4297 -20,7 78 50
503 5257 -27,6 68 60
500 5300 -27,9 67 60
448 6083 -32,9 65 60
417 6581 -37,3 65 60
400 6870 -39,9 64 60
397 6921 -40,5 66 60
316 8443 -51,3 67 68
6. Рассчитывают геопотенциальную высоту полета ВС относительно уровня метеостанции . С помощью барометрической формулы Лапласа:
= 18400. (1 + ях) \ё , (9)
^ВС
где - геопотенциальная высота ВС относительно уровня метеостанции;
.Роге - давление на уровне метеостанции;
Рвс - давление в точке полета ВС;
Я - коэффициент, равный 1/273;
t - средняя температура слоя атмосферы.
Средняя температура слоя t рассчитывается по формуле:
Х0 + Хи
х =
2
где ¿о - температура на уровне метеостанции; tя - температура на высоте полета ВС.
7. Рассчитывают геометрическую высоту ВС относительно метеостанции. Перевод высоты из геопотенциальной в геометрическую осуществляется в соответствии с документом ГОСТ 4401-81 «Атмосфера стандартная, параметры». Геопотенциальная высота ВС (H) и его геометрическая высота (h) преобразуются согласно (5) и (6).
8. Рассчитывают геометрическую высоту полета воздушного судна относительно уровня моря (Нвсгм):
тт гм т, . т,гм /11ч
hвс = hl + "2 (11)
где hi - превышение метеостанции над уровнем моря;
h2rM - геометрическая высота от воздушного судна относительно ВПП.
9. Выполняют сравнение полученных высот Нвсгм и Hads-b. Данные АЗН-В о геометрической высоте представляют собой значения превышения над эллипсоидом в системе координат WGS84. С использованием существующих баз данных волны геоида [Calculation of Gravity..., s.a.; EGM96..., s.a.; Evaluation of the EGM96..., s.a.] выполняется приведение к уровню MSL значений геометрической высоты, получаемой из сообщения АЗН-В. При этом модуль разности вычисленных значений не должен превышать максимально допустимой погрешности измерений:
hвс - hads-b <^макс._доп (12)
где Hads-b - приведенная к уровню MSL геометрическая высота ВС, полученная из сообщения ADS-B.
Обоснование значения максимально допустимой погрешности Амакс._доп
требует проведения отдельного исследования, поэтому в данной работе примем, что отклонение между высотами |нВСГМ -HADS-^ не должно превышать 25% от
нормы вертикального эшелонирования 300 м. Таким образом, Амакс._доп принимается равной 75 м.
Если разность высот |нВСГМ - HADS-^ меньше Амакс._доп, в таком случае
геометрическую высоту, полученную из сообщения АЗН-В, считают достоверной, а данные АЗН-В подтвержденными: их можно применять без необходимости сверки с данными других средств наблюдения. В противном случае данные о геометрической высоте сообщения АЗН-В считаются недостоверными, таким образом информации АЗН-В доверять нельзя и необходимо либо применять иные средства наблюдения, либо процедурные методы УВД.
Результаты
Выполним процедуру подтверждения данных АЗН-В с использованием разработанной методики. В данном примере используются прогнозные значения температуры по высотам (для наглядности алгоритм коррекции температуры по данным метеозонда вынесен отдельным пунктом). Необходимо отметить, что коррекция данных температуры выполняется при поступлении данных зондирования и эти данные имеют больший, в сравнении с прогнозными
значениями, приоритет. При определении температуры в точке полета ВС используется линейная интерполяция.
Рассмотрим случай полета ВС местных воздушных линий на эшелоне 060 (высота 6000 футов или 1850 м). В качестве исходных данных использовались реальные данные, полученные от наземной станции АЗН-В НС-1А производства АО «ВНИИРА», установленной на аэродроме Мезень, и данные от метеослужбы. Данные АЗН-В, полученные от бортового ответчика, включают геометрическую и барометрическую высоты. От метеослужбы получены данные QFE, температуры на уровне метеостанции и прогнозные значения температуры на различных высотах. Значения высоты превышения ВПП получены из сборника AIP. Так, для ВПП аэродрома Мезень высота превышения ВПП Нрев составляет 16 м. Значения исходных данных представлены в табл. 2 и табл. 3.
Таблица 2 - Исходные данные, полученные от наземной станции АЗН-В НС-1А производства АО «ВНИИРА», установленной на аэродроме Мезень
Описание параметра, единицы измерения Обозначение Значение параметра
Геометрическая высота, передаваемая в составе сообщения АЗН-В, м Hgeom 1630,22
Барометрическая высота, передаваемая в составе сообщения АЗН-В, фут Hbaro 6025
Таблица 3 - Исходные данные, полученные от метеослужбы аэродрома Мезень
Описание параметра, единицы измерения Обозначение Значение параметра
Давление ОТЕ, Па Pqfe 99848
Температура на уровне метеостанции, °С t0 -4
Прогнозное значение для высоты 600 м, °С tmz600 -11
Прогнозное значение для высоты 1500 м, °С tmz1500 -10
Прогнозное значение для высоты 3000 м, °С tmz3000 -17
Выполним шаги расчета согласно разработанной методике:
1. Высота превышения ВПП аэродрома Мезень: Нпрев=16 м.
2. Фактические значения давления и температуры на уровне метеостанции аэродрома Мезень: ^е = 99848 Па, ^ = -4°С.
3. Барометрическая высота, передаваемая в сообщении АЗН-В: Ньаго= 6025 фут. Учитывая, что 1 метр равняется 3,281 фут, выразим высоту в метрах: Ньаго_ш = 6025 / 3,281 = 1836,42 м.
4. Атмосферное давление на высоте полета ВС находим согласно выражению (7). Учитывая, что в стандартной атмосфере давление на уровне моря Ро составляет 101325 Па, температура поверхностного слоя стандартной атмосферы То равна 288,15 К, градиент температуры Ь составляет 0,0065 К/м,
газовая постоянная Я равна 287,0531 ^^, стандартное ускорение свободного
кг ■ К
падения gc составляет 9,8067 м/с2, получим:
р = Р-е
=101325-
ь-к
Р0-
т0 - ь-и
т
ь-к
288,15 - 0,0065-1836,42
9,8067
0,0065-287,0531
= 81123 Па.
288,15
5. Определим температуру на высоте полета ВС с использованием температуры на уровне метеостанции: ^ = -4°С и прогнозных значений температуры: -11°С, -10°С и -17°С для высот 600, 1500 и 3000 м соответственно. Прогнозное значение температуры ¿т^шег для высоты полета Нвс = 1630,22 м может быть получено с использованием линейной интерполяции по двум значениям температур на высотах Н1500 = 1500 м и Н3000 = 3000 м (рис. 3):
((тг1500 - ^3000)(НВС -Н1500) _
£
тгЫег ¿пк1500
(Н3000 - Н1500 )
(-10 + 17)-(1630,22-1500)
= -10-----^--¿ = -10,6 С.
(3000 -1500)
Температура.
Рисунок 3 - К задаче определения прогнозного значения температуры ¿т^шег для
высоты полета ВС 1630,22 м
Среднее значение температуры слоя составит:
-7,3 °С
Следует отметить, что при наличии данных метеозонда, полученных для данного аэродрома или района (см. табл. 1), необходимо произвести коррекцию температур для заданных высот.
6. Геопотенциальная высота полета ВС относительно уровня метеостанции, найденная с применением барометрической формулы Лапласа (9), составит:
кШ = 18400-(1+ К) ^ = 18400■( 1 --^к99848 = 1615,2 м.
2 РВС V 273 81123
7. Геометрическая высота, найденная с помощью выражения (6), равна:
гН 6356767-1615,2
к =-=-= 1615,6 м.
г -Н 6356767 -1615,2
8. Высота превышения ВПП аэродрома Мезень составляет 16 м, следовательно, геометрическая высота полета ВС относительно уровня моря согласно выражению (11) равна:
НВМ = к + кГ™ = 16 +1615,6 = 1631,6 м
9. Для сравнения рассчитанной геометрической высоты полета ВС и высоты, передаваемой в сообщении АЗН-В (равной 1630,22 м), необходимо привести значение геометрической высоты к уровню МБЬ согласно выражениям (1) или (2). Значение волны геоида N для рассматриваемого случая полета ВС в районе аэродрома Мезень согласно модели ЕОМ96 [12-14] составляет 34,5 м, при этом поверхность эллипсоида проходит над геоидом, следовательно, применяя выражение (2), получим приведенное к уровню МБЬ значение высоты полета:
Н = к - N = 1630,22 - 34,5 = 1595,72 м
Отсюда, модуль разности высот Амакс._доп согласно (12) составит:
^макс._доп
и гм тг нвс - hads-b
= |1631,6 -1595,72 = 35,88
м.
Этот результат меньше установленного порога 75 м. Таким образом, в рассмотренном примере данные АЗН-В подтверждаются.
Дискуссия
Проблеме совершенствования АЗН-В посвящено достаточно большое количество работ, включающих анализ характеристик линий передачи данных, методов обеспечения защиты от помех и преднамеренного вмешательства, повышения точности и достоверности передаваемой информации. Так, в работе [Сравнение пропускной., 2020] произведено сравнение пропускной способности различных линий передачи данных, применяемых для АЗН-В. В большинстве работ анализируется линия 1090ЕБ, являющаяся стандартом для многих стран мира, в том числе и для России. В работе [Косьянчук и др., 2019] выполнен анализ проблемы низкой защищенности линии 1090ЕБ, а также приведены возможные варианты решения этой проблемы. В работе [Алипов и др., 2019] выполнен сравнительный анализ информации, полученной от станции АЗН-В и радиолокатора, проанализированы невязки рассматриваемых данных и было установлено, что для московского региона порядка 30% ВС имеют дополнительную среднюю ошибку измерения координат 75 м.
Задача подтверждения данных АЗН-В без применения дополнительных средств наблюдения рассмотрена в работах [Калинцев и др., 2021; Плясовских и др., 2019; Плясовских и др., 2020], предложены методы алгоритмического подтверждения данных АЗН-В для обеспечения наземного движения на рабочей площади аэродрома и в воздушном пространстве. Методы, представленные в работах [Плясовских и др., 2019; Плясовских и др., 2020], основаны на накоплении измерений координат ВС с последующим статистическим анализом полученных данных. Необходимо отметить, что в перечисленных работах речь идет об ошибках и способах подтверждения данных АЗН-В в горизонтальной плоскости.
Применительно к высотам в работе [Лебедев и др., 2015] предложен новый метод мониторинга вертикального эшелонирования ВС на основе использования АЗН-В, в предложенном методе выполняется попарное сравнение передаваемых данных АЗН для встречных ВС.
В работе [Ali et al., 2019] проанализированы характеристики барометрических и геометрических данных АЗН-В. Проведен анализ показателей качества (параметров целостности высоты), дополняющих данные о высоте ВС.
Таким образом, можно сделать вывод, что применительно к АЗН-В ведутся работы по исследованию параметров и применению данных о вертикальном положении ВС как в РФ, так и за рубежом. Проблема подтверждения данных АЗН-В, в том числе и в вертикальной плоскости, является актуальной.
Разработанная в статье методика основана на использовании данных геометрической и барометрической высот и параметрах реальной атмосферы при подтверждении данных АЗН-В. Основной вклад в погрешности вычислений вносит температура. В методике предлагается использовать температуру, полученную методом интерполяции прогнозных данных температуры, передаваемых метеослужбой в блоке прогнозов по площадям (GAMET). Для коррекции и повышения точности температуры предлагается корректировать прогнозные значения по данным метеозондов.
Дальнейшие исследования планируется направить на повышение точности определения температуры воздуха на высоте полета ВС. Другим важным направлением является определение допустимого отклонения высот доп.
при превышении которого диспетчер УВД должен получить уведомление о недостоверности получаемых данных АЗН-В. Также в дальнейшей работе планируется оценить точность разработанной методики перевода барометрической высоты ВС в геометрическую путем сравнения с данными многопозиционной системы наблюдения, как более точного средства наблюдения.
Заключение
В статье разработана методика подтверждения данных АЗН-В путем оценки барометрической и геометрической высот полета ВС. При выполнении расчетов для получения температуры использовались значения, полученные интерполяцией прогнозных значений с возможностью коррекции по данным
метеозонда. Основной вклад в погрешность определения высоты полета ВС для реальной атмосферы вносит температура. Для коррекции прогнозных значений температуры на высотах выполнения полетов предложено использовать данные метеозонда.
В качестве примера применения методики, в работе выполнены расчеты с использованием реальных данных о высотах полета ВС, получаемых от станции АЗН-В НС-1А производства АО «ВНИИРА», установленной на аэродроме Мезень, а также данные, получаемые от метеослужбы. Для эшелона 060 (6050 футов) была рассчитана разность высот, составившая 35,88 м.
При выполнении расчета было установлено, что данные АЗН-В с использованием предложенной методики подтверждаются, поскольку разность высот не превосходит принятого допустимого уровня 75 м.
Предложенная методика может быть применена для подтверждения данных геометрической высоты АЗН-В ВС при выполнении полетов МВЛ.
Библиографический список
Алипов И. В. Моделирование «тонкой» структуры вероятностного распределения невязок в показаниях АЗН-В и ВОРЛ. Некомпенсированная задержка / И. В. Алипов, В. Л. Кузнецов // Научный вестник Московского государственного технического университета гражданской авиации. 2019. Т. 22. №. 4. С. 8-20. DOI 10.26467/2079-0619-2019-22-4-8-20 Калинцев А. С. Подтверждение данных АЗН-В в аэродромной зоне методом стробирования /
A. С. Калинцев, Е. А. Рубцов, А. П. Плясовских // T-Comm: Телекоммуникации и транспорт. 2021. Т. 15. № 7. С. 39-49. DOI 10.36724/2072-8735-2021-15-7-39-49.
Косьянчук В. В. Обзор основных путей повышения безопасности системы АЗН-В /
B. В. Косьянчук, Н. И. Сельвесюк, Р. Р. Хамматов // Научный вестник Московского государственного технического университета гражданской авиации. 2019. Т. 22. №. 1. С. 3950. DOI 10.26467/2079-0619-2019-22-1-39-50
Лебедев Б. В. Исследование метода контроля вертикального эшелонирования воздушных судов на основе использования АЗН-В / Б. В. Лебедев, В. В. Соломенцев, А. Н. Стратиенко // Научный вестник Московского государственного технического университета гражданской авиации. 2015. №. 213 (3). С. 135-140.
Плясовских А. П. Метод оценки достоверности информации АЗН-в в системе наблюдения и контроля аэродромного движения / А. П. Плясовских, Е. А. Рубцов // Вестник Санкт-Петербургского государственного университета гражданской авиации. 2019. № 3(24). С. 90102.
Плясовских А. П. Теоретическое обоснование подтверждения достоверности информации о местоположении объекта на рабочей площади аэродрома / А. П. Плясовских, Е. А. Рубцов // T-Comm: Телекоммуникации и транспорт. 2020. Т. 14. № 3. С. 32-40. DOI 10.36724/2072-87352020-14-3-32-40.
Сравнение пропускной способности систем автоматического зависимого наблюдения-вещания / А. Мирошниченко, И. Татарчук, Э. Фальков [и др.] // Первая миля. 2020. №. 3(88).
C. 24-29. DOI: 10.22184/2070-8963.2020.88.3.24.29.
Aeronautical Surveillance Manual. Doc. 9924 AN/474. ICAO, 2020. 372 p.
Ali B. S. A study on geometric and barometric altitude data in automatic dependent surveillance broadcast (ADS-B) messages / B. S. Ali, N. A. Taib // The Journal of Navigation, 2019. V. 72. №. 5. pp. 1140-1158.
Calculation of Gravity Field Functionals on User-Defined Points // International Centre for Global Earth Models // [Электронный ресурс]. URL: http://icgem.gfz-potsdam.de/calcpoints (дата обращения: 29.01.2023).
Department of Atmospheric Science // University of Wyoming // [Электронный ресурс]. 2023. URL: http://weather.uwyo.edu/upperair/sounding.html (дата обращения: 29.01.2023). EGM96. The NASA GSFC and NIMA Joint Geopotential Model // The Onio state university // [Электронный ресурс]. URL: https://cddis.nasa.gov/926/egm96/egm96.html (дата обращения: 29.01.2023).
Evaluation of the EGM96 Model of the Geopotential in the United States // National Geodetic Survey // [Электронный ресурс]. URL: https://www.ngs.noaa.gov/PUBS_LIB/ egm96.html (дата обращения: 29.01.2023).
Global air navigation plan 2016-2030. Doc. 9750 AN/963, Fifth Edition. ICAO, 2016. 142 p. Manual of the ICAO standard atmosphere. Doc. 7488/3, 1993. 305 p.
Minimum operational performance standards for 1090 MHz extended squitter automatic dependent surveillance - broadcast (ADS-B) and traffic information services - broadcast (TIS-B). RTCA DO-260B, 2009. 1410 p.
Taib N. A. An Analysis of Geometric Altitude Data in ADS-B Messages / N. A. Taib, B. S. Ali // Proceedings of the 2016 International Technical Meeting of The Institute of Navigation, 2016. pp. 697-704.
Use International Standard Atmosphere model // Mathworks // [Электронный ресурс]. URL: https://www.mathworks.com/help/aerotbx/ug/atmosisa.html (дата обращения: 29.01.2023).
References
Aeronautical Surveillance Manual. Doc. 9924 AN/474. ICAO, 2020. 372 p.
AliB. S., Taib N. A. (2019). A study on geometric and barometric altitude data in automatic dependent surveillance broadcast (ADS-B) messages. The Journal of Navigation. 72(5): 1140-1158. Alipov I. V., Kuznetsov V. L. (2019). Modeling of the "thin" structure probability discrepancy distribution in the indications of ADS-B and SSR. Uncompensated latency. Civil Aviation High Technologies. 22(4): 8-20. (in Russian)
Calculation of Gravity Field Functionals on User-Defined Points // International Centre for Global Earth Models // [Electronic resource]. URL: http://icgem.gfz-potsdam.de/calcpoints (accessed 29.01.2023).
Department of Atmospheric Science // University of Wyoming // [Электронный ресурс]. URL: http://weather.uwyo.edu/upperair/sounding.html (accessed 29.01.2023).
EGM96. The NASA GSFC and NIMA Joint Geopotential Model // The Onio state university // [Electronic resource]. URL: https://cddis.nasa.gov/926/egm96/egm96.html (accessed 29.01.2023). Evaluation of the EGM96 Model of the Geopotential in the United States // National Geodetic Survey // [Electronic resource]. URL: https://www.ngs.noaa.gov/PUBS_LIB/egm96.html (accessed 29.01.2023).
Global air navigation plan 2016-2030. Doc. 9750 AN/963, Fifth Edition. ICAO, 2016. 142 p. Kalintsev A. S., Rubtsov E. A., Plyasovskih A. P. (2021). Confirmation of ADS-B data in the aerodrome traffic zone by gating method. T-Comm. 15(7): 39-49. DOI 10.36724/2072-8735-202115-7-39-49. (in Russian)
Kosianchuk V. V., SelvesiukN. I., KhammatovR. R. (2019). An overview of the main ways to improve the ADS-B system security. Civil Aviation High Technologies. 22(1): 39-50. (in Russian) Lebedev B. V., Solomentchev V. V., Stratienko A. N. (2015). Analysis method of control vertical separation of aircraft on the base of ADS-B. Civil Aviation High Technologies. 213(3): 135-140. (in Russian)
Manual of the ICAO standard atmosphere. Doc. 7488/3, 1993. 305 p.
Minimum operational performance standards for 1090 MHz extended squitter automatic dependent surveillance - broadcast (ADS-B) and traffic information services - broadcast (TIS-B). RTCA DO-260B, 2009. 1410 p.
Miroshnichenko A., Tatarchuk I., Falkov E. [and others]. (2020). Comparison of the throughput of automatic dependent surveillance-broadcast systems. First mile. 3: 24-29. DOI: 10.22184/20708963.2020.88.3.24.29. (in Russian)
Plyasovskih А. Р., Rubtsov Е. А. (2019). Reliability estimation method of ADS-B information for surface movement guidance and control system. Vestnik Saint Petersburg State University of Civil Aviation. 3(24): 90-102. (in Russian)
Plyasovskih А. Р., Rubtsov Е. А. (2020). Theoretical substantiation of confirmation of the validity of information about the location of the object on the work area of the aerodrome. T-Comm. 14(3): 3240. DOI 10.36724/2072-8735-2020-14-3-32-40. (in Russian)
Taib N. A., Ali B. S. (2016). An Analysis of Geometric Altitude Data in ADS-B Messages. Proceedings of the 2016 International Technical Meeting of The Institute of Navigation. 697-704. Use International Standard Atmosphere model // Mathworks // [Electronic resource]. URL: https://www.mathworks.com/help/aerotbx/ug/atmosisa.html (accessed 29.01.2023).