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

УДК 621.316.99

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

Канд. техн. наук ГЕРАСИМОВИЧ Д. А., магистр техн. наук ДЕРЮГИНА Е. А.

Белорусский национальный технический университет, ООО «Электротехническая компания "ЭКНИС"»

Проектирование и комплексная оценка заземляющих устройств электроустановок представляют собой серьезную научно-техническую проблему. Теоретические и экспериментальные исследования, выполненные в России, Китае, странах Европы и Северной Америки, охватывают широкий круг задач по оценке заземляющих устройств [1-6].

Строгие математические методы решения задачи об электрическом поле и основных параметрах заземлителей разработаны лишь для немногих частных случаев простых заземляющих устройств, например для некоторых тел вращения с осью, перпендикулярной поверхности земли: полусферы, вытянутого и сплюснутого эллипсоида, тора [2]. Применительно к сложным заземлителям широкое распространение получили метод наведенного потенциала [2], а также конечно-разностные методы [2, 4]. Параллельно развивались приближенные инженерные методы расчета - метод коэффициентов использования, метод обобщенных параметров и др. [2]. К сожалению, большинство указанных исследований основывается на модели потенциала постоянного тока. В реальных условиях в заземляющих устройствах протекают переменные токи промышленной частоты и импульсные токи, вызванные разрядом молнии. Для учета отмеченных факторов необходимо использовать математическую модель на основе уравнений Максвелла [7]:

Гт - dD - dB rot H = 5 +—; rot E =--; (1)

dt dt

div D = p; div B = 0; div 5 = 0,

где H и E - векторы напряженности магнитного и электрического полей, соответственно; 5 - вектор плотности тока проводимости; B - то же индукции магнитного поля; D - то же электрического смещения; p - плотность электрических зарядов.

Ниже рассматривается метод расчета параметров электромагнитного поля уединенного бесконечно длинного стержневого заземлителя трубчатого сечения (рис. 1) при стекании в землю переменного тока промышленной частоты.

Данная модель представляет собой двухслойную среду, состоящую из слоя 1 заземлителя с внутренним и внешним радиусами R1 и R2 соответственно и слоя 2 собственно земли. Модель заземлителя круглого сечения можно получить как частный случай при R1 ^ 0.

Рис. 1. Геометрическая модель уединенного стержневого заземлителя

При стекании переменного тока промышленной частоты плотность тока проводимости существенно превышает плотность тока смещения

(

с>>

дЮ

\

как для слоя заземлителя, так и для слоя земли. С учетом отме-

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

го* Нгт = сгт; го* Е1т = - ]&в1т;

Шу Д„ = 0; ШУ 5,.„ = 0;

(2)

Вгт т ; Сгт УгЕгт ,

где г = 1, 2 - индекс номера слоя; ц,г-, уг- - магнитная проницаемость и электрическая проводимость г-го слоя соответственно; ю - угловая частота переменного тока.

Дифференциальные уравнения (2) для каждого слоя можно разрешить

относительно любой из компонент поля (напряженности магнитного Н

или электрического Е полей). Для принятой цилиндрической системы координат с учетом симметрии системы относительно угловой координаты система (2) сводится к уравнению [8]

+ <3^=а2к (г, г,, г = 1, 2, (3)

г дг дг2 &

где а2 = jюyiц,г-; фгт - может принимать значения комплексных амплитуд векторов напряженности электрического Егт и магнитного полей Нт или плотности тока с гт.

Решение уравнения (3) осуществляется относительно вектора плотности тока с,т -г,г) [9], который имеет две проекции: радиальную с¡тг и вертикальную ст. При этом уравнение (3) следует дополнить предельными условиями на поверхности земли и заземлителя, а также условиями на границе раздела сред заземлителя и земли.

