Методические указания МУК 4.3.1167-02 "Определение плотности потока энергии электромагнитного поля в местах размещения радиосредств, работающих в диапазоне частот 300 МГц - 300 ГГц" (утв. Главным государственным санитарным врачом РФ 7 октября 2002 г.) (документ не действует)

4.3. Методы контроля. Физические факторы

Методические указания МУК 4.3.1167-02
"Определение плотности потока энергии электромагнитного поля в местах размещения радиосредств, работающих в диапазоне частот 300 МГц - 300 ГГц"
(утв. Главным государственным санитарным врачом РФ 7 октября 2002 г.)

Дата введения: с момента утверждения

Область применения

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

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

Документ введен взамен МУК 4.3.680-97. Отличается от прежнего документа тем, что распространяется на более широкий класс антенн, содержит рекомендации по учету решетчатой структуры антенного рефлектора, влияния земли и крыши на величину плотности потока энергии в расчетной точке.

Не распространяется на радиосредства базовых станций систем подвижной связи и станций телевизионного вещания.

1. Радиосредства

На частотах выше 300 МГц работают различные радиосредства: радиорелейные системы передачи прямой видимости (РРСП ПВ), тропосферные радиорелейные системы передачи (ТРРСП), спутниковые системы передачи (ССП), радиолокационные станции (РЛС).

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

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

2. Расчет плотности потока энергии вблизи параболических антенн с круглой апертурой

Общие положения. Значение плотности потока энергии (ППЭ) апертурной антенны в произвольной точке пространства (в рамках энергетического подхода к решению задачи) определяется по формуле:

                   П     = П + П   + П   + П  , где                (2.1)

                    Сумма   а   обл   диф   пр

П     - апертурная составляющая;

 а

П     - составляющая, определяемая излучением облучателя;

 обл

П     - составляющая, обусловленная токами, протекающими  вблизи  кромки

 диф    зеркала (дифракционными токами);

П     - составляющая,  возникающая  в  результате   прохождения  энергии

 пр     сквозь  основное зеркало  антенны,  если  оно  имеет  решетчатую

        структуру.

Все пространство вблизи антенны условно делится на ряд характерных областей, которые ввиду симметрии показаны на рис. 2.1 только в секторе углов 0 <= Тэта <= пи.

Значение ППЭ в области I определяется апертурной составляющей П_а и составляющей облучателя П_обл.

В области II - заднем полупространстве антенны - ППЭ определяется составляющей П_диф. Если зеркало выполнено в виде решетчатой структуры, то к дифракционному полю добавится поле, прошедшее в область II сквозь ячейки решетки - составляющая П_пр.

В области III необходимо учитывать составляющие П_обл и П_диф. Область III в заднем полупространстве существует только для длиннофокусных антенн, когда Пси <= пи/2.

В области IV ППЭ определяется в основном составляющими П_обл и П_диф, но следует учитывать и П_а (особенно вблизи границы раздела областей I и IV).

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

Используемые допущения:

- амплитудное распределение поля по апертуре задается в виде "параболы на пьедестале":

                                              2

                                           2r

                f(r) = 0,316 + 0,684 [1 - (---) ], где             (2.2)

                                            d

     r - текущее значения координаты на диаметре апертуры,

     d - диаметр апертуры;

- облучатель и антенна имеют характеристики направленности с осевой симметрией относительно направлений их максимального излучения;

- характеристика направленности облучателя вне сектора углов перехвата энергии основным зеркалом считается неизменной и равной 0,316 по напряженности поля (по мощности 0,1);

- апертура имеет затенение, характеризуемое коэффициентом затенения d_T/d = 0,1 (d_T - диаметр "теневого диска", d - диаметр апертуры).

Расчет плотности потока энергии в области I. Плотность потока энергии в расчетной точке М (рис.2.2) представляется в виде двух составляющих:

                         П     = П + П   , где

                          Сумма   а   обл

П     - апертурная составляющая ППЭ (рис. 2.2а);

 а

П     - составляющая   ППЭ,  определяемая   непосредственно   излучением

 обл    облучателя (рис. 2.2, б).

В предположении осевой симметрии характеристик направленности апертуры и облучателя составляющие ППЭ имеют вид:

                         P       2    2

                  П = ------- D B (R)F (Тэта,R), Вт/м2             (2.3)

                   a        2  0

                      4 пи R

                          P         2

                 П   = ------- D   F   (гамма), Вт/м2, где         (2.4)

                  обл        2  обл обл

                       4 пи R

Р           - мощность, излучаемая антенной, Вт;

D           - КНД  антенны  в  направлении  максимального  излучения   в

 0            дальней зоне (величина безразмерная);

 2

B (R)       - функция, учитывающая зависимость КНД от расстояния;

 2

