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

884612727 На правах рукописи

Довбыш Владимир Николаевич

ТЕХНОЛОГИЯ РЕГИОНАЛЬНОГО КОНТРОЛЯ ПРИРОДНОЙ СРЕДЫ ПО ФАКТОРУ ЭЛЕКТРОМАГНИТНОГО ИЗЛУЧЕНИЯ ОБЪЕКТОВ ЭНЕРГЕТИЧЕСКИХ СИСТЕМ

Специальность 05.11.13 Приборы и методы контроля природной среды, веществ, материалов и изделий

Автореферат диссертации на соискание ученой степени доктора технических наук

Самара-2010

1 8 НОЯ 2010

004612727

Работа выполнена в Государственном образовательном учреждении высшего профессионального образования «Поволжский Государственный Университет телекоммуникаций и информатики» (ГОУ ВПО ПГУТИ)

Научный консультант

доктор технических наук, профессор, Сподобаев Юрий Михайлович.

Официальные оппоненты:

доктор технических наук, профессор, Седельников Юрий Евгеньевич,

профессор кафедры Радиоэлектронных и телекоммуникационных систем Казанского государственного технического университета им. А.Н. Туполева;

доктор технических наук, профессор, Корнилов Владимир Юрьевич,

профессор кафедры Электропривода и автоматизации промышленных установок и технологических компонентов Казанского государственного энергетического университета;

доктор технических наук, профессор Стеблев Юрий Иванович,

заведующий кафедрой Автоматизации производственных процессов Самарского государственного технического университета.

Ведущая организация

Институт проблем управления сложными системами Российской академии наук (г. Самара)

Защита диссертации состоится 26 ноября 2010 г. в 14 часов на заседании диссертационного Совета Д 212.079.04 в ГОУ ВПО «Казанский государственный технический университет им. А.Н. Туполева» по адресу: 420111, г. Казань, ул. К. Маркса, д. 10.

С диссертацией можно ознакомиться в библиотеке Казанского государственного технического университета им. А.Н. Туполева, с авторефератом - на сайте ВАК РФ: http://vak.ed.gov.ru/ru/announccmcnts_l/technikal_scienccs/

Автореферат разослан «_»

2010 г.

Ученый секретарь диссертационного совета

Седов С.С.

ОБЩАЯ ХАРАКТЕРИСТИКА РАБОТЫ

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

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

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

Человеческое общество в процессе своей деятельности воздействует на различные компоненты природной среды: атмосферу, гидросферу и литосферу. В настоящее время такие воздействия приобретают глобальный характер, затрагивая все континенты нашей планеты. Анализ значительного количества публикаций геоэкологической тематики отечественного и зарубежного происхождения, позволяет отметить то обстоятельство, что среди всего многообразия видов антропогенного загрязнения природной среды, можно выделить один специфический вид энергетического загрязнения, а именно, электромагнитный - вид, влияние которого можно одновременно отнести практически ко всем компонентам природной среды одновременно. Отмеченные вопросы относятся к специфической области знания — «электромагнитной экологии» и традиционно проблемы с ними связанные решаются при помощи электромагнитного мониторинга природной среды, который состоит:

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

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

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

- разработки мероприятий и рекомендаций по защите от ЭМП и нормализации электромагнитной обстановки.

Современному этапу развития человеческого общества, соответствует региональный характер развития всех инфраструктур, в том числе и связанных с потреблением и излучением электромагнитной энергии. Действительно, все энергоемкие объекты промышленного производства и культурной деятельности человека, в значительной степени, сосредоточены в локальных областях, поэто-

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

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

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

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

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

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

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

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

Особая народнохозяйственная значимость проблем контроля среды пребывания человека по электромагнитному фактору подчеркнута Решением Коллегии Федеральной службы Роспотребнадзора от 27 ноября 2009 г., в котором определены задачи и приоритетные направления по детальному изучению электромагнитной обстановки в населенных пунктах, инвентаризации источников электромагнитных излучений, предусматривающие планомерное снижение электромагнитной нагрузки на население и создание баз данных.

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

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

Так методологические основы технологий регионально контроля природной среды развиты в работах Булгакова Н.Г., Левича А.П. и Максимова В.Н.

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

ческих средств стали работы Шередько Е.Ю., Сподобаева Ю.М., Кубанова В.П., Маслова О.Н., Бузова АЛ., Романова В.А., Казанского JT.C. В этих работах предложены и обоснованы подходы к расчетному прогнозированию электромагнитной остановки вблизи широкого класса излучающих технических средств и их комплексов, сформулированы подходы к системному электромагнитному мониторингу. Основные результаты этих работ, подтвержденные многочисленными экспериментальными исследованиями, нашли отражение в нормативно-методических документах, утвержденных государственными органами санитарно-эпидемиологического надзора. Очевидно, что общие подходы и ряд математических методов, примененные названными авторами для электромагнитного мониторинга комплексов излучающих технических средств, вполне применимы и для решения базовых задач комплексного анализа электромагнитной обстановки вблизи энергетического оборудования.

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

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

- группу сосредоточенных (локальных) технических средств — силовые трансформаторные установки, тяговые подстанции электротранспорта, распределительные пункты системы энергоснабжения.

Названные технические средства, в основном, являются источниками либо статического ЭМП, либо ЭМП ПЧ. Для случая ЭМП ПЧ выполняется условие квазистационарности, поэтому при расчетном прогнозировании электромагнитной обстановки задача может быть сформулирована аналогично статическому случаю. Иными словами, моделирование для целей расчета электрического и магнитного полей может производиться раздельно. Вопросы расчетов стационарных ЭМП электроустановок подробно рассмотрены в работах Колечицкого Е.С., Меликова H.A., Филиппова A.A., Тозони О.В., Демирчана К.С., Abou-Seada M.S. и Nasser Е., Silvester Р. и Chari M., a так же других авторов.

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

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

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

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

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

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

Многочисленность влияющих факторов и сложность структуры первичных токов диктует целесообразность применения для расчета численных электродинамических методов, ориентированных на непосредственное решение уравнений Максвелла таких, как метод конечных элементов (МКЭ) (Демирчан К.С., Abou-Seada M.S. и Nasser Е., Silvester Р. и Chari М.).

Основной специфической особенностью распределительных объектов сетевой иерархии является то, что их конструкции специфичны и не унифицирова-

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

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

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

Воздушные ЛЭП и силовые трансформаторы, помимо ЭМП ПЧ, являются источниками высокочастотных полей. Причинами данных явлений являются, соответственно, коронирование проводов и частичные разряды в обмотках. Оценка уровней ЭМП, создаваемых коронными разрядами, возможна в силу простоты конфигурации первичных токов, при помощи приближенного подхода, развитого в работах Казакова В.Н., а также Кима К.С., Лелевкина В.М., Токарева A.B., ЮдановаВЛ.

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

Как отмечалось выше, данные об электромагнитной обстановке в масштабах региона представляют собой массивы данных, значительного объема. При этом весьма существенным обстоятельством, отличающим такие данные от результатов, получаемых в иных формах экологического мониторинга, является «генетическая» привязанность к географическим координатам. Действительно, области «сильного поля» энергетического оборудования, очевидно оказываются локализованными вблизи мест расположения технических средств - источников, а результирующая электромагнитная обстановка образует сложную пространственную картину, привязанную к рельефу местности.

Иными словами, неотъемлемой частью технологии регионального контроля природной среды по фактору ЭМП, является частная технология визуализации и графической обработки геоэкологической информации. Применение геоинформационных технологий при региональном экологическом контроле в последние годы стало стандартным решением. Вопросам, связанным с геоэкологическим картографированием, повещены работы Берлянта A.M., Полякова М.М., Прогу-ловой Т.Б., Мясоедова Б.В., а также других авторов. Общие подходы и примененные способы представления векторных и матричных данных на электронных картах могут быть с успехом применены для реализации целей, поставленных в настоящей диссертации.

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

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

В диссертационной работе решаются следующие задачи:

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

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

3. Исследование электромагнитных полей реальных объектов, функционирующих в Самарской области.

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

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

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

7. Построение фрагментов геоинформационной системы электромагнитной безопасности энергетической инфраструктуры Самарской области.

Методы исследования. В работе использованы методы математического моделирования, аналитический аппарат вычислительной электродинамики, методы решения дифференциальных уравнений в частных производных, аналитической геометрии, численные методы, математический аппарат теории ИУ. Визуализация результатов произведена с использованием ГИС-технологий.

Все результаты получены с использованием вычислительных алгоритмов, реализованных на ЭВМ на языке Fortran и в среде Matlab. Фрагменты электронных карт получены при помощи программного пакета ГИС-Карта-2000.

Научная новизна исследований заключается в следующем:

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

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

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

4. Получены новые результаты численного анализа электромагнитной обстановки вблизи реальных объектов региональной энергетической системы.

5. Впервые, на основе комплекса разработанных подходов, представлены результаты регионального контроля природной среды по фактору электромагнитного излучения.

Практическая значимость результатов работы заключается в том, что:

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

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

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

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

Реализация результатов работы.

Созданная в диссертационной работе технология контроля состояния природной среды, методики и алгоритмы анализа электромагнитной обстановки использованы при подготовке отчетных материалов при выполнении работ «Организация и проведение мониторинга электромагнитного излучения по Самарской области» в рамках государственного контракта, заключенного между ПГУТИ, Фондом социально-экологической реабилитации Самарской области и Министерством природных ресурсов и охраны окружающей среды Самарской области, при поддержке управления Росприроднадзора по Самарской области (ГК №1/11 от 13.02.2006 г.)

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

Научно-прикладные результаты работы внедрены при выполнении работ по теме «Разработка методов прогнозирования и визуализации электромагнитной обстановки элементов энергосистемы городского округа Новокуйбышевска», 10

выполненных по договору 1/8 от 10.04.2008 г. между ПГУТИ и Целевым бюджетным экологическим фондом г. Новокуйбышевска.

Принципы и методики моделирования электромагнитных полей, создаваемых силовыми трансформаторами, а так же некоторые результаты расчетов использованы при проектировании, размещении и эксплуатации трансформаторных подстанций выпускаемых ЗАО «Электрощит - Самара».

Методика и алгоритмы моделирования ЭМП трансформаторных подстанций и ЛЭП использованы при проектировании ряда энергетических объектов, строительство которых проходит в г. Самара и Самарской области (заказчик -ЗАО «Самарский Электропроект», 2005-2008 г. г.).

Разработанные автором методики комплексного анализа электромагнитной обстановки в регионе, а так же результаты исследования ЭМП реальных объектов внедрены в учебный процесс и практику дипломного проектирования в ПГУТИ на кафедре «Электродинамики и антенн», а так же в Самарском региональном телекоммуникационном трейнинг-центре.

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

Апробация результатов работы н публикации.

Основные результаты по теме диссертационного исследования докладывались на IX и X международных конгрессах «Актуальные проблемы экологии человека» (Самара, 2004 и 2005 г., соответственно), X, XI, XII, XIII и XIV Всероссийских научно-технических конференциях ПГАТИ (Самара, 2003, 2004, 2005, 2006 и 2007 г., соответственно), VI и VIII Международных научно-технических конференциях «Проблемы техники и технологии телекоммуникаций» (Самара, 2005 г. и Уфа 2007 г.), XXXIII Всероссийском семинаре «Актуальные вопросы охраны окружающей среды. Программное обеспечение экологов» (Москва, 2005 г.), IV и VI Международных симпозиумах по электромагнитной совместимости и электромагнитной экологии (Санкт-Петербург, 2005 и 2009 г.).

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

Публикации

По тематике диссертационных исследований автором (лично и в соавторстве) опубликовано 49 печатных работ.

Основные научные и прикладные результаты опубликованы в 15 статьях в периодических научных изданиях (из них 11 в журналах, входящих в список ВАК) и в 31 публикациях в форме материалов российских и международных конференциях и семинарах.

Автором опубликованы две монографии (лично и в соавторстве) и один нормативно-методический документ регионального значения.

На защиту выносятся:

1. Технология регионального контроля природной среды по фактору электромагнитного излучения объектов энергетических систем, как методологически единый комплекс частных методик и алгоритмов, предназначенных для ЭМП, созда-

II

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

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

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

4. Новые результаты анализа электромагнитной обстановки вблизи ряда реальных объектов системы энергоснабжения.

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

Объем и структура работы.

Диссертационная работа состоит из введения, пяти разделов, заключения, списка литературы и приложений. Работа содержит 305 страниц, включая 112 рисунков и 5 приложений. Список литературы содержит 200 наименований.

СОДЕРЖАНИЕ РАБОТЫ

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

