Трехмерное моделирование акустического поля в скважине, возбуждаемого источником в пласте Текст научной статьи по специальности «Физика»

УДК 534-18

Н. В. Мутовкин1'2, Д. Н. Михайлов2, И. Л. Софропов1'2

1 Московский физико-технический институт (национальный исследовательский университет) 2 Московский научно-исследовательский центр Шлюмберже

Трехмерное моделирование акустического поля в скважине, возбуждаемого источником в пласте

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

Ключевые слова: пассивная акустическая шумометрия, фильтрационный шум, акустический резонанс, метод спектральных элементов.

N. V. Mutovkin1'2, D. N. MikhailotP, I. L. Sofronov1'2

1 Moscow Institute of Physics and Technology 2Schlumberger Moscow Research Center

Three-dimensional modeling of acoustic noise in wellbore induced by the source in reservoir

We consider three-dimensional modeling of wellbore acoustic fields induced by the fluid flow noise in the reservoir. The goal is to develop highly accurate mathematical models to improve the quality of interpretation of acoustic noise logging data for oil and gas producing wells. Using the spectral element method, we simulate acoustic fields in wellbore, analyze their spatial frequency patterns and the fine structure of the resulting 3D resonance modes. In the axisymmetric case we compare the results with the previously found solution simulated by the two-dimensional finite-difference algorithm.

Key words: passive acoustic noise logging, noise induced by fluid flow in reservoir, acoustic resonance, spectral element method.

1. Введение

Измерение интенсивности акустического шума в скважине является одним из действенных методов анализа технического состояния скважины и определения гидродинамических свойств коллекторов. Анализ спектральных характеристик акустических шумов, вызванных фильтрацией флюида в породе пласта (фильтрационные шумы) [1], применяется для обнаружения отдающих и принимающих интервалов пласта и определения типа поступающего в ствол скважины флюида. Спектр фильтрационного шума лежит в полосе 2-20 кГц, что соответствует экспериментальным данным [2]. Развитие математических методов для

© Мутовкин Н. В., Михайлов Д. Н., Софронов И. Л., 2020

(с) Федеральное государственное автономное образовательное учреждение высшего образования

«Московский физико-технический институт (национальный исследовательский университет)», 2020

повышения качества и надежности интерпретации данных шумометрии является актуальной задачей.

В [3] описана численная модель, позволяющая рассчитывать акустические поля, генерируемые фильтрационным потоком в пласте, в предположении наличия вращательной симметрии. По результатам моделирования в спектрах были обнаружены достаточно широкие интервалы резонансных частот, где амплитуда акустического поля существенно возрастает. Был проведен анализ чувствительности этих резонансных мод к геометрии скважины, свойствам пород-коллекторов и фазовому составу флюида, что позволило начать разработку новых подходов интерпретации [4,5].

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

2. Постановка задачи

Рассматривается трехмерная постановка задачи в декартовой системе координат (х, у, х). Область моделирования П ПредСтавляет собой прямоугольный параллелепипед {|х| ^ Ьх/2, 1у1 ^ Ьу/2, ^ Ьг/2}, где ЬХ,ЬУ— длины сторон по осям ж, у и ^ соответственно. Модель включает в себя цилиндрическую скважину радиуса г^ц, заполненную флюидом. Ее ось параллельна оси ^ и ПрОХОдИТ через центр параллелепипеда. Скважину окружает пласт^ С0СТ0Я1цИ^ в 0бщем случае из нескольких ПОперечных слоев ТОдщиной г = 1,... , где N — к0личестВ0 ело ев, ^^ х Н\ = Ьг. В некоторых из слоев пласта располагаются источники акустического шума. На данном этапе разработки модели каждый слой считается сплошным, изотропным и изотермическим.

