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

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

ВАК РФ 25.00.16, Горнопромышленная и нефтегазопромысловая геология, геофизика, маркшейдерское дело и геометрия недр

Содержание диссертации, доктора технических наук, Беклемишев, Андрей Борисович

Глава 1. Общая характеристика работы

Глава 2. Критический обзор известных подходов, история вопроса, основные термины и определения, концепция исследований

Состав и устойчивость микросейсм, попытки поиска причин микросейсмического движения

Методы оценки СЭ

Выводы

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

Глава 4 Теоретический аспект. «Безсобытийные» модели и механизмы динамики неоднородных сред. Эффективная параметризация, мультидисциплинарная верификация оценок и их интегральных распределений

Понятие о передаточных функциях, способах их оценки и свойствах

Традиционный анализ А ЧХ(со/к)и ФЧХ(со/к)

Дополнения к традиционному анализу

Движение поверхности среды и концепция обратной связи

Упрощенный анализ устойчивости и переходных процессов в системе

Определения и формализмы пространственного движения

Параметры поляризации микросейсм

Энергетические соотношения

Инструментальная индикация дисперсионного типа среды

Эффективная параметризация

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

Глава 5. Требования к методике и технике комплексных полевых исследований, включая скрининг и мониторинг, их реализация, в том числе в новых модификациях трехкомпонентного приемника (ДВЗ-024) и цифровой измерительной системы

Трехкомпонентные приемники естественной акустической эмиссии 89 Развитие и современное состояние полевых систем для регистрации, векторных оценок и измерений параметров акустической эмиссии 95 Мультидисциплинарная инструментальная составляющая полевых измерений: прецизионные грави- и гелиеметрия

Выводы

Глава 6. Фактическая коллекция собранных данных, ее статистическая представительность, инженерная и геологическая содержательность

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

Анализ точности измерений

Многопараметровые оценки

Выводы

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

Инструментальные оценки разлома Topesti-Prisaca (TP fault)

Результаты по локализации и активности корового (меле- и внутриплитного) разлома Naruja-Negrilesti-Soveja (NNS fault), включая его верхнюю унаследованную часть в олигоцен-миоценовой формации; их анализ и интерпретация

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

Изменения релаксируемой энергии микросейсм как оценок, контролирующих степень деградации конструкции и/или геосреды

Выводы

Глава 9. Ретроспективный анализ структурной эволюции геологических объектов по результатам инструментальной фиксации их текущего состояния и среднесрочный прогноз будущей эволюции: качественной трансформации и скорости

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

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

Выводы

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

Введение Диссертация по наукам о земле, на тему "Научно-методологические и экспериментальные основы технологии изучения пассивных сейсмических явлений"

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

ФРАНК СЛ.,

1997)

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

Актуальность проблемы. Управление государственным фондом недр, реализуемое через изучение глубинного строения и эндогенной активности территории России, разработки рациональной стратегии недропользования и оценки связанных с использованием недр и освоением территорий геологических и экономических рисков, -первый и важнейший аспект деятельности МПР РФ (Положение об МПР РФ, утвержденное Постановлением Правительства РФ № 726 от 25.09.2000 г.). Следовательно, необходима современная и надежная база знаний, корректных представлений о геологической среде, ее естественной и/или вынужденной эволюции.

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

• причины локально ярких возмущений,

• доминирующие механизмы естественных движений (упругие и/или вязкие),

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

• структуру и состояние среды, степень их устойчивости,

• направление и скорость структурной эволюции (процессов текущей деструкции или ре-консолидации). Этот пункт, по крайней мере, на сейсмоопасных территориях, может быть жизненно важным показателем устойчивого развития. И очевидные альтернативные условия последнего: либо плавная, монотонная и легко предсказуемая эволюция некоторых ключевых параметров среды (например, структурной прочности, усталости), либо заблаговременное, но сегодня проблематичное, слабо формализованное и отчасти спорное предсказание места, времени, силы и последствий экстремального события (Соболев Г.А., 1987; Соболев Г.А., 1997).

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

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

2) используемые в большинстве случаев приемники с максимальной чувствительностью к вертикальным движениям адаптированы исключительно к линейно поляризованным одноосным колебаниям,

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

4) оценки интенсивности (кинетической энергии), как правило, сложно поляризованного волнового движения поверхности, ненадежны и часто в различной степени ошибочно занижены, t< )

5) мощные зондирующие импульсы «активной» сейсморазведки способны изменять исходную структуру геосреды и, особенно, ВЧР (по крайней мере, вблизи источника), ее проницаемость и гидрогеологический режим, а также прочность встроенных инженерных сооружений и конструкций, то есть необратимо модифицирует среду как объект исследования (Адушкин В.В. и Спивак А.А., 1993),

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

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

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

• трехкомпонентные (ЗС) высокочувствительные (на 2-3 порядка превосходящие индустриальный стандарт) измерения микросейсмического ускорения идентичными широкополосными (от 0,1 до 300 Гц) каналами с полным динамическим диапазоном близким к 100 дб,

• продолжительные (до 60С), статистически представительные, синхронизированные временные ряды одноточечных измерений в рамках однократного скрининга территории и/или прерывистого мониторинга (многократный скрининг в идентичных климатических, сезонных условиях).

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

- интенсивность микросейсмического вектора в заданном временном интервале, диапазон и характер текущих изменений этого параметра;

- интервалы внутри временных рядов с преимущественно линейно поляризованным или с преимущественно плоско поляризованным движением, чтобы различить случаи доминирования объемных волн глубинного генезиса или поверхностных волн;

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

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

При развитии отечественной версии технологии изучения пассивных сейсмических явлений удалось отказаться от широко применяемой «событийной концепции», когда оценки вырабатываются в связи с реакцией на определенное заметное событие: землетрясение или взрыв. Этот отказ носил сначала гипотетический характер и основывался только на обнаруженном непрерывном характере фоновых микросейсм. Однако, теперь весь многолетний опыт пассивной сейсморазведки, коллекции ЗС данных, собранные на геологически разнообразных объектах (различные и обширные территории Русской платформы, Воронежского кристаллического массива и Балтийского щита, ЮВ Карпат и ЮВ Предкарпатья), эффективная параметризация и результаты анализа поддерживают нижеследующий тезис.

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

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

- через определение потенциально наиболее эффективных направлений и параметров для приложений ЗС микросейсмических измерений;

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

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

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

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

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

• методики одноточечного оперативного сбора ЗС данных микросейсмических измерений (скрининг), обновления этих данных (прерывистый мониторинг) и увязки с данными высокоточной грави-, наклоно-, гелиеметрии и других методов исследования в различной геологической обстановке;

• количественные подходы к оценке текущего состояния и динамики погребенныз зон деструкции и разломов (особенно «не залеченных», активных), а также к геодинамическому прогнозу для территорий, включающих встроенные важные объекты инфраструктуры (АЭС, магистральные трубопроводы, хранилища радиоактивных отходов, мосты и т.п.), базирующемуся на пространственно-временных распределениях эффективных спектрально-корреляционных параметров, производных из векторных измерений.

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

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

Ранее (Gutenberg В., 1936; Prescott W.H., 1982; Sammis C.G. et al., 1987) экспериментально показано, что резкое увеличение ущерба, наносимого землетрясением, приурочено к близким, импедансно контрастным и часто скрытым латеральным зонам неоднородности (разломы, трещины, карстовые полости, каверны - деструктивные зоны; крутые кромки аллювиальных долин - «топографические» зоны). Максимальный ущерб локализован в сравнительно узкой полосе, в границах относительно мягкой акустически области. Положение и характеристики источника, траектории распространения возмущения имеют при этом подчиненное значение. Определяющие резонанс контактные условия (импедансный контраст > 3) играют главную роль. Разумеется, внутри таких полос успешность проектов индустриальной сейсморазведки и высокоточной гравиметрии является, по меньшей мере, спорной. Сейсмические риски здесь высоки, и, при необходимости практического использования территории, требуется изучить их пространственное распределение детально и количественно (Соболев Г.Л., 1997). Предложенные и обоснованные в работе карты добротности среды (в версии нововведенного FI-параметра как оценки динамики правых спектральных срезов для полного микросейсмического вектора) поэтому вполне эффективны и рациональны для идентификации и оконтуривания названных небезопасных (с высоким уровнем добротности) полос - резонаторов.

