Научная статья на тему 'Влияние стрелы провеса провода воздушной линии электропередач на параметры ЭМП'

Влияние стрелы провеса провода воздушной линии электропередач на параметры ЭМП Текст научной статьи по специальности «Физика»

CC BY
2011
161
i Надоели баннеры? Вы всегда можете отключить рекламу.

Аннотация научной статьи по физике, автор научной работы — Белицын И. В., Макаров А. В.

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

i Надоели баннеры? Вы всегда можете отключить рекламу.
iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.
i Надоели баннеры? Вы всегда можете отключить рекламу.

THE INFLUENCE OF THE OVERHEAD LINE SLACK ON ELECTROMAGENTIC FIELD PARAMETERS

The influence of the overhead line slack on the value of the electric field intensity is shown, the technique of intensity calculation is cited, recommendations on definition quantity of the flights, necessary for maintenance of sufficient accuracy of calculations, are given.

Текст научной работы на тему «Влияние стрелы провеса провода воздушной линии электропередач на параметры ЭМП»

УДК 621.311:537.871.3

ВЛИЯНИЕ СТРЕЛЫ ПРОВЕСА ПРОВОДА ВОЗДУШНОЙ ЛИНИИ ЭЛЕКТРОПЕРЕДАЧ

НА ПАРАМЕТРЫ ЭМП

И.В. Белицын, А.В. Макаров

Алтайский государственный технический университет им. И.И. Ползунова, г. Барнаул E-mail: [email protected]

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

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

Одним из факторов, влияющих на точность расчетов, является наиболее точное математическое описание «геометрии» ВЛЭП. Поэтому необходимо учитывать, что провод, закрепленный в двух точках на опоре, не является горизонтальным проводником. Существующие способы учета провисания не имеют под собой научного обоснования. Наиболее распространенным является прием, применяемый как за рубежом [1], так и у нас [2], когда реальные провода и тросы заменяют бесконечно длинными горизонтальными прямыми, расположенные на высоте Иср=И-2/3/, где Н-высота подвеса; / - стрела провеса. Непонятно, почему для расчета параметра ЭМП провод располагается на высоте центра тяжести провода.

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

y = y ch

J ^max

Уm

(1)

Рис. 1. Провисание идеально гибкой нити

Уравнение (1) можно разложить в функциональный ряд Маклорена

y = y ch

max

ym

= ym

1 +1

211 у 1

+ 4

411 У

(2)

. (2«)! I Ут

Из свойств гибкой нити следует:

• для низшей точки провисания провода или троса;

= УтахГ, (3)

где ст0 - механическое напряжение в низшей точке провода или троса Н/мм2; у - удельная нагрузка провода или троса Н/(м.мм2).

• в точках подвеса провода или троса А и В

°Л = (Утах + /Л )7=°0 + /лУ;

= (Утах + /Б )У = °0 + /зГ, (4)

где /А, /в - стрела провеса относительно точек подвеса А и В соответственно. Если опоры установлены на идеально ровной местности, то /А=/в=/ и называется стрелой провеса, тогда (4) можно переписать как

= = (Утах + Г)У =°0 + /7 ■

Как видно из рис. 1, стрела провеса

/ = Н - Утах, (5)

где Н - высота опоры, на которой крепиться провод или трос, или с учетом (1)

f = УmaxCh

Уm

У m

(6)

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

Поскольку при расстановке опор по горизонтальной плоскости, ур. (1-6), координаты точек

подвеса хАВ=1/2, где I - расстояние между опорами, воспользовавшись (2) и (3), получим

У

сЬ

1 + -

I

2а0

I

2ст0

+4

-и=

I

2оо

(2п)!

I

2ст„

-1

I =

х л

2! 22с

-у-

± Л

4! 24о

1

(2п)! 22

(7)

Итак, математической моделью для определения геометрического распределения заряда может служить выражение (1) или (2); для обеспечения достаточной точности допустимо ограничиться тремя членами ряда, а для пролетов менее 700 м достаточно первых двух членов

У ~ У,

+1

тах 2 у

у та

(8)

1 - 4 *

(9)

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

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

Построение кривой 1 шаблона выполняется по уравнению параболы:

У = кх\ (10)

где х - переменная величина е[-1гб/2; 4^/2]; к -постоянная шаблона, зависящая от физико-механических характеристик провода и приведенной

нагрузки на провод. Выражение (10) представляет собой выражение, аналогичное (8).

Координата точки максимального провеса ушш зависит от стрелы провеса / Стрела провеса, в свою очередь, определяется физико-механическими свойствами самого провода, действующими на него механическими нагрузками и расстоянием между опорами. Для ее практического определения достаточно двух первых членов выражения (7), а для пролетов менее 700 м достаточно одного члена, ур. (9):

Рис. 2. Шаблон (общий вид) для размещения опор по профилю. Кривая: 1) провисания провода; 2) смещена по вертикали на величину габарита, «габаритная» кривая; 3) смещена относительно первой на активную высоту опоры, «земляная» кривая

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

