Обеспечение работоспособности акустических расходомеров вне пределов оптимальных значений расхода жидкости

Полный текст

Введение * Почти все применяемые на практике акустические расходомеры работают в ультразвуковом диапазоне частот и поэтому называются ультразвуковыми. Ультразвуковой расходомер - это устройство, прямым назначением которого является измерение акустических эффектов, возникающих при движении вещества, расход которого необходимо измерить [1]. После сборки счетчика-расходомера производитель проводит ряд стендовых испытаний при различных расходах и давлении рабочего тела в магистрали. Результаты стендовых испытаний заносят в технический паспорт датчика. После этого его можно использовать как в промышленной отрасли и частном секторе, так и в научной области. Однако даже если калибровка расходомера проводилась изготовителем, то с целью повышения точности и достоверности результатов измерения следует дополнительно проводить калибровку счетчика-расходомера непосредственно после окончательного его монтажа в магистрали установки [2]. Некоторые счетчики-расходомеры способны оптимально работать при определенных расчетных режимах расхода рабочего тела, установленных производителем. Поэтому при очень малых или при очень больших расходах присутствует большая погрешность измерений [3; 4]. Значение ЭВМ весьма существенно не только для повышения точности измерения расхода, но и для облегчения эксплуатации приборов и сокращения обслуживающего ее персонала [5]. Целью данного исследования является разработка системы, позволяющей повысить точность измерения и расширить границы оптимальной работы счетчика-расходомера. В статье рассмотрен метод поверки акустического счетчика-расходомера АС-001-15 производства ОАО «Завод “Старорусприбор”». 1. Схема и принцип действия расходомера Счетчик содержит проточную часть (1) и электронный блок (2) (рис. 1). Проточная часть выполнена в виде полого цилиндра особой конструкции, внутри которого установлены пьезодатчики (4), одновременно выполняющие роль излучателей и приемников ультразвуковых волн и системы зеркал (5), направляющих ультразвуковую волну под определенным углом к направлению движения среды. Электронный блок (2) соединен при помощи стоек (3) с проточной частью (1). Рис. 1. Принципиальная схема счетчика-расходомера ФС-0012-15 [Figure 1. Layout diagram of acoustic flow meter FS-0012-15] Работа счетчика основана на измерении разности времени прохождения импульсов ультразвукового сигнала по направлению движения потока рабочего тела и против его потока в трубопроводе [6]. Сигналы от пьезодатчиков поступают в вычислитель счетчика, который пересчитывает эти сигналы в объемный расход. При определении расхода учитывают количество импульсов, поступающих от пьезодатчиков, и их вес (цену). Идея системы исходит из время-импульсного принципа действия акустического счетчика-расходомера. В вычислителе счетчика производится расчет мгновенного объемного расхода по формуле Q = n·K, (1) где Q - расход жидкости; n - количество импульсов, поданных преобразователем в секунду; К - цена (вес) импульса. 2. Постановка проблемы При проведении эксперимента было установлено, что в зависимости от режима движения потока в трубопроводе меняется как количество импульсов, так и цена одного импульса [7]. Количество импульсов в секунду изменяется прямо пропорционально расходу жидкости. Очевидно, что погрешность расчета расхода пропорциональна неточности определения корректирующего коэффициента [8]. Отсюда можно предположить, что возникновение погрешностей при работе на нерасчетных режимах является следствием того, что расходомер запрограммирован на расчет расхода, учитывающий только одно определенное значение цены импульса К. Это объясняется тем, что в вычислителе отсутствуют устройства для сравнения количества импульсов, поданных преобразователем, а также не установлены диапазоны цены импульса, зависящие от режима потока жидкости. Устранить данный недостаток позволяет специальное ПО на базе программы SCADAZetView. Для работы системы счетчик должен иметь выход на компьютер. В данном случае расходомер имеет выход на компьютер через встроенный интерфейс RS-485. Во время работы пьезодатчики счетчика-расходомера подают сигналы на АЦП-модуль, которые регистрируются на компьютере. Полученное компьютером значение является текущим количеством импульсов n, подаваемых счетчиком за единицу времени. Для правильной работы счетчика необходимо также знать значение цены импульса К, зависящей от режима потока жидкости. Данные значения могут быть получены только по результатам проведения поверки. 3. Методика поверки Для проверки данного предположения был выбран метод поверки с помощью массового расхода [9]. С целью минимизации погрешностей при измерении массы и времени были приняты минимальная масса жидкости, равная 4 кг, и интервал времени поверки - 3 минуты, а также усредненные значения количества импульсов n, подаваемых счетчиком за единицу времени. Для определения цены импульса была использована программа, состоящая из следующих компонентов: o счетчика количества импульсов за определенный промежуток времени; o счетчика текущего количества импульсов, подаваемых счетчиком за секунду; o счетчика среднего значения количества импульсов, подаваемых счетчиком за секунду; o секундомера с погрешностью ±0,05 с. Для повышения точности измерения и высокой стабильности в счетчиках малых диаметров - Ду 15, 25, 32, 40 используется так называемое осевое зондирование с практически полным «озвучиванием» всего сечения потока (достигается использованием акустического луча диаметром близким к диаметру измерительного участка счетчика), что позволяет осуществить полное усреднение эпюры скоростей потока воды, и, как следствие, для этих счетчиков предъявляются менее жесткие требования к прямым участкам трубопроводов [6]. Поэтому при проведении эксперимента эпюра распределения скоростей в трубопроводе не учитывается в связи с малым диаметром трубопровода (15 мм). 1. Средства поверки: o персональная ЭВМ, соответствующая минимальным системным требованиям, с собранной схемой на базе программы ZetView [10]; o АЦП модуль; o поверочная емкость объемом до 5 л; o электронные весы с погрешностью ±0,001 кг. 1. Условия поверки: o установившийся режим движения потока в трубопроводе [11; 12]; o изменения расхода в процессе поверки не должно превышать ±5 % от установившегося значения; o отсутствие внешних электрических и магнитных полей, влияющих на работу счетчика. 1. Подготовка к поверке: o подготовка к работе проверяемых счетчиков и измерительного оборудования; o проверка правильности собранной схемы поверки счетчика; o проверка герметичности установки; o после установки и подачи давления счетчик перед поверкой следует выдержать в течении 5 минут [6]. За эталонное значение принят расход, рассчитанный при помощи поверки. 1. Обработка и анализ данных поверки Во время проведения поверки одновременно с наполнением поверочной емкости начинают работу счетчик суммы импульсов и счетчик среднего значения количества импульсов, подаваемых счетчиком за 1 секунду. С помощью секундомера отсчитывается время проведения поверки. По истечении времени подача жидкости в емкость поверки прекращается, вместе с этим одновременно счетчик количества импульсов и счетчик среднего значения количества импульсов, подаваемых счетчиком за 1 секунду, останавливаются. Результаты заносятся в протокол измерений и обрабатываются (таблица). После этого строятся графики зависимости цены импульса от расхода измеряемой жидкости (рис. 2) и зависимости погрешностей измерений счетчика и измерений, сделанных расчетным методом относительно эталонного расхода (рис. 3). Из результатов эксперимента следует подтверждение предположения о зависимости цены импульса от расхода измеряемой жидкости. Результаты эксперимента принято считать удовлетворительными при погрешности ±2 %. В ходе проведения эксперимента максимальное отклонение расчетного расхода от экспериментального не превысило 1 %. Появление двусторонних случайных погрешностей часто связывают с воздействием большого числа факторов на процесс измерений, влияние каждого из которых не является существенным по сравнению с остальными факторами [13]. После того, как зависимость цены импульса от величины расхода была доказана, необходимо было создать метод расчета расхода, опираясь на результаты поверки. Используя полученную зависимость, возможно минимизировать погрешности расхода жидкости. Из графика (рис. 2.) видно, что значение цены импульса зависит от значения расхода, а следовательно, и от количества импульсов в секунду, согласно (1). Также из графика на рис. 2. видно, что при расходе от 134 л/ч и выше значения находятся в пределах требуемой точности измерения (1 %), поэтому было принято упрощение, что на этом участке изменение цены импульса не происходит. Протокол измерений [Table. Test log] Таблица № [No.] mэкс, г [mexpr, gr] T, с [T, s] Nимп [Nimp] nэкс, имп/с [nexpr, imp/s] K, мл/имп [K, ml/imp] Qэкс, л/ч [Qexpr, L/h] Qрас, л/ч [Qcalc, L/h] Qпок.сч, л/ч [Qflmet, L/h] Δр, % [Δcalc, %] Δс, % [Δflmet, %] 1 4016 240 903 3,8 4,447 60,2 60,8 53,0 1,00 12,02 2 5403 180 1302 7,2 4,150 108,1 107,6 99,5 0,46 7,92 3 7978 180 2013 11,1 3,963 159,6 158,4 151,5 0,75 5,05 4 10539 180 2655 14,7 3,969 210,8 210,1 203,0 0,34 3,69 5 13117 180 3277 18,2 4,003 262,3 262,3 253,0 0,03 3,56 6 15511 180 3887 21,7 3,990 310,2 311,7 300,5 0,49 3,13 7 17950 180 4574 25,4 3,924 359,0 358,8 348,5 0,04 2,92 8 20797 180 5271 29,3 3,946 415,9 416,2 403,0 0,06 3,11 9 23207 180 5900 32,8 3,933 464,1 464,5 450,0 0,07 3,05 Рис. 2. Зависимость цены импульса К от расхода измеряемой жидкости Qэкс [Figure 2. Dependence of the pulse value K on the liquid flow Qexpr] Рис. 3. Зависимость погрешностей измерений от расхода жидкости [Figure 3. Dependence of the measuring accuracy on the liquid flow] 2. Определение функциональной зависимости Для создания рабочей программы в конструкторском приложении ZetView необходимо представить график в виде математической зависимости. Поэтому были рассмотрены три варианта описания графика функцией: o множество линейных функций; o полиноминальная функция 6-го порядка; o комбинация линейной и квадратичной функций. Исходя из анализа описанных выше вариантов, наиболее оптимальной была выбрана комбинация линейной и квадратичной функций, так как при одинаковой погрешности всех трех функций данный метод имел наименьшее количество граничных условий, что максимально упрощало расчет в пределах допустимой погрешности. Для нахождения переменных, входящих в функцию, был применен метод решения систем уравнений при помощи обратной матрицы. Для этого, основываясь на данных из таблицы результатов поверки, были составлены матрицы уравнений А, матрицы свободных членов В и матрицы неизвестных членов Х. Перемножая обратную матрицу А на матрицу В, была получена матрица значений, которые подставляются в формулу расчета расхода. В ходе вычислений были получены следующие зависимости: § для изменения цены импульса в диапазоне от 7,2 до 9,4 имп/с: 2 K = 0,0066 · nэкс - 0,1688 · nэкс + 5,0223; (2) чениям и не превышает 1 %, что свидетельствует о работоспособности модели. § для изменения цены импульса в диапазоне от 0 до 7,2 имп/с: К = -0,0875 · nэкс + 4,78; (3) § для изменения цены импульса в диапазоне свыше 9,4 имп/с: К = 3,961 мл/имп. Таким образом, каждому участку графика соответствует своя зависимость. Результат моделирования представлен на рис. 5. Рис. 6. Сравнение погрешностей счетчика, расчета и модели [Figure 6. Comparative analysis of the measuring accuracy of the meter, calculation and mathematical model] Рис. 4. Комбинация линейной и квадратичной функций [Figure 4. Combined line function and quadratic function] Рис. 5. Результат моделирования [Figure 5. Simulation data] Из полученных зависимостей (2) и (3) были составлены условия расчета цены импульса К и расхода Q, после чего был выполнен расчет с применением данной модели. В результате были построены графики зависимости погрешностей измерений счетчика, показаний, сделанных расчетным методом, и показаний, выполненных с помощью модели относительно эталонного расхода для наглядного их сравнения (рис. 6.). Как видно из графика на рис. 6. погрешность расчета при помощи модели соответствует допустимым зна- 3. Проверка работоспособности системы Для подтверждения полученных данных был проведен эксперимент, порядок проведения которого описан в разделе «Методика поверки». Однако в уже имеющуюся программу поверки в среде ZetView были добавлены компоненты, которые позволили вносить необходимую коррекцию при измерении расхода в реальном времени. Пример рабочих элементов приложения для одного из участков диапазона измерения показан на рис. 8. Рис. 7. Сравнение погрешностей счетчика и разработанной системы [Figure 7. Comparative analysis of the measuring accuracy of acoustic meter and created system] По результатам проведенных экспериментов были получены данные, которые позволили определить отклонение значений расхода с учетом примененной рабочей программы от эталонных значений в реальном времени. После этого были построены графики зависимости погрешностей измерений счетчика и измерений, сделанных расчетным методом при помощи модели относительно эталонного расхода. Как видно из графика, представленного на рис. 7, отклонение расчетного расхода от экспериментального не превысило 1 %. Рис. 8. Рабочие элементы участка диапазона измерения [Figure 8. Elements of the measurement range] Заключение В ходе проведения экспериментального исследования были достигнуты следующие результаты: 1. подтверждена поставленная в начале эксперимента теория; 2. проведен сбор статистических данных об отклонениях измерений счетчика-расходомера при различных режимах потока; 3. проведен сбор статистических данных значений цены импульса К с целью формирования диапазонов для работы счетчика-расходомера при различных режимах потока; 4. разработаны алгоритмы коррекции измерения расхода; 5. подтверждена работоспособность системы. Можно с уверенностью заявить, что данная разработка позволяет расширить диапазон оптимальных значений расхода для акустических расходомеров, а также минимизировать погрешности измерения во время работы в любом диапазоне, заявленном производителем в паспорте прибора.

Об авторах

Михаил Владимирович Лапин

Российский университет дружбы народов

Автор, ответственный за переписку.
Email: alo14@yandex.ru

студент 2 курса магистратуры, направление «Энергетическое машиностроение», Инженерная академия

Российская Федерация, 117198, Москва, ул. Миклухо-Маклая, 6

Андрей Олегович Гринин

Российский университет дружбы народов

Email: alo14@yandex.ru

студент 2 курса магистратуры, направление «Энергетическое машиностроение», Инженерная академия

Российская Федерация, 117198, Москва, ул. Миклухо-Маклая, 6

Кирилл Владимирович Шкарин

Российский университет дружбы народов

Email: alo14@yandex.ru

ассистент департамента машиностроения и приборостроения Инженерной академии

Российская Федерация, 117198, Москва, ул. Миклухо-Маклая, 6

Камбиз Сури

Российский университет дружбы народов

Email: alo14@yandex.ru

ассистент департамента машиностроения и приборостроения Инженерной академии

Российская Федерация, 117198, Москва, ул. Миклухо-Маклая, 6

Список литературы