МУК 4.3.678-97 «Определение уровней напряжений, наведенных электромагнитными полями на проводящие элементы зданий и сооружений в зоне действия мощных источников радиоизлучений»

ГОСУДАРСТВЕННАЯ СИСТЕМА САНИТАРНО-ЭПИДЕМИОЛОГИЧЕСКОГО
НОРМИРОВАНИЯ РОССИЙСКОЙ ФЕДЕРАЦИИ

4.3. МЕТОДЫ КОНТРОЛЯ. ФИЗИЧЕСКИЕ ФАКТОРЫ.

ОПРЕДЕЛЕНИЕ УРОВНЕЙ
НАПРЯЖЕНИЙ, НАВЕДЕННЫХ
ЭЛЕКТРОМАГНИТНЫМИ ПОЛЯМИ
НА ПРОВОДЯЩИЕ ЭЛЕМЕНТЫ ЗДАНИЙ
И СООРУЖЕНИЙ В ЗОНЕ ДЕЙСТВИЯ
МОЩНЫХ ИСТОЧНИКОВ
РАДИОИЗЛУЧЕНИЙ

МЕТОДИЧЕСКИЕ УКАЗАНИЯ

МУК 4.3.678-97

Минздрав России

Москва

1998

1. Разработаны Бузовым А.Л., Романовым В.А., Казанским Л.С., Кольчугиным Ю.И., Юдиным В.В. (Самарский отраслевой научно-исследовательский институт радио Министерства связи Российской Федерации).

2. Представлены Госкомсвязи России письмом от 27.05.97 № НТУОТ-1/058.

Одобрены к утверждению Комиссией по государственному санитарно-эпидемиологическому нормированию при Минздраве России.

3. Утверждены и введены в действие Главным государственным санитарным врачом Российской Федерации от 6 ноября 1997 г.

4. Введены впервые.

СОДЕРЖАНИЕ

УТВЕРЖДАЮ

Главный государственный

санитарный врач Российской Федерации

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

06 ноября 1997 г.

МУК 4.3.678 - 97

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

4.3. МЕТОДЫ КОНТРОЛЯ. ФИЗИЧЕСКИЕ ФАКТОРЫ.

ОПРЕДЕЛЕНИЕ УРОВНЕЙ
НАПРЯЖЕНИЙ, НАВЕДЕННЫХ
ЭЛЕКТРОМАГНИТНЫМИ ПОЛЯМИ
НА ПРОВОДЯЩИЕ ЭЛЕМЕНТЫ ЗДАНИЙ
И СООРУЖЕНИЙ В ЗОНЕ ДЕЙСТВИЯ
МОЩНЫХ ИСТОЧНИКОВ
РАДИОИЗЛУЧЕНИЙ

Методические указания

1. ОБЛАСТЬ ПРИМЕНЕНИЯ

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

2. СУЩНОСТЬ МЕТОДА

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

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

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

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

3. ОСНОВНЫЕ ПОЛОЖЕНИЯ МЕТОДИКИ РАСЧЕТА НАВЕДЕННЫХ НАПРЯЖЕНИЙ

Расчет выполняется в следующей последовательности:

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

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

- составляется и решается СЛАУ - аналог интегрального уравнения; найденные коэффициенты разложения совместно с базисными функциями восстанавливают токовую функцию, аппроксимирующую истинное распределение тока;

- по найденному распределению тока с учетом поля излучения антенн радиостанций рассчитываются наведенные напряжения.

3.1. РАСЧЕТ ПРОВОЛОЧНОЙ МОДЕЛИ ОБЪЕКТА

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

Исследуемый объект представляется в виде системы цилиндрических проводников. Совокупность осей проводников представляет собой контур L¢ (не обязательно гладкий и непрерывный), на котором определяется токовая функция I (l¢), где l¢ - криволинейная координата, отсчитываемая вдоль L¢ (каждой точке на проводниках модели должно однозначно соответствовать некоторое значение l¢). Для наложения граничных условий строится контур L, который представляет собой контур L¢, перенесенный на поверхность проводников (L и L¢ нигде не пересекаются и не соприкасаются). На контуре L¢ выделяются N коротких отрезков - сегментов, как показано на рис. 1. Каждый k-й сегмент определяется тремя точками: l¢ = ak - начало, l¢ = bk - средняя точка, l¢ = ck - конец. Отрезки [ak, bk] и [bk, ck] - соответственно 1-е и 2-е плечи сегмента (k-го) - могут не лежать на одной прямой и иметь разную длину. Соседние сегменты частично перекрываются: средняя точка k-го сегмента bk совпадает с концом (k-1)-го ck-1 и началом (k+1)-го ak+1 сегментов. Электрические соединения описываются введением дополнительных сегментов, плечи каждого из которых лежат на разных проводах, как показано на рис. 2.