В первом разделе «Разработка электродинамических моделей распределительных элементов региональной энергетической инфраструктуры» проведен общий анализ объекта исследования, предложена классификация источников по критериям пространственной локализации и характеристикам излучаемого ЭМП. Также выявлены особенности соответствующих электродинамических задач, обоснованы общие подходы к электродинамическому моделированию источников высоковольтных ЛЭП. Приняты ряд допущений для моделирования и указаны границы применимости моделей.

В разделе показано, что основными техническими средствами входящими в состав региональной энергосистемы - «поставщиками» электромагнитной энергии ПЧ, существенно влияющей на общую электромагнитную обстановку в регионе, являются высоковольтные ЛЭП, силовые трансформаторные установки, сети питания наземного и подземного электротранспорта. При этом контролируемыми параметрами ЭМП, создаваемых элементами энергетических систем, при оценке их воздействия на окружающую среду являются:

- напряженность электрического поля Е в В/м;

- напряженность магнитного поля, Я в А/и.

В настоящее время основными документами, которые определяют предельно допустимые уровни поля от воздушных ЛЭП и других объектов, являются СН № 2971-84, СанПиН 2.1.2.1002-00 СанПиН 2.2.4.1191-03 и ГН 2.1.8/2.2.4.2262-07. В принципе этих документов достаточно, чтобы оценивать электромагнитную обстановку, создаваемую элементами энергетических систем. 12

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

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

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

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

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

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

Рассмотрим модель прямолинейного участка цепи электроснабжения с точки зрения вычисления электрического поля. Поскольку напряжение в сети не зависит от нагрузки, электрическое поле также оказывается независимым от потребляемого тока. ЛЭП, конфигурация проводов которой соответствует типовой опоре У-35, размещенная в декартовой системе координат, показана на рис.1.

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

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

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

(1)

где Иф-класс напряжения ЛЭП, А = 120° - фазовый сдвиг, j-мнимая единица, [¡/ - угол, образованный векторами Et и Е\, косинус которого легко

R. R4Й,2

определяется из рис.1 следующим образом: 2 cos у = — +----.

Я, R, R,R,

<z

-Ч:

Рис.1. К расчету электрического поля высоковольтной ЛЭП

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

Напряженность магнитного поля воздушной ЛЭП определяется аналогичным

образом:

/•е

Л.-1)Д

1

1

V (Й;)2 •(«;)

- eos ^,

(2)

— От .....1* .... 2» • — •4»

/ \ "Í--9т

\

/ N

■1

V/

На рис.2 приведены результаты расчета электрического поля ЛЭП 35 кВ, размещенной на опоре У-35-1. Расчет проведен в направлении перпендикулярном коридору прохождения ЛЭП.

Слияние кривых изменения поля, построенных для различных высот над поверхностью Земли, в одну линию, на значительном удалении от ЛЭП обусловлено преобладанием линейного размера х, по сравнению с размером г в тре-

Рис.2. Результат расчета электрического поля ЛЭП 35 кВ

' угольнике расстоянии,

что, в свою очередь приводит к тому, что величина 2 со?, у/ оказывается чрезвычайно малой и не оказывает существенного влияния на вид зависимости поля от расстояния.

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

Основным источников высокочастотного излучения ЛЭП является электрическая корона, возникающая на токонесущих проводниках вследствие резко выраженной неоднородности электрического поля вблизи последних. Часто при коронном разряде проводники окружены характерным свечением, также получившим название короны, или коронирующего слоя.

Распространенной формулой для расчета потерь на корону на переменном напряжении является эмпирическая формула Пика для одиночного провода:

Р - ~~—(f + 25)

ику

•1(Г5, Вт/км,

(3)

где 5 - относительная плотность воздуха, / - частота, Гц, г0 - радиус одиночного провода, см, 5 - расстояние между проводами, см, {/к - напряжение возникновения короны, определяемое выражением:

ик~2\,2-б-1п — тгт2,кВ

(4)

где т1 - коэффициент гладкости провода (для идеально гладкого провода т1 = 1, для реального витого провода /я, = 0.85-0.92, согласно эмпирическим данным, приведенным в литературе), т2 - коэффициент погоды.

Используя формулу (4) получено выражение для напряженности электрического поля короны в стационарном режиме:

1 [У

£ = -

(5)

2лК1аа \ Лд

где иа - электропроводность атмосферы при данных погодных условиях, Е0 -электрическое сопротивление цилиндра, коаксиального проводнику, длиной /. При расчетах, как правило, рассматривают несколько характерных случаев: хорошая погода (т2 = 1), дождь, включая мокрый снег и изморось (т2 =2), сухой снег (т2 = 3 ),иней и наледь на проводах (т2 = 4).

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

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

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

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

Систему уравнений Максвелла для квазистационарного случая преобразуется в систему уравнений второго порядка:

—

■ Е&.Я-})

- ЩГ,**!)

\ Ztr.nl)

----------- -4. :

- Г'"--

Рис.3. Распределение напряженности электрического поля коронного разряда токонесущего проводника ЛЭП 35 кВ для различных погодных условий

У2£ = ёгас1 —, -У2Я = г<* }. (6)

£0£

Данную систему дополняется условиями, накладываемыми на искомое решение на границе области тока.

Для электрического поля это условия следующего вида:

1 3

£Й0£[-£Й0£2 =—[Я0£1]-[Й0£2] = 0. (7)

е0 дБ

Условия для магнитного поля имеют вид:

Мп0Нх-МИ0Н2 =0, [п0Н1]-[п0Н2] = 10-^-, (8)

01

где Щ - вектор нормали к границе области тока, у - заряд, распределенный на границе, / -тангенциальная составляющая тока, на краевой линии элемента, о/ - элемент контура, нормального линиям тока, Г0 - орт линий тока.

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

В общем случае в обобщенных координатах уравнения Лапласа (6) могут быть представлены в виде (использование обобщенных координат целесообразно ввиду неопределенности геометрических особенностей некоторых задач, встречающихся в разделах 2 и 3, решаемых также при помощи данной реализации МКЭ):

(9)

дд ОГ]

где 1]) и /?(£, т]) - линейные функции или операторы, зависящие от типа выбранной системы координат; £,т] - обобщенные координаты; Ф( £,?/)- обобщенное решение уравнения (компоненты векторов в или Н^ -

известная функция, определяемая сторонними источниками.

Решение аппроксимируется линейной функцией, специфической для

каждого элемента, что можно представить как приближенное решение:

= . (ю)

п=1 т=1

где, так называемая, функция формы глобальной системы координат

имеет смысл базисной функции для разложения (10); Ф,, Ф2, Ф3 - узловые

значения искомой величины; Фт т^р)-принимают различные значе-

ния при различных индексах т; М- максимальный номер элемента.

Для определения значений коэффициентов разложения целесообразно использование какого-либо проекционного метода, например, метода Галеркина.

При аппроксимации решения выражением (10) невязка решения определяется выражением:

e = a{$,rj)

82Ф(4,г])

д2Ф(4,П)

(11)

а/72

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

j4">-*-dD = 0, (12)

D

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

Краевые условия для уравнения (11) в общем виде:

Ж. (6 Л) - 1)Ф2 (6 п) = 7'), (13)

где координаты на границе раздела двух элементов; и 8(^,1]) -

функции, определяющие ориентацию векторов поля относительно границы.

Численная реализация МКЭ требует совместного решения системы линейных уравнений вида (12), количество уравнений в которой равно количеству узлов в модели.

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

дуемой ЛЭП - вторая очередь набережной реки Волга. Класс напряжения ЛЭП - 110 кВ, предполагаемая номинальная загрузка - 100 MB А. Рассмотрен традиционный способ прокладки подземных ЛЭП с заземлением внешнего проводника силового кабеля. При таком способе оболочки кабелей соединены с общими контурами заземления на обоих оконечных пунктах магистрали. Сопротивление шины заземления, согласно норма-

и сии «¿j wo они iuüu —■ r\ e

r. »• тивам составляет не более 0.5

Рис.4. Распределение напряженности электрического поля q^ ЛЭП 110 кВ над поверхностью Земли на различных высотных

отметках Результаты расчета напря-

> --Се .............0.5 в ----------

: \ ! \

!N\

\

.....

•'х

женности электрического поля на примере подземной ЛЭП 110 кВ с заземлением приведены на рис.4.

Корректность результатов, получаемых при помощи МКЭ подтверждена сравнением расчетов поля воздушной ЛЭП с расчетами, проведенными по формулам (1) и (2). Данные результаты приведены в тексте диссертации.

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

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

Все методики расчетного прогнозирования электромагнитной обстановки, разработанные в настоящем разделе, реализованы в виде программы для ЭВМ в среде МаИаЬ.

Во втором разделе «Разработка электродинамических моделей силовых трансформаторных подстанций» выявлены и учтены особенности возникающей при этом электродинамической задачи. В частности, для уменьшения затрат времени на вычисления при реализации МКЭ размерность задачи целесообразно понизить. Нетрудно заметить, что анализируемые устройства имеют либо осевую (обмотки), либо плоскостную симметрию (трансформатор в целом), причем эта симметрия распространяется и на конфигурацию сторонних источников.

Е.У/т

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

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

Результаты расчета электрического поля в помещении типовой трансформаторной подстанции, содержащей два трансформатора, и помещениях второго этажа приведены на рис.6. При расчете предполагалось, что один из трансформаторов работает с нагрузкой 160 кВА, а второй с нагрузкой 50 кВА Картина распределения напряженности электрического поля построена для вертикального сечения модели. Из приведенных результатов видно, что электрическое поле концентрируется вблизи высоковольтных вводов и корпусов трансформаторов. При этом с самой трансформаторной уровни поля достаточно высоки (до 50 кВ/м).

Однако на втором этаже и выше уровень напряженности электрического поля невелик и значительно ниже предельно допустимых уровней и составляет примерно 0.2 кВ/м. Значительное ослабление электрического поля обусловлено экранирующим действием стен помещения, имеющего заземленную арматуру. Магнитное поле рассчитывалось на первом, и втором этажах в вертикальном сечении помещения. При расчете использовались следующие сведения о конструкциях трансформаторов: число витков трансформатора: высокого напряжения (ВН) - 2800 витков, низкого напряжения (НН) - 120 витков; средняя толщина обмотки трансформатора 160кВА:

Рис.6. Распределение электрического поля в трансформаторной подстанции

ВН - 320мм, НН -470мм. На рис.7 приведена картина распределения магнитного поля в вертикальном сечении.

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

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

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

21

Рис.7. Распределение магнитного поля в вертикальном сечении расчетного объема

писанному СанПиН 2.2.4.1191-03.

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

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

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

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

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

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

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

Моделирование ЭМП электротранспорта проведено при помощи метода конечных элементов.

Корпус электрического транспортного средства находится в непосредственной близости от токоведущих частей контактной сети. В связи с этим, очевидно, окажутся значительной величины вторичные поля, созданные вихревыми токами, наведенными в проводящих элементах кузова. Учет вихревых токов требует модификации исходных уравнений Лапласа.

Источниками вторичного поля являются токи, наведенные в толще металла переменным электрическим полем, созданным токами контактной сети. 22

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

Таким образом, следует подвергнуть модификации лишь уравнение относительно электрического вектора, исходя из предпосылки, что rot Ё * 0, следующим образом:

со scot

1-

tg S

rot j ,

(14)

где tg 8 - тангенс угла диэлектрических потерь материала корпуса.

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

На рис.8 приведены результаты расчета электрического поля, создаваемого цепью питания троллейбуса в присутствии вагона типа ЗиУ-9, оснащенного

1.1760е«ЮЗ

i.OSSfie* i.es26«+C03

J . *) J7S."I<j j

1.57ВЭе*СОЗ «.£ S««<

