Свидетельство:
О регистрации средства массовой информации: "Предотвращение аварий зданий и сооружений".
Номер: №ФС77-35253
Выдано: Федеральная служба по надзору в сфере связи и массовых коммуникаций
Дата: от 16.02.2009 г.
Форма распространения: электронное периодическое издание
Язык: русский
Учредитель: ООО "ВЕЛД"
Инженерные сооружения по конструкции являются замкнутым объёмом, и распространение сейсмических колебаний в них приводит к формированию стоячих волн. Ещё в девятнадцатом веке была доказана общеизвестная теорема [6, 8, 9] о независимости решения по стоячим волнам от начальных условий. Для инженерных сооружений это означает, что стоячие волны, формирующиеся в них, не зависят от источника колебаний. Для поля стоячих волн в здании не важно, где размещён источник колебаний и какого типа источник используется.
Опираясь на данную теорему, можно сказать, что детальность и точность исследования стоячих волн в инженерных сооружениях должны быть одинаковы независимо от того, изучаем мы микросейсмические колебания или применяем высокостабильные по излучению искусственные источники колебаний. До последнего времени в экспериментах по изучению зданий и сооружений измерения на микросейсмических колебаниях обосновано считались менее точными, чем исследования с искусственными источниками. Дело в том, что принципиальные возможности метода и практическая реализация не всегда совпадают.
В данной работе обосновывается методика наблюдений за микросейсмами в инженерном сооружении, создаётся математическая модель волнового поля, строятся алгоритмы, позволяющие выделить стоячие волны в чистом виде и преобразовывать их из разновремённых наблюдений в одновременные записи стоячих волн на плотной системе наблюдений в инженерном сооружении. Рассматриваются возможности качественной интерпретации детальной информации о стоячих волнах для диагностики физического состояния и определения сейсмостойкости сооружений, а также затрагиваются проблемы разработки методов количественной интерпретации детальных наблюдений.
Колебания для каждой из нормальных мод описываются выражением [х]
, (1)
где – геометрическая форма
стоячей волны;
– фаза колебаний в
заданной точке для конкретной нормальной моды;
– частота собственных
колебаний.
Геометрическая форма стоячей волны, фаза и частота собственных колебаний зависят от строения инженерного объекта, от физических параметров материала, из которого изготовлено сооружение. При детальном измерении в объёме инженерного сооружения колебаний на каждой из нормальных мод мы можем изучить строение объекта и выявить возникающие в нём дефекты и ослабленные места. Такие исследования можно отнести к сейсмической дефектоскопии инженерных сооружений.
Стандартная оценка сейсмостойкости зданий основана на определении усиления колебаний в здании на частотах наиболее интенсивных нормальных мод [7]. Детальное изучение стоячих волн в инженерном сооружении – это не только определение дефектов в здании, но и новый уровень точности и детальности в оценке сейсмостойкости.
Из уже упомянутой теоремы вытекает важное свойство когерентности во времени колебаний в стоячих волнах. По стоячим волнам линейная связь между записями колебаний двух точек не зависит от времени. И эта связь выполняется как для амплитуд, так и для фаз стоячих волн. Отмеченное свойство когерентности во времени мы будем использовать для создания методики детального изучения характеристик стоячих волн в инженерных сооружениях.
Для детального изучения стоячих волн в инженерном сооружении необходимо решить два актуальных вопроса. Если мы установим сейсмоприёмник в здании, то будем регистрировать сумму стоячих и бегущих волн. Первый вопрос: как выделить из регистрируемых колебаний стоячие волны? Детальное изучение стоячих волн в инженерном сооружении требует плотных систем наблюдений. При обследовании крупных сооружений потребуются сотни и даже тысячи точек с одновременными записями стоячих волн. Второй вопрос: как обойтись для изучения стоячих волн в инженерном сооружении на плотной системе наблюдений, измерениями с малоканальной аппаратурой?
Основу решения поставленных задач составит метод восстановления когерентных составляющих сейсмических волновых полей [3], позволяющий выделять когерентные волны из экспериментально наблюдаемых волновых полей. В упомянутой работе введено понятие когерентность по параметру. Дело в том, что, когда мы говорим о линейной связи между сейсмическими записями, то независимой она может быть от любого параметра, например, по измерению в пространстве, по элементам системы наблюдений, по площади, по энергии и т.п. Когерентность по времени – это частный случай когерентных сигналов.
Для детального обследования здания с использованием микросейсм идеально реализовать плотную систему одновременной регистрации колебаний, но это не представляется возможным. Рассмотрим иную систему наблюдений и её возможности. Одновременная регистрация колебаний здания под воздействием микросейсм ведётся в опорной и i-й точках (группе точек), затем i-я точка (группа точек) меняет свое положение и вновь проводится регистрация сейсмических колебаний одновременно с опорной точкой. Такими наблюдениями можно детально покрыть исследуемый объект с малоканальной аппаратурой.
Задача состоит в том, как преобразовать разновремённые наблюдения в разных точках здания в одновременную запись стоячих волн на всей системе наблюдений. Для решения этой задачи необходима математическая модель колебаний. Представим колебания каждой точки как сумму двух составляющих. Первая составляющая – стоячие волны, вторая – бегущие волны.
Заложим в основу модели сейсмических колебаний инженерных сооружений, следующие предположения:
1. При воздействии микросейсмических колебаний на здание (инженерное сооружение) отличия формирующихся в нём стоячих волн, зарегистрированных в двух произвольных точках, описываются линейной системой, характеристика которой не зависит от времени.
2. При воздействии микросейсмических колебаний на здание отличия бегущих волн, зарегистрированных в двух произвольных точках, не могут быть описаны линейной системой, характеристика которой не зависит от времени.
3. Линейные связи в изменениях стоячих волн, регистрируемых в двух точках, существуют для каждой пары одноимённых компонент регистрируемых колебаний.
Третье предположение существенно упрощает поставленную задачу. Вместо того, чтобы работать с векторными сигналами с девятью характеристиками связи, мы будем работать с данными регистрации по каждой из компонент колебаний независимо друг от друга и обойдёмся тремя характеристиками связи.
Модель связи колебаний, регистрируемых в двух разных точках инженерного сооружения (здания), примет вид
,
,
(2)
где – импульсная характеристика
линейной системы, описывающей связь между одновременными записями стоячих
волн в точках 0 и i на обследуемом объекте;
– бегущие волны;
– запись стоячей волны
в опорной точке.
Подход к обработке инженерно-сейсмологических данных рассматриваемый в данной работе является развитием работ [2-5].
Для модели процесса колебаний в здании есть возможность получить одновременные записи стоячих волн из разновремённых, последовательных наблюдений с опорной точкой. Процедура обработки в таком случае сводится к следующим операциям.
1. Нахождение частотных характеристик линейных
систем .
2. Запись или формирование независимой реализации процесса колебаний опорной точки при сейсмическом воздействии на исследуемый объект.
3. Пересчёт стоячих волн из опорной точки с использованием
во все точки обследуемого
объекта.
Ключевым вопросом рассматриваемой схемы является
задача определения с необходимой точностью
по записям микросейсмических колебаний. Рассматривая модель сигналов, регистрируемых
в двух точках здания на одинаковых компонентах, можно видеть, что задача определения
сводится к определению
характеристики линейной системы по сигналам на её входе и выходе, зарегистрированных
на фоне шумов. Фактически необходимо определить характеристику фильтра, обеспечивающего
наилучший пересчёт стоячих волн из опорной точки в i-ю.
Будем искать характеристику
в виде оптимального
фильтра Винера [1], преобразующего сигнал
в запись стоячей волны
в точке i. Для этого необходимо минимизировать по
математическое ожидание
квадрата ошибки между сигналом
, пропущенным через
фильтр, и сигналом
, являющимся выходом
линейной системы [1]:
(3)
Взяв производную по от данного выражения
и приравняв её к нулю, получим
Учитывая, что – взаимная корреляционная
функция между сигналами
и
, а
– автокорреляционная
функция колебаний, зарегистрированных в опорной точке, получим
.
(4)
Осуществив преобразование Фурье и перейдя в частотную область, получим
(5)
В знаменателе выражения (5) стоит усреднённый
квадрат модуля спектра колебаний в опорной точке. Решив, каким образом осуществляется
усреднение, можно осуществить его вычисление. В числителе выражения (5) стоит
усреднённый взаимный спектр колебаний в опорной точке и стоячей волны в i-й
точке. экспериментально не
измеряется. Приняв, что
и
являются стационарными,
случайными процессами, некоррелированными друг с другом и с записями стоячих
волн в опорной и i-й точках, легко доказывается, что
[1]
(6)
Выражение (6) позволяет подставить в числитель формулы (5) усреднённый взаимный спектр колебаний, зарегистрированных одновременно в опорной и i-й точках.
Усреднение спектров в данной ситуации можно осуществить, разбивая реализацию колебаний, записанных одновременно, на n блоков. Разбив на неперекрывающиеся блоки запись, получаем множество реализаций, по которым можно осуществить усреднение. Формула для расчёта фильтра Винера, обеспечивающего пересчёт колебаний из опорной точки в i-ю, примет вид
. (7)
Отличия частотной характеристики фильтра Винера от характеристики линейной системы, описывающей различия стоячих волн в двух точках обследуемого объекта, описываются смещённостью оценки, которая определяется отклонением математического ожидания определяемой частотной характеристики от истинного её значения. Можно показать, что
. (8)
В соответствии с (8) частотная характеристика фильтра является смещённой оценкой характеристики линейного фильтра, описывающего связь колебаний двух точек обследуемого объекта. Смещённость оценки выражается в регуляризации фильтра по соотношению шум/сигнал. Другими словами, фильтр запирается на частотах с низким соотношением сигнал/шум. Отметим, что оценка фазовой характеристики фильтра не смещена.
Формула (7) позволяет вести расчёт характеристики
фильтра с погрешностью, которая зависит от следующих параметров наблюдений:
шаг дискретизации по времени – , длина единичного
блока – Т, количество блоков в записи – n. Выбор первых двух параметров не вызывает затруднений. Шаг
дискретизации увязывается с частотным диапазоном, в котором изучается реакция
объекта на сейсмическое воздействие. Длина единичного блока связана с необходимым
разрешением спектрального анализа
. Сложнее вопрос о
количестве блоков. Расчёт по формуле (7) – лишь некоторая оценка характеристики
фильтра, погрешность которой зависит от количества блоков n
и соотношения энергии шумов и полезных сигналов в модели (2).
Для оценки количества блоков, необходимых для
обеспечения заданной погрешности фильтра, полезен спектр когерентности
, который можно рассчитать
по записям микросейсмических колебаний в двух точках обследуемого объекта
по формуле
,
(9)
где – спектр автокорреляционной
функции колебаний, записанных в опорной точке;
– тоже в i-й точке.
Применяя модель (2), получим:
,
,
.
Эти выражения получены в предположении о независимости реализаций шумов на входе и выходе линейной системы и их некоррелированности с полезными сигналами. Подставляя полученные соотношения в (9), получим
,
где ,
– отношения средневзвешенных
квадратов шум/полезный сигнал в опорной и i-й точках.
Из формулы (10) следует, что спектр когерентности
зависит только от отношения энергии бегущих волн к энергии стоячих волн. Спектр
когерентности может быть рассчитан по одновременным записям в любой паре точек
обследуемого объекта. Значения лежат в диапазоне
0-1. Единица соответствует случаю, когда в двух точках регистрируются только
стоячие волны. При
в двух точках здания
регистрируются только бегущие волны.
Экспериментами доказано [2-5], что значения спектра когерентности на собственных частотах инженерных сооружений высоки 0.8-0.99, а в промежутке между ними соответствуют значениям 0.1-0.3.
Полученный алгоритм расчёта характеристики фильтров является симбиозом фильтра Винера, свойства когерентности стоячих волн и методов математической статистики. Математическая статистика позволяет не только строить фильтры для выделения когерентных во времени стоячих волн, но и определять точность полученного волнового поля через дисперсию оценки характеристик. Подробное описание алгоритмов оценки точности можно найти в работе [4]. Среднеквадратическая ошибка фазовой характеристики фильтра Винера для пересчёта стоячих волн из опорной точки в i-ю будет равна
.
Данная формула позволяет рассчитать ошибку фазовой характеристики в зависимости от частоты в радианах. Относительная ошибка амплитудной характеристики фильтра равна
,
где – среднеквадратическое
отклонение оценки частотной характеристики фильтра.
Из этих формул заданная точность и число блоков могут быть увязаны после вычисления значения спектра когерентности на частоте собственных колебаний здания. Количество блоков определяет длину реализации, записанную в каждой точке обследуемого объекта, и существенно влияет на производительность работ. Оптимальная производительность достигается последовательной регистрацией колебаний в точках объекта при длине реализации 5-10 минут и существенно снижается при длинах реализаций, измеряемых часами. При высокой когерентности колебаний 0.8-0.99, наблюдаемой на инженерных сооружениях на частотах нормальных мод, удаётся достигнуть погрешности порядка 5% и менее для пересчёта стоячих волн и производительности обследования зданий – за несколько дней, а такого объекта, как плотина Саяно-Шушенской ГЭС, – за две недели (при детальности обследования через 5-10 м).
Пример восстановления поля стоячих волн дан на рис.1. На рисунке приведены разновременные наблюдения в одной из галерей Саяно-Шушенской ГЭС (рис.1,а), трассы пересчитанных в эти же точки системы наблюдений колебаний (рис.1,б), результат двухмерной фильтрации пересчитанных трасс (рис.1,в). Вместо нерегулярных колебаний мы получили регулярное волновое поле, в котором явно видны отражения от бортов плотины.
Рис.1. Пример пересчёта стоячих волн из опорной точки
в другие
точки системы наблюдений:
а – разновремённые записи микросейсмических колебаний в галерее
плотины Саяно-Шушенской ГЭС (отметка 521 м);
б – результат пересчёта стоячих волн из опорной точки в фрагмент
системы наблюдений, соответствующий а;
в – двумерная фильтрация б
По этим записям начинается детальное обследование. Первый вопрос обследования – это обнаружение и изучение нормальных мод колебаний инженерного сооружения. На рис.2 дан монтаж спектров трасс вдоль одной из галерей плотины Саяно-Шушенской ГЭС (радиальные колебания). По спектрам модуляция амплитуды от трассы к трассе позволяет сразу же определять кратность обнаруженной моды. Для таких объектов, как плотины ГЭС, предложенная методика позволяет выявлять и изучать особенности распределения амплитуд и фаз стоячих волн по пространству на нескольких десятках нормальных мод.
Рис.2. Карты амплитудных спектров радиальной скорости
смещения
пересчитанных колебаний в одну из галерей Саяно-Шушенской ГЭС
(отметка 521 м). Длина секции 15 м:
а – для уровня воды в водохранилище 539 м;
б – для уровня воды в водохранилище 500 м
Выделенные стоячие волны легко отфильтровываются и строятся карты амплитуд и фаз колебаний по каждой из нормальных мод в обследуемом объекте. По ним изучаются особенности каждой из стоячих волн. При обследовании зданий основная информация содержится в искажениях геометрической формы стоячей волны, вызванных физическим состоянием конструкций. Анализ искажений позволяет проводить диагностику физического состояния зданий и выявлять элементы наиболее вероятных разрушений при сейсмических воздействиях.
Заключение
1. Разработаны алгоритмы восстановления одновременных записей стоячих волн в инженерных сооружениях по разновремённым измерениям микросейсмических колебаний с записями опорных колебаний в одной или нескольких точках, а также алгоритм оценки точности, в основе которых заложены методы Винеровской фильтрации, свойство когерентности стоячих волн во времени и методы статистического оценивания характеристик.
2. Экспериментально доказаны высокая точность и детальность восстановления стоячих волн в инженерных сооружениях как по пространственным измерениям, так и по шкале частот колебаний.
3. Созданная методика детального изучения стоячих волн в инженерных сооружениях является основой диагностики физического состояния зданий и сооружений и изучения их сейсмостойкости на уровне элементов конструкций.
Библиографический список
1. Бендат Дж., Пирсол А. Измерение и анализ случайных процессов. М.: Мир, 1971, 408 с.
2. Еманов А.Ф. О применении вибраторов для определения сейсмостойкости зданий и в микросейсморайонировании. Геология и геофизика, 1995, №7. С.87-92.
3. Еманов А.Ф., Селезнёв В.С. Пересчёт колебаний фильтрами Винера как основа универсального метода обработки сейсмических волн.// Проблемы сейсмологии ІІІ тысячелетия: Материалы международной геофиз. конф., г.Новосибирск, 15-19 сент. 2003г.-Новосибирск: Издательство СО РАН, 2003. – С.207-231.
4. Еманов А.Ф., Селезнёв В.С., Бах А.А. и др. Пересчёт стоячих волн при детальных инженерно-сейсмологических исследованиях / Геология и геофизика, 2002. №2. С.192-207.
5. Еманов А.Ф., Бах А.А., Данилов И.А. и др. Детальные инженерно-сейсмологические исследования плотины Саяно-Шушенской ГЭС. – Вестник Красноярской государственной архитектурно-строительной академии: Сб. науч. тр. Вып.6/ Под ред. В.Д. Наделяева. – Красноярск: КрасГАСА, 2003. C.86-108.
6. Крауфорд Ф. Волны. // Берклиевский курс физики. – Т.3. М.: Наука. 1976, 526 с.
7. Медведев С.В. Инженерная сейсмология. – М.: Изд-во литературы по строительству, архитектуре и строительным материалам. 1962. – 283 с.
8. Фейнман Р., Лейтон Р., Сэндс М. Фейнмановские лекции по физике. Т.3. Излучение. Волны. Кванты// – М.: Мир, 1976. 496 с.
9. Пейн Г. Физика колебаний и волн. – М.: Мир, 1979. 387 с.