В геофизике есть задача, которая звучит как оксюморон: нужно смоделировать сейсмический отклик от трещины, ширина которой составляет микроны, а длина волны — дециметры. Если подходить в лоб — моделировать трещину как отдельный слой с низким импедансом, — то сетка должна разрешать эту микроскопическую толщину. Шаг по времени падает вместе с шагом сетки, и расчет становится невозможным даже на современных GPU.
Выход — не моделировать трещину как объем, а заменить ее граничным условием.
Линейная модель скольжения (Linear Slip Model, LSM) делает именно это. Трещина становится интерфейсом нулевой толщины, на котором вектор напряжений непрерывен, а вектор смещений может испытывать скачок. Этот скачок линейно связан с напряжением через матрицу жесткости K:
\[
\mathbf{t} = K \Delta \mathbf{u}
\]
При K = 0 получаем свободную поверхность; при очень больших K — практически идеальный контакт. Все промежуточные значения описывают реальные трещины с разной степенью сцепления.
Но как заставить спектральный метод работать с таким разрывным условием?
В классическом SEM на гексаэдральной сетке узлы на границе общие. А здесь нужно, чтобы узел с одной стороны трещины и узел с другой стороны существовали независимо — физически они находятся в одной точке, но численно имеют разные смещения. Решение: строим сетку так, чтобы граница трещины проходила по поверхностям элементов, а затем "разлепляем" узлы, дублируя их на двух берегах.
В слабой формулировке метода Галеркина появляется дополнительный граничный интеграл:
\[
\int_{S_{\text{трещины}}} (\mathbf{K} \Delta \mathbf{u}) N_m \, ds
\]
Ключевое наблюдение: благодаря свойствам ортогональности базисных функций Гаусса–Лобатто–Лежандра этот интеграл дает ненулевой вклад только для одинаковых индексов базисных функций. Вычислительная сложность — линейна по числу узлов на трещине. И, что еще важнее, схема остается полностью явной: не нужно решать систему на каждом шаге.
Валидация
Проверяли на двух тестах. Первый — релаксация ступеньки скорости на интерфейсе. Аналитическое решение для разности скоростей на трещине имеет вид экспоненциальной релаксации с частотой ν = k(Z₁+Z₂)/(Z₁Z₂). Численный расчет на 3D-модели с 7-м порядком аппроксимации совпал с аналитикой с хорошей точностью.
Второй тест — плоская волна, падающая на трещину. Для нормального падения коэффициент отражения и прохождения имеют простой частотно-зависимый вид:
\[
R(\omega) = \frac{Z_1 - Z_2 - i\omega Z_1Z_2/k}{Z_1 + Z_2 - i\omega Z_1Z_2/k}
\]
Расчет показал, что при 7-м порядке и шаге сетки, равном длине волны, численные коэффициенты отражения и прохождения совпадают с теорией вплоть до 900 Гц. Для P-волны допустимый частотный диапазон выше — потому что длина волны больше.
Реальная задача: CO₂ и микротрещины в цементе
Применяли метод к задаче, важной для хранения CO₂. В обсаженной скважине цементный камень может иметь микротрещины (ширина — микроны), через которые газ может мигрировать десятилетиями. Моделировали акустический каротаж: монопольный источник с центральной частотой 10 кГц, 21 кольцо приемников, 8 азимутальных датчиков на каждом кольце.
Исследовали два типа трещин:
- Азимутальные — идут по окружности вокруг скважины
- Радиальные — расходятся от скважины наружу
Варьировали жесткость k от 10⁹ до 10¹² Па/м и угловой размер от 45° до 180°.
Результат оказался нетривиальным. Радиальные трещины не дали заметного отклика на сейсмограмме — при любых параметрах. Азимутальные — напротив, сильно меняли амплитуды псевдо-рэлеевской (pR) и стоунлиевской (ST) волн, хотя на первое вступление (P-волну) влияли слабо.
На полярных диаграммах (рис. 12 в статье) видно, как "октагон" амплитуд вытягивается в сторону трещины. Причем увеличение жесткости и увеличение углового размера дают разные паттерны. Это открывает возможность обратной задачи: по данным акустического каротажа оценивать параметры микротрещины.
Что в итоге
SEM + LSM на GPU — это инструмент, который позволяет моделировать микротрещины там, где другие методы упираются в вычислительную стену. Полностью явная схема, линейная сложность по граничным узлам, высокий порядок аппроксимации — и в итоге возможность за приемлемое время получить сейсмограммы, чувствительные к параметрам трещины. А это уже база для создания датасетов и, возможно, для применения машинного обучения в оценке целостности цементного камня.