S.Glü8e> i . 11(3**C01 в ,5112eK00 3.77S6e*ÍOO

блзэгв-eoi

1.8430е-СШ . 1.3ü04e-C01 1 3 .Oí 37г-СО; 1.1Ш1-С02

5.C3 68e-C03 3.8157»-СОЭ l.Sllle-СОЗ

6.ВОЗЭя-СЗ :-CW

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

двумя тяговыми двигателями, суммарной мощностью 260 кВт. Троллейбус оснащен штанговым токоприемником.

Приведенные результаты подтверждают, что уровни электрического поля внутри салона троллейбуса достаточно высоки (0,2 кВ/м), и в ряде случаев приближаются к ПДУ, установленному для населения.

Относительно структуры поля отметим то обстоятельство, что области «сильного поля» и по электрической и по магнитной компонентам локализованы вблизи проводов контактной сети и практически не «касаются» салона. Уровни поля вне салона, оказываются небольшими и уже на расстоянии 1... 1,5 м от расположения контактной сети не превышают ПДУ. Присутствие вагона очевидным образом приводит к увели-

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

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

В разделе 4 «Экспериментальные исследования электромагнитных полей энергетической инфраструктуры региона» приведены результаты экспериментальных исследований ЭМП различных объектов системы энергоснабжения Самары и Самарской области.

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

Измерения проводились для ряда ЛЭП, эксплуатирующихся на территории Самарский области -ЛЭП-500 кВ, ЛЭП-200 кВ, ЛЭП-100 кВ, а также для некоторых трансформаторных подстанций. На рис.12 и 13 приведены некоторые результаты измерений.

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

Рис. 12. Электрическое поле ЛЭП 110 кВ

Рис. 13. Магнитное поле трансформаторной подстанции на высоте 0.5 м от плоскости пола

В разделе 5 «Систематизация подходов к региональному контролю природной среды по фактору электромагнитного излучения объектов энергетических систем» проведена разработка технологического алгоритма контроля. Предложен и детально обоснован геоинформационный подход к визуализации электромагнитной обстановки в регионе. Приведены скриншоты некоторых фрагментов ГИС электромагнитной безопасности региональной энергосистемы Самарской области.

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

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

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

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

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

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

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

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

В качестве примера на рис.14 и рис.15 приведены фрагменты электронной карты электромагнитного мониторинга г. Самары.

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

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

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

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

В заключении сформулированы основные научные и научно-практические результаты работы.

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

В том числе получены следующие научные результаты.

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

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

- группа локальных (сосредоточенных) источников.

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

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

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

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

26

Рис. 14. Слой карты с географическими координатами трансформаторных подстанций

Рис.15. Фрагмент карты с нанесенными ЛЭП, трансформаторными подстанциями и создаваемыми ими ЭМП

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

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

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

Все результаты исследований перечисленных объектов были в различное время получены в рамках хоздоговорных НИР и использованы в проектной документации на объекты.

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

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

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

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

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

7. Построены фрагменты геоннформационной системы электромагнитной безопасности энергетической инфраструктуры Самарской области.

В результате проведенных работ создано четыре слоя электронной карты:

- воздушные линии электропередач;

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

- трансформаторные подстанции;

- электромагнитные поля трансформаторных подстанций.

Таким образом, цели, поставленные в диссертационном исследовании, достигнуты полностью.

В приложениях приведены некоторые алгоритмы и результаты расчетов, не включенные в основную часть: в приложении 1 - выражения для функций формы и преобразования координат в методе конечных элементов; в приложении 2 -некоторые математические выкладки, результаты которых использованы в диссертации; приложение 3 содержит справочные данные по типовым опорам ЛЭП; в приложении 4 - формулы, используемые при расчете высокочастотных полей трансформаторов, а приложение 5 содержит результаты расчетов, не включенные в основную часть. В приложении 6 приведены акты внедрения результатов диссертационной работы.

СПИСОК ПУБЛИКАЦИЙ Монографии

1. Доб6ыш В.Н. Региональный контроль электромагнитных полей объектов энергетических систем: монография - Самара: Офорт, 2008. - 130 с.

2. Электромагнитная безопасность элементов энергетических систем: монография / Дов-быш В.Н., Маслов М.Ю., Сподобаев Ю.М. - Самара: «ИПК «Содружество», 2009. -198 с.

Нормативно-методические документы

3. Электромагнитные поля в окружающей среде. Расчет электромагнитных полей распределительных и оконечных устройств сетей энергоснабжения: методические указания / Довбыш В.Н., Маслов М.Ю., Ружников В.А. [и др.]. - Самара: ООО «САМБР», 2005. -57 с.

Статьи в журналах, рекомендованных ВАК

4. Довбыш В.Н. Моделирование и расчет электромагнитных источников тепла при индукционном нагреве осесимметричных тел сложной формы // Вестник Самарского государственного технического университета. Серия «Технические науки» - 2002. -Вып.14. -С. 181-184.

5. Довбыш В.Н. Оптимизация нестационарных режимов непрерывного индукционного нагрева в условиях ступенчатой вариации входных параметров изделий // Вестник Самарского государственного технического университета. Серия «Технические науки» -2002.-Вып.15-С. 18-23.

6. Электродинамическая модель трансформатора с сердечником, выполненным из материала с нелинейными свойствами / Довбыш В.Н., Маслов М.Ю., Ситником C.B. [и др.] // Инфокоммуникационные технологии - 2004. - Том 2. - №3. - С. 72-76.

7. Довбыш В.Н., Маслов М.Ю., Сарокваша О.Ю. Экологический мониторинг состава почвы в зоне размещения высоковольтной ЛЭП // Вестник Самарского государственного университета. Естественнонаучная серия. 2006. - №4 (44). - С. 136-148.

8. Доб6ыш В.Н. Анализ помехового излучения линии электропередач // Радиотехника. 2006. - №10. - С.69-71.

9. Довбыш В.Н., Сивков B.C., Сподобаев Ю.М. Визуализация электромагнитной обстановки, создаваемой телекоммуникационными техническими средствами, расположенными на больших территориях // Антенны. 2006. - №10 (113). - С.58-62.

10. Довбыш В.Н., Сивков B.C. Цифровая электромагнитная модель местности // Ин-фокоммуникационные технологии, 2007. - Т. 5. -№1. -С.85-88.

11. Довбыш В.Н. Электродинамическое моделирование цепей питания электротранспорта // Известия Самарского научного центра российской академии наук. Специальный выпуск «XIII конгресс Экология и здоровье человека». Самара, 2008. - Т. 1. — С.141-146.

12. Довбыш В.Н. Экологический анализ использования природных ресурсов и оценка состояния окружающей среды Самарской области // Известия Самарского научного центра Российской Академии наук, 2009. - Т. 11. - № 13. - С.472-476.

13. Довбыш В.Н. Экологический мониторинг элементов региональной энергетической инфраструктуры по фактору электромагнитного излучения // Известия Самарского научного центра Российской Академии наук, 2009. - Т. 11. - №13. - С.476-482.

14. Довбыш В.Н., Сподобаев Ю.М. Оценка радиочастотных помех, создаваемых короной высоковольтной линии электропередач // Телекоммуникации и транспорт. Спецвыпуск «Технологии информационного общества». Часть III. 2009. - С.77-79.

Публикации в прочих изданиях

15. Довбыш В.Н. Электромагнитные поля - фактор антропогенного загрязнения окружающей природной среды // Труды VII Всероссийского конгресса «Экология и здоровье человека». - Самара, 2001. - С.55-57.

16. Антропогенные загрязнения территории и их влияние на состояние здоровья населения / Довбыш В.Н., Котельников Г.П., Самыкина [и др.] // Труды VII Всероссийского конгресса «Экология и здоровье человека». - Самара, 2001. - С. 90-94.

17. Довбыш В.Н. Идентификация процесса непрерывного индукционного нагрева изделий при утилизации взрывателей // Труды двенадцатой межвузовской конференции «Математическое моделирование и краевые задачи»: ч.2. — Самара: Сам! ТУ, 2002. - С.31- 34.

18. Довбыш В.Н. Модели а алгоритмы оптимального управления индукционными нагревателями непрерывного действия в установках утилизации // Труды двенадцатой межвузовской конференции «Математическое моделирование и краевые задачи»: ч.2. - Самара: СамГТУ, 2002. - С.34- 37.

19. Данилушкин А.И., Довбыш В.Н. Проблемы экологии при утилизации взрывателей II Труды седьмой Международной конференции «Окружающая среда для нас и будущих поколений». - Самара: СамГТУ, 2002. - С.58- 59.

20. Абакумов A.M., Довбыш В.Н., Зимин JI.C. Энергосбережение и экология в электротехнологии // Труды седьмой Международной конференции «Окружающая среда для нас и будущих поколений». - Самара, СамГТУ, 2002. -С.13-14.

21. Довбыш В.Н. Проблемы и тенденции состояния природной среды Самарской области // Труды VIII Международного конгресса «Актуальные проблемы экологии человека». - Самара 2002. - С.64- 68.

22. Численный анализ электромагнитного поля силовых трансформаторов в реальных условиях размещения / Довбыш В.Н., Маслов М.Ю., Сподобаев Ю.М. // Ежегодник РНКЗНИ 2003. Сб. трудов. - М.: Изд-во АЛАНА, 2004. - С. 151- 163.

23. Довбыш В.Н. Практика управления отходами на территории Самарской области // Труды II Всероссийской научно-практической конференции «Процессы, технологии и оборудование для переработки отходов и вторичного сырья. Полигоны по захоронению отходов». - Самара, 2003. - С. 10- 13.

24. Довбыш В.Н. Проблемы электромагнитной экологии энергетических систем // Материалы IX Всероссийского конгресса «Экология и здоровье человека». Самара, 2004. - С. 89- 92.

25. Довбыш В.Н., Сподобаев Ю.М. Электромагнитный мониторинг энергетических систем // Материалы V Международной научно-технической конференции «Проблемы техника и технологии телекоммуникаций». - Самара, 2004. — С. 13—14.

26. Моделирование электромагнитных полей встроенной трансформаторной подстанции жилого дома / Довбыш В.Н., Маслов М.Ю., Сподобаев Ю.М. И Материалы V Международной научно-технической конференции «Проблемы техники и технологии телекоммуникаций». - Самара, 2004. - С.201 -202.

27. Встроенные трансформаторные подстанции - проблемы электромагнитного мониторинга / Довбыш В.Н., Маслов М.Ю., Сподобаев Ю.М. // Материалы IV Международного симпозиума по электромагнитной совместимости и электромагнитной экологии. -Санкт-Петербург. Июнь 2005. - С.25-30.

28. Довбыш В.Н. Электромагнитные поля и здоровье человека // Материалы X Конгресса «Экология и здоровье человека». - Самара, 2005. - С.88-91.

29. Довбыш В.Н., Сивков B.C. Применение геоинформационных технологий в системах электромагнитного мониторинга // Материалы X Конгресса «Экология и здоровье человека». - Самара, 2005. -С.91-93.

30. Довбыш В.Н., Сподобаев Ю.М. Проблемы мониторинга электромагнитного загрязнения окружающей среды // ХХХ1П Всероссийский семинар «Актуальные вопросы охраны окружающей среды. Программное обеспечение экологов, НПО «Логус». - Москва, 2005. - С.7-8.

31. Электромагнитный мониторинг энергетических систем / Довбыш В.Н., Маслов М.Ю., Сподобаев Ю.М. // Академия Энергетики, 2006. - Xsl (09). - С.4-7.

32. Довбыш В.Н. Принципы электромагнитного мониторинга энергетических систем // Материалы XIII Юбилейная Российская научная конференции профессорско-преподавательского состава, научных сотрудников и аспирантов. - Самара: ПГАТИ. -С. 150.

33. Довбыш В.Н., Маслов М.Ю. Анализ излучения коронного разряда высоковольтных линий передачи // Материалы VII Международной научно-технической конференции «Проблемы техники и технологии телекоммуникаций». - Самара, 2006. - С.270-272.

34. Довбыш В.Н., Сподобаев Ю.М. Проблемы электромагнитного мониторинга мегаполисов // Материалы VII Международной научно-технической конференции «Проблемы техники и технологии телекоммуникаций». - Самара, 2006. - С.272-274.

35. Довбыш В.Н. Современные методы оценки и прогнозирования экологической обстановки и масштабов природопользования И Труды XI Всероссийского конгресса «Экология и здоровье человека» на тему: «Проблемы выживания человека в техногенной среде современных городов». - Самара, 2006. - С.80-85.

36. Анализ влияния электромагнитных полей промышленной частоты на биохимический состав почвы в зоне прохождения высоковольтной ЛЭП / Довбыш В.Н., Маслов М.Ю., Сарокваша О.Ю. [и др.] // Ежегодник РНКЗНИ: сб.трудов, 2006. - М.: Изд-во АЛАНА, 2007. - С.111-126.

37. Довбыш В.Н., Сподобаев Ю.М. Уровни и источники электромагнитного загрязнения Самары // XIV Российская научная конференция профессорско-преподавательского состава, научных работников и аспирантов. - Самара: ПГАТИ, 2007. - С. 110-111.

38. Довбыш В.Н. Оценка экологической обстановки и масштабов природопользования с использованием ГИС // XIV Российская научная конф. профессорско-преподавательского состава, научных работников и аспирантов. - Самара: ПГАТИ, 2007. - С. 108.

39. Довбыш В.Н., Семаков ДМ. Расчет электромагнитных полей первичных объектов энергоснабжения области // XTV Российская научная конференция профессорско-

преподавательского состава, научных работников и аспирантов. - Самара: ПГАТИ, 2007. -С.110-111.

40. Довбыш В.Н., Сарохваша О.Ю. Эффект хронического воздействия электромагнитного излучения ЛЭП-110 кВ на эколого-биохимические показатели почвы // Материалы IV Международного симпозиума по электромагнитной совместимости и электромагнитной экологии. - Санкт-Петербург, 2007 — С.364-367.

41. Довбыш В.Н., Сивков B.C. формирование пространственных данных дан геоинформационной системы электромагнитной безопасности // Материалы VIIIМНТК «Проблемы техники и технологий телекоммуникаций». - Уфа, 2007. - С. 191-194.

42. Довбыш В.Н., Сподобаев МЛО. Современные методы оценки и прогнозирования экологической обстановки // Материалы VIII МНТК «Проблемы техники и технологий телекоммуникаций». - Уфа, 2007.- С. 196-199.

43. Довбыш В.Н., Маслов М.Ю., Сподобаев Ю.М. Метод оценки электромагнитной обстановки в помещениях со встроенными распределительными пунктами систем энергоснабжения // Материалы VIII МНТК «Проблемы техники и технологий телекоммуникаций». - Уфа, 2007. - С.207-210.

44. Довбыш В.Н., Маслов М.Ю., Сарокваша О.Ю., Сподобаев Ю.М. Анализ влияния электромагнитных полей промышленной частоты на биохимический состав почвы в зоне прохождения высоковольтной ЛЭП // Ежегодник РНКЗНИ: сб. трудов, 2007 г. - М.: Изд-во АЛАНА, 2007. -С.111-155.

45. Довбыш В.Н., Сивков B.C., Сподобаев Ю.М. Визуальное представление электромагнитной обстановки на территории мегаполиса // Труды XIII международной научно-технической конференции «Радиолокация, навигация, связь», 2007 -Т2. - С.2036-2042.

46. Довбыш В.Н., Маелов М.Ю., Сподобаев Ю.М. Региональный мониторинг электромагнитных полей создаваемых элементами энергетических систем // Труды VII международного научно-практического семинара (¿Проблемы электромагнитной экологии в науке, технике и образовании». - Ульяновск, 2008. - С.22-25.

47. Довбыш В.Н. Расчет электромагнитных полей, создаваемых цепями питания электротранспорта II Материалы VIII Международного симпозиума по электромагнитной совместимости и электромагнитной экологии - ЭМС-2009. - Санкт-Петербург, 2009. - С.27-30.

48. Довбыш В.Н., Маслов М.Ю., Сподобаев Ю.М. Электромагнитная безопасность элементов энергетических систем // Материалы VIII Международного симпозиума по электромагнитной совместимости н электромагнитной экологии - ЭМС-2009. - Санкт-Петербург, 2009. - С27-30.

1 49. Довбыш В.Н., Маслов М.Ю., Сподобаев Ю.М. Энергетические системы: электромагнитный мониторинг и безопасность Л Материалы Второго международного экологического конгресса - ELPIT-2009. Тольятти, 2009. - С.46-50.

Государственное образовательное учреждение высшего профессионального

образования «Поволжский государственный университет телекоммуникаций и

информатики» 443010, г. Самара, ул. Льва Толстого 23.

Отпечатано фотоспособом в соответствии с материалами, представленными

__заказчиком_

Подписано в печать 6.09.10 г. Формат 60x84''/16 Бумага писчая№ 1. Гарнитура Тайме.

_Заказ 741. Печать оперативная .Усл. печ, л. 1.84. Тираж 100 экз._

Отпечатано в издательстве учебной и научной литературы Поволжского государственного университета телекоммуникаций и информатики 443090, г. Самара, Московское шоссе 77. т. (846) 228-00-44

ВВЕДЕНИЕ.

1. РАЗРАБОТКА ЭЛЕКТРОДИНАМИЧЕСКИХ МОДЕЛЕЙ РАСПРЕДЕЛИТЕЛЬНЫХ ЭЛЕМЕНТОВ РЕГИОНАЛЬНОЙ ЭНЕРГЕТИЧЕСКОЙ ИНФРАСТРУКТУРЫ

1.1. Классификация технических средств систем энергоснабжения с точки зрения пространственной локализации и характеристик создаваемых электромагнитных полей.

1.2. Электродинамическое моделирование воздушных линий электропередач

1.2.1. Выбор и обоснование подходов к электродинамическому моделированию

1.2.2. Вывод выражений для компонент векторов электрического и магнитного полей.

1.2.3. Учет разветвленного характера воздушных линий.

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

1.2.5. Расчет высокочастотных полей высоковольтных линий.

1.3. Электродинамическое моделирование подземных линий электропередач

1.3.1. Анализ поля частично экранированной многопроводной линии при помощи метода конечных элементов.

1.3.2. Моделирование линий электропередач и оценка корректности конечноэлементных моделей.

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

1.4. Анализ электромагнитных полей силовых распределительных пунктов системы энергоснабжения.

1.4.1. Моделирование силовых распределительных пунктов расчетно-экспериментальным методом.

1.4.2. Выражения для компонент векторов электромагнитного поля расчетно-экспериментальной модели.

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

1.5. Выводы по разделу

2. РАЗРАБОТКА ЭЛЕКТРОДИНАМИЧЕСКИХ МОДЕЛЕЙ СИЛОВЫХ ТРАНСФОРМАТОРНЫХ ПОДСТАНЦИЙ.

2.1. Моделирование силового электрооборудования с учетом условий размещения при помощи метода конечных элементов.

2.1.1. Исходные уравнения и формулировка соответствующих краевых задач.

2.1.2. Выбор и обоснование вида аппроксимирующего элемента и порядка интерполяционного полинома.

2.2. Построение модели трансформатора и трансформаторной подстанции

2.3. Вывод выражений для коэффициентов интерполяционных полиномов методом Галеркина.

2.4. Вывод выражения для баланса энергий и оценка точности вычислений

2.5. Учет особенностей структуры матричного уравнения, возникающего в методе конечных элементов.

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

2.7. Расчет высокочастотных полей силовых трансформаторов.

2.7.1. Общие замечания.

2.7.2. Построение проволочной модели силового трансформатора.

2.7.3. Применение метода ИУ к электродинамическому моделированию силовых трансформаторов.

2.7.4. Численное решение интегрального уравнения.

2.7.5. Результаты расчета характеристик излучения высоковольтного силового трансформатора.

2.8. Выводы по разделу 2.

3. ЭЛЕКТРОДИНАМИЧЕСКОЕ МОДЕЛИРОВАНИЕ ЦЕПЕЙ ПИТАНИЯ ЭЛЕКТРОТРАНСПОРТА.

3.1. Варианты технических решений и конструктивного исполнения систем питания и энергоснабжения электротранспорта.

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

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

3.3.1. Электромагнитные поля системы питания троллейбуса.

3.3.2. Электромагнитные поля системы питания трамвая.

3.3.3. Электромагнитные поля системы питания метрополитена.

3.4. Выводы по разделу 3.

4. ЭКСПЕРИМЕНТАЛЬНЫЕ ИССЛЕДОВАНИЯ ЭЛЕКТРОМАГНИТНЫХ ПОЛЕЙ ЭНЕРГЕТИЧЕСКОЙ ИНФРАСТРУКТУРЫ РЕГИОНА.

4.1. Методика проведения эксперимента и обработки результатов.

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

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

4.5. Выводы по разделу 4.

5. СИТЕМАТИЗАЦИЯ ПОДХОДОВ К РЕГИОНАЛЬНОМУ КОНТРОЛЮ ПРИРОДНОЙ СРЕДЫ ПО ФАКТОРУ ЭЛЕКТРОМАГНИТНОГО ИЗЛУЧЕНИЯ ОБЪЕКТОВ ЭНЕРГЕТИЧЕСКИХ СИСТЕМ.

5.1. Алгоритм комплексного анализа электромагнитной обстановки в регионе

5.2. Применение геоинформационных технологий к представлению результатов электромагнитного мониторинга.

5.3. Результаты электромагнитного мониторинга энергетической системы города Самара и Самарской области.

5.4. Выводы по разделу

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

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

Многообразие переплетающихся между собой процессов, связывающих человеческий организм и среду его обитания, требуют комплексной оценки последствий как непреднамеренного воздействия на окружающую среду, так и целенаправленного преобразования природы [2, 199]. Поэтому решение современных экологических проблем немыслимы без участия практически всех областей научного знания и отраслей техники.

Анализ народнохозяйственных планов развития большинства стран планеты, проведённый в соответствии с программой ООН по проблемам окружающей среды [17, 20, 29, 30, 82, 95, 150, 182], показывает, что интегральные экономические показатели государств, а так же их тенденции на ближайшие десятилетия, по многим позициям не обеспечены имеющимися на планете ресурсами. В связи с этим важным обстоятельством, выбор экологической стратегии и формирование на ее основе общественного экологического сознания является на сегодняшний день необходимым, хотя и не достаточным, условием для устойчивого (квазиустойчивого) развития человеческого общества.

Традиционно к задачам экологии относят следующие направления:

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

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

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

При этом очевидно, что достаточно устойчивое состояние экосистем возможно только при непрерывном системном контроле состояния природной среды по всем существенным факторам антропогенного воздействия.

Человеческое общество в процессе своей деятельности, направленной, прежде всего, на создание комфортных условий существования воздействует на различные компоненты природной среды: атмосферу, гидросферу и литосферу. В настоящее время такие воздействия приобретают глобальный характер, затрагивая все, без исключения, континенты нашей планеты. Анализ значительного количества публикаций геоэкологической тематики отечественного [17, 20, 3032, 35, 44, 99, 103, 105, 150, 157, 171, 187, 193] и зарубежного [200, 204, 205, 216, 217] происхождения, позволяет отметить то обстоятельство, что среди всего многообразия видов антропогенного загрязнения природной среды, можно выделить один специфический вид энергетического загрязнения, а именно, электромагнитный — вид, влияние которого можно одновременно отнести практически ко всем компонентам природной среды одновременно.

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

И в быту, и в процессе трудовой деятельности человека окружают разнообразные технические средства, создающие электромагнитные поля (ЭМП), которые обладают различными пространственно-временными характеристиками. Причем для одних технических средств генерация электромагнитной энергии является специфической особенностью, диктуемой их функциональным назначением, а для других - напротив, побочным явлением. Однако, в обоих случаях генерируемые поля являются активным фактором загрязнения окружающей среды. Отмеченные вопросы относятся к специфической области знания -«электромагнитной экологии» [26, 44, 129, 131, 132, 141, 182-185, 187, 193, 202-204, 217] и традиционно проблемы с ними связанные решаются при помощи электромагнитного мониторинга включающего в себя:

- расчетное прогнозирование электромагнитных полей, что весьма важно для стадий разработки, проектирования и размещения технических средств, являющихся источниками ЭМП [24, 129, 138, 182-185];

- инструментальный контроль электромагнитной обстановки на стадии эксплуатации объектов и их комплексов [44, 80, 128, 131, 141, 217];

- разработку мероприятий и рекомендаций по защите от ЭМП и нормализации электромагнитной обстановки [182, 185, 187].

Современному этапу развития человеческого общества, соответствует региональный характер развития всех инфраструктур [142, 150, 179], в том числе и связанных с потреблением и излучением электромагнитной энергии. Действительно, все энергоемкие объекты промышленного производства и культурной деятельности человека, в значительной степени, сосредоточены в локальных областях, поэтому требованиям более или менее полного контроля над состоянием природной среды по электромагнитному фактору, может отвечать только система, основанная на регионально-ориентированном подходе.

Следует отметить, что в задачах электромагнитного мониторинга зачастую необходим анализ поля в непосредственной близости от технических средств с учетом реальных особенностей их размещения и наличия других технических средств и материальных тел [24, 182-185], что накладывает известные трудности на корректное проведение эксперимента и воспроизводимость результатов, получаемых эмпирическим путем. Это обстоятельство, очевидно, и обуславливает повышенный, особенно в последние десятилетия, интерес к созданию методик расчетного прогнозирования электромагнитной обстановки. Инструментальные же методы традиционно используются при проверке корректности расчетных методик, а так же в случаях, когда получение исчерпывающей информации об объекте исследования, необходимой для построения корректной теоретической модели, невозможно. Данное обстоятельство существенно отличает электромагнитную экологию от других экологических направлений, в которых для оценки состояния природной среды используют в основном данные, получаемые в результате экспериментов.

Вопросы расчетного прогнозирования ЭМП излучающих технических средств телекоммуникаций достаточно хорошо изучены как в нашей стране [2525, 179, 182-185], так и за рубежом [197, 198, 206, 207, 214]. Вопросам же загрязнения окружающей среды и контроля экологической безопасности по фактору электромагнитных полей промышленной частоты в рамках решения общих проблем электромагнитной экологии уделялось явно недостаточное внимание. Существование такого «белого пятна» в электромагнитной экологии можно оправдать только тем обстоятельством, что строительство объектов энергоснабжения различного назначения и увеличение их характерных энергетических нагрузок приняло широкомасштабный характер именно в последние годы, и только сегодня эти проблемы приобрели особенную актуальность и социальную значимость.

Ежегодно в регионах России вводятся в эксплуатацию новые и реконструируются существующие объекты и технические средства системы энергоснабжения. Этот процесс принял в последнее десятилетие лавинообразный характер. Отличительной особенностью современного этапа развития региональных энергетических инфраструктур России является многократный- рост мощностей, характерных для типичных циклов жизни регионов и мегаполисов.

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

Известно [86, 132, 185, 204], что энергетическое оборудование, в частности, линии электропередач, сети питания и тяговые подстанции электротранспорта, силовые трансформаторы, силовые распределительные пункты создают электромагнитные поля промышленной частоты, которые вносят существенный, а зачастую и определяющий вклад в общую электромагнитную обстановку на селитебных территориях [2, 171, 183, 185, 217].

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

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

В то же время следует отметить, что в то время как для излучающих технических средств телекоммуникаций существует развитая система санитарной паспортизации, для энергетического оборудования подобной системы не существует, а контроль электромагнитной обстановки в настоящее время проводится эпизодически при помощи методик, лишь в некоторых случаях имеющих статус отраслевых стандартов [20, 132, 150]. Так же, для оценки состояния природной среды по факторам электромагнитного излучения телекоммуникационного оборудования в нашей стране создан ряд автоматизированных программных комплексов анализа электромагнитной обстановки [184, 185]. Для энергосистем подобных программных средств не существует.

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

Особая народнохозяйственная значимость проблем контроля среды пребывания человека по электромагнитному фактору подчеркнута Решением Коллегии Федеральной службы Роспотребнадзора от 27 ноября 2009 г., в котором определены задачи и приоритетные направления по детальному изучению электромагнитной обстановки в населенных пунктах, инвентаризации источников электромагнитных излучений, предусматривающие планомерное снижение электромагнитной нагрузки на население и создание баз данных.

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

Состояние вопроса в рассматриваемой области характеризуется следующими основными достижениями.

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

Так методологические основы технологий регионально контроля природной среды развиты в работах Булгакова Н.Г., Левича А.П. и Максимова В.Н. [29, 135, 199, 200].

Общие вопросы электромагнитного мониторинга антропогенных воздействий достаточно хорошо освещены как в отечественной [2, 12, 17, 20, 30-32, 35, 51, 97, 99, 103, 106, 124, 142, 150, 157, 159, 168, 171, 179, 181, 182 и др.], так и зарубежной [200, 202, 203, 205, 216, 217] литературе. Общие подходы, к сбору и обработке экологической и геоэкологической информации, разработанные и представленные в указанных работах могут быть с успехом использованы для целей настоящего диссертационного исследования.

Проблемы, методы и средства численного анализа технических средств, являющихся источниками ЭМП, в том числе и тех, для которых данное свойство не продиктовано функциональным назначением (так называемых нетрадиционных источников излучения), также достаточно полно освещены в следующей литературе [4, 5, 9, 18, 21, 23-25, 27, 41, 42, 49, 50, 101, 111, 116-120, 122, 138, 146, 147, 184, 185, 190, 194, 195, 196-198, 201,206-208,211-214].

Следует отметить, что большинство подзадач, сходных с возникающими в диссертационной работе, решено в рамках традиционной электромагнитной экологии [23-27, 97, 128, 129, 138-141, 182-185, 187, 217]. При этом большинство авторов используют методы математического моделирования, хорошо себя зарекомендовавшие при решении задач вычислительной электродинамики и теории антенн [3-5, 9, 18, 21, 23, 27, 36, 40-42, 49, 50, 101, 102, 110, 111, 116-120, 122, 137, 138, 146-149, 151-156, 160, 161, 170, 180, 184-186, 195, 196198, 201,206-208,211-215].

Исследования биологического действия электромагнитных полей промышленной частоты, выполненные в нашей стране в 60-70х годах прошлого века [45, 132], ориентировались в основном на действие электрической составляющей, поскольку экспериментальным путем значимого биологического действия магнитной составляющей при типичных уровнях не было обнаружено. В 70-х годах для населения по электрическому полю промышленной частоты были введены нормативы, являющиеся и по настоящее время одними из самых жестких в мире [43]. Впоследствии, в 80-е годы XX столетия, на основании массовых эпидемиологических обследований населения, проживающего в условиях облучения магнитными полями высоковольтных линий, как безопасный уровень для условий продолжительного облучения, не приводящий к онкологическим заболеваниям, независимо друг от друга шведскими и американскими специалистами была рекомендована величина плотности потока магнитной индукции [204, 217]. Этот зарубежный стандарт и послужил основой для введения норматива по магнитному полю в нашей стране. Данными сведениями, практически исчерпываются современные достижения в области анализа биологического воздействия ЭМП энергетического оборудования. Тем не менее, в качестве критериев оценки электромагнитной обстановки, в настоящей работе могут выступать именно результаты упомянутых исследований.

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

Однако большинство работ, выпущенных ранее 70-х годов, посвящены в основном, проблемам биологического воздействия и гигиенического нормирования. Фундаментальными же в области расчета ЭМП различных технических средств стали работы Шередько Е.Ю., Сподобаева Ю.М., Кубанова В.П., Маслова О.Н., Бузова A.JL, Романова В.А., Казанского Л.С. и др. [19, 23, 28, 74-77, 118-120, 150156]. В этих работах предложены и обоснованы подходы к расчетному прогнозированию электромагнитной остановки вблизи широкого класса излучающих технических средств и их комплексов, сформулированы подходы к системному электромагнитному мониторингу. Основные результаты этих работ, подтвержденные многочисленными экспериментальными исследованиями, нашли отражение в нормативно-методических документах, утвержденных государственными органами санитарно-эпидемиологического надзора [140-144]. Так же на основе разработанных методик создан известный программный комплекс [185].

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

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

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

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

- разработка методик измерения ЭМП, для оценки корректности предложенных теоретических моделей;

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

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

Предварительная оценка качественного состава источников ЭМП, входящих в энергетическую инфраструктуру региона позволяет сделать ряд заключений. Так, основными техническими средствами - «поставщиками» электромагнитной энергии, существенно влияющей на общую электромагнитную обстановку в масштабах региона, являются линии электропередач (ЛЭП), распределительные пункты системы энергоснабжения, силовые трансформаторные подстанции, а так же силовые установки и сети питания наземного и подземного электротранспорта и т.п. [217]. В связи с этим можно утверждать, что все рассматриваемые технические средства можно подразделить на две основные группы, исходя из особенностей их пространственной локализации, следующим образом:

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

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

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

Вопросы расчетов ЭМП электроустановок подробно рассмотрены в работах Колечицкого Е.С. [116 - 120], Меликова H.A. [102], Филиппова A.A. [120], Тозони О.В. [190], Демирчана К.С. [49], Abou-Seada M.S. и Nasser Е. [196], Silvester Р. и Chari М. [214], а так же других авторов.

Известно, что любая электродинамическая проблема, в конечном счете представляет собой краевую электродинамическую задачу для системы дифференциальных уравнений в частных производных - уравнений Максвелла, или иных, из них выводимых - однородных или неоднородных уравнений д'Аламбера и Гельмгольца, уравнений Пуассона или Лапласа. Таким образом, принципиальная возможность расчета ЭМП любым из известных методов возникает в случаях, в которых может быть сформулирована (явно, или неявно) соответствующая полная система граничных условий [216, 156]. Из них, применительно, к задачам, поставленным в диссертационном исследовании, очевидно, необходимо рассмотреть следующие. Первое и самое распространенное граничное условие состоит в том, что должны быть известны первичные потенциалы и токи во всех проводниках (на их поверхностях или в их объемах, в зависимости от конкретной постановки задачи), образующих модель. В некоторых задачах граничное условие сводится к тому, что значение потенциала или тока проводника подлежит определению при условии, что известен его полный заряд. Так, например, в конструкциях, содержащих промежуточные экраны, не имеющие гальванического контакта с токоведущими частями устройства, заряд экранов обычно полагают равным нулю. Указанные граничные условия могут быть заданы в различных комбинациях. В соответствующей литературе принято разделение статических и стационарных полей на плоские (двумерные), плоскомеридианные (псевдодвумерные) и трехмерные [49, 117, 214]. Такой способ классификации может быть распространен и на соответствующие расчетные модели и, в значительной степени, определить выбор адекватных, для решения конкретной задачи, расчетных методов. Так, поля, образованные объемными проводниками и содержащие тонкие незамкнутые витки тока, достаточно эффективно могут быть определены в результате решения задачи в двумерной псевдодвумерной) формулировке (плоские и плоскомеридианные поля) [117, 214]. Расчет трехмерных полей таких проводников хотя и не представляет принципиальных математических трудностей, как правило, производится со значительно большими затратами машинного времени.

Практически все существующие на сегодняшний день методы расчета статических и стационарных полей вполне обеспечивают принципиальную возможность расчета потенциалов и напряженностей полей, впоследствии успешно верифицируемых. Однако, при применении различных методов оказывается существенно различным объем информации, получаемой при решении задачи, и, кроме того, различен объем вычислений, необходимых для получения численных значений искомых функций [214]. В связи с данным обстоятельством, при выборе адекватного расчетного метода, важно учесть, что при анализе поля наибольший интерес, как правило, представляет сравнительно небольшая область пространства вблизи источника. В большинстве рассматриваемых в настоящей работе задач, это область вблизи точки, в которой искомая функция принимает максимальное значение - так называемая область «сильного поля». Такой подход, позволяет ограничить область анализа до сравнительно небольших размеров, снижая при этом общую ресурсоемкость задачи, в смысле потребности производительности центрального процессора и объема оперативной памяти ЭВМ.

Решение любой задачи, по расчету ЭМП может производиться аналитически или при помощи численных методов вычислительной электродинамики [36, 154]. Содержание известных аналитических методов расчета статических и стационарных полей подробно изложено в [146, 153]. Их эффективное применение возможно в тех случаях, когда форма проводников, несущих первичные токи может быть достаточно простым образом представлена в какой либо системе координат (декартовой, цилиндрической и т.п.). Применительно к задачам, поставленным в настоящей диссертационной работе, аналитический подход целесообразен при расчете поля распределенных технических средств — линий электропередач, цепей питания электротранспорта. При этом искомые выражения для компонент векторов поля могут быть получены их известных интегралов уравнений Пуассона и Лапласа. Исключения составляют задачи анализа ЭМП подземных ЛЭП, а так же цепей питания электротранспорта в непосредственной близости транспортного средства. В первом случае на структуру и уровни поля оказывают существенное влияние реальные условия размещения, во втором случае - расположенный вблизи токоведущих частей проводящий корпус транспортного средства. В указанных задачах целесообразно уточнение решения при помощи какого-либо численного метода.

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

В настоящее время в вычислительной электродинамике наибольшее распространение получили следующие численные методы [49, 207,214]: метод сеток, или метод конечных разностей (МКР) [40, 129, 190]; вариационные методы (ВМ) [111, 122]; метод конечных элементов (МКЭ) [38, 49, 50, 102, 145, 156, 178, 207, 214]; метод интегральных уравнений (МИУ) [21, 27, 36, 41, 116-120, 152, 170, 194, 195, 197,198,206,208,211,212]; метод эквивалентных зарядов (МЭЗ) [119,190,196].

Применительно к статическим и стационарным задачам, упомянутыми методами осуществляется численное решение уравнений Пуассона или Лапласа или эквивалентных этим уравнениям задач. Численная процедура всех этих методов сводится к составлению и решению системы линейных уравнений [36, 195, 207]. Все необходимые характеристики поля в дальнейшем могут быть вычислены на основании решения соответствующей системы линейных уравнений. Различные методы различаются между собой способом составления такой системы, видом и размерностью матрицы коэффициентов этой системы, а также способом учета граничных условий.

Названные численные методы расчета ЭМП поля обладают различными функциональными возможностями, для их применения нужны разные способы подготовки исходных данных, различается для них и время счета. Следует отметить, что задача объективного выбора расчетного метода, наиболее пригодного для решения определенного класса задач электродинамики, до настоящего времени не может считаться окончательно решенной. Такое положение объясняется целым рядом причин. Одной из них является отсутствие общепризнанных критериев сравнения эффективности методов для однотипных задач [117]. Второй причиной служит тот факт, что один и тот же метод может быть реализован различным образом. И, наконец, невозможно утверждать, что возможности какого-либо метода на сегодняшний день раскрыты исчерпывающим образом, и каждый из них имеет свои перспективы развития.

Рассмотрим различные методы вычислительной электродинамики с точки зрения их функциональных возможностей и целесообразности применения для решения поставленных в диссертации задач. В этом плане наиболее универсальными представляются МКР и МКЭ. Оба этих метода позволяют описывать поверхности границ раздела сред весьма совершенно, так как в них присутствует принципиальная возможность аналитического описания любых, сколь угодно сложных поверхностей. Учет объемного заряда и наличия внешнего поля обоими методами производится весьма эффективно [116, 207]. При этом собственно характер его распределения не играет существенной роли - он может быть задан в виде точечных, объемных или поверхностных распределений. Точно так же внешнее поле может быть описано аналитической функцией произвольного вида, или задано таблично.

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

Методы решения подобных электродинамических задач можно разделить на три группы. К первой группе можно отнести прямые методы, подробно описанные в работах Говоркова В.А. [40] и монографии Миролюбова H.H. и др. [146], основанные на интегрировании тока вдоль криволинейного контура, проходящего по осям всех токоведущих частей, и предполагающие получение решения в замкнутой аналитической форме. Применительно к настоящему случаю данный подход малоприменим, поскольку витки обмоток имеют сложную форму, специфическую для каждого устройства, а так же различную для каждого слоя намотки, что затрудняет вычисление возникающих при этом интегралов. Кроме того, практически невозможным становится учет нелинейных элементов конструкции.

Методы, которые мы отнесем ко второй группе, предполагают сведение исходной краевой задачи к интегральному уравнению (ИУ) (системе уравнений). Решение ряда подобных задач указанными методами продемонстрировано в монографиях Тозони О.В. [190], Колечицкого К.С. [216], а так же Никольского В.В. и Никольской Т.Н. [155]. Возможных реализаций МИУ, применительно к решению поставленной задачи, возможно множество, при этом искомой функцией может выступать и магнитный поток в магнитопроводе (эквивалентный магнитный ток) и тангенциальные компоненты векторов электрического и магнитного поля. Методы ИУ могут быть ориентированы как на аналитическое, так и на численное решение, и позволяют учитывать нелинейности. С точки зрения задач поставленных в диссертации, основным недостатком МИУ является их неуниверсальность в смысле конфигурации токоведущих обмоток, т.е. для каждого конструктивного элемента необходимо получение отдельного уравнения. Причем свойства этих уравнений могут быть различными из-за геометрических особенностей входящих в них подынтегральных выражений. Подобный недостаток имеют и прямые методы, относимые к первой группе.

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

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

Как отмечалось выше, численное решение дифференциальных уравнений удобно проводить либо МКР, либо МКЭ. Характерной особенностью обоих методов решения дифференциальных уравнений является весьма значительная вычислительная трудоемкость. Обзор и сравнительная характеристика разностных и конечно-элементных методов решения электродинамических задач дана в работе [207]. Отмеченная работа является обобщением опыта коллектива авторов, занимавшихся разработкой универсальных программных пакетов электродинамического моделирования, таких как XFDTD, производства компании Remcom - HFSS и Maxwell, производства фирмы Ansoft. В работе отмечается преимущество метода конечных элементов при решении стационарных задач как более гибкого и экономичного. Действительно, практически все программные комплексы, ориентированные на решение статических и стационарных задач, используют именно данный метод.

Метод конечных элементов, традиционно применяемый для решения дифференциальных уравнений в задачах сопротивления материалов, строительной механики [38, 39, 102, 156], теории упругости и теплопроводности [102, 113, 178], в электродинамике появился сравнительно недавно. Причина в существенной ресурсоемкости многомерных задач, к коим относятся задачи электродинамики.

Сущность данного метода состоит в том [156], что анализируемая часть пространства разбивается на подобласти, в пределах которых решение аппроксимируется функциями- специального вида. Сшивание отдельных функций на границах обеспечивается граничными условиями.

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

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

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

- распределительные устройства сетей онергоснабэ!сения — объекты систем энергоснабжения, осуществляющие распределение электрической энергии, вырабатываемой первичными источниками (ГЭС, АЭС и т.д.), между потребителями всех уровней сетевой иерархии, кроме конечных - высоковольтные линии электропередач (ВЛЭП, ЛЭП); класс напряжения которых не ниже 6 кВ; потребительские цепи низкого напряжения; линии питания электротранспорта;

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

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

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

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

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

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

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

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

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

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

Так для высоковольтных линий электропередач, электромагнитное поле которых создается протяженными участками многопроводных линий, целесообразно применять метод расчета, основанный на использовании интегралов уравнений Максвелла, известных в замкнутой форме. Влияние подстилающей поверхности (поверхности Земли) на структуру и уровни ЭМП при этом учитываются введением соответствующих зеркальных изображений первичных токов [204]. Использование такого относительно простого подхода к построению электродинамических моделей оправдано тем обстоятельством, что ЛЭП содержат преимущественно линейные коллинеарные токи.

Электрическое поле ЛЭП практически не зависит от нагрузки и полностью определяется напряжением. Магнитное же поле, напротив существенно зависит от нагрузки и при расчете необходимо учитывать изменяющийся характер последней. Для оценки экологической безопасности ЛЭП необходим расчет магнитного поля с учетом максимальной за год средней нагрузке. Стандартные расчеты ЭМП ВЛЭП целесообразно проводить с учетом последующего сравнения результатов с критериями принятыми в СанПиН [174-177] следующим образом:

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

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

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

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

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

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

Электрическое поле трансформаторных подстанций практически полностью экранируется арматурой помещения подстанции [217], поэтому расчет электрического поля необходимо проводить только внутри технического помещения. Магнитное поле свободностоящих (внешних) трансформаторных подстанций рассчитывается в горизонтальном сечении на высоте 2 м от поверхности Земли до расстояний, превышающих размеры зоны, в которой контролируемые параметры оказываются больше своих предельно допустимых значений, регламентированных СанПиН [174-177]. Магнитное поле встроенных трансформаторных подстанций рассчитывается в техническом помещении, а так же в примыкающих помещениях на первом и далее этажах в горизонтальных сечениях, в плоскостях отстоящих на 0.5 и 1.5 м от уровня поля.

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

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

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

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

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

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

Учет влияния материальных тел, при анализе указанных выше исключительных случаев следует проводить отдельно, руководствуясь индивидуальными соображениями.

Следует также отметить, что воздушные линии электропередач и силовые трансформаторы, помимо ЭМП ПЧ, являются источниками высокочастотных полей. Причинами данных явлений являются, соответственно, коронирование проводов [109] и частичные разряды в обмотках [169].

Оценка уровней электромагнитных полей, создаваемых короной ЛЭП, возможна в силу простоты конфигурации первичных токов, при помощи приближенного подхода, развитого в работах Казакова В.Н. [109], а так же Кима К.С., Лелевкина В.М., Токарева A.B., Юданова В.А. [115].

В случае с силовыми трансформаторами, реальные устройства могут быть представлены для моделирования в виде сетки проводников, с сосредоточенным возбуждением. Расчет распределения тока на проводниках такой модели требует решения одномерного интегрального уравнения (ИУ). Подобные задачи широко известны в соответствующих областях теории антенн. Применительно к целям, поставленным в диссертационной работе, наиболее удобным и перспективным представляется подход на основе так называемого тонкопроволочного приближения [27, 36, 208]. Действительно, сетка, аппроксимирующая реальное техническое средство, состоит из проводников, размеры и форма поперечного сечения которых не имеют значения, а методы решения электродинамических задач на основе тонкопроволочного приближения достаточно хорошо изучены и относительно легко алгоритмизируются.

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

Как отмечалось выше, данные об электромагнитной обстановке в масштабах региона представляют собой массивы данных, значительного объема. При этом весьма существенным обстоятельством, отличающим такие данные от результатов, получаемых в иных формах экологического мониторинга, является «генетическая» привязанность к географическим координатам. Действительно, области «сильного поля» энергетического оборудования, очевидно оказываются локализованными вблизи мест расположения технических средств - источников, а результирующая электромагнитная обстановка образует сложную пространственную картину, привязанную к рельефу местности.

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

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

Берлянта A.M. [16], Полякова М.М. [157], Прогуловой Т.Б. [168], Мясоедова Б.В. [179], а так же ряда других авторов. Общие подходы и примененные способы представления векторных и матричных данных на электронных картах могут быть с успехом применены для реализации целей, поставленных в настоящей диссертации.

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

Для достижения поставленной цели в настоящей диссертационной работе выполнена следующая программа исследований.

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

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

3. Исследование электромагнитных полей реальных объектов, функционирующих в Самарской области.

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

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

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

7. Построение фрагментов геоинформационной системы электромагнитной безопасности энергетической инфраструктуры Самарской области.

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

5.4. Выводы по разделу 5

В данном разделе проведен анализ критериев оценки электромагнитной обстановки в регионе с точки зрения целей электромагнитной безопасности.

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

Рассмотрены принципы визуализации данных экологического мониторинга с использованием геоинформационных технологи.

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

Результаты, приведенные в настоящем разделе, опубликованы в [61, 63, 67, 68-70, 73, 76, 79, 80, 82, 84, 85, 88, 89-92, 95-97]. я о Л

К) к м о Й к я я я р 9 я я н я о ч о я о » аз О я тз

