Оценка параметров наведенных токов, протекающих в заземлителях при грозовых разрядах

УДК 621.316.925

ОЦЕНКА ПАРАМЕТРОВ НАВЕДЕННЫХ ТОКОВ, ПРОТЕКАЮЩИХ В ЗАЗЕМЛИТЕЛЯХ ПРИ ГРОЗОВЫХ РАЗРЯДАХ

СЛЫШАЛОВ В.К., д-р техн. наук, СОРОКИН А.Ф., СЛЫШАЛОВ А.В., КИСЕЛЕВА Ю.А., кандидаты техн. наук, КАНДАЛОВ Ю.В., асп.

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

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

EVALUATING INDUCED CURRENT PARAMETERS IN EARTH LEADS DURING LIGHTNING DISCHARGES

V.K. SLYSHALOV, Doctor of Engineering, A.F. SOROKIN, A.V. SLYSHALOV,

Yu.A. KISELEVa, Candidates of Engineering, Yu.V. KANDALOV, Post-Gradute Student

The authors consider methods of evaluating induced current impulse parameters when lightning strikes the ground near an overhead power transmission line. The latter is presented by a system of wires and grounded cables parallel to each other and the ground surface. The authors determine the current values in the earth lead which connects the groundwire of the power transmission line and cables.

Key words: impulse current, lightning, mathematic model, grounded object.

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

Решение осуществляется при следующих допущениях: лидерный канал считается заряженным равномерно; ток молнии является заданной функцией времени /м() скорость перемещения фронта нейтрализации зарядов лидерного канала пропорциональна току молнии; для системы заземления выполняется неравенство Тз<<Тм, Тз, ТМ - характеристические значения расчетных интервалов времени для заземлителя и грозового разряда соответственно (рис. 1).

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

2Ф < 2 < Л0 канал заряжен при прохождении лидера молнии, поэтому тк = тп, ток на этом участке полагаем равным нулю. В пределах первого участка протекает ток главного разряда /м@), замыкающийся в окружающем канал пространстве током смещения (рис. 1, пунктирные линии) и током проводимости в грунте.

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

ных, необходимых для его проведения, используются, наряду с геометрией системы «канал разряда -провода - тросы» (рис. 1), физические характеристики молнии: ток iM(t) и скорость движения фронта нейтрализации Vo(t), причем допустимо принять

iм

Для построения расчетной математической модели представляем нейтрализацию зарядов лидера как процесс заполнения канала зарядом с плотностью тс= -тп. Заряд вдвигается в канал лидера со скоростью Vo(t).

Количественное соотношение между наведенным током и движущимся зарядом dq = -тЛбг (рис. 1) позволяет установить теорема Шокли-Рамо [1-3]. Для выбранного направления iM(t) , т.е. при

i(t) = -(qj - наведенный заряд троса), эта dt

теорема в форме [2] имеет вид di (t) = Т л dzVo(t )(eot + E'ot ) +

+{

dE-

где 8q—- - плотность тока смещения, созданного dt

движущимся элементарным зарядом dq = -тr/dz ;

EOT, EOt - напряженности фиктивного электростатического поля, создаваемого в точке расположения заряда dq системой проводов и тросов при условии, что наблюдаемый трос (рис. 1) имеет безразмерный потенциал фт = 1 , а провода и остальные тросы заземлены (ф1 = ф2 = ... = фп = 0), т.е. это

поле условных линейных зарядов ТТ, т1....Тп и их

изображений ТТ =-тТ , Ti=-t1 , , размерность Кл

которых dim Т =

м  В

р

Рис. 1. Основные параметры системы провода-тросы и канала грозового разряда

а величина определяется при решении системы уравнений

- ат„т п

1 — а-

ттТ

+ аТііі +... + (

к7п 1п ’

0 = аи Ту + аі і'гі +... + аіпТп,

0 — а пт Ту + а п-1 Т +... + с

кпіЧ а п

атт, аіт, ... апп - потенциальные коэффициенты.

Дальнейшее рассмотрение методики определения і(Ї) выполним для уединенного троса (рис. 1) в двух вариантах: без учета и с учетом запаздывания сброса зарядов тросом относительно нейтрализации зарядов канала. И затем обобщим методику на систему провода-тросы. Без учета запаздывания имеем

2пє0

1п

Еот =

2т

Ттг0

2пє0Г г |п2Ьт

г = ^ $ + ( - Ьт )2,

Ео т =

г' 1п

2Ьт

Л )2;

Уф

1п

2НТ/

* - Ьт

$ +

Л

где т - радиус троса.

Задача о поле тока смещения как в лидер-ной стадии, так и в стадии главного разряда была решена в работе [4], поэтому подынтегральное выражение во втором слагаемом (2) считаем известным и вычисляем ток /@), обусловленный движением заряда тсбг . Фиксируя момент наблюдения t, определяем положение фронта нейтрализации

гф

0

и, полагая выполненным условие divVФ (?) = 0 , интегрируем (2) по г от 0 до гф^). В итоге получаем

