1. Trang chủ
  2. » Tất cả

Luận văn tốt nghiệp tìm hiểu một số quy luật kết cấu cơ bản lâm phần keo lai (acacia hybrid), phục vụ công tác điều tra rừng

60 3 0

Đang tải... (xem toàn văn)

Tài liệu hạn chế xem trước, để xem đầy đủ mời bạn chọn Tải xuống

THÔNG TIN TÀI LIỆU

PHẦN ĐẶT VẤN ĐỀ 1.1 TÍNH CẤP THIẾT CỦA VẤN ĐỀ NGHIÊN CỨU Trong năm gần nghiên cứu lai giống sử dụng giống lai mối quan tâm nhà chọn giống Nông, Lâm nghiệp đặc điểm bật giống lai có ưu lai đời lai F1 Chính lồi Keo, lai nhân tạo (thụ phấn có kiểm sốt) Giáo sư Lê Đình Khả, Kỹ sư Nguyễn Việt Cường cộng tác viên tiến hành năm 1997 - 2000 tạo số tổ hợp lai thuận nghịch lai trở lại làm nguồn vật liệu để khảo nghiệm có đánh giá bước đầu dựa khả sinh trưởng tiêu quan trọng định suất dễ dàng quan sát để xác định biểu ưu lai Keo tràm (Acacia auriculiformis A.Cunn ae Benth), Keo tai tượng (A.mangium Willd) gần giống lai tự nhiên hai loài (gọi tắt Keo lai - Acacia hybrid) trở thành loài đưa vào trồng rừng đại trà số loài cấu trồng Chương trình, Dự án trồng rừng nước ta, đặc biệt trồng rừng sản xuất cung cấp nguyên liệu thô cho ngành công nghiệp giấy, công nghiệp ván nhân tạo (ván dăm, ván ép, ván dán ) Keo lai Nhà khoa học Nhà kinh doanh đánh giá lồi có nhiều triển vọng việc tạo nên vùng nguyên liệu gỗ tập trung cho công nghiệp Kết nghiên cứu khảo nghiệm trồng rừng thử nghiệm địa phương nước bước đầu cho thấy khả sinh trưởng tăng trưởng Keo lai cao so với hai giống bố mẹ Keo tai tượng (Acacia mangium) Keo tràm (Acacia auriculiformis) Keo lai khẳng định loài ưu việt so với lồi Keo khác, là: Biên độ sinh thái rộng, thích ứng với nhiều điều kiện lập địa khác nhau, có khả chịu đựng khô hạn, tăng trưởng nhanh vùng đất nghèo dinh dưỡng Ngoài ra, Keo lai loài sinh trưởng nhanh, chu kỳ kinh doanh ngắn nên kinh doanh Keo lai nhanh chóng tạo nguồn nguyên liệu đem lại lợi ích cho đơn vị kinh doanh lâm nghiệp cho người sản xuất Để nâng cao suất, chất lượng hiệu của công tác trồng rừng kinh doanh rừng Keo lai việc nghiên cứu nắm bắt quy luật khách quan tồn đời sống lâm phần, nghiên cứu động thái cấu trúc sinh trưởng làm sở để dự doán sản lượng rừng Keo lai thời điểm điều tra, kinh doanh rừng khác cần thiết Ngoài xây dựng công cụ, bảng biểu chuyên dụng, phục vụ cho công tác điều tra, thống kê dự tính, dự báo sản lượng rừng, đề xuất biện pháp kỹ thuật lâm sinh kịp thời cho giai đoạn sinh trưởng Keo lai lâm phần Keo lai Qua tham khảo cơng trình nghiên cứu nước gần cho thấy: Quảng Trị địa phương có diện tích trồng rừng Keo lai tương đối lớn khu vực miền Trung, cơng trình nghiên cứu chuyên sâu lĩnh vực điều tra lâm sinh cho Keo lai tản mạn khiêm tốn Cơng trình nghiên cứu chun sâu lĩnh vực điều tra lâm sinh Nguyễn Trọng Bình (2005), tác giả nghiên cứu cấu trúc sinh trưởng tăng trưởng rừng Keo lai địa bàn rộng lớn rải khắp địa phương nước thu thành công bước đầu đáp ứng cho công tác kinh doanh rừng Keo lai nước ta Tuy nhiên, địa bàn nghiên cứu rộng, việc áp dụng cụ thể hóa cho địa phương cụ thể, đặc biệt khu vực Quảng trị nhiều vấn đề cần phải nghiên cứu Gần có Trần Xuân Dưỡng (2007) nghiên cứu cách tương đối hệ thống vấn đề cấu trúc, sản lượng phân chia cấp suất lập biểu cấp đất tạm thời cho rừng Keo lai tỉnh Quảng Trị Bên cạnh thành công bước đầu, chưa tìm thấy cơng trình nghiên cứu động thái cấu trúc rừng trồng Keo lai, nghiên cứu biến đổi theo thời gian trình sinh trưởng bình quân lâm phần, lập biểu cấp suất (hay cấp đất tạm thời) cho rừng Keo lai Cơng trình dừng lại việc nghiên cứu cấu trúc lâm phần trạng thái tĩnh, chưa thiết lập mơ hình dự báo sản lượng, chưa xét đến vận động biến đổi lâm phần theo giai đoạn tuổi khác cho đơn vị phân chia (từng cấp đất) Vì khó khăn việc dự tính dự báo sản lượng nói chung cho cấp suất nói riêng, tiến tới lập biểu q trình sinh trưởng cho Keo lai khu vực Để đáp ứng u cầu địi hỏi đó, chúng tơi tiến hành triển khai nghiên cứu đề tài: “Tìm hiểu số quy luật kết cấu lâm phần Keo lai (Acacia hybrid), phục vụ công tác điều tra rừng” Nghiên cứu đề tài nhằm góp phần giải yêu cầu thực tiễn sản xuất nay, thông qua nghiên cứu số phương pháp thực nghiệm thích hợp, phát quy luật cấu trúc bản, làm sở khoa học để dự tính, dự báo sản lượng rừng, đề xuất biện pháp kỹ thuật tác động hợp lý, phục vụ công tác kinh doanh nuôi dưỡng rừng keo lai nói riêng lồi trồng rừng đại trà khác địa phương khu vực miền Trung 1.2 MỤC ĐÍCH NGHIÊN CỨU Góp phần làm sáng tỏ thêm sở lý luận nghiên cứu cấu trúc, lâm phần, phục vụ công tác dự báo sản lượng Đề xuất biện pháp tác động hợp lý nhằm nâng cao suất, chất lượng hiệu công tác trồng rừng kinh doanh rừng Keo lai địa phương Bổ sung vào hệ thống quy luật kết cấu lâm phần rừng trồng lồi tuổi nói chung lâm phần Keo lai nói riêng PHẦN TỔNG QUAN VẤN ĐỀ NGHIÊN CỨU Nghiên cứu thiết lập mơ hình cấu trúc sinh trưởng rừng nhiều tác giả giới đề cập từ năm đầu kỷ XX Những nghiên cứu có xu hướng từ định tính đến định lượng quy luật tự nhiên, góp phần giải nhiều vấn đề kinh doanh rừng Để thiết lập mơ hình dự tính, dự báo sản lượng rừng việc nghiên cứu Sinh trưởng rừng lâm phần đánh giá nghiên cứu trọng tâm, có tính chất tảng để xây dựng mơ hình dự đốn sản lượng Định hướng nghiên cứu cấu trúc sản lượng rừng Nhà khoa học khái qt lại dạng mơ hình toán học từ đơn giản đến phức tạp nhằm định lượng quy luật tự nhiên, nhờ giải nhiều vấn đề kinh doanh rừng, đặc biệt lĩnh vực lập biểu chuyên dụng, phục vụ cho cơng tác điều tra, dự đốn sản lượng, đề xuất hệ thống biện pháp kinh doanh, nuôi dưỡng rừng cho đối tượng cụ thể Điểm qua số cơng trình nghiên cứu nước giới có liên quan tới nội dung nghiên cứu đề tài 2.1 TRÊN THẾ GIỚI 2.1.1 Nghiên cứu định lượng quy luật cấu trúc lâm phần Nghiên cứu quy luật cấu trúc nhằm xây dựng dạng cấu trúc phổ biến dạng tối ưu theo quan điểm kinh tế, nghĩa kiểu cấu trúc cho suất gỗ cao nhất, chất lượng gỗ phù hợp với nhu cầu sử dụng gỗ bảo vệ môi trường Dựa sở nghiên cứu quy luật cấu trúc với phát triển cách mạng khoa học kỹ thuật Nhà khoa học thu nhiều thành tựu khả quan 2.1.1.1 Nghiên cứu định lượng quy luật cấu trúc đường kính thân rừng Quy luật phân bố số theo cỡ kính (N/D) quy luật cấu trúc quan trọng hệ thống quy luật kết cấu lâm phần, nghiên cứu đầy đủ từ đầu kỷ XX Quy luật phân bố số theo cỡ đường kính biểu thị nhiều hình thức khác phân bố thực nghiệm N/D, phân bố số theo cỡ tự nhiên, biểu đồ đặc trưng thống kê Để nghiên cứu mô tả quy luật này, tác giả tìm phương trình tốn học nhiều dạng phân bố xác suất khác Các tác giả sau người xây dựng quy luật này: Veize (1880), Vimmenauer (1890, 1918), Schiffel (1898, 1899, 1902), Tretchiakov (1921, 1927, 1934, 1965), J Tuirin (1923, 1927, 1931, 1945), Moiseenko (1930, 1958), A noutchin (1931, 1936, 1954), Moiseev (1966, 1969, 1971), Prodan (1961, 1965) Các hàm toán học tác giả sử dụng để mô quy luật như: Hàm Hyperbol, họ đường cong Pearson, họ đường cong Poisson, hàm Chalier (kiểu A), Chalier (kiểu B) để xây dựng phân bố kinh nghiệm số theo đường kính N/D Ngồi nghiên cứu quy luật trạng thái tĩnh, tác giả sâu nghiên cứu biến đổi quy luật phân bố số theo thời gian theo biến số hay đại lượng có liên quan tới yếu tố thời gian mà điều tra rừng gọi động thái cấu trúc rừng + Hàm phân bố chuẩn Lôgarit Bliss, C, i ; Reinker, K, A (1964) xác lập tham số a, M, S phân bố chuẩn Lơgarit với đường kính bình qn theo dạng Lôgarit hai chiều: (2.1) (2.2) (2.3) + Hàm Weibull Clutter, J.L Allison, B.J (1973) dùng đường kính bình qn cộng, sai tiêu chuẩn đường kính đường kính nhỏ để tính tham số phân bố Weibull với giả thiết đại lượng có quan hệ với tuổi, mật độ lâm phần Quá trình biến đổi phân bố N/D theo tuổi, ngồi phụ thuộc vào sinh trưởng đường kính cịn chịu ảnh hưởng sâu sắc q trình tỉa thưa Từ Preussner đề nghị mơ hình tỉa thưa sở quan niệm biến đổi phân bố đường kính q trình xác định, nghĩa tổng hợp hai mơ hình: Mơ hình tỉa thưa mơ hình tăng trưởng đường kính tác giả sử dụng hàm: (2.4) Với: (2.5) (2.6) Trong đó: : Phần trăm số tỉa thưa theo cỡ kính i : Đường kính trung bình cỡ kính i : Đường kính nhỏ : Tham số : Các đại lượng biểu thị loại tỉa thưa : Tỷ lệ phần trăm chặt : Tuổi Hàm dùng xác định phân bố N/D phận tỉa thưa Để xác định phân bố này, cần phải biết phân bố N/D trước tỉa thưa, tuổi tỷ lệ chặt Số cịn lại sau tỉa thưa cỡ kính tính hiệu số số trước tỉa thưa số tỉa thưa Đối với mơ hình tăng trưởng, tác giả sử dụng hàm: (2.7) Với: Zi: Tăng trưởng đường kính cỡ kính i khoảng thời gian từ t đến t + t di: Đường kính trung bình cỡ kính i thời điểm t d: Đường kính trung bình cộng thời điểm t p.(t +t): Suất tăng trưởng đường kính a: Tham số phương trình (2.8) Do tăng trưởng số định chuyển dịch từ cỡ kính thấp đến cỡ kính cao Số xác định theo công thức hệ số chuyển cấp: (2.9) Hệ số phân thành hai phận f f2, f1 biểu thị phần nguyên f2 biểu thị phần thập phân Từ đó, số cỡ kính j thời điểm t chuyển lên cỡ kính i i + thời điểm t + t xác định sau: (2.10) (2.11) Trong đó: (2.12) Ngoài ra, hàm số thường sử dụng để xây dựng dãy phân bố kinh nghiệm số theo đường kính Nhà khoa học sử dụng như: + Hàm Beta Bennet F.A (1969) dùng phân bố Beta xác định đại lượng đường kính nhỏ (dm), đường kính lớn (dM) thơng qua phương trình tương quan kép với mật độ (N), tuổi (A) cấp đất (S) sau: (2.13) (2.14) Burkhart (1974) Strub (1972) tính tốn tham số d m, dM, phân bố Beta theo dạng phương trình: (2.15) (2.16) (2.17) (2.18) Với ho chiều cao tầng trội; A tuổi; N mật độ lâm phần + Hàm Gamma Roemisch, K (1975) nghiên cứu khả dùng hàm Gamma mô biến đổi phân bố đường kính rừng theo tuổi, xác lập quan hệ tham số Beta với tuổi, đường kính trung bình, chiều cao tầng trội khẳng định quan hệ tham số Beta với chiều cao tầng trội chặt chẽ Lembeke, Knapp Dittmar sử dụng phân bố Gamma với tham số thơng qua phương trình biểu thị mối tương quan với tuổi chiều cao tầng trội (2.19) (2.20) (2.21) Ngồi hàm tốn học đây: Hàm Mayer, hàm Hyperbol, hàm Poisson, hàm Charlier, hàm Logarit chuẩn, họ Pearson, đề cập nhiều nghiên cứu quy luật cấu trúc đường kính lâm phần cho đối tượng rừng khác nhau, đặc biệt rừng tự nhiên mà cơng trình khơng có điều kiện đề cập Một số tác giả khác: Suzuki (1971), Preussner.K (1974), Bock.W Diener.W (1972) lại nghiên cứu theo xu hướng khác với quan điểm đường kính rừng đại lượng ngẫu nhiên phụ thuộc vào thời gian coi trình biến đổi phân bố đường kính theo tuổi q trình ngẫu nhiên Q trình biểu thị tập hợp đại lượng ngẫu nhiên (Xt) với thời gian t lấy khoảng thời gian Nếu trị số đường kính thời điểm t phụ thuộc vào trị số thời điểm t - q trình Markov Nếu Xt = X có nghĩa q trình thời điểm t có dạng X Nếu tập hợp trạng thái xảy q trình Markov đếm chuỗi Markov, tức trị số t ứng với số tự nhiên Dùng hàm hàm khác để xây dựng dãy phân bố kinh nghiệm N/D phụ thuộc vào kinh nghiệm tác giả chất quy luật điều tra đo đạc Một dãy phân bố kinh nghiệm phù hợp cho dạng hàm số, phù hợp cho nhiều hàm số mức xác suất khác Từ nghiên cứu định lượng cấu trúc N/D đề cập cho thấy: - Các nghiên cứu phân bố số theo đường kính ứng dụng thường dựa vào dãy tần số lý thuyết - Các hàm toán học sử dụng để mô đa dạng phong phú - Xu hướng chung tìm hàm tốn học thích hợp, xác định tham số phân bố N/D hàm tương quan trực tiếp gián tuổi, thiết lập q trình ngẫu nhiên Ngồi ra, mô tả biến đổi phân bố N/D trình xác định sở quan niệm động thái phân bố N/D kết trình sinh trưởng trình tỉa thưa 2.1.1.2 Nghiên cứu quy luật quan hệ chiều cao với đường kính thân Quy luật tương quan H/D quy luật quan trọng hệ thống quy luật cấu trúc lâm phần Từ kết nghiên cứu nhiều tác giả cho thấy, chiều cao tương ứng với cỡ kính cho trước ln tăng theo tuổi, kết sinh trưởng tự nhiên lâm phần Trong cỡ kính xác định, cấp tuổi khác có thuộc cấp sinh trưởng khác Cấp sinh trưởng giảm tuổi lâm phần tăng lên dẫn đến tỷ lệ H/D tăng theo tuổi Từ đường cong quan hệ H/D thay đổi hình dạng ln dịch chuyển phía tuổi lâm phần tăng lên Vagui, A.B (1955) khẳng định "Đường cong chiều cao thay đổi dịch chuyển lên phía tuổi tăng lên" Tiurin.D.V (1927) phát hiện tượng ông xác lập đường cong chiều cao cấp tuổi khác Prodan.M (1944) nghiên cứu kiểu rừng “Plenterwal” kết luận đường cong chiều cao khơng bị thay đổi vị trí cỡ kính định Prodan.M (1965) lại phát độ dốc đường cong chiều cao có chiều hướng giảm dần tuổi tăng lên Curtis.R.O mô quan hệ chiều cao với đường kính tuổi theo dạng phương trình: (2.22) Tiếp theo nắn theo đường định kỳ năm tương ứng với định kỳ kiểm kê tài nguyên rừng Lĩnh Sam, tuổi định phương trình s l: (2.23) 10 Các phơng trình phơng trình phi tuyến cần thiết phải trải qua bớc tuyến tính hóa thông qua việc đặt ẩn số phụ để đa dạng phơng trình tuyến tính lớp Sử dụng số liệu điều tra đo cao lâm phần phơng pháp bình phơng bé để thiết lập phơng trình Kết lập phơng trình tơng quan H/D cho lâm phần đợc trình bày bảng Bng 6: Kt qu thăm dị dạng phương trình tốn mơ tả mối tương quan chiều cao với đường kính thân tiêu chuẩn Ơtc DPT pt lập   H=6.79+0.43D H=-6.69+12.06lgD   lgH=0.62+0.43lgD   H=6.88+0.43D 2 H=0.07+11.4lgD   lgH=0.65+0.41lgD   H=4.15+0.71D H=-5.91+17.38lgD   lgH=0.37+0.68lgD   H=4.71+0.56D H=-3.62+14.02lgD   lgH=0.41+0.60lgD   H=5.26+0.50D H=-3.14+13.49lgD   lgH=0.44+0.57lgD   H=5.15+0.54D 10 H=-3.11+13.88lgD   lgH=0.44+0.59lgD   H=5.15+0.54D 11 H=-3.38+13.60lgD   lgH=0.52+0.53lgD   H=6.27+0.50D 12 H=-1.48+12.86lgD   lgH=0.54+0.51lgD   H=6.10+0.50D 13 H=-1.38+12.65lgD R 0.85 0.84 0.84 0.90 0.87 0.88 0.79 0.80 0.81 0.78 0.80 0.81 0.83 0.87 0.87 0.70 0.73 0.75 0.91 0.92 0.92 0.67 0.69 0.70 0.83 0.85 Sy Ta Tb(Tr) n T05(K) 1.168 13.72 13.34 72 1.99 1.195 -0.64 12.91 72 1.99 0.043 15.80 12.85 72 1.99 0.907 18.05 16.62 70 2.00 0.996 0.08 14.76 70 2.00 0.034 21.51 15.38 70 2.00 1.302 6.00 11.76 84 1.99 1.286 -3.91 11.99 84 1.99 0.049 6.50 12.30 84 1.99 1.004 8.17 11.04 81 1.99 0.951 -2.98 12.04 81 1.99 0.040 8.02 12.34 81 1.99 0.909 11.81 13.48 81 1.99 0.816 -3.43 15.63 81 1.99 0.034 11.46 15.74 81 1.99 1.358 7.09 8.65 80 1.99 1.306 -2.01 9.35 80 1.99 0.051 7.20 10.16 80 1.99 0.790 16.93 18.70 77 1.99 0.749 -3.39 19.96 77 1.99 0.029 18.70 19.67 77 1.99 1.214 9.27 7.95 80 1.99 1.176 -0.96 8.52 80 1.99 0.046 8.98 8.64 80 1.99 0.988 13.78 12.83 76 1.99 0.948 -1.44 13.60 76 1.99 46     14     15     16     17     18     19     20     21     22     23     24     25     26   3 3 3 3 3 3 3 lgH=0.55+0.49lgD H=6.06+0.51D H=-1.69+13.03lgD lgH=0.55+0.50lgD H=5.75+0.53D H=-1.71+12.87lgD lgH=0.54+0.501gD H=6.00+0.27D H=-4.21+12.05lgD lgH=0.44+0.471gD H=5.80+0.27D H=-5.69+13.05lgD lgH=0.38+0.511gD H=6.80+0.22D H=-1.53+9.88lgD lgH=0.56+0.371gD H=7.81+0.17D H=1.10+7.89lgD lgH=0.66+0.30lgD H=7.58+0.19D H=-0.10+8.82lgD lgH=0.61+0.33lgD H=7.70+0.18D H=0.09+8.03lgD lgH=0.65+0.31lgD H=7.78+0.19D H=0.20+8.55lgD lgH=0.63+0.33lgD H=6.26+0.25D H=-3.45+11.01lgD lgH=0.5+0.43lgD H=7.03+0.21D H=-0.87+9.43lgD lgH=0.58+0.36lgD H=6.75+0.24D H=-2.29+10.67lgD lgH=0.54+0.40lgD H=7.54+0.18D H=0.41+8.37lgD lgH=0.63+0.32lgD 0.85 0.79 0.80 0.80 0.85 0.86 0.86 0.86 0.86 0.86 0.84 0.85 0.85 0.86 0.84 0.84 0.85 0.84 0.84 0.83 0.83 0.84 0.88 0.87 0.87 0.86 0.85 0.85 0.90 0.89 0.89 0.84 0.82 0.82 0.93 0.91 0.92 0.83 0.82 0.82 0.036 1.085 1.061 0.041 0.883 0.858 0.033 0.640 0.640 0.025 0.697 0.675 0.026 0.545 0.585 0.022 0.443 0.461 0.017 0.465 0.460 0.017 0.405 0.422 0.016 0.444 0.454 0.017 0.375 0.395 0.015 0.595 0.619 0.024 0.412 0.448 0.016 0.491 0.501 0.019 15.14 12.25 -1.49 12.63 13.59 -1.85 15.08 12.50 -3.16 8.51 10.42 -3.70 6.51 17.59 -1.28 12.74 24.42 1.15 18.40 19.95 -0.09 14.99 27.41 1.06 20.44 22.98 0.20 16.94 18.91 -2.99 12.60 16.89 -0.72 12.64 23.64 -2.58 17.05 20.94 0.39 15.43 13.94 11.26 11.65 11.72 13.67 14.21 14.46 11.90 11.90 11.86 10.84 11.35 11.57 11.97 10.83 11.02 11.38 10.79 10.88 10.40 10.55 10.69 13.16 12.47 12.56 11.80 11.44 11.54 14.77 13.87 13.95 10.94 10.34 10.36 17.31 15.68 16.36 10.46 10.15 10.12 76 77 77 77 76 76 76 51 51 51 53 53 53 52 52 52 52 52 52 51 51 51 53 53 53 52 52 52 53 53 53 53 53 53 52 52 52 52 52 52 1.99 1.99 1.99 1.99 1.99 1.99 1.99 2.01 2.01 2.01 2.01 2.01 2.01 2.01 2.01 2.01 2.01 2.01 2.01 2.01 2.01 2.01 2.01 2.01 2.01 2.01 2.01 2.01 2.01 2.01 2.01 2.01 2.01 2.01 2.01 2.01 2.01 2.01 2.01 2.01 47 Qua kết thu đợc bảng nhận thấy rằng: Mối quan hệ đờng kính chiều cao thân luôn tồn mức độ tơng quan chặt chẽ Dựa vào nguyên tắc chọn phơng trình thích hợp phơng trình có hệ số tơng quan cao đng thời có sai số bé nh mức độ tồn tham số hồi quy phơng trình lập sở số liệu thực nghiệm, đà định chọn phơng trình hay nói cách khác để mô tả quy luật tơng quan chiều cao với đờng kính thân cho lâm phần, tốt sử dụng dạng phơng trình Lôgarít chiều thích hợp LgH = a + blgD Bảng 7, tổng hợp kết chọn lựa phơng trình tốt mô tả mối quan hệ chiều cao với đờng kính cho lâm phần Bng 7: Tổng hợp kết lập phân tích hồi quy phương trình tương quan H/D theo dạng phương trình LgH =a+blgD tiêu chuẩn Ơtc pt lập R Sy Ta Tb(Tr) n T05(K) lgH=0.62+0.43lgD 0.84 0.043 15.80 12.85 72 1.99 lgH=0.65+0.41lgD 0.88 0.034 21.51 15.38 70 2.00 lgH=0.59+0.45lgD 0.87 0.037 16.77 14.49 67 2.00 lgH=0.66+0.40lgD 0.89 0.032 22.53 15.49 68 2.00 lgH=0.61+0.43lgD 0.90 0.032 21.45 16.97 69 2.00 lgH=0.48+0.56lgD 0.76 0.047 8.67 10.53 83 1.99 lgH=0.37+0.68lgD 0.81 0.049 6.50 12.30 84 1.99 lgH=0.41+0.60lgD 0.81 0.040 8.02 12.34 81 1.99 lgH=0.44+0.57lgD 0.87 0.034 11.46 15.74 81 1.99 10 lgH=0.44+0.59lgD 0.75 0.051 7.20 10.16 80 1.99 11 lgH=0.52+0.53lgD 0.92 0.029 18.70 19.67 77 1.99 12 lgH=0.54+0.51lgD 0.70 0.046 8.98 8.64 80 1.99 13 lgH=0.55+0.49lgD 0.85 0.036 15.14 13.94 76 1.99 17 lgH=0.38+0.511gD 0.85 0.026 6.51 11.57 53 2.01 18 lgH=0.56+0.371gD 0.84 0.022 12.74 11.02 52 2.01 19 lgH=0.66+0.30lgD 0.84 0.017 18.40 10.88 52 2.01 20 lgH=0.61+0.33lgD 0.84 0.017 14.99 10.69 51 2.01 21 lgH=0.65+0.31lgD 0.87 0.016 20.44 12.56 53 2.01 22 lgH=0.63+0.33lgD 0.85 0.017 16.94 11.54 52 2.01 23 lgH=0.5+0.43lgD 0.89 0.015 12.60 13.95 53 2.01 24 lgH=0.58+0.36lgD 0.82 0.024 12.64 10.36 53 2.01 48 25 lgH=0.54+0.40lgD 0.92 26 lgH=0.63+0.32lgD 0.82 KÕt qu¶ b¶ng cho thÊy: 0.016 17.05 0.019 15.43 16.36 10.12 52 52 2.01 2.01 C¸c tham sè håi quy tính toán đợc phơng trình có khác nhau, diện tích điều tra quy trình tính toán Hệ số tơng quan tính toán đợc phơng trình lớn (thấp 0,70 rơi vào lâm phần số 12 cao 0,92 rơi vào lâm phần 25) chứng tỏ chiều cao thân đờng kính ngang ngực có tơng quan đồng biến chặt chẽ Sai số phơng trình tính đợc (S, S% P%) nhỏ phạm vi sai số cho phép công tác nghiên cứu Vì vậy, phơng trình lập đợc có độ xác cao Hầu hết phơng trình lập đợc có trị tuyệt đối ta, tb tr lớn t05 tra bảng, chứng tỏ phơng trình lập đợc tồn mẫu mà lâm phần Nh việc sử dụng dạng phơng trình lôgarít chiều để mô tả quy luật tơng quan chiều cao với đờng kính thân rừng hoàn toàn hợp lý có sở khoa học Để xem xét khả thành lập phơng trình bình quân chung mô tả mối quan hệ chiều cao với đờng kính thân cây, tiến hành kiểm tra 26 phơng trình lập đợc bảng Kết kiểm tra nhất, cho bảng 8: Bng 8: Kt kiểm tra phương trình tương quan chiều cao với đường kính thân PT bi Sbi W=1/Sbi2 0.43 0.033 918.2736 0.41 0.026 1433.537 0.45 0.031 1044.245 17 0.51 0.044 512.7664 18 0.37 0.034 866.3856 19 0.3 0.028 1313.834 20 0.33 0.031 1016.228 Wbi*bi Wbi*bi2 394.8577 169.7888 587.7502 240.9776 469.9102 211.4596 261.5109 133.3705 320.5627 118.6082 394.1503 118.2451 335.3553 110.6673 49 21 0.31 0.025 1661.54 515.0773 22 0.33 0.028 1234.141 407.2666 23 0.43 0.031 1067.682 459.1032 24 0.36 0.035 814.4063 293.1863 25 0.4 0.024 1677.656 671.0626 26 0.32 0.032 994.0279 318.0889 Tổng cộng 24250.69 10237.93 Kết ti bng cho thấy: Khi bình phơng 159.6739 134.398 197.4144 105.5471 268.425 101.7885 4503.234 tính toán đợc ( n2=181.08) lớn nhiều so với bình phơng tra bảng với xác suất 95% (n2=37.6) Vì vậy, cha có sở khoa học để thành lập phơng trình bình quân chung, biểu diễn mối quan hệ chiều cao thân với đờng kính thân cho đối tợng nghiên cứu Để nghiên cứu mối quan hệ này, cần thiết lô rừng hay lâm phần nghiên cứu cần lập phơng trình tơng quan H/D theo dạng phơng trình LgH=a+blgD với tham số a, b đợc xác định từ tài liÖu thùc nghiÖm 5.2.3 Quy luật tương quan đường kớnh tỏn vi ng kớnh ngang ngc Để thấy đợc toàn cảnh tình hình sinh trởng lâm phần, bên cạnh việc nghiên cứu quy luật tơng quan chiều cao với đờng kính thân cây, tiến hành nghiên cứu mối quan hệ đờng kính tán với đờng kính ngang ngực Khác với quy luật tơng quan H/D, để mô tả quy luật tơng quan Dt/D13 qua khảo cứu công trình nghiên cứu nớc lĩnh vực hầu hết tác giả sử dụng dạng phơng trình quan hÖ: Dt = a + b.D13 KÕ thõa kinh nghiÖm đà tiến hành mô tả quy luật thông qua việc lập phân tích hồi quy sở nguồn số liệu thu thập đợc từ lâm phần khác Kết đợc tổng hợp bảng 9: 50 Bảng 9: Kết lập phân tích hồi quy phương trình tương quan đường kính tán với đường kính thân ƠTC 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 pt lập Dt1=1.02+0.11D13 Dt2=1.17+0.09D13 Dt3=1.32+0.09D13 Dt4=0.89+0.11D13 Dt5=1.16+0.10D13 Dt6=0.58+0.17D13 Dt7=-0.14+0.22D13 Dt8=0.53+0.16D13 Dt9=0.14+0.19D13 Dt10=0.64+0.16D13 Dt11=0.81+0.13D13 Dt12=1.21+0.08D13 Dt13=0.73+0.14D13 Dt14=0.65+0.14D13 Dt15=0.54+0.51D13 Dt16=2.46+0.15D13 Dt17=3.01+0.13D13 Dt18=3.12+0.13D13 Dt19=3.52+0.11D13 Dt20=3.86+0.09D13 Dt21=3.12+0.12D13 Dt22=3.01+0.13D13 Dt23=1.4+0.23D13 Dt24=2.77+0.14D13 Dt25=2.70+0.14D13 Dt26=3.08+0.19D13 R 0.84 0.66 0.71 0.80 0.75 0.66 0.83 0.70 0.84 0.63 0.82 0.52 0.83 0.81 0.86 0.79 0.75 0.77 0.73 0.68 0.74 0.74 0.94 0.79 0.80 0.73 Sy 0.299 0.434 0.391 0.341 0.373 0.250 0.358 0.370 0.345 0.480 0.294 0.297 0.269 0.274 0.033 0.471 0.436 0.446 0.419 0.376 0.472 0.469 0.266 0.476 0.440 0.477 Ta Tb(Tr) n T05(K) 8.09 13.13 72 1.99 6.44 7.25 70 2.00 7.62 8.11 67 2.00 5.95 10.82 68 2.00 7.39 9.20 69 2.00 2.32 7.82 83 1.99 -0.72 13.40 84 1.99 2.50 8.70 81 1.99 0.84 13.83 81 1.99 2.50 7.11 80 1.99 6.57 12.26 77 1.99 7.30 5.33 80 1.99 6.07 12.94 76 1.99 5.24 12.17 77 1.99 15.08 14.46 76 1.99 6.98 9.00 51 2.01 8.63 8.07 53 2.01 9.85 8.67 52 2.01 11.65 7.56 52 2.01 12.56 6.47 51 2.01 9.54 7.85 53 2.01 8.75 7.72 52 2.01 5.96 19.35 53 2.01 8.32 9.19 53 2.01 8.86 9.39 52 2.01 8.81 7.46 52 2.01 Qua kết bảng nhận thấy rằng: Đờng kính tán có mối quan hệ với đờng kính thân mức độ chặt chẽ với hệ số tơng quan r tính đợc biến động từ 0.52 (phơng trình lõm phn 12) đến 0,94 (phơng trình lõm phn 24), chứng tỏ loại đờng kính tán có mối quan hệ tơng đối chặt chẽ với đờng kính ngang ngực Mặt khác, sai số phơng trình nhỏ Vì cho phép sử dụng phơng trình ®Ĩ tÝnh to¸n tỉng diƯn tÝch 51 t¸n rõng sÏ cho kết với độ xác đảm bảo yêu cầu công tác điều tra rừng Trị tuyệt đối tb tr lớn nhiều so với t05 tra bảng điều chứng tỏ phơng trình tơng quan lập đợc tồn cao mẫu quan sát mà tổng thể Từ kết nghiên cứu đề nghị sử dụng phơng trình để mô tả quy luật tơng quan Dt/D13 cho rừng Keo lai hợp lý cho độ xác cao Để xem xét khả thành lập phơng trình bình quân chung mô tả mối quan hệ ng kớnh tỏn cõy với đờng kính thân cây, tiến hành kiểm tra 26 phơng trình lập đợc bảng trªn thơng qua tiêu chuẩn phù hợp Pearson người Anh xut Kết kiểm tra nhất, cho bảng 10: Bng 10: Kt qu kim tra thun phương trình tương quan đường kính tán với đường kính thân PT bi Sbi W=1/Sbi2 0.11 0.008 14722.44 0.09 0.012 6649.7193 0.09 0.011 7699.6076 0.11 0.010 9649.8717 0.1 0.011 8897.8103 0.17 0.453 4.8628109 0.22 0.224 19.963874 0.16 0.019 2870.0983 0.19 0.014 5068.0306 10 0.16 0.022 1981.6481 11 0.13 0.011 8668.7297 12 0.08 0.015 4196.2928 13 0.14 0.011 8858.7675 Wbi*bi Wbi*bi2 1619.468 178.1415 598.4747 53.86273 692.9647 62.36682 1061.486 116.7634 889.781 88.9781 0.826678 0.140535 4.392052 0.966252 459.2157 73.47452 962.9258 182.9559 317.0637 50.73019 1126.935 146.5015 335.7034 26.85627 1240.227 173.6318 52 14 0.14 0.012 7558.7397 15 0.51 0.035 819.04498 16 0.15 0.017 3653.7751 17 0.13 0.016 4114.2313 18 0.13 0.015 4351.7113 19 0.11 0.014 4804.8005 20 0.09 0.015 4749.3441 21 0.12 0.016 4005.0216 22 0.13 0.017 3505.3423 23 0.23 0.012 7012.7542 24 0.14 0.016 4115.4465 25 0.14 0.015 4667.9779 26 0.19 0.017 3395.5615 Tổng cộng 136041.59 KÕt qu¶ bảng 10 cho thÊy: 1058.224 417.7129 548.0663 534.8501 565.7225 528.5281 427.441 480.6026 455.6945 1612.933 576.1625 653.5169 645.1567 17814.08 148.1513 213.0336 82.20994 69.53051 73.54392 58.13809 38.46969 57.67231 59.24029 370.9747 80.66275 91.49237 122.5798 2621.069 Khi bình phơng tính toán đợc ( n2=288.4) lớn nhiều so với bình phơng tra bảng với xác suất 95% (n2=37.7) Vì vậy, cha có sở khoa học để thành lập phơng trình bình quân chung, biểu diễn mối quan hệ ng kớnh tỏn với đờng kính thân cho đối tợng nghiên cứu Để nghiên cứu mối quan hệ này, cần thiết lô rừng hay lâm phần nghiên cứu cần lập phơng trình tơng quan Dt/D1.3 theo dạng phơng trình Dt=a+bD13 với tham số a, b đợc xác định từ tài liệu thực nghiệm 5.3 MỘT SỐ ỨNG DỤNG KẾT QUẢ NGHIÊN CỨU PHỤC VỤ CƠNG TÁC ĐIỀU TRA RỪNG 5.3.1 ỈÏng dủng quy lût cáúu trục âỉåìng kênh Lám pháưn Quy lût cáúu trục lám pháưn nọi chung, cáúu trục âỉåìng kênh (hay quy lût phán bäú säú cáy theo âỉåìng kênh) nọi riãng, l cå såí khoa hc cho cạc phỉång phạp thäúng kã, dỉû 53 âoạn trỉỵ lỉåüng, sn lỉåüng v toạn cạc chè tiãu kü thût kinh doanh, âiãưu chãú rỉìng Tỉì kãút qu nghiãn cỉïu quy lût ny, cho phẹp xạc âënh cạc nhán täú âiãưu tra cå bn ca lám pháưn åí thåìi âiãøm hiãûn tải nhỉ: Máût âäü (N), täøng tiãút diãûn ngang (G), trỉỵ lỉåüng (M), caùc loaỷi õổồỡng kờnh bỗnh quỏn nhổ : D, D g ,Dgo Nghiãn cæïu quy luáût phán bäú N/D cn cọ nghéa quan trng kinh doanh, ni dỉåỵng rỉìng Tỉì kãút qu nghiãn cỉïu phán bäú N/D, biãút âỉåüc trảng thại lám pháưn âiãưu tra, giụp cho viƯc âiãưu chènh máût âäü theo hỉåïng cọ låüi Tọm lải: Quy lût phán bäú säú cáy theo âỉåìng kênh laỡ mọỹt nhổợng luỏỷt cỏỳu truùc coù nhỏỳt cuớa lỏm phỏửn Vỗ vỏỷy, noù laỡ nọỹi dung chờnh âiãưu tra lám pháưn 5.3.2 ỈÏng dủng quy lût tỉång quan giỉỵa chiãưu cao v âỉåìng kênh thán cáy Âáy cng l mäüt quy lût cáúu trục cå bn vaỡ quan troỹng cuớa Lỏm phỏửn Vióỷc nghión cổùu tỗm hiãøu v nàõm vỉỵng quy lût ny cọ nghéa låïn cäng tạc âiãưu tra, kinh doanh v ni dỉåỵng rỉìng Trỉåïc hãút chiãưu cao cng l mét nhỉỵng nhán täú cáúu thnh thãø têch thán cáy v trỉỵ lỉåüng lám pháưn, khäng thãø thiãúu âỉåüc sỉí dủng cạc biãøu, phủc vủ cäng tạc âiãưu tra v kinh doanh rỉìng Thäng qua quy lût ny, kãút håüp våïi cạc quy lût N/D cho phẹp xạc âënh cạc nhán täú âiãưu tra cå bn ca lám pháưn, phủc vủ cäng tạc nghiãn cỉïu tiãúp theo Tỉång quan H/D cn l cå såí âãø phán chia cáúp nàng sút ca rỉìng 54 Nh×n chung: Trõ chiỊu cao b×nh quân cộng chiều cao bình quân Lorêy, loại chiều cao bình quân đợc tính thông qua đờng kính bình quân tơng ứng phơng trình tơng quan H/D lập cho lô rừng, cách thay đờng kính bình quân vào phơng trình tơng quan H/D Tng hp ứng dụng kết nghiên cứu quy luật N/D, H/D Dt/D13 tính tốn nhân tố điều tra lâm phần trình bày bảng 11 Bảng 11: Ứng dụng tổng hợp quy luật cấu trúc lâm phần, xác định số nhân tố điều tra lâm phần ÔTC N/ha 1440 1400 1340 1360 1380 1660 1680 1620 1620 10 1600 11 1540 12 1600 14 1540 15 1520 16 1020 17 1060 18 1040 19 1040 20 1020 21 1060 22 1040 23 1060 24 1060 25 1040 26 1040 G/ha 29.3457 26.4098 26.7149 25.2030 26.2097 18.6954 19.2876 18.0546 20.2075 18.1438 16.8724 16.3650 16.0692 16.2411 39.6030 45.2813 38.8227 38.5068 38.1302 38.4491 36.9627 34.7006 41.1590 38.2007 36.9998 Dg Hg Vcbq 16.11 13.8 0.1460 15.50 13.7 0.1348 15.93 13.5 0.1402 15.36 13.6 0.1314 15.55 13.3 0.1309 11.98 12.1 0.0710 12.09 12.8 0.0762 11.91 11.4 0.0659 12.60 11.7 0.0757 12.02 11.9 0.0704 11.81 12.3 0.0698 11.41 12.0 0.0638 11.53 12.0 0.0654 11.66 11.8 0.0658 22.23 16.3 0.3301 23.32 12.0 0.2656 21.80 11.4 0.2205 21.71 11.5 0.2216 21.82 11.3 0.2190 21.49 11.6 0.2181 21.27 11.7 0.2162 20.42 11.6 0.1969 22.24 11.6 0.2345 21.63 11.9 0.2265 21.28 11.3 0.2100 M/ha St/ha 210.1834 8669.01 188.7130 7107.87 187.8463 7840.37 178.6658 7716.34 180.6908 7944.50 117.9211 8129.90 128.0452 8396.27 106.7105 7846.48 122.6863 8840.37 112.6711 7894.55 107.5255 6566.41 102.1390 7078.33 100.6618 6365.70 100.0139 6388.24 336.7143 10491.61 281.5155 11622.25 229.2702 10439.04 230.4480 9715.93 111.7023 10600.98 231.1558 10742.87 224.8886 10381.38 208.7602 9952.03 248.5540 11327.87 235.5507 10643.04 218.3707 10383.68 55 5.3.3 Xác định loại tiêu chuẩn theo lý thuyết tiêu chuẩn thực địa Trong thực tế lý thuyết thực tế nh Vì vậy, để xác định tiêu chuẩn thực tế, ô mẫu nghiên cứu, điều tra rừng cho phép sai số xác định kích thớc đờng kính chiều cao vút 10% Tuỳ theo mục đích nghiên cứu khác mà lựa chọn bình quân khác Trên sở bình quân, tiến hành chặt hạ giải tích, cho phép tính toán đợc loại lợng tăng trởng, nghiên cứu quy luật sinh trởng cho đại lợng đờng kính, chiều cao, thể tích nữa, cho phép xác định tốc độ tăng trởng tối đa xác định chu kỳ kinh doanh cho loài Tuy nhiên, khuôn khổ báo cáo có hạn mục tiêu nghiên cứu nh điều kiện không cho phép, nghiên cứu này, đa số ứng dụng quan trọng quy luật cấy trúc, đại diện cho hệ thống quy luật cấu trúc lâm phần để từ định hớng cho nghiên cứu tiÕp theo 56 PHẦN KÃÚT LÛN CHUNG, TÄƯN TẢI V KIÃÚN NGHË 6.1 KÃÚT LÛN CHUNG Xút phạt tỉì thỉûc tiãùn khạch quan ca sn xút lám nghiãûp, chụng täi nghiãn cỉïu âãư ti: “Tìm hiểu số quy luật kết cấu lâm phần Keo lai phục vụ cơng tác điều tra rừng” Våïi nhỉỵng kãút qu thu âỉåüc, chụng täi xin rụt mäüt säú kãút luỏỷn sau õỏy: 6.1.1 Các quy luật cấu trúc lâm phần nhìn chung tuân thủ theo quy luật chung lâm phần loài, tuổi nớc ta, là: * Phân bố số theo cỡ đờng kính có dạng đờng cong lệch trái, lệch trái đối xứng Tuy nhiên cá biệt có lâm phần ®ang chun sang giai ®o¹n lƯch hay cã xu híng lệch phải, Độ lệch phân bố nhìn chung không tuân theo quy luật * Tơng quan chiều cao với đờng kính thân cây, luôn tồn mức độ chặt chẽ mức ý nghĩa 95% Qua nghiên cứu, đà xác định đợc dạng phơng trình toán học thích hợp để mô tả mối quan hệ này, chiều cao với ®êng kÝnh, thĨ: Sư dơng quan hƯ l«garit chiều: LgH = a+bLgD thích hợp cho độ xác cao Qua phân tích hồi quy kiểm tra tính thích ứng câc kết nghiên cứu đề nghị sử dụng vào công tác điều tra nghiên cứu khoa học địa bàn nghiên cøu * Đỉåìng kênh tạn cáy våïi âỉåìng kênh thán cáy tải vë trê âäü cao 1,3m ln ln täưn taûi mäúi quan hãû tuyãún mäüt låïp: Dt=a+bD13 6.1.2 Ba quy lût chung ca lám pháưn âãưu chëu sỉû chi phäúi täøng håüp ca bäún úu täú âọ l: Âiãưu kiãûn láûp 57 âëa, úu täú âëa phỉång, âàûc âiãøm củ thãø ca lám pháưn v máût âäü hiãûn cn 6.1.3 Trãn cå såí cạc quy lût ny â đề xuất đợc hớng ứng dụng kết vào công tác điều tra rừng để tính toán xaùc õởnh âỉåüc cạc nhán täú âiãưu tra cå bn cho cạc lỏm phỏửn, laỡm cồ sồớ cho nghiên cứu chuyên sâu nhằm phục vụ cho công tác trồng rõng cịng nh sư dơng rõng 6.2 TÄƯN TẢI Tỉì nhỉỵng kãút qu m âãư ti âảt âỉåüc, chụng täi cng nháûn tháúy mäüt säú màût täưn tải sau âáy: 6.2.1 Trỉåïc hãút ngưn säú liãûu nghiãn cỉïu ca âãư ti cn hản chãú c vãư khu vỉûc nghiãn cỉïu cng âäü tøi nghiãn cỉïu 6.2.2 Viãûc nghiãn cỉïu xỏy dổỷng mọ hỗnh hỏửu nhổ chuớ yóỳu õi vaỡo hóỷ quaớ cuọỳi cuỡng cuớa quaù trỗnh sinh trổồớng vaỡ tàng trỉåíng cáy rỉìng, chỉa cọ âiãưu kiãûn nghiãn cỉïu sỉû nh hỉåíng ca cạc úu täú cnh 6.2.3 Âiãưu kiãûn ỉïng dủng kãút qu nghiãn cỉïu vo thỉûc tiãùn cn hản chãú Thỉûc cháút åí âáy chè måïi l nhỉỵng thnh cäng ban âáưu, cho nãn viãûc kiãøm tra tờnh thờch ổùng cuớa mọ hỗnh lỏỷp coỡn cọ nhỉỵng màût täưn tải nháút âënh 6.3 KIÃÚN NGHË Tỉì nhỉỵng kãút qu v täưn tải nãu trãn, chụng täi xin coï mäüt säú kiãún nghë sau âáy, âãø âãư ti hon chènh hån 6.3.1 Cáưn tiãúp tủc nghiãn cỉïu hon thiãûn âãư ti åí mỉïc âäü sáu v räüng hån c vãư phảm vi âëa bn nghiãn cỉïu cng âäü tøi nghiãn cỉïu 6.3.2 Ngoi nhỉỵng kãút qu mang cháút hãû qu ca cäng tạc träưng rỉìng l kêch thỉåïc cáy rỉìng tải thåìi âiãøm âiãưu tra cáưn nghiãn cỉïu bäø sung thãm cạc úu täú 58 cọ liãn quan khạc biãøu thë cho úu täú mäi trỉåìng nhỉ: nh hỉåíng ca cáúp âáút hay cáúp nàng sút ca lám pháưn, úu täú láûp âëa Nhỉỵng váún âãư ny âi hi phi cọ sỉû gii têch cáy rỉìng, phán têch v âạnh giạ theo quan âiãøm âäüng tỉïc l nghiãn cỉïu sỉû biãún âäøi ca lâm phần theo thồỡi gian mồùi õaớm baớo õuùng baớn chỏỳt cuớa rổỡng, vỗ cỏy rổỡng cuợng nhổ rổỡng cỏy l nhỉỵng âån vë sinh váût hc hon chènh ln ln åí trảng thại váûn âäüng v phạt triãøn 59 TI LIÃÛU THAM KHO [1]- Lê Đình Khả ( 1999), Nghiên cứu sử dụng giống lai tự nhiên Keo tai tượng Keo tràm Việt Nam, Nhà xuất Nông nghiệp [2]- Lưu Bá Thịnh, Phạm Văn Tuấn ( 1998), Báo cáo khoa học lâm nghiệp khảo nghiệm dịng vơ tính Keo lai Đơng Nam Bộ, Hội nghị tỉnh Đông Nam Bộ [3]- Lưu Bá Thịnh (1999), Báo cáo khoa học kết khảo nghiệm dịng vơ tính Keo lai tự nhiên tuyển chọn Đông Nam Bộ, Trung tâm KHSXLN Đông Nam Bộ [4]- Nguyễn Hải Tuất (1982), Thống kê toán học lâm nghiệp, Nhà xuất Nông nghiệp, Hà Nội [5]- Nguyễn Hải Tuất, Ngô Kim Khôi ( 1996), Xử lý thống kê kết nghiên cứu thực nghiệm Nơng nghiệp máy vi tính, nhà xuất nơng nghiệp, Hà Nội 60

Ngày đăng: 29/03/2023, 16:04

Xem thêm:

w