Вход
Купить систему ГАРАНТ
Получить демо-доступ
Узнать стоимость
Информационный банк
Подобрать комплект
Семинары
Методические указания МУК 4.3.678-97 "Определение уровней напряжений, наведенных электромагнитными полями на проводящие элементы зданий и сооружений в зоне действия мощных источников радиоизлучений" (утв. Главным государственным санитарным врачом РФ 6 ноября 1997 г.)
Методические указания МУК 4.3.678 - 97
"Определение уровней напряжений, наведенных электромагнитными полями на проводящие элементы зданий и сооружений в зоне действия мощных источников радиоизлучений"
(утв. Главным государственным санитарным врачом РФ 6 ноября 1997 г.)
Дата введения: с момента утверждения
Содержание
1. Область применения
2. Сущность метода
3. Основные положения методики расчета наведенных напряжений
3.1. Расчет проволочной модели объекта
3.2. Расчет поля излучения антенн радиостанций
3.3. Расчет распределения тока, наведенного на проводящие элементы
объекта
3.4. Расчет наведенных напряжений
4. Методика измерений уровней наведенных напряжений
4.1. Состав измерительного комплекса
4.2. Методика измерений
1. Область применения
Методические указания составлены в помощь инженерам органов и учреждений санитарно-эпидемиологической службы, инженерно-техническим работникам, проектным организациям средств связи с целью обеспечения предупредительного санитарного надзора за источниками излучения кило-(НЧ), гекто-(СЧ) и декаметрового (ВЧ) диапазонов технических средств радиовещания и радиосвязи, а также для прогнозирования уровней напряжений, наведенных электромагнитными полями на проводящие элементы зданий и сооружений в зоне действия мощных источников радиоизлучений этих диапазонов.
2. Сущность метода
Расчетный метод основан на предварительном расчете распределения тока, наведенного на проводящие элементы (металлоконструкции) объекта (здания или иного сооружения) полем излучения антенн радиостанций.
Задача о нахождении наведенного тока решается как внешняя электродинамическая задача методом интегрального уравнения в тонкопроволочном приближении. Исследуемый объект представляется в виде системы цилиндрических проводников. В соответствии с тонкопроволочным приближением, во-первых, считается, что ток течет по оси проводника (т.е. поверхностная плотность тока как векторная функция пространственных координат заменяется линейным осевым током - скалярной функцией криволинейной координаты, отсчитываемой вдоль оси проводника), а во-вторых, на поверхности проводника учитывается только тангенциальная составляющая электрического поля, параллельная его оси.
Для решения интегрального уравнения используется разновидность метода моментов - метод сшивания в дискретных точках при кусочно-синусоидальном базисе разложения искомой токовой функции. Интегральное уравнение имеет смысл граничного условия на поверхности проводника, которое выражается равенством нулю тангенциальной составляющей полного электрического поля - суперпозиции поля излучения антенн радиостанций и поля, создаваемого излучением наведенных токов. В рамках данного метода граничное условие накладывается в дискретных точках - точках сшивания. Интегральное уравнение сводится к системе линейных алгебраических уравнений (СЛАУ), решаемой относительно коэффициентов разложения токовой функции по кусочно-синусоидальному базису. Число уравнений (неизвестных) СЛАУ, а также число базисных функций равно числу точек сшивания, размещенных вдоль проводников объекта.
Поле излучения радиостанций рассматривается как стороннее (не зависящее от параметров объекта) и рассчитывается (в точках сшивания) с учетом электрофизических параметров реальной почвы, экранирующего и поглощающего действия стен зданий, а также подстилающей поверхности (нижнего перекрытия). Наведенные напряжения рассчитываются как результат действия полей излучения антенн радиостанций и наведенных токов.
3. Основные положения методики расчета наведенных напряжений
Расчет выполняется в следующей последовательности:
- рассчитываются геометрические параметры проволочной модели исследуемого объекта (пространственные координаты точек сшивания и сегментов - коротких отрезков проводников, на которых определяются соответствующие базисные функции);
- в точках сшивания рассчитываются тангенциальные составляющие стороннего поля - поля излучения антенн радиостанций с учетом электрофизических параметров реальной почвы, влияния стен здания и подстилающей поверхности;
- составляется и решается СЛАУ - аналог интегрального уравнения; найденные коэффициенты разложения совместно с базисными функциями восстанавливают токовую функцию, аппроксимирующую истинное распределение тока;
- по найденному распределению тока с учетом поля излучения антенн радиостанций рассчитываются наведенные напряжения.
3.1. Расчет проволочной модели объекта
Задается основная декартова система пространственных координат. Относительно исследуемого объекта введенную систему целесообразно ориентировать таким образом, чтобы оси абсцисс и ординат лежали в плоскости земной поверхности, ось аппликат была направлена вертикально вверх.
Исследуемый объект представляется в виде системы цилиндрических проводников. Совокупность осей проводников представляет собой контур L' (не обязательно гладкий и непрерывный), на котором определяется токовая функция I(l'), где l' - криволинейная координата, отсчитываемая вдоль L' (каждой точке на проводниках модели должно однозначно соответствовать некоторое значение l'). Для наложения граничных условий строится контур L, который представляет собой контур L', перенесенный на поверхность проводников (L и L' нигде не пересекаются и не соприкасаются). На контуре L' выделяются N коротких отрезков - сегментов, как показано на рис.1. Каждый k-й сегмент определяется тремя точками: l'=аk - начало, l'=bk - средняя точка, l'=сk - конец. Отрезки [аk, bk] и [bk, ck] - соответственно 1-е и 2-е плечи сегмента (k-го) - могут не лежать на одной прямой и иметь разную длину. Соседние сегменты частично перекрываются: средняя точка k-го сегмента bk совпадает с концом (k-1)-ro сk-1, и началом (k+1)-ro ak+1, сегментов. Электрические соединения описываются введением дополнительных сегментов, плечи каждого из которых лежат на разных проводах, как показано на рис.2.
Длина плеча сегмента не должна превышать 10% длины волны, отношение этой длины к радиусу провода должно находиться в пределах 3...10.
Каждому сегменту соответствует одна точка сшивания. В качестве точек сшивания берутся точки контура L, наименее удаленные от средних точек соответствующих сегментов.
В результате расчетов на данном этапе должны быть определены пространственные координаты (в основной системе) сегментов и точек сшивания.
3.2. Расчет поля излучения антенн радиостанций
Поле излучения антенн радиостанций в некоторой точке наблюдения вычисляется следующим образом.
Антенны рассматриваемых диапазонов представляют собой системы линейных проводников (вибраторов, проводников рефлекторов антенн СГД, проводников ромбических антенн). Каждый проводник представляется как совокупность элементарных электрических вибраторов (ЭЭВ).
Ток каждого ЭЭВ полностью определяется входным сопротивлением, мощностью передатчика и законом распределения тока. В качестве последнего следует задаться синусоидальным распределением тока по проводникам вибраторных антенн и режимом бегущей волны тока в проводниках ромбических антенн. Поле в точке наблюдения находится как суперпозиция полей, создаваемых всеми ЭЭВ антенны. Задача, таким образом, сводится к нахождению поля одного ЭЭВ.
Рассматриваются два типа ЭЭВ - горизонтальные и вертикальные. Практически все многообразие конструкций антенн НЧ, СЧ и ВЧ диапазонов исчерпывается данными типами ЭЭВ. Ниже для указания горизонтального и вертикального ЭЭВ будут использоваться надстрочные индексы "г" и "в", соответственно.
Горизонтальный ЭЭВ ориентируется вдоль оси ОХ вспомогательной декартовой системы (отличающейся от основной только поворотом на некоторый угол альфа относительно оси аппликат). Высота ЭЭВ относительно плоскости XOY (которая совпадает с поверхностью земли) h. Компоненты разложения по базисным ортам вспомогательной системы электрического поля в некоторой точке наблюдения с координатами х, у, z определяются выражениями:
M ex e(jk1R1) M ey e(jk1R1) M ez e(jk1R1)
Ex(г)=---------------, Ey(г)=--------------, Ez(г)=--------------,
R1 R1 R1 (3.1)
где М=jk1 Il ЭЭВ - дипольный момент;
k1 = 2пи/ламбда - волновое число для воздуха;
ламбда - длина волны, м;
I - амплитуда тока ЭЭВ, А;
lээв - длина ЭЭВ, м;
R1=(х(2)+y(2)+(z-h)(2))(1/2) - расстояние между точкой наблюдения и центром ЭВВ.
Размещение сегментов на контуре L'
|L'
|
(Сk+1) |
--------------------------+---------------------------------\
| |
| |
| |
(bk+1) | |(k+1)-й
----------------Ck--------+--------------\ }сегмент
| | |
| | |
| | |
(Сk-1=ak+1) | | |
-----------------bk-------+--------------}k-й сегмент ------|
| | |
| | |
| | |
(bk-1) | | |(k-1)-й
----------------ak--------+--------------/ }сегмент
| |
| |
| |
(аk+1) | |
--------------------------+---------------------------------/
|
|
|
|
Рис.1
Описание электрических соединений в проволочной модели
Провод 2 --| |
|
| | |
|------Крайний сегм.пров.2
| | |
|
| | |
-
| ||| |
|||
| ||| |
|||
| ||| |
|||
Провод 1 | ||| |
- - - - -/- - - - - - - - - / |||
||| |
|||
Крайний сегм.пров.1 ||| |
------------------------ -
| |
|
- - - - - - - - - -- - - - - - \ | |
| |-----Крайний сегм.пров.3
| |
| |
Провод 3---| |
| |
Рис.2
Функции ex, ey, ez в (3.1) имеют вид:
ex=erx cos фи - efx sin фи;
ey=erx sin фи - efx cos фи;
ez=esx cos фи, (3.2)
где
erx=[1+arr-ДЕЛЬТА(1+brr+2i сигма bs) - 2дельта(2) ДЕЛЬТА I(дельта) -
- 2ОМЕГА U(дельта) + 2ОМЕГА(U)(1/дельта)] cos фи,
efx={-[1+a'rr-ДЕЛЬТА-ДЕЛЬТАb'r + 2ДЕЛЬТАI(1/дельта) -
- 2ОМЕГА U(1/дельта) + 2 ОМЕГА U (дельта)]} sin фи, (3.3)
esx=ars - ДЕЛЬТА brs + 2iдельта ДЕЛЬТА U (дельта).
Вертикальный ЭЭВ ориентируется вдоль оси OZ вспомогательной цилиндрической системы и расположен на высоте h над плоскостью z=0, которая совпадает с поверхностью земли. Компоненты разложения по базисным ортам цилиндрической системы электрического поля в некоторой точке наблюдения с координатами р, z (в силу симметрии задачи азимут фи может быть произвольным) определяются выражениями:
Er(B) = (M er e(jk1 R1))/R,
Ez(B) = (M ez e(jk1 R1))/R, (3.4)
где функции еr и еz имеют вид:
er=ars+ДЕЛЬТА brs - 2i дельта U(дельта), (3.5)
ez=1+ass+ДЕЛЬТА(-1+bss-2i дельта bs) +
+ 2ДЕЛЬТА(1-дельта(2)) I(дельта).
Параметры, входящие в (3.3-3.5), вычисляются по формулам:
arr=a cos(2)тета' - b sin(2) тета - 1,
brr=a cos(2)тета` - b sin(2) тета` - 1,
ars=-(a+b)sin тета cos тета,
brs=-(a`+b`)sin тета` cos тета`,
ass=a sin(2)тета - b cos(2) тета - 1,
bss=a sin(2)тета` - b cos(2) тета` - 1,
ar=i(1+i/k1R1) sin тета, a`r=ar/k1r,
br=i(1+i/k1R2) sin тета`, b`r=br/k1r,
bs=id` cos тета, d`=1+i/k1R2, (3.6)
ОМЕГА=дельта(2)ДЕЛЬТА/k1r(1-дельта(2)),
ДЕЛЬТА=R1 exp[-ik(R1-R2)]/R2,
дельта=1/(эпсилон+i60 ламбда сигма + 1)(1/2),
а=1+i/k1R1 - 1/k1(2)R1(2), a`=1+i/k1R2 - 1/k1(2)R2(2),
b=2i(1+i/k1R1)/k1R1, b`=2i(1+i/k1R2)/k1R2,
cos тета=(z-h)/R1, cos тета`=(z-h)/R2.
где r, R2 - геометрические параметры задачи, выражающиеся формулами:
R2=(x(2)+y(2)+(z+h)(2))(1/2), r=(x(2)+(y(2))(1/2). (3.7)
Вспомогательные функции I(дельта) и U(дельта) выражаются через
функцию ослабления у (z,r):
I(дельта)=y(z,r), U(дельта)=r(i-1/k1R1)y(z,r)/R1. (3.8)
Для вычисления функции ослабления, в которую входит интеграл вероятности от комплексного аргумента, используются сходящиеся и асимптотические разложения:
y(z,r)=1+j кв.корень[пи S0 E(-s)] - 2 кв.корень[SS0] x
n (-2S)(v)
x сумма ---------, |S|<12;
v=0 (2v+1)!!
n
y(z,r)=1-j кв.корень[S0/S] сумма(2v-1)!!/(2S)(v),
v=0
|S|>=12 и Im кв.корень[S]>=0;
y(z,r)=1+2j кв.корень[пи S0 e(-s)] - кв.корень[S0/S] x
n
x сумма(2v-1)!!/(2S)(v), |S|>=12 и Im кв.корень[S]<0;
v=0
где S0=ik1 R2 дельта(2)/2 (r/R2)(2), S=S0[1+(z+h)/дельта R2](2).
После вычисления поля ЭЭВ как вектора в соответствующей вспомогательной системе координат должны рассчитываться компоненты его разложения по ортам основной системы. Поле излучения Еиз (вектор в основной системе) находится как векторная сумма полей отдельных ЭЭВ.
3.3. Расчет распределения тока, наведенного
на проводящие элементы объекта
Распределение тока определяется как токовая функция I(l') полностью определяемая коэффициентами разложения по кусочно-синусоидальному базису I1, I2, I3,... IN (N - число сегментов).
Коэффициенты I1, I2, I3,... IN находятся решением СЛАУ:
N
сумма Zik Ik = Ei, i=1, 2, ... N, (3.10)
k=1
где Zik - коэффициенты, образующие квадратную матрицу импедансов;
Еik - свободные члены СЛАУ.
Некоторый коэффициент 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-я компоненты поля (при ek=1 А) в точке сшивания поочередно от каждого плеча сегмента. Формулы для z-й и рo-й компонент поля, создаваемого в точке сшивания одним из плеч сегмента (в соответствующей цилиндрической системе, - см. рис.3) имеют вид:
exp(-j бета r1) exp(-j бета r0)
Ez=+-j30[--------------- - cos(бета l) --------------- -
r1 r0
- z sin (бетаl)(1/(бета r0(3)) + j/r0(2)) exp(-j бета r0)], (3.12)
exp(-j бета r1) exp(-j бета r0)
Epo=+-j 30[- --------------- - cos(бета l) --------------- -
po r1(z-1)(-1) r0 z(-1)
exp(-j бета r0)
- sin (бетаl)(1-z(2)/r0(2) - j бета z(2)/r0)--------------- ].
бета r0
(3.13)
где r1 - расстояние до точки наблюдения от начала (конца) сегмента, м;
r0 - расстояние до точки наблюдения от средней точки сегмента, м;
бета=2пи/ламбда - волновое число;
ламбда - длина волны, м;
l - длина рассматриваемого плеча, м;
z и р - цилиндрические координаты точки наблюдения (соответственно аппликата и проекция радиус-вектора точки на плоскость z=0, м).
Знак "+" в (3.11, 3.12) соответствует 1-му плечу сегмента, знак "-" - 2-му.
К выводу формул для расчета элементов матрицы импедансов
(О(ц) Z1(ц)) и (О(ц) Z2(ц)) - оси аппликат вспомогательных цилиндрических систем для 1-го и 2-го плеч k-гo сегмента, соответственно 1z1, 1p1, 1z2, 1p2 - единичные базисные векторы вспомогательных цилиндрических систем
-1z1, 1po2
|
|
|
i-я точка |
сшивания +------------ 1р1, 1z2
О(ц)
/=======>========== --------------- Z2(ц)
|| 2-е плечо
||
||
|| к-й сегмент
||
||1-е
||плечо
||
||
||
||
||
||
||
|
|
|
|
|
Z1(ц)
Рис.3
Пусть z- и рo-компоненты поля рассчитаны для обоих плеч k-го сегмента, т.е. получены 4 числа. Обозначим их Еm,k, m = 1, 2, 3, 4. Каждой m-й компоненте в исходной основной системе координат соответствует единичный вектор 1'm,k. С учетом этих обозначений формула для Z'ik имеет вид:
4
Zik=сумма(1i, 1'm,k) Em,k, (3.14)
m=1
где 1i=(ci-ai)/|ci-ai| - единичный тангенциальный вектор в i-й точке сшивания (сi и аi - радиус-векторы конца и начала i-го сегмента).
Величина Z"ik вычисляется аналогичным образом для зеркального изображения k-го сегмента.
Свободный член i-го уравнения СЛАУ (соответствующего i-й точке сшивания) Еi рассчитывается по формуле:
Еi=-(1i,Еиз)*Рст (3.15)
где Еиз - поле излучения радиостанций в i-й точке сшивания, В/м, вычисленное методом, описанным в п.3.2.
Рст - коэффициент, учитывающий влияние стен здания (для железобетонных стен Рст=0,316, во всех остальных случаях Рст=1).
3.4. Расчет наведенных напряжений
Наведенное напряжение между некоторыми точками с радиус-векторами v(0)1 и v(0)2 вычисляется по формуле:
N
U=-интеграл[сумма Z(0)k Ik + (1(0), Eиз) Рст] dl(0) (3.16)
L0 k=1
где L0 - контур интегрирования - прямолинейный отрезок, соединяющий точки v(0)1 и v(0)2;
l(0) - переменная интегрирования - криволинейная координата, отсчитываемая вдоль L(0);
v(0) - радиус-вектор точки l(0) (в основной системе);
1(0)=(v(0)2-v(0)1)/|v(0)2-v(0)1|.
Еиз - поле излучения радиостанций в точке v(0), В/м, вычисленное методом, описанным в п.3.2.
Коэффициенты Z в (3.16) вычисляются также как и коэффициенты матрицы импедансов (п.3.3), только вместо точки сшивания берется точка I(0) (с радиус-вектором v(0)).
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 - номер точки измерения). Все показания прибора фиксируются.
Схема измерений наведенных напряжений
Металлическая часть здания или другого сооружения
-- -- -- --------------------------------------------------- -- -- ---
-- -- -- --------------|------------------------------------ -- -- ---
|-Высокочастотный щуп
/---\
/-------\ |
| | |
| /---/ |
| | | |
| | \---/
| | ||
Шина-| | ||
соедините- | | ||
льная | | || /----------------\
| | || | Прибор ВЗ-48А |
| | |\----| |
| | \-----| |
| \---\ | |
| | | |
/---------------------------------------------------------\
| |
\---------------------------------------------------------/
|
Металлический лист,
имитирующий тело человека
---------------------------------------------------------------------
| |
Токопроводящий пол --------
или земная поверхность ------
--
Рис.4
В третьем случае определяется относительное среднеквадратическое значение как функция криволинейной координаты х, отсчитываемой вдольметаллоконструкции:
N 10Vi(x) 0,5
V (x) = [сумма --------] (4.1)
s i=l 10
где i - номер радиостанции; N - общее число радиостанций.
Максимумы функции Vs (х) будут точками, в которых необходимо выполнить измерение напряжения. Эти точки фиксируются, фиксируются также показания прибора в этих точках Vij, нормированные к максимальному измеренному уровню Vm.
4.2.2. В найденных точках максимумов напряжения производится его измерение прибором ВЗ-48А согласно схеме рис.4., предварительно удалив покрытия (например, краску) с поверхности металла в этих точках. При измерениях необходимо щуп прибора плавно приближать к точке касания, контролируя показания прибора. Если измеренное значение оказывается равным 300 В до момента касания, измерение в данной точке прекращается, и делается вывод о превышении в ней уровня 300 В. В противном случае осуществляется касание щупом поверхности металла и фиксируется показание прибора - среднеквадратическое значение U(c)j, В (j- номер точки измерения).
4.2.3. Если радиостанция одна, измерения прекращаются. Значение наведенного напряжения для j-й точки наблюдения берется равным U(c)j, В. Если радиостанций несколько для каждой j-й точки измерения вычисляются значения наведенных напряжений Uij, В по формулам:
k-1 N
Ukj = U(c)j [сумма 10((Vij-Vkj)/10) + сумма 10((Vij-Vkj)/10)](-1/2)
i=1 i=k+1
Uij = Ukj 10((Vij-Vkj)/20), i=1, 2, ... k-1, k+1, ... N, (4.2)
где k - номер радиостанции, уровень сигнала которой превышает все остальные.
4.2.4. При проведении измерений прибором ВЗ-48А по схеме, показанной на рис.4, необходимо обеспечить надежный электрический контакт между соединительной шиной, корпусом пробника прибора и металлическим листом МЛ.
4.2.5. Результаты измерений и расчетов заносятся в протокол.
|
Главный государственный санитарный врач |
Г.Г.Онищенко |
Откройте актуальную версию документа прямо сейчас
Если вы являетесь пользователем интернет-версии системы ГАРАНТ, вы можете открыть этот документ прямо сейчас или запросить по Горячей линии в системе.
Методические указания МУК 4.3.678-97 "Определение уровней напряжений, наведенных электромагнитными полями на проводящие элементы зданий и сооружений в зоне действия мощных источников радиоизлучений" (утв. Главным государственным санитарным врачом РФ 6 ноября 1997 г.)
Текст Методических указаний приводится по официальному изданию Минздрава РФ (Москва, 1998 г.)