Научная статья на тему 'Моделирование процессов в системе датчик - объект при сканировании подводных коммуникаций'

Моделирование процессов в системе датчик - объект при сканировании подводных коммуникаций Текст научной статьи по специальности «Медицинские технологии»

CC BY
76
18
i Надоели баннеры? Вы всегда можете отключить рекламу.
Ключевые слова
ПОДВОДНЫЙ АППАРАТ / ИНСПЕКТИРОВАНИЕ / ДАТЧИК ОБЪЕКТ / ВЕКТОР / UNDERWATER VEHICLE / MONITORING / SENSOR OBJECT / VECTOR

Аннотация научной статьи по медицинским технологиям, автор научной работы — Черный Сергей Григорьевич, Жиленков Антон Александрович

Смоделированы процессы, описывающие взаимодействие сканирующего датчика с элементом линии подводной коммуникации. Данная задача актуальна при построении систем автоматизированного инспектирования подводных объектов с использованием автономных необитаемых подводных аппаратов. Получена простая модель, позволяющая описать структурные свойства материала объекта или повреждений на его поверхности. Математическая модель взаимодействия объекта и сканирующего датчика представлена в виде суммы двух рядов. Показано, что в задачах контроля локальной структуры поверхности объекта возможно представление математической модели в форме ряда Котельникова. Полученные результаты имеют высокую практическую значимость при решении задач прокладки специализированных линий и магистралей на территории Крымской акватории Азово-Черноморского бассейна и регулярного мониторинга их состояния

i Надоели баннеры? Вы всегда можете отключить рекламу.

Похожие темы научных работ по медицинским технологиям , автор научной работы — Черный Сергей Григорьевич, Жиленков Антон Александрович

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.
i Надоели баннеры? Вы всегда можете отключить рекламу.

Process modeling in sensor - object system when scanning underwater communications

The study focuses on the modeled processes describing the interaction of the scanning sensor with the element of underwater communication. This task is relevant in constructing the systems for automated monitoring of underwater objects by means of unmanned underwater vehicles. A simple model has been obtained that makes it possible to describe the structural properties of the object’s material or damage on its surface. The mathematical model of interaction between the object and the scanning sensor is presented as a sum of two series. The study shows that in problems of controlling the local structure of the object’s surface, it is possible to present the mathematical model in Kotelnikov expanded form. The obtained results are of high practical importance in solving the problems of laying specialized lines and highways in the territory of the Crimean water area of the Azov-Black Sea basin and regular monitoring their condition

Текст научной работы на тему «Моделирование процессов в системе датчик - объект при сканировании подводных коммуникаций»

УДК 519.977.54:004.94

С. Г. Черный, А. А. Жиленков

Моделирование процессов в системе датчик - объект при сканировании подводных коммуникаций

Смоделированы процессы, описывающие взаимодействие сканирующего датчика с элементом линии подводной коммуникации. Данная задача актуальна при построении систем автоматизированного инспектирования подводных объектов с использованием автономных необитаемых подводных аппаратов. Получена простая модель, позволяющая описать структурные свойства материала объекта или повреждений на его поверхности. Математическая модель взаимодействия объекта и сканирующего датчика представлена в виде суммы двух рядов. Показано, что в задачах контроля локальной структуры поверхности объекта возможно представление математической модели в форме ряда Котельникова. Полученные результаты имеют высокую практическую значимость при решении задач прокладки специализированных линий и магистралей на территории Крымской акватории Азово-Черноморского бассейна и регулярного мониторинга их состояния.

Ключевые слова: подводный аппарат, инспектирование, датчик - объект, вектор.

о см

см

О!

<

I

