Lăm học HÀM Đ ộ THON VÀ SẢN LƯỢNG THÂN CÂY TRÀM Ở KHU vực THẠNH HÓA THUỘC TỈNH LONG AN Nguyễn Văn Thêm Hội Khoa học & Kỹ thuật Lâm nghiệp TP Hồ Chí Minh https://doi.Org/10.55250/jo.vnuf.2022.4.055-064 TĨM TẮT Các hàm độ thon thân sản lượng gỗ công cụ để xây dựng biểu thể tích sản lượng rừng Trong nghiên cứu này, hàm độ thon thân tràm xây dựng từ 104 mẫu cấp đường kính từ -ỉ-16 cm Hàm độ thon thích hợp kiểm định từ 12 hàm dự tuyển Hàm sản lượng gỗ kiểm định từ hàm dự tuyến Các hàm độ thon sản lượng gỗ thích họp chọn theo tiêu chuẩn tổng sai lệch bình phương nhỏ Kết nghiên cứu cho thấy hàm Them 2022-1 Them 2022-2 hai hàm thích hợp để ước lượng độ thon hình dạng thân tràm Thể tích thân đứng thể tích gỗ sản phẩm mức cá thể rừng tràm xây dựng hàm V = a + b (D2H) + c (D tì) So với thể tích thân vỏ (100%), thể tích thân khơng vỏ, thể tích gỗ sản phẩm vỏ khơng vỏ, thể tích vỏ tương ứng 66,0%, 99,6%, 62,7% 39,0% Từ khoá: Độ thon thân, hàm độ thon thân, hàm sàn lượng gỗ, hình dạng thân, rừng tràm ĐẶT VẤN ĐỀ Độ thon thân gỗ biểu thị suy giảm đường kính thân từ gốc đến Hình dạng thân đường cong biểu diễn độ thon thân từ gốc đến Các nhà lâm học sử dụng hàm độ thon thân để phân tích ảnh hưởng lập địa phương thức lâm sinh đến gỗ quần thụ (Nguyễn Văn Thêm, 2002) Trong điều tra rừng, hàm độ thon thân sử dụng để ước lượng không đường kính vị trí khác thân, m cịn thể tích phân đoạn gồ với chiều dài bất kỳ, thể tích thân dửng, thể tích gỗ sản phẩm sinh khối thân (Vũ Tiến Hĩnh, 2012) Hàm thể tích thân đứng hàm sản lượng gỗ công cụ để xây dựng biểu thể tích biểu sản lượng rừng (Ngọc Lung Đào Công Khanh, 1999; Vũ Tiến Hỉnh, 2012) Trên giới, nhiều tác giả (Muhainve, 1999; Kozak, 2004; Sharma and Zhang, 2004; Fonweban et al., 2011; Lee et al., 2017; Tang et al., 2017) phát triển hàm độ thon hàm thể tích cho lồi gỗ khác Ở Việt Nam, Ngọc Lung Đào Cơng Khanh (1999) xây dựng hàm thể tích thông ba (Pỉnus keysỉa ex Royle Gordon) khu vực Tây Nguyên Vũ Tiến Hĩnh (2012) xây hàm thể tích hàm độ thon nhiều loài gỗ tự nhiên Việt Nam Rừng tràm (Melaleuca caịuputỉ Powell) đóng vai trị to lớn kinh tế, mơi trường quốc phịng Việc thống kê trữ lượng gỗ đứng trữ lượng gỗ sản phẩm loại rừng đòi hỏi phải có biểu thể tích biểu gỗ sản phẩm Thế nhưng, thiếu hàm thể tích hàm độ thon thân tràm Điều gây khó khăn cho điều tra rừng quản lý rừng Bởi hình dạng độ thon thân thay đổi tùy theo loài gỗ, nên xây dựng hàm độ thon thân tràm từ hàm sẵn có phát triển hàm vấn đề cần đặt Các hàm độ thon phải đảm bảo yêu cầu cải thiện sai số ước lượng độ thon Mặt khác, độ xác hàm độ thon không phụ thuộc vào dung lượng mẫu, dạng hàm số lượng biến dự đốn, mà cịn vào phương pháp xây dựng hàm hồi quy Nghiên cứu tập trung trả lời hai câu hỏi Một độ thon hình dạng thân tràm ước lượng dự đoán dạng hàm nào? Hai sản lượng gỗ thân tràm ước lượng dự đoán dạng hàm nào? Hai câu hỏi trả lời cách kiểm định sai lệch hàm dự tuyển Xuất phát từ nhũng vấn đề đặt đây, mục tiêu nghiên cứu xây dựng hàm độ thon thể tích để thống kê trữ lượng gỗ rừng tràm PHƯƠNG PHAP NGHIÊN cứu 2.1 Đối tượng vị trí khu vực nghiên cứu Đối tượng nghiên cứu rừng tràm ưồng loài từ - 10 tuổi Mật độ ừồng rừng TẠP CHÍ KHOA HỌC VÀ CƠNG NGHỆ LÂM NGHIỆP SỐ - 2022 55 Lăm học ban đầu 10.000 - 20.000 cây/ha số liệu xác 0,1 cm Đoạn đo chiều dài (Ln, độ thon thân ừàm thu thập khu cm) đường kính đáy (D n , cm) Thể tích gỗ vực Thạnh Hóa thuộc tỉnh Long An Khu vực sản phẩm đo từ gốc đến D cv > cm nghiên cứu mang đặc tính chung khí hậu Bảng tổng họp D, H, Dcv, Dov đường nhiệt đới gió mùa Hàng năm khí hậu phân kính vỏ (Kvo, cm) trung bình 104 mẫu chia thành hai mùa mưa khô rõ rệt Mùa sử dụng để xây dựng hàm độ thon thân mưa kéo dài tháng từ tháng đến tháng 11, cịn mùa khơ từ tháng 12 năm trước đến tháng 2.3 Phương pháp xử lý số liệu năm sau Lượng mưa trung bình năm Từ 109 mẫu, sử dụng 104 để xây 2.000,0 mm Độ ẩm khơng khí trung bình 80% dựng hàm Dhcv Dhov, lại Độ cao từ 1,5 - 2,0 m so với mặt nước biển mẫu cấp D = - 18 cm sử dụng để Đất phèn phát triển phù sa kiểm tra khả ứng dụng hàm độ thon Các hàm độ thon thân tràm 2.2 Phương pháp thu thập số liệu Độ thon thân vỏ (D hcv, cm) độ thon kiểm định từ 12 hàm dự tuyển (Bảng 2); hàm + tác giả đề xuất nghiên thân không vỏ (Dhov, cm) mức cá thể cứu Các hàm + 12 ký hiệu tương rừng tràm từ - 10 tuổi (A, năm) xác ứng Them22-1, Them22-2, Them22-3, định từ 109 mẫu (N, cây) thuộc cấp D = + 16 cm Những mẫu chặt hạ cách Them22-4, Them22-5, Kozak88, Kozak01, Kozak02, Muhairwe99, Lee03, Sharmamặt đất 30 cm Tất mẫu xác định chiều dài toàn thân (H, m) thước dây với Zhang04, Sharma-Parton09 Ở hàm - 12, D độ xác 1,0 cm Thân phân chia (cm) = đường kính thân ngang ngực; Dh (cm) thành phân đoạn với chiều dài (L) 100 cm; = đường kính vị trí khác đoạn có chiều dài (Ln) 100 cm thân; H (m) = chiều cao toàn thân; h (m) = Đường kính thân vỏ đường kính thân chiều cao từ gốc đến vị trí khác khơng vỏ vị trí ngang ngực, đường kính thân; Y = (h/H); X i = (1 - Y 1/2)/(l - p 1/2); thân vỏ đường kính thân khơng vỏ đầu X j = (1 - Y 1/4)/(l - 0,011/4); Xk = Q /(l - p 1/3); Q nhỏ (tương ứng D ic v D iov) đầu lớn = - Y 1/3; p = điển uốn - (1,3/H); u = (H h)/(H -1,3); z = (h/.l,3) (tương ứng D 2CV D 2ov) phân đoạn xác định^bằng thước kẹp với độ Bảng Đặc trưng thống kê đếi với mẫu 56 A (năm) N (cây) (1) 10 A (năm) 10 D (cm) Min Max (2) (3) (4) (5) (6) 26 28 17 12 10 N (cây) 26 28 17 12 10 3,7 4,7 6,4 8,2 10,0 11,9 14,4 H(m) 6,6 7,9 9,2 10,5 13,0 14,9 14,7 3,9 5,3 7,4 9,0 10,6 12,8 16,0 Max 6,9 9,6 11,6 15,6 15,6 16,3 16,1 0,18 0,35 0,58 0,51 0,49 0,68 0,97 SEE 0,2 1,0 1,1 1,9 2,1 1,0 0,9 3,4 4,0 5,4 7,4 9,2 11,0 13,0 Min 6,5 6,3 5,4 7,2 10,1 11,9 13,8 SEE cv% (7) 5,0 7,4 9,0 6,2 4,9 5,7 6,7 cv% TẠP CHÍ KHOA HỌC VÀ CƠNG NGHỆ LÂM NGHIỆP SỐ - 2022 2,5 12,3 11,8 18,1 16,3 6,8 5,9 Lâm học TT Bảng Những hàm độ thon dự tuyển nghiên cứu Các hàm độ thon dự tuyển: Tác giả (D3H)a2(a3 - V Ỹ ) (a4Y + a5Ỵ A + a' u + a7(D /H » Them (2022-1) Dh = Dh = ai(D2H)a2(a3 - V Ỹ ) (a4Y + a5ỴA2 + a6 U + a7 (D /H )) Them (2022-2) Dh = ai(D3H)a2(a3 - Them (2022-3) Dh = ai(D2*H)a2(a3 - Dh = ai(D3H)a2(a3 - V Ỹ ) (a4Y + a5YA2 + a6/exp(D/H) + a7(D /H » Dh = a iD a2a 3DX (4ỴA2 + a5Ln(Y + ° ’001) + a6YAl/2 + a7exp(Y) + a8(D/H)) V Ỹ ) (a4Y + a5Ỵ A + a u + a e x P(D /H )) V Ỹ ) (a4Y + a5yA2 + a u + a7exP(D/H)) D ^ X / 33 + a4/exP(Z) + a5D AX J + a6X JAZ) Dh = 10 11 12 Dh = ãiD^ạ - V Ỹ ) (a3ỴA2 + a4Y + a5) Dh2 = aiD2UZ(2 ■ +a3Y + a4Ỵ A 2)) (Dh/D) = aiUY(a2 + a3Y + a4Ỵ A 2> ịu a iD a2H a3X k(a4YA4 + a5/exP