Настройка NI USRP для контроля целостности бетонных конструкций

Настройка платформы NI USRP (Universal Software Radio Peripheral) для контроля целостности бетонных конструкций представляет собой нетривиальную инженерную задачу, лежащую на стыке радиолокационного зондирования подповерхностных сред и программно-определяемого радио (SDR). Бетонные конструкции — мосты, плотины, туннели, здания — подвержены образованию трещин, пустот, коррозии арматуры и расслоению. Для неразрушающего контроля таких объектов традиционно используются георадары (GPR) и ультразвуковые системы. NI USRP, будучи гибкой SDR-платформой, может быть сконфигурирован как маломощный георадар с возможностью программного управления всеми параметрами сигнала — от несущей частоты до формы импульса. Однако достижение приемлемой глубины зондирования (до 50–100 см в бетоне) и пространственного разрешения (порядка нескольких сантиметров) требует тщательной настройки как аппаратной части (выбор модели USRP, антенн, усилителей), так и цифровой обработки сигналов. В данном ответе рассматриваются шаги по конфигурации NI USRP для контроля бетона, включая выбор частотного плана, синхронизацию приёмо-передающих каналов, использование открытого программного обеспечения (GNU Radio) и алгоритмы визуализации дефектов. Внимание уделяется также практическим ограничениям: сильному затуханию электромагнитных волн в бетоне, помехам от арматурной сетки и необходимости калибровки.

Бетон представляет собой гетерогенную диэлектрическую среду с относительной диэлектрической проницаемостью ε_r от 5 до 9 (в зависимости от влажности и марки) и проводимостью до 0,1 См/м. Электромагнитная волна, распространяясь в бетоне, затухает по экспоненте с коэффициентом затухания α ≈ (σ/2)·√(μ₀/ε₀ε_r), что для частот выше 1 ГГц приводит к глубине зондирования менее 20 см. Поэтому для контроля целостности на глубинах 30–100 см применяют частоты от 100 до 900 МГц. NI USRP, особенно модели серий B200 и N210, работают в диапазоне от 70 МГц до 6 ГГц, что перекрывает требуемый диапазон. Принцип действия: USRP генерирует зондирующий сигнал (например, линейно-частотно-модулированный импульс или гауссов моноцикл) через передающую антенну; отражённый от границ раздела (бетон-воздух, бетон-арматура, бетон-пустота) сигнал принимается второй антенной или той же при работе в импульсном режиме с циркулятором. Анализируя время задержки и амплитуду отражений, можно построить разрез внутренней структуры. USRP позволяет менять форму сигнала программно, что выгодно отличает его от специализированных георадаров с фиксированными параметрами. Однако для бетона критичны две проблемы: сильное отражение от поверхности (до -14 дБ) и рассеяние на щебне. Их решение достигается выбором согласованной антенны и последующей фильтрацией.

Для контроля бетона подходят следующие модели NI USRP:

• USRP B200/B210 (USB 3.0, частота 70 МГц – 6 ГГц, полоса до 56 МГц). B210 имеет два независимых канала приёма-передачи, что удобно для дифференциальных измерений. Однако выходная мощность передатчика ограничена 10–15 дБм (около 30 мВт), что недостаточно для зондирования бетона на глубину более 40 см без внешнего усилителя.
• USRP N210 (Gigabit Ethernet, частота 0–6 ГГц, полоса до 40 МГц, опционально внешний GPSDO для синхронизации). Большая динамика (16 бит АЦП/ЦАП) полезна для приёма слабых отражений. Мощность также мала — до 20 дБм.
• USRP E320 (встроенный ARM, полоса до 100 МГц, два канала, мощность до 20 дБм). Плюс — возможность работы автономно, но цена выше.

