Статья: Методы нечеткой логики в георадарных задачах

Методы нечеткой логики в георадарных задачах

Сугак В.Г. 1, Михайлюк Е.А. 2, Дубовицкий А.Н. 1, Маматов Е.М. 3, Маматов М.Е. 3

1) ООО «Геоэксперт»

2) Белгородский государственный технологический университет им. В.Г. Шухова

3) Белгородский государственный национальный исследовательский университет

Аннотация

В статье представлен опыт обработки результатов георадарного зондирования мелового карьера методами нечеткой логики для выделения таких признаков сигналов, которые трудно осуществлять стандартными методами в силу «размытости» или нечеткости этих признаков. Данная методика позволяет выделить дополнительную полезную информацию, прогнозировать качество и более точно оценить состояние подповерхностной структуры участка, в данном случае, степень белизны мела. В основе предлагаемого сценария лежит разработка моделей и методов выделения параметров георадарных сигналов, в частности, их фазочастотных спектров, позволяющих установить причинно-следственные связи, влияющие на качество принятия решения при наличии совокупности мешающих факторов, на основе использования алгебры нечеткой логики.

Ключевые слова: георадар, фазочастотный спектр, меловой карьер, нечеткая логика, меловой карьер.

Abstract

Fuzzy logic methods in georadar applications

Vladimir G. Sugak 1, Ekaterina A. Mikhajlyuk 2, Alexander N. Dubovitskiy 1, Evgeny M. Mamatov 3, Mikhail E. Mamatov 3

1 «Geoexpert», LLC

2 Belgorod State Technological University named after V.G. Shukhov

3 Belgorod National Research University

A theoretical justification for the development of a scenario for the subsurface structure state estimation of a chalk quarry site based on the results of GPR sounding is proposed, which allows predicting its quality, namely, the degree of whiteness. The proposed scenario is based on the development of models and methods for identifying parameters of GPR signals, in particular, their pha se-frequency spectra, which make it possible to establish causal relationships that affect the quality of decision making in the presence of a combination of interfering factors. In this paper, for the first time, data on GPR sounding of a chalk quarry were obtained, showing the presence of stable signs for the selection of chalk areas with increased whiteness. Such a feature is the phase fluctuations of the reflected signals obtained at the cross-section with respect to the radiated polarization of the received signals. This requires a special antenna, in which the radiation of signals is carried out on the main polarization, and the reception is orthogonal to the radiated one. To obtain a clear conclusion about the probability of attributing the chalk layer on the studied area of the chalk quarry to the categories of high, medium or ordinary whiteness, it is proposed to use an algorithm based on the theory of fuzzy logic.

Key words: georadar, phase-frequency spectrum, chalk quarry, fuzzy logic, chalk pit.

Введение

Развитие информационных технологий на современном этапе позволяет извлекать существенно больше информации из данных измерений, полученных при проведении всевозможных изысканий в геофизике, чем это было ранее во времена, когда соответствующая аппаратура, методика проведения измерений и интерпретации результатов только разрабатывались.

Необходимость выделения дополнительной полезной информации из данных геофизических измерений диктует необходимость усложнять методы обработки получаемых сигналов для выделения в них таких признаков, обусловленных воздействием исследуемых объектов, которые трудно осуществлять стандартными методами в силу их «размытости» или нечеткости [Замятин, 2006].

В частности, в такой области геофизики, как георадарное зондирование, часто возникает необходимость оценивать возмущения параметров отраженных сигналов, вызванные воздействием исследуемых объектов, находящихся в подповерхностном пространстве, которые являются либо слабо выраженными, либо сопровождаются сопутствующими дополнительными факторами, связанными с воздействием различного рода аномалий структуры грунта [Золотарев, 2017, 2].

Рассмотрим указанные особенности обработки и выделения полезной информации на примере применения георадарного зондирования для обнаружения аномалий, связанных со степенью белизны природного карбоната кальция или мела, в меловом карьере «Мелстром», г. Белгород.

