Luận án tiến sĩ lâm nghiệp nghiên cứu một số biện pháp kỹ thuật trồng rừng thâm canh keo lá tràm cung cấp gỗ lớn ở vùng đông bắc bộ

BỘ GIÁO DỤC VÀ ĐÀO TẠO BỘ NÔNG NGHIỆP VÀ PTNT VIỆN KHOA HỌC LÂM NGHIỆP VIỆT NAM ==================== PH M NH SÂM NGHIÊN CỨU MỘT SỐ BIỆN PHÁP KỸ THUẬT TRỒNG RỪNG THÂM CANH KEO LÁ TRÀM CUNG CẤP GỖ LỚN Ở VÙNG ÔNG BẮC BỘ LUẬN ÁN TIẾN SĨ LÂM NGHIỆP HÀ NỘI - 2022 BỘ GIÁO DỤC VÀ ĐÀO TẠO BỘ NÔNG NGHIỆP VÀ PTNT VIỆN KHOA HỌC LÂM NGHIỆP VIỆT NAM ==================== PH M NH SÂM NGHIÊN CỨU MỘT SỐ BIỆN PHÁP KỸ THUẬT TRỒNG RỪNG THÂM CANH KEO LÁ TRÀM CUNG CẤP GỖ LỚN Ở VÙNG ÔNG BẮC BỘ N n o t o: Lâm sinh M n n : 9.62.02.05 LUẬN ÁN TIẾN SĨ LÂM NGHIỆP NG ỜI H ỚNG D N KHOA HỌC: PGS TS NGU ỄN HU S N HÀ NỘI - 2022 i LỜI CAM OAN Tôi xin cam đoan cơng trình nghiên cứu khoa học thân tơi, cơng trình đƣợc hồn thành theo chƣơng trình đào tạo tiến sĩ kho Viện Khoa học Lâm nghiệp Việt Nam, thời gian từ năm 16 đến 22 C c số liệu kết nghiên cứu trình bày luận n trung thực Luận n kế thừa phần số liệu đề tài nghiên cứu khoa học công nghệ cấp Bộ “ ” PGS.TS Nguyễn Huy Sơn làm chủ nhiệm, cộng t c viên chính, trực tiếp tham gia thiết kế thí nghiệm, thu thập số liệu rừng trồng Keo tràm có, xử lý số liệu, viết báo cáo nội dung nghiên cứu vùng Đông Bắc Bộ C c số liệu kế thừa đề tài đƣợc chủ nhiệm đề tài thành viên tham gia cho phép sử dụng công bố luận n Hà Nội, ngày 15 tháng năm 2022 N iên cứu sin P m n S m ii LỜI CẢM N Lời đầu tiên, t c giả xin bày tỏ lòng biết ơn sâu sắc đến ngƣời hƣớng dẫn khoa học PGS.TS Nguyễn Huy Sơn dành nhiều thời gian công sức cho việc hƣớng dẫn giúp đỡ nghiên cứu sinh hoàn thành luận n T c giả xin chân thành cảm ơn Ban Gi m đốc, Ban Khoa học, Đào tạo Hợp t c quốc tế - Viện Khoa học Lâm nghiệp Việt Nam; Ban lãnh đạo - Viện Nghiên cứu Lâm sinh quan tâm, đạo tạo điều kiện thuận lợi qu trình học tập, nghiên cứu hoàn thiện luận n T c giả c ng xin trân trọng cảm ơn nhóm thực đề tài Ban lãnh đạo Trƣờng Cao đẳng Nông lâm Đông Bắc h trợ t c giả suốt qu trình thực thí nghiệm thu thập số liệu trƣờng T c giả xin gửi lời cảm ơn sâu sắc tới c c thầy phản biện thành viên c c hội đồng cấp Bộ mơn cấp Cơ sở góp ý cho Luận n hoàn thiện T c giả c ng xin chân thành cảm ơn gia đình bạn bè đồng nghiệp h trợ động viên để t c giả hồn thành luận n Luận n khơng tr nh khỏi thiếu sót, t c giả mong nhận đƣợc ý kiến đóng góp c c nhà khoa học đồng nghiệp Hà Nội, ngày 15 tháng năm 2022 N iên cứu sin P m n S m iii MỤC LỤC LỜI CAM ĐOAN i LỜI CẢM ƠN ii MỤC LỤC iii DANH MỤC CÁC KÝ HIỆU VÀ TỪ VIẾT TẮT vi DANH MỤC BẢNG vii DANH MỤC HÌNH ix PHẦN MỞ ĐẦU Sự cần thiết luận n Mục tiêu nghiên cứu Ý nghĩa luận n Những đóng góp Đối tƣợng, địa điểm giới hạn nghiên cứu Cấu trúc luận n Chƣơng TỔNG QUAN VẤN ĐỀ NGHIÊN CỨU 1.1 Kh i qu t chung vấn đề nghiên cứu 1.1.1 Đặc điểm phân bố sinh th i loài Keo l tràm Lịch sử gây trồng Keo l tràm Việt Nam Một số kh i niệm kỹ thuật trồng rừng thâm canh 1.1.4 Quy định rừng trồng g lớn 1.2 Trên giới Nghiên cứu chọn giống Keo l tràm Nghiên cứu trồng rừng thâm canh Keo l tràm 11 Trong nƣớc 15 Nghiên cứu chọn giống Keo l tràm 15 1.3.2 Nghiên cứu trồng rừng thâm canh Keo l tràm 20 Thảo luận chung 26 Chƣơng NỘI DUNG VÀ PHƢƠNG PHÁP NGHIÊN CỨU 29 Nội dung nghiên cứu 29 iv Đ nh gi thực trạng số mơ hình rừng trồng Keo l tràm có tiềm cung cấp g lớn vùng Đơng Bắc Bộ số tỉnh lân cận 29 Nghiên cứu x c định số giống Keo l tràm thích hợp trồng rừng vùng Đông Bắc Bộ 29 Nghiên cứu ảnh hƣởng số biện ph p kỹ thuật trồng rừng thâm canh đến sinh trƣởng suất g Keo l tràm vùng Đông Bắc Bộ 29 2.1.4 Đề xuất số biện ph p kỹ thuật trồng rừng thâm canh g lớn Keo l tràm vùng Đông Bắc Bộ 29 Phƣơng ph p nghiên cứu 30 Quan điểm c ch tiếp cận 30 Vật liệu địa điểm nghiên cứu 32 Phƣơng ph p nghiên cứu 32 Đặc điểm điều kiện tự nhiên vùng Đông Bắc Bộ 42 Chƣơng KẾT QUẢ NGHIÊN CỨU VÀ THẢO LUẬN 44 Đ nh gi thực trạng số mơ hình rừng trồng Keo l tràm có tiềm cung cấp g lớn vùng Đông Bắc Bộ số tỉnh lân cận 44 Quy mơ số mơ hình điều tra rừng trồng Keo l tràm 44 Một số đặc điểm lập địa c c mơ hình Keo l tràm 45 3.1.3 Một số biện ph p kỹ thuật trồng rừng g lớn Keo l tràm 50 .4 Sinh trƣởng suất g số mơ hình rừng trồng Keo l tràm .55 Nghiên cứu x c định số giống Keo l tràm thích hợp trồng rừng vùng Đông Bắc Bộ 59 Đ nh gi sinh trƣởng số dịng vơ tính Keo l tràm 59 3.2.2 Đ nh gi suất g số dịng vơ tính Keo l tràm 65 3.2.3 Đ nh gi chất lƣợng thân số dịng vơ tính Keo l tràm 67 3.3 Nghiên cứu ảnh hƣởng số biện ph p kỹ thuật trồng rừng thâm canh đến sinh trƣởng suất g Keo l tràm vùng Đông Bắc Bộ 69 3.3 Ảnh hƣởng biện ph p xử lý thực bì đến sinh trƣởng suất g rừng trồng Keo l tràm 69 v 3.3 Ảnh hƣởng biện ph p làm đất đến sinh trƣởng suất g Keo l tràm 74 3.3.3 Ảnh hƣởng mật độ trồng đến sinh trƣởng suất g rừng trồng Keo l tràm 78 3.3.4 Ảnh hƣởng phân bón đến sinh trƣởng suất rừng trồng Keo l tràm 82 3.3.5 Ảnh hƣởng biện ph p chăm sóc đến suất chất lƣợng g Keo l tràm 90 3.3.6 Đặc điểm cấu trúc rừng trồng Keo l tràm giai đoạn năm tuổi giải ph p tỉa thƣa nuôi dƣỡng rừng g lớn 94 3.4 Đề xuất số biện ph p kỹ thuật trồng rừng thâm canh g lớn Keo l tràm vùng Đông Bắc Bộ 99 Chƣơng KẾT LUẬN – TỒN TẠI – KHUYẾN NGH 103 Kết luận 103 Tồn 104 4.3 Khuyến nghị 105 DANH MỤC CÁC CƠNG TRÌNH KHOA HỌC ĐÃ CÔNG BỐ LIÊN QUAN ĐẾN LUẬN ÁN 106 TÀI LIỆU THAM KHẢO 107 Tiếng Việt 107 Tiếng nƣớc 112 PHỤ LỤC 118 vi DANH MỤC CÁC KÝ HIỆU VÀ TỪ VIẾT TẮT Từ vi t tắt N ầ CT Cơng thức CTTN Cơng thức thí nghiệm Dt Đƣờng kính t n D1.3 Đƣờng kính ngang ngực ĐBB Đông Bắc Bộ F.sig X c suất F (Fisher) tính tốn G Tiết diện ngang Hdc Chiều cao dƣới cành Hvn Chiều cao vút KLT Keo tràm M Trữ lƣợng Nht Mật độ NN&PTNT Nông nghiệp Phát triển nông thôn Ntr Mật độ trồng ban đầu OTC Ô tiêu chuẩn SD Hệ số biến động đƣờng kính ngang ngực SH Hệ số biến động chiều cao vút TB Trung bình (gi trị trung bình c c tiêu đo đếm) TBKN Trung bình khảo nghiệm TBKT Tiến kỹ thuật TCVN Tiêu chuẩn Việt Nam VLHCSKT Vật liệu hữu sau khai th c TLS Tỷ lệ sống VS Phân hữu vi sinh D Tăng trƣởng bình quân hàng năm đƣờng kính H Tăng trƣởng bình qn hàng năm chiều cao vii DANH MỤC BẢNG Bảng Sinh trƣởng tăng trƣởng bình quân Keo l tràm 12 Bảng Kết khảo nghiệm chọn xuất xứ số vùng sinh th i 16 Bảng C c giống Keo l tràm đƣợc công nhận giống TBKT 18 Bảng Thông tin c c giống Keo l tràm đƣợc chọn để khảo nghiệm 32 Bảng Tổng hợp số lƣợng OTC điều tra 34 Bảng Tổng hợp c c mơ hình điều tra rừng trồng Keo l tràm vùng Đông Bắc Bộ số tỉnh lân cận 44 Bảng Đặc điểm lập địa nơi điều tra số mơ hình trồng Keo l tràm vùng Đông Bắc Bộ số tỉnh lân cận 47 Bảng 3.3 Đặc điểm đất số mơ hình trồng Keo l tràm vùng Đơng Bắc Bộ số tỉnh lân cận 49 Bảng 3.4 Tổng hợp lịch sử rừng trồng số biện ph p kỹ thuật p dụng trồng rừng Keo l tràm vùng Đông Bắc Bộ c c tỉnh lân cận 53 Bảng 3.5 Đặc điểm sinh trƣởng số mơ hình Keo l tràm vùng Đông Bắc Bộ số tỉnh lân cận 57 Bảng 3.6 Sinh trƣởng số dòng vơ tính Keo l tràm năm tuổi ng Bí, Quảng Ninh 63 Bảng 3.7 Năng suất g số dịng vơ tính Keo l tràm năm tuổi ng Bí, Quảng Ninh 67 Bảng 3.8 Chỉ tiêu chất lƣợng thân số dịng Keo l tràm năm tuổi ng Bí, Quảng Ninh 68 Bảng 3.9 Ảnh hƣởng biện ph p xử lý thực bì đến sinh trƣởng Keo l tràm năm tuổi ng Bí, Quảng Ninh 73 Bảng Ảnh hƣởng biện ph p làm đất đến sinh trƣởng Keo l tràm năm tuổi ng Bí, Quảng Ninh 76 viii Bảng Ảnh hƣởng mật độ đến tỷ lệ sống sinh trƣởng Keo l tràm năm tuổi ng Bí, Quảng Ninh 80 Bảng Ảnh hƣởng phân bón lót đến sinh trƣởng năm tuổi ng Bí, Quảng Ninh 83 Bảng 3 Kết phân tích N, P, K mẫu l Keo l tràm năm tuổi ng Bí, Quảng Ninh 84 Bảng Ảnh hƣởng phân bón thúc đến tỷ lệ sống sinh trƣởng Keo l tràm năm tuổi ng Bí, Quảng Ninh 87 Bảng Ảnh hƣởng biện ph p chăm sóc đến sinh trƣởng Keo l tràm năm tuổi ng Bí, Quảng Ninh 93 Bảng C c đặc trƣng cấu trúc rừng Keo l tràm năm tuổi, mật độ 1.660 cây/ha ng Bí, Quảng Ninh 94 Bảng C c đặc trƣng cấu trúc rừng Keo l tràm năm tuổi, mật độ ng Bí, Quảng Ninh 97 23 Cơng thức thí nghiệm (CTTN) TLS D1,3 Sd Hvn Sh Hdc Shdc Dt Sdt (%) (cm) (%) (m) (%) (m) (%) (m) (%) 0,6kg P2O5 + 0,1kg K2O CT3.2: 0,2kg NPK 88,8 6,1 24,7 6,2 19,5 2,7 21,0 85,5 6,2 26,2 6,7 23,8 2,6 18,8 86,1 5,8 26,4 6,8 23,4 2,5 25,4 86,0 5,5 19,3 5,9 23,5 2,4 26,9 CT4.1: 0,8kg P2O5 + 0,1kg K2O CT4.2: 0,2kg NPK CT5 Khơng bón (ĐC) Sig F 0,012 0,026 0,045 0,038 ΔM M (m3/ha) (m3/ /năm) 24 Cơng thức thí nghiệm (CTTN) ă TLS D1,3 Sd Hvn Sh Hdc Shdc Dt Sdt (%) (cm) (%) (m) (%) (m) (%) (m) (%) ΔM M (m3/ha) (m3/ /năm) ổi (6/2016-6/2019) CT1.1: 0,2kg P2O5 + 0,1kg 88,8 10,0 10,7 12,3 5,3 6,4 9,2 2,7 18,2 54,1 18,0 85,1 9,3 15,2 11,9 4,7 6,4 9,9 2,9 20,4 45,6 15,2 86,1 10,0 11,5 12,6 3,1 6,6 9,8 2,9 21,1 55,7 18,5 85,1 9,8 18,8 12,1 6,3 6,3 10,1 3,0 22,1 51,6 17,2 87,0 10,0 13,9 12,3 7,4 6,6 9,4 3,0 19,1 53,6 17,8 K2O CT1.2: 0,2kg NPK CT2.1: 0,4kg P2O5 + 0,1kg K2O CT2.2: 0,2kg NPK CT3.1: 0,6kg P2O5 25 Cơng thức thí nghiệm (CTTN) TLS D1,3 Sd Hvn Sh Hdc Shdc Dt Sdt (%) (cm) (%) (m) (%) (m) (%) (m) (%) ΔM M (m3/ha) (m3/ /năm) + 0,1kg K2O CT3.2: 0,2kg NPK 88,8 9,7 17,9 12,1 7,2 6,6 9,5 2,8 20,5 50,1 16,7 83,3 9,8 14,4 12,4 4,5 6,5 8,6 2,8 18,0 52,6 17,5 83,3 9,2 19,0 12,3 5,6 6,8 9,4 2,8 25,0 45,4 15,1 9,3 15,8 11,7 6,4 6,5 7,1 2,7 24,6 45,1 15,0 CT4.1: 0,8kg P2O5 + 0,1kg K2O CT4.2: 0,2kg NPK CT5 Khơng bón 84,2 (ĐC) Sig F ă 0,0245 ổi (6/2016-6/2020) 0,034 0,058 0,048 26 Cơng thức thí nghiệm (CTTN) TLS D1,3 Sd Hvn Sh Hdc Shdc Dt Sdt (%) (cm) (%) (m) (%) (m) (%) (m) (%) ΔM M (m3/ha) (m3/ /năm) CT1.1: 0,2kg P2O5 + 0,1kg 88,8 11,9 14,1 14,3 10,6 8,0 85,1 10,7 16,5 13,4 9,3 7,6 86,1 12,5 16,3 14,4 8,7 8,1 85,1 11,7 17,2 14,0 9,5 7,9 9,5 3,2 21,5 78,5 19,6 3,4 23,4 56,7 14,1 3,4 21,5 84,0 21,0 3,5 22,2 71,0 17,7 3,5 18,3 74,0 18,5 K2O CT1.2: 0,2kg NPK 10,1 CT2.1: 0,4kg P2O5 + 0,1kg 10,0 K2O CT2.2: 0,2kg NPK CT3.1: 0,6kg P2O5 + 0,1kg 10,3 9,6 87,0 11,9 16,4 13,8 10,0 7,9 27 Cơng thức thí nghiệm (CTTN) TLS D1,3 Sd Hvn Sh Hdc Shdc Dt Sdt (%) (cm) (%) (m) (%) (m) (%) (m) (%) ΔM M (m3/ha) (m3/ /năm) K2O CT3.2: 0,2kg NPK 88,8 11,2 15,8 13,4 9,6 7,7 10,0 CT4.1: 3,3 22,1 65,8 16,4 3,3 19,6 68,2 17,0 3,3 18,3 61,8 15,4 3,2 18,7 49,0 12,2 107,2 21,4 8,7 0,8kg P2O5 + 0,1kg 83,3 11,6 17,6 14,0 8,2 7,9 83,3 11,3 16,4 13,3 11,1 7,6 K2O CT4.2: 0,2kg NPK CT5 9,2 7,3 Khơng bón 84,2 10,2 18,6 12,8 10,2 7,3 (ĐC) Sig F ă CT1.1: 0,038 0,016 0,013 0,027 ổi (6/2016-6/2021) 88,8 13,2 14,0 15,9 10,8 9,3 9,8 3,6 21,3 28 Cơng thức thí nghiệm (CTTN) TLS D1,3 Sd Hvn Sh Hdc Shdc Dt Sdt (%) (cm) (%) (m) (%) (m) (%) (m) (%) ΔM M (m3/ha) (m3/ /năm) 0,2kg P2O5 + 0,1kg K2O CT1.2: 0,2kg NPK 85,1 12,0 16,8 15,0 9,3 8,8 86,1 13,7 16,5 15,9 8,5 9,2 10,2 3,6 23,3 80,1 16,0 3,8 21,8 112,8 22,5 3,9 22,5 97,9 19,5 3,9 18,5 101,4 20,2 CT2.1: 0,4kg P2O5 + 0,1kg 10,1 K2O CT2.2: 0,2kg NPK 85,1 13,0 17,3 15,6 9,4 9,1 87,0 13,2 16,2 15,3 10,3 9,0 10,1 CT3.1: 0,6kg P2O5 + 0,1kg K2O 9,4 29 Cơng thức thí nghiệm (CTTN) CT3.2: 0,2kg NPK TLS D1,3 Sd Hvn Sh Hdc Shdc Dt Sdt (%) (cm) (%) (m) (%) (m) (%) (m) (%) 88,8 12,6 15,9 15,1 9,8 8,9 83,3 12,9 17,5 15,6 7,1 9,0 83,3 12,6 16,7 14,9 11,5 8,8 84,2 11,5 18,7 14,4 10,4 8,5 10,3 ΔM M (m3/ha) (m3/ /năm) 3,6 22,2 92,2 18,4 3,6 19,4 93,8 18,7 3,6 18,5 86,6 17,3 3,5 20,2 69,6 13,9 CT4.1: 0,8kg P2O5 + 0,1kg 8,5 K2O CT4.2: 0,2kg NPK 9,5 CT5 Khơng bón 7,0 (ĐC) Sig F 0,000 0,000 0,000 0,037 30 Phụ lục 3.5 Ảnh hƣởng biện ph p chăm sóc đến sinh trƣởng KLT Cơng thức thí nghiệm M TLS D1,3 Sd Hvn Sh Hdc Shdc Dt Sdt (%) (cm) (%) (m) (%) (m) (%) (m) (%) (CTTN) ă ổi (6/2016-6/2017) CT Chăm sóc lần; 93,3 1,3 44,0 1,3 36,1 0,7 42,0 96,3 1,3 41,3 1,4 30,7 0,7 41,2 95,3 1,1 37,1 1,3 35,1 0,6 39,3 0,5kg NPK CT Chăm sóc lần; 0,5kg NPK CT3 Chăm sóc lần; không phân Sig F ă 0,236 0,559 0,191 ổi (6/2016-6/2018) CT Chăm sóc lần; 88,8 4,5 41,1 6,0 19,5 2,5 22,2 (m3/ha) (m3/ ΔM /năm) 31 Công thức thí nghiệm M TLS D1,3 Sd Hvn Sh Hdc Shdc Dt Sdt (%) (cm) (%) (m) (%) (m) (%) (m) (%) 85 5,2 34,1 6,1 20,4 2,6 16,7 86,6 4,9 41,3 6,1 22,3 2,5 29,0 (m3/ha) (m3/ ΔM /năm) (CTTN) 0,5kg NPK CT Chăm sóc lần; 0,5kg NPK CT3 Chăm sóc lần; khơng phân Sig F ă 0,190 0,880 0,590 ổi (6/2016-6/2019) CT Chăm sóc lần; 88,8 9,5 15 11,8 12,5 6,8 11,0 2,6 13,1 47,0 15,6 84,2 9,7 13,8 11,7 13,4 6,8 10,0 3,0 12,8 48,5 16,1 0,5kg NPK CT Chăm sóc lần; 0,5kg NPK 32 Cơng thức thí nghiệm M TLS D1,3 Sd Hvn Sh Hdc Shdc Dt Sdt (%) (cm) (%) (m) (%) (m) (%) (m) (%) (m3/ha) (m3/ ΔM /năm) (CTTN) CT3 Chăm sóc lần; 86,1 9,4 15,0 11,7 8,8 6,8 11,2 2,9 10,0 45,9 15,3 không phân Sig F ă 0,004 0,049 0,052 0,169 ổi (6/2016-6/2020) CT Chăm sóc lần; 88,8 11,5 17,6 14,4 13,0 8,1 12,0 3,4 19,8 73,6 18,4 84,2 11,7 12,2 14,2 11,7 8,0 11,3 3,5 16,3 70,9 17,7 3,4 21,2 69,2 17,3 0,5kg NPK CT Chăm sóc lần; 0,5kg NPK CT3 Chăm sóc lần; 12,5 86,1 11,4 12,4 14,2 8,6 8,0 không phân Sig F 0,247 0,760 0,650 0,760 33 Công thức thí nghiệm M TLS D1,3 Sd Hvn Sh Hdc Shdc Dt Sdt (%) (cm) (%) (m) (%) (m) (%) (m) (%) (m3/ha) (m3/ ΔM /năm) (CTTN) ă ổi (6/2016-6/2021) CT Chăm sóc lần; 88,8 12,8 11,0 15,9 9,5 9,2 5,3 3,6 8,1 100,7 20,1 84,2 12,9 12,4 15,7 11,5 9,1 5,4 3,6 8,5 96,2 19,2 86,1 12,7 15,4 15,6 12,2 9,1 5,5 3,5 9,0 94,4 18,8 0,5kg NPK CT Chăm sóc lần; 0,5kg NPK CT3 Chăm sóc lần; khơng phân Sig F 0,24 0,76 0,62 0,76 34 BIỂU O ẾM RỪNG TRỒNG Địa điểm: …………………………………………………………………………………………………………………… Số hiệu ô:……….……….…Tiểu khu:……….………….… Khoảnh……… ………Lô………… …………… Tọa độ X:……………………………… ………… Tọa độ Y:…………………………….………………… Độ dốc:……………………… ……….………….Độ cao:…………………………………… ………………… Mật độ trồng: …………… …….… Cây/ha; Độ tàn che:……… %; Độ che phủ rừng ……… % Biện ph p KTLS t c động:…………………………………………………………………………………………… Kích thƣớc tiêu chuẩn:…………………………… …m2 TT 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 Loài D 1,3 (cm) Hvn (m) Hdc (m) Ngày đo ………………………………… Dt (m) Chất lượng (A,B,C) Ghi 35 MỘT SỐ H NH ẢNH T I HIỆN TR ỜNG NGHIÊN CỨU iều tr iện tr n mô n trồn Keo tr m Điều tra mô hình KLT Tun Quang Điều tra mơ hình KLT Hải Phòng 36 T u t ập số liệu 37 Mơ hình t í n iệm