Влияние ближнего разряда молнии на линию передачи

УДК 51.73

Ю. М. Григорьев, М. Н. Орлова ВЛИЯНИЕ БЛИЖНЕГО РАЗРЯДА МОЛНИИ НА ЛИНИЮ ПЕРЕДАЧИ1

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

Ключевые слова: молния, линия передачи, многолетняя мерзлота, индуцированное напряжение.

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

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

В предыдущих работах авторов аналитически и численно решены задачи о ВТН в предположении о мгновенном характере разряде молнии, рассмотрены случаи разряда молнии в землю и между двумя облаками [1; 2]. В нашей работе [3] представлена математическая модель ВТН, учитывающая зависимость тока молнии от времени и аналитически решена задача о ВТН. В данной статье приведен другой вид решения для случая произвольных погонных параметров Я, С, Ь, О, проведены новые численные расчеты.

Математическая модель. Рассмотрим линию передач с одним проводником. Ось х направим вдоль этого

бесконечного проводника. Проводник характеризуется распределенными параметрами Я, С, Ь, О - сопротивлением, емкостью, индуктивностью и коэффициентом утечки, рассчитанными на единицу длины. Грозовое облако моделируем точечным зарядом Q, расположенным на высоте И над кабелем и на высоте И +1 над проводящим полупространством (проводящий слой земли), параметр I описывает толщину многолетнемерзлого слоя грунта. Следовательно, проводник находится в поле заряда и его электростатического изображения (рис. 1).

Рис. 1. Геометрия модели

Проводник считаем заземленной в бесконечно удаленной точке, потенциал Земли у = 0 , т. е. индуцированные на проводнике заряды создают такой потенциал, который компенсирует потенциал, создаваемый двумя точечными зарядами Q и -Q .

Пусть в момент времени t = 0 облако начинает разряжаться, зависимость заряда облака от времени выражается функцией Q(t), t > 0, тогда ток молнии равняется dQ(t)/ dt. При уменьшении заряда облака на dQ за время dt те заряды на линии, которые компенсировали потенциал заряда dQ (и его изображение - dQ), стано-

1 Работа поддержана грантом № 06-08-96020 конкурса РФФИ-Дальний Восток.

вятся свободными и начнут «разбегаться», возбуждая ВТН в линии. Другими словами, разряд грозового облака эквивалентен действию некоторого источника зарядов в линии. За время М на элементе ёх этот источник вырабатывает количество зарядов dq1, равное

dq1 = Сёхёи = С/

где

/ (х) =

1

ґ

4пє„

1

1

Л

\1 к2 + х2

-у/(2/ + И)2

X + Си, + Ои = С/ (х) ^, ё,

и(х,0) = 0,і(х, 0) = 0.

и( х, ,) =

1

2аС

I е

1

2аС

х

С - х

С - х ) -X,,----------------------- I е

С - х

2аС