Конфигурация источника рассматривается в двух вариантах: кольцевой и секториаль-ный. Кольцевой источник предполагается осесимметричным и занимает область по радиусу в околоскважинной зоне пласта, соответствующую кольцу г^ц ^ г ^ г^ц + Аг8ГС (рис. 1), где г = д/х2 + у2 — радиальная координата, Аг8Гс — радиус зоны ИСТОчника. Секториальный источник локализован в секторе кольца того же радиуса и занимает область в\ ^ в ^ #2 (рис. 2), где в — азимутальная координата, #2 — границы углового сектора.

Пласт 1

со

^ I Я5 | СО щж Пласт 2 ьГ*

О Пласт 3

ги>е11 (Дг5ГС| X

1х/2

а)

б)

Рис. 1. Типичная геометрия модели с кольцевым источником и тремя слоями в пласте: а) продольный разрез модели, б) поперечный разрез слоя

Распространение волн в пласте моделируется при помощи линейной системы уравнений упругости, записанной в форме волнового уравнения Навье:

рад*и(8) -V- Т = ¥31С(х,у,г,г),

(1)

где Ь — время; — вектор смещений частиц упругой среды; р8 — плотность упругой среды; Е8ГС — объемный источник, отвечающий за генерацию акустического шума. Тензор Т — отклонение от стационарного значения тензора напряжений упругой среды, определяющееся по закону Гука:

Т = Л (V ■ и(в)) I + 2С ■ 1/2 ^+ ,

где А, С — параметры Ламе для упругой среды, I — единичный тензор второго ранга.

Рис. 2. Геометрия модели с секториальным источником: поперечный разрез слоя Распространение волн в скважине описывается уравнением акустики:

р/ д— к/ = 0,

(2)

Где р — скалярный потенциал, определяющий вектор смещения частиц флюида и(/) = V(p/pf] к/ — модуль всестороннего сжатия флюида; р/ — плотность флюида. В условиях па интерфейсе скважина-порода используется выражение р — для откло-

нения от стационарного значения давления флюида.

Источник звуковых колебаний Е8ГС занимает область в пласте, где происходит приток флюида в скважину. Его интенсивность предполагается неизменной по вертикали и азимуту и обратно пропорциональной кубу радиальной координаты в силу снижения скорости фильтрации по мере удаления от скважины. На достаточном расстоянии от скважины, а также в слоях с низким притоком флюида интенсивность источника задается равной нулю. Таким образом, выражение для источника записывается в виде

ЕбГС У) - ег

Мг™е\\/г)33(Ь), Г^п ^ г ^ + ДГзгс, Zí < 2 < Zí + Щ 0, еке,

Где ег — единичный радиальный вектор цилиндрической системы координат, задающий направление действия источника; А — амплитуда источника; Б(Ь) — сигнал, генерируемый источником; Zí — нижняя граница слоя пласта.

На внешних границах расчетной области при х — ±Ьх/2,шли у — ±ЬУ/2,шли г — /2 используются простейшие неотражающие граничные условия [6]. На интерфейсе между скважиной и пластом заданы следующие условия:

Тгг — р, ТГг = 0;

) — ).

(3)

Здесь Тгг, Тг х, и^^ — компоненты соответствующих величин в цилиндрической системе координат.

3. Метод спектральных элементов

Для численного решения рассматриваемой гиперболической системы (1) - (2) с заданными условиями на интерфейсе (3), однородными начальными условиями и неотражающими граничными условиями применяется метод спектральных элементов [7|, реализованный в открытом коде SPECFEM3D Cartesian [8].

Дадим вариационную формулировку задачи (1) - (3), умножив скалярно (1) на произвольную вектор-функцию w(x,y,z), а (2) — на произвольную скалярную функцию w(x, у, z) и проинтегрировав оба уравнения п о областям us (упругая сред а) и Uf (жидкая среда) по частям:

J psw • df u(s) dw -J w • (T • n) ds + J Vw : T du = J Fsrc • w du;

J pf wd^pdu — J kf w n •V^ds + J kf Vw •V^du = 0.

Uf duf Uf

(4)

Расчетная область П = ш8 иразбивается на некоторое количество пе объемных конечных элементов Пе, е = 1,... ,пе} П = 1 Пе, имеющих форму шестигранников в качестве объемных элементов (рис. 3). Граница расчетной области разбивается на щ поверхностных конечных элементов, которые представляют собой четырехугольники дЬ =