Н - И, = к\

к = 4

Н - И

1 габ

(11)

Постоянную шаблона в выражении (10) можно определить, анализируя выражение (8), приняв в нем УмхтК, тогда

У = И +

1 х

2 И, "

Сравнивая (10) и (12) находим

к =

1

2И , •

(12)

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

(13)

Для дальнейших расчетов можно использовать постоянную шаблона, рассчитанную как по (11), так и по (13).

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

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

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

а

У

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

Напряженность от заряда с элементарного участка провода

на 0Х

СЕ+=-^- г, 4пе0г

(14)

где е0=8,85.10-12 Ф/м - электрическая постоянная; йд=тсИ - элементарный заряд; г - расстояние от точечного заряда до произвольной точки Ы(Х,У^

(рис. 3) г =Л/ (X - х)2 + (У - у )2 + (г - 7 )2.

СЕ = -

4пе0г'

-г .

(15)

С1 =

008 в '

(16)

008 в =

лА+^в'

(17)

Угловой коэффициент касательной tgв = СУ ,

т. е. с учетом (10), получим:

tgв = 2кх. (18)

Будем считать, что заряд равномерно распределен вдоль проводов ВЛЭП с линейной плотностью т, тогда, учитывая выражения (14-18), получим

Сц =т/1 + (2кх)2 Сх. (19)

Проекции напряженности от всего провода найдем интегрированием по всей длине провода:

Е =_ь_

х 4пее0

ф + (2кх)2 (X - х)

((X -х)2 + (У -кх2 -кг)2 + (г -7)1)2 ф + (2 кх)2 (X - х)

Сх;

(20)

((X - х)2 + (У + кх2 + кг)2 + (I - 2)2)

2\1

на 0У

Е =_ь_

У 4пее0

ф + (2кх)2 (У - кх2 - кг) ((X - х)2 + (У - кх2 - кг )2 + (I - 2)1)2 ф + (2кх)2 (У + кх2 + кг)

Сх;

((X -х)2 + (У + кх2 + кг)2 + (I - г)2)2

(21)

на 0Z

Рис 3. К определению напряженности от одного пролета в пространстве

Напряженность от зеркального изображения элементарного заряда

Сд

Е = —ъ_

1 4пее0

Ф + (2кх)2 (г - 7)

((X - х)2+(У - кх2 - к )2 + (г - 2)1)2 + (2кх)2 (г - 7) ((X - х)2+(У+кх2 + кг )2+(г - 7 )2)2

Сх.

(22)

где г' - расстояние от изображения точечного заряда до произвольной точки М,

г' = ^(X - х)2 + (У + у)2 + (г - 7)2.

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

Сх

Данное выражение действительно для любых точек провода, поскольку у параболы нет точек, касательная к которым составляла бы угол 90°, т. к. ве(-90°; 90°) то

1_

Модуль напряженности в произвольной точке

Е (М) = 7 Е2 + Е2 + Е2. (23)

По полученным зависимостям (19-23) построим характерные зависимости напряженности электрического поля от координат точки наблюдения, при этом в качестве исходных данных примем длину пролета 4^=300 м, линейную плотность заряда на проводе т=3,72 мкКл/м, провод подвешен на высоте Н=25 м, нормированный габарит кг=8 м.

На рис. 4 показаны зависимости напряженности электростатического поля ДДсоп^сош^ от расстояния вдоль одного пролета Х, на постоянной высоте У=1,8 м, в зависимости от расстояния до плоскости провеса провода Д при условии, что стрела провеса вертикальна и ее координата ¿=0.

На рис. 5 показаны зависимости напряженности электрического поля от расстояния до плоскости провода Z на постоянной высоте 1,8 м, в зависимости от расстояния вдоль одного пролета X, при условии, что стрела провеса вертикальна и ее координата ¿=0.

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

• напряженность поля максимальна непосредственно в плоскости нахождения провода;

• на расстоянии, приблизительно равным трети длины габаритного пролета от провода, максимальная напряженность уменьшается на порядок;

х

X

X

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

E _ т С_H + kx2_

n 2пее0 С ((H + kx2)2 + (X - x)2)3

2

Для построения зависимостей напряженности электрического поля на поверхности земли будем использовать те же исходные данные, соответствующие ВЛЭП напряжением 500 кВ, а именно: линейная плотность заряда на проводе т=3,72 мкКл/м, габаритный пролет 300 м, &=7,56.10-4 м-1. Поскольку каждый пролет создает электрическое поле, то результирующая напряженность определяется согласно принципу суперпозиции

E (X) =

2тг°

N L

<Х i

i=-N -L 2

H + kx2

((H + kx2)2 + (X +i- 1габ -x)2)3

Рис. 4. Напряженность электрического поля вдоль одного пролета в зависимости от расстояния до плоскости провода на высоте 1,8 м. 1) 1=0 м; 2) 1=10 м; 3) 1=20 м; 4) 1=30 м; 5) 1=40 м; 6) 1=50 м

