Tổng quan nghiên cứu

Hạ du sông Hồng là trung tâm kinh tế, chính trị và xã hội đặc biệt quan trọng của Việt Nam, bao gồm 10 tỉnh và thành phố với tổng diện tích tự nhiên 14.956,9 km2 cùng quy mô dân số trên 19,2 triệu người, đóng góp khoảng 24% GDP cả nước. Khu vực này tiếp nhận nguồn nước mặt dồi dào từ thượng lưu đổ về trạm thủy văn Sơn Tây với tổng lượng trung bình hàng năm vượt 83 tỷ m3, kết hợp lượng mưa nội vùng đạt xấp xỉ 1.700 mm/năm. Mặc dù vậy, quy luật phân bố dòng chảy tự nhiên bị chi phối mạnh mẽ bởi hoàn lưu gió mùa: hơn 65% tổng lượng dòng chảy tập trung trong 4 tháng mùa lũ từ tháng 6 đến tháng 9, trong khi mùa kiệt kéo dài từ tháng 11 đến tháng 5 năm sau chỉ chiếm chưa đầy 35% lượng nước.

Trong những năm gần đây, sự suy thoái dòng chảy mùa cạn kết hợp tác động của biến đổi khí hậu và khai thác thượng nguồn đã dẫn tới hiện tượng hạ thấp mực nước lịch sử. Điển hình vào mùa kiệt năm 2009-2010, mực nước tại trạm thủy văn Hà Nội tụt sâu xuống mức kỷ lục 0,1m, gây tê liệt mạng lưới giao thông thủy và đe dọa an ninh nguồn nước tưới cho hơn 1,17 triệu ha gieo trồng cây lương thực.

Nghiên cứu tập trung giải quyết vấn đề bất đồng bộ pha cạn giữa các nhánh sông chính (sông Đà, sông Thao, sông Lô) và xác định các tổ hợp cạn bất lợi nhất ảnh hưởng trực tiếp tới dòng chính hạ du. Mục tiêu cốt lõi của luận văn là xây dựng cơ sở khoa học định lượng thông qua chuỗi số liệu quan trắc dài hạn 52 năm (1960-2011) và mô hình hóa thủy lực, từ đó đánh giá mức độ tổn thương của 5 hệ thống thủy nông liên tỉnh cùng 31 tuyến đường thủy nội địa, làm tiền đề cho công tác quy hoạch điều tiết nguồn nước tối ưu.

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 xây dựng trên nền tảng lý thuyết thủy động lực học dòng chảy sông ngòi và phương pháp thống kê xác suất thủy văn chuyên sâu. Mô hình toán sử dụng hệ phương trình Saint-Venant một chiều mô tả dòng chảy không ổn định trong mạng lưới sông hở, kết hợp sơ đồ sai phân hữu hạn 6 điểm ẩn Abbott để giải số trị hệ phương trình vi phân phi tuyến:

  1. Phương trình liên tục biểu diễn sự bảo toàn khối lượng nước qua mặt cắt ngang lòng dẫn theo thời gian và không gian.
  2. Phương trình động lượng mô tả sự biến đổi gia tốc, gradient mực nước, ma sát đáy lòng dẫn và lực cản dòng chảy.

Các chỉ tiêu xác suất thống kê được chuẩn hóa để phân định cấp độ dòng chảy kiệt: thời đoạn cạn kiệt chung khi lưu lượng trung bình nhỏ hơn giá trị ứng với tần suất bảo đảm P = 66%, và thời đoạn đặc biệt cạn khi lưu lượng nhỏ hơn tần suất P = 85%. Trạng thái đồng bộ dòng chảy cạn được phân loại theo tiêu chí: rất đồng bộ khi toàn bộ 3 nhánh thượng lưu xuất hiện dòng chảy kiệt cùng thời điểm, khá đồng bộ khi có 2 trong 3 nhánh trùng pha, và không đồng bộ khi các pha kiệt phân tán. Đồng thời, nghiên cứu lồng ghép lý thuyết động lực học thủy triều vịnh Bắc Bộ với chế độ nhật triều có biên độ truyền sâu vào nội địa tới 150 km trong mùa khô.

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

Nguồn dữ liệu quan trắc phục vụ phân tích bao gồm chuỗi số liệu lưu lượng và mực nước trung bình ngày, trung bình tháng liên tục suốt 52 năm (giai đoạn 1960-2011) tại 4 trạm quan trắc thủy văn xương sống: trạm Hòa Bình (nhánh sông Đà), trạm Yên Bái (nhánh sông Thao), trạm Vụ Quang (nhánh sông Lô) và trạm Sơn Tây - Hà Nội (dòng chính sông Hồng). Cỡ mẫu gồm 624 tháng quan trắc được thu thập theo phương pháp thống kê liên tục toàn diện, bảo đảm tính đại diện tự nhiên trước và sau khi xây dựng các công trình hồ chứa quy mô lớn ở thượng nguồn.

