СИЛЬНЫЕ ДВИЖЕНИЯ В ЭПИЦЕНТРЕ ЧУЙСКОГО ЗЕМЛЕТРЯСЕНИЯ КАК РЕЗУЛЬТАТ ВЫХОДА УДАРНОЙ ВОЛНЫ НА ПОВЕРХНОСТЬ ЗЕМЛИ

В.В. Кузнецов

ИКИР ДВО РАН, п.Паратунка, Камчатская обл., vvk@ikir.kamchatka.ru


Приводятся обобщенные результаты предварительного изучения сильных движений Чуйского землетрясения, полученные двумя экспедициями, организованными Институтом физики Земли и Горно-Алтайским университетом, во время которых изучены особенности сейсмотектонических разрывов, обвалов и оползней, исследованы случаи вспучивания и проседания почвы, а так же данные по излиянию воды и грязи на поверхность и образованию грифонов. Предлагается объяснение особенностей сильных движений грунта в рамках ударно-волновой модели землетрясения. С этой целью представлена информация о поведении твердого вещества при выходе на его поверхность ударной волны, о возникновении при этом волны разрежения и взаимодействие этих волн, результатом которого является образование разрывов на поверхности, вспучивание поверхности твердого тела, возникновение сложной картины сильных движений поверхностного слоя за счет взаимодействия возникших волн. Получены оценки эффектов, обсуждаются пути и перспективы изучения уникального природного явления - очага Чуйского землетрясения.


Чуйское землетрясение на Горном Алтае (М = 7.5, 03/09/27)

На юге Горного Алтая на территории Кош-Агачского района 27 сентября 2003 г. в 11час. 33 мин. по Гринвичу произошло сильное землетрясение с магнитудой 7.5 (по оценке ГС СОРАН и 7.3 - по оценке ГС РАН), которое произвело значительные разрушения в селах Бельтир (8 баллов), Чаган-Узун, Ортолык, Курай, Акташ (6-7 баллов). Очаг землетрясения располагался на территории Северо-Чуйского хребта, Чуйской и Курайской межгорных впадин и разделяющего их Чаган-Узунского приподнятого блока.

Землетрясение сопровождалось многочисленными афтершоками. Наиболее сильные из них два: первый - 27.09.03 в 18 час. 52 мин. по Гринвичу с магнитудой М = 6.4 и второй - 1.10.03 в 1ч. 03 мин.28 сек. по Гринвичу с магнитудой 7.0. Главный удар афтершока 1 октября произошел в районе пос. Акташ, где отмечались повреждения зданий 7-ми, 8-ми бальной интенсивности.

Магнитуды, координаты основного толчка и двух афтершоков, а так же глубины гипоцентров по данным Международной Сейсмологической сети USGS, National Earthquake Information Center, соответственно: 1) 2003/09/27 11:33 M 7.5 Z= 18 км 50.06N 87.73E; 2) 2003/09/27 18:52 M 6.3 Z= 33 км 50.08N 87.81E; 3) 2003/10/01 01:03 M 6.7 Z= 10 км 50.19N 87.67E. Наибольший интерес для нас представляют точные параметры Главного удара. По предварительным данным Геофизической службы СО РАН, его координаты: 49.8N 88.3E, что значительно отличается от данных международной сети.

Согласно уточненным данным, координаты гипоцентра Главного удара: 50.09N 87.98E, глубина 18 км. Так же уточнены глубины гипоцентров сильных афтершоков, которые не превышают глубины основного толчка, т.е. меньше или равны 18 км.

Область первичных, сильных разрушений поверхности Земли (область выхода очага на поверхность) можно изобразить в плане как прямоугольник со сторонами 5 20 км2, причем, его 20-ти километровая сторона проецируется на разлом Кускуннур. Зная координаты гипоцентра и размер области афтершоков, очаг землетрясения можно представить в виде параллелепипеда размером 5 20 18 км3, длинная грань которого ориентирована вдоль разлома Кускуннур, а его боковые поверхности наклонены к поверхности земли на угол примерно 65o в юго-западном направлении. По-видимому, главный удар Чуйского землетрясения пришелся именно под этим углом к поверхности. Вероятно, два главных афтершока можно рассматривать как землетрясения, которые внесли свой "вклад" в сильные движения грунта. К сожалению, оценить направления главных ударов этих землетрясений с желаемой точностью, вероятнее всего, не удастся. Судя по характеру афтершоков и особенностям разрушений (со слов очевидцев) направления главных ударов этих землетрясений отличались от главного толчка. Их наиболее вероятные направления: второго, - на северо-восток, в направлении с. Чаган-Узун; третьего, - на северо-запад, вдоль по разлому Кускуннур, в направлении альплагеря Актру (см. рис. 1).


Рис. 1. Карта распространения сейсмодислокаций [1]: 1 - сейсмотектонические разрывы; 2 - обвалы и осыпи; 3 - разжижения грунта; 4 - оползни; 5 - сейсмогравитационные трещины отрыва и оседания; 6 - разломы; 7 - озерные и ледниковые отложения.


В верхнем правом углу схема: стрелками отмечено направление главного удара землетрясения и область очага, обозначенная черной линией, показывающей расположение гипоцентров афтершоков. Прерывистая линия показывает разлом Кускуннур.


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

Участники экспедиции в эпицентральную зону Чуйского землетрясения от Горно-Алтайского университета (время работы - первые числа октября) обнаружили трещины уже в с. Курай, а также по дороге на Актру. В Курае трещина шириной до 10 см и длиной порядка 100 м образовалась в районе болота на окраине села. По дороге на Актру трещины наблюдались обычно вдоль крутых откосов, ширина их была 1-2 см, длина не превышала 5 м. Похожие трещины, но большей протяженности наблюдались на Чуйском тракте в районе Чаган-Узуна. На въезде в с. Бельтир в результате оползня было порвано дорожное полотно (этот случай широко известен благодаря снимкам в прессе). Одна из трещин прошла через с. Бельтир, вдоль ул. Москанова, значительно повредив дома, которые на нее попали. Ширина трещины около 1 м, глубина проседания грунта - от 0.5 м до 1 м. Были обследованы несколько трещин протяженностью порядка 1 км на плато Ада-Кель. Со слов местных жителей известно, что в эпицентре землетрясения (урочище Узук) образовалось две большие трещины, идущие в направлении Карагема и Кускуннура.

Наиболее подробно следы сильных движений (макросейсмичность) изучены совместной группой Института физики Земли и Алтайского университета. Результаты предварительного изучения сильных движений Чуйского землетрясения 2003 представлены в [1]. В дальнейшем, в этом разделе будем в основном цитировать эту работу, иногда привлекая неопубликованные данные очевидцев. Нарушения рельефа, по оценке авторов [1], образуют овал длиной 70 км и шириной около 15 км. Этот овал вытянут вдоль простирания разрыва (СЗ-ЮВ) и хорошо коррелирует с ориентировкой длиной оси овального облака эпицентров афтершоков. Концентрация вторичных сейсмодислокаций меняется, затухая по мере удаления от линейной вытянутой области выхода очага на земную поверхность. Обсуждая итоги предварительных исследований сильных движений Чуйского землетрясения, будем следовать авторам [1], которые выделяют четыре вида нарушений: 1) сейсмотектонические разрывы, 2) обвалы и оползни, 3) вспучивание и проседание почвы и 4) излияние воды и грязи на поверхность и образование грифонов.

