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

УДК 621.316.99

Н.В. Коровкин, С.Л. Шишигин

РАСЧЕТНЫЕ МЕТОДЫ В ТЕОРИИ ЗАЗЕМЛЕНИЯ

Заземляющие устройства (ЗУ) рабочего, защитного, молниезащитного и помехозащитного заземления играют важную роль в обеспечении надежности, электробезопасности и электромагнитной совместимости объектов электроэнергетики. Поэтому к ним предъявляются высокие требования по точности и адекватности как расчетов на этапе проектирования, так и измерений на этапе эксплуатации [1]. ЗУ служит для растекания токов короткого замыкания (КЗ), молнии, импульсных помех и состоит из заземлителей — проводников в земле и в воздухе, включая систему молниезащиты, металл окон-струкции зданий и оборудования, грозозащитные тросы воздушных линий, экраны кабельных линий и т. д. Как правило, геоэлектрическая структура земли — многослойная, горизонтально слоистая, включающая локальные неоднородности произвольной формы. Расчету подлежит сопротивление ЗУ, распределение его токов, напряжение прикосновения, шаговое напряжение, распределение потенциала, напряженности электрического и магнитного поля. Таким образом, имеем сложную, теоретически и практически важную задачу. Показать современный уровень ее решения — цель настоящей работы.

Выбор математической модели ЗУ. Наиболее точное решение задачи дает модель заземлителя, основанная, например, на решении уравнений Максвелла с использованием программы FDTD. Но дискретизация расчетного объема потребует огромного числа узлов, поскольку линейные размеры ЗУ на 3—5 порядков превышают сечение проводников, которые определяют шаг сетки. Приходится также искусственно замкнуть расчетную область для задания граничных условий. Цепная модель ЗУ, которая может быть использована в программе ЕМТР, существенно проще, но адекватно учесть все электромагнитные связи между элементами, тем более — структуру земли и экранирующие эффекты, здесь невозможно. Кроме того, результатом расчета будут интегральные параметры — напряжения и токи элементов, но не дифференциальные — напряженности электромагнитного поля, необходимые в задачах электромагнитной совместимости (ЭМС).

Задачи расчета ЗУ по постановке, методам и требуемым результатам являются цепно-по-левыми, а для их решения требуются две взаимосвязанные модели — полевая и цепная (рис. 1). Полевая модель позволяет рассчитать

б)

°у1] Су1

"ист

ту1

МУ Оуу Суу

1 "I

-у11

2

¿1 + Ыи чт

у2

у22

щ I

°у3] Су33

3 п + 1

ь2+м22

ту3

у22

'■Су33 Ц°у33

Рис. 1. Полевая (а) и цепная (б) модели стержневого заземлителя

2

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

Цепно-полевая модель ЗУ. Заземлитель дробится на элементы длиной I <Х/10 = 7105 р//

(где X—длина электромагнитной волны частотой /в проводящей среде с удельным сопротивлением р), что позволяет проводить расчет электромагнитных параметров элементов в статическом приближении. Предполагается, что источники поля — электрические заряды, а также стекающие и продольные токи сосредоточены на осях проводников круглого сечения. Внутреннее активное сопротивление и индуктивность элементов определяются с учетом поверхностного эффекта и описываются диагональными матрицами г и К Электрические и магнитные связи между элементами описываются полностью заполненными матрицами собственных и взаимных проводимостей G, емкостей С, индуктив-ностей M. Матрица проводимостей растекания тока получается обращением матрицы сопротивлений — G = R-1.

Элемент Яу матрицы R, имеющий смысл взаимного сопротивления в однородной среде с удельным сопротивлением р, определяется как отношение потенциала в средней точке 1-го элемента к току у-го элемента:

Яу = Я( р, д, I) =

4п /

71п

(д +1 - р)1 + \д +1 - р||/|

(д - р)1 + \д - р\\1\

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

Собственное сопротивление стержня диаметром ё определяется как

Яа =-р-1 11 2п 111

1п

2 +ё2

\ I \> ё.

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

мене р на 1/е) сопротивлениям элементов Яу в однородной земле.

Взаимная индуктивность элементов равна Му = Чу/у где Чу = /4 4 « Ау\ — потокосце-

пление 1-го стержня с током у-го стержня, которое определяется интегрированием векторного потенциала по длине 1-го элемента (стандарт МЭК 60050—121). Используя обозначения рис. 1, получим

му =

4п |/,.|

1п

(д+1у - р)1у + д+у - р Ы

(д - р)1у + д- рЩ

Для перехода к цепной схеме ЗУ (рис. 2, б) выполняется преобразование матриц G и С, определенных в средних точках элементов (рис. 2, а), в узловые матрицы Gy и СУ (рис. 2, б) из условия неизменности стекающего тока элементов. Топология продольных ветвей схемы (рис. 2, б) описывается стандартной матрицей соединений А и вводимой матрицей В (Ьу = = \а1у\/2). Тогда [2] искомые узловые матрицы будут определяться так: Gy = В^-ВТ, СУ = В-С-ВТ.

