CC BY
97
УДК 630.52:587/588
DOI 10.21685/2307-5538-2017-3-6
С. П. Санников, В. В. Побединский, А. А. Побединский
МОДЕЛЬ РАССЕИВАНИЯ РАДИОВОЛН В ЛЕСНОЙ СРЕДЕ
MODEL SCATTERING OF RADIO WAVES IN A FOREST
Аннотаци я. Актуальность и цели. Объектом исследования являлись участки лесной среды с различными лесонасаждениями. Предметом исследований являлись закономерности изменения величины рассеяния радиоволн в процессе радиочастотного мониторинга леса в зависимости от характеристик лесной среды, конструктивных параметров аппаратуры, параметров электромагнитного излучения и климатических условий. Цель работы заключалась в разработке модели рассеивания радиоволн в лесной среде. Материалы и методы. Для проведения исследований в качестве научных методов использованы методы математического моделирования, теории информации и передачи сигналов, методы расчета электромагнитных полей в лесных массивах, теория эксперимента, математической статистики. Для экспериментов были разработаны методики и оборудование для лабораторных и экспериментальных исследований, которое включает датчики (метки) RFID, передатчик, антенны трех типов, приемник, компьютер и специально разработанное программное обеспечение для обработки принимаемой информации. Результаты. Результатами являются предложенные расчетная модель распространения прямой волны в лесной среде, геометрическая расчетная модель распространения луча радиоволны в лесу и модель рассеивания радиоволн в лесной среде. Полученные результаты необходимы для проектирования систем радиочастотного мониторинга лесного фонда. Выводы. Полученная модель рассеивания радиоволн в лесной среде учитывает все основные влияющие параметры, дает достаточно точные результаты, адекватность модели подтверждена результатами экспериментальных исследований, поэтому результаты могут быть рекомендованы для использования при проектировании систем радиочастотного мониторинга леса.
Ab s t r a c t. Background. The subject of the study were sections of the forest environment with different forest plantations. The subject of the research were the regularities of the change in the magnitude of the scattering of radio waves in the process of radio frequency monitoring of the forest, depending on the characteristics of the forest environment, the design parameters of the equipment, the parameters of electromagnetic radiation and climatic conditions. The aim of the work was to develop a model for scattering radio waves in the forest environment. Materials and methods. To conduct research as scientific methods, mathematical modeling methods, information theory and signal transmission, methods for calculating electromagnetic fields in forest areas, theory of experiment, mathematical statistics were used. For the experiments, techniques and equipment for laboratory and experimental studies were developed that includes RFID sensors, transmitters, three types of antennas, a receiver, a computer and specially developed software for processing received information. Results. The results are the proposed design model for propagation of a direct wave in a forest environment, a geometric design model for the propagation of a radio wave beam in a forest, and a model for dispersing radio waves in a forest environment. The obtained results are necessary for the design of radio-frequency monitoring systems of the forest fund. Conclusions. The received model of radio wave scattering in the forest environment takes into account all the main influencing parame-
S. P. Sannikov, V. V. Pobedinsky, A. A. Pobedinskiy
ters, gives fairly accurate results, the adequacy of the model is confirmed by the results of experimental studies, therefore the results can be recommended for use in the design of radio frequency monitoring systems of the forest.
Ключевые слова: радиочастотный мониторинг лесного фонда; рассеивание радиоволн в лесной среде; комплексная диэлектрическая проницаемость участка леса; параметры лесной среды.
Key words: radio frequency monitoring the timber fund; the diffusing radio waves in timber ambience; complex dielectric constant of the forest area; parameters of the timber ambience.
В настоящее время для решения глобальной проблемы сохранения лесов приняты соответствующие программные документы на государственном уровне, в том числе «Основы государственной политики в области использования, охраны, защиты и воспроизводства лесов в Российской Федерации на период до 2030 года», которые предусматривают создание новых дистанционных систем наземного, авиационного и космического мониторинга пожарной опасности, защиты от незаконных рубок и использование инновационных информационных технологий. Разработанная [1, 2] с этой целью принципиально новая система (рис. 1) в виде сети устройств для непрерывного радиочастотного мониторинга лесного фонда в качестве одной из исходных величин использует значение рассеивания радиоволн в лесной среде. В наиболее подробных по этой теме исследованиях [3] указанный параметр получен на моделях деревьев с абсолютно гладкой поверхностью, геометрически правильной цилиндрической формы и равномерно распределенных по площади. Такая идеализация может давать значительные погрешности работы и приводить к функциональной непригодности всей системы радиочастотного мониторинга, которая должна учитывать перемещение лесоматериалов в лесу с точностью до одного бревна.
Рис. 1. Схема сети радиочастотного мониторинга лесного фонда: КРГО-1 - КРГО-4 - датчики; Р - мощность сигнала; Ж - влажность; Т - температура; п - количество деревьев; Ь - расстояние между датчиками; V - объемная доля /-го компонента лесной среды; а - константа вида лесного массива; ек - комплексная диэлектрическая проницаемость
Таким образом, определилась цель настоящих исследований, которая заключалась в разработке модели рассеивания радиоволн в лесной среде. Для достижения цели решались следующие задачи:
- обоснование перечня параметров, необходимых для создания адекватной модели;
- разработка геометрической расчетной модели распространения луча радиоволны в лесу;
- разработка математической модели в общем виде рассеивания радиоволн в лесной
Введение
<
среде;
- разработка методики и аппаратуры для экспериментальных исследований процессов радиочастотного мониторинга лесной среды;
- выполнение натурных экспериментальных исследований в лесной среде и оценка адекватности теоретической модели.
Для проведения исследований в качестве научных методов использованы методы математического моделирования, теории информации и передачи сигналов, методы расчета электромагнитных полей в лесных массивах, теория эксперимента, математической статистики.
Расчетная схема представлена на рис. 2, где на среднее поле прямой волны влияют диаметры стволов деревьев ё, диэлектрическая проницаемость древесины ед и остальных элементов лесной среды, хвои, листьев, воздушного пространства вокруг деревьев, поэтому в расчетах будем использовать усредненные значения параметров растительности:
- диаметр ствола деревьев ё;
- высота деревьев к;
- расстояние источника электромагнитных волн до точки измерения поля, т.е. приемной антенны ресивера I;
- длина волны X;
- амплитуда энергии источника электромагнитных волн и или вектор рассеяния поля бесконечным цилиндром и(х);
- электрическая проводимость древесины о;
- диэлектрическая проницаемость среды и древесины е;
- распространение электромагнитной волны происходит в лесной среде со средней постоянной густотой деревьев О.
Рис. 2. Расчетная модель распространения прямой волны в лесной среде: Н, Ь - координаты модели (высота, расстояние); Тх, Ях - трансмиттер, приемник; и(х) - энергия источника электромагнитных волн, которая ослабляется на пути Ь; - координата приемной антенны с диаметром ё; га - энергии электромагнитной волны в точке приема
Деревья на трассе распространения потока электромагнитной волны располагаются дискретно и случайным образом, а общую среднюю густоту деревьев О можно подсчитать по формуле
О 104N / (,)
О =--, шт./га, (1)
где N. - количество деревьев на трассе распространения потока электромагнитной волны, шт.; 5 - площадь вдоль трассы распространения потока электромагнитной волны с шириной 1 м, м2.
Приведенные параметры рассматриваются в работе [3] на моделях деревьев с абсолютно гладкой поверхностью цилиндрической формы и распределенных на модельной территории равномерно.
Для получения зависимости для реальных условий введем понятие комплексной диэлектрической проницаемости лесной среды е = Дед, ел, ех, ек, 8т, ез, ев). Здесь индексы при переменных е указывают на принадлежность к следующим элементам лесной среды:
- д - древесине ствола дерева;
- л - листьям;
- х - хвое;
- к - кроне деревьев и кустов;
- т - траве;
- з - земле (почве);
- в - воздуху.
Среднее поле рассеивания электромагнитных волн в лесной среде W (ld), где деревья представляют собой объекты со случайным расположением и случайной неоднородностью, зависит от множества параметров, влияющих на когерентность прямой волны I (ld). Рассеивание волн в такой случайной дискретной лесной среде подробно исследовано А. Исимару [4], но имеются ограничения, которые в нашем случае являются недостатками. Так, при рассмотрении автор использовал модели в виде плоских прямоугольных поверхностей, которые не сочетаются с формой стволов деревьев. Автор использовал обозначение функции ra (xa) энергии электромагнитной волны в точке приема а на радиусе распространения волны. Напряженность электромагнитного поля ультракоротких волн (УКВ) в точке приема зависит от протяженности l с встречающими элементами на пути с диаметром d. Отсюда следует, что функция ra (xa) в произвольной точке приема а связана с величиной рассеивания радиоволн W (ld)
на протяженности ld до приемной антенны в лесу с определенным количеством деревьев на пути, определяемым густотой деревьев G (1):
Fa (xa ) ^ W (ld ). (2)
Тогда для вывода уравнения расчета потенциала U (ra) [В, Вт, дБм] электромагнитной волны в точке а на пути распространения радиосигнала в лесном массиве применим теорию многократного рассеивания радиоволн [4]:
N
U(Fa) = £EQ¥s,), (3)
s=1
где Е - величина напряженности электромагнитного поля в направлении ra на пути распространения [В, Вт, дБм]; - параметр, влияющий на прохождение радиоволн в массиве лесной среды определенным индексом дерева s, для /-го измерения.
Параметр Y S имеет прямую зависимость с протяженностью канала передачи данных.
С учетом вышесказанного модель в общем виде величины рассеивания радиоволн W (ld) в лесной среде выглядит так:
W(l ) =
kl(< E, J,S,a,e*,G >
(4)
где к - волновое число среды, м 1; I - протяженность от источника радиоволны до точки наблюдения, м; 3 - плотность тока в антенне, А/м2; - площадь поперечного сечения дерева, м2; е - комплексная диэлектрическая проницаемость среды, Ф/м; О - густота деревьев в лесной среде, шт./м2; о - удельная проводимость среды (волновое сопротивление), См/м.
Известно, что радиоволны распространяются равномерно во всех направлениях от источника сигнала. Для определения величины рассеивания радиоволн Ж (1Л) на протяженности
и в точке а расположения приемной антенны в лесу рассмотрим один луч в направлении антенны ресивера га, который выразим в векторной форме с углом раствора а. Для исследования
этого вопроса разработана геометрическая расчетная модель распространения луча радиоволны в лесу, ее графическое представление показано на рис. 3.
Рис. 3. Геометрическая расчетная модель распространения луча радиоволны в лесу: U, U(x) - энергия источника электромагнитных радиоволн и вектор ее распространения; ra - радиус-вектор расстояния до антенны ресивера (приемника); d - размер антенны; ld - координата приемной антенны; L - протяженность распространения радиоволн; b - длина антенны; а - угол раствора луча
Размер d антенны ресивера (приемника) ra выбирается кратным длине волны в сторону увеличения и в сторону уменьшения, например: 2Х; X; ЛХ; и т.д. Тогда плотность энергии радиоволн в точке приемной антенны ресивера с координатой ld зависит от волнового числа среды k и площади антенны S:
S = wG = EH,
(5)
где V - плотность энергии электромагнитной волны, которая зависит от электрического и магнитного поля V = + wм; Е, Н - электрическое и магнитное поле электромагнитной волны; О - коэффициент пропорциональности, густота деревьев в лесной среде.
Вектор плотности потока электромагнитной энергии радиоволны принято называть вектором Умова - Пайнтинга [5], и с его учетом формула (5) записывается в следующем виде:
S = [ E, H ],
или в дифференциальной форме:
dW dt
+ divS = 0.
(6)
(7)
Из курса физики известно, что электрическое wэ и магнитное wм поля связаны между собой и распространяются в противоположных плоскостях, зависят от электрической е и магнитной ц проницаемости среды. Тогда запишем общую плотность энергии электромагнитной волны с учетом сказанного:
w = Wэ + wM = (So s E2)/2 + (ц Ц H2)/2,
(8)
где индекс 0 означает начальную проницаемость.
Частота электромагнитных волн в экспериментах по настоящей работе была принята 0,9 и 2,4 ГГц. Поскольку длина волны составляет несколько десятков миллиметров (точнее 330 и 135 мм соответственно), а расстояния до объекта наблюдения от источника сигнала больше 1 м, то магнитной проницаемостью в дальнейшем можно пренебречь в силу ее ничтожно малой величины.
Из уравнения (8) можно предположить, что энергия распространения электромагнитной энергии зависит от проводимости лесной среды. С другой стороны, лес является неоднородной анизотропной средой с хаотично расположенными элементами лесной среды, их размеры по-разному влияют на распространение радиоволн.
Для учета этих факторов введем понятие комплексной проницаемости и проводимости леса. Тогда для расчета эффективной проводимости леса оэ на участке от передатчика до приемника можно получить выражение
=Ё а, (9)
где ог - проводимость г элемента леса на пути распространения радиоволны; п - количество элементов леса на пути распространения радиоволны (деревья, подрост, кусты и пр.).
Под элементами лесной среды, которые влияют на распространение электромагнитной энергии от передатчика, будем понимать:
- стволы деревьев;
- ветки деревьев и кустарников;
- листья;
- траву;
- крону деревьев как интегрированный элемент дерева;
- воздушное пространство в лесу.
Все перечисленные элементы лесной среды имеют объем со свойствами, влияющими на распространение радиоволны. Для каждого из элементов лесной среды значения можно усреднить и измерить отдельно, например в лаборатории. По данным [5], высокочастотная энергия электромагнитных волн распространяется по поверхности материального тела. В данном случае все элементы леса являются отдельно расположенными в пространстве, имеют различные размеры, и можно сделать вывод, что они дискретны. Для простоты объяснения примем воздушное пространство в лесу как сплошную среду. Эффективная проводимость леса является комплексной величиной, поэтому по правилу амплитуд в математике запишем выражение для ее расчета:
аэ = а' + га", (10)
где о', о" - действительная и мнимая части комплексного числа; г - мнимая единица, равная V-!
Рассмотрим среду распространения электромагнитной волны (рис. 2) с условной электрической проводимостью о = 0. Пусть распространяется электромагнитная волна в направлении, задаваемом единичным вектором т. Тогда можно записать волновое уравнение в одномерном пространстве х = тг, которое является произведением единичного вектора т и радиус-вектора расстояния до точки наблюдения г для сферического распространения электромагнитных волн [5]:
Э 2 Ё = с2 Э 2 Ё
'ч -ч 2 '
Эг ец Эх2
где Е - вектор напряженности электромагнитного поля в лесной среде; с - скорость света; е, ц - диэлектрическая и магнитная проницаемости лесной среды; х - координата в направлении распространения электромагнитной энергии.
Уравнение (11) описывает две волны в направлении ± т с некоторой фазовой скоростью V = ю/к:
с2 с
— = = V . (12)
ец ^ец
Тогда покажем напряженность электромагнитного поля для скалярного (круглые скобки) и векторного (квадратные скобки) произведения членов волнового уравнения:
Э
ёыЁ = — (тЁ),
ЭЭх (13)
го(Ё = — [тЁ ]. Эх
Воспользуемся уравнением Максвелла с нулевой проводимостью в векторной форме [6]:
9 ,-ЕЛ дЙ
—(тЕ) = -|10|1—, 9х т
9 (тН) 9Е
—(тН) = —£о^——, 9х 9t
•л
—[тЕ ] = о,
9х
(14)
9х
[тн ] = о.
После совместного решения уравнений (14) и (11) получим связь между полями в электромагнитной волне:
Е = — М: [тН ]. V еое
(15)
Отношение магнитной проницаемости к диэлектрической проницаемости обозначим
символом X = и назовем волновым сопротивлением лесной среды [5]. Величина Z ука-
V еое
зывает на характеристику среды, которая состоит из двух частей. Соответственно отношение
/— - это волновое сопротивление свободного пространства, а отношение ./— - волновое \ео Vе
сопротивление среды, в которой распространяется радиоволна.
Перепишем уравнение (15) для гармонически меняющейся энергии поля в лесу в любой ее точке:
Е(х, 0 = Ео(х) ехр(-/Ю),
(16)
где х, t - координата и время; ю - круговая частота.
Из уравнения (16) видно, что энергия поля меняется от положения координаты, времени и частоты [7].
Теперь для уравнения (11) избавимся от производных по времени. С этой целью для решения уравнения Гельмгольца воспользуемся методикой, предложенной В. Б. Ивановым [5]. Условие решения уравнения с гармонической зависимостью от начальной пространственной координаты Ео ~ ехр(/кх) определяется выражением
<2Е с2 с12Е
< е— <Лх
или после соответствующих преобразований [5]:
2
<2Е ю2 (
<х2 с2
Л
е +1-
е ою у
—Ео = о .
(17)
(18)
Из уравнения (18) получим дисперсионное соотношение с учетом характеристики лесной среды, называемое в теории волн волновым числом среды [5]:
к2 =
ю
2
е +1-
еою у
ю2 * — = —г-—е . с
(19)
Здесь проводимость тождественна эффективной проводимости лесной среды (1о), или оэ = о, и является комплексной величиной, тогда и диэлектрическая проницаемость лесной среды становится комплексной благодаря проводимости:
9
G *
S + i-= S
(20)
которую по закону амплитуд можно записать в следующем виде:
* /,■ /Г
е = е + /е .
(21)
С учетом вышеизложенного анализа распространения электромагнитных волн в лесной среде, где на ее параметры влияют физические свойства, такие как проводимость, диэлектрическая проницаемость, рассмотренные в формулах (6)-(21), перепишем формулу через волновое уравнение (1):
G = ke^/t(E),
(22)
где t - тензор поля распространения электромагнитной энергии во времени.
Вернемся к выражению (4), которое перепишем с учетом уравнений (5), (19), (21) и модели распространения луча (рис. 3) и получим следующее соотношение для расчета величины рассеивания радиоволн в отдельно взятой точке лесной среды:
W ( x) = exp \_ikx (1 + GSa (S-1) )
(2З)
где к - волновое число среды (определено выше); х - координата расположения приемного устройства; О - коэффициент пропорциональности, густота деревьев в лесной среде; 8а - площадь лесной среды, занимаемая антенной Тх или Кх; е - диэлектрическая проницаемость среды.
Перепишем уравнение (23) для лесной среды с учетом модели распространения луча (рис. 3), заменив координату а расстоянием от источника сигнала до приемника, а также перейдя от площади антенны к ее диаметру и получим модель рассеивания радиоволн в лесной среде в общем виде:
W (id ) = exp
ikd
1 +
G nd2
(S* -1)
(24)
где I - протяженность от источника радиоволны до точки измерения; й - диаметр антенны (луча электромагнитных волн); е - комплексная диэлектрическая проницаемость лесной среды в точке приема сигнала.
Для подтверждения адекватности модели проведены экспериментальные исследования радиочастотного мониторинга лесной среды. Разработанный экспериментальный комплекс изображен на рис. 4. Фрагменты проведения экспериментальных исследований показаны на рис. 5. Результаты эксперимента показаны на рис. 6.
Рис. 4. Структурная схема установки сети беспроводных ИРГО-устройств: 1 - компьютер; 2 - анализатор спектра; 3 - ХВее модуль (приемник); 4 - ХЬее-модуль (передатчик); 5 - ИРГО-датчик; 6, 7 - таймер
а) б)
Рис. 5. Выполнение экспериментальных исследований: а - установка датчиков; б - установка приемника
Рис. 6. Результаты экспериментальных исследований: а - рассеивание радиоволн в сосновом лесу; б - рассеивание радиоволн в березовом лесу; 1 - экспериментальная зависимость; 2 - теоретическая зависимость
Заключение
Подытоживая сказанное, можно заключить, что предложенная модель учитывает все основные параметры лесной среды, конструкции радиочастотных устройств, климатические факторы и является достаточно адекватной, что подтверждено результатами, полученными экспериментально.
Библиографический список
1. Зависимость падения мощности сигнала от параметров лесной среды при радиочастотном мониторинге лесного фонда / С. П. Санников, В. В. Побединский, И. В. Бородулин, М. А. Черницын, Н. С. Кузьминов // Системы. Методы. Технологии. - 2016. - № 4. -С. 181-187.
2. Пат. 2492891 Российская Федерация, МПК А62С 37/00 С1. Система обнаружения лесного пожара / Лисиенко В. Г., Санников С. П. ; заявитель и патентообладатель Екатеринбург : Урал. гос. лесотехн. ун-т. - № 2012117204/12 ; заявл. 26.04.2012 ; опубл. 20.09.2013, Бюл. № 26. - 3 с.
3. Магазинникова, А. Л. Статистические методы расчета УКВ полей в лесных районах : автореф. дис. ... канд. физ.-мат. наук : 01.04.03 / Магазинникова А. Л. - Томск, 1998. -23 с.
4. Исимару, А. Распространение и рассеивание волн в случайно-неоднородных средах : в 2 т. / А. Исимару. - М. : Мир, 1981.
5. Иванов, В. Б. Теория волн : курс лекций / В. Б. Иванов. - Иркутск : Иркут. ун-т, 2ооб. -С. 76.
6. Лавров, В. М. Теория электромагнитного поля и основы распространения радиоволн / В. М. Лавров ; под. ред. И. Г. Кляцкина. - М. : Связь, 1964. - 368 с.
7. Яковлев, О. И. Распространение радиоволн : учебник / О. И. Яковлев, В. П. Якубов, В. П. Урядов, А. Г. Павельев ; под ред. О. И. Яковлева. - М : ЛеНАНД, 2оо9. - 496 с.
Санников Сергей Петрович
кандидат технических наук, доцент,
кафедра автоматизации
производственных процессов,
Уральский государственный
лесотехнический университет
(Россия, г. Екатеринбург, ул. Сибирский тракт, 37)
E-mail: ssp-54@mail.ru
Sannikov Sergey Petrovich
candidate of technical sciences, associate professor,
sub-department of automation
of production processes,
Ural State Forestry Technical University
(37 Sibirskiy trakt street, Ekaterinburg, Russia)
Побединский Владимир Викторович
доктор технических наук, профессор, кафедра сервиса и технической эксплуатации, Уральский государственный лесотехнический университет (Россия, г. Екатеринбург, ул. Сибирский тракт, 37) E-mail: pobed@e1.ru
Pobedinskiy Vladimir Viktorovich
doctor of technical sciences, professor, sub-department of service and technical operation, Ural State Forestry Technical University (37 Sibirskiy trakt street, Ekaterinburg, Russia)
Побединский Андрей Анатольевич
соискатель,
кафедра автоматизации
производственных процессов,
Уральский государственный
лесотехнический университета
(Россия, г. Екатеринбург, ул. Сибирский тракт, 37)
E-mail: pobed@e1.ru
Pobedinskiy Andrey Anatol'evich
applicant,
sub-department of automation production processes,
Ural State Forestry Technical University
(37 Sibirskiy trakt street, Ekaterinburg, Russia)
УДК 630.52:587/588 Санников, С. П.
Модель рассеивания радиоволн в лесной среде / С. П.Санников, В. В. Побединский, А. А. Побединский // Измерение. Мониторинг. Управление. Контроль. - 2017. - № 3 (21). - С. 39-48. БО! 10.21685/2307-5538-2017-3-6.