Сейсмотектонические разрывы, вспучивание и проседание почвы.

В ходе вспарывания земной коры, очаг землетрясения вышел на поверхность в виде системы первичных сейсморазрывов, прослеженной на протяжении 20 км в бассейне р. Чаган-Узун на восточном погружении Северо-Чуйского хребта (рис. 2). Первичный характер разрывов характеризуется их морфологией и то, что эти трещины представляют собой систему параллельных рвов, секущих различные формы рельефа. Эта система вытянута в северо-западном направлении на протяжении 20 км.

Сейсмотектонические разрывы представляют собой эшелонированную систему трещин сжатия и расширения, укладывающуюся в линию СЗ простирания. Такая структура сейсмотектонического разрыва характеризует подвижку в очаге как сдвиг в горизонтальной плоскости. Строение зоны сейсморазрывов несколько отличается для различных участков. Так, например, на водоразделе рек Чаган и Елангаш разрыв разветвляется на две основных трещины, которые затем смыкаются, представляясь трещинами СЗ ориентировки (рис. 2-а). По западной ветви наблюдается правостороннее смещение на 0.5 м, по восточной, - смещение достигает 0.2 м. Блок, заключенный между этими рвами, опущен на 0.2-0.3 м и рассечен многочисленными диагональными трещинами. Поперечные разрывы перемежаются вспучиваниями и наползанием почвы высотой 0.3-0.7 м. Картина разрушений в долине реки Талдуры несколько отличается. Здесь алювиальные отложения рассечены рвами шириной до 2 м и длиной до 50 м. и закрытыми трещинами с амплитудой правого сдвига 0.1-0.2 м. Трещины сопровождаются линейными зонами с вздыбленной на высоту до 0.5 м почвой. Наиболее впечатляющие сейсмотектонические разрывы деформировали земную поверхность в широкой седловине на водоразделе рек Талдура-Кускуннур. Ширина седловины составляет около 4 км, сложена она отложениями морены среднего плейстоцена. Зона сейсморазрывов здесь имеет субширотное простирание. Сами сейсморазрывы выстраиваются в эшелонированную систему зияющих рвов, смыкающихся линейными зонами сжатия (рис 2-б). Сейсморвы в этой районе принимают колоссальные масштабы, характеризуясь шириной до 10 м, длиной до 300 м и глубиной более 30 м. По отдельным зияющим разрывам амплитуда правого сдвига составляет 0.5 м. Зоны сжатия здесь представляют собой валы вспучивания высотой до 2 м и длиной до 50 м и надвиги, сдваивающие слой почвы. Сейсморазрывы разветвляются на две ветви и зажимают между собой просевшие блоки (микрограбены). На восточном борту седловины ширина такого микрограбена составляет не более 50 м, а амплитуда опускания около 0.5 м. На западном склоне просевший блок имеет ширину до 500 м и опущен с амплитудой до 2.5 м. Авторы отмечают, что в зоне эпицентра Чуйского землетрясения находятся следы предыдущих землетрясений. Землетрясения, происходившие ранее, так же как современные, были следствием разрядки одного и того же очага. Чуйское землетрясение, по всей видимости, не является уникальным, - оно лишь несколько "обновило" следы прошлых, возможно, более сильных землетрясений. Этот факт однозначно демонстрирует, что очаг сильного землетрясения оказывается "привязанным" к определенной геологической структуре.

Обвалы и оползни.

Одно из наиболее эффектных и впечатляющих последствий Чуйского землетрясения (по мнению Гвоздарева) находится в урочище Арха-Узюк, где произошло нарушение рыхлого осадочного чехла на склоне, ослабленном вечной мерзлотой. После землетрясения в этом районе сошел оползень шириной около 1 км, вся его поверхность при этом была разбита сетью трещин шириной 1-3 м и глубиной до 10 м и более. Рядом с оползнем образовались ряд трещин сходных размеров. Несколько трещин прошло вдоль по хребту. В дальнейшем возможно сползание блоков, на которые эти трещины разбили склон. Само название места (в переводе - "оборванный лес") говорит о том, что подобное явление там наблюдается не в первый раз. По-видимому, самый масштабный след землетрясения - крупный оползень на правом берегу реки Талдура, приурочен к зоне выхода очага на поверхность, причем стенка отрыва оползневых масс точно ложится в зону распространения первичных сейсморазрывов, один из которых, по мнению авторов, мог служить "заколом" на поверхности склона, инициировав при этом, сход оползня. Оползень отсел по стенке отрыва высотой около 150 м и продвинулся в долину на расстояние не меньше 100 м. Объем оползневой массы по приблизительным подсчетам составляет около 30 млн. м3. Ширина оползневого тела более 500, длина - около 700 м. Подобные оползни меньших масштабов образовались на склонах долин рек Чаган и Чаган-Узун.


Рис. 2. Слева: (a) структура сейсморазрыва на водоразделе рек Чаган и Елангаш: 1 - закрытые трещины и разрывы, 2 - валы вспучивания, 3 - сейсморвы, 4 - амплитуда (м) и направления смещения. Большая стрелка - направление главного удара, параллельные стрелки - направления подвижек почвы. Справа вверху: (b) структура сейсморазрыва на водоразделе р. Талдура-Кускуннур. 1 - сейсморвы; 2 - бровки валов вспучивания; 3 - контуры валов вспучивания; 4 - оси валов вспучивания. Справа внизу: (c, d) разрезы по А-А и В-В, звездочка - гипоцентр, стрелка - направление главного удара, а - фронт УВ.


Основной толчок и самые сильные афтершоки землетрясения сопровождались массовыми обвалами и осыпными явлениями. Часто обвалы наблюдались на склонах, с углом наклона, менее угла естественного откоса (25o ). В отложениях осыпей долины р. Актру наблюдается аномальное распределение обломков разной величины: в верхней части большинства осыпных конусов сконцентрированы наиболее крупные обломки, что свидетельствует о развитии импульсных обвальных явлений.

Грифоны.

Одним из удивительных явлений, сопутствующих этому землетрясению являются выбросы и разливы разжиженных песка и грязи в виде грифонов или грязевых вулканчиков. Приурочены они в основном к низменным водонасыщенным рыхлым породам (глины, суглинки, пески), способным при динамических воздействиях, в частности при землетрясениях 7 и более баллов, "разжижаться" до текучего состояния и изливаться на поверхность. При этом образуются грифоны в виде конических или песчанно-глинистых "вулканчиков" воронкообразной формы. По оценке, приводимой Гвоздаревым, объем воды, выброшенной каждым из обследованных ими грифонов, составляет около 100 куб. м. Во время экспедиции, этой группой, только с машины, на дороге Курай - Актру было замечено 8 таких грифонов, хотя их было, конечно, значительно больше. Авторы [1] описывают грифоны, в которых наблюдается выброс глины на поверхность в виде слоя толщиной не менее одного метра и занимаемой площадью около 50 м2. Объем песчанно-глинистой грязевой массы такого грифона, без учета ушедшей воды, значительно больше, он превышает 1000 куб. м. Количество грифонов в эпицентральной зоне никто не подсчитывал, по общему впечатлению, грифонов очень много. Известно, что грифоны бывают не только грязевые, а так же пылевые и песочно-пылевые "вулканчики", которые неоднократно наблюдались во время Ташкентского землетрясения 1966 г.

