СПОСОБ АНАЛИЗА ПОЛЕТНЫХ ДАННЫХ
METHOD OF FLIGHT DATA ANALYSIS

Область техники, к которой относится изобретение
Настоящее изобретение относится к способу анализа набора полетных данных, зарегистрированных в течение по меньшей мере одного полета по меньшей мере одного воздушного судна.
Уровень техники
Нормативные требования к обслуживанию техники и движению воздушных судов определяют стандарты, обеспечение соответствия которым авиакомпаниями необходимо для обеспечения максимального уровня безопасности пользователей.
Для оптимизации процедур обслуживания техники авиакомпании устанавливают системы анализа полетных данных.
Известны системы анализа полетных данных, такие как FDM (Flight Data Monitoring, контроль полетных данных) или же FOQA (Flight Operational Quality Assurance, операционный контроль качества полета). При использовании таких систем воздушное судно оборудуют регистратором полетных данных. Такой регистратор может представлять собой «черный ящик» или специализированный регистратор, например, типа ACMS (Aircraft Condition Monitoring System, система контроля состояния воздушного судна).
Такие системы позволяют авиакомпании восстанавливать подробную картину полета по регулярно регистрируемым значениям полетных данных, получаемым в ходе каждого полета каждого из ее воздушных судов.
Для этого такие системы выявляют заранее определенные события, происходящие в течение полета, после чего специалист анализирует такие события, указывающие на возникновение в ходе полета технических неисправностей или на несоблюдение правил или условий, определенных правилами полета, тем самым обеспечивая раннее предупреждение о возможных происшествиях или авариях.
Применение таких технологий требует заблаговременного определения правил выявления событий, которые обычно определяют как превышение пороговых значений одного или более параметров, что может вызвать срабатывание предупреждения, наличие которого позволяет специалисту выделить данный полет для более подробного анализа.
Недостаток таких технологий состоит в отсутствии возможности выявления единичных событий, не подпадающих под заранее определенные правила, что может привести к пропуску детектирования полета с отклонениями от нормы.
Раскрытие изобретения
Задача, на решение которой направлено настоящее изобретение, состоит в обеспечении возможности выявления полетов с отклонениями от нормы без необходимости задания правил выявления.
i
i
Для решения поставленной задачи в соответствии с изобретением предлагается способ анализа полетных данных, зарегистрированных в течение N полетов по меньшей мере одного воздушного судна, причем данные группируют по полету i в сигнатурный вектор X полета размерности d, компоненты которого соответствуют данным, зарегистрированным в течение указанного полета i воздушного судна, так что полет i задан указанным сигнатурным вектором X, причем способ включает этапы, на которых:
i
i
- производят Гауссовский покомпонентный анализ энтропии ядра полетных сигнатур X для получения области нормальных полетов и классифицируют полетные сигнатуры X по их расстоянию до указанной области;
i
i
- определяют для каждого полета i уровень z, отклонения от нормы, определенный расстоянием полетной сигнатуры X относительно области нормальных полетов;
i
- выявляют по меньшей мере один полет с отклонением от нормы в зависимости от уровня z отклонения от нормы.
В оптимальном варианте осуществления изобретения предусмотрены следующие дополнительные признаки, которые могут быть использованы по отдельности или в любом возможном техническом сочетании:
- для каждого выявленного полета с отклонением от нормы определяют модельный полет, ближайший к выявленному полету с отклонением от нормы и расположенный в пределах области нормальных полетов;
- параметры выявленного полета с отклонением от нормы сравнивают с параметрами определенного таким образом модельного полета с целью выявления по меньшей мере одного параметра полета с отклонением от нормы, который вызвал отклонение от нормы указанного выявленного полета с отклонением от нормы;
- Гауссовский покомпонентный анализ энтропии включает этапы, на которых:
i
- определяют матрицу К сходства размерности N×N, компоненты которой выражают в численном виде близость между двумя полетными сигнатурами X;
1
N
1
N
i
i
i
- производят разложение матрицы К сходства на собственные векторы, получая N собственных векторов а,…,a и N собственных значений λ,…,λ, таких, что ∀i=1,…,N K⋅a=λ⋅a;
i
- определяют для каждого собственного вектора соответствующий ему коэффициент энтропии Y;
m
i
- отбирают такое подмножество собственных векторов , в котором сумма значений энтропии Y больше, чем доля суммы N значений энтропии Y;
- коэффициент энтропии определяют по формуле
- уровень отклонения от нормы полета i определяют по формуле
- компоненты матрицы К сходства определяют по формуле
2
где σ - заранее определенный параметр сглаживания;
2
- параметр σ сглаживания определяют следующим образом:
- определяют матрицу D размерности N×N, соответствующую расстоянию между двумя полетами и заданную следующим образом: ;
- все компоненты каждого столбца располагают в порядке возрастания, получая матрицу D’;
- отбирают первые k строк полученной таким образом матрицы D’, а оставшиеся строки удаляют, в результате чего получают матрицу D’ размерности k×N;
1
N
- определяют среднее каждого из столбцов матрицы D’, получая таким образом N значений y,…,y;
1
N
1
N
- определяют абсолютное медианное отклонение значений y,...,y, определенное как , где med=median{y,…,y};
2
1
N
- определяют параметр σ сглаживания по абсолютному медианному отклонению значений y,…,y в соответствии с формулой .
i
i
i
i
- данные полета i группируют в матрицу F размерности Т×Р, где Т - количество элементов данных, зарегистрированных в течение полета i, а Р - количество зарегистрированных параметров, причем компоненты вектора X представляют собой столбцы матрицы F один за другим, так что вектор X имеет размерность d=T×P и определен формулой
.
i
i
i
i
- данные полета i группируют в матрицу F размерности Т×Р, где T - количество элементов данных, зарегистрированных в течение полета i, а Р - количество зарегистрированных параметров, причем компоненты вектора X представляют собой столбцы матрицы F, из которых отобраны n<T записей  параметров, так что вектор X имеет размерность d=n×Р и определен формулой
i
i
i
- данные полета i группируют в матрицу F размерности Т×Р, где T - количество элементов данных, зарегистрированных в течение полета i, а Р - количество зарегистрированных параметров, причем компоненты вектора X представляют собой среднее значение, дисперсию, минимальное значение и максимальное значение параметра по всем записям данного параметра, так что вектор X определен формулой
Решение по изобретению обладает многочисленными преимуществами.
Выявление полетов с отклонением от нормы в соответствии с изобретением может быть осуществлено автоматически, не требуя вмешательства специалиста для такого выявления.
В соответствии с изобретением осуществлено статистическое выявление с учетом параметров. Таким образом, происходит автоматическое обучение системы.
Изобретение обеспечивает возможность выявления неожиданных проблем.
Изобретение обеспечивает возможность выявления проблем, порожденных вкладами нескольких параметров.
Способ по изобретению может быть использован неспециалистом в области статистики или углубленного анализа данных.
Кроме того, настоящее изобретение позволяет упростить диагностику полетов с отклонением от нормы, так как обеспечивает возможность создания эталонного полета, модельного полета для каждого из выявленных полетов с отклонением от нормы.
Краткое описание чертежей
Другие особенности, задачи и преимущества настоящего изобретения станут ясны из нижеследующего описания, приведенного исключительно в иллюстративных целях и не накладывающего каких-либо ограничений, содержащего ссылки на прилагаемые чертежи. На чертежах:
- фиг. 1а и 1b схематически иллюстрируют этапы способа по одному из вариантов осуществления изобретения;
- фиг. 2 иллюстрирует представление полетных данных по одному из вариантов осуществления изобретения.
Аналогичные элементы обозначены на всех чертежах одинаковыми ссылочными номерами.
Осуществление изобретения
На предварительном этапе (подробно не описан) способа анализа полетных данных получают полетные данные, зарегистрированные в течение нескольких полетов по меньшей мере одного воздушного судна.
Такие полетные данные соответствуют параметрам полета, например, таким как скорость, высота, температура и т.д.
Затем, на первом этапе Е1 для каждого полета i следует извлечь из всех данных, зарегистрированных в течение полета i, данные, характеризующие полет i, или сигнатуру каждого полета i.
0
0
0
Для N≥1 полетов существует некоторое количество Р параметров, зарегистрированных в течение времени Т полета. Например, в некоторых регистраторах это количество может составлять Р=2000 параметров.
Подразумевается, что регистрация/выборка параметров производятся с одной и той же частотой; при необходимости к данным могут быть применены хорошо известные технологии передискретизации и линейной аппроксимации.
0
0
Для извлечения на этапе Е1 данных, характеризующих полет i, в соответствии с одним из вариантов осуществления выполняется ограничение до одной из стадий полета (например, стадии посадки) с ограничением по времени Т<Т и оставляется лишь ограниченное количество Р<Р релевантных параметров для анализа данных данного полета i. Например, могут быть установлены ограничения Т=1000 секунд и Р=30 параметров.
i
В нижеследующем описании значение каждого из элементов данных полета i обозначено следующим образом: F, (t=5, p=1), т.е. значение первого параметра, зарегистрированное на пятой секунде (если периодичность выборки измеряется в секундах) для полета i.
i
После извлечения на этапе Е1 данные группируются на этапе Е2 в сигнатурный вектор X полета i, имеющий размерность d.
i
Таким образом, каждый сигнатурный вектор X содержит набор релевантной информации, относящейся к полету i.
Группирование, производимое на этапе Е2, может быть осуществлено в соответствии с одним из трех вариантов.
В соответствии с первым вариантом группирование данных в вектор на этапе Е2 происходит с сохранением всех данных. Такой подход называют исчерпывающим.
i
i
В соответствии с данным вариантом осуществления для получения вектора X следует взять матрицу F и поставить ее столбцы один за другим.
i
Получаемый при этом вектор X имеет размерность d=Т×Р и определен следующим образом:
.
i
Например, при T=1000 и Р=30 получают вектор X размерности 30000 (а не матрицу размерности 1000×30).
Группирование в соответствии с описываемым первым вариантом осуществления может быть полезно, например, при отсутствии значительного количества подробной информации и позволяет легко выявить любые отклонения от маршрута полета.
i
В соответствии со вторым вариантом осуществления при группировании в вектор X на этапе Е2 происходит фокусирование на значениях параметров в конкретные моменты, так называемых моментальных снимков. Такие конкретные моменты соответствуют определенным рабочим операциям; например, на стадии посадки моментальные снимки могут быть произведены, начиная с высоты в 10000 футов через каждые 1000 футов снижения или же в моменты определенных событий - выпуска посадочного шасси и т.п.
Если определены n моментов моментальных снимков, то для полета i отбирают данные моментальных снимков . Следует отметить, что в разных полетах такие моменты могут не совпадать по времени.
i
В таком случае вектор X имеет размерность d=n×Р и определен следующим образом:
Второй вариант осуществления позволяет получить лучшие результаты, чем первый. Действительно, в описываемом втором варианте осуществления, поскольку T>>n, объем обрабатываемых данных значительно меньше, чем в первом варианте осуществления.
i
В соответствии с третьим вариантом осуществления группирование в вектор X на этапе Е2 состоит в «обобщении» каждой «кривой» каждого из параметров с использованием малого набора из k значений, где k ≈5 или 10: такие k значений могут, например, соответствовать среднему, дисперсии, максимальному значению, минимальному значению каждого из зарегистрированных параметров.
i
Такие k значений затем комбинируются для всех параметров с получением сигнатурного вектора X размерности d=k×Р при k≈4, определяемого по формуле
i
Например, если Р=30 параметров, сигнатурный вектор X имеет размерность d=120.
i
В конце этапа Е2 группирования как результат имеются N≥1 сигнатурных векторов X размерности d, компоненты которых соответствуют релевантным данным, зарегистрированным во время указанного полета/воздушного судна.
i
Дополнительно может быть предусмотрен необязательный этап Е3 нормализации каждого из компонентов каждого сигнатурного вектора X.
i
Такая нормализация может быть полезна, если компоненты каждого сигнатурного вектора X имеют величины разных порядков.
В результате нормализации на этапе Е3 в наборах данных для всех полетов предпочтительно получают компоненты с нулевым средним значением и стандартным отклонением, равным 1.
i
В нижеследующем описании символ X обозначает сигнатурный вектор полета i, компоненты которого могут быть или не быть нормализованы.
i
Для определения того, отклоняется ли полет от нормы, сигнатурные векторы X проходят несколько этапов обработки.
i
i
Для этого на этапе Е4 производят Гауссовский покомпонентный анализ энтропии ядра сигнатурных векторов X для определения области Е нормальных полетов и классифицируют сигнатурные векторы X по их расстоянию до указанной области.
i
Другими словами, все полетные сигнатуры X группируют так, чтобы демаркировать область Е нормальных полетов.
Анализ, производимый на этапе Е4, описан, например, в публикации R. Jenssen: "Kernel entropy component analysis", IEEE transactions on pattern analysis and machine intelligence, vol. 32, n°5, pages 847-860, May 2010.
i
Для осуществления такого анализа на этапе Е41 определяют матрицу К сходства размерности N×N, компоненты которой выражают в численном виде близость между двумя сигнатурными векторами X (т.е. между двумя полетами).
Матрица К известна под названием матрицы Грама.
2
Матрицу К получают из матрицы D расстояний, к которой применяют функцию Гаусса с заранее определенным параметром σ сглаживания (определение этого параметра описано ниже).
Матрица К определена следующим выражением:
i
2
Параметр сглаживания получают путем вычисления расстояний между каждым из сигнатурных векторов X и его k ближайшими соседями, где k - целое число в диапазоне от 4 до 10. Значение k зависит от числа N полетов и размерности d. Для определения параметра σ сглаживания на этапе Е42 определяют матрицу D размерности N×N, соответствующую расстоянию между двумя полетами и определенную следующим образом: .
Затем элементы каждого из столбцов матрицы D располагают в порядке возрастания, получая матрицу D’, такую, что .
Затем отбирают первые k строк полученной таким образом матрицы D’, а оставшиеся строки удаляют, в результате чего получают матрицу D’ размерности k×N.
1
N
После этого определяют среднее каждого из столбцов матрицы D’, получая таким образом N значений y,…,y, которые по существу выражают среднее расстояние между каждым из векторов и его к ближайшими соседями.
1
N
1
N
1
N
1
N
Затем удаляют аберрантные значения y,…,y. Для этого сначала определяют медианное значение med набора y,…,y, которое выражают в виде med=median{y,…,y}. Затем определяют абсолютное медианное отклонение mad значений y,…,y, определенное как медианное значение набора , записываемое в виде mad=median .
2
1
N
i
Тогда параметр σ сглаживания можно определить по абсолютному медианному отклонению значений y,…,y в соответствии с формулой . Другими словами, все значения y, превышающие med+20⋅mad, исключают из рассмотрения, а наибольшее из оставшихся значений принимают за параметр сглаживания.
1
N
1
N
j
i
i
После получения матрицы К на этапе Е43 осуществляют разложение матрицы К сходства на собственные векторы, получая при этом N собственных векторов a,…,a и N собственных значений λ,…,λ, таких, что ∀i=1,…,N K⋅a=λ⋅a.
i
i
i,j
i
i
N
Следует отметить, что каждый из собственных векторов a имеет размерность N:a∈ℜ, а символ a обозначает j-й компонент вектора a. Величины λ также называют коэффициентами энергии.
i
По полученным векторам и собственным значениям на этапе Е44 для каждого собственного вектора определяют его коэффициент энтропии Y по следующей формуле
Коэффициенты энтропии представляют собой критерий релевантности, позволяющий отобрать из набора сигнатурных векторов всех полетов лишь наиболее релевантные данные.
m
i
А именно на этапе Е45 отбирают такое подмножество собственных векторов , сумма соответствующих значений энтропии Y больше, чем доля суммы N значений энтропии Y.
Такая доля предпочтительно составляет от 75 до 95% и предпочтительно равна 90%.
i
Отбор такого подмножества позволяет, возвращаясь к параметрам каждого из векторов X, соответствующих отобранным собственным векторам, определить область Е нормальных полетов.
На фиг. 2 представлен простой пример, в котором использованы два параметра, что позволяет представить все полеты на двумерном графике и определить область Е нормальных полетов в виде окружности, окружающей кластер данных.
i
Затем на этапе Е5 для каждого полета i определяют уровень z, отклонения от нормы, определяемый расстоянием полета X относительно области Е нормальных полетов.
Уровень отклонения от нормы определяют по формуле ,
k
k
где a - отобранные собственные векторы, а λ - соответствующие им собственные значения.
Уровень отклонения от нормы имеет значения в диапазоне от 0 до 1. Чем ближе его значение к 1, тем более отклоняющимся от нормы считают данный полет.
Представленные на фиг. 2 полеты, расположенные вдали от области Е нормальных полетов, имеют высокий уровень отклонения от нормы.
В соответствии с предпочтительным вариантом осуществления изобретения полеты i, уровень отклонения от нормы которых превышает 0,99, являются отклоняющимися от нормы с высокой вероятностью, а полеты i, уровень аномальности которых превышает 0,999, считают отклоняющимися от нормы практически наверняка.
i
Затем на этапе Е6 в зависимости от уровня z отклонения от нормы определяют, является ли по меньшей мере один полет отклоняющимся от нормы.
Уровень отклонения от нормы может быть отображен для рассмотрения аналитиком.
Следует отметить, что для определения того, отклоняется ли полет от нормы, в соответствии с данным способом критерии, соответствие которым определяет нормальный полет, а несоответствие которым, соответственно, свидетельствует об отклонении от нормы, определяют с использованием имеющихся данных. Таким образом, необходимость в предварительном определении правил выявления отсутствует.
В случае выявления полета с отклонением от нормы для него на этапе Е7 определяют условный эталонный полет, называемый «модельным полетом», ближайший к данному полету с отклонением от нормы, но находящийся при этом в области нормальных полетов, причем модельный полет может быть отличным от сигнатурного полета.
Другими словами, модельный полет определяют как теоретический полет, расположенный в области нормальных полетов, ближайший к выявленному полету с отклонением от нормы.
0
В частности, определяют сигнатуру Z модельного полета.
Уровень z отклонения от нормы для модельного полета составляет менее 0,99 (или менее 0,95, если необходимо гарантированно обеспечить нормальность такого полета).
Таким образом, модельный полет представляет собой проекцию выявленного полета с отклонением от нормы на область нормальных полетов. В частности, он представляет собой ближайшую такую проекцию.
Модельный полет, сформированный в соответствии со способом по изобретению, потенциально не зависит от проанализированных полетов.
0
i
Модельный полет с сигнатурой Z, соответствующий выявленному полету с отклонением от нормы с сигнатурой X, формируют при помощи процедуры нелинейной условной оптимизации. Такая процедура оптимизации для заданного постоянного порогового значения seuil может быть выражена следующей формулой:
Следует отметить, что вышеописанная нормализация на единицу позволяет обеспечить минимизацию числа компонентов, различных в данном полете с отклонением от нормы и соответствующем ему модельном полете.
0
Процедура оптимизации требует вычисления величины уровня z(Z) для получения сигнатуры Z модельного полета.
d
z
Для этого для всех Z∈ℜ определяют следующий вектор k:
d
Тогда функция z, определяющая уровень любого вектора, принадлежащего ℜ, определена следующим образом:
1
m
1
m
где a,…,a - m ранее определенных собственных векторов (которые охватывают 90% суммарной энтропии), а λ,…,λ - их соответствующие собственные значения.
Если полет, соответствующий параметрам, обозначенным на фиг. 2 как 20, отклоняется от нормы, то ему соответствует модельный полет, параметры которого обозначены как 20’.
Это позволяет сравнить на этапе Е8 параметры выявленного полета с отклонением от нормы с параметрами определенного таким образом модельного полета с целью выявления по меньшей мере одного параметра полета с отклонением от нормы, который вызвал отклонение от нормы указанного выявленного полета с отклонением от нормы.