Tổng quan nghiên cứu

Lưu vực sông Nhuệ - sông Đáy giữ vị trí chiến lược đặc biệt đối với sự phát triển kinh tế - xã hội và an ninh nguồn nước của Thủ đô Hà Nội cũng như toàn bộ khu vực đồng bằng sông Hồng. Với tổng diện tích lưu vực sông Đáy trên địa bàn Hà Nội đạt 1900 km² và lưu vực sông Nhuệ khoảng 603 km², khu vực này hiện đang phục vụ cho hơn 10,77 triệu người với mật độ dân số trung bình lên đến 1405 người/km², cao gấp 5,5 lần bình quân cả nước. Tuy nhiên, sự gia tăng dân số đô thị cùng quá trình công nghiệp hóa mạnh mẽ đang đặt hệ thống tài nguyên nước trước những áp lực nặng nề.

Vấn đề nghiên cứu trọng tâm là sự tác động đa chiều của biến đổi khí hậu đến các hiện tượng thủy văn cực đoan, bao gồm lũ lụt khốc liệt trong mùa mưa và hạn hán, suy kiệt dòng chảy trong mùa khô. Mục tiêu cụ thể của luận văn là lượng hóa mức độ biến động của các đặc trưng cực trị dòng chảy (lưu lượng đỉnh lũ lớn nhất và dòng chảy kiệt nhất) trên lưu vực sông Nhuệ - Đáy thuộc Hà Nội trong giai đoạn tương lai 2010–2049 theo hai kịch bản phát thải A1B và A2 của Ban Liên chính phủ về Biến đổi khí hậu, so sánh với giai đoạn nền 1970–1999.

Nghiên cứu mang ý nghĩa khoa học và thực tiễn sâu sắc khi cung cấp hệ thống cơ sở dữ liệu định lượng chính xác cho các nhà quản lý. Các chỉ số dự báo về biên độ lũ, dòng chảy kiệt và tần suất cực trị đóng vai trò là căn cứ kỹ thuật then chốt để quy hoạch hệ thống thoát lũ sông Đáy, vận hành cống Liên Mạc và đập Đáy, đồng thời giảm thiểu thiệt hại kinh tế ước tính hàng trăm tỷ đồng do thiên tai ngập úng và ô nhiễm nước mặt gây ra cho thành phố Hà Nội.

Cơ sở lý thuyết và phương pháp nghiên cứu

Khung lý thuyết áp dụng

Nghiên cứu được xây dựng trên nền tảng lý thuyết biến đổi khí hậu toàn cầu của Ban Liên chính phủ về Biến đổi khí hậu (IPCC), kết hợp lý thuyết mô hình hóa thủy văn mưa - dòng chảy trong điều kiện biến đổi mặt đệm và khí quyển. Khung phân tích sử dụng hai kịch bản phát thải khí nhà kính chính: kịch bản A1B (phát thải trung bình, thể hiện sự phát triển công nghệ cân bằng và có xác suất xảy ra cao nhất) và kịch bản A2 (phát thải cao, mô tả kịch bản tiêu cực với tốc độ tăng trưởng phát thải liên tục).

Các khái niệm chuyên ngành then chốt bao gồm: cực trị dòng chảy (dòng chảy đỉnh lũ lớn nhất năm và lưu lượng dòng chảy tháng kiệt nhất), hạ quy mô khí hậu (climate downscaling), đường cong thời khoảng dòng chảy (Flow Duration Curve - FDC) và hệ số hiệu quả mô hình Nash-Sutcliffe (NSE). Về mặt thủy văn, mô hình NAM thuộc bộ công cụ MIKE 11 của Viện Thủy lực Đan Mạch (DHI) được lựa chọn. Đây là mô hình tất định, thông số tập trung hoạt động dựa trên cơ chế cân bằng ẩm liên tục giữa bốn tầng trữ nước tương tác: bể chứa tuyết tan, bể chứa bề mặt (với thông số dung tích trữ tối đa Umax), bể chứa tầng rễ cây (dung tích trữ tối đa Lmax) và hệ thống bể chứa nước ngầm hai tầng.

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

