Промышленное освоение месторождений криолитозоны показало необходимость контроля влажности горных пород, которая является весьма важной характеристикой, т.к. влияет на основные физико-механические свойства пород. Знание распределения влажности в горных породах криолитозоны необходимо как при эксплуатационной разведке месторождений, так и при контроле устойчивости грунтов оснований горнотехнических сооружений. В работах [1–3] обоснована и предложена методика георадиолокационного мониторинга криогенного состояния горных пород, позволяющая дистанционно определять влажность пород по результатам анализа данных разносезонных георадиолокационных измерений. Методика ориентирована на изучение влажности горизонтально-слоистых сред и не учитывает наличие аномалий георадиолокационных разрезов, связанных с развитием криогенных процессов в грунтах. Выявление подобных аномалий возможно с помощью специализированных алгоритмов обработки данных георадиолокации.
Для создания соответствующего программно-методического комплекса георадиолокационного мониторинга влажности горных пород необходима разработка алгоритмического обеспечения, включающего модели неоднородностей мерзлых горных пород, связанных с криогенными процессами, алгоритмы их выявления, а также процедуры для автоматизированной обработки, хранения и визуализации данных о выявленных неоднородностях.
Аномалии георадиолокационных разрезов, связанные с криогенными процессами
При разработке алгоритмического обеспечения необходимо основываться на моделях геологических разрезов, включающих изучаемые криогенные процессы (рис. 1), и особенностях распространения в них электромагнитных волн. Одним из основных электрофизических свойств горных пород является комплексная диэлектрическая проницаемость (КДП): . Мнимая часть
характеризует поглощение электромагнитной энергии в среде, а действительная часть определяет скорость распространения электромагнитных волн в среде и их коэффициент отражения (k) на границе двух сред:
.
Опыт проведения георадиолокационных исследований показывает, что основными аномалиями георадиолокационных разрезов, связанными с криогенными процессами, являются (рис. 2):
а) переотражение сигналов на участках повышенной влажности;
б) высокие значения амплитуд сигналов в зонах пучения пород;
в) скачкообразное изменение времени задержки сигналов, характерное для пустот;
г) разрывы осей синфазности сигналов с изменением времени задержки сигналов в зонах просадочных деформаций;
д) «хаотичные» сигналы на участках разуплотненных пород;
е) оси синфазности высокоамплитудных сигналов, характерные для границ сброса поверхностных вод или миграции надмерзлотных вод сезонно-талого слоя.
При исследовании участков повышенной влажности методом георадиолокации на границе мерзлых () и талых увлажненных (
) горных пород коэффициент отражения составляет около 0,4 и, соответственно, 40 % энергии отражается вверх. Далее происходит цикличное переотражение электромагнитной волны в мерзлом слое горных пород, которое формирует на радарограмме кратные по времени регистрации сигналы (рис. 2, а).
Рис. 1. Модели геологических разрезов, формирующие аномальные волновые картины в данных георадиолокации
Рис. 2. Фрагменты георадиолокационных разрезов, пересекающих зоны развития криогенных процессов
Зоны пучинистых горных пород в мерзлом состоянии на радарограмме представляются в виде линзообразных объектов (рис. 2, б). При георадиолокационном зондировании пустот оси синфазности сигналов, отраженных от нижней границы полости, имеют локальный характер и аномально малые значения времени регистрации (сигналы «прыгают» вверх) из-за высокой скорости распространения электромагнитных волн в воздухе (рис. 2, в). Зона просадки горных пород состоит из горизонтального участка просадки поверхности пород, в пределах которого просадка достигает максимальной величины, и двух боковых криволинейных участков (в виде круто наклонных осей синфазности), на которых просадка изменяется от максимальной величины до нуля (рис. 2, г, д). Зоны разуплотнения образуются в результате выноса водным потоком части горных пород определенной фракции, что приводит к образованию искусственной трещиноватости и нарушениям текстуры пород.
На радарограмме зоны разуплотнения отображаются в виде «хаотичных» сигналов, не имеющих выраженных осей синфазности (рис. 2, е). Миграция надмерзлотных вод сезонно-талого слоя происходит по кровле многолетнемерзлых пород, представляющих собой водоупор. Коэффициент отражения электромагнитных волн (k) на границе обводненных и мерзлых пород очень высок, вследствие чего георадиолокационные сигналы, отраженные от этой границы, имеют высокие амплитудные значения (рис. 2, ж). В данном случае кратного переотражения сигналов не происходит из-за значительного поглощения электромагнитной энергии в талых породах ( [4]).
Статистические методы оценки данных георадиолокационных исследований
Для автоматизированной обработки георадиолокационных данных возможно использование таких статистических методов, как дисперсионный и корреляционный анализ. Опыт георадиолокационных зондирований неоднородного массива мерзлых пород и натурные эксперименты показали, что локальные изменения свойств пород проявляются на георадиолокационном разрезе зонами с хаотичным распределением осей синфазности сигналов [5]. Появление подобных сигналов приводит к увеличению такой статистической характеристики, как дисперсия (DХ), которая определяет отклонение амплитуд сигналов от среднего ожидаемого значения. Анализ особенностей характера изменения значения дисперсии амплитуд сигналов по разрезу позволяет определять тип геологической неоднородности: относительно сглаженное изменение дисперсии определяет слоистую неоднородную структуру, резкие изменения дисперсии соответствуют нарушенной структуре с многочисленными локальными неоднородностями [6].
Корреляционный анализ является оптимальным методом обнаружения сигнала известной формы на фоне белого шума [7]. Данный метод хорошо зарекомендовал себя при прослеживании осей синфазности георадиолокационных сигналов. При обработке данных, полученных в непрерывном режиме съемки, коэффициент корреляции соседних сигналов, отраженных от одной геологической границы, составляет 0,95–0,99, что позволяет с высокой точностью прослеживать отражающие границы даже в случаях резкого изменения положения на временной оси [8].
Перечисленные выше аномалии георадиолокационных разрезов, связанные с криогенными процессами, могут быть картированы в автоматическом режиме на основе использования указанных статистических методов. В зонах с переотражениями сигналов в георадиолокационных трассах присутствуют сигналы одинаковой формы, на кратных временах задержки и с убывающими амплитудными значениями (закономерность убывания зависит от электрофизических свойств исследуемых горных пород). В этом случае в результате вычисления коэффициента корреляции искомого сигнала-образца по всей длине анализируемой трассы будет присутствовать несколько локальных максимумов, проверка времени задержки которых позволит выявить их кратность и картировать зоны повышенной влажности горных пород. Если в результатах расчета коэффициента корреляции количество максимумов равно 2, то после проведения опережающей проверки последующих трасс (для исключения случайных выбросов) возможно прослеживание обеих осей синфазности. После анализа радиофизических характеристик сигналов временной промежуток между выявленными границами может быть отмечен как пучинистые породы. В случае, когда после расчета коэффициента корреляции количество максимумов равно 1 и время задержки максимума значительно отличается, то после проведения опережающей проверки последующих трасс до возвращения границы на прежний уровень, этот участок следует отметить как пустоту.
Просадочные деформации горных пород, проявляющиеся на радарограммах в виде круто наклонных границ, могут быть выявлены посредством вычитания из каждой трассы разреза соседней трассы (данные берутся по абсолютным значениям). Это исключит из георадиолокационного разреза горизонтальные элементы и оставит высокоамплитудные сигналы, имеющие наклонные оси синфазности. Полученный разрез будет содержать наибольшие значения в местах просадки. Для более точного определения участка просадки необходимо рассчитать мощность георадиолокационных сигналов вдоль разреза.
Зоны разуплотнения грунтов на георадиолокационном разрезе отображаются высокоамплитудными, «хаотичными» сигналами, эти участки автоматически можно выделить по высоким значениям дисперсии амплитудных значений сигналов. Так как на подобных участках, как правило, появляются нерегулярные сигналы, которые сильно отличаются по своим радиофизическим характеристикам от сигналов, полученных при зондировании слоистых сред.
Границы талых и мерзлых горных пород, по которым осуществляется сброс поверхностных вод или миграция надмерзлотных вод сезонно-талого слоя, отличаются высокими значениями амплитуд. Для выявления таких границ необходимо отслеживать повышенные амплитудные значения при автоматизированном выявлении осей синфазности георадиолокационных сигналов с помощью расчета коэффициента корреляции.
Алгоритм обработки данных георадиолокационного мониторинга
В благоприятных условиях георадиолокационная съемка может быть оперативно проведена по длинным (в несколько километров) профилям, а при использовании транспортных средств длина профиля ограничена только характеристиками аппаратного и программного обеспечения георадара. Большой объем получаемых георадиолокационных данных требует рационального подхода к созданию алгоритма их обработки. Для выявления указанных выше аномалий георадиолокационного разреза предлагается проводить обработку по алгоритму, представленному на рис. 3.
В начале обработки производится ввод исходных данных:
1. Количество трасс в георадиолокационном разрезе (КТ).
2. Сигнал-образец.
3. Амплитуда сигнала образца (А).
4. Время задержки сигнала-образца (ВЗС).
5. Номер обрабатываемой трассы (N). Выполнение первых двух процедур позволяет выявить аномалии г) и д), для определения наличия аномалий остальных типов производится последовательный расчет коэффициента корреляции сигнала-образца со всеми трассами анализируемого георадиолокационного разреза.
Рис. 3. Блок-схема алгоритма обработки данных аномалий георадиолокационных разрезов при дистанционном мониторинге влажности горных пород
Далее анализируется количество максимальных совпадений (КМ), т.е. производится поиск сигналов максимально похожих на образец, при этом регистрируется время задержки для каждого максимума (ВЗМ) и амплитудное значение сигнала (Ат). Наличие остальных аномалий определяется на основе анализа параметра КМ, в соответствии с выявленными амплитудно-временными особенностями георадиолокационных сигналов.
Предлагаемый алгоритм обработки аномалий георадиолокационного разреза дополняет ранее разработанное программно-методическое обеспечение [5], как решение подзадачи «Картирование неоднородностей».
Заключение
В результате проведенного анализа особенностей амплитудно-временных характеристик сигналов в зонах аномалий георадиолокационных разрезов установлены признаки их выявления и интерпретации, что позволило в сочетании с существующим программно-методическим обеспечением разработать алгоритмическое обеспечение георадиолокационного мониторинга влажности горных пород криолитозоны. Разработанный алгоритм основан на использовании статистических характеристик (дисперсии и корреляции), расчет которых не требует больших вычислительных мощностей и может быть эффективно использован при обработке больших объемов данных.