Предельные условия для данной задачи неизвестны. Их нахождение для аналогичной задачи расчета электромагнитного поля в многослойной структуре плоского контактного соединения подробно рассмотрено в [9]. Для этого устанавливаются зависимости радиального и вертикального токов по координатам г и г. На основании этих зависимостей определяется распределение плотности тока. Рассмотрим расчет токораспределения отдельно для радиальной и вертикальной плотностей тока.

Расчет распределения радиальной составляющей плотности тока. Уравнение (3) относительно радиальной составляющей плотности тока принимает вид

+ + = а?с,т>, г), , = 1, 2. (4)

г дг дг2 дг

Для его решения установим зависимость радиального тока 5 ттг, приходящегося на единицу угловой координаты, от координаты г. При этом считаем известным ток 1т, стекающий в заземлитель. Его распределение (в данном случае вертикальная составляющая плотности тока) на поверхности при г = 0; г е [, Л2 ] определяется как для уединенного проводника трубчатого сечения [8]

с т = а То(х)Н0) (х1) - ^ (х1 )Н0)(х) и (п т (5)

С1тг (г,0) = 7а1-(2)'--27-Н1тф(R2,0), (5)

■Г0-х1)Н02) Ы-!0(Х2)Н02) -Х1)

где 10 -х), 10 -х) - функция Бесселя первого рода и ее первая производная,

соответственно; Н02)(х), Н02) -х) - функция Бесселя третьего рода (вторая функция Ганкеля) и ее первая производная соответственно; х = 7'а1г;

X] = 7'а1Л1; х2 = 7'а1 Я2; -Я2 0) = —-т--угловая составляющая векто-

2%Я2

ра напряженности магнитного поля на внешней поверхности трубчатого проводника.

На поверхности земли очевидно, что

С2тг -г,0) = 0, г е[»), -6)

а при г ^ да также

С1ш -г,да) = 0, г е [0,да); , = 1, 2. -7)

На основании (5)-(7) определяются зависимости радиального тока:

г г

(г) = 2п | (¿т (г, 0)- а1иг (г, да)) Мг = | то1)яг (г, 0)^г, г е [«, Я2 ];

12тг (г) = 1 т , г е[ К2, да) ,

а также тока, приходящегося на единицу угловой координаты:

' I т (Х1 )Н2)(х)-ВДНЧъ)

5. (г) =

гтг V '

1гтг (г)

2пг

2п«2 Мх^Н^Ы-^ЫН2^) I

, 1 = 1; г е[Я„ «2];

(9)

2пг'

1 = 2; г е [«2, да).

Полученные таким образом зависимости иллюстрирует рис. 2.

К К г

Рис. 2. Зависимость тока 5. от координаты г

Зная зависимости радиального тока, представим решение (4) как сумму одного из частных решений с«гтг (г, г), определяемого по радиальной составляющей тока в заземлителе, и невязочного решения АсЫг (г, г), обусловленного требованием непрерывности плотности тока с1тг (г, г) на границе раздела слоев земли и заземлителя [9]:

С1тг г) =С*тг г) + АС 1тг (г, ■ (10)

Нахождение решений вида (10) в [9] основано на применении рядов Фурье по тригонометрическим функциям. В данном случае аналогичные ряды получаются по функциям Бесселя. Исследования показали, что для рассматриваемой задачи ряды сходятся достаточно медленно, а их дифференцирование дает расходящийся ряд. Поэтому для определения решений сЧтг (г, г) и Ас 1тг (г, г) был использован метод конечных разностей. Рассматриваемая область при этом заменяется разностной сеткой (рис. 3) с узлами (к, Г). В направлении координаты г область г е (-да,0] ограничивается конечной областью г е[-Н ,0] из условия

Н > 4тах /—2—, 1 = 1, 2, (11)

V

т. е. в направлении оси г распределение электромагнитного поля рассматривается по крайней мере на четырех глубинах проникновения плоской электромагнитной волны. В остальной области г е (—да, - Н) с достаточной степенью точности с гт (г, г) = 0.

о

к = 0

к = Й1 (й2)

1 = 0 (щ)

/ = п1 («1 + п2)

Агг < гтг /, к—

/Аz¡ с гтг1- -1,т «гтгг/, С гтг/+1, к

' 1 < гт Г ' / к+1

г = К1 (К)

Рис. 3. Разностная сетка для г-го слоя (в скобках приведены обозначения для слоя земли)

Аналогично в направлении оси г ограничимся рассмотрением области г е [К, К3 ], где К3 выбирается с учетом (3) и (9) из условия

Кз >>,

2

|ЮУ 2^2

(12)

В остальной области г е (К3,да) распределение плотности тока рассчитывается по общеизвестным выражениям для распространения плоской электромагнитной волны.

Таким образом, узлам сетки (к, /) соответствуют точки сечения г1 = Агг/; гк = Аг{к. Шаги разбиения для каждого слоя по координатам г и г соответственно Агг = Кг+1—К- и Агг = —Н, где « и Иг - количество точек разбие-

пг Иг

ния слоя по координатам г и г соответственно, г = 1, 2 (рис. 3). Для обоих слоев: / = 0,1,...,п1,...,п1 + п2 и к = 0,1,...,Иг.

Для принятой разностной сетки первые и вторые производные по осям г и г в (4) аппроксимируются центральными разностными соотношениями [10] и уравнение (4) для внутренних узлов принимает вид

1 <гтг/+1,к — <гтг/—1,к , <гтг/+1,к — 2Сгтг/ ,к + <гтг/—1,к ,

2Аг

*гтг/ ,к+1 2Сгтг/,к + <гтг/ ,к—1 2 •

-= а с

Аг,

гтг/ ,к

Аг-2

, / = 0,1,...,«!,..., «1 + «2; (13)

к = 0,1,..., ^; г = 1,2.

г

г = - Н

Рассмотрим решение этого уравнения относительно составляющих à*imrl,k и Даimrl,k. Для нахождения частного решения à*imrl,k приведем (13)

к виду

ail,kà*imrl,k-1 + cil,kà*imrl,k + bil,kà*imrl,k+1 = dil,kà*imrl-l,k + fil,kà*imrl+l,k , (14)

где

l = 0,1,..., «!,..., n + n2; k = 0,1,..., hi,

Az2 2 2

aii,k =1; bii,k =1; cii,k = -2 тт- 2;

Ar

d,.,,, =

Az,

2 f 1 1 \ K„2 f

il,k i Ar

i 4

2r Ar-

; fil,k =

Az

Ar;

1 1

i 4

2r Ar

Очевидно, что уравнение (14) удовлетворяет условиям хорошей обусловленности относительно (7 w, к(jw к и (7*lmrl, к+1 [10]

\aa,k| + \ba,k| ^\ca,k| (15)

при любых соотношениях между шагами А^ и Azi.

Его решение осуществляется методом прогонки по слоям координаты z (к = 0,1,..., hi) в сочетании с методом итераций. Прогоночные коэффициенты в данном случае определяются выражениями:

L = -1 •

il,к ~ + Т '

Cil,k + Lil ,к-1

(16)

Fil,к - Kil,к-1 Cil ,к + Lil,k-1

Kil,k =

где ^П,к = ,кС*шг1 -1,к + кС*шг1+1,к •

По известным прогоночным коэффициентам рассчитываются величины плотности тока в узлах расчетной сетки

С*шг1 ,к = Ц1,кС*шг1 ,к+1 + КП к • (17)

Для использования (16), (17) необходимо знать значения плотности тока на внешних узлах, т. е. предельные условия. Кроме того, неизвестны значения плотности тока сЩтН-1 к и с*1шг1+1 к на смежных слоях по координате г по отношению к рассматриваемому слою. Нахождение этих величин осуществляется методом итераций. Для этого задается первоначальное (нулевое) приближение распределения плотности тока. Последующее математическое моделирование показало, что наименьшее количество итераций получается, если за нулевое приближение принять распределение плоской электромагнитной волны по координате г. Зависимость по координате г при этом принимается вида (9), т. е. нулевое приближение:

- 0 _

° * 'тг1, к ~

а,4 еа'гк Щ 1}Н(2)(Х,) - г}Н(2)(Х,) _

2пК2 ^ 1)Н(2)(х2) - ^(х2)Н(2)(Х1)' Х1 _ ^' ' _Г ^ а/ (18)

' т : _ .

2пГ

■е~''к, г _ 2; г1 е [Я2, <х>).

Зная нулевое приближение, рассчитываем новое приближение распределения плотности тока. При этом необходимо учесть, что для внешних узлов 1 _ 0; к _ 0,1,..., И' и 1 _ 0,1,..., п1,..., щ + п2; к = И; 6Птг1 ,к _ 0, а для

узлов 1 _ щ + п2; к _ 0,1,..., И; о.*'тг1к _о:°'тг1к. Расчет плотностей тока на

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

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

Д'тг1

и

| ё*тг1 (19)

Величина Д'тг1 сравнивается с известным током (9), а полученное токо-распределение уточняем по соотношению

_ Д'тг(Г) • (20)

а*'тг1 ,к ~ д а*'тг1 ,к ■

'тг1

Таким образом окончательно определяется первое приближение для плотности тока о*.'тг1 к. Это приближение сравнивают с нулевым и определяют максимальную величину невязки плотностей тока на двух смежных итерациях

1^1 I

"'"■<е, (21)

\^*'тг1,к ^*'тг1 ,к I

тах

(®*'тг1 ,к )

- Н

где етах - максимальная величина невязки; е - предварительно заданная достаточно малая величина; 0, 1 - верхний индекс соответствует номеру итерации.

Если условие (21) не выполняется, то рассчитывается следующее приближение на основании (16), (17), (19), (20). Итерационные циклы повторяются до выполнения условия (21).

После определения частного решения ЬЧтЛ к производим расчет невя-

зочного решения А<5'тг1к. Оно обусловлено тем, что частное решение

°*шг1 к не удовлетворяет условию непрерывности нормальных составляющих плотности тока на границе раздела слоев:

а1мг(Д2, г) = 62тг(Я,, 2). (22)

На основании (10) и (22) можно записать

Д°1тг (Я2 , 2) - Д°2тг (Я2 , = °*2тг (Я2 , - °*1тг (Я2 , 4 (23)

Для получения величин невязочных составляющих плотности тока на границе раздела используем дифференциальное уравнение непрерывности плотности тока сгт = 0 и равенство тангенциальных к границе раздела составляющих напряженности электрического поля, обусловленных величинами невязочных составляющих [11]. В принятой цилиндрической системе координат эти уравнения приводят к соотношению

+ дДа^Я 2)\ г^ДЯ^) , дД^тгЯ г)V

Г I г_дг ) = Г I Г_дг )

Ъ У 2 '

Так как проводимость заземлителя существенно превышает проводимость земли (у1 >> у2), то из (24) следует, что Дст1тг (Я2,2) >> Да2тг (Я2,2). Тогда (23) принимает вид:

Д^1тг С^ 2) = ^*2тг 2)-^*1тг С^ 2У;

(25)

ДС2тг (Я2, 2) = 0.

Таким образом, невязочные составляющие плотности тока появляются только в слое с хорошей проводимостью (в данном случае в заземлителе). Данный вывод хорошо согласуется с результатами [11].

Для нахождения невязочных составляющих (13) целесообразно привести к виду

а11 ,к Д°1тг1 -1,к + С11 ,к Д^1тгг ,к + Ь11 ,к Д°1тг1+1,к = ¿11 ,к ДЬ1тг1 ,к-1 + /и ,к ДЬ1тг1 ,к+1, (26)

I = 0, 1, ..., «1; к = 0, 1, ..., ¿1,

где

= ___L и =__1___1_ = А. + А. + 2

^ = 2г, Дг1 Дг12 ; Л1',к = 2г, Дг1 Дг12; ^ = Дг12 Д^2 1;

и = А. / = А.

^ = Д22; /11,к = Д22.

Его решение в целом аналогично решению (14). Только в данном случае метод прогонки применяется по слоям в направлении координаты г (I = 0, 1, ..., «1).

Прогоночные коэффициенты определяются по формулам:

L __Ь11 ,k_

М1,k _ L

С11 ,k + а11 ,kL1l+1,к

^ а11,кК11+1,к Ки ,к _■

С11 ,k + а11 ,kL1l+1,k где

^11,k _ d1l,kЛа 1тг1 ,k-1 + ^,kЛа1тг1,k+1,

а невязочные составляющие плотности тока

Ла 1mrl,k _ Ец,кЛа 1mrl-1,к + К11,k • (28)

Значения плотности тока на внешних узлах:

Ла 1тп, к _0;

(29)

Л^а1тгп1,к а*2тщ,к а*1тгп^к •

При к _ Ь1 имеем Лег 1тг1 к _ 0. Для расчета плотности тока в узлах к _ 0

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

В качестве первоначального (нулевого) приближения невязочных составляющих плотности тока целесообразно принять

Ла0 1к _Ла 1 к ^^ " ^1)) . (30)

1тг1 ,к 1тгп ,к л / / ^ г» \\ ^ '

1 8п(а1 (К2 - К1))

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

Расчет распределения вертикальной составляющей плотности тока. Для вертикальной составляющей плотности тока (3) имеет вид

+ + _ а2атг(г,г), / = 1, 2. (31)

г дг дг дг

Зависимость распределения вертикального тока по координате г (рис. 4) определяем по уже найденным распределениям горизонтального тока:

г

¡т(г) _ 1т 2пЯ2*2тг(^2,г)dг, ге[-И,0];

0

¡2тг (г) _ 0.

Рис. 4. Зависимость распределения вертикального тока /1тг по координате 7

Из условия равенства тангенциальных составляющих напряженности электрического поля и так как у1 >> у2 следует, что можно пренебречь

вертикальной составляющей плотности тока в земле.

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

а11, к -1, к + С11, к , к + Ь11 ,к +1, к = И11, к , к-1 + /и, к , к+1,

(33)

I - 0, 1, ..., «1; к - 0, 1, ..., ¿1,

где

1 1 11 2 2 2

- - -; К,, ------; с,,, -—- +—т + а'

11Л~ 2гг Ат1 Аг12 ' 1и ~ 2гг Аг1 А^2' Як " Аг12 А72 1'

и = -1- / = Й1',к = А72; ■/1',к = А^2.

Значения плотности тока на внешних узлах I - 0,1,..., «1, к - 0 (на основании (5))

• = а ^0 (Хг )Н02)' (Х1) - 10(Х1)Н02)(Хг) К

