АЛГОРИТМ РЕШЕНИЯ ПРЯМОЙ КИНЕМАТИЧЕСКОЙ ЗАДАЧИ СЕЙСМИКИ В ТРЕХМЕРНЫХ НЕОДНОРОДНЫХ ИЗОТРОПНЫХ СРЕДАХ
Решение прямой кинематической задачи сейсмики - важный этап обработки данных полевых измерений в сейсмологии и сейсморазведке, например, при решении обратной кинематической задачи сейсмики методом нелинейной лучевой сейсмотомографии [1], в методе пространственно-временной миграции Кирхгофа и др. К настоящему моменту существует несколько подходов к решению данной задачи. К наиболее популярным из них относятся: 1) методы, основанные на решении двухточечной краевой задачи; 2) метод волновых фронтов; 3) разностные методы решения уравнения эйконала и др. [2, 3]. Предлагаемый в работе метод относится к методам лучевого трассирования [4] и преодолевает некоторые ограничения, присущие вышеперечисленным методам. Основой подхода, предложенного в данной работе, является итерационный метод Ньютона - Канторовича. Алгоритм трехмерной пристрелки прост в реализации и допускает «почти» линейное распараллеливание. Вопросы сходимости метода Ньютона в применении к поставленной задаче в данной работе не рассматриваются. Также не рассматриваются вопросы применимости метода при наличии каустик, зон тени, дифракции.