Методические указания МУК 4.3.3831-22 "Расчетное прогнозирование уровней электромагнитных полей на судах" (утв. Федеральной службой по надзору в сфере защиты прав потребителей и благополучия человека 26 декабря 2022 г.)

Методические указания МУК 4.3.3831-22
"Расчетное прогнозирование уровней электромагнитных полей на судах"
(утв. Федеральной службой по надзору в сфере защиты прав потребителей и благополучия человека 26 декабря 2022 г.)

ББК 28.071+51.22

Р24

ISBN 978-5-7508-2041-2

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

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

1.1. Настоящие методические указания (далее - МУК) устанавливают порядок расчета уровней магнитных полей (далее - МП) частотой 50 Гц, создаваемых электрооборудованием в судовых помещениях, и электромагнитных полей (далее - ЭМП) радиочастотного диапазона (далее - РЧ), создаваемых антеннами средств радиосвязи и навигационными радиолокационными станциями (далее - НРЛС) на открытых палубах судов различного назначения.

1.2. МУК не распространяются на электрические поля (далее - ЭП) частотой 50 Гц, создаваемые электрооборудованием, на носимую и аварийную аппаратуру.

1.3. МУК применяются при:

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

- выполнении расчетов уровней МП и ЭМП на стадии проектирования судов, при модернизации электрооборудования и радиоэлектронных средств (РЭС).

1.4. МУК предназначены для органов и организаций Федеральной службы по надзору в сфере защиты прав потребителей и благополучия человека, а также могут быть использованы проектными организациями и институтами, занимающимися проектированием и модернизацией судов, судовладельцами.

1.5. МУК носят рекомендательный характер.

II. Общие положения

2.1. Расчет уровней МП частотой 50 Гц и ЭМП РЧ проводится для определения соответствия параметров ЭМП гигиеническим нормативам ¹, установленным для плавательных средств и морских сооружений, с целью разработки конструктивных и организационных мероприятий по защите экипажа от неблагоприятного воздействия физических факторов.

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

¹_{Таблица 5.74 главы V СанПиН 1.2.3685-21 "Гигиенические нормативы и требования к обеспечению безопасности и (или) безвредности для человека факторов среды обитания", утвержденных постановлением Главного государственного санитарного врача Российской Федерации от 28.01.2021 N 2 (зарегистрировано Минюстом России 29.01.2021, регистрационный N 62296).}

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

2.2. Расчетное прогнозирование проводится:

- при проектировании судна (является составной частью процесса проектирования судна);

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

- при реконструкции и изменении архитектуры надстроек судна.

2.3. Расчет уровней напряженности ЭМП, создаваемых средствами радиосвязи в средневолновом (далее - СВ), коротковолновом (далее - КВ) ультракоротковолновом (далее - УКВ) диапазонах, выполняется методом граничных элементов с учетом архитектурно-конструктивных особенностей судна, технических характеристик средств радиосвязи и вторичного излучения от ближайших к антенне металлических конструкций.

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

2.5. Определение уровней МП частотой 50 Гц выполняется на высотах 0,5, 1,0, 1,7 м от поверхности палубы, от источника ЭМП начиная от 0,2 м и далее с шагом 0,5 м. Расчеты проводятся для различных расстояний от источника до установления предельно допустимого уровня (далее - ПДУ).

III. Расчет уровней напряженности магнитного поля частотой 50 Гц, создаваемого судовым электрооборудованием

3.1. Основные положения методики расчета МП частотой 50 Гц.

3.1.1. К основным источникам МП частотой 50 Гц на судах относятся генераторы, распределительные щиты, электродвигатели, трансформаторы, статические преобразователи.

3.1.2. Исходными данными для расчета уровней магнитного поля частотой 50 Гц судового электрооборудования являются:

- дипольные магнитные моменты электрооборудования частотой 50 Гц;

- номинальная мощность (для генераторов и электродвигателей);

- число пар полюсов (для генераторов и электродвигателей);

- полная мощность (для трансформаторов);

- номинальный ток (для распределительных щитов);

- номинальный выпрямленный ток (для статических преобразователей);

- координаты геометрических центров источников.

3.1.3. При расчете магнитного поля частотой 50 Гц, каждый источник аппроксимируется точечным дипольным источником с определенным дипольным магнитным моментом, расположенным в геометрическом центре источника.

3.1.4. Если источник составлен из нескольких модулей, то при расчете каждый модуль рекомендуется представить как отдельный дипольный источник со своим дипольным магнитным моментом, расположенным в геометрическом центре модуля.