1 т&1 ,к J{X1 ~ г, '

■Г0(Х1)н02) (Х2) - ^(Х2)н02) (Х1)2п^2

(34)

Хг = /«Л; г1 е [ ^ ^2 ].

Прогоночные коэффициенты и плотности тока на каждой итерации определяем по выражениям (I - 0,1,..., п1):

т =_-Ь11 ,к_;

мг ,к = т '

с1г ,к + а1г ,кТ1г+1,к

(35)

^ ^1г,к - а1г,кК1г+1,к К1г,к

с1г ,к + а1г ,кТТи+1,к где ^1г, к = И11, к , к-1 + /и, к , к+1;

<1тг/, к = А/,к<1тг/—1, к + ^1/, к . (36)

Первоначальное (нулевое) приближение распределения вертикальной составляющей плотности тока принимаем

•0 ■ ^(Х/ )Н0)(Х1) — ^(Х1)Н0)(Х/) Ктг (гк) „ _ Г п Г>1 /-¡7^

<1тг/,к = 7а1-—'-—'---, , г К2 ]. (37)

■Г0(х1)н02) (Х2)—Т0(Х2)Н02) (Х1) 2пК2

Для каждого нового приближения определяем зависимость тока от координаты г

К2

11тгк

| 2пК2<1тг/,к^г. (38)