F (Тэта,R)  - нормированная  характеристика направленности  антенны  по

              мощности;

(Тэта,R)    - сферические координаты расчетной точки;

D           - КНД   облучателя  в  направлении  максимального излучения;

 обл

 2

F   (гамма) - нормированная характеристика направленности облучателя  по

 обл          мощности (угол гамма = 180° - Тэта).

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

Вводятся переменные: u - обобщенная угловая координата, х - относительное расстояние:

                      u = (пи d sin Тэта)/лямбда,                  (2.5)

                             x = R/R  , где                        (2.6)

                                    гр

d               - диаметр апертуры, м;

лямбда          - длина волны, м;

       2

R  = 2d /лямбда - граничное расстояние, начиная с которого можно считать,

 гр                что расчетная точка находится в дальней зоне.

С учетом введения обобщенных координат выражение (2.3) принимает вид:

                                 2    2

                         P лямбда    B (x)  2

                     П = --------- D ----- F (u,x), Вт/м2          (2.7)

                      а          4  0  2

                          16 пи d     х

Перевод размерности ППЭ Вт/м2 в мкВт/см2 осуществляется в (2.7) умножением на 100.

Переход от абсолютных значений величины ППЭ к относительным (децибелам относительно 1 мкВт/см2) осуществляется по формуле:

                         2

                 P лямбда                        B(x)

       П = 10 lg ---------  + 10 lg  D   + 20 lg[----] + 20 lg F(u,x)  + 3, дБ, (2.8)

        a            4                0           x

                    d

В(х)   - функция,   учитывающая   изменение   КНД   в   зависимости   от

         относительного расстояния;

F(u,x) - нормированная   характеристика    направленности   апертуры   в

         обобщенных координатах u, х.

С учетом принятых допущений выражение (2.4) при переходе к относительным значениям (децибелам относительно 1 мкВт/см2) приводится к виду:

                             P

              П   = 10 lg[-------] + 10 lg[D   ] + 10, дБ          (2.9)

               обл              2           обл

                          4 пи R

Аналитическое выражение функции В(х)/х для круглой апертуры с амплитудным распределением типа (2.2) имеет вид:

В области х < 0,105 функция (2.10) сильно осциллирующая, а в области x > 0,105 - изменяется монотонно. Осциллирующую часть функции следует заменить огибающей ее максимумов.

                                                                B(x)

     На  рис. П1.1  (приложение 1)  приведена  функция  20  lg[-----]. В

                                                                 х

                             B(x)

области х > 1 функция 20 lg [----] = - 20 lg x.

                              x

На рис. 2.3 показана круглая апертура с центральным затенением (затенение моделируется отсутствием элементов Гюйгенса в центре апертуры - светлый круг). Поле в точке М является суперпозицией полей элементов Гюйгенса dS, расположенных на поверхности раскрыва (апертуре).

Напряженность электрического поля, создаваемая в расчетной точке М всей совокупностью элементов Гюйгенса, находящихся в апертуре, определяется следующей формулой:

В (2.12) геометрические параметры Тэта_s, фи_s, r_s являются функциями Тэта и R.

Нормированная характеристика направленности апертуры имеет вид:

                                   |E(Тэта,R)|

                      F(Тэта,R)= ----------------                 (2.14)

                                  max|E(Тэта,R)|

В терминах координат u, x направленные свойства апертуры характеризуются функцией F(u,x). Процесс расчета каждой функции F(u,x) требует значительных вычислительных затрат, которые быстро растут с увеличением отношения d/лямбда. Функции F(u,х) сильно осциллирующие, поэтому в практических расчетах ППЭ следует использовать их гарантированные огибающие. Для удобства практических расчетов гарантированные огибающие табулированы (приложение 1, таблицы 1.1 и 1.2). При значениях х > 1, что соответствует дальней зоне, необходимо пользоваться огибающими для х = 1.

Значение КНД облучателя рассчитывается по формуле:

- характеристика направленности облучателя, реализующая амплитудное распределение вида (2.2) и обеспечивающая уровень 0,316 вне сектора углов перехвата энергии зеркалом (Пси_0 <= гамма <= пи).

Термин "облучатель" следует понимать, как собственно облучатель в однозеркальной антенне, так и систему "облучатель-вспомогательное зеркало" в двухзеркальной антенне.

График зависимости D_обл.дБ = 10 lg D_обл как функции аргумента Пси_0 для усредненной модели антенн приведен на рис. П1.2 (приложение 1).

Постановку задачи при расчете ППЭ в области I вблизи антенны с круглой апертурой иллюстрирует рис. 2.4, на котором центр апертуры - это О, а его высота над землей - Н_А. Угол альфа характеризует отклонение направления максимального излучения от плоскости горизонта. Ось Y системы координат XYZ совмещена с проекцией направления максимального излучения на плоскость XOY. Горизонтальная плоскость, на которой определяется ППЭ, находится на высоте Н_T над землей. Расчетная точка задается либо через координаты ро_м (расстояние до точки М вдоль поверхности земли) и фи_м, либо через y и x.