Отметим важную геологическую особенность района, которую необходимо учитывать при построении модели сильных движений, и которая состоит в том, что сильные движения в основном происходят в слое грунта толщиной не более 20-30м, а нижней границей слоя является верхняя граница вечной мерзлоты. На это обратил внимание автора А.Ю.Гвоздарев. Если наличие слоя вечной мерзлоты в районе эпицентра Чуйского землетрясения подтвердится, - это должно значительно отличать сильные движения Алтайских землетрясений, например, от макросейсмичности Калифорнийских.


Общие закономерности при выходе ударной волны на поверхность Земли

Основные соотношения для ударных волн (УВ) [2-4]:

Закон сохранения массы:

ρoD = ρn(D - un) или un +εD = 0,        (1)

где ρn - плотность вещества на фронте волны, ρo - невозмущенная плотность, D - скорость распространения фронта волны, un - скорость течения среды за фронтом ударной волны.

Закон сохранения импульса:

σn = pn - po = ρounD,        (2)

где po - начальное давление среды (до прохождения ударной волны), pn - давление на фронте ударной волны.

Уравнение сохранения энергии имеет классический вид:

En - Eo = 1/2 ( pn · pn) (vo - vn),        (3)

где Eo и En - внутренняя энергия до фронта ударной волны и за её фронтом, vo и vn, соответственно, - удельные объемы (v = 1/ ρ).

В неподвижной среде (uo = 0) соотношения скоростей на фронте ударной волны представляются в виде:

D = [ρ/ρo (p - po) / (ρ - ρo )]1/2,
u = [1/ρ ρo (p - po) (ρ - ρo )]1/2.         (4)

Уравнения (4) справедливы для сухого грунта. Если грунт водонасыщенный, то D и u, после подстановки уравнения состояния водонасыщенного грунта в виде [2]:

ρ = ρo {a1(p/po)-1/k1 + a2 [k2(p - po)/ρ2c22 + 1]-1/k2 + a3 [k3(p - po)/ρ3c23 + 1]-1/k3}-1

где a1,a2,a3 - содержание компонент (по объему): воздуха, воды и твердого грунта, соответственно, (a1 + a2 + a3 = 0), k1, k2, k3 - показатели изоэнтропы соответствующих компонентов; будут определяться следующим образом:

D2 = (p - po)/ρo {1 - a1(p/po)-1/k1 - a2 [k2(p - po)/ρ2c22 + 1]-1/k2 - a3 [k3(p - po)/ ρ3c23 + 1]-1/k3}-1;


u2 = (p - po)/ ρo {1 - a1(p/po)-1/k1 - a2 [k2(p - po)/ ρ2c22 + 1]-1/k2 - a3 [k3(p - po)/ ρ3c23 + 1]-1/k3}.         (5)

Пористость грунта, иначе, суммарная концентрация воздуха и воды ( a1 + a2) являются определяющими при оценке величин скоростей D и u. Увеличение концентрации защемленного воздуха снижает величину скорости ударной волны D и приводит к увеличению массовой скорости течения грунта u.

Экспериментально неоднократно отмечалось, что при прохождении УВ через среду происходит "размыв" фронта (толщиной l). Принято считать, что этот эффект связан с переходом среды в квазижидкое состояние, характеризующееся наличием у среды динамической вязкости [5]. Например, в преcсованном порошке соли NaCl, при прохождении ударной волны (u = 103 м/с, l = 10-3 м, τ = 10-6 с) размыв фронта соответствовал вязкости = 104 - 105 Пз. Вязкость сплошного материала каменной соли имеет порядок 1014 - 1015 Пз, что на десять порядков меньше вязкости, наблюдаемой в ударных экспериментах. Можно оценить величину вязкости по параметрам колебательного неустойчивого характера сильных движений почвы при Нортриджском землетрясении 1994 г. При этом вертикальный размер пространства, на котором наблюдался колебательный процесс в эпицентре землетрясения: а ≈ 1 м, величина напряжения в главном толчке землетрясения: σ ≈ 100 бар, а скорость колебаний частиц грунта v достигала 1 м/с. Знание этих параметров позволяет оценить величину вязкости "жидкой" среды: = σ a/v ≈ 105 - 106 Пз. Полученные оценки вязкости показывает, что вещество грунта в момент землетрясения в его эпицентре находится в "жидком" состоянии. Это обстоятельство позволяет воспользоваться формулами Эйлера, Навье-Стокса и Пуазейля для течения жидкости, соответственно:

dv/dt = F - 1/ ρ grad p,               (6)
dv/dt = F - 1/ ρ grad p + χdv,    (7)
v = ∫ u dS,                                  (8)
     S

где F - внешняя сила, χ=/ρ - кинематическая вязкость.

Скорости волн. Воспользуемся подходом, развитым для оценки скоростей перемещений частиц упругой среды под действием силы, приложенной в очаге землетрясения. Этот подход был предложен Кейлис-Бороком и развит Матуямой ([6] с. 40). Скорость перемещения в горизонтальной плоскости:

u = 1/4 ρ [2/x2 (φ - Ψ) + 2Ψ],     (9)

где φ и Ψ сферические волны вида: φ = (1/r) F [t - (r/Vp)]; Ψ = (1/r) F [t - (r/Vs)]. Эти волны распространяются на расстояние r со скоростями Vp и Vs, соответственно [6]. Функция F(t) связана с силой f(t): f(t) = d2F(t)/dt2. Развитый Матуямой подход имеет значение для нашей модели в том плане, что речь идет не прямо о Р- и S-волнах, а о волнах, имеющих скорости, близкие к скоростям распространения Р- и S-волн.

По-видимому, некоторые авторы допускают ошибку, утверждая, что воздействие на грунт именно S-волны - приводит к разрушениям в гипоцентре землетрясения. Согласно представлений Матуямы, правильнее говорить о волнах, двигающихся со скоростями, близкими (или равными) скоростям Р- и S-волн. Как известно из физики ударных волн, сильная ударная волна распространяется со скоростью D, близкой к скорости Р- волны. Соответственно, чем слабее ударная волна, тем меньше её скорость распространения, и, в пределе, она может быть равной скорости распространения S-волны: D → Vs.


Особенности распространения ударных волн в грунтах.

Приведем некоторые данные по особенностям распространения ударных (взрывных) волн в грунтах [2, 5]. При камуфлетном (подземном) взрыве наблюдается так называемое геометрическое подобие (скейлинг) по корню кубическому из энергии взрыва Е [5]. Данные по многим параметрам движения среды при взрыве удается свести к универсальным зависимостям, если использовать в качестве линейного масштаба комбинацию (E/ρc2)1/3. Большое количество проведенных экспериментов по изучению давления в ударной волне от расстояния и количества взрывчатого вещества подтверждает зависимость [2]:

p = po(C1/3/R),      (10)

где po - начальное давление в атм., R - расстояние от взрыва то точки измерения в метрах, C - вес взрывчатого вещества в кг.

