Phƣơng pháp nghiên cứu

Một phần của tài liệu (LUẬN văn THẠC sĩ) nghiên cứu đa dạng thành phần loài thân mềm hai mảnh vỏ thuộc phân lớp pteriomorphia trong hệ sinh thái rừng ngập mặn tại khu vực huyện tiên yên, tỉnh quảng ninh (Trang 36)

3. Nội dung nghiên cứu

2.2. Phƣơng pháp nghiên cứu

2.2.1. hương pháp thu thập s liệu

Thu thập tài liệu liên quan tới vấn đề nghiên cứu

+ Thu thập số liệu về điều kiện tự nhiên, kinh tế xã hội... của huyện Tiên Yên từ UBND huyện Tiên Yên, tỉnh Quảng Ninh.

+ Điều tra, khảo sát ngƣời dân sống xung quanh khu vực RNM nghiên cứu.

2.2.2. hương pháp nghi n u tài iệu

Phƣơng pháp nghiên cứu tài liệu: Từ các công trình nghiên cứu, sách báo tài liệu xuất bản về lớp thân mềm hai mảnh vỏ. Kế thừa các thông tin về điều kiện tự nhiên xã hội và các công trình đã từng nghiên cứu về lớp thân mềm hai mảnh vỏ nói chung và phân lớp Pteriomorphia nói riêng.

Phƣơng pháp tham vấn chuyên gia: Tham vấn những chuyên gia có những hiểu biết nhất định tại địa điểm nghiên cứu để có thể tìm hiểu, đánh giá khách quan đƣợc về điều kiện tự nhiên, xã hội, môi trƣờng ở khu vực nghiên cứu, và những

chuyên gia đã từng nghiên cứu về lớp thân mềm hai mảnh vỏ. Tìm hiểu và xin ý kiến những ngƣời đã khảo sát trƣớc, cụ thể:

- Tham khảo về cách thức thu mẫu, địa điểm thu mẫu dƣới sự giúp đỡ và hƣớng dẫn của PGS.TS. Hoàng Ngọc Khắc và các cán bộ thuộc UBND thị trấn Tiên Yên, huyện Tiên Yên, tỉnh Quảng Ninh.

- Thực hiện thu mẫu đợt 2 và xin ý kiến dƣới sự giúp đỡ và hƣớng dẫn của NCS. Phạm Thanh Bình để tìm hiểu them thông tin về RNM và các loài thân mềm hai mảnh vỏ phân lớp Pteriomorphia tại KVNC.

2.2.3. hương pháp thu mẫu ngoài th a

- Vật dụng, thiết bị:

Dụng cụ và thiết bị phục vụ dùng trong quá trình thu mẫu, định loại và phân tích bao gồm: xẻng nhỏ xúc đất, túi nylon (túi polyethylene) đựng mẫu, máy định vị GPS, túi bấm, thùng xốp, xô nhựa, giấy ghi nhãn, thƣớc dây, máy ảnh, bút, sổ và dụng cụ phân tích mẫu gồm kính lúp, thƣớc palme, panh, hộp nhựa nhỏ có nắp kín để lƣu trữ mẫu. Hóa chất dùng trong nghiên cứu gồm cồn 900

để bảo quản mẫu. - Vị trí lấy mẫu: Theo sơ đồ hình 2.1 và bảng 2.1

a. Phƣơng pháp thu mẫu

- Thu mẫu định lƣợng:

+ Mẫu định lƣợng đƣợc tính trên diện tích 1 m². Tổng số ô định lƣợng là 5 ô ở mỗi sinh cảnh. Các mẫu thu đƣợc đƣợc để trong các lọ nhựa, túi nhựa theo kích thƣớc. Giá trị của mẫu định lƣợng cho biết mật độ, sự phong phú về số lƣợng hoặc sự đa dạng về thành phần loài của KVNC.

+ Sau khi xác định đƣợc vị trí cần thu mẫu, dùng thƣớc dây xác định ô tiêu chuẩn theo diện tích ở trên, thu tất cả các mẫu có trong ô đó. Mẫu thu đƣợc ở mỗi ô vuông cho vào một túi nilon có đề nhãn. Nhãn ghi các thông tin: Địa điểm, thời gian, tọa độ, sinh cảnh, đặc điểm thảm thực vật… hoặc các lƣu ý cần thiết khác. - Thu mẫu định tính: Tiến hành thu mẫu tại các địa điểm khác nhau nhằm mục đích thu thành phần các loài thân mềm hai mảnh vỏ thuộc phân lớp Pteriomorphia.