д п = и£= 1д пь.

а)

б)

Рис. 3. Разбиение расчетной области на элементы: а) вид сверху, б) вид сбоку

Описание процедуры интегрирования слагаемых в системе (4) по элементам приведено в приложении. В результате система (4) сводится к системе обыкновенных дифференциальных уравнений относительно значений искомых полей в узлах СЬЬ-сеткн в каждом элементе, см. приложение. Полученная система уравнений может быть записана следующим образом:

М(в) и + С (з)и + К (8)и = ^;

М(/) Ф + С(/) Ф + К(/) Ф = 0, ^

где и — вектор, содержащий значения и(в) во всех узлах сетки, Ф — вектор значений р в узлах сетки. Матрицы М(з), С(з\ К(в) возникают при интегрировании первого, второго и третьего слагаемого соответственно в первом уравнении системы (4), Р — вектор, полученный при интегрировании правой части в (4). Матрицы М(?\ С(/\ К(*) выводятся аналогично из второго уравнения в (4).

Благодаря использованию сетки СЬЬ матрицы М оказываются диагональными, что упрощает явную разностную схему для интегрирования по времени. Мы применяем схему Ньюмарка типа предиктор-корректор второго порядка аппроксимации по времени.

Аппроксимация решения задачи по пространству в методе спектральных элементов определяется двумя параметрами: АИ — характерный размер элемента и щ — максимальная степень интерполяционных многочленов Лагранжа, используемых для представления функции на элементе.

В данной работе при моделировании применялся базис из многочленов Лагранжа с П[ = 4. Соответственно на каждый объемный элемент приходится по 125 узлов сетки СЬЬ, а на каждый поверхностный — 25 узлов сетки СЬЬ (5 по каждому из направлений).

4. Результаты моделирования для случая кольцевого источника 4.1. Параметры моделирования

Рассматривается постановка задачи, представленная на рис. 1, с тремя слоями в пласте (Ж = 3). Геометрические параметры модели следующие: Н2 = И = 0.24 м, = 0.12 м, Аг8ГС = 0.24 м, Ьх = Ьу = Ь = 2.88 м, = 6.24 м, Нг = Н3 = Н = 3 м.

Скважина заполнена газом с параметрами, соответствующими метану при давлении 2 ■ 107 Па и температуре 343 К. Порода представляет собой известняк и моделируется линейной упругой средой. Упругие свойства газа и породы представлены в табл. 1.

Таблица1

Свойства газа и породы

Газ (метан) Порода

Плотность (р), кг/м3 123 2125

Скорость продольной волны (ьр), м/с 535 3953

Скорость поперечной волны (и3), м/с 0 2603

Спектр сигнала S(t), генерируемого источником, является равномерным в диапазоне от 2 до 6 кГц с шагом дискретизации в 50 Гц. Амплитуда источника ^^на 5 ■ 10-4 Па/м.

Построенная конечно-элементная сетка содержит порядка 48 тысяч элементов. «Минимальный» и «максимальный» линейный размер элемента, оцениваемый через кубический корень из объема элемента в скважине, составляют 0.04 м и 0.05 м соответственно. В пласте линейный размер элемента варьируется от 0.04 м вблизи скважины до 0.19 м на внешней границе расчетной области.

Минимальная длина волны Amin гармонического по времени акустического поля в скважине, соответствующая максимальной частоте источника, равна 0.09 м. Таким образом, в скважине на длину волны в среднем приходится около двух конечных элементов. Относительная погрешность решения на данной сетке в Ь2-норме составляет 2% в скважине, что было оценено по дополнительным расчетам при проверке сеточной сходимости.

Время моделирования 100 мс было выбрано из условия выхода на периодический режим колебаний в скважине. Количество шагов по времени составляет примерно 250 000.

4.2. Анализ акустического поля в скважине

В [3,4] представлены результаты исследований спектров акустического давления в скважине в осесимметричной постановке задачи, полученные посредством моделирования двумерным кодом на основе конечно-разностной схемы [9]. В спектрах были обнаружены интервалы резонансных частот, на которых происходит существенное усиление сигнала.

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

fi = Jres 2п \

\rwen J \XZJ

,1 = 1,2,

(6)

где с — скорость звуковой волны в газе, — корни производной функции Бесселя 70(ж), А^ — длина волн, распространяющихся вдоль оси г.

Для случая трехмерного моделирования распределение спектра акустического давления вдоль ствола скважины (называемого далее спектрограммой) представлено на рис. 4. Распределение взято для координаты г = 0.01 м.

Рис. 4. Спектрограмма акустического давления для случая кольцевого источника

Для анализа осесимметричных решений, полученных 2Б- и ЗБ-кодами, рассмотрим спектры в точке г = 0, представленные на рис. 5а. Видно, что резонансы, соответствующие ЗБ-моделированию, смещены вправо.

а)

б)

Рис. 5. Спектры акустического поля в точке г = 0 г = 0.01 м, полученные двумерным и трехмерным моделированием: а) случай описанной окружности, б) случай вписанной окружности в многоугольник, имитирующий скважину на ЗБ-сетке

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