Длина плеча сегмента не должна превышать 10 % длины волны, отношение этой длины к радиусу провода должно находиться в пределах 3 … 10.

Каждому сегменту соответствует одна точка сшивания. В качестве точек сшивания берутся точки контура L, наименее удаленные от средних точек соответствующих сегментов.

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

3.2. РАСЧЕТ ПОЛЯ ИЗЛУЧЕНИЯ АНТЕНН РАДИОСТАНЦИЙ

Поле излучения антенн радиостанций в некоторой точке наблюдения вычисляется следующим образом.

Антенны рассматриваемых диапазонов представляют собой системы линейных проводников (вибраторов, проводников рефлекторов антенн СГД, проводников ромбических антенн). Каждый проводник представляется как совокупность элементарных электрических вибраторов (ЭЭВ).

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

Рассматриваются два типа ЭЭВ - горизонтальные и вертикальные. Практически все многообразие конструкций антенн НЧ, СЧ и ВЧ диапазонов исчерпывается данными типами ЭЭВ. Ниже для указания горизонтального и вертикального ЭЭВ будут использоваться надстрочные индексы «г» и «в», соответственно

Горизонтальный ЭЭВ ориентируется вдоль оси ОХ вспомогательной декартовой системы (отличающейся от основной только поворотом на некоторый угол α относительно оси аппликат). Высота ЭЭВ относительно плоскости XOY (которая совпадает с поверхностью земли) h. Компоненты разложения по базисным ортам вспомогательной системы электрического поля в некоторой точке наблюдения с координатами х, у, z определяются выражениями:

                          (3.1)

где М = jk1Ilээв - дипольный момент;

k1 = 2p/l - волновое число для воздуха;

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

I - амплитуда тока ЭЭВ, А;

lээв - длина ЭЭВ, м;

R1 = (х2 + у2 + (z - h)2)1/2 - расстояние между точкой наблюдения и центром ЭЭВ.

Размещение сегментов на контуре L

Рис. 1

Описание электрических соединений в проволочной модели

Рис. 2

Функции ех, еу, еz в (3.1) имеют вид:

ex = erxcosj - efxsinj;

ey = erxcosj + efxsinj;                                                   (3.2)

ez = esxcosj,

где

erx = [1 + arr - D(1 + brr + 2idbs) - 2d2DI(d) - 2WU(d) + 2WU(1/d)]cosj,

efx = {-[1 + a¢rr - D -Db¢r + 2DI(1/d) - 2WU(1/d) + 2WU(d)]}sinj,                  (3.3)

esx = ars - Dbrs + 2idDU(d).

Вертикальный ЭЭВ ориентируется вдоль оси OZ вспомогательной цилиндрической системы и расположен на высоте h над плоскостью z = 0, которая совпадает с поверхностью земли. Компоненты разложения по базисным ортам цилиндрической системы электрического поля в некоторой точке наблюдения с координатами r, z (в силу симметрии задачи азимут j может быть произвольным) определяются выражениями:

                                          (3.4)

где функции еr и еz имеют вид:

er = ars +Dbrs -2idU(d).

ez = 1 + ass + D(-1 + bss - 2idbs) + 2D(1-d2)I(d).                           (3.5)

Параметры, входящие в (3.3 - 3.5), вычисляются по формулам:

arr = acos2q - bsin2q - 1,

brr = acos2q¢ - bsin2q¢ - 1,

ars = -(a + b)sinq×cosq,

brs = -(a¢ + b¢)sinq¢cosq¢,

ass = asin2q - bcos2q - 1,

bss = asin2q¢ - bcos2q¢ - 1,

ar = i(1 + i/k1R1)sinq,                 a¢r = ar/k1r,

br = i(1 + i/k1R2)sinq¢,                b¢r = br/k1r,

bs = id¢cosq,                       d¢ = 1 + i/k1R2,                                                                      (3.6)

W = d2D/k1r(1 - d2),

D = R1exp[-ik(R1 - R2)]/R2,

d = 1/(e + i60ls + 1)1/2,

a = 1 + i/k1R1 - 1/k12R12,                  a¢ = 1 + i/k1R2 - 1/k12R22,