- Thu lƣợm động vật ở các bãi triều khi thủy triều xuống.

- Mua mẫu tại các chợ, thuyền của ngƣời dân sau khi đi đánh bắt hải sản về. - Ký hiệu các điểm thu mẫu

Bảng 2. 1. Địa điểm, tọa độ các xã lấy mẫu tại KVNC

STT Địa điểm Kí hiệu viết tắt Tọa độ Đặc điểm sinh cảnh 1 Đông Hải ĐH 21°18'48.97"Bắc 107°30'55.03"Đông Rừng tự nhiên, nền đáy, đất cát bùn, phù sa, đá sỏi… 2 Đồng Rui ĐR 21°14'2.14"Bắc 107°23'29.40"Đông Rừng tự nhiên gần cửa sông, đất cát bùn 3 Hải Lạng HL 21°16'20.67"Bắc 107°22'24.56"Đông Rừng trồng, rừng tái sinh, nhiều các giá thể vô sinh 4 Thác Bƣởi TB 21°18'35.57"Bắc 107°26'1.25"Đông Rừng tự nhiên 5 Thôn Thƣợng TT 21°14'36.76"Bắc 107°23'54.06"Đông Rừng tự nhiên chiếm diện tích lớn 6 Thủy Cơ TC 21°19'5.54"Bắc 107°25'51.69"Đông Rừng trồng, nền đáy bùn, cát 7 Tiên Lãng TL 21°18'57.15"Bắc 107°25'30.45"Đông Rừng trồng, nền đáy bùn, cát 8 Thôn 4 T4 21°14'37.90"Bắc 107°26'10.82"Đông Rừng tự nhiên chiếm diện tích lớn 2.2.4. Phương pháp xử lý mẫu - Rửa sạch mẫu.

- Luộc mẫu thu đƣợc sau đó tiến hành tách thân mềm ra khỏi vỏ để tiến hành phân tích mẫu.

2.2.5. Phương pháp phân tích mẫu

- Mô tả mẫu theo thứ tự bao gồm:

+ Tên loài: Tài liệu công bố gốc. + Synonym: Tên gọi khác (nếu có).

+ Mẫu vật: Phân tích để định loại.

+ Phân bố: KVNC và khu vực khác ở Việt Nam.

+ Nhận xét: Rút ra một số nhận xét: nơi đã thu đƣợc loài, kích thƣớc các cá thể của loài so với mô tả gốc các cá thể cùng loài thu đƣợc ở các khu vực lân cận…

- Định loại xác định tên loài: Dựa vào hình thái thái vỏ, các đặc điểm kích thƣớc, màu sắc, các đặc điểm khác biệt.

2.2.6 Phương pháp xử lý số liệu

Dữ liệu mẫu vật thu đƣợc tính toán bằng các công cụ để xác định mật độ cá thể. Các chỉ số đa dạng sinh học của loài thân mềm hai mảnh vỏ gồm:

+ Mật độ (Số cá thể /m2):

Mật độ: V (Số cá thể /m2

) = Σn ΣS

Trong đó:

Σn: Là tổng số cá thể trong các ô nghiên cứu.

ΣS: Là tổng diện tích các ô nghiên cứu.

+ Độ phong phú của loài (P%): Đƣợc tính theo công thức của Kreds (1989)

P% = Σn

ΣS x 100%

+ Độ đa dạng của loài (H’): Đƣợc tính theo chỉ số Shannon – Weiner

H’= ∑

Trong đó, H’ : độ đa dạng của loài

ni: số lƣợng cá thể của loài thứ i N: Tổng số cá thể của tất cả các loài

+Tần số xuất hiện (độ thƣờng gặp): đƣợc tính bằng công thức của Sharma (2003)

C’ = Pp x 100%

Trong đó:

C’: Là tần số xuất hiện (độ thƣờng gặp).

p: Là số lƣợng các địa điểm thu mẫu có loài xuất hiện.

P: Là tổng số các địa điểm thu mẫu khi nghiên cứu.

Đánh giá tần số xuất hiện theo giá trị của C’:

Loài thƣờng gặp C’> 50%,