і

1п

*УУ|п>/Ь +

%

| (2от + Е0т) єо "дт .

і

т

т

р дБ 5

В этом выражении є0 — — 8 - плотность тока ді

смещения, обусловленного перемещением заряда в канале. Его составляющие по осям т и р (рис. 1) цилиндрических координат имеют вид [4]

* - гФ

S7 =--

4п

ф

(z - zф )2

( + ^ )2

sp =-

4п

1 1

p2 +(z - Zф )2 32 P2 +( + Zф )2

где с =

Уг~

/ VS0^0

В [4] приводится уравнение для линий 8; его решение и анализ полученных результатов выполнены в [5].

Таким образом, на основе приведенных формул для 82, 8р может быть вычислено второе слагаемое в (7), т.е. рассчитана форма импульса наведенного тока /А).

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

поля вблизи троса, где напряженность Еот значительна. Однако в этой области плотность тока смещения допустимо считать постоянной по условию гт << Ь. Получающийся приближенный (оценочный) результат для () будет соответствовать формулам теоремы Шокли-Рамо [1, 3], учитывающим только ток переноса.

Итак, определяем параметры импульса наведенного тока /А), используя только первое слагаемое в (7). Предварительно заметим следующее: заряд лидерного канала может быть как положительным, так и отрицательным, тогда как скорость

продвижения фронта Vф ()= всегда положительна: Vф ()> 0 . Поэтому при тп < 0 имеем, согласно формулам (1), (7), для принятых положительных направлений этих токов /м(0 < 0, /^) > 0, т.е. в этом случае в заземлитель поступает положительный токовый импульс, что соответствует оттоку с троса в землю положительного наведенного заряда дт > 0. При тл > 0 картина явления - обратная.

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

di п ь

-----= 0 , решение которого дает гф = Лт и соответ-

dZф

ственно

M

MAX

ад ьт+(2hT)2

In

2T

b

{8)

где hT - решение уравнения

ьт

Ьт = J Vф (u)dи.

{9)

0

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

В качестве иллюстрации применения предложенной методики и оценки параметров ожидаемых импульсов наведенного тока выполним расчет i(t) при следующих исходных данных: радиус троса - 5-10'3 м; варианты высоты подвеса - hn = 10 м, Лт2 = 1 м, hT3 = 0,2 м; варианты расстояния до точки удара молнии - b1 = 20 м, b2 = 100 м; ток молнии задан в форме биэкспоненциального импульса, соответствующего 50 %-му току молнии, для которого, согласно [5], имеем iM max ~ 37 кА, Vo max = 0,31 -3108 м/с, тл = 410-4 Кл/м. Форму этого импульса описывает выражение

f _ t/ - t/ \

M

(t )| = im

т2

- e

/т1

{10)

где т2 = 46,7 мкс; п = — = 0,0214 1/мкс; т = 1,77 мкс; т2

1

5 =— = 0,565 1/мкс; im = 43,77 кА. т1

Максимальное значение тока имеем при tm = 6,02 мкс. Скорость продвижения фронта нейтрализации, м/с, определяется как

- e

VФ (t) = m(

т л тл Соответственно, по {6) получаем

).

{11)

zф (t) = 1,094 • 102

гф^т) ~ 430 м.

Длительность процесса нейтрализации, т.е. длительность импульса наведенного тока, определяется заданной высотой облака Л0 = 2000 м (рис. 1). Найденное значение этого параметра - И = 26,5 мкс.

Расчеты кривой наведенного тока /ф были выполнены по первому слагаемому выражения (7) и формулам (8) и (9). Полученные результаты иллюстрирует рис. 2 для значений Ь1 и Ь2, указанных выше.

т

З

а)

б)

Рис. 2. Кривая наведенного тока Щ рассчитанная по первому слагаемому (7): а - при Ь1 = 20 м и различных вариантах высоты подвеса троса ЛТ1..3, приведенных выше; б - при Ь2 = 20 м и различных вариантах высоты подвеса троса ЛТ1..3, приведенных выше

Значение времени, за которое наведенный ток достигает максимума, вычисляется по формуле (9). Для заданных высот подвеса троса Лт имеем ^ = 0,62 мкс, ^2 = 0,22 мкс, ^3 = 0,11 мкс.

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

£з =

и, соответственно, с учетом формул (2), (5) di (£) = ТЛ^Ф (£ - £З ) (от + Е0т),

