Đang tải... (xem toàn văn)
Tạp chí Khí tượng Thủy văn 2022, EME4, 85 93; doi 10 36335/VNJHM 2022(EME4) 85 93 http //tapchikttv vn/ Bài báo khoa học Ứng dụng phân tích thứ bậc (AHP) và GIS đánh giá nguy cơ trượt lở đất ở huyện L[.]
Bài báo khoa học Ứng dụng phân tích thứ bậc (AHP) GIS đánh giá nguy trượt lở đất huyện Lộc Ninh, tỉnh Bình Phước Nguyễn Huy Anh1*, Nguyễn Trung Kiên1, Trần Mỹ Hảo1 Đại học Tài ngun Mơi trường TP Hồ Chí Minh; huyanhunre@gmail.com; tranhao_tnmt@yahoo.com *Tác giả liên hệ: anhnh@hcmunre.edu.com; Tel.: +84–905210473 Ban Biên tập nhận bài: 5/2/2022; Ngày phản biện xong: 16/3/2022; Ngày đăng bài: 25/4/2022 Tóm tắt: Nghiên cứu sử dụng kết hợp phương pháp phân tích thứ bậc Saaty công nghệ GIS để xây dựng đồ nguy trượt lở đất địa bàn huyện Lộc Ninh, tỉnh Bình Phước Kết nghiên cứu chia nguy trượt lở đất thành cấp bao gồm: vùng có nguy trượt lở thấp diện tích 24.718,81ha chiếm 28,95% diện tích tự nhiên, vùng có nguy thấp 28.889,21ha chiếm 33,83%, vùng nguy trượt lở trung bình 5.751,43 chiếm 6,74%, vùng có nguy trượt lở cao có diện tích 20.660,45 chiếm 24,19% cao 5.375,09ha chiếm 6.29% Từ khóa: Trượt lở đất; AHP; GIS; Viễn thám; Huyện Lộc Ninh Đặt vấn đề Trượt lở đất (TLĐ) loại hình thiên tai phổ biến giới Việt Nam Ba phần lãnh thổ Việt Nam thuộc khu vực miền núi, có địa hình sườn dốc cao, hoạt động phát triển kinh tế–xã hội (KT–XH) chưa qui hoạch hợp lý, nên tượng TLĐ, lũ bùn đá lũ quét thường xảy Những năm gần đây, loại hình thiên tai xảy với tần suất cường độ ngày tăng, gây thiệt hại nghiêm trọng, điển hình tỉnh miền núi phía Bắc, dun hải miền Trung Tây Nguyên [1] Trượt lở đất dạng tai biến vừa có tính chất tiềm ẩn vừa có tính chất hiểm họa, gây nhiều thiệt hại lớn, khó khăn bất lợi cho sống sinh hoạt người dân [2] Tai biến trượt lở khiến mơi trường cảnh quan bị hủy hoại Do đó, việc nghiên cứu thành lập đồ nguy tai biến trượt lở nhu cầu cần thiết Lộc Ninh, huyện vùng sâu thuộc tỉnh Bình Phước có địa hình đồi núi chiếm 1/3 diện tích khu vực, phân bố hầu hết mở rộng phía đơng huyện Lộc Ninh, nơi cao có độ cao 235 m [3] Trong năm gần diện tích rừng ngày suy giảm, nạn khai thác đất trái phép gia tăng, vào mua mưa nguy xảy trượt lở đất vùng đồi núi cao Nhằm giảm thiểu thiệt hại tai biến trượt lở đất gây cần phải có nghiên cứu đánh giá xác lập khu vực có nguy trượt lở khác Từ sở cho quy hoạch phịng tránh nguy trượt lở đất địa bàn huyện Lộc Ninh Trên giới nghiên cứu trượt lở đất thực từ lâu, Việt Nam năm gần có nhiều cơng trình nghiên cứu đánh giá tai biến trượt lở đất Ngoài nghiên cứu theo hướng truyền thống, hướng ứng dụng Viễn thám (RS) công nghệ GIS, RS nghiên cứu trượt lở quan tâm nhiều với nghiên cứu [1–2, 4–7]…, nghiên cứu sử dung phương pháp phân tích thứ bậc Saaty (AHP) để tính tốn trọng số cho yếu tố thể vai trò chúng nguy trượt lở đất địa bàn nghiên cứu GIS RS sử dụng cơng cụ quan trọng để tính tốn, tích hợp yếu tố thành phần để xây dựng đồ nguy trượt lở đất Tạp chí Khí tượng Thủy văn 2022, EME4, 85-93; doi:10.36335/VNJHM.2022(EME4).85-93 http://tapchikttv.vn/ Tạp chí Khí tượng Thủy văn 2022, EME4, 85-93; doi:10.36335/VNJHM.2022(EME4).85-93 86 Tại địa bàn huyện Lộc Ninh chưa có cơng trình nghiên cứu liên quan đến tai biến trượt lở đât Mặc dù chưa có hậu lớn gây TLĐ địa bàn nghiên cứu, nhiên việc xây dựng đồ phân vùng nguy xảy TLĐ cần thiết, hỗ trợ quyền địa phương công tác quy hoạch quản lý tài ngun, mơi trường phịng tránh thiên tai có TLĐ, đặt biệt vùng đồi núi Lãnh thổ phương pháp nghiên cứu 2.1 Khái quát lãnh thổ Lộc Ninh có diện tích tự nhiên 853,95 km² với địa hình đồi chủ yếu, độ cao trung bình từ 70 m đến 220 m Địa bàn huyện nằm vùng khí hậu nhiệt đới gió mùa cận xích đạo, gió bão có nhiệt cao, nhiệt độ trung bình khoảng 26–27°C, giống vùng khác khu vực Đơng Nam Lượng mưa bình quân hàng năm biến động từ 2.045–2.325 mm Mùa mưa diễn từ tháng 5–11, chiếm 85–90% tổng lượng mưa năm, tháng có lượng mưa lớn 376 mm (tháng 7) Mùa khô từ cuối tháng 11 đến đầu tháng năm sau, lượng mưa chiếm 10–15% tổng lượng mưa năm, tháng có lượng mưa tháng 2, Nằm vùng có lượng xạ dồi dào, tổng số nắng năm từ 2.400 – 2.500 Số nắng trung bình ngày từ 6,2 – 6,6 Thời gian nắng nhiều vào tháng 1, 2, 3, thời gian nắng vào tháng 7, 8, [8] 2.2 Dữ liệu nghiên cứu Các liệu sử dụng cho nghiên cứu đánh giá nguy trượt lở đất huyện Lộc Ninh gồm có: (1) Các tài liệu điều kiện tự nhiên, kinh tế–xã hội huyện Lộc Ninh tỉnh Bình Phước năm 2020, thu thập từ UBND huyện phịng Tài ngun mơi trường huyện Lộc Ninh; (2) Bản đồ địa hình file số với tỷ lệ 1:50.000 thu thập sở TNMT tỉnh Bình Phước, đồ địa hình dạng file số sử dụng giai đoạn xây dựng đồ phân bố mạng lưới thủy văn, tính tốn mật độ chia cắt ngang địa hình, kết hợp với mơ hình số độ cao (DEM) tính tốn độ dốc, độ chia cắt sâu, hướng sườn; (3) Mơ hình số độ cao (DEM) độ phân giải 12,5 m thu thập từ website https://search.asf.alaska.edu [9]; (4) Ảnh vệ tinh Landsat độ phân giải 30 m thu thập ngày 28/01/2020 sử dụng để xây dựng liệu lớp phủ thực vật khu vực nghiên cứu, thông số ảnh Landsat thể bảng 1, kênh 3, 4, sử dụng tính tốn số thực vật (NDVI) [3]; (4) Bản đồ đất, trạng sử dụng đất, liệu mưa thu thập từ Sở Tài Nguyên Môi Trường tỉnh Bình Phước Bảng Thơng số ảnh vệ tinh landsat Vệ tinh LDCM Landsat8 cảm OLI TRIs Kênh/Band Bước sóng (micrometers) Độ phân giải (meters) Kênh – Costal aerosol (xanh tím) – quan sát vùng nước ven bờ hạt mịn Kênh – Blue (xanh lơ) Kênh – Green (xanh lục) Kênh – Red (Đỏ) Kênh – Near Infrered (NIR) – (Cận hồng ngoại) Kênh – SWIR1 (Hồng ngoại sóng ngắn 1) Kênh – SWIR2 (Hồng ngoại sóng ngắn 2) Kênh – Panochromatic (Tồn sắc) Kênh – Cirrus (Phát mật độ, độ dày mây ti) Kênh 10 – Thermal Infrered (NIR) (Hồng ngoại nhiệt 1) 0,433–0,453 30 0,450–0,515 0,525–0,600 0,630–0,680 0,845–0,885 1,560–1,660 2,100–2,300 0,500–0,600 1,360–1,390 10,3–11,3 30 30 30 30 30 30 15 30 30 Kênh 11 – Thermal Infrered (NIR) (Hồng ngoại nhiệt 2) 11,5–12,5 30 Hội nghị khoa học toàn quốc “Chuyển đổi số công nghệ số Khoa học Trái đất, Mỏ Môi trường” (EME 2021) Tạp chí Khí tượng Thủy văn 2022, EME4, 85-93; doi:10.36335/VNJHM.2022(EME4).85-93 87 2.3 Phương pháp nghiên cứu Phương pháp khảo sát thực địa: Trong nghiên cứu này, nhóm tác giả triển khai nhiều đợt khảo sát Trong đợt khảo sát địa bàn nghiên cứu thực nhiệm vụ định vị điểm trượt, vết trượt trạng, khứ ghi chép đặc điểm tính chất khối trượt, điểm trượt, vết trượt đồng thời tiến hành đo vẽ quy mô điểm, khối trượt để hỗ trợ tính tốn thể tích, khối lượng đất đá trượt Phương pháp sử dụng ảnh viễn thám: Bài báo sử dụng nguồn ảnh thu nhận từ vệ tinh Landsat năm 2020 với độ phân giải 30 m Ảnh vệ tinh sử dụng nghiên cứu để xác định điểm trượt khứ, xây dựng đồ trạng lớp phủ thảm thực vật theo số thực vật (NDVI) Bằng cách sử dụng ảnh Landsat tải miễn phí từ https://earthexplorer.usgs.gov, ảnh thu nhận ngày 28/01/2020 (công thức 1) (1) 𝑁𝐷𝑉𝐼 = Trong 𝑃 ,𝑃 giá trị phản xạ phổ kênh cận hồng ngoại kênh đỏ Khi đó, NDVI > 0,5 pixel coi hoàn toàn bao phủ thực vật (đối tượng thực vật), độ che phủ thực vật Nếu NDVI ≤ 0,2 pixel coi hoàn toàn bao phủ thổ nhưỡng (đối tượng thổ nhưỡng), độ che phủ thực vật Nếu 0,2 < NDVI ≤ 0,5 độ che phủ thực vật xác định theo công thức (4) Hình Sơ đồ nghiên cứu Phương pháp sử dụng GIS: Đây phương pháp quan trọng mang lại hiệu cao trình nghiên cứu, GIS thiết kế đồ chuyên đề, CSDL phục vụ cho đánh giá nguy TLĐ Ngoài modul GIS phân tích nội suy khơng gian tạo sở cho việc xác định, phân vùng nguy TLĐ Trong báo này, nhóm tác giả sử dụng phần mềm hãng ESRI AcrGIS 10.2 (phiên dùng thử) để tổng hợp lớp thông tin thành phần độ dốc, lượng mưa, thảm thực vật, … theo công thức (2) để xây dựng đồ nguy TLĐ vùng nghiên cứu Từ liệu thu thập được, xây dựng từ đồ thành phần dạng raster tích hợp thành lập đồ nguy trượt lở theo công thức (1) [2, 4, 5, 8]: LSI = ∑ 𝑊𝑗𝑋𝑖𝑗 (2) Trong LSI (Landslide Susceptibility Index) số nhạy cảm trượt lở đất; Wj trọng số nhân tố thứ j; Xij điểm số lớp thứ i nhân tố gây trượt j Sau tính tốn tiến hành kiển tra xem xét tỉ lệ quán (Consistency Ratio) Nếu CR ≤ 0,1 (10%) kết chấp nhận đánh giá chuyên gia tương đối quán Ngược lại, CR > 10% đánh giá không quán, phán đốn có phần ngẫu nhiên cần tiến hành đánh giá, xem xét lại Cơng thức (3) tính tỷ số quán [1–2, 5]: Hội nghị khoa học tồn quốc “Chuyển đổi số cơng nghệ số Khoa học Trái đất, Mỏ Môi trường” (EME 2021) Tạp chí Khí tượng Thủy văn 2022, EME4, 85-93; doi:10.36335/VNJHM.2022(EME4).85-93 88 𝐶𝑅 = (3) Trong CI số quán (Consistency Index); RI số ngẫu nhiên (Random Index) xác định theo bảng giá trị có sẵn; Trong bài, tỷ số quán lớp thành phần CR = 0,1 (10%) nên kết chấp thuận Phương pháp so sánh cặp (AHP – Analytical Hiearchy Process): Phương pháp Analytical Hierarchy Process (AHP) Saaty dựa vào việc xây dựng ma trận “so sánh cặp” để lựa chọn trọng số khai thác sử dụng triệt để Để định lượng hóa mức độ nhạy cảm phản ánh nguy TLĐ, tiến hành tích hợp tiêu theo cơng thức (2) (theo Patrono, et al., 1995) [1–2, 4–5] Bảng Thang tỉ lệ so sánh theo cặp (theo Saaty) [1–2] quan trọng 1/9 Ít cực độ 1/7 vơ quan trọng 1/5 Rất 1/3 Ít vừa Quan trọng Quan trọng Quan trọng nhiều Quan trọng nhiều Cực kỳ quan trọng Ma trận so sánh cặp (AHP): Để phân cấp hai tiêu chuẩn Saaty (1970) phát triển loại ma trận đặc biệt gọi ma trận so sánh cặp [1–2, 4–5] Những ma trận đặc biệt sử dụng để liên kết tiêu chuẩn đánh giá theo thứ tự thang phân loại Trong trình nghiên cứu tham vấn ý kiến số chuyên gia (8 chuyên gia) liên quan đến lĩnh vực nghiên cứu địa chất, tai biến môi trường, GIS nhằm thu thập thông tin ý kiến đánh giá so sánh yếu tố đưa vào tính tốn TLĐ, từ xây dựng ma trận trọng số Bảng Ma trận so sánh cặp yếu tố gây trượt lở Lớp thành phần Độ dốc Lượng mưa (a) 1/3 1/3 1/3 1/6 1/9 (a) (b) (c) (d) (e) (f) (b) 1/5 1/3 1/4 1/7 Hướng sườn (c) 1/2 1/2 1/5 Đứt gãy sông suối (d) 3 1/6 1/5 Thổ nhưỡng (e) 6 1/3 Lớp phủ (f) 5 Việc đánh giá mức độ nhạy cảm nhân tố thang điểm biểu thị ưu tiên chúng cách thích đáng q trình trượt lở đất Trong đó, độ dốc sườn nhân tố quan trọng bậc ảnh hưởng đến trượt lở đất, cho đểm Các nhân tố lại, ứng với mức độ nhạy cảm đến khả gây tai biến trượt lở đất khác trình bày bảng Để tính tốn trọng số cho thành phần sử dụng vector ngun lý Eigen (có thể tính tốn gần vector nguyên lý Eigen cách chia giá trị cột cho tổng giá trị cột) để thiết lập ma trận lúc giá trị trung bình hàng trọng số cho nhân tố (Bảng 4) [2, 6] Bảng Ma trận xác định trọng số yếu tố Lớp thành phần Độ dốc (a) Lượng mưa (b) Hướng sườn (c) Đứt gãy (d) Thổ nhưỡng (e) Lớp phủ (f) Trọng số (a) 0,44 0,61 0,29 0,32 0,31 0,30 0,38 (b) 0,15 0,20 0,49 0,32 0,21 0,23 0,27 (c) 0,15 0,04 0,10 0,21 0,10 0,17 0,13 Hội nghị khoa học toàn quốc “Chuyển đổi số công nghệ số Khoa học Trái đất, Mỏ Mơi trường” (EME 2021) Tạp chí Khí tượng Thủy văn 2022, EME4, 85-93; doi:10.36335/VNJHM.2022(EME4).85-93 Độ dốc Lớp thành phần (a) Lượng mưa (b) Hướng sườn (c) Đứt gãy (d) Thổ nhưỡng (e) 89 Lớp phủ (f) Trọng số (d) 0,15 0,07 0,05 0,11 0,31 0,17 0,14 (e) 0,7 0,05 0,05 0,02 0,05 0,10 0,06 (f) 0,05 0,03 0,02 0,02 0,02 0,03 0,03 Tổng 1,00 Kết thảo luận 3.1 Vai trò yếu tố Yếu tố độ dốc: Nghiên cứu sử dụng liệu DEM sau sử dụng cơng cụ Slope phần mềm ArcGIS tính tốn độ dốc, kết xây dựng đồ độ dốc (Bảng Hình 2b) Bảng Phân cấp ảnh hưởng độ dốc Cấp ảnh hưởng Độ dốc (0) Diện tích (ha) Tỷ lệ (%) 0–5 31339,33 36,70 – 15 49199,23 57,61 15 – 30 4744.20 5,55 30 – 45 110.69 0,13 >45 1,55 0,002 85395 100,00 Tổng Yếu tố hướng sườn: Hướng sườn liên quan đến khả tiếp nhận ánh sáng mặt trời, gió lượng mưa; từ ảnh hưởng gián tiếp đến nhân tố khác gây trượt lở đất độ ẩm đất, thảm phủ, tầng dày đất Sử dụng DEM cơng cụ Flow direction (hướng dịng chảy) ArcGIS để có hướng phơi sườn phục vụ cho u cầu (hình 3a) Hình (a) Mơ hình số độ cao DEM; (b) Bản đồ độ dốc Yếu tố lớp phủ: Thảm phủ thực vật rừng có vai trò lớp phủ che chắn hạn chế xói mịn đất, rễ thực vật làm tăng độ kết dính đất (Bảng 6) Bảng Phân cấp ảnh hưởng lớp phủ Cấp ảnh hưởng Lớp phủ Diện tích (ha) Tỷ lệ (%) Mặt nước 684,99 0,80 Rừng 6699,97 7,85 Cây lâu năm 10959,21 12,83 Trảng cỏ, bụi, 53159,32 62,25 Hội nghị khoa học tồn quốc “Chuyển đổi số cơng nghệ số Khoa học Trái đất, Mỏ Môi trường” (EME 2021) Tạp chí Khí tượng Thủy văn 2022, EME4, 85-93; doi:10.36335/VNJHM.2022(EME4).85-93 90 Cấp ảnh hưởng Lớp phủ Diện tích (ha) Tỷ lệ (%) Khu dân cư 13364,02 15,65 Đất trống 527,49 0,62 85395 100 Tổng Bằng phương pháp phân loại không giám sát phần mềm ENVI chia lãnh thổ nghiên cứu theo loại lớp phủ (mặt nước, rừng, lâu năm, trảng cỏ bụi, nơng nghiệp, dân cư, đất khơng có lớp phủ) (Hình 3b) Hình (a) Bản đồ hướng sườn; (b) Bản đồ thảm thực vật Thổ nhưỡng: Các loại đất chủ yếu địa bàn huyện đất xám phù sa cổ, đất dốc tụ, đất nâu đỏ đá bazan, đất nâu vàng bazan, đất xám gley, đất nâu vàng phù sa cổ có diện tích lớn 36.627,78 chiếm 42,89% đất xám gley có diện tích nhỏ 278,24 chiếm 0,33% (Bảng Hình 4a) Bảng Phân cấp yếu tố thổ nhưỡng Cấp ảnh hưởng Loại đất Diện tích (ha) Tỷ lệ (%) 2 Tổng Đất dốc tụ Đất nâu đỏ đá bazan Đất nâu vàng đá bazan Đất nâu vàng phù sa cổ Đất xám gley Đất xám phù sa cổ 3593,23 25212,06 8690,00 36627,78 278,24 10993,68 85395 4,21 29,52 10,18 42,89 0,33 12,87 100 Lượng mưa: Mùa mưa mùa phát sinh sạt lở đất mạnh nhất, thời gian mưa dài ngày lượng mưa lớn dịng chảy mặt phát triển mạnh làm gia tăng tượng sạt lở bờ sơng Bên cạnh đó, nước mưa làm tăng tải trọng khối đất đá sườn dốc làm giảm liên kết chúng gây nên trượt lở đất Dựa vào số liệu trạm đo mưa Phước Long, Đồng Xồi (Bình Phước) Tây Ninh, Gị Dầu (Tây Ninh) cơng cụ nội suy phần mềm ArcMap thành lập đồ mưa trung bình năm (Hình 4b) Bảng Phân cấp ảnh hưởng đứt gãy sông suối Cấp ảnh hưởng Khoảng cách (km) Diện tích (ha) Tỷ lệ (%) Tổng 28620,78 13434,90 10150,10 8720,95 24468,27 85395 33,52 15,73 11,89 10,21 28,65 100 Hội nghị khoa học toàn quốc “Chuyển đổi số công nghệ số Khoa học Trái đất, Mỏ Mơi trường” (EME 2021) Tạp chí Khí tượng Thủy văn 2022, EME4, 85-93; doi:10.36335/VNJHM.2022(EME4).85-93 91 Hình (a) Bản đồ thổ nhưỡng; (b) Bản đồ lượng mưa; (c) Bản đồ mật độ đứt gãy; (d) Bản đồ nguy trượt lở đất Mật độ đứt gãy: Sau số hoá, chuyển đổi đổi tách lớp liệu mặt nước thành lập lớp liệu sông suối khu vực nghiên cứu (Hình 4c) 3.2 Đánh giá nguy trượt lở đất Đánh giá nguy trượt lở đất dựa đồ nguy trượt lở đất khu vực nghiên cứu, đồ xây dựng dựa việc tích hợp tiêu (bản đồ thành phần) raster với độ phân giải pixel = 12.5×12.5 m Dùng Raster calculator phần mềm ArcGIS Chỉ số nguy trượt lở đất (LSI) tính sau [1–2, 5–7] 𝐿𝑆𝐼 = 0.38 ∗ 𝑆𝐷 + 0.27 ∗ 𝑆𝐷 + 0.13 ∗ 𝑆𝐷 + 0.14 ∗ 𝑆𝐷 + 0.06 ∗ 𝑆𝐷 + 0.03 ∗ 𝑆𝐷 Sau tính tốn ta thành lập đồ nguy trượt lở đất khu vực nghiên cứu tiến hành phân cấp nguy trượt lở đất ta đồ nguy trượt lở khu vực nghiên cứu cách sử dụng công thức 4: _ _ ∆ = (4) làm bước nhảy ta chia giá trị thành khoảng gồm cấp độ thấp, Sử dụng ∆ thấp, trung bình, cao cao Từ bảng hình 4d cho thấy vùng có nguy trượt lở thấp thấp chiếm phần lớn diện tích khu vực Trong khu vực có nguy trượt lở đất thấp chiếm 24.719,81 (chiếm 28,95%), nguy trượt lở thấp 28.889,21 (chiếm 33,83%), nguy trượt lở trung bình 5.751,43 (chiếm 6,74%), khu vực có nguy trượt lở cao 20.660,45 (chiếm 24,19%) cao 5.375,09 (chiếm 6,29%) Bảng Diện tích cấp trượt lở STT Khả trượt lở Diện tích (ha) Phần trăm (%) Rất thấp Thấp Trung Bình 24718,81 28889,21 5751,43 28,95 33,83 6,74 Hội nghị khoa học tồn quốc “Chuyển đổi số cơng nghệ số Khoa học Trái đất, Mỏ Môi trường” (EME 2021) Tạp chí Khí tượng Thủy văn 2022, EME4, 85-93; doi:10.36335/VNJHM.2022(EME4).85-93 STT Khả trượt lở Diện tích (ha) Phần trăm (%) Tổng Cao Rất cao 20660,45 5375,09 85.395 24,19 6,29 100,00 92 Kết luận Việc sử dụng GIS phương pháp thứ bậc (AHP) đánh giá nguy trượt lở đất cho phép đánh giá cách hiệu quả, khách mang tính đại diện Trong nghiên cứu này, có yếu tố sử dụng Tuy nhiên, tác động yếu tố ảnh hưởng đến q trình trượt lở đất khác nhau, thứ nguyên yếu tố không giống Việc kết hợp yếu tố không đồng thứ nguyên kể ưu điểm phương pháp Nghiên cứu xác định khu vực có nguy xảy TLĐ xây dựng đồ nguy trượt lở đất, sở để huyện Lộc Ninh xây dựng kế hoạch giải pháp phòng tránh Do hạn chế thời gian, nên nghiên cứu sử dụng 06 yếu tố chính, nghiên cứu cần bổ sung thêm yếu tố khác để xác định kết có tính thực tế Đóng góp tác giả: Xây dựng ý tưởng nghiên cứu: N.H.A., N.T.K.; Điều tra, khảo sát, phân tích số liệu: N.T.K., T.M.H., N.V.A.; Viết thảo báo: N.H.A., N.T.K.; Chỉnh sửa báo: N.H.A Lời cam đoan: Tập thể tác giả cam đoan báo cơng trình nghiên cứu tập thể tác giả, chưa công bố đâu, không chép từ nghiên cứu trước đây; khơng có tranh chấp lợi ích nhóm tác giả Tài liệu tham khảo Tuấn, T.A.; Dần, N.T Nghiên cứu nhạy cảm phân vùng nguy trượt – lở đất khu vực hồ thuỷ điện Sơn La theo phương pháp phân tích cấp bậc Saaty Tạp chí Khoa học Trái đất 2012, 34(3), 223–232 Anh, N.H Tích hợp GIS, viễn thám phương pháp phân tích thứ bậc Saaty (AHP) nghiên cứu trượt lở đất huyện Phú Lộc, tỉnh Thừa Thiên Huế Tạp chí khoa học, Đại học Huế – chuyên san Khoa học Trái đất Môi trường 2015, 13(4) Hội khảo sát địa chất Hoa Kỳ Online avaliable: http://earthexplorer.usgs.gov Hà, T.T Nghiên cứu địa mạo phục vụ giảm nhẹ thiệt hại tai biến trượt lở đất, lũ bùn đá tỉnh Lào Cai Luận án Tiến sỹ, Đại học Quốc gia Hà Nội, 2010 Thạch, N.N Áp dụng viễn thám hệ thông tin địa lý (GIS) để nghiên cứu dự báo tai biến thiên nhiên tỉnh Hịa Bình Báo cáo tổng kết đề tài Nghiên cứu Khoa học, Đại học Quốc gia Hà Nội, 2002 Ngọc, Đ.M.; Thùy, Đ.T.; Đức, Đ.M Ứng dụng GIS phương pháp phân tích thứ bậc (AHP) thành lập đồ nguy trượt lở huyện Xín Mần, tỉnh Hà Giang, Việt Nam Tạp chí Khoa học, Đại học Quốc gia Hà Nội: chuyên san Khoa học Trái đất Môi trường 2016, 32(2S), 206–216 Hà, Đ.T.; Trang, B.T.T.; Thành, N.K Nghiên cứu xây dựng đồ nguy sạt lở đất cho huyện Văn Yên, tỉnh Yên Bái Tạp chí Khoa học học công nghệ Việt Nam tập, 2019, 61(3) UBND tỉnh Bình Phước Địa chí Bình Phước, tập NXB Chính trị Quốc gia – thật, Hà Nội, 2015 Alaska Satellite Facility, Geophysical Institute, University of Alaska Fairbanks Dữ liệu DEM Online avaliable : https://search.asf.alaska.edu Hội nghị khoa học tồn quốc “Chuyển đổi số cơng nghệ số Khoa học Trái đất, Mỏ Môi trường” (EME 2021) Tạp chí Khí tượng Thủy văn 2022, EME4, 85-93; doi:10.36335/VNJHM.2022(EME4).85-93 93 Studying landslide hazards in Loc Ninh district, Binh Phuoc province, using GIS and the analytic hierarchy process (AHP) Nguyen Huy Anh1*, Nguyen Trung Kien1, Tran My Hao1 Ho Chi Minh University of Natural Resources and Environment; anhnh@hcmunre.edu.vn; tranhao_tnmt@yahoo.com Abstract: The goal of this research is to use AHP and GIS technologies to create a landslide risk map in the Loc Ninh district of Binh Phuoc province According to the findings, landslide risk is divided into five levels: very low, with an area of 24,718,81 accounting for 28.95 percent; low risk, with an area of 28,889.21 hectares accounting for 62.78 percent of the total area of Loc Ninh district; moderate risk, with 5,751.43 hectares accounting for 6.74 percent; high risk, with 20,660.45 hectares accounting for 24.19 percent; very high risk, with 5,375 hectares Keywords: Landslide; Analytic Hierarchy Process; GIS; Remote sensing; Loc Ninh District Hội nghị khoa học tồn quốc “Chuyển đổi số cơng nghệ số Khoa học Trái đất, Mỏ Môi trường” (EME 2021)