(Luận văn) xây dựng cơ sở khoa học cho việc điều tra thể tích thân cây từ kích thước gốc chặt của một số loài cây ở rừng tự nhiên vùng bắc trung bộ

Bộ giáo dục đào tạo Bộ nông nghiệp PTNT Tr-ờng đại học lâm nghiệp Vũ QUốC PHòNG lu Xây dựng sở khoa học cho việc điều tra thĨ tÝch an n va th©n c©y tõ kÝch th-íc gốc chặt số loài p ie gh tn to rừng tự nhiên vùng bắc trung w d oa nl Chuyên ngành : Lâm học Mã số : 60.62.60 u nf va an lu ll Luận văn thạc sĩ khoa học Lâm nghiệp oi m z at nh z NGƯỜI HƯỚNG DẪN KHOA HỌC @ m co l gm GS.TS VŨ TIẾN HINH an Lu n va Hµ néi - 2011 ac th si i LỜI CẢM ƠN Với mục tiêu góp phần hồn thiện sở khoa học thực tiễn vấn đề xác định thể tích gỗ thân từ kích thước gốc chặt số loài khai thác phổ biến Việt Nam, tiến hành thực đề tài: “Xây dựng sở khoa học cho việc điều tra thể tích thân từ kích thước gốc chặt số loài rừng tự nhiên vùng Bắc Trung bộ” Sau thời gian thực hiện, hướng dẫn tận tình GS.TS lu an Vũ Tiến Hinh, với giúp đỡ thầy cô giáo bạn đồng n va nghiệp, đến đề tài hoàn thành tn to Nhân dịp này, cho phép tơi bày tỏ lịng cảm ơn chân thành tới GS.TS gh Vũ Tiến Hinh, thầy cô môn Điều tra quy hoạch rừng – Trường p ie Đại học Lâm nghiệp bạn đồng nghiệp giúp đỡ tơi hồn thành luận văn w oa nl Mặc dù có cố gắng, lực thân thời d gian nghiên cứu hạn chế, nên kết đạt đề tài chắn lu va an khơng tránh khỏi thiếu sót, mong quan tâm, đóng góp ý u nf kiến thầy cô giáo, bạn đồng nghiệp người quan tâm đến ll vấn đề để luận văn hồn thiện m oi Tơi xin cam đoan, số liệu thu thập, kết tính tốn luận văn z at nh hoàn toàn trung thực trích dẫn rõ ràng z Tơi xin chân thành cảm ơn! @ gm Hà Nội, tháng 10 năm 2011 m co l Tác giả an Lu Vũ Quốc Phòng n va ac th si ii MỤC LỤC Trang phụ bìa Trang Lời cảm ơn i Mục lục ii Danh mục ký hiệu ivv Danh mục bảng v ĐẶT VẤN ĐỀ lu Chương 1: TỔNG QUAN VẤN ĐỀ NGHIÊN CỨU an va 1.1 Trên giới .2 n 1.2 Trong nước .7 gh tn to Chương 2: MỤC TIÊU, ĐỐI TƯỢNG, PHẠM VI, NỘI DUNG VÀ PHƯƠNG ie PHÁP NGHIÊN CỨU 15 p 2.1 Mục tiêu nghiên cứu 15 nl w 2.2 Đối tượng nghiên cứu 15 oa 2.3 Phạm vi nghiên cứu .15 d 2.4 Nội dung nghiên cứu: 15 lu va an 2.4.1 Xác định đặc trưng thống kê đường kính chiều cao gốc chặt 15 u nf 2.4.2 Xác lập quan hệ thể tích thân với đường kính chiều cao gốc ll chặt 16 m oi 2.4.3 Xác định thể tích thân phương trình thể tích biểu thể tích z at nh lập sẵn 16 2.4.4 Chọn phương pháp xác định thể tích thân thơng qua kích thước gốc z @ chặt 16 gm 2.5 Phương pháp nghiên cứu .16 m co l 2.5.1 Quan điểm phương pháp luận 16 2.5.2 Phương pháp thu thập số liệu 17 an Lu 2.5.3 Phương pháp xử lý số liệu 17 n va ac th si iii Chương 3: KẾT QUẢ NGHIÊN CỨU VÀ THẢO LUẬN 25 3.1 Khái quát số liệu nghiên cứu 25 3.2 Đặc điểm kích thước gốc chặt 27 3.3 Xác lập quan hệ thể tích thân với đường kính gốc chặt chiều cao gốc chặt .31 3.4 Xác định thể tích thân sở sử dụng phương trình thể tích biểu thể tích lập sẵn .37 3.4.1 Kiểm tra phụ thuộc đường kính ngang ngực D1.3 vào đường kính gốc chặt Dgc, chiều cao gốc chặt Hgc, đường kính gốc D0 37 lu an 3.4.2 Xác định quan hệ đường kính ngang ngực D1.3 với đường kính gốc va chặt Dgc đường kính gốc D0 39 n 3.4.3 Xác lập quan hệ đường kính ngang ngực D1.3 với đường kính gốc to 3.4.4 Kiểm tra phụ thuộc chiều cao vút Hvn vào đường kính gốc p ie gh tn chặt Dgc chiều cao gốc chặt Hgc 43 chặt Dgc, chiều cao gốc chặt Hgc, đường kính gốc D0 45 nl w 3.4.5 Xác lập quan hệ chiều cao vút Hvn với đường kính gốc chặt d oa Dgc chiều cao gốc chặt Hgc 47 lu 3.4.6 Xác định thể tích gỗ thân bị chặt thơng qua đường va an kính ngang ngực D 1.3 chiều cao vút H 49 u nf 3.4.7 Xác định thể tích gỗ thân sở sử dụng biểu thể tích ll lập sẵn 55 m oi 3.5 Chọn phương pháp xác định thể tích gỗ thân thơng qua kích thước z at nh gốc chặt 57 KẾT LUẬN - TỒN TẠI - KIẾN NGHỊ 60 z gm @ Kết luận .60 Tồn .61 l Kiến nghị 62 an Lu PHỤ LỤC m co TÀI LIỆU THAM KHẢO n va ac th si iv DANH MỤC CÁC KÝ HIỆU D1.3 : Đường kính thân vị trí 1.3m Dgc : Đường kính gốc chặt (đường kính mặt cắt gốc chặt) D0 : Đường kính gốc vị trí sát mặt đất f : Là hình số (độ tròn đầy) Hvn : Chiều cao vút Hgc : Chiều cao gốc chặt khai thác gỗ l : Chiều dài xúc gỗ lu an Ntt: Số tính tốn n va Nkt: Số kiểm tra tn to PP1 : Phương pháp xác định tương quan V với Dgc Hgc PP3 : Phương pháp tra biểu thể tích hai nhân tố D1.3 Hvn p ie gh PP2 : Phương pháp xác định tương quan V với D1.3 Hvn w R2 : Hệ số xác định oa nl R : Hệ số tương quan d S% : Hệ số biến động lu va an S : Sai tiêu chuẩn u nf Sbi : Phương sai hệ số hồi quy bi ll V : Thể tích gỗ thân m oi Vt : Thể tích thực tính theo cơng thức kép tiết diện bình quân m co l gm ∆ v%: Sai số thể tích trung bình @  : Hệ số thon z Xbq : Giá trị bình quân z at nh Vlt : Thể tích theo lý thuyết ∆%ΣV: Sai số tổng thể tích tập hợp thân an Lu ∆ v% : Sai số tương đối cá lẻ n va ac th si v DANH MỤC CÁC BẢNG Tên bảng TT Trang 3.1: Khái quát số liệu nghiên cứu 26 3.2: Tổng hợp đặc trưng thống kê chiều cao gốc chặt 28 3.4: Kết phân tích tương quan V với Dgc Hgc 32 3.5: Phương trình tương quan V với Dgc Hgc chọn 36 3.6: Kiểm tra phụ thuộc D1.3 vào Dgc, Hgc, D0 38 lu an 3.7: Tương quan D1.3 với Dgc D0 40 va n 3.8: Kiểm tra phương trình tương quan D1.3 Dgc 42 tn to 3.9: Kết tính tương quan D1.3 với Dgc Hgc 44 ie gh 3.10: Kết kiểm tra phụ thuộc Hvn vào Dgc, Hgc, D0 46 p 3.11: Kết tính tương quan Hvn với Dgc Hgc 48 nl w 3.12: Kết phân tích tương quan V với Hvn D1.3 50 d oa 3.13: Phương trình tương quan V với D1.3 Hvn chọn 54 an lu 3.14: Kết tính sai số thể tích thân xác định từ biểu thể tích lập theo tổ f01 56 ll u nf va 3.15: Kết xác định sai số thể tích phương pháp 58 oi m z at nh z m co l gm @ an Lu n va ac th si ĐẶT VẤN ĐỀ Rừng tự nhiên Việt Nam nói chung vùng Bắc Trung nói riêng có nhiều lồi với nhiều hình dạng kích thước khác Sự đa dạng phong phú có ý nghĩa quan trọng việc phát triển kinh tế, bảo tồn đa dạng sinh học trì ổn định hệ sinh thái Tuy nhiên, việc khai thác rừng tự nhiên Bắc Trung năm vừa qua diễn biến phức tạp, đặc biệt vụ phá rừng với quy mô lớn nhiều lu loài gỗ quý quan chức phát cịn lại vết tích an va gốc chặt Điều gây nhiều khó khăn cho cơng tác điều tra, kiểm kê, n đánh giá lượng gỗ bị hàng năm, đòi hỏi phải tốn nhiều công sức tiền gh tn to Trong năm gần đây, Việt Nam có nhiều nghiên cứu lập p ie biểu thể tích cho lồi chủ yếu Tuy nhiên, phần lớn nghiên cứu tập trung vào đối tượng rừng trồng mà có điều kiện quan tâm đến oa nl w lồi nhóm lồi rừng tự nhiên Có dừng lại d việc lập biểu thể tích hai nhân tố (thường đường kính D1.3 chiều cao an lu Hvn), hai nhân tố phải xác định gián tiếp từ kích thước gốc chặt u nf va Xuất phát từ yêu cầu thực tế nói trên, tiến hành thực đề tài: “Xây dựng sở khoa học cho việc điều tra thể tích thân từ kích ll oi m thước gốc chặt số loài rừng tự nhiên vùng Bắc Trung bộ” z at nh Kết đề tài góp phần vào việc đề xuất phương pháp xác định thể tích gỗ thân cho số lồi chủ yếu rừng tự nhiên vùng Bắc z m co l thành so với phương pháp khác gm @ Trung bộ, đồng thời đảm bảo độ xác cần thiết giảm chi phí giá an Lu n va ac th si Chương 1: TỔNG QUAN VẤN ĐỀ NGHIÊN CỨU Mục đích đề tài lựa chọn phương pháp xác định thể tích gỗ thân thơng qua kích thước gốc chặt mà đảm bảo độ xác cần thiết, phần tổng quan chủ yếu đề cập đến phương pháp xác định thể tích có mơ hình dự đốn mối quan hệ kích thước gốc chặt với nhân tố tham gia cấu thành thể tích gỗ thân D1.3, Hvn , làm sở cho việc định hướng lựa chọn phương pháp nghiên lu cứu an n va 1.1 Trên giới tn to Các phương pháp tính thể tích ie gh Cho đến nay, giới có nhiều cơng trình nghiên cứu liên quan đến p lập biểu thể tích thân thiết lập mối quan hệ nhân tố dễ xác định với nl w nhân tố cấu thành thể tích thân cây, giới hạn cho phép đề tài đề cập d oa đến số kết tiêu biểu có liên quan đến nội dung nghiên cứu an lu Ngay từ cuối kỷ XIX nhà lâm học biết sử dụng công thức u nf va hình học (viên trụ, paraboloid bậc cụt, đơn tiết diện giữa, đơn tiết diện bình quân, Simpson, Hostfeild,…) để đo tính thể tích xúc gỗ sản phẩm cá lẻ ll m oi Sang đầu kỷ XX xuất nghiên cứu điều tra z at nh gỗ tròn Trước hết tác giả xây dựng loại biểu thể tích hình viên trụ, thực chất bảng tính sẵn để tiện cho việc áp dụng cơng thức hình z l số vị trí xúc gỗ gm @ học nêu Tuy nhiên, loại biểu khó sử dụng phải biết đường kính m co Giai đoạn 1906 – 1908 Cruidener, giám đốc sở lâm nghiệp Hồng gia có chứa phần bạnh gốc khơng gồm bạnh gốc) an Lu (nước Nga) lập biểu thể tích gỗ trịn cho lồi (phân biệt thành gỗ tròn n va ac th si Do biểu ban đầu lập theo phương pháp thực nghiệm, dung lượng mẫu có hạn, nên cịn mắc sai số lớn, đặc biệt cỡ cực đoan Mendeleep D.I(1899), Belanovxki I.G(1917) Wimmenauer K(1918) đặt mục tiêu xác định hình dạng đường sinh biểu thị phương trình tốn học, xem đường kính (Y) hàm chiều cao (x): Y = F(x) đề nghị biểu thị hàm số phương trình bậc hai, bậc ba bậc bốn (Theo Đồng Sĩ Hiền 1974) Theo (Skindele LeMay, 2006; Đồng Sĩ Hiền, 1974; Husch et al, lu an 2003), mơ hình tốn học thể tích thân xem hàm n va biến độc lập: đường kính, chiều cao hình số [26] Nó biểu thị tn to dạng phương trình: (1.1) Petrovxki V.S (1963, 1964) Liên Xô cũ, biểu thị quan hệ đường p ie gh V = F(D, H, f) nl w kính lấy vị trí với khoảng cách (L) từ đường kính đến gốc d oa phương trình Parabol sau: lu va an X2 = 2.P.(y - h) ll u nf Trong đó: + P thông số tiêu đỉnh đường sinh oi m +X, y toạ độ Parabol, h chiều dài thân bớt 1m H z at nh Khi thể tích thân tính theo công thức: (1.2) z V =   X dl   M d 052 H m co l gm Trong M tuỳ thuộc vào lồi @ an Lu n va ac th si Turxki (theo Anoutchin [1971]) dùng phương pháp biểu đồ để nắn số liệu Cruidener hiệu chỉnh thành biểu lấy tên tác giả thừa nhận làm tiêu chuẩn quốc gia sử dụng Liên xô cũ ngày Tuy nhiên, sai số phương pháp biều đồ cao so với loại biểu thể tích hình viên trụ trước Ở cộng hòa Séc Korsum (theo Anoutchin [1971]) cho thể tích gỗ trịn quan hệ chặt chẽ với chiều dài sản phẩm theo dạng phương trình: v  k  lm (1.3) lu an Và xét nhân tố đường kính thì: va n v  k  lm  d n (1.4) to gh tn Nghiên cứu thực nghiệm Korsum kết luận k, m, n khác p ie tính tốn cho đối tượng khác nên việc ứng dụng tương quan w tương đối khó khăn vào giai đoạn kỷ XX oa nl Vấn đề nghiên cứu hình dạng thân phục vụ cho việc lập biểu thể d tích thực ý vào năm cuối thể kỷ XX Các tiêu lu va an tác giả đặc biệt quan tâm độ thon tuyệt đối, độ thon bình quân u nf (Anoutchin [1971]), độ thon tương đối (Zakharov [1967]) hình số gỗ trịn ll (Dementiev) oi m z at nh Anoutchin dựa vào tài liệu 4000 súc gỗ trịn xác định độ thon bình qn phụ thuộc chặt chẽ vào đường kính đầu sản phẩm theo phương z trình: (1.5) l gm @ s  0,39  0,021d m co Đồng thời xác định s dao động từ 0.77 đến 1.87, bình quân = 0.96 Hệ thay đổi từ 15 đền 55cm an Lu số biến động từ 26 đến 47% bình quân 38%, tương ứng với đường kính n va ac th si 52 TT 21 22 23 lu an n va 24 Loài gh tn to Dạng PT 2.21 2.22 Táu nước 2.23 2.24 2.21 2.22 Trâm Móc 2.23 2.24 2.21 2.22 Trám Trắng 2.23 2.24 2.21 2.22 Trường Sâng 2.23 2.24 2.21 2.22 Trường Vải 2.23 2.24 2.21 2.22 Kiền Kiền 2.23 2.24 2.21 2.22 Chò nâu 2.23 2.24 2.21 2.22 Dầu rái 2.23 2.24 2.21 2.22 Trám Chủa 2.23 2.24 2.21 2.22 Ươi 2.23 2.24 25 p ie d oa nl w 26 ll oi m z at nh z m co l gm @ 30 u nf 29 va 28 an lu 27 Phương trình V1 = 0.00013614*D^1.5676*H^1.1606 V2 = 0.1392+0.00003596*D^2*H V3 = -0.3209+0.0421*H+0.00003031*D^2*H V4 = 0.00009497*(D^2*H)^0.9170 V1 = 0.00004635*D^2.2404*H^0.6154 V2 = -0.0627+0.00004085*D^2*H V3 = 0.6892-0.0479*H+0.00004305*D^2*H V4 = 0.00002633*(D^2*H)^1.0373 V1=0.000126*D^1.9509*H^0.6134 V2=0.2075+0.0000297*D^2*H V3=0.8066-0.0361*H+0.0000312*D^2*H V4=0.0000794*(D^2*H)^0.9203 V1= 0.00001997*D^2.1278*H^1.0571 V2 = -0.1260 + 0.00004178*D^2*H V3 =-0.1376+0.0008*H+0.00004172*D^2*H V4 = 0.00001988*(D^2*H)^1.0625 V1 = 0.00009601*D^1.7803*H^0.9888 V2 = 0.0952+0.00003715*D^2*H V3 = 0.1340-0.0029*H+0.00003742*D^2*H V4 = 0.00009967*(D^2*H)^0.9130 V1 = 0.000040*D^2.1185*H^0.8203 V2 = 0.0077+0.0000362*D^2*H V3 = 0.5626-0.0285*H+0.0000381*D^2*H V4= 0.0000392*(D^2*H)^0.9930 V1 = 0.000067*D^1.7389*H^1.1305 V2 = 0.7201+0.000029*D^2*H V3 = -2.9782+0.1507*H+0.000025*D^2*H V4= 0.000097*(D^2*H)^0.9111 V1 = 0.00022*D^1.9969*H^0.4505 V2 = 0.7376+0.0000276*D^2*H V3 = 0.3233+0.0158*H+0.0000272*D^2*H V4= 0.0000982*(D^2*H)^0.9079 V1 = 0.000096*D^1.8138*H^0.9195 V2 = 0.5033+0.000029*D^2*H V3 = -0.5779+0.0522*H+0.000027*D^2*H V4= 0.0000969*(D^2*H)^0.9095 V1 = 0.000046*D^2.0381*H^0.8582 V2 = -0.0582+0.0000339*D^2*H V3 = 0.7804-0.03657*H+0.0000358*D^2*H V4= 0.0000391*(D^2*H)^0.9861 an Lu R2 0.92 0.90 0.92 0.90 0.94 0.94 0.95 0.93 0.93 0.92 0.92 0.93 0.94 0.94 0.94 0.94 0.98 0.98 0.98 0.98 0.97 0.95 0.95 0.97 0.99 0.98 0.97 0.98 0.99 0.99 0.99 0.98 0.99 0.98 0.98 0.99 0.96 0.95 0.95 0.96 R 0.96 0.95 0.96 0.95 0.97 0.97 0.97 0.96 0.96 0.96 0.96 0.96 0.97 0.97 0.97 0.97 0.99 0.99 0.99 0.99 0.98 0.98 0.98 0.98 Sig 4E-21 1E-20 4E-21 7E-21 5E-18 2E-18 2E-21 5E-13 4E-22 2E-22 4E-21 3E-23 5E-23 3E-24 8E-23 2E-24 1E-32 2E-34 9E-33 4E-34 1E-23 2E-22 5E-21 5E-25 Ntt 40 40 40 40 40 40 40 40 40 40 40 40 40 40 40 40 40 40 40 40 34 34 34 34 0.99 5E-19 31 0.99 5E-20 31 0.99 4E-19 31 0.99 2E-20 31 0.99 0.99 0.99 0.99 0.99 0.99 0.99 0.99 0.98 0.97 0.98 0.98 2E-30 5E-33 3E-31 4E-31 6E-30 9E-29 1E-27 9E-32 5E-21 5E-21 9E-20 2E-22 35 35 35 35 35 35 35 35 32 32 32 32 n va ac th si 53 Kết nghiên cứu cho thấy: - Giữa thể tích gỗ thân với đường kính D1.3 Hvn lồi ln tồn mối quan hệ với hệ số xác định dao động từ mức chặt đến chặt (0.83 ÷ 0.99) - Xác suất kiểm tra tồn hệ số R2 nhỏ 0.05, chứng tỏ thực tồn hệ số R2 tổng thể (R2 > 0) Qua phân tích tương quan V với D1.3 Hvn cho thấy dạng phương trình áp dụng thực tế để xác định thể tích gỗ lu an thân cho lồi Tuy nhiên, để có sở lựa chọn phương trình phù hợp va n nhất, đề tài tiến hành tính sai số xác định thể tích theo dạng phương trình tn to cho loài từ kiểm tra Kết tính sai số dạng Qua phụ biểu 10 cho thấy phương trình có số lần mắc sai số (-) p ie gh phương trình tổng hợp phụ biểu 10 oa nl w (+) cân nhau, tỷ lệ 666/690, chứng tỏ phương trình khơng mắc sai số hệ thống Số lần sai số > 20% chiếm 6% tổng số lần kiểm tra d an lu tính chung cho lồi, sai số chấp nhận điều tra thể u nf va tích riêng lẻ Sai số thể tích trung bình 8.6% sai số tổng thể tích kiểm tra chung cho dạng phương trình 3.5%, với sai số đáp ll oi m ứng độ xác cơng tác điều tra rừng z at nh Để thuận tiện cho công tác điều tra tính tốn, lồi cần chọn dạng phương trình tốt nhất, đảm bảo độ tin cậy Căn vào mức z @  gm độ phức tạp phương trình, sai số thể tích trung bình  v % , sai số tổng thể l tích ∆%ΣV, hệ số xác định R2 chênh lệch sai số âm, dương m co phương trình lồi, đề tài chọn phương trình biểu thị tốt mối an Lu quan hệ thể tích V với D1.3 Hvn Kết cụ thể sau: n va ac th si 54 Bảng 3.13: Phương trình tương quan V với D1.3 Hvn chọn lu an n va Loài 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 Bộp Chẹo Dẻ đỏ Dẻ Trắng Nang Táu Muối Vạng Trứng Ngát Dẻ Bốp Dẻ Cau Gội Tẻ Vàng Tâm Bời Lời Bộp Vàng Cây Chủa Cây Trín Gội Nếp Lim Xanh Ràng Ràng Re Đá Táu Nước Trâm Móc Trám Trắng Trường Sâng Trường Vải Kiền Kiền Chò nâu Dầu rái Trám Chủa Cây Ươi p ie gh tn to TT Phương trình d oa nl w V2 = -0.1691+0.00003614*D^2*H V2 = 0.3006+0.00002452*D^2*H V2 = 0.0627+0.00002821*D^2*H V2=-0.0747+0.00003081*D^2*H V2 = 0.0532+0.00003390*D^2*H V2 = -0.3036+0.00004027*D^2*H V2=-0.0226+0.00002812*D^2*H V2 = 0.0827+0.00002543*D^2*H V2 = 0.1146+0.0000314*D^2*H V2 = 0.0471+0.00002648*D^2*H V2 = -0.1206+0.00003213*D^2*H V2 = -0.0962+0.0000343*D^2*H V2 = 0.0365+0.00003672*D^2*H V2=0.0335+0.00003989*D^2*H V2 = 0.3857+0.00002999*D^2*H V2 = 1.2261+0.00002322*D^2*H V2 = 0.1745+0.00003461*D^2*H V2 = 0.2119+0.00004132*D^2*H V2 = 0.1700+0.00003794*D^2*H V2 = -0.02064+0.00004321*D^2*H V2 = 0.1392+0.00003596*D^2*H V2 = -0.0627+0.00004085*D^2*H V2=0.2075+0.0000297*D^2*H V2 = -0.1260 + 0.00004178*D^2*H V2 = 0.0952+0.00003715*D^2*H V2 = 0.0077+0.0000362*D^2*H V2 = 0.7201+0.000029*D^2*H V2 = 0.7376+0.0000276*D^2*H V2 = 0.5033+0.000029*D^2*H V2 = -0.0582+0.0000339*D^2*H ll u nf va an lu oi m z at nh R2 Nkt % Số lần sai số  ∆%  % max + ΣV v