Nghiên cứu khả năng hấp thụ co2 của rừng trồng keo tai tượng acacia mangium tại huyện định hóa tỉnh thái nguyên

ĐẠI HỌC THÁI NGUYÊN TRƢỜNG ĐẠI HỌC SƢ PHẠM HÀ THỊ DIỆU LINH NGHIÊN CỨU KHẢ NĂNG HẤP THỤ CO2 CỦA RỪNG TRỒNG KEO TAI TƢỢNG (Acacia mangium) TẠI HUYỆN ĐỊNH HÓA TỈNH THÁI NGUYÊN Chuyên ngành: Sinh thái học Mã số: 60.42.01.20 LUẬN VĂN THẠC SỸ SINH HỌC NGƢỜI HƢỚNG DẪN KHOA HỌC: TS Nguyễn Thế Hƣng PGS.TS Nguyễn Hồng Quang Thái Nguyên – Năm 2013 Số hóa Trung tâm Học liệu – Đại học Thái Nguyên http://www.lrc-tnu.edu.vn i LỜI CAM ĐOAN Tôi xin cam đoan cơng trình nghiên cứu riêng tơi với hƣớng dẫn khoa học TS Nguyễn Thế Hưng PGS.TS Hoàng Ngọc Quang (trƣờng Đại học Tài Nguyên Môi trƣờng Hà Nội) Các số liệu, kết nghiên cứu nêu luận văn hoàn toàn trung thực Nếu sai tơi chịu hồn tồn trách nhiệm Thái Ngun, ngày 10 tháng 04 năm 2013 Tác giả Hà Thị Diệu Linh Số hóa Trung tâm Học liệu – Đại học Thái Nguyên http://www.lrc-tnu.edu.vn ii LỜI CẢM ƠN Luận văn đƣợc hoàn thành trƣờng Đại học Sƣ phạm Thái Nguyên theo chƣơng trình đào tạo cao học Sinh học, chuyên ngành Sinh thái học, khoá 19 (2011 – 2013) Trong q trình thực hồn thành luận văn, tác giả nhận đƣợc quan tâm giúp đỡ Ban giám hiệu, khoa Sau đại học thầy, cô giáo Trƣờng Đại học Sƣ phạm, bạn bè đồng nghiệp cán địa phƣơng nơi tác giả thực nghiên cứu Nhân dịp này, tác giả xin chân thành cảm ơn giúp đỡ hiệu Trƣớc tiên, tác giả xin đặc biệt cảm ơn TS Nguyễn Thế Hưng PGS.TS Hoàng Ngọc Quang – ngƣời hƣớng dẫn khoa học, tận tình hƣớng dẫn, giúp đỡ tác giả trình thực luận văn Tác giả xin tỏ lòng cám ơn sâu sắc đến Ban giám hiệu trƣờng DHSP_ TN, Khoa sau đại học ĐHTN tạo điều kiện thuận lợi cho tác giả trình học tập nhƣ hoàn thành luận văn thạc sĩ Xin cảm ơn giúp đỡ nhiệt tình quan, ban ngành địa bàn tỉnh Thái Nguyên, xã số hộ dân trồng rừng địa bàn nghiên cứu tạo điều kiện giúp đỡ tác giả việc thu thập số liệu ngoại nghiệp để thực luận văn Xin trân trọng cảm ơn! Thái Nguyên, tháng năm 2013 Tác giả Số hóa Trung tâm Học liệu – Đại học Thái Nguyên http://www.lrc-tnu.edu.vn iii MỤC LỤC LỜI CAM ĐOAN i LỜI CẢM ƠN ii MỤC LỤC iii DANH MỤC HÌNH iv DANH MỤC BẢNG v DANH MỤC CÁC TỪ VIẾT TẮT iv MỞ ĐẦU Chƣơng TỔNG QUAN VẤN ĐỀ NGHIÊN CỨU 1.1 Trên giới 1.1.1 Nghiên cứu khả hấp thụ cacbon rừng 1.1.2 Nghiên cứu Keo tai tƣợng 1.2 Ở Việt Nam 1.2.1 Nghiên cứu khả hấp thụ carbon rừng 1.2.2 Nghiên cứu Keo tai tƣợng 11 Chƣơng MỤC TIÊU - ĐỐI TƢỢNG - NỘI DUNG VÀ PHƢƠNG PHÁP NGHIÊN CỨU 13 2.1 Mục tiêu nghiên cứu 13 2.2 Đối tƣợng nghiên cứu 13 2.3 Phạm vi nghiên cứu 13 2.4 Nội dung nghiên cứu 13 2.5 Cách tiếp cận phƣơng pháp nghiên cứu 14 2.5.1 Quan điểm cách tiếp cận đề tài 14 Chƣơng ĐẶC ĐIỂM TỰ NHIÊN, KINH TẾ- XÃ HỘI KHU VỰC NGHIÊN CỨU 22 Số hóa Trung tâm Học liệu – Đại học Thái Nguyên http://www.lrc-tnu.edu.vn 3.1 Điều kiện tự nhiên 22 3.1.1 Vị trí địa lý 22 3.1.2 Địa hình 22 3.1.3 Khí hậu, thủy văn 23 3.1.4 Thổ nhƣỡng 24 3.1.5 Hiện trạng rừng đất rừng 24 3.2 Điều kiện kinh tế- xã hội 25 3.2.1 Kinh tế 25 3.2.2 Xã hội 25 3.2.3 Cơ sở hạ tầng 26 Chƣơng KẾT QUẢ NGHIÊN CỨU VÀ THẢO LUẬN 28 4.1 Nghiên cứu sinh khối rừng trồng Keo tai tƣợng 28 4.1.1 Sinh khối cá thể Keo tai tƣợng 28 4.1.2 Sinh khối bụi, thảm tƣơi 32 4.1.3 Sinh khối tƣơi vật rơi rụng rừng trồng Keo tai tƣợng 33 4.1.4 Tổng sinh khối phần mặt đất toàn lâm phần 34 4.2 Nghiên cứu khả tích lũy carbon sinh khối phần mặt đất rừng trồng Keo tai tƣợng độ tuổi khác 40 4.2.1 Lƣợng carbon tích lũy cá thể Keo tai tƣợng 40 4.3 Lƣợng khí CO2 đƣợc hấp thụ để tích lũy sinh khối phần mặt đất rừng trồng Keo tai tƣợng 45 4.4 Nghiên cứu mối quan hệ sinh khối, lƣợng carbon tích lũy với nhân tố điều tra rừng chủ yếu 46 4.4.1 Mối quan hệ sinh khối khô Keo tai tƣợng với nhân tố điều tra 46 4.4.3 Mối quan hệ lƣợng carbon hấp thụ toàn lâm phần với nhân tố điều tra lâm phần 48 Số hóa Trung tâm Học liệu – Đại học Thái Nguyên http://www.lrc-tnu.edu.vn 4.5 Đề xuất số ứng dụng xác định lƣợng carbon tích lũy sinh khối phần mặt đất rừng trồng Keo tai tƣợng 49 4.5.1 Đề xuất số ứng dụng xác định lƣợng carbon đƣợc tích lũy cá thể Keo tai tƣợng dựa vào nhân tố điều tra 49 4.5.2 Đề xuất ứng dụng xác định carbon đƣợc tích lũy sinh khối phần mặt đất tồn lâm phần thơng qua nhân tố điều tra 49 KẾT LUẬN VÀ KHUYẾN NGHỊ 50 TÀI LIỆU THAM KHẢO 52 Phụ lục Số hóa Trung tâm Học liệu – Đại học Thái Ngun http://www.lrc-tnu.edu.vn iv DANH MỤC HÌNH Hình 2.1 : Sơ đồ bƣớc tiến hành nghiên cứu đề tài 15 Hình 2.2: Sơ đồ bố trí tiêu chuẩn 16 Hình 4.1: Cấu trúc sinh khối cá thể keo tai tƣợng trồng loài xã Đồng Thịnh 29 Hình 4.2: Cấu trúc sinh khối cá thể keo tai tƣợng trồng loài xã Bảo Linh 30 Hình 4.3: Cấu trúc sinh khối tƣơi phần mặt đất trồng Keo tai tƣợng xã Đồng Thịnh xã Bảo Linh 35 Số hóa Trung tâm Học liệu – Đại học Thái Nguyên http://www.lrc-tnu.edu.vn v DANH MỤC BẢNG Bảng 1: Tổng hợp số liệu phục vụ cho đề tài 21 Bảng 3.1: Thống kê diện tích đất đai khu vực nghiên cứu 24 Bảng 4.1: Cấu trúc sinh khối tƣơi phần mặt đất cá thể Keo tai tƣợng rừng trồng xã Đồng Thịnh độ tuổi khác 28 Bảng 4.2: Cấu trúc sinh khối tƣơi phần mặt đất cá thể Keo tai tƣợng rừng trồng xã Bảo Linh độ tuổi khác 28 Bảng 4.3: Tỷ lệ sinh khối khô với sinh khối tƣơi Keo tai tƣợng 31 Bảng 4.4 Mối quan hệ tổng sinh khối khô với tổng sinh khối tƣơi Keo tai tƣợng 32 Bảng 4.5: Sinh khối tƣơi bụi, thảm tƣơi rừng Keo tai tƣợng xã Đồng Thịnh xã Bảo Linh (tấn/ha) 33 Bảng 4.6: Sinh khối vật rơi rụng rừng Keo tai tƣợng xã Đồng Thịnh xã Bảo Linh (tấn/ha) 33 Bảng 4.7: Sinh khối tƣơi phần mặt đất lâm phần rừng Keo 34 tai tƣợng trồng loài theo độ tuổi khác xã Bảo Linh 34 Bảng 4.8: Sinh khối tƣơi phần mặt đất lâm phần rừng Keo tai tƣợng trồng loài theo độ tuổi khác xã Đồng Thịnh 34 Bảng 4.9: Sinh khối khô phần mặt đất bụi, thảm tƣơi rừng trồng Keo tai tƣợng độ tuổi khác xã Đồng Thịnh xã Bảo Linh 37 Bảng 4.10: Sinh khối khô phần mặt đất vật rơi rụng rừng trồng Keotai tƣợng độ tuổi khác xã Đồng Thịnh xã Bảo Linh 37 Bảng 4.11: Sinh khối khô phần mặt đất Keo tai tƣợng rừng trồng xã Đồng Thịnh 38 Số hóa Trung tâm Học liệu – Đại học Thái Nguyên http://www.lrc-tnu.edu.vn Bảng 4.12: Sinh khối khô phần mặt đất Keo tai tƣợng rừng trồng xã Bảo Linh 38 Bảng 4.13: Sinh khối khô phần mặt đất rừng trồng Keo tai tƣợng độ tuổi khác xã Đồng Thịnh 39 Bảng 4.14: Sinh khối khô phần mặt đất rừng trồng Keo tai tƣợng độ tuổi khác xã Bảo Linh 40 Bảng 4.15: Kết phân tích hàm lƣợng carbon tích lũy Keo 40 tai tƣợng (%) 40 Bảng 4.16: Hàm lƣợng carbon tích lũy sinh khối phần mặt đât bụi thảm tƣơi vật rơi rụng rừng trồng Keo tai tƣợng (%) 41 Bảng 4.17: Lƣợng carbon tích lũy sinh khối phần mặt đất rừng trồng Keo tai tƣợng xã Đồng Thịnh 42 Bảng 4.18: Lƣợng carbon tích lũy sinh khối phần mặt đất rừng trồng Keo tai tƣợng xã Bảo Linh 42 Bảng 4.19: Mối quan hệ tổng sinh khối khô cá thể Keo tai tƣợng 44 Bảng 4.20: Mối quan hệ tổng sinh khối tƣơi cá thể với đƣờng kính ngang ngực (D1.3) 44 Bảng 4.21: Lƣợng CO2 đƣợc hấp thu trình tạo sinh khối phần mặt đất rừng trồng Keo tai tƣợng xã Đồng Thịnh xã Bảo Linh 46 Bảng 4.22: Mối quan hệ sinh khối khô Keo tai tƣợng với nhân tố điều tra lâm phần độ tuổi khác 46 4.4.2 Mối quan hệ tổng lƣợng carbon tích lũy cá thể với nhân tố điều tra (D1.3) 47 Bảng 4.23 Mối quan hệ tổng lƣợng carbon tích lũy cá thể Keo tai tƣợng với nhân tố điều tra (D1.3) 47 Bảng 4.24: Mối quan hệ lƣợng carbon hấp thụ toàn lâm phần với nhân tố điều tra lâm phần 48 Số hóa Trung tâm Học liệu – Đại học Thái Nguyên http://www.lrc-tnu.edu.vn vi DANH MỤC CÁC TỪ VIẾT TẮT Kí hiệu TT Ghi viết tắt D1.3 Đƣờng kính ngang ngực Hvn Chiều cao vút Lt Chiều cao tán Dt Đƣờng kính tán Hdc Chiều cao dƣới cành N Mật độ MiCB Carbon tích lũy phận i bụi MCB/ha Carbon tích lũy bụi/ha MVRR Carbon tích lũy vật rơi rụng 10 MVRR/ha Carbon tích lũy vật rơi rụng/ 11 MiCCT Carbon tích lũy phận i cá thể 12 A Tuổi 13 CDM Cơ chế phát triển 14 UNFCCC Cơng ƣớc chống biến đổi khí hậu tồn cầu (United Framework Convention Climate Change) 15 IPCC Uỷ ban liên Chính Phủ biến đổi khí hậu ( the intergovermental Panel on Climate Change ) 16 AR-CDM Trồng rừng theo chế phát triển 17 R Hệ số tƣơng quan 18 S Thể tích 19 TDW Tổng khối lƣợng khơ 20 SFW Khối lƣợng tƣơi Số hóa Trung tâm Học liệu – Đại học Thái Nguyên http://www.lrc-tnu.edu.vn 70 b ANOVA Model Sum of Squares Regression Mean Square 123.561 123.561 14.847 4.949 138.408 Residual Total df F Sig 24.966 015 a a Predictors: (Constant), D1.3 b Dependent Variable: Pt Coefficients a Standardized Unstandardized Coefficients Model B (Constant) D1.3 Coefficients Std Error Beta 122.635 13.016 5.368 1.074 t 945 Sig 9.422 003 4.997 015 a Dependent Variable: Pt Model Summary and Parameter Estimates Dependent Variable:Pt Model Summary Equation Linear R Square 893 F 24.966 df1 Parameter Estimates df2 Sig Constant 015 122.635 b1 5.368 The independent variable is D1.3 Số hóa Trung tâm Học liệu – Đại học Thái Nguyên http://www.lrc-tnu.edu.vn 71 Phụ lục : Mối quan hệ tổng sinh khối khô với sinh khối tươi cá thể Keo tai tượng Tuổi Model Summary Model R Adjusted R Std Error of the Square Estimate R Square 954 a 909 864 19585 a Predictors: (Constant), Pt b ANOVA Model Sum of Squares Df Mean Square Regression 771 771 Residual 077 038 Total 847 F Sig 20.095 046 a a Predictors: (Constant), Pt b Dependent Variable: Pk Coefficients a Standardized Unstandardized Coefficients Model B (Constant) Pt Std Error 14.700 3.063 240 053 Coefficients Beta t 954 Sig 4.798 041 4.483 046 a Dependent Variable: Pk Số hóa Trung tâm Học liệu – Đại học Thái Nguyên http://www.lrc-tnu.edu.vn 72 Tuổi 4: Model Summary Model R 966 Adjusted R Std Error of the Square Estimate R Square a 933 899 46368 a Predictors: (Constant), Pt b ANOVA Model Sum of Squares Regression Residual Total df Mean Square 5.980 5.980 430 215 6.410 F 27.814 Sig .034 a Predictors: (Constant), Pt b Dependent Variable: Pk Số hóa Trung tâm Học liệu – Đại học Thái Nguyên http://www.lrc-tnu.edu.vn a 73 Coefficients a Standardized Unstandardized Coefficients Model B (Constant) Coefficients Std Error Beta -75.050 26.433 1.000 190 Pt t 966 Sig -2.839 105 5.274 034 a Dependent Variable: Pk Tuổi 5: Model Summary Model R 983 R Square a Adjusted R Std Error of the Square Estimate 967 935 96333 a Predictors: (Constant), Pt b ANOVA Model Sum of Squares Regression Residual df Mean Square 27.412 27.412 928 928 Số hóa Trung tâm Học liệu – Đại học Thái Nguyên F 29.538 Sig .116 http://www.lrc-tnu.edu.vn a 74 Total 28.340 a Predictors: (Constant), Pt b Dependent Variable: Pk Coefficients a Standardized Unstandardized Coefficients Model B (Constant) Pt Std Error -1.006 16.625 478 088 Coefficients Beta t 983 Sig -.061 962 5.435 116 a Dependent Variable: Pk Số hóa Trung tâm Học liệu – Đại học Thái Nguyên http://www.lrc-tnu.edu.vn 75 Phụ lục 9: Mối quan hệ tổng lượng carbon tích lũy cá thể với nhân tố điều tra Tuổi 3: Variables Entered/Removed Variables Variables Entered Removed Model D1.3 b Method a Enter a All requested variables entered b Dependent Variable: Cz Model Summary Model R Adjusted R Std Error of the Square Estimate R Square 941 a 886 847 06177 a Predictors: (Constant), D1.3 b ANOVA Model Sum of Squares df Mean Square Regression 089 089 Residual 011 004 Total 100 F Sig 23.211 017 a a Predictors: (Constant), D1.3 b Dependent Variable: Cz Coefficients a Standardized Unstandardized Coefficients Model B (Constant) D1.3 Std Error 3.689 335 211 044 Số hóa Trung tâm Học liệu – Đại học Thái Nguyên Coefficients Beta t 941 Sig 10.996 002 4.818 017 http://www.lrc-tnu.edu.vn 76 Coefficients a Standardized Unstandardized Coefficients Model B (Constant) Coefficients Std Error Beta 3.689 335 211 044 D1.3 t Sig 10.996 002 4.818 017 941 a Dependent Variable: Cz Model Summary and Parameter Estimates Dependent Variable:Cz Model Summary Equation R Square Linear 886 F df1 23.211 Parameter Estimates df2 Sig Constant 017 3.689 b1 211 The independent variable is D1.3 Tuổi 4: Model Summary Model R 987 R Square a Adjusted R Std Error of the Square Estimate 974 Số hóa Trung tâm Học liệu – Đại học Thái Nguyên 966 02922 http://www.lrc-tnu.edu.vn 77 Model Summary Model R Adjusted R Std Error of the Square Estimate R Square 987 a 974 966 02922 a Predictors: (Constant), D1.3 b ANOVA Model Sum of Squares df Mean Square Regression 097 097 Residual 003 001 Total 100 F Sig 114.083 002 a a Predictors: (Constant), D1.3 b Dependent Variable: Cz Coefficients a Standardized Unstandardized Coefficients Model B (Constant) Coefficients Std Error Beta 5.973 144 132 012 D1.3 t 987 Sig 41.592 000 10.681 002 a Dependent Variable: Cz Model Summary and Parameter Estimates Dependent Variable:Cz Model Summary Equation Linear R Square 974 F 114.083 df1 Parameter Estimates df2 Sig Constant 002 5.973 b1 132 The independent variable is D1.3 Số hóa Trung tâm Học liệu – Đại học Thái Nguyên http://www.lrc-tnu.edu.vn 78 Tuổi 5: Model Summary Model R 958 Adjusted R Std Error of the Square Estimate R Square a 918 891 05218 a Predictors: (Constant), D1.3 b ANOVA Model Sum of Squares df Mean Square Regression 092 092 Residual 008 003 Total 100 F 33.725 Sig .010 a Predictors: (Constant), D1.3 b Dependent Variable: Cz Số hóa Trung tâm Học liệu – Đại học Thái Nguyên http://www.lrc-tnu.edu.vn a 79 Coefficients a Standardized Unstandardized Coefficients Model B (Constant) Coefficients Std Error Beta 7.484 314 151 026 D1.3 t 958 Sig 23.874 000 5.807 010 a Dependent Variable: Cz Model Summary and Parameter Estimates Dependent Variable:Cz Model Summary Equation Linear R Square 918 F 33.725 df1 Parameter Estimates df2 Sig Constant 010 7.484 b1 151 The independent variable is D1.3 Số hóa Trung tâm Học liệu – Đại học Thái Nguyên http://www.lrc-tnu.edu.vn 80 Phụ lục 10: Mối quan hệ tổng lƣợng carbon tích lũy cá thể với nhân tố điều tra Tuổi 3: Variables Entered/Removed Variables Variables Entered Removed Model Pzk b Method a Enter a All requested variables entered b Dependent Variable: Cz Model Summary Model R 941 Adjusted R Std Error of the Square Estimate R Square a 886 847 27568 a Predictors: (Constant), Pzk b ANOVA Model Sum of Squares Regression Residual Total df Mean Square 1.764 1.764 228 076 1.992 F 23.211 Sig .017 a Predictors: (Constant), Pzk b Dependent Variable: Cz Số hóa Trung tâm Học liệu – Đại học Thái Nguyên http://www.lrc-tnu.edu.vn a 81 Coefficients a Standardized Unstandardized Coefficients Model B (Constant) Pzk Coefficients Std Error Beta -9.420 3.751 4.200 872 t 941 Sig -2.512 087 4.818 017 a Dependent Variable: Cz Model Summary and Parameter Estimates Dependent Variable:Cz Model Summary Equation Linear R Square 886 F 23.211 df1 Parameter Estimates df2 Sig Constant 017 -9.420 b1 4.200 The independent variable is Pzk Số hóa Trung tâm Học liệu – Đại học Thái Nguyên http://www.lrc-tnu.edu.vn 82 Tuổi 4: Model Summary Model R Adjusted R Std Error of the Square Estimate R Square 987 a 974 966 21909 a Predictors: (Constant), Pzk b ANOVA Model Sum of Squares Regression Residual Total df Mean Square 5.476 5.476 144 048 5.620 F Sig 114.083 002 a a Predictors: (Constant), Pzk b Dependent Variable: Cz Coefficients a Standardized Unstandardized Coefficients Model B (Constant) Pzk Coefficients Std Error -50.300 5.890 7.400 693 Beta t 987 Sig -8.540 003 10.681 002 a Dependent Variable: Cz Model Summary and Parameter Estimates Dependent Variable:Cz Model Summary Equatio n Linear R Square 974 F 114.083 df1 Parameter Estimates df2 Số hóa Trung tâm Học liệu – Đại học Thái Nguyên Sig Constant 002 -50.300 b1 7.400 http://www.lrc-tnu.edu.vn 83 Model Summary and Parameter Estimates Dependent Variable:Cz Model Summary Equatio n R Square Linear 974 F df1 114.083 Parameter Estimates df2 Sig Constant 002 -50.300 b1 7.400 The independent variable is Pzk Tuổi 5; Model Summary Model R 958 Adjusted R Std Error of the Square Estimate R Square a 918 891 33216 a Predictors: (Constant), Pzk b ANOVA Model Sum of Squares Regression Residual Total df Mean Square 3.721 3.721 331 110 4.052 Số hóa Trung tâm Học liệu – Đại học Thái Nguyên F 33.725 Sig .010 http://www.lrc-tnu.edu.vn a 84 b ANOVA Model Sum of Squares Regression Residual Total df Mean Square 3.721 3.721 331 110 4.052 F Sig 33.725 010 a a Predictors: (Constant), Pzk b Dependent Variable: Cz Coefficients a Standardized Unstandardized Coefficients Model B (Constant) Std Error Beta -25.370 6.619 6.100 1.050 Pzk Coefficients t 958 Sig -3.833 031 5.807 010 a Dependent Variable: Cz Model Summary and Parameter Estimates Dependent Variable:Cz Model Summary Equation Linear R Square 918 F 33.725 df1 Parameter Estimates df2 Sig Constant 010 -25.370 b1 6.100 The independent variable is Pzk Số hóa Trung tâm Học liệu – Đại học Thái Nguyên http://www.lrc-tnu.edu.vn ... tƣợng nghiên cứu Rừng Keo tai tƣợng loài theo giai đoạn từ 3, 4, tuổi huyện Định Hóa tỉnh Thái Nguyên 2.3 Phạm vi nghiên cứu -Về địa bàn nghiên cứu: Nghiên cứu khả hấp thụ carbon rừng trồng Keo tai. .. khả hấp thụ carbon rừng trồng Keo tai tƣợng lồi huyện Định Hóa tỉnh Thái Nguyên Đề xuất số ứng dụng việc xác định lƣợng carbon hấp thụ rừng trồng Keo tai tƣợng lồi huyện Định Hóa tỉnh Thái Ngun... 4, huyện Định Hóa Tĩnh Thái Nguyên Nghiên cứu khả hấp thụ carbon số quần xã rừng Keo tai tƣợng trồng loài tuổi 3, 4, huyện Định Hóa Tĩnh Thái Nguyên Số hóa Trung tâm Học liệu – Đại học Thái Nguyên