AC TRUNG THONG KE KHI TNG PHÂN CẤP NGUY CƠ CHÁY RỪNG BẰNG HÀM TÁCH BIỆT Nguyễn Văn Thêm(*) Tel 08 8 975 537; 0918 204 950 – Fax 84 8 8961707 (*) Khoa Lâm nghiệp, ĐHNL Tp HCM SUMMARY This paper present[.]
PHÂN CẤP NGUY CƠ CHÁY RỪNG BẰNG HÀM TÁCH BIỆT Nguyễn Văn Thêm(*) Tel: 08.8.975.537; 0918.204.950 – Fax: 84.8.8961707 (*) Khoa Lâm nghiệp, ĐHNL Tp.HCM SUMMARY This paper presents results of a forest fire classification which is based on the discriminant funtions The results show that the discriminant funtions with five variables (atmosphere temperature, atmosphere moisture, rain, drew point temperature and wind) can be used for forest fire grading purposes ĐẶT VẤN ĐỀ Lửa rừng nhân tố môi trường có tác dụng nhiều mặt Trước hết, lửa rừng làm biến đổi môi trường sâu sắc nhanh Lửa nhà lâm nghiệp kiểm sốt có số tác dụng tốt điều chỉnh thành phần vật chất hữu sàn rừng, tiêu diệt mầm bệnh, vật ký sinh côn trùng gây hại cho gỗ, loại bỏ sinh khối khô (cành nhánh, thân mục), làm giảm tích tụ vật chất hữu sàn rừng, cải thiện điều kiện tái sinh rừng hoạt động hệ động vật rừng Nhưng bên cạnh ưu điểm đó, lửa tự nhiên cháy rừng vô ý thức người nhân tố mơi trường có ảnh hưởng xấu không đến hiệu kinh doanh rừng, mà cịn đến tính mạng tài sản người Ở nước ta năm xảy cháy rừng với mức độ khốc liệt khác Một thống kê chưa đầy đủ cháy rừng cho thấy, từ năm 1963 – 2000 nước ta xảy 45.921 vụ cháy rừng, thiêu hủy 632.049 [2,3] Để chủ động phòng chống cháy rừng hạn chế hậu xấu cháy rừng gây ra, nhà quản lý rừng phải thực dự báo nguy cháy rừng Hiện cấp dự báo nguy cháy rừng ngắn hạn dài hạn nước ta xây dựng sở áp dụng tiêu khí tượng tổng hợp (P) V G Nexterov số ngày khô hạn Phạm Ngọc Hưng [2] Ưu điểm phương pháp đơn giản, dễ tính tốn Nhưng phương pháp có nhiều hạn chế: (1) chưa tính hết nhân tố ảnh hưởng đến cháy rừng; (2) để tính cấp cháy, hàng ngày dự báo viên phải đo đạc thông tin trường vào lúc 13 giờ; (3) cần phải có dụng cụ máy móc quan trắc khí tượng; (4) số P phân chia rộng, P thay đổi vài đơn vị cấp cháy thay đổi; (5) dự báo viên phải theo dõi liên tục ngày có mưa hay khơng mưa Với mong muốn góp phần tìm kiếm phương pháp dự báo nguy cháy rừng có hiệu cao, báo giới thiệu kết nghiên cứu phương pháp dự báo nguy cháy rừng dựa hàm tách biệt (Discriminant Functions) ĐỐI TƯỢNG VÀ PHƯƠNG PHÁP NGHIÊN CỨU Đối tượng nghiên cứu điều kiện khí hậu diễn biến thời tiết hàng ngày năm (1996 – 2002) khu vực Đà Lạt - tỉnh Lâm Đồng Thời gian theo dõi diễn biến thời tiết hàng ngày từ tháng 12 năm trước đến tháng năm sau Đây thời kỳ xác định có nguy cháy rừng cao Để xây dựng hàm dự đoán nguy cháy rừng sử dụng tiêu khí tượng phản ánh thời tiết hàng ngày tháng có nguy cháy rừng cao (từ tháng 12 năm trước đến tháng năm sau) – nhiệt độ trung bình khơng khí (T, oC), độ ẩm tương đối trung bình khơng khí (Rh,%), lượng mưa trung bình (M, mm), tốc độ gió trung bình (G, m/s) nhiệt độ điểm sương lúc 13 (D, mb) Các tài liệu thời tiết hàng ngày thu thập qua năm, 05/12/1996 đến ngày 31/05/2002; tổng số 992 ngày Ngoài ra, để phân loại sơ cấp nguy cháy rừng hàng ngày sử dụng tài liệu phân cấp cháy rừng tính theo tiêu tổng hợp (P) Nexterov Tất thơng tin khí tượng cấp cháy rừng hàng ngày theo tiêu tổng hợp Nexterov Chi Cục Kiểm Lâm tỉnh Lâm Đồng cung cấp [1] Từ tiêu khí tượng cấp cháy rừng hàng ngày khu vực Đà Lạt, xây dựng hàm tách biệt cho cấp nguy cháy rừng theo dạng: F(k) = a*T(k) + b*Rh(k) + c*M(k) + d*D + e*G; T(k), Rh(k), M(k), D(k) G(k) tương ứng biến hàm tách biệt thứ k hay cấp cháy rừng thứ k (k = 1, 2…5); a, b, c, d e hệ số hàm tách biệt Các hệ số hàm tách biệt xác định theo thủ tục lập nhóm phần mềm thống kê SPSS 10.0 (Statistical Products for Social Services) KẾT QUẢ VÀ THẢO LUẬN Từ cấp nguy cháy rừng phân loại theo tiêu tổng hợp Nexterov, thủ tục phân tích tách biệt cấp cháy theo nhân tố (T, Rh, M, G D, mb) nhận kết phân loại bảng hình Sơ đồ phân bố cấp cháy rừng hình vẽ từ điểm số ngày tương ứng với hai hàm hợp quy theo phương Từ hình dễ dàng nhận thấy ranh giới cấp cháy thiếu rõ ràng Phân tích số liệu bảng thấy 47,4% hay 471/992 ngày rơi vào cấp cháy I; cấp cháy II III, IV V tương ứng 28,3 24,3% số ngày Như vậy, sử dụng hàm tách biệt để phân loại cấp cháy rừng theo tiêu tổng hợp P Nexterov trung bình có 49,2% số ngày phân loại xác vào cấp cháy, 51,8% số ngày bị phân loại nhầm (bảng 2) Theo dõi số liệu bảng thấy tượng phân loại nhầm xảy tất cấp cháy; nhiều thuộc cấp cháy II - IV Sở dĩ có tượng phân loại nhầm số nguyên nhân sau đây: (1) Chỉ tiêu tổng hợp P tính từ biến nhiệt độ khơng khí; (2) Chỉ tiêu P phân cấp rộng, (3) P thay đổi vài đơn vị từ cấp cháy chuyển sang cấp cháy khác, điều kiện gây cháy rừng khơng có khác biệt rõ rệt Phân tích cho thấy, dựa vào tiêu tổng hợp P Nexterov để phân chia cấp nguy cháy rừng nhiều ngày có nguy cháy rừng giống bị phân cấp nhầm Rõ ràng để nâng cao độ xác việc phân chia cấp nguy cháy rừng, hệ thống phân cấp cháy phải bao gồm nhiều yếu tố gây cháy rừng Bảng Đặc trưng khí tượng từ tháng 12 năm trước đến tháng năm sau tương ứng với cấp nguy cháy rừng theo tiêu tổng hợp Nexterov Cấp cháy Đặc trưng N Trung bình Sai tiêu chuẩn I V% Xmin Xmax Xmax - Xmin N Trung bình Sai tiêu chuẩn II V% Xmin Xmax Xmax - Xmin N Trung bình Sai tiêu chuẩn III V% Xmin Xmax Xmax - Xmin N Trung bình Sai tiêu chuẩn IV V% Xmin Xmax Xmax - Xmin N Trung bình Sai tiêu chuẩn V V% Xmin Xmax Xmax - Xmin N Trung bình Sai tiêu chuẩn Tổng cộng V% Xmin Xmax Xmax - Xmin T0C 471 22,1 2,3 10,4 13 28 15 166 23,1 2,2 9,5 16 28,4 12,4 115 23,1 1,8 7,8 17 27 10 74 23,7 1,6 6,8 20 27 167 24,7 2,1 8,5 20 30 10 992 23 2,4 10,4 13 30 17 Am (%) 471 73,7 13,9 18,9 25 100 75 166 60,6 13,5 22,3 26 91 65 115 61,1 14,5 23,7 30 91 61 74 54,5 13,6 25 24 82 58 167 46 10,3 22,4 26 81 55 992 64 16,9 26,4 24 100 76 Mưa (mm) 470 5,3 10,7 201,9 80 80 166 0,4 1,5 375 10 10 115 0,3 1,3 433,3 8 74 0,3 1,7 566,7 13 13 167 0,1 1000 11 11 992 2,7 7,8 288,9 80 80 D (mb) 471 8,1 3,9 48,1 24,3 24,3 166 11,5 43,5 26,5 24,5 115 11,3 5,2 46 2,3 23,5 21,2 74 13,7 4,8 35 4,4 24,2 19,8 167 17,3 4,9 28,3 4,4 31,5 27,1 992 11 5,6 50,9 31,5 31,5 Gió (m/s) 471 7,2 69,4 37 35 166 6,7 4,7 70,1 35 33 115 6,8 4,9 72,1 30 28 74 6,7 3,2 47,8 18 16 167 6,2 3,1 50 18 16 992 6,9 4,6 66,7 37 35 Bảng Kết phân loại cấp cháy rừng theo tiêu P Nexterova Cấp cháy I II Tần số III IV V I II % III IV V Cấp cháy Tổng I II III IV V 323 22 75 32 18 470 41 19 47 28 31 166 38 12 23 19 23 115 22 13 24 74 14 31 110 167 68,7 4,7 16 6,8 3,8 100 24,7 11,4 28,3 16,9 18,7 100 33 10,4 20 16,5 20 100 9,5 10,8 29,7 17,6 32,4 100 1,8 5,4 8,4 18,6 65,9 100 a 49,2% số ngày ban đầu phân loại xác vào cấp cháy Hàm Để nhận hàm dự đốn nguy cháy rừng có hiệu hơn, tác giả sử dụng biến dự đốn nhiệt độ khơng khí (T, 0C), độ ẩm khơng khí (Rh, %), lượng mưa ngày (M, mm), nhiệt độ điểm sương (D, mb) tốc độ gió (G, m/s) Canonical Discriminant Bằng thủ tục phân tích tách biệt 10 cấp nguy cháy rừng theo biến phân loại nhận kết phân loại cháy rừng bảng – Nexterov hình Số liệu bảng Group Centroids đặc trưng khí tượng từ tháng 12 5 năm trước đến tháng năm sau 4 23 tương ứng với cấp nguy Fu nct cháy rừng phân chia theo ion -2 2 hàm tách biệt với biến số Số -4 liệu bảng kết kiểm -4 -2 định ngang trung bình Hàm Function biến số dự đoán cấp Hình Đồ thị biểu diễn cấp nguy cháy rừng nguy cháy rừng Từ cho theo phân cấp Nexterov thấy biến phân loại có khác biệt lớn cấp cháy Bảng hệ số hàm phân loại hợp quy chuẩn hố Qua nhận thấy biến có ý nghĩa lớn hàm phân loại nhiệt độ, độ ẩm nhiệt độ điểm sương Số liệu bảng báo cáo hệ số hàm phân loại hợp quy chưa chuẩn hoá Các hệ số sử dụng để tính điểm số ngày thuộc hàm hợp quy khác Để biết điểm số ngày tương ứng với hàm hợp quy, dự báo viên thay biến dự đoán vào hàm tương ứng Điểm số trung bình cấp cháy ghi lại bảng Từ số liệu bảng cho thấy trọng tâm hàm hợp quy khác lớn nhóm Số liệu bảng hệ số hàm phân loại tuyến tính Fisher Chúng sử dụng để tính khoảng cách khác cực đại cấp cháy rừng Lưu ý rằng, sử dụng số liệu bảng để tính khoảng cách khác cực đại cấp cháy 100% số ngày phân loại xác vào cấp nguy cháy rừng Số ngày phân bố vào cấp cháy dẫn bảng 10 Bảng Đặc trưng khí tượng từ tháng 12 năm trước đến tháng năm sau tương ứng với cấp cháy rừng phân loại theo hàm tách biệt biến số Cấp cháy I II III IV V Tổng Đặc trưng N Trung bình Sai tiêu chuẩn V% Xmin Xmax Xmax - Xmin N Trung bình Sai tiêu chuẩn V% Xmin Xmax Xmax - Xmin N Trung bình Sai tiêu chuẩn V% Xmin Xmax Xmax - Xmin N Trung bình Sai tiêu chuẩn V% Xmin Xmax Xmax - Xmin N Trung bình Sai tiêu chuẩn V% Xmin Xmax Xmax - Xmin N Trung bình Sai tiêu chuẩn V% Xmin Xmax Xmax - Xmin T0C 191,0 20,6 1,9 9,2 13,0 25,0 12,0 222,0 23,3 1,9 8,2 16,0 28,0 12,0 277,0 22,7 1,8 7,9 15,0 26,0 11,0 184,0 23,7 1,9 8,0 18,0 29,0 11,0 118,0 25,7 2,0 7,8 20,0 30,0 10,0 992,0 23,0 2,4 10,4 13,0 30,0 17,0 Am (%) 191,0 87,5 7,5 8,6 69,0 100,0 31,0 222,0 61,6 3,9 6,3 53,0 73,0 20,0 277,0 71,7 4,2 5,9 59,0 82,0 23,0 184,0 48,9 4,0 8,2 39,0 58,0 19,0 118,0 35,6 4,5 12,6 24,0 42,0 18,0 992,0 63,9 16,9 26,4 24,0 100,0 76,0 Mưa (mm) 191,0 8,7 13,5 155,2 0,0 80,0 80,0 222,0 1,7 5,8 341,2 0,0 58,0 58,0 277,0 1,9 5,1 268,4 0,0 38,0 38,0 184,0 0,6 3,5 583,3 0,0 37,0 37,0 118,0 0,0 0,0 0,0 0,0 0,0 992,0 2,7 7,8 288,9 0,0 80,0 80,0 D (mb) 191,0 5,0 3,1 62,0 0,0 11,5 11,5 222,0 11,2 2,1 18,8 4,2 17,2 13,0 277,0 7,8 1,7 21,8 2,8 12,8 10,0 184,0 15,0 2,4 16,0 8,9 21,8 12,9 118,0 21,5 3,1 14,4 14,0 31,5 17,5 992,0 11,0 5,6 50,9 0,0 31,5 31,5 Bảng Kiểm định ngang trung bình nhóm Biến dự đốn Wilks' Lambda F df1 T (độ C) 0,626 147,7 Rh (%) 0,084 2684,4 Mưa (mm) 0,854 42,3 D (mb) 0,185 1087,5 Gió (m/s) 0,969 7,8 df2 987 987 987 987 987 Sig 0,000 0,000 0,000 0,000 0,000 Gió (m/s) 191,0 8,2 6,5 79,3 2,0 37,0 35,0 222,0 6,2 3,5 56,5 2,0 22,0 20,0 277,0 7,2 4,7 65,3 2,0 35,0 33,0 184,0 6,4 3,7 57,8 2,0 25,0 23,0 118,0 5,8 2,6 44,8 2,0 18,0 16,0 992,0 6,9 4,6 66,7 2,0 37,0 35,0 Bảng Các hệ số hàm phân loại hợp quy chuẩn hoá Biến phân loại T (0C) Rh (%) Mưa (mm) D (mb) Gió (m/s) Hàm phân loại hợp quy chuẩn hoá 0,730 -0,762 0,850 -0,115 1,058 0,451 0,450 0,021 -0,152 0,236 -0,333 -0,352 -0,628 1,025 0,378 0,085 0,520 0,080 0,266 0,872 Bảng Các hệ số hàm phân loại hợp quy chưa chuẩn hoá Biến phân loại T (độ C) Rh (%) Mưa (mm) D (mb) Gió (m/s) (Constant) Hàm phân loại hợp quy chưa chuẩn hoá 0,387 -0,404 0,451 -0,061 0,215 0,092 0,091 0,004 -0,021 0,033 -0,046 -0,049 -0,259 0,424 0,156 0,035 0,115 0,018 0,059 0,193 -20,529 -1,435 -18,198 -0,449 Bảng Các hệ số hàm trung tâm nhóm Cấp cháy I II III IV V 5,72 -0,48 2,43 -4,01 -7,81 Hàm hợp quy trung tâm nhóm 0,81 -0,05 0,00 -0,33 -0,03 -0,11 -0,53 0,15 0,05 -0,04 -0,35 0,05 0,64 0,33 0,01 Bảng Kiểm định hàm hợp quy Kiểm định hàm Wilks' Lambda df P χ2 1- 0,039 3207,550 20 0,000 2- 0,747 287,789 12 0,000 3-4 0,955 45,613 0,000 0,996 3,662 0,160 Từ số liệu bảng 10 cho thấy số ngày rơi vào cấp cháy I 191 ngày hay 19,3%; cấp cháy II III, IV V tương ứng 50,3 30,4% So sánh số liệu bảng với số liệu bảng 10 thấy đặc trưng khí tượng cấp cháy có thay đổi lớn Điều xảy ngày có điều kiện khí tượng tương đồng với phân loại vào cấp cháy Canonical Discriminant Functions Them12 Group Centroids Ungrouped Cases Hàm Function 2 -2 -4 -20 -10 Function 10 Hàm Hình Đồ thị phân chia cấp nguy cháy rừng theo hàm phân loại tuyến tính Fisher Bảng Các hệ số hàm phân loại cấp nguy cháy rừng Biến dự đoán T (độ C) Rh (%) Mưa (mm) D (mb) Gió (m/s) (Constant) Cấp nguy cháy rừng I II III IV V 18,061 16,134 17,416 14,497 13,061 7,391 5,955 6,580 5,194 4,500 -1,051 -0,954 -1,037 -0,863 -0,791 -3,415 -2,288 -3,096 -1,297 8,366E-02 4,459 3,707 4,080 3,317 2,924 -515,994 -370,884 -436,683 -301,021 -258,997 a Hàm phân loại tuyến tính Fisher Kết báo cáo bảng 10 100% số ngày phân loại xác vào cấp nguy cháy rừng Ranh giới cấp xác định rõ ràng (hình 2) Do đó, số liệu bảng sử dụng để phân chia cấp nguy cháy rừng cho khu vực Đà Lạt tỉnh Lâm Đồng Từ số liệu bảng xây dựng hàm phân loại tuyến tính Fisher để phân chia cấp nguy cháy rừng sau: Cấp I = 18,061T + 7,391Rh – 1,051M – 3,415D + 4,459G – 515,994 (1) Cấp II = 16,134T + 5,955Rh – 0,954M – 2,288D + 3,707G –370,884 (2) Cấp III = 17,416T + 6,580Rh – 1,037M – 3,096D + 4,080G – 436,683 (3) Cấp IV = 14,497T + 5,194Rh – 0,863M – 1,297D + 3,317G – 258,997 (4) Cấp V = 13,061 + 4,500Rh – 0,791M +8,366D + 2,924G – 258,997 (5) Lưu ý rằng, hàm – sử dụng để dự đoán ngày thuộc vào cấp cháy Về ý nghĩa thực tiễn, ứng dụng hàm - để phân cấp nguy cháy rừng ngày đó, trước hết dự báo viên cần đo đạc biến dự đoán (T, Rh, M, D G) ngày Tiếp theo biến số dự đốn vào hàm phân loại tính giá trị khoảng cách cực đại cho hàm Hàm nhận giá trị khoảng cách lớn cho biết ngày thuộc cấp cháy Cách tính tốn cấp nguy cháy rừng cho ngày cụ thể dẫn bảng 11 Bảng 10 Kết phân loại số ngày vào cấp cháy rừnga Cấp cháy Tần số % 5 Số ngày dự đoán vào Tổng cấp nguy cháy rừng 191 191 222 222 277 277 184 184 118 118 100 100 100 100 100 100 100 100 100 100 a 100% số ngày phân loại xác vào cấp cháy Bảng 11 Kết dự đoán cấp nguy cháy rừng TT Đặc trưng khí tượng T Kết hàm phân loại Mưa D Gió 72 8,2 442,0 442,6 447,8 466,8 391,8 IV 24/3/2002 20,0 100 20 576,7 539,3 561,1 543,0 445,2 I 23/3/2002 23,4 Cấp cháy Rh Cần nhận thấy rằng, giống tiêu tổng hợp Nexterov, hàm phân loại nguy cháy rừng tính cho thời điểm (ngày tháng) năm khác nhận giá trị khác Thật vậy, ngày 01 tháng 02 năm 1997, 1998 2002 thuộc cấp nguy cháy V (cực kỳ nguy hiểm), vào ngày năm 1999, 2000 2001 lại thuộc cấp cháy II – III (nguy cháy thấp đến trung bình) Tình xảy khí hậu yếu tố thời tiết có biến động từ năm qua năm khác Chính thế, bảng phân cấp nguy cháy rừng cho ngày tháng thời kỳ có nguy cháy rừng cần xây dựng từ giá trị cấp cháy trung bình nhiều năm Theo đó, từ số liệu cấp cháy năm (1996 – 2002) tác giả xây dựng bảng dự báo nguy cháy rừng hàng ngày cho khu vực Đà Lạt tỉnh Lâm Đồng Kết ghi lại bảng 12 Bảng 11 Cấp nguy cháy rừng khu vực Đà Lạt – tỉnh Lâm Đồng Ngày tháng Tháng có nguy cháy rừng 12 Nex(*) Dis(**) Nex Dis Nex Dis Nex Dis Nex Dis Nex Dis (1) (2) (3) (4) (5) (6) (7) (8) (9) (10) (11) (12) (13) 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 2 2 2 2 1 1 1 1 1 1 1 2 2 2 2 2 2 2 2 1 1 1 1 1 1 1 2 2 2 2 2 2 3 2 2 3 3 3 3 3 3 3 4 4 4 2 2 3 2 2 3 3 3 3 3 3 3 4 4 4 4 4 3 3 3 3 4 4 4 3 3 3 3 4 4 3 3 3 3 4 4 4 3 3 3 3 4 4 4 4 4 4 3 3 3 3 3 3 2 2 2 4 4 4 4 4 3 3 3 3 3 3 2 2 2 2 2 2 2 2 2 1 1 1 1 1 1 1 1 2 2 2 2 2 2 1 1 1 1 1 1 1 1 1 1 2 2 2 2 1 1 1 1 1 1 1 1 1 1 1 2 2 2 2 1 1 1 1 1 1 1 1 1 1 (*) Nex = phân cấp nguy cháy rừng theo tiêu tổng hợp Nexterov; (**) Dis = phân cấp nguy cháy rừng theo hàm phân loại tuyến tính Fisher Trong bảng 12, trị số cột chẵn cấp nguy cháy rừng trung bình tính theo tiêu tổng hợp Nexterov, số liệu cột lẻ tính theo hàm phân loại tuyến tính Fisher Quan sát số liệu bảng 12 dễ dàng nhận thấy cấp cháy trung bình vào thời điểm tính theo tiêu tổng hợp Nexterov theo hàm phân loại tuyến tính Fisher khơng có khác biệt Từ số liệu bảng 12, dự báo viên dự đoán nguy cháy rừng xảy tuần tháng có nguy cháy rừng So với phương pháp Nexterov, phương pháp dự đoán cấp nguy cháy rừng hàng ngày theo hàm tách biệt tính tốn dễ dàng chủ động Một ưu điểm lớn hàm tách biệt dự báo viên tính cấp cháy cho ngày bất kỳ, yếu tố khí tượng khơng cập nhật liên tục hàng ngày Vì lý đó, tác giả khuyến nghị đơn vị lâm nghiệp khu vực Đà Lạt tỉnh Lâm Đồng nơi khác nên phân loại cấp nguy cháy rừng theo hàm phân loại tuyến tính Fisher KẾT LUẬN Nghiên cứu chứng tỏ rằng, hàm phân loại tuyến tính Fisher xây dựng nhiều biến số sử dụng vào mục đích dự báo nguy cháy rừng Để xây dựng hàm phân loại nguy cháy rừng dùng vào mục đích dự báo cấp cháy rừng cần phải có thơng tin khí tượng nhiều năm Tác giả khuyến nghị đơn vị lâm nghiệp khu vực Đà Lạt tỉnh Lâm Đồng sử dụng kết nghiên cứu phối hợp với tiêu tổng hợp Nexterov để dự báo cấp nguy cháy rừng TÀI LIỆU THAM KHẢO CHÍNH Chi Cục Kiểm Lâm Nhân Dân tỉnh Lâm Đồng, 2002 Số liệu khí tượng năm 1996 – 2002 Tài liệu lưu hành nội Phạm Ngọc Hưng, 2001 Thiên tai khô hạn cháy rừng giải pháp phòng cháy chữa cháy rừng Việt Nam, Nxb Nơng Nghiệp, Hà Nội Đặng Vũ Cẩn, Hồng Kim Ngũ ctv, 1992 Quản lý bảo vệ rừng Tủ sách Trường Đại Học Lâm Nghiệp, Xuân Mai - Hà Tây Nguyễn Văn Thêm, 2002 Sinh thái rừng, Nxb Nông Nghiệp, Chi nhánh TP Hồ Chí Minh Nguyễn Văn Thêm, 2002 Hướng dẫn sử dụng phần mềm SPSS 10.0 để xử lý thông tin nông – lâm nghiệp Sách lưu hành nội 10 ... phân chia cấp nguy cháy rừng theo hàm phân loại tuyến tính Fisher Bảng Các hệ số hàm phân loại cấp nguy cháy rừng Biến dự đoán T (độ C) Rh (%) Mưa (mm) D (mb) Gió (m/s) (Constant) Cấp nguy cháy rừng... tiêu tổng hợp P Nexterov để phân chia cấp nguy cháy rừng nhiều ngày có nguy cháy rừng giống bị phân cấp nhầm Rõ ràng để nâng cao độ xác việc phân chia cấp nguy cháy rừng, hệ thống phân cấp cháy... Nexterov, hàm phân loại nguy cháy rừng tính cho thời điểm (ngày tháng) năm khác nhận giá trị khác Thật vậy, ngày 01 tháng 02 năm 1997, 1998 2002 thuộc cấp nguy cháy V (cực kỳ nguy hiểm), vào ngày