ISSN 0236-235X (P)
ISSN 2311-2735 (E)

Journal influence

Higher Attestation Commission (VAK) - К1 quartile
Russian Science Citation Index (RSCI)

Bookmark

Next issue

4
Publication date:
09 December 2024

The system of statistical analysis and control of hydraulic unit vibration stability

Date of submission article: 20.11.2017
UDC: 519.248:658.518.5
The article was published in issue no. № 3, 2018 [ pp. 620-625 ]
Abstract:The system of vibration monitoring of a hydraulic unit includes indications of a distributed network of relative and absolute vibration intelligent sensors, as well as sensors for measuring the shaft rotation speed. These sensors are designed to work as a part of the hydraulic unit control system. The sensors are combined into one local network to transfer data on vibration parameters to the control system. The network operation mode during stationary operation of the unit can use multivariate statistical control algorithms in order to estimate vibration stability. These algorithms are similar to those used for controlling multiparameter technological processes. Such approach ensures early detection of process abnormalities using statistical methods and timely taking of control actions. Based on the conducted research and result analysis of vibration monitoring of the hydroelectric unit in the Krasnopolanskaya hydroelectric power station, there is a developed software system that provides monitoring of the vibration stability. The system in-cludes three subsystems: for visualization and grouping indicators, as well as for monitoring independent and correlated groups of indicators by stability criteria of the average level and the scattering the corresponding processes. In order to control independent in-dicators, the authors use Shewhart charts for mean values and standard deviations. When controlling correlated indicators, Hotelling charts are used to monitor the mean process level and a generalized variance chart for monitoring multidimensional scattering. The system ensures timely detection of possible abnormalities and prevents overranging of controlled indicators.
Аннотация:Система вибромониторинга гидроагрегата включает показания распределенной сети интеллектуальных датчиков относительной и абсолютной вибраций, а также датчиков измерения скорости вращения вала. Эти датчики предназначены для работы в составе системы управления гидроагрегатом. Датчики объединяются в единую локальную сеть для передачи данных о параметрах вибраций в систему управления. В режиме работы в сети при условии стационарной работы агрегата для оценки стабильности вибраций могут быть использованы алгоритмы многомерного статистического контроля, аналогичные применяемым при контроле многопараметрических технологических процессов. Такой подход обеспечивает раннее выявление нарушений процесса статистическими методами и своевременное принятие управляющих воздействий. На основе проведенного исследования и анализа результатов вибромониторинга гидроагрегата на Краснополянской ГЭС разработана программная система, обеспечивающая контроль стабильности вибраций. Система включает три подсистемы: для визуализации и группировки показателей, для контроля независимых и коррелированных групп показателей по критериям стабильности среднего уровня и для рассеяния соответствующих процессов. Для контроля независимых показателей используются карты Шухарта для средних значений и стандартных отклонений. При контроле коррелированных показателей применяются карты Хотеллинга для мониторинга среднего уровня процесса и карта обобщенной дисперсии для мониторинга многомерного рассеяния. Система позволяет своевременно выявлять возможные нарушения и обеспечивает предотвращение выхода контролируемых показателей за пределы допустимых значений.
Authors: Klyachkin, V.N. (v_kl@mail.ru) - Ulyanovsk State Technical University (Professor), Ulyanovsk, Russia, Ph.D, Yu.E. Kuvayskova (v_kl@mail.ru ) - Ulyanovsk State Technical University (Associate Professor), Ulyanovsk, Russia, Ph.D, Ivanova A.V. (tiki.mikck@yandex.ru) - NPO RusBITex (Junior Researcher), Tver, Russia
Keywords: vibration, statistical control, hotelling’s chart, generalized variance, diagnosis of violation
Page views: 6106
PDF version article
Full issue in PDF (29.03Mb)

Font size:       Font:

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

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

Характеристики вибраций в различных зонах гидроагрегата могут быть взаимосвязаны и образуют систему временных рядов. Можно построить соответствующие математические модели как в режиме запуска агрегата, так и в режиме работы в сети. На базе этих моделей возможны прогнозирование изменений характеристик агрегата и обнаружение нарушений процесса до того, как контроли- руемые параметры превысили предельные значения [2, 3].

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

Постановка задачи

