APPROACHES TO THE REPRESENTATION OF THE TEMPERATURE DEPENDENCE OF CABLE LINES IN THE CHANNEL ON THEIR LOADING IN THE FORM OF A PARAMETRIC MODEL
- Authors: Trufanova N.M1, Kazakov A.V1, Kukharchuk I.B1
- Affiliations:
- Perm National Research Polytechnic University
- Issue: No 40 (2021)
- Pages: 61-75
- Section: Articles
- URL: https://ered.pstu.ru/index.php/elinf/article/view/2415
- DOI: https://doi.org/10.15593/2224-9397/2021.4.04
- Cite item
Abstract
The cable lines used to power electrical installations differ in designs, operating conditions and placement. During the operation of cable structures, certain difficulties arise in monitoring the thermal condition of structural elements of cable lines. Timely determination of the maximum temperatures inside the cable and areas of possible overheating is necessary to maintain the operability of power supply systems. Therefore, today an urgent task is to develop methods for fast and low-cost determination of the thermal state of cable elements. Object of research: four cable lines located in an underground cable channel with specified ranges of operating current loads. The purpose of the study: to develop a methodology for operational forecasting of technological parameters of cable lines operation for various load conditions. Tasks: development of a mathematical model replacing field experiments; conducting a series of experiments and a set of statistical data; development of a parametric model based on the data obtained. Modeling and experimental part: at the first stage of the study, the real cable channel was replaced by its equivalent mathematical model based on basic conservation equations. The numerical implementation of the model was carried out using the Ansys complex. At the next stage, statistical data were collected using this model, which made it possible to switch to a polynomial parametric model reflecting the dependence of heating of cable line elements on current loads. The model takes into account the mutual influence of lines in the channel and environmental conditions. An assessment of the compliance of the results obtained using the parametric model with experimental data was carried out, which showed the operability of the proposed technique. Significance of the results obtained: the developed two-stage research methodology is applicable to the study of similar objects in the analysis of thermal processes.
Full Text
Эксплуатация современного оборудования подразумевает наличие обширной инфраструктурной сети энергоснабжения. При этом целесообразным становится агрегирование систем питания различного оборудования и отдельных независимых объектов в общие коммутационные линии. Одним из возможных способов организации подобных линий является кабельный канал - совокупность отделенных друг от друга кабельных систем, расположенных на заданном расстоянии друг от друга, зачастую в замкнутом объеме воздушного (или засыпного) пространства [1]. Как правило, в подобных каналах предусмотрена возможность обслуживания и замены отдельных элементов. При передаче электрической энергии ее часть безвозвратно теряется, рассеивается в конечном итоге в виде тепла. Это тепло нагревает как элементы самой кабельной системы, так и окружающие конструкционные элементы и соседние кабельные линии. Максимальная температура кабельной линии обычно ограничена допустимыми температурами эксплуатации ее отдельных элементов. Зачастую наименее стойкой к нагреву элементов конструкции кабеля является его изоляция, представляющая для современных кабелей полимерный материал [2, 3]. При этом в силу ненулевого теплового сопротивления материала изоляции наибольшая температура обычно наблюдается у поверхности токопроводящей жилы, являющейся основным источником нагрева для всего кабеля в целом. В процессе эксплуатации линии измерить данный перегрев затруднительно без нарушения целостности конструкции кабеля, поэтому ограничиваются замером температуры внешней оболочки всего кабеля [4, 5]. Повышение температуры окружающей линию среды ухудшает теплоотвод с кабельной линии, поэтому для недопущения ее перегрева необходимо снижать токовую нагрузку. Поскольку в канале может быть расположено несколько кабельных линий, появляются возможности регулирования внутриканальной температуры путем варьирования токовых нагрузок отдельных потребителей, если такое допускается условиями энергоснабжения [6]. Оценка перегревов внутри конструкции кабельной линии - весьма нетривиальная задача, так как необходимо учесть довольно широкий ряд влияющих на процесс факторов [7-10]. К тому же, как было показано ранее, искомый фактор ограничения нагрузочной способности линии не поддается непосредственному наблюдению в процессе эксплуатации, а может быть оценен только косвенно. Однако при наличии большого количества источников тепла внутри кабельного канала измеряемая температура поверхности кабеля уже не так достоверно описывает состояние его внутренних элементов [11-13]. Для решения вышеописанной задачи может быть предложена методика, заключающаяся в комбинировании двух разноуровневых математических моделей рассматриваемого процесса. Модель первого уровня с достаточной степенью достоверности и детализации описывает конструкцию всех линий, расположенных в кабельном канале, их взаимное расположение, физические и реологические свойства конструкционных материалов и воздуха внутри канала, а также свойства грунта (или иной среды), в котором проложен кабельный канал. Помимо теплофизической составляющей такая модель должна описывать и некоторые электромагнитные явления [14], так как в зависимости от тока линии тепловые потери, вызываемые наводимым электромагнитным полем, будут меняться. Модель первого уровня достаточно громоздка, однако позволяет с необходимой для инженерного расчета степенью точности сымитировать конструкцию реального кабельного канала [15]. Эта модель является базой для проведения серии численных экспериментов, необходимых для детерминации параметров второй более простой модели канала. Имитационная модель второго уровня представляет собой уже абстрагированную от геометрии, свойств, конструкции и прочих параметров реального кабельного канала систему. Это необходимо для возможности быстрой оперативной симуляции такой модели. Подобные параметрические модели часто применяются в системах автоматического и автоматизированного управления разнообразными технологическими процессами и признаны достаточно эффективными [16, 17]. Помимо всего прочего они позволяют быстро решать оптимизационные задачи, разграничивая в режиме реального времени приоритеты загрузки линий, диапазон возможного регулирования нагрузок и многие другие параметры. На основе решений, полученных по этим же моделям, возможна выработка априорных превентивных решений для различных вероятных ситуаций, возникающих на линиях при их реальной эксплуатации [18, 19]. Это позволяет разработчикам таких линий заранее проработать возможные нештатные ситуации и ещё на этапе проектирования устранить или минимизировать их последствия. Существенным недостатком предлагаемого подхода является достаточно жёсткая привязка второй параметрической модели к конкретному кабельному каналу, так как любая модернизация (замена кабеля в линии на другой, добавление новых линий и т.п.) приведет к необходимости построения или доработки модели первого уровня и проведения новой серии экспериментов. Параметры модели второго уровня при этом изменятся, но общий подход к решению всей задачи останется. Для проверки вышеописанной методики был рассмотрен кабельный канал, конструкция которого схематично представлена на рис. 1. Канал содержит четыре кабельные линии, расположенные в полиэтиленовых трубах, размещенных в бетонном монолите. Каждая линия состоит из трех однофазных кабелей с изоляцией из сшитого полиэтилена и алюминиевой жилой с площадью сечения 150 мм2. Рис. 1. Конструкция подземного кабельного канала Температурное состояние всех элементов канала определяется путем совместного решения задач электромагнетизма и тепломассопереноса. Полная постановка задачи описана в [20]. Задача может быть решена различными методами, в том числе методом конечных объемов, реализуемым пакетом Ansys [21]. В результате численной реализации задачи была получена совокупность данных, включающая температуры, скорости и другие физико-реологические параметры. Были определены точки максимального нагрева элементов конструкции кабельных линий. Дополнительно были проведены исследования для определения достаточного количества элементов модели и счетных итераций [22]. Для определения параметров второй модели необходимо спланировать и провести ряд экспериментов для различных токовых нагрузок. Параметрическая модель второго уровня для данного объекта представляет собой совокупность 4 систем уравнений, связывающих максимальную температуру в изоляции линии с собственной нагрузкой и нагрузкой 3 остальных кабельных линий. При этом в силу особенностей модели первого уровня, автоматически будут учтены все дополнительные источники тепла, граничные условия, а также свойства материалов. Каждый набор уравнений дополняется допустимым диапазоном изменения тока в конкретной линии. В параметрическом виде максимальная температура изоляции для первой линии запишется так: , где T0 - начальная температура линии (293,15 К); - температура в самой нагретой точке изоляции первой линии; - приращение температуры, вызываемое протеканием тока по данной линии; - приращения температур, вызываемые влиянием остальных трех линий на первую; I1-I4 - рабочие токи 1-4 линий (параметры модели). Для трех оставшихся линий можно сформулировать аналогичные выражения. Таким образом, полную модель второго уровня, отражающую зависимость максимальных температур изоляции кабельных линий, для рассматриваемой структуры канала можно представить системой параметрических уравнений: Для практического применения данного метода необходимо дополнить систему ограничивающим параметром. Для рассматриваемого процесса таким параметром будет выступать температура на изоляции , для которой необходимо сформулировать ограничения где - предельно допустимая рабочая температура изоляции. Систему в общем виде с учётом ограничений можно представить следующим образом: где i - номер кабельной линии, температура которой определяется; j - номер линии, влияние которой учитывается этим слагаемым; N - общее количество линий. Каждое из представляет собой полином второй степени, сформированный по результатам решения модели первого уровня. Так, для линии 1 полином для был построен по результатам серии численных экспериментов с изменяющимся от нуля до допустимого тока линии с некоторым запасом, необходимым для более точной аппроксимации полиномом. При этом токи остальных линий нулевые. Аналогичным образом были проведены эксперименты и определены приращения температуры первой линии от остальных линий . В табл. 1 указаны экспериментальные значения приращений температуры первой линии, рассчитанные для различных вариантов загрузки кабельного канала. Начальная температура линий принималась равной 293,15 К. Таблица 1 Значения приращений температур линии номер 1, К Номер влияющей линии Рабочий ток, А 50 100 150 200 250 300 350 1 2,49 9,63 21,06 36,71 55,21 77,64 103,76 2 0,76 2,9 7,39 12,88 19,39 27,23 36,04 3 0,78 3,52 7,7 13,43 20,15 28,3 37,88 4 0,65 3,02 6,6 11,52 17,29 24,27 32,49 Аналогичные результаты были получены для остальных линий канала. Для формализации аналитической зависимости температуры от тока был выбран полином вида , где y - приращение температуры, x - ток линии, a и b - коэффициенты, вычисленные по экспериментальным данным. Таким образом, приращение температуры в общем виде будет определяться выражением: Для рассматриваемой задачи расчетные величины коэффициентов, полученные по экспериментальным данным, представлены в табл. 2. При реализации параметрической модели второго уровня были получены числовые значения температур изоляции кабельных линий. Для этой модели время счета и требования к вычислительным ресурсам гораздо меньшие, чем у модели первого уровня. Для верификации модели второго уровня было проведено сравнение T1-T4, рассчитанных для различных токовых нагрузок, с тепловым распределением, полученным из модели первого уровня при тех же условиях. Таблица 2 Значения коэффициентов полинома Приращение температуры Коэффициент полинома а b×103 Линия 1 ∆T11 1,255 0,847 ∆T12 -0,142 0,305 ∆T13 0,433 0,309 ∆T14 0,367 0,265 Линия 2 ∆T21 -0,142 0,305 ∆T22 1,255 0,847 ∆T23 0,367 0,265 ∆T24 0,433 0,309 Линия 3 ∆T31 0,468 0,301 ∆T32 0,051 0,261 ∆T33 1,342 0,880 ∆T34 -0,142 0,305 Линия 4 ∆T41 0,051 0,261 ∆T42 0,468 0,301 ∆T43 -0,142 0,305 ∆T44 1,342 0,880 На рис. 2 представлено тепловое распределение в сечении канала для токовых нагрузок I1 = 50 А, I2 = 100 А, I3 = 150 А, I4 = 200 А, построенное по результатам расчета модели первого уровня. Рис. 2. Тепловое распределение в сечении канала Для тех же токов был проведен расчет с применением модели второго уровня. Результаты расчетов двух моделей с оценкой отклонения приведены в табл. 3. Таблица 3 Результаты сравнения моделей Максимальная температура изоляции Т1 Т2 Т3 Т4 Рассчитанная по первой модели, К 317,280 322,086 330,482 341,309 Рассчитанная по второй модели, К 317,134 321,977 329,853 340,217 Абсолютная погрешность Δ, К 0,146 0,109 0,629 1,092 Максимальное отклонение не превышает 1,1 К, что достаточно для инженерных расчетов. Заключение Предлагаемая двухуровневая реализация задачи не только упрощает расчёт, делая возможным оперативное решение параметрической модели без привлечения больших вычислительных ресурсов, например на микроконтроллере, но и позволяет решать оптимизационные задачи организации пропускной способности конкретного канала с учетом динамически меняющихся условий. Подобный подход и предлагаемые модели могут с успехом применяться не только для подземных кабельных линий, но и для иных задач схожего класса, например, для коллективных трубопроводных шахт, каналов сплит-систем кондиционирования и т.п.About the authors
N. M Trufanova
Perm National Research Polytechnic University
A. V Kazakov
Perm National Research Polytechnic University
I. B Kukharchuk
Perm National Research Polytechnic University
References
- Халитов В. Кабели со СПЭ-изоляцией. Расчет блочной канализации // Новости электротехники. - 2017. - № 5(107)-6(108). - С. 72-78.
- Багаутдинов И.З., Кувшинов Н.Е. Преимущества применения кабелей с изоляцией из сшитого полиэтилена // Инновационная наука. - 2016. - № 3-3. - С. 51-53.
- Щебенюк Л.А., Антонец Т.Ю. Исследования потерь в изоляции высоковольтных силовых кабелей с полимерной изоляцией // Электротехника и электромеханика. - 2016. - № 4. - С. 58-62.
- Удовиченко О.В. Температурный мониторинг кабельных линий высокого напряжения на основе кабелей с изоляцией из сшитого полиэтилена // Линии электропередачи 2008: проектирование, строительство опыт эксплуатации и научно-технический прогресс: материалы III Рос. науч.-практ. конф. с междунар. участ. - Новосибирск, 2008. - С. 301-304.
- Real Time Monitoring of Power Cables by Fibre Optic Technologies. Tests, Applications and Outlook / G.J. Anders, J.-M. Braun, A. Downes John, N. Fujimoto, M-H. Luton, S. Rizzetto // 6th International Conference on Insulated Power Cables (JiCable'03). - Paris, 2003.
- Кухарчук И.Б., Труфанова Н.М. Управление распределением электрической энергии в кабельном канале // Электротехника. - 2019. - № 11. - С. 2-7. doi: 10.3103/S1068371219110099
- Программа моделирования температуры провода и потерь мощности на основе учета режимных и атмосферных факторов / А.Б. Баламетов, Э.Д. Халилов, М.П. Байрамов [и др.] // Программные продукты и системы. - 2018. - Т. 31. - № 2. - С. 396-402. doi: 10.15827/0236-235X.031.2.396-402
- Короткевич М.А., Азаров С.Н. Оценка воздействия кабельных линий электропередачи на окружающую среду // Энергетика. Изв. высш. учеб. заведений и энерг. объединений СНГ. - 2019. - Т. 62. - № 5. - С. 422-432. doi: 10.21122/1029-7448-2019-62-5-422-432
- Зализный Д.И., Широков О.Г. Адаптивная математическая модель тепловых процессов трехжильного силового кабеля // Вестник Гомель. гос. техн. ун-та им. П.О. Сухого. Электротехника и энергетика. - 2014. - № 2. - С. 51-63.
- Baazzim M.S., Al-Saud M.S., El-Kady M.A.Comparison of Finite-Element and IEC Methods for Cable Thermal Analysis under Various Operating Environments // International Journal of Electrical, Computer, Energetic, Electronic and Communication Engineering. - 2014. - Vol. 8. - № 3. - P. 484-489.
- Бирюлин В.И., Куделина Д.В., Горлов А.Н. Расчет температуры нагрева изоляции кабеля с учетом влияния рядом проложенных кабелей // Вестник Казан. гос. энергетич. ун-та. - 2019. - № 2 (42). - С. 56-64.
- Лебедев В.Д., Зайцев Е.С. Расчет температуры жилы однофазного высоковольтного кабеля с изоляцией из сшитого полиэтилена в режиме реального времени // Вестник Иванов. гос. энергетич. ун-та. - 2015. - № 4. - С. 11-16. doi: 10.17588/2072-2672.2015.4.011-016
- Математическое моделирование магнитодинамических процессов при определении потерь в металлических экранах силовых кабелей / А.Г. Щербинин, Е.А. Бородина, Е.В. Субботин [и др.] // Научно-технический вестник Поволжья. - 2020. - № 12. - С. 246-248.
- Навалихина Е.Ю., Труфанова Н.М. Математическое моделирование тепловых и электромагнитных процессов при определении допустимых токовых нагрузок кабельных линий // Известия Томск. политехн. ун-та. - 2014. - Т. 325. - № 4. - С. 82-90.
- Труфанова Н.М., Кухарчук И.Б. Оценка работоспособности кабельного канала на основе численного моделирования процессов термодинамики // Вестник Пермского национального исследовательского политехнического университета. Электротехника, информационные технологии, системы управления. - 2020. - № 35. - С. 30-42. doi: 10.15593/2224-9397/2020.3.02
- Казаков А.В. Численное моделирование охлаждающей системы радиоэлектронного устройства // Вестник Пермского национального исследовательского политехнического университета. Электротехника, информационные технологии, системы управления. - 2020. - № 34. - С. 128-144.
- Terlych A.E., Shcherbinin A.G. Controlling polymer-insulation extrusion process in production of cables // Russian Electrical Engineering. - 2013. - 84 (11). - P. 617-621.
- Мониторинг температуры СПЭ изоляции кабельных линий для определения остаточного ресурса / В.Н. Пугач, Д.А. Поляков, К.И. Никитин [и др.] // Динамика систем, механизмов и машин. - 2017. - Т. 5. - № 3. - С. 84-92. doi: 10.25206/2310-9793-2017-5-3-84-92
- Кучеров Ю.Н. Проблемы обеспечения безопасности потребителей и объектов электроэнергетики при нарушениях работы энергосистемы // Энергетик. - 2007. - № 8. - С. 8-12.
- Труфанова Н.М., Кухарчук И.Б., Феофилова Н.В. Расчет теплового поля кабельного канала с учетом тепловых потерь в экранах кабелей // Вестник Пермского национального исследовательского политехнического университета. Электротехника, информационные технологии, системы управления. - 2018. - № 28. - С. 179-193.
- Коржов А.В. Оценка теплового режима работы изоляции в 2D-модели звена «кабель 6(10) кВ - грунт» в Ansys с учетом подвижек грунта и уставок устройств релейной защиты // Вестник Южно-Урал. гос. ун-та. Сер. Энергетика. - 2013. - Т. 13. - № 1. - С. 39-45.
- Труфанова Н.М., Навалихина Е.Ю. Управление токовой нагрузкой линий в кабельном канале с учетом сложного теплообмена // Известия Самар. науч. центра Рос. акад. наук. - 2012. - Т. 14. - № 4(5). - С. 1318-1321.
Statistics
Views
Abstract - 91
PDF (Russian) - 63
Refbacks
- There are currently no refbacks.