Đánh giá sinh trưởng của một số loài cây trồng nguyên liệu giấy tại công ty lâm nghiệp đoan hùng huyện đoan hùng tỉnh phú thọ

i LỜI CAM ĐOAN Tôi cam đoan, công trình nghiên cứu riêng tơi Các số liệu, kết nêu luận văn trung thực chƣa đƣợc cơng bố cơng trình nghiên cứu khác Nếu nội dung nghiên cứu tơi trùng lặp với cơng trình nghiên cứu cơng bố, tơi xin hồn tồn chịu trách nhiệm tuân thủ kết luận đánh giá luận văn Hội đồng khoa học Hà Nội, ngày 30 tháng 11 năm 2017 Người cam đoan Nguyễn Đức Nghĩa ii LỜI CẢM ƠN Đề tài luận văn “Đánh giá sinh trưởng số loài trồng nguyên liệu giấy Công ty Lâm nghiệp Đoan Hùng huyện Đoan Hùng, tỉnh Phú Thọ” đƣợc hồn thành theo chƣơng trình đào tạo Thạc sỹ trƣờng Đại học Lâm nghiệp Việt Nam Trong trình học tập thực luận văn, tác giả nhận đƣợc quan tâm, giúp đỡ tạo điều kiện thuận lợi Ban giám hiệu Trƣờng Đại học Lâm nghiệp; Phòng đào tạo sau đại học;các thầy giáo, cô giáo Trƣờng Đại học Lâm nghiệp; anh, chị, em, bạn bè đồng nghiệp tạo điều kiện thuận lợi cho thân tác giả trình thực đề tài Nhân dịp này, tác giả xin đƣợc bày tỏ lòng biết ơn sâu sắc trƣớc quan tâm giúp đỡ quý báu Đặc biệt, tác giả xin đƣợc bày tỏ lịng biết ơn chân thành tới TS Trần Việt Hà, ngƣời thầy hƣớng dẫn nhiệt tình, truyền đạt kinh nghiệm quý giá, ý tƣởng nghiên cứu khoa học giúp tác giả hoàn thành luận văn Mặc dù cố gắng nỗ lực, nhƣng kinh nghiệm nghiên cứu chƣa nhiều, đặc biệt hạn chế mặt thời gian trình nghiên cứu nên luận văn khơng tránh khỏi thiếu sót định Tác giả mong nhận đƣợc góp ý thầy cô giáo bạn bè đồng nghiệp luận văn đƣợc hoàn chỉnh Xin chân thành cảm ơn! Hà Nội, ngày 30 tháng 11 năm 2017 Tác giả Nguyễn Đức Nghĩa iii MỤC LỤC Trang Trang phụ bìa LỜI CAM ĐOAN i LỜI CẢM ƠN ii MỤC LỤC iii DANH MỤC CÁC TỪ VIẾT TẮT vi DANH MỤC CÁC BẢNG vii DANH MỤC CÁC HÌNH viii ĐẶT VẤN ĐỀ Chƣơng TỔNG QUAN VẤN ĐỀ NGHIÊN CỨU 1.1 Nghiên cứu giới 1.1.1 Nghiên cứu Keo tai tƣợng 1.1.2 Nghiên cứu Keo lai 1.1.3 Các nghiên cứu Bạch đàn 1.2 Nghiên cứu Việt Nam 1.2.1 Những nghiên cứu trồng rừng nguyên liệu 1.2.2 Những nghiên cứu Keo tai tƣợng 1.2.3 Nghiên cứu Keo lai 1.2.3 Nghiên cứu Bạch đàn 11 1.3 Nhận xét chung 13 Chƣơng MỤC TIÊU, NỘI DUNG, PHƢƠNG PHÁP NGHIÊN CỨU 15 2.1 Mục tiêu nghiên cứu 15 2.2 Đối tƣợng phạm vi nghiên cứu 15 2.3 Nội dung nghiên cứu 15 2.4 Phƣơng pháp nghiên cứu 15 2.4.1 Phƣơng pháp kế thừa 15 iv 2.4.2 Phƣơng pháp điều tra ngoại nghiệp 16 2.4.3 Phƣơng pháp xử lý số liệu 17 Chƣơng ĐIỀU KIỆN TỰ NHIÊN, KINH TẾ - XÃ HỘI KHU VỰC NGHIÊN CỨU 24 3.1 Điều kiện tự nhiên 24 3.1.1 Vị trí địa lý 24 3.1.2 Địa hình địa mạo 24 3.1.3 Khí hậu, thủy văn 24 3.1.4 Đặc điểm đất 25 3.2 Điều kiện kinh tế-xã hội 25 3.2.1 Sản xuất lâm nghiệp 25 3.2.2 Dịch vụ 25 3.2.3 Lao động 25 3.2.4 Cơ sở hạ tầng 26 3.3 Lƣợc sử trồng rừng 26 3.3.1 Loài Keo tai tƣợng , keo lai 26 3.3.2 Bạch đàn 29 Chƣơng KẾT QUẢ NGHIÊN CỨU 32 4.1 Tình hình sinh trƣởng lâm phần 32 4.1.1 Kiểm tra OTC 32 4.1.2 Sinh trƣởng đƣờng kính D1,3 33 4.1.3 Sinh trƣởng chiều cao vút 37 4.2 Một số quy luật kết cấu lâm phần 40 4.2.1 Quy luật phân bố số theo cỡ đƣờng kính 40 4.2.2 Quy luật phân bố số theo cấp chiều cao 42 4.2.3 Tƣơng quan chiều cao Hvn đƣờng kính D1,3 45 4.3 Trữ lƣợng tăng trƣởng lâm phần 47 v 4.4 Đánh giá hiệu kinh tế 48 4.4.1 Xác định chi phí đầu tƣ cho 1ha rừng trồng 49 4.4.2 Xác định thu nhập cho 1ha rừng trồng 50 4.4.3 Xác định hiệu kinh tế cho 1ha trồng rừng 51 4.5 Đề xuất số giải pháp phát triển rừng trồng sản xuất khu vực nghiên cứu 52 KẾT LUẬN, TỒN TẠI, KHUYẾN NGHỊ 54 5.1 Kết luận 54 5.2 Tồn 54 5.3 Khuyến nghị 55 TÀI LIỆU THAM KHẢO vi DANH MỤC CÁC TỪ VIẾT TẮT Nguyên nghĩa Ký hiệu Bộ NN&PTNT D1,3 Bộ Nông nghiệp Phát triển nơng thơn Đƣờng kính thân vị trí 1,3 m Dt Đƣờng kính chiếu tán Fl Phân bố lý thuyết fi Phân bố thực nghiệm Hvn Chiều cao vút Hdc Chiều cao dƣới cành OTC Ô tiêu chuẩn QLSD Quản lý sử dụng SXKD Sản xuất kinh doanh vii DANH MỤC CÁC BẢNG STT Tên bảng Trang 4.1 Kết kiểm tra OTC 32 4.2 Sinh trƣởng đƣờng kính D1,3 bình qn 33 4.3 Tăng trƣởng bình qn đƣờng kính 1,3 36 4.4 Sinh trƣởng chiều cao bình quân 37 4.5 Tăng trƣởng bình quân chiều cao vút 39 4.6 Kết qua kiểm định phân bố N/D1,3 theo phân bố Weibull 40 4.7 Kết qua kiểm định phân bố N/Hvn theo phân bố Weibull 43 4.8 Tƣơng quan H/D hệ số phƣơng trình hồi quy 46 4.9 Phƣơng trình tƣơng quan Hvn/D1,3 47 4.10 Tăng trƣởng thể tích trữ lƣợng lâm phần 47 4.11 Chi phí đầu tƣ cho rừng trồng 49 4.12 Chi phí khai thác vận chuyển (đ/m3) 50 4.13 Thu nhập cho rừng trồng 50 4.14 Tổng hợp tiêu đánh giá hiệu kinh tế 51 viii DANH MỤC CÁC HÌNH STT Tên hình Trang 4.1 Keo tai tƣợng tuổi 34 4.2 Keo lai tuổi 34 4.3 Bạch đàn tuổi 35 4.4 4.5 Biểu đồ sinh trƣởng D1,3 Keo tai tƣợng, Keo lai Bạch đàn tuổi 5; 6; Biểu đồ sinh trƣởng Hvn Keo tai tƣợng, Keo lai Bạch đàn tuổi 5; 6; 35 38 4.6 Biểu đồ phân bố N/D1,3 Keo tai tƣợng 41 4.7 Biểu đồ phân bố N/D1,3 Keo lai 41 4.8 Biểu đồ phân bố N/D1,3 Bạch đàn 42 4.9 Biểu đồ phân bố N/Hvn Keo tai tƣợng 44 4.10 Biểu đồ phân bố N/Hvn Keo lai 44 4.11 Biểu đồ phân bố N/Hvn Bạch đàn 45 ĐẶT VẤN ĐỀ Trong bối cảnh nguồn tài nguyên thiên nhiên bị suy thoái, đặc biệt tài nguyên rừng, năm gần Việt Nam hạn chế mở cửa khu rừng tự nhiên dẫn đến tình trạng lƣợng cung tài ngun khơng đáp ứng đủ nhu cầu sử dụng ngồi nƣớc Từ đó, rừng trồng trở thành đối tƣợng đƣợc quan tâm đầu tƣ phát triển Tập trung phát triển vùng rừng trồng nguyên liệu để cung cấp gỗ lớn, nhỏ cho khu công nghiệp, nhà máy, xƣởng chế biến gỗ, ngồi cịn phục vụ cho nhu cầu nƣớc nƣớc mà mục tiêu, phƣơng hƣớng giải nƣớc ta Tập đoàn trồng rừng sản xuất gồm loài nhƣ Keo tai tƣợng, Keo lai, Bạch đàn, Bồ đề, Thông, Mỡ, Lát hoa, Trám, Xoan, Giổi, Lim xanh, Dẻ, Chò chỉ, Sao đen, Trong đó, Keo Bạch đàn hai lồi chủ lực với tỷ lệ gây trồng 80% diện tích Đây lồi mọc nhanh, sinh trƣởng tốt, phù hợp với nhiều lập địa khác nhau, chu kỳ kinh doanh ngắn – năm – 15 năm tùy mục đích kinh doanh Từ đó, nhiều công ty, địa phƣơng, đơn vị chủ rừng lựa chọn trồng để phát triển trồng rừng kinh tế Tuy nhiên, trồng rừng kinh tế chƣa đƣợc đồng hóa dẫn đến suất chất lƣợng rừng thu đƣợc chƣa cao Đây vấn đề xảy nhiều nơi cần có nghiên cứu, đề xuất nhằm giải tình trạng Cơng ty Lâm nghiệp Đoan Hùng nằm địa bàn huyện Đoan Hùng tỉnh Phú Thọ đơn vị hạch toán Tổng công ty Giấy Việt Nam; đƣợc giao nhiệm vụ chủ yếu thực trồng, chăm sóc, bảo vệ, khai thác loài trồng rừng nguyên liệu giấy đáp ứng nhu cầu nguyên liệu tiêu kế hoạch Tổng Cơng ty giao Các lồi đƣợc gây trồng chủ yếu là: Keo tai tƣợng, Keo lai, Bạch đàn Năm 2010 Công ty đƣợc cấp chứng quản lý rừng bền vững FSC cho 1.918,93 rừng trồng nguyên liệu giấy, sản lƣợng bình quân đạt 70m3/ha Trong năm tới, Công ty đề mục tiêu ổn định nâng cao suất rừng trồng, phấn đấu đạt đƣợc chứng FSC toàn diện tích Trƣớc thực tiễn tác giả tiến hành nghiên cứu đề tài: “ Đánh giá sinh trưởng số lồi trồng ngun liệu giấy Cơng ty Lâm nghiệp Đoan Hùng huyện Đoan Hùng, tỉnh Phú Thọ” nhằm đánh giá sinh trƣởng hiệu kinh tế số lâm phần công ty Phụ biểu 09-c: Kiểm tra giả thuyết luật phân bố N/D1,3 Bạch đàn tuổi D7 Xi- Xi+ ft Xi Xiα Ft.Xiα Pi fl=n.Pi fl gộp fi gộp (fl - fi)^2/fl 1 7.0 0.131 16.066 16.066 7.00 5.116 10 28 252.0 0.298 36.667 36.667 28.00 2.049 12 59 25 1475.0 0.288 35.368 35.368 59.00 15.790 14 19 49 931.0 0.177 21.798 21.798 19.00 0.359 16 10 9 81 729.0 0.076 9.383 12.303 10.00 0.431 18 10 12 11 121 121.0 0.024 2.920 ∑ 123 3515 k=1 λ2 23.745 λ2(0,05) 3.841 Phụ biểu 10-a: Kiểm tra giả thuyết luật phân bố N/Hvn Bạch đàn tuổi Xi Xiα Ft.Xiα fl=n.Pi fl gộp fi gộp H5 Xi- Xi+ ft 14 12 0.5 0.25 0.127 16.086 16.086 12 1.038 15 35 1.5 2.25 78.75 0.292 37.032 37.032 35 0.112 16 42 2.5 6.25 262.5 0.286 36.345 36.345 42 0.88 17 22 3.5 12.25 269.5 0.181 22.991 22.991 22 0.043 18 16 4.5 20.25 324 0.081 10.247 10.247 16 3.231 ∑ 127 Pi 937.75 k=2 (fl - fi)^2/fl λ2 5.302 λ2(0,05) 5.991 Phụ biểu 10-b: Kiểm tra giả thuyết luật phân bố N/Hvn Bạch đàn tuổi H6 Xi- Xi+ ft Xi Xiα Ft.Xiα 14 0.5 0.25 15 20 1.5 2.25 16 43 2.5 17 18 ∑ fi gộp (fl - fi)^2/fl 0.01 12.21 12.21 5.518 45 0.24 29.82 29.82 20 3.235 6.25 268.75 0.27 32.96 32.96 43 3.06 34 3.5 12.25 416.5 0.20 24.91 24.91 34 3.314 22 4.5 20.25 445.5 0.11 14.08 14.08 22.0 4.452 λ2 19.578 λ2(0,05) 5.991 123 Pi fl=n.Pi fl gộp 1176.75 k=2 Phụ biểu 10-c: Kiểm tra giả thuyết luật phân bố N/Hvn Bạch đàn tuổi H7 Xi- Xi+ ft Xi Xiα Ft.Xiα Pi fl=n.Pi fl gộp fi gộp (fl - fi)^2/fl 14 0.5 0.25 0.5 0.08 9.557 9.557 5.975 15 13 1.5 2.25 29.25 0.2 24.441 24.441 13 5.356 16 33 2.5 6.25 206.25 0.24 29.604 29.604 33 0.389 17 43 3.5 12.25 526.75 0.21 25.678 25.678 43 11.685 18 21 4.5 20.25 425.25 0.14 17.436 17.436 21 0.728 19 11 5.5 30.25 332.75 0.08 11 0.206 λ2 23.406 λ2(0,05) 7.815 123 9.595 9.595 1520.75 k=3 Phụ biểu 11: Kết xác định phƣơng trình tƣơng quan H/D Keo lai Model Summary and Parameter Estimates Dependent Variable:H5 Model Summary Parameter Estimates R Equation Square Linear ,529 153,021 136 ,000 6,531 ,483 Logarithmic ,560 173,140 136 ,000 -,539 5,201 ,575 183,979 136 ,000 ,585 95,096 135 ,000 1,365 1,225 ,000 Compound ,549 165,419 136 ,000 7,420 1,043 Power 192,506 136 ,000 3,992 ,454 Inverse Cubic ,586 F df1 df2 Sig Constant 16,816 b1 b2 b3 52,932 -,002 The independent variable is D5 Model Summary and Parameter Estimates Dependent Variable:H6 Model Summary Parameter Estimates R Equation Square F df1 df2 Sig Constant Linear ,469 116,581 132 ,000 Logarithmic ,498 131,018 132 ,000 b1 6,962 ,489 -1,734 5,891 b2 b3 Inverse - ,520 143,026 132 ,000 18,682 ,537 75,997 131 ,000 -5,417 2,610 -,089 ,000 Compound ,487 125,303 132 ,000 7,992 1,040 Power 143,580 132 ,000 3,975 ,472 Cubic ,521 68,505 The independent variable is D6 Model Summary and Parameter Estimates Dependent Variable:H7 Model Summary Parameter Estimates R Equation Square Linear ,394 85,776 132 ,000 Logarithmic ,419 95,248 132 ,000 -1,777 6,023 ,439 103,420 132 ,000 19,426 ,459 55,497 131 ,000 -6,725 2,779 -,092 ,000 Compound ,414 93,215 132 ,000 8,452 1,037 105,091 132 ,000 4,165 ,462 Inverse Cubic Power ,443 F df1 The independent variable is D7 df2 Sig Constant b1 b2 b3 7,435 ,475 74,063 Phụ biểu 12: Kết xác định phƣơng trình tƣơng quan H/D Keo tai tƣợng Model Summary and Parameter Estimates Dependent Variable:H5 Model Summary Parameter Estimates R Equation Square Linear ,417 98,706 138 ,000 Logarithmic ,442 109,382 138 ,000 Inverse ,457 116,243 138 ,000 Cubic ,455 57,078 137 ,000 Compound ,422 100,932 138 ,000 7,280 Power 115,355 138 ,000 4,109 ,430 ,455 F df1 df2 Sig Constant b1 b2 b3 6,558 ,448 ,373 4,673 15,827 -46,036 ,610 1,626 -,056 ,000 1,042 The independent variable is D5 Model Summary and Parameter Estimates Dependent Variable:H6 Model Summary Parameter Estimates R Equation Square F df1 df2 Sig Constant b1 Linear ,374 79,922 134 ,000 7,575 ,415 Logarithmic ,383 83,257 134 ,000 1,038 Inverse ,385 83,716 134 ,000 Cubic ,386 41,757 133 ,000 b2 b3 ,000 ,000 4,663 16,838 -50,018 4,918 ,777 Compound ,381 82,649 134 ,000 8,329 Power 87,327 134 ,000 4,880 ,380 ,395 1,034 The independent variable is D6 Model Summary and Parameter Estimates Dependent Variable:H7 Model Summary Parameter Estimates R Equation Square Linear ,387 84,514 134 ,000 Logarithmic ,394 87,024 134 ,000 Inverse ,392 86,273 134 ,000 Cubic ,397 43,703 133 ,000 5,154 ,766 Compound ,393 86,793 134 ,000 8,631 Power 90,516 134 ,000 4,853 ,394 ,403 F df1 The independent variable is D7 df2 Sig Constant b1 b2 b3 ,000 ,000 7,764 ,431 ,344 5,096 17,878 -57,667 1,034 Phụ biểu 13: Kết xác định phƣơng trình tƣơng quan H/D Bạch đàn Model Summary and Parameter Estimates Dependent Variable:H5 Model Summary Parameter Estimates R Equation Linear Square Constan F df1 df2 Sig ,846 686,201 125 ,000 ,857 748,320 125 ,000 Inverse ,837 643,606 125 ,000 Cubic ,861 382,983 124 ,000 ,828 603,494 125 ,000 ,848 697,494 125 ,000 Logarithm ic Compoun d Power t b1 b2 b3 10,089 ,555 2,894 5,582 21,075 -52,364 8,119 ,850 10,929 ,000 ,000 1,036 6,864 ,359 The independent variable is D5 Model Summary and Parameter Estimates Dependent Variable:H6 Model Summary Parameter Estimates R Equation Linear Logarithm ic Inverse Square Constan F df1 df2 Sig ,819 548,481 121 ,000 ,833 605,672 121 ,000 ,821 556,323 121 ,000 t b1 10,599 ,516 2,924 5,607 21,649 -57,318 b2 b3 Cubic Compoun d Power ,838 309,841 120 ,000 ,803 493,192 121 ,000 ,824 566,422 121 ,000 8,172 ,854 11,414 ,000 ,000 1,033 7,061 ,350 The independent variable is D6 Model Summary and Parameter Estimates Dependent Variable:H7 Model Summary Parameter Estimates R Equation Linear Square Constan F df1 df2 Sig ,772 409,130 121 ,000 ,788 449,310 121 ,000 Inverse ,781 431,393 121 ,000 Cubic ,791 227,145 120 ,000 ,755 371,902 121 ,000 ,777 421,034 121 ,000 Logarithm ic Compoun d Power The independent variable is D7 t b1 b2 b3 10,899 ,497 3,026 5,599 21,944 -59,584 8,292 ,847 11,690 1,031 7,198 ,344 ,000 ,000 Phụ biểu 14: Chi phí đầu tƣ 1ha rừng Keo STT Hạng mục I Chi phí trực tiếp Trồng rừng '' Trồng rừng '' Cây giống '' Phân bón II III ĐVT đ/ha '' Chăm sóc '' Chăm sóc năm '' Chăm sóc năm '' Chăm sóc năm '' Bảo vệ rừng '' Bảo vệ rừng năm '' Bảo vệ rừng năm '' Bảo vệ rừng năm '' Bảo vệ rừng năm Chi phí gián tiếp BHXH,YT