У Ь2Т +( - Иг)2

VфU -

Ш2Ит/

I

Z - Ьт

dz.

Ьт2 +(z - Ьт) Ьт2 +( + Ь

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

заряды, для которых имеем t > tз; верхний предел в интеграле (14) вычисляется по уравнению

-|>/ьт+((-ьт?"

I Vф (о)/и, £ > ьТ + Ь2,

с

которое для тока молнии, заданного в форме (10), преобразуется к виду

zt =^Д

тл 1п

1- ехр

1- ехр

£-

точке

с

а кривая /(^ начинается

£ = tЗ (0) = —Ьт + И2 . Поскольку из процесса формирования /^) выпадает участок канала < г < гФ^), получаемая кривая наведенного тока в системе отсчета £' = £ - З (0) = £ - Ьт + И2 мажорируется най-

денной выше зависимостью (), даваемой формулой (7). Результаты расчетов по (14), (15) приведе-

ны на

Сравнение кривых /(§ и /(Р) на рис. 2, 3 позволяет сделать вывод о том, что для оценочных расчетов импульсов наведенного тока допустимо использовать сравнительно простую формулу (7) (первое слагаемое), рассматривая данную зависимость с учетом сдвига во времени для переменной

£' = £ -

При наличии нескольких заземленных тросов (тросовые экраны [6], защитные тросы воздушных линий электропередачи и другие устройства с заземленными объектами иных типов) расчет наведенных токов проводится для каждого из них в отдельности; ток в заземлителе находится путем суммирования с учетом запаздывания tз1, tз2, и т.д. Соответствующие расчетные соотношения запишем, полагая, что в системе (рис. 1) имеется т заземленных тросов и п - т проводов. Составляющая фиктивного электростатического поля, параллельная

вектору скорости Vф (£), вычисляется при определе-

нии тока с к-го заземленного троса (к = 1, 2, формуле

т) по

п М Е(к) = у %

Е01 ~ А п

I 2пе0

i=1'

z -

Ь2

в которой плотности условных зарядов т0/ находятся при решении системы уравнений (3) для случая 9^ = 1, ф,- ^ = 0 .

0

х

Л

а)

б)

Рис. 3. Кривая наведенного тока / (£'), рассчитанная с

учетом запаздывания: а - Ь1 = 20 м и различных вариантах высоты подвеса троса ЛТ1..3, приведенных выше; б - при Ь2 = 20 м и различных вариантах высоты подвеса троса ЛТ1..3, приведенных выше

Выполняя преобразования без учета запаздывания, получаем

^ . п 7x2

0 /=1

=ч« (2ф),

где ф0к^(ф) - потенциал фиктивного электростатического поля, созданного зарядами * в точке

I = Хф, р = 0 (фронт волны нейтрализации).

Для тока /(к>А), стекающего с к-го троса в за-землитель, согласно (7), имеем формулу

)(£ ).-ТМ) V ,<« .

2"° '=1 +

= ТлVф (£)ф0к}[Zф (£)].

Поскольку ф0к^(ф) как электростатический

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

ниях 2ф, рф, достигаемых фронтом к фиксированному моменту времени t.

Расчет по выражению (18) выполняется для всех заземленных тросов, найденные токи ^(О, к =1, 2, ...т, позволяют найти суммарный ток, протекающий в заземлителе:

,

к=1 V с У

В качестве иллюстрации применения изложенной методики рассчитаем наведенный ток, стекающий в заземлитель с защитных тросов ЛЭП-500 кВ при ударе молнии в землю в середине пролета линии. Расчет поясняет рис. 4; геометрические параметры системы провода-тросы соответствуют данным, взятым из справочника [7].

Для этой линии имеем:

• провод ЭхЛС 400/51: гп = 1Э,8-10'3 м, радиус расщепления б = 0,4 м;

• тросы С-70: гТ = 5,5-10'3 м, dТ = 16,2 м;

• в середине пролета Лп = 16,6 м, ЛТ = 24,2 м;

• расстояние от троса до места удара молнии: Ь1 = 20 м;

• параметры импульса соответствуют формуле (10).

При решении задачи последовательно выполняем следующие расчеты:

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

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

Результаты расчетов приведены в таблице,

для безразмерных величин т (к) = Т' зуемых в формулах (16)-(18).

'2пеп

исполь-

Относительные значения условных зарядов

/ 1 2 3 4 5 6

т? -0,119 0,00869 0,00941 0,00933 0,00941 0,00653

т*(2) 0,00869 -0,11949 0,0024 0,00238 0,0024 0,00653

/ 6 7 8 9 10 11

т*((1) 0,00653 0,00647 0,00653 0,0024 0,00238 0,0024

Т*/2) 0,00653 0,00647 0,00653 0,00941 0,00933 0,00941

0

Рис. 4. Геометрические параметры системы провода-трос 2. Поскольку радиус расщепления при расчетах /(к) (£) по вы-

Ь1 - I dП - dг

ражению (18) допустимо использовать заряды фаз

Т*(к) = Т*(к) + *(к) + *(к) ;

% А 3 +% 4 + т 5 ;

Т*(к) = Т*(к) + Т*(к) .

1 С ~ 1 9 + 1 10

Т*(к) = Т*(К) + Т*(К) + Т*(к) ;

% В ~т 6 +т7 +т8 ;

т Ц. Согласно данным таблицы,

имеем Т

Т*(2) =

А) = 0,02815, Т*В> = 0,01952, т*(1) = 0,00719; А = 0,00719, т*В2) = 0,01952, т*С° = 0,02815.

