Исследование распределённых элементов систем молниезащит на основе вычислительных экспериментов Текст научной статьи по специальности «Электротехника, электронная техника, информационные технологии»

УДК 621.398 З622

ИССЛЕДОВАНИЕ РАСПРЕДЕЛЁННЫХ ЭЛЕМЕНТОВ СИСТЕМ МОЛНИЕЗАЩИТ НА ОСНОВЕ ВЫЧИСЛИТЕЛЬНЫХ ЭКСПЕРИМЕНТОВ

© 2010 А Н. Потапенко1, А.И. Штифанов2, Т.А. Потапенко1

1 Белгородский государственный технологический университет им.В.Г. Шухова 2 Белгородский государственный университет

Поступила в редакцию 22.11.2010

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

Ключевые слова: распределённые молниеотводы, математическая модель, вычислительный эксперимент

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

Потапенко Анатолий Николаевич, кандидат технических наук, профессор кафедры электротехники и автоматики. E-mail: [email protected] Штифанов Андрей Иванович, кандидат технических наук, доцент кафедры информатики и вычислительной техники. E-mail: [email protected] Потапенко Татьяна Анатольевна, аспирантка. Email: [email protected]

применяют распределённые стержневые молниеотводы (СМ) [1] или распределённые тросовые молниеотводы [2].

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

Схема модели подстанции в виде территории В0 с размерами /1х/2 и с двенадцатью плотно расположенными объектами х11 - х34 при их максимальной высоте И с расположением СМ (молниеотводы отмечены (*)) показана на рис. 1. При этом могут быть использованы СМ с расположенными в их верхней части изолированными ИЭУ, основанными на новых физических принципах их действия [3], а именно, с возможностью возбуждения ЭДС в ИЭУ при grad у. Схема для проведения вычислительных экспериментов с расположением СМ для защиты объектов (при разрезе вертикальной плоскостью АБ) показана на рис. 2 с возможностью определения поля распределения потенциала и электрической функции потока напряженности поля исследуемых объектов.

Постановка задачи. На основе анализа известных методов решения внешних краевых задач принимается метод инверсии для численного решения (МИЧР) этих краевых задач с применением уравнения эллиптического типа

С учётом схемы на рис. 2 в вертикальной плоскости области В\(х,£) определяется электрическое поле потенциала у относительно проводящей поверхности земли на основе уравнения Лапласа с учётом подхода [4]:

д V д V

~дхГ ~д1Г

0

; (х^) е Вх(х,7). (1)

Граничные условия для этой задачи следующие:

- на проводящей границе А2 (нижняя часть области В^х^) п. 2 - п. 28):

Рис. 1. Схема расположения СМ для защиты объектов

у(ад)=0, (х^) е А2;

(2)

Рис. 2. Схема для вычислительных экспериментов с учётом СМ - на границе полубесконечной области Ада: уИ=0; (3)

- на границе А1 (верхняя часть исследуемой области В^х^) от п.1 до п.29):

у(х^„) = уа, (х,г„) е А1;

(4)

где В1(х,2) - расчетная область, ограниченная А1, А2, Ада, причем - условная граница на бесконечности слева и справа от исследуемой

области. Для выполнения расчетов в области В1(х,£) введены некоторые внутренние границы А3, А4 (от поз. 1 до поз. 2 и соответственно от поз. 28 до поз. 29, т.е. дополнительно введённые границы для решения задачи с использованием МИЧР для внешних краевых задач [4]), необходимые для выделения некоторой условно ограниченной области ВА1(х,£). Принимаем, что часть выделенной области в виде ВА1(х,£), входящей в состав В1(х,£), является основной, так как в ней задаются границы