b = 2i(1 + i/k1R1)/k1R1,                    b¢ = 2i(1 + i/k1R2)/k1R2,

cosq = (z - h) /R1,                            cosq¢ = (z - h)/R2.

где r, R2 - геометрические параметры задачи, выражающиеся формулами:

R2 = (x2 + у2 + (z + h)2)1/2,       r = (х2 + у2)1/2.                                 (3.7)

Вспомогательные функции I(d) и U(d) выражаются через функцию ослабления у(z, r):

I(d) = у(z, r),                      U(d) = r(i - 1/k1R,)у(z,r)/R1.                               (3.8)

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

            (3.9)

где S0 = ik1R2d2/2(r/R2)2,                             S = S0[1 + (z + h)/dR2]2.

После вычисления поля ЭЭВ как вектора в соответствующей вспомогательной системе координат должны рассчитываться компоненты его разложения по ортам основной системы. Поле излучения Еиз (вектор в основной системе) находится как векторная сумма полей отдельных ЭЭВ.

3.3. РАСЧЕТ РАСПРЕДЕЛЕНИЯ ТОКА, НАВЕДЕННОГО НА ПРОВОДЯЩИЕ ЭЛЕМЕНТЫ ОБЪЕКТА

Распределение тока определяется как токовая функция I(l¢) полностью определяемая коэффициентами разложения по кусочно-синусоидальному базису I1, I2, I3, … IN (N - число сегментов).

Коэффициенты I1, I2, … IN находятся решением СЛАУ:

                                                  (3.10)

где Zik - коэффициенты, образующие квадратную матрицу импедансов;

Eik - свободные члены СЛАУ.

Некоторый коэффициент Zik выражает связь между k-м сегментом и i-й точкой сшивания (физический смысл - тангенциальное поле в i-й точке сшивания при Ik = 1 А) и вычисляется по формуле:

Zik = Z¢ik + kzZ"ik,                                                               (3.11)

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

kz = 1, если объект содержит нижнее (на уровне земной поверхности) железобетонное перекрытие, в противном случае kz = 0.

Величины Z¢ik и Z"ik рассчитываются следующим образом.

На рис. 3 показаны произвольные i-я точка сшивания и k-й сегмент, с каждым плечом которого связана соответствующая вспомогательная цилиндрическая система координат. Сначала вычисляются z-я и r-я компоненты поля (при Ik = 1 А) в точке сшивания поочередно от каждого плеча сегмента. Формулы для z-й и r-й компонент поля, создаваемого в точке сшивания одним из плеч сегмента (в соответствующей цилиндрической системе, - см. рис. 3) имеют вид:

                      (3.12)

                   (3.13)

где r1 - расстояние до точки наблюдения от начала (конца) сегмента, м;

r0 - расстояние до точки наблюдения от средней точки сегмента, м;

b = 2p/l - волновое число;

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

l - длина рассматриваемого плеча, м;

z и r - цилиндрические координаты точки наблюдения (соответственно аппликата и проекция радиус-вектора точки на плоскость z = 0, м).

Знак «+» в (3.11, 3.12) соответствует 1-му плечу сегмента, знак «-» -2-му.

К выводу формул для расчета элементов матрицы импедансов

(ц)Z1(ц)) и (О(ц)Z2(ц)) - оси аппликат вспомогательных цилиндрических систем для 1-го и 2-го плеч k-го сегмента, соответственно 1z1, 1r1, 1z2, 1r2 - единичные базисные векторы вспомогательных цилиндрических систем

Рис. 3

Пусть z- и r-компоненты поля рассчитаны для обоих плеч k-го сегмента, т.е. получены 4 числа. Обозначим их Em,k, m = 1, 2, 3, 4. Каждой m-й компоненте в исходной основной системе координат соответствует единичный вектор 1¢m,k. С учетом этих обозначений формула для Z¢ik имеет вид:

                                                      (3.14)

где 1i = (сi - ai)/|ci - аi| - единичный тангенциальный вектор в i-й точке сшивания (сi и аi - радиус-векторы конца и начала i-го сегмента).

Величина Z²ik вычисляется аналогичным образом для зеркального изображения k-го сегмента.

Свободный член i-го уравнения СЛАУ (соответствующего i-й точке сшивания) Еi рассчитывается по формуле:

Ei = (1i, Eиз)×Рст                                                             (3.15)

где Еиз - поле излучения радиостанций в i-й точке сшивания, В/м, вычисленное методом, описанным в п. 3.2.

