Сейсмологические методы используют для регистрации упругие колебания, вызываемые землетрясениями и техногенными воздействиями на геологическую среду, и составляют основу контроля напряженно-деформированного состояния горных пород [4].
За инструментальный период наблюдений с 1963 г. (≈ 60 лет) на территории северо-восточного замыкания системы разломов Тан-Лу зарегистрировано около 650 землетрясений магнитудой M ≥ 2, то есть ежегодно регистрируется около 10 событий, или одно событие в месяц, при которых выделяется энергия порядка 109 Дж.
Высокий сейсмический потенциал данной области определяется сильными землетрясениями, которые произошли здесь в конце XIX и начале XX веков: (23.08.1888 г., M 5,5, 134° В.Д., 50° С.Ш.; 23.12.1914 г., M 6, 139,8° В.Д., 52,5° С.Ш.; 11.03.1924 г., M 5,6, 139,7° В.Д., 51,2° С.Ш.). Имеются сведения о трех исторических землетрясениях с магнитудами M > 6 (900 г., M 6,3, 134,15° В.Д., 50,85° С.Ш.; 1490 г., M 6,5, 134,22° В.Д., 50,53° С.Ш.; 1500 г., M 6,1, 134,22° В.Д., 50,53° С.Ш.). За инструментальный период зарегистрировано также одно сейсмическое событие с магнитудой M = 5,5 и координатами эпицентра (135,1° В.Д.; 51,1° С.Ш.), произошедшее 29.08.1970 г.
Высокая сейсмическая активность в историческом прошлом на юго-западе области в виде 8-балльной зоны и современная активность на востоке, в виде 8–9-балльной меридиональной зоны, поставили ряд актуальных задач по определению реальной сейсмической опасности в ее центральной части. Согласно современным сейсмотектоническим представлениям северо-восточная область системы разломов Тан-Лу относится к зоне 7–8-балльных сотрясений. Детальная схема активных разломов района исследований и прилегающих территорий, составленная коллективом авторов [1], показана на рис. 1, а.
За период наблюдений с 1970 г. за сейсмичностью региона данную область можно представить в виде области сейсмического затишья. Однако это может быть связано с достаточно редкой сетью стационарных сейсмических станций [3], расположением станций в населенных пунктах с высоким уровнем промышленных шумов, что приводит к многочисленным пропускам сейсмических событий.
Одним из способов преодоления указанных трудностей является организация сезонных полевых сейсмологических наблюдений в пунктах, удаленных от промышленных объектов, с использованием современной цифровой широкополосной аппаратуры.
В данной работе для анализа использованы результаты регистрации слабых землетрясений вблизи озера Удыль [4], на северо-восточном замыкании системы разломов Тан-Лу. Выбор места регистрации и современной аппаратуры позволил получить новые данные о распределении эпицентров слабых землетрясений, на основании которых построена предварительная модель затухания упругих волн в системе активных разломов.
Сейсмические наблюдения проводились в течение 10 дней с 23 июля по 1 августа 2014 г. вблизи оз. Удыль (Ульчский р-н Хабаровского края). Используемое оборудование включало цифровой регистратор с непрерывной записью REF TEK-130 и широкополосный сейсмометр REF TEK 150-120.
За время наблюдений в радиусе 360 км было зафиксировано 21 землетрясение, из них 10 сейсмических событий в радиусе 200 км от пункта регистрации (рис. 1, б). Для всех землетрясений были определены значения локальной магнитуды, а для 20 из них также координаты эпицентра по данным одной станции. Для этого были определены значения эпицентрального расстояния и азимута. Четыре землетрясения были достаточно сильными и были уверенно зарегистрированы также на одной или двух стационарных сейсмических станциях, расположенных в п. Ванино, п. Чегдомын, г. Хабаровске и в п. Горный. Это позволило проверить точность определения координат эпицентра по одной станции, которая составила порядка 40 км, что может считаться приемлемым результатом.
Результаты сейсмологических наблюдений
Уровень микросейсмического шума в районе пункта наблюдений позволил регистрировать сейсмические события с магнитудой M < 2. На рис. 2 приведены кривые спектральной плотности мощности (СПМ) для X-канала (ориентированного на север), соответствующие значениям вероятности 5, 50 и 95 % за все время наблюдений.
а б
Рис. 1. а – схема активных разломов района исследований и прилегающих территорий; б – эпицентры зарегистрированных землетрясений и основные тектонические нарушения района исследования
Рис. 2. Спектральные плотности мощности для X-канала, соответствующие значениям вероятности 5, 50 и 95 %, а также модели шума NLNM и NHNM
Также на рис. 2 приведены модели низкого (NLNM) и высокого (NHNM) шума, согласно общепринятой в настоящее время модели из работы [9]. В диапазоне периодов регистрации слабых землетрясений (0,06–0,5 с, частоты 2–16 Гц) медианный уровень микросейсмического шума на 18–25 дБ выше модели NLNM и значительно ниже модели NHNM (необходимо отметить, что указанные модели ограничены минимальным периодом 0,1 с).
Из-за большой удаленности пункта наблюдений от промышленных и транспортных объектов, да и в целом от населенных пунктов, антропогенная составляющая в микросейсмическом шуме практически отсутствовала. В интересующем нас диапазоне периодов наиболее существенным фактором, влияющим на уровень шума, по-видимому, являлся ветер и вызванные им волновые явления на поверхности озера. Кроме того, влияние могут оказывать и другие метеорологические параметры. Метеорологических измерений непосредственно вблизи пункта наблюдений не велось, ближайшая метеостанция находилась в п. Богородское (40 км к северо-востоку). Данные по силе ветра, температуре воздуха и атмосферному давлению, полученные на этой метеостанции, были взяты из архива погоды в Богородском, с сайта http://rp5.ru/Архив_погоды_в_Богородском (01.08.2016) и приведены на рис. 3, б. Дискретность указанных измерений составляет 3 часа.
Было проанализировано изменение во времени среднего значения СПМ на частотах 2–16 Гц (рис. 3, а). Для этого были рассчитаны значения СПМ через каждые две минуты. Полученные значения для устранения коротких случайных флуктуаций, в т.ч. вызванных влиянием сейсмических волн слабых землетрясений, были отфильтрованы скользящим медианным фильтром. Каждый отсчет был заменен на медианное значение 41 соседнего отсчета (т.е. 20 до и 20 после текущего). Таким образом, из ряда значений были исключены вариации короче 80 минут.
Из рис. 3 следует, что некоторые временные отрезки повышенного уровня шума совпадали с повышением силы ветра. Наиболее ярко это проявляется в диапазонах 204,5–208,5 и 212–213 суток. Проведенный корреляционный анализ показал наличие в среднем умеренной корреляционной зависимости между средним значением СПМ и скоростью ветра. Нормализованный коэффициент корреляции составил 0,396, а вероятность, что эта корреляция случайна – 0,06 %.
В то же время значимой корреляции между уровнем шума, с одной стороны, и температурой воздуха и атмосферным давлением, с другой стороны, не обнаружено. Полученные коэффициенты корреляции (от 0,11 до 0,16 по модулю) соответствуют очень низкому уровню корреляции, а их надёжность не слишком высока – от 64 до 83 %.
На рис. 4 показана сейсмограмма одного из землетрясений с магнитудой М = 1,1.
Результаты исследования и их обсуждение
Одной из задач, решаемых высокоточной сейсмометрией, является определение затухания сейсмических волн в различных азимутах и определение анизотропии среды. Для определения анизотропии применялся анализ отношений амплитуд P и S-волн. Подразумевалось, что на близких расстояниях (до 300 км) объемные волны от коровых землетрясений регистрируются как прямые, распространяющиеся по законам геометрической оптики. Затухание поперечной S-волны существенным образом зависит от того, вдоль или поперек разломных структур она распространяется. В первом случае затухание будет значительно меньше, чем во втором, поскольку поперечная волна вдоль разломов распространяется как в волноводе.
Соотношение амплитуд P и S-волн было определено для 19 землетрясений из 21. Амплитуды измерялись в смещениях, в частотном диапазоне 1–20 Гц или более узком, если на некоторых частотах интенсивность шума превышала интенсивность сейсмических колебаний. Полученные значения нанесены на радиальную диаграмму (рис. 5, а). Соотношения амплитуд, полученные в каждом квадранте, были усреднены, после чего был построен эллипс, проходящий через полученные четыре точки (рис. 5, б).
а
б
Рис. 3. а – изменение во времени среднего значения спектральной плотности мощности на частотах 2–16 Гц (X-канал); б – графики изменения скорости ветра, температуры воздуха и атмосферного давления на метеостанции в п. Богородское
Рис. 4. Сейсмограмма землетрясения с магнитудой М = 1,1, произошедшего 23.07.2014 г. в 22:48:05 UTC. Сверху вниз показаны вертикальный, северный и восточный каналы. Вертикальными линиями отмечены вступления P и S-волн
а б
Рис. 5. Отношения амплитуд P и S-волн: а – радиальная диаграмма; б – радиальная диаграмма с усреднением по квадрантам (звёздочки, соединенные прямыми линиями) и описанным вокруг полученного четырехугольника эллипсом
Ориентация большой оси полученного эллипса соответствует направлению ЮЗ-СВ. Это хорошо согласуется с пространственной ориентировкой активных разломов районе оз. Удыль (рис. 1, а). Основные тектонические нарушения в этом районе – Лимурчанский и Удыльский разломы характеризуются соответственно меридиональным и северо-восточным простиранием.
Привлечение материалов по блоковой структуре земной коры, распределению аномалий геофизических полей, сейсмотектоническим реконструкциям полей напряжений по фокальным механизмам и линеаментному анализу цифровых моделей рельефа показало, что северо-восточное простирание главной оси эллипса в среднем азимуте 45° согласуется с пространственной ориентацией сейсмоактивных границ доменов [6], азимутами максимальных деформаций [5], ориентацией оси главных напряжений сжатия [7] и направлением систем линеаментов [10]. Новые данные по ориентации осей сжатия-расширения [8] и разработанная тектоническая кластерная модель по распределению сильных землетрясений [10] позволяют сделать вывод, что применение описанной методики вполне допустимо для построения модели анизотропии на расстояниях эпицентров землетрясений до 300 км от пункта регистрации при достаточной представительности количества сейсмических событий. Проведение данных работ будет способствовать уточнению сейсмической опасности, что становится актуальным при интенсивном развитии инфраструктуры Дальневосточного региона [2].
Авторы благодарны д.г.-м.н. А.Н. Диденко за методическую помощь в организации полевых наблюдений. Исследование выполнено при финансовой поддержке РФФИ в рамках научного проекта № 16-05-00097а.