Заданными параметрами считаются: Р - мощность передатчика радиосредства, Вт; лямбда - длина волны, м; d - диаметр апертуры, м; D_0 - КНД антенны в направлении максимального излучения в дальней зоне; 2 Пси_0 - угол раскрыва антенны (угол перехвата энергии облучателя зеркалом).

В области I порядок расчета следующий:

1. Вычисляется расстояние от центра апертуры до расчетной точки М:

                                     2   2      2

                     R = кв.корень((Н - Н ) + pо )                (2.17)

                      М              А   Т      М

2. Рассчитывается угол между направлением максимального излучения и направлением линии "центр апертуры - расчетная точка М":

     Тэта = arсcos{[pо cos фи cos альфа - (Н - H ) sin альфа]/R } (2.18)

         М            М                     А   Т              М

3. Определяется граничное расстояние R_гр:

                                   2

                            R  = 2d /лямбда                       (2.19)

                             гр

4. Вычисляются параметры u, x по формулам:

                u = (пи d sin Тэта)/лямбда и х = R/R              (2.20)

                                                    гр

                                          В(х)

      5. Находится значение функции 20 lg[----], дБ (по огибающей).

                                           х

6. Вычисляется значение функции 20 lg F(u,x), дБ(по огибающей).

7. По формуле (2.8) рассчитывается значение составляющей П_а:

                        2

                Р лямбда                       В(х)

      П = 10 lg ---------] + 10 lg D   + 20 lg[----] + 20 lg F(u,х) + 3, дБ (2.21)

       a             4              0           x

                    d

8. Определяется:

                          D   , дБ = 10 lg D                      (2.22)

                           обл              обл

9. По формуле (2.9) находится значение составляющей П_обл:

                            P

             П   = 10 lg -------  + 10 lg  D     + 10, дБ         (2.23)

              обл              2            обл

                         4 пи R

10. Значение суммарной ППЭ рассчитывается по формуле:

                        П [дБ]     П   [дБ]

                         a          обл

                        ------     --------

                          10          10

              П     = 10       + 10        , мкВт/см2             (2.24)

               Сумма

Методика и порядок расчета имеют ограничение по минимальному удалению расчетной точки от центра апертуры - R_min >= d/2, что соответствует значению

                               ^      ^   лямбда

                          х >= х, где х = ------.

                                            4d

Расчет ППЭ для относительных расстояний х < х(^) выполнятся с помощью интерполяции. Для этого сначала по формуле (2.21) вычисляется П(х(^)) - величина ППЭ при х = х(^). Далее определяется П_s - усредненная величина ППЭ на апертуре по формуле:

                               2                       2

         П = 10 1g [400 P/(пи d 0,65)] ~ 10 lg [200 P/d ], дБ    (2.25)

          s

При этом интерполяционная формула имеет вид:

                                    ^

                              П - П(х)

                    1    ^     s        ^

                   -- [П(х) + -------- (х - х)]

                   10            ^

                                 x

         П (х) = 10                            , мкВт/см2         (2.26)

          a

Повторяя операции, предусмотренные порядком расчета ППЭ, последовательно для различных точек выбранного азимута фи = const определяется характер изменения ППЭ вдоль данного азимута, а также отыскивается точка, где значение ППЭ равно предельно допустимому уровню. Совокупность подобных точек на других азимутальных линиях, проведенных, например, через 1°, определяет границу соответствующей санитарной зоны в области I.

Расчет плотности потока энергии в области II. В области II поле обусловлено дифракцией электромагнитных волн на кромке параболического зеркала. Однако область II имеет три характерных подобласти (рис. 2.5), которые образуются после проведения двух касательных (1 и 2) в точках А и Б. Из подобласти II-а видна вся кромка антенны, из подобласти II-б видна часть кромки, из подобласти II-в кромка вообще не видна.

Исходными данными для расчета являются следующие параметры радиосредства: мощность излучения Р, Вт; длина волны лямбда, м; диаметр антенны d, м; половина угла раскрыва зеркала Пси_0; коэффициент использования поверхности зеркала k_un; альфа - уровень напряженности электрического поля на кромке зеркала. Постановка задачи иллюстрируется на рис. 2.6.

Сферические составляющие напряженности электрического поля в подобласти II-а имеют вид:

                       пи d                   пи d

              E    = интеграл dE    ; E   = интеграл dE           (2.27)

               Тэта     0       Тэта   фи      0       фи

В этих формулах: Е_0, определяется в соответствии с (2.36), альфа = 0,316, коэффициенты дифракции рассчитываются по формулам (2.28)...(2.32).