Loài ít gặp 25% ≤ C’≤ 50%, Loài ngẫu nhiên C’< 25%.

+ Độ gần gũi về thành phần loài:

Để phân tích mối quan hệ về thành phần loài thân mềm hai mảnh vỏ giữa KVNC và khu vực khác, tôi sử dụng chỉ số tƣơng đồng S (Sorensen, 1948). Chỉ số này đƣợc tính theo công thức sau:

S = 2C A+B

Trong đó: S là chỉ số tƣơng đồng; A, B là tổng số loài của hai khu hệ cần so sánh; C là số loài trùng nhau của hai khu hệ.

CHƢƠNG 3: KẾT QUẢ NGHIÊN CỨU VÀ THẢO LUẬN 3.1. Danh lục các loài thuộc phân lớp Pteriomorphia trong KVNC.

Thành phần loài Thân mềm hai mảnh vỏ thuộc phân lớp Pteriomorphia ở vùng rừng ngập mặn ven biển ngập triều, qua quá trình thu mẫu và phân loại mẫu sau 2 đợt đã thu đƣợc 96 mẫu, trong đó có 6 mẫu định tính và 90 mẫu định lƣợng. Tổng số thu đƣợc 1839 cá thể. Bƣớc đầu đã phát hiện ở KVNC có 23 loài thuộc 12 giống: Anadara, Estellacar, Crassotrea, Saccostrea, Anomia, Placuna, Brachidontes, Xenostrobus, Perna, Modiolus, Septifer Isognomon. Đƣợc xếp vào 7 họ; Arcidae, Noetiidae, Ostreidae, Anomiidae, Placunidae, Mytilidae, Isognomonidae và 5 bộ; Arcoida, Mytiloida, Ostreoida, Pectinoida và Pterioida. Cấu trúc thành phần loài thuộc phân lớp Pteriomorphia trong bảng 3.1.

Bảng 3. 1. Danh lục thành phần loài thuộc phân lớp Pteriomorphia tại KVNC

STT Tên khoa học Tên Việt Nam Số lƣợng cá thể Mẫu Mẫu ĐL ĐT Bộ - ARCOIDA Họ - Arcidae Họ sò

1 Anadara antiquata (Linnaeus, 1758) Sò quéo 58 + + 2 Anadara granosa (Linnaeus, 1758) Sò huyết 110 + + 3 Anadara subcrenata (Lienschke, 1869) Sò lông 18 + + 4 Anadara nodifera (Martens, 1860) Sò nứa 14 +

Họ - Noetiidae

5 Estellacar olivacea (Reeve, 1844) Sò hình thang 25 + +

Bộ - MYTILOIDA

Họ - Mytilidae Họ vẹm

6 Brachidontes curvatus, Dunker, 1857 Quéo 143 + + 7 Brachidontes emarginatus (Reeve, 1858) Quéo 186 + + 8 Brachidontes senhousei (Benson, 1842) Quéo 66 + + 9 Xenostrobus atrata (Lischke, 1871) Dòm đen 308 + + 10 Perna viridis (Linnaeus, 1758) Vẹm vỏ xanh 51 + + 11 Modiolus philippinarum (Hanley, 1843) Quéo 55 + +

12 Septifer virgatus (Wiegmann, 1837) Quéo 77 + +

Bộ - OSTREOIDA

Họ - Ostreidae Họ Hàu

13 Crassostrea ariakensis (Fujita, 1913) Hàu 125 + + 14 Crassotrea gigas (Thunberg, 1793) Hàu 110 + + 15 Crassotrea rivularis (Gould,1864) Hàu cửa sông 162 + + 16 Saccostrea cucullata (Born, 1778) Hàu lá 50 + + 17 Saccostrea glomerata (Gould, 1850) Hàu tròn 29 + + 18 Saccostrea mordax (Gould, 1850) Hàu 60 + + 19 Saccostrea pestigris (Hanley, 1846) Hàu 20 + +

Bộ - PECTINOIDA

Họ - Anomiidae Họ điệp cánh

20 Anomia aenigmatica (Holten, 1803) Điệp lá 31 + +

21 Anomia cyteum (Gray, 1850) 73 + +

Họ - Placunidae

22 Placuna placenta (Linnaeus, 1758) Điệp giấy 62 + +

Bộ - PTERIOIDA

