CC BY
37
УДК 681.78
МАТЕМАТИЧЕСКАЯ МОДЕЛЬ И ФУНКЦИЯ ПРЕОБРАЗОВАНИЯ ЛАЗЕРНО-ТЕЛЕВИЗИОННОГО МЕТОДА КОНТРОЛЯ УРОВНЯ ЖИДКОСТИ
© 2014 Б.В. Скворцов, С.А. Борминский, Д.И. Блинов, Д.М. Живоносновская
Самарский государственный аэрокосмический университет имени академика С.П. Королёва (национальный исследовательский университет)
Поступила в редакцию 21.11.2014
В настоящей статье предлагается теоретическое обоснование и разработка математической модели лазерно-телевизионного метода измерения уровня по меткам, создаваемым на контролируемой поверхности одним или группой лазерных лучей с последующей обработкой изображения. Ключевые слова: математическая модель, функция преобразования, лазерно-телевизионный метод, контроль уровня жидкости, лазерное излучение.
В настоящее время во всём мире наблюдается рост количества измерительных устройств, в которых применяется лазерное излучение. Измерения, проводимые с помощью лазерного излучения, содержат в себе ряд преимуществ. Использование лазерного излучения позволяет выполнять измерения с высокой чувствительностью, поскольку вследствие высокой энергетической яркости лазерного луча улучшаются энергетические соотношения в приёмных устройствах. Также, важным фактором в пользу применения лазерного излучения является то, что лазерный луч представляет собой направленное электромагнитное излучение, а скорость распространения электромагнитной волны мало зависит от окружающих факторов, таких как плотность среды, температура, влажность, давление. Известен ряд устройств [1, 2, 3, 4] для измерения уровня жидкости с помощью лазерного луча (одного или нескольких).
Выделим основные параметры измерительного лазерного луча, зависящие от контролируемого уровня жидкости:
-время распространения сигнала от поверхности до приёмника.
-амплитуда и фаза пришедшего на фотоприёмник сигнала.
-спектральный состав сигнала. -позиция отражённого сигнала на изображении, регистрируемом фоточувствительной матрицей приёмника (определение данного па-
Скворцов Борис Владимирович, доктор технических наук, профессор кафедры «Электротехника», научный руководитель НИЛ «Аналитические приборы и системы». E-mail: [email protected].
Борминский Сергей Анатольевич, кандидат технических наук, доцент кафедры «Электротехника». E-mail: [email protected].
Блинов Дмитрий Игоревич, аспирант кафедры «Электротехника». E-mail: [email protected]. Живоносновская Дарья Михайловна, студентка. E-mail: [email protected]
раметра с последующей обработкой используется в лазерно-телевизионном методе контроля уровня жидкости).
В данной статье предложены несколько вариантов реализации лазерно-телевизионного метода контроля уровня жидкости, а именно:
- по одной лазерной метке;
- по двум лазерным меткам, расстояние между которыми заведомо известно;
- по N лазерным меткам, расстояние между которыми заведомо известно.
Для пояснения рассматриваемого метода рассмотрим лазерно-телевизионный уровнемер, представленный на рис. 1. Данный уровнемер содержит: телевизионную камеру 1, цифровое вычислительное устройство 2, источники лазерного излучения 3, 4, 5, 6, блок управления 7, блок индикации 8. По команде блока управления активируются источники лазерного излучения, на контролируемой поверхности появляются метки лазерных лучей, а фотоприёмная матрица регистрирует изображение контролируемой поверхности, которое поступает в цифровое вычислительное устройство, где по формулам определяется искомое расстояние Н.
Сущность измерительной процедуры в самом общем случае иллюстрируется рис. 1, 2 и состоит в том, что на поверхности жидкости, уровень которой контролируется, с помощью лазерных лучей создаются метки А1, А2,... А .. Ап.
Каждая из этих меток имеет координаты в пространстве А.(Х. Известно, что все эти
точки принадлежат одной плоскости, уравнение которой можно определить через три произвольные точки Аг , Аг , Аг+1 по формуле [5].
X - X,-X - X,
X, . -X, .
г+1 г-1
Y - Y-
Y - Y
Y - Y
-1 г+1 * г-1
Z - Zt -Z - Z,
Z . - Z, .
г+1 г-1
= 0. (1)
Рис. 1. Структурная схема лазерно-телевизионного уровнемера
Выражение (1) преобразуется к виду
ОХ + ОУ + Тг + Р = 0 , (2)
где
О = у - у-1)(+1 - г --1) - (у+1 - у-1)(г , - г .-1),
Т = (X, - X,-1 )(у+1 - у-1) - (X,+1 - X,-1 )(7. -у-1),
д = -((X, - x¡_l)(z¡+l - г-) - (xм -- гн)),
(3)
р = -Gx.. - оу.. - тг,.
1-1 ^ 1-1 1-1
Предположим, что расстояние между всеми точками известно и определяется по формуле
Ь ^^/(Xn-X)г+(У™-У)г+(Zn-Z;)2=.
= SX+bYГ+h¿I. . (4)
\ ш ш т к 7
Контролируемая поверхность снимается те-ле(фото)камерой.
Выбранные точки имеют координаты на изображении
4 -1 (м-; п-1X 4 (М; пX 4 ++1(М+1; п+1).
Из законов геометрической оптики, по правилам сохранения пропорций, при съёмке со стороны X (центр линзы находится на оси X) между координатами точек в пространстве X, У,2 и на фотоматрице р., П справедливы следующие соотношения [6] (при г > / )
X/ у/
п (5)
Рис. 2. Иллюстрация к математической модели лазерно-телевизионного метода
контроля уровня жидкости: 1 - фотоматрица, 2 - оптическая система, 3 - контролируемая поверхность
хности до главной плоскости оптической системы, которое связано с контролируемым уровнем Н и другими конструктивными параметрами формулой (см. рис. 1)
Г=1 + ;- И = Г ... - Н , (6) где гт = Ь + В - высота установки центра оптической системы над нулевым уровнем жидкости, то есть при Н = 0, г = гт.
Тогда выражения, связывающие координаты точек поверхности и их изображения на фотоматрице будут иметь вид
М =■
X/
п =■
у/
". (7)
гт-н-/ гт-н-/
Рассмотрим варианты реализации лазерно-телевизионного метода контроля уровня жидкости.
1. Контроль по одной метке. Это возможно, если метка создаётся лазерным лучом, строго параллельным главной оптической оси телекамеры и перпендикулярным контролируемой поверхности. При этом фокусное расстояние известно. В этом случае при смещении уровня вдоль оси X координаты метки X. и У. не изменяются. Тогда уровень и координаты метки будут связаны соотношениями
М =-
п =-
г-/ г-/'
где f -внешнее фокусное расстояние объектива, для стандартных объективов это величина известная, либо определяется экспериментально. Здесь X - расстояние от контролируемой повер-
М =
X/
г -н-/
т
п =
у/
г - н - /
т
(8)
X
н=гт - / ■ (1+М М
Н = 1т-/■ (1 + У). (10)
п
Из формул (8), (9) видно, что измерения можно проводить по любой координате метки. Для проведения процедуры измерений нужно при известных значениях Хг, Уг, /, задаться начальным значением координаты метки
Мо ="
X/
По =-
у/
(11)
2т - Н0 -/ 2т " Н0 " /'
При изменении уровня координаты метки А (м, ) на изображении будут смещаться по фотоматрице (рис. 3).
п I
Птк Пш
Пок По!
® №ок Утк /7
Рис. 4. Перемещение меток на фотоматрице при изменении уровня
Уровень и координаты меток связаны соотношениями (8). Тогда уровень можно найти из системы уравнений
М1 =
Мп =
X/
2 - Н - /
т
X/
2т - Н - /
П =
У/
Пп =
2 - Н - /
т
У /
п
2т - Н - /
(13)
V пп = (Хп - X,)2 + (Уп - У,)2
Рис. 3. Смещение координаты метки на фотоматрице в зависимости от уровня
При этом каждому значению уровня будет соответствовать точка на фотоматрице. Единственным условием применения метода является захват телекамерой метки на контролируемой поверхности при минимальном и максимальном значениях уровня, текущее значение которого определяется по формулам (9), (10). Очевидно, что измерения будут более точными, если мы сделаем отсчёт по отдельным координатам, но результат усредним
Н= + (12)
где НМ,Нп - уровень, определённый по формулам (9), (10) через соответствующие координаты. Контроль уровня по одной метке требует безусловного знания конструкционных параметров - фокусного расстояния / и высоты что не всегда возможно.
2. Контроль уровня по двум меткам, расстояние между которыми известно. Это возможно, если метки создаются лазерными лучами, строго параллельным главной оптической оси телекамеры и ортогональными контролируемой поверхности. В этом случае также при смещении уровня вдоль оси X координаты X ,Уг меток не изменяются, X. - X = 0. В этом случае при изменении
7 г п ^
уровня будет изменяться расстояние между метками, (рис. 4).
где м,п,,мп;пп-текущие координаты на изображенииг-ой и п-ой меток соответственно. Из(13) получаем
(X - X) / (У - У) /
Мп-м = \ п '" , п-п =- п ,Л/ '
2т - Н - /
2т- Н - /
X - X , =
п
(м -мру* - Н - /) / '
У - У = (Пп -П )(2т - Н - /)
п
/
I ,„ =
(2 т - Н - /) /
^(Мп -М)2 + (Пп -П)2, (14)
Отсюда контролируемый уровень определяется из выражения
Н = 2т - / [1 + ■ ь;п ] = гя - / (1+), (15)
\1(Мп -М) + (Пп -п) Ь'п
где ь\п = V(мп - м,)г + (пп - п,)2 - расстояние между метками на фотоматрице.
Выражение (15) является функцией преобразования прибора, позволяющей определять уровень по координатам метки.
3. Контроль по N меткам, расстояние между которыми известно. Метки создаются лазерными лучами, параллельным главной оптической оси телекамеры и ортогональными контролируемой поверхности. Для трёх меток г, к, п система (13) будет иметь вид
М-1 = М =
Xl-l/
у-1/
гт-н-/,п' 'гт-н-/ X/ п у/
гт - н - / гт - н - /
у+1/
М = )См/ п = ^ гя -н-Гп гт-н -/
2 , /V
(16)
Ь,-1, , = (X,-1 - X у + (у-1 - у г у Ьl+l, = (X,+1 - X ,)2 + (у+1 - у )2 Ь2 ,-1,,+1 = (X1+1 - X,-1)2 + (у+1 - у,-1)2.
Здесь значения Ь .., Ь . ., Ь . .- считаются
^ г-1, г ' г, г+1 ' г-1, г+1
известными. Координаты меток на изображении являются данными измерительной процедуры. Целью решения является определение уровня Н. Система (16) путём последовательных подстановок сводится к виду
Ь ,-1,, =
Ь.
г - н - /) /
(гт - н - /)
л/(М-1 -М,)2 + (п,-1 -п)2
/
7(М+1 - М,)2 + (п+1 -п) (17)
Ь 1-1,1+1
(гт - н - /) /
л/(М+1 -М)2 + (п+1 -п)2
Система (17) является вырожденной, так как размеры отрезков линейно связаны одним и тем же уравнением (9). По каждому из уравнений можно найти уровень Н, но результат усреднить
н = гт - / [1+
Ь, +
>/(М+1 -М)2 + (п+1 -п):
1 ^
н = — У н.
N £ ^
-I=гт -/(1+7^) .(18)
(19)
Усреднённый контроль уровня по нескольким
меткам повышает точность измерений, в том числе компенсирует оптические искажения, если метки разбросаны по противоположным границам оптической системы.
Выражения (15)-(18) пригодны при горизонтальном положении контролируемого уровня, когда он ортогонален оптической оси. Реально оптическая ось может быть не ортогональна контролируемой поверхности, что имеет место быть при колебаниях резервуара, при этом уровень остаётся горизонтальным. Рассмотрим метод, позволяющий исключить влияние на результат измерений колебаний уровня относительно горизонта.
Выберем за ноль системы координат центральную точку объектива и оптическая ось объектива чтобы совпадала с осью 2. При выборе за начало координат других точек отсчета математическое описание процедуры формирования изображения сильно усложнится. Метки Аг на контролируемой поверхности формируются параллельными лазерными лучами, ортогональными к контролируемой поверхности (рис. 1). Расстояние между метками известно только при горизонтальной поверхности. При колебаниях уровня относительно оптической оси телекамеры оно изменяется. Из рис. 1 видно, что координаты Хг,¥гметок Аг(Х;У;2г) не зависят от уровня Н. Если за начало координат принять центр оптической системы, то их можно считать известными, определяемыми конструкционными размерами прибора. От уровня зависят только координаты 2г. Причём расстояние от каждой метки до поверхности будет разным и определяться через конструкционный параметр 2т по формуле
2 =2- Н, (20)
г т г
где Н-уровень г-ойметки относительно горизонтального дна.
Возьмём четыре точки А,Б,С,Б, координаты ХА, УЛ, ХБ, УБ, ХС, УС, Хв, Ув которых будем считать известными (рис. 5).
а) б)
Рис. 5. Расположение лазерных меток: а - на контролируемой поверхности, б - на фотоприёмной матрице телевизионной камеры
Метки принадлежат одной плоскости, поэтому удовлетворяют уравнениям [5]
Рх- ХА + Ру-УА + Рг-2А = 1 Р - Х„ + Р-У0 + Р-20 = 1
Рх- Хс + Ру-Ус + Р2-2С = 1 Рх' Хв + Ру-У0 + Рг-20 = 1.
(21)
Здесь коэффициенты Рх, Р, Р, 2А, 2В, 2С, 2В, определяющие расположение плоскости в пространстве неизвестны. Однако известно [5], что при такой форме записи уравнения плоскости
Р2 =
1
(22)
= + + 2с + 2В ) .
(24)
Обратимся к проекции меток на матрицу (рис. 5, б). Точки проекции будут иметь координаты 4ма;Г|А ), в(мв;пв), п), пъ) • Учитывая (5) для координат проекций каждой точки можно записать
Ма =
Мб =
Мс =
2Л - Г
Хв/
2В -/
Хс/
2с - /
Хъ/
Па =
Па =
Па =
Па ='
?А /
2а - / , 2в - /, 2с - /,
Уъ/
2ъ-Г 'А 2В-/. (25)
Оптическая система может состоять не из одной линзы, однако главное её назначение - пропорциональная передача изображения. Поэтому соотношения (24) универсальны и применимы для любых сколь угодно сложных оптических систем.
Параметры ма , Па , Мв, Пв, М С 11с ' Г-Б' 1 \ъ являются данными измерительной процедуры, определяются по изображению и считаются известными.
Из совокупности соотношений (21) - (25) мож-
но набрать нужное количество уравнений для определения неизвестных параметров, среди которых нас интересуют 20 , 2А , 2В , 2С, 2В , которые связаны с контролируемым уровнем по формуле (20), (23). Для контроля уровня по четырём точкам всего нам необходимо 9 уравнений. Исключив путём взаимных подстановок из уравнения (21) координату 2, получим систему из 4-х уравнений.
где 2Ц - отрезок, отсекаемый плоскостью контролируемой поверхности по оси 2 от начала координат, он и определяет расстояние от центра оптической системы до контролируемой поверхности и связан с искомым уровнем по формуле Н = 2 - 2 . (23)
т ц 4 7
Учитывая, что оптические метки симметричны относительно начала координат, можно записать также
Рх-Ха + Ру-Уа + р
Р-Х. + Р-У„ + Р
Рх^с + Ру'Ус + Р;
Рх-Хъ+Ру-Ув + Р;
(Ма + ха ) /
Ма
(Мв + Хв ) /
Мв
(Мс + Хс ) /
Мс
(Мъ ,+ Хъ ) /
= 1
= 1
= 1
= 1
(26)
Мъ
Неизвестными и искомыми в этой системе являются Рх, Ру, Р;, / .
Система нелинейная. Несмотря на то, система имеет четыре неизвестных, она не имеет решения, так как является вырожденной относительно фокусного расстояния. Для решения проблемы введём обозначение
Я = Рг/ =
/
(27)
' Ц
Тогда система (25) сводится к трём уравнениям
Рх 'ХА + ру-Уа + Я-(1 + ^) = 1
Ма
Х
(28)
Рх-Хв + Ру-Ув + Я-(1 + ) = 1 Мв
Рх-Хс + РуУс + Я-(1 + —) = 1 Мс
Это линейная система, в которой нас интересует Q, имеет прямое решение по формулам Крамера.
Я =
А
Я
А
(29)
где
А =
Х„ У. 1 +
1 + Ха
1 + Ма Х в Ха УА 1
Мв , А Я = Хв Уа 1
1 + Хс Мс Хс Ус 1
. (30)
Дальше при известном фокусном расстоянии из (26) получаем
Z ц = f
ц Q •
(31)
Отметим, что при известном фокусном расстоянии величина 20 может быть найдена из соотношения (24), которое с учётом (25),можно переписать в виде
X. X„ XГ Xr
■). (32)
г ц = 0,25/ (4
МА МВ Мс Ми
Искомый уровень нефтепродукта определяется по формуле (23).
Выражения (23), (29) - (32) определяют уровень колеблющейся относительно горизонта поверхности по координатам полученных на фотоматрице оптических меток и в совокупности определяют функцию преобразования прибора. Результаты, полученные по формулам (31),(32) можно усреднить. Выражение (31) предпочтительнее, так как оно не накладывает ограничений на симметричность лазерных меток относительно центра оптической системы и позволяет контролировать уровень по трём произвольным меткам, созданными на поверхности лазерными лучами, параллельными главной оптической оси.
Таким образом, общий алгоритм решения задачи определения уровня жидкости по четырём оптическим меткам в случаях горизонтальных колебаний уровня состоит в следующем:
1. Для конкретных значений координат меток ХА,УА,ХВ,УВ, ХС,УС решаем систему уравнений (28), из которой находим параметр &
2. По выражениям (31) - (32) находим искомый уровень.
3. При необходимости результат усредняем
Описанные в статье оптические методы позволяют повысить точность измерения уровня жидкости, обеспечить контроль за наклоном резервуара, что дает возможность применять данные разработки на железнодорожном, автомобильном транспорте, а также на танкерном флоте.
СПИСОК ЛИТЕРАТУРЫ
1. Патент на полезную модель №115886, МПК G 01 F 1/86.Устройство для измерения уровня жидкости в резервуарах / Б.В. Скворцов, Д.И. Блинов, А.В. Солнцева, С.А. Борминский; заявитель и патентообладатель Самарский государственный аэрокосмический университет имени академика С.П. Королёва. -№2011151843/28; заявл.19.12.2011; опубл. 10.05.2012.
2. Патент на полезную модель №135121, МПК G 01 F 1/86. Лазерно - акустический уровнемер. / Б.В. Скворцов, Д.И. Блинов, А.В. Солнцева, С.А. Борминский; заявитель и патентообладатель Самарский государственный аэрокосмический университет имени академика С.П. Королёва; опубл. 27.10.2013.
3. Патент на изобретение №2332644, МПК G 01 F 23/ 292.Способ и устройство для измерения уровня жидкости / Леонидова З.Е., Мукаев Р.Ю.; заявитель и патентообладатель Уфимский государственный авиационный технический университет. №2007111271/ 28; заявл. 27.03.2007; опубл. 27.08.2008.
4. Патент на изобретение №2338163, МПК G 01 F 23/ 292. Способ и устройство для измерения уровня жидкости / Леонидова З.Е., Мукаев Р.Ю.; заявитель и патентообладатель Уфимский государственный авиационный технический университет. №2007121235/ 28; заявл. 06.06.2007; опубл. 10.11.2008.
5. Выгодский М.А. Справочник по высшей математике. М.: Наука, 1973. 872 с.
6. Метод многомерных тестовых объектов в оптических измерительных системах / В.Н. Нестеров, В.М. Мухин, А.В. Мещанов. Самара: СНЦ РАН, 2013. 224 с.
MATH MODEL AND TRANSFORMATION FUNCTION OF THE LASER-TELEVISION LIQUID LEVEL CONTROL METHOD
© 2014 B.V. Skvortsov, S.A. Borminsky, D.I. Blinov, D.M. Zhivonosnovskaya
Samara State Aerospace University named after Academician S.P. Korolyov (National Research University)
The present article contains theoretical substantiation and development of the laser-television liquid level control method by marks, which are generated on the controllable surface by one ore more laser beams with subsequent image processing.
Key words: math model, transformation function, laser-television method, liquid level control, laser radiation.
Boris Skvortsov, Doctor of Technical Sciences, Professor at the Electrical Engineering Department, Scientific Supervisor of Scientific Research Laboratory "Analytical Devices and Systems". E-mail: [email protected]
Sergey Borminsky, Candidate of Technical Sciences, Associate Professor at the Electrical Engineering Department. E-mail: [email protected]
Dmitrij Blinov, Postgraduate Student at the Electrical Engineering Department. E-mail: [email protected]. Darya Zhivonosnovskaya, Student. E-mail: [email protected]
