В статье предложен метод расчета трехмерных сейсмических волновых полей для анизотропной упругости. Авторами разработаны устойчивые конечно-разностные схемы до 8-го порядка точности внутри области и до 4-го на границе. Для расчета сложных геологических моделей с топографией используются адаптивные криволинейные сетки. Сравнение численных и аналитических решений задачи Лэмба демонстрирует повышение точности моделирования и значительную экономию компьютерной памяти при повышении порядка аппроксимации.