Для оценки эффективного радиуса скважины reff, при котором резонансы совпадают, был проведен расчет с помощью двумерного кода, в котором радиус скважины был задан равным радиусу, вписанной в многоугольник окружности Г[п = rweii cos (^/16) = 0.118 м. Сравнение спектров показано на рис. 56. Из него видно, что сдвиг резонансов практически незаметен и форма спектров схожа между собой за некоторым отличием вблизи резонансов. Таким образом, reff для трехмерного случая можно приближенно оценить значением rweii cos (п/п) где п — число сторон многоугольника, имитирующего скважину в ЗБ-сетке, и использовать его для построения сетки в скважине.

5. Результаты моделирования для случая секториального источника

Рассмотрим конфигурацию с сскториальным источником, представленную на рис. 2. Геометрические параметры модели, физические свойства среды, а также параметры моделирования описаны в разделе 4.1 (при этом reff = 0.118 м). Угловой размер источника составляет 0 ^ в ^ к/8, что соответствует сектору 16-угольника.

Проанализируем спектрограммы акустического давления на оси скважины (рис. 6).

а)

б)

Рис. 6. Спектрограмма акустического давления на оси скважины: а) кольцевой источник, б) сскториальный источник

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

Ниже приведены спектрограммы на удалении от оси скважины, для координаты г = 0.06 м, для случая кольцевого (рис. 7) и секториального (рис. 8) источников. Из них видно, что в случае секториального источника появляется большое количество дополнительных резонансов, вызванных возбуждением азимутальных мод, амплитуда которых зависит от угла в.

Рис. 7. Спектрограмма акустического давления для координаты г = 0.06 м для случая кольцевого источника

Формула (6) для оценки резонансных частот /г'ез легко обобщается для этого случая:

rl,m _ "

J res

2п\

ш2

^weliy \Xz J

+ ( — ) ,l,m _ 1, 2,...,

где — корни производной функции Бесселя J'm (х).

а)

б)

Рис. 8. Спектрограммы акустического давления для координаты г = 0.06 м для случая сектори-ального источника: а) при азимуте в = 0, б) при азимуте в = -ж/2

В табл. 2 приведены значения резонансов (7) для данной конфигурации модели. Из нее видно, что некоторые резонансные моды (например, при т = 0 и т = 2) расположены очень близко друг от друга.

Т а б л и ц а 2

Значения резонансов для рассматриваемой скважины

Порядок функции Бееееля Номер корня производной функции Бееееля

1 = 1 1 = 2 1 = 3 1 = 4 1 = 5 1 = 6 1 = 7

т = 0 2722 Гц 4985 Гц 7228 Гц 9466 Гц 11702 Гц 13 937 Гц 16171 Гц

