Слайд 1Сейсмический мониторинг месторождений углеводородов
Семестровый курс
Лекция 1
Слайд 2Содержание курса
Вопросы, решаемые пассивным сейсмическим мониторингом на месторождениях углеводородов. Физика
очага землетрясения. Определение энергии землетрясения. Определение положения сейсмического источника.
Модели
очага. Силовые модели источника. Уравнение нодальных линий. Сейсмический режим. Закон повторяемости Гуттенберга-Рихтера. Графики Беньофа. Фрактальность пространственного распределения очагов землетрясений.
Связь микросейсмичности с флюидодинамическими явлениями. Использование регистрации микросейсмических явлений для определения положения трещин гидроразрыва пласта.
Слайд 3Содержание курса
Определение локальных характеристик напряженного состояния по данным инструментальных измерений
и натурных индикаторов:
определение напряжений по механизмам землетрясений;
определение напряжений по
измерениям в скважинах (гидроразрыв, радиусометрия скважин);
определение палеонапряжений по трещиноватости осадочных горных пород.
Определение полей напряжений по экспериментально определяемым ориентациям главных напряжений. Понятие траекторий главных напряжений (ТГН). Запись уравнений равновесия в форме Ламе-Максвелла; гиперболичность системы уравнений равновесия при заданных ТГН. Определение полей напряжений по заданным ТГН в литосфере произвольной реологии. Уравнения плоской задачи теории упругости, формулы Колосова-Мусхелишвили. Определение голоморфной функции по ее аргументу, оператор Шварца. Определение полей напряжений по заданным ТГН в упругой литосфере.
Слайд 4Содержание курса
Обратная задача сейсмологии в томографической постановке. Межскважинная томография на
проходящих волнах. Общие положения межскважинной сейсмики. Использование математического аппарата томографии.
Области применения. Осложняющие эффекты: искривление лучей и анизотропия.
Профильная сейсмическая томография. Томография, использующая форму записи. Томографическое восстановление сейсмических источников. Поверхностно-волновая томография. Дифракционная томография. Эмиссионная томография
Слайд 5Содержание курса
Применение методов нелинейной динамики для анализа рядов сейсмических наблюдений.
Представление сейсмического режима в фазовом пространстве. Динамические системы и их
устойчивость. Аттракторы, их типы и физический смысл. Фрактальная размерность аттракторов. Бифуркации. Инерциальные многообразия. Параметры порядка. Теорема Такенса. Метод Грасбергера – Прокачи оценки размерности аттракторов. Самоорганизованная критичность. Метод клеточных автоматов.
Методы анализа пространственных структур сейсмичности. Метод «ближайшего соседа». Выявление кластеров событий. Энергетический критерий близости. Пространственно – временная взаимосвязь сейсмических событий.
Практическое применение пассивного сейсмического мониторинга для поиска и оптимизации разработки месторождений углеводородов. Требования к системе сейсмического мониторинга. Выявление активных тектонических разломов. Взаимосвязь между сейсмической активностью и параметрами эксплуатации месторождения нефти. Опыт анализа микросейсмов при поиске месторождений газа.
Слайд 6Лабораторные практикумы
Определение скорости прохождения упругих волн через образцы с различной
пористостью.
Определение положения импульсных волновых источников.
Изучение связи между пористостью и проницаемостью.
Изучение
связи между поровым давлением и микросейсмической эмиссией.
Слайд 7Вопросы, решаемые пассивным сейсмическим мониторингом на месторождениях углеводородов.
Физика очага
землетрясения.
Определение энергии землетрясения.
Определение положения сейсмического источника.
Лекция 1
Слайд 8Пассивный микросейсмический мониторинг на месторождениях углеводородов
1970-е – первые опыты ведения
микросейсмического (МС) мониторинга на месторождениях УВ
1980-90-е годы – ряд проектов
по МС мониторингу, в основном в США, крупнейший – Cotton Valley Consortium Project в 1997, начало коммерческих проектов
2000 – е годы – бурный рост проектов по МС мониторингу, в основном в США, около половины – на Barnett Shale, тысячи ГРП с применением МС мониторинга
Гистограмма изменения результатов запросов по ключевому слову «микросейсмичность» в библиотеках SPE и SEG
Слайд 9Пример регистрации микросейсмических событий при проведении серии гидроразрывов пласта на
разных глубинах в двух скважинах:
а) проекция гипоцентров микросейсмических событий
на поверхность земли;
б) трехмерное изображение гипоцентров микросейсмических разрывов вблизи скважины 8-11
Геофоны расположены в скважине
Связь микросейсмичности с трещинами ГРП
Слайд 10Изменение порового давления сопровождается появлением сейсмических событий, обычно не ощущающихся
на поверхности.
По данным о распространении «облака» индуцированных микросейсмических событий
можно оценить проницаемость коллекторов.
Определение проницаемости по данным микросейсмического мониторинга
Слайд 11Морские буровые платформы – один из наиболее подходящих технологических носителей
системы комплексного многоуровневого мониторинга
Слайд 12Сеть сейсмического мониторинга месторождения Валхал
Life of field seismic array
Extensive
use of Advanced Technology
More than 120 km seismic cables
Covers 45 sq. Km.
2500 4-component sensors
Flexible and cost effective acquisition strategy
Seismic source to be operated by standby vessel
Direct access to data
Transmission of data to shore by fiber cable
Слайд 13Определение положения трещины гидроразрыва пласта
Определение положения фронта заводнения пласта
Определение положения
активных разломов
Определение пространственных и временных параметров естественных деформационных процессов, сопровождающихся
изменением НДС и проницаемости пласта
Определение опасного нарастания техногенной сейсмичности
Различие естественной и техногенной сейсмичности
Вопросы, решаемые пассивным сейсмическим мониторингом на месторождениях углеводородов.
Слайд 14Физика очага землетрясения. Модели очага. Силовые модели сейсмического источника.
Модель упругой
отдачи Рейда
Модель прерывистого скольжения
Новообразованные разломы
Фазовые переходы
Роль флюидов
Слайд 15Физика очага землетрясения. Модели очага. Силовые модели сейсмического источника.
Трещина отрыва
Силовой
диполь – пара сил с моментом
Силовой квадруполь – две пары
сил без момента
Три пары сил
Слайд 16Определение механизмов очагов землетрясений
Квадратичное распределение знаков первых вступлений P-волн
Нодальные линии
Распределение
знаков первых вступлений для трещин отрыва и сдвига.
Распределение знаков первых
вступлений S-волн для диполя с моментом и двух диполей без моментов
Сетка Вульфа
Слайд 17Сетка Вульфа (стереографическая сетка) - проекция меридианов и параллелей сферической
поверхности на плоскость основного меридиана.
Окружность сетки называют ее ОСНОВНЫМ МЕРИДИАНОМ
Точки,
в которых сходятся все меридианы, называются ПОЛЮСАМИ СЕТКИ
Диаметр , проходящий через полюса сетки, называется ОСЬЮ СЕТКИ
Диаметр , перпендикулярный к оси сетки, называется ЭКВАТОРОМ СЕТКИ
Методика построения сетки Вульфа
Построение линий меридианов
Линия меридиана представляет собой дугу окружности, которая проходит через полюса и точку А, лежащую на пересечении экватора и прямой, соединяющей один из полюсов с точкой на основном меридиане, отстоящую от другого меридиана на величину, равную долготе искомого меридиана.
Слайд 18Построение меридиан и параллелей на сетке Вульфа
Слайд 20Характеристики землетрясений
Эпицентр – проекция на поверхность очага землетрясения.
Гипоцентр – положения
очага землетрясения.
Сила (интенсивность) землетрясения – характеристика воздействия сейсмических волн в
баллах. XII-балльная шкала Рихтера (модифицированная шкала Меркалли).
Слайд 21Характеристики землетрясений
Магнитуда землетрясения – характеристика сейсмической энергии в очаге землетрясения,
определяется по записям сейсмических волн. По определению
где A – максимальная
амплитуда в мм записи сейсмометра с собственным периодом колебаний 0.8с, усилением 2800 и коэффициентом затухания 50:1.
Различают магнитуду, определяемую по поверхностным волнам, по объемным волнам, локальную магнитуду.
Слайд 22Характеристики землетрясений
A – амплитуда в мкм, T – период в
сек, Δ – эпицентральное расстояние в град., h – глубина
очага
Энергетический класс – логарифм сейсмической энергии в очаге землетрясения, выраженной в джоулях.
Слайд 23Параметры сейсмического источника
W – объемная плотность упругой энергии, θ –
дилатансия, λ,μ – коэффициенты Лямэ, ε – деформация
E=VW
V – объем
очага lgV=a+bM
L – длина разлома V~L3 M=c+dlgL
Слайд 24Параметры сейсмического источника
Сейсмическая эффективность – отношение излученной в виде сейсмических
волн энергии к общей энергии деформации, обычно около 1%
Скорость распространения
разлома (скорость вспарывания)
Сброшенные напряжения
Сейсмический момент
Слайд 25Сейсмологические методы
Уравнение движения упругой среды
Плоские упругие волны в неограниченной изотропной
среде
Слайд 26Типы волн
Объемные волны
Продольные
(p – волны, primary waves, pressure waves)
Поперечные
(s –волны, secondary waves, shear waves)
Слайд 27Разделение упругих волн на независимо распространяющиеся продольные и поперечные волны
возможно и в случае произвольной (не плоской) волны. Для этого
представим решение уравнения в виде суммы двух частей:
одна из которых представляет решение для волны, в которой не происходит изменения объема, так что
а для другой выполняется условие
Подставим сумму в уравнение движения и применим к обеим частям уравнения операции rot и div, получим два уравнения одинакового вида
Слайд 28Для монохроматической волны вектор смещения можно записать в виде
где
u0 – функция только координат, ω – циклическая частота колебаний
в волне. Тогда
- волновые векторы продольной и
поперечной волн
Слайд 29Поверхностные волны Рэлея
Пусть плоская монхроматическая волна распространяется вдоль поверхности
На поверхности
тела
Т.к. вектор нормали направлен по оси z, получаем
откуда по закону
Гука
Второе из уравнений, в силу независимости всех величин от y, дает
получаем
Т.е. вектор смещения лежит в плоскости, перпендикулярной к поверхности и проходящей через направление распространения волны.
Слайд 30Поверхностные волны Лява
При наличии твердого слоя, расположенного на твердом полупространстве,
возникают волны Лява. Эти волны – чисто поперечные: в них
имеется только одна компонента смещения y, упругая деформация в волне Лява представляет собой чистый сдвиг. Смещения в слое (индекс 1) и в полупространстве (индекс 2) описываются выражениями
смещения в слое распределены по косинусу, а в полупространстве экспоненциально убывают с глубиной
Слайд 31Отражение и преломление сейсмических волн
(Закон Снеллиуса)
Слайд 32Годограф сейсмических волн
Годографом сейсмической волны называется график зависимости времени пробега
волны от источника до приемника волны от эпицентрального расстояния.
Эпицентральное
расстояние – это угол с вершиной в центре шара, которым изображается Земля, и сторонами, образованными радиус-векторами источника и приемника
Слайд 33Определение положения сейсмического источника.
По разности времен прихода P и
S волн
По разности времен первых вступлений на нескольких станциях
Слайд 34Литература
Ландау Л.Д., Лифшиц Е.М. Теоретическая физика, т.VII. Теория упругости.
Родионов В.Н.,
Сизов И.А., Цветков В.М. Основы геомеханики.
Браун Д., Массет А. Недоступная
Земля.
Болт Б. В глубинах Земли.
Жарков В.Н., Трубицин В.П. Физика планетных недр.
Рихтер Ч. Элементарная сейсмология.