Для широкого набора различных грунтов, как водонасыщенных, так и сухих, используется следующая формула:

p = K(C1/3/R) µ ,     (11)

здесь K и µ - коэффициенты, величины которых приведены в Таблице 1.

Для того чтобы перейти от веса взрывчатого вещества к энергии и применить формулу (10) в нашей ситуации, воспользуемся соотношением между энергией и весом тротила. Известно, что удельная кинетическая энергия тротила равна 4 кДж/г. Положим, что сейсмическая энергия и энергии взрывчатых веществ эквивалентны. Энергия землетрясения с магнитудой М = 7.5 равна 1023 эрг, (1013 кДж), что эквивалентно 1012 г (одна мегатонна) тротила, тогда параметр C1/3 ≈ 103 кг-1/3. Давление р в очаге (на фронте ударной волны) Чуйского землетрясения (М = 7.5), по-видимому, может достигать 1 кбар. (Для сравнения, давление на фронте УВ в момент Нортриджского землетрясения М = 6.7, Е = 1022 эрг, порядка 100 бар). Давление на фронте ударной волны (волны сжатия) на расстоянии 650 км от очага, в Новосибирске, если принять самые высокие значения коэффициента К из Таблицы 1, оказывается порядка 1 бара. Очевидно, что эта оценка, по-видимому, наиболее высокая из всех возможных. Если, к примеру, параметр К будет меньше, а µ - больше, то давление р будет меньше чем 1 бар (1 атмосфера) на несколько порядков.

Величина смещения грунта во время главного толчка Чуйского землетрясения (по данным, полученным А.Г.Филиной на Новосибирской сейсмостанции), составляет примерно 50 микрон. Зная величину смещения, можно оценить деформацию ε, производимую приходящей волной: ε ≈ (1 - 5)10-3. Согласно [2] величина напряжения σ водонасыщенного грунта, соответствующая ε = 0.1, составляет примерно 500 бар и эта величина возрастает линейно с ростом ε: σ ε. Эта закономерность дает основание оценить σ = р ≈ 5 бар при ε ≈ 10-3. Аналогичная зависимость (σ ε) для вакуумированного песка приведена в [Н], откуда следуют примерно те же цифры для величины σ. Скорость звука в таком грунте ≈ 400 м/с. По-видимому, реальная величина σ меньше, однако эти оценки показывают (при различных способах их получения) один порядок величин напряжения в волне. Скорость распространения волны, амплитуда которой составляет 5 бар, от очага Алтайского землетрясения до сейсмостанции Новосибирск, расположенной примерно в 650 км от очага равна 3.6 км/с. Это намного больше, чем местная (в грунте) скорость звука. Возникает вопрос, что это за волна (ударная или нет) приходит к сейсмостанции, расположенной так далеко (650 км) от очага землетрясения? Согласно представлениям, развиваемым в [2], волна такой амплитуды ударной быть не должна. Однако, если производная d2σ/dε2 > 0, волну можно считать ударной вне зависимости от величины σ [Л стр. 52]. Если условие d2σ/dε2 > 0 не выполняется, - волна считается плоской волной сжатия, которая распространяется в упруго-пластической среде.

Таблица 1


Характеристика грунта К µ
Водонасыщенный песок с α= 0 600 1.05
Водонасыщенный песок с α= 5 10-4 450 1.5
Водонасыщенный песок с α= 10-2 250 2.0
Водонасыщенный песок с α= 4 10-2 45 2.5
Неводонасыщенный песок с ρ= 1.6 - 1.7 г/см3 15 2.8
Неводонасыщенный песок с ρ= 1.5 - 1.6 г/см3 7.5 3.0
Неводонасыщенный песок с ρ= 1.4 - 1.5 г/см3 2.5 3.5

Явления на поверхности Земли, связанные с выходом на неё ударной волны.

Акцентируем внимание на процессах, сопровождающих выход не слишком сильной ударной волны на поверхность (см. рис. 3). Выход наклонной к линии горизонта ударной волны изображен на рис. 3-1. Здесь фронт ударной волны, шириной а, двигающийся со скоростью со, подходит под углом α. к земной поверхности, вещество которой обладает плотностью ρ. Распространение ударной волны сопровождается сравнительно медленным движением сплошной среды с массовой скоростью u. Ударная волна отражается на участке АВ и уходит вниз со скоростью со. При отражении возникает быстро затухающая волна разгрузки (разрежения, растяжения), движение которой сопровождается перемещением вещества со скоростью uо. Действие на грунт силы, связанной с перемещением массы с вертикальными (вверх) компонентами этих скоростей приводят к разуплотнению, расширению грунта и подъему его относительно начального уровня на высоту порядка метра. На участке АВ может возникнуть система сил, направленных друг навстречу другу соответственно тому, как направлены горизонтальные компоненты массовых скоростей: u - uo (рис. 3-4), либо u' - uo' (рис. 3-6).

На рис. 3-2 изображено распределение плотности до (слева) и после выхода нормально падающей ударной волны (справа). Видно, что после выхода ударной волны и отражения её от поверхности, вещество с плотностью ρ1 (меньшей ρo) находится выше начального уровня. Аналогичная картина представлена на рис. 3-3, где изображены скорости в нормально падающей ударной волне со скоростью co и массовые скорости u и u' ≈ 2u до выхода ударной волны на поверхность (слева) и после её выхода на поверхность. На рисунках 3-5 и 3-7 изображены разломные диаграммы, соответствующие ситуациям, изображенным на рисунках 3-4 и 3-6. Темные секторы диаграмм показывают области сжатия, светлые - растяжения.

По данным сейсмологической сети ИОЦ ГС РАН определен механизм основного толчка Чуйского землетрясения, который практически совпадает с диаграммой, изображенной на рис. 3-5, что может говорить в пользу предлагаемой модели.


Рис. 3. Выход ударной волны (шириной а) на поверхность Земли под углом α к линии горизонта. Пояснения в тексте.



Угол выхода ударной волны на поверхность

Гипоцентр Чуйского землетрясения по уточненным данным расположен на глубине 18 км. "Ось зоны активности" афтершоков практически совпадает с разломом Кускуннур. Область первичных разрушений отстоит от этого разлома примерно на 5-7 км. Это может означать, что главный удар землетрясения (его очаг) выходит на поверхность под углом порядка 65 градусов. Рассмотрим, к каким выводам можно прийти, если известна величина этого угла.

