Все выпуски

В работе представлен подробный и систематический анализ сходимости метода конечных элементов при моделировании распространения упругих волн в анизотропных средах, обладающих ортотропной симметрией тензора упругости. Пространственная дискретизация основана на применении метода конечных элементов с квадратичными базисными функциями $P_2$ на неструктурированных тетраэдральных сетках, что обеспечивает гибкость при работе со сложными геометриями. Для интегрирования по времени используется схема Ньюмарка с параметрами $\beta = 0,302$, $\gamma = 0,6$; при этом выбор именно этой схемы не оказывает влияния на пространственную сходимость, которая является предметом настоящего исследования.

Анализ сходимости выполнен на тестовой задаче с гладким локализованным начальным импульсом в полупространстве с граничными условиями свободной поверхности. Начальные условия подобраны таким образом, что решение остается гладким, а импульс не достигает границ расчетной области до момента сравнения решений. Это позволяет исключить влияние аппроксимации граничных условий на исследуемую сходимость и сосредоточиться исключительно на пространственной составляющей погрешности. Серия вычислительных расчетов проведена на шести последовательно сгущающихся неструктурированных сетках с числомтет раэдров от 39 728 до 372 904. Сравнение производилось с эталоннымрешением , полученным на существенно более мелкой сетке, что обеспечивает достаточную точность референсных данных. Для корректного вычисления нормы разности решения на всех сетках проецировались на единую грубую пространственную сетку.

Установлено, что порядок сходимости в норме $L_2$ составляет 2,314 при использовании медианной длины ребра тетраэдров в качестве характерного размера сетки и 2,878 при использовании обратного кубического корня из числа элементов. Полученные значения свидетельствуют о хорошей сходимости реализованной численной схемы и подтверждают корректность предложенного подхода.

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

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

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

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