Criteria for identification of the form and functional properties of the soil profile in georadiolocated studies by using georadar “Loza-B”

Бесплатный доступ

The valuable criteria for identification of the soil profile form and its some functional properties are considered as determined by interpretation of radarograms in the course of georadiolocated profiling of different territories using the geophysical device “Loza-B”. Apart from a standard set of procedures for digital processing of radarograms and constructing a temporary profile of changes in the reflected impulse amplitude along the direction of its sounding and transformation with depth, it is suggested to use the information on changing the dynamic characteristics of the wave field in homogenous media. This information involves the data about splash of the wave field, determination of diffraction points (synphasic axes) corresponding to horizons and objects under study, calculation of the frequency characteristics for distribution of reflected signals in the soil profile. A simple homogenous layer has a periodical sequence of picks in the amplitude spectrum within the depth of this layer. Changes in the thickness of the soil horizon are calculated by analyzing the local values of amplitude maximum of reflected impulses in the radarogram. The procedure of soil verification based upon the comparison of real soil profiles with standard ones within the framework of a definite classifier is conducted by statistic methods for calculating the polynomial functions of trend values of the signal amplitude according to its duration. The statistical processing includes calculation of average values, standard deviations and the distribution frequency of signal amplitude connected with the layer thickness and the amount of signal oscillations. The indices of Furie’s spectral transformation (frequency, period, spectral density) along with trend characteristics of polynomial modulations are considered as the major criteria for the two-staged identification of soils. In spite of attractive pattern of Furie’s spectral transformation the calculation of trend models should be testified by data obtained for soil profiles in field.

Еще

Amplitude of reflected impulse, georadar, identification, soil, radarogram, spectral density, trend

Короткий адрес: https://sciup.org/14313623

IDR: 14313623

Текст научной статьи Criteria for identification of the form and functional properties of the soil profile in georadiolocated studies by using georadar “Loza-B”

Принцип действия георадара основан на зондировании электромагнитными импульсами земной поверхности на различную глубину и восстановлении картины раздела сред с различной диэлектрической проницаемостью по отраженному сигналу (Владов и др., 2001; Янковский, 2005; Любушкин, 2007; Старовойтов, 2008). Этот подход можно рассматривать в качестве одного из направлений гносеологического познания природных объектов, в том числе и почвенного покрова, и формирования картины пространственной организации территории. Механизм такой организации, ее векторная составляющая важны для реализации создаваемых управленческих моделей природных объектов. Именно фактор универсальности, малой затратности, высокой скорости проведения полевых исследований может явиться приоритетом в использовании данного направления в реализации различных экономических, экологических, природоохранных проектов страны.

За последние десятилетия георадиолокация как неразрушающий метод успешно используется при изучении подповерхностных сред. Специалистами накоплен огромный положительный опыт как в проведении полевых работ, так и в приемах обработки и интерпретации георадиолокационных данных. Области применения георадарных обследований обширны – геология, строительство, археология, экология и др. (Методические..., 2003).

В почвенных исследованиях таких работ также становится все больше (Золотая и др., 2004; Калинкевич и др., 2007; Пягай и др., 2009; Fisher, 2013; Белобров и др., 2013; Воронин и др., 2013). Это связано со спецификой функционирования почвенного тела, как целостного образования с определенным набором генетиче- ских горизонтов, системой почвенных режимов и плодородием, но развивающегося в рамках однородной геологической среды (почвообразующей породы). Это обстоятельство вносит существенные коррективы в обработку радарограмм.

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

Практически во всех программах обработки георадарных данных предусмотрен стандартный набор приемов цифровой обработки сигналов: частотная фильтрация, вычитание среднего, сглаживание, выделение огибающей сигналов, настройка параметров визуализации и т.п., а в некоторых реализована возможность автоматизированного выделения и анализа поля обратного рассеяния и построения разреза по целому ряду атрибутов волнового поля и параметров исследуемой среды, вычисленных на основе данных атрибутов (Денисов и др., 2010).

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

Цель работы – разработать технологию интерпретации данных георадарного профилирования для оценки строения, свойств и состояния почв и почвенного покрова.

МЕТОДИКА ИССЛЕДОВАНИЙ

