1. Trang chủ
  2. » Luận Văn - Báo Cáo

NGHIÊN CỨU TÁC ĐỘNG CỦA BIẾN ĐỔI KHÍ HẬU ĐẾN CỰC TRỊ DÒNG CHẢY TRÊN LƯU VỰC SÔNG NHUỆ ĐÁY THUỘC THÀNH PHỐ HÀ NỘI

85 710 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

Thông tin cơ bản

Định dạng
Số trang 85
Dung lượng 5,9 MB

Nội dung

ĐẠI HỌC QUỐC GIA HÀ NỘI TRƯỜNG ĐẠI HỌC KHOA HỌC TỰ NHIÊN Nguyễn Ý Như NGHIÊN CỨU TÁC ĐỘNG CỦA BIẾN ĐỔI KHÍ HẬU ĐẾN CỰC TRỊ DÒNG CHẢY TRÊN LƯU VỰC SÔNG NHUỆ ĐÁY THUỘC THÀNH PHỐ HÀ NỘI Chuyên ngành: Thủy văn học Mã số: 60 44 90 LUẬN VĂN THẠC SĨ KHOA HỌC NGƯỜI HƯỚNG DẪN KHOA HỌC: PGS.TS Nguyễn Thanh Sơn Hà Nội - 2011 LỜI CẢM ƠN Luận văn hoàn thành Khoa Khí tượng Thủy văn Hải dương học, Trường Đại học Khoa học Tự nhiên, Đại học Quốc gia Hà Nội hướng dẫn khoa học PGS TS Nguyễn Thanh Sơn Tác giả xin bày tỏ lòng biết ơn chân thành tới thầy, người đã hết lòng động viên, tận tình giúp đỡ quan tâm tới bước nghiên cứu học viên Để thực luận văn, tác giả hỗ trợ mặt tài đề tài cấp Đại học Quốc gia mã số QGTD.10.06, giúp đỡ thời gian, điều kiện nghiên cứu thuận lợi từ thầy cô Bộ môn Thủy văn, thầy cô giáo, đồng nghiệp khoa bè bạn trình học tập nghiên cứu Nguyễn Ý Như MỤC LỤC LỜI CẢM ƠN .2 MỤC LỤC DANH MỤC HÌNH DANH MỤC BẢNG DANH MỤC CÁC KÝ HIỆU VIẾT TẮT MỞ ĐẦU 10 Chương 11 TỔNG QUAN 11 1.1 TÌNH HÌNH NGHIÊN CỨU TÀI NGUYÊN NƯỚC TRÊN THẾ GIỚI LIÊN QUAN TỚI BIẾN ĐỔI KHÍ HẬU 11 1.2 TÌNH HÌNH NGHIÊN CỨU TRONG NƯỚC 14 1.3 MỘT SỐ THẢO LUẬN 15 1.4 ĐIỀU KIỆN ĐỊA LÝ TỰ NHIÊN – KINH TẾ XÃ HỘI LƯU VỰC SÔNG NHUỆ ĐÁY THUỘC ĐỊA PHẬN THÀNH PHỐ HÀ NỘI 16 1.4.1 Điều kiện địa lý tự nhiên 16 1.4.2 Đặc điểm kinh tế - xã hội .21 1.4.3 Hiện trạng tài nguyên nước 22 1.4.4 Một số nghiên cứu tiêu biểu lưu vực sông Nhuệ - Đáy .23 Chương 25 LỰA CHỌN KỊCH BẢN BIẾN ĐỔI KHÍ HẬU VÀ MÔ HÌNH MÔ PHỎNG KHÍ HẬU – DÒNG CHẢY 25 2.1 KỊCH BẢN BIỂN ĐỔI KHÍ HẬU CHO KHU VỰC NGHIÊN CỨU 25 2.1.1 Lựa chọn kịch biến đổi khí hậu .25 2.1.2 Một số công cụ sử dụng xây dựng kịch .26 2.1.3 Phương pháp nội suy 27 2.1.4 Hệ số chỉnh sai .28 2.2 CÁC MÔ HÌNH KHÍ HẬU – DÒNG CHẢY 30 2.2.1 Giới thiệu số mô hình khí hậu – dòng chảy 30 2.2.2 Nhận xét lựa chọn mô hình .36 2.3 MÔ HÌNH THỦY VĂN CHO LƯU VỰC NGHIÊN CỨU 37 2.3.1 Cấu trúc mô hình NAM 37 2.3.2 Các yếu tố ảnh hưởng đến dòng chảy mô hình NAM .39 2.3.3 Các thông số mô hình NAM 42 2.3.4 Điều kiện ban đầu mô hình 43 2.3.5 Hàm mục tiêu .43 Chương 44 ĐÁNH GIÁ BIẾN ĐỘNG CỰC TRỊ DÒNG CHẢY DƯỚI TÁC ĐỘNG CỦA BIẾN ĐỔI KHÍ HẬU .44 3.1 CƠ SỞ DỮ LIỆU 44 3.1.1 Số liệu đầu vào mô hình NAM 44 3.1.2 Số liệu sử dụng đánh giá biến đổi 45 3.2 ÁP DỤNG MÔ HÌNH CHO KHU VỰC NGHIÊN CỨU 47 3.2.1 Chỉ tiêu đánh giá hoạt động mô hình 48 3.2.2 Hiệu chỉnh kiểm định mô hình mưa dòng chảy NAM 49 3.3 ĐÁNH GIÁ BIẾN ĐỘNG CỰC TRỊ DÒNG CHẢY 51 3.3.1 Biến động đặc trưng dòng chảy lũ 51 3.3.2 Biến động đặc trưng dòng chảy kiệt 64 .76 .76 KẾT LUẬN .78 TÀI LIỆU THAM KHẢO 82 DANH MỤC HÌNH DANH MỤC BẢNG DANH MỤC CÁC KÝ HIỆU VIẾT TẮT A1B, A2 Các kịch phát thải khí nhà kính IPCC ANN Artificial Neural Network (Mạng thần kinh nhân tạo) AQUASIM A Computer Program for the Identification and Simulation of Aquatic Systems (Phần mềm sinh thái) BĐKH Biến đổi hậu BOD Biochemical oxygen demand (Nhu cầu oxy sinh hóa) COD Chemical oxygen demand (Nhu cầu oxy hóa học) DEM Digital Elevation Model (Mô hình độ cao số hóa) DHI Danish Hydraulic Institute (Viện nghiên cứu thủy lực Đan Mạch) ECHAM European Centre Hamburg Model (Mô hình khí hậu toàn cầu Trung tâm châu Âu Hamburg) ESRI Environmental Systems Research Institute (Viện nghiên cứu hệ thống môi trường) FDC Flow duration curve (Đường cong thời khoảng dòng chảy) GCM Global Climate Model (Mô hình khí hậu toàn cầu) HBV Hydrologiska Byråns Vattenbalansavdelning (Mô hình cân nước thủy văn) HEC-HMS Hydrologic Engineering Center’s Hydrologic Modelling System (Hệ thống mô hình hóa thủy văn trung tâm kỹ thuật thủy văn) HIRLAM High Resolution Limited Area Model (Mô hình khu vực hạn chế độ phân giải cao) Hydro-BEAM Hydrological River Basin Environment Assessment Model (Mô hình đánh giá môi trường lưu vực sông) IHMS Integrated Hydrological Modelling System hệ thống mô hình thủy văn kết hợp IPCC Intergovernmental Panel on Climate Change (Ban Liên phủ Biến đổi khí hậu) KNK Khí nhà kính KT – XH Kinh tế - Xã hội MIKE – SHE Système Hydrologique Européen (Mô hình hệ thống thủy văn Châu Âu) NAM Nedbør - Afstrømnings – Models (Mô hình mưa – dòng chảy) NASIM Niederschlag – Abfluss Simulation Model (Mô hình tính mưa – dòng chảy ND1 – ND5 Tiểu lưu vực thuộc Lưu vực sông Nhuệ - Đáy ND Nhuệ - Đáy NNK người khác (chỉ đồng tác giả công trình, báo …) NSE Nash–Sutcliffe efficiency (hệ số Nash–Sutcliffe) NWSRFS National Weather Service River Forecasting System (Hệ thống dự báo thời tiết quốc gia) OPYC Ocean General Circulation Model (Mô hình hoàn lưu chung đại dương) QUAL2E Enhanced Stream Water Quality Models (Mô hình chất lượng nước) RegCM REGional Climate Model (Mô hình khí hậu khu vực ICTP) SAC – SMA Sacramento Soil Moisture Accounting (Mô hình tính toán hàm lượng ẩm đất) SCS Soil Conservation Service (Phương pháp bảo tồn đất) SWAT Soil and Water Assessment Tool KH KTTV& MT Khoa học Khi tượng Thủy văn Môi trường WHO World Health Organization MỞ ĐẦU Lưu vực sông Nhuệ - sông Đáy lưu vực lớn, có vị trí địa lý đặc biệt, đóng vai trò quan trọng kinh tế nước nói chung, vùng đồng sông Hồng nói riêng Sông Nhuệ sông Đáy hai sông cung cấp nguồn nước quan trọng cho sản xuất nông nghiệp, công nghiệp dân sinh cho cộng đồng dân cư, chịu áp lực mạnh mẽ gia tăng dân số, trình đô thị hoá, hoạt động kinh tế xã hội (KT – XH) diễn lưu vực Trước yêu cầu lớn đặt với nguồn nước lưu vực để đáp ứng nhu cầu sử dụng nước ngày cao số lượng phát triển dân sinh KT - XH đòi hỏi nguồn nước để trì hệ sinh thái, hạn chế ô nhiễm nguồn nước, xu tài nguyên nước suy giảm chất lượng, tác động yếu tố tự nhiên hoạt động người tác động biến đổi khí hậu Đối với lưu vực vấn đề đánh giá tác động biến đổi khí hậu đến cực trị dòng chảy lưu vực sông Nhuệ Đáy thuộc địa bàn thành phố Hà Nội nhằm tạo tiền đề cho việc xây dựng giải pháp giảm thiểu tác hại biến đổi khí hậu, giúp nhà quản lý, nhà hoạch định sách xác định chiến lược phát triển kinh tế bền vững đảm bảo an sinh xã hội vấn đề có ý nghĩa khoa học thực tiễn Luận văn bố cục thành chương, mở đầu, kết luận, tài liệu tham khảo sau: • Chương Tổng quan • Chương Lựa chọn kịch biến đổi khí hậu mô hình mô khí hậu – dòng chảy • Chương Đánh giá biến động dòng chảy cực đoan tác động biến đổi khí hậu 10 giá trị dòng chảy nhỏ dòng chảy trung bình đồng thời tăng giá trị dòng chảy lớn trung bình, xu biến đổi diễn biến mạnh theo kịch A2 so với kịch A1B Dòng chảy vượt 90% thường sử dụng số hướng dẫn mà điều khiển công việc xử lý nước thải đánh giá cấp phép cho hoạt động bơm hút nước Sự giảm tần suất xuất Q90 nghĩa tiềm cho lưu lượng khai thác giảm Một điểm khác cần quan tâm tốc độ thay đổi cường độ vượt giá trị Q90 phân vị khác thay đổi không đáng kể tác động kịch bản, tác động kịch A2 gây biến đổi lớn Mặc dù xu hướng đối lập tìm thấy, xu hướng giảm cường độ vượt phân vị lớn 50 xu hướng ngược lại, tăng với phân vị nhỏ 50 (bảng 6) Bảng Thay đổi cường độ vượt Qx Thời đoạn Giai đoạn A1B A2 Phân vị 90 90 82 75 75 71 50 52 52 20 23 29 Kết cho thấy lưu vực nghiên cứu, lượng dòng chảy nhỏ yêu cầu để trì dòng chảy kiệt, nước ngầm để trì dòng chảy kiệt lưu vực nhỏ dẫn đến biến đổi dòng chảy kiệt lớn; đồng thời số cho thấy tác động biến đổi khí hậu không đáng kể dòng chảy kiệt 3.3.2.4 Tần suất dòng chảy kiệt Không giống với đường cong thời khoảng dòng chảy, thể tỉ lệ thời gian giá trị dòng chảy bị vượt quá, đường cong tần suất dòng chảy kiệt thể tỉ lệ năm dòng chảy bị vượt Các số tần suất dòng chảy kiệt sử dụng rộng rãi nghiên cứu hạn, thiết kế hệ thống cấp nước, tính toán lượng nước mặt khai thác đảm bảo an toàn, phân loại tiểm dòng chảy cho khả pha loãng, điều phối lượng chất thải vào dòng chảy, trì lưu lượng định sông việc phân tích cần thiết 71 Đường tần suất dòng chảy kiệt xây dựng sở chuỗi số liệu dòng chảy tháng nhỏ năm Một lý việc lựa chọn bước thời gian tháng Nhuệ Đáy nằm khu vực nhiệt đới ẩm gió mùa, dòng chảy năm không phản ánh trận kiệt khắc nghiệt nhất, mùa hè kiệt bất thường mùa mưa, dòng chảy trung bình năm rõ sai lệch so với bình thường với dòng chảy năm thời gian xuất kiệt xác định dẫn đến kết luận tác động bị giới hạn (Bonacci, 1993) [19], dòng chảy tháng khoảng thời gian thích hợp cho nghiên cứu hạn nông nghiệp, cấp nước nước ngầm chứa nhiều thông tin chi tiết năm, có khoảng thời gian đủ dài để loại tượng có ý nghĩa gọi kiệt cục Các hàm phân bố thường sử dụng tài liệu liên quan đến dòng chảy kiệt, phân bố Weibull, Gumbel, Pearson III, phân bố log chuẩn, kiểm tra để lựa chọn hàm phân bố lý thuyết phù hợp với chuỗi dòng chảy nhỏ năm với bước thời gian tháng lưu vực Nhuệ Đáy, với phương pháp đánh giá dựa vào kiểm tra đồ thị Kết cho thấy phân bố log – chuẩn nên sử dụng mô tả thống kê dòng chảy kiệt lưu vực Nhuệ Đáy cho thời kỳ kịch biến đổi khí hậu A1B, A2 (hình 14a, b, c) Kết cho thấy biến đổi tần suất dòng chảy tháng kiệt tương tự với dòng chảy tháng lớn xét phạm vi không gian, xu hướng biến đổi mạnh tác động kịch A2 Tuy nhiên, điểm khác biệt so với tần suất lũ cực trị thay đổi nhỏ (Bảng 7), diễn biến có khả giá trị cực trị dòng chảy kiệt thấp Cụ thể cường độ dòng chảy mùa kiệt với tần suất xuất lại 20 năm tăng 8.7% 17.4% tương ứng với kịch A1B, A2, kịch A1B A2 dẫn đến tăng 21.17%, 59.7% tính toán cường độ dòng chảy kiệt với p 95% 72 Hình 14a So sánh mức độ phù hợp đường tần suất lý luận đường thực nghiệm dòng chảy tháng kiệt – lưu vực ND1 – Giai đoạn Hình 14b So sánh mức độ phù hợp đường tần suất lý luận đường thực nghiệm dòng chảy tháng kiệt – lưu vực ND1 – kịch A1B 73 Hình 14c So sánh mức độ phù hợp đường tần suất lý luận đường thực nghiệm dòng chảy tháng kiệt – lưu vực ND1 – kịch A2 Bảng Tần suất dòng chảy tháng kiệt Subbasin P 5% 10% 25% 50% 75% 90% 95% Subbasin P 5% 10% 25% 50% 75% 90% 95% ND1 A1B 1.75 2.15 3.03 4.42 6.46 9.09 11.16 ND4 Baseline A1B 1.27 1.07 1.54 1.23 1.95 1.47 3.04 2.04 4.35 2.67 6.00 3.39 7.28 3.91 Baseline 1.61 1.95 2.69 3.85 5.50 7.60 9.21 A2 1.89 2.37 3.46 5.27 8.02 11.72 14.71 A2 1.01 1.26 1.82 2.74 4.13 5.98 7.47 ND2 A1B 0.65 0.76 0.99 1.34 1.79 2.34 2.75 ND5 Baseline A1B 0.91 0.65 1.06 0.75 1.29 0.93 1.86 1.20 2.50 1.53 3.27 1.92 3.83 2.19 Baseline 0.64 0.75 0.99 1.34 1.83 2.41 2.84 74 A2 0.68 0.82 1.10 1.53 2.12 2.86 3.41 A2 0.58 0.72 1.04 1.57 2.37 3.43 4.28 Baseline 0.8 0.99 1.42 2.11 3.15 4.52 5.61 ND3 A1B 0.70 0.83 1.11 1.54 2.13 2.85 3.40 A2 0.75 0.94 1.38 2.11 3.24 4.75 5.97 3.2.3.5 Đặc trưng thiếu hụt Cả đường cong thời khoảng dòng chảy đường tần suất dòng chảy kiệt không cung cấp thông tin độ dài giai đoạn liên tục giá trị dòng chảy quan tâm Phương pháp mô tả không đưa số lượng thiếu hụt dòng chảy, mà xây dựng giai đoạn dòng chảy kiệt liên tục Lý thuyết thời đoạn thiếu hụt sử dụng để đánh giá khả hồ chứa cần thiết đề cung cấp lượng xác định Hình 15a Dòng chảy kiệt đặc trưng khô hạn năm 1977 Một cách tiếp cận phổ biến sử dụng khái niệm ngưỡng (Yevjevich, 1967) [39] Một thời đoạn phạm vi dòng chảy kiệt số ngày dòng chảy trì dòng chảy ngưỡng xác định Giá trị ngưỡng cho dòng chảy thường niên dải lưu lượng với xác suất vượt từ 70 – 90% từ đường cong FDC, nghiên cứu lựa chọn giá trị Q90 dòng chảy tháng nhiều năm lưu vực (5.5 m3/s với lưu vực Nhuệ Đáy) ngưỡng lưu lượng thấp để trì môi trường sinh thái bình thường cho lưu vực Ba đặc trưng xem xét lý thuyết thời đoạn khoảng thời gian, tính nghiêm trọng (tổng thiếu hụt tích lũy) 75 cường độ Tổng lượng dòng chảy thiếu hụt giai đoạn Q def xác định tổng (Qthreshold-Qday) cho tất ngày Qday < Qthreshold từ chuỗi giá trị lưu lượng trung bình ngày Qday T3=149 ngày T1=85 ngày Hình 3.15b Dòng chảy kiệt đặc trưng khô hạn năm 2040 76 Dòng chảy năm khô 30 năm toàn lưu vực giảm mạnh 56.4%, 28.27% theo kịch A1B A2, dòng chảy năm ẩm tăng 5.57% toàn tác động kịch A1B, tăng 15.3% tác động kịch A2, có khác biệt tiểu lưu vực Dựa vào tính toán dòng chảy, năm 1977 xác định năm khô giai đoạn 1971 -1999, với dòng chảy năm 355 m3/s Trong giai đoạn 2010 – 2049 năm 2040 năm khô với dòng chảy năm đạt 87.79 m3/s Vì chuỗi dòng chảy ngày năm 1977 2040 phân tích chi tiết Hình 15a cho thấy lưu lượng ngưỡng Q90 năm 1977 xuất vài đợt diễn khoảng thời đoạn ngắn, đợt dài kéo dài tuần, thể tích thiếu hụt tích lũy nhỏ Năm 2040, tranh hoàn toàn khác, mức độ hạn tăng nghiêm trọng với nhiều đợt kéo dài với đợt kéo dài 85 ngày từ 21/5 đến 13/8/2040, đợt kéo dài từ 29/9 đến 12/11 (45 ngày), đợt kéo dài 149 ngày từ 3/12/2040 đến 30/4/2041 với tổng lượng nước thiếu hụt đợt tương ứng 185.3 m3, 141.7 m3 235.35 m3 tương ứng (hình 15b) 77 KẾT LUẬN 1/ Qua tổng quan tìm hiểu phương pháp đánh giá tài nguyên nước tác động biến đổi khí hậu phương pháp hiệu chỉnh số liệu khí tượng, từ lựa chọn cách tính toán thủy văn offline, nhận định phương pháp thực tiễn hiệu nghiên cứu biến đổi khí hậu đến tài nguyên nước nay, lựa chọn phương pháp hạ quy mô thống kê kết hợp với kết hạ quy mô động lực thông qua hệ số hiệu sai nhằm cải thiện kết đầu mô hình khí hậu Qua nghiên cứu đặc điểm địa lý tự nhiên lưu vực tình hình phát triển kinh tế xã hội sông Nhuệ Đáy thấy khu vực đóng vai trò vùng kinh tế trọng điểm , tập trung đông dân dẫn đến tăng nhu cầu nước, trình đô thị hóa diễn mạnh, Hà Nội có khả phải chịu tác động nghiêm trọng biến đổi khí hậu Một điểm cần ý tượng úng ngập Hà Nội thường xuyên xảy ra, nhiều trận mưa chí với cường độ 50 mm/giờ gây úng ngập nhiều khu vực, gây hậu nghiêm trọng nhiều ngành kinh tế quốc dân 2/ Đã tìm hiểu kịch biến đổi khí hậu mô hình toán mô khí hậu – dòng chảy, đặc biệt mô hình NAM, mô hình lựa chọn làm công cụ khảo sát kịch biến đổi khí hậu A1B A2 lưu vực sông Nhuệ Đáy Kết cho thấy mô hình có khả phân tích tác động biến đổi khí hậu cực trị đại lượng trung bình dòng chảy kịch sử dụng nhân tố khí hậu mô hình hóa 3/ Mô hình NAM khảo sát tác động biến đổi khí hậu thông qua kịch cho kết sau: Tác động biến đổi khí hậu có tính phân kỳ mạnh theo không gian Ở tiểu lưu vực thượng lưu, dạng biểu đồ thủy văn biến đổi đáng kể, dòng chảy lũ tập trung vào tháng từ tháng VII đến tháng X, tháng đỉnh lũ dịch chuyển 78 cuối năm, tháng X Dòng chảy lũ cường độ đỉnh lũ thể xu tăng trung bình 6% 16% so với Trong dòng chảy lũ yếu tố mang tính cực trị biến đổi lớn, thể qua cường độ thay đổi dòng chảy lũ tăng 6% dòng chảy tháng tăng 11% dòng chảy tháng lớn tăng 16% Phân tích tần suất dòng chảy lũ cho thấy kiện lũ biến đổi lớn, diễn thường xuyên với cường độ trận lũ với tần suất 1% 1200 m 3/s 2800 m3/s (so với 698 m3/s tại) tác động kịch A1B A2, cường độ lũ với thời gian lặp lại 100 năm có khả xuất lại với tần suất 20 năm tương lai Kết luận phù hợp với xu biến đổi theo nghiên cứu IPCC Trên tiểu lưu vực hạ lưu, đặc trưng lũ thể xu hướng biến đổi đối lập tác động kịch BDKH Không giống với kịch A2 gây tác động tương tự thượng lưu, kịch A1B có thiên hướng làm giảm rủi ro lũ so với tại, với dòng chảy lũ cường độ tháng đỉnh lũ giảm 5%, 21% tương ứng, cường độ lũ tần suất 1% giảm xuống 342 m3/s (hiện 528 m3/s), độ lặp lại trận lũ 10 năm tăng lên 40 năm Kết phân tích dòng chảy kiệt cho thấy tranh hoàn toàn tương tự dòng chảy lũ khác cường độ, điều cho thấy lưu vực Nhuệ Đáy thuộc hệ thống sông Hồng - Thái Bình lại biến đổi theo quy luật riêng tác động biến đổi khí hậu Dòng chảy kiệt thể tác động phân kỳ theo không gian Mặc dù có xen kẽ chu kỳ nhiều nước nước 10 năm dòng chảy thể xu biến đổi chính, tăng thượng lưu giảm hạ lưu Với trung binh 25% tăng, 8% giảm thượng hạ lưu tương ứng Cực trị kiệt có xu hướng dao động thiên giảm so với đại lượng trung bình kiệt, với khoảng biến đổi trung bình từ +23% đến -6% 79 Hàm phân bố log chuẩn đề nghị sử dụng phân tích tần suất cực trị dòng chảy lưu vực Nhuệ Đáy thể mức độ phù hợp tốt so với hàm phân bố xét đến (Weibull, Gumbel, Pearson III) Đặc trưng thiếu hụt năm khô giai đoạn tương lai xem xét cho thấy năm 2040 xuất năm kiệt khoảng từ 1970 đến 2050 Xảy đợt kiệt liên tục với đợt dài 149 ngày với tổng lượng nước thiết hụt khoảng 600 m3 so với giá trị ngưỡng Q90 lấy từ đường cong thời khoảng dòng chảy tháng giai đoạn Một cách tổng quát thấy rõ tác động biến đổi khí hậu trội hẳn với biến đổi theo không gian, biến đổi dòng chảy thể phân kỳ theo không gian biến đổi theo thời gian chúng diễn theo xu tương tự với biên dao động tăng dòng lũ kiệt khác cường độ Từ khía cạnh quản lý lũ, kết cho thấy biến động lớn, đại diện khoảng chênh lệch kết từ kịch Trong hướng thay đổi biết đến toàn lưu vực, cường độ thay đổi biến đổi đáng kể theo độ phân giải thời gian xem xét Các kết dường phân tán lưu vực sông Đáy sông Nhuệ, biểu qua khoảng biến động bé tính toán biến đổi khí hậu, với khả khoảng biến đổi nhỏ chế độ dòng chảy tương lai lưu vực Với số từ biến đổi khí hậu cho thấy xu hướng có khả tiếp tục biến đổi lớn tương lai, sách quy hoạch nên đặt quy trình để thực Cần thiết phải có kế hoạch quản lý lũ thời gian dài để kết hợp yếu tố biến đổi khí hậu chiến lược đánh giá rủi ro bao gồm thay đổi khí hậu khoảng biến động 4/ Một số hướng nghiên cứu Phải thừa nhận nghiên cứu tác động biến đổi khí hậu nhiều bất định không mô hình hóa thủy văn mả đánh giá tác động tương lai 80 Đối với mô hình hóa thủy văn, luận văn cố gắng loại bỏ bất định thông qua việc hiệu chỉnh, kiểm định mô hình thủy văn tính toán nhằm mục đích tối thiểu hóa sai số từ mô hình, biến đổi trình mưa – dòng chảy diễn điều kiện biến đổi tương lai, việc thay đổi thông số mô hình theo thời gian cần thiết Những thay đổi không đáng kể so với cường độ chế độ lũ hoàn toàn khác với dòng chảy kiệt việc tìm quy luật thay đổi góp phần cung cấp kết đánh giá xác, có độ tin cậy cao Đối với khía cạnh mô hình hóa khí hậu đánh giá tác động, độ phân giải thời gian không gian trình thủy văn khác so với đầu mô hình khí hậu, bước đầu khắc phục vấn đề thực luận văn nhiên hạn chế mặt thời gian, hệ thống máy móc, kết luận văn dừng lại mức độ định Do việc cải thiện kết đầu mô hình khí hậu cải thiện cầu nối thủy văn khí hậu vấn đề cần nghiên cứu nhiều, có hướng tiếp cận để khắc phục vấn đề này, từ khía cạnh thủy văn xây dựng mô đun thủy văn tích hợp vào mô hình khí hậu sử dụng mưa dạng lưới độ phân giải cao, hai từ khía cạnh khí tượng nghiên cứu phương pháp chuyển đổi liệu khí tượng sang liệu đáp ứng yêu cầu mô hình thủy văn với độ tin cậy cao 81 TÀI LIỆU THAM KHẢO Tiếng Việt Bộ Tài nguyên Môi trường, 2009 Kịch biến đổi khí hậu, nước biển dâng cho Việt Nam, Hà Nội Nguyễn Văn Cư, 2005 Báo cáo tổng kết đề tài cấp nhà nước: Xây dựng đề án tổng thể bảo vệ môi trường lưu vực sông Nhuệ sông Đáy – Hà Nội, Lưu trữ Viện Địa lý Hồ Thị Minh Hà, 2008 Nghiên cứu khả mô mùa yếu tố khí tượng lãnh thổ Việt Nam phương pháp thủy động thống kê, Luận án tiến sỹ khí tượng học, Trường Đại học KHTN Hà Nội Văn Thị Hằng 2010, Đánh giá tác động biến đổi khí hậu đến biến động tài nguyên nước lưu vực sông Nhuệ - Đáy thuộc thành phố Hà Nội, Luận văn thạc sỹ khoa học, Trường Đại học KNTN Hà Nội Vũ Văn Minh, Nguyễn Hoàng Minh, Trần Hồng Thái, 2011 Đánh giá tác động Biến đổi khí hậu đến dòng chảy lũ lưu vực sông Hồng - Thái Bình Tuyển tập Báo cáo Khoa học lần thứ XIII, tr 72 – 78 Nguyễn Ý Như, 2009 Ứng dụng mô hình SWAT nghiên cứu ảnh hưởng biến đổi khí hậu sử dụng đất đến dòng chảy sông Bến Hải Khóa luận tốt nghiệp, Trường Đại học Khoa học Tự nhiên Nguyễn Ý Như, Lê Văn Linh, Nguyễn Thanh Sơn, Trần Ngọc Anh, 2011 Nghiên cứu tác động biến đổi khí hậu đến dòng chảy lưu vực sông Đáy địa bàn thành phố Hà Nội Tạp chí Khoa học ĐHQGHN, Khoa học Tự nhiên Công nghệ Tập 27, số 1S, tr 192 Nguyễn Ý Như, Nguyễn Thanh Sơn, 2011 Biến đổi dòng chảy kiệt bối cảnh Biến đổi khí hậu lưu vực sông Nhuệ Đáy Tạp chí Khoa học ĐHQGHN, Khoa học Tự nhiên Công nghệ Tập 27, số 3S, tr 39 – 43 Nguyễn Ý Như, Nguyễn Thanh Sơn, Trần Ngọc Anh, Nguyễn Quang Trung, 2011 The potential impacts of climate change on flood flow in Nhue – Day 82 river basin The second International MAHASRI/HyARC Workshop, August 2224, 2011, Nha Trang, Vietnam 10 11 12 13 14 15 16 17 Lê Vũ Việt Phong, 2006 Nghiên cứu áp dụng mô hình toán MIKE 11 tính toán chất lượng nước sông Nhuệ sông Đáy Khóa luận tốt nghiệp trường Đại học Thủy lợi Nguyễn Thanh Sơn, 2008 Nghiên cứu mô trình mưa – dòng chảy phục vụ sử dụng hợp lý tài nguyên nước đất số lưu vực sông thượng nguồn miền Trung Luận án tiến sỹ Địa lý, Trường Đại học Khoa học Tự nhiên Hà Nội Nguyễn Thanh Sơn, Ngô Chí Tuấn, Văn Thị Hằng, Nguyễn Ý Như, 2011 Ảnh hưởng biến đổi khí hậu đến biến đổi tài nguyên nước lưu vực sông Nhuệ Đáy Tạp chí Khoa học ĐHQGHN, Khoa học Tự nhiên Công nghệ Tập 27, số 1S, tr 218 Nguyễn Thanh Sơn, Nguyễn Ý Như, Trần Ngọc Anh, Lê Thị Hường, 2011 Khảo sát trạng tài nguyên nước lưu vực sông Nhuệ - Đáy, Tạp chí Khoa học ĐHQGHN, Khoa học Tự nhiên Công nghệ Tập 27, số1S, tr 227 Nguyen Thanh Son, Nguyen Y Như, 2009, Applying SWAT model to simulate streamflow in BenHai River Basin in response to climate change scenarios Journal of Science, Earth Sciences, VNU, V25, No3, tr 161 Trung tâm Tư vấn Khí tượng Thủy văn Môi trường, Viện KH KTTV & MT, 2010 Đánh giá tác động biến đổi khí hậu lên tài nguyên nước biện pháp thích ứng – Lưu vực sông Hồng – Thái Bình, Lưu trữ Viện KH KTTV & MT Viện Khoa học Khí tượng Thủy văn & Môi trường, 2010 Tác động biến đổi khí hậu lên tài nguyên nước biện pháp thích ứng, Lưu trữ Viện KH KTTV & MT Trần Thanh Xuân, Hoàng Minh Tuyển, Lê Tuấn Nghĩa, Lương Hữu Dũng, 2011 Tác động Biến đổi khí hậu đến dòng chảy sông Tuyển tập Báo cáo Khoa học lần thứ XIII, tr 146 – 153 Tiếng Anh 83 18 19 20 21 22 23 24 Andersen H.E., Kronvang B., Larsen S.E., Hoffmann C.C., Jensen T.S, Rasmussen E.K., 2006 Climate change impacts on hydrology and nutrients in a Danish lowland river basin Science of the Total Environment 365, pp 223 – 237 Bonacci, O., Tadic, Z., Trninic, D., 1992 Effects of dams and reservoirs on the hydrological characteristics of the lower Drava River Regul Rivers Res Manag (4), pp 349–357 Cooley, K.R 1990: Effects of CO2 – induced climate changes on snowpack and streamflow Hydrological Sciences Journal 35, pp 511 – 22 D Labat, Y Godderis, J L Probst, 2004 Evidence for global runoff increase related to climate warming Advances in Water Resources, 27, pp 631 – 642 DHI, 2007 Mike 11 – User’s Manual Feddes, R.A., M Menenti, and P Kabat, 1989 Modeling the soil water and surface energy balance in relation to climate models, European coordination group on land – surface processes, hydrology, Dessertification in Europe, Barcelona, pp 21 Hayhoe K.A., 2010 A standardized framework for evaluating the skill of regional climate downscaling techniques Doctor of Philosophy in Atmospheric, Sciences in the Graduate College of the University of Illinois at Urbana – Champaign 25 IPCC (2007), Climate Change 2007: Synthesis Report 26 IPCC (2001), Special Report on Emissions Scenarios 27 28 29 Gleick, P.H., 1987a Regional hydrologic consequences of increases in atmospheric CO2 and other trace gases, Climate Change 10, tr 137 – 61 Kalvova, J., Nemesova, I., 1997 Projections of climate change for the Czech Republic Climate Change 36, tr 41 – 64 Kim U., Kaluarachchi J J., Smakhtin V U., 2008 Climate Change Impacts on Hydrology and Water Resources of the Upper Blue Nile River Basin, Ethiopia, Colombo, Sri Lanka: International Water Management Institute pp 27 (IWMI Research Report 126) 84 30 31 32 33 34 35 36 37 38 39 Mark Maslin, 2004 GLOBAL WARMING: A Very Short Introduction, OXFORD University Press Michal Jenicek, 2007 Rainfall-runoff modelling in small and middle-large catchments – an overview Nemec, J and Schaake, J., 1982 Sensitivity of water resources systems to climate variation Hydrological Sciences Journal 2, pp 327 – 43 S.L.Neitsch, J.G Arnold, J.R.Kiniry, J.R.Williams, 2001 Soil and water assessment tool user’s manual, USDA_ARS Publications Nash, J E and J V Sutcliffe, 1970 River flow forecasting through conceptual models part I — A discussion of principles, Journal of Hydrology, 10 (3), 282– 290 Sapkota M., Hamagochi T., Kojiri T., 2010 Geostatistical bias correction of super high resolution GCM outputs under climate change and its application to runoff simulations in Red river basin Proceedings of the fifth conference of APHW Conference, Labor and Social Publisher US Army Corps of Engineers (2001), Hydrology Model System HEC-HMS Users’ Manual Wigley, T.M.L, Jones, P.D., 1987 England and Wales precipitations: a discussion of recent changes in variability and an update to 1985 J Climatol 7, 231 – 246 Z.X Xu, F.F Zhao, J.Y Li, 2008 Response of streamflow to climate change in the headwater catchment of the Yellow River basin Quaternary International xxx, pp – 14 Yevjevich, V., 1972 Probability and Statistics in Hydrology Water Resources Publications, Fort Collins, CO, USA, pp 302 85 [...]... Lưu vực sông Đáy Lưu vực sông Đáy trên địa bàn thành phố Hà Nội với chiều dài 114 km Các chi lưu của sông Đáy: sông Tích, sông Bùi, sông Thanh Hà Nói chung 85% lượng dòng chảy trên lưu vực sông Đáy trên địa bàn thành phố Hà Nội có nguồn gốc từ sông Hồng chuyển sang, chỉ 15% còn lại bắt nguồn từ lưu vực Trong mùa mưa, mực nước và lưu lượng các sông suối lớn thay đổi nhanh, tốc độ dòng chảy đạt từ 2-... đổi khí hậu đến cực trị dòng chảy trên lưu vực sông Nhuệ Đáy thuộc địa bàn Thành phố Hà Nội đã được hình thành, có ý nghĩa khoa học và thực tiễn sâu sắc Tiếp theo, luận văn đi vào tổng quan các đặc điểm địa lý tự nhiên và tài nguyên nước của lưu vực sông Nhuệ - Đáy trong phần tiếp theo 1.4 ĐIỀU KIỆN ĐỊA LÝ TỰ NHIÊN – KINH TẾ XÃ HỘI LƯU VỰC SÔNG NHUỆ ĐÁY THUỘC ĐỊA PHẬN THÀNH PHỐ HÀ NỘI 1.4.1 Điều kiện... xã Hà Đông và chảy vào sông Đáy tại thành phố Phủ Lý Nước sông Tô Lịch thường xuyên xả vào sông Nhuệ với lưu lượng trung bình từ 11-17m3/s, lưu lượng cực đại đạt 30m3/s Các sông chính trong lưu vực: sông Nhuệ, sông Thanh Hà, sông Tích, sông Bùi và Sông Tô Lịch là nhánh chính của sông Nhuệ, nhận nước từ sông Lừ, Kim Ngưu, Sét Ngoài 2 sông chính ra có rất nhiều các chi lưu như: sông Tích, sông Bùi, sông. .. dòng chảy trung bình 4/ Các nghiên cứu trong nước hiện chỉ tập trung vào những lưu vực có quy mô lớn hệ thống sông Hồng – Thái Bình hay sông Mê Công mà chưa có những nghiên cứu cho các lưu vực nhỏ, như lưu vực sông Nhuệ - Đáy, một lưu vực đóng vai trò đặc biệt với Thủ đô Hà Nội, nơi tập trung đông dân cư và quá trình đô thị hóa diễn ra mạnh Đề tài luận văn Nghiên cứu tác động của Biến đổi khí hậu đến. .. chảy thẳng vào sông Đáy Diện tích lưu vực là 271 km2, sông dài 40 km, chiều rộng trung bình lưu vực 9 km Chế độ thủy văn lưu vực sông Đáy không những chịu ảnh hưởng của các yếu tố mặt đệm trên bề mặt lưu vực, các yếu tố khí hậu mà còn phụ thuộc vào chế độ dòng chảy của nước sông Hồng và các sông khác Vì thế mà chế độ thủy văn ở đây rất phức tạp và có sự khác nhau nhất định giữa các đoạn sông Dòng chảy. .. hướng giảm dòng chảy trung bình năm và tăng lượng thiếu hụt tài nguyên nước trên lưu vực nghiên cứu, tuy nhiên biến động thủy văn tương ứng với mỗi số liệu biến đổi GCMs tương đối lớn Kim U & nnk (2008) [29] đã đánh giá tác động của biến đổi khí hậu đối với cả chế độ thủy văn và tài nguyên nước trên lưu vực thượng lưu sông Blue Nile ở Ethiopia, khu vực có dữ liệu quan trắc hạn chế Nghiên cứu này cũng... Mực nước và lưu lượng lớn nhất năm có có khả năng xuất hiện trong tháng VII, VIII, hoặc IX, nhưng phổ biến vào tháng VIII Lưu vực sông Nhuệ Sông Nhuệ là một con sông nhỏ dài 74 km, nhận nước từ sông Hồng và lượng mưa trên toàn lưu vực khống chế của sông Nhuệ cuối cùng đổ vào sông Đáy qua cống Lương Cổ ở khu vực thành phố Phủ Lý Về mùa kiệt cống Liên Mạc luôn mở để lấy nước sông Hồng vào sông Nhuệ, còn... cống La Khê, Ngoại Độ…Các sông này thường phải đóng lại khi có phân lũ trong thời gian dài, ngắn tùy thuộc vào thời gian lũ Sông Đáy có vị trí rất quan trọng, nó vừa là đường thoát nước chính của sông Hồng, vừa là đường tiêu lũ của bản thân lưu vực sông Đáy 1.4.2 Đặc điểm kinh tế - xã hội Phân bố dân cư Dân số trên lưu vực sông Nhuệ - sông Đáy trên lưu vực nghiên cứu ước tính đến năm 2009 là 10,77 triệu... ảnh hưởng đến rất nhiều lĩnh vực trong đó có tài nguyên nước Trong khoảng 10 – 15 năm qua đã có nhiều nhà thủy văn trên thế 11 giới nghiên cứu tác động của biến đổi khí hậu đối với tài nguyên nước Trong những nghiên cứu này vận dụng nhiều cách tiếp cận các mô hình khác nhau Dù là theo cách tiếp cận nào thì mục tiêu chính của các hoạt động nghiên cứu tài nguyên nước liên quan đến biến đổi khí hậu là nhằm... đập Đáy nước sông Hồng không thường xuyên vào sông Đáy qua cửa đập Đáy trừ những năm phân lũ, vì vậy phần đầu nguồn sông (từ km 0 đến Ba Thá dài 71km) sông Đáy coi như đoạn sông chết Lượng nước để nuôi sông Đáy chủ yếu là do các sông nhánh, quan trọng nhất là sông Tích Sông Nhuệ lấy nước từ sông Hồng qua cống Liên Mạc để tưới cho hệ thống thủy nông Đan Hoài Sông Nhuệ còn tiêu nước cho thành phố Hà Nội,

Ngày đăng: 18/06/2016, 15:34

Nguồn tham khảo

Tài liệu tham khảo Loại Chi tiết
1. Bộ Tài nguyên và Môi trường, 2009. Kịch bản biến đổi khí hậu, nước biển dâng cho Việt Nam, Hà Nội Sách, tạp chí
Tiêu đề: Kịch bản biến đổi khí hậu, nước biển dângcho Việt Nam
2. Nguyễn Văn Cư, 2005. Báo cáo tổng kết đề tài cấp nhà nước: Xây dựng đề án tổng thể bảo vệ môi trường lưu vực sông Nhuệ và sông Đáy – Hà Nội, Lưu trữ Viện Địa lý Sách, tạp chí
Tiêu đề: Báo cáo tổng kết đề tài cấp nhà nước: Xây dựng đề ántổng thể bảo vệ môi trường lưu vực sông Nhuệ và sông Đáy – Hà Nội
3. Hồ Thị Minh Hà, 2008. Nghiên cứu khả năng mô phỏng mùa các yếu tố khí tượng trên lãnh thổ Việt Nam bằng phương pháp thủy động và thống kê, Luận án tiến sỹ khí tượng học, Trường Đại học KHTN Hà Nội Sách, tạp chí
Tiêu đề: Nghiên cứu khả năng mô phỏng mùa các yếu tố khítượng trên lãnh thổ Việt Nam bằng phương pháp thủy động và thống kê
4. Văn Thị Hằng 2010, Đánh giá tác động của biến đổi khí hậu đến biến động tài nguyên nước lưu vực sông Nhuệ - Đáy thuộc thành phố Hà Nội, Luận văn thạc sỹ khoa học, Trường Đại học KNTN Hà Nội Sách, tạp chí
Tiêu đề: Đánh giá tác động của biến đổi khí hậu đến biến động tàinguyên nước lưu vực sông Nhuệ - Đáy thuộc thành phố Hà Nội
5. Vũ Văn Minh, Nguyễn Hoàng Minh, Trần Hồng Thái, 2011. Đánh giá tác động của Biến đổi khí hậu đến dòng chảy lũ lưu vực sông Hồng - Thái Bình. Tuyển tập Báo cáo Khoa học lần thứ XIII, tr. 72 – 78 Sách, tạp chí
Tiêu đề: Tuyểntập Báo cáo Khoa học lần thứ XIII
6. Nguyễn Ý Như, 2009. Ứng dụng mô hình SWAT nghiên cứu ảnh hưởng của biến đổi khí hậu và sử dụng đất đến dòng chảy sông Bến Hải . Khóa luận tốt nghiệp, Trường Đại học Khoa học Tự nhiên Sách, tạp chí
Tiêu đề: Ứng dụng mô hình SWAT nghiên cứu ảnh hưởng củabiến đổi khí hậu và sử dụng đất đến dòng chảy sông Bến Hải
7. Nguyễn Ý Như, Lê Văn Linh, Nguyễn Thanh Sơn, Trần Ngọc Anh, 2011.Nghiên cứu tác động của biến đổi khí hậu đến dòng chảy lưu vực sông Đáy trên địa bàn thành phố Hà Nội. Tạp chí Khoa học ĐHQGHN, Khoa học Tự nhiên và Công nghệ Tập 27, số 1S, tr. 192 Sách, tạp chí
Tiêu đề: Tạp chí Khoa học ĐHQGHN, Khoa học Tự nhiên vàCông nghệ
8. Nguyễn Ý Như, Nguyễn Thanh Sơn, 2011. Biến đổi dòng chảy kiệt trong bối cảnh Biến đổi khí hậu trên lưu vực sông Nhuệ Đáy. Tạp chí Khoa học ĐHQGHN, Khoa học Tự nhiên và Công nghệ Tập 27, số 3S, tr. 39 – 43 Sách, tạp chí
Tiêu đề: Tạp chí Khoa họcĐHQGHN, Khoa học Tự nhiên và Công nghệ
10. Lê Vũ Việt Phong, 2006. Nghiên cứu áp dụng mô hình toán MIKE 11 tính toán chất lượng nước sông Nhuệ và sông Đáy. Khóa luận tốt nghiệp trường Đại học Thủy lợi Sách, tạp chí
Tiêu đề: Nghiên cứu áp dụng mô hình toán MIKE 11 tính toánchất lượng nước sông Nhuệ và sông Đáy
11. Nguyễn Thanh Sơn, 2008. Nghiên cứu mô phỏng quá trình mưa – dòng chảy phục vụ sử dụng hợp lý tài nguyên nước và đất một số lưu vực sông thượng nguồn miền Trung. Luận án tiến sỹ Địa lý, Trường Đại học Khoa học Tự nhiên Hà Nội Sách, tạp chí
Tiêu đề: Nghiên cứu mô phỏng quá trình mưa – dòng chảyphục vụ sử dụng hợp lý tài nguyên nước và đất một số lưu vực sông thượngnguồn miền Trung
12. Nguyễn Thanh Sơn, Ngô Chí Tuấn, Văn Thị Hằng, Nguyễn Ý Như, 2011. Ảnh hưởng của biến đổi khí hậu đến biến đổi tài nguyên nước lưu vực sông Nhuệ - Đáy. Tạp chí Khoa học ĐHQGHN, Khoa học Tự nhiên và Công nghệ Tập 27, số 1S, tr. 218 Sách, tạp chí
Tiêu đề: Tạp chí Khoa học ĐHQGHN, Khoa học Tự nhiên và Công nghệ
13. Nguyễn Thanh Sơn, Nguyễn Ý Như, Trần Ngọc Anh, Lê Thị Hường, 2011.Khảo sát hiện trạng tài nguyên nước lưu vực sông Nhuệ - Đáy, Tạp chí Khoa học ĐHQGHN, Khoa học Tự nhiên và Công nghệ Tập 27, số1S, tr. 227 Sách, tạp chí
Tiêu đề: Tạp chí Khoahọc ĐHQGHN, Khoa học Tự nhiên và Công nghệ
14. Nguyen Thanh Son, Nguyen Y Như, 2009, Applying SWAT model to simulate streamflow in BenHai River Basin in response to climate change scenarios.Journal of Science, Earth Sciences, VNU, V25, No3, tr. 161 Sách, tạp chí
Tiêu đề: Journal of Science, Earth Sciences
15. Trung tâm Tư vấn Khí tượng Thủy văn và Môi trường, Viện KH KTTV &amp; MT, 2010. Đánh giá tác động của biến đổi khí hậu lên tài nguyên nước và các biện pháp thích ứng – Lưu vực sông Hồng – Thái Bình, Lưu trữ Viện KH KTTV &amp;MT Sách, tạp chí
Tiêu đề: Đánh giá tác động của biến đổi khí hậu lên tài nguyên nước và các biệnpháp thích ứng – Lưu vực sông Hồng – Thái Bình
16. Viện Khoa học Khí tượng Thủy văn &amp; Môi trường, 2010. Tác động của biến đổi khí hậu lên tài nguyên nước và các biện pháp thích ứng, Lưu trữ Viện KH KTTV &amp; MT Sách, tạp chí
Tiêu đề: Tác động của biến đổikhí hậu lên tài nguyên nước và các biện pháp thích ứng
17. Trần Thanh Xuân, Hoàng Minh Tuyển, Lê Tuấn Nghĩa, Lương Hữu Dũng, 2011. Tác động của Biến đổi khí hậu đến dòng chảy trong sông Tuyển tập Báo cáo Khoa học lần thứ XIII, tr. 146 – 153.Tiếng Anh Sách, tạp chí
Tiêu đề: Tác động của Biến đổi khí hậu đến dòng chảy trongsông Tuyển tập Báo cáo Khoa học lần thứ XIII
9. Nguyễn Ý Như, Nguyễn Thanh Sơn, Trần Ngọc Anh, Nguyễn Quang Trung, 2011. The potential impacts of climate change on flood flow in Nhue – Day Khác

TÀI LIỆU CÙNG NGƯỜI DÙNG

TÀI LIỆU LIÊN QUAN

w