( ( X/о

К (, - т) -

(х - О2

а1

к

(, - Т)2 - ^

к* (, - т) -

(х - С)2

ёт +

с - х

2аС

к* (, - т) -

к

(, - т)2 - (х - ^

^1

V

к* (, - т) -

(х - с)2

(х- о2

ёт +

Кае

а

X/ ,

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

Аналогично [4] выведем систему телеграфных уравнений с учетом действия такого источника зарядов. Применяя закон Ома к участку длиной ёх, получаем

-ихёх = Кёх + Ы,ёх.

Количество зарядов, притекающих в элемент проводника ёх за время ё, равно

dq2 = (і (х,,) - і( х + ёх,, ))ё, = -і^хйї, за это время ё, на участке ёх за счет действия указанного выше «генератора» дополнительно образуется количество зарядов dq1. Количество заряда dq3, которое идет на зарядку элемента ёх за счет электроемкости равно dq3 = С(и(х,, + ё,) - и(х,,))ёх = Си^хё,, количество заряда dq4, теряющегося за счет несовершенства изоляции:

dq4 = Оиёхё,.

Согласно закону сохранения заряда, dq1 + dq2 = dq3 + dq4.

Подставляя сюда выражения для зарядов, получаем

С/ (x)dxdQ - ixdxdt = Си^хё, + Оиёхё,.

Из этого соотношения и уравнения (3) получаем искомую неоднородную систему телеграфных уравнений и, добавляя нулевые начальные условия, получаем задачу Коши, описывающую ВТН в линии передачи, с учетом зависимости тока молнии от времени:

их + Ы, + Кі = 0, -¥ > х >¥, , > 0,

а

- а

е~Хт *(С,т)/ 0

К (, - т) -

(х - О2

а1

____-X,

&>,

(

і(х,,) = --

<| е Х (*(х + а(, - т), т) - *(х - а(, - т), т))ёт +

ке

Ь а

2а

'іёх I

х-С

X /1

- а(, - х) (

(

(,-Х)2 - <хх-£

У(С,х);

к* (, -х) -

(х -С)2

X = -

к =

2ІС

1ЯС - “I а2 = -±_

2ЬС ' ЬС Рассмотрим неискажающую линию передачи, т. е. когда погонные параметры удовлетворяют соотношению ЯС = ЬО . В этом случае

ч = Я

4 = - ~Ь ’ ^ = 0, и формулы для величин ВТН существенно упрощаются:

и( х,,) =

1

2аС

, , +

С - х

I е

Аналитическое решение. Аналитическое решение задачи (4), описывающей ВТН, получаем в следующем виде:

і( х,,) = --

а 0

С - х 1 -Xі , ,---------- I е

а

-Хх /

1

2аС

Г е “ (*(х + а(, - т), т) -1 - *( х - а(, - т), т))ёт

Прямоугольный импульс тока молнии. Рассмотрим случай, когда ток молнии имеет вид прямоугольного импульса с длительностью t0, с амплитудой 10 = Q0 / ^ , где Q0 - первоначальный заряд облака, в этом случае заряд облака выражается функцией

Q

где

(1, при , > 0,

^ (0, при , < 0

- единичная функция Хевисайда. Тогда

> =-.,(,0 -.),

откуда

0

х

У

х

0

0

- а,

а

У

¥

) С00

*0

/(х)п(*0 -* )•

и(х *) = -0С I / П [ *0 - * -^Ч ек а «С -

2аС '} , 1 п 1

2аС

I /(Оп

С - х 1 -^

*0 - * + ---------------| е

«с -

С - х

.ее

2а*0

| «С I /(Ое-

[

1т [ XI0 1? -

( (

к

(* - т)2 -

2 (х - С)2

-л

к* (* - т) -

2 (X - С)2

(X - О2

п(*0 - т)Жт -

2а/„

1

(-1)І0 ч

( (

к

(* - т)2 -

2 (X - С)2

Ч(* - т)2 -

(X - С)2

П(*0 - т)«т -

ЯС00 ае

-и * х + а

2*0

-{ Ж т I е^ / (С)/ 0

(

0 0 х - а

К\(* - т)2 - ^ '

V

а

п(*0 - т)«С

У

/ (х, *) =

00 Сае~

2*0

*

Ге~ (у(х + а(* - т)) -Л -у(х - а(* - т))) п(*0 - т)«т -

7^-4 ^ - X* * х + а

к00Се Г Жт Г е-X,

2а*0 і

х/і

х - С

е * - т) і (* - т)2 - (^ а

1(*- т)2 - (х - О2 2

¥(С,т)>

В случае неискажающей линии передач решение (7), (8) будет иметь вид

и(х, *) = - I

I /(Оп

( /<0ц (*0-'- £-£) а

2а*г ^ -іі

0 х - а* - X---------------

0 е а «С -

0

2а*0

надущельем), Я = 10 3 Ом/м, С = 10~10 Ф/м, І = 10-7 Гн/м, 1

Таким образом, решение (5) задачи (4) будет иметь следующий вид:

£ - х ^ -’— а

ех + а1 / 9» \ Г — х

0 - г г — ^ -^

О = 10— .

Ом • м

Результаты расчетов в разные моменты времени для безразмерной силы тока и напряжения при разряде молнии в землю приведены на рис. 2, 3.

1 [к1

для прямоугольного импульса тока молнии

Іі [х ;

Аналогично преобразуется формула (6) для силы тока:

1

( , ) 2*0 Г -¥(х - а(Г - т))) п(*о - т)А.

Численные расчеты. Численные расчеты проведены для случая прямоугольного импульса тока молнии с длительностью /0 = 50 мкс и амплитудой J0 = 20 кА , Q0 = 1 Кл при И = 1 500 м, I = 250 м (сухой грунт над вечной мерзлотой или скальный грунт, воздушная линия

Рис. 3. Волна напряжения для прямоугольного импульса тока молнии

Расчеты показывают, что максимальное значение силы тока в линии при приведенных выше параметрах достигается в точках, расположенных на значительном расстоянии. Сравним полученные результаты с результатами расчетов по модели с мгновенным разрядом молнии [5]. Распределения силы тока в линии при достижении ее максимального значения, рассчитанные по этим двум моделям при одинаковых параметрах, приведены на рис. 4

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

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

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

а

х

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

i [==]

10 го зо

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

0.31 -Нгнов tmntp разряд

0.25І

0. l.j

0. A.

0.05 - с уче том Ер ексен

зо -20 -ю ю го зо

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

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

1. Григорьев, Ю. М. Исследование влияние электромагнитного воздействия на кабельные линии / Ю. М. Григорьев, В. В. Наумов, П. И. Николаев // Физика высокоширотной ионосферы и распространение электромагнитных волн. Якутск, 1988.

2. Григорьев, Ю. М. Математическая модель грозового перенапряжения в линии передачи при разряде молнии между двумя облаками / Ю. М. Григорьев, М. Н. Орлова // Динамика спл. среды. Новосибирск, 2004. Вып. 22. С. 53-56.

3. Григорьев, Ю. М. Расчет волны тока индуцированной разрядом молнии в кабельных линиях / Ю. М. Григорьев, С. Н. Еремеев, В. В. Наумов, А. А. Семенов // X Междунар. Вроцлав. симпозиум по электромагнит. совместимости (26-29 июня 1990, г. Вроцлав). Вроцлав, 1990, С. 247-252.

4. Тихонов, А. Н. Уравнения математической физики / А. Н. Тихонов, А. А. Самарский. М. : Наука, 1977.

5. Григорьев, Ю. М. Математические модели грозо-

вых перенапряжений в линиях передач / Ю. М. Григорьев, М. Н. Орлова // VII Лаврентьевские чтения : науч. конф. (7-11 апр. 2003, г. Якутск) : сб. ст. Якутск, 2003. Т. 1.

С. 29-33.

6. Rachidi, F. Calculation methods of the horizontal component of lightning return stroke electric fields / F. Rachidi, М. Ianoz, С. А. Nucci, С. Mazzetti // XI Intern. Wroclaw Symposium and Exhibition on Electromagnetic compatibility (September 2-4 1992, Wroclaw). Wroclaw, 1992. P. 452-456.

7. Yoshihiro, B. Voltages induced on an overhead wire by lightning strikes to a nearby tall grounded object / B. Yoshirio, V. A. Rakov // IEEE transaction on electromagnetic compatibility. 2006. Febr. Vol. 48. N° 1. P. 212-224.

Yu. M. Grigoriev, M. N. Orlova INFLUENCE OF A LIGHTNING ON A TRANSMISSION LINE

Analytical and numerical solutions of a current and voltage wave problem in a transmission line while a lightning discharge are given. The presented model takes into account the time dependence of the lightning current.

Keywords: lightning, transmission line, permafrost, induced voltage.

2б