Как выбрать частоту дискретизации для сбора вибрационных данных: практическое руководство

Частота дискретизации (Fs) — первый параметр, который задаёт предел того, что вы вообще сможете увидеть в спектре вибрации. Поставьте её слишком низко — и дефекты подшипников, резонансы корпуса, зубчатые зацепления просто «свернутся» в низкочастотный шум. Поставьте излишне высоко — и загоните хранилище и каналы связи в терабайты бесполезных отсчётов, не выиграв ни герца диагностической информации.

Главный ориентир: частота дискретизации должна быть не просто «больше частоты интересующего дефекта», а обеспечивать надёжное разделение полезного сигнала и его зеркальных копий (алиасов) после прохождения через антиалиасинг-фильтр. На практике это значит выбор коэффициента перекрытия относительно границы полосы пропускания датчика и фильтра, а не относительно теоретического Найквиста.

Физика процесса: почему теорема Найквиста — это только начало

Теорема Найквиста–Шеннона гласит: для восстановления сигнала без потерь частота дискретизации должна быть как минимум в два раза выше максимальной частоты в сигнале (Fs ≥ 2·Fmax). В учебниках это выглядит как чёткая граница. В реальном датчике и АЦП граница размыта.

Любой реальный антиалиасинг-фильтр (AAF) имеет полосу пропускания, полосу подавления и переходную полосу. Если вы ставите Fs = 2·Fmax, то частота Найквиста (Fs/2) совпадает с Fmax. Но фильтр не может идеально отсечь всё выше Fmax — у него есть спад. Компоненты в переходной полосе попадут в диапазон 0…Fs/2 как ложные низкочастотные артефакты.

Поэтому инженерный запас задаётся коэффициентом k = Fs / (2·Fmax). Типичные значения:

• k = 2.56 — стандарт для многих промышленных анализаторов (например, линейка CSI/Emerson, Bruel & Kjaer). Даёт комфортную переходную полосу для фильтров 8-го порядка.
• k = 3.2 – 4.0 — консервативный выбор для критических машин, когда цена пропущенного дефекта высока.
• k = 2.0 – 2.2 — допустимо только при использовании фильтров очень высокого порядка (12+) и полной уверенности, что выше Fmax сигнала нет ничего.

Правило «Fs = 2.56 × Fmax» — не догма, а результат компромисса между стоимостью АЦП 1990-х и качеством доступных аналоговых фильтров. Современные sigma-delta АЦП с цифровыми децимационными фильтрами позволяют работать ближе к теоретическому пределу, но настройку всё равно делает инженер, а не маркетинг даташита.

Что задаёт Fmax: не «частота вращения», а физика дефектов

Частая ошибка новичков: ставить Fmax = 10× или 20× от частоты вращения (1×RPM). Это работает для дисбаланса и невыравнивания, но проваливается на подшипниках и резонансах.

Реальный Fmax определяется самым высокочастотным признаком, который вы хотите ловить:

Источник вибрации	Типичный частотный диапазон	Что диктует Fmax
Дисбаланс, невыравнивание, ослабление креплений	0.5–10× RPM	10× RPM обычно достаточно
Зубчатые зацепления (GMF)	Z×RPM ± боковые полосы	Число зубьев × RPM + 3–5 боковых полос
Роллационные подшипники (BPFO, BPFI, BSF, FTF)	0.3–5 kHz (низкие гармоники), до 10–20 kHz (ударные компоненты)	3–5-я гармоника BPFI/BPFO для раннего стадия; ударный резонанс Körper для очень ранней стадии
Резонансы корпуса/фундамента	1–10 kHz (часто 2–5 kHz)	Первый–второй резонанс, модулируемый дефектом
Электрические дефекты двигателей (эксцентричность ротора, дефекты обмотки)	2×f_line (100/120 Hz) и гармоники, боковые полосы вокруг 1×RPM	До 1–2 kHz обычно достаточно
Кавитация насосов, турбулентность	Широкополосный шум 1–20 kHz	Зависит от задачи: детектирование факта — 5 kHz; диагностика стадии — 10–20 kHz