W1 + 4k2 x2 dx.

(25)

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

I Е, кВ/н

Рис. 5. Напряженность электрического поля от расстояния до плоскости провода, в зависимости от расстояния вдоль одного пролета на высоте 1,8 м. 1) Х=0 м; 2) Х=30 м; 3) Х=60 м; 4) Х=90 м; 5) Х=120 м; 6) Х=150 м

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

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

Напряженность, создаваемая всем проводом (24), будет иметь только нормальную составляющую

W1 + 4k2 x2 dx.

(24)

-uo -iao -so -eo -so о so eo so iao uo

Рис. 6. Продольная эпюра распределения напряженности от ВЛ 500 кВ при учете пролетов: 1) пяти; 2) трех; 3) одного

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

E — E

-1°0%' (26)

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

эт

где Е - напряженность, вычисленная по (25) для N<10 пролетов.

На рис. 7 показана относительная погрешность определения напряженности поля (26) в зависимости от расстояния от оси симметрии при варьировании количества смежных пролетов.

Анализируя рис. 7, можно сделать выводы:

• учитывая влияние одного пролета, минимальная погрешность 0,28 % при расчете поля наблюдается возле оси симметрии; максимальная погрешность достигает 50 %, наблюдается возле координаты точки подвеса провода;

• для инженерных расчетов при точности «5 % можно учитывать один пролет на 75 % его длины;

Рис. 7. Относительная погрешность при учете пролетов: а) одного; б) трех (крив. 1), пяти (крив. 2)

• при учете влияния двух смежных пролетов погрешность становится пренебрежимо малой «0,12 %, поэтому во всех расчетах, связанных с определением напряженности электрического поля и его моделированием, без потери точности можно учитывать только три пролета, которые в дальнейшем будем называть расчетными.

Зависимость напряженности электрического поля от координаты (25) - это некоторая периодическая четная функция с периодом, равным длине габаритного пролета. Любую периодическую функцию f(x)=E(X в интервале (-/, /), где /=1аб/2, можно разложить в ряд Фурье [7]:

f (x) = ^ + Х

an cos

nnx l '

(27)

+0,194cos3nX - 0,056cos4nX+...

150

150

(28)

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

E = 2-

2пе0 HCp

(29)

Среднюю высоту подвеса провода определим, приравняв (29) к постоянной составляющей а0/2:

Я„ = 2-^-.

применяя численные методы вычисления интегралов, определяем для данного случая а0=22,412, а1=5,431, а2=0,139, а3=0,194, а4=-0,056.

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

Е (X) = 11,206 + 5,431008 +ОДЗ^оэ2^ +

Поскольку напряженность изменяется в зависимости от координаты Х и известно среднее значение на всем протяжении под пролетом, равное постоянному члену ЕС1=а0/2 в выражении (27), то можно определить среднюю или эффективную высоту под-

В приведенном численном примере Нср«12 м, это на 12,5 % меньше средней высоты подвеса проводов, рассчитанной по формуле Н,=Н-2/3у=13,7 м. Очевидно, что средняя высота подвеса провода будет различна при определении Е на различных высотах от поверхности земли, т. к. с увеличением высоты (координата У), напряженность поля увеличивается, при прочих равных условиях.

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

n=1

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

1. Kumari M.N., Kumar O.R., Nambu-diri P.V.V., Srinivasan K.N. Computation of Electrical Environmental Effects ofTransmission Lines // In: Proc. 11th Int. Symp. on High Voltage Eng. - London, 1999.

2. Справочник по проектированию линий электропередачи / Под ред. М.А. Реута, С.С. Рокотяна. - 2-е изд., перераб. и доп. - М.: Энергия, 1980. - 295 с.

3. Проектирование линий электропередачи сверхвысокого напряжения / Под ред. Г.Н. Александрова, Л.Л. Петерсона. - Л.: Энергоатомиздат, 1983. - 368 с.

4. Справочник по электрическим сетям 0,4-35 кВ и 110-1150 кВ. Т. 2 / Под ред. И.Т. Горюнова, А.А. Любимова. - М.: Папирус ПРО, 2003. - 622 с.

5. Правила устройства электроустановок. Раздел 2. Передача электроэнергии. Главы 2.4, 2.5. - М.: Изд-во НЦ ЭНАС, 2003. - 160 с.

6. Бессонов Л.А. Теоретические основы электротехники. Электромагнитное поле: учебник. - М.: Гардарики, 2003. - 317 с.

7. Бермант А.Ф. Краткий курс математического анализа / А.Ф. Бермант, И.Г. Араманович. - М.: Лань, 2006. - 720 с.

Поступила 24.11.2006 г.

i Надоели баннеры? Вы всегда можете отключить рекламу.