Поскольку мощности встроенных усилителей недостаточно, необходим внешний широкополосный усилитель мощности (например, ZX60-33LN-S+ с усилением 20 дБ и выходной мощностью до 30 дБм). При этом следует помнить о допустимой мощности на входе приёмника USRP (максимум -15 дБм для защиты). Поэтому применяют дуплексер или переключатель приёма-передачи, а также встроенный аттенюатор. Антенны — ключевой элемент: для диапазона 200–900 МГц используют дипольные антенны с нагрузочными резисторами (Vivaldi-антенны или спиральные). Для контактного сканирования поверхности бетона подходят экранированные щелевые антенны с ферритовым поглотителем. Рекомендуемая поляризация — линейная, параллельная направлению арматуры для минимизации помех.

Основной режим работы USRP для контроля бетона — стробоскопический (импульсный) либо непрерывный ЛЧМ (FMCW). Для бетона с типичной ε_r=7 скорость волны v = c/√ε_r ≈ 113 млн м/с. Разрешение по глубине Δd = v / (2B), где B — полоса частот. При полосе 200 МГц разрешение составляет около 28 см — этого недостаточно для выявления тонких трещин. Необходима полоса минимум 500 МГц, что даёт разрешение ~11 см. Однако при полосе выше 1 ГГц глубина падает до 10-15 см. Компромисс — ступенчатое зондирование на двух центральных частотах: сначала 300 МГц для поиска крупных пустот, затем 900 МГц для детализации.

Настройка через UHD (Universal Hardware Driver) выполняется следующими параметрами:

• Центральная частота (f_c) — выбирается из условия минимального поглощения. Для бетона оптимальны 400–600 МГц (минимальный коэффициент затухания ~4 дБ/м).
• Частота дискретизации (samp_rate) — должна быть в 2-3 раза выше полосы сигнала. Например, для полосы 400 МГц установите samp_rate=1.25e6 (коэффициент интерполяции).
• Форма волны — гауссов моноцикл с длительностью 1–2 нс (генерируется как производная гауссова импульса). Для FMCW используется пилообразный сигнал с девиацией частоты.
• Усиление передатчика (Tx gain) — выставляется максимальное значение (например, 30 дБ для B200), но с учётом внешнего усилителя снижается, чтобы не перегрузить приёмник. Начальное значение 15 дБ.
• Усиление приёмника (Rx gain) — подбирается экспериментально по уровню шума. Для слабых отражений полезно включить LNA (до 30 дБ).

Важно: излучаемая мощность ограничена регуляторными нормами (макс. 30 дБм для ISM-диапазонов). В России диапазон 433 МГц разрешён без лицензии при мощности до 10 мВт, поэтому для полевых испытаний используйте частоты 2400 МГц (Wi-Fi) — но глубина зондирования на 2,4 ГГц в бетоне менее 5 см. Лучше работать в экранированной лаборатории или получить экспериментальное разрешение.

Основная экосистема для работы с USRP — GNU Radio с пакетом gr-uhd. Альтернативно — LabVIEW с модулем NI-USRP, но GNU Radio предпочтительнее из-за открытости и гибкости цифровой обработки. Настройка включает следующие этапы:

1. Установка драйверов UHD (версия 4.0 и выше) и прошивки для конкретной модели. Проверка через команду uhd_find_devices и uhd_usrp_probe.
2. Создание графа обработки в GNU Radio Companion. Блоки: Signal Source (например, "Chirp" для FMCW) → UHD: USRP Sink. Приёмный тракт: UHD: USRP Source → Complex to Mag → File Sink или QT GUI Time Sink.
3. В параметрах USRP Sink задаются: "Samp Rate", "Center Freq", "Gain", "Antenna" (выберите TX/RX). Для FMCW необходимо выставить "Command" = "lo_offset" для синхронизации гетеродинов.
4. Для режима зондирования с одним USRP (полудуплекс) используют переключение приёма-передачи через GPIO или внешний PIN-диод. Удобнее использовать два USRP (один только на передачу, другой на приём) — тогда синхронизация по 10 МГц и PPS.
5. В GNU Radio также добавляют блоки для формирования зондирующего импульса: например, "Random Source" с коротким гауссовым шумом, пропущенным через "FIR Filter" с гауссовой импульсной характеристикой.