Цепная модель ЗУ позволяет определить потенциалы и токи (стекающие и продольные) элементов методами теории цепей. При гармонических воздействиях используется метод узловых потенциалов. Методы расчета переходных процессов, возникающих при импульсных воздействиях, рассмотрены далее. По найденным стекающим токам стержней в модели (рис. 2, а) определяется распределение потенциала и напряженности электрического поля, по найденным продольным токам — напряженность магнитного поля.

р = Р1

Яу

Му

\

1 = у

д = ду т = ту

Рис. 2. Определение взаимных параметров у-х элементов с координатами д = ду = (Хд, Уд, гд)Т,

т = ту = (Хт, Ут, *т) Т> 1 = 1у = т-д,

р = р( = (хр, Ур, ?-р)Т

а

у

Методы расчета импульсных процессов ЗУ.

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

Частотный метод (ЧМ) с искусственной периодизацией сигнала. Ограничим длительность наблюдения за переходным процессом величиной 'пп = (3—10) где Тх — длительность переднего фронта волны. Продолжим импульс тока /(') вне расчетного интервала (' > 'пп) этой же функцией, но с обращением знака (рис. 3). Получим периодическую функцию /(') с периодом 2'пп, обладающую нечетной симметрией:

VСпи) -1 (' - 'пп), 'пп < ' < 2^пп.

Это позволяет исключить из гармонического ряда четные гармоники:

/(')=А + X Ак сов(к+Фк), ®1 = П'пп.

к=1,3,5...

Расчет коэффициентов ряда проводится с помощью быстрого преобразования Фурье (БПФ).

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

предлагается метод дискретных комплексных схем. Индуктивность моделируется сопротивлением ..Ь, конденсатор — проводимостью .С, где

^ = (2 + л/2/)/И , И — длина шага по времени. Начальные условия п-го шага учитываются источниками Еп = Ь/п для катушки и 1П = Сип для конденсатора. Проводим расчет комплексной схемы. Переход от изображения к оригиналу проводится по формуле / (') = Re((5V2/-2)¥(.))/И.

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

Учет частотнозависимых сопротивлений в МДС. Обычно эта задача решается с использованием эквивалентных схем с близкими частотными характеристиками, но можно предложить более эффективный способ [2].

Пусть задано сопротивление г(/ю) частотно-зависимого элемента или его операторный аналог Z(s). Проинтегрируем Z(.y) в пространстве изображений и, перейдя к оригиналу, получим переходное сопротивление г(') = Ь-1^(.)/,у|, связывающее между собой напряжение и ток в виде интеграла Дюамеля:

(

и(') = 2(0/(0) + ¡г(' - х)/'(х) ёх, /(0) = 0. (*)

о

Рис. 3. Импульс тока (1) при единичной амплитуде импульса, периодизация (2), АЧХ (3)

Дискретная форма записи этого интеграла на сетке с узлами 7п = (п— 1)к, п = 1, ...(И + 1) при кусочно-постоянной аппроксимации производной тока дает

¿п+1

Я,

ип+1 = Е

¿т+1 ¿т

т=1

к

т+1

} z(tn+1 - х) ёх =

= Е ('т+1 - ) Я

т=1

где (с учетом обозначения к = п — т + 1 и подстановки у = ^+1 — х) дискретное переходное сопротивление к-го интервала равно

1 кк

Як = - | г(у)ёу, к = 1,...N. к

(к-1) к

1 пк

Яп = - | г(у) ёу =

(п-1) к

/ ТЙЛ

ип+1

Интегрируя в пространстве изображений, вместо z(t) имеем ¿1(7) = Х-12(,у)/,у2]. Тогда Як = [г1(кк) - г1(кк - к)]/к, к = 1,... N, где ^(7) —

реакция цепи на воздействие тока единичной ступеньки ¿(7) = 7. Выделяя первое слагаемое ип+1, получим соотношение

п

Мп+1 = Я1^п+1 - Е (Яп-т+1 - Яп-т+2) ¿т =

т=2

= Я1гп+1 - Еп ,

которому соответствует схема (рис. 4).

Таким образом, сопротивление г(/ю) полностью описывается дискретными переходными сопротивлениями в шаговых алгоритмах.

В качестве примера найдем дискретное переходное сопротивление стального стержня с операторным сопротивлением г (я) = (// 2ла)^/ / у, где / — длина, а — радиус стержня. Переходное сопротивление стержня ¿(7) = г^г (*)/ ^ ] =

= (// 2ла)Л/ц7уП7. Тогда дискретные переходные сопротивления

- 1 (уШ-4п-\), п = 1,...N а пу]пк

зависят лишь от номера п шага.

Учет нелинейных элементов. К ним относится собственная проводимость элемента ЗУ

Рис. 4. Дискретная модель сопротивления ¿(ую)

с учетом искрообразования в земле (в работе используется модель Е.Я. Рябковой) и с учетом намагничивания железа (используется универсальная кривая Л.Р. Неймана). Нелинейность не вносит принципиальных изменений в работу шаговых алгоритмов — параметры нелинейных элементов при расчете переходных процессов принимаются кусочно-постоянными и равными значению в начале каждого шага. В большинстве случаев насыщение стальных стержней зазем-лителя мало сказывается на результатах расчета, поэтому рекомендуется выбирать постоянное значение ^/у стержней по универсальной кривой Л.Р. Неймана, а для оценки погрешности проводить расчет дважды — с минимальным и максимальным значениями ^/у.

Сопротивление заземлителя при импульсных воздействиях. Допустим, что в результате расчета или измерений известны временные функции напряжения и(7) и тока ¿(7) на входе зазем-лителя (рассматриваемого как пассивный двухполюсник). Требуется найти топологию, параметры и входное сопротивление двухполюсника.

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

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

т

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

Подойдем к определению импульсного сопротивления заземлителя, опираясь на методы теоретической электротехники. Заменим импульсы тока и напряжения эквивалентными синусоидами с частотой / = 0,25/ Тх (эквивалентная частота) при неизменной амплитуде и ограничимся расчетным интервалом К7\. Тогда комплексное сопротивление заземлителя 2 = и/1 = \2\е](? , а модуль сопротивления

1Д7х)| = ит/1т полностью совпадает с импульсным сопротивлением заземлителя ги = |Д7\)|. Таким образом, импульсное сопротивление имеет смысл модуля сопротивления на эквивалентной частоте импульса, поэтому является функцией длительности фронта волны 2и(7\). Расчет зависимости zи(T]) можно проводить на импульсе с линейным фронтом, либо по частотной характеристике модуля сопротивления \2()\ при / = 0,25/7\. Таким образом, понятие импульсного сопротивления заземлителя ги( 7\) обосновано с помощью метод эквивалентных синусоид.

Этот приближенный параметр, зависящий от фронта импульса тока, не рекомендуется применять при t > Тх.

Мгновенное сопротивление заземлителя г(/) = и(И)/1(И) — следующий формальный параметр, призванный показать изменение сопротивления заземлителя во времени и используемый рядом исследователей для подбора схем замещения импульсных заземлителей. Введение этого параметра не имеет под собой теоретической основы для ^ХС-цепей, поскольку не учитывает запасы энергии в электрической цепи, поэтому ограничимся лишь констатацией факта его существования.

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

Для нахождения 1(1) воспользуемся дискретной формой записи интеграла Дюамеля на временной сетке с узлами 1п = (п— 1)к, п = 1, ... N + 1), где N — число интервалов длиной к.