Анализ и обработка данных при георадиолокационных исследованиях, проводимых геофизическим прибором “Лоза-В”, осуществляется в специализированной программе Крот (Резников и др., 2000; Пягай и др., 2009; Пягай и др., 2013). Алгоритм обработки включает: установку параметров георадарных профилей (блок 2), цифровую фильтрацию данных, которая реализована в двух опциях: вертикальной и горизонтальной одномерной фильтрации (блок 3), вычисление средней скорости распространения волны в грунте (блок 4), ручную прорисовку границ слоев с разной диэлектрической проницаемостью. Интерпретация результатов георадарного зондирования завершается преобразованием временной шкалы радарограммы в глубинную (Методические..., 2003). Результаты зондирования в волновых формах при открытии в программе Крот могут представляться в цветных или в серых тонах картинках. Палитра из 128 цветов позволяет получить наглядную картину амплитудно-фазовых характеристик георадар-ного разреза. Такая предварительная обработка радарограмм характерна для абсолютно всех георадиолокационных исследований независимо от целей и объекта.

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

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

Моделирование поверхности полученных показателей осуществляется методом Inverse Distance, power-2 в программе Voxler 3.

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

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

Эти методы и критерии идентификации почв предполагают пространственно-геометрическую оценку расположения почвенных контуров, возможность применения ГИС-технологий при 3D-визуализации георадарных данных, использование математических и статистических методов оценки варьирования идентификационных георадарных показателей свойств почв и установление связей их с урожайностью сельскохозяйственных культур (Дмитриев и др., 1974; Любушкин, 2007).

РЕЗУЛЬТАТЫ И ОБСУЖДЕНИЕ

Критерии идентификации гранулометрического состава грунтов и их влажности. Глубина проникновения электромагнитных волн определяется электрическим сопротивлением и диэлектрической проницаемостью грунта. Для определения истинной глубины залегания границ раздела сред и объектов необходимо знать скорость распространения волны в соответствующей среде. Влажность и содержание минеральных солей ухудшают условия распространения электромагнитных волн и снижают максимальную глубину зондирования. Диэлектрическая проницаемость грунтов находится в тесной функциональной связи с влажностью и гранулометрическим составом пород (Денисов и др., 2010). Обобщение этих данных позволило составить номограмму изменения скоростей электромагнитного импульса для различных категорий влажности грунтов разного гранулометрического состава (рис. 1).

Диэлектрическая проницаемость (е)

Рис. 1. Изменение скорости импульса от диэлектрической проницаемости грунтов разного гранулометрического состава по категориям влажности (1 – максимальная гигроскопичность; 2 – наименьшая влагоем-кость).

Расчетные показатели диэлектрической проницаемости грунтов и скоростей распространения импульса в подземной среде являются критериями идентификации грунтов по гранулометрическому составу и влажности. Визуализация этих зависимостей на радарограммах сопровождается появлением характерных рисунков (рис. 2).

Электромагнитное профилирование и выделение локальных объектов искусственного происхождения. Георадар позволяет обнаруживать методами импульсной радиолокации объекты по значениям проницаемости и проводимости. Этими свойствами обладают почти все объекты искусственного происхождения в грунте: погребенные подземные сооружения (фундаменты, подвалы, ходы, захоронения), отдельные предметы, размеры которых превышают 5–10 см.

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

В качестве примера можно привести материалы георадарно-го профилирования территории косимой степи ЦентральноЧерноземного государственного природного биосферного заповедника им. В.В. Алехина, Курская область .

Пески           Суглинки             Глины

Рис. 2. Стандартный рисунок визуализации параметров отраженных импульсов, характерный для грунтов разного гранулометрического состава.

Рис. 3 . Поверхность локального объекта на территории косимой степи биосферного заповедника им. В.В. Алехина (Курская область), выделенного по осям синфазности георадарного профилирования.

На нескольких радарограммах на различной глубине четко выделяются оси синфазности (рис. 3). Соединение точек этих осей позволило выделить поверхность локального объекта размерностью 100 × 20 м, находящегося на глубине от 2 до 5 м.

Анализ отраженных импульсов волновых функций по трассам георадарного профилирования. Основная информация о среде заключена в амплитуде отраженного импульса, его полярности, времени задержки и ширине осцилляций. Реальная радарограмма в общем случае состоит из суперпозиции многих отраженных знакопеременных сигналов, которые могут накладываться друг на друга. Максимальные значения амплитуды отраженных сигналов, имеющие локальный характер распределения на временной оси радаро-граммы, являются критериями идентификации выделения границ почвенных горизонтов. Простой однородный слой имеет периодическую последовательность пиков в амплитудном спектре суммарного отражения от кровли до подошвы этого пласта. Границы горизонтов проведены по локальным амплитудным максимумам на временном отрезке радарограммы (табл. 1).