Mô hình thủy lực một chiều MIKE 11-HD do Viện Thủy động lực học Đan Mạch phát triển được lựa chọn làm công cụ tính toán chính. Lý do lựa chọn mô hình này là nhờ khả năng diễn toán xuất sắc các mạng lưới sông dạng mắt lưới phức tạp, chịu tác động đồng thời của dòng chảy kiệt từ thượng lưu và dao động thủy triều mạnh ở cửa sông hạ du.

Mô hình được thiết lập với sơ đồ mặt mạng sông hoàn chỉnh, tích hợp hàng trăm mặt cắt địa hình thực đo dọc các trục sông chính và sông nhánh. Quá trình hiệu chỉnh và kiểm định thông số nhám Manning được thực hiện công phu dựa trên số liệu thực đo các mùa cạn lịch sử 1998-1999, 2003-2004 và 2004-2005 tại các trạm Sơn Tây, Hà Nội, Thượng Cát, Hưng Yên. Kết quả kiểm định đạt chỉ số hiệu quả Nash-Sutcliffe dao động từ 0,85 đến 0,92, khẳng định độ chính xác và tin cậy rất cao của mô hình mô phỏng.

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

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

Phân tích chuỗi số liệu 52 năm chỉ ra rằng sự đồng bộ về thời gian xuất hiện dòng chảy kiệt giữa các nhánh thượng nguồn không diễn ra thường xuyên mà phân hóa rõ rệt. Nhánh sông Đà và sông Thao có tính đồng pha cao nhất với khoảng 70% số năm xuất hiện dòng chảy 3 tháng nhỏ nhất (3Tmin) và dòng chảy 1 tháng nhỏ nhất (Tmin) đồng thời gian. Khi xét tổng thể cả 3 nhánh sông Đà, Thao và Lô, trạng thái khá đồng bộ chiếm tỷ lệ 54,9%, trạng thái rất đồng bộ chiếm 35,29%, và chỉ có 9,8% số năm thuộc trạng thái hoàn toàn không đồng bộ.

Về mức độ cạn kiệt, tổ hợp cạn bất lợi nhất trong lịch sử quan trắc được xác định là mùa cạn các năm 2004-2005 và đặc biệt là năm 2009-2010. Trong năm 2009-2010, cả 4 vị trí kiểm soát Hòa Bình, Yên Bái, Vụ Quang và Sơn Tây đều rơi vào cấp độ đặc biệt cạn ứng với tần suất vượt P = 85%, khiến lưu lượng dòng chảy mùa kiệt suy giảm trên 40% so với trung bình nhiều năm.

Tác động của tổ hợp cạn bất lợi này đến hệ thống thủy nông là cực kỳ nghiêm trọng. Tại vùng hạ du sông Thái Bình, diện tích được tưới chủ động chỉ đạt 57% so với yêu cầu, tồn tại 57.878 ha đất canh tác bị thiếu hụt nguồn nước do các cống đầu mối không thể tự chảy. Vùng hữu sông Hồng chỉ đáp ứng được 64% diện tích tưới thiết kế (đạt 222.089 ha trên tổng số 350.000 ha đất canh tác), các công trình trọng điểm như cống Liên Mạc, cống Cẩm Đình bị phơi đáy hoặc hụt cột nước tĩnh.

Trên mạng lưới giao thông thủy, tổ hợp cạn khắc nghiệt làm suy giảm mớn nước nghiêm trọng trên toàn bộ 31 tuyến đường thủy nội địa quy định tại Thông tư 36/2012/TT-BGTVT. Điển hình tuyến vận tải số 4 (Cửa Đáy - Ninh Bình) độ sâu luồng hạ thấp dưới 0,8m, tuyến sông Đuống chỉ duy trì độ sâu 1,5-1,8m, khiến hơn 45% đội tàu vận tải có trọng tải từ 300 đến 600 DWT phải giảm tải từ 30% đến 50% hoặc buộc phải ngừng hoạt động để tránh mắc cạn.

Thảo luận kết quả

Nguyên nhân chính dẫn tới tình trạng kiệt quệ nguồn nước hạ du bắt nguồn từ sự kết hợp bất lợi giữa suy giảm lượng mưa mùa đông (lượng mưa tháng 12 đến tháng 2 chỉ đạt dưới 30 mm) và hiện tượng tụt sâu đáy sông Hồng do xói lòng dẫn nhân tạo hạ lưu.

