Структура и динамика геофизических полей и сейсмических процессов в блоковой модели земной коры

Трофименко, Сергей Владимирович

Бесплатный автореферат и диссертация по наукам о земле на тему
Структура и динамика геофизических полей и сейсмических процессов в блоковой модели земной коры
ВАК РФ 25.00.10, Геофизика, геофизические методы поисков полезных ископаемых

Автореферат диссертации по теме "Структура и динамика геофизических полей и сейсмических процессов в блоковой модели земной коры"

4856273

ТРОФИМЕНКО СЕРГЕЙ ВЛАДИМИРОВИЧ

СТРУКТУРА И ДИНАМИКА ГЕОФИЗИЧЕСКИХ ПОЛЕЙ И СЕЙСМИЧЕСКИХ ПРОЦЕССОВ В БЛОКОВОЙ МОДЕЛИ ЗЕМНОЙ КОРЫ

25.00.10 - геофизика и геофизические методы поисков полезных ископаемых

АВТОРЕФЕРАТ диссертации на соискание ученой степени доктора геолого-минералогических наук

2 4 20//

Томск-2011

4856273

Работа выполнена в Техническом институте (филиале) федерального государственного автономного образовательного учреждения высшего профессионального образования "Северо-Восточный федеральный университет имени М.К. Аммосова"

Научный консультант: доктор технических наук

Гриб Николай Николаевич

Официальные оппоненты: доктор физико-математических наук

Викулин Александр Васильевич

доктор геолого-минералогических наук Дмитриев Александр Георгиевич

доктор геолого-минералогических наук Исаев Валерий Иванович

Ведущая организация: Институт физики Земли им. О.Ю. Шмидта РАН (г. Москва)

Защита состоится «27» апреля 2011 г. в 15— на заседании совета по защите докторских и кандидатских диссертаций Д 212.269.12 при Государственном образовательном учреждении высшего профессионального образования «Национальный исследовательский Томский политехнический университет» по адресу: 634050, г. Томск, пр. Ленина, 2, строение 5, (20 корп. ТПУ), ауд.504.

С диссертацией можно ознакомиться в Научно-технической библиотеке ГОУ ВПО «НИ ТПУ» (634050, г. Томск, ул. Белинского, д.55).

Автореферат разослан «26 » января 2011г.

Ученый секретарь совета, доктор геолого-минералогических наук

А.А. Поцелуев

ОБЩАЯ ХАРАКТЕРИСТИКА РАБОТЫ

Актуальность темы.

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

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

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

Исторически, начиная с классических работ Г.А. Гамбурцева по прогнозу землетрясений, изучение сейсмотектонических процессов проводилось в контексте комплексных геолого-геофизических исследований в сейсмоактивных регионах с развитием методов натурных полевых регистрирующих систем геофизических процессов и методов моделирования физики очага землетрясения и природы сейсмичности в целом. Природа сейсмичности рассматривалась во всех основополагающих работах по физике Земли. В различных подходах данные вопросы изложены в ставших теперь классическими монографиях [Джеффрис, 1960; Гуттенберг, 1963; В.А. Магницкий, 1965; Стейси, 1972; Ботг, 1974; Буллен, 1978; Цубои, 1982; Эйби, 1982; Гир, Шах, 1988]. Новые подходы к анализу сейсмичности изложены в монографиях [Садовский, Писарснко, 1991; Садовский, 2004; Соболев, 1993; Быков, 2000; Викулин, 2003; Долгих, 2004; Ребецкий, 2007; Гуфельд, 2007; Хаин, Халилов, 2008].

Существующие представления об источниках аномалий геофизических полей предвестников землетрясений [Гохберг и др., 1988; Электромагнитные предвестники ..., 1982; Виноградов, 1989] базируются на физических [Физика очага ...,1975; Физические основания ...,1970] моделях очага землетрясения. Традиционный путь решения прогнозных задач - анализ корреляционных связей между аномальными проявлениями в физических полях и пространственным распределением, механизмами и динамикой очагов землетрясений с привлечением геоморфологических, геологических, тектонических и космических критериев сейсмичности. За практически 50-ти летний период разработки данного на-

правления исследований наиболее детально изучено влияние катастрофических землетрясений на динамику геофизической среды. Не исследованными остаются вопросы формирования аномальных полей вследствие миграционных волн и временных циклов сейсмической активности.

За рамками давдых исследований остаются также вопросы по комплексной интерпретации и моделированию результатов мониторинга геологической среды, структурированию сейсмического процесса, по взаимодействию сейсмических зон в пределах переходных областей, по анализу многообразных сейсмологических связей и геофизических полей. Открытыми остаются вопросы по изучению состояния геологической (деформируемой) среды в зонах активных разломов геофизическими методами в современный (го-лоценовый) период.

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

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

Объектом данного исследования являются физические поля магнитной, гравитационной и электромагнитной природы, иродуцируегиые переходной зоной взаимодействия Амурской, Евроазиатской и Охотоморской литосферных плит.

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

построение пространственных моделей распределений аномалий геофизических полей методом структурного статистического анализа;

изучение пространственно-временных характеристик сейсмических процессов;

разработку комплексной геодинамической модели формирования геофизических полей и сейсмических процессов контроля и прогноза сейсмической активности.

были поставлены и решены следующие основные задачи:

1. проведение статистических исследований по распределению азимутов линейных элементов аномалий физических полей гравитационной и магнитной природы с выделением преобладающих азимутов, построение моделей (систем) по распределению плотно-

сти индикаторов линейных элементов аномалий и поиск закономерностей их пространственной структуры.

2. установление подобия выделенных систем структурным тектоническим элементам исследуемой территории.

3. построение статистических моделей пространственно - временных распреде-

лений слабых землетрясений с энергией Е < 10Дж и изучение их закономерностей.

4. анализ динамических характеристик сейсмического режима сильных землетрясс-

ний Е > 10 Дж и рассеянной сейсмичности за инструментальный период. Проверка выявленных закономерностей на очагах сильных землетрясений Байкальской рифговой зоны, Сейсмического пояса Черского, Тихоокеанской островной дуги (п-ов Камчатка, о. Сахалин), обобщение материалов модельных построений и вьивление общих закономерностей статистических моделей сейсмического процесса.

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

6. проведение мониторинговых геофизических исследований в зонах активных разломов для оценки их современной тектонической активности и изучение прогностических эффектов с разработкой технологии моделирования геофизических полей и сейсмических процессов переходных зон и построения геодинамической модели сейсмического процесса.

Методы исследований. Системный анализ и научное обобщение результатов теоретических исследований взаимосвязи геофизических полей и сейсмических процессов; обработка статистических данных, построение прогнозных статистических зависимостей прогнозирования состояния геофизической среды; моделирование сейсмических процессов в зонах взаимодействия различных сейсмогенерирующих зон; выявление и систематизация вариаций геофизических полей в зонах активных разломов в периоды сильных землетрясений в условиях изменяющихся геодинамических условий.

Научные положения, защищаемые в диссертации:

1. Пространственное распределение плотности линейных элементов аномалий гравитационного и магнитного полей (индикаторов) обнаруживает высокую степень упорядоченности в виде ромбических структур (систем доменов) и является следствием блокового строение геофизической среды.

2. Проявления статистических закономерностей сейсмичности в характерных чувствительных зонах обусловлены влиянием общих возмущающих факторов вследствие взаимодействия физических полей в системе «Солнце-Земля-Луна» и идентично для быстрых геодинамических процессов в виде афтершоковых последовательностей и медленных стационарных сейсмических процессов.

3. Зона релаксации энергии упругих деформаций вследствие сильных землетрясений охватывает область активных разломов (по сейсмологическому признаку), в виде устойчивой геодинамической системы не зависимо от механизма очага землетрясения, энергии землетрясения, начиная с некоторой пороговой величины и положения очага внутри области. Длительность процессов релаксации пропорциональна энергии землетрясения.

4. Экспериментально установлены закономерности проявления аномалий физических полей, сформированных динамикой сейсмотектонических процессов переходной зоны взаимодействия Амурской, Евроазиатской и Охотоморской литосферных плит. Показано, что источниками аномалий являются геофизические процессы тектонической природы, а характерные особенности вариаций геофизических полей обусловлены реологией земной коры в пункте регистрации и отражают блоковую структуру литосферы.

Научная новизна работы:

В рамках учения М.А. Садовского о блочном строении земной коры разработана комплексная геодинамическая модель сейсмического процесса в стационарной и не ста-

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

• методами статистического анализа распределершй физических полей гравитационной и магнитной природы определен линейный размер неоднородности геофизической среды (домена). Впервые показано, что доменные структуры с линейными размерами от 120 до 150км образуют взаимосогласованные пространственные решетки, тем самым проведено обобщение теоретической модели формирования транспрессионных структур Вуд-кока - Фишера (Woodcock N.II., Fisher М, 1986.);

• развиты теоретические модели сейсмических процессов в переходных зонах ли-тосферных плит. Впервые показано что пространственное распределение очагов землетрясений формируется в векторном поле внешних сил, а временное распределение землетрясений в течение суток идентично для быстрых геодинамических процессов в виде аф-тершоковых последовательностей и медленных (стационарных) сейсмических процессов, тем самым доказан принцип согласованности самоподобия пространственной структуры земной коры и автомоделъности сейсмического процесса;

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

Фактический материал. В работе используются данные многолетних наблюдений вариаций геофизических полей на геодинамических полигонах. Для построения статистических моделей геофизических полей использовались карты м-ба 1:500 ООО. Основная часть исследуемых данных получена Проблемной научно-исследовательской лабораторией комплексных геофизических методов поиска предвестников землетрясений Якутского госуниверситета (ПНИЛЗ ЯГУ) в процессе работы по госбюджетным тематикам, а также при проведении детальных геофизических исследований зон активных разломов.

Личный вклад автора. В основу диссертации положен материал, полученный автором при моделировании геофизических полей и сейсмических процессов. Всеми перечисленными вопросами автор занимался в качестве ответственного исполнителя и исполнителя госбюджетных и хоздоговорных НИР Проблемной научно-исследовательской лаборатории комплексных геофизических методов поиска предвестников землетрясений Якутского госуниверситета (ПНИЛЗ ЯГУ) и Технического института (филиала) СевероВосточного федерального университета им. М.К. Аммосова в г. Нерюнгри (ТИ(ф) СВФУ) и непосредственно участвуя в проведении полевых геофизических работ в течение 1980 -2009гг.

При содействии заведующего ПНИЛЗ ЯГУ A.C. Стативы автором организованы и проведены деформационные и др. виды геофизических исследований в Олекминском эпи-центральном поле катастрофических землетрясений. В период Южно-Якутского землетрясения 1989г. с мапштудой М=6.6, при исследовании геофизических полей-предвестников землетрясений автором установлены закономерности формирования вариаций геофизических полей вследствие сейсмотектонических процессов и выделены характерные периоды вариаций сопоставимые с периодами сейсмических процессов.

Совместно с H.H. Гриб организованы мониторинговые геофизические исследования на базе ТИ(ф) СВФУ. Изучение вариаций электромагнитного излучения - импульсного электромагнитного поля земли (ЭМИ-ИЭМПЗ) в период повышенной сейсмической активности 2005-2007гг. в сопоставлении с результатами 1989г. позволило автору сформулировать основные положения методики и технологии мониторинговых наблюдений в

блоковой модели земной коры и определить ограничения на возможности геофизических методов прогнозирования землетрясений.

При содействии В.М. Никитина автором организованы и проведены детальные геофизические исследования зон активных разломов. На основе всей совокупности данных, включая пространственное распределение эпицентров рассеянной и очаговой сейсмичности в сопоставлении с азимутами активных разломов, детальные геофизические и мониторинговые исследования, разработана физико-динамическая модель сейсмической активности в блочной модели строения земной коры

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

Научное и практическое значение.

Полученные результаты моделирования дополняют систему знаний в области современных фундамею-альных исследований геодинамических процессов и прогноза землетрясений геофизическими методами мониторинга геологической среды, вследствие чего, результаты проведенных исследований могут быть использованы для уточнения и прогноза дифференцированных зон ВОЗ.

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

Предложенная модифицированная модель Беньофа может быть использована в сейсмологических отчетах для среднестатистического прогноза периодов повышенной сейсмической активности.

Геолого-геофизические материалы автора по центральной части Алданского щита вошли в отчет института Физики Земли РАН и использованы при проектировании трубопроводной системы ВСЮ в части уточнения сейсмической опасности.

Модельные построения по изучению пространственной неоднородности геофизической среды методами статистического анализа распределении линейных элементов аномалий гравитационного и магнитного полей могут иметь практическую значимость для уточнения и прогноза дифференцированных зон сейсмической опасности, т.к. выделенные азимуты и границы ромбических структур (систем доменов) коррелируют с пространственным распределением очагов землетрясений.

Апробация работы и публикации.

Результаты диссертационной работы докладывались на следующих конференциях: пятые геофизические чтения имени В. В. Федынского 27 февраля - 01 марта 2003 года-Москва, центр ГЕОН,(2СЮЗ); международной конференции «Проблемы сейсмологии 1П-го тысячелетия», Новосибирск (2003); 2-м Международном симпозиуме «Активный геофизический мониторинг литосферы Земли», Новосибирск (2005); международной сейсмологическая школа «Современные методы обработки и интерпретации сейсмологических данных», Обнинск (2006, 2008, 2009, 2010); IV-VI Всероссийских симпозиумах «Сейсмоа-кустика переходных зон», Владивосток, ТОЙ ДВО РАН (2005, 2007,2009); Всероссийской научно-практической конференции «Сейсмичность Южно-Якутского региона и прилегающих территорий», Нерюнгри (2005); IV-V Международной конференция «Солнечно -земные связи и предвестники землетрясений», с. Паратунка, Камчатской обл. (2007, 2010); Всероссийском совещании с международным участием «Проблемы современной сейсмологии и геодинамики Центральной и Восточной Азии», Иркутск ИЗК СО РАН (2007); международной научно-практическая конференция «Южная Якутия - новый этап индустри-

алыюго развития», Нерюнгри (2007); IV Международном симпозиуме «Геодинамика внутриконтинентальных орогенов и геоэкологические проблемы», Бишкек (2008); Всероссийской конференции «Тсктонофизика и актуальные вопросы наук о Земле», Москва -ИФЗ РАН - 20-25 октября 2008; Всероссийской конференции «Тектоника и глубинное строение востока Азии» Хабаровск: ИТиГ им. Ю.А. Косыгина ДВО РАН -20 -23 января 2009; Всероссийском совещании с участием приглашенных исследователей из других стран «Разломообразование и сейсмичность в литосфере: тектонофизические концепции и следствия», Иркутск ИЗК СО РАН (2009); на второй научно-технической конференции «Проблемы комплексного геофизического мониторинга Дальнего Востока России» 11-17 октября 2009г. г.Петропавловск-Камчатский; на научно-практической конференции «Сейсмические исследования земной коры» (Пузыревские чтения) Новосибирск: Институт нефтегазовой геологин и геофизики им. A.A. Трофимука СО РАН 22-25 ноября 2009; на научном симпозиуме «Кайнозойский континентальный рифтогенез» Иркутск ЙЗК СО РАН (2010); на научном симпозиуме «Проблемы сейсмичности и современной геодинамики Дальнего Востока и Восточной Сибири» Хабаровск: ИТиГ им. Ю.А. Косыгина ДВО РАН -1 -4 июня 2010.

Структура и объем работы. Диссертация состоит из введения, пяти разделов основного текста, заключения, изложенных на 240 страницах машинописного текста, содержит 91 рисунок, 4 таблицы, список литературы из 247 наименований и 1 приложения.

Благодарности. Автор выражает благодарность д.т.н. Грибу H.H., д.г.-м.н. Никитину В. М., д.г.-м.н. Имаеву В. С., за многолетнее содействие научному исследованию, а A.C. Стативе за сотрудничество в процессе работы, д.г.-м.н. Рогожину Е.П., поддержавшему работу на завершающем этапе исследований; к.г.-м.н. Овсюченко А.Н. за плодотворное сотрудничество при интерпретации материалов геолого - геофизических исследований.

КРАТКОЕ СОДЕРЖАНИЕ РАБОТЫ

В ведеиии изложены актуальность проблемы, цели и задачи исследований. Сформулированы основные защищаемые положения и определена общая структура диссертации.

Глава 1. Модели геофизической среды, вариаций физических нолей сейсмотектонической природы и сейсмических процессов

Первая глава носит обзорно-аналитический характер. Здесь приводится анализ состояния вопросов и подходов к изучению сейсмотектонических процессов и физических полей и закономерностей их проявления.

Модели геофизической среды и сейсмических процессов

Современные взгляды на сейсмический процесс, как на один из видов геотектонических процессов, развиваются па принципах автомоделыюсти, использующих признаки перехода от случайного, хаотического распределения сейсмических событий к их пространственной и временной самоорганизации, упорядоченности [Пригожин, Стенгерс, 1989; Кузнецов, 2000,2008]. Признаки самоорганизации должны проявляться в различных формах детерминированных проявлений аномалий сейсмического режима, с эффектами цикличности и периодичности сейсмической активности, миграции очагов землетрясений и т.д. Разрушения самоорганизованных систем, т.е. обратный переход к хаосу и разупоря-доченности, происходит в результате влияния внешних дестабилизирующих факторов на геофизическую среду. Методика анализа сейсмического процесса в блоковой среде [Садовский, Писаренко, 1991], развивающегося во времени и состоящего из последовательности различных этапов, связанных друг с другом и характеризуемых определенной временной последовательностью, предопределяется тем, что землетрясения являются одним из проявлений самоорганизующегося процесса в системе «литосфера - физические поля -сейсмический процесс». Геофизическая среда приспосабливается в процессе энергетического обмена к внешней силовой нагрузке, самостоятельно изменяя свою структуру, в результате чего, составные части среды могут консолидироваться в единое целое и разру-

шаться [Садовский, Писаренко, 1991]. Процессы приспособления среды происходят вблизи от некоторого положения равновесия, в соответствии со средним состоянием ее энергоемкости, и определяются глобальным взаимодействием тектонических плит. Следствием данного подхода должна быть долговременная пространственная упорядоченность сейсмичности, характеризующаяся местоположением очагов землетрясений, отражающая структуры современного поля деформаций. Основные положения учения о блочном строении земной коры сформулированы М.А. Садовским [Садовский, 2004] и заключаются в следующем:

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

2. Среда не только поглощает, но и, находясь в метастабильном состоянии, постоянно излучает энергию - сейсмическую, электрическую, магнитную, радиационную и другие. В системе «литосфера - океан - атмосфера» происходит постоянный энергетический обмен вследствие энергетической насыщенности среды.

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

4. Изменчивость физических свойств среды во времени происходит главным образом в результате изменения деформационных процессов вследствие перераспределения силовых нагрузок эндогенной, экзогенной и космической природы, и связи с этим, реакция физических параметров геологической среды отображается в геофизическом пространстве - времени во всех геосферах Земли.

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

Катастрофы порождают практически необратимые изменения среды, которые описываются глубиной воздействия, его объемом и другими пространственно-временными характеристиками. При этом в неравновесной сложной системе, какой является земная кора, возникают уникальные события и флуктуации, способствующие возникновению этих событий, а также происходит расширение масштабов системы [Николаев, 2003].

Различные модели блоковой структуры и неоднородпостей литосферы, проясняющие природу движений блоков, источников сейсмического процесса, сейсмического режима литосферы и сильных землетрясений изложены в работах [Гуфельд, 2009; Левин и др., 2007,2009, Ребецкий, 2008].

В качестве другого нового метода к анализу сейсмичности использован вероятностный подход, на основании которого впервые в континуальной форме созданы электронные слои карт периодов повторяемости сейсмических сотрясений с различной интенсивностью для всей территории Северной Евразии, охватывающей территорию Россию, всех стран СНГ и сопредельные сейсмоактивные регионы. Основой для данного вида исследований явилась электронная база сейсмологических данных, созданная в 1991-1997 гг. авторами [Уломов, 1999; Уломов, Шумилина, 1999] при выполнении работ по общему сейсмическому районированию России - ОСР-97. Вероятностный анализ сейсмической опасности сочетает в себе альтернативные модели сейсмических очагов, периоды повторяемости землетрясений, зависимости затухания сейсмического эффекта с расстоянием, а также целый ряд неопределенностей, обусловленных как нечетким знанием тех или иных параметров, так и случайным характером возникновения самих сейсмических событий.

Модели геофизических полей сейсмотектонической природы и проблема прогноза землетрясений

