Математическая модель и функция преобразования лазерно-телевизионного метода контроля уровня жидкости

Автор: Скворцов Борис Владимирович, Борминский Сергей Анатольевич, Блинов Дмитрий Игоревич, Живоносновская Дарья Михайловна

Журнал: Известия Самарского научного центра Российской академии наук @izvestiya-ssc

Рубрика: Физика и электроника

Статья в выпуске: 6-1 т.16, 2014 года.

Бесплатный доступ

В настоящей статье предлагается теоретическое обоснование и разработка математической модели лазерно-телевизионного метода измерения уровня по меткам, создаваемым на контролируемой поверхности одним или группой лазерных лучей с последующей обработкой изображения.

Математическая модель, функция преобразования, лазерно-телевизионный метод, контроль уровня жидкости, лазерное излучение

Короткий адрес: https://sciup.org/148203536

IDR: 148203536

Текст научной статьи Математическая модель и функция преобразования лазерно-телевизионного метода контроля уровня жидкости

В настоящее время во всём мире наблюдается рост количества измерительных устройств, в которых применяется лазерное излучение. Измерения, проводимые с помощью лазерного излучения, содержат в себе ряд преимуществ. Использование лазерного излучения позволяет выполнять измерения с высокой чувствительностью, поскольку вследствие высокой энергетической яркости лазерного луча улучшаются энергетические соотношения в приёмных устройствах. Также, важным фактором в пользу применения лазерного излучения является то, что лазерный луч представляет собой направленное электромагнитное излучение, а скорость распространения электромагнитной волны мало зависит от окружающих факторов, таких как плотность среды, температура, влажность, давление. Известен ряд устройств [1, 2, 3, 4] для измерения уровня жидкости с помощью лазерного луча (одного или нескольких).

Выделим основные параметры измерительного лазерного луча, зависящие от контролируемого уровня жидкости:

-время распространения сигнала от поверхности до приёмника.

-амплитуда и фаза пришедшего на фотоприёмник сигнала.

-спектральный состав сигнала.

раметра с последующей обработкой используется в лазерно-телевизионном методе контроля уровня жидкости).

В данной статье предложены несколько вариантов реализации лазерно-телевизионного метода контроля уровня жидкости, а именно:

  • -    по одной лазерной метке;

  • -    по двум лазерным меткам, расстояние между которыми заведомо известно;

  • -    по N лазерным меткам, расстояние между которыми заведомо известно.

Для пояснения рассматриваемого метода рассмотрим лазерно-телевизионный уровнемер, представленный на рис. 1. Данный уровнемер содержит: телевизионную камеру 1, цифровое вычислительное устройство 2, источники лазерного излучения 3, 4, 5, 6, блок управления 7, блок индикации 8. По команде блока управления активируются источники лазерного излучения, на контролируемой поверхности появляются метки лазерных лучей, а фотоприёмная матрица регистрирует изображение контролируемой поверхности, которое поступает в цифровое вычислительное устройство, где по формулам определяется искомое расстояние H .

Сущность измерительной процедуры в самом общем случае иллюстрируется рис. 1, 2 и состоит в том, что на поверхности жидкости, уровень которой контролируется, с помощью лазерных лучей создаются метки A1 , A2 ,… Ai ,… An .

Каждая из этих меток имеет координаты в пространстве Ai ( Xi ; Yi ; Zi ) . Известно, что все эти точки принадлежат одной плоскости, уравнение которой можно определить через три произвольные точки Ai-1 , Ai , Ai+1 по формуле [5].

X-X i-1

Y - Y i -1

Z

-

Z

i -1

X i

-

X

I +1

X i -1

-

X i-1

Y i - Y i -1

Y + i - Y -i

Z i

-

Z + 1

Zi -1

-

Z i -1

= 0 . (1)

Рис. 1. Структурная схема лазерно-телевизионного уровнемера

Рис. 2. Иллюстрация к математической модели лазерно-телевизионного метода контроля уровня жидкости:

1 – фотоматрица, 2 – оптическая система, 3 – контролируемая поверхность

