CC BY
54
Георадиолокационный метод определения физических параметров конструкционных слоев автомобильных и железных дорог
З.Б. Хакиев, В.Л. Шаповалов, А.А. Кругликов, В.А. Явна
Для решения задачи стабилизации автомобильных и железных дорог необходимо располагать постоянной и объективной информацией о состоянии конструкционных слоев. Для решения этой задачи используется метод георадиолокации [1], позволяющий получать непрерывную информацию от оборудования, размещенного на движущихся специализированных измерительных средствах [2, 3].
Вместе с этим метод георадиолокации, как неразрушающий геофизический метод определения свойств грунтов, необходимо использовать в комплексе с другими геофизическими или лабораторными методами. Наиболее достоверную информацию о состоянии грунта удается получить в случаях, когда комплекс включает пенетрацию или бурение с отбором проб грунтов для дальнейшего лабораторного изучения.
Скоростной режим георадиолокационных обследований протяженных участков автомобильных и железных дорог позволяет получать большие объемы информации, требующие интерпретации высокого качества. Естественно, применять комплексы геофизических методов, в объемах, диктуемых протяженностью сети автомобильных и железных дорог, нерационально в связи с требуемыми материальными затратами. Это обстоятельство повышает актуальность исследований по расширению круга задач, решаемых методом георадиолокации вне рамок геофизических комплексов. Данная работа посвящена разработке метода оценки показателя преломления и удельной проводимости конструкционных слоев автомобильных и железных дорог при регулярных георадиолокационных обследованиях. Эти физические параметры необходимы для диагностики
состояния конструкционных слоев [4-6], в частности для определения их толщины и влажности.
Теоретические модели при обработке радарограмм
При диагностике протяженных участков дорог используются широкополосные георадары импульсного типа [7]. Импульс электромагнитного излучения частично рассеивается грунтовой средой и отражается от границ ее слоев. Это отражение может регистрироваться в виде зависимости напряжения на элементах приемного тракта георадара от времени. Обычно эту зависимость называют трассой [1].
Ограничимся рассмотрением плоской электромагнитной волны с вектором напряженности электрического поля, определенным в точке с координатой г в виде [8]:
Ю 7 ■ ■
—пкт го(—т-() го(—т)
Е = Е0в с е с = Е0е'рге с ,
где Е0 - амплитуда волны, р - коэффициент затухания волны, о - круговая
частота, п - показатель преломления, с - скорость света в вакууме. Учет временной протяженности импульса, а также зависимостей р и п от частоты может выполняться в более строгих моделях Фурье-преобразованием полученных результатов.
Для простоты рассмотрим случай нормального падения электромагнитного излучения на границы изотропных грунтовых слоев. Это приближение хорошо соответствует случаю обследования грунтовых слоев, когда антенный блок георадара оснащен рупором из проводящих материалов [9] и поднят над поверхностью грунта. В этом случае можно не учитывать поляризацию излучения при описании взаимодействия с границами раздела грунтовых слоев.
Система уравнений для определения показателя преломления первого грунтового слоя
Рассмотрим случай падения плоской электромагнитной волны на границу грунтового слоя из воздуха. Будем отсчитывать координату грунтового разреза и номера точек радарограммы от верхней границы первого грунтового слоя.
В случае, когда грунтовый слой обладает проводимостью, его показатель преломления может быть представлен в виде:
П = п (1 + щ), где 12 = -1 и индексом указан номер слоя.
Напряженность электрического поля сигнала (Е), отраженного от верхней поверхности грунтового слоя определяется соотношением [8]:
Е = Еп
Л
П 2(1 + щ2) + 1-2щ
П 2(1 + щ2) +1 + 2щ
Учтем, что коэффициент затухания равен [8] р = — пщ. Обозначим
с
пЩ =~ А =Хх = —
(1)
тогда
Е = Е„
1
(П - 1)2 +%!
2 '
(П + 1)2 + ^
Обозначим Е' напряженность электрического поля сигнала, отраженного верхней поверхностью грунта при изменении его влажности. Определим отношение величин сигналов, отраженных от нижней границы слоя грунта при изменении влажности грунта: Е'
А =
Е
Соответствующее уравнение преломления приобретает вид:
А =
(«;-1)2 +Х1 (п +1)2+Х12
'(«;+1)2+х1 \1 (п -1)2+Х1
для определения показателей
(2)
с
2
С другой стороны, изменение влажности грунта приводит к изменению скорости распространения электромагнитного сигнала. В приближении нормального падения электромагнитного излучения на поверхность грунтового слоя пройденный путь (И) до границы раздела со вторым грунтовым слоем определится соотношениями: тДш, с
К = —с (3)
п
тДш | с или К = -;- ,
п
где т - временной интервал между соседними точками радарограммы, Дш1 и
Дш1 - число точек трассы радарограммы между границами первого грунтового слоя при разной влажности. Приравнивая величины И, получим второе уравнение, связывающее диэлектрические проницаемости сред при разных значениях влажности:
(4)
п = п
(Дш ЧДш1 у
Система уравнений (1, 2, 4) позволяет определить показатели преломления грунтового слоя при изменении его влажности. При этом необходимо учесть то обстоятельство, что коэффициент затухания р в (1) далее вычисляется с использованием преобразования Гильберта от радарограмм, построенных в шкале глубин, которая в свою очередь зависит от показателя преломления грунтового слоя. Поэтому решение системы (1, 2, 4) должно осуществляться самосогласованным методом.
Расчет удельной проводимости
Значение коэффициента поглощения и соотношение, связывающее показатель преломления и диэлектрическую проницаемость в проводящих средах определим согласно [5] в СИ:
2 4 2
2 , Ц0С 0
ш 1
И 2 I V ш2
И = —^8 + ^)8 +
2 4 2
2 Ц О С 0
ш
где ш - центральная круговая частота излучения; с - скорость света в среде; 8- диэлектрическая постоянная среды; ц0 - магнитная постоянная;
а- удельная проводимость среды.
В используемом диапазоне частот величину а можно считать малой [9, 10], тогда:
8
2 4 2
1 -Ш 8
(
1 +
Ц с 4а2^
V
2ш 8
22
и система уравнений для определения удельной проводимости принимает вид:
2 рл/в
а
Ц ос
а
2ш 8
2 4
Ц о с
(И2 -8)
(5)
(6)
1
Система уравнений для определения диэлектрической проницаемости второго грунтового слоя
Электромагнитное излучение, распространяющееся в грунтовом слое, будет испытывать отражение от его нижней границы. Рассчитаем величину сигнала, отраженного нижней границей грунтового слоя. Исходим из того, что на грунт нормально падает электромагнитное излучение с модулем вектора напряженности электрического поля, равным Е0. Тогда, после частичного отражения от верхней границы слоя, в среде будет распространяться излучение с напряженностью электрического поля, равной:
Е = Ео
'\(пх +1)2 +Х12
При распространении в грунте, излучение будет затухать, частично отражаться на границе раздела со следующим грунтовым слоем (показатель преломления щ ). Отраженный сигнал испытает затухание в первой грунтовой среде и частично отразиться от верхней границы первого грунтового слоя. Таким отражением в приемную антенну попадет излучение с напряженностью электрического поля, равной:
Е = Е
(П1-1)2 +х1 '(П +1)2 +х12
2
2
-2рЛ (п2- П1)2 +х2 +х1 -2х2x1 }(п2 + П1)2 +х 22 +х12 + 2х 2 х1
где толщину первого грунтового слоя можно определить согласно (3), а показатели преломления и удельные проводимости считаются известными после решения систем уравнений (1, 2, 4) и (5, 6) соответственно.
Определим отношение величин сигналов, отраженных от нижней
Е'
границы слоя грунта при изменении влажности грунта А = —. Обозначим
Е
( -П2 ■ - 2
(
1-
(П1-1)2 +Х12 1(П1 +1)2 +Х12
I I Т
(П1-1)2 +Х1
' (п1 +1)2 + ХГ
2
,-2 РА
,-2 р'А
тогда:
9 2 2 I ' ' 9 ' 2 ' 2 ''
л „ I (п2 - п) +х2 +х1 -2х2х1 _ /(п2- п1) +х2 +х1 -2х2x1
А2\2 2 2 0 = 1. ' Ч? '2 .2 . .
(п2 + п1) +х 2 +х1 + 2х 2x1 \(П2 + щ)2 +х 2 +х1 + 2х 2x1
(7)
Вычисление показателей преломления второго грунтового слоя можно выполнить решением системы уравнений, включающей уравнения (1, 4, 7). Для определения удельной проводимости второго грунтового слоя можно использовать систему уравнений (5, 6).
Расчет электрофизических параметров третьего и последующих грунтовых слоев можно выполнять по аналогии с расчетом для второго слоя.
1
е
а =
2
Применение разработанной методики к обследованию земляного полотна железнодорожного пути
Для исследования влияния сезонных (погодных) факторов на конструкционные слои железнодорожного пути в декабре 2009 г., мае и сентябре 2010 г проведены георадиолокационные измерения на перегоне Туапсе - Сочи Северо-Кавказской железной дороги. Измерения выполнены в осевом сечении железнодорожного пути с использованием экранированного антенного блока с центральной частотой 400 МГц, размещенного на специализированном вагоне. Высота подвеса антенного блока над поверхностью балласта составляет 0.35 м, что обеспечивает значения угла падения луча, равным примерно 150. Развертка выбрана равной 48 нс. Трасса радарограммы содержала N=512 точек. Точность привязки радарограмм к железнодорожному полотну обеспечена системами спутниковой навигации. Фрагмент результатов обследования (участок между 1928 и 1929 км) приведен на рис. 1.
1928 кт 1929 1928 кт 1929
1 2
Рис. 1 - Радарограммы после совмещения характерных особенностей.
Время обследования: 1 - декабрь 2009 г., 2 - сентябрь 2010 г.
Для определения значений Ат - числа точек трассы между границами грунтовых слоев при разной влажности, необходимых для составления уравнения (4), радарограммы выделенного фрагмента подвергнуты предварительной обработке. Она включала удаление верхней части, не относящаяся к грунтовой среде. Полученные результаты использованы для определения усредненной на выделенном интервале трассы. Для подавления систематических помех, выполнено вычитание усредненной трассы из трасс радарограммы. Для позиционирования границ раздела грунтов к полученному результату применено преобразование Гильберта [11], позволившее определить огибающую к трассе (Р(0).
Для иллюстрации решения уравнений (1, 2, 4), (5, 6) и (1, 7, 4) из обрабатываемых радарограмм выделен фрагмент, относящийся к выбранному участку железнодорожного пути. Размеры участка ограничены так, чтобы границы раздела слоев располагались практически горизонтально. На рис. 2 приведены зависимости Р() = Р(т- т) как функции число точек измерения т после усреднения обработанных трасс. Соответствующие толщины конструкционных слоев сведены в таблицу 1.
Ат,
Ат
Рис. 2 - Радарограммы после предварительной обработки. Время обследования: 1 - декабрь 2009 г., 2 - сентябрь 2010 г.
Для определения коэффициентов А и р , необходимых для решения уравнений (2) и (5), при обработке радарограмм не выполнялось вычитание средней трассы. Значения А для верхней и нижней границ балластного слоя
приведены в таблице 1.
Таблица №1
Исходные параметры и результаты расчетов
№ Параметры балластного слоя Декабрь 2009 г. Сентябрь 2010г.
1 Значение преобразования Гильберта трасс на верхней границе, отн.ед. 7243 6016
2 Параметр А для формулы (1) на верхней границе - 1.20
3 Д^ - видимая толщина (в точках трассы при расстоянии между ними г=0.0967 нс) 115 98
4 Коэффициент затухания р сигнала в балластном слое, -1 м 2.15 1.69
5 Показатель преломления балластного слоя 2,45 2,09
6 Диэлектрическая проницаемость балластного слоя 5.81 4.23
7 Удельная проводимость балластного слоя, См/м 0.027 0.018
№ Параметр подбалластного основания Декабрь 2009 г. Сентябрь 2010г.
8 Значение преобразования Гильберта трасс на верхней границе, отн.ед. 5344 3277
9 Параметр А для формулы (1) на нижней границе балластного слоя - 1.63
10 Дда2 - видимая толщина подбалластного основания (в точках трассы при расстоянии между ними г =0.0967 нс) 110 118
11 Коэффициент затухания рсигнала в верхней части подбалластного основания, м- 2.63 1.93
12 Показатель преломления подбалластного основания 3,20 2,99
13 Диэлектрическая проницаемость подбалластного основания 9,89 8.79
14 Удельная проводимость подбалластного основания, См/м 0.044 0.030
При вычислении коэффициентов затухания необходимо учитывать то, что отраженный от границы грунтового слоя сигнал имеет пространственную протяженность и может существенным образом влиять на крутизну функции Р^) = Р(г • п / с) вблизи границ слоя. Расчет коэффициента затухания желательно осуществлять экспоненциальной аппроксимацией
Р(г) = Р(г0)е рг, начиная с точки г0, в которой значение Р(г) рассеянного
сигнала существенно больше аналогичного значения отраженного поверхностью сигнала. Очевидно, что на выбор значения г0 оказывает влияние соотношение показателей преломления слоев, образующих границу, пространственная протяженность импульса электромагнитного излучения, толщина грунтового слоя.
В данной работе величина г0 определена отступлением от границы слоя на 0.15 м. Исследования формы импульса позволяют установить, что на этом расстоянии от максимума (г=0) на частоте излучения 1700 МГц
обстоятельство обосновывает очевидное требование использовать коротковолновое излучение для исследования балластного слоя, толщина которого близка к 0.5 м.
отношение
равно 0.05, а на частоте 400 Мгц - 0.3. Это
1 -
о-
0
7!о ' 1^5 ' 2$ ' 2,5 г, т
0,0
0,1
0,2
г, т
Рис. 2 - Экспоненциальная аппроксимация в первом слое.
Рис. 3 - Экспоненциальная аппроксимация во втором слое. Время обследования:
1 - декабрь 2009 г.,
2 - сентябрь 2010 г.
Время обследования:
1 - декабрь 2009 г.,
2 - сентябрь 2010 г.
Результат аппроксимации зависимостей Р(г) приведены на рис. 2 и 3 для балластного слоя и верхней части подбалластного основания. Соответствующие значения р приведены в таблице 1.
Расчет показателей преломления балластного слоя и верхней части подбалластного основания
Для оценки влияния сезонных изменений на показатели преломления конструкционных слоев железнодорожного пути на выделенном участке железнодорожного пути решена система уравнений (1, 2, 4). Требуемые для решения параметры взяты из таблицы 1. При расчетах использована центральная частота импульса, определенная как максимальная частота Фурье разложения соответствующих радарограмм. Ее значение оказалась равным 330 МГц. Результаты параметрического решения системы уравнений для показателей преломления в декабре 2009 г. и сентябре 2010 г. приведены: для балластного слоя на рис. 4, для верхней части подбалластного основания - на рис. 5. Результаты расчета соответствующих показателей преломления представлены таблице 1.
Рис. 4 - Расчет показателей преломления балластного слоя:
1-декабрь 2009 г. (п'), 2 - сентябрь 2010 г (п).
Рис. 5 - Расчет показателей преломления подбалластного основания: 1-декабрь 2009 г. (п'), 2 - сентябрь 2010 г (п).
Видно, что изменение диэлектрической проницаемости в течение года может превышать 20%.
Расчет удельной проводимости и диэлектрической проницаемости балластного слоя и верхней части подбалластного основания
Для оценки влияния сезонных изменений на удельную проводимость и диэлектрическую проницаемость балластного слоя и верхней части верхней части земляного полотна железнодорожного пути решена система уравнений состоящей из (5) и (6). Необходимые для решения параметры приведены в таблице 1.
Результаты параметрического решения системы уравнений для диэлектрических проницаемостей и удельных проводимостей в декабре 2009 г. и сентябре 2010 г. приведены: для балластного слоя на рис. 4, 5 и для верхней части подбалластного основания - на рис. 6, 7. Соответствующие результаты расчетов приведены в таблице 1.
Рис. 6 - Расчет удельной проводимости и диэлектрической
Рис. 7 - Расчет удельной проводимости и диэлектрической
проницаемости балластного слоя в проницаемости балластного слоя в сентябре 2010 г. декабре 2009 г.
Рис. 8 - Расчет удельной проводимости и диэлектрической проницаемости подбалластного слоя в сентябре 2010 г.
Рис. 9 - Расчет удельной проводимости и диэлектрической проницаемости подбалластного слоя в декабре 2009 г.
Из приведенных данных можно сделать вывод о том, что на рассмотренном участке пути сезонное изменение удельной проводимости балластного слоя и подбалластного основания может составлять до 50%.
Выводы
Разработан метод определения сезонного изменения электрофизических свойств грунтовых слоев с использованием георадиолокации, который можно использовать без данных комплекса геофизических методов. Предложенный метод включает методику обработки георадиолокационных данных и теоретические методики определения:
- показателя преломления конструкционных слоев по величине отраженного сигнала и его «кажущейся» толщине;
- удельной проводимости грунтовых слоев по коэффициенту затухания и диэлектрической проницаемости;
Применение указанного метода к обследованию участка железнодорожного пути позволяет оценивать:
- диэлектрическую проницаемость балластного материала, необходимую для корректного профилирования нижней границы балластной призмы;
- степень загрязненности балластного материала;
- качество работы водоотводных сооружений (канавы, дренажи и пр.).
Работа выполнена при поддержке Министерства Образования и Науки РФ, Соглашение № 8403 от 24 августа 2012 г.
Литература:
1. Старовойтов, А.В. Интерпретация георадиолокационных данных. [Монография]: Учебное пособие. - М.: Изд-во МГУ, 2008. - 192 с.
2. Диагностический комплекс автоматизированной оценки состояния технических объектов инфраструктуры, связанных с обеспечением безопасности движения поездов «ИНТЕГРАЛ». [Электронный ресурс] -
Режим доступа: http://www.tvema.ru (доступ свободный) - Загл. с экрана. -Яз. рус.
3. Low-cost measurement of ballast quality, trackbed integrity and maintenance efficiency [Электронный ресурс] - Режим доступа: http://www.zeticarail.com (доступ свободный) - Загл. с экрана. - Яз. англ.
4. Явна, В.А. Определение эффективности георадиолокационного метода в грунтах с различными электрофизическими свойствами [Электронный ресурс] / В.А. Явна, З.Б. Хакиев, К.Ю. Кислица, М.С. Рассудов // «Инженерный вестник Дона», 2012, №2. - Режим доступа: http://www.ivdon.ru/magazine/archive/n2y2012/872 (доступ свободный) - Загл. с экрана. - Яз. рус.
5. Кругликов, А.А. Интеллектуальные системы мониторинга высоких железнодорожных насыпей [Электронный ресурс] / А.А. Кругликов, Г. И. Лазоренко, В.Л. Шаповалов, З. Б. Хакиев, В. А. Явна // «Инженерный вестник Дона», 2012, №3. - Режим доступа: http://www.ivdon.ru/magazine/archive/n3y2012/899 (доступ свободный) - Загл. с экрана. - Яз. рус.
6. Berthelot, С. "Ground-Penetrating Radar Evaluation of Moisture and Frost across Typical Saskatchewan Road Soils / C. Berthelot, D. Podborochynski, T. Saarenketo, B. Marjerison, C. Prang // Advances in Civil Engineering, vol. 2010, Article ID 416190, 9 pages, 2010. doi:10.1155/2010/416190
7. Шаповалов, В.Л. Георадиолокационная диагностика железнодорожного пути в скоростном режиме [Текст] / B.JI. Шаповалов, В.В. Ковдус, З.Б. Хакиев, А.А. Зарифьян, В.А. Явна, В.В. Помозов, А.В. Дудник, Б.А. Наумов // Труды третьей научно-технической конференции с международным участием «Современные проблемы. проектирования, строительства и эксплуатации железнодорожного пути», Москва, 2006. С. 58-60.
8. Борн, М. Основы оптики [Монография] / М. Борн, Э. Вольф. - 2-е изд., перераб. и доп. - М.: Наука, 1973, - 716 с.
9. Topp, G. C. Electromagnetic determination of soil water content: Measurements in coaxial transmission lines.
Water Resources Research // G. C. Topp, J. L. Davis, A. P. Annan, Volume 16, Issue 3, January 1980, pages 574-582.
10. Khakiev, Z. Efficiency evaluation of ground-penetrating radar by the results of measurement of dielectric properties of soils [Текст] / Z. Khakiev, K. Kislitsa, V. Yavna // Journal of Applied Physics, DOI: 10.1063/1.4770470, Vol. 112, Issue 12.
11. Гахов, Ф.Д. Краевые задачи [Монография] // Ф.Д. Гахов. - М. : ГИФМЛ. - 1958. - 544 с.
