УДК 621.31.048.015
doi: 10.20998/2074-272X.2017.3.07
А.Ю. Чернухин
ВЛИЯНИЕ ПАРАМЕТРОВ КОРОННОГО РАЗРЯДА НА ЭФФЕКТИВНОСТЬ ЭЛЕМЕНТОВ СИСТЕМ МОЛНИЕЗАЩИТЫ
Наведено результаты теоретичных та експериментальних дослгджень процес1в формування об'емного заряду в областг вершини стрижньового блискавкоприймача Розглянуто особливостг формування стримерноТ короны на стрижньових блискавкоприймачах з рнними конфкурацмями вершини в електричному пол грозовоТ хмари. Встановлено залежностг сили струму корони вд напруженостг електричного поля та висоти для кожного вар1анта стрижня Показано, що при цьому часов1 параметри шпульав струму стримеру змтюються не больше нгж на 30%, а заряд стримеру /, як наслгдок, сила струму корони, зб/льшуються з ростом напруженостг електричного поля. Запропоновано метод вим1рювання швидкостг просування стримеру. Метод базуеться на одночасному вижрювант напруги та струму у розрядному пром1-жку. Результати пропонуеться врахувати при сертифжаци блнскавкоприймачгв. Бiбл. 10, табл. 4, рис. 9. Ключовi слова: стрижньовш блискавкоприймач, iмпульсна корона, напружешсть електричного поля грозово!" хмари, струм корони, швидккть стримеру.
Приведены результаты теоретических и экспериментальных исследований процессов формирования объемного заряда в области вершины стержневого молниеприемника. Рассмотрены особенности формирования стримерной короны на стержневых молниеприемниках с различными конфигурациями вершины в электрическом поле грозового облака. Установлены зависимости силы тока короны от напряженности электрического поля и высоты для каждого варианта стержня. Показано, что при этом временные параметры импульсов тока стримера меняются не более чем на 30%, а заряд стримера, и как следствие, сила тока короны увеличиваются по мере роста напряженности электрического поля. Предложено метод измерения скорости продвижения стримера. Метод основан на одновременном измерении напряжения и тока на разрядном промежутке. Результаты предлагается учесть при сертификации молние-приемников. Библ. 10, табл. 4, рис. 9.
Ключевые слова: стержневой молниеприемник, импульсная корона, напряженность электрического поля грозового облака, ток короны, скорость стримера.
Введение. Основой систем молниезащиты зданий и сооружений от прямых ударов молнии являются металлические конструкции в виде стержневых, тросовых и сеточных молниеприемников. Требования к устройству таких систем регламентируются стандартами IEC 62305-1:2010 и IEC 62305-3:2010. В отличие от упомянутых выше систем, которые условно можно назвать «пассивными», в мире на протяжении последнего десятилетия предпринимаются попытки создания активных устройств, обеспечивающих существенное увеличение размеров зоны защиты, по сравнению с зоной защиты классического молниепри-емника Франклина (далее ПСМ - пассивный стержневой молниеприемник). К числу таких устройств относятся, так называемые «Early streamer emission air terminals» (ESE), которые, по утверждению разработчиков, обеспечивают более быстрое по сравнению с ПСМ создание встречного стримера, способствующего перехвату молнии. Декларируемый радиус защиты ESE молниеприемников прямо пропорционален времени опережения [1]. Основой такого подхода является многократно подтвержденный экспериментально факт формирования встречных лидеров с металлических объектов как заземленных так и не имеющих контакта с землей. Однако не все особенности и условия образования встречного лидера достаточно изучены. В частности, вероятность встреч канала нисходящей молнии с восходящими лидерами носит вероятностный характер. Кроме того, скорость продвижения встречного восходящего лидера зависит от потенциала головки канала молнии, который также имеет случайное значение.
На рынке также предлагаются различные варианты устройств, названных рассеивателями
(Dissipater), обеспечивающих снижение вероятности удара молнии в объект за счет создания вокруг вершины объекта объема воздуха, насыщенного заряженными частицами. Этот эффект обеспечивает снижение напряженности электрического поля над вершиной объекта, в результате чего встречный лидер возникает при больших значениях напряженности электрического поля, формируемого каналом молнии, что приводит к снижению вероятности поражения объекта молнией.
В Научно-исследовательском и проектно-конструкторском институте «Молния» Национального технического университета «Харьковский политехнический институт» (НИПКИ «Молния» НТУ «ХПИ») за последние 10 лет проведены исследования более 20 типов различных образцов ESE молниепри-емников и рассеивателей практически всех компаний, представленных на мировом рынке. В 2007 году были проведены сравнительные испытаний ESE терминала компании ERICO (США) в лабораториях Технического университета Валенсии (Испания) и НТУ «ХПИ» (Украина). Результаты испытаний показали, что положения первой редакции (1995 года) стандарта Франции [1] не достаточные для достоверной оценки размеров защитной зоны ESE молниеприемников. В ходе обсуждения результатов было доказано [2], что необходимо внести дополнительные требования на величину дисперсии времени опережения и установить параметры эталонного ПСМ. В новой редакции (2011 г.) стандарта [1] эти рекомендации учтены.
Значительная часть научного сообщества, объединенного в рамках международной конференции ICLP (International Conference Lightning Protection)
© А.Ю. Чернухин
категорически отрицает научную обоснованность положений стандарта [1]. Основными аспектами, не имеющими должного экспериментального подтверждения, являются такие:
- оценка радиуса зоны защиты ESE молниеприем-ника осуществляется путем умножения времени опережения (определяемого как среднеарифметическое значение по 100 разрядам) на скорость продвижения встречного лидера с молниеприемника, которая установлена равной 106 м-с-1. Экспериментальные значения скорости, полученные с помощью скоростных видеокамер, лежат в диапазоне (104 - 106) м-с-1;
- отсутствие документально подтвержденного экспериментами большего радиуса защиты по сравнению с ПСМ. Неоднократные совместные испытания на моделирующих установках не подтверждали этого факта.
Другая часть сообщества, в которую входят в основном производители и дистрибьюторы рассматриваемых устройств, объединена в рамках Международной ассоциации молниезащиты ILPA (International Lightning Protection Association), которая проводит собственные научные симпозиумы.
Возникшее противоречие имеет не только научный аспект, но обуславливает важные проблемы практического применения новых устройств. Ошибочное определение размеров зоны защиты обуславливает снижение вероятности перехвата молнии системой молниезащиты, что может иметь негативные последствия. Тем не менее, на сегодняшний день, стандарт [1] введен в действие в ряде стран, в том числе Франции, Испании, Казахстане, Латвии, которые фактически легализовали использование ESE молниеприемников с учетом расчетной зоны защиты. Настойчивые попытки предпринимаются внедрить стандарт [3] в Украине. Сложившиеся обстоятельства побудили НИПКИ «Молния» НТУ «ХПИ» и функционирующий на его базе Технический комитет Украины по стандартизации в сфере электромагнитной совместимости (ТК 22) инициировать исследования с целью разобраться в физике сопутствующих явлений.
Основой принципа действия обоих типов новых устройств являются процессы коронного разряда с металлических стержней и игл. Поэтому, научное обоснование реальных защитных свойств таких устройств требует детального исследования физики коронного разряда в условиях нахождения изделий в электрическом поле грозового облака.
Анализ последних исследований и публикаций. Исследования процессов коронного разряда проводятся на протяжении многих десятилетий. Однако, они не прекращаются и по сей день. Причина заключается в том, что процессы имеют вероятностный характер, а полнота и достоверность результатов исследований в значительной мере зависит от используемых средств измерения электрофизических величин и скорости видео регистрации. Естественно, что информация, полученная в последние годы, точнее за счет качественно нового уровня упомянутой выше техники.
Принцип действия ESE терминала основан на предположении, что возникновение стримера, переходящего во встречный к каналу нисходящей молнии лидер, с такого терминала происходит ранее, чем с обычного стержневого молниеприемника [1]. Следовательно, величина времени опережения является основной технической характеристикой ESE терминала. При этом, восходящий стример возникает на фоне коронного разряда, являющегося неотъемлемой частью процесса. Поэтому, решение рассматриваемой проблемы в нашем институте начато с исследования закономерностей развития стримерной короны и перехода ее во встречный лидер [2].
В определенной степени появление короны негативно влияет на защитные свойства молниеприемни-ка, т.к. ее наличие препятствует развитию встречного лидера. Вопросам изучения процесса коронного разряда посвящено большое число работ, в том числе [3-6]. Для зажигания короны величина напряженности электрического поля на вершине коронирующего электрода должна превышать определенное критическое значение (Ес). Впервые это значение было получено эмпирическим путем Пиком для электрода цилиндрической формы:
Ec = 29,8S
1 + -
0.3
(1)
(б Я )1/2
где: Ес - значение критической напряженности, кВ/см; д = N/N0; N N - плотность газа при существующих и нормальных условиях соответственно, Я - радиус электрода, см.
В работе [5] представлена формула (2) для электрода сферической геометрии радиусом Я
0.54
Ec = 27,8
1 + -
(2)
(Я)1/2
Сравнение результатов оценки по этим формулам для электродов радиусом Я<1 см, для которых формулы справедливы, представлено в табл. 1.
Таблица 1
Значения критической напряженности электрического поля, вычисленные по формуле Пика и Базеляна [5]
R, см 0,1 0,5 1,0
Ес, кВ/см 58,5 42,6 39,0
Ес [5], кВ/см 75,1 48,9 42,8
Очевидно, что расхождение между значениями критической напряженности возрастает по мере уменьшения радиуса электрода.
В промежутках с большой степенью неоднородности электрического поля плазменная область, возникающая вследствие коронного разряда, может проникнуть в зону с низкой напряженностью только в виде тонкого канала - стримера. Когда стример попадает в область слабого поля, его скорость замедляется, он может остановиться. Такой незавершенный процесс называется импульсной или стримерной короной. Распространение продолжается до тех пор, пока напряженность электрического поля больше
минимально допустимой. Среднее значение напряженности электрического поля Es вдоль стримерного канала при положительной полярности находится в диапазоне 450 кВ-м-1 - 500 кВ-м-1 [3]. В работах Уотерса (Waters, 1987) и Галлимберти (Gallimberti, 1979) показано, что увеличение длины стимера L прямо пропорционально росту напряжения U до тех пор, пока L мало по сравнению с длиной разрядного промежутка. Коэффициент пропорциональности равен, соответственно упомянутым авторам: 0,145 см-кВ-1 и 0,152 см-кВ-1.
По мнению авторов работы [5] тупые молниеотводы более эффективны при перехвате молний, чем острые. Напротив, в работе [7] утверждается, что, меняя радиус молниеотвода в достаточно широких пределах (практически значимых), нельзя повлиять на величину и распределение объемного заряда короны в его окрестности. Это противоречие также требует экспериментальной проверки.
Граничной скоростью нарастания фронта импульса для квазистационарного режима, когда можно пренебречь изменением напряжения за время развития стримеров считается величина 50 кВ-мкс-1. Согласно формулы (2) из работы [7], следует, что для линейного роста напряженности электрического поля, ток короны нарастает линейно, а при стабилизации Е0 ток снижается во времени. Данные экспериментальных исследований [6] этот факт не подтверждают. Если крутизна фронта выше граничной, стример продвигается в течение всего времени роста напряжения на промежутке, при этом скорость его нарастает.
Сильная стримерная вспышка инжектирует в промежуток так много заряда, что поле на всем аноде падает гораздо ниже порога ионизации, в результате чего в развитии разряда наступает пауза [7]. Если напряжение на промежутке не меняется или меняется медленно, время паузы может быть большим - около времени дрейфа ионов на расстояние, сравнимое с радиусом анода (при ra ~ 10 см и средней напряженности электрического поля 5 кВ-см-1 A ~ 10-3 с). Это явление характерно для стержня с острой конической вершиной, потому что элемент, который инжектирует заряд, имеет размер точки. Следовательно, плоская вершина стержня должна иметь преимущества при создании стабильной последовательности стримеров.
Цель статьи - исследование процессов формирования объемного заряда в области вершины стержневого молниеприемника в условиях действия электрического поля грозового облака для оценки эффективности элементов систем молниезащиты.
Результаты теоретических исследований. Очевидно, что молниеприемник, связанный с грунтом проводником имеет нулевой потенциал. При приближении грозового облака на нем наводится электрический заряд, поверхностная плотность которого обусловлена напряженностью электрического поля на вершине в области наибольшей неоднородности. При определенных условиях, описанных выше, возникает коронный разряд и стримеры, в последствии перерастающие во встречный к молнии лидер.
Быстродействие молниеприемника зависит от степени ионизации окружающего воздушного пространства. Если ионизация активная, образуется облако заряженных частиц, осложняющих прорастание встречного лидера. На этом эффекте основаны конструкции рассеивателей, представляющие собой комбинации тонких проводников. Если ионизация слабая, отсутствуют заряды, необходимые для формирования лавин. Значит, можно предположить, что существует конструкция вершины молниеприемни-ка, которая обеспечит ток короны оптимальный для минимизации времени формирования встречного лидера. Этот факт важен для выбора эталонного образца молниеприемника, необходимость определения которого была инициирована в работе [2] и подтверждена стандартом [1].
С целью исследования начальной фазы процесса образования стационарной короны на вершине стержневого молниеприемника, выбрана математическая модель максимально приближенная к реальным условиям при проведении испытаний по стандарту [1]. Электростатическое поле образовано двумя круглыми дискообразными проводящими пластинами, диаметром 10 м каждая. Одной из пластин (верхней) присвоен потенциал 10 кВ, другая - заземлена (и = 0). Расстояние между пластинами модели изменяется в диапазоне от 3 м до 10м. Значение потенциала выбрано таким из соображений уровня напряженности электрического поля Е0 > 1 кВ-м-1, характеризующей предгрозовую обстановку. Высота стержневого молниеприемника принята равной 2 м. Математическая модель задачи такова. В воздушном пространстве электрическое поле является потенциальным, а потенциал удовлетворяет уравнению Лапласа. Граничные условия задачи заданы равенством нулю потенциала на нижней пластине и стержне; положительный потенциал 10 кВ на верхней пластине; остальные границы соответствуют условию непрерывности потенциала.
Цель моделирования - определение такого соотношения параметров стержневого молниеприемника, при котором достигается оптимальное соотношение максимальной напряженности поля около вершины и объема этой области (напряженного объема), благоприятное для стримерного процесса.
Расчет напряженности электрического поля на поверхности заостренного стержня проведен вдоль образующей конуса его вершины. Степень остроты вершины стержня характеризуется коэффициентом заостренности (к), который равен отношению высоты острия к радиусу стержня. Рассмотрены варианты, когда коэффициент принимает значения 1; 2; 3 и 4 для стержня радиусом 0,05 м. Результаты расчета напряженности электрического поля на вершине острия заостренного стержня высотой 2 м, находящего в электрическом поле напряженностью 2 кВ-м-1 представлены в табл. 2. Там же даны расчетные значения коэффициента К усиления напряженности поля в рассматриваемой точке.
Таблица 2
Напряженность электрического поля на вершине острия заостренного стержня высотой 2 м
к 1 2 3 4
Е, кВ-м-1 220 600 1020 1250
К 110 300 510 625
Результаты решения задачи для различных вариантов высоты стержня (в диапазоне от 1 м до 8 м), поперечного сечения (круг, квадрат, многоугольники), его размера (от 10 мм до 30 мм) и формы вершины (плоскость, острие) представлены в работе [8]. Результаты моделирования показали, что максимальное значение напряженности электрического поля на поверхности стержня увеличивается прямо пропорционально высоте стержня. В частности, для стержня квадратного сечения, эта зависимость описывается формулой (3) с погрешностью < 5%.
Ет (к) = 26,7 • Е0 • к , (3)
где Ет(к) - максимальное значение напряженности Е-поля на ребре вершины квадратного стержня с плоской вершиной, В-м-1; Е0 - напряженность электрического поля, в котором размещен стержень, В-м-1; к - высота стержня, м.
С целью проверки результатов, проведено сравнение с результатами расчетов, полученных другими исследователям: по формулам, приведенным в работе [4] для заземленной половины сфероида во внешнем электрическом поле и численного метода работы [9]. Рассчитано напряженность Е-поля на вершине стержня радиусом 0,05 м с вершиной в виде полусферы. Входные параметры модели таковы: Н = 5 м, к = 2 м, и = 104 В. По нашим оценкам максимальное значение напряженности электрического поля на вершине стержня равно 7,2-104 В •м-1, то есть в 36 раз больше чем Е0 = 2-103 В •м- . При использовании графика 10 из работы [4] для сфероида с соотношением высоты к радиусу равным 40, коэффициент усиления определяется близко к 30. В результате численного решения в работе [9] результат близок к 35. Учитывая некоторую разницу между геометрическими фигурами, используемыми в расчетах, и не слишком точную шкалу на графике в работе [4], есть основания считать достоверными результаты применяемого нами численного метода достоверными. Иллюстрация напряженного объема вокруг стержня квадратного сечения представлена на рис. 1.
Рис. 1. Напряженный объем Е-поля вокруг стержня квадратного сечения (модель повернута на 90° влево)
С целью обобщения результатов исследования процесса формирования коронного разряда со стержневых молниеприемников, проведено компьютерное моделирование на основе многофакторного плана эксперимента. В качестве функции отклика выбрано значение напряженного объема пространства вокруг вершины стержня. Напряженный объем - объем
пространства, в котором напряженность внешнего электрического поля превышает 30 кВ-м-1. На основании анализа результатов экспериментальных исследований отобраны количественные и качественные факторы, которые в полной мере определяют значение функции отклика. Эти факторы таковы: напряженность электрического поля, длина стержня, форма поперечного сечения, размеры сечения, форма вершины. Априори, предположено, что уровни значимости этих факторов сопоставимы. Из опыта известно, что функция отклика имеет несколько максимумов, зависящих от формы сечения стержня. Поэтому, выбор нулевой точки модели связан с формой поперечного сечения стержня. Построенные ранее по экспериментальным данным функциональные зависимости показывают, что они являются аналитическими и позволяют задаться конкретным видом зависимости целевой функции от каждого их выбранных факторов. Это позволило ограничиться двух уровневым факторным планом типа 2* (к - число факторов). Для обеспечения возможности проверки адекватности выбранной математической модели и правильности определения констант, использовался ненасыщенный план (число экспериментов превышает число неизвестных констант в модели). Проведен полный факторный эксперимент, на основании которого построена целевая функция отклика.
Величина напряженного объема, определяемая по установленной функциональной зависимости, сопоставлена со значением числа импульсов стримеров, полученных экспериментальным путем. Сопоставление позволило определить необходимый минимальный напряженный объем, при котором начинается процесс образования коронного разряда. Кроме того, установлено значение длины начальных стимеров, в зависимости от значений факторов эксперимента.
Основные результаты теоретических исследований таковы:
- поскольку зона повышенной напряженности локализована около ребер стержня, «напряженный объем» увеличивается с ростом периметра стержня;
- максимальное значение напряженности электрического поля и величины напряженного объема на поверхности молниеприемника, при прочих равных условиях, достигается при использовании стержня квадратного сечения с плоской вершиной. Такой стержень предлагается в качестве нового эталонного для стандарта [1].
Результаты экспериментальных исследований. При некотором значении напряженности электрического поля, в которое помещен молниеприем-ник, на фоне «тихой» короны, ток которой составляет не более сотен микроампер, возникают стримерные вспышки, характеризующиеся током в десятки миллиампер. Только стримерная вспышка при определенных условиях может преобразоваться в лидер. Поэтому, определение критического значения напряженности электрического поля, при котором возникают стримерные вспышки, и характер их поведения при изменении уровня напряженности электрического поля является важной задачей. Моделирование
процессов осуществлено на высоковольтном испытательном стенде ВВС-1.2 НИПКИ «Молния» НТУ «ХПИ». Для оценки быстродействия конкретного молниеприемника предлагается использовать значения параметров стримеров, которые возникают при размещении его в электрическом поле.
Чтобы выявить зависимость между параметрами стримеров и защитными свойствами молниеприемни-ка, проведено комплекс исследований, включая определение характеристик стримеров при постоянном и импульсном напряжении. Исследования на постоянном напряжении позволяют оценить поведение мол-ниеприемника при приближении грозового облака,
другая же часть исследований связана с изучением процесса возникновения стримеров при воздействии импульсного электрического поля, сопровождающего прорастающий канал молнии.
Исследования проведены на стержневых молни-еприемниках, размещавшихся между двумя параллельными металлическими плоскостями. Размеры плоскостей: нижняя - 4x6,5 м, верхняя - 3,6x5,2 м. Схема испытательной установки для исследования характеристик тока короны при постоянном и изменяющемся напряжении и ее внешний вид приведены на рис. 2 и 3 соответственно.
Рис. 2. Схема испытательной установки ВВС-1.2 (С1 = 0,381 мкФ; С2 = 0,385 мкФ; С3 = 0,4 мкФ; Ю1 = 300 МОм; Ю2 = 30 кОм; Ю3 = 60 кОм; Ю4 = 510 кОм; Ю5, Ю6 = 60 кОм; = 75 Ом; Ю = 75 Ом; Тр - трансформатор ИОМ 100/25)
Рис. 3. Внешний вид испытательной установки ВВС-1.2
В данном варианте нижняя плоскость заземлена. Верхняя плоскость потенциальная - на нее подается высокое напряжение отрицательной полярности. Расстояние (S) между вершиной стержня и потенциальным верхним электродом варьировалось в диапазоне от 2,5 м до 0,5 м
При этом, начальное значение напряженности электрического поля в промежутке не изменялось. Результаты экспериментальных исследований подробно представлены в работах [6-8]. Обобщение этих результатов позволяет сделать следующие выводы. Типовая осциллограмма силы тока одиночного стримера представлена на рис. 4. Заряд, содержащийся в стримере примерно 5,85'10-9 Кл. Для сравнения, в работе [5] измерение заряда стримера альтернативным методом дало среднее значение 5'10-10 Кл. Учитывая различие конфигураций вершин стержней и высот стержней совпадение результатов удовлетворительное.
Тек Ли »Acq Complete М Pos: 160.0ns MEASURE
Г '''+ :
;
; 1 ^ ;
■ \ ■ ,,,, ,,,, ,,,, ,,,, 4
:
:
CH1 Fall Time 23.60ns
CH1 Neg Width 35.40ns
CHI Min -¿1.12 V
CH2 None
CH2 Min -3.60V
CH1 1.00VBW CH2 5.00 VEu M 50.0ns
CHI \ -3.61V
-♦-Острие 10мм ^»-Квадрат 12x12мм |
Рис. 5. Зависимость частоты следования стримеров с заостренного (диаметр 10 мм, коэффициент заостренности к = 3) и квадратного (12х 12 мм2) стержней от напряженности электрического поля. Длина каждого стержня 3,4 м
Результаты статистической обработки экспериментальных данных представлены в табл. 3 и 4 для заостренного и квадратного стержней соответственно. Следует отметить, что временные параметры (форма) тока единичного стримера для всех рассмотренных случаев изменяются в пределах ±30% от средне арифметических значений. Существенные изменения наблюдаются только для пикового значения импульса, следовательно, величины заряда стримера. При каждом напряжении было произведено не менее 50 измерений. Представленные в таблицах значения среднеквадратичного отклонения (СКО), однозначно показывают существенно меньший разброс значений для стержня квадратного сечения по сравнению со стержнем с заостренной вершиной. Это свойство характерно для стержня любой длины в рассмотренных пределах. Не случайно, стержни квадратного сечения определены в качестве элементов разрядного промежутка для эталона величины высокого напряжения (стандарт 1ЕС 60062).
Таблица 3
Результаты измерения частоты следования стримеров для заостренного стержня круглого сечения (диаметр 10 мм, коэффициент заостренности к = 3, высота 3,4 м)
U„, Ее, Число стримеров СКО
кВ кВ-м-1 в секунду, N
14,5 2,9 9,16 1,31
20 4 276 7,02
25 5 636 12,58
30 6 78 15,97
50 10 136,6 59,53
75 15 1460 145,60
100 20 2656 322,20
120 24 3940 523,70
Рис. 4. Осциллограмма единичного стримера со стержня квадратного сечения 12x12 мм2 высотой стержня 1,7 м, находящего в электрическом поле напряженностью 30 кВ-м-1 (сигнал снят с резистора номиналом 75 Ом, заряд в импульсе 58,5 10-10 Кл)
Зависимости среднего арифметического числа стримерных вспышек от напряженности электрического поля для заостренного и квадратного стержней представлены на рис. 5.
Из представленных результатов видно, что у стержня с заостренной вершиной есть «мертвая» зона (в диапазоне от 6 кВ-м-1 до 10 В-м-1). Эффект обуславливается точечным характером коронирующего элемента. Наличие такого эффекта объясняет упомянутое выше мнение авторов работы [5], что тупые молниеотводы более эффективны при перехвате молний, чем острые. Важно отметить, что для стержня квадратного сечения такой эффект не наблюдается.
Таблица 4
Результаты измерений частоты следования стримеров для квадратного стержня (12x12 мм2, высота 3,4 м)
U0, кВ Е0, кВ-м 1 Число стримеров в секунду, N СКО
47 9,4 14,7 1,40
60 12 116 4,81
75 15 258 4,74
100 20 1302 20,54
120 24 1720 23,85
150 30 3064 32,77
180 36 4930 93,15
Зависимость величина силы тока коронного разряда от напряженности электрического поля, для описанных выше стержней, представлена рис. 6. Из рис. 6 видно, что при напряженности электрического поля большей 15 кВ-м-1 сила тока с квадратного стержня
больше, чем с заостренного, не смотря на то, что число стримеров с заостренного стержня примерно в 2 раза больше (см. рис. 5). Причина такого факта обуславливается большим зарядом каждого стримера со стержня квадратного сечения.
-КруГПИЙ. 1C1IM. «ТРЫЛ -Ш-КвадрД! 12*12. ппхкнй
Рис. 6. Зависимость силы тока короны от напряженности
электрического поля для заостренного (диаметр 10 мм, коэффициент заостренности к = 3) и квадратного стержней.
Длина каждого стержня 3,4 м
Из представленных экспериментальных результатов следуют такие выводы:
- заостренный стержень (диаметр 10 мм, к = 3, высота 3,4 м) начинает реагировать с напряженности электрического поля 3 кВ-м-1, а стержень квадратного сечения (12x12 мм2) той же высоты не ранее чем с 9 кВ-м-1;
- коронный разряд с заостренного стержня в диапазоне напряженности Е-поля от 6 кВ-м-1 до 10 кВ-м-1 прекращается, что в полной мере соответствует процитированным выше во введении результатам работы [5];
- при напряженности Е-поля, превышающей 10 кВ-м-1 число стримеров с заостренного стержня примерно вдвое больше, чем с квадратного;
- среднее пиковое значение силы тока отдельного стримера выше для квадратного стержня, например, при напряженности Е ~ 10 кВ-м-1это значение составляет 5 мА для заостренного стержня и 60 мА для квадратного стержня;
- значение среднеквадратичного отклонения, при прочих одинаковых условиях для случая квадратного стержня меньше, чем для случая заостренного стержня.
Альтернативный метод измерения скорости стримера. Традиционным методом измерения скорости стримеров-лидеров является сверхскоростная видеосъемка. Для реализации этого метода требуется наличие сверхскоростной видеокамеры, например, типа FASTCAM SAS (с опцией до 106 кадров в секунду) и специального режима освещенности исследуемого разрядного промежутка. Ограничением реализации этого метода является не только высокая стоимость камеры, но и фиксированный сектор обзора, что затрудняет эффективное использование метода при исследовании длинных (более 3 м) разрядных промежутков.
Важным параметром предразрядных процессов есть скорость встречного стримера, формирующегося с заземленного объекта. Чем выше усредненное значение скорости, тем выше вероятность попадания молнии в объект. Известно [10], что скорость стримера зависит от напряженности электрического поля в разрядном промежутке. Так, в области, где напряженность Е поля превышает 3-106 В-м-1, скорость стримера может достигать 107 м/с. Однако, такая напряженность, а, следовательно, и скорость, возможны только в стримерной зоне головки лидера молнии. В работе [10] указано, что измеренная средняя скорость лидера 3,36-105 м/с, минимальная - 8-104 м/с, а максимальная - 2,6-106 м/с. Доказано, что существует минимальная скорость стримера, равная в воздухе при нормальных условиях 103 м/с. При меньших значениях стример останавливается. Следовательно, диапазон изменения исследуемого параметра (скорости стримера-лидера) составляет более четырех порядков.
Известно [10], что волна ионизации формируется с вершины стержня и распространяется в межэлектродном промежутке по направлению к противоположному электроду. Как правило, в воздушных промежутках длиной до 5 м (характерных для большинства лабораторных условий), стример - лидер имеет одноканальную структуру без значительных ветвлений. Вокруг канала стримера формируется слабо проводящий чехол. Принятая в рамках модель стримера базируется на следующих допущениях:
- канал критического стримера одиночный (не ветвящийся);
- диаметр канала с учетом чехла постоянный;
- плотность зарядов в канале однородная.
При этих допущениях, длина стримера прямо пропорциональна заряду в канале (с учетом чехла). Интегральная величина заряда в канале пропорциональна току стримера и времени по формуле (4). Фактически, это площадь под кривой ДО, представленной на осциллограмме (рис. 7)
Q = jI(t)dt,
(4)
где Q - величина заряда в канале; T - время протекания процесса, I(t) - зависимость силы тока от времени с начала процесса.
Тек JL • Acq Complete W Pos: 462.?jus HORIZONTAL
CH1 5.00V CH2 5.00V W 250ди CHI \ -220V
CH2 vertic-al position -0.36 divs (-1.80V)
Рис. 7. Осциллограмма процесса распространения критического стримера
Обозначим длину канала стримера Ьк. Очевидно, что в момент касания противоположного электрода Ьк равна длине воздушного промежутка Ь, а заряд стримера имеет значение 2. Этот момент фиксируется по резкому провалу напряжения и(г), приложенного к воздушному промежутку.
Следует отметить, что момент замыкания промежутка стримером нельзя определять по максимуму тока стримера, поскольку, максимальное значение силы тока достигается после замыкания промежутка и определяется величиной разрядного напряжения в разрядном контуре.
Следовательно, результат эксперимента дает такие исходные данные:
- общий заряд канала (0;
- интервал времени (Т), в течение которого стример пересек воздушный промежуток;
- длина воздушного промежутка (Ь).
По этим данным можно легко определить значение средней на промежутке скорости стримера по очевидной формуле (5).
Ул = ЬТ-1. (5)
Для оценки мгновенных значений скорости стримера в различных точках пространства воздушного промежутка, воспользуемся указанными выше предположениями о характере распространения стримера. Канал стримера, в виде цилиндра длиной Ь, определенного радиуса (значение которого для рассматриваемой задачи не существенно) разбивается на N равновеликих частей. На каждом п-том участке величина заряда (2п) оценивается по формуле (6).
2п = а/N. (6)
Дальнейшие действия иллюстрирует рис. 8, на котором изображена типовая кривая 1(г) в увеличенном масштабе. Определяется площадь под кривой 1(г), дающая величину заряда 2. Эта площадь разбивается на N равновеликих частей. Алгоритм разбиения для некоторых частных случаев будет представлен далее. Отыскиваются соответствующие значения гк. Скорость стримера (Уп) на любом, сколь угодно малом, участке определяется по формуле (7).
Уп = Ь( N -Мп)-1, (7)
где Ып - величина я-ного интервала времени.
с постоянной скоростью. Такой вариант носит гипотетический характер, поскольку в экспериментах не наблюдается.
В качестве первого приближения, кривую зависимости 1(г) можно аппроксимировать прямой, проходящей под определенным углом наклона (к). В этом случае, функциональная зависимость силы тока от времени описывается формулой (8)
1 (г) = к • г. (8)
После подстановки (8) в формулу (4), для общего заряда стримера получим (9)
2 = 0,5к • Т2. (9)
Момент времени (?1) прохождения стримером
первого отрезка (совпадает с интервалом времени) отыскивается по формуле (10)
. (10)
Для вычисления последующих значений гп несложно получить рекуррентную формулу (11)
п =
1п
( Т 2 ^
N
V У
+
(11)
Характер кривой на осциллограмме (рис. 7, нижняя кривая) более точно можно описать экспоненциальной зависимостью. При этом, удобнее рассматривать зависимость в варианте инверсного времени, без учета знака, что позволяет ввести функцию (12):
I (г) = ехр(-а • г), (12)
где а - коэффициент, определяющий скорость затухания функции.
Максимальное значение силы тока положим равным единице, поскольку его абсолютное значение не влияет на дальнейшие результаты. По кривой 1(г) без труда определяется значение параметра а. Для этого достаточно найти момент времени (ге) при котором ток достигает значение 0,368. Параметр а = ге-1.
Рекуррентное соотношение для отыскания значений времени (гп) по предыдущему значению (гп-1), описывается формулой (13).
гп = И ~
ехр(-аЯп-1) - (-1 (1 - ехр(-аТ))
N
. (13)
Рис. 8. Принцип разбиения площади под кривой
Результаты расчетной оценки. Рассмотрим частные случаи функциональной зависимости силы тока стримера от времени распространения. Если сила тока не зависит от времени, то стример распространяется
Это соотношение позволяет легко отыскать все значения гп с любой заданной точностью, определяемой числом разбиений N. Следует еще раз указать, что в этом случае отсчет отрезков времени осуществляется от максимума силы тока к началу.
В представленном на рис. 7 примере, расстояние от вершины молниеприемника до верхнего потенциального электрода составляло 1,2 м. Время развития процесса прорастания стримера, определяемое как интервал от точки расхождения кривых до момента резкого уменьшения напряжения, составляет 9 мкс. Выбрано число интервалов равное 10. Значения соответствующих интервалов времени, рассчитаны с использованием формулы (13). На основании полученных данных, учитывая, что каждый отрезок длины стримера равен 0,12 м, найдена пространственная зависимость скорости, представленная графиком на рис. 9.
В рассмотренном примере скорость стримера изменялась от 1,8-104 м/с до 1,1-106 м/с. Среднее значение скорости стримера на воздушном промежутке составляет 1,3-105 м/с. Рассчитанные значения скорости стримера хорошо согласуются с данными, полученными другими методами [10].
1.Е+07
1.E+Q6
г
О 0,2 0,4 0,6 0,8 1 1,2
iS, М
Рис. 9. График изменения скорости стримера от пройденного пути
Предложенный алгоритм можно описать так:
- определяется интегральное значение заряда Q (площадь под кривой зависимости ДО);
- выбирается число разбиений, из соображений достаточной точности;
- отыскиваются значения интервалов времени, для которых площади под кривой равновеликие и равны Ж-ой части Q;
- скорость на конкретном пространственном промежутке определяется как частное от деления длины промежутка на время его пересечения.
Метод может дополнить процедуру аттестации молниеприемников по стандарту [1].
Выводы. Экспериментально определены функциональные зависимости частоты следования стримеров с металлических стержней длиной 3,4 м круглого сечения с заостренной вершиной и квадратного сечения с плоской вершиной.
Установлено, что частота следования стримеров для случая стержня квадратного сечения имеет детерминированный характер: с ростом напряженности электрического поля частота следования монотонно возрастает с коэффициентом близким к 0,25 с-1-В-1-м.
Значение СКО, при прочих одинаковых условиях для случая квадратного стержня в несколько раз меньше, чем для случая заостренного стержня.
Установлено, что коронный разряд с заостренного стержня в диапазоне напряженности электрического поля от 6 кВ/м до 10 кВ/м прекращается, а частота следования c ростом напряженности достигает 6 кГц.
На основании полученных результатов исследования особенностей процесса формирования стри-мерной короны с вершины стержневого молниепри-емника высотой более 1 м предлагается принять в качестве эталонного молниеприемника для проведения испытаний ESE терминалов в соответствии со стандартом [1] стержень квадратного сечения (12x12 мм2) с плоской вершиной высотой 1 м.
Сертификацию молниеприемников ESE типа рекомендуется начинать со сравнительных испытаний со стандартным молниеприемником, путем размещения их обоих в рабочем объеме испытательного стенда и регистрации частоты разрядов в каждый из них при подаче на промежуток коммутационного импульса напряжения с длительностью фронта более 100 мкс.
Предложен метод оценки скорости продвижения критического стримера, основанный на синхронном измерении напряжения на разрядном промежутке и силы тока стримера с молниеприемника.
На основании совокупности полученных результатов, рекомендуется стандарт [1] не внедрять в качестве национального стандарта Украины до полного введения в требования стандарта научно обоснованных данных.
Результаты получены в рамках выполнения научно-исследовательской работы (ГР №0115U000611), которая выполнялась по заказу Министерства образования и науки Украины в 2015-2016 годах.
Автор благодарит П.Н. Мельникова и С.П. Ша-ламова за помощь при проведении экспериментальных и теоретических исследований, а научного руководителя В.В. Князева - за ценные рекомендации и консультации.
СПИСОК ЛИТЕРАТУРЫ
1. NFC 17-102:2011 (Ed.2) Protection contre la foudre -Systèmes de protection contre la foudre à dispositif d'amorçage. Union Technique de l'Electricité (UTE), 2011. 82 p.
2. Князев В.В., Лесной И.П., Чернухин А.Ю. Результаты исследования параметров активных молниеприемников и рассеи-вателей // Вюник НТУ «ХП1». - 2008. - №21. - С. 78-87.
3. Ortéga P., Heilbronner F., Rühling F., Díaz R., Rodière M. Charge-voltage relationship of the first impulse corona in long air gaps // Journal of Physics D: Applied Physics. - 2005. -vol.38. - no.13. - pp. 2215-2226. doi: 10.1088/00223727/38/13/021.
4. Moore C.B., Rison W., Mathis J., Aulich G. Lightning rod improvement studies // Journal of Applied Meteorology. - 2000. vol.39. - no.5. - pp. 593-609. doi: 10.1175/1520-0450-39.5.593.
5. Базелян Э. М. Эффект объемного заряда короны в мол-ниезащите // Труды IV Российской конференции по молние-защите. НПО «Стример»: Санкт-Петербург, 2014. - С. 1-16.
6. Чернухин А.Ю. Результаты измерения параметров импульсной короны со стержневых молниеириемников // Вю-ник НТУ «ХП1». - 2014. - №50. - С. 155-160.
7. Чернухин А.Ю., Князев В.В. Особенности стримерной короны со стержневых молниеприемников // Вюник НТУ «ХП1». - 2015. - №20. - С. 149-155.
S. Князев В.В., Чернухин А.Ю. Эффекты коронирования в условиях предразрядной обстановки // Труды IV Российской конференции по молниезащите. НПО «Стример»: Санкт-Петербург, 2014. - С. 62-70.
9. Резинкина М.М. Моделирование электрических полей при наличии стержней со скругленными вершинами // Журнал технической физики. - 2015. - Т.85. - №3. - С. 21-27. doi: 10.1134/S1063784215030238.
10. Lightning Protection / Edited by Vernon Cooray // Ins.Eng.Tech., London, UK. 2010. - 1070 р.
REFERENCES
1. NFC 17-102:2011 (Ed.2) Protection contre la foudre -Systèmes de protection contre la foudre à dispositif d'amorçage. Union Technique de l'Electricité (UTE), 2011. 82 p.
2. Kniazev V.V., Lesnoi I.P., Chemukhin A.Iu. Investigation results of the ESE-terminal and dissipaters' parameters. Bulletin of NTU«KhPI», 2008, no.21, pp. 78-87. (Rus).
3. Ortega P., Heilbronner F., Ruhling F., Diaz R., Rodiere M. Charge-voltage relationship of the first impulse corona in long air gaps. Journal of Physics D: Applied Physics, 2005, vol.38, no.13, pp. 2215-2226. doi: 10.1088/0022-3727/38/13/021.
4. Moore C.B., Rison W., Mathis J., Aulich G. Lightning rod improvement studies. Journal of Applied Meteorology, 2000, vol.39, no.5, pp. 593-609. doi: 10.1175/1520-0450-39.5.593.
5. Bazelian E.M. The effect of the space charge of the corona in lightning protection. Trudy IVRossiiskoi konferentsii po mol-niezashchite [Proceedings of the IV Russian Conference on Lightning Protection]. NPO «Strimer», Sankt-Peterburg, 2014, pp. 1-16. (Rus).
6. Chernukhin A.Iu. Parameters of pulse corona on metal roads in strength electric field. Bulletin of NTU «KhPI», 2014, no.50, pp. 155-160. (Rus).
7. Chernukhin A.Iu., Kniazev V.V. A features of streamer crown are from cored air terminals. Bulletin of NTU «KhPI», 2015, no.20, pp. 149-155. (Rus).
8. Kniazev V.V., Chernukhin A.Iu. Effects of corona in pre-discharge conditions. Trudy IV Rossiiskoi konferentsii po mol-niezashchite [Proceedings of the IV Russian Conference on Lightning Protection]. NPO «Strimer», Sankt-Peterburg, 2014, pp.62-70. (Rus).
9. Rezinkina M.M. Simulation of electric fields in the presence of rods with rounded upper ends. Technical Physics, 2015, vol.60, no.3, pp. 337-343. doi: 10.1134/S1063784215030238.
10. Vernon Cooray. Lightning Protection. Ins.Eng.Tech., London, UK, 2010. 1070 р.
Поступила (received) 24.04.2017
Чернухин Александр Юрьевич, м.н.с., НИПКИ «Молния»
Национальный технический университет «Харьковский политехнический институт», 61013, Харьков, ул. Шевченко, 47,
тел/phone +38 057 7076292, e-mail: [email protected]
A.Yu. Chernukhin
Scientific-&-Research Planning-&-Design Institute «Molniya» National Technical University «Kharkiv Polytechnic Institute», 47, Shevchenko Str., Kharkiv, 61013, Ukraine. Influence of coronary discharge parameters on the efficiency of lightning protection system elements.
Purpose. Investigation of the formation of space charge in the region of the apex of lighting rod, under the action of the electric field of a thunderstorm cloud, to evaluate the efficiency of elements of lighting protection systems. Methodology. We have applied the mathematical simulation of electromagnetic field distribution on the top of metal rod by different forms. As a mathematical apparatus, we use the finite element method. We considered two forms of the rod section: round and square. The round (cylindrical) rod has a sharp apex. The square rod has a flat top. The experimental study investigates the features of corona discharge formation. A high-voltage test equipment is created an electric field. Experiments carried out using a configuration consisting a potential plate and vertical rod electrode on grounded plate. The electric field strength varied from 1 kVm to 100 kVm-1. This range corresponds to the thunderstorm condition. Results. We have obtained a correlation between the corona current and the strength of the electric field for various shapes of the rod top. The results of experimental studies confirmed the correctness of the conclusions of theoretical estimates. We show that the time parameters of streamer current pulses vary by no more than 30%, but a streamer charge increase to three time with increasing electric field strength. We proposed and applied a method for measuring the velocity of motion of a streamer in the discharge gap. As a result, it is established that the streamer speed is nonlinear in time. For a discharge gap of 1.2 m, the speed varies from 1.810'4 m/s to 1.110 m/s. Originality. For the first time, we have carried out a complex of studies of corona discharge parameters from lightning protector rod to apply for the certification procedure of ESE terminals. Practical value. Based on the set of obtained results, it is obviously that the standard NF C 17-102:2011 will not be introduced as a national standard of Ukraine before full introduction of scientifically justified data will be include into the requirements of the standard. References 10, tables 4, figures 9. Key words: rod lightning terminal, pulse corona, electric field intensity of thunderstorm cloud, corona current, streamer speed.