Российский фонд
фундаментальных
исследований

Физический факультет
МГУ им. М.В.Ломоносова
 

Вычислительные методы и программирование. 2014. 15, № 3

 

Неклюдов Д.А., Сильвестров И.Ю., Чеверда В.А. «Итерационный метод решения трехмерного уравнения Гельмгольца с "почти аналитическим" предобусловливателем для моделирования акустических волновых полей в задачах сейсморазведки» Вычислительные методы и программирование, 15, № 3, с. 514-529 (2014)

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

Вычислительные методы и программирование, 15, № 3, с. 514-529 (2014) | Рубрики: 04.12 12.06