CC BY
20УДК 550.348.435
СПЕКТРАЛЬНЫЕ И ДИНАМИЧЕСКИЕ ПАРАМЕТРЫ ОЧАГОВ ЗЕМЛЕТРЯСЕНИЙ КРЫМА 2017 ГОДА
Пустовитенко Б. Г., Эреджепов Э. Э., Бондарь М. Н.
Институт сейсмологии и геодинамики ФГАОУ ВО «Крымский федеральный университет имени В. И. Вернадского», Симферополь, Республика Крым, Россия
E-mail: bpustovitenko@mail.ru
Приведены спектральные и динамические параметры очагов (Мо, ro, Act, 8, г|ст, Астг, и, Eu и Mw) десяти землетрясений Крыма 2017 года в диапазоне энергетических классов Кп =6.5^10.9, восстановленные по спектрам записей продольных и поперечных сейсмических волн, зарегистрированных цифровыми региональными сейсмическими станциями. Всего получено 78 станционных амплитудных спектров, положенных в основу расчета динамических параметров очагов с использованием теоретической дислокационной модели Брюна. Дается анализ полученных результатов.
Ключевые слова: амплитудный спектр, модель Брюна, спектральная плотность, угловая частота, сейсмический момент, радиус дислокации, сброшенное и кажущееся напряжения, подвижка по разрыву, радиационное трение.
ВВЕДЕНИЕ
В 2017 г. региональной сетью станций Крыма зарегистрировано 89 местных землетрясений, часть из которых относится к представительным для региона с Кп>8.0 [1]. Инструментальные записи сейсмических колебаний послужили основой не только определения основных параметров землетрясений (время, место, энергетический класс и магнитуда), но и дополнительных характеристик источника излучения: спектральных и динамических параметров очагов.
К динамическим параметрам очага относятся: скалярный сейсмический момент М0, радиус круговой дислокации Г0, сброшенное напряжение Аст, величина деформации сдвига 8, кажущееся напряжение пст, величина радиационного трения Астг, средняя подвижка по разрыву U (или величина дислокации), Еи - энергия образования дислокации в очаге. Для восстановления этих параметров используются количественные связи, установленные между спектральными характеристиками амплитудных спектров объемных волн, зарегистрированных на сейсмических станциях и существующими теоретическими моделями источников.
В данной работе, как и в предыдущие годы [2-4], в качестве теоретического источника использована дислокационная модель Брюна [5], в рамках которой хорошо описываются экспериментальные данные по спектрам записей местных землетрясений Крыма [6].
1. ИСХОДНЫЕ ДАННЫЕ
Для расчета спектров и определения на их основе динамических параметров очагов выбраны все наиболее значимые сейсмические события 2017 года с Кп>9.0, а также более слабые толчки, произошедшие в непосредственной близости к Южному
69
берегу Крыма и в менее изученных районах региона: Севастопольском (№1) и Степном Крыму (№6). По такому критерию первоначально были отобраны и просмотрены сейсмограммы большинства землетрясений из этих районов по всем региональным станциям с оценкой качества записей, пригодных для расчета спектров. В результате анализа для решения обратной задачи восстановления динамических параметров очагов выбраны сейсмограммы пяти сейсмических станций («Алушта», «Севастополь» «Симферополь», «Судак», «Ялта») десяти землетрясений с энергетическими классами Кп=6.5-10.9 (Рис. 1). Основные параметры изученных землетрясений по данным сводной обработки в Крыму приведены в таблице 1.
Рис. 1. Карта эпицентров землетрясений Крымско-Черноморского региона за 2017 г., для которых восстановлены динамические параметры: 1 - энергетический класс Кп [7]; 2 - глубина очага; 3 - сейсмическая станции; 4 - граница и № районов. Цифры рядом с эпицентрами - номера землетрясений в соответствии с Табл. 1.
Наибольшее количество рассмотренных землетрясений (3 события), относится к Керченско-Анапской зоне очагов (район №5), по два - к Севастопольскому и Ялтинскому (районы №1 и №2) , по одному к Алуштинскому (район №3), Степному Крыму (район №6) и к Черноморской впадине (район №9).
Таблица 1.
70
Основные параметры землетрясений Крыма за 2017 год, для которых восстановлены динамические параметры очагов
№ Дата, t0, Эпицентр h, Магнитуда Район
Э м ч мин с ф°, N Г, E км Мс Mw/n, MSH Кп
Табл. 3
1 2 3 4 5 6 7 8 9 10 11
1 27.01 19 05 25.9 44.35 34.37 29 3.0/1 8.1 Ялтинский (№ 2)
2 16.06 18 17 34.8 44.75 36.74 9 3.6 4.0/9 3.7 10.9 Керченско-Анапский (№ 5)
3 05.07 20 12 41.4 43.94 34.10 47 2.5 3.5/6 3.1 9.0 Ялтинский (№ 2)
4 09.07 20 37 12.1 44.62 34.97 24 2.5/2 6.5 Алуштинский (№ 3)
5 13.07 12 49 40.5 44.81 32.42 7 2.4/2 7.4 Севастопольский (№ 1)
6 08.08 04 08 50.5 44.68 36.77 8 3.1 3.3/8 3.4 9.9 Керченско-Анапский (№ 5)
7 08.08 05 32 20.3 44.59 36.77 5 3.4/8 3.6 10.2 Керченско-Анапский (№ 5)
8 30.08 20 24 40.9 45.30 34.19 25 2.4/3 7.7 Степной Крым (№ 6)
9 02.10 13 58 34.0 43.05 36.44 30 3.5/4 3.3 9.6 Черноморская впадина(№ 9)
10 12.11 12 58 29.2 44.86 32.29 9 2.9 3.6/5 3.4 9.8 Севастопольский (№ 1)
Примечание. Параметры землетрясений в графах 2-7, 9-11 даны по данным сводной обработки в Крыму; значения Mw/n - из Табл. 2
Для расчета амплитудных спектров и восстановления по ним динамических параметров очагов использовано 16 записей продольных (P) и 62 - поперечных (5) волн на сейсмических станциях «Алушта» (ALU), «Севастополь» (SEV), «Симферополь» (SIM), «Судак» (SUDU) и «Ялта» (YAL). Для этих станций в 2017 году была проведена калибровка сейсмографов и уверенно определены амплитудно-частотные характеристики.
Наибольшее количество спектров получено по записям сейсмических станций «Алушта» и «Судак». Процент участия этих же станций в общей оценке динамических параметров отдельных землетрясений также наибольший, соответственно 90% и 80%. Меньше всего для анализа привлекались записи станции «Ялта» (50%), особенно для относительно слабых землетрясений. Большая часть динамических параметров очагов землетрясений восстановлена по характеристикам амплитудных спектров поперечных волн (5), амплитуды которых практически всегда в несколько раз превышают амплитуды продольных волн (P). Для анализа отобраны только записи с четкими фазами объемных P- и 5-волн с превышением полезного сигнала над фоном помех в два и более раза. Примеры таких записей даны на Рис. 2.
Как и в предыдущие годы [2-4] не использованы для анализа записи сейсмостанции «Феодосия», которая расположена на мергелистых глинах с пониженными прочностными свойствами подстилающих пород и повышенным затуханием высоких частот. Сейсмические колебания от местных землетрясений на этой станции существенно осложнены микросейсмическим шумом и не отвечают критерию выбора записей для анализа. Сейсмограммы новой станции «Тарханкут» также пока не привлечены для исследования.
71
а
б
Рис. 2. Примеры записей землетрясений 2017 г. по составляющим Ъ, N-S, E-W (сверху вниз): а - 9 июля с Кп = 6.5 (№4 по Табл. 1); б - 30 августа с Кп =7.7 (№8 по Табл. 1) соответственно на станциях «Судак», Д=30 км и «Симферополь», Д=44 км.
2. МЕТОДИКА ОБРАБОТКИ И ИНТЕРПРЕТАЦИЯ АМПЛИТУДНЫХ СПЕКТРОВ
Станционные спектры поперечных волн (5) рассчитывались по полному вектору колебаний на двух горизонтальных компонентах и (Е^), а
продольных (Р) - по вертикальной составляющей (Ъ). Относительная длительность
72
т исследуемого участка записи принята равной интервалу времени от момента вступления S- и Р-волн до времени спада их максимальных амплитуд Атах на уровень 7з Атах [6]. Спектры рассчитаны стандартным методом быстрого преобразования Фурье. Всего по записям объемных волн, отобранных для анализа десяти землетрясений, рассчитано 78 амплитудных спектров по разным типам волн и составляющим записи на сейсмограмме. Примеры спектров даны на Рис. 3.
ю
10
10
10
О, М- I
10
1 1 аьц N ;
: : 1 : \ А1
; ..... 1 У.1 . 1 1
г ;; .......;; \
10
10
ю1 *,Гц
.6 Г2, м- с
«Г1 10° юЧГц
27 января в 19ь05т25.98 (№ 1)
16 июня в 18ь17т34.85 (№ 2)
ю
10
10
О,»..
8Ши, N
5 июля в 20ь12т41.4 5 _(№ 3).
9 июля в 20ь37т12.85 (№4)_
13 июля в 12М9т40.58 (№5)_
73
-в О. м • с
10° 101Г,Гц
гл
10
10
10
10
м- <
10
7 Г Т Г Ч ГГЦ УЛ1д Е
„з ;,
ю-1 10° ю1 Г, Гц
8 августа в 04ь08т50.55 (№6)
8 августа в 05ь32т20.35 (№7)
О, м-
а,м.
30 августа в 20ь24т40.98 _(№8)_
2 октября в 13ь58т34.08 _(№9)_
12 ноября в 12ь58т29.28 _(№10)_
Рис. 3. Примеры амплитудных спектров объемных сейсмических волн землетрясений Крыма за 2017 г. по записям региональных цифровых сейсмостанций и их аппроксимация в рамках теоретической модели Брюна. Номера и даты землетрясений соответствуют таковым в таблице 1.
Согласно теоретической дислокационной модели Брюна (ю-2) [5] амплитудный спектр объемной волны характеризуется тремя основными параметрами: спектральной плотностью О0(х, /) в длиннопериодной части (при /^0), угловой частотой/0 (ю 0) и углом наклона у—2. Из приведенных примеров (Рис. 3) видно, что спектры реальных записей землетрясений Крыма за 2017 г. уверенно аппроксимируются в рамках выбранной теоретической модели (ю-2).
74
3. РАСЧЕТ ДИНАМИЧЕСКИХ ПАРАМЕТРОВ ОЧАГА
Для восстановления очаговых параметров землетрясений в рамках дислокационной модели Брюна [5] используются две главные характеристики амплитудных спектров: спектральная плотность Qо, которая пропорциональна скалярному сейсмическому моменту Мо (1) и угловая частота f0, прямо связанная с размером дислокации го (2).
Mo = Q о(4лр-V3)/R 0Ф-G (А, h) • C (и) • Sm(f) (1)
где Qо- максимальное значение спектральной плотности при и ^ 0; V - скорость распространения объемной волны; р - плотность пород в окрестности очага, Я0ф -направленность излучения из очага на станцию регистрации; G(A, h) - поправка за геометрическое расхождение; С(ю) - частотная характеристика среды под станцией, Smf - поправка за неупругое затухание в мантии.
Методика учета всех этих факторов при переходе от станционного спектра к спектру источника подробно изложена в работах [6, 8, 9] и не изменилась по сравнению с предыдущими годами.
Для дислокационной модели Брюна с разрывом в виде круга, радиус дислокации го вычисляется по угловой частоте fо по формуле:
r = 2.34-^— (2)
о 2nfü У }
По найденным значениям сейсмического момента Мо и размеров разрыва го определены средние величины напряжений в очаге: сброшенное До и кажущееся По, средняя подвижка по разрыву U по формулам из работ [1о, 11]:
Л 7M о _ и • ES - M о Ас =-3 т]а=——- U =-^
16го M о и •л го
где и - модуль сдвига; п - сейсмический коэффициент полезного действия; Es-энергия источника (lg Es=11.8+1.5Mlh).
Величина радиационного трения Acr, связанного с шероховатостью поверхности разрыва в очаге, определяется по формуле из [Ю]:
Acr =1/2Ас -Г)С (3)
Оценка энергии Еи выполнена по дислокации в очаге и сброшенному напряжению [5]:
1 — 2 Eu = — AcU •лго 2 о
С использованием уравнения из теории упругости, связывающего деформацию с напряжением, оценена упругая деформация сдвига в окрестности очага:
75
Да £ --
М
Моментная магнитуда рассчитана по индивидуальным оценкам сейсмических моментов М0 с использованием формулы Канамори [12]:
Ш =2/3-(^ М0+7)-10.7
где М0 имеет размерность Н-м.
Вычисление средних значений параметров и стандартных отклонений для конкретного землетрясения выполнено по станционным определениям с учетом логнормального закона распределения величин [6]. Станционные параметры радиационного трения Даг получились отрицательными или знакопеременными, в связи с этим их среднее значение вычислено по среднегеометрическим для данного очага напряжениям Да и п^ по формуле (3). В этом случае графа «стандартное отклонение» в таблице 2 не заполнялась. Средняя величина моментной магнитуды Мн> определена как среднее арифметическое с соответствующей погрешностью.
Результаты расчета станционных и средних значений спектральных и динамических параметров очагов землетрясений за 2017 год даны в таблице 2. Для каждой станции указаны эпицентральное расстояние А, км, тип использованной волны (Р, 5) и составляющая записи, где ^+Е) означает полный вектор колебаний по N+S и E-W.
4. ОБСУЖДЕНИЕ РЕЗУЛЬТАТОВ
Количественные параметры амплитудного станционного спектра зависят от энергии в очаге, направленности излучения на станцию регистрации, свойств глубинной среды на пути распространения сейсмических волн и под станцией регистрации, а также от ряда факторов, которые трудно учесть. Поэтому для надежных оценок очаговых параметров землетрясений Крыма за 2017 год использованы записи как продольных (Р), так и поперечных (5) волн четырех региональных станций. В осреднении очаговых параметров в большинстве случаев участвовало от четырех до девяти станционных определений. Только для одного относительно слабого сейсмического толчка (№1), с энергетическим классом Кп =8.1, произошедшего в Ялтинском районе, динамические параметры рассчитаны по 5-волне на ближайшей к очагу станции «Алушта».
Наилучшая сходимость станционных определений получена для радиуса круговой дислокации (5г0 = 0.01^0.04), величина которого в явном виде не зависит от направленности излучения энергии из очага и условий среды под станцией регистрации. Стандартные отклонения по другим параметрам для большинства землетрясений 8Х<0.2. Наибольший разброс станционных определений, отмечен для слабого сейсмического события (Кп =7.7), произошедшего 30.08.2017 г. в районе Степного Крыма (№8). Вероятно, это связано с разной направленностью очагового излучения Reф (формула 1) разных типов волн на станцию регистрации «Алушта». Так значения сейсмического момента М0 по спектрам Р- и 5-волн получились
76
различающимися на порядок (Табл. 2), что при близких значениях радиуса круговой дислокации го повлекло за собой к разбросу данных по другим динамическим параметрам на порядок и более.
Таблица 2.
Спектральные и динамические параметры очагов землетрясений Крыма за 2017 год
Станция на нл о т Составляющая < о N О о а о N © о а £ 3 £ гя ~ о ^ о 1Л 0 1 о ш 3 сч ~ О 'а £ О Ь р О о <1 00 ^ 0 1 * 5
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16
№ 1. 27 января, ¿0=19 ч 05 мин 25.9 с; ф = 44.35°, Х= 34. 37°; h = 29 км; Ял =8.1
АЦП | S |№Е| 37 | 10.32 | 3.85 | 4.59 |0.37 |4.04 | 13.5 | 0.36 | 0.83 | 1.2 | 3.1 | 3.0
№ 2. 16 июня, ¿0=18 ч 17 мин 34,8 с; ф = 44.72°, X = 36. 14°; h = 9 км; Ял =10.9
АЦи S Е 185 1.8 1.5 263 0.89 16.1 53.8 3.5 2.55 5.51 708 4.22
SEV р Z 243 0.12 2.4 116 0.74 12.3 41 2.22 5.79 0.36 238 3.98
SEV-1 S N+E 243 0.56 1.68 117 0.8 10.1 33.6 1.95 5.74 0.7 196 3.98
SEV S N+E 243 0.8 1.6 176 0.84 12.4 41.4 2.53 4.02 0.22 346 4.09
SUDU р Z 139 0.18 2.4 102 0.93 5.51 18.4 1.24 6.61 -3.86 93.2 3.94
YAL р Z 206 0.1 2.45 84.6 0.91 4.9 16.3 1.08 7.93 -5.5 68.9 3.89
YAL S N+E 206 0.54 1.6 98.7 0.84 7.34 24.5 1.49 6.8 -3.13 121 3.93
SIM р Z 208 0.12 2.41 102 0.93 5.51 18.4 1.24 6.61 -3.86 93.2 3.94
SIM S N+E 208 0.54 1.6 98.7 0.84 7.34 24.5 1.49 6.8 -3.13 121 3.93
X- среднее значение 120.56 0.86 8.36 27.88 1.73 5.6 -1.42 166.8 3.99
стандартное отклонение 0.05 0.01 0.06 0.06 0.06 0.05 0.11 0.07
№ 3. 5 июля, ¿0=20 ч 12 мин 41.4 с; р = 43.94°, X = 34. 10°; h = 47 км; Ял =9.0
АЦи р Z 86 0.02 3.1 7.38 0.99 0.34 56.3 0.04 5.76 -5.59 0.21 3.18
АЦи S №Е 86 0.91 2.17 51.9 0.76 5.27 8.78 0.48 0.82 1.82 22.8 3.75
SEV S N+E 75 0.14 2.15 10.9 0.76 1.08 1.79 0.1 3.9 -3.37 0.97 3.29
SUDU р Z 127 0.05 3.0 37.7 1.02 1.56 2.61 0.2 1.13 -0.34 4.95 3.65
SUDU S N+E 127 0.32 2.2 37.3 0.74 3.95 6.58 0.36 1.14 0.84 12.3 3.65
YAL S N+E 61 0.34 2.2 19.0 0.74 2.01 3.35 0.2 2.24 -1.23 3.18 3.46
X-среднее значение 21.94 0.83 1.7 6.09 0.17 1.94 -1.09 3.11 3.5
стандартное отклонение 0.14 0.03 0.17 0.22 0.16 0.14 0.3 0.19
№ 4. 9 июля. ¿0=20 ч 37 мин 12.1 с; р = 44.62°, X = 34. 97°; h = 24 км; Ял =6.5
АЦП S №Е 46 0.03 3.6 0.61 0.37 0.52 1.73 0.05 0.27 -0.015 0.053 2.46
SUDU S N 30 0.05 3.8 0.76 0.35 0.76 2.53 0.065 0.22 0.15 0.096 2.52
X-среднее значение 0.68 0.36 0.63 2.1 0.06 0.24 0.08 0.07 2.49
стандартное отклонение 0.05 0.01 0.08 0.08 0.06 0.04 0.13 0.03
№ 5. 13июля, ¿0=12 ч 49 мин 40.5 с; ф = 44.81°, X = 32. 42°; h = 7 км; Ял =7.4
SEV S N+E 104 0.012 3.1 0.35 0.38 0.27 0.88 0.025 2.75 -2.62 0.015 2.3
SUDU S N+E 203 0.015 3.15 0.84 0.32 0.68 2.26 0.062 1.13 -0.79 0.095 2.55
77
Продолжение таблицы
1 |2|3|4| 5 |6| 7 8 9 10 11 12 13 14 15 16
S №Е 187 0.45 2.2
X-среднее значение 0.54 0.35 0.43 1.41 0.04 1.76 -1.55 0.04 2.42
стандартное отклонение 0.19 0.04 0.2 0.2 0.2 0.19 0.4 0.12
№ 6. 8 августа, ¿¡=04 ч 08 мин 50.5 с; (р = 44.68°, X = 36. 77°; h = 8км; ЯЛ =9.9
ALU S N+E 187 0.45 2.2 32.9 0.61 6.36 21.2 0.94 2.57 0.61 34.9 3.61
SEV Р Z 245 0.02 4.0 10.0 0.58 2.23 7.42 0.31 8.44 -7.32 3.72 3.27
SEV S №Е 245 0.14 2.3 13.4 0.58 2.96 9.87 0.42 6.31 -4.83 6.62 3.35
SIM Р Z 211 0.025 3.1 10.8 0.75 1.12 3.73 0.2 7.81 -7.25 2.02 3.29
SIM S N+E 211 0.14 2.1 11.6 0.64 1.95 6.49 0.3 7.3 -6.3 3.76 3.31
SUDU Р Z 141 0.02 4.0 5.79 0.58 1.28 4.28 0.18 14.6 -14.0 1.24 3.11
SUDU S N+E 141 0.4 2.2 22.1 0.61 4.27 14.2 0.63 3.82 -1.69 15.8 3.5
YAL S №Е 208 0.09 2.9 9.77 0.46 4.33 14.4 0.48 8.66 -6.5 7.04 3.26
X-среднее значение 12.73 0.6 2.61 8.7 0.38 6.64 -5.34 5.55 3.34
стандартное отклонение 0.08 0.02 0.09 0.09 0.09 0.08 0.17 0.11
№ 7. 8 августа, ¿0=05 ч32 мин 20.3 с; р = 44.59°, X = 36. 77°; h = 5км; Кп =10.2
ALU S N+E 188 0.5 2.8 31.0 0.45 14.7 48.9 1.61 5.44 1.89 75.7 3.6
SEV Р Z 245 0.03 3.95 13.4 0.57 3.23 10.8 0.44 12.6 -11.0 7.19 3.35
SEV S N+E 245 0.16 3.0 13.0 0.42 7.54 25.1 0.77 13.0 -9.24 16.3 3.35
SIM S N+E 214 0.16 2.75 11.3 0.46 5.08 16.9 0.57 14.9 -12.4 9.58 3.31
SUDU Р Z 144 0.05 4.05 13.1 0.55 3.41 11.4 0.46 12.9 -11.2 7.45 3.35
SUDU S N+E 144 0.43 2.8 20.5 0.45 9.7 32.3 1.06 8.23 -3.38 33.1 3.48
YAL Р Z 208 0.02 0.04 7.54 0.55 1.96 6.54 0.26 22.4 -21.4 2.47 3.19
YAL S N+E 208 0.15 3.0 10.3 0.42 5.99 20.0 0.61 16.4 -13.4 10.6 3.28
X-среднее значение 13.74 0.48 5.38 17.95 0.63 12.28 -9.59 12.37 3.36
стандартное отклонение 0.07 0.02 0.1 0.1 0.09 0.07 0.16 0.06
№ 8. 30 августа, ¿0=20 ч 24 мин 40.9 с; ч р = 45.30°, X = 34.19°; h = 25 км; Кп =7.7
ALU Р Z 75 0.0012 5.8 0.13 0.4 0.08 0.29 0.01 10.5 -10.5 0.02 2.01
ALU S №Е 75 0.07 3.1 1.42 0.43 0.77 2.56 0.08 0.09 -0.06 0.18 2.7
SIM S N+E 44 0.09 3.05 0.9 0.44 0.46 1.55 0.05 1.49 -1.26 0.07 2.57
X-среднее значение 0.55 0.42 0.3 1.05 0.03 1.12 -0.97 0.06 2.43
стандартное отклонение 0.32 0.01 0.3 0.29 0.27 0.6 0.28 0.28
№ 9. 2 октября, ¿0=13 ч 58 мин 34.0 с; р = 43.05°, X = 36.44°; h = 30 км; Кп =9.6
ALU Р Z 245 0.01 3.0 12.1 0.82 0.96 3.2 0.19 4.96 -4.48 1.93 3.33
YAL S N 244 0.056 2.3 12.9 0.62 2.42 8.08 0.36 4.64 -3.43 5.21 3.34
SUDU S N+E 236 0.18 2.0 39.8 0.71 4.92 16.4 0.84 1.5 0.95 32.6 3.67
SIM S №Е 282 0.085 2.1 22.5 0.67 3.21 10.7 0.53 2.67 -1.06 12.0 3.5
X-среднее значение 19.34 0.7 2.46 8.21 0.42 3.1 -1.87 7.92 3.46
стандартное отклонение 0.12 0.03 0.15 0.15 0.14 0.12 0.26 0.15
№ 10. 12 ноября, ¿0=12 ч 58 мин 29.2 с; р = 44.86°, X = 32.29°, h = 9 км; Кп =9.8
ALU S N+E 168 0.71 2.1 32.8 0.57 7.85 26.2 1.08 2.58 1.34 42.9 3.61
SEV S N+E 115 0.56 2.3 18.1 0.52 5.7 19.0 0.72 4.87 -1.95 17.2 3.44
SIM S N+E 144 0.7 1.9 28.0 0.63 4.97 16.6 0.76 3.02 -0.5 23.2 3.57
SUDU Р Z 213 0.1 3.3 38.7 0.68 5.45 18.2 0.89 2.19 0.5 35.2 3.66
SUDU S №Е 213 0.7 1.9 28.0 0.63 4.97 16.6 0.76 3.02 -0.53 23.2 3.57
X-среднее значение 28.26 0.06 5.7 19.03 0.83 3.02 -0.17 26.86 3.48
стандартное отклонение 0.05 0.02 0.04 0.04 0.03 0.06 0.07 0.15
78
Наибольшие значения динамических параметров (Мо, го, Ас, в, па, й, Ей и Mw) получены для землетрясения 16 июня с Кп =10.9, произошедшего в восточной части региона (район 5) на глубине h = 9 км (№2), а наименьшие - для слабых толчков 9 июля с Кп = 6.5 (№4) и 30 августа с Кп = 7.7 (№8). В целом, с ростом энергетического уровня землетрясений количественные значения динамических параметров очагов
увеличивались в пределах разброса измеряемых величин. Радиационное трение Асг для исследованных землетрясений имело отрицательное значение, указывая на сложное скольжение разрыва в очаге за счет неравномерного распределения прочностных свойств глубинной среды в очаговых зонах. Во всем диапазоне энергий средняя величина сброшенных напряжений не превысила Да = 8-105 Па (8 бар), а кажущихся напряжений па < 13 -105 Па (13 бар).
Анализ индивидуальных станционных определений по данным за 2017 год показал, что значения сейсмических моментов М0 и, соответственно, моментной магнитуды Мц>, восстановленные по записям продольных волн Р в большинстве случаях получились существенно ниже, чем по записям поперечных 5": среднее по всем станциям значение отношений k =Мо 5 / Мо Р получилось равным 2.5 при максимальном 11 по данным станции «Алушта» для очага землетрясения №8.
Сравнение полученных в 2017 г. динамических параметров очагов со средними их долговременными величинами [13], как и в предыдущие годы [2-4], проведено только для сейсмического момента М0 и радиуса круговой дислокации Г0 (Рис. 4). Эти параметры определяются непосредственно по характеристикам спектров записей сейсмических колебаний на станциях регистрации по формулам 1 и 2, а остальные параметры вычисляются по их значениям. Для сравнения использованы зависимости М0 (Кп) и Г0 (Кп) из [13]:
^М0 = 0.645(±0.027)Кп + 15.142(±0.271), р=0.99 ^ Г0 = 0.112(±0.011)Кп - 1.293(±0.107), р=0.93,
где р - коэффициент корреляции.
Рис. 4. Сравнение динамических параметров очагов землетрясений Крыма за 2017 г.: а - сейсмического момента М0; б - радиуса круговой дислокации Г0 с долговременными зависимостями Мо(Кп) и го(Кп) [13]. Пунктиром обозначены пределы погрешностей долговременных зависимостей.
79
Полученные в 2017 г. средние значения сейсмических моментов и радиусов круговой дислокации (Рис. 4) попадают в доверительный интервал средних долговременных параметров. Наибольшее отличие от средней аппроксимирующей прямой Mo(Ru) отмечено для землетрясений 9 июля (№4) с Кп = 6.5 и №7 с Кп =10.2. В первом случае Mo получилось завышенным, а во втором - заниженным.
ЗАКЛЮЧЕНИЕ
По результатам сейсмического мониторинга в 2017 г, база данных по очаговым параметрам местных землетрясений Крыма дополнена 47 станционными определениями для 10 сейсмических событий в диапазоне энергетических классов Кп=6.5^10.9, произошедших в различных районах региона на эпицентральных расстояниях от 30 км до 245 км от станций регистрации. Для достоверных оценок спектральных и динамических параметров очагов землетрясений за 2017 год были использованы четкие записи объемных сейсмических волн, не осложненные микросейсмическими помехами на сейсмических станциях с надежными амплитудно-частотными характеристиками регистрирующей аппаратуры. Средние значения динамических параметров очагов в большинстве случаев оценены по группе станций и по разным типам волн, что обеспечило в основном малую величину показателя степени рассеяния индивидуальных оценок. В этой связи, количественные оценки очаговых параметров землетрясений за 2017 год можно отнести к категории надежных для дальнейшего их использования в научных и научно-прикладных задачах, в том числе при моделировании прогнозных сейсмических воздействий в соответствии с методикой [14, 15].
Список литературы
1. Свидлова В. А., Сыкчина З. Н., Бондарь М. Н., Бойко В. А. Сейсмичность Крыма в 2016 году // Ученые записки Крымского федерального университета им. В. И. Вернадского. Серия «География. Геология». 2016. Том 3(69). N4. С. 7-50.
2. Пустовитенко Б. Г., Мержей Е. А., Пустовитенко А. А. Калинюк И. В. Очаговые параметры землетрясений Крыма 2013 года // Сейсмологический бюллетень Украины за 2013 г. Севастополь: НПЦ «ЭКОСИ Гидрофизика», 2014. С. 12-21.
3. Пустовитенко Б. Г., Эреджепов Э. Э. Спектральные и динамические параметры землетрясений Крыма 2015 года // Ученые записки Таврического национального университета им. В. И. Вернадского; серия «География. Геология». 2016. Том 2(68). N4. 2016. С. 53-68.
4. Пустовитенко Б. Г., Эреджепов Э. Э. Очаговые параметры землетрясений Крымско-Черноморского региона 2016 года // Ученые записки Таврического национального университета им. В. И. Вернадского; серия «География. Геология». 2017. Том 3(69). N4. С. 51-69.
5. Brune I. V. Tectonic stress and the spectra of seismic shear waves from earthquakes // J. Geophys. Res. 1970. V. 75, no. 26, pp. 4997-5009.
6. Пустовитенко Б. Г., Пантелеева Т. А. Спектральные и очаговые параметры землетрясений Крыма. Киев: Наукова думка, 1990. 251 c.
7. Пустовитенко Б. Г., Кульчицкий В. Е. Об энергетической оценке землетрясений Крымско-
Черноморского региона // Магнитуда и энергетическая классификация землетрясений. М.: ИФЗ АН СССР. 1974. Т. 2. С. 113-125.
80
8. Аптекман Ж. Я., Белавина Ю. Ф., Захарова А. И., Зобин В. М., Коган С. Я., Корчагина О. А., Москвина А. Г., Поликарпова Л. А., Чепкунас Л. С. Спектры Р-волн в задаче определения динамических параметров очагов землетрясений. Переход от станционного спектра к очаговому и расчет динамических параметров очага // Вулканология и сейсмология. 1989. N2. C. 66-79.
9. Аптекман Ж. Я., Дараган С. К., Долгополов В. В., Захарова А. И., Зобин В. М., Коган С. Я., Корчагина О. А., Москвина А. Г., Поликарпова Л. А., Чепкунас Л. С. Спектры Р-волн в задаче определения динамических параметров очагов землетрясений. Унификация исходных данных и процедуры расчета амплитудных спектров // Вулканология и сейсмология. 1985. N2. С. 60-70.
10. Костров Б. В. Механика очага тектонического землетрясения. М.: Наука, 1975. 179 с.
11. Кейлис-Борок В. И. Исследование источников, приближенно эквивалентных очагам землетрясений // Труды Геофизического Института АН СССР. 1959. N9(136). C. 20-42.
12. Hanks T. C., Kanamori H. A. Moment magnitude scale // J. Geophys. Res. 1979, v. 84, no. 35, pp. 2348-2350.
13. Пустовитенко Б. Г., Пустовитенко А. А., Капитанова С. А., Поречнова Е. И. Пространственные особенности очаговых параметров землетрясений Крыма. Сейсмичность Северной Евразии. Обнинск: ГС РАН, 2008. С. 238-242.
14. Пустовитенко Б. Г., Кульчицкий В. Е., Калинюк И. В., Мержей Е. А. Сравнительные характеристики акселерограмм, смоделированных для разных сеймотектонических условий // Сейсмостойкое строительство. Безопасность сооружений. 2012. N6. C. 20-23.
15. Пустовитенко Б. Г. Мержей Е. А. К прогнозу ожидаемых сейсмических воздействий от сильных землетрясений Крыма // Доповщ Национально! Академй Наук Украши. 2014. N5. C.102-109.
SPECTRAL AND DYNAMIC PARAMETERS OF THE FOCI OF 2017 CRIMEAN EARTHQUAKES
Pustovitenko B. G., Eredzhepov E. E., Bondar M. N.
Institute of seismology and geodynamics FSAEI HE «Of the V. I. Vernadsky Crimean Federal University», Simferopol, Republic of Crimea, Russia E-mail: bpustovitenko@mail.ru
The spectral and dynamic parameters of foci (М0, Г0, Ac, s, no, Acr, ü, Eu and Mw) of ten 2017 Crimea earthquakes in the range of energy classes Кп=6.5^10.9 are given. Focal parameters reconstructed from the recordings of longitudinal (P) and transverse (5) seismic wave records recorded by digital regional seismic stations: «Alushta» (ALU), «Sevastopol» (SEV), «Simferopol» (SIM), «Sudak» (SUDU), «Yalta» (YAL). A total of 78 station amplitude spectra were calculated. To approximate the spectra and calculate the dynamic parameters of the foci, the Brune dislocation model (ю-2) was used. For reliable estimates of the average values of the focal parameters of each earthquake, in most cases, from four to nine station definitions are used. Only for one relatively weak seismic shock on January 27 with the energy class Кп = 8.1, which occurred in the Yalta zone, the dynamic parameters are calculated using the 5-wave at the closest to the source of the «Alushta» station.
The average values of seismic moments М0 and radii of circular dislocations Г0 obtained in 2017 corresponded to confidence intervals of average long-term parameters. The greatest difference from the average approximating straight line Mci(Kn) was observed for
81
earthquakes on July 9 with An = 6.5 and on August 8 with An = 10.2. In the first case, Mo turned out to be overestimated, and in the second - underestimated. The largest values of dynamic parameters (Mo, ro, Ac, s, n^, u, Eu and Mw)) were obtained for the June 16 earthquake with An = 10.9, which occurred in the eastern part of the region at a depth of h = 9 km, and the smallest - for weak shocks on July 9 with An = 6.5 and on August 30 with An = 7.7. In general, with an increase in the energy level of earthquakes, the quantitative values of the dynamic parameters of foci increased within the range of measured values. Radiation friction Acr for the studied earthquakes had a negative value, indicating a complex slip of the gap in the source due to the uneven distribution of the strength properties of the deep environment in the source zones. In the whole energy range, the average value of the stress drop did not exceed Ao= 8-105 Pa (8 bar), and the apparent stresses n^ < 13 •lO5 Pa (13 bar).
Analysis of individual station definitions from the data for 2017 showed that the values of the seismic moments M0 and the moment magnitude Mw from the records of longitudinal waves P in most cases turned out to be significantly less than from the records of transverse S-waves: the mean value of the ratios k = M0S / M0 P for all stations, it turned out to be 2.5 with a maximum k =11 according to the «Alushta» station data for the earthquake focus of August 30 with a An = 7.7.
Keywords: amplitude spectrum, Brune model, spectral density, the angular frequency, seismic moment, radius of dislocation, stress drop, displacement along a discontinuity, radiation friction.
1. References
2. Svidlova V. A., Sykchina Z. N., Bondar' M. N., Bojko V. A. Sejsmichnost' Kryma v 2016 godu // Uchenye zapiski Krymskogo federal'nogo universiteta im. V. I. Vernadskogo. Geografiya. Geologiya. 2017, T. 3(69), no. 4, pp. 7-50 (in Russian).
3. Pustovitenko B. G., Merzhej E. A., Pustovitenko A. A. Kalinyuk I. V. Ochagovye parametry zemletryasenij Kryma 2013 goda // Sejsmologicheskij byulleten' Ukrainy za 2013 g. Sevastopol': NPC «EHKOSI Gidrofizika». 2014, pp. 12-21 (in Russian).
4. Pustovitenko B. G., Ehredzhepov E. E. Spektral'nye i dinamicheskie parametry zemletryasenij Kryma 2015 goda // Uchenye zapiski Tavricheskogo nacional'nogo universiteta im. V. I. Vernadskogo. Geografiya. Geologiya. 2016. T. 2(68), no. 4, 2016, pp. 53-68 (in Russian).
5. Pustovitenko B. G., Ehredzhepov E. E. Ochagovye parametry zemletryasenij Krymsko-Chernomorskogo regiona 2016 goda // Uchenye zapiski Tavricheskogo nacional'nogo universiteta im. V. I. Vernadskogo. Geografiya. Geologiya. 2017. T. 3(69), no. 4. pp. 51-69 (in Russian).
6. Brune I. V. Tectonic stress and the spectra of seismic shear waves from earthquakes // J. Geophys. Res. 1970. V. 75, no. 26, pp. 4997-5009 (in English).
7. Pustovitenko B. G., Panteleeva T. A. Spektral'nye i ochagovye parametry zemletryasenij Kryma. Kiev: Naukova dumka, 1990. 251 p. (in Russian).
8. Pustovitenko B. G., Kul'chickij V. E. Ob ehnergeticheskoj ocenke zemletryasenij Krymsko-Chernomorskogo regiona // Magnituda i ehnergeticheskaya klassifikaciya zemletryasenij. M. : IFZ AN SSSR. 1974. T. 2, pp. 113-125 (in Russian).
9. Aptekman Zh. Ya., Belavina Yu. F., Zaharova A. I., Zobin V. M., Kogan S. Ya., Korchagina O. A., Moskvina A. G., Polikarpova L. A., Chepkunas L. S. Spektry P-voln v zadache opredeleniya dinamicheskih parametrov ochagov zemletryasenij. Perekhod ot stancionnogo spektra k ochagovomu i raschet dinamicheskih parametrov ochaga // Vulkanologiya i sejsmologiya. 1989, no. 2, pp. 66-79 (in Russian).
82
10. Aptekman Zh. Ya., Daragan S. K., Dolgopolov V. V., Zaharova A. I., Zobin V. M., Kogan S. Ya., Korchagina O. A., Moskvina A. G., Polikarpova L. A., Chepkunas L. S. Spektry P-voln v zadache opredeleniya dinamicheskih parametrov ochagov zemletryasenij. Unifikaciya iskhodnyh dannyh i procedury rascheta amplitudnyh spektrov // Vulkanologiya i sejsmologiya. 1985, no. 2, pp. 60-70 (in Russian).
11. Kostrov B. V. Mekhanika ochaga tektonicheskogo zemletryaseniya. M.: Nauka, 1975. 179 p. (in Russian).
12. Kejlis-Borok V. I. Issledovanie istochnikov, priblizhenno ehkvivalentnyh ochagam zemletryasenij // Trudy Geofizicheskogo Instituta AN SSSR. 1959, no. 9(136), pp. 20-42 (in Russian).
13. Hanks T. C., Kanamori H. A. Moment magnitude scale // J. Geophys. Res. 1979, v. 84, no. 35, pp. 2348-2350 (in English).
14. Pustovitenko B. G., Pustovitenko A. A., Kapitanova S. A., Porechnova E. I. Prostranstvennye osobennosti ochagovyh parametrov zemletryasenij Kryma. Sejsmichnost' Severnoj Evrazii. Obninsk: GS RAN. 2008, pp. 238-242 (in Russian).
15. Pustovitenko B. G., Kul'chickij V. E., Kalinyuk I. V., Merzhej E. A. Sravnitel'nye harakteristiki akselerogramm, smodelirovannyh dlya raznyh sejmotektonicheskih uslovij // Sejsmostojkoe stroitel'stvo. Bezopasnost' sooruzhenij. 2012, no. 6, pp. 20-23 (in Russian).
16. Pustovitenko B. G. Merzhej E. A. K prognozu ozhidaemyh sejsmicheskih vozdejstvij ot sil'nyh zemletryasenij Kryma // Dopovidi Nacional'noi Akademii Nauk Ukraini. 2014, no. 5, pp. 102-109 (in Russian).
83