Выражение (1) преобразуется к виду

GX + QY + TZ + P = 0 ,      (2)

где

G = (Yi - Yi-,)(Zi+1 - Z-_x) - (Yi+1 - Yi-,)(Zi - ZM), t = (x, - X-1)(Y+1 - Y-1) - (Xw - Xi-1 )(Y - Y-1),

Q = - ( ( x - X -1 )( Z i +1 - Z i -1 ) -( - м - X -1 )( Z , - Z -1 ) ) ,

P = - GX i - 1 - QY - 1 - TZ i - 1

Предположим, что расстояние между всеми точками известно и определяется по формуле

L n = V ( - . - X ) + ( Y - Y )2 + ( Z . - Z )2 =

= л / A - 2 + A Y 2 + A Z 2 . n n n

Контролируемая поверхность снимается те-ле(фото)камерой.

Выбранные точки имеют координаты на изображении

А - 1 ( ц / - 1 ; n г - 1 )’ А ( М / ; n ii ), А + 1 (M1 / + 1 ; H + 1 ) .

Из законов геометрической оптики, по правилам сохранения пропорций, при съёмке со стороны Z (центр линзы находится на оси Z ) между координатами точек в пространстве X , Y , Z и на фотоматрице Ц, П справедливы следующие соотношения [6] (при Z f )

Ц =

Xf

Z - f

п =

Yf Z - f

где f –внешнее фокусное расстояние объектива, для стандартных объективов это величина известная, либо определяется экспериментально. Здесь Z – расстояние от контролируемой повер-

хности до главной плоскости оптической системы, которое связано с контролируемым уровнем Н и другими конструктивными параметрами формулой (см. рис. 1)

,        (6)

где Z m = L + B - высота установки центра оптической системы над нулевым уровнем жидкости, то есть при H = 0, Z = Z m .

Тогда выражения, связывающие координаты точек поверхности и их изображения на фотоматрице будут иметь вид

Ц = Xf          п = Yf .

Zm - H - f ’          2, - H - f ■ (7)

Рассмотрим варианты реализации лазернотелевизионного метода контроля уровня жидкости.

1. Контроль по одной метке. Это возможно, если метка создаётся лазерным лучом, строго параллельным главной оптической оси телекамеры и перпендикулярным контролируемой поверхности. При этом фокусное расстояние известно. В этом случае при смещении уровня вдоль оси Z координаты метки Xi и Yi не изменяются. Тогда уровень и координаты метки будут связаны соотношениями

Xif                     Yif

Ц Zm - H- f ’     п Zm - H- f (8)

mm

H = Z m - f (1 + ),           (9)

μ

„ . Y

H = Z m - f (1 + Y ) .          (10)

n

Из формул (8), (9) видно, что измерения можно проводить по любой координате метки. Для проведения процедуры измерений нужно при известных значениях Xi , Yi , f, задаться начальным значением координаты метки

Xif

Yif

П о Z - H -f" (11)

Z m H 0 f

При изменении уровня координаты метки A i ( ц i ; n i ) на изображении будут смещаться по фотоматрице (рис. 3).

Рис. 4. Перемещение меток на фотоматрице при изменении уровня

Рис. 3. Смещение координаты метки на фотоматрице в зависимости от уровня

Уровень и координаты меток связаны соотношениями (8). Тогда уровень можно найти из системы уравнений

xf n = Yf

Z - H - f , 11 Z - H - f mm

<

M n

Z m

X n f

-

H - f ’

П п

Y n f

Z m - H - f