Nguồn dữ liệu nghiên cứu bao gồm chuỗi số liệu quan trắc khí tượng thủy văn theo bước thời gian ngày trong 30 năm giai đoạn nền (1970–1999) và chuỗi quan trắc liên tục 1971–1978 tại 7 trạm đo đại diện: Ba Thá, Lâm Sơn, Sơn Tây, Hà Đông, Hưng Thi, Phủ Lý và Láng. Phương pháp chọn mẫu số liệu là chọn chuỗi liên tục, đồng bộ và có tính đại diện cao nhất về không gian mặt đệm. Dữ liệu địa hình được trích xuất từ mô hình số độ cao DEM kích thước ô lưới 30x30m kết hợp công cụ AVSWAT để phân chia lưu vực nghiên cứu thành 4 tiểu lưu vực cơ sở đồng nhất về đặc tính dòng chảy.

Phương pháp phân tích tích hợp mô hình khí hậu khu vực RegCM3 độ phân giải ngang 36km với phương pháp hạ quy mô thống kê thông qua hệ số hiệu sai tỷ lệ theo từng tháng và thuật toán nội suy song tuyến tính để chuyển đổi dữ liệu từ dạng ô lưới về các điểm trạm đo thực tế. Lý do lựa chọn tổ hợp phương pháp này là vì mô hình RegCM3 nắm bắt tốt các biến đổi vật lý quy mô khu vực, trong khi phương pháp hạ quy mô thống kê giúp loại bỏ sai số hệ thống giữa mô hình khí hậu và số liệu quan trắc thực địa. Mô hình NAM được lựa chọn nhờ cấu trúc gọn gàng, khả năng mô phỏng liên tục dòng chảy ngày vượt trội và phù hợp tối ưu với điều kiện dữ liệu khí tượng thủy văn tại lưu vực sông Nhuệ - Đáy so với các mô hình phân bố phức tạp đòi hỏi quá nhiều thông số địa chất thổ nhưỡng.

Kết quả nghiên cứu và thảo luận

Những phát hiện chính

Quá trình hiệu chỉnh và kiểm định mô hình mưa - dòng chảy NAM tại trạm thủy văn Ba Thá (tiểu lưu vực ND1) cho thấy mức độ chính xác rất cao. Trong giai đoạn hiệu chỉnh (1972–1974), hệ số hiệu quả Nash-Sutcliffe đạt 0,85 với sai số tổng lượng nước tích lũy chỉ 0,3% và hệ số tương quan đạt 0,811. Trong giai đoạn kiểm định độc lập (1976–1978), chỉ số Nash-Sutcliffe duy trì ở mức 0,86, sai số tổng lượng nước tích lũy là 0,4% và sai số lưu lượng đỉnh lũ dưới mức ước tính 8,23%.

Dòng chảy mùa lũ có xu hướng gia tăng rõ rệt về cả cường độ lẫn tính ác liệt trong giai đoạn 2010–2049. Dưới tác động của kịch bản phát thải A1B và A2, lưu lượng dòng chảy 4 tháng mùa lũ chính (chiếm 70–80% tổng lượng dòng chảy năm) và dòng chảy ngày cực đại đều tăng mạnh. Cụ thể, lưu lượng đỉnh lũ cực đại 3 ngày lớn nhất tại tiểu lưu vực Ba Thá tiếp nhận nước từ sông Tích và sườn đông núi Ba Vì (nơi có lượng mưa năm cao nhất đạt từ 2200 đến 4000 mm) có xu hướng gia tăng từ 10% đến 25% so với thời kỳ nền.