Со времени первых работ Ф. Рейда [Яе1(1, 1910, 1911] было разработано множество моделей по физике и динамике очага землетрясения и сейсмического процесса в целом. Совершенствование подходов к моделированию сейсмичности, основанных на теории взаимодействия всех геосфер Земли [Хаин, Короновский, 2007], модели блокового строения геофизической среды [Садовский, 2004], волновой природы сейсмического процесса [Быков, 2000, 2005; Викулин, 2003, 2008], порождает новую проблему, связанную с возможностью проявления вариаций физических полей вследствие геодинамических процессов (геофизических полей). Пространственное распределение геофизических полей и их динамические (временные) аномалии в виде геофизических предвестников землетрясений в процессе формирования и разрушения консолидационной зоны на разных временных масштабах не могут быть обусловлены только локальными причинами. Традиционные представления о возможности локального выделения аномалии предвестника от формирующегося очага землетрясения должны трансформироваться в неопределенность пространственного положения, соразмерного области консолидации, длины волны медленных деформационных волн и взаимодействием геофизических сред в переходных зонах. Отсюда, по видимому, следуют многочисленные ошибочные прогнозы в определении места, времени и силы землетрясения.

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

Формирование геофизических полей в результате деформационных процессов возможно при сохранении структуры поля деформаций в течение длительных интервалов времени. Не привлекая геологические и тектонические причины деформаций можно предположить возможность структурирования аномалий геофизических полей в результате самоорганизации согласно разработанным моделям блоковой геологической [Пейве, 1961] и геофизической [Садовский, 2004] сред. В данной концепции глобальная общепланетарная самоорганизация нашла свое отражение в теории ротационного структурообра-зования разломной тектоники [Тяпкин, Гонтаренко, 1990] и, как следствие, в возможности формирования устойчивых разнонаправленных систем распределений геофизических полей.

Физическая интерпретация изменения геофизических параметров среды предложена Т. Нагата в конце 60-х годов 20 столетия [Ыада1а Т., 1968, 1972]. Результаты данного направления исследований обобщены в работе [Сковородкин, 1985]. Полный анализ исследований наблюдений электромагнитного излучения (ЭМИ) приведен в работе [Гохберг, 1998]. Изучение неприливных изменений силы тяжести и их связи с современными сейсмотектоническими процессами представлено многими коллективами авторов и рассматривалось на ежегодных совещаниях под руководством Ю.Д. Буланже [Буланже, 1983]. Деформационные и наклономерные наблюдения на геодинамических полигонах с особой тщательностью были проанализированы в период Спитакского землетрясения [Деформационные ..., 1989].

Сейсмотектоническая позиция изучаемой территории и формулировка проблемы исследований

В сейсмическом отношении территория исследований охватывает сейсмические пояса Северо-Восточной Азии, сейсмичность которой обусловлена интенсивным геодинамическим взаимодействием нескольких крупных литосферных плит - Евроазиатской

(ЕАЗ), Амурской (AM). Тихоокеанской (ТИХ), Северо-Американской (САМ) и Охотомор-ской (ОХО) [Имаев и др., 2000] (рис.1).

Очаги землетрясений различной величины и глубины гипоцентров Ш 0 - 70 км

71 - 300 км *» 301 - 700 км

Границы питосферных ппит т ш Оси конвергентных сейсмоактивных структур - ■*•.• ' * ? :Л

зоны субдукцни и их реликты на континенте . - ' '•У?' 7

Направление миграции сейсын' в 22-х летних циклах

120' 50' 140'

Рис.Х. Границы литосферных плиг и сейсмические пояса Северо-востока Азии (с использованием материалов B.C.. Имаева).

Литосферных плиты (EA3-EU - Евроазиатская, AM - Амурская, ОХО-ОК - Охотоморская, CAM-NA - Севере - Американская, ТИХ-РА - Тихоокеанская); показано распределение очагов землетрясений в пределах сейсмоактивных зон: I - Верхояно-Колымская, П- Байкало-Становая (Па - восточное замыкание БРЗ, 116 Олекмо - Становая зона), Ш - Сахалино-Японская, IV - Арктическая рифтовая, V - Курило-Камчагасая

Наиболее активными являются конвергентные структуры литосферы. Они достаточно упорядочены по своим размерам и представлены дугообразными границами между ли-тосферными плитами, расположенными по периферии Тихого океана в виде зон субдукции [Уломов, 2007].

В тектоническом отношении область моделирования объединяет восточное замыкание Байкальской рифтовой зоны и Олекмо-Становое орогеническое поднятие [Имаев, 2000]. При моделировании сейсмических процессов привлечены материалы по Верхояно-Колымскому региону, Арктической рифтовой зоне, часть Курило-Камчагской зоны в пределах (56° - 58°СШ, 160° - 165° ВД), часть Сахалино-Японской островной дуги в пределах о. Сахалин.

Байкальская рифтовая зона - уникальный сейсмогеодинамический регион мира. Впадина озера представлена тремя сейсмоактивными котловинами - южной, средней и северной. Аналогичная зональность свойственна и проявлению сейсмичности восточнее озера, вплоть до р. Олекмы. Восточнее БРЗ Олекмо-Становая сейсмоактивная зона трассирует границу между Евроазиатской и Китайской литосферными плитами (некоторые исследо-

ватели выделяют еще промежуточную, меньшую по площади, Амурскую плиту). Среди относительно недавних крупных землетрясений - Муйское (1957 г., М-7.7,10—10 баллов).

Рассматриваемая часть Верхояно-Колымского региона соответствует сейсмическому поясу Черского, протягивающемуся в юго-восточном направлении от устья р. Лены к побережью Охотского моря, Северной Камчатке и Командорским островам. Самые сильные из известных в пределах пояса землетрясений - два Булунские (1927 г., М=6.8 и 10=9 баллов каждое) в низовьях р. Лены и Артыкское (1971 г., М=7.1, 1сг9 баллов) - у границы Якутии с Магаданской областью.

Арктическая рифтовая зона является северо-западным продолжением сейсмоактивной структуры Верхояно-Колымского региона, уходящей узкой полосой в Северный Ледовитый океан и соединяющейся на западе с аналогичной рифтовой зоной Срединно-Атлантического хребта. На шельфе моря Лаптевых в 1909 г. и 1964 г. произошли два землетрясения с магнитудой М=б.8.

Курило-Камчатская зона является классическим примером субдукции Тихоокеанской литосферной плиты под материк. Она протягивается вдоль восточного побережья Камчатки, Курильских островов и о-ва Хоккайдо. Здесь возникают самые крупные в Северной Евразии землетрясения с М более 8.0 и сейсмическим эффектом 1о=Ю баллов и выше. Структура зоны четко прослеживается по расположению очагов в плане и на глубине. Протяженность ее вдоль дуги около 2500 км, по глубине - свыше 650 км, толщина -около 70 км, угол наклона к горизонту - до 50°.

Сахалин представляет собой северное продолжение Сахалино-Японской островной дуги и трассирует границу Охотоморской, Амурской и Евразиатской литосферных плит. Катастрофическое Нефгегорское землетрясение (1995 г., М=7.5) было самым разрушительным из известных за все время на территории России. Погибло более 2000 чел. В результате полностью ликвидирован рабочий поселок Нефтегорск. Крупным сейсмическим событием было Углегорское землетрясение в южной части острова (2000 г., М=7.0), что подтвердило повышенную сейсмическую опасность Сахалина.

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

Интерпретация результатов экспериментальных данных в системе прогнозирования землетрясений и быстрых геодинамических процессов в литосфере осуществляется в рамках классической теории линейных однородных сред и существенно отстает от теории взаимодействия геосфер [Садовский, 2004; Викулин, 2003] и разрабатываемых волновых [Быков, 2005] и вихревых [Викулин, 2008] моделей геосфер. Новые подходы и концепции формирования физических полей в геофизической среде обусловлены нелинейными свойствами самоорганизованных, энергонасыщенных геосфер и их постоянном взаимодействии.

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

Глава 2. Структурная статистическая модель распределения физических полей гравимагцитной природы

Предпосылки для построения статистических пространственных моделей физических натек

Цубои Т. [Цубои, 1982] при изучении изостазии Северо - Американского континента методом разложений в двойной ряд Фурье аномалий силы тяжести в редукции Буге установлены пространственные параметры периодичности структур рельефа L = 140 км вдоль широты.

Авторы [С.И. Шерман и др., 2005], изучая закономерности строения блоковой среды по дискретному распределению средних размеров блоков (L) для регионов с различными режимами геодинамического развитая (Камчатка, Тянь-Шань, Байкальская рифтовая зона и сопредельные территории), установили полимодальный характер распределения параметра L.

Общей выделенной закономерностью данных построений явилось наличие локального максимума среднестатистических размеров блоков протяженностью L = (110 - 130) км. Для БРЗ характерно наличие максимума Ц =\/2L = 60км.

Исследования афтершоковых областей сильных землетрясений, происшедших в пределах северо-западной окраины Тихого океана в XX столетии [Викулин, 1987, 1990] показывают, что существует такое значение магнитуды землетрясений М^, для которых взаимное расположение афтершоковых областей имеет следующую особенность: афтер-шоковые области землетрясений с М > М^ имеют тенденцию не перекрывать друг друга, в то время как афтершоковые области землетрясений с М < Msp практически полностью расположены внутри афтершоковых областей более сильных событий. Линейный размер не пересекающих зон составляет L = 150 - 200км.

Изучению пространственного размещения аномалий магнитного и гравитационного полей литосферы посвящены многочисленные публикации, обобщение результатов которых и их тектоническая интерпретация приведеда в монографиях: (Малышев, 1977; Парфенов и др., 1985; Фотиади и др., 1988, Абрамов, 1993; Стогний, В.В. Стогний, 2005).

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

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

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

• ступени в уровнях полей (используются осевые линии линейных гравитационных и магнитных ступеней);

• локальные аномалии полей обоих знаков (используется большая ось эллиптических и квазиэллиптических аномалий);

• нарушения регулярного поведения изоаномал гравитационного и магнитного полей (используются линейные границы областей с разным рисунком этих полей).

Для построения статистических моделей распределения линейных элементов геофизических полей использовались электронные карты аномалий гравитационного и магнитного полей в масштабе 1 ;500 ООО в координатах 119 - 135° ВД и 55 - 60° СШ. При сканировании карг м - ба 1:500 000 на стандартный лист А4 помещается участок (фрагмент) размером 1.5x1.5°. Координаты всех видимых линейных элементов (признаков индикаторов разломов) по всем фрагментам рассматриваемого поля заносятся в электронные таблицы Excel. По заданным формулам вычисляется длина каждого линейного элемента (отрезка) и его азимут, общая длина отрезков в заданном азимуте и значения частот индикаторов, нормированные на максимальную величину отрезка по каждому полю с учетом округления до целых. Посредством модулей «описательная статистика, гистограмма, скользящее среднее и др.» встроенного пакета «Анализ данных» производятся стандартные статистические процедуры для оценки значимости выделенных аномальных уровней, оценка погрешностей и построения графиков. Данный алгоритм допускает полную автоматизацию процесса предварительной обработки данных и построения карт линейных элементов по любым заданным азимутам индикаторов. Выделение преобладающих азимутов линейных элементов аномалий по описанной методике производилось методом 2-х кратного сглаживания исходных данных в скользящем окне в 5°. В результате обработки структуры аномалий были выделены максимумы частот распределения азимутов, по которым составлены пары взаимоортогональных направлений: (3-273°), (26-296°), (35-305°), (40-310°), (64-334°), (70-340°), (84-354°). Средняя расчетная погрешность определения максимумов для выделенных направлений составила (± 3° - ± 5°).

Анализ плотностиых неоднородностей пространственных распределений линейных индикаторов аномалий. На рис.2 представлено пространственное распределение линейных элементов аномалий геофизических полей (признаков индикаторов разломов) выделенных максимумов (26-296°) в границах 120 - 132° ВД и 56 - 59° СШ. По всей площади условные осевые линии увеличенной плотности линейных индикаторов образуют регулярную и квазирегулярную сеть, в виде ромбических структур (доменов) со стороной 120-150км, условно названных доменами первого порядка.

Р»с.2. Пространственное расположение линейных элементов аномалий гравитационного, магнитного полей для выделенного максимума (26296°). 2.03 - 2.07 - зоны увеличенной плотность индикаторов северо-западного направления, 2.04-2.10 -северо-восточного. Окружность - узел ортогональной системы

Построение систем индикаторов азимутами (3-273°). (35-305°), (40-310°), (64-334°), (70-340 ), (84-354°) приводят к следующим результатам: домены систем (26 и 296°) и (35 и 305°) представлены на 75-85% и проявляются в северо-восточной части Алданского щита. По всей площади проявляются системы (3 и 273°), (40 и 310°), (64 и 334°). (84 и 354°). Причем системы (40 и 310°) и (64 и 334°) обладают наибольшей плотностью индикаторов разломов. Все системы имеют регулярную сетку (ромбические структуры - домены) с расстоянием между осевыми линиями увеличенной плотности индикаторов порядка L = (120 - 150) ± 10км, что можно считать пространственным линейным размером аномалий геофизических полей 1 порядка или линейным параметром 1-го порядка неоднородности геофизической среды. Доменами первого порядка условно названы ромбические структуры со стороной L = (120 - 150) ± 10 км. Сопоставление результатов исследований показывает, что ромбические структуры сохраняют принцип самоподобия, т.е. в доменах 1-го порядка наблюдаются структуры 2-го порядка со стороной L, = V~L.

Пространственный параметр L = 140км размера средней формы рельефа автор [Цу-бои, 1982] интерпретирует как граничное значение, начиная с которого горные хребты не могут удерживаться упругостью земной коры, а должны находиться в изостатическом равновесии.

Наличие локального максимума на пространственных периодах L = (110 - 130) и Lj = HL = 60км [Шерман и др.. 2005] в контексте данной работы соответствует характерным размерам доменов 1-го и 2-го порядка соответственно.

Характерный размер домена L = (120 - 150) ± 10 км, полученный методом статистического анализа геофизических полей, согласуется с пространственными параметрами блоковой делимости литосферы для различных регионов Земли. Природа характерного размера может быть связана с региональностью изостазии. а отклонения от среднего значения данного параметра, выявленные для отдельных регионов - особенностями геодинамической обстановки.

Ромбические структуры (домены) образуются как транспресионные в условиях горизонтального сжатия, теоретически модель формирования которых рассмотрена Вудкоком, Фишером [Woodcock, Fisher, 1986]. На рис. 3 представлены системы: ортогональная (84354°) в виде первичной или основной сетки разрывов и диагональная (40-310°), которая могла быть сформирована как сетка разрывов при внедрении индентора. Причем, теоретически азимут простирания должен быть 45°. Так как основная система имеет азимут 84° (-6°), то динамическая модель Р. Вудкока адекватно согласуется с геометрической моделью блоковой делимости Алданского щита (45° - 6° = 39° ~ 40°). В дополнение к модели Р. Вудкока разработана модель, в которой формируются структуры в азимутах (26 - 296°), в виде диагональной взаимоортогональной сетки на первичной (основной) системе разрывов первого порядка.

//7 \

Теоретические значения углов

относительно относительно

ортогональной системы (34-354),

системы ■6 град

26,56935 20,56505

.33,69007 27.69007

45 39

56,30993 50,33993

63,43495 57.43495

71.56505 65,56505

Рис.З. Расчетные значения и взаимное расположение осевых линий увеличенной плотности линейных элементов геофизических полей (индикаторов) ортогональных систем

Система (64 - 334°) являясь дополнительной к системе (26 - 296°), т.к. 26° + 64° = 90°, также формируется как подчиненная система. Близкие системы по азимутам простираний (26-296°) и (35-305°), (84-354°) и (3-273°), (64-334°) и (70-340°) оказываются пространственно совмещенными системами. Продолжение осевых линий систем (26-296") и (32-302°)установлено, что обе системы совмещены по линиям 55.5°, 59° - 59.5° СШ. Такое расположение двух систем может служить косвенным подтверждением ротационной природы структурообразования, разработанной К.Ф. Тяпкиным [Тяпкин, 1990]. С другой стороны северо - западный фланг двух систем с координатами -120° - 122°х58° может бьггь центром вращения Алдано - Станового мегаблока в условиях смены геодинамических условий сжатия - расширения. Аналогичное распределение имеют системы (64-334°) и (70340°), (3-273°) и (84-354°).

Выводы по главе 2

1. Методами статистического анализа азимутов линейных элементов аномалий, гравитационных и магнитных ступеней и других геофизических признаков разломов, выделено сет преимущественных направлений, линейных индикаторов разломов азимутами (3-273°), (26-296°), (35-305°), (40-310°), (64-334°), (70-340°), (84-354°), соответственно.

2. Пространственная привязка индикаторов в выделенных направлениях позволила установить наличие линейных зон их повышенной плотности, которые укладываются в определенные, строго ориентированные в пространстве системы. Каждая система образует доменные структуры с линейными размерами Ь от 120 до 150км

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

4. Сравнение результатов данного исследования с работами предшественников по тектонической интерпретации геофизических аномалий свидетельствует о тектонической природе выявленных закономерностей. Все ориентировки разломов на региональном и локальном уровнях, выделенные раннее различными авторами, подтверждены статистическим анализом распределений линейных элементов аномалий геофизических полей. Дополнительно к известным ориентировкам выделены статистические значимые азимуты, ортогональные основным направлениям.

5. Сопоставление систем индикаторов разломов с особенностями тектонического плана территории показало, что каждая из выделенных систем укладывается в определенные геолого - тектонические структуры. Это позволило установить время заложение систем (64-334°), (70-340°) не позднее раннего протерозоя. Субмеридионалыюе расположение узлов выделенных систем в пределах восточного фланга Алданского щита может свидетельствовать о более раннем (докембрийском) заложении близких по простиранию (сопряженных) систем (3-273°) и (84-354 ) в виде первичной сетки разломов. Системы (26296°), (35-305°) и (40-310°) могли быть последовательно сформированы на первичной сетке разломов.

6. Наличие одной (первичной) ортогональной системы разломов позволяет объяснить происхождение всех существующих систем Алданского щита. Сопоставление результатов данного исследования с аналогичными исследованиями по многим регионам позволило предположить, что, выделенные (основные) азимуты простирания и пространственная периодичность носят глобальный всепланетарный характер.

Глава 3. Статистические пространственно-временные модели сейсмичности

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

водится моделирование сейсмичности в различных пространственно - временных масштабах по трем направлениям: статистический анализ сейсмического процесса в модели временного ряда, поиск периодических компонент сейсмического процесса и изучение пространственной структуры сейсмичности с построением моделей сильных землетрясений.

Для анализа сейсмичности использованы каталоги землетрясений ГС СО РАН и ЯФ ГС СО РАН. Представительность землетрясений за исследуемый период во времени и пространстве, параметры сейсмичности, описание механизмов землетрясений и пространственная структура сейсмичности изложены в Приложении 1 в соответствие с опубликованными материалами в открытой печати. Проверка статистических гипотез и оценивание параметров модели проведено с использованием классических методов математической статистики.

Изучение динамики распределения энергетических параметров сейсмичности Олекмо-Становои сейсмической зоны

Изменение энергии и количества землетрясений во времени. Для моделирования сейсмического процесса использовалась база сейсмологической информации о землетрясениях в Южной Якутии, содержащая данные с 1969 г. по 2005 г. Сейсмичность в рассматриваемом районе анализировалась на основе инструментальных наблюдений, главным образом, над слабыми землетрясениями, число которых за период наблюдений составило более 15 тыс. Информация о землетрясениях содержала следующие данные: год, месяц, число, время, координаты и К= 1вЕ.

Анализ распределения энергии и количества землетрясений по месяцам и часам в пределах суток приводит к следующим выводам: максимум высвобождаемой энергии за исследуемый период времени приходится на июнь и ноябрь; по времени в сутках отмечается три пика в 8, 11, 14 часов (по Гринвичу).

Распределение энергии по площади приблизительно постоянно и не имеет никакой связи с количеством землетрясений. Среднегодовое значение высвобожденной энергии может рассматриваться как нормальный сейсмический фон. Если обозначить через т? -число землетрясений в с энергетическими классами К]в. период регистрации р, то суммарную энергию землетрясений, отнесенную к периоду р, можно выразить в виде уравнения:

>12 „ (3.1)

где г1 и г2 нижний и верхний предел регистрируемых энергетических классов без учета сильных землетрясений, К] = ДК-], ДК=1. Пусть Г длительность анализируемого периода (1 год в данном случае), тогда

(3.2)

Ч = IУь^РУ10^ = согШ = с1-

Ех(Р)=ЕРр,Р = 1,2,3... (3.3)

Так как Т ■ р имеет размерность текущего времени, то формулу (3.3) можно переписать в виде

Ег (/) = ££■/ = С,/, (3.4)

