Источник: http://tt.pstu.ru/mnp10/mnp10/s2/liskov.htm
Для выявления и исследования оседаний земной поверхности используется радарная интерферометрия. Этот метод хорошо подходит для решения задач по определению небольших относительных смещений земной поверхности с применением радиоволн с короткой длиной волны (сантиметрового порядка). В настоящее время для выявления оседаний подработанных территорий используется техника интерферометрического анализа спутниковых съемок. В данной статье представлен анализ оседаний территории г. Березники методом площадной интерферометрии.
Ключевые слова: мониторинг, радар, радарная интерферометрия, деформации земной поверхности.
Интенсивная разработка рудных и нефтегазовых месторождений, промышленно-гражданская застройка территорий, наличие многочисленных оползнеопасных мест, наличие разломных тектонических структур и ряда отдельных геодинамических и геоморфологических факторов могут оказать негативное влияние на состояние горного массива и земной поверхности. Антропогенное влияние на недра и земную поверхность может сопровождаться как обычными проседаниями земной поверхности, так и аномальными деформациями и проявлениями сейсмической активности, что может вызвать значительный экономический и экологический ущерб, особенно в районах промышленно-гражданской застроенной территории. В этой связи, возникает необходимость проведения мониторинга за состоянием поверхности и горного массива. На сегодняшний день существуют несколько методов мониторинга для точного и достоверного выявления оседаний земной поверхности.
Традиционно для выполнения такого рода работ используется высокоточное нивелирование и GPS измерения. Данные технологии геодезических измерений позволяет проводить наблюдения на миллиметровом уровне точности. В последние время все более широкое применение находит метод дифференциальной интерферометрии SAR, который также позволяет выявлять проседания также высоком - миллиметровом уровне точности. Технология дифференциальной интерферометрии SAR использует методы радиолокационной съемки земной поверхности с космических спутников, на которых установлены антенны с синтезированной апертурой (SAR). Основной информацией, получаемой в результате радарной съемки, являются интенсивность и фаза (временная задержка сигнала). Повторная радарная съемка позволяет определить разность фаз, обусловленную, например, сдвижением земной поверхности. Такие смещения могут быть определены в результате обработки радарной съемки с высокой точностью.
Преимущества такого способа заключается в том, что можно вести наблюдение за всей площадью исследуемого участка (площадью до 3000 кв. км), а не за конкретными объектами или линиями. Экономически данный метод тоже эффективен, так как исключается непосредственное участие человека при выполнении дорогостоящих геодезических полевых работ. Получение радарных данных с космических спутников не занимает много времени (максимум 6 дней), что в свою очередь позволяет получать картину деформаций гораздо оперативнее. Так как спутник вращается по повторяющимся орбитам, период наблюдений может колебаться от 11 до 45 дней, в зависимости от типа спутника.
В рамках Национального Образовательного Проекта (НОК), выигранного Пермским государственным техническим университетом, был приобретен современный программный комплекс «GAMMA» (фирма «GAMMA REMOTE SENSING AG», Швейцария) для обработки данных радарной космической съемки. Специалисты ПГТУ прошли обучение и в настоящее время полностью освоили работу с данным комплексом. В качестве исследовательских проектов были выполнены обработки спутниковых данных по территории г.Березники, а также на Астраханском и Оренбургском геодинамических полигонах. Полученные результаты свидетельствуют о полной работоспособности метода.
Верхнекамское месторождение калийных солей. Город Березники.
В рамках научно-исследовательского проекта Пермский государственный технический университет выполнил исследования по определению вертикальных смещений земной поверхности г.Березники дистанционным методом радарной интерферометрии. Как известно в 2007 году в результате прорыва воды был затоплен и первый рудник БКПРУ-1. В связи с разработкой и последующим затоплением первого рудника в некоторых частях города зафиксированы серьезные деформации земной поверхности, а, соответственно, и промышленных и жилых зданий. В связи с этим мониторинг развивающихся деформаций территории города стал чрезвычайно актуальной задачей.
В этом проекте были использованы данные со спутника ENVISAT, предоставленные Европейским космическим агентством ESA. Данный спутник имеет радар с синтезированной апертурой (SAR). Средние значения деформаций земной поверхности и их история развития являлись основным интерферометрическим результатом.
Для исследования была использована серия повторных наблюдений со спутника ENVISAT. Общее количество сцен составило 28 штук. Они охватывали период времени с начала 2004 по конец 2008 года. В качестве дополнительной информации для подтверждения полученных интерферометрических результатов использовались данные нивелирования, предоставленные ОАО «Уралкалий», а также данные GPS измерений на Верхнекамском глобальном геодинамическом полигоне, созданном Пермским государственным техническим университетом.
В результате сложнейшей обработки и интерпретации радарных снимков была получена деформационная картина на территорию города Березники. Для сравнительного анализа полученных интерферометрических данных был выбран небольшой участок, который охватывал центральную часть города, где были зафиксированы значительные смещения земной поверхности. На данном участке выполняются ежегодные инструментальные наблюдения за осадками силами маркшейдерской службы ОАО «Уралкалий». Высокоточное нивелирование выполняется по профильным линиям, которые представляют собой сеть грунтовых реперов и стенных марок. Всего для анализа было выбрано четыре профильные линии, лежащие в зоне деформаций. Помимо этого, на каждой профильной линии было отмечено по одному реперу для изучения их истории развития смещений. Положение профильных линий и выбранных реперов показано на рис. 1. Здесь же показан фрагмент карты деформаций земной поверхности г. Березники.
Как видно на рис. 1 на участке сформировались две мульды оседаний, которые из-за близкого расположения слились в одну общую зону деформаций. За четыре года наблюдений максимальные оседания на этой территории достигли значений в 20-25 см.
Для сравнения данных нивелирования и результатов интерферометрического анализа по каждой профильной линии были построены графики вертикальных смещений (на рис. 2 приведен один график по профильной линии XXVII). Рассматриваемый период — четыре года: с июня 2004 по июль 2008 года. И как видно на графике, интерферометрические результаты практически полностью совпадают с инструментальными.
С помощью точечного интерферометрического анализа была также получена информация о том, как развивались деформационные процессы в каждой точке за наблюдаемый период времени. На рисунке 2а приведены данные о развитии смещений для одного из реперов (их местоположение представлено на рис. 1). Анализируя этот график (рис. 2б), можно сделать вывод, что смещения реперов имеют линейную закономерность с незначительными отклонениями, а, соответственно, и весь участок с 2004 по конец 2008 года оседал равномерно без изменения скоростей деформаций, даже после полного затопления в 2007 году рудника. Незначительные расхождения интерферометрических результатов с данными нивелирования могут быть связаны с остаточным влиянием атмосферы или с ошибками развертывания фазы. В целом же, полученные по результатам дистанционного зондирования данные о смещениях, очень хорошо согласуются с результатами инструментальных наблюдений.
Полученная информация об оседаниях земной поверхности позволила восстановить историю развития деформаций не только в конкретной точке, но и всего участка. Так на рис. 3 приведена пространственная информация о развитии деформаций на всем участке за период с мая 2007 по май 2008 года.
Интерферометрический анализ радарных сцен на территорию города Березники дал очень интересные результаты об активности деформаций, связанных с разработкой первого рудника и дальнейшим его затоплением, и их пространственном распределении. Высокая корреляция результатов анализа с данными инструментальных наблюдений говорит о потенциальной возможности использования метода радарной интерферометрии для контроля смещений земной поверхности, или хотя бы для прогнозирования деформаций и выявление потенциально опасных зон, в которых, в дальнейшем, необходимо заложить наблюдательные станции. Применение точечного анализа позволило получить информацию не только о пространственном распределении оседаний в городе, но и во времени. Анализ развития смещений во времени позволил судить о скоростях изменения деформаций на участке.
Представленная в статье техника особенно хорошо себя показала в застроенной части города, чего не скажешь о периферии, где степень декорреляции очень велика из-за наличия лесной растительности. Малое количество точек-отражателей на участке промзоны первого рудника связано еще и с относительно высокими смещениями поверхности, возникшими после затопления рудника. Использование радаров с меньшим интервалом повторных съемок, например, TERRASAR-X, COSMO-SKYMED позволит осуществлять наблюдение даже за такими участками, где скорость оседаний превышает 50-70 см/год.
Таким образом, интерферометрический анализ радарных сцен техногенно нагруженных территорий, а также территорий с высокой плотностью застройки промышленными и гражданскими объектами дает очень важные и практические результаты об активности деформаций и их пространственном распределении. Применение точечного анализа позволяет получить информацию не только о пространственном распределении деформаций, но и во времени. Открывается возможность оперативного мониторинга деформационных процессов больших территорий. После выявления на основе спутниковой интерферометрии наиболее опасных мест для предупреждения возможных аварий можно выполнять локальный геодезический мониторинг данных участков традиционными геодезическими методами.
1. Cumming, I. G., and Wong, F. H., Digital processing of synthetic aperture radar data. Norwood, MA: Artech House, Inc., 2005.
2. Curlander, J. C., and McDonough, R. N., Synthetic Aperture Radar: Systems and Signal Processing. New York: John Wiley & Sons, Inc., 1991.
3. Elachi, C., Spaceborne Radar Remote Sensing: Applications and Techniques. New York: IEEE Press, 255, 1987. Franceschetti, G., and Schirinzi, G., "A SAR Processor Based on Two-Dimensional FFT Codes," IEEE Transactions on Aerospace and Electronic Systems, vol. 26, 2, pp. 356-365, 1990.
4. ESA TM-19, InSAR Principles: Guidelines for SAR Interferometry Processing and Interpretation. Noordwijk: ESTEC, 2007.
5. Raney, R. K., "Radar Fundamentals: Technical Perspective," in Principles and Applications of Imaging Radar, vol. 2, Manual of Remote Sensing, Henderson, F. M., and Lewis, A. J., Eds., 3rd ed. New York: John Wiley & Sons, pp. 9-130, 1998.
6. Rees W.G., Physical principles of Remote sensing. Cambridge University Press, 2001.
7. Бакулев П.А., Радиолокационные системы. — М: Радиотехника, 2004.
8. Гужов В.И. Ильиных С.П., Компьютерная интерферометрия: Учеб. пособие. — Новосибирск: Изд-во НГТУ, 2004.
9. Кашников Ю.А., Ашихмин С.Г., Механика горных пород при разработке месторождений углеводородного сырья. — М.: ООО "Недра-Бизнесцентр", 2007.