В настоящем исследовании изучалось влияние пространственного разрешения цифровой модели рельефа (ЦМР) на достоверность оценки оползневой опасности в Шапе провинции Лаокай. Для решения поставленной задачи был использован модифицированный метод анализа иерархий (MAHP) при двух разных разрешениях ЦМР - 30 м (MAHP30) и 10 м (MAHP10). Обе модели учитывали следующие факторы оползнеобразования: 1) крутизна склонов; 2) состав грунтов; 3) расстояние от активных разломов; 4) вертикальное расчленение рельефа; 5) горизонтальное расчленение рельефа; 6) обводненность массивов; 7) количество осадков; 8) тип растительности. В результате район исследования был разделен на пять зон по оползневой опасности: очень низкая, низкая, средняя, высокая и очень высокая. Достоверность карт определялась с использованием обучающей выборки (91 оползень - обучающая выборка и 39 оползней – данные для проверки). Выполненные исследования показали, что модель MAHP10 является более достоверной в сравнении с моделью MAHP30. Таким образом, качество оценки оползневой опасности имеет прямую корреляцию с пространственным разрешением цифровой модели рельефа.
Качество имеющихся исходных данных играет важную роль в достоверности результатов моделирования оползневой опасности. Одним из наиболее значимых параметров для выполнения анализа является информация о рельефе местности [16].
Цифровая модель рельефа (ЦМР) может быть создана с использованием различных методов, включая наземную съемку, аэрофотограмметрию, LiDAR (англ. Light Detection and Ranging) и InSAR (англ. Interferometric Synthetic Aperture Radar) съемки. Качество ЦМР может быть оценено на основе ее пространственного разрешения. Пространственное разрешение – величина, определяющая количество точек (элементов растрового изображения) на единицу площади (или единицу длины). Разрешение растровой структуры данных может быть охарактеризовано в терминах размера ячеек сетки. Выбор оптимального размера зависит от целей исследований и масштаба анализа [6].
Лаокай – одна из северных провинций Вьетнама, наиболее подверженная оползневым процессам [15]. Расположенный на ее северо-западе район Шапа (Рисунок 2.1) находится примерно в 38 км от города Лаокай и в 376 км от Ханоя. Он покрывает площадь 675,8 км2 и имеет высоту над уровнем моря от 150 м до более чем 3000 м. Район Шапа столкнулся с проблемой активизации оползневых и эрозионных процессов вследствие сельскохозяйственного освоения и развития туризма [7].
Рис.1. Расположение района исследования
Оценка оползневой опасности является важнейшим компонентом национальной стратегии предотвращения и уменьшения последствий стихийных бедствий. Главной задачей при проведении анализа оползневой опасности является определение условий оползнеобразования, под которыми следует понимать всю совокупность природных и антропогенных факторов, нарушающих равновесие масс горных пород [3]. Для каждой территории набор таких условий может существенно отличаться. Разработка карты оползневой опасности зависит от сложной взаимосвязи между проявлениями оползней и оползневыми факторами. Переменные, которые способствовали составлению карт, были выбраны на основе обзора литературы и подробных полевых исследований [8]. Обе модели учитывали следующие факторы оползнеобразования: 1) крутизна склонов; 2) состав грунтов; 3) расстояние от активных разломов; 4) вертикальное расчленение рельефа; 5) горизонтальное расчленение рельефа; 6) обводненность массивов; 7) количество осадков; 8) тип растительности.
Важнейшим элементом, необходимым для оценки достоверности выполненной оценки оползневой опасности является информация о фактическом распределении существующих проявлений оползневых процессов, которая отображается на «карте инвентаризации оползней» [9]. При составлении карт инвентаризации (LIM) внимание фокусируется на определении границ оползней и игнорирует особенности оползневых деформаций. Наряду с отображением пространственного распределения оползней карта инвентаризации может включать другие виды информации, такие как: геометрические характеристики оползня (размеры, площадь, глубина захвата склонового массива оползневыми деформациями), особенности геологического строения (литология, структура, свойства грунтов) и гидрогеологических условий [14].
В районе исследования было выявлено в общей сложности 130 оползней, при этом наименьшая площадь оползня составила 395,6 м2, а самая большая площадь оползня – 70447,4 м2. Большинство оползней было зафиксировано в коммуне Чунгчай, в центре района Шапа, в коммунах Лаочай, Таван и Банхо (Рис. 1).
Метод анализа иерархий (англ. Analytic Hierarchy Process – AHP), предложенный Т. Л. Саати [5, 12], предоставляет экспертам относительно простой и эффективный способ измерения объективных и субъективных факторов посредством попарных относительных сравнений и вычисления соответствующих приоритетов шкалы отношений. Он является примером полуколичественного подхода – многокритериального анализа [13].
Ранее выполненные исследования [4, 10, 11] показали, что основным недостатком классической формулировки AHP является необходимость использования экспертных оценок, что вносит определенный субъективизм в результат оценки региональной оползневой опасности.
В настоящем исследовании был использован
модифицированный метод анализа иерархий (MAHP) [1, 2, 4] при двух разных разрешениях ЦМР — 30 м (MAHP30) и 10 м (
MAHP10). В исследовании Нгуена Ч. К. [3], модель MAHP30 использовалась для оценки оползневой опасности в северо-западном регионе провинции Лаокай, который включал район Батсат, район Шапа и город Лаокай. Модель MAHP10 была использована авторами для оценки оползневой опасности в районе Шапа [1]. Основным критерием для определения значимости фактора в MAHP является форма функции распределения выявленных оползней по информационным классам рассматриваемого фактора [2].
В результате были получены карты, на которых район исследования был разделен на пять зон по оползневой опасности: очень низкая, низкая, средняя, высокая и очень высокая.
Достоверность карт определялась с использованием обучающей выборки (91 оползень — обучающая выборка и 39 оползней – данные для проверки).
Зоны очень высокой и очень низкой оползневой опасности – две зоны, которые вызывают наибольший интерес для анализа, поскольку категории неустойчивости/устойчивости для них выражены наиболее контрастно. Уровень достоверности, полученный для этих двух зон имеет определяющее значение для оценки точности карт оползневой восприимчивости и опасности (Рис. 2 и Рис. 3).
Рис. 2. Карта оползневой восприимчивости (a), карта оползневой опасности (б) с применением модели MAHP30
Рис. 3. Карта оползневой восприимчивости (a), карта оползневой опасности (б) с применением модели MAHP10
Распределение зон оползневой опасности по площади приведено на рис. 4 и в таблице 1. Площадь зоны очень низкой оползневой опасности составляет 30,26% на карте, построенной на основе ЦМР, с разрешением 30м. (MAHP30) и 24,86%, на карте, построенной на основе ЦМР, с разрешением 10м. (MAHP10). Зона очень высокой оползневой опасности занимает 12,3% (MAHP 30) и 14,73% (MAHP10) соответственно.
Рис. 4. Распределение зон оползневой опасности по площади с использованием моделей MAHP30 и MAHP10
Таблица . Распределения оползней по классам оползневой опасности
|
Модель |
|
Класс |
Очень низкая |
Низкая |
Средняя |
Высокая |
Очень высокая |
|
MAHP30 |
% площадь класса |
30,26 |
24,92 |
17,46 |
15,06 |
12,30 |
|
|
% площадь оползня |
Обучение |
32,14 |
17,53 |
10,23 |
13,94 |
26,16 |
|
|
Проверка |
8,60 |
37,66 |
25,26 |
17,21 |
11,27 |
||
|
MAHP10 |
% площадь класса |
24,86 |
20,29 |
19,11 |
21,01 |
14,73 |
|
|
% площадь оползня |
Обучение |
6,38 |
19,42 |
24,46 |
21,15 |
28,59 |
|
|
Проверка |
3,76 |
13,63 |
19,80 |
41,85 |
20,96 |
||
Чтобы оценить, какая модель является более достоверной, было изучено распределение оползней в пределах, выделенных по степени опасности зон с использованием результатов по обучающей выборке и данных проверки (Рис. 5 и Рис. 6).
Рис. 5. Сравнение распределения процента площади оползней с использованием обучающей выборки
Рис. 6. Сравнение распределения процента площади оползней с использованием данных проверки
AUC ROC (площадь под кривой ошибок) анализ так же показал, что при использовании ЦМР более высокого разрешения, достоверность модели увеличивается (Рис. 7) и, таким образом, можно сделать вывод, что качество оценки оползневой опасности имеет прямую корреляцию с пространственным разрешением цифровой модели рельефа.
Рис. 7. Кривая успеха (a) и прогнозирования (b) моделей MAHP30 и MAHP10 по кривой ROC
Заключение
Выполненный анализ показал, что модель MAHP10 является более достоверной (площадь оползней в зоне очень высокой опасности, полученная оползня с использованием обучающей выборки – 28,59%, площадь оползней, полученная в зоне очень высокой опасности с использованием данных проверки — 20,96%), в сравнении с моделью MAHP30 (площадь оползней в зоне очень высокой опасности, полученная оползня с использованием обучающей выборки – 26,16%, площадь оползней, полученная в зоне очень высокой опасности с использованием данных проверки – 11,27%). Аналогичные результаты получены и для зоны с очень низкой опасностью. В модели MAHP10 ее площадь составила 6,38% (для обучающей выборки) и 3,76% (для данных проверки), в модели MAHP10 – 32,14% и 8,6% соответственно.