Ngược lại, dòng chảy mùa kiệt thể hiện sự suy giảm nghiêm trọng. Trong 7 tháng mùa cạn (vốn chỉ đóng góp 20–25% tổng lượng nước cả năm), lưu lượng dòng chảy tháng kiệt nhất theo kịch bản A1B và A2 giảm từ 8% đến 18% qua các thập niên dự báo. Mô phỏng năm điển hình khô hạn 2040 cho thấy mức độ cạn kiệt sâu hơn 20% so với năm hạn lịch sử 1977. Không gian biến động cho thấy vùng hữu ngạn sông Đáy chịu sự gia tăng lũ đột biến, trong khi dòng chính sông Nhuệ dài 74 km đối mặt với nguy cơ cạn dòng cơ sở, làm mất khả năng tự làm sạch khi tiếp nhận nước thải từ sông Tô Lịch với lưu lượng xả từ 11 đến 17 m³/s.

Thảo luận kết quả

Nguyên nhân chính dẫn đến sự phân hóa cực đoan của dòng chảy bắt nguồn từ sự mất cân bằng trong cán cân nhiệt - ẩm. Nhiệt độ không khí trung bình tăng đẩy lượng bốc hơi tiềm năng vượt mức 900–1000 mm/năm, làm cạn kiệt dung tích ẩm tầng rễ cây (Lmax) vào mùa khô. Đồng thời, sự gia tăng cường độ các trận mưa rào nhiệt đới tập trung trong thời gian ngắn khiến lượng nước mưa vượt ngưỡng thấm, chuyển hóa trực tiếp thành dòng chảy tràn trên bề mặt.

So sánh với các nghiên cứu quy mô lớn trên hệ thống sông Hồng - Thái Bình của các tác giả trong nước, kết quả của luận văn hoàn toàn nhất quán về xu thế lũ tăng từ 5% đến 22% và kiệt giảm sâu, nhưng nghiên cứu này đã đi sâu chi tiết hóa đến từng tiểu lưu vực và từng dạng cực trị dòng chảy ngày. Trong báo cáo kết quả, dữ liệu biến động được biểu diễn trực quan và sinh động thông qua hệ thống bảng tần suất Pearson loại III, biểu đồ đường quá trình dòng chảy mô phỏng (hydrograph) so sánh chuỗi thực đo với tính toán, cùng đồ thị đường cong thời khoảng dòng chảy (FDC) phản ánh rõ nét sự dịch chuyển phân bố lưu lượng kiệt giữa giai đoạn nền và hai kịch bản tương lai.

Đề xuất và khuyến nghị

Nhằm chủ động thích ứng với các tác động tiêu cực của biến đổi khí hậu đến chế độ thủy văn cực trị, luận văn đề xuất 4 nhóm giải pháp chiến lược:

  1. Nâng cấp và số hóa mạng lưới quan trắc cảnh báo sớm thiên tai: Triển khai lắp đặt bổ sung 5 đến 10 trạm đo mưa tự động và trạm quan trắc mực nước số hóa trên các nhánh sông Tích, sông Bùi và hạ lưu sông Đáy. Mục tiêu là nâng cao độ chính xác dự báo lũ trước 12–24 giờ, hoàn thành trước năm 2030 dưới sự chủ trì của Sở Nông nghiệp và Phát triển nông thôn Hà Nội phối hợp cùng Trung tâm Dự báo Khí tượng Thủy văn Quốc gia.

  2. Hoàn thiện quy trình vận hành liên hồ chứa và công trình phân lũ: Ban hành quy chế điều tiết mới cho đập Đáy, cống Liên Mạc và cống Lương Cổ nhằm chủ động cắt giảm từ 15% đến 20% lưu lượng đỉnh lũ sông Hồng chuyển sang khi có báo động khẩn cấp, đồng thời duy trì lưu lượng xả tối thiểu mùa kiệt trên 15 m³/s để duy trì dòng chảy môi trường trong giai đoạn 2025–2030 do Ủy ban nhân dân thành phố Hà Nội và Bộ Nông nghiệp và Phát triển nông thôn thực hiện.

  3. Bảo vệ và tái tạo thảm phủ thực vật đầu nguồn: Tăng cường khoanh nuôi, trồng mới và phục hồi tối thiểu 160,2 km² diện tích rừng phòng hộ tự nhiên và cải tạo 105,64 km² rừng nghèo tại vùng đồi núi Ba Vì và Kim Bôi trước năm 2035. Mục tiêu là nâng cao khả năng giữ nước mặt đệm, giảm hệ số dòng chảy tràn bề mặt (CQOF) xuống dưới mức 0,35 do Chi cục Kiểm lâm Hà Nội và Hòa Bình phụ trách.

  4. Xây dựng công trình kiểm soát ô nhiễm và trữ nước mùa kiệt: Đẩy nhanh tiến độ thi công hệ thống cống gom và các nhà máy xử lý nước thải tập trung dọc lưu vực sông Tô Lịch và sông Nhuệ, giảm 80% nồng độ các chất ô nhiễm hữu cơ trước năm 2030, do Sở Tài nguyên và Môi trường thành phố Hà Nội làm đầu mối triển khai.