т = 1 1308 Гц 3788 Гц 6065 Гц 8317 Гц 10 560 Гц 12 800 Гц 15 037 Гц

т = 2 2170 Гц 4765 Гц 7083 Гц 9357 Гц 11615 Гц 13 864 Гц 16108 Гц

т = 3 2985 Гц 5695 Гц 8061 Гц 10 363 Гц 12 639 Гц 14 901 Гц 17155 Гц

т = 4 3778 Гц 6595 Гц 9010 Гц 11342 Гц 13 639 Гц 15 916 Гц 18181 Гц

Уширение интервалов резонансов, обычно наблюдаемое в системе скважина-порода, ведет к тому, что близко расположенные резонанеы накладываются друг на друга в спектре, что видно на рис. 8 в районе 5 кГц.

6. Заключение

Для исследования акустических шумов, порождаемых фильтрационным потоком из пласта в скважину, построена трехмерная численная модель. В ней используется метод спектральных элементов, реализованный в открытом коде SPECFEM3D Cartesian.

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

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

Авторы благодарны Шлюмберже за разрешение опубликовать представленные результаты.

Приложение. Дискретизация модели методом спектральных элементов

Здесь представлен переход от системы (4) к (5).

Для каждого объемного элемента Ое задается отображение (х,у,г) ^ (£> С) в параметрический элемент в форме куба (поверхностные элементы отображаются в квадрат), —1 ^ С ^ 1 —1 ^ V ^ 1 —1 ^ С ^ 1- Предполагается, что отображение является линейным по каждому из направлений, т.е. рассматриваются элементы первого порядка. Для задания отображения используются 8 угловых точек ха, а = 1,..., 8, шестигранника (для поверхностных элементов - 4 угловые точки). Соотношение, связывающее произвольную точку х в Ое с произвольной точкой (£, г/, () в параметрической ячейке, может быть записано в виде

8

х ) = Е ^ (Ьъс) ха,

а= 1

где N0,^, г/, () = 1а(£)1р(С) = 0,1 — произведения линейных интерполяционных

многочленов Лагранжа.

Для решения задачи (4) все интегралы по расчетной области О и по границе расчетной области дО разбиваются па подынтегралы по объемным и поверхностным элементам соответственно. Функции на элементах представляются интерполяционными многочленами Лагранжа Ьа(£) высокого порядка щ, определенными на узлах Гаусса-Лобатто-Лежандра (далее СЬЬ) а = 0,..., щ, которые являются корнями уравнения вида

(1 — £2) Р'а, (0=0,

где Рщ — многочлен Лежандра порядка щ. Произвольная скалярная функция /, а также

компоненты векторной функции, представляются па элементе Ое в виде

щ