3.2. Определение дипольных магнитных моментов электрооборудования производится по номинальной мощности или номинальному току электрооборудования с использованием формул, представленных в табл. 1.

Таблица 1

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

Тип электрооборудования	Формула для расчета
Генераторы и электродвигатели с одной парой полюсов	M = P
Генераторы и электродвигатели с более чем одной парой полюсов	М = 0,36 Р sin /p
Трансформаторы	М = 0,36 Q
Распределительные щиты	М = 0,33 I
Статические преобразователи	М = 0,26 I

где: М - дипольный магнитный момент в А·м ²;

Р - номинальная мощность в кВт;

р - число пар полюсов;

Q - полная мощность в кВ·А;

I - номинальный ток в А.

3.3. Расчет магнитного поля дипольного источника.

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

Расстояние между центром источника и расчетной точкой вычисляется по формуле (1):

,

(1)

где: r - расстояние между центром источника и расчетной точкой, м;

x ₀, y ₀, z ₀ - координаты точки расположения источника, м;

x ₁, y ₁, z ₁ - координаты расчетной точки, в которой рассчитывается поле, м.

Напряженность магнитного поля одного точечного дипольного источника рассчитывается по формулам (2 - 4):

;

(2)

;

(3)

,

(4)

где: H _x, H _y, H _z - проекции вектора напряженности магнитного поля на оси х, у, z декартовой системы координат, А/м;

M _x, M _y, M _z - проекции вектора дипольного магнитного момента на оси х, у, z А/м ²;

x ₀, y ₀, z ₀ - координаты точки расположения источника, м;

x ₁, y ₁, z ₁ - координаты расчетной точки, в которой рассчитывается поле, м.

Если известна абсолютная величина М дипольного магнитного момента источника, но неизвестно направление вектора дипольного магнитного момента, то необходимо выполнить три варианта расчета для дипольных магнитных моментов с проекциями, формула (5):

;

;

,

(5)

где: М - абсолютная величина дипольного магнитного момента, А·м ²;

M _x, M _y, M _z - проекции вектора дипольного магнитного момента на оси x, y, z, А·м ².

Для оценки уровня магнитного поля рекомендуется использовать максимальное из трех полученных значений.

Модуль напряженности магнитного поля источника вычисляется по формуле (6):

,

(6)

где: Н - модуль напряженности магнитного поля источника, А/м;

n - число точечных дипольных источников, на которые заменялся при расчете реальный источник;

Н _ix, Н _iy, Н _iz - проекции вектора напряженности магнитного поля i-того точечного дипольного источника на оси х, у, z, вычисленные по формулам (2 - 4), А/м.

IV. Расчет уровней напряженности электрических и магнитных полей судовых средств радиосвязи СВ, КВ, и УКВ диапазонов

4.1. Судовые средства радиосвязи работают в СВ 0,03 - 3 МГц, КВ 3 - 30 МГц и УКВ 30 - 300 МГц диапазонах. При эксплуатации средств радиосвязи используются различные типы антенн - штыревые, дипольные, коллинеарные дипольные, Т-образные, Г-образные и другие.

4.2. Исходными данными для расчета уровней ЭП и МП судовых средств радиосвязи являются:

- рабочие диапазоны частот;

- максимальные мощности;

- коэффициенты полезного действия антенн;

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

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

- данные о расположении палуб и надстроек судна (координаты и размеры).

4.3. Расчеты уровней ЭП и МП выполняются в расчетных точках на различных расстояниях от антенны до получения значения ЭМП, не превышающего ПДУ, что определяет границы зон с повышенными уровнями ЭП и МП.

4.4. Расчеты уровней ЭП и МП выполняются отдельно для различных частот f. Если антенна предназначена для излучения в рабочем диапазоне частот, формула (7):

,

(7)

где: f _мин - минимальная частота диапазона, Гц;

f _макс - максимальная частота диапазона, Гц.

Если отношение максимальной частоты диапазона к минимальной частоте диапазона удовлетворяет соотношению, формула (8):

,

(8)

расчет может выполняться для одной произвольной частоты диапазона.

Если условие формулы (8) не выполняется, необходимо выбрать набор частот f ₁, f ₂, f _N где f ₁ = f _min; f _N = f _max, так, чтобы каждая следующая частота не более чем в 1,4 раза превосходила предыдущую, и выполнить расчет для каждой из этих частот.