где константа С, имеет смысл видимой (регистрируемой) скорости диссипативных процессов {Дж/год}. Равенство Е? = С[ (формула (3.2)) при закономерном изменении т^(р) означает, что за период изменений происходит перераспределение энергии между низкоэнергетическими и высокоэнергетическими классами землетрясений.

На уровне 10 Дж распределение землетрясений во времени имеет период около 2,5 лет, что немного меньше, чем период повторяемости для землетрясений с энергетическим классом К = 13 для всей сейсмической зоны

Статистические оценки повторяемости землетрясений в модели стационарного Пуассоновского потока событий

Одним из параметров сейсмического режима является изменение графика повторяемости (угла наклона) у и сейсмической активности А во времени. Для изучения динамики сейсмического процесса был использован каталог землетрясений Южной Якутии с 1969 по 1987 гг. Общий объем выборки 4060 значений. График повторяемости рассчитывался по 10-ти летним циклам со смещением на один год. В каждой серии расчетов использовано от 1500 до 2500 событий.

Из расчетов исключены данные 7-го энергетического класса, не являющегося представительным в данном каталоге (г!=8). Верхний предел статистически значимых данных равен 12 (г2= 12). Интегральный параметр сейсмической активности, рассчитанный по формуле. ■ ■ .

, 1 1 „ 1 1 -'^г (3-5)

равен 6-106 Дж/(км"Тод). В формуле (3.5) N] - средне статистическое количество зарегистрированных землетрясений в классе ) е (г1. г шах)за 10-ти летний цикл, Т10 - период измерений (год* 10), 5«200(50) - площадь (км2) ОСЗ, = 1000км2 стандартный размер площади при.расчете сейсмической активности, Е} =10А'(Дж), Ef - суммарная энергия землетрясений ^ -го энергетического класса. По средним значениям величин повторяемости землетрясений рассматриваемых энергетических классов получено уравнение линейной регрессии в логарифмическом масштабе по десятилетним циклам. Это позволило составить ряд соотношений, для прогнозирования сейсмической активности в разлитых вариантах регистрируемых классов землетрясений (экстраполяцией в низкоэнергетическую область землетрясений)::

1В(Л') = 6.81-0.6А\ (3.6)

Л' = 410-ехр(-1.3811(А - А'6)) (3.7)

Е~ =4-108ехр(0.9215(А', -А',)). (3-8)

В различных вариантах аппроксимации величина энергии, которая рассеивается в виде землетрясений в течение года, с учетом микроземлетрясений (от Л=1), составляет ^=1012-10125Дж.

Верхний предел статистического энергетического параметра Етм = 1012 5 Дж можно определить как предельную энергетическую насыщенность упругой геологической среды в стационарной фазе развития сейсмического процесса. В результате расчетов было установлено изменение угла наклона графика повторяемости у во времени по закону

у(/)=-0,0019/2 + 0,0128^+0,5428. (3-9)

где начало цикла » = 1 относится к середине интервала 1969-1978гг.

Наклон [рафика повторяемости для Южной Якутии оказался зависящим от времени в период подготовки землетрясения 1989г. и, следовательно, данный параметр сейсмического процесса может быть использован при анализе скрытых периодов.

Уменьшение наклона графика повторяемости при сохранении суммарной энергии землетрясений (3.2) означает, что в данный период происходило уменьшение количества слабых землетрясений и увеличение количества событий высших энергетических классов. Этот факт может свидетельствовать о сложности формальной интерпретации статистиче-

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

Проявление землетрясений и их фор-афтершоков на фоне стационарного сейсмического процесса. Вопросы исследования закономерностей пространственно-временных параметров сейсмической активности Олекмо - Становой зоны (ОСЗ) связаны с задачами возможного влияния динамических характеристик сейсмо - тектонических процессов на формирование непериодических изменений геофизических полей. Расчеты показали, что суммарная среднегодовая энергия, высвобождаемая в течение года на изучаемой сейсмогеаерирующей области, остается постоянной (соотношение (3.2)) в пределах энергетических классов К = 12 - 12.5 (К = ^Е, где Е в [Дж]). Следовательно сильное землетрясение по анергии должно превышать порог с К " 12.5. Это с одной стороны. С другой стороны само сильное землетрясение должно менять закономерности развития сейсмического процесса.

Рассмотрение миграций эпицентров землетрясений за период с 05.1988 по 08.1990г. показало, что после сейсмического события с К=16.5 все события сгруппировались в области очага. То есть в период действия афтершоков появился единый энергетический центр выделения сейсмической энергии. Сопоставление времени действия энергетического центра с энергией землетрясений позволило составить уравнение регрессии для периода релаксации: Ь§(Е) = 1.51^(Т) +• 12.63, где Е в [Дж], Т в [суг]. Если обозначить через 1^(Ео) = Ко = 12.63, то формул)' можно получить ряд соотношений:

(3.10)

Г = Т0\О 15 , Г и Г0(£7£0)а667 , Г = Т„ 10А'-М°, М0 = 5.1

Для анализа сейсмической активности было предложено [Трофимснко, 2007] рассматривать временные ряды изменения суммарной энергии землетрясений в интегральном виде, график которого можно описать уравнением линейной регрессии

где иа - скорость дис-сипативных процессов, п = 0,1,2,... Т - текущее время в единицах периода усреднения (рис.4) по сути, представляет собой модифицированный график Беньофа.

Детальный анализ динамики изменения суммарной энергии по всей ОСЗ, позволил выявить ряд закономерностей развития сейсмического процесса. В период 1981 - 1989гг. угол наклона уменьшился (уменьшилось количество высвобождаемой энергии упругих деформаций) и в 1989г. произошло Южно-Якутское землетрясение. В период 1999-2005гг наблюдается аналогичное изменение и как следствие ноябрьское землетрясение 2005г с М=5.7. В период 1991-1996гг. отмечен всплеск сейсмической активности, сопровождавшейся землетрясениями 11-12 классов. Значимого события в ОСЗ не произошло. Однако в БРЗ зарегистрированы события с М=5 26.04.1994т и М=6.1 21.08.1994, а на восточном об-

у = 0.0263х + 12.641 И2 = 0.9836

8 13

2 §

й-л

о 12.5

19Е4 1965 1974 1979 . 1984 1989 1994 199В 2004

Рис.4. Среднегодовое выделение энергии землетрясений с накоплением по всей ОСЗ (каталог ИЗК). 1 - периоды повышенной сеГкашческоЯ опасности со средней еидимой скоростью диссипации энергии < 10" Дис/год. 2 - периоды повышенной сейсмической активности с Ц, > 1012 Дж/год

рамлении Амурской плиты произошло землетрясение М=7.5 - 27.05.1995. Если рассмотреть график, начиная с 1978г., то можно обратить внимание на то, что время накопления энергии упругих деформаций (уменьшенный наклон) и период повышенной разгрузки среды имеют одинаковую длительность.

Во-вторых, в отличие от энергетических (магнитудных) представлений Беньофа данный график (рис. 4) носит более регулярный характер. Это связано с тем, что в расчет энергии принимаются данные без сильных землетрясений и их афгершоков. Особые точки графика (точки перегиба) связаны с сейсмическими событиями и, в данном случае, могут быть проинтерпретированы следующим образом. После 1967г. (Тас-Юряхское землетрясение) скорость диссипативных процессов оставалась постоянной до 1981г. Величина наклона графика практически совпадает с трендовой прямой, что можно принять за нормальную скорость диссипативных процессов иё. Более детальное рассмотрение приводит к следующему: 1968-1971гг. - ол =4.45х10п Дж/год. В конце периода землетрясение с энергией Е=1015 Дж. Период 1981 - 1989гг с о^ =5.16x10" Дж/год закончился землетрясением с энергией Е=1016"6 Дж. В конце периода 1999-2005гг. с ил =5.16x10" Дж/год в течение месяца произошло два землетрясения общей энергией Е=1014 8 Дж. Окончание периода 1991-1996гг. с =2.93x10" Дж/год отмечено всплеском сейсмической активности, сопровождавшегося землетрясениями 11-12 классов. В периоды повышенной сейсмической активности с видимой скоростью диссипативных процессов ил от 1.1х1012 до ц,

2.2х Ю12 Дж/год не произошло ни одного значимого события в ОСЗ.

Видимая скорость диссипативных процессов вследствие упругих деформаций, проявляющаяся в виде изменений скорости высвобождаемой энергии землетрясений, может служить индикатором готовящегося сильного землетрясения. Если рассмотреть график, начиная с 1978г., то можно обратить внимание на то, что величина временного интервала накопления энергии упругих деформаций (её вариативной составляющей) с и^ <1 .Ох 10й Дж/год больше либо равна периоду повышенной разгрузки среды с >1.0х10и Дж/год.

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

Изучение периодических компонент сейсмического процесса.

Исследование стационарности сейсмического процесса сформированной области афтершоков сильного землетрясения М=7.0. Для исследования параметров сейсмического режима был выбран район Тас-Юряхского землетрясения 1967г с М=7, в пределах юго-западной части ОССЗ, для чего использован каталог землетрясений, насчитывающий более 4 тыс. событий. Для расчетов координаты отсчитывались от геометрического центра полигона: (56.8 с.ш., 121.14в.д.).

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

1. среднее значение широты во времени меняется с периодом порядка 22-лст. С 1971 по 1983гг. отмечен северный дрейф энергетической области. Выделяются квазициклы с периодом от 4 до 5 лет; график дисперсии показывает увеличение абсолютного значения, что означает расширение области и содержит квазипериодичность.

2. среднее значение долготы во времени меняется с периодичностью 11 лет; за это же время область увеличилась по долготе (максимум 1982-1983гг.). Общий период равен, по-видимому, 22 года;

3. среднее значение энергетического класса имеет минимум в 1978 - 1979гг. Разброс значений дисперсии при закономерном уменьшении общей энергонасыщенности области означает перераспределение энергии между землетрясениями различных классов.

Рис.5. Динамика параметров сейсмического процесса области Tac - Юряхского землетрясения N1=7. А - изменение широты эпицентров; Б - изменение энергетических классов; В - изменение долготы эпицентров; Г - изменение суммарной энершн; Д - пространственное распределение и вектора смещения эшщетров землетрясении исследуемой области; Е - изменение количества землетрясений

4. отмечается абсолютный минимум энергетической активности области в 1988г. (рис. 5Г, Е). предшествующий Южно Якутскому землетрясению М=6.6 1989г.

На основании проведенного анализа можно заключить, что суммарная энергия землетрясений для данного участка не является стационарным процессом и может быть использована в прогнозных целях. Однако причина не стационарности может быть связана, во первых, с длиннопериодной миграцией сейсмической активности за пределы анализируемой области и, во вторых, с циклическими процессами, обусловленными процессами во всей анализируемой области.

Закономерные изменения параметров сейсмичности области Tac - Юряхского землетрясения позволяют определить динамические характеристики сейсмического процесса: эпоха 1971 - 1972гг. является началом цикла сейсмической активности длительностью 22 года; эпоха 1983 - 1984гг. экстремальный период развития сейсмического процесса, от-

стоящий на 11 лет от начала цикла; 1981 - 1984гг. период неустойчивого периодического развития сейсмического процесса (рис. 5Е); в течение выделенного полупериода сформировалось поле распределений эпицентров землетрясений в виде ромбической структуры; скорость миграции эпицентров землетрясений вдоль долготы (изменение широты) равна ч& =4 км/год. Направление дрейфа показано на рис. 5Д; в короткопериодных 3-х и 5-ти летних циклах скорость миграции эпицентров землетрясений вдоль долготы (изменение широты)равна =10-30 км/год, вдоль широты- «10 км/год.

Изучение суточных и годовых распределений сейсмических событий.

Распределение количества землетрясений в течение суток. Если принять независимость распределения очагов землетрясений, то логично было бы предположить, что за длительный период наблюдений средняя величина энергии и количества землетрясений в течение суток должны стремиться к равномерному распределению, как равновероятные события. Статистические оценки аномальных уровней ряда за 15-ти летний период с 1970 по 1985гг показали наличие трёх пиков в 8, 11, 14 часов по Гринвичу в течение суток. В результате анализа, который был проведен за 1982 -1994 гг. и отдельно за 1989г., т. е. для афтершоков Южно Якутского землетрясения, оказалось, что имеются закономерности распределения количества событий во времени в течение суток как для всей ОССЗ так и для отдельного сейсмического домена. Основной максимум событий приходится на 13 -15 часов иТ.

Было изучено пространственное распределение эпицентров землетрясений в особых точках суточной активности [Трофименко, 2008]. Для центрального участка ОСЗ (зона Нб, рис.1), с известным пространственным положением активных разломов [Рогожин и др., 2007], выделено 3 максимума суточной сейсмической активности. Землетрясения, локализованные вблизи каждого максимума, укладываются в выделенные активные структуры, что дало основание полагать, чго активные разломы не только формируют поле сейсмичности, но и определяют время происхождения (в статистическом смысле) землетрясений.

Аналогичные особенности в распределениях эпицентров землетрясений обнаружены и для других сейсмоактивных зон (глава 5), что дало основание предположить, что поле сейсмичности на энергетическом уровне 10® Дж и ниже формируется в векторном поле внешних сил.

Распределение количества землетрясений в течение годх Анализ распределения энергии и количества землетрясений по месяцам в течение года по статистике распределения землетрясений за 1970-1985гг показал, что максимум высвобождаемой энергии за исследуемый период времени приходится на июнь и ноябрь. Максимумы количества землетрясений с учетом 7-го энергетического класса укладываются в интервалы «март-май», «июнь-август», «октябрь-декабрь» и аналогичны статистическим распределениям представительных классов (9-11 класса). Сопоставление аномалий годовой сейсмичности с изменениями угловой скорости вращения Земли указывает на наличие корреляции двух процессов. В сезонных колебаниях отклонений длительности суток (5Р) отмечается ряд закономерностей: скорость вращения Земли бывает наименьшей б апреле и ноябре, а наибольшей - в январе и июле. Январский максимум значительно меньше июльского. Разность между минимальной величиной отклонения длительности земных суток от эталонных в июле и максимальной в апреле или ноябре составляет 0,001с [Сидоренков, 2004].

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

Первые результаты по выявлению периодических явлений в сейсмическом процессе были получены при изучении динамики афтершоковой области Тас-Юряхского землетрясения 1967г. Был установлен длиннопериодный тренд миграции широты и долготы эпи-

центров и определен вектор смещения. Для суммарной энергии был выделен полупериод длительностью 11 лет (рис.5).

Исследованием миграционных потоков сейсмичности для всей Олекмо - Становой зоны были установлены вектора смещения энергетических центров, периоды изменения сейсмической активности для различных энергетических уровней сейсмичности. Основным периодом изменения сейсмической активности в обоих случаях был определен 22-х летний цикл. Для выявления скрытых периодичностей сейсмического процесса Олекмо -Становой зоны использовался комплексный подход на основе применения различных математических методов, включающих, как традиционные статистические методы [Любу-шин и др., 1998], так и разрабатываемые методы кластерного и др. виды анализа, методика которых разработана применительно к сейсмическому процессу. Результаты расчетов показали, что максимум 22-х летнего цикла приходится на 1983-1984 годы.

Кластерный подход к анализу данных основан на изучении пространственно - временных характеристик дискретных временных рядов в виде определения свойств объекта по некоторой совокупности данных, объединенных по некоторому признаку. Афтершоко-вый сейсмический процесс можно выделить в отдельный кластер по признаку порядкового группирования [Сташевский, 2000] методика расчетов которых основана на вычислении парного корреляционного интеграла [Гусев, 2005] и для которого справедлив закон Омори [Оьшка М., 1987].

Кластерный подход позволяет отбраковывать сейсмические события из общего потока событий, не участвующие в процессах определенной закономерности. Это дает возможность более дифференцированно (и, следовательно, более точно) выделять скрытые периодичности. На рис.1 представлены поля распределений эпицентров землетрясений переходных областей литосферных плит, где показаны направления миграции сейсмической активности. Начало сейсмического цикла для ОСЗ относится к 1972г с максимальным смещением северу в 1983г. В области Нб (рис.1) миграция эпицентров происходит вкрест основных активизированных разломов, имеющих субширотное простирание. В области 11а вектор смещения контролируется системой скрещивающихся разломов меридионального и субширотного простирания.

Третий способ построения кластеров основан пространственной привязке сейсмического процесса к активным тектоническим структурам. В этом случае область динамического влияния разлома рассматривается в виде пространственного элемента, любые сейсмические события в котором относятся к его геометрическому центру. Последовательная активизация сейсмичности в виде обобщенного порядкового группирования [Гусев, 2005] в отдельных кластерах позволяет выделять периоды активизации событий по сейсмическим событиям различных энергетических классов.

Кластерный анализ динамики сейсмичности позволил установить ряд закономерностей: в миграционных процессах участвуют не все землетрясения; землетрясения одного кластера образуют последовательные временные цепочки - порядковое группирование. Это осложняет интерпретацию результатов моделирования. Тем не менее, установленным фактом является северный дрейф эпицентров с 1972 по 1983гг [Трофименко, 2007]. Данный вывод согласуется с результатами, в представленном авторами [Хаин, Халилов, 2008] обзоре но другим регионам. Суммарный эффект от суточных, годовых и в среднем 22-х летних циклов сейсмического режима составляет до 80% всех землетрясений в стационарной фазе сейсмического процесса.

Изучение пространственной структуры сейсмичности и динамические модели сейсмического процесса

Пространственная структура поля сейсмичиости. Для уточнения места возможного землетрясения была построена карта эпицентров землетрясений, на которую вынесены эпицентры землетрясений с размером, пропорциональным классу точности определения очага, т.е. эпицентру присваивается площадь, пропорциональная погрешности определения очага. Фактически это означает, что тем выше класс землетрясения, тем меньшим

значком он фиксируется (инверсное построение). В результате такого построения были прослежены некоторые закономерности: акгивизированые структуры имеют азимуты простирания от 40° до 45° и ортогональные им; сгущения эпицентров землетрясений расположены в узлах ортогональной системы. Крупнейшие сейсмические события либо формируют данную систему, либо формируются в вершинах ромбовидных структур. В последнем случае имеется возможность прогноза мест будущих событий: в пограничных местах скрещивающихся структур. Построение карты афтершоков Южно - Якутского землетрясения по данной методике приводит к аналогичным результатам. Основное событие произошло в вершине ромбической структуры. 99% афтершоков ограничено этой структурой. На рис 4.Д пространственное распределение очагов землетрясений за 1968 - 1994гг. вблизи эпицентра Тас - Юряхского землетрясения также имеют ромбовидную структуру, которая была также сформирована в процессе афтершоковой деятельности. Основное событие произошло на восточной оконечности данной структуры.

Геофизические процессы в переходных областях сейсмоактивных зон. Для установления закономерностей сейсмического процесса переходной области к рассмотре-

переходной зоны от БЗР к ОС.З

Звездочками обозначена трассировка миграции сейсмической активности в пределах западного фланга Алданского тцита.

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

Вследствие этого можно не учитывая реологию определить среднюю скорость распространения деформационного процесса, для которой получена оценка Vj = 6 км'сут. Для быстрой компоненты (событие 1970г.) скорость деформационного процесса составила 40 км/сут.

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

фланге Олекмо - Становой сейсмической зоны в области сопряженных систем разломов Тан Лу, также существенно меняющих сейсмический режим ОСЗ.

Модельные оценки баланса энергии очага и расчет коэффициента полезного действия землетрясения. Проявление динамики сейсмических процессов в геофизических полях означает, что не вся энергия упругих деформаций выделяется в виде сейсмических волн основного события.

Часть запассшюй энергии преобразуется в другие виды тепловой, магнитной, электромагнитной и др. природы. Результаты геофизического мониторинга очаговой зоны Южно - Якутского землетрясения показали, что часть энергии упругих деформаций ре-лаксировалась 4 главными афтершоками (ДЕа = 2-1016 Дж). Около трех тысяч афтершоков высвободили энергию ДЕ3 = 3000(108 -г Ю10) » З Ю13). В электромагнитное излучение преобразовалась часть энергии, равная ДЕ< = 1,2 (109 - 10й) Дж. Однако большая часть энергии, по-видимому, преобразовалась в тепловую энергию. Расчетная температура в очаге оказывается равной 0.092-п, где п в см (температура на 1 см сдвига/с), т.е. порядка ДГ = 300°С, а общая выделившаяся теплота - 1.6101бДж, что по порядку совпадает с расчетной теплотой принятой атмосферой земли. Таким образом, самые общие построения показывают, что оценка КПД сейсмической энергии очага (КПД землетрясения), рассчитанная по формуле

где Е„ - - суммарная, а Е* - кинетическая энергия соответствующей группы

/Г=1

волн, ДII1, = Ер + Еа +Е„г + ЗЕ, - общая энергия, для Южно - Якутского землетрясения

не превышает 0.1%., соответственно Тас - Юряхского в данных предположениях - 1%.

Выводы по главе 3

Моделирование сейсмических процессов Олекмо - Становой сейсмической зоны (ОСЗ), проведенные расчеты и графические представления позволяют сделать следующие выводы:

• Под элементарной сейсмогенерирующей областью будем понимать область, обладающую свойством постоянства выделяющейся сейсмической энергии за некоторый (фиксированный, обычно ежегодный) интервал времени. В расчет принимается минимальная область, обладающая данным свойством. При расчетах суммарной энергии принимаются к расчету все сейсмические события за исключением сильных землетрясений и их афтершоков. Дополнительным ограничением для формирования массивов данных служит представительность энергетических классов, которая должна быть одинаковой для всей анализируемой области. Среднегодовая энергия, высвобождаемая в данной области без учета влияния энергии сильных землетрясений, принимается за нормальный сейсмический фон Еф. Любое сейсмическое событие в данной области с энергией меньше Еф не нарушает стационарность сейсмического процесса. При этом установлена предельная энергетическая насыщенность упругой геологической среды в стационарной фазе развития сейсмического процесса, равная верхнему пределу статистического энергетического параметра Е^ = 10125 Дж. В течение периода изменений происходит перераспределение энергии между низкоэнергетическими и высокоэнергетическими классами землетрясений.

. Рассчитан наклон графика повторяемости по 10-ти летним циклам. Для Олекмо -Становой зоны данный параметр оказался зависящим от времени в период подготовки Южно - Якутского землетрясения 1989г. и, следовательно, данный параметр сейсмического процесса может быть использован для долгосрочного прогнозирования катастрофических землетрясений.

о Изучение миграции очагов землетрясений области Тас - Юряхского землетрясения позволило установить эпоху начала нового 22-х летнего сейсмического цикла, длительность периода и направленность миграции эпицентров землетрясений.

. Землетрясения с энергией более 1013Дж формируют динамическую сейсмическую брешь, пространственные параметры которой охватывают единую структуру геодинамической системы. После ощутимого (значимого) сейсмического события вся область представляет собой зону сейсмического затишья в течение некоторого времени, определяемого магнитудой (энергией) землетрясения и скоростью диссипативных процессов.

• Пространственно-временной анализ изменения широты и долготы эпицентров землетрясений как очаговой так и рассеянной компонент сейсмичности показывает, что миграция эпицентров (энергетических центров) землетрясений происходит со скоростью в среднем равной 10-15км/год. Более динамичной является широтная миграция (периодическое изменение долготы) на фоне длиннопериодных изменений широты на энергетическом уровне 108 - 109Дж. Для землетрясений с энергией Ю10Дж выделен полный цикл длительностью около 22 лет. Скорость миграции сейсмической активности также составляет 12 - 15км/год.

• Перед началом периода повышенной сейсмической активности средняя видимая скорость диссипации энергии для ОСЗ уменьшается до У^ < 1012 Дж/год. На графике суммарной ежегодной энергии землетрясений начало периода повышенной сейсмической активности фиксируется в виде точки перегиба кривой. Длительность периода повышенной сейсмической опасности (Уд < 1012 Дж/год) равна или превышает длительность периода повышенной сейсмической акгивности (У^ 1012 Дж/год).

■ Сейсмический процесс Олекмо - Становой сейсмической зоны развивается под действием региональных напряжений плейттектонической природы, на что указывает одновременное повышение сейсмической активности северо-западной и юго-восточной флангов Амурской плиты 1971-1973, 1995 и 2005 - 2007 гг.

• Количество землетрясений по часам в течение суток распределено не равномерно и имеет характерные участки (пики), отличающиеся для центральной и западной частей ОСЗ. Пространственное распределение эпицентров землетрясений в характерных точках суточной активности показывает наличие упорядоченных структур, укладывающихся в систему активных разломов. В рамках аддитивной модели установлено, что в трех циклах сейсмической активности (суточных, годовых и 22-х летних) участвуют до 80% от общего числа землетрясений. В случае смены геодинамической обстановки «сжатие - расширение» на изучаемой территории может произойти землетрясение с магнитудой М>8, при этом КПД Тас-Юряхского землетрясения составляет около 1%.

Глава 4. Геофизические исследования переходных зон взаимодействия лито-сферных плит