Khi thảo luận kết quả mô phỏng, dữ liệu thủy lực được biểu diễn trực quan qua biểu đồ đường quá trình mực nước tại các vị trí đầu mối và biểu đồ không gian phân bố độ sâu cột nước dọc tuyến sông. Biểu đồ tương quan phân tán giữa các trạm thượng nguồn với trạm Sơn Tây cho thấy hệ số tương quan rất cao, đạt R = 0,78 đến 0,89 đối với dòng chảy năm và R = 0,63 đến 0,82 đối với dòng chảy mùa cạn. Điều này chứng minh sự phụ thuộc mật thiết của nguồn nước hạ du vào sự điều tiết của các hồ chứa thượng nguồn. Các phát hiện này hoàn toàn tương thích với các báo cáo khí tượng thủy văn gần đây, khẳng định tính tất yếu của việc phải tính toán tổ hợp bất lợi thay vì chỉ xét riêng lẻ từng nhánh sông độc lập.

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

  1. Vận hành điều tiết xả nước liên hồ chứa thượng lưu: Bộ Nông nghiệp và Phát triển Nông thôn phối hợp cùng Tập đoàn Điện lực Việt Nam thiết lập quy chế xả nước tăng cường từ các hồ thủy điện Hòa Bình, Tuyên Quang và Thác Bà. Mục tiêu duy trì lưu lượng xả tối thiểu tại Sơn Tây trên 1.400 m3/s và mực nước tại Hà Nội đạt cao trình từ 2,2m trở lên trong các đợt đổ ải vụ Đông Xuân (khoảng 15-20 ngày tập trung vào tháng 1 và tháng 2 hàng năm).

  2. Hiện đại hóa và cải tạo hạ tầng công trình thủy lợi: Ủy ban nhân dân 10 tỉnh đồng bằng sông Hồng chỉ đạo các đơn vị quản lý khai thác công trình thủy lợi triển khai hạ thấp cao trình đặt ống hút, chuyển đổi sang máy bơm chìm hoặc bơm trục đứng công suất lớn tại hơn 1.100 trạm bơm vùng tả và hữu sông Hồng. Mục tiêu nâng tỷ lệ diện tích tưới chủ động từ mức 64% hiện nay lên trên 85% trước năm 2030, đồng thời kiên cố hóa 80% hệ thống kênh mương nội đồng để cắt giảm tỷ lệ thất thoát nước từ 35% xuống dưới 15%.

  3. Nạo vét luồng lạch và phân luồng giao thông thủy: Cục Đường thủy nội địa Việt Nam chủ trì kế hoạch nạo vét định kỳ hàng năm tại các điểm nghẽn cạn trọng yếu trên sông Đuống, sông Luộc, sông Đáy và sông Kinh Thầy theo đúng tiêu chuẩn kỹ thuật cấp I và cấp II. Đảm bảo duy trì chiều sâu luồng chạy tàu mùa kiệt tối thiểu từ 2,5m đến 3,5m, đáp ứng an toàn cho phương tiện vận tải từ 600 DWT lưu thông liên tục 24/7.

  4. Xây dựng hệ thống quan trắc và cảnh báo sớm hạn kiệt: Trung tâm Dự báo Khí tượng Thủy văn Quốc gia đầu tư lắp đặt mạng lưới trạm đo mực nước tự động truyền dữ liệu thời gian thực tại 15 điểm nút phân lưu hạ du. Thiết lập mô hình dự báo tổ hợp dòng chảy cạn trước từ 30 đến 45 ngày nhằm cung cấp dữ liệu số phục vụ điều hành sản xuất nông nghiệp và giao thông thủy.

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

  1. Cơ quan quản lý nhà nước về tài nguyên nước và giao thông: Bộ Tài nguyên và Môi trường, Tổng cục Thủy lợi, Cục Đường thủy nội địa Việt Nam sử dụng kết quả nghiên cứu làm căn cứ khoa học trực tiếp để rà soát, điều chỉnh Quy trình vận hành liên hồ chứa trên lưu vực sông Hồng và quy hoạch mạng lưới đường thủy nội địa quốc gia.

  2. Doanh nghiệp và công ty khai thác công trình thủy nông: Các Công ty TNHH MTV Khai thác Công trình Thủy lợi (như hệ thống Bắc Hưng Hải, Bắc Đuống, Sông Nhuệ, Liễn Sơn) áp dụng số liệu dự báo mực nước kiệt để xây dựng phương án vận hành trạm bơm dã chiến, cải tạo bể hút và lên lịch lấy nước đổ ải chuẩn xác.

  3. Các đơn vị vận tải thủy và doanh nghiệp logistics: Các hiệp hội vận tải thủy, chủ tàu thuyền và đơn vị khai thác cảng sông sử dụng các bản đồ cảnh báo vị trí bất lợi và độ sâu luồng chạy tàu để lập kế hoạch xếp dỡ tải trọng an toàn, hạn chế tối đa rủi ro tổn thất do tàu mắc cạn mùa khô.

  4. Viện nghiên cứu, trường đại học và học viên chuyên ngành: Giảng viên, nghiên cứu sinh và học viên cao học thuộc các ngành Thủy văn học, Kỹ thuật Tài nguyên nước, Biến đổi khí hậu khai thác phương pháp luận kết hợp mô hình MIKE 11-HD và lý thuyết tổ hợp xác suất để mở rộng nghiên cứu cho các lưu vực sông khác trên toàn quốc.

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

  1. Tại sao dòng chảy mùa cạn hạ du sông Hồng ngày càng khốc liệt dù đã có hệ thống hồ chứa lớn ở thượng nguồn? Hiện tượng này do tác động kép: sự bất đồng bộ dòng chảy tự nhiên giữa các nhánh sông kết hợp hiện tượng hạ thấp lòng dẫn sông Hồng (tụt đáy sông từ 1,5m đến 2,0m) khiến mực nước tại các trạm hạ du giảm sâu dù lưu lượng xả vẫn duy trì.

  2. Mùa cạn năm nào trong chuỗi quan trắc được đánh giá là bất lợi và khắc nghiệt nhất? Mùa cạn năm 2009-2010 là năm bất lợi nhất khi cả 3 nhánh sông Đà, Thao, Lô cùng rơi vào trạng thái đặc biệt cạn ứng với tần suất P = 85%, khiến mực nước tại Hà Nội xuống mức thấp kỷ lục 0,1m.

  3. Mô hình MIKE 11-HD được kiểm định đạt mức độ tin cậy như thế nào trong nghiên cứu? Mô hình được kiểm định qua 3 mùa cạn thực tế với chuỗi số liệu đo đạc chi tiết tại 4 trạm thủy văn, đạt chỉ số hiệu quả Nash-Sutcliffe từ 0,85 đến 0,92, chứng minh khả năng mô phỏng diễn biến thủy lực chính xác vượt trội.

  4. Hệ thống thủy lợi nào tại hạ du sông Hồng chịu mức độ tổn thương nặng nề nhất? Vùng hạ du sông Thái Bình và vùng hữu sông Hồng chịu tổn thương lớn nhất, trong đó hạ du sông Thái Bình có tới 57.878 ha chưa được tưới chủ động và chỉ đáp ứng được 57% yêu cầu cấp nước do phụ thuộc dòng chảy phân lưu.

  5. Luận văn đã phân tích mạng lưới giao thông thủy dựa trên cơ sở pháp lý và kỹ thuật nào? Nghiên cứu căn cứ vào Thông tư 36/2012/TT-BGTVT của Bộ Giao thông vận tải, phân loại chi tiết cấp kỹ thuật của 31 tuyến đường thủy nội địa trên 25 sông chính và định lượng chính xác các vị trí cạn bất lợi trên từng tuyến.

