BỘ GIÁO DỤC VÀ ĐÀO TẠO BỘ XÂY DỰNG VIỆN KHOA HỌC VÀ CÔNG NGHỆ XÂY DỰNG  NGUYỄN QUANG HUY NGHIÊN CỨU ĐÁNH GIÁ NGUY CƠ TRƢỢT ĐẤT VÀ LUẬN CHỨNG HỆ THỐNG QUAN TRẮC PHỤC VỤ CẢNH BÁO TAI BIẾN TRƢỢT ĐẤT CHO KHU VỰC TÂY NAM TỈNH HÀ GIANG LUẬN ÁN TIẾN SĨ ĐỊA CHẤT CÔNG TRÌNH Mã số: 62.44.65.01 HÀ NỘI, 2017 BỘ GIÁO DỤC VÀ ĐÀO TẠO BỘ XÂY DỰNG VIỆN KHOA HỌC VÀ CÔNG NGHỆ XÂY DỰNG  NGUYỄN QUANG HUY NGHIÊN CỨU ĐÁNH GIÁ NGUY CƠ TRƢỢT ĐẤT VÀ LUẬN CHỨNG HỆ THỐNG QUAN TRẮC PHỤC VỤ CẢNH BÁO TAI BIẾN TRƢỢT ĐẤT CHO KHU VỰC TÂY NAM TỈNH HÀ GIANG LUẬN ÁN TIẾN SĨ CHUYÊN NGÀNH ĐỊA CHẤT CÔNG TRÌNH Mã số : 62.44.65.01 NGƢỜI HƢỚNG DẪN KHOA HỌC: PGS.TSKH Trần Mạnh Liểu PGS.TS Đoàn Thế Tƣờng HÀ NỘI, 2017 LỜI CAM ĐOAN Tôi xin cam đoan công trình nghiên cứu riêng Các kết nghiên cứu luận án tính toán đƣa hoàn toàn trung thực chƣa đƣợc công bố công trình khác LỜI CÁM ƠN Luận án đƣợc hoàn thành Viện Khoa học Công nghệ xây dựng dƣới hƣớng dẫn khoa học PGS.TSKH Trần Mạnh Liểu PGS.TS Đoàn Thế Tƣờng Tôi xin bày tỏ lòng biết ơn sâu sắc thày giáo hƣớng dẫn, ngƣời gợi mở tƣ nghiên cứu có đóng góp to lớn trình thực hoàn thiện luận án Tôi xin chân thành cám ơn tổ chức, cán thuộc Viện Khoa học Công nghệ xây dựng, Trung tâm nghiên cứu đô thị - ĐHQGHN, Viện Khoa học Địa chất & Khoáng sản, Sở Khoa học Công nghệ tỉnh Hà Giang, Công ty tƣ vấn lập Xí nghiệp Xây dựng Chuyển giao Công nghệ Môi trƣờng tạo điều kiện thuận lợi, quan tâm, giúp đỡ nhiệt tình suốt trình thực luận án Xin trân trọng cám ơn chuyên gia, thày cô, giáo, nhà khoa học tham gia hội thảo mở rộng, hội đồng đánh giá bƣớc thực luận án, đóng góp ý kiến quý báu, tâm huyết để hoàn thiện luận án Tất giúp đỡ nêu trên, ghi nhớ trân trọng mang theo suốt trình học tập, công tác nghiên cứu Nghiên cứu sinh Nguyễn Quang Huy MỤC LỤC PHẦN MỞ ĐẦU 1 Tính cấp thiết đề tài, lĩnh vực nghiên cứu Mục tiêu nghiên cứu Đối tƣợng phạm vi nghiên cứu Nội dung nghiên cứu Phƣơng pháp nghiên cứu Luận điểm bảo vệ Những điểm khoa học Cơ sở tài liệu Cấu trúc luận án CHƢƠNG I TỔNG QUAN CÁC VẤN ĐỀ NGHIÊN CỨU VỀ TRƢỢT ĐẤT ĐÁ 1.1 Khái niệm trƣợt đất đá Phân loại trƣợt đất đá 1.3 Các yếu tố nguyên nhân, điều kiện gây trƣợt đất 13 1.3.1 Nguyên nhân gây trƣợt đất 13 1.3.2 Điều kiện thành tạo trƣợt 16 1.4 Vai trò (tỷ trọng tham gia) yếu tố điều kiện thành tạo nguyên nhân gây trƣợt 17 1.5 Tình hình nghiên cứu đánh giá dự báo nguy trƣợt đất hệ thống quan trắc phục vụ cảnh báo tai biến trƣợt đất giới 18 1.6 Tình hình nghiên cứu đánh giá dự báo nguy trƣợt đất Việt Nam số phƣơng pháp phân vùng nguy tai biến trƣợt đất đƣợc áp dụng 24 1.6.1 Các hƣớng nghiên cứu đƣợc triển khai 24 1.6.2 Một số mô hình tính toán đánh giá nguy trƣợt đất đƣợc áp dụng .27 CHƢƠNG II ĐẶC ĐIỂM ĐIỀU KIỆN HỆ TỰ NHIÊN – KINH TẾ - XÃ HỘI KHU VỰC TÂY NAM TỈNH HÀ GIANG 45 2.1.Vị trí địa lý 45 2.2 Đặc điểm khí hậu 46 2.3 Đặc điểm địa hình, độ đốc 49 2.3.1.Yếu tố địa hình 49 2.3.2 Yếu tố độ dốc 53 2.4 Đặc điểm địa chất, thạch học 56 2.5 Vỏ phong hóa 59 2.6 Hiện trạng sử dụng đất khu vực nghiên cứu 61 2.7 Hệ thống đƣờng giao thông 63 2.8 Đặc điểm đứt gẫy, kiến tạo 64 2.9 Hoạt động kinh tế - xã hội 67 2.10 Hiện trạng trƣợt đất khu vực nghiên cứu 68 2.11 Đặc điểm quy luật phát sinh, phát triển trƣợt đất 73 CHƢƠNG III ĐÁNH GIÁ DỰ BÁO NGUY CƠ TRƢỢT ĐẤT 85 3.1 Lựa chọn phƣơng pháp đánh giá dự báo nguy trƣợt đất cấp khu vực 85 3.2 Đánh giá nguy trƣợt đất khu vực Tây Nam tỉnh Hà Giang 85 3.2.1 Cơ sở liệu phƣơng pháp xây dựng đồ thành phần phân lớp yếu tố gây trƣợt xác định trọng số tƣơng ứng lớp 85 3.2.2 Kết tính toán, xây dựng đồ thành phần phân lớp yếu tố gây trƣợt trọng số tham gia lớp phân chia tƣơng ứng 89 3.2.3 Xác định trọng số yếu tố gây trƣợt đất (Wj) 104 3.2.4 Bản đồ đánh giá nguy trƣợt đất 104 CHƢƠNG IV: 108 LUẬN CHỨNG HỆ THỐNG QUAN TRẮC PHỤC VỤ CẢNH BÁO TAI BIẾN TRƢỢT ĐẤT KHU VỰC NGHIÊN CỨU 108 4.1 Tổng quan hệ thống quan trắc phục vụ cảnh báo tai biến trƣợt đất 108 4.1.1 Phân cấp hệ thống quan trắc 108 4.1.2 Mục tiêu nhiệm vụ hệ thống quan trắc 108 4.1.3 Quan trắc (giám sát) trạng trƣợt đất toàn khu vực 112 4.1.4 Mạng lƣới quan trắc đo vẽ điểm trƣợt 112 4.1.5 Sơ đồ nguyên tắc tổ chức hệ thống quan trắc trƣợt đất 114 4.2 Luận chứng hệ thống quan trắc cho khu vực nghiên cứu 115 4.2.1 Cấu trúc hệ thống quan trắc 115 4.2.2 Mạng lƣới tuyến, trạm quan trắc 117 4.2.3 Tổng hợp khối lƣợng quan trắc 119 4.2.4 Phƣơng pháp, thiết bị quan trắc 119 4.2.4.1 Quan trắc phân bố không gian khối trƣợt 120 4.2.4.2 Quan trắc đo vẽ đặc điểm khối trƣợt yếu tố điều kiện gây trƣợt chủ yếu 122 4.2.5 Quan trắc yếu tố tác động, biến đổi nhanh 123 4.3 Định hƣớng mô hình cảnh báo trƣợt đất khu vực nghiên cứu từ liệu quan trắc 125 4.3.1 Điều chỉnh trọng số mô hình số thống kê tích hợp đa biến 125 4.3.2 Sử dụng mô hình định lƣợng tiêu tích hợp yếu tố điều kiện, nguyên nhân gây trƣợt đất để dự báo phân vùng nguy trƣợt đất 125 4.3.3 Cảnh báo nhanh quy mô, cƣờng độ (thể tích) khối trƣợt xuất thông qua quan trắc yếu tố gây trƣợt chủ yếu hàm hồi quy 126 4.3.4 Mô hình cảnh báo, dự báo ngƣỡng phát triển trƣợt đất theo yếu tố có tỷ trọng gây trƣợt lớn, biến đổi nhanh hàm hồi quy 128 KẾT LUẬN 130 TÀI LIỆU THAM KHẢO 133 DANH MỤC HÌNH VẼ Hình 1: Khu vực nghiên cứu (giới hạn phạm vi đường viền màu đỏ) Hình 2: Ảnh số khối trƣợt khu vực nghiên cứu Hình 1.1: Đặc điểm cấu trúc khối trƣợt 12 Hình 1.2: Sơ đồ biểu diễn lực tác động lên khối trƣợt (theo V.D Lomtadze) 14 Hình 1.3: Phân loại dự báo nguy trƣợt đất Kunxen V.V; Postoev G.P, Gylakian K.A có bổ sung (theo tác giả) 22 Hình 1.4: Sơ đồ thể liệu đầu vào kết đầu Sinmap 28 Hình 1.5: Sơ đồ ổn định mái dốc 29 Hình 1.6: Sơ đồ mái dốc chiếu chiếu đứng 31 Hình 1.7: Quy trình áp dụng Mô hình Chỉ số Thống kê dự báo khả trƣợt đất ArcGIS 35 Hình 2.1: Bản đồ mô hình số độ cao khu vực nghiên cứu 52 Hình 2.2: Mô hình số độ dốc khu vực nghiên cứu 55 Hình 2.3: Bản đồ phân bố địa chất thạch học khu vực nghiên cứu 58 Hình 2.4: Bản đồ phân bố đứt gãy khu vực nghiên cứu 65 Hình 2.5: Biểu đồ phân loại điểm trƣợt theo quy mô trƣợt đất 70 Hình 2.6: So sánh quy mô trƣợt đất sau chu kỳ 12 năm (2003-2015) 72 Hình 2.7: Bản đồ trạng trƣợt đất khu vực nghiên cứu 73 Hình 2.8: Khối trƣợt chảy khu vực Bản Nà Lống, xã Bản Díu, XM 74 Hình 2.9: Mô tả khối trƣợt chảy khu vực Bản Nà Lống, xã Bản Díu, Xín Mần 75 Hình 2.10: Khối trƣợt hỗn hợp Bản Nà Lống, xã Bản Díu, Xín Mần 76 Hình 2.11: Mô tả khối trƣợt hỗn hợp Cẩm Phu, xã Bản Phùng, HSP 77 Hình 2.12: Trƣợt chảy khu vực Kết Thành, Hồ Thầu, HSP 79 Hình 2.13: Trƣợt chảy khu vực Bản Péo, HSP 80 Hình 2.14: Trƣợt chảy khu vực Bản Chè, Tân Tiến, HSP 81 Hình 2.15: Trƣợt hỗn hợp khu vực Pố Lồ, HSP 82 Hình 2.16: Trƣợt hỗn hợp khu vực km 85+588, XM 83 Hình 3.1: Các yếu tố đƣa vào tính toán nguy trƣợt đất 87 Hình 3.2: Bảng phân loại cao độ địa hình theo phƣơng pháp Natural Break 90 Hình 3.3: Bản đồ phân lớp cao độ địa hình giá trị trọng số tƣơng ứng 91 Hình 3.4: Bảng phân loại độ dốc theo phƣơng pháp Nature Break 92 Hình 3.5: Bản đồ độ dốc giá trị trọng số tƣơng ứng 93 Hình 3.6: Bản đồ hƣớng dốc giá trị trọng số tƣơng ứng 94 Hình 3.7: Bản đồ mật độ phân cắt ngang giá trị trọng số tƣơng ứng 95 Hình 3.8: Bản đồ mật độ phân cắt sâu giá trị trọng số tƣơng ứng 96 Hình 3.9: Bản đồ khoảng cách đến đƣờng giao thông trọng số tƣơng ứng 97 Hình 3.10 Bản đồ phân bố khoảng cách đến đứt gãy trọng số tƣơng ứng 98 Hình 3.11: Bản đồ phân bố địa chất thạch học trọng số tƣơng ứng 100 Hình 3.12 Bản đồ phân bố lƣợng mƣa trọng số tƣơng ứng 101 Hình 3.13: Bản đồ cấu sử dụng đất trọng số tƣơng ứng 102 Hình 3.14: Bản đồ chiều dày vỏ phong hóa giá trị trọng số tƣơng ứng 103 Hình 3.15: Bản đồ nguy trƣợt đất khu vực nghiên cứu 106 Hình 3.16: Bản đồ nguy trƣợt đất khu vực điểm trƣợt 107 Hình 4.1: Sơ đồ nguyên tắc tổ chức hệ thống quan trắc trƣợt đất 114 Hình 4.2: Bản đồ phân bố tuyến trạm quan trắc 118 Hình 4.3: Nguyên tắc hoạt động ống ngắm đo khoảng cách 122 Hình 4.4: Diễn biến thay đổi yếu tố mƣa khu vực nghiên cứu 130 DANH MỤC BẢNG BIỂU Bảng 1.1: Hệ thống phân loại trƣợt đất đá theo Lomtadze V.D (1970) Bảng 1.2: Hệ thống phân loại trƣợt đất theo Varnes (1978) 11 Bảng 1.3: Phân loại mức độ ổn định sƣờn theo SI 32 Bảng 1.4: Chỉ tiêu Saaty so sánh cặp đôi yếu tố 41 Bảng 1.5: Ma trận kết so sánh tầm quan trọng yếu tố 41 Bảng 1.6: Đánh giá điểm số mức độ ảnh hƣởng yếu tố 42 hình thành phát triển tai biến 42 Bảng 2.1: Lƣợng mƣa tháng năm trạm Hà Giang 47 Bảng 2.2: Lƣợng mƣa tháng năm trạm Bắc Quang 47 Bảng 2.3: Lƣợng mƣa tháng năm trạm Hoàng Su Phì 48 Bảng 2.4: Tỷ lệ diện tích loại đất dốc khu vực huyện 53 Bảng 2.5: Hiện trạng cấu sử dụng đất 61 (nguồn:[1], [7]) 61 Bảng 2.6: Phân cấp mật độ đƣờng giao thông 64 Bảng 2.7: Đặc điểm phân bố mật độ đứt gãy mật độ điểm trƣợt tƣơng ứng khu vực nghiên cứu 66 Bảng 2.8: Quy mô trƣợt đất khu vực cấp xã, thị trấn 70 Bảng 3.1: Phân lớp cao độ địa hình giá trị trọng số tƣơng ứng 91 Bảng 3.2: Phân loại độ dốc giá trị trọng số tƣơng ứng 93 Bảng 3.3: Phân lớp hƣớng dốc giá trị trọng số tƣơng ứng 94 Bảng 3.4: Phân lớp mật độ phân cắt ngang giá trị trọng số tƣơng ứng 95 Bảng 3.5: Phân lớp mật độ phân cắt sâu giá trị trọng số tƣơng tứng 96 Bảng 3.6: Phân lớp khoảng cách đến đƣờng giao thông trọng số lớp 97 Bảng 3.7: Phân lớp nhóm khoảng cách đến đứt gãy trọng số tƣơng ứng 98 Bảng 3.8: Phân lớp thành phần thạch học trọng số tƣơng ứng 99 Bảng 3.9 Phân lớp lƣợng mƣa trọng số tƣơng ứng 100 Bảng 3.10: Phân lớp sử dụng đất giá trị trọng số 102 Bảng 3.11: Phân lớp chiều dày vỏ phong hóa giá trị trọng số tƣơng ứng 103 Bảng 3.12: Kết tính toán trọng số yếu tố 104 Bảng 4.1: Số phân lớp yếu tố phát sinh, phát triển trƣợt đất 117 Bảng 4.2: Tổng hợp khối lƣợng, chu kỳ quan trắc theo hệ thống quan trắc thành phần theo năm 119 Bảng 4.3: Thống kê lƣợng mƣa trung bình đo trạm khí tƣợng 128 (Advanced Land Observing Satellite-2 or DAICHI-2) hãng JAXA (The Japan Aerospace Exploration Agency) Phiên nâng cấp từ đầu thu PALSAR – đầu thu kênh L bƣớc sóng 23,6m, lắp vệ tinh ALOS với qu đạo bay thời gian xác để theo dõi biến dạng chậm bề mặt trái đất ALOS đóng góp vào việc đo vẽ, quan sát khu vực, theo dõi tai biến trƣợt đất khảo sát tài nguyên từ phóng lên trời năm 2006 Các ƣu điểm đầu thu kênh L so với kênh C kênh X nhƣ sau: Việc xuyên sâu qua khu vực có thực vật bao phủ làm giảm tƣơng quan theo thời gian, cho phép tăng thời gian thu chụp cặp ảnh; bƣớc sóng dài làm tăng khoảng đƣờng sở tới hạn, điều giúp tăng số lƣợng ảnh sử dụng đƣợc Vệ tinh ALOS-2 giúp đạt đƣợc thành công quan trắc trƣợt đất với khả đƣợc cải tiến đầu thu PALSAR2 có đặc tính đặc biệt cho phép hoạt động nhiều phƣơng thức quan sát với độ phân giải khác Ảnh vệ tinh dùng để quan trắc trƣợt đất cho khu vực nghiên cứu có độ phân giải cao, đƣợc chụp theo phƣơng thức Spotlight, độ phân giải 1m (azimuth) x 3m (range), kích thƣớc vùng chụp 25x25km Quá trình phân tích vị trí trƣợt đất khu vực nghiên cứu ảnh viễn thám song song tiến hành Với số liệu ảnh viễn thám đƣợc sử dụng ảnh vệ tinh với độ phân giải 3m đồ địa hình tỷ lệ 1:50.000 Bản đồ địa hình đƣợc số hóa với khoảng cách đƣờng đồng mức 10m Sau sở đồ địa hình đƣợc gán giá trị độ cao đƣờng đồng mức, thuật toán nội suy đƣợc sử dụng để nội suy tạo mô hình số độ cao (cao độ địa hình) khu vực nghiên cứu Để hiển thị phân tích ảnh vệ tinh dƣới kính lập thể (Stereo), ảnh vệ tinh khu vực nghiên cứu đƣợc nắn chỉnh chồng chập lên DEM để tạo ảnh khu vực nghiên cứu theo không gian chiều (hình 4.4) Trên sở ảnh 3D, tất vị trí trƣợt đất đƣợc phân tích, giải đoán, đánh dấu mắt thƣờng không gian chiều thông qua đặc điểm địa hình, dạng kiến trúc dòng chảy, xói mòn bề mặt, tông ảnh, thảm phủ trạng sử dụng đất 121 4.2.4.2 Quan trắc đo vẽ đặc điểm khối trƣợt yếu tố điều kiện gây trƣợt chủ yếu Các yếu tố thể tích, diện tích khối trƣợt (đặc điểm khối trƣợt), khoảng cách đến đƣờng giao thông, khoảng cách đến đứt gãy, độ cao địa hình, độ dốc, hƣớng dốc, phân cắt ngang, phân cắt sâu, thành phần thạch học, chiều dày vỏ phong hóa (các yếu tố điều kiện gây trƣợt chủ yếu) đƣợc quan trắc đo thiết bị đo quan sát, ghi chép thực địa Tại điểm khảo sát, tiến hành định vị tọa độ địa lý điểm khảo sát độ cao địa hình thiết bị GPS Garmin 62 (độ xác ± 3m) GPS 72 (độ xác ± 5m) Đo thể tích khối trƣợt, chiều dày vỏ phong hóa xuất lộ mặt trƣợt ống ngắm đo khoảng cách Nikon Laser 550A S Khoảng cách đo tối đa 1000m góc đo tối đa 89o độ phóng đại lên tới lần Với chiều dài dƣới 30m dùng thƣớc dây để đo Hình 4.3: Nguyên tắc hoạt động ống ngắm đo khoảng cách 122 Quan trắc mái dốc, độ dốc, hƣơng dốc khu vực nghiên cứu máy đo Nikon Laser 550A S Quan trắc mức độ phân cắt ngang, phân cắt sâu, khoảng cách đến đứt gãy, đến đƣờng giao thông giải đoán ảnh vệ tinh độ phân giải cao, đồ địa hình tỷ lệ 1/50.000 kết hợp máy đo Nikon Laser 550A S Quan trắc chiều dày vỏ phong hóa thông qua giải đoán đồ địa chất tỷ lệ 1/50.000 Mô tả thành phần vật chất khối trƣợt, đặc điểm, chế trƣợt đất thực tế kết hợp vấn ngƣời dân để biết thêm thông tin xung quanh khối trƣợt, đƣa bảng thông tin khối trƣợt lân cận khu vực quan trắc Nhận định khối trƣợt theo bảng phân loại trƣợt đất Lomtazde (1979) để xếp loại phân cấp khối trƣợt 4.2.5 Quan trắc yếu tố tác động, biến đổi nhanh 4.2.5.1 Thiết bị đo lƣợng mƣa Yếu tố lƣợng mƣa đƣợc quan trắc thông qua việc thu thập liệu lƣợng mƣa từ trạm khí tƣợng có khu vực nghiên cứu (Bắc Quang, Hoàng Su Phì, Xín Mần) kết hợp với thiết bị đo lƣợng mƣa đƣợc lắp đặt trạm quan trắc Tại trạm quan trắc, máy đo mƣa đƣợc lắp đặt hệ thống ghi đo tự động mô-đem điện thoại kết nối máy tính PC đặt từ xa với hệ thống ghi đo tự động Máy đo mƣa TE525 MM hãng Texas Electronics chế tao M loại máy đo mƣa hoạt động theo nguyên lý cốc lật điện từ Nƣớc mƣa đƣợc hứng qua phễu có đƣờng kính 25.4 mm xuống cốc chứa đƣợc hiệu chuẩn tƣơng ứng với lƣợng mƣa 0.1 mm Khi nƣớc mƣa đƣợc hứng qua phễu chảy xuống đầy cốc, cốc lật, làm bật công tắc điện từ, gây xung điện, xung điện đƣợc ghi đo máy ghi đo tự động 4.2.5.2 Thiết bị đo thay đổi cấu sử dụng đất Việc kết hợp ảnh vệ tinh độ phân giải cao với công nghệ laser scanner tăng độ phân giải không gian thời gian ảnh chụp bề mặt khu 123 vực nghiên cứu Công nghệ có tính linh hoạt cao, tốc độ thu thập xử lý số liệu nhanh chóng so với phƣơng pháp viễn thám truyền thống Với công nghệ này, tốc độ dịch chuyển bề mặt đƣợc xác định sở so sánh số liệu mô hình số bề mặt đƣợc xác định từ ảnh chụp ghi nhận thời điểm khác Các công tác tiến hành tự động thủ công Công nghệ sử dụng UAV (Thiết bị bay chụp điều khiển từ xa) cho phép phân tích chi tiết hình ảnh bề mặt phạm vi vùng nghiên cứu nhỏ (cấp xã,thị trấn) cho ta biết đặc điểm đất loại đất diện tích loại đất, từ tích hợp hình ảnh cho phép thành lập đồ trạng cấu sử dụng đất phạm vi toàn khu vực nghiên cứu 4.2.5.3 Thiết bị đo thay đổi mực nƣớc ngầm Thiết bị đo mực nƣớc ngầm Well Sounder 2010 PRO hãng Enoscientific – USA sản xuất Thiết bị đo sử dụng sóng âm để đo khoảng cách từ miệng giếng đến bề mặt nƣớc ngầm Khi đo, đặt đầu đo vào lỗ nắp giếng bật máy, vài giây, giá trị đo mực nƣớc ngầm đƣợc hiển thị Ngoài ra, máy hiển thị thông tin khác nhƣ: tổng lƣợng nƣớc giếng, tốc độ thay đổi, nhiệt độ điện áp pin Đặc điểm thiết bị là: tự động bù nhiệt; tính toán khoảng cách mực nƣớc tĩnh mực nƣớc động, tốc độ phục hồi; lấy số liệu lập tức; ghi số liệu cho phép ghi tới 25 triệu giá trị đo với thông tin thời gian đo, mã số giếng, nhiệt độ thông tin khác; ghi số liệu tự động, đặt thời khoảng lấy mẫu từ giây đến 60 phút; ghi số liệu thủ công: lần bấm phím lần lấy mẫu Số liệu ghi đƣợc đƣợc xem hình tải máy tính qua cổng USB RS232 Đơn vị đo: ft, in, m; độ phân giải: cm; độ xác: cm; dải đo: ~ 600 m Thời gian cập nhật: giây 150 m / giây 600 m 124 4.3 Định hƣớng mô hình cảnh báo trƣợt đất khu vực nghiên cứu từ liệu quan trắc 4.3.1 Điều chỉnh trọng số mô hình số thống kê tích hợp đa biến Để tính toán định lƣợng trọng số yếu tố điều kiện, nguyên nhân gây trƣợt, luận án sử dụng công thức kết nghiên cứu Trung tâm Nghiên cứu đô thị - Đại học Quốc gia Hà Nội, cụ thể Wj  MaxDij n  MaxDij (4-2) j 1 n: Số lƣợng yếu tố gây trƣợt đất khu vực nghiên cứu MaxDij : Mật độ trƣợt lớn phân lớp i yếu tố j Hệ thống quan trắc cung cấp liệu quan trắc thay đổi theo chu kỳ quan trắc Đối với liệu thuộc 01 chu kỳ quan trắc có đƣợc 01 giá trị MaxDij Wj Tùy kết quan trắc thực tế, giá trị trọng số yếu tố thay đổi dẫn đến việc điều chỉnh mạng lƣới quan trắc cho phù hợp, cụ thể: phân bố lại tuyến quan trắc cho phù hợp với đồ phân vùng cƣờng độ tác động yếu tố chiếm tỷ trọng lớn, phân bố lại trạm quan trắc cho phù hợp với đồ phân vùng nguy trƣợt đất đƣợc xác định lại sau sử dụng kết liệu quan trắc 4.3.2 Sử dụng mô hình định lƣợng tiêu tích hợp yếu tố điều kiện, nguyên nhân gây trƣợt đất để dự báo phân vùng nguy trƣợt đất Khu vực nghiên cứu đƣợc chia làm ô lƣới sở với kích thƣớc 20x20m (01 pixel) Chỉ tiêu tích hợp yếu tố phát sinh, phát triển trƣợt đất đƣợc tính toán cho tất ô lƣới sở mạng tính toán, sau xây dựng mô hình trƣờng biến đổi dƣới dạng đƣờng đẳng trị tiêu tích hợp Đó sở để tiến hành phân vùng lãnh thổ khả phát triển trƣợt đất Chỉ tiêu tích hợp yếu tố điều kiện, nguyên nhân gây trƣợt đất trƣợt đất đƣợc tính toán theo công thức sau: 125 (4-3) Trong đó: gi tỷ trọng yếu tố điều kiện, nguyên nhân gây trƣợt đất thứ i, tham số định lƣợng yếu tố điều kiện, nguyên nhân gây trƣợt đất thứ i đƣợc quan trắc Tỷ trọng (gi) yếu tố điều kiện, nguyên nhân gây trƣợt đất đƣợc xác định sở mối tƣơng quan yếu tố điều kiện, nguyên nhân gây trƣợt vói với đặc điểm trƣợt đất khu vực nghiên cứu (thông qua liệu quan trắc) Tham số định lƣợng yếu tố ( ), đƣợc điều chỉnh thông qua liệu quan trắc cƣờng độ tác động yếu tố điều kiện, nguyên nhân gây trƣợt Nhƣ vậy, dựa vào liệu quan trắc ta xây dựng đồ dự báo nguy trƣợt đất dƣới dạng đƣờng đẳng trị tiêu tích hợp 4.3.3 Cảnh báo nhanh quy mô, cƣờng độ (thể tích) khối trƣợt xuất thông qua quan trắc yếu tố gây trƣợt chủ yếu hàm hồi quy Xây dựng hàm hồi quy liên hệ thể tích khối trƣợt với yếu tố gây trƣợt chủ yếu Xác định yếu tố gây trƣợt chủ yếu dựa việc xác định định lƣợng tỷ trọng tham gia yếu tố gây trƣợt sở thống kê số liệu đo vẽ thể tích khối trƣợt yếu tố gây trƣợt tƣơng ứng - Tính toán hệ số tƣơng quan cặp đôi tất yếu tố gây trƣợt đƣợc xem xét (rij) yếu tố gây trƣợt với thể tích khối trƣợt, xây dựng ma trận chúng r21 r31 … … rp1 r12 r32……… rp ∆ = r13 r23 1… … rp … ……………… r1p r2p………… 126 - Lập hệ phƣơng trình: r1 y  1   r21    p rp1 r2 y  1 r12      p rp … rpy  1 r1 p   r2 p    p - Tính toán hệ số tiêu chuẩn βi Với  i  i ,  Hệ số i đƣợc xác định từ ma trận cách thay cột tƣơng ứng i hệ số tự phƣơng trình rij hệ số tƣơng quan yếu tố điều kiện tai biến thứ i j, r iy hệ số tƣơng quan yếu tố điều kiện thứ i hàm mục tiêu F TB (thể tích khối trƣợt) - Tính toán hệ số tƣơng quan nhiều chiều R p R =  i ry i (4-4) i 1 Hệ số tƣơng quan nhiều chiều cho phép xem xét tham số phát triển tai biến tham gia phân vùng có hợp lý hay không Thực tế hệ số tƣơng quan nhiều chiều R >0.7 tham số lựa chon đủ, hệ số tƣơng quan nhiều chiều nhỏ chắn việc xác định tham số phát triển tai biến thiếu số tham số quan trọng - Tính toán tỷ trọng tham số phát triển tai biến theo công thức sau: gi   i riy p r i 1 (4-5) i iy Mức độ quan trọng yếu tố xác định theo tỷ trọng chúng, tỷ trọng lớn, yếu tố quan trọng Phụ thuộc vào giá trị tỷ trọng 127 yếu tố gây trƣợt đƣợc tính toán khu vực nghiên cứu, xây dựng hàm hồi quy liên hệ thể tích khối trƣợt yếu tố chiếm tỷ trọng lớn Hàm hồi quy đƣa số liệu cảnh báo nhanh quy mô thể tích khối trƣợt xuất theo yếu tố chiếm tỷ trọng gây trƣợt lớn 4.3.4 Mô hình cảnh báo, dự báo ngƣỡng phát triển trƣợt đất theo yếu tố có tỷ trọng gây trƣợt lớn, biến đổi nhanh hàm hồi quy Căn bảng trọng số yếu tố định trƣợt đất (bảng , nhận thấy yếu tổ chiếm tỷ trọng gây trƣợt lớn lƣợng mƣa, khoảng cách đến đƣờng giao thông, Cao độ địa hình, trạng sử dụng đất Thực tế khảo sát cho thấy, khu vực nghiên cứu, năm qua đặc điểm độ cao địa hình, trạng sử dụng đất, hệ thống đƣờng giao thông biến đổi chậm Yếu tố lƣợng mƣa đƣợc ghi nhận yếu tố chiếm tỷ trọng gây trƣợt lớn biến động nhanh thất thƣờng qua năm Bảng 4.3: Thống kê lƣợng mƣa trung bình đo trạm khí tƣợng (nguồn:[40]) (Đơn vị: mm) Trạm Hà Giang Bắc Quang Năm Tổng 1991 191,8 50,4 400,2 1992 159,5 37,6 1993 215,9 1994 TB toàn khu vực Tổng Max 107,7 133,0 35,1 241,7 64,4 375,5 77,2 119,7 29,5 218,2 48,1 59,0 496,8 111,1 113,4 31,0 275,4 67,0 217,6 69,4 277,3 89,5 152,1 44,0 215,7 67,6 1995 224,7 50,1 447,8 87,3 133,9 33,2 268,8 56,9 1996 191,9 57,8 392,7 109,3 149,5 46,4 244,7 71,2 1997 177,4 45,5 498,2 117,2 117,0 29,3 264,2 64,0 1998 171,6 46,0 479,9 118,6 127,2 35,5 259,6 66,7 1999 230,6 55,7 515,4 138,1 121,4 30,4 289,1 74,7 2000 176,3 57,9 376,7 96,9 122,3 32,9 225,1 62,6 2001 188,5 47,4 353,7 96,0 111,9 36,3 218,0 59,9 Max Tổng Hoàng Su Phì Max 128 Tổng Max 2002 213,0 54,3 458,5 122,4 157,7 45,9 276,4 74,2 2003 220,9 52,5 439,6 96,2 130,6 35,6 263,7 61,4 2004 194,5 50,9 285,7 72,3 126,6 36,9 202,3 53,4 2005 204,8 59,9 306,0 80,2 136,7 36,0 215,8 58,7 2006 122,1 35,0 418,6 141,5 110,6 26,5 217,1 67,7 2007 246,2 48,2 226,9 63,8 135,2 33,6 202,8 48,5 2008 254,3 63,8 367,6 83,8 198,8 58,8 273,6 68,8 2009 195,1 57,5 351,5 87,4 112,0 40,6 219,5 61,8 2010 262,3 104,2 322,7 85,6 150,4 38,4 245,1 76,1 2011 241,8 60,3 361,5 95,1 199,0 52,3 267,4 69,2 2012 195,3 48,9 343,2 90,3 119,5 32,3 219,3 57,2 2013 212,5 53,8 327,9 88,6 150,0 40,5 230,1 61,0 2014 150,7 37,7 265,1 69,7 96,1 27,5 170,6 45,0 Tổng: Giá trị lượng mưa trung bình theo tháng (mm) Max: Giá trị lượng mưa max theo ngày (mm) 350,0 300,0 250,0 200,0 289,1 276,4 275,4 268,8 273,6 267,4 264,2 263,7 259,6 245,1 244,7 241,7 230,1 225,1 219,5 219,3 218,2 215,7 218,0 217,1 215,8 202,8 202,3 170,6 150,0 Max 100,0 50,0 Tổng 76,1 74,7 74,2 64,4 67,067,656,971,264,066,7 62,659,9 61,4 58,767,7 68,861,8 69,257,261,0 53,4 48,5 48,1 45,0 1991 1992 1993 1994 1995 1996 1997 1998 1999 2000 2001 2002 2003 2004 2005 2006 2007 2008 2009 2010 2011 2012 2013 2014 0,0 129 Hình 4.4: Diễn biến thay đổi yếu tố mƣa khu vực nghiên cứu Biểu đồ cho ta thấy liên quan mật thiết thay đổi lƣợng mƣa thay đổi quy mô trƣợt đất Năm 1999, 2002, 2008, 2011 lƣợng mƣa tăng đột biến khu vực, thời điểm số lƣợng khối trƣợt gia tăng cách đáng kể, đặc biệt khối trƣợt lớn Nhƣ vậy, kết tính toán trọng số yếu tố liệu thống kê qua năm cho thấy lƣợng mƣa yếu tố chiếm tỷ trọng gây trƣợt lớn, biến đổi nhanh có tác động trực tiếp đến tai biến trƣợt đất khu vực nghiên cứu Mô hình cảnh báo nhanh tai biến trƣợt đất đƣợc xây dựng số liệu quan trắc biến đổi lƣợng mƣa tần suất xuất trƣợt đất khu vực nghiên cứu Đồ thị phân tán tần suất xuất khối trƣợt biến đổi lƣợng mƣa trạm quan trắc cho ta phƣơng trình hồi quy tuyến tính Phƣơng trình đƣợc sử dụng để cảnh báo ngƣỡng biển đổi lƣợng mƣa gây trƣợt cho khu vực ngƣỡng biến đổi lƣợng mƣa gây trƣợt ạt, quy mô lớn toàn khu vực nghiên cứu Từ đó, đƣa khuyến cáo để quyền địa phƣơng có giải pháp chủ động ứng phó, giảm thiểu thiệt hại trƣợt đất gây KẾT LUẬN Trên sở nghiên cứu đánh giá nguy trƣợt đất luận chứng hệ thống quan trắc phục vụ cảnh báo tai biến trƣợt đất cho khu vực Tây Nam tỉnh Hà Giang, rút số kết luận sau: Khu vực Tây Nam Hà Giang, trƣợt đất xảy với quy mô rộng khắp huyện Hoàng Su Phì, Xín Mần Quang Bình, nhiều khối trƣợt tích lớn (>100.000m3), mật độ tập trung cao ( trung bình 01 khối/3,27 km2) Khu vực nghiên cứu tồn kiểu trƣợt đất chủ yếu trƣợt chảy trƣợt hỗn hợp nằm xen kẽ quy luật phân bố 130 Trƣợt đất xuất tập trung vào mùa mƣa, đặc biệt tăng số lƣợng, quy mô tần xuất vào năm có lƣợng mƣa thay đổi đột biến (2002, 2006, 2008) Phƣơng pháp số thống kê tích hợp đa biến tính toán dự báo nguy trƣợt đất cấp khu vực phản ánh đƣợc tính đa dạng phức tạp yếu tố phát sinh, phát triển trƣợt đất, áp dụng tốt để dự báo nguy trƣợt đất cho khu vực Tây Nam tỉnh Hà Giang Nguy trƣợt đất khu vực Tây Nam tỉnh Hà Giang đƣợc định chủ yếu 11 yếu tố quan trọng với trọng số (vai trò) tham gia tƣơng ứng, bao gồm: Hiện trạng sử dụng đất: 12,9%; cao độ địa hình:12%; khoảng cách đến đƣờng giao thông: 11,7%; lƣợng mƣa: 9,5%; mật độ phân cắt sâu: 9,2%; địa chất: 7,9%; độ dốc: 7,7; hƣớng dốc: 7,4%; khoảng cách đến đứt gãy: 7,3%; chiều dày vỏ phong hóa: 7,1%; mật độ phân cắt ngang: 6,8%) Nguy trƣợt đất khu vực nghiên cứu đƣợc chia thành cấp: cao chiếm 15,78% diện tích toàn khu vực nghiên cứu, cao -35,59%, trung bình 33,37%, thấp- 10,74%, thấp - 4,52% Hệ thống quan trắc phục vụ cảnh báo tai biến trƣợt đất khu vực Tây Nam tỉnh Hà Giang đƣợc xây dựng sở đồ đánh giá nguy trƣợt đất đồ giá trị trọng số phân lớp yếu tố, bao gồm 03 phụ hệ thống quan trắc thành phần: - Hệ thống quan trắc phân bố không gian khối trƣợt toàn khu vực nghiên cứu phƣơng pháp phân tích ảnh vệ tinh kết hợp thị sát trƣờng: phục vụ xây dựng đồ trạng trƣợt đất - Hệ thống lƣới trạm quan trắc đo vẽ đặc điểm khối trƣợt yếu tố điều kiện chủ yếu: phục vụ thống kê xác định quan hệ trạng đặc điểm trƣợt đất với yếu tố điều kiện gây trƣợt làm sở xác định trọng số 12 yếu tố, bao gồm 07 tuyến, 50 trạm đo vẽ khảo sát - Hệ thống quan trắc yếu tố tác động biến đổi nhanh (nguyên nhân gây trƣợt): phục vụ cảnh báo tai biến trƣợt đất, bao gồm 13 trạm quan trắc tự 131 động lƣợng mƣa phân tích ảnh vệ tinh kết hợp thị sát trƣờng để xác định biến động sử dụng đất Dữ liệu thu đƣợc thông qua hệ thống quan trắc đƣợc sử dụng để phục vụ cảnh báo tai biến trƣợt đất khu vực nghiên cứu thông qua mô hình cảnh báo sau: - Điều chỉnh trọng số yếu tố điều kiện, nguyên nhân gây trƣợt - Phân vùng nguy trƣợt đất khu vực (sử dụng mô hình định lƣợng tiêu tích hợp yếu tố điều kiện nguyên nhân gây trƣợt) - Cảnh báo nhanh quy mô, cƣờng độ thể tích khối trƣợt xuất thông qua quan trắc yếu tố gây trƣợt chủ yếu - Cảnh báo, dự báo ngƣỡng phát triển tai biến trƣợt đất toàn khu vực nghiên cứu theo yếu tố có tỷ trọng gây trƣợt lớn, biến đổi nhanh Bản đồ dự báo nguy trƣợt đất luận chứng hệ thống quan trắc phục vụ cảnh báo tai biến trƣợt đất khu vực Tây Nam tỉnh Hà Giang có đầy đủ sở khoa học thực tiễn, làm để triển khai đề án quy hoạch, khai thác bền vững lãnh thổ đề xuất giải pháp giảm thiểu thiệt hại trƣợt đất gây khu vực vùng Tây Nam tỉnh Hà Giang 132 TÀI LIỆU THAM KHẢO 10 11 12 13 14 Công ty tƣ vấn lập Xí nghiệp Xây dựng Chuyển giao Công nghệ Môi trƣờng (2013), Báo cáo điều tra, đánh giá mức độ tác động biến đổi khí hậu đến lũ quét, lũ ống, trượt, sạt lở đất địa bàn huyện Yên Minh, Hoàng Su Phì, Xín Mần, Quang Bình-tỉnh Hà Giang xây dựng biện pháp ứng phó với biến đổi khí hậu Ragozin A.L (2000), Đánh giá quản lý rủi ro thiên nhiên, NXB ANKIN, Matxcova (Tiếng Nga) Đỗ Ngọc Ánh (2015), Nghiên cứu phân vùng cảnh báo trượt lở thị xã Mường Lay, tỉnh Điện Biên điều kiện vận hành hồ chứa thuỷ điện Sơn La đề xuất giải pháp ứng phó, giảm thiểu hậu quả, Báo cáo tổng kết đề tài cấp Bộ Viện Thủy điện lƣợng tái tạo, Hà Nội Hồ Việt Cƣờng (2011), Báo cáo lập đồ nguy lũ quét sạt lở đất địa bàn huyện Đắc Tô, Đắc Glei Tu Mơ rông tỉnh Kon Tum, Báo cáo tổng kết đề tài, CTCP Tƣ vấn thiết kế xây dựng thƣơng mại Thăng Long, Hà Nội Lomtadze V D (1979), Địa chất công trình, địa chất động lực công trình, NXB Đại học Trung học chuyên nghiệp, HN (Bản dịch tiếng Việt) Lê Mục Đích (2001), Kinh nghiệm phòng tránh kiểm soát tai biến địa chất, Nxb Xây dựng, Hà Nội (Dịch từ tiếng Trung Quốc) Sở KHCN Hà Giang (2003, 2004), Báo cáo Hiện trạng, quy hoạch sử dụng đất đến năm 2010, (bản đồ tỷ lệ: 1/10 000 kèm theo) xã thuộc huyện Yên Minh, Hoàng Su Phì Xín Mần, Hà Giang, Niên giám thống kê 2003; 2004 Trần Trọng Huệ (2010), Báo cáo kết KHCN đề tài “Nghiên cứu đánh giá dự báo chi tiết tượng trượt – lở xây dựng giải pháp phòng chống cho Thị trấn Cốc Pài, huyện Xín Mần, tỉnh Lào Cai, Báo cáo tổng kết đề tài cấp Nhà nƣớc mã số KC.08.33/06-10, Viện KH&CN VN, HN Lê Quốc Hùng (2014), Điều tra, đánh giá phân vùng cảnh báo nguy trượt lở đất đá vùng miền núi Việt Nam, Viện Khoa học Địa chất & Khoáng sản Nguyễn Xuân Huyên (2009), Đánh giá nguy đề xuất giải pháp phòng chống, phòng tránh tai biến trượt lở đất thành phố Đà Nẵng (20072008), Viện Hàn Lâm KHCN Việt Nam, Hà Nội Đặng Quang Khang (2011), Đánh giá rủi ro tai biến trượt lở (lấy ví dụ khu vực đèo gió, huyện Ngân Sơn, tỉnh Bắc Cạn), Luận văn thạc sĩ Địa chất học Trƣờng Đại học KHTN – ĐHQGHN, Hà Nội Trần Mạnh Liểu (2007), "Cơ sở tiếp cận thệ thống đánh giá dự báo tổng hợp tai biến địa chất", Tạp chí Địa kỹ thuật, số 2/2007 Trần Mạnh Liểu (2007), "Phƣơng pháp phân vùng dự báo khả phát triển tai biến địa chất theo tiêu tích hợp yếu tố phát triển tai biến", Tạp chí xây dựng số 9/2007 Trần Mạnh Liểu (2013), Một vài phương pháp đánh giá định tính định lượng vai trò yếu tố hình thành phát triển tai biến địa chất, Trung tâm Nghiên cứu đô thị - ĐHQGHN, HN 133 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 Nguyễn Đức Lý (2011), Nghiên cứu trình dịch chuyển trọng lực đất đá sườn dốc, mái dốc tuyến đường giao thông Tây Quảng Bình đề xuất giải pháp phòng chống, Luận án TS, HN Trần Mạnh Liểu nnk (2011), Báo cáo đề tài “Nghiên cứu đánh giá rủi ro Dự báo nguy trượt lở thị xã Bắc Kạn”, Trung tâm Nghiên cứu đô thị ĐHQGHN, HN Nguyễn Thị Thanh Nhàn (2014), Nghiên cứu tượng dịch chuyển đất đá sườn dốc vùng đồi núi Quảng Trị - Thừa Thiên Huế, đề xuất phương pháp dự báo phòng chống phù hợp, Luận án Tiến s , Trƣờng Đại học Mỏ địa chất, Hà Nội Xoigu S K Osipova V I (2003), Những tai biến thiên nhiên Nga (Tiếng Nga) (Tập 1, 2, 3, 4, 5, 6), Nxb firma “KPYK”, Moskva Xôigu S K Osipova V I (2003), Những tai biến thiên nhiên Nga (Tiếng Nga) (Tập 1, 2, 3, 4, 5, 6), NXB Firma “KPYK”, Moskva Ngô Thị Phƣợng (2000), Xác định vùng có nguy trượt lở đe doạ trực tiếp đến môi trường sống tính mạng nhân dân sau trận lũ quét đêm 7/6/2001 Trùng Khánh, Hạ Lang (Cao Bằng) làm sở khoa học cho phòng tránh hữu hiệu, Viện Địa chất, Viện Hàn lâm KH&CN Việt Nam, Hà Nội Hệ đồ Địa chất môi trƣờng Trung Quốc (1992), Hƣớng dẫn sử dụng đồ phân loại, phân bố trƣợt lở đá lở Trung Quốc, chủ biên, (Bản dịch Viện địa chất), NXB Bản đồ, Trung Quốc, Bắc Kinh Đinh Văn Toàn (2006), Nghiên cứu phân vùng dự báo trượt lở đất, lũ quéy tỉnh Hoà Bình giải pháp phòng tránh giảm nhẹ thiệt hại (2005-2006) , Viện Địa chất, Viện Hàn lâm KHCN Việt Nam, Hà Nội Cotlov F V (1978), Những thay đổi địa chất môi trường ảnh hưởng hoạt động người (Tiếng Nga)Nxb “Nedra”, Moskva Osipov V.I (1999), Các trình tai biến nguồn gốc nội sinh, NXB Geos, Matxcova (Tiếng Nga) Nguyễn Trọng Yêm (2011), Những đặc điểm, nguyên nhân T-L, LQ-LBĐ Yên Minh, Hoàng Su Phì, Xín Mần, đề xuất giải pháp phòng chống thích hợp cho địa phương, Đề tài NCKH cấp Bộ KC-08-01 & KC-0801BS, Hà Nội Nguyễn Trọng Yêm (2011), Những đặc điểm, nguyên nhân T-L, LQ-LBĐ Yên Minh, Hoàng Su Phì, Xín Mần, đề xuất giải pháp phòng chống thích hợp cho địa phương, Đề tài NCKH cấp Bộ KC.08.01 & KC.08.01BS Viện Địa chất, Viện Hàn lâm KHCN Việt Nam, Hà Nội Fabbri A.G Chung C.F, Van Westen C.J (1995), Multivariate regression analysis for landslide hazard zonation In Geographical infomation systems in assessing natural hazards, sous la dir De A Carrara, F Guzzetti p.107-133 Netherland: Kluwer Academic Publishers Cruden D.M and Varnes D.J (1996), Landslide types and processes In A K Turner and R L Schuster (eds) Landslides - Investigation and Mitigation, National Academy Press, Transportation Reseach Board Special Report 247, Washington D.C pp 36-75 Varnes D.J (1978), Landslide types and processes in R L schuster and R J Krizek (eds), Landslide Analysis and control pp 11-33 Special Report 176 134 30 31 32 33 34 35 36 37 38 39 40 41 42 43 Ellott C Spiker and Paula L Gori (2000), National Landslide Hazards Mitigation Strategy, Open file report 00-450 Department of the interior US Geological Survey ltd (JRC) Japan Radio Co (2005), Natural disaster prevention system, Debris flow Monitoring system Arora M K Kanungo D P , Sarkar S and Gupta R P (2009), Landslide Susceptibility Zonation (LSZ) Mapping - A Review, Journal of South Asia Disaster Studies.,Vol No June 2009 Zhoua C.H Lana H.X , Wangb L.J , Zhangc H.Y , Lic R.H (2004), Landslide hazard spatial analysis and prediction using GIS in the Xiaojiang watershed, Yunnan, China United Nations (2001), Guidelines for reducing flood losses Krizek R.J (ed) Schucter R.L (1978), Landslides: Analysis and Control, NAS, NRC, Ja Special report N-176 P.R China The Ministry of Geology and Mineral Resources (1991), Geological hazards in China and their prevention and control, Geological Publishing house, Beijing, China UNEP (1998), Landslides and Mudflows, Vol 1&2 Unesco, Moscow Lomtadze V.D (1977), Engineering Geology Engineering Geodynamic, Nedra Publishing house,Leningrad, p.495 Jerusalimckaya E.N Bondaric G.K., Yarg L.A (1981), "Nguyên tắc tổ chức giám sát cấp khu vực hệ thống tự nhiên k thuật môi trƣờng địa chất ", Tạp chí Địa chất, số 3/2009 Trung tâm khí tƣợng thủy văn Quốc gia (2014), Thống kê lượng mưa trạm đo khí tượng địa phận tỉnh Hà Giang Nguyễn Xuân Khang (2016), Phát triển công nghệ đánh giá rủi ro trượt đất dọc theo tuyến giao thông Việt Nam, Dự án hợp tác k thuật JICA - JST, Viện Khoa học Giao thông Vận tải - Bộ Giao thông Voogd.H (1983), " Multicriteria Evaluation for Urban and Regional Planning", University of Groningen, The Netherlands Cees van Westen (1997), "statistical landslide hazard analysis", ITC Publication, Eschede, Netherland, tr 73-84 135 ... hội khu vực Tây Nam tỉnh Hà Giang Chƣơng 3: Đánh giá nguy trƣợt đất khu vực Tây Nam tỉnh Hà Giang Chƣơng 4: Luận chứng hệ thống quan trắc phục vụ cảnh báo tai biến trƣợt đất khu vực Tây Nam tỉnh. .. sinh khu vực nghiên cứu tác động đến trình trƣợt đất - Đánh giá nguy trƣợt đất khu vực Tây Nam tỉnh Hà Giang - Luận chứng hệ thống quan trắc phục vụ cảnh báo tai biến trƣợt đất khu vực nghiên cứu: ... 108 LUẬN CHỨNG HỆ THỐNG QUAN TRẮC PHỤC VỤ CẢNH BÁO TAI BIẾN TRƢỢT ĐẤT KHU VỰC NGHIÊN CỨU 108 4.1 Tổng quan hệ thống quan trắc phục vụ cảnh báo tai biến trƣợt đất 108 4.1.1 Phân cấp hệ thống