Bài viết Những hàm thể tích đối với thân cây Tràm (Melaleuca cajuputi Powell) ở khu vực Tây Nam Bộ được nghiên cứu với mục tiêu là xây dựng các hàm thể tích thân cây đứng và thể tích gỗ sản phẩm ở múc các thể của rừng Tràm.
9 Trường Đại học Nơng Lâm TP Hồ Chí Minh Stem volume functions for Melaleuca cajuputi trees in Southwestern Region Them V Nguyen Forest Science and Technology Association of Ho Chi Minh City, Vietnam ARTICLE INFO ABSTRACT Research Paper The objective of this study was to construct standing stem and commercial timber volume functions at the individual tree level of Melaleuca cajuputi plantations The insidebark and outsidebark stem volume functions were constructed from 56 sample trees at the diameter classes of to 16 cm The appropriate volume functions are tested from five candidate functions The applicability of the volume functions were tested from 10 trees that were not included in the construction of the volume functions The research results showed that the function V = a + b(D2 * H)c was a suitable function to build the insidebark stem volume function of Melaleuca cajuputi The function V = a + b(D2 H) + c(Dd He ) was a suitable function to build the outsidebark stem and commercial timber volume functions of Melaleuca cajuputi The volume functions gave errors less than 5.0% Compared with the insidebark stem volume, the average ratio for the outsidebark stem volume, insidebark and outsidebark commercial timber volume, and bark volume was 65.7%, 95.2%, 60.6% and 34.3%, respectively Received: November 19, 2021 Revised: April 27, 2022 Accepted: April 28, 2022 Keywords Bark volume Insidebark stem volume Outsidebark stem volume Melaleuca cajuputi plantation Stem volume function Corresponding author Nguyen Van Them Email: nvthem2009@gmail.com Cited as: Nguyen, T V (2022) Stem volume functions for Melaleuca cajuputi trees in Southwestern Region The Journal of Agriculture and Development 21(2), 9-16 www.jad.hcmuaf.edu.vn Tạp chí Nơng nghiệp Phát triển 21(2) 10 Trường Đại học Nơng Lâm TP Hồ Chí Minh Những hàm thể tích thân Tràm (Melaleuca cajuputi Powell) khu vực Tây Nam Bộ Nguyễn Văn Thêm Hội Khoa Học Kỹ Thuật Lâm Nghiệp TP.HCM, TP Hồ Chí Minh THƠNG TIN BÀI BÁO TĨM TẮT Bài báo khoa học Mục tiêu nghiên cứu xây dựng hàm thể tích thân đứng thể tích gỗ sản phẩm mức cá thể rừng Tràm Các hàm thể tích xây dựng từ 56 mẫu cấp đường kính từ đến 16 cm Các hàm thể tích thích hợp kiểm định từ hàm dự tuyển Khả ứng dụng hàm thể tích kiểm định từ 10 khơng tham gia xây dựng hàm thể tích Kết nghiên cứu cho thấy, hàm V = a + b(D2 * H)c hàm thích hợp để xây dựng hàm thể tích thân vỏ mức cá thể Hàm V = a + b(D2 H) + c(Dd He ) hàm thích hợp để xây dựng hàm thể tích thân khơng vỏ thể tích gỗ sản phẩm mức cá thể Các hàm thể tích nhận sai số nhỏ 5,0% So với thể tích thân vỏ, tỷ lệ trung bình thể tích thân khơng vỏ, thể tích gỗ sản phẩm vỏ không vỏ, tỷ lệ vỏ tương ứng 65,7%, 95,2%, 60,6% 34,3% Ngày nhận: 19/11/2021 Ngày chỉnh sửa: 27/04/2022 Ngày chấp nhận: 28/04/2022 Từ khóa Hàm thể tích thân Rừng Tràm Thể tích thân vỏ Thể tích thân khơng vỏ Thể tích vỏ Tác giả liên hệ Nguyễn Văn Thêm Email: nvthem2009@gmail.com Đặt Vấn Đề Ước lượng xác thể tích thân nhiệm vụ quan trọng lâm học điều tra rừng (Nguyen, 2002; Nguyen, 2005; Vu, 2012) Trữ lượng gỗ rừng trồng loài đồng tuổi ướng lượng cách kết hợp hàm thể tích thân với hàm mật độ quần thụ Hàm thể tích thân kết hợp với tỷ trọng gỗ sở cho việc xác định sinh khối gỗ rừng Thể tích thân gỗ thay đổi khơng theo đườgn kính, chiều cao hình dạng thân, mà cịn theo điều kiện khí hậu, địa hình, đất hoạt động lâm sinh Trong điều tra rừng, hàm độ thon (Dh, cm) thể tích thân (V, m3 ) đượ xây dựng theo hàm với nhiều biến dự đoán khác Khi thống kê tài nguyên gỗ phạm vi không gian hẹp, hàm thể tích thân lồi gỗ nhóm lồi gỗ bao gồm yếu tố dự đốn, thường đường kính thân ngang ngực (D, cm) Khi thống kê tài nguyên gõ phạm Tạp chí Nơng nghiệp Phát triển 21(2) vi khơng gian rộng lớn yêu cầu độ xác cao, hàm thể tích bao gồm hai yếu tố dự đốn, thường D chiều cao tồn thân (H, m) Thể tích thân thay đổi theo hình dạng thân (F) Vì thế, hàm thể tích bao gồm ba yếu tố dự đoán (D, H, F) Một số tác giả (Sherrill & ctv., 2011; Adekunle & ctv., 2013; Lee & ctv., 2017) xây dựng hàm thể tích thân gỗ theo ba biến dự đoán (Tuổi = A, năm; D H) Nói chung, xây dựng hàm thể tích thây theo hay nhiều biến dự đoán yêu cầu lý thuyết thực hành Trước nhiều nhà lâm học điều tra rừng xây dựng hàm thể tích loài gỗ khác Việt Nam (Nguyen & Dao, 1988; Nguyen, 2005; Vu, 2012) Hiện cịn thiếu hàm ước lượng thể tích thân vỏ thể tích thân khơng vỏ mức cá thể rừng Tràm khu vực Tây Nam Bộ Mặt khác, độ xác hàm thể tích khơng phụ thuộc vào dạng hàm số lượng biến dự đốn, mà cịn vào phương pháp xây dựng www.jad.hcmuaf.edu.vn Trường Đại học Nông Lâm TP Hồ Chí Minh 11 hàm hồi quy Xuất phát từ vấn đề đặt dựng hàm thể tích đây, mục tiêu nghiên cứu xây dựng hàm thể tích thân đứng thể tích 2.3 Phương pháp xử lý số liệu gỗ sản phẩm múc thể rừng Tràm Các Các hàm thể tích (VCV , VOV , VSPCV hàm thể tích xây dựng dựa giả thuyết VSPOV , m3 /cây) xây dựng từ 56 mẫu, V = f(D, H, D2 * H, Db *Hb ) lại 10 mẫu sử dụng để kiểm tra Vật Liệu Phương Pháp Nghiên Cứu khả ứng dụng hàm thể tích Hai đại lượng VCV VOV xác 2.1 Đối tượng nghiên cứu định theo công thức 1; Vj = thể tích j, Vij = thể tích phân đoạn i j, Vnj = Đối tượng nghiên cứu rừng Tràm trồng thể tích đoạn j Thể tích phân lồi Mật độ trồng rừng ban đầu 10.000 - đoạn (Vij ) xác định theo cơng thức hình nón 20.000 cây/ha Số liệu thể tích thân Tràm cụt (Cơng thức 2); L = 100 cm (Chiều thu thập khu vực Thạnh Hóa thuộc tỉnh dài phân đoạn), D1 (cm) = đường kính Long An Khu vực nghiên cứu mang đặc tính đầu lớn vỏ không vỏ phân đoạn, chung khí hậu nhiệt đới gió mùa Hàng năm D2 (cm) = đường kính đầu nhỏ vỏ khơng khí hậu phân chia thành hai mùa mưa khơ vỏ phân đoạn Thể tích đoạn rõ rệt Mùa mưa kéo dài tháng từ tháng đến (Vnj ) tính theo cơng thức hình nón tháng 11, cịn mùa khơ từ tháng 12 năm trước đến (Công thức 3); gn (tiết diện đáy đoạn tháng năm sau Lượng mưa trung bình năm ngọn) = 0,00007854 D2n với Dn = đường kính đáy 2.000 mm Độ ẩm khơng khí trung bình 80% Độ đoạn vỏ không vỏ, Ln = chiều dài cao từ 1,5 - 2,0 m so với mặt biển Đất phèn phát đoạn Sau xây dựng hàm VCV , VOV , triển phù sa VSPCV VSPOV theo hàm dự tuyển (4) - (8); V (m3 ) biến phụ thuộc, biến dự 2.2 Phương pháp thu thập số liệu đoán D H Các hàm nhiều tác giả sử dụng ( Sherrill & ctv., 2011; Vu, 2012; Giả định hình dạng (F) thể tích (V, m3 ) Lee & ctv., 2017) Tràm có kích thước đường kính (D, cm) chiều cao (H, m) Theo VJ = Vij + Vnj (1) giả định này, thể tích thân đứng vỏ (VCV , m3 ) không vỏ (VOV , m3 ), thể tích gỗ sản phẩm vỏ (VSPCV , m3 ) không vỏ (VSPOV , m3 ) πL (D2 + D22 + D1 D2 ) (2) Vij (m3 ) = Tràm xác định từ 66 mẫu cấp D = ∗ ∗ 10000 ö 16 cm Những mẫu chặt hạ vị trí 30 cm kể từ gốc Tất mẫu xác Vnj = (1/3)gnj Ln (3) định chiều dài toàn thân (H, m) thước dây với độ xác 1,0 cm Thân phân chia thành phân đoạn với chiều dài (L) 100 cm; V = a(D2 H) (4) đoạn có chiều dài (Ln ) 100 cm Đường kính thân vỏ (DCV , cm) khơng vỏ V = a(D2 H)b (5) (DOV , cm) vị trí ngang ngực, đường kính đầu lớn vỏ không vỏ (tương ứng D1CV D1OV ), đường kính đầu nhỏ vỏ khơng vỏ (tương ứng V = a + b(D2 H) (6) D2CV , D2OV ) phân đoạn xác định thước Panme với độ xác 0,1 mm Hiệu V = a + b(D2 H)c (7) số DCV DOV đường kính vỏ (DVo , cm), cịn (DCVo /2) độ dày vỏ Đoạn đo chiều dài (Ln, cm) đường kính đáy (Dn , cm) V = a + b(D2 H) + c(Dd He ) (8) Thể tích gỗ sản phẩm vỏ khơng vỏ (tương ứng VSPCV , VSPOV , m3 ) đo từ gốc đến DCV Các hệ số hồi quy thống kê sai lệch ≥ cm Bảng tổng hợp đặc trưng thống kê kích hàm dự tuyển (4) - (8) xác định thước mẫu sử dụng để xây phương pháp hồi quy tương quan phi tuyến tính www.jad.hcmuaf.edu.vn Tạp chí Nơng nghiệp Phát triển 21(2) 12 Trường Đại học Nơng Lâm TP Hồ Chí Minh Bảng Đặc trưng thống kê kích thước mẫu rừng Tràm Chỉ tiêu (1) D (cm) H (m) VCV (m3 ) VOV (m3 ) VSPCV (m3 ) VSPOV (m3 ) N (cây) (2) 56 56 56 56 56 56 Trung bình (3) 7,9 10,0 0,02944 0,01890 0,02879 0,01847 Marquartz Mức độ chặt chẽ mối quan hệ V với D H đánh giá theo hệ số xác định (R2 ; Công thức 9) Sai lệch hàm hồi quy so với số liệu thực tế đánh giá theo tiêu chuẩn: (1) Tổng sai lệch bình phương (SSR; Cơng thức 10); (2) Sai số chuẩn ước lượng (SEE; Công thức 11); (4) Sai số tuyệt đối trung bình (MAE; Cơng thức 12) sai số tuyệt đối trung bình theo phần trăm (MAPE; Cơng thức 13); (6) Sai lệch trung bình (ME; Cơng thức 14) sai lệch trung bình theo phần trăm (MPE; Công thức 15) Ở công thức (9) – (15), Vi VJ = giá trị biến V thực nghiệm ước lượng; VBq = giá trị trung bình biến V thực nghiệm; n = dung lượng quan sát; p = số lượng tham số hàm thể tích Hệ số R2 cho biết mức độ chặt chẽ mối quan hệ V với biến dự đoán D H Tiêu chuẩn SSR sử dụng để đánh giá tổng sai lệch hàm quy so với số liệu thực tế Các tiêu chuẩn SEE, MAE, MAPE, ME MPE sử dụng để đánh giá sai số hàm hồi quy Mục đích phân tích hồi quy xác định hàm ước lượng V với sai lệch nhỏ Theo mục đích này, trước hết phân tích so sánh hệ số R2 thống kê sai lệch (SSR, SEE, MAE, MAPE, ME, MPE) hàm dự tuyển Sau chọn hàm thể tích thích hợp theo tiêu chuẩn SSRMin Khả ứng dụng hàm thể tích kiểm định từ 10 mẫu; sai lệch hàm thể tích so với thể tích thực tế đánh giá theo tiêu chuẩn MPE Các bước phân tích hồi quy tương quan thực phần mềm thống kê STATGRAPHICS Centurion XV.I 15.1.02 n (Vi − VJ ) i=1 ∗ 100 (9) R2 = − n 2 (Vi − Vbq ) i=1 n (Vi − VJ ) SSR = i=1 Tạp chí Nơng nghiệp Phát triển 21(2) Min (4) 3,7 6,3 0,00435 0,00268 0,00300 0,00179 Max (5) 16,0 16,3 0,10871 0,07606 0,10820 0,07574 ± SEE (6) 3,1 2,8 0,02806 0,02019 0,02811 0,02023 SEE = MAE = MAPE = CV(%) (7) 39,1 28,4 95,3 106,8 97,7 109,5 SSR n-p Vi − VJ n MAE ∗ 100 Vi ME = (Vi − VJ ) MPE = Vi − V J ∗ 100 Vi (11) (12) (13) (14) (15) Kết Quả Thảo Luận 3.1 Hàm thể tích thân Tràm Các hàm ước lượng VCV (Hàm 4a ö 8a) VOV (Hàm 4b ö 8b) Tràm ghi lại Bảng Các hàm VCV VOV tồn mức ý nghĩa cao (P < 0,01) Hệ số xác đinh (R2 ) nhận giá trị cao hàm 7a 8b (tương ứng R2 = 96,61% 92,96%), thấp hàm 4a (R2 = 94,92%) 4b (R2 = 92,59%) Đối với VCV , giá trị MAPEMin MAPEMax tương ứng hàm 8a 4a (tương ứng 10,6% 17,8%) Đối với VOV , giá trị MAPEMin MAPEMax tương ứng hàm 6b 7b (tương ứng 15,8% 20,4%) Tất hàm ước lượng VCV nhận sai số hệ thống (ME); hàm 5a 8a mang giá trị âm, hàm khác nhận giá trị dương Đối với hàm ước lượng VOV , giá trị ME hàm 5b mang giá trị âm, hàm khác nhận giá trị dương Khi sử dụng hàm dự tuyển để ước lượng VCV , hàm 7a nhận SSRMin (0,00147), SSRMax hàm 4a (0,0022) Khi sử dụng hàm dự tuyển để ước lượng VOV , hàm 8b nhận SSRMin (10) (0,00158), cịn SSRMax hàm 4b (0,00166) Từ phân tích tương quan sai lệch www.jad.hcmuaf.edu.vn 13 Trường Đại học Nơng Lâm TP Hồ Chí Minh Bảng Các hàm ước lượng thể tích thân vỏ Tràm Thống kê (1) a b c d e Pα R2 SSR SEE MAE MAPE ME MPE (4a) (2) 0,000032 < 0,01 94,92 0,00220 0,00632 0,00407 17,8 0,00205 14,8 Các hàm thể tích thân vỏ (5a) (6a) (7a) (3) (4) (5) 0,000093 0,003845 -0,002348 0,860660 0,000030 0,000142 0,811226 < 0,01 96,55 0,00150 0,00526 0,00322 10,7 - 0,00021 -3,4 < 0,01 95,94 0,00176 0,00571 0,00332 12,9 0,00000 -6,1 < 0,01 96,61 0,00147 0,00527 0,00337 12,8 0,00000 0,6 (8a) (6) 0,000865 0,000016 0,000022 0,837038 1,992220 < 0,01 96,50 0,00152 0,00530 0,00314 10,6 -0,00012 -3,2 Bảng Các hàm ước lượng thể tích thân khơng vỏ Tràm Các hàm thể tích thân vỏ (5b) (6b) (7b) (3) (4) (5) 0,000034 0,000799 -0,001550 0,000021 0,000053 0,887867 < 0,01 92,59 0,00166 0,00550 0,00344 17,9 0,00042 3,7 < 0,01 92,83 0,00161 0,00546 0,00341 16,7 -0,00019 -4,8 hàm dự tuyển cho thấy, theo tiêu chuẩn SRMin , hàm 7a (Bảng 2; Hình 1) hàm 8b (Bảng 3; Hình 1) hai hàm thích hợp để ước lượng VCV VOV theo D H Tràm 3.2 Hàm thể tích gỗ sản phẩm Các hệ số hồi quy thống kê sai lệch hàm ước lượng VSPCV VSPOV Tràm dẫn Bảng Tất hàm ước lượng VSPCV (Hàm 4c ö 8c; Bảng 4) VSPOV (Hàm 4d ö 8d; Bảng 5) tồn mức ý nghĩa cao (P < 0,01) Hệ số xác định (R2 ) nhận giá trị cao hàm 8c 8d (tương ứng R2 = 96,82% 93,13%), thấp hàm 4c www.jad.hcmuaf.edu.vn < 0,01 92,67 0,00164 0,00552 0,00335 15,8 0,00000 -3,8 < 0,01 92,90 0,00159 0,00548 0,00357 20,4 0,00000 0,9 (8b) (6) 0,000372 0,000009 0,000004 0,842667 2,568450 < 0,01 92,96 0,00158 0,00541 0,00334 15,9 0,00000 -3,3 VCV = f(D, H) VOV = f(D, H) VOV (m3 /cây) (1) a b c d e Pα R2 SSR SEE MAE MAPE ME MPE (4b) (2) 0,000021 0,940304 VCV (m3 /cây) Thống kê (a) (b) Hình Đồ thị biểu diễn mối quan hệ thể tích thân vỏ (a) khơng vỏ (b) với đường kính ngang ngực chiều cao toàn thân Tràm 4d (tương ứng R2 = 95,24% 92,67%) Đối với VSPCV , giá trị MAPEMin MAPEMax tương Tạp chí Nơng nghiệp Phát triển 21(2) 14 Trường Đại học Nông Lâm TP Hồ Chí Minh Bảng Các hàm ước lượng thể tích gỗ sản phẩm vỏ Tràm Thống kê (1) a b c d e R2 Pα SSR SEE MAE MAPE ME MPE (4c) (2) 0,000031 95,24 < 0,01 0,00207 0,00614 0,00374 15,0 0,00167 11,3 Các hàm thể tích thân vỏ (5c) (6c) (7c) (3) (4) (5) 0,000082 0,003149 -0,003273 0,875241 0,000030 0,000148 0,805999 96,50 < 0,01 0,00152 0,00531 0,00325 12,6 -0,00031 -6,1 95,91 < 0,01 0,00178 0,00573 0,00336 15,2 0,00000 -8,0 96,62 < 0,01 0,00147 0,00527 0,00334 13,2 0,00000 0,7 (8c) (6) -0,003243 -0,000019 0,000063 1,182740 1,716990 96,82 < 0,01 0,00138 0,00506 0,00327 13,4 0,00001 1,3 Bảng Các hàm ước lượng thể tích gỗ sản phẩm không vỏ Tràm Thống kê (1) a b c d e R2 Pα SSR SEE MAE MAPE ME MPE (4d) (2) 0,000021 92,67 < 0,01 0,00165 0,00547 0,00334 16,4 0,00018 -1,26 Các hàm thể tích thân vỏ (5d) (6d) (7d) (3) (4) (5) 0,000030 0,000337 -0,002136 0,955370 0,000021 0,000055 0,882869 92,81 < 0,01 0,00162 0,00547 0,00340 18,0 -0,00027 -8,07 92,69 < 0,01 0,00165 0,00552 0,00335 17,2 0,00000 -4,89 92,94 < 0,01 0,00159 0,00547 0,00356 21,6 0,00000 1,50 (8d) (6) -0,001376 0,000006 0,000008 0,705896 2,511580 93,13 < 0,01 0,00155 0,00535 0,00343 19,0 -0,00004 -0,52 tương ứng hàm 4d 7d (tương ứng 16,4% 21,6%) Tất hàm ước lượng VSPCV nhận sai số hệ thống (ME); hàm 5c mang giá trị âm, hàm khác mang giá trị dương Đối với hàm ước lượng VSPOV , giá trị ME hàm 5d 8d mang giá trị âm, hàm khác nhận giá trị dương Khi ước lượng VSPCV , hàm 8c nhận SSRMin (0,00138), SSRMax (a) (b) hàm 4c (0,00207) Khi ước lượng VSPOV , hàm Hình Đồ thị biểu diễn mối quan hệ thể tích 8d nhận SSRMin (0,00155), cịn SSRMax hàm gỗ sản phẩm vỏ (a) khơng vỏ (b) với đường kính 4d (0,00165) Từ phân tích tương quan ngang ngực chiều cao toàn thân Tràm sai lệch hàm dự tuyển cho thấy, theo tiêu chuẩn SRMin , hàm 8c (Bảng 4; Hình 2) hàm 8d (Bảng 5; Hình 2) hai hàm thích hợp để ước ứng hàm 5c 6c (tương ứng 12,6% 15,2%) lượng V SPCV VSPOV theo D H Đối với VSPOV , giá trị MAPEMin MAPEMax VSPOV = f(D, H) VSPOV (m3 /cây) VSPCV (m3 /cây) VSPCV = f(D, H) Tạp chí Nơng nghiệp Phát triển 21(2) www.jad.hcmuaf.edu.vn 15 Trường Đại học Nông Lâm TP Hồ Chí Minh Bảng Kiểm định sai lệch hàm thể tích thân vỏ thể tích thân khơng vỏ TT (1) 10 D (cm) H (m) (2) (3) 4,5 6,8 5,4 8,5 6,5 8,0 7,5 7,5 8,0 9,8 8,5 8,5 10,5 10,4 12,5 12,7 14,5 13,5 15,5 14,5 Trung bình Thể tích thân Thực tế Hàm 7a (4) (5) 0,0060 0,0054 0,0109 0,0101 0,0124 0,0136 0,0159 0,0168 0,0205 0,0240 0,0229 0,0236 0,0369 0,0407 0,0687 0,0648 0,0963 0,0874 0,0947 0,1036 vỏ MPE (6) 10,4 7,4 -10,1 -5,7 -17,4 -3,0 -10,4 5,6 9,2 -9,5 -2,3% Thể tích Thực tế (7) 0,0033 0,0061 0,0068 0,0090 0,0138 0,0131 0,0226 0,0508 0,0729 0,0617 thân không vỏ Hàm 8b MPE (8) (9) 0,0034 -4,1 0,0064 -5,3 0,0072 -5,1 0,0077 14,1 0,0136 1,1 0,0114 13,4 0,0217 4,0 0,0397 21,9 0,0543 25,5 0,0678 -10,0 5,5% Bảng Thể tích thân trung bình mức cá thể rừng Tràm D (cm) 10 12 14 16 H (m) 6,6 8,4 10,2 11,8 13,4 14,8 16,0 VCV (m3 ) 0,00386 0,01222 0,02481 0,04183 0,06326 0,08890 0,11843 VOV (m3 ) 0,00286 0,00710 0,01458 0,02599 0,04162 0,06143 0,08507 VSPCV (m3 ) 0,00307 0,01134 0,02413 0,04163 0,06339 0,08852 0,11585 VSPOV (m3 ) 0,00172 0,00644 0,01444 0,02625 0,04196 0,06134 0,08387 Vo (m3 ) 0,00100 0,00513 0,01023 0,01585 0,02164 0,02747 0,03337 Tràm ứng tốt tiêu chuẩn SSRMin MPEMin Vì thế, bốn hàm sử dụng để ước lượng VCV , 3.3 Kiểm định khả ứng dụng VOV , VSPCV vaà VSPOV mức cá thể hàm thể tích thân Tràm rừng Tràm; cấp D = - 16 cm Từ phân tích mục 3.1 cho thấy, hai hàm (7a) (8b) nhận SSRMin MPEMin Sai lệch VCV từ hàm (7a) VOV từ hàm (8b) so với số liệu 10 kiểm tra ghi lại Bảng Ở Bảng 6, hai cột tương ứng đường kính ngang ngực vỏ chiều cao toàn thân; cột tương ứng VCV VOV thực tế; cột tương ứng VCV VOV ước lượng tương ứng từ hàm 7a hàm 8b; cột sai số ước lượng VCV VOV theo phần trăm Từ số liệu Bảng cho thấy, so với số liệu kiểm tra, hàm (7a) nhận sai số âm 6/10 cây, hàm (8b) nhận sai số âm 4/10 Hàm (7a) nhận sai số trung bình -2,3% Hàm (8b) nhận sai số trung bình 5,5% Từ số liệu Bảng - cho thấy, bốn hàm (7a), (8b), (8c) (8d) nhận MPE < 5,0% Thông thường, sai số thống kê thể tích thân gỗ dao động từ ±5 ±10% Từ phân tích cho thấy, bốn hàm (7a), (8b), (8c) (8d) đáp www.jad.hcmuaf.edu.vn 3.4 Biểu thể tích thân thể tích gỗ sản phẩm rừng Tràm Biểu thể tích thân vỏ khơng vỏ trung bình mức cá thể rừng Tràm xây dựng tương ứng từ hàm 7a 8b Biểu thể tích gỗ sản phẩm vỏ khơng vỏ xây dựng tương ứng từ hàm 8c 8d (Bảng 7) Bảng ghi lại tỷ lệ thành phần thể tích thân Tràm Từ cho thấy, so với VCV , tỷ lệ trung bình bốn thành phần (VOV , VSPCV , VSPOV VVo ) tương ứng 65,7%, 95,2%, 60,6% 34,3% Kết Luận Nghiên cứu xây dựng hàm thể tích thân đứng vỏ khơng vỏ, hàm thể tích gỗ sản phẩm vỏ khơng vỏ mức cá thể rừng Tràm Hàm V = a + b(D2 * H)c hàm thích hợp để xây dựng hàm thể tích thân Tạp chí Nơng nghiệp Phát triển 21(2) 16 Trường Đại học Nông Lâm TP Hồ Chí Minh Bảng Tỷ lệ thành phần thể tích thân Tràm D (cm) H (m) 6,6 8,4 10,2 10 11,8 12 13,4 14 14,8 16 16,0 Trung bình VCV (%) 100 100 100 100 100 100 100 100 VOV (%) 74,1 58,1 58,8 62,1 65,8 69,1 71,8 65,7 vỏ mức cá thể rừng Tràm Hàm V = a + b(D2 H) + c(Dd He ) hàm thích hợp để xây dựng hàm thể tích thân khơng vỏ thể tích gỗ sản phẩm mức cá thể rừng Tràm Các hàm thể tích nhận sai số nhỏ 5,0% So với thể tích thân vỏ, tỷ lệ trung bình thể tích thân khơng vỏ, thể tích gỗ sản phẩm vỏ không vỏ, tỷ lệ vỏ tương ứng 65,7%, 95,2%, 60,6% 34,3% Trong thực tế, sử dụng hàm thể tích để ước lượng thể tích thân đứng vỏ khơng vỏ, thể tích gỗ sản phẩm vỏ khơng vỏ mức cá thể rừng Tràm Lời Cam Đoan Tôi cam đoan báo thực khơng có mâu thuẫn kết nghiên cứu với nhà nghiên cứu khác VSPCV (%) 79,5 92,8 97,3 99,5 100,2 99,6 97,8 95,2 VSPOV (%) 44,64 52,68 58,21 62,74 66,33 69,00 70,82 60,6 Vo (%) 25,9 41,9 41,2 37,9 34,2 30,9 28,2 34,3 Nguyen, N L., & Dao, C K (1999) Research on growth and yield of planted forests (Applied to Pinus keysia Royle ex Gordon) forests in Vietnam Ho Chi Minh City, Vietnam: Agricultural Publishing House Nguyen, T B (2005) Scheduling growth and temporary yield of pure Acacia hybrid forest Journal of Agriculture and Rural Development 13, 91-95 Nguyen, V T (2002) Forest ecology Ho Chi Minh City, Vietnam: Agricultural Publishing House Sherrill J R., Bullock, B P., Mullin, T J., McKeand, S E., & Purnell, R C (2011) Total and merchantable stem volume equations for midrotation loblolly pine (Pinus taeda L.) Southern Journal of Applied Forestry 35(3), 105-108 https://doi.org/10 1093/sjaf/35.3.105 Vu, T H (2012) Method of constructing the volume tables for standing trees in natural forests in Vietnam Ho Chi Minh City, Vietnam: Agricultural Publishing House Lời Cám Ơn Tác giả xin chân thành cảm ơn Bộ môn Lâm sinh – Khoa Lâm nghiệp – Trường Đại học Nông Lâm TP.HCM cổ vũ cho tác giả hoàn thành báo Tài Liệu Tham Khảo (References) Adekunle, V A J., Nair, K N., Srivastava, A K., & Singh, N K (2013) Models and form factors for stand volume estimation in natural forest ecosystems: a case study of Katarniaghat Wildlife Sanctuary (KGWS), Bahraich District, India Journal of Forestry Research 24(2), 217-228 https://doi.org/10.1007/ s11676-013-0347-8 Lee, D., Seo, Y., & Choi, J (2017) Estimation and validation of volume equations for Pinus densiflora, Pinus koraiensis, and Larix kaempferi in South Korea Forest Science and Technology 13(2), 77-82 https://doi.org/10.1080/21580103.2017.1315963 Tạp chí Nơng nghiệp Phát triển 21(2) www.jad.hcmuaf.edu.vn ... ) hàm thích hợp để xây dựng hàm thể tích thân khơng vỏ thể tích gỗ sản phẩm mức cá thể Các hàm thể tích nhận sai số nhỏ 5,0% So với thể tích thân vỏ, tỷ lệ trung bình thể tích thân khơng vỏ, thể. .. dựng hàm thể tích lồi gỗ khác Việt Nam (Nguyen & Dao, 1988; Nguyen, 2005; Vu, 2012) Hiện thiếu hàm ước lượng thể tích thân vỏ thể tích thân khơng vỏ mức cá thể rừng Tràm khu vực Tây Nam Bộ Mặt... rừng Tràm Hàm V = a + b(D2 H) + c(Dd He ) hàm thích hợp để xây dựng hàm thể tích thân khơng vỏ thể tích gỗ sản phẩm mức cá thể rừng Tràm Các hàm thể tích nhận sai số nhỏ 5,0% So với thể tích thân