Для подшипников ключевой параметр — частота резонанса самого подшипника (bearing natural frequency), обычно 2–8 kHz для катящихся тел и 5–15 kHz для колец. Удар от дефекта возбуждает этот резонанс. Если Fmax ниже резонанса — вы видите только низкочастотную модуляцию (обволачку), теряя информацию об амплитуде ударов. Для раннего обнаружения (стадия 1–2 по ISO 17359) часто нужно Fmax ≥ 5–10 kHz, а для трендинга известного дефекта — достаточно 2–3 kHz.

Практический алгоритм выбора Fs

1. Определите целевую полосу диагностики Fmax. Исходя из типа оборудования, критичности и стадии дефекта, которую нужно ловить. Для общего мониторинга насосов/вентиляторов/компрессоров — часто 2–5 kHz. Для бумажных машин, турбин, высокоскоростных шпинделей — 10–20 kHz.
2. Проверьте полосу пропускания датчика. Пьезоакселерометр с резонансом 25 kHz и полосой ±1 dB до 10 kHz не даст чистый сигнал на 20 kHz. Fmax не может превышать полосу датчика.
3. Узнайте порядок и тип антиалиасинг-фильтра в вашем приборе. Если это аналоговый фильтр 4-го порядка — нужен k ≥ 3.2. Если 8-го — k = 2.56 комфортен. Если sigma-delta с цифровым фильтром — смотрите datasheet: часто там уже заложен k = 2.0–2.2 при указанной Fmax.
4. Рассчитайте Fs = k × 2 × Fmax. Округлите до стандартного значения прибора (1024, 2048, 5120, 10240, 20480, 25600, 51200 Гц).
5. Оцените нагрузку на хранение и передачу. Один канал 16-бит, Fs = 25.6 кГц = 51.2 кБ/с = 4.4 ГБ/сутки. Десять таких каналов — 44 ГБ/сутки. Нужно ли непрерывное записывание сырого сигнала или достаточно спектров/параметров?
6. Проведите тестовый замер. Снимите спектр на выбранной Fs и на удвоенной. Если разницы в диагностически важных пиках нет — можно понизить. Если появляются новые пики выше Fmax/2 — значит, алиасинг или реальный сигнал выше порога.

Специфика разных типов оборудования

Электродвигатели (до 3000 об/мин, 50/60 Гц)

Основные дефекты: подшипники, эксцентричность, нарушение симметрии обмотки. Fmax = 2–5 kHz покрывает 90% задач. Fs = 10.24–25.6 кГц. Для высоковольтных машин с катушечными роторами иногда нужен Fmax до 10 kHz (резонансы статора).

Зубчатые редукторы

GMF = Z × RPM. Боковые полосы от модуляции несущей дефектом зубьев расширяют спектр. Правило: Fmax ≥ GMF + 5×RPM (для модуляции) или ≥ 3×GMF (для анализа формы зуба). При Z=50, 1500 об/мин GMF = 1250 Гц → Fmax ≈ 4–5 kHz. При высокоскоростных редукторах (30 000 об/мин, Z=30) GMF = 15 kHz → Fmax ≥ 30–50 kHz, нужны специализированные датчики и АЦП.

Роллационные подшипники: ловушка «низкочастотного мониторинга»

Частота дефекта внешнего кольца BPFO ≈ 0.4 × n × RPM (n — число тел). При 1500 об/мин, n=8 → BPFO ≈ 80 Гц. Новички ставят Fmax = 1 kHz («ведь дефект на 80 Гц»). Но физика удара: импульс возбуждает резонанс подшипника на 3–8 kHz. Амплитуда этого резонанса растёт раньше, чем амплитуда 1×BPFO. Если Fmax < резонанса — вы обнаружите дефект на 30–50% позже. Компромисс: Fmax = 5 kHz для раннего обнаружения, 2.5 kHz для трендинга известного дефекта.

Турбины, компрессоры, высокоскоростные шпиндели (>10 000 об/мин)

1×RPM уже в сотнях Гц. Подшипники скольжения дают субсинхронные компоненты (0.4–0.48×RPM). Резонансы лопаток/колес — в диапазоне 1–10 kHz. Fmax обычно 10–20 kHz, Fs = 51.2–102.4 кГц. Требуются датчики с резонансом >50 kHz и кабелями с низкой ёмкостью.