Переход к декартовым составляющим дифракционного поля осуществляется по формулам (2.39), а к сферическим компонентам ППЭ по формулам (2.38).

В секторе углов подобласти II-в следует принять, что П_диф = 0.

Расчет плотности потока энергии в области III. Если зеркало антенны длиннофокусное (Пси < 90°), то в области III ППЭ имеет две составляющие. Одна определяется излучением облучателя, другая - дифракцией на части кромки:

                         П      = П   + П                         (2.42)

                          Сумма    обл   диф

Составляющая П_обл определяется по формуле (2.23) с последующим переходом к мкВт/см2, а составляющая П_диф также, как для подобласти II-б.

В случае короткофокусной антенны (Пси < 90°) в области III ППЭ имеет составляющие П_диф и П_а. Составляющая П_диф рассчитывается также, как для подобласти II-а - по формулам (2.28) ... (2.38).

Расчет плотности потока энергии в области IV. В области IV ППЭ определяется в основном составляющими П_обл и П_диф, поэтому значение ППЭ в этой области формально определяется по формуле (2.42) с добавлением составляющей П_а. При этом дифракционная составляющая поля определяется так же, как это сделано в подобласти II-а - по формулам (2.28) ... (2.38).

Расчет плотности потока энергии в области V. В этой области ППЭ следует определять следующим образом:

                         П      = П   + П

                          Сумма    обл   диф

При этом составляющая П_а рассчитывается по формулам (2.25) и (2.26), после того как предварительно геометрически определена граница раздела между областями IV и V (значение х(^)) на выбранном азимутальном направлении.

Алгоритм определения областей и подобластей для расчетной точки и примеры расчета ППЭ приведены в приложении 2 - рис. П1.3 и П1.4.

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

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

Квадратная апертура. При анализе квадратной апертуры (рис. 3.1) используются допущения:

- распределение амплитуды поля в одной из главных плоскостей - "косинус на пьедестале":

                                             пи тау

                   f(тау) = 0316 + 0,684 cos ------, где          (3.1)

                                               а

а - сторона квадрата, тау - текущее значение координаты апертуры в одной из главных плоскостей;

- облучатель и антенна имеют характеристики направленности с осевой симметрией относительно направлений их максимального излучения;

- значение характеристики направленности облучателя вне сектора углов перехвата энергии основным зеркалом равно 0,316 (по напряженности поля).

Общая расчетная формула для определения значения ППЭ имеет вид (2.1). Физический смысл отдельных составляющих прежний. Составляющие П_а и П_обл в децибелах относительно 1 мкВт/см2 будут иметь вид:

                        2

                P лямбда                      B(x)

      П = 10 lg --------  + 10 lg D   + 20 lg ----  + 20 lg F(u,x) + 3, дБ,  (3.2)

       a            4              0           x

                   a

                         P

          П   = 10 lg -------  + 10 lg  D     + 10, дБ, где      (3.3)

           обл              2            обл

                      4 пи R

Р                          - мощность, излучаемая антенной, Вт;

лямбда                     - длина волны, м;

а                          - сторона квадрата (апертуры антенны), м;

D                          - КНД  антенны   в  направлении максимального

 0                           излучения в волновой зоне;

В(х)                       - функция,    учитывающая   изменение     КНД

                             квадратной   апертуры   в  зависимости   от

                             относительного расстояния;

F(u,x)                     - нормированная  ХН  квадратной  апертуры   в

                             координатах u, х;

Тэта, R                    - сферические координаты расчетной точки;

u = (пи а sin Тэта)/лямбда - обобщенная координата угла;

х = R/R                    - относительное расстояние;

       гр

       2

R  = 2а /лямбда            - граничное расстояние.

 гр

D                          - КНД  облучателя в направлении максимального

 обл                         излучения  (величина  безразмерная). График

                             зависимости    D_обл,дБ = 10 lg D_обл   как

                             функции  аргумента Пси = 0 для  усредненной

                             модели   антенн   приведен   на   рис. П1.2

                             (приложение 1).

Вычисление функции F(u,x) сводится к расчету характеристики направленности линейного синфазного источника (рис. 3.2) с распределением амплитуды тока, совпадающем с распределением амплитуды поля вида (3.1).

Значение напряженности поля в расчетной точке определяется выражением:

                            a

                           ---

                            2

                     Е = интеграл dE(Тэта,R), где                 (3.4)

                            a

                           ---

                            2

                                                 -jkr

                                f(тау)          e

          dE (Тэта,R) = i 60 пи ------ cos тэта ------ d тау       (3.5)

                                лямбда             r

В (3.5) f(тау) определяется распределением поля, а геометрические параметры тета, r являются функциями Тэта, R.

