Для оценки состояния сложных промышленных объектов активно используют мониторинг – постоянное наблюдение за техническим состоянием конструкции или агрегата. Развитие и распространение систем мониторинга связано со старением оборудования и необходимостью продления срока его эксплуатации, а их успешное функционирование обеспечивается высоким техническим уровнем средств неразрушающего контроля, развитием их алгоритмического и программного обеспечения [1]. Возможность отслеживать возникновение и развитие трещин, разломов и иных дефектов позволяет планировать ремонтные работы или профилактическое обслуживание, а также предотвращать аварийные ситуации.
Наиболее эффективным для мониторинга особо опасных объектов является метод акустической эмиссии (АЭ), под которой понимается явление возникновения и распространения упругих колебаний (акустических волн) в различных процессах, например, при деформа- ции напряженного материала, истечении газов, жидкостей, горении, взрыве и др. [1]. Количественно АЭ – критерий целостности материала, который определяется звуковым излучением материала при его контрольной нагрузке. Эффект АЭ может использоваться для определения образования дефектов на начальной стадии разрушения конструкции и слежения за характером образования и развития дефектов в материале всего объекта в целом.
Основа АЭ-мониторинга – обнаружение и превращение упругих волн в электрический сигнал. Анализ этих сигналов дает ценную информацию о наличии и происхождении дефектов в материале [2]. Одной из проблем, возникающих при проведении АЭ-мониторинга в задаче диагностики состояния контролируемого объекта, является большой объем разнородной диагностической информации: сигналы АЭ, их параметры, показания датчиков температуры, давления, влажности и др. [3]. Помимо этого, при проведении мониторинга контроли- руемый объект работает в штатном режиме, но вибрация оборудования, течение продукта (например, жидкость, нефть и др.) в трубопроводе, механические воздействия (например, ослабленные болты или незакрепленные разъемы, которые перемещаются при воздействии ветра, дождь, летящие объекты, пыль) и другие воздействия создают акустические помехи, которые несут в себе заведомо ложную информацию, что затрудняет проведение мониторинга. Для исключения сигналов-помех необходимо решать проблемы идентификации сигналов и фильтрации данных.
Некоторые реальные подходы компенсации влияния фоновой помехи включают в себя изготовление специальных датчиков с электронными фильтрами для блокирования шумов [4], принимая во внимание расстановку датчиков как можно дальше от источника шума и электронную фильтрацию (с помощью либо учета времени появления сигнала, либо различия в спектральном составе истинных сигналов АЭ и фоновых шумов) [3]. Однако даже применение указанных подходов не всегда и не полностью исключает появление шумовых сигналов в данных, предоставляемых эксперту для дальнейшей диагностики. Помимо этого, эксперту зачастую необходимо выбирать из большого числа осциллограмм определенные виды сигналов для анализа дефекта на рассматриваемом объекте.
Для помощи эксперту в решении задачи идентификации сигналов могут применяться методы искусственного интеллекта, в частности, прецедентный подход. Каждый сигнал можно рассматривать как отдельный прецедент, и тогда, имея информацию о ранее выявленных сигналах, можно без активного привлечения эксперта идентифицировать определенные сигналы для их последующего анализа экспертом [5].
Представление информации о сигналах АЭ в виде прецедентов
В прецедент предлагается включить описание ситуации, решение для данной ситуации и информацию о результате применения решения: CASE = (Situation, Solution, Result), где Situation – ситуация, к которой применяется данный прецедент; Solution – решение (например, диагноз и рекомендации пользователю); Result – результат применения решения, который может включать список выполненных действий, дополнительные комментарии и ссылки на другие прецеденты, а при необходимости обоснование выбора данного решения и возможные альтернативы [6].
Полученный сигнал содержит достаточно много данных, являющихся промежуточными и не несущих в себе важную информацию. Поэтому отображение сигнала можно упрощенно представить в виде ломаной из экстремумов осциллограммы. Экстремумом будет считаться точка с максимальной/минимальной амплитудой в ε-окрестности, которая определяется длиной сигнала.
По упрощенной версии сигнала строится последовательность возрастаний и спадов сигнала с фиксированием значений относительной амплитуды и диапазона времени. Другими словами, строится последовательность r1, r2, …, rK, каждый член которой ri представлен четверкой {T, δA, tн, tк}, где T – тип (рост или спад), δA ∈ [–1; 1] – относительная амплитуда на этом участке, tн, tк ∈ [0; 1] – относительные временные границы диапазона времени.
Далее полученная последовательность проходит процедуру группировки по следующим правилам:
- пара возрастания и спада сигнала, действующая в пределах £ e, считается пиком;
- пара возрастания и спада сигнала, действующая в пределах > e, считается колебанием;
- идущие подряд N > 1 колебаний с близкими значениями относительной амплитуды группируются (объединяются в одно с суммарным диапазоном времени).
После группировки образуется аналогичная последовательность S1, S2, …, SL, где Si представляется такой же четверкой {T, δA, tн, tк}, но только T ∈ {пик; колебания}.
Формирование данной последовательности позволяет значительно сократить информацию о сигнале, что, в свою очередь, ведет к существенному уменьшению временных затрат на поиск похожих сигналов в базе знаний (БЗ).
В качестве иллюстрации приведем следующий пример: для сигнала (рис. 1а) длиной в 2 000 точек была получена последовательность из 30 четверок (рис. 1б). В графическом отображении сформированной последовательности (рис. 1в) серым цветом показаны колебания, черным – пики.
Полученную последовательность можно сохранить в БЗ как прецедент для дальнейшего Рис. 1. Пример преобразования данных осциллограммы в прецедент Fig. 1. An example of converting oscillogram data into a use case сравнения с имеющимися прецедентами в БЗ, в которой могут храниться прецеденты как с сигналами-помехами, так и с актуальными сигналами.
Заметим, что прецедент необязательно содержит полную информацию о сигнале (например, могут отсутствовать конкретные данные осциллограммы сигнала), но он должен содержать информацию, необходимую для экспертной оценки актуальности сигнала [7].
Прецедентный подход для идентификации сигналов АЭ
CBR-методы поиска решения включают четыре основных этапа, образующих так называемый CBR-цикл: извлечение, повторное использование, адаптация и сохранение прецедента [8]. Существует целый ряд методов извлечения прецедентов, например, метод ближайшего соседа, на основе деревьев решений, метод извлечения с учетом применимости прецедента и др. [6].
Для определения сходства текущего сигнала с прецедентами из БЗ предлагается использовать метод ближайшего соседа (nearest neighbor) [9] – наиболее используемый метод сравнения и извлечения прецедентов. Он позволяет достаточно просто вычислить степень сходства текущей проблемной ситуации и пре- цедентов из БЗ. Для определения степени сход- ства на множестве параметров, используемых для описания прецедентов и текущей ситуации, вводится некоторая метрика. Далее в соответствии с выбранной метрикой определяется расстояние от целевой точки, соответствующей текущей проблемной ситуации, до точек, представляющих прецеденты из БЗ, и выбирается точка, ближайшая к целевой.
Эффективность метода ближайшего соседа во многом зависит от выбора метрики (меры сходства). Две последовательности могут различаться как самими членами, так и их количеством, поэтому для сравнения последовательностей колебаний выбрана метрика на основе расстояния Левенштейна [9].
Расстояние Левенштейна для последовательностей S1 и S2 (длиной M и N соответственно) можно подсчитать по формуле d(S1, S2) = D(M, N), где
где S1, S2 – входные последовательности; m(a, b) = 0, если a = b, в противном случае m(a, b) = 1.
Между новым прецедентом T и сохранен- ными в базе прецедентами C в выбранной мет- рике определяются расстояния dCT. Для опреде- ления значения степени сходства Sim(C, T) необходимо найти максимальное расстояние dmax в выбранной метрике, используя границы диапазонов соответствующих параметров [7].
Оценка по методу ближайшего соседа: , где dCT – расстояние между текущим сигналом и прецедентом из БП; dmax – максимальное расстояние в выбранной метрике.
При определении класса сигнала вместо конкретных чисел эксперт сравнивает образы сигналов, в частности, наличие колебаний и их амплитуду. Предложенный подход к построению последовательности колебаний для нового прецедента и поиск похожего прецедента в БЗ позволяют ускорить и упростить работу эксперта при идентификации сигналов во время анализа файла данных осциллограмм.
Архитектура прецедентного модуля для идентификации сигналов АЭ
Разработанный CBR-модуль (рис. 2) является надстройкой, которую в виде отдельной библиотеки можно подключать к програм- мным системам соответствующего назначения, реализованным на различных языках про- граммирования. Модуль может получать от ос- новной программы буфер с данными осцилло- граммы, который далее поступает в блок извлечения осциллограмм и проходит стадии CBR-цикла.
Из данных осциллограммы формируется прецедент, далее полученный прецедент проходит стадии извлечения и сохранения в БЗ. Сохраненные прецеденты проходят процедуру маркировки в зависимости от заданных пользователем классов. Однако с ростом числа сохраненных прецедентов в БЗ появляется проблема снижения быстродействия системы [10], для решения которой применяется блок сокращения числа прецедентов.
CBR-модуль решает задачу маркировки и фильтрации файла данных осциллограмм для облегчения и ускорения работы эксперта (оператора). После анализа содержимого файла данных с помощью интерфейса модуль пе- редает основной программе номера интересующих (маркированных положительно) осциллограмм, что позволяет оператору быстро перемещаться по файлу с данными АЭ-мониторинга.
Взаимодействие между разработанным модулем и основной программой
Разработан набор сообщений, с помощью которых идет общение между основной про- граммой АЭ-мониторинга и реализованной надстройкой. Для начала работы CBR-модуля необходимо загрузить конфигурацию с подобранными и помеченными экспертом осциллограммами. После этого модуль готов к приему сообщений с осциллограммами для идентификации. Сообщения должны содержать в себе следующую информацию: размер буфера (в байтах), размер единицы данных (в байтах), номер осциллограммы, начало блока точек осциллограммы (опционально), количество точек осциллограммы (опционально).
Если последние два пункта содержат в себе нулевое значение, то из буфера извлекается вся осциллограмма, в противном случае извлекается указанный набор точек. Извлеченный набор сохраняется как массив вещественных чисел (double) и проходит обработку прецедента, реализуя CBR-цикл. Соответствующий алгоритм представлен на рисунке 3. После обработки в БЗ CBR-модуля сохраняются сформированная последовательность, результат маркировки и полученный номер осциллограммы.
Модуль может принимать сообщения, в которых будут указаны классы осциллограмм, требующие просмотра и изучения оператором. В этом случае происходит настройка приоритетных для отображения классов маркированных осциллограмм.
Номера осциллограммы используются для быстрого перемещения по имеющимся в основной программе данным. Необходимо отправить сообщение, в котором будут содержаться текущий номер осциллограммы основного проекта и направление поиска ближайшей подходящей по маркировке осциллограммы. Таким образом, пользователь не ограничен в перемещении по файлу данных.
В случае, когда пользователя не устраивает результат, полученный CBR-модулем, можно отправить сообщение об изменении класса маркировки.
Программная реализация прецедентного модуля для идентификации сигналов АЭ
Предложенный прецедентный подход для идентификации сигналов АЭ на основе метода ближайшего соседа с использованием метрики Левенштейна был реализован в программном модуле для системы анализа файлов, содержащих данные АЭ-мониторинга (рис. 4), на языке C# (.NET 4.5) [11] в среде Microsoft Visual Studio 2012 под операционную систему MS Windows.
Пользователь задает настройки работы CBR-модуля: устанавливает виды уже из-вестных (классифицированных экспертом) сигналов и выбирает необходимые для отоб-ражения классы сигналов. После этого CBR-модуль выполняет процедуру маркировки файла данных и предоставляет пользовате-лю упрощенное перемещение по содержи-мому файла, позволяющее просматривать только заданные пользователем классы сигналов. При желании можно совмещать обычный и упрощенный просмотры для детального анализа файла данных.
Работа модуля была проверена экспертами на тестовых данных, полученных при АЭ-мониторинге металлических конструкций (более 20 файлов измерений), результаты подтвердили эффективность применения прецедентного подхода для идентификации сигналов, получаемых во время АЭ-мониторинга сложных технических объектов. Применение реализованного CBR-модуля позволяет значительно сократить время анализа данных благодаря полученной на основе прецедентов маркировке содержимого файла (более 106 осциллограмм).
Заключение
В работе исследована возможность применения методов поиска решения на основе прецедентов для идентификации сигналов АЭ-мониторинга сложных технических объектов. Для извлечения прецедентов из БЗ (базы прецедентов) предложено использовать метод ближайшего соседа и расстояние Левенштейна в качестве меры сходства. Разработан метод представления сигналов АЭ-мониторинга как прецедента. На основе предложенного подхода реализован прецедентный модуль для идентификации сигналов АЭ-мониторинга. Работа модуля была проверена на тестовых данных, полученных при АЭ-мониторинге металлических конструкций. Полученные результаты по- казали эффективность предложенного подхода и его реализации в виде CBR-модуля.
Работа выполнена при финансовой поддержке проектов РФФИ №№ 18-01-00459, 18-29-03088, 18-51-00007.
Литература
1. Барат В.А., Алякритский А.Л. Метод статистической обработки данных акустико-эмиссионного мониторинга на примере реактора гидроочистки Мозырского НП // В мире НК. 2008. № 4. С. 52–55.
2. Терентьев Д.А. Идентификация сигналов акустической эмиссии при помощи частотно-временного анализа // В мире НК. 2013. № 2. С. 51–55.
3. Kwon O.-Y., Ono K. Acoustic emission characterization of the deformation and fracture of an sic-reinforced aluminum matrix composite. JAE, 1990, vol. 9, pp. 123–129.
4. Кузьмин А.Н., Журавлев Д.Б., Филип- пов С.Ю. Коррозия – приговор или диагноз? К вопросу технической диагностики тепловых сетей // Технадзор. 2009. № 3. С. 32–35.
5. Alekhin R.V., Varshavsky P.R., Ereme- ev A.P., Kozhevnikov A.V. Application of the case-based reasoning approach for identification of acoustic-emission control signals of complex technical objects. Proc. 3rd Russian-Pacific Conf. on Computer Technology and Applications (RPC), 2018, pp. 28–31. DOI: 10.1109/RPC.2018.8482124.
6. Еремеев А.П., Варшавский П.Р., Куриленко И.Е. Моделирование временных зависимостей в интеллектуальных системах поддержки принятия решений на основе прецедентов // Information Technologies and Knowledge. 2012. Т. 6. № 3. С. 227–239 (рус.).
7. Варшавский П.Р., Алехин Р.В., Кожевников А.В. Разработка прецедентного модуля для идентификации сигналов при акустико-эмисси- онном контроле сложных технических объектов // КИИ-2018: тр. XVI Национальн. конф. 2018. Т. 2. С. 222–229.
8. Aamodt A., Plaza E. Case-based reasoning: foundational issues, methodological variations, and system approaches. Artificial Intelligence Communications. 1994, vol. 7, no. 1, pp. 39–59.
9. Рассел С., Норвиг П. Искусственный интеллект: современный подход; [пер. с англ.]. М.: Вильямс, 2007. 1408 с.
10. Варшавский П.Р., Ар Кар Мьо, Шунке- вич Д.В. Применение методов классификации и кластеризации для повышения эффективности работы прецедентных систем // Программные продукты и системы. 2017. № 4. С. 625–631.
11. Albahari J., Albahari B. C# 4.0 in a Nutshell. The Definitive Reference. O'Reilly Media, Inc., 2010, 1056 p.
References
- Barat V.A., Alyakritski A.L. Statistical data processing of acoustic emission monitoring: Mozyrski NPZ hydrotreater case study. NDT World. 2008, no. 4, pp. 52–55 (in Russ.).
- Terentev D.A. Time frequency analysis for identification of acoustic emission signals. NDT World. 2013, no. 2, pp. 51–55 (in Russ.).
- Kwon O.-Y., Ono K. Acoustic emission characterization of the deformation and fracture of an sic-reinforced aluminum matrix composite. JAE. 1990, vol. 9, pp. 123–129.
- Kuzmin A.N., Zhuravlev D.B., Filippov S.Yu. Corrosion – sentence or diagnosis? On the issue of technical diagnostics of heat networks. Engineering Supervision. 2009, no. 3, pp. 32–35 (in Russ.).
- Alekhin R.V., Varshavsky P.R., Eremeev A.P., Kozhevnikov A.V. Application of the case-based reasoning approach for identification of acoustic-emission control signals of complex technical objects. Proc. 3rd Russ.-Pacific Conf. on Computer Technology and Applications (RPC). 2018, pp. 28–31 (in Russ.). DOI: 10.1109/RPC.2018.8482124.
- Eremeev A.P., Varshavsky P.R., Kurilenko I.E. Modeling the time-dependent in intellectual systems of decision making support on the basis of precedents. Information Technologies & Knowledge. 2012, vol. 6, no. 3, pp. 227–239 (in Russ.).
- Varshavsky P.R., Alekhin R.V., Kozhevnikov A.V. Development of the precedent module for identification of signals with acoustic-emission control of complex technical objects. Proc. 16th Russ. Conf. (RCAI 2018). Moscow, 2018, no. 2, pp. 222–229 (in Russ.).
- Aamodt A., Plaza E. Case-based reasoning: foundational issues, methodological variations, and system approaches. Artificial Intelligence Communications. IOS Press, 1994, vol. 7, no. 1, pp. 39–59.
- Russell S., Norvig P. Artificial Intelligence: Modern Approach. Prentice Hall Publ., 2 ed., 2002, 1132 p. (Russ. ed.: 2nd ed., Moscow, Villiams Publ., 2006, 1408 p.).
- Varshavsky P.R., Ar Kar Myo, Shunkevich D.V. Classification and clustering methods for improving efficiency of case-based systems. Software & Systems. 2017, no. 4, pp. 625–631 (in Russ.).
- Albahari J., Albahari B. C# 4.0 in a Nutshell. The Definitive Reference. O'Reilly Media, 2010, 1056 p.