CC BY
48
УДК 621.039.564.5
АЛГОРИТМЫ ОБРАБОТКИ ИНФОРМАЦИИ ПРИ АКУСТИЧЕСКОМ БЕСКОНТАКТНОМ ПОИСКЕ ПРОТЕЧЕК НА ВЕРХНЕМ БЛОКЕ РЕАКТОРА ВВЭР-1000
Канд. техн. наук, доц. БОЛТЕНКОВ В. А., доктора техн. наук, профессора КОРОЛЁВ А. В., MAКСИМОВ M. В., канд. техн. наук, доц. МАСЛОВ О. В.
Одесский национальный политехнический университет
В повышении безопасной и надежной эксплуатации АЭС важную роль играют системы диагностики течей теплоносителя. Особую опасность представляют малые течи, в частности во фланцевых соединениях, которые не фиксируются штатными средствами, но приводят к повреждению основного оборудования АЭС. В этом плане весьма актуально внедрение на АЭС Украины эффективных систем обнаружения и локализации течей на раннем этапе их развития, а именно акустических бесконтактных систем течеискания [1]. Принцип действия таких систем основан на регистрации звукового сигнала, возникающего при истечении теплоносителя через дефект, пространственно разнесенной системой акустических датчиков (микрофонов). Поскольку при истечении теплоносителя, находящегося на рабочем режиме энергоблока при температуре до 320 °С и давлении до 16 МПа, происходят достаточно сложные процессы звукообразования, акустический сигнал от течи является широкополосным и занимает частотную полосу 2-45 кГц при общем уровне звукового давления до 95 дБ. По сравнению с другими типами систем течеискания акустические бесконтактные системы имеют ряд преимуществ: дистанционность, долговечность, живучесть, возможность определения места течи с высокой точностью, хотя их существенный недостаток - невысокая чувствительность по расходу теплоносителя.
Основными функциями системы течеискания являются: обнаружение течи, ее локализация с максимально возможной точностью, оценка расхода течи, прогноз развития течи, классификация типа течи (щелевая/фланцевая). Акустические бесконтактные системы течеискания обеспечивают все перечисленные функции, главная из которых - локализация течи, т. е. оценивание ее координат.
Эксплуатация акустических бесконтактных систем течеискания осложняется рядом факторов. Это, в частности, низкие отношения сигнал/шум (ОСШ), связанные с интенсивными фоновыми акустическими шумами технологических агрегатов и узлов реакторной установки (РУ), наличие многократных переотражений звуковых волн, создаваемых течью в контролируемых помещениях АЭС, повышенная температура в локальных областях истечения теплоносителя, вызывающая особенности распространения звуковых волн. Все перечисленные сложности требуют разработки новых эффективных алгоритмов обработки сигналов от течей. Главные требования к алгоритмам вытекают из того, что течь должна быть зафиксирована в кратчайшие сроки и диагностическая информация должна быть
выдана с максимально возможной точностью, что позволяет принять правильное решение о действиях персонала АЭС при обнаружении течи. Поэтому разработанные алгоритмы исследованы с точки зрения их быстродействия и точности диагностической информации.
В статье приводятся результаты моделирования алгоритмов обработки реальных сигналов от течей, зарегистрированных на экспериментальном теплофизическом стенде НИЛ «Атомспецавтоматика» ОНПУ [2]. Специальный теплофизический стенд позволяет имитировать сигналы истечения теплоносителя от различных видов дефектов в элементах оборудования АЭС на рабочих давлениях и температурах, соответствующих работе энергоблока на мощности. Течи через дефекты в трубопроводах (щельные течи) моделировались на калиброванных щелях, а фланцевые течи - на полномасштабных макетах фланцев верхнего блока (ВБ) РУ ВВЭР-1000. Генерируемый при модельной течи акустический сигнал регистрировался в цифровом виде группой из восьми пространственно разнесенных измерительных конденсаторных микрофонов с улучшенными характеристиками (рабочая амплитудно-частотная характеристика - 20-40000 Гц при неравномерности не более 5 дБ).
Алгоритм локализации (оценивания координат течи) основывается на обработке сигналов, зарегистрированных каждой четверкой микрофонных датчиков. Составим систему гиперболических уравнений:
д/(х0 - Х1)2 + (Уо - Л)2 + (*0 - ->1 (Х0 - Х2)2 + (Уо - У2)2 + (20 - ^ = СТ12;
\1(Хо - X )2 + (Уо - У )2 + (^0 - ^ )2 -(Х0 - Х3)2 + (У0 - Уз)2 + (^0 - гъ)2 = СТ13; (1)
>/(х0 - Х1)2 + (У0 - У1)2 + (^0 - )2 - у!(Х0 - Х4)2 + (У0 - У4)2 + (*0 - ^ = СТ14,
где (х0,у0,z0) - координаты источника звука; (хк,ук,zk), к = 1,..., 4 - то же четырех микрофонов из всей сети акустических сенсоров; т7]-, 7 = 1, ] = 2, 3, 4 - соответствующие оценки разности времен приходов
(РВП) на микрофоны с номерами 7, у; с - скорость звука для условий измерения.
Уравнения (1) описывают гиперболоиды вращения. Решение системы уравнений дает оценку местоположения источника (Х0, у0, ¿0).
Техника оценивания РВП такова. Для каждой пары микрофонов в скользящем временном окне длительностью Т оценивается обобщенная взаимно корреляционная функция (ВКФ) Щ(т). Положение максимума ВКФ позволяет оценить разность времен прихода (РВП) Ту на 7-й и у-й микрофоны
Ту = ащшах Щ (т), (2)
теТ
после чего уравнение (1) решается численными методами.
Типичный вид ВКФ для шумового сигнала от фланцевой течи приведен на рис. 1.
Номер отсчета
Рис. 1. Типичный вид ВКФ для сигнала фланцевой течи: частота дискретизации - 125 кГц
В условиях помещений АЭС наблюдается многократное переотражение звуковых сигналов от стен помещения и находящегося в нем оборудования (реверберация). Это вызывает искажение формы ВКФ и дает смешанную оценку ВКФ, что в свою очередь приводит к ошибке в оценке координат
(хо, уо, 2о)-
Так, моделирование с конфигурацией микрофонной сети, приведенной на рис. 2 (система контроля течей для ВБ), показало, что при двухкратном отражении точность определения координат течи определяется со среднеквадратичной погрешностью 0,7-1,0 м (по всем четверкам датчиков, зарегистрировавших сигналы от течи).
Рис. 2. Конфигурация микрофонной системы для течеискания на верхнем блоке РУ ВВЭР-1000
Для повышения точности был применен метод областей локальной когерентности [3]. Точностью оценивания значений РВП для каждой пары микрофонов т.., которая в свою очередь определяется степенью когерентности сигналов, регистрируемых каждой парой микрофонов. Количественной мерой когерентности двух случайных сигналов х ) и у(?), принимаемых парой микрофонов, является комплексная функция когерентности, равная взаимному спектру мощности Сху (/), деленному на корень квадратный из произведения собственных спектров мощности этих сигналов:
У ху ( / ) =
^ (/)
^ХХ (/(/) '
(3)
где / - частота.
Обычно используется квадрат модуля когерентности (КМК) у 2у (/). По
физическому смыслу как функция когерентности, так и КМК - это коэффициент корреляции на каждой частоте / из анализируемого частотного диапазона. На практике удобнее использовать КМК как функцию, по определению нормированную к единице: в частности, если у2ху(/) = 1, то на всех частотах исследуемого диапазона сигналы считаются полностью когерентными, при у2у (/) = 0 - сигналы некогерентны. На практике при обработке широкополосных сигналов 0 <у2ху (/) < 1, причем порогом когерентности, выше которого когерентность считается существенной для широкополосных сигналов, применяется величина 0,75-0,8 [3]. В каждом скользящем окне строится КМК, и для обработки по (2) используются только те частотные области, в которых КМК выше порогового значения. На рис. 3 показаны две зоны локальной когерентности для пары сигналов, зарегистрированных в условиях сильной реверберации. Применение метода областей локальной когерентности позволило уменьшить погрешность оценивания координат течи до 0,3-0,5 м.
КМК 1,0 0,9 0,8 0,7 0,6 0,5 0,4 0,3 0,2 0,1 0
0,5 1,0 1,5 2,0 2,5 3,0 Частота, Гц х104
Рис. 3. Функция квадрата модуля когерентности с двумя выраженными областями локальной когерентности
Другой проблемой, снижающей точность оценки координат течи, с которой мы столкнулись, является смещение (систематическая ошибка), достигающее в отдельных случаях 15-18 %. Установлено, что смещение оценки координат источника вызвано тем, что в отдельных точках зоны контроля температура воздуха при истечении теплоносителя достигала 175 °С и более. Это вызывало ошибку в оценивании координат путем решения системы (1) из-за существенного изменения скорости звука с с ростом температуры. Температурная зависимость скорости звука описывается соотношением
с = 331,3 + 0^, м/с. (4)
После установки в зоне контроля термодатчиков и изменения реальной температуры воздуха при истечении перегретого теплоносителя рассчитывалось ее усредненное значение по зоне контроля t0. Уравнение (4) со значением t0 = t0 добавлялось четвертым в систему уравнений (1), которая решалась методом Ньютона - Рафсона с учетом скорости звука с для реальных условий. Результаты расчетов показали, что смещение оценки координат протечки (х0,у0,£0) уменьшилось до 1,5-3,0 %.
Разработаны алгоритмы классификации фланцевых и щелевых течей. Проблема такой классификации весьма актуальна в практическом плане, так как возникновение фланцевой и щелевой течей - события различной сложности, требующие разных оперативных действий. Классификация типов течей основана на оценке фрактальных показателей акустических сигналов от течей [4], в частности показателя Херста. Установлено, что показатель Херста для зарегистрированных сигналов от фланцевых течей лежит в пределах 0,62-0,64, а для сигналов от течи через микродефект -в пределах 0,69-0,71. Это позволяет достаточно простыми методами классификации распознавать тип течи. Так, вероятности правильного распознавания шума щелевой течи и шума фланцевой течи составляют соответственно 0,97 и 0,93 при ошибке распознавания второго рода (отнесение объекта к чужому классу) 0,03 и 0,04 соответственно [5]. Кроме того, установлено, что для щелевой течи степень фрактальной динамики процесса звукообразования заметно выше, чем для фланцевой течи. По временному поведению фрактальных показателей можно строить прогноз расхода течи, поскольку чем выше показатель Херста процесса, тем более процессу присуща нарастающая динамика истечения теплоносителя.
В Ы В О Д
Все представленные алгоритмы протестированы на реальных сигналах от щелевых и фланцевых течей, зарегистрированных экспериментально. Алгоритмы построены в расчете на минимальные вычислительные затраты, что позволяет повысить оперативность систем акустического бесконтактного течеискания, и как следствие - обнаруживать течи в оборудовании АЭС на ранней стадии их проявления. Повышение точности диагностической информации облегчает принятие правильного решения при обнаружении течи оператором энергоблока.
Л И Т Е Р А Т У Р А
1. Б о л т е н к о в, В. А. Исследование принципов построения автоматизированной системы мониторинга протечек теплоносителя на верхнем блоке реакторной установки АЭС / В. А. Болтенков, В. И. Верпета, А. Н. Калашников // Вестник Харьковск. гос. политехн. унта, 1999. - Вып. 73. - С. 13-18.
2. К о р о л ё в, А. В. Экспериментальное исследование акустических сигналов при истечении теплоносителя через фланцевые соединения верхнего блока реактора ВВЭР-1000 / А. В. Королёв, А. Н. Литвин, М. В. Максимов // Ядерная и радиационная безопасность. -2001. - № 4. - С. 73-80.
3. Б о л т е н к о в, В. А. Локально-когерентная обработка широкополосных сигналов в системах локализации источника звука / В. А. Болтенков // Збiрник наукових праць Одеського ордена Лешна iнституту сухопутних вшськ. - 2006. - № 12. - С. 21-26.
4. Б о л т е н к о в, В. А. Системы мониторинга течей в оборудовании АЭС: технологии обработки диагностической информации / В. А. Болтенков // 1нтелектуальш системи прий-няття ршень та прикладт аспекти шформацшних технологш: матерiали науково-практично! конференцп. Херсон, 2006. - Т. 1. - С. 45-48.
5. Б о л т е н к о в, В. А. Методика расчета критериев эффективности распознающих систем / В. А. Болтенков, О. Б. Максимова, В. О. Давыдов // Холодильна технжа i технологи. - 2003. - № 1 (81). - С. 83-86.
Представлена кафедрой АЭС Поступила 03.03.2009
УДК 621.1.016: 536.242
О ВЛИЯНИИ ГЕОМЕТРИИ КОРИДОРНОГО ТРУБНОГО ПУЧКА НА ЕГО ТЕПЛООТДАЧУ В ПОПЕРЕЧНОМ ПОТОКЕ*
Канд. техн. наук, доц. СМОЛЯК А. А.
Могилевский государственный университет продовольствия
Теплоотдача пучков труб в поперечном потоке к настоящему времени достаточно хорошо исследована многими авторами. По современным представлениям, при поперечном обтекании коридорного пучка трубки предыдущего ряда закрывают (затеняют) от потока трубки последующего ряда. За обтекаемой трубкой образуется застойная вихревая зона. Движение в этой зоне в значительной мере зависит от относительного продольного шага с2 = S2/d. При уменьшении продольного шага вихревое движение тормозится. Это приводит к более слабому омыванию кормовой части трубки, более сильному затенению последующей трубки и, следовательно, к уменьшению средней теплоотдачи от поверхности трубки в пучке.
Основной поток жидкости в коридорном пучке течет между трубками рядов без изменения направления. Приближенно этот поток можно рассматривать как движущийся в плоском канале. Ширина канала является определяющим размером для характеристики движения и теплообмена. Поэтому и величина поперечного шага также должна влиять на теплоотдачу коридорного пучка труб. Но такое влияние значительно слабее. По этой причине в уравнениях подобия для коридорных пучков чаще всего предлагается учитывать только влияние продольного относительного шага.
* Публикуется в порядке обсуждения.