Họ - Isognomonidae

23 Isognomon ephippium (Linnaeus, 1758) 6 +

Tổng 1839

Ghi hú: ĐL = Đ nh ượng, ĐT: = Đ nh tính.

Trong hệ sinh thái rừng ngập mặn cửa sông ven biển có số lƣợng loài, mức độ đa dạng của các giống và họ trong bảng trên cho thấy số lƣợng loài ở có sự khác nhau rõ rệt. Họ Ostreidae và Mytilidae chiếm số lƣợng lớn với 7 loài, Arcidae có 4 loài, họ Anomiidae có 2 loài, họ Noetiidae, Placunidae và Isognomonidae có 1 loài tại KVNC.

3.2. Cấu trúc thành phần loài phân lớp Pteriomorphia trong KVNC

Về bậc bộ:

Trong hệ sinh thái rừng ngập mặn KVNC có 5 bộ là: Arcoida, Ostreoida, Mytiloida, Pectinoida và Pterioida. Trong đó bộ Mytiloida có 1 họ (chiếm 14,29% tổng số họ), 5 giống (chiếm 41,67% tổng số giống) và 7 loài (chiếm 30,44% tổng số loài); Bộ Ostreoida có 1 họ (chiếm 14,29% tổng số họ), 2 giống (chiếm 16,67% tổng số giống) và 7 loài (chiếm 30,44% tổng số loài); Bộ Arcoida có 2 họ (chiếm

28,57% tổng số họ), 2 giống (chiếm 16,67% tổng số giống) và 5 loài (chiếm 21,74% tổng số loài); Bộ Pectinoida có 2 họ (chiếm 28,57% tổng số họ), 2 giống (chiếm 16,67% tổng số giống) và 3 loài (chiếm 13,04% tổng số loài); Bộ Pterioida có 1 họ (chiếm 14,29% tổng số họ), 1 giống (chiếm 8,33% tổng số giống) và 1 loài (chiếm 4,35% tổng số loài).

Hình 3. 1.Đa dạng về các bộ tại KVNC

Xét về số lƣợng cá thể của các bộ, bộ Mytiloida chiếm tổng số cá thể nhiều nhất với 932 cá thể chiếm 50,68% tổng số cá thể thu đƣợc. Bộ Ostreoida có 556 cá thể chiếm 19,58% tổng số cá thể thu đƣợc. Bộ Arcoida có 225 cá thể chiếm 12,24% tổng số cá thể thu đƣợc. Bộ Pectinoida chiếm có 166 cá thể chiếm 9,03% tổng số cá thể thu đƣợc. Bộ Pterioida có 6 cá thể chiếm 0,33% tổng số cá thể thu đƣợc.

0 1 2 3 4 5 6 7 8

Arcoida Mytiloida Ostreoida Pectinoida Pterioida

Loài Giống Họ 50,68% 30,23% 12,24% 9,03% 0,33% Mytiloida Ostreoida Arcoida Pectinoida Pterioida

Về bậc họ:

Hệ sinh thái rừng ngập mặn tại KVNC có tổng số 7 họ phát hiện đƣợc ở khu vực nghiên cứu thì họ Mytilidae có số giống nhiều nhất là 5 giống chiếm 41,67% tổng số giống. Họ Ostreidae là 2 giống, chiếm 16,67%. Các họ còn lại chỉ có 1 giống, chiếm 8,33% tổng số giống.

Hình 3. 3. Tỉ lệ % các họ thuộc phân lớp Pteriomorphia tại KVNC

Xét theo số lƣợng loài của các họ, số lƣợng loài phát hiện đƣợc ở khu vực nghiên cứu khá phong phú về thành phần loài 23 loài, trong đó họ Mytilidae (họ Vẹm) và Ostreidae (họ Hàu) có số lƣợng loài nhƣ nhau với 7 loài, chiếm 30,43% tổng số loài. Họ Arcidae có 4 loài, chiếm 17,39% tổng số loài. Họ Anomiidae có 2 loài, chiếm 8,7% tổng số loài. Họ Isognomonidae, Noetiidae và Placunidae có 1 loài, chiếm 4,35% tổng số loài.

Bảng 3. 2. Tỉ lệ % số cá thể, giống trong các họ phân lớp Pteriomorphia tại KVNC