Đối tượng nên tham khảo luận văn

Nội dung và kết quả nghiên cứu của luận văn mang lại giá trị ứng dụng thiết thực cho 4 nhóm đối tượng chính:

  1. Cơ quan quản lý nhà nước và hoạch định chính sách: Điển hình là Sở Tài nguyên và Môi trường, Sở Nông nghiệp và Phát triển nông thôn Hà Nội cùng các cơ quan ban ngành liên quan. Tài liệu cung cấp căn cứ khoa học trực tiếp để xây dựng Kế hoạch hành động ứng phó với biến đổi khí hậu và quy hoạch tích hợp tài nguyên nước toàn lưu vực.

  2. Đơn vị tư vấn, thiết kế công trình thủy lợi và hạ tầng đô thị: Các kỹ sư xây dựng thủy lợi, giao thông, cấp thoát nước có thể khai thác các thông số tính toán thủy văn cực trị, lưu lượng lũ thiết kế tương ứng với các chu kỳ lặp để xác định khẩu độ cầu cống, cao trình đê kè và công suất trạm bơm tiêu úng dọc sông Nhuệ - Đáy.

  3. Giảng viên, nghiên cứu sinh và học viên chuyên ngành: Các nhà khoa học trong lĩnh vực Thủy văn học, Biến đổi khí hậu, Quản lý tài nguyên nước và Môi trường có thể sử dụng luận văn như tài liệu tham khảo mẫu mực về quy trình liên kết mô hình khí hậu RegCM3, phương pháp hạ quy mô thống kê và mô hình thủy văn NAM.

  4. Các tổ chức phi chính phủ và cộng đồng bảo vệ môi trường: Cung cấp bức tranh toàn cảnh về nguy cơ suy thoái chất lượng nước mặt trên 74 km dòng chính sông Nhuệ và 114 km sông Đáy, hỗ trợ xây dựng các đề án phục hồi hệ sinh thái sông ngòi đô thị.

Câu hỏi thường gặp

Tại sao nghiên cứu lại lựa chọn hai kịch bản phát thải A1B và A2 thay vì các kịch bản khác? Theo khuyến nghị của Bộ Tài nguyên và Môi trường và nhóm nghiên cứu khí hậu, kịch bản A1B mô tả sự phát triển cân bằng công nghệ có tính khả thi cao nhất trong thế kỷ 21, còn A2 phản ánh kịch bản phát thải cao và tiêu cực nhất. Việc kết hợp hai kịch bản này giúp định hình đầy đủ dải biến động từ mức trung bình đến khắc nghiệt nhất của cực trị dòng chảy trong 40 năm tới.

Độ tin cậy của mô hình thủy văn NAM trên lưu vực sông Nhuệ - Đáy được chứng minh như thế nào? Mô hình NAM đạt độ tin cậy vượt trội thông qua chỉ số Nash-Sutcliffe đạt 0,85 trong giai đoạn hiệu chỉnh (1972–1974) và 0,86 trong giai đoạn kiểm định (1976–1978) tại trạm Ba Thá. Sai số cân bằng nước tích lũy rất nhỏ, chỉ từ 0,3% đến 0,4%, chứng minh bộ thông số tối ưu phản ánh chính xác cơ chế dòng chảy thực tế.

