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

Методика. Суть метода дифференциальной РСА интерферометрии состоит в обработке пар радарных изображений, снятых в два разных момента времени, между которыми произошла деформация. Разность фаз сигналов, полученных при съемке, отображается с помощью интерферограмм. Монохромная или цветная шкала интерферограммы характеризует величину и направление сдвига, произошедшего при деформации. В нашем исследовании применялись данные с двух спутников: ENVISAT (инструмент ASAR, С-диапазон) и ALOS (инструмент PALSAR, L-диапазон). При длине волны 5.6 см (С-диапазон) фазовый набег 2π на интерферограмме порождается изменением дальности по направлению к спутнику на 2.8 см; при длине волны 23.6 см (L-диапазон) – на 11.8 cм. Разности дальностей в направлении на точку съемки пересчитываются в истинные изменения положения земной поверхности. Преимуществом метода перед другими способами измерения деформаций является высокоточное измерение смещений земной поверхности по всей площади объекта. Существуют ограничения метода, такие, как временная декорреляция, происходящая из-за изменения радиофизических свойств поверхности, и влияние атмосферных и ионосферных явлений. Ограничения и методы их преодоления детально описаны в [Hanssen, 2001].

Геодинамическая обстановка. Район наших исследований находится в пределах Верхнеангарско-Муйской междувпадинной перемычки северо-восточного сегмента Байкальской рифтовой системы, в зоне динамического влияния Муяканского и Перевального разломов (рис. 1). Высокий уровень сейсмичности подтверждает современную геодинамическую активность района исследований. Эпицентральное поле с магнитудами землетрясений 2.6<ML<6.2 наблюдается с 1962 по 2012 гг. (по данным Байкальского филиала Геофизической службы Сибирского отделения Российской академии наук (рис. 2)).

В исследуемом районе трасса БАМ проходит вдоль Муяканской структуры почти на всем ее протяжении и пересекает отдельные ее дислокации. Таким образом, инженерно-сейсмогеологические условия трассы на этом участке неблагоприятны [Solonenko, Mandelbaum, 1985]. Региональный Перевальный разлом простирается в северо-западном направлении. Выраженность и влияние разлома на геодинамическую обстановку наразных его отрезках различна. На исследуемом участке разлом проявлен в поле силы тяжести в виде гравитационной ступени небольшой интенсивности [San’kov et al., 1991].

Интерпретация результатов дифференциальной РСА интерферометрии. Нами произведена интерферометрическая обработка снимков ENVISAT и ALOS/PALSAR. Из-за сильной заболоченности долин и расчлененного рельефа на всех интерферограммах по снимкам ENVISAT велико влияние временной декорреляции. Однако, несмотря на это, на одних и тех же участках на трех парах снимков (13.01.2004–08.03.2005, 16.05.2004–18.09.2005 и 17.02.2004–17.01.2006) обнаружены четыре области высокой когерентности. Для анализа выбрана пара с базой интерферометра, равной 3 м (13.01.2004–08.03.2005) (рис. 3). При такой базе влияние рельефа на интерферометрическую разность фаз исчезающе мало.

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

Линейно локализованные деформации приурочены непосредственно к зоне сместителя разлома (рис. 3), прослеживающегося в основании уступа Муяканского хребта. Рассмотрение размеров аномалий, их протяженности и характера пространственного распределения позволяет исключить экзогенную природу деформаций. Сейсмическая активность зоны разлома в течение 2004–2006 гг. весьма незначительна. Следовательно, мы не можем объяснить генезис этих деформаций накоплением и реализацией упругих напряжений в виде сейсмических событий. Одним из видов деформаций, фиксируемых геодезическими методами в зонах активных разломов, как в подвижных областях, так и на платформах, являются параметрические деформации, впервые выделенные и описанные Ю.О. Кузьминым и называемые им суперинтенсивными [Kuzmin, 2013]. Скорости параметрических деформаций достигают больших величин – до 10–5. Чаще всего, они не сопровождаются сейсмичностью и локализуются непосредственно в зонах разломов. Автор связывает происхождение этих деформаций с изменениями параметров среды внутри зоны разлома, вблизи ее главного сместителя.

Площадные деформации, обнаруженные вблизи окончания Муяканского разлома (рис. 3) на СВ замыкании Муяканской впадины и в западном борту Улан-Макитской впадины, также не могут быть прямо объяснены сейсмогенными движениями. Картина распределения площадных деформаций соответствует распределению напряжений на окончании левостороннего сдвига в упругой среде [Osokina, 2010]. На восточном окончании северного крыла разлома располагается область растяжения и фиксируется относительное опускание земной поверхности, а на юго-восточной части разлома располагается область сжатия, где фиксируется относительное поднятие земной поверхности. Эти данные согласуются с данными геодинамических исследований [San’kov et al., 2000], по которым Муяканский разлом наряду с вертикальной имеет левостороннюю сдвиговую компоненту смещений.

