Анализ инструментальных записей схода ледника Копка по данный покапьной сети сейсмических наблюдений
В.Б.ЗаапишвишН.И. Невская; РЖМаребое3
СХОД ЛЕДНИКА КОЛКА.
ОБРАБОТКА СЕЙСМИЧЕСКИХ ДАННЫХ
Локальная сеть сейсмических цифровых станций на территории Республики Северная Осетия-Алания была организована Геофизическим центром экспериментальной диагностики ВНЦ РАН (ГФЦЭД) в 1998 году. На всех пунктах наблюдения установлена регистрирующая аппаратура «Альфа-Геон». В 2002 году функционировало 6 сейсмических пунктов наблюдений: во Владикавказе, Ардоне, Чиколе, Заманкуле, Фиагдоне и Цее (соответственно коды станций VLD, ARD, CI К, ZAM, FIG, CEY). Кроме перечисленных сейсмических станций, следует отметить установку в 2003 г. станции в Дигории (DIG). Крайне необходимая сейсмическая станция в пункте Кармадон (KAR) была включена в сеть 8 сентября 2003 г, после почти двухлетнего бездействия. В настоящее время в сеть объединено 8 сейсмических станций [1 ].
Для обработки данных сейсмических наблюдений по методике Центра «ГЕОН» применяется комплекс компьютерных программ, позволяющих выполнять
Г_"'| '/üi-SC/r
*S. rn ,:.W 43 v22 Л Ii. НЧЮм
Cl". ¿WihAJ'l I «У' г*
I 44,-tiF гЛ Г,I ,4lM н
(31-37) i'l. HMhmMkffJ
Mjtlli
in ui.
С I '-1' L' Y. 111111
r Дпмщнй
I i
J.
j IJ t >L L.k
hoftii
кзшк MiiHia
EIil,iiiiik k.kiti3
Рис. 1. Расположение сейсмических станций, зарегистри-_ровавших сход ледника Колка._
локализацию очагов землетрясении, а также определять напряженно-деформированное состояние массива, т.е. механизмы очагов землетрясений.
20 сентября 2002 г. действующими на тот период станциями указанной локальной сети (Фиагдон, Заман-кул, Чикола, Владикавказ 31, Владикавказ 37, Ардон) был зарегистрирован сход ледника Колка (рис. 1-6). Обе станции Владикавказ были расположены на одном участке [2]. Необходимо отметить, что время срабатывания станций сети определялось особенностями настройки. Локальная сеть наблюдений создавалась в основном для регистрации сейсмических событий, которые, как правило, в регионе характеризуются значительно меньшей длительностью.
Необходимо отметить, что величина регистрируемого минимального уровня воздействия, учитывая процесс затухания сейсмической энергии, фактически является уровнем трнггерного включения регистраторов РСС «Альфа-Геон» станций сейсмических наблюдений локальной сети ГФЦЭД. Другими словами, сейсмические станции работают в ждущем режиме и включаются при достижении амплитудой движения грунта определенного уровня. При этом, очевидно, на чувствительность срабатывания станций заметное влияние может оказывать и спектральный состав приходящих колебаний (рис. 7).
Анализ амплитудных спектров Фурье показывает, что с удалением от места схода ледника или относительно «исходных» записей ст. Фиагдон спектральный состав колебаний становится низкочастотным (рис. 7). Особенно заметно поглощается высокочастотная составляющая в районе расположения ст. Чикола (рис. 5). Это соответствует хорошо известному факту прямой пропорциональности поглощения частоте или даже квадрату частоты колебаний. С другой стороны, на записи ст. Заманкул (рис. 6) наблюдается факт сдвига спектра в высокочастотную область спектра. Указанные особенности могут быть объяснены рядом причин, которые более корректно реализовать после получения большей информации.
В то же время на основе работы с фондовыми материалами было установлено, что грунтами ос-
' В.Б. Заалишвили - д.ф.-м.н., директор ГФЦЭД ВНЦ РАН и РСО-А.
2 НИ. Невская - зав. лаб. ГФЦЭД ВНЦ РАН и РСО-А.
3 А.К. Харебов - аспирант, ГФЦЭД ВНЦ РАН и РСО-А,
С-.Ф-141И.+ 4 11Л1 ли 11:Н 55 !+: 11:32
..... , ,.
1г
Л 1 л '
и
II к л 1-Я Л1 Г! С 311
Рис. 2. Запись схода ледника Колка станцией Фиагдон*
59
более близкой к месту расположению ледниковой массы до схода. Следовательно, в качестве исходных рассматриваются записи указанной станции. Очевидно, мы можем видеть только продолжение воздействия, которое включило станции (само запаздывание, очевидно, не превышает нескольких секунд). Ясно видно, что воздействие имеет отчетливую вертикальную составляющую. Но, при этом, амплитуда вертикальной составляющей раз в пять меньше горизонтальных. Другими словами, наблюдается преобладающее горизонтальное движение ледника в направлении N8 и Е\¥.
Таким образом, началу отсчета времени движения каменно-ледовой массы соответствует отметка (1). Необходимо отметить, что дальше на записях отсут-
Рис. 3. Запись схода ледника Колка станцией Владикавказ 37. нования ст.Заманкул являются конгломераты, перекрытые слоем элювия мощностью 2-3 м [3]. Необходимо уточнить грунтовые условия с помощью бурения скважин.
С целью корреляции инструментальных записей с процессом схода ледника были построены графики зависимости выделенной энергии от времени. Для более четкого выделения экстремумов было проведено сглаживание данных и фильтрация по высокочастотным составляющим (рис. 9-10).
Привязку к местности можно проводить по вектору наибольшего выделения энергии, а также по изменению угла между компонентами N8 и Е\\'. На всех записях направление преобладающего выделения энергии соответ ствует компоненте Север-Юг (N8), что подтверждает преобладающее движение ледовых масс (рис.. 9-10).
В соответствие со схемой транзита ледника Колка (рис. 8) первое отклонение или само включение сейсмических станций можно связать с обвалом или обвалами ледовых масс висячего ледника с г.Джимарай-хох на верхнюю часть ледника Колка (1), второе - с соударением ледово-каменной массы о правый склон транзита после прохождения ледника Майли (2), третье - с ударом ледовой массы по т.н. «Кармадонским воротам» (3).
Для изучения особенностей движения ледника Колка наиболее корректно сопоставление всех трех составляющих инструментальной трехкомпонент-ной записи ст. Фиагдон (рис. 2, начало записи), наи-
............. Л31*Ы*»11:Н:1! 4)1:11
II) 1Ц
Рис. 4. Запис ь схода ледника Колка станцией Владикавказ 31.
с- ......Jiil.HU 14.11 и К I! II
■1 1 ......«I
I И I и РШР1 .....*
ш Н) III
Рис. 5. Запись схода ледника Колка станцией Чикола.
С1.11-Н.РТГ-3«.«.»Ф1 1«; К-Р1 I? ч II
Рис. 6. Запись схода ледника Колка станцией Заманкул.
' Составляющие - 1Л) - вертикальная, N8 - Север-Юг, Е\У - Восток-Запад.
» С Н 4
№ 3
СТ. Ф И4ГЦ □ Н
£
с с
4
У ч Л ►»„'■ЧА1
л 1 у \ ■ м^КРСО
«Ч-—_ -л »-43ЕС-
10 1! Ч|Е 191». Г и
и
а
4 т. Впацнкжвка} 31
и 3
1 /
У '_¡1— у^_
Л" г' и л г я « VI ____
10 I & м
Ч.чг 1 йт А. Г 4
ствует какое-либо движение. Анализ особенностей различных составляющих записей на включившихся станциях и механизма процесса схода показывает, что это может быть обусловлено равномерным движением ледово-каменной массы (может быть, даже в виде скольжения), которое не регистрируется из-за отсутствия ускорения. При этом начальный пик, хорошо видный на вертикальной составляющей, есть именно начало схода. Таким образом, начало движения процесса схода ледника Колка есть ни что иное, как пик с отметкой (1).
Точку (2) - прохождение через ледник Майли с резким поворотом - мы выделяем на основе анализа записей горизонтальной составляющей т. к. вследствие удара ледниковых масс преобладающий вектор энергии в данном месте должен быть направлен на восток. При этом мы имеем торможение к источнику, но благодаря гравитационному ускорению равнодействующий вектор направлен вниз по его транзиту и руслу р. Геналдон.
Наибольшее выделение энергии наблюдается при торможении ледника в точке с отметкой (3). На рас. 9 ясно виден соответствующий пик на состав-
« т. Владикавказ 3 Т
Л 1 Н 1>11№
\ А 1 г ■ ни ¡¿ГГ
V ■ ни 1Е ЕС
ш
ч *е г i г я . Ги
ГТ ЧИЯПЛ*
\ в- а »1- ■
Л / Вы в л геТШ
кЛ'/«^
I и 14
Ч л г га и, Г и
.'II
с т. Эаадантьш
Л ч «•.!£+■• АЛ?
г № ) 1 г V Ч ■ лпв 1Ш
АЛ* 1/ к ь-ламДО
10 И
Члг 1 А. Г 4
2Ъ
Рис. 7. Амплитудный спектр Фурье инструментальных записей схода ледника Колка.
ляющейДалее наблюдается движение следующей за фронтальной частью ледово-каменной массы.
Далее оценим скорость движения ледниковой массы. Учитывая, что расстояние до Майлийского поворота от середины (середины инерционной массы)
■ if¿parte >* «■1» *з кто
.1^.1» Hki.i|>
Рис. 8. Схема схода ледника Колка.
сформированного ледника равно 4 110 м и время прохождения -80 с, средняя скорость на участке от метки 1 до метки 2 равна Vj -! 80 км/ч, и на участке от метки 2 до -3 при расстоянии, равном 13 450 м, и времени прохождения 150 с средняя скорость составит Г, =320 км/ч. В дальнейшем расстояние и соответственной скорости будут уточняться с уче-
том рекомендции геологов.
Исследование графиков показывает, что во время торможения «Кармадонскими воротами», где выделение энергии наибольшее, вначале произошел вертикальный удар (рис. 9, составляющая 1X0), очевидно, обусловленный падением передней части или фронта массы, а затем - горизонтальный N8 (основное соударение). Указанный вертикальный удар произошел на расстоянии около 400 м до главного удара. Это можно объяснить выполаживанием в этом месте ложа р. Ге-налдон и т.п. Далее, через период времени около 4 с произошел удар о «Кар-мадонские ворота» с преобладающим движением на север. На составляющей МВ ст. Фиагдон после основного удара отчетливо прослеживаются последующие волны. Достаточно большие амплитуды колебаний обусловлены: остановкой основной массы в районе «Кармадонских ворот», движением «хвоста» ледового потока, сходом грязе-селевого потока в направлении на север и т. п.
Для оценки времени прохождения ле-дово-каменной массы ледника по данным различных станций измеряли величину времени между соответствующими метками на исходной (Фиагдонской) записи и выделяли их на других записях. Следует отметить, что особенности выделенных метками точек надежно прослежи-
ваются. Такой анализ с конца записей, в частности, показал, что есть станции (Владикавказ 37, Чикола и Заманкул), где начало процесса практически утеряно. Это обусловлено поздним срабатыванием триггеров удаленных станций из-за малого амплитудного уровня и других причин. В то же время благодаря этому была получена возможность регистрации продолжения движения уже после выключения других вовремя включившихся станций. Особо необходимо отметить неудовлетворительно организованную службу времени на всех станциях локальной сети на период схода ледника. Это значительно усложнило обработку данных и повысило величину погрешности. Кроме того, сеть характеризовалась на период схода ледника неудовлетворительной подборкой параметров срабатывания триггеров, нередко приводящей к отсутствию на записях землетрясений первых вступлений и т. д. [4].
СОБЫТИЕ, ПРЕДВАРЯЮЩЕЕ СХОД ЛЕДНИКА
Анализ данных, зарегистрированных сетью станций ГФЦЭД, показывает, что за 12 минут до схода ледника имело место событие, зарегистрированное только двумя станциями: Чикола и Владикавказ 31 и одной из четырех сейсмических станций «Института Гидропроект», расположенных на расстоянии около 36-40 км к западу-юго-западу от Геналдонского ущелья.
ы». (г «ufri¿к..
j < 1
1 г г-,,-/4 \ Л л/ \/ Ч
1 хм 1 V и V
4D IEE II] ■ ршл. с
»: ж ::s
Рис. 9. Зависимость энергии от времени процесса схода, ст. Фиагдон.
ТО*1 4
ъ
■ lli:Lci.n.n
О")
э 1 — — —L-J МШ Jl "IV Б
— — =г ± Л £ =Н
*>
* л *
г а>
t *
о
О
liU, pfff
l!i' l+j № Ь»^^ г
A .Ж ДО
V
» Ш Ш
Рис. 10. Зависимость составляющей ТО-вектора выделения энергии
от времени.
При этом указанная станция АЦРСС № 11 располагалась на «узком мысу с крутыми обрывами к р. Ар-дон и р. Баддон на высоте 1624 м» [5] и, таким образом, была наиболее чувствительной [5].
Станция в Фиагдоне, расположенная (наиболее близко) в 18 км к северо-западу от Геналдонскош ущелья, не включилась на указанное событие. Обе включившиеся станции Г'ФЦЭД Чикола и Владикавказ 31, очевидно, находились достаточно близко от источника колебаний (рас. 1). На это указывает; в частности, сопоставимость значений амплитуд по всем составляющим на записях обеих станций. Такое бывает при почти сферическом источнике, когда горизонтальные составляющие, как правило более низкочастотные и, отсюда, менее поглощаемые средой, не успевают на относительно малом расстоянии от источника превысить амплитуду вертикальных колебаний.
Анализ спектрального состава зарегистрированных колебаний (рис. 12) показывает, что в спектре вертикальных колебаний ст. Чикола наблюдаются два почти равных по амплитуде пика, разнесенных на частотном диапазоне: один на частоте колебаний 1 Гц и второй - 8 Гц. При этом на спектрах обоих горизонтальных колебаний наблюдаются максимумы на частоте 1 и 5 Гц. В спектрах всех составляющих записей ст. Владикавказ 31 наблюдаются максимумы только на частоте 5 Гц. В то же время максимум спектра колебаний ст. АЦРСС № 11 «Института Гидропроект» составляет 0,81 Гц.
На основании вышеизложенного и других соображений это событие, на наш взгляд, является, скорее всего, технологическим взрывом. К сожалению, мы имеем только одну составляющую записи «Института Гидропроект». Это обстоятельство и отсутствие данных амплитудного уровня колебаний не позволяет установить предполагаемый эпицентр события. Мы надеемся, что вскоре, несмотря на существующие организационные сложности, эта задача будет решена. Строго говоря, нельзя также полностью исключить, что это событие было локальным землетрясением. Но для этого необходимо полностью проанализировать записи, полученные за указанный период сейсмическими станциями из сопредельных республике районов.
ПРЕДВАРИТЕЛЬНАЯ РЕКОНСТРУКЦИЯ СОБЫТИЯ СХОДА ЛЕДНИКА
Предварительная, на настоящее время, реконструкция события схода ледника на основе недостаточно полного материала выглядит следующим образом.
В течение последних 1-2 месяцев с г. Джимарай-хох могли иметь место неоднократные обвалы, обусловленные ослаблением вывет-релой части горных пород и несколькими местными землетрясениями с магнитудами М^3,9 и М=5,1, интенсивность которых составляла 3-4 балла в районе ледника Колка. Ледник быстро набирал критическую массу. Могли иметь место ряд сопутствующих факторов (обводнение под ледником, подогрев и т.п.). На это указывает почти полное отсутствие трения и высокая скорость схода. Очередное событие, зарегистрированное за 12 минут до начала процесса схода (взрыв/землетрясение) двумя станциями ГФЦЭД и одной - «Института Гидропроект», очевидно, могло раскачать тело висячего ледника. Произошел обвал. Ледник Колка, достигший критической массы под воздействием очередного обвала (именно он обусловил включение сейсмических станций локальной сети ГФЦЭД, но не был записан), был буквально сдвинут со своего ложа (за несколько секунд до метки 1 на рис. 9-10). Далее ледник двинулся (заскользил) под воздействием гравитационной силы вниз по не особенно крутой поверхности со скоростью 180 км/ч. После удара о поворот у ледника Майли (метка 2 на рис. 9-10) ледово-каменная масса понеслась вниз со средней скоростью 320 км/ч на водной (может быть - водно-воздушной) подушке. Иначе трудно объяснить такую скорость. Удар о скалы «Кармадонских ворот» выделен меткой 3
(рас. 9-10). За 400 м до указанного удара имел место вертикальный удар, очевидно, обусловленный частичным обрушением фронтальной части ледо-во-каменной массы при выходе на равнину.
Для большей глубины анализа и обоснованности выводов необходимо:
• получение данных из сопредельных республике районов;
• уточнение данных грунтовых условий расположения сейсмических станций ГФЦЭД;
• анализ всех данных в их взаимосвязи.
ВЫВОДЫ
1. 20 сентября 2002 г. действующими станциями локальной сети Геофизического центра экспериментальной диагностики (Фиагдон, Заманкул, Чикола, Владикавказ 31, Владикавказ 37, Ар дон) был зарегистрирован сход ледника Колка. Запись ст.Ардон, крайне неудачная, подлежит дальнейшей обработке и в работе не рассматривалась.
2. Получены уникальные данные цифровых записей схода ледника. Обработка и анализ данных показывает их соответствие основным этапам предполагаемого реконструированно- Рис. 11. го сценария.
3. Проведены работы по установлению грунтовых условий сейсмических станций локальной сети ГФЦЭД по фондовым материалам. В дальнейшем предполагается детальное изучение грунтовых усло-
Инструментальная запись предшествующего события на станциях ГФЦЭД.
вий с помощью бурения скважин. Это позволит уточнить скоростную модель, принятую ранее для расчетов параметров землетрясений, и влияние грунтовых условий на волновые формы записей.
о ю 20 м « 50 еи та го
Чтила ?[ 1-3 11:11 Ни - 1Ь:&Т:П1
0 10 Ж ЗС а:. № 50 та .Ц| Са-кн <■. Е
ет. Вшаниш: 21 ЗКИ ЛИТ 13:31:]!! - 13:1Г:ЗЬ
а-
1
1 Е
Е
5-<
ст. ЧикАпа
и |№И
и
/и
■ 0 15
ЧА.етлгя, Гц
20 25
ет. Владикавказ ц
1 1
J 1
> у * Ш1>М
Ш 15
Ча ггл г.1 Гц
20 25
Рис. 12. Амплитудные спектры Фурье предшествующего события.
ЮМ 4
№ 3
ч ее*
J»
i ujk
4. С целью локализации событий с юга предполагается установить сейсмическую станцию в с. Казбеги. Чтобы оконтурить Казбекский вулканический центр, необходима установка сейсмических станций южнее и восточнее г. Казбек.
5. Необходимо получить записи событий схода ледника из сопредельных республике районов и государств. Это позволит значительно уточнить сценарий происшедшего события.
6. Необходимо ускорить обработку инструментальных данных, сведенных в специальную базу данных за период, предшествующий и последующий катастрофическому сходу ледника.
7. На основе изучения и анализа дополнительных данных с сопредельных территорий уточнить наличие-отсутствие технологических взрывов за указанный период на территории республики и за ее пределами и возможность их локализации.
8. Построить механизмы очагов наиболее представительных событий в исследуемый период с помощью адаптации уже полученных передовых алгоритмов и программ, и реализуемых при помощи коллег из центральных научных учреждений РАН.
9. Особо следует отметить неудовлетворительно организованную службу времени на всех станциях локальной сети ГФЦЭД на период схода ледника. Это значительно усложнило обработку данных и повысило величину погрешности. Необходимо срочно улучшить ситуацию.
10. При действенной государственной поддержке на основе результатов исследования будут разработаны эталонные сценарии опасных проявлений рутинного функционирования природного процесса в районе Казбекского вулканического центра при совокупном воздействии внешних и внутренних факторов, а также в других районах республики.
Кпнппнннтл |Ю, er. JiihMia
KnunDHiHTj U^ci.Diii#iiaiii]'31
Нешпонент? вд.рт Чнндоэ
Ü-JÖ
ч JP- ».в
i
к в*
■и
G.Z
П.'5
4il
л*. Ч н н*м
40
КОЯГ'ОНйчТй 0Y. ei.
П III
Güt □ ДО
> ft >
•7. Ii LH /V 1
щ £|Д1 yv v- и Щ
L'.'jD - V - -
Ii
M
bfcjüfja l
Ь-fll'NH К
Рис.13. Записи колебаний предшествующего события.
11. Анализ вышеизложенного показывает, что исследования необходимо продолжить. В частности, необходимо провести полевые исследования в районе предполагаемых эпицентров сейсмических событий.
Литература
1. Заалишвили В.Б., Невская Н.И., Канунов Э.А., Ходов В.Э., Громова H.A. Развитие существующей локальной сети сейсмических наблюдений на территории РСО-А и совершенствование обработки данных // Отчет Геофизического центра экспериментальной диагностики, 2003. 37с.
2. Заалишвили В.Б., Невская Н.И. Сбор и анализ сей-смогралт землетрясений (с определением их эпицентров), произошедших в сентябре и до сентября 2002г.// Комплексные геолого-геофизические исследования в Кар-мадонском ущелье Республики Северная Осетия Алания /июль-август-сентябрь 2003г. Док. 41 с.
3. Особенности геологического строения и гидрогеологических условий участков расположения пунктов контрольного сейсмического наблюдения (ПКН) на территории РСО-А // Отчет Геофизического центра экспериментальной диагностики, 2003. 54с.
4. Протокол стажировки Невской Н.И. в ГС РАН. Центральная опытно-методическая экспедиция геофизической службы РАН, 2002. Док. 4 с.
5. Сейсмические записи движения Кармадонского ледника. ЦСГЕНЭО-филиал АО «Институт Гидропроект» // О:ивег\ОЪ]ес1Хагата^\ПолМат092002Карл1 JIedHUK.doc. Док. 6 с.