Контроль границ областей высоких резонансных свойств и добротности среды - не самоцель, поскольку им дополняется основа для эффективного планирования кондиционных геофизических исследований на конкретных территориях. Как результат такого дополнения растет надежность комплексных прогнозов и заключений относительно строения и эволюции геологической формации. Характер и скорость эволюции оказались разными: для древних платформ и бассейнов отмечена только медленная деструкция; для относительно молодых и активных формаций - чередование быстрых деструктивных процессов, связанных с мощными землетрясениями, и относительно медленных (но заметных в пятилетних повторных наблюдениях) конструктивных процессов «залечивания» или ре-консолидации. В новой эффективной технологии предусмотрена также количественная оценка и сопоставление других динамических параметров. Эти параметры - интенсивность, скорость изменения интенсивности, спектральное отношение горизонтальной и вертикальной составляющих вектора (HVSR-параметр), частотная привязка и уровни шести последовательно убывающих доминант спектра, отношение коэффициентов поляризации (Kp/Ki-параметр, где Кр и К| - коэффициенты плоской и линейной поляризации соответственно), скорость изменения параметра Kp/Ki, и степень корреляционной связи флуктуирующих оценок (АКФ и ВКФ соответственно).

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

Математическая модель для составляющих обычно слабого колебательного процесса вдоль каждой из трех ортогональных осей, представлена (на примере оси X) в форме линейного неоднородного дифференциального уравнения второго порядка с постоянными коэффициентами (упругости - к и динамической вязкости - h): d2X2(t)/dt2+2hdX2(t)/dt+k2X2(t)=k2X1(t) (Bl), где Xlj2(t) - временные ряды возмущения и реакции соответственно (в смещениях). Вещественные решения (В1) в частотной области (как передаточные функции среды - ПФС): nOC(co/k)=ModSX2(t)/ModSXi (t)={ [ 1 -(co/k)2]2+4 (co/k)2(h/k)2}I/2 (В2), где Mod SXi,2(t) - модули комплексных спектров для временных рядов Xij2(t) соответственно, со - круговая частота, h/k = относительный коэффициент демпфирования, экспериментально подтверждены как эффективные для описания (покомпонентно и средневзвешенно) поля реальных «безсобытийных» микросейсм. Модель допускает разные значения коэффициентов по координатным осям и, таким образом, подходит неоднородным и анизотропным, вязким и упругим средам.

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

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

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

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

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

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

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

1. выявление геодинамически наиболее потенциально опасных по максимальному градиенту параметра (СЭ) участков активных разломных зон, в пределах которых недопустимо размещение особо важных объектов инфраструктуры;

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

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

Фактический материал, методы исследования и личный вклад автора. Измерения, сохраненные в виде временных рядов оценок, составляют проанализированную базу данных, относящуюся исключительно к «безсобытийным» микросейсмам. ЗС измерения выполнены с использованием серийных и опытных образцов акселерометров от четырех различных производителей и адаптированного ко всем этим акселерометрам опытного образца цифровой полевой регистрирующей системы «Вибротестер-МС». Функцией системы является, в частности, предобработка зарегистрированных временных рядов. Интервал выборки — 2 мС. Несколько тысяч продолжительных микросейсмических выборок скрининга собрано и сохраняется, начиная с 1994 г., в виде электронного архива, доступного для просмотра и углубленной обработки с помощью вновь разработанных и стандартных программных средств. Эти выборки относятся к различным и обширным регионам древних формаций Русской платформы, Воронежского и Балтийского щитов, а также относительно молодых, активных, продолжающих формироваться образований Карпатского орогена и Предкарпатского краевого прогиба. При этом детальность наблюдений варьирует от 1 м до нескольких км. Географическая привязка измерений выполнена на основе опубликованных карт, планов и чертежей встроенных сооружений, а также на основе GPS данных фиксированной точности. Повторные микросейсмические измерения (прерывистый мониторинг, TLM) выполнены с интервалами продолжительностью от нескольких месяцев до 5 лет в сходных климатических и сезонных условиях, в дневное время. Известно, что достоверность однометодных оценок принципиально невелика для многофакторных природных явлений и моделей, поэтому ЗС микросейсмические данные в большинстве случаев сопровождены и заверены прецизионными и синхронизированными данными высокоточной грави-, наклоно- и гелиеметрии. Использована также доступная скважинная и другая геологическая информация. В результате реализована мультидисциплинарная верификация содержательности и представительности исходных микросейсмических данных, динамики производных параметров и итоговой интерпретации.

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

Предлагаемые результаты исследований, составляющие основу для инновационной технологии, были получены, начиная с 90-х годов по настоящее время, исключительно диссертантом или при его непосредственном участии в итоге плановой и/или инициативной научно-производственной и организационной деятельности в ГФУП ВНИИГеофизика МПР РФ и в ИФЗ РАН.

Апробация работы. Результаты, включенные в диссертацию, доложены в течение последних десяти лет на следующих всероссийских и международных конгрессах, конференциях, симпозиумах, чтениях: Новые методики и технологии геофизических исследований на основе комплексирования методов и взаимодействия полей, Геофиз. конф. ЕАГО и ВНИИГеофизика, Москва, 2005; 4th Congress of the Balkan Geophysical Society,

Bucharest, 9-12 October 2005; Workshop on preparation of large scale fault zone map of Tehran area and guidelines for civil construction, Tehran, Iran, 2004; Междунар. геофиз. конф, Москва, 2003; Междунар. конф. "Научное наследие академика Г. А. Гамбурцева и современная геофизика", Москва, 2003; Междунар. конф. "Живая тектоника", Москва, 2003; Четвертые геофиз. чтения им. В.В. Федынского, Москва, 2002; Конф. и выставка по подпрограмме "Топливо и энергетика" Программы "Научные исследования Высшей школы по приоритетным направлениям науки и техники", Москва, 2001; Междунар. геофиз. конф. "300 лет горно - геологической службе России", Санкт-Петербург, 2000; XXI11 General Assembly of EGS, Nice, 1998; Third research symposium "Developing and managing turbidite reservoirs. Case histories and experiences", Almeria, Spain, 1998; Cairo' 98 Africa/Middle East Second International Geophysical Conference and Exp., 1998; XXI1 General Assembly of EGS, Vienna, 1997; 59 EAGE Conference, Geneva, 1997; Istanbul' 97 International Geophysical Conference and Exp., 1997; XXI General Assembly of EGS, Hague, 1996; 30 International Geological Congress, Beijing, China, 1996; Междунар. геофиз. конф.и выставка, Санкт-Петербург, 1995.

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

Публикации.

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

1. Динамика реального («безеобытийного») микросейсмического вектора наилучшим образом определяется в спектре микросейсм через средневзвешенные распределения эмпирических оценок трех текущих однокомпонентных (по осям X,Y,Z) относительных коэффициентов демпфирования - эффективных, статистически устойчивых и надежных параметров идентификации ведущего механизма движения (в подавляющем большинстве случаев имеющего упруго-вязкую, в более редких случаях - вязкую и еще более редко - реликтовую квазиупругую природу), что позволяет I производить инструментальное районирование среды по типу преимущественного механизма движения.

2. Более высокая крутизна правого («высокочастотного») из срезов зарегистрированного микросейсмического спектра соответствует случаю вязкоупругой или квазиупругой среды с нормальной дисперсией. Более высокая- крутизна левого («низкочастотного») из срезов зарегистрированного микросейсмического спектра соответствует случаю вязкой (часто водонасыщенной) среды с аномальной дисперсией. Независимая оценка дисперсионного параметра увеличивает надежность идентификации механизма движения.

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

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

Исследования проведены в соответствии с планами НИР по Федеральным целевым программам и проектам фундаментальных и прикладных исследований МПР РФ и РАН: «Экология и природные ресурсы России, 2002-2010 гг.»; «Сейсмобезопасность территории России, 2002-2010 гг.»; «Развитие федеральной системы сейсмологических наблюдений и прогноза землетрясений, 1995-2000 гг.». С середины 90-х годов инициативные проекты НИР и экспедиции поддержаны грантами Российского фонда фундаментальных исследований

РФФИ): № 98-05-65087, 1998-2000 гг. и др., приборные НИР - грантом Минобразования и науки РФ.

Объем и структура работы. Диссертация включает в себя 13 разделов: введение (1), девять глав основной части (2-10), в том числе три главы (2-4), характеризующие состояние проблемы, заключение, выводы и рекомендации (11), список литературы из 146 наименований (12) и сведения об авторе (13). Полный объем диссертации 233 страниц: 176-текст; 63 - иллюстрации (59) и таблицы (7).

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

Выводы