Kết luận

  • Luận văn đã phân tích toàn diện chuỗi số liệu thủy văn 52 năm (1960-2011), làm sáng tỏ quy luật đồng bộ dòng chảy cạn giữa 3 nhánh sông thượng nguồn với dòng chính hạ du sông Hồng.
  • Xác định chính xác tổ hợp cạn bất lợi cực đoan nhất mùa kiệt 2009-2010 với mực nước Hà Nội tụt đáy 0,1m làm chuẩn tính toán thủy lực.
  • Ứng dụng thành công mô hình số trị MIKE 11-HD đạt hệ số Nash trên 0,85, định lượng hóa mức độ suy giảm khả năng cấp nước của 5 vùng thủy lợi liên tỉnh.
  • Đánh giá chi tiết sự suy giảm độ sâu cột nước chạy tàu trên 31 tuyến đường thủy nội địa, chỉ ra các điểm nghẽn nghiêm trọng cản trở tàu trên 300 DWT.
  • Đề xuất hệ thống giải pháp công trình và phi công trình đồng bộ, tích hợp giữa điều tiết liên hồ chứa, cải tạo trạm bơm và nạo vét luồng lạch.

Đóng góp khoa học cốt lõi của luận văn là thiết lập khung phương pháp luận kết hợp thống kê xác suất tổ hợp với mô hình thủy lực hiện đại, giải quyết căn cơ bài toán quản lý nguồn nước trong điều kiện cạn kiệt. Kế hoạch tiếp theo trong giai đoạn 2025-2030 cần tiếp tục mở rộng mô phỏng với các kịch bản biến đổi khí hậu cực đoan và sự hình thành của các bậc thủy điện mới. Các cơ quan quản lý và đơn vị khai thác cần nhanh chóng ứng dụng các kết quả tính toán này vào công tác điều hành xả nước mùa khô và quy hoạch giao thông thủy bền vững.