Сетевое издание
Современные проблемы науки и образования
ISSN 2070-7428
"Перечень" ВАК
ИФ РИНЦ = 1,006

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

Петрова Е.В. 1 Гиршин С.С. 1 Ляшков А.А. 1 Бигун А.Я. 1
1 Федеральное государственное бюджетное образовательное учреждение высшего профессионального образования «Омский государственный технический университет»
Рассмотрен новый подход к расчету уравнения теплового баланса воздушных линий электропередачи при условии вынужденной конвекции, основанный на приближенном аналитическом решении уравнения нагрева. Уравнение теплового баланса после преобразований представлено как алгебраическое уравнение четвертой степени относительно абсолютной температуры поверхности провода. В предлагаемом подходе решение полученного алгебраического уравнения четвертой степени осуществляется при помощи метода Фиррари. При решении данного уравнения методом Фиррари получается четыре корня уравнений, но только один соответствует реальному тепловому режиму. Сравнительные расчеты показали, что значения температуры, полученные предлагаемым методом и путем численного решения исходного уравнения, совпадают друг от друга. Предложенный метод может использоваться при расчете теплового режима линии в условиях вынужденной конвекции и изменяющейся нагрузки, а также для расчета потерь энергии.
уравнение теплового баланса
потери энергии
изолированный провод
температура провода
вынужденная конвекция
1. Воротницкий В.Э., Туркина О.В. Оценка погрешностей расчета переменных потерь электроэнергии в ВЛ из-за неучета метеоусловий / В.Э. Воротницкий, О.В. Туркина // Электрические станции. – 2008. – №10. – С. 42–49.
2. Герасименко А.А. Комплексный учет режимно-атмосферных факторов в расчете активного сопротивления и потерь электроэнергии в ЛЭП / А.А. Герасименко, И.В. Шульгин, Г.С. Тимофеев // Оптимизация режимов работы электрических систем: Межвуз. сб-к научных трудов. – Красноярск, 2008. – С. 188–206.
3. Гиршин С.С. Упрощение уравнений теплового баланса воздушных линий электропередачи в задачах расчета потерь энергии / С. С. Гиршин [и др.] // Омский научный вестник. – 2013. – № 1. – С. 148–151.
4. Зарудский Г.К. Уточнение выражений для расчета температуры проводов воздушных линий электропередачи сверхвысокого напряжения / Г.К. Зарудский, С.Ю. Сыромятников // Вестник МЭИ. – 2008. – № 2. – С. 37–42.
5. Курс высшей алгебры. / А.Г. Курош. – М.: Наука, 1968. – 431 с.
6. Левченко И.И. Нагрузочная способность и мониторинг воздушных линий электропередачи в экстремальных погодных условиях / И.И. Левченко, Е.И. Сацук // Электричество. – 2008. – № 4. – С. 2–8.
7. Математическая модель расчета потерь мощности в изолированных проводах с учетом температуры / С. С. Гиршин [и др.] // Омский научный вестник. – 2009. – № 3. – С. 176–179.
8. Оценка дополнительных потерь мощности от снижения качества электрической энергии в элементах систем электроснабжения / Долингер С.Ю. [и др.]. //Омский научный вестник. – 2013. – № 2. – С. 178–183.
9. Схематические решения активной фильтрации кривой тока в четырехпроводной трехфазной сети для обеспечения качества электрической энергии / В. Н. Горюнов [и др.] // Омский научный вестник. – 2011. – № 3. – C. 214–217.
10. Уточнение метода расчета температуры провода при постоянной нагрузке с учетом климатических факторов / С. С. Гиршин [и др.]; Омский гос. техн. ун-т. – Омск, 2010. – 23 с. – Деп. в ВИНИТИ 08.04.2010, N198-В2010.

Практически важными задачами в электроэнергетике являются улучшение качества электрической энергии и анализ потерь в линиях электропередачи [7-10]. Решение второй задачи предполагает использование уравнений теплового баланса. Уравнения теплового баланса проводов в общем случае нелинейны (в разных источниках приводятся несколько различающиеся варианты этих уравнений, неодинаковость которых обусловлена принимаемыми допущениями) и поэтому решаются численными [2, 6] или приближенными аналитическими [1, 3, 4] методами. Те и другие имеют известные недостатки: при использовании численных методов затрудняется анализ результатов, а приближенные аналитические методы дают достаточную точность лишь в ограниченном диапазоне исходных данных.

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

1. Расчетные формулы

Уравнение теплового баланса провода воздушной линии может быть записано в виде [3]:

(1)

где ΔP0 – потери активной мощности на единицу длины, вычисленные при сопротивлении, приведенном к 0 ºC; α – температурный коэффициент сопротивления; Θвнеш и Θокр – температуры поверхности провода и окружающей среды в ºC; dпр – диаметр провода; Sиз – тепловое сопротивление изоляции на единицу длины провода; αвын – коэффициент теплоотдачи вынужденной конвекцией; εп – коэффициент черноты поверхности провода для инфракрасного излучения; C0 = 5,67·10-8 Вт/(м2·К4) – постоянная излучения абсолютно черного тела; Tвнеш и Tокр – абсолютные температуры поверхности провода и окружающей среды; As – поглощательная способность поверхности провода для солнечного излучения; qсолн – плотность потока солнечной радиации на провод.

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

