Nghiên cứu ảnh hưởng của mật độ và kích cỡ thả nuôi đến sinh trưởng, tỷ lệ sống của nghêu (meretrix lyrata sowerby, 1851) giai đoạn nuôi thương phẩm tại vùng bãi triều, kiến thụy, hải phòng

67 6 0
Nghiên cứu ảnh hưởng của mật độ và kích cỡ thả nuôi đến sinh trưởng, tỷ lệ sống của nghêu (meretrix lyrata sowerby, 1851) giai đoạn nuôi thương phẩm tại vùng bãi triều, kiến thụy, hải phòng

Đ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

BỘ GIÁO DỤC VÀ ĐÀO TẠO TRƯỜNG ĐẠI HỌC NHA TRANG NGÔ XUÂN BA NGHIÊN CỨU ẢNH HƯỞNG CỦA MẬT ĐỘ VÀ KÍCH CỠ THẢ NI ĐẾN SINH TRƯỞNG, TỶ LỆ SỐNG CỦA NGHÊU (Meretrix lyrata SOWERBY, 1851) GIAI ĐOẠN NUÔI THƯƠNG PHẨM TẠI VÙNG BÃI TRIỀU, KIẾN THUỴ, HẢI PHỊNG LUẬN VĂN THẠC SĨ Khánh Hồ - 2015 BỘ GIÁO DỤC VÀ ĐÀO TẠO TRƯỜNG ĐẠI HỌC NHA TRANG NGÔ XUÂN BA NGHIÊN CỨU ẢNH HƯỞNG CỦA MẬT ĐỘ VÀ KÍCH CỠ THẢ NI ĐẾN SINH TRƯỞNG, TỶ LỆ SỐNG CỦA NGHÊU (Meretrix lyrata SOWERBY, 1851) GIAI ĐOẠN NUÔI THƯƠNG PHẨM TẠI VÙNG BÃI TRIỀU, KIẾN THUỴ, HẢI PHÒNG LUẬN VĂN THẠC SĨ Ngành: Nuôi trồng thuỷ sản Mã số: 60 62 03 01 NGƯỜI HƯỚNG DẪN KHOA HỌC TS NGUYỄN TẤN SỸ CHỦ TỊCH HỘI ĐỒNG KHOA SAU ĐẠI HỌC PGS.TS LẠI VĂN HÙNG HOÀNG HÀ GIANG Khánh Hồ - 2015 i LỜI CAM ĐOAN Tơi xin cam đoan rằng, số liệu kết nghiên cứu luận văn trung thực chưa sử dụng để bảo vệ học vị Tôi xin cam đoan, giúp đỡ cho việc thực luận văn cảm ơn, thông tin trích dẫn luận văn rõ nguồn gốc Tác giả Ngô Xuân Ba ii LỜI CẢM ƠN Xin trân trọng gửi lời cảm ơn tới Trường Đại học Nha Trang, Ban lãnh đạo Sở Nông nghiệp & Phát triển nông thôn Hải Phòng, Phòng Nuôi trồng thủy sản - Sở NN&PTNT Hải Phòng tập thể lãnh đạo, cán Dự án ICA (Dự án Tác động biến đổi khí hậu giải pháp thích ứng an tồn sinh học cho ni trồng thủy sản miền Bắc) giúp đỡ tơi hồn thành khố học Tơi xin dành biết ơn sâu sắc tới thầy cô truyền thụ cho kiến thức nhất, đặc biệt giáo viên hướng dẫn - TS Nguyễn Tấn Sỹ tận tình bảo giúp đỡ tơi q trình thực luận văn này; Đồng thời tơi xin gửi lời cảm ơn đến cán bộ, nhân viên Sở NN&PTNT Hải Phòng, phịng Nơng nghiệp Phát triển nơng thôn huyện Kiến Thuỵ, tạo điều kiện thuận lợi thời gian học tập thực luận văn Tôi xin gởi lời cảm ơn chân thành đến phòng Đào tạo sau Đại học - Viện Nghiên cứu Hải sản, Khoa sau đại học & Viện Nuôi trồng thủy sản - Trường Đại học Nha Trang giúp đỡ trình học tập thực luận văn tốt nghiệp hoàn tất thủ tục bảo vệ luận văn Lời cảm ơn chân thành xin dành cho gia đình bạn bè, đồng nghiệp giúp đỡ, động viên, cổ vũ tơi q trình học tập cơng tác Hải Phịng, ngày tháng Tác giả Ngơ Xn Ba năm 2014 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, CHỮ VIẾT TẮT v DANH MỤC CÁC BẢNG vi DANH MỤC CÁC HÌNH VẼ, ĐỒ THỊ vii MỞ ĐẦU CHƯƠNG - TỔNG QUAN TÀI LIỆU 1.1 Một số đặc điểm sinh học đối tượng nghiên cứu 1.1.1 Đặc điểm phân loại 1.1.2 Đặc điểm cấu tạo 1.1.3 Phân bố nghêu 1.1.4 Đặc điểm dinh dưỡng .6 1.1.5 Đặc điểm sinh trưởng .6 1.1.6 Đặc điểm sinh sản 1.2 Tình hình nghiên cứu, ni nghêu giới 1.2.1 Tình hình nuôi nghêu 1.2.2 Nghiên cứu ảnh hưởng yếu tố môi trường 1.2.3 Khả tích tụ chất nhiễm, độc tố nghêu .10 1.2.4 Nghiên cứu ảnh hưởng cố, suy thối nhiễm mơi trường gây chết nghêu 11 1.2.5 Giải pháp phòng tránh, giảm thiểu ảnh hưởng môi trường tới nghêu nuôi 11 1.2.6 Nghiên cứu sản xuất giống, nuôi quản lý 12 1.2.7 Nghiên cứu dịch bệnh 12 1.2.8 Nghiên cứu mật độ, kích cỡ thả nuôi tới sinh trưởng nghêu .12 1.3 Tình hình nghiên cứu, ni nghêu ở Việt Nam 13 1.3.1 Tình hình ni 13 1.3.2 Nghiên cứu nguồn lợi, phân bố 13 1.3.3 Nghiên cứu dinh dưỡng, sinh trưởng 14 1.3.4 Nghiên cứu đặc điểm sinh học, sản xuất giống 14 1.3.5 Nghiên cứu ảnh hưởng yếu tố môi trường 15 1.3.6 Nghiên cứu giải pháp phòng tránh, giảm thiểu ảnh hưởng môi trường tới nghêu nuôi 16 iv 1.3.7 Nghiên cứu dịch bệnh 16 1.3.8 Nghiên cứu ảnh hưởng mật độ ni, kích thước thả giống đến sinh trưởng, tỷ lệ sống nghêu 17 1.3.9 Thực trạng nuôi nghêu Hải Phòng 17 CHƯƠNG - VẬT LIỆU VÀ PHƯƠNG PHÁP NGHIÊN CỨU 20 2.1 Địa điểm, thời gian, đối tượng nghiên cứu 20 2.2 Phương pháp nghiên cứu 20 2.2.1 Sơ đồ khối nội dung nghiên cứu 20 2.2.2 Phương pháp bố trí thí nghiệm 21 2.3 Phương pháp thu thập, phân tích xử lý số liệu 22 2.3.1 Phương pháp xác định thông số môi trường 22 2.3.2 Phương pháp xác định tiêu theo dõi 22 2.3.3 Phương pháp xử lý số liệu 23 CHƯƠNG - KẾT QUẢ NGHIÊN CỨU VÀ THẢO LUẬN 26 3.1 Ảnh hưởng mật độ ni, kích cỡ thả ni đến tăng trưởng, tỷ lệ sống nghêu 26 3.1.1 Biến động yếu tố mơi trường thí nghiệm 26 3.1.2 Tốc độ tăng trưởng nghêu 30 3.1.3 Tỷ lệ sống nghêu 34 3.2 Đánh giá hiệu kinh tế 35 KẾT LUẬN VÀ KIẾN NGHỊ 38 Kết luận 38 Kiến nghị 38 TÀI LIỆU THAM KHẢO 39 v DANH MỤC CÁC KÝ HIỆU, CHỮ VIẾT TẮT ctv Cộng tác viên cm Centimeter DO Hàm lượng ơxy hồ tan DLG Tăng trưởng chiều dài tuyệt đối DWG Tăng trưởng khối lượng tuyệt đối ĐBSCL Đồng sông Cửu Long g gam GHCP Giới hạn cho phép L Chiều dài M Mét Mm Milimet mg miligam NT Nghiệm thức QCVN Quy chuẩn Việt Nam S‰ Độ mặn SGRW Tốc độ tăng trưởng tương đối khối lượng SGRL Tốc độ tăng trưởng tương đối chiều dài Tr đ Triệu đồng W Khối lượng vi DANH MỤC CÁC BẢNG Trang Bảng 1.1: Hiện trạng nuôi nghêu địa bàn thành phố Hải Phòng 17 Bảng 3.1: Điều kiện môi trường thời gian thí nghiệm 27 Bảng 3.2: Tốc độ tăng trưởng tuyệt đối chiều dài (mm/ngày) nghêu nuôi 32 Bảng 3.3: Tốc độ tăng trưởng tuyệt đối khối lượng (g/ngày) nghêu nuôi 33 Bảng 3.4: Tỷ lệ sống nghêu ở nghiệm thức thí nghiệm 35 Bảng 3.5: Hạch toán hiệu kinh tế nghêu nuôi ở nghiệm thức 36 vii DANH MỤC CÁC HÌNH VẼ, ĐỒ THỊ Trang Hình 1.1: Đặc điểm cấu tạo nghêu Hình 1.2: Cấu tạo nghêu Hình 2.1: Sơ đồ khối nội dung nghiên cứu 20 Hình 2.2: Bố trí thí nghiệm ni 21 Hình 2.3: Đo chiều dài, khối lượng nghêu 22 Hình 3.1: Sự biến động nhiệt độ thời gian thí nghiệm 28 Hình 3.2: Biến động pH thời gian thí nghiệm 29 Hình 3.3: Biến động độ mặn thời gian thí nghiệm 29 Hình 3.4: Biến động hàm lượng ơxy hồ tan thời gian thí nghiệm .30 Hình 3.5: Tương quan chiều dài khối lượng nghêu 30 Hình 3.6: Tốc độ tăng trưởng chiều dài đặc trưng nghêu (SRGL) 31 Hình 3.7: Tốc độ tăng trưởng khối lượng đặc trưng nghêu (SRGW) 33 MỞ ĐẦU Động vật thân mềm ngành có số lượng lồi lớn, khoảng 130.000 loài, chiếm 10% tổng số loài động vật trái đất, đứng thứ hai sau ngành giáp xác (với 800.000 loài) [3] Những năm gần đây, nghề nuôi động vật thân mềm trọng quan tâm Nghêu Bến Tre Meretrix lyrata động vật thân mềm có giá trị dinh dưỡng cao, chứa khoảng 56% protein tính theo trọng lượng khơ [15] Nghêu Bến Tre sinh trưởng nhanh, có giá trị kinh tế cao đối tượng thuỷ sản nuôi chủ lực tỉnh thành ven biển Nghêu Bến Tre phân bố tự nhiên Đồng sông Cửu Long Những năm gần nghêu M lyrata di nhập vùng triều ven biển tỉnh miền Bắc (Thái Bình, Hải Phòng, Nam Định ) để nuôi Xác định nuôi trồng thủy sản mạnh hướng chiến lược ngành thủy sản, Hải Phòng chủ trương phát triển nuôi trồng thủy sản thành ngành sản xuất hàng hóa lớn với đối tượng ni mới, có giá trị kinh tế Trong năm gần đây, nông dân vùng ven biển Hải Phòng đầu tư vào nuôi nhuyễn thể ở vùng bãi triều, chủ yếu nghêu Meretrix lyrata Kiến Thụy huyện đồng ven biển thành phố Hải Phòng, có xã Đại Hợp với diện tích rừng ngập mặn 500 ha, diện tích cồn cát Quần Mục phía ngồi rừng ngập mặn 1.000 ha, có điều kiện thuận lợi để hình thành vùng nuôi nghêu Bãi triều nơi vùng hạ triều, thời gian ngập nước dài, đáy cát pha bùn, đáy mềm, xốp, sinh vật phù du phong phú Hiện nay, diện tích quy mơ đầu tư ngày mở rộng phát triển lớn thị trường lại có nhu cầu sử dụng sản phẩm ngày gia tăng, không tiêu thụ nội địa mà còn xuất sang Trung Quốc, Đài Loan Tuy nhiên, nghêu nuôi vùng triều đối tượng chịu nhiều ảnh hưởng tác động biến đổi khí hậu Sự tăng nhiệt độ bề mặt nước biển dẫn đến giảm lượng dưỡng chất nước, giảm tốc độ tăng trưởng sản lượng Các mầm bệnh ký sinh trùng phát triển với gia tăng nhiệt độ Những năm gần xảy tượng nghêu chết hàng loạt, có kết hợp yếu tố bất lợi như: thời tiết nắng nóng, nhiệt độ độ mặn cao, mật độ thả nuôi cao Mật độ thả nuôi dày dẫn đến thiếu thức ăn, đồng thời nguyên nhân gia tăng tác nhân gây bệnh nghêu PHỤ LỤC Thí nghiệm kích cỡ * Tỷ lệ sống 1.1 Ảnh hưởng tới sinh trưởng chiều dài Between-Subjects Factors Value Label N Giong 300 12 2 400 80 120 160 200 12 6 6 MD Descriptive Statistics Dependent Variable:CD Std Deviation Giong MD Mean 300 80 062667 0007234 120 059967 0006506 160 057033 0005859 200 052567 0008083 058058 067600 062833 060133 056833 061850 065133 061400 058583 054700 059954 0039560 0010536 0005132 0004041 0004041 0041542 0028204 0016553 0017566 0024058 0044146 12 3 3 12 6 6 24 Total 400 80 120 160 200 Total Total 80 120 160 200 Total N Tests of Between-Subjects Effects Dependent Variable:CD Type III Sum of Squares Source Mean Square df Corrected 000a 6.299E-5 Model Intercept 086 086 Giong 8.626E-5 8.626E-5 MD 000 000 Giong * MD 4.295E-6 1.432E-6 Error 7.287E-6 16 4.554E-7 Total 087 24 Corrected 000 23 Total a R Squared = 984 (Adjusted R Squared = 977) Estimated Marginal Means F 138.320 000 1.894E5 189.410 256.467 3.143 000 000 000 054 Grand Mean Dependent Variable:CD 95% Confidence Interval Std Error Mean 060 Lower Bound 000 Upper Bound 060 060 Giong Dependent Variable:CD 95% Confidence Interval Giong Mean 300 400 058 062 Std Error 000 000 Lower Bound 058 061 Sig Upper Bound 058 062 MD Dependent Variable:CD 95% Confidence Interval MD Std Error Mean 80 120 160 200 065 061 059 055 Lower Bound 000 000 000 000 Upper Bound 065 061 058 054 066 062 059 055 Giong * MD Dependent Variable:CD 95% Confidence Interval Giong MD 300 400 Mean Std Error Lower Bound Upper Bound 80 063 000 062 063 120 060 000 059 061 160 057 000 056 058 200 80 120 160 200 053 068 063 060 057 000 000 000 000 000 052 067 062 059 056 053 068 064 061 058 Oneway Descriptives CD 95% Confidence Interval for Mean Std Mean Deviation Std Error N 3 3 3 062667 059967 057033 052567 067600 062833 0007234 0006506 0005859 0008083 0010536 0005132 0004177 0003756 0003383 0004667 0006083 0002963 Lower Bound Upper Bound 060870 058350 055578 050559 064983 061559 064464 061583 058489 054575 070217 064108 Minimum Maximum 0622 0593 0566 0517 0666 0624 0635 0606 0577 0533 0687 0634 Tota l 060133 0004041 0002333 056833 0004041 0002333 059129 055829 061137 057837 0597 0566 0605 0573 24 059954 0044146 0009011 058090 061818 0517 0687 Test of Homogeneity of Variances CD Levene Statistic df1 672 df2 Sig 16 693 ANOVA CD Sum of Squares Between Groups Within Groups Total Post Hoc Tests Mean Square df 000 000 000 16 23 F Sig .000 138.320 000 000 Homogeneous Subsets CD Duncan Subset for alpha = 0.05 NT Sig N 052567 056833 057033 059967 060133 062667 062833 067600 1.000 721 766 766 1.000 1.2 Ảnh hưởng tới tỷ lệ sống Between-Subjects Factors Value Label N Giong 300 12 2 400 80 120 160 200 12 6 6 MD Descriptive Statistics Dependent Variable:KL Std Deviation Giong MD Mean 300 80 039467 0010786 120 037333 0005859 160 034367 0005508 200 031367 0004509 035633 041400 039700 037300 034700 038275 040433 038517 035833 033033 036954 0032503 0006245 0004583 0003606 0003464 0026677 0013201 0013790 0016597 0018608 0032057 12 3 3 12 6 6 24 Total 400 80 120 160 200 Total Total 80 120 160 200 Total N Tests of Between-Subjects Effects Dependent Variable:KL Type III Sum of Squares Source Mean Square df F Corrected 000a 3.295E-5 92.054 Model Intercept 033 033 9.157E4 Giong 4.187E-5 4.187E-5 116.984 MD 000 6.235E-5 174.204 Giong * MD 1.711E-6 5.704E-7 1.594 Error 5.727E-6 16 3.579E-7 Total 033 24 Corrected 000 23 Total a R Squared = 976 (Adjusted R Squared = 965) Estimated Marginal Means Grand Mean Dependent Variable:KL 95% Confidence Interval Mean Std Error 037 Lower Bound 000 Upper Bound 037 037 Giong Dependent Variable:KL 95% Confidence Interval Giong Mean 300 400 036 038 Std Error 000 000 Lower Bound 035 038 Upper Bound 036 039 Sig .000 000 000 000 230 MD Dependent Variable:KL 95% Confidence Interval MD Std Error Mean 80 120 160 200 040 039 036 033 Lower Bound 000 000 000 000 Upper Bound 040 038 035 033 041 039 036 034 Giong * MD Dependent Variable:KL 95% Confidence Interval Giong MD 300 400 Mean Std Error Lower Bound Upper Bound 80 039 000 039 040 120 037 000 037 038 160 034 000 034 035 200 80 120 160 200 031 041 040 037 035 000 000 000 000 000 031 041 039 037 034 032 042 040 038 035 Oneway Descriptives KL 95% Confidence Interval for Mean Std Mean Deviation Std Error N 3 3 3 039467 037333 034367 031367 041400 039700 037300 0010786 0005859 0005508 0004509 0006245 0004583 0003606 0006227 0003383 0003180 0002603 0003606 0002646 0002082 Lower Bound 036787 035878 032999 030247 039849 038562 036404 Upper Bound Minimum Maximum 042146 038789 035735 032487 042951 040838 038196 0387 0369 0340 0309 0409 0393 0369 0407 0380 0350 0318 0421 0402 0376 Tota l 034700 0003464 0002000 033839 035561 0345 0351 24 036954 0032057 0006544 035601 038308 0309 0421 Test of Homogeneity of Variances KL Levene Statistic df1 df2 1.825 Sig 16 151 ANOVA KL Sum of Squares Between Groups Within Groups Total Mean Square df F 000 000 000 000 16 23 000 Sig 92.054 000 Post Hoc Tests Homogeneous Subsets KL Duncan Subset for alpha = 0.05 NT Sig N 031367 034367 034700 037300 037333 039467 039700 041400 1.000 505 946 639 1.000 1.3 Ảnh hưởng tỷ lệ sống Univariate Analysis of Variance Between-Subjects Factors Value Label N Giong 300 12 2 400 80 120 160 200 12 6 6 MD Descriptive Statistics Dependent Variable:TLS Std Deviation Giong MD Mean 300 80 91.9167 2.02073 120 89.6667 2.12890 160 87.9600 2.66923 200 84.9667 1.53080 88.6275 86.0833 84.0533 81.7100 78.3333 82.5450 89.0000 86.8600 84.8350 81.6500 85.5862 3.21084 2.50416 4.32835 3.28288 1.33167 3.98637 3.78814 4.33124 4.34507 3.85318 4.70979 12 3 3 12 6 6 24 Total 400 80 120 160 200 Total Total 80 120 160 200 Total N Tests of Between-Subjects Effects Dependent Variable:TLS Source Type III Sum of Squares Mean Square df F Sig Corrected 398.909a 56.987 8.194 Model Intercept 175800.149 175800.149 2.528E4 Giong 221.981 221.981 31.917 MD 176.007 58.669 8.436 Giong * MD 921 307 044 Error 111.280 16 6.955 Total 176310.337 24 Corrected 510.188 23 Total a R Squared = 782 (Adjusted R Squared = 686) Estimated Marginal Means Grand Mean Dependent Variable:TLS 95% Confidence Interval Mean Std Error 85.586 Lower Bound 538 84.445 Upper Bound 86.727 Giong Dependent Variable:TLS 95% Confidence Interval Giong Mean 300 400 88.628 82.545 Std Error Lower Bound 761 761 87.014 80.931 Upper Bound 90.241 84.159 MD Dependent Variable:TLS 95% Confidence Interval MD Mean 80 120 89.000 86.860 Std Error 1.077 1.077 Lower Bound 86.718 84.578 Upper Bound 91.282 89.142 000 000 000 001 987 160 200 84.835 81.650 1.077 1.077 82.553 79.368 87.117 83.932 Giong * MD Dependent Variable:TLS 95% Confidence Interval Std Error Lower Bound Upper Bound Giong MD Mean 300 80 91.917 1.523 88.689 95.144 120 89.667 1.523 86.439 92.894 160 87.960 1.523 84.732 91.188 200 80 120 160 200 84.967 86.083 84.053 81.710 78.333 1.523 1.523 1.523 1.523 1.523 81.739 82.856 80.826 78.482 75.106 88.194 89.311 87.281 84.938 81.561 400 Oneway Descriptives TLS 95% Confidence Interval for Mean N Std Mean Deviation Std Error 91.9167 2.02073 1.16667 89.6667 2.12890 1.22912 3 87.9600 2.66923 1.54108 84.9667 1.53080 88380 86.0833 2.50416 1.44578 84.0533 4.32835 2.49898 81.7100 3.28288 1.89537 78.3333 1.33167 76884 Tota 24 85.5862 4.70979 96138 l Test of Homogeneity of Variances TLS Lower Bound 86.8969 84.3782 81.3293 81.1640 79.8626 73.3011 73.5549 75.0253 Upper Bound Minimum Maximum 96.9364 94.9552 94.5907 88.7694 92.3040 94.8056 89.8651 81.6414 89.75 87.83 86.00 83.80 83.50 79.33 79.75 77.20 93.75 92.00 91.00 86.70 88.50 87.83 85.50 79.80 83.5975 87.5750 77.20 93.75 Levene Statistic df1 1.215 df2 Sig 16 350 ANOVA TLS Sum of Squares Between Groups Within Groups Total Mean Square df 398.909 56.987 111.280 510.188 16 23 6.955 F Sig 8.194 000 Post Hoc Tests Homogeneous Subsets TLS Duncan NT Sig N Subset for alpha = 0.05 78.3333 81.7100 81.7100 84.0533 84.9667 86.0833 3 136 079 84.0533 84.9667 84.9667 86.0833 86.0833 87.9600 87.9600 87.9600 89.6667 89.6667 91.9167 113 060 100 Một số hình ảnh trình thí nghiệm Vùng ni ngao Quần Mục Chịi canh ngao Chuẩn bị bố trí thí nghiệm ... tốt nghiệp ? ?Nghiên cứu ảnh hưởng mật độ kích cỡ thả ni đến sinh trưởng, tỷ lệ sống nghêu (Meretrix lyrata Sowerby, 1851) giai đoạn nuôi thương phẩm vùng bãi triều, Kiến Thụy, Hải Phòng? ?? Mục tiêu... hưởng mật độ kích cỡ thả ni đến sinh trưởng, tỷ lệ sống nghêu (Meretrix lyrata Sowerby, 1851) giai đoạn nuôi thương phẩm vùng bãi triều, Kiến Thụy, Hải Phòng Thí nghiệm mật độ kích cỡ: yếu... VÀ ĐÀO TẠO TRƯỜNG ĐẠI HỌC NHA TRANG NGÔ XUÂN BA NGHIÊN CỨU ẢNH HƯỞNG CỦA MẬT ĐỘ VÀ KÍCH CỠ THẢ NI ĐẾN SINH TRƯỞNG, TỶ LỆ SỐNG CỦA NGHÊU (Meretrix lyrata SOWERBY, 1851) GIAI ĐOẠN NUÔI THƯƠNG PHẨM

Ngày đăng: 02/06/2021, 22:15

Tài liệu cùng người dùng

Tài liệu liên quan