TẠP CHÍ KHOA HỌC, Đại học Huế, Tập 74B, Số 5, (2012), 143-155 143 XÂY DỰNG BẢN ĐỒ NGUY CƠ TRƯỢT LỞ ĐẤT TỈNH QUẢNG TRỊ BẰNG PHƯƠNG PHÁP TÍCH HỢP MÔ HÌNH PHÂN TÍCH THỨ BẬC (AHP) VÀO GIS Nguyễn Thám 1 , Nguyễn Đăng Độ 1 , Uông Đình Khanh 2 1 Trường Đại học Sư phạm, Đại học Huế 2 Viện Địa lý, Viện KHCN Việt Nam Tóm tắt. Quảng Trị là một tỉnh có diện tích không lớn, nhưng với địa hình đồi núi chiếm trên 70% diện tích, lượng mưa lớn và tập trung, lớp phủ rừng bị tàn phá do chiến tranh cùng với các hoạt động dân sinh đã làm cho quá trình trượt lở đất xảy ra ngày càng tăng về quy mô và cường độ. Phương pháp tích hợp mô hình phân tích thứ bậc (AHP) vào GIS để thành lập bản đồ nguy cơ trượt lở đất (TLĐ) là cách tiếp cận mang tính định lượng và có độ tin cậy cao. Bài báo sử dụng mô hình AHP để xác định trọng số của các yếu tố ảnh hưởng đến TLĐ kết hợp với chức năng phân tích không gian của GIS đã xây dựng được bản đồ nguy cơ TLĐ tỉnh Quảng Trị tỷ lệ 1/100.000. 1. Mở đầu Quảng Trị là một tỉnh nghèo, có diện tích 4.746,4km 2 , với trên 70% diện tích là đồi núi, mức độ chia cắt sâu và độ dốc của địa hình lớn, mạng lưới sông suối dày đặc, khoảng 0,8 - 1 km/km 2 , lòng sông hẹp, nhiều thác ghềnh, sông ngắn và dốc. Sự tương tác giữa hoàn lưu gió mùa với đặc điểm hệ thống sơn văn đã tạo cho Quảng Trị một lượng mưa dồi dào, trung bình 2000 - 2700mm/năm, mưa tập trung theo mùa (mùa mưa chiếm 80 - 90% lượng mưa cả năm), kết hợp với đặc điểm địa chất phức tạp, mức độ biến động lớn của thảm phủ thực vật đã làm cho hiện tượng trượt lở đất xảy ra với cường độ, quy mô và tần suất ngày càng cao, gây thiệt hại lớn về người và tài sản của nhân dân, nhất là vùng miền núi. Ngày nay, việc áp dụng công nghệ GIS trong nghiên cứu TLĐ đang trở nên phổ biến và chiếm ưu thế. Các dữ liệu không gian liên quan có thể được mô hình hóa như DEM, lượng mưa, phân cắt sâu, phân cắt ngang, thạch học, hiện trạng sử dụng đất, dòng chảy sông suối, kết hợp với việc xác định trọng số cho từng yếu tố ảnh hưởng đến TLĐ bằng phương pháp mô hình phân tích thứ bậc (AHP) để thành lập bản đồ nguy cơ TLĐ đã cho kết quả đáng tin cậy. 2. Cơ sở lý thuyết của phương pháp 2.1. Lựa chọn các nhân tố tác động gây nên quá trình trượt lở đất Chúng ta biết rằng, có rất nhiều nhân tố tác động lên môi trường gây nên hiện tượng trượt lở đất. Căn cứ vào điều kiện cụ thể của lãnh thổ nghiên cứu, chúng tôi chỉ 144 Xây dựng bản đồ nguy cơ trượt lở đất tỉnh quảng trị… lựa chọn 7 nhân tố đặc trưng, mang tính tiên quyết và bao trùm lên các nhân tố khác để xây dựng bản đồ nguy cơ TLĐ, bao gồm: độ dốc, lượng mưa trung bình năm, địa chất (thạch học), mật độ đứt gãy, phân cắt sâu, phân cắt ngang và hiện trạng sử dụng đất, còn các nhân tố khác chỉ được phân tích định tính và mang tính chất tham khảo. 2.2. Sử dụng AHP để xác định trọng số của các nhân tố ảnh hưởng đến trượt lở đất Có rất nhiều nhân tố tác động đến quá trình TLĐ, tuy nhiên vai trò của chúng là không hoàn toàn giống nhau. Vì vậy, việc xác định trọng số cho mỗi nhân tố này là rất cần thiết. Phương pháp phân tích thứ bậc (AHP) (hay còn gọi là phương pháp mô hình trọng số) là một phương pháp bán định lượng. Nội dung của phương pháp bao gồm việc xây dựng một hệ thống các cặp ma trận so sánh giữa các yếu tố khác nhau cho trượt lở đất. Cách tiếp cận này có thể được mô tả như là sự phân bậc tầm quan trọng của các nhân tố gây nên TLĐ, mỗi nhân tố được so sánh với các nhân tố khác để xác định tầm quan trọng của chúng đối với TLĐ [6], [7] (Bảng 1). Bảng 1. Mức độ tác động của các nhân tố phát sinh trượt lở đất Nhân tố Độ dốc Lượng mưa trung bình năm Địa chất (thạch học) Phân cắt sâu Phân cắt ngang Mật độ đứt gãy Hiện trạng sử dụng đất Điểm 9 5 4 4 3 2 1 Việc đánh giá mức độ nhạy cảm của các nhân tố trên thang điểm biểu thị sự ưu tiên của chúng một cách thích đáng đối với quá trình TLĐ. Trong đó, độ dốc sườn là nhân tố quan trọng bậc nhất ảnh hưởng đến TLĐ, do đó được cho 9 điểm. Các nhân tố còn lại, ứng với mức độ nhạy cảm đến khả năng gây ra tai biến TLĐ khác nhau được trình bày trên bảng 1. Trọng số của các nhân tố ảnh hưởng được xác định thông qua việc lập ma trận so sánh tương quan giữa các nhân tố và tính được trọng số tương ứng của từng nhân tố theo bảng 2. Bảng 2 . Ma trận xác định trọng số của các nhân tố ảnh hưởng đến TLĐ Các nhân tố Độ dốc Lượng mưa trung bình năm Mật độ đứt gãy Địa chất (thạch học) Phân cắt sâu Phân cắt ngang Hiện trạng sử dụng đất Trọng số (A) (B) (C) (D) (E) (F) (G) (A) 0.326 0.335 0.333 0.332 0.309 0.333 0.333 0.328 NGUYỄN THÁM, NGUYỄN ĐĂNG ĐỘ, UÔNG ĐÌNH KHANH 145 (B) 0.261 0.259 0.259 0.258 0.240 0.259 0.259 0.256 (C) 0.196 0.184 0.185 0.188 0.175 0.185 0.185 0.185 (D) 0.107 0.111 0.111 0.110 0.103 0.111 0.111 0.109 € 0.035 0.036 0.037 0.036 0.103 0.037 0.037 0.042 (F) 0.035 0.036 0.037 0.036 0.034 0.037 0.037 0.036 (G) 0.035 0.036 0.037 0.036 0.034 0.037 0.037 0.036 2.3. Phương pháp tích hợp AHP vào GIS để xây dựng bản đồ nguy cơ trượt lở đất Hệ thống tin địa lý (GIS) cho phép xây dựng các phân tích không gian, quản lý, tích hợp và chồng ghép các lớp thông tin. Mô hình phân tích thứ bậc sẽ hỗ trợ cho hệ thống thông tin địa lý, tổng hợp các thông tin, gán các trọng số phù hợp nhất cho các yếu tố đã được lựa chọn. Sau khi đã phân cấp và tính trọng số của các chỉ tiêu thì việc tích hợp chúng sẽ cho ta chỉ số nhạy cảm trượt lở đất. Để định lượng hóa mức độ nhạy cảm phản ánh nguy cơ trượt lở đất, bài báo tiến hành tích hợp các chỉ tiêu theo công thức (1) (theo Patrono, et al., 1995):    n j WjXijLSI 1 (1) Trong đó: LSI (Landslide Susceptibility Index): là chỉ số nhạy cảm trượt lở đất Wj: là trọng số của nhân tố thứ j Xij: là điểm số của lớp thứ i trong nhân tố gây trượt j Việc tích hợp AHP vào GIS thông qua công thức (1) được chúng tôi liên kết và tính toán bằng công cụ Raster Calculator của phần mềm ArcGis 9.3. Hình 1. Quy trình xây dựng bản đồ nguy cơ trượt lở đất tỉnh Quảng Trị Cơ sở dữ liệu các bản đồ chuyên đề: Bản đồ địa chất, bản đồ đẳng trị mưa, bản đồ hiện trạng sử dụng đất Cơ sở dữ liệu bản đồ địa hình Chuyên gia BĐ lượng mưa TBN BĐ thạch học BĐ mật độ đứt gãy Bản đồ HTSDĐ Mô hình DEM BĐ độ dốc BĐ phân cắt sâu BĐ phân cắt ngang Mô hình trọng số GTTS theo lượng mưa TBN GTTS theo thạch học GTTS theo mật độ đứt gãy GTTS theo HTSDĐ GTTS theo độ dốc GTTS theo phân cắt sâu GTTS theo phân cắt ngang Tích hợp GIS Bản đồ nguy cơ TLĐ 146 Xây dựng bản đồ nguy cơ trượt lở đất tỉnh quảng trị… 3. Đánh giá vai trò của các nhân tố chủ yếu tác động đến quá trình trượt lở đất 3.1. Nhân tố độ dốc Độ dốc địa hình có vai trò quyết định tới sự hình thành và phát triển trượt lở. Khi góc dốc càng lớn thì mức độ ổn định của sườn càng nhỏ và ngược lại khi độ dốc bằng không thì sẽ không có trượt lở. Thống kê độ dốc khu vực tỉnh Quảng Trị cho thấy giá trị độ dốc biến thiên từ 0-58 o . Địa hình có độ dốc > 15 o chiếm đến 29.74%, trong đó những khu vực có độ dốc từ 25-35 o chiếm 8.02% và khu vực có độ dốc > 35 o chỉ chiếm 0.78%. Bản đồ độ dốc sau khi được thành lập và thống kê được chúng tôi chia thành 5 cấp giá trị tương ứng với mức độ tác động của nó đến tai biến TLĐ khác nhau (Bảng 3). Bảng 3. Phân cấp ảnh hưởng của nhân tố độ dốc đến quá trình TLĐ Cấp ảnh hưởng Nguy cơ TLĐ Độ dốc ( o ) Diện tích (ha) Diện tích (%) Cấp 1 Rất thấp < 8 248.764 52,42 Cấp 2 Thấp 8-15 84.655 17,84 Cấp 3 Trung bình 15-25 99.396 20,94 Cấp 4 Cao 25-35 38.049 8,02 Cấp 5 Rất cao > 35 3.706 0,78 Toàn tỉnh 474.570 100 3.2. Nhân tố lượng mưa trung bình năm Lượng mưa là thông số rất quan trọng quyết định đến quá trình TLĐ. Cường độ TLĐ gia tăng tỷ lệ thuận với lượng mưa, đặc biệt là với cường độ mưa trận. Lượng mưa trung bình năm của tỉnh Quảng Trị là khá lớn, trong khu vực toàn tỉnh giá trị này dao động từ 2151,87-2679,99 mm/năm. Bản đồ lượng mưa trung bình năm được thành lập bằng việc nội suy giá trị lượng mưa trung bình năm của các trạm đo mưa và đường đẳng trị mưa trung bình năm. Sau khi thống kê các giá trị nội suy cho toàn tỉnh, chúng tôi chia thành 5 cấp giá trị tương ứng với mức độ tác động của nó đến tai biến trượt lở đất khác nhau (Bảng 4). Bảng 4. Phân cấp ảnh hưởng của nhân tố lượng mưa trung bình năm đến quá trình TLĐ Cấp ảnh hưởng Nguy cơ TLĐ Lượng mưa TBN (mm/năm) Diện tích (ha) Diện tích (%) Cấp 1 Rất thấp < 2.200 5.826 1,23 Cấp 2 Thấp 2.200-2.300 4.687 0,99 Cấp 3 Trung bình 2.300-2.400 96.697 20,38 NGUYỄN THÁM, NGUYỄN ĐĂNG ĐỘ, UÔNG ĐÌNH KHANH 147 Cấp 4 Cao 2.400-2.500 123.844 26,10 Cấp 5 Rất cao > 2.500 243.516 51,31 Toàn tỉnh 474.570 100 3.3. Nhân tố địa chất (thạch học) Ảnh hưởng của điều kiện địa chất, kiến tạo được coi là một nhân tố cơ bản gây ra quá trình TLĐ, đặc biệt thành phần thạch học là một trong những nhân tố quan trọng nhất ảnh hưởng đến sự ổn định của sườn. Các đá có độ bền thấp dễ có xu hướng phong hoá thành các vật liệu kém bền vững. Trong nhóm này bao gồm sét, đá phiến, một số tuf núi lửa và các đá có chứa các khoáng vật dạng tấm yếu như mica (các đá biến chất). Trên cơ sở bản đồ địa chất tỉnh Quảng Trị tỉ lệ 1:100.000 do Lê Đức An thành lập [1], trong đó các thuộc tính về nguồn gốc, tuổi, thành phần thạch học và đặc biệt là thành phần vật chất của các thành tạo địa chất đã được chúng tôi tích hợp để xây dựng bản đồ phân cấp ảnh hưởng của thành phần thạch học đến TLĐ với 5 cấp như sau (Bảng 5). Bảng 5. Phân cấp ảnh hưởng của nhân tố địa chất (thạch học) đến quá trình TLĐ Cấp ảnh hưởng Nguy cơ TLĐ Địa chất (thạch học) Diện tích (ha) Diện tích (%) Cấp 1 Rất thấp Nhóm carbonat 3.523 0,74 Cấp 2 Thấp Nhóm đá phun trào 42.460 8,95 Cấp 3 Trung bình Nhóm đá xâm nhập 69.356 14,61 Cấp 4 Cao Nhóm đá trầm tích lục nguyên và biến chất 248.264 52,31 Cấp 5 Rất cao Nhóm trầm tích bở rời 110.967 23,38 Toàn tỉnh 474.570 100 3.4. Nhân tố mật độ đứt gãy Mật độ đứt gãy là một nhân tố quan trọng gây nên tai biến TLĐ, chúng ảnh hưởng trực tiếp đến mật độ và quy mô điểm trượt. Khi phân tích các khối trượt đơn lẻ, nhân tố này được đặc trưng bằng mức độ dập vỡ, nứt nẻ của đất đá. Nhưng khi xem xét trên phạm vi một lãnh thổ thì nhân tố này có ý nghĩa lớn đối với quá trình TLĐ. Trong cùng một loại đất đá thì TLĐ dễ phát sinh ở những đới dập vỡ, nứt nẻ vì ở đây đất đá thường dễ bị phong hóa, dễ bị bão hoà nước nên có độ bền chống cắt thấp. Mức độ dập vỡ, nứt nẻ của đất đá thường là do các quá trình phá huỷ kiến tạo như các đứt gãy kiến tạo, các đới tiếp xúc,… và các quá trình phong hóa. Việc xây dựng chỉ tiêu mật độ đứt gãy được nội suy từ bản đồ đứt gãy kiến tạo 148 Xây dựng bản đồ nguy cơ trượt lở đất tỉnh quảng trị… của khu vực với sự hỗ trợ của công cụ GIS. Kết quả nội suy được thể hiện theo 5 cấp giá trị tương ứng với 5 cấp ảnh hưởng của nó đến quá trình TLĐ (Bảng 6). Bảng 6. Phân cấp ảnh hưởng của nhân tố mật độ đứt gãy đến quá trình TLĐ Cấp ảnh hưởng Nguy cơ TLĐ Mật độ đứt gãy (km/km 2 ) Diện tích (ha) Diện tích (%) Cấp 1 Rất thấp < 1 156.287 32,93 Cấp 2 Thấp 1-2 62.294 13,13 Cấp 3 Trung bình 2-3 62.465 13,16 Cấp 4 Cao 3-4 50.299 10,60 Cấp 5 Rất cao > 4 143.226 30,18 Toàn tỉnh 474.570 100 3.5. Nhân tố phân cắt sâu Độ cao tương đối của địa hình là biên độ dao động về độ cao của bề mặt đất, nghĩa là độ chênh cao tương đối giữa đỉnh các địa hình dương với đáy của các địa hình âm gần nhất. Nhân tố này thể hiện vai trò thế năng của địa hình. Khi độ cao tương đối càng lớn thì thế năng địa hình càng cao và ngược lại, điều này đã thúc đẩy quá trình dịch chuyển của đất đá xảy ra mạnh hơn và động năng va đập của đất đá thể hiện tính khốc liệt rõ nét hơn. Giá trị độ cao tương đối của địa hình thông thường được thể hiện qua chỉ số phân cắt sâu của nó. Khi tính toán trên GIS chỉ tiêu này được xác định bằng việc tính độ chênh cao địa hình (mét) tại mỗi ô lưới vuông có diện tích là 1km 2 trên bản đồ địa hình tỉ lệ 1:100.000. Cũng như các nhân tố khác quyết định đến khả năng gây ra tai biến TLĐ, lớp thông tin bản đồ phân cắt sâu địa hình được tính toán, thống kê và chia thành 5 cấp giá trị: < 10; 10-20; 20-100; 100-300 và > 300m/km 2 tương ứng với 5 cấp ảnh hưởng đến quá trình TLĐ (Bảng 7). Bảng 7. Phân cấp ảnh hưởng của nhân tố phân cắt sâu đến quá trình TLĐ Cấp ảnh hưởng Nguy cơ TLĐ Phân cắt sâu (m/km 2 ) Diện tích (ha) Diện tích (%) Cấp 1 Rất thấp < 10 225.211 47,46 Cấp 2 Thấp 10-20 96.950 20,43 Cấp 3 Trung bình 20-100 133.859 28,21 Cấp 4 Cao 100-300 17.938 3,78 Cấp 5 Rất cao > 300 612 0,13 Toàn tỉnh 474,570 100 NGUYỄN THÁM, NGUYỄN ĐĂNG ĐỘ, UÔNG ĐÌNH KHANH 149 3.6. Nhân tố phân cắt ngang Hệ thống sông suối là bức tranh thể hiện kết quả sự phân cắt địa hình dưới tác động của dòng chảy. Nước trên bề mặt địa hình rất nhạy cảm và linh động với sự thay đổi của địa hình. Vì thế, nó cũng phản ánh phần nào chế độ kiến tạo của khu vực mà cụ thể là nhiều hệ sông suối được hình thành từ các hệ thống đứt gãy. Mật độ phân cắt ngang thể hiện sự phân cắt theo chiều ngang của địa hình, là thông số xác định gián tiếp nguy cơ xảy ra trượt lở đất. Mật độ phân cắt ngang địa hình được hiểu là tổng độ dài tất cả các rãnh xâm thực, khe xói (dòng chảy tạm thời), sông suối (dòng chảy thường xuyên) trên một diện tích nhất định nào đó (thường là 1km 2 ). Mật độ chia cắt ngang được tính bằng công thức (2): L = S l  (km/km 2 ) (2) Trong đó: L - là giá trị mật độ chia cắt ngang; l (km) - là tổng chiều dài các đường đáy được xác định trên một đơn vị diện tích S (km 2 ). Khu vực tỉnh Quảng Trị có mức độ phân cắt ngang phức tạp, nơi có giá trị cao không chỉ tập trung ở vùng đồi và đồng bằng ven biển mà còn ở các thung lũng lớn, các bồn trũng lớn phần trung tâm và phía Tây của tỉnh. Giá trị này thay đổi từ 0m/km 2 đến 4.86km/km 2 , được chia làm 5 cấp: < 1; 1-2; 2-3; 3-4 và > 4km/km 2 tương ứng với 5 cấp ảnh hưởng của nó đến quá trình trượt lở đất (Bảng 8). Bảng 8. Phân cấp ảnh hưởng của nhân tố phân cắt ngang đến quá trình TLĐ Cấp ảnh hưởng Nguy cơ TLĐ Phân cắt ngang (km/km 2 ) Diện tích (ha) Diện tích (%) Cấp 1 Rất thấp < 1 80.728 17,01 Cấp 2 Thấp 1-2 299.000 63,00 Cấp 3 Trung bình 2-3 78.232 16,48 Cấp 4 Cao 3-4 13.889 2,93 Cấp 5 Rất cao > 4 2.721 0,57 Toàn tỉnh 474.570 100 3.7. Nhân tố hiện trạng sử dụng đất Trong nghiên cứu, phân tích các nhân tố phát sinh TLĐ thì hiện trạng sử dụng đất được xem xét ở khía cạnh ảnh hưởng của lớp phủ thực vật. Các tính chất ăn sâu, ăn ngang của rễ, mật độ lớp phủ là những thông số quan trọng trong đánh giá ảnh hưởng của lớp phủ thực vật đối với tai biến trượt lở. Thảm thực vật rừng rậm thường xanh 150 Xây dựng bản đồ nguy cơ trượt lở đất tỉnh quảng trị… thường giữ cho địa hình ổn định hơn các kiểu thảm thực vật khác. Tính ổn định của địa hình còn tỷ lệ thuận với mật độ che phủ của lớp phủ thực vật. Tuỳ thuộc vào kiểu rừng và mức độ che phủ rừng mà tác động của chúng đối với quá trình TLĐ khác nhau. Hoạt động chặt phá, đốt rừng làm nương rẫy càng tăng nguy cơ TLĐ. Hiện trạng sử dụng đất của tỉnh Quảng Trị được chia thành 5 cấp ảnh hưởng đến quá trình TLĐ (Bảng 9). Bảng 9. Phân cấp ảnh hưởng của nhân tố hiện trạng sử dụng đất đến quá trình TLĐ Cấp ảnh hưởng Nguy cơ TLĐ Hiện trạng sử dụng đất Diện tích (ha) Diện tích (%) Cấp 1 Rất thấp Đất mặt nước (sông, hồ,…) 8.408 1,77 Cấp 2 Thấp Đất rừng trồng, rừng tự nhiên 72.126 15,20 Cấp 3 Trung bình Đất trồng cây công nghiệp lâu năm 10.094 2,13 Cấp 4 Cao Đất nông nghiệp, đồng cỏ, trảng cây bụi 327.612 69,03 Cấp 5 Rất cao Đất thổ cư, đất màu, bãi cát 56.331 11,87 Toàn tỉnh 474.570 100 4. Xây dựng bản đồ nguy cơ tai biến trượt lở đất tỉnh Quảng Trị Bản đồ nguy cơ trượt lở đất được xây dựng trên cơ sở của các phép phân tích không gian trong phần mềm ArcGIS. Các bản đồ nhân tố thành phần sau khi được phân cấp ảnh hưởng đến trượt lở đất, xác định trọng số tương ứng, được tích hợp tuyến tính theo công thức (3): LSI = 0,328*A + 0,256*B + 0,185*C + 0,109*D + 0,042*E + 0,036*F + 0,036*G (3) Trong đó: LSI: là chỉ số nhạy cảm trượt lở đất A: Nhân tố độ dốc B: Nhân tố lượng mưa trung bình năm C: Nhân tố mật độ đứt gãy D: Nhân tố địa chất (thạch học) E: Nhân tố phân cắt sâu F: Nhân tố phân cắt ngang G: Nhân tố hiện trạng sử dụng đất Kết quả thu được bản đồ nguy cơ TLĐ với các giá trị khác nhau trên mỗi pixel. Vì vậy, ở dạng nguyên thủy nó chưa đặc trưng cho một bản đồ cảnh báo nguy cơ TLĐ. Để hình thành bản đồ cảnh báo nguy cơ TLĐ chúng tôi tiến hành phân cấp lại thành 5 cấp nguy cơ tương ứng: rất thấp, thấp, trung bình, cao và rất cao. Ngưỡng giá trị để NGUYỄN THÁM, NGUYỄN ĐĂNG ĐỘ, UÔNG ĐÌNH KHANH 151 phân cấp bản đồ nguy cơ TLĐ được lựa chọn sau khi thực hiện xử lý thống kê giá trị trong phần mềm ArcGIS, kết quả cho ra đường cong tích lũy xác suất với các thông số như sau: Giá trị điểm nhỏ nhất (LSI min ) = 1,426 Giá trị điểm lớn nhất (LSI max ) = 8,556 Khoảng cách điểm giữa các cấp được xác định theo công thức: n LSILSI TLĐ minmax   = 426,1 5 426,1556,8   Kết quả phân cấp bản đồ nguy cơ TLĐ giá trị số được phân ra 5 cấp nguy cơ trượt lở đất với diện tích các cấp tương ứng cho toàn vùng nghiên cứu thể hiện tại bảng 10. Bảng 10. Phân cấp diện tích lãnh thổ Quảng Trị theo nguy cơ trượt lở đất Cấp ảnh hưởng Nguy cơ TLĐ Khoảng cách điểm Diện tích (ha) Diện tích (%) Cấp 1 Rất thấp 1,426 – 2,852 71.204,36 15,00 Cấp 2 Thấp 2,853 – 4,279 158.653,19 33,43 Cấp 3 Trung bình 4,280 – 5,705 128.167,60 27,01 Cấp 4 Cao 5,706 – 7,131 108.427,38 22,85 Cấp 5 Rất cao 7,132 – 8,556 8.117,47 1,71 Toàn tỉnh 474.570 100 5. Đánh giá nguy cơ tai biến trượt lở đất tỉnh Quảng Trị Kết quả thống kê trong bảng 10 cho thấy, khu vực có nguy cơ TLĐ cao và rất cao chiếm 24,56% (116.544,85ha) toàn tỉnh, trong đó cấp nguy cơ TLĐ rất cao chỉ chiếm 1.71% (8.117,47ha). Nguy cơ TLĐ ở cấp cao và rất cao tập trung lớn nhất trong huyện Đa Krông (chiếm 12,51% diện tích toàn tỉnh) và huyện Hướng Hoá (chiếm 7,93%). Các huyện Gio Linh (1,54%), huyện Vĩnh Linh (1,31%) có cấp nguy cơ TLĐ cao và rất cao chiếm tỷ lệ thấp. Các huyện còn lại giá trị TLĐ rất thấp, đặc biệt là thành phố Đông Hà và thị xã Quảng Trị thì gần như không xảy ra. Nguy cơ TLĐ ở cấp trung bình cũng tập trung chủ yếu ở huyện Đa Krông (chiếm 10,37% diện tích toàn tỉnh), huyện Hướng Hóa (8,76%); các huyện khác tỷ lệ này đều nhỏ. 152 Xây dựng bản đồ nguy cơ trượt lở đất tỉnh quảng trị… Hình 2. Bản đồ nguy cơ tai biến trượt lở đất tỉnh Quảng Trị Trên cơ sở bản đồ nguy cơ trượt TLĐ tỉnh Quảng Trị, bài báo tiến hành xem xét và thống kê trong từng đơn vị hành chính cấp xã của 4 huyện có nguy cơ tai biến TLĐ cao và rất cao như sau: - Huyện Vĩnh Linh: + Nguy cơ TLĐ ở cấp cao và rất cao: tập trung chủ yếu trong xã Vĩnh Ô (4.163,18ha), xã Vĩnh Hà (737,96ha) là các xã miền núi phía Tây của huyện; các xã Vĩnh Hiền (474,72ha), Vĩnh Hoà (292,74ha), Vĩnh Thạch (187,03ha), Vĩnh Kim (186,42ha) cũng có nguy cơ TLĐ cao tập trung ở phần đỉnh của khối basalt Vĩnh Linh. + Nguy cơ TLĐ ở cấp Trung bình: lớn nhất ở xã Vĩnh Hà (5.714,18ha), tiếp đến là xã Vĩnh Ô (4.154,43ha), xã Vĩnh Chấp (524,06ha). - Huyện Gio Linh: [...]... huyện Cam Lộ Như vậy, kết quả xây dựng bản đồ nguy cơ trượt lở đất đã phản ánh chính xác phạm vi phân bố theo không gian của các vùng có nguy cơ xảy ra trượt lở ở tỉnh Quảng Trị Các vùng được đánh giá có nguy cơ TLĐ cao và rất cao rơi vào những xã đã xảy ra trượt lở mạnh trong những năm gần đây 7 Kết luận - Tích hợp mô hình phân tích thứ bậc vào GIS để xây dựng bản đồ nguy cơ TLĐ là hướng tiếp cận mới... thông tin sử dụng để xây dựng bản đồ nguy cơ TLĐ tỉnh Quảng Trị tỷ lệ 1/100.000 bằng phương pháp sử dụng mô hình AHP chạy trong môi trường GIS khá đồng bộ cả về nội dung và tỷ lệ bản đồ cho phép đảm bảo độ tin cậy Kết quả đối sánh với hiện trạng TLĐ của tỉnh Quảng Trị năm 2009 - 2010 là minh chứng cho mức độ chính xác của bản đồ nguy cơ TLĐ tỉnh Quảng Trị - Bản đồ nguy cơ TLĐ tỉnh Quảng Trị tỷ lệ 1/100.000... ra trượt lở không đáng kể (Hình 3) 100 80 60 40 20 0 Số điểm Trượt lở Vĩnh Linh Gio Linh Đắkrông Hướng Hóa 15 7 74 81 Hình 3 Số điểm xảy ra trượt lở đất ở một số xã vùng đồi núi tỉnh Quảng Trị năm 2009-2010 Qua hình 3 cho thấy, trên địa bàn tỉnh Quảng Trị năm 2009 - 2010 có 177 điểm trượt lở Trong đó, huyện Hướng Hóa có số điểm trượt lở nhiều nhất (81 điểm, xảy ra ở xã Xây dựng bản đồ nguy cơ trượt lở. .. tai biến trượt lở đất, lũ bùn đá tỉnh Lào Cai, Luận án tiến sĩ, ĐHKHTN-ĐHQGHN, Hà Nội, 2009 NGUY N THÁM, NGUY N ĐĂNG ĐỘ, UÔNG ĐÌNH KHANH 155 [4] Nguy n Thám, Phan Văn Trung, Đánh giá hiện trạng và nguy cơ trượt lở đất dọc theo hành lang các tuyến giao thông quan trọng trên địa bàn tỉnh Quảng Trị bằng phương pháp đa chỉ tiêu, Tạp chí Khoa học Trường ĐHSP Hà Nội, số 56(3), (2011), 133-141 [5] Nguy n Thám,... (2.897,34ha),…, thấp nhất ở xã Mò Ó (786,45ha) + Nguy cơ TLĐ ở cấp trung bình: lớn nhất ở xã Tà Long (6.365,1ha), Hướng Hiệp (5.958,74ha), A Bung (5.336,28ha), Húc Nghì (5.157,68ha),…và thấp nhất ở xã A Ngo (1.393,07ha) 6 Hiện trạng trượt lở đất tỉnh Quảng Trị Kết quả khảo sát thực trạng trượt lở vùng đồi núi tỉnh Quảng Trị năm 2009 2010 của nhóm tác giả cho thấy, trượt lở đất xảy ra chủ yếu ở các huyện Vĩnh Linh,.. .NGUY N THÁM, NGUY N ĐĂNG ĐỘ, UÔNG ĐÌNH KHANH 153 + Nguy cơ TLĐ ở cấp cao và rất cao: tập trung chủ yếu ở xã miền núi Linh Thượng (6.622,82ha) và 2 xã vùng đồi Trung Sơn (133,54ha), Hải Thái (99,23ha) + Nguy cơ TLĐ ở cấp trung bình: lớn nhất ở xã Linh Thượng (7.888,17ha) và 3 xã vùng đồi Gio An (432,8ha), Vĩnh Trường (386,55ha), Trung Sơn (80,27ha) - Huyện Hướng Hóa: + Nguy cơ TLĐ ở cấp... cấp nguy cơ: rất thấp, thấp, trung bình, cao, rất cao và đánh giá cho các huyện, cũng như một số xã trong các huyện có nguy cơ TLĐ cao Kết quả đánh giá đã cho thấy: các huyện có nguy cơ xảy ra tai biến TLĐ cao và rất cao tập trung chủ yếu huyện Đa Krông, Hướng Hoá, Vĩnh Linh và Gio Linh Các xã có nguy cơ xảy ra tai biến TLĐ cao và rất cao: Hướng Sơn, Hướng Lập, Hướng Linh, Hướng Phùng, xã Húc, Tân Hợp. .. Phúc, Ba Nang của huyện Đa Krông; Vĩnh Ô, Vĩnh Hà của huyện Vĩnh Linh và Linh Thượng, Trung Sơn của huyện Gio Linh Kết quả này khá phù hợp với thực tế xảy ra TLĐ của tỉnh Quảng Trị hiện nay TÀI LIỆU THAM KHẢO [1] Lê Đức An, Uông Đình Khanh, Địa mạo và địa chất tỉnh Quảng Trị, Nxb Khoa học Tự nhiên và Công nghệ, Hà Nội, 2007 [2] Đào Đình Bắc và nnk., Vấn đề cảnh báo-dự báo tai biến thiên nhiên đảm bảo... xảy ra ở xã Xây dựng bản đồ nguy cơ trượt lở đất tỉnh quảng trị 154 Ba Tầng, Tân Thành, Húc, Hướng Phùng, Hướng Linh, Hướng Lập và Hướng Sơn) Huyện Đa Krông đã xảy ra 74 điểm trượt lở, tập trung ở xã Mò Ó, Ba Lòng, Triệu Nguy n, Tà Rụt, A Ngo, Ba Nang, Hải Phúc, Đắkrông, A Vao, Húc Nghì, A Bung, Hướng Hiệp và Tà Long Huyện Vĩnh Linh đã xảy ra 15 điểm trượt lở, chủ yếu là ở các xã Vĩnh Thủy, Vĩnh Khê,... (1.995,26ha), Tân Hợp (1.545,63ha) là các xã có địa hình núi trung bình và độ dốc cao + Nguy cơ TLĐ ở cấp trung bình: lớn nhất ở xã Hướng Lập (7.504,78ha), tiếp đến là các xã Hướng Sơn (6.908,48ha), Hướng Phùng (6.643,91ha), Hướng Linh (4.452,6ha), xã Húc (3.054,03ha), Ba Tầng (2.949,2ha), Hướng Lộc (2.448,95ha), Tân Thành (1.858,8ha), Tân Hợp (1.308,36ha) - Huyện Đa Krông: + Nguy cơ TLĐ ở cấp cao và . XÂY DỰNG BẢN ĐỒ NGUY CƠ TRƯỢT LỞ ĐẤT TỈNH QUẢNG TRỊ BẰNG PHƯƠNG PHÁP TÍCH HỢP MÔ HÌNH PHÂN TÍCH THỨ BẬC (AHP) VÀO GIS Nguy n Thám 1 , Nguy n. trình trượt lở đất xảy ra ngày càng tăng về quy mô và cường độ. Phương pháp tích hợp mô hình phân tích thứ bậc (AHP) vào GIS để thành lập bản đồ nguy cơ trượt