Расчет напряженностей ЭП и МП судовых средств радиосвязи методом граничных элементов

4.5. Расчет уровней напряженности ЭП и МП, создаваемых антеннами средств радиосвязи, выполняется на основе численного решения интегрального уравнения электрического поля (ИУЭП), формула (9):

,

(9)

где: i - комплексная единица;

- 2f - круговая частота (где f - частота излучения антенны, Гц), рад/с;

 ₀ - 1,257·10 ^-6 Гн/м - магнитная постоянная;

 ₀ - 8,854·10 ^-12 Ф/м - электрическая постоянная;

S - совокупность поверхностей всех проводников;

r - точка наблюдения на поверхности S;

r' - точка интегрирования на поверхности S;

I(r') - вектор поверхностной плотности электрического тока в точке, А/м.

Функция Грина, формула (10):

,

(10)

где: - волновое число, |r-r'| - расстояние между точками r и r';

grad _s - оператор поверхностного градиента;

div _s - оператор поверхностной дивергенции;

- вектор напряженности электрического поля источника в точке r, В/м.

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

4.6. Уравнение (10) сводится к системе линейных алгебраических уравнений (далее - СЛАУ) с помощью метода граничных элементов. Для этого вектор поверхностной плотности электрического тока аппроксимируется суммой, формула (11):

,

(11)

где: - вектор поверхностной плотности электрического тока в точке r, А/м;

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

- векторные базисные функции, м ^-1;

N - общее число базисных функций.

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

 _мин/10,

(12)

где:  _мин - минимальная длина волны, м.

Минимальная длина волны вычисляется по формуле (13):

 _мин = с/f _макс,

(13)

где: с 3·10 ⁸ м/с - скорость света.

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

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

,

(14)

где: и - треугольные элементы, имеющие общую сторону;

+ и - площади треугольных элементов и , м ²;

и - радиус-векторы из вершин треугольных элементов и, противолежащих общей стороне и , в точку r (рис. 1), м.

Рис. 1. К определению базисных функций на треугольных элементах

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

(15)

где: и - цилиндрические элементы, соприкасающиеся торцами;

и - радиусы цилиндрических элементов и , м;

и - длины цилиндрических элементов и , м;

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

Рис. 2. К определению базисных функций на цилиндрических элементах

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

4.8. Для всех пар индексов m к n (1 m N, 1 n N) вычисляются коэффициенты СЛАУ а _mn по формуле (16):

,

(16)

где: - скалярное произведение векторов и .

Все коэффициенты b _m правой части СЛАУ задаются равными 0, за исключением коэффициента b _m с индексом m = n _пит, который полагается равным 1.

4.9. Решается СЛАУ, формула (17):

.

(17)

Результатом решения системы (17) являются значения коэффициентов I _n.

4.10. Для каждой расчетной точки r _рас напряженности электрического и магнитного полей вычисляются по формулам (18):

,

(18)

где: E(r _рас) - напряженность электрического поля в расчетной точке r _рас, В/м;

H(r _рас) - напряженность электрического поля в расчетной точке r _рас, А/м;

k _мощ - безразмерный коэффициент, вычисляемый по формуле (19):

,

(19)

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

- комплексный коэффициент с индексом n _пит, найденный в результате решения СЛАУ (17), А;

Re - вещественная часть комплексного коэффициента , А.

Максимальная мощность, излучаемая антенной, связана с максимальной мощностью P _макс, подводимой к антенне, соотношением (20):

.

(20)

где: - коэффициент полезного действия антенны.

Если коэффициент полезного действия антенны неизвестен, то рекомендуется производить расчет для = 1, при этом P _изл = P _макс.

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

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

5.1. Принципы расчета плотности потока энергии (далее - ППЭ) навигационных радиолокационных станций.

5.1.1. При оценке навигационных радиолокационных станций рассчитывается средняя за период следования импульсов ППЭ, которая определяется исходя из предполагаемого распределения ЭМП в апертуре прямоугольной антенны НРЛС. Распределение ЭМП в апертуре выбирается на основе характеристик диаграммы направленности НРЛС.

5.1.2. Расчет уровней ППЭ производится только для расчетных точек, находящихся в прямой видимости относительно НРЛС (то есть на прямой между антенной и расчетной точкой отсутствуют экранирующие или отражающие конструкции).

5.1.3. Исходными данными для расчета уровней ППЭ НРЛС являются:

- импульсная мощность НРЛС;

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

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

- коэффициент усиления антенны;

- длина волны излучения;

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

- размер апертуры антенны в вертикальной плоскости;

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

- координаты центра апертуры антенны.

5.2. Расчет ППЭ навигационных радиолокационных станций.

Средняя излучаемая мощность НРЛС определяется по формуле (21):

,

(21)

где: P _ср - средняя излучаемая мощность, Вт;

P _имп - импульсная мощность НРЛС, Вт;

- длительность импульсов, с;

f _имп - частота следования импульсов, Гц.

Расчет ППЭ рекомендуется проводить для условий эксплуатации навигационной радиолокационной станции в режиме максимальной дальности действия, что соответствует работе НРЛС на максимальной шкале дальности.

Когда антенна направлена в сторону расчетной точки ППЭ (расчетная точка и ось излучения антенны находятся в одной вертикальной плоскости), взаимное расположение антенны и расчетной точки описываются расстоянием R от центра апертуры антенны до расчетной точки и углом в вертикальной плоскости между осью излучения антенны и направлением на расчетную точку, которые вычисляются по формулам (22 - 23):

;

(22)

,

(23)

где: R - расстояние от центра апертуры антенны до расчетной точки, м;

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

x ₀, y ₀, z ₀ - координаты центра апертуры антенны, м;

x ₁, y ₁, z ₁ - координаты расчетной точки, м.

ППЭ в расчетной точке определяется по формуле (24):

,

(24)

где: G - коэффициент усиления антенны (безразмерная величина);

B(d) - функция, описывающая изменение ППЭ на оси излучения антенны, (безразмерная величина);

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

d ₁, d ₂, u - безразмерные параметры, которые рассчитываются по формулам (25 - 27):

;

(25)

;

(26)

,

(27)

где: R - расстояние от центра апертуры антенны до расчетной точки, м;

- длина волны излучения, м;

L - размер апертуры антенны в горизонтальной плоскости, м;

H - размер апертуры антенны в вертикальной плоскости, м.

Функции B(d) и F(u,d) задаются таблично и зависят от распределения ЭМП в апертуре антенны. Вид распределения ЭМП в апертуре антенны рекомендуется выбирать в зависимости от уровня первого бокового лепестка диаграммы направленности антенны, параметры которого указаны в технической документации к НРЛС. Если уровень первого бокового лепестка более - 20 дБ, то рекомендуется выбирать равномерное распределение ЭМП по апертуре. Если уровень первого бокового лепестка от 20 до 30 дБ, рекомендуется выбирать косинусоидальное распределение ЭМП по апертуре. Если уровень первого бокового лепестка менее 30 дБ, рекомендуется выбирать распределение "косинус-квадрат". В приложении 1 к настоящим МУК приведены значения функций B(d) (табл. 1.1) и F(u,d) (табл. 1.2 - 1.4) для равномерного распределения ЭМП по апертуре, косинусоидального распределения и распределения "косинус-квадрат". При значениях d > 1 рекомендуется использовать значения функций B(d) и F(u,d) для d = 1.

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

6.1. Методика расчета ППЭ судовых станций спутниковой связи (далее - ССС).

6.1.1. Расчет уровней ППЭ ССС проводится для максимального излучения антенны (только для положения антенны, при котором ось излучения антенны направлена в расчетную точку).

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

6.1.3. Расчет уровней ППЭ производится только для расчетных точек, находящихся в прямой видимости относительно ССС (то есть на прямой между антенной и расчетной точкой отсутствуют экранирующие или отражающие конструкции).

6.1.4. Исходными данными для расчета уровней ППЭ являются:

- максимальная мощность ССС;

- коэффициент усиления антенны;

- длина волны излучения;

- радиус апертуры антенны;

- координаты центра апертуры антенны.

6.2. Расчет ППЭ ССС.

Расстояние между центром апертуры антенны и расчетной точкой вычисляется по формуле (28):

,

(28)

где: R - расстояние от центра апертуры антенны до расчетной точки, м;

x ₀, y ₀, z ₀ - координаты центра антенны, м;

x ₁, y ₁, z ₁ - координаты расчетной точки, м.

ППЭ в расчетной точке определяется по формуле (29):

,

(29)

где: R - расстояние между центром апертуры антенны и расчетной точкой, м;

Р - максимальная мощность ССС, Вт;

G - коэффициент усиления антенны (безразмерная величина);

B(d) - функция, описывающая изменение ППЭ на оси излучения антенны (безразмерная величина);

d - безразмерный параметр, который рассчитывается по формуле (30):

,

(30)

где: R - расстояние между центром апертуры антенны и расчетной точкой, м;

- длина волны излучения, м;

а - радиус апертуры антенны, м.

В приложении 1 к настоящим МУК (таблица 1.5) приведены значения функций B(d). Они получены в предположении, что ЭМП в апертуре пропорционально функции распределения, которая задается формулой (31):

,

(31)

где: f(r) - функция распределения ЭМП по апертуре (безразмерная величина);

r - расстояние от точки апертуры до центра апертуры, м;

а - радиус апертуры антенны, м.

При значениях d > 1 рекомендуется использовать значение функции B(d) для d = 1.

VII. Мероприятия по защите членов экипажа от облучения ЭМП

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

7.2. Обеспечение защиты включает в себя конструктивные, организационные мероприятия.

7.3. Конструктивные мероприятия:

- рациональное размещение излучающего оборудования;

- применение защитных экранов.

7.4. Организационные мероприятия:

- применение леерных ограждений антенн;

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

Библиографические ссылки

1. Федеральный закон от 30.03.1999 N 52-ФЗ "О санитарно-эпидемиологическом благополучии населения".

2. Федеральный закон от 28.12.2013 N 412-ФЗ "Об аккредитации в национальной системе аккредитации".

3. Федеральный закон от 26.06.2008 N 102-ФЗ "Об обеспечении единства измерений".

4. СанПиН 1.2.3685-21 "Гигиенические нормативы и требования к обеспечению безопасности и (или) безвредности для человека факторов среды обитания".

5. РД 5Р.8903-96 "Аппаратура радиосвязи и радиолокации. Порядок выполнения работ по защите личного состава судов от облучения ЭМП".

6. РД 5Р.8713-93 "Аппаратура радиосвязи и радиолокации. Методы оценки электромагнитных полей и средства защиты личного состава судов от облучения".

7. Никитина В.Н., Ляшко Г.Г., Калинина Н.И. Анализ современного состояния электромагнитной безопасности экипажа на судах ледокольного флота//Гигиена и санитария. 2018. N 12. С. 1210 - 1214.

8. Ерохин В.Л. Арктическое судостроение и морская техника России: состояние, потребности, перспективы развития И Маркетинг и логистика. 2019. N 1 (21). С. 12 - 31.

9. Никитина В.Н., Дубровская Е.Н., Калинина Н.И., Ляшко Г.Г. Расчетное прогнозирование уровней электромагнитных полей радиочастотного диапазона на открытых палубах современных судов рыбопромыслового флота // Материалы I Национального конгресса с международным участием по экологии человека, гигиене и медицине окружающей среды "СЫСИНСКИЕ ЧТЕНИЯ - 2020" Сборник тезисов, Москва, 19 - 20 ноября 2020 года. Москва: Федеральное государственное бюджетное учреждение "Центр стратегического планирования и управления медико-биологическими рисками здоровью", 2020. С. 259 - 262.

10. Разлетова А.Б., Никитина В.Н., Ляшко Г.Г., Калинина Н.И., Дубровская Е.Н., Котенко П.К. Гигиеническая оценка особенностей трудовой деятельности экипажа в условиях Арктики и нормативное обеспечение охраны здоровья моряков: проспективное исследование // Морская медицина. 2022. Т. 8. N 3. С. 88 - 99.

11. Lange S., Sabath F., Baganz N. Radiation hazard assessment of shipbome RF communication // International Symposium on Electromagnetic Compatibility-EMC EUROPE. - IEEE, 2012.

12. Sanchez L.F. et al. HF broadband antenna design for shipboard communications: Simulation and measurements // Measurement. 2016. Vol. 89. P. 13 - 20.

13. Sanchez L.F. et al. Radiation hazards to personnel from non-ionizing fields of broadband HF systems onboard a vessel: Measurement and simulation // Measurement. 2018. Vol. 115. P. 223 - 232.

14. Matkovic A., Sarolic A. Near Electric Field Around the Shipboard Navigational Radar Patch Array Antenna: Comparison to Human Exposure Limits // 2019 International Conference on Software, Telecommunications and Computer Networks (SoftCOM). - IEEE, 2019.

15. Spandonidis C.C., Arvaniti K. Numerical modeling of the exposure on radiofrequency radiation of marine mammal observers during their shift: a case study // Int. J. Saf. Secur. Eng. 2020. Vol. 10. P. 753 - 758.