/ (х ))= £ /ьа (О ^ (п)ь, (с),

«,,3,7=0

где — значение функции / в узлах СЬЬ.

Интегрирование внутри элемента производится при помощи квадратуры Гаусса-Лобатто-Лежандра:

щ

«,,3,7=0

е

где а = 0— веса интегрирования, ассоциированные с узлами GLL, J"^1 — якобиан отображения элемента Qe в соответствующую параметрическую ячейку, вычисленный в узлах GLL. Поскольку часть узлов может принадлежать нескольким элементам, которые соприкасаются друг с другом по граням, ребрам или вершинам, вклады таких узлов в систему уравнений суммируются с весами по всем элементам, которым этот узел принадлежит, при этом на границе упругой среды с жидкостью учитываются условия (3), см. [8].

Таким образом, после применения вышеописанной процедуры интегрирования ко всем слагаемым в системе (4), в том числе к слагаемым, содержащим взятие градиента, см. [7], задача сводится к решению системы обыкновенных дифференциальных уравнений относительно значений искомых полей во всех точках GLL по всем элементам (объемным и поверхностным).

Литература

1. Ипатов А.И., Кременецкий М.И. Геофизический и гидродинамический контроль разработки месторождений углеводородов. Москва : НИЦ «Регулярная и хаотическая динамика», 2006.

2. Sergeev S.I., Ryzhikov N.I., Mikhailov D.N. Laboratory investigation of sound induced by gas flow in porous media //J. Pet. Sci. Eng. 2019. V. 172. P. 654-661.

3. Мутовкин H.B., Михайлов Д. I/.. Софронов И.Л. Моделирование акустических полей, генерируемых фильтрационным потоком в околоскважинной зоне // Математическое моделирование. 2019. Т. 31, вып. 6. С. 95-107.

4. Mutovkin N.V., Mikhailov D.N., Sofronov I.L. Analysis of Modeling of Acoustic Fields Excited by the Flow Noise in the Formation Inflow Zones // Geomodel 2018 — 20th Conference on Oil and Gas Geological Exploration and Development. 2018. (in Russian).

5. Мутовкин H.B., Михайлов Д. I/.. Софронов И.Л. Оценка изменения фазового состава флюида в скважине на основе данных пассивной акустической шумометрии // Российская нефтегазовая техническая конференция SPE. 2019.

6. Clayton R., Engquist В. Absorbing boundary conditions for acoustic and elastic wave equations 11 Bull. Seismol. Soc. Am. 1977. V. 67, N 6. P. 1529-01540.

7. Komatitsch D., Tromp J. Introduction to the spectral element method for three-dimensional seismic wave propagation // Geophvs. J. Int. 1999. V. 139. P. 806-822.

8. Komatitsch В., Vilotte J.-P., Tromp J. SPECFEM3D Cartesian v2.0.2 [software] [Electronic resource] // Computational Infrastructure for Geodvnamics. 2012. URL: https://geodvnamics.org/cig/software/specfem3d/ (accessed: 04.06.2019).

9. Plyushchenkov B.D., Turchaninov V.I. Acoustic Logging Modeling by Refined Biot's Equations // Int. J. Mod. Phvs. C. 2000. V. 11, N 2. P. 365-397.

References

1. Ipatov A.I., Kremenetskiy M.I. Geophysical and hvdrodvnamic control of hydrocarbon field development. Moscow : NITS «Regular and chaotic dynamics», 2006. (in Russian).

2. Sergeev S.I., Ryzhikov N.I., Mikhailov D.N. Laboratory investigation of sound induced by gas flow in porous media. J. Pet. Sci. Eng. 2019. V. 172. P. 654-661.

3. Mutovkin N.V., Mikhailov D.N., Sofronov I.L. Modeling of acoustic fields, induced by the flow noise in the near-wellbore zone. Matem. Mod. 2019. V. 31, N 6. P. 95-107. (in Russian).

4. Mutovkin N. V., Mikhailov D.N., Sofronov I.L. Analysis of Modeling of Acoustic Fields Excited by the Flow Noise in the Formation Inflow Zones. Geomodel 2018 — 20th Conference on Oil and Gas Geological Exploration and Development. 2018. (in Russian).

5. Mutovkin N. V., Mikhailov D.N., Sofronov I.L. Estimation of Fluid Phase Composition Variation Along the Wellbore by Analyzing Passive Acoustic Logging Data. SPE Russian Petroleum Technology Conference. 2019.

6. Clayton R., Engquist B. Absorbing boundary conditions for acoustic and elastic wave equations. Bull. Seismol. Soc. Am. 1977. V. 67, N 6. P. 1529-1540.

7. Komatitsch D., Tromp J. Introduction to the spectral element method for three-dimensional seismic wave propagation. Geophvs. J. Int. 1999. V. 139. P. 806-822.

8. Komatitsch D., Vilotte J.-P., Tromp J. SPECFEM3D Cartesian v2.0.2 [software] [Electronic resource]. Computational Infrastructure for Geodvnamics. 2012. URL: https://geodvnamics.org/cig/software/specfem3d/ (accessed: 04.06.2019).

9. Plyushchenkov B.D., Turchaninov V.I. Acoustic Logging Modeling by Refined Biot's Equations. Int. J. Mod. Phvs. C. 2000. V. 11, N 2. P. 365-397.

Поступим в редакцию 18.06.2019