Biến đổi khí hậu làm gia tăng nguy cơ lũ lụt trên lưu vực sông Nhuệ - Đáy ra sao? Lượng mưa cực trị có xu hướng tập trung mạnh tại các tâm mưa lớn như núi Ba Vì với lượng mưa từ 2200 đến 4000 mm. Dưới kịch bản A1B và A2, lưu lượng đỉnh lũ lớn nhất và dòng chảy mùa lũ (chiếm 70–80% lượng nước cả năm) gia tăng từ 10% đến 25%, làm tăng tần suất xuất hiện các con lũ lớn có biên độ ngập từ 1 đến 4m.

Tình trạng kiệt quệ dòng chảy mùa khô trong tương lai diễn biến nghiêm trọng như thế nào? Trong 7 tháng mùa cạn, lưu lượng dòng chảy vốn chỉ chiếm 20–25% cả năm sẽ tiếp tục giảm sâu từ 8% đến 18%. Điển hình trong năm khô hạn mô phỏng 2040, dòng chảy kiệt sụt giảm trên 20%, kết hợp với lượng nước thải đô thị xả thường xuyên 11–17 m³/s từ sông Tô Lịch sẽ làm gia tăng nguy cơ cạn kiệt và ô nhiễm môi trường nước.

Quy trình xử lý chênh lệch độ phân giải giữa mô hình khí hậu và mô hình thủy văn được thực hiện như thế nào? Tác giả đã chuyển đổi dữ liệu dạng ô lưới 36km từ mô hình RegCM3 về 7 trạm đo khí tượng thông qua phương pháp nội suy song tuyến tính. Tiếp đó, thuật toán hạ quy mô thống kê dựa trên hệ số hiệu sai tỷ lệ theo tháng được áp dụng trên chuỗi số liệu nền 30 năm (1970–1999) để hiệu chỉnh hoàn toàn độ lệch chuẩn trước khi nạp vào mô hình NAM.

Kết luận

Luận văn đã giải quyết trọn vẹn và toàn diện bài toán đánh giá tác động của biến đổi khí hậu đến các đặc trưng cực trị dòng chảy trên lưu vực sông Nhuệ - Đáy thuộc địa bàn Hà Nội. Dưới đây là 5 đóng góp học thuật và thực tiễn cốt lõi của công trình:

  • Xây dựng thành công chuỗi phương pháp luận tích hợp chặt chẽ giữa mô hình khí hậu khu vực RegCM3, phương pháp hạ quy mô thống kê và mô hình thủy văn NAM với độ chính xác Nash-Sutcliffe đạt trên 0,85.
  • Lượng hóa rõ nét bức tranh gia tăng dòng chảy mùa lũ và đỉnh lũ cực đại từ 10% đến 25% cho giai đoạn 2010–2049, cảnh báo nguy cơ ngập lụt gia tăng tại các vùng phân lũ sông Đáy.
  • Xác định mức độ suy thoái dòng chảy tháng kiệt nhất giảm từ 8% đến 18%, cung cấp bằng chứng thực nghiệm về nguy cơ khủng hoảng thiếu nước và suy giảm khả năng tự làm sạch của sông Nhuệ.
  • Thiết lập hệ thống cơ sở dữ liệu thủy văn dự báo theo bước thời gian ngày cho 4 tiểu lưu vực cơ sở phục vụ trực tiếp công tác thiết kế công trình và quy hoạch đê điều của Hà Nội.
  • Đề xuất lộ trình 4 giải pháp công trình và phi công trình đồng bộ với các mục tiêu chỉ số cụ thể cần hoàn thành trong giai đoạn 2025–2035.

Các cơ quan quản lý, viện nghiên cứu và đơn vị tư vấn hạ tầng cần nhanh chóng ứng dụng bộ thông số và kết quả mô phỏng thủy văn của luận văn vào việc cập nhật quy hoạch thoát lũ và kế hoạch hành động thích ứng biến đổi khí hậu của Thủ đô Hà Nội ngay trong giai đoạn hiện nay.