Алиасинг: как узнать, что вы его подцепили

Алиасинг (подмена частот) возникает, когда компоненты выше Fs/2 отражаются в полосу 0…Fs/2. Признаки:

• Пики в спектре, которые не соответствуют никаким известным частотам машины (не 1×RPM, не GMF, не BPFO/BPFI, не электрические).
• Пики, которые «перескакивают» при смене Fs: если на Fs=10 кГц пик на 3.2 кГц, а на Fs=20 кГц он уходит на 6.8 кГц — это алиас от реальной частоты 13.2 кГц (Fs — 6.8 = 3.2).
• Шумовой пол, поднимающийся к Fs/2 — признак того, что антиалиасинг-фильтр не справляется или его нет.

Проверка: снимите тот же сигнал на двух разных Fs (например, 10 и 20 кГц). Наложите спектры. Реальные пики останутся на месте, алиасы сдвинутся.

Разрешение спектра vs частота дискретизации: не путайте

Частота дискретизации задаёт диапазон частот (0…Fs/2). Количество отсчётов N (размер БПФ) задаёт разрешение Δf = Fs/N.

Хотите видеть боковые полосы вокруг GMF при шаге 1 Гц? При Fmax=5 kHz (Fs=12.8 кГц) нужно N = 12 800 точек. При Fmax=20 kHz (Fs=51.2 кГц) — уже 51 200 точек. Больше Fs — больше вычислений и памяти для того же разрешения.

Практический совет: если задача — разделение близких пиков (боковые полосы, биения), подбирайте N под нужное Δf, а не гоняйте за Fs. Для подшипников часто достаточно Δf = 1–2 Гц (N = 8192–16384 при Fs=12.8–25.6 кГц).

Современные подходы: адаптивная и мультирейтовая дискретизация

Классическая схема: один Fs для всего. Современные мониторинговые системы (WirelessHART, IO-Link Wireless, проприетарные IIoT-шлюзы) используют гибридные стратегии:

• Низкочастотный непрерывный поток (Fs = 1–2 кГц, Fmax = 500–1000 Гц) — для трендинга RMS, 1×RPM, 2×RPM, температуры. Малый трафик, года хранения.
• Высокочастотные «снимки» по расписанию или триггеру (Fs = 25–100 кГц, длительность 1–10 с) — для спектрального анализа подшипников, резонансов, зацеплений. Передаются по событию или раз в сутки/смену.
• Order tracking (отслеживание порядков) — передискретизация по углу вращения, а не по времени. Устраняет размывание спектра при нестабильной скорости. Требует тахогенератора или точной оценки RPM из вибрации. Эффективно на переменных скоростях (ветряки, приводы с ЧРП).

Такой подход даёт и экономию трафика, и диагностическую глубину. Но требует умного шлюза, который решает, когда запускать высокочастотный захват (по превышению RMS, по календарю, по команде аналитика).

Типичные ошибки и их последствия

Ошибка	Последствие	Как избежать
Fs = 2×Fmax без учёта переходной полосы фильтра	Алиасинг из зоны затухания фильтра, ложные пики, потеря доверия к данным	Использовать k ≥ 2.56 для аналоговых фильтров 8-го порядка; проверять datasheet АЦП
Fmax = 10×RPM для всего подряд	Пропуск резонансов подшипников, GMF, кавитации	Определять Fmax по физике дефектов конкретной машины
Игнорирование полосы датчика	Сигнал выше резонанса датчика искажён, амплитуда некорректна	Fs/2 ≤ 0.8×F_рез_датчика (правило большого пальца для пьезоакселерометров)
Слишком высокий Fs «на всякий случай»	Перегруз каналов связи, затраты на хранение, дольше БПФ, сложнее анализ	Мультирейтовая стратегия: базовый поток + высокочастотные снимки
Отсутствие проверки на алиасинг при вводе в эксплуатацию	Месяцы сбора «мусорных» данных	Тестовый замер на 2×Fs при пуске, сравнение спектров
Путаница между Fs и разрешением Δf	Неразрешенные боковые полосы, невидимая модуляция	Подбирать N (размер БПФ) отдельно от Fs; Δf = Fs/N

Чек-лист перед настройкой системы сбора