Как следует из ударно-волновой модели, в сейсмически активной, напряженной геологической среде кроме литостатической нагрузки должна быть еще одна (или несколько) других напряжений, которые, взаимодействуя со средой, приводят её в некоторое, названное нами, - когерентное состояние. Если среда находится в таком состоянии, то в ней может развиться процесс лавинного образования трещин и генерации при этом акустических волн, которые, взаимодействуя с трещинами, помогают им раскрыться в едином процессе. В результате, в среде может возникнуть разрыв характеристик среды (например, плотности) и образоваться ударная волна. Угол выхода ударной волны на поверхность не случаен. Он определяется рядом обстоятельств, связанных с характером нагрузок и реологией среды: Во-первых, это литостатическая нагрузка, величина которой пропорциональна глубине. Она приводит к возникновению трещин, углы раскрытия которых относительно вектора напряжений, направлены, согласно закона Кулона-Мора (см. рис. 4-а). В среде, напряженной только литостатическим давлением, углы раскрытия трещин сильно различаются по глубине. В такой ситуации когерентной структуры возникнуть не может. К литостатической нагрузке "добавим" напряжение, возникающее в результате приложения силы, действующей в горизонтальном направлении. Если распределение нагрузки окажется таким, что в результате сложения литостатической и боковой, углы раскрытия трещин окажутся одинаковыми и направленными примерно под углом в 45 градусов, то и ударная волна подойдет к поверхности Земли снизу под этим углом к горизонту (рис. 4-b). Этот угол в сейсмологии принято считать углом выхода очага на поверхность Земли. Представим себе, что в формировании "угла когерентности" принимает участие напряжение, действующее под углом к поверхности Земли. Это напряжение может быть, связано, например, с влиянием разломов (рис. 4-с). В этом случае угол выхода ударной волны оказывается больше чем 45o и приближается к углу 90o. На рис. 4-d изображена ситуация, в которой когерентность среды достигается путем небольшой добавки, обязанной, например, действию атмосферных возмущений, или влиянию лунных приливов и пр. Эту ситуацию принято считать триггерным эффектом.


Рис. 4. Изменение угла образования трещин при вариации характера нагрузки на среду: а - вертикальная, литостатическая; b - литостатическая + боковая; c - литостатическая + "под углом"; d - то же, что "b" + дополнительная, внешняя. Дальний порядок и пространственная когерентность возникают в случаях: "b", "c" и "d". 1 - направление действия нагрузки; 2 - положение трещины и ее направление; 3 - направление пространственной когерентности.


Взаимодействие волн, возникающих при выходе ударной волны на поверхность.

Рассмотрим волновую картину процессов формирования в недрах Земли УВ, движение ее к дневной поверхности, возникновение ВР и взаимодействие волн. Воспользуемся приемом из практики изучения УВ. Изобразим рассматриваемые процессы на xt-диаграмме, где x - координата, а t - время. Обозначим: xо, tо - координату (глубину гипоцентра) и время формирования ударной волны в геологическом теле (см. рис. 5). Предположим, что в земной коре, на глубине между x = 0 и x = xо, и, далее, до x1, находится геологическое "тело", обладающее определенными уникальными акустическими свойствами: во-первых, это должно быть прочное тело, способное к накоплению упругой энергии и активному образованию трещин; во-вторых, оно должно быть изотропным. В таком теле не должно быть неоднородностей, у которых заметно различаются между собой плотность и скорость звука ( ρс). На таких неоднородностях происходит рассеяние (расщепление) УВ и генерация объемных волн. Характер нагрузки (напряжений) тела (в-третьих) должен быть таким, чтобы в нем могла возникнуть пространственно когерентная структура микротрещин или, иначе, - возникнуть (на языке физики когерентных структур) дальний порядок. В-четвертых, в этом теле должны раскрываться трещины, участвующие в процессе, - вполне определенного размера, не более 100 - 200 микрон, которые потом, после события, могли бы "залечиться" так, чтобы тело могло быть готово к очередному событию. (Это условие не исключает присутствие трещин других размеров). При выполнении этих требований в т.н. сейсмоактивной среде такого геологического тела возможно возникновение акустического стимулированного сверхизлучения (режима лавинного образования трещин) и возбуждение сильной акустической, центрированной волны, а затем и ударной волны.


Рис. 5. Волновая картина формирования ударной волны, волны разрежения и их взаимодействия в геологическом теле и на его поверхности. 1 - геологическое тело; 2 - земная кора; 3 - поверхность; 4 - фронт УВ; 5 - фронт волны разрежения (ВР); 6 - кривая распределения скоростей Р-волн по глубине.


В точке (x1, t1) УВ выходит на границу между геологическим телом и рыхлым грунтом, а в точке (x2, t2) - на земную (свободную) поверхность. В каждой из этих точек происходит отражение УВ и возникновение волн разрежения (ВР), "двигающихся" в обратном (по отношению к УВ) направлении. Длительность импульса (удара) УВ конечна (обозначим ее Δ t). За это время УВ охватывает пространство на земной поверхности Δx = Vx × Δt - в точке (x3, t3), где так же может возникнуть ВР. Таким образом, на поверхности Земли, в области эпицентра землетрясения, в пространстве (x1 - x3), - возникает ансамбль волн: ударных и разрежения. Эти волны, взаимодействуя между собой и с окружающим пространством, вызывают комплекс сложных волновых и гидродинамических явлений, приводящих, в конечном счете, к разрывам поверхности, разуплотнениям, вспучиваниям и проседаниям грунта, излиянию на поверхность жидкой фазы (мокрого песка и глины) в виде грифонов и т.п. На рис. 5 черным кружком (точка С) отмечено место пересечения ВР, в котором могут возникнуть разрывы сплошной среды. Рассмотрим эти явления подробнее.


Взаимодействие ударной волны с преградой, приводящее к возникновению отколов и разрушений поверхности

Взаимодействие волн. Рассмотрим механизмы взаимодействия волн, возникающих при выходе ударной волны на "свободную" поверхность твердого тела. Наиболее полно этот вопрос рассмотрен в книге [7], которой будем придерживаться при изложении материала этого раздела.


Рис. 6. Возможные типы распада разрыва [7]. Давления Р1 и Р2 - на фронтах УВ и ВР, давление РF - в теле, после прохождения по нему этих волн. Стрелки - направления распространения волн.


Как известно, ударная волна, это волна сжатия, на очень узком фронте которой термодинамические характеристики претерпевают разрыв. Эта волна, при прохождении через твердое тело приводит к его уплотнению и, при расщеплении на неоднородностях тела, - к генерации и излучению объемных волн. Ударная волна отличается от обычных акустических (сейсмических) волн ещё и тем, что вслед за фронтом УВ происходит движение (течение) вещества твердого тела, причем массовая скорость этого течения значительно меньше скорости УВ, и, чем интенсивнее ударная волна, тем выше массовая скорость. Совсем другая ситуация возникает, например, при выходе ударной волны на поверхность твердого тела. При этом образуется волна разрежения, двигающаяся в сторону, обратную движению УВ. Образование волн, идущих в обе стороны от "свободной" поверхности принято называть распадом разрыва. Разрыв в начальных условиях распространения УВ действительно исчезает, так как образуются два других разрыва. В книге Забабахина [7] рассмотрены четыре возможных типа распада разрыва, из которых нас будет интересовать лишь четвертый. Все типы распада разрыва изображены на рис. 4 (сверху вниз): 1) Образование двух ударных волн. Такой тип распада встречается при выходе УВ из более "мягкого" вещества в более "жесткое". 2) Образование ударной волны и волны разрежения. Это происходит при выходе УВ из более жесткого вещества в более мягкое. Массовая скорость частиц более жесткого вещества удваивается и в нем начинает распространяться волна разрежения. 3) Образование двух волн разрежения. Они возникают в том случае, когда, строго говоря, разрыва в начальных условиях нет, а волна разрежения приходит на границу двух сред, до этого сжатых и покоящихся. На Земле такая ситуация может возникнуть при выходе ударной волны на поверхность и образовании волны разрежения, которая приходит на границу двух сред, например, грунта, лежащего на гранитном основании, либо рыхлого грунта на границе вечной мерзлоты и т.п. Волна разрежения является фактически волной плавного изменения состояния и не носит характера ударной волны, т.к. фронт такой волны имеет конечную ширину. 4) Образование двух волн разрежения и разрыва сплошного вещества. Такое движение возможно, как и в ситуации 3, но при слабом начальном сжатии. При интенсивной волне разрежения элементы вещества начинают двигаться друг от друга в противоположных направлениях, т.е. происходит откол, - разрушение ранее целого и компактного вещества.