Детальные геофизические исследования в зонах активных разломов. Одной из основных задач современных геофизических исследований в сейсмически активных регионах является изучение геологических предпосылок возникновения землетрясений разной силы и составление региональных прогнозных карт сейсмической опасности. Несмотря на определенные успехи, достигнутые в этой области за последние десятилетия [Геодинамика ...1985; Живая тектоника ..., 1966; Имаев, 2000; Николаев и др., 1982; 2000], корректность используемых критериев сейсмичности остается до сих пор остро дискуссионной. Морфоструктурный и геодинамический анализ показывает, что активное развитие морфострукгур, характеризующихся различной степенью дифференциации и интенсивности тектонических движений, приводит к накоплению напряжений в крупных объемах земной коры [Солоненко, 1966]. Разрядка этих напряжений в виде сильных землетрясений (М=6 и более), хотя и возможна в пределах активно развивающихся блоков, обычно связана С подвижками по зонам активных разломов. Это является решающим при определении сейсмического потенциала не только структурных швов, но и морфострукгур, развитие которых они контролируют. Таким образом, в задачах по детальной оценке степени сейсмической опасности для данной территории одно из важнейших мест занимает геологический анализ активности зон сочленения главных морфоструктурных элементов на протяжении последних нескольких тысяч лет. Актуальность их диктуется нарастающей

динамикой природных и техногенных катастроф и интенсификацией освоения природных ресурсов. На рис.7 представлена карта детальных геофизических работ.

Рис. 7. Активные разломы и сейсмичность Южной Якутии (Овсюченко, Трофименко и др., 2006). Условные обозначения: 1.2-акгивные разломы (1-основные; 2-второстепенные): 1-Кабактинский; 2-Нижне-Нерюнгринский; З-Беркакнтсыш; 4-Суннапщо-Лар&ттск1ш; 5-Южно-Якутсиш; б-Вчжне-Гонаыский; 7-Северо-Становои; 8-Южно-Становой (северная ветвь): 9-Южно-Становои (основная ветвь); 3-мезозойские отложения Южно-Якутской системы впадин; 4-места геофизических исследований активных разломов; 5-пункта геодинамического профиля

Основные результаты геофизических исследований заключались в следующем [Овсюченко и др. 2009].

Южно-Якутский разлом. По данным геофизических исследований (ВЭЗ) в подножии приразломного уступа отмечается появление полого падающих под уступ зон с пониженными кажущимися сопротивлениями, отождествляемых с тектоническими нарушениями. Геоэлектрический разрез на данном участке отражает строение подножия уступа в виде фрагментарно горизонтально слоистой среды (рис.8).

Рис. 8. Электроразведочный профиль, пройденный на меридиональном отрезке зоны ЮжноЯкутского разлома ъ левом боргу р. Горбылях.

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

Беркакитский разлом. По данным геофизических исследований (ВЭЗ) в зоне разлома выделяется падающая на север под углом 70-60° зона повышенного кажущегося сопротивления. В подножии уступа зона выходит практически под дневную поверхность. Ниже по склону расположена еще одна подобная аномалия, выделяющая также и на сейсмическом разрезе. Здесь происходит резкое увеличение мощности приповерхностного высоко-омного (низкоскоростного) слоя. Таким образом, амплитуда вертикального смещения в предполагаемых рыхлых отложениях составляет 8-10 м.

Нижне-Нерюнгринская зона разломов. По данным геофизических исследований (ВЭЗ) разлом выделяется вертикальными зонами резкого увеличения кажущегося сопротивления на условных глубинах от 60 до 40м. С северным приподнятым бортом разлома коррелирует нижняя граница приповерхностного слоя, которая в данном интервале испытывает резкое вертикальное смещение с амплитудой до 20 м и ассоциируется с тектоническим уступом.

Северо-Становая разломная зона. Зона разлома исследована геофизическими методами в районе руч. Половинка. По данным сейсморазведки (КМПВ) наблюдается интерференция волн, что связывается с резким воронкообразным (двусторонним) погружением преломляющей границы, ассоциируемым с зоной дробления выходящей в чехол молодых отложений. Зона дробления выделяется также и по результатам электроразведочных работ.

Южно-Становая система разломов. По данным электроразведочных исследований (ВЭЗ), выполненных в зоне разлома, огранивающего с севера Иенгрскую впадину, выделяется резкая, крутопадающая на юг зона пониженного кажущегося сопротивления. К югу от этой зоны появляется слоистый разрез, выявляемый по появлению слоев с повышенным сопротивлением. В северном крыле разлома слой с повышенным сопротивлением залегает в верхней части разреза. По результатам тектонофизического моделирования на основе повторных геодезических измерений установлено, что описываемый участок расположен на границе участков близгоризонтального сжатия и растяжения. Кроме того, здесь была выделена локальная положительная аномалия скоростей вертикальных движений. Причем южнее северного ограничения Иенгрской впадины наблюдается стабильная отрицательная аномалия скоростей. Севернее в эпоху 1967-1972 гг. выделялась зона отрицательных значений скоростей. В 80-х годах значения скоростей были близки к нулевым. Следовательно, на описываемом участке, характеризующемся рассеянием скоростей вертикальных движений реперов, можно предположить современную тектоническую активность структур, разделяющих блоки с разными скоростями движения. В целом, характер горизонтальных смещений реперов свидетельствует о субмеридиональном горизонтальном сжатии в пределах изучаемого участка. По данным элекгроразведочных исследований (ВЭЗ), выполненных в зоне разлома, огранивающего с севера Верхне-Тимптонскую впадину, выделяется контрастная, крутопадающая зона пониженного кажущегося сопротивления. Этот участок характеризуется относительно низкоомным включением в разрыв вы-сокоомных градиентных зон. На одном участке резко (на 70 м) увеличивается мощность образований с повышенным сопротивлением, которые можно интерпретировать как рыхлые, водопроницаемые отложения пролювиального происхождения.

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

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

Целевые объекты - разрывные нарушения - однозначно выделяются на геоэлектрических разрезах контрастными аномалиями понижения сопротивления пород. Исключение составляет профиль №5, на котором разрывные нарушения выделяются крайне неоднозначно.

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

Зоны дробления и разуплотнения пород характерные для разрывных нарушении, выделяются на сейсмических разрезах в виде характерных воронок. Показательным в этом отношении является разрез по профилю 4 - в дополнение к основным методам по этому профилю были проведены георадарные исследования, результаты которых, в виде волнового разреза, совмещены с основными данными. На разрезе совпадают абсолютно все аномалии - как от рззрывных нарушений, так и обусловленные иными геологическими факторами - и это несмотря на полную независимость исследований и разные физические осноеы методов. Многоактная активизация устанавливается для системы Южно-Якутских разломов.

Комплексная интерпретация результатов геофизических исследований зон активных разломов. Оценка сейсмической опасности данного региона проведена с использованием комплексных сейсмотектонических исследований, которые позволили в региональном масштабе выявлять взаимосвязи сейсмичности с особенностями геолого-структурного и тектонического плана территории. Сейсмический цикл для 8-ми бальных, землетрясений определен в 200 лет [Козмин, Имаев, 2000]

Историческая сейсмическая активность изучаемого региона подтверждается детальными палеосейсмо геологическими исследованиями. Наличие многочисленных дислокаций, свидетельствуют о том, что в Южной Якутии, на протяжении голоцена, не раз происходили сильные землетрясения. Полученные радиоуглеродные датировки палеопочв, погребенных в пределах сейсмодислокаций, позволяют в первом приближении установить возраст некоторых древних сейсмических событий. По этим данным, наиболее уверенно выделяются два палеоземлетрясения [Овсюченко, Трофименко и др., 2007, 2009].

Первое такое событие, вероятно, произошло не ранее рубежа 1540 ± 200 лет назад, согласно радиоуглеродным датам или интервала 1280-1638 лет назад, согласно калиброванному возрасту палеопочв. Очаг землетрясения увязывается с тектоническим узлом, образованным сопряженной системой Южно-Якутского, Суннагино-Ларбинского, СевероСтанового (северная ветвь) разломов.

Второе событие, возможно, произошло около 200 лет назад. Амплитуда вертикального смещения по разрыву в зоне Южно-Якутского разлома при этом составила 10-15 см. Столь небольшая амплитуда смещения позволят предположить, что Южно-Якутский разлом при этом землетрясении был активизирован пассивным образом, т.к. следов выхода очага события на поверхность в зоне разлома обнаружено не было. По имеющимся данным, очаг землетрясения не был связан и с зоной Нижне-Нерюнгрииского разлома, т.к. последняя подвижка по этому нарушению произошла не ранее 950±60 лет назад.

Результаты комплексных исследований стали основанием для существенного уточнения имеющихся сведений о сейсмотектонике региона. В результате применения палео-сейсмогеологического подхода к изучению активных разломов были намечены конкретные геологические структуры, ответственные за реализацию сейсмического потенциала. Это, в свою очередь, позволило построить дифференцированную карту зон возможных очагов землетрясений (ВОЗ) Южной Якутии.

На основании модельных построений авторов [Имаев, Козмин, 2000] и. учитывая, что сильное Учурское землетрясение произошло в 1895г., с учетом радиоуглеродных датировок центральная часть Олекмо-Становой сейсмической зоны находится в середине сейсмического цикла.

Выводы по главе 4

Впервые проведенные геофизические исследования зон активных разломов в переходной области взаимодействия литосферных плит (Амурской и Евроазиатской) показали высокую эффективность данного метода при комплексной оценке энергетического потенциала тектонических структур

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

Результаты геолого-геофизических исследований свидетельствуют о неоднократной тектонической активизации центральной части ОСЗ. Данные результаты являются существенным дополнением к характеристике активных разломов, определенных по сейсмологическому признаку. Некоторая осторожность интерпретации результатов связана с повсеместным распространением мерзлых пород, физические свойства которых изменяются в течение года [Фролов, 1998]. Поведение грунтов при этом определяется не только приповерхностными геологическими условиями, но и более глубинными неоднородностями среды [Джурик и др., 2000].

Глава 5. Модели физических полей и сейсмических процессов в блочной структуре земной коры

Дистанционный геофизический мониторинг активны! тектонических структур

Геофизический мониторинг геологической среды - одно из основных направлений геодинамических исследований, который позволяет дистанционно отслеживать процессы в земной коре. Результаты изучения физических полей, продуцируемых литосферой в связи с сейсмичностью, показывают, что изменение геодинамических условий в переходных зонах литосферных плит отражается в вариациях магнитного, гравитационного, электрического и электромагнитного полей. Масштабы проявления аномалий физических полей варьируются от локальных, при изменении сейсмического режима в близи пунктов наблюдений, до глобальных, при формировании очагов катастрофических землетрясений (глава 1).

Изучение вариаций физических полей в связи с сейсмичностью показывает, что пространственный период проявления аномалий намного превосходит область, охватываемую экспериментальной базой. В качестве классического примера приводится временной ход разности среднесуточных значений магнитного склонения на двух японских станциях, удаленных от эпицентра Ханкайского землетрясения 1946 г. на 60 и 600 км. По данным Като [Ка1о, ШавЫго, 1948] разность вначале - перед землетрясением, примерно на 1 минуту возросла, потом резко уменьшилась на 4 минуты и, далее, в течение восьми месяцев шло восстановление до нормального уровня. Большая длительность заключительной фазы эффекта, по мнению Ногато Т.[ Иа§а1а, 1969]. противоречит пьезомагнитной модели наблюдавшейся вариации, но эффект может быть объяснен как результата движения земной коры или как проявление электрокинетических явлений, возникающих из-за диффузии грунтовых вод. Отклик геомагнитного поля на изменения сейсмической активности был обнаружен и в его вековых вариациях, обусловленных дрейфом главного магнитного поля Земли [Б^Иа, 1965; Тагтт, 1968, 1976; УатагаИ, 1970]. Установлено, что области с ано-

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

В работах [Брейс и др., 1976; Гусева, Добровольский, 1984; Гохберг и др., 1S83; Добровольский, 1984, 1986; Зубков, 1983, 2002, Мячкин, 1974; Сидорин, 1979, 1972] приведены результаты расчета горизонтальной составляющей АН и экспериментальные значения абсолютной величины ДТ геомагнитного ноля по данным наблюдений. Значения ДТ представлены в зависимости ог эпицснтрального расстояния, измеряемого в едииицах радиуса очага землетрясения. Следует заметить, что Нагата [Nagata, 1972] для объяснения сейсмомагнитного эффекта привлекает иные подходы, что, скорее всего, подчеркивает необходимость учета конкретной геологической обстановки для выбора тех или иных моделей сейсмомагнитных эффектов.

Пионерные работы A.A. Воробьева, его теоретические разработки, лабораторные и полевые эксперименты по изучению сейсмоэлектромагнигных явлений позволили установить физическую природу электромагнитного излучения (ЭМИ) горных пород в их естественном залегании в условиях изменяющегося деформационного процесса земных недр [Воробьев, 1971, 1980]. В этой связи, представляют несомненный интерес результаты натурных наблюдений за аномальными вариациями импульсного электромагнитного поля Земли (ЭМИ - ИЭМПЗ) в сейсмоактивных регионах [Соболев ГА., Демин, 1980].

Полный анализ исследований наблюдений электромагнитного поля (ЭМИ) приведен в работах [Гершензон и др., 1986; Гохберг и др., 1983, 1986, 1988]. Первые исследования аномальных возмущений в ионосфере было изложено Я.Г. Бирфельдом в 1974 г. Кривые вариации по ЭМИ имеют квазисинусоидальный характер с максимумом интенсивности излучений, приходящимися па дневное время суток и минимумом интенсивности в ночное время. Характерная цикличность нарушается за несколько часов до реализации наиболее сильных землетрясений. Анализ материала позволяет сделать вывод о возможности выхода на поверхность электромагнитных эмиссий и дает авторам основание утверждать, что существует принципиальная возможность обнаружения предвестников электромагнитных аномалий. Аномалии ЭМИ проявляются на больших площадях и расстояниях от очага землетрясения. Причем, в отличие от аномалий магнитного поля, в данном случае отсутствует корреляция между амплитудой аномалии и параметрами очага.