Для обработки информации и разработки методики анализа и контроля стабильности вибраций использовались результаты записи показателей в системе управления гидроагрегатом на Краснополянской ГЭС (Краснодарский край) в режиме работы в сети, при котором контролируемые показатели должны быть стабильными. Проводился мониторинг биения валов и вибраций гидроагрегата с периодичностью два измерения в секунду. Процесс определяется десятью показателями: Х1–Х2 – вибрации нижнего генераторного подшипника верхнего бьефа и на правом берегу; Х3–Х4 – соответствующие вибрации верхнего генераторного подшипника; Х5–Х6 – бой вала гидротурбины, нижний бьеф и правый берег соответственно; Х7–Х8 – бой вала гидрогенератора; Х9–Х10 – вибрации крышки гидротурбины.

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

Для выявления участков установившейся (стационарной) работы гидроагрегата вибрации усреднялись по времени и визуализировались. На рисунке 2 в качестве примера показаны вибрации при усреднении данных по 10 сек. Начальная часть участка стабильной работы (25–30 усредненных данных) используется для оценки характеристик процесса: среднего уровня, рассеяния, корреляций, положения контрольных границ.

Статистический контроль многопараметри- ческого процесса предполагает выделение из множества контролируемых показателей групп кор- релированных параметров: для них необходимо применение многомерных методов, основанных на алгоритме Хотеллинга [5, 6]. Кроме того, статистический контроль ориентирован на нормально распределенные данные, однако, имея в виду большой объем измерений и усреднение данных, в условиях рассматриваемой задачи это обстоятельство может не учитываться.

Предварительный анализ данных

	Х1	Х2	Х3	Х4	Х5	Х6	Х7	Х8	Х9	Х10Х1	1									Х2	-0,01	1								Х3	-0,14	-0,03	1							Х4	-0,18	-0,01	0,10	1						Х5	-0,01	-0,02	0,15	-0,17	1					Х6	0,03	0,38	-0,02	0,10	0,12	1				Х7	-0,07	0,10	-0,20	0,35	-0,81	-0,07	1			Х8	0,03	0,44	-0,12	0,22	-0,13	0,61	0,31	1		Х9	-0,06	0,14	0,05	0,18	-0,24	0,09	0,31	0,20	1	Х10	-0,01	0,04	-0,08	0,18	-0,10	0,05	0,22	0,09	0,21	1Рис. 3. Корреляционная матрица (выделены значимые корреляции)Fig. 3. A correlation matrix (significant correlations are highlighted)Для анализа взаимосвязей между контролируемыми показателями строилась корреляционная матрица и проверялась значимость корреляций по критерию Стьюдента. Корреляционная матрица показана на рисунке 3, значимые корреляции выделены жирным шрифтом. Имеют место значимая отрицательная корреляция r = – 0,81 между показателями Х5 и Х7 (боем валов гидрогенератора и гидротурбины на нижнем бьефе) и положительная корреляция r = 0,61 между Х6 и Х8 (аналогичными показателями на правом берегу). Следует отметить, что корреляции зависят от того, как проводится усреднение данных (на рисунке приведены данные при усреднении по 10 сек.).

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

1.    Независимые показатели Х1–Х4 и Х9–Х10. Для контроля этих шести показателей могут быть использованы стандартные карты Шухарта для средних значений и стандартных отклонений.

2.    Попарно коррелированные показатели Х5–Х7, Х6–Х8. Эти две пары показателей следует контролировать по стабильности среднего уровня с использованием карт Хотеллинга, а по контролю рассеяния – карт обобщенной дисперсии.

Статистический контроль независимых показателей

Контрольные карты Шухарта для показателя Х1 (вибрации нижнего генераторного подшипника верхнего бьефа), построенные в системе Statistica на этапе анализа для 30 выборок [7], показаны на рисунке 4. Видно, что процесс стабилен: нет выходов контролируемых значений за границы ни на одной из карт. Полученные значения среднего уровня процесса, среднего размаха и контрольных границ могут быть использованы для дальнейшего мониторинга процесса по этому показателю. Подобные карты были построены и для других независимых показателей.