STT Họ Cá thể Loài Giống N % n % n % 1 Arcidae 200 10,88 4 17,39 1 8,33 2 Noetiidae 25 1,36 1 4,35 1 8,33 3 Mytilidae 886 48,18 7 30,43 5 41,67 4 Ostreidae 556 30,23 7 30,43 2 16,67 5 Anomiidae 104 5,66 2 8,70 1 8,33 6 Placunidae 62 3,37 1 4,35 1 8,33 7 Isognomonidae 6 0,33 1 4,35 1 8,33 Tổng 1839 100,00 23 100,00 12 100,00 31% 31% 17% 9% 4% 4% 4% Mytilidae Ostreidae Arcidae Anomiidae Isognomoni dae Placunidae

Hình 3. 4. Số lƣợng cá thể tại các họ thuộc phân lớp Pteriomorphia tại KVNC

Xét về số lƣợng cá thể của từng họ Mytilidae là họ chiếm ƣu thế với 886 cá thể, chiếm 48,18%. Họ Ostreidae với 556 cá thể, chiếm 30,23% tổng số cá thể. Họ Arcidae có 200 cá thể, chiếm 10,88%. Họ Anomiidae có 104 cá thể, chiếm 5,66%. Họ Placunidae chỉ có 62 cá thể, chiếm tỷ lệ với 3,37%. Họ Noetiidae có 18 cá thể, chiếm 0,95%. Họ Isognomonidae chỉ có 6 cá thể, chiếm tỷ lệ thấp nhất với 0,33%.

Về bậc loài và phân loài:

- Về độ đa dạng loài:

Qua qua quá trình thu mẫu 2 đợt và phân tích đã xác định đƣợc 23 loài thân mềm hai mảnh vỏ thuộc phân lớp Pteriomorphia, chứng tỏ khu vực này có độ đa dạng khá cao (H = 0,922). .000% 10.000% 20.000% 30.000% 40.000% 50.000% 60.000%

- Về độ số lƣợng loài:

Với tổng số 23 loài trong KVNC cho thấy khu hệ thân mềm 2 mảnh vỏ KVNC có số lƣợng loài tƣơng đối cao. Kết quả quá trình nghiên cứu đã bổ sung cho khu vực rừng ngập mặn Tiên Yên, tỉnh Quảng Ninh 23 loài thân mềm hai mảnh vỏ thuộc phân lớp Pteriomorphia. Thân mềm hai mảnh vỏ ở khu vực nghiên cứu không chỉ đa dạng về thành phần loài mà còn đa dạng về kích thƣớc và hình thái.

- Về kích thƣớc:

Thành phần loài của thân mềm hai mảnh vỏ thuộc phân lớp Pteriomorphia tại KVNC bao gồm đầy đủ nhóm có kích thƣớc từ bé đến rất lớn:

+ Nhóm có kích thƣớc lớn (lớn hơn 100mm): Gồm các loài 01 loài thuộc họ Mytilidae: Perna viridis. 7 loài thuộc họ Ostreidae: Crassostrea ariakensis,

Crassotrea gigas, Crassotrea rivularis, Saccostrea cucullata, Saccostrea glomerata, Saccostrea mordax, Saccostrea pestigris. 1 loài thuộc họ Isognomonidae: Isognomon ephippium.

+ Nhóm có kích thƣớc trung bình (từ 50mm - 100mm): Gồm 04 loài thuộc họ Arcidae gồm: Anadara antiquata, Anadara granosa, Anadara subcrenata, Anadara nodifera. 01 loài thuộc họ Mytilidae: Xenostrobus atrata. 02 loài thuộc họ Anomiidae: Anomia aenigmatica, Anomia cyteum. 01 loài thuộc họ Placunidae:

Placuna placenta.

+ Nhóm có kích thƣớc nhỏ (dƣới 50mm): Gồm các 01 loài thuộc họ Noetiidae:

Estellacar olivacea. 05 loài thuộc họ Mytilidae: Brachidontes curvatus, Brachidontes emarginatus, Brachidontes senhousei, Modiolus philippinarum, Septifer virgatus.

- Về mật độ cá thể của loài:

Mật độ trung bình của thân mềm hai mảnh vỏ thuộc phân lớp Pteriomorphia trong 90 mẫu ĐL thuộc các sinh cảnh trong KVNC ở mức khá cao (V = 12,27 con/m2). Trong đó Xenostrobus atrata là loài chiếm ƣu thế nhất với 2,48 con/m2, kế đến là Brachidontes emarginatus với 1,2 con/m2, tiếp theo là Brachidontes curvatus