Смыкание откола. При соблюдении некоторых условий образовавшийся откол (разрыв сплошности материала) может сомкнуться. Для этого в материале должна возникнуть волна сжатия, т.е. ударная волна, которая может остановить процесс откола. Условия возникновения УВ могут быть связаны как с внутренними причинами, т.е. особенностями распространения волны разгрузки, так и внешними, - заключающимися в особенностях волновой обстановки вокруг материала, претерпевающего откол. Возможно, подобное явление, т.е. раскрытие трещины и, следующее за ним, смыкание её, наблюдал очевидец (пастух) во время Чуйского землетрясения в районе с Чибит. По словам очевидца, скорости раскрытия трещины и последующего её смыкания были сравнительно невелики, а сам процесс хорошо различался во времени. Опытные сейсмологи утверждают, что подобные случаи встречались и раньше, например, на Памире. Для объяснения этих фактов необходимо в рамках модели сильных движений получить оценки времени и скорости образования различимых глазом сейсморазрывов, которые должны быть сравнимы со скоростями прихода к сейсморазрывам ударных волн.

Влияние прочности материала на скорость откола. Обсуждая проблему откола (разрыва), предполагалось, что для проявления этого явления достаточно того, чтобы плотность вещества после выхода ударной волны на поверхность оказалась ниже исходной. Действительно, это не совсем верно. Для начала откола необходимо, чтобы плотность оказалась меньше заданной, заранее определенной. Это связано с тем, что вещество, через которое распространяется УВ, обладает прочностью σр.

О прочности грунта. Известно, что теоретическая прочность на разрыв твердого тела σр ≈ 0.1 Е, где Е - модуль Юнга. Учитывая справочные данные, скорости сейсмических волн в грунте и его плотность, можно с уверенностью считать, что величина Е ≈ 1 кбар. Фактически прочность меньше чем σp примерно в 100 - 1000 раз. Для грубых, качественных оценок примем σp = 1 - 10 бар, или 106 - 107 дин/см2, что соответствует давлению вышележащих слоев грунта р =σр = ρgh на глубине h порядка нескольких метров (обычная для землетрясений глубина разрыва поверхности).

Скорость сейсморазрыва. Воспользуемся подходом оценок эффектов, связанных с отколом, развитых в [7]. Будем полагать, что возникновение первичных трещин сейсморазрывов в грунте идентично появлению откола в твердом теле при воздействии на него ударных волн. В этом случае можно написать формулу для величины скорости раздвижения трещины U:

U = 1/ρo × P/x Δt,

где P/x - градиент давления, Δt - время процесса. Очевидно, что P/x = σр/δx, где σр - растягивающее напряжение, при котором происходит откол, δx - толщина слоя откола (ширина трещины сейсморазрыва). Полагая Δt - время выхода ударной волны на поверхность до момента образования откола: Δt = δ x/C, C - скорость звука, получаем выражение для скорости U:

U = σрo С,      (13)

что позволяет оценить порядок величины скорости (в нашей интерпретации) раскрытия трещины-разлома U. Скорость оказывается порядка нескольких метров в секунду. По-видимому, это несколько заниженная оценка, так как, со слов очевидцев, раскрытие трещины происходит, как говорится, "на глазах у свидетеля". По-видимому, человеческий глаз способен различить скорость порядка 100 м/с. Это потребует корректировки (увеличения в 10 раз) величины σр, что вполне правдоподобно. Важно, что скорость раскрытия трещины прямо связана с прочностью материала, что так же вполне правдоподобно.


Модели сильных движений Чуйского землетрясения

Рассмотренные нами явления, объединенные в группу, которая носит название сильные движения (макросейсмичность) при землетрясении, можно разделить по природе механизма. В первую очередь, это ударно-волновой механизм. О физике этого механизма в основном и шла речь выше. Однако некоторые сопутствующие эффекты имеют явно гидродинамический характер. Сюда с очевидностью можно отнести образование грифонов, образование валов вспучивания как своеобразных волн, типа волн на воде, правосторонние сдвиги по простиранию образовавшихся трещин, когда правый борт трещины сдвигается к югу, а левый, - к северу и т.п. Ударно-волновой механизм включает такие вопросы как, образование волн разрежения в результате выхода ударной волны на свободную поверхность и последующее взаимодействие этих волн между собой. Так как в большинстве случаев процессы взаимодействия волн происходят в сравнительно тонком слое грунта, лежащем на гранитном основании или на слое вечной мерзлоты, то, рассматривая эту задачу можно иметь в виду, что особенности распространения и взаимодействия волн в пластине подходят под описание т.н. волн Лэмба.

Гидродинамические задачи базируются на двух главных моментах, следующих из теории ударных волн. Во-первых, это учет того, что при распространении ударных волн вещество, через которое распространяется УВ, движется по направлению волны с массовой скоростью, значительно меньшей, чем скорость распространения ударной волны и волны разряжения. Это движение (течение) вещества направлено вдоль полиниям сжатия-растяжения. Во-вторых, следует учитывать, что сильная ударная волна при распространении через среду приводит её в "квазижидкое" состояние обладающее конечной величиной вязкости.

Кратко рассмотрим модели сильных движений в контексте физики ударно-волновой модели. Как отмечалось выше, среди множества сопутствующих землетрясению процессов авторы [1] выделили четыре вида нарушений:

1) сейсмотектонические разрывы, 2) обвалы и оползни, 3) вспучивание и проседание почвы и 4) излияние воды и грязи на поверхность связанное с образованием грифонов. Практически все виды нарушений нами рассмотрены выше при описании процессов, сопутствующих выходу ударной волны на поверхность среды. Тем не менее, дадим краткое содержание каждого из них.

Сейсмотектонические разрывы. Причины разрывов - взаимодействие двух и более волн разрежения, возникающих при выходе ударной волны на поверхность.

Обвалы и оползни. Эти явления возникают вследствие "подрезания склона сейсмотектоническими разрывами. Иногда такие явления называются сейсмогравитационными.

Вспучивание почвы. Возникает в поверхностном слое при отражении ударной волны от свободной поверхности, в то время как вещество ещё продолжает двигаться как бы вслед за УВ с массовой скоростью, которая в такой ситуации - удваивается.

Проседания почвы. Возможны две причины, или их комбинация. 1) Неравномерное вспучивание, там где меньше "вспучилось" - это и есть проседания; 2) Ударная волна, проходя по твердой среде, уплотняет её. При выходе на поверхность, возникает волна разрежения, которая разуплотняет верхний слой среды. В принципе, в силу различных геометрических, реологических и пр. обстоятельств, возможна ситуация, в которой волна разрежения, в некоторой области пространства, - не образовалась. В этом месте вместо вспучивания, будет наблюдаться проседание почвы.