и корректируем плотности тока по соотношению

1 (г )

• ' _ 1тг У к/ • Поч

и1тг/, к ~ 1 и1тг/, к' ^^

1тгк

Итерации повторяем до тех пор, пока не выполнится условие вида (21) для вертикальных составляющих плотности тока

_

У]'7? ^ 1 УИ Х' I

< в. (40)

CT1mzl ,k °1mzl ,k

max

(^lmzl,k )

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

В Ы В О Д Ы

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

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

Л И Т Е Р А Т У Р А

1. Р я б к о в а, Е. Я. Заземления в установках высокого напряжения / Е. Я. Рябкова. -М.: Энергия, 1978. - 224 с.

2. Б у р г с д о р ф, В. В. Заземляющие устройства электроустановок / В. В. Бургсдорф, А. И. Якобс. - М.: Энергоатомиздат, 1987. - 400 с.

3. T a s k Force F0A of the IEEE/PES Substations West Coast Subcommittee. Current North American assessment and refurbishment practices of substations grounding systems // IEEE Trans. on Power Del. - 2005. - Vol. 20, № 3. - P. 1886-1889.

4. T a k a h a s h i, T. Calculation of earth resistance for deep-driven rod in multi-layer earth structure / T. Takahashi, T. Kawase // IEEE Trans. Power Del. -1990. - Vol. 6, № 1. -P. 608-614.