Пример конфигурации UHD для работы на частоте 500 МГц с полосой 200 МГц (импульс длительностью 5 нс):
uhd_usrp_probe --args="type=b200"
uhd_usrp_set_tx_freq 500e6
uhd_usrp_set_rx_freq 500e6
uhd_usrp_set_tx_gain 20
uhd_usrp_set_rx_gain 25
Затем запуск потока через Python API: tx_streamer = usrp.get_tx_stream(), rx_streamer = usrp.get_rx_stream().

Для получения временных задержек отражений с точностью до наносекунд (что соответствует сантиметрам глубины) необходима жесткая синхронизация опорного генератора USRP и момента начала приёма. В одноплатном варианте (один USRP с переключением TX/RX) внутренняя задержка переключения вносит неопределённость. Лучшее решение — два USRP, объединённых внешним опорным сигналом 10 МГц и сигналом PPS (pulse per second) от GPS-дисциплинированного генератора (например, GPSDO от Jackson Labs). Также можно использовать выход триггера с функции "Burst mode" в UHD.

Шаги синхронизации:

• Подключите выход 10 МГц от внешнего генератора к разъёмам "REF IN" обоих USRP. Установите параметр clock_source = external.
• Подайте PPS на вход "PPS IN". В программном коде перед началом передачи дождитесь фронта PPS и запустите буферизированную передачу строго по этому фронту, используя команду uhd::stream_cmd_t::STREAM_MODE_START_CONTINUOUS с временной меткой.
• Для измерения абсолютного времени задержки отражения от поверхности (t0) необходимо зарегистрировать момент излучения. Это делается путём подачи копии передаваемого сигнала на дополнительный вход приёмника через калиброванную линию задержки.
• В GNU Radio реализована система "Tag Propagation": метки времени добавляются к потоку через блок "UHD: Tagged Stream". Используйте блок "UHD: USRP Source" с опцией "Enable Timestamps".

Без надёжной синхронизации разница фаз между передачей и приёмом будет случайной, что сделает невозможным когерентное накопление сигналов для повышения отношения сигнал/шум (SNR). Для бетона с большим затуханием когерентное накопление необходимо (усреднение 64–256 измерений).

Сырой сигнал с выхода USRP содержит прямую волну (передатчик-приёмник по воздуху), отражение от поверхности бетона, многолучевые помехи и шум. Полезная информация — отражения от арматуры, трещин и пустот — имеет амплитуду на 40–60 дБ ниже поверхности. Типовой пайплайн обработки в GNU Radio или Python (после сохранения в файл):

1. Вычитание фонового сигнала: запись профиля без объекта (на открытом воздухе) и последующее вычитание из каждого A-скана.
2. Фильтрация нижних частот: полосовой фильтр Баттерворта 5-го порядка с полосой, соответствующей зондирующему сигналу.
3. Огибающая (преобразование Гильберта): для импульсных сигналов это даёт аналитический сигнал, амплитуда которого огибающая показывает мощность отражения.
4. Компенсация затухания: умножение огибающей на экспоненциальный множитель exp(2αd), где d — глубина, α — коэффициент затухания (оценивается калибровкой).
5. Миграция (алгоритм Кирхгофа или f-k миграция) — перенос отражённых сигналов в истинные положения подповерхностных объектов, особенно важна для наклонных трещин.
6. Построение B-скана: матрица A-сканов (время-амплитуда) по горизонтальной координате. Для визуализации используют псевдоцвета или серую шкалу.

Пример программного кода (Python с библиотекой numpy и scipy) для одного A-скана:

import numpy as np from scipy.signal import hilbert def process_ascan(signal, fs, f0, attenuation_coeff): filtered = bandpass_filter(signal, f0-100e6, f0+100e6, fs) envelope = np.abs(hilbert(filtered)) t = np.arange(len(envelope))/fs depth = t * (3e8/2/np.sqrt(7)) # для εr=7 comp = envelope * np.exp(2*attenuation_coeff*depth) return depth, comp

Для улучшения SNR применяют "stacking" (накопление нескольких измерений в одной точке). NI usrp позволяет делать до 2000 A-сканов в секунду при частоте дискретизации 1 МГц.

Перед началом контроля бетонной конструкции необходима калибровка на эталонном блоке с известными дефектами (например, просверленные отверстия диаметром 10 мм на глубинах 5, 10, 20 см). Алгоритм калибровки:

• Измерить время прихода отражения от поверхности (t_surface) — это ноль глубины.
• По времени отражения от металлической пластины на известной глубине d_plate вычислить эффективную скорость v = 2d_plate / t_plate.
• Уточнить коэффициент затухания α, сравнивая амплитуду от пластины с теоретическим законом A(d) = A0 * exp(-2αd).
• Настроить усиление приёмника так, чтобы отражение от поверхности не насыщало АЦП (обычно уровень -10 дБ от полной шкалы).

Помехи от арматурной сетки проявляются как регулярные гиперболы на B-скане. Для их подавления применяют:

• Фильтрацию в волновой области (f-k фильтр), удаляющую наклонные компоненты, соответствующие скорости распространения в бетоне.
• Вычитание синтетической модели: зная шаг арматуры (например, 20 см), генерируют синтетический радарограмм и вычитают его.
• Использование двух поляризаций (HH и VV) — трещины проявляются по-разному.

Также необходимо экранирование антенн от внешних радиочастотных помех (сотовые вышки, Wi-Fi). Для этого антенны помещают в алюминиевый корпус с ферритовой обкладкой.

Ниже приведены рекомендуемые параметры USRP для контроля бетонных конструкций с различной глубиной и требуемым разрешением. Значения получены эмпирически для сухого бетона марки M300 (εr=6,5, затухание 3 дБ/м на 400 МГц).

Примечание: разрешение по глубине рассчитано как v/(2B) с v=1.13e8 м/с. Реальное разрешение может быть хуже из-за шумов. Мощность указана на выходе антенны с учётом потерь в кабелях (≈2 дБ). Для глубины более 100 см требуется использовать частоты ниже 200 МГц, но для NI USRP B200 минимальная частота 70 МГц — подходит. Однако антенны для таких частот становятся громоздкими (длина диполя >1 м).

При настройке NI USRP для контроля целостности бетона инженер сталкивается с несколькими фундаментальными ограничениями:

• Динамический диапазон: Отражение от поверхности (через воздух) может быть на 60 дБ сильнее, чем от глубокой пустоты. 12-битные АЦП USRP (кроме N210 с 14 битами) дают динамику ~72 дБ, чего хватает с запасом при правильном подборе усиления. Рекомендуется использовать цифровое подавление прямой волны (адаптивный фильтр).
• Джиттер синхронизации: даже с GPSDO остаточная нестабильность фазы составляет 1-2 нс, что соответствует ошибке глубины 2-3 см. Для прецизионных работ нужен внешний синтезатор с фазовой автоподстройкой (например, ADF4351).
• Влияние влажности бетона: при повышении влажности с 2% до 10% затухание растёт в 5 раз. Необходимо перед сканированием измерить влажность и скорректировать коэффициент α в программной обработке.
• Ограничения по скорости сканирования: для получения B-скана длиной 2 м с шагом 2 см потребуется 100 A-сканов. При накоплении 64 измерений на точку и частоте повторения 500 Гц общее время составит около 13 секунд. Это допустимо для статического контроля, но не для динамического. Повысить скорость можно снижением накоплений до 16, жертвуя SNR.
• Электромагнитная совместимость: мощные импульсы USRP могут создавать помехи для другого измерительного оборудования. Рекомендуется работать в режиме "burst" с малым коэффициентом заполнения (скважность 1:100).

Для преодоления этих ограничений часто используют гибридную схему: NI USRP генерирует зондирующий сигнал, а оцифровка отражений производится внешним высокоскоростным АЦП (например, 14 бит, 1 ГГц), синхронизированным с USRP по PPS. Такой подход увеличивает стоимость, но позволяет достичь глубины до 2 м в бетоне. В целом, правильно настроенная система на базе USRP способна выявлять трещины шириной от 1 мм на глубине до 30 см, пустоты диаметром от 5 см на глубине до 70 см и отслоения арматуры. Это делает её полезным инструментом для неразрушающего контроля в строительной диагностике.