Crassostrea ariakensis với 1,19 con/m2, Crassotrea rivularis với 1,06 con/m2. Các loài còn lại xuất hiện với mật độ thấp từ dƣới 1 con/m2, góp phần làm phong phú thêm về thành phần loài Thân mềm hai mảnh vỏ ở khu vực nghiên cứu.

Mật độ trung bình của các loài ở phân lớp Pteriomorphia trong KVNC trên mỗi ô nghiên cứu (1m2) cũng có sự khác nhau. Trong mỗi ô nghiên cứu thì mật độ của loài Crassostrea ariakensis ở ô nghiên cứu TT- 1 là cao nhất có v = 55 con/m2, tiếp theo là Brachidontes emarginatus với v = 33 con/m2 ở ô nghiên cứu TT-15,

Xenostrobus atrata với v = 32 con/m2 ở ô nghiên cứu TT-9. - Về độ phong phú:

Bảng 3. 3. Độ phong phú của các loài thuộc phân lớp Pteriomorphia tại KVNC

STT Tên khoa học Độ phong phú (P%) Ghi

chú

Họ - Arcidae ĐH ĐR HL TL KVNC

1 Anadara antiquata (Linnaeus, 1758) 0,09 3,53 3,62

2 Anadara granosa (Linnaeus, 1758) 0,27 0,63 0,18 4,44 5,62

3 Anadara subcrenata (Lienschke, 1869) 0,09 0,27 0,36

4 Anadara nodifera (Martens, 1860)

Họ - Noetiidae

5 Estellacar olivacea (Reeve, 1844) 0,09 0,63 0,72

Họ - Mytilidae

6 Brachidontes curvatus, Dunker, 1857 1,90 3,17 3,99 0,63 9,69

7 Brachidontes emarginatus (Reeve, 1858) 4,98 4,80 9,78

8 Brachidontes senhousei (Benson, 1842) 0,45 1,27 1,72

9 Xenostrobus atrata (Lischke, 1871) 6,07 9,51 2,08 2,54 20,20

10 Perna viridis (Linnaeus, 1758) 0,09 1,18 0,45 1,72

11 Modiolus philippinarum (Hanley, 1843) 0,54 0,36 0,63 1,27 2,81

12 Septifer virgatus (Wiegmann, 1837) 0,18 0,45 1,45 2,08

Họ - Ostreidae

13 Crassostrea ariakensis (Fujita, 1913) 0,27 7,25 0,82 1,36 9,69

14 Crassotrea gigas (Thunberg, 1793) 1,45 1,09 3,08 5,89

15 Crassotrea rivularis (Gould,1864) 3,71 3,44 0,72 0,72 8,61

16 Saccostrea cucullata (Born, 1778) 2,99 0,27 0,09 3,35

17 Saccostrea glomerata (Gould, 1850) 0,27 1,90 0,18 0,27 2,63

18 Saccostrea mordax (Gould, 1850) 0,82 1,99 0,63 3,44

Họ - Anomiidae

20 Anomia aenigmatica (Holten, 1803) 0,18 0,54 0,45 0,54 1,72

21 Anomia cyteum (Gray, 1850) 1,18 1,45 1,27 0,72 4,62

Họ -Placunidae

22 Placuna placenta (Linnaeus, 1758) 0,18 0,36 0,27 0,18 1,00

Họ - Isognomonidae

23 Isognomon ephippium (Linnaeus, 1758) 0,18 0,18

Tổng 100%

Ghi hú: ĐH = Đông H i, ĐR = Đồng Rui, HL = H i Lạng, TL = Ti n Lãng.

Nhận xét:

+ Ở Đông Hải: Loài có số lƣợng cá thể nhiều nhất là Xenostrobus atrata với

Một phần của tài liệu (LUẬN văn THẠC sĩ) nghiên cứu đa dạng thành phần loài thân mềm hai mảnh vỏ thuộc phân lớp pteriomorphia trong hệ sinh thái rừng ngập mặn tại khu vực huyện tiên yên, tỉnh quảng ninh (Trang 36)

Tải bản đầy đủ (PDF)

(106 trang)