5. G u e m e s, J. A. Method for calculating the ground resistance of grounding grids using FEM / J.A. Guemes, F.E. Hernando // IEEE Trans. on Power Del. - 2004. - Vol. 19, № 2. -P. 595-600.

6. П о л е в а я и цепная модели волновых процессов в протяженном заземлителе / В. К. Слышалов [и др.] // Вестник ИГЭУ. - 2006. - Вып. 2. - С. 50-58.

7. Т а м м, И. Е. Основы теории электричества / И. Е. Тамм. - 9-е изд., испр. - М.: Наука. Гл. ред. физ.-мат. лит., 1976. - 616 с.

8. Б е с с о н о в, Л. А. Теоретические основы электротехники: в 3 ч. / Л. А. Бессонов // 4-е изд. - М.: Высш. шк., 1964. - 751 с.

9. Г е р а с и м о в и ч, А. Н. Обобщенная модель электромагнитных процессов в плоских многослойных контактных соединениях / А. Н. Герасимович, Д. А. Герасимович, Г. В. Яковлев // Энергетика... (Изв. высш. учеб. заведений и энерг. объединений СНГ). -2001. - № 3. - С. 37-47.

10. С а м а р с к и й, А. А. Теория разностных схем / А. А. Самарский. - 3-е изд., испр. -М.: Наука, 1989. - 616 с.

11. Э л е к т р о м а г н и т н ы е процессы в слоистых проводниках и структурах контактных соединений / А. Н. Герасимович [и др.] // Энергетика... (Изв. высш. учеб. заведений и энерг. объединений СНГ). - 2001. - № 2. - С. 11-25.

Представлена кафедрой

электрических станций Поступила 06.06.2009

УДК 621.311

УПРОЩЕННЫЙ РАСЧЕТ ВЕЛИЧИНЫ ТЕХНОЛОГИЧЕСКОГО РАСХОДА ЭЛЕКТРИЧЕСКОЙ ЭНЕРГИИ НА ЕЕ ТРАНСПОРТИРОВКУ

Канд. техн. наук, доц. ПАВЛОВЕЦ В. В.

Белорусский национальный технический университет

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