Афонина Е. В., Вишневский Д. М., Горшкалев С. Б., Карстен В. В., Лисица В. В., Чеверда В. А. Институт нефтегазовой геологии и геофизики СО РАН, Новосибирск.

Презентация:



Advertisements
Похожие презентации
Баюк И.О. *, Рыжков В.И. ** * ИФЗ РАН, ** РГУ нефти и газа Использование дипольного акустического каротажа для оценки параметров пор и трещин карбонатных.
Advertisements

Институт вычислительной математики и математической геофизики СО РАН Рис. 1. Численные снимки волнового поля в фиксированный момент времени для сложнопостроенной.
Новосибирский государственный университет Механико-математический факультет Кафедра вычислительных систем Численное моделирования распространения упругих.
Мортиков Е.В. 2 4 апреля 2014 г. НИВЦ МГУ М. В. Ломоносова Лаборатория суперкомпьютерного моделирования природно - климатических процессов ЧИСЛЕННОЕ МОДЕЛИРОВАНИЕ.
Численное моделирование задач геологии и сейсморазведки. Панкратов С. А., Московский Физико-Технический Институт (ГУ) Вычислительная физика: алгоритмы,
Расчеты низкоскоростного режима развития детонации ВВ Бахрах С.М., Володина Н.А., Кузьмицкий И.В., Леонтьев М.Н., Циберев К.В. РФЯЦ-ВНИИЭФ ИТМФ, Саров.
Проект 17: Алгоритмическое и программное обеспечение для моделирования деформации микроразрушенных и пористых сред на многопроцессорных вычислительных.
Разработка вычислительных методов и алгоритмов для численного решения динамических задач механики сплошных сред. Московский физико-технический институт.
1 Параллельный алгоритм расчета трехмерного поля давления при моделировании пространственных теплогидравлических процессов Ю.В. Юдов, А.В. Владимиров ФГУП.
М ОДЕЛИРОВАНИЕ ГАЗОВОГО ПОТОКА ЗА СОПЛОМ Кондаков В.Г. Якутск, ноября 2011 г.
Моделирование динамических задач МДТТ с помощью параллельной версии сеточно- характеристического метода Васюков Алексей Викторович, МФТИ, аспирант Научный.
В.И. Исаев Дисциплина «Интерпретация данных ГИС» АКУСТИЧЕСКИЙ МЕТОД 1.
Параллельные алгоритмы для симплициального подразделения области с итерационным измельчением вблизи границы Кафедра параллельных алгоритмов Математико-Механический.
Геометрическое моделирование трехмерных объектов..
Параллельная реализация экономичных методов параболических задач.
Система видеоуправления Кафедра УиИТС, МГИЭМ (Московского Государственного Института Электроники и Математики) Докладчики: Николаева М.В. Сидоров С.М.
Смешанная модель параллельных вычислений OpenMP&MPI в программе газовой динамики Быков А.Н., Жданов А.С. (РФЯЦ-ВНИИЭФ, Россия) 17 мая 2013 г.
Адаптация комплекса программ M2DGD для работы на МВС с использованием среды параллельного программирования OST Павлухин Павел Научный руководитель: Меньшов.
1 ФГУП «НАУЧНО-ИССЛЕДОВАТЕЛЬСКИЙ ТЕХНОЛОГИЧЕСКИЙ ИНСТИТУТ им. А.П. Александрова» Ю.В. ЮДОВ DIRECT NUMERICAL SIMULATION DNS 5-я международная научно-техническая.
Основные теоремы теории очага землетрясения. Тензор сейсмического момента. Лекция 4.
Транксрипт:

Афонина Е. В., Вишневский Д. М., Горшкалев С. Б., Карстен В. В., Лисица В. В., Чеверда В. А. Институт нефтегазовой геологии и геофизики СО РАН, Новосибирск МОДЕЛИРОВАНИЕ ВОЛНОВЫХ ПРОЦЕССОВ В АНИЗОТРОПНОЫХ УПУГИХ СРЕДАХ

Актуальность Анизотропия в сейсмике: Трещиноватые кабонатные коллекторы; Топкослоистые пачки; Предварительные напряжения (в окрестности скважин, вблизи соляных тел и пр.). Анизотропия может выступать как косвенный признак: Ориентации трещин (AVO, AVA analysis); Концентрации трещин; Состава флюида в трещиноватом пространстве и пр.

Выбор метода

Почему неприменима стандартная схема на сдвинутых сетках

Схема Лебедва

Схема Лебедева и схема на повернутых сетках Схема Лебедева, как минимум, на 30% экономичнее схемы на повернутых сетках!

Параллельная реализация 1D1D 2D 3D Трехмерная декомпозиция расчетной области; Неблокирующие процедуры обмена: Isend, Irecv; Параллельные I/O процедуры: MPI I/O

Эффективность и ускорение

Упрощенная модель Юрубчено-Тохомской зоны Модель - горизонтально-слоистая среда с двумя трансверсально-изотропными слоями, направления осей симметрии в которых горизонтальны и составляют с осью Х углы 30° и 60° м. Дискретизация 1 метр (10 точек на длину волны), Общий объем RAM – 4 Tb!

Упрощенная модель Юрубчено-Тохомской зоны

Упрощенная модель Юрубчено-Тохомской зоны, сейсмограмма VSP Расщепление S-волн Контрастный обмен S-волн

Акустический каротаж Анизотропная вмещающая среда; Цилиндрическая система координат; Применение схемы Лебедева; Параллельная реализация на основе декомпозиции расчетной области; Периодическое азимутальное измельчение сетки для компенсации увеличения шага сетки; Основные особенности:

Акустический каротаж

Гибридный алгоритм SSGS LS Увеличение требований на вычислительные ресурсы в пять раз!

Гибридный алгоритм SSGS LS

Гибридный алгоритм Схема Лебедева Стандартная схема на сдвинутых сетках Группа процессоров для LS Группа процессоров для SSGS

Гибридный алгоритм

Потребности в вычислительных ресурсах 15 км. h=2.5 м. 20 км., h=5 м. Число ядер Изотропная упругая модель (3+9) 12 параметров 4000 Анизотропная упругость (22+36) 58 параметров Гибридный алгоритм 20% анизотропии 6500 Анизотропная вязкоупругость (64+84) 148 параметров 45000

Заключение Разработаны и реализованы алгоритмы моделирования волновых процессов в анизотропных средах для задач поверхностной сейсмики, VSP, cross-well, акустического каротажа. параллельная реализация на основе расщепления по пространственным подобластям; используются неблокирующие процедуры Isend, Irecv эффективность порядка 90 % Разработан гибридный алгоритм моделирования волновых полей в средах, содержащих анизотропные включения, что позволяет экономить до 80% вычислительных ресурсов.

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

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