-1 ю ш Я о о о я о ЯС я о 1=1 о н ю я я я я

Х= 5 450,90 ш У= 9 207,51 т

ССКИИ

Карта 2000 - [D:\ASPIRANT\GEOINFGOROD\maps\домики.МАР] домики" 1:100 000 Листов :1

Объектов: 7527610 (отображено I выделено)

1 : 500 [домики (объектов: 75277)

0 Файл Правка Вид Поиск Задачи База Масштаб Параметры Окно Помощь

В диссертационной работе получены следующие научные результаты.

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

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

- группа локальных (сосредоточенных) источников.

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

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

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

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

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

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

7. Построены фрагменты геоинформационной системы электромагнитной безопасности энергетической инфраструктуры Самарской области.

В результате проведенных работ создано четыре слоя электронной карты:

- воздушные линии электропередач;

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

- трансформаторные подстанции;

- электромагнитные поля трансформаторных подстанций.

Таким образом, цели, поставленные в диссертационном исследовании, достигнуты полностью.

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

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

Итоговые результаты экологического контроля состояния природной среды по фактору ЭМП, создаваемого элементами региональной энергосистемы, представлены в виде слоев электронной карты, вошедших в ГИС ЕЦКО.

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

1. Абакумов A.M., Довбыш В.Н., Зимин Л.С. Энергосбережение и экология в электротехнологии. // Труды VII Международной конференции «Окружающая среда для нас и будущих поколений». Самара, СамГТУ, 2002, С. 13-14.

2. Артемов И.Е. Влияние современной антропогенной деятельности на природную среду и климатические особенности Урала / И. Е. Артемов // Использование и охрана природ, ресурсов в России. 2006, №1. - С.98-104.

3. Айзенберг Г.З. Коротковолновые антенны. -М.: Связьиздат, 1962. -815 с.

4. Александров Г.Н. Сверхвысокие напряжения. Л.: Энергия, Ленинградское отделение, 1973. - 180 с.

5. Аполлонский С.М. Внешние электромагнитные поля электрооборудования и средства их снижения. СПб.: Безопасность, 2001. - 620 с.

6. Атабеков В.Б., Крюков В.И. Городские электрические сети. М.: Строй-издат, 1987.-384 с.

7. Атаманова И.Г., Левин В.И. Уравнения математической физики. М.: Наука, 1964.-365 с.

8. Баге К., Вилсон Е. Численные методы анализа и метод конечных элементов. М.: Стройиздат, 1982. - 250 с.

9. Базуткин В.В., Кадомская К.П., Колечицкий Е.С. и др. Физико-математические основы техники и электрофизики высоких напряжений. / Под ред. К.П. Кадомской. М.: Энергоатомиздат, 1995. - 416 с.

10. Бартенев О.В. Visual Fortran: новые возможности. М.: Диалог МИФИ, 1999.-232 с.

11. Бартенев О.В. Современный Fortran. M.: МИФИ, 2000. - 446 с.

12. Башкин В.Н. Управление экологическим риском. МГУ Им. М.В. Ломоносова, Науч.-исслед. и проект.-изыскат. инст-т экологии города. - М.: Науч. мир, 2005. - 367 с.

13. Бахвалов Н.С., Жидков Н.П., Кобельков Г.М. Численные методы. -М.: Наука, 2000. 350 с.

14. Беклемишев Д.В. Курс аналитической геометрии и линейной алгебры. -М.: Наука, 1976. -320 с.

15. Берлянт A.M., Тикунов B.C. Геоинформационные системы: Сб. переводных статей М.: Картгеоцентр-Геодезиздат, 1994. - 180 с.

16. Битюкова, В.Р. Экологическая безопасность как фактор качества жизни и социального развития регионов России // Межрегиональный экологический форум «Человек. Экология. Здоровье». Барнаул, 2004. - С. 161-163.

17. Блохин Ю.В., Журавлев З.Н., Ярославский Э.Н. К расчету электростатических полей методом эквивалентных зарядов Электричество, 1980, №2, С.26-31.

18. Богородицкий Н.П., Пасынков В.В., Тареев Б.М. Электротехнические материалы. Ленинград: Энергоатомиздат, 1985. 304 с.

19. Борисов Б.М. К вопросу об оценке состояния здоровья населения в условиях антропогенного загрязнения окружающей среды. // Экология промышленного производства. 1999. №1. С.3-6.

20. Брауде Л.Г. Использование сетчатых моделей для расчета входных сопротивлений самолетных антенн ДКМВ.// Труды НИИР. № 3, 1989. - С.79-82.

21. Бронштейн И.Н., Семендяев К.А. Справочник по математике для инженеров и учащихся втузов. М.: Наука. Гл. ред. физ.-мат. литературы, 1981. - 652 с.

22. Бузов А.Л., Казанский Л.С., Минкин М.А., Юдин В.В. Принципы моделирования антенно-фидерного устройства как сложной пространственной структуры обобщенными LC-цепями. // Физика волновых процессов и радиотехнические системы. 1, № 4, 1998. - С.38-41.

23. Бузов А.Л., Маслов М.Ю. Моделирование электромагнитных полей, возникающих за счёт антенного эффекта технических средств в закрытых помещениях. // Антенны и электродинамика СВЧ. № 7, 2002. - С.9-12.

24. Бузов А.Л., Сподобаев Ю.М. Электромагнитная экология. Основные понятия и нормативная база. М.: Радио и связь, 1999. - 78 с.

25. Бузов А.Л., Сподобаев Ю.М., Филиппов Д.В., Юдин В.В./ Под ред. В.В. Юдина. Электродинамические методы анализа проволочных антенн. М.: Радио и связь, 2000. - 153 с.

26. Бузов А.Л., Сподобаев Ю.М., Юдин В.В. Электромагнитные поля и волны. Термины и определения. Справочное пособие. Самара: СОНИИР, 1999.-70 с.

27. Булгаков Н.Г. Технология регионального контроля природной среды по данным биологического и физико-химического мониторинга. // Диссертация на соискание ученой степени доктора биологических наук. Москва, МГУ, 2003.-270 с.

28. Быстрова А.К. Экология и капиталистический город. М.: Н., 1980.- 173 с.

29. Бялобок С. Регулирование загрязнения атмосферы. // Загрязнение воздуха и жизнь растений. Л., 1988. С500-531.

30. Ветров В.В., Хрупачев А.Г. Метод оценки и прогнозирования влияния вредных техногенных факторов на продолжительность жизни человека // Вестник новых медицинских технологий 1998. 5. №3-4. С. 15-17.

31. Владимиров B.C. Уравнения математической физики. М.: Наука, 1981.-250 с.

32. Волков Е.А. Численные методы. М.: Наука, 1982. - 200 с.

33. Воробейчик Е.Л. Садыков О.Ф., Фарафонтов М.Г. Экологическое нормирование техногенных загрязнений. Екатеринбург: Наука, 1994. 280 с.

34. Вычислительные методы в электродинамике. Под ред. Р. Митры./ Пер. с англ. Под. ред. Э.Л. Бурштейна. М.: Мир, 1977. - 487 с.

35. Градштейн И.С., Рыжик И.М. Таблицы интегралов, сумм, рядов и произведений. М.: Наука, 1971. - 512 с.

36. Галлагер Р. Метод конечных элементов. Основы. М.: Мир, 1984.-428 с.

37. Гахов Ф.Д. Краевые задачи. М.: Наука, 1977. - 640 с.

38. Говорков В.А. Электрические и магнитные поля. М., Госэнергоиз-дат, 1960.-463 с.

39. Гринберг Г.А., Избранные вопросы математической теории электрических и магнитных явлений. М.: Изд-во АН СССР, 1948. - 727 с.

40. Гринберг Г.А., Пименов Ю.В. К вопросу о дифракции электромагнитных волн на бесконечно тонких идеально проводящих экранах.// Журнал теоретической и экспериментальной физики, т. XXVII, вып. 10, 1957. с. 2326- 2339.

41. ГОСТ 12.1.002 84. Электрические поля промышленной частоты. Допустимые уровни напряженности и требования к проведению контроля на рабочих местах.

42. Григорьев Ю.Г., Степанов B.C., Григорьев O.A., Меркулов A.B. Электромагнитная безопасность человека. Справочно-информационное издание. -М.: Российский национальный комитет по защите от неионизирующего излучения, 1999.- 151 с.

43. Гудериан Р. Загрязнение воздушной среды. М.: Мир, 1979. 200 с.

44. Данилушкин А.И., Довбыш В.Н. Проблемы экологии при утилизации взрывателей. // Труды VII Международной конференции «Окружающая среда для нас и будущих поколений». Самара, СамГТУ, 20021, С.58-59.

45. Дезин A.A. Общие вопросы теории граничных задач. М.: Наука, 1980.- 120 с.

46. Демидович Б.П., Марон И.А., Шувалова Э.З. Численные методы анализа./ 3-е, перераб. М.: Наука, 1967.-368с.

47. Довбыш В.Н., Котельников Г.П., Самыкина Л.Н., Федосейкина И.В. Антропогенные загрязнения территории и их влияние на состояние здоровья населения // Труды VII Всероссийского конгресса «Экология и здоровье человека», Самара 2001. -С.90-94.

48. Довбыш В.Н. Идентификация процесса непрерывного индукционного нагрева изделий при утилизации взрывателей // Труды двенадцатой межвузовской конференции «Математическое моделирование и краевые задачи».Ч.2. Самара: СамГТУ, 2002. -С.31-34.

49. Довбыш В.Н. Моделирование и расчет электромагнитных источников тепла при индукционном нагреве осесимметричных тел сложной формы // «Вестник Самарского государственного технического университета».Серия «Технические науки». Вып. 14, 2002. -С. 181-184.

50. Данилушкин А.И., Довбыш В.Н. Проблемы экологии при утилизации взрывателей // Труды седьмой Международной конферении «Окружающая среда для нас и будущих поколений». Самара, СамГТУ, 2002. -С.58-59.

51. Абакумов A.M., Довбыш В.Н., Зимин JI.C. Энергосбережение и экология в электротехнологии // Труды седьмой Международной конферении «Окружающая среда для нас и будущих поколений». Самара, СамГТУ, 2002. -С. 13-14.

52. Довбыш В.H. Проблемы и тенденции состояния природной среды Самарской области // Труды VIII Международного конгресса «Актуальные проблемы экологии человека», Самара 2002. -С.64-68.

53. Довбыш В.Н., Маслов М.Ю., Сподобаев Ю.М. Численный анализ электромагнитного поля силовых трансформаторов в реальных условиях размещения //ЕжегодникРНКЗНИ2003. Сб. трудов.-М.: Изд-во АЛАНА, 2004. -С. 151-163.

54. Довбыш В.Н. Проблемы электромагнитной экологии энергетических систем // Материалы IX Всероссийского конгресса «Экология и здоровье человека». Самара, 2004. -С.89-92.

55. Довбыш В.Н., Маслов М.Ю., Ситникова C.B., Сподобаев Ю.М. Электродинамическая модель трансформатора с сердечником, выполненным из материала с нелинейными свойствами // Инфокоммуникационные технологии, Том 2, №3, 2004. -С. 72-76.

56. Довбыш В.Н., Сподобаев Ю.М. Электромагнитный мониторинг энергетических систем // Материалы V Международной научно-технической конференции «Проблемы техники и технологии телекоммуникаций». Самара, 2004,.-С. 13-14.

57. Довбыш В.Н. Электромагнитные поля' и здоровье человека // Материалы X Конгресса «Экология и здоровье человека», Самара, 2005. -С.88-91.

58. Довбыш В.Н., Сивков B.C. Применение геоинформационных технологий в системах электромагнитного мониторинга // Материалы X Конгресса «Экология и здоровье человека», Самара 2005. -С.91-93.

59. Довбыш В.Н., Маслов М.Ю., Сподобаев Ю.М. Электромагнитный мониторинг энергетических систем // Академия Энергетики. №1 (09), 2006. -С.4-7.

60. Довбыш В.Н. Принципы электромагнитного мониторинга энергетических систем // Материалы XIII Юбилейная Российская научная конференции профессорско-преподавательского состава, научных сотрудников и аспирантов. 30.01 -4.02.2006. Самара, ПГАТИ. -С.150.

61. Довбыш В.Н., Маслов М.Ю., Сарокваша О.Ю. Экологический мониторинг состава почвы в зоне размещения высоковольтной ЛЭП // Вестник Самарского государственного университета. Естественнонаучная серия. 2006. №4(44). -С. 136-148.

62. Довбыш В.Н., Сивков B.C., Сподобаев Ю.М. Визуализация электромагнитной обстановки, создаваемой телекоммуникационными техническими средствами, расположенными на больших территориях // Антенны, №10 (113), 2006. -С.58-62.

63. Довбыш В.Н. Анализ помехового излучения линии электропередач // Радиотехника, №10, 2006. -С.69-71.

64. Довбыш В.Н., Маслов М'.Ю. Анализ излучения коронного разряда высоковольтных линий передачи // Материалы VII Международной научно-технической конференции «Проблемы техники и технологии телекоммуникаций». Самара, 2006. -С.270-272.

65. Довбыш В.Н., Сподобаев Ю.М. Проблемы электромагнитного мониторинга мегаполисов // Материалы VII Международной научно-технической конференции «Проблемы техники и технологии телекоммуникаций». Самара, 2006. -С.272-274.

66. Довбыш В.Н., Сподобаев Ю.М. Уровни и источники электромагнитного загрязнения Самары // XIV Российская научная конференция профессорско-преподавательского состава, научных работников и аспирантов. Самара, ПГАТИ, 2007.-С.110-111.

67. Довбыш В.Н. Оценка экологической обстановки и масштабов природопользования с использованием ГИС // XIV Российская научная конф. профессорско-преподавательского состава, научных работников и аспирантов. Самара, ПГАТИ, 2007. -С. 108.

68. Довбыш В.Н., Семаков Л.М. Расчет электромагнитных полей первичных объектов энергоснабжения области // XIV Российская научная конф. профессорско-преподавательского состава, научных работников и аспирантов. Самара, ПГАТИ, 2007. -С. 110-111.

69. Довбыш В.Н., Сивков B.C. Цифровая электромагнитная модель местности // Инфокоммуникационные технологии, том 5, №1, 2007. -С.85-88.

70. Довбыш В.Н., Сивков B.C. Формирование пространственных данных для геоинформационной системы электромагнитной безопасности // Материалы VIII МНТК «Проблемы техники и технологий телекоммуникаций». Уфа, 2007.-С.191-194.

71. Довбыш В.Н., Сподобаев М.Ю. Современные методы оценки и прогнозирования экологической обстановки // Материалы VIII МНТК «Проблемы техники и технологий телекоммуникаций». Уфа, 2007. -С. 196-199.

72. Довбыш В.Н., Сивков B.C., Сподобаев Ю.М. Визуальное представление электромагнитной обстановки на территории мегаполиса // Труды XIII международной научно-технической конференции «Радиолокация, навигация, связь». Том 2. 2007. -С.2036-2042.

73. Довбыш В.Н. Региональный контроль электромагнитных полей объектов энергетических систем (монография). Офорт, Самара, 2008. 130 с.

74. Довбыш В.Н. Электродинамическое моделирование цепей питания электротранспорта // Известия Самарского научного центра российской академии наук. Специальный выпуск «XIII конгресс Экология и здоровье человека» Том 1. Самара, 2008. -С.141-146.

75. Довбыш В.Н., Маслов М.Ю., Сподобаев Ю.М. Электромагнитная безопасность элементов энергетических систем (монография). «ИПК «Содружество», Самара, 2009. 198 с.

76. Довбыш В.Н. Расчет электромагнитных полей, создаваемых цепями питания электротранспорта // Материалы VIII Международного симпозиума по электромагнитной совместимости и электромагнитной экологии ЭМС-2009. Санкт-Петербург, 2009. -С.27-30.

77. Довбыш В.Н., Маслов М.Ю., Сподобаев Ю.М. Электромагнитная безопасность элементов энергетических систем // Материалы VIII Международного симпозиума по электромагнитной совместимости и электромагнитной экологии ЭМС-2009. Санкт-Петербург, 2009. -С.27-30.

78. Довбыш В.Н., Маслов М.Ю., Сподобаев Ю.М. Энергетические системы: электромагнитный мониторинг и безопасность // Материалы Второго международного экологического конгресса ELPIT-2009. Тольятти, 24-27.09. 2009 г.

79. Довбыш В.Н. Экологический анализ использования природных ресурсов и оценка состояния окружающей среды Самарской области // Известия Самарского научного центра Российской Академии наук. Том 11, №13, 2009. -С.472-476.

80. Довбыш В.Н. Экологический мониторинг элементов региональной энергетической инфраструктуры по фактору электромагнитного излучения // Известия Самарского научного центра Российской Академии наук. Том 11, №13,2009. -С.476-482.

81. Ефимов H.B. Высшая геометрия. М., Физматгиз, 1961. 580 с.

82. Жигальский O.A. Экологическое нормирование антропогнных нагрузок // Материалы III междунар. конф. «Освоение Севера и пробл. рекультивации». Сыктыкар, 1997. С.73-75.

83. Заворыкин В.М. и другие. Численные методы. М.: Просвещение, 1980.-250 с.

84. Захаров Е.В., Пименов Ю.В. Численный анализ дифракции радиоволн. М.: Радио и связь, 1982. - 184 с.

85. Зенкевич О.С. Метод конечных элементов в технике. М.:Мир, 1975-541с.

86. Иванов В.А. Система управления промышленной безопасностью. // Экология и промышленность России. 2005. - №6. - Спец. вып.: Природоохранная деятельность химического предприятия.

87. Иванов В.К. О некорректно поставленных задачах. // Математический сборник. №2 (61), 1963. - С.75-79.

88. Израэль Ю.А. Экология и контроль состояния природной среды. М.: Гидрометеоиздат, 1984. 435 с.

89. Икрамов К.Д. Численные решения матричных уравнений. — М.: Наука, 1984. 300 с.

90. Инженерные расчеты на ЭВМ. Справочное пособие. / Под ред. В.А. Троицкого. JL: Машиностроение, Ленингр. отд., 1979. - 288 с.

91. Иоссель Ю.Я., Качанов Э.С, Струнский М.Г. Расчет электрической емкости. 2-е изд. - Л.: Энергия, Ленингр. отд-ние, 1981. - 282 с.

92. Казаков В.Н. Уравнение коронного разряда. // Цветные металлы. -№6 2003 г. С.44 -47.

93. Кайданов Ф.Г., Перельман Л.С. Экранирующие навесы для защиты персонала подстанции сверхвысокого напряжения. Электрические станции, 1975, №1.- С.75-78.

94. Корягина Е.Е., Коськин O.A. Электрооборудование трамваев и троллейбусов. Учебник для техникумов городского транспорта. М.: изд. Транспорт, 1982. 296 с.

95. Капцов H.A. Коронный разряд. М.: ОГИЗ Гостехиздат, 1947 г. - 256 с.

96. Краскович и др. Численные методы в инженерных исследованиях. -Киев: В/ш, 1986.-320 с.

97. Краснов M.JL, Киселёв А.И., Макаренко Г.И. Интегральные уравнения. М.: Наука, 1968. - 192 с.

98. Крылов В.А., Юченкова Т.В. Защита от электромагнитных излучений. М.: Сов. Радио, 1972. - 130 с.

99. Кубанов В.П., Маслов О.Н., Сподобаев Ю.М. Электромагнитная экспертиза независимость и компетентность. // Телекоммуникационное поле регионов, № 3, 1999. - С. 22-25.

100. Кудрявцев Л.Д. Курс математического анализа. В 3 т., Т. 3, — 2-е изд., переработ, и доп. М.: Высшая школа, 1988. - 352 с.

101. Кузнецов А.Н. Биофизика электромагнитных воздействий (основы дозиметрии). М.: Энергоатомиздат, 1994. - 80 с.

102. Кузнецов А.Н. Биофизика низкочастотных электромагнитных воздействий. Учебное пособие. М.: МФТИ, 1994. - 90 с.

103. Кунцман С.Е. Численные методы. М.: Наука, 1979. - 350 с.

104. Ландау Л.Д., Лифшиц Е.М. Теоретическая физика. Т2. Теория поля. М.: Наука, 1988.-509 с.

105. Левич А.П., Булгаков Н.Г., Максимов В.Н. Теоретические и методические основы технологии регионального контроля природной среды по данным экологического мониторинга. М.: Наука, 2004. - 250 с.

106. Мак-Кракен Д., Дорн У. Численные методы и программирование на ФОРТРАНе. М.: Мир, 1969. - 583 с.

107. Марков Г.Т., Чаплин А.Ф. Возбуждение электромагнитных волн. / 2-е изд., перераб. и доп. М.: Радио и связь, 1983. - 296 с.

108. Маслов М.Ю. Численный анализ электромагнитной обстановки в офисном помещении. // Вестник СОНИИР №1, 2004. С. 162 - 168.

109. Маслов О.Н. Вероятностное моделирование и нормирование уровней электромагнитного фона.// Труды Межд. Ак. Связи, № 2(6), 1998. С. 12-16.

110. Маслов О.Н. Вероятностное моделирование последствий непороговых электромагнитных воздействий. // Физика волновых процессов и радиотехнические системы, Т.1, № 4, 1998. С.30-34.

111. Маслов О.Н. Электромагнитная безопасность радиоэлектронных средств. Серия изданий «Связь и бизнес». М.: МЦНТИ, 2000. - 82 с.

112. Макеева Е.Д. Становление и развитие природоохранного и экологического движения в Самарской обл. // Дисс. на соиск. уч. ст. канд. ист. наук. С амарский государственный педагогический университет. — Самара, 2002—238 с.

113. Мамихин C.B. Компьютеризация экологических исследований. //Вест. Российского университета дружбы народов. Экология и безопасность жизнедеятельности. 2000. №4. С. 150-157.

114. Маткад 6.0. Финансовые, инженерные и научные расчёты. / Пер с англ. MathSoft Inc. 101 Main Street. Cambridge. Massachusetts, 02142USA- 698 c.I

115. Методы граничных элементов. Перевод с англ. / Бреббия К., Телес Ж., Вроубел Л. М., Мир, 1987. 524 с.

116. Миролюбов H.H., Костенко М.В., Левинштейн М.Л., Тиходеев H.H. Методы расчета электростатических полей. М.: Высшая школа, 1963. 415 с.

117. Моисеев В. Н. Способ вычисления электрического поля на заряженной поверхности, обладающей осевой симметрией. Труды МЭИ, 1979, вып. 412.-С. 97-100.

118. Молчанов А.П., Занадворнов П.Н. Курс электротехники и радиотехники. М.: Наука, глав. ред. физ.-мат. лит., 1969. - 480 с.

119. Назаров В.Е., Рунов A.B., Подининогин В.Е. Численное решение задач об основных характеристиках и параметрах сложных проволочных антенн // Радиотехника и электроника, Вып. 6. — Минск: Вышейшая школа, 1976. С.153-158.

120. Пименов Ю.В., Вольман В.И.', Муравцов А.Д. Техническая электродинамика. М.: Радио и связь, 2000. - 450 с.

121. Писсанецки С. Технология разреженных матриц. / Пер. с англ. Ик-рамова Х.Д. и Капорина И.Е., под ред. Икрамова Х.Д. М., Мир, 1988. 411 с.

122. Попков В. И. Электропередачи сверхвысокого напряжения, в кн.: Наука и человечество, (т. 6), М., 1967 г. 230 с.

123. Правила технической эксплуатации электроустановок потребителей. Раздел 1.-М.: НЦЭНАС, 2003.-304 с.

124. Правила технической эксплуатации электроустановок потребителей. Раздел 2. Передача электроэнергии. М.: НЦ ЭНАС, 2003. - 152 с.

125. Правила устройства электроустановок. НЦ ЭНАС, 2003. - 299 с.

126. Предварительные итоги деятельности Федеральной службы надзора в сфере природопользования в 2005 году. // Использование и охрана природ, ресурсов в России. 2006. - № 1. - С.42-44.

127. Прогулова Т.Б. Геоинформационные системы при подготовке специалистов в области недропользования. // Дисс. на соиск. уч. ст. к.т. н. / Московский государственный геолого-разведочный университет. М., 2003. - 100 с.

128. Разевиг Д.В., Соколова М.В. Расчет начальных и разрядных напряжений газовых промежутков. М.: Энергия, 1977. - 199 с.

129. Резвых К. А. Расчет электростатических полей в аппаратуре высокого напряжения. -М.: Энергия, 1967. 121 с.

130. Ряполов С.И. Обобщённый метод численного решения задач Коши. / Под ред. Баринова. М-во обороны, 1975. - 80 с.

131. Самарский A.A. Теория разностных схем. М., Наука, 1977. - 656 с.

132. Abou-Seada M. S., Nasser E. Digital Computer Calculation of the Electric Potential of a Rod Gap. Proc. IEEE, 1968, vol. 56, № 5. - P.813-820.

133. Antenna Modeling Program Engineering Manual, MB Associates Report No. MB-R-74/62, 1994. - 85 p.

134. Antenna Modeling Program Supplementary Computer Program Manual (AMP2), MB Associates Report No. MB-R-75/4, 1975. - 72 p.

135. Bulgakov N.G. Determination analysis as a method for diagnostics of ecosystem condition. // Environmetrics. 2003. P. 122-135.

136. Calow P. Ecological risk assessment: Risk for what? How do we decide? //Ecotoxicol. and Environ. Safety. 1998. P.15-18.

137. Comparative Study of three Methods for Computing Electric Fields / Beasley B. A., Pickles J. H., D'Amico G. et. al. -Proc. IEE, 1979, vol. 126, №1. P.126-134.

138. David O. Carpenter. Sinerik Ayrapetyan «Biological Effects of Electric and Magnetic Fields» Academic press, 1994. 369 p.

139. Electromagnetic Fields (300 Hz to 300 GHz). Environmental health criteria; 137, Geneva: WHO, 1993. P.290.

140. Handbook of Biological Effects of Electromagnetic Fields edited by Charles Polk, Elliot Postow. New York, London, Tokyo: CRC Press, 1996. - 618 p.

141. Hannon B. Ecosystem flow analysis //Can. Bull. Fish. Aq. Sc. 1985.-P.97-118.

142. Harrington R.F. Field Computation by Moment Method. Macmillan, New York, 1968.- 150 p.

143. Hoorfar A., Jamnejad V. Electromagnetic Modeling and Analysis of Wireless Communication Antennas.// IEEE Microwave Magazine. Mar. 2003/ P.51-67.

144. Numerical Electromagnetic Code (NEC-1) Part I: NEC Program Description Theory, Lawrence Livermore Laboratory, 1977. 85 p.

145. Lytle R.J., Lager D.L. Numerical Evaluation of Sommerfeld Integrals, UCRL-51688, Lawrence Livermore Laboratory, CA, October 23, 1974. 170 p.

146. Maximov V.N., Bulgakov N.G., Levich A.P. Quantitative methods of ecological control: Diagnostics, standardization, and prediction // Environmental indices: Systems Analysis Approach. London: EOLSS Publishers, 1999. 363 p.

147. Poggio A.J. and Miller E.K. Integral Equation Solutions of Three-Dimensional Scattering Problems, Chapt. IV in Computer Techniques for Electromagnetics, edited by R.Mittra, Pergamon Press, New York, 1973. P. 125 - 132.

148. Poggio A.J. and Adams R.W. Approximations for Terms Related to the Kernel in Thin-Wire Integral Equation, UCRL-51985, Lawrence Livermore Laboratory, CA, December 19, 1975. 52 p.

149. Shannon, C.E., Weaver W. The Mathematical Theory of Communication. Urbana, Illinois: Univ. of Illinois Press, 1949 350 p.

150. Silvester P., Chari M. Finite element solution of saturate magnetic field problems. IEEE Trans. Power Appar. Syst., 1970, vol. 89, #7. p. 1642 - 1651.

151. Thide B. Electromagnetic Field Theory. Swedish Institute of Space Physics; Department of Astronomy and Space Physics Uppsala University, Sweden; School of Mathematics and Systems Engineering Vaxjo University, Sweden. 2002.-191 p.

152. Urban ecology/ / edited: J. Breuste, H. Feldmann, O. Uhlmann. Berlin et al.: Springer, 1998. - XVIII. - 714 p.

153. William Ralpf Bennett, jr. Health and Low-Frequency Electromagnetic Fields. Yale University, 1994. 185 p.