Использование мела в производстве в качестве наполнителя позволяет существенно снизить себестоимость готовой продукции и улучшить ее качество [Гурова, Коваленко]. Уровень белизны природного карбоната кальция, используемого за рубежом в качестве пигмента, составляет 90-95 % (по Эльрефо). Для наполнителя бумаги применяют природный карбонат кальция с белизной 85-95 % и выше. В бумажной промышленности используется измельченный природный карбонат кальция высокой степени дисперсности [Маскайкин, Федотов, Кирюшин, Масляев, 2019].

Поэтому контроль степени белизны мела при разработке меловых карьеров является очень актуальной и важной задачей.

Методика

нечеткая логика георадарный

В работе впервые оценивается возможность использования специализированного георадара для определения качества мела, а именно степени его белизны, в слое от поверхности выработки до глубины в несколько метров на примере мелового карьера «Мелстром», г. Белгород. Данный меловой карьер является типовым для многих регионов России [Сергеев, Сидорова, 1950, Мощанский, 1975, Sugak, Dubovitsky, 2018]. В качестве георадара использовался экспериментальный макет радиолокатора подповерхностного зондирования со ступенчатым изменением несущей частоты в диапазоне 150-300 мГц, разработанный компанией «Геоэксперт», г. Белгород [Sugak, Sugak, 2010]. Антенна георадара позволяет излучать сигналы с выбранной поляризацией и принимать поочередно сигналы на согласованной и ортогональной к излучаемой поляризациях. В указанной антенне использовано электронное переключение двух взаимно перпендикулярных щелей треугольной формы, прорезанных в металлическом корпусе с помощью электронного коммутатора, осуществляющего согласование волнового импеданса антенн с кабелем и переключение щелей антенны на прием или передачу. В результате антенна позволяет принимать и после обработки отраженных сигналов получать их полную поляризационную матрицу.

Так как в литературе отсутствует какая-либо информация об опыте применения георадара для решения указанной задачи, то методика изысканий и методы обработки сигналов разрабатывались практически с нуля.

Изыскания проводились на нескольких площадках мелового карьера «Мелстром», г. Белгород (рис. 1) размерами порядка 20х10 м, которые были отобраны его сотрудниками по результатам лабораторных измерений. Всего было выбрано 4 площадки (2 участка с мелом повышенной белизны, 1 участок с мелом средней белизны и 1 участок с обычным мелом), представленные на карте карьера (рис. 1).

Рис. 1. Карта мелового карьера

Fig. 1. Chalk pit map

Данные зондирования, полученные на указанных выше участках, сравнивались между собой и с данными, полученными в лаборатории мелового карьера по образцам, взятым в каждом исследуемом участке, для выявления характерных признаков в георадарных сигналах, которые могли бы в будущем лечь в основу технологии выявления мела повышенной белизны.

На каждой площадке зондирование проводилось в среднем по 4-м параллельным профилям движения георадара. В результате получались радиолокационные изображения сечения грунта по каждому из профилей. В дальнейшем для анализа использовались радиолокационные изображения фазовой структуры сигналов (рис. 1) по каждому из указанных профилей, так как только на них удавалось обнаружить признаки, позволяющие отличить участки с мелом повышенной белизны от других участков с нормальной и обычной белизной.

Оценка возможности применения специализированного георадара для определения степени белизны мела

Наибольшие отличительные признаки присутствия на исследуемом участке мела повышенной белизны выделяются на изображениях фазовой структуры отраженных сигналов (фазо-частотных спектров) - и только на ортогональной к излученной поляризации этих сигналов. Изображения фазовой структуры сигналов описывают поведение изолиний фазо-частотного спектра отраженных сигналов [Яхъяева, 2012], в которых ось ординат соответствует глубине от поверхности грунта, ось абсцисс - дистанции, пройденной георадаром при движении по поверхности вдоль выбранного направления, а ряд квазипараллельных кривых - изолиний фазовой структуры сигналов, соответствующих пересечению фазо-частотного спектра отраженных сигналов нулевого значения фазы. По расстоянию между изолиниями и их флуктуациями можно судить о физических свойствах среды зондирования и аномалий в ней [Яхъяева, 2012].

Рис. 2. Изображения фазовых структур на участке с мелом обычной белизны (левый рисунок) и с мелом повышенной белизны (правый рисунок)

Fig. 2. Images of phase structures in an area with chalk of normal whiteness (left figure) and with chalk of increased whiteness (right figure)