При мониторинге процесса для оценки стабильности вибраций, кроме выходов точек за границы соответствующих доверительных интервалов, использовались и другие критерии, связанные с наличием так называемых неслучайных структур: наличие 9 точек подряд по одну сторону от центральной линии, 6 возрастающих или убывающих точек подряд (тренд), 14 попеременно возрастающих и убывающих точек и других [4].

Многомерный контроль среднего уровня процесса

Статистический контроль коррелированных показателей, как и независимых, может проводиться по среднему уровню вибраций и их рассеянию. Для контроля среднего уровня многопараметрического процесса используется алгоритм Хотеллинга, а многомерное рассеяние характеризуется обобщенной дисперсией – определителем ковариационной матрицы [6, 8].

На этапе анализа процесса в условиях его стабильности по результатам исследования 30 выборок для каждой группы коррелированных показа- телей (в рассматриваемом процессе таких групп две: Х5–Х7 и Х6–Х8) оцениваются следующие показатели:

– вектор средних m0 = (m1, ..., mp)T,

mj =                               (1)

где р – количество контролируемых показателей; m – количество выборок; n – объем выборки; xijt – результат i-го наблюдения по j-му показателю в t-й выборке;

– ковариационная матрица, компоненты которой sjk = ,    (2)

j, k = 1, …, p                               

(определитель этой матрицы используется как оценка целевой обобщенной дисперсии |Σ0|);

– контрольные границы для карт Хотеллинга и обобщенной дисперсии.

Положение контрольной границы карты Хотеллинга определяется формулой

Tkp2 =         (3)

где  – квантиль F-распределения Фишера с числами степеней свободы в числителе k1, в знаменателе – k2 порядка 1 – α, a – значение уровня значимости (вероятность ложной тревоги).

Верхняя и нижняя контрольные границы карты обобщенной дисперсии:

.                                    (4)

Коэффициенты b1 и b2 вычисляются по формулам:

 ,                               (5)

. (6)

Если значение на нижней границе окажется меньше нуля, она принимается равной нулю.

Найденные характеристики используются при последующем мониторинге вибраций. На рисунке 5 показана карта Хотеллинга при мониторинге группы показателей Х5–Х7: на данном этапе наблюдений имеет место выброс за контрольную границу в выборке 13.

Возможные нарушения процесса анализировались также на наличие неслучайных структур (тренд, цикличность, приближение к оси абсцисс или контрольной границе, резкие скачки), а также с использованием предупреждающей границы [9]. Были выявлены резкие скачки в выборках 1–2, 10–11, 12–13–14, 19–20 и 22–23–24, свидетельствующие о возможном нарушении процесса.

Заметим, что эти скачки, как и выход точки за контрольную границу, свидетельствуют о нарушении статистической стабильности среднего уровня процесса по группе показателей Х5–Х7, но не о выходе соответствующих характеристик за пределы допустимых значений. Принятие необходимых мер управленческого воздействия обеспечит стабилизацию процесса.

Одновременно проводился и мониторинг многомерного рассеяния по этим же группам показателей. На рисунке 6 показана карта обобщенной дисперсии, полученная для той же группы показателей Х5–Х7 в то же время, что и карта Хотеллинга (рис. 5). Эта карта свидетельствует о стабильности процесса по критерию многомерного рассеяния.

Структура системы статистического анализа и контроля стабильности вибраций гидроагрегата

Система (структурная схема представлена на рисунке 7) предназначена для анализа и контроля вибраций гидроагрегата при его стационарной работе в сети с целью своевременного выявления нарушений стабильности и предотвращения выхода вибраций за допустимые границы. Система включает три программы.