Правостороннее смещение. На рис. 2 можно видеть, что направление трещин сейсмотектонических разрывов не совпадают с валами перемещения почвы. Углы расхождения могут достигать 45o. Несмотря на то, что оба явления возникают как результат выхода ударной волны на поверхность грунта, в первом случае картина обязана взаимодействию сильных волн: ударной и разрежения. Во втором случае, имеет место гидродинамический перенос, связанный с тем, что касательные напряжения, согласно закона Кулона-Мора, ориентированы под углом к действующей на тело нагрузке.

Образование грифонов. Напряжения, возникающие в грунте в процессе воздействия на него ударной волны определяются его реологией: чем суше грунт, тем в большей степени преобладают горизонтальные, касательные к поверхности грунта нагрузки над - вертикальными, нормальными к его поверхности. Соответственно, амплитуды перемещений частиц грунта, скорости и ускорения перемещений при воздействии на них ударной волны так же зависят от степени "сухости" грунта. Поэтому, во влагонасыщенном грунте возникают грифоны, а в сухом грунте - разрывы и образование трещин. Если вместо грунта, на поверхности находится монолитная порода, основное направление силы - вертикальное, в этом случае наблюдаются отколы. В сейсмологической практике известны случаи, когда в момент землетрясения откалывается часть скалы, "подпрыгивает", иногда разворачивается при этом, и остается на месте. Совсем другая физика образования грифона, которая, скорее, сродни физике, объясняющей возникновение "султанов" на воде при подводном взрыве [8]. И в том, и в другом случаях основным механизмом является схлопывание пузыря. Ударная волна выходит на поверхность и возникшее при этом давление выталкивает некоторое количество воды с массовой скоростью u. На поверхности воды образуется воронка в форме впадины.


Энергетика разрушений в эпицентре Алтайского землетрясения

Рассмотрим два различных подхода к оценке энергии разрушений в эпицентре землетрясения: один, - с точки зрения теории упругости, другой - гидродинамики. Обратим внимание на то, что в преобладающем большинстве работ по физике землетрясения принято считать потери на разрушения очень малыми, - менее одного процента от полной энергии землетрясения. Выясним, насколько справедливо это утверждение в случае Чуйского землетрясения.

Сейсмическая энергия. В механике землетрясения принято считать [6], что сейсмическая волна, достигая поверхности Земли, вызывает смещение земной поверхности x:

x = a cos (2 t/To),     (14)

тогда скорость движения грунта:

V = - (2ao/To) sin (2t/To).     (15)

Удельная энергия:

e = (ρ/4) (2t/To)2,

Кинетическая энергия, излучаемая очагом:

E = 33 h2 c to ρ (ao /To)2.    (16)

Здесь ао и То - амплитуда и период волны, h - глубина очага, с - скорость распространения волны, t = nTo - длина цуга волн.

Оценим величину (ao/To), принимая: с = 3,6 км/с, h = 18 км, ρ= 2.7 г/см3 и Е = 1023 эрг. Величина (ao/To) , иначе, скорость движения грунта, V ≈(5 - 10) см с-1. По-видимому, эту оценку можно принять как некоторое усреднение, т.к. в эпицентре 7.5 бального землетрясения скорость движения грунта значительно (в 10 - 20 и более раз) выше. Тем не менее, эта оценка полезна и в дальнейшем, обсуждая параметры сильных движений в эпицентре Чуйского землетрясения, будем ориентироваться на неё.

Характерно, что в формуле для Е (16) нет расстояния до места прихода излучения от очага. По-видимому, это имеет глубокий смысл. Например, если оценить величину σ =ρ сV, где с - скорость звука, то оказывается, что напряжение σ соответствует примерно 1 бар. Как было оценено выше, напряжение в волне в 1 бар было зарегистрировано на сейсмостанции в Новосибирске (4 балла). Однако, главный толчок Чуйского землетрясения в Барнауле вызвал сотрясение всего в 2 балла, хотя Барнаул в два раза ближе к эпицентру, чем Новосибирск.

Энергетика сильных движений. Воспользуемся формулой: p =ρDu, где примем: D = 3.5 км/с, u = 1м/с, р = 100 бар. Как отмечалось выше, эти параметры измерены при наблюдении Нортриджского землетрясения в Калифорнии и примерно соответствуют его энергетике (М = 6.7). Энергия Чуйского землетрясения (М = 7.5) примерно на порядок выше, следовательно, величины D и u, соответственно, выше.

Предположим невероятный случай, что вся энергия землетрясения Е пошла, например, на разуплотнение грунта на поверхности S = 10 × 10 км2 под действием волны разрежения. В этом случае глубина слоя разуплотнения грунта была бы равна E/pS =10 метрам. Эта оценка показывает масштаб явления. Реально наблюдаемые нарушения грунта в эпицентре Чуйского землетрясения были рассмотрены нами выше. Обратимся к рисунку 2 (рис. 11 [1]), на котором изображена структура сейсморазрывов и зон вспучиваний на водоразделе рек Талдура (на некоторых картах - Талтура) - Кускуннур в эпицентральной зоне разрушений Чуйского землетрясения. Грубая оценка суммарной длины сейсморазрывов, изображенных на этом рисунке, показывает, что это примерно 5 км. Если принять глубину разрывов 5 м, а ширину - 20 м, то объем, который подвергся разрушению, равен 5 105 м3. На эту "операцию" было затрачено примерно 5 1019 эрг. Оценим объем грунта, подвергшегося вспучиванию. Оказывается, что грубая оценка величины этого объема равна объему трещин. Возможно, это совпадение, но, возможно, -это реализация закона сохранения массы в ударной волне. В любом случае, суммарная энергия, потраченная на получение разрывов и вспучиваний, составляет примерно 1020 эрг на один из многих участков земной коры в эпицентре землетрясения. Будем, к примеру, считать, что таких участков 10, тогда получается, что на реализацию этого типа сильных движений и разрушений в эпицентре землетрясения "уходит" примерно 1 % от суммарной энергии. Пока не учитывалась часть энергии, которая "пошла" на работу грифонов, подвижку среды как квазижидкости и образование оползней.

Оценим, для сравнения, порядок кинетической энергии, израсходованной на образование грифона. Воспользуемся формулой Пуазейля для оценки объемного расхода жидкости V(см3/с), изливающейся через круглое отверстие радиусом R, длиной "трубы" l, вязкостью "жидкости" и под давлением Δp:

V = ( R4 /8 l) Δp      (17)

Подставляя R = 102 см, = 105 -106 Пз, l = 102 см, Δp = 108 дин/см2, получаем: V = 108 см3/с. Учитывая, что время "работы" грифона не более 10 с, масса (объем ΔV) выброшенного материала может достигать: m = 1000 т, что согласуется с приведенной выше оценкой. Линейная скорость переноса массы m порядка v = 104 см/с, что позволяет сделать грубую оценку кинетической энергии работающего грифона: Е = 1017 эрг. Величину Е можно оценить и по формуле: Е = Δp ΔV. Подставляя: Δp = 108 бар, ΔV =109 см3, получаем ту же величину Е. Даже если допустить, что во время Чуйского землетрясения "работало" 100 таких грифонов, их суммарная энергия (Е = 1019 эрг) оказывается значительно меньше (в сто раз), чем энергия, которая пошла на образование сейсморазрывов.