Представлены полученные диссертантом с коллегами оригинальные и свежие данные о динамике полного вектора микросейсмических движений вблизи двух коровых разломов: в ЮВ Карпатах (TPV) и на Мезийской плите (IMF), а также некоторые результаты их анализа в рамках первичного скрининга и двух фаз прерывистого мониторинга (TLM), разделенных пятилетним промежутком (lapse). Конечно, еще до постановки измерений была вполне ясной множественность причин, вызывающих микросейсмы. A priori ясной являлась также множественность доминирующих волновых типов и распределенных во времени и пространстве возможных траекторий распространения микросейсмической энергии в сложной среде.

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

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

Время от времени, особенно после крупных событий с катастрофическими последствиями, подобных цунами 26.12.2004 в ЮВ Азии, активизируется широкий интерес к проблемам геодинамики, предвестникам событий и прогнозированию последних. Этот интерес, разумеется, оправдан, но, чтобы быть адекватным проблемам, он очевидно должен быть непрерывным или, по крайней мере, периодическим. Учет эволюции принятых геодинамических ЗС микросейсмических параметров, реализуемый в пространственных координатах при скрининге и в пространственно-временных координатах при TLM, представляется в этом смысле содержательным и перспективным. Многочисленные эмпирические данные, производные параметры и результаты их интерпретации свидетельствуют о чередовании эпох (ритмах) низкой и высокой микросейсмической активности. Трактовка этих эпох через механизмы накопления и релаксации (сброса) избыточной энергии в среде - тривиальна. Однако, несомненно, что само по себе признание доминирования тех или других механизмов и даже факта их чередования без надежной оценки характера (скорости) переходного процесса в конкретной геологической формации, учитывающей ее предысторию, мало продуктивно в отношении даже вероятностного прогнозирования даты, силы и места события (Beklemishev А.В. et al., 2004; 2005).

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

Инструментальная идентификация действующих в геологической среде и на ее поверхности фактических механизмов нагружения - релаксации, количественная их характеристика, предусмотрены только в отечественных, совместимых с TLM, модификациях пассивной сейсморазведки. Как результат, эти легкие и экономичные модификации приобретают способность заблаговременно (до реализации проекта) прогнозировать качество данных индустриальной сейсморазведки и высокоточной гравиметрии, используя прежде всего карты добротности (FI). На этапе проектирования эти карты полезны для оптимизации пространственной привязки позиций сейсмоприемников и профилей к упругим высокодобротным участкам, а опорных гравиметрических сетей - к вязким. Безусловным эффектом такого прогноза качества являются улучшенная увязка (сопоставимость) мультидисциплинарных данных и большая надежность комплексной интерпретации. Следствиями являются: уточнение геометрии и режима погребенных каналов энерго- и массопереноса, флюидо- и газонасыщения этих каналов и степени современной активности глубинных разломов и вообще деструктивных зон; коррекция прогнозов локальных сейсмических рисков и рекомендаций инженерного характера.

На рис. 7.4 было представлено состояние на момент публикации (Беклемишев А.Б. и др., 2004) текущей деструкции осадочного чехла в значениях FI параметра во внутренней и СЗ периферической частях Московской синеклизы (Тверская область). Правый сегмент этого распределенияя, не более 10%, соответствующий высоким значениям FI, характеризует не измененную (со времени заложения) или слабо измененную консолидированную геологическую формацию. Левый сегмент распределения, около 20%, соответствующий предельно низким значениям FI, характеризует полностью деструктурированную, деконсолидированную, фрагментированную формацию в погребенных разломных зонах. Наиболее представительный центральный сегмент распределения, около 70%, соответствующий промежуточным значениям FI, имеет ясно выраженный максимум, более 12%, который указывает наиболее вероятное состояние текущей'деструкции формации на дату измерений. Очевидна асимметрия центрального сегмента относительно максимума: более вероятна формация, пока слабо затронутая деструктивными процессами. В итоге имеет место следующая интерпретация: деструктивная геологическая история древнего бассейна является в значительной мере незавершенной, то есть не продвинувшейся повсеместно до минимальных значений FI параметра. При этом реликтовый статус, под которым понимается эмпирическое предельно высокое значение FI, практически повсюду утрачен, но скорость деструктивной эволюции существенно различна для разных фрагментов бассейна. Действительно, эта скорость максимальна для фрагментов бассейна, соответствующих левой периферии распределения и убывает по мере смещения направо. В правом предельном случае эта скорость близка к нулю.

На рис. 10.1 (TLM) и 10.2 (скрининг) представлены карты, отображающие зоны деструкции и близреликтового статуса осадочного чехла бассейна в диапазонах значений FI параметра, соотнесенных с вышеназванными сегментами распределения (Беклемишев А.Б. и др., 2004).

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

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

По крайней мере, три участка (рис. 10.1) вблизи Каблуково, Юрьево-Девичье и, особенно, Харитоново являются рискованными как потенциальные носители СЭ. Здесь кинетическая энергия микросейсм экстремально высока во всех трех TLM фазах. Характерно, что СЭ прогнозированы в позициях, где переход от предельно вязкой к предельно упругой среде

РИС. 10.1 Зоны завершенной и текущей деструкции, а также близреликтового статуса осадочной формации бассейна, представленные в форме карты в диапазонах значений FI параметра. Результат трехфазного TLM [15]. ы о

Зоны деструкции осадочного чехла по данным ЗС мнкросейсмнческнх исследований, 2003 г.

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

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

На рис. 10.3 представлен гидрогеологический разрез от Ижорского плато на ЮВ до Финского залива вдоль 80 км участка воздымающегося в СЗ направлении Балтийского щита (Ar-Pri-г). Здесь в 20 км прибрежной полосе в течение последнего десятилетия детально изучены ЗС микросейсмические движения и производные параметры с тем, чтобы в Предглинтовой низменности оконтурить основные погребенные каналы массопереноса, связанные с разгрузкой водоносных горизонтов. Самый верхний из них, четвертичный (Q), разгружается через трещины в нижележащие водоносные горизонты в карбонатах ордовика (O1-2) и в породах кембрия, а также через родники у кромки Ижорского плато. Грунтовые воды из этой естественной и сложной напорной системы разгружаются непосредственно в залив и, частично, в нижележащий Ломоносовский водоносный горизонт. Таким образом, область разгрузки содержит сложно организованную в погребенном пространстве систему проводящих каналов и перемычек - водоупоров Лонтовской свиты, Котлинского горизонта и покрышки Гдовского комплекса.

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

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

Водонасыщенные породы канала (в предельном случае - пульпы) по своим акустическим характеристикам (крутизне срезов АЧХ, дисперсионным свойствам и др.) отличаются от

Я2£>

У слоёные c<fopia чек г1 st и - УЗодскосныг-с гсри^оиггг. чегк£ер?7-1Ысу>чых стложенигс tS)- Медоносный zopujCKm opgcSukckusc сгпложенкгс

- Ясуочоскый icpujcj-vn. верхнею -нижнего Кембрия

- ЛоитоёасЫя с$ита €±Ы н.ижне-10 кежория - fogoyrtop fcfffl- Ломоносовский £одсксснъо1 горизонт -Cjlm

- Хотминскый zepujenm позднего npomepejest £R3kt-Sc^ovnop

- Шрекрыбающие Sccjcyncpu гдо£сксго #соохссноьо комплекса:

Тдс^ский. So^CHCCHbcw горизонт

- Стрелбнинский. бс^оноскый горизонт. zgcSckoic комплекса.

- Ярхей -протерозойский бодонхеный zopujou-m 'AR-P2, 2

- Урс&нь педз&пнызе cpgcS-икскочо водоносного горизонта.

PR^cl j--* 1~ Урссеш? подъемных: £сд лсиконссоёского ёсаснссксго -гсри?сн-т.сс

I ? . — Ыстпочник

РЯС. 10.3 Гидрогеологический разрез: Ижорское плато - Финский залив (см. текст). вмещающих пород подобно тому, как водонасыщенные породы внутриразломной зоны отличаются от консолидированной формации вне разлома.

Действительно, на рис. 10.4 а (Беклемишев А.Б., 2003) представлена детальная карта значений FIw параметра по 109 позициям скрининга. Интенсивной окраской здесь выделены перемычки - водоупоры, обладающие относительно высокой упругостью и, соответственно, значениями FIw параметра в пределах ,05-,2. Основные каналы подземного массопереноса в западном направлении (к урезу залива), оставленные незакрашенными, характеризуются относительно высокими вязкостями и относительными коэффициентами демпфирования и, наоборот, низкими значениями FIw параметра. Слабо окрашены полупроницаемые области с эмпирически оцененными промежуточными значениями FIw параметра.