Программа для предварительного анализа показаний интеллектуальных датчиков контроля вибраций дает возможность усреднить и визуализировать показания, оценить коррелированность между показаниями различных датчиков и разделить их на группы независимых и коррелированных с целью дальнейшего статистического контроля. Интерфейс программы на этапе усреднения данных (см. http://www.swsys.ru/uploaded/image/2018-3/ 2018-3-dop/19.jpg), графики вибраций (рис. 2) построены в этой же программе.

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

Программа контроля стабильности вибраций по коррелированным показателям предназначена для анализа и мониторинга с применением карты Хотеллинга для среднего уровня многомерного процесса и карты обобщенной дисперсии. При этом на этапе анализа оцениваются соответствующие характеристики по формулам (1)–(6), используемые при мониторинге коррелированных показателей. При мониторинге процесса также выявляются структуры специального вида, кроме того, используется предупреждающая граница для выявления нарушений по карте Хотеллинга [9]. Интерфейс программы на этапах выбора параметров и расчетов по алгоритмам Хотеллинга и обобщенной дис- персии представлен на рисунках (см. http://www. swsys.ru/uploaded/image/2018-3/2018-3-dop/20.jpg, http://www.swsys.ru/uploaded/image/2018-3/2018-3-dop/21.jpg), а построенные в программе контрольные карты – на рисунках 5 и 6.

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

Заключение

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

Работа выполнена при поддержке гранта РФФИ № 16-48-732002.

Литература

1.     Владиславлев Л.А. Вибрация гидроагрегатов гидроэлектрических станций. М.: Энергия, 1972. 153 с.

2.     Клячкин В.Н., Кувайскова Ю.Е., Алешина А.А. Моделирование вибраций гидроагрегата на основе адаптивных динамических регрессий // Автоматизация. Современные технологии. 2014. № 1. С. 30–34.

3.     Кувайскова Ю.Е., Алешина А.А. Повышение эффективности системы управления техническими объектами при использовании адаптивного динамического регрессионного моделирования временных рядов // Автоматизация процессов управления. 2013. № 4. С. 77–83.

4. Уилер Д., Чамберс Д. Статистическое управление процессами. Оптимизация бизнеса с использованием контрольных карт Шухарта. М.: Альпина Бизнес Букс, 2009. 409 c.

5.     Ryan T.P. Statistical methods for quality improvement. NY, John Wiley and Sons, 2011, 687 р.

6.     Montgomery D.C. Introduction to statistical quality control. NY, John Wiley and Sons, 2009, 754 р.

7.     Боровиков В. STATISTICA: искусство анализа данных. Для профессионалов. СПб: Питер, 2001. 656 с.

8.     Bersimis S., Psarakis S., Panaretos J. Multivariate statistical process control charts: an overview. Quality and Reliability Eng. Intern., 2007, vol. 23, рр. 517–543.

9.     Клячкин В.Н., Кравцов Ю.А. Обнаружение нарушений при многомерном статистическом контроле технологического процесса // Программные продукты и системы. 2016. № 3. С. 192–197.

References

  1. Vladislavlev L.A. Vibration of hydroelectric units in hydroelectric stations. Moscow, Energiya Publ., 1972, 153 p.
  2. Klyachkin V.N., Kuvayskova Yu.E., Aleshina A.A. Modeling of vibrations of a hydrounit on the basis of adaptive dynamic regressions. Automation. Modern Technologies. 2014, no. 1, pp. 30–34 (in Russ.).
  3. Kuvayskova Yu.E., Aleshina A.A. Increasing the efficiency of technical object control system with the use of adaptive dynamic regression modeling of the time series. Automation of Control Processes. 2013, no. 4, pp. 77–83 (in Russ.).
  4. Wheeler D.J., Chambers D.S. Understanding Statistical Process Control. SPC Press, 1992, 406 p. (Russ. ed.: Moscow, AlpinaBiznesBuks Publ., 2009, 409 p.).
  5. Ryan T.P. Statistical Methods for Quality Improvement. NY, John Wiley and Sons Publ., 2011, 687 р.
  6. Montgomery D.C. Introduction to Statistical Quality Control. NY, John Wiley and Sons, 2009, 754 р.
  7. Borovikov V. STATISTICA: the Art of Data Analysis. For Professionals. St. Petersburg, Piter Publ., 2001, 656 p.
  8. Bersimis S., Psarakis S., Panaretos J. Multivariate Statistical Process Control Charts: An Overview. Quality and reliability Engeneering Intern. 2007, vol. 23, рр. 517–543.
  9. Klyachkin V.N., Kravtsov Yu.A. Irregularities in multivariate statistical control of a technological process. Software & Systems. 2016, no. 3, pp. 192–197 (in Russ.).

Permanent link:
http://swsys.ru/index.php?page=article&id=4508&lang=en
Print version
Full issue in PDF (29.03Mb)
The article was published in issue no. № 3, 2018 [ pp. 620-625 ]

Perhaps, you might be interested in the following articles of similar topics: