← Ко всем статьям
Альтернативный алгоритм статистической обработки данных УЗК свай — метод межквартильного размаха
УЗК свайМетод IQRГОСТ Р 71039-2023МИТ УЗК 1-2К/1-3К/1-4КДефектоскопия свай

Альтернативный алгоритм статистической обработки данных при УЗК бетона свай

Метод межквартильного размаха (IQR) для ультразвукового контроля сплошности бетона свай — универсальный алгоритм для приборов линейки МИТ УЗК (1-2К, 1-3К, 1-4К) и их аналогов. Совместим с ГОСТ Р 71039-2023.

Денис Воробьев
Денис Воробьев
Начальник лаборатории ООО «МАТТЕСТ»
26 июл. 2026 г.
12 мин чтения
Кратко. Альтернативный алгоритм статистической обработки данных при ультразвуковом контроле сплошности бетона свай построен на методе межквартильного размаха (IQR) и применим для любых приборов УЗК свай, включая линейку МИТ УЗК (1-2К, 1-3К, 1-4К) и её аналоги. Для каждой акустической линии фиксируются глубина, время первого вступления, скорость и затухание. Массив скоростей Vi делится квартилями C1 и C2 на внутренние (C ± 1.5·ΔC) и наружные (C ± 2·ΔC) барьеры: значения за наружными считаются аномальными, за внутренними — выбросами. Средняя (нормальная) скорость Vср принимается равной медиане очищенного массива. Отклонение ΔVi и затухание γij оцениваются по критериям ГОСТ Р 71039-2023: аномалия считается значительной при снижении скорости >25% или 15–25% при затухании >12 дБ. Алгоритм рекомендован для любых испытательных лабораторий вне зависимости от применяемого оборудования.

О чём эта статья

Испытательные строительные лаборатории, а также экспертные центры сталкиваются с проблемой статистической обработки данных при ультразвуковом контроле бетона свай. МИТ УЗК от М-Инструмент предоставляет решение этой проблемы путём применения линейки приборов МИТ УЗК-1-2К, МИТ УЗК-1-3К и МИТ УЗК-1-4К (см. Рис. 1) совместно с программным комплексом МИТ УЗК.

В нашей первой статье «Алгоритм статистической обработки данных УЗК свай в МИТ УЗК МОНОЛИТ» подробно изложен алгоритм, реализованный в программном комплексе МИТ УЗК на базе JGJ 106-2014 (Китай). Здесь же мы постарались изложить общий альтернативный алгоритм, применимый для всех приборов по ультразвуковому контролю свай — как приборов МИТ УЗК, так и их аналогов.

Испытатели, не работающие с приборами линейки МИТ УЗК, сталкиваются с проблемой отсутствия понятного и достаточно простого алгоритма статистической обработки своих результатов измерений. Данная статья — ответ на этот запрос.

Линейка приборов МИТ УЗК: 1-2К двухканальный, 1-3К трёхканальный, 1-4К четырёхканальный
Рис. 1. Линейка приборов МИТ УЗК: 1-2К (двухканальный), 1-3К (трёхканальный), 1-4К (четырёхканальный)

Альтернативный алгоритм: параметры измерения

Данный алгоритм обработки и оценки полученных результатов применим как для данных, полученных приборами МИТ УЗК, так и другими специализированными приборами по ультразвуковому контролю свай.

Для каждой акустической линии i (поперечного сечения профиля j) определяются и фиксируются следующие параметры:

Схема расположения преобразователей при параллельном сканировании j-го профиля УЗК свай
Рис. 2. Схема расположения преобразователей при параллельном сканировании j-го профиля

i — номер акустической линии, который присваивается последовательно снизу вверх (или сверху вниз) для каждого измерительного профиля j.

Шаг 1. Время первого вступления

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

Уровень первого вступления устанавливается обычно в диапазоне от 20% до 80% (см. Рис. 3).

Сигнал МИТ УЗК: ширина диапазона измерения уровня первого вступления 20%
Рис. 3а. Ширина диапазона измерения 20%
Сигнал МИТ УЗК: ширина диапазона измерения уровня первого вступления 80%
Рис. 3б. Ширина диапазона измерения 80%. Регулировка ширины диапазона измерения уровня первого вступления в МИТ УЗК на развёртке сигнала.

Шаг 2. Скорость распространения ультразвукового сигнала

Скорость Vi, м/с, распространения ультразвукового сигнала для i-ой акустической линии (i-ого поперечного сечения) j-го профиля рассчитывается прибором ультразвукового контроля свай автоматически по формуле (1):

Vi(j) = Ljti · 1000 (1)

где:

Панель информации о профиле сваи в МИТ УЗК: корректировка значений расстояния между трубами доступа
Рис. 4. Корректировка значений расстояния между трубами доступа в МИТ УЗК на панели информации о профиле сваи

Шаг 3. Средняя (нормальная) скорость — метод межквартильного размаха

С целью выявления резко выделяющихся (аномальных, выпадающих) значений из общего массива скоростей ультразвуковых волн для j-го профиля используется метод расчёта диаграммы размаха (box plot / IQR).

3.1. Расчёт первого C1 и третьего C2 квартилей. Первый и третий квартили равны соответственно 25-му (1/4) и 75-му (3/4) процентилю.

3.2. Межквартильный (интерквартильный) размах ΔCj:

ΔCj = C2C1 (2)

3.3. Расчёт шага hj:

hj = 1,5 · ΔCj (3)

3.4. Внутренние барьеры β1 и β2:

β1 = C1hj (4)
β2 = C2 + hj (5)

3.5. Наружные барьеры B1 и B2:

B1 = C1 − 2hj (6)
B2 = C2 + 2hj (7)

Значения скоростей Vi i-ой акустической линии j-го измерительного профиля, выходящих за наружные барьеры, принимаются как аномальные значения; значения скоростей Vi, выходящих за внутренние барьеры, считаются выбросами.

Аномальные значения и выбросы значительно влияют на результаты статистического расчёта, следовательно, должны быть отброшены из дальнейшего расчёта. Наиболее робастной оценкой закона распределения является медиана Me [1].

3.6. Расчёт медианы Me очищенного массива скоростей.

Нечётное количество значений скоростей Vi:

Me = V(n+1)/2 (8)

Чётное количество значений скоростей Vi:

Me = (Vn/2 + Vn/2 + 1) / 2 (9)

Полученное значение медианы будет соответствовать средней (нормальной) скорости Vср на j-ом профиле измерений, т. е.:

Vср = Me (10)

В случае, если исходный массив значений скоростей не имеет выбросов и аномальных значений, то в качестве средней (нормальной) скорости профиля принимаем медиану всего массива значений.

Шаг 4. Отклонение скоростей от среднего значения

Рассчитываются отклонения ΔVi всех значений скорости Vi на j-ом профиле измерений от значения средней скорости Vср в % по формуле:

ΔVi = VсрVi · 100 (11)

где:

Шаг 5. Затухание сигнала

В приборах линейки МИТ УЗК предусмотрена функция автоматического определения значения амплитуды сигнала Ai i-ой акустической линии (см. Рис. 5).

Развёртка сигнала и массив данных j-го профиля с автоматически определёнными значениями времени и амплитуды в программном комплексе МИТ УЗК
Рис. 5. Развёртка сигнала и массив данных j-го профиля с автоматически определёнными значениями времени ti и амплитуды Ai программным комплексом МИТ УЗК

Значение затухания γij для каждой i-ой акустической линии j-го профиля рассчитывается по формуле:

γi(j) = −20 · lg (Ai(j) / A0) (12)

где:

1
Примечание 1. Оценка затухания для приборов, отличных от МИТ УЗК.

В случае, если применяются приборы, отличные от МИТ УЗК, то у испытателей иногда возникают некоторые сложности в определении значения амплитуды сигналов. Если в таком приборе предусмотрен режим дефектоскопа — т. е. графический режим отображения результатов измерения для каждой акустической линии (поперечного сечения), то значение амплитуды полученного сигнала оценивается:

  • по величине амплитуды сигнала относительно высоты экрана в %;
  • по величине фронта первого вступления Фi (В/мкс), учитывающему напряжение и время первого вступления. Величина фронта первого вступления пропорциональна величине амплитуды первого сигнала.

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

γi(j) = −20 · lg (Кi(j) / К0) (13)

где:

  • Фij — значение фронта первого вступления i-ого поперечного сечения, измеряемое в процентах высоты экрана прибора;
  • Ф0 — максимальное значение фронта первого вступления, равное максимальному значению, зарегистрированному для j-го профиля.
2
Примечание 2. Ограничения приборов без измерения энергии сигнала.

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

Методика анализа полученных данных