Пути и перспективы дальнейших исследований очага Чуйского землетрясения.

В то время, когда готовилась эта статья (декабрь 2003) в Калифорнии произошло землетрясение, погибли люди, разрушены дома. И это в том районе, где сосредоточена самая совершенная аппаратура, где ведется постоянный мониторинг всех возможных геофизических параметров, огромное количество геофизиков получают немалые деньги на прогноз землетрясения. А его нет, да и не может быть. Ещё в 1988 г. после Спитакского землетрясения, когда у меня впервые возник интерес к физике этого загадочного явления и пришла идея о возможности представить очаг землетрясения как источник генерации ударной волны, а механизм её возникновения, - как сверхизлучение звука когерентно раскрывающимися трещинами [9], в голову пришла мысль, которую я впервые опубликовал в 1990 в препринте, а затем повторял практически во всех своих статьях и книгах, касающихся физики этого явления: "Землетрясение, уже начавшись, само не знает, каким получится". Этой фразой можно определить отношение к прогнозу землетрясения. Получается так, что точный прогноз места, времени и силы землетрясения не возможен принципиально.

Какая же альтернатива? И есть ли она? Очень многие серьезные ученые и инженеры полагают, что в настоящее время у человечества существуют всего два пути: один, это хорошо строить и быть готовым к тому, что такое событие может произойти в любое время, и второй, - заниматься физикой явления. Причем оба пути должны развиваться одновременно, быть тесно связаны, а не противопоставляться друг другу.

Чуйское землетрясение представляет огромные возможности для реализации второго пути - изучения физики этого явления. Кое-что уже сделано, однако, это очень малая толика из того, что следовало бы сделать в будущем. Группой сотрудников Алтае-Саянской Опытно-Методической Сейсмологической Экспедиции (АСОМСЭ) Геофизической службы СО РАН определены глубины гипоцентров, координаты афтершоков, механизмы землетрясений и пр. Объединенной группой сотрудников ИФЗ (проф. Е.А.Рогожин) и Горно-Алтайского Государственного Университета (С.Г.Платонова) проведено предварительное изучение морфологии эпицентра этого землетрясения. Группой сотрудников Горного института УрО РАН (г. Пермь) (Р.А.Дягилев) проведены наблюдения микросейсмичности на сейсмостанции Акташ в период афтершоковой активости в начале ноября 2003 г. Там же, и в эти же сроки, сотрудником АСОМСЭ (С.Ю.Хомутов) проведены наблюдения за вариациями геомагнитного поля, инфразвуком и изменениями атмосферного электрического поля. (В этих наблюдениях принимал участие и автор). В настоящее время результаты наблюдений обрабатываются, можно надеяться, что полученные данные позволят хоть не намного приблизиться к пониманию физики землетрясения, а так же подтвердить, или опровергнуть ударно-волновую модель, суть которой частично изложена в этой статье.

Что следовало бы сделать в будущем? Отвечу на этот вопрос очень кратко, без особых обоснований:

  1. Продолжить наблюдения за микросейсмичностью, инфразвуком, атмосферным электрическим полем. Наблюдения за вариациями геомагнитного поля могут быть продолжены, но они, по всей видимости, бесперспективны, т.к. магнитометр на эффекте Оверхаузера, обладающий нулевой чувствительностью к колебаниям почвы, не зафиксировал ни одного сигнала, связанного с точками. Наблюдения велись в эпицентре землетрясения - на сейсмостанции Акташ.
  2. Наладить GPS-измерения за подвижностью грунта по обеим сторонам разлома Кускуннур.
  3. На ближайших к эпицентру землетрясения сейсмостанциях организовать регистрацию низкочастотного электромагнитного излучения. Аналогичные наблюдения вести в геофизической обсерватории Новосибирска. Это даст возможность отделить сигналы из зоны очага от излучения от внешних источников.
  4. Провести электромагнитное зондирование в районе разлома Кускуннур вдоль и поперек него. Это может быть как МТЗ (при выборе диапазона пульсаций, обеспечивающих проникновение зондирующего сигнала до глубин 20 км) или с помощью внешнего источника тока.
  5. Провести зондирование очага с помощью сейсмовибратора по трассам вдоль и поперек разлома Кускуннур. При выполнении п.п. 4 и 5 следует искать аномалии, связанные с очагом землетрясения.
  6. Изучить афтершоки, выяснить объем очага, определить, анализируя закон Омори, являются ли источниками афтершоков главный удар или в этом принимают участие и два других сильных землетрясения.
  7. Выяснить, действительно ли очаг Чуйского землетрясения, это долгоживущее образование, в котором немалую роль играет разлом Кускуннур. Для этого необходимо провести сравнение современных разрывов с палеонарушениями.
  8. Оценить глубину слоя вечной мерзлоты и выяснить его роль в особенностях сильных движений.

Благодарности. Благодарю сотрудников Геофизической службы СО РАН (АСОМСЭ) С.Ю.Хомутова - за предоставление возможности изучения последствий землетрясения и афтершоковой активности в близи его эпицентра и А.Г.Филину - за многочисленные консультации и обсуждения представленного в этой работе фактического материала. Благодарю сотрудника Горно-Алтайского Университета А.Ю.Гвоздарева, участника экспедиции, за информацию об особенностях сильных движений Чуйского землетрясения. Хочу особо отметить авторов статьи [1], опубликованной в Вестнике, в которой наиболее полно освещены явления сильных движений, имевших место при Чуйском землетрясении 27 сентября 2003 г.

Литература

  1. Геодаков А.Р., Овсюченко А.Н., Платонова С.Г., Рогожин Е.А. Материалы предварительного изучения сильного землетрясения 2003 г. на Горном Алтае. Электронный научно-информационный журнал "Вестник отделения наук о Земле РАН" N 1 (21) 2003.
  2. Ляхов Г.М., Покровский Г.И. Взрывные волны в грунтах. М.: ГНТИ. 1962 г. 103 с.
  3. Станюкович К.П. Неустановившиеся движения сплошной среды. М.: Наука. 1971. 854 с.
  4. Зельдович Я.Б., Райзер Ю.П. Физика ударных волн и высокотемпературных гидродинамических явлений. М.: Наука. 1966.
  5. Николаевский В.Н. Динамическая прочность и скорость разрушения. Сб. Механика. Удар, взрыв и разрушение. 26. М.: Мир. 1981. С. 166-203.
  6. Касахара К. Механика землетрясений. М.: Мир. 1985. 264 с.
  7. Забабахин Е.И. Некоторые вопросы газодинамики взрыва. Снежинск: РФЯЦ-ВНИИТФ. 1997.
  8. Лаврентьев М.А., Шабат Б.В. Проблемы гидродинамики и их математические модели. М.: Наука. 1973
  9. Кузнецов В.В. Физика земных катастрофических явлений. Новосибирск. Наука. 1992; Ударно-волновая модель очага землетрясения. Науки о Земле. М.: Вузовская книга. 2001. С.47-62; Ударно-волновая модель землетрясения. Вулканология и Сейсмология, принято в печать 2004. Физика горячей Земли. www.uiggm.nsc.ru/~kuz/site.htm