Lập biểu thể tích, cấp đất và sản lượng cho thông ba lá (pinus kesiya royle ex gordon) ở hà giang

BỘ GIÁO DỤC VÀ ĐÀO TẠO BỘ NÔNG NGHIỆP VÀ PTNT TRƯỜNG ĐẠI HỌC LÂM NGHIỆP PHẠM QUANG TUYẾN LẬP BIỂU THỂ TÍCH, CẤP ĐẤT VÀ SẢN LƯỢNG CHO THÔNG BA LÁ (PINUS KESIYA ROYLE EX.GORDON) Ở HÀ GIANG CHUYÊN NGÀNH LÂM HỌC MÃ SỐ: 60.62.60 LUẬN VĂN THẠC SỸ KHOA HỌC LÂM NGHIỆP HƯỚNG DẪN KHOA HỌC: PGS.TS TRẦN VĂN CON TS PHAN MINH SÁNG Hà Nội, 2010 BỘ GIÁO DỤC VÀ ĐÀO TẠO BỘ NÔNG NGHIỆP VÀ PTNT TRƯỜNG ĐẠI HỌC LÂM NGHIỆP PHẠM QUANG TUYẾN LẬP BIỂU THỂ TÍCH, CẤP ĐẤT VÀ SẢN LƯỢNG CHO THÔNG BA LÁ (PINUS KESIYA ROYLE EX.GORDON) Ở HÀ GIANG CHUYÊN NGÀNH LÂM HỌC MÃ SỐ: 60.62.60 LUẬN VĂN THẠC SỸ KHOA HỌC LÂM NGHIỆP Hà Nội, 2010 i LỜI CẢM ƠN Để hoàn thành chương trình đào tạo Cao học Lâm nghiệp khoá học 20082010, đồng ý Khoa sau đại học - Trường Đại học Lâm nghiệp, thực đề tài tốt nghiệp: Tên đề tài: “Lập biểu thể tích, cấp đất sản lượng cho Thông ba (Pinus kesiya Royle ex.Gordon) Hà Giang” Sau thời gian tiến hành làm đề tài tốt nghiệp đến luận văn hoàn thành Cho phép bày tỏ lòng biết ơn chân thành đến thầy giáo, đặc biệt PGS.TS Trần Văn Con, TS Phan Minh Sáng tận tình hướng dẫn, bảo giúp đỡ hoàn thành luận văn Tôi xin bày tỏ lòng biết ơn chân thành đến TS Đào Công Khanh, GS.TS Vũ Tiến Hinh, TS Phạm Ngọc Giao Thầy giáo môn Điều tra - Quy hoạch rừng cho Tôi ý kiến đóng góp quý báu Cuối xin bày tỏ lòng biết ơn đến tất bạn bè đồng nghiệp người thân giúp đỡ có luận văn Tác giả vui lòng nhận góp ý, bổ sung bạn đọc để luận văn hoàn chỉnh Tôi xin chân thành cảm ơn ! Hà Nội, tháng 09 năm 2010 Tác giả Phạm Quang Tuyến ĐẶT VẤN ĐỀ Bảng biểu điều tra kinh doanh rừng có vai trò vị trí đặc biệt sản xuất lâm nghiệp không khả ứng dụng: tính thuận tiện, tính khoa học, tính xác mà khả dự báo sản lượng, suất rừng tương lai Với vai trò to lớn đó, bảng biểu điều tra nhiều nước tiên tiến châu Âu, châu Mỹ vận dụng vào việc kinh doanh rừng từ kỷ XIX Trong vừa năm qua, thực chương trình trồng triệu rừng, diện tích rừng trồng nước ta ngày tăng, cấu loài trồng ngày nhiều như: Thông bá lá, Thông mã vĩ, Thông nhựa, Keo tai tượng, Keo tràm, Keo lai, Mỡ, Bồ đề, Sa Mộc, Tếch, Quế, Bạch đàn Uro, hàng loạt bảng biểu điều tra xây dựng nhằm đáp ứng thực tiễn sản xuất Tuy nhiên, điều kiện thời gian, kinh phí hạn hẹp nên số bảng biểu xây dựng áp dụng cho vùng, địa phương chương trình mà chưa xây dựng cho phạm vi toàn quốc Trong nghiên cứu xây dựng bảng biểu điều tra có điểm quan trọng loài vùng sinh thái lại có kiểu hình dạng riêng, kiểu sinh trưởng riêng, cấp suất khác nhau, Do đó, nghiên cứu lập biểu điều tra, việc thu thập số liệu phải đại diện cho vùng sinh thái, kiểu sinh trưởng đặc biệt kết nghiên cứu phải kiểm tra với số liệu độc lập không tham gia vào trình lập biểu Bảng biểu lập cho Thông bá Tây Nguyên hệ thống bảng biểu xây dựng đầy đủ cho rừng trồng loại áp dụng kinh doanh rừng Thông ba Tây Nguyên Nhưng ví dụ điển hình cho việc xây dựng bảng biểu dừng lại vùng sinh thái mà chưa có kết kiểm tra ứng dụng cho vùng sinh thái đặc trưng khác Hà Giang nằm vùng nhiệt đới gió mùa miền núi cao, khí hậu Hà Giang mang đặc điểm vùng núi Việt Bắc – Hoàng Liên Sơn mát, lạnh nhiều sương mù, nhiệt độ trung bình năm khoảng 21,60C - 23,9 C Điều kiện kiện tự nhiên Hà Giang thích hợp cho Thông ba sinh trưởng, phát triển mở rộng kinh doanh rừng Thông ba Theo thống kê sơ huyện Hoàng Su Phì tỉnh Hà Giang có khoảng 7000 rừng Thông ba (Báo cáo Ban quản lý rừng phòng hộ Hoàng Su Phì) Diện tích trồng với mục đích trước rừng phòng hộ nên không khai thác Nhưng số năm gần đây, theo chủ trương nhà nước chuyển đổi số rừng phòng hộ thành rừng sản xuất việc kinh doanh rừng Thông ba Hà Giang gặp phải khó khăn nghiên cứu quy trình hướng dẫn cụ thể Từ đòi hỏi thực tế sản xuất, cần phải có biểu sản lượng để phục vụ việc kinh doanh rừng Thông ba địa phương Nhưng có khó khăn việc sử dụng bảng biểu lập áp dụng cho Hà Giang điều kiện sinh thái hai vùng hoàn toàn khác xa Điều làm cho điều kiện sinh trưởng Thông ba hai vùng khác dẫn đến có sai số lớn điều tra Đây câu hỏi đặt thực tiễn cần phải giải Trước yêu cầu sản xuất việc kinh doanh rừng Thông ba Hà Giang vấn đề đặt là: Các biểu điều tra rừng Thông ba lập cho vùng Tây Nguyên sử dụng cho rừng Thông ba Hà Giang không? Nếu sử dụng độ xác bao nhiêu? Hay phải phải xây dựng số bảng biểu phù hợp với điều kiện sinh thái địa phương? Vì vậy, để sử dụng biểu Thông ba lập cho vùng Tây Nguyên ứng dụng hiệu cho tỉnh Hà Giang địa phương miền núi phía Bắc, cần phải có số liệu kiểm tra, điều chỉnh biểu xây dựng sở kế thừa kết nghiên cứu biểu điều tra rừng Thông ba Tây Nguyên Từ yêu cầu thực tiễn trên, tiến hành thực đề tài: “Lập biểu thể tích, cấp đất sản lượng cho Thông ba (Pinus kesiya Royle ex.Gordon) Hà Giang” Chương TỔNG QUAN VỀ VẤN ĐỀ NGHIÊN CỨU 1.1 Ngoài nước 1.1.1 Nghiên cứu lập bảng biểu điều tra rừng giới Nghiên cứu sinh trưởng, suất rừng nội dung môn điều tra sản lượng rừng (Growth an yield study) Bộ môn khoa học đời từ kỷ XVIII, gắn liền với tên tuổi Oettelt (1765) Paulsen (1787) Vào thời điểm nhu cầu buôn bán gỗ, rừng dẫn tới đời môn khoa học Trong giai đoạn phát triển đầu tiên, sản lượng rừng chưa gắn với quan điểm sinh thái học Việc nghiên cứu hoàn thiện phát triển thêm bước (Trịnh Đức Huy, 1988) Bảng biểu điều tra rừng ngày có bước tiến quan trọng việc ứng dụng vào việc kinh doanh rừng Việc ứng dụng bảng biểu: thể tích, thương phẩm, biểu cấp đất, vào kinh doanh rừng công cụ đắc lực khó thay hoạt động sản xuất lâm nghiệp tính tiện lợi khoa học tính toán Ngày nay, với phát triển công nghệ thông tin, công nghệ phần mềm, tính toán bảng biểu điều tra lập trình nhằm tính toán nhanh hơn, sử dụng tiện lợi hơn; thông số, tiêu cần tính toán xây dựng nhằm phát huy tối đa vai trò sản xuất lâm nghiệp (Vũ Tiến Hinh., 2003) 1.1.1.1 Biểu thể tích Nghiên cứu hình dạng thân biểu thể tích tác giả số nước Châu Âu nghiên cứu từ kỷ XIX (Van Laar & Akca, 2007) Tuy nhiên, bảng biểu giai đoạn dừng lại việc nghiên cứu cá lẻ thực nghiệm tính toán trực tiếp Các nghiên cứu hoàn thiện đến đầu kỷ XX, nhiều tác giả vận dụng hàm toán học vào việc tính toán biểu thể tích dựa vào tương quan đường kính thể tích (Van Laar & Akca, 2007) Nghiên cứu bảng biểu nhà nghiên cứu hệ thống hóa sâu vào số loài cụ thể nhiều nước Châu Âu, Châu Mỹ, Châu Phi, Châu Á (Van Laar Akca, 2007; Vũ Tiến Hinh Phạm Ngọc Giao, 1997) Ngày nay, với phát triển công nghệ máy tính, nhà nghiên cứu sản lượng xây dựng phần mềm máy tính thành chương trình để lập biểu thể tích cho nhiều loài (Avery & Burkhart, 1994; Husch et al., 2003; Van Laar & Akca, 2007) Đây tiến vượt bậc việc nghiên cứu lập biểu thể tích, việc sử dụng phần mềm máy tính giúp cho việc tính toán xác hơn, tiện lợi so với loại bảng biểu cũ 1.1.1.2 Nghiên cứu sinh trưởng Các nhà nghiên cứu rằng, đời sống rừng, lâm phần rừng chia thành giai đoạn: non, sào, trung niên, thành thục, thành thục (Belov, 1983) giai đoạn khác rừng có sức sinh trưởng phát triển khác Burkhat (1982) tổng kết có ba hướng nghiên cứu sinh trưởng cho rừng đồng tuổi (Trịnh Đức Huy., 1988) sau: * Tiếp cận toàn lâm phần (Whole stand): Theo hướng này, ta nghiên cứu mối quan hệ tiêu sinh trưởng lâm phần trữ lượng/ha, sản lượng/ha… với nhân tố tạo thành tuổi lâm phần, chiều cao, đường kính bình quân lâm phần, tổng diện ngang/ha, mật độ, số cấp đất …rất quan tâm Đầu kỷ XX, có biểu sản lượng xây dựng phương pháp đồ thị (Missellaneous Publication 50 - dẫn theo Burkhart,1982) Phương trình tương quan nhiều biến số Markiney Chaiken (1959) sử dụng để lập biểu sản lượng cho rừng thông Lobelly (dẫn theo Burkhart, 1982) * Tiếp cận theo phân bố kích cỡ ((Size class distribution): Việc nghiên cứu phân phối trữ lượng rừng theo kích cỡ, thường áp dụng cho cỡ đường kính Hướng nghiên cứu phát triển mạnh từ năm 60 trở lại với nhiều tác giả tiêu biểu như: Cluter (1967, 1971), Lekhart (1971, 1972), Burkhart (1974)… Trên cỡ kính lấy giá trị chiều cao, thể tích, số cây… quy luật tương quan hay thống kê trị số trung bình sử dụng mẫu biểu lập sẵn Vấn đề cốt lõi phương pháp nghiên cứu xác định hàm phân phối mật độ cho cỡ đường kính trình vận động theo thời gian rừng * Tiếp cận theo cá lẻ (Invidual tree): Một số tác giả tập trung vào nghiên cứu quy luật sinh trưởng đại diện lâm phần như: trung bình, sinh trưởng tốt nhất, từ đó, ta thấy có hai cách nghiên cứu là: xem sinh trưởng độc lập với (Distance independent) có xem xét với bạn lâm phần (Distance dependent) Trong trường hợp thứ hai, gắn số cạnh tranh (Conpettition index) coi số yếu tố ảnh hưởng đến sinh trưởng (Trịnh Đức Huy., 1988) Theo Pretzsch (2001, 2009), mối quan hệ để xây dựng nên biểu sản lượng, là: 1, Chiều cao lâm phần = f (tuổi lâm phần) 2, Tổng trữ lượng lâm phần = f (chiều cao lâm phần) 3, Tổng trữ lượng lâm phần = f (tuổi lâm phần) Các hàm sinh trưởng thông dụng phù hợp cho mô sinh trưởng rừng, rừng cho (Zeide, 1993): Hossfeld IV: y  t c /(b  t c / a) Gompertz: y  aebe (1.1)  bt (1.2) Logistic: y  a /(1  cebt ) (1.3) Schumacher: y  aeb / t (1.4) Korf: y  aebt c c (1.5) Trong y = biến lâm phần (chiều cao, đường kính thể tích); t = tuổi; y’ = tăng trưởng; a, b, c d = tham số phương trình; ln = hàm logarit tự nhiên Phương pháp xây dựng mô hình sinh trưởng, sản lượng giới 71 Bảng 4-19: Biểu trình sinh trưởng cấp đất II Bộ phận nuôi dưỡng Bộ phận tổng hợp Tia thua A N Hg Dg G M PV% ∆M ZM N Hgc Dgc Gc Mc G M ZM ∆M PM % 1242 6,92 10,20 10,13 26,73 26,21 4,45 7,01 99 5,47 7,45 0,43 0,97 10,56 27,70 7,14 4,62 25,79 1160 8,09 11,20 11,42 34,51 22,56 4,93 7,79 82 6,39 8,18 0,43 1,11 11,85 35,63 7,93 5,09 22,25 1090 9,21 12,20 12,74 43,07 19,86 5,38 8,55 70 7,28 8,92 0,44 1,26 13,18 44,33 8,70 5,54 19,63 1029 10,29 13,21 14,09 52,40 17,80 5,82 9,33 61 8,13 9,65 0,44 1,42 14,54 53,82 9,48 5,98 17,62 10 975 11,33 14,23 15,49 62,51 16,18 6,25 10,12 54 8,95 10,40 0,46 1,58 15,95 64,10 10,28 6,41 16,04 11 927 12,33 15,26 16,95 73,45 14,89 6,68 10,93 48 9,74 11,15 0,47 1,76 17,42 75,21 11,11 6,84 14,77 12 883 13,29 16,32 18,46 85,23 13,82 7,10 11,78 44 10,50 11,93 0,49 1,95 18,95 87,17 11,97 7,26 13,73 13 842 14,22 17,41 20,04 97,90 12,94 7,53 12,67 40 11,24 12,72 0,51 2,15 20,56 100,04 12,87 7,70 12,86 14 805 15,13 18,53 21,70 111,51 12,20 7,96 13,61 37 11,95 13,54 0,54 2,36 22,24 113,87 13,82 8,13 12,14 15 771 16,00 19,68 23,44 126,11 11,58 8,41 14,60 35 12,64 14,38 0,56 2,59 24,00 128,70 14,83 8,58 11,53 16 735 16,92 20,99 25,45 143,27 11,98 8,95 17,16 35 13,37 15,34 0,65 3,14 26,10 146,41 17,71 9,15 12,09 17 703 17,82 22,36 27,58 161,82 11,46 9,52 18,55 33 14,08 16,34 0,69 3,47 28,27 165,29 18,88 9,72 11,42 18 667 18,82 24,02 30,21 185,05 12,55 10,28 23,23 35 14,87 17,55 0,86 4,50 31,06 189,55 24,26 10,53 12,80 19 633 19,82 25,83 33,14 211,38 12,46 11,13 26,34 34 15,66 18,87 0,96 5,24 34,10 216,63 27,08 11,40 12,50 20 600 20,82 27,82 36,44 241,36 12,42 12,07 29,98 33 16,45 20,33 1,08 6,13 37,52 247,49 30,87 12,37 12,47 72 Bảng 4-20: Biểu trình sinh trưởng cấp đất III Bộ phận nuôi dưỡng A Tỉa thưa Bộ phận tổng hợp N Hg Dg G M PV% ∆M ZM N Hgc Dgc Gc Mc G M ZM ∆M PM% 1.403 5,03 8,62 8,2 16,3 27,14 2,71 4,42 104 3,97 6,30 0,33 0,56 8,5 16,8 4,5 2,8 26,7 1.318 5,96 9,39 9,1 21,1 22,82 3,01 4,81 85 4,71 6,86 0,31 0,62 9,4 21,7 4,9 3,1 22,5 1.247 6,85 10,13 10,1 26,2 19,65 3,28 5,16 71 5,41 7,40 0,30 0,68 10,4 26,9 5,2 3,4 19,4 1.187 7,69 10,85 11,0 31,7 17,23 3,52 5,46 61 6,07 7,93 0,30 0,74 11,3 32,5 5,5 3,6 17,0 10 1.134 8,49 11,55 11,9 37,4 15,33 3,74 5,74 53 6,71 8,44 0,30 0,80 12,2 38,2 5,8 3,8 15,2 11 1.087 9,26 12,24 12,8 43,4 13,80 3,95 6,00 47 7,31 8,95 0,29 0,85 13,1 44,3 6,1 4,0 13,7 12 1.046 9,99 12,92 13,7 49,7 12,55 4,14 6,24 42 7,89 9,44 0,29 0,91 14,0 50,6 6,3 4,2 12,4 13 1.008 10,69 13,59 14,6 56,1 11,51 4,32 6,46 38 8,44 9,93 0,29 0,96 14,9 57,1 6,5 4,4 11,4 14 974 11,36 14,26 15,5 62,8 10,64 4,49 6,68 34 8,97 10,42 0,29 1,01 15,8 63,8 6,7 4,6 10,5 15 942 12,00 14,92 16,5 69,7 9,89 4,65 6,89 31 9,48 10,90 0,29 1,06 16,8 70,8 6,9 4,7 9,8 16 913 12,61 15,57 17,4 76,8 9,24 4,80 7,10 29 9,96 11,38 0,29 1,11 17,7 77,9 7,2 4,9 9,2 17 887 13,20 16,22 18,3 84,1 8,68 4,95 7,30 27 10,43 11,86 0,30 1,17 18,6 85,3 7,4 5,0 8,6 18 862 13,77 16,88 19,3 91,6 8,19 5,09 7,50 25 10,88 12,33 0,30 1,22 19,6 92,8 7,6 5,2 8,1 19 838 14,32 17,53 20,2 99,3 7,75 5,23 7,70 23 11,31 12,81 0,30 1,27 20,5 100,6 7,8 5,3 7,7 20 816 14,85 18,18 21,18 107,2 7,36 5,36 7,90 22 11,733 13,28 0,30 1,32 21,5 108,5 7,95 5,427 7,3 73 Chương KẾT LUẬN – TỒN TẠI – KHUYẾN NGHỊ 5.1 Kết luận Từ kết nghiên cứu, phân tích tiêu thống kê mô tả qui luật sinh trưởng Thông ba Hà Giang, rút kết luận sau: 5.1.1 Kết nghiên cứu ứng dụng bảng biểu Thông ba Hà Giang Do đặc điểm sinh thái vùng Tây Nguyên Hà Giang khác nên loài Thông ba sinh trưởng hai vùng có sai khác đáng kể Vì biểu thể tích biểu cấp đất lập cho Thông ba Tây Nguyên chưa thể áp dụng cho Thông ba Hà Giang - Sai số lớn biểu thể tích hai nhân tố kiểm tra biểu 46,98%, sai số bình quân 15,60% Do đó, cần lập biểu thể tích hai nhân tố cho Thông ba Hà Giang - Kiểm tra biểu cấp đất cho thấy đường sinh trưởng chiều cao tầng trội H0 biểu cấp đất lập cho Thông ba Tây Nguyên chưa phản ánh hết quy luật sinh học Thông ba Hà Giang Từ suy biểu sản lượng lập không phù hợp biểu sản lượng xây dựng dựa cấp đất lập cho biểu cấp đất 5.1.2 Kết xây dựng số biểu Thông ba cho địa phương Hà Giang 5.1.2 Kết nghiên cứu số mối quan hệ nhân tố điều tra lâm phần Trong lâm phần tương quan D/H mô hàm S: H = e3,673-17,724/D1.3 (5.1) Tương quan giưa Hg/H0 mô tả hàm Power (Y = b1*Xb2) : H0 = 1,707*Hg0,851 (5.2) 74 5.1.2.2 Kết lập biểu thể tích Thông ba Đã lập biểu thể tích có vỏ biểu thể tích không vỏ cho Thông ba lập cho tỉnh Hà Giang phương trình trương quan thể tích với đường kính chiều cao (bảng 4-9) xây dựng hướng dẫn sử dụng biểu thể tích Vcvo = 0,00006439(D2H)0,942 (5.3) Vkvo = -0,012 + 0,921*Vcvo (5.4) 5.1.2.3 Kết nghiên cứu lập biểu cấp đất Biểu cấp đất lập cho Thông ba Hà Giang có cấp đất theo phương pháp dựa vào suất tăng trưởng chiều cao tiết diện (Ph g) với số cấp đất Si tuổi 15, thông qua tương quan H0 với Hg lập biểu cấp đất theo H0 (Bảng 4-17) Xây dựng hướng dẫn xác định biểu cấp đất Thông ba Hà Giang Phương trình suất tăng trưởng chiều cao cấp đất lập sau: Cấp I: Ph = 177,723*A-1,325 (5.5) Cấp II: Ph = 177,255*A-1,274 (5.6) Cấp III: Ph = 193,985*A-1,335 (5.7) 5.1.2.4 Kết xây dựng mô hình lập biểu sản lượng Đề tài thiết lập mô hình sinh trưởng vận dụng mô hình sinh trưởng để lập biểu sản lượng sau: (1) Mô hình mật độ N N = 6685*Dg-0,725 (5.8) (2) Mô hình đường kính bình quân Dg = e(1,774 + 0,075.Hg) (5.9) (3) Mô hình trữ lượng M = 0,000423* A0,1578 N0,9602.Dg2,0154 Biểu sản lượng lập theo biểu ( 4.18; 4.19; 4.20) (5.10) 75 5.2 Tồn Do Thông ba Hoàng Su Phì lâm phần thuộc cấp tuổi I (1-5 tuổi), rừng trồng theo thời kỳ, không liên tục Nên số liệu tính toán không đại diện phản ánh hết quy luật sinh học theo tuổi nhân tố điều tra lâm phần Số lượng giải tích ô tiêu chuẩn, cấp tuổi cấp đất chưa nhiều nên kết nghiên cứu chưa phản ánh hết quy luật sinh trưởng lâm phần Rừng trồng nằm gần với khu dân cư, người dân thường vào rừng trích nhựa lấy củi dẫn tới số lâm phần không phản ánh chất sinh học, sinh thái học cá thể lâm phần Biến dạng trình trích nhựa lâm phần tuổi lớn nguyên nhân khiến hình dạng thân bị thay đổi, sinh trưởng đường kính bị ảnh hưởng 5.3 Kiến nghị Đề nghị có nghiên cứu cho Thông ba tuổi cao mở rộng nghiên cứu với số tỉnh miền núi phía Bắc Đề tài đề xuất bổ sung thêm số liệu cấp tuổi thiếu số vùng có điều kiện lập địa tương tự (Yên Bái) để có kết khách quan xác việc ứng dụng biểu Thông ba Hà Giang Miền Bắc Đề xuất thiết lập ô tiêu chuẩn định vị theo dõi, bảo vệ ô tiêu chuẩn định vị để mô tả xác chất sinh học lâm phần 76 Tài liệu tham khảo Tiếng Việt Nguyễn Trọng Bình, Vũ Thế Hồng Hoàng Xuân Y (2003), Lập biểu sinh trưởng sản lượng tạm thời cho rừng keo lai trồng loài, Trường Đại học Lâm nghiệp Nguyễn Trọng Bình, (2003), Lập biểu sinh trưởng sản lượng tạm thời cho rừng Keo lai trồng loài Trần Văn Con (2008), Báo cáo kết kiểm tra biểu thể tích, biểu cấp đất loài trồng vùng dự án trồng rừng KfW1 KfW3 Bắc Giang, Lạng Sơn Đông Triều - Quảng Ninh, Viện Khoa học Lâm nghiệp Việt Nam Lê Huy Cường, Trần Quốc Dũng, Trần Bạch Đằng Mưai Kỳ Vinh (2001), Tổng hợp hoàn thiện loại biểu số loài trồng rừng Việt Nam, Viện Điều tra Quy hoạch rừng Võ Đại Hải, Đặng Thịnh Triều, Nguyễn Hoàng Tiệp, Nguyễn Văn Bích Đặng Thái Dương (2009), Năng suất sinh khối hấp thụ carbon số dạng rừng trồng chủ yếu Việt Nam, Nhà xuất Nông nghiệp Võ Đại Hải, (2007), Điều tra đánh giá suất sinh trưởng loài trồng rừng chủ yếu dạng lập địa, làm sở đề xuất biện pháp lâm sinh cho trồng rừng phục vụ mục tiêu kinh doanh gỗ lớn cho xuất Viện Khoa học Lâm nghiệp Việt Nam Đồng Sĩ Hiền (1974), Lập biểu thể tích biểu độ thon đứng cho rừng Việt Nam, Nhà xuất Khoa học Kỹ thuật, Hà Nội Vũ Tiến Hinh (1998), Sản lượng rừng, Nhà xuất Nông nghiệp, Hà Nội, Vũ Tiến Hinh, Nguyễn Trọng Bình Hoàng Xuân Y (1996), Lập biểu trình sinh trưởng sản lượng rừng trồng keo tràm, Đại học Lâm nghiệp 77 10 Vũ Tiến Hinh Phạm Ngọc Giao (1997), Điều tra rừng, Nhà xuất Nông nghiệp, Hà Nội 11 Vũ Tiến Hinh, Hoàng Xuân Y, Phạm Ngọc Giao, Nguyễn Thị Bảo Lâm Nguyễn Trọng Bình (2000), Lập biểu sinh trưởng sản lượng cho ba loài cây: Sa mộc (Cunninghamia lanceolata - Hook), Thông đuôi ngựa (Pinus mưassosiana - Lamb), Mỡ (Mưanglietia glauca) tỉnh phía Bắc Đông bắc Việt Nam, Trường Đại học Lâm nghiệp 12 Vũ Tiến Hinh, Phạm Xuân Hoàn, Phạm Ngọc Giao Ngô Kim Khôi (2002), Lập biểu sản lượng rừng Quế Văn Yên - Yên Bái, Trường Đại học Lâm nghiệp 13 Vũ Tiến Hinh (2003), Sản lượng rừng, Nhà xuất Nông nghiệp, Hà Nội 14 Bảo Huy (1995), Dự đoán sản lượng rừng Tếch Đắk Lắk, Tạp chí Lâm nghiệp Số 15 Bảo Huy Đào Công Khanh (2008), Biểu sản lượng rừng trồng Trám trắng tỉnh Lạng Sơn, Bắc Giang, Quảng Ninh, Dự án "Các biện pháp đào tạo hỗ trợ quản lý rừng bền vững tỉnh Bắc Giang Lạng Sơn, KfW - MƯARD" 16 Trịnh Đức Huy, (1988), Dự đoán trữ lượng rừng suất gỗ đất trồng Bồ đề (Styrax tonkinensis) loài tuổi vùng trung tâm ẩm Bắc Việt Nam, Luận án PTS KHVN, Viện Khoa học Lâm nghiệp Việt Nam 17 Vũ Tiến Hưng, (2006), Nghiên cứu phương pháp hiệu chỉnh biểu trình sinh trưởng để xác định số tiêu sản lượng cho lâm phần Mỡ (Manglietia Glauca) Sa Mộc (Cunninghamia lenceolata) số tỉnh phía Bắc, Luận văn thạc sĩ khoa học Lâm nghiệp, Trường đại học Lâm nghiệp, Hà Tây 18 Đào Công Khanh, Bảo Huy, Đặng Văn Thuyết, Phan Minh Sáng, Bùi Thanh Hằng, Hoàng Văn Thắng Nguyễn Thanh Đạm (2001), Lập biểu trình sinh trưởng sản lượng cho rừng trồng loài Bạch đàn urophylla (Eucalyptus urophylla), Keo tai tượng (Accacia mưangium), Tếch (Tectona 78 grandis), Thông nhựa (Pinus merkusii) kiểm tra biểu sản lượng loài Đước (Rhizophora apiculata) Tràm (Mưalaleuca cajuputi), Viện Khoa học Lâm nghiệp 19 Nguyễn Ngọc Lung (1987a), Bàn lý thuyết chủ động điều khiển mật độ rừng theo mục tiêu điều chế, Tạp chí Lâm nghiệp 20 Nguyễn Ngọc Lung (1987b), Mô hình hóa trình sinh trưởng loài mọc nhanh để dự đoán sản lượng, Tạp chí Lâm nghiệp 21 Nguyễn Ngọc Lung Đào Công Khanh (1999), Nghiên cứu tăng trưởng sản lượng rừng trồng (áp dụng cho rừng Thông ba Việt Nam), Nhà xuất Nông nghiệp, TP.Hồ Chí Minh 22 Vũ Nhâm (1988), Lập biểu sản phẩm thương phẩm cho rừng Thông đuôi ngựa (Pinus mưassosiana - Lamb) kinh doanh gỗ mỏ vùng Đông bắc - Việt Nam, Viện Nghiên cứu Lâm nghiệp 23 Nguyễn Thị Tú Oanh (2002), Thiết lập số mô hình sinh trưởng sản lượng Keo Lai (Acacia mangium & Acacia auriculiformis), Luận văn thạc sỹ Lâm nghiệp, Hà Tây 24 Vũ Đình Phương (1975), Cơ sở xác định mật độ trồng phương thức tỉa thưa kinh doanh rừng Bồ đề trồng, Tập san Lâm nghiệp 25 Phan Minh Sáng Lưu Cảnh Trung (2006), Hấp thụ Các bon lâm nghiệp, Trong: Nguyễn Ngọc Bình (Ed,), Cẩm nang ngành Lâm nghiệp, Nhà xuất Nông nghiệp 26 Khúc Đình Thành (2003), Lập biểu sinh trưởng sản phẩm rừng Keo tai tượng (Accacia mưangium) kinh doanh gỗ trụ mỏ vùng Đông Bắc - Việt Nam, Trường Cao đẳng Nông Lâm Đông Bắc 27 Nguyễn Hải Tuất, Vũ Tiến Hinh Ngô Kim Khôi (2006), Phân tích Thống kê Lâm nghiệp, Nhà Xuất Nông Nghiệp, Hà Nội 28 Tổng công ty Giấy Việt Nam (2009), Quyết định giá mua gỗ nguyên liệu thô từ tháng 6/2009 79 29 Viện Điều tra Quy hoạch rừng (1995), Sổ tay Điều tra Quy hoạch rừng, Nhà xuất Nông nghiệp, Hà Nội 30 Viện Khoa học Lâm nghiệp Việt Nam, (2005), Dự án "Điều tra đánh giá xác định tập đoàn trồng rừng sản xuất có hiệu dạng lập địa chủ yếu vùng kinh tế lâm nghiệp toàn quốc", Hà Nội 31 Vụ KHCN & CLSP - Bộ Nông nghiệp PTNT (2003), Tiêu chuẩn ngành 04 - TCN - 66 - 2003: Biểu điều tra kinh doanh rừng trồng 14 loài chủ yếu, Nhà xuất Nông nghiệp, Hà Nội Tiếng nước Alder D (1980) Yield prediction, In: Forest volume estimation and yield prediction, FAO, Rome Amaro A, Reed D, Soares P (2003) Modelling Forest Systems, CABI Publishing; Assmann,E,, (1970) The principles of Forest yiel study, Pergamon Press 1970 (translation by Gardiner S,N) 506 trang; Avery TE, Burkhart HE (1994) Forest Measurements (4th edition), McGraw-Hill, Inc Evans J,, (1974) Som aspects of growth of Pinus patula in Swaziland, Commonwealth Forestry Review 53, 57-62pp Evans J,, (1992) Plantation Forestry in the tropics, Clarendon Press-Oford, Grafen A, Hails R (2002) Modern Statistics for the Life Sciences, Oxford University Press Husch B, Beers TW, Kershaw JA (2003) Forest Mensuration (4th edition), John Wiley & Son, Inc Jayaraman K (1999) A statistical manual for Forestry Research, FAO 10 Johnson JB, Omland KS (2004) Model selection in ecology and evolution, Trends in Ecology and Evolution, 19, 101-108 11 Monserud RA (2003) Evaluating Forest Models in a Sustainable Forest Management Context, Forest Biometry, Modelling and Information Sciences, 1, 35-47 80 12 Pretzsch H (2001) Models for Pure and Mixed Forests, In: The Forests Handbook, (ed Evan J), pp, 210-228, Blackwell Science 13 Pretzsch H (2009) Forest Dynamics, Growth and Yield, Springer 14 Pretzsch H, Grote R, Reineking B, Rotzer T, Seifert S (2008) Models for Forest Management: A European Perspective, Annals of Botany, 101, 10651087 15 Phan Minh Sáng, Schmidt, S, Lamb, D, (2008), Optimal rotation for timber production and carbon sequestration of Eucalytpus urophylla and Pinus merkusii plantations in Vietnam, Trong, Linking Forest Inventory and Optimization – LIFO 2008, Munich, Germany 16 Phan Minh Sáng (2008), Carbon sequestration and soil fertility of tropical tree plantations and secondary forests in Vietnam, Department of Botany, University of Queensland 17 Phan Minh Sáng, Yates, D,, Erskine, P,, Lamb, D, Schmidt, S, (2009), Simulating productivity of commercial Acacia mangium plantations using the 3-PG model (under review), Forest Ecology and Management 18 Van Laar A, Akca A (2007) Forest Mensuration, Springer 19 Vanclay JK (1999) Modelling Forest Growth and Yield, CAB International 20 Vanclay JK, Skovsgaard JP (1997) Evaluating of Forest Growth Models, Ecological Modelling, 98, 1-12 21 Veiga, V,P,1964a: “Five years growth and yield of residual forest in Benguet, Mountain Provice” Bur, For, Philippnies Res, Note No 68” 22 Wenk G, Forstliche Ertraskunde, T,U, Dresden 1977 (300 trang) 23 Wenk G, Holzmesskunde, TU Dresen 1973, 250 trang 24 Zeide B (1993), Analysis of Growth Equations, Forest Science, 39, 594 – 616 81 ii MỤC LỤC Trang phụ bìa Lời cảm ơn ………………………………………………………………………….i Mục lục …………………………………………………………………………… ii Danh mục bảng ……………………………………………………………… iv Danh mục hình ………………………………………………………………… v ĐẶT VẤN ĐỀ .1 Chương .3 TỔNG QUAN VỀ VẤN ĐỀ NGHIÊN CỨU .3 1.1 Ngoài nước 1.1.1 Nghiên cứu lập bảng biểu điều tra rừng giới 1.1.2 Phương pháp kiểm tra bảng biểu .8 1.2 Trong nước 1.2.1 Nghiên cứu lập bảng biểu điều tra rừng Việt Nam 1.2.2 Phương pháp kiểm tra bảng biểu 16 1.2.3 Bảng biểu điều tra cho kinh doanh rừng Thông ba Việt Nam 17 Chương 21 MỤC TIÊU, ĐỐI TƯỢNG, PHẠM VI, NỘI DUNG VÀ 21 PHƯƠNG PHÁP NGHIÊN CỨU .21 2.1 Mục tiêu nghiên cứu 21 2.2 Đối tượng phạm vi nghiên cứu 21 2.2.1 Đối tượng nghiên cứu 21 2.2.2 Phạm vi nghiên cứu 21 2.3 Nội dung nghiên cứu 22 2.3.1 Nghiên cứu khả ứng dụng số biểu điều tra Thông ba lập cho vùng Tây Nguyên Hà Giang 22 2.3.2 Nghiên cứu số mối quan hệ nhân tố điều tra lâm phần .22 2.3.3 Bước đầu nghiên cứu xây dựng số bảng, biểu điều tra rừng Thông ba Hà Giang 22 2.4 Phương pháp nghiên cứu 23 2.4.1 Phương pháp điều tra 23 82 2.4.2 Phương pháp xử lý số liệu .25 Chương 32 ĐIỀU KIỆN TỰ NHIÊN KHU VỰC NGHIÊN CỨU .32 3.1 Vị trí địa lý 32 3.2 Địa hình, địa mạo đất đai 32 3.3 Khí hậu .33 3.4 Chế độ thuỷ văn 33 Chương 34 KẾT QUẢ VÀ PHÂN TÍCH KẾT QUẢ 34 4.1 Nghiên cứu khả ứng dụng số biểu điều tra Thông ba lập cho vùng Tây Nguyên Hà Giang 34 4.1.1 Kiểm nghiệm khả ứng dụng biểu thể tích 34 4.1.2 Kiểm kiểm nghiệm khả ứng dụng biểu cấp đất 36 4.1.3 Kiểm tra biểu sản lượng 40 4.2 Nghiên cứu số mối quan hệ nhân tố điều tra lâm phần 40 4.2.1 Nghiên cứu quan hệ đường kính chiều cao 40 4.2.2 Nghiên cứu mối quan hệ chiều cao theo tiết diện bình quân (Hg) chiều cao tầng trội (H0) 42 4.2.3 Nghiên cứu quan hệ thể tích có vỏ (Vcvo ) thể tích không vỏ (Vkvo) .43 4.3 Bước đầu nghiên cứu xây dựng số bảng, biểu điều tra rừng Thông ba Hà Giang 45 4.3.1 Nghiên cứu lập biểu thể tích hai nhân tố 45 4.3.2 Nghiên cứu cứu lập biểu cấp đất .49 4.3.3 Nghiên cứu lập biểu sản lượng .64 Chương 73 KẾT LUẬN – TỒN TẠI – KHUYẾN NGHỊ 73 5.1 Kết luận 73 5.2 Tồn 75 5.3 Kiến nghị 75 Tài liệu tham khảo .76 Phụ lục 83 iv DANH MỤC CÁC BẢNG TT Tên Trang Bảng 2-1: Tổng hợp dung lượng mẫu điều tra 25 Bảng 4-1 Bảng kết kiểm tra biểu thể tích cá lẻ 35 Bảng 4-2: Bảng tổng hợp kiểm tra dạng liên hệ D/H 40 Bảng 4-3: Bảng kiểm tra tồn tham số phương trình S .41 Bảng 4-4: Bảng kiểm tra tồn dạng liên hệ Hg với H0 .42 Bảng 4-5: Kết kiểm tra dạng liên hệ Vkvo/Vcvo .43 Bảng 4-6: Kết kiểm tra tồn tham số phương trình 44 Biểu 4-7: Kết tính toán phương trình thử nghiệm để lập biểu 45 Bảng 4-8: Bảng kiểm tra sai số biểu thể tích Vcvo theo phương trình 46 Bảng 4-9: Biểu thể tích nhân tố Thông ba trồng Hà Giang 48 Bảng 4-10: Kết kiểm tra phương trình xác lập 51 Bảng 4-11 : Biểu cấp đất lập theo suất tăng trưởng chiều cao (Ph) 51 Bảng 4-12: Kết xác định tham số phương trình sinh trưởng .53 Bảng 4-13 : Phương trình cấp đất theo a chung, b thay đổi .55 Bảng 4-14: Biểu đồ cấp đất tính theo phương pháp a chung b thay đổi 55 Bảng 4-15: Biểu cấp đất lập theo phương pháp Affill .57 Bảng 4-16: Kiểm tra hệ số hồi quy đường cong cấp đất đường cong sinh trưởng chiều cao thực nghiệm 61 Bảng 4-17: Biểu cấp đất H0 lập theo suất tăng trưởng chiều cao Ph .62 Bảng 4-18: Biểu trình sinh trưởng cấp đất I 70 Bảng 4-19: Biểu trình sinh trưởng cấp đất II 71 Bảng 4-20: Biểu trình sinh trưởng cấp đất III 72 84 v DANH MỤC CÁC MỤC HÌNH TT Tên hình Trang Hình 4-1: Biểu đồ đám mây điểm kết kiểm tra sai số thể tích 35 Hình 4-3: Biểu đồ kiểm tra biểu cấp đất II 38 Hình 4-4: Biểu đồ kiểm tra biểu cấp đất III .39 Hình 4-3: Biểu đồ mô tả đám mây điểm quan hệ D/H 41 Hình 4-4 : Mô tả tương quan Vkvo/Vcv .44 Hình 4-7: Biểu đồ đám mây điểm chiều cao thực nghiệm Hg 49 Hình 4-8 : Biểu đồ cấp đất lập theo suất tăng trưởng chiều cao Ph 52 Hình 4-9: Mô sinh trưởng chiều cao Hg hàm toán học 53 Hình 4-10: Biểu đồ cấp đất theo phương pháp a chung, b thay đổi 56 Hình 4-11: Biểu đồ cấp đất lập theo phương pháp Affill 58 Hình 4-12: Biểu đồ kiểm tra cấp đất lập theo Ph .59 Hình 4-13: Biểu đồ kiểm tra cấp đất lập theo phương pháp a chung, b thay đổi 60 Hình 4-14: Biểu đồ kiểm tra cấp đất lập theo phương pháp Affill 60 Hình 4-15: Biểu đồ cấp đât H0 lập theo Ph 63 85 NHỮNG KÝ HIỆU SỬ DỤNG TRONG LUẬN VĂN A : Tuổi lâm phần, tuổi D1.3 : Đường kính ngang ngực, đườn kính thân đo độ cao 1,3m Dt : Đường kính tán G : Tổng tiết diện ngang Hvn : Chiều cao vút H0 : Chiều cao tầng ưu Hg : Chiều cao có tiết diện bình quân HF : Hình cao N : Mật độ M : Trữ lượng lâm phần ∆M : Tăng trưởng bình quân chung trữ lượng ZM : Tăng trưởng thường xuyên hàng năm trữ lượng ∆% : Sai số tương đối kv : Hệ số tỉa thưa theo trữ lượng kg : Hệ số tỉa thưa theo tiết diện kh : Hệ số tỉa thưa theo chiều cao n : Dung lượng quan sát R2 : Hệ số xác định r : Hệ số tương quan St : Tổng diện tích tán lâm phần t05 Tiêu chuẩn t Student ttinh : Tiêu chuẩn tính kiểm tra tồn tham số phương trình Sig : Tính xác xuất khoảng tin cậy tiêu chuẩn kiểm tra ÔTC : Ô tiêu chuẩn ... dụng bảng biểu điều tra, kinh doanh rừng Thông ba lập Tây Nguyên cho rừng Thông ba Hà Giang - Xây dựng biểu thể tích, biểu cấp đất, biểu sản lượng kinh doanh rừng Thông ba tỉnh Hà Giang 2.2 Đối... nghiệp: Tên đề tài: Lập biểu thể tích, cấp đất sản lượng cho Thông ba (Pinus kesiya Royle ex. Gordon) Hà Giang Sau thời gian tiến hành làm đề tài tốt nghiệp đến luận văn hoàn thành Cho phép bày tỏ... biểu điều tra rừng Thông ba Tây Nguyên Từ yêu cầu thực tiễn trên, tiến hành thực đề tài: Lập biểu thể tích, cấp đất sản lượng cho Thông ba (Pinus kesiya Royle ex. Gordon) Hà Giang Chương TỔNG