исследуемых объектов от п. 1 до п. 29. В рассматриваемом случае ВА1(х,£) имеет форму прямоугольника. Для реализации метода МИЧР необходима дополнительная область В А1(х,г) с учетом вида основной области ВА1(х^) для того, чтобы часть области В^х^), находящаяся между границами А3, А4 и Ада, отобразилась на область В А1(х^). При этом в центре области В А1(х,г), присваивается значение, совпадающее со значением потенциала у на бесконечности. При этом считаем, что области ВА1(х^) и В А1(х^) соприкасаются по границам А3, А4 [4]. В области ВА1(х^) параметры следующие: к - максимальная высота объекта; Н1 - высота от поверхности земли до верхней границы заряженного облака; Нт - высота СМ, рассчитанная в соотношении не менее 0,2к [5]. Вычислительные эксперименты выполнялись на основе метода конечных разностей с применением универсальной прямоугольной равномерной сетки с учётом особенностей исследуемых задач, которые проводились в

прямоугольной области на базе ускоренного метода Либмана [6].

При исследованиях распределенных СМ систем молниезащит определяются поля потенциала у и эквипотенциальных поверхностей в виде уравнений y,=const, а также поверхности равных напряженностей электростатического поля в виде уравнений Et = const. Расчёты ведутся в безразмерном виде, а именно, потенциалы yai =уг/уА; где уА, - базовые значения потенциалов на границах исследуемой области, а компоненты напряженности поля в виде E*x= - ду */dx*; E*y= - д у */ду*. Аналогично записываются безразмерные параметры для других величин. При вычислительных экспериментах для СМ за основу принимается исходное внешнее однородное поля с напряженностью E0=const как в [7]. Типовые результаты расчета согласно уравнениям (1) - (4) для схемы с СМ при разрезе в вертикальной плоскости АБ (рис. 2), проходящей через объекты подстанции (х21 - х24), показаны на рис. 3 в виде уравнений Et = const.

Рис. 3. Результаты расчета в виде уравнений Et = const при H*m0 = Hm/h=1,22

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

поля в виде уравнений Е^со^ с учётом того, что СМ имели относительную величину Н т=1,66, показаны на рис. 4.

Рис. 4. Результаты расчета в виде уравнений E,=const при H m1=1,66

Сравнительный анализ полученных результатов, представленных на рис. 3 и рис. 4 показал следующее:

1. Максимальная напряжённость Етэ на верхних концах распределенных СМ при Н то=1,22 с изменениями относительных величин от 938,36 до 1077,13 увеличивается при Нт1=1,66 с изменениями относительных величин от 1404,96 до 1438,62, т.е. при увеличении высоты молниеотводов Н т в 1,36 раза в среднем пропорционально увеличивается и максимальная напряжённость поля Етэ* с 1007,74 до 1420,87, т.е. в 1,4 раза.

2. Для технологических объектов подстанции при увеличении высоты СМ от Н*т0 до Н т\ наблюдается следующее:

- с высотой объектов меньшей к и с увеличением Н т напряженность поля Етц в среднем уменьшается с 168,30 до 134,87, т.е. 1,25 раза;

- с высотой объектов равной к и с увеличением Н т напряженность поля Ет}2 в среднем уменьшается с 557,40 до 382,13, т.е. в 1,46 раз.

Следует отметить, что при вычислительных экспериментах с учётом относительного шага ЛЕ «33,66 при расчётах распределения поверхностей равных напряженностей электрического поля было установлено, что максимальная напряжённость поля Етэ находится на верхних концах крайних СМ систем мол-ниезащит (см. рис. 1) в сравнении центральным СМ

Вычислительные эксперименты (рис. 5) по определению напряженности поля для одного СМ с расположенным на нём изолированного ИЭУ, основанного на новом физическом принципе действия, а именно, с возможностью возбуждения ЭДС в сравнении для СМ без ИЭУ показали возможность увеличения напряженности поля в пределах до 20% [8] и это частично совпадает с полученными результатами в этой работе.

Е

21-1

о-1—.—,—*—,—.—~ : т •■—. - ■!— I

1,00 1,05 1,10 1.15 1,20 1,25 1,30

Рис. 5. Распределение Е* = flHL*) по вертикали: 1 - СМ без ИЭУ; 2 - СМ с ИЭУ

