CC BY
21
УДК 621.316.925
ПОЛЕВАЯ И ЦЕПНАЯ МОДЕЛИ ВОЛНОВЫХ ПРОЦЕССОВ В ПРОТЯЖЕННОМ ЗАЗЕМЛИТЕЛЕ*
СЛЫШАЛОВ В.К., д-р техн. наук; ГОЛОВ П.В. асп., КИСЕЛЕВА Ю.А., ТИМОФЕЕВА И.В., инженеры
Рассмотрены вопросы математического моделирования волновых процессов в протяженных зазем-лителях, проложенных в однородной среде, при прохождении по ним синусоидального тока. Моделирование осуществляется на основе описания и расчета волнового поля заземлителя и построения модели в форме цепи с распределенными параметрами, адекватной полевым процессам. Полученные результаты обобщены на случай слоистой среды при прокладке заземлителя перпендикулярно слоям, как это имеет место в системах скважинных зеземлителей.
Ключевые слова: математическое моделирование волновых процессов, электромагнитное поле, заземли-
тель
FIELD AND CHAIN MODELS OF WAVE PROCESSES IN DISTANT GROUND
SLYSHALOV V.K., Ph.D., GOLOV P.V., postgraduate, KISELYOVA Yu.A., eng., TIMOFEEVA I.V., eng.
The article deals with the wave processes simulation in distant ground made in uniform environment while carrying sinusoidal current. The modeling is carried out on the basis of wave field ground description and calculation and model construction in the form of a circuit with distributed parameters adequate to field processes. The results are summed up in case of a flaky environment during laying the ground perpendicular to the layers as it happens in borehole ground systems.
Key words: wave processes simulation, electromagnetic field, ground.
Введение. Задача расчета волнового поля протяженных естественных заземлителей (оболочки кабелей, водопроводные и обсадные трубы, арматура, рельсовые пути и другие т.п. объекты) является базовой при оценке их эксплуатационных характеристик: входного сопротивления, рабочей длины, токов и напряжений при подключении к системам искусственных заземлителей. Интерес представляют и волновые параметры: постоянные распространения, волновое сопротивление, скорость распространения волны, коэффициенты отражения и преломления - знание которых необходимо для разработки и совершенствования дистанционных волновых методов определения мест повреждения кабельных линий, при исследовании распределений токов и напряжений в заземленных дискретно объектах, таких, как нефте- и газопроводы, проложенные на поверхности земли, кабели и трубы в туннелях и на эстакадах, защитные тросы ЛЭП. Расчет волнового поля и перечисленных параметров и характеристик обязателен и при проектировании искусственных протяженных заземлителей, среди которых следует выделить скважинные зазем-лители [1], непосредственно относящиеся к тематике настоящей работы.
Постановку общей задачи расчета волновых параметров протяженного заземлителя иллюстрирует рис.1, на котором схематично показан путь протекания тока молнии от молниеотвода до скважинного зазем-лителя, обеспечивающего растекание тока в земле. Все обозначенные на рис. 1 участки заземлителя, исключая участок 3, характеризуемый сосредоточенным сопротивлением, рассматриваем как протяженные, выполненные из цилиндрических стержней и труб.
стержневой протяженный заземлитель, проложенный в воздухе вблизи поверхности земли; 3 - участок перехода из воздуха в землю; 4 - аналогичен 2, но проложен в земле; 5 -скважинный заземлитель, 5' и 5" - его участки, находящиеся в слоях земли с разными электрическими свойствами; -концевое сосредоточенное сопротивление
В данной статье исследуется электромагнитное поле прямого протяженного заземлителя в форме цилиндрической трубы (рис. 2), проложенной в неограниченном проводящем пространстве (участки 5' и 5" на рис. 1). Как частный случай задачи рассматривается поле стержневого цилиндрического заземлителя в неограниченной проводящей среде и воздухе (участок 1 на рис. 1). Из решений для векторов поля устанавливаются параметры цепи с распределенными параметрами.
Электромагнитное поле и волновые параметры трубы и стержня. Полагаем, что труба проложена в неограниченном пространстве, заполненном средой, электрические свойства которой характеризу-
* Изложенные в статье результаты получены при финансовой поддержке в форме гранта Министерства образования Российской Федерации (шифр гранта Т02-01.5-1186).
ются удельным сопротивлением р2, диэлектрической и магнитной проницаемостью 8а2 = 8280, ца2 = цо. Внутренняя область трубы заполнена средой с параметрами р2, 8а2, Ца2 = цо. По заземлителю протекает синусоидальный ток, амплитуда и фаза которого изменяется по длине трубы за счет оттока в среду, и поэтому комплексное значение тока является функцией координаты х: I = 1(х).
|1(х)+1з(х)=1(х)
Рис. 2. Координатная система и геометрические характеристики заземлителя: 1, 2, 3 - расчетные области
Описанный процесс, как это следует из уравнений Максвелла [2], является волновым, причем все компоненты комплексов напряженностей электрического и магнитного полей Е,Н удовлетворяют уравнению Гельмгольца
Ду, + = 0, (1)
где у, - любая из составляющих напряженностей в среде , = 1,2,3 (рис. 2);
kf - га2 gaiца(- = га2 | g
ai — J — |Hai-тар,- 1
Классическая методика решения уравнения (1) [2-5] предполагает экспоненциальное изменение у,- по координате х с показателем экспоненты ±ух, у = а + _/'Р -постоянная распространения, инвариантная относительно номера среды ,; а и р - коэффициенты затухания и фазы; знаки «-», «+» соответствуют прямой и обратной волнам.
Полагая, что заземлитель имеет неограниченную длину, будем рассматривать в полевой модели только движение прямой волны, а процессы отражения и преломления волн, обусловленные конечной длиной участков заземлителя и послойным изменением удельного сопротивления среды, разберем на основе цепной модели.
С учетом указанного выше допущения, формаль-
но означающего для прямой волны, что
д_
dx
-Y.
уравнение (1) относительно составляющей Ёх, в координатной системе рис. 2 принимает вид
1d r dr
(
dix, dr
\
■mfÈxi -
m2 - k2
■y
(2)
Под Ёх, понимается комплексная амплитуда продольной составляющей напряженности Е,-, т.е.
Ех(г, х) = Ех1 (г, 0)е/УЕ(г,0) • е. (3)
В аналогичной форме представлены и другие составляющие векторов поля, причем из уравнений Максвелла для них имеем, согласно [6], формулы:
д
Е i
дх
mr
дЕЕх,
dr
\
■ H - — i raga Ёж
' ui j 2 a-
m,2 dr
(4)
Общая форма решений уравнения (2) известна, поэтому, опуская промежуточные выкладки, запишем выражения для Ёх, Ёг, Ни в областях 1, 2 (рис. 2): поле в заземлителе Г < г < г2) -
ЕЕ xi = ^i^o(kir ) + B1Ho(k1r ),
Er 1 -1 d?1 J0 (kir ) + f dB1 H0 (kir ),
k1 dx k1 dx
Hu1 = —г-1—(V ) + S1H0 (k1r )]■ k1P1
поле в грунте (r2 < r < м) -
(5)
Ex 2 = C2H0(m2r ),
ЕЕr 2 =
1 dC2
m dx
H0 (m2r ),
Hu2 = —"
1
-C2H0 (m2r ).
(6)
m2P2
В выражениях (5), (6) и0(кг), Н0(тг) - функции Бесселя и Ханкеля. При записи (5), (6) учтено, что в теле заземлителя, где р1 ~ 10-7 Омм, допустимо принять [6]:
2 2 2 m( - k( + у2
• k2 - — j
P1
g1a = g1a — J'
1
4-
1
(7)
rap1 rap1
Для грунта при рз = 102-105 Омм и £2 ~ 10 [7] второе из соотношений (7) выполняется вплоть до частот ю ~ 106-107 1/с, т.е. и для импульсных токов молний, условно характеризуемых на фронте разряда частотами ~ 105 1/с. Следовательно, для грунта
2 2 2 m2 = k2 + у2
- j ^ + у2, Р2
(8)
Поскольку некоторые участки протяженных за-землителей расположены в воздушной среде, одним из вариантов (6) являются решения уравнения (2) при Р2 и £2а ~ £0. В этом случае:
m| -ra2g0^0 + у2, ga -g0;
Ex 2 = C2H0(m2r ), Er 2 =— dri ®H0 (m2r ),
m2 dx
(9)
(10)
Hu2 =— j^o C2H0(m2r). m2
Во внутренней области трубы (0 < r < Г1)
Ex3 = D3J0(m2r ), E з =— ^г3 J0 (m3r),
m3 dx
Hu3 = -
1
m3P3
D3J0 (mir).
(11)
В интересующем нас случае, когда труба является скважинным заземлителем, внутренняя область заполняется глиняным раствором [1] с рз = 50-100 Омм, т.е. электромагнитный параметр тз определяется выражением (8), в котором р2 следует заменить на р3, и в силу инвариантности постоянной у вычисляется через т2по формуле
"Р2
— (12)
2 2 P3
m3 = m2 - jюц^-^3-
Р3 -Р2
Граничные условия задаются на поверхностях r = r1 и r = r2 в форме равенств:
Exi( ri) = Ex3(ri); Нui(ri) = Hu3(ri);
Exite) = Ex2(r2); H^) = H^).
(i3)
Условие «на бесконечности» реализуется выбором функции Ханкеля H0(m2r) = H§\m2r), при котором для положительного знака у корня m2 = +yjk2 + у2 обеспечиваются условие Jm m2 > 0
и, соответственно, равенство lim Ex2(m2r) = 0 . За-
r ^^
данным также является ток в заземлителе, что дает возможность записать еще одно условие однозначности, используя закон полного тока
/ = i + 3 =
2п
= [ H^foMи =--C2H0 (m2r2)2nr2
m2P3
(i4)
Таким образом, записано пять условий однозначности - (13), (14), что при четырех константах интегрирования позволяет определить электромагнитный параметр т2, а через него - постоянную распространения у2 = т2 - к| .
Для определения т2 воспользуемся получающимися из (13) равенствами:
Е х1(1) = Ех 3(г1) ; Е х1 (г2) = Ех 2(^2)
подставляя в которые граничные значения напряжен-ностей, вычисляемые по (5), (6), (11), записываем уравнения
в
J0(kiri) + ^г1 ^(kn)
л
J0(m3r,)
i
kiPi
J0 (Vi) + Bl H0 (kiri)
i
m3P3
J0 (m3ri)
(i5)
в
J0(kir2) + А H^V) "i
H0(m2r2)
i
kiPi
J0 (V) + B H0 (V) Ai
H0(m2r2)' (16)
m2P3
образующие совместно с равенством (12) систему,
В
позволяющую найти отношение — и параметр т2. По
Ai
в
отношению находится и ток /3 в сердечнике трубы. Для этого опять воспользуемся законом полного тока. Имеем:
Hui(ri)2nri = /3 = i
kiPi
Ai
B,
J0 (V) + A H0 (kiri) Ai
Hui(r2)2nr2 = / = i
kiPi
Ai
Bi
J0(V2) + ^H0(kir2)
2nri;
2nr2,
Для отношения токов получаем выражение
Bi
J0 (kiri) + AL H0 (kiri) _A_
Bi
(i7)
1 ^ 4 С^) + В "0 (к/2)
А,
При расчете поля скважинного заземлителя величина параметра т2 может быть оценена по соотношению т2 ~ к2, значение внешнего радиуса трубы Г2 не превышает, согласно [1], 0,3 м. Вычисления показывают, что выполняются неравенства |т2г2| << 1,
\т3г1\ << 1, в силу чего является допустимым использование при решении уравнений (15) и (16) приближенных формул [8]:
^31) „ 2 "0(т2Г2) я т2г21п Хт2к, (18)
J0 (m3ri) m3ri H0 (m2r2)
2 j
где у' = i,78i07.
Решая (i5) с учетом (i8), получаем J0(kiri)
Bi = J0 (kiri)
z+
J0 (kiri)
Ai H0(kiri) z + Hski)'
(i9)
H0 (kiri)
,P3 1
где С = 2- .
Р1 к1г1
При исследованиях полей реальных заземлите-лей выражение (19) может быть существенно упрощено. Поясним этот момент, полагая, что по заземлите-лю протекает ток промышленной частоты. В этом случае для медного заземлителя (р1 ~ 1,72-10"8 Ом м, щ =1)получаем
kiM = • I-
j
Pi
= i5i,46• e~j45 i/м,
для железного (pi ~ iü"' kiM = 1,987-103е'45°1/м.
Омм, щ = 1000)
Приняв рз ~ См = 3,839 108еj4S
м, находим
- 102 Омм, г^ = 0,2 ', Сж = 0,503107е '45°. При возрастании частоты до 5000 1/с значения См, Сж уменьшаются всего на порядок, т.е. и при импульсных токах для любых металлических заземлите-лей имеем \С\>>1.
Поскольку для достаточно больших по модулю значений
формула (19) преобразуется к виду В ^(к^) =_ ^1(к1Г-|) А" Н0(к1г1) Н1(к1г1>.
(20)
3
При подстановке (20) в (17) получаем -3
'i
i 0 , т.е.
допустимо принять, что ток и поле во внутренней области скважинного заземлителя при заполнении ее любой слабопроводящей средой отсутствуют.
Подставляя (20) в (16), получаем, с учетом второй из формул (18), уравнение для определения электромагнитного параметра т2
к1Р1 ^/2) _ НШ^
2ПГ2 Ji(kf2) - ^г?1-Hi(kir2)
m|P2|n Y m2r2
(21)
2п 2 \
Н^кг)"
Левая часть этого уравнения дает комплексное сопротивление трубы на единицу длины
¿ОГ
k P J0(k/2) -Нгтггнкгг)
, kiPi___Hi(kiri)_
|nri Ji(kir2) -^^Hi(kir2) ' n i 2 Hi(kiri) n i 2У
(22)
в чем можно убедиться, вводя в (22) функции Кельвина: Jn (|k|r,V-7) = bern (\k\r,) + jbein (|k|r,);
Hn (( ) = hern (( r,) + jhein flk|r,),
где n - 0, i; i - i, 2. После этого получаемое выражение легко преобразуется в формулу, приведенную в справочнике [9].
Выполняя в (2i) предельный переход при ri ^ 0 , получаем уравнение для определения параметра m2 стержневого заземлителя, проложенного в неограниченной среде:
kiPi Jj(kia) m2"
____ m2P2 |n Y m2a
Ina Ji(kia) In 2j '
(2з)
где a - радиус стержня; ние на единицу длины
= kiPi J0<kia)
^ос -
его комплексное сопротивле-
(24)
2па ^1(к1а)
Аналогичным образом формула (22) преобразуется в (24) и при достаточно больших значениях частоты, когда эквивалентная глубина проникновения
А-
2Pi
<< r2 - ri
Трансцендентные уравнения (2i), (23) решаются методом последовательных приближений [2, 6].
Расчеты были выполнены для ряда частот от f = 50 1/с до f = 20000 1/с, причем была установлена возможность использования при /1 > 8000 1/с формул
а-^Т2 (a/HP!^)-2 -
Р-^М. ^Ji7(P2Sa2ra)-2 +
il =
1P2 ' 1P2 '
(15)
соответствующих полю плоской электромагнитной волны, распространяющейся в среде с конечной проводимостью [13]. На подробном анализе этих результатов в данной статье мы не останавливаемся.
В табл. 1, 2 приведены некоторые результаты, позволяющие конкретизировать значения физических характеристик стержневых заземлителей. Результаты даны для стержня Г2= а = 5-10"3 м, материал - медь и железо; удельное сопротивление среды - в первом варианте Р2= 102 Омм, во втором р2 = 103 Омм, расчетные значения частоты - / = 50 1/с, /2 = 5000 1/с.
Таблица 1. Стержневой заземлитель в грунте. Частота тока в заземлителе Ъ = 50 1/с
Медь PI=1,7I-10-8 Омм Mi=i p2=i02 Омм mi = i.075i0-3-- j-3.I78-i0-5 Y = i.6i3i0-3+ + j-i.I0I-i0-3
p2=i03 Омм m2 = 3.I48i0-4-- j-7.97-i0-6 Y = 5.045 i0-4+ + j-3.86i-i0-4
Железо p2=i02 Омм mi = 5.937i0-3+ + j-i.9I8-i0-3 Y = 6.039 i0-3+ + ji.III-i0-3
|Ji=i000 p2=i03 Омм mi = i.776i0-3+ + j-5.89i-i0-4 Y = i.8i0i0-3+ + j-6.868-i0-4
Таблица I. Стержневой заземлитель в грунте. Частота тока в заземлителе f2 = 5000 1/c
Медь pi=i,7Ii0-8 Омм Mi=i р2 = i0! Омм m2 = I.08Ii0-3+ + j-6.7I-i0-4 Y = 0.0i4+ + j0.0i4
р2 = i03 Омм mi = 6.I59i0-4+ + ji.057 i0-4 Y = 4.477 i0-3+ + j-4.438-i0-3
Железо pi=i0-7 Омм Mi=1000 р2 = i02 Омм m2 = 0.0I0+ + j-6.570-i0-3 Y = 0.0I3+ + j0.0i4
р2 = i03 Омм mi = 5.804 i0-3+ + jI.00Ii0-3 Y = 7.038 i0-3+ + j-4.456-i0-3
С аналогичными результатами для кабельной оболочки (г1 = 2-10"2 м, г2 = 2,3-10"2 м, р1 = 10-7 Омм, Р2 = 102 Омм, £2 = 10, / = 50 1/с) можно ознакомиться по [6].
Для участка заземлителя в воздухе (участок 1 на рис. 1) в уравнении (23) изменяется его правая часть
1
путем замены активной проводимости — на емкост-
Рз
ную j®е0 . В итоге получаем
kiPi Jp(kia) = -_m|_ln Y'mia Ina Ji(kia) In^ra I j
(26)
Результаты расчета для этого случая приведены в табл. 3; геометрические и электрические параметры стержня оставлены прежними.
Таблица 3. Стержневой заземлитель в воздухе. Частота тока в заземлителе f1 = 50 1/c, f2 = 5000 1/c
ние
Медь pi=i,72i0-8 Омм Mi=i fi m2 = 3.042 i0-7+ + j-3.i44-i0-7 Y = 9.073 i0-8+ + j-i.054-i0-6
f2 m2 = 3.949 i0-6+ + j-8.585-i0-6 Y = 3.228 i0-7+ + j-i.050-i0-4
Железо pi=i0-7 Омм |Ji=i000 fi m2 = i.027 i0-6+ + j-2.279-i0-6 Y = 9.452 i0-7+ + j-2.476-i0-6
f2 m2 = 3.439 i0-5+ + j-7.950-i0-5 Y = 2.i25i0-5+ + j-i.287-i0-4
После определения параметра т2 вычисляются постоянные интегрирования А1, В1, С2. По (14) получаем
С2 =-/
m2P2 ~ j m2P2 j 2nr2H0 (m2r2) 4
кроме того, по (13), (20) и (5) имеем
A^V)
1 -
Ji(kfi) ^0(^/2)
= С2Н)(т2Г2).
(27)
(28)
Оставшаяся константа В? находится по соотношению (20). Эти же формулы при Г1 ^ 0, = а дают постоянные интегрирования для задачи о поле стержня. На основе выполненного расчета электромагнитного поля заземлителя определяются его интегральные параметры: волновое сопротивление ¿в, переходное сопротивление ¿п и сопротивление заземления (входное сопротивление) 2вх, являющиеся совместно с параметрами т2 и у = а + ¡в базовыми при построении цепной модели, адекватной по своим характеристикам полевым процессам в заземлителе и окружающей его среде.
Интегральные параметры заземлителей. Полевая и цепная модели. Дадим определение напряжениям, характеризующим работу заземлителя, и запишем для них расчетные формулы. Для пояснений используем рис. 3, на котором схематично изображен трехфазный силовой кабель в режиме короткого замыкания одной из его жил на проводящую оболочку. Полагаем, что кабель имеет значительную длину, позволяющую считать его бесконечным и погруженным в проводящую среду неограниченного объема.
х
иС^Х)^^
Ü(1)(0)-U(1) (x)
P2, S2
Рис. 3. Электрические характеристики оболочки силового кабеля при работе его в качестве заземлителя
Для заземлителя, с которого в проводящую среду стекает ток, необходимо в общем случае различать
два напряжения: и (1)(х) и и (2)(х) (рис. 3), обусловленных продольным и поперечным токами. Напряже-
U (i)(x), измеряемое
вдоль заземлителя, является
характеристикой заземлителя как токоведущего элемента и в случае оболочки кабеля, например, определяет в различных сечениях величины напряжений между жилами кабеля и оболочкой, а значит, и емкостные токи в промежутках жила-оболочка [6]. Для вычисления и (1)( х) имеем формулу
ад ад
и (1)( х) = | Ех1 (£, Г2) С1 £ = | Ех1 (0, Г2) е£ =
= - Exi (°> Г2 К
Y
yx
= U (i)(0)e
-yx
(29)
которая с учетом соотношения
— = Exi (x-r2 ) = Exi (x-r2 )
-0 2nr2HHui(x, Г2) /(x) может быть записана в виде
U (i)( x) = /(0)e-yx-° = /(x)=.
Y Y
Напряжение U(2)(x), измеряемое в поперечном
относительно заземлителя направлении (координата r на Рис.3) между заземлителем и виртуальной цилиндрической поверхностью S3 радиуса г» является характеристикой процесса стекания тока в среду. Оно же является напряжением прямой волны, сопровождающей волну тока в заземлителе. Для вычисления
U(2)(x) имеем формулу
гад <*> .
U(2)(x) = lim fEEr2(x,г)dr = f — (-Y)C2H0(m2r)dr =
Г»^ад r2 Г2 m2
= C2x) H0 (т2Г2).
Подставляя в (30) выражение (27), получаем
(30)
= Y m2P2H0 (т2Г2 ) v) = Y—0 /,
(3i)
и(2)(х) = -У7Н01"2'2 Кх) = ^/(х).
т| 2пг2Н0 (т2г2) т2
Таким образом, напряжения и (1)( х), и (2)( х) различны, что обусловлено изменениями во времени магнитного поля, созданного токами в заземлителе и грунте. Действительно, в соответствии с законом электромагнитной индукции можем записать (см. рис. 3) соотношение
и(1)(х) - и(2)(х) =
% Гад
jro lim ff^0^2Hu2(x,г)d%dr = -jroö(x).
(32)
Гад^ад xr2
P2
(33)
Интеграл в (32) легко вычисляется, и в итоге с учетом равенств
2 =_ т2Р2Н0 (т2Г2 ) к = _ I — 2пг2Н0 (т2г2) ' 3 находим
¿0к22 . _ } гоФ (х) = =Т/(х). ут2
х
x
Этот же результат получаем, вычисляя и(1)(х) - и(2)(х) по формулам (29), (31).
Волновое сопротивление заземлителя, равное по определению отношению напряжения и тока прямой волны, вычисляется по формуле
и (2)( х) ч^о
Z
/(*)
m2
(34)
и в любом сечении заземлителя равно входному сопротивлению или сопротивлению заземления
и(2)(0) и (2)(х) уТо
— /(0)
/(х)
m.
2 S Zb .
(35)
Переходное сопротивление, характеризующее стекание тока с заземлителя в грунт, следует определить как погонный параметр, равный отношению напряжения и(2)(х) к току утечки с единицы длины за-
землителя. Вычисляя его, получаем с учетом формул (6), (27) значение
/о(х) = 2п • Г2 —ЕГ2(х,г2) = У • /(х) = .
Рз дх
Zn
Соответственно
U (2)( х) U (2)( х)
Y
2
m22
(36)
— /о( х) у/( х) т.е. качество заземления зависит не только от удельного сопротивления грунта, но и от материала зазем-лителя и его геометрических параметров.
Подставляя в (33) параметр т2 из соотношения (36), получаем с учетом равенства
Н0(т2г2) = Н1(т2г2) уравнение для переходного сопротивления Тп в критериальной форме П Но(П2)
2п Н,(П2)
(37)
где П =;
Р2
П2 = Г2
Критерий П, содержащий только заданные электрические и геометрические характеристики заземлителя и среды, может быть назван критерием заземлителя. Критерий П2 определяется при решении уравнения (37) и является критерием волнового процесса работы заземлителя. По найденному значению П2 определяются Zn , Zb , Y, тг, U(1)(х),
U(2)(х), Ф(х) и вектора напряженностей электрического и магнитного полей. В силу известных теорем о подобии [10] при выполнении для группы заземлителей условия П1 = idem волновые процессы в них будут подобны, т.е. перечисленные выше характеристики и параметры всех заземлителей будут вычисляться исходя из единого для всей группы значения критерия П2 .
Если воспользоваться приближенной формулой (18), можно преобразовать (37) в следующее удобное для практических расчетов уравнение:
z n=
J
Z 0/
■ = РЗ In. 2п
1.123 J
Z 0,
(38)
В табл. 4, 5 приведены результаты расчетов Тп, Твв для медного и железного стержней при рассмотренных выше значениях Г2, Р2, f и данных табл. 1, 2. Расчет выполнен по соотношению (36).
Таблица 4. Стержневой заземлитель в грунте. Частота тока в заземлителе 11 = 50 1/с
Медь р1=1,7210-8 Омм М1=1 Р2=102 Омм Zn = 196.607-е i7'447" Zb = 0.396-е Н4152°
Р2=103 Омм Zn = 2155-е №766° Zb = 1.369-е Н4193°
Железо Р1=10-7 Омм |J1=1000 Р2=102 Омм Zn = 168.157-е №832° Zb = 1.082-е i27032°
Р2=103 Омм Zn = 1872-е №103° Zb = 3.625-е i26'876"
Таблица 5. Стержневой заземлитель в грунте. Частота тока в заземлителе f2 = 5000 1/c
Медь р1=1,7210-8 Омм |1=1 Р2=102 Омм Zn = 184.737-е №225' Zb = 3.684-е i5a945°
Р2=103 Омм Zn = 2037-е i5'619" Zb = 12.842-е i50'367"
Железо р1=10-7 Омм |1=1000 Р2=102 Омм Zn = 149.270-е i7'638" Zb = 4.045е i38-899i
Р2=103 Омм Zn = 1684-е №723° Zb = 14.026-е i39063°
Для этих же стержней в воздухе значения Тп, Тв приведены в табл. 6.
Таблица 6. Стержневой заземлитель в воздухе. Частота тока в заземлителе 1 = 50 1/с, 12 = 5000 1/с
Медь p1=1,7210-8 Омм J1=1 f1 Zn = 1.149-109-е -i'87-805° Zb = 1.216103е -iZ723°
f2 Zn = 9.725 1 06е -i88546° Zb = 1.021103еi'1-278'
Железо Pi=10-7 Омм |1=1000 f1 Zn = 1.049 1 09е -i88'676" Zb = 2.779 103е -im574°
f2 Zn = 8.457 106е -i88416° Zb = 1.103103е -i7J92"
Для практических расчетов протяженных заземлителей в [11] предлагается использовать упрощенный вариант уравнения (38), введенный в рассмотрение Е.Д. Зунде [12]. В наших обозначениях это уравнение имеет вид
Rn = рЗ in. n 2п
1.12
Zn
(39)
где ^ = КеТп ; То = \То\,
и предполагает наличие у переходного сопротивления только вещественной составляющей, что допустимо в оценочных расчетах для заземлителя в грунте и совершенно неприемлемо для заземлителя в воздухе, где ТП (см. табл. 6) имеет практически чисто емкостный характер.
Решения уравнения Гельмгольца (1) для напря-женностей электрического и магнитного полей (5), (6), (11); уравнения (21), (23) для расчета постоянной распространения и выводимые из решений для поля формулы (29), (3о), (34), (35), (36), по которым опре-
деляются интегральные характеристики и параметры волновых процессов в заземлителе и среде, образуют полевую модель работы протяженного неэквипотенциального заземлителя в проводящей или изолирующей среде при синусоидальном токе в нем. Ключевым для этой модели является расчет постоянной распространения по уравнениям (21), (23) или уравнениям (37), (38). На основе этой модели целесообразно рассматривать те задачи, в которых изучение режимов работы заземлителя сопрягается с исследованиями полевых процессов (искрообразование и пробой в сыпучих грунтах, разогрев и спекание грунта, корони-рование заземлителя и т.д.). Если же достаточным
является знание распределений и(х), /(х), /3(х), в
том числе и для отраженных и преломленных волн, целесообразно воспользоваться цепной моделью, под которой понимается совокупность уравнений и формул для определения перечисленных интегральных характеристик работы заземлителя. Поскольку для заземлителя имеем
и(2)(х) = и(2)(0}е~ух, /(х) = /(0)е~ух =
и(2)(х)
¿в в
(40)
одной из возможных моделей является однородная цепь с распределенными параметрами, уравнения которой имеют своими решениями зависимости (40). Конфигурация такой цепи, допускающая и физическую реализацию, показана на рис. 4.
1(0)
£ом
2ом
1_Г'(0)
'(х)
£ом
Уом
У|2|(0)
1 Уом
Ук
Рис. 4. Конфигурация цепной модели заземлителя
Уравнения установившегося режима при синусоидальных токах и напряжениях записываются для нее в форме [3]
ди(2\х)
бх ¿К х) бх
=1ом1( х ь = X оми (2)( х)
(41)
где ¿оом , Уом - первичные погонные параметры.
Вводя волновые (вторичные) параметры - постоянную распространения YМ и волновое сопротивление
¿ВМ:
У м
-^¿ОМ УОМ ;
¿ВМ ГОМ
У
ОМ
(42)
можем записать решения системы (41) при отсутствии отраженных волн в форме
и (2)(х
" (43)
и(2)(х) - и(2)(0)е~™х, /(х) - /(0)е~™х - -
-ВМ
Условиями адекватности процессов в заземлите-ле и модели являются, согласно (40), (43), равенства
Ум - У , 1вм - ¿в , (44)
позволяющие определить первичные параметры модели через параметры заземлителя и среды. Исходными в этом расчете являются сопротивление зазем-лителя ¿о (формулы (22)-(24)) и параметр т2, вычисленный по уравнениям (21), (23), или переходное со-
противление ¿п, определяемое по уравнению (38). Используя совместно с (42), (44) формулы (34), (36), находим
1
(45)
22 2 г -у22 У -_
¿ОМ - 2 Ь-0 - ' Ь-П ' уОМ - ^ - ^ т| ¿0 ¿п
Применение в уравнениях (41) этих параметров обеспечивает тождественность распределений по длине заземлителя токов и напряжений в исходной полевой задаче и ее цепной модели. Соответственно, недопустимым является нарушение этих условий, например, использование в уравнениях (41) сопротивления ¿о вместо ¿ом, как это сделано в работе [11 ], без оценки возникающей при этом погрешности. Поскольку ключевым расчетным моментом при использовании уравнений (41) является, как и в полевой модели, определение электромагнитного параметра т2 или переходного сопротивления ¿п, то при рассмотрении бесконечно длинного заземлителя в однородной среде следует предпочесть полевую модель. Цепная модель имеет преимущества перед полевой при исследовании волновых процессов в заземлителях, проложенных в многослойных средах или средах с непрерывным изменением удельного сопротивления р-р(х), для которых имеется развитая методика составления и решения уравнений типа (41) [13]. Приведем основные результаты, относящиеся к расчету процессов в заземлителе по уравнениям (41) в случае его прокладки в многослойной среде (участок 5 на рис. 1).
Полагаем, что имеется / = 1,2, ... п слоев, для каждого из которых по изложенной выше методике определены волновые параметры у/, ¿в/. Задан ток на входе / и комплексное сопротивление К (рис. 1). Искомыми являются сопротивление заземления ¿вх и распределения тока и напряжения \(х), и(х) по длине заземлителя.
Напряжение и/ (х) и ток ii(x) (/ = 1,2, ... п) формируются при наложении прямой и обратной волн: и31 (х,) - йпр; + йоБР/ - А/е-У/х/ + А2е~У/К/;
А/
//(х,)- 1пр1 _ 1оБР1 -^е_'/х/
(46)
¿81
¿вВ/
Координаты х, в выражении (46) отсчитываются от начала слоя «/»; константы А?, А2/ связывают коэффициент отражения - комплексное число [3]
п
К.ВХ/+1 _ ¿в/ - 2у/е/
- ±^вХ/+1 -вI е ■ ■
отр/ - у-е ' (47)
ZВХ/+1 + ¿в/
где zвх|+1 - входное сопротивление заземлителя в «/+1» слое, равное при / = п сосредоточенному концевому сопротивлению тК;
А2/ - А/
^/+1 _ Z в / _ 21/1 /
- п,
отр/
А/.
-ВХ/+1
в
(48)
Подставляя (48) в (46), получаем для формулу
- ОдЮ - 7 1 + ПоТР/ (49)
1вх/ - /,(0) - 1вх/ 1 _ п^ , (49)
которая при / = 1 дает входное сопротивление зазем-лителя
1 + пг
(50)
¿вх -
'отр/
1_п
отр/
Коэффициент отражения в (5о) вычисляется путем последовательного, начиная с / = п, применения формулы (47), причем находим
_ Zк - ZBn e -2ln 'n
Пптрп ---п—e
'ОТРп затем пОТРп-1
ZK + ZBn
ZBXn _ ZBn
1+n,
ОТРп
1 - n
ОТРп
ZBXn - ZBn-1 ^ 2Yn-ltn-1
BXn
-1 _ Z
BXn Bn-1
1+n,
Bn-1
ОТРп-1
1 - n
ОТРп-1
и т.д. до получения значений ПОТР1 и ТВХ .
Аналогичным образом, но в обратном порядке, начиная с тока на входе I = 11(0), последовательно при учете граничного условия 1,(И) = //+1(0)определяются постоянные Лц по формуле
А/ _
ZBili (0)
1 - п
ОТР/
постоянные A2/ по (48) и функции ii(x), Üi (x) согласно
выражениям (46).
Вопросы расчета волновых процессов в заземли-телях при р = р(х), а также процессов в заземлителях, проложенных в грунте или воздухе вблизи границы раздела этих сред, будут рассмотрены авторами в последующих публикациях.
Заключение
1. Предложенный в статье метод расчета волновых процессов в протяженных заземлителях, являющийся вполне строгим для заземлителей кругового и трубчатого сечений, допустимо применять как приближенный при расчете волновых параметров заземлителей типа уголка, полосы, швеллера и других практических сечений, используя в формуле (38) эквивалентный радиус изучаемого заземлителя.
2. Обобщение полученных решений на задачи с током в сердечнике i3 ï 0 позволяет рассмотреть волновые процессы в биметаллических проводах и за-землителях, окруженных цилиндрическим слоем грунта с измененными свойствами.
3. Представляется перспективным использование результатов, относящихся к полю и характеристикам заземлителей при синусоидальном токе, для исследования импульсных режимов заземления и, в частности, для расчета параметров и поля лидерного канала и канала главного разряда молнии.
4. Подобие волновых процессов в заземлите-лях, характеризуемое пропорциональностью волно-
вых параметров m2
(2) _
_ kmm(1) _■
Г (1)
!L_ m(1)
.(2) m2 ,
обеспечива-
ется
выполнением 2
критериального
условия
П1 _ Z°rL _ /dem .
P2
5. При физическом моделировании процессов в заземлителях и заземляющих системах, при исследовании этих процессов с помощью математических моделей в форме цепей с распределенными параметрами необходимым условием адекватности волновых процессов в реальных заземлителях и моделях является использование модельных погонных параметров Zoм , Уом, вычисляемых по формулам (45) и уравнению (38) или уравнениям (21), (23).
Список литературы
1. Максименко Н.Н., Агеев Г.Г. Расчет скважинных заземлителей // Электричество. - 1970. - № 6. - С. 44-47.
2. Зоммерфельд А. Электродинамика. - М.: ИЛ, 1958.
3. Гринберг Г.А., Бонштедт Б.Э. Основы точной теории волнового поля линии передачи // ЖТФ. - 1954. -T. XXIV. -Вып. 1. - С. 67-95.
4. Ломоносов В.Ю. Электромагнитное поле двухпроводной линии: Сб. науч. тр. ИЭИ. Вып. IV. - М.-Л.: ГЭИ, 1951. -С. 223-247.
5. Костенко М.В., Перельман Л.С., Шкарин Ю.П. Волновые процессы и электрические помехи в многопроводных линиях высокого напряжения. - М.: Энергия, 1973. - 272 с.
6. Слышалов В.К., Гречухина И.В. Расчет электромагнитного поля кабельной оболочки в режиме однофазного замыкания // Вестник ИГЭУ. - 2001. - Вып. 2. - С. 5-8.
7. Рябкова Е.А. Заземления в установках высокого напряжения. - М.: Энергия, 1978. - 224 с.
8. Янке Е., Эмде Ф., Леш Ф. Специальные функции, формулы, графики, таблицы. - М.: Наука, 1977. - 342 с.
9. Калантаров П.Л., Цейтлин Л.А. Расчет индуктивно-стей. Справочная книга. - Л.: Энергоатомиздат, 1986. - 488 с.
10. Веников В.А. Теория подобия и моделирование применительно к задачам электроэнергетики. - М.: Высш. шк., 1966. - 487 с.
11. Карякин Р.Н. Электромагнитные процессы в протяженных заземлителях в неоднородных структурах // Электричество. - 1996. - № 7. - С. 43-51.
12. Sunde E.D. Currents and Potentials Along Leaky Ground-Return Conductors // Electrical Engineering. - 1936. -December (12).
13. Шимони К. Теоретическая электротехника. М., 1964. - 773 с.
Слышалов Владимир Константинович,
ГОУВПО «Ивановский государственный энергетический университет имени В. И. Ленина», доктор технических наук, профессор кафедры электрических систем, телефон (4932) 26-99-21, e-mail: [email protected]
Голов Швел Валерьевич,
ГОУВПО «Ивановский государственный энергетический университет имени В. И. Ленина», аспирант кафедры электрических систем. телефон (4932) 26-99-21, e-mail: [email protected]
Киселева Юлия Александровна, ООО «РЗА-сервис»,
кандидат технических наук, ведущий инженер, e-mail: [email protected]
Тимофеева Ирина Владимировна,
ГОУВПО «Ивановский государственный энергетический университет имени В. И. Ленина», Инженер вычислительного центра электроэнергетического факультета, телефон (4932) 26-99-21.