На основе полученных предварительных результатов, не вдаваясь в детали обработки сигналов, так как этот материал доступен в большей степени специалистам в области георадарного зондирования, удалось сделать следующие выводы: во-первых, наблюдаются отличия в фазовых структурах сигналов, проявляющиеся в том, что наличие фракций примесей материалов, отличных от чистого мела на участках с мелом обычной белизны, приводят к увеличению флуктуаций фазы сигналов (флуктуации глубины изолиний) в разных элементах объема мела по мере движения георадара вдоль выбранного профиля. Причем этот эффект проявляется на кроссовой к излученной поляризации отраженных сигналов, что требует применения специальных антенн, о которых сказано выше. Вывод иллюстрирует рис. 2, на котором приведены изолинии фазовой структуры отраженных сигналов (рисунки представляют изображения, полученные на 2 -х разных профилях). По мере роста глубины наблюдается увеличение флуктуаций фазы отраженных сигналов, соответствующих второму участку с мелом обычного качества. Ось абсцисс соответствует точкам зондирования вдоль профиля движения георадара (всего было пройдено расстояние около 15-20 м). Видно, что флуктуации фазы сигналов накапливаются с глубиной и отличия начинают проявляться уже на глубине около 2-3 м. Во-вторых, указанный эффект наблюдался практически по всем профилям движения георадара на выбранных участках, что позволяет характеризовать его как устойчивый признак, по которому можно обнаруживать и отличать участки карьера с мелом разной степени белизны. Среднеквадратичные значения флуктуаций фазы, соответствующие третьей по счету изолинии фазовой структуры сигналов от поверхности по всем профилям движения на участке с мелом повышенной белизны (левые гистограммы) и участке с обычной белизной (правые гистограммы), приведены в таблице 1. Видно, что наблюдаются существенно большие (на 2030 %) значения среднеквадратичных отклонений, соответствующих участку с мелом обычной белизны. Этот эффект можно использовать для разработки алгоритма автоматического обнаружения участков с мелом повышенной белизны без использования экспертных оценок, проводимых оператором. В таблице 2 приведены наибольшие среднеквадратичные значения флуктуаций фазы 4-х первых от поверхности изолиний, начиная с глубины 1,3 м. Видно, что в целом наблюдаются существенно большие флуктуации этой фазы на участках с мелом обычной и средней белизны. Однако существуют отдельные профили на участках мела с повышенной белизной, на которых этот вывод не подтверждается. Является ли это следствием того, что на этих профилях наблюдаются существенные отклонения указанных признаков от предполагаемых или на этих участках существуют только локальные зоны с мелом повышенной белизны, пока выяснить не удалось вследствие отсутствия априорной информации.

Таблица 1. Среднеквадратичные значений флуктуаций третьей изолинии фазовой структуры сигналов по глубине

Table 1. Root mean square values of fluctuations of the third isoline of the phase structure of signals over depth

В-третьих, наблюдается также стабильное небольшое, порядка 10-20 см, смещение по глубине изолиний фазовой структуры на кроссовой к излучаемой поляризации сигнала на участке с мелом средней и обычной белизны в сторону больших значений по сравнению с участком с мелом повышенной белизны. Этот результат также говорит в пользу модели, в которой предполагается, что в структуре мела обычной белизны присутствуют примеси, которые слегка увеличивают диэлектрическую проницаемость, что и приводит к смещению изолиний фазовой структуры по глубине на соответствующем участке. На рис. 2 показан этот эффект и указаны изолинии, на которых максимально наблюдается это смещение. Кроме того, видно также, что флуктуации фазы изолинии на участке с мелом обычной белизны больше соответствующих флуктуаций на участке с мелом повышенной белизны.

Таблица 2. Наибольшие среднеквадратичные значения флуктуаций фазы 4-х первых изолиний

Table 2. The largest root mean square values of phase fluctuations of the first 4 isolines

Таким образом, результаты зондирований хотя и показывают положительные тенденции в появлении таких характеристик сигналов, которые можно использовать для идентификации нужных аномалий структуры мела по глубине, но их проявления являются нечеткими. Выделять такие структурные признаки классическими методами обработки сигналов, основанными, например, на применении метода максимального правдоподобия, не представляется возможным.