(0 та

0 ^

СО та

1

о.

3

и <и со

см ■ч-ю

с?

см ■ч-ю см

Введение

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

Для выполнения нормативных требований для визуального исследования линий подводных коммуникаций сегодня привлекаются водолазы или гораздо реже подводные аппараты с телеуправлением. В то же время в ведущих странах мира наблюдается тенденция к разработке эффективных методов и средств подобного контроля с помощью автономных необитаемых подводных аппаратов [1-3]. Теоретически использование подобных аппара-

- © Черный С. Г., Жиленков А. А., 2017

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

В научной литературе достаточно полно освещено решение задач автоматизированного поиска и распознавания подводных коммуникаций необитаемыми аппаратами [4-6], но информация о разработках в области автоматизированной инспекции целостности протяженных подводных объектов практически отсутствует. Известны буквально единицы аппаратов, обладающих подобной функцией [7-9]. На основании анализа доступных публикаций можно сделать вывод, что наибольших успехов разработчикам удалось достичь в области поиска и инспекции по электромагнитному излучению посредством соответствующих бортовых датчиков [9, 10]. Однако применение данного метода весьма ограничено узким кругом задач и имеет ряд существенных недостатков. Основными недостатками являются обеспечение обнаружения в ходе инспекции только мест разрыва запитанных электролиний и необходимость наличия последних в составе инспектируемой линии [7, 8].

В рамках решения задачи автоматизации инспектирования линий подводных коммуникаций, связывающих Республику Крым с материковой частью России, в данной статье исследуются проблемы распознавания дефектов на поверхностях элементов инспектируемых линий. В частности, проводится анализ изображений, получаемых с помощью телевизионных датчиков, подводных акустических сканеров (типа BlueViе^ BV 5ххх, Echoscope и т. п.) и их комплексирования.

Цели настоящей работы - построение модели взаимодействия датчика и объекта инспектирования, поиск на ее базе эффективных методов выделения основных особенностей структуры в изображении объекта, идентификации скрытых периодичностей, неравно-мерностей, т. е. решение задачи выделения существенных особенностей структуры изображения, характерных для элементов подводных коммуникаций. Математическое описание взаимодействия датчика и объекта контроля

Опишем систему объект - датчик в случае контроля поверхности объекта по его изображению в оптическом диапазоне (рис. 1).

f(x\,y\)

H{TUÖ

л

и(х, у) = Л /(х1, л)Я(х - х, у1 - у)(х(ух.

ан

Из-за расхождения области определения модели и окна теряется возможность применения свертки в виде

и(х, у) = Ц /(дц, У1)Н(х - х, у - у)йМу1.

ан

Действительно, отсчет координат модели и окна датчика не совпадают, что приводит к нарушению симметричности свертки. Данная особенность учтена в работе [11] и практически игнорируется в более поздних исследованиях. Это свойство присуще описаниям для телевизионных систем [12]. В этом случае для замены переменных является целесообразным переход к координатной системе X, У и переносной системе п, В результате сигнал от подвижного окна описывается уравнением

и(х, = \\ /(п + х, % + у)И(п, (1)

ан

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

и(х, ^) = /(х, ^) ® Н(х, у) = = Л /(п + х, \ + у)Н(п, пй\ =

S н(п Ъ

mn (am,n C0s п(П + X)'

n=0 n=0

(2)

x,

Рис. 1. Координатная система для подвижного окна: х, у - координаты переноса; /(х, у) - импульсная переходная функция объекта; Н - чувствительность датчика, по физическому смыслу совпадающая с локальной импульсной переходной характеристикой

В этом случае реакция и(х, у) определяется сверткой весовой функции материала объекта и функции распределения чувствительности датчика Н(х, у) по формуле:

хcosm (£, + y) + ... + bm n sin n(n + x) x x cos m (£, + y) + cm n cos n(n + x) x x sin m (£, + y) + ... ... + dmn sinn(n + x)sinm(£, + y))]dnd£..

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

f (x, y) = JJ H(n, $ X

aH

ю ю

XXX(am,n cosn(n +x)cosm(^ + y) +

и=0 я=0

+ bm n sin n(n + x)cos m(t, + y) + + cm n cos n(n + x)sin m(% + y) + + dm n sin n(n + x) sin m(% + y))] dnd

(3)

ra

ra 5 <D

H

H

Cd oo

fa

О

со О

I I

E

S S

0 ° л я

fe о 43 43

01 Ol

ti s er о о

"ы S о

3

43

s s

сг

о Sc H

Cd

О

оо S

0

*

1

о о H

er □

о

ti ^

X

s H er

0

s

01

s

о

Sc

s

о

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

43

ti s

ЕС

и

s

s

Ol

0 ti

01 ti er

Cd

s ti

s

OO

s S

о ^

N

I

о

Cd

er ti

fa

Cd

EC

Ol

EC

s

Sc

s s

s H сг о о о H fa

Cd о

S Я Й «

о S

fC

S

H

о

Cd

3

43 fa ti CT

S о

Sc

S

я

43

й а

О w

о 9 •р

Cd

S ti Ol

fa S

s §

ä m

H

E о S CD

43 fa t СГ

S о

Sc

£ 43

fa

fD 43 S

S Я S ti 9

+ cß

P S J о о

cß g

О

о

S

iTrs

J

P

'1 +

о о

cß S J

P

^

о о

s

iTrs

J

Ч

+

о о

s

iTrs

J

ч

ES S J о о

cß g

CO t—' • ES

S

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

iTr« о

0

1 +

О

о

cß S J

P

^

СО t—» •

P

0

1 +

Л

S

СО >—» •

S +

cß P

X

iTrs +

Ч

P S J о о

cß g

+ О

о

cß S J

P

g

СО >—» •

0

1 +

+

о о

cß S J

О

о

cß g

О

о

S

iTrs cß

4

ES

5 J

ES

^

о о

S if4« Cß

4

<Ггt J^rt О о

о

3 ъ о о

cß S

5s

+

P

Д

iTrs +

4

о

о о

cß S J

О

о

ä 8 s. S

Ез a

S J g

о

0

1

P S J

P

g

СО t—» •

P

S

iTrs о

0

1 +

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

Cß g

P S J

P

I,

X Р

S

iTrs о

0

1 +

О

0

S if4« Cß

P

1

I

P S J о о

cß g

СО >—» •

S if4« Cß

4

о

0

5 J

P

^

СО t—» •

P

S

iTrs cß

P

1

о-

S' s J о о

cß g

О

о

S

iTr« о

0

1 +

о о

cß S J

ЕЗ

g

о о

S

iTr« о

0

1

S

со » •

ЕЗ

s

3

+

о о

Д

iTr« +

Ч^ II

3' s J о о

cß g

+ О

о

cß S J

ЕЗ

g

о о

cß S J

О

о

cß g

ЕЗ S J

ЕЗ

g

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

3

Ol

о О cß S J

О

о

cß g

iTr«

ES

I +

ES S J

ES

^

CO t—» • es

S

iTr« cß

4

о о

cß S J

О

о

cß g

О

о

S

iTr« о

0

1

ЕЗ S J

ЕЗ

g

о о

S

iTr« о

0

1

о о

cß S

3

+

о о

iTr« +

4

о о

cß S J

О

о

cß g

ЕЗ S J

ЕЗ

g

о о

iTr« о о

cß '1 Cß

ЕЗ

S if4« Cß

ЕЗ

'1

Обозначив элементы, связанные с характеристиками датчика как компоненты, касающиеся соответствующих элементов спектральной матрицы модели материала, получаем:

К,п = JJ H(n, £)cos щ cos mtdndЪ,

ан

hbm,n = JJ H(n, £)sin nn cos mt,dnd, (5)

aH

hcm,n = JJ H(n, ^)cos nn sin mt,dnd%,

aH

hm,« = JJ H(n, %)sin nn sin m%dnd%.

aH

Подстановка (5) в уравнение (2) значительно упрощает модель реакции системы:

го го

Uу) = ^^^m,n (hm,nam,n + hm,nbm,n +

m=0 n=0

+ К >ncm>n + hi >ndm>n )c°s nx cos my +

+(—hb a + ha h — hd c + hc d ) X

'V m,n m,n m,n m,n 4,n m>n lm,n m,nJ

x sin nx cos my + (6)

+(—hc a — hd b + ha c + hb d ^ X

1 V "m,^m,n 4,num,n 1 "m.n^m.n 1 'lm,n m,nj

x cos nx sin my +

+(hd a — hc b — hb c + ha d ) X

V m,^wm,n 'm,^m,n ím,nc'm,n ím,^^ím,n/

x sin nx sin my.

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

В уравнении (6) компоненты структуры на изображении инспектируемой поверхности взаимодействуют с компонентами модели датчика, что обусловливает использование более высокой размерности. Введем векторы с компонентами спектральных характеристик датчика h и материала объекта a:

Ma =

Mb =

'1 0 0 0Л

0 10 0

0 0 10

0 0 0 1^

0 -10 0

V

1 0 0 0 0 0

0 0 0 -1 1 0

(8)

J

Mc =

^0 0 -1 0Л

0 0 0 -1

10 0 0

0 10 0

Md =

V

0 0 0 1 0 0 -10 0 -10 0 10 0 0

J

Данные матрицы в качестве элементов фактически содержат матрицы растяжения (сжатия) и матрицы поворота

Ma =

Mb =

'1 0 0 1

0

V

^0 -1 1 0

ha ' am,n

hb рm,n hc a — f m,n b ^ m,n c ^m,n

hd 1dm,n

Mc =

(7)

Сформируем матрицы, описывающие последовательность взаимодействия элементов модели датчика и изображения объекта [13]:

Md =

0 0

1 0 0 1

0

0 -1 1 0

0

1 0 0 1У

0

0 -1 1 0

-1 0

0 -1

(9)

0 1 -1 0

га

га

s ф

0

о см

см

OI

<

I

о та

s

о ^

со та г о. ф

£

и

V

со

см ■ci-io

с?

см ■ci-io см

(П (П

Матрица модели взаимодействия датчика и объекта в этом случае имеет вид:

(

Su —

m,n

^ ch m,n

am,nMbh m,n am,nMdhm,n J

(10)

Полученная модель позволяет сохранить форму изображения отклика как результат взаимодействия двух основных элементов - объекта и модели датчика в виде

го го

пу) .

т=0 п=0

К тому же значительно упрощается описание контроля поврежденных участков объекта. Благодаря тому, что выражение (8) имеет простое ядро, введение модели поврежденного материала объекта сводится к простым операциям над матрицами.

При использовании матриц модели деформированного материала и векторов возбужденной модели (9) можно определить реакцию системы объект - датчик (10) с помощью спектральных матриц

я =

" а + " иm,n am,n

bmn

m+ , Я m,n

^m,n ^'mm

d+ ^ m,n dm ,n

(12)

S"+ _

Си- _

^m,n

( Я + n---ah m,n a+m nMbh m Л n

у К n^^ch m,n а m nMd h m n J

am n---ah m,n am n-bh m Л n

a m nMd h m n J

(13)

Далее сигнал датчика в общем случае описывается уравнением

U(X У) = ХХ^шп{Ф+r. Sm+nФ^у) -

т=0п=0 го го

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

(14)

m=0 n=0

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

векторы как сумму собственного и совмещенного векторов

a+ m,n + d+mn

1 b+ m,n + c+ ~ m,n

2 с+ m,n +bm ,n

d+ ы m,n + a+m,n _

am,n —d— " **m,n

1 b~ ^ m,n —с— ^m,n

2 с — ^m,n —b- ^m,n

dm,n —a— Mm,n

1 + я

(15)

m,n'

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

С учетом заявленных технических характеристик сканирующих систем, использование которых возможно для рассматриваемых задач инспектирования, рассмотрим один из основных случаев - сканирование зондирующим лучом. В данном случае сканирование производится с помощью лучевого акустического сканера BlueView BV 5000, использующего элемент с Гауссовским распределением энергии и эффективным радиусом р(х, у) = = exp (-а2( х2 + у2) ).

Компоненты спектральной матрицы элемента разбиения имеют вид: ha = 1/ 2exp х х (-(m2 + n2) / 4а2). Эта спектральная матрица включает следующие элементы:

Ann = 1/2namn exp(-(m2 + n2)/4a2),

Bmn = 1 / 2nbmn exp(-(m2 + n2) /4a2), (16) X4 = XS = SIGN Y( XY),

Dmn = 1 / 2%dmn exp (-(m2 + n2)/4a2).

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

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

изображения поверхности объекта от заранее заданных.

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

fmn (x У) = ат>п cos nx cos my + + bm n sin nx cos my + cm n cos nx sin my +

+ dm n sin nx sin my = rm neJ(nx+my) + (17)

+ r eJ (nx-my) i r eJ (-nx+my) I ~'-m,n^ ' 'm,-nc' '

+ r e - J ( nx+my)

-m,-n

Учитывая, что в данном выражении (17) имеем отсчеты по узлам решетки со смещени-

3

ем по координатам на — п, запишем:

го го

f (x, y) = Z Z vJ(nx+my) +

m=0 n=0

гого

+Z Z r-m,nej(nx-my) +

m=0 n=0

гого

+Z Z m-nej(-nx+m) +

m=0 n=0

(18)

гого

+ Z Z "л> =

т=0 п=0

/\х, у) + /\х, у + А у) + Г(х + Ах, у) + /\х + Ах, у + Ау).

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

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

и(х, у) = Н /*(х, у) + /\х, у + А у ) + / •( х + А х, у) + (19) +1 Ч х + А х, у + А у )}.

Приняв, что сканируемая область является прямоугольником для равномерной чувствительности преобразователя в окне, можно записать:

U (X, J) = ZZ fm

sin (n(x / h - n))

-ZZ fm

n(x / h - n) sin(n(y /1 - m)) n(y /1 - m)

sin (n(x / h - n))

m+A,n

n(x / h - n)

sin(n(y /1 - m)) n(y /1 - m)

sin(n(x / h - n))

(20)

-ZZ fm

n+А

n(x / h - n) sin(n(y /1 - m))

-ZZ fm

n(y /1 - m)

sin (n(x / h - n))

m+А ,n +А

n(x / h - n) sin(n(y /1 - m))

X-.

n(y /1 - m)

Следовательно, сигнал сканирования формируется из четырех компонентов, представляющих собой ряд Котельникова. Отсчеты берутся в точках (x, y), (x + A, y), (x, y + A), (x + A, y + А) (рис. 2).

x + A,y + A

х,у + А

,х + А, у

Рис. 2. Координаты отсчетов функции

Таким образом, сохранение чувствительности модели к топологическим свойствам поверхности объекта достигается контролем сигнала в четырех точках области сканирования. Учитывая, что спектр сигнала ограничен, можно утверждать, что использование в качестве модели формулы (20) обеспечивает воспроизведение сигнала. Также можно сделать важный вывод: для контроля без поте- —

га

ри информации о структурных особенностях | поверхности объекта достаточно иметь окно § 2^2 пикселя, сравнимое с областью сканиро- £ вания (табл. 1).

х

х

Таблица 1

Минимальный размер окна сканирования

о см

см

О!

<

I

о га

г

о ^

со га г о.

V

3

и <и со

см ■ч-ю

с?

см ■ч-ю см

(П (П

/(У + А) /(X + А, ^ + А)

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

/(х, у) /(х + А, у)

Из анализа (14) следует, что линейно независимыми являются только четыре компоненты спектральных векторов. Для контроля в общем случае необходимо иметь достаточно простое окно, дополненное компонентами задержки (табл. 2).

Таблица 2

Структура окна сканирования

/(У + А) /(х + А, ^ + А)

/(х, у) /(х + А, у)

/(У - А) /(х + А, ^ - А)

Существенным является необходимость определения исходного шага в окне. Данный вопрос решается просто, если принять, что при контроле поверхности объекта достаточно проанализировать только первую гармонику. Отсчеты берутся в точках максимумов и нулей функции. Таким образом, достаточно просто построить модель, сохраняющую основные топологические свойства поверхности объекта контроля за максимумами функции и ее производных. Проверка адекватности математической модели Проверим сходимость модели при описании поверхности объекта. Использование при построении модели ряда Фурье позволяет заключить, что поскольку элементы сканируемого изображения физически не реализуют разрывов и структура изображения поверхности периодическая, то общая сходимость гарантирована. Важным остается вопрос сложности модели, поэтому необходимо оценить количество элементов ряда для получения модели с требуемой точностью.

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

1

Ц 1у

Модель первого приближения обеспечивает только описание основных компонентов, при этом средняя квадратическая погрешность достигает 30 %. При приближении двумя членами ряда погрешность снижается до 18 %. Как известно, увеличение числа членов ряда приводит к снижению ошибки. Однако, как показал анализ, переход к приближению с учетом пяти гармоник обеспечивает снижение ошибки до 4,8 %, при дальнейшем увеличении количества гармоник снижение ошибки замедляется так, что при использовании всего спектра удается получить ошибку не менее 2,7 %. Таким образом, использование при построении модели частот выше пятой гармоники можно признать нецелесообразным.

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

4ЬХЬ

| | I(х, y)dxdy. (21)

х У -Цх-Цу

Рис. 3. Изображение поверхности подводного кабеля

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

В спектре, кроме нулевой гармоники, доминируют только две составляющие, связанные с чередованием областей сканирования 1 и существованием возмущений структуры 2. Собственно, полный пространственный спектр изображения исследуемого образца достаточно сложный (рис. 5).

"0,0 _

Рис. 4. Одномерный спектр образца: Рис. 5. Пространственный

1 - область сканирования; 2 - возмущения структуры спектр изображения образца

а

у=2Ле+ОО4*х10-3.2е+ОО4*х9+1.9е+ОО5*х8+2.7е+ОО5*х7-6.8е+ОО5*хб-6 -8.2е+005*х5+1Ле+006*х4+1Ле+006*х3-7.9е+005*х2-4.9е+005*х+1.5е+005

Рис. 6. Строка (а) и столбец (б) матрицы автокорреляции

о

см

>s

ф

IX <

I

о (О

ц

о ^

ш

(О X

а

ф j

X

о *

S X I-

о ф

ш

см

ю о

см

ю см

О) О)

Реализация модели была выполнена в среде MATLAB. Исходный кадр изображения образца обрабатывался с выделением контуров. Выделение контуров изображения позволило подчеркнуть особенности периодичности структуры поверхности объекта контроля. Наличие корреляционного максимума позволяет установить основные частоты структуры материала. На рис. 6 приведены строка и столбец матрицы автокорреляции. Строка и столбец корреляционной матрицы являются автокорреляционными функциями для фиксированных координат.

Определить основные частоты структуры, являющиеся периодами структуры изображения, можно и с корреляционными функциями. Однако для корреляции характерно обострение не только основных частот, но и выделение скрытых периодичностей, что усложняет задачу определения основных параметров структуры модели [14-17]. Для получения значения основных частот модели найдем энергетический спектр изображения и его корреляционное поле (рис. 7).

Рис. 7. Энергетический спектр изображения образца (а) и корреляционное поле образца (б)

Рис. 8. Свертка модели поверхности инспектируемого объекта и эталонного образца

Определенные по энергетическим спектрам частоты позволяют построить эталонную полигармоническую модель поверхности инспектируемого объекта (рис. 8). Заключение

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

• математическая модель взаимодействия объекта (линии подводных коммуникаций) и сканирующего датчика (рассмотрен пример обработки изображения, полученного акустическим сканирующим датчиком) представляется в виде суммы двух рядов;

• компоненты спектральных векторов являются определенными только при отсутствии деформации структуры материала объекта;

• в задачах контроля локальной структуры поверхности объекта можно представлять математическую модель в форме ряда Котель-никова;

• для сохранения в модели топологических свойств объекта отсчеты необходимо выполнять в максимумах функции и ее производных.

Настоящие исследования имеют высокую практическую значимость ввиду прокладки специализированных линий и магистралей на территории Крымской акватории Азово-Черноморского бассейна и необходимости регулярного мониторинга их состояния. Список литературы

1. Matsumoto S., Ito Y. Real-time vision-based tracking of submarine-cables for AUV/ROV // OCEANS'95 MTS/IEEE. October, 1995. Vol. 3. Pp. 1997-2002.

2. Ortiz A., Simó M., Oliver G. A vision system for an underwater cable tracker // Machine Vision and Applications. July, 2002. Vol. 13, Iss. 3. Pp. 129-140.

3. Rives P., Borrelly ./../Underwater pipe inspection task using visual servoing techniques // Proceedings of the IEEE International Conference on Intelligent Robotic Systems. September, 1997. Vol. 7-11. Pp. 63-68.

4. Fukasawa T., Hozumi S., Morita M., Oketani T. Dissolved methane sensor for methane leakage monitoring in methane hydrate production // 0CEANS'06 MTS/IEEE. September, 2006. Vol. 18-21. CD-ROM.

5. Tracking // OCEANS'03 MTS/IEEE, Biloxi. September, 2003. Vol. 5. Pp. 2651-2658.

6. Evans J., Patryn P., Privat B., Johnson N., Capus C. AUTOTRACKER. Autonomous inspection - Capabilities and lessons learned in offshore operations // OCEANS'09 MTS/IEEE, Biloxi. October, 2009. Vol. 26-29.

7. Dubberley J.R., Gendron M.L. Synthetic aperture sonar low frequency vs. high frequency automatic contact generation // OCEANS'09 MTS/IEEE, Biloxi. October, 2009. Vol. 26-29. Pp. 1-4.

8. Copros T., Scourzic D. Pipeline inspection latest results from the Alistar // International Conference on Subsea Technologies, Saint Petersburg, Russia, June 2007.

9. Conte G., Zanoli S., Perdon A., Tascini G., Zingaretti P. Automatic analysis of visual data in submarine pipeline inspection // OCEANS' 96 MTS/IEEE. 1996. Vol. 3. Pp. 1213-1219.

10. Прэтт У. Цифровая обработка изображений. В 2 кн. Кн. 1. М.: Мир, 1982. 312 с.

11. Радионов В. А. Математическая модель текстильного материала и оптикоэлектронно-го преобразователя в потоке лучистой энергии // Известия высших учебных заведений. Технология текстильной промышленности. 1982. № 3. C. 82-87.

12. Системы технического зрения (принципиальные основы, аппаратное и математическое обеспечение) / А. Н. Писаревский,

А. Ф. Чернявский, Г. К. Афанасьев и др. ; под ред. А. Н. Писаревского, А. Ф. Чернявского. Л.: Машиностроение. Ленингр. отд-ние, 1988. 424 с.

13. Рыфтин К. В. Телевизионные системы. М.: Радио, Ленинград, 1968. 210 с.

14. Chernyi S., Zhilenkov A. Modeling of complex structures for the ship's power complex using XILINX system // Transport and Telecommunication. 2015. Vol. 16, no. 1. Pp. 73-82.

15. Chernyi S. Use of information intelligent components for the analysis of complex processes of marine energy systems // Transport and Telecommunication. 2016. Vol. 17, no. 3. Pp. 202-211.

16. Nyrkov A., Sokolov S., Zhilenkov A., Chernyi S., Mamunts D. Identification and tracking problems in qualimetry inspections in distributed control systems of drilling platforms // IEEE NW Russia Young Researchers in Electrical and Electronic Engineering Conference, Saint Petersburg, Russia, 2016. Pp. 641-645.

DOI 10.1109/EIConRusNW.2016.7448265

17. Жиленков А. А., Черный С. Г. Модель взаимодействия телевизионного датчика и объекта в задачах автоматизации инспектирования состояния подводных коммуникаций // Инженерная физика. 2016. № 4. С. 43-49.

Поступила 28.01.17

Черный Сергей Григорьевич - кандидат технических наук, заведующий кафедрой «Электрооборудование судов и автоматизации производства» Керченского государственного морского технологического университета, г. Керчь, доцент кафедры «Комплексное обеспечение информационной безопасности» Государственного университета морского и речного флота имени С. О. Макарова, г. Санкт-Петербург. Область научных интересов: экспертные системы, морские буровые платформы.

Жиленков Антон Александрович - кандидат технических наук, доцент кафедры систем управления и информатики Санкт-Петербургского национального исследовательского университета информационных технологий, механики и оптики, г. Санкт-Петербург. Область научных интересов: энергетические системы.

| MaTeMaTMKa |

Process modeling in sensor - object system when scanning underwater communications

The study focuses on the modeled processes describing the interaction of the scanning sensor with the element of underwater communication. This task is relevant in constructing the systems for automated monitoring of underwater objects by means of unmanned underwater vehicles. A simple model has been obtained that makes it possible to describe the structural properties of the object's material or damage on its surface. The mathematical model of interaction between the object and the scanning sensor is presented as a sum of two series. The study shows that in problems of controlling the local structure of the object's surface, it is possible to present the mathematical model in Kotelnikov expanded form. The obtained results are of high practical importance in solving the problems of laying specialized lines and highways in the territory of the Crimean water area of the Azov-Black Sea basin and regular monitoring their condition Keywords: underwater vehicle, monitoring, sensor - object, vector.

Cherny Sergey Grigorevich - Candidate of Engineering Sciences, Head of the Department of Electrical Equipment of Vessels and Manufacturing Automation, Kerch State Maritime Technological University, Kerch, Associate Professor, Department of Integrated Information Security, State University of Marine and River Fleet named after S. O. Makarov, Saint Petersburg.

Science research interests: expert systems and offshore drilling platforms.

Zhilenkov Anton Aleksandrovich - Candidate of Engineering Sciences, Associate Professor, Department of Control Systems and Informatics, Saint Petersburg National Research University of Information Technologies, Mechanics and Optics, Saint Petersburg. Science research interests: energy systems.

i Надоели баннеры? Вы всегда можете отключить рекламу.