L 2m = ( X n - X J2 + ( Y n - Y )

При этом каждому значению уровня будет соответствовать точка на фотоматрице. Единственным условием применения метода является захват телекамерой метки на контролируемой поверхности при минимальном и максимальном значениях уровня, текущее значение которого определяется по формулам (9), (10). Очевидно, что измерения будут более точными, если мы сделаем отсчёт по отдельным координатам, но результат усредним

, (12) где H M , H — уровень, определённый по формулам (9), (10) через соответствующие координаты. Контроль уровня по одной метке требует безусловного знания конструкционных параметров – фокусного расстояния f и высоты Zm, что не всегда возможно.

2. Контроль уровня по двум меткам, расстояние между которыми известно. Это возможно, если метки создаются лазерными лучами, строго параллельным главной оптической оси телекамеры и ортогональными контролируемой поверхности. В этом случае также при смещении уровня вдоль оси Z координаты Хi ,Yi меток не изменяются, Zi Zn = 0. В этом случае при изменении уровня будет изменяться расстояние между метками, (рис. 4).

где M^M n ; П —текущие координаты на изображении i -ой и n -ой меток соответственно. Из (13) получаем

M n - Mi

( X n - X ) f Zm - H - f ’

X n - X 1 =

П п - n =

( Y n - Y i ) f ; Z m - H - f ’

( M n - M 1 )( Z m - H - f )

f

,

Y - Y = ( П п - П 1 )( Z m - H - f ) n                            f ,

L n = ( Z m - H " f ) V ( M n - M i )2 + ( n n - П 1 )2, (14)

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

m

^^^^^B

f [1 +

,         Ln           ] = Z„

7( m - M i )2 + ( n - n )2

- f (1 + L n ) , (15) L n

где L n = V ( M n - M i J2 + ( n n - n i J2 — расстояние между метками на фотоматрице.

Выражение (15) является функцией преобразования прибора, позволяющей определять уро-

вень по координатам метки.

3. Контроль по N меткам, расстояние между которыми известно. Метки создаются лазерными лучами, параллельным главной оптической оси телекамеры и ортогональными контролируемой поверхности. Для трёх меток i, k, n система (13) будет иметь вид

I         X. ff            Y ff a i = ———, n = ——

- 1 Z m - H - f '* Z m - H - f

X f        Yf

A ----i---, n - —i---- i Zm - H - f    Zm - H - f ц - Xi+i f   n -   Y+i f

‘+1 Z m - H - f , ^n Z m - H - f

L 2 , _v- ( X -1 - X i )2 + ( Y -1 - Y )2

L 2 , +v - ( X i +1 - X i )2 + ( Y +1 - Y .)2

, L2 i -1, i +1 - ( X 1+1 - X i -1 )2 + ( Y +1 - Y -1 )2 .

Здесь значения Li-1, i , Li , i+1 , Li-1 , i+1 – считаются известными. Координаты меток на изображении являются данными измерительной процедуры. Целью решения является определение уровня Н . Система (16) путём последовательных подстановок сводится к виду

L - 1 ,.- ^Z m ff ^ ( A - 1 - A ) 2 + ( n - 1 - n ) 2 L,, + 1 - ^ Z m - H^ f ) V ( A + 1 - a )2 + ( n + 1 - n )2 (17) L i - 1, i + 1 - ( Zm H f ) V ( A + 1 - A )2 + ( n + 1 - n ) 2

Система (17) является вырожденной, так как размеры отрезков линейно связаны одним и тем же уравнением (9). По каждому из уравнений можно найти уровень Hi , но результат усреднить H - Z m - f [1 + ,       Li + 1        2 ] - Z m - f (1 + L ^) .(18)

V ( A +1 - A ) + ( n + 1 - n )                L i , i +1

H - N- ^ H,.          (19)

Усреднённый контроль уровня по нескольким меткам повышает точность измерений, в том числе компенсирует оптические искажения, если метки разбросаны по противоположным границам оптической системы.

Выражения (15)–(18) пригодны при горизонтальном положении контролируемого уровня, когда он ортогонален оптической оси. Реально оптическая ось может быть не ортогональна контролируемой поверхности, что имеет место быть при колебаниях резервуара, при этом уровень остаётся горизонтальным. Рассмотрим метод, позволяющий исключить влияние на результат измерений колебаний уровня относительно горизонта.

Выберем за ноль системы координат центральную точку объектива и оптическая ось объектива чтобы совпадала с осью Z . При выборе за начало координат других точек отсчета математическое описание процедуры формирования изображения сильно усложнится. Метки Аi на контролируемой поверхности формируются параллельными лазерными лучами, ортогональными к контролируемой поверхности (рис. 1). Расстояние между метками известно только при горизонтальной поверхности. При колебаниях уровня относительно оптической оси телекамеры оно изменяется. Из рис. 1 видно, что координаты Xi , Yi меток Аi(X ; Yi ; Zi) не зависят от уровня Н . Если за начало координат принять центр оптической системы, то их можно считать известными, определяемыми конструкционными размерами прибора. От уровня зависят только координаты Zi . Причём расстояние от каждой метки до поверхности будет разным и определяться через конструкционный параметр Zm по формуле

Z i =Z m – H i , (20) где Нi – уровень i -ойметки относительно горизонтального дна.

Возьмём четыре точки A,B,C,D, координаты XA , YA , XB , YB , XC , YC , XD, YD которых будем считать известными (рис. 5).

а)

б)

Рис. 5. Расположение лазерных меток: а – на контролируемой поверхности, б – на фотоприёмной матрице телевизионной камеры

Метки принадлежат одной плоскости, поэтому удовлетворяют уравнениям [5]

' P x - a + P y Y a + P z- Z a = 1 P x - B + P y Y B + P z" Z b = 1 P x X c + P y Y c + P z- Z c = 1 . P x - d + P y Y d + P z" Z d = 1 .

Здесь коэффициенты Px, Py, Pz , ZA , ZB , ZC , ZD , определяющие расположение плоскости в пространстве неизвестны. Однако известно [5], что при такой форме записи уравнения плоскости

P Z

Z Ц

где ZЦ – отрезок, отсекаемый плоскостью контролируемой поверхности по оси Z от начала координат, он и определяет расстояние от центра оптической системы до контролируемой поверхности и связан с искомым уровнем по формуле

H = Zm Zц . (23)

Учитывая, что оптические метки симметричны относительно начала координат, можно записать также но набрать нужное количество уравнений для определения неизвестных параметров, среди которых нас интересуют Z0, ZA , ZB , ZC , ZD , которые связаны с контролируемым уровнем по формуле (20), (23). Для контроля уровня по четырём точкам всего нам необходимо 9 уравнений. Исключив путём взаимных подстановок из уравнения (21) координату Z, получим систему из 4-х уравнений.

P XA + P Ya + P ( M a + X A ) f = 1

xAyAz

M a

P x X B + P y Y B + P z ( M b + X B ) f = 1

M b            (26)

P XC + P Y C + P ( M c + X C ) f = 1

Mc p .X +p .у +P ■(Md + XD)f = 1

  • 1    x л D + 1 y J D + 1 z                  1

Неизвестными и искомыми в этой системе являются Px , P y , Pz , f .

Система нелинейная. Несмотря на то, система имеет четыре неизвестных, она не имеет решения, так как является вырожденной относительно фокусного расстояния. Для решения проблемы введём обозначение

Z ц = 4( Z a + Z b + Z c + Z d ) .   (24)

Обратимся к проекции меток на матрицу (рис. 5, б). Точки проекции будут иметь коорди-

на™ А ( ц A ; n A ) , B ( ц B ; n B ) , C ( Ц с ; П с ) -u D ( h d ; П D ) • Учи тывая (5) для координат проекций каждой точки

-

f

Q = Pzf = ^.       (27)

Тогда система (25) сводится к трём уравнениям

можно записать

M a =

XАf

z а - Г

П а =

YАf

P x X a + P y Y a + Q (1 + - A ) = 1 M a

M b

Xf B

Zb f ’

B

M c Zc " f

M d   Z d " f

Z а - f ,

П а =

Yf B

Z„ f , B

Y C f

П А   Z c f

П а =

YDf

Z D - f .

Оптическая система может состоять не из одной линзы, однако главное её назначение – пропорциональная передача изображения. Поэтому соотношения (24) универсальны и применимы для любых сколь угодно сложных оптических систем.

Параметры Ц A , П A , Ц B , n B , Ц C , П C , Ц D , П D являются данными измерительной процедуры, определяются по изображению и считаются известными.

Из совокупности соотношений (21) – (25) мож-

j P x - B + P y Y B + Q (1 + - B ) = 1      (28)

M b

P x - C + P y Y c + Q (1 + - C ) = 1 M c

Это линейная система, в которой нас интересует Q , имеет прямое решение по формулам Кра-

мера.

Q =

Л Q

А

где

X A

Y A

+ - A

M a

А = -

X C

Y A

Y C

X B

M b

X

+ — C

M c

Y

, А Q

X B

X C

Y A

Y C

. (30)

Дальше при известном фокусном расстоянии из (26) получаем

f

Z Ц - Q.             (31)

Отметим, что при известном фокусном расстоянии величина Z0 может быть найдена из соотношения (24), которое с учётом (25),можно переписать в виде

X X X X

Z Ц - 0,25 f (4 +   + - ^ + - ^ + -D-) . (32)

μAμBμCμD

Искомый уровень нефтепродукта определяется по формуле (23).

Выражения (23), (29) – (32) определяют уровень колеблющейся относительно горизонта поверхности по координатам полученных на фотоматрице оптических меток и в совокупности определяют функцию преобразования прибора. Результаты, полученные по формулам (31),(32) можно усреднить. Выражение (31) предпочтительнее, так как оно не накладывает ограничений на симметричность лазерных меток относительно центра оптической системы и позволяет контролировать уровень по трём произвольным меткам, созданными на поверхности лазерными лучами, параллельными главной оптической оси.

Таким образом, общий алгоритм решения задачи определения уровня жидкости по четырём оптическим меткам в случаях горизонтальных колебаний уровня состоит в следующем:

  • 1.    Для конкретных значений координат меток XA,YA,XB,YB, XC,YC решаем систему уравнений (28), из которой находим параметр Q .

  • 2.    По выражениям (31) – (32) находим искомый уровень.

  • 3.    При необходимости результат усредняем

Описанные в статье оптические методы позволяют повысить точность измерения уровня жидкости, обеспечить контроль за наклоном резервуара, что дает возможность применять данные разработки на железнодорожном, автомобильном транспорте, а также на танкерном флоте.

Список литературы Математическая модель и функция преобразования лазерно-телевизионного метода контроля уровня жидкости

  • Патент на полезную модель №115886, МПК G 01 F 1/86.Устройство для измерения уровня жидкости в резервуарах/Б.В. Скворцов, Д.И. Блинов, А.В. Солнцева, С.А. Борминский; заявитель и патентообладатель Самарский государственный аэрокосмический университет имени академика С.П. Королёва. -№2011151843/28; заявл. 19.12.2011; опубл. 10.05.2012.
  • Патент на полезную модель №135121, МПК G 01 F 1/86. Лазерно-акустический уровнемер./Б.В. Скворцов, Д.И. Блинов, А.В. Солнцева, С.А. Борминский; заявитель и патентообладатель Самарский государственный аэрокосмический университет имени академика С.П. Королёва; опубл. 27.10.2013.
  • Патент на изобретение №2332644, МПК G 01 F 23/292.Способ и устройство для измерения уровня жидкости/Леонидова З.Е., Мукаев Р.Ю.; заявитель и патентообладатель Уфимский государственный авиационный технический университет. №2007111271/28; заявл. 27.03.2007; опубл. 27.08.2008.
  • Патент на изобретение №2338163, МПК G 01 F 23/292. Способ и устройство для измерения уровня жидкости/Леонидова З.Е., Мукаев Р.Ю.; заявитель и патентообладатель Уфимский государственный авиационный технический университет. №2007121235/28; заявл. 06.06.2007; опубл. 10.11.2008.
  • Выгодский М.А. Справочник по высшей математике. М.: Наука, 1973. 872 с.
  • Метод многомерных тестовых объектов в оптических измерительных системах/В.Н. Нестеров, В.М. Мухин, А.В. Мещанов. Самара: СНЦ РАН, 2013. 224 c.
Еще
Статья научная