Наряду с деформациями эндогенной природы, в зонах активных разломов Верхнеангарско-Муйской междувпадинной перемычки нами были обнаружены экзогенные деформации. На участке Байкало-Амурской железнодорожной магистрали в нескольких километрах от Северо-Муйского тоннеля (ст. Казанкан, п. Северомуйск, координаты 56.1N 113.8E) с начала 90-х годов прошлого столетия развивается оползневой процесс. Согласно данным нивелирования и геодезических GPS-измерений, проведенных на железной дороге в 2002– 2003 годах «Иркутскжелдорпроектом», железнодорожное полотно сдвигалось в зоне действия оползня, вдоль склона со скоростью до 2 см/мес. Развитие оползня приводит к риску разрушения железнодорожного пути и крушения поездов. Происхождение оползня связывается исследователями с процессом деградации многолетнемерзлых пород и подрезанием склона дорогой, а также чрезмерной нагрузкой на склон [Trzhtsinsky et al., 2004]. Нельзя также недооценивать влияние уровня сейсмической и тектонической активности района на разрушение пород и развитие оползней, поскольку Казанканский оползень находится в зоне влияния сейсмически активных Муяканского и Перевального разломов.

Для изучения оползневого процесса в районе ст. Казанкан мы использовали снимки с восходящего витка ALOS/PALSAR. Несмотря на зимнее время съемки и разницу между съемками в два года, хорошую корреляцию и результаты показала пара радарных снимков ALOS/PALSAR, снятых 17.01.2009 и 12.01.2007. Местоположения деформаций, наблюдаемых на интерферограммах, точно совпадают с расположением мостов (рис. 4, 5, 6), где по данным предыдущих исследований [Trzhtsinsky et al., 2004] обнаружены значительные подвижки. У второго и третьего моста, ниже железнодорожного полотна, обнаруживаются области, испытывающие воздымание, что, по-видимому, связано с характером развития оползня. Модель развития оползня по геологическим данным описывается как выдавливание по цилиндрической плоскости скольжения с отрицательным знаком смещения в области цирка и положительным – во фронтальной части [Trzhtsinsky et al., 2004].

На исследуемой интерферограмме наибольшее изменение дальности наблюдается ниже железнодорожных путей у второго моста. За двухлетний период, прошедший между измерениями, изменение дальности по направлению от поверхности к радару составило 6 см. Кроме того, вблизи мостов наблюдаются области с опусканием земной поверхности, связанные, по-видимому, с процессами солифлюкции. По данным РСА интерферометрии уточнено место положения оползневого цирка между вторым и третьим мостом (рис. 6). Ранее, вследствие недостатка данных, положение этого цирка отмечалось ниже по склону. Выявлено распространение процесса оползания, происходящее на склоне между первым и вторым мостом, которое не фиксировалось ранее. В будущем возможно объединение цирков и, как следствие, усиление оползневого процесса. Из полученных данных можно сделать вывод о том, что проведенные мероприятия по укреплению склона и железнодорожной насыпи не привели к остановке оползневого процесса, что говорит об острой необходимости продолжения укрепительных работ и принятия мер по изменению трассы на данном участке магистрали.

Заключение. Таким образом, исследование деформаций в районе Верхнеангарско-Муйской междувпадинной перемычки северо-восточного сегмента Байкальской рифтовой системы с применением метода дифференциальной РСА интерферометрии позволило получить новые данные о характере современных смещений земной поверхности различного происхождения. В зоне динамического влияния Муяканского разлома обнаружены два типа деформаций, которые можно определить как тектонические линейно локализованные и площадные. Наблюдаемые деформации, находящиеся непосредственно в зоне разлома, можно отнести к параметрическим, связанным с изменением состояния среды в зоне разлома. Картина распределения площадных деформаций соответствует распределению напряжений на окончании левостороннего сдвига в упругой среде.

Первые РСА-исследования, проведенные в районе п. Северомуйск, показали, что снимки ALOS/PALSAR применимы для детального изучения оползневого процесса, происходящего в районе ст. Казанкан БАМ ВСЖД. Оползневой процесс продолжает свое развитие, несмотря на проведенные укрепления. В настоящее время происходит объединение цирков на участке между первым и вторым мостом, что может привести к усилению процесса.

Таким образом, деформации в зоне Муяканского разлома могут представлять опасность для линейных сооружений Байкало-Амурской магистрали. Метод РСА интерферометрии имеет большой потенциал как для изучения уже произошедших событий в районе исследований, так и для мониторинга развивающихся процессов. В настоящее время РСА ALOS/PALSAR закончил свою работу на орбите, однако с 2014 года Японским аэрокосмическим агентством JAXA запущены новые усовершенствованные спутники РСА миссии ALOS2/PALSAR2, работающие на той же длине волны L-диапазона. Мы планируем продолжить наши исследования с использованием новых данных.