Journal influence
Bookmark
Next issue
A program of modeling wire temperature and power losses based on operation and atmospheric factors
Abstract:In modern conditions, due to the growth of power consumption and noticeable climatic changes, the interest in calculating admissible mechanical and current loads on the aerial wires has been increased. The maximum use of overhead line capacity can be achieved with reliable information on the line state. The accuracy of calculating electric power losses in overhead power transmission lines might be increased by determining active resistance of wires taking into account the operating current flowing along the lines, the ambient air temperature, wind speed and solar radiation heat. In these conditions, it is very important to have software for evaluating and continuously monitoring wire temperature under various weather conditions, to reliably determine the permissible current load and to be able, if unacceptable line overload occurs, including repair modes, to carry out measures for unloading the remaining work of an aerial line. Modern programs for calculating power systems modes do not include thermal balance equations of aerial lines. This leads to the impossibility of taking into account their temperature conditions due to the mismatch of wire active resistance to their temperature. The authors have developed an algorithm and a program for calculating the specific active resistance of aerial lines wires taking into account the air temperature, working current, wind speed and solar radiation. They also carried out a quantitative assessment of the load current influence, ambient temperature, solar radiation and wind speed on the active resistance of aerial wires. There are the results of calculation of limited current loads for aerial wires. Based on the calculations to evaluate the effect of load current, ambient temperature, solar radiation and wind speed on active resistance of aerial wires, it is established that when calculating power losses for wires without considering the resistance temperature dependence, relative errors can reach 26 % or more, which is not permissible. Therefore, it is required to automatically detect the wire temperature in order to improve the accuracy of the calculation of active power and energy losses.
Аннотация:В современных условиях в связи с ростом электропотребления и заметными климатическими изменениями повысился интерес к расчетам допустимых механических и токовых нагрузок на провода воздушных линий. Максимальное использование их пропускной способности может быть достигнуто при наличии достоверной информации о состоянии линии. Повысить точность расчета потерь электроэнергии в воздушных линиях электропередачи можно, определив активное сопротивление проводов с учетом протекающего по линиям рабочего тока, температуры окружающего воздуха, скорости ветра и теплоты солнечного излучения. В этих условиях очень важно иметь ПО для оценки и непрерывного контроля температуры провода при различных погодных условиях, достоверно определять допустимую токовую нагрузку и иметь возможность при возникновении недопустимой перегрузки линии, в том числе в ремонтных режимах, выполнять мероприятия по разгрузке остающихся в работе воздушных линий. Современные программы расчета режимов энергосистем не включают в себя уравнения тепловых балансов проводов воздушных линий. Это приводит к невозможности учета их температурного режима из-за несоответствия активных сопротивлений проводов их температуре. Авторами данной статьи разработаны алгоритм и программа расчета удельного актив¬ного сопротивления проводов воздушных линий с учетом температуры воздуха, рабочего тока, скорости ветра и солнечной радиации. Проведена количественная оценка влияния тока нагрузки, температуры окружающей среды, солнечной радиации и скорости ветра на активное сопротивление проводов воздушных линий. Приведены результаты расчета предельных токовых нагрузок для проводов с активным сопротивлением. На основе проведенных расчетов по оценке влияния тока нагрузки, температуры окружающей среды, солнечной радиации и скорости ветра на активное сопротивление проводов воздушных линий установлено, что при расчете потерь электроэнергии для проводов без учета температурной зависимости сопротивления относительные погрешности могут достигать 26 % и более, что недопустимо. Поэтому требуется автоматическая регистрация температуры провода для повышения точности расчета потерь активной мощности и энергии.
Authors: Balametov A.B. (balametov.azniie@gmail.com) - Azerbaijan Scientific-Research and Design-Prospecting Power Engineering Institute, Azerbaijan State Oil and Industry University (Professor, Head of a Scientific Direction), Baku, Ph.D, Khalilov E.D. (elmanxalilov2010@mail.ru) - Azerbaijan Scientific-Research and Design-Prospecting Power Engineering Institute (Head of Department), Baku, Ph.D, M.P. Bayramov (balametov.azniie@gmail.com) - Azerbaijan State University of Oil and Industry (Senior Lecturer), Baku, Ph.D, K.A. Agakhanova (balametov.azniie@gmail.com) - Sumgait State University (Senior Laboratory Assistant), Sumgait, | |
Keywords: aerial line, heat balance equation, operating current, weather conditions, wire temperature, numerical simulation, the automated information system |
|
Page views: 10371 |
PDF version article Full issue in PDF (19.46Mb) |
В связи с увеличением энергопотребления нагрузки на воздушные линии (ВЛ) электропередачи растут довольно высокими темпами, а сетевое строительство пока явно отстает. Максимальное использование пропускной способности ВЛ может быть достигнуто при наличии достоверной информации об их состоянии, а также актуальных данных о температуре проводов и плотности протекающего по ним тока. Одним из направлений повышения точности расчета переменных потерь электроэнергии в ВЛ является определение активных сопротивлений проводов с учетом протекающего по линиям рабочего тока, температуры окружающего воздуха, скорости ветра и теплоты солнечного излучения. В соответствии с [1–6], при определении технологических потерь электроэнергии по электросетям необходимо учитывать температуру провода, влияющую на величину активного сопротивления. Температура провода, в свою очередь, зависит от режимных и климатических факторов: плотности тока, температуры окружающего воздуха tв, силы и направления ветра, интенсивности солнечной радиации и т.д. [1–14]. Удельное активное сопротивление проводов ВЛ
где R20 – удельное активное сопротивление при температуре провода 20 °С, Ом/км; α = 0.004 – температурный коэффициент электрического сопротивления сталеалюминиевых проводов, 1/град; tnp. – температура провода, °С. Температура провода в формуле (1), как правило, неизвестна. Она зависит от большого числа факторов и прежде всего от протекающего по проводу электрического тока, температуры окружающего воздуха и скорости ветра. В традиционных расчетах переменных потерь электроэнергии активные сопротивления принимаются по справочникам, в которых они приведены к температуре 20 °С. В лучшем случае активные сопротивления пересчитываются с учетом температуры воздуха tв. Цель авторов данной статьи – разработка алгоритма и оценка влияния тока нагрузки, темпе- ратуры провода, солнечной радиации на актив- ное сопротивление проводов при скорости ветра 0.5÷10 м/с, а также определение диапазонов погрешностей расчета годовых переменных потерь электроэнергии. Зависимость температуры провода от параметров внешней среды может быть представлена в виде уравнения теплового баланса провода, коэффициенты которого для различных случаев впервые были получены Бургсдорфом [1]. Уравнение теплового баланса для установившегося теплового режима записывается следующим образом:
где I – ток линии, A; R20 – удельное активное сопротивление при температуре провода 20 °С, Ом/м; tв – температура воздуха, °С; aк, aл – коэффициенты теплоотдачи провода при конвективном и лучистом теплообмене, Вт/(м2∙°С); Wc – теплота солнечного излучения, поглощаемая 1 м провода в единицу времени, Вт; dпр – диаметр провода, м. В уравнении (2) левая часть представляет собой поглощаемую теплоту солнечного излучения и теплоту от токовой нагрузки, а правая часть состоит из суммы конвективных потерь и потерь на радиацию. Коэффициент теплоотдачи лучеиспусканием определяется по формуле Стефана–Больцмана [5, 6, 9]:
где ε – степень черноты поверхности провода для окисленного алюминия, равная 0.13 отн. ед.; С0 – коэффициент излучения абсолютно черного тела, равный 5.67∙10-8 Вт/м2; S – площадь поверхности провода, м2. Потери мощности при теплоотдаче конвекцией [2] определяем по формуле
где aк – коэффициент теплоотдачи конвекцией, Вт/(м2∙°С); tрад. – температура нагрева солнечной радиацией, °С. Коэффициент теплоотдачи конвекцией [3] определяем по формуле
где kv – коэффициент, учитывающий влияние угла направления ветра к оси ВЛ, равный 0.5; v – ско- рость ветра, м/с; а – коэффициент температуропроводности воздуха, равный 18.8∙10-6 м2/с; λв – теплопроводность воздуха, равная 0.0244 Вт/(м∙°С). Из уравнений (2–5) для тока [5, 6] можно получить
Получение зависимости температуры провода от тока нагрузки, температуры воздуха и скорости ветра в явном виде из уравнения (6) не удается. В связи с этим для решения уравнения (6) относительно температуры провода и получения аналитических зависимостей от тока нагрузки, температуры воздуха и скорости ветра применяется следующий алгоритм. 1. Получение зависимости тока от температуры провода, скорости ветра и температуры воздуха для заданного провода по формуле (6) для конкретных значений tв, tрад., ν в табличном виде. 2. Аппроксимация полученных зависимостей в областях температур от –40 °С до +40 °С полиномом третьей степени:
В связи с зависимостью коэффициента теплоотдачи от температуры провода температуру неизолированных проводов можно вычислить по формуле (6) и методом последовательных приближений на основе уравнения
Здесь к – номер итерации. По формуле (8) можно вычислить температуру провода при известном токе и различных погодных условиях. Для определения теплоты солнечного излучения, поглощаемого проводом, существует несколько моделей [3–7]. Наиболее полно эти вопросы рассмотрены в [5–8]. Теплота солнечного излучения определяется уравнением
где eп = 0,6 – коэффициент поглощения провода; kН – коэффициент, учитывающий влияние высоты над уровнем моря; Wc – интенсивность суммарной радиации (прямой и отраженной), Вт/м2; yc – активный угол наклона солнечных лучей. Из уравнений теплового баланса (6) провода следует, что минимальное допустимое значение температуры tдоп. состоит из трех составляющих: tдоп. = tI + tв + tрад., (10) где tI – составляющая от нагрева провода током нагрузки за счет потерь энергии в активном сопротивлении провода. Температура tрад. зависит от интенсивности солнечной радиации, высоты и плотности облаков. Интенсивность радиации солнца изменяется в те- чение года. По его состоянию в сентябре от 14 ч 30 мин. до 15 ч 30 мин. при синем небе и отсутствии ветра для проводов АС в [3] получено уравнение
Здесь Крад. = 92.0375 °С/м. В период с 20 до 24 часов и с 24 до 7 часов принимается tрад. = 0. Для проводов разных марок для учета солнечной радиации в июне для времени с 11 до 13 ч (наибольший нагрев) при расчете в формулу (11) вводятся коэффициенты Кt = 1.15, Кt' = 1 при синем небе и Кt' = 0,5 при сером небе и видимом солнце. Определение допустимой температуры провода. В настоящее время длительно допустимая температура провода для нормального режима принята равной 70 °С. Для утяжеленного режима температура провода выше 70 °С, но ниже максимально допустимой по условиям механической прочности 90 °С и допустимой по условию сохранения габаритов. Аварийный режим по токовой перегрузке – установившаяся температура провода превышает 100 °С или допустимую по условию сохранения габаритов или гололедная нагрузка превышает допустимую. Максимально допустимый ток провода [6] определяется из выражения
С ростом скорости ветра V (и ак) значение Iдоп. увеличивается, однако для повышения надежности допустимый ток вычисляют для наименьшей нормативной скорости ветра (0.5 м/с) и его продольном направлении (kV = 0,5). Для моделирования изменения температуры провода во времени уравнение теплового баланса записывается в дифференциальной форме и расчет ведется методами численного интегрирования. Время нарастания температуры провода. Температурная характеристика переходного процесса может быть рассчитана в соответствии с методом, описанным в [8, 9]. Для нагрева провода от температуры t1 с током I1 до t2 с током I1 требуется время, определяемое минутами. В программе реализован следующий алгоритм нарастания температуры провода в зависимости от времени. Постоянная времени рассчитывается по формуле Характер нарастания температуры провода до установившейся температуры определяется переходным процессом по формуле
В Азербайджанском научно-исследовательском и проектно-изыскательском институте энергетики разработан алгоритм, реализованный в виде про- граммы моделирования температуры провода с учетом тока нагрузки ВЛ и влияния атмосферных факторов. Уравнения (1–14) составляют основу комплекса алгоритмов и программ данного института. Блок-схема программы и экранные формы приведены на рисунках 1 и 2. Программа использовалась при моделировании технических потерь электроэнергии на Азербайджанской ЭЭС и мониторинга ВЛ 110 кВ для определения допустимой токовой нагрузки по условию передачи мощности электрической станции Шимал-1 по пяти линиям 110 кВ и введения ограничений по мощности и нашла практическое применение. С использованием программы проведены рас- четы для провода АС 240/39 с удельным активным сопротивлением R20 = 0.121 Ом/км и диаметром провода d = 21.6 мм. Получены зависимости тока от температуры провода, скорости ветра и температуры воздуха по формуле (7) при ν = 0,5 м/с и tрад. = 0 °С в области температур от –40 до +40 °С. После подстановки численных значений R20, d и ν в формулу (7) получаем значения токовых нагрузок провода при заданной температуре воздуха и температуре провода, которые представлены в таблице 1. Аппроксимация кубической зависимости тока от температуры провода, скорости ветра и температуры воздуха имеет достаточную для практики точность. В таблице 2 представлены коэффициенты кубической зависимости из формулы (7) для расчета температуры провода АС 150/24 при ν = 0.5 м/с и tрад. = 0 °С, полученные по разработанной программе. Представление графика кусочно-нелинейной аппроксимации делением на две части tв ≤ 0 и tв > 0 и применение для аппроксимации полинома второй степени имеют определенные удобства и практически приемлемую точность по сравнению с результатами полиномов более чем во второй степени. К примеру, tв = –40÷0 °С, Iпр.(-40÷0) = –0.164x2 + + 2.3118x + 490.58 и tв = 0÷90 °С, Iпр.(0÷90) = = –0.0149x2 + 4.1767x + 498.45. Сравнение с полиномом 3-й степени показывает, что СКО имеет значение 4.35 %, тогда как кусочно-нелинейная аппроксимация полинома вто- рой степени делением на две части имеет лучшую точность: первая часть имеет СКО 3 %, а вторая – 0.21 %. По формуле (7) для коэффициентов таблицы 2 рассчитана температура провода в зависимости от тока в диапазоне 110–350 А для температуры воздуха от –40 °С до +40 °С. График зависимости допустимого тока провода АС 240/39 при температуре воздуха –40 °С, 0 °С и 40 °С и скорости ветра ν = 0.5 м/с приведен на рисунке 3.
Для получения зависимостей активного сопротивления от тока нагрузки, температуры воздуха и скорости ветра подставим температуру провода из таблицы 2 в формулу (1). На рисунке 4 представлен график зависимос- ти удельного активного сопротивления провода АС 240/39 от тока при температуре воздуха –40 °С, скорости ветра 0.5 м/сек. и tрад. = 0 °С. Зависимость температуры провода от скорости ветра и тока провода АС150 при tв = 20 °С и tрад. = = 0 °С показана на рисунке 5. Погрешности расчета сопротивлений линий электропередачи могут быть достаточно большими и составлять от –8 до 24 % и более. В таблице 3 представлены результаты расчетов суточных переменных потерь электроэнергии по формуле (10) с учетом диапазона токов для проводов АС 150/24. Таблица 3 Переменные потери электроэнергии провода АС-150/24 Table 3 Variable power losses of AC-150/24 wire
Рассчитана переходная температура провода АС 240/39 по параметрам погоды: v = 0.5 м/сек., tв = 40 °С, Wc = 900 Вт/м. При этом происходит изменение тока от начального Iнач. = 70 A до конечного Iкон. = 100 A. Постоянная времени рассчитана по формуле (13): Результаты В современных условиях эксплуатации электрических сетей важное значение имеют ПО для мониторинга температуры провода с учетом погод- ных условий, возможность повышения пропускной способности линии и выполнение мероприятий по разгрузке для недопущения перегрузки линий. Авторами разработаны алгоритм и программа расчета удельного активного сопротивления приводов BJI с учетом температуры воздуха, рабочего тока, скорости ветра и солнечной радиации. Проведена количественная оценка влияния тока нагрузки, температуры окружающей среды, солнечной радиации и скорости ветра на активное сопротивление проводов ВЛ. Алгоритм и программа протестированы на примере разных марок проводов АС 70/11, АС 95/16, АС 120/19 и АС 150/24, АС 240/39. Приведены результаты расчета зависимости удельного активного сопротивления от температуры воздуха, рабочего тока, солнечной радиации, предельных токовых нагрузок, характера нарастания температуры провода до установившейся температуры, опреде- ляющейся переходным процессом. Установлено, что при расчете потерь электроэнергии для проводов без учета температурной зависимости сопротивления относительные погрешности могут достигать 26 % и более, что недопустимо. Поэтому в условиях эксплуатации необходима автоматическая регистрация температуры провода для повышения точности расчета потерь активной мощности и энергии. Литература 1. Бургсдорф В.В., Никитина Л.Г. Определение допустимых токов нагрузки воздушных линий электропередачи по току их проводов // Электричество. 1989. № 11. С. 1–8. 2. Зарудский Г.К., Зиннер Л.Э., Сыромятников С.Ю. Расчет температуры проводов воздушных линий электропередачи СВН на основе метода критериального планирования эксперимента // Вестн. МЭИ. 1997. № 12. C. 85–90. 3. Никифоров Е.П. Предельно допустимые токовые нагрузки на провода действующих ВЛ с учетом нагрева проводов солнечной радиацией // Электрические станции. 2006. № 7. С. 56–59. 4. Воротницкий В.Э., Туркина О.В. Оценка погрешностей расчета переменных потерь электроэнергии в ВЛ из-за неучета метеоусловий // Электрические станции. 2008. № 10. С. 42–49. 5. Левченко И.И., Сацук Е.И. Нагрузочная способность воздушных линий электропередачи в экстремальных погодных условиях // Электричество. 2008. № 4. С. 2–8. 6. Фигурнов Е.П., Жарков Ю.И., Петрова Т.Е. Уточненная методика расчета нагрева проводов воздушных линий электропередачи // Электрические станции. 2013. № 9. С. 54–59. 7. Guide for Thermal Rating Calculations of Overhead Lines. WG B2.43, Tech. Rep. 601, CIGRE, 2014, 95 p. 8. Механошин Б.И., Шкапцов В.А., Васильев Ю.А. Повышение эффективности использования существующих ВЛ на основе анализа их технического состояния и данных мониторинга температуры проводов // Электро. 2007. № 6. С. 37–41. 9. Герасименко А.А., Тимофеев Г.С., Тихонович А.В. Учет схемно-режимных и атмосферных факторов при расчете технологических потерь электроэнергии в распределительных сетях // Журн. Сибирского федер. ун-та. Сер.: Техника и технологии. 2008. Т. 1. № 2. С. 188–206. 10. Петрова Е.В., Бубенчиков А.А., Кириченко Н.В., Птицына Е.В. Разработка комплекса программ сравнения методов расчета потерь электрической энергии в воздушных линиях электроэнергетических систем с учетом режимных и климатических факторов // Омский науч. вестн. 2011. № 3. С. 221–225. 11. Баламетов А.Б., Байрамов М.П. Моделирование температуры провода для расчета потерь электроэнергии воздушных линий // Проблемы энергетики. 2013. № 2. С. 4–12. 12. Cecchi V., Leger A.St., Karen M., Nwankpa C.O. Incorporating temperature variations into transmission-line models. IEEE Transactions on Power Delivery, 2011, vol. 26, i. 4, pp. 2189–2196. 13. Mo Y., Zhou X., Wang Ya., Liang L. Study on operating status of overhead transmission lines based on wind speed variation. Progress In Electromagnetics Research M, 2017, vol. 60, pp. 111–120. 14. Rahman M., Kiesau M., and Cecchi V. Investigating the impacts of conductor temperature on power handling capabilities of transmission lines using a multi-segment line model. Proc. SoutheastCon, 2017, pp. 1–7. |
Permanent link: http://swsys.ru/index.php?page=article&id=4473&lang=en |
Print version Full issue in PDF (19.46Mb) |
The article was published in issue no. № 2, 2018 [ pp. 396-402 ] |
Perhaps, you might be interested in the following articles of similar topics:
- Разработка прецедентного модуля для идентификации сигналов при акустико-эмиссионном мониторинге сложных технических объектов
- Возможности мониторинга динамики развития проекта в интеллектуальном проектном репозитарии
- Разработка графической оболочки для параллельных расчетов на базе платформы OpenFOAM
- Технологические процессы с гибкими связями типа «звезда» для интерактивных программных систем
- Метод оценки технического уровня изделия по мониторингу рекламаций в системе управления качеством производства
Back to the list of articles