Численное моделирование задач геологии и сейсморазведки. Панкратов С. А., Московский Физико-Технический Институт (ГУ) Вычислительная физика: алгоритмы, методы и результаты. г. Таруса. Октябрь 19-21, 2011

Содержание Численное моделирование в сейсмологии Численный метод Практические задачи Моделирование гетерогенных сред Кластеры трещин и каверн Слои с трещинами Многослойные среды Пористые среды

Численное моделирование в сейсмологии Исследование свойств породы Выявление закономерностей Реализация модели осредненной породы Обратные задачи численного моделирования

Численный метод Сеточно-характеристический методыю Учитывает физические свойства процесса (распространение волн вдоль характеристик). Корректная обработка граничных условий. Схемы 1 и 2 порядка. Гибридная схема. Монотонный метод.

Математическая модель ρ – плотность λ, μ – упругие коэффициенты Ляме v – скорость σ – тензор упругости ε – тензор деформаций Уравнение движения Закон Гука Связь скорости и деформаций

Сеточно-характеристический метод Расщепление по направлениям Гиперболическая система уравнений

Вычислительная сетка Треугольная Локальное измельчение

Практические задачи Единичная трещина или каверна Кластеры трещин и кавверн Многослойная порода с трещинами Многослоная порода Пористая среда

Газонасыщенная трещина

Экранирование Поле скоростей

Газонасыщенная трещина Поле скоростей Отращенная продольная и поперечная волны

Флюидонасащенная трещина

Поле скоростей Отраженная поперечная волна(без продольной) Слабое экранирование

Флюидонасыщенная каверна

Прошедшая волна Отраженная продольная волна Отраженная поперечная волна Поле скоростей

Отражение от нерегулярности

Кластер из каверн

Комплексный кластер

Сравнеие с осредненной моделью Осредненная средаКластер из трещин и каверн

Зависимость амлитуды отклика от количества трещин Вертикальные трещины Наклонный фронт падающей волны Постоянная длина резервуара Количество трещин: 3-21

Резервуар с флюидонасыщенными трещинами

Резервуар с флюидонасыщенными трещинами Количество трещин

Анизотропия отраженных волн

Приповерхностный взрыв

Поперечные волныПродольные волныВолна Рэлея

Двухслоная среда Конус Маха

Многослоная среда (5 слоев) Приемник

Сравнение многослоной породы с осредненным слоем. Вертикальная компонента скорости Отраженные волны от контактов

Многослойная среда (10 слоев) Приемник

Вертикальная компонента скорости на приемнике First layer 23 4 Incident wave

Пористая среда

Пористая среда (10% )

Пористая среда (40%)

Сравнение сред с различной пористостью 10%40%

Спасибо за внимание!