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