Chương 3. KẾT QUẢ NGHIÊN CỨU VÀ THẢO LUẬN
3.5. Dự báo một số đặc trƣng của VRR
3.5.3. Dự báo khối lượng VRR phân hủy
3.5.3.1. Xác định các tham số của phương trình dự báo khối lượng VRR phân hủy Tham số của các phương trình được lựa chọn để dự báo lượng VRR phân hủy đƣợc tổng hợp tại bảng 3.38.
Bảng 3.38. Xác định tham số của phương trình dự báo lượng VRR phân hủy
Trạng thái TTV
Lƣợng tồn dƣ (X) tại thời điểm
Các thời điểm dự báo
phân hủy
R^2
Tham số
Sig
a b
Rừng giàu
T1
T1 0,64 251,64 0,23 2,1E-14
T2 0,96 40,21 0,41 9,4E-43
T3 0,99 4,26 0,57 5,1E-66
T4 0,96 -115,57 0,55 1,3E-40
T2
T1 0,60 355,13 0,27 4,9E-13
T2 0,95 107,19 0,48 3,4E-38
T3 0,99 66,28 0,67 5,9E-56
T4 0,97 -128,52 0,66 5,1E-46
T3
T1 0,66 262,44 0,39 6,9E-15
T2 0,98 74,41 0,67 9,2E-49
T3 0,997 92,14 0,94 1,6E-72
T4 0,93 56,90 0,89 4,6E-35
T4
T1 0,82 14,29 0,30 4,1E-23
T2 0,99 161,26 0,46 4,9E-57
T3 0,95 390,60 0,63 1,4E-39
T4 0,78 738,46 0,56 3,5E-20
Rừng Trung bình
T1
T1 0,71 115,76 0,30 2,8E-17
T2 0,98 45,26 0,46 3,0E-52
T3 0,995 112,87 0,57 1,4E-68
T4 0,96 127,19 0,48 5,7E-43
T2
T1 0,74 9,52 0,37 1,2E-18
T2 0,99 -43,34 0,57 1,6E-56
T3 0,999 4,66 0,71 7,5E-92
T4 0,96 45,09 0,59 4,6E-43
T3
T1 0,59 265,02 0,45 7,3E-13
T2 0,88 125,13 0,72 2,6E-28
T3 0,94 64,51 0,93 8,8E-38
T4 0,97 -42,44 0,80 4,8E-44
T4
T1 0,69 40,11 0,35 1,8E-16
T2 0,87 115,50 0,51 8,1E-28
T3 0,91 136,91 0,65 9,2E-32
T4 0,89 113,32 0,55 1,1E-29
Rừng nghèo T1 T1 0,69 134,34 0,34 1,4E-16
Trạng thái TTV
Lƣợng tồn dƣ (X) tại thời điểm
Các thời điểm dự báo
phân hủy
R^2
Tham số
Sig
a b
T2 0,89 20,69 0,52 9,9E-30
T3 0,98 -24,10 0,73 1,7E-51
T4 0,87 179,63 0,58 6,0E-27
T2
T1 0,66 152,60 0,44 2,7E-15
T2 0,83 76,29 0,66 2,6E-24
T3 0,94 9,48 0,94 1,1E-36
T4 0,94 -3,43 0,80 8,2E-37
T3
T1 0,48 710,56 0,55 1,7E-07
T2 0,54 788,86 0,87 2,7E-11
T3 0,70 796,49 1,34 7,2E-17
T4 0,99 145,82 1,35 2,3E-57
T4
T1 0,84 -12,31 0,56 3,9E-25
T2 0,99 -55,61 0,81 1,7E-55
T3 0,99 49,12 1,10 1,8E-63
T4 0,72 545,07 0,79 1,7E-17
Từ bảng 3.38 nhận thấy, hệ số xác định R2 của phương trình Y=a+ b*X dao động từ 0,48 đến 0,99. Nhƣ vậy, có thể khẳng định, quan hệ giữa lƣợng VRR bổ sung với VRR tồn dư tồn tại ở mức tương quan vừa phải đến rất chặt theo dạng phương trình Y=a+ b*X. Giá trị Sig của tất cả các phương trình đều nhỏ hơn 0,05 cho thấy phương trình các tham số của phương trình đều tồn tại và đó thể hiện sự tồn tại của phương trình.
Hình 3.26. Quan hệ giữa khối lƣợng VRR phân hủy với khối lƣợng VRR tồn dƣ 3.5.3.2. Xác định sai số của phương trình dự báo khối lượng VRR bổ sung
Dựa vào các tham số của các phương trình dự báo khối lượng VRR phân hủy và theo các công thức tính sai số. Xác định các loại sai số của phương trình dự báo lƣợng VRR phân hủy đƣợc tổng hợp tại bảng 3.39.
Từ bảng 3.39 có thể rút ra một số nhận xét về sai số của phương trình dự báo khối lƣợng VRR bổ sung tại các thời điểm khác nhau trong năm nhƣ sau:
Tỷ lệ ô kiểm tra mắc sai số âm trung bình là 43,2% và sai số dương là 56,8%.
Sai số lớn nhất mắc phải từ 1,8% đến 72,4%, trung bình 29,2%.
Sai số bình quân về khối lƣợng dao động từ 1,8% đến 30,6%, bình quân là 11,3%,
Sai số quân phương về khối lượng dao động từ 2,2% đến 39,6%, bình quân là 13,3%.
Hệ số chính xác từ 0,6% đến 10,2%, bình quân là 3,4%.
Sai số tổng thể về khối lƣợng từ - 7,9% đến 4,3%, bình quân là - 0,4%.
Bảng 3.39. Sai số xác định lƣợng VRR phân hủy
Trạng thái TTV
Tồn dƣ (X) tại
thời điểm
Các thời điểm dự báo
phân hủy
Các loại sai số
sai số (-) sai số (+) ∆Max ∆TB Sq P% ∆%(∑PH)
Rừng giàu
T1
T1 5 10 67,8 32,6 39,3 10,2 2,5
T2 9 6 15,1 9,2 10,5 2,7 1,6
T3 9 6 5,2 3,5 3,8 1,0 0,7
T4 6 9 15,9 9,5 11,3 2,9 -1,3
T2
T1 5 10 77,9 34,4 42,3 10,9 2,4
T2 9 6 19,8 10,8 12,7 3,3 1,7
T3 9 6 8,8 5,1 5,9 1,5 0,9
T4 6 9 12,6 7,6 8,9 2,3 -1,1
T3
T1 5 10 62,0 31,5 36,3 9,4 1,6
T2 5 10 10,9 6,6 7,4 1,9 0,6
T3 9 6 4,3 2,0 2,5 0,6 -0,4
T4 6 9 15,7 12,3 13,5 3,5 -2,6
T4
T1 5 10 23,8 20,6 21,5 5,6 0,4
T2 6 9 9,5 4,8 5,8 1,5 -1,5
T3 6 9 17,8 10,0 11,8 3,0 -2,6
T4 6 9 49,7 22,7 27,3 7,0 -5,2
Rừng Trung bình
T1
T1 5 10 63,6 32,9 40,2 10,4 4,3
T2 5 10 9,8 6,4 7,3 1,9 0,8
T3 8 7 4,5 2,6 3,0 0,8 -0,8
T4 6 9 13,7 9,0 9,9 2,6 -2,5
T2
T1 5 10 55,7 31,4 39,3 10,2 5,0
T2 9 6 9,5 5,8 6,6 1,7 1,4
T3 6 9 1,8 1,3 1,5 0,4 -0,3
T4 6 9 13,0 9,0 10,4 2,7 -2,0
T3
T1 9 6 40,5 43,1 56,2 14,5 7,0
T2 9 6 31,1 16,4 19,3 5,0 3,9
T3 4 11 15,7 9,4 11,3 2,9 2,5
Trạng thái TTV
Tồn dƣ (X) tại
thời điểm
Các thời điểm dự báo
phân hủy
Các loại sai số
sai số (-) sai số (+) ∆Max ∆TB Sq P% ∆%(∑PH)
T4 4 11 11,2 7,1 8,4 2,2 1,0
T4
T1 9 6 71,5 38,7 48,8 12,6 7,4
T2 4 11 24,9 15,6 17,6 4,5 3,5
T3 4 11 17,4 12,6 13,6 3,5 1,9
T4 8 7 27,8 13,1 16,9 4,4 0,2
Rừng nghèo
T1
T1 5 10 48,3 28,3 31,0 8,0 2,3
T2 5 10 27,8 16,4 20,2 5,2 3,0
T3 9 6 11,0 6,9 8,0 2,1 1,5
T4 6 9 22,5 15,6 17,0 4,4 -4,7
T2
T1 5 10 60,2 30,7 36,0 9,3 3,7
T2 9 6 41,5 21,9 26,6 6,9 4,3
T3 9 6 20,0 12,2 14,1 3,6 3,0
T4 6 9 14,1 10,6 11,4 2,9 -2,8
T3
T1 6 9 64,8 41,9 58,5 12,1 3,5
T2 7 8 58,7 36,3 48,0 11,4 4,7
T3 7 8 72,4 27,2 34,2 8,8 4,2
T4 4 11 9,9 4,3 5,5 1,4 0,2
T4
T1 5 10 21,9 18,0 18,9 4,9 0,8
T2 5 10 9,4 5,9 6,9 1,8 1,0
T3 10 5 7,3 4,1 4,7 1,2 -0,7
T4 6 9 48,2 24,2 27,9 7,2 -7,1
Trung bình 6 9 29,2 16,3 19,6 5,2 1,0
Min 4 5 1,8 1,3 1,5 0,4 -7,1
Max 10 11 72,4 43,1 58,5 14,5 7,4
Nhận xét
Các phương trình trên sẽ góp phần dự báo được một số đặc trưng của VRR nhƣ: khối lƣợng VRR bổ sung, tồn dƣ, phân hủy tại các thời điểm khác nhau trong năm. Tuy nhiên, việc dự báo trên chỉ căn cứ vào khối lƣợng VRR tồn dƣ hiện có
trên nền rừng để dự báo cho các trạng thái rừng giàu, rừng trung bình và rừng nghèo, nên kết quả dự báo sẽ có sai số sử dụng cao do có sự biến động lớn của khối lƣợng VRR tồn dƣ hiện có trên nền rừng. Mặt khác, các khái niệm về "rừng giàu",
"rừng trung bình", "rừng nghèo" có sự biến động rất lớn, nên kết quả dự báo cho từng trạng thái rừng sẽ có sai số lớn. Mặc dù vậy, đây là phương pháp dự báo nhanh, đơn giản, dễ làm, dễ thực hiện. Phương pháp này được sử dụng khi cần dự báo nhanh và không yêu cầu cao về độ chính xác.
Để hạn chế những nhƣợc điểm của nội dung dự báo trên, luận án tiếp tục nghiên cứu và đưa ra phương pháp dự báo một số đặc trưng của VRR với sự tham gia của các nhân tố có ảnh hưởng tới VRR, gồm: i) khối lượng VRR tồn dư tại thời điểm dự báo Wtdư; ii) trữ lượng của rừng tại thời điểm dự báo M.
Về cả lý luận và thực tiễn, đây là những nhân tố có ảnh hưởng quan trọng đến VRR. Kết quả đƣợc tổng hợp tại bảng 3.40.
Bảng 3.40. Phương trình dự báo khối lượng VRR bổ sung và phân hủy Thời điểm Phương trình dự báo khối lượng VRR bổ sung R2 Sig
T1 YWBS= - 0,8184 + 0,4638*Wtdƣ + 0,0063*M 0,85 4E-05 T2 YWBS= 0,9239 + 0,0830*Wtdƣ + 0,0080*M 0,80 2E-04 T3 YWBS= 0,6384 + 0,3426*Wtdƣ + 0,0082*M 0,86 2E-05 T4 YWBS= 1,2591 + 0,4798*Wtdƣ + 0,0163*M 0,90 2E-06 Thời điểm Phương trình dự báo khối lượng VRR phân hủy R2 Sig
T1 YWPH= 1,0413 + 0,0088*Wtdƣ + 0,0143*M 0,86 2E-02 T2 YWPH= 0,9475 + 0,2509*Wtdƣ + 0,0091M 0,87 9E-06 T3 YWPH= 0,9239 + 0,6047*Wtdƣ + 0,0079*M 0,89 5E-06 T4 YWPH= 1,4727 + 0,2413*Wtdƣ + 0,0140*M 0,80 2E-04 Bảng 3.40 cho thấy, quan hệ F (Wtdư, M) có hệ số tương quan R từ 0,66 đến 0,90, nên quan hệ giữa khối lƣợng VRR với các nhân tố (Wtdƣ, M) ở mức chặt. Giá trị Sig đều nhỏ hơn 0,05, nên các tham số của phương trình đều tồn tại và qua đó khẳng định sự tồn tại của phương trình.
Lƣợng VRR bổ sung và phân hủy đƣợc dự báo cho từng thời điểm khác nhau trong năm. Khi biết đƣợc khối lƣợng VRR tại từng thời điểm sẽ xác định đƣợc tổng
khối lượng VRR của một năm tùy theo mục đích của người sử dụng. Từ kết quả xác định tham số của phương trình dự đoán khối lượng VRR bổ sung và phân hủy đã xây dựng bảng tra khối lƣợng VRR bổ sung và phân hủy ở các thời điểm trong năm đƣợc thể hiện ở bảng 3.41.
Bảng 3.41. Bảng tra khối lƣợng VRR Dữ liệu đầu vào Dữ liệu đầu ra (tấn/ha) Wtdƣ (tấn/ha) M
(m3/ha)
Tra khối lƣợng VRR bổ sung Tra khối lƣợng VRR phân hủy
T1 T2 T3 T4 T1 T2 T3 T4 T1 T2 T3 T4
11 9,5 7 9 30 4,47 1,95 3,28 6,07 1,57 3,60 5,39 4,07
... ... ... ... ... ... ... ... ... ... ... ... ...
11 9,5 7 9 50 4,60 2,11 3,45 6,39 1,85 3,79 5,55 4,35 11 9,5 7 9 55 4,63 2,15 3,49 6,47 1,93 3,83 5,59 4,42
... ... ... ... ... ... ... ... ... ... ... ... ...
11 9,5 7 9 100 4,91 2,51 3,86 7,21 2,57 4,24 5,95 5,05 11 9,5 7 9 105 4,95 2,55 3,90 7,29 2,64 4,29 5,99 5,12 11 9,5 7 9 110 4,98 2,60 3,94 7,37 2,71 4,33 6,03 5,19
... ... ... ... ... ... ... ... ... ... ... ... ...
11 9,5 7 9 250 5,86 3,72 5,10 9,66 4,72 5,61 7,14 7,15 11 9,5 7 9 255 5,89 3,76 5,14 9,74 4,79 5,65 7,18 7,22 11 9,5 7 9 260 5,92 3,80 5,18 9,82 4,86 5,70 7,22 7,29
... ... ... ... ... ... ... ... ... ... ... ... ...
11 9,5 7 9 300 6,18 4,12 5,51 10,47 5,44 6,06 7,53 7,85 11 9,5 7 9 305 6,21 4,16 5,55 10,55 5,51 6,11 7,57 7,92 11 9,5 7 9 310 6,24 4,20 5,59 10,63 5,58 6,15 7,61 7,99
... ... ... ... ... ... ... ... ... ... ... ... ...
12 10 9 11 30 4,85 2,03 3,83 7,03 1,57 3,85 6,36 4,55 12 10 9 11 35 4,88 2,07 3,87 7,12 1,65 3,90 6,40 4,62 12 10 9 11 40 4,91 2,11 3,91 7,20 1,72 3,94 6,44 4,69
... ... ... ... ... ... ... ... ... ... ... ... ...
12 10 9 11 75 5,13 2,40 4,20 7,77 2,22 4,26 6,72 5,18 12 10 9 11 80 5,17 2,44 4,24 7,85 2,29 4,31 6,76 5,25 12 10 9 11 85 5,20 2,48 4,28 7,93 2,36 4,35 6,80 5,32
... ... ... ... ... ... ... ... ... ... ... ... ...
Dữ liệu đầu vào Dữ liệu đầu ra (tấn/ha) Wtdƣ (tấn/ha) M
(m3/ha)
Tra khối lƣợng VRR bổ sung Tra khối lƣợng VRR phân hủy
T1 T2 T3 T4 T1 T2 T3 T4 T1 T2 T3 T4
12 10 9 11 130 5,48 2,84 4,66 8,66 3,01 4,76 7,15 5,95 12 10 9 11 135 5,51 2,88 4,70 8,75 3,08 4,81 7,19 6,03 12 10 9 11 140 5,54 2,92 4,74 8,83 3,15 4,85 7,23 6,10
... ... ... ... ... ... ... ... ... ... ... ... ...
12 10 9 11 190 5,86 3,32 5,15 9,64 3,87 5,31 7,63 6,80 12 10 9 11 195 5,89 3,36 5,19 9,73 3,94 5,35 7,67 6,87 12 10 9 11 200 5,92 3,40 5,23 9,81 4,01 5,40 7,71 6,94
... ... ... ... ... ... ... ... ... ... ... ... ...
12 10 9 11 260 6,30 3,88 5,73 10,79 4,87 5,94 8,18 7,78 12 10 9 11 265 6,33 3,92 5,77 10,87 4,94 5,99 8,22 7,85 12 10 9 11 270 6,36 3,96 5,81 10,95 5,01 6,03 8,26 7,92
... ... ... ... ... ... ... ... ... ... ... ... ...
12 10 9 11 330 6,74 4,44 6,30 11,93 5,87 6,58 8,74 8,76 12 10 9 11 335 6,77 4,48 6,35 12,01 5,95 6,63 8,78 8,83 12 10 9 11 340 6,81 4,52 6,39 12,09 6,02 6,67 8,82 8,90
... ... ... ... ... ... ... ... ... ... ... ... ...
12 10 9 11 385 7,09 4,88 6,76 12,82 6,66 7,08 9,17 9,53 12 10 9 11 390 7,12 4,92 6,80 12,91 6,73 7,13 9,21 9,60 12 10 9 11 395 7,15 4,96 6,84 12,99 6,81 7,17 9,25 9,67
... ... ... ... ... ... ... ... ... ... ... ... ...
14 11 9 12 55 6,08 2,30 4,27 8,03 1,95 4,27 6,97 5,20 14 11 9 12 60 6,12 2,34 4,31 8,11 2,03 4,31 7,01 5,27 14 11 9 12 65 6,15 2,38 4,35 8,19 2,10 4,36 7,05 5,34
... ... ... ... ... ... ... ... ... ... ... ... ...
14 11 9 12 135 6,59 2,94 4,93 9,33 3,10 4,99 7,60 6,32 14 11 9 12 140 6,62 2,98 4,97 9,41 3,17 5,04 7,64 6,39 14 11 9 12 145 6,65 3,02 5,01 9,50 3,24 5,09 7,68 6,46
... ... ... ... ... ... ... ... ... ... ... ... ...
14 11 9 12 210 7,06 3,54 5,55 10,56 4,18 5,68 8,20 7,37 14 11 9 12 215 7,09 3,58 5,59 10,64 4,25 5,72 8,24 7,44 14 11 9 12 220 7,12 3,62 5,63 10,72 4,32 5,77 8,27 7,51
Dữ liệu đầu vào Dữ liệu đầu ra (tấn/ha) Wtdƣ (tấn/ha) M
(m3/ha)
Tra khối lƣợng VRR bổ sung Tra khối lƣợng VRR phân hủy
T1 T2 T3 T4 T1 T2 T3 T4 T1 T2 T3 T4
... ... ... ... ... ... ... ... ... ... ... ... ...
14 11 9 12 310 7,69 4,34 6,37 12,19 5,61 6,59 8,99 8,78 14 11 9 12 315 7,72 4,38 6,41 12,27 5,68 6,63 9,03 8,85 14 11 9 12 320 7,76 4,42 6,45 12,35 5,75 6,68 9,07 8,92
... ... ... ... ... ... ... ... ... ... ... ... ...
14 11 9 12 385 8,17 4,95 6,99 13,41 6,68 7,27 9,58 9,83 14 11 9 12 390 8,20 4,99 7,03 13,49 6,76 7,31 9,62 9,90 14 11 9 12 395 8,23 5,03 7,07 13,57 6,83 7,36 9,66 9,97
... ... ... ... ... ... ... ... ... ... ... ... ...
16 14 10 15 35 7,05 2,36 4,49 9,08 1,69 4,76 7,49 5,61 16 14 10 15 40 7,08 2,40 4,53 9,16 1,76 4,80 7,53 5,68 16 14 10 15 45 7,11 2,44 4,57 9,24 1,83 4,85 7,57 5,75
... ... ... ... ... ... ... ... ... ... ... ... ...
16 14 10 15 105 7,49 2,92 5,07 10,22 2,69 5,39 8,04 6,59 16 14 10 15 110 7,52 2,96 5,11 10,30 2,76 5,44 8,08 6,66 16 14 10 15 115 7,56 3,00 5,15 10,38 2,83 5,48 8,12 6,73
... ... ... ... ... ... ... ... ... ... ... ... ...
16 14 10 15 200 8,09 3,68 5,85 11,77 4,05 6,26 8,79 7,92 16 14 10 15 205 8,12 3,72 5,89 11,85 4,12 6,30 8,83 7,99 16 14 10 15 210 8,16 3,76 5,93 11,93 4,20 6,35 8,87 8,06
... ... ... ... ... ... ... ... ... ... ... ... ...
16 14 10 15 295 8,69 4,45 6,63 13,32 5,41 7,12 9,55 9,26 16 14 10 15 300 8,72 4,49 6,67 13,40 5,49 7,17 9,59 9,33 16 14 10 15 305 8,75 4,53 6,71 13,48 5,56 7,21 9,63 9,40
... ... ... ... ... ... ... ... ... ... ... ... ...
16 14 10 15 375 9,20 5,09 7,29 14,62 6,56 7,85 10,18 10,38 16 14 10 15 380 9,23 5,13 7,33 14,71 6,63 7,90 10,22 10,45 16 14 10 15 385 9,26 5,17 7,37 14,79 6,70 7,94 10,26 10,52
... ... ... ... ... ... ... ... ... ... ... ... ...
Trên đây là bảng sử dụng để tra khối lƣợng VRR bổ sung và phân hủy, cách sử dụng bảng tra này nhƣ sau:
Các số liệu về khối lƣợng VRR tồn dƣ nếu đƣợc đo đếm từ tháng 1 đến tháng 3 (ứng với thời điểm T1); từ tháng 4 đến tháng 6 (ứng với thời điểm T2); từ tháng 7 đến tháng 9 (ứng với thời điểm T3); từ tháng 10 đến tháng 12 (ứng với thời điểm T4).
Ví dụ: M=100 m3/ha.
+ Nếu đo tại thời điểm T1, giá trị Wtdư=11 tấn/ha, tra giá trị tương ứng khối lượng VRR bổ sung là 4,91 tấn/ha và khối lượng VRR phân hủy là 2,57 tấn/ha.
+ Nếu đo tại thời điểm T2, giá trị Wtdư=9,5 tấn/ha, tra giá trị tương ứng khối lượng VRR bổ sung là 2,51 tấn/ha và khối lượng VRR phân hủy là 4,24 tấn/ha.
Tương tự tra tại thời điểm T3 và T4 tùy theo người sử dụng cần xác định khối lượng VRR tồn dư tại thời điểm nào trong năm.
Nhận xét chung
Các phương trình tại bảng 3.40 và hệ thống bảng tra 3.41 sẽ giúp người sử dụng xác định đƣợc khối lƣợng VRR bổ sung và phân hủy tại từng thời điểm khác nhau trong năm. Tuy nhiên, phương pháp này phức tạp hơn vì ngoài việc xác định khối lƣợng VRR tồn dƣ (Wtdƣ) cần phải xác định thêm các nhân tố khác nhƣ: trữ lƣợng, đường kính tán, chiều dài tán. Nhưng theo phương pháp này thì độ chính xác sẽ cao hơn vì hạn chế đƣợc những sai số do biến động của khối lƣợng VRR tồn dƣ gây nên và phản ánh tốt sự biến động của các trạng thái thảm thực vật. Vì vậy, tùy theo mức độ chính xác yêu cầu và mục đích của người sử dụng mà có thể sử dụng một trong hai phương pháp xác định khối lượng VRR như đã nêu ở trên.