Сравнительный анализ на основе полевых испытаний, выполненных сотрудниками французской компании INDELEC, для устройств типа Prevectron 2 (СМ с ИЭУ) и одиночного стержня показал преимущества молниеотводов этого типа [3]. Однако в литературных источниках идут весьма острые дискуссии относительно эффективности СМ с применением ИЭУ (молниеотвод типа Prevectron) [9]. Выводы:

1. Исследования, связанные с определением напряженности поля СМ с изолированными ИЭУ, расположенными в их верхней части, с учётом схемы на рис. 1, показали некоторое совпадение результатов с данными [8]. Это связано с тем, что на верхних концах распределенных СМ существует значительный диапазон изменения относительных величин максимальной напряженности поля Етэ , который зависит от места расположения СМ по отношению к объектам подстанции, а это влияет на определение эффективности молниеотводов с применением ИЭУ.

2. Из анализа полученных результатов для распределённых СМ на территории подстанции следует, что при увеличении высоты молниеотводов H т в 1,36 раза максимальная напряженность поля Етэ на верхних концах распределенных СМ систем молниезащит увеличивается, но при этом на объектах подстанции в их верхних угловых зонах напряженность поля Emi снижается. Следовательно, чем больше относительная высота СМ по отношению объектам, тем больше уменьшение Emi на этих объектах, а это важно для проектирования распределённых элементов систем молниезащит.

СПИСОК ЛИТЕРАТУРЫ:

1. Юриков, П.А. Защита электростанций и подстанций 3-500 кВ от прямых ударов молнии / П.А. Юриков. - М.: Энергоиздат, 1982. 88 с.

2. http:// www.siemens.com/hvds (High Voltage Direct Current Transmission - Proven Technology for Power Exchange).

3. http:// www.prevectron.ru.

4. Дыльков, М.И. Метод инверсии для численного решения внешних краевых задач для уравнений эллиптического типа: дис. канд. физ.-мат. наук: 05.13.18 / М.И. Дыльков. - Белгород, 2004. 140 с.

5. Базелян, Э.М. Физика молнии и молниезащиты / Э.М. Базелян. - М.: Физматлит, 2001. 320 с.

6. Сипайлов, Г.А. Электрические машины (специальный курс) / Г.А. Сипайлов, Е.В. Кононенко, К.А. Хорьков. - М.: Высш. шк., 1987. 287 с.

7. Резинкина, М.М. Расчет трехмерных электрических полей в системах, содержащих тонкие проволоки // Электричество. 2005. №1. С. 44-49.

8. Потапенко, А.Н. Численное моделирование электрических полей в системах «электрод - поверхность земли» для элементов молниезащит / А.Н. Потапенко, Е.А. Канунникова, М.И. Дыль-

ков // Известия вузов. Проблемы энергетики. 2008. №11-12. С. 72-78. 9. Базелян, Э.М. Молниезащита высоких сооружений // Известия РАН. Энергетика. 2005. №3. С. 55-74.

RESEARCH OF THE ARRANGED ELEMENTS IN LIGHTNING PROTECTION SYSTEMS ON THE BASIS OF COMPUTING EXPERIMENTS

© 2010 A.N. Potapenko1, A.I. Shtifanov2, T.A. Potapenko1

1 Belgorod State Technological University named after V.G. Shuhov 2 Belgorod State University

Efficiency of the arranged elements of lightning protection systems in the form of lightning rods, including in view the new physical principles of operation, for substations of electric power systems on the basis of mathematical model by means of computing experiments is examined. Results of researches are presented.

Key words: arranged lightning rods, mathematical model, computing experiment

Anatoliy Potapenko, Candidate of Technical Sciences, Professor at the Department of Electrical Engineering and Automation. E-mail: [email protected]

Andrey Shtifanov, Candidate of Technical Sciences, Associate Professor at the Department of Computer Science and Computing Machinery. E-mail: [email protected]

Tatiana Potapenko, Post-graduate Student. E-mail: [email protected]