3. Ý NGHĨA KHOA HỌC VÀ THỰC TIỄN
2.2.5. xuất một số biện pháp PCCCR phù hợp trên địa bàn nghiên cứu
2.3. PHƯƠNG PHÁP NGHIÊN CỨU
2.3.1. Phương pháp thu thập số liệu
+ Phương pháp thu thập số liệu thứ cấp
- Số liệu điều kiện tự nhiên, kinh tế xã hội huyện A Lưới.
- Số liệu nhiệt độ, độ ẩm, lượng mưa từ các trạm quan trắc trong nhiều năm của tỉnh Thừa Thiên Huế.
- Các khu vực thường xảy ra cháy rừng trên địa bàn và những phương án PCCR hàng năm của tỉnh.
+ Phương pháp thu thập số liệu sơ cấp
- Thu thập các loại bản đồ số: Bản đồ hiện trạng, bản đồ địa hình, bản đồ 3 loại rừng, bản đồ phân bố lượng mưa, bản đồ hành chính.
- Tư liệu ảnh viễn thám (Landsat) năm 2005, 2010, 2015, ảnh vệ tinh Spot.
2.3.2. Phương pháp xử lý số liệu
Số liệu sau khi thu thập ở các nguồn khác nhau, sau đó tiến hành tổng hợp và đồng bộ hóa và đưa vào phần mềm ENVI, Ecognition Developper, ArcGis và Mapinfo để xử lý.
2.3.3. Phương pháp xây dựng Hệ thống hỗ trợ quyết định
Trình tự các bước phát triển Hệ thống hỗ trợ quyết định trong nghiên cứu lửa rừng được thể hiện qua sơ đồ sau.
đồ địa hình đồ hiện trạng hành chính vệ tinh điều tra kê thủy văn trạng rừng kê khí hậu
Raster Vector Xây dựng cơ sở dữ liệu nghiên cứu lửa rừng
CSDL không gian CSDL thuộc tính Hệ thống hỗ trợ quyết định (HTHTQD)
Giao diện người dùng
Chức năng Ứng dụng Thiết lập bản đồ chuyên đề Truy vấn và cập nhật dữ liệu Kết xuất dữ liệu Thu thập dữ liệu Lưu trữ dữ liệu Hiển thị lớp bản đồ Phân tích dữ liệu
Lập kế hoạch quản lý và giám sát lửa rừng
Xây dựng bản đồ phân cấp nguy cơ cháy rừng
Vị trí tọa độ các điểm cháy và hệ thống cảnh báo cháy rừng Đánh giá biến động lớp thảm thực vật
Thành lập bản đồ chuyên đề Nhập dữ liệu
Hình 2.1. Trình tự các bước phát triển HTHTQĐ quản lý lửa rừng
2.3.3.1. Phương pháp xây dựng bản đồ bằng phần mềm GIS
Phần mềm GIS cung cấp các chức năng và các công cụ cần thiết để lưu giữ, phân tích và hiển thị thông tin địa lý. Các thành phần chính trong phần mềm GIS là:
+ Công cụ nhập và thao tác trên các thông tin địa lý. + Hệ quản trị cơ sở dưc liệu (DBMS).
+ Công cụ hỗ trợ hỏi đáp phân tích và hiển thị địa lý.
+ Giao diện đồ họa người – máy để truy cập các công cụ dễ dàng.
Hình 2.2. Sơ đồ tổng quan làm bản đồ bằng GIS
Ứng dụng các chức năng cơ bản của GIS để xác định được các trọng số các yếu tố nghiên cứu và xây dựng được các lớp thuộc tính cho các lớp bản đồ: lớp bản đồ băng nhiệt, lớp bản đồ mật độ sông suối, lớp bản đồ độ cao, lớp độ dốc, lớp hướng dốc, bản đồ vùng trọng điểm cháy, bản đồ quy hoạch công trình PCCCR.
2.3.3.2. Phương pháp xây dựng bản đồ hiện trạng rừng từ ảnh vệ tinh Landsat
Tiến hành nghiên cứu giải đoán ảnh bằng các phần mềm, các công cụ đã lựa chọn để đưa ra được kết quả giải đoán tốt nhất. Nghiên cứu đã sử dụng phương pháp phân loại có kiểm định theo xác suất cực đại (Maximum likelihood) (Nguyễn Văn Lợi 2007) và kiểm định thực địa. Đánh giá kết quả giải đoán và chỉnh sửa bản đồ thành quả.
Số liêu ngoại nghiệp Bản đồ Tư liệu ảnh viễn thám
Bàn phím Bản đồ số hóa Máy quét Băng đĩa từ
GIS
BỘ XỬ LÝ TRUNG TÂM
Màn hình
Trình tự các bước được thể hiện qua sơ đồ sau:
Hình 2.3. Quy trình giải đoán ảnh viễn thám bằng phần mềm ENVI
Chồng lớp Nắn ảnh
Cắt ảnh theo khu vực nghiên cứu
Lựa chọn phương pháp phân loại ảnh phù hợp
- Phân loại và giải đoán ảnh - Xây dựng mẫu phân loại Bản đồ 3 loại
rừng
Lấy mẫu phân loại ở hiện trường
Kiểm tra mẫu phân loại
Phân loại ảnh Ảnh sau phân loại
năm 2005
Ảnh sau phân loại năm 2010
Ảnh sau phân loai năm 2015
Kiểm tra hiệu chỉnh sau phân loại
Bản đồ che phủ năm 2005
Bản đồ che phủ năm 2010 Bản đồ che phủ năm 2015
Bản đồ biến động hiện trạng tài nguyên rừng
2.3.3.3. Phương pháp xây dựng bản đồ chỉ số thực vật (NDVI)
Chỉ số thực vật (NDVI) là một đại lượng thay thế về số lượng thực vật và điều kiện sống. Chỉ số này liên kết với đặc điểm độ che phủ của thực vật như là sinh khối, chỉ số diện tích lá và phần trăm thực phủ.
Chỉ số thực vật (NDVI) được xác định dựa trên sự phản xạ khác nhau của thực vật thể hiện giữa kênh phổ khả biến và kênh phổ cận hồng ngoại, dùng để biểu thị mức độ tập trung của thực vật trên mặt đất.
Chỉ số thực vật được tính theo công thức sau: NDVI = (BIR - BR)/( BIR + BR) (2.1)
Trong đó: NDVI là chỉ số thực vật
BIR là kênh hồng ngoại BR là kênh đỏ
Giá trị của chỉ số thực vật là dãy số từ -1 đến +1. Nếu giá trị NDVI càng cao thì khu vực có độ che phủ thực vật cao. Nếu giá trị NDVI thấp thì khu vực đó độ che phủ thấp. Nếu giá trị NDVI âm cho thấy khu vực không có thực vật.
Từ tài liệu thu thập được như: Ảnh vệ tinh, bản đồ hiện trạng, bản đồ nền của huyện A Lưới, tiến hành xây dựng bản đồ chỉ số thực vật thông qua phần mềm Ecognition và Mapinfo. Trình tự các bước được thể hiện ở sơ đồ sau:
Hình 2.4. Sơ đồ các bước xây dựng bản đồ chỉ số thực vật.
Bản đồ cơ sở Tư liệu ảnh vệ tinh
Lớp bản đồ hiện trạng rừng Lớp bản đồ nền Lớp bản đồ segmentation Chồng lớp, nội suy BẢN ĐỒ CHỈ SỐ THỰC VẬT (NDVI)
2.3.3.3. Phương pháp xây dựng bản đồ khí hậu (nhiệt độ không khí, ẩm độ không khí, độ ẩm bão hòa, mưa, chỉ số khô hạn)
Khí hậu được xem là một nhân tố có ảnh hưởng lớn đến cháy rừng bao gồm các yếu tố nhiệt độ không khí, ẩm độ không khí, mưa, chỉ số khô hạn…Sau khi tiến hành thu thập số liệu các yếu tố trên, tiến hành xây dựng các bản đồ đơn tính và chồng xếp các lớp bản đồ dựa trên cơ sở phân tích không gian của phần mềm Arcgis. Trình tự các bước được thể hiện ở hình 2.5.
.
Hình 2.5. Sơ đồ các bước xây dựng bản đồ khí hậu.
2.3.3.4. Phương pháp xây dựng bản đồ thủy văn
Yếu tố thủy văn cũng ảnh hưởng tới cháy rừng. Mật độ sông ngòi dày thưa, dài ngắn đều ảnh hưởng khác nhau tới quá trình cháy rừng cũng như công tác PCCCR.
Lớp lượng mưa Lớp ẩm độ không khí Lớp độ ẩm bão hòa Lớp nhiệt độ Bản đồ hành chính Lớp chỉ số khô hạn Nội suy Lớp bản đồ nhiệt độ Lớp bản đồ chỉ số khô hạn Lớp bản đồ ranh giới Nội suy Chồng lớp, tính điểm BẢN ĐỒ KHÍ HẬU
Xây dựng bản đồ thủy văn dựa vào tài liệu quan trắc thu thập được và mô hình số độ cao (DEM) của huyện A Lưới, sử dụng công cụ trong phần mềm Arcgis hiệu chỉnh chồng xếp các lớp bản đồ đơn tính để tạo ra bản đồ thành quả. Trình tự các bước được thể hiện ở hình 2.6.
Hình 2.6. Sơ đồ các bước xây dựng bản đồ thủy văn.
2.3.3.5. Phương pháp xây dựng bản đồ tiếp cận khu rừng từ mạng lưới đường khu dân cư
Mạng lưới đường giao thông, khu dân cư gần rừng là những nhân tố quan trọng liên quan đến lửa rừng kể cả về mặt tích cực lẫn tiêu cực. Người dân đốt nương làm rẫy, khai thác rừng trái phép vô ý làm cháy rừng, đốt thực bì, hun khói bắt ong, đốt dọn vệ sinh ven đường, hệ thống đường giao thông kém đi lại khó khăn…đó là những mặt tiêu cực gây ra cháy rừng cũng như công tác phòng chống cháy. Ngược lại mạng lưới đường giao thông thuận lợi, nâng cao ý thức của người dân gần rừng là những mặt tích cực trong công tác PCCCR. Từ những số liệu thu thập được kết hợp với phần
Lớp ranh giới Lớp đường bình độ Lớp điểm độ cao Bản đồ thủy văn Lớp bề mặt DEM Lớp bản đồ lưu vực sông suối Lớp bản đồ mạng lưới sông suối Nội suy Chồng lớp, tính điểm BẢN ĐỒ THỦY VĂN
mềm chuyên dụng Arcgis, tiến hành hiệu chỉnh và chồng xếp các lớp bản đồ đơn tính để tạo ra bản đồ tiếp cận khu rừng, được thể hiện ở hình 2.7.
Hình 2.7. Sơ đồ các bước xây dựng bản đồ tiếp cận khu rừng.
2.4. PHƯƠNG PHÁP NGHIÊN CỨU CÁC ỨNG DỤNG CHÍNH CỦA HTHTQĐ GIS 2.4.1. Phương pháp đánh giá biến động diện tích của lớp thảm thực vật giai đoạn 2.4.1. Phương pháp đánh giá biến động diện tích của lớp thảm thực vật giai đoạn 2005 - 2010 và gia đoạn 2010 - 2015
Dữ liệu sau khi được thu thập từ các nguồn, phân tích, xử lý dữ liệu, tiến hành giải đoán, thành lập bản đồ hiện trạng các năm 2005,2010 và 2015 được thể hiện ở sơ đồ 2.8. Tiến hành chồng xếp bản đồ theo từng giai đoạn 2005-2010 và 2010-2015 Tổng kết và đánh giá kết quả. Toàn bộ phương pháp đánh giá biến động được thể hiện ở hình sau. Bản đồ hiện trạng rừng năm 2005 Bản đồ hiện trạng rừng năm 2010 Bản đồ hiện trạng rừng năm 2015 Chồng xếp, xử lý Kết quả biến động rừng theo từng giai đoạn
Bản đồ hành chính Tư liệu ảnh Landsat năm 2014
Lớp bản đồ phân bố dân cư Lớp bản đồ mạng lưới đường giao thông Lớp bản đồ ranh giới hành chính Lớp trạng thái rừng Chồng lớp, tính điểm BẢN ĐỒ TIẾP CẬN KHU RỪNG TỪ MẠNG LƯỚI ĐƯỜNG, KHU DÂN CƯ
Hình 2.8. Sơ đồ các bước xây dựng thống kê biến đổi hiện trạng rừng.
2.4.2. Phương pháp xây dựng bản đồ phân cấp nguy cơ cháy rừng
2.4.2.1. Ứng dụng AHP để xác định trọng số các nhân tố ảnh hưởng đến nguy cơ cháy
AHP dựa vào 3 nguyên tắc: (i) Phân tích vấn đề ra quyết định (thiết lập thứ bậc), (ii) Đánh giá so sánh các thành phần, (iii) Tổng hợp các độ ưu tiên.
Để tính toán mức độ ưu tiên giữa các chỉ tiêu, giả sử ta có Xn chỉ tiêu cần giả định thì một ma trận được giả thuyết như sau:
Bảng 2.1. Ma trận so sánh cặp đôi tầm quan trọng giữa các nhân tố
X1 X2 … Xn X1 a11 a12 … a1n X2 a21 a22 … a2n
… … … … …
Xn an1 an2 … ann
Để điền vào ma trận trên, sử dụng thang đánh giá từ 1- 9 của Saaty trong đánh giá các nhân tố ảnh hưởng đến nguy cơ cháy rừng ở vùng nghiên cứu như sau:
Bảng 2.2: Thang so sánh cặp nhân tố ảnh hưởng đến nguy cơ cháy rừng
Điểm giá trị Mức độ quan trọng
1 Hai nhân tố có tầm quan trọng như nhau 3 Nhân tố này quan trọng hơn nhân tố kia
5 Nhân tố này quan trọng hơn nhiều so với tiêu chí kia 7 Nhân tố này quan trọng hơn rất nhiều so với nhân tố kia 9 Nhân tố này tuyệt đối quan trọng hơn nhân tố kia
2,4,6,8 Khoảng trung gian ở giữa các mức độ trên
Sử dụng ma trận so sánh cặp đôi tính toán các trọng số của các nhân tố ảnh hưởng đến nguy cơ cháy rừng thông qua các đại lượng sau đây:
Wi=[(ai1)x(ai2)x...x(aim)]1/m (2.2) Nwi= I (2.3)
Dựa trên kết quả tính toán theo phương trình (2.2) và (2.3) ta có ma trận 2 như sau:
Bảng 2.3.Ma trận trọng số các nhân tố ảnh hưởng đến nguy cơ cháy rừng
X1 X2 X3 X4 Xn X1 w11 w12 w13 w14 w1n X2 w21 w22 w23 w24 w2n X3 w31 w32 w33 w34 w3n X4 w41 w42 w43 w44 w4n Xn Wn1 Wn2 Wn3 Wn4 Wnn
Ghi chú: X1, X2… Xn là các nhân tố ảnh hưởng đến nguy cơ cháy rừng.
Tác giả AHP( ông Saaty), đưa ra khái niệm “xung khắc”. Nếu tỉ số “ xung khắc” đạt mức <=10% thì đánh giá trong Bảng 2.3 chấp nhận được và ngược lại thì người đánh giá phải đánh giá và tính toán lại Bảng 2.1. Tính tỉ số nhất quán (consistency ratio – CR):
CR= (2.4)
CI: chỉ số nhất quán (Consistency Index)
RI: chỉ số ngẫu nhiên (Random Index). RI được xác định : RI = 1.845 x (1 - ) (2.5)
CI = (2.6)
λmax : giá trị riêng của ma trận so sánh n : số nhân tố
Sử dụng phương pháp phân tích thứ bậc AHP để xác định trọng số của các nhân tố ảnh hưởng chính đến phân vùng nhạy cảm cháy rừng bao gồm khí hậu, thủy văn, yếu tố xã hội và trạng thái thực vật. Năm nhân tố ảnh hưởng chính này được chia ra làm 8 nhân tố phụ và được xác định theo Bảng 2.4.
Bảng 2.4. Nhân tố chính và nhân tố phụ xác định STT Nhân tố chính Nhân tố phụ 1 Khí hậu Chỉ số khô hạn Nhiệt độ 2 Thảm thực vật rừng Trạng thái Chỉ số thực vật (NDVI) 3 Thủy văn
4 Điều kiện kinh tế - xã hội Khoảng cách từ khu dân cư
Khoảng cách từ đường giao thông
5 Địa hình Đai cao
Độ đốc
Theo phương pháp AHP tổng các giá trị của các nhân tố sẽ bằng 1 và tất cả các nhân tố sẽ có một trọng số tương ứng với mức độ quan trọng của nó, trọng số của các tiêu chí chính và phụ có giá trị từ 0 (kém quan trọng nhất) đến 1 (quan trọng nhất). Giá trị phù hợp tiềm năng nguy cơ cháy rừng được lựa chọn theo phương pháp cho điểm. Điểm xác định có giá trị từ thấp đến cao: giá trị 1 (ảnh hưởng thấp nhất), giá trị 9 (ảnh hưởng cao).
Bảng 2.5. Điểm phân cấp nguy cơ cháy rừng theo các lớp nhân tố đầu vào được lựa chọn
STT Lớp nhân tố đầu vào Khoảng chỉ tiêu Điểm Phân cấp rủi ro Nguy cơ cháy
1 Thảm thực vật rừng
Rừng trồng thông, keo, cây
bụi thường xanh 9 Cấp 5 Rất nguy hiểm Rừng phục hồi, lau lách, cỏ
tranh, cây bụi liền kề. 7 Cấp 4 Nguy hiểm Rừng thứ sinh nghèo, cây
STT Lớp nhân tố
đầu vào Khoảng chỉ tiêu Điểm
Phân cấp
rủi ro Nguy cơ cháy
Trảng cỏ thấp, thưa trên đất trơ sỏi đá, cây trồng nông nghiệp.
3 Cấp 2 Ít nguy hiểm
Rừng trung bình,rừng giàu
đất khác. 1 Cấp 1 Không nguy hiểm NDVI (< 0.1) 1 Cấp 1 Không nguy hiểm NDVI (0.1 - 0.2) 3 Cấp 2 Ít nguy hiểm
NDVI (0.2 - 0.3) 5 Cấp 3 Nguy hiểm trung bình NDVI (0.3 - 0.4) 7 Cấp 4 Nguy hiểm
NDVI (> 0.4) 9 Cấp 5 Rất nguy hiểm
2 Nhiệt độ (độ C)
≥ 26 9 Cấp 5 Rất nguy hiểm 24-26 7 Cấp 4 Nguy hiểm
22 độ - 24 5 Cấp 3 Nguy hiểm trung bình 19 độ - 22 3 Cấp 2 Ít nguy hiểm
17-19 độ 1 Cấp 1 Không nguy hiểm
3 Khí hậu
Rất ẩm 1 Cấp 1 Không nguy hiểm Ẩm 3 Cấp 2 Ít nguy hiểm
Hơi khô 5 Cấp 3 Nguy hiểm trung bình
Khô 7 Cấp 4 Nguy hiểm
Rất khô 9 Cấp 5 Rất nguy hiểm
4 Khoảng cách từ khu dân cư đến rừng (m) < 500 9 Cấp 5 Rất nguy hiểm 500 - 1.000 7 Cấp 4 Nguy hiểm
1.000 - 1.500 5 Cấp 3 Nguy hiểm trung bình 1.500 -2000 3 Cấp 2 Ít nguy hiểm
STT Lớp nhân tố
đầu vào Khoảng chỉ tiêu Điểm
Phân cấp
rủi ro Nguy cơ cháy
≥ 2000 1 Cấp 1 Không nguy hiểm
5 Khoảng cách từ đường giao thông đến rừng (m) < 500 9 Cấp 5 Rất nguy hiểm 500 - 1.000 7 Cấp 4 Nguy hiểm
1.000 - 1.500 5 Cấp 3 Nguy hiểm trung bình 1.500 - 2000 3 Cấp 2 Ít nguy hiểm
≥ 2000 1 Cấp 1 Không nguy hiểm
6 Thủy văn (km)
< 0-0.5 1 Cấp1 Không nguy hiểm 0.5-1 3 Cấp 2 Ít nguy hiểm
1-1.5 5 Cấp 3 Nguy hiểm trung bình 1.5-2 7 Cấp 4 Nguy hiểm ≥ 2 9 Cấp5 Rất nguy hiểm 7 Độ cao tuyệt đối (m) < 400 3 Cấp 2 Ít nguy hiểm 400 – 600 9 Cấp 1 Rất nguy hiểm 600 – 800 5 Cấp 3 Nguy hiểm
800 – 1000 7 Cấp 4 Nguy hiểm trung bình ≥ 1000 1 Cấp5 Không nguy hiểm
8 Độ dốc (độ)
0 – 8 1 Cấp 1 Không nguy hiểm 8 – 15 3 Cấp 2 Ít nguy hiểm
15 – 25 5 Cấp 3 Nguy hiểm trung bình