На карте отчетливо видно, что погребенные каналы в северной части территории промышленного предприятия, периметр которого отмечен пунктиром, в отличие от каналов в южной части всюду блокированы или имеют весьма широкие (до 300 и более метров) полупроницаемые барьеры. В южной части территории, наоборот, имеют место относительно тонкие (до 100 м по латерали) полупроницаемые барьеры, а вблизи точки 14 идентифицирован единственный канал (русло) массопереноса шириной около 50 метров.

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

Район работ (экспедиционный проект поддержан РФФИ, 09.2005 г.) представляет собой прямоугольник 30x40 км2, примыкающий к Воронежскому водохранилищу с юга, и включает 66 позиций ЗС измерений. Модуль «безсобытийного» микросейсмического вектора характеризуется здесь следующим сравнительно узким интервалом энергий: -5,6 < lg Е < -4,97. Эмпирические оценки параметров занимают интервалы: ,03 < FI < ,18 и 1,3 < HVSR < 4. Скважиной 7801 (вблизи позиции 1-19) с забоем на глубине 500 м вскрыты кровля коренных кристаллических пород на глубине 75 м, представленная эффузивами преимущественно основного состава, и покровные песчано-глинистые осадочные отложения девона, мела, палеогена, неогена и четвертичных с включениями карбонатных слоев (до двух) мощностью до 10 м. Ближайшие к поверхности грунтовые воды обнаружены на глубине 3-12 м. До

Координатная сетка со стороной ячейки - SO и.; О контур здания; надаеииый трубопровод; ж/д ветка: шоссе; аномально диспергирующая среда; аномально низкие (до 7 Гц.} доминирующие компоненты АЧС ЗС микроускорений;

Параметр деструктуриэации среды (Fl-lg dt/iav) в пределах от -0.35 (30) до +0.20 (77): в диапазоне от +0,05 до +0,20 в диапазоне от-0.08 до +0.05 Г 97? - .83" ;

I 1,'Bfi V I ч7 \

•101.ой I ч.б2 .во -J—г». 1—'"Ч ,■. у"96 , * '«.56 * 64 . , V -г 72

Ш I 81 Ц

К"81 Л» \

ЩЗЩ: «Ы V ла Ч? Л

ШЗ "81 ь9 bJ* fi5 asa

•. - : •• • шшттщ штштШШ&ш г.

РИС. 10.4 а Основные каналы подземного массопереноса, полупроницаемые перемычки - барьеры и водоупоры поверх кристаллических (Ar-Pri-2) пород Балтийского щита, представленные в форме детальной карты в диапазонах значений FIw параметра [17]. Результат скрининга.

Малышево Шилово

Гремячье

Масловка г. N5T32.9941 \ Е39Ч 6.859' Ф

Боево

• 049

116

0 т Колодезный 1

Россошь i и

НОВОВОРОНЕЖ

336^

FI <,1

10<FK15 ШШЯ FI>J5

11 6 Фскважина

