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