Примем, что производная тока в пределах каж-

г \

дого интервала постоянна

= const

t =tn

на-

J

чальные условия нулевые. Тогда

Un+1 = Z

n ; —; m+i lm+1 lm

m=1

J z(tn+1 - x) dx =

X ('т+1 'т ) 2и-т+1 2п'2 + Х ('т+1 'т ) 2п-т+1' т=1 т=2

откуда получим рекуррентную формулу для определения дискретных значений переходного сопротивления 2п = 2(1п):

21 = и2/

п

2п = ип+1 - X 2п-т+1('т+1 - 'т ) / г2;

т=2 _

п = 2,... N; /1 = и1 = 0.

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

Программная реализация. Рассмотренные модели и методы расчета ЗУ реализованы в программе ZYM как приложение к AutoCAD. Использование AutoCAD в качестве геометрического процессора (требование проектировщиков) позволяет моделировать сложные объекты (рис. 5) и представлять результаты расчета в виде 2D-, BD-графиков с анимацией динамики. Сеточная структура стен и металлоконструкций ЗРУ (см. рис. 5) создает объемный электромагнитный экран, существенно снижающий уровень напряженности магнитного поля в местах расположения микропроцессорных устройств.

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

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

1. Дьяков, А.Ф. Электромагнитная совместимость в электроэнергетике и электротехнике [Текст] / А.Ф. Дьяков, Б.К. Максимов, Р.К. Борисов [и др.]; под ред. А.Ф. Дьякова.— М.: Энергоатомиздат, 2003.— 768 с.

2. Шишигин, С.Л. Математические модели и методы расчета заземляющих устройств [Текст] / С.Л. Шишигин // Электричество.— 2010. № 1.— С. 16— 23.

3. Демирчян, К.С. Теоретические основы электротехники [Текст]: Учебник для вузов / К.С. Демирчян, Л.Р. Нейман, Н.В. Коровкин.— 5-е изд. Том 1, 2. — СПб.: Питер, 2009.

4. Влах, И. Машинные методы анализа и проектирования электронных схем [Текст] / И. Влах, К. Сингхал / Пер. с англ.— М.: Радио и связь, 1988.— 560 с.

m