Сторожевое ([3 область СЭ

РИС. 10.4 б Обобщение результатов скрининга оценок FI параметра по участку Воронежского кристаллического массива (ARu, PRi), перекрытому осадочными юродами (D, К, Pg, N, Q) с суммарной мощностью до 70 м [6] глубины 100м обнаружены два изолированных водоносных горизонта. Деструктивные изменения в разломных зонах вблизи позиции 1-19 выражены в форме разрыва и ступенчатого смещения кровли (подошвы) жесткого основания и в облекающих формах песчано-глинистых слоев, лесса. Состав осадочного покрова, альтитуды и геометрия границ отдельных его составляющих весьма неустойчивы даже на дистанциях, не превышающих 3 км. Альтитуды поверхности наблюдений имеют выраженную тенденцию роста с юга на север (от 150 до 205 м). Для правобережья р. Дон характерно развитие 15° оползневых склонов и оползней нескольких возрастных генераций площадью до нескольких квадратных км, существенно усложняющих оценку устойчивости формации. В теле оползней обнаруживаются очаги новейших подвижек. База оползней — на десятки метров ниже уровня реки. В попытке установить причины, механизмы и потенциальные последствия активизации оползней были изучены относительная микросейсмическая активность и ее связь с разрывными нарушениями и с гидрогеологическими характеристиками. Вследствие существования нескольких плоскостей легкого скольжения практически равновозможны альтернативные сценарии разделения оползней на самостоятельные блоки.

Каков же полезный информационный прирост? Во-первых, вблизи скважин 116 и 7801 с почти одинаковыми суммарными мощностями покровных отложений получены следующие также почти совпадающие оценки HVSR параметра: 2,2 и 2,1. Грубая интерпретация, этого факта позволяет заключить, что относительные изменения мощности покровных отложений и HVSR параметра коррелированы и часто почти совпадают по уровню. Во-вторых, вблизи двух подозреваемых разломных. зон обнаружены (как и ожидалось согласно теоретическим и модельным предпосылкам) аномально низкие значения FI параметра: ,07; ,09. В то же время над консолидированной формацией, подтвержденной скважинными данными, оценка этого параметра почти вдвое выше (,14). Здесь необходимо учитывать, что дистанция до границы: консолидированный блок - разломная зона имеет порядок длины волны (до 3 км), доминирующей в составе микросейсм. Таким образом, FI параметр как индикатор зоны глубинного разлома эффективен. В третьих, на фоне общей тенденции спада уровня в направлении на юг вблизи всех ранее установленных разломных зон имеют место локальные минимумы оценок параметра lg Е, составляющие по отношению к оценкам в бортовых зонах соответственно: ,4-,7. Это дает определенные, основания для утверждений об однотипности разломных зон, а также для распознавания участков в границах разломных зон, связанных с оползнями и с потенциально опасными СЭ (Антонов Ю.В. и др., 2005). Повышение геодинамической активности изменяет водонасыщение тела оползня: обычно оно локально растет. Вблизи зон активной геодинамики формируются границы новых блоков (и/или изменяется геометрия границ уже существующих), трещины отрыва, то есть имеет место процесс подготовки потенциальных срывов оползневых масс.

11. Заключение, выводы и рекомендации

Теоретически вынужденное движение поверхности геологической формации тесно связано с известным видом неоднородного дифференциального уравнения второго порядка с постоянными коэффициентами (к - упругости и h - динамической вязкости). Тремя такими уравнениями (в общем случае с различными к и h) характеризуется разложение вектора микросейсмического движения на три ортогональные компоненты по координатным осям. Как результат, покомпонентные микросейсмические спектры (амплитудно-частотные и энергетические) существенно различаются в четырех-пяти октавной области относительных частот со/к > ,5. Область выраженного спектрального максимума надежно устанавливается по доминирующему периоду АКФ для ряда оценок кинетической энергии (lg Е). На спектры, сохраняющие в целом форму, продиктованную заданным видом уравнения движения, накладываются возмущения, особенно существенные в области со/к>1, связанные с реальной интерференцией волн вблизи акустически контрастных границ внутри среды и поверхности. Теоретические предпосылки и модели нашли многократное и серьезное практическое подтверждение при анализе продолжительных ЗС микросейсмических выборок, полученных, начиная с 1994 года, в разнообразных сейсмогеологических условиях платформ, бассейнов, кристаллического щита, краевого прогиба, орогена, в России и в Румынии. Эти же предпосылки и модели использованы при анализе динамики реальных рукотворных конструкций - мостов и многоэтажных зданий.

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

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

2. разработаны основы и развиты (и/или адаптированы) прикладные программные пакеты обработки и интерпретации ЗС микросейсмических данных "и производных параметров для идентификации деструктивных процессов в реальных геологических средах и рукотворных конструкциях, для обнаружения погребенных разломов и зон деструктуризации, а также изучения их современного состояния и активности, для идентификации реальных сред с принципиально различными механизмами отклика на возмущение: преимущественно упругим, преимущественно вязким или вязкоупругим,*

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

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

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

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

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

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

Содержание аэрозолей в атмосфере, например, сегодня и 150 лет назад характеризуется двумя следующими точными цифрами: 356 и 280 частиц на миллион соответственно. Среднее содержание в атмосфере долгоживущей двуокиси углерода сравнительно с доиндустриальным периодом (завершившимся 250 лет назад) характеризуется сегодня превышением на 31% и продолжающимся экспоненциальным ростом. Скорость этого роста координирована с глобальным потеплением, составляющим 0,7 град Цельсия за последние сто лет (по независимым отечественным данным 0,45±0,15 град, К.Я. Кондратьев) и с прогнозируемым 2-6 градусным потеплением в следующем столетии, а также с ростом аллергических и легочных заболеваний, с сокращением видового разнообразия флоры и фауны, особенно в Северном полушарии, с сокращением эффективных площадей для земледелия, с ростом количества и мощности экстремальных событий. Прогнозы обещают удвоение названного превышения в течение ближайших 50-100 лет (Беклемишев А.Б., 2003).

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

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

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

Но эти сценарии и прогнозы учитывают только три из четырех названных составляющих. Первая же из них и, повидимому, наиболее существенная, остается пока вне внимания экспертов. Это так, поскольку ни один из экспертов (J. Shinn, К. Trenberth, R. Repetto из Национального центра атмосферных исследований, L. С. Gerhard из геологической службы штата Канзас, США) не приводит данных относительно доли естественных парниковых газов в суммарном тренде ВПГ (по информации из AAPG Explorer, 2001-2002).

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

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

Нельзя сказать, что данные о естественном тренде ВПГ полностью отсутствуют или нет инструментов для его оценки, хотя бы косвенной или качественной. Через геохимический и тепловой мониторинг (time lapse monitoring-TLM), в том числе использующий средства газовой хроматографии над глубинными разломами, идентифицированными и позиционно прослеженными на земной поверхности прецизионными геофизическими и GPS методами, через измерения плотности вертикальных электрических токов в атмосфере как функции эндогенного аэрозольного и/или радонового потока (норма - единицы мкА/кв. км) - эти данные могут быть получены повсеместно. Разлом же является, как известно, наиболее газопроницаемым и теплопроводным (благодаря прежде всего незатрудненной циркуляции грунтовых вод) литосферным объектом.

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

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

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

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

Библиография Диссертация по наукам о земле, доктора технических наук, Беклемишев, Андрей Борисович, Москва

1. Аббасов Т.А., Беклемишев А.Б. и др. Устройство для накапливания сейсмических сигналов, Авт. свид. СССР № 354380, Бюлл. изобр. № 30, 1972

2. Адушкин В.В., Спивак А.А. Геомеханика крупномасштабных взрывов, Недра, М., 1993, 319с.

3. Алампиев В.П., Беклемишев А.Б. Адресация элементов в жидкокристаллической матрице, Приборы и системы управления, № 8, изд. Машиностроение, М., 1974, с. 44-45

4. Александров С.И. Поляризационный анализ сейсмических волн, М., ОИФЗ РАН, 1999

5. Александров С.И., Гамбурцев А.Г., Кузнецов В.В. Использование параметров поляризации сейсмических волн при изучении неоднородных сред, Докл. АН СССР, т. 250, № 5, 1980

6. Арутюнов С.Л., Давыдов В.Ф., Кузнецов О.Л., Графов Б.М., Сиротинский Ю.В. Явление генерации инфразвуковых волн углеводородной залежью, диплом № 109, Научные открытия, Сб. кратких описаний, вып. 1, М., 1999

7. Ахулков С.Е., Бархоткин В.А., Беклемишев А.Б., Брехов Р.С., Крыликов Н.О. Реализация требований к методике и технике полевых микросейсмических измерений, Оборонный комплекс научно-техническому прогрессу России, № 3, М., 2005, с. 60-63

8. Ахулков С.Е., Беклемишев А.Б., Крыликов Н.О. и др. Прибор для оперативных полевых исследований вибросейсмического излучения, Конверсия, вып. 11, изд. «Изана», М., 1995, с. 51-52

9. Ахулков С.Е., Беклемишев А.Б., Клинов В.В., Крыликов Н.О. и др. Вибротестер для оперативного полевого контроля параметров сейсмических виброисточников, Электронная промышленность, вып. 7, изд. "ЦНИИЭлектроника", М., 1995, с. 75-76

10. Беклемишев А.Б. Геодинамика двух зон коровых разломов (IMF и TPVF) Восточной Румынии по данным ЗС микросейсмического скрининга 1999 и 2004 годов: сравнительные оценки и тенденции, Геоинформатика, № 3, М., 2005, с. 43-48

11. Беклемишев А.Б. Прогноз потенциальной эффективности сейсморазведки через опережающий ЗС микросейсмический скрининг, Геоинформатика, № 1, М., 2005, с. 35-39

12. Беклемишев А.Б. Роль параметров поляризации в технологии пассивной сейсморазведки, Известия ВУЗ-ов: «Геология и разведка», № 5, М., 2005, с. 37-42

13. Беклемишев А.Б., Кудинов В.Г. Диагностика и дефектоскопия Земли как природной тепловой машины, Проблемы машиностроения и надежности машин, № , М., 2005 (в печати)

14. Беклемишев А.Б., Потапов О.А., Ларин Г.В., Буданов В.Г. Технология пассивной сейсморазведки при оценке состояния бассейнов, глубинных разломов и их эволюции, Приборы и системы разведочной геофизики, 02, изд. СО ЕАГО, Саратов, 2004, с. 24-29

15. Беклемишев А.Б. Штрихи к портрету выдающейся личности, Сб. «М.К. Полшков Воспоминания близких, учеников, коллег, друзей», изд. ВНИИГеофизика, ЕАГО, М., 2004,с. 177-181

16. Беклемишев А.Б., Потапов О. А., Буданов В.Г. Сейсмогравитационный метод: современный статус, Расширенные тезисы докладов Международной геофизической конференции, М., 2003

17. Беклемишев А.Б. Инструментальное изучение погребенных разломов, Доклады Конференции "Строение, живая тектоника и дислокации платформ и их горно-складчатых обрамлений", ИДГ и ИФЗ РАН, М., 2003

18. Беклемишев А.Б. Локализация и активность разломов юго восточных Карпат в свете новых данных ЗС микросейсмического мониторинга и многопараметровой интерпретации, Геофизический вестник, № 6, 2002, с. 8-14

19. Беклемишев А.Б., Потапов О.А. Инструментальное прослеживание геофизическими методами погребенных разломов и зон деструктуризации, Разведка и охрана недр, № 10, 2002, с. 59-63

20. Беклемишев А.Б., Потапов О.А. Оценка возможностей экстремальных событий на действующих и проектируемых трубопроводах, Доклады 4-ых геофизических чтений им. В.В. Федынского, Центр "ГЕОН", М., 2002, с. 104

21. Беклемишев А.Б., Буданов В.Г., Верба M.JL, Локшин Б.С. Изучение активных внутриплитных разломов земной коры с целью снижения экологических угроз, Разведка и охрана недр, № 7-8, 1998, с. 70-71

22. Беклемишев А.Б. и др. К оценке электрооптических параметров матричных жидкокристаллических индикаторов, Электронная техника, сер. «Электровакуумные и газоразрядные приборы», вып. 2(107), 1985, с. 59-64

23. Беклемишев А.Б. и др. Развитие электроэрозионной технологии формирования приводящего рисунка на подложках жидкокристаллических матриц, Электронная техника, сер. «Электровакуумные и газоразрядные приборы», вып. 1(106), 1985, с. 31-33

24. Беклемишев А.Б. Алексеев М.И. Применение электроэрозии при • изготовлении жидкокристаллических матриц, Электронная техника, сер. 4, № 1, 1979, с. 87-90

25. Беклемишев А.Б. Устройство визуализации сейсмозаписей на основе применения жидкокристаллических сред, Сб. «Голография и оптическая обработка информации в геологии и геофизике» под ред. проф. С.Б. Гуревича, Ленинград, 1979, с. 144-152

26. Беклемишев А.Б. и др. Применение методов корреляционного анализа для обработки данных сейсморазведки, Обзор, сер. «Региональная, разведочная и промысловая геофизика», изд. ОНТИ-ВИЭМС, М., 1974, с. 1-58

27. Беклемишев А.Б. Исследование разрешающей способности корреляционной сейсмической аппаратуры, Канд. диссерт., фонды ВНИИГеофизика, М., 1970

28. Беклемишев А.Б., Рапопорт М.Б., Шнеерсон М.Б. и др. Способ обработки сейсмических данных, Авт. свид. СССР № 284333, Бюлл. изобр. № 32, 1970

29. Беклемишев А.Б. и др. Об определении направления на источник возбуждения упругих колебаний, Разведочная геофизика, вып. 4, Недра, М., 1965, с. 3-11

30. Быков И.А. Определение параметров поляризации сейсмических волн, Вопросы динамической теории распространения сейсмических волн, вып. XXIV, Наука, Л., 1984,с. 196-212

31. Воюцкий B.C., Беклемишев А.Б., Барсов В.В. Способ получения спектров сейсмических сигналов, Авт. свид. СССР № 258637, Бюлл. изобр. № 1, 1970

32. Гамбурцев А.Г. и др. Атлас временных вариаций природных, антропогенных и социальных процессов, т. 3,2002, 652 с.

33. Геофизики России, Информ биогр. сб., изд. ЕАГО, М., 2001

34. Голицын Б.Б. Лекции по сейсмометрии, Петербург, 1912

35. Гольдин С.В. Линейное преобразование сейсмических сигналов, Недра, М., 1974

36. Гольцман Ф.М. Статистические модели интерпретации, Наука, М., 1971

37. Гурвич И.И. Сейсмическая разведка, Изд. нефтяной и горно-топливной лит., М., 1960, 504 с.

38. Гурвич И.И. Об отражениях от тонких пластов в сейсморазведке, Прикл. геофизика, вып. 9, Гостоптехиздат, М., 1952

39. Добрев Т.Б., Щукин Ю.К. Геофизические поля и сейсмичность восточной части Карпато-Балканского региона, Наука, М.,1974

40. Дьяконов Б.П. и др. Результаты исследований сейсмоакустической эмиссии пород в Уральской и Кольской сверхглубоких скважинах, Расширенные тезисы докладов Международной геофизической конференции, М., 2003

41. Епинатьева A.M. О суммарных многократных отраженных волнах, Изв. АН СССР, сер. геофиз., № 8, 1959

42. Журков С.Н. и др. Концентрационный критерий объемного разрушения твердых тел, сб. «Физические процессы в очагах землетрясений», Наука, М., 1978

43. Ивакин Б.Н. Головные, проходящие и другие волны в случае тонкого твердого слоя в жидкости, Тр. Геофиз. инст. АН СССР, № 35 (162), 1956

44. Инструкция по забивке свай вблизи зданий и сооружений ВСН 358-76, Минмонтажспецстрой СССР, М., 1976

45. Канасевич Э.Р. Анализ временных последовательностей в геофизике, Недра, М., 1985, 400 с.

46. Кейлис-Борок В.И. Интерференционные поверхностные волны, Академиздат, М., 1960

47. Козлов Е.А. Распознавание и подавление многократных волн в сейсморазведке, Недра, М., 1982

48. Кристенсен Р. Введение в теорию вязко-упругости, Мир, М., 1974

49. Кудинов В.Г. Отработка методов и средств оценки виброакустических свойств роторных и зубчатых механизмов, Доклады Пятой Международной конференции по проблемам вибрации, ICOVP-2001, ИМАШ, М., 2001

50. Кудинов В.Г., Генкин М.Д. Пьезоэлектрический датчик, Авт. свид. СССР № 542955, Бюлл. изобр. № 19, 1972

51. Лугвин В.Г. Элементы современной низкочастотной электроники, Энергия, М-Л., 1964, 88 с.

52. Максимов Л.С., Шейнин И.С. Измерение вибраций сооружений, Л., Стройиздат, 1974

53. Минорский В.П. Сборник задач по высшей математике, ГИТТЛ, М., 1955

54. Михальцев А.В., Мушин И.А., Погожев В.М. Обработка динамических параметров в сейсморазведке, Недра, М., 1990., 250 с.

55. Невский М.В. Исследование современных геодинамических процессов сейсмическими методами, Диссертация, ОИФЗ, 1996

56. Никитин А.А., Тархов А.Г. Комплексирование геофизических методов, М.,1982

57. Никитин А.А. Статистические методы выделения геофизических аномалий, Недра, М., 1979,342 с.

58. Николаев А.В. Динамические процессы в геофизической среде, Наука, М., 1994

59. Николаев А.В. Физические основы сейсмического метода. Нетрадиционная геофизика, М.,1991

60. Николаев А.В. Сейсмические волновые поля в зонах разломов, М., 1978.

61. Николаевский В.Н. Механика пористых и трещиноватых сред, Недра, М., 1984

62. Облогина Т.И. О дифрагированных сейсмических волнах, Прикладная геофизика, вып. 20, Гостоптехиздат, М., 1958

63. Орлов А.Н. Введение в теорию дефектов в кристаллах, Высшая школа, М, 1983, 144 с.

64. Пасечник И.П. Сравнение результатов теоретического и экспериментального исследования резонансных явлений в системе почва сейсмограф, Изв. АН СССР, сер. геофиз., № 4, 1952

65. Полшков М.К., Беклемишев А.Б., Барсов В.В. Применение серийной корреляционной аппаратуры для оценок спектральных характеристик сейсмических колебаний, Прикладная геофизика, вып. 69, Недра, 1972, с. 35-45

66. Потапов О.А., Безрук И. А., Беклемишев А.Б. Новые геофизические технологии, основанные на взаимодействии полей, Тезисы докладов 4-х геофизических чтений им. В.В. Федынского, изд. Центра ГЕОН, М., 2002, с. 26

67. Потапов О.А., Беклемишев А.Б., Будагов А.Г., Гродзенский В.А., Шехтман Г.А. Новые направления в развитии наземной, скважинной и дистанционной сейсморазведки, Прикладная геофизика, вып. 131, Недра, М., 1994, с. 41-56

68. Пронин А.П. Применение гелиевого метода при геоэкологических исследованиях и охране недр, изд. Геоинформмарк, 1997, 72 с.

69. Пронин А.П., Башорин В.Н. Современная флюидная активность на Русской платформе: экологические аспекты, Сб. "Геоэкологические исследования и охрана недр", № 2, изд. Геоинформмарк, М., 1996, с. 3-7

70. Пузырев Н.Н. Методы и объекты сейсмических исследований. Изд. НИЦ ОИГГМ СО РАН, Новосибирск, 1997

71. Работнов Ю.Н. Элементы наследственной механики твердых тел, Наука, М., 1977

72. Ребиндер П.А. и др. Поверхностные явления в твердых телах в процессах их деформации и деструкции, Успехи физических наук, т. 108, вып.1, 1972

73. Саваренский Е.Ф., Кирнос Д.П. Элементы сейсмологии и сейсмометрии, Гостеоретиздат, М., 1955

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

75. Свешников А.А. Прикладные методы теории случайных функций, Судпромгиз, Л., 1961, 252 с.

76. Смирнов Н.В., Дунин Барковский И.В. Краткий курс математической статистики для технических приложений, Физматгиз, М., 1959, 436 с.

77. Соболев Г.А. Оценка сейсмической опасности и сейсмического риска, М.,1997

78. Соболев Г.А. Физические основы прогнозирования разрушения горных пород при землетрясениях, М., 1987

79. Солодилов Л.Н. Центр ГЕОН 25 лет глубинных сейсмических исследований, Разведка и охрана недр, № Ю, М., 1994

80. Страхов В.Н. Вклад Г.А. Гамбурцева в теорию интерпретации гравитационных и магнитных аномалий, Сб. "Развитие идей Г.А. Гамбурцева в геофизике", Наука, М., 1982,с. 17-89

81. Урупов А.К., Кондратович Ю.В. Основы теории комбинированной обработки спектров сейсмических записей для определения свойств геологического разреза, Прикладная геофизика, изд. ВИЭМС, М., 1976

82. Франк C.JI. Сб. "Русская религиозная антропология", т. 11, изд. МДА, 1997

83. Харкевич А.А. Спектры и анализ, Гостехтеоретиздат, М., 1956., 235 с.

84. Харкевич А.А. Теория электроакустических аппаратов, ОНТИ, М., 1940

85. Харкевич А.А. О применении критерия квадратичной погрешности линейных искажений, ЖТФ, вып. 5, 1937

86. Яницкий И.Н. Гелиевая съемка, Недра, М., 1979, 96 с.

87. Aki K. Local site effects on strong ground motion, Proc. of Earthq. Engineering and Soil Dynamics 11, GT Div/ASCE, Park City, Utah, 1988

88. Aki K. and Richards P.G. Quantitative seismology: Theory and Methods, Freeman and Co, 1980

89. Aki K. Space and time spectra of stationary stochastic waves with special reference to microtremors, Bull. Earthq. Res. Inst., 35, 1957, p. 415-456

90. Angelier I. Tectonic analysis of fault slip data sets, Journ. Geophys. Res. 89, В 7, 1984, p.5835-5848

91. Angot A. Complements de mathematiques, Paris, 1957, 772 p.

92. Bard P.Y. Microtremor measurements: a tool for site effect estimation? The effects of surface geology on seismic motion, Eds: Irikura et al., Balkema, Rotterdam, 1999, p. 1251-1279

93. Bard P.Y. and Bouchon M. The seismic response of sediment filled valleys, (a) Part 1, the case of incident SH waves, Bull. Seism. Soc. Am. 70, 1980, p. 1263-1286; (b) Part 2, 1985, p. 519541

94. Bartlett M.S. Smoothing periodograms from time series with continuous spectra, Nature, 161, 1948

95. Beklemishev A.B. et al. Active deep faults of the Vrancea seismogenic area (Romania): results of recent multidiscipline measurements, Article is submitted for publication to the Revue Roumaine de Geophysique, Bucharest, 2000, p. 1-10

96. Beklemishev A.B. et al. Gravity microseismic approach to prediction of buried active faults, Annales Geophysicae Supplement, vol. 16, Proc. of XXI11 General Assembly of EGS, Nice, 1998

97. Beklemishev A.B. et al. Microtremor spectrum dynamics, geochemical and gravity data in the monitoring of the intraplate fracture zones, Annales Geophysicae, Supplement to vol. 16, Proc. of XXI11 General Assembly of EGS, Nice, 1998

98. Beklemishev A.B. Natural three component vibrations. What good is it? , Extended Abstract Book of the Cairo' 98 Africa/Middle East Second International Geophysical Conf. and Exp., 1998

99. Beklemishev A.B. et al. Instrumental confirmation and using of reading increments at different type gravimeters induced by seismic events, Annales Geophysicae, Supplement 1 to vol. 15, Proc. of XXI1 General Assembly of EGS, Vienna, 1997, p. С 133

100. Beklemishev A.B. et al. Structural dynamic approach to complex study of Russian platform basement and sediments, Extended Abstract Book of the 59 th EAGE Conf., vol. 2, Geneva, 1997,p. D012

101. Beklemishev A.B. et al. Forecast of extreme events on operating and projecting pipelines, Extended Abstract Book of the Istanbul' 97 International Geophysical Conf. and Exp., 1997, p. 97

102. Besutiu L. Vrancea active seismic area: a continental unstable triple junction? Revue Roumaine de Geophysique, tome 45, Editura Academiei Romane, 2001

103. Boatwright J. et al. Ground motion amplification in the Marina district, Bull. Seism. Soc. Am., 81, 1991

104. Borcherdt R.D., Effects of local geology on ground motion near San Francisco Bay, Bull. Seism. Soc. Am., 60, 1970

105. Diaconescu M. et al. Distribution of earthquakes from Vrancea area, Supplement to Romanian Geophysics, International geophysical conference, Conference Volume, Bucharest, 2000

106. Dobrin M.B. et al. Surface and near surface waves in the Dalaware basin, Geophys., № 4, 1954

107. Field E.H., Jacob K.H. A comparison and test of various site response estimation techniques including three that are non reference site dependent, submitted to Bull. Seism. Soc. Am., 1994,p. 1127-1143

108. Friedrich A. et al. Ocean generated microseismic noise located with the Grafenberg array, Journ. of Seismology, 2, 1998, p. 47-64

109. Gephart J.W. and Forsyth D.W. An improved method for determining the regional stress tensor using earthquake focal mechanism data: Application to the San Fernando Earthquake Sequence, Journ. Geophys. Res., 89, В 11, 1984, p. 9305-9320

110. Gumbel E. J. Statistical Theory of Extreme Values and Some Practical Applications, National Bureau of Standards, Washington, 1954

111. Gutenberg B. The amplitudes of waves to be expected in seismic prospecting, Geophys., № 2, 1936

112. Irikura K. and Kawanaka T. Characteristics of microtremors on ground with discontinuous underground structure, Bull. Disaster Prev. Inst., Kyoto Univ., 30, 1980, p. 81-96

113. Kanai K. Engineering seismology, Univ. of Tokyo Press, 1983, 251 p.

114. Kelly D.W. The manual of free energy devices and systems, Burbank, CA, 1991

115. Korteveg D.J. and de Vries G. On the change of form of long waves advancing in a rectangular canal and on a new type of long stationary waves, Philosophical Magazine 39, 1895

116. Kudinov V.G. and Genkin M.D. Pickup for measuring vibration parameters of operating machinery parts, Pat. USA № 4021688, 1977

117. Landee R.W. et al. Electronic Designers' Handbook, 1957, 704 p.

118. Lermo J. and Chavez Garcia F.J. Site effect evaluation using spectral ratios with only one station, Bull. Seism. Soc. Am., 83,1993, p. 1574-1594

119. Mc Kenzie D.P. The relation between fault plane solutions for earthquakes and the directions of the principal stresses, Bull. Seism. Soc. Am., 59, 1969, p. 591-601

120. Metric Catalog Supplement № 400, Newport Corp., 1985, 111 p.

121. Moczo P. et al. Amplification and differential motion due to an antiplane 2d resonance in the sediment valleys embedded in a layer over the half space, Bull. Seism. Soc. Am., 60,1996,p. 1434-1446

122. Moczo P. and Bard P.Y. Wave diffraction, amplification and differential motion near strong lateral discontinuities, Bull. Seism. Soc. Am., 83, 1993, p. 85-106

123. Moriya H. et al. Precise estimation of acoustical emission source direction by spectral matrix analysis, Progress in Acoustic Emission, V, The Japanese Soc. for NDI, 1990, p. 244-251

124. Nakamura Y. A method for dynamic characteristics estimation of subsurface using microtremor on the ground surface, QR of R.T.R., 1989

125. Polonic G. and Zugravescu D. Geodynamic compartments and present day stress state of the Romanian territory, Revue Roumaine de Geophysique, 1999, 32 p.

126. Prescott W.H. Circumstances surrounding the preparation and suppression of a report on the 1868 California earthquake, Bull. Seism. Soc. Am., 72, 1982, p. 2389-2393

127. Sammis C.G., Aki K. and Kawase H. Damage patterns due to the Whittier Narrows earthquake of October 1987: a ring of destruction by SV waves at critical incidence, submitted to Bull. Seism. Soc. Am., 1987

128. Samson C. Matrix and Stokes vector representation of detector for polarized waveforms. Theory with some applications to teleseismic waves, Geophys. Journ. Astr. Soc., 51, 1977

129. Sanchez Sesma F.J. Diffraction of elastic waves by three - dimensional surface irregularities, Bull. Seism. Soc. Am., 75, 1985, p. 1435-1446

130. Seo K. et al. Study of site effects in Kobe area using microtremors, X World Conf. Earthq. Eng. Acapulco, Elsevier Science Ltd., 1996

131. Sharpe J.A. The production of elastic waves by explosion pressures, Geophys., № 2, 3, 1942

132. Stotter C. Numerical Modeling of Seismic Noise in Canonical Structures, Diplomarbeit, Univ. Wien, Formal- und Naturwissenschaftliche Facultat, Institut fur Meteorologie und Geophysik, Wien,2003

133. Основные результаты диссертации опубликованы в ведущих научных журналах и изданиях (см. вышеприведенные позиции 3, 8, 12-14, 21, 26-28) согласно Перечню и Дополнениям к Перечню, утвержденным Президиумом ВАК Минобразования России в2004 году.13. Об авторе

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

135. Диссертант автор и соавтор более чем 130 научных работ и изобретений, из которых более 95 опубликованы (более 40 по теме диссертации, см. раздел 12).

136. Перечень иллюстраций и подписей к иллюстрациям

137. РИС. 2.1 стр 25 Амплитуды смещения на свободной поверхности при падении плоской SV волны. Коэффициент Пуассона равен ,25 141.

138. РИС. 2.3 стр 26 Траектории частиц поверхностной волны Релея в функции глубины (а) и зависимости (б) нормированных амплитуд вертикальной и горизонтальной компонент смещения от глубины 43,124.

139. РИС. 2.4 стр 26 Пример кривых эллиптичности (H/V) для фундаментальной и четырех первых мод волн Релея на поверхности стратифицированного полупространства. Острые пики соответствуют исчезновению вертикальной, минимумы горизонтальной компоненты 99.

140. РИС. 2.6 стр 33 Эмпирические отношения H/V как функции частоты для: (а) мягких (осадочных) и (б) жестких (коренных, кристаллических) приповерхностных сред 99.

141. РИС. 3.1 стр 44 Влияние сейсмогеологических условий и условий установки сейсмоприемника на временную и спектральную реакцию 70.

142. РИС. 3.2 стр 44 Подошва (кружки) акустически мягкого осадочного чехла, перекрывающего погребенные коренные породы каменноугольного возраста, осложненные надвигом.

143. Выявлена посредством анализа микросейсм в 19 позициях (закрашенные треугольники). Данные бурения в 10 позициях (вертикальные отрезки) приведены для сравнения 125.

144. РИС. 3.3 стр 47 Импульсное возмущение (Е) и реакция среды (V) как соответствующие скорости вертикального смешения: а во временной области, б - в частотной области; в -передаточная функция системы в вертикальном направлении как V/E 134.

145. РИС. 4.1 стр 54 АЧХ (со/к) для различных значений относительного коэффициента демпфирования 131.

146. РИС. 4.2 стр 56 ФЧХ (со/к) для различных значений относительного коэффициента демпфирования 131.

147. РИС. 4.3 стр 56 АЧХ для системы с отрицательной обратной связью для пяти отмеченных значений (от 0 до -63) параметра Ар 131.

148. РИС. 4.5 стр 66 К описанию пространственной кривой (траектории) через тангенциальный, бинормальный и главный нормальный векторы. Единичные векторы для соответствующих направлений: x,p,v (см. текст).

149. РИС. 4.6 стр 66 Результат прямого Фурье-преобразования для бесконечного и ограниченного во времени (отрезка) гармонического процесса

150. РИС. 4.8 стр 72 Эмпирические "спектральная когерентность (а) и вклад (б)", а также прирост точности оценки азимута источника возмущения (в) в технологии 137.

151. РИС. 4.13 стр 86 Позиции, профили и обобщенные данные микросейсмических и отчасти гелиеметрических измерений 1999 года в пределах четырех (1-1V) блоков геодинамического полигона Vrancea (Румыния) на ранее установленной разломной сети (TPV, NNS, TP).

152. РИС 5.1 стр 91 Частотная зависимость чувствительности компонент акселерометра АПТ-1м.

153. РИС. 5.2 стр 93 Общий вид приборов типа ДВЗ-024 (а) и А 1632 (б).

154. РИС. 5.3 стр 93 Пьезоэлектрический (механоэлектрический) преобразователь мостового типа (б) и его эквивалентная схема (а).

155. РИС. 5.4 а стр 97 Полевая оперативная оценка ускорения плиты (УСК 1), инертной массы (УСК 2), силы при изменяющейся частоте (ЧСТ) развертки в графическом (а) и в табличном (б) виде (вывод обеспечен рулонным принтером в составе системы "Вибротестер").

156. РИС. 5.4 б стр 98 Полевая оперативная оценка фазового сдвига (ФАЗ ГРД) излучаемой гармоники относительно опорного (ОП) сигнала, стабильность которого контролируется, и коэффициента нелинейных искажений (КНИ %) в графическом (а) и в табличном (б) виде.

157. РИС. 5.5 стр 100 Общий вид полевой системы "Вибротестер МС"8.

158. РИС. 5.6 стр 102 Полевые распечатки синхронных микросейсмических временных рядов X,Y,Z компонент (а) и соответствующих АЧХ (б).

159. РИС. 6.1 стр 110 Распечатка экрана монитора с результатами предобработки.

160. РИС. 6.2 б стр 115 Структурно геологическая и тектоническая подосновы ЮВ периферии Карпатского орогена, включая Северо - Добруджинский ороген (см. текст).

161. РИС. 7.1 стр 127 Гистограммы (г-е) и АКФ (а-в, ж) временных рядов FI параметра.

162. РИС. 7.2 стр 130 Гистограмма (а) и АКФ (б) временного ряда lg Кр/К1 параметра.

163. РИС. 7.3 а стр 131 Эмпирические распределения позиций на оси частот главной и двух следующих (по значимости) доминант АЧХ (1-3) для X компоненты и соответствующие АКФ (4-6)

164. РИС. 7.3 б стр 132 То же для Y компоненты (главная и следующая по значимости доминанта АЧХ)

165. РИС. 7.3 в стр 133 То же для Z компоненты (главная и две следующих по значимости доминанты АЧХ).

166. РИС. 7.4 стр 137 Одномерный вариационный ряд FI параметра 16.

167. РИС 8.1 стр 145 Гистограммы и АКФ временных рядов FI параметра в консолидированной зоне.

168. РИС 8.2 стр 146 Гистограмма (а) и АКФ (б) временного ряда lg Кр/К1 параметра в консолидированной зоне.

169. РИС. 8.3 а стр 147 Эмпирические распределения позиций на оси частот главной и следующей (по значимости) доминанты АЧХ для X компоненты и соответствующие АКФ в консолидированной зоне

170. РИС. 8.3 б стр 149 То же для Y компоненты (главная и следующая по значимости доминанта АЧХ)

171. РИС. 8.3 в стр 150 То же для Z компоненты, для главной (1) и для следующей по значимости доминанты (2) АЧХ.

172. РИС. 8.4 стр 151 Распределения оценок FI и lg Кр/К1 параметров, соответствующие АКФ в позиции 1 консолидированной зоны по другую сторону TP разлома.

173. РИС. 8.5 стр 152 Распределения оценок FI параметра, соответствующие АКФ в позиции 2 консолидированной зоны по другую сторону TP разлома.

174. РИС. 8.6 стр 163 Результаты интерпретации данных микросейсмического скрининга 1999 года (белые линии): инструментально подтвержденная и пространственно скорректированная сеть глубинных разломов.

175. РИС. 8.7 стр 163 Эскиз двух однотипных мостовых конструкций (А и Б), профили и позиции ЗС микросейсмических измерений.

176. РИС. 8.8 стр 170 К оценке степени прочностной деградации мостовой конструкции через рассеяние Кр и К1 параметров.

177. РИС. 9.1 стр 177 Зависимость скорости дислокации от напряжения 69. РИС. 9.2 стр 177 Зависимость энергии тела с трещиной от ее длины. РИС. 9.3 стр 180 Карты разрушения для Ni и Мо [69].

178. РИС. 9.4 стр 183 Временные ритмы солнечной активности и даты 23 землетрясений (М>5) 117.

179. РИС. 9.5 стр 186 Временные ряды полной энергии микросейсмического движения (а, б) и соответствующие АКФ (в, г) в TLM фазах 1999 и 2004 года.

180. РИС. 9.6 стр 187 АКФ для временных рядов FIw параметра в TLM фазах 1999 и 2004 года (а, б); сравнение рассеяния синхронных оценок FIw и FI параметров (в) и сопоставление распределений оценок FI параметра для двух позиций в названных TLM фазах (г, д).

181. РИС. 9.7 стр 190 Временные ряды скорости нагружения и релаксации в TLM фазах 1999 и 2004 года в одной позиции Интрамезийского разлома IMF (а, в) и соответствующие АКФ (б, г).

182. РИС. 9.8 стр 192 Временные ряды полной энергии микросейсмического движения (а, в) и соответствующие АКФ (б, г) в TLM фазах 1999 и 2004 года. Разлом TPV.

183. РИС. 9.9 стр 193 Сопоставление временных рядов полной энергии микросейсмического движения (а, г) в двух позициях и соответствующие АКФ (б, в, д, е) в TLM фазах 1999 и 2004 года. Консолидированная зона на удалении одного км от разлома TPV.

184. РИС. 9.10 стр 195 Сравнительная динамика временных рядов FI параметра для трех позиций вблизи и на кромке разлома TPV в TLM фазах 1999 и-2004 года.

185. РИС. 9.11 стр 196 АКФ временных рядов оценок FIw параметра для тех же трех позиций у разлома TPV в TLM фазе 1999 года.

186. РИС. 10.1 стр 201 Зоны завершенной и текущей деструкции, а также близреликтового статуса осадочной формации бассейна, представленные в форме карты в диапазонах значений FI параметра. Результат трехфазного TLM 13.

187. РИС. 10.2 стр 202 Зоны завершенной и текущей деструкции, а также близреликтового статуса осадочной формации бассейна, представленные в форме карты в диапазонах значений FI параметра. Результат скрининга 13.

188. РИС. 10.3 стр 204 Гидрогеологический разрез: Ижорское плато Финский залив (см. текст).

189. РИС. 10.4 б стр 207 Обобщение результатов скрининга оценок FI параметра по участку Воронежского кристаллического массива (Ari^, Pri), перекрытому осадочными породами (D, К, Pg, N, Q) с суммарной мощностью до 70 м 6.1. Таблицы:

190. Стр. 61 Таблица 4.1 Прямое и обратное преобразования Лапласа

191. Стр. 61 Таблица 4.2 К понятию устойчивости системы

192. Стр. 92 Таблица 5.1 Технические характеристики ЗС акселерометра

193. Стр. 101 Таблица 5.2 Технические характеристики полевой МС системы

194. Стр. 112-114 Таблица 6.1 Состояние электронной коллекции МС данных

195. Стр. 124 Таблица 7.1 База МС данных: X компонента

196. Стр. 125 Таблица 7.2 База МС данных: Y компонента

197. Стр. 126 Таблица 7.3 База МС данных: Z компонента