Представлен численный метод решения двумерной обратной динамической задачи сейсмики для вязкоупругой изотропной среды. В качестве математической модели рассматривается система дифференциальных уравнений упругости для изотропных сред с памятью. Искомыми величинами являются смещение точек поверхности, функция памяти среды (ядро интегрального члена) и скорость распространения упругих поперечных волн в слабо горизонтально-неоднородной среде при воздействии на границу полупространства направленной мгновенной силы. Дополнительной информацией для решения обратной задачи является отклик смещения, измеренный на дневной поверхности. Метод основан на сведении обратной задачи к системе интегральных уравнений типа Вольтерра и их последовательной численной реализации. Приводится анализ результатов исследования и сравнение с аналитическим решением. Показано, что результаты находятся в удовлетворительном соответствии.