Нормированная характеристика направленности апертуры имеет вид:

                                   |E(Тэта,R)|

                     F(Тэта,R) = ---------------                   (3.6)

                                  max|E(Тэта,R)|

В терминах обобщенных координат (u,х) направленные свойства апертуры будут характеризоваться функцией F(u,x).

Функции F(u,x) сильно осциллирующие, поэтому в практических расчетах ППЭ следует использовать их гарантированные огибающие. Для удобства практических расчетов гарантированные огибающие табулированы (приложение 3, таблицы 3.1 и 3.2). При значениях х > 1, что соответствует дальней зоне, необходимо пользоваться огибающими для х = 1.

Аналитическое выражение функции В(х)/х для квадратной апертуры с амплитудным распределением типа "косинуса на пьедестале" имеет вид:

В области х <= 1 функция (3.7) осциллирующая, а в области х > 0,15 - изменяется монотонно. При расчете ППЭ осциллирующую часть функции следует заменить огибающей ее максимумов.

На рис. П3.1 (приложение 1) приведена функция 20 lg(B(x)/x). В области x > 1 функция 20 lg(B(x)/x) = -20 lgx.

Порядок расчета ППЭ в переднем полупространстве не отличается от приведенного для круглой апертуры.

Для расчета ППЭ в области заднего полупространства вводится понятие эквивалентной круглой апертуры. Диаметр эквивалентной апертуры определяется из условия равенства площадей квадратной и круглой апертур:

                         d = 2a/кв.корень(пи)                     (3.11)

                          э

Угол раскрыва эквивалентной апертуры определяется по формуле:

                                      d

                                       э

                       Пси = 2 arctg [---], где                   (3.12)

                          0           4 f

f - фокусное расстояние зеркала.

Прямоугольная апертура. Прямоугольная апертура со сторонами а и b показана на рис. 3.3.

Расчетная формула для апертурной составляющей ППЭ имеет вид:

Функции F(u,x), входящие в (3.13), вычисляются в соответствии с (3.4)...(3.6) с учетом перехода к обобщенным координатам u и х, а функции В(х)/х - по формуле (3.7) при х = x_1 и х = х_2, соответственно.

Расчет составляющей ППЭ от облучателя выполняется по формуле (3.3). КНД облучателя следует рассчитывать по формуле:

                      D   = кв.корень(D D ), где                  (3.15)

                       обл             1 2

D_1 и D_2 определяются по формулам (2.15) и (2.16) с учетом различных значений угла Пси_0 в главных плоскостях.

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

                                           ab

                          d = 2 кв.корень(----)                    (3.16)

                           э               пи

Угол раскрыва эквивалентной круглой апертуры определяется по формуле (3.12).

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

4. Расчет плотности потока энергии вблизи антенн типа параболический цилиндр и рупорных антенн

Антенна типа параболический цилиндр. Антенна имеет прямоугольную апертуру (рис. 4.1). Распределение амплитуды поля вдоль сторон апертуры в главных плоскостях XOZ и ZOY равномерное. Методика расчета ППЭ соответствует случаю прямоугольной апертуры при f(тау) = const. Огибающие F(u,x) для случая f(тау) = const отличаются от аналогичных кривых, соответствующих распределению (3.1) не существенно (единицы дБ). Поэтому в практических расчетах можно использовать данные, приведенные в табл. П3.1 и табл. П3.2.

Значение КНД облучателя рассчитывается по формуле:

                                     kL

               D    = ----------------------------------, где      (4.1)

                обл               cos(kL) - 2   sin kL

                      2 [Si(kL) + ----------- + ------]

                                     kL              2

                                                 (kL)

L - длина облучателя;

                                   2 пи

                              k = ------;

                                  лямбда

                        kL    sin x

            Si(kL) = интеграл ----- dx - интегральный синус.

                        0       x

В области заднего полупространства расчет ППЭ следует вести по формуле (3.3), приняв D_обл = D_0

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

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

                        Р       2             Вт

                 П = ------- D F (Тэта, фи), ----, где            (4.2)

                           2  p                2

                     4 пи R                   м

Р           - мощность, изучаемая антенной, Вт;

F(Тэта, фи) - характеристика направленности рупора;

D           - КНД рупора.

 р

Для пирамидального рупора в расчетах следует принять:

             F (Тэта, фи) = кв.корень(F (Тэта) F (Тэта)), где      (4.3)

                                       E        H

сомножители - это нормированные характеристики направленности рупорной антенны в Е и Н плоскостях.

Нормированная характеристика направленности имеет вид:

                                    f(Тэта)

                        F (Тэта) = -----------                      (4.6)

                         Е         max f(Тэта)

Для расчета характеристики направляемости в плоскости Н (фи = пи/2) сначала вычисляется функция:

Нормированная характеристика направленности имеет вид:

                                    f(Тэта)

                        F (Тэта) = -----------                      (4.9)

                         Н         max f(Тэта)

Значение КНД пирамидальной рупорной антенны рассчитывается по формуле:

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

Нормированная характеристика направленности имеет вид:

                                    f(Тэта)

                        F (Тэта) = -----------                     (4.19)

                                   max f(Тэта)

Значение КНД конической рупорной антенны рассчитывается по формуле:

                                           2

                                      r

                            D = 20 (------)                       (4.20)

                             p      лямбда

В области заднего полупространства расчет ППЭ выполнятся по формуле:

                             0,025 Р      Вт

                         П = ------- D , ----                     (4.21)

                                  2   р    2

                              пи R        м

Примеры расчетов ППЭ вблизи антенн типа параболический цилиндр, пирамидального и конического рупора приведены в приложении 4.

5. Расчет плотности потока энергии вблизи рупорно-параболической и перископической антенн

Конструкция типовой рупорно-параболической антенны (РПА) схематично представлена на рис. 5.1. Апертуру можно считать квадратной с размером 2,7 х 2,7 м2. Методика расчета ППЭ в переднем полупространстве (Y > 0) соответствует приведенной в разделе 3 для квадратной апертуры с распределением амплитуды поля (3.1) - "косинус на пьедестале". Составляющая П_а рассчитывается по формуле (3.2), составляющая П_обл - по формуле (3.3). При этом угол раскрыва рупора 2 Пси_0 = 35°. Учитывая, что РПА обладает хорошим защитным действием (уровень бокового излучения исключительно низок: почти во всем секторе углов он не превосходит - 60 ... - 70 дБ), дифракционную сооставляющую# и составляющую П_обл при расчете ППЭ в области заднего полупространства (Y < 0) не учитывать.

Перископические антенные системы (ПАС) обычно строятся по трехэлементной схеме (рис. 5.2) - первичный рупорный облучатель, нижнее зеркало и верхнее зеркало. Диаметр верхнего зеркала 3,9 м, нижнего 3,2 м. Диаметры апертур с учетом наклонного положения зеркал следует брать равными 0,7 реальных диаметров зеркал.

В общем случае значение ППЭ в произвольной точке пространства М определяется тремя составляющими, каждая из которых соответствует одной из трех антенн - А1, А2, A3 (рис. 5.3).

Антенна А1 рупорная. Ее вклад в значение ППЭ определяется в соответствии с методикой, изложенной в разделе 4. Антенны А2 и A3 имеют круглые апертуры. Их вклад в значение ППЭ определяется апертурными составляющими, которые рассчитываются по формуле (2.8).

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

В области Y < 0 следует ограничиться только составляющей от антенны А1, то есть не учитывать составляющие ППЭ, обусловленные дифракцией электромагнитных волн на зеркалах.

Примеры расчета ППЭ вблизи антенн РПА и ПАС приведены в приложении 5.

6. Учет влияния решетчатой структуры рефлектора

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

Значение составляющей П_пр в точке М, находящейся в области заднего полупространства, определяется по формуле:

                     2   Р         2

               П  = T ------- D   F   (гамма), Вт/м2, где          (6.1)

                пр          2  обл обл

                      4 пи R

Т           - коэффициент прохождения волны через  решетчатую  структуру

              (по полю);

Р           - мощность излучения облучателя, Вт;

D           - коэффициент направленного  действия  облучателя  (величина

 обл          безразмерная);

F   (гамма) - характеристика    направленности   облучателя    (величина

 обл          безразмерная);

R           - расстояние  от фазового  центра  облучателя  до  расчетной

              точки М, м.

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

Характеристика направленности облучателя в области углов 0 < гамма <= Пси_0 имеет вид:

Значение КНД облучателя рассчитывается по формуле (2.15).

Для случая выполнения рефлектора (зеркала) из сетки линейных проводников (рис. 6.2) формула для расчета коэффициента прохождения по полю имеет вид:

Функции Бесселя рассчитываются по формулам:

                    (2)              2       лямбда

                   H  (kpо) ~ 1 + i -- ln -----------              (6.4)

                    0               пи    1,781 пи ро

Ряд, стоящий в правой части (6.5), сходится достаточно быстро (можно ограничиться десятью членами).

В случае выполнения отражательного зеркала в виде поверхности со щелями (рис. 6.3а) при длине щелей, отвечающих условию (t > лямбда), считать, что полоски и щели имеют безграничную длину (рис. 6.3б).

Для практически важного случая d < (0,4...0,5) лямбда, формула расчета коэффициента прохождения имеет вид:

Для поверхности с круглыми отверстиями (рис. 6.4) и расстоянии между центрами отверстий d < (0,3...0,4) лямбда формула для расчета коэффициента прохождения имеет вид:

                                               t

                                3       - 1,6 ---

                          2 пи D               D

                    T = ----------- x 10         , где             (6.7)

                                  2

                        3 лямбда d

     d - расстояние между щелями,

     D - диаметр щели,

     t - толщина отражающей поверхности.

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

Пример расчета ППЭ с учетом решетчатой структуры рефлектора приведен в приложении 6.

7. Расчет плотности потока энергии вблизи вибраторных антенн

Расчет ППЭ вблизи вибраторных антенн осуществляется по формуле:

                             ->   ->

                   П = 50|Re{E  x H*}|, мкВт/см2, где              (7.1)

    ->   ->

    Е  и Н* - векторы   напряженности электрического (В/м) и магнитного,

              (А/м) полей, соответственно (символом "х" здесь  и  далее

              обозначается векторное произведение).

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

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

Интегральное уравнение Харрингтона имеет вид:

                                    ->     ->                           1        dG(l,l') dI(l')

      E (l) = i интеграл [омега мю (h (l), h (l')) G(l,l') I(l') + ------------- -------- ------] dl', где (7.2)

       0           L              0                                омега эпсилон    dl      dl'

                                                                                0

E (l)          - стороннее поле;

 0

L              - контур,   последовательно  проходящий   по  осям   всех

                 проводников - вдоль   данного   контура   отсчитываются

                 криволинейные координаты l и l' (см. рис. 7.1);

омега          - круговая частота;

мю             - магнитная постоянная;

  0

эпсилон        - диэлектрическая постоянная;

       0

->      ->

h (l) и h (l') - единичные   векторы,  тангенциальные  к  проводникам  в

                 точках l и l';

                        -i бета R(l,l')

                       e

             G(l,l') = ---------------- - функция Грина;

                         4 пи R(l,l')

            ->  -> 2    2    1/2

R(l,l') = {|r - r'|  + a (l)}    - расстояние      между        точками,

                                   соответствующими значениям l и l';

->  ->

r и r'                           - радиус-векторы    точек   l   и   l',

                                   соответственно;

а(l)                             - радиус проводника в точке l;

I(l')                            - искомая токовая функция.

Искомая токовая функция I(l') при выбранном базисе разложения определяется как линейная комбинация:

                                M

                      I(l') = Сумма {I B (l')}, где                (7.3)

                               m=1    m m

М      - число базисных функций, равное числу сегментов;

I      - коэффициент при m-ой базисной функции;

 m

В (l') - m-ая кусочно-синусоидальная базисная функция.

 m

Представление токовой функции в виде (7.3) сводит интегральное уравнение (7.2) к системе М линейных алгебраических уравнений относительно I_m:

                   М            N

                 Сумма Z  I = Сумма Z     I     , где              (7.4)

                  m=1   km m   i=1   k,M+i (вх)i

     k = 1, 2, ... М,

N     - число входов (активных вибраторов) антенны;

I     - заданные входные токи.

 (вх)i

Коэффициенты системы (7.4) рассчитываются по формуле:

Первые М сегментов вводятся собственно на проводниках антенны. Сегменты с номерами, превышающими М, вводятся в областях зазоров активных вибраторов, к которым подводятся питающие напряжения (токи). Соответствующие им коэффициенты Z_k,M+i, обусловливающие неоднородность системы (7.4), также рассчитываются по формуле (7.5).

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

Векторы напряженности электрического и магнитного полей в заданной точке наблюдения с радиус-вектором r(->) рассчитываются по формулам:

                      ->  M + N -> ->  M + N ->

                      Е = Сумма Е, H = Сумма H , где               (7.6)

                          m = 1  m     m = 1  m

     ->   ->

     E  и H  - векторы   напряженности   соответственно электрического  и

      m    m   магнитного   полей,   создаваемые в  точке наблюдения  k-m

               сегментом     (включая    сегменты    в   зазорах активных

               вибраторов).

Векторы Е(->)_m и Н(->)_m рассчитываются по формулам:

При m > M в фомулах (7.7) и (7.8) в качестве коэффициентов I_m берутся соответствующие заданные входные токи:

                        I = I        , m > N.

                         m   (вх) m-N

Пример расчета ППЭ вблизи вибраторной антенны приведен в приложении 7.

8. Учет влияния подстилающей поверхности (земли, крыши)

Учет влияния плоской безграничной поверхности сводится к применению двухлучевой модели (рис. 8.1).

В точке В плотность потока энергии рассчитывается по формуле:

           П = П (Тэта , фи , r ) + П (Тэта , фи , r ), где         (8.1)

                1     1    1   1     2     2    2   2

П (Тэта , фи , r ) - ППЭ, создаваемая прямой волной (луч АВ);

 1     1    1   1

П (Тэта , фи , r ) - ППЭ,  создаваемая  волной,   распространяющейся  по

 2     2    2   2    траектории АСВ;

                                    h + h          h - h

                                     1   2          1   2

                Тэта = Тэта + arctg ------ - arctg ------;          (8.2)

                    2      1          ро             ро

                                      2          2

                      r = кв.корень(ро + (h - h ) )                 (8.3)

                       1                   1   2

                                      2          2

                      r = кв.корень(ро + (h + h ) )                 (8.4)

                       2                   1   2

Обе составляющие рассчитываются по методике, соответствующей свободному пространству.

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

I - область прямых лучей, II - область интерференции прямых и отраженных лучей, III - область тени, в которой предполагается отсутствие поля.

Исходными данными для учета конечных размеров крыши вдоль азимута фи (рис. 8.3) являются: координаты фазового центра Ф{0, 0, z_Ф}, координаты расчетной точки М{х_М, у_М, z_M}, координаты углов крыши А{х_А, у_А, z_A}, В{х_В, у_В, z_B}.

Последовательность расчета:

1. Вычисляются значения некоторых вспомогательных параметров:

                    y - y        y (x - x ) - x (y - y )

                     B   A        A  B   A     A  B   A

               k  = ------, a  = -----------------------           (8.5)

                AB  x - x    AB          x - x

                     B   A                B   A

2. Определяются координаты точки К (точки пересечения кромки крыши АВ линии O кси):

                         пи    пи

      - для фи не равно ---, 3 ---

                         2      2

                           a                a  tg фи

                            AB               AB

                   x  = -----------,   у  = -----------;            (8.6)

                    K   tg фи - k       К    tg фи - k

                                 AB                 AB

                 пи     пи

      - для фи = ---, 3 ---

                  2      2

                              х = 0, у = а                         (8.7)

                               К      К   АВ

3. Рассчитываются расстояния:

                           2   2                        2   2

       l = OK = кв.корень(x + у ), l = ON = кв.корень (x + y )     (8.8)

        K                  K   K    M                   M   M

4. В соответствии с алгоритмом, приведенном на рис. 8.4, определяется принадлежность точки М к одной из трех областей.

5. Приведенный порядок расчета повторяется для всех выбранных значений фи, соответствующих положению линии O кси внутри сектора углов АОВ.

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

9. Инструментальный контроль уровней электромагнитных полей радиосредств

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

Измерения проводятся:

- на этапе предупредительного санитарного надзора (при приемке объекта в эксплуатацию);

- на этапах текущего санитарного надзора (при изменении ситуационных планов размещения антенн, технических характеристик или режимов работы объекта;

- после проведения защитных мероприятий, направленных на снижение уровней ЭМП;

- в порядке плановых контрольных измерений (не реже одного раза в год).

9.2. Подготовка к проведению инструментального контроля включает несколько этапов:

- согласование цели, времени и условий проведения измерений с заинтересованными организациями и предприятиями;

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

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

- изучение нормативной санитарно-гигиенической документации, относящейся к данному виду технических средств;

- обеспечение измерений дальности до точки измерений;

- подбор и подготовку необходимой измерительной аппаратуры;

- выполнение расчетного прогнозирования;

- определение необходимости использования средств индивидуальной защиты.

9.3. Для инструментального контроля используются измерители плотности потока энергии (ППЭ), к которым предъявляются следующие требования:

- измерительная антенна (зонд, датчик поля) не должна существенно искажать структуру измеряемого поля;

- измерительная антенна (зонд, датчик поля) должна иметь изотропную характеристику направленности;

- кабель, соединяющий антенну с измерительным прибором не должен обладать антенным эффектом;

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

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

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

9.5. Выбор трасс (маршрутов) измерений.

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

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

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

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

- наклон трассы должен быть минимальным по сравнению с наклоном всех возможных трасс в данном секторе;

- трасса должна быть доступной для пешего передвижения или для автотранспорта;

- протяженность трассы определяется на основе расчетного удаления границ СЗЗ и глубины зоны ограничения застройки (в 1,5 - 2 раза больше);

- точки (площадки) для проведения измерений следует выбирать с интервалом 5 - 10 м - при удалении до 100 м от излучающей антенны; 10-20 м - при удалении от 100 м до 300 м; 20-50 м - при удалении более 300 м.

9.6. Проведение измерений.

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

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

9.6.3. При измерении уровней ЭМП необходимо ориентировать измерительную антенну (зонд, датчик поля) на излучающую антенну и вращением ее вокруг своей оси добиться максимального значения.

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

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

9.6.6. По результатам измерений составляется протокол, который является неотъемлемой частью санитарного паспорта.

Главный государственный санитарный врач Российской Федерации - Первый заместитель Министра здравоохранения Российской Федерации

Г.Г. Онищенко