Критерии выделения ультразвуковых аномалий по отклонению значений скорости распространения ультразвуковых волн Vi от средней скорости Vср и по значениям затухания ультразвуковых волн αj для каждого j-го профиля соответствуют требованиям ГОСТ Р 71039-2023 и приведены на Рис. 6.

Критерии выделения ультразвуковых аномалий по ГОСТ Р 71039-2023: снижение скорости и значение затухания
Рис. 6. Критерии выделения ультразвуковых аномалий по ГОСТ Р 71039-2023

Значительная аномалия интерпретируется как нарушение сплошности бетона конструкции.

Аномалия считается значительной, если определено:

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

Если при выполнении измерений выявлены участки профиля со значительным снижением скорости (более 15% от среднего значения), то шаг сканирования уменьшается до 50 мм. При необходимости такие участки могут быть проверены повторно.

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

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

Заключение: совместимость с ГОСТ Р 71039-2023

Приведённый выше алгоритм статистической обработки и методика анализа полученных данных в целом соответствуют требованиям ГОСТ Р 71039-2023 и могут быть рекомендованы нами к применению испытательными лабораториями вне зависимости от применяемого оборудования — линейки МИТ УЗК или её аналогов.

Видеообзор прибора и программного комплекса

Работа с прибором МИТ УЗК-1-2К:

Работа с прибором МИТ УЗК-1-2К — превью видео Смотреть на YouTube

Подробнее о работе с программным комплексом МИТ УЗК:

Работа с программным комплексом МИТ УЗК — превью видео Смотреть на YouTube

Первая статья Дениса про алгоритм JGJ 106-2014 в МИТ УЗК МОНОЛИТ доступна здесь: Алгоритм статистической обработки данных УЗК свай в МИТ УЗК МОНОЛИТ.

Заказать линейку МИТ УЗК и программный комплекс

Линейка МИТ УЗК: 1-2К, 1-3К, 1-4К + программный комплекс. Автоматический расчёт времени первого вступления, скорости и затухания. Совместимость с ГОСТ Р 71039-2023. Приборы 2/3/4 канала для параллельного и наклонного сканирования.

Узнать цену и заказать →

Ответы на частые вопросы

Чем альтернативный алгоритм отличается от алгоритма в МИТ УЗК МОНОЛИТ?

Алгоритм МИТ УЗК МОНОЛИТ построен на JGJ 106-2014 (Китай): итерационное отбрасывание экстремумов + λ-коэффициент. Альтернативный алгоритм построен на классическом методе межквартильного размаха (IQR) и применим для любых приборов УЗК свай, а не только МИТ.

Для каких приборов подходит этот алгоритм?

Для любых специализированных приборов ультразвукового контроля сплошности бетона свай, включая линейку МИТ УЗК (1-2К, 1-3К, 1-4К) и их аналоги, при условии, что прибор фиксирует и время прохождения, и энергию/амплитуду сигнала.

Что делать, если прибор не измеряет амплитуду сигнала?

Такой прибор не может быть применим для полноценного контроля сплошности бетона свай. Если в приборе есть режим дефектоскопа (графический режим), то амплитуда может быть оценена по величине сигнала относительно высоты экрана в % или по фронту первого вступления Фi (В/мкс).

Как определяется значительная аномалия по ГОСТ Р 71039-2023?

Аномалия считается значительной, если снижение скорости более 25% при любых значениях затухания, либо снижение скорости 15–25% при затухании более 12 дБ. Если аномалия на одной глубине встречается как незначительная в двух и более продольных сечениях сваи — она также считается значительной.

Что делать при выявлении участка со снижением скорости более 15%?

Шаг сканирования уменьшается до 50 мм, участок при необходимости проверяется повторно. Анализ ведётся по каждому продольному профилю отдельно, вывод о сплошности сваи делается с учётом всех профилей.

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

  1. Хампель Ф., Рончетти Э., Рауссеу П., Штаэль В. Робастность в статистике. Подход на основании функций влияния. — М.: Мир, 1989.
  2. ГОСТ Р 71039-2023. Сваи буронабивные и «стены в грунте» траншейного и свайного типов. Межскважинный ультразвуковой метод контроля качества бетона.

Автор статьи: Денис Воробьев, начальник лаборатории ООО «МАТТЕСТ». Экспертный текст подготовлен для сайта партнёрской компании ООО «М-Инструмент» (МАТТЕСТ ИНСТРУМЕНТ).