3. Выполняем расчет токов /

стекающих с тросов в заземлитель. Согласно выра жениям (17), (18), имеем

/(к)(£) = -тV (£)Хт*<к)1п

ф )

/ = 1, 2, Л, В, С.

4. По формуле (19) определяем ток в зазем-

лителе:

/(£) = /(1) I£ - Ь.2 + Ь.2 "1 + /(2) {£ - Ь2 + Ь

£ > .д/Ь.2 + Ь. . с

Результаты расчетов, выполненные по методике, иллюстрирует рис. 5.

Рис. 5. Результаты расчета наведенного тока, стекающего в заземлитель с защитных тросов при ударе молнии в землю, в середине пролета линии ЛЭП-500 кВ: /1(Ц и /2(0 -токи, стекающие с тросов 1 и 2 (см. рис. 4) в заземлитель при ударе молнии; ^ЭЕетОД - суммарный ток, протекающий в заземлителе

Заключение

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

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

Применение рассмотренного варианта теоремы Шокли-Рамо в тех случаях, когда для многопроводной системы выполняется расчет только интегральных характеристик, является предпочтительным, по сравнению с расчетом распределений зарядов через поле методом интегральных уравнений [8].

Список литературы

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

2. Фельд Я.Н. Наведение токов движущимися зарядами // Докл. АН СССР. - 1953. - Т. 93. - № 3 - С. 447-450.

3. Герштейн Г.М. Моделирование полей методом электростатической индукции. - М.: Наука, 1970.

4. Слышалов А.В. Распределение тока смещения при грозовом разряде: Межвуз. сб. науч. тр. «Высоковольтные техника и электротехнология». Вып. 2. - Иваново, 1999.

5. Слышалов А.В. Разработка методов расчета механических сил и моментов, действующих на фазные провода и защитные тросы ЛЭП в поле лидера молнии: Дис... канд. техн. наук: 05.14.02; 05.14.12. - Иваново, 2001.

6. Барсуков А.И., Васильев А.В., Мисриха-нов М.Ш., Седунов В.Н., Токарский А.Ю. Применение тросовых экранов для снижения напряженности электрического поля ВЛ СВН.

7. Справочник по проектированию линий электропередачи / М.Б. Вязьменский, В.Х. Ишкин, К.П. Крюков и др.; Под ред. М.А. Реута и С.С. Рокотяна. 2-е изд. пере-раб. и доп. - М.: Энергия, 1980.

8. Слышалов В.К., Киселева Ю.А. Расчет распределений зарядов, наводимых в системе параллельных проводов полем лидера молнии: Межвуз. сб. науч. тр. «Высоковольтные техника и электротехнология». Вып. 3. -Иваново, 2003. - С. 91-96.

Слышалов Владимир Константинович,

ГОУВПО «Ивановский государственный энергетический университет имени В.И. Ленина», доктор технических наук, профессор кафедры электрических систем, телефон (4932) 26-99-21, zav_es@es.ispu.ru

Сорокин Александр Федорович,

ГОУВПО «Ивановский государственный энергетический университет имени В.И. Ленина», кандидат технических наук, доцент, декан электроэнергетического факультета, телефон (4932) 26-99-23, deaneef@eef.ispu.ru

Слышалов Андрей Владимирович,

ОАО «Ивэнерго»,

Начальник отдела, кандидат технических наук, zav_es@es.ispu.ru

Киселева Юлия Александровна,

ООО «РЗА-сервис».

ведущий инженер, кандидат технических наук, zav_es@es.ispu.ru

Кандалов Юрий Владимирович,

ГОУВПО «Ивановский государственный энергетический университет имени В.И. Ленина», аспирант,

телефон (4932) 26-99-21, zav_es@es.ispu.ru