(2)

Здесь введены следующие обозначения:

, (3)

. (4)

Видно, что выражение (2) представляет собой алгебраическое уравнение четвертой степени относительно температуры Tвнеш. Найдем его решение методом Феррари [5].

Слагаемое, содержащее третью степень неизвестного, отсутствует. Тогда первым шагом решения уравнения методом Феррари является преобразование с помощью вспомогательного параметра β [5]:

(5)

Параметр β выбирается таким образом, чтобы выражение в квадратных скобках было полным квадратом. Это соответствует условию

(6)

Решая кубическое уравнение (6) по формуле Кардано [5], получим следующее выражение для β:

. (7)

При таком значении β уравнение (5) сводится к двум квадратным уравнениям:

(8)

Корни уравнений (8):

(9)

(10)

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

1. Предположим, что температурная зависимость сопротивления отсутствует (α=0). Из формул (3), (4) при таком допущении следует, что A1>0, A0<0. Численные расчеты показывают, что такие знаки сохраняются и в реальных условиях при α≠0, по крайней мере при токах, соответствующих нормальным режимам эксплуатации линий.

2. Если A1>0, A0<0, то из формулы (7) видно, что β – действительное положительное число.

3. Так как A1>0 и β>0, то формула (10) дает комплексно-сопряженные корни, а меньший корень, определяемый по формуле (9), представляет собой действительное отрицательное число. Из этого следует, что реальному тепловому режиму линии соответствует больший из корней, определяемых по выражению (9). Таким образом, окончательная расчетная формула для температуры внешней поверхности провода имеет вид

(11)

Температура токоведущей части определяется по следующему выражению [10]:

(12)

Приведенные формулы справедливы как для изолированных, так и для неизолированных проводов воздушных линий. В последнем случае следует принимать Sиз=0.

2. Пример расчета

Пусть по проводу марки SAX-50 протекает ток 200 А. Исходные данные для расчета (табл. 1) примем те же, что были использованы в [3].

Таблица 1

Исходные данные для расчета температуры провода SAX-50

Наименование и обозначение параметра

Численное значение

Погонное активное сопротивление при 0 ºC r0

0,000663 Ом/м

Тепловое сопротивление изоляции на единицу длины Sиз

0,193566 м∙К/Вт

Температурный коэффициент сопротивления a

0,0043 ºC-1

Диаметр провода dпр

0,0127 м

Коэффициент теплоотдачи вынужденной конвекцией αвын

13,3764 Вт/(м2∙K)

Степень черноты поверхности провода εп

0,8

Поглощательная способность поверхности провода для солнечного излучения As

0,9

Температура окружающей среды Θокр

0 ºC

Плотность потока солнечной радиации qсолн

526,291 Вт/м2

Ток в проводе I

200 А

Расчет. Потери активной мощности на единицу длины при сопротивлении, приведенном к 0 ºC

Вт/м.

Коэффициенты A1, A0 по формулам (3), (4):

Определяем параметр β по формуле (7):

Температура внешней поверхности провода по формуле (11):

°C.

Температура токоведущей части по формуле (12):

°C.

Полученные значения совпадают с результатами численного решения уравнения теплового баланса, приведенными в [3].

Из приведенного расчета видно, что сделанные предположения о знаках коэффициентов для данного примера оказались верными. Однако следует заметить, что это будет справедливо не при любых исходных данных. Так, видно, что при некотором достаточно большом значении потерь мощности (тока) знаменатель во втором слагаемом формулы (3) и в первом слагаемом формулы (4) обращается в ноль. Поэтому коэффициенты A1, A0 при увеличении тока проходят через точки разрыва и меняют знаки. Однако это не свидетельствует о том, что данный метод несправедлив при больших токах, а отражает известный факт, что при больших токах установившегося теплового режима провода не существует. Это является следствием температурной зависимости сопротивления, которая проявляет себя как положительная обратная связь, препятствуя установления теплового равновесия.

Рецензенты:

Черемисин В.Т., д.т.н., профессор, директор Научно-исследовательского института энергосбережения на железнодорожном транспорте, заведующий кафедрой ФГБОУ ВПО «Омский государственный университет путей сообщения», г. Омск;

Кузнецов А.А., д.т.н., профессор, Федерального государственного бюджетного образовательного учреждения высшего профессионального образования «Омский государственный университет путей сообщения», г. Омск.


Библиографическая ссылка

Петрова Е.В., Гиршин С.С., Ляшков А.А., Бигун А.Я. АНАЛИТИЧЕСКОЕ РЕШЕНИЕ УРАВНЕНИЯ ТЕПЛОВОГО БАЛАНСА ПРОВОДА ВОЗДУШНОЙ ЛИНИИ В УСЛОВИЯХ ВЫНУЖДЕННОЙ КОНВЕКЦИИ // Современные проблемы науки и образования. – 2015. – № 1-1. ;
URL: https://science-education.ru/ru/article/view?id=18506 (дата обращения: 18.04.2024).

Предлагаем вашему вниманию журналы, издающиеся в издательстве «Академия Естествознания»
(Высокий импакт-фактор РИНЦ, тематика журналов охватывает все научные направления)

«Фундаментальные исследования» список ВАК ИФ РИНЦ = 1,674