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