Однако при знакопеременной функции по этому алгоритму трудно отличить сигнал, отраженный новым слоем, от “хвоста” предыдущего сигнала. Необходим некоторый инвариантный по

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

Почва,    Радаро- угодье,    грамма

Импульс радарограммы

Профиль почвы

горизонты

3D-визуализация гумусового горизонта

Черноземы южные, пашня A–Bk– Cк–D

Черноземы типичные, целина A– AB–Bк– BCk–Ck

отношению к различным функциям импульса способ обнаружения отраженных сигналов (Резников и др., 2000; Копейкин, 2012). Такой способ предусматривает учет поведения электромагнитной волны при затухании в однородной среде, которое может протекать плавно или неравномерно. Кроме того, затухающий сигнал, в соответствии с законами электродинамики, при отражении от границ меняет свою полярность, если волна переходит из среды с меньшей диэлектрической проницаемостью или проводимостью, в среду с большими значениями этих параметров, и не меняет полярности в противоположном случае. Отраженный от размытой границы сигнал “растягивается” по времени. Слабая проводимость приводит только к “затуханию” сигналов, в то время как сильная – к их дисперсионному искажению (Копейкин, 2012).

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

Длительность затухания определяется мощностью передатчика и физическими свойствами среды. Энергетический потенциал георадара “Лоза-В” – 128 дБ, излучаемая мощность – 50 мВт. Глубина зондирования в зависимости от гранулометрического состава грунта может колебаться от единиц до десятков метров (для влажных глинистых грунтов 3–7 м).

Учет и оценка этих особенностей поведения электромагнитного импульса при почвенных исследованиях обусловили применение в качестве идентификационного метода статистическую обработку значений амплитуды сигнала длительностью 512 нс. В качестве критерия идентификации почв используется среднее значение амплитуды импульса. Предусмотрена возможность оценки вариаций этого параметра, его дисперсии. По нашим данным, вариация амплитуды отраженных импульсов по трассам профилирования может достигать 40–90%.

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

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

Таблица 2. Статистическая обработка амплитуды отраженных сигналов по трассам георадарного профилирования различных угодий почв ЕЧР

Почвы, угодье, местоположение

Число трасс, n*

Среднее значение амплитуды, X

Стандартное отклонение, s

Ошибка среднего, m

Коэффициент вариации, v , %

Каштановые, целина, Волгоградская область

5

128

118

5

93

Черноземы южные, пашня, Саратовская область

5

149

058

3

39

Черноземы типичные, некосимая степь, Курская область

6

112

059

3

53

Черноземы типичные, лесополоса, Курская область

6

120

065

3

54

Черноземы типичные, пашня, Курская область

9

132

060

3

46

Черноземы обыкновенные, косимый участок, Воронежская область

6

111

066

3

60

Черноземы обыкновенные, пашня, Воронежская область

6

128

054

2

42

Черноземы обыкновенные, лесополоса, Воронежская область

6

114

077

4

68

Серые лесные, пашня, Тульская область

5

119

062

3

52

Дерновоподзолистые,  паш

ня, Московская область

5

137

061

3

45

* В трассе 512 значений амплитуды.

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

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

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

Такая оценка становится возможной при применении спектрального анализа, который позволяет разложить любой электромагнитный сигнал, создаваемый источником, на более простые составляющие. Процесс разложения сложного апериодического и периодического сигналов на простые гармонические составляющие называется анализом Фурье. Среди атрибутов спектрального преобразования выделяется параметр спектральной плотности сигнала (spectral dencity), т.е. его распределение в частотной области, состоящей из многих близких частот, которые вносят наибольший вклад в периодическое поведение всего ряда. Чтобы убрать случайные колебания, периодограмму сглаживают, например, методом взвешенного скользящего среднего.

Ниже приведены примеры преобразования Фурье отраженных импульсов по некоторым объектам с распространением черноземов (рис. 4). Расчеты проведены в программе Statistica.

некос имая степь

1E5

ЙПППП

4 иии

0,00   0,05   0,10   0,15   0,20   0,25   0,30   0,35   0,40   0,45   0,50

лесополоса, Курская область

1E5

0,00   0,05   0,10   0,15   0,20   0,25   0,30   0,35   0,40   0,45   0,50

частота

1E5

частота

