(Luận văn) đá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 lu an Hà Nội, ngày 30 tháng 11 năm 2017 n va Người cam đoan gh tn to p ie Nguyễn Đức Nghĩa d oa nl w ll u nf va an lu oi m z at nh z m co l gm @ an Lu n va ac th si 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 lu an nghiệp; anh, chị, em, bạn bè đồng nghiệp tạo điều kiện thuận lợi cho n va 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ỏ tn to 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 gh p ie 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 nl w Mặc dù cố gắng nỗ lực, nhƣng kinh nghiệm nghiên cứu chƣa d oa 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 an lu không tránh khỏi thiếu sót định Tác giả mong nhận đƣợc góp ý va thầy giáo bạn bè đồng nghiệp luận văn đƣợc hoàn chỉnh u nf Xin chân thành cảm ơn! ll Hà Nội, ngày 30 tháng 11 năm 2017 oi m Tác giả z at nh z m co l gm @ Nguyễn Đức Nghĩa an Lu n va ac th si 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 lu DANH MỤC CÁC HÌNH viii an ĐẶT VẤN ĐỀ va n Chƣơng TỔNG QUAN VẤN ĐỀ NGHIÊN CỨU gh tn to 1.1 Nghiên cứu giới p ie 1.1.1 Nghiên cứu Keo tai tƣợng 1.1.2 Nghiên cứu Keo lai nl w 1.1.3 Các nghiên cứu Bạch đàn d oa 1.2 Nghiên cứu Việt Nam an lu 1.2.1 Những nghiên cứu trồng rừng nguyên liệu u nf va 1.2.2 Những nghiên cứu Keo tai tƣợng 1.2.3 Nghiên cứu Keo lai ll oi m 1.2.3 Nghiên cứu Bạch đàn 11 z at nh 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 z @ 2.1 Mục tiêu nghiên cứu 15 l gm 2.2 Đối tƣợng phạm vi nghiên cứu 15 m co 2.3 Nội dung nghiên cứu 15 2.4 Phƣơng pháp nghiên cứu 15 an Lu 2.4.1 Phƣơng pháp kế thừa 15 n va ac th si 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 lu 3.1.4 Đặc điểm đất 25 an 3.2 Điều kiện kinh tế-xã hội 25 va n 3.2.1 Sản xuất lâm nghiệp 25 gh tn to 3.2.2 Dịch vụ 25 ie 3.2.3 Lao động 25 p 3.2.4 Cơ sở hạ tầng 26 nl w 3.3 Lƣợc sử trồng rừng 26 d oa 3.3.1 Loài Keo tai tƣợng , keo lai 26 an lu 3.3.2 Bạch đàn 29 u nf va Chƣơng KẾT QUẢ NGHIÊN CỨU 32 4.1 Tình hình sinh trƣởng lâm phần 32 ll oi m 4.1.1 Kiểm tra OTC 32 z at nh 4.1.2 Sinh trƣởng đƣờng kính D1,3 33 4.1.3 Sinh trƣởng chiều cao vút 37 z @ 4.2 Một số quy luật kết cấu lâm phần 40 l gm 4.2.1 Quy luật phân bố số theo cỡ đƣờng kính 40 m co 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 an Lu 4.3 Trữ lƣợng tăng trƣởng lâm phần 47 n va ac th si 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 lu 5.2 Tồn 54 an 5.3 Khuyến nghị 55 va n TÀI LIỆU THAM KHẢO p ie gh tn to d oa nl w ll u nf va an lu oi m z at nh z m co l gm @ an Lu n va ac th si vi DANH MỤC CÁC TỪ VIẾT TẮT Nguyên nghĩa Ký hiệu Bộ NN&PTNT Bộ Nông nghiệp Phát triển nơng thơn Đƣờng kính thân vị trí 1,3 m D1,3 Dt Đƣờng kính chiếu tán Fl Phân bố lý thuyết fi Phân bố thực nghiệm lu Hvn Chiều cao vút Hdc Chiều cao dƣới cành an n va OTC Ô tiêu chuẩn Quản lý sử dụng SXKD Sản xuất kinh doanh p ie gh tn to QLSD d oa nl w ll u nf va an lu oi m z at nh z m co l gm @ an Lu n va ac th si vii DANH MỤC CÁC BẢNG STT Tên bảng Trang lu an n va 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 gh tn to 4.1 4.11 Chi phí đầu tƣ cho rừng trồng p ie 49 4.12 Chi phí khai thác vận chuyển (đ/m3) w 50 50 oa nl 4.13 Thu nhập cho rừng trồng 4.14 Tổng hợp tiêu đánh giá hiệu kinh tế d 51 ll u nf va an lu oi m z at nh z m co l gm @ an Lu n va ac th si 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 35 Bạch đàn tuổi 5; 6; Biểu đồ sinh trƣởng Hvn Keo tai tƣợng, Keo lai 38 lu Bạch đàn tuổi 5; 6; an n va 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 Biểu đồ phân bố N/Hvn Keo tai tƣợng 44 gh tn to 4.6 p ie 4.9 4.10 Biểu đồ phân bố N/Hvn Keo lai w 44 4.11 Biểu đồ phân bố N/Hvn Bạch đàn d oa nl 45 ll u nf va an lu oi m z at nh z m co l gm @ an Lu n va ac th si ĐẶ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 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 lu mà mục tiêu, phƣơng hƣớng giải nƣớc ta an Tập đoàn trồng rừng sản xuất gồm loài nhƣ Keo tai tƣợng, va n Keo lai, Bạch đàn, Bồ đề, Thông, Mỡ, Lát hoa, Trám, Xoan, Giổi, Lim xanh, gh tn to Dẻ, Chò chỉ, Sao đen, Trong đó, Keo Bạch đàn hai lồi chủ lực với tỷ ie lệ gây trồng 80% diện tích Đây lồi mọc nhanh, sinh trƣởng tốt, p phù hợp với nhiều lập địa khác nhau, chu kỳ kinh doanh ngắn – năm nl w – 15 năm tùy mục đích kinh doanh Từ đó, nhiều cơng ty, địa phƣơng, đơn d oa vị chủ rừng lựa chọn trồng để phát triển trồng rừng kinh tế Tuy an lu nhiên, trồng rừng kinh tế chƣa đƣợc đồng hóa dẫn đến suất u nf va 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 ll oi m Công ty Lâm nghiệp Đoan Hùng nằm địa bàn huyện Đoan Hùng z at nh tỉnh Phú Thọ đơn vị hạch tố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 z trồng rừng nguyên liệu giấy đáp ứng nhu cầu nguyên liệu tiêu kế hoạch @ l gm Tổng Cơng ty giao Các lồi đƣợc gây trồng chủ yếu là: Keo tai m co 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 an Lu bình quân đạt 70m3/ha Trong năm tới, Công ty đề mục tiêu ổn n va ac th si định nâng cao suất rừng trồng, phấn đấu đạt đƣợc chứng FSC tồ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ố 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ọ” nhằm đánh giá sinh trƣởng hiệu kinh tế số lâm phần công ty lu an n va p ie gh tn to d oa nl w ll u nf va an lu oi m z at nh z m co l gm @ an Lu n va ac th si 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 lu ∑ 123 3515 an 23.745 λ2(0,05) 3.841 n va k=1 λ2 gh tn to Phụ biểu 10-a: Kiểm tra giả thuyết luật phân bố N/Hvn Bạch đàn tuổi H5 Xi- Xi+ Xi Xiα Ft.Xiα fl=n.Pi fl gộp Pi fi gộp (fl - fi)^2/fl p ie ft 12 0.5 0.25 0.127 16.086 16.086 12 1.038 35 1.5 2.25 78.75 0.292 37.032 37.032 35 0.112 oa 6.25 262.5 0.286 36.345 36.345 42 0.88 22 3.5 12.25 269.5 0.181 22.991 22.991 22 0.043 va 15 16 42 2.5 17 18 16 4.5 20.25 324 0.081 10.247 10.247 16 3.231 nl w d an lu 937.75 ll 127 u nf ∑ 14 m oi k=2 λ2 5.302 λ2(0,05) 5.991 z at nh z m co l gm @ an Lu n va ac th si 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 fl=n.Pi fl gộp Pi 1176.75 lu k=2 an n va tn to 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+ Xi Xiα Ft.Xiα Pi fl=n.Pi fl gộp fi gộp (fl - fi)^2/fl 0.5 0.25 0.5 0.08 9.557 9.557 5.975 1.5 2.25 29.25 0.2 24.441 24.441 13 5.356 6.25 206.25 0.24 29.604 29.604 33 0.389 526.75 0.21 25.678 25.678 43 11.685 p ie gh 14 ft 13 16 33 nl 2.5 17 43 3.5 12.25 d 18 21 4.5 20.25 an 425.25 0.14 17.436 17.436 21 0.728 19 11 5.5 30.25 va 332.75 0.08 11 0.206 λ2 23.406 λ2(0,05) 7.815 w 15 oa lu u nf 9.595 1520.75 ll 123 9.595 m oi k=3 z at nh z m co l gm @ an Lu n va ac th si 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 lu an va 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 192,506 136 ,000 3,992 ,454 n Inverse F df1 df2 Sig Constant gh tn to 16,816 b2 b3 52,932 -,002 p ie Cubic b1 ,586 oa nl w Power d The independent variable is D5 an lu u nf va Model Summary and Parameter Estimates Dependent ll oi m Variable:H6 Parameter Estimates z at nh Model Summary R z Square F df1 df2 Sig Constant Linear ,469 116,581 132 ,000 Logarithmic ,498 131,018 132 ,000 b1 b2 b3 6,962 ,489 m co l gm @ Equation -1,734 5,891 an Lu n va ac th si 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 lu Dependent an Variable:H7 va n Model Summary Parameter Estimates tn to R gh Square F df1 df2 Sig Constant b1 b2 b3 p ie Equation Linear ,394 132 ,000 7,435 ,475 95,248 132 ,000 -1,777 6,023 103,420 132 ,000 19,426 131 ,000 -6,725 2,779 -,092 ,000 132 ,000 8,452 1,037 132 ,000 4,165 ,462 nl w 85,776 Cubic u nf va an ,439 lu Inverse d oa Logarithmic ,419 55,497 Compound ,414 93,215 74,063 ll ,459 - oi m ,443 105,091 The independent variable is D7 z at nh Power z m co l gm @ an Lu n va ac th si 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 lu an n va Square Linear ,417 98,706 138 ,000 Logarithmic ,442 109,382 138 ,000 Inverse ,457 116,243 138 ,000 ,455 57,078 137 ,000 Compound ,422 100,932 138 ,000 7,280 115,355 138 ,000 4,109 ,430 tn to Equation p ie gh Cubic ,455 df1 df2 Sig Constant b1 b2 b3 6,558 ,448 ,373 4,673 15,827 -46,036 ,610 1,626 -,056 ,000 1,042 nl w Power F d oa The independent variable is D5 an lu Model Summary and Parameter Estimates va Dependent ll u nf Variable:H6 Parameter Estimates oi m Model Summary F z at nh R df1 df2 Sig Constant ,374 79,922 134 ,000 Logarithmic ,383 83,257 134 ,000 Inverse ,385 83,716 134 ,000 Cubic ,386 41,757 133 ,000 gm @ b2 b3 ,000 ,000 7,575 ,415 1,038 4,663 16,838 -50,018 m co Linear b1 l Square z Equation 4,918 ,777 an Lu n va ac th si 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 lu an Square F df1 df2 Sig Constant Linear ,387 84,514 134 ,000 Logarithmic ,394 87,024 134 ,000 b1 b2 b3 ,000 ,000 7,764 ,431 n va Equation 5,096 Inverse ,392 86,273 134 ,000 Cubic ie ,397 43,703 133 ,000 5,154 ,766 Compound ,393 86,793 134 ,000 8,631 90,516 134 ,000 4,853 ,394 p gh tn to ,344 nl w ,403 1,034 d oa Power 17,878 -57,667 ll u nf va an lu The independent variable is D7 oi m z at nh z m co l gm @ an Lu n va ac th si 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 Constan Square df1 df2 Sig 686,201 125 ,000 ,857 748,320 125 ,000 Inverse ,837 643,606 125 ,000 tn lu ,846 to Linear F Cubic ,861 382,983 124 ,000 603,494 125 ,000 697,494 125 ,000 an Logarithm n va ic t b1 b2 b3 10,089 ,555 2,894 5,582 21,075 -52,364 ,000 ,000 ie gh 8,119 ,850 p Compoun ,828 10,929 1,036 ,848 oa Power nl w d 6,864 ,359 d The independent variable is D5 an lu Model Summary and Parameter Estimates va u nf Dependent ll Variable:H6 oi m R Sig t 548,481 121 ,000 ,833 605,672 121 ,000 ,821 556,323 121 ,000 b2 b3 10,599 ,516 2,924 5,607 an Lu ,819 b1 m co Inverse df2 l ic df1 gm Logarithm F @ Linear Square Constan z Equation Parameter Estimates z at nh Model Summary 21,649 -57,318 n va ac th si 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 lu an Model Summary Parameter Estimates n va R Square ,772 df1 df2 Sig 409,130 121 ,000 449,310 121 ,000 431,393 121 ,000 120 ,000 121 ,000 121 ,000 t b1 b2 b3 10,899 ,497 ie gh Linear F ll tn to Equation Constan p Logarithm ,788 3,026 5,599 ,791 371,902 ,777 421,034 21,944 -59,584 8,292 ,847 11,690 ,000 ,000 1,031 u nf 7,198 ,344 oi z at nh The independent variable is D7 m Power ,755 va d an Compoun 227,145 lu Cubic d ,781 oa Inverse nl w ic z m co l gm @ an Lu n va ac th si Phụ biểu 14: Chi phí đầu tƣ 1ha rừng Keo STT Hạng mục I Chi phí trực tiếp ĐVT đ/ha Trồng rừng '' Trồng rừng '' Cây giống '' '' Phân bón Chăm sóc năm '' Chăm sóc năm '' Chăm sóc năm '' '' n va '' an Chăm sóc ie lu gh tn to p Khối lƣợng Tính cho 20,764,118.30 10,787,744.29 107,423 66.98 7,195,192.54 400 1834.00 733,600.00 4,290 666.50 2,858,951.75 8,257,606.01 107,423 45.10 4,844,777.30 107,423 27.27 2,929,425.21 107,423 4.50 483,403.50 1,718,768.00 oa nl w Bảo vệ rừng Đơn giá '' d Bảo vệ rừng năm lu an '' va Bảo vệ rừng năm u nf '' Bảo vệ rừng năm 4.00 429,692.00 107,423 4.00 429,692.00 107,423 4.00 429,692.00 107,423 4.00 429,692.00 9,460,391.63 ll 107,423 '' oi đ/ha z at nh '' 2,589,058.63 z '' đ/ha 871,333.00 30,224,509.93 m co Thiết kế lâm sinh Tổng chi phí 6,000,000.00 l '' gm Quản lý công ty @ III m II Bảo vệ rừng năm Chi phí gián tiếp BHXH,YT