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