Изучение неприливных изменений силы тяжести и показывает их высокую связь с современными сейсмотектоническими процессами [Буланже, 1983]. Отмечено, что неоднородность земной коры может играть существенную роль в весьма быстрых современных процессах, происходящих в земных недрах. Планетарные вариации силы тяжести оценивались при условии, что они обусловлены «таянием» земной коры. Расчетный эффект (систавляет около 0,004 мГал в год. Вулканогенные и сейсмогенные изменения Аg обычно связывают с вертикальными перемещениями земных блоков или с изменением влагонасыщенности горных пород. Выявление указанных эффектов сопряжено с рядом методических трудностей, обусловленных большой чувствительностью гравиметров к изменению внешней температуры, атмосферного давления, а также, к микросейсмам и вибрациям. Отмечено наличие взаимосвязанных изменений Ag и современных движений земной коры [Курскеев и др., 1981]. Однако, анализ показывает, что вертикальные микро движения сами по себе практически не могуг обеспечить наблюдаемые изменения &g. В подобных случаях гравитационный эффект может бьггь обусловлен различием изменений состояния и свойств горных пород (фактически изменением плотности). Из других причин изменений силы тяжести можно выделить эффекты влияния на величину силы тяжести аномальных масс атмосферы [Бровар, 1980; Гриднев и др., 1986].

Собственные исследования автора на геодинамическом полигоне Южной Якутии, проложенном в крест Становых разломов (рис. 8), показали, что аномалии силы тяжести не связаны с конкретным сейсмическим событием, а могут быть обусловлены региональными причинами сейсмотектонической природы. Аномалии силы тяжести выявлены на отдельных пунктах профиля протяженностью в 120км. Расположение пунктов регистра-

ции показано на рис. 8 [Трофименко, 1990; Гриб, Трофименко, 2003]. Очаг землетрясения 1983г расположен в 80км севернее п. Чульман. Для отдельных пунктов выделены аномалии различной формы и длительности с практически одинаковой амплитудой. На основании данных исследований было предположено, что в период формирования очага землетрясения исследуемая область находилась под действием горизонтальных сил субмери-диопального сжатия. В данном случае форма аномалии будет определяться физическими процессами в зонах разломов, амплитуда аномалии будет зависеть от величины избыточных напряжений сжатия и не будет зависеть от энергии землетрясения [Трофименко и др., 1995].

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

Интерпретация вариаций ЭМИ в блоковой модели земной коры. Анализ взглядов на природу аномальных проявлений физических полей в сейсмотектоническом контексте свидетельствует, что аномалии глобальной геофизической среды (глава 1) наиболее интенсивно продуцируются в зонах динамического влияния разломов, к которым приурочены эпицентры землетрясений. Максимальным уровнем сейсмической опасности характеризуются зрелые деструктивные поля зоны динамического влияния генеральных разломов. Именно эти области земной коры, наделённые высоким сейсмическим потенциалом (глава 1, приложение Ц следует считать источниками аномалий физических шлей. В этой, связи, важной задачей становится изучение пространственных закономерностей распределения аномалий физических полей.

Исследование вариаций физических полей в Южной Якутии показало, что наиболее чувствительными методами геофизического мониторинга геологической среды являются поля электромагнитной (ЭМИ) природы [Трофименко и др., 1993, 2007]. За время наблюдений отмечено 3 периода повышенной сейсмической активности в ОСЗ: в 1989 г. и 20052007 гг. Аномальные изменения в электромагнитных полях, предваряющие землетрясения с энергией более 1013Дж в эти периоды, были ассоциированы с сейсмотектоническими процессами в земной коре [Трофименко 2007, 2009]. Произошедшее землетрясение 20.04.89 г. с энергией О16,5 Дж позволило сопоставить аномальные суточные вариации ЭМИ с динамикой афтершоковых последовательностей (рис. 9) [Южно - Якутское..., 1993; Афтершоки Южно -Якутского..., 1993]. Изменение напряженности ЭМИ в моменты, предваряющие главное событие 20 апреля 1989 г. и его основные афтершоки 24 и 29 апреля, а также 7 и 17 мая проявилось в виде увеличения амплитуды сигналов в течение нескольких суток до события. Погрешность прогноза времени события составила +2 суток без указания на характер самого землетрясения: магнитуды, энергетического класса,

интенсивности и места события [Трофименко, 2005]. В течение второго цикла измерений в 2004-2007 гг. в ОСЗ отмечено два периода повышения сейсмической активности. В нервом периоде произошло два землетрясения с энергией Е = 1015,7 Дж - 10.11.2005 и Е

1014'8 Дж 11.12.2005, во втором - также два землетрясения с энергией Е = 1014'6 Дж - 03.03.2007 и Е = 1014-2 Дж -13.03.2007.

Во всех четырех случаях имели место аномальные суточные вариации ИЭМГ13. в каждом из двух пунктов измерений [Трофимеико, 2007].

Детальное рассмотрение временных изменений геофизических полей вследствие сейсмотектонических процессов привело к необходимости решения двух задач: поиск корреляционных зависимостей «аномалии ЭМИ - землетрясение» и определение радиуса (сферы) влияния тектонических процессов на формирование геофизических предвестников землетрясений. Сложности интерпретации, которые при этом возникают, связаны с тем, что групп аномальных сигналов ЭМИ-ИЭМПЗ больше чем землетрясений и, следовательно, при наличии явных предвестников землетрясений нет возможности детерминированного прогноза событий. Это указывает на то, что аномалии ЭМИ в условиях современной активной геолого-тектонической структуры не являются достаточным (100%) признаком возможного повышения сейсмической активности. Вероятность наступления события (отношение числа прогнозных аномалий к общему числу зарегистрированных аномалий) оценивается как или 75%. Привлекая модель сейсмического процесса в блоковой среде (глава 1, [Садовский, 2004]), можно предположить, что аномалии ЭМИ должны проявляться дважды: в период консолидации блоков на начальной стадии формирования очага землетрясения и в период потери устойчивости консолидированных систем блоков, т.е. в периоды наибольшего упруго-кинематического взаимодействия между блоками. Таким образом, аномалии ЭМИ, выделенные в 1989 г., могут быть разложены по составляющим: афгершоку 28 апреля соответствует две группы аномалий с 27 апреля по 3 мая; афтершоку 7 мая - две группы аномалий с 4 по 11 мая; афтершоку 17 мая - с 12 по 25 мая. Характер изменения суточных аномалий ИЭМПЗ в 2005 и 2007 годах аналогичен вариациям ЭМИ 1989г.

В работе [Гохберг М.Б. и др. 1991] показано, что даже при трансформации всей энергии очага землетрясения с Е=\()и Дж в электромагнитную энергию напряженность ЭМИ на дневной поверхности будет меньше фоновой составляющей. Следовательно, регистрируемые аномалии являются следствием приповерхностного взаимодействия блоков земной коры при их относительном перемещении вблизи пункта наблюдения и определяются физическими свойствами горных пород, вещественным составом, характером движения по разлому.

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

Рис. 9. Изменения напряженности ЭМИ в апреле-мае 1989 г. в период Южно-Якутского землетрясения и его афтершоков.

Стрелками отмечены моменты сейсмических событий с указанием энергетических классов (К=1§Е, где Е-энертя землетрясения в Дж), скобками выделены серии аномалий с указанием их порядковых номеров

ластях. В работе [Трофименко, 2007] показано, что в период активизации сейсмической активности в Олекмо - Становой сейсмической зоне в апреле-мае 1989 г. аналогичные изменения происходили как в Сейсмическом поясе Черского, так и на Камчатке, а в 19711973гг. - на Сахалине. Таким образом, в апреле-мае 1989 г. мог существовать единый источник, как повышенной сейсмической активности, так и непериодических изменений в геофизических полях. Протяженность зоны активизации, по-видимому, составила более 3000км.

Полевые исследования вариаций электромагнитного излучения (ЭМИ) показали высокую эффективность данного метода геофизических наблюдений для контроля динамики литосферных плит. В частности перед всеми сейсмическими событиями в пределах обрамления Амурской плиты (рис.1) наблюдались аномальные изменения в суточных вариациях ЭМИ в виде отдельных серий, что свидетельствует об их сейсмотектонической природе.

В рамках теории (модели) сейсмического процесса в блоковой среде сделано предположение, что аномалии ЭМИ должны появляться дважды: перед землетрясением и после него, к тому же не обязательно симметрично относительно сейсмического события, так как скорость накопления энергии упругих деформаций и скорость разгрузки среды не одинаковы. Аномальные амплитуды суточных вариаций могут регистрироваться в течение периода не стационарности сейсмического процесса. Результаты мониторинга свидетельствуют о том, что активные тектонические структуры развиваются как единая геодинамическая система, взаимодействующая с окружающими ОСЗ сейсмогенерирующими зонами.

Обобщение статистических моделей на сейсмические зоны северо-востока Азии

Исследование пространственно - временных характеристик землетрясений сейсмических поясов Северо-востока Азии для различных геодинамических условий развития деформационных процессов установлен ряд общих закономерностей, на основании которых выделены области, чувствительные к быстро изменяющимся геодинамическим условиям. Исследование динамики сейсмичности сейсмических поясов Северо-востока Азии показало идентичность распределений количества землетрясений в течение суток и года (рис.10).

Зама На Зона IIa

Рис. 10. Аномалии суточной и годовой активности сейсмических зон Северо-востока Азии

Для областей скрещивающихся тектонических структур в суточных распределениях отмечается 2-3 максимума активности и 1-2 максимума в годовой. По результатам определения максимумов суточной активности для северной границы Амурской плиты в предположении плоской волны установлена фазовая скорость смещения максимума сейсмической суточной активности. Истинное долготное время для выделенных областей составляет: БРЗ - 7.4ч, ОСЗ - 8ч, Сахалин - 9.47ч. Максимумы суточной активности соответственно - 18, 19 и 20ч Ш'. Отсюда истинное время максимумов 25.4, 27 и 29.5 часа ЬТ. Разница между максимумами на Сахалине и БРЗ составляет 2 часа, т.е. относительно нулевого меридиана они формируются одновременно. Отсюда фазовая скорость смещения максимума суточной активности не может бьгть меньше чем ЮООкм/час.

В годовых аномалиях отмечается наличие максимума в ноябре не зависимо от долготы исследуемого участка. Для Арктической зоны один из максимумов активности смещен на январь февраль.

Динамическая модель формирования геофизических нолей и процессов и проблема прогнозирования землетрясений.

В модели блокового строения литосферы в системе гравитационного взаимодействия «Земля - Солнце» в суточных аномалиях сейсмического режима должно проявляться от двух до четырех максимумов сейсмической активности: в максимумах горизонтальной и вертикальной составляющих гравитациошюго потенциала. Причем в максимуме горизонтальной составляющей потенциала за счет увеличения межблокового взаимодействия максимум суточной сейсмичности должен проявляться в утренние и вечерние часы в зонах активизированных разломов с вертикальными межблоковыми шовными зонами (рис.11). Максимумы сейсмической активности в периоды экстремальных значений вертикальной составляющей гравитационного потенциала будут проявляться при наличии наклонных межблоковых шовных зон также за счет увеличения кинематического взаимодействия блоков. В данной модели прогнозируется отсутствие суточных аномалий на линейных субширотных структурах.

Изменение годовой активности в данной модели возможно при формировании ансамбля блоков и при наличии наклонных межблоковых шовных зон. Исходя из инерционной природы годовых циклов сейсмического режима, следовало бы ожидать 4 максимума сейсмической активности. Отсутствие 3-го и 4-го максимума может быть связано с недостаточной статистикой исходных данных, вследствие ограниченности ряда наблюдений либо недоучетом низших энергетических классов (К<7).

50 . 100

1 I I з I I'

/ У S У,

Рис. 11. Динамическая модели земной коры во взаимодействии в системе «Земля - Солнце».

1 - блоки литосферы в модели М.А. Садовского [Садовский, 2004]; 2 - положение фазовой границы в модели Б.В. Левина [Левин и др., 2007]; 3 - шовные межблоковые зоны; 4 - зона эндогенной активизации в результате дегазации [Летников, Дорогокупец, 2001; Гуфельц, 2009]; 5 - направление горизонтальной деформации е положение Солнца (А), соответствующее первому максимуму суточной активности сейсмическою процесса; б - направление деформации в положение Солнца (Б), соответствующее второму максимуму суточной активности сейсмического процесса при наклонном падении шовных межблоковых зон; 7 -межблокоше границы второго порядаа

Моделирование геофизических полей и процессов методами статистического анализа, геофизические исследования и пассивный геофизический мониторинг геологической среды позволяют сформулировать ряд положений методики и технологии организации геофизических исследований переходных зон литосферных плит для контроля геодинамического режима литосферы и прогноза периодов сейсмической активности.

Геофизическим мониторингом геологической среды установлен пространственный параметр влияния сейсмотектонических процессов на формирование аномалий физических полей в пределах активных тектонических структур северной границы Амурской плиты.

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

Не зависимость размеров динамических брешей после сильных землетрясений от положения очага землетрясения (глава 3) внутри области означает, что при наличии геофизических признаков готовящегося землетрясения, не определенность положения очага землетрясения эквивалентна области, в пределах которых снимается видимая часть тектонических напряжений.

Не зависимость размеров геодинамической системы после землетрясений от энергии, начиная с некоторой пороговой величины энергии, даюг возможность прогнозировать именно пороговое (минимальное) значение энергии предстоящего землетрясения (глава 2,4).

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

Исследованием вариаций ЭМИ в зонах активных разломов различных регионах показано наличие двух максимумов суточной активности. Сопоставление данных результатов с изменением суточной активности сейсмического режима (глава 3) может свидетельствовать о едином источнике возмущений для обоих процессов. Таким источником может быть не стабильность суточного вращения Земли. Косвенным подтверждением такой связи может служить результат сопоставления динамики годовой сейсмичности (раздел 3.2) и изменением угловой скорости вращения Земли [Сидоренков, 2001].

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

ЗАКЛЮЧЕНИЕ

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

геодинамических полигонах Южной Якутии, в непосредственной близости от Олёкмин-ского эпяцентрального поля катастрофических землетрясений с магнитудой М=6.5 - 7. В результате проведенной интерпретации суточных изменений электромагнитного излучения в период Южно-Якутского землетрясения обоснована модель аномалий в блоковой структуре литосферы. Данный результат стал стимулом дальнейших исследований по структурированию геофизических полей и сейсмических процессов.

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

Установлено подобие сейсмических процессов для всех сейсмических поясов северо-востока Азии в суточных, годовых и, в среднем, 22-х летних циклах сейсмической активности. Этот результат позволил построить физико-геодинамическую модель сейсмического процесса, адекватно отображающую аномальные и циклические изменения сейсмической активности, а также, направление и скорость пространственно миграции эпицентров слабых землетрясений с энергией Е< 1013Дж.

На основе анализа динамики энергетических параметров сейсмичности произведены оценки скоростей релаксационных процессов в стационарной фазе развития сейсмического процесса и определено пороговое значение скорости, при которой возможно накопление энергии упругих деформаций с последующей разрядкой в виде сильного землетрясения. Данные исследования позволили вывести количественные соотношения между периодами повышения и понижения сейсмической активности и длительности релаксационных процессов в виде афтершоковой деятельности вследствие сильных землетрясений.

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

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

1. Взаимодействие различных геологических оболочек Земли свидетельствует о многоплановом характере реакции одной среды при изменении параметров другой. Геофизические процессы в блоковой модели геофизической среды в каждой из составляющих оболочек глобальной геологической среды обуславливают неоднозначность их взаимодействия. Это в свою очередь может приводить к не равновесному хаотическому состоянию всей геологической среды. Литосфера продуцирует физические поля, которые накладываются на аномальные изменения параметров взаимодействующих систем. Детерминировано - стохастические процессы во взаимодействующих геосферах обуславливает слабую предсказуемость развития сейсмотектонических процессов.

2. Изучение пространственного распределения геофизических полей методом статистического анализа распределений линейных индикаторов аномалий физических полей магнитной и гравитационной природы позволило обобщить результаты геофизических исследований литосферы. На основе рассмотренных материалов подтверждена глобальная общепланетарная структурированность литосферы. Выделенные ортогональные системы линейных индикаторов аномалий физических полей гравимагнитной природы образуют пространственные решетки (домены) с линейными размерами (120-150)±10 км, условно названными системами первого порядка, а домены более высоких рангов образуют решетки с периодами, соответствующими 1/2, 1/4 и т.д. этой величины.

3. Ромбические структуры (домены) в азимутах (40 и 310°) образуются как транспресионные в условиях горизонтального сжатия (модель Вудкока-Фишера). В дополнение к данной модели разработана модель, в которой формируются вложенные структуры во всех выделенных азимутах в виде диагональных взаимоортогональных сеток на первичной (основной) ортогональной системе первого порядка. Взаимное про-странсгвегагое расположение бимодального (3-223° и 84-254°) и тримодального (26-296°, 35-305° и 40-310°) простираний линейных элементов аномалий гравимагнитных полей может свидетельствовать о том, что изменения палеотектонических обстановок и соответствующих ориентировок векторов деформаций на изучаемой территории происходили не единовременно и резко, а постепенно и плавно, что может найти объяснение в неравномерном (циклическом) вращении, или "косой субдукции" Амурской плиты.

4. Изучена пространственная структура и временная динамика сейсмического режима ОСЗ. Получены параметры изменения сейсмической активности как для всей ОСЗ, так и для отдельного афтершокового домена. Установлено, что скорости миграции очагов землетрясений имеют один и то г же порядок.

5. На основании анализа распределений поля сейсмичности всей ОСЗ в периоды сильных землетрясений получены аналитические зависимости для длительности действия афтсршоков.

6. Энергетический подход на основе модифицированной методики Беньофа моделирования сейсмотектонических процессов позволил установить наиболее вероятные периоды возникновения сильных землетрясений. Перед началом периода повышенной сейсмической активности средняя видимая скорость диссипации энергии для ОСЗ уменьшается до < 1012 Дж/год. Длительность периода повышенной сейсмической опасности (ц,< 1012 Дж/год) равна или превышает длительность периода повышенной сейсмической активности (иа > 1012 Дж/год). Сейсмический процесс Олекмо - Становой сейсмической зоны развивается под действием региональных напряжений плейттектошгческой природы, а видимая скорость диссипации энергии о1/ в отдельных частях ОСЗ определяется упруго пластическими свойствами горных пород.

7. Распределение количества землетрясений по часам в течение суток не равномерно и имеет характерные участки, отличающиеся для центральной и западной частей ОСЗ. Пространственное распределение эпицентров землетрясений в характерных точках суточной активности показывает наличие упорядоченных структур, укладывающихся в систему активных разломов.

8. Сопоставлены структурные элементы систем разломов с пространственным полем сейсмичности. На основаши данного исследования установлено их геометрического подобие, как на региональном уровне для всей ОСЗ, так и на локальном уровне по полю афтершоков сильных землетрясений.

9. На основании результатов модельных построений, детальных геофизических работ в зонах активных разломов, мониторинговых наблюдений разработана комплексная геофизическая модель сейсмического процесса ОСЗ, предложена методика анализа сейсмического процесса, установлена неоднократная тектоническая активизация активных тектонических структур, современная активность которых подтверждается мониторинговыми геофизическими наблюдениями, причем разломы различной кинематики активизируются как единая геодинамическая система. По комплексному анализу показана возможность получения состоятельной оценки сейсмической опасности и построения дифференцированных зон возможных очагов землетрясений.

разломов установлена последовательная активизация в течение суток (в статистическом смысле), для сдвиговых деформаций выявлены ромбические структуры, характерные для присдвиговых бассейнов (pull-apart basin); для всех исследуемых зон выделен 22-х летний период миграции сейсмичности с максимумом широты в 1983г., определено пороговое значение энергии землетрясений, меняющее структуру потока сейсмических событий, величина которого находится в пределах 1012 5 - 1013 Дж; наличие динамических брешей затишья выявлено для всех землетрясений с энергией более 1013 Дж, при этом размеры области затишья не зависят от энергии землетрясений и определяют динамическую систему, объединенную общим деформационным процессом плейттектонической природы.

11. На основе результатов анализа пространственного распределения геофизических полей и сейсмических процессов, моделирования результатов геофизического мониторинга геологической среды разработана геодинамическая модель сейсмического процесса в блоковой структуре литосферы. В рамках построенной модели возможна интерпретация основных результатов геофизических исследований.

ОПУБЛИКОВАННЫЕ РАБОТЫ ПО ТЕМЕ ДИССЕРТАЦИИ Статьи в журналах из Перечня ВАК

1. Гриб H.H., Трофимепко C.B. Моделирование техногенной сейсмичности в зоне открытых горных работ Нерюшринского разреза Республики Саха (Якутия) // Горный информационно-аналитический бюллетень/Региональное приложение ЯКУТИЯ. Выпуск -4, 2005. Издательство МГГУ. С. 153-166.

2. Гриб H.H., Трофимепко C.B. Сейсмичность Олекмо-Становой сейсмической зоны инициированная ядерными взрывами // Горный информационно-аналитический бюллетень/Региональное приложение ЯКУТИЯ. Выпуск - 4, 2005. Издательство МГГУ. С. 167-175.

3. Трофимепко C.B. Оценка энергии возможного землетрясения Олекмо-Становой зоны // Горный информационно-аналитический бюллетень/Региональное приложение ЯКУТИЯ. Выпуск - 3,2006.-М.: Издательство МГГУ, С. 149-154.

4. Статива A.C., Трофимепко C.B. Методика и технология выделения систем разломов Алданского щита // Горный информационно-аналитический бюллетень/Региональное приложение ЯКУТИЯ. Выпуск - 3, 2006.-М.: Издательство МГГУ, С. 155-166.

5. Статива A.C., Трофимепко C.B., Гриб H.H. Анализ существующих представлений о тектонике Алданского щита // Горный информационно-аналитический бюллетень/Региональное приложение ЯКУТИЯ. Выпуск - 3, 2006.-М.: Издательство МГГУ, С. 167-188.

6. Трофимепко C.B., Гриб H.H. О возможности прогноза эпицентра сильного землетрясения по структуре сейсмического поля Олекмо-Становой зоны // Горный информационно-аналитический бюллетень/Региональное приложение ЯКУТИЯ. Выпуск - 3, 2006,-М.: Издательство МГГУ, С. 189-192.

7. Имаев B.C., Козьмин Б.М., Гриб H.H., Трофимепко C.B. Анализ геодинамических моделей сейсмических поясов Якутии // Горный информационно-аналитический бюллетень/Региональное приложение ЯКУТИЯ. Выпуск - 4, 2006.-М.: Издательство МГГУ, С. 165-172.

8. Гриб H.H., Трофимепко C.B. О возможности долгосрочного прогноза времени сильного землетрясения по динамике сейсмического процесса Олекмо-Становой зоны // Горный информационно-аналитический бюллетень/Региональное приложение ЯКУТИЯ. Выпуск -4, 2006.-М.: Издательство МГГУ, С. 173-178.

9. Статива A.C., Трофименко C.B. Ортогональные системы разломов Алданского щита и их геологическая позиция // Горный информационно-аналитический бюллс-

тень/Регионалыгое приложение ЯКУТИЯ. Выпуск - 4, 2006.-М.: Издательство МГГУ, С. 179-187.

10. Статива A.C., Трофименко C.B., Имаев B.C. Пространственное распределение индикаторов систем разломов Алданского щита // Горный информационно-аналитический бюллетень/Региональное приложение ЯКУТИЯ. Выпуск - 4, 2006.-М.: Издательство МГГУ, С. 188-196.

11. Гриб H.H., Трофименко C.B. Регистрация электромагнитного излучения (ЭМИ) в период землетрясений // Горный информационно-аналитический бюллетень/Региональное приложение ЯКУТИЯ. Выпуск - 4, 2006.-М.: Издательство МГГУ, С. 197-204.

12. Трофименко C.B. Геофизические поля и сейсмичность Южной Якутии // Горный информационно-аналитический бюллетень/Отдельный выпуск 3, Якутия 1, - М.: Издательство МГГУ, 2007,- С. 188-195.

13. Трофименко C.B. Закономерные изменения параметров сейсмического процесса Олекмо-Становой зоны (ОСЗ) // Горный информационно-аналитический бюллетень/Отдельный выпуск3, Якутия 1. - М.: Издательство МГГУ, 2007 - С. 196-202.

14. Трофименко C.B. Статистические оценки сейсмического режима области Тас-Юряхского землетрясения Олекмо-Становой зоны (ОСЗ) // Горный информационно-аналитический бюллетень/Отдельный выпуск 3, Якутия 1. - М.: Издательство МГГУ, 2007,-С. 213-216.

15. Трофименко C.B., Имаев B.C. Направленность научно исследовательских работ по изучению закономерностей развития сейсмотектонических процессов Олекмо-Становой зоны (ОСЗ) // Горный информационно-аналитический бюллетень/Отдельный выпуск 3, Якутия 1.- М.: Издательство МГГУ, 2007,- С. 203-207.

16. Трофименко C.B. Проявление землетрясений на фоне стационарного сейсмического процесса Олекмо-Становой зоны (ОСЗ) // Горный информационно-аналитический бюллетень/Отдельный выпуск 3, Якутия 1- М.: Издательство МГГУ, 2007,- С. 208-212.

17. Е.А.Рогожин, А.Н.Овсюченко, С.В.Трофименко и др. Сейсмотектоника зоны сочленения структур Байкальской рифтовой зоны и орогенного поднятия Станового хребта.- Геофизические исследования.-М.: Изд-во ИФЗ, 2007. вып. 8, с.81-116.

18. Трофименко C.B., Гриб H.H., Никитин В.М. Геодинамические процессы Олекмо-Становой зоны и вариации геофизических полей.- Наука и образование - №4 (48), 2007,- С .44-49.

19. Карасев П.С., Овсюченхо А.Н., Мараханов A.B., Трофименко C.B. Активные тектонические нарушения участка Алдан-Нагорный нефтепроводной системы Восточная Сибирь - Тихий океан,- Нефтяное хозяйство, 2008. № 9, с. 80-84.

20. А.Н. Овсюченко, С.В.Трофименко, Мараханов A.B., П.С., Карасев, Е.А. Рогожин Очаговые зоны сильных землетрясений Южной Якутии,- Физика Земли. 2009. №2. С. 15-33.

21. А.Н. Овсюченко, C.B. Трофименко, Мараханов A.B., П.С., Карасев, Е.А. Рогожин, B.C. Имаев, В. М. Никитин, Н.Н.Гриб Детальные геолого-геофизические исследования зон активных разломов и сейсмическая опасность Южно-Якутского региона.- Тихоокеанская геология, 2009, том 28, №4, с. 55-74.

22. А. Н. Овсюченко, С. В. Трофименко, А. В. Мараханов, П. С. Карасев, Е. А. Рогожин Сейсмотектоника переходной области от Байкальской рифтовой зоны к орогенному поднятию Станового хребта. - Геотектоника, 2010, № 1, с. 29-51.

23. Трофименко С. В. Тектоническая интерпретация статистической модели распределений азимутов аномалий гравимагнитных полей Алданского щита,- Тихоокеанская геология, 2010, том 29, №3, с. 64-77.

24. Трофименко С. В. Моделирование пространственно - временной цикличности сейсмического процесса на примере Олекмо - Становой зоны.- Наука и образование,- №1, 2010-С.54-58.

25. Трофименко C.B. Статистические пространственно - временные модели сейсмичности. - Вестник ЯГУ, 2010, том 7, №2. - С. 48-55.

Другие публикации

26. Трофименко C.B., Лозовский JI.A. Регистрация лунно - солнечных вариаций донным гравиметром ГДК в автономном режиме //Геофизические исследования в Якутии: Сб. науч. тр. - Якутск: ЯГУ, 1990. - С. 49-53.

27. Трофименко C.B. Высокоточные гравиметрические наблюдения при условии низких температур //Геофизические исследования в Якутии: Сб. науч. тр. - Якутск: Я1"У, 1990. - С. 49-53.

28. Козьмин Б.М. Голенецкий С.И. Николаев В.В., ...,Трофимеико C.B. и др. Южно

- Якутское землетрясение 20 апреля 1989г. - Южно - Сахалинск: ИМГиГ ДВО РАН, 1993. -45с.

29. Козьмин Б.М. Голенецкий С.И. Николаев В.В., ....Трофименко C.B. и др. Аф-тершоки Южно - Якутского землетрясения 20 апреля 1989г. - Южно - Сахалинск: ИМГиГ ДВО РАН, 1993. - 66с.

30. Трофименко C.B., Широков И.А., Анохина K.M. Результаты земноприливных наблюдений на территории Южной Якутии //Геофизические исследования в Якутии: Сб. науч. тр. - Якутск ЯГУ, 1995. - С. 55-63.

31. Трофименко C.B., Статива A.C. Методика и техника стационарных геофизических исследований для целей прогноза землетрясений //Материалы III региональной научно практической конференции, посвященной 10 легию ТИ(ф) ЯГУ, Нерюнгри: 2002. -С.73-78.

32. Гриб H.H., Трофименко C.B. Возможности прогноза землетрясений по результатам геофизических исследований в Южной Якутии в период Южно - Якутского землетрясения 1989 г. //Пятые геофизические чтения имени В. В. Федынского 27 февраля - 01 марта 2003 года. - М.: ГЕОН, 2003. - С. 67-68.

33. Гриб H.H., Трофименко C.B. Неприливные изменения силы тяжести в зонах влияния современных активизированных разломов //Проблемы сейсмологии Ш-го тысячелетия: Материалы международной конференции (Новосибирск 15-19 сентября 2003 г.).

- Новосибирск: Изд-во СО РАН, 2003. - С. 271-274.

34. Трофименко C.B. Модели геофизических процессов в период Южно-Якутского землетрясения 1989 г. //Сейсмичность Южно-Якутского региона и прилегающих территорий: Материалы всероссийской научно-практической конференции 24-27 октября 2005 г.

- Нерюнгри, изд-во ЯГУ, 2005,-С.99-104.

35. Трофименко C.B. Термодинамическая модель Южно-Якутского очага землетрясения // Сейсмичность Южно-Якутского региона и прилегающих территорий: Материалы всероссийской научно-практической конференции 24-27 октября 2005 г. - Нерюнгри, изд-во ЯГУ, 2005,-С. 163-165.

36. Трофименко С. В., Гриб H. Н. Временные макропроявления сейсмичности Олекмо-Становой сейсмической зоны (ОССЗ) //Активный геофизический мониторинг' литосферы земли: Материалы 2-го международного симпозиума 12-16 сентября 2005 г., Академгородок, Новосибирск. - Новосибирск: Изд-во СО РАН, 2005.-С.381-385.

37. Трофименко С. В., Гриб H.H., Никитин В.М., Муллаяров В.А. Результаты наблюдений за вариациями естественного электромагнитного поля Земле в Южной Якутии // Материалы IV международной конференции 14-17 августа 2007 г., «Солнечно-земные связи и предвестники землетрясений»,- с. Паратунка, Камчатской обл.-Петропавловск-Камчатский.-Изд-во ИКИР ДВО РАН, 2007.-С.453-458.

38. Трофименко C.B. Статистические модели сейсмического режима Олекмо-Становой зоны (ОСЗ) // Физика геосфер: Материалы пятого всероссийского симпозиума.-Владивосток.: Изд-во Дальнаука, 2007,- С.218-225,

39. Трофимснко C.B. Проявление землетрясений и их фор-афтершоков на фоне стационарного сейсмического процесса // Материалы всероссийского совещания с международным участием 18-24 сентября 2007г. «Проблемы современной сейсмологии и геодинамики Центральной и Восточной Азии».-Иркутск,- Изд-во ИЗК СО РАН, 2007,- в 2-томах.-Т.2.-С.171-175.

40. Трофименко С. В., Никитин В.М., Статива A.C., Гриб H.H. Системы разломов Алданского щита в связи с сейсмичностью Южной Якутии // Проблемы современной сейсмологии и геодинамики Центральной и Восточной Азии: Материалы всероссийского совещания с международным участием 18-24 сентября 2007г. - Иркутск.: Изд-во ИЗК СО РАН, 2007- В 2-х томах.-Т.2.-С. 176-179.

41. Овсюченко А.Н., Мараханов A.B., Карасев П.С., Рогожин Е.А., Трофименко C.B., Никитин В.М. Зоны возможных очагов землетрясений и сейсмическая опасность Южно-Якутского региона // Южная Якутия - новый этап индустриального развития: Материалы международной научно-практической конференции 24-26 октября 2007, г. Не-рюнгри, Россия. — Нерюнгри,: изд-во Технического института, 2007.- в 2-х томах.-Т.2.-С.46-58.

42. Трофименко C.B. Тектонические модели прогнозируемых очагов землетрясений //Южная Якутия - новый этап индустриального развития: Материалы международной научно-практической конференции 24-26 октября 2007, г. Нерюнгри, Россия. - Нерюнгри,: изд-во Технического института, 2007,- в 2-хтомах.-Т.2.-С.88-91.

43. Трофименко C.B., Гриб H.H., Никитин В.М. Вариации электромагнитного поля как отражение сейсмотектонических процессов Олекма-Становой зоны// Известия Томского политехнического университета. Том314, №1, 2009. «Науки о Земле» - С. 48 -53.

44. Трофименко C.B. Геоморфологические признаки модели поля сейсмичности Олекмо - Становой зоны // Геодинамика внутиконтинентальных орогенов и геоэкологические проблемы: Сборник материалов 4 Международного Симпозиума, 15-20 июня 2008 г.: Москва - Бишкек; 2009, С. 241-245.

45. Трофименко C.B., Гриб H.H., Никитин В.М. Вариации электромагнитного поля как отражение сейсмотектонических процессов Олекмо-Становой зоны // Геодинамика внутиконтинентальных орогенов и геоэкологические проблемы: Сборник материалов 4 Международного Симпозиума, 15-20июня2008 г.: Москва-Бишкек; 2009, С. 137-141.

46. Трофименко C.B. Статистические модели пространственных распределений аномалий гравимагнитных полей и их отражение в структурно - тектоническом строении Алданского щита // Тектоника и глубинное строение востока Азии:VI Косыгинские чтения: доклады всероссийской конференции, 20 - 23 января 2—9, г. Хабаровск/под. ред. А.Н. Диденко, A.A. Степашко. - Хабаровск: ИТиГ им. Ю.А. Косыгина ДВО РАН, 2009 -С. 136-139.

47. Трофименко C.B., Гриб H.H., Никитин В.М. Особенности сейсмического режима переходной области восточного замыкания Байкальской рифтовой и западного фланга Олекмо - Становой зон // Современные методы обработки и интерпретации сейсмологических данных. Материалы Международной сейсмологической школы. - Обнинск: ГС РАН, 2009. -С.195-199.

48. Трофименко C.B. Термическая модель Южно-Якутского землетрясения // Шестой всероссийский симпозиум «Физика геосфер», 7-11 сентября 2009г., Владивосток.-ТОИ ДВО РАН, 2009.-С.250-255.

49. Трофименко C.B., Гриб H.H., Никитин В.М. Электромагнитный мониторинг переходной зоны сочленения Алданского щита и Становой складчатой системы // Шестой всероссийский симпозиум «Физика геосфер», 7-11 сентября 2009г., Владивосток.- ТОЙ ДВО РАН, 2009.-С.256-259.

50. Трофименко С. В. Динамика сейсмического режима Олекмо-Становой сейсмической зоны // «Тектонофизика и актуальные вопросы наук о земле. К 40-летию создания

M.B. Гзовским лаборатории тектонофизики в ИФЗ РАН». Материалы конференции. М.: Изд-во ИФЗ РАН. 2009. Т.2. С.403-410.

51. Трофименко С. В. Предварительный анализ статистической модели аномалий гравимагнитных полей Алданского щита.- Вестник Технического института (филиала) Якутского государственного университета,- Выпуск 4.-Нерюнгри: изд-во Технического института, 2009,- С.86-96.

52. Трофименко С. В. Тектонофизические особенности проявления сейсмического режима Олекмо - Становой зоны - Вестник Технического институ та (филиала) Якутского государственного университета- Выпуск 4.-Нерюнгри: изд-во Технического института, 2009.- С.97-101.

53. Трофименко С.В., Гриб H.H., Никитин В.М. 22-х летние циклы сейсмической активности сейсмических поясов Северо-востока Азии //Проблемы сейсмичности и современной геодинамики Дальнего Востока и Восточ! :i Сибири. Доклады научного симпозиума 1-4 июня 2010, г. Хабаровск /под.ред.В.Г. Быкова, А.Н. Дидекко.- Хабаровск: ИТИГ им. IO.A. Косыгина ДВО РАН, 2010.-С.291-294.

54. Трофименко С.В. Геофизический мониторинг геологической среды для прогнозирования сейсмической опасности //Проблемы сейсмичности и современной геодинамики Дальнего Востока и Восточной Сибири. Доклады научного симпозиума 1-4 июня 2010, г. Хабаровск /под.ред.В.Г. Быкова, А.Н. Диденко.- Хабаровск: ИТИГ им. Ю.А. Косыгина ДВО РАН, 2010.-С.295-298.

С.В. Трофименко

Автореферат диссертации на соискание ученой степени доктора геолого-минералогических наук

Подписано в печать 23.12.2010. Формат 60x84/16. Бумага тип. №2. Гарнитура «Тайме». Печать офсетная. Печ. п. 2,75. Тираж 120 экз. Заказ 04 . Издательство ТИ (ф) СВФУ, 678960, г. Нерюнгри, ул. Кравченко, 16

Отпечатано в ООО «Печатный дом». I. Нерюнгри, ул. Южно-Якутская, 22.

Содержание диссертации, доктора геолого-минералогических наук, Трофименко, Сергей Владимирович

ВВЕДЕНИЕ

ГЛАВА 1. Модели геофизической среды, вариаций физических полей сейсмотектонической природы и сейсмических процессов

1.1. Модели геофизической среды и сейсмических процессов

1.2. Модели геофизических полей сейсмотектонической природы и 15 проблема прогноза землетрясений

1.3. Сейсмотектоническая позиция изучаемой территории

1.4. Обсуждения и выводы по главе 1 42'

ГЛАВА 2.

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

2.1 Предпосылки для построения статистических пространствен- 46 ных моделей геофизических полей земной коры

2.1.1 Пространственная периодичность распределений геофизи- 46 ческих полей и зон очаговой сейсмической активности

2.1.2 Структурно-геологическая изученность территории иссле- 51 дований

2.1.3. Изучение строения земной коры по аномалиям геофизи- 66 ческих полей

2.2 Методика и технология построения статистических систем ано- 71 малий геофизических полей

2.3 Анализ статистики распределений азимутов линейных индика- 74 торов аномалий физических полей

2.4 Анализ пространственного распределения плотностных неод- 78 нородностей линейных индикаторов аномалий

2.5 Геолого-геофизическая. интерпретация статистической модели 82 пространственной структуры аномалий магнитного и гравитационного полей

2.6 Прогнозирование сейсмоактивных зон на основе статистиче- 87 ской модели геофизических полей

2.7 Пространственное распределение плотностных неоднородно- 89* стей линейных индикаторов аномалий как отражение деформационных процессов

2.9 Обсуждения и выводы по главе

ГЛАВА 3. Статистические пространственно-временные модели сейсмичности

3.1 Изучение динамики распределения энергетических параметров^ сейсмичности Олекмо-Становой сейсмической-зоны

3.1.1. Изменение энергии и количества землетрясений-во времени.

3.1.2. Статистические оценки повторяемости землетрясений в модели стационарного Пуассоновского потока событий

3.2 Оценка стационарности сейсмического процесса сформированной области афтершоков сильного землетрясения М=7.

3.3 Изучение суточных и годовых распределений сейсмических событий

ГЛАВА 4.

ГЛАВА 5.

3.3.1. Распределение количества землетрясений по часам в течение суток

3.3.2. Пространственное распределение эпицентров землетрясений в максимумах суточной активности

Проявление землетрясений и их фор-афтершоков на фоне стационарного сейсмического процесса

Особенности 22-двух летних циклов сейсмической активности Изучение пространственной структуры очаговой и рассеянной сейсмичности Олекмо-Становой сейсмической зоны Геофизические процессы в переходных областях сейсмоактивных зон

Модельные оценки баланса энергии очага и расчет коэффициента полезного действия землетрясения Выводы по главе

Геофизические исследования переходных зон взаимодействия литосферных плит

Детальные геофизические исследования в зонах активных разломов

4.1.1. Исходные данные для производства детальных геофизических работ

4.1.2. Методика электроразведочных работ

4.1.3. Методика сейсморазведочных работ

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

Комплексная интерпретация результатов геофизических исследований зон активных разломов Обсуждения и выводы по главе

Изучение физических полей в связи сейсмичностью в переходных зонах литосферных плит

Дистанционный геофизический мониторинг активных тектонических структур

5.1.1 Вариации геофизических полей в сейсмоактивных зонах 5.1.2. Интерпретация результатов наблюдений вариаций электромагнитного излучения (ЭМИ) в блоковой модели земной коры

Обобщение статистических моделей на сейсмические зоны северо-востока Азии

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

Введение Диссертация по наукам о земле, на тему "Структура и динамика геофизических полей и сейсмических процессов в блоковой модели земной коры"

Диссертационная работа «Структура и динамика геофизических полей и сейсмических процессов в блоковой модели земной коры» посвящена изучению общих, региональных и локальных характеристик геофизических полей как реакцию на сейсмотектонические процессы литосферы.

Актуальность темы.

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

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

Природа сейсмичности рассматривалась во всех основополагающих работах по физике Земли. В различных подходах данные вопросы изложены в ставших теперь классическими монографиях [Джеффрис, 1960; Гуттенберг, 1963; В.А. Магницкий, 1965; Дарвин, 1965; Стейси, 1972; Ботг, 1974; Буллен, 1978; Цубои, 1982; Эйби, 1982; Касахара, 1985; Ризниченко, 1985; Гир, Шах, 1988]. Новые подходы к анализу сейсмичности изложены в монографиях [Садовский, Писаренко, 1991; Садовский, 2004; Соболев, 1993; Быков, 2000; Викулин, 2003; Арефьев, 2003; Долгих, 2004; Ребецкий, 2007; Гуфельд, 2007; Хаин, Хали-лов, 2008].

Существующие представления об источниках аномалий геофизических полей предвестников землетрясений [Гохберг и др., 1988; Аксенов, 1982; Воробьев, 1971, 1980; Электромагнитные предвестники ., 1982; Виноградов, 1989] базируются на физических [Физика очага .,1975; Физические основания .,1970; Шамина, 1981, Соболев, 1980, 1989] моделях очага землетрясения.

Традиционный путь решения прогнозных задач — анализ корреляционных связей между аномальными проявлениями в физических полях и пространственным распределением, механизмами и динамикой очагов землетрясений с привлечением геоморфологических, геологических, тектонических и космических критериев сейсмичности.

За практически 50-ти летний период разработки данного направления исследований наиболее детально изучено влияние катастрофических землетрясений на динамику геофизической среды. Не исследованными остаются вопросы формирования аномальных полей вследствие миграционных волн и временных циклов сейсмической активности.

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

Открытыми остаются вопросы по изучению состояния геологической (деформируемой) среды в зонах активных разломов геофизическими методами в современный (го-лоценовый) период.

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

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

Объектом данного исследования являются физические поля магнитной, гравитационной и электромагнитной природы, продуцируемые переходной зоной взаимодействия Амурской, Евроазиатской и Охотоморской литосферных плит.

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

1 .проведение статистических исследований по распределению азимутов линейных элементов аномалий физических полей гравитационной и магнитной природы с выделением преобладающих азимутов, построение моделей (систем) по распределению плотности индикаторов линейных элементов аномалий и поиск закономерностей их пространственной структуры.

2.установление подобия выделенных систем структурным тектоническим элементам исследуемой территории.

3.построение статистических моделей пространственно - временных распреде

12 лений слабых землетрясений с энергией Е < 10 Дж и изучение их закономерностей.

4.анализ динамических характеристик сейсмического режима сильных землетрясе-12 ний Е > 10 Дж и рассеянной сейсмичности за инструментальный период. Проверка выявленных закономерностей на очагах сильных землетрясений Байкальской рифтовой зоны, Сейсмического пояса Черского, Тихоокеанской островной дуги (п-ов Камчатка, о. Сахалин), обобщение материалов модельных построений и выявление общих закономерностей статистических моделей сейсмического процесса.

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

6.проведение мониторинговых геофизических исследований в зонах активных разломов для оценки их современной тектонической активности и изучение прогностических эффектов с разработкой технологии моделирования геофизических полей и сейсмических процессов переходных зон и построения геодинамической модели сейсмического процесса.

Научные положения, защищаемые в диссертации:

1. Пространственное распределение плотности линейных элементов аномалий гравитационного и магнитного полей (индикаторов) обнаруживает высокую степень упорядоченности в виде ромбических структур (систем доменов) с устойчивыми линейными размерами и является следствием блокового строение геофизической среды.

2.Проявления статистических закономерностей сейсмичности в характерных чувствительных зонах в суточных, годичных и в среднем 22-х летних циклах обусловлены влиянием общих возмущающих факторов вследствие взаимодействия физических полей в системе «Солнце-Земля-Луна» и идентично для быстрых геодинамических процессов в виде афтершоковых последовательностей и медленных стационарных сейсмических процессов.

З.Зона релаксации энергии упругих деформаций вследствие сильных землетрясений охватывает область активных разломов (по сейсмологическому признаку), в виде устойчивой геодинамической системы не зависимо от механизма очага землетрясения, энергии землетрясения, начиная с некоторой пороговой величины и положения очага внутри области. Длительность процессов релаксации пропорциональна энергии землетрясения.

Научная новизна работы:

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

• методами статистического анализа распределений физических полей гравитационной и магнитной природы определен линейный размер неоднородности геофизической среды (домена). Впервые показано, что доменные структуры с линейными размерами от 120 до 150км образуют взаимосогласованные пространственные решетки, тем самым про1 ведено обобщение теоретической модели формирования транспрессионных структур Вуд-кока - Фишера (Woodcock N.H., Fisher М, 1986.); i

• развиты теоретические модели сейсмических процессов в переходных зонах ли-I тосферных плит. Впервые показано что пространственное распределение очагов земле-< трясений формируется в векторном поле внешних сил, а временное распределение землетрясений в течение суток идентично для быстрых геодинамических процессов в виде аф-тершоковых последовательностей и медленных (стационарных) сейсмических процессов, тем самым доказан принцип согласованности самоподобия пространственной структуры земной коры и автомодельности сейсмического процесса;

При содействии заведующего ПНИЛЗ ЯГУ A.C. Стативы автором организованы и проведены деформационные и др. виды геофизических исследований в Олекминском эпи-центральном поле катастрофических землетрясений. В период Южно-Якутского землетрясения 1989г. с магнитудой М-6.6 при исследовании геофизических полей-предвестников землетрясений автором установлены закономерности формирования вариаций геофизических полей вследствие сейсмотектонических процессов и выделены характерные периоды вариаций сопоставимые с периодами сейсмических процессов.

Научное и практическое значение.

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

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

Геолого-геофизические материалы автора по центральной части Алданского щита вошли в отчет института Физики Земли РАН и использованы при проектировании трубопроводной системы ВСТО в части уточнения сейсмической опасности.

Модельные построения по изучению пространственной неоднородности геофизической среды методами статистического анализа распределений линейных элементов аномалий гравитационного и магнитного полей могут иметь практическую значимость для уточнения и прогноза дифференцированных зон сейсмической опасности, т.к. выделенные азимуты и границы ромбических структур (систем доменов) коррелируют с пространственным распределением очагов землетрясений.

Апробация работы и публикации.

Результаты диссертационной работы докладывались на следующих конференциях: пятые геофизические чтения имени В. В. Федынского 27 февраля - 01 марта 2003 года.-Москва, центр ГЕОН,(2003); международной конференции «Проблемы сейсмологии Ш-го тысячелетия», Новосибирск (2003); 2-м Международном симпозиуме «Активный геофизический мониторинг литосферы Земли», Новосибирск (2005); международной сейсмологическая школа «Современные методы обработки и интерпретации сейсмологических данных», Обнинск (2006, 2008, 2009, 2010); IV-VI Всероссийских симпозиумах «Сейсмоа-кустика переходных зон», Владивосток, ТОЙ ДВО РАН (2005, 2007, 2009); Всероссийской научно-практической конференции «Сейсмичность Южно-Якутского региона и прилегающих территорий», Нерюнгри (2005); IV-V Международной конференция «Солнечно -земные связи и предвестники землетрясений», с. Паратунка, Камчатской обл. (2007, 2010); Всероссийском совещании с международным участием «Проблемы современной сейсмологии и геодинамики Центральной и Восточной Азии», Иркутск ИЗК СО РАН (2007); международной научно-практическая конференция «Южная Якутия - новый этап индустриального развития», Нерюнгри (2007); IV Международном симпозиуме «Геодинамика внутриконтинентальных орогенов и геоэкологические проблемы», Бишкек (2008); Всероссийской конференции «Тектонофизика и актуальные вопросы наук о Земле», Москва -ИФЗ РАН - 20-25 октября 2008; Всероссийской конференции «Тектоника и глубинное строение востока Азии» Хабаровск: ИТиГ им. Ю.А. Косыгина ДВО РАН - 20 - 23 января 2009; Всероссийском совещании с участием приглашенных исследователей из других стран «Разломообразование и сейсмичность в литосфере: тектонофизические концепции и следствия», Иркутск ИЗК СО РАН (2009); на второй научно-технической конференции

Проблемы комплексного геофизического мониторинга Дальнего Востока России» 11-17 октября 2009г. г.Петропавловск-Камчатский; на научно-практической конференции «Сейсмические исследования земной коры» (Пузыревские чтения) Новосибирск: Институт нефтегазовой геологии и геофизики им. A.A. Трофимука СО РАН 22-25 ноября 2009; на научном симпозиуме «Кайнозойский континентальный рифтогенез» Иркутск ИЗК СО РАН (2010); на научном симпозиуме «Проблемы сейсмичности и современной геодинамики Дальнего Востока и Восточной Сибири» Хабаровск: ИТиГ им. Ю.А. Косыгина ДВО РАН - 1-4 июня 2010.

Заключение Диссертация по теме "Геофизика, геофизические методы поисков полезных ископаемых", Трофименко, Сергей Владимирович

ЗАКЛЮЧЕНИЕ

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

В работе приведены результаты моделирования геофизических полей и сейсмических процессов, базирующихся на экспериментальном материале, полученном при исследовании процессов подготовки землетрясений при помощи геофизического мониторинга геологической среды. Основная доля экспериментальных исследований выполнена на геодинамических полигонах Южной Якутии, в непосредственной близости от Олёкмин-ского эпицентрального поля катастрофических землетрясений с магнитудой М=6.5 - 7. В результате проведенной интерпретации суточных изменений электромагнитного излучения в период Южно-Якутского землетрясения обоснована модель аномалий в блоковой структуре литосферы. Данный результат стал стимулом дальнейших исследований по структурированию геофизических полей и сейсмических процессов.

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

35-305° и 40-310°) простираний линейных элементов аномалий гравимагнитных полей может свидетельствовать о том, что изменения палеотектонических обстановок и соответствующих ориентировок векторов деформаций на изучаемой территории происходили не единовременно и резко, а постепенно и плавно, что может найти объяснение в неравномерном (циклическом) вращении, или "косой субдукции" Амурской плиты.

5. На основании анализа распределений поля сейсмичности всей ОСЗ в периоды сильных землетрясений получены аналитические зависимости для длительности действия афтершоков.

19 шается до оа < 10 Дж/год. Длительность периода повышенной сейсмической опасности

I 9 и{,< 10 Дж/год) равна или превышает длительность периода повышенной сейсмической

19 активности (ис1 > 10 Дж/год). Сейсмический процесс Олекмо - Становой сейсмической зоны развивается под действием региональных напряжений плейттектонической природы, а видимая скорость диссипации энергии иа в отдельных частях ОСЗ определяется упруго пластическими свойствами горных пород.

8. Сопоставлены структурные элементы систем разломов с пространственным полем сейсмичности. На основании данного исследования установлено их геометрического подобие, как на региональном уровне для всей ОСЗ, так и на локальном уровне по полю афтершоков сильных землетрясений.

10. Выявленные закономерности сейсмического процесса ОСЗ подтверждены результатами моделирования сейсмического режима Байкальской рифтовой зоны, сейсмического пояса Черского, отдельных областей Тихоокеанской островной дуги. Для всех сильных землетрясений определена длительность афтершоковых последовательностей, укладывающаяся в единое феноменологическое уравнение; для узлов скрещивающихся активных разломов установлена последовательная активизация в течение суток (в статистическом смысле), для сдвиговых деформаций выявлены ромбические структуры, характерные для присдвиговых бассейнов (pull-apart basin); для всех исследуемых зон выделен 22-х летний период миграции сейсмичности с максимумом широты в 1983г., определено пороговое значение энергии землетрясений, меняющее структуру потока сейсмических событий, величина которого находится в пределах 1012'5- 1013 Дж; наличие динамических брешей за

1 ^ тишья выявлено для всех землетрясений с энергией более 10 Дж, при этом размеры области затишья не зависят от энергии землетрясений и определяют динамическую систему, объединенную общим деформационным процессом плейттектонической природы.

Библиография Диссертация по наукам о земле, доктора геолого-минералогических наук, Трофименко, Сергей Владимирович, Нерюнгри

1. Аксенов В.В. Интерпретация электромагнитных вариаций. М.: Наука, 1982.264с.

2. Алексеев A.B., Долгих Г.И. Сейсмоакустика переходных зон //Четвертый Всероссийский симпозиум «Сейсмоакустика переходных зон». Материалы докладов. Владивосток: Изд-во Дальневост. Ун-та, 2005.- с.3-4.

3. Арефьев С.С. Эпицентральные сейсмологические исследования. М.: ИКЦ «Академкнига», 2003. 376 с.

4. Балакина JI.M., Введенская A.B., Голубева Н.В., Мишарина J1.A., Широкова Е.И. /Поле упругих напряжений Земли и механизм очагов землетрясений. М.: Наука. - 1972. -191 с.

5. Билибина Т.В. К геологии и металлогении Алданского щита //Вопросы генезиса и закономерности размещения эндогенных месторождений. М.: Наука, 1966. - С. 233-256.

6. Билибина Т.В., Крутихевская З.А., Титов В.К. Региональные магнитные аномалии и тектоника Алданского щита //Советская геология. 1986. - №1. - С. 106-115.

7. Билибина Т.В., Дашкова А.Д., Донаков В.И. Петрология щелочного вулканогенно-интрузивного комплекса Алданского щита. Л.:Недра, 1967. - 263 с.

8. Ботт М. Внутреннее строение Земли. М.: Мир, 1974. 376 с.

9. Брейс В. и др. Две модели объяснения предвестников землетрясений.- В кн.: Сборник советско-американских работ по прогнозу землетрясений.-Душанбе,1976.-Т.1.-к.2.-С.9-21.

10. Бриджмен П. Исследование больших пластических деформаций и разрыва.-М.: Изд-во иностранной литературы.- 1955.-444с.

11. Бровар B.B. Влияние изменений атмосферного давления на силу тяжести //Повторные гравиметрические наблюдения.- М.: Из-во ВНИИгеофизики,1980.-С.27-33.

12. Быков В. Г. Нелинейные волновые процессы в геологических средах. Владивосток: Дальнаука, 2000. 190 с.

13. Быков В. Г. Деформационные волны Земли: концепция, наблюдения и модели // Геология и геофизика. 2005. Т. 46. № 11. С. 1176-1190.

14. Буланже Ю.Д. Неприливные изменения силы тяжести//Повторные гравиметрические наблюдения.-М.:Из-во ВНИИгеофизики.-1983.

15. Буланже Ю. Д. Некоторые результаты изучения неприливных изменений силы тяжести. В кн. Проблемы расширения и пульсации Земли. М.: Наука, 1984, с.73-84.

16. Буллен К.Е. Плотность Земли. М.: Мир, 1978. 442 с.

17. Викулин A.B.www.kscnet.ru/monograhp/ivs/vikulin/index.html

18. Викулин A.B., Журавлев В.И. Статистический анализ распределения Курило-Камчатских землетрясений во времени и в пространстве. Петропавловск-Камчатский: ИВ ДВНЦ АН СССР. 1987. Препринт. 28 с.

19. Викулин A.B. Особенности распределения землетрясений северо-западной части Тихого океана / Препринт. Петропавловск-Камчатский: ИВ ДВО АН СССР. 1987. 54 с.

20. Викулин A.B. Феноменологическая волновая модель сейсмического процесса // Докл АН СССР. 1990. Т.310. № 4. С.821-824.

21. Викулин A.B. Прогноз времени сильнейших землетрясений у берегов Камчатки и Северных Курил // Вулканология и сейсмология. 1992. № 1. С.62-69.

22. Викулин А. В. Физика волнового сейсмического процесса. Петропавловск-Камчатский: КОМСП ГС РАН, 2003. 150 с.

23. Викулин А. В. Энергия и момент силы упругого ротационного поля геофизической среды // Геология и геофизика. 2008. Т. 49. № 6. С. 559-570.

24. Витинский Ю. И. Цикличность и прогнозы солнечной активности., JL: Наука, 1973,458 с.

25. Воробьев A.A. Физические условия залегания вещества в земных недрах. Томск: Изд-во Томского политех, ин-та, 1971. - Ч. 1. - 270 с.

26. Воробьев A.A. Равновесие и преобразование видов энергии в недрах. Томск: Изд-во Томского ун-та, 1980. - 211 с.

27. Гамбурцев Г.А., В.В.Белоусов Прогноз землетрясений. Г.А. Гамбурцев. Избранные труды. М.: Изд. АН СССР. 1960. С. 436-443.

28. Гамбурцев А.Г. Сейсмический мониторинг. М.: Наука, 1992. 200 с.

29. Гамбурцев Г.А. Воспоминания, очерки, статьи. М.: ОИФЗ РАН. 1998. 270с.

30. Геодинамика Олекмо-Становой сейсмической зоны / Л.М.Парфенов, Б.М.Козьмин, В.С.Имаев и др. Якутск: ЯФ СО АН СССР, 1985. 136 с.

31. Геология и сейсмичность зоны БАМ. Новосибирск. Сибирское отделение: Наука, 1983.- 188 с.

32. Геологическая карта СССР масштаба 1:200000. Литы: 0-51-XXIX, XXXV (авторы И. А. Петрова, О. А. Янополец). Госкомгеология РС(Я), 2004.

33. Геологическая карта Южной Якутии масштаба 1:500 000 /Под ред. Л.М.Реутова. -Якутск, 1976.

34. Гершензон Н.И. и др. Об источниках электромагнитного излучения, предваряющего сейсмическое событие // Прогноз землетрясений.-Душанбе:Дониш,1986.-№7.-С.54-62.

35. Гир ДЖ., Шах X. Зыбкая твердь. М., Мир, 1988. 220 с.

36. Гольдин C.B., Тимофеев В.Ю., М. ван Раумбеке, Ардюков Д.Г., Лаврентьев М.Е., Се-дусов Р.Г. Приливная модуляция слабой сейсмичности для южной части Сибири //Физическая мезомеханика,- 2008.-№11(4).-С. 81-93.

37. Горьковатый H.H., Левицкий Л.С., Тайдакова Т.А. и др. О корреляции графиков угловой скорости вращения Земли и модуля ее временной производной с чистотой землетрясений в зависимости от магнитуды / Физика Земли. 1994, № 10. - С. 33-38.

38. Горьковатый H.H., Трапезников Ю.А., Фридман A.M. О глобальной составляющей сейсмического процесса и её связи с наблюденными особенностями вращения Земли // Докл. РАН. 1994а. - Т. 338, №4. - С. 525 - 527.

39. Гохберг М.Б. и др. Процессы подготовки, признаки и предвестники коровых землетрясений.- Известия АН СССР, Физика Земли, 1983.-№2.-С.59-67.

40. Гохберг М.Б. О связи аномальных ЭМИ с деформациями при подготовке землетрясений // Прогноз землетрясений.-Душанбе.:Дониш,1986.-№7.-С.288.

41. Гохберг М.Б., Моргунов В.А., Матвеев И.В. О наблюдении аномального электромагнитного излучения в эпицентральной зоне землетрясений.-Изв.АН СССР, Физика Земли, 1986.-№8.-С.95-97.

42. Гохберг М.Б., Моргунов В.А., Похотелов O.A. Сейсмоэлектромагнитные явления.-М.:Наука.-1988.-174с.

43. Грачев А.Ф., Деменицкая P.M., Карасик A.M. Проблемы связи Момского континентального рифта со структурой срединно-океанического хребта Гаккеля, Геофизические методы разведки в Арктике, Л., НИИГА, 1973, вып. 8, С. 56-75.

44. Гриднев Д.Г. и др. Региональные приливные исследования в Азиатской части СССР.- В кн.: Современная геодинамика литосферы Сибири ( Сб.научн.тр. ).- Новосибирск: ИГиГ, 1986.-С.63-79.

45. Гришкян Р.И. О механизме образования позднедокембрийских разрывных нарушений в кристаллических толщах центральной части Алданского щита //Геотектоника. -1968. №4. - С. 136-139.

46. Гришкян Р.И. Новые данные о строении и развитии системы докембрийских разломов решетчатого типа на Алданском щите //Вопросы тектоники докембрия континентов. М.: Наука, 1970. - С. 120-131.

47. Гутенберг Б. Физика земных недр. М.: ИЛ, 1963. 264 с.

48. Гусев A.A. Прогноз землетрясений по статистике сейсмичности // Сейсмичность и сейсмический прогноз, свойства верхней мантии и их связь с вулканизмом на Камчатке. Новосибирск: Наука, 1974. С.109-119.

49. Гусев A.A. Многомасштабное порядковое группирование в последовательности землетрясений Земли //Физика Земли, 2005, № 10, с. 30-45.

50. Гусева Е.В., Добровольский И.П. Вертикальные движения и магнитные вариации перед землетрясением.-Физика Земли,1984.-№1.-С.49-53.

51. Гуфельд И.Л. Сейсмический процесс. Физико-химические аспекты. Королев: ЦНИИМАШ. 2007.160 с.

52. Данилкин Н. Прогноз землетрясений из космоса http://www.novosti-kosmonavtiki.ru/content/numbers/238/26.shtml.- электронный ресурс

53. Дарвин Д.Г. Приливы и родственные им явления в солнечной системе. М.: Мир, 1965.207 с.

54. Джалилов Н. С., Ораевский В. Н. К теории 11-летнего цикла солнечной активности. Известия Академии Наук. Серия физическая. 2000, т.64, №9, с. 1793 1798.

55. Джеффрис Г. Земля ее происхождение история и строение. М.: Изд-во иностранной литературы,- i960.- 486с.

56. Джурик В. И., Дреннов А. Д., Басов А. Д. Прогноз сейсмических воздействий в условиях криолитозоны. Новосибирск: Издательство СО РАН, 2000. 272 с.

57. Дзевановский Ю.К., Миронюк E.JL, Лагздина Г.Ю. История геологического развития Алданского щита //Тектоника Сибири, т.З. -М.: Наука, 1970. С. 132-142.

58. Дзевановский Ю.К. Мезозойские гранитоиды хр. Станового и их структурное положение //Инф. сб. ВСЕГЕИ. 1959. - № 7. - С. 55-60.

59. Дзевановский Ю.К. Архейский метаморфический комплекс Алданского щита // Тр. Межд. совещ. по стратиграфии Сибири. Стратиграфия докембрия.- М.- 1958.- С. 37-42.

60. Добровольский И.П. Механика подготовки тектонического землетрясения.-М.: ИФЗ АН СССР,1984.-189с.

61. Добровольский И.П. Некоторые вопросы изучения процесса подготовки землетрясения.-В кн.: Прогноз землетрясений.-Душанбе: Дониш,1986.-№7.-С.7-15.

62. Долгих Г.И. Исследование волновых полей океана и литосферы лазерно-интерференционными методами. Владивосток: Дальнаука. 2004 г. 160 с.

63. Емельяненко A.C. Структурные границы Становой зоны тектоно-магматической активизации /Тектоника активизированных областей. -Чита, 1982. С. 135-137.

64. Зубков С. И. Предвестники землетрясений. М.: ОИФЗ РАН.2002. 140 с.

65. Зубков С.К. О зависимости времени возникновения и радиуса зоны проявления.-Физика Земли,1983.-С.101-106.

66. Имаев B.C. Тектонические критерии сейсмичности Южной Якутии. М.: Наука, 1986.- 127 с.

67. Имаев B.C., Имаева Л.П., Козьмин Б.М. Активные разломы и сейсмотектоника Северо-Восточной Якутии, Якутск, ЯНЦ СО АН СССР, 1990, 138 с.

68. Имаев В. С., Имаева Л. П., Козьмин Б. М. Сейсмотектоника Якутии. М.: ГЕОС.-2000,- 226с.

69. Имаев B.C., Трофименко C.B., Гриб H.H. и др. Разломная тектоника и геодинамика в моделях очаговых зон сильных землетрясений Южной Якутии. Томск: Изд-во Томского политехнического университета.- 2007. 274с.

70. Казанский В.И., Иванов И.Б. Два этапа метаморфизма в раз ломах Центрально-Алданского района//Докл. АН СССР. 1968. - т. 181. № 1. - С. 169-170.

71. Казанский В.И. О внутреннем строении архейских разломовЦентрально-Алданского района//Геология рудных месторождений. 1965. - T. VII, № 2. - С. 63-79.

72. Карасев П.С., Овсюченко А.Н., Мараханов A.B., Трофименко C.B. Активные нек-тонические нарушения участка Алдан-Нагорный нефтепроводной системы Восточная Сибирь — Тихий океан.- Нефтяное хозяйство, 2008. № 9, с. 80-84.

73. Калинин Ю.Д. Солнечная обусловленность изменения длины суток и сейсмической активности. Красноярск: Инст. Физ. СО АН ССР, 1974, 23 с.

74. Касахара К. Механика землетрясений. М.: Мир, 1985. - 260с.

75. Кейлис-Борок В.И. Динамика литосферы и прогноз землетрясений // Природа, 1989. N2. С.10-18.

76. Кейлис-Борок В.И. Повторный сильный толчок землетрясений: прогноз возможен. Наука в России. 1992. №1. С.60-63.

77. Ключевский A.B., Демьянович В.М. Напряженно-деформированное состояние литосферы северо-восточного фланга Байкальского региона по данным о сейсмических моментах землетрясений // Вулканология и сейсмология. 2006. № 2. С.65-78.

78. Козьмин Б.М. Сильные землетрясения Якутии (макросейсмические и инструментальные данные за 1851-1974 гг.) // Сейсмическое районирование Якутии и сопредельных территорий. — Якутск, 1975. С. 5-19.

79. Козьмин Б.М. Механизм очагов сильных южно-якутских землетрясений. // В кн.: Сейсмические и сейсмогеологические исследования на центральном участке БАМа. — Якутск: ЯФ СО АН СССР, 1978. С. 58-71.

80. Козьмин Б.М. Сейсмические пояса Якутии и механизмы очагов их землетрясений. -М.: Наука, 1984.-125 с.

81. Козьмин Б.М., Голенецкий С.И., Николаев В.В. и др. Южно-Якутское землетрясение 20 апреля 1989 года,- Южно-Сахалинск: ИМГиГ ДВО РАН, 1992,- 45 с

82. Козьмин Б.М., Голенецкий С.И., Николаев В.В. и др. Афтершоки Южно-Якутского землетрясения 20 апреля 1989 года.- Южно-Сахалинск: ИМГиГ ДВО РАН, 1993.- 66 с.

83. Кокс Д., П. Льюис Статистический анализ последовательности событий. М.: Мир, 1969.-312с.

84. Красный Л.И. Геологическая картина мира и проблемы тектонической систематики //Труды ВСЕГЕИ. 1972. - 188 с.

85. Кропоткин П. Н. Трапезников Ю. А. Вариации угловой скорости вращения Земли, колебаний полюса и скорости дрейфа геомагнитного поля и их возможная связь с геотектоническими процессами. Изв. АН СССР, Сер. Геол., №14, 1963, с.32-50.

86. Кропоткин П.Н. Возможная роль космических факторов в геотектонике// Геотектоника, 1970. N2. С.30-76.

87. Кузнецов В. В. Физика земных катастрофических явлений. 1992. Новосибирск Наука. 96 с.

88. Кузнецов В.В. Физика горячей Земли. Новосибирск. 2000. 365 с.

89. Кузнецов В. В. Введение в физику горячей Земли. Петропавловск-Камчатский: Изд-во КамГУ им. Витуса Беринга, 2008. 367 с.

90. Курскеев A.K. и др. Вариации гравитационного поля в зонах глубинных разломов Алма-Атинского полигона // Повторные гравиметрические наблюдения.-М.:Из-во ВНИИ-геофизики.-1981.

91. Левин Б.В., Родкин М.В., Сасорова Е.В. О возможной природе сейсмической границы на глубине 70км. Докл РАН. 2007. Т. 414, №1. С.101-104.

92. Летников Ф.А., Дорогокупец П.И. К вопросу о роли суперглубинных флюидных систем земного ядра в эндогенных геологических процессах //ДАН,2001. Т. 378, №4. С535-537.

93. Логачев H.A., Борняков С.А., Шерман С.И. О механизме формирования Байкальской рифтовой зоны по результатам физического моделирования. // Докл. РАН. 2000. Т. 373, № 3. С.388-390.

94. Любушин A.A., Писаренко В.Ф., Ружич В.В., Буддо В.Ю. Выделение периодично-стей в сейсмическом режиме //Вулканология и сейсмология. 1998. №1. - С. 62 — 76.

95. Магницкий В. А. Внутреннее строение и физика Земли. М.: Недра, 1965. 380 с.

96. Максимов Е.П., Угрюмов А.Н. Геологическое строение центральной части Алданского щита //Геология и петрология докембрия Алданского щита, М.: Наука, 1966. - С. 51-38.

97. Максимов Е.П. Опыт Формационного анализа мезозойских магматических образований Алданского щита //Известия АН СССР. Сер. геол. 1975. - №4. - С. 16-32.

98. Малышев Ю.Ф. Геофизические исследования докембрия Алданского щита. М.: Наука, 1977, 128 с.

99. Минкин Л.М, Закономерности размещения и возраст мезозойских интрузий Алдана //Материалы по геологии и полезным ископаемым Якутской АССР. Вып.2. - М.: Госгео-лтехиздат, 1960. - С. 3-22.

100. Моги К. Предсказание землетрясений (перевод с англ.). М.: Мир. 1988. 382 с.

101. Мокроусов В.А. Тектоника нижнего докембрия Алданского щита //Тектоника Сибири, т.З. М.: Наука, 1970. - С. 150-156.

102. Мокшанцев К.Б., Горнштейн Д.К., Гусев Г.С., Деньгин Э.В., Штех Г.И. Тектоническое строение Якутской ССР. М.: Наука, 1964. 240 с.

103. Моргунов В.А. Реальности прогноза землетрясений. Физика Земли. 1999. №1. С.7991.

104. Мухин В.А. Зоны повышенной тектонической активности в западной части Алданского щита и связь с ними хрусталеносных жил /Тр.Всес. научно-исслед. ин-та пьезокварца, т.7, 1962.

105. Мячкин В.И. и др. Лабораторные и теоретические исследования процесса подготовки землетрясений.-Изв.АН СССР,Физика Земли, 1974.-№10.-С. 107-112.

106. Николаев В.В., Семенов Р.Ф., Семенова В.Г., Солоненко В.П. Сейсмотектоника, вулканы и сейсмическое районирование хребта Станового. Новосибирск: Наука, 1982. 150 с.

107. Николаев A.B. Реакция сильных землетрясений на фазы земных приливов / Физика Земли. 1996. - № 11. - С. 49-58.

108. Овсюченко А.Н., Трофименко C.B., Мараханов A.B., П.С., Карасев, Е.А. Рогожин Очаговые зоны сильных землетрясений Южной Якутии.- Физика Земли. 2009. №2. С. 1533.

109. Овсюченко А.Н., Трофименко C.B., Мараханов A.B., П.С., Карасев, Е.А. Рогожин,

110. B.C. Имаев, В. М. Никитин, Н.Н.Гриб Детальные геолого-геофизические исследования зон активных разломов и сейсмическая опасность Южно-Якутского региона.- Тихоокеанская геология, 2009, том 28, №4, с. 55-74.

111. Овсюченко А. Н., Трофименко С. В., А. В. Мараханов, П. С. Карасев, Е. А. Рогожин Сейсмотектоника переходной области от Байкальской рифтовой зоны к орогенному поднятию Станового хребта. Геотектоника, 2010, № 1, с. 29-51.

112. Павловский Е.В. О специфике стиля тектонического развития земной коры в раннем докембрии /Труды Вост.-Сиб. геол. ин-та. 1952. - Вып.5. М.: Изд-во АН СССР, 1962. -С. 77-108.

113. Парфенов Л.М., Натальин Б.А., Попеко Л.И. Тектоника восточной части Монголо-Охотской складчатой системы и природа сопряженных с ней зон мезозойской тектоно-магматической активизации Алдано-Становой области. Тихоокеанская геол., 1983, № 2,1. C. 26-33.

114. Парфенов Л.М., Козьмин Б.М., Имаев В.С и др. Геодинамика Олекмо-Становой сейсмической зоны. Якутск: ЯФ СО АН СССР, 1985. - 136 с.

115. Парфенов Л.М., Козьмин Б.М., Гриненко О.В. и др., Сейсмичность и геодинамика Восточной Сибири, Современная тектоническая активность Земли и проблемы сейсмичности, М, 1986, С 26-27.

116. Пейве A.B. Тектоника и магматизм // Изв. АН СССР, сер. геологическая. 1961. № 3. С. 36-54.

117. Попов Н.В., Смелов А.П. Метаморфические формации Алданского щита. Геология и геофизика, 1996, т. 37, № 1, с. 148-161.

118. Пригожин И., Стенгерс И. Порядок из хаоса. М.: Прогресс, 1989. 241 с.

119. Разломы и горизонтальные движения платформенных областей СССР. М.: Наука, 1977. - 143 с.

120. Разломная тектоника территории Якутской АССР. Якутск, ЯФСО АНСССР, 1976,173 с.

121. Рао С.Р. Линейные статистические методы и их применение. — М.: Наука. — 1968. —548с.

122. Ребецкий Ю.Л. Тектонические напряжения и прочность природных горных массивов. Научное издание. М.: ИКЦ «Академкнига», 2007. - 406с.

123. Ризниченко Ю.В. Метод суммирования землетрясений для изучения сейсмической активности // Изв. АН СССР. сер. геофиз. 1964. № 7. С.969-977.

124. Ризниченко Ю.В. Размеры очага корового землетрясения и сейсмический момент // Исследования по физике землетрясений. М.: Наука. 1976. С.9-27.

125. Ризниченко Ю.В. Проблемы сейсмологии. Избранные труды. М.: Наука. 1985. 408с.

126. Садовский М. А. О естественной кусковатости горных пород. Доклады АН СССР, 1979, т. 247, № 4, с. 829.

127. Садовский М. А., Волховитинов Л. Г., Писаренко В. Ф. О свойствах дискретности горных пород.-Препринт ИФЗ АН СССР, 1981, № 1.

128. Садовский М.А., Волховитинов Л.Г., Писаренко В.Ф. Деформирование геофизической среды и сейсмический процесс. М.: Наука, 1987. 100 с.

129. Садовский М. А. Избранные труды. Геофизика и физика взрыва. М.: Наука, 2004.440 с.

130. Садовский М.А., Писаренко В. Ф. Сейсмический процесс в блоковой среде,- М.: Наука.- 1991.-96с.

131. Салман А.Г., Шилин Б.В. Сейсмическая активность: взгляд из космоса. Природа. 1989. №12. С. 55-58.

132. Салтыков В.А., Иванов В.В., Кугаенко Ю.А. Воздействие земных приливов на сейсмичность перед землетрясением 13 ноября 1993 года Mw=7.0 (Камчатка) // Физика Земли. 2004. №7. С.25-34.

133. Салтыков В.А. О воздействии земных приливов на сейсмические процессы //Проблемы сейсмичности Дальнего Востока / под ред. Викулина A.B. Петропавловск-Камчатский, 2000. С. 12-21.

134. Саньков В.А., Леви К.Г., Лухнев A.B., Мирошниченко А.И. Современные движения литосферных блоков Центральной Азии по данным GPS-геодезии // Актуальные вопросы современной геодинамики Центральной Азии. Новосибирск: Изд-во СО РАН, 2005. С.165-179.

135. Саркисов Ю.М., Ермаков Б.В. Положение разломов Байкало-Амурского региона в общей регматической сети Востока СССР //Разломыи эндогенное оруденение Байкало-Амурского региона. М.: Наука,1982. - С. 15-21.

136. Сергин С.Я., Сергин В.Я. Взаимодействие литосферы и климата как одна из причин возникновения геологических циклов. Бюлл. МОИП. Отд. Геол. Т.62, вып.2, 1987, с. 3-17.

137. Сергин С.Я., Сергин В.Я. Природа глобальных геологических циклов. Системный подход. М., Наука, 1993 123 с.

138. Сидоренков Н.С. Межгодовые колебания системы атмосфера-океан-Земля. Физика, № 25.-1998.

139. Сидоренков Н.С. http://www. 1 september.ru/ru/fiz/2003/01 /поО 1 1 .htm электронныйресурс

140. Сидорин А.Я. Зависимость времени проявления предвестника землетрясений от эпицентрального расстояния,- ДАН СССР,1979.-Т.245.-№4.-С.225-228.

141. Сидорин А.Я. Суточная периодичность землетрясений и ее сезонные изменения //Сейсмические приборы. 2009. Т. 45, № 4, с. 69-84

142. Сидорин А.Г. Предвестники землетрясений. М.: Наука, 1992. 191 с.

143. Сковородкин Ю.П. Изучение тектонических процессов методами магнитометрии.-Из-во ИФЗ АН СССР.-1985.-197с.

144. Соболев Г.А., Демин В.М. Механоэлектрические явления в Земле. М.: Наука, 1980.-215 с.

145. Соболев Г.А. Проблема прогноза землетрясений. Природа. 1989. №12. С.47-55.

146. Соболев Г.А. Основы прогноза землетрясений. М.: Наука, 1993. 312 с.

147. Соболев Г.А. Стадии подготовки сильных камчатских землетрясений // Вулканология и сейсмология. 1999. -№4—5. С. 63-72.

148. Соломатина Э.К. Предвидеть земную бурю. Наука в СССР. 1990. №3. С. 5-13.

149. Солоненко В.П. Определение эпицентральных зон землетрясений по геологическим признакам // Изв. АН СССР. Сер. геол. 1962. № 11. С. 58-74.

150. Солоненко В.П. Сейсмогеология, геофизические поля Монголо-Охотского сейсмического пояса и проблемы прогноза землетрясений / Основные проблемы сейсмотектоники. М.: Наука. - 1986. - С. 171-177.

151. Солоненко A.B., Штейман Е.А. 1994. Самоподобие поля сейсмичности Байкальского рифта // Доклады РАН. Т. 337, № 2. С. 253-257.

152. Статива A.C., Трофименко C.B. Методика и технология выделения систем разломов Алданского щита // Горный информационно-аналитический бюллетень/Региональное приложение ЯКУТИЯ. Выпуск 3, 2006. Издательство МГГУ. С. 155-166.

153. Статива A.C., Трофименко C.B. Построение систем разломов Алданского щита по геофизическим данным //Горный информационно-аналитический бюллетень/Региональное приложение ЯКУТИЯ. Выпуск 3, 2006. Издательство МГГУ. С. 193196.

154. Статива А. С., Трофименко С. В. Ортогональные системы разломов Алданского щита и их геологическая позиция //Горный информационно-аналитический бюллетень/Региональное приложение ЯКУТИЯ. Выпуск 4, 2005. Издательство МГГУ.

155. Стаховский И.Р. Временной и пространственно-временной мультимасштабный анализ сейсмичности в период подготовки Рачинского землетрясения // Физика Земли. 2000. №4. С. 41-47.

156. Стейси Ф. Физика Земли. М.: Мир, 1972. 344 с.

157. Стогний В.В., Стогний Г.А. Тектоническая раслоенность Алдано-Станового геоблока. Новосибирск: Наука, 1997. 151 с.

158. Стогний Г.А., Стогний В.В. Геофизические поля восточной части СевероАзиатского кратона. Якутск: ГУП НИПК «Сахаполиграфиздат», 2005. - 174 с.

159. Суворов В.Д., Корнилова З.А. Глубинное строение Алданского щита по данным сейсмологии близких землетрясений // Геол. и Геоф.1985.- N 2.- С. 86-89.

160. Сытинский А.Д. О влиянии солнечной активности на сейсмичность Земли. ДАН СССР. - 1973. -Т. 209, № 15. - С. 1078-1081.

161. Таймазов Д. Г. О возможности прогноза землетрясений по скачкообразным изменениям деформационного поля. Тез. докл. межд. симпоз. "Мониторинг и прогнозирование чрезвычайных ситуаций". Махачкала, 1997. С. 25-26.

162. Тектоника, геодинамика и металлогения территории Республики Саха (Якутия) / Отв. ред. J1.M. Парфенов, М.И. Кузьмин. М.: МАИК "Наука/Интерпериодика", 2001. -571 с.

163. Тимофеев Д.А. Некоторые общие вопросы неотектоники и морфотектоники на примере геоморфологической истории Южной Якутии // Проблемы геоморфологии и неотектоники орогенных областей Сибири и Дальнего Востока.- Новосибирск.- 1968.- Т. 2.-С. 235-244.

164. Трофименко C.B. Высокоточные гравиметрические наблюдения при условии низких температур // Геофизические исследования в Якутии,- Якутск: ЯГУ, 1990.-С.54-60.

165. Трофименко C.B., Широков И.А., Анохина К. М.Результаты земноприливных наблюдений на территории Южной Якутии // Геофизические исследования в Якутии,-Якутск: ЯГУ, 1995.-С.55-63.

166. Гриб H.H., Трофименко C.B. Сейсмичность Олекмо-Становой сейсмической зоны инициированная ядерными взрывами // Горный информационно-аналитический бюллетень/Региональное приложение ЯКУТИЯ. Выпуск 4, 2005. Издательство МГГУ. С. 167175.

167. Трофименко C.B. Проявление землетрясений на фоне стационарного сейсмического процесса Олекмо-Становой зоны (ОСЗ) //Горный информационно-аналитический бюллетень. Отдельный выпуск 3. Якутия 1. М.: Изд-во МГГУ, 2007. - С. 208-212.

168. Тяпкин К.Ф., Гонтаренко В.Н. Системы разломов Украинского щита. Киев: Наук, думка, 1990. - 184 с.

169. Уломов В.И. Волны сейсмогеодинамической активизации и долгосрочный прогноз землетрясений // Физика Земли. 1993. № 4. С. 43-53.

170. Уломов В.И., Полякова Т.П., Медведева Н.С. О долгосрочном прогнозе сильных землетрясений в Центральной Азии и в Черноморско-Каспийском регионе // Физика Земли. №4. 2002. С. 31-47.

171. Уломов В.И. От сейсмического районирования к прогнозу землетрясений // Сб. статей "Научное наследие Г.А. Гамбурцева и современная геофизика: Развитие идей. Воспоминания". М.: ОИФЗ РАН, 2003. С. 189-202.

172. Уломов В.И. Сейсмичность // Национальный атлас России. Том 2. Природа. Экология. 2007. С. 56-57.

173. Федотов С.А. О закономерностях распределения сильных землетрясений Камчатки, Курильских островов и северо-восточной Японии // Труды ИФЗ АН СССР. 1965. № 203 (36). С.66-93.

174. Федотов С.А. О сейсмическом цикле, возможности количественного сейсмического районирования и долгосрочном сейсмическом прогнозе // Сейсмическое районирование СССР. М.: Наука. 1968. С.121-150.

175. Федотов С.А. О сейсмичности области очага катастрофического Итурупского землетрясения 6.11.1958 г. и сейсмическом прогнозе // Изв. АН СССР. Физ. Земли. 1969. № 1. С.3-12.

176. Федотов С.А., Чернышев С.Д., Чернышева Г.В., Викулин A.B. Уточнение границ очагов землетрясений с M > ТА, свойств сейсмического цикла и долгосрочного сейсмического прогноза для Курило-Камчатской дуги //Вулканология и сейсмология. 1980. .№ 6. С.52-67.

177. Федотов С.А., Чернышев С.Д. 20 лет долгосрочного сейсмического прогноза для Курило-Камчатской дуги: достоверность в 1981-1985 гг., в целом за 1965-1985 гг. и прогноз на 1986-1990 гг. // Вулканология и сейсмология. 1987. № 6. С.93-109.

178. Фирстов П.П. Мониторинг объемной активности подпочвенного радона на Пара-тунской геотермальной системе в 1997-1998 гг. с целью поиска предвестников сильных землетрясений Камчатки // Вулканология и сейсмология, 1999. N 6. С.22-32.

179. Фотиади Э.Э., Есиков Н.П., Бочаров Г.В. и др. Тектонофизическое районирование зоны сочленения Алданского щита и Станового нагорья по данным геодезических измерений // Геология и геофизика. 1982. №3. С.3-15.

180. Фотиади Э.Э., Шарловская JI.A., Захарова Т.Л. Строение земной коры и верхней мантии региона БАМ. Новосибирск: Наука, 1988. - 103 с.

181. Фролов А. Д. Электрические и упругие свойства мерзлых пород и льдов. М.: ОНТИ ПНЦ РАН, 1998. 515 с.

182. Фрумкин И.М. Стратиграфия и тектоника архея восточней части Алданского щита. Якутск: ЯГУ, 1971.-272 с.

183. Хаин В.Е., Короновский Н.В. Планета Земля от ядра до ионосферы. М.: КДУ, 2007.244 с.

184. Хаин В.Е., Халилов Э.Н. Пространственно-временные закономерности сейсмической и вулканической активности. Бургас, SWB, 2008, 304 с.

185. Халилов Э.Н. Выявление некоторых особенностей вулканической активности по оценкам функций спектральных плотностей. Вулканология и сейсмология. М.: № 3, 1987, с. 104-107.

186. Хирт Дж., Ломе И. Теория дислокаций. М.: Атомиздат, 1972. - 600с.

187. Хотина Е.Б., Жежель О.Н. Кайнозой бассейна Алдана // Кайнозойский седименто-генез и структурная геоморфология СССР. Л.: Недра, 1987. С.29-33.

188. Хренов П.М. Скрытые зоны глубинных разломов складчатых областей и платформ (на примерз юга Восточной Сибири) //Глубинные разломы юга Восточней Сибири и их металлогеническое значение. -М.: Наука, 1971. С. 5-38.

189. Худяков Г.И. Принципиальные основы морфотектонических исследований (на примере изучения южной материковой части Советского Дальнего Востока).: Автореф. Дис. докт. геол.-мин. наук. Новосибирск, 1974, 34 с

190. Цубои Т. Гравитационное поле Земли.- М.: Мир, 1982. 286с.

191. Черкасов Р.Ф. Архей Алданского щита. М.: Наука, 1979. - 160 с.

192. Шамина О.Г. Модельные исследования физики очага землетрясения. М.: Наука, 1981.-192с.

193. Шерман С.И., Днепровский Ю.И. Новая карта полей напряжений Байкальской рифтовой зоны по геолого-структурным данным // Докл. РАН. 1986. Т. 287, № 4. С.943-947.

194. Шерман С.И., Леви К.Г. Трансформные разломы Байкальской рифтовой зоны и сейсмичность ее флангов // Тектоника и сейсмичность континентальных рифтовых зон. -М.: Наука, 1978.- С. 7-18.

195. Шерман С.И., Семинский К.Ж., Черемных A.B. Разломно блоковая тектоника Центральной Азии: опыт тектонофизического анализа //Актуальные вопросы современной геодинамики Центральной Азии. - Новосибирск: Изд-во СО РАН, 2005. - 297 с.

196. Широков В.А. Влияние 19-летнего лунного прилива на возникновение больших камчатских извержений и землетрясений и их долгосрочный прогноз // Геологические и геофизические данные о БТТИ 1975-1976 гг. М.: Наука, 1978. С. 164-170.

197. Широков В. А. Некоторые вопросы методики комплексного прогноза побочных извержений вулкана Ключевского // Вулканология и сейсмология. N6. 1985. С.40-55.

198. Широков В.А., Кузьмин Ю.Д. О суточной цикличности вулканических землетрясений в Курило-Камчатской зоне // Вулканология и сейсмология. 1988. N5. С.92-100.

199. Шулейкин В.В. Взаимодействие звеньев в системе «Океан-Атмосфера-Материки».- Природа, № 10.-1971.

200. Эйби Дж. А. Землетрясения. М.: Недра, 1982. - 264с.

201. Юшманов В.В. Типизация мезозайнозойских тектоно-магматических комплексов Алдано-Станового региона /Тектоника активизированных областей. Чита, 1982. - С. 7275.

202. Янополец О.А., Петрова И.А., Скопич С.И. и др. Государственная геологическая карта РФ масштаба 1:200 ООО. Серия Алданская. Листы 0-51-XXIX, 0-51-XXXY: Объясн. зап. / Филиал ГУП PC (Я) "Алдангеология" "Южноякутская экспедиция". СПб., 2004.

203. Aki К. Earthquake generating stress in Japan for the years 1961 to 1963 obtained by smoothing the first motion radiation patterns. Bull. Earthquakes Res. Inst., Univ. Tokyo. 44, №2, 1966.

204. Anderson E.M. The dynamics of faulting and dyke formations with applications to Britain.- Edinburgh.- London.- Oliver and Boyd.- 1951.- 206 p.

205. Ashida M. Faultu premise. Sciences. 1996. Vol. 36, №5. P. 15-19.

206. Byerlee J.D. Friction of Rocks //Pure and applied geophysics. 1978. V 116. - P. 615626.

207. Doser D.I. Faulting within the eastern Baikal rift as characterized by earthquake studies // Tectonophysics. 1991. Vol. 196. P.109-139.

208. FujitaN. The magnetic distucbance accompaning the Nigata earthquake,J.Geod.-v. 11.-p.8.-1965y.

209. Fujita K., Stone D.B., Layer P.W., Parfenov L.M. and Koz'min B.M., Cooperative program helps decipher tectonics of Northeastern Russia, Eos, v.78, 1997, p.245, 252-253.

210. Hamaguchi H., Hasegava A. An investigation on the aftershocks of the Tokachi-oki earthquake of 1968. 2. Statistical study on time distribution // Sci. Repts Tohoku Univ. 1970 Sep. 5. 20. N3. P.l 19-133.

211. Kato V., Utashiro S. On the changes of the terres trial magnetic field accompanyng the GreatNankaido earthquake of 1946, Sci, Rep.Tohoku Univ.Ser.5, №l.-p.40.-1948y.

212. Kelleher J., Sykes L., Oliver J. Possible criteria for predicting earthquake locations and their application to major plate boundaries of the Pacific and the Caribbean // J. Geophys. Res. 1973. N 14. P.2547-2585.

213. Mackey K.G., Fujita K., Gunbina L.V., Kovalev V.N., Imaev V.S., Koz'min B.M. and Imaeva L.P., Seismicity of the Bering Strait region: evidance for a Bering block: Geology, 1997, 64(6), p. 467-471.

214. Meyer K.B., Olsson R., Kulhanek O. High-velocity migration of large earthquakes along the Azores-Iran plate boundary segment // Pure and Appl. Geophys. 1985. 122. N 6. P.831-847.

215. Mogi K. Migration of seismic activity // Bull, of the Eathquake Res. Inst. 1968 a. V.46. P.53-74.

216. Mogi K. Development of aftershock areas of great earthquakes // Bull, of the Earthquake Res. Inst. 1968 b. V.46. P.l75-203.

217. Mogi K. Two kinds of seismic gaps // Pure and Appl. Geophys. 1979. V.117. P.l 1721186.

218. Nagata T. Tectonomagnetism, Jnt, Assoc. Geomag. Aeron. Bull., №27, p. 12, 1969y.

219. Nagata T. Application of tectonomagnetism to earthquake phenomena, Tectonophis., №14.- p.263.-1972y.

220. Otsuka M. Подобие последовательностей землетрясений. 4.8. Описание афтершо-ковой активности по закону Омори // J. Seismol. Soc. Jap. 1987. 40. N 1. Р.69-75 (яп.).

221. Page R. Aftershocks and micro aftershocks of the Great Alaska earthquake of 1964 // Bull. Seismol. Soc. Amer. 1968. 58. N 3. P.l 131-1168.

222. Reid HF.The California earthquake of April 18 1906. V.2 // The mechanics of the earthquake. The Carnegie Inst. Washington, 1910.

223. Reid HF. The elastic-rebound theory of earthquakes. University of California Publ. Geol. Sci. 1911. V. 6. P. 413-444.

224. Savostin L.A., Karasik A.M., Recent plate tectonics of the Arctic basin and Northeastern Asia, Tectonophysics, 1981, vol. 74, P.l 11-145.

225. Sykes L. Aftershock zones of great earthquakes, seismicity gaps, and earthquake prediction for Alaska and the Aleutians // J. Geophys. Res. 1971. 76. N 32. P.8021-8041.

226. Thanassoulas, C., Klentos, V., 2003. Seismic potential map of Greece, calculated by the application of the "Lithospheric energy flow model". Open File Report A.4403, Inst. Of Geol. Min. Expl. IGME, Greece.

227. Tazima M. Accurasy of resent magnetic survey and a locally anomalies behaviour of the geomagnetic secular variation in Japan // Bull Geograph,Earthq.Res.Jnat.-v.l3.-p.21.- 1968y.

228. Tazima M.H.,Mizuna and Tanaka M. Geomagnetic secular change anomaly in Japan // J.Geomag.and Geoelectr.-v.28.-p.69.-1976y.

229. Vamazaki V., Rikitake T. Local anomalies changes in the geomagnetic field at Matsu-shira // Bull.Earthq.Res.Jnst.-v.48.-p.637.-1970y.

230. Woodcock N.H., Fisher M. Strike slip duplexes //Journ. Of Struct. Geology. 1986. Vol. 8, N7. P. 725-735.1. Фондовая литература

231. Штех Г.И., Трофименко С.В. Прогноз землетрясений.- Отчет.- Якутск, ЦНТИ, фонд ЯГУ,1982. 95с.

232. Штех Г.И., Трофименко С.В. и др. Провести исследования по выявлению предвестников землетрясений в зоне Южно Якутского территориально промышленного комплекса.- Отчет.-Якутск, ЦНТИ, 1983. 178с.

233. Штех Г.И., Трофименко С.В. и др. Режимные исследования гравитационного, магнитного и теллурического полей в зоне Становика.- Отчет.- Якутск, ЦНТИ, 1985. 444с.

234. Штех Г.И., Трофименко С.В. и др. Создать комплексные прогностические полигоны и станции по поиску новых и изучению известных прогностических эффектов и выявлению наиболее информативных из них.- Отчет.- Якутск, ЦНТИ, 1985. 372с.

235. Статива А.С., Трофименко С.В. и др. Провести исследования по получению и накоплению геофизической информации на прогностических полигонах в зоне Центрального и Южно — Якутского ТПК.- Отчет.- Якутск, ЦНТИ, 1987. 128с.

236. Статива А.С., Трофименко С.В. и др. Изучение Вариаций геофизических полей и деформаций земной поверхности на территории Южной Якутии,- Отчет.- Якутск, ЦНТИ, 1990. 122с.

237. Статива А.С., Трофименко С.В. и др. Проведение режимных исследований геофизических полей а Олекминской эпицентральной зоне,- Отчет.- Якутск, ЦНТИ, 1991. 218с.

238. Статива А.С., Трофименко С.В. и др. Результаты и задачи исследований магнито-теллурического поля на Алданском щите.- Отчет.- Якутск, ЦНТИ, 1991. 140с.

Информация о работе

Трофименко, Сергей Владимирович
доктора геолого-минералогических наук
Нерюнгри, 2011
ВАК 25.00.10

Диссертация

Структура и динамика геофизических полей и сейсмических процессов в блоковой модели земной коры - тема диссертации по наукам о земле, скачайте бесплатно

Автореферат

Похожие работы