частота

пашня, Ку рская область

1E5

ЙПППП

0,00   0,05   0,10   0,15   0,20   0,25   0,30   0,35   0,40   0,45   0,50

частота

Рис. 4. Графики изменения спектральной плотности амплитуды отраженных сигналов по трассам георадарного профилирования черноземов различных сельскохозяйственных угодий.

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

Рис. 5. Изменение средних значений спектральной плотности отраженных импульсов по трассам георадарного профилирования черноземов (1 – обыкновенных, Воронежской области; 2 – типичных Курской области) по видам угодий.

Статистическая обработка выявила возможность использования показателя спектральной плотности в качестве критерия при идентификации почв (рис. 5).

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

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

Такие параметры применяются в качестве поправочных коэффициентов к реальным или расчетным значениям амплитуды для нивелирования эффекта дифракции импульса. Поправочные коэффициенты изменяются в пределах 0.1–1.0.

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

Метод полиномиальных функций предусматривает расчет тренда изменений амплитуды. Тренд, или тенденция f ( t ), представляет собой вероятностную устойчивую закономерность, наблюдаемую в течение определенного периода времени. При ге-орадарном профилировании общая длительность сигнала составляет 512 нс. Обычно тренд (тенденция) описывается с помощью той или иной функции f ( t ), как правило, монотонной аргументом которой является время. Тренды могут быть описаны различными уравнениями: линейными, логарифмическими, степенными и т.д. В нашем случае рассматриваются полиномиальные функции первого и второго порядков.

Фактический тип тренда устанавливают на основе подбора его функциональной модели статистическими методами либо сглаживанием исходного временного ряда. Сглаживание временного ряда, т.е. замена фактических данных расчетными значениями, имеющими меньшую колеблемость, чем исходные данные, является простым методом выявления тенденции развития. Соответствующее преобразование называется фильтрованием. Расчет проводится в программах: Крот, Excel, Statistica. При подстановке в модель (уравнение тренда) соответствующего значения фактора времени проводится точечный прогноз. Прогнозирование методом экстраполяции базируется на предположениях затухания амплитуды импульса, описывающегося плавной кривой, и учета случай-

Бюллетень Почвенного института им. В.В. Докучаева. 2015. Вып. 80. ности, которая позволяет оценить вероятность отклонения от закономерного развития.

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

Обычно в георадарном профиле отмечаются тренды одной полярности сигнала, но могут встречаться профили с чередованием разнознаковых трендов. Например, георадарный профиль балки в зоне распространения каштановых почв Волгоградской области. На южном склоне балки гумусовый горизонт смыт. На поверхность выходят карбонатные легкоглинистые отложения. Нижняя часть склона и днище балки заполнены делювиальными отложениями. Эта часть балки отмыта от солей. Не вскипает. Скопления гипса отмечаются на глубине более 3 м от поверхности склонов разной экспозиции. Исходные данные георадарного профилирования отфильтрованы, рассчитаны тренды изменения амплитуды по полиномам второй степени для 11-ти трасс. Пять из них, имеющие восходящий вид тренда, преобразованы в нисходящий. Введены поправочные коэффициенты на дифракцию сигнала рассматриваемой толщи через отношение средних значений спектральной плотности сигнала каждой трассы к максимальным их параметрам на данном профиле (табл. 3). Поправочные коэффициенты изменяются в пределах 1.0–0.4.

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

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

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

Время, nc

Расстояние между трассами, м

0

30

50

70

80

100

y

0.86 y*

y

0.40 y

y

0.66 y

y

1.00 y

y

0.82 y

y

0.88 y

001

129

111

164

66

158

104

092

092

150

123

102

90

040

115

099

149

60

132

087

118

118

120

098

107

94

080

100

086

132

53

115

076

131

131

095

078

105

92

120

086

074

118

47

100

066

136

136

075

062

096

84

160

072

062

102

41

095

063

130

130

060

049

084

74

200

060

052

090

36

095

063

112

112

050

041

065

57

240

045

039

078

31

100

066

085

085

045

037

040

35

256

040

034

072

29

104

069

072

072

042

034

030

26

* С поправочным коэффициентом на дифракцию отраженного импульса.

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

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

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

Для визуализации ее в программе Voxler по координатам строится модель с окраской значений амплитуд. Программа позволяет представить данные в плоскости и 3D-измерении.

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

ЗАКЛЮЧЕНИЕ

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

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

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

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

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

Статья научная