Рст - коэффициент, учитывающий влияние стен здания (для железобетонных стен Рст = 0,316, во всех остальных случаях Рст = 1).

3.4. РАСЧЕТ НАВЕДЕННЫХ НАПРЯЖЕНИЙ

Наведенное напряжение между некоторыми точками с радиус-векторами v01, и v02 вычисляется по формуле:

                                           (3.16)

где L0 - контур интегрирования - прямолинейный отрезок, соединяющий точки v01, и v02;

l0 - переменная интегрирования - криволинейная координата, отсчитываемая вдоль L0;

v0 - радиус-вектор точки l0 (в основной системе);

10 = (v02 - v01)/|v02 - v01|.

Еиз - поле излучения радиостанций в точке v0, В/м, вычисленное методом, описанным в п. 3.2.

Коэффициенты Z в (3.16) вычисляются также как и коэффициенты матрицы импедансов (п. 3.3), только вместо точки сшивания берется точка l0 (с радиус-вектором v0).

4. МЕТОДИКА ИЗМЕРЕНИЙ УРОВНЕЙ НАВЕДЕННЫХ НАПРЯЖЕНИЙ

4.1. СОСТАВ ИЗМЕРИТЕЛЬНОГО КОМПЛЕКСА

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

- селективный микровольтметр SMV 11;

- антенна для измерения напряженности поля «ОРТ» (разработана Самарским отраслевым НИИ радио);

- милливольтметр ВЗ-48А (широкополосный измеритель среднеквадратического значения напряжения произвольной формы в диапазоне 10 Гц … 50 МГц);

- металлический лист МЛ круговой формы диаметром 1 м, имитирующий тело человека как емкость относительно проводящего пола или реальной земли;

- соединительная шина, обеспечивающая соединение корпуса пробника прибора ВЗ-48А с МЛ;

- диэлектрическая штанга, осуществляющая надежную изоляцию тела оператора от металлоконструкции (штанга прикрепляется к корпусу пробника прибора ВЗ-48А).

4.2. МЕТОДИКА ИЗМЕРЕНИЙ

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

4.2.1. Комплектом в составе SMV 11 и антенны «ОРТ» (измеритель напряженности поля) выполняются измерения напряженности поля вблизи металлоконструкции (в местах, доступных для прикосновения). При этом возможны два случая:

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

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

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

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

Во втором случае измерителем напряженности поля находятся точки максимумов поля (потенциала) на частоте станции с превалирующим уровнем сигнала. Фиксируются показания прибора в этих точках, наибольшее значение Vm принимается равным 0 дБ (остальные нормируются по этому уровню). Затем в этих же точках определяются относительные (нормированные к Vm) уровни сигнала на частотах остальных радиостанций Vij, дБ (i - номер радиостанции, j - номер точки измерения). Все показания прибора фиксируются.

Схема измерений наведенных напряжений

Рис. 4

В третьем случае определяется относительное среднеквадратическое значение как функция криволинейной координаты х, отсчитываемой вдоль металлоконструкции:

                                                     (4.1)

где i - номер радиостанции; N - общее число радиостанций.

Максимумы функции Vs(x) будут точками, в которых необходимо выполнить измерение напряжения. Эти точки фиксируются, фиксируются также показания прибора в этих точках Vij, нормированные к максимальному измеренному уровню Vm.

4.2.2. В найденных точках максимумов напряжения производится его измерение прибором ВЗ-48А согласно схеме рис. 4, предварительно удалив покрытия (например, краску) с поверхности металла в этих точках. При измерениях необходимо щуп прибора плавно приближать к точке касания, контролируя показания прибора. Если измеренное значение оказывается равным 300 В до момента касания, измерение в данной точке прекращается, и делается вывод о превышении в ней уровня 300 В. В противном случае осуществляется касание щупом поверхности металла и фиксируется показание прибора - среднеквадратическое значение Ucj, В (j - номер точки измерения).

4.2.3. Если радиостанция одна, измерения прекращаются. Значение наведенного напряжения для j-й точки наблюдения берется равным Ucj, В. Если радиостанций несколько для каждой j-й точки измерения вычисляются значения наведенных напряжений Uij, В по формулам:

                                (4.2)

где k - номер радиостанции, уровень сигнала которой превышает все остальные.

4.2.4. При проведении измерений прибором ВЗ-48А по схеме, показанной на рис. 4, необходимо обеспечить надежный электрический контакт между соединительной шиной, корпусом пробника прибора и металлическим листом МЛ.

4.2.5. Результаты измерений и расчетов заносятся в протокол.