Технологии раннего мониторинга технических объектов с применением оценок плотности распределения помехи
Технологии мониторинга технических объектов дают возможность по комбинациям оценок среднего квадратического отклонения, математического ожидания, плотности распределения, ее максимума и точек перегиба помехи каждого зашумленного сигнала в отдельности, а также всех поступающих сигналов вместе определ...
Gespeichert in:
Datum: | 2019 |
---|---|
Hauptverfasser: | , |
Format: | Artikel |
Sprache: | Russian |
Veröffentlicht: |
Інститут кібернетики ім. В.М. Глушкова НАН України
2019
|
Schriftenreihe: | Проблемы управления и информатики |
Schlagworte: | |
Online Zugang: | http://dspace.nbuv.gov.ua/handle/123456789/180838 |
Tags: |
Tag hinzufügen
Keine Tags, Fügen Sie den ersten Tag hinzu!
|
Назва журналу: | Digital Library of Periodicals of National Academy of Sciences of Ukraine |
Zitieren: | Технологии раннего мониторинга технических объектов с применением оценок плотности распределения помехи / Т.А. Алиев, Н.Ф. Мусаева // Проблемы управления и информатики. — 2019. — № 5. — С. 114-124. — Бібліогр.: 20 назв. — рос. |
Institution
Digital Library of Periodicals of National Academy of Sciences of Ukraineid |
irk-123456789-180838 |
---|---|
record_format |
dspace |
spelling |
irk-123456789-1808382021-10-21T01:26:26Z Технологии раннего мониторинга технических объектов с применением оценок плотности распределения помехи Алиев, Т.А. Мусаева, Н.Ф. Технические средства для измерений и управления Технологии мониторинга технических объектов дают возможность по комбинациям оценок среднего квадратического отклонения, математического ожидания, плотности распределения, ее максимума и точек перегиба помехи каждого зашумленного сигнала в отдельности, а также всех поступающих сигналов вместе определить начальный период зарождения и динамику развития неисправности, причем для каждого отдельно промышленного объекта применение предложенной технологии следует применять с учетом специфики исследуемого объекта. Показано, що технологія моніторингу технічних об’єктів із застосуванням комбінацій оцінок середнього квадратичного відхилення, математичного очікування, щільності розподілу, її максимуму і точок перегину завади кожного зашумленого сигналу окремо, а також всіх сигналів, що надходять, дають можливість визначити початковий період зародження несправності і динаміку її розвитку. Розглянуто специфіку застосування щільності розподілу завади зашумленого сигналу при відмінному від нуля математичному очікуванні для раннього моніторингу аварій різних технічних об’єктів. It is shown that the technologies for monitoring technical objects using combinations of estimates of the mean square deviation, mathematical expectation, distribution density, its maximum and inflection points of each noisy signal separately and of all incoming signals together make it possible to determine the initial period of initiation of malfunction and the dynamics of its development. The specifics of using the noise distribution density of a noisy signal with a non-zero mathematical expectation for early monitoring of accidents of various technical objects are considered. 2019 Article Технологии раннего мониторинга технических объектов с применением оценок плотности распределения помехи / Т.А. Алиев, Н.Ф. Мусаева // Проблемы управления и информатики. — 2019. — № 5. — С. 114-124. — Бібліогр.: 20 назв. — рос. 0572-2691 http://dspace.nbuv.gov.ua/handle/123456789/180838 519.216 ru Проблемы управления и информатики Інститут кібернетики ім. В.М. Глушкова НАН України |
institution |
Digital Library of Periodicals of National Academy of Sciences of Ukraine |
collection |
DSpace DC |
language |
Russian |
topic |
Технические средства для измерений и управления Технические средства для измерений и управления |
spellingShingle |
Технические средства для измерений и управления Технические средства для измерений и управления Алиев, Т.А. Мусаева, Н.Ф. Технологии раннего мониторинга технических объектов с применением оценок плотности распределения помехи Проблемы управления и информатики |
description |
Технологии мониторинга технических объектов дают возможность по комбинациям оценок среднего квадратического отклонения, математического ожидания, плотности распределения, ее максимума и точек перегиба помехи каждого зашумленного сигнала в отдельности, а также всех поступающих сигналов вместе определить начальный период зарождения и динамику развития неисправности, причем для каждого отдельно промышленного объекта применение предложенной технологии следует применять с учетом специфики исследуемого объекта. |
format |
Article |
author |
Алиев, Т.А. Мусаева, Н.Ф. |
author_facet |
Алиев, Т.А. Мусаева, Н.Ф. |
author_sort |
Алиев, Т.А. |
title |
Технологии раннего мониторинга технических объектов с применением оценок плотности распределения помехи |
title_short |
Технологии раннего мониторинга технических объектов с применением оценок плотности распределения помехи |
title_full |
Технологии раннего мониторинга технических объектов с применением оценок плотности распределения помехи |
title_fullStr |
Технологии раннего мониторинга технических объектов с применением оценок плотности распределения помехи |
title_full_unstemmed |
Технологии раннего мониторинга технических объектов с применением оценок плотности распределения помехи |
title_sort |
технологии раннего мониторинга технических объектов с применением оценок плотности распределения помехи |
publisher |
Інститут кібернетики ім. В.М. Глушкова НАН України |
publishDate |
2019 |
topic_facet |
Технические средства для измерений и управления |
url |
http://dspace.nbuv.gov.ua/handle/123456789/180838 |
citation_txt |
Технологии раннего мониторинга технических объектов с применением оценок плотности распределения помехи / Т.А. Алиев, Н.Ф. Мусаева // Проблемы управления и информатики. — 2019. — № 5. — С. 114-124. — Бібліогр.: 20 назв. — рос. |
series |
Проблемы управления и информатики |
work_keys_str_mv |
AT alievta tehnologiirannegomonitoringatehničeskihobʺektovsprimeneniemocenokplotnostiraspredeleniâpomehi AT musaevanf tehnologiirannegomonitoringatehničeskihobʺektovsprimeneniemocenokplotnostiraspredeleniâpomehi |
first_indexed |
2025-07-15T21:10:53Z |
last_indexed |
2025-07-15T21:10:53Z |
_version_ |
1837748826230751232 |
fulltext |
© Т.А. АЛИЕВ, Н.Ф. МУСАЕВА, 2019
114 ISSN 0572-2691
УДК 519.216
Т.А. Алиев, Н.Ф. Мусаева
ТЕХНОЛОГИИ РАННЕГО МОНИТОРИНГА
ТЕХНИЧЕСКИХ ОБЪЕКТОВ С ПРИМЕНЕНИЕМ
ОЦЕНОК ПЛОТНОСТИ РАСПРЕДЕЛЕНИЯ ПОМЕХИ
Ключевые слова: помеха с отличным от нуля математическим ожиданием,
зашумленный сигнал, плотность распределения, система контроля, диагностика
мониторинга.
Введение
Известно, что чрезвычайные ситуации возникают в результате негодности
производственных фондов, ветшания технологического оборудования, износа
систем противоаварийной защиты, нарушения производственного процесса,
отсутствия непрерывного контроля за состоянием технологического процесса,
нарушения технологии производства и т.д. [1]. К авариям приводят поврежде-
ния, поломки, коррозии, выход из строя агрегатов, двигателей, оборудования,
машин, станков, механизмов, разрушение строительных сооружений, систем
жизнеобеспечения зданий и т.д.
Например, основными причинами, приводящими к разрушению строи-
тельной конструкции, могут быть визуально невидимые дефекты, которые по-
являются даже при относительно слабых, но частых сейсмических воздействи-
ях на высотные здания, большепролетные конструкции, спортивные сооруже-
ния, строительные конструкции, энергетические объекты, промышленные
объекты и т.д. Кроме того, основными причинами аварий строительных объ-
ектов и конструкций являются также внешние техногенные (неравномерные
деформации основания) и природные воздействия (снег, ветер, вода, темпера-
тура и т.п.); превышение эксплуатационных нагрузок и износ; ошибки при
проектировании; экстремальные воздействия в результате пожара, взрыва и
т.д. [1].
При бурении основными видами аварий являются прихваты, поломка и от-
винчивание бурильных труб, поломка в скважине долот и турбобуров, падение
бурильного инструмента и других предметов в скважину. Подобные аварии могут
возникнуть при резком изменении гидравлического давления в скважине в ре-
зультате выброса, нарушения целостности ствола скважины, заклинивания бу-
рильной и обсадной колонн и т.д.
Аварии на объектах теплоснабжения, к которым относятся теплотрассы, ко-
тельные установки, когенерационные станции и т.д., происходят в результате по-
вреждений, коррозий, дефектов и т.д. в трубопроводах, тепловых сетях, котлах и
общекотельном оборудовании и т.д. [1].
Для того чтобы в вышеуказанных случаях исключить перерастание мел-
ких повреждений в крупномасштабную чрезвычайную ситуацию, необходимо
своевременно выявить возникшие неполадки на ранней стадии и провести
профилактические и ремонтные работы [1]. Для этого на технических объек-
тах совершенствуются технологические процессы, повышается надежность
технологического оборудования и систем, применяются конструкторские и
Работа выполнена при финансовой поддержке Научного фонда Государственной нефтяной компании
Азербайджана SOCAR.
Международный научно-технический журнал
«Проблемы управления и информатики», 2019, № 5 115
технологические решения, используется высококачественное сырье, материа-
лы, комплектующие изделия, используются и совершенствуются системы тех-
нологического мониторинга, контроля, технической диагностики, безаварий-
ной остановки производства, проводится локализация и подавление аварий-
ных ситуаций и т.д. Это особенно важно, например, на объектах нефтяной,
химической и газовой промышленности; на предприятиях, которые произво-
дят, используют, хранят и транспортируют взрыво- и огнеопасные вещества и
материалы, ядовитые вещества [1].
Однако перечисленные мероприятия и работы не дают возможность выявить
дефекты на начальном периоде их появления, когда они еще находятся в скрытом
виде. Поэтому для раннего обнаружения техногенных ситуаций и мониторинга
аварий необходимо разработать алгоритмы обнаружения раннего периода зарож-
дения дефекта, определить динамику развития неисправности, своевременно опо-
вестить обслуживающий персонал о необходимости остановки и проведения про-
филактических ремонтных работ.
В последние годы для решения этих задач разработаны методы, основанные
на оценке текущего технического состояния производственных объектов с ис-
пользованием характеристик помехи [2–15].
В данной работе предлагаются технологии раннего мониторинга техниче-
ских объектов на основе анализа оценок плотности распределения помехи,
имеющей отличное от нуля математическое ожидание. С помощью этой ин-
формации система мониторинга может формировать информацию о целесооб-
разности проведения необходимых мероприятий на раннем этапе возникнове-
ния неисправности.
1. Постановка задачи
Пусть от датчика, размещенного в зоне действия влияющих на объект факторов,
поступает зашумленный цифровой сигнал ).(tG При стабильном состоянии объекта
сигнал )(tG состоит из полезного сигнала )(tX и аддитивной помехи ),(1 t воз-
никшей в результате воздействия внешних факторов: )(tG = )(tX + ).(1 t В системе
мониторинга, контроля, диагностики, прогноза, управления, идентификации и
т.д. при появлении износа, коррозии, трещин, поломок, повреждений, выходов
из строя агрегатов, двигателей, механизмов, разрушений строительных соору-
жений, систем жизнеобеспечения зданий и т.д. также возникает помеха ).(2 t
Тогда зашумленный сигнал )(tG можно представить в виде )(tG = )(tX + ),(t
где )(t = )(1 t + ).(2 t
Сигналы )(tX , ),(t )(tG являются случайными стационарными эргодиче-
скими процессами, и помеху )(t невозможно выделить из ).(tG Для зашумленно-
го сигнала )(tG можно вычислить математическое ожидание ,Gm дисперсию ,GD
среднее квадратическое отклонение ,G корреляционную функцию )(GGR [16–9]:
N
i
G tiG
N
m
1
),(
1
,])([
1
1
2
N
i
GG mtiG
N
D ,GG D (1)
N
i
GG tiGtiG
N
R
1
),)(()(
1
)(
где t — шаг дискретизации, = 0, ,t ,2 t … — временной сдвиг.
116 ISSN 0572-2691
При этом сигнал )(tX оценивает текущее состояние исследуемого процесса.
Априори известно, что помеха )(t имеет нормальное распределение ),;( mN
и нулевое среднее .0m Кроме того, помеха )(t не коррелирована с полезным
сигналом, причем в момент зарождения дефекта время корреляции шума )(t
меньше времени корреляции полезного сигнала.
Поскольку стационарная случайная помеха )(t — эргодическая, то ее мате-
матическое ожидание m и среднее квадратическое отклонение имеют одно
и то же значение для любого из случайных сигналов, входящих в совокупность.
Поэтому функцию плотности нормального распределения )(),;( NmN
помехи представим
.
2
1
)(
2
2
2
)(
m
eN (2)
При этом следует отметить, что дисперсия D помехи )(t , входящая
в формулу (2), состоит из дисперсии
1
D помехи )(1 t от внешних факторов и
дисперсии
2
D помехи )(2 t в результате образования неисправности: D =
=
1
D + .
2
D Предполагается, что дисперсия
1
D помехи )(1 t от внешних
факторов постоянна или же изменения внешних воздействий поддаются кон-
тролю. Дисперсия
2
D помехи )(2 t , возникшей в результате дефекта, не
поддается контролю. Поэтому по виду и динамике изменения функции плот-
ности распределения )(N помехи )(t можно судить об изменениях, проис-
ходящих в техническом состоянии исследуемого объекта, так как различным
степеням неисправности соответствуют различные виды кривой )(N . Поэто-
му ниже предлагаются технологии раннего мониторинга различных конкрет-
ных технических объектов на основе анализа оценок плотности распределения
)(N помехи )(t с отличным от нуля математическим ожиданием.
2. Алгоритмы вычисления оценок плотности распределения помехи
с отличным от нуля математическим ожиданием для определения раннего
периода и динамики развития неисправностей технических объектов
Из формулы (2) очевидно, что для определения плотности распределения поме-
хи )(t необходимо вычислить среднее квадратическое отклонение D и
математическое ожидание m . Так как помеха )(t распределена по нормальному
закону с нулевым средним ,0m то задача сводится к вычислению параметра и
построению доверительного интервала для оценки математического ожидания.
Вычислим сначала дисперсию D помехи, используя выражение (1)
,)))(())(())(()((
1
)(
1
N
i
GG titiXtitiX
N
R
,)()( XmtiXtiX
N
i
X tiX
N
m
1
)(
1
— математическое ожидание ).(tX
Учитывая, что сигнал )(tX и помеха )(t некоррелированы, т.е.
,0))(()(
1
1
N
i
titiX
N
,0))(()(
1
1
N
i
tiXti
N
Международный научно-технический журнал
«Проблемы управления и информатики», 2019, № 5 117
можно написать [2–15]:
)(GGR = )(XXR + ),(R (3)
где
N
i
XX tiXtiX
N
R
1
))(()(
1
)(
— корреляционная функция ),(tX )(R
N
i
titi
N 1
))(()(
1
— корреляционная функция помехи )(t .
При этом на практике для таких инфранизкочастотных медленно протекающих
технологических процессов: нефтепереработка, нефтехимия, когда t значи-
тельно мало по сравнению со временем наблюдения ,T помеха )(t формируется
из высокочастотных спектров в результате возникновения таких неисправностей,
как износ, коррозия, нагарообразование и т.д. и имеет более высокий спектр, чем
сама полезная составляющая )(tX . Значение же полезной составляющей за вре-
мя t не успевает измениться, и )( ttX совпадает со значением ),(tX т.е.
)( ttX = )(tX . (4)
Это равенство выполняется для случаев, когда T составляет, например, 10–20 ча-
сов, а t — секунды или минуты (в зависимости от специфики исследуемого
процесса). В этом случае шаг дискретизации t выбирается исходя из конечного
времени корреляции cor помехи )(t с полезным сигналом. Очевидно, что такое
строгое равенство справедливо не для всех реальных процессов, а для таких, как
нефтепереработка, нефтехимия. Для остальных технологических процессов допу-
стимо приближенное равенство. Тогда для указанных объектов при выполнении
условия (4) отношение
)0(
)(
XX
XX
R
tR
равно единице, т.е. [16]:
).0()( XXXX RtR
В то же время в силу того, что для случайной помехи шаг дискретизации t
выбран, исходя из конечного времени корреляции cor помехи, то )(R можно
представить [2–15]
.при0
,0при)0(
)(
t
R
R (5)
Поэтому, если по формуле (1) вычислить оценки корреляционной функции
)(GGR зашумленного сигнала при =0 и достаточно малом по сравнению со
временем наблюдения T временном интервале tcor и найти разницу между
этими оценками, получим
)0(GGR – )( tRGG = )0(XXR + )0(R – )( tRXX – ).( tR (6)
Учитывая условия (5), (6), а также то, что оценки автокорреляционных функ-
ций соответственно полезного сигнала )(tX и помехи )(t при нулевом времен-
ном сдвиге =0 представляют собой дисперсии соответственно полезного сигна-
ла и помехи: )0(XXR = ,XD )0(R = ,D получаем
D = )0( GGR – ).( tRGG (7)
В работах [2–15] была выведена формула вычисления дисперсии помехи для
более общего случая, когда шаг дискретизации t выбран, исходя из частот-
118 ISSN 0572-2691
ной полосы спектра помехи )(t , т.е. ),2/(1 ft где f — частота среза по-
мехи (герц), а также имеется корреляция между полезным сигналом )(tX и
помехой )(t :
D = )0( GGR – 2 )( tRGG + ).2( tRGG (8)
После вычисления среднего квадратического отклонения
помехи можно
определить доверительный интервал для математического ожидания 0m по-
мехи. Известно, что доверительный интервал для оценки математического ожида-
ния при известном среднем квадратическом отклонении вычисляется по выраже-
нию [19]
N
zm
N
zm pp ; , (9)
где pz можно найти, задавая определенную доверительную вероятность p
);(1 z например, для вероятности p = 0,95 имеем ;96,195,0 z N — объем
выборки.
Тогда с учетом выражений (7)–(9) и условия 0m можно вычислить уточ-
ненную оценку математического ожидания
m помехи:
.
N
zm p
(10)
Принимая во внимание (2), (7)–(10), предлагается вычислять плотность рас-
пределения )(N нормально распределенной помехи )(t по выражению
,
2
1
))((
2
2
)(2
))((
mi
eiN (11)
где дискретные значения функции плотности распределения )(N помехи )(t
вычисляются в интервале ,3
m т.е. при .3)(3
mtm
Максимум функции плотности нормального распределения ))((max iN поме-
хи, который находится в точке ,m определяется по выражению
2
1
)(max mN
случая.общегодля
))2()(2)0((2
1
случая,частногодля
))()0((2
1
tRtRR
tRR
GGGGGG
GGGG
(12)
Точки перегиба с координатами
e
m
2
1
; и
e
m
2
1
;
вычисляются по выражениям:
— для первой и второй точек по оси абсцисс:
Международный научно-технический журнал
«Проблемы управления и информатики», 2019, № 5 119
1A
случая;общегодля))2()(2)0((
))()0((
случая,частногодля))()0((
))()0((
tRtRR
N
tRR
z
tRR
N
tRR
z
GGGGGG
GGGG
p
GGGG
GGGG
p
(13)
2A
случая;общегодля))2()(2)0((
))()0((
случая,частногодля))()0((
))()0((
tRtRR
N
tRR
z
tRR
N
tRR
z
GGGGGG
GGGG
p
GGGG
GGGG
p
(14)
— для первой и второй точек по оси ординат:
случая.общегодля
))2()(2)0((2
1
случая,частногодля
))()0((2
1
etRtRR
etRR
O
GGGGGG
GGGG (15)
3. Технологии раннего мониторинга технических объектов
на основе анализа оценок плотности распределения помехи
Для обнаружения появления дефекта на ранней стадии необходимо создать
банк данных, состоящий из оценок характеристик помехи, полученных в различ-
ные моменты времени ,it ,3,2,1i Для обработки этих данных предлагается
следующий алгоритм.
1. В момент времени ,1t когда объект находится в нормальном режиме
эксплуатации, по выражениям (7)–(15) вычисляются значения характеристик
помехи ,
1
t
D ,
1
t
,
1
t
m )),((1 iNt ),(max,1
mNt ,1 1
tA ,2 1
tA .1
tO Полученные
значения заносятся в банк данных в виде информативных признаков: 1tS
],21)())(([ 111max,11111
ttttt OEAEAEmNiNmD
ttt
где отмеча-
ется, что в техническом состоянии объекта нет неисправностей. Данные значения
фиксируются как характеристики помехи от внешнего воздействия и предполага-
ется, что внешние помехи находятся в неизменном состоянии.
2. Затем аналогичные характеристики помехи ,
2
t
D ,
2
t
,
2
t
m )),((2 iNt
),(max,2
mNt ,1 2
tAE ,2 2
tAE
2tOE вычисляются в момент времени 2t и также
заносятся в банк данных информативных признаков .2tS
Проводится анализ оценок помехи, полученных в моменты времени 1t
и ,2t и делаются выводы. Если
1t
D =
2t
D ,
1t
= ,
2
t
1t
m = ,
2
t
m ))((1 iNt =
= ))((2 iNt , )0(max,1
tN = ),0(max,2
tN
11tAE = ,1 2
tAE
12tAE = ,2 2
tAE
1tOE = ,2
tOE
то это означает, что техническое состояние объекта стабильное и неисправности
не наблюдаются.
Если
1t
D ≠ ,
2
t
D
1t
≠ ,
2
t
1t
m ≠ ,
2
t
m ))((1 iNt ≠ )),((2 iNt )0(max,1
tN ≠
≠ ),0(max,2
tN
11tAE ≠ ,1 2
tAE
12tAE ≠ ,2 2
tAE
1tOE ≠ ,2
tOE то имеется неисправ-
120 ISSN 0572-2691
ность. Причем если
1t
D >
2t
D ,
1t
>
2t
, то это означает, что имеется дефект,
его следует выявить и провести ремонтные работы. Если
1t
D ≫ ,
2
t
D
1t
≫
≫ ,
2
t
то это означает, что появился дефект, который интенсивно развивается.
Тогда для вектора информативных признаков 2tS устанавливается соответствие
между значениями характеристик помехи и степенью этого вида дефекта.
3. Затем характеристики помехи ,
3
t
D ,
3
t
,
3
t
m )),((3 iNt ),(max,3
mNt
,1 3
tAE ,2 3
tAE
3tOE вычисляются в момент времени ;3t полученные значения
также заносятся в банк данных информативных признаков .3tS Если
3t
D ≫
≫ ,
2
t
D
3t
≫ ,
2
t
то это означает, что повреждение интенсивно развивается.
Тогда в векторе информативных признаков этим значениям характеристик поме-
хи ставится в соответствие состояние «быстро развивающаяся степень дефекта».
В этом случае необходимо остановить работу оборудования во избежание воз-
никновения аварии.
4. Затем характеристики помехи вычисляются в моменты времени ,4t ,5t …
…, nt , и для каждого момента времени составляется банк информативных при-
знаков ,4tS ,5tS …, .tnS При этом для каждого из этих векторов устанавливается
соответствие между значениями характеристик помехи и степенью развития дан-
ного вида дефекта. Составляется матрица информативных признаков динамики
развития данного вида неисправности технического объекта:
tntntnttnn
ttttt
ttttt
ttttt
OAAmNNmD
OAAmNNmD
OEAEAEmNNmD
OAAmNNmD
S
tntn
ttt
ttt
ttt
21)()(
21)()(
21)()(
21)()(
max,3
333max,33
222max,22
111max,11
333
222
111
. (16)
5. Для других видов дефектов составляются аналогичные матрицы информа-
тивных признаков динамики развития неисправностей технических объектов. Это
дает возможность обнаружить момент возникновения дефекта и провести ранний
мониторинг аварий.
4. Примеры применения технологии раннего мониторинга технических
объектов на основе анализа оценок плотности распределения помехи
4.1. Технологии составления noise-матрицы раннего сейсмического мо-
ниторинга строительных объектов. Проведение сейсмометрического монито-
ринга строительных объектов — чрезвычайно важная задача, так как даже при от-
носительно слабых, но частых сейсмических воздействиях могут появиться визу-
ально невидимые дефекты, приводящие к разрушению конструкции. Ниже пред-
лагается следующий алгоритм раннего мониторинга. Составляется noise-матрица
информативных признаков (16). Затем после каждого маломощного землетрясе-
ния начинается поиск возникшего дефекта и устанавливается соответствие между
значениями характеристик помехи и видом дефекта. Одновременно в матрице
влияния сейсмических воздействий на техническое состояние строительного объ-
Международный научно-технический журнал
«Проблемы управления и информатики», 2019, № 5 121
екта выделяется i-я строка и отмечается тип и степень тяжести повреждения. Для
этого к noise-матрице (16) добавляют еще два столбца, первый соответствует виду
дефекта, а второй указывает на степень его повреждения. Виды дефекта записы-
ваются в закодированном виде: 1 — дефект несущей конструкции, в результате
которого теряется прочность и устойчивость; 2 — дефект ограждающей кон-
струкции, в результате которого сооружение ослабевает; 3 — дефект второсте-
пенного элемента конструкции, в результате которого снижаются эксплуатацион-
ные характеристики сооружений; 4 — дефекты механического происхождения,
например изменения грунтового основания, механические воздействия; 5 — де-
фекты физико-химического характера, например от коррозии, влажности, повы-
шенной температуры и т.д. Таким образом, получают расширенную noise-мат-
рицу технического состояния строительного объекта. При этом одному и тому же
дефекту могут соответствовать различные степени повреждения: 1 — незначи-
тельная, 2 — слабая, 3 — средняя, 4 — сильная, 5 — полное разрушение:
.
5
4
3
2
1
5
4
3
2
1
1
1
1
1
1
21)()(
21)()(
21)()(
21)()(
21)()(
21)()(
21)()(
21)()(
21)()(
21)()(
4,54,4,4,
3,43,3,3,
2,32,2,2,
1,21,1,1,
11111
55555
44444
33333
22222
11111
4,4,4,
3,3,3,
2,2,2,
1,1,1,
111
555
444
333
222
111
n
n
n
n
n
OAAmNNmD
OAAmNNmD
OAAmNNmD
OAAmNNmD
OAAmNNmD
OAAmNNmD
OAAmNNmD
OAAmNNmD
OAAmNNmD
OAAmNNmD
nttntntnt
nttntntnt
nttntntnt
nttntntnt
ttttt
ttttt
ttttt
ttttt
ttttt
ttttt
ntntnt
ntntnt
ntntnt
ntntnt
ttt
ttt
ttt
ttt
ttt
ttt
(17)
После каждого маломощного землетрясения проводится анализ технического
состояния строительного сооружения, и система мониторинга выдает сообщение
о мерах, которые необходимо предпринять в данном конкретном случае: 1 — ре-
монт не требуется; 2 — необходим текущий ремонт; 3 — провести капитальный
ремонт; 4 — заменить поврежденную строительную конструкцию; 5 — разобрать
и усилить строительную конструкцию.
4.2. Технологии составления noise-таблицы начала аварий в системе кон-
троля буровых установок. В настоящее время для управления процессом бурения
создаются системы контроля установки (СКУ). Однако существующие техноло-
гии контроля не обеспечивают бурильщика информацией, позволяющей выявить
начальный скрытый период возникновения аварий. Чаще всего этот момент уста-
навливается буровым мастером интуитивно.
Ниже предлагается технология определения ранней стадии зарождения де-
фекта с помощью характеристик (7)–(15) помехи. Пусть от датчиков, установлен-
ных на столе ротора, наверху буровой вышки и середине каркаса станка, посту-
пают вибрационные сигналы ),(1 tGV ),(2 tVG ),(3 tVG которые возникли от кру-
тильной, осевой и боковой вибраций, и передаются в компьютер. В компьютере
122 ISSN 0572-2691
по формулам (7)–(15) для каждого сигнала вычисляются характеристики помехи,
строится график плотности распределения помехи, которые затем выдаются на
экран: ,T мин: ,
t
D ,
t
,
t
m )),(( iNt ),(max,
mNt ,1tA ,2tA .tO
При изменении технического состояния буровой установки меняются зна-
чения среднего квадратического отклонения помехи и ее математическое ожи-
дание. Естественно, что при увеличении значения среднего квадратического
отклонения помехи максимум плотности распределения уменьшается, кривая
приобретает более плоскую форму. При изменении значения математического
ожидания помехи график плотности распределения помехи перемещается
вдоль оси абсцисс. При изменении формы кривой распределения меняются
координаты точек перегиба. В зависимости от значений характеристик помехи
для каждого отдельного зашумленного сигнала система сначала формирует
предварительное решение о состоянии буровой установки. Затем, анализируя
значения характеристик помехи для всех зашумленных сигналов в совокупно-
сти, система выносит окончательное решение. При этом для всех сигналов вы-
даются графики кривых плотности распределения с отмеченными точками пе-
региба и максимумами. Эти графики сравниваются с графиками, полученными
при нормальном функционировании объекта. Этим самым система оповещает
о произошедших изменениях и дает возможность буровому мастеру более
наглядно в этом убедиться. Одновременно система ставит в известность буро-
вого мастера о степени повреждения: 1 — повреждение несущественное и воз-
растает неинтенсивно и неравномерно; 2 — повреждение неощутимое и разви-
вается с небольшой интенсивностью; 3 — повреждение ощутимое и развива-
ется ускоренно; 4 — повреждение сильное и характеризуется интенсивным
развитием.
4.3. Технологии выявления ранней стадии возникновения неисправно-
стей объектов теплоснабжения. Известно, что одна из основных подсистем
энергетики — теплоснабжение. Для эффективного и надежного функционирова-
ния объектов теплоснабжения часто разрабатываются и внедряются системы
автоматизации, которые осуществляют контроль и сигнализацию отклонения
параметров, контроль за техническим состоянием трубопроводов, оперативное
оповещение о режимах работы тепловых сетей, мониторинг и контроль состо-
яния оборудования котлов и общекотельного оборудования, планирование ре-
монтных работ и оперативное реагирование на возникающие неполадки и т.д.
без остановки работы самого объекта на капитальный ремонт. Разработанные
алгоритмы вычисления характеристик помехи (7)–(15) рекомендуется также
применить для выявления ранней стадии возникновения неисправностей объ-
ектов теплоснабжения. Это особенно важно, например, при транспортировке
газа или другого горючего вещества.
Заключение
Технологии мониторинга технических объектов дают возможность по
комбинациям оценок среднего квадратического отклонения, математического
ожидания, плотности распределения, ее максимума и точек перегиба помехи
каждого зашумленного сигнала в отдельности, а также всех поступающих
сигналов вместе определить начальный период зарождения и динамику раз-
вития неисправности, причем для каждого отдельно промышленного объекта
применение предложенной технологии следует применять с учетом специфи-
ки исследуемого объекта [2–15, 20].
Международный научно-технический журнал
«Проблемы управления и информатики», 2019, № 5 123
Т.А. Алієв, Н.Ф. Мусаєва
ТЕХНОЛОГІЇ РАННЬОГО МОНІТОРИНГУ
ТЕХНІЧНИХ ОБ’ЄКТІВ ІЗ ЗАСТОСУВАННЯМ
ОЦІНОК ЩІЛЬНОСТІ РОЗПОДІЛЕНИХ ЗАВАД
Показано, що технологія моніторингу технічних об’єктів із застосуванням
комбінацій оцінок середнього квадратичного відхилення, математичного
очікування, щільності розподілу, її максимуму і точок перегину завади
кожного зашумленого сигналу окремо, а також всіх сигналів, що надхо-
дять, дають можливість визначити початковий період зародження неспра-
вності і динаміку її розвитку. Розглянуто специфіку застосування щільно-
сті розподілу завади зашумленого сигналу при відмінному від нуля мате-
матичному очікуванні для раннього моніторингу аварій різних технічних
об’єктів.
Ключові слова: завада з відмінним від нуля математичним очікуванням,
зашумлений сигнал, щільність розподілу, система контролю, діагностика
моніторингу.
T.A. Aliev, N.F. Musaeva
TECHNOLOGIES FOR EARLY
MONITORING OF TECHNICAL OBJECTS
USING THE ESTIMATES OF NOISE
DISTRIBUTION DENSITY
It is shown that the technologies for monitoring technical objects using combi-
nations of estimates of the mean square deviation, mathematical expectation,
distribution density, its maximum and inflection points of each noisy signal
separately and of all incoming signals together make it possible to determine
the initial period of initiation of malfunction and the dynamics of its develop-
ment. The specifics of using the noise distribution density of a noisy signal
with a non-zero mathematical expectation for early monitoring of accidents of
various technical objects are considered.
Keywords: Noise with non-zero mathematical expectation, noisy signal, distribution
density, control, diagnostics and monitoring system.
1. Aliev T. Noise control of the beginning and development dynamics of accidents. Switzerland :
Springer. 2019. 201 p. DOI: 10.1007/978-3-030-12512-7
2. Aliev T.A. Digital noise monitoring of defect origin. New York : Springer. 2007. 223 p. DOI:
10.1007/978-0-387-71754-8
3. Aliev T.A., Musaeva N.F., Suleymanova M.T. Density function of noise distribution as an indica-
tor for identifying the degree of fault growth in sucker rod pumping unit (SRPU). Journal of Au-
tomation and Information Sciences. 2017. 49, N 4. P. 1–11. DOI: 10.1615/JAutomatInf
Scien.v49.i4.10
4. Position-binary technology for processing the signals of a load for identification of the technical
state of deep-well sucker-rod pumps. T.A. Aliev, O.G. Nusratov, G.A. Guluev, As.G. Rzaev,
F.G. Pashaev, M.G. Rezvan. Measurement Techniques. 2018. 61, N 9. P. 885–890. https://doi.
org/10.1007/s11018-018-1519-x
5. Алиев Т.А., Мусаева Н.Ф., Сулейманова М.Т., Газызаде Б.И. Чувствительные алгоритмы
выявления степени развития неисправности штанговой глубинной насосной установки.
Мехатроника, автоматизация, управление. 2017. 18, № 2. С. 94–102. DOI: 10.17587/mau.
18.91-102
6. Aliev T.A., Musaeva N.F. An algorithm for eliminating microerrors of noise in the solution of
statistical dynamics problems. Automation and remote control. 1998. 59 (2), N 5. P. 679–688.
https://link.springer.com/article/10.1007/s11018-018-1519-x
https://link.springer.com/article/10.1007/s11018-018-1519-x
https://doi.org/10.17587/mau.18.91-102
https://doi.org/10.17587/mau.18.91-102
http://www.refdoc.fr/?traduire=en&FormRechercher=submit&FormRechercher_Txt_Recherche_name_attr=auteursNom:%20(ALIEV)
http://www.refdoc.fr/?traduire=en&FormRechercher=submit&FormRechercher_Txt_Recherche_name_attr=auteursNom:%20(MUSAEVA)
124 ISSN 0572-2691
7. Aliev T.A., Musaeva N.F., Suleymanova M.T., Gazizade B.I. Analytic representation of the den-
sity function of normal distribution of noise. Journal of Automation and Information Sciences.
2015. 47(8), N 4. P. 24–40. DOI: 10.1615/JAutomatInfScien.v47.i8.30
8. Aliev T.A., Musaeva N.F., Suleymanova M.T., Gazizade B.I. Technology for calculating the pa-
rameters of the density function of normal distribution of the useful component in a noisy process.
Journal of Automation and Information Sciences. 2016. 48, N 4. P. 35–55. DOI: 10.1615/
JAutomatInfScien.v48.i4.50
9. Aliev T.A., Musaeva N.F., Gazizade B.I. Algorithms of building a model of the noisy process by
correction of the law of its distribution. Journal of Automation and Information Sciences. 2017.
49, N 9. P. 61–75. DOI: 10.1615/ JAutomatInfScien.v49.i9.50
10. Aliev T.A., Musaeva N.F., Suleymanova M.T. Algorithms for Indicating the Beginning of
Accidents Based on the Estimate of the Density Distribution Function of the Noise of
Technological Parameters. Automatic Control and Computer Science. 2018. 52, N 3. P. 231–242.
DOI: 10.3103/S0146411618030021
11. Musaeva N.F. Robust method of estimation with «contaminated» coarse errors. Automatic Cont
rol and Computer Sciences. 2003. 37, N 6. P. 50–63. https:// elibrary.ru/contents.asp?
id=33405883
12. Aliev T.A., Musaeva N.F. Statistical identification with error balancing. Journal of computer and
systems sciences international. 1996. 34, N 5. P. 119–124
13. Aliev T.A., Musaeva N.F. Algorithms for improving adequacy of statistical identification. Jour-
nal of computer and systems sciences International, 1997. 36, N 3. P. 363–369. https://www.
tib.eu/en/search/id/olc%3A1518633188/Algorithms-for-Improving-Adequacy-of-Statistical
14. Aliev T.A., Musaeva N.F., Gazizade B.I. Algorithm of application of high-order moments of the
useful component as a diagnostic indicator of changes in the technical state. Journal of Automa-
tion and Information Sciences. 2018. 50, N 11. P. 29-43. https:// DOI: 10.1615/JAuto-
matInfScien.v50.i11.30
15. Aliev T.A., Musaeva N.F., Gazizade B.I. Algorithms for calculating high-order moments of the
noise of noisy signals. Journal of Automation and Information Sciences. 2018. 50, N 6. P. 1–13.
https:// DOI: 10.1615/JAutomatInfScien.v50.i6.10.
16. Техническая кибернетика. Кн. 2. Под ред. В.В. Солодовникова. М. : Машиностроение,
1967. 496 с.
17. Вентцель Е.С., Овчаров Л.А. Теория случайных процессов и ее инженерные приложения.
5-е изд. М. : КНОРУС, 2013. 448 с.
18. Пугачев В.С. Теория вероятностей и математическая статистика. 2-е изд. М. : Физматлит,
2002. 398 с.
19. Иванова В.М., Калинина В.Н., Нешумова Л.А., Решетникова И.О. Математическая стати-
стика. М. : Высшая шк., 1975. 398 с.
20. Abbasov A.M., Маmеdоva М.H., Orujov G.H.Aliyev H.B. Synthesis of the methods of subjective
knowledge representations in problems of fuzzy pattern recognition. Mechatronics. 2001. N 11.
P. 439–449.
Получено 18.06.2019
http://apps.webofknowledge.com/full_record.do?product=WOS&search_mode=GeneralSearch&qid=1&SID=F6arLwWxoPxgiiUzNY1&page=2&doc=13
http://apps.webofknowledge.com/full_record.do?product=WOS&search_mode=GeneralSearch&qid=1&SID=F6arLwWxoPxgiiUzNY1&page=2&doc=12
Journal%20of%20computer%20and%20systems%20sciences%20International,%201997.%2036,%20N 3.%20P.%20363–369%20
Journal%20of%20computer%20and%20systems%20sciences%20International,%201997.%2036,%20N 3.%20P.%20363–369%20
|