• Известен ли тип машины, режим вращения (const/variable), критичность?
• Какие дефекты являются приоритетными для обнаружения (ранняя стадия подшипников, зубья, резонансы)?
• Какова полоса пропускания установленных датчиков (±3 dB, резонанс)?
• Какой антиалиасинг-фильтр стоит в АЦП/модуле (порядок, тип, коэффициент k в datasheet)?
• Какое разрешение спектра (Δf) нужно для разделения близких пиков?
• Каков бюджет трафика и хранилища (непрерывный поток или по событию)?
• Есть ли тахогенератор для order tracking при переменной скорости?
• Проведена ли верификация на алиасинг при пуске (замер на 2×Fs)?

Сценарии выбора: быстрые ориентиры

Сценарий А: Стандартный насос/вентилятор/компрессор до 3600 об/мин, мониторинг состояния подшипников и общей вибрации.

Fmax = 5 кГц (покрывает 3–5 гармоники BPFO/BPFI и резонанс корпуса). Датчик IEPE, полоса 10 кГц. АЦП с фильтром 8-го порядка. k = 2.56. Fs = 2.56 × 2 × 5 = 25.6 кГц. Стандартное значение приборов. Δf = 1 Гц → N = 25 600 (ближайшее 2^n = 32 768). Хороший базовый выбор.

Сценарий Б: Тот же насос, но беспроводной датчик с ограниченным трафиком (LoRaWAN, BLE).

Непрерывный поток: Fs = 1.024 кГц (Fmax = 400 Гц) — только RMS, 1×, 2×RPM, тренд. Раз в сутки/смену: «снимок» Fs = 25.6 кГц, 4 с — полный спектр подшипников. Шлюз решает, когда слать снимок (по превышению RMS или по расписанию).

Сценарий В: Редуктор высокой скорости, GMF = 12 кГц, нужна диагностика зубьев.

Fmax ≥ 3×GMF = 36 кГц (для формы зуба) или ≥ GMF + 5×RPM = 14 кГц (для модуляции). Датчик с резонансом >70 кГц. АЦП sigma-delta, k = 2.0. Fs = 80–100 кГц. Нужен кабель низкой ёмкости, синхронизация по тахо. Order tracking обязателен, если скорость плавает >1%.

Сценарий Г: Турбина/генератор 3000 об/мин, подшипники скольжения, субсинхронная вибрация.

Интересует 0.1–0.5×RPM (5–25 Гц) и резонансы педестала до 2 кГц. Fmax = 2–3 кГц. Fs = 7.68–10.24 кГц. Важно: высокое разрешение на низких частотах (Δf ≤ 0.125 Гц → N ≥ 65 536 при Fs=8192 Гц). Длинные окна БПФ, возможно, нужен zoom-FFT (зум-спектр) вокруг 0.42×RPM.

Что проверить на дате обращения

Параметры АЦП, порядок фильтров, стандартные значения Fs — зафиксированы в документации конкретного прибора. Полосы датчиков — в каталогах производителей (PCB, Wilcoxon, B&K, Metrix, Vibro-Meter, локальные бренды). Нормы по вибрации (ISO 10816-3, ISO 20816, ISO 17359, API 670, GOST ISO 10816) задают зоны тревоги в полосе 10–1000 Гц или 10–5000 Гц, но не диктуют Fs. Выбор Fmax за пределами нормы — инженерное решение, зафиксируйте его в процедуре мониторинга.

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

Главный принцип и следующий шаг

Частота дискретизации — не число «побольше для надёжности», а результат цепочки: физика дефекта → целевая Fmax → полоса датчика → порядок антиалиасинг-фильтра → коэффициент k → Fs → проверка на алиасинг. Пропуск любого звена даёт либо слепоту к дефектам, либо переплату за данные, либо ложные тревоги.

Следующий шаг: возьмите паспорт датчика и datasheet модуля АЦП, запишите их реальные полосы и параметры фильтра, примените алгоритм из раздела «Практический алгоритм выбора Fs» к своим трём самым критичным машинам. Сравните получившиеся Fs с текущими настройками. Если расхождение >30% — перенастройте и сделайте контрольный замер на удвоенной частоте. Это даст больше уверенности, чем любые общие рекомендации.